[发明专利]一种地震叠前反演的方法及装置

说明书

本发明提供一种地震叠前反演的方法及装置，所述方法包括：通过测井数据正演模拟，确定页岩气地震数据中噪声信号的噪声类型，其中，所述噪声类型至少包括以下之一：高斯、亚高斯和超高斯；在第一目标函数中加入L2范数和L4范数；基于设定规则自适应调节所述噪声信号在所述目标函数中L2范数和L4范数之间的权重参数；根据所述权重参数进行叠前同步反演确定反演结果。反演结果能够较好的展示出明显的碳酸盐岩与页岩分界面，岩性区分清楚，反演结果与地质情况吻合度较高。

技术领域

本发明涉及地震反演技术领域，尤其涉及一种地震叠前反演的方法及装置。

背景技术

在求解地震反演问题的过程中，通常需要加入一定的先验信息。这种先验信息有的是为了减少反演的不适定性，有的是体现了地质学家的先验知识，有的是为了提高反演算法对地震数据噪声的压制效果等。

然而，现有的反演方法，大都假设噪声服从高斯分布。但由于地震数据噪声来源的多样性，噪声很难严格满足高斯分布的先验假设，噪声先验假设的不合理性会影响反演结果的精度。

发明内容

本发明实施例提供一种地震叠前反演的方法及装置，用以解决现有的反演方法，反演结果的精度不高的问题。

第一方面，本发明实施例提供一种地震叠前反演的方法，包括：

通过测井数据正演模拟，确定页岩气地震数据中噪声信号的噪声类型，其中，所述噪声类型至少包括以下之一：高斯、亚高斯和超高斯；

在第一目标函数中加入L2范数和L4范数；

基于设定规则自适应调节所述噪声信号在所述目标函数中L2范数和L4范数之间的权重参数；

根据所述权重参数进行叠前同步反演确定反演结果。

在一可能的实施例中，所述确定页岩气地震数据中噪声的噪声类型，包括：

根据峭度确定噪声信号的噪声类型。

在一可能的实施例中，所述峭度为：

其中，为均值，σ_i为标准差，N为信号采样长度。

在一可能的实施例中，所述设定规则包括：

所述权重参数为反演误差的函数；

或，

当噪声信号的类型为高斯或是超高斯时，权重参数Υ(n_b)≈0，采用L2范数；当噪声信号的类型为亚高斯时，权重参数Υ(n_b)≈1，采用L4范数；

其中，所述权重参数Υ(n_b)与所述峭度χ(n_b)之间符合如下公式：

在一可能的实施例中，所述根据所述权重参数进行叠前同步反演确定反演结果，包括：

根据所述权重参数确定第二目标函数；

基于第二目标函数进行叠前同步反演确定反演结果；

所述第二目标函数为：

其中，向量d为M×1维向量，表示地震叠前道集；G为M×N维矩阵，表示正演模型矩阵；m为N×1维向量。

在一可能的实施例中，所述方法还包括：

采用迭代残差n_k代替真实的噪声数据误差n_b，重新确定权重参数Υ(n_b)；