精河6.6级地震前精河水管仪、伸缩仪异常特征分析1

张智慧 李 娜 邢喜民



精河6.6级地震前精河水管仪、伸缩仪异常特征分析1

张智慧1）李 娜2）邢喜民2）

1）河南省工业科技学校，河南新乡 453000 2）新疆维吾尔自治区地震局，乌鲁木齐 830011

2017年8月9日精河发生6.6级地震。本文首先利用相关、回归分析研究了气温对水管仪、伸缩仪的影响，并进行了干扰的定量剔除；其次，通过回归残差分析研究地震前水管仪、伸缩仪的异常特征，认为精河6.6级地震前存在3项异常：①2016年5月22日起，水管仪NS分量S倾幅度0.53"；②2017年3月24日起，EW分量E倾幅度0.28"；③2016年9月10日起，伸缩仪NS分量出现压缩现象。

相关系数 气温 回归分析 水管仪 伸缩仪

引言

由构造运动引起的地应力和地倾斜变化，在大地震发生前可能发生异常变化。由于台站建立在地壳表层，位于岩石圈与大气圈、水圈、生物圈的交界面上，其观测不可避免地受到诸多因素影响（中国地震局监测预报司，2008）。许多台站能观测到地应力（变）和倾斜的年度变化，且幅度因台站条件不同而各异，造成这种随季节变化的因素包括气温、气压、降水、地下水位的变化等，温度引起热弹性变形也被用来解释地应力和地形变的年变化（曹建玲等，2005）。Berge（1975）曾对正弦变化的沿水平地表传播的热波动引起的地应变和地倾斜给出了理论解。曹建玲等（2005）利用有限元方法考虑在平行山脊的地形条件下，气温年变化造成的应力和倾斜的幅度及其随深度的分布特征。就目前的认识水平，外界干扰与地震前兆的变化形态有时难以区分（王梅等，2004）。于桂春等（2011）仅定性地分析了气温与阿合奇地震台地倾斜变化关系。邢喜民等（2016）研究了温度、水位、气压对乌什体应变的影响，利用分段回归研究气象因素对精河、库尔勒水平摆的影响（邢喜民等，2015），然而并没有研究定量剔除气象因素之后的形变观测与周边强震之间的关系。

由于地壳构造的复杂性和震源区的不可直观性，地震的孕育和发生、其成因和机制等问题，至今尚无完满的解答（杨龙，2015）。众多大地震的震例研究表明，前兆异常会或多或少地出现在地震孕育过程的不同阶段（吴中海等，2013）。根据中国地震台网测定，2017年8月9日07时27分，新疆博尔塔拉蒙古自治州精河县发生6.6级地震（44.27°N，82.89°E）。本文首先尝试利用相关、回归分析研究气温对精河伸缩仪、水管仪的影响特征，并将其所造成的干扰定量剔除；其次，通过对回归残差进行分析，研究精河6.6级地震前精河伸缩仪、水管仪的异常特征，以期推动地震数值预测方面的研究。

1 台站、仪器概况

精河地震台（44.62°N，82.97°E）位于新疆维吾尔自治区博尔塔拉蒙古自治州精河县境内，海拔387.3m，处于博罗科努、准格尔南缘、科古琴三大断裂带交汇区，该断裂以北为准噶尔中新生代凹陷和围绕准噶尔古老地块形成的准噶尔-北天山古生代褶皱系；断裂以南为塔里木地台北缘的多旋回地槽在晚古生代末期形成的天山褶皱系。历史上该区曾多次发生6.0级以上的中强地震。精河地震台与精河6.6级地震分布如图1所示。

图1 精河地震台与精河6.6级地震分布

精河台观测山洞东西向为主洞，进深约90m，南北洞长约12m，宽约2.2m，洞内年温差0.5℃，日温差小于0.1℃，精河FSQ水管仪、SSY-II型伸缩仪于1990年投入使用，水管、伸缩北南向长9.04m，东西向长29.98m（朱令人，2002）。2011年台站进行了数字化改造，仪器观测精度高，稳定性好，抗干扰能力强。精河水管、伸缩台站观测环境、观测条件均达到观测规范要求。

2 气温对水管仪、伸缩仪影响的定量剔除

为研究精河水管仪EW向、伸缩仪NS向分量与气温之间的关系，选取2013年以来水管仪、伸缩仪和气温的观测数据为研究对象。由日均值曲线图（图2）可以看出，精河水管仪EW向具有趋势E倾、伸缩仪NS向具有趋势压缩的变化，而精河的气温却没有与之一致的趋势变化。根据干扰异常识别与排除的4个“相关性”原则（车用太等，2011），认为精河水管仪EW向、伸缩仪NS向的趋势性变化与气温无关。选用最佳逼近多项式对精河水管仪EW向、伸缩仪NS向进行拟合，去掉2分量观测资料趋势，结果见图3。

（a）水管仪NS向；（b）水管仪EW向；（c）伸缩仪NS向；（d）伸缩仪EW向；（e）气温

（a）水管仪EW向；（b）伸缩仪NS向

通过研究精河水管仪、伸缩仪各分量与气温的相关系数和滞后天数之间的关系，得到了相关系数与滞后天数的计算结果（表1）及关系曲线（图4）。同时，为判定观测数据与气温的相关特征，绘制了其与气温的散点图，如图5所示。

表1 精河水管仪、伸缩仪各分量和气温的相关系数与滞后天数

（a）水管仪NS向；（b）水管仪EW向；（c）伸缩仪NS向；（d）伸缩仪EW向

（a）水管仪NS分量；（b）水管仪EW分量

（c）水管仪EW分量去倾后；（d）伸缩仪NS分量；（e）伸缩仪NS分量去倾后；（f）伸缩仪EW分量

由表1及图5可知，为定量剔除气温对精河水管仪、伸缩仪各分量的影响，以气温作为自变量，选用适合的回归模型进行回归分析。为了方便记录，记定点形变观测数据为，气温为（℃），残差为，并通过可决系数评价各回归方程的拟合效果，其结果见表2。

表2 精河水管仪、伸缩仪与气温回归分析结果

通过回归分析得到的回归残差序列见图6。由上述分析可知，气温是影响精河水管仪、伸缩仪年变的主要影响因素，气温与水管仪、伸缩仪各分量具有准线性关系，且水管仪NS、EW向，伸缩仪NS、EW向分别滞后气温天数约为33、92、82、73天。

（a）水管仪NS向；（b）水管仪EW向；（c）伸缩仪NS向；（d）伸缩仪EW向

3 精河6.6级地震前的异常特征分析

由图6可以看出，在剔除气温对精河水管仪、伸缩仪的影响之后，2013年起水管仪NS分量持续N倾，并于2016年5月22日出现转向S倾的异常，S倾幅度达0.53"，精河6.6级地震发生后异常恢复；2014年沙湾5.0级地震前，水管仪EW分量出现E倾变化，并于2015年9月由E倾转为W倾，之后发生了新源5.0级地震；在精河6.6级地震前，2017年3月24日水管仪再次出现E倾异常变化，至精河地震发生时E倾幅度达0.28"，地震发生后持续E倾，2017年9月29日E倾结束，异常恢复。此外，自2013年起精河伸缩仪NS分量变化较平稳，而自2016年9月10日起出现拉张的异常变化，2017年5月20日拉张速率减慢，精河6.6级地震发生后出现压缩变化，异常恢复；2017年以来伸缩仪EW向虽有拉张迹象，但异常并不明显。由精河6.6级地震构造应力场图（图7）可以看出，精河地震台附近区域的构造应力场水平最大主应力优势方向为近NS向，区域压应力轴为近NS向，这说明精河6.6级地震前伸缩仪NS分量的趋势压缩异常是可靠的。

图7 精河6.6级地震构造应力场图（引自中国地震局地壳应力研究所）

4 结论与讨论

形变台站（测点）作为多输入、单输出（观测值）的系统，在观测值序列中不仅包含着来自地球内部的地形变信息，还包含着各种影响因素的信息。观测值序列既然是多种信息的综合，必然也是可分离的（中国地震局监测预报司，2008）。

本文利用相关、回归分析研究了气温对精河水管仪、伸缩仪的影响特征，并尝试定量剔除气温对其影响，结果表明精河水管仪、伸缩仪观测数据的趋势性变化与气温无关；气温与精河水管仪、伸缩仪各分量具有准线性关系；气温是影响精河水管仪、伸缩仪各分量年变形态的主要影响因素；精河水管仪NS、EW向，伸缩仪NS、EW向滞后气温天数分别为33、92、82、73天；另外，对回归残差序列进行了分析，在精河6.6级地震前出现水管仪NS分量S倾0.53"，EW分量E倾0.28"以及2016年9月10日起伸缩仪的压缩变化等异常变化。通过对比图2、图6可知，剔除气象因素对形变观测数据的影响后，震前异常特征与幅度将更接近由地震引起的异常变化。今后，可通过对大量震例进行总结，尝试建立异常幅度、异常持续时间等异常特征与地震关系，以期推动地震数值预测的发展。

曹建玲，石耀霖，2005．地表温度年变化对地应力和地倾斜的影响．中国科学院研究生院学报，22（3）：303—308．

车用太，鱼金子，刘成龙等，2011．判别地下水异常的干扰性与前兆性的原则及其应用实例．地震学报，33（6）：800—808．

王梅，李峰，孔向阳等，2004．数字化形变观测干扰识别．大地测量与地球动力学，24（1）：94—98．

吴中海，赵根模，2013．地震预报现状及相关问题综述．地质通报，32（10）：1493—1512．

邢喜民，张涛，2015．分段回归在剔除精河、库尔勒水平摆气象因素影响的探索．地震工程学报，37（2）：623—628．

邢喜民，孙吉泽，张治广，2016．乌什台钻孔体应变年频段及日频段信息的影响因素探索，大地测量与地球动力学，36（5）：456—459．

杨龙，2015．钢筋沥青隔震装置的位移控制研究及工程应用．长沙：湖南大学．

于桂春，王在华，刘辉等，2011．阿合奇地震台地倾斜观测与气温变化关系．地震地磁观测与研究，32（6）：143—149．

中国地震局监测预报司，2008．地形变测量．北京：地震出版社．

朱令人，2002．新疆通志-第十一卷-地震志．乌鲁木齐：新疆人民出版社．

Berger J., 1975. A note on thermoelastic strains and tilts. Journal of Geophysical Research, 80(2): 274—277

Analysis of Anomalous Characteristics of Jinghe Water Tube Tilt and Extensometer Temperature before the JingheS6.6 Earthquake

Zhang Zhihui1), Li Na2)and Xing Ximin2)

1) Henan Industrial Technology Institute, Xinxiang 453000, Henan, China 2) Earthquake Agency of Xinjiang Uygur Autonomous Region, Urumqi 830011, China

TheS6.6 earthquake occurred on August 9, 2017 in Jinghe. In order to investigate the anomaly characteristics from water tube tiltmeter and extensometer in Jinghe station, we analyzed the effects of air temperature on water tube tiltmeter and extensometer in Jinghe by using the correlation and regression analysis then rejected the interference of air temperature quantitatively. The analysis of regression residuals shows that there existed three abnormalities before the earthquake: ① S tilt amplitude 0.53" in NS component from water tube tiltmeter since May 22, 2016; ② EW component tilt amplitude 0.28" from March 24, 2017; ③From September 10, 2016, the NS component of the telescopic instrument appears to be compressed.

Correlation coefficient; Air temperature; Regression analysis; Water tube tiltmeter; Extensometer

张智慧，李娜，邢喜民，2018．精河6.6级地震前精河水管仪、伸缩仪异常特征分析．震灾防御技术，13（3）：671—678．

10.11899/zzfy20180318

中国地震局地震科技星火计划（XH19049），中国地震局2018年度地震监测、预测、科研三结合课题项目（CEA- JC/3JH-173103）

2017-12-20

张智慧，男，生于1982年。讲师。主要从事数学的实际应用与教学研究。E-mail：632821511@qq.com

邢喜民，男，生于1980年。高级工程师。主要从事大地测量、定点形变数据分析方面研究。 E-mail：442635690@ qq.com