基于最小二乘梯度更新速度模型的地震波全波形反演方法
【技术领域】
[0001] 本发明涉及一种地震波全波形反演方法，特别是关于一种基于最小二乘梯度更新 速度模型的地震波全波形反演方法。
【背景技术】
[0002] 地震波场反演方法很多，如相位反演成像、振幅反演成像以及全波形反演成像。相 位反演和振幅反演利用的是波场运动学信息，反演出的速度分布是一种平滑估计值，它是 真实速度模型的低频分量，反映的是大尺度的速度结构。基于波动方程理论的全波形反演 综合利用了地震记录中振幅、走时和相位等完整的波场信息，通过拟合，全波形反演通常利 用给定时窗内的理论波场记录与实测波场记录差值的二次泛函以及额外的地质约束条件 建立优化目标函数，利用迭代线性化反演方法或完全非线性反演法求解符合条件的模型， 定量提取地下介质的弹性参数，进而为深部大尺度构造演化分析，为勘探地震成像及速度 建模等方面提供可靠依据。
[0003]目前，全波形反演的方式大致可以分为时间域全波形反演与频率域全波形反演。 频率域全波形反演理论与时间域全波形反演理论有着良好的对应关系。两种反演方法的观 测数据表现形式不同造成正演波场计算方式的差别，而反演思想完全一致，与时间域相比， 频率域波场是相互解耦的。不同频率数据对异常体反映能力不一，根据需要既可以利用部 分频段数据反演，也可以使用全频段数据同时进行反演。
[0004] 全波形反演可利用逆时偏移产生的成像结果近似来更新速度模型，基于逆时偏移 的全波形反演获得低频成像不再需要炮集的低频信息且具有很高的保真度。但是，反演使 用大量低频的信息来更新背景速度，而偏移使用高频信息来勾画地下构造的比较精确的边 界，由于反演和偏移的目标在不同的频带，因此，在实际处理中，由于噪声和模型假象，高保 真度很难达到，且常规全波形反演流程中，每次的全波形反演的非线性迭代过程中，只利用 一次的逆时偏移的成像结果作为速度模型更新梯度的方向，这样容易导致更新梯度计算的 更新量不足，且收敛速度较慢。
【发明内容】
[0005] 针对上述问题，本发明的目的是提供一种可有效提高速度模型更新梯度的准确性 且可快速完成速度模型更新的基于最小二乘梯度更新速度模型的地震波全波形反演方法。
[0006] 为实现上述目的，本发明采取以下技术方案：一种基于最小二乘梯度更新速度模 型的地震波全波形反演方法，包括以下步骤：
[0007] 1)在时间域内，地震波场的声波方程为：
[0008]
[0009]其中，u是波场向量，s是震源向量，x和t分别是空间和时间坐标参数，M为质量 矩阵，A为刚度矩阵；
[0010] 2)构建初始速度模型，设置速度模型更新迭代次数N与速度模型允许的最小误差 值e;
[0011] 3)构造观测波场数据和计算波场数据的波场误差矢量Ad;
[0012] 4)利用最小二乘流程和波场误差矢量Ad，构造目标函数C(mk)，目标函数C(mk)的 表达式为：
[0013] C(mk) = 1/2AdTAd；
[0014] 式中，T表示转置；
[0015] mk为第k次更新的速度模型，更新公式为：
[0016] mk=mk !+Am；
[0017] 式中，mki为第k-1次更新的速度模型；Am为扰动模型；
[0018] 5)对目标函数进行计算获得地震波全波形反演的标准方程为：
[0019]
[0020] 式中，mn为模型参数；
[0021] 6)引入速度模型更新梯度方向和更新量步长，扰动模型Am表示为：
[0022] Am=agk；
[0023] 式中，gk为速度模型更新梯度方向，a为更新量步长；
[0024] 7)利用最小二乘的方法求解地震波全波形反演的速度模型更新梯度方向gk，其公 式为：
[0025]
[0026] 8)在速度模型更新梯度方向gkl对更新量步长a进行插值计算求取更新量步长 a;
[0027] 9)对速度模型进行更新得：
[0028] mk=mkl+agk;
[0029] 当|agk| <e时或速度模型更新次数达到速度模型更新迭代次数N时，速度模 型更新结束；否则，进入步骤3)。
[0030] 所述步骤3)中，波场误差矢量Ad的计算公式如下：
[0031] Ad=dobs-dcal；
[0032] 式中，(1^为地震观测波场数据，d为在迭代过程中的计算波场数据；
[0033] 计算波场数据(1。31的计算公式如下：
[0034]
[0035] 式中，浪为计算波场数据dMl和波场向量u之间的关系函数。
[0036] 所述步骤5)中，地震波全波形反演的标准方程具体计算过程包括以下步骤：
[0037] (1)对目标函数C(mk)进行二阶泰勒-拉格朗日展开，得到以下形式：
[0038]
[0039] 式中，j、l分别为二维数据的位置参数；M为正整数；m,、!^分别为更新模型变量mk: 在二维参数模型中某一个方向的分量；
[0040] (2)对步骤（1)中目标函数C(mkAm)的展开式以模型参数叫为变量进行求导： [0041 ]
[0042] 当目标函数C(mk)的导数;火/%为零时，目标函数C(叫）得到极值，此时，扰 动模型Am的表达式为：
[0043]
[0044] 扰动模型Am的表达式即为地震波全波形反演的标准方程。
[0045] 所述步骤8)中，更新量步长a满足伍尔夫条件：
[0046]
[0047]
[0048] 式中，cJPc2均为固定参数，其中ClG (〇，l)，c2G (Cl，l)，VC为目标函数C的 梯度。
[0049] 本发明由于采取以上技术方案，其具有以下优点：1、本发明由于采用迭代非线性 反演方法求解符合条件的模型，为深部大尺度构造演化分析、勘探地震成像及速度建模等 方面提供可靠依据。2、本发明由于采用在传统全波形反演的计算流程基础上，在每次的全 波形反演的非线性迭代过程中，引入最小二乘逆时偏移的流程，通过最小二乘的约束和多 次迭代修正，有效地提高了速度模型更新梯度准确性，可快速完成速度模型的更新。综上所 述，本发明可广泛应用于地震波全波形反演方法中。
【附图说明】
[0050] 图1是本发明的流程示意图。
【具体实施方式】
[0051] 下面结合附图和实施例对本发明进行详细的描述。
[0052] 本发明提供一种基于最小二乘梯度更新速度模型的地震波全波形反演方法，包括 以下步骤：
[0053] 1)在时间域内，地震波场的声波方程为：
[0054]
(1)
[0055] 其中，u是波场向量，s是震源向量，x和t分别是空间和时间坐标参数，M为质量 矩阵，A为刚度矩阵。
[0056] 2)构建初始速度模型，设置速度模型更新迭代次数N与速度模型允许的最小误差 值e〇
[0057] 3)构造观测波场数据和计算波场数据的波场误差矢量Ad，波场误差矢量Ad的 计算公式如下：
[0058] Ad=dobs-dcal (2)
[0059] 式中，(1^为地震观测波场数据，dMl为在迭代过程中的计算波场数据。
[0060] 其中，根据波场向量u计算波场数据cL1:
[006