具有脆性形变和蠕动变形的复合单断层系模型中的应力、滑动与地震

摘要

对弹性半空间中的几个含有不连续断层系的三维弹性连续介质模型,在沿单元垂直走滑断层发生的滑动过程进行了数值模拟。模型的几何形状和所使用的边界条件大致与圣安德烈斯断层中部的情形相符合。数值模拟综合考虑了脆性形变和蠕动变形。断层带的总形变率为蠕变率和摩擦滑动率之和。脆性断层的性质由破裂段(数值单元)的不同地震应力降的分布来给定。所假设的分布描述了两种理想化的情形,它们分别对应于不同断层滑动的极限状态:(1)以大尺度范围为特征的强无序状态,表示不成熟的断层带及扩展空间域。(2)相对有序的状态,以小尺度范围为特征,表示成熟的完全滑动断层。假定断层的蠕变特性在所有的情形中都相同,并且是用蠕滑速率与应力的幂次律关系来刻画的。利用这种脆性蠕变过程和模型参数可以得到一个在12.5km 深度含有“脆性-延性”过渡带的应力-深度剖面,以及一个在1857年地震破裂西北65km 处存在有“脆性-蠕变”过渡带的沿走向的应力剖面。模拟的震源空间分布在统计意义上与观测数据相类似。结果表明,断层强烈不均匀性的特征尺度范围会在断层系统的地震响应中显著表现出来。若尺度范围小,则会使震级-频度分布接近于特征地震分布,在时间上会导致如地震空区说所预言的大地震的准周期分布。而大尺度范围的地震活动性符合古登堡-里克特的震级-频度关系以及大地震的随机或是成簇的时间分布特征。数值模拟的结果表明,各种形式的震级-频度关系以及地震时间分布的统计特征都可以通过表征断层带不均匀性的尺度范围统一起来。这可以通过借助一个给定的断层带或是岩石层的结构特征给出一个清晰的物理解释,并且得到地震及断层观测资料的支持。在一些模拟中,小地震的震级-频度关系曲线明显低于自相似的古登堡-里克特关系曲线。这说明在构造加载过程中,小地震通过平滑应力的长波成分,为大地震的发生创造了条件。这一过程是通过与小地震的大量断裂有关的短波应力的消失过程来完成的。大地震的大尺度破裂的发生正好导致与上述长、短波应力的平滑与消失变化相反的过程。