Artificial Seismic Source Field Research on the Impact of the Number and Layout of Stations on the Microseismic Location Error of Mines

A site experiment is performed herein within a 100m range using a high-frequency structure activity monitor to explore the impact of different factors on the microseismic source location and analyze the range of influence of the velocity model, number of stations, and array surface on the seismic source location. Moreover, the impact of wave velocity, velocity-free location algorithm, and position of the seismic source on the microseismic location error of mines is discussed by establishing the ideal theoretical model of the wave velocity location and with particle swarm optimization. -e impact of the number of stations and tables on the location precision is also explored by using the microseismic signals produced by the artificial seismic source. -e results show that, for the location model containing the velocity, the velocity error would greatly affect the location precision, and the velocityfree algorithm receives good location results.-e location result is more satisfactory when the seismic source point falls in between array envelope lines. -e seismic source location precision is in direct proportion to the number of stations. According to the experiment, within a 100m range, when the number of stations is over 12, the effect does not significantly grow with the increase of stations; the number of tables affects the location precision; and the multitable location effect is significantly superior to the single-table effect. -e research shows that the optimal station density is 0.0192%, and the appropriate sensor layout to form a multitable monitoring network may effectively enhance the microseismic source precision of mines through the selection of a velocity-free location model. On the contrary, the number of stations can be reduced on the premise of the allowable error of the seismic source location, which may effectively reduce the monitoring cost.


Introduction
e shallow coal resource in China has already been exhausted in recent years; hence, the shift to deep coal mining is inevitable.e average mining depth of mines in China has already reached 700 m, and the number of deep shafts would increase year by year.However, a more complicated and dangerous stress environment would lead to more dynamic disasters in coal mines.Mine earthquake is one of these primary disasters, which usually give rise to a series of secondary disasters and result in a tremendous loss of personnel and property.erefore, mine earthquake has become a challenge for both Chinese and foreign scholars and institutions [1,2].e microseismic monitoring system is widely applied in all coal mines [3][4][5][6][7].To date, the microseismic monitoring technology has made great progress in research and application in the world, for example, in analyzing focal mechanisms [8], simulating three-dimensional (3D) geological structures [9,10], estimating rock shear-wave velocity [11], and developing microseismic data processing software packages [12].However, a prominent unresolved problem in the seismic monitoring is the large location error.Consequently, studies on improving the microseismic location precision would be of great theoretical significance and practical value.
Most of the current methods of microseismic location have been developed from earthquake locations.Accordingly, many scholars conducted studies on the location algorithm and how to improve the location precision.Several research directions on the decrease of the location error, such as the location algorithm [13], velocity model, station layout [14], and the accuracy of a seismic instrument [15,16], have been presented.urber used the nonlinear Newton iteration method for location for the first time.Zhao and Zeng [17] then introduced the simplex method by relying on the rapid development of the computer.However, the simplex method has defects; for instance, it may easily fall in local error.Sabbione et al. [18] proposed an adaptive filtering method using the apex-shifted parabolic Radon transform to denoise downhole microseismic data.Dong [19] and Li et al. [20] proposed the velocity-free location method after considering the impact of the velocity error on the location precision.Afterward, Dong et al. [21] proposed the comprehensive location method based on a threedimensional analysis of the acoustic emission and microseismic source in an unknown velocity system.Zhou et al. [22] proposed the microseismic source location method based on the small-region signal pickup.Increasing the number of stations is considered an effective method of improving the location precision; hence, Gong et al. [23] established the principle of determining the optimal number of channels by aiming at how to determine a proper number of stations.Gong et al. [24] believed that six channels could guarantee the location precision of the seismic source and proposed the anisotropic location method caused by the consideration of the anisotropic error.Lurka and Swanson [25] constructed a model according to the inversion principle and evolutionary algorithm for a deep coal mine in Western United States.
is author then modified and improved the seismic location precision with a new target function model, and as a result, the seismic location precision was evidently enhanced compared with the constant velocity model.e station layout would affect the location precision.Kovaleva et al. [26] developed a comprehensive workflow to optimize the microseismic acquisition design and obtain optimum depth and location for deploying the borehole geophone arrays.Chen [27] studied the impact of a substation network on the location precision.After determining the optimal number of channels, Feng and Lan [28] conducted the optimal channel combiner fixed-point test to enhance the microseismic location precision.Jia and Li [29] proposed the reduction of the location error by optimizing the station layout.Ruigrok et al. [30] designed and applied a temporary array of 38 seismic stations.I. Sanina et al. [31] designed and applied a two-dimensional small-aperture seismic array.Meanwhile, Hu and Li [32] reduced the location error by accurately picking up the velocity of the P wave.Chen et al. [33] applied the ranking of the P wave to estimate the location error.
In conclusion, substantial research results have already been achieved for the microseismic source location.However, the abovementioned research mainly aims at the algorithm of a limited number of stations and station layout, and few comprehensive studies focused on the impact of comprehensive factors in a high-density station layout on the location error.In this study, an in situ microseismic monitoring experiment is conducted within a 100 m range for the first time with the high-frequency structure activity monitor to improve the microseismic location precision.Moreover, the impact of the location target function, velocity error, seismic source position, number of stations, and number of tables on the location precision is analyzed and studied by applying the theoretical analysis and field experimentation.

Microseismic Monitoring Experiment of an Artificial Seismic Source Site
Microseismic monitoring has become one of the effective methods of dynamic disaster control for all kinds of mines across the country, but location precision has always been a difficult problem to resolve.High-density station monitoring is considered an effective solution, but research on this topic is seldom performed in mines.Small-range microseismic monitoring is performed in this research.Moreover, a microseismic experiment site is built for the high-density station microseismic monitoring experiment, which can provide a theoretical basis for mine earthquake research.

Introduction to the Experimental Region.
A gold mine in Fuxin was selected for the small-range microseismic monitoring experiment.e geographic coordinates of the mine are as follows: east longitude, 121 °43′04″-121 °43′06″, and north latitude, 41 °53′03″-41 °53′04″.e surface soil of the coal mine is mild clay, and rocks in the mine mainly include biotite plagioclase mylonite and felsic mylonite, which are solid and of good integrity.e mines are explored through vertical and inclined shafts in a single cage with a balanced weight.Coal mining currently has four midlevels, namely, mid-223 m, mid-180 m, mid-140 m, and mid-100 m.Among which, the first (mid-223 m) has already been emptied and abandoned.e top-down horizontal sublevel mining method is adopted according to the status of mine production and the hoisting condition of the ore body.Underground mining is being adopted at the present, and vertical shaft-blind inclined shaft is employed for codevelopment based on the full utilization of the original development system.A vertical shaft is employed for transportation from the surface to 180 m below the ground, while an inclined shaft is adopted from 180 m to 140 m and 180 m to 100 m plane.e experiment section is the plane of 140 m, 180 m, and 100 m.

Monitoring Instrument.
Experiment monitoring includes sensors, cables, a 12-channel data hub, and a 48channel data recorder for constructing the physical deformation site (Figure 1).e sensors used in the experiment adopt a direction-free setting.Sensors feature high sensitivity; hence, the data acquisition frequency can be set as 100,000 Hz at the maximum.While recording the microblasting signal in the experiment, the sampling frequency is set to 25 kHz.

2
Advances in Civil Engineering 2.3.Array Layout.Sensors were mainly laid out in three planes and two inclined shafts.e vertical size of the distribution region was 80 m, while the horizontal size was between 100 m and 200 m.e array mainly consisted of 33 sensors.Table 1 shows their location information.Figures 2-4 depict the planar layout of the sensor.Figure 5 shows the space layout of the sensors.

Microblasting Seismic Source.
ree microblasting experiments were performed, and the blasting location was set in three operation planes.A monitoring station was set near the artificial seismic source to precisely collect the seismic moment.e arrival time difference of this station and the other monitoring stations was calculated and became the travel time of the seismic signal spent in reaching the detection station.Table 2 presents the artificial seismic source coordinates, amount of explosive, and the signal waveform recorded.Figure 5 illustrates the spatial arrangement of the artificial microseismic source.
e monitoring result showed that only the No. 2 blasting at the 140 m plane and the No. 3 blasting at the 100 m plane were well received.e No. 1 seismic signal at the 180 m plane and the effective acceptance amount of stations were quite dissatisfactory and thus abandoned.e analysis implied that this finding may be caused by the remote seismic location, the small amount of explosive, and the dense transportation channels, which affected the signal transmission from the artificial seismic source to all the sensors.Figure 6 shows the signal received by all stations after the second and third blasting.e picked up signals were filtered.Subsequently, onset time picking was then conducted, which was not elaborated here because of the length limit.
e equivalent velocity of the wave transmission in the medium in the theoretical model is simplified as V � 5600 m/s to make it convenient to explore the impact of different factors on the location error.e central point O of the monitoring region is (500, 500, 500).e location effect worked out for the seismic source within and without the region is studied.After the artificial seismic source coordinates and the velocity are given, the arrival time from the artificial seismic source to the sensors is worked out by supposing that the artificial seismic source is stirred at 0 moment and taken as the true value.e location calculation is then performed with the arrival time.e coordinates of seismic source points 1 to 11 are (500, 500, 500), (400, 400, 400), . . ., (−500, −500, −500).Artificial seismic source points 1-6 are within the envelope line of the monitoring region (i.e., AG line), while points 7-11 are beyond the envelope line (i.e., extending line of AG). e arrival time from the artificial seismic source points to all sensors is 0.1 ms as the minimum unit in the round-off method to simulate the error caused by the onset time picking.Figure 7 shows the model.

Location Algorithm.
Particle swarm optimization is an emerging evolutionary algorithm, which constantly collects the direction and speed of search according to the flight course and information transfer between swamps.e search process is mainly accomplished by relying on the interactions and the mutual influence of particles featuring easy implementation, fast convergence, and high precision.
e free flight of "particle swarm" in the solution space can perfectly resolve the problem of the final solution being the local optimal solution.e updated versions of the speed and the location of particle i are shown in the following equations: where w is the inertia weight; c 1 and c 2 are the integrals ranging between [2,4] known as the learning or accelerated factor; r 1 and r 2 are the random numbers between [0, 1] known as the maximum inertia weight and minimum inertia weight, respectively; d � 1, 2, 3, . . ., D, where D stands for Advances in Civil Engineering the dimension of the solution vector; V id , X id , P id , and P gd are the flight speed of particle i in the dth dimension space, location, optimal solution searched, and optimal solution searched by the swamp.e location was worked out through the MATLAB PSO toolbox with the following PSO parameter settings: learning factor: c 1 � c 2 � 2; inertia weight: r 1 � 0.9, r 2 � 0.4; number of particles: 20, w � 1; and end condition: maximum number of iterations � 20,000 or ε � 1e − 25.When the velocity is known, V � 5600 m/s and D � 3.Meanwhile, D � 4 when the location is worked out with the unknown velocity.

Impacts of the Location Methods on the Location.
is study located the two target functions of the known wave velocity and the unknown wave velocity through the ion swarm algorithm to explore the effects of different location methods on the location error.Equation ( 2) is the target function for solving the known wave velocity localization, and equation ( 3) is the target function for solving the unknown wave velocity localization; Figure 8 shows the location result error: Figure 8(a) illustrates that with the known wave velocity, the errors in the x, y, and z directions tended to increase with the source point away from the center of the monitoring station envelope.e error presented a jump change when the unknown wave velocity was used for the location, with the location error of each source point getting farther and farther away from the center point of the monitoring station.e location error was smaller in the area.e error outside the area also increased when the distance increased and became greater than the location error of the known wave velocity.Figure 8(b) depicts a large difference of the location error between the known and unknown wave velocities, which cannot re ect the actual geologic wave velocity.e in uences of the two location methods on the location error were not obvious when the wave velocity error was not considered.However, the variation trend of the location error was di erent.

In uences of the Wave Velocity Error on the Location.
e P-wave velocity of the rock samples in the experimental area was rst required to be determined by an ultrasonic velocity meter based on the known wave velocity model for the location.Di erences between the given wave velocity and the actual wave velocity were observed.erefore, the errors that are ±2%, ±5%, ±10%, and ±15% are applied to disturb the original wave velocity V 5600 m/s and take it as the true value of wave velocity.A disturbance error greater than 0 implied that the actual wave velocity was less than V; otherwise, the actual wave velocity was greater than V. e arrival time was calculated and taken as the initial value of the location calculation together with V 5600 m/s to verify the in uences of the wave velocity error on the location.Figure 9 shows the location distance error of each arti cial seismic source under di erent wave velocity error disturbances.
Figure 9 also depicts that the wave velocity error has a great in uence on the location performance when the known wave velocity was used for the location.e location error of each source point increased in line with the increase of the wave velocity error when the actual wave velocity was less than that used for the location.Meanwhile, the location error of Advances in Civil Engineering each source point inside the area increased in line with the increase of the wave velocity error when the actual wave velocity was greater than that used for the location.e location error of each source point outside the area was greater than that inside the area that approximately showed an exponential increase with the increase of the wave velocity error.e error brought by a lower location wave velocity was generally larger with a greater in uence on the location.

In uences of the Source Position on the Location.
Figure 10 clearly shows that the location error of each source point inside and outside the monitoring area was signi cantly di erent.
e location error increased in a straight line with the enlarged distance from the central point inside the area when the actual wave velocity was greater than location wave velocity, and the disturbed value was less than −5%.e change rate of the error suddenly became higher from source point 6 (i.e., the demarcation point between the inside and the outside) and then became lower.Its position of becoming smaller became farther and farther away with the decrease of the disturbed value.In addition, the change rate of the error eventually tended to be stable.e disturbed value was greater than 5% when the actual wave velocity was less than the location wave velocity.
e location error rst increased and then decreased with the enlarged distance from the central point inside the area and gradually increased again with a more enlarged distance from the demarcation point.e location error was very small when the disturbed value was ±2%. is error slowly increased in a straight line inside the area, but remained stable outside the area.e in uences of the source position on the location are suggested to be interfered by the wave velocity error.With the great error in the wave velocity, the location error changes could be divided into the three following areas from the central point of the area (Figure 10): near enlarged, mutation, and far enlarged areas.e change trend of the location error changed within the mutation area.Furthermore, the range and the changing characteristics of the mutation area depended on the wave velocity error.

Analysis on the Influences of the Quantity of Station Planes and Stations on the Location Error
e particle swarm algorithm was applied to undertake the source location to the microexplosion performed on the site.e location parameters are as follows: location area coordinates X [8800, 9000], Y [6100, 6400], and Z [95,185].e e ective arrival time was adopted for the location.
e ultrasonic testing indoor showed that the P-wave velocity of the rock samples of the experimental diggings was 6357 m/s.e unknown wave velocity was adopted for the microexplosion source location to avoid the in uences of the wave velocity error on the location performance.

In uences of the Quantity of Station Planes on the
Location.Di erent station sensors, including 180 m plane, 140 m plane, 100 m plane, inclined shaft (long and short inclined shafts), and some plane sensors, were successively applied for the location to analyze the high-density array for the location performance of microseismic monitoring.e unknown wave velocity location was adopted for a better location performance (Figure 11).When the single station plane was used for the location, the source location was more accurately realized only when the source was located on or nearby the monitoring plane.
e single location of the monitoring plane far away from the source would have an extremely great error, especially the vertical location error.When the sensor of the 140 m plane was adopted for the location of the 100 m plane, its vertical error unexpectedly 6 Advances in Civil Engineering reached 51 m, and the vertical distance between the two planes was only 40 m. e location performance was obviously greater than that of the single plane when two monitoring planes were applied for the location.However, this value was still less than that in the location with all the sensors used. is nding illustrates that the increase of the sensor "array plane" could largely reduce the location error with a positive and stimulating in uence on the microseismic location.

In uences of the Station Quantity on the Location.
Four sensors were selected from all the sensors in the array to verify the impact of the monitoring array density on the location performance and avoid locating in a single plane.e sensors at the other positions were successively superimposed to select di erent sensor quantities and perform an unknown wave velocity location for the microblasting on the 140 m plane.ree sets of data were selected to take the mean value for the location.Figure 12 shows the resulting location error.e left gure depicts that the sensor quantity plays a positive role in improving the location accuracy of the source.e location errors became increasingly smaller as the sensor quantity increased.e derivation was taken for the location error to obtain the variation trend curve of the error.As shown in the right side of Figure 12, the location error decreased with the increase of the sensor quantity when the sensor quantity was between 4 and 14. e value change of the derivative then became slow but was still less than 0, thereby indicating that the increase of the sensor quantity can improve the location performance, which proves that the high-density array is bene cial to control the microseismic monitoring accuracy.e source location error of the whole array was approximately 2.6 m.
e mathematical relationship between the location error ΔL and the sensor quantity n was obtained as shown in equation ( 4) based on the tting of the location error data in the gure.e correction determination coe cient Adj.R 2 0.98972, which shows that the tting result is very close to the real situation.
rough this equation, an exponential relationship can be found between both: ΔL 251.999 × e −(n/2.343)+ 3.166. (4) In the above calculation, the P-wave velocities of the 140 m plane and the 100 m plane obtained from the   Advances in Civil Engineering inversion of the whole array were 5898 m/s and 5903 m/s, respectively, which were rather di erent from the results of the indoor experiment.In addition, compared with the model analysis results, the location error of the actual monitoring was larger and related to the complexity of a geological structure.e location, size, and direction of various discontinuity surfaces, such as fractures, joints, and faults in the rock mass, all in uence the wave propagation in the medium.Inevitably, errors were found in the arrival time required for the location because the time was acquired from the waveform recorded in the noise environment through phase identi cation.On the one hand, the high-density array microseismic monitoring can extend the quantity of the arrival time data to reduce the signi cant impact of the arrival time error.On the other hand, the high-density station, especially when being three-dimensionally laid, decreases the travel time distance from the source to the various sensors and reduces the di erence between space and time and the theoretical calculation of the travel time data caused by the inhomogeneity of the rock mass.In summary, bene tting from the two abovementioned points, the high-density station can e ectively improve the location performance of the microseismic source of the mine.

Conclusions
e following conclusions are obtained from this study: (1) eoretical modeling was used herein to analyze the target function with the wave velocity, nonwave velocity, and source location performance.Consequently, the wave velocity error had a great in uence   8 Advances in Civil Engineering on the location result in the target function with the wave velocity, and such in uence exponentially increased.e positioning error was small when the source point was at the envelope center of the monitoring array.e error was larger when the source point was outside the envelope.
(2) e in uence of the station density on the location error was analyzed by a eld experiment.e experimental results indicated that the in uence of the station density can be divided into two parts: station quantity and station plane quantity.e location accuracy signi cantly increased with the increase of the station quantity when the station plane quantity was certain.When station quantity was more than 12, such an increase was no longer obvious.When the station quantity is certain, the increase of the station plane quantity can signi cantly improve the location accuracy.e accuracy of the multistation plane was higher than that of the single-station plane.
e best station density was 0.0192%.

Figure 8 :
Figure 8: (a) Location error distribution under the known wave velocity.(b) Location error distribution under the unknown wave velocity.

Figure 9 :
Figure 9: Location distance error distribution of di erent wave velocity errors.

Figure 10 :
Figure 10: Location distance error distribution of di erent positions.

Figure 11 :
Figure 11: Relationship between the location error and the sensor position.(a) 100 m plane seismic source location error.(b) 140 m plane seismic source location error.

Figure 12 :
Figure 12: Relationship between the location error and the number of stations.

Table 1 :
Sensor location information.