Method for improving monitoring capability of borehole-surface micro-seismic monitoring system

Abstract

A method for improving a monitoring capability of a borehole-surface micro-seismic monitoring system includes selecting multiple candidate points for installing surface wireless sensors to form a natural-number-coded candidate point set and combining a fixed number of candidate points randomly selected from the candidate point set with an underground installed sensor set to form a borehole-surface micro-seismic monitoring network; carrying out multiple random selections until a certain scale of borehole-surface micro-seismic monitoring network deployment plans are generated; establishing an evaluation model for a monitoring capability of each borehole-surface micro-seismic monitoring network deployment plan according to a propagation relation equation between a micro-seismic energy and a first-arrival peak amplitude of a P-wave, forming an initial population; determining an optimal borehole-surface micro-seismic monitoring network deployment plan through a genetic algorithm; and determining an optimal surface wireless sensor deployment plan that significantly improves the monitoring capability of the borehole-surface micro-seismic monitoring network deployment plan.

Claims

1. A method for improving a monitoring capability of a borehole-surface micro-seismic monitoring system comprising the following steps: (1) selecting multiple candidate points for installing surface wireless sensors to form a natural-number-coded candidate point set S={1,2,3,4,5, . . . , n}; (2) combining a fixed number m of candidate points randomly selected from the candidate point set S formed in step (1) with an underground installed sensor set U to form a borehole-surface micro-seismic monitoring network deployment plan G.sub.v=[S.sub.2.sup.1 S.sub.4.sup.2 . . . S.sub.n−2.sup.m U.sub.1 U.sub.2 . . . U.sub.k]; wherein, S.sub.2.sup.1 denotes a candidate point that is a first candidate point randomly selected from the candidate point set S and is a second candidate point in the candidate point set S; similarly, S.sub.n−2.sup.m denotes a candidate point that is an m-th candidate point randomly selected from the candidate point set S and is an (n−2)-th candidate point in the candidate point set S; and k denotes a number of underground sensors; (3) repeating step (2) until v=p borehole-surface micro-seismic monitoring network deployment plans are generated to form ap-scale deployment plan set G: $G = [\begin{matrix} G_{1} \\ G_{2} \\ .Math. \\ G_{p} \end{matrix}]$ (4) forming, by each borehole-surface micro-seismic monitoring network deployment plan G.sub.v generated in step (3), an initial population Gen: 401) determining, according to a micro-seismic signal acquired by the underground sensor, a propagation relation equation between a micro-seismic energy E and a first-arrival peak amplitude f of a P-wave: $f = E α_{1} \frac{1}{r} e^{- α_{2} r}$ wherein, α.sub.1 denotes an amplitude-energy ratio coefficient; α.sub.2 denotes an attenuation coefficient; and r denotes a distance from a micro-seismic source to the underground sensor; 402) forming a three-dimensional (3D) equidistant grid model comprising $floor (\frac{X \max - X \min}{d x} + 1) \times floor (\frac{Y \max - Y \min}{d y} + 1) \times floor (\frac{Z \max - Z \min}{d z} + 1)$ grids with an X-direction spacing dx, a Y-direction spacing dy, and a Z-direction spacing dz for a mining and production area defined by [Xmin, Xmax], [Ymin, Ymax], and [Zmin, Zmax], wherein m.sub.1, n.sub.1, and p.sub.1 denote a number of X-direction grids, a number of Y-direction grids, and a number of Z-direction grids, respectively; 403) rewriting the propagation relation equation determined in step 401) to obtain $E = \frac{f r}{α_{1} e^{- α_{2^{r}}}};$ and calculating a minimum micro-seismic energy E.sub.i,j,k.sup.min to trigger the borehole-surface micro-seismic monitoring network deployment plan G.sub.v to record the micro-seismic signal, at a point (X.sub.i, Y.sub.j, Z.sub.k) in the 3D equidistant grid model formed in step 402); wherein, i∈1,2, . . . , m1; j∈1,2, . . . , n1; k∈1,2, . . . , p1; v∈1,2, . . . , p; 404) establishing, according to step 403), an evaluation model for a monitoring capability Q.sub.v of the borehole-surface micro-seismic monitoring network deployment plan G.sub.v: $Q_{v} = \frac{{.Math.}_{i = 1}^{m_{1}} {.Math.}_{j = 1}^{n_{1}} {.Math.}_{k}^{p_{1}} E_{i, j, k}^{\min}}{m_{1} n_{1} p_{1}}$ 405) forming the initial population Gen: $Gen = [\begin{matrix} G_{1} & Q_{1} \\ G_{2} & Q_{2} \\ .Math. & .Math. \\ G_{p} & Q_{p} \end{matrix}]$ (5) determining, according to the initial population formed in step (4), an optimal borehole-surface micro-seismic monitoring network deployment plan through a genetic algorithm, and determining an optimal surface wireless sensor deployment plan.

2. The method for improving the monitoring capability of the borehole-surface micro-seismic monitoring system according to claim 1, wherein in step (1), the multiple candidate points selected for installing the surface wireless sensors satisfy the following conditions: the candidate points cooperate with the underground sensors to surround the mining and production area; the candidate points have a distance of no more than 2,000 m from the mining and production area; the candidate points avoid a waterlogged area, a surface water system, a highway facility, and a noisy place; and the candidate points provide strong fourth-generation/fifth-generation (4G/5G) network signals.

3. The method for improving the monitoring capability of the borehole-surface micro-seismic monitoring system according to claim 1, wherein step 401) comprises: determining α.sub.1 and α.sub.2 as follows: manually marking the first-arrival peak amplitude f of the P-wave recorded by each underground sensor; calculating a source position and micro-seismic energy E of multiple micro-seismic signals with different energy levels; calculating a distance r from the micro-seismic source to the underground sensor; and determining α.sub.1 and α.sub.2 by a nonlinear least squares (NLS) method.

4. The method for improving the monitoring capability of the borehole-surface micro-seismic monitoring system according to claim 2, wherein step 401) comprises: determining α.sub.1 and α.sub.2 as follows: manually marking the first-arrival peak amplitude f of the P-wave recorded by each underground sensor; calculating a source position and micro-seismic energy E of multiple micro-seismic signals with different energy levels; calculating a distance r from the micro-seismic source to the underground sensor; and determining α.sub.1 and α.sub.2 by a nonlinear least squares (NLS) method.

5. The method for improving the monitoring capability of the borehole-surface micro-seismic monitoring system according to claim 3, wherein in step 403), the calculating a minimum micro-seismic energy E.sub.i,j,k.sup.min to trigger the borehole-surface micro-seismic monitoring network deployment plan G.sub.v to record the micro-seismic signal, at a point (X.sub.i, Y.sub.j, Z.sub.k) comprises: 40301: determining, according to a micro-seismic positioning principle based on a first arrival time of the P-wave, that the micro-seismic monitoring system is triggered to record the micro-seismic signal when the first-arrival peak amplitude f of the P-wave received by at least four sensors is greater than or equal to three times an ambient noise level NL; 40302: calculating a distance r.sub.l from the point (X.sub.i,Y.sub.j, Z.sub.k) to each sensor in the borehole-surface micro-seismic monitoring network deployment plan G.sub.v; determining, according to step 40301, a first-arrival peak amplitude f.sub.l, of the P-wave required to trigger each sensor; and back-calculating, according to the propagation relation equation between the micro-seismic energy E and the first-arrival peak amplitude f of the P-wave determined in step 401, a micro-seismic energy E.sub.i,j,k.sup.l required to trigger each sensor: $E_{i, j, k}^{l} = \frac{f_{l} r_{l}}{α_{1} e^{- α_{2^{r} l}}},$ $wherein, l = 1, 2, .Math., m + k$ 40303: sorting, according to step 40301, the micro-seismic energy E.sub.i,j,k.sup.l calculated in step 40302 in an ascending order; and selecting a fourth micro-seismic energy after the sorting as the minimum micro-seismic energy E.sub.i,j,k.sup.min to trigger the borehole-surface micro-seismic monitoring network deployment plan to record the micro-seismic signal.

6. The method for improving the monitoring capability of the borehole-surface micro-seismic monitoring system according to claim 5, wherein in step 40301, the ambient noise level NL comprises a surface ambient noise level NL.sub.s monitored by the surface sensor installed on a surface and an underground ambient noise level NL.sub.u monitored by the underground sensor installed in an underground roadway.

7. The method for improving the monitoring capability of the borehole-surface micro-seismic monitoring system according to claim 6, wherein in step (5), the genetic algorithm sets a generation number of not less than 100; and the genetic algorithm carries out mutation operation through a mixture of adjacent gene mutation, gene insertion mutation, gene exchange mutation, three-point gene exchange mutation, and two-point inversion mutation, and carries out crossover operation through a mixture of partially mapped crossover, cycle crossover operator, edge recombination crossover, linear sequential crossover, ordered crossover operator, and uniform crossover.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0032] FIG. 1 is a flowchart of a method for improving a monitoring capability of a borehole-surface micro-seismic monitoring system according to the present disclosure;

[0033] FIG. 2 shows a comparison between original data and fitted data of a first-arrival peak amplitude f of a P-wave according to an embodiment;

[0034] FIG. 3 shows a three-dimensional (3D) grid model formed according to the embodiment; and

[0035] FIG. 4 shows a process of determining an optimal individual in an initial population according to an evaluation value of a monitoring capability based on a generation number according to the embodiment.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0036] The present disclosure is further described below with reference to the accompanying drawings.

[0037] As shown in FIG. 1, a method for improving a monitoring capability of a borehole-surface micro-seismic monitoring system includes the following steps: [0038] (1) Multiple candidate points for installing surface wireless sensors are selected to form a natural-number-coded candidate point set S={1,2,3,4,5, . . . , n}. [0039] (2) A fixed number, namely m, of candidate points randomly selected from the candidate point set S formed in step (1) are combined with underground installed sensor set U to form a borehole-surface micro-seismic monitoring network deployment plan G.sub.v=[S.sub.2.sup.1 S.sub.4.sup.2 . . . S.sub.n−2.sup.m U.sub.1 U.sub.2 . . . U.sub.k]. [0040] S.sub.2.sup.1 denotes a candidate point that is a first candidate point randomly selected from the candidate point set S and is a second candidate point in the candidate point set S; similarly, S.sub.n−2.sup.m denotes a candidate point that is an m-th candidate point randomly selected from the candidate point set S and is an (n−2)-th candidate point in the candidate point set S; and k denotes a number of underground sensors. [0041] (3) Step (2) is repeated until v=p borehole-surface micro-seismic monitoring network deployment plans are generated to form ap-scale deployment plan set G:

[00008] $G = [\begin{matrix} G_{1} \\ G_{2} \\ .Math. \\ G_{p} \end{matrix}]$ [0042] (4) Each borehole-surface micro-seismic monitoring network deployment plan G.sub.v generated in step (3) is combined to form an initial population Gen. [0043] 401) According to a micro-seismic signal acquired by the underground sensor, a propagation relation equation between micro-seismic energy E and first-arrival peak amplitude f of a P-wave is determined:

[00009] $f = E α_{1} \frac{1}{r} e^{- α_{2} r}$ [0044] where, α.sub.1 denotes an amplitude-energy ratio coefficient; α.sub.2 denotes an attenuation coefficient; and r denotes a distance from a micro-seismic source to the underground sensor. [0045] 402) A three-dimensional (3D) equidistant grid model including

[00010] $m_{1} \times n_{1} \times p_{1} = floor (\frac{X \max - X \min}{dx} + 1) \times floor (\frac{Y \max - Y \min}{dy} + 1) \times floor (\frac{Z \max - Z \min}{dz} + 1)$ [0046] grids with an X-direction spacing dx, a Y-direction spacing dy, and a Z-direction spacing dz is formed for a mining and production area defined by [Xmin,Xmax], [Ymin,Ymax], and [Zmin,Zmax], where m.sub.1, n.sub.1, and p.sub.1 denote a number of X-direction grids, a number of Y-direction grids, and a number of Z-direction grids, respectively. [0047] 403) The propagation relation equation determined in step 401) is rewritten to obtain

[00011] $E = \frac{fr}{α_{1} e^{- α_{2} r}};$ [0048] and minimum micro-seismic energy E.sub.i,j,k.sup.min to trigger the borehole-surface micro-seismic monitoring network deployment plan G.sub.v to record the micro-seismic signal, at point (X.sub.i, Y.sub.j, Z.sub.k) in the 3D equidistant grid model formed in step 402) is calculated.

[0049] i∈1,2, . . . , m1; j∈1,2, . . . , n1; k∈1,2, . . . , p1, v∈1,2, . . . , p. [0050] 404) According to step 403), an evaluation model for a monitoring capability Q.sub.v of the borehole-surface micro-seismic monitoring network deployment plan G.sub.v is established.

[00012] $Q_{v} = \frac{{.Math.}_{i = 1}^{m_{1}} {.Math.}_{j = 1}^{n_{1}} {.Math.}_{k = 1}^{p_{1}} E_{i, j, k}^{\min}}{m_{1} n_{1} p_{1}}$ [0051] 405) The initial population Gen is formed as follows:

[00013] $Gen = [\begin{matrix} G_{1} & Q_{1} \\ G_{2} & Q_{2} \\ .Math. & .Math. \\ G_{p} & Q_{p} \end{matrix}]$ [0052] (5) According to the initial population formed in step (4), an optimal borehole-surface micro-seismic monitoring network deployment plan is determined through a genetic algorithm and an optimal surface wireless sensor deployment plan is determined.

[0053] Further, in step (1), the multiple candidate points selected for installing the surface wireless sensors satisfy the following conditions: the candidate points cooperate with the underground sensors to surround the mining and production area; the candidate points have a distance of no more than 2,000 m from the mining and production area; the candidate points avoid a waterlogged area, a surface water system, a highway facility, and a noisy place; and the candidate points provide strong fourth-generation/fifth-generation (4G/5G) network signals.

[0054] Further, in step 401), α.sub.1 and α.sub.2 are determined as follows. The first-arrival peak amplitude f of the P-wave recorded by each underground sensor is manually marked and a source position and micro-seismic energy E of multiple micro-seismic signals with different energy levels are calculated. A distance r from the micro-seismic source to the underground sensor is calculated. α.sub.1 and α.sub.2 are determined by a nonlinear least squares (NLS) method.

[0055] Further, in step 403), the minimum micro-seismic energy E.sub.i,j,k.sup.min to trigger the borehole-surface micro-seismic monitoring network deployment plan G.sub.v to record the micro-seismic signal, at a point (X.sub.i, Y.sub.j, Z.sub.k) in the grid model is calculated as follows: [0056] 40301. According to a micro-seismic positioning principle based on a first arrival time of the P-wave, it is determined that the micro-seismic monitoring system is triggered to record the micro-seismic signal when the first-arrival peak amplitude f of the P-wave received by at least four sensors is greater than or equal to three times an ambient noise level NL. [0057] 40302. A distance r.sub.1 from the point (X.sub.i, Y.sub.j, Z.sub.k)) to each sensor in the borehole-surface micro-seismic monitoring network deployment plan G.sub.v is calculated. According to step 40301, a first-arrival peak amplitude f.sub.l of the P-wave required to trigger each sensor is determined. According to the propagation relation equation between the micro-seismic energy E and the first-arrival peak amplitude f of the P-wave determined in step 401, micro-seismic energy E.sub.i,j,k.sup.l required to trigger each sensor is back-calculated.

[00014] $E_{i, j, k}^{l} = \frac{f_{l} r_{l}}{α_{1} e^{- α_{2} r_{l}}},$ [0058] where, l=1,2, . . . , m+k [0059] 40303. According to step 40301, the micro-seismic energy E.sub.i,j,k.sup.l calculated in step 40302 is sorted in an ascending order. A fourth micro-seismic energy after the sorting is selected as the minimum micro-seismic energy E.sub.i,j,k.sup.min to trigger the borehole-surface micro-seismic monitoring network deployment plan to record the micro-seismic signal.

[0060] Further, in step 40301, the ambient noise level NL includes a surface ambient noise level NL.sub.s, monitored by the surface sensor installed on a surface and an underground ambient noise level NL.sub.u monitored by the underground sensor installed in an underground roadway.

[0061] Further, in step (5), the genetic algorithm sets a generation number of not less than 100; and the genetic algorithm carries out mutation operation through a mixture of adjacent gene mutation, gene insertion mutation, gene exchange mutation, three-point gene exchange mutation, and two-point inversion mutation, and carries out crossover operation through a mixture of partially mapped crossover, cycle crossover operator, edge recombination crossover, linear sequential crossover, ordered crossover operator, and uniform crossover.

Embodiment

[0062] (1) Multiple candidate points for installing surface wireless sensors are selected to form a natural-number-coded candidate point set S={1, 2, 3, 4, 5, . . . , 72}. The coordinates of each candidate point in the set are provided in the following table.

TABLE-US-00001 SN x/m y/m z/m 1 19381846 4321678 1400 2 19381846 4321878 1400 3 19381846 4322078 1400 4 19381846 4322278 1400 5 19381846 4322478 1400 6 19381846 4322678 1400 7 19381846 4322878 1400 8 19381846 4323078 1400 9 19382046 4321678 1400 10 19382046 4321878 1400 11 19382046 4322078 1400 12 19382046 4322278 1400 13 19382046 4322478 1400 14 19382046 4322678 1400 15 19382046 4322878 1400 16 19382046 4323078 1400 17 19382246 4321678 1400 18 19382246 4321878 1400 19 19382246 4322078 1400 20 19382246 4322278 1400 21 19382246 4322478 1400 22 19382246 4322678 1400 23 19382246 4322878 1400 24 19382246 4323078 1400 25 19382446 4321678 1400 26 19382446 4321878 1400 27 19382446 4322078 1400 28 19382446 4322278 1400 29 19382446 4322478 1400 30 19382446 4322678 1400 31 19382446 4322878 1400 32 19382446 4323078 1400 33 19382646 4321678 1400 34 19382646 4321878 1400 35 19382646 4322078 1400 36 19382646 4322278 1400 37 19382646 4322478 1400 38 19382646 4322678 1400 39 19382646 4322878 1400 40 19382646 4323078 1400 41 19382846 4321678 1400 42 19382846 4321878 1400 43 19382846 4322078 1400 44 19382846 4322278 1400 45 19382846 4322478 1400 46 19382846 4322678 1400 47 19382846 4322878 1400 48 19382846 4323078 1400 49 19383046 4321678 1400 50 19383046 4321878 1400 51 19383046 4322078 1400 52 19383046 4322278 1400 53 19383046 4322478 1400 54 19383046 4322678 1400 55 19383046 4322878 1400 56 19383046 4323078 1400 57 19383246 4321678 1400 58 19383246 4321878 1400 59 19383246 4322078 1400 60 19383246 4322278 1400 61 19383246 4322478 1400 62 19383246 4322678 1400 63 19383246 4322878 1400 64 19383246 4323078 1400 65 19383446 4321678 1400 66 19383446 4321878 1400 67 19383446 4322078 1400 68 19383446 4322278 1400 69 19383446 4322478 1400 70 19383446 4322678 1400 71 19383446 4322878 1400 72 19383446 4323078 1400 [0063] (2) A fixed number, namely 5, of candidate points randomly selected from the candidate point set S formed in step (1) are combined with an underground installed sensor set U to form a borehole-surface micro-seismic monitoring network deployment plan G.sub.1.

G.sub.1=[S.sub.2.sup.1S.sub.4.sup.2S.sub.20.sup.3S.sub.34.sup.4S.sub.69.sup.5U.sub.1U.sub.2U.sub.3U.sub.4]

[0064] In the plan, S.sub.2.sup.1 denotes a candidate point that is a 1.sup.st candidate point randomly selected from the candidate point set S and is a 2.sup.nd candidate point in the candidate point set S; similarly, S.sub.69.sup.5 denotes a candidate point that is a 5.sup.th candidate point randomly selected from the candidate point set S and is a 69.sup.th candidate point in the candidate point set S; and k=4 denotes a number of underground sensors. The coordinates of the sensors are provided in the following table.

TABLE-US-00002 U x/m y/m z/m U.sub.1 19382780.65 4322747.45 693.80 U.sub.2 19382493.17 4322201.48 696.80 U.sub.3 19382492.86 4322681.22 692.80 U.sub.4 19382780.95 4322183.56 697.00 [0065] (3) Step (2) is repeated until v=100 borehole-surface micro-seismic monitoring network deployment plans are generated to form a p=100-scale deployment plan set G:

[00015] $G = [\begin{matrix} G_{1} \\ G_{2} \\ .Math. \\ .Math. \\ G_{100} \end{matrix}] = [\begin{matrix} S_{2}^{1} & S_{4}^{2} & S_{2 0}^{3} & S_{3 4}^{4} & S_{6 9}^{5} & U_{1} & U_{2} & U_{3} & U_{4} \\ S_{1}^{1} & S_{8}^{2} & S_{3 7}^{3} & S_{6 5}^{4} & S_{7 2}^{5} & U_{1} & U_{2} & U_{3} & U_{4} \\ .Math. \\ .Math. \\ S_{5}^{1} & S_{3 3}^{2} & S_{3 7}^{3} & S_{4 0}^{4} & S_{6 8}^{5} & U_{1} & U_{2} & U_{3} & U_{4} \end{matrix}]$ [0066] (4) Each borehole-surface micro-seismic monitoring network deployment plan G.sub.v generated in step (3) is combined to form an initial population Gen. [0067] 401) According to a micro-seismic signal acquired by the underground sensor, a propagation relation equation between a micro-seismic energy E and a first-arrival peak amplitude f of a P-wave is determined as follows:

[00016] $f = E α_{1} \frac{1}{r} e^{- α_{2} r}$ [0068] α.sub.1 denotes an amplitude-energy ratio coefficient; α.sub.2 denotes an attenuation coefficient; and r denotes a distance from a micro-seismic source to the underground sensor. α.sub.1 and α.sub.2 are determined as follows. Three micro-seismic signals of different energy levels are selected and their micro-seismic energy and source position are calculated as follows:

TABLE-US-00003 Micro- Micro- seismic Source position seismic signal x/m y/m z/m energy E/J 1 19382734.46 4322782.88 727.02 2.60E+2 2 19382591.66 4322777.12 746.85 7.74E+3 3 19382648.81 4322816.84 752.25 4.10E+4

[0069] The first-arrival peak amplitude f of the P-wave recorded by the underground sensors is manually marked and a distance r from the micro-seismic source to the underground sensor is calculated as follows:

TABLE-US-00004 First-arrival amplitude Distance from source Micro-seismic SN of P-wave f/(m/s) to sensor r/m energy E/J 1 2.9281E−05 67.025 2.60E+2 2 2.9587E−06 264.341 2.60E+2 3 8.4109E−07 630.206 2.60E+2 4 1.70725E−07 753.718 2.60E+2 5 1.00089E−04 147.918 7.74E+3 6 4.61004E−05 198.524 7.74E+3 7 1.18401E−05 586.146 7.74E+3 8 3.15171E−06 625.003 7.74E+3 9 1.48796E−04 160.041 4.10E+4 10 7.20168E−05 215.052 4.10E+4 11 1.92005E−05 637.155 4.10E+4 12 6.50699E−06 649.274 4.10E+4

[0070] The data are brought into the propagation relation equation between the micro-seismic energy E and the first-arrival peak amplitude f of a P-wave.

[00017] $f = E α_{1} \frac{1}{r} e^{- α_{2} r}$

[0071] According to the NLS method, it is determined that α.sub.1 and α.sub.2 are 7.28352×10.sup.−7 and 0.0012836 respectively and a comparison between the original data and fitted data of the first-arrival peak amplitude f of the P-wave is plotted, as shown in FIG. 2. [0072] 402) As shown in FIG. 3, a 3D equidistant grid model including m.sub.1×n.sub.1×p.sub.1=14×10×5 grids with an X-direction spacing dx=50, a Y-direction spacing dy=50, and a Z-direction spacing dz=50 for a mining and production area defined by [Xmin=19382312,Xmax=19382962], [Ymin=4322587,Ymax=4323037], and [Zmin=600,Zmax=800]. [0073] 403) According to the propagation relation equation determined in step 401), the micro-seismic energy is:

[00018] $E = \frac{f r}{7.2 8 3 5 2 \times 1 0^{- 7} e^{- 0.0 0 1 2 8 3 6 r}},$

[0074] The minimum micro-seismic energy E.sub.i,j,k.sup.min to trigger the borehole-surface micro-seismic monitoring network deployment plan G.sub.v to record the micro-seismic signal, at the point (X.sub.i, Y.sub.j, Z.sub.k) in the 3D equidistant grid model formed in step 402) is calculated, i∈1,2, . . . , 14; j∈1,2, . . . , 10; k∈1,2, . . . , 5; v∈1,2, . . . 100. The minimum micro-seismic energy E.sub.1,1,1.sup.min to trigger the borehole-surface micro-seismic monitoring network deployment plan G.sub.2=G.sub.2=[S.sub.1.sup.1S.sub.8.sup.2S.sub.37.sup.3S.sub.65.sup.4S.sub.72.sup.5U.sub.1U.sub.2U.sub.3U.sub.4] to record the micro-seismic signal at the point [X.sub.1=19382312, Y.sub.1=4322587, Z.sub.1=600] is calculated as follows: [0075] 40301. According to a micro-seismic positioning principle based on a first arrival time of the P-wave, it is determined that the micro-seismic monitoring system is triggered to record the micro-seismic signal when the first-arrival peak amplitude f of the P-wave received by at least four sensors is greater than or equal to three times an ambient noise level NL. The ambient noise level NL includes a surface ambient noise level NL.sub.s=2.0×10.sup.−8 m/s monitored by the surface sensor installed on a surface and an underground ambient noise level NL.sub.u=5.0×10.sup.−7 m/s monitored by the underground sensor installed in an underground roadway. Specifically, the surface sensor satisfies f≥3×2.0×10.sup.−8 m/s=6×10.sup.−8 m/s, and the underground sensor satisfies f≥3×5.0×10.sup.−7 m/s=1.5×10.sup.−6 m/s. [0076] 40302. A distance r={1297, 1048, 874, 1659, 1472, 504, 437, 224, 626} from the point [X.sub.1=19382312, Y.sub.1=4322857, Z.sub.1=600] to each sensor in the borehole-surface micro-seismic monitoring network deployment plan G.sub.2 is calculated. According to step 40301, a first-arrival peak amplitude of the P-wave required to trigger each sensor is determined as follows:

f={6×10.sup.−8,6×10.sup.−8,6×10.sup.−8,6×10.sup.−8,1.5×10.sup.−6,1.5×10.sup.−6,1.5×10.sup.−6,1.5×10 .sup.−6} [0077] 40303. According to the propagation relation equation

[00019] $\frac{f_{l} \times r_{l}}{7.2 8 3 5 2 \times 1 0^{- 7} \times e^{- 0.0 0 1 2 836 \times r_{l}}}$

between the micro-seismic energy E and the first-arrival peak amplitude f of the P-wave determined in step 401, a micro-seismic energy required to trigger each sensor is back-calculated as follows:

[0078] E.sub.1,1,1.sup.l={565.2, 331.4, 220.9, 1149.4, 802.4, 1983.2, 1576.1, 615.2, 2880.8}.The micro-seismic energy calculated is sorted in an ascending order. A fourth micro-seismic energy in {220.9, 331.4, 565.2, 615.2, 802.4 , 1149.4, 1576.1, 1983.2, 2880.8} is selected as the minimum micro-seismic energy E.sub.1,1,1.sup.min=615.2 to trigger the borehole-surface micro-seismic monitoring network deployment plan to record the micro-seismic signal.

[0079] Steps 40302 and 40303 are repeated until the E.sub.i,j,k.sup.min at all points in the grid model is calculated. [0080] 404) According to step 403), an evaluation model for a monitoring capability Q.sub.v of the borehole-surface micro-seismic monitoring network deployment plan G.sub.v is established.

[00020] $Q_{v} = \frac{{.Math.}_{i = 1}^{m_{1}} {.Math.}_{j = 1}^{n_{1}} {.Math.}_{k}^{p_{1}} E_{i, j, k}^{\min}}{m_{1} n_{1} p_{1}}$

[0081] For example, the monitoring capability of the borehole-surface micro-seismic monitoring network deployment plan G.sub.2 is calculated as follows:

[00021] $Q_{2} = \frac{{.Math.}_{i = 1}^{1 4} {.Math.}_{j = 1}^{1 0} {.Math.}_{k = 1}^{5} E_{i, j, k}^{\min}}{1 4 \times 1 0 \times 5} = 6 3 8.6 7$

[0082] The monitoring capabilities of the borehole-surface micro-seismic monitoring network deployment plan set G are calculated as follows:

[00022] $Q = [\begin{matrix} 5 0 3.6 6 \\ 6 3 8.6 7 \\ .Math. \\ 486.99 \end{matrix}]$ [0083] 405) The initial population Gen is formed as follows:

[00023] $Gen = [\begin{matrix} S_{2}^{1} & S_{4}^{2} & S_{20}^{3} & S_{34}^{4} & S_{69}^{5} & U_{1} & U_{2} & U_{3} & U_{4} & 503.66 \\ S_{1}^{1} & S_{8}^{2} & S_{37}^{3} & S_{65}^{4} & S_{72}^{5} & U_{1} & U_{2} & U_{3} & U_{4} & 638.67 \\ .Math. & .Math. \\ S_{5}^{1} & S_{33}^{2} & S_{37}^{3} & S_{40}^{4} & S_{68}^{5} & U_{1} & U_{2} & U_{3} & U_{4} & 486.99 \end{matrix}]$ [0084] (5) As shown in FIG. 4, an optimal individual of the initial population Gen generated in step 405) is determined by the genetic algorithm. The minimum Q acquired by the evaluation model is 179.9. The genetic algorithm sets the generation number to be 200. In each generation, the genetic algorithm carries out mutation operation through a mixture of adjacent gene mutation, gene insertion mutation, gene exchange mutation, three-point gene exchange mutation, and two-point inversion mutation, and carries out crossover operation through a mixture of partially mapped crossover, cycle crossover operator, edge recombination crossover, linear sequential crossover, ordered crossover operator, and uniform crossover. Finally, the present disclosure determines an optimal surface wireless sensor deployment plan that significantly improves the monitoring capability of the borehole-surface micro-seismic monitoring network deployment plan, as shown in the following table.

TABLE-US-00005 S x/m y/m z/m S.sub.39 19382646 4322878 1400 S.sub.38 19382646 4322678 1400 S.sub.31 19382446 4322878 1400 S.sub.47 19382846 4322878 1400 S.sub.30 19382446 4322678 1400

Method for improving monitoring capability of borehole-surface micro-seismic monitoring system

Assignee

Inventors

Cpc classification

Classification Explorer

G06Q50/02

PHYSICS

Classification Explorer

G01V1/288

PHYSICS

Classification Explorer

G06F30/27

PHYSICS

Classification Explorer

G06F2111/06

PHYSICS

International classification

Classification Explorer

G06F30/27

PHYSICS

Classification Explorer

G01V1/28

PHYSICS

Classification Explorer

G06Q50/02

PHYSICS

Abstract

Claims

Description