[发明专利]一种地震资料的反演方法

说明书

本发明提供一种地震资料的反演方法。基于Lagrange乘子法和Tikhonov正则化方法构造光滑泛函；建立双参数模型函数根据双滤子正则化算法计算双参数模型函数中的参数t_n和b_n；使用上述双参数模型函数逼近光滑泛函；根据修正的广义偏差原理对上述双参数模型函数进行求解；循环计算直至参数满足预设条件。该方法能够兼顾数据噪音和子波算子本身估计不准的影响。

技术领域

本发明涉领域地震勘探领域，具体涉及一种应用于地震波弹性阻抗反演中的双参数正则化方法。

背景技术

传统的地震波阻抗反演是基于P波在地下界面垂直入社的假设。但是，当AOV效应存在时，远偏移距资料反演得到的反射系数信息将明显不同于垂直入射情况，因此，弹性阻抗反演问题，需要充分提取叠钱地震资料中丰富的信息，将远偏移距资料与测井数据、岩心测量数据等直接相关联，有效降低声波阻抗油气检测的多解性，使得反演结果更为准确。

基于卷积模型的反演方法，直接求解离散线性算子方程表达式相关的最小二乘问题会来带不稳定的计算效果，需要使用正则化方法并辅之以适当的最优化技巧。另一方面，由于地震资料数据量大，信噪比差，资料误差水平难以准确估计，弹性阻抗反演需要借助后验正则化参数选取准则有效的确定最优化正则参数。

基于弹性阻抗反演原理，建立的Tikhonov正则化反演模型，在实际生产中，采用Tikhonov正则化方法求解弹性阻抗问题时，估算的子波可能不准确，即算子存在扰动的情况。

发明内容

为解决上述技术问题，本发明提出了一种应用于地震波弹性阻抗反演中的双参数正则化方法。与传统正则化方法的不同之处在于，传统方法在估计正则参数时只考虑数据噪音的影响，本发明兼顾了数据噪音和子波算子本身估计不准的影响，提出双参数正则化方法。

该方法具体步骤如下：

(1)基于Lagrange乘子法和Tikhonov正则化方法构造光滑泛函；

(2)建立双参数模型函数其中α是正则参数；

(3)根据双滤子正则化算法计算双参数模型函数中的参数t和b；

(4)使用上述双参数模型函数逼近所述光滑泛函；根据修正的广义偏差原理对上述双参数模型函数进行求解；

(5)循环计算直至参数满足预设条件。

其中，步骤(1)中构造光滑泛函的具体步骤如下：

构造光滑泛函Jα[R，δ，h]＝||W_hR-d_δ||²+αR||²，并考虑其极值问题令η＝(δ，h)。

按照正则化方法的思想，构造函数α＝α(η)，使得当η→0时，由上式求得的极小值满足R*为原问题的精确解。

记

令M_η(α)＝Jα[R_η(α)，η]＝||W_hR_η(α)-d_δ||²+α||R_η(α)²||(1)，

记θ_η(α)＝||R_η(α)||²，β_η(α)＝||W_hR_η(α)-d_δ||²，其中，R是反射系数，W是地震子波算子，h是地震子波算子的扰动误差，W_h是带有扰动误差的地震子波算子，d是地震波数据，δ是数据中的干扰噪音，d_δ是存在噪音干扰的地震数据，α是正则参数。