[发明专利]基于多信息约束反演的品质因子确定方法、装置及系统

权利要求书

1.基于多信息约束反演的品质因子确定方法，其特征在于，包括：

获取叠前地震数据，所述叠前地震数据包括时空几何结构信息、时空域振幅衰减信息、时空域相位衰减信息；

基于岩石物理理论中子波衰减的物理信息，建立叠前地震数据、品质因子、反射系数之间的关系模型；

基于所述关系模型，反演获得多道反射系数；

利用预设稀疏准则函数计算所述多道反射系数的稀疏度；

根据固有衰减与反射系数稀疏度的关系，确定品质因子。

2.如权利要求1所述的方法，其特征在于，所述基于所述关系模型，反演获得多道反射系数，包括：

根据预设衰减模型以及时变子波与时不变子波关系，获得时变地震子波；

基于时变地震子波与反动校正算子，获得时变混合反动校正算子；

基于所述关系模型、时变混合反动校正算子以及叠前地震数据，获得反射系数。

3.如权利要求2所述的方法，其特征在于，所述基于所述关系模型、时变混合反动校正算子以及叠前地震数据，获得反射系数，包括：

利用交替方向法对公式求解，获得反射系数；

其中，表示R_i的l_2,1范数值的最小值，表示地震数据匹配项，σ表示数据误差或噪声能量，S表示未经校正的地震数据，i表示扫描的品质因子的个数，G_i表示第i个品质因子对应的时变混合反动校正算子，R_i表示第i个品质因子对应的反射系数矩阵，下标F表示弗罗贝尼乌斯范数。

4.如权利要求1所述的方法，其特征在于，所述利用预设准则函数计算所述多道反射系数的稀疏度，包括：

利用l_2,p范数计算所述多道反射系数的稀疏度；或者，

选取窗函数，利用所述窗函数对多道反射系数进行滑动开窗处理；

基于l_2,p范数计算窗内准则函数值，将所述窗内准则函数值作为所述多道反射系数的稀疏度。

5.如权利要求4所述的方法，其特征在于，所述根据固有衰减与反射系数稀疏度的关系，确定品质因子，包括：

当利用l_2,p范数计算所述多道反射系数的稀疏度时，将所述稀疏度中最小值对应的固有衰减确定为品质因子；

当将窗内准则函数值作为所述多道反射系数的稀疏度时，基于所述准则函数值构建品质因子谱图，将所述品质因子谱图中最大值对应的固有衰减确定为品质因子。

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

基于所述品质因子，对地震波的吸收衰减进行补偿处理。

7.基于多信息约束反演的品质因子确定装置，其特征在于，包括：

地震数据获取模块，用于获取叠前地震数据，所述叠前地震数据包括时空几何结构信息、时空域振幅衰减信息、时空域相位衰减信息；

关系模型建立模块，用于基于岩石物理理论中子波衰减的物理信息，建立叠前地震数据、品质因子、反射系数之间的关系模型；

反射系数获得模块，用于基于所述关系模型，反演获得多道反射系数；

稀疏度计算模块，用于利用预设稀疏准则函数计算所述多道反射系数的稀疏度；

品质因子确定模块，用于根据固有衰减与反射系数稀疏度的关系，确定品质因子。

8.如权利要求7所述的装置，其特征在于，所述反射系数获得模块，包括：

时变地震子波获得单元，用于根据预设衰减模型以及时变子波与时不变子波关系，获得时变地震子波；

时变混合反动校正算子获得单元，用于基于时变地震子波与反动校正算子，获得时变混合反动校正算子；

反射系数获得单元，用于基于所述关系模型、时变混合反动校正算子以及叠前地震数据，获得反射系数。