CN107515419A - The method of estimation and device of rock stability - Google Patents
The method of estimation and device of rock stability Download PDFInfo
- Publication number
- CN107515419A CN107515419A CN201710740277.1A CN201710740277A CN107515419A CN 107515419 A CN107515419 A CN 107515419A CN 201710740277 A CN201710740277 A CN 201710740277A CN 107515419 A CN107515419 A CN 107515419A
- Authority
- CN
- China
- Prior art keywords
- microseismic event
- energy
- working face
- event
- distribution
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 143
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000009826 distribution Methods 0.000 claims abstract description 113
- 230000002902 bimodal effect Effects 0.000 claims abstract description 31
- 230000002123 temporal effect Effects 0.000 claims abstract description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 239000003381 stabilizer Substances 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 description 20
- 238000012544 monitoring process Methods 0.000 description 17
- 238000005065 mining Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 10
- 230000006378 damage Effects 0.000 description 9
- 238000006073 displacement reaction Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 239000003245 coal Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 239000000463 material Substances 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 206010044565 Tremor Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000004080 punching Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Emergency Management (AREA)
- Business, Economics & Management (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Alarm Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of method of estimation and device of rock stability.This method includes：Obtain the temporal information and strength information of microseismic event；According to the temporal information and strength information of the microseismic event, the intensity distribution of the microseismic event is determined；Intensity distribution based on bimodal distribution model and the microseismic event, determine the Evolution of the high-energy microseismic event and low energy microseismic event in the microseismic event；According to the Evolution of the high-energy microseismic event and low energy microseismic event, the rock stability of working face is estimated.Method provided by the invention, the estimation of rock stability is realized, improve the accuracy of estimation, ensure Mine Safety in Production.
Description
Technical field
The present embodiments relate to the method for estimation and device in mining technique field, more particularly to rock stability.
Background technology
During coal mining, with the continuous extension of mine, under the influence of mine pressure and mining pressure,
Rock mass can destroy at any time, cause the loss of personnel and property, therefore the stability of rock mass is estimated, can ensure ore deposit
The safety in production on mountain, reduce the generation of mine accident.
At present, the method for the estimation of the stability of rock mass mainly has：Numerical simulation and displacement monitoring method.Clearly as
The complexity of rock mass materials so that people are difficult thoroughly to recognize rock mass materials this structures, and the anisotropism of rock mass materials, makes in addition
The mechanics parameter for obtaining rock mass materials is difficult to choose, and this stability for all estimating rock mass to numerical simulation brings difficulty.Displacement is supervised
Survey method has important as a kind of main field monitoring method for being used to estimate rock stability to estimation rock stability
Meaning, it is mainly using rock mass before destruction, and based on small deformation, rate of deformation is smaller, the stage Rock Slide Stability；Utilize rock
Body is close to when destroying, and rate of deformation can increased, and the stage rock mass is gradually from stable state to instability status transition；Rock mass
After destruction, rock mass deformation increases rapidly, the stage Instability of Rock Body.It is substantially from the angle of displacement, and research rock mass is from starting to become
Shape to unstability overall process, so as to realize the judgement to rock stability.Due to the complexity of geological conditions, so that rock mass is in
Existing different rock mass deformation rule, even occurs rock mass sometimes under conditions of very little deformation directly from stable state to unstability shape
State changes, and is thus difficult with the displacement monitoring method of rock mass to judge rock stability.
Therefore, when estimating rock stability there is problems with prior art, it is difficult to realize to the rock in mining active process
The accurate estimation of body stability.
The content of the invention
The present invention provides a kind of method of estimation and device of rock stability, to solve accurately estimate in the prior art
The problem of counting the rock stability in mining active process.
On the one hand, the present invention provides a kind of method of estimation of rock stability, and methods described includes：
Obtain the temporal information and strength information of microseismic event；
According to the temporal information and strength information of the microseismic event, the intensity distribution of the microseismic event is determined；
Intensity distribution based on bimodal distribution model and the microseismic event, determines the high-energy in the microseismic event
The Evolution of microseismic event and low energy microseismic event；
According to the Evolution of the high-energy microseismic event and low energy microseismic event, the Rock Slide Stability of working face is estimated
Property.
On the other hand, the present invention provides a kind of estimation unit of rock stability, and described device includes：
Acquisition module, obtain the temporal information and strength information of microseismic event；
Estimation module, according to the temporal information and strength information of the microseismic event, determine the intensity of the microseismic event
Distribution；Intensity distribution based on bimodal distribution model and the microseismic event, determine that the high-energy in the microseismic event is micro-
The Evolution of shake event and low energy microseismic event；According to the evolution of the high-energy microseismic event and low energy microseismic event
Rule, estimate the rock stability of working face.
The method of estimation and device of rock stability provided by the invention, pass through the microseismic event of collecting work face rock mass
Temporal information and strength information, according to the temporal information and strength information of microseismic event, determine the Energy distribution of microseismic event, base
In bimodal distribution model and microseismic event Energy distribution, determine that the evolution of high-energy microseismic event and low energy microseismic event is advised
Rule, according to the Evolution of the high-energy microseismic event and low energy microseismic event, estimate the stability of working face rock mass, it is real
Show the estimation of rock stability, improved the accuracy of estimation, ensure Mine Safety in Production.
Brief description of the drawings
Accompanying drawing herein is merged in specification and forms the part of this specification, shows the implementation for meeting the disclosure
Example, and be used to together with specification to explain the principle of the disclosure.
Fig. 1 is a kind of method of estimation schematic flow sheet of rock stability provided in an embodiment of the present invention；
Fig. 2 is Microseismic monitoring system site layout project schematic diagram provided in an embodiment of the present invention；
Fig. 3 is the method for estimation schematic flow sheet of another rock stability provided in an embodiment of the present invention；
Fig. 4 is that schematic diagram is arranged at coal mine work area underground work station in Dong Jia rivers provided in an embodiment of the present invention；
Fig. 5 is newly-increased microseismic event number-energy profile on April 9 provided in an embodiment of the present invention to October 10；
Fig. 6 is newly-increased microseismic event number-moment magnitude distribution map on April 9 provided in an embodiment of the present invention to October 10；
Fig. 7 moves towards distribution for newly-increased microseismic event number on April 9 provided in an embodiment of the present invention to October 10 along working face
Figure；
Fig. 8 is distributed for newly-increased microseismic event number on April 9 provided in an embodiment of the present invention to October 10 along working depth
Figure；
Fig. 9 is a kind of estimation unit schematic diagram of rock stability provided in an embodiment of the present invention.
Pass through above-mentioned accompanying drawing, it has been shown that the clear and definite embodiment of the disclosure, will hereinafter be described in more detail.These accompanying drawings
It is not intended to limit the scope of disclosure design by any mode with word description, but is by reference to specific embodiment
Those skilled in the art illustrate the concept of the disclosure.
Embodiment
Here exemplary embodiment will be illustrated in detail, its example is illustrated in the accompanying drawings.Following description is related to
During accompanying drawing, unless otherwise indicated, the same numbers in different accompanying drawings represent same or analogous key element.Following exemplary embodiment
Described in embodiment do not represent all embodiments consistent with the disclosure.On the contrary, they be only with it is such as appended
The example of the consistent apparatus and method of some aspects be described in detail in claims, the disclosure.
How to be solved to the technical scheme of technical scheme and the application with specifically embodiment below above-mentioned
Technical problem is described in detail.These specific embodiments can be combined with each other below, for same or analogous concept
Or process may repeat no more in certain embodiments.Below in conjunction with accompanying drawing, embodiments of the invention are described.
In embodiments of the present invention, when rock test piece is in low-stress state, now rock test piece sound emission is with small energy
Based on event, the yardstick of micro- destruction in rock test piece is much smaller than the critical slope langth of rock test piece, and rock test piece is more stable；When
Rock test piece is in high-stress state, and now it is rapid to send out energy acoustie emission event big for rock test piece sound, micro- broken in rock test piece
Bad yardstick is close to or greater than the critical slope langth of rock test piece, rock test piece less stable, close to state of rupture.Therefore,
The stability of rock test piece and big energy events generation rule are closely bound up.
Each vibrations that micro rupture is formed during Coal Exploitation is dynamic are referred to as a microseismic event, microseismic event bimodal distribution
It is due to (low energy microseismic event and high-energy microseism thing caused by the mixing changed at random caused by two kinds of different phenomenons
Part), the low energy microseismic event in the bimodal distribution model of microseismic event is by micro- destruction generation less than critical dimension, high energy
Amount microseismic event is by micro- destruction generation more than critical dimension.The bimodal distribution rule of microseismic event is with loaded load
Increase, a kind of inevitable phenomenon that the small yardstick micro rupture in rock mass inside changes to large scale micro rupture.Therefore bimodal distribution mould is utilized
The Evolution of the big energy microseismic event of type analysis and small energy microseismic event, it is feasible that estimation is carried out to rock stability.
Embodiment one
Fig. 1 is a kind of flow chart of the embodiment 1 of the method for estimation of rock stability provided by the invention, as shown in figure 1,
The method of the present embodiment includes：
S110, the temporal information and strength information for obtaining microseismic event.
In the present embodiment, Microseismic monitoring system is established on working face to be measured, can be spaced very short time constantly monitor it is micro-
Shake event, or, can be with work to be measured in 24 hours uninterrupted continuous monitoring recovery process in order to improve the accuracy of estimation
The microseismic event in face, obtain time, the strength information of microseismic event.
Wherein, working face is the first production scene of coal, is the job site of direct exploitation mineral or rock, with adopting
Driving degree and to depths move.The Microseismic monitoring system being monitored to working face, as shown in Fig. 2 the microseism including underground passes
Sensor, underground work station are data acquisition substation, the data storage server and data processing server of ground matrix section.Microseism passes
Sensor typically can be that quality test of pile is used to gather Wave data, and the quantity of microseismic sensors can be according to specific operating mode
It is configured, such as could be arranged to more than 6, wherein at least 4 microseismic sensors can accurately determines microseismic event
Hypocentral location, remaining sensor can correct the precision of the microseismic event position data of first to fourth sensor acquisition.Cause
This, when a microseismic event occurs, each sensor can collect Wave data in synchronization, be carried out to Wave data
During processing, the Wave data reported of each sensor can be associated according to temporal information, so as to a microseism thing
The spatial positional information that part occurs is confirmed, it is achieved thereby that microseismic event information accurately obtains.
Further, microseismic sensors arrangement basic principle can be：Track lane and return airway are moving towards direction along colliery
It is in interlaced arrangement with depth direction；It is arranged in non-producing side；Microseismic sensors spacing is in 100m, office in track lane or return airway
Portion's position sensor spacing is encrypted to 60m.Microseismic sensors are arranged in hard rock, are close to palisades；Microseismic sensors are in level
Arranged with target area to be monitored staggeredly, is surround in vertical direction.Treated it should be noted that microseismic sensors are typically located at
In the both sides tunnel for detecting working face, working face to be detected is surrounded so as to the positioning precision in working face to be detected
Highest.In addition, microseismic sensors are arranged in hard rock by the way of shallow punching, therefore embodiments of the invention realize non-connect
Touch, the information of microseismic event is obtained at a distance.
Specifically, using the microseismic event in Microseismic monitoring system monitoring face, include gathering the elastic wave of microseismic event
Information, time, the intensity letter of microseismic event are being obtained by the data processing software on the data processing server of earth's surface
Before breath, microseismic event can also be judged, and the elastic wave to collecting is filtered processing.It should be noted that
The energy and moment magnitude of microseismic event are generally smaller, in various embodiments of the present invention, it is necessary to the microseismic event energy of monitoring
It may range from 10-3-10-7Joule, so as to avoid monitoring other energy are larger caused by recovery activity shock event,
Ensure the accuracy of acquisition microseismic event.
Further, time, the strength information of microseismic event are obtained using data processing software.Wherein, microseismic event
Strength information can include microseismic event energy, microseismic event moment magnitude.
S120, temporal information and strength information according to the microseismic event, determine the intensity distribution of the microseismic event.
In the present embodiment, it is that the temporal information of periodic statistics microseismic event and strength information obtain the strong of microseismic event
Degree distribution, specifically, cycle can be one day, one week or one month, and the present invention is not limited herein.
Alternatively, in the present embodiment, it was the cycle with one month, statistics increases microseismic event number and energy newly in one month,
It is determined that increasing microseismic event number-Energy distribution newly in this month, obtain and increase microseismic event number-energy point in this month newly
Butut.
In the present embodiment, measurement period can flexibly be set, therefore when rock mass is in compared with stable state, can be set
In the longer cycle, when rock mass rock mass occurs and unstable state occurred, the shorter cycle is could be arranged to, ensure that rock mass is steady
The real-time of qualitative estimation.
S130, the intensity distribution based on bimodal distribution model and the microseismic event, are determined in the microseismic event
The Evolution of high-energy microseismic event and low energy microseismic event.
In the present embodiment, after the intensity distribution based on bimodal distribution model and the microseismic event, it is also necessary to really
Determine the critical strength of high-energy microseismic event and low energy microseismic event.
Specifically, the bimodal model of microseismic event is due to two kinds of different phenomenons, i.e. high-energy microseismic event and low energy
Caused by the mixing changed at random caused by amount microseismic event, two kinds of peak values of appearance on distribution map are embodied in, together
When there is also low ebb, therefore by finding separation be assured that bimodal critical value.The intensity distribution of microseismic event is divided into
Low energy microseismic event and high-energy microseismic event, in the case where working face is in stable, high-energy microseismic event is less, because
This on newly-increased microseismic event number-energy profile (or newly-increased microseismic event number-moment magnitude distribution map), high energy peak
Newly-increased microseismic event number is less, can be ignored, now distribution map appears as unimodal.But when working face appearance is unstable
During situation, newly-increased high-energy microseismic event number increases, and the situation of bimodal distribution occurs.Therefore, monitoring that rock mass is micro- first
Shake intensity distribution rule meets the cycle of bimodal distribution model, it is possible to determines high-energy microseismic event and low energy microseismic event
Critical strength.
Carried out according to the critical strength can of high-energy microseismic event and low energy microseismic event for microseismic event
Classification, the microseismic event that intensity is more than or equal to critical value are defined as high-energy microseismic event, and intensity is less than the micro- of critical strength
Shake event is defined as low energy microseismic event.
Therefore, according to the critical strength of high-energy microseismic event and low energy microseismic event, it is possible to judge in the cycle
Each microseismic event be high-energy microseismic event or low energy microseismic event, to obtain high-energy microseismic event and low energy
The Evolution of microseismic event.
S140, the Evolution according to the high-energy microseismic event and low energy microseismic event, estimate the rock of working face
Body stability.
Specifically, the Evolution of the high-energy microseismic event and low energy microseismic event can include following feelings substantially
Condition：
The peak value of energy first for increasing microseismic event in cycle newly is higher than the second peak value, and the rock mass of working face, which is in, stablizes shape
State；
The energy variation for increasing microseismic event in cycle newly is less than the second peak value for the first peak value, and working face rock mass is in shakiness
Determine state.
Wherein, the range of energy distribution of the first peak value is less than critical strength.The range of energy distribution of second peak value is more than
Critical strength.
A kind of method of estimation of the rock stability present embodiments provided, it is micro- by three-dimensional arrangement around working face
Shake sensor so that at a distance to the accurate acquisition of microseismic event information；And by flexibly setting monitoring cycle microseism
The real-time of event rock stability estimation, using based on bimodal distribution model, analysis high-energy microseismic event and energy microseism
The Evolution of event, so as to estimate the stability of working face rock mass exactly.
Embodiment two
Fig. 3 is a kind of flow chart of the embodiment 2 of the method for estimation of rock stability provided by the invention, as shown in figure 3,
The method of the present embodiment on the basis of the embodiment 1 further to working face in the stability of each region rock mass estimate,
Therefore methods described also includes：
S150, obtain high-energy microseismic event spatial information.
Specifically, high-energy microseismic event spatial information gathers elastic wave when microseismic event occurs by microseismic sensors
Information, gathered by underground work station to data processing server, using data processing software to micro- on data processing server
The spatial information of shake event is obtained.
Specifically, the spatial information of high-energy microseismic event be high-energy microseismic event three-dimensional location coordinates, the three-dimensional
Marker coordinates are obtained in data processing software by the legal positions of Geiger, and specific method is：By from a given initial point
(testing site) approaches final result by iteration.Iteration each time, it is all based on least square method and calculates a modification vector Δ θ
(Δ x, Δ y, Δ z, Δ t), vectorial Δ θ is added in the result (testing site) of last iteration, obtains a new testing site, so
After judge whether this new testing site meets to require, if meeting to require, this point coordinates is required hypocentral location；It is if discontented
It is sufficient then continue iteration.
Further, after the three-dimensional location coordinates of high-energy microseismic event are determined using Geiger methods, at data
Manage the optimization in software and calculate theoretical method under the existing microseismic sensors arrangement condition of working face, high-energy microseismic event is three-dimensional
Position coordinates error optimizes analysis, to obtain the three-dimensional location coordinates of more accurate high-energy microseismic event.
S160, according to the high-energy microseismic event spatial information, determine the space distribution rule of high-energy microseismic event.
Specifically, after the spatial information of high-energy microseismic event is got, determine to increase high-energy microseism thing in the cycle newly
Part increases the high-energy microseismic event number vertically regularity of distribution newly along the working face horizontal direction regularity of distribution and in the cycle.
Alternatively, according to the spatial information of high-energy microseismic event, determine that increase event number in the cycle newly moves towards along working face
The horizontal distribution rule in direction, obtain microseismic event and move towards (horizontal direction) distribution map along working face；According to high-energy microseism thing
The spatial information of part, microseismic event is obtained along working depth direction (vertical direction) distribution map.
S170, the space distribution rule according to the high-energy microseismic event, estimate the rock in each region in the working face
Body stability.
Specifically, the regularity of distribution and the relation in horizontal 3rd area of the working face according to microseismic event in the horizontal direction, it is determined that
The stability of the upward rock mass of working face reclaimed water square；Three bands are erected with working face according to the regularity of distribution of the microseismic event in vertical direction
Relation, determine the stability of rock mass in vertical direction in working face；According to rock stability in horizontal direction and Vertical Square
Upward rock stability, determine the rock stability in each region in working face.
Wherein, horizontal 3rd area of working face refer to after being produced in coal seam that original original stress equilibrium state is by broken around goaf
Bad, overlying rock moves to form horizontal 3rd area along direction of propulsion, is that rib supports the zone of influence, again absciss layer area, compacting zone respectively, its
The middle rib support zone of influence is supported by wall, vertically moves small, is moved horizontally big；Absciss layer area, under vertical displacement greatly it is upper it is small formed from
Layer；Again compacting zone, it is small big under vertical displacement, it is compacted again.It is violent with absciss layer area strata movement degree in 3rd area, its work
Feature, shield is pushed away on face roof influences the balanced structure very big, area's disrupted bed formation is engaged, and can mitigate work face pressure
Power, balanced structure unstability, working face will bear large-scale pressure effect.
Working face erects three bands and refers to after being produced in coal seam that original original stress equilibrium state is destroyed around goaf, on
Rock stratum moves along aspect vertical direction, produces perpendicular three bands, is caving zone, fissure zone, warp damage, wherein fissure zone respectively
Absciss layer, stretching crack are produced on caving zone, between rock stratum, global configuration is preferable, and the wide displacement of scope is small, general high 15-
70m；Warp damage, absciss layer and fracture do not occur substantially for rock stratum, but large area is slowly sunk.The formation of perpendicular three bands is basic
On successively carry out from bottom to top.
The present embodiment calculates theoretical algorithm using optimization, the spatial information of accurate high-energy microseismic event is obtained, by height
The distribution in the spatial information of energy microseismic event horizontal with working face 3rd area, perpendicular three bands is combined, and is adopted so as to further have estimated
During each region rock mass of working face stability, realize the estimation of specific region rock stability.
Embodiment three
The reality of stability evaluation of rock mass of the embodiment of the present invention three with the working face of Dong Jia rivers colliery 22517 in mining active process
Example illustrates the method for estimation of rock stability：
S110, the temporal information and strength information for obtaining microseismic event.
First, according to the concrete condition of the working face of Dong Jia rivers colliery 22517, the type and sensing of microseismic sensors are determined
The particular location of device arrangement.The face length of Dong Jia rivers colliery 22517 is in 1200m or so, and width is in 180m or so, therefore microseism
Sensor selects geophone, and the response frequency scope of geophone is 15Hz to 1000Hz, sensitivity 43.3V/m/
sec.Underground work station by the elastic wave signal that geophone collects by sending on the data processing server of earth's surface,
Schematic diagram is arranged at the working face underground work station of Dong Jia rivers colliery 22517, as shown in figure 4, the working face of Dong Jia rivers colliery 22517 has 5
Individual underground work station, be respectively positioned at two underground work stations of transportation roadway and positioned at three underground work stations in track lane, its
Two work stations of middle transportation roadway, which are realized, transmits the signal that the geophone in transportation roadway collects to earth's surface service
On device, three work stations in track lane are realized to transmit the signal that the geophone in track lane collects to earth's surface and taken
It is engaged on device；In addition, two underground work stations of transportation roadway are first connected to by optical fiber on three underground work stations in track lane, then
Transmitted a signal to by the 1# underground works station in track lane on ground list server.It should be noted that the arrangement of microseismic sensors
Position needs to combine monitoring purpose and geological condition further determines that.In the present embodiment, geophone is at least 6, its
Middle at least four is used for gathering the elastic wave of microseismic event.
Secondly, using the microseismic event of working face to be measured in the uninterrupted continuous monitoring recovery process of Microseismic monitoring system 24h,
Obtain time, the strength information of each microseismic event.What is directly collected due to geophone is elastic wave information, therefore is adopted
The elastic wave information collected is handled with ESG data processing softwares.For the microseismic event elastic wave information collected
It is filtered；Judge whether the energy of microseismic event is in 10 afterwards-3—10-7In the range of joule, to ensure to obtain microseismic event
Accuracy, however, it is determined that be microseismic event, then the temporal information and intensity of microseismic event are then determined by ESG data processing softwares
Information.
S120, temporal information and strength information according to the microseismic event, determine the intensity distribution of the microseismic event.
In the present embodiment, it was the cycle with one month, the statistics working face of Dong Jia rivers colliery 22517 monthly increases microseismic event number-energy newly
(moment magnitude) is distributed, and draws 2015.4.9-2015.10.10 as shown in Figure 5 and monthly increases microseismic event number-Energy distribution newly
Figure, likewise, drawing 2015.4.9-2015.10.10 as shown in Figure 6 monthly increases microseismic event number newly.
S130, the intensity distribution based on bimodal distribution model and the microseismic event, are determined in the microseismic event
The Evolution of high-energy microseismic event and low energy microseismic event.With 2015.4.9-2015.10.10 shown in Fig. 5 monthly
Exemplified by newly-increased microseismic event number-energy profile, the distribution meets bimodal distribution model, in Figure 5, the sheet that April 9 counted
The moon, there is obvious bimodal distribution in newly-increased microseismic event number-Energy distribution, i.e. the first peak value is higher than the second peak value, its
In the first peak value be range of energy distribution corresponding to low energy microseismic event in 0-300J, the second peak value is high-energy microseism thing
Range of energy distribution corresponding to part is in 300-30000J, it follows that in the present embodiment, critical strength 300J.
It should be noted that in the present embodiment, the microseismic event intensity distributions of 22517 working faces in first month just
Through there is bimodal distribution, it is possible to determine critical strength in first month, in other cases, microseismic event is former
Individual lunar rock body high energy concentrating amount microseismic event frequency is less, mainly based on low energy microseismic event, therefore microseismic event intensity
Distribution is constantly in unimodal state in the past few months, can not now determine critical strength；But elapse over time, work as high energy
Amount microseismic event frequency increases, when bimodal distribution occurs in microseismic event intensity distribution, it is possible to critical strength is determined,
Therefore, when microseismic event intensity distribution bimodal distribution occurs first, it is possible to determine critical strength.
After 2015.4.9 Energy distribution determines the critical strength of 22517 working faces, continue to count
2015.4.9-2015.10.10 monthly increases microseismic event number-Energy distribution newly, it may be determined that 2015.4.9- as shown in Figure 5
2015.10.10 this 7 months microseismic event energy evolution rules.
S140, the Evolution according to the high-energy microseismic event and low energy microseismic event, estimate the rock of working face
Body stability.As shown in figure 5, in the newly-increased microseismic event number-energy profile counted 2015.4.9 days, the first peak value is higher than
Second peak value, illustrate that now rock mass is in stable state；2015.5.8 in the newly-increased microseismic event number-energy profile counted,
It is slightly above the first peak value that second peak value is varied, the newly-increased microseismic event monthly counted from 2015.5.8-2015.10.10
As can be seen that being slightly above the first peak value until the peaks of 2015.10.10 second by the peak values of 2015.5.8 second in number-energy profile
For value more than one times of height of the first peak value, i.e., the energy for monthly increasing microseismic event newly is distributed mainly on second in bimodal distribution model
Corresponding to peak value in energy range, illustrate that unstable state has occurred in rock mass in the process.
S150, obtain high-energy microseismic event spatial information.ESG data processing softwares are using Geiger to microseismic event
Focus is positioned, and final result is approached by iteration especially by from a given initial point (testing site).Each time
Iteration, is all based on least square method and calculates modification vector Δ θ (Δ x, Δ y, Δ z, Δ t), vectorial Δ θ being added to last time and changed
In the result (testing site) in generation, a new testing site is obtained, then judges whether this new testing site meets to require, if full
Foot requires that this point coordinates is required hypocentral location；Continue iteration if being unsatisfactory for.
Further, after the three-dimensional location coordinates of high-energy microseismic event are determined using Geiger methods, at data
Manage the optimization in software and calculate theoretical method under the existing microseismic sensors arrangement condition of working face, high-energy microseismic event is three-dimensional
Position coordinates error optimizes analysis, to obtain the three-dimensional location coordinates of more accurate high-energy microseismic event.
S160, according to the high-energy microseismic event spatial information, determine the space distribution rule of high-energy microseismic event.
High-energy microseismic event number and mining are monthly increased newly in 2015.4.9-2015.10.10 according to the working face of Dong Jia rivers colliery 22517
Progress, determine that 2015.4.9-2015.10.10 monthly increases high-energy microseismic event number newly and moves towards the regularity of distribution along working face, such as scheme
Shown in 7, April 9, statistics was monthly about 780 in working face position corresponding to newly-increased high-energy microseismic event number along work EDS maps
Rice, its microseismic event number peak value illustrate to release energy in the region along 750 meters of working face more；May 8, statistics monthly newly increased
Working face position corresponding to energy microseismic event number along work EDS maps is about 745 meters, and its microseismic event number peak value explanation exists
Released energy along 700 meters of region of working face more；June 4, statistics monthly increased high-energy microseismic event number newly along working face point
Working face position is about 700 meters corresponding to cloth, and its microseismic event number peak value illustrates to discharge energy in the region along 650 meters of working face
Measure more；July 9, statistics was monthly about working face position corresponding to newly-increased high-energy microseismic event number along work EDS maps
675 meters, its microseismic event number peak value illustrates to release energy in the region along 650 meters of working face more；August statistics on the 3rd is monthly new
Working face position corresponding to increasing energy microseismic event number along work EDS maps is about 650 meters, and its microseismic event number peak value is said
The bright region along 650 meters of working face releases energy more；September statistics on the 11st monthly increases high-energy microseismic event number newly along work
Working face position is about 610 meters corresponding to EDS maps, and its microseismic event number peak value illustrates to release along the region of 650 meters of working face
Exoergic amount is more；It is big to count monthly working face position corresponding to newly-increased high-energy microseismic event number along work EDS maps on October 10
General is 575 meters, and its microseismic event number peak value illustrates to release energy in the region along 650 meters of working face more.
Similarly, according to the working face of Dong Jia rivers colliery 22517, in 2015.4.9-2015.10.10, monthly newly-increased high-energy is micro-
Event number and mining progress are shaken, it is deep along working face to determine that 2015.4.9-2015.10.10 monthly increases high-energy microseismic event number newly
The regularity of distribution is spent, as shown in Figure 8.
S170, the space distribution rule according to the high-energy microseismic event, estimate the rock in each region in the working face
Body stability.According to the regularity of distribution of microseismic event in the horizontal direction and the relation in horizontal 3rd area of the working face, working face is determined
The stability of the upward rock mass of reclaimed water square；According to the regularity of distribution of the microseismic event in vertical direction and the pass of perpendicular three bands of working face
System, determine the stability of rock mass in vertical direction in working face；According in rock stability in horizontal direction and vertical direction
Rock stability, determine the rock stability in each region in working face.
Specifically, in the horizontal direction, as shown in fig. 7, high-energy microseism in 2015.4.9-2015.10.10 this seven months
Event number moves towards the regularity of distribution along working face and is distributed mainly on 650 meters of working face monthly to increase releasing energy for microseismic event newly,
In the areas of Ji Heng tri- in front of the support pressure peak of the rib support zone of influence 50 meters to rock stratum absciss layer area.
Further, in vertical direction, as shown in figure 8, high-energy is micro- in 2015.4.9-2015.10.10 this seven months
Shake event number is distributed mainly on working face absolute altitude along the working depth regularity of distribution monthly to increase releasing energy for microseismic event newly
460 meters of position, in addition, the absolute altitude of common working face is 250 meters, therefore monthly increase releasing energy main point for microseismic event newly
It is distributed in the bottom of warp damage in perpendicular three band.
Further, with reference to rock stability in rock stability in horizontal direction and vertical direction, i.e., in Dong Jiahe
On the working face of colliery 22517, horizontal direction is 50 meters in front of the support pressure peak of the rib support zone of influence to rock stratum absciss layer area,
Vertical direction is the bottom of warp damage, and the rock stability in this region is poor, close to state of rupture.
A kind of method of estimation of the rock stability present embodiments provided, it is micro- by three-dimensional arrangement around working face
Shake sensor so that at a distance to the accurate acquisition of microseismic event information；And by flexibly setting monitoring cycle microseism
The real-time of event rock stability estimation, using based on bimodal distribution model, analysis high-energy microseismic event and energy microseism
The Evolution of event, so as to estimate the stability of working face rock mass exactly.And the present embodiment is calculated using optimization
Theoretical algorithm, the spatial information of accurate high-energy microseismic event is obtained, by the spatial information of high-energy microseismic event and work
The distribution in the areas of Mian Heng tri-, perpendicular three bands combines, so as to further have estimated the stabilization of each region rock mass of working face in mining active process
Property, realize the estimation to specific region rock stability.
Example IV
Fig. 9 is a kind of schematic diagram of the estimation unit of rock stability provided by the invention, as shown in figure 9, the present embodiment
Device include：
Acquisition module, obtain the temporal information and strength information of microseismic event；
Estimation module, according to the temporal information and strength information of the microseismic event, determine the intensity of the microseismic event
Distribution；Intensity distribution based on bimodal distribution model and the microseismic event, determine that the high-energy in the microseismic event is micro-
The Evolution of shake event and low energy microseismic event；According to the evolution of the high-energy microseismic event and low energy microseismic event
Rule, estimate the rock stability of working face.
The estimation module is specifically used for：Intensity distribution based on bimodal distribution model and microseismic event, determines high energy
Measure the critical strength of microseismic event and low energy microseismic event；
According to the critical strength, the Evolution of the high-energy microseismic event and low energy microseismic event is determined.
Face stabilizer is estimated according to the Evolution of high-energy microseismic event and low energy microseismic event, if high-energy
The Evolution of microseismic event and low energy microseismic event is：
The first peak value for increasing the energy of microseismic event in cycle newly is higher than the second peak value, and the rock mass of the working face is in steady
Determine state；
The energy variation for increasing microseismic event in cycle newly is less than the second peak value for the first peak value, at the rock mass of the working face
In unstable state；
The range of energy distribution of first peak value is less than the critical strength, the Energy distribution model of second peak value
Enclose and be more than the critical strength.
Alternatively, the acquisition module, it is additionally operable to obtain high-energy microseismic event spatial information；
Alternatively, the estimation module, it is additionally operable to, according to the high-energy microseismic event spatial information, determine that high-energy is micro-
The space distribution rule of shake event；According to the space distribution rule of the high-energy microseismic event, estimate each in the working face
The rock stability in region.
Alternatively, the regularity of distribution and the relation in horizontal 3rd area of the working face according to the microseismic event in the horizontal direction,
Determine the stability of the upward rock mass of the working face reclaimed water square；According to the microseismic event vertical direction the regularity of distribution with
The working face erects the relation of three bands, determines the stability of rock mass in vertical direction in the working face；According to the level side
Rock stability in upward rock stability and the vertical direction, determine the Rock Slide Stability in each region in the working face
Property.
The device of embodiment illustrated in fig. 9 accordingly can be used for performing rock stability in embodiment of the method shown in Fig. 1 and Fig. 2
The step of method of estimation performs, its implementing principle and technical effect is similar, and here is omitted.
A kind of rock stability estimation unit that the present embodiment provides, obtain three-dimensional arrangement microseism around working face and pass
The microseismic event information of sensor collection；Using based on bimodal distribution model, analysis high-energy microseismic event and energy microseismic event
Evolution, so as to estimate the stability of working face rock mass exactly, in addition the present embodiment provide device can obtain
The spatial information of accurate high-energy microseismic event is taken, by the spatial information of high-energy microseismic event and horizontal 3rd area of working face, is erected
The distributions of three bands combines, and so as to further have estimated the stability of each region rock mass of working face in mining active process, realizes pair
The estimation of specific region rock stability.
One of ordinary skill in the art will appreciate that：Realizing all or part of step of above-mentioned each method embodiment can lead to
The related hardware of programmed instruction is crossed to complete.Foregoing program can be stored in a computer read/write memory medium.The journey
Sequence upon execution, execution the step of including above-mentioned each method embodiment；And foregoing storage medium includes：ROM, RAM, magnetic disc or
Person's CD etc. is various can be with the medium of store program codes.
Those skilled in the art will readily occur to the disclosure its after considering specification and putting into practice invention disclosed herein
Its embodiment.It is contemplated that cover any modification, purposes or the adaptations of the disclosure, these modifications, purposes or
Person's adaptations follow the general principle of the disclosure and including the undocumented common knowledges in the art of the disclosure
Or conventional techniques.Description and embodiments are considered only as exemplary, and the true scope of the disclosure and spirit are by following
Claims are pointed out.
It should be appreciated that the precision architecture that the disclosure is not limited to be described above and is shown in the drawings, and
And various modifications and changes can be being carried out without departing from the scope.The scope of the present disclosure is only limited by appended claims
System.
Claims (10)
- A kind of 1. method of estimation of rock stability, it is characterised in that including：Obtain the temporal information and strength information of microseismic event；According to the temporal information and strength information of the microseismic event, the intensity distribution of the microseismic event is determined；Intensity distribution based on bimodal distribution model and the microseismic event, determine the high-energy microseism in the microseismic event The Evolution of event and low energy microseismic event；According to the Evolution of the high-energy microseismic event and low energy microseismic event, the rock stability of working face is estimated.
- 2. according to the method for claim 1, it is characterised in that methods described also includes：Obtain high-energy microseismic event spatial information；According to the high-energy microseismic event spatial information, the space distribution rule of high-energy microseismic event is determined；According to the space distribution rule of the high-energy microseismic event, the rock stability in each region in the working face is estimated.
- 3. according to the method for claim 1, it is characterised in that described to be based on bimodal distribution model and the microseismic event Intensity distribution, determine the Evolution of high-energy microseismic event and low energy microseismic event, specifically include：Intensity distribution based on bimodal distribution model and microseismic event, determines high-energy microseismic event and low energy microseismic event Critical strength；According to the critical strength, the Evolution of the high-energy microseismic event and low energy microseismic event is determined.
- 4. according to the method for claim 3, it is characterised in that described micro- according to the high-energy microseismic event and low energy The Evolution of shake event, estimate the rock stability of working face, specifically include：If the Evolution of the high-energy microseismic event and low energy microseismic event is：The first peak value for increasing the energy of microseismic event in cycle newly is higher than the second peak value, and the rock mass of the working face, which is in, stablizes shape State；The energy variation for increasing microseismic event in cycle newly is less than the second peak value for the first peak value, and the rock mass of the working face is in not Stable state；The range of energy distribution of first peak value is less than the critical strength, and the range of energy distribution of second peak value is big In the critical strength.
- 5. according to the method for claim 2, it is characterised in that advised according to the spatial distribution of the high-energy microseismic event Rule, estimate the rock stability in each region in the working face, including：According to the regularity of distribution of the microseismic event in the horizontal direction and the relation in horizontal 3rd area of the working face, the work is determined The stability of the upward rock mass of face reclaimed water square；The relation of three bands is erected in the regularity of distribution of vertical direction and the working face according to the microseismic event, determines the work In face in vertical direction rock mass stability；According to rock stability in rock stability in the horizontal direction and the vertical direction, determine in the working face The rock stability in each region.
- A kind of 6. estimation unit of rock stability, it is characterised in that including：Acquisition module, for obtaining the temporal information and strength information of microseismic event；Estimation module, for the temporal information and strength information according to the microseismic event, determine the intensity of the microseismic event Distribution；Intensity distribution based on bimodal distribution model and the microseismic event, determine that the high-energy in the microseismic event is micro- The Evolution of shake event and low energy microseismic event；According to the evolution of the high-energy microseismic event and low energy microseismic event Rule, estimate the rock stability of working face.
- 7. device according to claim 6, it is characterised in thatThe acquisition module, it is additionally operable to obtain high-energy microseismic event spatial information；The estimation module, it is additionally operable to, according to the high-energy microseismic event spatial information, determine the sky of high-energy microseismic event Between the regularity of distribution；According to the space distribution rule of the high-energy microseismic event, the rock mass in each region in the working face is estimated Stability.
- 8. device according to claim 6, it is characterised in that the estimation module is specifically used for：Intensity distribution based on bimodal distribution model and microseismic event, determines high-energy microseismic event and low energy microseismic event Critical strength；According to the critical strength, the Evolution of the high-energy microseismic event and low energy microseismic event is determined.
- 9. device according to claim 8, it is characterised in that the estimation module is specifically used for：Face stabilizer is estimated according to the Evolution of high-energy microseismic event and low energy microseismic event, if high-energy microseism The Evolution of event and low energy microseismic event is：The first peak value for increasing the energy of microseismic event in cycle newly is higher than the second peak value, and the rock mass of the working face, which is in, stablizes shape State；The energy variation for increasing microseismic event in cycle newly is less than the second peak value for the first peak value, and the rock mass of the working face is in not Stable state；The range of energy distribution of first peak value is less than the critical strength, and the range of energy distribution of second peak value is big In the critical strength.
- 10. device according to claim 7, it is characterised in that the estimation module is specifically used for：According to the regularity of distribution of the microseismic event in the horizontal direction and the relation in horizontal 3rd area of the working face, the work is determined The stability of the upward rock mass of face reclaimed water square；The relation of three bands is erected in the regularity of distribution of vertical direction and the working face according to the microseismic event, determines the work In face in vertical direction rock mass stability；According to rock stability in rock stability in the horizontal direction and the vertical direction, determine in the working face The rock stability in each region.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710740277.1A CN107515419B (en) | 2017-08-25 | 2017-08-25 | The estimation method and device of rock stability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710740277.1A CN107515419B (en) | 2017-08-25 | 2017-08-25 | The estimation method and device of rock stability |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107515419A true CN107515419A (en) | 2017-12-26 |
| CN107515419B CN107515419B (en) | 2019-02-15 |
Family
ID=60724255
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710740277.1A Active CN107515419B (en) | 2017-08-25 | 2017-08-25 | The estimation method and device of rock stability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107515419B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110910613A (en) * | 2019-12-10 | 2020-03-24 | 大连理工大学 | Rock micro-seismic wireless monitoring, receiving and early warning system |
| CN111323812A (en) * | 2020-04-09 | 2020-06-23 | 南方科技大学 | Method and system for monitoring stability of landslide body based on micro-fracture |
| CN111396128A (en) * | 2020-02-28 | 2020-07-10 | 北京科技大学 | Mining surrounding rock fracture sliding starting condition and instability process analysis method and system |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105807312A (en) * | 2016-03-15 | 2016-07-27 | 大连理工大学 | Coal mine roof rock mass vertical zone determination method based on micro-seismic monitoring |
| CN106054243A (en) * | 2016-06-03 | 2016-10-26 | 河南理工大学 | Rockburst multi-index prediction method based on micro-seismic monitoring |
| WO2017007464A1 (en) * | 2015-07-08 | 2017-01-12 | Halliburton Energy Services, Inc. | Improved fracture matching for completion operations |
| CN106443761A (en) * | 2016-10-09 | 2017-02-22 | 华北科技学院 | Full-frequency band and wide-frequency domain shake monitoring system for mine earthquakes and rock bursts |
| CN106646591A (en) * | 2017-01-17 | 2017-05-10 | 克拉玛依市海晟达石油科技有限公司 | Oil layer fireflood front end micro-seismic monitoring system and method |
| CN106869910A (en) * | 2017-02-08 | 2017-06-20 | 平安煤炭开采工程技术研究院有限责任公司 | Crack recognition methods and system |
-
2017
- 2017-08-25 CN CN201710740277.1A patent/CN107515419B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017007464A1 (en) * | 2015-07-08 | 2017-01-12 | Halliburton Energy Services, Inc. | Improved fracture matching for completion operations |
| CN105807312A (en) * | 2016-03-15 | 2016-07-27 | 大连理工大学 | Coal mine roof rock mass vertical zone determination method based on micro-seismic monitoring |
| CN106054243A (en) * | 2016-06-03 | 2016-10-26 | 河南理工大学 | Rockburst multi-index prediction method based on micro-seismic monitoring |
| CN106443761A (en) * | 2016-10-09 | 2017-02-22 | 华北科技学院 | Full-frequency band and wide-frequency domain shake monitoring system for mine earthquakes and rock bursts |
| CN106646591A (en) * | 2017-01-17 | 2017-05-10 | 克拉玛依市海晟达石油科技有限公司 | Oil layer fireflood front end micro-seismic monitoring system and method |
| CN106869910A (en) * | 2017-02-08 | 2017-06-20 | 平安煤炭开采工程技术研究院有限责任公司 | Crack recognition methods and system |
Non-Patent Citations (2)
| Title |
|---|
| 于群 等: "基于微震监测的锦屏二级水电站深埋隧洞岩爆孕育过程分析", 《岩土工程学报》 * |
| 唐礼忠: "深井矿山地震活动与岩爆监测及预测研究", 《中国博士学位论文全文数据库 工程科技I辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110910613A (en) * | 2019-12-10 | 2020-03-24 | 大连理工大学 | Rock micro-seismic wireless monitoring, receiving and early warning system |
| CN111396128A (en) * | 2020-02-28 | 2020-07-10 | 北京科技大学 | Mining surrounding rock fracture sliding starting condition and instability process analysis method and system |
| CN111323812A (en) * | 2020-04-09 | 2020-06-23 | 南方科技大学 | Method and system for monitoring stability of landslide body based on micro-fracture |
| CN111323812B (en) * | 2020-04-09 | 2022-03-11 | 南方科技大学 | Method and system for monitoring stability of landslide body based on micro-fracture |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107515419B (en) | 2019-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Konicek et al. | Stress changes and seismicity monitoring of hard coal longwall mining in high rockburst risk areas | |
| Xu et al. | Microseismic monitoring and stability evaluation for the large scale underground caverns at the Houziyan hydropower station in Southwest China | |
| Zhao et al. | Study of the rock mass failure process and mechanisms during the transformation from open-pit to underground mining based on microseismic monitoring | |
| Guo et al. | An integrated approach to study of strata behaviour and gas flow dynamics and its application | |
| CN105785471A (en) | Impact danger evaluation method of mine pre-exploiting coal seam | |
| CN105807312A (en) | Coal mine roof rock mass vertical zone determination method based on micro-seismic monitoring | |
| CN107515419A (en) | The method of estimation and device of rock stability | |
| CN110159347A (en) | A kind of deep high stress tight roof stope dynamic disaster monitoring and pre-alarming method | |
| Liu et al. | Trapezoidal collapse model to calculate the height of the overburden collapse zone in coal seam mining: An example from Guo'Jiahe Coal Mine, Western China | |
| CN111695790A (en) | Mining method for security pillar | |
| Deng et al. | Experimental study on the interrelation of multiple mechanical parameters in overburden rock caving process during coal mining in longwall panel | |
| Froese et al. | 25 years of movement monitoring on South Peak, Turtle Mountain: understanding the hazard | |
| Nordlund | Deep hard rock mining and rock mechanics challenges | |
| Bossi et al. | Capabilities of continuous and discontinuous modelling of a complex, structurally controlled landslide | |
| Aydin et al. | Geomechanical characterization, 3-D optical monitoring and numerical modeling in Kirkgecit-1 tunnel, Turkey | |
| CN110133748A (en) | A kind of method of integrated alkalic-metasomatism type uranium deposit deep metallogenetic data | |
| CN105589104A (en) | Method for determining influence of fracture structure on rock burst | |
| Vinoth et al. | Slope stability monitoring by quantification and behavior of microseismic events in an opencast coal mine | |
| Kieffer et al. | Continuous real‐time slope monitoring of the Ingelsberg in Bad Hofgastein, Austria | |
| Ren et al. | Strata movement and fracture evolution characteristics in adjacent seam mining | |
| Dyagilev | Induced seismicity and seismic hazard in mining regions | |
| Kajzar | Geodetic and seismological observations applied for investigation of subsidence formation in the CSM mine (Czech Republic) | |
| CN112035931A (en) | Numerical simulation research method and device for roadway surrounding rock damage under dynamic load effect | |
| Ding et al. | Numerical simulations of how staged dewatering and mining influence surface subsidence | |
| CN110823425B (en) | Construction method and system for constructing nonmetal underwater water-sealed oil depot by using underwater natural rock mass |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 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 |