[发明专利]一种异常初至时间修正方法

权利要求书

1.一种异常初至时间修正方法，其特征在于，包括以下步骤：

A、根据采集的地震数据计算每道每一个采样点的能量比值，得到所述道的能量比值曲线；

B、从所述道的能量比值曲线中得到最大值对应的采样点，并将采样点作为相应地震道初至时间的样点；

C、利用滑动时窗对初至时间线性拟合，根据初至时间与拟合线的距离确定可靠初至时间；

D、计算可靠初至时间的正负波形面积、两个面积的比值、比值的均值以及初至波峰点能量的均值；

E、计算异常初至道在预设的时窗内每个子波的正负波形面积的比值和每个子波的波峰点能量；

F、根据上述得到的数据计算异常初至道在预设的时窗内每个子波的综合评价值，确定异常初至时间的正确初至时间；

异常初至道在预设的时窗内每个子波的综合评价公式如下：

式中：S为子波的综合评价值，最小值对应的子波为第k₁道的初至时间，和分别为第k₁道预设时窗内某个子波的波峰点能量和正负波形面积比值，和分别为与第k₁-1道相邻的l道可靠初至波峰点能量的均值和正负波形面积比值的均值，ε₁和ε₂为加权系数，ε₁+ε₂＝1，l为预设的可靠初至时间的道数；

G、滑动着对所有异常初至时间进行计算修正。

2.根据权利要求1所述的一种异常初至时间修正方法，其特征在于，所述步骤A中采集的地震数据利用预设的能量比值公式计算每道每一个采样点的能量比值，其预设的能量比值公式如下：

式中：R_i(r)代表第i道第r个采样点的能量比值，w为计算时窗的点数，i为地震道道号，1≤i≤m，r为对应的第i道的采样点点号，1≤r≤n，n为一道的采样点数，α为稳定因子，β为大于0的整数。

3.根据权利要求1所述的一种异常初至时间修正方法，其特征在于，所述步骤C的具体过程为：

a、以第一道初至时间为起点，滑动时窗长度为L，滑动步长为1，逐道滑动着对拾取的初至时间线性拟合，计算滑动时窗内每道初至时间与拟合直线的时间差，若滑动时窗内某道大于给定的阈值δ，则为异常初至，否则为可靠初至；

b、以第一道初至时间为起点，确定第k₁道为第一个异常初至时间，逐道滑动着向后寻找，直到第k₂道为异常初至，且第k₂+1道为可靠初至时，将第k₁至k₂道之间的初至时间标记为异常初至段，对所有拾取的异常初至时间进行标记。

4.根据权利要求1所述的一种异常初至时间修正方法，其特征在于，所述步骤D中可靠初至时间的正负波形面积比值的计算公式如下：

式中：为第k₁-d道可靠初至时间的正负波形面积比值，和分别为第k₁-d道初至时间的正波形面积和负波形面积，且d为大于等于1小于等于l的正整数，l为预设的可靠初至时间的道数。

5.根据权利要求1所述的一种异常初至时间修正方法，其特征在于，所述步骤D中比值的均值的计算公式如下：

式中：为与第k₁-1道相邻的可靠初至时间的正负波形面积比值的均值，为第k₁-d道可靠初至时间的正负波形面积比值，d为大于等于1小于等于l的正整数，l为预设的可靠初至时间的道数。

6.根据权利要求1所述的一种异常初至时间修正方法，其特征在于，所述步骤D中初至波峰点能量的均值的计算公式如下：

式中：为与第k₁-1道相邻的l道可靠初至波峰点能量的均值，为第k₁-d道可靠初至波峰点的能量，且d为大于等于1小于等于l的正整数，l为预设的可靠初至时间的道数。