评“紫坪铺水库造成孔隙弹性耦合变化及其对2008年汶川地震触发作用”二维模拟的局限性

2014年10月《地球物理学报》出版的陶玮等(2014)"紫坪铺水库造成孔隙弹性耦合变化及其对2008年汶川地震触发作用"一文基于完全耦合孔隙弹性理论,采用二维有限单元法模拟了紫坪铺水库蓄水造成的库区应力场和孔隙压力场的变化.该文模拟计算原则和思路是正确的,但仍然采用二维模型是一个不足.真实三维模型中水库是有限的一定面积内的载荷;而二维计算模型,却把水库载荷作为无限长的带状载荷,夸大了水库的作用.如果考虑三维模型,该研究得到的震源处库仑应力增长2~25kPa的结论可能只能解读为0.7~8.0kPa.对于紫坪铺水库是否对汶川地震有触发作用,目前各个研究组估计的库仑应力变化恰好处在是否可以触发的边缘.取得进一步认识的关键有赖于对地下结构、状态和物性认识的突破性进展.     

红河断裂带库仑应力演化及未来地震危险性估计

现今地震活动性显示红河断裂带强震主要发生在断裂带的北段,中段长期存在地震空区.为了更好地理解红河断裂带不同段落地震活动性差异和评估未来潜在地震危险性,本研究基于分层半无限空间黏弹性地球模型计算了自1833年嵩明M8.0地震以来,红河断裂带上及其周边共25次强震由于同震应力阶变、震后黏滞松弛和震间构造应力加载的综合作用,红河断裂带上库仑应力变化的演化过程.结果表明,在近180年红河断裂带南北两段得到加载,中段始终位于应力影区,这可能部分解释红河断裂带南北段地震相对活跃、中段长期存在的地震平静现象.假设未来三十余年该地区不再发生大震,考虑震后和震间作用,红河断裂带中段应力影区仍然存在,但范围在缩小;洱源附近、大理至大斗门以北地区、元江以南地区应力增加超过0.1 MPa,可能仍是地震潜在危险区段.     

卷积神经网络余震自动识别及在中国地震台网地震检测中的应用

正随着地震网络在全球范围大规模部署,每天产生的大量连续波形记录对数据处理的新方法提出了要求。传统的波形拾取算法都有一定的局限性,很难同时适用于复杂区域的不同类型地震事件的自动识别与到时挑选。机器学习和深度学习神经网络方法由于具有较强的泛化能力,近些年来受到越来越多关注。本研究使用汶川地震的10 000多个余震和专家手动挑取的P、     

国医大师熊继柏辨治胃脘痛学术特点及组方配伍规律研究

目的：分析和研究国医大师熊继柏辨治胃脘痛的学术特点及组方配伍规律。方法：搜集整理学术论著、期刊论文、报纸文章中熊老治疗胃脘痛的方剂,分析其辨治特点,并利用中医传承计算平台进行用药频次分析、关联规则分析,利用复杂系统熵聚类方法分析其核心组合,形成新处方以研究其组方配伍规律。结果：共纳入处方72首,涉及中药100味,多为理气药、补虚药,总用药频次875次。中药药性以温、寒、平为主,五味以苦、辛、甘为主,多入脾、胃经,肺、肝经次之。药物频次、组合和关联规则分析显示多以甘草、陈皮、木香、砂仁、枳实、白芍、厚朴、半夏八味药的两味至五味配伍为主。结论：熊老治疗胃脘痛辨证首重分虚实,治疗注重行气止痛,处方以甘草、陈皮、木香、砂仁、枳实、白芍、厚朴、半夏为核心,通过不同配伍达到治疗目的。     

大数据时代的地球动力学数值计算方法回顾与展望

随着大数据时代的到来,计算地球动力学数值计算方法体系更加完善.本文系统地回顾了传统数值模拟方法在计算地球动力学领域的应用进展,包括:有限差分法、有限单元法、谱方法和谱元法;并对近年来一些新发展的算法和应用前景进行了综述,如:不连续Galerkin法、小波方法和格子玻尔兹曼方法等.本综述有助于读者以整体视角了解地球动力学数值计算方法的发展脉络,并对大数据时代下研究适应日益丰富的数据和新算法提供有益参考.     

贝叶斯同化重力反演方法构建龙门山地壳密度模型

重力异常对地壳横向密度变化敏感,而无约束重力反演得到的密度模型其垂向分辨能力往往不理想.为了改善反演结果的垂向分辨率,本文参考已有先验分层模型,基于贝叶斯原理,提出了一种重震联合反演的新策略,可实现多种参考模型和复杂加权参数条件下的最大后验概率估计.理论模型测试结果表明,对于深度加权、多参考模型约束等多种问题,本文提出的新方法都可以稳健地获得最优化的模型参数.本文同时以中国地震科学台阵在龙门山地区及周边的一维接收函数分层模型和地震层析成像结果为参考,通过此方法对该区的重力异常进行反演,获得了该区的高精度三维密度结构,其水平分辨率优于10 km,垂直分辨率优于5 km.结合四条通过汶川和芦山地震震中的剖面进行分析后发现,反演得到的密度结构模型在过强震震源区位置横向变形显著,其揭示的分层地壳结构和变形模式与地表已知断裂构造具有相关性.本文提出的重震联合反演新策略,可为研究潜在强震风险源区的地壳结构和物性特征提供有效的科技方法支撑.     

2015年智利Illapel M_W8.3地震同震效应及其对南美大陆地震危险性影响

基于高性能并行有限元方法,建立含地表地形和Moho面起伏的大规模非均匀椭球地球模型,计算了2015年9月17日智利Illapel M_W8.3地震同震效应,同时讨论2010年智利Maule M_W8.8地震和2014年智利Iquique M_W8.1地震对2015年智利Illapel地震的加载作用,进一步分析这3个地震共同作用对南美大陆及周围断层的地震危险性影响.结果表明:对此次智利Illapel地震CEA(中国地震局地球物理研究所)滑动模型计算的地表最大水平位移约为3 m,USGS(美国地调局)模型的最大水平位移约6 m.特大地震影响范围广,南美大陆接近一半区域同震水平位移量级达到0.5 mm;受3次地震共同作用,智利2010年Maule地震南部俯冲带、2015年Illapel地震和2014年Iquique地震之间俯冲带,以及2014年Iquique地震以北俯冲带,受到10 kPa东西向压应力加载;2010年Maule地震和2014年Iquique地震共同加速了2015年Illapel地震滑动;1960年Valdivia M_W9.5地震破裂区、2014年Iquique地震和1995年Antofagasta地震之间区域以及2010年Maule地震和2015年Illapel地震的未破裂区重合处库仑应力变化均超过10 kPa,有可能加速强震的发生.     

2016年11月13日新西兰凯库拉(Kaikoura)MW7.8地震同震位移和应力数值分析

2016年11月13日新西兰南岛北端凯库拉（Kaikoura）发生了M_W7.8大地震,造成了强烈的地表变形并引发大面积滑坡和海啸的发生.基于美国地质调查局（USGS）断层滑动模型,建立全球同震横向不均匀并行椭球型地球模型,计算了此次新西兰凯库拉大地震产生的同震形变和应力及库仑应力变化.初步计算结果表明：新西兰凯库拉M_W7.8地震造成断层上盘东北向抬升,下盘西南俯冲;引起发震区域同震位移较大,从凯库拉到坎贝拉（Campbell）以及首都惠灵顿（Wellington）整体上东北向抬升,最大同震水平位移1.2 m,垂直位移1.1 m.此次大地震释放了发震断层上积累的压应力,但增加了发震断层两端的挤压力;同时,同震剪应力变化增加了NE-SW向断层发生右旋滑动的危险性;采用此次地震发震断层参数计算得出的最大库仑应力变化增加区域集中在发震断层两端,可达到MPa量级.当分别采用新西兰北岛Awatere断裂系和南岛Wellington断裂系参数计算库仑应力变化时,发现新西兰北岛和南岛震中以南区域的库仑应力均增加,可触发部分余震的发生.     

2010智利Maule特大地震的同震效应

2010年智利Maule大地震是人类有观测记录以来第六大地震,基于高性能并行有限元方法,建立含地表地形和Moho面起伏的大规模非均匀椭球地球模型,计算了该地震同震效应,并根据库仑应力变化分析周围断层地震活动性和主震对余震的触发关系.结果表明:对于2010年Maule地震,地球曲率和层状结构对同震水平位移影响不可忽略,如果不采用球面分层模型,在低纬度远场可能会造成位移和应变量值均被低估,在高纬度远场可能会造成位移被夸大而应变被低估.如果不考虑介质横向不均匀和地形,同震应力降被显著低估,近场同震东西向应力变化误差最高95.4%,南北向应力变化误差最高90.8%.Illapel地区库仑应力变化约10 kPa,相当于15年构造积累,Maule地震加速了2015年Illapel地震发生.以最优破裂方向投影得到库仑应力变化,经统计70.9%余震落在库仑应力变化量值大于10 kPa区域.本次特大地震触发了Pichilemu附近M_w6.9和M_w7.0两次正断型余震.