CN101000378A - 利用全波列、偶极横波测井资料确定气层的方法 - Google Patents
利用全波列、偶极横波测井资料确定气层的方法 Download PDFInfo
- Publication number
- CN101000378A CN101000378A CN 200610000623 CN200610000623A CN101000378A CN 101000378 A CN101000378 A CN 101000378A CN 200610000623 CN200610000623 CN 200610000623 CN 200610000623 A CN200610000623 A CN 200610000623A CN 101000378 A CN101000378 A CN 101000378A
- Authority
- CN
- China
- Prior art keywords
- gas
- time difference
- wave
- saturation
- compressional wave
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及碎屑岩常规、非常规气藏勘探开发中利用横纵波时差测井资料确定低孔、低渗、低阻气层的方法。根据已知的岩石物理特征建立勘探地区含气储层多种含气饱和度的横—纵波转换交会图,然后利用声波资料横波时差确定多种含气饱和度情况下纵波时差,实际测量纵波时差并与利用声波资料横波时差确定多种含气饱和度情况下纵波时差比较,对应交会图中多种含气饱和度确定储层的含气性。本发明确定疑难气层是利用声波横纵波转换模型,回避了疑难气层的电性、放射性测井响应特征不明显,地提高了测井确定气层的能力。同时,也提高了全波列、偶极横波资料的应用,增强了解决地质问题的能力。
Description
技术领域
本发明涉及碎屑岩常规、非常规气藏勘探开发中利用横纵波时差测井资料确定低孔、低渗、低阻气层的方法。
技术背景
随着气田勘探开发的深入,确定碎屑岩疑难气层的问题日益突出,重点表现在低孔、低渗和低阻气层的确定和识别上。测井识别气层基本方法是根据含气储层的电性和放射性特征变化。常规典型含气储层表现为电阻率升高,声波时差增加,中子孔隙度降低和密度孔隙度升高等测井特征,中子孔隙度-密度孔隙度曲线组合差异明显。但对低孔、低渗和低阻气层,由于储层物性差,骨架对测井信息的贡献远大于孔隙部分的贡献,增加了测井资料识别孔隙流体的困难,电阻率与围岩接近,中子-密度交会差异基本重合,常规资料识别困难,这是测井确定(寻找)气层的难点。
目前为解决这个问题,只能借助于声波全波列、偶极横波、核磁共振等成像测井技术。核磁共振测井对储层物性要求较高,限制了它的应用范围。全波列、偶极横波等测井采集技术提供了丰富信息,但没有真正地深掘这些资料背后的有效信息,因此,来自于全波列和偶极横波测井系列采集的全波列、偶极横波基础测井资料一直也没有取得明显的效果。
发明内容
本发明目的在于借助于声波全波列、偶极横波等成像测井技术,利用测井采集技术提供的丰富信息,提供一种利用全波列、偶极横波测井资料提高确定疑难气层的方法。
本发明提供以下技术方案。
利用全波列、偶极横波测井资料确定气层的方法,根据已知的岩石物理特征建立勘探地区含气储层多种含气饱和度的横-纵波转换交会图,然后利用声波资料横波时差确定多种含气饱和度情况下纵波时差,实际测量纵波时差并与利用声波资料横波时差确定多种含气饱和度情况下纵波时差比较,对应交会图中多种含气饱和度确定储层的含气性。
本发明还提供以下技术方案。
建立多种含气饱和度情况下储层横-纵波转换交会图是:用常规的方法对岩心进行气水两相条件下驱替,在多种状态采用超声脉冲透射法测量岩心纵横波速度,所述的多种状态的含气饱和度覆盖了从完全含水到高饱和气的范围,为0-80%。
本发明采用以下步骤:
a.选择勘探目标地区的含气储层岩心进行确定岩石物理性,包括采用常规的方法对岩心进行洗油洗盐预处理,测量孔隙度和渗透率参数,然后加压饱和对岩心进行气水两相条件下驱替,在0、20%、50%和80%四种不同含气饱和度状态采用超声脉冲透射法测量岩心纵横波速度和时差,纵横波速度采用以下公式计算:
,式中:V:纵波、横波速度,L为样品长度,T为声波传播时间,T0为***零延时;饱和度的计算公式:Sw=(G-Gma)/(Gw-Gma),式中:Sw:含水饱和度,G为岩样当前重量,Gma为岩样完全干燥重量,Gw为完全饱和水岩样重量;
b.所述的不同饱和度覆盖从完全含水到高饱和气的范围,采用常规的方法建立完全含水、低含气饱和度、高含气饱和度几种情况下的横-纵波时差交会图,对其数据进行拟合,得到不同含气饱和度情况下横-纵波时差的关系图;
c.根据不同含气饱和度情况下横-纵波时差的关系,采用常规的测井方法对勘探目标区气井处理成图或资料,再提取纵波时差曲线与计算的多种含气饱和度横-纵波时差的曲线比较,当提取纵波时差曲线位于多种含气饱和度横-纵波时差的曲线之间,或与含气饱和度横-纵波时差的曲线重合时,确定相应的实际储层含气量。
所述的多种含气饱和度为0、20%、50%和80%。
本发明采用的技术方法是以含气储层岩石物理性质为基础,通过测量分析岩心声波纵横波时差(速度)关系,建立了0、20%、50%和80%四种含气饱和度情况下横波纵波的交会图,从而实现不同含气饱和度情况下从横波时差到纵波时差的转换。在实际应用中,利用从全波列、偶极横波测井资料提取出的横波时差,根据上述模型计算得到四种含气饱和度的纵波时差,将提取出的纵波时差与计算的多种含气饱和度情况下纵波时差比较,由于含气饱和度升高引起纵波时差降低,计算纵波时差覆盖了从完全含水到高饱和气的范围,实际纵波时差必然位于计算的多种纵波时差之中,它们之间相对差异大小就可以定量确定储层的含气性。这种方法避开了疑难气层的电性和中子、密度的响应特征不明显难题,利用含气层的声波响应规律,是一种确定疑难气层的有效方法。
本发明确定疑难气层是利用声波横纵波转换模型,回避了疑难气层的电性、放射性测井响应特征不明显,地提高了测井确定气层的能力。同时,也提高了全波列、偶极横波资料的应用,增强了解决地质问题的能力。
附图说明
附图1横-纵波时差在不同含气饱和度下交会图;
附图2利用横-纵波时差转换模型对疑难气层处理成果图。
图1中横坐标DTS是是测量的横波时差(单位:μs/m),纵坐标DTC是测量的纵波时差(单位:μs/m)。图中菱形、叉形、方形和三角形分别代表0、20%、50%和80%四种含气饱和度下测量的数据点;对应线分别代表根据0、20%、50%和80%四种含气饱和度测量数据点拟合出的数据线。
图2是利用偶极横波资料,根据横-纵波时差转换模型对低阻气层的处理成果图。第五道中,DTC为从偶极横波资料中提取出的纵波时差,DTC0、DTC2、DTC8分别为根据横-纵波时差转换模型计算得出的0、20%和80%含气饱和度的纵波时差。充填部分为左边DTC,右边为DTC0。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
1、通过含气储层岩心计算不同饱和度情况下的纵横波时差
首先选择研究地区的含气储层岩心进行岩石物理性质实验。对岩心进行洗油洗盐等基本预处理,测量孔隙度和渗透率等基本参数,然而加压饱和对岩心进行气水两相条件下驱替实验。依据行业标准“SY/T 6351-1998岩石声波特性的实验室测定方法”,在0、20%、50%和80%四种不同含气饱和度状态采用超声脉冲透射法测量岩心纵横波速度(速度的倒数即得到时差)。
纵横波速度的计算公式为:
式中:V:纵波、横波速度;
L为样品长度;
T为声波传播时间;
T0为***零延时。
饱和度的计算公式:
Sw=(G-Gma)/(Gw-Gma) (2)
式中:Sw:含水饱和度;
G为岩样当前重量;
Gma为岩样完全干燥重量;
Gw为完全饱和水岩样重量。
2、建立横-纵波时差转换交会图
第一步确定了四种含气饱和度为0、20%、50%和80%,饱和度覆盖了从完全含水到高饱和气的范围。建立完全含水、低含气饱和度、高含气饱和度等几种情况下的横-纵波时差交会图,见图1。对其数据进行拟合,得到0、20%、50%和80%四种含气饱和度情况下横-纵波时差的关系。
DTC0=a0ln(DTS)-b0 (3)
DTC2=a2ln(DTS)-b2 (4)
DTC5=a5ln(DTS)-b5 (5)
DTC8=a8ln(DTS)-b8 (6)
式中:DTS:横波时差(μs/m);
DTC0、DTC2、DTC5、DTC8:分别为0、20%、50%和80%四种含气饱和度情况下纵波时差(μs/m);
a0、a2、a5、a8、a0、b0、b2、b5、b8:分别为拟合公式的系数。
3、对疑难气层进行确定(寻找)
以一实例来说明此方法确定(寻找)疑难气层的步骤,见图2。
(1)地区基本情况
图2显示的是对某地区一口井疑难气层处理成果图,上面含气储层岩心取自该地区。该地区是一个典型低孔、低渗岩性气藏。储层岩性以石英砂岩和岩屑砂岩为主。此井气层可明显分为两类,一类是正常气层,电性和放射性特征明显,如图2中3350-3358m储层;另一类是低阻气藏,它与围岩电性特征很接近不能区分,同时中子-密度曲线基本重合,常规测井资料不能识别疑难气层,如图2中334l-3349m储层。
(2)图2中曲线说明
左边第一道为测井仪器测量的曲线,其中:绿线为仪器测量的自然电位曲线,蓝线为井径曲线,红线为自然伽马曲线;
左边第二道为曲线深度道;
左边第三道为测井仪器测量的曲线,其中:绿线为微球曲线,蓝线为浅侧向曲线,红线为深侧向曲线;
左边第四道为测井仪器测量的曲线,其中:绿线为中子曲线,蓝线为声波曲线,红线为密度曲线;
左边第五道中,DTC为全波列或偶极横波测井测量的纵波时差曲线,蓝线DTC0为根据上述方法计算得到的100%含水的纵波时差曲线,绿线DTC2为根据上述方法计算得到的20%含气饱和度的纵波时差曲线,紫线DTC8为根据上述方法计算得到的80%含气饱和度的纵波时差曲线。
(3)确定(寻找)疑难气层方法
从全波列、偶极横波测井资料提取出纵横波时差,根据上述四种交会计算得到四种含气饱和度的纵波时差,图中显示了三种饱和度的纵波时差,分为DTC0、DTC2和DTC8。将提取出的纵波时差DTC与计算的多种含气饱和度纵波时差比较,从图上左边第五道可以看出,在3341-3349m疑难含气储层DTC位于DTC0和DTC8之间,说明实际储层含气,含气饱和度在0-80%之间,显示明显。因此通过全波列、偶极横波测井的纵波时差DTC与上述方法计算得到的四种饱和度相对差异大小就可以定性到半定量确定(寻找)疑难含气储层。
Claims (4)
1、一种利用全波列、偶极横波测井资料确定气层的方法,其特征是:根据已知的岩石物理特征建立勘探地区含气储层多种含气饱和度的横-纵波转换交会图,然后利用声波资料横波时差确定多种含气饱和度情况下纵波时差,实际测量纵波时差并与利用声波资料横波时差确定多种含气饱和度情况下纵波时差比较,对应交会图中多种含气饱和度确定储层的含气性。
2、根据权利要求1所述的一种利用全波列、偶极横波测井资料确定气层的方法,其特征是,建立多种含气饱和度情况下储层横-纵波转换交会图是:用常规的方法对岩心进行气水两相条件下驱替,在多种状态采用超声脉冲透射法测量岩心纵横波速度,所述的多种状态的含气饱和度覆盖了从完全含水到高饱和气的范围,为0-80%。
3、根据权利要求1或2所述的一种利用全波列、偶极横波测井资料确定气层的方法,其特征是采用以下步骤:
a.选择勘探目标地区的含气储层岩心进行确定岩石物理性,包括采用常规的方法对岩心进行洗油洗盐预处理,测量孔隙度和渗透率参数,然后加压饱和对岩心进行气水两相条件下驱替,在0、20%、50%和80%四种不同含气饱和度状态采用超声脉冲透射法测量岩心纵横波速度和时差,纵横波速度采用以下公式计算:
,式中:V:纵波、横波速度,L为样品长度,T为声波传播时间,T0为***零延时;饱和度的计算公式:Sw=(G-Gma)/(Gw-Gma),式中:Sw:含水饱和度,G为岩样当前重量,Gma为岩样完全干燥重量,Gw为完全饱和水岩样重量;
b.所述的不同饱和度覆盖从完全含水到高饱和气的范围,采用常规的方法建立完全含水、低含气饱和度、高含气饱和度几种情况下的横-纵波时差交会图,对其数据进行拟合,得到不同含气饱和度情况下横-纵波时差的关系图;
c.根据不同含气饱和度情况下横-纵波时差的关系,采用常规的测井方法对勘探目标区气井处理成图,再提取纵波时差曲线与计算的多种含气饱和度横-纵波时差的曲线比较,当提取纵波时差曲线位于多种含气饱和度横-纵波时差的曲线之间,或与含气饱和度横-纵波时差的曲线重合时,确定相应的实际储层含气量。
4、根据权利要求3所述的一种利用全波列、偶极横波测井资料确定气层的方法,其特征是:所述的多种含气饱和度为0、20%、50%和80%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100006234A CN100429529C (zh) | 2006-01-10 | 2006-01-10 | 利用全波列、偶极横波测井资料确定气层的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100006234A CN100429529C (zh) | 2006-01-10 | 2006-01-10 | 利用全波列、偶极横波测井资料确定气层的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101000378A true CN101000378A (zh) | 2007-07-18 |
CN100429529C CN100429529C (zh) | 2008-10-29 |
Family
ID=38692428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100006234A Active CN100429529C (zh) | 2006-01-10 | 2006-01-10 | 利用全波列、偶极横波测井资料确定气层的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100429529C (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101158724B (zh) * | 2007-09-14 | 2010-06-09 | 中国石油集团西北地质研究所 | 基于偶极小波的储层厚度预测方法 |
CN101520517B (zh) * | 2008-02-25 | 2011-06-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种能准确评价碎屑岩盆地含油气目标的方法 |
CN102508299A (zh) * | 2011-09-29 | 2012-06-20 | 中国石油大学(华东) | 在井中发射并接收偶极横波的远探测方法 |
CN102721984A (zh) * | 2012-05-28 | 2012-10-10 | 中国海洋石油总公司 | 一种预测碎屑岩隐蔽储层的方法 |
CN103630939A (zh) * | 2013-11-06 | 2014-03-12 | 中国石油天然气集团公司 | 一种气层识别评价方法 |
CN104514553A (zh) * | 2013-09-29 | 2015-04-15 | 中国石油化工股份有限公司 | 一种岩屑砂岩层中气层的识别方法及其应用 |
CN105549088A (zh) * | 2014-10-29 | 2016-05-04 | 中国石油天然气股份有限公司 | 裂缝性致密砂岩中气层的识别方法和装置 |
CN106842303A (zh) * | 2016-12-30 | 2017-06-13 | 中国石油天然气股份有限公司 | 含气饱和度确定方法及装置 |
CN107238521A (zh) * | 2017-06-07 | 2017-10-10 | 中国石油大学(华东) | 一种冻土带水合物岩心的制作方法 |
CN112014878A (zh) * | 2019-05-30 | 2020-12-01 | 中石化石油工程技术服务有限公司 | 一种气水层识别方法 |
CN112065361A (zh) * | 2020-09-16 | 2020-12-11 | 成都理工大学 | 一种基于声波衰减的致密储层含气饱和度确定方法 |
CN112965106A (zh) * | 2021-01-21 | 2021-06-15 | 中国地质大学(北京) | 利用纵横波速度的奇异性识别天然气水合物赋存类型方法 |
CN113253350A (zh) * | 2021-05-25 | 2021-08-13 | 电子科技大学 | 一种基于联合字典的孔隙度反演方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7382135B2 (en) * | 2003-05-22 | 2008-06-03 | Schlumberger Technology Corporation | Directional electromagnetic wave resistivity apparatus and method |
-
2006
- 2006-01-10 CN CNB2006100006234A patent/CN100429529C/zh active Active
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101158724B (zh) * | 2007-09-14 | 2010-06-09 | 中国石油集团西北地质研究所 | 基于偶极小波的储层厚度预测方法 |
CN101520517B (zh) * | 2008-02-25 | 2011-06-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种能准确评价碎屑岩盆地含油气目标的方法 |
CN102508299A (zh) * | 2011-09-29 | 2012-06-20 | 中国石油大学(华东) | 在井中发射并接收偶极横波的远探测方法 |
CN102721984B (zh) * | 2012-05-28 | 2014-07-23 | 中国海洋石油总公司 | 一种预测碎屑岩隐蔽储层的方法 |
CN102721984A (zh) * | 2012-05-28 | 2012-10-10 | 中国海洋石油总公司 | 一种预测碎屑岩隐蔽储层的方法 |
CN104514553A (zh) * | 2013-09-29 | 2015-04-15 | 中国石油化工股份有限公司 | 一种岩屑砂岩层中气层的识别方法及其应用 |
CN104514553B (zh) * | 2013-09-29 | 2017-08-25 | 中国石油化工股份有限公司 | 一种岩屑砂岩层中气层的识别方法及其应用 |
CN103630939A (zh) * | 2013-11-06 | 2014-03-12 | 中国石油天然气集团公司 | 一种气层识别评价方法 |
CN103630939B (zh) * | 2013-11-06 | 2018-07-13 | 中国石油天然气集团公司 | 一种气层识别评价方法 |
CN105549088A (zh) * | 2014-10-29 | 2016-05-04 | 中国石油天然气股份有限公司 | 裂缝性致密砂岩中气层的识别方法和装置 |
CN105549088B (zh) * | 2014-10-29 | 2018-01-05 | 中国石油天然气股份有限公司 | 裂缝性致密砂岩中气层的识别方法和装置 |
CN106842303B (zh) * | 2016-12-30 | 2019-07-05 | 中国石油天然气股份有限公司 | 含气饱和度确定方法及装置 |
CN106842303A (zh) * | 2016-12-30 | 2017-06-13 | 中国石油天然气股份有限公司 | 含气饱和度确定方法及装置 |
CN107238521A (zh) * | 2017-06-07 | 2017-10-10 | 中国石油大学(华东) | 一种冻土带水合物岩心的制作方法 |
CN107238521B (zh) * | 2017-06-07 | 2019-09-10 | 中国石油大学(华东) | 一种冻土带水合物岩心的制作方法 |
CN112014878A (zh) * | 2019-05-30 | 2020-12-01 | 中石化石油工程技术服务有限公司 | 一种气水层识别方法 |
CN112014878B (zh) * | 2019-05-30 | 2023-08-18 | 中石化石油工程技术服务有限公司 | 一种气水层识别方法 |
CN112065361A (zh) * | 2020-09-16 | 2020-12-11 | 成都理工大学 | 一种基于声波衰减的致密储层含气饱和度确定方法 |
CN112065361B (zh) * | 2020-09-16 | 2021-03-12 | 成都理工大学 | 一种基于声波衰减的致密储层含气饱和度确定方法 |
CN112965106A (zh) * | 2021-01-21 | 2021-06-15 | 中国地质大学(北京) | 利用纵横波速度的奇异性识别天然气水合物赋存类型方法 |
CN113253350A (zh) * | 2021-05-25 | 2021-08-13 | 电子科技大学 | 一种基于联合字典的孔隙度反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100429529C (zh) | 2008-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100429529C (zh) | 利用全波列、偶极横波测井资料确定气层的方法 | |
Conyers | Ground-penetrating radar for archaeology | |
Illman | Hydraulic tomography offers improved imaging of heterogeneity in fractured rocks | |
CN105426612A (zh) | 一种地层组分最优化确定方法及装置 | |
CN103163553B (zh) | 基于多重孔隙介质模型的地震烃类检测方法及装置 | |
CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
CN103675907A (zh) | 一种基于岩相约束的avo反演油气检测方法 | |
CN107229076B (zh) | 一种基于测井资料进行温度响应特征分析的方法 | |
Barrio et al. | CO2 migration monitoring methodology in the shallow subsurface: Lessons learned from the CO2FIELDLAB project | |
Li et al. | Exploration methods for late Quaternary shallow biogenic gas reservoirs in the Hangzhou Bay area, eastern China | |
Chen et al. | Geochemical characterization using geophysical data and Markov Chain Monte Carlo methods: A case study at the South Oyster bacterial transport site in Virginia | |
CN106285623B (zh) | 确定总有机碳含量的方法及*** | |
Michael et al. | The South West Hub In-Situ Laboratory–a facility for CO2 injection testing and monitoring in a fault zone | |
CN110018527B (zh) | 幔源油气动态油气藏勘探方法 | |
Flinchum et al. | Identifying recharge under subtle ephemeral features in a flat-lying semi-arid region using a combined geophysical approach | |
Hyndman et al. | Hydrogeophysical case studies at the local scale: The saturated zone | |
Young et al. | 3-D dip filtering and coherence applied to GPR data: A study | |
George et al. | Estimation of aquifer hydraulic parameters via complementing surfacial geophysical measurement by laboratory measurements on the aquifer core samples | |
Roach et al. | Applied geophysics for cover thickness mapping in the southern Thomson Orogen | |
Arts et al. | Results of geophysical monitoring over a “leaking” natural analogue site in Italy | |
CN106353813A (zh) | 基于阵列声波测井的流体性质识别方法 | |
Popp et al. | Assessment of shallow subsurface characterisation with non-invasive geophysical methods at the intermediate hill-slope scale | |
Oldenborger et al. | Radar determination of the spatial structure of hydraulic conductivity | |
Alhumimidi et al. | Imaging fracture distributions of the Al-Khuff Formation outcrops using GPR and ERT geophysical techniques, Al-Qassim area, Saudi Arabia | |
Tian et al. | Prediction study of hydrocarbon reservoir based on time–frequency domain electromagnetic technique taking Ili Basin as an example |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |