Numerical investigation of heating effect on the earthquake faulting based on the Chester-Higgs model
-
摘要: 速率和状态相依赖的摩擦定律是本文采用的重要定律。结合Chester-Higgs摩擦模型和McKenzie-Brune摩擦生热模型,在一维弹簧-滑块-断层近似模型下,利用四阶变步长的Dormand-Prince算法,研究探讨了断层摩擦生热对断层演化的影响。结果表明:与忽略温度影响的情形相比,摩擦生热造成的温度上升可导致断层滑移时刻的略微提前,并伴随着摩擦系数和状态变量的下降,同时也使得断层的滑移量和应力降略有减小,而滑移速率有所增大;另外,在考虑温度影响时,有效正应力和临界滑移距离也会影响断层的演化过程,断层上的有效正应力越大,断层失稳时刻越提前,温度上升越明显;断层的临界滑移距离越大,断层失稳时刻则越迟,温度上升越显著,但当临界滑移距离超过5 cm时,具有不同临界滑移距离的断层,失稳时的温度则基本保持一致。
-
关键词:
- Chester-Higgs摩擦模型 /
- McKenzie-Brune摩擦生热模型 /
- 一维弹簧-滑块-段层模型 /
- 摩擦生热 /
- 有效正应力 /
- 临界滑移距离
Abstract: Rate- and state-dependent friction (RSF) law is an empirical law derived from labo-ratory experiments related to rock friction. RSF law has been used to quantitatively describe complex fault friction processes. Currently, it has emerged as the theoretical basis for the study of seismogenesis and earthquake faulting. With a combination of the Chester-Higgs friction model and the McKenzie-Brune frictional heat generation model, in this study we have investi-gated the effect of frictional heating process on the fault temporal evolution based on a spring-slider-fault system subjected to a rate- and state-dependent friction law. The system equations are solved efficiently by Dormand-Prince method with adaptive steps. The results show that, compared with the case in which the temperature effect is neglected (unheated fault), the rise of temperature caused by frictional heating can lead to a slight time advance of fault instability, accompanied by abrupt decreases of the friction coefficient and state variable, respectively. In the case when the temperature effect is taken into consideration (heated fault), the slip and stress drop on the fault are slightly smaller than that on the unheated fault, while the slip rate becomes larger. In addition, the effective normal stress and critical slip distance can also affect the fault temporal evolution. The greater the effective normal stress on the heated fault is, the earlier the fault instability occurs, accompanied with higher temperature rising. The larger the critical slip distance of the heated fault is, the later the fault instability occurs with a significant temperature increase. However, when the critical slip distance is larger than 5 cm, the peak temperatures are almost the same when the fault is unstable. -
引言
近年来,从强震之间的应力交互角度开展强震间相互作用的研究,已成为震源物理学的一个热点课题(Harris,1998;Perfettini et al,1999 ;Stein,1999;King,Cocco,2001;Steacy et al,2005 )。通过计算强震序列中前一次事件对后续事件的库仑应力改变量,可以定量地判断前者对后者的触发/延迟作用,帮助研究人员理解强震序列的孕育过程(Wan,Shen,2010;Wang et al,2014 ;Jia et al,2018 )。同时,大地震的发生将显著影响其周缘地区的应力场状态,并可能造成周围断层提前或推迟发生地震,影响区域地震危险性(Parsons et al,2008 ;Toda et al,2008 ;Hu et al,2017 )。此外,以统计学的方法研究区域地震活动性的时空演化规律已成为刻画、了解地震事件之间的因果关系及区域构造活动的有力手段(Zhuang et al,2005 ;Peng et al,2012 ;Jia et al,2014 ,2018)。自1970年我国国家地震台网建立以来,地震目录越来越完备,这为统计分析地震活动特征提供了数据基础。在速率-状态本构关系框架下(Dieterich,1979;Ruina,1983),强震引起的库仑应力变化会直接导致地震活动性的时空变化,从而将库仑应力变化和地震活动性变化从理论上衔接起来,为强震之间的触发机制以及未来区域地震危险性的评估提供了有效方法(Aron,Hardebeck,2009;Cocco et al,2010 ;Toda,Enescu,2011;Ishibe et al,2015 )。
青藏高原是我国最主要的地震活动区之一,尤其是青藏高原中部的巴颜喀拉块体周缘,近20年来发生了一系列强震事件,显示了其构造运动正处于活跃期(程佳等,2011;闻学泽等,2011;Jia et al,2012 ;徐锡伟等,2017)。自本世纪初开始,巴颜喀拉块体周缘相继发生了2001年可可西里MW7.8地震、2008年于田MW7.1地震、2008年汶川MW7.9地震、2010年玉树MW6.9地震、2013年芦山MW6.6地震、2014年于田MW6.9地震和2017年九寨沟MW6.5地震(图1)。这7次强震在时间和空间上的密集发生,造成了巨大的人员伤亡和财产损失,引起了学术界和公众的普遍关注。从动力学的角度看,这些强震的发震机制均与巴颜喀拉块体东南向的运动相吻合,可能是块体构造运动的产物(Jia et al,2012 ;徐锡伟等,2017)。另一方面,一些研究表明,强震间的触发/延迟作用也可能主导了某些强震的发生(Jia et al,2014 ,2018)。因此,深入地细致地研究该强震序列之间可能存在的触发/延迟机制,有助于理解该强震序列的孕震机理和巴颜喀拉块体及周缘构造区的强震危险性。
图 1 巴颜喀拉块体构造背景及7次强震事件(活动断层数据引自张培震等,2003)Figure 1. The tectonic background of the Bayan Har block and the seven major earthquakes(The active faults data is from Zhang et al,2003 )多项研究表明,部分强震之间存在应力交互现象,并且强震的发生会对周围其它发震断层产生触发或抑制的作用,影响余震丛集、强震孕育等地震活动(Stein et al,1992 ;Harris,1998;Stein,1999;King,Cocco,2001;Pollitz et al,2003 ;Steacy et al,2005 )。King等(1994)研究了强震引起的同震库仑应力改变,并运用于1992年MW7.3 Landers地震,成功地解释了Landers地震的余震空间分布,即大多数余震发生在库仑应力增加的区域,只有少量地震发生在库仑应力减小的区域。Freed和Lin (2001)首次使用黏弹性应力触发假说解释了1992年Landers MW7.3地震对1999年Hector Mine MW7.1地震的延迟触发作用,使得黏弹性应力触发机制在地震触发中的作用被广泛关注。1999年土耳其MW7.4 Izmit地震提高了其后续MW7.1 Duzce地震的震源破裂起始点的黏弹性应力,增加量达70% (Hearn et al,2002 )。上述研究表明各个时空尺度的地震触发已成为一种普遍存在的现象,可为地震孕震机理提供更多的信息。
另一方面,使用统计学的方法分析地震活动的时空特征,进而判断地震之间的触发关系,也成为地震触发关系研究的热点内容之一。其中,结合了大森-宇津(Omori-Utsu)定律和古登堡-里查德(Gutenberg-Richter)定律的传染型余震序列(epidemic type aftershock sequence,简写为ETAS模型)(Ogata,1988,1998;Zhuang et al,2002 ;Ogata,Zhuang,2006),因其提供了一种有效区分背景地震和丛集(触发)地震的方法,而被广泛应用于不同构造背景的地震活动性研究中(Zhuang et al,2005 ;蒋长胜,庄建仓,2010;Peng et al,2012 ;Jia et al,2014 ,2018)。Dieterich (1994)从理论上将地震活动性与库仑应力改变直接联系起来,证实了两者之间存在定量的正相关关系,从而为结合地震活动性分析和库仑应力改变来探究地震触发关系提供了理论依据。因此,同时考虑地震活动性变化和库仑应力改变,能够互相印证地甄别地震之间的触发关系,有助于深入理解区域构造活动和地震危险性。
巴颜喀拉块体周缘强震序列之间是否存在触发关系是科学界普遍关注的问题。尤其是关于2008年汶川地震与2013年芦山地震及2017年九寨沟地震之间,以及2008年和2014年两次于田地震之间的关系,科学界存在不同的看法。刁桂苓等(2010)研究发现,于田地区的地震活动性在汶川地震发生之后显著增强,进而推断2008年于田地震和汶川地震之间可能存在触发关系。而关于芦山地震是否是汶川地震的余震这一科学问题,学术界仍存在争论:既有学者认为是余震的(陈运泰等,2013;王为民等,2013),也有研究推断是构造事件的(杜方等,2013;Xu et al,2013 );同时,一些研究人员认为芦山地震的发生既有汶川地震的触发作用,也有构造应力加载的贡献(Jia et al,2014 ;Wang et al,2014 );Jia等(2014)定量分析了汶川地震对芦山地区背景地震活动性的改变,认为芦山地震是汶川地震余震的概率是62%;Wang等(2014)以构造应力和库仑应力改变的比例为判断依据,推断芦山地震是汶川地震余震的概率是85%。针对2017年九寨沟地震与汶川地震的关系,也存在不同的声音:一些学者认为前者是后者触发引起的(单斌等,2017;Hu et al,2017 ),另一些研究则表明后者延迟了前者的发震时间(Wang et al,2014 ;Jia et al,2018 )。关于两次于田地震之间的关系,多项研究均表明,2008年于田地震所引起的静态和黏弹性应力的增加提前了2014年于田地震的发生(Li et al,2015b ;Wang et al,2017 )。程佳等(2011)考虑黏弹性应力作用,计算了玛尼地震、可可西里地震和汶川地震所形成的同震和震后形变场的变化过程和特征,认为各地震之间的影响很微弱。以上针对巴颜喀拉块体周缘强震之间触发关系的研究,对理解该强震序列的发震机理是构造运动还是触发作用,或者两者皆有,均具有重要的科学意义。
本文拟对自本世纪初起发生在巴颜喀拉块体边界的7次强震之间的触发关系进行系统研究(表1)。以库仑应力改变量和地震活动性的改变为依据,甄别地震之间的交互作用。根据发震时间和震中距离的不同,本文对汶川地震−芦山地震−九寨沟地震、两次于田地震进行单独研究,将其余强震作为一个序列进行分析讨论,再根据这些地震之间的触发关系,推断巴颜喀拉块体现今的构造运动和未来的地震危险性。
表 1 巴颜喀拉块体周缘7次强震的震源机制解Table 1. The focal mechanism solutions of seven major earthquakes occurred around the Bayan Har block序号 发震时间
年-月-日东经/° 北纬/° MW 走向/° 倾角/° 滑动角/° 深度/km 1 2001−11−14 92.91 35.80 7.8 94 51 −12 15.0 2 2008−03−21 81.37 35.43 7.1 203 52 −74 12.0 3 2008−05−12 104.10 31.44 7.9 231 35 138 12.8 4 2010−04−14 96.79 33.05 6.9 300 88 23 15.7 5 2013−04−20 103.12 30.22 6.6 212 42 100 21.9 6 2014−02−12 82.57 36.22 6.9 242 86 −5 18.3 7 2017−08−08 103.90 33.21 6.5 150 78 −13 14.9 注:数据参考Ekström等 (2012)。 1. 库仑应力改变
由于各强震间发震时刻的间隔(数年)远大于动态应力触发的特征时间(数小时至数周),因此本文不考虑该强震序列的动态应力改变,仅关注同震时的静态应力改变和震后的黏弹性应力改变。本研究根据巴颜喀拉块体深部结构的研究结果,将地球模型简化为包含弹性上地壳、黏弹性下地壳和黏弹性上地幔的层状模型,其中下地壳和上地幔的粘性系数分别为1×1018 Pa·s和1×1020 Pa·s,模型具体细节请参考Xiong等(2010)和Jia等(2012)。由于汶川地震与芦山地震、九寨沟地震以及两次于田地震之间的空间距离较近,且相关研究较多(Wang et al,2014 ;赵立波等,2016;单斌等,2017)。为了方便对比研究,分别将汶川地震对芦山地震和九寨沟地震的影响,芦山地震对九寨沟地震的影响,以及2008年于田地震对2014年于田地震的影响进行单独分析。同时定量地计算可可西里地震、2008年于田地震、汶川地震和玉树地震之间的库仑应力改变量。
1.1 库仑应力改变的计算方法
库仑应力改变的计算公式为(Robinson,2005;Jia et al,2014 )
${\rm d}{\rm{CFS}} = {\rm d}{\tau _{\rm{s}}} + \mu \left( {{\rm d}{\tau _{\rm{n}}} - \left( {\frac{\beta }{3}} \right)\sum {{\tau _{ii}}} } \right),$
(1) 式中,dτs为剪切应力改变量,μ为摩擦系数,dτn为正应力改变量(拉张为正),β为Skempton系数。相比于部分研究使用视摩擦系数综合孔隙压的变化,我们更倾向于使用Skempton系数来更准确地刻画孔隙压力的变化,因为前者可能导致库仑应力计算结果不准确(Beeler et al,2000 )。Skempton系数根据实验室和野外观测的结果取0.5,而摩擦系数因不同发震断层所处的环境不同而取不同数值,取值范围介于0—1.0之间 (Robinson,2005;Jia et al,2014 )。使用PSGRN/PSCMP程序(Wang et al,2006 ),给定地震滑动分布模型,接受断层的发震机制,即可根据黏弹性位错模型计算出强震引起的静态及黏弹性应力改变量。
1.2 可可西里地震、2008年于田地震、汶川地震和玉树地震之间的库仑应力改变
可可西里地震、2008年于田地震和汶川地震的破裂模型分别参考Lasserre等(2005)、张国宏等(2011)以及Ji和Hayes (2008)的结果。摩擦系数和Skempton系数根据实验室和野外观察的结果,分别选取0.75和0.5 (Robinson,2005)。库仑应力的计算结果列于表2。其中由于2008年于田地震和汶川地震的震中相距非常远(2 500 km),而同震应力改变与两次地震的震中距离呈二次反比的关系,因此这两次地震之间的同震应力改变可以忽略;同时这两次地震的发震时间间隔较短(52天),远小于震后黏弹性应力释放发挥作用的特征时间,因此震后黏弹性应力改变也可以忽略,所以表2中并未列出这两次地震之间的库仑应力改变量。从表2可见,可可西里地震对后续2008年于田、汶川和玉树地震的库仑应力改变非常小,远小于可能引起地震活动性变化的阈值0.01 MPa (Reasenberg,Simpson,1992;Stein et al,1992 )。类似的,2008年于田地震和汶川地震对玉树地震的库仑应力改变量也非常小。由上述结果可以推断,这4次地震事件之间并无触发作用。
表 2 巴颜喀拉强震序列在其后续地震发震位置所引起的库仑应力变化Table 2. The dCFS results induced by the strong earthquake sequence occurred around the Bayan Har block on following earthquakes2001年可可西里地震 2010年玉树地震 同震应力改变
/MPa震后应力改变
/MPa同震应力改变
/MPa震后应力改变
/MPa于田地震 7.00×10−5 1.30×10−4 7.28×10−7 8.00×10−7 汶川地震 1.14×10−5 1.20×10−5 −1.32×10−4 −2.50×10−4 玉树地震 −7.15×10−4 −1.35×10−3 注:震后应力改变指同震和震后黏弹性松弛共同导致的结果。 1.3 汶川地震、芦山地震和九寨沟地震之间的库仑应力改变
为了讨论汶川地震对芦山地震和九寨沟地震的影响,以及芦山地震对九寨沟地震的影响,本文使用不同的滑动分布模型和摩擦系数计算强震之间的库仑应力改变量,计算结果如图2所示。汶川地震的滑动模型使用Ji和Hayes (2008)以及Wang等(2011)的反演结果,芦山地震使用Jiang等(2014)和Zhang等(2014)的反演结果。由图2a可知,分别使用Wang等(2011)、Ji和Hayes (2008)的破裂模型,汶川地震对芦山地震的同震库仑应力改变分别达到了0.006—0.009 MPa和0.010—0.012 MPa,而总体的库仑应力改变(同震与震后库仑应力改变之和)分别为0.031—0.034 MPa和0.041—0.049 MPa。显然,这种程度的库仑应力改变量表明汶川地震对芦山地震具有一定的触发作用,也与其它研究结果相似(Parsons et al,2008 ;Toda et al,2008 ;单斌等,2009;Wan,Shen,2010;Jia et al,2014 ;Wang et al,2014 )。由图2b可知,除了使用Wang等(2011)的模型在μ=0.8的条件下,库仑应力改变为正值外,在μ=0.4以及使用Ji和Hayes (2008)的模型时,库仑应力改变量均为负值。这是由库仑应力改变的计算结果对摩擦系数有较高的敏感度造成的,也能解释为何针对汶川地震对九寨沟地震的库仑应力改变,不同研究可能得到相反的结果(Parsons et al,2008 ;Toda et al,2008 ;Luo,Liu,2010;Wan,Shen,2010;Nalbant,McCloskey,2011;Wang et al,2014 ;Hu et al,2017 )。通常,断层面上摩擦系数的精确值很难获得,因而仅假设其取值范围介于0—1.0之间。Li等(2015a)通过汶川地震断裂带科学钻探项目直接观测到断层面上的温度数据,推断龙门山断裂带处于低摩擦环境(μ<0.2)。因此,使用低摩擦系数计算得到的负库仑应力改变量更为可信,表明汶川地震可能对九寨沟地震具有一定的延迟作用。由图2c可知,芦山地震对九寨沟地震的库仑应力改变量非常微弱,约为10−6 MPa,可以推断芦山地震对九寨沟地震并无影响。
图 2 强震之间的库仑应力改变随时间的变化(a) 汶川地震对芦山地震的库仑应力改变;(b) 汶川地震对九寨沟地震的库仑应力改变;(c) 芦山地震对九寨沟地震的库仑应力改变Figure 2. Coulomb stress changes dCFS of strong earthquakes versus time(a) The dCFS induced by the Wenchuan earthquake on the Lushan earthquake;(b) The dCFS induced by the Wenchuan earthquake on the Jiuzhaigou earthquake;(c) The dCFS induced by the Lushan earthquake on the Jiuzhaigou earthquake1.4 两次于田地震之间的库仑应力改变
使用同样的方法计算两次于田地震之间的库仑应力改变量(图3),解析深度为10 km。由图3可知,2008年于田地震对2014年于田地震的库仑应力改变量从同震时的约−0.001 MPa至2014年于田地震前的约−0.003 MPa,虽然是负值,但其强度较小。该结果与前人的研究结果并不完全一致(Li et al,2015b ;王辉等,2016),但各项计算结果均显示库仑应力改变量非常小(赵立波等,2016;Wang et al,2017 ),表明这两次地震事件之间触发关系很微弱。
2. 地震活动性改变
2.1 ETAS模型和随机除丛法
ETAS模型结合了古登堡-里查德定律、大森-宇津定律及地震发生率与震级的关系等统计学规律,将某一地区某一时段内大于某震级的地震活动性描述为点过程,并将整体地震活动性看作不随时间变化的背景活动性和余震活动性之和。最重要的是,ETAS模型中的每次地震(包括余震)均有能力触发其余震,这种触发能力与其震级大小有关。ETAS模型以概率的形式来定义每次地震,即每次地震均有一定概率是其之前地震的余震,余下的概率作为背景地震的概率。然而在现实中无法实际观测到每次地震与其余震之间的关系,而只能记录到地震的时间序列。但ETAS模型能够通过对某一地区的地震目录建模,拟合其模型参数,从而分辨出各地震的余震概率和背景概率。
ETAS模型的数学含义为:给定时刻 t之前的地震目录,在单位时间内,位置(x,y)处震级大于完备震级Mc的地震数目,其表达式为(Ogata,1998;Zhuang et al,2005 )
$\lambda \left( {t,x,y} \right) = \mu \left( {x,y} \right) + \sum\limits_{\left( {i:{t_i} < t} \right)} {\kappa \left( {{M_i}} \right)g\left( {t - {t_i}} \right)f\left( {x - {x_i},y - {y_i};{M_i}} \right)} ,$
(2) $\kappa \left( M_i \right) = A{{\rm{e}}^{\alpha \left( {M_i - {M_{\rm{c}}}} \right)}},\quad M_i {\text{≥}} {M_{\rm{c}}},$
(3) $g\left( t \right) = \frac{{p - 1}}{c}{\left( {1 + \frac{t}{c}} \right)^{ - p}},$
(4) $f\left( {x,y;M_i} \right) = \frac{{q - 1}}{{{\rm{\pi }}D{{\rm{e}}^{\gamma \left( {M_i - {M_c}} \right)}}}}{\left[ {1 + \frac{{{x^2} + {y^2}}}{{D{{\rm{e}}^{\gamma \left( {M_i - {M_{\rm{c}}}} \right)}}}}} \right]^{ - q}},$
(5) 式中:μ(x,y)为背景地震活动率,假设其不随时间变化;κ(Mi)为地震发生率,即震级为Mi的主震能够触发的余震数目;g(t)为大森-宇津定律的概率密度函数;f(x,y;Mi)为空间概率密度函数;A,α,p,c,q,D,γ均为模型参数,表示给定区域地震活动性的统计特征。给定完备的地震目录和合适的模型参数初始值,通过最大似然法即可估计出各模型参数的值。
2.2 ETAS模型拟合
将整个巴颜喀拉块体看作一个整体进行ETAS模型拟合,空间范围为(29°N—37°N,79°E—105°E),时间范围为1980年1月1日至2017年8月27日。前人的研究表明,强震发生之后存在小地震无法被记录到的现象(Enescu et al,2007 ;Peng et al,2012 ),导致地震目录不完备,进而引起ETAS模型拟合参数的不稳定(Utsu et al,1995 ;Ogata,1998)。这一问题可以通过选用高于完备震级的地震目录来解决。由于观测到汶川地震之后大量小震(M≤4.0)缺失,因此本研究选取的完备震级为4.0 (Jia et al,2014 ,2018)。利用最大似然法估计ETAS模型参数,最终结果为μ=0.754,A=0.474,c=0.025,α=0.839,p=1.134,D=2.848×10−2,q=1.765,γ=0.46。
2.3 背景地震活动性变化
得到ETAS模型参数之后,可以计算地震事件j处背景地震活动率与整体地震活动率的比值,即地震事件j的背景概率,其表达式为
${\varphi _j} = \frac{{\mu \left( {{x_j},{y_j}} \right)}}{{\lambda \left( {{t_j},{x_j},{y_j}} \right)}}.$
(6) 为了描述背景地震活动率随时间的变化,Zhuang等(2005)定义了累积背景地震概率S(t)为
$S(t) = \sum\limits_{{t_i} < t} {{\varphi _i}} .$
(7) 对于构造运动稳定的地区,背景地震活动性不随时间变化,因此S(t)曲线应为随时间稳定上升的直线。当S(t)曲线发生突然变化时,预示着该区域背景地震活动性的突变。
根据活动构造与地震活动性空间分布的特征,在7次强震周围分别选取6个研究区(图4),展示其各自的地震活动性随时间的变化规律。由于2008年与2014年两次于田地震的震中位置较近,故将这两次地震的发震区域作为同一个研究区来分析地震活动性的变化。
图5给出了各个研究区的总体地震活动性和背景地震概率累积曲线,可以根据后者随时间的变化推断各地区构造活动的变化,进而推断该研究区被强震影响的可能性。从图5a可以看出,2001年可可西里地震之后,该地区背景地震活动性显著升高,表明该地区构造运动活跃;从图5b可以看出,于田地区背景地震活动性在每次强震发生时均无明显变化,表明其它地震对两次于田地震并无明显影响,两次于田地震之间也无显著触发/延迟关系。同时也可以看出,两次于田地震也并未提升该地区的背景地震活动性;类似的,汶川地区背景地震活动性在可可西里地震时并未出现显著变化,而2008年汶川地震显著提高了该地区的背景地震活动性(图5c);玉树地区在可可西里地震、2008年于田地震和汶川地震前后背景地震活动性也均未发生显著变化(图5d),指示了玉树地区并未受这3次地震事件的影响;而芦山地区的背景地震活动性也仅仅在2008年汶川地震之后有明显升高,表明汶川地震对芦山地震有触发作用(图5e);九寨沟地区在汶川地震后背景地震活动性出现明显降低,而在其它地震发生时则无显著变化(图5f),表明汶川地震对九寨沟地区地震活动性有压制作用。
综合库仑应力改变和地震活动性变化的结果可以看出,各强震的库仑应力变化结果与其造成的背景地震活动性变化相一致。在该强震序列中,汶川地震对芦山地震和九寨沟地震均有显著的库仑应力改变以及背景地震活动性变化,而其它地震之间均无明显的触发/延迟作用以及显著的背景地震活动性变化。
3. 讨论与结论
强震之间的应力交互是一个复杂的物理过程。由于计算参数的不确定性,例如缺乏约束的摩擦系数或滑移分布模型的不一致性等,往往很难准确地计算出强震引起的库仑应力改变。此外,观测到的地震活动性不仅包含来源于构造运动的背景地震事件,还包含了强震触发引起的触发地震。利用ETAS模型能够有效地区分出背景地震和触发地震,从而为分析背景地震活动性的时空变化提供基础。库仑应力改变与地震活动性改变呈正相关,通过速率-状态摩擦率有机地联系在一起,为研究强震之间的触发机制提供有力的证据。
总体来说,本文针对巴颜喀拉块体周缘的强震序列,从库仑应力改变和地震活动性这两个角度进行探讨分析,试图辨别这一备受关注的强震序列是构造运动的产物还是触发作用的结果,并进一步讨论了巴颜喀拉块体未来的地震危险性。库仑应力改变和地震活动性分析的结果均表明,汶川地震对芦山地震的触发作用和对九寨沟地震的延迟作用较为明显,展现出局部地区也存在强震之间互相影响的特征。除此之外,其它地震之间均无明显的触发或延迟作用。这表明该强震序列的孕震机制主要来源于巴颜喀拉块体不断向东南挤压的构造运动。因此,可以推断巴颜喀拉块体目前的构造运动较为活跃,该区域未来的地震危险性不可忽视。同时,强震触发/延迟周围活动断层的破裂也是值得关注的现象。
中国地震局地球物理研究所蒋长胜研究员提供了整理和修正后的地震目录,德国波茨坦地学研究中心汪荣江研究员提供了PSGRN/PSCMP的程序,日本数理统计研究所庄建仓教授提供了ETAS模型和随机除丛法的程序,两位匿名评审专家为本文的撰写提出了宝贵的意见和建议,作者在此一并表示衷心的感谢。
-
图 3 小时间尺度下失稳前后断层模型的演化
(a) 滑移速率
${\dot \delta}$ 演化图;(b) 摩擦系数μ演化图;(c) 状态变量Θ演化图;(d) 温度T演化图Figure 3. Simulated time histories of fault evolution from small time scale
The evolution of system around the onset of instability is shown. Figs. (a),(b),(c) and (d) show the evolutions of slip rate
${\dot \delta}$ ,frictional coefficient μ,state variable Θ and temperature change T,respectively图 4
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
(a) 摩擦系数-滑动速率相图;(b) 摩擦系数-位移相图;(c) 位移-滑动速率相图
Figure 4.
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
图 6 不同有效正应力对应的断层模型的演化(小时间尺度)
(a) 滑移速率
${\dot \delta}$ 演化图;(b) 摩擦系数μ演化图;(c) 状态变量Θ演化图;(d) 温度T演化图Figure 6. Simulated time histories of fault evolution corresponding to different effective normal stresses (small time scale)
Figs. (a),(b),(c) and (d) show the evolutions of slip rate
${\dot \delta}$ ,frictional coefficient μ,state variable Θ and temperature change T,respectively图 8 不同临界滑移距离所对应的断层模型的演化(大时间尺度)
(a) 滑移速率
${\dot \delta}$ 演化图;(b) 摩擦系数μ演化图;(c) 状态变量Θ演化图;(d) 温度T演化图Figure 8. Simulated time histories of fault evolution corresponding to different critical slip distances (large time scale)
Figs. (a),(b),(c) and (d) show the evolutions of slip rate
${\dot \delta}$ ,frictional coefficient μ,state variable Θ and temperature change T,respectively表 1 一维弹簧-滑块-断层模型参数的定义和取值
Table 1 Definitions and values of parameters for 1D spring-slider-fault model
参数 含义 取值 参数 含义 取值 a 直接影响系数 0.012 ${\dot \delta }$ u失稳速率 1 m/s b 演化影响系数 0.017 T* 参考温度 550 K σ 有效正应力 100—1 000 MPa Tini 初始温度 550 K Dc 临界滑移距离 0.01—10 cm Qa 直接影响的表面激活能 105 J/mol k/kc 断层刚度与临界刚度的比值 ≈1 Qb 演化影响的表面激活能 105 J/mol μ0 T=T*时以 ${\dot \delta} $ =${\dot \delta }$ *稳定滑动的摩擦系数0.6 R 气体常数 8.314 μini 初始摩擦系数 0.623 P 介质密度 2 600 kg/m3 ${\dot \delta} $ 0远场加载速率 3.5 cm/a c 介质比热 1 000 J/(kg·K) ${\dot \delta }$ *参考速率 3.5 cm/a κ 介质固体热扩散系数 10−6 m2/s ${\dot \delta }$ ini初始速率 0.035 cm/a -
陈运泰. 2003. 地震参数: 数字地震学在地震预测中的应用[M]. 北京: 地震出版社: 7–11. Chen Y T. 2003. Earthquake Parameters: Application of Digital Seismology in Earthquake Prediction[M]. Beijing: Seismological Press: 7–11 (in Chinese).
李世愚, 和泰名, 尹祥础. 2015. 岩石断裂力学[M]. 北京: 科学出版社: 1−208. Li S Y, He T M, Yin X C. 2015. Fracture Mechanics of Rock[M]. Beijing: Science Press: 1−208.
姚路,马胜利. 2013. 断层同震滑动的实验模拟:岩石高速摩擦实验的意义、方法与研究进展[J]. 地球物理学进展,28(2):607–623 Yao L,Ma S L. 2013. Experimental simulation of coseismic fault sliding:Significance,technological methods and research progress of high-velocity frictional experiments[J]. Progress in Geophysics,28(2):607–623 (in Chinese)
姚路. 2014. 龙门山断裂带断层泥中速-高速摩擦性质的实验研究[J]. 国际地震动态,8:41–42 doi: 10.3969/j.issn.0235-4975.2014.02.010 Yao L. 2014. Experimental study on mid-high speed tribological properties of fault mud in Longmenshan fault zone[J]. Recent Development in World Seismology,8:41–42 (in Chinese)
姚路, 马胜利, 嶋本, 利彦. 2016. 利用岩石高速摩擦实验认识地震断层滑动的物理化学过程: 现状与展望[C]//2016中国地球科学联合学术年会论文集(二十五)——专题48: 地震震源物理研究进展、专题49: 利用人工震源探测地下介质结构及其变化. 北京: 中国地球物理学会, 中国地震学会, 全国岩石学与地球动力学研讨会组委会, 中国地质学会构造地质学与地球动力学专业委员会, 中国地质学会区域地质与成矿专业委员会: 1357. Yao L, Ma S L, Dao B, Li Y. 2016. Recognition of physicochemical process of seismic fault slip by high speed friction experiment of rocks: Status and prospect[C]// The Paper Collection of 2016 Annual Meeting of Chinese Geoscience Union (Twenty-fifth)-Topic 48: Advances in Seismic Source Physics, Topic 49: Detection of Underground Medium Structure and Change by Artificial Source. Beijing: Chinese Geophysical Society, Seismological Society of China, Organizing Committee of National Symposium on Petrology and Geodynamics, Professional Committee of Tectonic Geology and Geodynamics of the Chinese Geological Society, Regional Geology and Mineralization Committee of the Chinese Geological Society: 1357.
曾融生. 1984. 固体地球物理学导论[M]. 北京: 科学出版社: 351–363. Zeng R S. 1984. Introduction of Solid Geophysics[M]. Beijing: Science Press: 351–363 (in Chinese).
Ampuero J P,Rubin A M. 2008. Earthquake nucleation on rate and state faults:Aging and slip laws[J]. J Geophys Res,113(B1):B01302
Barbot S,Lapusta N,Avouac J P. 2012. Under the hood of the earthquake machine:Toward predictive modeling of the seismic cycle[J]. Science,336(6082):707–710 doi: 10.1126/science.1218796
Beeler N M,Lockner D A. 2003. Why earthquakes correlate weakly with the solid Earth tides:Effects of periodic stress on the rate and probability of earthquake occurrence[J].J Geophys Res,108(B8):2391 doi: 10.1029/2001JB001518
Bhattacharya P,Rubin A M. 2014. Frictional response to velocity steps and 1-D fault nucleation under a state evolution law with stressing-rate dependence[J]. J Geophys Res,119(3):2272–2304 doi: 10.1002/2013JB010671
Bizzarri A,Cocco M,Andrews D J,Boschi E. 2001. Solving the dynamic rupture problem with different numerical approaches and constitutive laws[J]. Geophys J Int,144(3):656–678 doi: 10.1046/j.1365-246x.2001.01363.x
Bizzarri A. 2010. Determination of the Temperature Field due to Frictional Heating on a Sliding Interface[R]. Istituto Nazionale di: 1–16.
Bizzarri A. 2011. Temperature variations of constitutive parameters can significantly affect the fault dynamics[J]. Earth Planet Sci Lett,306(3/4):272–278
Blanpied M L,Tullis T E,Weeks J D. 1998. Effects of slip,slip rate,and shear heating on the friction of granite[J]. J Geophys Res,103(B1):489–511 doi: 10.1029/97JB02480
Byerlee J D. 1970. The mechanics of stick-slip[J]. Tectonophysics,9(5):475–486 doi: 10.1016/0040-1951(70)90059-4
Byerlee J D. 1978. Friction of rocks[J]. Pure Appl Geophys,116(4):615–626
Chester F M,Higgs N G. 1992. Multimechanism friction constitutive model for ultrafine quartz gouge at hypocentral conditions[J]. J Geophys Res,97(B2):1859–1870 doi: 10.1029/91JB02349
Chester F M. 1994. Effects of temperature on friction:Constitutive equations and experiments with quartz gouge[J]. J Geophys Res,99(B4):7247–7261 doi: 10.1029/93JB03110
Dieterich J H. 1979. Modeling of rock friction:1. Experimental results and constitutive equations[J]. J Geophys Res,84(B5):2161–2168 doi: 10.1029/JB084iB05p02161
Dieterich J H. 1992. Earthquake nucleation on faults with rate- and state-dependent strength[J]. Tectonophysics,211(1/4):115–134
Dieterich J H. 1994. A constitutive law for rate of earthquake production and its application to earthquake clustering[J]. J Geophys Res,99(B2):2601–2618 doi: 10.1029/93JB02581
Gu J C,Rice J R,Ruina A L,Tse S T. 1984. Slip motion and stability of a single degree of freedom elastic system with rate and state dependent friction[J]. J Mech Phys Solids,32(3):167–196 doi: 10.1016/0022-5096(84)90007-3
Gu Y J,Wong T F. 1991. Effects of loading velocity,stiffness,and inertia on the dynamics of a single degree of freedom spring-slider system[J]. J Geophys Res,96(B13):21677–21691 doi: 10.1029/91JB02271
Hatano T. 2015. Rate and state friction law as derived from atomistic processes at asperities[J/OL]. Geophys J Int [2017−09−21]. https://arxiv.org/pdf/1512.05078.pdf.
He C R,Wong T F,Beeler N M. 2003. Scaling of stress drop with recurrence interval and loading velocity for laboratory-derived fault strength relations[J]. J Geophys Res,108(B1):2037
He C R,Luo L,Hao Q M,Zhou Y S. 2013. Velocity-weakening behavior of plagioclase and pyroxene gouges and stabilizing effect of small amounts of quartz under hydrothermal conditions[J]. J Geophys Res,118(7):3408–3430 doi: 10.1002/jgrb.50280
Kame N,Fujita S,Nakatani M,Kusakabe T. 2013. Effects of a revised rate- and state-dependent friction law on aftershock triggering model[J]. Tectonophysics,600:187–195 doi: 10.1016/j.tecto.2012.11.028
Kaneko Y,Avouac J P,Lapusta N. 2010. Towards inferring earthquake patterns from geodetic observations of interseismic coupling[J]. Nat Geosci,3(5):363–369 doi: 10.1038/ngeo843
Kato N. 2001. Effect of frictional heating on pre-seismic sliding:A numerical simulation using a rate-,state- and temperature-dependent friction law[J]. Geophys J Int,147(1):183–188 doi: 10.1046/j.0956-540x.2001.01531.x
Kimura T. 2009. On dormand-prince method[EB/OL]. [2009−09−24]. http://depa.fquim.unam.mx/amyd/archivero/DormandPrince_19856.pdf.
King G C P,Cocco M. 2001. Fault interaction by elastic stress changes:New clues from earthquake sequences[J]. Adv Geophys,44:1–38 doi: 10.1016/S0065-2687(00)80006-0
Lu Z,He C R. 2014. Frictional behavior of simulated biotite fault gouge under hydrothermal conditions[J]. Tectonophysics,622:62–80 doi: 10.1016/j.tecto.2014.03.002
Marone C. 1998. Laboratory-derived friction laws and their application to seismic faulting[J]. Annu Rev Earth Plant Sci,26:643–696 doi: 10.1146/annurev.earth.26.1.643
McKenzie D,Brune J N. 1972. Melting on fault planes during large earthquakes[J]. Geophys J Int,29(1):65–78 doi: 10.1111/j.1365-246X.1972.tb06152.x
Perfettini H,Avouac J P. 2004a. Postseismic relaxation driven by brittle creep:A possible mechanism to reconcile geodetic mea-surements and the decay rate of aftershocks,application to the Chi-Chi earthquake,Taiwan[J]. J Geophys Res,109(B2):B02304
Perfettini H,Avouac J P. 2004b. Stress transfer and strain rate variations during the seismic cycle[J]. J Geophys Res,109(B6):B06402
Raleigh C B,Healy J H,Bredehoeft J D. 1976. An experiment in earthquake control at Rangely,Colorado[J]. Science,191(4233):1230–1237 doi: 10.1126/science.191.4233.1230
Rice J R,Tse S T. 1986. Dynamic motion of a single degree of freedom system following a rate and state dependent friction law[J]. J Geophys Res,91(B1):521–530 doi: 10.1029/JB091iB01p00521
Rice J R. 2006. Heating and weakening of faults during earthquake slip[J]. J Geophys Res,111:B05311
Ruina A. 1983. Slip instability and state variable friction laws[J]. J Geophys Res,88(B12):10359–10370 doi: 10.1029/JB088iB12p10359
Scholz C H. 1998. Earthquakes and friction laws[J]. Nature,391(6662):37–42 doi: 10.1038/34097
Scholz C H. 2002. The Mechanics of Earthquakes and Faulting[M]. 2nd ed. New York: Cambridge University Press: 351–412.
Segall P. 2010. Earthquake and Volcano Deformation[M]. Princeton: Princeton University Press: 332–369.
Stein R S. 1999. The role of stress transfer in earthquake occurrence[J]. Nature,402(6762):605–609 doi: 10.1038/45144
Stein S, Wysession M. 2003. An Introduction to Seismology, Earthquakes, and Earth Structure[M]. Malden: Blackwell Publishing: 215–217.
-
期刊类型引用(1)
1. 黄长生,侯保全,易秤云,李龙,张胜男,周耘,Waseem Akram,王芳婷. 赣南地热水形成条件分析与赣县区地热水勘查靶区圈定. 华南地质. 2021(01): 64-74 . 百度学术
其他类型引用(2)