近场强地震动数值模拟的简化计算方法 - 地震地质.pdfVIP

第２８卷　第４期 地　震　地　质 Ｖｏｌ．２８，Ｎｏ．４ ２００６年 １２月 ＳＥＩＳＭＯＬＯＧＹＡＮＤＧＥＯＬＯＧＹ Ｄｅｃ．，２００６ 近场强地震动数值模拟的简化计算方法 １，２） ３） ４） 张冬丽 　陶夏新 　周正华 １）中国地震局地质研究所，北京　１０００２９ ２）中国地震局地震预测研究所兰州科技创新基地，兰州　７３００００ ３）哈尔滨工业大学，哈尔滨　１５００８３ ４）中国地震局工程力学研究所，哈尔滨　１５００８３ 摘　要　　近场强地震动除受场地条件的影响外，还受到震源破裂面上子源的空间分布特点、子源破 裂先后顺序的强烈控制，基于数值格林函数法的近场强地震动数值模拟方法可以综合考虑震源、传播 途径及局部场地条件的影响，对计算过程进行合理简化，分２步完成地震动模拟：第 １步，在介质均匀 区采用矩张量的解析解计算所有子源在盖层底面的位移，形成下一步有限元计算的输入场；第 ２步， 在盖层介质不均匀区，结合局部人工透射边界技术，采用时、空解耦的波动显式有限元方法计算地表 强地震动。在有限断层模型中，采用具有９个力偶的等效地震矩张量表达断层产状、滑动方向等的影 响，采用 Ｂｒｕｎｅ模型定义各子源的滑动时间函数，描述滑动的时、空不均匀分布特征，从而细化震源模 型。通过对 Ｎｏｒｔｈｒｉｄｇｅ地震中４个基岩台站地表地震动的模拟结果和强震记录，验证了此简化计算方 法的可行性。 关键词　　近场强地震动　数值模拟　数值格林函数法　Ｎｏｒｔｈｒｉｄｇｅ地震 中图分类号：Ｐ３１５．８ 文献标识码：Ａ 文章编号：０２５３－４９６７（２００６）０４－０６１２－１１ ０　引言 地震动是由震源、地震波从震源至场地的传播途径和局部场地条件三者控制的。强地震动 数值模拟过程则可以归纳为震源模型的建立、地震波在地壳介质中的传播和局部场地条件的影 响３个环节。基于这３个环节的近场强地震动研究一直是工程地震学的研究核心，数值模拟作 为近场强地震动研究的重要手段，近年来随着强震记录的丰富和迅速发展，大区域、复杂介质及 地形的三维数值模拟和预测研究已在国内外广泛展开，其模拟方法和预测成果在工程及地震危 害性评价中得到了广泛应用（Ｇｒａｖｅｓ，２０００；Ｔａｏｅｔａｌ．，２００２）。目前，正在开展的大城市活断层 探测与地震危险性评价项目中的地震危害性评价是项目可否达到预期目标的重要因素，也是为 城市建设和发展提供基础保证的重点所在。但是，由于近场强地震模拟耗费的计算机资源较 大，往往只是在大型计算机或并行计算机上才能实现，且计算过程中所需要的基础参数也较为 复杂。因此，本文在数值格林函数方法的基础上，采用有限元法，对计算方法进行了合理简化， 为在单机上实现近场强地震动数值模拟的全部过程提供了理论依据和技术方法。 〔收稿日期〕　２００６－０５－２５收稿，２００６－０８－２１改回。 〔基金项目〕　国家发展与改革委员会发改投资“城市活断层试验探测”项 目、国家自然科学基金 （ 、地震科学联合基金（Ａ０７０２７）、中国地震局“十五”重大建设项目（Ｊ１０５Ｂ１６）共同资助。 ４期 张冬丽等：近场强地震动数值模拟的简化计算方法 ６１３ １　震源细节是影响近场强地震动特征的重要因素 在震源、传播途径和局部场地条件３个环节中，近场地震动受震源的影响最为显著，而后者 是目前人类了解得最少的。震源的研究主要依靠地震学方法，震源反演技术在２０世纪８０年代 有了迅速发展（姚振兴等，１９９４；Ｗａｌｄ，１９９６）。由于格林函数在高频段的估计精度不够和高频 振荡波形反演等方面的困难，反演结果仅在频率低于 １Ｈｚ的范围内取得了很好的效果，在大多 数工程结构十分关注的高于 １Ｈｚ的频段内效果却不好。随后，在有条件的地区借助经验格林函 数法（亦称小震合大震）改善了高