工程6（2020）827研究应用地球物理学-文章三维含孔结构的无速度MS/AE源定位方法董龙军a，胡庆春a，童晓杰a，刘有芳ba中南大学资源与安全工程学院，湖南长沙410083b美国科罗拉多矿业学院地球物理系，科罗拉多州戈尔登，邮编80401阿提奇莱因福奥文章历史记录：收到2019年2019年11月15日修订2019年12月16日接受2020年1月11日在线关键词：微震源声发射三维含孔结构的无速度定位法A B S T R A C T微震/声发射（MS/AE）源定位方法是复杂结构潜在危险源预测与控制的关键。然而，现有的定位方法对不规则结构和预先测量的速度引起的定位误差知之甚少针对三维复杂含孔结构的高精度定位要求，提出了一种无速度的MS/AE源定位方法。该算法采用等距网格点搜索路径，避免了人工重复训练，引入了A* 搜索算法，并利用网格点适应不规则孔洞的复杂结构。它还利用了无速度源定位方法。为了验证该方法的有效性，对尺寸为10cm× 10 cm× 10 cm的立方体混凝土试件进行了它被剪掉了放入一个/6 cm× 10 cm的圆柱形空间中根据到达时间，经典的盖革方法和所提出的方法被用来定位断铅源。结果表明，该方法的定位误差为1.20 cm，小于Geiger法的2.02 cm。因此，该方法可以有效地定位源在复杂的三维结构与孔，并达到更高的精度要求。©2020 THE COUNTORS.Elsevier LTD代表中国工程院出版，高等教育出版社有限公司。这是一篇CC BY-NC-ND许可下的开放获取文章（http://creativecommons.org/licenses/by-nc-nd/4.0/）。1. 介绍随着浅层矿产资源的逐渐枯竭和地下空间的开发利用，地下工程逐渐向地球深部发展由于深部巨大的地应力和构造应力，大量的能量在深部岩体中积累在地下工程开挖过程中，岩体中储存的能量在卸荷过程中瞬时释放，往往导致岩爆灾害。这种高应力引起的灾害往往会对岩体工程活动造成极大的破坏和人员伤亡。岩爆灾害涉及岩体力学演化过程中的诸多复杂因素，难以用传统的岩石力学理论进行分析，是当今世界性难题之一。微地震/声发射震源监测技术是监测和预报冲击地压的有效手段[1]。目前已广泛应用于煤矿安全监测中，*通讯作者。电子邮件地址：lj. csu.edu.cn（L. Dong）。许多深井和高应力矿井的地压[2此外，加拿大、南非等国也建立了多个国家级矿山MS/AE监测网络系统[7，8]。MS/AE监测技术也得到了迅速发展，并在各个领域得到应用，如航空、军工、桥梁结构、隧道工程、精密制造等[9MS/AE源定位是MS/AE监测中最经典、最基本的问题之一[13]。1912年，Geiger[14]首先提出了一种定位方法。因此，提高MS/AE源定位的精度和效率成为MS/AE监测技术领域的重要研究内容。随着大量定位方法的引入，震源定位的误差和定位事件所需的时间也在减少。经典的震源定位方法有：Ingladahttps://doi.org/10.1016/j.eng.2019.12.0162095-8099/©2020 THE COMEORS.由爱思唯尔有限公司代表中国工程院和高等教育出版社有限公司出版。这是一篇基于CC BY-NC-ND许可证的开放获取文章（http://creativecommons.org/licenses/by-nc-nd/4.0/）。可在ScienceDirect上获得目录列表工程杂志主页：www.elsevier.com/locate/eng828L. Dong等人/工程6（2020）827该方法消除了测量区域与实际MS/AE区域速度差异的影响，大大提高了定位精度。节省了因预先测速而造成的人员、时间和经济成本，比传统方法更方便实用在他们的初步研究之后，Dong等人开发了几种方法并取得了良好的结果，包括三维综合分析定位[24，25]，多步定位[26]以及数值和分析协同定位[27]。随着研究的深入，一些学者开始关注物体的几何不规则性对定位的影响。通过使用训练数据，Park等人。[9]，Baxter等人。[28]和Eaton等人。[29]建立了一个数据集来评估源位置。从理论上讲，这类方法可以应用于物体表面的任何源，但对块结构内部的源的适应性较差目前的此类研究的实验对象都是平面结构，而且这些定位方法非常耗时。虽然Gollob等人[30]考虑了几何形状对定位的影响，但他们没有讨论由所用速度与实际波速度为了提高三维含孔结构中声源的定位精度，提出了一种基于A* 搜索算法的三维含孔结构无速度MS/AE声源定位方法。称为VFH（velocity-free for hole-containing structure）。该算法采用等距网格点搜索路径，避免了人工的重复训练，引入了A* 搜索算法，利用网格点来适应具有不规则孔洞的复杂结构2. 方法VFH采用改进的A* 搜索算法和无速度测量的定位方法对MS/AE源进行定位该方法分为4个部分：首先对需要定位的目标进行网格划分，用0和1表示目标的形状;其次获取传感器接收到的MS/AE事件产生的P波到达信号;然后采用A* 搜索算法搜索各传感器到各网格点的最快波形路径;最后引入最小偏离量D计算MS/AE源的 图 1示出了VFH的定位过程的流程图。本节介绍该方法的各个部分2.1. 确定初始环境在定位区域中确定空区域的几何形状和具体位置。并根据空区情况和位置精度要求总体上图1.一、VFH定位过程的流程图Sk：接收信号的第k个传感器;Dxyz：点Pxyz与未知MS/AE源P0的偏离程度;Dtlm：两个不同传感器Sl和Sm的实际时间差;Dtlm：两个不同站l和m之间的行程时间差。0xyzðÞxyzxyzxyzð Þð Þð ÞL. Dong等/工程6（2020）827-834829网格越密，定位精度越高，计算量也会成倍增加。然而，当网格划分足够密集时，分区并不能明显提高定位精度建立与网格节点大小相同的零矩阵M，并将网格节点对应的矩阵索引位置x;y;z标记为1。网格节点形成一个集合，该集合用作在后续节点中搜索最快波形路径的起点。假设P波在周围非空区域的传播速度是一个固定的未知值，并表示为C. 对于未知源P0，将其位置坐标设置为x0;y0;z0，以及激励到t0的初始时间。2.2. 搜索最快的波形路径以集合中的每一个网格点Pxyz作为潜在源激发位置，追踪最短路径，得到理论解。最短路径Lk从采样点到第k个传感器。如果图二、传统的A* 搜索算法（a）当前节点连接到采样点P xyz在空白区域，L k 被认为（b）图2（a）的主视图;（c）通过传统A* 搜索算法搜索的路径。为1. VFH使用改进的A* 搜索算法跟踪最短路径Lk。2.2.1. A* 搜索算法A* 搜索算法是由斯坦福研究所（现在的SRI International）的Peter Hart，Nils Nilsson和Bertram Raphael在1968年提出的[31，32]。在计算机科学中，A* 搜索算法作为Dijkstra算法的一种扩展，由于其高效性而被广泛用于路径查找和图遍历A* 搜索算法通过使用启发式算法来指导搜索，从而获得了更好的性能。在其主循环的每次迭代中，A* 搜索算法基于函数fm选择节点，由两个函数dm和hm组成。A * 搜索算法在每一步选择具有最小值fmfm d mhm 1其中n表示所选路径的当前节点;d m用于计算从起始节点到节点m的确切最快路径的距离;h m是一个启发式函数，可用于估计从节点n到目标的最小距离[33]。在本文中，两个节点之间的欧几里德距离作为估计的最小距离。A* 搜索算法在从起始点扩展到目标点或没有路径可以扩展时停止搜索。在实际的三维应用中，传统的A* 搜索算法采用中心点搜索路径，一般只考虑相邻层的26个节点来选择下一个节点，如图1和图2所示。 2（a）和（b）。在"L“型区域中，通过使用传统的A* 搜索算法来跟踪最短路径，并且所从图中可以发现，在追踪最短路径时存在两个不合理的地方：①得到的路径具有尖锐的锯齿状。这是由于传统A* 搜索算法的局限性造成的②路径的节点是网格的中心，这意味着传感器也要放置在网格的中心，这在实际应用中是不适用的。2.2.2. 改进的A* 搜索算法为了更有效地跟踪最短路径，对A* 搜索算法进行了改进，使其按网格点进行搜索，如图1所示。3.第三章。这允许搜索路径的节点位于网格节点上，这也意味着传感器可以附着到对象的表面节点。更符合实际情况。图三.当前栅格节点和关联栅格节点之间的潜在方向。(a)电流节点连接到相邻的一层;（b）电流节点连接到相邻的两层;（c）电流节点连接到相邻的三层。为了避免搜索路径的尖锐锯齿形状，在节点和其周围的节点之间建立有效连接。在传统的A* 搜索算法中，一个节点扩展到相邻层中的26个节点，如图3（a）所示。这意味着当前节点只能向26个方向扩展与当前节点相连的相邻层节点越多，可供路径节点选择的方向就越多，导致搜索路径的准确性越来越高，由于节点间方向的对称性，只绘制了总方向数的八分之一进行搜索。图图3（b）和（c）示出了网格连接到周围的两层（124个节点）和三层（342个节点）。节点的个数Zi与表示为：Zi2i13-1 2在向外扩展的过程中，节点形成的某些方向是重复的，因此可以忽略。通过去除这些网格点可以降低计算为了减少搜索的误差，增加了计算量，为每个网格节点提供更多的连接。DRð Þ公司简介830L。Dong等人/工程6（2020）827采用已划分网格的块体模型和条带模型，探索i层的合理取值，如图4所 示 。 假 设 源 在 点 O 处 被 激 发 ， 并 且 波 从 点 K 到 达 点 K（K^A;B;C;D;E;F;G;H;A0;B0;C0;D0;E0;F0;G0;H0）。O，形成路径LK。根据几何关系，将波浪的实际最快路径距离DR记录到表1中。利用第i层的模型跟踪路径LK（1/41; 2; 3; 4; 5; 6;：：）。搜索到的路径距离DSi和DR之间的相对误差E可以表示为Ei¼. DSi-1×100%3所得路径的相对误差E记录在表1中，并选择每个模型的最大路径误差Ei-max改进的A* 搜索算法的时间复杂度O为正图四、 从源节点到邻居节点的理想路径K（K^A;B;· ··;H0）。与扩展节点Z的数量有关，即Oi Z i。因此，可以使用节点的数量Z来近似计算量O。 图 5显示了数量之间的关系层i和Ei-max以及关联节点的数量Z i。通过不同的测试发现，当拓扑层数达到3时，理论路径的相对误差小于3%。相对误差受传感器尺寸或系统误差的影响。这意味着该算法能够以较小的计算量增长率基本满足定位要求。因此，本文将利用A* 搜索算法来建立节点与周围三层之间的连接。为了确定路径的具体位置，在用A* 搜索算法（或Dijkstra算法）确定最短路径的长度后，还需要进行反向搜索.因此，本文增加了在搜索路径时记录当前节点对应位置上的前一个节点坐标的功能，这也有助于避免反向搜索，提高计算效率。2.3. 收集到达m传感器安装在样品的不同位置，以捕获MS/AE的到达每个传感器都位于一个网格上图五.层数、相对误差和关联节点的数量。表1实际路径距离与搜索路径距离之间的相对误差路径符号DR使用第i层模型搜索路径的相对误差E1/11/21/3四分之一五分之一1/61/71/8LA11.311.760.650.350.200.120.070.030LB11.363.240.920.3000000LC11.494.310.730.090.090.020.020.020LD11.704.880000000LE12.004.811.010.110.110.020.020.020LF12.374.011.400.4900000LG12.812.411.100.630.380.230.130.060LH13.3000000000LA01.7300000000LB02.4511.540000000LC03.3212.534.01000000LD04.2411.544.881.7400000LE05.2010.314.882.320.890000LF06.169.214.622.471.270.51000LG07.148.274.312.451.420.770.3200LH08.127.484.012.371.460.890.500.210Ei-max-12.534.882.471.460.890.500.2100ðÞ0xyzxyzxyz0xyz×××××xyzxyz0L. Dong等/工程6（2020）827-834831点，以便知道位置。对于三维模型，有五个未知数（P波波速、MS/AE源坐标x0;y0;z0、激发初始时间t0），因此传感器数量m应是大于或等于5的整数为第k个传感器Sk接收到信号，坐标记为xk;yk;zk，接收到的P波信号的初至时间为tk。计算两个不同的实际时间差传感器Sl和Sm，由Dtlm表示。2.4. 源位置对于由点Pxyz产生的源，时间tk对应于最短传播路径Lk之间源和第k个传感器除以速度C两个不同的站l和m之间的时间差对应于行进时间差Dtlm。根据Dtlm与Dtlm之差的平方，引入Dxyz来描述点Pxyz与未知MS/AE源P0的偏离程度，表示为D¼X.Dtlm-Dtlm24其中，当源落入空区时，Dxyz¼ 1。每个网格点将获得相应的Dxyz值。点Pxyz与未知源P0的偏离程度随Dxyz值的增大而增大。因此，对应于最小Dxyz值的坐标可以被认为是MS/AE源位置的坐标。3. 结果和讨论通过对空心柱砂浆结构的破铅试验，验证了该方法的有效性。 立方体混凝土试件，的大小为10cm10厘米10厘米被切成一个圆柱形的空的空间尺寸为/6 cm × 10 cm，如图所示。 6（a）. 有10个6立方相同大小的网格块，以形成模型时，双方立方体的长度为1 mm。为了加快计算速度，试件被分成25 25 25立方体，尺寸为4 mm长，如图所示。 6（b）.空白区域的立方体标记为未通过，其他立方体标记为通过。将六个MS/AE传感器连接到测试件，并且在传感器和测试件之间添加偶联剂以获得更好的声学耦合。传感器的位置位于网格节点上（图1）。 6（c）），坐标位置如表2所示。为了便于计算，对模型坐标进行变换，使得节点矩阵M的索引与坐标位置一一对应。数据采集使用40 dB阈值和5 MHz采样率。在试件的不同位置进行引线断裂试验，每个位置（见附录A表S1）进行两次。传感器记录MS/AE事件产生的P波到达（表S1）。搜索每个潜在源和传感器之间的路径L S，并计算路径的距离，如图所示。 6（d）. 然后，根据传感器记录的到达，计算每个潜在源的偏差值Dxyz样品上的源位置的坐标X;Y;Z到的最小Dxyz是确定并且转换为位置坐标。 的定位结果（c）确定传感器的坐标（单位：cm）;（d）用改进的A* 搜索算法搜索到的从节点到传感器的路径（单位：cm）。表2传感器的坐标传感器数量对应坐标（cm）Xyz10.80.810.029.20.810.030.89.210.048.002.0510.08.02.062.010.02.0VFH见表3。同时，执行TD的结果并在表3中示出。VFH和TD之间的定位结果与实际的引线断裂点，以及误差记录如表3所示。从表中可以看出，TD的最大定位误差为4.5 cm，这比VFH的定位误差大得多。为了便于观察，图。 在图7（a）中，圆的大小表示源定位结果的误差的大小。可以清楚地看到，TD的圆比VFH的圆大得多图图7（b）表示两种方法的定位误差E的箱形图。在每个方框上，中心标记、方框底边和方框顶边分别表示中位数、第25百分位数和第75百分位数框上方和下方的须线显示最极端的数据点，离群值使用“+”符号单独绘制如图7（b）所示，通过使用VFH获得的中值源定位误差约为1.0 cm，而通过使用TD获得的中值定位误差约为1.9 cm。显然，VFH的定位效果更好。图7（b）的其他参数也可以验证该观点。根据表中各导线断裂试验的定位误差，VFH的平均定位误差为1.20 cm，TD的平均定位误差为2.02 cm。VFH的平均定位精度比TD高近40%因此，我们认为，图六、与VFH的本地化过程（a）样品;（b）建模和832L. Dong等人/工程6（2020）827表3引线断裂点的定位结果及误差编号VFH TD定位结果（cm）定位误差（cm）定位结果（cm）定位误差（cm）y z4x4y4z E1.9 10.0 1.7-0.1 0 1.71.7 10.0 1.7-0.3 0 1.71.2 9.2 0.8-0.8-0.8 1.41.1 9.4 1.2-0.9-0.6 1.61.1 10.0-0.2-0.9 0 0.90.4 0.4-0.8 1.0 7.2 1.7 10.0 1.2-0.3 0 1.3电话：+86-10- 8888888传真：+86-10 - 88888888电话：+86-10- 8888888传真：+86-10 - 888888880.4-0.4 0 0.6 2.9 4.0 9.8 0.9 0-0.2 0.90 0 0 0 1.9 4.2 9.1-0.1 0.2-0.9 0.9电话：+86-10-8888888传真：+86-10 - 88888888电话：+86-10 - 8888888传真：+86-10 - 888888880.4 0 0 0.4 2.8 5.8 10.0 0.8-0.2 0 0.8电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-10 - 8888888传真：+86-10 - 888888880.4-0.4-1.2 1.3 0.5 7.9 8.2-1.5-0.1-1.8 2.31.2-0.8 0 1.4 2.1 7.6 10.0 0.1-0.4 0 0.4电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-10 - 8888888传真：+86-10 - 8888888电话：+86-10 - 88888888传真：+86-10 - 88888888电话：+86-10 - 8888888传真：+86-10 - 888888881.6 1.6 2.0 3.0 3.4 2.3 10.0 1.4 2.3 2.0 3.3电话：+86-10 -8888888传真：+86-10 - 888888880.8 1.2 1.6 2.2 0.9 0 5.3-1.1 0-0.7 1.3电话：+86-021 - 8888888传真：+86-021 - 8888888电话：+86-10 -8888888传真：+86-10 - 88888888电话：+86-021 - 8888888传真：+86-021 - 8888888电话：+86-021 - 8888888传真：+86-021 - 8888888电话：+86-0511 - 8888888传真：+86-0511 - 88888880.4 0 0.8 0.9 3.4-0.1 3.6-0.6-0.1-0.4 0.7电话：+86-10 - 88888888传真：+86-10 - 88888888电话：+86-0571 - 8888888传真：+86-0571 - 8888888电话：+86-10 - 88888888传真：+86-10 - 88888888电话：+86-021 - 88888888传真：+86-021 - 88888888电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-21 - 8888888传真：+86-21 - 88888888电话：+86-10- 8888888传真：+86-10 - 88888888电话：+86-10- 8888888传真：+86-10 - 88888888电话：+86-0511 - 8888888传真：电话：+86-0511 - 8888888传真：+86-0511 - 8888888电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-0511 - 8888888传真：+86-0511 - 8888888电话：+86-10 - 88888888传真：+86-10 - 88888888电话：+86-10 - 88888888传真：+86-10 - 88888888电话：+86-021 - 88888888传真：+86-021 - 8888888电话：+86-021 - 88888888传真：+86-021 - 8888888电话：+86-21 - 88888888传真：+86-21 - 88888888电话：+86-21 - 6666666传真：+86-21 - 66666666电话：+86-021 - 8888888传真：+86-021 - 88888880.4 0 1.6 1.6 8.7 10.1 10.0 0.7 0.1 2.0 2.1电话：+86-0511 - 8888888传真：+86-0511 - 8888888电话：+86-021 - 8888888传真：+86-021 - 8888888电话：+86-0511 - 8888888传真：+86-0511 - 8888888电话：+86-10 - 8888888传真：+86-10 - 88888888电话：+86-100.4 0 1.2 1.3 6.9 10.1 9.7 0.9 0.1 1.7 1.9电话：+86-10- 8888888传真： +86-10 - 88888888电话：+86-10- 8888888传真： +86-10 - 888888880.4-0.2 0.8 2.2 4.8 10.1 10.0 0.8 0.1 2.0 2.2电话：+86-10 - 8888888传真：+86-10 - 88888888电话：021 - 88888888传真：021 - 888888880 0 0 0 3.1 2.1 8.1-0.1-0.7 0.1 0.7电话：+86-0511 - 8888888传真：+86-0511 - 8888888电话：+86-0511 - 8888888传真：+86-0511 - 8888888电话：+86-0511 - 8888888传真：+86-0511 - 8888888（接下页）Xyz13.62.410.023.62.410.034.01.210.044.41.210.056.02.49.266.42.49.276.82.810.086.82.810.092.43.610.0102.04.010.0117.22.810.0127.23.210.0132.46.010.0142.05.610.0158.05.210.0168.05.610.0172.47.68.8183.27.210.0193.68.010.0203.68.010.0216.07.69.6226.07.69.6237.67.69.6247.67.610.0252.81.210.0263.61.610.0272.81.67.2282.81.27.6292.40.45.2302.81.25.6314.40.89.2324.40.89.6334.004.8344.404.8356.01.68.8366.01.28.0375.61.26.4386.00.86.4396.41.68.0406.82.88.4418.42.86.4426.82.86.84310.02.04.44410.02.04.4459.24.08.0469.64.08.0479.24.04.44810.04.04.8498.86.48.0509.26.48.05110.06.46.85210.06.06.4537.66.03.6547.66.03.6558.410.09.6568.410.09.6578.010.07.2588.010.07.6598.010.04.8607.28.84.0615.68.08.8626.410.09.2635.69.64.8645.69.64.8654.48.08.8664.08.08.4673.22.88.4683.22.88.0693.22.86.8703.22.86.8713.22.84.84x4y4zEX1.6 0.4 01.63.71.6 0.4 01.63.70 - 0.8 00.84.80.4 - 0.8 00.95.20 0.4-0.80.95.8L. Dong等/工程6（2020）827-834833表3（续）编号VFH TD定位结果（cm）定位误差（cm）定位结果（cm）定位误差（cm）x y z4x4y4z E x y z4x4y4z E2019 - 10 -15 00： 00：0073 7.2 6.8 8.0 0 0 0 10.1 8.9 9.3 2.9 2.1 1.3 3.874 7.2 6.8 8.4 0 0 0.4 0.4 10.0 8.7 10.0 2.8 1.9 2.0 3.975 7.2 6.8 6.4 0 0 0.4 0.4 8.0 7.3 6.4 0.8 0.5 0.4 1.076 7.2 6.8 6.8 0 0 0.8 0.8 8.0 7.3 7.1 0.8 0.5 1.1 1.477 7.2 6.8 4.0 0 0 0 7.5 6.3 3.7 0.3-0.5-0.3 0.778 7.2 6.8 5.2 0 0 1.2 1.2 7.8 7.4 4.9 0.6 0.6 0.9 1.2平均误差-0.18 0.05 0.44 1.20 - 0.09 0.31 0.50 2.02图7.第一次会议。两种方法的声发射源定位结果及误差分析（a）结果和误差的可视化;（b）定位误差的箱形图在复杂的三维结构中，VFH的定位精度比TD有很大提高。4. 结论本文提出了一种基于A* 搜索算法的不需要预先测量速度的定位方法。采用等距网格点搜索路径，避免了人工的重复训练，对于空间不规则的三维复杂结构，改进了A* 搜索算法，采用网格点搜索反映弹性波实际传播的路径。此外，在计算中将弹性波的波速视为未知量，以减小监测过程中波速变化引起的定位误差为了评价VFH的准确定性和有效性，在空心立方体混凝土结构上进行了铅断裂试验。定位结果表明，VFH的平均定位误差为1.20 cm，平均定位精度比TD提高了40%。结果表明，该方法能有效地适应几何形状不规则的三维复杂结构。该方法克服了现有定位方法考虑几何影响的缺点，适用于三维含孔结构。本文的方法也适用于声发射定位无损检测等其他领域。确认作者希望感谢国家自然科学基金委员会（51822407和51774327）、湖南省自然科学基金（2018JJ1037）、中南大学创新驱动项目（2020CX 014）。遵守道德操守准则董龙军、胡庆春、童晓杰和刘有芳声明，他们没有利益冲突或财务冲突需要披露。附录A.补充数据本文的补充数据可以在https://doi.org/10.1016/j.eng.2019.12.016上找到。引用[1] 冯X，刘J，陈B，肖Y，冯G，张F. 深部金属矿山岩爆监测预警与控制。 Engineering2017;3（4）：538-45.[2] Milev AM，Spottiswoode SM，Rorke AJ，Finnie GJ.地下隧道模拟岩爆的地震监测。《南非金属研究所杂志》2001;101（5）：253-60.[3] Urbancic TI，Trifu C.加拿大矿山地震监测技术的最新进展。应用地球物理学 报2000;45（4）：225-37.[4] 王宏，葛明.高水平应力石灰石矿声发射/微地震震源定位分析J Rock Mech MinSci2008;45（5）：720-8.[5] 平田A，龟冈Y，平野T.隧道开挖过程中基于声发射监测的岩爆安全管理。RockMechRock Eng 2007;40（6）：563-76.[6] Li L，Tan J，Wood DA，Zhao Z，Becker D，Lyu Q，et al.非常规致密油气藏水力压裂诱发地震监测的现状综述。《燃料》 2019;242：195-210。[7] 通用汽车公司有效的矿山微震监测。国际煤炭地质杂志2005;64（1- 2）：44-56.[8] 达勒姆河减轻南非深层硬岩矿山岩爆风险 上一篇：提取科学：一个世纪的834L. Dong等人/工程6（2020）827采矿研究。丹佛：采矿、冶金和勘探协会;2010. p. 156比71[9] Park B，Sohn H，Olson SE，DeSimio MP，Brown KS，Derriso MM.使用基于激光的时间反转在复杂结构中的冲击定位。结构健康监测2012;11（5）：577-88。[10] Marantidis C，Van Way CB，Kudva JN.军用飞机机载智能结构健康监测系统中的声发射传感。In：Sirkis JS，editor.会议录第2191卷，智能结构和材料1994年：智能传感，处理和仪器; 1994年2月13日至18日;奥兰多，佛罗里达州，美国; 1994年。p.2191号[11] Feng G，Feng X，Chen B，Xiao Y，Yu Y.隧道岩爆发展过程动态预警的微地震方法。Rock Mech RockEng2015;48（5）：2061-76.[12] 冯G，冯X，陈B，肖Y，刘G，张W.锦屏二级水电站深埋隧洞贯通过程中的微震活动特征。Int J Geomech 2020;20（2）：04019163.[13] 郑洁，宋国，孙晓，文玲，李芳。矿山微地震震源定位研究进展与展望。工程2018;4（5）：653-60。[14] 盖革湖仅由地震到达时间确定震中的概率方法。Bull St Louis Univ1912;8：60-71.[15] Inglada V. Die berechnung der uncoordinated eines nahbebens aus denditrittszeiten der in einingen benachbarten stationen aufgezeichneten P-oderP-wellen.Gerlands Beitr Geophys 1928;19：73-98. 德语.[16] 放大图片作者：Leighton F，Blake W.岩石噪声源定位技术。 华盛顿特区：美国内政部情报局，矿业局，1970年。[17] Leighton FW，Duvall WI.改进岩石噪声源定位技术的最小二乘法。华盛顿特区：美国矿业局，1972年。报告编号：BM -R I -7626[18] ThurberCH. 非 线 性 地 震定 位 ： 理 论 与 实 例 。 Bull SeismolSoc Am 1985;75（3）：779-90.[19] 唐湾用计算机确定地震参数的一般方法。《地震学报》1979;1（2）：186-96.中文.[20] PruggerA，Gendzwill D. 微震定位：一种利用单纯形步进法的非线性方法。Bull Seismol Soc Am 1988;78（2）：799-815.[21] Waldhauser F，Ellsworth WL.双差地震定位算法：方法及其在加利福尼亚州北海沃德断层的应用。Bull Seismol Soc Am2000;90（6）：1353-68.[22] 董丽，李新，唐丽，龚飞.无测速微震震源定位的数学函数和参数。中国岩石机械工程杂志2011;30（10）：2057-67. 中文.[23] 李新，董丽. 两种不测量声速的四传感器声发射源定位方法的比较。Sens Lett2011;9（5）：2025-9.[24] 董丽，舒伟，李新，韩刚，邹伟.未知速度挖掘系统中传感器网络源定位的三维综合解析解。IEEE Access2017;5（99）：11337-51.[25] 董丽，李晓，马军，唐丽. 未知速度系声发射/微震源三维解析综合解。 中国岩石机械工程杂志2017;36：186-97. 中文.[26] 董立，舒伟，韩刚，李晓，王军。建筑物内信息物理系统的一种多步变窄速度间隔源定位方法。IEEEAccess 2017;5：20207-19.[27] 董丽，邹伟，李晓，舒伟，王正. 地下采矿岩体结构中微地震/声发射源解析迭代协同定位方法。工程机械2019;210：95-112。[28] BaxterG，Pullin R，Holford KM，Evans SL. 用于声发射的Δ T源定位。机械系统信号处理2007;21（3）：1512-20.[29] EatonMJ，Pullin R，Holford KM. 用Delta T标测法定位复合材料中的声发射源。 Compos Appl Sci Manuf 2012;43（6）：856-63.[30] [10]李文辉，李文辉，李文辉.基于非均匀速度模型的多段路径分析在复杂传播介质声发射源定位中的应用。超声2017;74：48-61。[31] 放大图片作者：HartPE，Nilsson NJ，拉斐尔B. 启发式确定最小费用路径的形式基础。 IEEE Trans Sys Sci Cybern 1968;4（2）：100-7.[32] 放大图片作者：Hart PE，Nilsson NJ，拉斐尔B.修正最小费用路径的启发式确定的形式基础。 ACM Sigart Bull 1972;37（37）：28-9.[33] Oommen BJ ， Rueda LG. A formal analysis of why heuristic functions work.ArtifIntell 2005;164（1-2）：1-22.
我的内容管理
展开
