[发明专利]智能化时变盲反褶积宽频处理方法及装置

权利要求书

1.一种智能化时变盲反褶积宽频处理方法，其特征在于，包括：

从衰减地震数据中提取出初始地震子波；

随机产生包含多个地层品质因子Q值的Q值群，包括：

随机产生每个地层品质因子Q值的初始位置和初始速度；

确定所述Q值群的初始全局最好位置，每个地层品质因子Q值的初始历史最好位置；

根据所述初始地震子波和所述Q值群建立时变子波矩阵群；

根据所述衰减地震数据和所述时变子波矩阵群，获得多个反射系数；

建立适应度函数；

根据所述适应度函数和所述多个反射系数，获得多个适应值，其中，一个适应值对应一个地层品质因子Q值；

根据所述多个适应值确定最优地层品质因子Q值。

2.如权利要求1所述的智能化时变盲反褶积宽频处理方法，其特征在于，根据所述衰减地震数据和所述时变子波矩阵群，获得多个反射系数，包括：

根据所述衰减地震数据和所述时变子波矩阵群，运用稀疏贝叶斯反演算法进行时变盲反褶积处理，获得多个反射系数。

3.如权利要求1所述的智能化时变盲反褶积宽频处理方法，其特征在于，建立适应度函数，包括：

选择l_0.1范数作为评价准则建立适应度函数。

4.如权利要求1所述的智能化时变盲反褶积宽频处理方法，其特征在于，根据所述多个适应值确定最优地层品质因子Q值，包括：

将所述多个适应值两两之间进行比较，获得最小适应值；

将所述最小适应值与预设适应值阈值进行比较，当最小适应值小于预设适应值阈值时，所述最小适应值对应的地层品质因子Q值为最优地层品质因子Q值。

5.如权利要求4所述的智能化时变盲反褶积宽频处理方法，其特征在于，还包括：

当最小适应值大于预设适应值阈值时，更新所述Q值群；

根据更新后的所述Q值群重新确定最优地层品质因子Q值。

6.如权利要求5所述的智能化时变盲反褶积宽频处理方法，其特征在于，当最小适应值大于预设适应值阈值时，更新所述Q值群，包括：

根据所述每个地层品质因子Q值的初始位置和初始速度、所述Q值群的初始全局最好位置、每个地层品质因子Q值的初始历史最好位置，更新所述Q值群。

7.如权利要求5或6所述的智能化时变盲反褶积宽频处理方法，其特征在于，还包括：

当更新次数达到预设次数上限和/或最小适应值小于预设适应值阈值时，结束更新。

8.如权利要求1所述的智能化时变盲反褶积宽频处理方法，其特征在于，根据初始地震子波和地层品质因子Q值群建立时变子波矩阵群，包括：

G＝F^HWA(Q)；

其中，G表示时变子波矩阵群；F表示傅里叶算子；W表示初始地震子波矩阵；A(Q)表示依赖于Q值的时变衰减矩阵；^H表示转置共轭。