CN108549103A - 饱和多孔介质叠后地震烃类检测方法 - Google Patents
饱和多孔介质叠后地震烃类检测方法 Download PDFInfo
- Publication number
- CN108549103A CN108549103A CN201810283307.5A CN201810283307A CN108549103A CN 108549103 A CN108549103 A CN 108549103A CN 201810283307 A CN201810283307 A CN 201810283307A CN 108549103 A CN108549103 A CN 108549103A
- Authority
- CN
- China
- Prior art keywords
- seismic
- oil
- amplitude
- porous medium
- gas
- 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.)
- Pending
Links
- 239000004215 Carbon black (E152) Substances 0.000 title claims abstract description 23
- 238000001514 detection method Methods 0.000 title claims abstract description 23
- 229930195733 hydrocarbon Natural products 0.000 title claims abstract description 23
- 150000002430 hydrocarbons Chemical class 0.000 title claims abstract description 20
- 229920006395 saturated elastomer Polymers 0.000 title claims abstract description 14
- 239000012530 fluid Substances 0.000 claims abstract description 43
- 230000000694 effects Effects 0.000 claims abstract description 29
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 12
- 238000001228 spectrum Methods 0.000 claims abstract description 12
- 238000000265 homogenisation Methods 0.000 claims abstract description 8
- 230000002159 abnormal effect Effects 0.000 claims abstract description 4
- 238000004458 analytical method Methods 0.000 claims description 16
- 238000000034 method Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 description 10
- 239000011435 rock Substances 0.000 description 7
- 239000011148 porous material Substances 0.000 description 6
- 125000001183 hydrocarbyl group Chemical group 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种饱和多孔介质叠后地震烃类检测方法，针对目的层段寻找含油气井和非含油气井出现振幅异常的敏感频率段，再提取新的敏感频带地震体，在振幅均一化的基础上进行差值运算，得到流体活动属性体。本发明将一个针对饱和多孔介质的实验室基础理论属性表达式转换成易实现的叠后地震数据烃类检测方法，且利用子波分解重构生成的敏感频段地震体所求的流体活动属性相较于频谱分解的单频体稳定性更好。
Description
技术领域
本发明涉及一种饱和多孔介质叠后地震烃类检测方法，属于地球物理方法在储层预测及油气检测中的应用范畴。
背景技术
烃类检测技术在油气勘探中具有重要作用，目前国内外学者对于烃类检测的研究主要包括以下几个方面，如叠前AVO油气预测技术、瞬时谱分析探测低频阴影油气检测技术、单频成像油气预测技术、多波多分量岩石物理参数油气预测技术、利用地震子波数据体进行油气预测技术、通过地震数据结构分析进行油气预测技术、多相介质油气检测技术等。根据应用的数据类型总体可以归纳为叠前和叠后两类，前者理论基础是AVO技术，需要全波列测井资料，叠前道集保幅性直接影响预测精度；后者更多的是根据地震波通过含烃储层前后地震波场即传播能量的变化寻找油气存在的信息，分析数据量相对较少，运算速度快，适用性更广。
地震信号的低频成分和储层流体具有很强的相关性，根据BIOT双相介质理论，若多孔介质的孔隙单元相互连通，地震波在含流体多孔介质中传播时，流体和固体的振动相互作用、藕合，致使孔隙中流体在从孔隙空间中产生流动（比如从一个孔隙向周围较大的孔隙空间流动），从而引起流体和固体颗粒的相对运动，导致地震波的振幅衰减。而这种振幅衰减为进行烃类检测创造了客观条件，但该理论如何应用于碳酸盐岩薄气藏预测是一个非常值得探索的方向。
发明内容
本发明的目的在于克服现有技术存在的上述问题，提供一种饱和多孔介质叠后地震烃类检测方法。本发明将一个针对饱和多孔介质的实验室基础理论属性表达式转换成易实现的叠后地震数据烃类检测方法，且利用子波分解重构生成的敏感频段地震体所求的流体活动属性相较于频谱分解的单频体稳定性更好。
为实现上述目的，本发明采用的技术方案如下：
一种饱和多孔介质叠后地震烃类检测方法，其特征在于：针对目的层段寻找含油气井和非含油气井出现振幅异常的敏感频率段，再提取新的敏感频带地震体，在振幅均一化的基础上进行差值运算，得到流体活动属性体。
所述方法具体包括如下步骤：
a、输入原始地震数据及目的层位数据，并根据储层厚度确定分析时窗；
b、在分析时窗中选取含油气井和非含油气井进行频谱特征差异性分析，确定振幅差异对应的敏感频率段；
c、应用子波分解重构方法提取新的敏感频带地震体；
d、对重构的敏感频带地震体进行振幅均一化后的差值计算得到流体活动属性体；
e、根据计算得到流体活动属性体，对目的层段提取平面属性，根据提取的平面属性值得出检测结果。
所述步骤d中，通过如下公式计算得到流体活动属性体：
式中：
F——流体活动属性体（无量纲）；
A——瞬时振幅；
——地震波角频率。
所述步骤e中，平面属性值为0.5—1表征含油气区域。
采用本发明的优点在于：
1、本发明对流体活动属性公式进行了进一步的推导，并在其在实际应用中的方法提出了改进，并在求取流体活动属性时利用子波分解重构生成的敏感频段地震体代替频谱分解的单频体，将一个实验室公式推导为一种具有理论基础且易于实现的叠后地震烃类检测工具，提高了该方法在烃类检测中的有效性和稳定性。
2、本发明在多个区块取得了良好效果。以苏里格气田某区块马家沟组岩溶风化壳气藏为例，该气藏储层厚度薄，仅1~3米，储集空间以次生溶蚀孔洞和裂缝为主，主要发育在距离马家沟组0~10米层段内，上覆砂泥岩和煤层与储层段碳酸盐岩形成一套强反射，储层段和非储层在常规地震剖面上差异不明显，烃类预测难度大。在本发明基础上对该区块南部主力产层马五1 3段16口井的试气结果和流体预测属性预测情况相比对，共有14口井吻合，吻合度达到87.5%，说明该方法的适用性和稳定性得到了极大提升，适于推广应用。
附图说明
图1为本发明研究区含油气井和干井的目的层段频谱特征差异性分析图；
图2为由含油气层段厚度决定分析时窗的大小示意图；
图3a 为10~15HZ重构剖面及目的层位图；
图3b 为15~20HZ重构剖面及目的层位图；
图4为利用能量均一化后的振幅差求取的频谱能量衰减斜度剖面图。
具体实施方式
实施例1
本发明针对研究区目的层段应用频谱特征分析技术寻找含油气井和非含油气井出现振幅异常的敏感频率段，再利用子波分解重构技术提取新的敏感频带地震体，在振幅均一化的基础上进行差值运算，得到流体活动属性体。
所述方法具体包括如下步骤：
a、输入原始地震数据及目的层位数据，并根据储层厚度确定分析时窗；
b、在分析时窗中选取含油气井和非含油气井进行频谱特征差异性分析，确定振幅差异对应的敏感频率段；
c、应用子波分解重构方法提取新的敏感频带地震体；
d、对重构的敏感频带地震体进行振幅均一化后的差值计算得到流体活动属性体；
e、根据计算得到流体活动属性体，对目的层段提取平面属性，根据提取的平面属性值得出检测结果。
所述步骤d中，通过如下公式计算得到流体活动属性体：
式中：
F——流体活动属性体（无量纲）；
A——瞬时振幅；
——地震波角频率。
所述步骤e中，平面属性值为0.5—1表征含油气区域。
实施例2
本实施例结合附图对本发明做进一步说明。
图1为研究区含油气井和干井的目的层段频谱特征差异性分析图。可以看到原始地震资料目的层段主频约为30HZ，低频段10HZ~20HZ振幅频谱变化斜率最大，能量衰减最为明显，对应敏感频率带。
图2为输入的原始地震剖面和目的层位，由含油气层段厚度决定分析时窗的大小。
图3a和图3b为利用子波分解与重构技术提取的敏感频率段重构剖面，即对应15~20HZ频段地震体，对应10~15HZ频段地震体A（10~15HZ）。
图4为利用能量均一化后的和A（10~15HZ）振幅差求取的频谱能量衰减斜度剖面，即流体活动属性剖面，可以发现含油气层段和干层在流体活动性剖面上表现差异明显，含油气层段对应流体活动性强区。
实施例3
本实施例对本发明中流体活动属性体计算公式的推导得出进行说明。
在流体和岩石特性在合理范围时，无量纲参数在低频（<1KHz）时很小，弹性介质和饱和流体多孔介质分界面上的地震反射系数R与地震波角频率ω存在以下关系：
(1)
其中，和是实系数，表征岩石和流体的力学特性；，是储层渗透率，是流体粘滞系数，是储层岩石体密度，当=0时，反射系数绝对值取得最大值。
将(1)式对频率求导，得到：
(2)
令，它是孔隙流体和岩石骨架的弹性性质的复函数；
设储层流体活动属性，即储层渗透率和流体粘滞系数之比；
则(2)式可表示为：
(3)
则流体活动属性F可推导为：
(4)
从上式可见，储层流体活动属性F和地震角频率为时的反射系数的一阶导数成正比，对地震信号进行时频分解后，其单频瞬时谱振幅能准确刻画该频率对应的地震反射能量，在实际计算中，可用瞬时谱振幅A(ω)代替相应频率处的反射系数R，则有：
(5)
低频域饱和流体多孔介质储集层中流体的活动性近似与储集层渗透率、流体密度与流体黏度比值的函数成正比，且与地震反射振幅对地震反射频率偏导的绝对值成正比。
对式(5)右侧进行推导，流体活动性和低频域的两个频率体振幅差值梯度成正比：
(6)。
Claims (4)
1.一种饱和多孔介质叠后地震烃类检测方法，其特征在于：针对目的层段寻找含油气井和非含油气井出现振幅异常的敏感频率段，再提取新的敏感频带地震体，在振幅均一化的基础上进行差值运算，得到流体活动属性体。
2.根据权利要求1所述的饱和多孔介质叠后地震烃类检测方法，其特征在于：所述方法具体包括如下步骤：
a、输入原始地震数据及目的层位数据，并根据储层厚度确定分析时窗；
b、在分析时窗中选取含油气井和非含油气井进行频谱特征差异性分析，确定振幅差异对应的敏感频率段；
c、应用子波分解重构方法提取新的敏感频带地震体；
d、对重构的敏感频带地震体进行振幅均一化后的差值计算得到流体活动属性体；
e、根据计算得到流体活动属性体，对目的层段提取平面属性，根据提取的平面属性值得出检测结果。
3.根据权利要求2所述的饱和多孔介质叠后地震烃类检测方法，其特征在于：所述步骤d中，通过如下公式计算得到流体活动属性体：
式中：
F——流体活动属性体（无量纲）；
A——瞬时振幅；
——地震波角频率。
4.根据权利要求2或3所述的饱和多孔介质叠后地震烃类检测方法，其特征在于：所述步骤e中，平面属性值为0.5—1表征含油气区域。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810283307.5A CN108549103A (zh) | 2018-04-02 | 2018-04-02 | 饱和多孔介质叠后地震烃类检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810283307.5A CN108549103A (zh) | 2018-04-02 | 2018-04-02 | 饱和多孔介质叠后地震烃类检测方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108549103A true CN108549103A (zh) | 2018-09-18 |
Family
ID=63513741
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810283307.5A Pending CN108549103A (zh) | 2018-04-02 | 2018-04-02 | 饱和多孔介质叠后地震烃类检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108549103A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109581498A (zh) * | 2018-11-08 | 2019-04-05 | 中国石油天然气集团有限公司 | 储层厚度分布确定方法、系统、设备、可读介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106483564A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 一种利用地震低频信息进行流体识别的方法 |
| CN107132573A (zh) * | 2017-02-28 | 2017-09-05 | 中国石油化工股份有限公司 | 一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法 |
| US20170292377A1 (en) * | 2016-04-11 | 2017-10-12 | Schlumberger Technolgoy Corporation | Methods for determining saturation parameters of a formation |
| CN107450101A (zh) * | 2017-07-21 | 2017-12-08 | 杰奥世博（北京）技术有限公司 | 基于多子波的地震数据分频方法 |
-
2018
- 2018-04-02 CN CN201810283307.5A patent/CN108549103A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106483564A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 一种利用地震低频信息进行流体识别的方法 |
| US20170292377A1 (en) * | 2016-04-11 | 2017-10-12 | Schlumberger Technolgoy Corporation | Methods for determining saturation parameters of a formation |
| CN107132573A (zh) * | 2017-02-28 | 2017-09-05 | 中国石油化工股份有限公司 | 一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法 |
| CN107450101A (zh) * | 2017-07-21 | 2017-12-08 | 杰奥世博（北京）技术有限公司 | 基于多子波的地震数据分频方法 |
Non-Patent Citations (2)
| Title |
|---|
| 杨静 等: "GeoEast流体活动性属性技术在S井区储层描述中的应用", 《石油地球物理勘探》 * |
| 汲生珍 等: "子波分解与重构技术在储层预测中的应用", 《石油天然气学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109581498A (zh) * | 2018-11-08 | 2019-04-05 | 中国石油天然气集团有限公司 | 储层厚度分布确定方法、系统、设备、可读介质 |
| CN109581498B (zh) * | 2018-11-08 | 2020-07-10 | 中国石油天然气集团有限公司 | 储层厚度分布确定方法、系统、设备、可读介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102809762B (zh) | 基于全频带地震信息挖掘的储层成像技术 | |
| CN103728659B (zh) | 一种提高地下岩溶探测精度的方法 | |
| CN104142519B (zh) | 一种泥岩裂缝油藏预测方法 | |
| CN110348135B (zh) | 一种随钻声波测井评价地层渗透率的方法 | |
| CN112746837A (zh) | 基于分布式光纤传感的页岩油藏勘探数据采集系统及方法 | |
| CN105134189A (zh) | 一种LogGMIR测井地应力识别优质储层的方法 | |
| CN106526670A (zh) | 一种碎屑岩储层中地震属性砂体空间分布描述及评价的方法 | |
| CN105093313B (zh) | 一种岩溶型油藏单井油气产能预测方法及装置 | |
| CN103424776A (zh) | 一种碳酸盐岩油气储层裂缝地震检测方法 | |
| CN111399056B (zh) | 一种基于划分方位角滤波预测裂缝强度的方法 | |
| CN103235339A (zh) | 一种时频分解地震流体识别方法 | |
| CN110456412B (zh) | 一种基于叠后地震数据识别碳酸盐岩储层流体饱和度的方法 | |
| CN104007462B (zh) | 基于衰减各向异性的裂缝预测方法 | |
| CN105931125B (zh) | 一种致密油分段多簇体积压裂水平井产量预测方法 | |
| CN103412327A (zh) | 一种裂缝性储层的粘弹性参数提取方法 | |
| CN111060986B (zh) | 一种地层压力预测方法及岩性油藏评价方法 | |
| CN214576965U (zh) | 基于分布式光纤传感的页岩油藏勘探数据采集系统 | |
| CN101852864B (zh) | 一种利用地表一致性统计频谱分析技术处理海量地震数据的方法 | |
| CN104880730A (zh) | 基于Synchrosqueezing变换的地震资料时频分析和衰减估计方法 | |
| CN103777245A (zh) | 基于地震资料的油气成藏条件定量评价方法 | |
| CN108549103A (zh) | 饱和多孔介质叠后地震烃类检测方法 | |
| CN117452491A (zh) | 一种用于复杂山地地表条件下煤系气藏特征识别的组合勘探方法 | |
| CN116027416A (zh) | 一种砂砾岩体岩性圈闭识别方法 | |
| CN104007465A (zh) | 基于弹性反演的纯纵波数据刻画储层发育细节方法 | |
| Li et al. | Predicting the distribution of thin bed reservoirs by broad frequency band seismic |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180918 |