用二相介质有限元方法对共和70级地震前兆的数值模拟-地震工程学报.PDFVIP

第 19 卷　第3 期 西　北　地　震　学　报 Vol.19　No.3 1997 年9 月 NORTHWESTERN SEISMOLOGICAL JOURNAL Sept.1997 用二相介质有限元方法对共和7.0 级 地震前兆的数值模拟 梁子斌　钱家栋 (国家地震局兰州地震研究所, 兰州　730000) 夏玉胜 (青海省地震局, 西宁　810001) 赵和云　蔡红卫　张慧 (国家地震局兰州地震研究所, 兰州　730000) 摘要　利用固-液二相介质平面应变问题的有 元程序对共和7.0 级地震的前兆 现象的时空分布进行了数值模拟.研究中不仅考虑了固相的非线性、硬化及膨胀等效 应, 也充分考虑了水对前兆的影响, 因此模拟结果与实际的前兆分布特征符合较好. 主题词:　地震前兆　数值模拟　有限元分析　共和地震 1　前言 共和地震 西北地区近20 年来发生的最大的一次地震, 破坏性也最大, 因而在西北地区 具有一定的代表性, 对该次地震进行深入研究 有一定意义的.共和地震发生在青海共和盆地 内, 其发震断层的走向为北15°西, 倾向北东, 倾角为75°～ 80°, 断层两侧显示出强烈的挤压现 象.震源机制解也表明该次地震以倾滑为主.据对航磁资料的分析, 在震中附近14 ～30 m 深 度上存在一盆形强磁性高密度体, 界面形状为南陡北缓, 这一块体对共和地震的孕育过程起着 重要作用. 对地壳这样复杂的多相介质, 不仅要考虑它的非线性力学性质(弹、塑性及膨胀等), 还要 考虑它内部多相物质共存及相互作用的影响.Rudni i and Rice 利用简单的微观模型给出了多 [ 1] 相介质的本构关系 , 这一模型不仅考虑了介质的非线性性质, 也同时考虑了流体的运移情 况, 因而更符合大地介质的真实情况.牛志仁等据此利用平面应力问题采用有限差分方法对唐 ① 山7.8 级地震过程进行了数值模拟, 模拟结果与实际地震前兆有较好的对应 .虽然共和地震 与唐山地震的地震序列均为主余震型, 但震源机制不同, 共和地震以倾滑逆冲为主.基于此, 本 文采用平面应变问题对共和地震的孕震过程进行数值模拟, 探讨共和地震的前兆在时空上的 分布规律. 地震科学联合基金资助课题(编号:95170). 现在国家地震局分析预报中心工作. 第一作者简介:梁子斌, 男, 1963 年8 月生, 助理研究员, 现从事地震预报地电方法理论研究工作. ① 牛志仁, 等.华北北部地震孕育及前兆机理的数值模拟研究.1995. 　　　　　　　　　　　　　　　　　　西　北　地　震　学　报　　　　　　　　　　　　　　　　　　第 19 卷 70 2 　地质模型、介质参数及边界条件的选取 地下流体基本包括水、油和气三种物质, 它们存在于岩石的孔隙中, 本文只考虑水的作用, 就 介质孔隙中只有水的二相介质情形.由多相介质的本构关系、介质连续性条件及达西定律 可直接得到介质刚度方程.本文拟采用有限元方法进行求解, 通过求解线性方程组, 从而得到 问题的解.具体实现参见附录. 依据共和地震的震源机制解, 主断层 走向为NW 、倾向为SW 的逆冲断层, 其主压应力 方向为NE, 而且附近的几条较大断层或破裂带也呈NW 向展布, 与发震断层近于平行, 参见 图1.因此本文选取沿NE 向通过震源的一条总长为300 m 的剖面作为研究区域, 如图1 中的 虚线AB 所示, 其深度取莫霍面以上, 从SW 到NE 剖面深度由深变浅.如果不考虑与剖面垂直 方向(NW 向)上的变化, 则可视为一个平面应变问题.将研究区域分为5 种不同介