周江林　沈　萍

1　中國(guó)地震局地球物理研究所,北京市民族大學(xué)南路5號(hào)，100081 2　北京白家疃地球科學(xué)國(guó)家野外科學(xué)觀測(cè)研究站，北京市白家疃，100095

?

北京地震臺(tái)gPhone重力儀對(duì)尼泊爾M8.1強(qiáng)震所激發(fā)地球自由振蕩的探測(cè)

周江林1，2沈萍1

1中國(guó)地震局地球物理研究所,北京市民族大學(xué)南路5號(hào)，100081 2北京白家疃地球科學(xué)國(guó)家野外科學(xué)觀測(cè)研究站，北京市白家疃，100095

利用北京地震臺(tái)gPhone重力儀觀測(cè)資料，采用功率譜密度估計(jì)方法，檢測(cè)2015-04-25尼泊爾M8.1地震激發(fā)的0S0～0S60之間的地球球型自由振蕩，并與地球初步參考模型(PREM)的理論自由振蕩周期進(jìn)行對(duì)比。結(jié)果表明，0S2～0S56基頻自由振蕩的實(shí)測(cè)周期值與PREM 模型的理論值基本一致，證實(shí)了gPhone重力儀能有效檢測(cè)出地球自由振蕩信號(hào)。

地球自由振蕩；gPhone重力儀；功率譜密度；PREM模型；北京地震臺(tái)

地球自由振蕩現(xiàn)象[1-7]的研究與檢測(cè)主要包括地球自由振蕩頻率和相應(yīng)譜峰品質(zhì)因子兩個(gè)方面[8]。這些工作通常都是基于地震技術(shù)開展，用重力技術(shù)檢測(cè)地球自由振蕩的研究較少。gPhone重力儀是目前監(jiān)測(cè)地表垂直形變和重力場(chǎng)變化最精密的儀器之一，具有很寬的動(dòng)態(tài)線性測(cè)量范圍和極低的噪聲水平，能夠檢測(cè)從幾s到幾a周期尺度的地球重力場(chǎng)信號(hào)，具備檢測(cè)球型地球自由振蕩的觀測(cè)能力。

北京地震臺(tái)觀測(cè)山洞位于京郊燕山地震帶上，觀測(cè)室上覆蓋層厚度約為25 m，進(jìn)深約為110 m，采用雙層被覆，至儀器安裝處有5道密封門，洞溫約為12.7 ℃，日氣溫變化<0.05 ℃,年氣溫變化<0.5 ℃，濕度約為70%。2007-05洞室內(nèi)安裝gPhone重力儀，觀測(cè)資料穩(wěn)定、可靠、連續(xù)。2015-04-25尼泊爾發(fā)生8.1級(jí)大地震，北京地震臺(tái)gPhone重力儀清晰地記錄到了這次地震，為提取地球球型自由振蕩信息提供了良好的數(shù)據(jù)基礎(chǔ)。

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

采用改進(jìn)的平均周期圖法求取隨機(jī)信號(hào)的功率譜密度估計(jì)值，并使用快速傅里葉變換算法(FFT)提高運(yùn)算效率。采用功率譜密度估計(jì)的方法分析地球球型自由振蕩，其定義為記錄到的波形資料的自相關(guān)函數(shù)Rn的Fourier變換[9]：

(1)

式中，K=0，1，…，N-1。自相關(guān)函數(shù)Rn可以表示為：

(2)

式中，n=-(N-2)，-(N-1)，…，-1，0，1，…，N-2，N-1，N為所用數(shù)據(jù)的數(shù)目，x為重力數(shù)據(jù)，SK離散值為功率譜密度值。

為消除數(shù)據(jù)不能無(wú)限長(zhǎng)而必須加窗對(duì)功率譜密度估計(jì)的影響，采用Hanning窗來(lái)抑制旁瓣，突出主瓣[9]。

因?yàn)榇蟮卣鸺ぐl(fā)的地球自由振蕩通常能夠持續(xù)數(shù)d的時(shí)間，所以選取數(shù)據(jù)長(zhǎng)度為2015-04-25～30共144 h。本文計(jì)算重力資料采用分鐘值進(jìn)行分析，故N=8 640。在自由振蕩發(fā)生的幾d內(nèi)，不考慮地球重力場(chǎng)長(zhǎng)期變化的影響。

采用功率譜密度估計(jì)值的計(jì)算方法，得到北京地震臺(tái)gPhone重力儀對(duì)尼泊爾地震所激發(fā)球型振蕩觀測(cè)值，并與PREM理論值進(jìn)行對(duì)比。

2　地球自由振蕩的檢測(cè)

北京地震臺(tái)gPhone重力儀觀測(cè)數(shù)據(jù)對(duì)尼泊爾8.1級(jí)地震后的功率譜密度估計(jì)譜見圖1(a)。為將檢測(cè)到的地球自由振蕩信號(hào)清楚地分辨并標(biāo)識(shí)出來(lái)，依據(jù)地球自由振蕩信號(hào)在頻域的分布特點(diǎn)，把功率譜密度估計(jì)結(jié)果分成0.20～1.80 mHz、1.75～2.95 mHz、2.90～3.90 mHz、3.85～4.85 mHz和4.75～6.55 mHz共5個(gè)頻段進(jìn)行分析(圖1(b)～(f))。為避免因分段而可能造成的檢測(cè)信號(hào)丟失，使相鄰頻段保持0.05～0.10 mHz頻域的重疊。

圖1　北京地震臺(tái)gPhone重力儀檢測(cè)到的地球自由振蕩信號(hào)Fig.1　Distribution of the signals in frequency domain of earth’s free oscillations detected by the gPhone gravimeter at Beijing seismic station

計(jì)算結(jié)果清晰地檢測(cè)到了0S0～0S60之間絕大部分的基頻振型自由振蕩系列。參照PREM模型[10]理論計(jì)算值，在各個(gè)譜峰對(duì)應(yīng)的頻率上標(biāo)識(shí)出相應(yīng)的基頻振型(圖1(b)～(f)虛線頂部)，以方便對(duì)比觀測(cè)到的自由振蕩和PREM模型的理論自由振蕩周期。

圖1(b)為0.20～1.80 mHz頻段內(nèi)檢測(cè)到的0S2～0S10基型振蕩。整個(gè)頻段無(wú)較大干擾，各個(gè)振型附近噪音很小，信噪比高。0S0振型觀測(cè)效果不明顯，不能清晰檢測(cè)。雖然0S2與0S3振型相對(duì)于其他振型觀測(cè)效果不明顯，但由于背景噪聲很微弱，仍能識(shí)別出來(lái)。由萬(wàn)永革等[11]的工作可知，0S2與0S3振型存在譜線分裂的情況，且地震一般難以激發(fā)低頻率的低階振型，所以檢測(cè)的信號(hào)較弱。

圖1(c)為1.75～2.95 mHz頻段內(nèi)檢測(cè)到的0S11～0S20基型振蕩。除0S12振型觀測(cè)不明顯外，該頻段整體噪聲較低，信噪比較高，觀測(cè)的振型附近無(wú)較大干擾，觀測(cè)效果清晰。除0S11振型觀測(cè)頻率值與理論值稍有偏差外，其他振型觀測(cè)頻率值與理論值吻合非常好。

圖1(d)為2.90～3.90 mHz頻段內(nèi)檢測(cè)到的0S21～0S30振型。除0S24振型觀測(cè)不明顯外，其余振型檢測(cè)效果良好。0S23、0S25觀測(cè)頻率值與理論值稍有偏差，其他振型吻合非常好。

圖1(e)為3.85～4.85 mHz頻段內(nèi)檢測(cè)到的0S31～0S41振型。該頻段的噪聲較大，0S32、0S33、0S38、0S39、0S40振型未形成突出的譜峰，觀測(cè)效果不明顯。但0S31、0S34、0S35、0S36、0S37、0S41譜線相對(duì)清晰，檢測(cè)效果很好，觀測(cè)頻率值與理倫值偏差較小。

圖1(f)為4.75～6.55 mHz頻段內(nèi)檢測(cè)到的0S41～0S60振型。0S48、0S49振型觀測(cè)效果不明顯；0S50、0S51、0S52振型譜峰不突出，但仍能分辨出來(lái)；0S41～0S47、0S53～0S56振型觀測(cè)效果良好。

由圖1(f)看出，從0S57振型開始，gPhone重力儀已不能檢測(cè)到地球自由振蕩信號(hào)，可能是因?yàn)楦哳l自由振蕩信號(hào)的衰減較快。此外，由于地球是一個(gè)有限的整體，從特定點(diǎn)源發(fā)射的瞬態(tài)波最終還是會(huì)產(chǎn)生干涉，當(dāng)傳播長(zhǎng)度是波長(zhǎng)的整數(shù)倍時(shí)，相互干擾的振型也會(huì)出現(xiàn)。因此，在時(shí)長(zhǎng)幾d(72 h)的波譜尾部(大于4×10-3Hz)出現(xiàn)的較為離散的譜峰與地球自由振蕩周期有關(guān)，譜峰的分離是頻率相差不大的振型干涉所致，所以檢測(cè)級(jí)數(shù)較大的諧頻振型譜峰較為困難[12]。

根據(jù)功率譜密度估計(jì)結(jié)果，統(tǒng)計(jì)北京地震臺(tái)gPhone重力儀檢測(cè)到的尼泊爾M8.1強(qiáng)震激發(fā)的球型基頻振蕩0S0～0S60頻率值，并將結(jié)果與PREM模型[10]理論值進(jìn)行對(duì)比(表1)。由表1可以看出，0S4、0S5、0S9、0S10、0S20、0S21、0S28、0S42共8個(gè)振型的觀測(cè)頻率與理論值偏差近似為0。0S2、0S3因?yàn)槭堑皖l的低階振型，又由于存在譜線分裂的情況，觀測(cè)頻率值與理論值稍有偏差，分別為0.65%與-0.85%。0S0、0S12、0S24、0S32、0S33、0S38、0S39、0S40、0S48、0S49等10個(gè)振型不能與周圍噪聲區(qū)分出來(lái)，未形成突出的譜峰，檢測(cè)效果不明顯，從0S57振型開始的較高頻率的自由振蕩信號(hào)已不能檢測(cè)出來(lái)。絕大部分振型的觀測(cè)周期與理論偏差都在0.20%以內(nèi)，總體平均偏差更是在0.10%以內(nèi)。這說(shuō)明本文對(duì)0S0～0S60基頻自由振蕩的檢測(cè)結(jié)果與PREM模型的理論值基本一致，所以北京地震臺(tái)gPhone重力儀對(duì)地球自由振蕩的檢測(cè)是比較成功的。

表1　北京地震臺(tái)gPhone重力儀球型振蕩觀測(cè)值、PREM模型理論值及其偏差

3　結(jié)　語(yǔ)

運(yùn)用北京地震臺(tái)gPhone重力儀觀測(cè)資料，成功對(duì)尼泊爾M8.1地震激發(fā)的球型自由振蕩進(jìn)行了提取，除0S0、0S12、0S24、0S32、0S33、0S38、0S39、0S40、0S48、0S49等10個(gè)振型觀測(cè)效果不明顯外，準(zhǔn)確檢測(cè)到了0S2～0S56基頻振型自由振蕩。與PREM的理論自由振蕩頻率進(jìn)行對(duì)比，結(jié)果表明，實(shí)測(cè)的頻率值與PREM理論值互相吻合，表明了本文方法的有效性，同時(shí)也反映了北京地震臺(tái)gPhone重力儀優(yōu)秀的觀測(cè)質(zhì)量。

[1]張少泉．地球物理學(xué)概論[M]．北京：地震出版社，2007(Zhang Shaoquan.An Introduction to Geophysics[M].Beijing: Seismological Press,2007)

[2]Lapwood E R，Usami T．Free Oscillations of the Earth[M]．Cambridge:Cambridge University Press，1981

[3]唐磊，邱澤華，闞寶樣．中國(guó)鉆孔體應(yīng)變臺(tái)網(wǎng)觀測(cè)到的地球球型振蕩[J]．大地測(cè)量與地球動(dòng)力學(xué)，2007，27(6)：37-44(Tang Lei,Qiu Zehua,Kan Baoxiang.Earth’s Spheroidal Oscillations Observed with China Borehole Dilatometer Network[J]．Journal of Geodesy and Geodynamics，2007，27(6)：37-44)

[4]嚴(yán)珍珍，張懷，楊長(zhǎng)春，等．地震激發(fā)的地球自由振蕩研究綜述[J]．地球物理學(xué)進(jìn)展，2008，23(3)：686-693(Yan Zhenzhen,Zhang Huai,Yang Changchun,et al. A Review of the Earth’s Free Oscillations Excited by Earthquakes[J].Progress in Geophysics, 2008，23(3)：686-693)

[5]于海英，朱元清，壽海濤．印度洋9．0級(jí)大地震激發(fā)的地球球型振蕩和環(huán)型振蕩[J]．地球物理學(xué)進(jìn)展，2006，21(3)：706-716(Yu Haiying, Zhu Yuanqing, Shou Haitao.Toroidal and Spheroidal Mode of the Earth’s Free Oscilation Excited by the Indian Ocean Earthquake[J]．Progress in Geophysics，2006，21(3)：706-716)

[6]徐曉楓，萬(wàn)永革，王惠琳．由中國(guó)CDSN臺(tái)網(wǎng)檢測(cè)到的汶川地震所激發(fā)的地球球型自由振蕩[J]．地震，2010，30(1)：36-49(Xu Xiaofeng,Wan Yongge,Wang Huilin.Spheroidal Mode of the Earth Free Oscillations Excitated by Wenchuan

Earthquake Recorded by China Digital Seismograph Network[J]．Earthquake，2010，30(1)：36-49)

[7]邱澤華，馬瑾，池順良，等．鉆孔差應(yīng)變儀觀測(cè)的蘇門答臘大地震激發(fā)的地球環(huán)型自由振蕩[J]．地球物理學(xué)報(bào)，2007，50(3)：797-805(Qiu Zehua,Ma Jin,Chi Shunliang,et al.Earth’s Free Torsional Oscillations of the Great Sumatra Earthquake Observed with Borehole Shear Straiumeter[J]．Chinese Journal of Geophysics，2007，50(3)：797-805)

[8]孫和平， 雷湘鄂.重力技術(shù)在地球自由振蕩檢測(cè)研究中的應(yīng)用[J]．地球物理學(xué)進(jìn)展,2005,20(4):974-979(Sun Heping,Lei Xiang’e.Application of the Gravity Technique on Detection of the Earth’s Free Oscillations[J]．Progress in Geophysics,2005,20(4):974-979)

[9]萬(wàn)永革．?dāng)?shù)字信號(hào)處理的MATLAB實(shí)現(xiàn)[M].北京：科學(xué)出版社，2007(Wan Yongge. Digital Signal Processing of MATLAB [M]. Beijing: Science Press, 2007)

[10] Dziewonski A M，Anderson D L．Preliminary Reference Earth Model(PREM)[J]．Phys Earth Planet Inter,198l，25(4)：297-356

[11] 萬(wàn)永革，盛書中，周公威,等．中國(guó)數(shù)字地震臺(tái)網(wǎng)記錄的蘇門答臘-安達(dá)曼地震激發(fā)的地球球型自由振蕩的檢測(cè)[J]．地震學(xué)報(bào)，2007，29(4)：369-380(Wan Yongge,Sheng Shuzhong,Zhou Gongwei,et al.Spheroidal Oscillations of the Earth Stimulated by the Sumatra-Andaman Earthquake with CDSN Data[J]．Acta Seismologica Sinica，2007，29(4)：369-380)

[12] 王林松，陳超，梁青，等.東日本大地震重力信號(hào)同震響應(yīng)[J]．武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2012,37(11):1 348-1 351(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)

Foundation support:Special Fund for Basic Scientific Research of Central Public Research Institutes, No.DQJB14B08.

About the first author:ZHOU Jianglin,senior engineer,majors in earthquake monitoring and prediction,E-mail:zhou.jiang.lin@163.com.

Detecting the Earth Free Oscillations Excited by the Nepal M8.1 Earthquake by Using gPhone Gravimeter at Beijing Seismic Station

ZHOUJianglin1, 2SHENPing1

1Institute of Geophysics, CEA, 5 South-Minzudaxue Road,Beijing 100081,China 2Beijing Baijiatuan Earth Science State Geophysical Observatory,Baijiatuan, Beijing 100095, China

Using the gPhone gravimeter digital observation data recorded at Beijing seismic station, we detect spherical modes of0S0-0S60of the earth free oscillation stimulated by the M8.1 earthquake that occurred on April 25,2015 in Nepal. In our process, the method of power spectral density(PSD) is used．In comparison with the periods of free oscillations of the PREM model, the observed periods coincide well with the theoretical values.The results verify that gPhone gravimeter can detect the earth’s free oscillation.

the earth’s free oscillation;gPhone gravimeter;power spectral density;preliminary reference earth model (PREM); Beijing seismic station

2015-09-08

周江林，高級(jí)工程師，主要從事地震監(jiān)測(cè)預(yù)報(bào)研究，E-mail:zhou.jiang.lin@163.com。

10.14075/j.jgg.2016.09.017

1671-5942(2016)09-0829-04

P315

A

項(xiàng)目來(lái)源：中央級(jí)公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)(DQJB14B08)。