CN110118991A - A kind of mining induced stress appraisal procedure based on microseism damage reconstruct - Google Patents
A kind of mining induced stress appraisal procedure based on microseism damage reconstruct Download PDFInfo
- Publication number
- CN110118991A CN110118991A CN201910404955.6A CN201910404955A CN110118991A CN 110118991 A CN110118991 A CN 110118991A CN 201910404955 A CN201910404955 A CN 201910404955A CN 110118991 A CN110118991 A CN 110118991A
- Authority
- CN
- China
- Prior art keywords
- induced stress
- grid
- microseism
- grid node
- damage
- 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
- 238000005065 mining Methods 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000009826 distribution Methods 0.000 claims abstract description 35
- 239000003245 coal Substances 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 238000009825 accumulation Methods 0.000 claims abstract description 8
- 239000011435 rock Substances 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 10
- 238000005553 drilling Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 230000035939 shock Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005670 electromagnetic radiation Effects 0.000 description 2
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/25—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons
- G01L1/255—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons using acoustic waves, or acoustic emission
-
- 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
-
- 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/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of mining induced stress appraisal procedures based on microseism damage reconstruct, are suitable for field of mine safety micro seismic monitoring field and use.It specifically includes: grid dividing being carried out to assessment area, and calculates the accumulated deformation energy ε of each grid node using accumulation methodEiWith microseismic event number NiOr the loaded elapsed-time standards Δ t of coal petrographyi, find out the maximum accumulated deformation energy numerical value max { ε of assessment areaEi, and calculating average accumulated deformation can εF, can be obtained mining induced stress distribution on the basis of calculating and obtaining damaging parameter D distribution.The Parameters Calculation that this method physical mechanics meaning is obvious, is related to is clear, universality and strong operability, is suitable for programming and realizes, application feasibility it is good, it can be achieved that in the mining process of coal seam mining induced stress approximate real time inverting.
Description
Technical field
The present invention relates to a kind of mine mining induced stress appraisal procedures, are particularly suitable for the base in mine safety micro seismic monitoring field
In the mining induced stress appraisal procedure of microseism damage reconstruct.
Background technique
Mining induced stress is the one kind of primitive stress after by the redistribution of digging disturbing influence in the country rock body of down-hole mining space
This life stress.The distribution characteristics for understanding the stress is working face design, stops adopting line design, Coal Pillar Design, safety and Protection, branch
Indispensable one of important evidence and work in the mining production process such as shield design.
Since digging disturbance in coal seam is a dynamic process, also the moment generates dynamic change to mining induced stress, finally leads
This "black box" process is caused to be difficult to directly a wide range of observation and quantitative description.The method for being currently used primarily in mining induced stress observation
Have: drilling hole stress monitoring method, electromagnetic radiation method, drilling cuttings method, micro-seismic method, shock wave topographic imaging by velocity method etc..Wherein, drilling is answered
Power monitoring method is a kind of method of direct observation, but its observation for being limited only to tunnel coal wall superficial part small range region；Electromagnetic radiation method
It is indirect assessment method with drilling cuttings method, is a kind of indirect relation progress according between electromagnetic parameter, coal powder quantity of bore and mining induced stress
The method of assessment, and such method also cannot achieve a wide range of observation, while be influenced more obviously by ambient noise；It is micro-
Shake method induces the tool of coal petrography micro rupture event as a kind of strong effective monitoring digging disturbance, supervises in mine safety at present
Survey field is widely applied, and the shock wave topographic imaging by velocity technology of especially its extension can realize a wide range of detection, together
When be based on microseismic event spatial distribution, microseism frequency Density Distribution and energy density distribution, can qualitative evaluation mining induced stress indirectly
Coverage.However, shock wave topographic imaging by velocity technology is to be distributed to calculate mining induced stress distribution indirectly according to velocity of longitudinal wave,
And due to needing a certain number of microseismic events as inverting initial data in its calculating process, thus inevitably
Cause to lead to certain error calculated because of the unreasonable hypothesis that velocity of longitudinal wave must be set to constant during inverting, simultaneously should
Technology calculation amount is generally bigger, is difficult to realize real time inversion；Based on microseismic event spatial distribution, microseism frequency Density Distribution with
The appraisal procedure of energy density distribution, although can infer mining induced stress point by the distribution of reflection mining induced fissure to a certain extent
Cloth, but it lacks apparent physical mechanics association.Therefore, gone out based on micro seismic monitoring data reconstruction a kind of with physical mechanics meaning
The method of justice and the approximate real time a wide range of assessment mining induced stress of energy, has very important practical value and realistic meaning.
Summary of the invention
Goal of the invention: shortcoming in view of the above technology, provide it is a kind of based on microseism damage reconstruct mining induced stress comment
Estimate method, with specific reference to the micro seismic monitoring data obtained in real time, synchronous calculate obtains mining induced stress distribution, realizes coal seam digging
The approximate real time inverting of mining induced stress in journey.
Technical solution: to achieve the above object, the mining induced stress appraisal procedure of the invention based on microseism damage reconstruct is first
The Damage Parameter of stope is first reconstructed according to microseism parameter；It is then based on damage mechanics and stope stress is obtained by association Damage Parameter
Distribution, and then mining induced stress field distribution is obtained,
Specific step is as follows:
Assessment area is carried out grid dividing and forms grid dividing figure by a, with the corresponding statistics circle of each grid node for region
Statistical window calculates the accumulated deformation energy ε of each statistical regions using accumulation methodEiWith microseismic event number NiOr the loaded experience of coal petrography
Time Δ tiNumerical value as each grid node；
B traverses the ordered series of numbers accumulated deformation energy ε in grid dividing figureEiFind its maximum value max { εEi, and calculate assessment
The average accumulated deformation in region can εE；
C utilizes formula:Calculate corresponding damage at each grid node in grid dividing figure
Parameter Di, D in formulaiFor the corresponding Damage Parameter numerical value of i-th of grid node；
D calculates the mining induced stress numerical value at each grid node:
When working face mining speed approximation is at the uniform velocity stable, strain-time mode is preferentially used:
σi=E·αt·Δti·(1-Di)
When working face mining speed is unstable, approximation uses strain-microseism frequency mode:
σi=E·αN·Ni·(1-Di)
In formula: σiFor mining induced stress numerical value corresponding at i-th of grid node；E is elasticity modulus；αtFor strain-time system
Number；αNFor strain-microseism frequency coefficient, interpolation finally is carried out to the mining induced stress numerical value at each grid node, can be obtained and comment
Estimate the mining induced stress space distribution information in region, finally obtain the stress envelope of tested region using distributed intelligence, is mine
Foundation is instructed in safe design offer.
It with s is grid dividing spacing in grid dividing figure, r is statistics sliding radius, to avoid in statistics slipping
It omits individual microseismic events and result is caused to be distorted, it is as follows that the two meets relationship:Specific calculating process are as follows: with each net
The corresponding statistics circle of lattice node is range statistics window, and the accumulated deformation energy ε of each statistical regions is calculated using accumulation methodEiWith it is micro-
Shake event number NiOr the loaded elapsed-time standards Δ t of coal petrographyiAs the numerical value of each grid node, its calculation formula is:
Δti=tiN-ti1
In formula: εEiIndicate the corresponding accumulated deformation energy for counting circle region of i-th of grid node；NiIndicate i-th of grid
The corresponding microseismic event number for counting circle region of node；EijIndicate that corresponding j-th for counting circle region of i-th of grid node is micro-
The energy of shake event；ΔtiIndicate the corresponding loaded elapsed-time standards for counting circle region of i-th of grid node；tiNIndicate i-th of net
Lattice node is corresponding to count the time that the last one microseismic event occurs in circle region；ti1Indicate the corresponding statistics of i-th of grid node
The time that first microseismic event occurs in circle region.
Average accumulated deformation can εFCalculation formula are as follows:In formula: max { εEiIt is that assessment area is maximum
Accumulated deformation energy numerical value；DcFor Damage Parameter numerical value corresponding under complete faulted condition, 0.95 is chosen here.
The utility model has the advantages that the parameter that mining induced stress calculation formula physical mechanics meaning of the present invention is obvious, formula is related to
Clear, universality and strong operability are calculated, is suitable for programming and realizes, application feasibility is good；Number is surveyed involved in calculating process
According to the microseism data using a wide range of real-time monitoring of mine, timeliness is high, can a wide range of approximate real time inverting coal seam mining process
In mining induced stress distribution, while can realize daily monitoring and warning, there is very important practical value and realistic meaning.
Detailed description of the invention
Fig. 1 is that the present invention is based on the mining induced stress distribution schematic diagrams of the mining induced stress appraisal procedure of microseism damage reconstruct；
Fig. 2 is that the present invention is based on the grid dividing schematic diagrames of the mining induced stress appraisal procedure of microseism damage reconstruct；
Fig. 3 is microseismic event spatial distribution map；
Fig. 4 is the accumulated deformation energy spatial distribution map based on microseism Parameters Calculation；
Fig. 5 is the loaded elapsed-time standards spatial distribution map of coal petrography based on microseism Parameters Calculation；
Fig. 6 is the Damage Parameter distribution map based on microseism Parameters Calculation；
Fig. 7 is the mining induced stress distribution map based on microseism damage reconstruct；
Specific embodiment
Further description is made to the present invention with reference to the accompanying drawing.
With down-hole coal bed exploitation, mining induced stress as shown in Figure 1 will be formed in coal and rock in front of working face and is distributed
(ABCD), longitudinal overlying strata including elastic region (AB), plastic zone (BC), post-peak softening area (CD), are respectively corresponded spatially
Warp damage, fissure zone, caving zone.From damage mechanics angle analysis, the coal petrography material of D point to mined out region is to reach
Complete faulted condition, corresponding accumulation microseismic event distribution density also can reach maximum in this region；
Mining induced stress appraisal procedure based on microseism damage reconstruct of the invention, step are as follows: first according to microseism parameter weight
The Damage Parameter of structure stope；It is then based on damage mechanics and stope stress distribution is obtained by association Damage Parameter, and then adopted
Dynamical field distribution；
Specific steps are as follows:
A. assessment area progress grid dividing is formed into grid dividing figure as shown in Figure 2, s is between grid dividing in figure
It is statistics sliding radius away from, r, causes result to be distorted to avoid omitting individual microseismic events in statistics slipping, the two is full
Sufficient relationship is as follows:Its specific calculating process are as follows: with the corresponding statistics circle of each grid node for range statistics window, adopt
The accumulated deformation energy ε of each statistical regions is calculated with accumulation methodEIWith microseismic event number NIOr the loaded elapsed-time standards Δ t of coal petrographyIMake
For the numerical value of each grid node, wherein εeiFor calculating Damage Parameter numerical value (see step c), NiWith Δ tiIt is respectively used to calculate
Under strain-microseism frequency mode and strain-time mode mining induced stress numerical value (see step d), its calculation formula is:
Δti=tiN-ti1
In formula: εEiIndicate the corresponding accumulated deformation energy for counting circle region of i-th of grid node；NiIndicate i-th of grid
The corresponding microseismic event number for counting circle region of node；EijIndicate that corresponding j-th for counting circle region of i-th of grid node is micro-
The energy of shake event；ΔtiIndicate the corresponding loaded elapsed-time standards for counting circle region of i-th of grid node；tiNIndicate i-th of net
Lattice node is corresponding to count the time that the last one microseismic event occurs in circle region；ti1Indicate the corresponding statistics of i-th of grid node
The time that first microseismic event occurs in circle region.
B. the ordered series of numbers accumulated deformation energy ε in grid dividing figure is traversedEiFind its maximum value max { εEi, and calculate assessment
The average accumulated deformation in region can εF；In formula: max { εEiIt is the maximum accumulated deformation energy numerical value of assessment area；DcIt is complete
Corresponding Damage Parameter numerical value, chooses 0.95 here under faulted condition；
C. corresponding damaging parameter D at each grid node is calculatedi:
In formula: DiFor the corresponding Damage Parameter numerical value of i-th of grid node；
D. the mining induced stress numerical value at each grid node is calculated:
When working face mining speed approximation is at the uniform velocity stable, strain-time mode is preferentially used:
σi=E αt·Δti·(1-Di)
When working face mining speed is unstable, approximation uses strain-microseism frequency mode:
σi=E αN·Ni·(1-Di)
In formula: σiFor mining induced stress numerical value corresponding at i-th of grid node；E is elasticity modulus；αtFor strain-time system
Number；αNFor strain-microseism frequency coefficient.Interpolation finally is carried out to the mining induced stress numerical value at each grid node, can be obtained and comment
Estimate the mining induced stress spatial distribution in region.
Embodiment one:
The micro seismic monitoring data that certain coal mine work area recovery phase is chosen in instance analysis are analyzed, since the working face is opened
It is at the uniform velocity stable to adopt speed approximation, average daily footage is 1.2m, and therefore, present invention calculating is finally opened up by taking strain-time mode as an example
Open explanation.Implement the present invention according to inventive concept:
(1) as shown in Figure 3 according to microseismic event spatial distribution；Three-dimensional grid division is carried out to assessment area, is taken between grid
It is 10m away from s, statistics sliding radius r is 30m, and calculates the accumulated deformation energy ε at each grid node using accumulation methodEiWith
The loaded elapsed-time standards Δ t of coal petrographyi, interpolation computing method is then used, you can get it such as Fig. 4 accumulated deformation energy spatial distribution and figure
Accumulated deformation energy spatial distribution map shown in 5；
Fig. 5 is the loaded elapsed-time standards spatial distribution map of coal petrography based on microseism Parameters Calculation
(2) grid node sequence ε is traversedEi, obtaining cumulative maximum deformation can max { εEiNumerical value be 9141.3678, according to
This calculates average accumulated deformation can εFNumerical value be 3051.4635.
(3) by εEiAnd εFSubstitute into formulaThe Damage Parameter numerical value at each grid node is obtained,
Then show that Damage Parameter distribution is as shown in Figure 6 by interpolation calculation.
(4) by actual elastic modulus E=9GPa and factor alphat=0.000026 substitutes into σi=E αt·Δti·(1-Di),
The mining induced stress numerical value at each grid node is obtained, then show that mining induced stress distribution is as shown in Figure 7 by interpolation calculation.
Example shows that clear Parameters Calculation of the present invention, universality and strong operability, Simultaneous Inversion are calculated
Mining induced stress be reasonably distributed, effect preferably, it can be achieved that in the mining process of coal seam mining induced stress approximate real time inverting.
Claims (3)
1. a kind of mining induced stress appraisal procedure based on microseism damage reconstruct, it is characterised in that adopted first according to the reconstruct of microseism parameter
The Damage Parameter of field；It is then based on damage mechanics and stope stress distribution is obtained by association Damage Parameter, and then obtain to adopt and answer
Force distribution, the specific steps are as follows:
Assessment area is carried out grid dividing and forms grid dividing figure by a, with the corresponding statistics circle of each grid node for range statistics
Window calculates the accumulated deformation energy ε of each statistical regions using accumulation methodEiWith microseismic event number NiOr the loaded elapsed-time standards of coal petrography
ΔtiNumerical value as each grid node；
B traverses the ordered series of numbers accumulated deformation energy ε in grid dividing figureEiFind its maximum value max { εEi, and calculate assessment area
Average accumulated deformation can εF；
C utilizes formula:Corresponding Damage Parameter at each grid node in calculating grid dividing figure
Di, D in formulaiFor the corresponding Damage Parameter numerical value of i-th of grid node；
D calculates the mining induced stress numerical value at each grid node:
When working face mining speed approximation is at the uniform velocity stable, strain-time mode is preferentially used:
σi=E αt·Δti·(1-Di)
When working face mining speed is unstable, approximation uses strain-microseism frequency mode:
σi=E αN·Ni·(1-Di)
In formula: σiFor mining induced stress numerical value corresponding at i-th of grid node；E is elasticity modulus；αtFor strain-time coefficient；
αNFor strain-microseism frequency coefficient, interpolation finally is carried out to the mining induced stress numerical value at each grid node, can be obtained assessment area
The mining induced stress space distribution information in domain finally obtains the stress envelope of tested region using distributed intelligence, is mine safety
Design provides and instructs foundation.
2. the mining induced stress appraisal procedure according to claim 1 based on microseism damage reconstruct, it is characterised in that: in grid
It with s is grid dividing spacing in division figure, r is statistics sliding radius, to avoid omitting individual microseism things in statistics slipping
Part and cause result to be distorted, it is as follows that the two meets relationship:Specific calculating process are as follows: with the corresponding system of each grid node
Meter circle is range statistics window, and the accumulated deformation energy ε of each statistical regions is calculated using accumulation methodEiWith microseismic event number NiOr coal
The loaded elapsed-time standards Δ t of rockiAs the numerical value of each grid node, its calculation formula is:
Δti=tiN-ti1
In formula: εEiIndicate the corresponding accumulated deformation energy for counting circle region of i-th of grid node；NiIndicate i-th of grid node
The corresponding microseismic event number for counting circle region；EiiIndicate corresponding j-th of the microseism thing for counting circle region of i-th of grid node
The energy of part；ΔtiIndicate the corresponding loaded elapsed-time standards for counting circle region of i-th of grid node；tiNIndicate i-th of grid section
Point is corresponding to count the time that the last one microseismic event occurs in circle region；ti1It indicates that i-th of grid node is corresponding and counts circle area
The time that first microseismic event occurs in domain.
3. the mining induced stress appraisal procedure according to claim 1 based on microseism damage reconstruct, it is characterised in that: average tired
Product deformation can εFCalculation formula are as follows:In formula: max { εEiIt is the maximum accumulated deformation energy number of assessment area
Value；DcFor Damage Parameter numerical value corresponding under complete faulted condition, 0.95 is chosen here.
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910404955.6A CN110118991B (en) | 2019-05-16 | 2019-05-16 | Mining induced stress assessment method based on microseismic damage reconstruction |
| PCT/CN2020/088156 WO2020228546A1 (en) | 2019-05-16 | 2020-04-30 | Mining-induced stress assessment method based on microseismic damage reconstruction |
| AU2020275806A AU2020275806B2 (en) | 2019-05-16 | 2020-04-30 | Mining-induced stress assessment method based on microseismic damage reconstruction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910404955.6A CN110118991B (en) | 2019-05-16 | 2019-05-16 | Mining induced stress assessment method based on microseismic damage reconstruction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110118991A true CN110118991A (en) | 2019-08-13 |
| CN110118991B CN110118991B (en) | 2020-06-23 |
Family
ID=67522497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910404955.6A Active CN110118991B (en) | 2019-05-16 | 2019-05-16 | Mining induced stress assessment method based on microseismic damage reconstruction |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN110118991B (en) |
| AU (1) | AU2020275806B2 (en) |
| WO (1) | WO2020228546A1 (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110244354A (en) * | 2019-07-11 | 2019-09-17 | 东北大学 | A kind of metal mine mining disturbance stress field quantifies dynamic playback method |
| WO2020228546A1 (en) * | 2019-05-16 | 2020-11-19 | 中国矿业大学 | Mining-induced stress assessment method based on microseismic damage reconstruction |
| CN112377257A (en) * | 2020-10-26 | 2021-02-19 | 中国矿业大学 | Working face mining advance influence range determining method based on microseismic monitoring |
| CN114935513A (en) * | 2022-07-26 | 2022-08-23 | 四川中水成勘院工程物探检测有限公司 | Method for predicting generation and expansion of concrete dam body crack based on microseismic signal characteristics |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112324506B (en) * | 2020-11-20 | 2024-05-14 | 上海大屯能源股份有限公司江苏分公司 | Dynamic early warning method for preventing rock burst of coal mine based on microseism |
| CN113312742B (en) * | 2021-03-15 | 2023-10-10 | 中国再保险(集团)股份有限公司 | Annular space grid data structure and construction and retrieval method and device thereof |
| CN113482720A (en) * | 2021-07-02 | 2021-10-08 | 中煤第三建设（集团）有限责任公司 | Mechanized operation line construction process under rock burst condition |
| CN114810211B (en) * | 2021-11-26 | 2023-03-10 | 中国矿业大学 | Rock burst danger prediction method based on mine seismic group shock wave energy attenuation characteristics |
| CN115014613B (en) * | 2022-06-28 | 2023-05-23 | 中国科学院武汉岩土力学研究所 | Monitoring method for surrounding rock stress and deformation of coal mine tunnel |
| CN115826037B (en) * | 2022-11-02 | 2023-12-12 | 北京国信安科技术有限公司 | Evaluation method for monitoring and early warning capability of microseism monitoring and early warning system in mine goaf |
| CN116540299B (en) * | 2023-07-05 | 2023-09-26 | 煤炭科学研究总院有限公司 | Early warning method based on microseismic energy accumulation tendency for coal mine scene |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7460436B2 (en) * | 2005-12-05 | 2008-12-02 | The Board Of Trustees Of The Leland Stanford Junior University | Apparatus and method for hydraulic fracture imaging by joint inversion of deformation and seismicity |
| CN103306722A (en) * | 2013-06-21 | 2013-09-18 | 中国矿业大学 | Detection evaluation method for microearthquake multi-dimensional information integration area of impact danger zone |
| CN104331533A (en) * | 2014-09-18 | 2015-02-04 | 山东科技大学 | Poor-sealing drill hole mining response value simulation method |
| CN105422170A (en) * | 2015-11-02 | 2016-03-23 | 安徽理工大学 | Grouting reinforcement treatment method for mid-depth goaf below building foundation |
| CN105607127A (en) * | 2016-01-27 | 2016-05-25 | 中国矿业大学 | Microquake multi-parameter early warning method of rock burst in high-stress concentration area |
| US9389326B2 (en) * | 2011-03-23 | 2016-07-12 | Global Ambient Seismic, Inc. | Methods, systems and devices for near-well fracture monitoring using tomographic fracture imaging techniques |
| CN106285782A (en) * | 2016-08-30 | 2017-01-04 | 中国矿业大学(北京) | Bump method for early warning under a kind of complicated geological Environmental effect and system |
| CN106873029A (en) * | 2017-01-19 | 2017-06-20 | 秦福亮 | A kind of determination method on coal and gas outburst index and its critical condition |
| US10013800B1 (en) * | 2009-10-23 | 2018-07-03 | Emerson Paradigm Holding Llc | Systems and methods for coordinated editing of seismic data in dual model |
| CN109063257A (en) * | 2018-07-02 | 2018-12-21 | 山东科技大学 | A kind of coal and rock subregion water filling seepage flow-damage-stress coupling method for numerical simulation |
| CN109447837A (en) * | 2018-11-15 | 2019-03-08 | 国家能源投资集团有限责任公司 | A kind of risk assessment method of mining area bump |
| CN109581492A (en) * | 2017-09-29 | 2019-04-05 | 中国石油化工股份有限公司 | Petrophysical parameter calculation method and system based on Simulating Seismic Wave |
| CN109597124A (en) * | 2018-12-29 | 2019-04-09 | 煤炭科学技术研究院有限公司 | A kind of fine detection method of working face stress concentration region channel wave seismic |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8902707B2 (en) * | 2007-04-09 | 2014-12-02 | Baker Hughes Incorporated | Analysis of uncertainty of hypocenter location using the combination of a VSP and a subsurface array |
| US9176245B2 (en) * | 2009-11-25 | 2015-11-03 | Halliburton Energy Services, Inc. | Refining information on subterranean fractures |
| CA2945467C (en) * | 2014-05-23 | 2019-12-17 | Halliburton Energy Services, Inc. | Enhancing reservoir characterization using real-time srv and fracture evolution parameters |
| CN104100297B (en) * | 2014-07-31 | 2017-01-25 | 煤炭科学技术研究院有限公司 | Self-vibrating type micro-quake monitoring system and self-vibrating type micro-quake monitoring method |
| CN105891874B (en) * | 2016-06-30 | 2017-12-12 | 大连理工大学 | One kind adopts coal and rock gushing water micro seismic monitoring method |
| CN106646607B (en) * | 2016-12-22 | 2018-11-27 | 中国矿业大学 | A kind of adaptive unequal spacing Meshing Method improving CT resolution of inversion and efficiency |
| CN110118991B (en) * | 2019-05-16 | 2020-06-23 | 中国矿业大学 | Mining induced stress assessment method based on microseismic damage reconstruction |
-
2019
- 2019-05-16 CN CN201910404955.6A patent/CN110118991B/en active Active
-
2020
- 2020-04-30 WO PCT/CN2020/088156 patent/WO2020228546A1/en active Application Filing
- 2020-04-30 AU AU2020275806A patent/AU2020275806B2/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7460436B2 (en) * | 2005-12-05 | 2008-12-02 | The Board Of Trustees Of The Leland Stanford Junior University | Apparatus and method for hydraulic fracture imaging by joint inversion of deformation and seismicity |
| US10013800B1 (en) * | 2009-10-23 | 2018-07-03 | Emerson Paradigm Holding Llc | Systems and methods for coordinated editing of seismic data in dual model |
| US9389326B2 (en) * | 2011-03-23 | 2016-07-12 | Global Ambient Seismic, Inc. | Methods, systems and devices for near-well fracture monitoring using tomographic fracture imaging techniques |
| CN103306722A (en) * | 2013-06-21 | 2013-09-18 | 中国矿业大学 | Detection evaluation method for microearthquake multi-dimensional information integration area of impact danger zone |
| CN104331533A (en) * | 2014-09-18 | 2015-02-04 | 山东科技大学 | Poor-sealing drill hole mining response value simulation method |
| CN105422170A (en) * | 2015-11-02 | 2016-03-23 | 安徽理工大学 | Grouting reinforcement treatment method for mid-depth goaf below building foundation |
| CN105607127A (en) * | 2016-01-27 | 2016-05-25 | 中国矿业大学 | Microquake multi-parameter early warning method of rock burst in high-stress concentration area |
| CN106285782A (en) * | 2016-08-30 | 2017-01-04 | 中国矿业大学(北京) | Bump method for early warning under a kind of complicated geological Environmental effect and system |
| CN106873029A (en) * | 2017-01-19 | 2017-06-20 | 秦福亮 | A kind of determination method on coal and gas outburst index and its critical condition |
| CN109581492A (en) * | 2017-09-29 | 2019-04-05 | 中国石油化工股份有限公司 | Petrophysical parameter calculation method and system based on Simulating Seismic Wave |
| CN109063257A (en) * | 2018-07-02 | 2018-12-21 | 山东科技大学 | A kind of coal and rock subregion water filling seepage flow-damage-stress coupling method for numerical simulation |
| CN109447837A (en) * | 2018-11-15 | 2019-03-08 | 国家能源投资集团有限责任公司 | A kind of risk assessment method of mining area bump |
| CN109597124A (en) * | 2018-12-29 | 2019-04-09 | 煤炭科学技术研究院有限公司 | A kind of fine detection method of working face stress concentration region channel wave seismic |
Non-Patent Citations (1)
| Title |
|---|
| 窦林名 等: "煤矿冲击矿压动静载的"应力场-震动波场"监测预警技术", 《岩石力学与工程学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020228546A1 (en) * | 2019-05-16 | 2020-11-19 | 中国矿业大学 | Mining-induced stress assessment method based on microseismic damage reconstruction |
| CN110244354A (en) * | 2019-07-11 | 2019-09-17 | 东北大学 | A kind of metal mine mining disturbance stress field quantifies dynamic playback method |
| CN112377257A (en) * | 2020-10-26 | 2021-02-19 | 中国矿业大学 | Working face mining advance influence range determining method based on microseismic monitoring |
| CN114935513A (en) * | 2022-07-26 | 2022-08-23 | 四川中水成勘院工程物探检测有限公司 | Method for predicting generation and expansion of concrete dam body crack based on microseismic signal characteristics |
Also Published As
| Publication number | Publication date |
|---|---|
| AU2020275806B2 (en) | 2022-11-24 |
| WO2020228546A1 (en) | 2020-11-19 |
| AU2020275806A1 (en) | 2021-02-04 |
| CN110118991B (en) | 2020-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110118991A (en) | A kind of mining induced stress appraisal procedure based on microseism damage reconstruct | |
| Zhou et al. | Seepage channel development in the crown pillar: Insights from induced microseismicity | |
| Xu et al. | Microseismic monitoring and stability analysis of the left bank slope in Jinping first stage hydropower station in southwestern China | |
| 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. | The analysis of rock damage process based on the microseismic monitoring and numerical simulations | |
| Xu et al. | The dynamic evaluation of rock slope stability considering the effects of microseismic damage | |
| Hosseini | Evaluation of the rockburst potential in longwall coal mining using passive seismic velocity tomography and image subtraction technique | |
| Feng et al. | Sectional velocity model for microseismic source location in tunnels | |
| Zhang et al. | The role of seismic triggering in a deep-seated mudstone landslide, China: historical reconstruction and mechanism analysis | |
| Xu et al. | Excavation-induced microseismicity: microseismic monitoring and numerical simulation | |
| Hosseini et al. | Passive seismic velocity tomography on longwall mining panel based on simultaneous iterative reconstructive technique (SIRT) | |
| He et al. | Integrated rockburst early warning model based on fuzzy comprehensive evaluation method | |
| CN106156432A (en) | Deformation of tunnel based on three-dimensional geological information assessment and support design method and system | |
| Jiang et al. | A novel method for automatic identification of rock fracture signals in microseismic monitoring | |
| Jahandideh et al. | Inference of rock flow and mechanical properties from injection-induced microseismic events during geologic CO2 storage | |
| CN112901158A (en) | Hydraulic fracture length prediction method and fracture network modeling method and device | |
| Xiao et al. | Stability analysis of surrounding rock mass in underground powerhouse considering damage effect of microfractures | |
| Maxwell et al. | A comparison between controlled source and passive source seismic velocity images | |
| Mreyen et al. | Dynamic numerical modelling of co-seismic landslides using the 3D distinct element method: Insights from the Balta rockslide (Romania) | |
| Viegas et al. | Mapping cave front growth utilising the collective behaviour of seismicity and velocity fields | |
| Zhang et al. | Improvement of microseismic source location during cavern excavation in faulted rock mass using fast marching method | |
| Molka | Tomographic imaging associated with a Mw 2.6 fault-slip event in a Deep Nickel Mine | |
| CN109239777A (en) | A method of it is developed using joint inversion method detection deformation coal | |
| Qian et al. | Microseismic activity characteristics and range evaluation of hydraulic fracturing in coal seam | |
| Kerr | Applications of double-difference tomography for a deep hard rock mine |
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 |