断层自发破裂传播的数值模拟:速率状态摩擦定律本构参数的影响 被引量：6

1

作者 张慧茜 张海明 +1 位作者 盖增喜 黄清华 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2014年第10期3285-3295,共11页

地震往往受控于滑动面的摩擦性质,这种摩擦性质可以由速率状态摩擦定律较好地描述.速率状态摩擦定律中的本构参数a和b与动态摩擦系数相关,从而影响着同震位移与剪切应力的时空演化.本文在前人工作的基础上,采用三维边界积分方程法模拟... 地震往往受控于滑动面的摩擦性质,这种摩擦性质可以由速率状态摩擦定律较好地描述.速率状态摩擦定律中的本构参数a和b与动态摩擦系数相关,从而影响着同震位移与剪切应力的时空演化.本文在前人工作的基础上,采用三维边界积分方程法模拟速率状态摩擦定律控制下均匀全空间中平面断层的自发破裂传播过程,并详细讨论了a和b对滑动速率、剪切应力和破裂传播速度的影响.数值结果表明a和b的不同取值将导致不同的破裂行为,b-a的值越大,断层越不稳定,这种不稳定性有利于裂纹的产生与扩展.但滑动速率的时空分布不只依赖b-a,而且还与a和b的具体取值有关,断层面上滑动速率峰值与剪切破裂强度均随着a的减小而增大,随着b的增大而增大.相关结果有助于加深对断层自发破裂传播的认识. 展开更多

关键词 震源动力学 自发破裂 速率状态摩擦定律 边界积分方程法

在线阅读 下载PDF 职称材料

基于速率-状态摩擦定律的地震二维准动力学数值模拟:浅层正应力变化对断层演化的影响 被引量：3

2

作者 申文豪 薛腾飞 +1 位作者 张景发 杨建钦 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2018年第8期3198-3210,共13页

本文基于结合速率-状态摩擦定律(RSF)的二维准动力学数值模型,以半空间垂直走滑断层为研究对象,通过比较两种正应力随深度变化模型的模拟结果,研究了浅层正应力变化对断层演化参数、地震孕育过程、震后滑移传播等方面的影响.结果显示,... 本文基于结合速率-状态摩擦定律(RSF)的二维准动力学数值模型,以半空间垂直走滑断层为研究对象,通过比较两种正应力随深度变化模型的模拟结果,研究了浅层正应力变化对断层演化参数、地震孕育过程、震后滑移传播等方面的影响.结果显示,我们的数值模型在给定模型参数和约束条件下,能够完整模拟出地震周期中震间、震前、同震以及震后多个特征阶段.常数正应力模型下,动态破裂在浅层速率强化区停止,而在浅层变化正应力模型下动态破裂可以传播至自由表面,导致浅层更高的最大滑移速率和同震滑移量.两种模型下的地震矩、地震周期、平均应力降和震后滑移传播等差别不明显.两种滑移模型的傅氏振幅谱与理论K-2模型傅氏振幅谱均符合较好,且浅层变化正应力模型下的拐角波数值高于常数正应力模型,说明两种模型均符合地震同震滑移模型的运动学特征,并且浅层变化正应力模型下最终应该产生高于常数正应力模型的高频强地面运动水平.我们认为选用不同的模型参数对最终结果存在显著影响,应当根据具体问题来选择模型参数,这样才能在保证结果准确前提下有效提高计算效率. 展开更多

关键词 速率-状态摩擦定律 准动力学数值模拟 断层演化 地震参数

在线阅读 下载PDF 职称材料

基于差分Ricatti方程的Dieterich地震活动性模型构建及其应用

3

作者 尚园程 史保平 李丹 《地球物理学报》 北大核心 2025年第5期1856-1865,共10页

地震活动性模型(Dieterich模型)对研究地震成因机制及其时空动态发展过程具有重要意义.然而,传统Dieterich模型仅适用于剪切应力按照对数、线性函数规律加载/卸载以及由同震滑移导致的剪切应力阶跃等条件,面对更复杂、更贴近真实情况的... 地震活动性模型(Dieterich模型)对研究地震成因机制及其时空动态发展过程具有重要意义.然而,传统Dieterich模型仅适用于剪切应力按照对数、线性函数规律加载/卸载以及由同震滑移导致的剪切应力阶跃等条件,面对更复杂、更贴近真实情况的力学状态,例如断层上Coulomb应力率函数为谐函数或多段分段函数等函数形式时,传统模型不再适用.为此,本文提出一种新的数学处理方法,构建利用差分法处理后的Ricatti方程形式Dieterich地震活动性模型.该模型相较于传统模型在模拟、探究复杂形态Coulomb应力率函数对地震发生频率特征影响时更有优势,能够更细致再现震源区域内Coulomb应力加载速率随时间演化的细节,提升Dieterich地震活动性理论在应对复杂问题中的适用性和普适性.运用本文改进的模型对2013年前郭震群和2017年长宁诱发地震事件进行模拟与分析,揭示前郭震群阶段震源区域内Coulomb应力加载速率多次变动,辨识出长宁诱发地震期间Coulomb应力加载速率为负值的特殊时段,验证了负向Coulomb应力加载速率对地震活动频率的抑制作用.此外,该模型为深部油气和干热岩勘探开发过程中可能诱发的地震现象提供了理论依据和量化手段,有助于提高资源开发效率和降低开发成本,进一步推动勘探理论的深入发展. 展开更多

关键词 速率-状态摩擦定律 Dieterich模型 地震活动性 诱发地震

在线阅读 下载PDF 职称材料

利用应力加载速率变化与地震活动性模型探讨2013年吉林前郭震群成因的力学机制 被引量：2

4

作者 尚园程 史保平 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2022年第1期26-36,共11页

理解震群的成因机制是地震科学的重要挑战.本文基于速率-状态摩擦定律,利用Dieterich断层失稳模型探讨了流体运移、慢滑移和震群触发三者之间的关系,研究结果表明:(1)断层内部或周边流体流动所导致的孔隙压力升高可诱发断层发生慢滑或... 理解震群的成因机制是地震科学的重要挑战.本文基于速率-状态摩擦定律,利用Dieterich断层失稳模型探讨了流体运移、慢滑移和震群触发三者之间的关系,研究结果表明:(1)断层内部或周边流体流动所导致的孔隙压力升高可诱发断层发生慢滑或非地震滑移,慢滑或非地震滑移可使区域内临近断层上库仑应力加载速率升高,进而导致震群产生;(2)基于Dieterich地震活动性模型所给出的慢滑或非地震滑移触发的震群的地震发生率表达式可定量描述时间域内震群发生率上升、平稳以及快速下降的过程;(3)活动断层可能在震群期间重复失稳.针对2013年前郭震群的分析表明,前郭震群序列的产生与序列中b值(古登堡-里克特定律)的降低可能同整个震源区域内库仑应力加载速率的突然上升相关,而孔隙压力上升所诱发的断层慢滑或非地震滑移则可能为导致该区域内库仑应力加载速率上升的关键原因之一. 展开更多

关键词 速率-状态摩擦定律 Dieterich模型 震群 慢滑移 流体运移 2013年前郭震群

在线阅读 下载PDF 职称材料

龙门山断裂带上“破裂空段”的发震概率研究

5

作者 宋金 蒋海昆 +1 位作者 孟令媛 臧阳 《中国地震》 北大核心 2017年第2期219-228,共10页

本文采用分层粘弹性介质模型计算了汶川地震对芦山震中产生的库仑应力加载的影响,进而结合Dieterich(1994)提出的速率状态摩擦定律给出芦山附近区域6级地震累积发震概率随时间的变化。结果显示,2013年芦山7.0级地震时其累积发震概率达1... 本文采用分层粘弹性介质模型计算了汶川地震对芦山震中产生的库仑应力加载的影响,进而结合Dieterich(1994)提出的速率状态摩擦定律给出芦山附近区域6级地震累积发震概率随时间的变化。结果显示,2013年芦山7.0级地震时其累积发震概率达18%,说明汶川地震产生的应力扰动加速了芦山地震的发生。本文还计算了汶川、芦山2次地震对其间"破裂空段"处产生的累积库仑应力扰动的影响,结合背景地震发生率,给出了"破裂空段"处6级地震累积发震概率变化。虽然计算结果可能受到大邑地震、介质模型参数的选取和背景地震发生概率等因素影响而存在一定误差,但"破裂空段"在2次强震应力加载下累积发震概率是不断增大的,因此我们认为"破裂空段"处发生中强地震的紧迫性不断增强。 展开更多

关键词 库仑应力 速率状态摩擦定律 累积发震概率

在线阅读 下载PDF 职称材料

汶川地震对芦山地震及周边断层发震概率的影响 被引量：11

6

作者 刘博研 史保平 雷建设 《地震学报》 CSCD 北大核心 2013年第5期642-651,1,共10页

大震后区域静态库仑应力变化常常被用于解释区域地震活动性速率的变化、主震断层外余震的发生以及即将失稳断层的地震发生概率的变化.2013年4月20日芦山MS7.0地震的发生重新引起了对2008年5月12日汶川MS8.0大地震的热议.利用含(滑移)速... 大震后区域静态库仑应力变化常常被用于解释区域地震活动性速率的变化、主震断层外余震的发生以及即将失稳断层的地震发生概率的变化.2013年4月20日芦山MS7.0地震的发生重新引起了对2008年5月12日汶川MS8.0大地震的热议.利用含(滑移)速率和状态的摩擦定律,结合汶川大地震前后的地震活动性水平,定量化计算了汶川地震后雅安地区发震概率的变化,并着重解释了芦山地震发震的可能根源.此外,还对库仑应力明显增加的鲜水河断层和熊坡断层进行了发震概率的定量化计算,计算结果与中国地震台网中心的地震目录基本符合.鲜水河断层从汶川地震后至今近5年来未发生M>6.0地震,而M>6.0的发震概率已约为60%;熊坡断层自汶川地震以来尚未发生M>4.0地震,芦山地震后M>4.0的发震概率已接近90%.所以,我们认为鲜水河断层附近将成为M>6.0地震的重点防范地区,熊坡地区将来仍旧存在发生中强地震的危险性. 展开更多

关键词 含（滑移）速率和状态的摩擦定律 地震发生概率 库仑应力变化

在线阅读 下载PDF 职称材料

2008年和2014年于田地震对周边断层发震概率的影响 被引量：7

7

作者 刘博研 史保平 雷建设 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2015年第2期463-473,共11页

大震后区域静态库仑应力变化直接影响地震活动性速率的变化、主震断层外余震和即将失稳断层的发震概率的变化.利用滑移速率和状态相依赖的摩擦定律,结合2008年3月21日于田地震前后的地震活动性水平,定量计算了2008年于田地震后该地区周... 大震后区域静态库仑应力变化直接影响地震活动性速率的变化、主震断层外余震和即将失稳断层的发震概率的变化.利用滑移速率和状态相依赖的摩擦定律,结合2008年3月21日于田地震前后的地震活动性水平,定量计算了2008年于田地震后该地区周边断层发震概率的变化,着重解释了2014年于田地震发震的可能根源.此外,本文还对库仑应力明显变化的周边三条断层进行了发震概率的定量计算.贡嘎错断裂中段、贡嘎错断裂西南段和康西瓦断裂中段分别经历了发震概率先降后升、先升后升和先降后降两个阶段,充分显示了库仑应力的细微变化造成的周边断层的危险性的变化.这三条断裂发生7.0级以上地震的发震概率超越95%均需要500年左右;贡嘎错断裂西南段发生中强地震的可能性较大,而康西瓦断裂中段活跃度较低. 展开更多

关键词 滑移速率和状态相依赖的摩擦定律 发震概率 库仑应力变化

在线阅读 下载PDF 职称材料

M_S8.1昆仑山口西地震和M_S8.0汶川地震余震序列的时空分布特征和持续时间的对比 被引量：9

8

作者 刘博研 史保平 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2012年第6期1942-1951,共10页

2001年MS8.1昆仑山口西地震和2008年MS8.0汶川地震发生在同一构造单元,但其余震序列无论在个数、空间分布,还是持续时间上都表现了显著的差别.余震通常由主震区域内背景场地震活动性受到的扰动所引起,这样的扰动则来自于主震造成的应力... 2001年MS8.1昆仑山口西地震和2008年MS8.0汶川地震发生在同一构造单元,但其余震序列无论在个数、空间分布,还是持续时间上都表现了显著的差别.余震通常由主震区域内背景场地震活动性受到的扰动所引起,这样的扰动则来自于主震造成的应力场状态的变化.本文从滑移速率和状态相依赖的摩擦定律(Rate-and State-Dependent Friction Law)出发,结合区域主震前后的地震活动性资料,定量地估算了这两个大地震后余震序列可能的持续时间,并对不同模型所得的结果进行了比较和对比.结果表明,汶川地震余震持续时间约为昆仑山口西地震余震持续时间的20倍,这是由于昆仑山口西地震和汶川地震余震序列的个数和持续时间不仅与地震成核过程的状态变化有关,还与作用在断层面上的正应力σN和剪应力加载速率τ.的大小有关.主震前后剪应力速率τ.的差别导致了在相同大小应力扰动ΔCFS之后的余震的活动性变化率的明显不同,导致了所触发的余震的个数和余震序列的持续时间的巨大差别.通过对昆仑山口西地震和汶川地震余震序列的时空分布特征和持续时间的定量化认识,可以为地震灾害定量评估提供合理和有益的物理参数. 展开更多

关键词 应力降 滑移速率和状态相依赖的摩擦定律 库仑失稳准则

在线阅读 下载PDF 职称材料

汶川地震同震过程中断层间的相互作用及子事件触发 被引量：1

9

作者 刘博研 史保平 《地震地质》 EI CSCD 北大核心 2015年第2期357-373,共17页

2008年5月12日四川省汶川县发生了MS8.0地震,造成了重大人员伤亡和经济、财产损失。主震发生时,龙门山前山断层从都江堰北部开始错动,且一直破裂到地表。文中利用映秀-北川断裂和灌县-江油断裂2条断层模型,着重研究和深入探讨了灌县-江... 2008年5月12日四川省汶川县发生了MS8.0地震,造成了重大人员伤亡和经济、财产损失。主震发生时,龙门山前山断层从都江堰北部开始错动,且一直破裂到地表。文中利用映秀-北川断裂和灌县-江油断裂2条断层模型,着重研究和深入探讨了灌县-江油断裂的触发机制以及错动时间。基于滑移速率和状态相依赖的摩擦定律,推导出在剪切应力变化率是否为0的2种前提条件下的地震触发新破裂的时间延迟,并将结果用于研究汶川地震是否为同震触发的时间过程上。通过与绵竹清平和汶川卧龙台站的三分量强地面运动观测资料计算出的时间延迟对比,得到了比较一致的结果,发现汶川地震中映秀-北川断裂和灌县-江油断裂并不是同时错动的,而是映秀-北川破裂产生的动态地震波(应力波)触发了灌县-江油断裂,并造成了几秒到十几秒的触发延迟。 展开更多

关键词 汶川地震 地震触发 滑移速率和状态相依赖的摩擦定律 时间延迟

在线阅读 下载PDF 职称材料

基于应力变化的余震预测研究综述

10

作者 贾若 蒋海昆 《中国地震》 北大核心 2013年第3期293-305,共13页

本文综述了近期基于应力变化的余震预测研究进展。从静态库仑应力触发角度对余震活动进行研究主要有两方面内容:一是集中于对主震触发余震机制的讨论;二是通过库仑应力变化分布区域与实际余震空间分布区域进行定性对比分析,并认为库仑... 本文综述了近期基于应力变化的余震预测研究进展。从静态库仑应力触发角度对余震活动进行研究主要有两方面内容:一是集中于对主震触发余震机制的讨论;二是通过库仑应力变化分布区域与实际余震空间分布区域进行定性对比分析,并认为库仑应力变化为正的区域有利于余震的发生,而这正是余震活动与应力变化高度相关的重要证据。基于应力变化的余震预测研究的一项重要内容是建立在速率-状态依从摩擦定律基础之上的余震活动率研究,即将余震活动作为一系列独立的成核事件,从断层滑动触发地震的角度出发,结合静态库仑应力计算得到余震活动率,进而得到余震数量。余震活动是与断层面上的剪切应力变化相关联的,余震活动率受应力扰动幅度、应力加载速率、断层基本物理特性及背景地震活动率等的影响。考虑到"余震激发余震"的特点,ETAS模型与静态库仑应力变化相结合的研究方法提供了一条统计学与物理学方法(应力变化)相结合、能更完整地刻画余震活动特征的可能途径。 展开更多

关键词 余震 速率-状态依从摩擦定律 静态库仑应力变化 余震活动率 背景地震活动率

在线阅读 下载PDF 职称材料

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

11

作者 解孟雨 史保平 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2016年第2期593-605,共13页

基于大量岩石力学实验,Dieterich和其他研究者(Dieterich,1978;Ruina,1983)首先提出了描述岩石摩擦过程的速率-状态摩擦定律(R-S摩擦定律).如今R-S摩擦定律已成为研究地震成核等地震演化机制的有效手段.Dieterich(1992,1994)最早提出了... 基于大量岩石力学实验,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模型则不能完整地给出受到静态剪应力扰动后断层失稳时间提前或推后的估计值. 展开更多

关键词 速率-状态摩擦定律 静态剪应力扰动 Coulomb应力失稳模型 Dieterich模型 失稳时间提前量 失稳时间推后量

在线阅读 下载PDF 职称材料

断层演化过程中的周期和非周期地震滑移及非地震滑移——一维断层模型分析 被引量：3

12

作者 解孟雨 史保平 《地震学报》 CSCD 北大核心 2016年第4期590-608,658,共19页

本文首先根据Dieterich和Ruina提出的含速率和状态的摩擦定律(Dieterich-Ruina定律),基于一维弹簧-滑块模型推导了地震复发周期的解析表达式,然后将该近似解与数值模拟结果以及Barbot等的相关研究进行了对比分析.此外,本文还利用数值模... 本文首先根据Dieterich和Ruina提出的含速率和状态的摩擦定律(Dieterich-Ruina定律),基于一维弹簧-滑块模型推导了地震复发周期的解析表达式,然后将该近似解与数值模拟结果以及Barbot等的相关研究进行了对比分析.此外,本文还利用数值模拟与理论分析研究了断层周期和非周期演化的力学成因机制以及非地震滑移形成的另类力学机制,并讨论了一维弹簧-滑块模型的优点及其局限性.结果表明:1震后滑移和自加速/成核阶段的持续时间在整个演化过程中不能被忽略;2在修正后的复发周期模型中,复发周期的长短除了与断层特征尺度、作用于断层面上的有效正应力和远场加载速率相关外,还受Dieterich-Ruina定律中摩擦参数的取值以及临界滑移距离的影响;3当给定各个物理参数和几何参数时,目前所得到的解析近似解可以很好地估计地震的复发周期,其相对误差可小于5%;4在断层演化过程中,施加剪切应力加载会产生非周期的地震滑移,而在自加速/成核阶段后期或震后滑移阶段早期,施加较大的剪切应力加载,则会出现非地震滑移. 展开更多

关键词 含速率和状态的摩擦定律 地震循环 地震复发周期 非地震滑移

在线阅读 下载PDF 职称材料

基于Dieterich地震发生率模型分析前震成因的力学机制 被引量：4

13

作者 尚园程 史保平 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2020年第9期3355-3369,共15页

前震是地震前兆观测中的一个重要物理量,前震的研究对地震预测的发展尤为重要,前震序列的典型特征包括加速发生的小震与古登堡-里克特定律(Gutenberg-Richter Law)中的b值变小.本文在Dieterich提出的地震发生率模型的基础上,探讨前震发... 前震是地震前兆观测中的一个重要物理量,前震的研究对地震预测的发展尤为重要,前震序列的典型特征包括加速发生的小震与古登堡-里克特定律(Gutenberg-Richter Law)中的b值变小.本文在Dieterich提出的地震发生率模型的基础上,探讨前震发生的力学成因机制及前震-主震-余震时域演化特征.模型分析结果表明,前震震源区断层及其断层周边剪应力加载速率的变化可能是前震发生及地震发生率变化的一个关键原因.由于前震震源区次级断层之间剪应力加载速率受来自主震成核过程中断层自加速滑移的影响升高,从而可导致这些次级断裂的加速失稳,即加速前震的发生.当震源区内由前震所产生的静态应力扰动不可忽略时,应力扰动和主震成核的共同作用也可对后续前震发生率产生影响.当正向应力扰动出现时,后续的前震序列的地震发生率会出现陡增,随后其地震发生率逐渐下降.而当负向应力扰动出现时,地震发生率会出现陡降,然后再次逐渐上升.基于Kostrov模型,本文得到了剪应力加载速率与古登堡-里克特定律中b值的关系式,结果表明前震序列中b值的减小与前震区内的剪应力加载速率的上升有关. 展开更多

关键词 速率-状态摩擦定律 Dieterich模型 前震 B值 台湾花莲地震

在线阅读 下载PDF 职称材料