CN103630939A - 一种气层识别评价方法 - Google Patents
一种气层识别评价方法 Download PDFInfo
- Publication number
- CN103630939A CN103630939A CN201310548161.XA CN201310548161A CN103630939A CN 103630939 A CN103630939 A CN 103630939A CN 201310548161 A CN201310548161 A CN 201310548161A CN 103630939 A CN103630939 A CN 103630939A
- Authority
- CN
- China
- Prior art keywords
- fluid
- dry
- pore space
- rho
- well
- 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
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开一种气层识别评价方法。所述方法,包括如下步骤:根据研究区域内砂岩储层的不同饱和状态下岩心的声学参数测试结果,计算其孔隙空间模量值,确定不同流体的孔隙空间模量值的参数范围,作为研究区域的气层识别的界限值;对研究区域内的待识别井进行密度、阵列声波测井,按照阵列声波资料处理的标准流程进行计算分析,绘制出单口井的纵波时差曲线及横波时差曲线,通过密度、纵波时差、横波时差曲线计算出待识别井砂岩储层的孔隙空间模量值;将待识别井砂岩储层的孔隙空间模量值与研究区域的气层识别的界限值进行比较,判定所述待识别井砂岩储层的流体性质类型。本发明能够准确反映储层孔隙中的流体特性,提高了流体性质判识的准确率。
Description
技术领域
本发明涉及石油勘探技术领域,特别涉及一种气层识别评价方法。
背景技术
低孔隙度和低渗透率使得储层中的流体对测井响应的贡献大大减小,再加上流体性质和岩性的复杂情况,造成部分测井信息对地层的反映失真,增加了测井评价的难度,因此在油田勘探开发初期低孔低渗气层往往被遗漏。
目前应用常规测井资料进行气层识别评价多采用响应特征法、三孔隙度重叠/差值/比值法、孔隙度与电阻率交会法、束缚水饱和度法、视地层水电阻率法等,这些方法对致密砂岩气层的判识不准确。应用偶极声波测井是偶极横波成像测井仪通过利用单极和偶极测量的组合在硬地层和软地层中都能测得基础的声波数据以及纵、横波和斯通利波等的特征参数,当孔隙介质含气时,其纵波速度明显下降,纵波时差增大,但横波速度下降不明显,纵、横波速度比值(Vp/Vs)将减小,因此可根据测得的声学特征参数来判别致密砂岩储层的流体性质。但目前应用偶极声波测井资料进行气层识别的常规方法多采用交会图、差比值、重叠显示、基线值差异等直观定性解释方法,判识参数选择单一,且没有充分挖掘测井资料中的声学特征信息,致使储层流体识别中误判率较高,使得致密砂岩气层进一步解释评价的不确定性增大。
发明内容
本发明的目的在于提供一种能够对低孔低渗致密砂岩储层的流体性质进行准确判识的气层识别评价方法。
为了达到上述目的,本发明采用的技术方案为:
一种气层识别评价方法,包括如下步骤:
根据研究区域内砂岩储层的不同饱和状态下岩心的声学参数测试结果,计算其孔隙空间模量值,确定不同流体的孔隙空间模量值的参数范围,作为所述研究区域的气层识别的界限值;
对所述研究区域内的待识别井进行密度、阵列声波测井,按照阵列声波资料处理的标准流程进行计算分析,绘制出单口井的纵波时差曲线及横波时差曲线,通过密度、纵波时差、横波时差曲线计算出待识别井砂岩储层的孔隙空间模量值;
将待识别井砂岩储层的孔隙空间模量值与所述研究区域的气层识别的界限值进行比较,判定所述待识别井砂岩储层的流体性质类型。
上述方案中,所述根据研究区域内砂岩储层的不同饱和状态下岩心的声学参数测试结果,计算其孔隙空间模量值,具体包括如下步骤:
将致密砂岩储层看作各向同性的理想介质,根据波动理论,岩石的纵波、横波速度可以表示为:
将饱和流体岩石看作由干燥岩石和流体两部分组成,由Gassmann理论推导出的饱和流体岩石与干燥岩石的体积模量关系式如下,
Ksat=Kdry+β2M (8)
式中Ksat为饱和流体岩石的体积模量;Kdry为干燥岩石的体积模量;β为当孔隙压力为常数时的流体体积变化与岩石体积变化之比;M为地层模量,M是在体积不变的情况下,使流体进入地层的压力。
定义f为孔隙空间模量值,与Biot-Gassmann方程相结合,得到:
由于横波参数基本不受孔隙流体性质的影响,因此干燥岩石与饱和流体岩石的剪切模量基本不变,即
μsat=μdry (10)
将公式(2)、(4)代入公式(1),得到:
由以上公式可以推导出:
其中,干岩样的纵横波速度比值可以通过实验室测试得到。
与现有技术相比,本发明技术方案产生的有益效果为:
本发明通过对储层综合模量值的岩石骨架响应校正,得到反映孔隙中流体性质的体积模量值,根据其分布范围确定储层的流体性质类型。该孔隙空间模量值能够准确反映储层孔隙中的流体特性,提高了流体性质判识的准确率。
附图说明
图1为本发明实施例提供的一种气层识别评价方法的流程图;
图2为本发明实施例中某地区地区干岩样纵横波速度比统计图;
图3为本发明实施例中不同饱和岩样的孔隙空间模量值分布范围统计图
图4为本发明实施例中某井实际资料的处理解释成果图。
具体实施方式
下面结合附图和实施例对本发明的技术方案做进一步说明。
本发明实施例提供一种气层识别评价方法,包括如下步骤:
步骤110,根据研究区域内砂岩储层的不同饱和状态下岩心的声学参数测试结果,计算其孔隙空间模量值,确定不同流体的孔隙空间模量值的参数范围,作为所述研究区域的气层识别的界限值;
步骤120,对所述研究区域内的待识别井进行密度、阵列声波测井,按照阵列声波资料处理的标准流程进行计算分析,绘制出单口井的纵波时差曲线及横波时差曲线,通过密度、纵波时差、横波时差曲线计算出待识别井砂岩储层的孔隙空间模量值;
步骤130,将待识别井砂岩储层的孔隙空间模量值与所述研究区域的气层识别的界限值进行比较,判定所述待识别井砂岩储层的流体性质类型。
本实施例中,根据研究区域内砂岩储层的不同饱和状态下岩心的声学参数测试结果,计算其孔隙空间模量值,具体包括如下步骤:
将致密砂岩储层看作各向同性的理想介质,根据波动理论,岩石的纵波、横波速度可以表示为:
将饱和流体岩石看作由干燥岩石和流体两部分组成,由Gassmann理论推导出的饱和流体岩石与干燥岩石的体积模量关系式如下,
Ksat=Kdry+β2M (14)
式中Ksat为饱和流体岩石的体积模量;Kdry为干燥岩石的体积模量;β为当孔隙压力为常数时的流体体积变化与岩石体积变化之比;M为地层模量,M是在体积不变的情况下,使流体进入地层的压力。
定义f为孔隙空间模量值,与Biot-Gassmann方程相结合,得到:
由于横波参数基本不受孔隙流体性质的影响,因此干燥岩石与饱和流体岩石的剪切模量基本不变,即
μsat=μdry (16)
将公式(2)、(4)代入公式(1),得到:
由以上公式可以推导出:
其中,干岩样的纵横波速度比值可以通过实验室测试得到。
实施例1:
选取某地区为研究的目标区块,该地区盒8、山西组地层是主力产气层。对该地区的126块干岩样,按照《岩石声波特性的实验室测定(SY/T6351-1998)》标准流程进行实验及计算,测量每块岩心水饱和、干燥条件下的纵波速度、横波速度、密度值,求得干岩样的纵波时差、横波时差。
根据实验测试数据统计结果得出,该地区干岩样的纵横波速度比值为1.53,根据公式(6)计算得到C值为2.34,图1为某地区干岩样纵横波速度比统计图。
根据研究区域砂岩储层的气层识别的界限值的计算方法,计算出实验样品点的孔隙空间模量值,图2为绘制的饱和水岩样和干燥岩样的孔隙空间模量值分布范围统计图。从该图中可以看出,两种饱和状态下岩样的孔隙空间模量值存在不同的分布范围,当岩样为水饱和时,其孔隙空间模量值介于10~25GPa之间;当岩样干燥时,其孔隙空间模量值小于5GPa。由此可以确定,对于研究区内的砂岩储层,当计算出的孔隙空间模量值小于5GPa时便认为储层含气。
实际生产应用中,对区域内的井进行密度、阵列声波测井,按照阵列声波资料处理的标准流程进行计算分析,绘制出单口井的纵波时差曲线及横波时差曲线;通过密度、纵波时差、横波时差曲线计算出实际砂岩储层的孔隙空间模量曲线,将孔隙空间模量值小于5GPa的储层段划分为气层。
实际工作中,孔隙空间模量的计算通过编写程序实现,图3为某井实际资料的处理解释成果图。井段2895~2897.4m的孔隙空间模量值较上下围岩储层出现明显减小的趋势,根据地层产气的实际情况,判定该层段为主要的含气性储层段;后期生产中,对2895~2899m井段进行压裂施工测试,测试结果证实该层段出气12068m3/d,验证了该方法的有效性。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种气层识别评价方法,其特征在于,包括如下步骤:
根据研究区域内砂岩储层的不同饱和状态下岩心的声学参数测试结果,计算其孔隙空间模量值,确定不同流体的孔隙空间模量值的参数范围,作为所述研究区域的气层识别的界限值;
对所述研究区域内的待识别井进行密度、阵列声波测井,按照阵列声波资料处理的标准流程进行计算分析,绘制出单口井的纵波时差曲线及横波时差曲线,通过密度、纵波时差、横波时差曲线计算出待识别井砂岩储层的孔隙空间模量值;
将待识别井砂岩储层的孔隙空间模量值与所述研究区域的气层识别的界限值进行比较,判定所述待识别井砂岩储层的流体性质类型。
2.如权利要求1所述的气层识别评价方法,其特征在于,所述根据研究区域内砂岩储层的不同饱和状态下岩心的声学参数测试结果,计算其孔隙空间模量值,具体包括如下步骤:
将致密砂岩储层看作各向同性的理想介质,根据波动理论,岩石的纵波、横波速度可以表示为:
将饱和流体岩石看作由干燥岩石和流体两部分组成,由Gassmann理论推导出的饱和流体岩石与干燥岩石的体积模量关系式如下,
Ksat=Kdry+β2M (2)
式中Ksat为饱和流体岩石的体积模量;Kdry为干燥岩石的体积模量;β为当孔隙压力为常数时的流体体积变化与岩石体积变化之比;M为地层模量,M是在体积不变的情况下,使流体进入地层的压力;
定义f为孔隙空间模量值,与Biot-Gassmann方程相结合,得到:
由于横波参数基本不受孔隙流体性质的影响,因此干燥岩石与饱和流体岩石的剪切模量基本不变,即
μsat=μdry (4)
将公式(2)、(4)代入公式(1),得到:
由以上公式可以推导出:
其中,干岩样的纵横波速度比值可以通过实验室测试得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310548161.XA CN103630939B (zh) | 2013-11-06 | 2013-11-06 | 一种气层识别评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310548161.XA CN103630939B (zh) | 2013-11-06 | 2013-11-06 | 一种气层识别评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103630939A true CN103630939A (zh) | 2014-03-12 |
CN103630939B CN103630939B (zh) | 2018-07-13 |
Family
ID=50212187
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310548161.XA Active CN103630939B (zh) | 2013-11-06 | 2013-11-06 | 一种气层识别评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103630939B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105003257A (zh) * | 2015-08-07 | 2015-10-28 | 中国海洋石油总公司 | 一种定性识别高温高压甲烷气层与二氧化碳气层的方法 |
CN105549088A (zh) * | 2014-10-29 | 2016-05-04 | 中国石油天然气股份有限公司 | 裂缝性致密砂岩中气层的识别方法和装置 |
CN105589110A (zh) * | 2015-12-10 | 2016-05-18 | 中国石油集团川庆钻探工程有限公司 | 孔隙流体压缩系数判别致密砂岩气层的方法 |
CN110133724A (zh) * | 2019-05-06 | 2019-08-16 | 中国石油天然气集团有限公司 | 一种应用阵列声波测井数据进行气层识别的方法 |
CN111077589A (zh) * | 2018-10-19 | 2020-04-28 | 中国石油天然气股份有限公司 | 干岩石纵横波速度比的确定方法和装置 |
CN111271057A (zh) * | 2020-03-25 | 2020-06-12 | 西北大学 | 一种测井识别气层方法 |
CN111625750A (zh) * | 2019-02-28 | 2020-09-04 | 中国石油化工股份有限公司 | 一种针对致密砂岩储层气层与气水同层的识别方法 |
CN111781663A (zh) * | 2020-07-30 | 2020-10-16 | 陕西延长石油(集团)有限责任公司 | 一种快速判断砂岩储层流体性质的方法及模型 |
CN112133376A (zh) * | 2020-08-27 | 2020-12-25 | 中国石油天然气集团有限公司 | 一种饱含水地层体积模量的流体识别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4399525A (en) * | 1979-10-05 | 1983-08-16 | Chevron Research Company | Method for interpreting well log records to yield indications of gas/oil in an earth formation such as a sandstone, limestone, or dolostone |
CN101000378A (zh) * | 2006-01-10 | 2007-07-18 | 中国石油天然气集团公司 | 利用全波列、偶极横波测井资料确定气层的方法 |
CN101634717A (zh) * | 2009-08-26 | 2010-01-27 | 中国石油大学(华东) | 基于测井和叠前道集地震数据的精细横波阻抗求取技术 |
CN102012526A (zh) * | 2010-09-09 | 2011-04-13 | 四川德阳西德电器有限公司 | 采用电阻率资料进行储层流体类型判别的方法 |
US8547788B2 (en) * | 2010-05-17 | 2013-10-01 | Schlumberger Technology Corporation | Methods for making acoustic anisotropy logging while drilling measurements |
-
2013
- 2013-11-06 CN CN201310548161.XA patent/CN103630939B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4399525A (en) * | 1979-10-05 | 1983-08-16 | Chevron Research Company | Method for interpreting well log records to yield indications of gas/oil in an earth formation such as a sandstone, limestone, or dolostone |
CN101000378A (zh) * | 2006-01-10 | 2007-07-18 | 中国石油天然气集团公司 | 利用全波列、偶极横波测井资料确定气层的方法 |
CN101634717A (zh) * | 2009-08-26 | 2010-01-27 | 中国石油大学(华东) | 基于测井和叠前道集地震数据的精细横波阻抗求取技术 |
US8547788B2 (en) * | 2010-05-17 | 2013-10-01 | Schlumberger Technology Corporation | Methods for making acoustic anisotropy logging while drilling measurements |
CN102012526A (zh) * | 2010-09-09 | 2011-04-13 | 四川德阳西德电器有限公司 | 采用电阻率资料进行储层流体类型判别的方法 |
Non-Patent Citations (1)
Title |
---|
李山生 等: ""利用弹性参数识别气、水层"", 《西南石油大学学报(自然科学版)》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105549088A (zh) * | 2014-10-29 | 2016-05-04 | 中国石油天然气股份有限公司 | 裂缝性致密砂岩中气层的识别方法和装置 |
CN105549088B (zh) * | 2014-10-29 | 2018-01-05 | 中国石油天然气股份有限公司 | 裂缝性致密砂岩中气层的识别方法和装置 |
CN105003257A (zh) * | 2015-08-07 | 2015-10-28 | 中国海洋石油总公司 | 一种定性识别高温高压甲烷气层与二氧化碳气层的方法 |
CN105589110A (zh) * | 2015-12-10 | 2016-05-18 | 中国石油集团川庆钻探工程有限公司 | 孔隙流体压缩系数判别致密砂岩气层的方法 |
CN111077589A (zh) * | 2018-10-19 | 2020-04-28 | 中国石油天然气股份有限公司 | 干岩石纵横波速度比的确定方法和装置 |
CN111625750A (zh) * | 2019-02-28 | 2020-09-04 | 中国石油化工股份有限公司 | 一种针对致密砂岩储层气层与气水同层的识别方法 |
CN111625750B (zh) * | 2019-02-28 | 2023-04-28 | 中国石油化工股份有限公司 | 一种针对致密砂岩储层气层与气水同层的识别方法 |
CN110133724A (zh) * | 2019-05-06 | 2019-08-16 | 中国石油天然气集团有限公司 | 一种应用阵列声波测井数据进行气层识别的方法 |
CN111271057A (zh) * | 2020-03-25 | 2020-06-12 | 西北大学 | 一种测井识别气层方法 |
CN111271057B (zh) * | 2020-03-25 | 2023-05-16 | 中国石油大学(北京) | 一种测井识别气层方法 |
CN111781663A (zh) * | 2020-07-30 | 2020-10-16 | 陕西延长石油(集团)有限责任公司 | 一种快速判断砂岩储层流体性质的方法及模型 |
CN112133376A (zh) * | 2020-08-27 | 2020-12-25 | 中国石油天然气集团有限公司 | 一种饱含水地层体积模量的流体识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103630939B (zh) | 2018-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103630939A (zh) | 一种气层识别评价方法 | |
US11098565B2 (en) | Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves | |
CN103513270B (zh) | 一种基于岩石声学特性的气层识别评价方法及装置 | |
CN103603659B (zh) | 利用常规测井资料识别复杂孔隙结构储层流体类型的方法 | |
CN112526107B (zh) | 裂缝型致密砂岩储层甜点识别与定量表征的方法 | |
Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
CN110320139A (zh) | 缝洞型储层裂缝孔隙度定量评价方法及*** | |
CN104252007A (zh) | 一种相容性岩石物理建模方法 | |
CN111271057A (zh) | 一种测井识别气层方法 | |
CN107831540A (zh) | 储层物性参数直接提取新方法 | |
CN107829731A (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
CN103792575B (zh) | 一种去骨架影响的致密砂岩气层识别评价方法 | |
CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
US6374186B1 (en) | Method for overpressure detection from compressional-and- shear-wave data | |
CN116641702A (zh) | 一种页岩气藏甜点区识别方法 | |
CN104484573B (zh) | 一种确定地层刚性系数的方法 | |
CN113703052B (zh) | 海上中、高孔渗砂岩压裂效果预评估方法 | |
FANG et al. | Quantitative Method of Classification and Discrimination of a Porous Carbonate Reservoir Integrating K‐means Clustering and Bayesian Theory | |
CN111594155B (zh) | 一种致密复杂岩性储层气层、水层识别方法 | |
CN114859407A (zh) | 火山岩储层声学特征参数的确定方法及装置 | |
CN110297264B (zh) | 一种低渗气藏薄储层甜点地震预测方法 | |
CN107764697A (zh) | 基于孔隙介质渐进方程非线性反演的含气性检测方法 | |
CN113960666B (zh) | 一种基于叠前弹性参数的基岩气层识别方法及装置 | |
CN108931813B (zh) | 一种识别礁盖、礁核、礁基的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |