1. 引言
许多年来,国内外科学家公认的是地震前兆具有较大的复杂性,而地震预报必须要依据某种前兆来进行,所以首先我们就要加强地震前兆理论的研究和实践。地震前兆异常的研究是地震预测预报的基础,对地震预报工作具有重要的指导作用。众多学者通过运用前兆观测资料对诸多震例开展了震前异常特征分析,取得了许多有益的研究成果。邱永平等 [1] 对宁波地震台与台湾三次中强震处于近南北轴线的特殊位置进行了同震响应探讨,发现在震前垂直摆EW分量出现阶跃异常现象;蒋薇等 [2] 对通海地震台垂直摆倾斜仪观测资料在宁蒗5.7级地震前出现的异常进行了分析,结果表明通海台垂直摆倾斜仪在震前出现了明显的中期、短临异常,有明显的趋势转折、破年变、矢量反向、速率超二倍均方差等现象;柯昌安等 [3] 对宁强5.7级地震前陕西定点形变的异常现象进行了分析,阐述了临震前陕西形变观测资料出现了短临异常现象,异常形态趋势转折;李杰等 [4] 对濮阳ML4.6地震前泰安台形变进行异常特征分析,泰安台倾斜仪出现了显著的短、临前兆异常,地震发生在异常达到最大值前后,且震中异常背向初期地壳的倾斜方向。
地形变又称地震形变。在地震发生时或地震前后,由于构造运动或地震构造力和地震波振动作用,地表产生变形的现象。通过地形变观测可监测到震前出现的前驱波、慢地震或蠕变事件等现象 [5] [6]。垂直摆倾斜仪就是地形变测量手段和仪器的其中一种 [7]。
会昌地震台距离台湾地区约530 km,台站垂直摆倾斜仪对台湾中强震反应特别明显。会昌地震台垂直摆倾斜仪安装在山洞内,它能够有效地排除来自地表的各种干扰,取得可信度较高的观测数据。本文对垂直摆2010年至2021年的观测资料与台湾地区Ms ≥ 6.5级地震之间的对应关系及异常图像进行分析,探讨其观测曲线异常特征与地震的关系 [8]。
2. 观测环境
会昌地震台设在江西省赣州市会昌县文武坝镇岚山公园内,距主城区5 km,有简易公路到达台站。台站风景秀丽,无明显干扰,是地震观测的良好场地,具有优质的地质观测条件,台基岩性为加里东期花岗岩,台基坚硬,背景噪声极低,周围无厂矿与输、变电设施,满足综合台站观测规范的要求 [9]。
会昌地震台山洞总长约85米。主巷道安装了三道密封门,洞内温度常年稳定在20℃左右。2007年3月垂直摆倾斜仪布设在台站山洞内,摆房距洞口深约60米,与山洞过道由一扇密封门隔离。垂直摆观测墩直接浇筑在基岩上,随着时间的推移,观测墩已同下面基岩全面融合,周边有隔震槽 [10]。
3. 异常分析
3.1. 资料选取
会昌地震台垂直摆倾斜仪周边无人为干扰,具有较好的观测环境,同时近年来仪器运行稳定,故障断记少,产出了大量连续稳定的观测资料 [11] [12]。本文选用2010年1月至2021年12月的观测资料,以台湾地区Ms ≥ 6.5级地震为研究对象,对此期间垂直摆观测到的震例进行震前异常曲线和矢量图异常特征分析,总结震例异常的特点、形态,并运用一阶差分求速度值和二阶差分求加速度值提取异常 [13] [14] [15]。通过中国地震局地震数据中心,检索了2010年1月1日至2021年12月31日台湾地区Ms ≥ 6.5级地震,共计9个地震。详见表1。

Table 1. Earthquakes with Ms ≥ 6.5 in Taiwan from January 1, 2010 to December 31, 2021
表1. 2010年1月1日至2021年12月31日台湾地区Ms ≥ 6.5级地震
3.2. 震前观测曲线异常分析
2018年2月4日台湾花莲6.6级、2月6日台湾花莲6.8级和2019年4月18日台湾花莲6.7级三个地震发生前后因仪器故障,无法分析(如图1,图2)。对表1中其余6个地震震前一个月内的垂直摆观测资料进行研究分析后发现,除2013年10月31日6.8级地震震前异常不明显外,垂直摆北南分量观测曲线在其余地震震前均出现了不同形态的异常。2010年3月4日台湾高雄6.7级地震前垂直摆北南分量在2月13日至2月16日出现了小幅突跳、波浪形抖动等异常图像(如图3);2010年4月26日台湾东南6.7级地震前在4月9日至11日出现了褶皱、锯齿型连续错动、掉格型突跳等异常图像(如图4);2013年6月2日台湾南投6.7级地震前在5月15日至21日出现了褶皱、固体潮畸变波动等异常图像(如图5);2013年10月31日台湾花莲6.8级地震前一个月内垂直摆北南分量未发现明显异常(如图6);2015年4月20日台湾花莲6.7级地震前在4月9日至17日出现了褶皱、掉格型突跳、波浪形抖动等异常图像(如图7);2016年2月6日台湾高雄6.7级地震前在1月15日至17日出现了掉格型突跳、波浪形抖动等异常图像(如图8)。

Figure 1. Vertical pendulum failure from December 19, 2017 to February 3, 2018
图1. 2017年12月19日至2018年2月3日垂直摆故障

Figure 2. Vertical pendulum failure from March 27 to May 20, 2019
图2. 2019年3月27日至5月20日垂直摆故障

Figure 3. Anomalies before the Ms6.7 Kaohsiung, Taiwan earthquake on March 4, 2010
图3. 2010年3月4日台湾高雄6.7级地震震前异常曲线图

Figure 4. Anomalies before the Ms6.7 southeast Taiwan earthquake on April 26, 2010
图4. 2010年4月26日台湾东南6.7级地震震前异常曲线图

Figure 5. Anomalous curve before the Ms6.7 Nantou Taiwan earthquake on June 2, 2013
图5. 2013年6月2日台湾南投6.7级地震震前异常曲线图

Figure 6. Graph of the earthquake before Nantou Ms6.7 earthquake in Hualien, Taiwan, October 31, 2013
图6. 2013年10月31日台湾花莲6.8级地震震前曲线图

Figure 7. Anomalous curve before the Ms6.7 Hualien Taiwan earthquake on April 20, 2015
图7. 2015年4月20日台湾花莲6.7级地震震前异常曲线图

Figure 8. Anomalous curve before the Ms6.7 Kaohsiung, Taiwan earthquake on February 6, 2016
图8. 2016年2月6日台湾高雄6.7级地震震前异常曲线图
3.3. 震前矢量图异常分析
对表1中2010年3月4日台湾高雄6.7级地震、2010年4月26日台湾东南6.7级地震、2013年6月2日台湾南投6.7级地震、2015年4月20日台湾花莲6.7级地震和2016年2月6日台湾高雄6.7级地震进行震前一个月日均值矢量图研究分析,我们发现在震前日均值矢量图中呈现出临震反向、速率加大变化、拐弯、打结等异常特征(如图9~13),观测资料异常时段仪器运行正常,未发现人为或自然干扰。
3.4. 差分法异常分析
在地震研究中,通常要排除各种干扰,突出前兆信息 [16] [17]。单分量测值的一阶差分和二阶差分反映的是该分量形变速度和加速度的变化,是抑制较长周期变化,突出前兆信息的有效方法 [18] [19] [20] [21]。运用一阶差分和二阶差分分别绘制了表1中2010年3月4日台湾高雄6.7级地震、2010年4月26日台湾东南6.7级地震、2013年6月2日台湾南投6.7级地震、2015年4月20日台湾花莲6.7级地震和2016年2月6日台湾高雄6.7级地震震前速度和加速度值曲线(如图14~18),从各图中可以看出通过一阶差分和二阶差分滤除曲线中较长周期变化后,地震震前均出现了不同程度的高频异常。

Figure 9. Vector image of the earthquake before the magnitude Ms6.7 earthquake in Kaohsiung, Taiwan, March 4, 2010
图9. 2010年3月4日台湾高雄6.7级地震震前矢量图

Figure 10. Vector image of the earthquake before the magnitude Ms6.7 earthquake southeast of Taiwan on April 26, 2010
图10. 2010年4月26日台湾东南6.7级地震震前矢量图

Figure 11. Vector image of the earthquake before the magnitude Ms6.7 earthquake in Nantou, Taiwan, June 2, 2013
图11. 2013年6月2日台湾南投6.7级地震震前矢量图

Figure 12. Vector image of the earthquake before the magnitude Ms6.7 earthquake in Hualien, Taiwan, April 20, 2015
图12. 2015年4月20日台湾花莲6.7级地震震前矢量图

Figure 13. Vector image of the earthquake before the Ms6.7 Kaohsiung, Taiwan, February 6, 2016
图13. 2016年2月6日台湾高雄6.7级地震震前矢量图

Figure 14. Abnormal curves of velocity and acceleration before the Ms6.7 Kaohsiung, Taiwan earthquake on March 4, 2010
图14. 2010年3月4日台湾高雄6.7级地震震前速度、加速度异常曲线图

Figure 15. Abnormal curves of velocity and acceleration before the Ms6.7 southeast Taiwan earthquake on April 26, 2010
图15. 2010年4月26日台湾东南6.7级地震震前速度、加速度异常曲线图

Figure 16. Abnormal curves of velocity and acceleration before the Ms6.7 Nantou Taiwan earthquake on June 2, 2013
图16. 2013年6月2日台湾南投6.7级地震震前速度、加速度异常曲线图

Figure 17. Abnormal curves of velocity and acceleration before the Ms6.7 Hualien Taiwan earthquake on April 20, 2015
图17. 2015年4月20日台湾花莲6.7级地震震前速度、加速度异常曲线图

Figure 18. Abnormal curves of velocity and acceleration before the Ms6.7 Kaohsiung, Taiwan earthquake on Feb 6, 2016
图18. 2016年2月6日台湾高雄6.7级地震震前速度、加速度异常曲线图
4. 结论与讨论
综上所述,会昌地震台垂直摆十多年的观测资料表明会昌台山洞基岩条件好,产出的形变观测资料精度较高,有利于开展地震前兆异常研究。本文通过原始形态、矢量图分析和差分法异常分析三种方法提取的异常特征与台湾地震有良好对应关系,但这也属震后分析总结,说明地震预报研究的难度很大。
1) 会昌地震台观测环境较好,垂直摆观测资料连续稳定,具备良好的映震能力,但对于每个台湾地震震前异常表现也存在一定差异,这也充分说明了地震前兆异常的复杂性,仅为台站地震监测及临震预报提供一种参考。
2) 会昌台垂直摆北南分量在台湾地区Ms ≥ 6.5级地震除2013年10月31日台湾花莲6.8级地震震前异常不明显外,其他地震震前一个月内均出现了不同程度的褶皱、锯齿型连续错动、台阶状、掉格型突跳或波浪型抖动等异常图像。同时在震前日均值矢量图中也呈现出临震反向、速率加大变化、拐弯、打结等异常特征。
3) 通过一阶差分和二阶差分滤除曲线中较长周期变化后,地震震前速度和加速度值曲线均出现了不同程度的高频异常,震前异常对应好,特征明显,易于识别。
基金项目
江西九江扬子块体东部地球动力学野外科学观测研究站开放基金资助项目(OGYB202006)。
NOTES
*通讯作者。