CN104236404A - 一种基于***振动测试的现场快速测定毫秒***延时精度的方法 - Google Patents

一种基于***振动测试的现场快速测定毫秒***延时精度的方法 Download PDF

Info

Publication number
CN104236404A
CN104236404A CN201410471008.6A CN201410471008A CN104236404A CN 104236404 A CN104236404 A CN 104236404A CN 201410471008 A CN201410471008 A CN 201410471008A CN 104236404 A CN104236404 A CN 104236404A
Authority
CN
China
Prior art keywords
time
detonator
explosion
big gun
wave
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.)
Granted
Application number
CN201410471008.6A
Other languages
English (en)
Other versions
CN104236404B (zh
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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201410471008.6A priority Critical patent/CN104236404B/zh
Publication of CN104236404A publication Critical patent/CN104236404A/zh
Application granted granted Critical
Publication of CN104236404B publication Critical patent/CN104236404B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于***振动测试的现场快速测定毫秒***延时精度的方法,包括如下步骤:⑴布置炮孔;⑵钻孔、装药完成以后,用导爆索将所有炮孔内的***同时点燃同时起爆;⑶接收***传来的爆轰波,并记录波形数据;⑷利用小波分析方法将实测波形按从低到高不同的频段进行分解,并求得低频段波形的时能密度曲线,以获得每一段***振动到达的精确时间,扣除***应力波在岩石中传播所用时间,得到***的精确延期时间。其优点在于:在一个起爆网络中同时布置多个段号的***,可同时测得每个***的延期时间。在***施工现场中易于应用,能够实时动态地检验毫秒***的精度,根据试验场地的不同地形,可灵活布置起爆网络。

Description

一种基于***振动测试的现场快速测定毫秒***延时精度的方法
技术领域
本发明涉及基于记录***振动到达时间来批量、快速测定毫秒***的延期时间的方法,具体的说是一种基于***振动测试的现场快速测定毫秒***延时精度的方法。
背景技术
随着我国西南地区水电开发的进一步深入以及矿产资源开采深度的日益增加,地下工程将越来越多地应用在国民经济基本建设各个领域,而且大多数的边坡工程涉及到地下洞室的***开挖。近年来,随着工程实践的不断发展,***技术和***振动控制措施取得了巨大的进步,但岩土介质的复杂性给***控制技术的研究带来了较大的困难,如何在复杂的周边环境和施工条件下,使保护对象所受影响在工程许可范围内,快速经济的进行地下洞室的***开挖仍是工程难点。
现在的钻爆工艺已经在飞速地发展,目前较为成熟的有光面***、微差***、预裂***和洞室***等。利用毫秒量级间隔,实现按顺序起爆的方法称为微差***。微差***是一种毫秒级的延期***,使各药包造成的能量场相互影响而产生一系列良好的***效果,是使用最为广泛的***技术,常见于地下掘进***和露天台阶深孔***之中。所以***的精度对***开挖的结果有着至关重要的影响,有必要对开挖所用的***精度进行测试。
已有的***精度测量方法根据测量原理大致可以分为三类:第一类是传统导爆管***精度的测定方法——声光电法,原理是将导爆管***放入声音传导管中,以引爆***产生的光信号为记录开始标志,以******产生的声音信号为记录关闭标志,导爆管***的延期时间就是它们的时间差减去爆轰波在导爆管内的传播时间和声波在传导管中的传播时间,导爆管***延期时间在测定实施过程中,是以光开声关来记录的。第二类则采用特殊的设备——防爆箱来进行试验,原理是将待测试***放在特殊容器中,用一发***作为激发***,激发***起爆后喷出的火焰被光信号接收器接收后,测时仪启动;当被测的***被导爆索引爆后,其喷出的火焰被光讯号接收器接收,测时仪停止,测时仪记录的时间即是被测的某一导爆管***的延期时间。第三类是测震法,原理是通过***振动传感器和记录仪记录导爆管***的***振动波形,经过分析振动波形,找出波形中前后各有一个最大波峰或波谷,这两个波峰或波谷对应的时间就是引爆导爆管和导爆管******时刻,以这两个时刻的时差减去爆轰波在导爆管中传播的时间和声波在传导管中的传播时间,就是导爆管***的延期时间。
另外,***振动号分析是研究***振动控制的基础,传统的信号分析是建立在傅立叶变换(Fourier Transform)的基础之上的,由于傅立叶分析使用的是一种全局变换,即要么完全在时域,要么完全在频域,无法表达信号的时——频局部化性质,而这种性质恰恰是非平稳信号最根本、最关键的性质。因此,傅氏变换从本质上而言只适用于平稳信号而不适用于非平稳信号。
发明内容
本发明的目的是针对现有的技术状况,以***产生的应力波在岩石中传播为依据,提供一种基于***振动测试的现场快速测定毫秒***延时精度的方法,该方法利用常规的***作业,操作方便可靠,效率高,能够在***施工过程中随时测量***的精度。
本发明的核心思想是:根据***应力波在岩石中传播的原理,选择一块岩石完整、节理裂隙不发育的平整场地,钻设按一定顺序排列的系列炮孔,每孔装1-2节***,然后联网起爆,并在不同的爆心距处布置高精度振动检测设备记录振动信号;再利用适合非平稳信号分析的小波分析方法对实测振动波形按从低到高不同的频段进行分解,并求得低频段波形的时能密度曲线,以获得每一段***振动到达的精确时间,扣除***应力波在岩石中传播所用时间,得到***的精确延期时间。在一个起爆网络中同时布置多个段号的***,可同时测得每个***的延期时间。
一种基于***振动测试的现场快速测定毫秒***延时精度的方法,包括如下步骤:
步骤1,布置炮孔,各炮孔所激发的应力波传播路线不重叠;
步骤2,钻孔、装药完成以后,在已有的钻爆参数基础上,用导爆索将所有炮孔内的***同时点燃同时起爆;
步骤3,在不同的爆心距处布置多个振动监测仪,接收***传来的爆轰波,并记录波形数据;
步骤4,利用小波分析方法将实测波形按从低到高不同的频段进行分解,并求得低频段波形的时能密度曲线,以获得每一段***振动到达的精确时间,扣除***应力波在岩石中传播所用时间,得到***的精确延期时间;
步骤4具体包括如下子步骤:
4.1)根据所得到的波形文件,确定***开始时间,记为t0
4.2)找出随后每个波峰开始时刻,记为Tn
4.3)测出每个炮孔中***与监测器的距离Sn
4.4)用测得的数据根据公式(1)计算出各个***的延期时间
t n = T n - t 0 - S n v - - - ( 1 )
式中:t0—开始记录时间,实际t0取0;
tn—***的延期时间,ms;
Tn—波峰开始时刻,ms;
v—***应力波在岩石中的传播速度,m/ms;
Sn—爆心距。
布置炮孔的方式为双排布置或者环形布置或者布置成四边形、五边形。
***应力波在岩石中的传播速度由公式(2)确定:
v = S t - - - ( 2 )
式中:v—***应力波在岩石中的传播速度;
t—应力波在岩石中传播时间,ms;
S—爆心距,m。
对任意能量有限的函数f(t),关于ψ(t)的连续小波变换定义如下:
W f ( a , b ) = < f , &psi; a , b > = | a | - 1 2 &Integral; R f ( t ) &psi; &OverBar; ( t - b a ) dt - - - ( 3 )
上式中的a、b分别为伸缩因子和平移因子。根据Moyal内积定理,有下式成立:
1 C &psi; &Integral; R da a 2 &Integral; R | W f ( a , b ) | 2 ab = &Integral; R | f ( t ) | 2 dt - - - ( 4 )
上式表明,小波变换幅度平方的积分和被分析信号的能量成正比。根据能量密度的概念,式(4)可以写成:
&Integral; R | f ( t ) | 2 dt &Integral; R E ( b ) db - - - ( 5 )
式中:
E ( b ) = 1 C &psi; &Integral; R 1 a 2 | W f ( a , b ) | 2 da - - - ( 6 )
小波变换中,尺度a在一定意义上对应于频率ω,因此(6)式给出了信号所有频带的能量随时间b的分布情况,称为时能密度函数。实际应用中,可以通过改变上式的积分上、下限,使积分区间落在待分析信号的某频率范围内,从而得到该频带内信号能量密度随时间的分布特征。
如将某次***作为一个***来考虑,则每一段***的起爆就是向***输入能量的过程,每一段***起爆必然引起***内能量密度的改变。因此,可以适当选取积分上、下限,按(6)式计算单段围岩振动信号主频段内的能量密度,并画出其时-能密度图。根据图中出现的突峰位置即可得到各微差段***荷载所诱发的振动的到达时刻,从而得各段***的实际起爆延迟时间。
本发明一种基于***振动测试的现场快速测定毫秒***延时精度的方法的优点是:在一个起爆网络中同时布置多个段号的***,可同时测得每个***的延期时间。较传统的***精度测量方法操作简单,效率高,而且分析结果可靠、安全,成本低廉。在***施工现场中易于应用,能够实时动态地检验毫秒***的精度,根据试验场地的不同地形,可灵活布置起爆网络。
附图说明
图1为测试岩石中纵波速度时炮孔与监测器的位置示意图。
图2为***信号滤波图。
图3为时能密度图。
图4为炮孔双排布置时炮孔与监测器的位置示意图。
图5为炮孔双排布置时监测器记录的振动波形图。
图6为炮孔环形布置时炮孔与监测器的位置示意图。
图7为炮孔环形布置时监测器记录的振动波形图。
具体实施方式
下面结合附图,对本发明进行进一步说明。
实施例一
确定测试场地的纵波速度。在测试场地中央布置一个炮孔,孔内装填段号为MS1的***,在距离炮孔不同距离处布置三个***振动监测仪A、B、C,用来接收***传来的应力波,监测仪之间的夹角为120°。引爆***,测出应力波到达监测设备的时间t以及***与监测仪的距离S,由于MS1段的***延期时间为0ms,利用公式(2)即可求出***应力波在岩石中的传播速度。
在某宽阔的平整的完整无裂隙岩体中,布置两排炮孔,排距为5m,试验时每排炮孔取5个;垂直钻孔,孔径取70mm,孔深为1m,炮孔之间的距离为4m。
在两排炮孔之间布置三个振动监测传感器,监测器之间的距离为2m。将炮孔接成如图4所示的由MS1、MS3、MS5、MS7、MS9、MS11、MS13、MS15、MS17、MS19以及***组成的起爆网络。
按上述联网进行起爆,记录***发生时刻,此后导爆管***由于延期时间不同而依次***,***产生的应力波传至监测点被监测设备接收。将测得的数据进行处理,最终得到的波形图象如图5所示,图象中包含多个波段,理论上波段个数与***数量一致。
用滤波的方法分离出低频信号,再用小波分析及时能密度方法确定每个波段中波峰开始时刻,根据图中出现的突峰位置即可得到各微差段***荷载所诱发的振动的到达时刻,从而得各段***的实际起爆延迟时间。
具体实施过程如下:
1)确定***开始时间,记为t0
2)找出随后每个波峰开始时刻,记为Tn
3)测出每个炮孔中***与监测器的距离Sn
4)用测得的数据根据公式(1)计算出各个***的延期时间。
为了得到更加准确的结果,我们在测试时,在炮孔不同位置布置3个监测器,分别测出每个***的延期时间,对结果进行比较分析。
实施例二
确定测试场地的纵波速度。在测试场地中央布置一个炮孔,孔内装填段号为MS1的***,在距离炮孔不同距离处布置三个***振动监测仪A、B、C,用来接收***传来的应力波,监测仪之间的夹角为120°。引爆***,测出应力波到达监测设备的时间t以及***与监测仪的距离S,由于MS1段的***延期时间为0ms,利用公式(2)即可求出***应力波在岩石中的传播速度。
在某宽阔的平整的完整无裂隙岩体中,布置一圈环形炮孔,初次试验时炮孔取9个;垂直钻孔,孔径取70mm,孔深为1m,炮孔之间的夹角为22.5°。
在炮孔连线所在的圆形的圆心处布置一个振动监测传感器,记为A;同时,在圆心附近布置另外两个监测器,分别记为B、C,监测器之间的距离为2m。将炮孔接成如图6所示的由MS1、MS3、MS5、MS7、MS9、MS11、MS13、MS15、MS17以及***组成的起爆网络。
按上述联网进行起爆,记录***发生时刻,此后导爆管***由于延期时间不同而依次***,***产生的应力波传至监测点被监测设备接收。最终得到的波形图象如图7所示,图象中包含多个波段,理论上波段个数与***数量一致。
用滤波的方法分离出低频信号,再用小波分析及时能密度方法确定每个波段中波峰开始时刻,记为Tn,测出每个导爆管***与监测器之间的距离Sn,根据公式(1)计算出各个***的延期时间。
为了得到更加准确的结果,我们在测试时,在炮孔不同位置布置3个监测器,分别测出每个***的延期时间,对结果进行比较分析。

Claims (2)

1.一种基于***振动测试的现场快速测定毫秒***延时精度的方法,其特征在于包括如下步骤:
(1)布置炮孔,各炮孔所激发的应力波传播路线不重叠;
(2)钻孔、装药完成以后,在已有的钻爆参数基础上,用导爆索将所有炮孔内的***同时点燃同时起爆;
(3)在不同的爆心距处布置多个振动监测仪,接收***传来的爆轰波,并记录波形数据;
(4)利用小波分析方法将实测波形按从低到高不同的频段进行分解,并求得低频段波形的时能密度曲线,以获得每一段***振动到达的精确时间,扣除***应力波在岩石中传播所用时间,得到***的精确延期时间;具体包括如下子步骤:
(4.1)根据所得到的波形文件,确定***开始时间,记为t0
(4.2)找出随后每个波峰开始时刻,记为Tn
(4.3)测出每个炮孔中***与监测器的距离Sn
(4.4)用测得的数据根据公式(1)计算出各个***的延期时间
t n = T n - t 0 - S n v - - - ( 1 )
式中:t0—开始记录时间,实际t0取0;
tn—***的延期时间,ms;
Tn—波峰开始时刻,ms;
v—***应力波在岩石中的传播速度,m/ms;
Sn—爆心距。
2.如权利要求1所述的基于***振动测试的现场快速测定毫秒***延时精度的方法,其特征在于:所述布置炮孔的方式为双排布置或者环形布置或者布置成四边形、五边形。
CN201410471008.6A 2014-09-16 2014-09-16 一种基于***振动测试的现场快速测定毫秒***延时精度的方法 Active CN104236404B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410471008.6A CN104236404B (zh) 2014-09-16 2014-09-16 一种基于***振动测试的现场快速测定毫秒***延时精度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410471008.6A CN104236404B (zh) 2014-09-16 2014-09-16 一种基于***振动测试的现场快速测定毫秒***延时精度的方法

Publications (2)

Publication Number Publication Date
CN104236404A true CN104236404A (zh) 2014-12-24
CN104236404B CN104236404B (zh) 2016-03-02

Family

ID=52224963

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410471008.6A Active CN104236404B (zh) 2014-09-16 2014-09-16 一种基于***振动测试的现场快速测定毫秒***延时精度的方法

Country Status (1)

Country Link
CN (1) CN104236404B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106092311A (zh) * 2016-06-08 2016-11-09 武汉大学 一种基于实测高频振动***振动信号的起爆网路质量评价方法
CN108008349A (zh) * 2017-12-07 2018-05-08 西安近代化学研究所 爆心定位方法
CN109827481A (zh) * 2019-01-18 2019-05-31 华北科技学院 一种微差干扰降振***中最优微差时间确定方法及装置
CN112130207A (zh) * 2020-09-25 2020-12-25 中国科学院武汉岩土力学研究所 一种基于球形装药条件下由地面振动计算地下振动的方法
WO2022198840A1 (zh) * 2021-03-23 2022-09-29 山东科技大学 模拟隧道工程***的脉冲声源及***噪声测试装置和方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
PL202555B1 (pl) * 2004-01-23 2009-07-31 O & Sacute Rodek Badawczo Rozw Sposób i urządzenie do testowania zapalnika czasowo-ciśnieniowego
CN102927862A (zh) * 2012-10-25 2013-02-13 长江水利委员会长江科学院 一种快速测量***延时精度的方法
CN103134398A (zh) * 2013-02-08 2013-06-05 前进民爆股份有限公司 一种力学测量数码电子***起爆延时的方法及测时仪

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
PL202555B1 (pl) * 2004-01-23 2009-07-31 O & Sacute Rodek Badawczo Rozw Sposób i urządzenie do testowania zapalnika czasowo-ciśnieniowego
CN102927862A (zh) * 2012-10-25 2013-02-13 长江水利委员会长江科学院 一种快速测量***延时精度的方法
CN103134398A (zh) * 2013-02-08 2013-06-05 前进民爆股份有限公司 一种力学测量数码电子***起爆延时的方法及测时仪

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘翼等: "***振动法测导爆管***延期时间", 《***器材》, vol. 35, no. 01, 28 February 2006 (2006-02-28), pages 15 - 18 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106092311A (zh) * 2016-06-08 2016-11-09 武汉大学 一种基于实测高频振动***振动信号的起爆网路质量评价方法
CN106092311B (zh) * 2016-06-08 2019-01-04 武汉大学 一种基于实测高频振动***振动信号的起爆网路质量评价方法
CN108008349A (zh) * 2017-12-07 2018-05-08 西安近代化学研究所 爆心定位方法
CN109827481A (zh) * 2019-01-18 2019-05-31 华北科技学院 一种微差干扰降振***中最优微差时间确定方法及装置
CN109827481B (zh) * 2019-01-18 2021-06-25 华北科技学院 一种微差干扰降振***中最优微差时间确定方法及装置
CN112130207A (zh) * 2020-09-25 2020-12-25 中国科学院武汉岩土力学研究所 一种基于球形装药条件下由地面振动计算地下振动的方法
CN112130207B (zh) * 2020-09-25 2021-07-20 中国科学院武汉岩土力学研究所 一种基于球形装药条件下由地面振动计算地下振动的方法
WO2022198840A1 (zh) * 2021-03-23 2022-09-29 山东科技大学 模拟隧道工程***的脉冲声源及***噪声测试装置和方法

Also Published As

Publication number Publication date
CN104236404B (zh) 2016-03-02

Similar Documents

Publication Publication Date Title
CN104236404B (zh) 一种基于***振动测试的现场快速测定毫秒***延时精度的方法
CN102829899B (zh) 一种适用于深埋圆形隧洞的围岩应力快速测算方法
Bogdanoff Vibration measurements in the damage zone in tunnel blasting
CN107238538B (zh) 弱***诱导的应变型岩爆现场模拟试验方法
CN104197795B (zh) 一种基于***声压测试的现场***延时精度测试方法
CN102829678B (zh) 一种错相减震的***方法
US7370513B2 (en) System and method for monitoring features of a blast
CN103344156B (zh) ***工程中的盲炮识别方法
CN101858715B (zh) 一种***中药包发生拒爆的识别及其拒爆药包的定位方法
Leidig et al. Quantification of ground vibration differences from well-confined single-hole explosions with variable velocity of detonation
Hinzen Comparison of seismic and explosive energy in five smooth blasting test rounds
Simioni et al. Field measurements of snowpack response to explosive loading
CN105865280A (zh) 一种优化设计岩石匹配现场混装乳化***的方法
CN103256874A (zh) 利用***振动波检测盲炮的方法
CN105973456A (zh) 一种现场测量深孔******与岩石匹配关系的方法
CN114577386B (zh) 炮孔内电子***所受真实冲击波荷载测试方法和测试***
CN107401958A (zh) 一种利用电子芯片******振动波分析识别盲炮的方法
RU2407984C2 (ru) Способ определения относительной работоспособности промышленного взрывчатого вещества
Battison et al. Reducing crest loss at Barrick Cowal gold mine
Bilgin et al. Assessment of explosive performance by detonation velocity measurements in Turkey
CN105928601B (zh) 一种测量露天矿山爆区下方台阶岩体声速的方法
Yang et al. An integrated approach of signature hole vibration monitoring and modeling for quarry vibration control
Smith et al. Detonation wave propagation in underground mine entries
Garaliu-Busoi et al. Ensuring the seismic protection of the overground objectives in the neighboring area of industrial cement producers quarries
Pytlik et al. Determining the timing accuracy of electronic detonators using different methods

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant