[发明专利]一种高维地震数据规则化方法

权利要求书

1.一种高维地震数据规则化方法，其特征在于：所述方法包括：

(1)输入原始三维炮道集；

(2)利用高维傅里叶变换将时间空间域的所述原始三维炮道集变换到频率空间域中，然后进行分选得到五维炮道集；

(3)对一个频率片进行最优化求解，得到无假频频谱；

(4)判断是否达到迭代门槛值，如果是，则转入步骤(5)，如果否，则返回步骤(3)；

(5)利用高维傅里叶反变换将所述无假频频谱变换到时间空间域，即得到规则炮道集数据；

(6)输出所述规则炮道集数据,

所述步骤(3)是采用共轭梯度算法求解下面的公式实现的：

$\begin{array}{c}\tilde{m}=m_{prior}+{(A^H{C}_D^{-1}A+C_M^{-1})}^{-1}{A}^H{C}_D^{-1}(d_{obs}-{\mathrm{Am}}_{prior})\\ =m_{prior}+C_M{A}^H{({\mathrm{AC}}_M{A}^H+C_D)}^{-1}(d_{obs}-{\mathrm{Am}}_{prior})\end{array}$

其中，是无假频的规则数据频谱，即无假频频谱，m_prior是先验模型，C_M为数据的Radon谱，C_D为一般约束条件，d_obs为观测数据，A表示采样算子矩阵，A^H为A的共轭转置矩阵，m_prior＝0。