测井信息约束的波形反演方法与流程

技术特征：

1.测井信息约束的波形反演方法，其特征是在正则化方法约束时间域全波形反演过程中，利用偏移成像结果中的大尺度地质构造，迭代控制测井数据的空间插值；

加入正则化约束项的目标泛函为：

$E\left(v\right)=\frac{1}{2}{\[u\left(v\right)-u_{obs}\]}^T\[u\left(v\right)-u_{obs}\]+\frac{1}{2}\mathrm{\λ}{(v-v_{prior})}^T{W}^{T}W(v-v_{prior})---\left(5\right)$

目标泛函的梯度公式为：

$\frac{\∂E}{\∂v}={\[\frac{\∂u\left(v\right)}{\∂v}\]}^T\[u\left(v\right)-u_{obs}\]+{\mathrm{\λW}}^{T}W(v-v_{prior})---\left(6\right)$

其中，E(v)为目标泛函，u(v)和u_obs分别是地震合成记录和地震实际观测数据，λ是正则化约束参数，v是反演速度，v_prior是先验速度，W是平滑矩阵，T代表矩阵转置，代表偏导数。

2.根据权利要求1所述的测井信息约束的波形反演方法，其特征是由于上式(6)中不好求取，利用伴随状态法间接得到，通过共轭梯度方法实现时间域的全波形反演计算；W平滑矩阵，取一阶或者二阶差分矩阵。

3.根据权利要求1所述的测井信息约束的波形反演方法，其特征是实现过程如下：

(1)根据给定一个初始速度模型，并对测井数据空间插值，构建准确的先验速度模型；

(2)利用交错网格有限差分方法合成地震记录，运用合成的地震记录与实际输入地震记录做差，计算剩余波场；

(3)建立测井数据约束正则化反演方程，求取速度更新量，更新速度模型，井位置处速度网格不更新；通过地震正演得到的正传波场和反向传播波场并在测井数据的约束下计算梯度方向；

(4)对梯度方向进行修正，变成共轭梯度方向，先给出试探步长，然后利 用线性搜索或者抛物线插值方法求解最佳步长，在时间域中运用共轭梯度的方法实现全波形的计算；

(5)通过迭代公式更新速度参数判断是否符合迭代收敛条件，如果不满足条件，当前得到的速度模型作为输入，返回步骤(2)，如果满足，迭代终止，输出结果。