青藏高原及邻区未来地震活动性趋势数值分析          下载免费PDF全文

董培育  石耀霖  程惠红  乔学军 《地球物理学报》2020,63(3):1155-1169

区域中长期地震危险性数值分析研究,需要对其初始构造应力场有所了解,但目前以及未来一段时期内仍无法直接观测到深部孕震层区域的应力场状况.本文首先基于岩石库仑-摩尔破裂准则,利用青藏高原及邻区百年历史范围内的强震信息,来反演估算该区域的初始应力场.然后,考虑区域构造应力加载及强震造成的应力扰动共同作用,重现了历史强震的发展过程.然而对于初始应力场的反演估算,本文仅能给出区域其上下限的极限值,并不能唯一确定.因此,采用Monte Carlo随机法,进行大量独立的随机试验计算,生成数千种有差异的区域初始应力场模型,且保证每种模型都能令历史强震有序发生,但未来应力场演化过程不尽相同.最后,将数千种模型在未来时间段内的危险性预测结果集成为数理统计结果,据此给出了区域未来的地震危险性概率分布图.初步结果显示未来强震危险性概率较高地区集中在巴颜喀拉块体边界及鲜水河断裂带地区.  相似文献   

2001年昆仑山口西M_S8.1地震对周围断层的应力影响数值分析          下载免费PDF全文

任天翔  程惠红  张贝  石耀霖 《地球物理学报》2018,61(12):4838-4850

大地震的发生往往会引起周围区域形变场和应力场变化，且对临近断层上的应力状态也有影响.2001年11月4日，昆仑山口西发生了半个世纪以来中国最大的M_S8.1级地震.本文基于已有的滑动模型，建立了三维含地形高程的横向不均匀性椭球型地球有限元模型，采用等效体力方法，分析了此次M_S8.1地震产生的全球同震位移和应力场变化.与解析方法相比，该模型考虑了地形、Moho面起伏和地球介质横向不均匀性；与一般的有限元数值模拟相比，该模型考虑了地球曲率和椭率，合理地规避了有限块体模型假定边界位移为零所引入的误差.计算得出同震位移与GPS观测数据可以很好地吻合.据库仑破裂应力准则和震源参数，计算得出昆仑山口西M_S8.1地震的发生造成了汶川、芦山、改则和当雄地震的发震断层上库仑应力增加，对这些地震的发生起促进作用；而造成玉树和德令哈地震发震断层上的库仑应力变化为负值，在一定程度上抑制了这些断层的地震活动性.此外，计算结果显示地球地形高程、介质非均匀性和椭率对昆仑山口西M_S8.1地震同震变化计算有一定的影响，其中地形和椭率造成的同震位移场相对误差约10%.  相似文献   

历史大地震断层滑动模型的建立及其对同震数值计算的影响——以1920年宁夏海原M_S8.5大地震为例          下载免费PDF全文

邓园浩  程惠红  张贝  张怀  石耀霖 《地球物理学报》2018,61(3):975-987

大地震的发生会引起区域位移场和应力场发生变化，进而改变区域内及临近断层的应力状态和地震活动性.目前，研究学者可据已有的断层滑动模型来计算分析大地震同震应力变化，同时采用库仑应力触发理论来进一步分析震后余震分布和断层危险性.然而，历史上曾经发生过不少大地震，例如，1920年的海原M_S8.5大地震，是全球范围内少见的特大地震之一.局限于无确切的地震台站地震波等资料，前人在研究历史地震的影响时往往给出一些简单的断层滑动模型，将断层面上错动量视为均匀分布.为更准确地了解历史地震对后续地震的影响，基于前人研究和一般地震滑动形态分布规律及地震反射剖面等资料，以海原M_S8.5大地震为例，探讨了如何建立海原大地震断层滑动模型，并分别搭建了简单断层滑动模型和复杂断层滑动模型的全球同震横向不均匀并行椭球型地球模型.通过对海原M_S8.5地震的同震位移场和应力场的计算，发现采用复杂断层滑动模型比简单断层滑动模型地表位错分布更切合实际.同时，进一步计算和分析了此次大地震对青藏高原东北缘近100年历史地震和周围断层的应力触发作用，得出断层滑动模型对同震计算结果的影响集中在发震断层附近而对远场影响较小.  相似文献   

长短时记忆神经网络在中期地震预报中的探索--以川滇地区为例     

李林芳  石耀霖  程术 《地球物理学报》2022,65(1):12-25

利用机器学习方法对地震活动大数据进行挖掘,识别出一些过去认识不到的异常,提高地震预测的准确性,是一个非常具有挑战性的科学问题.本文基于川滇部分地区(24°N—32°N,98°E—106°E)地震目录,采用滑动的时空窗口,选取16个反映地震时空强度分布特征的地震预测因子,建立了长短时记忆(LSTM)神经网络,对研究区域9个子区块未来一年的最大地震震级进行预测.通过设置训练集∶测试集=8∶2和训练集∶测试集=7∶3两个不同的训练测试模型,对过往发生的地震进行了回溯性预报.结果表明:训练集∶测试集=7∶3模型能够利用1970年1月至2004年9月的地震目录进行学习,成功回溯性预报2008年汶川地震;训练集∶测试集=8∶2模型利用1970年1月至2009年5月资料进行训练,回溯性预报2010—2019年间6级以上地震的R评分为0.407,回溯性7级地震预报时准确率高达92.31%.本文还探讨性给出预测意见:2022年2月前研究区西部、中部、东部、西南部存在发生5.1~5.3级地震的潜在危险性.  相似文献   

巴颜喀拉块体东缘形变及九寨沟地震孕震环境数值分析   总被引：6，自引：3，他引：3       下载免费PDF全文

庞亚瑾  程惠红  张怀  石耀霖 《地球物理学报》2017,60(10):4046-4055

2017年8月8日四川发生九寨沟M7.0地震,是继2008年汶川M8.0地震后发生在巴颜喀拉块体东部的又一强震.现今GPS速度观测数据显示,2008年汶川地震前后的1999-2007年和2011-2016年两个时间段内巴颜喀拉块体东部地表速度场存在明显的差异.本文以实际GPS速度观测资料为约束,构建三维有限元地球动力学模型,分别计算分析了两个时段内震源区及周边现今地壳形变、弹性应变能和应力积累特征,进一步探讨汶川地震的发生对九寨沟地区变形及应力的影响.数值模拟结果显示,汶川地震之后(2011-2016年)巴颜喀拉块体东部的地壳形变、弹性应变能积累及应力积累速率均明显大于震前,增加量值达1.5-3倍;九寨沟地震发震断裂上库仑应力增长率在1999-2007年约为0.7 kPa·a~(-1),2011-2016年间增至1.2 kPa·a~(-1).上述结果表明,现今巴颜喀拉块体东部地壳应力积累过程有利于左旋走滑型九寨沟地震的发生,汶川地震的发生调整了区域应力状态,加速了九寨沟地震的孕育过程.  相似文献   

利用GRACE观测数据研究苏门答腊区域的黏滞性结构   总被引：1，自引：0，他引：1  

王武星  石耀霖  孙文科  张晶 《中国科学:地球科学》2011,(6):773-783

本文尝试利用卫星重力观测资料和震后黏弹性松弛理论研究苏门答腊地区的区域流变结构,为更好地认识区域地球动力学环境提供依据．利用GRACE卫星重力资料,计算了2004年苏门答腊Mw9．3地震的同震及震后的重力变化．计算中使用平滑半径为500km的高斯滤波器．结果显示苏门答腊Mw9．3地震破裂东侧以陆地为主的上盘同震下沉量很大,造成约9x10＾-8mS＾－2的重力下降阶变,而西侧处于海水下的下盘重力同震上升较小约2x10＾－8mS＾－2,但其震后上升较快．流变结构对岩石的变形有很大的影响,是地球动力学数值模拟取得可靠科学结果的基础．本文尝试了利用卫星重力变化时间序列来反演苏门答腊地区的黏滞性结构;即基于GRACE时变重力场,利用自重力、黏弹性、平面分层模型模拟了该地震的同震和震后变形,并将获得的空间固定点的重力变化与GRACE重力场及点位时间序列进行比较,估计该地区的黏滞性系数在1．0x10＾18PaS的量级,且断层两侧的流变参数存在差异．最后结合苏门答腊区域的构造特点讨论了黏滞性系数的影响因素．  相似文献   

高分辨瞬时振幅和瞬时相位在薄互煤层识别中的应用   总被引：2，自引：1，他引：1  

于捷  苗广文  冯晅  刘财  石耀霖  李鹏 《地球物理学进展》2011,26(4):1386-1392

东北地区煤田中的煤层一般是一个典型的低速薄层,可采煤厚约1～8 m,这使得常规地震方法勘探煤层受到很大的限制.但由于煤层与围岩组成的是一个有规律的薄互层组合,所以通过研究薄互层的地震波波形、振幅、能量分配关系等特征来进行识别成为地震勘探煤层的一种可能.而高分辨三瞬技术能够提供精细化的波形、振幅等信息,因此本文将高分辨三...  相似文献   

河北文安5.1级地震前后重力场动态变化   总被引：1，自引：0，他引：1  

王武星  石耀霖  张晶  史彦华 《地震》2009,29(2)

分析了2006年7月4日文安5.1级地震前后文霸流动重力测网的观测结果,讨论了垂直形变和浅层地下水位变化引起的重力改变,并对重力资料进行了高程改正.研究结果认为,文安地震前后,测网出现了明显的异常变化,测区的重力异常以重力上升为主,且上升幅度大的区域主要在测网南部;地震震中在震前形成重力高梯度带,2005年10月至2006年4月短期内的重力变化表现最为明显;测区内的正负重力变化的分界线,地震前后出现逐渐偏转的特点,且在2007年10月逐渐恢复.  相似文献   

虚拟川滇——基于千万网格并行有限元计算的区域强震演化过程数值模型设计和构建   总被引：1，自引：0，他引：1  

张怀  吴忠良  张东宁  刘杰  王辉  严珍珍  石耀霖 《中国科学D辑》2009,(3)

利用合理的数值模型,对区域性的强震活动性分析,是近些年来地震学研究领域的一个热点问题.但无论从时间还是空间角度来讲,区域性强震演化的理论与数值模拟模型都是多尺度问题,对该问题的充分研究需要借助于大规模并行数值模拟技术.侧重于基于千万有限元网格技术,采用大规模并行有限元计算手段,对该问题的理论建模与并行数值算法实现,以及如何应用到对川滇地区区域强震演化过程的模拟进行了具体论述.本研究的关键目标是为中长期地震危险性估计提供一个大规模数值预报模拟平台系统.  相似文献   

汶川大地震地震波传播的谱元法数值模拟研究   总被引：3，自引：0，他引：3  

严珍珍  张怀  杨长春  石耀霖 《中国科学D辑》2009,39(4):393-402

汶川大地震震中位于震害活跃的龙门山断裂带上,其强度超过1976年唐山大地震．地震激发的地震波可以用来研究震源破裂过程及地球内部的层次结构．从弹性波理论基础出发,采用AK135理论地球模型,考虑地表地形、地球介质衰减及地球椭率等特性,利用谱元法结合高性能并行计算,分别对人工点源和复合源所激发的汶川大地震地震波的全球传播过程进行数值模拟重现．数值模拟结果显示了地震波在地球表面的传播形态．通过比较发现,复合源数值模拟结果较点源数值模拟结果更能体现汶川地震震源破裂过程的时空特性．另外,将两种震源的数值模拟结果分别与实际观测台站记录波形资料进行拟合对比,进一步定性地认识到多个点源组成的复合源的数值模拟结果较单个点源的数值模拟结果与实际观测资料的拟合程度更好,证明了汶川大地震破裂过程为多阶段破裂组成的一个复合破裂过程．  相似文献