超導(dǎo)重力儀背景噪聲水平分析

江　穎　劉子維　張苗苗　郝洪濤　李　輝

1　中國(guó)地震局地震研究所(地震大地測(cè)量重點(diǎn)實(shí)驗(yàn)室)，武漢市洪山側(cè)路40號(hào)，430071　2　中國(guó)科學(xué)院測(cè)量與地球物理研究所大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，武漢市徐東大街340號(hào)，430077

?

超導(dǎo)重力儀背景噪聲水平分析

江穎1,2劉子維1張苗苗2郝洪濤1李輝1

1中國(guó)地震局地震研究所(地震大地測(cè)量重點(diǎn)實(shí)驗(yàn)室)，武漢市洪山側(cè)路40號(hào)，4300712中國(guó)科學(xué)院測(cè)量與地球物理研究所大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，武漢市徐東大街340號(hào)，430077

摘要：基于拉薩超導(dǎo)重力儀(SG)4 a的連續(xù)觀測(cè)數(shù)據(jù)計(jì)算其在地震頻段(200～600 s)的背景噪聲水平，并研究氣壓、采樣率等對(duì)功率譜密度(PSD)的影響。分析比較全球地球動(dòng)力學(xué)計(jì)劃(GGP)公布的2013年全球SG在地震頻段和亞地震頻段的背景噪聲水平，結(jié)果發(fā)現(xiàn)：1)在頻率低于10-3Hz時(shí)，潮汐和氣壓改正可以很好地改進(jìn)重力PSD；2)在頻率低于0.5×10-3Hz時(shí)，超導(dǎo)重力儀的噪聲水平明顯優(yōu)于地震儀，表明SG更適合研究長(zhǎng)周期地震和亞地震模；3)全球SG地震頻段噪聲水平的范圍為0.180～1.964，亞地震頻段噪聲水平的范圍為1.860～3.853。

關(guān)鍵詞：超導(dǎo)重力儀；噪聲水平；新低噪聲模型；功率譜密度

超導(dǎo)重力儀(SG)[1-4]背景噪聲的研究有助于臺(tái)站選址、儀器調(diào)整及判斷儀器在近期是否具有檢測(cè)地震簡(jiǎn)正模的能力。利用背景噪聲水平結(jié)合連續(xù)觀測(cè)數(shù)據(jù)可以確定全球地球參數(shù)，如自由核章動(dòng)和錢德勒擺動(dòng)、地球內(nèi)核平動(dòng)振蕩引起的微小重力變化[3,5]。利用觀測(cè)序列的功率譜密度(PSD)能夠很好地分析SG的噪聲水平。Peterson[6]基于全球75個(gè)地震儀臺(tái)站功率譜密度的下包絡(luò)線得到新低噪聲模型(NLNM)，作為估計(jì)站點(diǎn)傳感器質(zhì)量的參考水平。Banka[7]認(rèn)為背景噪聲水平可以從PSD中提取出來，將一個(gè)窄窗在200～600 s簡(jiǎn)正模頻段的平均功率譜密度稱為地震噪聲水平(SNM)。利用該方法對(duì)SG、彈簧重力儀、絕對(duì)重力儀和寬頻帶地震儀的背景噪聲水平進(jìn)行比較發(fā)現(xiàn)，最好的超導(dǎo)儀在1 mHz以下的頻段要比其他儀器安靜[8-12]。

本文利用拉薩SG連續(xù)4 a的觀測(cè)記錄，研究氣壓、采樣率等對(duì)功率譜密度的影響，分別計(jì)算了拉薩超導(dǎo)重力儀4 a連續(xù)觀測(cè)數(shù)據(jù)在地震頻段(200～600 s)的背景噪聲水平；基于相同方法，分析比較了地球動(dòng)力學(xué)計(jì)劃(GGP)公布的2013年全球SG分別在地震頻段和亞地震頻段(1～6 h)的背景噪聲水平。

1　數(shù)據(jù)處理

對(duì)SG在此期間每年每天的重力和氣壓數(shù)據(jù)分別進(jìn)行預(yù)處理。首先對(duì)氣壓數(shù)據(jù)進(jìn)行預(yù)處理，去除突跳和階躍等；然后基于彈性參考地球模型計(jì)算得到合成潮，并在觀測(cè)數(shù)據(jù)中去掉合成潮汐信號(hào)[13]；去掉氣壓對(duì)重力信號(hào)的影響，大氣導(dǎo)納值取-0.3 μGal/mbar。

為確定一個(gè)時(shí)間序列的頻率內(nèi)容，需要進(jìn)行譜分析[14]。如果f(j)為信號(hào)，j=1，…，N，采樣率是Δt，那么振幅譜s(k)和功率譜密度PSD為[1]：

(1)

PSD=s2Δt/N

(2)

(3)

其中，k1和k2分別是頻率的范圍，M為k1和k2之間的樣本數(shù)量。對(duì)經(jīng)過以上處理的每天的重力殘差信號(hào)計(jì)算均方根值(RMS)，并選出RMS最小的5 d，用來代表該年背景噪聲最小的5 d。傅氏分析后計(jì)算最平靜5 d的平均功率譜，計(jì)算地震頻段的平均功率譜密度和SNM：

SNM=log10(meanPSD)+0.5

(4)

其中，PSD單位為(nm/s2)2/Hz。SNM作為儀器噪聲水平的一個(gè)品質(zhì)因子，可以很好地對(duì)比儀器在200～600 s頻段內(nèi)的表現(xiàn)。同樣，在利用式(4)計(jì)算亞地震噪聲水平(SSNM)時(shí)，應(yīng)選取相應(yīng)頻段范圍的數(shù)據(jù)進(jìn)行計(jì)算。

2　拉薩超導(dǎo)重力儀地震頻段噪聲水平分析

2.1日重力數(shù)據(jù)的PSD

為研究不同程度的預(yù)處理對(duì)SG數(shù)據(jù)PSD結(jié)果的影響，取拉薩臺(tái)2010-01-14重力記錄進(jìn)行分析，其中采樣率為1 s。

從圖1可以看出，原始重力信號(hào)的振幅約為1 600 nm/s2，經(jīng)過潮汐和氣壓改正后，振幅約為40 nm/s2；扣除9階多項(xiàng)式后，振幅約為8 nm/s2。潮汐和大氣改正很大程度上可以改進(jìn)重力數(shù)據(jù)，扣除9階多項(xiàng)式，得到更精準(zhǔn)的重力殘差?；谝陨?種不同的數(shù)據(jù)預(yù)處理結(jié)果，分別計(jì)算相應(yīng)的重力PSD，結(jié)果如圖2。其中，2(a)為原始重力觀測(cè)數(shù)據(jù)的PSD；2(b)為去掉潮汐和大氣影響的重力PSD；2(c)為扣掉9階多項(xiàng)式的重力PSD；NLNM為新低噪聲水平。

圖1　拉薩臺(tái)站2010-01-14連續(xù)重力觀測(cè)的預(yù)處理結(jié)果Fig.1　The gravity observation of the Lhasa SG on Jan.14 2010

由圖2可見，每經(jīng)過一個(gè)處理步驟，噪聲水平都有一定程度的降低。潮汐和大氣改正可以很好地改進(jìn)重力數(shù)據(jù)的PSD，尤其在頻率低于10-3Hz時(shí)，對(duì)地震頻段內(nèi)重力PSD的影響比對(duì)亞地震頻段的影響小。因此，利用重力數(shù)據(jù)檢測(cè)低頻信號(hào)時(shí)，如0S2、0S3、2S1、1S1，必須進(jìn)行嚴(yán)密的潮汐、大氣以及9階多項(xiàng)式改正，以降低噪聲水平。

圖2　拉薩臺(tái)站2010-01-14經(jīng)過不同程度改正的重力PSDFig.2　The gravity PSD of observation after correction of the Lhasa SG on Jan.14 2010

2.2采樣率對(duì)PSD的影響

為了研究采樣率對(duì)PSD的影響，將拉薩超導(dǎo)重力儀2010年全年采樣率為1 s和1 min的連續(xù)重力數(shù)據(jù)分別進(jìn)行氣壓和潮汐改正，并扣掉9階多項(xiàng)式，計(jì)算每天重力殘差的均方根值(RMS)，然后選出RMS最小的5 d。對(duì)2010年連續(xù)觀測(cè)的重力數(shù)據(jù)進(jìn)行分析，雖然由于采樣率不同，最平靜的5 d并不重合，但是其RMS結(jié)果非常相近。分別計(jì)算秒采樣和分鐘采樣的重力PSD最平靜5 d 的平均振幅譜及其功率譜，如圖3(其中TNM為潮汐頻段范圍)。結(jié)果發(fā)現(xiàn)，在秒采樣的重力PSD中，在頻率為8×10-3Hz處有一個(gè)非常明顯的超導(dǎo)共振模，這是由于一個(gè)球型水平旋轉(zhuǎn)振動(dòng)耦合進(jìn)了垂直力和反饋系統(tǒng)。另外可以明顯看出，分鐘采樣的重力PSD在高于1.67 mHz處急劇降低，對(duì)低頻的影響較小，即不同的采樣率對(duì)地震頻段的重力PSD影響不大。計(jì)算地震頻段的平均PSD和SNM發(fā)現(xiàn)，不同采樣率的重力PSD在地震頻段的SNM差別比較小，秒采樣數(shù)據(jù)得到的SNM為1.231，分鐘采樣數(shù)據(jù)得到的SNM為0.904。為了便于計(jì)算及與其他SG作比較[9-10]，下文均采用分鐘采樣數(shù)據(jù)。

圖3　拉薩超導(dǎo)臺(tái)站2010年最平靜5 d的重力PSDFig.3　The gravity PSD of the 5 quietest days in 2010 of the Lhasa SG

2.3噪聲水平計(jì)算

采用相同的方法計(jì)算拉薩超導(dǎo)重力儀2010～2013年共4 a連續(xù)重力的地震頻段噪聲水平，結(jié)果列于表1，相應(yīng)的重力PSD如圖4。可以明顯看出，頻率低于0.5×10-3時(shí)，SG的噪聲水平優(yōu)于地震儀，說明超導(dǎo)重力儀更適合研究長(zhǎng)周期地震和亞地震模。4 a的SNM范圍為0.904～1.044，4 a中最平靜5 d得到的SNM為0.903，說明拉薩臺(tái)站屬于低噪聲臺(tái)站。

表1　拉薩超導(dǎo)重力儀2010～2013年的背景噪聲水平

圖4　拉薩超導(dǎo)重力儀2010～2013年重力PSDFig.4　The gravity PSD of 2010～2013 of the Lhasa SG

3　全球SG噪聲水平分析

圖5　全球SG地震頻段噪聲水平Fig.5　The SNM of global SGs

圖6　全球SG亞地震頻段噪聲水平Fig.6　The SSNM of global SGs

根據(jù)GGP網(wǎng)站更新的2013年全球SG重力數(shù)據(jù)，計(jì)算15個(gè)SG分別在地震頻段和亞地震頻段的背景噪聲水平，包括中國(guó)的拉薩(LH)、武漢(WH)和武漢大學(xué)的iGrave，美國(guó)的Apache Point(AP),德國(guó)的Schiltach(BF)、Bad Homburg(BH)、Moxa(MO)和Wettzell(WE)，加拿大的Cantley(CA)，意大利的Medicina(MC)，芬蘭的Metsahovi(ME)，捷克的Pecny(PE)，南非的Sutherland(SU)，智利的Concepcion(TC)和西班牙的Yebes(YS)。2009年以后，新布設(shè)的超導(dǎo)重力儀包括AP、BF、LH和YS。其中，WH和iGrave重力儀在地震頻段和亞地震頻段的地震噪聲水平已經(jīng)由Zhang等[15]給出。本文計(jì)算以上超導(dǎo)重力儀在地震頻段和亞地震頻段的噪聲水平，結(jié)果如圖5和圖6。由圖可知，SNM的范圍為0.180(BF)～1.964(TC)，其中拉薩SG、BH超導(dǎo)重力儀和iGrave重力儀的地震頻段噪聲水平基本一致，優(yōu)于全球其他10臺(tái)SG。德國(guó)新增加的BF超導(dǎo)重力儀在地震頻段噪聲水平最低。SSNM的范圍為1.860(WH)～3.853(TC)，其中中國(guó)的3臺(tái)超導(dǎo)重力儀在亞地震頻段的噪聲水平最低。相比于Rosat等的結(jié)果，武漢SG搬離了市區(qū)，很大程度上降低了其噪聲水平。拉薩SG的SSNM優(yōu)于全球其他12個(gè)SG。全球SG的背景噪聲水平分析為全球動(dòng)力學(xué)信號(hào)的探測(cè)提供了數(shù)據(jù)選擇依據(jù)。選用噪聲水平較低的SG觀測(cè)數(shù)據(jù)，結(jié)合有效的疊積方法，可以進(jìn)一步壓制噪聲水平，有利于提取微弱信號(hào)。

4　結(jié)　語

本文利用拉薩超導(dǎo)重力儀連續(xù)4 a的觀測(cè)記錄研究氣壓、采樣率等對(duì)功率譜密度的影響，分別計(jì)算其在地震頻段和亞地震頻段的背景噪聲水平。結(jié)果表明，潮汐和大氣改正在頻率低于10-3Hz時(shí)可以很好地改進(jìn)重力PSD，對(duì)地震噪聲頻段內(nèi)重力PSD的影響要比對(duì)亞地震頻段內(nèi)的小。基于秒采樣重力觀測(cè)計(jì)算得到的SNM大于基于分鐘采樣的結(jié)果，分別為1.231和0.904。在頻率低于0.5×10-3Hz時(shí)，超導(dǎo)重力儀的噪聲水平明顯優(yōu)于地震儀，說明超導(dǎo)儀更適合研究長(zhǎng)周期地震和亞地震模。拉薩超導(dǎo)重力儀4 a的SNM范圍為0.904～1.044，4 a中最平靜5 d得到的SNM為0.903，拉薩臺(tái)站屬于低噪聲臺(tái)站。分析比較了GGP公布的2013年全球超導(dǎo)重力儀分別在地震頻段和亞地震頻段的背景噪聲水平，SNM的范圍為0.180～1.964，SSNM的范圍為1.860 ～3.853。武漢SG搬離了市區(qū)，很大程度上降低了其噪聲水平；拉薩SG在潮汐頻段的噪聲水平優(yōu)于全球其他SG。

基于以上結(jié)果，利用重力數(shù)據(jù)檢測(cè)低頻信號(hào)時(shí)，必須進(jìn)行嚴(yán)密的潮汐和大氣改正，以降低噪聲水平。相比較而言，超導(dǎo)重力儀具有低頻性、穩(wěn)定性，在探測(cè)地球自由振蕩信號(hào)方面，尤其是低頻信號(hào)方面具有非常大的優(yōu)勢(shì)。全球SG的背景噪聲水平分析為全球動(dòng)力學(xué)信號(hào)的探測(cè)提供了數(shù)據(jù)選擇依據(jù)，選用噪聲水平較低的SG觀測(cè)數(shù)據(jù)，結(jié)合有效的疊積方法，有利于提高信號(hào)的信噪比。

致謝：感謝中國(guó)科學(xué)院測(cè)量與地球物理研究所提供拉薩超導(dǎo)重力儀觀測(cè)數(shù)據(jù)，全球地球動(dòng)力學(xué)計(jì)劃提供超導(dǎo)重力儀觀測(cè)數(shù)據(jù)。

參考文獻(xiàn)

[1]Banka D, Crossley D. Noise Levels of Superconducting Gravimeters at Seismic Frequencies [J]. Geophys J Int, 1999,139(1): 87-97

[2]孫和平, 許厚澤.國(guó)際地球動(dòng)力學(xué)合作項(xiàng)目的實(shí)施與展望[J]. 地球科學(xué)進(jìn)展, 1997,12(2):152-157(Sun Heping, Xu Houze. Execution and Prospect for the Global Geodynamics Project Cooperation[J]. Advance in Earth Science, 1997,12(2):152-157)

[3]徐建橋, 陳曉東, 周江存, 等.拉薩重力潮汐變化特征[J].科學(xué)通報(bào), 2012,57(22):2 094-2 101(Xu Jianqiao, Chen Xiaodong, Zhou Jiangcun, et al. Characteristics of Tidal Gravity Changes in Lhasa, Tibet,China[J]. Chinese Science Bulletin, 2012, 57(20): 2 586-2 594)

[4]Hinderer J, Amalvict M, Crossley D,et al. Tides, Earthquakes and Ground Noise as Seen by the Absolute Gravimeter FG5 and Its Superspring; Comparison with a Superconducting Gravimeter and a Broadband Seismometer[J]. Metrologia, 2002,39(5): 495

[5]Jiang Y, Hu X G, Liu C L, et al. Constraining Focal Mechanism of Lushan Earthquake by Observations of the Earth’s Free Oscillation[J]. Science China:Earth Sciences, 2014,57(9):2 064-2 070

[6]Peterson J. Observations and Modeling of Seismic Background Noise[R]. US Department of Interior, Geological Survey, Albuquerque, New Mexico,1993

[7]Banka D. Noise Levels of Superconducting Gravimeters at Seismic Frequencies [D]. GDMB-Information-GesellschaftmbH, Clausthal, Germany, 1997

[8]Rosat S, Hinderer J. Noise Levels of Superconducting Gravimeters: Updated Comparison and Time Stability[J]. Bulletin of the Seismological Society of America, 2011,101(3): 1 233-1 241

[9]Rosat S, Hinderer J, Crossley D. A Comparison of the Seismic Noise Levels at Various GGP Stations[J]. Bull Inf MaréesTerrestres, 2002,135: 10 689-10 700

[10]Rosat S, Hinderer J, Crossley D,et al. The Search for the Slichter Mode: Comparison of Noise Levels of Superconducting Gravimeters and Investigation of a Stacking Method[J]. Physics of the Earth and Planetary Interiors, 2003,140(1): 183-202

[11]Freybourger M, Hinderer J, Trampert J. Comparative Study of Superconducting Gravimeters and Broadband Seismometers STS-1/Z in Seismic and Subseismic Frequency Bands[J]. Physics of the Earth and Planetary Interiors, 1997,101(3): 203-217

[12]Camp M. Measuring Seismic Normal Modes with the GWR C021 Superconducting Gravimeter[J]. Physics of the Earth and Planetary Interiors, 1999,116(1): 81-92

[13]Xi Q W. The Precision of the Development of the Tidal Generating Potential and Some Explanatory Notes[J].Maréesterrestres,1989,105: 7 396-7 404

[14]Bronstein I N, Semendjajew K A. Taschenbuch der Mathematik[M]. Verlag Harri Deutsch: Thun und Frankfurt am Main,1985

[15]Zhang M M, Xu J Q, Sun H P, et al. Comparison of Noise Levels of the New iGrav-007 Superconducting Gravimeter and the SG-065 Superconducting Gravimeter in Wuhan (China)[J]. Bull Inf Mar Terr, 2014, 148: 11 987-12 000

Foundation support:Special Fund for Earthquake Research of CEA, No. 201308009; Director Fund of Institute of Seismology, CEA, No. IS201526225; National Natural Science Foundation of China, No.41404064; Open Fund of State Key Laboratory of Geodesy and Earth’s Dynamics, No. SKLGED2015-1-2-EZ.

About the first author:JIANG Ying, PhD, majors in gravity data analysis and its application. E-mail:jiangyingchen@126.com.

收稿日期：2015-08-31

第一作者簡(jiǎn)介：江穎，博士，助理研究員，主要從事重力數(shù)據(jù)分析及應(yīng)用研究，E-mail: jiangyingchen@126.com。 通訊作者：劉子維，副研究員，主要從事重力觀測(cè)技術(shù)和數(shù)字信號(hào)處理研究，E-mail:Liuzw.99@gmail.com。

DOI：10.14075/j.jgg.2016.08.008

文章編號(hào)：1671-5942(2016)08-0689-05

中圖分類號(hào)：P312

文獻(xiàn)標(biāo)識(shí)碼：A

Corresponding author:LIU Ziwei,associate researcher,majors in gravity observation technique and digital signal processing, E-mail: Liuzw.99@gmail.com.

Research on Background Noise Level of Global Superconducting Gravimeter

JIANGYing1,2LIUZiwei1ZHANGMiaomiao2HAOHongtao1LIHui1

1Key Laboratory of Earthquake Geodesy, Institute of Seismology, CEA, 40 Hongshance Road,Wuhan 430071, China2State Key Laboratory of Geodesy and Earth’s Dynamics, Institute of Geodesy and Geophysics, CAS,340 Xudong Street,Wuhan 430077, China

Abstract:The background noise level in seismic frequency band (200-600 s) based on four years continuous observation of superconducting gravimeter (SG) at Lhasa station is investigated, and the influence of power spectral density (PSD) caused by atmosphere and sampling rate is also studied. The noise level of SG observation provided by global geodynamic project (GGP) for 2013 in seismic frequency band and sub-seismic frequency band are calculated, and they are further analyzed and compared. It is found that the PSD of gravity can be improved by tides and atmosphere correction when the frequency is below 10-3Hz. When the frequency is below 0.5×10-3Hz, the noise level of SG is obviously better than that of a seismometer. This means that SG is more suitable for researching long-period seismic and sub-seismic modes.The range of the seismic noise magnitude of global SG is 0.180-1.964, and the range of sub-seismic noise magnitude is 1.860-3.853.

Key words:superconducting gravimeter; noise level; NLNM; power spectral density

項(xiàng)目來源：中國(guó)地震局地震行業(yè)科研專項(xiàng)(201308009);中國(guó)地震局地震研究所所長(zhǎng)基金(IS201526225)；國(guó)家自然科學(xué)基金(41404064)；大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室開放基金(SKLGED2015-1-2-EZ)。