[发明专利]基于AVO异常类型约束的剩余相位校正方法

权利要求书

1.一种基于AVO异常类型约束的剩余相位校正方法，步骤包括：

(1)激发地震波并记录工区地震波，按常规地震资料处理流程对记录的地震资料进行前期基础处理包括数据解编、去噪、振幅补偿、速度分析、压制多次波、叠前时间偏移、常规动校正步骤，生成可使用的CMP道集或CRP道集；

(2)对输入的CMP/CRP道集进行傅里叶变换，将其由时间域转换到频率域，在频率域中进行振幅及相位计算，并对振幅谱进行去除噪音处理；

(3)对处理后的CMP/CRP道集进行相位动校正PMO，如下：

对于动校正后某地震道时间序列s(t)可写为：

其中i为地震道时间序列的序号，n为地震道的采样点数目，ti及ai分别为到达时和对应的振幅值；δ(t-ti)为δ函数，δ函数为一个广义函数，该函数在除了t-ti＝0以外的点取值都等于零，而其在整个定义域上的积分等于1；

s(t)的振幅谱为：

其中ω为角频率，式(2)对应的相位谱为：

其中ti，tj为所有到达时，(ti-tj)为所有到达时的时差；

(4)对经过相位动校正处理后的CMP/CRP道集进行相位统计及相位滤波；

(5)对CMP/CRP道集进行频率域PMO修正，对于具有第一、二类极性反转的AVO资料，采用相位谱条件约束，在做 PMO之前，用井的资料做AVO正演模型，分析含气砂岩的类型，方法如下：

A.远、中、近叠加道相位组合法：

其中Φa(ω)、Φb(ω)、Φc(ω)分别表示近、中、远道叠加相位谱，Wa、Wb、Wc是相对应的加权系数；

B.约束条件相位谱法：

定义为第j道条件约束相位谱，其中Φc为约束相位；

取Φc＝0、π/4、π、2π；其实，只有当实际地震道相位谱Φj(ω)与模型道相位谱Φm(ω)之差异大于π值时，使地震道相位产生反转；

设模型道的相位谱为Φm(ω)，则PMO算法为：

sj(t)＝FFT-1[Aj(ω)exp(iΦm(ω))] (5)

其中sj(t)为PMO算法处理后的CMP/CRP道集时间序列，j为CMP/CRP道集中的道序号，Aj(ω)为第j道的振幅谱，FFT-1表示为逆FFT过程；

定义：

其中k为CMP/CRP道集中的道序号，则PMO的模型道相位谱可取用：

Φm(ω)＝limj→JΦj(ω),J∈(1,n) (7)

(6)对完成频率域PMO修正的CMP/CRP道集进行反傅里叶变换，将其由频率域转换返还至时间域；

(7)将时间域的CMP/CRP道集直接输出得到时域PMO道集；

(8)将输出得到的时域PMO道集应用于AVO分析解释，进行油气检测。