[发明专利]一种混合全局优化算法的叠前地震多参数反演方法

摘要

本发明公开了一种混合全局优化算法的叠前地震多参数反演方法，该方法将粒子群算法和快速模拟退火算法有效结合，解决传统粒子群算法易不成熟收敛的问题，并在粒子群算法中添加基于三参数联合概率密度择优组合的多维学习项，克服叠前地震多参数同步反演的不稳定，可以同步而准确的获取纵波速度、横波速度和密度三参数反演结果。使用本发明方法对合成地震资料进行反演取得满意的结果，与使用传统粒子群算法的反演结果相比反演效果改善明显。

主权项

1.一种混合全局优化算法的叠前地震多参数反演方法，其特征在于，包括如下步骤：步骤1，确定粒子待反演参数的取值范围，设定粒子的数量、最大迭代次数、模拟退火冷却进度表、模拟退火学习项加权系数和多维学习项加权系数；步骤2，设置各个粒子的初始位置和初始速度，位置包括四个维度：纵波速度、横波速度、密度以及前述三个参数的联合概率密度；步骤3，反演第k次迭代时，根据粒子的位置计算粒子适应度，其中粒子适应度由观测地震记录与合成地震记录的误差和步骤2所述三个参数的先验约束项构成；所述粒子适应度的表达式为：其中，为第i个粒子适应度，P_i为第i个粒子位置，E(P_i)为P_i的期望，L为观测数据采样长度，θ为入射角，t为地震记录采样时间，W为震源子波，R为纵波反射系数，D为观测地震记录，σ为纵波速度V_Pi、横波速度V_Si和密度ρ_i这三个参数的协方差矩阵，λ₁、λ₂均为预设定系数，上标T表示矩阵的转置；步骤4，根据步骤3的粒子适应度，计算各粒子被选择作为模拟退火学习项引导粒子的接收概率，并基于各粒子的接收概率，通过轮盘选择法则确定模拟退火学习项引导粒子；步骤5，固定每个粒子的横波速度和密度，用所有粒子的纵波速度依次替换每个粒子的纵波速度，并利用联合概率密度公式计算所有粒子的纵波速度对应的联合概率密度，并找出最大联合概率密度；从每个粒子对应的最大联合概率密度中找出最大值及最大值所对应的纵波速度VPmax；按同样的方法，找到每个粒子对应的最大联合概率密度中最大值对应的横波速度VSmax、密度ρmax；根据VPmax、VSmax和ρmax，得到多维学习项引导粒子；步骤6，根据模拟退火学习项引导粒子和多维学习项引导粒子更新各粒子的速度和位置；步骤7，进行第k+1次迭代，重复步骤3至步骤6，直至达到模拟退火冷却进度表中的终止退火温度和最大迭代次数后，输出粒子适应度最优的粒子。