[发明专利]一种沿地层走向光滑约束的反射波波形反演方法

说明书

本发明公开了一种沿地层走向光滑约束的反射波波形反演方法，包括以下步骤：选取初始背景速度模型；设置反射轴模型为零；计算临时反射轴模型；沿地层走向的光滑约束的更新背景速度；判断是否终止反演。本发明公开的沿地层走向光滑约束的反射波波形反演方法，针对目前无约束的反射波波形反演反演过程不稳定等问题，根据地质规律引入沿地层走向的光滑约束，利用反射波波形反演的临时反射轴图像识别地层界面，并利用地层界面得到了地层的走向，及通过解各向异性扩散方程来实现沿地层走向的光滑约束等手段来提高反演的稳定性，消除了柱状和球状的假异常，从而提高反演的准确性。

技术领域

本发明涉及油气勘探技术领域，具体涉及一种沿地层走向光滑约束的反射波波形反演方法。

背景技术

地震波是了解地球内部复杂的结构和构造、寻找矿产和油气资源的重要信息来源。 地震波包括了反射波、折射波等多种类型。然而，在现有的观测技术和条件下，只有 反射波才到达油气储层的深度。目前，波形反演是最为先进的一种成像方法，结合反 射波数据和波形反演技术的反射波波形反演，因此形成了一种成像深层油气藏的重要 方法。然而，由于在有限的地表观测范围内，能采集到的反射波对目标体的照明角度 有限，从而导致反射波波形反演稳定性差，给出的反演结果准确性差，不能满足表征 定量储层参数的要求，进而限制了反射波波形反演的应用。

发明内容

本发明提供一种沿地层走向光滑约束的反射波波形反演方法，用以解决目前的反射波波形反演的稳定性较差的问题。

本发明提供一种沿地层走向光滑约束的反射波波形反演方法，包括以下步骤：

步骤A：选取初始背景速度模型v_b，设置终止误差∈，最大迭代次数k，设立最小 化目标函数φ的表达式为：

式1中的u_δ表示预测的反射波波场，R表示用于从反射波波场中提取检波点的记录的检波器算子，d表示实际观测到的记录，v_b表示背景速度，m_δ表示反射轴模型；

步骤B：设置反射轴模型m_δ＝0；

步骤C：计算临时反射轴模型m_δ；

步骤D：沿地层走向光滑约束的更新背景速度v_b；

步骤E：判断是否终止反演：若完成步骤D后，目标函数的值仍然大于∈并且迭 代次数小于k，重复步骤B至步骤D，直至满足终止条件。

进一步地，所述步骤C包括以下步骤：

步骤C1：求最小化目标函数关于反射轴模型m_δ的梯度；

步骤C2：根据步骤C1中梯度通过局部梯度法对模型m_δ进行更新。

进一步地，所述步骤D包括以下步骤：

步骤D1：求最小化目标函数关于背景速度v_b的梯度；

步骤D2：根据步骤D1中梯度通过局部梯度法对模型v_b进行更新；

步骤D3：背景速度沿地层光滑约束。

进一步地，所述步骤D3包括以下步骤：

步骤D31：通过步骤C产生的临时反射轴模型m_δ求解该模型图像的梯度

步骤D32：由梯度的三个分量形成张量T，所述张量T的表达式为：