CN107450098B - A kind of latent gushing water karst collapse col umn dynamic positioning method of seat earth - Google Patents
A kind of latent gushing water karst collapse col umn dynamic positioning method of seat earth Download PDFInfo
- Publication number
- CN107450098B CN107450098B CN201710657566.5A CN201710657566A CN107450098B CN 107450098 B CN107450098 B CN 107450098B CN 201710657566 A CN201710657566 A CN 201710657566A CN 107450098 B CN107450098 B CN 107450098B
- Authority
- CN
- China
- Prior art keywords
- value
- point
- monitoring
- trend
- positioning
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 41
- 238000012544 monitoring process Methods 0.000 claims abstract description 51
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 238000012937 correction Methods 0.000 claims description 9
- 230000003068 static effect Effects 0.000 claims description 9
- 238000005211 surface analysis Methods 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 claims description 7
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 230000000295 complement effect Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 239000011435 rock Substances 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000010287 polarization Effects 0.000 claims description 2
- 238000002360 preparation method Methods 0.000 claims 2
- 230000000644 propagated effect Effects 0.000 claims 1
- 238000005070 sampling Methods 0.000 claims 1
- 230000004807 localization Effects 0.000 abstract description 5
- 238000003384 imaging method Methods 0.000 abstract description 3
- 239000003245 coal Substances 0.000 description 7
- 238000005065 mining Methods 0.000 description 2
- 238000006424 Flood reaction Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000009172 bursting Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 239000011440 grout Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1429—Subsurface, e.g. in borehole or below weathering layer or mud line
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to a kind of localization methods, belong to technical field of geophysical exploration, and in particular to a kind of latent gushing water karst collapse col umn localization method of seat earth.This method acquires " secondary vibration " that seat earth hidden karst pillar is generated as high artesian conduit pipe by ground On Microseismic Monitoring Technique and is used as focus, pass through shortest path imaging, back analysis, so that it is determined that EXIT POINT that underground is lain concealed, dynamic, and then it is accurately positioned out hidden karst pillar position.
Description
Technical field
The present invention relates to a kind of localization methods, belong to technical field of geophysical exploration, and in particular to a kind of seat earth
The dynamic positioning method of latent gushing water karst collapse col umn.
Background technique
In recent years, watered-out well accident of dashing forward caused by the high pressure-bearing Ordovician karst water of seat earth presents some new features, i.e., difficult to understand
The top of grey Genesis of Karst Subsided Column is developed to 30~50m below seat earth, but and do not enter coal seam, therefore working face of coal seam is visited
Also the karst collapse col umn of small scale is difficult to find when all going well when survey, carrying out grout transformation to seat earth;And after the back production of coal seam,
The zone of fracture of seat earth, fissure zone are gradually linked up with potential, latent karst collapse col umn lateral, vertical, so that Gao Cheng
The Ordovician karst water of pressure is lagged from goaf and is discharged by bottom plate " two bands ", is ultimately caused and is flooded well, floods face even personnel casualty accidents
Occur.
Karst collapse col umn is a kind of very special isolated geologic body, and small with cross section, longitudinal extension is long, form is not advised
Then, the features such as filler difference is big, conventional ground probing means are difficult to control.Around the geophysical exploration about exhausted of karst collapse col umn, coal
The high-resolution three-dimension earthquake of ore mining area has played important function.Previous a large amount of three dimensional seismic data prospecting-mining ratios the result shows that:Ground
The accuracy rate that the 3-D seismics karst collapse col umn above for major diameter 30m in face is detected is about 50% or so, other small-sized are subside
The discovery degree of column is lower.In recent years, underground coal mine seam seismic exploration technology is close with its detection range, influence factor is few, divides
The unique advantages such as resolution height, can be realized to karst collapse col umn latent, major diameter 10m or more inside underground coal mine working face can
By detecting, and it is then helpless in the following karst collapse col umn of seat earth for lying concealed.
Therefore, after seat earth hidden karst pillar causes mine water inrush to flood well, how gushing water point is quickly and accurately determined
It is set to as the key technical problem of mine water damage accident emergency rescue.
Summary of the invention
Present invention is generally directed to the prior arts to be difficult to the technical issues of detecting to small-sized karst collapse col umn and hidden karst pillar, provides
A kind of seat earth lies concealed the dynamic positioning method of gushing water karst collapse col umn.This method acquires coal seam by ground On Microseismic Monitoring Technique
Bottom plate hidden karst pillar is used as focus as " the secondary vibration " that high artesian conduit pipe generates, be imaged by shortest path,
Back analysis, so that it is determined that EXIT POINT that underground is lain concealed, dynamic, and then it is accurately positioned out hidden karst pillar position.
Above-mentioned technical problem of the invention is mainly to be addressed by following technical proposals:
(1) arrange that optimal micro seismic monitoring monitors system.The planar range of the Microseismic monitoring system of surface deployment " can wrap
Enclose " EXIT POINT, and there are certain time differences for the first arrival time of different micro seismic monitoring points.
(2) it is arranged on each monitoring point reception device (wave detector);
(3) Microseismic monitoring system needs continuous uninterrupted recording;
(4) microseismic event is extracted inside collected record;
(5) coordinate (x in monitoring point is extracted from each microseismic eventi,yi,zi,ti) (i=1,2,3 ..., n);
(6) static correction is carried out to each monitoring point coordinate, is corrected to all receiving points all on a datum level, it is specific to walk
Suddenly it is:
Wherein Δ t is positive correcting value;ziTo monitor point height;z0For benchmark face elevation;V is velocity of wave.
(7) focus shortest path imaging method is used, gushing water range is finely quickly determined on time isopleth plan view.
The above method the specific steps are:
1) under the assumed condition based on underground uniform dielectric model
2) the first arrival time plan view for the same microseismic event that different micro seismic monitoring points receive is drawn;
3) residual error extraction is carried out to microseismic event first arrival time isogram using trend surface analysis technology, found out in plane
The position of time shortest point.
Trend surface analysis technology residual error extract the specific steps are:
1. the trend surface fitting method of first arrival time isopleth
If focus is to the first arrival time t of monitoring pointiIt is by its Trend value(ambient field), exceptional value uiWith interference value si(with
Machine noise) composition, i.e.,
Wherein, tiMonitor value when to walk,Trend value when to walk, uiExceptional value when to walk, siFor random noise.
Because random noise is suppressed well when microseismic event is extracted, above formula can be approximately:
If there is n group monitor value (x on groundi,yi,zi,ti) (i=1,2,3 ..., n), wherein (xi,yi,zi) it is that monitoring point is sat
Mark, tiMonitor value when to walk, due to data carry out static correction, so use P order polynomial it is carried out trend surface fitting method for:
Wherein, P is the number of trend surface;b0,b1,b2,...,bL-1For undetermined coefficient;L is the number of fitting coefficient, L=
(P+1)(P+2)/2。
By n to monitoring point (xi,yi) and it is corresponding walk when monitor value tiBring formula (4) into, it is former according to least square method valuation
Reason enables
According to extreme value theorem, to make Q reach minimum, large linear systems is resolved to (5) derivation, the L of (4) can be found out
A coefficient value, finds out fitting surface.
2. degree of fitting calculates
If ti,Monitor value and Trend value when to be walked on i-th of monitoring point,For its average value, total sum of squares of deviations
Wherein, Q2 numerical value is bigger, then fitting degree is higher.
3. significance test
The effect of trend surface fitting method can carry out complementary inspection using F distribution:
In formula, L is the number of multi-trend equation coefficient;N is monitoring points;The freedom degree of Q1 is n-L;Q2 from
By spending for L-1.
Under normal circumstances, significant property coefficient is divided into 0.05 and 0.01 two rank, for given insolation level α, can there is F
(L-1, n-L) distribution table finds critical value Fα;As F > FαWhen, then it is assumed that the situation of change of trend surface reacting dose is significant, it is on the contrary then
It is not significant.
4. residual error is extracted
After microseismic event carries out trend surface analysis, it is 0 that the abnormal Traveltime data separated, which has mathematical expectation, side
The smallest statistical property of difference, integrally tends to normal distribution.This is that the quantitative criterion that travel time anomaly threshold value divides can be determined as:
For abnormal point, time most short point as on first arrival time plan view.
This point is the monitoring point nearest apart from underground EXIT POINT floor projection.
(8) at the time of generation by the available rock mass micro rupture of inversion method, position and property, and according to micro rupture
Size, intensity, failure density, so definitely under latent projective water point spatial position.
The inversion method of seismic source location much has at present:Three circle intersection positioning modes, longitudinal and shear wave time difference method, the homotype wave time difference
The methods of method, polarization analysis positioning mode, Geiger seismic source positioning method, Geiger revised law, joint inversion.It can be according to different situations
Different inversion methods is selected to determine the spatial position of latent projective water point.
Therefore, the invention has the advantages that:It can be accurately positioned seat earth hidden karst pillar below, drawn a circle to approve rapidly
The range of gushing water karst collapse col umn, the position for determining projective water point.
Detailed description of the invention
Fig. 1 seat earth lies concealed gushing water karst collapse col umn localization method schematic diagram;
Fig. 2 T mine W6 water gush face point microseism positioning figure.
Specific embodiment
Below with reference to the embodiments and with reference to the accompanying drawing the technical solutions of the present invention will be further described.
Embodiment:
Prominent watered-out well occurs for T mine, is Ordovician karst water by hydrologic monitoring hole preliminary judgement water bursting source, thus it is speculated that water inrush channel is
Karst collapse col umn.When the high pressure-bearing Ordovician karst water gushing water of seat earth, projective water point can impact back and form " secondary focus ", using ground
On Microseismic Monitoring Technique when being nearby unfolded long to gushing water karst collapse col umn, dynamic micro seismic monitoring and quickly positions, it is determined that doubtful prominent at 2
Water spot.By drilling verification, the karst collapse col umn position that microseism determines differs 26m with probing control point, projective water point for seat earth with
Lower 77m, small-sized karst collapse col umn is lain concealed for 19m × 7m.
The present embodiment proposes a kind of localization method of latent gushing water karst collapse col umn of seat earth.Specifically include following steps:
(1) arrange that optimal micro seismic monitoring monitors system.In order to rapidly determine projective water point position, it is desirable that surface deployment it is micro-
The planar range of shake monitoring system " can surround " EXIT POINT, and the first arrival time of different micro seismic monitoring points there are it is certain when
Difference.First in the sparse laying micro seismic monitoring point of larger range, " coarse positioning " is carried out to latent EXIT POINT;Then, according to first circle
The microseismic event that fixed working face nearby occurs nearby carries out " fine positioning " to working face by encrypting micro seismic monitoring point;Most
Afterwards, on the basis of preceding positioning twice, micro seismic monitoring point is further encrypted, realizes " fine positioning " of latent projective water point;
(2) it is arranged on each monitoring point reception device (wave detector);
(3) Microseismic monitoring system needs continuous uninterrupted recording;
(4) microseismic event is extracted inside collected record;
(5) coordinate (x in monitoring point is extracted from each microseismic eventi,yi,zi,ti) (i=1,2,3 ..., n);
(6) static correction is carried out to each monitoring point coordinate, is corrected to all receiving points all on a datum level, it is specific to walk
Suddenly it is:
Wherein Δ t is positive correcting value;ziTo monitor point height;z0For benchmark face elevation;V is velocity of wave.
(7) focus shortest path imaging method, fine quickly positioning are used.
The above method the specific steps are:
1) it selects to may be assumed that the subsurface model for uniform dielectric;
2) the first arrival time plan view for the same microseismic event that different micro seismic monitoring points receive is drawn;
3) residual error extraction is carried out to microseismic event first arrival time isogram using trend surface analysis technology, found out in plane
The position of time shortest point.
Trend surface analysis technology residual error extract the specific steps are:
1. the trend surface fitting method of first arrival time isopleth
If focus is to the first arrival time t of monitoring pointiIt is by its Trend value(ambient field), exceptional value uiWith interference value si(with
Machine noise) composition, i.e.,
Wherein, tiMonitor value when to walk,Trend value when to walk, uiExceptional value when to walk, siFor random noise.
Because random noise is suppressed well when microseismic event is extracted, above formula can be approximately:
If there is n group monitor value (x on groundi,yi,zi,ti) (i=1,2,3 ..., n), wherein (xi,yi,zi) it is that monitoring point is sat
Mark, tiMonitor value when to walk, due to data carry out static correction, so use P order polynomial it is carried out trend surface fitting method for:
Wherein, P is the number of trend surface;b0,b1,b2,...,bL-1For undetermined coefficient;L is the number of fitting coefficient, L=
(P+1)(P+2)/2;
By n to monitoring point (xi,yi) and it is corresponding walk when monitor value tiBring formula (4) into, it is former according to least square method valuation
Reason enables
According to extreme value theorem, to make Q reach minimum, large linear systems is resolved to (5) derivation, the L of (4) can be found out
A coefficient value, finds out fitting surface.
2. degree of fitting calculates
If ti,Monitor value and Trend value when to be walked on i-th of monitoring point,For its average value, then total sum of squares of deviations
Q is:
Wherein, Q1 is residual sum of squares (RSS), and the more big then fitting degree of value is lower, and Q2 is regression sum of square, is worth bigger fitting journey
It spends higher.
3. significance test
The effect of trend surface fitting method can carry out complementary inspection using F distribution:
In formula, L is the number of multi-trend equation coefficient;N is monitoring points;The freedom degree of Q1 is n-L;Q2 from
By spending for L-1.
Under normal circumstances, significant property coefficient is divided into 0.05 and 0.01 two rank, for given insolation level α, can there is F
(L-1, n-L) distribution table finds critical value Fα;As F > FαWhen, then it is assumed that the situation of change of trend surface reacting dose is significant, it is on the contrary then
It is not significant.
4. residual error is extracted
After microseismic event carries out trend surface analysis, it is 0 that the abnormal Traveltime data separated, which has mathematical expectation, side
The smallest statistical property of difference, integrally tends to normal distribution.This is that the quantitative criterion that travel time anomaly threshold value divides can be determined as:
For abnormal point, time most short point as on first arrival time plan view.
This point is the monitoring point nearest apart from underground EXIT POINT floor projection.
(8) at the time of generation by the available rock mass micro rupture of inversion method, position and property, and according to micro rupture
Size, intensity, failure density, so definitely under latent projective water point spatial position.
Specific embodiment described herein is only an example for the spirit of the invention.The neck of technology belonging to the present invention
The technical staff in domain can make various modifications or additions to the described embodiments or replace by a similar method
In generation, however, it does not deviate from the spirit of the invention or beyond the scope of the appended claims.
Claims (5)
1. a kind of seat earth lies concealed gushing water karst collapse col umn dynamic positioning method, which is characterized in that including:
Monitor preparation process:In the microtremor observation system of surface deployment monitoring range covering projective water point;
Signal acquisition step:Using the microseismic signals of the secondary focus of microtremor observation system acquisition, the secondary focus is by height
Bearing bottom plate water slug roof is formed;
Signal analysis step:Analyze the microseismic signals so that definitely under latent projective water point position;
The signal analysis step specifically includes:Microseismic event is extracted in the record of continuous uninterrupted sampling;From each microseism
Coordinate (x in monitoring point is extracted in eventi,yi,zi,ti), wherein i=1,2,3 ... n;Static correction is carried out to monitoring point coordinate, is made
All receiving points are all on a datum level;It is imaged using focus shortest path, finely quickly on time isopleth plan view
Determine gushing water range;
It is described to be specifically included using focus shortest path image-forming step:Step 5.1, it selects to may be assumed that the ground lower die for uniform dielectric
Type;Step 5.2, the first arrival time plan view for the same microseismic event that different micro seismic monitoring points receive is drawn;Step 5.3 uses
Trend surface analysis technology carries out residual error extraction to first arrival time plan view, finds out the position of the shortest point of time in plane;Residual error
Extraction specifically includes:
Step 6.1, the trend surface fitting method of first arrival time isopleth;
If focus is to the first arrival time t of monitoring pointiIt is by its Trend valueExceptional value uiWith interference value siComposition, i.e.,
I=1 in formula, 2,3 ... n;
Wherein, tiMonitor value when to walk,Trend value when to walk, uiExceptional value when to walk;
If there is n group monitor value (x on groundi,yi,zi,ti) wherein i=1,2,3 ... n, wherein (xi,yi,zi) it is monitoring point coordinate,
tiMonitor value when to walk, due to data carry out static correction, so use P order polynomial it is carried out trend surface fitting method for:
Wherein, P is the number of trend surface;b0,b1,b2,...,bL-1For undetermined coefficient;L is the number of fitting coefficient, L=(P+1)
(P+2)/2;
By n to monitoring point (xi,yi,zi) and it is corresponding walk when monitor value tiIt brings above formula into, according to least square method valuation principle, enables
According to extreme value theorem, to make Q reach minimum, large linear systems is resolved to above formula derivation, L coefficient can be found out
Value, finds out fitting surface;
Step 6.2, degree of fitting calculates;
If ti,Monitor value and Trend value when to be walked on i-th of monitoring point,For its average value, total sum of squares of deviations
Wherein, Q2 numerical value is bigger, then fitting degree is higher;In formula,
Step 6.3, significance test
The effect of trend surface fitting method carries out complementary inspection using F distribution:
In formula, L ' is the number of multi-trend equation coefficient;N is monitoring points;The freedom degree of Q1 is n-L;The freedom of Q2
Degree is L-1;
For given insolation level α, table look-up to obtain critical value Fα;As F > FαWhen, then it is assumed that the variation feelings of trend surface reacting dose
Condition is significant, on the contrary then not significant;
Step 6.4, residual error is extracted;
After microseismic event carries out trend surface analysis, it is 0 that the abnormal Traveltime data separated, which has mathematical expectation, and variance is most
Small statistical property, integrally tends to normal distribution;This is that the quantitative criterion that travel time anomaly threshold value divides can be determined as:
For abnormal point, time most short point as on first arrival time plan view;
This point is the monitoring point nearest apart from underground EXIT POINT floor projection.
2. a kind of seat earth according to claim 1 lies concealed gushing water karst collapse col umn dynamic positioning method, it is characterised in that:Institute
It states in monitoring preparation process, first lies concealed EXIT POINT for coarse positioning in the sparse laying micro seismic monitoring point of larger range;Then, foundation
The microseismic event that the working face drawn a circle to approve when coarse positioning nearby occurs, encryption micro seismic monitoring point lie concealed EXIT POINT for fine positioning;Instead
Coarse positioning and fine positioning are carried out again until positioning accuracy is met the requirements.
3. a kind of seat earth according to claim 1 lies concealed gushing water karst collapse col umn dynamic positioning method, which is characterized in that logical
It crosses inversion method and obtains rock mass micro rupture at the time of occur, position and property, and according to the size of micro rupture, intensity, broken
Split density, so definitely under latent projective water point spatial position.
4. a kind of seat earth according to claim 1 lies concealed gushing water karst collapse col umn dynamic positioning method, which is characterized in that institute
Static correction step is stated to specifically include:
Step 4.1, datum elevation point Z is chosen0;Monitor the velocity of wave V that Earthquakes wave is propagated;
Step 4.2, the static correction value Δ t of each monitoring point, formula are calculated:
In formula, wherein Δ t is positive correcting value;ziTo monitor point height;z0For benchmark face elevation;V is velocity of wave;
Step 4.3, according to the static correction value of calculating, each monitoring point is corrected, makes it all on a datum level.
5. a kind of seat earth according to claim 1 lies concealed gushing water karst collapse col umn dynamic positioning method, which is characterized in that shake
Source positioning inversion method has:Three circle intersection positioning modes, longitudinal and shear wave time difference method, homotype wave time difference method, polarization analysis positioning mode,
Geiger seismic source positioning method, Geiger revised law, joint inversion method;Can be selected according to different situations different inversion methods come
Determine the spatial position of latent projective water point.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710657566.5A CN107450098B (en) | 2017-08-03 | 2017-08-03 | A kind of latent gushing water karst collapse col umn dynamic positioning method of seat earth |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710657566.5A CN107450098B (en) | 2017-08-03 | 2017-08-03 | A kind of latent gushing water karst collapse col umn dynamic positioning method of seat earth |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107450098A CN107450098A (en) | 2017-12-08 |
CN107450098B true CN107450098B (en) | 2018-11-23 |
Family
ID=60490460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710657566.5A Active CN107450098B (en) | 2017-08-03 | 2017-08-03 | A kind of latent gushing water karst collapse col umn dynamic positioning method of seat earth |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107450098B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109839663B (en) * | 2019-03-20 | 2020-04-10 | 山西山地物探技术有限公司 | Earthquake recognition method and device for hidden collapse column |
CN110633557B (en) * | 2019-10-30 | 2023-04-14 | 太原理工大学 | Identification method for favorable area of coal bed gas structure |
CN110907991B (en) * | 2019-12-11 | 2021-03-16 | 重庆大学 | Seismic source positioning method and system based on data field potential value and readable storage medium |
CN111577279B (en) * | 2020-05-18 | 2021-04-23 | 中国矿业大学 | Coal-geothermal water collaborative mining method based on collapse column water guide channel |
CN114017016B (en) * | 2021-11-04 | 2024-02-13 | 中煤科工集团西安研究院有限公司 | Classifying method for hidden water guide channels of coal seam floor |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105607121A (en) * | 2016-02-02 | 2016-05-25 | 中国矿业大学(北京) | Coal collapse column identification method and apparatus |
CN106437843A (en) * | 2016-08-30 | 2017-02-22 | 大连理工大学 | Coal mine bottom plate water guiding channel identification method based on micro-seismic monitoring |
CN106990435A (en) * | 2017-06-07 | 2017-07-28 | 中煤科工集团西安研究院有限公司 | It is a kind of to weaken the microseism localization method and device for relying on first break pickup precision |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9354336B2 (en) * | 2011-10-19 | 2016-05-31 | Global Ambient Seismic, Inc. | Microseismic data acquisition array and corresponding method |
-
2017
- 2017-08-03 CN CN201710657566.5A patent/CN107450098B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105607121A (en) * | 2016-02-02 | 2016-05-25 | 中国矿业大学(北京) | Coal collapse column identification method and apparatus |
CN106437843A (en) * | 2016-08-30 | 2017-02-22 | 大连理工大学 | Coal mine bottom plate water guiding channel identification method based on micro-seismic monitoring |
CN106990435A (en) * | 2017-06-07 | 2017-07-28 | 中煤科工集团西安研究院有限公司 | It is a kind of to weaken the microseism localization method and device for relying on first break pickup precision |
Non-Patent Citations (6)
Title |
---|
The Application and Prospect of Microseismic Technique in Coalmine;Si Qin et al.;《Procedia Environmental Sciences》;20121231;第218-224页 * |
微震监测技术在煤矿突水监测中的应用;白越等;《辽宁工程技术大学学报(自然科学版)》;20100831;第29卷(第4期);第549-552页 * |
煤层底板突水监测预警系统的开发及应用;靳德武等;《煤炭科学技术》;20111130;第39卷(第11期);第14-17页 * |
煤矿井下地震勘探技术应用现状与发展展望;程建远等;《勘探地球物理进展》;20090430;第32卷(第2期);第96-100、111页 * |
煤矿底板突水防治微震监测系统的建立及应用;王苏健等;《陕西省煤炭学会学术年会论文集(2016)》;20161231;第33-37页、附图1、4-8 * |
高精度微震监测技术在煤矿突水监测中的应用;姜福兴等;《岩石力学与工程学报》;20080930;第27卷(第9期);第1932-1938页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107450098A (en) | 2017-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107450098B (en) | A kind of latent gushing water karst collapse col umn dynamic positioning method of seat earth | |
Zhou et al. | Seepage channel development in the crown pillar: Insights from induced microseismicity | |
CN108930554B (en) | Well-ground-air combined monitoring method of the fracture of coal mine overlying strata and surface strata movement | |
Boncio et al. | Zoning surface rupture hazard along normal faults: Insight from the 2009 M w 6.3 L’Aquila, Central Italy, earthquake and other global earthquakes | |
CN102645674B (en) | Three-dimensional detection method for use before exploitation of working surface coal bed | |
CN103399354B (en) | The forecasting procedure of underground river tunnel geology | |
CN103389523B (en) | The forecast of tunnel geology and the method and system of construction | |
Hao et al. | Present-day crustal vertical motion around the Ordos block constrained by precise leveling and GPS data | |
CN102426384A (en) | Method for detecting underground goaf and karst distribution | |
CN103038670A (en) | Processing geophysical data | |
CN108415066A (en) | Tunnel construction geological disaster forecasting method | |
CN108104822B (en) | A kind of Tunnel Passing large size karst cavity construction geology disaster real-time predicting method | |
CN106597528A (en) | Tunnel geology three-dimensional seismic reflection and seismo-electric integrated advanced exploration device and method | |
Shynar et al. | Methodology of assessment and prediction of critical condition of natural-technical systems | |
CN103399356A (en) | Forecasting method and system for tunnel geology | |
CN108915766A (en) | A kind of working face deep concealed conduit pipe method for surveying | |
CN109059847A (en) | Overburden of the goaf destroys monitoring method | |
Ellenberger et al. | Three dimensional microseismic monitoring of a Utah longwall | |
Slaker et al. | Analysis of extensometer, photogrammetry and laser scanning monitoring techniques for measuring floor heave in an underground limestone mine | |
Nie et al. | Comprehensive ahead prospecting of tunnels in severely weathered rock mass environments with high water inrush risk: a case study in Shaanxi Province | |
Panzhin et al. | Satellite geodesy-aided geodynamic monitoring in mineral mining in the Urals | |
CN109521467A (en) | A kind of forward probe method based on projecting coal bed tunnel | |
CN112799129B (en) | Method and system for identifying water guide channel | |
Heasley et al. | An analysis of rock failure around a deep longwall using microseismics | |
Paine et al. | Assessing sinkhole potential at Wink and Daisetta using gravity and radar interferometry |
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 |