论文部分内容阅读
[摘 要]典型的大地震老震区出现某些特征性异常变化与外围发生破坏性地震有着密切的关系。本文重点对典型老震区地震窗活动反常预报效能进行探究。
[关键词]地震；活动反常；预报；老震区
中图分类号：P315.75 文献标识码：A 文章编号：1009-914X（2015）06-0190-01
0 引言
本文以海城地震窗的地震活动时间序列为例，提取该序列中的不规则随机部分，以一定的突变量为异常判据，定量地对华北地区6级强震和地震窗附近5级以上地震进行了预测效能的探索研究。海城地震窗位于中国东部。中国东部的政治、经济、文化发达，人口众多，特别是京、津、唐和沿海地区发生6级以上地震会造成巨大损失。因此，分析、判定和预测这些地区的地震形势非常重要。1975年海城MS7.3地震是我国实现破坏性地震成功预报的典型震例，这次地震本身及其相关的资料具有很高的科学研究价值。该地震在地震构造划分上位于一级活动地块的边界带上，是发生强震的主要构造部位。地震已过去36年，震区积累了丰富的余震记录，为地震形势和地震预测研究提供了宝贵资料。为此，本文研究了海城地震窗地震频次时间序列，分解、提取了异常并进行预报效能评价。对地震形势分析、判定和预测有科学参考价值。
1 海城地震窗范围确定及地震资料选取
地震窗是区域应力场与震源应力场变动效应的一种地震活动表现，我们可以把握地震窗的应力变化来实现对相关地震区的地震预测.地震窗的应力变化可以用其地震活动起伏来反映，但在研究时必须解决的关键问题是合理地确定好地震窗的边界范围，因为不同的范围统计得到的地震起伏结果不同。本文综合断裂构造背景和地震空间分布特征来确定海城地震窗范围。分析1975年海城MS7.3地震的地震断裂构造，北东向的营口—海城断裂、金州断裂和太子河断裂是主要控震构造，皮口断裂和碧流河断裂是规模较小的次要控震构造；北西向的大洋河断裂是发震构造。这些断裂都与海城MS7.3地震孕震发生有关。控震构造北东向展布，大致300km长；发震构造北西向展布，大致100km长。按照上述地震构造条件（控震构造和发震构造）的展布，结合海城地震以来的小震空间分布，本文使用这两个约束条件合理地勾画了海城地震窗的边界范围.采用该边界范围选取了1975-2011年2月的地震资料进行海城地震窗地震活动起伏及映震效能分析。采用上述方法确定的范围选取地震资料，还应考虑震级较小地震和震中定位精度对地震起伏的影响。选取资料时应做到在台网最低监测能力一致的基础上进行，即确定好最低震级下限，震级下限要由观测多年的台网监测能力来决定。震中定位精度依赖于诸多条件：台站分布的状况，震相识别及其到时读取的准确程度，以及地壳结构的认识程度等。
2 海城地震窗地震活动异常提取及对华北强震的预测效能
（1）异常提取方法设XT是一个时间观测数据序列，据加法分解原理则有式（1）包括长期趋势周期（循环和季节）和不规则随机两部分。Mt是关于时间的多项式函数，它是XT数据序列中包含长期趋势、循环周期、季节周期成分组合部分项，反映观测数据本身具有的内在规律和特性，是可以预测的。It为XT数据序列中随时间的不规则随机扰动变化项，是受误差和主观或客观方面不可预测之突然、偶然异常因素造成的。为了达到提取异常目的，应根据研究目的要求分离需要项。为此采用对XT数据序列建立自回归单整移动平均模型求取Mt。在得到Mt（T，C，S）的情况下，It即可自Xt中分离出来，用来进行异常研究分析。由于海城地震窗月频次时间序列的自相关函数不存在季节频度的波动循环.
（2）建立海城地震窗地震月频次时间序列的自回归单整移动平均模型由月频次图可以看出，海城地震窗月频次时间序列是非平稳的。首先数据变化幅度较大，故采取对数处理；其次数据序列存在上下趋势和某些周期性变化，经检验中增广迪基-福勒，由此采取一阶差分处理。经过上述处理，较大幅度和趋势的变化剔除，数据序列成为正态平稳化，满足建模要求。计算处理后数据的自相关函数和偏相关函数值，模型合理。使用该模型计算可分离得到海城地震窗包含参差和受突然或偶然等异常因素造成的。
（3）预测效能分析做海城地震窗震级ML≥1.8的1975-2011年2月地震月频次曲线，在获得解释模型的基础上，分离出不规则随机扰动项，本文称之为提取异常。对华北地区MS6.0和地震窗附近地区MS5.0以上地震做预测反应分析。经3种异常判据对比分析，以扰动值高于2倍方差线为异常判据.在研究时间内，应除去大震初期1975年高值异常。如把海城地震刚发生这次高值对应成1976年发生强震（和林格尔MS6.4、唐山MS7.8）的异常说服力不强，为此，本文做预测反应分析的效能评价时不包括大震初期高值异常。当发生强震群时，震前的异常也多次出现，如大同强震群；当发生单发或双震地震时，震前的异常往往出现1次或2次；1999年高值异常对应在附近地区（即地震窗本身）发生的岫岩MS5.4地震；其它的高值异常对应华北地区发生MS6.0以上地震；1998年张北MS6.2地震前虽未出现2σ异常，但有1.5σ的异常量，看来强震前的异常频次和幅度可能与地震的距离、构造走向有关，这需作更多地震应力窗的统计研究。在422个月的研究时间内共发生了11次地震。本文规定高值异常出现后，预报1年内对应地震，1年内不发生地震情况下，再延迟半年的对应期，在1年半的预报期发生地震为报准；在1年半的预报期未发生地震为虚报；发生地震前1年半内无高值异常为漏报。按照许绍燮（1989）的评分方法，进行了预测效能评价计算。结果显示：报准10次，虚报4次，漏报1次，R=0.497，具有97.5置信水平的R0值为0.323。因此，该方法具有较高信度和预报价值。
3 讨论与结论
本研究结果表明，利用具有较高监测能力的台网监测统计海城地震窗地震活动时间序列，有效地提取地震异常并对历史震例预报效能评价，是切实有效的震情预测方法。实践中，在地震构造区，可以联合几次大地震老震区窗口形成震情网络窗口，提取同步地震异常，能更加强预测的有效性。
参考文献
[1] 徐敬海，刘伟庆，邓民宪.GIS技术在城市防震减灾的应用[D].南京：东南大学，2002.
[2] 姚宝华，谢礼立，陶夏新.基于Web GIS的防震减灾系统研究及其实现方法[J].地震工程与工程震动，2010，23（l）：l-8.
[关键词]地震；活动反常；预报；老震区
中图分类号：P315.75 文献标识码：A 文章编号：1009-914X（2015）06-0190-01
0 引言
本文以海城地震窗的地震活动时间序列为例，提取该序列中的不规则随机部分，以一定的突变量为异常判据，定量地对华北地区6级强震和地震窗附近5级以上地震进行了预测效能的探索研究。海城地震窗位于中国东部。中国东部的政治、经济、文化发达，人口众多，特别是京、津、唐和沿海地区发生6级以上地震会造成巨大损失。因此，分析、判定和预测这些地区的地震形势非常重要。1975年海城MS7.3地震是我国实现破坏性地震成功预报的典型震例，这次地震本身及其相关的资料具有很高的科学研究价值。该地震在地震构造划分上位于一级活动地块的边界带上，是发生强震的主要构造部位。地震已过去36年，震区积累了丰富的余震记录，为地震形势和地震预测研究提供了宝贵资料。为此，本文研究了海城地震窗地震频次时间序列，分解、提取了异常并进行预报效能评价。对地震形势分析、判定和预测有科学参考价值。
1 海城地震窗范围确定及地震资料选取
地震窗是区域应力场与震源应力场变动效应的一种地震活动表现，我们可以把握地震窗的应力变化来实现对相关地震区的地震预测.地震窗的应力变化可以用其地震活动起伏来反映，但在研究时必须解决的关键问题是合理地确定好地震窗的边界范围，因为不同的范围统计得到的地震起伏结果不同。本文综合断裂构造背景和地震空间分布特征来确定海城地震窗范围。分析1975年海城MS7.3地震的地震断裂构造，北东向的营口—海城断裂、金州断裂和太子河断裂是主要控震构造，皮口断裂和碧流河断裂是规模较小的次要控震构造；北西向的大洋河断裂是发震构造。这些断裂都与海城MS7.3地震孕震发生有关。控震构造北东向展布，大致300km长；发震构造北西向展布，大致100km长。按照上述地震构造条件（控震构造和发震构造）的展布，结合海城地震以来的小震空间分布，本文使用这两个约束条件合理地勾画了海城地震窗的边界范围.采用该边界范围选取了1975-2011年2月的地震资料进行海城地震窗地震活动起伏及映震效能分析。采用上述方法确定的范围选取地震资料，还应考虑震级较小地震和震中定位精度对地震起伏的影响。选取资料时应做到在台网最低监测能力一致的基础上进行，即确定好最低震级下限，震级下限要由观测多年的台网监测能力来决定。震中定位精度依赖于诸多条件：台站分布的状况，震相识别及其到时读取的准确程度，以及地壳结构的认识程度等。
2 海城地震窗地震活动异常提取及对华北强震的预测效能
（1）异常提取方法设XT是一个时间观测数据序列，据加法分解原理则有式（1）包括长期趋势周期（循环和季节）和不规则随机两部分。Mt是关于时间的多项式函数，它是XT数据序列中包含长期趋势、循环周期、季节周期成分组合部分项，反映观测数据本身具有的内在规律和特性，是可以预测的。It为XT数据序列中随时间的不规则随机扰动变化项，是受误差和主观或客观方面不可预测之突然、偶然异常因素造成的。为了达到提取异常目的，应根据研究目的要求分离需要项。为此采用对XT数据序列建立自回归单整移动平均模型求取Mt。在得到Mt（T，C，S）的情况下，It即可自Xt中分离出来，用来进行异常研究分析。由于海城地震窗月频次时间序列的自相关函数不存在季节频度的波动循环.
（2）建立海城地震窗地震月频次时间序列的自回归单整移动平均模型由月频次图可以看出，海城地震窗月频次时间序列是非平稳的。首先数据变化幅度较大，故采取对数处理；其次数据序列存在上下趋势和某些周期性变化，经检验中增广迪基-福勒，由此采取一阶差分处理。经过上述处理，较大幅度和趋势的变化剔除，数据序列成为正态平稳化，满足建模要求。计算处理后数据的自相关函数和偏相关函数值，模型合理。使用该模型计算可分离得到海城地震窗包含参差和受突然或偶然等异常因素造成的。
（3）预测效能分析做海城地震窗震级ML≥1.8的1975-2011年2月地震月频次曲线，在获得解释模型的基础上，分离出不规则随机扰动项，本文称之为提取异常。对华北地区MS6.0和地震窗附近地区MS5.0以上地震做预测反应分析。经3种异常判据对比分析，以扰动值高于2倍方差线为异常判据.在研究时间内，应除去大震初期1975年高值异常。如把海城地震刚发生这次高值对应成1976年发生强震（和林格尔MS6.4、唐山MS7.8）的异常说服力不强，为此，本文做预测反应分析的效能评价时不包括大震初期高值异常。当发生强震群时，震前的异常也多次出现，如大同强震群；当发生单发或双震地震时，震前的异常往往出现1次或2次；1999年高值异常对应在附近地区（即地震窗本身）发生的岫岩MS5.4地震；其它的高值异常对应华北地区发生MS6.0以上地震；1998年张北MS6.2地震前虽未出现2σ异常，但有1.5σ的异常量，看来强震前的异常频次和幅度可能与地震的距离、构造走向有关，这需作更多地震应力窗的统计研究。在422个月的研究时间内共发生了11次地震。本文规定高值异常出现后，预报1年内对应地震，1年内不发生地震情况下，再延迟半年的对应期，在1年半的预报期发生地震为报准；在1年半的预报期未发生地震为虚报；发生地震前1年半内无高值异常为漏报。按照许绍燮（1989）的评分方法，进行了预测效能评价计算。结果显示：报准10次，虚报4次，漏报1次，R=0.497，具有97.5置信水平的R0值为0.323。因此，该方法具有较高信度和预报价值。
3 讨论与结论
本研究结果表明，利用具有较高监测能力的台网监测统计海城地震窗地震活动时间序列，有效地提取地震异常并对历史震例预报效能评价，是切实有效的震情预测方法。实践中，在地震构造区，可以联合几次大地震老震区窗口形成震情网络窗口，提取同步地震异常，能更加强预测的有效性。
参考文献
[1] 徐敬海，刘伟庆，邓民宪.GIS技术在城市防震减灾的应用[D].南京：东南大学，2002.
[2] 姚宝华，谢礼立，陶夏新.基于Web GIS的防震减灾系统研究及其实现方法[J].地震工程与工程震动，2010，23（l）：l-8.