从库仑破裂应力和余震分布角度探讨汶川地震和芦山地震的关系

本文以龙门山及周边地区为研究对象,考虑区域地质构造差异、主要活动断裂带、地表附加重力影响,建立能反映地表起伏和岩石圈分层结构的龙门山地区三维粘弹性有限元模型。以地壳水平运动速率观测值为约束条件重建研究区现今构造背景应力场,在此基础上分别模拟了汶川地震和芦山地震的发生机理。通过分析同震库仑破裂应力变化与余震空间分布的关系,探讨了2次地震主震对余震的触发作用以及汶川地震对芦山地震的影响。研究表明,汶川地震和芦山地震的余震大部分由其主震触发,汶川地震对芦山地震的余震有约6.78%的触发作用。汶川地震的同震库仑破裂应力在芦山地震主震位置的增加值约为0.016MPa,如果龙门山断裂带南段库仑破裂应力年累积速率按照0.4×10-3-0.6×10-3MPa·a-1计算,汶川地震使芦山地震提前了约27-40年。计算还表明汶川地震和芦山地震的发生使鲜水河断裂带南段和虎牙断裂的库仑破裂应力增加,这些断裂带在未来发生地震的可能性增加。     

地面VLF波穿透电离层的能量衰减变化

地面VLF人工源可以导致辐射带高能粒子沉降.为研究辐射带粒子加速和沉降机制,乃至实现对辐射带电子"人工控制"的设想,需要在精确计算得到VLF人工源在电离层激发电磁场能量分布的基础上,计算哨声波对辐射带电子的调制作用.以往计算VLF人工源在电离层中激发的电磁场能量分布多基于Cary给出的波导中VLF波衰减率和Helliwell提出的经典电离层吸收曲线,但近期研究表明,这些模型结果存在较为明显的误差.本研究建立了地面VLF信号穿透电离层传播的全波计算模型,将计算结果与DEMETER卫星记录的NWC通讯台激发的电磁响应进行对比.虽然模型没有考虑电离层参数水平不均匀性,但模型计算结果与卫星观测值也有较好的一致性.利用经过验证的全波模型计算了不同地磁参数、电离层参数情况下,不同辐射特性的地面VLF辐射在波导中的衰减和穿透电离层时D/E区的吸收值,探讨了上述参数对电磁波能量在波导和D/E区中衰减的影响规律.     

2008年5月12日汶川MS8.0地震动力学背景的数值模拟实验

2008年5月12日汶川M_S8.0地震发生在龙门山断裂带.本文基于龙门山断裂带的地质与地球物理研究结果,以及高精度的地形数据、大地热流测量数据,建立了以龙门山断裂带为主要研究对象的有限元模型;以GPS观测数据、构造应力场和震源破裂过程研究结果为约束,研究了此次强烈地震的动力学背景.模拟实验结果显示:在考虑青藏高原物质向东挤压流动的同时,青藏高原与四川盆地的地形差异和流变强度差异、断层摩擦强度差异和断层产状形式均对地震起始破裂的发生位置和断层错动形式有着重要影响.本文利用地球动力学的有限元软件模拟了汶川地震地表破裂在龙门山断裂带上传播的过程.     

Faults有限元程序在Windows平台上的移植

Faults有限元软件综合高程、热流、断层等多个地球物理场参数,相较于Ansys、Adina等商用有限元软件能更好的处理断层的摩擦机制,因此Faults有限元软件在地学模拟中有独特优势.目前使用的Faults软件基于IBM科学与工程数学库(ESSL),使用环境较为苛刻,使用不便,需要移植源代码.利用Linpack免费数...     

无人机快速获取地震灾情的应用探索

快速获取地震灾情对于应急救援决策和指挥是非常重要的。大地震造成通讯和交通设施破坏,使应急救援指挥部不能及时全面了解灾情,汶川地震已经给我们惨痛的教训。为此,我们开展了用无人机快速获取地震灾情的探索。     

2019年6月17日四川长宁6.0级地震震后应力演化与余震关系

2019年6月17日四川长宁县发生6.0级地震,该次地震余震活动频度高、强度大,其中超过5.0级的强余震就有4次,具有不同于以往6.0级地震的独特特征.余震活动与震后区域应力变化密切相关,为研究它们之间的关系,考虑区域主要活动构造、地表起伏和深部反演结果,建立长宁地区岩石圈三维黏弹性有限元模型.采用数值方法重建基本符合研究区GPS观测和最大水平主压应力方向测量结果的现今构造背景应力场.进而依次模拟了长宁6.0级地震和5.0级以上强余震序列.通过计算库仑破裂应力变化研究了震后应力演化与余震分布,以及主震和5.0级强余震序列之间的关系.研究表明,长宁6.0级地震的发生可能与区域内非构造加载因素有关,余震活动明显受震后区域应力变化的控制.长宁地震后,于滩-长宁背斜在10 km深度应力得到充分释放,库仑破裂应力明显减小;而在3 km深度库仑破裂应力明显增加,应力水平仍然较高.      

汶川Ms8.0地震动力学背景的数值模拟实验

2008年5月12日汶川Ms8.0地震发生在龙门山断裂带.地球物理研究者利用多种方法对此次地震的动力学成因及过程进行了研究,本文从研究断层错动行为的角度出发来探索影响汶川地震断层初始破裂及破裂的传播过程的动力学因素.     

巴颜喀拉块体东南缘应力演化及鲁甸地震自发破裂模拟研究

巴颜喀拉块体活动导致的昆仑—汶川地震系列,是目前中国大陆唯一一组7级以上地震序列,未来几年该块体仍是发生7级地震的主体地区。巴颜喀拉块体东缘从走滑运动向逆冲运动转换的构造动力学性质使得活动块体两个相邻边界带存在大地震序列的关联性。研究表明近期逐渐加速的块体东边界强震序列是北边界强震活动的响应,巴颜喀拉块体南边界的鲜水河断裂带同样以左旋走滑的方式向东边界挤压,但巴颜喀拉块体东南缘地震呈现以下两个特点,一是强震大都发生在巴颜喀拉块体东边界断裂带,东边界强震序列呈现加速状态;二是南边界鲜水河断裂带在1973年炉霍7.6级地震后,已有40多年没有发生7级以上地震,鲜有6.5级以上强震,在此次地震序列中表现出一种大震缺失的状态。因此,巴颜喀拉块体东南缘主要断裂带的地震危险性值得关注和研究。同震库仑破裂应力及震后粘弹性松弛在解释余震序列分布,分析强震序列相互作用关系及未来地震危险性分析等方面发挥了重要的作用。比如汶川地震之后,一些研究者利用弹性位错模型计算了汶川地震的同震库仑破裂应力变化,并对龙门山区域未来的地震危险性进行了估计。一些研究者分别计算了汶川地震和芦山地震后周边断层的同震库仑应力变化,以及汶川地震和芦山地震引起的同震和震后粘弹性松弛应力场变化,得到了汶川地震可能有效地促进了芦山地震的发生的结论。在应力演化过程中,震间构造应力加载对强震危险性的作用也非常重要,目前比较流行的计算这部分应力积累的做法是采用"负位错"的理论进行计算,这一方法虽然简洁,但假设了断层是完全闭锁的,可能会造成一定误差。本文拟以巴颜喀拉块体区域构造活动作为三维有限元模型的边界约束,综合研究区域介质速度、密度结构最新反演结果计算介质的弹性材料参数,从初始应力场出发模拟在构造应力和重力作用下的现今背景应力场,计算强震序列引起的区域内断裂带震前、震间、震后的应力演化,获得更为准确可靠的应力场结果。利用速率-状态关系计算库仑应力变化对研究区域背景地震发生率的影响,利用泊松概率过程计算研究区域在汶川地震和芦山地震震后未来10年6级以上强震的发震概率及研究区域主要断裂带的6级以上强震发震概率。除以上述准动力学的方法研究巴颜喀拉块体东南缘应力演化及地震危险性,本文还以鲁甸地震为例,研究鲁甸地震从成核区开始破裂,然后扩展直至终止的自发破裂动力学过程,探索影响地震自发破裂位错的主要因素。研究内容:(1)建立巴颜喀拉块体东南缘三维粘弹性有限元模型。(2)模拟巴颜喀拉块体东南缘主要断裂带上1900年以来7级以上强震序列的应力演化过程,并对强震序列对芦山地震的触发作用进行了研究,对比库仑应力场及等效应力场演化结果与背景地震活动的关系。(3)应用速率-状态关系研究库仑破裂应力对汶川地震、芦山地震震后研究区域的地震发生率的影响,应用泊松过程给出研究区域主要断层未来的强震发生概率。(4)应用可以描述断层摩擦机制的FAULTS有限元软件对芦山地震断层错动方式,初始破裂点进行了模拟,研究了汶川地震后,断层强度的变化对芦山地震的影响。(5)应用曲线有限差分方法建立鲁甸地震非平面断层自发破裂模型,模拟鲁甸地震在包谷垴—小河断裂上的自发破裂过程,分析影响鲁甸地震自发破裂过程的主要因素。本文的创新点:(1)以三维速度、密度反演结果建立能反映更真实地下介质的三维粘弹性有限元模型;(2)从初始应力场出发,进行数十万年的应力场演化模拟直至获得与现今构造应力场更为吻合的背景应力场用于地震序列模拟;(3)采用泊松概率模型计算地震发生概率时直接基于三维粘弹性有限元模型得到的库仑应力结果,充分考虑了同震、震间和震后效应的影响。同时,在考虑背景地震发生率时采用了更全面完整的地震目录,这可能使得对未来地震危险性的估计更准;(4)对鲁甸地震自发破裂过程进行了动力学模拟,将震间应力积累过程与同震地震破裂过程的模拟结合起来。     

太原盆地应力场演化与强震关系的数值模拟研究

近年来,太原盆地内尽管没有强震发生,小震却异常活跃,而研究太原盆地构造应力场的分布和演化规律的工作比较少.采用有限元数值方法,综合考虑区域地质构造差异、主要活动断裂带、地形起伏和岩石圈分层结构,引入深部速度结构,建立包含太原盆地在内的山西地区三维粘弹性模型.以GPS观测值和最大主压应力方向测量值为约束条件,重建研究区现今构造背景应力场.在此基础上,依次模拟了自公元512年以来太原盆地6级以上和山西地震带7级以上历史强震序列.计算结果显示,太原盆地内的3次6级以上历史地震均落在应力场增加值大于0.01 MPa的区域,盆地近年来小震空间分布与现今应力增加值大于0.01 MPa的区域有较好的对应关系.研究表明:山西地震带内历史强震序列和长期构造加载对太原盆地内3次历史强震均有促进作用,太原盆地的地震活动性明显受区域当前应力水平的控制.      

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

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