一种自适应方向增强的地质不连续异常检测方法

本发明涉及石油地震勘探领域，具体是一种自适应方向增强的地质不连续异常检测方法，该方法在压制地震数据特定空间方向噪声的同时，能够提取和突出储层结构及地质构造在特定方向上的不连续性异常信息，为储层预测与地质构造解释提供支持。

背景技术：

1、地震资料中包含丰富的具有各种方向特征的地质异常信息，提取这些地质异常及其方向信息是油气地球物理勘探的重要内容之一。然而，这些有效信息常受到地震资料采集过程中所含噪声的影响，因此，压制噪声的同时，保留地震资料中的地震异常信息及其方向特征十分重要。传统的基于高斯滤波的噪声压制方法，在数据各个方向上的滤波权重相同，缺乏在特定方向上压制噪声，突出相应方向有效信息的能力。针对这一问题，不少研究者做了大量工作，freeman(1991)提出了一种方向自适应滤波设计方法，这种方向滤波器由一个方向滤波器基集的元素通过插值得到，但却未能指出如何正确得到这个滤波器基集和正确的插值规则；van ginkel等(1997)提出去卷积的方法来提高高斯滤波的角分辨率，但是该方法用到计算复杂的傅立叶去卷积算法。geusebroek等(2003)提出了一种快速各向异性高斯滤波方法，这种滤波器通过在两个不同方向上选用不同的高斯尺度，从而使滤波器可以在去噪时可以较好保留图像边缘等重要信息，并利用高斯函数的可分解性将滤波器沿长轴和短轴方向分解为两个一维滤波器与图像卷积，使计算简化。

2、地震曲率属性是继相干技术后的又一种强大的地震几何属性分析方法。20世纪90年代中期以来，以二阶导数为基础的曲率属性在构造解释上得到了迅速发展和广泛应用，如检测地质的张开裂缝与高斯曲率的关系(lisle，1994)；ericsson等人(1998)论证了油气产量与曲率之间的关系，认为曲率非常有助于常规构造和地形解释、改善断层的成像质量；roberts(2001)详述了曲率属性的基本理论，提出了第一代曲率分析方法——层面曲率属性(surface curvature attribute)的计算和工作流程，表明曲率属性对提取断层和裂缝走向等构造几何特征十分有效，为曲率属性在地震资料构造解释中的推广和应用奠定了基础；hart(2002)研究了新墨西哥州西北部地层的走向曲率与张开裂缝有密切的关系；bergbauer等人(2003)采用kx-ky滤波计算了不同波长的曲率；al-dossary和marfurt(2006)实现了基于三维地震数据体的第二代曲率分析方法——体曲率属性(volumetriccurvature attribute)，它同时利用了地震资料的振幅信息和地震反射的几何特征信息，可以减少噪声和层位追踪的闭合差对曲率的影响。

3、地震数据去噪时，如果不能较好地保留地震异常信息及其方向特征，不利于提高储层地震检测与构造解释的精度。如果在地震数据去噪处理中，能根据解释的需要，保留和加强特定方向的地质异常、并在后续的地震几何属性分析中进一步提取和突出相应方向的不连续性异常信息，则能更精确可靠地检测储层结构及地质构造信息(如不同方向断层、裂缝带及其密度等指标)，这正是本发明专利的方法实现的工作。

技术实现思路

1、本发明的目的是提供一种自适应方向增强的地质不连续异常信息的检测方法，本发明的原理是根据最大方位角及裂缝走向，不断的更换各向异性高斯滤波器的方向，使得数据所受的异常干扰小的同时，保留的有效信息多。本发明方法能有效地压制与裂缝走向不同的干扰，突出地震数据的有效信息，实现特定方向的不连续信息的增强及高精度地质异常信息的提取，并可以根据储层与构造解释的需要，保留、提取和突出一个或者多个特定方向的不连续性异常信息，在提高地震资料信噪比的同时，能够更加准确可靠地检测地下断裂系统及储层的空间分布和内部结构。本发明方法包括以下主要步骤：

2、(1)输入三维地震数据集，采样周期为δt，设置窗长为2w+1个采样点，对输入地震数据集进行分时窗处理，设在时窗位置τ处，构建相应的时窗地震数据集s(τ+mδt,x,y)|m∈[-w,w]，m为窗采样点序号，(x,y)为平面坐标；

3、(2)利用如下各向异性高斯滤波器g：

4、

5、式中，

6、

7、(σx,σy)为各向异性高斯滤波器的调节参数，θ为方位角；在0°～180°的方位角范围内，选取k个方位角[θ1,θ2,θ3,...,θk]，对时窗位置τ处的时窗地震数据集s(τ+mδt,x,y)|m∈[-w,w]，如下式进行各个方位角的各向异性高斯滤波处理：

8、f(τ+mδt,x,y；θk)|m∈[-w,w],k∈[1,k]＝gθ(x,y；σx,σy,θk)s(τ+mδt,x,y)|m∈[-w,w]；

9、(3)拾取f(τ+mδt,x,y；θk)|m∈[-w,w],k∈[1,k]的最大值，并将其对应的方位角记为最大方位角θm，表示如下：

10、

11、(4)按下式计算时窗位置τ处的层段的最大方位角滤波数据z：

12、

13、(5)对得到的数据z(τ,x,y)进行l级量化，构成最大方位角滤波量化数据体zl(τ,x,y)，利用zl(τ,x,y)计算k个方位角[θ1,θ2,θ3,...,θk]的纹理参数：能量e(θk)，相关性c(θk)；

14、(6)利用全部方位角的纹理参数e(θk)和c(θk)，按照下式计算能量与相关性的权重系数，以权重系数最大时的方位角，作为不连续地质异常信息的最优方向，记为θo：

15、

16、(7)对数据z(τ,x,y)进行方位角为θo的各向异性高斯滤波处理，如下式：

17、fo(τ,x,y)＝gθ(x,y；σx,σy,θo)z(τ,x,y)；

18、(8)利用处理后的数据fo(τ,x,y)，按照下式计算，得到自适应方向增强的地质不连续异常信息检测结果：

19、

20、其中：

21、

技术特征：

1.一种自适应方向增强的地质不连续异常检测方法，包括以下步骤：

技术总结
本发明涉及石油地震勘探数据处理与解释领域，是一种通过地震信号的各向异性高斯迭代滤波、纹理参数计算、方向增强的不连续信息检测等处理，确定地质不连续信息的方向，实现自适应方向增强的地质不连续异常检测方法。它在压制地震数据中存在的不同空间方向噪声和背景干扰的同时，还能突出该各个方向的地质不连续异常和有效的地质构造特征，能更好地提取储层结构及地质构造的不连续性异常信息，可为精确的地质构造解释、断裂系统分布的高精度描述、储层预测等提供可靠的技术支持。

技术研发人员：陈学华,黎康毅,吕丙南,吴昊杰,王聪,王胜义,赵晨斐
受保护的技术使用者：成都理工大学
技术研发日：
技术公布日：2024/1/13