[发明专利]地震波数据的处理方法及装置

权利要求书

1.一种地震波数据的处理方法，包括：

将共炮域地震波数据转换为共偏移距地震波数据，所述共偏移距地震波数据包括压缩波数据和剪切波数据；

进行正演计算：输入初始压缩波速度模型和初始剪切波速度模型，确定计算数据函数；以及输入当前迭代后的压缩波速度模型和当前迭代后的剪切波速度模型，更新计算数据函数；

进行反演计算：构建目标函数，输入所述计算数据函数、压缩波数据和剪切波数据进行迭代，确定当前迭代后的压缩波速度模型和当前迭代后的剪切波速度模型；以及输入当前迭代后的压缩波速度模型和当前迭代后的剪切波速度模型，更新目标函数；以及

重复迭代，当迭代次数达到预设迭代次数时，迭代暂停，判断目标函数的结果是否达到预定结果，若达到该预定结果，则停止迭代；若未达到该预定结果，则继续迭代，直至达到该预定结果，输出最终压缩波速度模型和最终剪切波速度模型。

2.根据权利要求1所述的方法，其中，所述目标函数包括交叉梯度函数，其公式为||α_hh(m_p，m_s)||²，其中，||||²指求模平方，h(m_p，m_s)为压缩波速度模型m_p和剪切波速度模型m_s的交叉梯度，其公式为为梯度算子，α_h是交叉梯度的权重因子，所述压缩波速度模型包括初始压缩波速度模型、当前迭代后的压缩波速度模型和最终压缩波速度模型，剪切波速度模型包括初始剪切波速度模型、当前迭代后的剪切波速度模型和最终剪切波速度模型，且所述压缩波速度模型与所述剪切波速度模型一一对应。

3.根据权利要求2所述的方法，其中，所述目标函数还包括拟合残差函数和正则化函数，其中，所述拟合残差函数的公式为d_p为压缩波数据，d_s为剪切波数据，u为计算数据函数；正则化函数的公式为其中，L为正则化矩阵，α_p和α_s分别为L(m_p)和L(m_s)的权重因子。

4.根据权利要求3所述的方法，其中，所述目标函数的公式为：

5.根据权利要求4所述的方法，其中，当前迭代与前一次迭代的速度模型m的增加量其中，m包括所述压缩波速度模型m_p与所述剪切波速度模型m_s，Δm为下一次迭代的速度模型增加量，为φ(m)的偏导。

6.根据权利要求5所述的方法，其中，在进行正演计算之前，还包括步骤：通过时窗选取共偏移距地震波数据中的压缩波数据和剪切波数据；以及

采用数值计算所述预设地震波在传输的介质中作e的微小扰动，计算近似公式其中，m为介质的速度模型量，e大于0且远小于1。

7.根据权利要求5所述的方法，其中，所述计算数据函数的公式为：

其中，t为时间，k为波数，i为虚部，n为序数，为基坐标向量，

Y_n(r，θ)＝J_n(kr)e^inθ，J为以m_p和m_s为自变量的贝塞尔函数。