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

权利要求书

1.一种地震叠前反演的方法，其特征在于，包括：

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

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

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

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

2.根据权利要求1所述的方法，其特征在于，所述确定页岩气地震数据中噪声的噪声类型，包括：

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

3.根据权利要求2所述的方法，其特征在于，所述峭度为：

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

4.根据权利要求1所述的方法，其特征在于，所述设定规则包括：

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

或，

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

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

5.根据权利要求1所述的方法，其特征在于，所述根据所述权重参数进行叠前同步反演确定反演结果，包括：

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

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

所述第二目标函数为：

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

6.根据权利要求5所述的方法，其特征在于，所述方法还包括：

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

m_k＝m_k-1+α_k-1P_k-1

n_k＝d-Gm_k＝ΔGm_k+GΔm+ΔGΔm+n_b

其中，m_k表示第k次迭代结果，P_k-1表示第k次迭代搜索方向，α_k-1表示搜索步长，ΔG表示真实正演模型与使用的有偏正演模型之间的模型误差，Δm表示当前迭代模型m_k与真实模型参数m之间的反演误差；

所述第二目标函数还为：

7.一种地震叠前反演的装置，其特征在于，包括：

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

添加模块，用于在第一目标函数中加入L2范数和L4范数；

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

第三确定模块，用于根据所述权重参数进行叠前同步反演确定反演结果。

8.根据权利要求7所述的装置，其特征在于，所述第一确定模块，用于根据峭度确定所述噪声信号的噪声类型。

9.根据权利要求8所述的装置，其特征在于，所述峭度为：

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