# The Influence Mechanism and Optimization of the Sensor Network on the MS/AE Source Location

Authors: Nan Li; Maochen Ge; Enyuan Wang; Shaohua Zhang
Journal: Shock and Vibration
 (2020)

---

## Abstract

The sensor network layout is a key factor affecting the accuracy and stability of the microseismic/acoustic source (MS/AE) location. Based on the arrival-time-difference principle, the hyperbolic/hyperboloidal governing equations for the source location are derived. The nonuniformity geometrical characteristics of hyperbolic/hyperboloidal field for the source location are obtained. The sensor network does not induce any location errors; it only affects the source location accuracy by amplifying the existing errors in the input data during the source location process. Also, this amplication effect of the input data errors is characterized by nonuniformity because of the nonuniformity of the hyperbolic/hyperboloidal field. Furthermore, two basic effects, the geometrical spreading and the directional control, of the sensor network are investigated, and the three-dimensional space quantitative models of these two effects are established, respectively. The influence of the wave velocity error and arrival time error on the source location accuracy is analytically compared, and the propagation characteristics of these two types of errors during the source location process are revealed. The concepts of critical arrival-time difference and critical hyperbola/hyperboloid are proposed. Based on these two concepts, the monitoring area can be divided into two regions where the source location accuracy is controlled by the velocity error and the arrival time error, respectively. The concept of direction angle of paired sensors is proposed, and the relationship between the source location and the layout of four typical paired sensors is studied. Finally, the principles of sensor network optimization are determined.


---

## Body

## 1. Introduction

It is generally accepted today that rock emits seismic signals when it is stressed and/or deformed. The microseismic (MS)/acoustic emission (AE) technique refers to a geophysical monitoring technology that studies the internal structure and stability, fracture mechanism, and morphology of rock mass by using the seismic signals. In the recent years, the MS/AE technique has been widely used in monitoring and early warning of coal and rock dynamic disasters, hydraulic fracture (HF) monitoring, and stability monitoring of underground structures such as mines, tunnels, and natural gas, and petroleum storage caverns, as well as surface structures such as dams, slopes, and foundations [1–3]. A high accuracy source location is the first step in studies of the MS/AE technique. The factors that influence the source location accuracy may be divided into two types: input data errors and systematic bias [4, 5]. The source location accuracy is definitely affected if there are errors in the input data such as the velocity model and the arrival time. The influence of the systematic error is mainly caused by the sensor network layout. Some effects of the sensor network relative to location accuracy have been noted; for example, source location accuracy inside a sensor network is usually better than that outside this network. But, the general influence of the sensor network is still not clear.Because the errors associated with the input data are unavoidable, the realistic solution is to reduce the impact of these input errors. So, many advanced source location algorithms and optimization analysis methods are developed, such as Geiger’s method, genetic algorithms, the simplex method, the double-difference method, and the particle swarm optimization method [6–10]. On the other hand, the stability and accuracy of the source location are also dependent on the sensor network layout. The sensor network geometry and its mechanism in the control of source location are critical problems needing study in detail [3, 11]. Some people studied the source location efficiency for particular sensor network geometries, such as triangle, quadripartite, and hexagon [12–14]. Kijko proposed a general method to the distribution of seismic sensors based on the D-optimality theory. Also, he developed a procedure for obtaining the optimal spatial distribution of MS sensors in mines by using D-optimality theory [15–17]. The D-optimality method and its variants were widely applied in natural and induced earthquakes [18–22]. Based on the D-optimization theory, Gong et al. [4, 23] proposed a principle to choose the optimal channel numbers and created a model of the expected epicenter and hypocenter error by using the numerical emulation experiment. Finally, they built an optimal design and assessment system of the MS network. Mendecki developed a method for optimal design of mine seismic networks based on the C-optimality theory [24]. The research results of Jia and Li indicated that the source location accuracy depends upon the geometry between the sources and MS sensors when the wave velocity model is determined [25]. An optimization of the MS network was carried out in the Dongguashan Copper Mine by Tang et al. [21], and they designed some different schemes of the MS network. The optimum MS network is obtained by calculating the errors of source location and system sensitivities associated with these different networks. Wang et al. studied the relationship between the MS station coordinates deviation and source location accuracy, and results show that the source location errors increase with errors of the MS sensor station [26]. Jia et al. found that the impact of the sensor density on the source location accuracy can be divided into two parts: the number of sensors and sensor plane quantity [27]. Stabile et al. developed a method for evaluating network performance in surface and borehole seismic monitoring. This method can estimate the location and origin time errors of seismic events with different magnitudes by a specified network geometry [28]. According to Stabile’s method, Landro et al. evaluated the performance of sensor networks with varying the station density and geometry in the induced MS monitoring. The results show that the station density and the noise level represent the crucial parameters for the seismic network performances [29]. LI et al. studied the influence law of the MS/AE source location affected by the sensor network, velocity model, and arrival time by laboratory tests [5, 11]. They also studied two types of multiple solutions related to the sensor network for the source location if the arrival-time-difference principle is used [30].In summary, previous studies mainly focused on the optimization and evaluation of the specific sensor network. For example, the performance of a network is usually evaluated after installation by plotting the contour map of the location error in the monitoring area. However, compared to the importance of sensor network geometry, it is remarkable how little attention has been paid to the effect of the sensor network on the MS/AE source location. Therefore, based on the space geometry theory and arrival-time-difference principle, we first demonstrated that why and how the sensor network has such an important influence on the MS/AE source location. Then, propagation characteristics of input data errors, such as wave velocity error and arrival time error, are determined. Finally, the optimization principle of the sensor network layout is proposed.


## 2. Hyperbolic/Hyperboloidal Governing Equation of the Source Location and Its Geometrical Distribution Characteristics

### 2.1. The Hyperbolic/Hyperboloidal Governing Equation

SupposeSx0,y0,z0,t0 and Tixi,yi,zi,ti are the MS or AE source and the ith sensor, respectively. It is very difficult to identify S waves because the MS and AE monitoring area is small. So, P-wave velocity is the most commonly used velocity in the MS and AE source location. Also, the source location equation based on the arrival-time-difference principle can be expressed as(1)xi−x02+yi−y02+zi−z02=vPti−t0,where x0, y0, z0 and xi, yi, zi are the coordinates of the source and the ith sensor, respectively. t0 and ti are the origin time of the source and the observed arrival time of the P-wave at the ith sensor. vP is the P-wave velocity. Also, i=1,2,…,N, and here, N denotes the number of this equation.Let the observed arrival time of any two sensors,Ti and Tj, be ti and tj, and ti < tj. Subtracting the ith equation from the jth equation, which are all defined by system of equation (1), gives the source location-governing equation [30]:(2)xj−x02+yj−y02+zj−z02−xi−x02+yi−y02+zi−z02=vtj−ti.According to the geometry theory, equation (2) represents a hyperboloid with the sensors Ti and Tj as the focal points, which indicates that the path of the source locations determined by two sensors is a hyperboloid in space. For a two-dimensional source location, the source location-governing equation (2) reduces to(3)xj−x02+yj−y02−xi−x02+yi−y02=vptj−ti.Equation (3) represents a hyperbola, with sensor Ti and Tj as the foci, in a plane. According to analytical geometry, the difference of distances of any point on this hyperbola from its foci is a constant, say 2aij, which satisfies 2aij=vptj−ti. Therefore, the path of the potential source determined by two sensors in the plane is a hyperbola. In other words, any point on this hyperbola may be the potential source. So, the geometrical meaning of the arrival-time-difference source location method is looking for the intersection point of different hyperbolas. As shown in Figure 1, two hyperbolas generated by three sensors in a plane and the intersection point of hyperbola I and II is the source.Figure 1
MS/AE source determined by two hyperbolas on a plane.


### 2.2. Nonuniformity Geometrical Characteristics of the Hyperbolic Field of the Source Location

According to equation (3), the distance difference of any point on the hyperbola focusing on the sensors Ti and Tj from its foci is a constant which is equal to 2aij, and 2aij=vtj−ti. Thus, the position and the shape of a hyperbola just depend on the value of aij when the positions of the sensors (foci) are determined. Based on the geometry of a hyperbola and the arrival-time-difference source location principle, 2aij satisfies the following relationship:(4)0≤2aij≤2cij,where 2cij is the distance between the foci.As shown in Figure2, when the value of 2aij is equal to zero, the source locates on the perpendicular bisector (denotes as line L) of the line jointing the two sensors. When 2aij is 2cij, the source is on a focus or on the outside of the line jointing the foci. Also, a hyperbola field is formed (Figure 2) when the value of 2aij increases from 0 to 2cij. In general, the smaller the value of 2aij, the closer the hyperbola to the line L. Besides, the density and the directions of hyperbola in this field are different.Figure 2
Hyperbolic field associated with two MS sensorsTi and Tj.Based on the source location-governing equation, the source is always located somewhere on one of the hyperbolas whose value of the real axis is equal to2aij=vtj−ti. It means that the shape and position of a hyperbola that contains the source is affected by the parameters of arrival time and wave velocity. However, in practice, 2aij always deviates from its true value because there are always some errors in arrival time and wave velocity. If the input data error is assumed equal to ±δa, the calculating source will jump to a place on the adjacent hyperbola which controlled by the value of 2aij±δa.As shown in Figure2, there are five different sources in the hyperbolic field, and these five sources will incorrectly locate on the neighbouring hyperbolas if there are errors in the input data. For the same input error, say δa=0.2c, it clearly shows that the location errors are different of these five sources. Generally speaking, the error increases as the sources are located further from the origin point O. In order to get good location accuracy, the real source need be at a point where the hyperbola corresponding to 2aij and 2aij±δa are close enough. As a result, the location error caused by the “source jumping effect” associated with the sensor network can be controlled. Otherwise, the location error will be amplified. Therefore, the higher the density of hyperbolic field in the vicinity of the source is, the better the source location accuracy there will be. Unfortunately, the hyperbolic field associated with sensor network is nonuniform, and the density of hyperbolas varies from point to point in the hyperbolic field. It makes the source location accuracy heavily dependent on the relative position of the true source and the sensor network layout.From the abovementioned analysis, some important conclusions associated with sensor network layout and hyperbolic field are summarized as follows: the sensor network layout itself does not cause any location errors; it only amplifies input data errors already existed. The essence of this amplification in terms of hyperbolic geometry is that the calculated source is incorrectly located on an adjacent hyperbola rather than on the one associated with the true source. In addition, the effect of this amplification is not uniform, and this nonuniform amplification effect is due to the nonuniformity of the hyperbolic field. The results indicate that the location error is highly dependent on the position of the true source relative to the sensor network layout.


## 2.1. The Hyperbolic/Hyperboloidal Governing Equation

SupposeSx0,y0,z0,t0 and Tixi,yi,zi,ti are the MS or AE source and the ith sensor, respectively. It is very difficult to identify S waves because the MS and AE monitoring area is small. So, P-wave velocity is the most commonly used velocity in the MS and AE source location. Also, the source location equation based on the arrival-time-difference principle can be expressed as(1)xi−x02+yi−y02+zi−z02=vPti−t0,where x0, y0, z0 and xi, yi, zi are the coordinates of the source and the ith sensor, respectively. t0 and ti are the origin time of the source and the observed arrival time of the P-wave at the ith sensor. vP is the P-wave velocity. Also, i=1,2,…,N, and here, N denotes the number of this equation.Let the observed arrival time of any two sensors,Ti and Tj, be ti and tj, and ti < tj. Subtracting the ith equation from the jth equation, which are all defined by system of equation (1), gives the source location-governing equation [30]:(2)xj−x02+yj−y02+zj−z02−xi−x02+yi−y02+zi−z02=vtj−ti.According to the geometry theory, equation (2) represents a hyperboloid with the sensors Ti and Tj as the focal points, which indicates that the path of the source locations determined by two sensors is a hyperboloid in space. For a two-dimensional source location, the source location-governing equation (2) reduces to(3)xj−x02+yj−y02−xi−x02+yi−y02=vptj−ti.Equation (3) represents a hyperbola, with sensor Ti and Tj as the foci, in a plane. According to analytical geometry, the difference of distances of any point on this hyperbola from its foci is a constant, say 2aij, which satisfies 2aij=vptj−ti. Therefore, the path of the potential source determined by two sensors in the plane is a hyperbola. In other words, any point on this hyperbola may be the potential source. So, the geometrical meaning of the arrival-time-difference source location method is looking for the intersection point of different hyperbolas. As shown in Figure 1, two hyperbolas generated by three sensors in a plane and the intersection point of hyperbola I and II is the source.Figure 1
MS/AE source determined by two hyperbolas on a plane.


## 2.2. Nonuniformity Geometrical Characteristics of the Hyperbolic Field of the Source Location

According to equation (3), the distance difference of any point on the hyperbola focusing on the sensors Ti and Tj from its foci is a constant which is equal to 2aij, and 2aij=vtj−ti. Thus, the position and the shape of a hyperbola just depend on the value of aij when the positions of the sensors (foci) are determined. Based on the geometry of a hyperbola and the arrival-time-difference source location principle, 2aij satisfies the following relationship:(4)0≤2aij≤2cij,where 2cij is the distance between the foci.As shown in Figure2, when the value of 2aij is equal to zero, the source locates on the perpendicular bisector (denotes as line L) of the line jointing the two sensors. When 2aij is 2cij, the source is on a focus or on the outside of the line jointing the foci. Also, a hyperbola field is formed (Figure 2) when the value of 2aij increases from 0 to 2cij. In general, the smaller the value of 2aij, the closer the hyperbola to the line L. Besides, the density and the directions of hyperbola in this field are different.Figure 2
Hyperbolic field associated with two MS sensorsTi and Tj.Based on the source location-governing equation, the source is always located somewhere on one of the hyperbolas whose value of the real axis is equal to2aij=vtj−ti. It means that the shape and position of a hyperbola that contains the source is affected by the parameters of arrival time and wave velocity. However, in practice, 2aij always deviates from its true value because there are always some errors in arrival time and wave velocity. If the input data error is assumed equal to ±δa, the calculating source will jump to a place on the adjacent hyperbola which controlled by the value of 2aij±δa.As shown in Figure2, there are five different sources in the hyperbolic field, and these five sources will incorrectly locate on the neighbouring hyperbolas if there are errors in the input data. For the same input error, say δa=0.2c, it clearly shows that the location errors are different of these five sources. Generally speaking, the error increases as the sources are located further from the origin point O. In order to get good location accuracy, the real source need be at a point where the hyperbola corresponding to 2aij and 2aij±δa are close enough. As a result, the location error caused by the “source jumping effect” associated with the sensor network can be controlled. Otherwise, the location error will be amplified. Therefore, the higher the density of hyperbolic field in the vicinity of the source is, the better the source location accuracy there will be. Unfortunately, the hyperbolic field associated with sensor network is nonuniform, and the density of hyperbolas varies from point to point in the hyperbolic field. It makes the source location accuracy heavily dependent on the relative position of the true source and the sensor network layout.From the abovementioned analysis, some important conclusions associated with sensor network layout and hyperbolic field are summarized as follows: the sensor network layout itself does not cause any location errors; it only amplifies input data errors already existed. The essence of this amplification in terms of hyperbolic geometry is that the calculated source is incorrectly located on an adjacent hyperbola rather than on the one associated with the true source. In addition, the effect of this amplification is not uniform, and this nonuniform amplification effect is due to the nonuniformity of the hyperbolic field. The results indicate that the location error is highly dependent on the position of the true source relative to the sensor network layout.


## 3. The Mechanism of the Sensor Network in the Control of the Source Location

### 3.1. Geometrical Spreading Effect of the Sensor Network

Since the source location accuracy in the monitoring area is due to the nonuniformity of the associated hyperbolic field, as seen in Figure2, the line segment between the two sensors (or the foci) can be regarded as the original line where all hyperbolas are evenly spaced. But, the interval between the adjacent hyperbolas becomes larger and larger along these hyperbolas spread into the space. This effect associated with the nonuniform hyperbolic field is defined as the geometrical spreading effect. The density of hyperbolas in the quickest change direction is used to study the geometrical spreading effect quantitatively. Based on the geometry and mathematics theory, it is known that hyperbolas and ellipses are orthogonal when they are confocal. So, the quickest change direction of the hyperbola is the direction along its confocal ellipse.Based on the space geometry theory, the geometrical spreading quantitative formulas can be derived as follows.As shown in Figure3, sensors T1 and T2 are the foci of the hyperbola and ellipse. A1A2 is the real axis of the hyperbola with the length of 2a, and B1B2 is the imaginary axis which has the length of 2b. E1E2 is the long axis of the ellipse which has the length of 2e, and F1F2 is the minor axis with the length of 2f. The hyperbola and ellipse in this confocal system can be expressed by the following equations:(5)x2a2−y2b2=1,a2+b2=c2,(6)x2e2+y2f2=1,c2+f2=e2.Figure 3
The confocal ellipse and hyperbola with fociT1 and T2.As shown in Figure4, hyperboloid of two sheets and ellipsoid are generated by the rotation, about the X-axis, of corresponding hyperbola and ellipse located in the XY plane. The equations of the hyperboloid and ellipsoid are as follows:(7)x2a2−y2+z2b2=1,(8)x2e2+y2+z2f2=1.Figure 4
The confocal ellipsoid and hyperboloid hyperbola with fociT1 and T2.The trajectory of the intersecting points of the confocal hyperboloid and ellipsoid is a curve, and its equation can be expressed as follows:(9)y2+z2=b2f2c2,x=±aec.From mathematics, the change rate of the intersection curve (9) with the parameter a can be expressed as(10)drda=dxda2+dyda2+dzda2.Actually,dr/da is also the change rate of the hyperboloid along its confocal ellipsoid, which means the change rate of the source location accuracy. The parameter a is the vertex of the hyperboloid, which determines the position of the hyperboloid on the X-axis. Also, the term of dr/da can be interpreted as the distance that the hyperboloid will move along on the ellipsoid relative to the change of its position on the X-axis.Differentiating equation (9) with respect to a and substituting the partial derivatives into equation (10), we obtain the following equation:(11)drda=1ce2+a2f4y2+a2f4z2.Therefore, the density of hyperboloids at any point in the three-dimensional space can be expressed as(12)ρs=ce2+a2f4/y2+a2f4/z2.Equations (11) and (12) are the geometrical spreading quantitative formulas in the three-dimensional space. Equation (11) will be degraded as follows in the 2-dimensional plane:(13)drda=1+f2b2.According to the properties of the hyperbola, the larger value ofa, the further the hyperbola from the central line. So, the term of dr/da can be interpreted as the distance that the hyperboloid will move along on the ellipsoid relative to the change of its position on the X-axis.Equation (14) is given based on equation (5) and normalized by parameter c:(14)1b2=11−a2.The relationship between1/b2 and a is shown in Figure 5. As shown in Figure 5, 1/b2 increases with a, and this increase rate is accelerated when a approaches c (c = 1). So, the term 1/b2 represents the influence of the position of hyperbolas. Furthermore, the influence increases dramatically along with the hyperbolas located close to their foci. According to the discussion above, it is understood that the location errors on the hyperbolas (except the intersection point of hyperbolas and X-axis) will increase tremendously with the hyperbolas located close to their foci.Figure 5
The relationship between1/b2 and a (c = 1).As shown in Figure3, point K is the intersection of the confocal hyperbola and ellipse. Based on the geometry theory, the distance from point K to the origin point O can be calculated as follows:(15)d=x2+y2=a2+f2.According to equation (15), the value of d only depends on f2 once parameter a is determined. The relationship between distance d and parameter f is shown in Figure 6. As illustrated in Figure 6, the value of d increases with the increase of the parameter f, and d increases almost linearly with a slope of 1 when f is greater than 1.Figure 6
The relationship betweend and f (c = 1).Based on the abovementioned analysis, the change of any point in the hyperbolic field is determined by the distance from the origin and the position of its associated hyperbola. The elliptic equation (6) will be transformed into a straight line when f = 0, and this line is the X-axis. The change rate of hyperbolas along the X-axis is a constant and is equal to 1.Besides the effect of the decreasing density as the sources move away from the center of the sensor network, what is more important is that the density of hyperbolas is uneven spreading. Generally speaking, the density of hyperbolas in the center of the sensor network is higher than the density of those at the edge or outside of network. Thus, the location error of the sources in the center is much smaller than that of the sources at the edge or outside of the network. Based on the characteristic of the uneven geometrical spreading, the location error increases with a nonlinear trend as the sources move away from the center of the network. Also, when the sources are outside of the network, the location error increases sharply because the density of hyperbolas decreases sharply. The geometrical spreading effect shows that it is impossible to get an accurate location when a source is outside of the network. The uneven spreading effect reveals the importance of the relative position of a source and a sensor network.


### 3.2. Directional Control Effect of the Sensor Network

MS/AE monitoring tests in the field and laboratory both show that the source location errors in different directions are also different. For example, the location error in theZ-direction (vertical direction) is usually much larger than the errors in the X- and Y-direction (horizontal plane) during the MS monitoring process in underground mines. So, how the sensor network layout affects the source location in different directions? In the following part, the influence mechanism of the sensor network layout on source location in different directions is studied mathematically, which will explain why the accuracy of location in different directions is different.For the hyperboloid equation (7), based on the derivative rule of the implicit function, the change rate of any point in the X-direction due to the change of the position of a certain hyperboloid can be obtained as(16)2xa2dx/da−2ax2a4−y2+z2−1−2ac2−a22=0.The density of hyperboloids in theX-direction can be calculated as(17)ρsx=dadx=xaa4y2+z2/c2−a22+x2.Similarly, the density of hyperboloids in theY- and Z-direction can be expressed as follows:(18)ρsy=dady=−yc2−a2x2c2−a22/a3+ay2+z2,(19)ρsz=dadz=−zc2−a2x2+z2c2−a22/a3+y2a.Equations (19)–(21) are the quantitative solution formulae of the directional control effect. The equations of the density of hyperboloids in the Y- and Z-direction have similar expressions, but the equation of the density of hyperboloid in the X-direction is totally different. This is because there are two sensors in the X-direction which forms a pair of sensors, while there are no paired sensors on the Y- and Z-direction (Figure 4). The source location accuracy in the X-direction is much higher than that in the Y- and Z-direction.For a two-dimensional source location, it is very difficult to have a good location accuracy in theY-direction when all sensors are located on the X-axis. As shown in Figure 7, three sensors are deployed on the X-axis. According to the geometrical spreading effect, the location accuracy is closely related to the mesh size and distribution in the hyperbolic field. The smaller the mesh size is, the higher the location accuracy is. Based on the directional control effect, the source location accuracy in different directions can be estimated by the shape and the size of the mesh in the hyperbolic field. For the area between the sensors T1 and T3, except for the area I and II, most of the meshes are rhomboidal in shape with the long diagonal line almost parallel to the Y-axis, which clearly shows that the location accuracy in the X-direction is much higher than the accuracy in the Y-direction. In addition, the ratio of the length and the width of these meshes becomes larger and larger as the source moves away from the X-axis, which indicates that the location accuracy in the Y-direction decreases faster than that in the X-direction.Figure 7
Hyperbolic field associated with three sensors located in theX-axis.The direction control effect clearly shows that a source location with good accuracy in all directions can be obtained only when sensors are deployed on all directions. So, it is clear that a sensor network should not be “dimension-reduced.” For example, compared to the 3-dimensional underground mining space, most sensors in underground mines are deployed on a near horizontal plane, and these networks are vertical-dimension-reduced (Z-direction-reduced). As a result, the source location accuracy control in the direction (Z-direction) perpendicular to the horizontal plane will be very poor. This is why the vertical location errors are always larger than the horizontal location errors for MS monitoring projects in mines. Therefore, the sensors should be arranged in all directions (X-direction, Y-direction, and Z-direction) if the monitoring object is a 3-dimensional space structure.


### 3.3. Propagation Characteristics of the Velocity Error and Arrival Time Error

It is generally known that the AE/MS source location accuracy is mainly affected by the sensors coordinates and network layouts, wave velocity model, arrival time, and source location methods. Among the abovementioned influence factors, sensors coordinates, velocity model, and arrival time are the input parameters. According to the abovementioned discussion, the sensor network layout itself does not induce any errors during the source location process, and it merely transfers the input data errors (such as sensor coordinate, velocity, and arrival time errors) already existing. With the intelligent instruments, the accuracy of sensor coordinates can be as high as 0.2×10−6 mm. Thus, the error of coordinates is normally negligible. But, the velocity and arrival time errors are inevitable. Therefore, it is important to study how these velocity and arrival time errors transfer to the final location errors.According to the hyperbolic governing equation of the source location, the shape of the hyperbola of the source location just depends on the product of wave velocity and the difference of the arrival time. For convenience, this product is denoted as a signal element, namely,(20)g=2aij=vtj−ti.If there is total initial errorδg, the source will be wrongly located on the adjacent hyperbola with focal length d+δg. If the sensor coordinate error is negligible, then, δg is only affected by the arrival time error and velocity error. If there is only arrival time error δt, according to equation (20), the total initial error is(21)δgt=vδt.Equation (21) clearly demonstrates that δt is always amplified v times, where v is the wave velocity. Also, if the velocity is a constant, the total errors is not affected by the position of a source when the arrival time error is a primary source of initial error.If there is only wave velocity errorδv, it is clear that the total initial error is as follows:(22)δgv=δvdt,where dt is the arrival-time difference of two sensors and dt=tj−ti.It is clear from equation (22) that the total error affected by the velocity error is different from the arrival time error. The velocity error is amplified by dt. Because the arrival-time difference dt is highly affected by the hyperbola where the source is, the total initial error depends on the position of the source and only velocity errors are under consideration.The arrival-time difference is denoted as the “critical arrival-time difference (dtc)” when the total error caused by the arrival time error equals to total error caused by the velocity error, i.e., δgv=δgt. The hyperbola calculated by the critical arrival-time difference is called the critical hyperbola. As illustrated in Figure 8, the monitoring area is divided into two parts by the critical hyperbola. Based on equation (21), the location errors are mainly caused by the arrival time error in part I, but the velocity error has the major influence on source location accuracy in part II. The ratio of the area II and area I is denoted as error coefficient k. For a square area with a length of 1000 meters, the ratio of area II and I is about 12, which means that the source location error is mainly controlled by the velocity error in most monitoring areas.Figure 8
The areas controlled by the arrival time error and wave velocity error.The critical hyperbola divides the monitoring area into two regions. Generally speaking, the arrival time error has the major influence on source location in the central area. But, the wave velocity error becomes the leading factor that affects the source location error as the sources move away from the central area. Based on the concept of the critical arrival-time difference, the error coefficientk is getting smaller and smaller with the reduction of the distance between two sensors. For example, k will drop to 9, 5, 2.5, and 1.5 if the distance of two sensors reduces to 800, 500, 300, and 200 meters. So, the source location errors of most monitoring areas are mainly affected by the velocity errors in field MS monitor project because the distance of any two sensors is usually greater than 200 m. In contrast, the AE source location errors are mainly controlled by the arrival time error in the laboratory rock AE location tests.


## 3.1. Geometrical Spreading Effect of the Sensor Network

Since the source location accuracy in the monitoring area is due to the nonuniformity of the associated hyperbolic field, as seen in Figure2, the line segment between the two sensors (or the foci) can be regarded as the original line where all hyperbolas are evenly spaced. But, the interval between the adjacent hyperbolas becomes larger and larger along these hyperbolas spread into the space. This effect associated with the nonuniform hyperbolic field is defined as the geometrical spreading effect. The density of hyperbolas in the quickest change direction is used to study the geometrical spreading effect quantitatively. Based on the geometry and mathematics theory, it is known that hyperbolas and ellipses are orthogonal when they are confocal. So, the quickest change direction of the hyperbola is the direction along its confocal ellipse.Based on the space geometry theory, the geometrical spreading quantitative formulas can be derived as follows.As shown in Figure3, sensors T1 and T2 are the foci of the hyperbola and ellipse. A1A2 is the real axis of the hyperbola with the length of 2a, and B1B2 is the imaginary axis which has the length of 2b. E1E2 is the long axis of the ellipse which has the length of 2e, and F1F2 is the minor axis with the length of 2f. The hyperbola and ellipse in this confocal system can be expressed by the following equations:(5)x2a2−y2b2=1,a2+b2=c2,(6)x2e2+y2f2=1,c2+f2=e2.Figure 3
The confocal ellipse and hyperbola with fociT1 and T2.As shown in Figure4, hyperboloid of two sheets and ellipsoid are generated by the rotation, about the X-axis, of corresponding hyperbola and ellipse located in the XY plane. The equations of the hyperboloid and ellipsoid are as follows:(7)x2a2−y2+z2b2=1,(8)x2e2+y2+z2f2=1.Figure 4
The confocal ellipsoid and hyperboloid hyperbola with fociT1 and T2.The trajectory of the intersecting points of the confocal hyperboloid and ellipsoid is a curve, and its equation can be expressed as follows:(9)y2+z2=b2f2c2,x=±aec.From mathematics, the change rate of the intersection curve (9) with the parameter a can be expressed as(10)drda=dxda2+dyda2+dzda2.Actually,dr/da is also the change rate of the hyperboloid along its confocal ellipsoid, which means the change rate of the source location accuracy. The parameter a is the vertex of the hyperboloid, which determines the position of the hyperboloid on the X-axis. Also, the term of dr/da can be interpreted as the distance that the hyperboloid will move along on the ellipsoid relative to the change of its position on the X-axis.Differentiating equation (9) with respect to a and substituting the partial derivatives into equation (10), we obtain the following equation:(11)drda=1ce2+a2f4y2+a2f4z2.Therefore, the density of hyperboloids at any point in the three-dimensional space can be expressed as(12)ρs=ce2+a2f4/y2+a2f4/z2.Equations (11) and (12) are the geometrical spreading quantitative formulas in the three-dimensional space. Equation (11) will be degraded as follows in the 2-dimensional plane:(13)drda=1+f2b2.According to the properties of the hyperbola, the larger value ofa, the further the hyperbola from the central line. So, the term of dr/da can be interpreted as the distance that the hyperboloid will move along on the ellipsoid relative to the change of its position on the X-axis.Equation (14) is given based on equation (5) and normalized by parameter c:(14)1b2=11−a2.The relationship between1/b2 and a is shown in Figure 5. As shown in Figure 5, 1/b2 increases with a, and this increase rate is accelerated when a approaches c (c = 1). So, the term 1/b2 represents the influence of the position of hyperbolas. Furthermore, the influence increases dramatically along with the hyperbolas located close to their foci. According to the discussion above, it is understood that the location errors on the hyperbolas (except the intersection point of hyperbolas and X-axis) will increase tremendously with the hyperbolas located close to their foci.Figure 5
The relationship between1/b2 and a (c = 1).As shown in Figure3, point K is the intersection of the confocal hyperbola and ellipse. Based on the geometry theory, the distance from point K to the origin point O can be calculated as follows:(15)d=x2+y2=a2+f2.According to equation (15), the value of d only depends on f2 once parameter a is determined. The relationship between distance d and parameter f is shown in Figure 6. As illustrated in Figure 6, the value of d increases with the increase of the parameter f, and d increases almost linearly with a slope of 1 when f is greater than 1.Figure 6
The relationship betweend and f (c = 1).Based on the abovementioned analysis, the change of any point in the hyperbolic field is determined by the distance from the origin and the position of its associated hyperbola. The elliptic equation (6) will be transformed into a straight line when f = 0, and this line is the X-axis. The change rate of hyperbolas along the X-axis is a constant and is equal to 1.Besides the effect of the decreasing density as the sources move away from the center of the sensor network, what is more important is that the density of hyperbolas is uneven spreading. Generally speaking, the density of hyperbolas in the center of the sensor network is higher than the density of those at the edge or outside of network. Thus, the location error of the sources in the center is much smaller than that of the sources at the edge or outside of the network. Based on the characteristic of the uneven geometrical spreading, the location error increases with a nonlinear trend as the sources move away from the center of the network. Also, when the sources are outside of the network, the location error increases sharply because the density of hyperbolas decreases sharply. The geometrical spreading effect shows that it is impossible to get an accurate location when a source is outside of the network. The uneven spreading effect reveals the importance of the relative position of a source and a sensor network.


## 3.2. Directional Control Effect of the Sensor Network

MS/AE monitoring tests in the field and laboratory both show that the source location errors in different directions are also different. For example, the location error in theZ-direction (vertical direction) is usually much larger than the errors in the X- and Y-direction (horizontal plane) during the MS monitoring process in underground mines. So, how the sensor network layout affects the source location in different directions? In the following part, the influence mechanism of the sensor network layout on source location in different directions is studied mathematically, which will explain why the accuracy of location in different directions is different.For the hyperboloid equation (7), based on the derivative rule of the implicit function, the change rate of any point in the X-direction due to the change of the position of a certain hyperboloid can be obtained as(16)2xa2dx/da−2ax2a4−y2+z2−1−2ac2−a22=0.The density of hyperboloids in theX-direction can be calculated as(17)ρsx=dadx=xaa4y2+z2/c2−a22+x2.Similarly, the density of hyperboloids in theY- and Z-direction can be expressed as follows:(18)ρsy=dady=−yc2−a2x2c2−a22/a3+ay2+z2,(19)ρsz=dadz=−zc2−a2x2+z2c2−a22/a3+y2a.Equations (19)–(21) are the quantitative solution formulae of the directional control effect. The equations of the density of hyperboloids in the Y- and Z-direction have similar expressions, but the equation of the density of hyperboloid in the X-direction is totally different. This is because there are two sensors in the X-direction which forms a pair of sensors, while there are no paired sensors on the Y- and Z-direction (Figure 4). The source location accuracy in the X-direction is much higher than that in the Y- and Z-direction.For a two-dimensional source location, it is very difficult to have a good location accuracy in theY-direction when all sensors are located on the X-axis. As shown in Figure 7, three sensors are deployed on the X-axis. According to the geometrical spreading effect, the location accuracy is closely related to the mesh size and distribution in the hyperbolic field. The smaller the mesh size is, the higher the location accuracy is. Based on the directional control effect, the source location accuracy in different directions can be estimated by the shape and the size of the mesh in the hyperbolic field. For the area between the sensors T1 and T3, except for the area I and II, most of the meshes are rhomboidal in shape with the long diagonal line almost parallel to the Y-axis, which clearly shows that the location accuracy in the X-direction is much higher than the accuracy in the Y-direction. In addition, the ratio of the length and the width of these meshes becomes larger and larger as the source moves away from the X-axis, which indicates that the location accuracy in the Y-direction decreases faster than that in the X-direction.Figure 7
Hyperbolic field associated with three sensors located in theX-axis.The direction control effect clearly shows that a source location with good accuracy in all directions can be obtained only when sensors are deployed on all directions. So, it is clear that a sensor network should not be “dimension-reduced.” For example, compared to the 3-dimensional underground mining space, most sensors in underground mines are deployed on a near horizontal plane, and these networks are vertical-dimension-reduced (Z-direction-reduced). As a result, the source location accuracy control in the direction (Z-direction) perpendicular to the horizontal plane will be very poor. This is why the vertical location errors are always larger than the horizontal location errors for MS monitoring projects in mines. Therefore, the sensors should be arranged in all directions (X-direction, Y-direction, and Z-direction) if the monitoring object is a 3-dimensional space structure.


## 3.3. Propagation Characteristics of the Velocity Error and Arrival Time Error

It is generally known that the AE/MS source location accuracy is mainly affected by the sensors coordinates and network layouts, wave velocity model, arrival time, and source location methods. Among the abovementioned influence factors, sensors coordinates, velocity model, and arrival time are the input parameters. According to the abovementioned discussion, the sensor network layout itself does not induce any errors during the source location process, and it merely transfers the input data errors (such as sensor coordinate, velocity, and arrival time errors) already existing. With the intelligent instruments, the accuracy of sensor coordinates can be as high as 0.2×10−6 mm. Thus, the error of coordinates is normally negligible. But, the velocity and arrival time errors are inevitable. Therefore, it is important to study how these velocity and arrival time errors transfer to the final location errors.According to the hyperbolic governing equation of the source location, the shape of the hyperbola of the source location just depends on the product of wave velocity and the difference of the arrival time. For convenience, this product is denoted as a signal element, namely,(20)g=2aij=vtj−ti.If there is total initial errorδg, the source will be wrongly located on the adjacent hyperbola with focal length d+δg. If the sensor coordinate error is negligible, then, δg is only affected by the arrival time error and velocity error. If there is only arrival time error δt, according to equation (20), the total initial error is(21)δgt=vδt.Equation (21) clearly demonstrates that δt is always amplified v times, where v is the wave velocity. Also, if the velocity is a constant, the total errors is not affected by the position of a source when the arrival time error is a primary source of initial error.If there is only wave velocity errorδv, it is clear that the total initial error is as follows:(22)δgv=δvdt,where dt is the arrival-time difference of two sensors and dt=tj−ti.It is clear from equation (22) that the total error affected by the velocity error is different from the arrival time error. The velocity error is amplified by dt. Because the arrival-time difference dt is highly affected by the hyperbola where the source is, the total initial error depends on the position of the source and only velocity errors are under consideration.The arrival-time difference is denoted as the “critical arrival-time difference (dtc)” when the total error caused by the arrival time error equals to total error caused by the velocity error, i.e., δgv=δgt. The hyperbola calculated by the critical arrival-time difference is called the critical hyperbola. As illustrated in Figure 8, the monitoring area is divided into two parts by the critical hyperbola. Based on equation (21), the location errors are mainly caused by the arrival time error in part I, but the velocity error has the major influence on source location accuracy in part II. The ratio of the area II and area I is denoted as error coefficient k. For a square area with a length of 1000 meters, the ratio of area II and I is about 12, which means that the source location error is mainly controlled by the velocity error in most monitoring areas.Figure 8
The areas controlled by the arrival time error and wave velocity error.The critical hyperbola divides the monitoring area into two regions. Generally speaking, the arrival time error has the major influence on source location in the central area. But, the wave velocity error becomes the leading factor that affects the source location error as the sources move away from the central area. Based on the concept of the critical arrival-time difference, the error coefficientk is getting smaller and smaller with the reduction of the distance between two sensors. For example, k will drop to 9, 5, 2.5, and 1.5 if the distance of two sensors reduces to 800, 500, 300, and 200 meters. So, the source location errors of most monitoring areas are mainly affected by the velocity errors in field MS monitor project because the distance of any two sensors is usually greater than 200 m. In contrast, the AE source location errors are mainly controlled by the arrival time error in the laboratory rock AE location tests.


## 4. Laboratory Testing

### 4.1. Experimental System and Specimen Preparation

The experimental system (Figure9) consists of the 32-channel AE workstation form Physical Acoustics Company (PAC) and the 5058PR high-voltage pulse generator from Olympus Company. The 5058PR high-voltage pulse generator is designed for the ultrasonic detection of high attenuation material such as rock and coal, which can generate high-voltage ultrasonic signals. 100 V pulse signals are generated by 5058PR as the known AE sources in this test. The R15 AE sensors from PAC are used to record the AE signals. The resonant frequency of R15 is 150 kHz, and its operating frequency is from 150 to 400 kHz. The sampling rate is set as 10 MHz.Figure 9
Sketch of the experimental system of the sensor network and AE source location test.The granite and marble with high homogeneity were used and cut into a plate and cuboid. The thin plate specimen (GP-1) is granite, and its size is 400 mm × 250 mm × 50 mm (Figure11). The cuboid specimen (MP-2) is marble, and its size is 250 mm × 200 mm × 250 mm (Figure 10). There are five drill holes in the specimen MP-2, which were used to generate internal AE sources during the experiment. The P-wave velocity of the specimen GP-1 and MP-2 is 4472.9 m/s and 3192.6 m/s.Figure 10
The AE sensor network layout for the specimen GP-1, (a) is the sketch map, (b) is the physical photo.
(a)(b)


### 4.2. Experimental Program

Two different AE sensor networks were designed by using the GP-1 and MP-2 specimen to study the impacts of the sensor network on the source location. These two different sensor networks are marked as GP-1 sensor network source location and MP-2 sensor network source location experimental program. In order to enhance the coupling effect between the AE sensor and the specimen, the AE sensors are fixed on the specimen by using the tailor-made AE sensor mounting device, and the special coupling agent for AE test is used as well. In addition, the coupling quality of AE sensors is tested by the pencil lead breaking before the formal experiments.


#### 4.2.1. GP-1 Sensor Network Source Location Experimental Program

As shown in Figure11, the AE sensor network consists of 12 AE sensors in GP-1 experimental program. There are 29 AE pulse sources on the thin plate specimen. The coordinates of the sensors and sources are given in Tables 1 and 2.Table 1
The coordinates of sensors on the specimen GP-1.
SensorsX (mm)Y (mm)Z (mm)T110024928T229024828T3104030T4298027T510012050T61513250T72602950T830012250T925120550T1015021050T11012927T1239812028Table 2
The coordinates of the pulse sources on the specimen GP-1.
SourcesX (mm)Y (mm)Z (mm)S120012050S217016050S315012050S417018050S52308050S625012050S723016050S820021050S912016050S101208050S112003050S122808050S1328016050S147018050S154012050S16706050S173306050S1836010050S1933018050S2060030S21200030S22340030S2334025030S2420025030S256025030S26018030S2706030S2840018030S294006030


#### 4.2.2. MP-2 Sensor Network Source Location Experimental Program

As illustrated in Figure10, the AE sensor network consists of 18 AE sensors in MP-2 experimental program. According to experimental requirements, 5 AE sources are located in the drilling holes, and 16 AE sources are located on the surface of the cuboid. The coordinates of the sensors and sources are given in Tables 3 and 4.Figure 11
The AE sensor network layout for the specimen MP-2: (a) the sketch map and (b) the physical photo.
(a)(b)Table 3
The coordinates of sensors on the specimen MP-2.
SensorsX (mm)Y (mm)Z (mm)T13030250T231168250T322528250T4219165250T51620160T632028T7220018T824830220T9248112111T102483028T11221202220T12120202218T1321820230T142920120T150180130T170102190T17020132T1803030Table 4
The coordinates of the pulse sources on the specimen MP-2.
SourcesX (mm)Y (mm)Z (mm)S1120100170S25544190S31708770S4190130180S580160110S6100160250S740100250S816060250S9180140250S1019020080S11140200130S1230200160S13250130200S14250180120S1525050120S162100210S17120070S18200210S190170220S20030220S21015030


### 4.3. Analysis of the AE Source Location

The simplex method and L1-norm were used to calculate the AE source in the test [9]. The 29 AE source locations in GP-1 experimental program and 21 source locations in MP-2 experimental program were calculated. The location errors of each AE source in the GP-1 and MP-2 sensor network are shown in Figure 12.Figure 12
AE source location errors of the sample GP-1 and MP-2.As clearly shown in Figure12, the location errors are, in general, the smallest when the sources are located at the center of the sensor network. The location accuracy decreases with a nonlinear downtrend as the source moves away from the center of the network. For the test GP-1, the seven sources at the center of the network have the highest horizontal location accuracy, and the horizontal location errors are less than 8 mm. However, the location error increases to around 20 mm when the source is outside the hexagon of sensors on the upper surface of sample GP-1. The sources in the five boreholes of test MP-2 have the highest location accuracy, and their errors are only 2∼4 mm. But, the location error of the sources at the surface of the sample is up to 30 mm.According to the directional control effect, it is a dimension-reduced (Z-direction) sensor network in test GP-1, and as a result, it is very hard to have an accurate location in the Z-direction because all the sensors are arranged nearly horizontally. The source location results in Figure 12 clearly show that the location error of each source, whether inside or outside the network, in the Z-direction is much larger than the location error in the horizontal direction (X-direction or Y-direction) in the GP-1 test. Thus, the total location error is mainly affected by the Z-direction error. The horizontal location errors (X-direction and Y-direction location errors) are further studied to verify the geometrical spreading effect. The contour map of location errors on the X-direction and Y-direction is plotted in Figure 13. The solid blue circles represent sensors. The deeper the red color is, the higher the location accuracy is. The location accuracy rate can be represented by the density of the contour lines. Also, the higher the contour density, the faster the location accuracy decreases.Figure 13
The contour map of AE source location errors on theX- and Y-direction of the sample GP-1. (a) X-direction. (b) Y-direction.
(a)(b)As illustrated in Figure13, the density of the location error contour varies at different points. Usually, the density of the location error contour in the center area is smaller than that at the edge or outside of the sensor network, which means that the decline rate of the location error is getting faster and faster as the source moves away from the center of the sensor network. In other words, the location error increases with a nonlinear growth trend as the source moves away from the center of the network.The distribution law of the location error in theX-direction and Y-direction is different from each other. As shown in Figure 13, the location error in the X-direction is mainly affected by two subnetworks: a rectangular subnetwork composed of T1∼T4 sensors and a orthohexagonal subnetwork composed of T5∼T10 sensors. Also, the location error in the Y-direction is mainly affected by three subnetworks, namely, a orthohexagonal subnetwork composed of T5∼T10 sensors and two triangular subnetworks composed of T6, T10, T11 and T7, T9, T12. The location errors in both the X-direction and Y-direction of the sources located in the corners of the rectangle (black circle areas) are very large and out of control because these four corners are outside of the sensor network. So, we will get a more optimized network configuration if the sensors T1, T2, T3, and T4 need move to the four corners of the rectangle.As can be seen from Tables1 and 2, the distances from the edge to the center of the sensor network are about 80 cm and 50 cm in the X-direction and Y-direction, respectively, in the GP-1 test. As can be seen from Tables 3 and 4, the distances from the edge to the center of the sensor network are about 90 cm in all directions in the MP-2 test. The relationship between the location error in different directions and the distance from the source to the center of the network is shown in Figure 14. When the source is inside the sensor network, the location error in each direction increases slowly with the increase in distance from the center of the network. But, when the source is outside the sensor network, the location error in every direction increases sharply with the increase in distance. There is an inflection point that the increasing rate of the location error increases significantly when the distance exceeds this point. As shown in Figure 14, the inflection points are 50 cm and 80 cm in the X-direction and Y-direction, respectively, in the GP-1 test and are 90 cm in all directions in the MP-2 test. According to Tables 1∼4, this infection point is just located on the edge of the sensor network.Figure 14
Relationship between the location errors and the distance from the center of the network. (a) GP-1. (b) MP-2.
(a)(b)The wave velocity error is an important type of input data error during the source location process. Here, we assume different errors in the velocity model and relocate the sources by using velocities with different errors. The calculated results are illustrated in Figure15. The location errors increase with the increase in the velocity error. For the sources inside and outside the network, the location errors caused by the same velocity error are completely different. For example, if 10% error is involved in the velocity, the location error is about 20 mm when the sources are inside the network, but the location error will reach 300 mm when the sources are outside the network. In practice, the velocity error is usually within 10%, which means that the location errors caused by the velocity error are acceptable if the sources are located in the sensor network.Figure 15
Relationship between location errors and velocity errors. (a) Inside of the sensor network. (b) Outside of the sensor network.
(a)(b)


## 4.1. Experimental System and Specimen Preparation

The experimental system (Figure9) consists of the 32-channel AE workstation form Physical Acoustics Company (PAC) and the 5058PR high-voltage pulse generator from Olympus Company. The 5058PR high-voltage pulse generator is designed for the ultrasonic detection of high attenuation material such as rock and coal, which can generate high-voltage ultrasonic signals. 100 V pulse signals are generated by 5058PR as the known AE sources in this test. The R15 AE sensors from PAC are used to record the AE signals. The resonant frequency of R15 is 150 kHz, and its operating frequency is from 150 to 400 kHz. The sampling rate is set as 10 MHz.Figure 9
Sketch of the experimental system of the sensor network and AE source location test.The granite and marble with high homogeneity were used and cut into a plate and cuboid. The thin plate specimen (GP-1) is granite, and its size is 400 mm × 250 mm × 50 mm (Figure11). The cuboid specimen (MP-2) is marble, and its size is 250 mm × 200 mm × 250 mm (Figure 10). There are five drill holes in the specimen MP-2, which were used to generate internal AE sources during the experiment. The P-wave velocity of the specimen GP-1 and MP-2 is 4472.9 m/s and 3192.6 m/s.Figure 10
The AE sensor network layout for the specimen GP-1, (a) is the sketch map, (b) is the physical photo.
(a)(b)


## 4.2. Experimental Program

Two different AE sensor networks were designed by using the GP-1 and MP-2 specimen to study the impacts of the sensor network on the source location. These two different sensor networks are marked as GP-1 sensor network source location and MP-2 sensor network source location experimental program. In order to enhance the coupling effect between the AE sensor and the specimen, the AE sensors are fixed on the specimen by using the tailor-made AE sensor mounting device, and the special coupling agent for AE test is used as well. In addition, the coupling quality of AE sensors is tested by the pencil lead breaking before the formal experiments.


### 4.2.1. GP-1 Sensor Network Source Location Experimental Program

As shown in Figure11, the AE sensor network consists of 12 AE sensors in GP-1 experimental program. There are 29 AE pulse sources on the thin plate specimen. The coordinates of the sensors and sources are given in Tables 1 and 2.Table 1
The coordinates of sensors on the specimen GP-1.
SensorsX (mm)Y (mm)Z (mm)T110024928T229024828T3104030T4298027T510012050T61513250T72602950T830012250T925120550T1015021050T11012927T1239812028Table 2
The coordinates of the pulse sources on the specimen GP-1.
SourcesX (mm)Y (mm)Z (mm)S120012050S217016050S315012050S417018050S52308050S625012050S723016050S820021050S912016050S101208050S112003050S122808050S1328016050S147018050S154012050S16706050S173306050S1836010050S1933018050S2060030S21200030S22340030S2334025030S2420025030S256025030S26018030S2706030S2840018030S294006030


### 4.2.2. MP-2 Sensor Network Source Location Experimental Program

As illustrated in Figure10, the AE sensor network consists of 18 AE sensors in MP-2 experimental program. According to experimental requirements, 5 AE sources are located in the drilling holes, and 16 AE sources are located on the surface of the cuboid. The coordinates of the sensors and sources are given in Tables 3 and 4.Figure 11
The AE sensor network layout for the specimen MP-2: (a) the sketch map and (b) the physical photo.
(a)(b)Table 3
The coordinates of sensors on the specimen MP-2.
SensorsX (mm)Y (mm)Z (mm)T13030250T231168250T322528250T4219165250T51620160T632028T7220018T824830220T9248112111T102483028T11221202220T12120202218T1321820230T142920120T150180130T170102190T17020132T1803030Table 4
The coordinates of the pulse sources on the specimen MP-2.
SourcesX (mm)Y (mm)Z (mm)S1120100170S25544190S31708770S4190130180S580160110S6100160250S740100250S816060250S9180140250S1019020080S11140200130S1230200160S13250130200S14250180120S1525050120S162100210S17120070S18200210S190170220S20030220S21015030


## 4.2.1. GP-1 Sensor Network Source Location Experimental Program

As shown in Figure11, the AE sensor network consists of 12 AE sensors in GP-1 experimental program. There are 29 AE pulse sources on the thin plate specimen. The coordinates of the sensors and sources are given in Tables 1 and 2.Table 1
The coordinates of sensors on the specimen GP-1.
SensorsX (mm)Y (mm)Z (mm)T110024928T229024828T3104030T4298027T510012050T61513250T72602950T830012250T925120550T1015021050T11012927T1239812028Table 2
The coordinates of the pulse sources on the specimen GP-1.
SourcesX (mm)Y (mm)Z (mm)S120012050S217016050S315012050S417018050S52308050S625012050S723016050S820021050S912016050S101208050S112003050S122808050S1328016050S147018050S154012050S16706050S173306050S1836010050S1933018050S2060030S21200030S22340030S2334025030S2420025030S256025030S26018030S2706030S2840018030S294006030


## 4.2.2. MP-2 Sensor Network Source Location Experimental Program

As illustrated in Figure10, the AE sensor network consists of 18 AE sensors in MP-2 experimental program. According to experimental requirements, 5 AE sources are located in the drilling holes, and 16 AE sources are located on the surface of the cuboid. The coordinates of the sensors and sources are given in Tables 3 and 4.Figure 11
The AE sensor network layout for the specimen MP-2: (a) the sketch map and (b) the physical photo.
(a)(b)Table 3
The coordinates of sensors on the specimen MP-2.
SensorsX (mm)Y (mm)Z (mm)T13030250T231168250T322528250T4219165250T51620160T632028T7220018T824830220T9248112111T102483028T11221202220T12120202218T1321820230T142920120T150180130T170102190T17020132T1803030Table 4
The coordinates of the pulse sources on the specimen MP-2.
SourcesX (mm)Y (mm)Z (mm)S1120100170S25544190S31708770S4190130180S580160110S6100160250S740100250S816060250S9180140250S1019020080S11140200130S1230200160S13250130200S14250180120S1525050120S162100210S17120070S18200210S190170220S20030220S21015030


## 4.3. Analysis of the AE Source Location

The simplex method and L1-norm were used to calculate the AE source in the test [9]. The 29 AE source locations in GP-1 experimental program and 21 source locations in MP-2 experimental program were calculated. The location errors of each AE source in the GP-1 and MP-2 sensor network are shown in Figure 12.Figure 12
AE source location errors of the sample GP-1 and MP-2.As clearly shown in Figure12, the location errors are, in general, the smallest when the sources are located at the center of the sensor network. The location accuracy decreases with a nonlinear downtrend as the source moves away from the center of the network. For the test GP-1, the seven sources at the center of the network have the highest horizontal location accuracy, and the horizontal location errors are less than 8 mm. However, the location error increases to around 20 mm when the source is outside the hexagon of sensors on the upper surface of sample GP-1. The sources in the five boreholes of test MP-2 have the highest location accuracy, and their errors are only 2∼4 mm. But, the location error of the sources at the surface of the sample is up to 30 mm.According to the directional control effect, it is a dimension-reduced (Z-direction) sensor network in test GP-1, and as a result, it is very hard to have an accurate location in the Z-direction because all the sensors are arranged nearly horizontally. The source location results in Figure 12 clearly show that the location error of each source, whether inside or outside the network, in the Z-direction is much larger than the location error in the horizontal direction (X-direction or Y-direction) in the GP-1 test. Thus, the total location error is mainly affected by the Z-direction error. The horizontal location errors (X-direction and Y-direction location errors) are further studied to verify the geometrical spreading effect. The contour map of location errors on the X-direction and Y-direction is plotted in Figure 13. The solid blue circles represent sensors. The deeper the red color is, the higher the location accuracy is. The location accuracy rate can be represented by the density of the contour lines. Also, the higher the contour density, the faster the location accuracy decreases.Figure 13
The contour map of AE source location errors on theX- and Y-direction of the sample GP-1. (a) X-direction. (b) Y-direction.
(a)(b)As illustrated in Figure13, the density of the location error contour varies at different points. Usually, the density of the location error contour in the center area is smaller than that at the edge or outside of the sensor network, which means that the decline rate of the location error is getting faster and faster as the source moves away from the center of the sensor network. In other words, the location error increases with a nonlinear growth trend as the source moves away from the center of the network.The distribution law of the location error in theX-direction and Y-direction is different from each other. As shown in Figure 13, the location error in the X-direction is mainly affected by two subnetworks: a rectangular subnetwork composed of T1∼T4 sensors and a orthohexagonal subnetwork composed of T5∼T10 sensors. Also, the location error in the Y-direction is mainly affected by three subnetworks, namely, a orthohexagonal subnetwork composed of T5∼T10 sensors and two triangular subnetworks composed of T6, T10, T11 and T7, T9, T12. The location errors in both the X-direction and Y-direction of the sources located in the corners of the rectangle (black circle areas) are very large and out of control because these four corners are outside of the sensor network. So, we will get a more optimized network configuration if the sensors T1, T2, T3, and T4 need move to the four corners of the rectangle.As can be seen from Tables1 and 2, the distances from the edge to the center of the sensor network are about 80 cm and 50 cm in the X-direction and Y-direction, respectively, in the GP-1 test. As can be seen from Tables 3 and 4, the distances from the edge to the center of the sensor network are about 90 cm in all directions in the MP-2 test. The relationship between the location error in different directions and the distance from the source to the center of the network is shown in Figure 14. When the source is inside the sensor network, the location error in each direction increases slowly with the increase in distance from the center of the network. But, when the source is outside the sensor network, the location error in every direction increases sharply with the increase in distance. There is an inflection point that the increasing rate of the location error increases significantly when the distance exceeds this point. As shown in Figure 14, the inflection points are 50 cm and 80 cm in the X-direction and Y-direction, respectively, in the GP-1 test and are 90 cm in all directions in the MP-2 test. According to Tables 1∼4, this infection point is just located on the edge of the sensor network.Figure 14
Relationship between the location errors and the distance from the center of the network. (a) GP-1. (b) MP-2.
(a)(b)The wave velocity error is an important type of input data error during the source location process. Here, we assume different errors in the velocity model and relocate the sources by using velocities with different errors. The calculated results are illustrated in Figure15. The location errors increase with the increase in the velocity error. For the sources inside and outside the network, the location errors caused by the same velocity error are completely different. For example, if 10% error is involved in the velocity, the location error is about 20 mm when the sources are inside the network, but the location error will reach 300 mm when the sources are outside the network. In practice, the velocity error is usually within 10%, which means that the location errors caused by the velocity error are acceptable if the sources are located in the sensor network.Figure 15
Relationship between location errors and velocity errors. (a) Inside of the sensor network. (b) Outside of the sensor network.
(a)(b)


## 5. Discussion

The MS monitoring area is a semiopen three-dimensional space in underground mines. The spatial form and scale and rock properties of the monitoring area are very complicated. In addition, the installation position of the sensor is also limited by the field conditions and coupled with the economic considerations, and the number of sensors in the monitoring system cannot be increased indefinitely. Therefore, it is of great practical significance to establish an optimal sensor network with the least number of sensors. Some principles of optimizing sensor network are given based on the mechanism of sensor network in the control of source location.As shown in Figure16, it is assumed that the main MS monitoring area is around −500 m level. Four different pairs of MS sensors are illustrated in Figure 16, and the vertical drop of every paired sensors is 200 miters. The two sensors in Figures 16(a), 16(c), and 16(d) are located on both sides of the −500 m level, and the vertical distance of the two sensors to the −500 m level is equal. However, the difference is that the angle between the connection line of the two sensors and the horizontal plane in the three cases is 90°, α and β, respectively. Also, the relationship of the three angles satisfies the following equation: α<β<90°. According to the directional control effect, the case of Figure 16(a) has the highest source location accuracy, the case of Figure 16(c) has the second highest, and the case of Figure 16(d) has the lowest accuracy in the vertical direction. Here, we define the angle between the connection line of the two sensors and the horizontal plane as the direction angle. Thus, when the sensors are located on both sides of the monitoring plane, the larger the direction angle is, the higher is the location accuracy in the direction perpendicular to the monitoring plane.Figure 16
Four typical paired sensors at −500 m monitoring level.
(a)(b)(c)(d)Although the direction angles of Figures16(a) and 16(b) are the same, the relative positions of sensors and monitoring level are different. The monitoring plane vertically bisects the connection line between the two sensors in Figure 16(a). But, in Figure 16(b), the monitoring plane vertically passes through one third of the connection line between the two sensors. According to the geometrical spreading effect, the location accuracy in the vertical direction of the case in Figure 16(b) is less than that in Figure 16(a).In addition, the whole monitoring area should be considered when designing the sensor network. The sensor network should surround the entire monitoring area. As illustrated in Figure17, the whole MS monitoring area is the rectangle surrounded by points A, B, C, and D. Although the direction angle w>y and x > z in Figure 17, the sensor network in Figure 17(a) does not surround the whole monitoring area. The shadow part is located outside the sensor network, so the location accuracy of the shadow part cannot be guaranteed. However, the sensor network completely surrounds the whole monitoring area in Figure 17(b), so the sensor network layout of Figure 17(b) is more preferable in practice.Figure 17
Optimization of the MS monitoring sensor network layout.According to the abovementioned analysis, combined with the geometrical spreading effect, directional control effect, and the critical hyperbola, some principles of the sensor network are given as follows:(1)
The sensor network shall surround the main monitoring area in theX-, Y-, and Z-direction. Also, the sensors should avoid installing on the same straight line, plane, hyperbola, or hyperboloid at the same time.(2)
In consideration of the real monitoring environmental and economic factors, the number of sensors should be increased as much as possible. In addition, we make the sensor network form an irregular spatial layout because this network layout can increase the arrival-time difference of different sensors, thereby reducing the effect of the arrival time error on the source location and improving the source location accuracy.(3)
After the monitoring area and the total number of sensors are determined, the sensors located at the center and edge of the monitoring area should be balanced. In this way, not only the overall location accuracy of the entire monitoring area can be balanced but also the influence of the wave velocity error of the sources located on the edge area of the network and the influence of the arrival time error of the sources located in the center area of the network can be reduced at the same time.(4)
Some sensor subnetworks should be formed based on the spatial characteristics of the monitoring area.(5)
Since the establishment of MS monitoring system is a costly and time-consuming project, it is necessary to consider not only the current requirements but also the long-time monitoring requirements for the MS monitoring when the sensor network is designed.


## 6. Conclusions

(1)
Based on the arrival-time-difference principle, the hyperbolic/hyperboloidal governing equations of the source location are established. The path of the potential source determined by two sensors is a hyperbola on the plane and a hyperboloid in the space. The position and the shape of a hyperbola or hyperboloid where the source is located just depend on the arrival-time difference when the velocity model is determined. The sensor network does not induce any errors during the source location process; it only affects the source location accuracy by amplifying the existing errors of the input data such as arrival time and wave velocity. Also, the amplification effect of the input data errors is characterized by nonuniformity because of the nonuniformity of the hyperbolic/hyperboloidal field. This nonuniformity amplification effect is an important reason that the accuracy and stability of source location depend on the relative position of the true source and the sensor network.(2)
According to the governing equation of the source location and the nonuniformity geometrical characteristics of the hyperbolic/hyperboloidal field, the geometrical spreading effect and the directional control effect of the sensor network are revealed, and the three-dimensional quantitative models of these two effects are established. In addition, these two effects are verified by two source location experiments with different AE sensor networks. The geometrical spreading effect indicates that the location accuracy declines when a source moves away from the center of the network. The decline rate of the source location accuracy is small when the source is inside the sensor network. But, the decline rate drops sharply when the source is crossing the edge of the sensor network. Also, this decline rate has nonuniformity characteristics. The directional control effect shows that the source location accuracy of the same source in different directions is different. As the source moves away from the center of the network, there exists a direction that the location accuracy in this direction decreases faster than the other directions. In other words, the geometrical spreading effect in different directions is also different.(3)
The influence of the wave velocity error and the arrival time error on the source location accuracy is analytically compared, and the propagation characteristics of these two input data errors during the source location process are revealed. The concepts of the critical arrival-time difference and the critical hyperbola/hyperboloid are proposed. The critical hyperbola/hyperboloid can divide the monitoring area into two different areas where the source location errors are dominated by the arrival time error and the wave velocity error, respectively. Usually, the arrival time errors have the major influence in the central region. However, as the source moves away from the center area of the sensor network, the velocity error gradually becomes the major influence in the remaining area. This study shows that the velocity error is more significant for the large-scale MS monitoring in field. In contrast, the source location error is mainly controlled by the arrival time error in the small-scale laboratory rock AE location test.(4)
The concept of the paired sensors and the direction angle of paired sensors are proposed, and the relationship between the source location and the layout of four typical paired sensors is discussed. When the two sensors are located on both sides of the MS monitoring plane, the larger the direction angle, the higher the source location accuracy in the direction perpendicular to the monitoring plane. An optimal sensor network should surround the potential monitoring area in the three directions ofX, Y, and Z, that is, the monitoring area is best located inside the sensor network. Also, a sensor network should not be dimension-reduced. According to the spatial distribution characteristics of the monitoring area, if the total number of sensors is determined, a more reasonable sensor network will be formed by appropriately reducing the number of sensors in the center area and adding sensors to the periphery of the monitoring area. In addition, the sensor subnetwork should be formed to increase the interference resistance performance of the sensor network and source position.



---

*Source: 2651214-2020-09-07.xml*