• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

地下核爆诱发工程性地震实测数据分析与不可逆变形范围计算

徐天涵 李杰 王明洋 徐小辉 何雯静

徐天涵, 李杰, 王明洋, 徐小辉, 何雯静. 地下核爆诱发工程性地震实测数据分析与不可逆变形范围计算[J]. 爆炸与冲击, 2019, 39(12): 121101. doi: 10.11883/bzycj-2018-0505
引用本文: 徐天涵, 李杰, 王明洋, 徐小辉, 何雯静. 地下核爆诱发工程性地震实测数据分析与不可逆变形范围计算[J]. 爆炸与冲击, 2019, 39(12): 121101. doi: 10.11883/bzycj-2018-0505
XU Tianhan, LI Jie, WANG Mingyang, XU Xiaohui, HE Wenjing. Analysis of test data of underground nuclear explosions and calculation of irreversible deformation range[J]. Explosion And Shock Waves, 2019, 39(12): 121101. doi: 10.11883/bzycj-2018-0505
Citation: XU Tianhan, LI Jie, WANG Mingyang, XU Xiaohui, HE Wenjing. Analysis of test data of underground nuclear explosions and calculation of irreversible deformation range[J]. Explosion And Shock Waves, 2019, 39(12): 121101. doi: 10.11883/bzycj-2018-0505

地下核爆诱发工程性地震实测数据分析与不可逆变形范围计算

doi: 10.11883/bzycj-2018-0505
详细信息
    作者简介:

    徐天涵(1995- ),男,博士研究生,martinxu41@126.com

    通讯作者:

    李 杰(1981- ),男,博士,副教授,lijierf@163.com

  • 中图分类号: O383.1

Analysis of test data of underground nuclear explosions and calculation of irreversible deformation range

  • 摘要: 地下核试验瞬间释放的巨大能量引起地壳能量的连锁反应,产生诱发地震等地球物理现象。本文对20世纪前苏联与美国进行的地下核试验数据进行整理与归纳,给出地下核爆炸试验诱发工程地震的范围以及激活岩块大小等。根据实测数据,指出地下核爆炸诱发工程性地震的力学本质,利用理论公式计算了地下核爆炸产生不可逆位移的临界能量因子范围,为相关研究提供理论基础和场地效应试验数据。
  • 地壳岩体是具有复杂内部构造层次的介质,在长期的重力和地质构造应力的作用下,不断聚集能量,当积聚的变形超出了岩体的承受能力,就会产生突然的断裂或沿已有断裂面的错动,在构造运动中长期积累的能量得到迅速释放,从而形成地震(构造地震)[1]。研究表明,在地壳中的弹性应变能被释放掉之前,地壳中的原有断层都具有潜在活化的可能,任何影响断层稳定性的微小扰动,例如水库蓄水、采矿和地下核爆炸等,都存在诱发地震的可能性[2]。理解地震诱发的原因及其破坏性影响,是地震学和防护学研究的重点。

    早在20世纪50~60年代,美国、苏联的科学家就注意到,在地下核爆炸几天之后,有时会在几百甚至几千米外发生地震,释放的总能量超出爆炸本身千倍量级[3]。这一巨大的摧毁能力,在冷战的格局下,甚至被作为地震武器的设想加以关注。美国、俄罗斯等国家高度重视地震触发机制的研究,例如Adushkin等[4-6]、Rodinov等[7]、Kocharyan[8-11]、Nikolaev等[12]都对于地震的触发机制进行了广泛的实验与理论研究,并每年召开地球系统的地震触发机制学术会议。

    目前关于核爆炸触发地震的可能性已基本取得学界共识,但相关研究仍处于定性解释的阶段,远未达到机理清晰、模型准确的程度。例如:目前已知,核爆炸是否会诱发地震,以及诱发地震的大小取决于场地的地质背景。特殊的岩性组合、有利于地震发生的断层与不连续结构面的存在以及岩体应力状态等,都是诱发地震的必要条件。然而其充分条件是什么?在一定的地质条件和爆炸当量条件下,会否诱发地震?诱发地震的量级是多少?能否对修建于岩层中的地下工程产生影响?可否采取相应的工程措施以减少即将到来的工程灾害?······这些问题目前尚未有一致性的相关研究成果。研究的困难主要在于工程地质条件的复杂不可探测性以及试验样本数量的稀少。

    目前国际社会已全面禁止核试验,对核爆炸诱发地震的研究只能立足已有的核试验数据。据不完全统计,自1945年首次核试验至今,全球共进行2 000余次核试验,其中地下核试验1 500余次。但由于试验国家对数据的保密以及受当时落后测量条件的制约,公开的试验数据及其完整性十分有限,且这些数据分散于各种文献,十分不利于研究。

    基于上述理由,本文对各文献记录的地下核爆炸试验场地数据进行整理汇总,并归纳总结地下核爆炸试验诱发工程地震的范围、临界能量阈值、激活岩块大小等,为相关研究提供场地效应试验数据。

    1966年12月20日,美国在内华达试验场Pahute Mesa试验区深1 215 m处进行了GREELY核爆炸[13]。爆炸当量为825 kt,周边岩体介质为沸石凝灰岩。爆炸震级为6.1级,其最远地表错动距爆心投影点(爆心在地表的垂直投影)约5 200 m,超出爆心投影点约760 m的错动基本发生在先前存在的自然断层附近。距爆心投影点若干千米外的地表最大相对垂直位移为120 cm。爆炸产生最大裂缝长度为6.9 km,地表断层及位移分布见图1(图中数字单位为cm,下同)。

    图  1  GREELEY引起的主要地表错动
    Figure  1.  Ground motion induced by GREELEY

    GREELY核爆炸主要引起3条断层的水平相对位移。第1条断层F-1在距爆心投影点东侧1 200 m之外,在点A南部出现长约1 800 m的连续右旋水平相对位移,最大为15 cm,而点A北部3 500 m范围内只有垂直相对位移出现。其余两处为断层F-2与F-3,分别在爆心投影点北部1 700 m和3 000 m处。在断层F-2上,若干厘米的右旋水平相对位移以及垂直位移从爆心投影点北部3 700 m处一直向南延伸直至F-2在距爆心投影点北部约1 700 m处分叉成两条断层。分叉点南部的两条断层都出现了从几厘米至最大120 cm的垂直相对位移。断层F-3处仅出现若干厘米的右旋水平位移。

    1966年4月14日,美国在内华达试验场Pahute Mesa区深547 m处进行了DURYEA核爆炸[13]。爆炸当量为65 kt,介质为流纹岩。爆炸震级为5.1级,产生的地表错动最远距爆心投影点790 m,爆炸产生最大裂缝长度为1.1 km,地表断层及位移分布见图2

    图  2  DURYEA引起主要地表错动
    Figure  2.  Ground motion induced by DURYEA

    几乎所有距爆心投影点600 m以内的自然断层都发生错动。断层F-2与F-3上最大垂直相对位移为12 cm,F-2北部端点附近出现大约5 cm的左旋水平相对位移。最显著的运动发生在断层F-1上,虽然在该断层上可见的垂直相对位移为10~12 cm,但在一条东西走向穿过断层的等高线表现出的位移为18 cm。更重要的是,沿该等高线的位移随距断层距离而衰减的梯度比穿过的其他断层要小。位移平缓的变化意味着断层F-1的错动延伸至一定深度。

    1968年4月26日,美国于内华达试验场Pahute Mesa区深1 158 m处进行了BOXCAR核爆炸[13]。爆炸当量为1 200 kt,介质为沸石凝灰岩。爆炸产生的沿已有断层的破裂最远距爆心投影点6 100 m,激活最大裂缝长度约为8.0 km,爆炸震级为6.2级。

    爆炸引起最明显的错动在距爆心投影点东部915 m处的南北走向断层上(图3)。错动并不与断层位置完全一致,但其关系是明显的。沿该断层的错动长达7 930 m,最大垂直相对位移为1 m,位移距爆心投影点东南方1 070 m处。一般来说地表位移随距爆心投影点距离增加而减小,但在该爆心投影点南部2 590 m处出现30 cm垂直相对位移,而在北部2 750 m处出现45 cm垂直相对位移。在爆心投影点东部观察到大于10 cm的左旋水平相对位移,爆心投影点东北方2 290 m处观察到约6 cm的左旋水平相对位移。

    图  3  BOXCAR引起主要地表错动
    Figure  3.  Ground motion induced by BOXCAR

    Cloud等[14]还对地下爆炸引起的地表运动进行了测量。加速度测量结果见图4,横坐标为测点距爆心投影点距离D,纵坐标为加速度a与重力加速度g的比值。图5所示为爆炸引起地表振动的幅值,横坐标为测点距爆心投影点距离D,纵坐标为地表振动最大振幅A,测点分别位于岩石与冲积层之上。图6所示为地表振动振幅A与振动周期T的比值,该比值可以表征速度的大小。

    图  4  BOXCAR引起地表最大加速度
    Figure  4.  Maximum ground acceleration induced by BOXCAR
    图  5  BOXCAR引起地表运动幅值
    Figure  5.  Amplitude of ground motion induced by BOXCAR
    图  6  BOXCAR引起地表振动最大速度
    Figure  6.  Maximum ground velocity induced by BOXCAR

    1968年1月19日,美国在内华达Hot Creek山谷深975 m处进行FAULTLESS核爆炸[13]。观测到的爆炸产生最远的断裂带位于爆心投影点西侧2 750 m处,观测到最大垂直位移为3 cm。爆炸产生最大裂缝长度为4.1 km,爆炸震级为6.1级。

    爆炸产生的最显著的效应是在距爆心投影点3 200 m范围内的沿线性构造的错动。3处主要的地表错动分别沿线性构造A、B和C(图7)。其他错动几乎都不同程度地与剩下的线性构造相一致。

    图  7  FAULTLESS引起地表错动
    Figure  7.  Ground motion induced by FAULTLESS

    为方便描述,将地表错动以距爆心投影点600 m为界分为两组。

    1.4.1   600 m以内的地表错动

    在距爆心投影点600 m内的两条向东北延伸的主要断裂带形成了地堑的边界。在地堑内有另一条断裂带与其余两条平行,且爆心投影点位于其西北方约60 m处。地堑西北方的边界上最大垂直相对位移为2 m,到端点处降为0。该边界上也出现了30 cm的右旋水平相对位移。地堑东南边界的断裂带最大垂直位移为4.5 m,位于爆心投影点南部。从最大位移处向东北方向垂直位移逐渐且无规律地衰减,向西方向衰减较快。在该断裂带上,位于爆心投影点南部处存在1 m的右旋水平相对位移,该处东北方约760 m位置出现约60 cm的左旋水平相对位移。断裂带端点附近位移降至3 cm以内。在爆心投影点南部约60 m处的断裂带上有最大3 m的垂直相对位移,沿该断裂带到处存在左旋水平相对位移。

    1.4.2   600 m以外的地表错动

    距爆心投影点600 m外的主要断裂带包括:(1)地堑两侧边界的延伸部分;(2)爆心投影点东南部60 m处断裂带的延伸部分;(3)与线性构造A相关联的很长的断裂带系统;(4)在爆心投影点西南方约1 830 m处爆炸前存在的线性构造上的断裂带;(5)一处山前断裂带;(6)爆心投影点东南方向一条2 440 m长的断裂带。地堑西北边界和东南边界上的断裂带总长分别为1 980 m和2 744 m。地堑中间的断裂带总长为1 525 m。爆心投影点西南方1 830 m处的断裂带没有明显的垂直运动。爆心投影点600 m外最明显的断裂为线性构造A处断裂带系统。该线性构造为全新世时期由断层形成的位于冲积层上的断崖。内华达试验场中在该种类型断层附近的地下核爆炸都会引起其运动。在此次爆炸中垂直相对位移在端点处几乎不可见,在中间达到最大值60 cm。右旋水平位移在个别位置出现。大多数的位移出现在一条主要断裂带上,其余位移很少或没有的断裂带与主要断裂带平行或斜交,并形成宽300 m的区域。两条与A斜交的断裂带上最大垂直相对位移为30 cm。

    1968年12月19日,美国在内华达深1 402 m处进行了BENHAM爆炸[15]。岩体介质为凝灰岩,当量1 100 kt。最远地表位移出现在距爆心投影点11.6 km处。爆炸所引起最长裂缝为9.0 km,爆炸震级为6.1级。

    爆炸后的6个星期内,在距爆心投影点12.9 km范围出现10 000余次微弱余震,其中640处绘于图8内用并小叉表示。图中细线表示已知断层,粗线表示出现位移的断层。在爆心投影点附近地表位移大都出现在已有断层附近,其中最大垂直相对位移为0.46 m,出现在爆心投影点北部2.41 km与4.02 km处。其他地区会出现更大的位移,但在断层附近的垂直相对位移大都小于0.3 m。这些位移虽不连续,但能绵延长达数公里。

    图  8  BENHAM引发其中640处余震及断层位移分布
    Figure  8.  Location of 640 aftershocks and faults induced by BENHAM

    爆炸引起的地表加速度见图9[14],地表振动最大幅值与振动周期之比见图10,其坐标轴意义分别与图4图6相同。

    图  9  BENHAM引起地表最大加速度
    Figure  9.  Maximum ground acceleration induced by BENHAM
    图  10  BENHAM引起地表振动最大速度
    Figure  10.  Maximum ground velocity induced by BENHAM

    1969年10月2日,美国在位于阿留申群岛的Amchitka岛深1 219.2 m处进行了MILROW核爆炸,爆炸当量为1 000 kt,介质为枕熔岩。爆炸激活断层长度为0.3~8.0 km,沿断层垂直相对位移为1.0~1.2 m,水平相对位移为0.15 m[5]

    该群岛位于地球上地震最频繁的区域之一[15]。爆炸形成的空腔及其上方“烟囱体”在37 h后完全塌陷。12个地质构造活动被记录[16],其中10个发生在爆炸后41 d内,所有活动的震级均小于3级。图11为爆心经纬度坐标及附近地形,爆炸引发地质活动在图中用圆圈表示。

    图  11  MILROW与CANNIKIN位置及附近地形
    Figure  11.  Location of MILROW and CANNIKIN and nearby terrain

    1971年11月6日在Amchitka岛下方1 790.7 m处进行了CANNIKIN核爆炸,当量“小于5 000 kt”[15]。爆炸引起的空腔及其上方“烟囱体”在38 h内完全塌陷。爆炸后记录到22个构造活动,除了其中一个,其余所有都在爆炸后23 d内发生,所有活动震级均小于4级。

    朝鲜于2017年9月3日进行了NKNT 6[17]地下核爆炸,地点为朝鲜北部丰溪里试验场万塔山下。根据美国地质调查局与中国地震台网监测显示,此次爆炸震级约为6.3级,远大于朝鲜自2006年以来的前5次地下核试验。

    南洋理工大学利用德国航天局TerraSAR-X卫星监测了核爆炸前后地表的变化(图12)。监测显示,核爆炸导致万塔山较陡峭的西侧和南侧向外位移3~4 m,产生的不可逆位移范围约为9 km2。根据地表位移推算得此次核爆炸当量约为120~304 kt,爆源位置为129.078°E/41.300°N±50 m,海拔(1750±100) m,即万塔山顶下方(450±100) m处。

    图  12  朝鲜核爆前后地表变形场
    Figure  12.  Deformation field before and after NKNT 6

    除了上述被详细介绍的地下核爆炸试验,尚有许多未被详细描述的爆炸试验,其数据散见于各文献,集中汇总如下。

    表1[13]为在内华达试验场其他几次试验中爆炸规模与最大激活断层长度的统计。

    表  1  爆炸规模与最大断层长度
    Table  1.  Magnitude of explosions and maximum length of active faults
    爆炸名称当量/Mt震级/mb最大激活断层长度/km
    SCOTCH0.155.41.54
    CHARTREUSE0.075.11.66
    HALFBEAK0.305.93.06
    STINGER5.32.61
    LONGSHOT6.12.25
    下载: 导出CSV 
    | 显示表格

    此外,另有表4~6以及图14见下文。

    从上述地下核爆现场数据可以看出,地下核爆炸在距离爆心很远(r/Q1/3=(6001200)m/kt1/3r为距爆心距离,Q为爆炸当量)的距离上诱发了工程性地震。按地下核爆炸地冲击波的试验监测结果,地下核爆炸时,应力波强度随距爆心距离的增加按rn速率衰减[11]

    σ=ρcpv,v=A(r/Q1/3)n
    (1)

    式中:σ为波前的应力,v为质点速度,ρ为介质密度,cp为介质纵波波速,An为系数,可由试验检测得到。Radionov等[7]基于大量地下核爆炸试验,给出不同岩石中An统计值见表2

    表  2  不同岩石的A、n取值
    Table  2.  Values of A and n of different rocks
    岩石类型An
    花岗石(1.0~1.3)×1041.6~1.75
    盐岩(0.8~1.0)×1041.6
    凝灰岩(0.3~0.4×10)41.6
    下载: 导出CSV 
    | 显示表格

    对于100万吨量级地下核爆炸,当介质为花岗岩时,其在r/Q1/3=(6001200)m/kt1/3距离上的应力为1~3 MPa,远小于岩石的破裂强度。

    在如此小的扰动应力作用下,岩体介质本身不会产生破坏,介质中出现的不可逆现象与现实岩体的块状构造体系有关,这种破坏表现为软弱夹层两侧块体的相对不可逆位移与转动。在自然状态下,岩体中的地应力可分为自重应力与构造应力,其中自重应力是上部岩体在重力作用下对下部岩体产生的应力,而构造应力主要是由于板块之间的相互挤压,以及随之产生的地形、地势等因素引起的应力[18]。即使在深部岩体中构造应力也可以大于自重应力。在地应力作用下,深部岩体块体之间的软弱面存在内应力,但在摩擦力及粘结力作用下尚未发生位移,此时岩体介质存储有应变能。当有应力扰动经过时,会产生两种效应:为软弱面提供了附加剪应力;出现超低摩擦效应导致接触面摩擦力减小[19]。正是由于这两种效应的作用,使原本不会发生位移的软弱面发生了不可逆变形,同时介质应变能释放,若释放的应变能很大,则会形成工程性地震。块系构造介质的变形主要由块体的变形及块体之间软弱夹层的变形组成,因此在爆炸远区,扰动应力很小时岩体的变形主要集中于块体边界及强度软弱区,变形的大小与初始内应力的大小有关,初始内应力越大,则扰动释放的能量越多,产生更大的不可逆位移。根据岩体的具体块体系列及块体之间夹层强度大小,局部的破坏可能在103 m/kt3的距离上发生。已知的块系尺度规律表明,自然岩体中块系构造是普遍存在的,若某一岩体不包含地质断层,则其在更大的尺度上一定具有强度软弱带,并且不可逆变形的大小与发生破坏的块系尺度水平是对应的。

    在确定由地下核爆炸引起的岩体介质变形区范围时,必须考虑103 m/kt3距离上可能会出现的局部不可逆变形。

    不可逆变形区半径rd(以爆心为球心)的确定过程如下:把岩体理解成由若干线性尺寸为L的构造单元组成的整体。构造单元即为组成块系的单位块体,其可以在很大程度进行自由的空间运动。如果爆炸作用引起的变形不超过临界值ε,就不会产生不可逆变形。在块系介质中最大变形发生在软弱面附近,因此极限变形的尺寸可取ε=u0/L,这里u0为邻近块体面之间的相对位移。不可逆变形区可由下式确定[20]

    rdQ1/3(2ABπεcpL)1nmQ13(nm)
    (2)

    式中:Bm同样为场地试验确定的参数。

    为了获得具体数值,必须给出ε的数值。可以把现实的岩体看成是存在时间很长的非平稳系统,作为ε的最小估计值,可以取地壳在自然变形过程中达到的数值。根据潮汐运动及其他地质构造运动中积累的众多资料,ε可以取为(1~2)×10−5[11]

    也可以借助于介质中不可逆位移发生时的质点速度来确定不可逆位移的范围。王明洋等[21]分析了介质中质点速度与能量的关系,首次提出用无量纲能量因子确定岩块不可逆位移的边界范围:

    rdQ1/3=φkd12n,φ=(Acp)1n
    (3)
    kd=vcp=WMc2p
    (4)
    kd=ε2
    (5)

    式中:kd为特征能量因子,W为爆炸扰动能量,M为岩石块体质量。

    将地下核爆炸试验测得的地表质点运动速度代入式(4),求得特征能量因子kd后代入式(3)即可确定不可逆变形区半径,结合爆炸埋深利用勾股定理可求得地表的不可逆变形范围。

    kd的物理意义是在爆炸作用下单位体积岩石的能流密度,在不可逆变形边界处的kd即为岩体发生不可逆位移的阈值。在实际爆炸中,利用地表不可逆位移范围及爆炸埋深可求得rd,代入式(3)即可求得作为阈值的特征能量因子kd

    表3中的计算结果可得块体产生位移的临界能量因子范围为10−9~10−10,该结果与Kurlenya等[22-23]基于摆型波理论给出的范围一致,也与根据式(5)得到的结果一致。Adushkin等[4]认为,当kd≥ (1~4) × 10−9时,整个岩块进入准共振状态。

    表  3  不可逆位移能量因子计算
    Table  3.  Calculation of irreversible displacement
    核爆炸试验岩性当量/kt埋深/m地表不可逆范围/m不可逆位移半径/(m·kt−1/3)vkd
    Greely沸石凝灰岩 8251 215 5 200 5700.164.9×10−10
    Duryea流纹岩 65 547 790 2390.638.0×10−9
    Boxcar沸石凝灰岩1 2001 160 6 100 5840.154.6×10−10
    Benham凝灰岩1 1001 40013 0001 2600.041.6×10−10
    Milrow枕熔岩1 0001 219 8 000 8090.091.6×10−10
    下载: 导出CSV 
    | 显示表格

    根据前苏联场地试验数据[5],不可逆位移范围为650~1 400 m/kt1/3,将表3中的kd代入式(3),取花岗岩的参数,可得不可逆位移范围为658~1 270 m/kt1/3,与试验结果一致。

    爆炸引起不可逆位移与激活块体尺寸是相对应的,且根据式(2),不可逆位移范围也与之有关。为研究在距爆心不同距离上激活块体尺寸的规律,前苏联在长的地下坑道中进行了相应观测。

    表4[10]为前苏联地下爆炸试验中在坑道中测得的块体尺寸。

    表  4  爆炸激活块体尺寸
    Table  4.  Dimension of block that became more active due to explosive action during explosions
    爆炸当量/kt距爆心投影点距离/m激活块体尺寸/m爆炸当量/kt距爆心投影点距离/m激活块体尺寸/m
    2.3 60111.4 150 12
    12.511035 165 15
    1.41202517 430 80
    1110012 500 62
    3 9015 560 55
    10012 625 78
    11012 680 40
    1201178 800100
    13018 900100
    140151 020155
    1501054 750 80
    15510 800 70
    19010 900115
    1.4110161 000100
    125 81 100 90
    13514
    下载: 导出CSV 
    | 显示表格

    从试验资料可以看出,随着距爆心距离增大,激活块体尺寸有增大的趋势,即在远处有可能诱发大的断层位移。当大尺寸块体被激活时,其内部的小尺寸块体同样被激活。确定激活块体尺寸可以帮助在地下工程选址时避开有可能被激活的不连续构造。

    图13[5]为前苏联某测试场布置图。a为Kalba-Chingiz断层,I、II分别为3#、4#测量区,布置有测量仪器。编号为1至4的圆孔为核爆炸孔,试验编号分别为1352、1350、1348及1346。3#及4#测得数据见表5表6,其中D/Q1/3为地表距爆心投影点距离与爆炸当量开立方之比,S/Q1/3为块体位移与爆炸当量开立方之比。

    图  13  测试场布置图
    Figure  13.  Test site layout
    表  5  3#测量区数据
    Table  5.  Results of light distance measurements at Site 3#
    试验试验日期距离D /m(D·Q−1/3)/(m·kt−1/3)位移S/mm(S·Q−1/3)/(mm·kt−1/3)
    13481987.08.021 70039819.54.56
    13501988.09.143 500674122.31
    13461988.12.173 900890143.19
    13521989.07.085 7001 74382.45
    14101989.09.02~7 0004 00021.07
    下载: 导出CSV 
    | 显示表格
    表  6  4#测量区数据
    Table  6.  Results of light distance measurements at Site 4#
    试验编号试验日期距离D /m(D·Q−1/3)/(m·kt−1/3)位移S/mm(S·Q−1/3)/(mm·kt−1/3)
    13481987.08.023 7008666.661.56
    13501988.09.141 45027919.73.8
    13461988.12.176 6501 51851.14
    13521989.07.082 9008873.51.07
    下载: 导出CSV 
    | 显示表格

    图14[5]为地下核爆炸试验中标准化距离与激活块体尺寸之间的关系,数据来自前苏联与美国的地下核爆炸试验。横坐标为爆心投影点距离与爆炸当量开立方之比,纵坐标为激活块体尺寸与爆炸当量开立方之比。

    图  14  标准化距离与块体尺寸关系
    Figure  14.  Scaled size of blocks activated by a nuclear explosion

    激活块体尺寸可用下式计算[24]

    LQ1/3=A(DQ1/3)mn,A=2ABk2πcpςk1(1+k2)
    (6)

    式中:ς=102101为间隙填充物质的极限压缩系数,k1k2为Sadovsky块系构造理论中的岩体不变量。利用式(6)对图14中的数据进行拟合得到A0.105,mn1.04

    (1)在扰动应力远低于岩石介质强度的爆炸远区,岩体会发生工程性地震导致的不可逆位移,这种位移是由于自然岩体块系构造中的断层及块体间隙两侧岩体积累的内应力在扰动作用下释放而导致的,主要集中于块体边界及强度软弱区。

    (2)根据试验数据可以确定不可逆位移半径rd以及发生不可逆位移的特征能量因子kd的阈值。当特征能量因子大于阈值时将发生不可逆位移,计算得阈值为10−9~10−10

    (3)任何岩块都在某一尺度的块系构造内,应力扰动在不同距离上激活的块体尺寸不同,随距离增大,激活块体尺寸增大。

  • 图  1  GREELEY引起的主要地表错动

    Figure  1.  Ground motion induced by GREELEY

    图  2  DURYEA引起主要地表错动

    Figure  2.  Ground motion induced by DURYEA

    图  3  BOXCAR引起主要地表错动

    Figure  3.  Ground motion induced by BOXCAR

    图  4  BOXCAR引起地表最大加速度

    Figure  4.  Maximum ground acceleration induced by BOXCAR

    图  5  BOXCAR引起地表运动幅值

    Figure  5.  Amplitude of ground motion induced by BOXCAR

    图  6  BOXCAR引起地表振动最大速度

    Figure  6.  Maximum ground velocity induced by BOXCAR

    图  7  FAULTLESS引起地表错动

    Figure  7.  Ground motion induced by FAULTLESS

    图  8  BENHAM引发其中640处余震及断层位移分布

    Figure  8.  Location of 640 aftershocks and faults induced by BENHAM

    图  9  BENHAM引起地表最大加速度

    Figure  9.  Maximum ground acceleration induced by BENHAM

    图  10  BENHAM引起地表振动最大速度

    Figure  10.  Maximum ground velocity induced by BENHAM

    图  11  MILROW与CANNIKIN位置及附近地形

    Figure  11.  Location of MILROW and CANNIKIN and nearby terrain

    图  12  朝鲜核爆前后地表变形场

    Figure  12.  Deformation field before and after NKNT 6

    图  13  测试场布置图

    Figure  13.  Test site layout

    图  14  标准化距离与块体尺寸关系

    Figure  14.  Scaled size of blocks activated by a nuclear explosion

    表  1  爆炸规模与最大断层长度

    Table  1.   Magnitude of explosions and maximum length of active faults

    爆炸名称当量/Mt震级/mb最大激活断层长度/km
    SCOTCH0.155.41.54
    CHARTREUSE0.075.11.66
    HALFBEAK0.305.93.06
    STINGER5.32.61
    LONGSHOT6.12.25
    下载: 导出CSV

    表  2  不同岩石的A、n取值

    Table  2.   Values of A and n of different rocks

    岩石类型An
    花岗石(1.0~1.3)×1041.6~1.75
    盐岩(0.8~1.0)×1041.6
    凝灰岩(0.3~0.4×10)41.6
    下载: 导出CSV

    表  3  不可逆位移能量因子计算

    Table  3.   Calculation of irreversible displacement

    核爆炸试验岩性当量/kt埋深/m地表不可逆范围/m不可逆位移半径/(m·kt−1/3)vkd
    Greely沸石凝灰岩 8251 215 5 200 5700.164.9×10−10
    Duryea流纹岩 65 547 790 2390.638.0×10−9
    Boxcar沸石凝灰岩1 2001 160 6 100 5840.154.6×10−10
    Benham凝灰岩1 1001 40013 0001 2600.041.6×10−10
    Milrow枕熔岩1 0001 219 8 000 8090.091.6×10−10
    下载: 导出CSV

    表  4  爆炸激活块体尺寸

    Table  4.   Dimension of block that became more active due to explosive action during explosions

    爆炸当量/kt距爆心投影点距离/m激活块体尺寸/m爆炸当量/kt距爆心投影点距离/m激活块体尺寸/m
    2.3 60111.4 150 12
    12.511035 165 15
    1.41202517 430 80
    1110012 500 62
    3 9015 560 55
    10012 625 78
    11012 680 40
    1201178 800100
    13018 900100
    140151 020155
    1501054 750 80
    15510 800 70
    19010 900115
    1.4110161 000100
    125 81 100 90
    13514
    下载: 导出CSV

    表  5  3#测量区数据

    Table  5.   Results of light distance measurements at Site 3#

    试验试验日期距离D /m(D·Q−1/3)/(m·kt−1/3)位移S/mm(S·Q−1/3)/(mm·kt−1/3)
    13481987.08.021 70039819.54.56
    13501988.09.143 500674122.31
    13461988.12.173 900890143.19
    13521989.07.085 7001 74382.45
    14101989.09.02~7 0004 00021.07
    下载: 导出CSV

    表  6  4#测量区数据

    Table  6.   Results of light distance measurements at Site 4#

    试验编号试验日期距离D /m(D·Q−1/3)/(m·kt−1/3)位移S/mm(S·Q−1/3)/(mm·kt−1/3)
    13481987.08.023 7008666.661.56
    13501988.09.141 45027919.73.8
    13461988.12.176 6501 51851.14
    13521989.07.082 9008873.51.07
    下载: 导出CSV
  • [1] 徐开礼, 朱志澄. 构造地质学[M]. 2版. 北京: 地质出版社, 1989.
    [2] HILL D P, REASENBERG P A, MICHAEL A, et al. Seismicity remotely triggered by the magnitude 7.3 landers, California, earthquake [J]. Science, 1993, 260(5114): 1617–1623. DOI: 10.1126/science.260.5114.1617.
    [3] 钱七虎. 深部地下空间开发中的关键科学问题[R]. 南京: 解放军理工大学工程兵工程学院, 2004.
    [4] ADUSHKIN V V, OPARIN V N. From the alternating-sign explosion response of rocks to the pendulum waves in stressed geomedia: part iii [J]. Journal of Mining Science, 2014, 50(4): 623–645. DOI: 10.1134/S1062739114040024.
    [5] Adushkin V V, Spivak A. Underground explosions[R]. WGC-2015-03, 2015
    [6] ADUSHKIN V V, SPIVAK A A. Changes in properties of rock massifs due to underground nuclear explosions [J]. Combustion Explosion and Shock Waves, 2004, 40(6): 624–634. DOI: 10.1023/B:CESW.0000048263.34894.58.
    [7] RODINOV V, ADUSHKIN V, KOSTUCHENKO V, et al. Mechanical effect of an underground explosion [M]. Moscow: Nedra, 1971.
    [8] KOCHARYAN G G. Experimental investigations on the dynamic deformation of underground excavations in rock [C] // 2nd North American Rock Mechanics Symposium. USA: American Rock Mechanics Association, 1996.
    [9] KOCHARYAN G G, BRIGADIN I V, KARYAKIN A G, et al. Failure of underground workings in rock of block structure under the action of dynamic forces: I: experimental data on the failure mechanics of real rock under the action of powerful explosions [J]. Journal of Mining Science, 1994, 30(4): 370–378. DOI: 10.1007/BF02048182.
    [10] KOCHARYAN G G, SPIVAK A A. Movement of rock blocks during large-scale underground explosions: Part i: experimental data [J]. Journal of Mining Science, 2001, 37(1): 64–76. DOI: 10.1023/A:1016736919590.
    [11] KOCHARYAN G G, SPIVAK A A, BUDKOV A M. Movement of rock blocks during large-scale underground explosions: part ii: estimates by analytical models, numerical calculations, and comparative analysis of theoretical and experimental data [J]. Journal of Mining Science, 2001, 37(2): 149–168. DOI: 10.1023/A:1012327627277.
    [12] NIKOLAEV A V. Distant aftershocks of earthquakes and underground nuclear explosions [J]. Reports of Russian Academy of Science, 1999, 364(1): 110–113.
    [13] MCKEOWN F A, DICKEY D D. Fault displacements and motion related to nuclear explosions [J]. Bulletin of Engineering Geology and the Environment, 1969, 59(6): 2253–2269.
    [14] CLOUD W K, CARDER D S. Ground effects from the boxcar and benham nuclear explosions [J]. Bulletin of the Seismological Society of America, 1969, 59(6): 2371–2381.
    [15] GLASSTONE S, DOLAN PHILIP J. The effects of nuclear weapons [M]. USA: U.S. Atomc Energy Commission, 1977.
    [16] ENGDAHL E R. Seismic effects of the milrow and cannikin nuclear explosions [J]. Bulletin of Engineering Geology and the Environment, 1972, 62(6): 1411–1423.
    [17] WANG T, SHI Q, NIKKHOO M, et al. The rise, collapse, and compaction of Mt. Mantap from the 3 September 2017 North Korean nuclear test [J]. Science, 2018, 361(6398): 166–170. DOI: 10.1126/science.aar7230.
    [18] 陈宗基, 吴海青. 我国在复杂岩层中的巷道掘进:兼论构造应力与时间效应的重要性 [J]. 岩石力学与工程学报, 1988, 7(1): 1–001. DOI: 10.3321/j.issn:1000-324X.1999.03.022.

    CHEN Zongji, WU Haiqing. Yunnelling in complex rock formations in China: importance of geotectonic stress and time effects [J]. Chinese Journal of Rock Mechanics and Engineering, 1988, 7(1): 1–001. DOI: 10.3321/j.issn:1000-324X.1999.03.022.
    [19] KURLENYA M V, OPARIN V N. Problems of nonlinear geomechanics: part ii [J]. Journal of Mining Science, 2000, 36(4): 305–326. DOI: 10.1023/A:1026673105750.
    [20] 王明洋, 李杰, 邱艳宇, 等. 基于能量原理的大规模地下爆炸不可逆位移计算方法 [J]. 爆炸与冲击, 2017, 37(4): 685–691. DOI: 10.11883/1001-1455(2017)04-0685-07.

    WANG Mingyang, LI Jie, QIU Yanyu, et al. Calculation method for irreversible deformation region radius under large-scale underground explosion based on law of energy [J]. Explosion and Shock Waves, 2017, 37(4): 685–691. DOI: 10.11883/1001-1455(2017)04-0685-07.
    [21] 王明洋, 李杰. 地下核爆炸诱发工程性地震的计算原理及深地防护[J]. 岩石力学与工程学报, 已录用.

    WANG Mingyang, LI Jie. The calculation principle of engineering seismic effects induced by underground nuclear explosion and deep underground protection[J]. Chinese Journal of Rock Mechanics and Engineering, accepted.
    [22] KURLENYA M V, OPARIN V N, VOSTRIKOV V I. Pendulum-type waves: part ii: experimental methods and main results of physical modeling [J]. Journal of Mining Science, 1996, 32(4): 245–273. DOI: 10.1007/BF02046215.
    [23] KURLENYA M V, OPARIN V N, VOSTRIKOV V I, et al. Pendulum waves: part iii: data of on-site observations [J]. Journal of Mining Science, 1996, 32(5): 341–361. DOI: 10.1007/BF02046155.
    [24] 李杰, 陈伟, 施存程, 等. 基于块系构造的大规模地下爆炸不可逆位移计算方法 [J]. 爆炸与冲击, 2018, 38(6): 1271–1277. DOI: 10.11883/bzycj-2017-0201.

    LI Jie, CHEN Wei, SHI Cuncheng, et al. Calculation method of irreversible deformation region radius under large-scale underground explosion based on block hierarchical structure [J]. Explosion and Shock Waves, 2018, 38(6): 1271–1277. DOI: 10.11883/bzycj-2017-0201.
  • 期刊类型引用(4)

    1. 陆天然,段梦乔,李子怡,周连庆. 人工智能地震分类研究进展综述. 地球物理学进展. 2025(01): 25-47 . 百度学术
    2. Chengzhi Qi,Mingyang Wang,Gevorg Kocharyan,Artem Kunitskikh,Zefan Wang. Dynamically triggered seismicity on a tectonic scale: A review. Deep Underground Science and Engineering. 2024(01): 1-24 . 必应学术
    3. 李志浩,李杰,王明洋,蒋海明,吴红晓,潘越峰,陈伟,王德荣. 大当量爆炸地冲击毁伤效应的理论与试验研究Ⅰ:深埋洞室地冲击破坏的现场实测分析. 岩石力学与工程学报. 2022(05): 865-876 . 百度学术
    4. 李杰,王明洋,陈昊祥,徐天涵. 深部非线性岩石动力学的理论发展及应用. 中国科学:物理学 力学 天文学. 2020(02): 17-24 . 百度学术

    其他类型引用(0)

  • 加载中
图(14) / 表(6)
计量
  • 文章访问数:  5464
  • HTML全文浏览量:  1810
  • PDF下载量:  79
  • 被引次数: 4
出版历程
  • 收稿日期:  2018-12-17
  • 修回日期:  2019-01-19
  • 刊出日期:  2019-12-01

目录

/

返回文章
返回