广东形变观测历经了模拟和数字长时间的数据观测：期间模拟观测经过近30年的运行，在地球动力学及地震预报研究中均显示出不可替代的作用。N3水准仪、007水准仪、摆式倾斜仪、电子倾斜仪、水管倾斜仪、金属摆倾斜仪、整体倾斜仪等，水准仪的数据所显示的各测区的总体形变过程及局部形变变化特征，为全省地震趋势会商会提供了背景资料。倾斜仪能观测到固体潮，这些仪器可以观测到目前唯一能从理论上精确计算的形变（倾斜、重力、应变）固体潮，还有可能捕捉到不同形式的地震前兆信息。

广东省形变数字化观测自2007年陆续开始至今，也经历了10多年的时间。形变仪目前主要为地倾斜和地应变仪，包括“九五”、“十五”地倾斜5套（垂直摆倾斜仪和水管倾斜仪）和地应变仪3套（伸缩仪和分量式钻孔应变仪）。其中汕头4套，分别为DSQ水管倾斜仪、VP宽频带倾斜仪、SS-Y铟瓦棒伸缩仪和RZB-2分量式钻孔应变仪；信宜2套为VS垂直摆倾斜仪和RZB-2分量式钻孔应变仪；潮州1套VP宽频带倾斜仪；韶关1套VS垂直摆倾斜仪。形变数字化前兆仪器的使用，避免了观测人员到观测室换纸等干扰因素，减少了对观测环境的影响，同时也避免了人工读数造成的误差。因此观测资料总体质量应更有保证。

广东地区虽然是少震区，但在过去模拟时期的30多年的地震预报探索过程中。也观测到一些震前形变异常，如1986年的阳江5.0级地震和自1970年以来新丰江水库区4级以上地震。信宜台整体摆倾斜仪观测结果显示出倾斜量有异常变化；新丰江台的形变仪在新丰江水库区4级以上地震前后，有异常变化显示。通过岩石力学实验和理论研究以及大量的观测和震例总结，以往地形变模拟观测在地震预报中的作用已为地震学者所公认。2007年以后投入的数字化观测，资料已应用到地震监测中。在积累了多年的资料后，在做跟踪分析之余，有必要根据形变观测规范的标准，对观测数据进行全面的认识和评价，以便能更有效的在强地震的监测中发挥作用^[2]。

在整理多年的资料时发现，部分台站的可用资料长度仅3~4个月，计算m_γ、s_a的误差太大，因此用月序列主波潮汐参数的稳定性研究资料的稳定性。另外在观测精度分析前增加了Nakai检验，对资料进行预处理，并以残差值的分布评判资料稳定性。

观测数据内精度分析对数据的基本要求是，①进行Nakai检验的资料连续率在48 h以上，并剔除不连续的数据段；②为保证潮汐参数的准确性，必须排除Nakai检验中拟合残差太大的数据段。原计划统一使用2017年1月到2017年12月的资料进行分析，但在整理数据后发现，部分台站不连续的数据段较多，如韶关台的垂直摆倾斜仪、潮州台的宽频带倾斜仪等，为保证计算精度和准确性，选择至少连续稳定观测3个月以上的资料参加分析。各套仪器所选取的时间见表1、表2。

连续形变主要观测倾斜与应变。本文根据中国地震局制定的《地震及前兆数字观测技术规范—地壳形变观测（试行）》^[3]及中国地震局监测预报司编制的《地形变测量》为标准，以2017年观测数据为例，对广东省现行在观测的各形变仪进行分析。

《地形变测量》中提出，现有的倾斜和应变观测资料内精度指标共有两项内容：观测资料调和分析，观测相对噪声水平计算。而优秀台站即Ⅰ类台的标准是，地倾斜M2波中误差m_γ≤0.02、日均值相对噪声水平M₁≤0.02"；地应变M2波相对中误差s_a/a≤0.05、日均值相对噪声水平M_l≤0.05× 10^{-6 [4]}。

1.3.1 固体潮观测精度

（1）利用相对M2波潮汐因子的中误差和相对误差。检验观测资料的精度。优秀台站（Ⅰ类台）的标准是，地倾斜潮汐因子的中误差m_γ≤0.02。应变相对中误差s_a/a≤0.05。

（2）内在精度计算结果分析。使用中国地震前兆台网数据处理系统软件对各台站数字化形变观测资料进行了数据检验和评价分析。采用以倾斜潮汐因子γ值、线应变因子α等为代表的、既具有明确物理意义，又有较高的有效率的误差分析方法进行计算^[5]。地倾斜仪的倾斜潮汐M2波γ值中误差以及应变仪潮汐M2波σ值相对中误差计算精度结果见表1、表2，由表可见各台形变仪精度大多数达到优秀标准。

表1 地倾斜观测数据精度结果表Table 1 Precision results of clinometer observation data

表2 地应变观测数据精度结果表Table 2 Precision results of strainmeter observation data

各台形变仪误差结果见图1、图2。对地倾斜仪资料的分析表明，①水管倾斜仪的误差最小， m_γ≤0.02的月份有9个，整年误差均值是0.0214，最小误差值为0.0077；②同为VP宽频带倾斜仪，汕头台的误差相对较小， NS向较EW向好， m_γ≤0.02的月份NS向有8个，EW向无，整年误差均值是0.0323，最小误差值为0.0128。潮州台为市县台，从安装开始就不稳定，其误差较大，尤其EW向更甚， m_γ≤0.02的月份NS向有6个， EW向无，整年误差均值是0.0573，最小误差值为0.0148；VS垂直摆倾斜仪，韶关台的误差较信宜台大，可能是因为仪器老化，故障较多，导致2台的误差均较大。韶关台NS向较EW向差， m_γ≤0.02的月份NS向有5个，EW向11个，整年误差均值是0.0276，最小误差值为0.0081。信宜台主要是NS向较差， EW向很好， m_γ≤0.02的月份NS向有6个， EW向12个，整年误差均值是0.0167，最小误差值为EW向的0.00。

对洞体应变仪资料的分析来看，相对误差较小， s_a/a≤0.05的月份12个，整年误差均值是0.0268，最小误差值为0.007。说明潮汐观测的日变形态正常，稳定性好。

钻孔体应变资料的分析表明，汕头和信宜2个台站的相对误差均很小。s_a/a≤0.05的月份分别为11和12个，整年误差均值是0.0183和0.0131，最小误差值为0.0054、0.0024。

图1 各台倾斜仪M2波γ值中误差Fig.1 The RMS of γ value of M2 wave of each inclinometer

图2 各台应变仪M2波σ值相对中误差Fig.2 The relative RMS of σ value of M2 wave of each strain gauge

1.3.2 相对噪声水平

相对噪声水平M₁是用来判断观测资料长期稳定性的一项定量指标。优秀台站即Ⅰ类台此项标准是，地倾斜日均值相对噪声水平M₁≤0.02"；地应变日均值相对噪声水平M_l≤0.05×10^-6。采用中国地震前兆台网数据处理系统集成版软件提供的相对噪声水平计算方法，计算了广东省4个台站8套形变仪器连续3年的数据均值。计算结果见表3。结果表明，除市县管理的潮州台VP宽频带倾斜仪相对噪声水平较高，为0.0308"和0.0568"未达到优秀外，其余均为优秀以上，其中伸缩仪的最低，其次是2套钻孔应变仪。

表3 形变仪相对噪声水平表Table 3 Relative noise level of deformation instruments

数据连续是地震预报对观测数据的基本要求之一。对2017年度各台仪器的连续率进行统计，结果见表4可以看出，各仪器的连续率均达到99%以上，观测仪器的运行状态良好。

由漂移量可以看出仪器的稳定性和年变形态，漂移量是检验仪器是否正常工作的重要指标。由表5和图3可以看出：对于地倾斜仪，汕头台的1套水管倾斜仪的月漂移量在0.001～0.08之间；汕头台、潮州台、信宜台和韶关台的垂直摆倾斜仪，其NS向月漂移量在0.0002～0.61之间， EW向在0.0015～0.52之间。总体来说，垂直摆倾斜仪的月漂移量都不大，多数在0.2内。信宜和韶关台的月漂移量较大，信宜台NS向超出0.2范围内有3个，韶关1个，EW向超出0.2信宜、韶关台各1个，说明这2台仪器的稳定性较其他台略差。

对于地应变仪，汕头台的洞体应变仪的稳定良好，月漂移量在0.0003～0.03之间；信宜台、汕头台的钻孔应变仪NS向和EW向月漂移量，除了汕头台各有1个较大值外，分别为2.8和0.6，其余均在0.1范围内。

表4 2017年度各台形变仪连续率Table 4 The continuity rate of each deformation instrument in 2017

表5 2017年各台形变仪年零漂均值Table 5 Mean value of annual zero drift of deformation instruments in 2017

噪声水平指标与观测环境、观测系统噪声水平等直接关联，主要反映数据动态变化的稳定性水平^[6-7]。从表3结果来看，汕头台形变观测仪器安装在168 m的坑道内，干扰少，观测条件较好，各仪的相对噪声水平均较低。信宜台洞壁为花岗片麻岩，岩性完好。2套形变仪的相对噪声水平应变仪较小，垂直摆稍大。韶关台基条件稳定，观测条件较理想，形变仪噪声水平较低。潮州台噪声水平较高，该台是隶属潮州市地震局的地方台站。形变观测的山洞洞壁和洞顶为石砌结构，底部为水泥埕面，洞体长年潮湿，加上南方雷雨天气频繁，因而一定程度影响了记录的连续性、完整性、稳定性和观测质量^[8-9]。

图3 2017年倾斜、应变仪各月零漂值Fig.3 Monthly zero drift value of inclinometer and strainometer in 2017

潮汐变化是外部天体对于地面观测站相对位置变化造成的附加引力而引起的，而地震研究关注的是从实际观测资料中分离出与地震孕育有关的非潮汐因素所引起的倾斜、应变及重力变化^[10]。用Venedikov法可将潮汐波背景值剔除，较完整地保留细微的形变信息。M2波潮汐因子中误差是评定观测资料精度的主要标准，用相位滞后可以比较不同台站资料的稳定程度^[11-12]。在地震平静期，用上述方法也可反映出观测系统受干扰情况、自身运行情况及台址条件，从而对各台资料稳定性和可靠性进行评价

2.3.1 地倾斜仪

洞体地倾斜测量是在山洞或地下室内应用倾斜仪测定地面的倾斜变化。广东此类形变观测仪有水管倾斜仪和垂直摆倾斜仪。

倾斜仪γ因子中误差和M2波相位滞后误差见图1和图4。对DSQ水管倾斜仪潮汐因子中误差及相位滞后误差计算情况来看：γ因子中误差达到优秀标准的月份NS向10个较EW向7个多，说明观测精度NS向较EW向高。汕头、信宜VP宽频带倾斜仪NS向γ因子中误差达到优秀标准的月份较EW向多。信宜、韶关垂直摆倾斜仪则EW向γ因子中误差达到优秀标准的月份较NS向多，误差小的分量说明观测精度较高。

从M2波相位滞后的误差可以得出：DSQ水管倾斜仪、汕头VP宽频带倾斜仪和韶关VS垂直摆倾斜仪的EW向误差较小，表明固体潮波动小，数据记录稳定。而潮州VP宽频带倾斜仪和韶关VS垂直摆倾斜仪的NS向误差相对较大，表明固体潮波动大，数据不够稳定。

垂直摆倾斜仪m_γ年误差均值和M2波相位滞后的误差由小到大的排序为信宜、韶关、汕头、潮州和汕头、信宜、韶关、潮州。信宜台的m_γ误差最小，相位误差也较小，说明观测精度高，资料较为稳定。汕头台相位误差最小，m_γ误差略大，潮州台各项指标均较差，观测精度及稳定性均不理想。

图4 形变仪相位滞后误差Fig.4 Phase lag error of deformation instrument