DOI： 10.11817/j.issn.1672-7207.2020.08.024

二维非测速条件下声发射震源定位方法数值验证

吴顺川^{1, 2},张光¹,张诗淮²,郭沛²,储超群²

（1. 昆明理工大学 国土资源工程学院，云南 昆明， 650093；

2. 北京科技大学 金属矿山高效开采与安全教育部重点实验室，北京，100083）

摘要:基于离散元方法构建平板模型，研究应力波信号传播规律。在平板模型上布设震源点激发应力波作为主动震源，布设监测点记录应力波传播信息，采用3种典型监测点布设方式(等腰直角三角形、一般直角三角形和等边三角形)在非连续介质中进行二维非测速条件下定位试验研究，探讨应力波在三维FJM模型中的传播规律。研究结果表明：1) 互相关技术能够准确计算应力波信号到达相邻监测点的时延，互相关系数可用于评估应力波衰减规律；2) 数值模拟计算所得震源定位结果集中在实际震源周边，且定位误差较小；3) 采用正弦波激发方式所得3种布设方式计算震源位置较优点的比例分别为100.0%，100.0%和87.5%，采用雷克子波激发方式所得3种布设方式计算震源位置较优点的比例分别为100.0%，95.8%和95.8%；4) 在震源处施加集中力，各监测点P波位移场在施加力方向(z方向)上最大，在垂直于施加力方向(x和y方向)上无位移；平板上各监测点实际振幅与位移场远场项理论解相符。

关键词:颗粒流理论；应力波；互相关理论；二维震源定位；定位精度；场源效应

中图分类号:TU45            文献标志码:A             开放科学(资源服务)标识码(OSID)

文章编号:1672-7207（2020）08-2299-13

Numerical verification on two-dimensional acoustic emission source location without knowing velocity profile

WU　Shunchuan^{1, 2}, ZHANG　Guang¹, ZHANG　Shihuai², GUO　Pei², CHU　Chaoqun²

(1. Faculty of Land Resources Engineering, Kunming University of Science and Technology,Kunming 650093, China；

2. Key Laboratory of Ministry of Education for Efficient Mining and Safety of Metal Mine, University of Science and Technology Beijing, Beijing 100083, China)

Abstract: Based on the discrete element method, a plate model was constructed to study the propagation rule of stress wave signals. The source point was established to stimulate stress wave as the active source on the plate model, and the monitoring points were set up to record the stress wave propagation information. Three typical monitoring point layouts (isosceles right triangle, general right triangle and equilateral triangle) were adopted to carry out the source location test under the two-dimensional non-velocity measurement conditions, which was simulated in the discontinuous medium. Further, the propagation law of stress waves in the three-dimensional flat-joint model was explored. The results are as follows. 1) Time delay of the stress wave signals reaching the adjacent receivers can be obtained accurately through cross-correlation technique, and the cross-correlation coefficient analysis can be used to evaluate the attenuation law of stress wave. 2) The source locations obtained by numerical simulation are clustered around the actual sources, and the location error is small. 3) The proportions of the better source location calculated by three typical layouts which are simulated by sine wave are 100.0%, 100.0%, and 87.5% respectively, while those which are simulated by ricker wavelet are 100.0%, 95.8%, and 95.8%, respectively. 4) When a concentrated force is applied in the source, the P wave displacement field is the maximum along the axis of applied force (z direction), and there is no displacement in the direction perpendicular to the applied force(x and y direction). The actual amplitudes of receivers are well matched with the theoretical calculation results of the far-field terms.

Key words: particle flow theory; stress wave; cross-correlation theory; two-dimensional source location; location accuracy; source effects

岩石类材料在受内力或外力作用下发生变形或破裂，在破裂过程中会激发应力波。利用波在传播过程中的时间信息，由时空坐标关系可定位震源。吴顺川等^[1]提出了一种二维非测速条件下声发射震源定位方法，并采用花岗岩平板验证其适用性。声发射定位技术在结构检测、矿山破岩、油井致裂以及隧道掘进等诸多方面得到了广泛应用。定位方法是确定岩石内部缺陷以及破裂损伤位置的关键，国内外已有大量学者研究了定位方法。传统定位方法包括GEIGER法^[2]、单纯形法^[3]和网格搜索法^[4]等。TOBIAS^[5]提出了三角时差定位法；CIAMPA等^[6]提出了无需测量波速的震源定位方法；DONG等^[7-9]提出无需测速的微地震/声发射定位方法，并采用复杂结构模型和原位微震监测验证了定位效果；KUNDU等^[10-11]总结了不同定位方法在各向同性与各向异性材料中的定位差异，并提出了一种特定传感器布设方式的无需测速震源定位方法；YIN等^[12]基于新三角时差定位方法原理，提出了一种新的无需测速震源定位方法，其传感器布设方式为“Z”形。在定位方法的研究基础上，众多学者研究定位精度，LI等^[13]采用花岗岩和大理岩研究速度偏差对线定位和平面定位精度的影响；许江等^[14]采用砂岩试样研究了岩石尺寸效应与声发射定位精度影响因素。由于应力波的传播规律极其复杂，室内试验中波的传播衰减规律难以满足理论条件，数值模拟对应力波进行研究具有独特优势^[15-16]。目前，模拟应力波方法主要包含连续介质理论和非连续介质理论。岩体包含众多断层、节理和裂隙等地质结构，其主要特性为非连续性、非均匀性和各向异性。采用连续介质理论难以准确、客观地描述波在岩体中的传播规律。CUNDALL等^[17]提出的离散元法能够有效模拟岩体的力学特性，反映其非连续性。离散元法适用于计算岩体等非连续性介质的动力学问题，在岩土工程领域应用广泛，国内外众多学者采用离散元模型研究了岩石类材料应力波传播问题。TRENT等^[18]采用离散元法，从细观尺度上研究了介质的波动过程；HAZZARD等^[19]基于颗粒流理论构建砂岩模型，研究了荷载损伤下的波速演化规律；NABIPOUR^[20]应用颗粒流理论建立了六角形排布和方形排布颗粒模型，研究了颗粒施加应力波传播特性，验证了颗粒流模型中应力波传播的合理性；TOOMEY等^[21]应用颗粒流理论构建六角形排布颗粒模型，研究了应力波在大尺度规则排列颗粒模型中的传播规律；SADD等^[22]基于离散元理论研究了枝矢量、接触法向和轴向矢量分布与波速和振幅衰减的关系；张国凯等^[23]采用颗粒流程序建立9种不同各向异性模型，研究了波的传播规律，揭示孔隙率、枝矢量、配位数张量和刚度张量等细观参数对波的传播和衰减规律；张诗淮等^[24]基于颗粒流理论研究应力波在大尺度一维和二维模型中传播，在震源激发方式、颗粒黏结模型、数值弥散和边界条件等方面研究了P波在数值模型中的传播规律；郭易圆等^[25]采用面-面接触模型模拟纵波在一维岩体中的传播，得出了阻尼比、软弱夹层以及节理对波传播规律的影响；徐小敏等^[26]针对砂土等颗粒材料建立离散元模型，分析了激发频率、激发幅值、激发源和接收源尺寸和阻尼等因素对剪切波速的影响；吴顺川等^[27]基于颗粒流理论构建层状岩体模型，在震源定位计算中考虑了层状结构对应力波传播的影响。前人多基于颗粒流理论采用均匀颗粒构建一维和二维模型，在相对简易模型中研究应力波的传播规律，且未进行定量分析。传统定位方法多适用于各向同性介质，对于岩石类非连续性介质，由于内部结构较为复杂，采用传统单一波速模型进行定位会产生误差。吴顺川等^[1]提出的非测速条件下声发射震源定位方法无需预先测速，避免了波速误差对定位结果的影响。室内试验会存在较多的外界干扰因素和人为因素引起的试验误差，数值试验条件下应力波的传播规律更接近于理论条件下的传播规律，采用离散元方法研究应力波具有独特优势。本文基于颗粒流程序构建平板模型，在平板上布设震源点和监测点，采用互相关理论分析监测点的波形衰减规律，确定监测点时延，依据二维非测速震源定位方法确定震源位置并进行误差分析，并基于平板模型探究波动理论在平节理模型(FJM)的适用性。

1  二维非测速定位方法

1.1　二维非测速条件下震源定位方法

根据文献[1]可知，定位监测所用的簇内传感器布设方式可为任意三角形，如图1所示。当在平板上布设多簇传感器时，可对震源进行定位，图1中在平板上布设4簇传感器，传感器簇编号为C_i(i=1，2，3，4)，簇内传感器编号为S_i-j(j=1，2，3)，每簇传感器将S_i-1设定为参考传感器，计算簇内其他2个传感器与该参考传感器的时延分别为t_i-21和t_i-31。

图1　四簇传感器布设示意图^[1]

Fig. 1　Four sets of sensors on a plate^[1]

由文献[1]可知，震源A与参考传感器S_i-1连线的坐标方位角以及超声波信号在该路径的传播速度计算结果如下：

(1)

(2)

式中：为S_i-1S_i-2与x轴正方向夹角；为S_i-1S_i-3与x轴正方向夹角；为S_i-1S_i-2的距离；为S_i-1S_i-3的距离。

由式(1)和(2)可知，坐标方位角以及该传播路径上波速的计算结果准确性取决于t_i-21和t_i-31。依据式(1)求出一簇传感器的震源坐标方位角，2簇传感器联合求解即可确定一个震源坐标。

1.2　时延测量和波形互相关系数分析

数值试验中采用正弦波和雷克子波作为震源输入波形，在平板中设置不同监测点对震源施加的波动信息进行监测，通过计算同一簇监测点的时延确定震源坐标，确定时延的方法主要包括初达波法、峰-峰值法和互相关分析法^[26](cross-correlation，CCR)。采用初达波法和峰-峰值法确定监测点间时延，测量准确性会受到波形扭曲和弥散的影响；采用互相关方法不仅可以确定各监测点时延，还可利用互相关系数评估不同监测点波的传播质量及衰减规律。

对所监测得到的2组波形信号x(t)和y(t)进行处理，其互相关函数^[28-29]定义为

(3)

式中：T为波形时长；为时延。

在对其进行处理中，波形时长会选取固定值，互相关函数估计值为

(4)

当时，函数取最大值，可知波形信号x(t)和y(t)在此处互相关程度最高，其时延为。

对随机信号进行相关性分析时，如果分析整个时间序列信号，较大样本量会造成差异增大，导致互相关系数计算不准确，可采用时间序列延迟相关分析法提高相关系数计算的准确性^[30]。信号x(t)={x₁，，x_t，，x_n}和y(t)={y₁，…，y_t，…，y_n}在时延处相关系数计算公式为

(5)

式中：；

相关系数的取值范围为[-1，1]，越接近于1时，表明信号x(t)和y(t)线性关系越密切；接近于0时，表明信号x(t)和y(t)线性关系越弱。

2  震源定位颗粒流模拟

2.1　离散元模型构建

文献[1]采用花岗岩平板进行震源定位试验，并对非测速条件下声发射震源定位方法进行了验证，为便于与室内试验进行对比，基于颗粒流程序建立的平节理模型^[31-32](flat-joint model, FJM)与室内试验所用花岗岩平板宏观力学参数均一致。花岗岩力学性质如表1所示^[33]。在PFC3D中生成直径为50 mm和高为100 mm的圆柱试样进行单轴压缩试验与直接拉伸试验，调整模型细观参数，将模型的宏观力学性质与实验室结果相匹配，当计算结果与实验室结果相近时，可将该组细观参数应用于模型计算。

表1　花岗岩力学性质^[33]

Table 1　Mechanical properties of granite^[33]

本文采用FJM模型构建花岗岩平板，验证定位方法，表2所示为平节理模型细观参数，当采用表2所示细观参数匹配结果时，所得力学性质与室内试验中花岗岩试样基本一致(见表1)。采用该细观参数构建长为500 mm、宽为500 mm、厚为18 mm的平板模型(如图2所示)，模型中颗粒最小直径为3.0 mm，最大直径为4.5 mm，包含颗粒总数为101 764个。在平板上分别布设震源点和监测点颗粒，对震源颗粒施加脉冲信号，采集监测点应力波信号并进行分析，研究波形衰减以及震源定位结果。

在数值试验中，震源布置位置和监测点布设方式与室内试验布设方式均一致，监测点采用图3所示的I，II和III类3种布设方式，通过数值试验验证非测速条件下声发射震源定位方法的适用性与准确性。

表2　平节理模型细观参数

Table 2　Micro-parameters used in FJM

图2　颗粒流平板模型

Fig. 2　Plate model in particle flow code

图3　监测点布设方式示意图^[1]

Fig. 3　Different sensor arrangements on a plate^[1]

2.2　震源激发方式

对震源颗粒以集中力形式分别施加正弦波(Sine wave)和雷克子波(Ricker wave)脉冲信号，模拟声发射信号。脉冲激发频率为100 kHz^[19]，为保证模型的稳定，激发幅值不宜过大，在该模型中设定振幅为1.0×10^-8 N，每一步运行时长。在震源点(A)的z方向上分别施加正弦波和雷克子波脉冲，在平板上布设监测点记录应力波信息，分析其传播规律，计算各监测点时延进行定位计算。

2.3　应力波振幅衰减及互相关系数分析

在平板上布设多个监测点，研究正弦波和雷克子波在不同震中距的振幅衰减以及不同监测点间的相关性规律。在震源点A与监测点S_3-1(此为第10衰减监测点，S_att-10)之间等间距布设9个监测点(监测点编号为S_att-i，i=1~9)，记录各监测点z方向位移变化信息，如图4所示。

图4　平板模型监测点分布示意图

Fig. 4　receivers placed on the plate

图5所示为AS3-1传播路径上监测点应力波振幅与互相关系数变化规律。为便于对比分析，以第1个衰减监测点(S_att-1)为基准对各监测点振幅归一化处理，互相关系数变化规律为各监测点波形与第1个监测点波形进行互相关处理所得。由图5可见：正弦波和雷克子波初峰振幅随着与震源距离增大而不断降低，且在较近距离范围内会急剧下降，第3个监测点(S_att-3)振幅降低为第1个监测点的50%，第4~10监测点振幅逐渐减小，最终在第10监测点减小至第1监测点的20%以下。

对各监测点波形与第1个监测点波形进行互相关处理得到各监测点的互相关系数。正弦波在第2~6监测点与第1监测点的波形互相关系数缓慢降低，第7~10监测点互相关系数骤降，在第2~10监测点互相关系均大于0.7，可知其与第1监测点波形高度相关。雷克子波在第6~8监测点比第2~5监测点的互相关系数的降低速率快，第9和10监测点的互相关系数与第8监测点相比有所增大，在第2~10监测点互相关系数值均大于0.5，与第1监测点波形显著相关。

图5　AS3-1传播路径上监测点应力波振幅与互相关系数变化规律

Fig. 5　Variations of the stress wave amplitude and cross-correlation coefficient of the receivers on the AS3-1 propagation path

对于2种不同震源激发方式，振幅的衰减规律存在差异。随监测点与震源距离增大，正弦波和雷克子波的振幅均在不断减小，但雷克子波振幅衰减快于正弦波振幅衰减，在远离震源位置的各监测点上正弦波互相关系数大于雷克子波互相关系数。

2.4　震源定位计算及误差分析

2.4.1　时延测量　

由式(1)可知：声发射震源定位计算中，定位结果的准确性主要取决于不同监测点间时延计算。图6所示为一次试验中C₁簇监测点(S_1-1，S_1-2和S_1-3)所记录的应力波信息，3个监测点所采集的波形高度相似，可采用互相关方法对应力波进行处理计算监测点间时延。图7所示为S_1-1，S_1-2和S_1-3监测点所采集的应力波进行互相关处理所得结果，S_1-2和S_1-3分别对S_1-1监测点应力波进行计算可得各监测点间时延，互相关函数峰值处横坐标即为监测点间时延。

图6　不同监测点z方向位移波形

Fig. 6　Waveforms of z-direction displacement in monitoring points

图7　波形信号互相关处理结果

Fig. 7　Cross-correlation of the acoustic emission waveform

2.4.2　震源定位结果　

平板模型上布设4簇监测点(C₁，C₂，C₃和C₄)分别采集1号~4号震源应力波波形信息，依据式(4)对监测点波形信号与参考监测点(S_i-1)进行互相关计算，得到时延并根据式(1)求解震源定位结果。图8所示为采用花岗岩平板进行定位试验中一次定位结果计算示意图，在试验中，传感器布设方式为II类布设，对4号震源进行定位验证^[1]。在该定位试验中，每2簇传感器定位直线交点可确定1个计算定位点，采用4簇传感器进行震源定位，可计算得到6个定位结果。

图8　花岗岩平板4号震源定位计算示意图

Fig. 8　Location calculation of 4号 source on granite plate

图9和图10分别为对震源施加正弦波脉冲和雷克子波脉冲模拟声发射信号进行试验所得定位结果，其中，监测点布设方式采用图3所示的3种布设方式(I，II和III类)。分析定位结果精度时，引入较优点概念^[14]，即定位结果与实际震源距离相对误差小于20%的点为较优点，定位结果与实际震源距离相对误差和定位直线与实际震源距离相对误差统计见表3。

由图9可见：1~4号震源定位试验中定位结果大部分集中在实际震源附近，落在较优点范围内，但仍有部分定位结果误差较大。采用III类布设方式时，1号和4号震源定位试验中C₂-C₃簇监测点定位结果误差较大，2号震源定位试验中存在1个较大定位偏离点(图9(b)虚线椭圆标记处)。2号震源与C₁簇和C₄簇监测点连线线段之间距离较小，C₁簇传感器与2号震源连线的坐标方位角与2号震源与C₄簇传感器连线的坐标方位角相近，在定位计算所得的坐标方位角的微小误差会导致定位点的较大误差，该误差产生原因见文献[1]。

由图10可见：1号震源采用II类布设方式和4号震源采用III类布设方式所得定位结果偏离实际震源较远，超出较优点范围(即相对误差大于20%)，其余定位结果较为集中，均在较优点范围内。

由图9和10所示的震源定位结果和表3中定位结果相对误差统计可知，2种震源激发方式(正弦波和雷克子波)均较好地验证了该震源定位方法在非连续介质的适用性与准确性。

图9　正弦波激发方式震源定位结果

Fig. 9　Source location results by the type of sine wave excitation

图10　雷克子波激发方式震源定位结果

Fig. 10　Source location results by the type of ricker wave excitation

表3　数值试验平板定位结果相对误差统计

Table 3　Relative error of source location on the plate from the numerical simulation

表4所示为文献[1]中采用花岗岩平板进行室内试验所得定位结果相对误差统计结果，表5所示为数值模拟定位结果相对误差统计结果。在数值试验定位结果中，采用正弦波激发震源，3种监测点布设方式(I类、II类和III类)所得震源定位结果较优点比例分别为100.0%，100.0%和87.5%；采用雷克子波激发震源，3种监测点布设方式(I类、II类和III类)所得震源定位结果较优点比例分别为100.0%，95.8%和95.8%。与室内试验定位结果相比，数值试验定位结果相对误差较小，定位准确性较高。在定位试验中，室内试验中影响因素更为复杂，花岗岩平板中存在的固有缺陷会对定位结果造成误差，而数值试验干扰因素较少，应力波传播更接近理论条件，故数值试验定位结果较优点比例明显高于室内试验较优点比例。

数值试验中颗粒粒径及位置为随机生成，颗粒在空间上非均匀排布。定位计算时，监测颗粒z方向坐标差异会对定位结果带来一定误差，且该定位误差无法通过计算消除。由于该二维定位方法并未考虑监测点z方向影响，定位方法理论计算中仅包含x及y方向，监测颗粒在实际z方向上会存在差异，导致时延测量出现误差。

表4　室内试验平板定位结果较优点统计^[1]

Table 4　The better points counts from the laboratory test^[1]

表5　数值试验平板定位结果较优点统计

Table 5　The better points counts from the numerical simulation

3  P波辐射花样定量分析

文献[24]定性验证了应力波在接触模型中的合理性，由于接触模型无法有效反映岩石力学特性，基于FJM模型，进一步定量对比应力波辐射花样与理论解的差异。

模型中弹性波由点震源激发，根据弹性波动力学理论，其三维波动矢量方程及其格林函数解^[24,34]如下：

(6)

(7)

式中：(r，)为质点位移；r为辐射距离矢量；t为时间；c为波速；F(r，)为源项；为位移分量；i为位移方向(i=x，y和z)；j为作用力方向(j=x，y和z)；和分别为震源与监测点的连线矢量与坐标轴的方向余弦；r为震源与监测点距离；T(t)为震源激励函数；为P波波速；为狄拉克函数；为S波波速；为密度。

在平板模型震源点以力形式在z方向施加脉冲波形。式(7)中为近场项，为P波远场项，为S波远场项。在实践中震源机制研究中主要取决于位移场远场项，因此，在计算中可忽略近场项^[35-36]，并将远场项中的常数项取1，则有

(8)

(9)

式中：方向余弦为，，。

图11所示为P波和S波三维辐射花样。由图11可知：P波在平行于施加力方向上(z方向)，位移最大，而在垂直于施加力方向上(x和y方向)，位移为0；S波在垂直于施加力方向上(x和y方向)，位移最大，而在平行于施加力方向上(z方向)，位移为0。在x-z二维坐标系中，P波和S波辐射花样如图12所示，其位移分量如下：

图11　z方向集中力激发震源P波和S波震动位移场远场项三维空间辐射花样分布

Fig. 11　Three-dimensional representations of the far-field P- and S-wave radiation patterns generated by a point source in the z direction

，(10)

，(11)

P波和S波的振幅表达式为

，(12)

图12　z方向集中力激发震源P波和S波震动位移场远场项x-z平面辐射花样分布^[24]

Fig. 12　Two-dimensional representations of the far-field P-wave and S-wave radiation patterns generated by a point source in the z direction^[24]

图13所示为平板模型运行400 时步(2×10^-5 s)颗粒位移矢量图。在震源z方向施加雷克子波，观察震源周边颗粒位移变化，颗粒位移振动均为z方向。图13中A范围为P波波峰，B范围为波谷，波形由震源点向四周传播时，颗粒位移方向均与施加力方向(z方向)一致，在x和y方向上无位移。应力波在x-y平面上以圆周形式向外传播，平板边界为自由边界，未设置吸收边界，C范围为平板边界，故该范围处颗粒z方向位移明显大于平板内部颗粒位移，应力波传播至平板边界会进行反射，入射波与反射波在此处位移矢量叠加，该处位移明显大于其他传播方向位移，该模拟结果与张诗淮等^[24]研究的边界条件对应力波传播振幅影响规律一致。

图13　2×10^-5 s时刻平板模型颗粒位移矢量图

Fig. 13　Displacement vectors in the plate model at 2×10^-5 s

图14所示为震源A与衰减监测点S_att-5的位移监测波形信息图，为便于比较，将震源点与监测点位移波形信息振幅进行归一化处理，其中z方向位移分量显著大于x和y方向位移分量，与图11和图12中辐射花样相符。P波传播位移方向均平行于施加力方向(z方向)，在垂直于施加力方向上无位移。

图14　震源点与监测点(S_att-5) x，y和z方向振幅对比

Fig. 14　The amplitude of x, y, and z direction displacement from the source and receiver(S _att-5)

在AS_3-1传播路径上布设的各监测点与震源之间距离(r)和波长(λ)之间关系均满足r>>λ，可知在该研究中位移场主要取决于远场项^[36]，忽略近场项影响，由式(7)可得，各监测点P波在z方向位移如下：

(13)

式中：为由10个监测点所得的P波速度，取4 335.3 m/s，ρ=2 800 kg/m³，如表2所示。震源激励函数T(t)如下：

(14)

式中：A₀取1.0×10^-8 N，f取100 kHz。

若以集中力作为激发震源，式(13)的理论解如图15(a)所示。由图15(a)可见：在距离震源较近时，位移振幅随r增大迅速衰减；在距离震源较远时，衰减速率显著变缓，两者呈反比关系。进一步，将AS_3-1传播路径上布设的监测点位移振幅与式(13)理论解进行对比，如图15(b)所示，监测点实际振幅与理论解基本一致，证明了波动理论在FJM模型中的适用性与准确性。

图15　z方向位移振幅衰减性质

Fig. 15　Attenuation property of z-direction displacement

4  结论

1) 采用互相关技术处理声发射信号，可准确得到相邻监测点时延，根据时延反演震源定位结果；随着监测点与震源点间距离增大，各监测点振幅及互相关系数均不断减小。

2) 计算定位点均分布在实际震源周边，相对误差较小，验证了二维非测速下震源定位方法在非连续介质中的适用性。

3) 在正弦波激发方式下，等腰直角三角形，一般直角三角形和等边三角形3种监测点布设方式定位结果较优点比例分别为100.0%，100.0%和87.5%；在雷克子波激发方式下，3种监测点布设方式定位结果较优点比例分别为100.0%，95.8%和95.8%。

4) 应力波信号在平板模型x-y平面上以圆周形式向外传播，各监测点P波位移在z方向上最大，在x方向和y方向无位移。平板各监测点实际振幅与理论解基本一致，证明了波动理论在FJM模型中的适用性与准确性。

参考文献：

[1] 吴顺川, 张光, 张诗淮, 等. 二维非测速条件下声发射震源定位方法试验研究[J]. 岩石力学与工程学报, 2019, 38(1): 28-39.

WU Shunchuan, ZHANG Guang, ZHANG Shihuai, et al. Experimental research on two-dimensional acoustic emission source location without knowing the velocity profile[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(1): 28-39.

[2] GEIGER L. Probability method for the determination of earthquake epicenters from the arrival time only[J]. Bulletin of Saint Louis University, 1912, 8(1): 60-71.

[3] NELDER J A, MEAD R. A simplex method for function minimization[J]. The Computer Journal, 1965, 7(4): 308-313.

[4] SAMBRIDGE M S, KENNETT B L N. A novel method of hypocentre location[J]. Geophysical Journal International, 1986, 87(2): 679-697.

[5] TOBIAS A. Acoustic-emission source location in two dimensions by an array of three sensors[J]. Non-Destructive Testing, 1976, 9(1): 9-12.

[6] CIAMPA F, MEO M. A new algorithm for acoustic emission localization and flexural group velocity determination in anisotropic structures[J]. Composites Part A: Applied Science and Manufacturing, 2010, 41(12): 1777-1786.

[7] DONG Longjun, ZOU Wei, LI Xibing, et al. Collaborative localization method using analytical and iterative solutions for microseismic/acoustic emission sources in the rockmass structure for underground mining[J]. Engineering Fracture Mechanics, 2019, 210: 95-112.

[8] 董陇军, 李夕兵, 唐礼忠, 等. 无需预先测速的微震震源定位的数学形式及震源参数确定[J]. 岩石力学与工程学报, 2011, 30(10): 2057-2067.

DONG Longjun, LI Xibing, TANG Lizhong, et al. Mathematical functions and parameters for microseismic source location without pre-measuring speed[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(10): 2057-2067.

[9] HU Qingchun, DONG Longjun. Acoustic emission source location and experimental verification for two-dimensional irregular complex structure[J]. IEEE Sensors Journal, 2020, 20(5): 2679-2691.

[10] KUNDU T. Acoustic source localization[J]. Ultrasonics, 2014, 54(1): 25-38.

[11] KUNDU T, NAKATANI H, TAKEDA N. Acoustic source localization in anisotropic plates[J]. Ultrasonics, 2012, 52(6): 740-746.

[12] YIN Shenxin, CUI Zhiwen, KUNDU T. Acoustic source localization in anisotropic plates with “Z” shaped sensor clusters[J]. Ultrasonics, 2018, 84: 34-37.

[13] LI Qiyue, DONG Longjun, LI Xibing, et al. Effects of sonic speed on location accuracy of acoustic emission source in rocks[J]. Transactions of Nonferrous Metals Society of China, 2011, 21(12): 2719-2726.

[14] 许江, 李伊, 田傲雪, 等. 声发射定位精度尺寸效应的试验研究[J]. 岩石力学与工程学报, 2016, 35(S1): 2826-2835.

XU jiang, LIYI, TIAN Aoxue, et al. Experimental research on size effect of acoustic emission location accuracy[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(S1): 2826-2835.

[15] 田振农, 李世海, 肖南, 等. 应力波在一维节理岩体中传播规律的试验研究与数值模拟[J]. 岩石力学与工程学报, 2008, 27(S1): 2687-2693.

TIAN Zhennong, LI Shihai, XIAO Nan, et al. Experimental studies and numerical simulation of stress wave propagation in one-dimensional rock mass[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(S1): 2687-2693.

[16] 李建春, 李海波, MA Guowei, 等. 节理岩体的一维动态等效连续介质模型的研究[J]. 岩石力学与工程学报, 2010, 29(S2): 4063-4067.

LI Jianchun, LI Haibo, MA GUOWEI, et al. An equivalent 1d dynamic continuum model for rock mass with parallel joints[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(S2): 4063-4067.

[17] CUNDALL P A. A computer model for simulating progressive large scale movements in blocky rock systems[C]// Proceedings of the Symposium of International Society of Rock Mechanics, France, Nancy: 1971: 129-136.

[18] TRENT B C, MARGOLIN L G. A numerical laboratory for granular solids[J]. Engineering Computations, 1992, 9(2): 191-197.

[19] HAZZARD J F, YOUNG R P. Numerical investigation of induced cracking and seismic velocity changes in brittle rock[J]. Geophysical Research Letters, 2004, 31: L01604

[20] NABIPOUR A. Experimental and numerical study of ultrasonic monitoring of hydraulic fracture propagation[D]. Perth: Department of Petroleum Engineering in Curtin University, 2013: 46-58.

[21] TOOMEY A, BEAN C J. Numerical simulation of seismic waves using a discrete particle scheme[J]. Geophysical Journal International, 2000, 141(3): 595-604.

[22] SADD M H, GAO J, SHUKLA A. Numerical analysis of wave propagation through assemblies of elliptical particles[J]. Computers and Geotechnics, 1997, 20(3/4): 323-343.

[23] 张国凯, 李海波, 夏祥, 等. 岩石波速与损伤演化规律研究[J]. 岩石力学与工程学报, 2015, 34(11): 2270-2277.

ZHANG Guokai, LI Haibo, XIA Xiang, et al. Wave velocity and damage development of rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(11): 2270-2277.

[24] 张诗淮, 吴顺川, 陈子健. 低频动载应力波传播规律及颗粒流模拟方法研究[J]. 岩石力学与工程学报, 2016, 35(8): 1555-1568.

ZHANG Shihuai, WU Shunchuan, CHEN Zijian. Stress wave propagation under low frequency dynamic loading and simulation method with particle flow code[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(8): 1555-1568.

[25] 郭易圆, 李世海. 有限长岩柱中纵波传播规律的离散元数值分析[J]. 岩石力学与工程学报, 2002, 21(8): 1124-1129.

GUO Yiyuan, LI Shihai. Distinct element analysis on propagation characteristics of p-wave in rock pillar with finite length[J]. Chinese Journal of Rock Mechanics and Engineering, 2002, 21(8): 1124-1129.

[26] 徐小敏, 凌道盛, 黄博, 等. 离散元模拟中颗粒材料剪切波速的剪切振动确定方法[J]. 岩土工程学报, 2011, 33(9): 1462-1468.

XU Xiaomin, LING Daosheng, HUANG Bo, et al. Determination of shear wave velocity in granular materials by shear vibration within discrete element simulation[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(9): 1462-1468.

[27] 吴顺川, 陈子健, 张诗淮, 等. 缓倾地层微震定位算法探讨及其数值验证[J]. 岩土力学, 2018, 39(1): 297-307.

WU Shunchuan, CHENG Zijian, ZHANG Shihuai, et al. Microseismic location algorithm for gently inclined strata and its numerical verification[J]. Rock and Soil Mechanics, 2018, 39(1): 297-307.

[28] 李孟源, 尚振东, 蔡海潮, 等. 声发射检测及信号处理[M]. 北京: 科学出版社, 2010: 52-56.

LI Mengyuan, SHANG Zhendong, CAI Haichao, et al. Acoustic emission detection and signal processing[M]. Beijing: Science Press, 2010: 52-56.

[29] 李灏, 陈励军. 互相关法时延估计中的多径抑制方法[J]. 声学技术, 2013, 32(5): 421-425.

LI Hao, CHEN Lijun. The multipath effect eliminating technique for time delay estimation of cross-correlation method[J]. Technical Acoustics, 2013, 32(5): 421-425.

[30] SAKURAI Y, PAPADIMITRIOU S, FALOUTSOS C. Braid: stream mining through group lag correlations[C]//Proceedings of The 2005 ACM SIGMOD International Conference on Management of Data. New York: ACM, 2005: 599-610.

[31] POTYONDY D O. PFC3D flat joint contact model version1[R]. Minneapolis: Itasca Consulting Group, 2013: ICG7234-L.

[32] WU Shunchuan, XU Xueliang. A study of three intrinsic problems of the classic discrete element method using flat-joint model[J]. Rock Mechanics and Rock Engineering, 2016, 49(5): 1813-1830.

[33] 吴顺川, 郭沛, 张诗淮, 等. 基于巴西劈裂试验的花岗岩热损伤研究[J]. 岩石力学与工程学报, 2018, 37(S2): 3805-3816.

WU Shunchuan, GUO Pei, ZHANG Shihuai, et al. Study on thermal damage of granite based on Brazilian splitting test[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(S2): 3805-3816.

[34] PUJOL J. Elastic wave propagation and generation in seismology[M]. Cambridge, UK: Cambridge University Press, 2003: 278-315.

[35] 曹安业, 窦林名. 采场顶板破断型震源机制及其分析[J]. 岩石力学与工程学报, 2008, 27(S2): 3833-3839.

CAO Anye, DOU Linming. Analysis of focal mechanism caused by rupture of stope roof[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(S2): 3833-3839.

[36] 冯俊军, 王恩元, 沈荣喜, 等. 煤体压剪破裂震源模型及远场震动特征[J]. 中国矿业大学学报, 2016, 45(3): 483-489.

FENG Junjun, WANG Enyuan, SHEN Rongxi, et al. A coal shear fracture source model and its far-field seismic characteristics[J]. Journal of China University of Mining & Technology, 2016, 45(3): 483-489.

（编辑  秦明阳）

收稿日期： 2019 -11 -18; 修回日期： 2020 -01 -14

基金项目(Foundation item)：国家自然科学基金资助项目(51774020, 51934003); 云南省高校科技创新团队支持计划资助 (Projects(51774020, 51934003) supported by the National Natural Science Foundation of China; Project supported by Program for Innovative Research Team (Science and Technology) in University of Yunnan Province)

通信作者：张光，博士研究生，从事岩石力学试验和数值模拟研究。E-mail：zhangguang09@163.com