数值模拟静态应力扰动下的断层失稳:结果分析兼与Dieterich模型和Coulomb模型的对比

基于大量岩石力学实验,Dieterich和其他研究者(Dieterich, 1978; Ruina, 1983)首先提出了描述岩石摩擦过程的速率-状态摩擦定律(R-S摩擦定律).如今R-S摩擦定律已成为研究地震成核等地震演化机制的有效手段.Dieterich (1992, 1994)最早提出了描述受静态剪应力扰动后断层失稳时间提前或推后的余震触发机制的解析模型.现在Dieterich模型已经成为解释余震随时间衰减规律的Omori定律等地震观测现象的有力工具.与之相对应,广泛使用的Coulomb应力失稳模型也可以给出断层受到静态剪应力扰动后,断层失稳时间的提前和推后量.Dieterich模型和Coulomb应力失稳模型基于不同的物理方法,所以在进行地震危险性评估时,二者均有各自的局限性.本文利用R-S摩擦定律控制的1-D弹簧-滑块模型,模拟计算了理论地震循环以及在不同静态剪切应力扰动下,失稳时间的提前和推后量的变化情况,然后将计算得到的时间提前和推后量分别与Dieterich模型和Coulomb模型的相应计算结果进行了定量化对比和差异性分析,并给出了相应的解释.数值模拟结果显示,对于R-S摩擦定律在参数不变的条件下,断层模型失稳时间的提前和推后量的大小强烈依赖于静态剪应力扰动的大小和作用时间,而且绝对值相同的正、负向静态剪应力扰动造成的失稳时间的提前和推后量的变化情况并不完全一致.在震后松弛/滑移阶段和闭锁阶段,时间提前和推后量是常数,且随静态剪应力扰动绝对值的增大而增大,两者的比值接近于1.0,这与Dieterich模型和Coulomb模型的结果是一致的,相应的差值小于两模型结果的10%.而在自加速阶段,模拟计算结果则存在与Dieterich模型和Coulomb模型结果不同的特征.首先,在自加速阶段模拟计算结果均偏离Coulomb模型,而且时间提前和推后量的比值小于1.0,相异于Coulomb模型的论断.不过当受到正向静态剪应力扰动后,Dieterich模型的结果和模拟计算结果是一致的,最大相差量不超过Dieterich模型结果的7%,可接近0.对于负向静态剪应力扰动,当其绝对值较小时,Dieterich模型的结果很接近模拟计算结果,相差量不超过该结果的14%.但对于绝对值较大的情况,模拟计算结果远大于Dieterich模型的结果,最大可达Dieterich模型结果的35倍,这是由于负向静态剪应力扰动后使得Ω=δθ/dc>>1的条件不再成立,进而使得Dieterich模型不再成立.总的来说,与模拟计算相比Dieterich模型可以很好地描述1-D断层受扰后失稳时间提前和推后量的变化情况,并且可以体现出正、负静态剪切应力扰动后失稳时间提前量和推后量变化的差异性,而Coulomb模型则不能完整地给出受到静态剪应力扰动后断层失稳时间提前或推后的估计值.

Numerical simulation of fault instability due to an arbitrary static stress perturbation:A comparison with the Dieterich model and Coulomb failure model