CN106019371B - A kind of advanced qualitative forecast method of projecting coal bed tunnel craven fault - Google Patents

A kind of advanced qualitative forecast method of projecting coal bed tunnel craven fault Download PDF

Info

Publication number
CN106019371B
CN106019371B CN201610316770.6A CN201610316770A CN106019371B CN 106019371 B CN106019371 B CN 106019371B CN 201610316770 A CN201610316770 A CN 201610316770A CN 106019371 B CN106019371 B CN 106019371B
Authority
CN
China
Prior art keywords
point
receiving point
ray
driving face
max
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
Application number
CN201610316770.6A
Other languages
Chinese (zh)
Other versions
CN106019371A (en
Inventor
王勃
刘盛东
周福宝
黄兰英
章俊
郝家林
姜永虎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN201610316770.6A priority Critical patent/CN106019371B/en
Publication of CN106019371A publication Critical patent/CN106019371A/en
Application granted granted Critical
Publication of CN106019371B publication Critical patent/CN106019371B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/65Source localisation, e.g. faults, hypocenters or reservoirs

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a kind of projecting coal bed advanced qualitative forecast methods of tunnel craven fault, this method in front of driving face by being arranged shot point, above driving face top plate, one receiving point is respectively set below bottom plate, above the top plate of driving face rear, one receiving point is respectively set below bottom plate, vector seismic detector is recycled to obtain the X of each receiving point respectively, Y-component seismic signal, according to the strong energy feature of coal seam breakpoint diffracted wave received, determine X, Y-component maximum value amplitude time, main polarization direction is calculated according to the signal of the time, further determine that the top broken-point of tomography and lower breakpoint, position and the property of craven fault are determined with this.The present invention can complete the calculating of main polarization direction merely with the signal of peak swing time, convenient, fast, it is easy to accomplish.The present invention can accurate judgement breakpoint location and property, the advanced qualitative forecast of craven fault is completed, to ensureing that projecting coal bed tunnel safety driving has important practical significance.

Description

A kind of advanced qualitative forecast method of projecting coal bed tunnel craven fault
Technical field
The present invention relates to a kind of heading tomography forecasting procedure, specifically a kind of projecting coal bed tunnel craven fault is fixed in advance Property forecasting procedure, belongs to technical field of mine safety.
Background technology
Coal and gas prominent refers under pressure, being dished out from coal body to tunnel or stope space within the extremely short time The event of a large amount of coal and (or) gas, not only results in casualties, can also national wealth be caused to lose.
Geological structure is the principal element for inducing coal and gas prominent, with tomography, particularly craven fault (drop is less than 5m) It is the most typical.Coal in Huainan Mining Area shows that crossdrift and coal roadway tunneling work in mining area with gas outburst point tectonic geology statistics Coal and gas prominent in face has 71.8% to be happened near craven fault.The qualitative advanced prediction of craven fault is exactly to position of fault The prediction carried out with fault feature judges, is of great significance for preventing coal and gas prominent.
Underground craven fault detects means mainly probing, physical prospecting at present.It is less efficient using drilling method, and visited before geology Drill easy off-design track, especially dip angle of hole be 0~25 ° of angle of depression and vertical angle when, vertical direction deviate it is more serious, Craven fault is easily caused to fail to judge.Using the earthquake advanced detection technology in geophysical prospecting method judge breakpoint at least need one it is complete when Between the period echo-signal, however since craven fault size is smaller, coal seam breakpoint diffracted wave is difficult to form a symmetric form complete Periodic signal, be unable to accurate judgement craven fault coal seam breakpoint, it is more difficult to which just/inverse property of craven fault is identified.And it is small Just/inverse property of tomography is very notable to coal and gas prominent control action difference, statistics show coal and gas prominent usually with Reversed fault is related.
Invention content
The purpose of the present invention is to provide a kind of projecting coal bed advanced qualitative forecast methods of tunnel craven fault, pass through this method It being capable of accurate judgement position of fault and fault feature.
To achieve the above object, a kind of projecting coal bed advanced qualitative forecast method of tunnel craven fault of the present invention, including it is following Step:
A shot point is arranged in step 1 in front of driving face, and first is arranged above driving face top plate and connects Secondary destination is arranged below driving face bottom plate in sink, and third is arranged above driving face rear, top plate and receives The 4th receiving point is arranged in point below driving face rear, bottom plate;
Step 2 selects the first receiving point, X, Y-component seismic signal is obtained using vector seismic detector, according to roof rock The strong energy feature of coal seam breakpoint diffracted wave received determines X, Y-component maximum value amplitude time t1
Step 3 calculates t1Main polarization direction angle θ1, specifically calculating step is:
A. to X, Y, this two multicomponent seismics signal does Hilbert transformation
X (t), y (t) are respectively X, Y-component, symbol in formulaIndicate Hilbert transformation,
B. Hermitian matrix constructions are established
C (t)=M* (t) M (t)
In formula, M (t)=[hx (t) hy (t)], the complex conjugate transposition of symbol * representing matrixes;
C. the maximum eigenvalue λ of Hermitian matrixes is soughtmaxAnd its corresponding normalized feature vector (xmax, ymax);
d.Re (x in formulamax)、Re(ymax) it is respectively xmax、ymaxReal part;
Step 4, according to the first receiving point position and main polarization direction angle θ1, the first receiving point is drawn to the of craven fault One ray;
Step 5 selects other receiving points, and method draws secondary destination to craven fault respectively with step 2 to step 4 The second ray, third receiving point to craven fault third ray and the 4th receiving point to craven fault the 4th ray;
6th step, the intersection point of the first ray and third ray are the top broken-point of craven fault, the second ray and the 4th ray Intersection point is the lower breakpoint of craven fault;The upper and lower disk relationship in coal seam is determined according to top broken-point and lower breakpoint, judges fault properties.
The present invention eliminates the reliance on a complete cycle signal, merely with the letter of the peak swing time of coal seam breakpoint diffracted wave The calculating of main polarization direction number can be completed, it is convenient, fast, it is easy to accomplish.The present invention can accurate judgement craven fault position and Feature completes the advanced qualitative forecast to craven fault, to ensureing that projecting coal bed tunnel safety driving has important practical significance.
Description of the drawings
Fig. 1 is the method for the present invention schematic diagram;
Fig. 2 is that the present invention obtains vector seismic detector acquisition X, Y-component peak swing time method schematic diagram;
In figure, the 1, first receiving point, 2, secondary destination, 3, third receiving point, the 4, the 4th receiving point, 5, headwork Face, 6, shot point, 7, top plate, 8, bottom plate, 9, top broken-point, 10, lower breakpoint, the 11, first ray, the 12, second ray, 13, third Ray, the 14, the 4th ray, 15, craven fault.
Specific implementation mode
Below in conjunction with the accompanying drawings, the present invention is further described.
A kind of advanced qualitative forecast method of projecting coal bed tunnel craven fault, which is characterized in that include the following steps:
Step 1 a, as shown in Figure 1, shot point 6 is arranged in 5 front of driving face, in 5 top plate 7 of driving face Top be arranged the first receiving point 1, below 5 bottom plate 8 of driving face be arranged secondary destination 2,5 rear of driving face, 7 top setting third receiving point 3 of top plate, is arranged the 4th receiving point 4 below 5 rear of driving face, bottom plate 8;
Step 2 selects the first receiving point 1, as shown in Fig. 2, obtaining X, Y-component seismic signal, root using vector seismic detector The strong energy feature of coal seam breakpoint diffracted wave received according to roof rock determines X, Y-component maximum value amplitude time t1
Step 3 calculates t1Main polarization direction angle θ1, specifically calculating step is:
A. to X, Y, this two multicomponent seismics signal does Hilbert transformation
X (t), y (t) are respectively X, Y-component, symbol in formulaIndicate Hilbert transformation,
B. Hermitian matrix constructions are established
C (t)=M* (t) M (t)
In formula, M (t)=[hx (t) hy (t)], the complex conjugate transposition of symbol * representing matrixes;
C. the maximum eigenvalue λ of Hermitian matrixes is soughtmaxAnd its corresponding normalized feature vector (xmax, ymax);
d.Re (x in formulamax)、Re(ymax) it is respectively xmax、ymaxReal part;
Step 4, according to 1 position of the first receiving point and main polarization direction angle θ1, draw the first receiving point 1 and arrive craven fault 15 The first ray 11;
Step 5 selects other receiving points, and method draws secondary destination 2 and break to small respectively with step 2 to step 4 Second ray 12 of layer 15, the third ray 13 of third receiving point 3 to craven fault 15 and the 4th receiving point 4 are to the of craven fault 15 Four rays 14;
6th step, as shown in Figure 1, the intersection point of the first ray 11 and third ray 13 is coal seam top broken-point 9, the second ray 12 Intersection point with the 4th ray 14 is breakpoint 10 under coal seam;The upper and lower disk relationship in coal seam is determined according to top broken-point 9 and lower breakpoint 10, is sentenced Disconnected fault properties.Hanging wall shown in FIG. 1 rises, and is reversed fault.
The present invention eliminates the reliance on a complete cycle signal, merely with the letter of the peak swing time of coal seam breakpoint diffracted wave The calculating of main polarization direction number can be completed, it is convenient, fast, it is easy to accomplish.The present invention can accurate judgement craven fault position and Feature completes the advanced qualitative forecast to craven fault, to ensureing that projecting coal bed tunnel safety driving has important practical significance.
Preferably, at 5 three meters of front of heading driving face, first connects the setting of shot point 6 described in step 1 Sink 1 is arranged above the top plate 7 of driving face 5 at three meters, and secondary destination 2 is arranged under the bottom plate 8 of driving face 5 Place of three meters of side, the setting of third receiving point 3 above five meters of 5 rear of driving face (lefts Fig. 1 be rear), the top plate 7 at six meters, The setting of 4th receiving point 4 is below five meters of 5 rear of driving face, the bottom plate 8 at six meters.The selection of these points is both convenient for setting It sets, and the requirement of detection accuracy can be met.

Claims (1)

1. a kind of advanced qualitative forecast method of projecting coal bed tunnel craven fault, which is characterized in that include the following steps:
Step 1 arranges a shot point (6), on face driving face (5), top plate (7) in front of driving face (5) Side's the first receiving point of setting (1), the setting secondary destination (2) below face driving face (5), bottom plate (8), in heading driver Make face (5) rear, top plate (7) top setting third receiving point (3), the setting below driving face (5) rear, bottom plate (8) 4th receiving point (4);
Step 2 selects the first receiving point (1), X, Y-component seismic signal is obtained using vector seismic detector, according to the coal received The layer strong energy feature of breakpoint diffracted wave, determines X, Y-component maximum value amplitude time t1
Step 3 calculates t1Main polarization direction angle θ1, specifically calculating step is:
A. to X, Y, this two multicomponent seismics signal does Hilbert transformation
X (t), y (t) are respectively X, Y-component, symbol in formulaIndicate Hilbert transformation,
B. Hermitian matrix constructions are established
In formula, M (t)=[hx (t) hy (t)], the complex conjugate transposition of symbol * representing matrixes;
C. the maximum eigenvalue λ of Hermitian matrixes is soughtmaxAnd its corresponding normalized feature vector (xmax,ymax);
d.Re (x in formulamax)、Re(ymax) it is respectively xmax、ymaxReal part;
Step 4, according to the first receiving point (1) position and main polarization direction angle θ1, draw the first receiving point (1) and arrive craven fault (15) The first ray (11);
Step 5 selects other receiving points, and method draws secondary destination (2) and arrive craven fault respectively with step 2 to step 4 (15) the third ray (13) of the second ray (12), third receiving point (3) to craven fault (15) and the 4th receiving point (4) arrive small 4th ray (14) of tomography (15);
The intersection point of 6th step, the first ray (11) and third ray (13) is coal seam top broken-point (9), the second ray (12) and the 4th The intersection point of ray (14) is breakpoint (10) under coal seam;The upper and lower disk relationship in coal seam is determined according to top broken-point (9) and lower breakpoint (10), Judge fault properties;
Shot point described in step 1 (6) setting is in front of the heading driving face (5) at three meters, the first receiving point (1) It is arranged at three meters of top plate (7) top of face driving face (5), secondary destination (2) is arranged in face driving face (5) below bottom plate (8) at three meters, third receiving point (3) setting is six above six meters of driving face (5) rear, top plate (7) At rice, the setting of the 4th receiving point (4) is at six meters of driving face (5) rear, six meters of bottom plate (8) lower section.
CN201610316770.6A 2016-05-13 2016-05-13 A kind of advanced qualitative forecast method of projecting coal bed tunnel craven fault Active CN106019371B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610316770.6A CN106019371B (en) 2016-05-13 2016-05-13 A kind of advanced qualitative forecast method of projecting coal bed tunnel craven fault

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610316770.6A CN106019371B (en) 2016-05-13 2016-05-13 A kind of advanced qualitative forecast method of projecting coal bed tunnel craven fault

Publications (2)

Publication Number Publication Date
CN106019371A CN106019371A (en) 2016-10-12
CN106019371B true CN106019371B (en) 2018-09-14

Family

ID=57100740

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610316770.6A Active CN106019371B (en) 2016-05-13 2016-05-13 A kind of advanced qualitative forecast method of projecting coal bed tunnel craven fault

Country Status (1)

Country Link
CN (1) CN106019371B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106646662B (en) * 2016-11-10 2018-09-18 中国矿业大学(北京) The prediction technique and device in Gas Outburst region
CN109765606B (en) * 2019-01-28 2020-08-04 阳泉煤业(集团)有限责任公司 Method for detecting nature of hidden fault of stope face based on reflected trough wave
CN110531418B (en) * 2019-08-21 2020-11-20 徐州工程学院 Breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging
CN110531415B (en) * 2019-08-21 2020-10-30 徐州工程学院 Three-dimensional small fault advanced detection method utilizing influence of surrounding rock loosening ring
CN111025383B (en) * 2019-11-21 2021-09-24 徐州工程学院 Method for qualitatively judging water filling condition of tunnel front karst cave based on diffracted transverse waves
CN111236940B (en) * 2020-01-14 2021-06-01 山西晋城无烟煤矿业集团有限责任公司 Method for safely and efficiently passing reverse fault group on fully mechanized coal mining face
CN112859196B (en) * 2021-03-03 2022-05-24 中国石油大学(北京) Accurate identification method for broken layer breakpoint in shaft

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102053267A (en) * 2010-10-22 2011-05-11 中国石油化工股份有限公司 Method for separating VSP (vertical seismic profiling) wave field based on parametric inversion during seismic profile data processing
CN104166136A (en) * 2014-07-11 2014-11-26 河海大学 Interference subspace tracking-based high-efficiency self-adaptive monopulse angle measurement method

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8325559B2 (en) * 2010-08-27 2012-12-04 Board Of Regents Of The University Of Texas System Extracting SV shear data from P-wave marine data
CN102540254B (en) * 2010-12-29 2014-11-26 中国石油天然气集团公司 Method for determining effective excitation region based on lighting analysis
CN103235334B (en) * 2013-04-11 2015-06-10 中国矿业大学 Low drainage roadway fault-throw advanced-detection method based on seismic signals
CN103217719B (en) * 2013-04-11 2015-07-15 中国矿业大学 Method of advanced detection of breaking-loss wing coal seam of coal road based on single offset pair observation system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102053267A (en) * 2010-10-22 2011-05-11 中国石油化工股份有限公司 Method for separating VSP (vertical seismic profiling) wave field based on parametric inversion during seismic profile data processing
CN104166136A (en) * 2014-07-11 2014-11-26 河海大学 Interference subspace tracking-based high-efficiency self-adaptive monopulse angle measurement method

Also Published As

Publication number Publication date
CN106019371A (en) 2016-10-12

Similar Documents

Publication Publication Date Title
CN106019371B (en) A kind of advanced qualitative forecast method of projecting coal bed tunnel craven fault
Ellsworth et al. Triggering of the Pohang, Korea, earthquake (M w 5.5) by enhanced geothermal system stimulation
CN104166164B (en) The many wave reflection 3-D seismics forward probe methods of coal roadway tunneling geological structure three-component
CN105785471A (en) Impact danger evaluation method of mine pre-exploiting coal seam
CN101270666B (en) Tunnel influence geological condition section-expanding horizontal drilling estimation and prediction method
CN105334548A (en) Geological forecasting method for tunnel construction in karst area
CN103235333A (en) Coal roadway fault advanced-detection method based on reflecting channel wave signals
CN103217719B (en) Method of advanced detection of breaking-loss wing coal seam of coal road based on single offset pair observation system
CN202649483U (en) Electric field constraint method mine security type full mechanized excavation machine carried geological structure detection system
CN102213773A (en) Roadway multi-azimuth advance detection method
CN109736796A (en) A kind of intensification blasthole forward probe mechanism and its forecasting procedure
CN112965136A (en) Multi-stage advanced detection method for water-rich karst tunnel
Tian et al. Coalbed methane reservoir fracture evaluation through the novel passive microseismic survey and its implications on permeable and gas production
CN103235334B (en) Low drainage roadway fault-throw advanced-detection method based on seismic signals
CN103176214B (en) Electric field leash law coal peace type roadheader carries tectonic structure detection system and method thereof
CN103606019A (en) Mine goaf overlying stratum sedimentation dynamic prediction method based on time-space relationship
Meng et al. In situ investigation and numerical simulation of the failure depth of an inclined coal seam floor: a case study
Lin et al. Transient electromagnetic detection method in water-sealed underground storage caverns
Tan et al. Researches on floor stratum fracturing induced by antiprocedure mining underneath close-distance goaf
CN113126147A (en) Detection method for determining spatial form of hidden collapse column of rock roadway floor
Ptáček et al. The impact of Outer Western Carpathian nappe tectonics on the recent stress-strain state in the Upper Silesian Coal Basin (Moravosilesian Zone, Bohemian Massif)
CN108228920A (en) A kind of sandstone-type uranium mineralization with respect soil element uranium exception screening technique
Cheng et al. Rapidly locating a water–inrush collapse column in a seam floor: A case study
Shang et al. Advanced predication of geological anomalous body ahead of laneway using seismic tomography technique
CN110748341A (en) Method for detecting influence radius of controllable seismic source based on carbon dioxide blasting

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant