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

利用机器学习方法对地震活动大数据进行挖掘,识别出一些过去认识不到的异常,提高地震预测的准确性,是一个非常具有挑战性的科学问题.本文基于川滇部分地区(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级地震的潜在危险性.     

利用GRACE观测数据研究苏门答腊区域的黏滞性结构

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

高分辨瞬时振幅和瞬时相位在薄互煤层识别中的应用

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

河北文安5.1级地震前后重力场动态变化

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

虚拟川滇——基于千万网格并行有限元计算的区域强震演化过程数值模型设计和构建

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

汶川大地震地震波传播的谱元法数值模拟研究

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

华北地区水位下降是否会减缓气温上升——浅部地温影响的数值模拟分析

华北地区由于长期持续的地下水过量开采,导致了大面积地下水位大幅下降,引发地面塌陷、地下水质污染等一系列地质环境问题,这些现象早已为人们所熟知和关注.然而地下水位下降还会造成百米量级浅部地温及其梯度的变化,因此即使来自地球深部的大地热流密度没有变化,年度平均的从表浅部位通过地表实际传导进入大气的热流密度会减小,这是中外文献中尚未见讨论过的问题.我们通过数值模拟发现假定大地热流密度不变的条件下,华北数万平方公里地下水位下降会造成百米尺度内的地温降低,从而传入大气的热流密度降低40%以上,且会持续数百年以上的时间.这种长时间大范围的传导入大气的热流密度变化对环境会造成什么影响是一个十分值得关注的问题.这一预测在一定程度上得到了气象站地温观测数据的支持,但由于目前气象观测站只有3.2m深度范围内的地温资料,累计不超过5、60年,中间还有10余年的间断,而且表浅深度地温受地表多种因素的影响也较大,这些资料难以对我们关心的地下水位下降引起流入大气的热流密度变化这一问题提供直接确凿的数据来进行分析,因此今后有必要开展对地下数十乃至数百米地温进行持续精确的监测工作.     

固体潮应变各向异性研究

钻孔应变仪的观测过程中发现潮汐应变响应存在方向各向异性, 推测是受到台站附近断层影响导致的。 我们对此展开数值模拟研究, 讨论软弱断层对潮汐变形传播的影响。 我们的研究结果证实了软弱断层的存在确实能够产生能够观测到的各向异性现象。 同时我们还发现其各向异性变化对断层弹性参数较为敏感。 观测到的部分地震震前反映固体潮应变各向异性的潮汐因子"玫瑰图"变化有可能是由于台站附近断层在震前的性质变化引起的。     

致密砂岩水压瞬态致裂液化格子波尔兹曼及有限元数值研究

水压瞬态致裂液化在地应力测量、 地震破坏评估和机理研究、 油气(天然气、 页岩气)及地热资源开发等地学领域都具有重要的理论和应用价值。 自20世纪60年代至今在理论和实际应用方面取得了一些奠基性和开创性成果, 但鉴于问题复杂性, 瞬态致裂液化机理至今尚不清楚。 本研究应用格子波尔兹曼及有限元多孔介质流固耦合物理模型, 对地震波载荷作用下致密砂岩水压瞬态致裂液化过程进行数值模拟研究。 首先, 以鄂尔多斯盆地某油田延长组致密砂岩为例, 利用X射线CT断层成像技术, 应用基于量子力学第一性原理格子波尔兹曼方法, 建立致密砂岩数字岩芯模型。 进而, 推导格子波尔兹曼及有限元多孔介质流固耦合数值模型公式, 建立致密砂岩水压致裂液化物理模型。 最后, 应用地震波载荷边界和初始条件, 模拟瞬态致裂液化流固耦合过程, 讨论了地震波载荷幅值、 频率及作用时间对致密砂岩孔隙结构(孔隙度大小及连通性)、 致密砂岩破裂最大主应力之间关系, 得到了地震波作用下致密砂岩致裂液化准则。     

先存组构对各向异性岩石流变强度的影响

地壳深部岩石普遍存在变形组构,花岗质岩石中的变形组构不仅影响岩石强度,而且对后期变形具有显著控制作用。近年来,先存组构对各向异性岩石的流变强度影响成为高温高压实验研究的热点之一。文中对前人给出的各向异性岩石(包括云母片岩-片麻岩、石英-钙长石均匀混合体与层状组构样品)半脆性-塑性流变实验数据进行了重新整理与分析,结合作者开展的不同组构条件下花岗片麻岩与糜棱岩流变实验结果,讨论了先存组构对各向异性岩石流变强度的影响。实验数据表明:1)各向异性岩石的面理与最大主应力方向的角度是影响强度的主要因素。在半脆性破裂域,样品压缩方向垂直于面理(PER)和平行于面理(PAR)的强度基本相同,在压缩方向与面理呈30°夹角时,岩石破裂强度最小;在塑性流变域,垂直于面理方向的强度显著高于平行于面理方向的强度,当面理与最大主应力方向的角度为45°时,岩石强度最小。2)后期变形对原有组构的继承与改造程度,决定了各向异性岩石强度高低。3)样品中矿物的含量、分布与粒度对各向异性岩石强度有显著影响。理论模型预测结果与云母片岩实验结果比较吻合,但其他类型各向异性岩石的流变比理论模型结果要复杂得多。因此,进一步开展具有先存组构的各向异性岩石的流变实验,并将实验变形与实际地质条件下更为复杂的岩石变形进行对比分析,是认识各向异性岩石流变和变形机制最有效的方法。