基于Lebenberg_Marquardt的声发射源定位算法.pdf

Advances in Applied Mathematics 应用数学进展应用数学进展,2024,13(4),1327-1333 Published Online April 2024 in Hans.https:/www.hanspub.org/journal/aam https:/doi.org/10.12677/aam.2024.134122 文章引用文章引用:王子维,阎少宏.基于 Lebenberg_Marquardt 的声发射源定位算法J.应用数学进展,2024,13(4):1327-1333.DOI:10.12677/aam.2024.134122 基于基于Lebenberg_Marquardt的声发射源定位的声发射源定位 算法算法 王子维，阎少宏王子维，阎少宏*华北理工大学理学院，河北 唐山 收稿日期：2024年3月17日；录用日期：2024年4月11日；发布日期：2024年4月18日 摘摘 要要 声发射是一种寻声发射是一种寻找材料破裂点或地震发声源位置非常有效的方法，找材料破裂点或地震发声源位置非常有效的方法，在本次实验中介绍了声发射源定位在本次实验中介绍了声发射源定位中时差定位法的基本原理中时差定位法的基本原理，将定位问题转化为多元非线性方程组求解问题。根据问题特点，将传统的，将定位问题转化为多元非线性方程组求解问题。根据问题特点，将传统的Geiger定位算法与定位算法与Lebenberg_Marquardt算法相结合，通过改进算法相结合，通过改进Geiger算法的迭代方向与迭代步，极算法的迭代方向与迭代步，极大大地地减少了算法的减少了算法的时间复杂度，提高了时间复杂度，提高了算法的精确度。数值仿真算例表明：该方法能有效地算法的精确度。数值仿真算例表明：该方法能有效地解决了迭代解决了迭代时雅可比矩阵与其转置矩阵接近奇异时时雅可比矩阵与其转置矩阵接近奇异时dk过大的问题，保证了算法的收敛性并且提高过大的问题，保证了算法的收敛性并且提高了了迭代算法的收敛迭代算法的收敛速度。速度。关键词关键词 声发射定位，声发射定位，Lebenberg_Marquardt算法，算法，Geiger算法，迭代方向算法，迭代方向 Sound Source Localization Algorithm Based on Levenberg_Marquardt Ziwei Wang,Shaohong Yan*College of Science,North China University of Science and Technology,Tangshan Hebei Received:Mar.17th,2024;accepted:Apr.11th,2024;published:Apr.18th,2024 Abstract Acoustic emission is an extremely effective method for locating material fracture points or seismic sound sources.In this experiment,the basic principle of time difference localization in acoustic emission source localization is introduced,which transforms the localization problem into a mul-*通讯作者。王子维，阎少宏 DOI:10.12677/aam.2024.134122 1328 应用数学进展 ti-dimensional nonlinear equation solving problem.Based on the characteristics of the problem,the traditional Geiger localization algorithm is combined with the Levenberg_Marquardt algorithm by improving the iteration direction and step of the Geiger algorithm,greatly reducing the time complexity of the algorithm and improving its accuracy.Numerical simulation examples demon-strate that this method effectively solves the problem of excessive dk when the Jacobian matrix and its transpose matrix are close to singular during iteration,ensuring the convergence of the algo-rithm and improving the convergence speed of the iterative algorithm.Keywords Acoustic Emission Localization,Levenberg_Marquardt Algorithm,Geiger Algorithm,Iteration Direction Copyright 2024 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 声发射(acoustic emission,AE)1 2 3 4作为一种研究材料破裂重要手段，在岩石力学应用领域应用十分广泛。如何精确地找到材料破裂点或地震发声源的位置，是声发射中的一个重要问题。在传统的声发射问题中，比较常用的算法有：最小二乘法、Bayesian 定位方法、相对定位法、Geiger 定位方法和单纯形定位方法等。目前，国内外学者针对声发射定位方法开展了大量的研究工作，比较常用的有最小二乘法5，Geiger法6，单纯形法7以及牛顿法8等。其中，Geiger 法是一个经典的定位技术，其中使用到地震到时(即地震波到达不同地震站的时间)进行震源位置的估计。林峰9等研究人员提出了一种联合定位方法，将线性定位方法和 Geiger 定位方法相结合，有效解决了线性定位精度不高的问题，并对常规 Geiger 定位方法的初值选择进行了优化。同时，Li Nan 10等人提出了一种基于单纯形法的非线性微震震源定位方法，显著提高了震源定位的精度。此外，李会义11等人利用牛顿线性迭代法成功求解定位方程组，实现了微震定位的可靠监测。传统的 Geiger 方法可能会使用迭代重新定权最小二乘方法进行求解。该过程中，可能会对方程或模型的某些部分进行简化或假设以便计算。在这种迭代类中常会遇到迭代步中雅可比矩阵与其转置矩阵的乘积接近于奇异的情况，在这种情况下会直接导致迭代步过大，从而使误差过大。针对上述问题，本文引进了 Levenberg_Marquardt 算法12。作为一种结合了梯度下降和高斯牛顿方法的优化技术，可以用来提高 Geiger 定位方法的效度。利用 LM 算法，可以更准确地进行震源位置的求解，因为它在数据不是完全线性的情况下仍然有效，并且提供了一种智能化的方式来调整步长以平衡收敛速度和稳定性。不仅提高了算法本身的鲁棒性，而且很大程度上降低了算法的时间复杂度。2.声发射源三维时差定位方法声发射源三维时差定位方法 2.1.Geiger 算法算法 Geiger 算法是声发射源定位中常用的一种定位算法。该算法基于声发射信号在不同传感器之间的传播时间差(time difference of arrival，TDOA)来确定声源的位置。Geiger 算法的基本思想是通过将声源位置的坐标表示为未知变量，并将 TDOA 表示为已知常数，将定位问题转化为一个非线性方程组求解问题。Open AccessOpen Access王子维，阎少宏 DOI:10.12677/aam.2024.134122 1329 应用数学进展 该方程组的解即为声源的位置。算法的核心是迭代过程，其中使用的最小二乘法用于拟合实测和计算得到的 TDOA 值，以更新声源位置。在每次迭代中，根据当前的声源位置估计值和 TDOA 信息，计算残差，并构建目标函数。然后，使用数值优化算法(如 Levenberg_Marquardt 算法)求解目标函数，并得到更新的声源位置估计值。重复这一过程，直到迭代收敛或达到预定的停止条件。Geiger 算法在声发射源定位中具有较高的定位精度和稳定性，并且在实践中得到广泛应用。通过结合其他改进算法，如 Levenberg_Marquardt 算法，可以进一步提高算法的性能和收敛速度。Geiger 定位算法是从一个给定的初始点(),x y z t v通过迭代而接近最终结果。每一次迭代都基于最小二乘法计算一个修正向量(),xyztv ，把向量加到上次迭代的结果(试验点)上，得到一个新的试验点，然后判断这个新试验点是否满足要求；如果满足要求，此点坐标即所求震源位置，如果不满足，则继续迭代。根据多个传感器获得的到达时间所建立的式(1)所给出的固定方程的到声发射源位置坐标。()()()()222iiiixxyyzzv tt+=(1)式中：,iiix y z代表第 i 个传感器的坐标位置；,x y z代表声发射源坐标；v 代表声波波速；it代表 P 波到达第 i 个传感器的时间，t 代表发声源发声的时间。在(1)式中未知变量共有 5 个，理论上我们只要建立 5 个方程就可以求出声发射源位置坐标。或通过建立 5 个以上的方程得到超定方程组对声发射源进行定位。现在假设有 n(n 5)个传感器接收到有效信 号：()()()()222iiiixxyyzzv tt+=(1,2,3,in=)。即可建立需求解的非线性方程组。我们令：()()()()222iiiiifxxyyzzv tt=+。则我们需求解的非线性方程组为：12000nfff=其中()1,2,ifn为给定 n 维欧式空间nR中的实值函数。引进向量记号，令()()()()12nflflF lfl=xylztv =则上式可写成：()0F l=。若存在*lR，使得()*0F l=，则称*l是方程的解。考虑非线性方程组：()0F l=其中():nmF lRR是连续可微的及()F l连续可微。设()J l是()F l的 Jacobian 矩阵 王子维，阎少宏 DOI:10.12677/aam.2024.134122 1330 应用数学进展 ()()()()()()()()()()()()()()()()111111111111222222222222nnnnnnnnnnnnflflflflflxxyyzzvttxyztvRRRxxyyzzflflflflflvttRRRJ lxyztvxxyyzzflflflflflvttRRRxyztv=利用泰勒公式，我们有：()()()()()()20000F lF lJ llloll=+(2)一般0l充分接近*l，略去上式的高阶小量，因而可用线性方程组()()()()0000F lF lJ lll=+=(3)代替(1-1)，设(1-3)的解为1l，则()()11000llJ lF l=(4)()()11kkkkllJ lF l+=(5)计算迭代步为：1kkkdJF=(6)其中()kkFF l=，()kkJFx=。我们将问题(1-1)的求解转化为求解()421min2x RF l (7)可知 Geiger 算法是牛顿法在非线性最小二乘法问题上的应用，其搜索方向为：()1TTkkkkkdJ JJ F=(8)故问题(1-1)的求解转化为：()1TT1kkkkkkllJ JJ F l+=(9)2.2.Levenberg_Marquardt 算法算法 LM 算法是一种非线性优化方法，处于牛顿法和梯度下降法之间，对于过参数化问题不敏感，能有效处理冗余参数问题，从而降低代价函数陷入局部极小值的可能性。因此，LM 算法在计算机视觉等领域有着广泛的应用。在上述的最小二乘问题中，由于我们所省略的二阶项对 Geiger 算法在收敛存在上有一定的影响。为了得到更加有效的方法，在运用 Geiger 算法时加上搜索策略即：1kkkkxxd+=+(10)其中k为搜索因子，k利用 Wolfe-Powell 准则确定()()()()()()()22T22T1kkkkkkkkkkkkF ldF lFl dF ldF lFl d+(11)王子维，阎少宏 DOI:10.12677/aam.2024.134122 1331 应用数学进展 其中10,2。Geiger 算法对于绝大多数的三维声发射问题中的最小二乘是全局收敛的，但对于某些问题仍然可能收敛很慢。且进行迭代求解时当迭代步中TkkJ J接近于奇异时kd过大，而且当TkkJ J是奇异时，其迭代步是没有意义的。因此使用 LM 算法对 Geiger 算法进行修正，对于当前迭代点 LM 算法搜索方向为:()()1TTkkkkkkkdJ JIJ F=+(12)其中k为正参数，LM 是对 G-N 迭代步的修正。k是为了防止TkkJ J接近奇异时kd过大。而且当TkkJ J是奇异时。迭代步是没有定义的。LM 克服了迭代对矩阵()J l必须满秩的要求。在得到修正向量kd后，以1kkkkxxd+=+为新的实验点继续迭代，直到满足误差要求。3.基于基于 Levenberg_Marquardt 的的声发射源定位算法声发射源定位算法 3.1.数值仿真数值仿真 数据来源为康玉梅等所做的实验仿真。为了定位声发射源在试件中的位置，需要根据虚拟声发射源位置、传感器坐标和材料波速生成一组模拟到时数据，并通过处理这组数据实现声发射源的定位。为了考察定位结果的重复性，在声发射源位置生成了三组带有随机误差的到时数据。由于一定量的到时误差与传感器布局有关，因此会影响定位结果的准确性。假设试件尺寸为 100 mm 100 mm 200 mm，模拟试样的传感器布置和模拟声发射源的位置如图1 所示。以点 A 为原点，依此可确定传感器位置坐标，详见表 1。试验中的破裂点为点 M(30,40,50)和点 C(100,100,0)，破裂时间为 100 秒。Figure 1.Schematic diagram of the location of sensors and acoustic emission sources 图图 1.传感器和声发射源的位置示意图 Table 1.Coordinates of the locations of various sensors 表表 1.各个传感器的位置坐标 传感器 1 2 3 4 5 6 7 8 坐标(50,0,20)(100,50,20)(50,100,20)(0,50,20)(50,0,180)(100,50,180)(50,100,180)(0,50,180)表 2 为各个传感器的真实到时以及含随机误差的模拟到时数据。根据不同的时差定位方法，可以计王子维，阎少宏 DOI:10.12677/aam.2024.134122 1332 应用数学进展 算出不同的声发射源位置。Table 2.Actual arrival times of each sensor and data on arrival times with three sets of random errors 表表 2.各个传感器的真实到时和三组随机误差到时数据 定位点 传感器 模拟真实到时 含随机误差的模拟到时数据 第 1 组 第 2 组 第 3 组 M 1 100.0000108 100.0000107 100.0000111 100.0000108 2 100.0000154 100.0000154 100.0000163 100.0000157 3 100.0000140 100.0000144 100.0000144 100.0000149 4 100.0000087 100.0000080 100.0000094 100.0000092 5 100.0000275 100.0000270 100.0000282 100.0000265 6 100.0000296 100.0000296 100.0000305 100.0000296 7 100.0000290 100.0000295 100.0000288 100.0000289 8 100.0000268 100.0000274 100.0000268 100.0000268 C 1 100.0000227 100.0000231 100.0000221 100.0000234 2 100.0000108 100.0000101 100.0000106 100.0000111 3 100.0000108 100.0000104 100.0000099 100.0000098 4 100.0000227 100.0000236 100.0000225 100.0000233 5 100.0000424 100.0000424 100.0000426 100.0000423 6 100.0000374 100.0000382 100.0000383 100.0000367 7 100.0000374 100.0000374 100.0000374 100.0000374 8 100.0000424 100.0000428 100.0000424 100.0000425 3.2.算法的实现算法的实现 针对有 8 个传感器的实验环境中，分别使用最小二乘法、Geiger 法以及修正算法进行计算实验。表3 为计算出的各个定位点的位置、计算误差以及迭代次数。从表 3 中还可以看出，Levenberg_Marquardt算法对 Geiger 算法进行修正后明显减少了迭代步数，使收敛速度更快，这对于声发射实时监测的海量数据的处理是非常有意义的。Table 3.Actual arrival times of each sensor and data on arrival times with three sets of random errors 表表 3.各个传感器的真实到时和三组随机误差到时数据 定位方法 实际破裂点/mm 定位点 计算误差 迭代次数 最小二乘法 1(30,40,50)(34.81,37.24,57.17)2.780885*106-2(27.46,43.78,52.06)5.416332*107-3(29.39,43.01,47.79)1.078191*106-1(100,100,0)(95.35,94.03,0.98)1.017506*106-2(92.95,97.40,0.28)5.533379*107-3(101.90,99.59,0.66)6.573619*107-王子维，阎少宏 DOI:10.12677/aam.2024.134122 1333 应用数学进展 续表 Geiger 法 1(30,40,50)(28.64,40.50,50.66)3.115651*107 5 2(28.42,38.04,48.71)3.470470*107 5 3(31.74,44.01,51.04)1.509080*107 5 1(100,100,0)(104.78,102.91,0.50)4.293240*107 18 2(97.75,101.10,0.50)2.505926*107 8 3(104.60,109.18,0.35)3.196705*107 8 Levenberg_Marquardt算法 1(30,40,50)(31.45,37.45,49.23)3.032803*107 4 2(28.17,42.79,50.37)3.354007*107 4 3(30.10,38.86,51.08)1.573923*107 4 1(100,100,0)(104.72,102.81,0.51)4.293178*107 12 2(97.74,101.09,0.49)2.505863*107 6 3(104.61,109.17,0.35)3.196455*107 5 4.结语结语 本文主要介绍了三维声发射定位的原理，以及以 Gergir 算法为基础，使用 Levenberg_Marquardt 算法对其迭代步进行了修正，一方面不仅仅在精度上相较于普通的定位算法得到了提高，而且相较于传统最小二乘法的迭代速度也得到了提高，另一方面解决了迭代中容易出现的迭代步过大的问题。很好地降低了算法的时间复杂度，提高了算法本身的鲁棒性。参考文献参考文献 1 张宇,刘嘉成,冯舒,等.针对不连续带洞结构的声发射源定位方法研究J.仪器仪表学报,2023,44(11):282-289.https:/doi.org/10.19650/ki.cjsi.J2311591 2 王巍浩,杜西岗,许城,等.三轴加载下灰岩的力学特性及声发射特征J.矿业研究与开发,2024,44(2):129-137.3 胡杰,张延欢,刘洪亮,等.隧道围岩块体垮塌声发射及自振频率演化特征试验研究J/OL.中国公路学报:1-18.http:/ 黄冬梅,卢欢欢,潘鑫,等.单轴压缩条件下含孔洞裂隙红砂岩力学行为及声发射特征分析J/OL.矿业安全与环保:1-10.http:/ 何永红,靳鹏伟.小波稳健最小二乘估计的轨道误差校正算法及其在DEM中的应用J.大地测量与地球动力学,2020(12):40.6 刘金丽,高见.联合采用Inglada算法和Geiger算法确定地震精确位置J.防灾科技学院学报,2012,14(2):54-59.7 刘建坡,王洪勇,杨宇江,等.不同岩石声发射定位算法及其实验研究J.东北大学学报(自然科学版),2009,30(8):1193-1196.8 姜天琪,裴烁瑾.基于网格搜索-牛顿迭代法的微震震源定位算法J.矿业科学学报,2019,4(6):480-488.https:/doi.org/10.19606/ki.jmst.2019.06.002 9 林峰,李庶林,薛云亮,等.基于不同初值的微震源定位方法J.岩石力学与工程学报,2010,29(5):996-1002.10 李楠,王恩元,孙珍玉,等.基于 L1 范数统计的单纯形微震震源定位方法J.煤炭学报,2014,39(12):2431-2438.https:/doi.org/10.13225/ki.jccs.2013.1855 11 李会义,姜福兴,杨淑华.基于 Matlab 的岩层微地震破裂定位求解及其应用J.煤炭学报,2006(2):154-158.12 徐克彬,陈祖斌,刘玉海,等.基于 L-M 算法的微地震定位方法J.石油地球物理勘探,2018,53(4):765-769+790+653.https:/doi.org/10.13810/ki.issn.1000-7210.2018.04.013