【技术实现步骤摘要】
本专利技术涉及一种地震数据处理方法,具体为地震数据梯度信息不连续性边界检测方法。
技术介绍
地震数据中含有大量丰富的地层和岩性信息,这些信息很多都是储层中断层、河道、裂缝发育区以及其他一些不连续性结构特征的响应,但是利用常规的地震数据解释方法并不能够得到清晰直观的不连续性信息。2002年,Cohen和Coifman提出了基于局部结构熵的相干性算法。该算法把由地震数据估计的局部结构熵作为相干测定,先构造一个分析数据体,并将它分为4个子数据体,并分别生成4个列向量,利用4个列向量的互相关生成一个4×4的相关矩阵,将该矩阵的归一化道作为局部熵估计。该算法从计算上讲是非常有效的,因为它避免了协方差矩阵中大量本征值的计算,不足之处是没有考虑地下构造的倾角对局部结构熵估计的影响。2005年,Lu Wenkai和Li Yandong等提出了基于高阶累积量的相干估计方法(HOSC)。高阶累积量(HOS)成功用于信号处理,而相干分析的基本问题就是延迟估计,这正是HOS方法的特点。HOSC方法同时利用3个地震道来计算拥有零延迟相关的归一化4阶矩的二维切片,将二维切片上最大相关点作为相干估计。另外,Lu Wenkai提出了超级地震道(ST-supertrace)概念。先定义一个分析时窗,它包含分析点周围的许多地震道,然后重新排列多道为超道,最后在整个三维数据体中移动分析时窗,这样就把原始地震数据体转换为超级地震道数据体 ...
【技术保护点】
地震数据梯度信息不连续性边界检测方法,其特征在于:包括以下步骤:(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为子窗口 ...
【技术特征摘要】
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
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。