[发明专利]一种速度场建模方法

说明书

本申请实施例公开了一种速度场建模方法，所述方法包括：S1，确定多个射线参数和多个参考炮点；S2，根据观测数据和所述多个射线参数、多个参考炮点获取平面波数据；S3，确定目标函数和所述多个射线参数的计算顺序，其中，所述目标函数与所述平面波数据相关；S4，根据所述目标函数和所述多个射线参数的计算顺序，对初始模型进行迭代更新，生成速度场模型。本申请实施例通过选择多个参考炮点和射线参数，获得多个平面波数据，有效的提高了观测数据对地下模型的覆盖次数，提高了模型的准确度。此外，通过对射线参数按照绝对值由高到低逐次迭代反演，反演过程更为稳定，生成的速度场模型更为准确。

技术领域

本申请涉及地震资料处理方法，尤其涉及一种速度场建模方法。

背景技术

当前技术中，在对地下构造反演时，较为先进的是全波形反演方法。全波形反演是一个极强的非线性问题，其反演能力在很大程度上依赖于初始模型的准确性。在实际应用中，初始模型往往是通过速度分析或走时层析获得的，其精度有限。

现阶段低频数据采集能力较弱，以及地震波对地下介质照明强度较低(照明强度即地震波对模型的覆盖次数，覆盖次数越多，则照明强度越强)。低频成分的缺失导致初始模型的中低波数成分往往难以恢复，造成波数成分的不连续，进而影响了全波形反演的稳定性。而地震波照明强度低导致波形反演的非线性和多解性增强，导致反演难度增大。

基于此，我们需要一种更准确的速度场模型作为全波形反演的初始模型。

发明内容

本申请实施例提供一种速度场建模方法，用于解决如下问题：为全波形反演提供更准确的速度场模型。

基于此，本申请实施例提供一种速度场建模方法，所述方法包括：

S1，确定多个射线参数和多个参考炮点；

S2，根据观测数据和所述多个射线参数、多个参考炮点获取平面波数据；

S3，确定目标函数和所述多个射线参数的计算顺序，其中，所述目标函数与所述平面波数据相关；

S4，根据所述目标函数和所述多个射线参数的计算顺序，对初始模型进行迭代更新，生成速度场模型.

相较于现有技术而言，本申请实施例通过选择多个参考炮点，每个射线参数获得多个平面波数据，有效的提高了观测中对地下模型的覆盖次数，生成的速度场模型更为准确。此外，通过对射线参数按照绝对值由高到低逐次迭代反演，实现了对地下介质由低波数到高波数的多尺度反演，生成模型的过程更为稳定；利用平面波合成中的叠加效应，具有压制噪音的作用，对含噪音数据具有较高的实用性；利用波数与散射角的关系，降低了波形反演对低频数据的依赖性。

附图说明

图1为本实施例中反演时波数与平面波入射角之间的关系示意图；

图2为本实施例中不同射线参数时得到的反演结果，其中：(a)为数值实验采用的真实模型，(b)为入射角参数为0.8时的反演结果；(c)为射线参数为0.4时的反演结果；(d)为以(c)为初始模型的常规全波形反演结果。

图3为本发明实施例中当低频数据缺失时反演所得结果，其中：(a)为高通滤波后雷克子波的频谱；(b)为高通滤波后雷克子波的时间域波形；(c)为低频信息缺失情况下的反演结果；(d)为以(c)为初始模型的常规全波形反演在低频缺失下的反演结果；

图4为本发明实施例中添加噪音数据时反演所得结果：(a)为含噪音数据情况下的反演结果；(b)为以(a)为初始模型的全波形反演结果。；

图5为本发明实施例中地震数据随机缺少50炮的反演结果：(a)为数据缺失情况下的反演结果；(b)为以(a)为初始模型的全波形反演结果；