专利名称：爆破振动预测方法
技术领域：
本发明涉及爆破振动安全领域，尤其涉及一种可在爆前预测目标点爆破振动波形 的爆破振动预测方法。
背景技术：
随着我国基础建设和矿山开采的迅猛发展，爆破作为土岩开挖的一种必要手段所 面临的环境越来越复杂，如何采用科学、有效的措施最大限度保证爆破安全，尽量降低爆破 振动对周边环境的影响，成为顺利实施爆破施工、构筑和谐社会的一项迫切任务。由此，对 爆破振动的评估及分析也成为人们日益重视的问题。目前，通常采用萨道夫斯基公式Γ= K f 对爆破振动效果进行预测评估。式
I R J
中，V代表介质质点振动速度的峰值，单位为CHMS—1 ；R代表爆心距离，单位为m ；Q代表同一 时刻起爆的炸药量，单位为kg ;Κ、α分别代表与爆破条件、地形地质条件等有关的系数和 衰减指数。在工程应用中，将相应数值代入萨道夫斯基公式就可计算出介质质点振动速度 的峰值V，进而再利用峰值V的大小对爆破振动进行安全评估。但采用上述公式进行爆破振动预测评估存在以下缺陷1、上述公式中K和α的值都是根据爆破现场的实际情况靠经验确定的。若想得 到比较可信的K和α的值，通常需要利用该爆破现场的多次爆破振动测试数据、采用回归 分析法才能得到相应分析结果，但是通常回归分析的相关系数较低，采用上述公式计算爆 破振动速度峰值V，只能从概率原理上进行预测。因此，预测结果的可信度和准确度难以把握。2、公式中炸药量Q的值通常取为起爆网路中同段装药总量的最大值，R代表最大 药量爆破药包对应到预测目标点的距离。而实际工程应用中，由于起爆器材自身延期精度 的限制，同段的起爆器材无法精确同步起爆；同时，同段起爆的多个炮孔至目标点的距离也 各不相同。因此，采用这种方式确定Q和R的值与公式的基本假设有较大差异，得出的预测 结果自然准确度降低。3、通常评估爆破振动效果时除需考虑质点振动速度的峰值V外，还需考虑爆破振 动频率和爆破持续时间，只有综合考虑到这三方面的因素，才能得到比较可信的评估结果。 而上述公式仅仅给出了对质点振动速度的峰值V的预测方法，无法得知爆破振动频率和爆 破持续时间，因此采用上述公式对爆破振动进行分析评估不够全面。

发明内容
本发明的技术目的在于解决上述现有技术的缺陷，提供一种爆破振动预测方法， 该方法基于单孔爆破试验获得单孔爆破信息，再利用振动波叠加原理得到预测爆破振动波 形。采用该方法得到的振动波形的预测结果能更加客观地体现出场地的地形地质条件对爆 破振动的影响，并能得到预测目标点的完整振动波形。采用本发明的方法预测质点振速峰值比传统的萨道夫斯基经验公式更加准确可靠。本发明所提供的爆破振动预测方法是按以下步骤进行的首先，在爆破区域附近确定一基准爆破孔；其次，在位置坐标为(U(l，v0, r0)的预测目标点与上述基准爆破孔之间的连线上 选取至少一个参考测振点，并在预测目标点和各参考测振点处分别安置可同步触发的测振 仪；再次，起爆基准爆破孔，以获得预测目标点和各参考测振点的爆破振动波形；最后，根据基准爆破孔的爆破信息以及爆破区域内各炮孔的布孔信息，执行爆破 振动波形预测进程，输出对预测目标点的振动波形。采用本发明的技术方案对振动波形进行预测时，先选取一基准爆破孔用于获取基 于场地条件的振动波形，作为原始分析数据。起爆基准爆破孔后，根据基准爆破孔的爆破信 息以及爆破区域内各炮孔的布孔信息，执行爆破振动波形预测进程，即可得到预测目标点 的预测爆破振动波形。这样得到的振动波形能综合体现质点振动速度的峰值、爆破振动频 率以及爆破持续时间等评价爆破振动效果需要的信息，进而可对爆破振动效果进行更全面 的评估。除此之外，利用本发明的技术方案进行爆破振动评估时，只需在爆破现场进行一次 单孔爆破试验、依据在预测目标点和参考测振点分别测得的数据就可以得到可信的分析结 果，爆破现场操作简单，也有利于提高测振效率。本方法执行爆破振动波形预测进程对从各 测振点测得的数据进行分析处理，这一数据分析过程几乎不受人为因素的影响，进一步提 高了分析结果的准确性和可信性，也为更好地优化调整爆破设计提供了有效手段。本发明爆破振动预测方法中的爆破振动波形预测进程，依据基准爆破孔的爆破信 息以及爆破区域内各炮孔的布孔信息进行预测。其中，基准爆破孔的爆破信息包括基准爆 破孔的位置坐标(Χο,Κ,Ζο)和装药量％，炮孔的布孔信息包括炮孔的总数η、以及各炮孔的 位置坐标(Xi, Yi, Zi) (i = 1，2，...，η)、装药量= 1，2，...，η)和起爆延期时间、(i
1，2，...，η) ο作为本发明的一种优选方案，上述基准爆破孔的爆破信息还可进一步包括炮孔深 度Htl，各炮孔的布孔信息还可进一步包括炮孔深度Hi (i = 1，2，...，n)。考虑炮孔的深度， 更有利于计算结果的准确性。爆破振动波形预测进程可按照以下步骤进行步骤一，分别读取从各测振仪分别测得的各测振点的爆破振动数据，并计算得到 基准爆破孔在预测目标点的振动波形Mtci)和爆破区域当地的地震波波速Cv ；步骤二，读取基准爆破孔的爆破信息；读取各炮孔的布孔信息；步骤三，依据基准爆破孔的爆破信息计算地震波从基准爆破孔传播到预测目标点 所经过的距离Ltl ；依据各炮孔的布孔信息分别计算地震波从各炮孔传播到预测目标点将经 过的距离Li ； 步骤四，依据公式f(。)=
αιβο
V
C τ \α L0
ν A j
，依次计算各炮孔的爆破振动波
L -L
形f(ti)，其中，所述α为爆破振动衰减指数；依据公式乃=A+ 1 c °，依次计算各炮
V孔的爆破振动波形f (、)在预测目标点的偏移叠加延期时间Ti (i = 1，2，. . .，η)；步骤五，基于振动波叠加原理，在本发明中体现为爆破振动波叠加原理，依据公式
η
F(0 = Ef^)，计算得到预测目标点的振动波形F(t)； /=1步骤六，结束本爆破振动波形预测进程。上述爆破振动波形预测进程中，计算出各个炮孔的起爆振动波形Mti)以及各炮 孔的起爆振动波形在预测目标点的偏移叠加延期时间Ti后，就可依据振波叠加原理计算出 各个炮孔作用于预测目标点时振动波形的叠加效果，从而得以预测预测目标点处振动速度 的峰值、爆破振动频率和爆破持续时间等评价指标，而且还能得知峰值出现的时刻，也就是 爆破发生后、振动对预测目标点影响最明显的时刻。这就可以准确地把握爆破振动的过程， 更加全面地预测爆破振动效果，从而为优化爆破网路设计提供参考依据。上述地震波从基准爆破孔传播到预测目标点所经过的距离Ltl可采用公式
L0 =V(^-Mo)2 +(^o -V0)2 +( -r-ο if ^ ；地震波从各炮孔传播到预测目标点将经过的距离Li可采用公式 L1 = ^(Xi-U0)2 +(yt -V0)2+(Zj -r。)2 计算得到。进一步考虑基准爆破孔的炮孔深度Htl和各炮孔的炮孔深 度HiG = 1，2，. . .，η)对地震波传播距离的影响，则距离Ltl可采用公式
L0 =^x0 -U0)2+(y0 -V0)2 +(z0-r0 +ψ)2 计算得到，距离 Li 可采用公式
Ii=^xi-U0)2+Cyi -V0)2 +(Zi-r0计算得到。这样算得的距离就能更准确地反映
地震波传播的实际距离，从而更有利于提高预测结果的准确性。作为本发明中爆破振动波形预测进程的优选方案，在执行上述步骤一时，可采用 低通滤波的方式读取各测振仪测得的振动数据，这样有利于抑制高频噪音，消除环境对采 集数据的影响，从而有利于提高振动数据的准确性。处理原始振动数据时采用的低通滤波 阈值优选取为500Hz以内，以300Hz以内为最佳。


图1-1为本发明单孔爆破试验中基准爆破孔起爆后传播到图5所示测振点302处 的振动波形图；图1-2为采用本发明爆破振动预测方法得出的测振点302处的预测波形图；图1-3为实际爆破后测振点302处的实际波形图；图2-1为本发明单孔爆破试验中基准爆破孔起爆后传播到图5所示测振点301处 的振动波形图；图2-2为采用本发明爆破振动预测方法得出的测振点301处的预测波形图；图2-3为实际爆破后测振点301处的实际波形图；图3-1为本发明单孔爆破试验中基准爆破孔起爆后传播到图5所示预测目标点 300处的振动波形图3-2为采用本发明爆破振动预测方法得出的预测目标点300处的预测波形图；图3-3为实际爆破后预测目标点300处的实际波形图；图4为振波叠加原理的示意图；图5为采用本发明的方法进行单孔爆破试验时的实施示意图。
具体实施例方式下面结合附图和具体实施方式
对本发明的技术方案做进一步详细说明。本发明所提供的爆破振动预测方法是按以下步骤进行的，参见图5所示实施示意 图首先，在爆破区域100附近确定基准爆破孔200的位置以及该孔200的装药量。基 准爆破孔200位置的选取应在确保该孔200的起爆不对爆破区域100内各炮孔构成影响的 前提下，尽可能地接近爆破区域100。其次，在位置坐标为(U(l，v0, r0)的预测目标点300与上述基准爆破孔200之间的 连线上选取至少一个参考测振点，例如图5中所示的测振点301和测振点302。在预测目 标点300和各参考测振点处分别安置可同步触发的测振仪。参考测振点301/302应选取在 基准爆破孔200起爆产生的地震波传播到预测目标点300的必然路径上，并且参考测振点 301/302与基准爆破孔200之间的距离应使得参考测振点处的测振仪不被损坏。再次，起爆基准爆破孔200，获得预测目标点300、参考测振点301以及参考测振点 302处的爆破振动波形，分别如图3-1、图2-1和图1-1所示。最后，根据基准爆破孔200的爆破信息以及爆破区域100内各炮孔的布孔信息，执 行爆破振动波形预测进程，输出预测目标点300的振动波形。采用本发明的实现方式对振动波形进行预测时，先选取一基准爆破孔200用于获 取振动波形的原始分析数据。起爆基准爆破孔200后，结合基准爆破孔200的爆破信息以 及爆破区域100内各炮孔的布孔信息，执行爆破振动波形预测进程，即可得到预测目标点 的预测爆破振动波形。这样得到的振动波形能综合体现质点振动速度的峰值、爆破振动频 率以及爆破持续时间等评价爆破振动效果需要的信息，进而可对爆破振动效果进行更全面 的评估。除此之外，利用本发明的技术方案进行评估时，只需在爆破现场进行一次单孔爆破 试验、依据在预测目标点和参考测振点分别测得的数据就可以得到可信的分析结果，爆破 现场操作简单，也有利于提高测振效率。本方法执行爆破振动波形预测进程对从各测振点 测得的数据进行分析处理，这一数据分析过程几乎不受人为因素的影响，进一步提高了分 析结果的准确性和可信性，也为更好地优化调整爆破设计提供了有效手段。本发明爆破振动预测方法中的爆破振动波形预测进程，依据基准爆破孔200的爆 破信息以及爆破区域100内各炮孔(例如孔A105、孔A112、孔A113等)的布孔信息进行预 测。其中，基准爆破孔200的爆破信息包括基准爆破孔的位置坐标( ，％，Ζ(Ι)和装药量仏， 炮孔的布孔信息包括炮孔的总数η、以及各炮孔的位置坐标(Xi，yi，Zi) (i = 1，2，...，n)、装 药量QiG = 1,2,... ,η)和起爆延期时间、(1 = 1，2，...，η)。由于爆破施工时，炮孔实际 深度往往与设计值有一定偏差，因此，各炮孔的炮孔深度需完成打孔后逐孔实际测量方能 得到。在执行本发明中的爆破振动波形预测进程时，若炮孔深度与该炮孔距预测目标点的 距离相比小一个或多个数量级，或各个炮孔深度相差不大时，可以不考虑炮孔深度的影响，这样有利于简化现场操作步骤，提高工作效率，同时，仍然可以得到较现有技术更加准确可 靠的计算结果，满足精度要求。考虑炮孔深度对振动波形的影响，上述基准爆破孔200的爆破信息进一步包括炮 孔深度Htl，各炮孔的布孔信息进一步包括炮孔深度HiG = 1，2，...，n)，这样就更加周全地 考虑了基准爆破孔200的爆破信息以及各炮孔的布孔信息，从而更有利于计算结果的准确 性。爆破振动波形预测进程可按照以下步骤进行步骤一，分别读取各测振仪分别测得的各测振点的爆破振动数据，并计算得到基 准爆破孔200在预测目标点300的振动波形f (、)和爆破区域100当地的地震波波速Cv ；步骤二，读取基准爆破孔的爆破信息；读取各炮孔的布孔信息；步骤三，依据基准爆破孔的爆破信息计算地震波从基准爆破孔200传播到预测目 标点300所经过的距离Ltl ；依据各炮孔的布孔信息，分别计算地震波从各炮孔传播到预测 目标点将经过的距离L^i = 1,2,... ,η)；步骤四，依据公式
权利要求
1.一种爆破振动预测方法，其特征在于 首先，在爆破区域附近确定一基准爆破孔；其次，在位置坐标为(Utl, v0, r0)的所述预测目标点与所述基准爆破孔之间的连线上选 取至少一个参考测振点，并在所述预测目标点和各所述参考测振点处分别安置可同步触发 的测振仪；再次，起爆所述基准爆破孔，以获得所述预测目标点和各所述参考测振点的爆破振动 波形；最后，结合所述基准爆破孔的爆破信息以及所述爆破区域内各炮孔的布孔信息，执行 爆破振动波形预测进程，输出对所述预测目标点的所述振动波形进行预测的结果。
2.按照权利要求1所述的预测方法，其特征在于所述爆破信息包括所述基准爆破孔的位置坐标( ，10，Z0)和装药量仏； 所述布孔信息包括所述炮孔的总数n，以及各所述炮孔的位置坐标(Xi，yi Zi) (i = 1， 2，...，η)、装药量 & (i = 1,2, ... ,η)和起爆延期时间、(1 = 1，2，...，η)。
3.按照权利要求2所述的预测方法，其特征在于 所述爆破信息还包括所述基准爆破孔的炮孔深度Η0。
4.按照权利要求3所述的预测方法，其特征在于所述布孔信息还包括各所述炮孔的炮孔深度Hi (i = 1，2，. . .，η)。
5.按照权利要求1、2、3或4所述的预测方法，其特征在于 所述炮孔振动波形预测进程按照以下步骤进行，步骤一，分别读取各所述测振仪测得的各测振点的爆破振动数据，并计算得到所述基 准爆破孔在所述预测目标点的振动波形Mtci)和所述爆破区域的地震波波速Cv ； 步骤二，读取所述基准爆破孔的所述爆破信息；读取各所述炮孔的所述布孔信息； 步骤三，依据所述基准爆破孔的所述爆破信息，计算地震波从所述基准爆破孔传到所 述预测目标点所经过的距离Ltl ；依据各所述炮孔的所述布孔信息，分别计算地震波从各所 述炮孔传播到所述预测目标点将经过的距离Li ； 步骤四，依据公式
6.按照权利要求5所述的预测方法，其特征在于 所述爆破振动衰减指数α的值取为1 3。
7.按照权利要求5所述的预测方法，其特征在于所述地震波波速Cv的值根据地震波传播到任意两个测振点所经过的距离以及所述地震波在该距离传播所对应的传播时间计算。
8.按照权利要求5所述的预测方法，其特征在于所述距离Ltl采用公式L0=^(X0-U0)2 +(y0 -V0)2+(z0-r0f或者公式 L0 = j(x。-U0)2 +(y0-V0)2 +(Z0 -r0 +ψ)2 计算得到；所述距离Li采用公式L1 =^j(Xi-U0)2 +Cyi-V0)2 +(z,-r0)2或者公式L1 = Jixi-u0)2+ (y, -V0)2+ (z,-r0+ ψ)2 计算得到。
9.按照权利要求5所述的预测方法，其特征在于执行所述步骤一读取各所述测振仪测得的所述振动数据时，采用低通滤波的方式读取。
10.按照权利要求9所述的预测方法，其特征在于 所述低通滤波的阈值取为500Hz以内。
11.按照权利要求9所述的预测方法，其特征在于 所述低通滤波的阈值取为300Hz以内。
全文摘要
本发明提供一种爆破振动预测方法。首先确定一个基准爆破孔和多个测振点，并在各测振点处安置测振仪；起爆基准爆破孔以获得不同距离点的爆破振动波形，最后执行爆破振动波形预测进程。上述预测进程中，首先获得实际爆破振动测试数据，再依据预测目标点至各炮孔和基准爆破孔的距离，计算各炮孔在预测目标点的爆破振动波形和偏移叠加延期时间，基于振波叠加原理即可得到预测目标点处的振动波形。采用该方法得到的爆破振动预测结果真实地体现出场地的地形地质条件对爆破振动的影响，并能得到预测目标点的完整振动波形；预测波形的质点振动速度峰值与实测振速峰值很接近、出现的时刻相差也很小，分析结果准确可靠，从而可更好地指导爆破方案的调整。
文档编号G01H17/00GK102135445SQ20101022244
公开日2011年7月27日 申请日期2010年6月30日 优先权日2010年6月30日
发明者孟海利, 张乐, 杨年华, 薛里, 郭尧, 颜景龙 申请人:中国铁道科学研究院铁道建筑研究所, 北京北方邦杰科技发展有限公司