[发明专利]一种方位高阶相干振幅变化率裂缝检测方法

权利要求书

1.一种方位高阶相干振幅变化率裂缝检测方法，其特征在于：所述方法包括以下步骤：

(1)叠后地震资料处理：采集地震数据，并对其进行处理得到叠前地震资料，对所述叠前地震资料进行叠加处理得到叠后地震数据体；

(2)对步骤(1)得到的叠后地震数据体中的一个时间点上的宏面元的中心点进行不同方位宏面元地震道的相关校正；

(3)计算不同方位地震道高阶振幅变化率：将该中心点的地震道与步骤(2)得到的相关校正后的相邻点的地震道做高阶振幅变化率计算，得到该中心点在四个方向上的地震道的高阶振幅变化率；

(4)最大方位高阶振幅变化率提取：取步骤(3)得到的四个方向上的地震道的高阶振幅变化率中的最大值，即得到该中心点的地震道在该时间上的最大方位高阶振幅变化率；

(5)多次重复步骤(2)到步骤(4)获得该中心点在各个时间点上的最大方位高阶振幅变化率，即该中心点在整个地震道上的最大方位高阶振幅变化率；

(6)将整个地震工区其它各个点分别作为中心点，然后重复步骤(2)至(5)，就得到整个叠后地震数据体的最大方位高阶振幅变化率。

2.根据权利要求1所述的方位高阶相干振幅变化率裂缝检测方法，其特征在于：所述步骤(2)中采用5*5或3*3的宏面元。

3.根据权利要求2所述的方位高阶相干振幅变化率裂缝检测方法，其特征在于：具体如下：所述步骤(2)具体如下：

在叠后地震数据体中取一个时间点上的宏面元的中心点，采用下面的公式对该中心点的地震道与相邻点的地震道在一个时窗T内做互相关计算：

rXY=Σi=1N(Xi-X‾)(Yi-Y‾)Σi=1N(Xi-X‾)2Σi=1N(Yi-Y‾)2]]>

式中，X为中心点的变量，Y为其他每个点做互相关计算的变量，Xi、Yi为变量的标志值，i＝1至N，N表示一个时窗T内地震数据点的总个数，分别为变量的平均值。