快速检索        
  武汉大学学报·信息科学版  2015, Vol. 40 Issue (12): 1683-1689

文章信息

贾剑钢, 栾威, 申文斌
JIA Jiangang, LUAN Wei, SHEN Wenbin
iGrav-007超导重力仪的性能分析及对球型自由振荡模态0Sm的检测
iGrav-007 SG and Detection of the Spherical Free Oscillation Modes 0Sm
武汉大学学报·信息科学版, 2015, 40(12): 1683-1689
Geomatics and Information Science of Wuhan University, 2015, 40(12): 1683-1689
http://dx.doi.org/10.13203/j.whugis20140145

文章历史

收稿日期: 2014-02-07

iGrav-007超导重力仪的性能分析及对球型自由振荡模态0Sm的检测
贾剑钢 , 栾威, 申文斌    
武汉大学测绘学院, 湖北 武汉, 430079
摘要: iGrav超导重力仪是世界上最新型的便携式高精度重力仪,可提供目前为止最稳定的相对重力测量。基于Tsoft软件处理iGrav-007在武汉九峰台站连续5个月的观测数据,进行气压负荷、极移重力效应改正,分析重力残差和漂移等特性。通过频谱分析,检测观测期间2013年Mw6.6级芦山地震激发的球型自由振荡模态0Sm,检测结果明显,综合表明iGrav-007观测性能良好。
关键词: iGrav-007     大气重力导纳值     重力残差     自由振荡模态0Sm    
iGrav-007 SG and Detection of the Spherical Free Oscillation Modes 0Sm
JIA Jiangang, LUAN Wei, SHEN Wenbin     
School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China
First author: JIA Jiangang, Phd candidate, specializes in determination of earth gravity field. E-mail:jgjia@sgg.whu.edu.cn
Foundation support: The Major State Basic Research Development Program of China (973 Program), No.2013CB733305;the National Natural Science Foundation of China, Nos. 41128003,41174011,41210006,40974015; the Experiment Technology Project of Wuhan University, No.WHU-2015-SYJS-08.
Abstract: iGrav Superconducting Gravimeter is the newest portable gravimeter with high accuracy in the world, which provides the most stable relative gravity measurement. Based on Tsoft software, this paper processed the consecutive iGrav-007 observation data of 5 months at Wuhan Jiufeng Station. Correction of station air pression and variations of polar motion on the gravity field and analysis of the nature of gravity residual error and drifting was made. Through the spectrum analysis, detaction of the spherical free oscillation modes 0Sm caused by the 2013 Lushan Earthquake(Mw=6.6) was made. The result of detaction is evident. It comprehensively shows that iGrav-007 having a good performance on observation.
Key words: iGrav-007     atmospheric gravity admittance     gravity residual     free oscillation modes 0Sm    

地球重力场变化与各种物理过程(如潮汐、构造活动、地球自由振荡等)紧密联系,国内外许多学者都试图通过重力潮汐的观测来探索地球内部的奥秘并解释与之相关的地球物理现象[1, 2, 3]。1997年7月启动的“全球地球动力学计划”(Global Geodynamics Projects,GGP)为研究全球动力学现象,如地球固体潮、重力场与海洋大气耦合、地球自由核章动、地球自转变化等导致的重力效应,提供了全球超导重力观测资料[4, 5, 6]

超导重力仪被认为是精度最高,稳定性和连续性最好的重力仪,对重力变化的灵敏度极高 (在频段低于0.8 mHz时超过当前最好的宽频地震仪[7]),适合长周期信号的观测,在自由振荡等研究中优势明显[8, 9, 10]。美国GWR公司生产的iGrav超导重力仪是最新型的便携式相对重力仪,其基本原理与台站型超导重力仪(OSG)解释相同,体积更小更为便携。其标称精度在频域中为1 nanoGal,在时域中1 min平均值为0.05 μGal,零漂小于0.5 μGal/月[11]

利用超导重力仪获取高精度地球重力信号,需进行各项必要的改正及残差分析。本文基于Tsoft软件[12]处理武汉九峰台站连续5个月的超导重力观测数据,进行气压负荷、极移重力改正及尖峰、阶越、间断等干扰信号的修正;进一步分析重力残差及仪器的长期漂移特征,评估观测精度。采用功率谱密度估计(Spanos-Tratskas方法)[13],检测2013年Mw6.6级芦山地震激发的地球自由振荡球型模态0Sm并分析各频段上的信噪比。

1 观测数据的预处理和分析 1.1 观测资料及预处理方法

武汉大学的iGrav-007超导重力仪引进于2012年,次年3月安装于中科院测地所武汉九峰台站。由于观测显示安装前3个月仪器尚未完全稳定,因此本文取2013年6月25日~11月25日连续5个月的原始秒采样观测数据,并采用最小二乘滤波获得分采样数据。

利用Tsoft软件进行数据预处理的过程包括:(1) 根据台站潮汐参数构制重力合成潮,得到重力残差;(2) 利用Tsoft多种修正方法去除尖峰、阶跃和间断等干扰信号;(3) 由回归分析法计算大气重力导纳值,改正大气负荷效应;(4) 根据IERS提供的地球自转参数,改正极移重力效应。对经大气负荷和极移重力改正后的重力残差分别进行多项式拟合和滤波,去除仪器漂移和残余潮汐信号,得到最终重力残差 ε(t)(单位为nm·s-2):

式中,t为采样历元;yobs(t)为iGrav-007的重力潮汐观测值;ymod(t)为台站重力合成潮;dgP(t)为气压引起的台站重力变化;dgpole(t)为极移引起的台站重力变化; 表示拟合仪器漂移和残余潮汐信号的p次多项式,ai是多项式系数。

1.2 重力合成潮汐

台站重力合成潮ymod(t)可表示为:

式中,n为波群数;N为波群总数;δn和Δφn为第n个波群数的模型振幅因子和相位滞后;m为第n个波群的潮汐分量数,Mn为第n个波群的潮汐分量总数;Hnm、ωnmφnm为第n个波群第m个潮汐分量的模型振幅、角频率和初始相位[2]

本文采用台站提供的重力合成潮波群(共17个波群)参数,并置测站点的大地经、纬度λ=114.489 80°E,φ=30.515 87°N,仪器高程H=89.3 m,构制ymod(t),如图1(b)所示。

图 1 重力残差中的台站大气负荷和极移重力改正 Fig. 1 Correction of Air Pression and Polar Motion at Station on the Gravity Residuals
1.3 大气负荷改正

大气负荷改正的常用方法是将重力潮汐观测值减去重力合成潮获得重力残差yresd(t),如图1(c)所示,再将yresd(t)与台站气压观测值P(t)进行最小二乘回归求得大气重力导纳值C,乘积CP(t)即为大气负荷改正量。

在时域内获得大气重力导纳值,一元线性回归模型为:

气压引起的台站重力变化dgP(t)可表示为:

采用最小二乘线性迭代回归法,以残差标准差σ为中误差,3σ为限差,每次迭代时均去掉原始重力残差和气压观测值中的二次多项式拟合部分,当C的变化小于0.000 1 nm·s-2/mBar时停止迭代,实际迭代4次即可,得到台站大气重力导纳值C=(-2.480 7± 0.003 5) nm·s-2/mBar,斜距C1= 4 361.244 8 nm·s-2。由图1(d)可知,台站大气负荷的影响振幅约100 nm·s-2。经大气负荷改正后,重力残差得到明显改善,如图1(e)所示,可见大气负荷效应是重力残差中的主要成分。

1.4 极移重力效应改正

极移引起的地表任一点重力变化可根据该点坐标和极坐标的变化(从IERS提供的地球自转参数中得到)算得,故极移重力效应可作为已知量扣除[14]。极移引起的台站重力变化dgpole可表示为:

式中,地球平均半径 a=6 378 137.0 m;地球平均自转角速度 ω=7.292 115×10-5 rad/s;(x,y)为测站点的极坐标;大地经度λ、大地纬度φ与上文给出的值相同。

根据式(5)算得武汉九峰台站的极移重力效应,并取与图1(e)中残差相同的时间段。因IERS提供的数据采样率为1 d,采用三次样条插值获得分采样的极移重力,如图1(f)所示,极移影响的振幅约20 nm·s-2,较明显。

1.5 经大气和极移改正后的重力残差分析

对最终重力残差线性拟合得iGrav-007的漂移率为0.40 μGal/月,如图2(a)所示,较其标称值小。采用20阶多项式拟合重力残差,如图2(b)所示,得到拟合残差均方根RMS=0.423 μGal/min (经实验,采用更高阶多项式RMS变化均小于0.01 μGal)。该值大于其标称值0.050 μGal/min。究其原因,一是因台站背景噪声较大;二是大气负荷改正不彻底,未考虑中远和远区的气压影响[15]以及对重力潮汐观测影响的频率依赖,仅采用线性回归计算大气重力导纳值[15, 16, 17]

图 2 重力残差的线性拟合、20次多项式拟合及多项式拟合残差 Fig. 2 Linear Fitting, 20 Times Polynomial Fitting and Polynomial Residual Error on the Gravity Residuals

对最终重力残差分别进行最小二乘高通、带通、低通滤波,截止频率分别以合成潮M3波的最高频率3.081 254 cpd和MM波的最低频率0.020 885 cpd为界。此处认为比M3波频更高的部分由仪器的高频噪声引起,比MM波频更低的部 分为仪器漂移,中间部分主要为残余潮汐。故将重力残差分为短周期、潮汐频段、长周期三部分。由图3(a)可知,重力残差中存在微弱的残余潮汐信号,能量约在±3 nm·s-2之间,相对较小。以短周期部分扣除平均值后的均方根0.040 μGal/min衡量其观测精度,表明iGrav-007观测性能良好。

图 3 重力残差的最小二乘滤波 Fig. 3 Least Squares Filtering on the Gravity Residuals
2 iGrav-007对球型自由振荡模态0Sm的检测 2.1 自由振荡信号检测方法

2013年芦山地震在2.3~5.3 mHz频段激 发了微弱的球型简正模信号,为研究iGrav-007的 高频特性提供了条件。采用功率谱密度估计,检测该地震激发的地球自由振荡球型模态0Sm,分析各频段上的信噪比,并与PREM模型[18]理论预测值比较,评价iGrav-007的残差水平。芦山地震发生于2013年4月20日0:2:48(UTC),选取了地震后18 h的iGrav-007观测数据,如图4所示。

图 4 iGrav-007记录的芦山地震后18 h的原始观测数据/V Fig. 4 Hours Original Recording After Lushan Earthquake by iGrav-007 /V
2.2 观测结果分析

在0.28~1.80 mHz频段上,自由振荡信号受噪声干扰严重,0S0~0S10基型振荡无法有效检测。分析原因,一是这些频率的低阶振型受地震强度较小的限制,难以激发;二是地球自转和椭率可导致自由振荡信号的分裂和耦合[19, 20]

图5(a)图5(b)分别为iGrav-007在1.8~3.9 mHz和3.9~5.5 mHz频段检测到的芦山地震激发的自由振荡振型。数据长度为18 h,傅里叶变换采用Hanning窗。由图5可知iGrav-007清晰地检测到了0S11~0S48之间除0S16振型之外的振型基频系列。为对比观测到的自由振荡周期与PREM理论值,图5中用虚线指示PREM模型的理论频率值。0S16振型未形成突出的谱峰,主要是受附近的噪声影响较大所致。1.80~3.00 mHz频段信噪比较低,但功率谱密度能较好地反映基频球型振荡频率。在3.00~4.75 mHz频段检测到0S21~0S40的所有振型,信噪比均大于3,谱峰清晰,大部分基频振型观测值与PREM理论值偏差较小。其中0S30振型虽形成清晰的谱峰,但观测值与理论值频率偏差相对其他频段较大,达到20.0×10-6 Hz(约8倍的傅氏谱分析频率分辨率)。究其原因,一是由于该振型附近噪声较大,二是这一基频振型可能受到地球上地幔剪切波的干扰[21],且0S31振型处可能出现谱线分裂现象。在4.75~5.50 mHz频段检测到0S41~0S48振型,该频段附近噪声较低,谱峰清晰可见,与PREM理论值较吻合。

图 5 iGrav-007记录到的芦山地震激发的自由振荡信号傅里叶变换频谱 Fig. 5 FFT Transform Spectrum of Earth’s Free Oscillation Excited by Lushan Earthquake Recorded by iGrav-007

表1可知,1.80~5.50 mHz频段iGrav-007检测到的基频自由振荡信号和PREM理论值偏差仅3个基频振型偏差绝对值超过0.30%,偏差绝对值小于0.15%的检测振型数占总振型数的65.7%。

表 1 自由振荡模态0S110S48的基频观测结果 Tab. 1 The Observation of Basic Frequencies Free Oscillation Modes from 0S11 to 0S48
振型PREM理论预测值/mHziGrav-007观测值/mHz偏差/%
0S11 1.862 1.867 0.27
0S12 1.99 1.991 0.05
0S13 2.113 2.114 0.05
0S14 2.231 2.238 0.31
0S15 2.346 2.346 0.00
0S16 2.458 不明显 ****
0S17 2.567 2.567 0.00
0S18 2.673 2.67 -0.11
0S19 2.777 2.778 0.04
0S20 2.878 2.878 0.00
0S21 2.978 2.978 0.00
0S22 3.075 3.071 -0.13
0S23 3.171 3.172 0.03
0S24 3.266 3.272 0.18
0S25 3.359 3.359 0.00
0S26 3.452 3.451 -0.03
0S27 3.544 3.549 0.14
0S28 3.635 3.637 -0.05
0S29 3.725 3.719 -0.16
0S30 3.816 3.796 -0.52
0S31 3.905 3.92 0.38
0S32 3.995 3.997 0.05
0S33 4.084 4.09 0.15
0S34 4.174 4.167 -0.17
0S35 4.263 4.275 0.28
0S36 4.352 4.352 0.00
0S37 4.442 4.444 0.05
0S38 4.531 4.522 -0.20
0S39 4.62 4.614 -0.13
0S40 4.71 4.722 0.25
0S41 4.8 4.799 -0.02
0S42 4.889 4.892 0.06
0S43 4.979 4.985 0.12
0S44 5.069 5.062 -0.14
0S45 5.159 5.170 0.21
0S46 5.249 5.231 0.34
0S47 5.339 5.340 0.02
0S48 5.429 5.417 0.22
3 结 语

本文基于iGrav-007在武汉九峰台站连续5个月的观测数据,应用Tsoft软件进行大气负荷和极移重力效应改正,采用最小二乘线性迭代回归,在时域内得到武汉九峰台站大气重力导纳值为(-2.480 7±0.003 5) nm·s-2/mBar。利用线性拟合最终重力残差,得到iGrav-007的漂移率为0.40 μGal/月;利用多项式拟合重力残差,去除仪器长期漂移量和残余潮汐信号,得到最终重力残差均方根为0.423 μGal/ min;利用滤波,以重力残差短周期部分扣除平均值后的均方根0.040 μGal/min衡量观测精度,表明iGrav-007观测性能良好。

iGrav-007很好地记录了2013年Mw6.6芦山地震,检测出1.80~5.50 mHz频段内0S11~0S40振型中除0S16振型外的基频振型自由振荡系列。观测值与PREM理论值偏差绝对值超过0.30%的有3个基频振型,小于0.15%的检测振型数占总振型数的65.7%。

下一步将对iGrav-007固体潮观测数据作调和分析,综合考虑固体潮、气压、极移、日长等变化的影响,获得更准确的潮汐参数和频域内的大气重力导纳值。

致谢: 感谢中国科学院测量与地球物理研究所提供的武汉九峰台站潮汐参数。

参考文献
[1] Wang Dijin,Shen Wenbin. Development of Investigations in some Geoscientific Problems Using Superconducting Gravimeter Data[J]. Science of Surveying and Mapping,2009,34:10-11(王迪晋,申文斌.利用超导重力数据研究若干地球科学问题的进展[J].测绘科学,2009,34:10-11)
[2] Chen Xiaodong,Sun Heping.New Method for Pre- processing and Analyzing Tidal Gravity Observations [J]. Journal of Geodesy and Geodynamics, 2002,22:83-87 (陈晓东,孙和平.一种新的重力潮汐数据预处理和分析方法[J].大地测量与地球动力学,2002,22:83-87)
[3] Shen Wenbin,Liu Renli.Study of Detecting the Inner Core Super Rotation Using SG Data[J].Geomatics and Information Science of Wuhan University,2009,34(1):72-76(申文斌,刘任莉.利用超导重力数据探测内核超速旋转的研究[J].武汉大学学报·5信息科学版,2009,34(1):72-76)
[4] Crossley D, Hinderer J. Report of GGP Activities to Commission 3, Completing 10 Years for the Worldwide Network of Superconducting Gravimeters[J]. Observing Our Changing Earth,IAG Symposia, 2008,133:511-521
[5] Warburton R J,Pillai H,Reineman R C. Initial Results with the New GWR iGravTM Superconducting Gravity Meter[C]. IAG Symposium Proceedings, Saint Petersbuly, Russia, 2010
[6] Rosat S, Rogister Y, Crossley D, et al. A Search for the Slichter Triplet with Superconducting Gravimeters: Impact of the Density Jump at the Inner Core Boundary[J]. Journal of Geodynamics, 2006,41(1-3): 296-306
[7] Rosat S,Hinderer J, Crossley D, et al. Performance of Superconducting Gravimeters from Long Period Seismology to Tides[J]. Journal of Geodynamics, 2004, 38:461-476
[8] Widmer R. What can Superconducting Gravimeters Contribute to Normal Mode Seismology[J]. Bull.Seism. Soc.Am, 2003, 93(3):1 370-1 380
[9] Hu Xiaogang, Liu Lintao, Sun Heping, et al. Wavelet Filter Analysis of Splitting and Coupling of Seismic Normal Modes Below 1.5 mHz with Superconducting Gravimeter Records After the December 26th 2004 Sumatra Earthquake[J].Science in China Series D: Earth Sciences,2006,49:1 259-1 269
[10] Wang Linsong, Chen Chao, Liang Qing, et al. Coseismic Response to Gravity Signals of the Great Earthquake in East Japan[J]. Geomatics and Information Science of Wuhan University,2012,37(11):1 348- 1 351(王林松,陈超,梁青,等,东日本大地震重力信号同震响应[J].武汉大学学报·5信息科学版,2012,37(11): 1 348-1 351)
[11] GWR Instruments, Inc. iGrav007 User's Guide[OL]. http://www.gwrinstruments.com,2012
[12] Van C M,Vauterin P. Tsoft:Graphical and Interactive Software for the Analysis of Time Series and Earth Tides[J].Computers & Geosciences,2005,31:631- 640
[13] Spanos D,Failla G. Evolutionary Spectra Estimation Using Wavelets[J].Journal of Engineering Mechanics, 2004,130(8):952-960
[14] Banka D, Crossley D L. Noise Levers of Super- conducting Gravimeter at Seismic Frequencies[J]. Geophys.J.Int,1999,139(1):87-97
[15] Sun Heping, Ducarme B, Dehant V. Correction of the Atmosphere Pressure on Gravity Measurements Recorded by a Superconducting Gravimeter at Brussels[C]. The 12th Int. Sympos.On Earth Tides, Beijing, 1993
[16] Sun Heping, Xu Houze, Tao Guoxiang, et al. Accurate Determination of the Earth Tidal Parameters Using Superconducting Gravimeter Observation[J].Crustal Deformation and Earthquake,1997,17(4):17-25(孙和平,许厚泽, 陶国祥,等.用超导重力仪观测数据精密测定地球潮汐常数[J].地壳形变与地震, 1997,17(4):17-25)
[17] Crossley D, Jensen O, Hinderer J. Effective Barometric Admittance and Gravity Residual[J]. Phys. Earth Planet.Int, 1995,90(3-4):1 679-1 682
[18] Dziewonski A, Anderson D. Preliminary Reference Earth Model[J].Phys. Earth Planet.Int,1981,25:297-356
[19] Lei Xiange, Xu Houze, Sun Heping. Detection of Earth Free Oscillation by Superconducting Gravity Observation Data[J]. Chinese Science Bulletin,2002,47 (18):1 432-1 436(雷湘鄂,许厚泽,孙和平.利用超导重力观测资料检测地球自由振荡[J].科学通报, 2002,47(18):1 432- 1 436)
[20] Lei Xiange, Sun Heping, Xu Houze. Test and Discussion of the Earth Free Oscillations and Spectral Line Splitting of Sumatra Earthquake[J].Science China: Earth Sciences,2007,37(4):504-511(雷湘鄂,孙和平,许厚泽.苏门达腊大地震激发的地球自由振荡及其谱线分裂的检测与讨论[J].中国科学:地球科学, 2007,37(4):504-511)
[21] Ishii M, Tromp J, Dziewonski A M, et al. Joint Inversion of Normal Mode and Body Wave Data for Inner Core Anisotropy: Laterally Homogeneous Anisotropy[J].J.Geophys.Res,2002,107(B12): 2 379-2 395