CN109738940A - A kind of sound emission there are under the conditions of dead zone/microseismic event localization method - Google Patents
A kind of sound emission there are under the conditions of dead zone/microseismic event localization method Download PDFInfo
- Publication number
- CN109738940A CN109738940A CN201910047005.2A CN201910047005A CN109738940A CN 109738940 A CN109738940 A CN 109738940A CN 201910047005 A CN201910047005 A CN 201910047005A CN 109738940 A CN109738940 A CN 109738940A
- Authority
- CN
- China
- Prior art keywords
- sound emission
- dead zone
- unit
- microseismic event
- wave
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Abstract
A kind of sound emission there are under the conditions of dead zone/microseismic event localization method provided by the invention, it overcomes since the difficulties of positioning accuracy are not allowed and are influenced in velocity of wave input under the influence of rock excavation, it can be on the basis of considering that rock engineering digging process influences, in conjunction with the numerical simulation of stress wave propagation, it is more accurately determined sound emission/microseismic event spatial position in lithosome.Using established space status as known conditions, the lag time that dead zone reaches sensor to the variation of the change of propagation path and wave resulting from propagation distance, wave is quantified in conjunction with numerical simulation, is significantly reduced because dead zone has the sound emission/microseismic event position error induced.In addition, grid search and Error weight Y-factor method Y have been used in combination in location algorithm, the shortcomings that excessively relying on primary iteration value in iterative algorithm and not restraining not only had been avoided, but also has breached limitation of the mesh scale to positioning result.
Description
Technical field
The invention belongs to sound emission in mining/microseismic event localization method technical field, it is related to that a kind of there are dead zones
Under the conditions of sound emission/microseismic event localization method.
Background technique
The rapid growth of China's economy greatly have stimulated the exploitation of resource, the energy, and wherein most is related to rock engineering.
The work progress of rock engineering can cause the formation in lithosome internal rupture face and simultaneously with the propagation of stress wave, this phenomenon
Referred to as sound emission (frequency is high, amplitude is small) or microseism (frequency is low, amplitude is big).Sound emission/microseism is lithosome rupture process
Attendant phenomenon, physical mechanics behavior with lithosome has close relationship therefore can be by sound emission/micro seismic monitoring
The stress and level of breakage of information analyzed to infer lithosome, and then the de-stabilise of early warning and control lithosome destroys.
Sound emission/On Microseismic Monitoring Technique is gradually paid attention to by countries in the world in recent years, and a large amount of indoor rocks at home and abroad
It is applied in stone experiment and rock engineering.National Security of China supervision general bureau has put into effect in metal and nonmetal sub-terrain mines for 2010
Mountain must install the regulation of safe " six big systems ", and sound emission/On Microseismic Monitoring Technique conduct can meet ground in " six big systems "
Press monitoring and controlling system related request powerful oneself gradually started large-scale application in mine at home.
Sound emission/microseismic event positioning is one of core function of the technology, can carry out accurately positioning be sound emission/
The critical evaluation index whether Microseismic monitoring system effectively plays a role.In mining engineering, for economic benefit often more water
Flat, hot spot is exploited simultaneously, and since (high quality in a set of 8 channel monitors system to sound emission/Microseismic monitoring system fancy price
System is at 200,000 yuan or more) and target monitoring range (frequently including multiple dead zones or multiple levels) and monitoring effect (positioning accuracy
It is general to require in 10m or so) requirement, the laying interval of sensor is generally in 50-120m or so.Therefore, in a large amount of lithosomes
There are dead zone barrier between the rupture location and sensor in portion, the stress wave that lithosome rupture induces must be around dead zone
Sensor position can be propagated to be collected, propagation path and the change for propagating distance will lead to the change that wave propagates duration,
Ignoring this influence factor will lead to the increase of positioning result error, but the location algorithm being widely used at present does not account for dead zone
Influence to propagation path.
Summary of the invention
The object of the present invention is to provide a kind of sound emission there are under the conditions of dead zone/microseismic event localization methods, overcome
Since the difficulties of positioning accuracy are not allowed and are influenced in velocity of wave input under the influence of rock excavation, can be excavated considering rock engineering
On the basis of process influences, in conjunction with the numerical simulation of stress wave propagation, it is more accurately determined sound emission in lithosome/microseism thing
The spatial position of part.
The present invention provides a kind of sound emission there are under the conditions of dead zone/microseismic event localization method, includes the following steps:
The foundation and assignment of step 1, the numerical model containing dead zone: establishing 3D solid geological model, grid division and to from
The grid cell for belonging to different lithology carries out the assignment of physical and mechanical parameter;
Step 2, the numerical simulation of stress wave propagation: side is applied to the 3D solid geological model established in step 1
Boundary's condition and instantaneous step force, and simulate and calculate matrix when stress wave is walked;
Step 3, sound emission/microseismic event positioning: establishing arrival time difference matrix/vector according to matrix when walking, matching search and
Sound emission/microseismic event positioning.
Of the invention there are in sound emission under the conditions of dead zone/microseismic event localization method, three are established in the step 1
Dimension entity geological model specifically includes:
Step 1.1 is generated log sheet according to bore database and is adjusted by its true coordinate to corresponding space
In position, the top of same lithology, bottom plate are respectively connected with, form the log sheet of each lithology, passes through finger using setting-out order
The section of series of identical lithology is determined to create the 3D solid geological model of different lithology;
Step 1.2, the design size according to goaf, length establish simple cuboid dead zone model, if empty
The actual form in area is more complicated and greatly differs from each other with cuboid, then can use three-dimensional laser scanner cavity monitoring system pair
The form of dead zone is reconstructed, to establish dead zone model;
3D solid geological model and dead zone model are imported in numerical simulation software and do difference set boolean fortune by step 1.3
It calculates, i.e. 3D solid geological model-dead zone model, completes the foundation of the numerical model containing dead zone.
Of the invention there are in sound emission under the conditions of dead zone/microseismic event localization method, net is divided in the step 1
Lattice specifically:
Step 1.4, to containing dead zone numerical model carry out grid dividing, will be close to dead zone mesh scale divide compared with
It is small, and larger, the mesh scale around dead zone is divided further away from the mesh scale of dead zone are as follows:
In formula, LMFor mesh scale, unit m;VPFor spread speed of the P wave in rock mass, unit m/s;SfFor sound emission/
The sample frequency of Microseismic monitoring system, unit Hz;emaxFor the worst error for meeting positioning requirements, unit m.
Of the invention there are in sound emission under the conditions of dead zone/microseismic event localization method, object is carried out in the step 1
Manage the assignment of mechanics parameter specifically:
Step 1.5 carries out assignment to the grid cell in different lithology by numerical simulation software, due to most Numerical-Modes
Quasi- software can not directly assign grid cell value of wave speed, can be by assigning grid cell elasticity modulus and density come indirectly to unit
Value of wave speed is assigned, such as following formula of the relationship between velocity of wave and elasticity modulus, density indicates:
In formula, E is the elasticity modulus of lithosome, unit Pa;ρ is the density of lithosome, units/kg/m3。
Of the invention there are in sound emission under the conditions of dead zone/microseismic event localization method, the step 2 specifically:
Step 2.1 applies boundary condition: reaching sensor since sound emission/microseismic event positioning need to only obtain stress wave
At the time of position, so improving settlement efficiency to reduce calculation amount, the subsequent stress wave of excessive interference of reflected wave being avoided to reach
The sensor position moment accurately identifies, and the boundary of numerical model should be set as wholly transmissive boundary;
Step 2.2 applies instantaneous step force: using the node of each grid cell as representing the sound emission of the position/micro-
Focus Pi(i=l, 2 ..., N), position coordinates Xi, Yi, Zi, successively in each sound emission/microquake sources PiPlace applies a wink
When step force Fi(i=l, 2 ..., N), the time interval that instantaneous step force applies should ensure that previous instantaneous step force induced
Stress wave has propagated to all sensors, gives time interval according to the following formula:
In formula, Δ t is the time interval that instantaneous step force applies, unit s;L, W, H be respectively numerical model length and width,
Height, unit m;
Step 2.3, simulation calculate matrix when stress wave is walked: simulation calculates the stress wave induced by instantaneous step force in numerical value
Communication process in model, it is first determined sound emission/microseismic sensors corresponding coordinate s in numerical modeli(xi, yi, zi),
Wherein i=l, 2 ..., n, n are number of sensors, then determine the application moment of instantaneous step force and are lured by the instantaneous step force
At the time of the stress wave propagation of hair to each sensor position, then matrix when walking of stress wave to each sensor are as follows:
Tt=Tij=Mij-MPi
In formula, TtFor matrix when walking of stress wave to each sensor, unit s;TijFor sound emission/microquake sources PiIt induces
Stress wave propagation is to j sensor when walking, unit s;MijSound emission/microquake sources P is received for j sensoriInduce stress
At the time of wave, unit s;MPiFor sound emission/microseism source point PiAt the time of applying instantaneous step force, unit s.
Of the invention there are in sound emission under the conditions of dead zone/microseismic event localization method, sound is sent out in the step 3
/ microseismic event positioning is penetrated, specifically includes the following steps:
Step 3.1, the foundation for simulating arrival time difference matrix: matrix when walking obtained according to numerical simulation calculates successively from sound
Transmitting/microquake sources PiPropagate to the arrival time difference matrix of each sensor position:
Δ T=Tij-Timin
In formula, Δ T is arrival time difference matrix, unit s;TijFor sound emission/microquake sources PiThe stress wave propagation of induction is to No. j
When walking of sensor, TiminMatrix T when to walkijIn the smallest element in the i-th row, unit s.
The foundation of step 3.2, true then difference vector: it is acquired using the sound emission installed at the scene/Microseismic monitoring system
Wave data carries out manual or automatic then pickup work to the waveform that sensor receives, and calculates sensor according to the following formula
True then difference vector:
ΔTireal=Tireal-min(Tireal)
In formula, Δ TirealFor really then difference vector, unit s;TirealFor sensor i it is true then, unit s;
Step 3.3, matching search and sound emission/microseismic event positioning: will really then difference vector and simulation arrival time difference square
Battle array in row vector carry out matching search, if there is the sum of departure absolute value be 0 row vector, then directly determine sound emission/
The coordinate of microseismic event cell node coordinate corresponding with the row vector is identical, completes positioning;If it is not, determining and true
Then difference vector is closest, i.e. the smallest 4 row vectors of total deviation amount, i.e. sound emission/microquake sources Pi, the size of departure is made
Final anchor point is determined for weight coefficient:
In formula, c is that sound emission/microseismic event positions coordinate, unit m;ekIt is total between true arrival time difference and quasi- arrival time difference
Departure, unit s;P1, P2, P3, P4 are respectively and really then immediate 4 sound emission/microquake sources of difference vector.
A kind of sound emission there are under the conditions of dead zone/microseismic event localization method of the invention at least has beneficial below
Effect:
The method overcome due under the influence of rock excavation velocity of wave input it is inaccurate and influence the difficulties of positioning accuracy, energy
Enough on the basis of considering that rock engineering digging process influences, in conjunction with the numerical simulation of stress wave propagation, more it is accurately determined
Sound emission/microseismic event spatial position in lithosome.Using established space status as known conditions, in conjunction with numerical simulation
The variation of distance is propagated into the change of propagation path and wave resulting from dead zone, the lag time of wave arrival sensor adds
With quantization, significantly reduce because dead zone has the sound emission/microseismic event position error induced.In addition, joining in location algorithm
Conjunction has used grid search and Error weight Y-factor method Y, had both avoided and has excessively relied on primary iteration value in iterative algorithm and do not restrain
The shortcomings that, and breach limitation of the mesh scale to positioning result.
Detailed description of the invention
Fig. 1 is the numerical model containing dead zone established;
Fig. 2 is the numerical model after grid division;
Fig. 3 a is the numerical simulation result that stress wave does not reach dead zone;
Fig. 3 b is that stress wave reaches dead zone and generates the numerical simulation result of reflection;
Fig. 3 c is that stress wave will bypass dead zone, because variation diagram has occurred in the shape that dead zone influences wavefront;
Fig. 4 is certain actual measurement acoustic emission waveform schematic diagram.
Specific embodiment
The present invention is further introduced below in conjunction with the drawings and specific embodiments.
A kind of sound emission there are under the conditions of dead zone/microseismic event localization method of the invention, includes the following steps:
The foundation and assignment of step 1, the numerical model containing dead zone: establishing 3D solid geological model, grid division and to from
The grid cell for belonging to different lithology carries out the assignment of physical and mechanical parameter;
3D solid geological model is established in the step 1 to specifically include:
Step 1.1 is generated log sheet according to bore database and is adjusted by its true coordinate to corresponding space
In position, the top of same lithology, bottom plate are respectively connected with, form the log sheet of each lithology, passes through finger using setting-out order
The section of series of identical lithology is determined to create the 3D solid geological model of different lithology;
Step 1.2 establishes simple cuboid dead zone model according to the design size (length) in goaf, if empty
The actual form in area is more complicated and greatly differs from each other with cuboid, then can use three-dimensional laser scanner cavity monitoring system pair
The form of dead zone is reconstructed, to establish dead zone model;
3D solid geological model and dead zone model are imported in numerical simulation software and do difference set boolean fortune by step 1.3
It calculates, i.e. 3D solid geological model-dead zone model, completes the foundation of the numerical model containing dead zone.
Grid division in the step 1, specifically:
Step 1.4, to containing dead zone numerical model carry out grid dividing, will be close to dead zone mesh scale divide compared with
It is small, and larger, the mesh scale around dead zone is divided further away from the mesh scale of dead zone are as follows:
In formula, LMFor mesh scale, unit m;VPFor spread speed of the P wave in rock mass, unit m/s;SfFor sound emission/
The sample frequency of Microseismic monitoring system, unit Hz;emaxFor the worst error for meeting positioning requirements, unit m.
The assignment of physical and mechanical parameter is carried out in the step 1 specifically:
Step 1.5 carries out assignment to the grid cell in different lithology by numerical simulation software, due to most Numerical-Modes
Quasi- software can not directly assign grid cell value of wave speed, can be by assigning grid cell elasticity modulus and density come indirectly to unit
Value of wave speed is assigned, such as following formula of the relationship between velocity of wave and elasticity modulus, density indicates:
In formula, E is the elasticity modulus of lithosome, unit Pa;ρ is the density of lithosome, units/kg/m3。
Step 2, the numerical simulation of stress wave propagation: side is applied to the 3D solid geological model established in step 1
Boundary's condition and instantaneous step force, and simulate and calculate matrix when stress wave is walked, specifically:
Step 2.1 applies boundary condition: reaching sensor since sound emission/microseismic event positioning need to only obtain stress wave
At the time of position, so improving settlement efficiency to reduce calculation amount, the subsequent stress wave of excessive interference of reflected wave being avoided to reach
The sensor position moment accurately identifies, and the boundary of numerical model should be set as wholly transmissive boundary;
Step 2.2 applies instantaneous step force: using the node of each grid cell as representing the sound emission of the position/micro-
Focus Pi(i=l, 2 ..., N), position coordinates are (Xi, Yi, Zi), successively in each sound emission/microquake sources PiPlace applies primary
Instantaneous step force Fi(i=l, 2 ..., N), the time interval that instantaneous step force applies should ensure that previous instantaneous step force induces
Stress wave propagated to all sensors, give time interval according to the following formula:
In formula, Δ t is the time interval that instantaneous step force applies, unit s;L, W, H be respectively numerical model length and width,
Height, unit m;
Step 2.3, simulation calculate matrix when stress wave is walked: simulation calculates the stress wave induced by instantaneous step force in numerical value
Communication process in model, it is first determined sound emission/microseismic sensors corresponding coordinate s in numerical modeli(xi, yi, zi),
Wherein i=l, 2 ..., n, n are number of sensors, then determine the application moment of instantaneous step force and are lured by the instantaneous step force
At the time of the stress wave propagation of hair to each sensor position, then matrix when walking of stress wave to each sensor are as follows:
Tt=Tij=Mij-MPi
In formula, TtFor matrix when walking of stress wave to each sensor, unit s;TijFor sound emission/microquake sources PiIt induces
Stress wave propagation is to j sensor when walking, unit s;MijSound emission/microquake sources P is received for j sensoriInduce stress
At the time of wave, unit s;MPiFor sound emission/microseism source point PiAt the time of applying instantaneous step force, unit s.
Step 3, sound emission/microseismic event positioning: establishing arrival time difference matrix/vector according to matrix when walking, matching search and
Sound emission/microseismic event positioning, specifically includes the following steps:
Step 3.1, the foundation for simulating arrival time difference matrix: matrix when walking obtained according to numerical simulation calculates successively from sound
Transmitting/microquake sources PiPropagate to the arrival time difference matrix of each sensor position:
Δ T=Tij-Timin
In formula, Δ T is arrival time difference matrix, unit s;TijFor sound emission/microquake sources PiThe stress wave propagation of induction is to No. j
When walking of sensor, TiminMatrix T when to walkijIn the smallest element in the i-th row, unit s.
The foundation of step 3.2, true then difference vector: it is acquired using the sound emission installed at the scene/Microseismic monitoring system
Wave data carries out manual or automatic then pickup work to the waveform that sensor receives, and calculates sensor according to the following formula
True then difference vector:
ΔTireal=Tireal-min(Tireal)
In formula, Δ TirealFor really then difference vector, unit s;TirealFor sensor i it is true then, unit s;
Step 3.3, matching search and sound emission/microseismic event positioning: will really then difference vector and simulation arrival time difference square
Battle array in row vector carry out matching search, if there is the sum of departure absolute value be 0 row vector, then directly determine sound emission/
The coordinate of microseismic event cell node coordinate corresponding with the row vector is identical, completes positioning;If it is not, determining and true
Then 4 row vectors (sound emission/microquake sources P of the difference vector closest to (total deviation amount is minimum)i), using the size of departure as
Weight coefficient determines final anchor point:
In formula, c is that sound emission/microseismic event positions coordinate, unit m;ekIt is total between true arrival time difference and quasi- arrival time difference
Departure, unit s;P1, P2, P3, P4 are respectively and really then immediate 4 sound emission/microquake sources of difference vector.
Embodiment:
1, establish 3D solid geological model according to the actual situation, as shown in Figure 1, suppose there is 4 sensor S1, S2, S3,
S4, practical sound emission/microseismic event c, sensor location coordinates are as shown in table 1.
Table 1
2, grid dividing is carried out to the 3D solid geological model of Fig. 1, as shown in Figure 2.
3, physical and mechanical parameter value is assigned to grid cell, for the sake of letter, using the uniform dielectric of unified lithology in the present embodiment.
4, boundary condition is set and successively applies instantaneous step force, instantaneous step force application time interval on cell node
For 0.1s.
5, stress wave propagation is simulated, according to analog result calculating stress wave (as shown in Figure 3) to each sensor
When walking, in the case where there are 4 sensors, each it is applied to the instantaneous step force of cell node just and can be obtained 4 when walking, if there is N
A cell node just will form matrix when walking of the column of N × 4.With coordinate be (200,200,0) cell node for, knot when walking
Fruit is as shown in table 2.
Table 2
Sensor | When walking |
S1 | 0.0169706 |
S2 | 0.0346987 |
S3 | 0.0346987 |
S4 | 0.0463500 |
6, table 2 is by the available simulation arrival time difference matrix of value least when subtracting away, as shown in table 3.
Table 3
Sensor | Arrival time difference matrix |
S1 | 0 |
S2 | 0.0177281 |
S3 | 0.0177281 |
S4 | 0.0293794 |
7, Wave data is acquired using the sound emission installed at the scene/Microseismic monitoring system, as shown in figure 4, to sensor
The waveform received carries out artificial then pickup work, and calculates true then difference vector.
8, the row vector really then in difference vector and simulation arrival time difference matrix is subjected to matching search, discovery is really arrived
When difference vector and simulation arrival time difference matrix in coordinate be (200,200,0) cell node induce stress wave when difference vector it is consistent
(total deviation amount is 0), it is thus determined that the coordinate of practical sound emission/microseismic event is (200,200,0).
The foregoing is merely presently preferred embodiments of the present invention, the thought being not intended to limit the invention, all of the invention
Within spirit and principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (6)
1. a kind of sound emission there are under the conditions of dead zone/microseismic event localization method, which comprises the steps of:
The foundation and assignment of step 1, the numerical model containing dead zone: establishing 3D solid geological model, grid division and to being subordinated to
The grid cell of different lithology carries out the assignment of physical and mechanical parameter;
Step 2, the numerical simulation of stress wave propagation: perimeter strip is applied to the 3D solid geological model established in step 1
Part and instantaneous step force, and simulate and calculate matrix when stress wave is walked;
Step 3, sound emission/microseismic event positioning: arrival time difference matrix/vector, matching search and sound hair are established according to matrix when walking
Penetrate/microseismic event positioning.
2. as described in claim 1, there are the sound emission under the conditions of dead zone/microseismic event localization methods, which is characterized in that institute
It states and establishes 3D solid geological model in step 1 and specifically include:
Step 1.1 is generated log sheet according to bore database and is adjusted by its true coordinate to corresponding spatial position
In, the top of same lithology, bottom plate are respectively connected with, the log sheet of each lithology is formed, using setting-out order by specifying one
The section of the identical lithology of series creates the 3D solid geological model of different lithology;
Step 1.2, the design size according to goaf, length establish simple cuboid dead zone model, if dead zone
Actual form is more complicated and greatly differs from each other with cuboid, then can use three-dimensional laser scanner cavity monitoring system to dead zone
Form be reconstructed, to establish dead zone model;
3D solid geological model and dead zone model are imported in numerical simulation software and do difference set Boolean calculation by step 1.3, i.e.,
3D solid geological model-dead zone model completes the foundation of the numerical model containing dead zone.
3. as described in claim 1, there are the sound emission under the conditions of dead zone/microseismic event localization methods, which is characterized in that institute
State grid division in step 1 specifically:
Step 1.4 carries out grid dividing to the numerical model containing dead zone, will be close to the smaller of the mesh scale division of dead zone, and
Mesh scale further away from dead zone divides larger mesh scale around dead zone are as follows:
In formula, LMFor mesh scale, unit m;VPFor spread speed of the P wave in rock mass, unit m/s;SfFor sound emission/microseism
The sample frequency of monitoring system, unit Hz;emaxFor the worst error for meeting positioning requirements, unit m.
4. as described in claim 1, there are the sound emission under the conditions of dead zone/microseismic event localization methods, which is characterized in that institute
State the assignment that physical and mechanical parameter is carried out in step 1 specifically:
Step 1.5 carries out assignment to the grid cell in different lithology by numerical simulation software, since most numerical simulations are soft
Part can not directly assign grid cell value of wave speed, can assign wave to unit indirectly by assigning grid cell elasticity modulus and density
Speed value, the relationship such as following formula between velocity of wave and elasticity modulus, density indicate:
In formula, E is the elasticity modulus of lithosome, unit Pa;ρ is the density of lithosome, units/kg/m3。
5. as described in claim 1, there are the sound emission under the conditions of dead zone/microseismic event localization methods, which is characterized in that institute
State step 2 specifically:
Step 2.1 applies boundary condition: reaching sensor position since sound emission/microseismic event positioning need to only obtain stress wave
At the time of, so improving settlement efficiency to reduce calculation amount, the subsequent stress wave of excessive interference of reflected wave being avoided to reach sensing
The device position moment accurately identifies, and the boundary of numerical model should be set as wholly transmissive boundary;
Step 2.2 applies instantaneous step force: using the node of each grid cell as the sound emission/microquake sources P for representing the positioni
(i=l, 2 ..., N), position coordinates Xi, Yi, Zi, successively in each sound emission/microquake sources PiPlace applies primary instantaneous step
Power Fi(i=l, 2 ..., N), the time interval that instantaneous step force applies should ensure that the stress wave that previous instantaneous step force induces
All sensors have been propagated to, have given time interval according to the following formula:
In formula, Δ t is the time interval that instantaneous step force applies, unit s;L, W, H are respectively the length of numerical model, single
Position m;
Step 2.3, simulation calculate matrix when stress wave is walked: simulation calculates the stress wave induced by instantaneous step force in numerical model
In communication process, it is first determined sound emission/microseismic sensors corresponding coordinate s in numerical modeli(xi, yi, zi), wherein i
=l, 2 ..., n, n are number of sensors, then determine the application moment of instantaneous step force and are answered by what the instantaneous step force induced
At the time of Reeb propagates to each sensor position, then matrix when walking of stress wave to each sensor are as follows:
Tt=Tij=Mij-MPi
In formula, TtFor matrix when walking of stress wave to each sensor, unit s;TijFor sound emission/microquake sources PiThe stress of induction
Wave propagates to when walking of j sensor, unit s;MijSound emission/microquake sources P is received for j sensoriInduce stress wave
Moment, unit s;MPiFor sound emission/microseism source point PiAt the time of applying instantaneous step force, unit s.
6. as described in claim 1, there are the sound emission under the conditions of dead zone/microseismic event localization methods, which is characterized in that institute
Sound emission in step 3/microseismic event positioning is stated, specifically includes the following steps:
Step 3.1, simulate arrival time difference matrix foundation: matrix when walking obtained according to numerical simulation, calculate successively from sound emission/
Microquake sources PiPropagate to the arrival time difference matrix of each sensor position:
Δ T=Tij-Timin
In formula, Δ T is arrival time difference matrix, unit s;TijFor sound emission/microquake sources PiThe stress wave propagation of induction to No. j sense
When walking of device, TiminMatrix T when to walkijIn the smallest element in the i-th row, unit s.
The foundation of step 3.2, true then difference vector: waveform is acquired using the sound emission installed at the scene/Microseismic monitoring system
Data carry out manual or automatic then pickup work to the waveform that sensor receives, and calculate the true of sensor according to the following formula
Difference vector when actual arrival:
ΔTireal=Tireal-min(Tireal)
In formula, Δ TirealFor really then difference vector, unit s;TirealFor sensor i it is true then, unit s;
Step 3.3, matching search and sound emission/microseismic event positioning: will be really then in difference vector and simulation arrival time difference matrix
Row vector carry out matching search, if there is the sum of departure absolute value be 0 row vector, then directly determine sound emission/microseism
The coordinate of event cell node coordinate corresponding with the row vector is identical, completes positioning;If it is not, it is determining with it is true then
Difference vector is closest, i.e. the smallest 4 row vectors of total deviation amount, i.e. sound emission/microquake sources Pi, using the size of departure as power
Coefficient is weighed to determine final anchor point:
In formula, c is that sound emission/microseismic event positions coordinate, unit m;ekFor the total deviation between true arrival time difference and quasi- arrival time difference
Amount, unit s;P1, P2, P3, P4 are respectively and really then immediate 4 sound emission/microquake sources of difference vector.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910047005.2A CN109738940B (en) | 2019-01-18 | 2019-01-18 | Acoustic emission/microseismic event positioning method under condition of existing empty zone |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910047005.2A CN109738940B (en) | 2019-01-18 | 2019-01-18 | Acoustic emission/microseismic event positioning method under condition of existing empty zone |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109738940A true CN109738940A (en) | 2019-05-10 |
CN109738940B CN109738940B (en) | 2021-01-29 |
Family
ID=66365192
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910047005.2A Active CN109738940B (en) | 2019-01-18 | 2019-01-18 | Acoustic emission/microseismic event positioning method under condition of existing empty zone |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109738940B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110907897A (en) * | 2019-12-23 | 2020-03-24 | 鞍钢集团矿业有限公司 | Method for positioning acoustic emission source suitable for rock with hole |
CN111880220A (en) * | 2020-09-07 | 2020-11-03 | 中国科学院武汉岩土力学研究所 | Seismic source positioning method, device, equipment and storage medium |
CN113153430A (en) * | 2021-03-23 | 2021-07-23 | 中国矿业大学 | Roadway surrounding rock damage acoustic emission positioning and wave velocity imaging monitoring and catastrophe early warning method |
CN113219068A (en) * | 2021-05-14 | 2021-08-06 | 重庆大学 | Cylinder acoustic emission positioning method, system, terminal and readable storage medium based on analysis solution of group sensor |
CN113552536A (en) * | 2021-07-30 | 2021-10-26 | 重庆大学 | Acoustic emission/microseismic event positioning method, system, terminal and readable storage medium containing round hole structure |
CN116593295A (en) * | 2023-07-19 | 2023-08-15 | 北京科技大学 | Method and device for improving acoustic emission positioning precision by utilizing rock anisotropic wave velocity |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2550770C1 (en) * | 2014-08-27 | 2015-05-10 | Открытое акционерное общество "Татнефть" им. В.Д. Шашина | Method to determine geometric characteristics of hydraulic fracturing crack |
CN106094021A (en) * | 2016-06-01 | 2016-11-09 | 北京科技大学 | A kind of microseism focus method for rapidly positioning based on arrival time difference data base |
CN107884822A (en) * | 2017-11-13 | 2018-04-06 | 北京矿冶研究总院 | Method for improving positioning precision of mining micro-seismic source |
-
2019
- 2019-01-18 CN CN201910047005.2A patent/CN109738940B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2550770C1 (en) * | 2014-08-27 | 2015-05-10 | Открытое акционерное общество "Татнефть" им. В.Д. Шашина | Method to determine geometric characteristics of hydraulic fracturing crack |
CN106094021A (en) * | 2016-06-01 | 2016-11-09 | 北京科技大学 | A kind of microseism focus method for rapidly positioning based on arrival time difference data base |
CN107884822A (en) * | 2017-11-13 | 2018-04-06 | 北京矿冶研究总院 | Method for improving positioning precision of mining micro-seismic source |
Non-Patent Citations (2)
Title |
---|
LI QI-YUE ET AL.: "Effects of sonic speed on location accuracy of acoustic emission source in rocks", 《TRANSACTIONS OF NONFERROUS METALS SOCIETY OF CHINA》 * |
朱双江: "浅部采空区微震定位方法研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110907897A (en) * | 2019-12-23 | 2020-03-24 | 鞍钢集团矿业有限公司 | Method for positioning acoustic emission source suitable for rock with hole |
CN110907897B (en) * | 2019-12-23 | 2023-09-15 | 鞍钢集团矿业有限公司 | Acoustic emission source positioning method suitable for rock containing holes |
CN111880220A (en) * | 2020-09-07 | 2020-11-03 | 中国科学院武汉岩土力学研究所 | Seismic source positioning method, device, equipment and storage medium |
CN113153430A (en) * | 2021-03-23 | 2021-07-23 | 中国矿业大学 | Roadway surrounding rock damage acoustic emission positioning and wave velocity imaging monitoring and catastrophe early warning method |
CN113153430B (en) * | 2021-03-23 | 2022-02-08 | 中国矿业大学 | Roadway surrounding rock damage acoustic emission positioning and wave velocity imaging monitoring and catastrophe early warning method |
CN113219068A (en) * | 2021-05-14 | 2021-08-06 | 重庆大学 | Cylinder acoustic emission positioning method, system, terminal and readable storage medium based on analysis solution of group sensor |
CN113219068B (en) * | 2021-05-14 | 2022-03-22 | 重庆大学 | Cylinder acoustic emission positioning method, system, terminal and readable storage medium based on analysis solution of group sensor |
CN113552536A (en) * | 2021-07-30 | 2021-10-26 | 重庆大学 | Acoustic emission/microseismic event positioning method, system, terminal and readable storage medium containing round hole structure |
CN113552536B (en) * | 2021-07-30 | 2022-08-09 | 重庆大学 | Acoustic emission/microseismic event positioning method, system, terminal and readable storage medium containing round hole structure |
CN116593295A (en) * | 2023-07-19 | 2023-08-15 | 北京科技大学 | Method and device for improving acoustic emission positioning precision by utilizing rock anisotropic wave velocity |
CN116593295B (en) * | 2023-07-19 | 2023-10-03 | 北京科技大学 | Method and device for improving acoustic emission positioning precision by utilizing rock anisotropic wave velocity |
Also Published As
Publication number | Publication date |
---|---|
CN109738940B (en) | 2021-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109738940A (en) | A kind of sound emission there are under the conditions of dead zone/microseismic event localization method | |
Feng et al. | Sectional velocity model for microseismic source location in tunnels | |
CA2886778C (en) | Propagating fracture plane updates | |
CA2969101C (en) | Seismic elastic wave simulation for tilted transversely isotropic media using adaptive lebedev staggered grid | |
CN106094021B (en) | A kind of microseism focus method for rapidly positioning based on arrival time difference database | |
Dong et al. | Some developments and new insights for microseismic/acoustic emission source localization | |
CN105301639B (en) | The method and its device of double weighting chromatography inversion speeds when being travelled based on VSP | |
CN109492262A (en) | A method of utilizing numerical simulation analysis non-uniform Distribution crack Drift stability | |
CN106814391A (en) | Ground micro-seismic state event location method based on Fresnel zone tomographic inversion | |
CN106772577A (en) | Source inversion method based on microseism data and SPSA optimized algorithms | |
CN109375253A (en) | Ground motion parameter evaluation method based on whole seismic structure maximum credible earthquakes | |
CN106980716B (en) | Underground cavern stability analysis method based on random block | |
CN109632016A (en) | Rock And Soil adit digging and surrouding rock stress, strain monitoring experimental rig and its method | |
CN105487117A (en) | Three-dimensional earthquake observation system optimization method and apparatus | |
Kan et al. | Study on influencing factors and prediction of peak particle velocity induced by roof pre-split blasting in underground | |
Wang et al. | Mechanical parameter inversion in sandstone diversion tunnel and stability analysis during operation period | |
CN113189644B (en) | Microseismic source positioning method and system | |
Li et al. | Anchoring parameters optimization of tunnel surrounding rock based on particle swarm optimization | |
CN111830557A (en) | Artificial fracture complexity index obtaining method and system based on fracturing microseism | |
Bondur et al. | Connection between variations of the stress-strain state of the Earth’s crust and seismic activity: The example of Southern California | |
Wang | Structure-soil-structure interaction between underground structure and surface structure | |
Huang et al. | Relocation method of microseismic source in deep mines | |
Fälth et al. | Simulation of co-seismic secondary fracture displacements for different earthquake rupture scenarios at the proposed nuclear waste repository site in Forsmark | |
CN114047546A (en) | Crowd-sourcing spiral mine earthquake positioning method based on three-dimensional spatial joint arrangement of sensors | |
Bi et al. | Research on seismic input method of layered ground foundation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |