CN104197795B - 一种基于***声压测试的现场***延时精度测试方法 - Google Patents

一种基于***声压测试的现场***延时精度测试方法 Download PDF

Info

Publication number
CN104197795B
CN104197795B CN201410471316.9A CN201410471316A CN104197795B CN 104197795 B CN104197795 B CN 104197795B CN 201410471316 A CN201410471316 A CN 201410471316A CN 104197795 B CN104197795 B CN 104197795B
Authority
CN
China
Prior art keywords
detonator
time
wave
blast
sound 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.)
Active
Application number
CN201410471316.9A
Other languages
English (en)
Other versions
CN104197795A (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 CN201410471316.9A priority Critical patent/CN104197795B/zh
Publication of CN104197795A publication Critical patent/CN104197795A/zh
Application granted granted Critical
Publication of CN104197795B publication Critical patent/CN104197795B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于***声压测试的现场***延时精度测试方法,将***按串联或并联的方式布置在一个圆周上,并在圆心处不同高度布置多台声压传感器,起爆后获得各******激发的空气冲击波压力曲线或声压曲线;再利用适合非平稳信号分析的小波分析方法对实测波形进行分析,并求得波形的时能密度曲线,以获得每一段声波到达的精确时间,扣除声波在空气中传播所用时间,得到***的精确延期时间。其优点是:在一个起爆网络中同时布置多个段号的***,可同时测得每个***的延期时间。操作简单,效率高,而且分析结果可靠,成本低廉。在***施工现场中易于应用,能够实时动态地检验毫秒***的精度,根据试验场地的不同地形,可灵活布置起爆网络。

Description

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

Claims (1)

1.一种基于***声压测试的现场***延时精度测试方法,其特征在于包括如下步骤:
⑴布置***,各***所激发的声波传播路线不重叠;所述***的布置方式为8个***均布在圆周上,***之间为并联或者串联;
⑵在已有的钻爆参数基础上,用导爆索引爆所有的***;
⑶在不同的爆心距处布置多个声压监测仪,接收***传来的声波信号,并记录波形数据;
⑷用小波分析及时能密度方法确定每个波段中波峰开始时刻,根据突峰位置得到各微差段***所产生的声波的到达时刻,从而获得各段***的实际起爆延迟时间;具体包括如下子步骤:
(4.1)根据所得到的波形文件,确定***开始时间,记为T0
(4.2)找出随后每个波峰开始时刻,记为Tn
(4.3)测出每个炮孔中***与监测器的距离Sn
(4.4)用测得的数据根据下述公式计算出各个***的延期时间
当***串联布置时
t 1 = T 1 - T 0 - S 1 v ( n = 1 ) t n = T n - T 0 - S n v - Σ i = 1 n - 1 t i ( n > 1 ) - - - ( 1 )
式中:tn—***的延期时间,ms;
T0—记录开始时间,实际T0取0;
Tn—每段波峰开始时刻,ms;
v—声波在空气中的传播速度,取340m/s;
Sn—***到监测仪的距离;
当***并联布置时
t n = T n - T 0 - S n v - Δ t - - - ( 2 )
式中:tn—***的延期时间,ms;
T0—记录开始时间,实际T0取0;
Tn—每段波峰开始时刻,ms;
v—声波在空气中的传播速度,取340m/s;
Sn—***到监测仪的距离;
Δt—爆轰波在导爆管中的传播时间。
CN201410471316.9A 2014-09-16 2014-09-16 一种基于***声压测试的现场***延时精度测试方法 Active CN104197795B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410471316.9A CN104197795B (zh) 2014-09-16 2014-09-16 一种基于***声压测试的现场***延时精度测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410471316.9A CN104197795B (zh) 2014-09-16 2014-09-16 一种基于***声压测试的现场***延时精度测试方法

Publications (2)

Publication Number Publication Date
CN104197795A CN104197795A (zh) 2014-12-10
CN104197795B true CN104197795B (zh) 2016-08-17

Family

ID=52083123

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410471316.9A Active CN104197795B (zh) 2014-09-16 2014-09-16 一种基于***声压测试的现场***延时精度测试方法

Country Status (1)

Country Link
CN (1) CN104197795B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104457457B (zh) * 2014-11-27 2016-08-24 安徽理工大学 ***延期体延时测试新方法
CN108981507B (zh) * 2018-06-26 2020-10-02 湖北省水利水电规划勘测设计院 一种基于高速摄像机的***延时精度现场快速测量方法
CN113834392B (zh) * 2021-10-14 2023-01-10 郑州大学 一种基于水击波信号的电子***水下***威力测试装置
CN115406321A (zh) * 2022-09-26 2022-11-29 三峡大学 一种基于空气冲击波的电子***拒爆识别方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101988814A (zh) * 2009-03-30 2011-03-23 中水东北勘测设计研究有限责任公司 控制***中毫秒延期时间的检验方法
CN102927862A (zh) * 2012-10-25 2013-02-13 长江水利委员会长江科学院 一种快速测量***延时精度的方法
EP2407748A3 (de) * 2010-07-15 2014-12-31 TDW Gesellschaft für verteidigungstechnische Wirksysteme mbH Messverfahren für einen Penetratorzünder

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101988814A (zh) * 2009-03-30 2011-03-23 中水东北勘测设计研究有限责任公司 控制***中毫秒延期时间的检验方法
EP2407748A3 (de) * 2010-07-15 2014-12-31 TDW Gesellschaft für verteidigungstechnische Wirksysteme mbH Messverfahren für einen Penetratorzünder
CN102927862A (zh) * 2012-10-25 2013-02-13 长江水利委员会长江科学院 一种快速测量***延时精度的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"***振动法测导爆管***延期时间";刘翼等;《***器材》;20060228;第35卷(第1期);第16页 *

Also Published As

Publication number Publication date
CN104197795A (zh) 2014-12-10

Similar Documents

Publication Publication Date Title
CN104236404B (zh) 一种基于***振动测试的现场快速测定毫秒***延时精度的方法
CN104197795B (zh) 一种基于***声压测试的现场***延时精度测试方法
CN102829899B (zh) 一种适用于深埋圆形隧洞的围岩应力快速测算方法
CN102829678B (zh) 一种错相减震的***方法
CN106874586A (zh) 一种考虑地应力瞬态卸荷的***振动安全判定方法
CN104181159A (zh) 一种网格式爆轰波形扫描测试方法
Gitterman Secondary shock features for large surface explosions: Results from the Sayarim Military Range, Israel and other experiments
CN105928425A (zh) 一种导爆索爆速测量装置及测量方法
US7370513B2 (en) System and method for monitoring features of a blast
CN110578516B (zh) 一种超临界co2相变脉冲下岩体冲击致裂损伤测试方法
CN103344156B (zh) ***工程中的盲炮识别方法
CN108981507B (zh) 一种基于高速摄像机的***延时精度现场快速测量方法
JP5757474B2 (ja) 断層探査装置および断層探査方法
Simioni et al. Field measurements of snowpack response to explosive loading
Leidig et al. Quantification of ground vibration differences from well-confined single-hole explosions with variable velocity of detonation
JP7030502B2 (ja) 岩盤評価方法
CA2807885A1 (en) Seismic energy sources and methods of use
CN105973456A (zh) 一种现场测量深孔******与岩石匹配关系的方法
CN107401958A (zh) 一种利用电子芯片******振动波分析识别盲炮的方法
CN102927862A (zh) 一种快速测量***延时精度的方法
CN112945034A (zh) 基于水下冲击波实测波形的电子***延时精度测试方法
CN102495191A (zh) 一种埋地卧式油罐燃爆事故危害测试装置
GB1588693A (en) Method of monitoring underground processes
Yang et al. An integrated approach of signature hole vibration monitoring and modeling for quarry vibration control
CN105928601B (zh) 一种测量露天矿山爆区下方台阶岩体声速的方法

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