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

1

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

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

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

在线阅读 下载PDF 职称材料

摩擦生热对断层演化过程的影响:Ruina-和Chester-Higgs-模型的对比分析

2

作者 高雅琪 史保平 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2020年第6期2297-2310,共14页

利用速率-状态摩擦定律(Rate-and State-Dependent Friction Law:简称RSF定律),结合McKenzie-Brune摩擦生热模型,本文分别从Ruina提出的RSF定律和Chester-Higgs提出的RSF定律出发,通过一维弹簧-滑块模型,采用四阶变步长的Dormand-Princ... 利用速率-状态摩擦定律(Rate-and State-Dependent Friction Law:简称RSF定律),结合McKenzie-Brune摩擦生热模型,本文分别从Ruina提出的RSF定律和Chester-Higgs提出的RSF定律出发,通过一维弹簧-滑块模型,采用四阶变步长的Dormand-Prince算法,对断层演化过程进行了数值模拟,探讨了摩擦生热对断层演化进程的影响.模拟结果显示,与Ruina-模型相比,Chester-Higgs-模型在断层高速滑动时存有更大的摩擦强度,表明摩擦生热对断层具有一定的强化作用,且同临界滑移距离的取值相关.而且,Chester-Higgs-模型在失稳时的断层面温度远远低于Ruina-模型,表明摩擦生热在断层演化过程中能抑制断层面温度的剧烈升高,且正应力和临界滑移距离越大,两种模型的温差越为明显,而断层的刚度和尺度则对温度的影响很小.模拟两种模型周期演化过程的结果表明,在相同的初始条件下,Chester-Higgs-模型给出的断层失稳周期明显比Ruina-模型更短,说明摩擦生热对断层自身演化最显著的影响是较大地缩短了地震重复发生周期.当断层进入周期性演化后,Chester-Higgs-模型给出的摩擦强度大于Ruina-模型,且对第一次非周期性失稳的摩擦强度和剩余应力的继承性更好.另外,由ChesterHiggs-模型给出的静态应力降远小于Ruina-模型给出的结果,所对应的单个事件的滑移量也小于Ruina-模型. 展开更多

关键词 Ruina-模型 Chester-Higgs-模型 McKenzie-Brune摩擦生热模型 一维弹簧-滑块模型 摩擦生热 速率和状态相依赖的摩擦

在线阅读 下载PDF 职称材料

远场应变速率对非均质断层滑动成核影响研究 被引量：2

3

作者 王威 任青文 张丽峰 《水电能源科学》 2008年第3期50-53,共4页

为进一步探讨断层滑动成核机理,以无限体弹性位错理论为基础,结合速率—状态相关摩擦准则,采用数值方法研究了远场应变加载率对非均质断层成核的影响及各物理量在断层失稳过程中的变化。结果表明,对于非均质断层,远场应变加载率对断层... 为进一步探讨断层滑动成核机理,以无限体弹性位错理论为基础,结合速率—状态相关摩擦准则,采用数值方法研究了远场应变加载率对非均质断层成核的影响及各物理量在断层失稳过程中的变化。结果表明,对于非均质断层,远场应变加载率对断层成核的时间和后期的速率、摩擦应力均有影响。远场应变加载率越高,断层成核时间越早,成核后期速率越高,应力降越小,成核区的尺寸略有增加。 展开更多

关键词 速率-状态相关摩擦准则 断层 成核 弹性位错 非均质

在线阅读 下载PDF 职称材料

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

4

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

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

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

在线阅读 下载PDF 职称材料

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

5

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

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

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

在线阅读 下载PDF 职称材料

CO_(2)地质封存诱发断层活化和流体泄漏模型及数值研究 被引量：2

6

作者 蒋林桐 马天然 +4 位作者 李采 郭朝斌 何庆成 马晓冬 张诗淮 《地质学报》 EI CAS CSCD 北大核心 2024年第11期3418-3432,共15页

地下注入二氧化碳(CO_(2))可能通过改变孔隙压力和地下应力而引发地震。深入了解流体注入导致的断层再活化的基本机理对于评估与地下流动过程相关的地震及流体泄漏风险具有重要意义。基于此问题,本文推导了地层和断层内CO_(2)与水两相... 地下注入二氧化碳(CO_(2))可能通过改变孔隙压力和地下应力而引发地震。深入了解流体注入导致的断层再活化的基本机理对于评估与地下流动过程相关的地震及流体泄漏风险具有重要意义。基于此问题,本文推导了地层和断层内CO_(2)与水两相渗流控制方程,建立了模拟CO_(2)注入诱发断层活化及泄漏的全过程数值模型,重点分析了不同的断层属性和CO_(2)注入工况对断层活化及流体泄漏的影响。数值结果表明,在断层活化后,其渗透率显著提升,最大增幅可达106倍,将断层滑移区从原本的阻碍渗透的屏障状态转变为主导渗流通道。不同的初始切向渗透率对断层活化时的力学特征有显著影响,初始切向渗透率越大,断层滑移量越大。具有较强的剪切刚度的断层对注入诱发的活化表现出更强的抵抗能力。此外,CO_(2)注入速率影响了断层活化的地震和CO_(2)泄漏。低注入速率在相同的注入量下引起的断层滑移量更小,但伴随着更多的CO_(2)泄漏;而高注入速率下,断层滑移量更大,但CO_(2)泄漏量较少。 展开更多

关键词 CO_(2)地质封存 两相流 状态-速率摩擦本构 断层活化

在线阅读 下载PDF 职称材料

断层几何形态对阿拉斯加中南部俯冲带慢滑移特征的影响 被引量：2

7

作者 李昊天 周仕勇 《地震学报》 CSCD 北大核心 2019年第6期681-694,808,共15页

本文采用三种不同的俯冲带几何模型,在速率-状态依赖型摩擦律和准动态算法的框架下,对阿拉斯加库克湾的慢滑移事件进行了数值模拟,以探究断层几何形状对慢滑移特征的影响。结果表明:几何因素对慢滑移的时空演化有较大影响;慢滑移区域的... 本文采用三种不同的俯冲带几何模型,在速率-状态依赖型摩擦律和准动态算法的框架下,对阿拉斯加库克湾的慢滑移事件进行了数值模拟,以探究断层几何形状对慢滑移特征的影响。结果表明:几何因素对慢滑移的时空演化有较大影响;慢滑移区域的宽度对数值模拟的结果起着至关重要的作用;断层几何形态更平缓的区域将导致更大、更快的事件。这一结果有助于我们进一步了解慢滑移的成因以及断层几何形态对慢滑移时空演化的影响。 展开更多

关键词 慢滑移 阿拉斯加中南部俯冲带 断层几何效应 速率-状态依赖型摩擦律

在线阅读 下载PDF 职称材料

M_s8.0汶川地震对M_s7.0芦山地震成核失稳的影响Mw 被引量：6

8

作者 解朝娣 雷兴林 +3 位作者 吴小平 虎雄林 赵小艳 朱荣欢 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2014年第6期1825-1835,共11页

依据速率-状态依赖性摩擦本构关系,并结合中下地壳和上地幔的黏弹性松弛效应,以震后库仑破裂应力变化和同震动态库仑破裂应力变化的计算为基础,模拟了两种应力变化对芦山地震断层的失稳发震时间的影响,研究了2008年Ms8.0汶川地震与2013... 依据速率-状态依赖性摩擦本构关系,并结合中下地壳和上地幔的黏弹性松弛效应,以震后库仑破裂应力变化和同震动态库仑破裂应力变化的计算为基础,模拟了两种应力变化对芦山地震断层的失稳发震时间的影响,研究了2008年Ms8.0汶川地震与2013年Ms7.0芦山地震之间的触发关系.计算得到,汶川地震在芦山地震震源断层面上产生的动态应力变化的峰值为0.127 MPa;此外,经过近5年中下地壳和上地幔的黏弹性松弛效应,芦山地震震源断层面上受到的震后应力变化值为0.025 MPa.结果表明,芦山地震的震源断层在应力积累逐渐接近临界状态的某一特定时期内,受到了汶川地震产生的动态应力变化、静态应力变化以及黏弹性松弛效应造成应力变化的共同触发作用,且动态应力的延迟触发作用可能更为显著.最后对芦山地震之后研究区域的应力变化场进行了初步探讨. 展开更多

关键词 汶川地震 芦山地震 库仑破裂应力变化 速率-状态依赖性摩擦本构关系

在线阅读 下载PDF 职称材料

地震动态应力触发研究进展 被引量：17

9

作者 张彬 杨选辉 陆远忠 《西北地震学报》 CSCD 北大核心 2008年第3期298-303,共6页

本文从以下几方面阐述了近几年来地震动态应力触发的研究进展:地震动态应力触发问题的提出;地震动态应力触发存在的证据(包括地震频度图中探寻被触发震群和独立地震、数字滤波后新发现的被触发地震和形变仪器记录的同震变化和应力阶跃... 本文从以下几方面阐述了近几年来地震动态应力触发的研究进展:地震动态应力触发问题的提出;地震动态应力触发存在的证据(包括地震频度图中探寻被触发震群和独立地震、数字滤波后新发现的被触发地震和形变仪器记录的同震变化和应力阶跃等);地震动态触发的研究现状和机理等。还提出了目前动态应力触发研究中存在的问题,并对其近期的发展进行了展望。 展开更多

关键词 动态应力触发 库伦摩擦破裂 速率-状态模型

在线阅读 下载PDF 职称材料

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

10

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

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

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

在线阅读 下载PDF 职称材料

龙门山断层地震周期及其动力学过程模拟研究 被引量：5

11

作者 马林飞 陶玮 +2 位作者 张永 曾明会 郑茜 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2018年第5期1824-1839,共16页

在断层面上引入速率-状态相依摩擦本构关系、考虑铲形逆冲断层几何结构特征、断层下盘和上盘中下地壳及上地幔为黏弹性介质、上盘上地壳为弹塑性介质,本文用二维有限元动力学模型模拟了龙门山断层上大震准周期复发行为、分析了断层上地... 在断层面上引入速率-状态相依摩擦本构关系、考虑铲形逆冲断层几何结构特征、断层下盘和上盘中下地壳及上地幔为黏弹性介质、上盘上地壳为弹塑性介质,本文用二维有限元动力学模型模拟了龙门山断层上大震准周期复发行为、分析了断层上地震孕育位置、地震周期不同阶段的应力/应变场演化特征.不同于近垂直走滑断层上的地震周期行为,大陆铲形逆冲断层上的构造应力的积累和释放过程更复杂、有其独特性.我们得到如下认识:(1)铲形逆冲断层上的地震复发是准周期行为.(2)龙门山断层最大库仑应力位于断层17~20km深处,应力长期积累和同震释放都在此深度最大,说明地震会在此处孕育、发动.(3)在断层破裂的深部和浅部,同震滑动大小和构造应力释放大小并非同步,而是差异悬殊.(4)地震仅部分释放区域积累的应变能,断层上盘上地壳顶部和底部的褶皱、破裂等永久变形形式也是释放应变能的重要形式.(5)应变能密度增量的演化图像分为:震间、同震、震后期,清晰反应了龙门山断层附近的地震动力学过程.(6)地震发生除释放能量外,同时也对近断层的中地壳和断层底部有很大的应变能加载;这些加载,在震后期可能通过震后滑移、余震或中下地壳乃至上地幔的驰豫形变用几十年时间释放.以上对大陆内铲形逆冲断层上变形特征的了解,有助于我们在其地震周期行为中评估地震危险性. 展开更多

关键词 铲形逆冲断层 龙门山断层 地震周期行为 速率-状态相依摩擦本构关系 动力学模型

在线阅读 下载PDF 职称材料

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

12

作者 解孟雨 史保平 《地球物理学报》 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 职称材料

基于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 职称材料