[发明专利]一种基于最小二乘参数反演的傅里叶域地震数据重构方法

权利要求书

1.一种基于最小二乘参数反演的傅里叶域地震数据重构方法，其特征在于：所述方法基于最小二乘参数反演，对不含假频的非均匀采样的地震数据进行插值，从而重构地震数据；

其中，所述方法包括：

(1)对原始地震数据即共炮道集进行预处理，得到NMO校正后的数据D；

(2)对于NMO校正后的数据D，计算：b＝A^H·W·D，

其中A^H是A的共轭转置矩阵，W是加权值，它是一个与局部采样密度成反比的值，m_kx，m_ky为空间波数采样序号，Δk_x，Δk_y为空间波数间隔，ΔS_F＝(X_i+1-X_i)(Y_i+1-Y_i)，(X_i，Y_i)、(X_i+1，Y_i+1)为空间上两点，x_n表示地震数据所在的空间横坐标，y_n表示地震数据的所在的空间纵坐标；

(3)计算：H＝A^H·W·A+λ²·I，然后求其逆H^-1，其中σ是先验模型变差，c是常数；I是单位对角矩阵；

(4)由第(2)、第(3)步得到的b和H^-1计算

(5)计算插值重构结果：

2.根据权利要求1所述的基于最小二乘参数反演的傅里叶域地震数据重构方法，其特征在于：所述步骤(1)中的预处理包括：置道头、去噪、振幅补偿、静校正、反褶积、分选得到CMP道集并对CMP道集进行NMO校正。