北京地区地壳三维结构的综合建模-震灾防御技术.pdf

第 4 卷 第 1 期 震灾防御技术 Vol. 4, No. 1 2009 年 3 月 Technology for Earthquake Disaster Prevention Mar., 2009 谭毅培，王伟君，齐诚，陈棋福，2009. 北京地区地壳三维结构的综合建模. 震灾防御技术，4 （1）：24—39. 北京地区地壳三维结构的综合建模1 谭毅培 王伟君 齐 诚 陈棋福 （中国地震局地震预测研究所，北京 100036 ） 摘要 综合基础地理信息数据、数字化地质图、沉积层的地脉动探测结果、人工地震测深和宽角 折射等多种资料，在对数据进行分析和整理的基础上，利用建模工具软件 GoCAD，建立了北京 地区（东经 115.5°—117.5°，北纬 39°—41° ）地壳及上地幔顶部 P 波速度结构和地质界面模型。 在模型建立的过程中，分析得到了人工地震剖面数据的地学统计特性，并设计了一种适合离散光 滑插值（DSI ）约束点权重的光滑估计方差，对模型的插值结果进行了质量评估。分析了北京地 区地壳上地幔速度结构的特征，又将速度结构模型与层析成像结果进行对比，模型比较好地反映 了北京地区的地下结构。 关键词：北京地区 地壳上地幔顶部 三维 P 波速度模型 引言 探测获取地下精细结构信息是认识和了解我们人类赖以生存的地球，有效抗御包括地震 在内的自然灾害，开发利用地下油气矿产资源，保持社会经济持续发展的重要前提。现代信 息计算技术发展带来的高速计算能力，为基于强地震动数值模拟计算的地震危险性分析和抗 震设防规划，以及利用理论地震图和观测地震波形的差异来反演修正已有的地下结构模型， 提供了强大的技术支撑，同时也对综合应用已有的各种地球物理探测研究成果，相互补充、 解决冲突建立一个高分辨率的可信的地下模型提出了迫切的需求。 环绕北京的首都圈地区是我国的政治、经济和文化中心，也是我国地球物理勘探和地质 调查研究成果最丰富的地区之一。30 多年来，在以首都圈为中心的华北地区，不同学者利用 人工地震测深资料（刘昌铨等，1982；嘉世旭等，1995；2005 ；李松林等，2001 ；张成科等， 2002 ；段永红等，2002 ；魏文博等，2007 ）、地震层析成像（朱露培等，1990；孙若昧等，1995； 1996；张先康等，1996；于湘伟等，2003 ；王夫运等，2005 ；黄金莉等，2005 ；齐诚等，2006 ； 郑建常等，2008 ）等方法得到了各种不同区域，不同分辨率的速度结构模型，并应用于震源 的精确定位（刘福田等，1986），地壳形成和演化机制（张先康等，1996），强地震动模拟（高 1 基金项目 中国地震局地震预测研究所基本科研业务专项（编号：0207690229) 、国家国际科技合作项目（编号： 2006DFA21650 ）共同资助 [收稿日期] 2008-11-09 [作者简介] 谭毅培，男，生于 1982 年。现为中国地震局地震预测研究所研究生。主要研究方向为地壳结构和地震动模拟。 E-mial:oivertan921@ 1 期 谭毅培等：北京地区地壳三维结构的综合建模 25 孟潭等，2002 ；潘波等，2006 ；刘博研等，2007 ）等研究中。在此基础上，朱介寿等（1997） 对人工地震测深及天然地震面波体波三维层析反演数据进行统一处理，建立了中国及其邻区 大范围地球三维结构初始模型。张亚明等（2008 ）综合京津唐及邻区的８条宽角反射地震剖 面地壳速度结构，计算了该区域不同构造域上、中、下地壳平均速度和厚度。但由于人工勘 探的高成本及探测区域有限，而天然地震层析成像反演得到的速度结构分辨率低且明显依赖 于初始模型，以及地质调查、地球物理勘探和反演方法的不同与实施目的和时间的差异，至 今仍没有类似洛杉矶（Magistrale 等，1996；2000 ；Kohler 等，2003 ）那样建立一个基于多 种资料且得到普遍认可的首都圈地区三维模型。 为此，本研究试图综合基础地理信息、数字化地质图、地脉动的沉积探测