影响b值计算误差的Monte Carlo实验研究 Monte Carlo Experiments on the Influencing Factors of b Value Calculation Errors.pdfVIP

第33卷第4期 地 震 V01．33，No．4 2013年lO月 EARTHQUAKE ()ct．。2013 影响易值计算误差的MonteCarlo实验研究“ 刘善琪，李永兵，田会全，朱伯靖，石耀霖 (中国科学院计算地球动力学重点实验室。中国科学院大学，北京 100049) 摘要：本文假设地震的发生与时问的关系遵循稳态泊松模型，假设地震震级M服从参数为1佃 的指数分布，其中口一6／lo酽，在此基础上生成地震目录。用最小二乘法估算6值，分别从样本 数量、震级间隔、震级误差三方面分析6值误差的大小。统计分析结果表明，6值的标准偏差 随着样本数量的增多而减小；随着震级问隔的增大而减小；震级误差使6值的标准误差变大。 本文可以为利用真实地震目录求6值时选择适当的方法提供参考。 关键词：6值；Monte Carlo；最小二乘法；地震危险性分析 中图分类号：P315．7 文献标识码：A 文章编号：1000一3274(2013)04一0135—10 引言 古登堡和里克特提出全球的地震活动遵循关系logN一以一6M[1。，式中M为地震震级， N是相应震级的地震次数，盘、6为常数。由于当时可用资料有限，认为此关系式不适用M 8的地震。通常认为，n代表一个地区的地震活动性。2‘3j；地震活动水平可以认为是由无 震地动位移速率决定的[4]，“可以看成是这一速率的反应r3’引。6反映该地区大小地震的比 为6值变化取决于岩石所处的应力状态；我国学者的实验则表明，6值变化既与介质的均 匀度有关，又与岩石所处的应力状态和破裂方式有关一]。 研究发现，一般情况下，6值相对稳定，在1附近波动；对于一些地震序列，6值在大 震之前发生异常变化一1n1 3：，即大震发生前6值下降，大震发生后6值又上升。根据6值变 化可以进行地震危险性分析，因此6值的物理机制和变化规律受到重视。国内从20世纪 70年代就有人研究利用6值的变化预测地震，李全林等[141和马鸿庆一”3最早发表了系统研 究成果；近lo年来，国外一些学者系统研究了6值空间变化与地震发生位置的关系¨6~2…。 因为6值的重要性，所以更应注重其精度问题。由于6值是由地震目录用统计方法确 定的，所以与地震资料的完整性、可靠性、统计量的大小、样本的起始震级、取样间隔以 及时空范围等都有关系[21I，因此研究6值统计中的种种影响因素，对合理计算6值是至关 0407 *收稿日期：2013O卜11；修改回日期：2013 基金项目：地壳深部探测专项(sinoProbe07)和国家科技支撑计划项目(2011BAB03809) 作者简介：刘善琪(1987)，女，山东济宁人，在读博士，主要从事矿物、岩石物性模拟等方面的研究。 万方数据 地 震 33卷 重要的。A“221基于最大似然法研究了6值的标准误差和置信度；张建中等[2朝从地震震级 etal[2 M服从指数分布出发，讨论各种估计的标准误差及其求取方法的优劣；Shi 4]得出了 6值标准误差的经验公式；黄玮琼等心门从各种因素取值的合理性及6值在时空域中的稳定 性两方面考察了6值统计中的影响因素；秦长源[2阳研究了地震震级误差对6值的影响；任 雪梅等[263讨论了地震目录的完整性对6值计算的影响。 求6值必须取一定长度(包含了足够数量地震)的时间窗口，然而人们又认为6值是随 时间变化的，那么实际计算的6值的误差究竟是多大，不同时段的6值差异究竟是真实的 差异变化，还是随机的误差起伏，这就成了一个重要的问题。由于真实的地震序列难以排 除6值自身变化的影响，因此本文采用人工生成6值为常数的地震目录，在此基础上，采 用最小二乘法来估算6值，分别从样本数量、震级问隔、震级误差三方面分析6值误差的 大小，从而为利用真实地震目录求6值时选择适当的方法提供参考。 1