CN104570078B - 一种基于频率域倾角的相似性横向变化率的溶洞检测方法 - Google Patents
一种基于频率域倾角的相似性横向变化率的溶洞检测方法 Download PDFInfo
- Publication number
- CN104570078B CN104570078B CN201310491952.3A CN201310491952A CN104570078B CN 104570078 B CN104570078 B CN 104570078B CN 201310491952 A CN201310491952 A CN 201310491952A CN 104570078 B CN104570078 B CN 104570078B
- Authority
- CN
- China
- Prior art keywords
- mrow
- frequency domain
- seismic
- inclination angle
- similitude
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于频率域倾角的相似性横向变化率的溶洞检测方法,属于石油地球物理勘探与开发领域。本方法包括:第一步,对地震资料进行Hilbert变换,将时间域的地震资料转换至复数域,得到Hilbert地震道,然后计算Hilbert地震道的瞬时相位和频率域倾角;第二步,对第一步得到的频率域倾角进行相似性检测,寻找连续性较差的异常点,即利用地震波场信息计算各道之间的相关性,突出不相关的异常边界,据此识别细微的岩层横向不均一性和断裂特征,使得裂缝和溶洞的响应特征得到放大;第三步,开展频率域倾角的相似性横向变化的边缘检测,从而确定了溶洞型储层的发育区。
Description
技术领域
本发明属于石油地球物理勘探与开发领域,具体涉及一种基于频率域倾角的相似性横向变化率的溶洞检测方法,利用频率域倾角属性结果开展相似性横向变化率的边缘检测扫描,可应用于油田勘探中碳酸盐岩储层中溶洞的检测和识别。
背景技术
相似性横向变化率即地震道的相似性,目前较流行的算法为基于互相关和基于相似的算法,主要缺陷一是计算的只是邻近道与中心道的相关关系,不能考虑相邻多道间彼此在结构上的相互关系;二是不同方向计算的相关系数组合后有平均效应,计算道数越多平均效应越严重,降低了空间分辨能力及异常的连续性。
基于频率域倾角的相似性横向变化率技术,主要是通过Hilbert变换将时间域地震数据转换到频率域,再进行倾角扫描和相似性横向变换的计算,最后对计算再进行边缘检测分析,从而便于平面上溶洞的识别。模型正演和实际应用效果表明,开展基于频率域倾角的相似性横向变化率,提高了地震道的相似性对于溶洞型异常体的检测能力。
相似性横向变化率即地震道的相似性(相干),目前技术发展较为成熟,通过多年来在多个工区的的应用表明,该方法对于检测地震剖面上同向轴的异常点有一定的效果,但是当对应地质体为不同尺度裂缝和溶洞同时发育时,该检测方法能反映裂缝和溶洞的整体发育状况,但是无法区分出裂缝与溶洞的发育特征。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种基于频率域倾角的相似性横向变化率的溶洞检测方法,采用沿反射曲面的面主测线或横测线方向进行频率域倾角属性计算的方法,提高检测精度的同时,放大地震剖面上裂缝和溶洞产生的异常点,在此基础上开展相似性横向变化率的计算,使得溶洞的异常被区分显示出来。
本发明是通过以下技术方案实现的:
一种基于频率域倾角的相似性横向变化率的溶洞检测方法,包括:
第一步,对地震资料进行Hilbert变换,将时间域的地震资料转换至复数域,得到Hilbert地震道,然后计算Hilbert地震道的瞬时相位和频率域倾角;
第二步,对第一步得到的频率域倾角进行相似性检测,寻找连续性较差的异常点,即利用地震波场信息计算各道之间的相关性,突出不相关的异常边界,据此识别细微的岩层横向不均一性和断裂特征,使得裂缝和溶洞的响应特征得到放大;
第三步,开展频率域倾角的相似性横向变化的边缘检测,从而确定了溶洞型储层的发育区。
所述第一步中计算Hilbert地震道的瞬时相位和频率域倾角是这样实现的:
利用公式(1)和公式(2)计算瞬时相位和频率域倾角:
瞬时相位
频率域倾角
其中f(D,t)是Hilbert地震道的实部,g(D,t)是Hilbert地震道的虚部,Ph(D,t)是瞬时相位,D是Hilbert地震道;x和y分别指三维地震数据的主测线和联络线方向。
所述第二步是通过计算公式(3)实现的:
通过公式(3)计算多个相邻道的相似度:
其中,S(t,τ)为t时刻的相似度值,是t时刻不同时间延迟的瞬时相似度值s(t,τ)之和;N为计算相似度的时窗;M为计算相似度所采用的相邻地震道数目;fm(t)为第m个地震道的频率域倾角数据,由公式(2)计算得到;τ为时间延迟。公式(3)中的fm(t)就是公式(2)中的Dip。
所述N的大小根据地震频带确定,选取40-80ms,频带越宽,时窗越短。
所述第三步是通过计算公式(4)实现的:
其中,D’为叠后地震属性体数据,该数据体现了溶洞型储层的发育区,F为利用频率域倾角计算得到的相似度值,即式(3)的计算结果。公式(4)中的F就是公式(3)中的
与现有技术相比,本发明的有益效果是:
1)相对于传统的频率域倾角属性,本发明又进一步的开展了相似性横向变化率扫描,剔除了裂缝对振幅连续性的影响,更好的凸显了碳酸盐岩溶洞的异常响应。
2)从模型和实际数据的结果上可以看出,该方法有较好的计算结果,能凸显碳酸盐岩溶洞的异常响应,同时为确定碳酸盐岩溶洞型储层的发育区提供更有利的预测手段。
附图说明
图1是Hilbert变换。
图2是频率域瞬时平均倾角。
图3是基于频率域倾角的相似性横向变化的边缘检测,对于溶洞型储集体的展布刻画清晰。
图4是试验区地震剖面,溶洞储集体特征明显。
图5是本发明方法的步骤框图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
基于频率域扫描倾角属性的相似性横向变化率技术,首先在频率域进行倾角扫描,然后通过高精度相似性检测的方法寻找连续性较差的异常点,在此基础上,开展相似性横向变化率的边缘检测,剔除裂缝引起异常点,从而确定溶洞型储层的发育区。
如图5所示,其具体实现的步骤为:
(1)频率域倾角扫描
频率域地震倾角属性是利用瞬时相位谱来计算地震倾角的属性。首先对地震数据进行Hilbert变换,得到Hilbert地震道,使得地震道信息由时间域变换至频率域,然后对Hilbert地震道进行瞬时相位计算,最后通过瞬时相位在空间上的一阶导数,从而得到频率域倾角的扫描结果(如公式1、2所示)。
瞬时相位
频率域倾角
其中x和y分别指三维地震数据的主测线和联络线方向,f(D,t)是Hilbert地震道的实部,g(D,t)是Hilbert地震道的虚部,Ph(D,t)是瞬时相位,D是Hilbert地震道;即
瞬时相位独立于振幅,不受其影响,仅与地震波传播相位有关。频率域倾角是计算瞬时相位在空间上的一阶导数,可以反映地层或反射界面的倾角变化,可以称之为瞬时倾角。
(2)相似性检测
相似度的计算是由计算两个相邻道之间的相关系数发展而来的,可以计算多个相邻道的相似程度,其定义为:
式中S(t,τ)为t时刻的相似度值,是t时刻不同时间延迟的瞬时相似度值s(t,τ)之和;N为计算相似度的时窗,其大小根据地震频带而定,一般选取40-80ms,频带越宽,时窗越短;M为计算相似度所采用的相邻地震道数目,选取时要考虑地震信噪比、反射界面弯曲程度和横向连续性等,对于信噪比高、界面曲率小、连续性好的地震数据体,M值可取小一些;fm(t)为第m个地震道集数据;τ为时间延迟。相似度能反映各地震道之间的相似程度,相似度越大说明地层的连续性越好。
(3)边缘异常分析
地震的边缘检测技术,是用来检测地震剖面上连续性差的异常点。通过对频率域倾角开展高精度相似性检测(完成步骤(2)的相似性就完成了相似性检测),使得溶洞的响应特征得到放大,在此基础上,开展频率域倾角的相似性横向变化的边缘检测,从而确定了溶洞的发育位置。物理正演模型和实际资料的应用结果表明,该方法能有效的削弱了裂缝对溶洞型储层的影响,从而明确溶洞的发育位置,计算公式如下:
其中F为利用频率域倾角计算得到的相似度值,即式(3)的计算结果。
公式(4)得到的是一种叠后地震属性体数据,属于地震几何属性的一种,与以往不同的是,通常这里的F都选用地震道振幅谱,而在本方法里的F采用的是经过相似性计算的频率域倾角,它放大了溶洞在地震振幅谱上的响应特征,使得计算得到的地震属性体上的溶洞信息更加清晰易辨认。
基于频率域倾角的相似性横向变化率技术,采用沿反射曲面的面主测线或横测线方向进行频率域倾角属性计算的方法,提高了检测精度的同时,放大了地震剖面上裂缝和溶洞产生的异常点,在此基础上开展相似性横向变化率的计算,使得溶洞的异常被区分显示出来。
图1为塔河油田某试验区碳酸盐岩储层段时窗内的Hilbert变换平面展布图,通过Hilbert变换,地震振幅体从时间域转换到了复数域,然后对复数域地震数据开展角度扫描,得到图2,即图1对应时窗范围的频率域瞬时平均倾角,最后通过高精度相似性边缘检测的方法,寻找连续性较差的异常点的同时,压制了裂缝引起异常反应,从而确定溶洞型储层的发育区;图4为试验区过W1、W2、W3井的地震剖面,其中W1井为溶洞型储层钻遇井,W2井为裂缝型储层钻遇井,W3为非溶洞型储层钻遇井,W1和W2井在地震剖面上响应特征明显,W3井无明显响应,这与图3中对溶洞型储层的平面预测结果对应井点位置的响应一致,通过多口钻遇目的层段的井资料的分析发现,碳酸盐岩储层段的溶洞型储层发育区与图3中预测的溶洞平面展布有较高的吻合度。
实施例中的具体步骤如下:
第一步,首先是对地震资料进行Hilbert变换,将时间域地震资料转换至复数域,然后再进行横向、垂向求导,并计算瞬时倾角。
第二步,对倾角扫描结果开展高精度相似性检测,寻找连续性较差的异常点(通过公式(3)的计算实现),利用地震波场信息计算各道之间的相关性,突出不相关的异常边界(计算得到相关性的值就突出了不相关的异常边界),据此识别细微的岩层横向不均一性和断裂特征,使得裂缝和溶洞的响应特征得到放大。;
第三步,在此基础上,开展频率域倾角的相似性横向变化的边缘检测,从而确定了溶洞型储层的发育区,有效的削弱了裂缝对溶洞型储层的影响。
本发明方法目前只在碳酸盐岩储层的模型和实际数据中开展了测试,碎屑岩等其它类型储层未进一步研究。
本发明属于石油地球物理勘探与开发技术方法,主要是针对碳酸盐岩型储集体应用地震资料进行溶洞型储层检测的方法,具体方法技术首先是对地震资料在频率域进行倾角扫描,然后通过高精度相似性检测的方法,寻找连续性较差的异常点,在此基础上,开展相似性横向变化率的边缘检测,剔除裂缝引起异常点,从而确定溶洞型储层的发育区,通过分析发现预测结果与钻井吻合程度较高。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (3)
1.一种基于频率域倾角的相似性横向变化率的溶洞检测方法,其特征在于:所述方法包括:
第一步,对地震资料进行Hilbert变换,将时间域的地震资料转换至复数域,得到Hilbert地震道,然后计算Hilbert地震道的瞬时相位和频率域倾角;
第二步,对第一步得到的频率域倾角进行相似性检测,寻找连续性较差的异常点,即利用地震波场信息计算各道之间的相关性,突出不相关的异常边界,据此识别细微的岩层横向不均一性和断裂特征,使得裂缝和溶洞的响应特征得到放大;
第三步,开展频率域倾角的相似性横向变化的边缘检测,从而确定了溶洞型储层的发育区;
所述第一步中计算Hilbert地震道的瞬时相位和频率域倾角是这样实现的:
利用公式(1)和公式(2)计算瞬时相位和频率域倾角:
瞬时相位
频率域倾角
其中f(D,t)是Hilbert地震道的实部,g(D,t)是Hilbert地震道的虚部,Ph(D,t)是瞬时相位,D是所述Hilbert地震道;x和y分别指三维地震数据的主测线和联络线方向;
所述第二步是通过计算公式(3)实现的:
通过公式(3)计算多个相邻道的相似度:
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&tau;</mi>
<mo>=</mo>
<mo>-</mo>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msubsup>
<mi>f</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msubsup>
<mi>f</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&tau;</mi>
<mo>=</mo>
<mo>-</mo>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,s(t,τ)为t时刻的相似度值,是t时刻不同时间延迟的瞬时相似度值s(t,τ)之和;N为计算相似度的时窗;M为计算相似度所采用的相邻地震道数目;fm(t)为第m个地震道的频率域倾角数据,由公式(2)计算得到;τ为时间延迟。
2.根据权利要求1所述的基于频率域倾角的相似性横向变化率的溶洞检测方法,其特征在于:所述N的大小根据地震频带确定,选取40-80ms,频带越宽,时窗越短。
3.根据权利要求1所述的基于频率域倾角的相似性横向变化率的溶洞检测方法,其特征在于:所述第三步是通过计算公式(4)实现的:
<mrow>
<msup>
<mi>D</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>F</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>F</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,D’为叠后地震属性体数据,该数据体现了溶洞型储层的发育区,F为利用频率域倾角计算得到的相似度值,即式(3)的计算结果,即:
<mrow>
<mi>F</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&tau;</mi>
<mo>=</mo>
<mo>-</mo>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
1
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310491952.3A CN104570078B (zh) | 2013-10-18 | 2013-10-18 | 一种基于频率域倾角的相似性横向变化率的溶洞检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310491952.3A CN104570078B (zh) | 2013-10-18 | 2013-10-18 | 一种基于频率域倾角的相似性横向变化率的溶洞检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104570078A CN104570078A (zh) | 2015-04-29 |
CN104570078B true CN104570078B (zh) | 2017-10-27 |
Family
ID=53086592
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310491952.3A Active CN104570078B (zh) | 2013-10-18 | 2013-10-18 | 一种基于频率域倾角的相似性横向变化率的溶洞检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104570078B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105093300B (zh) * | 2015-07-27 | 2017-10-17 | 中国石油天然气股份有限公司 | 一种地质体边界识别方法及装置 |
CN107229068A (zh) * | 2016-03-24 | 2017-10-03 | 中国石油化工股份有限公司 | 用于识别勘探地球物理信号边界的方法和装置 |
CN106199710B (zh) * | 2016-06-29 | 2018-05-18 | 中国石油化工股份有限公司 | 基于混合倾角扫描振幅变化率的潜山储层地震识别方法 |
CN108008455B (zh) * | 2017-11-29 | 2018-12-25 | 广州市公路勘察设计有限公司 | 桥梁基础溶洞勘测方法和装置 |
CN112630831B (zh) * | 2019-10-08 | 2024-04-09 | 中国石油化工股份有限公司 | 碳酸盐岩溶洞纵向尺度计算方法及系统 |
CN114428321B (zh) * | 2020-10-28 | 2024-09-20 | 中国石油化工股份有限公司 | 一种缝洞型储层识别方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EG24691A (en) * | 2007-01-17 | 2010-05-17 | Pgs Geophysical As | Diagonal gather trace interpolation |
CN102012521B (zh) * | 2010-10-27 | 2013-03-20 | 中国石油化工股份有限公司 | 一种地震储层预测中叠前裂缝的检测方法 |
-
2013
- 2013-10-18 CN CN201310491952.3A patent/CN104570078B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104570078A (zh) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104570078B (zh) | 一种基于频率域倾角的相似性横向变化率的溶洞检测方法 | |
US9645268B2 (en) | Seismic orthogonal decomposition attribute | |
CN106094029B (zh) | 利用偏移距矢量片地震数据预测储层裂缝的方法 | |
CN103149589B (zh) | 火成岩油气勘探方法及装置 | |
CN103645503B (zh) | 一种三维时间域照明分析及振幅补偿方法 | |
Song et al. | Multi-waveform classification for seismic facies analysis | |
CN103733089B (zh) | 用于包括不确定性估计的地下表征的系统和方法 | |
Zimmer | Microseismic design studies | |
US20160356909A1 (en) | Determining p-wave azimuthal anisotropy from walkaround vsp with offset dependent slowness corrections | |
CN105005080A (zh) | 一种利用振幅比率属性识别地层圈闭尖灭线的方法 | |
CN105425299B (zh) | 确定地层裂缝分布的方法和装置 | |
CN112305591B (zh) | 隧道超前地质预报方法、计算机可读存储介质 | |
CN104950327A (zh) | 确定地面微地震观测系统的检波器的位置的方法 | |
CN104216013B (zh) | 基于宽方位角资料的c3相干体的方法 | |
CN110320576A (zh) | 超前地质预报方法、装置及存储介质 | |
CN103576191A (zh) | 一种采用地震属性识别断层的方法 | |
CN102565852B (zh) | 针对储层含油气性检测的角度域叠前偏移数据处理方法 | |
CN104199088A (zh) | 一种提取入射角道集的方法及系统 | |
CN106597547A (zh) | 薄储层地震准确描述方法 | |
CN112946743B (zh) | 区分储层类型的方法 | |
CN104459787A (zh) | 一种垂直接收阵列地震记录的速度分析方法 | |
CN105259578A (zh) | 基于垂直地震剖面资料的地层速度确定方法 | |
CN107339099B (zh) | 一种确定储层岩性的方法和装置 | |
CN105093300A (zh) | 一种地质体边界识别方法及装置 | |
CN110703329B (zh) | 基于弱振幅地震反射形成机制的岩性油藏边界确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |