[发明专利]一种数据驱动的地震信号零相位反褶积方法

摘要

本发明公开了一种数据驱动的地震信号零相位反褶积方法，包括：1)从地震记录中估计地震子波振幅谱；2)基于估计的振幅谱，对地震记录在频率域进行维纳滤波；3)对维纳滤波器做逆傅里叶变换，得到反褶积结果。本发明利用子波振幅谱估计方法，不依赖于模型参数，能够较好地对子波振幅谱进行有效地估计，具有好的鲁棒性。用该方法对地震记录进行零相位反褶积，可以收到了很好的效果，同时该方法也可以有效地处理非平稳地震数据。

主权项

1.一种数据驱动的地震信号零相位反褶积方法，其特征在于，包括以下步骤：1)从地震记录中估计地震子波振幅谱，采用COM方法，具体方法如下：其估计子波使用的是压缩映射算子，表达式为：P[f,x]＝cp(Fp(f,x))α(1‑Fp(f,x))β   (1)其中：f为待估计函数，x为函数的自变量，cp为归一化因子，Fp为待估计函数的泛函，α为参数，β为参数；其中，fp(t)为p次方，a为积分下界，b为积分上界；在实际应用上式估计地震子波振幅谱时，α,β需要估计，而p是事先设定的，p的选取不影响子波估计计算精度，影响子波估计计算效率；2)基于估计的振幅谱，对地震记录在频率域进行维纳滤波，其中，频率域维纳滤波的具体方法如下：频率域维纳滤波的标准形式为：其中，S(ω),W(ω),ε2分别为地震记录的频谱、地震子波的频谱和正则化参数；ω为角频率，ε为正则化因子；W*(ω)表示对子波谱取共轭；采用如下的改进形式：其中，θ的值由地震信号的噪声水平决定，当地震信号的噪声水平高于一般值时，θ的取值大于1；当地震信号的噪声水平低于一般值时，θ的取值小于1；3)对维纳滤波结果做逆傅里叶变换，得到反褶积结果。