CN108121010B - Based on the united underground dead face slot wave forward probe method and system in hole lane - Google Patents
Based on the united underground dead face slot wave forward probe method and system in hole lane Download PDFInfo
- Publication number
- CN108121010B CN108121010B CN201711391055.XA CN201711391055A CN108121010B CN 108121010 B CN108121010 B CN 108121010B CN 201711391055 A CN201711391055 A CN 201711391055A CN 108121010 B CN108121010 B CN 108121010B
- Authority
- CN
- China
- Prior art keywords
- wave
- slot wave
- slot
- tunnel
- forward probe
- 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
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Abstract
The present invention relates to a kind of slot wave forward probe method and system, belong to technical field of geophysical exploration, and in particular to one kind is based on the united underground dead face slot wave forward probe method and system in hole lane.The present invention, which utilizes, beats the hole lane combined detection scheme of multi-component and multilevel geophone in deep hole and placement hole in two lateral wall vertical working direction of tunnel, the detection angle of reflection slot wave is effectively expanded, so as to front anomalous body, the especially reflection slot wave image quality of the reflecting surface anomalous body vertical with tunnel improves dead face reflection slot wave forward probe precision.
Description
Technical field
The present invention relates to a kind of slot wave forward probe method and system, belong to technical field of geophysical exploration, specifically relate to
And it is a kind of based on the united underground dead face slot wave forward probe method and system in hole lane.
Background technique
For digging laneway front geological construction advanced detection technology be a detection bad environments, technical difficulty greatly,
Detection accuracy requires geological guarantee technologies that are high and requiring real-time.But this technology can not only prevent Gas Outburst, dash forward
The accidents such as water occur, moreover it is possible to geological information are provided for gas drainage under suction, advanced leting speeper, optimization tunneling program etc., thus improving
Accelerate digging speed in the case where safety guarantee.
In the prior art, mining progress is not influenced using big gun pick as the real-time advanced detection technology of focus, while driving
It can complete the task of real-time forward probe.
But due to driving face narrow space, Hydrographic General Line in the prior art is limited to linearly, so that receive
Wave field space wave number narrow limits, the bearing accuracy of reflection wave imaging are zero, so that slot wave image quality is poor, are influenced anti-
Penetrate slot wave forward probe precision.
Summary of the invention
Slot wave image quality present in the present invention mainly solution prior art is poor, reflects slot wave forward probe precision
Low technical problem provides a kind of based on the united underground dead face slot wave forward probe method and system in hole lane, the party
Method and system widen wavefield reception visual angle by via hole lane joint observation scheme, and use phase according to the wave field characteristics of acquisition data
The noise-removed technology and imaging technique answered, improve the precision of dead face channel wave seismic advanced prediction.
Above-mentioned technical problem of the invention is mainly to be addressed by following technical proposals:
One kind being based on the united underground dead face slot wave forward probe method in hole lane, comprising:
Epicenter excitation step utilizes epicenter excitation signal at the tunnel place of meeting head on；
Wavefield reception step, be utilized respectively be installed on tunnel meet head on place two lateral wall lateral apertures in multiple holes in wave detector and
Be installed on tunnel meet head on two lateral wall of rear tunnel wave detector receive vibration signal caused by wave field；
Wave field processing step deviates wave field record using scattering imaging method for handling the wave field received
Imaging obtains slot wave forward probe imaging results；
Geologic interpretation step, explained on migration imaging section in front of dead face and around existing for suspicious geology it is different
Chang Ti, and determine its position, scale and type.
Preferably, above-mentioned one kind is based on the united underground dead face slot wave forward probe method in hole lane, the wave field
Processing step includes Signal Pretreatment sub-step, for being decayed using polarization filtering to the random noise in data record,
The signal-to-noise ratio that useful signal especially reflects slot wave signal is improved, is specifically included:
(1) to the covariance matrix for calculating three component seismic data in timing window:
(2) according to formula | Q- λ I |=0, find out eigenvalue λ1、λ2And λ3With corresponding feature vector V1、V2And V3；
(3) polarization coefficient is calculatedAnd utilize polarization coefficient structure
Build out Polarization filter
(4) Polarization filter is respectively acting on three components, just obtains filtered record, filtered record filters out
The wave field of low polarizability, remains the wave field of high polarization degree.
(5) window when mobile, and repeat the above steps until record tail.
Preferably, above-mentioned one kind is based on the united underground dead face slot wave forward probe method in hole lane, the wave field
Processing step specifically includes slot wave compression sub-step, specially utilizes the multiple filtering technique high precision computation slot based on S-transformation
Wave number is composed according to frequency dispersion, picks up frequency dispersion parameter in frequency dispersion spectrum, and with frequency dispersion parameter designing inverse dispersion filter, then act on slot
Wave record, realizes the compression to the slot wave wave train.
Preferably, above-mentioned one kind is based on the united underground dead face slot wave forward probe method in hole lane, the wave field
Processing step specifically includes modulation function and determines sub-step, for according to the spatial position of excitation point, receiving point to reflection point and
Spread speed parameter is based on ray tracing technique, calculates Travel time t, is determined according to Travel time with the signal period
Mathematical statistics means are based on by window amplitude information when received three-component for the when window of unit, calculate particle polarization direction
With back wave directions of raysSpace angle θ, according to the polarization direction coefficient of angle calculation imaging slot wave type usedDetermine the modulation function of polarization offset imagingWherein α is given tuning parameter, when wave field is LOVE
When slot wave wave fieldThen it is if RayLeigh slot wave
Preferably, above-mentioned one kind is based on the united underground dead face slot wave forward probe method in hole lane, the wave field
Handling in sub-step includes migration imaging sub-step, specially acts on record using modulation function, obtains polarization filtering signal
X'(t)=X (t) * F (t), then using scattering imaging method to record X'(t) migration imaging is carried out, it is advanced to obtain final slot wave
Detection imaging result.
One kind being based on the united underground dead face slot wave forward probe system in hole lane, comprising:
Big gun digs focus, meets head on to locate positioned at tunnel, for exciting vibration signal；
Wave detector in hole is installed on the tunnel place of meeting head on two sides and helps in lateral aperture, and each lateral aperture installs more than two holes
Middle wave detector, for receiving wave field caused by vibration signal；
Tunnel wave detector, along tunnel meet head on two lateral wall of rear setting, composition two be parallel to each other along, for receive shake
Wave field caused by dynamic signal；
Signal processing apparatus deviates wave field record using scattering imaging method for handling the wave field received
Imaging, explained on migration imaging section in front of dead face and around existing for suspicious geological anomalous body, and determine its position
It sets, scale and type.
Preferably, above-mentioned one kind is based on the united underground dead face slot wave forward probe system in hole lane, in the hole
6 grades or more three-component geophones are laid in wave detector 10 meters of spacing settings in lateral aperture, every hole.
Preferably, above-mentioned one kind is based on the united underground dead face slot wave forward probe system in hole lane, the tunnel
24 or more three-component geophones are arranged in 10 meters of geophone separation settings, every survey line.
Therefore, compared with prior art, the present invention, which utilizes, beats deep hole and placement hole in two lateral wall vertical working direction of tunnel
The hole lane combined detection scheme of middle multi-component and multilevel geophone has effectively expanded the detection angle of reflection slot wave, so as to front exception
It is super to improve dead face reflection slot wave for the reflection slot wave image quality of body, especially the reflecting surface anomalous body vertical with tunnel
Preceding detection accuracy.
Detailed description of the invention
Attached drawing 1 is observation system structural schematic diagram of the invention；
Attached drawing 2 is tunnel geophone signal processing schematic of the invention.
Attached drawing 3 is geophone signal processing schematic in hole of the invention.
Attached drawing 4 is that wave detector compares the Effect on Detecting of the tomography perpendicular to tunnel in tunnel wave detector and hole of the invention
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:
The observation system of based on well the united underground coal mine dead face slot wave forward probe of the present embodiment is according to Fig. 1
Mode designs.
The signal that is excited using dynamite source during big gun pick of focal point first is as excitaton source, therefore the usual position of excitation point
It meets head on to locate in tunnel；Tunnel receiving point is distributed in tunnel and meets head on 5 meters of rears, and two lateral walls respectively lay two surveys line, and every distribution is more
A three-component receiving point (it is recommended that at least 24, preferably to identify wave field), is divided into 10 meters；It meets head on two at 5 meters in tunnel
Side is each to make a call to the lateral aperture that a depth is greater than 60 meters, and 10 meters of detector interval in hole, 6 grades of three-component geophones are laid in every hole；
Then three-component geophone is connect using data line with from note formula in-seam seismograph, is excited when receiving big gun pick
Direct wave caused by vibration signal out, refracted wave, back wave and other wave fields when encountering anomalous geologic body around；
Decayed using polarization filtering to the random noise in data record, improves useful signal and especially reflect slot wave
The signal-to-noise ratio of signal, specific algorithm are as follows:
1. to the covariance matrix for calculating three component seismic data in timing window:
Wherein, the covariance matrix formed in window when Q is: (Xi,Yi,Zi) be when window in tri- components of X, Y and Z sampling point,For when window in sampling point average amplitude value.
2. according to formula | Q- λ I |=0, find out eigenvalue λ1、λ2And λ3With corresponding feature vector V1、V2And V3；
3. calculating polarization coefficientAnd utilize polarization coefficient structure
Build out Polarization filter
4. Polarization filter is respectively acting on three components, filtered record is just obtained, filtered record filters out
The wave field of low polarizability, remains the wave field of high polarization degree.
5. window when moving, and repeat the above steps until record tail.
Then it utilizes the multiple filtering technique high precision computation slot wave data frequency dispersion based on S-transformation to compose, is picked up in frequency dispersion spectrum
It takes frequency dispersion parameter, and with frequency dispersion parameter designing inverse dispersion filter, then acts on slot wave record, realize the pressure to the slot wave wave train
Contracting；
According to excitation point, the spatial position of receiving point to reflection point and spread speed parameter, it is based on ray tracing technique, meter
Travel time t is calculated, the when window as unit of the signal period is determined according to Travel time, passes through window when received three-component
Amplitude information is based on mathematical statistics means, calculates particle polarization direction(the corresponding feature vector of maximum eigenvalue) and reflection
Wave ray directionSpace angle θ (such as Fig. 2 and Fig. 3), according to the polarization side of angle calculation imaging slot wave type used
To coefficientDetermine the modulation function of polarization offset imaging, such as following formula:
Wherein α is given tuning parameter.When wave field is LOVE slot wave wave field
Then it is if RayLeigh slot wave
Record is acted on using modulation function, obtains polarization filtering signal X'(t)=X (t) * F (t), then using scattering
Imaging method is to record X'(t) migration imaging is carried out, obtain final slot wave forward probe imaging results；
By geologic interpretation, explained on migration imaging section in front of dead face and around existing for suspicious geology it is different
Chang Ti, and determine its position, scale and type.
Compared with prior art, the present invention, which utilizes, beats multistage in deep hole and placement hole in two lateral wall vertical working direction of tunnel
The hole lane combined detection scheme of wave detector has effectively expanded the detection angle of reflection slot wave, so as to front anomalous body, especially
It is the reflection slot wave image quality and areas imaging (as shown in Figure 4) of the reflecting surface anomalous body vertical with tunnel, improves only head
Slot wave forward probe precision is reflected in tunnel.
Claims (7)
1. one kind is based on the united underground dead face slot wave forward probe method in hole lane characterized by comprising
Epicenter excitation step utilizes epicenter excitation signal at the tunnel place of meeting head on；
Wavefield reception step is utilized respectively and is installed on wave detector and installation in multiple holes that tunnel is met head in two lateral wall lateral apertures of place
In tunnel meet head on two lateral wall of rear tunnel wave detector receive vibration signal caused by wave field；
Wave field processing step records wave field using scattering imaging method and carries out migration imaging for handling the wave field received,
Obtain slot wave forward probe imaging results；
Geologic interpretation step, explained on migration imaging section in front of dead face and around existing for suspicious geologic anomaly
Body, and determine its position, scale and type；
Wherein,
The wave field processing step specifically includes modulation function and determines sub-step, for according to excitation point, receiving point to reflection point
Spatial position and spread speed parameter, be based on ray tracing technique, calculate Travel time t, it is true according to Travel time
The fixed when window as unit of the signal period is based on mathematical statistics means by window amplitude information when received three-component, calculates matter
Point polarization directionWith back wave directions of raysSpace angle θ, according to the pole of angle calculation imaging slot wave type used
Change direction coefficientDetermine the modulation function of polarization offset imagingWherein α is given tuning parameter, when
When wave field is LOVE slot wave wave fieldThen it is if RayLeigh slot wave
2. it is according to claim 1 a kind of based on the united underground dead face slot wave forward probe method in hole lane, it is special
Sign is that the wave field processing step includes Signal Pretreatment sub-step, for using polarization filtering in data record with
Machine noise is decayed, and is improved the signal-to-noise ratio that useful signal especially reflects slot wave signal, is specifically included:
(1) to the covariance matrix for calculating three component seismic data in timing window:
(2) according to formula | Q- λ I |=0, find out eigenvalue λ1、λ2And λ3With corresponding feature vector V1、V2And V3；
(3) polarization coefficient is calculatedAnd it is constructed using polarization coefficient
Polarization filter
(4) Polarization filter is respectively acting on three components, just obtains filtered record, filtered record has filtered out low
The wave field of polarizability remains the wave field of high polarization degree；
(5) window when mobile, and repeat the above steps until record tail.
3. it is according to claim 1 a kind of based on the united underground dead face slot wave forward probe method in hole lane, it is special
Sign is that the wave field processing step specifically includes slot wave compression sub-step, specially utilizes the multiple filtering skill based on S-transformation
Art high precision computation slot wave data frequency dispersion spectrum is picked up frequency dispersion parameter in frequency dispersion spectrum, and is filtered with frequency dispersion parameter designing inverse dispersion
Then device acts on slot wave record, realizes the compression to the slot wave wave train.
4. it is according to claim 1 a kind of based on the united underground dead face slot wave forward probe method in hole lane, it is special
Sign is, includes migration imaging sub-step in the wave field processing sub-step, specially acts on record using modulation function, obtain
To polarization filtering signal X'(t)=X (t) * F (t), then using scattering imaging method to record X'(t) migration imaging is carried out, it obtains
Obtain final slot wave forward probe imaging results.
5. one kind is based on the united underground dead face slot wave forward probe system in hole lane characterized by comprising
Big gun digs focus, meets head on to locate positioned at tunnel, for exciting vibration signal；
Wave detector in hole is installed on the tunnel place of meeting head on two sides and helps in lateral aperture, and each lateral aperture installs two or more Kong Zhongjian
Wave device, for receiving wave field caused by vibration signal；
Tunnel wave detector, along tunnel meet head on two lateral wall of rear setting, composition two be parallel to each other along, for receive vibration letter
Wave field caused by number；
Signal processing apparatus records wave field using scattering imaging method and carries out migration imaging for handling the wave field received,
Explained on migration imaging section in front of dead face and around existing for suspicious geological anomalous body, and determine its position, rule
Mould and type；And according to excitation point, the spatial position of receiving point to reflection point and spread speed parameter, it is based on ray tracing skill
Art calculates Travel time t, determines the when window as unit of the signal period according to Travel time, passes through received three points
Window amplitude information when amount is based on mathematical statistics means, calculates particle polarization directionWith back wave directions of raysSpace angle
θ, according to the polarization direction coefficient of angle calculation imaging slot wave type usedDetermine the modulation function of polarization offset imagingWherein α is given tuning parameter, when wave field is LOVE slot wave wave fieldIf
RayLeigh slot wave is then
6. it is according to claim 5 a kind of based on the united underground dead face slot wave forward probe system in hole lane, it is special
Sign is that wave detector is arranged for 10 meters of spacing in lateral aperture in the hole.
7. it is according to claim 5 a kind of based on the united underground dead face slot wave forward probe system in hole lane, it is special
Sign is that 24 or more three-component geophones are arranged in 10 meters of tunnel geophone separation settings, every survey line.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711391055.XA CN108121010B (en) | 2017-12-21 | 2017-12-21 | Based on the united underground dead face slot wave forward probe method and system in hole lane |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711391055.XA CN108121010B (en) | 2017-12-21 | 2017-12-21 | Based on the united underground dead face slot wave forward probe method and system in hole lane |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108121010A CN108121010A (en) | 2018-06-05 |
| CN108121010B true CN108121010B (en) | 2019-05-07 |
Family
ID=62230643
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711391055.XA Active CN108121010B (en) | 2017-12-21 | 2017-12-21 | Based on the united underground dead face slot wave forward probe method and system in hole lane |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108121010B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108693561B (en) * | 2018-06-14 | 2019-11-08 | 中煤科工集团西安研究院有限公司 | Coal mining seismic acquisition system and method based on wave detector subdivision array |
| CN111399030A (en) * | 2020-02-29 | 2020-07-10 | 山西晋煤集团技术研究院有限责任公司 | Mine roadway slot wave distributed advanced detection system and method |
| CN113075730A (en) * | 2021-03-23 | 2021-07-06 | 西安中地博睿探测科技有限公司 | Travel time information automatic extraction method in transmission channel wave velocity CT imaging |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104111479A (en) * | 2014-07-17 | 2014-10-22 | 河北煤炭科学研究院 | Underground coal mine in-excavating channel wave seism fore-detection method |
| CN104199110A (en) * | 2014-09-05 | 2014-12-10 | 河北煤炭科学研究院 | Method for channel wave earthquake stereoscopic advanced detection during underground coal mine support |
| CN106610503A (en) * | 2016-11-29 | 2017-05-03 | 北京中矿大地地球探测工程技术有限公司 | Omnidirectional slot wave seismic detection method in coal mine excavation process |
| CN107015273B (en) * | 2017-04-10 | 2018-10-19 | 中煤科工集团西安研究院有限公司 | A kind of coal seam Love type groove wave pressure contracting method and devices based on frequency domain |
-
2017
- 2017-12-21 CN CN201711391055.XA patent/CN108121010B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108121010A (en) | 2018-06-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR102369752B1 (en) | Detection of underground structures | |
| CN102027390B (en) | Jointly interpolating and deghosting seismic data | |
| US9874646B2 (en) | Seismic data processing | |
| Liu et al. | Three-dimensional seismic ahead-prospecting method and application in TBM tunneling | |
| US6092025A (en) | Hydrocarbon edge detection using seismic amplitude | |
| CN101551463B (en) | Noise suppression evaluation method for three-dimensional observation system | |
| US20120275264A1 (en) | Removing free-surface effects from seismic data acquired in a towed survey | |
| US9274239B2 (en) | Wavefield deghosting | |
| CN104360388B (en) | A kind of 3 D seismic observation system evaluation method | |
| CN108121010B (en) | Based on the united underground dead face slot wave forward probe method and system in hole lane | |
| CN101598811A (en) | A kind of method of calculating shot-static correction for two-dimensional vertical seismic section data | |
| Brodic et al. | Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging | |
| WO2010082131A2 (en) | Processing seismic data | |
| Dasios et al. | Seismic imaging of the shallow subsurface: shear‐wave case histories | |
| CN104570115A (en) | Surface wave attenuating method and device | |
| US20050162974A1 (en) | Resonance scattering seismic method | |
| Stucchi et al. | High‐resolution reflection seismic survey at the Patigno landslide, Northern Apennines, Italy | |
| CN110579799A (en) | Seismic acquisition observation method and system with equal travel time intervals | |
| Pugin et al. | Update on recent observations in multi-component seismic reflection profiling | |
| Mari et al. | Well seismic surveying | |
| Blias et al. | High frequency VSP methodology and its application to the detailed investigation of near-well space | |
| Xia et al. | Feasibility of detecting voids with Rayleigh-wave diffraction | |
| Daley et al. | Borehole seismic monitoring of injected CO2 at the Frio site | |
| Tomio et al. | High-frequency seismic measurements using a piezoelectric type accelerometer array | |
| Koedel et al. | Project report on seismic tomography data interpretation and conceptual model for integrating DAS into borehole seismic tomography surveying |
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 |