黄河三角洲凹陷区地震勘探振动特征分析.pdf

第37卷第3期 2020年9月 Vol. 37 No. 3 Sep. 2020 bMg do i 10. 3963/j. issn . 1001 -487X. 2020.03.007 黄河三角洲凹陷区地震勘探振动特征分析* 王贤就俺 刘美张伟“，髙載刘志扩,王 （1.青海物产爆破技术服务有限公司，西宁810000 ； ；2,武汉理工大学土木工程与建筑学院，武汉430070； ； 3. 长江科学院水利部岩土力学与工程重点实验室，武汉430010； ； 4. 中石化石油工程地球物理有限公司胜利分公司，东营257000） 摘要地震勘探的振动特性分析是工程安全研究的基础。针对黄河三角洲凹陷区地震勘探震源类型、药 柱埋深和药量的不同，采用现场获得地震勘探振动衰减规律和频率分布特征。研究表明勘探炸药埋深越 浅，震源近区地表餉振动越大，振动速度衰减更快，振动影响范围更小；一般地，竖向的振动速度为径向的 2 5倍,切向的510倍;随药量增大和埋深的增加，萨道夫斯基公式的a和K值均有所减小；机械震源振 动特性则受到基础硬化程度影响。爆破振动主频率的分布在三个方向差异性明显,径向主振频率10 Hz以 下占比82,切向表现为“两头分布多（＞50 Hz和＜20 Hz）、中间（20-50 Hz）分布少”的特点，竖向的高频 部分（＞30 Hz）占比达到68； ；机械震源的主频率在三个方向分布基本一致。 关键词地震勘探；振动衰减；频率分布；机械震源 中图分类号TU45 文献标识码A 文章编号1001 -487X（2020）03 - 0040 - 07 Analysis of Vibration Characteristics of Seismic Exploration in Subsidence Zone of Yellow River Delta WANG Xian-wv） ,RAO Yu23 ,L IU Mei-shan3,ZHANG Wei4, GAO Feng1 ,UU Zhi-fang3, WANG Kai1 （1. Qin gh a i Ma t er ia l s In du st r y Bl a st in g Tech n o l o gy Ser vice Limit ed Co r po r a t io n,Xin in g 810000,Ch in a； 2. Sch o o l o f Civil En gin eer in g a n d Ar ch it ect u r e, Wu h a n Un iver sit y o f Tech n o l o gy, Wu h a n 430070, Ch in a ； 3・ Key La bo r a t o r y o f Geo t ech n ica l Mech a n ics a n d En gin eer in g o f Min ist r y o f Wa t er Reso u r ces, Ya n gt ze River Scien t ific Resea r ch In st it u t e,Wu h a n 430010,Ch in a；4. Sh en gl i Br a n ch,Pet r o l eu m En gin eer in g Geo ph y sica l Limit ed Co r po r a t io n,SINOPEC,Do n gy in g 257000,Ch in a） Abstract Vibration characteristics analysis of seismic exploration blasting is the foundation of eng ineering safety research. Taking the difference of seismic source type, explosive depth and explosive charg e used in the subsidence zone of Yellow River delta into account, the attenuation law of blasting vibration and distribution characteristics of main vibration frequency were obtained throug h field test. The results showed that if the depth of charg e was reduced, the vibration of surface near the seismic source would be increased, the attenuation of vibration velocity would be fas ten ,and the influence rang e of blasting vibration would be smaller. I n g eneral, the vertical component direction had the maximum peak particle velocity（PPV） ,which was 2 5 times of the PPV in radial component direction and 5 10 times of the PPV in tang ential component direction. With the increase of quantity of charg e and exploration depth, the values of K and a in the Sadaovsk ula decreased. The vibration characteristics of mechanical source were af fiected by the deg ree of foundation hardening . The main frequency distribution of blasting vibration was obviously dib ferent in three directions,82 main vibration frequencies were below 10 Hz in radial component direction and 68 were over 30 Hz in vertical component direction. I n the tang ential direction, main vibration frequencies were mostly below 20 Hz and over 50 Hz,and little between 20 Hz and 50 Hz. The main vibration frequency distribution charac- 第37卷第3期王贤武，饶 宇,刘美山，等 黄河三角洲凹陷区地震勘探振动特征分析41 teristics of the mechanical source were almost consistent in three directions. Key Words seismic exploration ； ； vibration attenuation ； ； frequency distribution ； ； mechanical source 振动是地震勘探引起的主要危害之一,随着地 震勘探作业范围向城镇推进，振动所激化的勘探作 业和各类建（构）筑物保护之间的矛盾日益突出，是 影响地震勘探工作顺利开展的重要因素。以胜利油 田东部探区所处的黄河三角洲凹陷区为例，地震勘 探主要采用震源药柱深井爆炸（井炮）和地面机械 震源相结合所进行。为确保勘探资料的质量，检炮 点布置需要满足一定密度且需要尽量降低保护物影 响所带来的不规则分布。因此，地震勘探振动特性 的分析对保护物的安全控制和确保勘探顺利实施具 有重要意义。 爆破或爆炸的振动特性主要包括振动速度和振 动频率两个方面，众多的学者和工程技术人员对此 进行了研究。卢文波通过理论分析研究了主振动频 率的衰减规律和影响因素⑴2〕；胡建华修正了萨道 夫斯基公式并进行了多元回归分析⑶；严鹏研究了 爆破振动的衰减及其与损伤深度的关系⑷；陶铁军 建立了地震波能量的衰减公式⑸；王玉杰获得了地 下厂房的振动衰减规律同；饶宇从滤波角度探讨了 预裂爆破的频率分布特性E ；柴修伟探讨了炸药性 能对爆破振动衰减的影响囲。张加海通过对地震 勘探震源药柱爆炸振动的回归分析获得了位移、速 度和加速度的衰减规律⑼；李云伟根据胶莱地区勘 探振动的衰减规律确定了不同岩体下建筑物的安全 距离3〕。杨风威比较了台阶爆破和预裂爆破的振 动衰减参数和主振频率等特性不同⑴］；周俊汝研究 了主振频率的衰减及其波动特征「⑵。 当前对振动特性的研究多是以微差爆破为基 础,主要针对岩土爆破和拆除爆破等,而地震勘探的 振动特性因其规模较小往往被忽略。实际上，地震 勘探具有震源类型多（爆破震源药柱和机械震源 等）、保护物类别多、涉及到城镇等复杂作业环境、 爆炸频次高、振动能量利用率高（尽可能将爆炸能 量转换为振动能量）、寻求极限的安全距离以满足 勘探精度要求等特征，往往更加容易引起安全问题, 而地震勘探的振动特性研究是安全控制的基础。针 对黄河三角洲凹陷区勘探作业采用不同的震源类 收稿日期 2020-05-28 作者简介王贤武（（1979 -），男,工程师、学士，从事爆破安全相关的 工作， ，（ （E-mail）751034172 qq. como 通讯作者饶 宇（1990-），男,工程师、博士研究生，主要从事岩土 动力学方面的研究工作,（E-mail） raoyuray 163. como 基金项目国家自然科学基金青年科学基金项目（51609017） 型、震源药柱埋深和药量,通过现场试验和原位测试 方法，获得了震源药柱和机械震源的振动衰减规律 和主振频率分布特征,研究了埋深和药量对勘探振 动特征的影响，获得了不同类型震源的振动衰减规 律,为地震勘探的安全控制制定提供了依据。 1地震勘探振动衰减规律 1.1药柱爆炸振动衰减规律药柱爆炸振动衰减规律 本地区的地震勘探常采用的单井药量规格有 1 kg、2 kg、4 kg、6 kg和8 kg,同时，震源一般药柱位 于潜水面以下1 m、3 m、5 m和7 m0对于药柱爆炸 振动的衰减，目前常采用萨道夫斯基公式进行回归 分析。药量不同导致药柱长度、位置不同，从而对衰 减规律产生影响。同样地，震源药柱在潜水面以下 不同位置也导致衰减规律出现差异。因此,需要对 震源位置和药量多少对地震勘探振动速度衰减的影 响进行分析。 式中v为峰值质点振动速度,cm/s；p为比例药 量P Q1/3/R; Q为最大单响药量，kg；丘 心匚宀分为爆心距,m；K为与岩石性质、药柱 参数、场地条件等有关的衰减系数;a为振动衰减 指数。 现场选取了三个不同试验区进行勘探试验分 析,三个试验区的潜水面深度不相同。其中,1地区 不同方向的药柱爆炸振动衰减曲线如图1 图4 所示。 根据图］可知，在勘探药量相同情况下,随着井 深减小,炸药埋深越浅,在距离震源较近的区域（简 称近区，下同，具体范围在不同方向存在差异，一般 震源距在20 m以内）的质点振动速度峰值越大，并 且振动峰值的衰减速度也更快，从而出现“近区大 远区小”的情况，其中，以铅锤向和水平径向尤为明 显。炸药埋深越大，虽然近区振动小，但衰减慢，在 中远区振动速度反而更大。但当炮孔存在冲孔导致 能量外泄时则上述规律是不明显的。 分析图2 图4可知,随药量的增加,振动速度 峰值显著增大，而振动速度峰值在三个方向的数值 大小表现为竖向〉径向〉切向，竖向振速峰值一般 为径向的2 5倍，为切向的5 10倍。竖向为三个 方向中的主控方向。 42爆破2020年9月 S ・ 0 5 0 5 0 5 0 5 0 5 0 5 0 5 0 5 0 5 0 5 5 4 4 3 3 2 2 1 1 5 4 4 3 3 2 2 1 1 0 厂S 吕磐蛾哟赳蠻哑皑 ★水下7m --水下5m △ 水下3m -----水下l m ♦水下7m实测 水下5m实测 *水下3m实测 ■水下l m实测 2 cm/s 线 s UI。、粵赳蚩旺径 ♦水下7 m实测 水下5 m实测 *水下3 m实测 ・水下1 m实测 ・ 5 0 5 5 0 5 3 3 2 3 3 2 5 5 0 5 0 0 5 0 2 1 1 2 1 1 2,cm/s 线 20 40 60 80 1001200 20 40 60 80 100 120 050 100 150 200 震源距/m震源距/m震源距/m a 径向b切向c竖向 a Th e h o r izo n t a l r a dia lb Th e h o r izo n t a l t a n gen t ia lc Th e ver t ica l 图11地区不同震源深度振动衰减拟合曲线8 kg Fig . 1 PPV attenuation curves with different explosive depth in test area one8 kg 5 0 5 0 5 0 5 5 0 5 0 5 0 5 O O 3 3 2 2 1 1 3 3 2 2 1 1 3 3 3 3 3 8 3 3 3 3 3 8 下下下下下合 水水水水水综 ♦水下3 m 1 kg实测 水下3 m2kg实测 *水下3 m 4 kg实测 ■水下3 m6kg实测 v水下3 m 8 kg实测 2 cm/s 线 S U I O M s m 赳缎無底 00005050150 200 40 35 30 25 20 15 10 5 0 水下 5 m 1 kg 水下5 m 2 kg △水下5 m 4 kg --水下5 m 6 kg ♦水下5 m 1 kg实测 水下5m2kg实测 *水下5 m4kg实测 ■水下5 m 6 kg实测 v水下5 m 8 kg实测 20 40 60 80 100 120 震源距/m a 药柱位于水面以下3 m a Th e expl o sive is l o ca t ed t h r ee met er s bel o w t h e wa t er l evel 震源距/m b药柱位于水面以下5 m b Th e expl o sive is l o ca t ed five met er s bel o w t h e wa t er l evel 图2 1地区不同药量竖向振动衰减拟合曲线 Fig . 2 Attenuation curves of vertical PPV with different explosive quantity in test area one 3030 S・日。 、 罄赳營哑蹊 0 5 0 0 5 0 2 1 1 2 1 1 ♦水下3 m 1 kg实测 水下3 m 2 kg实测 *水下3 m 4 kg实测 “水下3 m 6 kg实测 2 cm/s 线 5 0 5 0 5 0 5 0 5 5 0 5 0 5 0 5 0 5 4 4 3 3 2 2 1 1 4 4 3 3 2 2 1 1 下下下下合 水水水水综 - - - ♦水下5 m 1 kg实测 水下5 m 2 kg实测 *水下5 m 4 kg实测 ■水下5 m 6 kg实测 kg kg kg kg kg kg kg kg 2 cm/s 线 杠薜 S・曙赳迤哑抵 0 50 100 150 200 震源距/m a 药柱位于水面以下3 m a Th e expl o sive is l o ca t ed t h r ee met er s bel o w t h e wa t er l evel 0 50 100 150 200 震源距/m b药柱位于水面以下5 m b Th e expl o sive is l o ca t ed five met er s bel o w t h e wa t er l evel 图3 2地区不同药量竖向振动衰减拟合曲线 Fig . 3 Attenuation curves of vertical PPV with different explosive quantity in test area two 将回归结果汇总至表1可知，药量和埋置深度 还影响到振动衰减指数Q和衰减系数K, 一般地，随 着药量增大和埋深的增加，衰减指数和衰减系数 K均有所减小。 12机械振动衰减规律机械振动衰减规律 对于地表可控震源而言,其产生的振动衰减规 律在数学表达式上与药柱爆炸振动类似,只是在药 柱爆炸振动衰减中,采用最大单段药量来表征输入 的振动能量,在地表可控震源中,选择振动出力进行 表征,相关的数学表达如式2所示。 八K罟“硏 ⑵ 式中J为保护对象所在地面质点振动速度, cm/s；p为比例出力,p二W1/n/R； W为振动出力值, 第37卷第3期王贤武，饶 宇,刘美山，等 黄河三角洲凹陷区地震勘探振动特征分析43 kPa；R为平板（重锤）中心至观测点的距离，称震源 距,m；K、a分别为平板（重锤）与保护对象间的地 形、地质条件有关的衰减系数和衰减指数;“为表征 振动出力和震源距指数关系的参数。 s o o o o o o o o o o o o o o o o 8 7 6 5 4 3 2 1 8 7 6 5 4 3 2 1 薜薜 S ♦水下3 m 1 kg实测 水下3 m2 kg实测 ★水下3 m4kg实测 ■水下3 m 6 kg实测 水下3m8kg实测 6060 S 2 cm/s 线 o o o o o o o o o o 5 4 3 2 1 5 4 3 2 1 ri 水下 5 m 1 kg 水下5 m 2 kg a水下5 m 4 kg 水下5 m 6 kg 水下5 m 8 kg 综合8 kg ♦水下5 m 1 kg实测 水下5 m2 kg实测 *水下5 m 4 kg实测 ■水下5 m 6 kg实测 v水下5 m 8 kg实测 2 cm/s 线 0 50 100 150 200 震源距/m a 药柱位于水面以下3 m a Th e expl o sive is l o ca t ed t h r ee met er s bel o w t h e wa t er l evel 0 50 100 150 200 震源距/m b药柱位于水面以下5 m b Th e expl o sive is l o ca t ed five met er s bel o w t h e wa t er l evel 图4 3地区不同药量竖向振动衰减拟合曲线 Fig . 4 Attenuation curves of vertical PPV with different explosive quantity in test area three 可控震源与药柱爆炸不同的是,药柱爆炸很大 部分能量实际上用来破碎挤压岩土而非转化为振 动，而可控震源的能量基本都用于引起岩土体的波 动。为便于采集资料的分析处理，一般同一地区震 源型号、震源参数相同。因此,将本地区的震源参数 （型号为WTC5141TZY车载可控震源、出力设定 70、扫描长度23 s,平均出力67）作为基准值1, 即式（2）中W取值为1。采用震源车分别在硬化地 表、一层砖混民房墙体和土质地基等部位不同距离 的位置进行振动试验，以此研究机械震源的衰减 规律。 分析图5可知，采用当前震源参数的可控震源 引起的质点振动速度峰值远小于1炖震源药柱爆 炸所引起的质点振动速度峰值，振动的衰减速度也 较爆炸振动慢，整体上可控震源引起的振动较小。 可控震源所引起建（构）筑物的振动小于硬化地表 振动，而硬化地表的振动则明显小于天然软弱土层 的振动。 一*一地表第1次试验 -一地表第2次试验 △地表第3次试验 ---民房墙面 。民房土基 地标试验综合 ♦地表1次实测 *地表2次实测 ■地表3次实测 ▼民房墙面实测 v民房土基实测 Qcm/s 线 4 2 0 8 6 4 2 4 2 0 8 6 4 2 1 1 1 1 1 X 1 X S S・ 2 ・ 2 磐 粤赳 哑底 喰地表第1次试验 ---地表第2次试验 A地表第3次试验 ---民房墙面 ，民房土基 *地标试验综合 ♦地表1次实测 *地表2次实测 ■地表3次实测 ▼民房墙面实测 v民房土基实测 0.5 cm/s线 12 10 8 6 4 2 0 5 10 15 20 25 30 震源距/m a 竖向 a Th e ver t ica l 0 5 10 15 20 25 30 震源距/m b径向 b Th e h o r izo n t a l r a dia l .5 .O .5 .O .5O .5 .O .5 .O .5O 2. 2.L L O. 2. 2.L L O. S ・ O ---* 地表第1次试验 --地表第2次试验 . A地表第3次试验 --民房墙面 ＞民房土基 - 地标试验综合 ♦地表1次实测 *地表2次实测 ■地表3次实测 ▼民房墙面实测 ＜民房土基实测 0.2 cm/s 线 5 10 15 20 25 30 震源距/m c切向 c Th e h o r izo n t a l t a n gen t ia l 图5可控震源引起不同对象的振动衰减拟合曲线 Fig . 5 PPV attenuation curves of different objects by vibroseis 可控震源一般以竖向振动最大，在混凝土基础、 已夯实或硬化的地面采用可控震源进行勘探，径向 和切向振动速度均显著小于铅锤向振动速度。 而在软土基础上，因为可控震源可能引起近区 软质土层的挤压变形，从而导致软质土基径向振动 大于竖向，这在震源近区尤为明显,但影响范围较 小,且径向速度衰减比铅锤向速度衰减更快；因此, 软质土层中，径向振动速度呈现“近区大、远区小、 衰减快”的现象。 13衰减参数统计衰减参数统计 根据2.1 2. 2节所述，采用震源药柱进行勘探 作业的振动主控方向为铅锤向，而采用地面可控机 械震源进行勘探作业则视地质条件而定,一般硬化 地表主控方向为铅锤向，而天然地表和软弱基础主 控方向则不同，在近区为径向，中远区为铅锤向。黄 河三角洲凹陷区地震勘探振动衰减参数统计如表1 44爆破2020年9月 所示。 表1黄河三角洲凹陷区地震勘探振动衰减参数统计 Table 1 Attenuation parameters of the Sadaovsk ula by seismic exploration in subsidence zone of Yellow River delta 震源类型 工况 主控方向 比例药量（出力）对数 I n p 回归方程常数 K a 相关系数 地PC 装药高度 水下1 m-3.83 〜-1.78365.281.500.97 水下3 m-4.15 〜-2.02338.841.450.96 S1水下5 m-4.24 〜-2.04350.601.360.96 水下7 m-4.03 〜-2.11169.871.060.99 不同深度累计-4.24--1.78271.801.320.93 水下1 m-3.15 〜-2.69314.601.300.96 水下3 m-4.35 〜-2.15776.111.620.99 震源药柱S2水下5 m铅锤向-3.49--2.25645.401.420.99 水下7 m-3.60--2.79180.601.210.90 不同深度累计-4.35 〜-2.15815.071.590.96 水下1 m-3.82--2.04392.211.150.98 水下3 m-4.28 〜-2.11621.601.450.94 S3水下5 m-4.50--1.94434.021.290.97 水下7 m-3.92--2.26598.661.390.97 不同深度累计-4.50--1.94562.331.390.95 三地区合计-4.50--1.78466.351.400.92 硬化地面铅锤向-3.35 〜-0.643.210.760.96 可控震源软弱土体 铅锤向-3.40--1.612.420.530.86 径向-3.40--1.617.110.860.93 1层砖混民房墙面铅锤向 3.40 〜-1.611.450.380.87 2主振频率的分布 2.1震源药柱的主振频率分布震源药柱的主振频率分布 振动的危害主要与振动速度大小、持续时间、频 率有关。振动速度越大，对建（构）筑物的危害越 大;主频率越接近于建（构）筑物的自振频率，易形 成共振现象,对建筑物的危害大。振动持续时间越 长,显然危害越大，但对于地震勘探,振动持续时间 基本固定。根据现场实测地震勘探振动监测所得的 峰值振速频率数据绘制频率分布图如图6。 a 径向 a Th e h o r izo n t a l r a dia l b切向 b Th e h o r izo n t a l t a n gen t ia l c竖向 c Th e ver t ica l 图6药柱爆炸振动不同方向主振频率分布 Fig . 6 The distribution of main vibration frequency by explosive seismic blasting in different direction 根据图6可知，井炮的主振频率在三个方向上 20皿）分布多冲间（20〜50皿）分布少”的特点， 存在较大的差异，径向的主振频率在10 Hz以下占 竖向与切向和径向又显著不同，高频部分 到82,而切向则表现为“两头（50 Hz和v （ 30 Hz）占比达到68 o 第37卷第3期王贤武，饶 宇,刘美山，等 黄河三角洲凹陷区地震勘探振动特征分析45 考虑到竖向为三个方向中的主控方向,竖向主 振频率高对保护物的安全性有利。但地震勘探爆破 的主振频率分布离散、特别是径向低频、超低频占比 高。仅以主振动及其频率对建筑进行安全评价存在 不足，未考虑“低振速且低频率”的情况。 2.2机械震源的主振频率分布机械震源的主振频率分布 图7为地面可控震源所引起三个方向的主振频 率区间分布图。 a 径向 a Th e h o r izo n t a l r a dia l b切向 b Th e h o r izo n t a l t a n gen t ia l c竖向 c Th e ver t ica l 图7可控震源不同方向主振频率分布 Fig . 7 The distribution of main vibration frequency by vibroseis in different direction 根据图7可知,可控震源的主振频率分布在三 个方向上相对固定，主振频率在10〜30 Hz之间占 到90以上，竖向主振频率在20 -30 Hz之间比例 要高于其他两个方向。机械震源主振频率高对建筑 物保护总体有利。 3结论 通过对黄河三角洲凹陷区现场实测的地震勘探 振动的衰减规律和频率分布特征进行了分析，获得 如下结论 1 地震勘探药柱埋深越浅，近区的地表振动 越大，而振动速度衰减也更快;药柱埋深越大，中远 区的振动速度更大，振动影响范围更广。 2 勘探爆破竖向质点振动速度峰值最大，一般 为径向的2〜5倍，为切向的5〜10倍;随药量增大和 埋深增加，衰减指数和衰减系数K均有所减小。 3 可控震源振动特性受地表硬化程度影响， 地表越软则近区引起的振动更大，但衰减速度也更 快;一般竖向振动速度最大，但在软土上，径向可产 生更大的振动。 4 震源药柱爆炸的径向主振频率10 Hz以下 占比82,切向主振频率基本分布在50 Hz以上和 20 Hz以下，竖向的主振频率在30 Hz以上占比达到 68。机械震源的主振频率分布在三个方向上较为 一致,主振频率在10 30 Hz之间占到90以上。 参考文献参考文献References [1]卢文波，张 乐,周俊汝，等爆破振动频率衰减机制 和衰减规律的理论分析[J] 爆破,2013,3021-6, 11. [1 ] LU Wen-bo, ZHANG Le, ZHOU Jun-ru, et al. Theoretical analysis on decay mechanism and law of blasting vibration frequency[ J]. Blasting,2013,30 2 1-6,11. in Chi nese [2]卢文波，周俊汝，陈 明，等.爆破振动主频衰减公式 研究[J].工程爆破,2015,2161-6,24. [2] LU Wen-bo, ZHOU Jun-ru, CHEN Ming, et al. Study on attenuation ula of dominant frequency of blasting vi bration [J]. Eng ineering Blasting,2015,216 1-6,24. in Chinese [3] 胡建华，尚俊龙,罗先伟，等.单孔爆破振动监测与衰 减规律多元线性化回归[J].振动与冲击，，2013, 3216 49-53. [3] HU Jian-hua,SHANG Jun-long ,LUO Xian-wei,et al. Mo nitoring of sing le-hole blasting vibration and detection of its attenuation law by using multiple linear reg ression[ J]. Journal of Vibration and Shock,2013 ,3216 49-53. in Chinese [4] 严 鹏，邹玉君,卢文波，等.基于爆破振动监测的岩 石边坡开挖损伤区预测[J]岩石力学与工程学报， 2016,353 538-548. [4] YAN Peng,ZOU Yu-jun,LU Wen-bo,et al. Predicting the damag e zone of rock slopes under blasting excavation based on vibration monitoring [ J ]. Chinese Journal of Rock Mechanics and Eng ineering ,2016,353 538-548. in Chinese [5] 陶铁军，汪旭光,池恩安，等.基于能量的爆破地震波 衰减公式[J].工程爆破,2015,21678-83. [5 ] TAO Tie-jun, WANG Xu-g uang, CHI En-an, et al. Attenu ation ula of blasting seismic wave based on energ y theory [J]. Eng ineering Blasting ,2015 ,216 78-83. in Chinese [6] 王玉杰，梁开水，田新邦.周宁水电站地下厂房开挖爆 破地震波衰减规律的研究[J].岩石力学与工程学报， 200522 4111-4114, 46爆破2020年9月 [6] WANG Yu-jie, LI ANG Kai-shui, TI AN Xin-bang . Study on redundant reg ulation of underg round dig g ing blasting vibration of zhouning hydropower station [ J ]・ ・ Chinese Journal of Rock Mechanics and Eng ineering, 2005 22 4111-4114. in Chinese [7] 饶 宇,夏元友湖英国，等预裂缝对爆破振动频谱分布 特征的影响[J] 振动与冲击,2017,367 191-198. [7] RAO Yu,XI A Yuan-you,HU Ying -g uo,et al. I nfluence of pre-splitting crack on spectrum distribution characteristics of blasting vibration [ J ] . Journal of Vibration and Shock, 2017,36刀 191-198. in Chinese [8] 柴修伟不同炸药性能对爆破振动效应的影响[J].爆 炸与冲击,2011,315548-552. [8] CHAI Xiu-wei. Effects of different explosive properties on blast vibration [ J ]. Explosion and Shock Waves, 2011, 315 548-552. in Chinese [9] 张加海，王敬绵,张光德，等.炸药震源激发参数与场 地质点振动响应关系探讨[J].油气地质与采收率， 2010,172 57-59,63,115. [9 ] ZHANG Jia-hai, WANG Jing -jin, ZHANG Guang -de, et al. Relationship between seismic source parameters and particle response [ J ]. Petroleum Geolog y and Recovery Efficiency,2010,172 57-59,63,115. in Chinese [10]李云伟.胶莱地区地震勘探安全距离研究探讨[J]・ ・ 工程地球物理学报,2012,95 587-594. [10] LI Yun-wei. The study of seismic exploration safe dis tance in Jiaol