[发明专利]一种反褶积方法及装置

权利要求书

1.一种反褶积方法，其特征在于，包括以下处理步骤：

S1：将所需处理的地震数据处理成叠后地震数据记录；

S2：构建包括反射系数序列在内的目标函数；

S3：在所述叠后地震数据记录中选取一道所需处理的地震资料，从所述构建的目标函数中计算出所述地震资料的反射系数序列的振幅和位置；

S4：读取下一道所需处理的地震资料，计算出所述下一道地震资料反射系数序列的振幅和位置；

S5：计算出所述叠后地震数据记录中所有地震道反射系数序列的振幅和位置，输出反射系数序列，完成所述叠后地震数据记录的反褶积处理。

2.如权利要求1所述的一种反褶积方法，其特征在于，所述构建目标函数包括：

S201：建立地震数据褶积模型；

S202：将所述地震数据褶积模型转换成矩阵形式；

S203：由所述矩阵形式构建目标函数；所述构建的目标函数为：

J_α＝||Aα-S||²

上式中，J_α为关于反射系数序列α的范数，A为一个N×M的矩阵，其表达式为其中，ω_t为地震子波，t＝1，2，…N，N为所述叠后地震数据记录的时间窗长度，M为预先设置的所述时间窗长度N内反射系数序列的个数，τ_M第M个反射系数的时延，s为t时刻所述叠后地震数据记录的值。

3.如权利要求2所述的一种反褶积方法，其特征在于，所述构建的目标函数为：

J＝||Aα-S||²+μ||Cα-ε||²＝J_α+μJ_c

上式中，μ为阻抗权系数，C为左下三角全为1的矩阵，ε为波阻抗。

4.如权利要求3所述的一种反褶积方法，其特征在于，所述构建的目标函数为：

J＝J_α+μJ_c+βJ_α

上式中，β为正则化的阻尼系数，其中，β＝β₀f₀，β₀为预设值，f₀＝max{(A^TA+μC^TC_ii)}，i为矩阵(A^TA+μC^TC)的下标。

5.如权利要求1-4任意一项所述的一种反褶积方法，其特征在于，从所述构建的目标函数中计算出所述地震资料的反射系数序列的振幅包括步骤：

S310：从所述叠后地震数据记录中估计出地震子波；

S320：将所述估计出的地震子波代入所述目标函数中，计算出所述反射系数序列的振幅；

其中，所述估计地震子波的方法包括步骤：

S311：建立所述叠后地震数据记录的褶积模型；

S312：将所述建立的褶积模型变换到复赛谱域；

S313：设计低通滤波器，在复赛谱域将所述建立的褶积模型的地震子波和反射系数序列进行分离，得到复赛谱域的地震子波；

S314：将所述复赛谱域的地震子波转换到时间域，得到时间域的地震子波。

6.如权利要求5所述的一种反褶积方法，其特征在于，所述地震子波的估计方法还包括步骤：

315：对所述估计出的地震子波相位进行校正。