中国大陆岩石层动力学数值模型的边界条件 - 地震学报.PDF

第 卷 第 期 地 震 学 报 ， ２１ ２ Ｖｏｌ．２１ Ｎｏ．２ 　 　 　 　 　 　 年 月 （ ） ＡＣＴＡＳＥＩＳＭＯＬＯＧＩＣＡＳＩＮＩＣＡ ， １９９９ ３ １３３ １３９ Ｍａｒ． １９９９ 　 　 ～ 　 中国大陆岩石层动力学数值模型 的边界条件 张东宁 许忠淮 　　 （中国北京 １０００８１中国地震局地球物理研究所） 摘要 利用粘弹性蠕变本构关系的有限元方法，建立了中国大陆岩石层应力场及 　 应变场的动力学模型；以中国大陆岩石层应力场的水平向最大主压应力方向的 最新研究成果，以及块体水平运动的方向及速率为约束条件，用试错法得到了 模型边界作用力的大小；同时分析了青藏高原重力势能对中国大陆岩石层运动 格局的影响．数值模拟试验结果显示，印度洋板块对欧亚大陆的挤压作用及岩 石层底部的突然剥离，是导致青藏高原目前快速抬升及正断层发育的动力源． 青藏高原物质的水平侧向扩展是形成中国大陆现代构造特征的主要动力因素． 模型东部所受到的边界挤压作用，主要来自太平洋板块的俯冲作用，菲律宾板 块对欧亚板块的挤压作用可能非常微弱． 关键词 中国大陆 应力场 位移速率 边界条件 有限元 　　 　 　 　 　 引言 一些学者认为，地壳应力、应变场的动态演化和迁移过程，导致了强震活动的成组性 （张国民，傅征祥， ； ， ， ） 利用数值模拟方法，可以体现出地壳中各主 １９８５ Ｋｅｒｒ １９９５ １９９６． 要因素的变化对应力应变场不均匀分布动态图象的影响，进而可以预测强震危险区 目前 ． 已经开展了探讨构造应力场演化过程与强震活动之间关系的三维数值模拟工作（梅世蓉， 梁北援， ；张东宁，曾融生， ；张东宁，许忠淮， ） 这些模拟工作考虑了活动 １９８９ １９９５ １９９５． 断裂带的几何与力学特征，利用地震学、地球物理学和大地测量学等观测研究结果（如地 震勘探剖面和地震层析成像研究），确定地质构造的数值模型，模拟出不同尺度区域内应 力应变场动态演化图象 用计算应力变化来预测未来地震的方法很吸引人（ ， ）， ． Ｋｅｒｒ １９９５ 但利用数值模拟方法分析岩石层动力学过程的基本前提之一，是选取合适的边界条件，以 得到与实际情况基本符合的背景应力场及应变场． 以原地应力测量、地震震源机制解等资料显示的应力场方向为约束条件，张东宁和高 龙生（ ）利用有限元计算推断出，中国大陆地震活动发生在统一的区域性的主应力场背 １９８９ 景下，该应力场是印度洋板块和太平洋板块等共同作用于东亚大陆