应用闪电连接过程优化算法反演多阶瑞雷波频散曲线.pdf

2023 年 8 月第 58 卷 第 4 期应用闪电连接过程优化算法反演多阶瑞雷波频散曲线付宇1，2，艾寒冰1，2，姚振岸*1，2，李红星1，2，田宵1，2，张杏棉1，2（1.江西省防震减灾与工程地质灾害探测工程研究中心（东华理工大学），江西南昌 330013；2.东华理工大学地球物理与测控技术学院，江西南昌 330013）摘要：在瑞雷波勘探中，频散曲线反演是获取地下地层信息的关键，具有多参数、多极值特征。传统的全局优化算法，如粒子群（PSO）算法、遗传算法（GA）等，进行频散曲线反演时存在收敛速度慢、求解精度低等问题。为此，将一种新颖的全局优化算法闪电连接过程优化（Lightning Attachment Procedure Optimization，LAPO）算法用于频散曲线反演。LAPO 算法受自然界中闪电下行先导和上迎先导连接过程的启发所提出，通过模拟自然界中闪电形成过程寻求最优解。该算法全局搜索和局部开发能力都强，无需进行参数调整。使用理论地质模型和实际数据进行测试，结果表明：与 PSO 算法相比，LAPO 算法具有收敛速度快、求解精度高、算法性能稳定的特点；LAPO 算法不仅适用于基阶频散曲线反演，而且适用于多阶频散曲线反演，是瑞雷波频散曲线定量解释的有效工具。关键词：瑞雷波，频散曲线反演，多阶，全局优化，闪电连接过程优化算法中图分类号：P631 文献标志码：A doi：10.13810/ki.issn.1000-7210.2023.04.008Inversion of multimode Rayleigh wave dispersion curves based on lightning attachment procedure optimizationFU Yu1，2，AI Hanbing1，2，YAO Zhen an1，2，LI Hongxing1，2，TIAN Xiao1，2，ZHANG Xingmian1，2（1.Engineering Research Center for Earthquake Disaster Prevention and Engineering Geological Disaster Detection of Jiangxi Province（East China University of Technology），Nanchang，Jiangxi 330013，China；2.School of Geophysics and Measurementcontrol Technology，East China University of Technology，Nanchang，Jiangxi 330013，China）Abstract：In Rayleigh wave exploration，inversion of picked dispersion curves is the key to retrieving subsurface stratigraphic information，which is characterized by multiple parameters and optimal values.Traditional global optimization methods such as particle swarm optimization（PSO）algorithm and genetic algorithm（GA）feature slow convergence speed and low precision.To this end，this paper adopts a novel global optimization algorithm，lightning attachment procedure optimization（LAPO）for dispersion curve inversion.The LAPO algorithm is inspired by the natural lightning formation process involving the process of connecting downward and upward leader movements.This algorithm can find the optimal solution by simulating the lightning formation process in nature and has strong global search and local exploitation abilities without parameter tuning.Tests carried out through the theoretical geological model and actual data show that compared with PSO，LAPO yields faster convergence speed，higher solution accuracy，and more stable performance.Additionally，it can be applied to the inversion of both fundamental and multimode dispersion curves，and is a valuable tool for the quantitative interpretation of Rayleigh dispersion curves.处理技术 文章编号：1000-7210（2023）04-0830-09*江西省南昌市经开区广兰大道 418号东华理工大学地球物理与测控技术学院，330013。Email：本文于 2022年 7月 14日收到，最终修改稿于 2023年 5月 9日收到。本项研究受国家自然基金项目“黏弹各向异性介质复杂震源机制解析与微地震响应高精度正演”（42004113）、江西省自然科学基金项目“复杂介质微地震响应无网格高精正演与震源机制解析”（20212BAB211003）、“基于人工智能的江西地区天然地震和非天然地震事件识别方法研究”（2022AB213047）、江西省大学生创新创业项目“探地雷达在植物根干探测中的应用研究”（S202210405012）、“基于机器学习的瑞雷面波反演方法研究”（S202210405013）、“南昌市地下空间微动探测方法研究”（S202210405016）和东华理工大学研究生创新基金项目“利用地脉动资料研究萍乡矿区地壳浅层三维横波速度结构”（YC2021S644）联合资助。第 58 卷 第 4 期付宇，等：应用闪电连接过程优化算法反演多阶瑞雷波频散曲线Keywords：Rayleigh wave，dispersion curve inversion，multimode，global optimization，lightning attachment proce dure optimization（LAPO）付宇，艾寒冰，姚振岸，等.应用闪电连接过程优化算法反演多阶瑞雷波频散曲线 J.石油地球物理勘探，2023，58（4）：830838.FU Yu，AI Hanbing，YAO Zhen an，et al.Inversion of multimode Rayleigh wave dispersion curves based on lightning attachment procedure optimization J.Oil Geophysical Prospecting，2023，58（4）：830838.0引言瑞雷波是纵波和横波相互干涉叠加形成的一种沿着地表或介质分界面传播的波。起初，瑞雷波只被当作噪声，研究人员主要是根据瑞雷波的特点减小其在地震中带来的危害或压制其在石油地震勘探中对有效信息的影响。随着研究的深入，人们发现瑞雷波在层状介质中传播会出现频散现象，于是便开始将瑞雷波应用于地下地层结构的研究1。与反射地震等常规勘探方法相比，瑞雷波勘探具有无损、施工便捷、浅层分辨率高、检测速度快等优点2，在土层划分3、工程无损检测4、地质灾害勘察5 等领域应用广泛。瑞雷波勘探可分为三个步骤：野外面波数据采集、频散曲线拾取和频散曲线反演。频散曲线的反演作为瑞雷波勘探的关键步骤，直接决定了所求地层信息的可靠性。目前，频散曲线的反演算法主要有两种：局部优化算法和全局优化算法。局部优化算法包括最小二乘法、LevenbergMarquardt（LM）算法、奥卡姆算法等。Dorman等6 首次利用阻尼最小二乘法反演频散曲线，得到的横波速度与折射波法一致，证实了瑞雷波勘探的有效性。Xia等7 将 LM 算法与奇异值分解技术相结合反演瑞雷波相速度，加速了解的收敛，提高了反演的计算效率和稳定性。艾东海等8 在 Constable等9 研究的基础上，将奥卡姆算法用于瑞雷波频散曲线反演，并且还在反演中引入自适应修改拉格朗日乘子光滑系数协调了模型分辨率与解的精度。局部优化算法因计算速度快而被广泛应用，但反演效果过度依赖初始模型10。另外，局部优化算法易陷入局部极值，且涉及偏导数计算，反演结果受雅可比矩阵精度的制约，因此应用受到限制。全局优化算法能避免初始模型选取和偏导数计算，也被应用于频散曲线反演，如遗传算法（GA）、模拟退火算法等。石耀霖等11 将 GA 用于频散曲线反演，通过分析初始搜索结果修改搜索范围，从而提高了搜索效率；Yamanaka等12 在应用 GA 对瑞雷波频散曲线进行反演时，加入精英筛选策略，加快了解的收敛；Dal Moro等13 将GA与后验概率密度估计相结合反演瑞雷波频散曲线，提高了解的精度。Martinez等14 利用模拟退火算法对瑞雷波的群速度和相速度进行反演，降低了对初始模型的依赖，但需经验控制退火过程的一些参数；Beaty等15 采用快速模拟退火算法反演了多阶频散曲线，并与基阶频散曲线反演结果进行了对比，结果表明多阶频散曲线反演的横波速度分辨率更高、不确定性更小、更可靠；Lu等16 基于模拟退火算法提出了热浴模拟退火算法，并与LM算法进行对比，结果表明热浴模拟退火算法更适于频散曲线反演。与局部优化算法相比，全局优化算法有不依赖初始模型、不易陷入局部最优解等优势，但运算时间长、收敛速度慢。为解决全局优化算法反演频散曲线存在的问题，有学者引入了其他全局优化算法，如洗牌蛙跳算法17、自适应混沌遗传粒子群算法18、蚱蜢算法19、自适应蜻蜓算法20 等，也有学者结合多种算法的优势进行反演，如遗传算法与 LM算法的结合21、模拟退火算法与阻尼最小二乘法结合22 等。闪电连接过程优化（Lightning Attachment Procedure Optimization，LAPO）算法模拟雷云向地放电过程，是一种全局优化算法，本文将其用于瑞雷波频散曲线的反演。为详尽评价 LAPO算法反演频散曲线的能力，首先利用不含噪声的模型数据验证了LAPO算法用于频散曲线反演的可行性；其次，往模型数据中加入10%的随机噪声测试了LAPO算法的稳定性；同时应用含噪和不含噪理论模型数据，对比分析了LAPO算法与粒子群（PSO）算法的反演性能；再次，利用理论模型数据测试了LAPO算法的多阶频散曲线反演能力；最后用美国怀俄明地区实际数据检验了LAPO算法实用能力。1LAPO 算法原理1.1原理描述LAPO算法是 Nematollahi等23 根据闪电连接这831石 油 地 球 物 理 勘 探2023 年一物理现象提出的一种自然启发式算法，模拟了雷云对地放电过程。雷云对地放电时，雷云先发展出下行先导，下行先导向下延伸，接近地面目标（高建筑物、高树木等）后，在目标上激发向上发展的上迎先导，当上迎先导与下行先导距离足够小，出现完全击穿，雷击就会产生。下行先导是指从雷云开始发展、向下延伸的先导通道；上迎先导是指从地面（高建筑物、高树木等）开始发展、向上延伸的先导通道24。描述该过程的数学模型可表示为以下五步。（1）初始化。LAPO算法首先需要通过随机方式在定义域中产生一个初始种群。种群的每个个体代表优化问题中的一个解，可以看作云与地面之间的候选点。候选点的产生可以描述为X(0)i=Xlbi+(Xubi-Xlbi)r（1）式中：Xubi、Xlbi分别表示第i个个体的上界和下界；r是0，1 区间的随机数。生成初始种群后根据目标函数计算各种群候选点的电场值或适应度值F(0)i=X(0)i（2）式中为目标函数。（2）确定闪电的下一跳。闪电下一跳的移动需要用到所有候选点的均值 Xave和均值点的适应度值Fave，即 Xave=meanX(t)1，X(t)2，Fave=(Xave)（3）式中t为迭代次数。闪电在传播过程中需在多个点移动，其移动时存在多个选择点。在LAPO算法中对该过程进行数学描述时，对种群中的某一点，考虑了所有其他点都可能成为下一个跳跃点。由于闪电是一个随机行为，对于候选点i，在种群中随机选择一个点j（ij），若点 j的电场值高于均值点的电场值（电场值越大，适应度值越小），则闪电向j点跳跃，即X(t)i，new=X(t)i+r1Xave+r2X(t)j（4）否则，闪电向另一方向移动，即X(t)i，new=X(t)i-r1Xave+r2X(t)j（5）式中：X(t)j为从种群中随机选取的一个点j（ij）；r1、r2都是 0，1 区间的随机数。（3）闪电分支消失。在该阶段，如果新的候选点的电场值高于旧的电场值，则新候选点所在的分支保留，反之，该分支消失。这一过程可以表示为X(t)i=X(t)i，new X(t)i，new=X(t)i （6）以上三步模拟了闪电连接过程的下行先导向地面移动的阶段。（4）上迎先导的移动。在这一阶段，所有的候选点都被认为是上迎先导，并向上移动。上迎先导的移动受到下行先导电荷分布的影响，下行先导的电荷呈指数分布。指数因子定义为S=1-ttmax exp(-ttmax)（7）式中tmax为最大迭代次数。上迎先导的下一候选点的选择过程可以描述为X(t+1)i=X(t)i，new+r S X(t)min-X(t)max（8）式中X(t)min和X(t)max分别为种群第t次迭代的最佳解和最差解，是式（6）计算的所有个体中适应度值的最小值和最大值所对应的点的位置。该步骤完成后若未满足收敛条件（如完成整个迭代过程），则返回步骤（2）。（5）连接点的确定。当上迎先导和下行先导相遇时，连接点确定，闪电过程停止。上迎先导和下行先导相遇这一过程可表达为算法满足收敛条件。1.2算法测试对于全局优化算法性能评估，人们一般采用各种标准函数进行测试。这些测试函数可分为单峰函数、多峰函数和复合函数。本文引用参考文献 23 中Rastrigrin 和 Ackley 两 个 多 峰 函 数 测 试 LAPO 和PSO算法的性能，评价算法探索特性和跳出局部最优能力。Rastrigin函数的表达式为F1(x)=j=1nx2j-10cos(2xj)+10（9）式中 n 为函数维度。测试时设置 n=30，搜索区间为-5.12，5.12 ，最优解位于0。Ackley 函数表达式为F2(x)=-20exp()-0.21nj=1nx2j-exp 1nj=1ncos(2xj)+20+e式中e为自然常数，取值为2.7182。测试时同样设置n=30，搜索区间为-32，32，最优解位于 0。二维情况下，两个函数形状如图1所示。（10）832第 58 卷 第 4 期付宇，等：应用闪电连接过程优化算法反演多阶瑞雷波频散曲线在算法对比时，需控制参数变量以合理对算法进行评估。在寻优时，LAPO算法迭代一次需求取两次目标函数，PSO算法迭代一次只求取一次目标函数，因此前者的种群数设置为 40，是后者种群数 80的一半。迭代次数都设为 200次。PSO 算法的其他参数见参考文献 25。由于两种算法都是基于粒子的随机运动寻优，需要对多次反演结果进行统计分析。参照文献 23，本文对两种方法进行了 30 次试验。30 次测试结果中所求得的最优函数值、平均值及标准差统计如表 1所示（表中数据均保留两位小数），可以看出，对于 Rastrigin和 Ackley两个多峰函数，LAPO算法优于 PSO 算法，说明 LAPO 算法具有更强的全局探索能力。由参考文献 23 中其他测试结果可知，LAPO算法的局部开发能力和收敛性能也强于 PSO算法。2理论模型测试在面波数据中，基阶面波能量强、易观测，在工程中应用最为广泛，因此，本文首先测试了 LAPO算法反演基阶频散曲线的能力。Xia等7 的研究表明，对瑞雷波频散曲线特征变化影响最大的是地层的横波速度和厚度，其余参数影响较小。为节省计算量，本文仅对横波速度和层厚进行反演，其他参数根据先验信息确定。理论模型测试搜索区域的上、下限设定为与真实值相差 50%。反演参数中 LAPO算法的种群数设置为10，PSO算法的种群数设置为20，迭代次数都设置为 50。为避免算法中随机因子对反演结果的影响，每个理论模型测试都进行 20次独立反演，且每次反演时初始模型都随机生成。最后将 20次的均值作为反演结果，将20次反演数据的标准差作为衡量算法稳定性的指标。瑞雷波频散曲线反演的本质就是通过算法搜索到一个最佳横波速度模型，使反演频散曲线与实测频散曲线的拟合误差达到最小。本文使用均方根误差作为目标函数，以评价二者的拟合程度。对于基阶频散曲线反演，目标函数为1=k=1N()vobsk-vcalk2N（11）式中：N为基阶频散曲线的频点数；vobs和vcal分别为观测和反演的瑞雷波相速度。对于多阶频散曲线反演，目标函数为2=1Ll=1Lk=1Nl()vobsl，k-vcall，k2Nl（12）式中：L为频散曲线条数；Nl为第l条频散曲线频点数。频散曲线的正演模拟采用的是快速标量传递算法26。2.1无噪基阶频散曲线反演首先使用模型 A的无噪声理论曲线对 LAPO和PSO算法进行反演测试。模型 A为含软夹层的四层模型，各层的纵波速度（VP）、横波速度（VS）、密度（）和厚度（h）参数如表2所示。当频率大于70 Hz后，该模型观测不到基阶频散曲线，因此频率范围设为 568 Hz，间隔为3 Hz。两种方法反演结果如图 2与表 3所示（表中数据均保留两位小数）。由图 2a可知，LAPO算法反演的图 1二维 Ackley（a）和 Rastrigin（b）函数形状表 1两种方法的标准函数测试结果统计函数名称AckleyRastrigrin理论极值00LAPO算法最优函数值9.50 10-90均值5.87 10-81.53标准差5.05 10-83.70PSO算法最优函数值2.8744.69均值4.6184.67标准差1.0722.44表 2模型 A参数及反演搜索范围层号1234模型参数VSm/s200150250400VPm/s6634978291327g/cm32.02.02.02.0hm235搜索范围VSm/s10030075225125375200600hm1.03.01.54.52.57.5833石 油 地 球 物 理 勘 探2023 年频散曲线几乎和理论模型频散曲线完全重合，PSO算法反演的频散曲线与理论模型频散曲线也拟合得很好。但图 2b 中 LAPO 算法反演的地层横波速度曲线与真实模型基本吻合，PSO 算法反演的横波速度曲线与真实模型却存在一定偏差。另外，LAPO算法反演结果的最大相对误差为 1.38%和最大标准 差 为 6.16，显 著 低 于 PSO 算 法 反 演 结 果 的11.75%和 62.94。因此，与 PSO 算法相比，LAPO算法的寻优能力更强，能搜索到适应度值更小的结果。2.2含噪基阶频散曲线反演实际地震数据不可避免地会包含噪声27，会降低反演结果的准确性，因此有必要检验算法的抗噪能力。在模型A的频散曲线中加入10%的随机噪声28，然后使用LAPO和PSO算法对加噪后的频散曲线进行反演（图3和表4）。本文使用的加噪方式为vnoise=v1+2(0.5-r)p（13）式中：vnoise表示加入噪声后的频散数据；v为加噪前的频散数据；p为加入噪声的程度。由图 3a 可以看出，在频散曲线中加入噪声后，LAPO 算法反演的频散曲线依然能很好地拟合理论频散曲线，PSO算法的拟合效果稍差。在图3b中，与PSO算法相比，LAPO反演的横波速度曲线更接近真实模型。与不含噪声理论频散曲线的反演结果相比，含噪频散曲线反演结果的相对误差和标准差大多有所增大，可见噪声会影响反演结果的准确性。与不含噪理论曲线的反演结果相比，LAPO算法含噪频散曲线反演的最大相对误差从1.38%变为5.44%，最大标准差从6.16变为9.51，增幅较小，整体变化也不大，可见 LAPO算法具有良好的稳定性。PSO算法在反演含噪声理论曲线时，与不含噪声时相比，反演结果的相对误差和标准差整体上有所增加，大于 LAPO 反演含噪声数据和不含噪声数据结果间的变化幅度。由此可知，LAPO 比 PSO 算法的抗噪能力更强、更稳定。从无噪和含噪理论曲线的LAPO和PSO算法各20次反演结果中，选取精度最高的一次进行收敛曲线对比（图4）。对于无噪理论频散曲线，在迭代的前、中期，两种算法收敛精度和速度相差不大；在中、后期PSO算法早熟，LAPO算法还在继续收敛；迭代结束后 LAPO算法的收敛精度更高（图 4a）。对于含噪理论频散曲线，多数情况下，LAPO算法的收敛速度和收敛精度高于优于 PSO算法。可见，与 PSO算法相比，LAPO算法具有更高的收敛速度和精度。表 3模型 A无噪声理论数据两种算法反演结果统计参数VS1/(ms-1)VS2/(ms-1)VS3/(ms-1)VS4/(ms-1)h1/mh2/mh3/m真实值200150250400235LAPO算法反演均值200.04150.13249.84399.911.983.045.00相对误差%0.020.090.070.020.951.380.08标准差3.472.453.716.160.130.330.27PSO算法反演均值197.54153.60263.27400.272.233.124.53相对误差%1.232.405.310.0711.753.929.41标准差28.1921.1862.947.100.901.211.74图 3模型 A含噪理论频散曲线两种方法反演结果对比（a）频散曲线；（b）横波速度曲线图 2模型 A无噪声理论数据两种算法反演结果对比（a）频散曲线；（b）横波速度曲线834第 58 卷 第 4 期付宇，等：应用闪电连接过程优化算法反演多阶瑞雷波频散曲线 图 4 模型 A无噪（a）和含噪（b）理论曲线的 两种方法反演迭代收敛曲线对比2.3无噪多阶频散曲线反演面波在一些特殊地层中传播时，可能出现高阶模式占主导的情况29。若联合反演基阶和高阶频散曲线，能增加对反演结果的约束，提升反演结果的精度。模型 B（表 5）为含硬夹层的四层模型，理论频散曲线的频率范围为598 Hz，频率间隔为3 Hz。模型 B的基阶和一阶频散曲线联合 LAPO算法反演结果如图5和表6所示。在图5a中，反演的基阶和一阶频散曲线与理论频散曲线均拟合较好。由图5b可以看出，反演的横波速度曲线与真实模型曲线非常接近。对比基阶模式频散曲线的反演结果与多阶频散曲线的联合反演结果（表6）可以明显看出，多阶频散曲线联合反演的模型参数相对误差和标准差整体上都有所降低。可见，利用多阶频散曲线联合反演的结果确实更准确、更可靠，也说明 LAPO算法不仅可用于基阶频散曲线的反演，也适用于多阶频散曲线的反演。表 5模型 B参数及反演搜索范围层号1234模型参数VSm/s160260200400VPm/s 531 862 6631327g/cm32.02.02.02.0hm235搜索范围VSm/s80240130390100300200600hm1.03.01.54.52.57.5图 5模型 B无噪声多阶理论频散曲线反演结果（a）频散曲线；（b）横波速度曲线2.4含噪多阶频散曲线反演向模型 B的多阶频散曲线中加入 10%的随机噪声，对 LAPO算法进行测试，结果如图 6和表 7所示。表 4模型 A含噪理论曲线两种算法反演结果统计参数VS1/(ms1)VS2/(ms1)VS3/(ms1)VS4/(ms1)h1/mh2/mh3/m真实值200150250400235LAPO反演均值198.91149.04249.25396.942.112.945.01相对误差/%0.540.640.300.765.441.930.19标准差4.773.536.449.510.230.290.45PSO反演均值207.10154.68278.46397.792.173.224.59相对误差/%3.553.1211.380.558.367.438.12标准差36.8925.5281.3211.210.781.282.17表 6模型 B基阶与多阶无噪声频散曲线 LAPO 算法反演结果统计对比参数VS1/(ms1)VS2/(ms1)VS3/(ms1)VS4/(ms1)h1/(m)h 2/(m)h 3/(m)真实值160260200400235基阶数据反演均值160.60260.82200.68401.352.032.985.05相对误差/%0.370.310.340.341.520.730.98标准差1.773.261.902.140.070.180.17多阶数据反演均值159.74259.37199.76399.362.002.945.01相对误差/%0.170.240.120.160.151.970.24标准差1.592.512.104.000.060.200.21835石 油 地 球 物 理 勘 探2023 年由图 6可知，无论是反演的频散曲线还是横波速度曲线都与真实曲线拟合得很好。与表 6相比，表 7中统计数据变化较小，表明LAPO算法用于含噪多阶频散曲线反演表现不错。表 7模型 B含噪多阶理论频散曲线反演结果统计参数VS1/(ms1)VS2/(ms1)VS3/(ms1)VS4/(ms1)h1/mh2/mh3/m真实值160260200400235反演均值160.23260.84200.66400.782.002.985.08相对误差/%0.140.320.330.190.010.521.55标准差1.933.152.214.490.100.170.213实际资料测试实际瑞雷面波数据采集自美国怀俄明地区30。由锤击震源激发，48个8 Hz 垂直检波器接收，最小炮检距和道间距都为 0.9 m。原始地震记录如图 7a所示。图7b为频散能量谱，从中提取基阶和二阶频散曲线用于反演。实际数据反演参数与理论模型测试保持一致。根据测井资料岩性解释将模型划分为 5层，估计的横波速度等其他物性参数如表8所示。表 8怀俄明地区面波反演搜索范围及模型参数设置31层号12345VSm/s100300100400100600200600200800hm15151515泊松比0.380.380.350.350.30g/cm32.02.02.02.02.0图8a为实际数据LAPO算法的基阶和多阶反演的频散曲线与实测曲线对比，三者很吻合。图 8b为LAPO 算法反演过程中目标函数随迭代次数的变化曲线，可以看出，使用多阶数据比只使用基阶数据的反演精度更高。多阶数据反演的横波速度模型 （图 8c 中红色实线）比基阶数据反演得到的横波速度模型（图 8c中黑色实线）更接近测井数据（图 8c中带菱形的灰色折线），尤其是在浅层，说明 LAPO 算法对于多阶瑞雷面波频散曲线反演具有良好的实用性。图 7美国怀俄明地区实际面波地震记录（a）和频散能量谱（b）31图 6模型 B含噪多阶理论频散曲线反演结果（a）频散曲线；（b）横波速度曲线836第 58 卷 第 4 期付宇，等：应用闪电连接过程优化算法反演多阶瑞雷波频散曲线4结束语本文将一种新颖的全局优化算法LAPO 算法用于瑞雷波频散曲线反演。在反演时，通过设置较大的模型参数搜索范围以模拟无先验信息或者仅有少量已知信息的情况。利用两个理论模型对 LAPO算法的性能进行了测试，最后再应用美国怀俄明地区的实测数据检验了LAPO的实用性。测试结果表明：LAPO算法能有效用于频散曲线反演，该算法在获得高精度地层模型参数的同时，还保证了反演过程收敛速度以及较强的稳定性，具有良好的发展潜力。对于瑞雷波频散曲线反演，LAPO算法比传统的粒子群算法，精度更高、收敛更快、稳定性更强。参 考 文 献1 邱新明,王赟,韦永祥,等.多分量面波研究进展J.石油物探,2021,60(1):4656.QIU Xinming,WANG Yun,WEI Yongxiang,et al.Advancements in multicomponent surface waves:A reviewJ.Geophysical Prospecting for Petroleum,2021,60(1):4656.2 伍敦仕,孙成禹,林美言.基于互相关相移的主动源地震面波频散成像方法J.地球物理学进展,2017,32(4):16931700.WU Dunshi,SUN Chengyu,LIN Meiyan.Active seismic surface wave dispersion imaging method based on crosscorrelation and phaseshiftingJ.Progress in Geophysics,2017,32(4):16931700.3 宫胜家.瞬态面波的勘探及应用J.中国煤炭地质,2008,20(增刊 1):8486.GONG Shengjia.Prospecting and application of transient ground rollJ.Coal Geology of China,2008,20(S1):8486.4 赵硕九,金钢燮,汪利民.瑞雷面波在多尺度勘探上的 研 究 现 状 与 展 望 J.勘 探 地 球 物 理 进 展,2010,33(1):110.ZHAO Shuojiu,JIN Gangxie,WANG Limin.Review of Rayleigh wave techniques in survey of different scales and its future prospectsJ.Progress in Exploration Geophysics,2010,33(1):110.5 许新刚,岳建华,李娟娟,等.面波勘查技术及在滑坡地质调查中的应用研究J.地球物理学进展,2016,31(3):13671372.XU Xingang,YUE Jianhua,LI Juanjuan,et al.Surface wave exploration technique and its application in landslide geology surveyJ.Progress in Geophysics,2016,31(3):13671372.6 DORMAN J,EWING M.Numerical inversion of seismic surface wave dispersion data and crustmantle structure in the New YorkPennsylvania areaJ.Journal of Geophysical Research,1962,67(13):52275241.7 XIA J,MILLER R D,PARK C B.Estimation of nearsurface shearwave velocity by inversion of Rayleigh wavesJ.Geophysics,1999,64(3):691700.8 艾东海,程庆群.低速软弱夹层二维横波速度结构的OCCAM 反演J.工程勘察,2009,37(4):8790.AI Donghai,CHENG Qingqun.Estimation of 2D Swave velocity section with low velocity layers by OCCAM algorithmJ.Geotechnical Investigation&Surveying,2009,37(4):8790.9 CONSTABLE S C,PARKER R L,CONSTABLE C G.Occam s inversion:A practical algorithm for gene rating smooth models from electromagnetic sounding dataJ.Geophysics,1987,52(3):289300.10 王一鸣,宋先海,张学强.瑞雷面波频散曲线的粒子群蚁群混合优化反演J.石油地球物理勘探,2022,57(2):303310,356.WANG Yiming,SONG Xiahai,ZHANG Xueqiang.Inversion of Rayleigh wave dispersion curves based on particle swarm and ant colony hybrid optimizationJ.Oil Geophysical Prospecting,2022,57(2):303310,356.11 石耀霖,金文.面波频散反演地球内部构造的遗传算法J.地球物理学报,1995,38(2):189198.SHI Yaolin,JIN Wen.Genetic algorithms inversion of lithospheric structure from surface wave dispersionJ.Chinese Journal of Geophysics,1995,38(2):189198.12 YAMANAKA H,ISHIDA H.Application of genetic algorithms to an inversion of surfacewave dispersion dataJ.Bulletin of the Seismological Society of America,1996,86(2):436444.13 DAL MORO G,PIPAN M,GABRIELLI P.Rayleigh wave dispersion curve inversion via genetic algorithms and marginal posterior probability density estimationJ.Journal of Applied Geophysics,2007,61(1):3955.14 MARTNEZ M D,LANA X,OLARTE J,et al.Inversion of Rayleigh wave phase and group velocities by simulated annealingJ.Physics of the Earth and Planetary Interiors,2000,122(12):317.15 BEATY K S,SCHMITT D R,SACCHI M.Simulated annealing inversion of multimode Rayleigh wave dispersion curves for geological structureJ.Geophysi图 8怀俄明地区基阶和多阶瑞雷波相速度反演结果对比（a）频散曲线；（b）收敛曲线；（c）横波速度曲线837石 油 地 球 物 理 勘 探2023 年cal Journal Int