鲜水河断裂带库仑应力演化与强震间关系

以鲜水河断裂带为研究区,首先验证了该断裂带上1893年以来M6.7以上地震的相互触发作用,然后采用更符合实际的分层黏弹介质模型研究强震震后黏滞松弛引起的库仑应力变化对后续地震的影响,并基于负位错理论计算鲜水河断裂带10个断层段的震间长期构造加载作用引起的断层上的应力积累.在此基础上,讨论同震、震后、震间效应引起的累积库仑应力变化与区域强震活动的关系,给出断层上库仑应力随时间的演化.结果表明,鲜水河断裂带上1893年以来发生的7次强震均在其前面一系列强震及构造应力加载的驱使下发生,同震、震后、震间三方面效应均引起了鲜水河断裂带不可忽略的库仑应力变化.对断层上的库仑应力状态的研究,可进一步为揭示地震的发生规律、找寻危险断层段提供线索.     

盲逆断层震后变形耦合机制研究——以2017年伊朗MW7.3地震为例

震后变形是地震周期中一个重要阶段,深入研究其产生机制有助于提高对地震演化过程、断层性质及地震危险性评估的认识.2017年伊朗M_W7.3地震InSAR观测资料丰富,为研究扎格罗斯造山带山前褶皱带盲逆冲断层震后变形机制提供了切入点.合理详细的震源模型及高精度的形变观测是开展震后余滑或黏弹松弛研究的首要前提,本文采用多视角InSAR资料联合远场波形数据反演该地震的破裂滑动分布,并解算得到伪三维地表位移,利用InSAR时序分析提取震源区主震后一年半内变形特征,显示震源区震后形变特征与同震位移场类似,LOS向形变速率最大值约8 cm·a^-1,震后变形预示震源区可能存在明显的震后应力调整现象.利用震后半年内形变资料约束的纯运动学余滑模型表明同时期释放能量为矩震级M_W6.7,进一步探讨应力驱动震后余滑及下地壳黏弹性松弛对震后变形的贡献,基于分层黏弹模型的模拟计算表明震后余滑和黏弹松弛效应的耦合模型可以更好地解释地表震后变形特征,其中震后余滑主要分布在同震破裂区的上方浅部,对主震后一年半内震后地表变形起主导作用,震源区壳下黏滞系数量级下...     

1990年青海共和7.0级地震震后垂直形变研究

根据1990年青海共和地震震后地表垂直形变,通过模型拟合得到了支配共和地区震后形变场时空演化的形变源及其力学机制.分析穿过断层的震前1期和震后6期水准数据,结果表明震后垂直形变具有以下特征:①震后震区上盘继续发生继承性的大幅度上升,其中震后头一年上升速率最大;②震后上升区范围显著,范围随时间变化不大,但较同震形变上升区范围增大;③震后相邻测站高差观测值的时间序列明显具有对数衰减特征或指数衰减特征,衰减特征时间分别为0.165年和1.344年.本文还发展了一个利用水准数据与连续介质位错模型研究震后形变机制的新方法.该方法用相邻水准点之间的原始高差观测值而非它们相对参照点的积分值来约束连续介质位错模型,可以有效减少误差累积带来的偏差并充分利用观测数据.利用这一方法的初步分析结果表明,断层震后滑移和介质黏弹性松弛共同导致了共和地震震后形变.前者表现为发生在断层面及其延伸部分的滑移,特别是位于主破裂上方沉积层内的滑移;后者则表现为下地壳与上地幔内的黏弹性松弛,黏滞系数为1020Pa.s量级.     

2010年MS7.1级玉树地震同震库仑应力变化以及对2011年MS5.2级囊谦地震的影响

2010年4月14日青海省玉树藏族自治州发生M_S7.1级地震.和传统的板内地震相比,玉树M_S7.1级地震的余震具有数量少、震级大的特点.研究玉树地震主震与余震之间的关系,对于我们了解余震的发震机理具有十分重要的参考价值.本文利用弹性位错理论和分层岩石圈模型,计算玉树地震引起的同震及震后黏弹松弛应力场变化,讨论M_S7.1级玉树地震对余震分布的影响以及与2011年囊谦M_S5.2级地震之间的触发关系.结果显示,玉树地震导致了四处明显的库仑应力增强的扇区,2010年4月13日至6月17日的870次M_L>1.0级余震主要分布于主震破裂面附近区域以及破裂面东北端的应力增强扇区.分析玉树地震对余震分布的影响时,有效摩擦系数以及计算深度的选取对计算结果的影响较小,是否考虑区域构造应力场的影响较大.考虑区域构造应力场时,占总数86.7%的余震位于库仑应力增强区,地震应力触发理论较好地解释了余震的分布.选取囊谦地震震源机制解的两个节面作为库仑应力计算中的接收断层参数,并且考虑不同黏滞系数下的玉树地震同震及震后黏弹松弛效应,模型计算结果均表明囊谦地震位于玉树地震所导致应力影区,仅依靠地震的静态、震后黏弹松弛应力触发理论,无法解释囊谦地震的发生,说明该次地震可能是一次独立的事件.     

鲁甸地震同震与震后形变、重力时空分布的理论模拟

基于位错理论,考虑重力和黏弹性的影响,在分层介质模型下计算鲁甸地震引起的同震、震后形变和重力变化.结果表明形变和重力的显著变化主要发生于断层在地表投影附近区域.同震形变场显示发震断层有明显的走滑性质.考虑黏弹松弛效应,随着时间的推移,震后形变和重力有了明显改变,同震效应为正的区域得到加强,为负区域进一步减弱.震后松弛效应的影响范围相比同震明显增加.在靠近断层的GPS观测台站处,计算了由黏弹松弛效应引起的震后形变和重力时间序列.震后松弛效应引起的重力变化在50年之后均达到同震水平,除了NJ13的纬向、垂直位移,NJ16的垂直位移,NJ15的径向位移,其余台站的所有震后形变都超过1mm.观测台站的震后重力和垂直位移时间序列在震后100年趋于稳定,纬向位移和经向位移在震后50年趋于稳定.     

M_s8.0汶川地震对M_s7.0芦山地震成核失稳的影响Mw

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

MS8.0汶川地震对MS7.0芦山地震成核失稳的影响

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

利用同震和震后位移数据联合反演2011年日本MW9.0地震同震断层滑动

基于黏弹性球体位错理论, 联合陆地和海底同震GPS数据以及日本本岛330个陆地GPS站点5~10年的震后数据, 反演了日本M_W9.0地震的断层滑动模型, 提升了断层滑动分布在细节上的合理性。 首先, 基于日本本岛330个陆地GPS站点震前2年和震后10年的连续观测数据, 获取了日本M_W9.0地震震后5~10年的年平均位移, 该时段的位移几乎完全由地幔黏弹性松弛效应引起; 接着, 利用黏弹性球体位错理论对震后5~10年的位移进行反复拟合, 确定了日本M_W9.0地震震源及周边地区的地幔黏滞性系数最优解(9.0×10¹⁸ Pa·s)。 然后, 联合同震和震后位移数据, 引入黏弹性位错格林函数, 反演了2011年日本M_W9.0地震的断层滑动分布。 结果表明, 该地震同震破裂的最大值达到了62.72 m, 同震滑动的总地震矩为4.48×10²² Nm, 相应的矩震级为M_W9.03。 由于黏弹性松弛效应引起的震后位移中包含了同震破裂的信息, 基于黏弹性球体地震位错理论, 联合同震和震后位移数据反演断层同震破裂, 有效提高了日本M_W9.0地震断层滑动分布的可靠性。 最后, 本文提出的反演方法为同震观测结果缺乏的大地震震后科考提供了理论支撑: 在大地震发生之后, 即使在同震期间没有足够的观测数据, 也可以在震后通过对震源区的加密观测积累的震后数据, 使用本文提出的反演方法优化同震断层滑动模型。     

从应力演化角度讨论龙门山区域7级强震序列对芦山地震的影响

计算包括同震静态库仑破裂应力变化以及震后黏滞松弛引起的应力变化等,可以更好地解释余震分布、地震序列等地震观测结果。在芦山地震之前,从1900年以来龙门山区域发生了4次7级以上地震,分析这几次地震的同震应力变化以及震后黏弹性松弛对芦山地震的产生的影响,芦山地震以后区域断裂带上的应力伴随强震如何演化,芦山地震与汶川地震的破裂空段呈现何种应力状态,探讨这些问题可能会为了解芦山地震震源处震前的应力状态及该区域未来地震风险评估提供一定的依据。文中采用有限元数值模拟方法,根据地质构造、速度、密度结构深部反演结果以及GPS及应力观测资料等,建立龙门山地区三维黏弹性有限元模型进行研究。模拟结果显示:龙门山断裂带南段及鲜水河断裂带的库仑应力年变化速率在研究区域中相对更高,这与研究区域的地震活动性一致。芦山地震的前4次地震,除叠溪地震外,康定、松潘、汶川等3次地震在芦山地震震中位置产生的同震库仑破裂应力变化大于0,表明这3次地震可能促进了芦山地震的发生,汶川地震的同震库仑破裂应力超过了0.01MPa,同震触发效应十分显著。震间的黏弹性松弛对芦山震源处起加载作用,从1900年以来这种持续的加载作用也超过0.01MPa,因此在模拟应力演化的时候,介质的黏弹性松弛效应不能被忽略。从库仑破裂应力的角度计算龙门山区域断层的应力演化,可以发现龙门山断裂带上汶川地震和芦山地震破裂的空段,在芦山地震之后仍然属于相对应力水平较高的区域。     

唐山地震后发震断层和周围地区的地壳形变

介绍了设于唐山地震发震断层地表裂缝的形变台站，在地震后18年间观测到的断层形变．结果表明，发震断层的垂直和水平形变都集中在大震之后的7年之内；1983年是震后形变的转折期．大面积水准显示的形变与断层形变强度定性相符，目前呈现出不易于积累应力的松弛状态，近期不大可能发生较强地震．      

研究地震和岩石破裂现象的非线性科学方法

本文以地震孕育和岩石破裂过程的阶段(包括裂纹系形成、变形局部化和断层形成、断层活动和松弛)为线索,阐述了非线性科学方法(包括分形几何学,自组织临界现象和混沌动力学)的应用,提出了地震孕育过程和岩石破裂各阶段的普适性特征、机理和要着重研究的课题.     

共轭地震与地壳介质内摩擦系数的测定

安德逊理论表明地震的共轭破裂现象可以用库仑准则解释,断层面夹角与地下岩石内摩擦系数之间存在简单的定量关系。基于本文提出的通过共轭地震计算介质内摩擦系数的方法,对国内外19个具有X形或不完整的X形（L形）共轭破裂特征的地震进行系统性分析,与震源机制、区域应力轴方向等资料进行对比,确定共轭断层面,进而确定共轭面夹角并计算岩石的内摩擦系数。研究结果进一步证实了天然地震的共轭破裂与实验室岩石破裂具有较好的一致性,给出了库仑准则与安德逊理论在大尺度下也成立的证据,表明通过共轭角确定地下介质的内摩擦系数是对地下介质特性认知的有益补充。      

应力时空演化揭示出的汶川地震前亚失稳过程

近年来从岩石实验研究针对应力水平进入峰值强度后至断层快速失稳前的短临阶段变形物理场的时空演化过程提出了亚失稳阶段的概念,该阶段的基本特征为应力由峰值前的积累为主转为峰值后稳定释放并加速至最终快速失稳,这一过程伴随着其他变形物理场的相应转变.把握该阶段应力场的转变特征对认识未来地震危险区及其变形阶段至关重要.以2008年汶川地震为例,运用基于区域多个地震震源机制反演定量应力张量的方法对强震前是否存在可观测的亚失稳应力状态进行了研究.获得的震源附近紫坪铺水库区域临震前应力状态时空演化过程显示,水库测震台网开始运行的2004年8月至2007年6月期间该区域归一化有效平均应力水平相对较低且均匀分布,随时间无明显变化;而2007年6月后至汶川主震发生前应力水平局部显著增强.基于不同时段的定量应力进一步获得了随时间演化的增量场,显示2007年6月后除了汶川附近和彭县附近为应力增强区,其他大部分区域表现为幅值相对较低的应力释放,表明研究区在2007年6月后应力水平明显弱化;而在2007年6月之前的应力增量场却显示区域应力随时间以积累为主.根据应力演化特征及分时段的应力增量场可推断孕震区在2007年中期达到强度极限,峰值前区域应力随时间平稳变化并以积累为主,而峰值后则转变为明显的应力释放.该现象显示汶川地震前区域应力存在由积累为主转为释放为主的可观测阶段,符合实验研究中断层由稳态进入亚失稳状态的应力演化特征.结合其他研究结果,推测汶川地震前震源断层及附近区域存在联动加速变形过程,与震前亚失稳应力状态相关.     

分形方法在地震序列类型早期判别中的应用

着重讨论了地震序列的空间分布随时间变化的特性，利用分形方法提取大震后余震的空间分维信息，从而达到早期判别地震序列类型的目的。以１９６６年以来定位精度较好、余震记录较全的１２个７级以上大震和大同６级强震序列为资料，计算了震后１个月内余震随时间变化的空间分维数。结果发现，不同类型的地震序列分维数不同。一般震后５－６天，主震型序列的分维数在１．１－１．２之间，多震型序列的分维数在１．３以上。利用前７天的     

岩石分形断裂的统计理论

本文研究了岩石分形断裂的统计理论.认为岩石断裂是岩石内部微裂缝成核、形成和传播的非平衡过程,导出了描述这个过程的微分方程.运用位错塞积成核机理,得到了岩石发生分形断裂的微裂缝大小的统计分布函数和平均长大速率.从理论上得出了岩石非弹性体积应变与微裂缝数目成正比,从而解释了实验室研究的结果.另外发现岩石的分维是决定岩石脆性的一个重要因素,分维值愈小,岩石愈脆.但岩石的分维值与岩石强度无明显关系.最后,从微裂缝大小的统计分布函数出发,导出了b值与岩石分维值的关系.     

地震断裂系统的分形和在地震预报中的应用

本文基于地震断层具有自相似结构,从分形概念出发,依据地震资料提出了适合于前震和余震的分形模型,在此基础上讨论了前震和余震的b值特征,提出了地震矩不均衡度的概念。本文还探讨了利用分数维和地震矩不均衡度进行地震预报的可能性。研究发现,正常情况下,地震矩不均衡度呈现低值,大震发生前出现高值异常,震后降低;在空间上,震前2—4年其异常范围大,震前一年异常区逐渐缩小,大震发生在异常集中区内或其边缘。     

Towards a new view of earthquake phenomena

Recent advances in the theory of fracture and fragmentation are reviewed. Empirical laws in seismology are interpreted from a fractal perspective, and earthquakes are viewed as a self-organized critical phenomenon (SOC). Earthquakes occur as an energy dissipation process in the earth's crust to which the tectonic energy is continuously input. The crust self-organizes into the critical state and the temporal and spatial fractal structure emerges naturally. Power-law relations known in seismology are the expression of the critical state of the crust. An SOC model for earthquakes, which explains the Gutenberg-Richter relation, the Omori's formula of aftershocks and the fractal distribution of hypocenters, is presented. A new view of earthquake phenomena shares a common standpoint with other disciplines to study natural complex phenomena with a unified theory.     

东昆仑活动断裂带大地震之间的黏弹性应力触发研究

对青藏高原北部东昆仑破裂带大地震之间的应力转移和断层相互作用进行研究. 考虑1937年以来沿此破裂带发生的5个M≥7的地震：1937年M7.5花石峡地震，1963年MS7.1都兰地震，1973年MS7.3玛尼地震，1997年MW7.5玛尼地震和2001年MW7.8可可西里地震，模拟了黏弹性成层介质中地震断层错动产生的应力演化过程，并计算了在后续地震破裂面上产生的库仑破裂应力变化. 结果表明，前面4个地震均造成2001年可可西里地震断层面上库仑破裂应力的增加，并且中地壳和下地壳的黏弹性松弛效应使得库仑破裂应力场随着时间的推移而逐渐加强. 在计算过程中定量估计了可可西里地震发生时前面4个地震同震形变和黏弹性松弛导致可可西里地震破裂面上库仑破裂应力变化之间的比值，发现前3个地震由黏弹性松弛造成的变化远远大于同震形变所造成的变化. 可可西里地震之后应力场的模拟表明东昆仑断层中段的东大滩－西大滩断层段（位于可可西里地震破裂以东及都兰地震以西）的库仑破裂应力显著增加，变化值达0.05～0.1 MPa，预示这一地区地震危险性的增加.     

山东胜利油田角07井注水地震序列的演化及其机制

1985年12月至1986年1月发生的山东省胜利油田角07井注水地震序列,经历了特有的演化过程,即全序列可分解为4个子序列,且子序列逐渐由震群型向主余型过渡.前期地震(最大地震之前)呈升级趋势,地震b值前期高于后期,地震时间结构的分数维D则前期低于后期,以及较大地震的lg△t-lgt前期较为紊乱,而后期线性相关等.计算表明,注水后岩石孔隙压力大于临界孔隙压力,从而诱发了地壳浅部岩石低初始应力的提前释放,而水的渗透作用则导致岩石孔隙压力随时间由高度非均匀分布向相对均匀分布转换.角07井地震序列的演化与这一物理过程有关.      

皎口水库地震精定位、震源机制及发震机理探讨

利用双差定位法对宁波皎口水库2009—2021年地震进行重新定位,通过HASH方法计算了较大地震的震源机制解,基于精确定位后的地震震源位置分布、深度分布以及总体震源机制解特征,结合区域地质构造背景、库区地质条件以及水库水位变化,对发震构造及孕震机理进行探讨。结果表明,皎口水库2009—2012年地震序列与2018—2021年地震序列孕震机理相似,区域性断裂长兴—奉化可能为发震断层;库水通过断裂破碎带下渗,孔隙水压力扩散及库水对岩石的软化弱化长时间作用,造成断层面及其深部某些薄弱部位变为不稳定区域,最后在库水变化导致的加卸荷应力场作用下触发地震。