[发明专利]一种提高高陡构造速度反演效率的全波形联合反演方法

说明书

技术领域

本发明属于石油地球物理勘探领域，具体涉及一种提高高陡构造速度反演效率的全波形联合反演方法。

背景技术

常规一次反射波勘探，受观测系统限制，对断层、盐丘侧翼等陡倾角地下构造照明不足，高陡构造精确成像困难。高陡构造产生的棱柱波信息，能够有效增加高陡构造的照明程度，改善高陡构造成像效果。因此充分利用棱柱波信息对于提高高陡构造的成像能力是非常必要的。

由于全波形反演高精度以及高分辨率的特点，使其成为速度建模的一种有力工具，逐渐成为研究的热点。全波形反演是一个非线性数据拟合的过程，通过减少观测数据与预测数据的之间的差值来更新参数模型，这个过程以迭代的方式重复下去，直到数据差值足够小为止。全波形反演计算量很大，计算效率问题成为全波形反演的研究重点之一。

发明内容

针对现有技术中存在的上述技术问题，本发明提出了一种提高高陡构造速度反演效率的全波形联合反演方法，能够降低对初始模型中高陡构造信息是否缺失的依赖性，提高对高陡构造的反演精度和效率。

为了实现上述目的，本发明采用如下技术方案：

一种提高高陡构造速度反演效率的全波形联合反演方法，包括如下步骤：

步骤1：输入初始速度场的速度、野外观测炮记录及震源子波，并建立观测系统；

步骤2：使用下式表示的成像公式获得成像结果：

$I=\underset{x_s}{\Σ}\underset{t}{\Σ}uw/u^2$

其中，I表示成像结果，u和w分别表示正传波场和反传波场，w＝L^*(R^*(p_obs))，L^*为反传算子，R^*表示将限定在检波器上的数据残差空间扩展到整个模型空间，p_obs为野外观测炮记录；

步骤3：将步骤2中的成像结果作为反射系数；

步骤4：应用线性正演模拟和线性波场反传，分别计算两种棱柱波prism1和prism2的梯度方向；

$g{\left(v\right)}_{prism1}=\frac{2}{v^3}\underset{x_s}{\Σ}\underset{t}{\Σ}\frac{{\∂}^{2}B\left(u\right)}{\∂t^2}{w}^{\′}$