基于优化Morlet小波的分频扩散滤波断层强化方法技术

本发明专利技术提供一种基于优化Morlet小波的分频扩散滤波断层强化方法，包括：将叠后地震数据转换为MATLAB可识别的数据格式；对地震数据进行频谱分析并提取子波；改进Morlet子波母函数，构建Morlet子波函数组合；利用优化Morlet子波进行小波分频处理；计算叠后地震数据的结构张量，进而构建基于结构张量迹的新断层因子；在结构张量迹的保护下，对分频地震数据进行各项异性扩散滤波；对分频扩散滤波结果进行小波反变换，并重构得到提高了信噪比和分辨率，并保护了断层信息的地震数据。该方法实现了地震数据去噪保边的同时提高其分辨率，从而提高地震数据低序级断层识别能力。

Fault Enhancement Method Based on Frequency Division Diffusion Filter Based on Optimized Morlet Wavelet

The invention provides a method of fault enhancement based on optimized Morlet wavelet by frequency division diffusion filtering, which includes: converting post-stack seismic data into identifiable data format of MATLAB; spectrum analysis of seismic data and extraction of wavelet; improving Morlet wavelet generating function to construct Morlet wavelet function combination; using optimized Morlet wavelet to perform wavelet frequency division processing; and calculating post-stack seismic data. Structural tensor is used to construct a new fault factor based on structural tensor trace. Under the protection of structural tensor trace, anisotropic diffusion filtering is applied to the frequency-division seismic data. The results of the frequency-division diffusion filtering are inversely transformed by wavelet transform, and the seismic data with improved signal-to-noise ratio and resolution are reconstructed, and the fault information is protected. This method realizes the denoising and edge preservation of seismic data and improves its resolution, thus improving the recognition ability of low-order faults in seismic data.

【技术实现步骤摘要】
基于优化Morlet小波的分频扩散滤波断层强化方法
本专利技术涉及油田开发
，特别是涉及到一种基于优化Morlet小波的分频扩散滤波断层强化方法。
技术介绍
在地震数据采集过程中，不可避免地会受到传感器灵敏度、野外噪声干扰和数模转换量化问题等因素的影响，导致地震数据信噪比低等问题，给后续构造解释带来很大困难。针对常规叠后地震数据噪声压制的问题，国内外学者开展了大量研究，提出了中值滤波、高斯滤波、时频域滤波等方法。上述方法在噪声压制方面效果较好，但由于无法识别地震反射的局部方向信息，在噪声压制的同时，会破坏边缘信息。为此，国内外学者将各向异性扩散滤波方法引入到地震数据滤波处理中。Perona等最先提出各项异性扩散方程(即P-M模型)，Fehmers等将该方法引入到地震资料的处理和解释中，Weickert和Lavialle等对P-M模型进行了改进，提出一致性增强扩散滤波(CED滤波)，通过使用扩散张量作为扩散系数，实现了真正意义上的非线性各项异性扩散滤波。此后，众多学者在CED滤波的基础上进行改进，并引入地震相干属性、非平稳相似系数等断层因子，提出诸如双向增强扩散滤波、构造导向扩散滤波等。上述扩散滤波方法的基本思想都是通过地震剖面局部图像的方向特征，判断是否存在断层等异常地质体边界，进而确定是否进行滤波，从而达到去噪的同时保护边界的目的。需要注意的是，滤波处理在提高信噪比的同时会降低地震数据的分辨率。为了解决上述问题，通过对分频地震数据进行各向异性扩散滤波处理，然后重构得到全频地震数据，可以一定程度提高滤波后地震数据的分辨率。目前，用于分频处理的主要是谱分解法、傅里叶变换、S变换和小波变换等。其中小波变换可以随着频率成分的改变而自动调整窗口形态，能够对信号进行局部化分析，被称为数学显微镜，广泛用于提高地震资料的分辨率。公开号为CN106054250A的专利申请变频分量扩散滤波融合的地震资料噪声消减方法中提到类似的实现过程：(1)输入三维叠后地震数据体；(2)将三维叠后地震数据体分解为多个变频分量数据体；(3)对变频分量数据体进行地震反射倾角导向扩散滤波处理，得到滤波后的变频分量地震数据体；(4)将滤波后的变频分量地震数据体进行加权重构与融合，获得压制了噪声，同时保留了有效的不连续异常新信息的三维地震数据体；(5)对步骤(4)得到的三维地震数据体求取地震体曲率属性；(6)输入目的层位数据，从上述三维数据体中抽取岩层切片，用于储层解释。基于分频扩散滤波的地震数据去噪方法的关键是获取高分辨率的分频地震数据，而利用小波变换进行分频处理的关键选取合适的小波母函数，其选取标准是小波母函数的波形和频谱与地震子波的波形与频谱相似，并满足容许条件。目前，Morlet小波具有良好的局部性能，且当Morlet小波主频大于1Hz时近似满足容许条件，被广泛用于地震数据的小波分频中，但传统的Morlet小波及前人改进后的Morlet子波均存在其波形和频谱无法同时匹配地震子波的波形和频谱的问题，从而导致分频结果分辨率较低。为此我们专利技术了一种新的基于优化Morlet小波的分频扩散滤波重构断层强化处理方法，解决了以上技术问题。
技术实现思路
本专利技术的目的是提供一种通过优化Morlet小波提高小波分频结果的分辨率及准确度，然后结合各项异性扩散滤波，实现地震数据去噪保边的同时提高其分辨率，从而提高地震数据低序级断层识别能力的基于优化Morlet小波的分频扩散滤波断层强化方法。本专利技术的目的可通过如下技术措施来实现：基于优化Morlet小波的分频扩散滤波断层强化方法，该基于优化Morlet小波的分频扩散滤波重构断层强化处理方法包括：步骤1，将叠后地震数据转换为MATLAB可识别的数据格式；步骤2，对地震数据进行频谱分析并提取子波；步骤3，参照地震子波的波形和频谱特征，改进Morlet子波母函数，构建Morlet子波函数组合；步骤4，利用优化Morlet子波进行小波分频处理；步骤5，计算叠后地震数据的结构张量，进而构建基于结构张量迹的新断层因子；步骤6，在结构张量迹的保护下，对分频地震数据进行各项异性扩散滤波；步骤7，对分频扩散滤波结果进行小波反变换，并重构得到提高了信噪比和分辨率，并保护了断层信息的地震数据。本专利技术的目的还可通过如下技术措施来实现：在步骤2中，利用快速傅里叶变换对地震数据进行频谱分析，并进行地震子波提取。在步骤3中，改进的Morlet子波母函数公式为：式中：fc为小波母函数的主频；t为时间。在改进Morlet小波母函数的基础上，构建不同尺度的Morlet小波函数组合：式中：fc为小波母函数的主频，a为尺度参数，该参数与分频数据主频fa的关系为：a＝fs*fc/fa，其中fs为地震数据采样频率。在步骤4中，利用优化Morlet小波开展叠后地震数据的小波分频处理，其小波变换公式为：式中：Wψ(a，b)为小波分解后的小波系数；ψ(t)为优化Morlet小波母函数，a为尺度因子，b为时移因子。在步骤5中，结构张量的计算公式如下所示：式中：fδ＝Gδ*f，Gδ为尺度为δ的高斯核函数，*为卷积运算符，f为地震数据；为梯度算子；T为矩阵转置符号；算子Sδ是在尺度δ下的结构张量，包含着与梯度一样的方向信息；在此基础上，对上述Sδ进一步利用尺度为ρ的高斯核函数Gρ进行滤波，得到另一尺度下的结构张量：Sρ＝Gρ*Sδ对于平稳变化的区域，其不同尺度滤波后的结构张量反映的方向信息应该一致，而对于非平稳变化的区域，受高斯核函数尺度参数的影响，其滤波后反映的方向信息可能发生变化；因此，基于结构张量特征值的变化构建断层因子：式中：Sρ和Sδ是经过两个不同尺度滤波的结构张量；Tr表示矩阵的迹，其值为矩阵特征值的和；由ε的计算公式可知，ε取值范围为0<ε<1，ε取值较大时，反映图像平稳变化区域，而ε取值较小，则反映图像边界区域。在步骤6中，在地层连续区域，进行扩散滤波，在断层附近，减弱或停止扩散滤波，从而保护断层信息；针对不同分频结果，利用步骤3构建的结构张量构建扩散张量；对Sρ进行特征值分解得出特征值μ1≥μ2和对应的特征向量(v1，v2)；扩散张量的特征向量取结构张量的特征向量，使得扩散张量能够反映图像梯度变化，扩散张量的特征值则由研究人员赋值及其结构张量的特征值计算得到：式中：α用来控制沿图像梯度变化较大方向的扩散强度，为了保护图像的结构信息，α通常设为接近0的正数，C为阈值，通常取值为1，m取值为1；在引入结构张量迹后，扩散张量D构建如下：式中：a、b、c是扩散张量为矩阵表示形式时的参数。最终，基于结构张量迹保护的各向异性扩散滤波计算公式如下：式中：div为散度计算符号，▽为梯度算子，f为地震数据，分别为地震数据在x方向、对角方向和y方向上的二阶偏导。在步骤7中，小波反变换公式为：式中：ψ(t)为优化Morlet小波母函数，t为时间，ω为圆频率，Wψ(a，b)为小波分解后的小波系数；a为尺度因子，b为时移因子。该基于优化Morlet小波的分频扩散滤波断层强化方法还包括，在步骤7之后，在完成断层强化处理后，将处理后数据输出为标准SGY格式的数据。本专利技术中的基于优化Morlet小波的分频扩散滤波重构断层强化处理方法，在前人基础上进一步改本文档来自技高网...

【技术保护点】
1.基于优化Morlet小波的分频扩散滤波断层强化方法，其特征在于，该基于优化Morlet小波的分频扩散滤波重构断层强化处理方法包括：步骤1，将叠后地震数据转换为MATLAB可识别的数据格式；步骤2，对地震数据进行频谱分析并提取子波；步骤3，参照地震子波的波形和频谱特征，改进Morlet子波母函数，构建Morlet子波函数组合；步骤4，利用优化Morlet子波进行小波分频处理；步骤5，计算叠后地震数据的结构张量，进而构建基于结构张量迹的新断层因子；步骤6，在结构张量迹的保护下，对分频地震数据进行各项异性扩散滤波；步骤7，对分频扩散滤波结果进行小波反变换，并重构得到提高了信噪比和分辨率，并保护了断层信息的地震数据。

【技术特征摘要】
1.基于优化Morlet小波的分频扩散滤波断层强化方法，其特征在于，该基于优化Morlet小波的分频扩散滤波重构断层强化处理方法包括：步骤1，将叠后地震数据转换为MATLAB可识别的数据格式；步骤2，对地震数据进行频谱分析并提取子波；步骤3，参照地震子波的波形和频谱特征，改进Morlet子波母函数，构建Morlet子波函数组合；步骤4，利用优化Morlet子波进行小波分频处理；步骤5，计算叠后地震数据的结构张量，进而构建基于结构张量迹的新断层因子；步骤6，在结构张量迹的保护下，对分频地震数据进行各项异性扩散滤波；步骤7，对分频扩散滤波结果进行小波反变换，并重构得到提高了信噪比和分辨率，并保护了断层信息的地震数据。2.根据权利要求1所述的基于优化Morlet小波的分频扩散滤波断层强化方法，其特征在于，在步骤2中，利用快速傅里叶变换对地震数据进行频谱分析，并进行地震子波提取。3.根据权利要求1所述的基于优化Morlet小波的分频扩散滤波断层强化方法，其特征在于，在步骤3中，改进的Morlet子波母函数公式为：式中：fc为小波母函数的主频；t为时间；在改进Morlet小波母函数的基础上，构建不同尺度的Morlet小波函数组合：式中：fc为小波母函数的主频，a为尺度参数，该参数与分频数据主频fa的关系为：a＝fs*fc/fa，其中fs为地震数据采样频率。4.根据权利要求1所述的基于优化Morlet小波的分频扩散滤波断层强化方法，其特征在于，在步骤4中，利用优化Morlet小波开展叠后地震数据的小波分频处理，其小波变换公式为：式中：Wψ(a，b)为小波分解后的小波系数；ψ(t)为优化Morlet小波母函数，a为尺度因子，b为时移因子。5.根据权利要求1所述的基于优化Morlet小波的分频扩散滤波断层强化方法，其特征在于，在步骤5中，结构张量的计算公式如下所示：式中：fδ＝Gδ*f，Gδ为尺度为δ的高斯核函数，*为卷积运算符，f为地震数据；为梯度算子；T为矩阵转置符号；算子Sδ是在尺度δ下的结构张量，包含着与梯度一样的方向信息；在此基础上，对上述Sδ进一步利用尺度为ρ的高斯...

【专利技术属性】
技术研发人员：周海廷，刘显太，杜玉山，王建，武刚，赵开连，黄超，房环环，谭琴辉，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司勘探开发研究院，
类型：发明
国别省市：山东,37