介质黏滞性质对余震活动影响的数值模拟-地震地质.PDF

第 ３７卷　第 １期 地　震　地　质 Ｖｏｌ．３７，Ｎｏ．１ ２０１５年 ３月 ＳＥＩＳＭＯＬＯＧＹＡＮＤＧＥＯＬＯＧＹ Ｍａｒ．，２０１５ － ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５３ ４９６７．２０１５．０１．００５ 介质黏滞性质对余震活动影响的数值模拟 １，２） ３） ３） ４） 曲均浩 　蒋海昆 　宋　金 　李　金 １）中国地震局地质研究所，北京　１０００２９ ２）山东省地震局，济南　２５００１４ ３）中国地震台网中心，北京　１０００４５ ４）新疆维吾尔自治区地震局，乌鲁木齐　８３００１１ 摘　要　　假定主震破裂面上残余凹凸体的破裂导致余震，残余凹凸体尺度符合分形分布，在破裂 面上随机分布，残余凹凸体破裂／失稳强度符合对数正态分布；以岩石圈下部粘弹介质震后应力松弛 作为导致残余凹凸体破裂的动力输入，加载应力按负指数形式衰减；以库仑破裂准则作为失稳判据， 考虑残余凹凸体破裂对附近区域的应力加载，在不同条件下模拟生成包含时、空、强三要素的人工 － “余震”序列。通过 ＧＲ关系及修改的大森公式，检测模拟输出与实际余震活动统计特征的吻合程 度，以此作为进一步调整模型参数的依据。在此基础上重点讨论介质黏滞性质对余震活动的影响。 结果显示：表征岩石物理特性的黏滞系数对余震活动持续时间、衰减快慢起重要的控制作用，黏滞 － 系数与衰减系数之间显示负指数的相关关系。余震序列 ＧＲ关系与黏滞系数关系不大，主要受控于 残余凹凸体的尺度分布，即与构造不均匀性关系更为密切。 关键词　　主震破裂面　余震活动　残余凹凸体　应力松弛　黏滞系数 ＋ － － － 中图分类号：Ｐ３１５．３３ 文献标识码：Ａ 文章编号２０１５ ０　引言 大地震发生后，在震源区及邻近区域会观测到大量的余震活动。余震活动满足 ２个重要 － 的统计关系，一个是描述地震震级、频度之间关系的 ＧＲ关系，另一个是描述余震随时间衰 减特征的修改的大森公式。介质状况和应力分布的不均匀性及其相互作用，是影响余震活动 的重要因素（Ｍｏｇｉ，１９６２ａ，ｂ；Ｓｃｈｏｌｚ，１９６８ａ，ｂ，ｃ），以此为基础，许多研究从不同角度来解释 余震活动机理：１）从介质不均匀性的角度，普遍接受的事实是，余震是残余凹凸体的进一步破 裂（Ｋａｎａｍｏｒｉ，１９８１；Ｋｉｓｓｌｉｎｇｅｒｅｔａｌ．，１９９１；Ｓｏｍｅｒｖｉｌｌｅｅｔａｌ．，１９９９）；２）从应力变化的角度，速 － 率 状态依从、亚临界裂纹扩展、蠕变或余滑模型等认为主震造成的断层应力变化是余震发生 的主要原因（Ｄｉｅｔｅｒｉｃｈ，１９９４；Ｈｕｃｅｔａｌ．，２００３；Ｋａｎａｍｏｒｉｅｔａｌ．，２００４），但也有研究者通过对实 际余震观测资料的研究，认为主震造成的应力变化并非是控制余震发生的主要原因（Ｇａｓｐｅｒｉｎｉ ｅｔａｌ．，１９８９；Ｈｅｌｍｓｔｅｔｔｅｒｅｔａｌ．，２００３），或者此种控制作用非常微弱（Ｐｅｎｇｅｔａｌ．，２００９）；３）主震 断层滑动引起的孔隙流体扩散，也被认为是余震发生的重要诱因，余震频次与孔隙压变化之间 － － － － 〔收稿日期〕稿回。 〔基金项目〕　山东省自然科学基金（ＺＲ２０１４ＤＱ０１９）和山东省地震局青年基金（ＪＪ１３０８Ｙ）共同资助。