五指山形变台定点形变观测资料内在精度分析

郭昱琴　曾涛　杨世杰　李冬雅

摘要本文对2016年1月—2021年12月连续6年五指山形变台定点形变观测资料内在精度进行统计分析，台站水管倾斜仪和铟瓦棒伸缩仪内在精度质量较好；宽频带倾斜仪两个分量内在精度年均达到优秀指标的概率仅50%，且自2020年1月出现观测数据内在精度质量持续下降，导致观测资料质量不高。提高数据内在高质量的重要条件是确保观测数据的连续率，因此，提升台站人员业务水平和保证必要的备机备件是当务之急。

关键词倾斜仪；应变仪；潮汐因子；内在精度

中图分类号： P315.72+5文献标识码： A文章编号：2096-7780（2023）06-0278-07

doi：10.19987/j.dzkxjz.2022-093

Internal accuracy analysis of fixed-point deformation observation data of Wuzhishan deformation station

Guo Yuqin1），Zeng Tao2），Yang Shijie2），Li Dongya1）

1） Hainan Earthquake Agency， Hainan Haikou 570203， China

2） Qiongzhong Central Station of Hainan Earthquake Agency， Hainan Qiongzhong 572900， China

Abstract Thispaperpresentsastatisticalanalysisoftheinternalaccuracyofthefixed-pointdeformation observation data of Wuzhishan deformation station from January 2016 to December 2021. The internal accuracy of the stations water pipe tiltmeter and indium tile rod extensometer is good. The probability that the internal accuracy of the two components of the broadband tiltmeter reaches the excellent index is only 50%. Since January 2020，the internal accuracy quality of the observation data has been decreasing，resulting in the low quality of the observation data. The important condition for improving the internal accuracy quality of data is to ensure the continuous rate of observation data. Therefore，it is urgent to improve the operational level of station personnel and to ensure the necessary spare parts.

Keywordstiltimeter; strain gauge; tide factor; internal accuracy

引言

有大量震例证明，孕震过程中形变会出现阶段性前兆特征，即形变速率变化、形变方向变化和岩石力学性质变化[1]。其中，前两项特征（形变速率和形变方向变化）可以从地倾斜数据直接获得，而第三项特征（岩石力学性质）则是间接通过固体潮振幅因子及波形包络形态描述获得[2-3]。一般情况下，台站日常监测资料是连续动态观测时间序列数据，序列数据中有多种有用信息，也有噪声影响。为更好地使用资料和保障观测系统正常运行，建立客观、准确评估观测资料序列数据内在精度的方法尤为重要。杨晓东等[4-7]对韩城台、宁波台、格尔木台等的倾斜观测和数字化垂直摆资料进行了对比研究；曹白伦等[8]对云县台水平摆和垂直摆资料进行了对比研究。然而，这些研究成果较多的是从单台套地形变仪器资料得出的，本文针对五指山地震台地倾斜（DSQ 水管倾斜仪和 VP 宽频带倾斜仪）及地应变（ SSY 铟瓦棒伸缩仪）3套形变仪器资料内在精度进行比较，评估台站地形变资料的真实性和分析影响数据质量的因素，为研究该区域地壳应力—应变状态的长期变化及规律提供支持和依据。

1 台站概况

五指山形变台位于五指山市冲山镇阿里山度假村旁太平山山麓半山腰处，地理位置（18.79°N，109.53°E），该区域二叠系花岗岩和三叠系花岗岩出露显著[9]。台站洞体基岩坚硬完整，致密均匀，周边观测环境较为稳定，四周无振动干扰源。“十五”期间，经技术改造，山洞架设了 D SQ 水管倾斜仪、 VP 宽频带倾斜仪和 S SY 铟瓦棒伸缩仪。其中，DSQ 水管倾斜仪和 SSY 铟瓦棒伸缩仪安装在主洞室内，两套仪器相同分量观测部件平行布设在同一槽体，监测该区域的大地倾斜和应变变化，VP 宽频带倾斜仪安装在主洞南侧耳洞室（图1）。3套仪器均能记录到完整的固体潮，观测数据连续性好，同震记录效应明显。

2 资料分析

固体潮观测值相对其理论的偏离通常采用基于最小二乘法的维尼迪柯夫（Venedikov）调和分析方法，求解各个波群的观测振幅与理论振幅之比，即潮汐因子。台站倾斜资料调和分析精度也采用Venedikov调和分析方法，对月尺度倾斜资料进行调和分析，以 M2波中误差 mγ作为评定倾斜观测资料质量的一项指标，规范要求Ι类倾斜台站观测 M2波中误差 m γ＜0.02。应变观测资料调和分析精度同样采用Venedikov调和分析方法，对月尺度的应变观测资料进行调和分析，以 M2波振幅因子的相对中误差σα/α指标作为评定应变资料质量的一項指标，规范要求Ι类应变台站观测资料σα/α＜0.05[10]。

定点形变观测资料内在精度是指 M2波潮汐因子均方差，该指标是用来评定固体潮汐（地倾斜、地应变）观测资料内在质量精度的一项定量化指标。本文采用了Venedikov调和分析法，分析台站观测资料的内在精度（地倾斜仪器精度由 M2波的中误差 mγ表示，地应变仪器精度由 M2波振幅因子相对中误差σα/α表示，文中统称为“内在精度”），也是当前中国地震局地壳形变学科组通常公布的月评内在精度结果。通过统计2016年1月—2021年12月五指山台定点形变观测内在精度分析，我们发现五指山台倾斜应变各仪器各分量内在精度变化规律各不相同（图2）。对同一套观测仪器而言，两个观测分量内在精度存在同步上升下降，如图2a 所示，水管倾斜仪2019年6月、2020年3月、8月的渐变现象；或两个观测分量内在精度变化不一致或单观测分量内在精度下降，如图2b 中洞体应变伸缩仪所描述的。

根据计算仪器观测数据的内在精度，我们再对五指山台3套仪器不同分量计算内在精度均值，然后，定义a为地倾斜或地应变仪器 NS 分量内在精度与 EW 分量内在精度之比，a值可较直观地显示出同一仪器两个分量内在精度的差异[11]，a值越趋近于1，仪器两个分量的内在精度越相近，表明数据内在质量越好，反之亦同。统计数据见表1。

由表1计算数据可知：

（1）水管倾斜仪、铟瓦棒伸缩仪两个分量内在精度较好，a值趋向于1，接近优秀指标；而宽频带倾斜仪两个分量内在精度年均达到优秀满分概率仅50%（图2c），出现自2020年1月地倾斜仪数据内在精度下降趋势，导致数据质量不高。

（2）同台站地倾斜类仪器 NS 分量内在精度均高于 EW 分量，水管倾斜仪和宽频带倾斜 EW 分量内在精度优于 NS 分量（图2）。

（3）宽频带倾斜仪 NS 分量和 EW 分量内在精度平均值差距较大，宽频带倾斜仪的a值为3.5072，铟瓦棒伸缩仪和水管倾斜仪 NS 分量和 EW 分量内在精度接近，a值分别为1.2424、1.0763（表1）。

3 影响内在精度因素探讨

3.1 仪器稳定性

仪器故障是台站数据完整率和连续率的重要影响因素。缺数或干扰太多进而影响数据质量。五指山台周边观测环境较为稳定，四周无振动干扰源，且观测硐室温度年变化为0.1℃。宽频带倾斜仪自安装运行以来，因仪器 NS 分量自身稳定性能差致使后续故障频发，尤其是漂移现象造成调零频繁和超量程。如图3所示，2021年1—12月宽频带倾斜仪数据曲线 NS 分量持续向南倾斜，EW 分量持续向东倾斜。NS 分量数据统计计算内在精度满足0.02评价指标的月份仅5个月，其中2、3、5、6、8、9、10月内在精度均未满足评价指标0.02（图4a），EW 分量要好些，满足0.02评价指标的月份达11个月，仅3月内在精度未满足评价指标0.02（图4b）。

3.2 人为因素

影响数据的人为因素主要有检修仪器过程、人行走进洞干扰、人为操作调零及标定等，表现为数据突跳、数据出错误或阶跃变化。如2020年5月水管倾斜仪主机供电模块故障，剔除产出的错误数据，出现缺数情况（图5），导致无法统计计算内在精度（图6）。另外还存在其他的人为因素，如台站人员运维水平有限，对台站仪器性能了解、故障排查检修能力不足，只能等待厂家维修，导致故障断记时间较长，直接影响数据内在精度。

3.3 同震影响

地震对台站3套仪器观测数据的影响表现上阶跃、下阶跃或阶跃式振荡重叠同震响应现象，在统计2016年1月—2021年12月连续6年记录数据中发现，伸缩仪对应7级以上地震才出现同震效应。而水管倾斜仪和 VP 垂直摆对应不同震中距的地震响应存在较大差异：

（1）震中距较小时，水管仪基本无同震响应，VP 垂直摆同震波记录清晰；

（2）震中距＞6000km 时，两种观测仪器记录应震能力相当[12]。这种响应主要以高频振荡衰减的形式叠加在所记录的固体潮曲线上，震级越大，同震阶变数量越多，震荡时间越长，数据内在精度会越差（图7）。如宽频带倾斜仪受同震干扰影响有2021年 7月和8月的数据内在精度与其他月份及同期的水管倾斜仪和伸缩仪相比均表现为能力明显降低（图4）。

4 结语

通过统计2016年1月—2021年12月五指山台定点形变观测内在精度变化情况，得出以下结论：

（1）水管倾斜仪和铟瓦棒伸缩仪两个分量内在精度较好，基本满足优秀评价指标0.02；宽频带倾斜仪两个分量内在精度年均达到优秀指标的概率仅50%（图2c）；自2020年1月内在精度持续下降，超量程故障频繁，观测资料质量不高，亟待需要对仪器进行全方面检修或更新改造。

（2）地倾斜仪与地应变观测仪观测数据均为 EW 分量内在精度优于 NS 分量，分析其原因可能是安裝仪器正交方位或夹角参数不准确导致不一致所致。同时，可能原因是山洞仪器长期运行老化性能变差，或山洞地质构造环境特殊，这方面尚未做进一步探讨。

（3）仪器稳定性、人为因素和同震干扰等因素都会影响数据质量，导致数据内在精度降低。因此，提高台站仪器运行率，保证记录数据连续完整性，是提高数据内在精度的一个重要技术环节。当然，客观干扰因素无法避免，如每年沿海台风、雷暴雨、地震或近场爆破等。另外，尽快提升台站人员熟悉仪器原理、性能业务水平和仪器故障排查技术能力势在必行，配置必要的备机备件也应当考虑。所以亟待解决的问题是台站人才培养、业务培训和备机备件。

参考文献

［1］刘晓霞，江在森，武艳强.唐山、汶川地震地壳形变特征综合研究[J].地球物理学进展，2018，33（3）：941-950

Liu X X，Jiang Z S，Wu Y Q. Characteristic of crustal deformation related to Tangshan and Wenchuan earthquakes[J]. Progress in Geophysics，2018，33（3）：941-950

［2］张英杰，江永建，李晓鹏，等.河北涉县台伸缩仪稳定性及资料与映震分析[J].地震工程学报，2018，40（增刊1）：112-116

Zhang Y J，Jiang Y J，Li X P，et al. Stability and reflecting earthquake analysis of the extensometer at Shexian Seismic Station[J]. China Earthquake Engineering Journal，2018，40（S1）：112-116

［3］史小平，闫万龙，龙艳.白银形变观测站 DSQ 型水管倾斜仪观测资料分析[J].地震地磁观测与研究，2014，35（3/4）：178-181

Shi X P，Yan W L，Long Y. Preliminary analysis on observational data of DSQ type water tube tilted instrument at Baiyin Seismic Station[J]. Seismological and Geomagnetic Observation and Research，2014，35（3/4）：178-181

［4］杨晓东，张晨蕾，吴明，等.陕西省 VP 型垂直摆倾斜仪观测对比分析[J].地震科学进展，2021，51（8）：371-375

Yang X D，Zhang C L，Wu M，et al. Comparative analysis of VP type vertical pendulum tiltmeter in Shaanxi Province[J]. Progress in Earthquake Sciences，2021，51（8）：371-375

［5］蘇维刚，冯丽丽，马茹莹，等.水平摆与垂直摆倾斜仪观测资料干扰特征对比[J].地震地磁观测与研究，2018，39（4）：149-153

Su W G，Feng L L，Ma R Y，et al. Comparative of disturbance characteristics from horizontal pendulum and vertical pendulum tiltmeter[J]. Seismological and Geomagnetic Observation and Research，2018，39（4）：149-153

［6］夏超德.强冷空气对宁波地震台垂直摆记录影响[J].地震地磁观测与研究，2016，37（3）：82-85

Xia C D. Influence of strong cold air on the observation of vertical pendulum tiltmeter at Ningbo Seismic Station[J]. Seismological and Geomagnetic Observation and Research，2016，37（3）：82-85

［7］杨晓东，韩晓飞，潘存英，等.韩城台模拟金属摆与数字石英摆的对比研究[J].通信电源技术，2016，33（2）：76-77

Yang X D，Han X F，Pan C Y，et al. A comparative study on the simulation of metal pendulum and digital quartz pendulum in Hancheng[J]. Telecom Power Technologies，2016，33（2）：76-77

［8］曹白伦，匡福江，段勇，等.云县台 SQ 水平摆与 VP 垂直摆观测资料对比分析[J].内陆地震，2020，34（2）：196-206

Cao B L，Kuang F J，Duan Y，et al. Comparative analysis of observational data between SQ horizontal pendulum tiltermeter and VP vertical pendulum tiltermeter at Yunxian Seismic Station[J]. Inland Earthquake，2020，34（2）：196-206

［9］张卓.海南五指山地区早二叠世（强）过铝花岗岩成因及构造意义[D].南昌：东华理工大学，2014

Zhang Z. Petrogenesis of the Early Permian（ strongly） peraluminous granites from Wuzhishan area，Hainan and its geotectonic implications[D]. Nanchang：East China University of Technology，2014

［10］中国地震局监测预报司.地形变观测[M].北京：地震出版社，2008：1-2

Monitoring and Forecasting Department of China Earthquake Administration. Terrain deformation observation[M]. Beijing： Seism- ological Press， 2008：1-2

［11］杨祺，杨玲英，方伟，等.云南地倾斜定点形变 M2波月潮汐因子精度浅析[J].地震地磁观测与研究，2020，41（4）：129-135

Yang Q，Yang L Y，Fang W，et al. Analysis of the accuracy of the M2 wave moon tidal factor using deformation observations at fixed-points in Yunnan[J]. Seismological and Geomagnetic Observation and Research，2020，41（4）：129-135

［12］郭昱琴，付国超，孙三健，等.五指山台 DSQ 水管倾斜仪和 VP 垂直摆观测数据质量对比分析[J].地震科学进展，2022，52（5）：232-237

Guo Y Q，Fu G C，Sun S J，et al. Analysis on the observation data of DSQ water pipe inclinometer and VP vertical pendulum in Wuzhishan station，Hainan Province[J]. Progress in Earthquake Sciences，2022，52（5）：232-237