Impact of complete Coulomb failure stress changes of the 2008 Xinjiang Yutian MS7.3 earthquake on the subsequent earthquakes
-
摘要: 在离散波数法基础上计算2008年3月21日新疆于田MS7.3地震造成的近场区域完全库仑应力变化, 分析该变化对余震发生所产生的影响, 得到了此次地震在2008年10月5日乌恰MS6.8地震震中处所产生的动态库仑破裂应力变化. 计算结果表明, 该地震近场区域库仑应力变化图像演化大概持续了60 s, 库仑应力变化对余震的触发率达到90%以上, 其中动态库仑应力变化图像更好地解释了余震的分布. 余震震中处的完全库仑应力变化计算结果表明, 其动态库仑应力变化远远大于静态库仑应力变化. 于田MS7.3地震在乌恰MS6.8地震震中处造成的最大动态库仑应力变化为0.12 MPa, 说明后者可能受到了于田MS7.3地震的动态应力触发作用, 但不显著;而静态库仑应力则对其影响很小.
-
关键词:
- 于田地震 /
- 完全库仑破裂应力变化 /
- 动态应力触发 /
- 静态应力触发
Abstract: Based on the discrete-wavenumber method, we calculated the complete Coulomb failure stress (CFS) changes of near-field region induced by 21 March 2008, Yutian, Xinjiang, MS7.3 earthquake, and analyzed their impact on the aftershocks. We also get the dynamic CFS changes of the Yutian earthquake produced at the epicentral location of the Wuqia MS6.8 earthquake on 5 October 2008. The calculated results show that the evolution of CFS change patterns lasted for almost 60 s, and trigger ratio of aftershocks by the main shock is above 90%. Dynamic CFS change patterns can better explain the distribution of aftershocks. Complete CFS changes results of aftershocks indicate that the dynamic CFS changes are far greater than static ones. The maximum dynamic CFS change at the epicenter of Wuqia earthquake is 0.12 MPa, suggesting that the Wuqia event may be triggered by dynamic Coulomb stresses induced by the Yutian event, but the triggering effect is unremarkable and static CFS changes have little influence on it. -
引言
北京时间2008年3月21日,新疆于田县境内发生MS7.3地震,震中位于东北向的阿尔金断裂与东西向的普鲁断裂交汇区,地质构造较为复杂. 地震时和田、 喀什等地震感强烈,震中烈度达到Ⅶ度. 这次地震造成了部分建筑物损毁,但是没有人员伤亡. 地震后又发生了一系列余震,截止到2008年8月31日,新疆维吾尔自治区地震局共测定出于田地震序列MS≥4.0余震60次(李志海等,2009). 同年10月5日,其西北800 km的新疆乌恰县境内发生MS6.8地震,该地震微观震中位于乌恰县境内,与吉尔吉斯斯坦、 塔吉克斯坦接壤,该地震使乌恰县境内部分乡镇遭受到不同程度的地震灾害,而宏观震中位于吉尔吉斯斯坦共和国的努拉村,造成了努拉村70余人死亡(陈建波等,2009).
这两次地震事件之后一些研究人员从不同学科角度对其进行了研究(尹光华等,2008;陈建波等,2009;李志海等,2009;王琼等,2009;Wang et al, 2009; Elliott et al,2010;唐明帅等,2010;万永革等,2010;Shan et al,2011;王凡等,2011;张国宏等,2011;徐锡伟等,2011;唐兰兰等,2012). 其中王琼等(2009)计算了这两次地震产生的静态库仑应力,分析了两次事件对MS3.0以上余震以及天山中段、 柯坪块体和喀什-乌恰交汇区地震活动性的影响(Wang et al, 2009). 万永革等(2010)计算了此次地震对周围断层的影响,并讨论了该地震及其周围已发生的地震出现正断层机制的原因. 一次强震的发生所导致的应力变化对之后区域地震活动的影响被称为“应力触发”,一般分为静态应力触发和动态应力触发. 前者主要指发震断层错动在近场产生的永久性静态应力变化持续地对后续地震活动的触发作用;后者则指强震破裂传播或激发出的地震波通过某一区域时产生的较强即时应力变化对后续地震的触发作用. 另外,也有部分研究人员提出了黏弹性应力触发(Pollitz,Sacks,1997,2002;Freed,Lin,1998,2001),指震后下地壳和上地幔黏弹性介质在震后松弛过程中造成地壳应力场的调整,改变接收断层所处的应力环境,影响地震断层的地震活动性(汪建军,2010). 对于近场而言,静态和动态库仑应力变化的作用是分不开的(吴小平等,2007). Kilb等(2002)提出了完全库仑破裂应力概念,同时包含静态和动态库仑破裂应力变化.
Bouchon和Aki(1977)提出的离散波数法,能够精确地求解出完全格林函数,其特点是能够计算全波场,同时包括动态和静态的贡献. 本文基于离散波数法获取于田MS7.3地震的理论地震图,分析该地震对近场区域余震的完全库仑应力触发影响,并计算于田地震在乌恰地震震中处产生的动态库仑破裂应力变化.
1. 活动构造及动力学背景
2008年3月21日于田MS7.3地震震中位于青藏高原北缘ENE向阿尔金断裂西南端与近EW向康西瓦断裂及贡嘎错断裂的交汇部位(图1). 断错地貌测量显示阿尔金断裂西段全新世左旋走滑速率可达9—17.5 mm/a(徐锡伟等,2003). 康西瓦断裂的西端终止于塔什库尔干谷地东部,其东端在于田阿羌乡南部与阿尔金断裂相连,近EW向延伸达700 km,断错地貌显示晚第四纪以来,康西瓦断裂的平均左旋走滑速率为8—12 mm/a(付碧宏等,2006). 贡嘎错断裂在较大的阿尔金左旋走滑运动与较小的喀喇昆仑右旋走滑运动的共同作用下,呈左旋张扭性质. 野外地质考察显示,此次于田MS7.3地震走向为358°,近SN向的发震断层属于阿尔金断裂西南端NE向张性区构造(徐锡伟等,2003,2011; 王峰等,2004; Cowgill et al,2009;刁桂苓等,2010). 帕米尔高原东北侧位处天山褶皱带、 西昆仑造山带和塔里木块体等3大地质构造单元的交接地带. 印度板块与欧亚板块的碰撞和持续汇聚使其陆内俯冲和变形作用非常强烈,并在推挤前缘形成了一系列壮观的北向凸出的弧型推覆构造,这些大型推覆构造东部为正向逆冲,向西渐变为斜冲-走滑大断裂,断层显示出东部右旋、 西部左旋的活动特征. 其中位于北部的卡兹克阿尔特弧形活动褶皱-逆断裂带又由近EW向的卡巴加特弧和NW向的乌帕尔弧两个次级弧型构造构成,是6级以上地震频发的场所. 2008年10月5日的新疆乌恰MS6.8地震就发生在卡兹克阿尔特断裂的中西段(Molnar,1988;Burtman,Molnar,1993; 陈杰等,1997;唐兰兰等,2012).
图 1 于田MS7.3、 乌恰MS6.8地震主震震中位置及区域构造图 F1:焉耆盆地南缘断裂;F2:兴地断裂;F3:北轮台断裂;F4:柯坪断裂;F5:卡兹克阿尔特断裂;F6:米牙断裂;F7:西昆仑断裂;F8:康西瓦断裂;F9:喀喇昆仑断裂;F10:普鲁断裂;F11:贡嘎错断裂;F12:阿尔金断裂;F13:龙木错—邦达错断裂Figure 1. Epicenters of Yutian MS7.3 and Wuqia MS6.8 earthquakes all well as regional tectonic settings F1: Yanqi basin southern fringe fault; F2: Xingdi fault; F3: North Luntai fault;F4: Kalpin fault;F5: Stewart Katz Yasar fault;F6: Miya fault;F7: West Kunlun fault;F8: Kangxiwa fault;F9: Karakorum fault;F10: Pulu fault;F11: Gonggar Co fault;F12: Altyn Tagh fault;F13: Lungmu Co- -Bangdag Co fault2. 方法及原理
Bouchon和Aki(1977)发展了离散水平波数法,并将其应用于层状介质中震源的近场波场研究. 本文在此进行简单解释. 二维情况下,无限均匀介质中的一个线源稳态辐射可以用一个柱状波来表示,或者等价于一系列均匀和非均匀平面波的连续叠加. 令x和z分别表示水平和垂直坐标轴,则位移或应力可以表示为
式中,ω为角频率,k为水平波数. 对其进行积分时,如果介质有限或垂向不均匀,积分就会存在极点和奇异点,此时对水平波数的积分就会变得复杂和困难. 对于这个问题,一个简单的解决方法就是用一系列轴上周期分布的复合源来代替上述式(1)中的单一源. 设定复合源为沿轴周期性分布的源阵列,并且每两个相邻源之间都具有相同的距离间隔L和确定的时间延迟相,则式(1)变为
式中,L为周期源的距离间隔,对上述周期函数进行傅里叶展开严格求和,则式(2)化为
式中,kn=(2π/L)n. 如果其中的序列收敛,则可用有限求和形式近似表示为
由式(1)推到式(4),将一个单一源问题转化成了包含无限周期源的复合源问题,大大化简了运算. 下一步是将已经在频率域中解决的复合源问题再恢复至单一源问题. 但是实际操作时,仅能计算出某些频率的函数G,因此使用离散傅里叶反变换来求得时间域的解. 这样,一方面要处理一个含无限时间响应的单一源问题(由于假定了无限数目的周期源),另一方面要使用离散傅里叶反变换得到一个有限持续时间为T=2π/Δω的信号. 这些实际上可通过复频率域中作傅里叶反变换实现:
式中,ωI表示角频率的虚部常数,且满足
式(6)保证了先前无限时间响应解在时间窗T内的衰减. 如果选择的L足够大,以至于在时间窗T内没有最邻近的源造成对接收点(x,z)的干扰,那么时间域中单一源的解f(x,z; t)便可以通过恢复频率域中复合源的解G(x,z; ω)而得到,即
将上述二维问题延伸至三维情况下,与式(4)相对应有
式中,knx=(2π/Lx)nx,kny=(2π/Ly)ny.
与二维情况一样,对式(8)进行傅里叶反变换,把在频率域中已经解决的问题转换回时间域,从而得到地震破裂产生的位移场(Bouchon,Aki,1977;Bouchon,2003). 离散波数法可以计算全波场,可以得到完全库仑破裂应力变化,同时包含动态和静态库仑应力变化(Kilb et al,2002).
3. 资料选取及模型建立
3.1 资料选取
本文选取哈佛大学给出的于田MS7.3地震的基本参数,并采用冀战波等(2014)反演得到的于田地震震源破裂模型进行动态应力触发作用的研究. 结合乌恰地震震后考察及发震构造性质(谭明等,2008;陈建波等,2009),考虑到体波资料能更多地反映地震的初始应力状况,本文选用P波初动解作为乌恰地震的震源机制解. 于田MS7.3地震的基本参数如表1所示,地震震源模型如表2所示.
表 1 于田MS7.3地震的基本参数Table 1. Basic parameters of Yutian MS7.3 earthquake表 2 于田MS7.3地震的震源模型Table 2. Source model of Yutian MS7.3 earthquake利用地震台阵可以较为准确地估算出于田MS7.3地震信号的后方位角和慢度. 唐明帅等(2010)利用和田地震台阵数据,应用宽频f-k方法,对于田地震序列(2008年3月21日—8月31日,共119次ML≥4.0余震)进行了重新定位. 本文基于其余震精定位结果分析于田地震对近场区域余震的完全库仑应力触发影响.
3.2 计算范围及震源模型
对于近场和远场应力触发,本文选用了不同的计算范围和震源模型. 在近场为了获得分辨率更高的模型,选取1.2°×1°(35°—36.2°N,81°—82°E)为计算范围,以0.12°×0.1°为步长共获得接收点即场点121个;对于于田MS7.3地震对乌恰MS6.8地震的远场触发研究,本文选取5°×9°(35°—40°N,73°—82°E)为计算范围,以0.5°×0.9°为步长,获得场点121个. 根据冀战波等(2014)对于田地震破裂过程的反演研究,设震源为一个以有限速度扩展的双侧破裂正断层,断层上沿埋深0.5 km,下沿埋深30.5 km,破裂传播速度为3.0 km/s,视摩擦系数μ'取为0.5. 根据上节所述理论,将震源断层分解成10×5个小单元面,相当于50个位错点源,滑动振幅为1.5 m. 本文所选用的介质模型参考了胥颐等(2006)的地壳速度结构(表3).
表 3 介质模型参数Table 3. Layered structure model4. 计算结果及分析
4.1 于田MS7.3地震对近场区域余震的完全库仑应力作用
本文采用位错点源格林函数(Bouchon,1981),利用Bouchon(2003)提供的位移理论地震图计算程序,对上述模型计算其位移理论地震图; 进行弹性动力学转换后得到了于田MS7.3地震断层破裂在周围介质中产生的完全库仑破裂应力,并将余震序列的精定位结果投影到库仑破裂应力变化的空间分布图上. 图2给出了于田MS7.3地震近场完全库仑应力变化的演化图像. 可以看到,地震破裂开始后完全库仑应力图像发展很快,从40.00 s之后库仑应力变化图像即趋于稳定,在60.39 s以后图像已基本稳定. 此时由破裂产生的地震波已超出近场研究范围,剩下的仅是由永久性位移造成的静态库仑应力变化. 大部分余震分布在库仑破裂应力为正的区域,仅少量分布在其边缘及负区. 因此,静态库仑应力变化对余震有较为可靠的触发作用.
根据震源时间函数(图3)可以看出,地震矩释放在13 s左右达到峰值,而库仑应力变化图像在发展到13.33 s时比60.39 s时有更多的余震落在了正区(几乎在90%以上). 因此,考虑到地震波的传播速度和时间,有理由怀疑这部分在13.33 s时落在正区、 而在60.39 s时落在了负区的余震是由动态库仑应力触发的. 尽管静态库仑应力对其产生了延迟触发作用,但似乎动态库仑应力起到了最终的决定性作用. 而对于大部分的余震,动态和静态库仑应力变化对其均起到了一定的触发作用,但究竟是哪种库仑应力变化起到主导作用值得进一步研究.
为了更清楚地看到动态库仑应力和静态库仑应力变化各自对余震的触发影响,本文计算提取了余震震中处的完全库仑应力变化. 在此以两次余震为例,一次是发生于2008年3月24日04时56分16秒的ML5.6地震(图2中余震1),震中为35.59°N、 81.31°E. 由于受本文模型精度的限制,选取了最近的接收点(35.60°N、 81.30°E)作为参考,最终的完全库仑应力变化时程如图4所示. 可以看出,最大完全库仑应力变化出现在14.90 s,达到2.117 MPa,远远高于动态触发阈值;静态库仑应力变化稳定于0.7 MPa左右,同样远远高于0.01 MPa的静态触发阈值. 因此,动态和静态库仑应力变化可能对本次余震均有较明显的触发作用.
另外一次余震是发生于2008年3月21日20时02分57秒的ML5.4地震(图2中余震2),震中为35.28°N、 81.31°E. 同样由于受本文模型精度的限制,选取了最近的接收点(35.24°N、 81.30°E)作为参考,最终的完全库仑应力变化时程如图5所示. 最终静态库仑应力变化稳定在-2.5 MPa左右,超过了有意义的静态触发阈值±0.01 MPa; 正的完全库仑应 力变化最大值出现在10.20 s,为0.425 7 MPa,也接近有意义的动态触发阈值. 这个接收点同时受到了静态库仑应力的延迟触发作用和动态库仑应力的促进触发作用,于是在主震发生14个小时后,此处发生了ML5.4余震. 同样的情况,在附近区域至少有4次ML≥4.0的震例.
在于田MS7.3地震发生后,66.6%的M≥5.0和24.0%的M4余震发生在主震发生的当日,35.7%的M≥5.0和40.3%的M4余震发生在主震发生的一周内,其它余震均发生在主震发生的一周以后. Kilb等(2002)研究了美国L and ers M7.3地震之后库仑应力的变化(ΔCFS和ΔCFS(t)峰值)与地震活动率改变的空间关系,认为应力变化能够大到足以改变断层及其周围介质的物理和化学性质,使其更加接近于破裂. 这种改变可能并不会导致断层立即破裂,而是推迟一段时间. 另外,定量的统计学实验显示,ΔCFS和ΔCFS(t)与地震活动率改变之间的对应关系要比地震活动率改变作为一种随机结果更好一些,大约高10%. 虽然这10%并不起眼,但是它具有重要的统计学意义. 而且,它已经是一个非常保守的统计学估计. 对于一个在4.5 km深度处的模型,静态库仑应力变化振幅阈值为0.001—0.5 MPa,而对于峰值动态库仑应力变化来说,其阈值要≥0.5 MPa;其最佳的触发阈值分别为0.1和4.0 MPa(Kilb et al, 2002).
4.2 2008年于田MS7.3地震对乌恰MS6.8地震的动态库仑应力作用
乌恰MS6.8地震发生在于田MS7.3地震的6个月之后,两次地震震中相隔800 km,这两次地震事件之间是否存在应力触发关系是本文关心的问题. 如果应力触发关系存在,那么是动态还是静态触发呢? 王琼等(2009)计算分析了于田地震事件在乌恰地震事件震中位置上产生的静态库仑应力变化,结果为0—1 kPa. 针对这个问题,本文采用离散波数法重新计算了两次地震的动态库仑应力变化,得到了于田地震在乌恰地震震中位置产生的动态库仑破裂应力时程变化(图6).
图6为乌恰MS6.8地震震中12 km深处的库仑应力变化. 于田MS7.3地震断层破裂在乌恰地震震中处造成的最大动态库仑应力变化为0.12 MPa,其正的最大峰值没有达到0.5 MPa的动态触发阈值,静态库仑应力变化在1 kPa以内. 该结果与王琼等(2009)结果基本一致,远远小于动态变化,也未达到0.01 MPa的静态触发阈值. 因此本文认为,乌恰MS6.8地震可能受到了于田MS7.3地震一定的动态应力触发作用,但不显著,而静态库仑应力则对其影响很小.
5. 讨论与结论
5.1 于田MS7.3地震对近场区域余震的触发作用
于田MS7.3地震在近场区域的完全库仑应力演化图像显示,库仑应力变化对余震的触发率达到90%以上. 整个图像演化大概持续了60 s. 从地震波的传播速度和在近场研究区的传播时间来看,与震源时间函数的反演结果有较好的一致性. 根据震源时间函数,地震矩释放在13 s左右达到峰值,而库仑应力变化图像在发展到13.33 s时比60.39 s时有更多的余震落在正区. 在13.33 s时,有1个余震序列落在负区,9个余震序列落在不确定区域,109个余震序列落在正区,约占全部余震序列的91%;在60.39 s时,有10个余震序列落在负区,10个余震序列落在不确定区域,99个余震序列落在正区,约占全部余震序列的83%.在这次地震中,动态库仑应力变化图像更好地解释了余震的分布.
余震震中处的完全库仑应力变化计算结果表明,动态库仑应力变化远远大于静态库仑应力变化. 2008年3月24日的ML5.6地震最大动态库仑应力变化和静态库仑应力变化均远远高于各自触发阈值. 因此,动态和静态库仑应力变化可能对其均有较明显的触发作用. 2008年3月21日的ML5.4余震的震中接收点同时受到了静态库仑应力的延迟触发作用和动态库仑应力的促进触发作用,在主震发生14个小时后,此处发生了ML5.4余震. 同样的情况,在附近区域至少有4个ML≥4.0的震例.
5.2 于田MS7.3地震对乌恰MS6.8地震的远程触发作用
于田地震断层破裂在乌恰地震震中造成的最大动态库仑应力变化为0.12 MPa,其正的最大峰值没有达到0.5 MPa的动态触发阈值,但是远远高于1 kPa的静态库仑应力变化. 本文认为,乌恰MS6.8地震可能受到了于田MS7.3地震一定的动态应力触发作用,但不显著,而静态库仑应力则对其影响很小.
在于田MS7.3地震中,动态库仑应力变化图像更好地解释了余震的分布,那么在余震的发生中,两种库仑应力变化各自扮演了什么样的角色? 哪种库仑应力变化起到了更为重要的作用? 希望通过更多的震例和进一步的研究来厘清此问题.
新疆维吾尔自治区地震局唐明帅为本文提供了余震精定位结果,审稿专家对本文提出了宝贵修改意见,在此深表感谢.
-
图 1 于田MS7.3、 乌恰MS6.8地震主震震中位置及区域构造图 F1:焉耆盆地南缘断裂;F2:兴地断裂;F3:北轮台断裂;F4:柯坪断裂;F5:卡兹克阿尔特断裂;F6:米牙断裂;F7:西昆仑断裂;F8:康西瓦断裂;F9:喀喇昆仑断裂;F10:普鲁断裂;F11:贡嘎错断裂;F12:阿尔金断裂;F13:龙木错—邦达错断裂
Figure 1. Epicenters of Yutian MS7.3 and Wuqia MS6.8 earthquakes all well as regional tectonic settings F1: Yanqi basin southern fringe fault; F2: Xingdi fault; F3: North Luntai fault;F4: Kalpin fault;F5: Stewart Katz Yasar fault;F6: Miya fault;F7: West Kunlun fault;F8: Kangxiwa fault;F9: Karakorum fault;F10: Pulu fault;F11: Gonggar Co fault;F12: Altyn Tagh fault;F13: Lungmu Co- -Bangdag Co fault
表 1 于田MS7.3地震的基本参数
Table 1 Basic parameters of Yutian MS7.3 earthquake
表 2 于田MS7.3地震的震源模型
Table 2 Source model of Yutian MS7.3 earthquake
表 3 介质模型参数
Table 3 Layered structure model
-
汪建军. 2010. 同震、 震后和震间应力触发[D]: 武汉: 武汉大学测绘学院: 3-4. Wang J J. 2010. Coseismic, Postseismic and Interseismic Stress Triggering[D]. Wuhan: School of Geodesy and Geomatics, Wuhan University: 3-4 (in Chinese).
-
期刊类型引用(5)
1. 刘容,陈强,杨莹辉,钟霞,苑雨. 2021年云南漾濞M_S6.4地震动态应力对后续余震活动的触发作用. 地震学报. 2023(01): 17-28 . 本站查看
2. 魏芸芸,王琼,陈昊,钱才,向元,张琳琳. 2015年12月7日塔吉克斯坦M_S7.4地震对远场小震活动的动态应力触发. 内陆地震. 2023(01): 1-11 . 百度学术
3. 蒋子琴,杨莹辉,陈强,徐倩,徐浪,黄小梅. 2017年伊朗M_W7.3地震的震后余滑分布及其对2018年M_W6.0地震的触发影响. 地震学报. 2022(03): 452-466 . 本站查看
4. 刘方斌,古丽孜娜提·依德热斯,曲均浩. 新疆地区区域应力场特征及强震间库仑应力变化探讨. 地震工程学报. 2021(05): 1024-1036 . 百度学术
5. 陈昊,王琼,苏金波,魏芸芸. 2015年尼泊尔M_S8.1地震在喀什—乌恰交汇区的动态触发活动. 地震工程学报. 2017(02): 253-261 . 百度学术
其他类型引用(1)