CN106772595A - A kind of method that detonator time delay is eliminated in shock wave Velocity Inversion for colliery - Google Patents
A kind of method that detonator time delay is eliminated in shock wave Velocity Inversion for colliery Download PDFInfo
- Publication number
- CN106772595A CN106772595A CN201611103024.5A CN201611103024A CN106772595A CN 106772595 A CN106772595 A CN 106772595A CN 201611103024 A CN201611103024 A CN 201611103024A CN 106772595 A CN106772595 A CN 106772595A
- Authority
- CN
- China
- Prior art keywords
- time delay
- shock wave
- shot point
- ripple
- eliminated
- 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.)
- Pending
Links
- 230000035939 shock Effects 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000000523 sample Substances 0.000 claims abstract description 32
- 239000003245 coal Substances 0.000 claims abstract description 13
- 230000008030 elimination Effects 0.000 claims description 7
- 238000003379 elimination reaction Methods 0.000 claims description 7
- 238000000265 homogenisation Methods 0.000 claims description 6
- 238000004880 explosion Methods 0.000 claims description 5
- 230000001934 delay Effects 0.000 claims 1
- 229940079593 drug Drugs 0.000 claims 1
- 239000003814 drug Substances 0.000 claims 1
- 239000011435 rock Substances 0.000 abstract description 7
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000005065 mining Methods 0.000 abstract description 2
- 230000000644 propagated effect Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 3
- 238000005422 blasting Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000012502 risk assessment 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/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- 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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
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)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses the method that detonator time delay in a kind of shock wave Velocity Inversion for colliery is eliminated, it is adaptable to the P ripple first arrival-times of coal mining and technical field of mine safety, first labeled analysis probes records shocking waveshape；Calculate each shot point to the distance popped one's head in；Detonator time delay is then obtained with distance relation using linear least squares fit P ripples again；Detonator time delay is finally eliminated, the P ripples corrected are then.This method can eliminate detonator time delay so as to calculate the real time that P ripples are propagated in coal and rock exactly, improve the accuracy and reliability of colliery shock wave Velocity Inversion, and its reliable results is applied widely, and is easy to computer programming.
Description
Technical field
The present invention relates to the method that detonator time delay in a kind of shock wave Velocity Inversion for colliery is eliminated, belong to coal mining
Used with technical field of mine safety.
Background technology
Colliery shock wave Velocity Inversion technology can be used for coal mine rock burst Risk Assessments and early warning, and it is built upon P ripples
Exist between velocity of wave and stress positive correlation basis on, therefore the technology by P ripples Velocity Inversion can on a large scale, high-resolution
Rate obtains the stress distribution feature inside coal and rock, and extensive use is achieved in coal mine rock burst Hazard Assessment warning aspect.
In vibrations operating process of the big gun as excitaton source, typically swash unminding working face side crossheading and install explosive and carry out explosion
Hair, opposite side crossheading lays probe and carries out the reception of shock wave signal, finally by shot point and the P ripple first arrival-times of receiving point
Labeled analysis, inverting obtains wave speed distribution so as to obtain working face stress state.Wherein, P ripples how accurately to be obtained in coal and rock
In propagation time vital effect is played to the accuracy and reliability of inversion result.But in practice of construction, by
Often obscure using the detonator of different hop counts in workman, and the time delay desired value of different hop count detonators is different, even same section
The characteristics of several detonators also meets normal distribution, this causes the P ripples first arrival-time that receiving transducer is recorded to be not fully shock wave
The detonator time delay in propagation time and same hop count in coal and rock is also different, and this is obviously not to vibrations velocity of wave high accuracy inverting
Profit.Therefore, it is the accuracy and reliability of raising shock wave Velocity Inversion result, it is necessary to propose a kind of for colliery vibrations
The method that detonator time delay is eliminated in ripple Velocity Inversion.
The content of the invention
Goal of the invention：Existing realized present in shock wave Velocity Inversion not by explosion the invention aims to overcome
Foot, there is provided a kind of method that detonator time delay is eliminated in shock wave Velocity Inversion for colliery, so as to improve the accurate of inversion result
Property and reliability.
Technical scheme：For achieving the above object, detonator time delay in the shock wave Velocity Inversion for colliery of the invention
The method of elimination, comprises the following steps：
A. the probe for receiving shock wave signal is installed in coal mine work area side crossheading, pacifies in proper order in opposite side crossheading
Blasting charge excites shock wave signal as shot point, and marks the shock wave P ripple first arrival-times of collection of popping one's head in after each shot point explosion
tpi；
B. each shot point is calculated to probe apart from r according to shot point and receiving transducer positioni, i is the numbering of probe in formula；
C. it is based on formula：Reached using each shot point and all received
The P wave ray linear fit methods of probe obtain the time delay t at each shot pointdAnd homogenization P ripple velocities of wave v；
D. the time delay t that each shot point reaches the P wave rays for all receiving probe is eliminatedd, obtain the amendment P ripples on every ray
Then t 'pi=tpi-td。
The shock wave P ripple first break pickup methods of the probes records are：Obtain to blow out using artificial or computer and excite
Effective vibration signal in shock wave P ripple first arrival-times as tpi。
The shot point is to probe apart from riUsing formula：Obtain, x in formulaiIt is suitable to opposite to pop one's head in
The distance of groove projected position and shot point, L is face width.
The time delay t calculated at each shot pointdAnd the approximating method that homogenization P ripple velocities of wave v is utilized is a linear most young waiter in a wineshop or an inn
Multiplication, its formula is as follows：
Wherein, n is the probe number of handmarking's P ripple first arrivals.
T ' when being adapted to after described elimination P ripple time delayspiUsing formula：
It is calculated, 1,2.....n is the numbering of probe in formula.
Beneficial effect：The present invention is prolonged using the detonator that linear least-squares linear fit method has been obtained at each shot point
When, the P ripple first arrival-times t of each probes records is precisely calculated after detonator time delay is eliminatedpi, effectively improve colliery shock wave ripple
The accuracy of fast inversion result.This method is simple, convenient, workable, and rapid effectively acquisition coal mine rock burst danger is commented
Valency and early warning, with wide applicability.
Brief description of the drawings
Fig. 1 is the method flow diagram that the present invention is eliminated for detonator time delay in the shock wave Velocity Inversion of colliery；
Fig. 2 is the colliery active inverting of the method for detonator time delay elimination in the shock wave Velocity Inversion for colliery of the invention
Schematic diagram；
Fig. 3 is that the detonator time delay of the method for detonator time delay elimination in the shock wave Velocity Inversion for colliery of the invention is eliminated
Preceding P ripples first arrival-time and the relation of distance；
Fig. 4 is the revised P ripples of the method for detonator time delay elimination in the shock wave Velocity Inversion for colliery of the invention
The relation of first arrival-time and distance.
Specific embodiment
Embodiment is described further below in conjunction with the accompanying drawings：
As shown in figure 1, the method that detonator time delay is eliminated in the shock wave Velocity Inversion for colliery of the invention, including it is as follows
Step：
A. the probe for receiving shock wave signal is installed in coal mine work area side crossheading, pacifies in proper order in opposite side crossheading
Blasting charge excites shock wave signal as shot point, and marks the shock wave P ripple first arrival-times of collection of popping one's head in after each shot point explosion
tpi；
B. each shot point is calculated to probe apart from r according to shot point and receiving transducer positioni, i is the numbering of probe in formula；
C. it is based on formula：Reached using each shot point and all received
The P wave ray linear fit methods of probe obtain the time delay t at each shot pointdAnd homogenization P ripple velocities of wave v；
D. the time delay t that each shot point reaches the P wave rays for all receiving probe is eliminatedd, obtain the amendment P ripples on every ray
Then t 'pi=tpi-td。
The shock wave P ripple first break pickup methods of the probes records are：Obtain to blow out using artificial or computer and excite
Effective vibration signal in P wave P ripple first arrival-times as tpi。
The shot point is to probe apart from riUsing formula：Obtain, x in formulaiIt is suitable to opposite to pop one's head in
The distance of groove projected position and shot point, L is face width.
The time delay t calculated at each shot pointdAnd the approximating method that homogenization P ripple velocities of wave v is utilized is a linear most young waiter in a wineshop or an inn
Multiplication, its formula is as follows：
Wherein, n is the probe number of handmarking's P ripple first arrivals.
T ' when being adapted to after described elimination P ripple time delayspiUsing formula：
It is calculated, 1,2.....n is the numbering of probe in formula.
Embodiment one：
Shock wave inverting detection coal body inside is excited to answer as shown in Fig. 2 used before certain working face extraction of colliery 3310 and blown out
Power distribution characteristics, shot point 31, track crossheading side arrangement observation station 31 are arranged in haulage gate side.To the waveform of probes records
Data eliminate detonator time delay using the inventive method, and implementation steps are as follows：
(1) wave file to probes records obtains P ripple first arrival-times t using handmarkingpi, (i=1,2 ... 31)；
(2) utilizeFormula calculates each shot point to probe apart from ri, (i=1,2 ... 31)；
Relation of the P ripples obtained according to step (1) and (2) then with distance is as shown in Figure 3, it can be clearly seen that two kinds of sections
The presence of number detonator, and due to the difference of time delay, velocity of wave has obvious grouping feature.
(3) according to following formula, the P ripple first arrival-times t marked using linear least squares fitpiWith apart from riRelation ask
The time delay t gone out at each shot pointdAnd P ripple velocities of wave v：
Wherein,
The detonator time delay and velocity of wave that the step is obtained at all shot points are as shown in table 1.
Detonator time delay and velocity of wave result of calculation at all shot points of table 1
(4) the detonator time delay at each shot point that basis (3) is obtained, the P ripple first arrival-times t of amendment is calculated using following formula
′pi, (i=1,2 ... 31)：
Revised P ripples first arrival-time is as shown in Figure 4 with the relation of distance, it can be seen that velocity of wave packet problem is no longer deposited
.
In order that those skilled in the art more fully understand the technical scheme in the present invention, above-described embodiment is to the present invention
Technical scheme clearly and completely described.Obviously, described embodiment is only a part of embodiment of the invention,
Rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art are not making creative labor
The every other embodiment obtained under the premise of dynamic, should all belong to the scope of protection of the invention.
Claims (5)
1. a kind of method that detonator time delay is eliminated in shock wave Velocity Inversion for colliery, it is characterised in that comprise the following steps：
A. the probe for receiving shock wave signal is installed in coal mine work area side crossheading, installs fried in proper order in opposite side crossheading
Medicine excites shock wave signal as shot point, and marks the shock wave P ripple first arrival-times t of collection of popping one's head in after each shot point explosionpi；
B. each shot point is calculated to probe apart from r according to shot point and receiving transducer positioni, i is the numbering of probe in formula；
C. it is based on formula：Whole receiving transducers are reached using each shot point
P wave ray linear fit methods obtain time delay t at each shot pointdAnd homogenization P ripple velocities of wave v；
D. the time delay t that each shot point reaches the P wave rays for all receiving probe is eliminatedd, obtain amendment P ripples on every ray then
t′pi=tpi-td。
2. the method that detonator time delay is eliminated in the shock wave Velocity Inversion of colliery is used for according to claim 1, it is characterised in that described
The shock wave P ripple first break pickup methods of probes records are：The effective vibrations letter blown out and excite is obtained using artificial or computer
Shock wave P ripple first arrival-times in number are used as tpi。
3. the method that detonator time delay is eliminated in the shock wave Velocity Inversion of colliery is used for according to claim 1, it is characterised in that：It is described
Shot point is to probe apart from riUsing formula：Obtain, x in formulaiFor pop one's head in opposite crossheading projected position with
The distance of shot point, L is face width.
4. the method that detonator time delay is eliminated in the shock wave Velocity Inversion of colliery is used for according to claim 1, it is characterised in that：
The time delay t calculated at each shot pointdAnd the approximating method that homogenization P ripple velocities of wave v is utilized is linear least square, it is public
Formula is as follows：
Wherein, n is the probe number of handmarking's P ripple first arrivals.
5. the method that detonator time delay is eliminated in the shock wave Velocity Inversion of colliery is used for according to claim 1, it is characterised in that：It is described
Elimination P ripple time delays after t ' when being adapted topiUsing formula：
It is calculated, 1,2 ... ..n are the numbering of probe in formula.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611103024.5A CN106772595A (en) | 2016-12-05 | 2016-12-05 | A kind of method that detonator time delay is eliminated in shock wave Velocity Inversion for colliery |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611103024.5A CN106772595A (en) | 2016-12-05 | 2016-12-05 | A kind of method that detonator time delay is eliminated in shock wave Velocity Inversion for colliery |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106772595A true CN106772595A (en) | 2017-05-31 |
Family
ID=58884616
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611103024.5A Pending CN106772595A (en) | 2016-12-05 | 2016-12-05 | A kind of method that detonator time delay is eliminated in shock wave Velocity Inversion for colliery |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106772595A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111413734A (en) * | 2020-04-10 | 2020-07-14 | 徐州弘毅科技发展有限公司 | Calculation method for testing propagation speed and arrival time of underground vibration wave |
| CN112180430A (en) * | 2020-09-23 | 2021-01-05 | 中国矿业大学 | Mine earthquake P wave first arrival identification method in presence of interference signals |
| CN116842311A (en) * | 2023-07-05 | 2023-10-03 | 中交一航局第五工程有限公司 | Blasting detonator delay precision determination method based on air shock wave test |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101132191A (en) * | 2007-10-15 | 2008-02-27 | 北京航空航天大学 | Baseband signal processing method for GNSS receiver |
| CN101201365A (en) * | 2007-07-13 | 2008-06-18 | 北京工业大学 | System and method for measuring and analyzing voltage frequency |
| CN101581789A (en) * | 2009-06-23 | 2009-11-18 | 刘盛东 | Mine working face inter-lane seismic wave CT detection method |
| CN101634590A (en) * | 2009-07-29 | 2010-01-27 | 武汉光迅科技股份有限公司 | Detection and predetermination calculational method for rapid variation light power |
| WO2014088440A1 (en) * | 2012-12-03 | 2014-06-12 | BOCHKAREV, Andrej Vadimovich | A method for applying physical fields of an apparatus in the horizontal end of an inclined well to productive hydrocarbon beds |
| CN104656155A (en) * | 2015-02-10 | 2015-05-27 | 华北科技学院 | Superconductive device capable of increasing exploration depth of nuclear magnetic resonance instrument in coal roadway and signal extraction method |
-
2016
- 2016-12-05 CN CN201611103024.5A patent/CN106772595A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201365A (en) * | 2007-07-13 | 2008-06-18 | 北京工业大学 | System and method for measuring and analyzing voltage frequency |
| CN101132191A (en) * | 2007-10-15 | 2008-02-27 | 北京航空航天大学 | Baseband signal processing method for GNSS receiver |
| CN101581789A (en) * | 2009-06-23 | 2009-11-18 | 刘盛东 | Mine working face inter-lane seismic wave CT detection method |
| CN101634590A (en) * | 2009-07-29 | 2010-01-27 | 武汉光迅科技股份有限公司 | Detection and predetermination calculational method for rapid variation light power |
| WO2014088440A1 (en) * | 2012-12-03 | 2014-06-12 | BOCHKAREV, Andrej Vadimovich | A method for applying physical fields of an apparatus in the horizontal end of an inclined well to productive hydrocarbon beds |
| CN104656155A (en) * | 2015-02-10 | 2015-05-27 | 华北科技学院 | Superconductive device capable of increasing exploration depth of nuclear magnetic resonance instrument in coal roadway and signal extraction method |
Non-Patent Citations (2)
| Title |
|---|
| 刘腾飞: "煤岩声发射信号处理与定位方法的研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
| 钱守一: "HHT在***振动信号处理中的应用研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111413734A (en) * | 2020-04-10 | 2020-07-14 | 徐州弘毅科技发展有限公司 | Calculation method for testing propagation speed and arrival time of underground vibration wave |
| CN111413734B (en) * | 2020-04-10 | 2023-03-10 | 徐州弘毅科技发展有限公司 | Calculation method for testing propagation speed and arrival time of underground vibration wave |
| CN112180430A (en) * | 2020-09-23 | 2021-01-05 | 中国矿业大学 | Mine earthquake P wave first arrival identification method in presence of interference signals |
| CN112180430B (en) * | 2020-09-23 | 2021-08-20 | 中国矿业大学 | Mine earthquake P wave first arrival identification method in presence of interference signals |
| CN116842311A (en) * | 2023-07-05 | 2023-10-03 | 中交一航局第五工程有限公司 | Blasting detonator delay precision determination method based on air shock wave test |
| CN116842311B (en) * | 2023-07-05 | 2024-03-08 | 中交一航局第五工程有限公司 | Blasting detonator delay precision determination method based on air shock wave test |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106772595A (en) | A kind of method that detonator time delay is eliminated in shock wave Velocity Inversion for colliery | |
| CN102135445B (en) | Blasting vibration predicting method | |
| CN102262220B (en) | Positioning method of micro-seismic source or acoustic emission source based on non-linear fitting | |
| Rigby et al. | Angle of incidence effects on far-field positive and negative phase blast parameters | |
| JP2015184144A5 (en) | ||
| CN112012797A (en) | Evaluation method for coal mine impact danger pressure relief effect | |
| Anderson | Blast monitoring: regulations, methods and control techniques | |
| Ouchterlony et al. | Monitoring of large open cut rounds by VOD, PPV and gas pressure measurements | |
| Paswan et al. | Blast vibration and fragmentation control at heavily jointed limestone mine | |
| JP2014222119A (en) | Method and apparatus for detecting detonation | |
| Simioni et al. | Field measurements of snowpack response to explosive loading | |
| CN113217109B (en) | Waveform completion method for rockburst mine microseismic monitoring system | |
| Richards | Elliptical airblast overpressure model | |
| Aloui et al. | Evaluation of ground vibrations and the effect of air blast in open-pit phosphate mines | |
| CN106443645A (en) | Ultrasonic ranging system, and echo value taking method and apparatus | |
| JP2000205794A (en) | Bullet position locator | |
| CN104502913B (en) | Method and device for measuring the distance between coal uncovering roadway and coal seam | |
| JPS61503048A (en) | Geophysical exploration methods for complex ore bodies | |
| Yang et al. | An integrated approach of signature hole vibration monitoring and modeling for quarry vibration control | |
| CN105445378A (en) | Computing method for reflecting overall splitting damage process of argillaceous siltstone | |
| JP7296250B2 (en) | Excavation surface geological evaluation method | |
| JP7315449B2 (en) | Excavation surface geological evaluation method | |
| CN103808803A (en) | Method for detecting completeness of concrete pile body by virtue of wave train figure event image characteristics | |
| CN116842311B (en) | Blasting detonator delay precision determination method based on air shock wave test | |
| Richards | Blast vibration wavefront reinforcement model |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170531 |