[发明专利]面波频散曲线的反演方法、装置及计算机可读存储介质

权利要求书

1.一种面波频散曲线的反演方法，其特征在于，包括：

获取构建的第一地质模型集中每个第一地质模型的面波频散曲线与待反演的面波频散曲线的拟合误差，其中，所述拟合误差用于表征第一地质模型的面波频散曲线与所述待反演的面波频散曲线的适合度，且拟合误差越小，适合度越高，拟合误差越大，适合度越低；

b.判断是否存在有任一拟合误差满足预设条件；

c.若否，则在第二地质模型集中，确定每个第二地质模型的所属群体，其中，所述每个第二地质模型的所属群体包括在对应第二地质模型相邻半径内的第二地质模型，所述相邻半径是根据对应第二地质模型的搜索空间的半径得到的，且所述第二地质模型集为剔除了拟合误差最大以及最小的第一地质模型后的第一地质模型集；

d.对每个第二地质模型的所属群体，进行行为权重的调整，以通过调整后的行为权重，更新每个第二地质模型在搜索空间中的飞行速度以及位置，其中，所述行为权重用于表征对对应第二地质模型在搜索空间中运动轨迹的影响程度，且所述行为权重包括聚集权重、避撞权重以及结队权重；

计算更新了飞行速度以及位置后的每个第二地质模型的面波频散曲线与所述待反演的面波频散曲线的拟合误差，并重复前述步骤b～d，直至有任一拟合误差满足预设条件为止；

将满足预设条件的任一拟合误差对应的第二地质模型，作为与所述待反演的面波频散曲线最匹配的地质模型；

对每个第二地质模型的所属群体，进行行为权重的调整，包括：

按照如下公式，调整每个第二地质模型所属群体的行为权重；

式中，和分别表示第i个第二地质模型的所属群体调整后的聚集权重、避撞权重以及结队权重，c_i、s_i以及a_i分别表示第i个第二地质模型的所属群体的初始聚集权重、初始避撞权重以及初始结队权重，b₁和b₂分别为常数；m为第i个第二地质模型的所属群体中包含的第二地质模型的个数，n为第二地质模型的总数量，rankfit1(i)表示将第i个第二地质模型的所属群体内所有第二地质模型的拟合误差从大到小进行排序后，第i个第二地质模型在其中的排名，rankfit2(i)表示在当前步骤中，将所有第二地质模型的拟合误差从大到小进行排序后，第i个第二地质模型在其中的排名。

2.如权利要求1所述的方法，其特征在于，通过调整后的行为权重，更新每个第二地质模型在搜索空间中的飞行速度以及位置，包括：

按照如下公式更新第二地质模型在搜索空间中的飞行速度；

式中，V_i^t+1为第i个第二地质模型更新后的飞行速度，t代表迭代次数，S_i、A_i、C_i、F_i以及E_i分别为第i个第二地质模型的避撞行为向量、结队行为向量、聚集行为向量、觅食行为向量以及避敌行为向量；f_i和e_i分别为第i个第二地质模型的觅食行为权重以及避敌行为权重，V_i^t代表第i个第二地质模型在第t次迭代时的飞行速度，w_i为第i个第二地质模型的飞行速度的权重；

按照如下公式更新第二地质模型在搜索空间中的位置；

式中，为第i个第二地质模型更新后的位置，为第i个第二地质模型第t次迭代时的位置。

3.如权利要求1所述的方法，其特征在于，在第二地质模型集中，确定每个第二地质模型的所属群体前，所述方法还包括：

标记拟合误差最小的第一地质模型以及拟合误差最大的第一地质模型在搜索空间中的位置；

判断拟合误差最小的第一地质模型的位置在连续k次迭代中是否未发生变化，其中，k为正整数；

若是，则对所述第一地质模型集进行群体划分，以便在群体划分后，调整群体的避撞权重。