五指山台水管倾斜仪北南向变化的异常判定_郭昱琴.pdf

郭昱琴，李盛，张慧，等.五指山台水管倾斜仪北南向变化的异常判定 J.地震科学进展,2023,53(2):77-85.doi:10.19987/j.dzkxjz.2022-077GuoYQ,LiS,ZhangH,etal.Abnormaldeterminationofnorth-southcomponentofwaterpipetiltmeteratWuzhishanstationJ.ProgressinEarthquakeSciences,2023,53(2):77-85.doi:10.19987/j.dzkxjz.2022-077五指山台水管倾斜仪北南向变化的异常判定*郭昱琴李盛张慧胡久常孙三健（海南省地震局，海南海口570203）摘要五指山台 DSQ 水管倾斜仪北南向自 2017 年 5 月持续南倾。排除观测系统、洞室温度、辅助观测资料情况和周围环境因素变化等多方面原因，利用观测资料内精度质量评价分析、数字化前兆异常识别等方法，结合 GPS 及重力场变化，综合认为，五指山地震台 DSQ 水管倾斜仪北南向南倾变化为正常区域背景变化趋势，非前兆异常现象。关键词水管倾斜；南倾；干扰分析；S 变换分析中图分类号：P315.63文献标识码：A文章编号:2096-7780(2023)02-0077-09doi：10.19987/j.dzkxjz.2022-077Abnormal determination of north-south component ofwater pipe tiltmeter at Wuzhishan stationGuo Yuqin,Li Sheng,Zhang Hui,Hu Jiuchang,Sun Sanjian(HainanEarthquakeAgency,HainanHaikou,570203,China)Abstract The north-south component of the DSQ water pipe tiltmeter at Wuzhishan station has been tiltingsouthwardcontinuouslysinceMay2017.WeeliminatevariousreasonssuchasthedeformationobservationsystemoftheWuzhishan station，the temperature of the cavern，the auxiliary observation data and the changes in surroundingenvironmentalfactors，andusethemethodsoftheevaluationandanalysisofthequalityofthegroundtiltobservationdataandtherecognitionofdigitalprecursoryanomalies.CombinedwithGPSandgravityfieldchanges，itisconcludedthatthesouthtiltofDSQwaterpipetiltmeteratWuzhishanstationmaybeanormalregionalbackgroundchangetrend，notaprecursoranomaly.Keywordspipetilt;southtilt;interferenceanalysis;S-transformanalysis引言形变观测信息中包含固体潮汐和非固体潮汐信息，固体潮汐是由于日、月和近地行星对地球的引力变化所导致的地球内部和表面的周期性形变1，非潮汐信息反应的则是地表倾斜和应变的变化。因此，倾斜观测可作为地震前兆定点形变观测的重要手段。前人针对形变观测进行了大量研究，陈德福和张锡令2在搜集整理了大量的倾斜异常图像后对倾斜异常形态进行了分类和特征分析；王梅等3、曹建玲和*收稿日期：2022-05-06；采用日期：2022-07-21。基金项目：海南省地震局科研课题（HDKJ202012）资助。通信作者：郭昱琴（1988），女，工程师，主要从事地震前兆监测工作。E-mail：。第53卷第2期地震科学进展Vol.53No.22023年2月ProgressinEarthquakeSciencesFeb.,2023石耀霖4论述了气压、雷电、地表温度和环境载荷等干扰对观测记录的影响及其机理。由于影响机理复杂，在地震预报实践中，如何提取较为可靠的前兆异常信息就显得尤为重要5-6。研究区域地下水环境复杂，随时间推移，由水体质量变化引起的孔隙压力扩散和渗流效应，对地形变将产生不同影响7。目前定点形变仪器种类较多，不同类型仪器的观测物理量和信号频段不同，捕捉到的前兆异常具有一定的差异性8-10。本文对五指山台水管倾斜仪北南向变化从观测系统、洞室温度、辅助观测资料情况和周围环境因素变化等多方面进行详细论证，排除前兆异常现象。将有助于对台站观测资料的理解和利用，进而为该区域地壳应力应变状态的长期变化研究提供宝贵资料。1台站基本情况海南岛以琼州海峡与华南大陆相隔，由北向南发育有王五文教、昌江琼海、尖峰吊罗和九所陵水等 4 条近东西向隐伏断裂，由西向东发育北南向戈枕断裂和白沙断裂11。五指山台位于五指山市冲山镇阿里山度假村旁太平山山麓半山腰（18.79N，109.53E）（图 1），距台站 10km 处有一条隐形平移断层，而台站正南方 25km 有尖峰吊罗深大断裂，该区域广泛出露二叠系花岗岩和三叠系花岗岩11。台站周围基岩裸露，观测硐室的台基为角闪黑云二长花岗岩，厚度较大，约 40m。基岩性坚硬完整、致密均匀，测点周边观测环境也较为稳定。台站为无人值守台站，年平均气温 22.4，年平均气压 977.9hPa，年平均降雨量 1690mm。硐室内共布设 3 套形变观测仪器，分别为 DSQ 水管倾斜仪、SSY 铟瓦棒伸缩仪和 VP 宽频带倾斜仪，均为“十五”数字化观测项目所选用的仪器。其中 DSQ 水管倾斜仪和 SSY 铟瓦棒伸缩仪布设在主洞室内，两者相同分量的仪器部件平行布设且共用同一槽体，VP 宽频带倾斜仪布设在主洞室南侧的侧洞室。气象三要素综合观测仪架设于距台站 2.4km 处的五指山市地震局，观测环境稳定。2异常判定分析2.1 观测数据异常变化情况2015 年 1 月2017 年 4 月，五指山台水管倾斜仪北南分量观测数据变化形态平稳。2017 年 5 月起，北南分量出现了持续南倾现象。从水管仪北南分量观测数据日均值（图 2）来看，往年年变幅度约700103，而 2017 年 5 月2019 年 5 月 1 日，下降幅度达 5892.796103。而同期东西分量观测数据表现正常，年变幅度在 400103范围内波动。2.2 观测系统排查观测资料出现异常后对仪器的主机、传感器、通信系统、供电线路、UPS、避雷设施等进行了检查，信号电压与工作电压均正常。对水管水质情况检查发现北南分量的钵体、水管及胶皮管均完好，水质也较为清澈透明，未发现有明显的微生物或其他杂质。对墩体进行检查，发现尽管有挡雨棚遮挡、塑料布覆盖及泡沫箱隔断，墩体依然较为潮湿，并在墩体表面伴有水滴积累。但学科专家现场查看后，认为本台洞室内和观测墩体的潮湿程度总体情况较好，对观测数据影响不大。2.3 环境干扰调查对周边环境调查发现，台站上游 1.3km 处有一大型水库太平山水库，台站西侧约 1km 处有一建筑工地水墨雨林项目。N108 E10911011120NN五指山形变台山洞洞体平面图洞体示意图铁门铁门防盗门比例尺 1:500防盗门防盗门DSQ和SSY墩体DSQ和SSY墩体船舱门船舱门山体墩体1918五指山台主要断裂高程/m2583500129图1五指山台站构造位置及仪器布设图Fig.1StructurelocationandinstrumentlayoutofWuzhishanstation78地震科学进展2023年太平山水库距五指山形变台 1.3km，位于台站东北方向，主要为城市供水，兼有防洪、旅游等功能。20 世纪 80 年代初建成，总容积约 142 万 m3，其蓄水主要来自降雨。水库观测站水文资料显示，水库水位呈冬低夏高变化。水库距台站的距离符合形变台站建设规范中倾斜、应变观测距干扰源最小“安静”距离的规定，且其库容量与五指山形变台水管倾斜仪北南分量观测数据没有明显相关性，因此，认为水管倾斜仪北南分量的缓慢南倾现象并不是水库库容量变化所导致的（图 3）。水墨雨林项目距五指山形变台站约 1km，处于台站正西方位，建筑工地位于 224 国道与甘工路交汇处附近。通过走访项目部建筑技术工程师及查阅施工日志得知，施工队于 2017 年 5 月初进驻工地，309形变量(NS向)/103形变量(EW向)/103年份年份3661701376196468201520162017201820192020201520162017201820192020图2水管倾斜仪整时值数据变化曲线（2015-01-012019-05-01）Fig.2Datachangecurveofwaterpipetiltmeter（2015-01-012019-05-01）2014-12-022015-12-022016-12-012017-12-012018-12-0104080120水库库容量/104 m3日期图3太平山水库库容量年变化Fig.3AnnualvariationofTaipingshanreservoircapacity第2期郭昱琴等：五指山台水管倾斜仪北南向变化的异常判定79开始打桩，并于 2017 年 712 月完成大楼主体建设。水墨雨林项目共有 5 栋楼在建，楼高为 89 层。根据此项目的方位（台站西侧）、规模以及距台站距离等判断，均不足以对仪器产生影响（图 4）。图4水墨雨林项目工地图片Fig.4ThepictureoftheInkRainForestProjectSite2.4 气候因素调查五指山台气象三要素观测仪 2017 年 9 月 26 日投入运行，对比分析 2017 年 9 月以来降雨量和水管仪北南向观测数据，发现降雨对水管仪观测数据无明显影响，南倾变化与降水不存在相关性（图 5）。334降雨量/mm2.201.10038377340201520162017年份201820192020形变量(NS向)/103图5五指山台降水量（蓝色实线）和水管仪北南分量观测曲线（黑色实线）Fig.5Precipitation（blue solid line）and NS componentchangecurveofwaterpipetiltmeter（blacksolidline）atWuzhishanstation2.5 观测数据的可靠性分析现 有 倾 斜 观 测 资 料 调 和 分 析 精 度 指 标 是 按Venedikov 调和分析方法，对月尺度的倾斜观测资料进行调和分析，以 M2波的中误差 m指标作为评定倾斜观测资料质量的第一项指标，规范要求 类倾斜台站观测 m0.02。根据中国地震局地壳形变学科组公布的月评内精度结果，统计 2015 年 1 月2019 年 5 月五指山台定点形变观测内精度可以看出：水管仪北南向自 2016 年 5 月后观测资料内精度满足规定优秀指标，观测质量较好。2017 年 5 月出现趋势转折后，M2波潮汐因子内精度未出现异常变化（图 6）。2.6 数字化前兆异常识别方法分析2.6.1 非潮汐变化分析方法日均值方法是固体潮数据非潮汐分析方法之一，即对一天的 24 个整点值取平均，它可以滤去97%的潮汐成分，剩下的是零点飘移和非潮汐信息12。从水管倾斜仪北南分量日均值及其差分，5 日均值及其差分发现异常突出均为仪器校准检修影响（图 7）。2.6.2 超限率方法利用小波分解和改进的超限率可以得到的某个带通信号，可以分为正常背景和异常两部分。通过对这两部分界线划分的研究，找到了客观、唯一确定超限率分析的阈值方法13。结合五指山台水管倾斜仪北南分量长期（4 年时长）的超限率计算结果，发现超限率数量具有显著增强的均为仪器故障及洞室施工影响造成的。影响时间为：2016 年 5 月 10 日至6 月 2 日前置放大器故障、2016 年 7 月份五指山台洞室改造在水管倾斜仪北端建筑施工、2017 年 6 月5 日至 22 日传感器故障、2019 年 3 月 11 日更换传感器原因影响（图 8）。2.6.3 S 变换方法变换方法属于时频分析类方法的一种，可以80地震科学进展2023年将一维的观测数据时间序列转换成二维的时频域信号，从而增加获取的信息量，为研究数据中各类频率信号的频谱动态演化提供便利。通过 S 变换方法曾有效地提取了姑咱台四分量钻孔应变在芦山 7.0 级地震前的高频应变畸变信号14。我们尝试对五指山台水管倾斜仪北南分量进行 S 变换，时频变换结果发现，在趋势转折时未显示高频畸变信号（图 9）。3GPS 及重力场结果3.1 GPS 应力场从单个 GPS 站点数据（图 10）来看，海南三亚（HISY）站点自 2015 年 11 月份以来呈现向西偏南方201520162017年份2018201900.020.040.060.080.10中误差 m图6五指山台水管仪北南向内精度结果Fig.6NScomponentprecisioncurveofWuzhishanwaterpipetiltmeter334年份36677000201520162017201820192020340年份36707000201520162017201820192020120日均值一阶差分5日均值一阶差分5日均值日均值年份水管加水水管加水129620152016201720182019184186020101151年份732972015201620172018201920202020形变量/103形变量/103一阶差分/103一阶差分/103图7五指山台水管倾斜仪北南分量差分曲线Fig.7NScomponentdifferencecurveofwaterpipetiltmeteratWuzhishanstation第2期郭昱琴等：五指山台水管倾斜仪北南向变化的异常判定81向偏移；海南琼中（QION）站点自 2010 年 6 月份以来呈现向北偏西方向偏移；海南海口（HIHK）站点则一直较为平稳。从 GPS 站点间基线数据来看，海南海口（HIHK）海南琼中（QION）虽表现为缓慢缩短挤压状态，但海南海口（HIHK）海南三亚（HISY）、海南琼中（QION）海南三亚（HISY）则均表现为快速拉伸状态，这表明整个海南岛主要处于一个地壳拉伸状态。地壳拉伸可能会引起地壳出现下降（沉），而不均匀下降可以引起地壳南倾，从而使五指山台长基线水管仪数据出现南倾变化现象。3.2 重力场资料分析从海南岛陆重力场动态的变化（图 11）来看：2015 年 9 月2016 年 9 月，海南岛陆呈现东北部正变化，南部和西部负变化的过程，零值线穿琼东南而过，形成重力场正负异常变化高梯度区；而 2016 年9 月2017 年 10 月，海南岛陆变化比较平稳，表现为南部正变化，北部负变化。五指山台位于海南南部，自 2015 年 9 月以来，该地区的重力场变化经历了先负后正的变化过程。研究结果表明，重力值减小，说明地表隆升或地下物质减少或是两者迭加效应；重力值增大，说明地表沉降或地下物质增多或是其迭加效应15。比较琼中台重力非潮汐变化和 GPS站基线变化，显示琼中台重力非潮汐变化 2017 年1 月以来处于小幅上升变化趋势，琼中 GPS 基线垂直分量年中开始处于缓慢下降过程，两者的变化均表明 2017 年以来琼中台地区地表小幅沉降；而同时间段 GPS 基线琼中三亚呈伸长趋势；结合这些变化特征，认为 2017 年以来海南南部地区均处于地表沉降的过程，且北部的沉降幅度小于南部的沉降幅度。以上也可佐证五指山台基线式水管仪因区域构造背景变化导致数据出现持续南倾变化。4综合分析及讨论通过对五指山形变台观测系统、洞室温度、辅助观测资料情况和周围环境因素变化等多方面调查分析，排除了这些影响因素导致水管倾斜仪北南向出60040020001000500数量超限率/个0201520162017年份20182019201520162017年份201820192.5倍均方差阈值强度超限率/1032497%1357图8五指山台水管倾斜仪北南分量超限率分析结果Fig.8AnalysisresultofNScomponentexceedinglimitrateofwaterpipetiltmeteratWuzhishanstation4213S变换幅值/1030201520162017年份20182019最大幅度平均幅度710周期信号/min幅值/10317603210图9五指山台水管倾斜仪北南分量 S 变换结果Fig.9S-transformresultsofNScomponentofwaterpipetiltmeteratWuzhishanstation82地震科学进展2023年现异常变化的可能。根据中国地震局地壳形变学科组公布的月评结果统计发现，五指山台水管仪北南向自 2016 年 6 月起内精度满足学科优秀指标，观测资料质量较好。且 2017 年 5 月出现趋势转折后，M2波潮汐因子内精度平稳，未出现异常变化。通过数字化前兆异常识别方法分析，未发现超出阈值及高频信号异常信息。通过区域 GPS 基线数据分析表明，海口至琼中201010（a）（b）（c）（e）（d）位移/mm0102011 2012 2013 2014N向年份年份年份年份年份年份年份年份年份年份年份年份HIHK点位时间序列(去线性)2015 2016 2017 2018 2019201010位移/mm0102011 2012 2013 2014QION点位时间序列(去线性)2015 2016 2017 2018 2019201010位移/mm0102011 2012 2013 2014E向2015 2016 2017 2018 2019201010位移/mm0102011 2012 2013 2014 2015 2016 2017 2018 2019201030位移/mm0302011 2012 2013 2014U向N向E向U向2015 2016 2017 2018 2019201030位移/mm0302011 2012 2013 2014 2015 2016 2017 2018 2019201010位移/mm0102011 2012 2013 2014N向HISY点位时间序列(去线性)2015 2016 2017 2018 201920101020位移/mm0102011 2012 2013 2014HIHK_HISYHIHK_QIONQION_HISY2015 2016 2017 2018 2019201010位移/mm0102011 2012 2013 2014E向2015 2016 2017 2018 2019201020位移/mm0202011 2012 2013 2014 2015 2016 2017 2018 2019201030位移/mm0302011 2012 2013 2014108 E22 N2018GNSSQIONHISYHIHK站点112U向2015 2016 2017 2018 2019201025位移/mm51552011 2012 2013 2014 2015 2016 2017 2018 2019图10GNSS 基线变化曲线（a，b，c，d）及站点分布图（e）Fig.10GNSSbaselinechangecurve（a，b，c，d）andsitedistributionmap（e）第2期郭昱琴等：五指山台水管倾斜仪北南向变化的异常判定83地区地壳处于轻微挤压状态，而琼中至三亚地区地壳处于拉伸状态，整个海南岛区域受到这样的不均衡应力变化，易产生地壳的不均衡沉降，水管倾斜仪北南分量持续南倾可能是地壳形变的真实反映。海南岛陆重力场变化显示，五指山形变台附近区域出现重力异常变化，认为 2017 年以来海南南部地区均处于地表沉降的过程，且北部的沉降幅度小于南部的沉降幅度，此结论与 GPS 资料结果呼应。综合分析认为，五指山台水管倾斜仪北南向变化为正常区域背景变化趋势，非地震前兆异常现象。水管倾斜仪是通过检测和记录两端钵体液面的变化来模拟地壳的实际运动状态的重要观测手段，做好监测水管倾斜仪的变化趋势与跟踪分析，对琼东南地区的震情跟踪工作具有一定帮助。参考文献许厚泽.固体地球潮汐 M.武汉：湖北科学技术出版社，2010：1-2XuHZ.SolidearthtidesM.Wuhan：HubeiScienceandTechnologyPress，2010：1-21陈德福，张锡令.气流扰动对自记水管倾斜仪影响的实验结果 J.地壳形变与地震，1987，7(4)：288-296ChenDF，ZhangXL.Experimentalresultsoftheinfluenceofairflowdisturbanceontheself-recordingwaterpipeinclinometerJ.CrustalDeformationandEarthquake，1987，7(4)：288-2962王梅，李峰，孔向阳.数字化形变观测干扰识别 J.大地测量与地球动力学，2004，24(1)：94-98WangM，LiF，KongXY.IdentificationofdisturbanceofdigitaldeformationobservationsJ.JournalofGeodesyandGeodynamics，2004，24(1)：94-983曹建玲，石耀霖.地表温度年变化对地应力和地倾斜的影响 J.中国科学院研究生院学报，2005(3)：303-308CaoJL，ShiYL.StressandtiltinducedbyannualvariationofsurfacetemperatureJ.JournaloftheGraduateSchooloftheChineseAcademyofSciences，2005(3)：303-3084海南省地震局监测处.地震观测资料评比办法及管理细则汇编 S.海南省地震局，2014：9-10MonitoringOfficeofHainanEarthquakeAgency.CompilationofseismicobservationdataappraisalmethodsandmanagementrulesS.HainanEarthquakeAgency，2014：9-105中国地震局监测预报司.地震前兆异常落实工作指南 M.北京：地震出版社，2000：25-31Monitoring and Forecasting Department of China Earthquake Administration.Guidelines for the implementation of earthquakeprecursoranomaliesM.Beijing：SeismologicalPress，2000：25-316朱冰清，赵黎明，张磊，等.水库库容变化对蓟县地震台 DSQ 型水管倾斜仪观测的影响 J.地震地磁观测与研究，2019，40(2)：108-114ZhuBQ，ZhaoLM，ZhangL，etal.EffectofreservoircapacitychangeonobservationofDSQtypewaterpipeinclinometeratJixianSeismicStationJ.SeismologicalandGeomagneticObservationandResearch，2019，40(2)：108-1147龚丽文，陈丽娟，吕品姬，等.黔江台水管仪与垂直摆预报效能对比分析 J.地震，2021，41(4)：168-179GongLW，ChenLJ，LPJ，etal.Comparativeanalysisofthepredictionefficiencybetweenwater-tubeandverticalpendulumtiltmetersinQianjiangstationJ.Earthquake，2021，41(4)：168-1798海口海口琼中2015-092016-09三亚三亚琼中2016-092017-10-40-40-40-40-20-20-20-2000000020202040108.6 E 109.0 109.4109.8110.2 110.6 111.0108.6 E 109.0109.4109.8 110.2 110.6 111.018.218.619.019.419.820.2N18.218.619.019.419.820.2N乌石九所牛漏那大那大永发永发白马井白马井大致坡大致坡1009080706050403020100102030405060708090100Gal1009080706050403020100102030405060708090100Gal-20-2000000乌石九所牛漏图11海南岛陆重力场动态变化Fig.11DynamicchangeoflandgravityfieldinHainanIsland84地震科学进展2023年郭昱琴，付国超，孙三健，等.五指山台 DSQ 水管倾斜仪和 VP 垂直摆观测数据质量对比分析 J.地震科学进展，2022，52(5)：232-237GuoYQ，FuGC，SunSJ，etal.AnalysisontheobservationdataofDSQwaterpipeinclinometerandVPverticalpenduluminWuzhishanStation,HainanProvinceJ.ProgressinEarthquakeSciences，2022，52(5)：232-2379胡炜，张小燕，藤赫.不同地震台站水管倾斜仪数据资料对比分析以呼和浩特和包头台为例 J.地震科学进展，2020，50(11)：15-20HuW，ZhangXY，TengH.Comparativeanalysisofwaterpipeinclinometerdataofdifferentseismicstations：TakingHohhotandBaotouasexamplesJ.ProgressinEarthquakeSciences，2020，50(11)：15-2010张卓.海南五指山地区早二叠世（强）过铝花岗岩成因及构造意义 D.南昌：东华理工大学，2014：9-10Zhang Z.Genesis and tectonic significance of Early Permian(strong)peraluminous granites in Wuzhishan area，HainanD.Nanchang：DonghuaUniversityofScienceandTechnology，2014：9-1011中国地震局监测预报司.地形变观测 M.北京：地震出版社，2008，320-343Monitoring and Forecasting Department of China Earthquake Administration.Crustal deformation measurementM.Beijing：SeismologicalPress，2008：320-34312邱泽华，唐磊，张宝红，等.用小波超限率分析提取宁陕台汶川地震体应变异常 J.地球物理学报，2012，55(2)：538-546QiuZH，TangL，ZhangBH，etal.ExtractinganomalyoftheWenchuanEarthquakefromthedilatometerrecordingatNSHbymeansofWavelet-OverrunRateAnalysisJ.ChineseJournalofGeophysics，2012，55(2)：538-54613刘琦，张晶，晏锐，等.高采样率四分量钻孔应变同震响应分析 J.中国地震，2013，29(1)：57-67Liu Q，Zhang J，Yan R，et al.The analysis of coseismic response of 4-component borehole strainmeters with high samplingrateJ.EarthquakeResearchinChina，2013，29(1)：57-6714李辉，申重阳，孙少安，等.中国大陆近期重力场动态变化图像 J.大地测量与地球动力学，2009，29(3)：1-10Li H，Shen C Y，Sun S A，et al.Dynamic gravity change in recent years in China continentJ.Journal of Geodesy andGeodynamics，2009，29(3)：1-1015第2期郭昱琴等：五指山台水管倾斜仪北南向变化的异常判定85