[发明专利]一种介质品质因子的逐次差分进化估计方法

摘要

本发明一种介质品质因子的逐次差分进化估计方法，包括如下步骤，步骤1，先对实际VSP道集进行预处理，去除不同炮子波的不一致性；步骤2，构建一个含M个个体的初始群体X；步骤3，按相邻两个检波器的距离将地下介质划分为不同地层，将每一道作为一个子成分，用L₁范数定义第k个个体的局部适应度函数，从而得到每道的局部适应度值；步骤4，根据实际VSP道集和每道的局部适应度，对初始群体X中的每个子成分进行优化更新；根据变异策略求取变异后的待估品质因子和速度；步骤5，重复步骤4，将整个群体更新完后根据全局适应度的值选择最优个体，从而得到地下介质的品质因子和速度。其整体的计算速度快，收敛性好，且不依赖于初值。

主权项

一种介质品质因子的逐次差分进化估计方法，其特征在于，包括如下步骤，步骤1，用WEPIF方法估计实际VSP道集的品质因子和速度前，先对实际VSP道集进行预处理，去除不同炮子波的不一致性；步骤2，构建一个含M个个体的初始群体X，表示为，X＝[x₁,x₂,…,x_M]；群体中的第k个个体表示为，x_k＝[x_k,1,x_k,2,…,x_k,N],k＝1,…,M；其中，N为总道数，M为地层数量，每个地层对应一组要进行估计的品质因子和速度；步骤3，按相邻两个检波器的距离将地下介质划分为不同地层，将每一道作为一个子成分，用L₁范数定义第k个个体的局部适应度函数为，$f_{k,j}=\underset{n=P-P_1}{\overset{P+P_2}{\Σ}}|\[d_{o,j}\left(n\right)-d_{k,j}\left(n\right)\]\·g\left(n\right)|;$从而得到每道的局部适应度值；式中，d_o,j(·)为第j道观测数据，d_k,j(·)为由第k个个体正演计算的第j道数据，P为观测数据直达波瞬时振幅谱的峰值位置，P₁、P₂分别为子波瞬时振幅的包络峰值左右两侧的第一个波谷位置，g(n)是高斯窗函数；步骤4，根据实际VSP道集和每道的局部适应度，对初始群体X中的每个子成分进行优化更新；进行正演模拟逐道的优化每个子成分，当所有个体中当前层和下一层的参数均选择最优值后，再进行层位循环，逐层反演参数；反演第j层的参数时，将群体X中第j列的所有变量进行升序排列，根据变异策略求取变异后的待估品质因子和速度；步骤5，重复步骤4，将整个群体更新完后根据全局适应度的值选择最优个体，从而得到地下介质的品质因子和速度。