Study on the influence of local mountainous topography to fault dynamic rupture
-
摘要: 基于曲线网格有限差分方法研究了垂直走滑断层在不同山体地形情况下的动力学破裂传播,分析并讨论了局部山体地形对断层破裂过程及相应地面地震动的影响,得到了各模型断层面的动力学破裂过程及相应的地表峰值速度特征。研究结果表明,山体地形尺度(山体高度及底部延展距离等)对断层动力学破裂过程影响较大,进而影响到相应的地面地震动分布。当山体地形处于自由地表上亚剪切向超剪切转换的位置附近时,山体地形会阻碍断层面上自由地表超剪切的产生。一般而言,对于具有一定埋深的断层,当山体地形底部延展距离一定时,山体高度越高,其对自由地表超剪切的阻碍程度越大;当山体高度一定时,地形底部延展距离越大,越会阻碍自由地表超剪切的产生,这种破裂过程的变化会导致相应地面地震动呈现不同特征的分布。此外,还探讨了断层破裂过程及相应地震动对成核区外初始剪切应力变化的响应,结果显示,当初始剪切应力较高时,高应力降引起的超剪切破裂会对断层破裂及相应的地震动分布起主导作用。Abstract: In this study, the curved grid finite-difference method was implemented to investigate the effect of local irregular topography on the dynamic rupture process of a vertical strike-slip fault and the resultant strong ground motions. The rupture propagation and ground motions were simulated with different irregular topography in a three-dimensional homogeneous half-space. Our results show that the scale of ridge topography including its height and bottom extension size has great impact on the dynamic rupture process, and then will affect the distribution of ground motions. The mountainous topography will obstruct the generation of super-shear induced by free surface when it is located near the subshear-to-supershear transition position on free surface. Generally, for the faults with a certain buried depth, with the same size of topography bottom extension, the higher the mountain height is, the stronger prevention it has on the generation of super-shear. In addition, when the mountain height is fixed, the larger extension of mountain bottom size has more obstacles to the generation of the super-shear induced by free surface. The variation of fault rupture process will make different distribution of ground motions. Furthermore, the response of dynamic rupture process and the corresponding ground motion to the change of initial shear stress outside the nucleation area was discussed. Our result shows that with the high initial shear stress, the super-shear induced by high stress drop also plays an important role in dynamic rupture and distribution of the resultant ground motion.
-
引言
井下地震观测可以减少地表的噪声干扰影响。井底背景噪声的频谱幅度约为地表的数百分之一,甚至千分之一(Takahashi,Hamada,1975),可有效地提高观测波形的信噪比(徐纪人,赵志新,2009)。井下地震仪可以在高噪声背景区域观测到更多的微小近震,近距离观测的地震波更多地保留了其中的高频成分,而高频成分正是解析地球内部物理、地壳和上地幔结构不可缺少的信息,对研究断裂活动、理解地质过程等都起着极为重要的作用(Aster,Shearer,1991;Yamauchi et al,2005)。井下观测地球物理状态的变化,有利于地震活动性的分析及震源机制和地震预报研究(Oye et al,2004)。目前井下地球物理观测研究已逐渐发展成为多学科的前沿性科学研究课题(Baisch et al,2002),其成果推进了现代地震学、地球物理学以及钻孔板块构造学等学科的发展(Zoback et al,2011)。由于强地面运动受到厚沉积物的显著影响,井下观测也是研究沉积平原地区强地面运动的有效途径(Zhao et al,2004)。我国井下地震观测研究始于20世纪70年代,至上世纪80年代中期已有百余个井下地震观测台,遍布全国(修济刚,1988)。
1. 先进的井下地震观测台网
近年来井下地震观测台站在数量及深度上都得到迅速发展,观测多使用国产地震仪,且观测质量优良。在我国1 000多个数字化测震台站中,井下观测台约200个(张明等,2019)。2001年建成的首都圈测震台网的107个台站中,58个观测台站配备有井下地震仪(张尉等,2009),已发展成为先进的井下地震观测台网,促进了我国地震事业的发展。图1给出了我国大陆地区表1中部分井下观测台站的分布。其分布与地震构造带关系密切,主要分布在我国东部沉积层覆盖比较厚的地区,目的在于监控地震活动及开展地震研究工作。增设井下观测台站,形成科学的全方位观测台网,可以提高地震观测研究能力。为了监视大断裂的构造活动和地震活动,我国在西部一些地区也增设了井下地震观测研究台站。
表 1 我国部分井下观测台站基本情况Table 1. Basic information of some underground observation stations in China台站名 井深/m 地震仪型 参考文献 台站名 井深/m 地震仪型 参考文献 首都圈台网 150—480 宽、甚频 短周期 刘渊源等(2 011) 上海东滩 421 FSS-3DBH 裴晓等(2 012) 北京市台网 宽、甚频 短周期 兰从欣等(2 005) 上海金泽 305 FSS-3DBH 首都圈白家疃 257 韦士忠和李玉萍(1 990) 上海上戏 370 FSS-3DBH 首都圈文安 266 上海南汇 280 JDF-2 首都圈东三旗 250 上海大新中学 375 FSS-3DBH 首都圈大兴 110 上海竹园 317 JDF-1 首都圈雄县 358 上海虹桥 651 768 首都圈龙门庄 447 上海八角厅 780 JDF-1 吉林松原 384 FSS-3DBH 韦庆海等(2 015) 上海张江 350 KS-2000M 裴晓等(2 013) 吉林松原 243 FSS-3DBH 陈闯等(2 022) 上海崇明 463 李伟等(2 013) 大庆新台 708 JD-2 韦庆海等(2 015) 江苏宝应 460 CMG-3TB 仇中阳等(2 014) 河北赵县 260 BBVS-60 郑德高等(2 018) 江苏高邮 440 CMG-3TB 河北唐海 480 短周期 郑德高等(2 018) 江苏淮安 315 JDF-2 河北涿县 320 768 李彦林和郑淑兰(1 989) 江苏涟水 400 CMG-3TB 河北邯郸 400 JD-2 张新东(2 002) 江苏射阳 380 CMG-3TB 河北肥乡 400 JD-2 江苏盐城 445 CMG-3TB 河北临漳 440 JD-2 江苏金湖 447 GL-S60B 宫杰等(2 019) 天津静海 371 768 赵惠君等(1 991) 江苏滨海 470 GL-S60B 天津芦台 276 768 江苏丹阳 175 GL-S60B 天津武清 450 768 江苏响水 410 GL-S60B 内蒙赤峰 90 GL-S120B 郭延杰等(2 020) 江苏高邮 458 GL-S60B 甘肃天水 337 短周期 蔡耐芳(1 990) 江苏建湖 443 GL-S60B 新疆喀什 283 GL-S60B 赵瑞胜等(2 021) 江苏启东 410 GL-S60B 山西太原 500 JD-2 张少泉等(1 988) 江苏盐城 436 GL-S60B 陕西定边 300 BBSV-60BH 李少睿等(2 016) 江苏泰兴 425 GL-S60B 宁夏灵武 248 JDF-2 江苏东台 458 GL-S60B 宁夏陶乐 245 JDF-2 江苏溧阳 83 CMG-DM24 mk3 胡米东(2 014) 河南安阳 393 FSS-3DBH 江苏盐城 445 CMG-DM24 mk3 河南清丰 308 FSS-3DBH 江苏南通 105 CMG-DM24 mk3 四川泸州 95 CMG-3TB 江苏大丰 366 短周期 徐元耀(1 994) 山东荷泽 370 768 周焕鹏(1 986) 江苏淮阴 325 井下摆 安徽六安 126 GL-S60B 石英杰等(2 021) 江苏海安 420 JD-2 安徽霍邱 150 GL-S60B 石英杰等(2 021) 云南昆明 2 02 GL-S60B 李雷等(2 018) 浙江景宁 68 FSS-3DBH 张明等(2 019) 云南昆明 452 JD-2 修济刚(1 988) 浙江北仑 86 FSS-3DBH 云南大寨 375 井下仪 王芳等(2 017) 浙江南麂岛 110 GL-60DBH 广东汕头 2 00 TBG-60B 郭德顺等(2 014) 上海普陀 564 768 叶世元和柳国华(1 987) 江苏东海 4050 宽频 短周期 Xu等(2 016) 上海海运 600 768 叶世元和柳国华(1 987) 1.1 先进的井下地震观测技术和观测仪器
我国井下台网的观测技术已日趋完善成熟。在井下密封、井下定向、方位补偿、电缆应力解除等观测技术领域也取得了令人满意的进步。而适合不同井孔孔径的锁固装置的成功研制及利用反馈控制技术调节井下长周期地震仪的零位漂移、振幅效应等技术,可保持井下地震仪长期稳定工作(高龙生,1986)。通过编码遥控装置调整井下的地震仪摆体的夹紧、放松状态以及标定等多种功能,可改进井下地震仪的工作状态(陈宏水,付迺珍,1987)。规范井下地震仪安装过程以保证其良好的观测能力(胡米东,2014);规范安装观测房和井口之间的防雷设施及线路布设,使井下地震台站稳定可靠地运行(陈吉锋等,2014)。根据条件选择落底安装井下地震仪,可以提高地震仪的观测能力(宫杰等,2019)。利用波形相关方法,分析地表与井下地震仪方位角检测以及一致性对比测试,使井下地震仪方位角检测精度优于4° (李少睿等,2016)。在缆线应力释放完成后固定井下地震仪,同地表地震仪进行数据相关性分析,以确定方位角的可靠性(张明等,2019)。
井下观测多使用我国研制的先进的短周期、宽频和甚宽频带地震观测仪器(刘瑞丰等,2008)。井下观测资料的分析方法也在不断精益求精。对比分析各种地震观测参数,研究人员提出了高效、可靠的识别与提取井下微地震信号的方法(李稳等,2016)。利用我国井下台网已经可以进行全频域、高精度、长期稳定的地震监测研究,井下台网也成为获取可靠资料的保障,为地震科学的发展提供了更坚实的基础。表1列出了本文涉及的井下观测台站的深度、观测地震仪的频段特性以及观测成果来源。
1.2 井下观测网提高了地震观测能力
上述先进的观测仪器和井下观测技术提高了地震观测能力和科学研究水平。地面背景噪声常导致地震波形很难被地面地震台识别,采用井下摆观测方式就可以有效抑制各种人为活动噪声的干扰。宽频地震仪安装在地下10 m的中微风化层基岩上,可有效减小噪声(郭德顺等,2014)。研究平原地面、井下与山区地震的纵、横波的运动学和动力学特征得知,井下记录比较接近设置于基岩上的地震台的记录。井下记录的地震波形的噪声功率均值较小(张阿瑶等,2018),可以更精确地记录到地脉动(裴晓等,2013)。上海张江井下台的背景噪声的均方根比地表台站要小1个数量级,而且近震、远震波形记录都很清晰。井下摆的放大倍数亦可达地面摆的十倍多(李彦林,郑淑兰,1989)。井下观测的数据精度也能提高1—2个数量级,并且井下观测系统的有效动态范围超出地面系统约10%—30%,能记录到的地震动信号的范围更宽(叶世元等,1988)。井下90 m深地震仪也能观测到比山洞中深处花岗岩台址上地震仪更远的同等大小的地震(郭延杰等,2020)。
东海台阵深井地震仪可以检测几十至几百千米外的微小地震的清晰波形,而地表记录却识别不出这些地震的迹象(表2)。图2示出了东海台阵的日常监测屏上不同深度地震仪显示的2019年9月22日山东烟台ML3.1地震纪录图。图中地表记录识别不出地震迹象,井下400 m深处记录可看出明显地震波形。而且,地下3 km深处记录也可看出清晰的地震波形,表明监控记录可靠。井下观测填补了在高噪声区域获得高精度地震波形的空白,在提高地震监测分辨能力及观测地震波精度等方面具有明显优势(徐元耀,1994;徐纪人,赵志新,2006)。
表 2 江苏东海垂直地震台阵地面与井下观测能力比较Table 2. Comparison of observation ability between surface and underground of vertical seismic array in Jiangsu DonghaiML 震中距
Δ/km地面观测宽频
带地震仪井下观测 井深/m 仪器 观测效果 −1.3 64 无法识别 2545 短周期 图像清晰 −0.5 71 无法识别 2545 短周期 图像清晰 0.8 97 无法识别 3500 短周期 图像清晰 3.1 492 无法识别 400 宽频 图像清晰 图 2 东海垂直地震台阵的地面及井下仪器记录的2019年9月22日烟台ML3.1地震G0:地表宽频地震仪;G1:400 m深宽频地震仪;G3:3 km深短周期地震仪Figure 2. The September 22,2019,Yantai ML3.1 earthquake recorded by surface and underground instruments of Donghaivertical seismic array G0:Surface broadband seismograph;G1:400 m deep broadband seismograph; G3:3 km deep short period seismometer2. 井下观测研究对地震科学研究发展的重要贡献
井下地震观测台网记录到的地震波形,噪声小且精度高,这些可靠的高质量资料较精确地反映了地震波路径及震源信息。在远离地面干扰环境中,长期观测地球内部的信息,所得结果客观真实。用其高分辨率记录测定地震参数,研究震源过程、地震成像和地球内部结构,可以改善研究成果的精度,有利于创新地震波理论研究,且目前已取得了可喜的新成果。
2.1 井下观测得到的精确地震参数推进了地震学发展
井下观测研究避免了地表噪声和场地效应的影响,可全方位布局地震观测台,这样增用井下观测资料得到的地震定位结果更加准确(仇中阳等,2014),进而可推进对地震序列活动及震源位置的研究,以及断层构造运动等的研究。分析井下摆与地面摆测定的地震参数差异成因,并给出两者之差的适当校正值(朱雅山,李永勤,1992),可提高地震观测成果的客观性与可靠性。
井下观测地震波形噪声小,有利于高精度测定地震P波和S波到时、振幅、周期等参数(赵惠君等,1991),推进了地震学的研究。上海南汇台用统计方法分析了井下记录地震尾波衰减与震级之间的关系,求得了垂直向地震波持续时间与对数震级Mv在近震范围内(Δ<5°)呈线性关系(李慧民等,1992)。用井下地震仪记录的最大振幅测定近震震级,误差分析表明,结果合理可靠(章纯等,1992)。吉林地区台站由井下地震仪得到的近震震级比地表地震仪的小0.12 (陈闯等,2022)。河北涿县井深350 m的短周期摆与地面同型摆记录的地震波形对比研究表明,地面S波两水平向最大地动位移平均为井下摆的4.7倍。地面摆垂直向最大地动位移平均为井下摆的2.7倍。因此,井下地震摆观测所得震级比地面摆平均低0.7。如采用波形振动持续时间测定震级,井下记录可得到与地面记录较为一致的结果。涿县地震台P波首波到时,地面地震仪摆要比井下摆滞后0.3 s (张寿康等,1986)。大庆井下地震台网不同深度各子台观测结果表明,近震震级偏差和校正与井下摆深度和介质性质等有关(袁卫红等,2012)。
新疆喀什麦盖提台站井下宽频地震仪记录的近、远震波形比地面的更为清晰。井下地震仪得到的近震和远震震级分别比地表地震仪得到的震级小0.518和0.025,且功率谱也相对较低(赵瑞胜等,2021)。昆明基准地震台山洞内地面上的宽频带地震仪观测结果ML比井下结果平均大0.18 (李雷等,2018)。宽频带地震仪井下记录的地震波的运动学参数、地震波走时比地面的小,统计得出地面地震台的动力学参数地动振幅大约是井下振幅的10倍(王俊国等,1988),与上述短周期地震仪研究结果类似。地震矩和矩震级观测也表明井下震级较地面要小(仇中阳等,2014)。这可能是由于地面记录存在较大的地震波振幅效应,某些情形下地面振幅相对于井下振幅的放大因子可能大于1。井下观测结果是相对客观而科学的。图3显示出部分井下与地面观测的同一地震事件的震级比较。图中红线为横轴与纵轴等值参考线,横轴和纵轴数值分别表示地面地震仪与井下地震仪观测到同一事件的震级。可看出井下震级明显小于地面震级,特别是对于M4.0以下的小地震(毛华锋等,2014;仇中阳等,2014,赵瑞胜等,2021)。目前井下观测资料在分析地震参数、进行地震速报及监测地震活动等方面日益发挥出重要作用(修济刚,1988),推进了深井地震学研究的发展。
2.2 促进了地震波速度结构的研究
低噪音的井下地震观测资料已用于地球内部结构研究,有益于建立高精度三维地震波速度模型(Lay et al,2020)。地表和井下台站波形资料的联合利用提高了监控区域地震波速研究精度,揭示了京津唐地区地震波速度结构明显的横向不均匀性。并且这些资料也被用于探索地震成因及地震预报,发现许多地震均发生在P波低速异常体与高速异常体的交界地带(于湘伟等,2003)。唐山震区的三次强震均发生在低速体与高速体之间(齐诚等,2006)。地震波速时间的变化显示2006年文安MS5.1大地震前一年,文安和唐山地区波速比明显出现了可信的前兆性降低过程(王林瑛等,2008)。
井下地震观测记录到了体波的地表反射震相。这一发现促进了地震学观测的发展。井下地震仪记录的地震波震相较之地面记录更为丰富(乌统昱,赵惠君,1992)。大庆台井下708 m地震仪记录到2005年7月25日林甸M5.1地震和2003年6月26日松原M2.3地震的地表反射震相。用井下接收的地面反射震相PP,SP等可以确定井下观测站与地面之间的速度值和界面埋深(林云松,1989)。静海台井下地震仪记录到周围微震波形的4个主要震相,即P,S,地面反射波PP和SS,其中SS震相的水平分量通常较强。松原台384 m深处井下摆记录到许多地方震的地表反射波(韦庆海等,2015)。用地表反射波资料求得松原台地表到井下地震仪之间的P波速度约为2 km/s (陈闯等,2022)。新疆喀什台站用井下资料求得麦盖提下方浅层的平均S波速度为0.59 km/s (赵瑞胜等,2021)。由此结果构成的地表沉积层速度模型可有效改善地震定位精度。利用反卷积方法识别近震入射波与地表反射波,可建立高精度浅层地震波速度模型(王芳等,2017;Riga et al,2019)。
井下摆反射波形分析是一种测量深达数百米S波速度的有效方法。分析首都圈台网中58个井下摆记录的100多次M≥3.0地震的直达波及对应地表反射波的波形可知,这些台站上方厚约300 m浅地表土层P波和S波平均速度分别约为1.6—2.0 km/s和0.34—0.48 km/s,波速比约为4.0—5.3。该区域P波速度和波速比均无显著横向变化。清晰的震相图也表明,土层没有明显的速度断层(刘渊源等,2011)。首都圈地表S波平均速度梯度为0.8 m/(s·m),浅部100 m的平均S波速小于300 m/s,深度500 m处S波速度为800 m/s (沈伟森等,2010)。井下摆近震记录研究表明,华北盆地近地表P波平均速度约为1.98 km/s,S波平均速度约为0.46 km/s (朱音杰等,2020)。上述地震波速度成果对计算地球内部精细结构、地震波传播在地球内部的衰减特征以及资源评估均具有重要意义,同时也可为研究场地效应提供直接证据(罗诚等,2018)。
2.3 丰富了震源及地球内部物理学研究成果
研究地震记录的背景噪声功率谱可知,井下地震台在全频范围内记录数据质量明显优于地表地震仪(石英杰等,2021)。井下台网高质量的资料为研究震源和地球内部物理学奠定了坚实的基础,取得了大量可信赖的高精度成果。
利用井下和地面观测资料求得了京津唐地区1 132次ML1.7—6.2地震的精确震源参数(于湘伟等,2003),促进了震源过程研究。用最小二乘法拟合北京地区震源参数与震级的定量表达式,证实拟合所使用资料的台网规模对震源参数结果影响程度的差异。场地响应对震源参数的计算结果影响也比较大(兰从欣等,2005)。江苏部分地震台井下观测波形拐角频率比地面观测的小0.397 Hz,井下观测的地震矩和矩震级分别比地面的小8.642×1014 N·m和0.32,断层错距和应力降分别比地面的小2.268 cm和2.033 MPa (毛华锋等,2014)。仇中阳等(2014)的结果显示,井下记录地震波拐角频率较地面约小0.28 Hz;地震矩和矩震级数值较地面分别小3.104×1014 N·m和0.32。井下测定的地震波谱高频较弱,拐角频率较低,地震矩、应力降和平均位错都比地面基岩测定的结果偏低,而震源等效半径则偏高,井下记录高频成份没有地面的丰富(韦士忠,李玉萍,1990)。接收函数正演计算和井下理论波形研究发现,井下地震仪记录P波波形垂向分量的频谱极小值随深度的增加而减小(郑德高等,2018)。首都圈西北部地区快剪切波平均偏振方向为NE方向,它暗示了该区域的主压应力水平方向;慢剪切波时间延迟急剧的梯度变化可能与地壳深部的温度变化有关(吴晶等,2007)。上述结果对于探索地震发生机制有参考意义。加之钻孔井下记录受场地效应的影响比地表地震仪小,可为研究震源提供更直接的约束(Shearer,Abercrombie,2021),其成果丰富了现代地震学研究宝库。
2.4 井下地震观测与强地面震动预测
强地面运动的振幅和持续时间受到厚沉积物速度结构的显著影响。研究地震波传播在沉积物中的非线性效应,以及近地表地震波衰减和场地效应对于精确预测强地面运动至关重要(Stephenson et al,2005)。
将短周期井下摆与地面仪记录对比可知,地表土层对地震波有放大作用。北京西部地区地表基岩台的场地响应总体上比东部平原区井下台的响应大,其低频端场地响应变化平稳。东部平原区井下台的场地响应也与频率有关(兰从欣等,2005)。品质因子Q值与频率的定量关系显示,台站场地响应在低频段幅度变化亦保持稳定(蔡耐芳,1990;卞真付等,2005)。将地面和井下观测场地效应资料用于波场理论计算,可减少波场解的不确定性,为工程地震提供更可靠的依据(张少泉等,1992)。
井下强震动观测台阵中加速度仪记录是研究强地面运动的重要方法(Hollender et al,2023)。地表与井下实测数据可有效地用于研究各种土层的地震响应(丁海平等,2014)。进一步反演结果指出,土层的S波和P波速度及厚度是影响场地响应的水平与竖向谱比的重要因素(李小军等,2020)。由井下基岩记录的地震波作为输入研究土层地震动响应可知,土层对地震波的放大效应与波形主频率有关,放大倍数最高可达约20倍。该结论对软土地区的烈度区划、建筑的抗震设计有现实意义(徐永林等,2002)。
3. 东海大陆深井垂直地震台阵
深井地震垂直台阵观测研究有益于地壳活动性的监测及地震成因的研究,开创了高精度地震学研究的新途径(Ma,2021)。台湾车笼埔断层带上的井下垂直地震台阵中,多台波形水平分量之间的尾波互相关函数的研究结果显示,该区地震各向异性呈现出季节性变化(Hung et al,2022)。
江苏东海大陆深孔地壳活动国家野外科学观测研究站在国际大陆科学钻探项目(ICDP)之一的中国大陆科学钻探的
5000 m深井中,建设了我国大陆第一个深井垂直地震台阵,即东海大陆深井垂直地震台阵(简称东海垂直地震台阵)。除设立地面观测地震仪之外,还在井下不同深度设置了4层地震观测仪。在相距深井约500 m的另一个450 m井中进行了综合地球物理观测。2008年汶川MS8.0地震后,在汶川地震区跨断层建立了6口200—3000 m不同深度的井下观测台,深入研究大地震发生机制,监测断层活动。东海大陆深井垂直地震台阵及其周围处于高噪声背景区域,致使地表地震仪记录中难于识别ML0.8以下地方震波形。深井地震台阵能够观测到比地面地震仪更多的零级或负震级微小地震的清晰波形,观测资料清晰可靠。井下不同深度地震台观测的地震序列的波形信噪比随时间变化的趋势基本一致。其地震波形基本达到了70 dB以上的高信噪比。在高噪声背景区域获得了高保真度的微小地震波形(徐纪人等,2022),提高了监测地壳活动和地震活动的能力。不同深度处的记录波形的信噪比可为井下观测台站建设提供科学参考依据。低噪声高质量的地震波形研究成果可以加深对地震运动学和动力学的理解;建立精细波速模型,有助于全面清晰地揭示地球内部结构,创新地震波传播理论研究(Ma,2021)。分析低噪声的微小地震至大地震的波相特性,可为震源和震源谱的研究提供直接约束,更深刻地解释震源破裂过程和机制。利用深井地震波形研究场地效应可以提高强地面运动预测水平。未来深井垂直地震台阵观测研究或可在高精度地震学研究中发挥巨大作用。
4. 总结和展望
本文研究表明,井下观测可以在高噪声背景区域近距离观测地震活动,应科学地布局观测台站,以取得高质量地震资料,进而推进地震科学的发展。用井下资料研究波速空间分布的非均匀性,可以缩小因地面台站台基的不同而导致的各台站间的P,S波的平均速度差异,提高所建立的区域波速模型的可信度和科学性。井下观测与地面观测的震级、矩震级比地面基岩测定的结果偏低(图3),其差异可能与井下仪器上层介质的波形非线性增幅效应与波的频率有关(Zhao et al,2004)。井下观测拐角频率低于地表可能与井下仪器上层介质对不同频率波分量吸收与放大有关(周焕鹏,1986)。场地响应对震源谱参数影响也较大,在低频域,绝大多数井下台的场地响应都大于1,高频端,绝大多数井下台的场地响应都小于1(兰从欣等,2005),这些因素都可能导致井下台站观测的震级、拐角频率都低于地面观测结果。上述差异的成因也具有多重复杂性(Abercrombie,1997),仍是今后科学探索的重要课题。井下地震观测开创了我国在地下深处观测研究地球内部的新途径,提高了地震观测能力,为地球科学研究作出了巨大贡献;对资源评估、地震预报、防震减灾也具有重大的现实意义(徐纪人等,2022)。本文未能引用全国的数百个井下观测台的全部累累成果,这些论文和成果可以在相关的文献中见到。
深井地震观测研究是现今国内外的前沿课题,提高其观测研究水平可实现高精度地震学研究,而高精度地球物理学正是创新未来的必然发展轨迹(Zeng,Yang,2021)。我国还将要沿安宁河断裂带建设12个井深300 m和两个井深
1000 m的综合地震观测研究井(王晓蕾等,2021);目前已经建成8个井深300 m和两个井深1000 m的综合地震观测台阵,其中湖南2000 m深井宽频地震观测站已取得重大成果。东海大陆深孔地壳活动国家野外科学观测研究站将继续增强汶川及东海地区井下综合地球物理长期观测研究工作。随着井下台网观测研究工作的发展,我国的高精度科学研究成果将更上一层楼。 -
图 1 三维断层及地形模型 (a)与断层面所在垂直剖面 (b)示意图
蓝色矩形表示埋深为200 m的垂直走滑断层,黄色方形表示成核区,灰色区域表示包含吸收边界的高强度区
Figure 1. Three dimensional fault and topography model (a) and vertical profile along fault plane (b)
The blue rectangle depicts the 3-D vertical strike-slip fault,the yellow square indicates the nucleation area and the grey are aindicates a high strength area including absorbing boundary
图 2 成核区外初始剪切应力为5 MPa时各模型断层面上的破裂起始时间等值线分布
图(a)为水平自由地表模型;图(b)−(f)为山体地形高度分别为400,600,800,1 000和1 100 m时的模型
Figure 2. Initial time contours of ruptures on fault plane for different models with initial shear stress 5 MPa outside the nucleation area
Fig. (a) is for flat free surface model;Figs. (b)−(f) are for Gaussian hill models with height 400,600,800,1 000 and 1 100 m,respectively
图 3 成核区外初始剪切应力为5 MPa时各模型断层面上峰值滑动速率分布
图(a)为水平自由地表模型;图(b)−(f)为山体地形高度分别为400,600,800,1 000和1 100 m时的模型
Figure 3. Distribution of peak slip rate on fault plane for different models with initial shear stress 5 MPa outside of the nucleation area
Fig. (a) flat free surface model;Figs. (b)−(f) are for Gaussian hill models with height 400,600,800 ,1 000 and 1 100 m,respectively
图 4 成核区外初始剪切应力为5 MPa时各模型地表峰值速度的平行断层面水平分量分布
图(a)为水平自由地表模型;图(b)−(f)为山体地形高度分别为400,600,800,1 000和1 100 m时的模型
Figure 4. Fault-parallel component of peak ground velocity distribution for different models with initial shear stress 5 MPa outsidethe nucleation area
Fig. (a) is for flat free surface model;Figs. (b)−(f) are for Gaussian hill models with height 400,600,800,1 000 and 1 100 m,respectively
图 5 成核区外初始剪切应力为5 MPa时各模型地表峰值速度的垂直断层面水平分量分布
图(a)为水平自由地表模型;图(b)−(f)为山体地形高度分别为400,600,800,1 000和1 100 m时的模型
Figure 5. Fault-normal component of peak ground velocity distribution for different models with initial shear stress 5 MPa outside the nucleation area
Fig. (a)is for flat free surface model;Figs. (b)−(f) are for Gaussian hill models with height 400,600,800,1 000 and 1 100 m,respectively
图 10 成核区外初始剪切应力为5.6 MPa时各模型断层面上破裂起始时间分布
图(a)为水平自由地表模型;图(b)−(f)为山体地形高度分别为400,600,800,1 000和1 100 m时的模型
Figure 10. Initial time contours rupture on fault plane for different models with initial shear stress 5.6 MPa outside the nucleation area
Fig. (a) is for flat free surface model;Figs. (b)−(f) are for Gaussian hill models with height 400,600,800,1 000and 1 100 m,respectively
图 11 成核区外初始剪切应力为5.6 MPa时各模型地表峰值速度的平行断层面水平分量分布
图(a)为水平自由地表模型;图(b)−(f)为山体地形高度分别为400 ,600,800,1 000和1 100 m时的模型
Figure 11. Fault-parallel component of peak ground velocity distribution for different models with initial shear stress 5.6 MPa outside of the nucleation area
Fig. (a) is for flat free surface model;Figs. (b)−(f) are for Gaussian hill models with height 400,600,800,1 000 and 1 100 m,respectively
图 12 成核区外初始剪切应力为5.6 MPa时各模型地表峰值速度的垂直断层面的水平分量分布
图(a)为水平自由地表模型;图(b)−(f)为山体地形高度分别为400,600,800,1 000和1 100 m时的模型
Figure 12. Distribution of fault-normal component of peak ground velocity for different models with initial shear stress 5.6 MPa outside the nucleation area
Fig. (a) is for flat free surface model;Figs. (b)−(f) are for Gaussian hill models with height 400,600,800,1 000 and 1 100 m,respectively
表 1 断层面的模型参数
Table 1 Model parameters on fault plane
区域 初始剪切应力σ0/MPa 剪切破裂强度σu/MPa 残余应力σr/MPa 临界滑动弱化距离Dc/m 断层成核区内 10.1 10.0 0 0 断层成核区外 5.0 10.0 0 0.2 断层外 5.0 200.0 0 20.0 -
胡进军. 2009. 近断层地震动方向性效应及超剪切破裂研究[D]. 哈尔滨: 中国地震局工程力学研究所: 178. Hu J J. 2009. Directivity Effect of Near-Fault Ground Motion and Super-Shear Rupture[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration: 178 (in Chinese).
王铭锋,郑傲,章文波. 2017. 局部山体地形对强地面运动的影响研究[J]. 地球物理学报,60(12):4655–4670. doi: 10.6038/cjg20171210 Wang M F,Zheng A,Zhang W B. 2017. Effect of local mountain topography on strong ground motion[J]. Chinese Journal Geophysics,60(12):4655–4670 (in Chinese).
张伟. 2006. 含起伏地形的三维非均匀介质中地震波传播的有限差分算法及其在强地面震动模拟中的应用[D]. 北京: 北京大学: 7–109. Zhang W. 2006. Finite Difference Seismic Wave Modelling in 3D Heterogeneous Media With Surface Topography and Its Implementation in Strong Ground Motion Study[D]. Beijing: Peking University: 7–109 (in Chinese).
Andrews D J. 1976a. Rupture propagation with finite stress in antiplane strain[J]. J Geophys Res,81(20):3575–3582. doi: 10.1029/JB081i020p03575
Andrews D J. 1976b. Rupture velocity of plane strain shear cracks[J]. J Geophys Res,81(32):5679–5687. doi: 10.1029/JB081i032p05679
Aochi H,Ide S. 2014. Ground motions characterized by a multi-scale heterogeneous earthquake model[J]. Earth Planets Space,66:42.
Burridge R. 1973. Admissible speeds for plane-strain self-similar shearcracks with friction but lacking cohesion[J]. J Geophys Res,35(4):439–455.
Das S,Aki K. 1977. A numerical study of two-dimensional spontaneous rupture propagation[J]. Geophys J Rastr Soc,50(3):643–668. doi: 10.1111/j.1365-246X.1977.tb01339.x
Day S M,Dalguer L A,Lapusta N,Liu Y. 2005. Comparison of finite difference and boundary integral solutions to three-dimensional spontaneous rupture[J]. J Geophys Res,110(B12):B12307. doi: 10.1029/2005JB003813
Dieterich J H. 1972. Time-dependent friction in rocks[J]. J Geophys Res,77(20):3690–3697. doi: 10.1029/JB077i020p03690
Dieterich J H. 1978. Preseismic fault slip and earthquake prediction[J]. J Geophys Res,83(B8):3940–3948. doi: 10.1029/JB083iB08p03940
Duan B C,Oglesby D D. 2007. Nonuniform prestress from prior earthquakes and the effect on dynamics of branched fault systems[J]. J Geophys Res,112(B5):B05308.
Dunham E M. 2007. Conditions governing the occurrence of supershear ruptures under slip-weakening friction[J]. J Geophys Res,112(B7):B07302.
Ely G P,Day S M,Minster J B. 2010. Dynamic rupture models for the southern San Andreas fault[J]. Bull Seismol Soc Am,100(1):131–150. doi: 10.1785/0120090187
Frankel A. 2004. Rupture process of the M7.9 Denali fault,Alaska,earthquake:Subevents,directivity,and scaling of high-frequency ground motions[J]. Bull Seismol Soc Am,94(6B):S234–S255. doi: 10.1785/0120040612
Harris R A,Archuleta R J,Day S M. 1991. Fault steps and the dynamic rupture process:2-D numerical simulations of a spontaneously propagating shear fracture[J]. Geophys Res Lett,18(5):893–896. doi: 10.1029/91GL01061
Huang H Q,Zhang Z G,Chen X F. 2018. Investigation of topographical effects on rupture dynamics and resultant ground motions[J]. Geophys J Int,212(1):311–323. doi: 10.1093/gji/ggx425
Ida Y. 1972. Cohesive force across the tip of a longitudinal-shear crack and Griffith’s specific surface energy[J]. J Geophys Res,77(20):3796–3805. doi: 10.1029/JB077i020p03796
Kaneko Y,Lapusta N. 2010. Supershear transition due to a free surface in 3-D simulations of spontaneous dynamic rupture on vertical strike-slip faults[J]. Tectonophysics,493(3/4):272–284.
Langer S,Olsen-Kettle L,Weatherley D. 2012. Identification of supershear transition mechanisms due to material contrast at bimaterial faults[J]. Geophys J Int,190(2):1169–1180. doi: 10.1111/gji.2012.190.issue-2
Lozos J C,Oglesby D D,Brune J N. 2013. The effects of fault stepovers on ground motion[J]. Bull Seismol Soc Am,103(3):1922–1934. doi: 10.1785/0120120223
Madariaga R. 1976. Dynamics of an expanding circular fault[J]. Bull Seismol Soc Am,66(3):639–666.
Madariaga R,Olsen K B,Archuleta R J. 1998. Modeling dynamic rupture in a 3D earthquake fault model[J]. Bull Seismol Soc Am,88(5):1182–1197.
Oglesby D D. 2008. Rupture termination and jump on parallel offset faults[J]. Bull Seismol Soc Am,98(1):440–447. doi: 10.1785/0120070163
Oglesby D D,Day S M. 2001. Fault geometry and the dynamics of the 1999 Chi-Chi (Taiwan) earthquake[J]. Bull Seismol Soc Am,91(5):1099–1111.
Oglesby D D,Archuleta R J. 2003. The three-dimensional dynamics of a nonplanar thrust fault[J]. Bull Seismol Soc Am,93(5):2222–2235. doi: 10.1785/0120020204
Ruina A. 1983. Slip instability and state variable friction laws[J]. J Geophys Res,88(B12):10359–10370. doi: 10.1029/JB088iB12p10359
Scholz C,Molnar P,Johnson T. 1972. Detailed studies of frictional sliding of granite and implications for the earthquake mechanism[J]. J Geophys Res,77(32):6392–6406. doi: 10.1029/JB077i032p06392
Wen Y Y,Oglesby D D,Duan B C,Ma K F. 2012. Dynamic rupture simulation of the 2008 MW7.9 Wenchuan earthquake with heterogeneous initial stress[J]. Bull Seismol Soc Am,102(4):1892–1898. doi: 10.1785/0120110153
Yamashita T. 1976. On the dynamical process of fault motion in the presence of friction and inhomogeneous initial stress partⅠ : Rupture propagation[J]. J Phys Earth,24:417–444. doi: 10.4294/jpe1952.24.417
Zhang H M,Chen X F. 2006a. Dynamic rupture on a planar fault in three-dimensional half-space:Ⅱ . Validations and numerical experiments[J]. Geophys J Int,167(2):917–932. doi: 10.1111/gji.2006.167.issue-2
Zhang W,Chen X F. 2006b. Traction image method for irregular free surface boundaries in finite difference seismic wave simulation[J]. Geophys J Int,167(1):337–353. doi: 10.1111/gji.2006.167.issue-1
Zhang W B,Iwata T,Irikura K. 2006. Dynamic simulation of a dipping fault using a three-dimensional finite difference method with nonuniform grid spacing[J]. J Geophys Res,111(B5):B05301.
Zhang W,Zhang Z G,Chen X F. 2012. Three-dimensional elastic wave numerical modelling in the presence of surface topography by a collocated-grid finite-difference method on curvilinear grids[J]. Geophys J Int,190(1):358–378. doi: 10.1111/gji.2012.190.issue-1
Zhang Z G,Zhang W,Chen X F. 2014. Three-dimensional curved grid finite-difference modelling for non-planar rupture dynamics[J]. Geophys J Int,199(2):860–879. doi: 10.1093/gji/ggu308