[发明专利]一种增强反Q值补偿稳定性的地震资料处理方法及系统

权利要求书

1.一种增强反Q值补偿稳定性的地震资料处理方法，其特征在于，包括：

步骤1：获取地震数据、反射系数和Q值之间的关系表达式；

步骤2：获得Lp范数约束下的反演目标函数；

步骤3：获得所述Lp范数约束下的反演目标函数的Q值；

步骤4：基于所述关系表达式和所述Lp范数约束下的反演目标函数的Q值，获得反射系数；

所述步骤2包括：

子步骤201：在稳态褶积模型中结合地层正Q滤波算子，推导非稳态数据时变褶积模型，建立时变子波与Q值的关系模型，并基于所述Q值建立时变子波库；

子步骤202：基于所述时变子波库，建立L_P范数约束下的反演目标函数；

所述Lp范数约束下的反演目标函数为：

其中，s为地震数据，r为反射系数，p为范数，G＝F^-1WA为时变子波矩阵，为傅里叶逆变换矩阵，

减矩阵，为衰减值，t₀为初始时间，ω₀为参考角频率，ω_d为d点处的角频率，γ为衰减项，Q为介质的品质因子，W为初始子波矩阵。

2.根据权利要求1所述的地震资料处理方法，其特征在于，所述步骤1包括：

子步骤101：基于地震数据建立波动方程，求取所述波动方程的平面波解析解；

子步骤102：将所述平面波解析解中的实速度替换为复速度，获取衰减介质下的平面波解析解；

子步骤103：对所述衰减介质下的平面波解析解进行离散化，获得所述地震数据、反射系数和Q值之间的关系表达式。

3.根据权利要求2所述的地震资料处理方法，其特征在于，所述衰减介质下的平面波解析解为：

其中，ω为角频率，r(t′)为反射系数，t′为传播时间，为地震子波，t为时间，为振幅衰减项，为相位校正项或子波整形项，ω₀为参考角频率，Q(t′)为介质的品质因子，γ为衰减项，

4.根据权利要求1所述的地震资料处理方法，其特征在于，所述地震数据、Q值和反射系数之间的关系表达式为：

s＝F^HWA(Q)r (2)

其中，为傅里叶算子，H为复数，t为时间，f为频率，为初始子波矩阵，w(f₁)为初始子波，为衰减矩阵，α为衰减值，s为地震数据，r为反射系数。

5.根据权利要求1所述的地震资料处理方法，其特征在于，所述L_P范数中的P的数值范围为0-1。

6.一种增强反Q值补偿稳定性的地震资料处理系统，其特征在于，该系统包括：

存储器，存储有计算机可执行指令；

处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：

步骤1：获取地震数据、反射系数和Q值之间的关系表达式；

步骤2：获得Lp范数约束下的反演目标函数；

步骤3：获得所述Lp范数约束下的反演目标函数的Q值；

步骤4：基于所述关系表达式和所述Lp范数约束下的反演目标函数的Q值，获得反射系数；

所述步骤2包括：

子步骤201：在稳态褶积模型中结合地层正Q滤波算子，推导非稳态数据时变褶积模型，建立时变子波与Q值的关系模型，并基于所述Q值建立时变子波库；

子步骤202：基于所述时变子波库，建立L_P范数约束下的反演目标函数；

所述Lp范数约束下的反演目标函数为：

其中，s为地震数据，r为反射系数，p为范数，G＝F^-1WA为时变子波矩阵，为傅里叶逆变换矩阵，

为衰减矩阵，为衰减值，t₀为初始时间，ω₀为参考角频率，ω_d为d点处的角频率，γ为衰减项，Q为介质的品质因子，W为初始子波矩阵。

7.根据权利要求6所述的地震资料处理系统，其特征在于，所述步骤1包括：

子步骤101：基于地震数据建立波动方程，求取所述波动方程的平面波解析解；

子步骤102：将所述平面波解析解中的实速度替换为复速度，获取衰减介质下的平面波解析解；

子步骤103：对所述衰减介质下的平面波解析解进行离散化，获得所述地震数据、反射系数和Q值之间的关系表达式。