利用烈度数据点估算中强历史地震震级与震中的方法-震灾防御技术.PDF

第 9 卷 第 1 期 震灾防御技术 Vol. 9 ，No. 1 2014 年 3 月 Technology for Earthquake Disaster Prevention Mar. ，2014 吴清，高孟潭，2014．利用烈度数据点估算中强历史地震震级与震中的方法研究. 震灾防御技术，9 （1）：12—28. 利用烈度数据点估算中强历史地震震级与 震中的方法研究1 吴 清 高孟潭 （中国地震局地球物理研究所，北京 100081 ） 摘要 本文基于地震烈度数据点，用椭圆参数方程对烈度点空间分布进行最小二乘拟合，得到各烈度 区椭圆烈度估计线，然后对烈度估计线进行统计回归，得到一组适用于青、甘、川、滇 4 省中强地震 的椭圆烈度分布模型I f (M ,R ) 。基于所建烈度分布模型，联立考虑中心点和方向性的椭圆数学方程， 代入全部烈度数据点估算地震震级和宏观震中。本文的试算震例表明了此方法的可行性，并进行了不 确定性分析。利用该方法对明清时期 4 次烈度点较少的中强历史地震参数进行了估算，所得结果表明 了此方法对历史地震有效。 关键词：烈度数据点 椭圆烈度分布模型 宏观震中 震级 历史地震 引言 历史地震是分析区域地震活动特征，预测地震活动趋势不可缺少的基础资料。工程抗震设计 工作中提出的地震危险性分析方法，其历史重演原则和地质构造类比原则也在很大程度上依赖于 历史地震。近年来，历史地震研究取得了丰硕的成果，但始终存在着史料不均衡、不完备，地震 参数误差较大等问题。即使相同的史料，不同专家的理解也可能产生分歧，不同版本地震目录所 确定的同一地震参数可以出现较大差异。而新一代地震区划图的编制，对历史地震参数精度、可 信度都有更高的要求。因此，有必要探讨估算历史地震参数的新思路和新方法。 地震烈度是指某一地区地面及各人工建筑物遭受一次地震影响的强烈程度（胡聿贤， 2006 ）。公元 1900 年以前，在没有仪器可以记录地震动各种物理参数的时候，地震烈度资料 就成为记录历史地震的唯一指标。历史文献中关于地震破坏的记录，以及震后烈度调查资料， 是某个具体地点的破坏程度或烈度记录，可以据其评定该地点的烈度值，称为烈度数据点。 烈度数据点从某种意义上说与地震台站有一定的相似处，仿照台站数据处理方法，用烈度数 据点对一些主要地震参数进行定量估计，理论上切实可行（谢毓寿，1991），更是现今研究历 1 基金项目 国家国际科技合作项目远东地区地磁场、重力场及深部构造观测与模型研究（2011DFB20210-04 ）与中国地震 局地球物理研究所基本科研业务专项“基于烈度数据点的历史强震（M ≥6.5 ）参数估计方法研究（023903 ）” 联合资助 ［收稿日期］ 2013-06-27 ［作者简介］ 吴清，女，生于 1985 年。中国地震局地球物理研究所在读博士生。主要从事历史地震与工程地震研究。E-mail ： wuqing908@ 1 期 吴清等：利用烈度数据点估算中强历史地震震级与震中的方法研究 13 史地震的重要手段。 基于 Peruzza （1992）提出的残差最小化方法，Bakun 等（1997）以烈度数据点为基础， 利用衰减关系通过网格搜索来试算历史地震震级和震中区域，这种方法曾应用于美国加州和 美国西北太平洋地区（Bakun 等，1997；Bakun ，1999；Bakun 等，2003 ；Bakun ，2006a ； Diane ，2009 ）、日本（Bakun ，2005 ）、德国（Hinzen 等，2001 ）、法国（Bakun 等，2006b ） 及我国华北地区（张扬等，2009 ）。但该方法是以圆烈度衰减关系为基础，没有考虑烈度衰减 的方向性问题，而且宏观震中与仪器震中有一定的混淆。Gasperini 等（1999）用一套基于烈 度数据点的 Boxer 方法来估算大震级历史地震的震源参数。但该方法首先将最高烈度点的几