三维数值模拟在华北地区现今构造变形分析中的应用.pdfVIP

第 36卷 第 1期 地 震 地 质 Vo1．36．No．1 2014年 3月 Mar．．2014 SEISM0L0GY AND GE0L0GY doi：10．3969／j．issn．0253—4967．2014．O1．012 三维数值模拟在华北地区 现今构造变形分析中的应用 胡勐乾 邓 志辉 陆远 忠 宋 键 路 雨u 朱秀云 ¨ 孙 锋 ¨ 1)环境保护部核与辐射安全中心，北京 100082 2)中国地震局地质研究所 ，北京 100029 3)中国地震局地壳应力研究所 ，北京 100085 4)青 岛市地震局 ．青岛 266034 摘 要 基于并行版 ANSYS数值模拟软件平台，根据华北地 区活动地块划分及活动断裂分布 ，结 合 GPS资料等 ，确定模型的几何边界范围是 99．8。一121．4。E，27．9。～42．3。N，模型的边界范围包括 华北活动地块 区的绝大部分以及周边青藏、西域、南华和东北亚地块 的部分地 区。模型单元平均边 长为25km，模型共划分为单元416582个 ，节点 582392个 。主要研究成果为 ：1)地壳运动速度模拟 与分析对 比：华北地 区地壳运动速率的整体分布呈现由东向西逐渐减小、自北 向南逐渐增大的特点， 数值模拟结果与GPS观测运动速度大部分较吻合；2)断层运动性质模拟与分析对 比：华北地区有限 元模型考虑 了华北地区已知所有晚更新世 以来 的活动断裂 ，模拟结果所得的断层运动性质与活动断 裂地质调查的运动性质符合得较好 ；3)应变场模拟与分析对 比：通过数值模拟，计算了1999--2004 年和 2004--2007年华北地 区地表第 1主应变和第 3主应变 (最大压应变 )的大小和方向，与前人利 用震源机制结果、GPS测量结果等资料反演得到的华北地 区主应变场较一致 。 关键词 有限元 数值模拟 活动断裂 华北地区 中图分类号 ：P315．2 文献标识码 ：A 文章编号 ：0253-4967(2014)0l_0148—18 O 引言 华北地区地质构造演化过程复杂 、构造变形独特 、地震活动性强 ，历来是众多学者研究的 焦点。不论在地质资料、活动断层的研究 ，或者是地球物理资料 、GPS观测的密度和时间上， 华北地区都 已有众多研究成果 (郭慧等 ，2011)。前人对华北地区构造变形的数值模拟研究取 得 了丰富的成果(邓志辉等，2011a，b)，但也存在一些科学问题需要进一步深入研究。 在数值模拟实验中，模拟结果的对 比是工作的重要环节之一 ，其用来检验数值模型的合理 性及数据的可靠性 ，这在很大程度上决定了数值模拟结果的可信度。遗憾 的是 ，这个环节常常 被人们所忽视，或者仅进行某一方面的结果对 比就认定数值模型的可靠性。 本研究在综合分析华北地区地质构造、地球物理资料 的基础上 ，利用高性能并行计算技 [收稿 日期] 2013-06-02收稿 ，2013-08—25改回。 [基金项 目] 国家科技重大专项(CAP1400)安全评审技术及独立验证子课题 (2010ZX060021016)资助。 1期 胡勐乾等 ：三维数值模拟在华北地 区现今构造变形分析rI1的心 f『J 149 术 ，建立华北地区的三维有限元模型，以高精度的GPS观测研究结果为运动边界 ，对华北地区 的运动速度场 、断层运动性质 、应变场三方面进行数值模拟结果对 比，从而验证华北地 区数值 模型的可靠性与合理性 。 l 华北地区三维有限元模型及边界条件 华北地区(106。～124。E，31。～42。N)，位于燕 山、阴山以南 ，秦岭、大别山以北，贺兰一六 盘断裂带以东一直到