广义极值模型在早期强余震预报中的应用.pdfVIP

第23卷 第6期 广东石油化工学院学报 V01．23 No．6 2013年 12月 JoumalofGuangdongUniversityofPetrochemicalTechnology Dec．20l3 广义极值模型在早期强余震预报中的应用 贾朝勇，潘玉荣 (蚌埠学院数理系，安徽蚌埠 233030) 摘要 ：一次主震发生后，往往会产生很多余震，其中强余震的危害性很大。以汶川余震震级资料为基础 ，对余震数据建立广 义极值模型，应用最大似然估计计算出相应参数，并进行了分布拟合检验。利用建立的广义极值模型预测汶川震后一一个月 内最大余震震级和发生的概率，发现与实际情况基本吻合，这说明了极值理论是早期强余震预报中一种行之有效的方法。 关键词：极值理论；广义极值模型；汶川地震；强余震 中图分类号：P315．8 文献标识码：A 文章编号：2095—2562(20l3)06—0065—03 0 引言 极值理论是专门研究很少发生，一旦发生却会产生巨大影响的随机变量极端变异性的建模及统计分 析方法。 目前，极值理论的应用已经深入到许多领域，如水文、气象、金融、地震等…。1945年，诺德魁斯特 (Nordquist)首次把极值理论(冈贝尔分布)用于地震震级统计中 J。之后，许多国内外学者应用和发展了 极值理论，考虑到地震震级有上界，陈培善、林邦慧对冈贝尔分布作出了修正(极值 I型修正式)，使其更加 符合震级分布的实际情况b]。张卫东等利用四种极值分布对郯一营、京一渤两个地震带进行了研究，发现 极值Ⅲ型和极值 I型修正式与实际观测资料符合的较好，他们还对极值分布模型中的参数估计 、检验方 法、应用范围和预测效果进行了探讨，拓展了极值分布在地震危险性分析中的应用范围 J。为避免人为性 和随意l生，钱小仕等采取了广义极值分布对台湾地区最大震级分布进行建模，得到最大震级超过7级的地 震理论发震次数与实际发震次数完全一致，取得了很好的效果 ，而且避免了极值模型选择不当带来的影 响 。 余震也是地震，在主震后的一段时期内余震震级分布规律具有相对稳定性，这主要是因为该地区的内 部环境基本不变(主要受到主震的影响)，因而本文将借鉴地震预报的方法对余震情况进行研究和预报。 事实上，余震活动是极复杂的自然现象，主震后各个地区的余震震级分布存在差异，所以，在模型的选择上 检验实际观测的最大余震震级序列来 自那种理论模型分布是至关重要。现以汶川余震震级资料为基础， 对余震危险}生进行分析，鉴于所有的极值模型都是广义极值模型特例，在此，我们使用广义极值分布对强 余震震级进行研究和预测，分析了震后一个月余震危险性。 1 模型的选择 设 ， ，…， 为独立同分布的随机变量序列，其分布函数为 ，()。称 Mn=max{X，X：…，Xn}为 样本极大值，它的分布称为极值分布。如果存在常数列 { 0}和 {b∈E}，以及一个非退化的分布函数 H，使得 收稿 日期：2013—05—30；修回日期：2013—06—28 基金项 目：安徽省高等学校优秀青年教师人才科研资助项 目(2012SQRL215，2010SQRL115)；蚌埠学院教研项 目(2013jyxn~) 作者简介：贾朝勇(1978一)，男，安徽郎溪人，硕士，讲师，主要从事概率统计的教学和极值理论应用方面的研究。 广东石油化工学院学报 2013每 · lim( ≤戈)：H()，∈R 成立，则 H()必属于下列三种分布族之一： I型分布 日1()=e一 ，∈R Ⅱ型分布 f0， ≤ H2(x)={ 00 【e一( )～ ， Ⅲ型分