CN110118991B - Mining induced stress assessment method based on microseismic damage reconstruction - Google Patents
Mining induced stress assessment method based on microseismic damage reconstruction Download PDFInfo
- Publication number
- CN110118991B CN110118991B CN201910404955.6A CN201910404955A CN110118991B CN 110118991 B CN110118991 B CN 110118991B CN 201910404955 A CN201910404955 A CN 201910404955A CN 110118991 B CN110118991 B CN 110118991B
- Authority
- CN
- China
- Prior art keywords
- microseismic
- mining
- grid node
- statistical
- 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.)
- Active
Links
Images
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. analysis, for interpretation, for correction
- 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. analysis, for interpretation, for correction
- 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
Abstract
The invention relates to a mining stress assessment method based on microseismic damage reconstruction, which is suitable for the field of mine safety and microseismic monitoring. The method specifically comprises the following steps: carrying out grid division on the evaluation area, and calculating the accumulated deformation energy epsilon of each grid node by adopting an accumulation methodEiAnd number of microseismic events NiOr the coal rock loading elapsed time delta tiFinding the maximum value of the accumulated deformation energy max [ epsilon ] in the evaluation areaEiAnd calculating the average accumulated deformation energy epsilonFAnd acquiring mining stress distribution on the basis of acquiring damage parameter D distribution by calculation. The method has obvious physical and mechanical significance, definite related parameter calculation, strong universality and operability, is suitable for programming realization, has good application feasibility, and can realize approximate real-time inversion of mining stress in the coal seam mining process.
Description
Technical Field
The invention relates to a mining induced stress assessment method, in particular to a mining induced stress assessment method based on microseismic damage reconstruction, which is suitable for the field of mine safety microseismic monitoring.
Background
The mining induced stress is one of the induced stresses in the surrounding rock mass of the underground mining space after the original stress is redistributed under the influence of mining disturbance. The understanding of the distribution characteristics of the stress is one of the indispensable important basis and work in the mining production processes such as working face design, stope line design, coal pillar design, safety protection design, support design and the like.
Because coal seam mining disturbance is a dynamic process, mining stress is dynamically changed all the time, and finally the black box process is difficult to directly observe and quantitatively describe in a large range. At present, methods mainly used for observing mining induced stress comprise: a drilling stress monitoring method, an electromagnetic radiation method, a drilling cutting method, a microseismic method, a vibration wave velocity tomography method and the like. The drilling stress monitoring method is a direct observation method, but is only limited to the observation of a small-range area of a shallow part of a roadway coal wall; the electromagnetic radiation method and the drilling cutting method are indirect evaluation methods, and are evaluation methods according to indirect relations among electromagnetic parameters, drilling cutting quantity and mining stress, and the methods cannot realize large-scale observation and are obviously influenced by surrounding environment noise; the microseism method is used as a powerful tool for effectively monitoring the coal rock micro-fracture event induced by mining disturbance, is widely applied to the field of mine safety monitoring at present, particularly the expanded vibration wave velocity tomography technology can realize large-range detection, and meanwhile, the influence range of mining stress can be indirectly and qualitatively evaluated on the basis of the space distribution, the microseism frequency density distribution and the energy density distribution of the microseism event. However, the seismic wave velocity tomography technology indirectly calculates mining stress distribution according to longitudinal wave velocity distribution, and in the calculation process, a certain number of microseismic events are required to be used as inversion original data, so that a certain calculation result error is inevitably caused due to the unreasonable assumption that the longitudinal wave velocity must be set as a constant in the inversion period, and meanwhile, the calculated amount of the technology is generally large, and real-time inversion is difficult to realize; although the mining stress distribution can be inferred to a certain extent by reflecting the mining fracture distribution based on the evaluation method of the microseismic event space distribution, the microseismic frequency density distribution and the energy density distribution, the method lacks obvious physical and mechanical association. Therefore, the method which has physical and mechanical significance and can approximate real-time large-range estimation of the mining stress is reconstructed based on the microseismic monitoring data, and has very important practical value and practical significance.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the defects of the technology, the mining stress assessment method based on the microseismic damage reconstruction is provided, and specifically, the mining stress distribution is synchronously calculated according to microseismic monitoring data acquired in real time, so that the approximate real-time inversion of the mining stress in the coal seam mining process is realized.
The technical scheme is as follows: in order to achieve the purpose, the mining induced stress assessment method based on microseismic damage reconstruction comprises the following steps of reconstructing damage parameters of a stope according to microseismic parameters; then acquiring stope stress distribution by correlating damage parameters based on damage mechanics, further acquiring mining induced stress field distribution,
the method comprises the following specific steps:
a, carrying out grid division on the evaluation area to form a grid division graph, taking a statistical circle corresponding to each grid node as an area statistical window, and calculating the accumulated deformation energy epsilon of each statistical area by adopting an accumulation methodEiAnd number of microseismic events NiOr the coal rock loading elapsed time delta tiAs the value of each grid node;
b traversing the series of cumulative deformation energies ε in the mesh-divided graphEiFind its maximum value max [ epsilon ]EiAnd calculating the average accumulated deformation energy epsilon of the evaluation areaE;
c using the formula:calculating corresponding damage parameters D at each grid node in grid division diagramiIn the formula DiThe damage parameter value corresponding to the ith grid node is obtained;
d, calculating the mining stress value at each grid node:
when the mining speed of the working face is approximately stable at a constant speed, a strain-time mode is preferentially adopted:
σi=E·αt·Δti·(1-Di)
when the mining speed of the working face is unstable, a strain-microseismic frequency mode is approximately adopted:
σi=E·αN·Ni·(1-Di)
in the formula: sigmaiIs the corresponding mining stress value at the ith grid node, E is the elastic modulus, αtIs the strain-time coefficient αNFinally, interpolation is carried out on the mining stress values at each grid node for strain-microseismic frequency coefficient, so as to obtain the mining stress spatial distribution information of the evaluation area,and finally, obtaining a stress distribution map of the measured area by using the distribution information, and providing a guidance basis for mine safety design.
In the gridding chart, s is a gridding interval, r is a statistical slip radius, and in order to avoid the result distortion caused by missing individual microseismic events in the statistical slip process, the two satisfy the following relation:the specific calculation process is as follows: taking the statistical circle corresponding to each grid node as a regional statistical window, and calculating the accumulated deformation energy epsilon of each statistical region by adopting an accumulation methodEiAnd number of microseismic events NiOr the coal rock loading elapsed time delta tiAs the numerical value of each grid node, the calculation formula is as follows:
Δti=tiN-ti1
in the formula: epsilonEiRepresenting the accumulated deformation energy of the statistical circle region corresponding to the ith grid node; n is a radical ofiRepresenting the number of microseismic events of the statistical circle region corresponding to the ith grid node; eijRepresenting the energy of the jth microseismic event of the ith grid node corresponding to the statistical circle region; Δ tiRepresenting the loading elapsed time of the statistical circle region corresponding to the ith grid node; t is tiNRepresenting the time when the ith grid node corresponds to the last microseismic event in the statistical circle region; t is ti1Indicating the time when the ith grid node corresponds to the occurrence of the first microseismic event in the statistical circle region.
Average cumulative deformation energy εFThe calculation formula of (2) is as follows:in the formula: max [ epsilon ]EiThe maximum accumulated deformation energy value of the evaluation area is obtained; dcFor the corresponding value of the damage parameter in the completely damaged state, 0.95 is selected here.
Has the advantages that: the mining induced stress calculation formula has obvious physical and mechanical significance, definite parameter calculation related to the formula, strong universality and operability, suitability for programming realization and good application feasibility; the actual measurement data related to the calculation process adopts microseismic data monitored in real time in a large range in a mine, has high timeliness, can approximately invert mining stress distribution in the coal seam mining process in real time in a large range, can realize daily monitoring and early warning, and has very important practical value and practical significance.
Drawings
FIG. 1 is a schematic view of mining stress distribution of a mining stress evaluation method based on microseismic damage reconstruction according to the present invention;
FIG. 2 is a schematic diagram of grid division of the microseismic damage reconstruction-based mining induced stress assessment method of the present invention;
FIG. 3 is a microseismic event spatial distribution plot;
FIG. 4 is a diagram of a spatial distribution of accumulated deformation energy based on microseismic parameter calculations;
FIG. 5 is a time-space distribution diagram of coal rock loading experience calculated based on microseismic parameters;
FIG. 6 is a distribution diagram of damage parameters calculated based on microseismic parameters;
FIG. 7 is a distribution diagram of mining stress based on microseismic damage reconstruction;
Detailed Description
The present invention will be further described with reference to the accompanying drawings.
With the mining of the underground coal bed, a mining induced stress distribution (ABCD) as shown in figure 1 is formed in the coal rock body in front of the working face, and the ABCD comprises an elastic region (AB), a plastic region (BC) and a post-peak strain softening region (CD) which respectively correspond to a bending subsidence zone, a fracture zone and a caving zone on a longitudinal overlying strata space. Analyzing from the angle of damage mechanics, the coal rock materials from the point D to the goaf area all reach a complete damage state, and the distribution density of the corresponding accumulated microseismic events can also reach the maximum in the area;
the invention discloses a microseismic damage reconstruction-based mining stress evaluation method, which comprises the following steps: firstly, reconstructing damage parameters of a stope according to the microseismic parameters; then acquiring stope stress distribution by correlating damage parameters based on damage mechanics, and further acquiring mining induced stress field distribution;
the method comprises the following specific steps:
a. the evaluation area is subjected to grid division to form a grid division diagram as shown in fig. 2, wherein s is a grid division distance, r is a statistical slip radius, and in order to avoid the result distortion caused by missing individual microseismic events in the statistical slip process, the two satisfy the following relation:the specific calculation process is as follows: taking the statistical circle corresponding to each grid node as a regional statistical window, and calculating the accumulated deformation energy epsilon of each statistical region by adopting an accumulation methodEIAnd number of microseismic events NIOr the coal rock loading elapsed time delta tIAs a value of each mesh node, where εeiFor calculating damage parameter values (see step c), NiAnd Δ tiRespectively used for calculating the mining stress value under the strain-microseismic frequency mode and the strain-time mode (see step d), and the calculation formula is as follows:
Δti=tiN-ti1
in the formula: epsilonEiRepresenting the accumulated deformation energy of the statistical circle region corresponding to the ith grid node; n is a radical ofiRepresenting the number of microseismic events of the statistical circle region corresponding to the ith grid node; eijRepresenting the energy of the jth microseismic event of the ith grid node corresponding to the statistical circle region; Δ tiRepresenting the loading elapsed time of the statistical circle region corresponding to the ith grid node; t is tiNRepresenting the time when the ith grid node corresponds to the last microseismic event in the statistical circle region; t is ti1Indicating the time when the ith grid node corresponds to the occurrence of the first microseismic event in the statistical circle region.
b. Cumulative deformation energy epsilon of number series in traversal grid division graphEiFind its maximum value max [ epsilon ]EiAnd calculating the average accumulated deformation energy epsilon of the evaluation areaF(ii) a In the formula: max [ epsilon ]EiThe maximum accumulated deformation energy value of the evaluation area is obtained; dcSelecting 0.95 as the corresponding damage parameter value in the complete damage state;
c. calculating corresponding damage parameters D at each grid nodei:
In the formula: diThe damage parameter value corresponding to the ith grid node is obtained;
d. calculating the mining stress value at each grid node:
when the mining speed of the working face is approximately stable at a constant speed, a strain-time mode is preferentially adopted:
σi=E·αt·Δti·(1-Di)
when the mining speed of the working face is unstable, a strain-microseismic frequency mode is approximately adopted:
σi=E·αN·Ni·(1-Di)
in the formula: sigmaiIs the corresponding mining stress value at the ith grid node, E is the elastic modulus, αtIs the strain-time coefficient αNIs the strain-microseismic frequency coefficient. And finally, interpolating the mining stress values at the nodes of each grid to obtain the mining stress spatial distribution of the evaluation area.
The first embodiment is as follows:
the microseismic monitoring data of a certain coal mine working face extraction stage is selected for analysis, and the mining speed of the working face is approximately constant and stable, and the average daily footage is 1.2m, so that the calculation of the method is finally explained by taking a strain-time mode as an example. The invention is implemented according to the inventive idea:
(1) the spatial distribution according to microseismic events is shown in figure 3; three-dimensional grid division is carried out on the evaluation area, the grid spacing s is 10m, the statistical slip radius r is 30m, and the accumulated deformation energy epsilon of each grid node is calculated by adopting an accumulation methodEiLoaded with coal rockLapse of time Δ tiThen, by adopting an interpolation calculation method, the spatial distribution of the accumulated deformation energy shown in fig. 4 and the spatial distribution map of the accumulated deformation energy shown in fig. 5 can be obtained;
FIG. 5 is a coal rock loading experience time-space distribution diagram calculated based on microseismic parameters
(2) Traverse the grid node sequence epsilonEiObtaining the maximum accumulated deformation energy max [ epsilon ]EiThe value of 9141.3678, the average cumulative deformation energy ε is calculatedFThe value of (d) is 3051.4635.
(3) Will epsilonEiAnd εFSubstitution formulaThe values of the damage parameters at the nodes of each grid are obtained, and then the distribution of the damage parameters is obtained through interpolation calculation as shown in fig. 6.
(4) The actual modulus of elasticity E is set to 9GPa and the coefficient αtSubstitution of 0.000026 into σi=E·αt·Δti·(1-Di) Acquiring mining stress values at the nodes of each grid, and calculating the mining stress distribution by interpolation as shown in FIG. 7.
The example shows that the parameter calculation related to the invention is clear, the universality and the operability are strong, the mining stress obtained by inversion calculation is reasonable in distribution and good in effect, and the approximate real-time inversion of the mining stress in the coal seam mining process can be realized.
Claims (3)
1. A mining induced stress assessment method based on microseism damage reconstruction is characterized in that firstly, damage parameters of a stope are reconstructed according to microseism parameters; and then acquiring stope stress distribution by correlating damage parameters based on damage mechanics, and further acquiring mining induced stress field distribution, wherein the method comprises the following specific steps of:
a, carrying out grid division on the evaluation area to form a grid division graph, taking a statistical circle corresponding to each grid node as an area statistical window, and calculating the accumulated deformation energy epsilon of each statistical area by adopting an accumulation methodEiAnd number of microseismic events NiOr the coal rock loading elapsed time delta tiAs a number of nodes of each meshA value;
b traversing the series of cumulative deformation energies ε in the mesh-divided graphEiFind its maximum value max [ epsilon ]EiAnd calculating the average accumulated deformation energy epsilon of the evaluation areaF;
c using the formula:calculating corresponding damage parameters D at each grid node in grid division diagramiIn the formula DiThe damage parameter value corresponding to the ith grid node is obtained;
d, calculating the mining stress value at each grid node:
when the mining speed of the working face is approximately stable at a constant speed, a strain-time mode is adopted:
σi=E·αt·Δti·(1-Di)
when the mining speed of the working face is unstable, adopting a strain-microseismic frequency mode:
σi=E·αN·Ni·(1-Di)
in the formula: sigmaiIs the corresponding mining stress value at the ith grid node, E is the elastic modulus, αtIs the strain-time coefficient αNAnd finally, interpolating the mining stress numerical value at each grid node for the strain-microseismic frequency coefficient to obtain the mining stress spatial distribution information of the evaluation area, and finally obtaining the stress distribution map of the measured area by using the distribution information to provide a guidance basis for mine safety design.
2. The microseismic damage reconstruction-based mining induced stress assessment method according to claim 1, characterized in that: in the gridding chart, s is a gridding interval, r is a statistical slip radius, and in order to avoid the result distortion caused by missing individual microseismic events in the statistical slip process, the two satisfy the following relation:the specific calculation process is as follows: to be provided withTaking the statistical circle corresponding to each grid node as a regional statistical window, and calculating the accumulated deformation energy epsilon of each statistical region by adopting an accumulation methodEiAnd number of microseismic events NiOr the coal rock loading elapsed time delta tiAs the numerical value of each grid node, the calculation formula is as follows:
Δti=tiN-ti1
in the formula: epsilonEiRepresenting the accumulated deformation energy of the statistical circle region corresponding to the ith grid node; n is a radical ofiRepresenting the number of microseismic events of the statistical circle region corresponding to the ith grid node; eijRepresenting the energy of the jth microseismic event of the ith grid node corresponding to the statistical circle region; Δ tiRepresenting the loading elapsed time of the statistical circle region corresponding to the ith grid node; t is tiNRepresenting the time when the ith grid node corresponds to the last microseismic event in the statistical circle region; t is ti1Indicating the time when the ith grid node corresponds to the occurrence of the first microseismic event in the statistical circle region.
3. The microseismic damage reconstruction-based mining induced stress assessment method according to claim 1, characterized in that: average cumulative deformation energy εFThe calculation formula of (2) is as follows:in the formula: max [ epsilon ]EiThe maximum accumulated deformation energy value of the evaluation area is obtained; dcFor the corresponding value of the damage parameter in the completely damaged state, 0.95 is selected 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 CN110118991A (en) | 2019-08-13 |
CN110118991B true 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) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110118991B (en) * | 2019-05-16 | 2020-06-23 | 中国矿业大学 | 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 |
CN112377257B (en) * | 2020-10-26 | 2021-08-20 | 中国矿业大学 | Working face mining advance influence range determining method based on microseismic monitoring |
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 |
CN114935513A (en) * | 2022-07-26 | 2022-08-23 | 四川中水成勘院工程物探检测有限公司 | Method for predicting generation and expansion of concrete dam body crack based on microseismic signal characteristics |
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 (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
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 |
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 (14)
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 |
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 |
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 |
CN103306722B (en) * | 2013-06-21 | 2015-08-05 | 中国矿业大学 | A kind of microseism multidimensional information general area Detection and assessment method of danger of burst district |
US10302791B2 (en) * | 2014-05-23 | 2019-05-28 | 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 |
CN105607127B (en) * | 2016-01-27 | 2018-01-02 | 中国矿业大学 | A kind of microseism many reference amounts method for early warning of high stress concentration zone bump |
CN105891874B (en) * | 2016-06-30 | 2017-12-12 | 大连理工大学 | One kind adopts coal and rock gushing water micro seismic monitoring method |
CN106285782B (en) * | 2016-08-30 | 2018-11-13 | 中国矿业大学(北京) | A kind of bump method for early warning under complicated geological environmental effect and system |
CN106646607B (en) * | 2016-12-22 | 2018-11-27 | 中国矿业大学 | A kind of adaptive unequal spacing Meshing Method improving CT resolution of inversion and efficiency |
CN106873029B (en) * | 2017-01-19 | 2020-02-07 | 秦福亮 | Method for determining coal and gas outburst indexes and critical states thereof |
CN109447837B (en) * | 2018-11-15 | 2022-05-20 | 国家能源投资集团有限责任公司 | Danger assessment method for rock burst in mining area |
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 (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10013800B1 (en) * | 2009-10-23 | 2018-07-03 | Emerson Paradigm Holding Llc | Systems and methods for coordinated editing of seismic data in dual model |
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 |
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 |
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 |
---|
煤矿冲击矿压动静载的"应力场-震动波场"监测预警技术;窦林名 等;《岩石力学与工程学报》;20170430;第36卷(第4期);第803-811页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110118991A (en) | 2019-08-13 |
WO2020228546A1 (en) | 2020-11-19 |
AU2020275806B2 (en) | 2022-11-24 |
AU2020275806A1 (en) | 2021-02-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110118991B (en) | Mining induced stress assessment method based on microseismic damage reconstruction | |
CN110109895B (en) | Surrounding rock grading combined prediction method suitable for TBM tunneling tunnel and application | |
Xiao et al. | ISRM suggested method for in situ microseismic monitoring of the fracturing process in rock masses | |
Dai et al. | Microseismic early warning of surrounding rock mass deformation in the underground powerhouse of the Houziyan hydropower station, China | |
Hosseini | Evaluation of the rockburst potential in longwall coal mining using passive seismic velocity tomography and image subtraction technique | |
CN114810211B (en) | Rock burst danger prediction method based on mine seismic group shock wave energy attenuation characteristics | |
Xu et al. | Comprehensive evaluation of the stability of the left-bank slope at the Baihetan hydropower station in southwest China | |
CN113390458B (en) | Method for judging damage degree of surrounding rock in blasting area | |
He et al. | Integrated rockburst early warning model based on fuzzy comprehensive evaluation method | |
Li et al. | Characteristics of microseismic b-value associated with rock mass large deformation in underground powerhouse caverns at different stress levels | |
Jiang et al. | A novel method for automatic identification of rock fracture signals in microseismic monitoring | |
CN114779330B (en) | Mining working face main fracture azimuth analysis and prediction method based on microseismic monitoring | |
Pettitt et al. | The assessment of damage around critical engineering structures using induced seismicity and ultrasonic techniques | |
CN110376643B (en) | Micro-seismic effect data processing method for jet grouting pile diameter detection | |
Feng et al. | Microseismic energy and intensity criterion of rockburst in deep TBM tunnels: A case study of the Neelum-Jhelum hydropower project | |
Mooney et al. | Influence of geological conditions on measured TBM vibration frequency | |
Bakeev et al. | Modeling of karst deformation and analysis of acoustic emission during sinkhole formation | |
Gangrade et al. | Investigating seismicity surrounding an excavation boundary in a highly stressed dipping underground limestone mine | |
Hardy Jr et al. | Study of microseismic activity associated with a longwall coal mining operation using a near-surface array | |
Wang et al. | Experimental study on stress evolution and microseismic signals under vibration conditions of coal during excavation and subsequent waiting time | |
ZHANG et al. | P-wave velocity tomography and acoustic emission characteristics of sandstone under uniaxial compression | |
Sivakumar et al. | Real time microseismic monitoring to study geomechanics of underground structures | |
CN117233840B (en) | Underground cavern full-space advanced geological prediction method based on seismic waves | |
CN114895359B (en) | DAS (data acquisition system) same-well monitoring real-time microseism effective event denoising method and system | |
Hidayat et al. | Analysis of Hypocenter Relocation by Using Double Difference Method in a Deep Underground Mining |
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 |