论文部分内容阅读
在区域地震危险性分析中,地震预测效能评价、参考模型的构建和区域地震预测策略研究是重要的三个基本问题。而在当前我国地震预测研究和预测预报业务工作体系的“转型升级”过程中,迫切需要构建不同时间尺度的、效能可评价的地震预测模型,开展地震预测效能评价和区域地震预测策略研究。与此相对应,目前国际上正在开展的以严格的统计检验、“竞赛式”预测模型构建等为主的“地震可预测性合作研究”(CSEP)计划,以及“可操作的地震预测”(OEF)等工作对我国相关工作具有重要的借鉴意义。本论文针对上述三个基本问题,借鉴国际经验并结合中国区域地震预测预报实际,开展适合我国区域地震活动特点的地震预测效能评价、参考模型构建和比较研究、区域地震预测策略研究,试图为未来构建新型的区域地震预测预报体系提供科学技术参考。针对上述三方面问题,分别选用了南北地震带(21.0°~43.0°N,97.0°~107.0°E)和“晋冀蒙交界”地区(37.0°~41.0°N,110.0°~116.0°E)为研究区,其中南北地震带为全球“地震可预测性合作研究”(CSEP)计划中国检验区(Chinese Testing Region),而晋冀蒙交界地区先后20余次被划入中国地震局年度重点地震危险区,是华北地区最重要的天然地震预测试验场。数据使用了中国地震台网提供的1970年以来的“全国统一正式目录”。在数据处理中,为获得“背景地震”(background sesimicity)使用了时空ETAS模型和随机除丛法进行了筛选。为选定晋冀蒙交界地区的地震序列目录,采用了区分地震序列与周围背景地震的时-空天然分布界限的“自然边界法”,即通过纬度-时间图、经度-时间图以及震中分布图相结合的方式。为确保地震序列的目录完整性,研究中还使用了“震级-序号”法确定每个地震序列的最小完整性震级(Mc)。针对预测效能评价及其技术问题,论文考虑到当前全球“地震可预测性合作研究”(CSEP)计划和“可操作的地震预测”(OEF)前沿动态,开展了 CSEP计划中通用的Molchan图表法、区域技能评分、N-test、T-test等地震预测检验技术的引入研究和应用。上述方法将分别解决概率预测与危险区划分等“二元”预测结果检验、连续多时段预测的整体评价、地震发生率的预测检验、不同预测模型的比较研究等问题。上述检验方法一方面在优化和构建预测模型中起到重要作用,另一方面可量化评价预测模型和预测结果,这些预测效能评分方法的引入将对解决国内相对单一方法的预测结果评价发挥重要作用。针对预测模型构建中所必须的参考模型问题,论文将全球CSEP计划中国检验区—南北地震带地区为研究区,构建适合该区域地震活动特点的效能可评价的参考模型。参照原理简单、物理意义明确的构建原则,选用简单平滑模型(Triple-S model)和相对强度(RI)模型作为参考模型,使得“参考模型”成为其它地震预测模型继续发展的“竞技标靶”,并在考虑随机除丛等影响因素情况下,与“地震矩加速释放”(AMR)预测模型进行了比较研究。研究中还使用了时空ETAS模型进行除丛处理后获得背景地震目录,用于AMR模型的预测效能检验。相关研究结果表明,构建的两种参考模型(Triple-S模型和RI模型)可适用南北地震带地区的中期尺度预测。“加速矩释放”(AMR)模型的实际预测效果低于参考模型,表明该模型当前的模型设计尚无法满足区域地震预测的需要。针对区域地震预测策略研究问题,论文利用时间序列ETAS模型和“瘦化算法”,对晋冀蒙交界地区开展了连续滑动预测试验及系统的预测效能评估。研究结果表明,ETAS模型和“瘦化算法”对该地区的余震发生率具有一定预测能力,预测“过少”的比例较低,但预测“过多”的比例较高;在主震发生后的早期阶段预测效果较好,但预测效能随着序列的持续时间显著下降;适度增加预测时间窗长可改善预测“过少”的情况,但对预测“过多”的情况改善不明显。此外,高阶余震激发能力偏弱,可能是晋冀蒙交界地区ETAS模型和“瘦化算法”预测效能偏低的原因之一。作为余震短期发生率预测的可能策略,建议在该地区的ETAS模型和“瘦化算法”预测应用中,采用3天的预测时间窗、仅对序列早期阶段进行应用,且重点关注“不低于”相应预测地震数目的底线思维,可能更具有现实意义。