地震数据梯度信息不连续性边界检测方法技术

本发明专利技术涉及一种地震数据处理方法，具体为地震数据梯度信息不连续性边界检测方法，输入三维地震数据，计算梯度向量体，建立二维滑动分析时窗，二维滑动分析时窗分为5个子窗口，计算滑动时窗内每个样点的梯度幅值；再计算子窗口的峰态值，比较每个子窗口的峰态值；将峰态值最大的那个子窗口内所有的梯度向量构成梯度张量矩阵，求取梯度结构张量矩阵的特征值，并将三个特征值进行排序，利用排序后的特征值大小构建梯度张量属性，将求取的梯度张量属性赋给位于二维滑动分析时窗内中心的分析点；返回计算下一个分析点的梯度张量属性值。该方法既考虑了时窗内位于不同区域位置数据，又考虑了时窗内数据本身的不稳定性，充分利用了各种信息。

【技术实现步骤摘要】

本专利技术涉及一种地震数据处理方法，具体为地震数据梯度信息不连续性边界检测方法。
技术介绍
地震数据中含有大量丰富的地层和岩性信息，这些信息很多都是储层中断层、河道、裂缝发育区以及其他一些不连续性结构特征的响应，但是利用常规的地震数据解释方法并不能够得到清晰直观的不连续性信息。2002年，Cohen和Coifman提出了基于局部结构熵的相干性算法。该算法把由地震数据估计的局部结构熵作为相干测定，先构造一个分析数据体，并将它分为4个子数据体，并分别生成4个列向量，利用4个列向量的互相关生成一个4×4的相关矩阵，将该矩阵的归一化道作为局部熵估计。该算法从计算上讲是非常有效的，因为它避免了协方差矩阵中大量本征值的计算，不足之处是没有考虑地下构造的倾角对局部结构熵估计的影响。2005年，Lu Wenkai和Li Yandong等提出了基于高阶累积量的相干估计方法(HOSC)。高阶累积量(HOS)成功用于信号处理，而相干分析的基本问题就是延迟估计，这正是HOS方法的特点。HOSC方法同时利用3个地震道来计算拥有零延迟相关的归一化4阶矩的二维切片，将二维切片上最大相关点作为相干估计。另外，Lu Wenkai提出了超级地震道(ST-supertrace)概念。先定义一个分析时窗，它包含分析点周围的许多地震道，然后重新排列多道为超道，最后在整个三维数据体中移动分析时窗，这样就把原始地震数据体转换为超级地震道数据体。ST数据体保持了原始数据体中所拥有的地层倾角信息，故在ST数据体中利用倾角扫描很容易获取最佳倾角估计。把ST与HOSC算法组合起来可以得到有效的相干估计算法ST-HOSC。该算法在检测不连续性信息方面比C1相干体和HOS相干体更好。以上两种现有的检测不连续性边缘的算法，都是在一个整体时窗内所有的地震数据进行分析处理，并没有考虑整体时窗内数据本身的分布情况和结构特征，而且没有考虑时窗内不同位置区域的数据对中心分析点的贡献。
技术实现思路
针对上述技术问题，本专利技术提供一种新的不连续性边界的检测方法。该不连续性边界检测算法从地震数据梯度信息的角度出发，在构建滑动时窗时既考虑时窗内位于不同区域位置数据对中心分析点的贡献程度，又考虑时窗内数据本身的不稳定性。具体的技术方案是：地震数据梯度信息不连续性边界检测方法，包括以下步骤：(1)输入三维地震数据，并利用式(1)～(3)计算地震数据在x，y，z三个方向上的梯度向量体； g x = f ( x i + 1 , y , z ) - f ( x i , y , z ) Δx - - - ( 1 ) ; ]]> g y = f ( x , y j + 1 , z ) - f ( x , y j , z ) Δy - - - ( 2 ) ; ]]> g z = f ( x , y , z k + 1 ) - f ( x , y , z k ) Δz - - - ( 3 ) ; ]]>上式中，f(x，y，z)为三维地震数据；Δx、Δy和Δz分别为地震数据在x，y，z三个方向上的面元大小；gx、gy和gz分别为地震数据在x、y、z三个方向上的梯度向量体；(2)建立二维滑动分析时窗，二维滑动分析时窗分为五个3×3的子窗口，分别位于二维滑动分析时窗的中心位置、左上位置、右上位置、左下位置和右下位置；(3)利用式(4)计算滑动时窗内每个样点的梯度幅值；再利用式(5)计算二维滑动分析时窗内5个子窗口的峰态值： g = g x 2 + g y 2 + g z 2 本文档来自技高网...

【技术保护点】
地震数据梯度信息不连续性边界检测方法，其特征在于：包括以下步骤：(1)输入三维地震数据，并利用式(1)～(3)计算地震数据在x，y，z三个方向上的梯度向量体；gx=f(xi+1,y,z)-f(xi,y,z)Δx---(1);]]>gy=f(x,yj+1,z)-f(x,yj,z)Δy---(2);]]>gz=f(x,y,zk+1)-f(x,y,zk)Δz---(3);]]>上式中，f(x，y，z)为三维地震数据；Δx、Δy和Δz分别为地震数据在x，y，z三个方向上的面元大小；gx、gy和gz分别为地震数据在x、y、z三个方向上的梯度向量体；(2)建立二维滑动分析时窗，二维滑动分析时窗分为五个3×3的子窗口，分别位于二维滑动分析时窗的中心位置、左上位置、右上位置、左下位置和右下位置；(3)利用式(4)计算滑动时窗内每个样点的梯度幅值；再利用式(5)计算二维滑动分析时窗内5个子窗口的峰态值：g=gx2+gy2+gz2---(4);]]>K=1JΣi=1J(gi-g‾)4---(5);]]>上式中，g为计算的梯度幅值；gi为子窗口内第i个样点的梯度幅值大小；为子窗口内所有梯度幅值的平均值；K和J分别为子窗口的峰态值和样点数；(4)比较二维滑动分析时窗内每个子窗口的峰态值；(5)将峰态值最大的那个子窗口内所有的梯度向量构成梯度张量矩阵，即式(6)；T‾=g·gT‾=1JΣj=1Jgx2Σj=1JgxgyΣj=1JgxgzΣj=1JgygxΣj=1Jgy2Σj=1JgygzΣj=1JgzgxΣj=1JgzgyΣj=1Igz2---(6)]]>上式中gx、gy和gz分别为子窗口内的地震数据在x、y、z三个方向上的梯度向量；J为子窗口内的样点个数；(6)利用QR分解法求取梯度结构张量矩阵的特征值，并将三个特征值进行排序，λ1＞λ2＞λ3；(7)利用排序后的特征值大小构建梯度张量属性，即式(7)；AG=λ11NΣi=1N(λi-λ‾)3---(7)]]>上式中λ1为最大特征值；N为特征值的个数；为所有特征值的平均值；(8)将求取的梯度张量属性赋给位于二维滑动分析时窗内中心的分析点；返回步骤(3)，滑动二维滑动分析时窗并计算下一个分析点的梯度张量属性值；(9)将整个三维地震数据依据以上过程计算完毕，得到最后的三维梯度张量属性体。...

【技术特征摘要】
1.地震数据梯度信息不连续性边界检测方法，其特征在于：包括以下步骤：
(1)输入三维地震数据，并利用式(1)～(3)计算地震数据在x，y，z
三个方向上的梯度向量体；
g x = f ( x i + 1 , y , z ) - f ( x i , y , z ) Δx - - - ( 1 ) ; ]]> g y = f ( x , y j + 1 , z ) - f ( x , y j , z ) Δy - - - ( 2 ) ; ]]> g z = f ( x , y , z k + 1 ) - f ( x , y , z k ) Δz - - - ( 3 ) ; ]]>上式中，f(x，y，z)为三维地震数据；Δx、Δy和Δz分别为地震数据
在x，y，z三个方向上的面元大小；gx、gy和gz分别为地震数据在x、y、z
三个方向上的梯度向量体；
(2)建立二维滑动分析时窗，二维滑动分析时窗分为五个3×3的子窗口，
分别位于二维滑动分析时窗的中心位置、左上位置、右上位置、左下位置和右下
位置；
(3)利用式(4)计算滑动时窗内每个样点的梯度幅值；再利用式(5)计
算二维滑动分析时窗内5个子窗口的峰态值：
g = g x 2 + g y 2 + g z 2 - - - ( 4 ) ; ]]> K = 1 J Σ i = 1 J ( g i - g ‾ ) 4 - - - ( 5 ) ; ]]>上式中，g为计算的梯度幅值；gi为子窗口内第i个样点的梯度幅值大小；
为子窗口内所有梯度幅值的平均值；K和J分别为子窗口的峰态值和样点数；
(4)比较二维滑动分析时窗内每个子窗口的峰态值；
(5)将峰态值最大的那个子窗口内所有的梯度向量构成梯度张量矩阵，即
式(6)；
T ‾ = g · g T ‾ = 1 J Σ j = 1 J g x 2 Σ j = ...

【专利技术属性】
技术研发人员：丁峰，彭达，尹成，范廷恩，罗浩然，刘伟，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：四川;51