地震研究

针对gPhone连续重力仪漂移修正困难的情况,收集2019年1月—2024年5月云龙地震台gPhone连续重力仪观测数据,采用流动重力数据控制的方法进行漂移修正,并与目前常用的小波滤波和多项式拟合修正方法进行比较,修正后的结果显示:① 5.4 a累计变化量从约1 374×10^-8 m/s²降至69×10^-8 m/s²; ②漂移修正结果与5阶多项式修正结果相似; ③模型评价结果优于其他2种方法; ④残差分布与5阶多项式修正结果接近。

It is difficult to correct the drift of the data produced by the continuous gravimeter,gPhone.To solve this problem,we collect the data measured by gPhone at Yunlong Station from January 2019 to May 2024,and do the drift correction through mobile gravity data control.Then we compare the de-drifted data with the data corrected with the wavelet filtering method and the polynomial fitting method.We find that after drift correction through mobile gravity data control,the gravity data's cumulative change in 5.4 years decreases from about 1 374 × 10^-8 m/s² to 69 × 10^-8 m/s².The de-drifted data are similar to the corrected data processed through 5th-order polynomial of the wavelet filtering method.The residual distribution through mobile gravity data control is close to the one of the data processed through 5th-order polynomial correction of the wavelet filtering method.The method of the drift correction through mobile gravity data control is evaluated to be superior to the other two methods.The drift correction method we proposed is practical,and it is a new approach for the drift correction of the continuous gravity data.

引言
1 台站及仪器简介
2 数据与处理方法
3 漂移修正方法比较
4 结论

图1 云龙台周边地区地质构造图<br/>Fig.1 Geological structures in the region whereYunlong Station is located

图1 云龙台周边地区地质构造图
Fig.1 Geological structures in the region whereYunlong Station is located

图2 云龙台重力M2波潮汐因子变化时间序列<br/>Fig.2 The time series of M2 wave tidal factors of gravity at Yunlong Station

图2 云龙台重力M2波潮汐因子变化时间序列
Fig.2 The time series of M2 wave tidal factors of gravity at Yunlong Station

图3 云龙台重力原始值及预处理值时间序列<br/>Fig.3 The time series of the original minute valuesof gravity recorded by Yunlong Station andtheir preprocessed ones

图3 云龙台重力原始值及预处理值时间序列
Fig.3 The time series of the original minute valuesof gravity recorded by Yunlong Station andtheir preprocessed ones

图4 云龙台与2个绝对重力台站联测<br/>Fig.4 The joint observation network composed ofYunlong Station and other two absolute gravitystations:Xiaguan and Lijiang

图4 云龙台与2个绝对重力台站联测
Fig.4 The joint observation network composed ofYunlong Station and other two absolute gravitystations:Xiaguan and Lijiang

表1 流动重力观测变化量与同时段连续重力观测变化量及差值<br/>Tab.1 The variation of the mobile gravity observation and the variation of the continuous gravity observation in the same period and the difference between them

表1 流动重力观测变化量与同时段连续重力观测变化量及差值
Tab.1 The variation of the mobile gravity observation and the variation of the continuous gravity observation in the same period and the difference between them

图5 基于流动重力控制的漂移修正分钟值时间序列<br/>Fig.5 The time series of minute values after driftcorrection based on mobile gravity control

图5 基于流动重力控制的漂移修正分钟值时间序列
Fig.5 The time series of minute values after driftcorrection based on mobile gravity control

图6 各漂移修正方法结果对比<br/>Fig.6 Comparison of the results from differentdrift correction methods

图6 各漂移修正方法结果对比
Fig.6 Comparison of the results from differentdrift correction methods

图7 去漂移后数据加回理论固体潮结果<br/>Fig.7 The theoretical solid tide after thede-drifted residuals were added back

图7 去漂移后数据加回理论固体潮结果
Fig.7 The theoretical solid tide after thede-drifted residuals were added back

表2 3种方法漂移修正后加回理论固体潮精度评定结果<br/>Tab.2 Evaluation of the theoretical solid tide after the residualsde-drifted with three methods were added back

表2 3种方法漂移修正后加回理论固体潮精度评定结果
Tab.2 Evaluation of the theoretical solid tide after the residualsde-drifted with three methods were added back

图8 各漂移修正方法残差频率分布<br/>Fig.8 Distribution of residuals coming from various drift-correction methods

图8 各漂移修正方法残差频率分布
Fig.8 Distribution of residuals coming from various drift-correction methods

表3 3种方法漂移修正结果的残差统计指标<br/>Tab.3 Statistical indicators of residuals of the resultscorrected through three drift-correction methods

表3 3种方法漂移修正结果的残差统计指标
Tab.3 Statistical indicators of residuals of the resultscorrected through three drift-correction methods

窦喜英,韦进,王恩利,等.2019.双重约束条件下的连续重力观测非线性漂移改正方法分析[J].地震工程学报,41(5):1251-1258.
Dou X Y,Wei J,Wang E L,et al.2019.Nonlinear drift correction method for continuous gravity observation under dual constraints[J].China Earthquake Engineering Journal,41(5):1251-1258.(in Chinese)
李响.2024.基于大数据平台的连续重力数据算法模块实现[D].三河:防灾科技学院.
Li X.2024.Implementation of continuous gravity data algorithm module based on big data platform[D].Sanhe:Institute of Disaster Prevention.(in Chinese)
李永生,周晨,张思萌,等.2025.利用Prophet模型进行地下水位异常识别初探[J].大地测量与地球动力学,45(3):313-318.
Li Y S,Zhou C,Zhang S M,et al.2025.Preliminary exploration of using the Prophet model for identifying groundwater level anomalies[J].Journal of Geodesy and Geodynamics,45(3):313-318.(in Chinese)
刘东,郝洪涛,王青华,等.2021.2021年云南漾濞MS6.4地震前重力变化[J].地震地质,43(5):1157-1170.
Liu D,Hao H T,Wang Q H,et al.2021.Gravity field change before the 2021 Yangbi MS6.4 earthquake,Yunnan[J].Seismology and Geology,43(5):1157-1170.(in Chinese)
马娟娟,王乐行,江颖,等.2024.利用回归模型对gPhone重力仪漂移改正研究[J/OL].(2024-10-29)[2024-11-07].大地测量与地球动力学,https://doi.org/10.14075/j.jgg.2024.09.419.
Ma J J,Wang L X,Jiang Y,et al.2024.Study on gPhone gravimeter drift correction using regression Model[J/OL].(2024-10-29)[2024-11-07].Journal of Geodesy and Geodynamics,https://doi.org/10.14075/j.jgg.2024.09.419.(in Chinese)
王乐行,刘子维,江颖,等.2025.gPhone重力仪漂移改正方法研究[J].大地测量与地球动力学,45(1):104-110.
Wang L X,Liu Z W,Jiang Y,et al.2024.Research on gPhone Gravimeter Drift Correction Method[J].Journal of Geodesy and Geodynamics,45(1):104-110.(in Chinese)
王林松,陈超,梁青,等.2012.东日本大地震重力信号同震响应[J].武汉大学学报(信息科学版),37(11):1348-1351.
Wang L S,Chen C,Liang Q,et al.2012.Coseismic response to gravity signals of the great earthquake in east Japan[J].Geomatics and Information Science of Wuhan University,37(11):1348-1351.(in Chinese)
韦进,刘高川,李辉,等.2012.弹簧式连续重力观测非线性漂移影响因素分析[J].大地测量与地球动力学,32(5):137-142.
Wei J,Liu G C,Li H,et al.2012.Analysis of non-linear drift of spring gravimeter[J].Journal of Geodesy and Geodynamics,32(5):137-142.(in Chinese)
韦进,刘子维,郝洪涛,等.2011.日本MS9.0地震前的连续重力观测异常[J].大地测量与地球动力学,31(2):9-11.
Wei J,Liu Z W,Hao H T,et al.2011.Gravity anomalies before MS9.0 earthquake in Japan observed by continuous relative gravimeter in China[J].Journal of Geodesy and Geodynamics,31(2):9-11.(in Chinese)
邢乐林,李辉,夏正超,等.2010.CG-5重力仪零漂特性研究[J].地震学报,32(3):369-373.
Xing L L,Li H,Xia Z C,et al.2010.Study on zero drift characteristics of CG-5 gravimeter[J].Acta Seismologica Sinica,32(3):369-373.(in Chinese)
羊锴,韦进.2015.基于连续重力观测数据应用平台的数据预处理[J].中国科技信息,11:72-73.
Yang K,Wei J.2015.The date preprocessing of application platform based on continuous gravity observation data[J].China Science and Technology Information,11:72-73.(in Chinese)
杨学慧,常玉巧,李燕玲,等.2016.用同相叠加方法识别定点形变数据中的降雨影响[J].大地测量与地球动力学,36(10):936-940.
Yang X H,Chang Y Q,Li Y L,et al.2016.Distinguishing the rainfall disturbing pattern in crustal deformation recording by in-phase stack of data[J].Journal of Geodesy and Geodynamics,36(10):936-940.(in Chinese)
周江林,沈萍,田鑫.2015.北京地震台gPhone重力仪同震响应特征分析[J].中国地震,31(3):553-561.
Zhou J L,Shen P,Tian X.2015.Analysis of the characteristics of co-seismic response of gPhone gravimeter at Beijing seismic station[J].Earthquake Research in China,31(3):553-561.(in Chinese)
周鹏飞,王艳霞,周汝良.2024.基于ICESat-2的横断山区NASA DEM高程精度评价与误差校正[J].西南林业大学学报(自然科学),44(5):136-146.
Zhou P F,Wang Y X,Zhou R L,et al.2024.Evaluation of NASA DEM elevation accuracy and error correction in the Hengduan Mountain region based on ICESat-2[J].Journal of Southwest Forestry University(Natural Science),44(5):136-146.(in Chinese)

备注

引言

1 台站及仪器简介

2 数据与处理方法

3 漂移修正方法比较

4 结论

期刊信息