楊躍文 任 佳 楊 星 楊玲英 和玉權(quán) 李麗紅

(1)云南省地震局麗江地震臺(tái),麗江 674100

2)河北省地震局張家口中心臺(tái),張家口 075000 3)云南省地震局,昆明 650224)

利用伸縮儀觀測(cè)資料檢測(cè)地球自由振蕩＊

楊躍文1)任 佳2)楊 星3)楊玲英3)和玉權(quán)1)李麗紅1)

(1)云南省地震局麗江地震臺(tái),麗江 674100

2)河北省地震局張家口中心臺(tái),張家口 075000 3)云南省地震局,昆明 650224)

運(yùn)用云南省云龍地震臺(tái)和山東省泰安地震臺(tái)伸縮儀的數(shù)字化觀測(cè)資料,利用功率譜密度估計(jì)方法,在沒(méi)有對(duì)資料進(jìn)行去固體潮處理的情況下,準(zhǔn)確地獲得了 2004年 12月 26日蘇門答臘 9.0級(jí)地震激發(fā)的0S3～0S26基頻球形自由振蕩,并與地球初步參考模型 PREM的理論自由振蕩周期進(jìn)行了對(duì)比,發(fā)現(xiàn)實(shí)測(cè)振蕩周期與 PRE M預(yù)測(cè)的振蕩周期基本相一致。

伸縮儀;蘇門答臘 9.0級(jí)地震;PREM模型;自由振蕩;球形振蕩

1 引言

地震振動(dòng)產(chǎn)生的地震波,其中體波直接通過(guò)地球傳播,而其他波則沿地球表面?zhèn)鞑ァＣ娌◤恼鹪聪蚋鞣较蛏涑?由于地球是球體,這些不同方向的波最終產(chǎn)生相互干涉,這種干涉方式叫“自由振蕩”。

地球自由振蕩的理論研究至少可追溯到 18世紀(jì)。1882年,Lamb[1]詳細(xì)討論了均勻球體的振型,得到通常稱為球型振蕩和環(huán)型振蕩的兩種振型。最早對(duì)地球自由振蕩的準(zhǔn)確觀測(cè)是在 1960年智利Ms8.3地震后,Benioff[2]和 Ness等[3]分別用 Isabella應(yīng)變儀和 LaCoste-Romberg重力儀實(shí)現(xiàn)的。兩套觀測(cè)結(jié)果非常吻合并與理論值相當(dāng)一致,確認(rèn)了長(zhǎng)周期自由振蕩的存在。至今,用于觀測(cè)地球自由振蕩的儀器主要包括：寬頻帶地震儀、基線應(yīng)變儀、重力儀、水管傾斜儀、垂直擺傾斜儀和水位儀[2,4-15]。2004年 12月 26日的蘇門答臘大地震,給人類造成了重大傷亡和損失,同時(shí)也為研究地球自由振蕩提供了難得的機(jī)會(huì)。本文將利用云南省云龍地震臺(tái)(后面簡(jiǎn)稱云龍臺(tái))3個(gè)方向和山東省泰安基準(zhǔn)地震臺(tái)(后面簡(jiǎn)稱泰安臺(tái))2個(gè)方向伸縮儀的數(shù)字化觀測(cè)資料提取蘇門答臘地震激發(fā)的地球球型自由振蕩。

2 資料

2004年 12月 26日印尼蘇門答臘島西部附近海域印度洋發(fā)生M9.0地震,震中位置為 3.9°N、95.9°E,并在印度洋多個(gè)周邊國(guó)家引起大規(guī)模海嘯,造成了約 1 300 km長(zhǎng)的余震區(qū)[16]。

云南省云龍臺(tái)和山東省泰安臺(tái)使用的伸縮儀是中國(guó)地震局武漢地震科學(xué)儀器研究院研制的 SS-Y型數(shù)字化伸縮儀,是目前我國(guó)地震傾斜潮汐觀測(cè)臺(tái)網(wǎng)中使用的主要觀測(cè)儀器。針對(duì)這次地震,筆者研究了該伸縮儀的觀測(cè)數(shù)據(jù),發(fā)現(xiàn)云龍臺(tái)三分量和泰安臺(tái)二分量均可以較好地識(shí)別出球型振蕩。圖 1是云龍臺(tái)三分量和泰安臺(tái)二分量伸縮儀記錄的蘇門答臘 9.0級(jí)大地震分鐘值曲線,此曲線沒(méi)有經(jīng)過(guò)任何消除干擾的處理,從圖 1中可清晰地看到固體潮、地震波和同震階躍。分析采用的數(shù)據(jù)云龍臺(tái)從地震發(fā)生后開始,提取了有明顯地震影響的 6天的資料,共8 640分鐘;泰安臺(tái)因?yàn)榈诹觳糠秩睌?shù),采用 5天的數(shù)據(jù),共 7 200分鐘。

圖 1 云龍臺(tái)和泰安臺(tái)伸縮儀記錄的印尼蘇門答臘大地震分鐘值曲線Fig 1 Minute-value curves recorded by the digital extensometer at Yunlong,Taian seis mostations excited by theM9.0 Sumatra-Andaman great earthquake on 26 December,2004

3 處理方法

研究采用通常功率譜密度估計(jì)方法提取地球自由振蕩[14-19],即

其中,自相關(guān)函數(shù) Rn為:

式中,N =8 640為所用數(shù)據(jù)的數(shù)目,x為臺(tái)站記錄的伸縮儀數(shù)據(jù)。Sk離散值為功率譜密度值。為了消除數(shù)據(jù)不能無(wú)限長(zhǎng)而必須加窗造成的對(duì)功率譜密度估計(jì)的影響,本文采用 Hanning窗抑制旁瓣,突出主瓣。

4 計(jì)算結(jié)果

云龍臺(tái)三分量和泰安臺(tái)二分量伸縮儀的功率譜密度 (PSD)估計(jì)如圖 2～4所示 (垂直虛線表示PREM模型給出在其頂點(diǎn)標(biāo)出振型的自由振蕩頻率值)。圖 2為 0.28～1.30 MHz頻段的功率譜密度估計(jì),在這一頻段可以看到,云龍臺(tái)三分量和泰安臺(tái)二分量伸縮儀的觀測(cè)數(shù)據(jù)可以清楚地檢測(cè)到球型振蕩的基型振蕩0S3、0S4、0S5、0S6、0S7,同時(shí),還檢測(cè)到諧頻球型振蕩2S1和環(huán)型振蕩0T5。“足球”式振蕩0S2振型雖然附近功率譜密度有一個(gè)峰值,但相對(duì)于周圍的噪聲水平并沒(méi)有明顯分開。作為初步研究結(jié)果,在后面的分析中,我們只用 PREM模型給出的自由振蕩的頻率值與云龍臺(tái)三分量和泰安臺(tái)二分量的功率譜密度進(jìn)行對(duì)比。

圖 2 0.28～1.3 MHz頻段地傾斜數(shù)據(jù)功率譜密度估計(jì)曲線Fig 2 Curves of power spectrum density esti mation of ground tilt data in the 0.28～1.33 MHz frequency band of the digital extensometer

圖 3為 1.3～2.4頻段范圍內(nèi)云龍臺(tái)三分量和泰安臺(tái)二分量的功率譜密度(相對(duì)值)與 PREM模型給出的球型振蕩頻率值的對(duì)應(yīng)情況?？梢钥闯?S8～0S15在三分量中均被清楚地檢測(cè)出來(lái)。

圖 4為 2.4～3.5 MHz頻段范圍內(nèi)云龍臺(tái)三分量和泰安臺(tái)二分量的功率譜密度(相對(duì)值)與 PREM模型給出的球型振蕩頻率值的對(duì)應(yīng)情況?？梢钥吹?雖然在這個(gè)頻段噪聲比較大,但云龍臺(tái)三分量和泰安臺(tái)二分量的觀測(cè)數(shù)據(jù)仍然清楚地檢測(cè)到球型基頻振蕩0S16～0S26。

圖 3 1.3～2.4 MHz頻段地傾斜數(shù)據(jù)功率譜密度估計(jì)曲線Fig 3 Curves of power spectrum density estimation of ground tilt data in the 1.33～2.4 MHz frequency band of the digital extensometer

圖 4 2.4～3.5 MHz頻段地傾斜數(shù)據(jù)功率譜密度估計(jì)曲線Fig 4 Curves of power spectrum density estimation of ground tilt data in the 2.4～3.5 MHz frequency band of the digital extensometer

通過(guò)圖 2～4得出：SS-Y型數(shù)字化伸縮儀的觀測(cè)資料可以提取到印尼蘇門答臘大地震激發(fā)的地球球型振蕩。

5 結(jié)論

1)云龍臺(tái)和泰安臺(tái) SS-Y型數(shù)字化伸縮儀的觀測(cè)資料可以提取出印尼蘇門答臘 9.0級(jí)大地震激發(fā)的球型自由振蕩,并能準(zhǔn)確地檢測(cè)到0S3～0S26的基型球型振蕩,同時(shí)還能檢測(cè)到諧頻球型振蕩2S1和環(huán)型振蕩0T5,并與地球初步參考模型 (PREM)的理論自由振蕩頻率的對(duì)比說(shuō)明,實(shí)測(cè)振蕩頻率與 PREM預(yù)測(cè)的振蕩頻率基本符合。

2)需要說(shuō)明的是,除了地震可以激發(fā)地球自由振蕩外,還有許多其他因素,如大氣擾動(dòng)[18,19]、隨季節(jié)變化的大氣壓力變化[20,21]也可引起地球自由振蕩。因此本文所檢測(cè)到的印尼蘇門答臘大地震激發(fā)的地球自由振蕩資料還可能有其他源的貢獻(xiàn)。

3)如何用觀測(cè)資料檢驗(yàn)伸縮儀觀測(cè)的頻響特性,還需要進(jìn)一步研究。

致謝 衷心感謝云南省地震監(jiān)測(cè)中心前兆臺(tái)網(wǎng)和山東省地震局泰安基準(zhǔn)臺(tái)提供數(shù)據(jù),以及盧雙玲、孟建國(guó)、萬(wàn)永革教授提供的幫助!

1 傅承義,陳運(yùn)泰,祁貴仲.地球物理學(xué)基礎(chǔ)[M].北京：科學(xué)出版社,1985.

2 Garland G D.Introduction to Geophysics(Mantle,Core and Crust).2nd ed[M].Toronto：W B Saunders Company, 2003.

3 Benioff H,Press F and Smith S.Excitation of the free oscillations of the earth by earthquakes[J].Geophys.Res,1961, 66(2)：605-619.

4 NessN R,Harrison C T and SliehterL B.Observation of the free oscillation of the earth[J].Geophys Res,1961,66：621 -629.

5 Alsop L E,Sutton G H and Ewing M.Free oscillation of Earth observed on strain and pendulum seismographs[J]. Geophy.Res,1961,66(2)：631-641.

6 Bogert B P.An observation of free oscillations of the Earth [J].Geophys.Res,1961,66(2)：643-646.

7 DziewonskiA M and Gillbert F.Observation of normalmodes from 84 recordings of the Alaskan earthquakes 0f 1964 March 28[J].Geophys.J.R.Astr.Soc,1972,27：393-446.

8 Bolt B A and Currie R G.Maximum entropy estimates of Earth torsional eigenperiods from 1960 trieste data[J].Geophys.J.R.Astr.Soc,1975,40：107-114.

9 Zadro M and Braitenberg C.Measurements and interpretations of tiltstrain gauges in seis mically active areas[J].Earth Science Reviews,1999,47：151-187.

10 呂永清,蔡亞先,周云耀.用 JCZ-1超寬頻帶數(shù)字地震儀觀測(cè)地球自由振蕩[J].地殼形變與地震,1997,(3)：103-105.

11 唐磊,邱澤華,闞寶祥.中國(guó)鉆孔體應(yīng)變臺(tái)網(wǎng)觀測(cè)到的地球自由振蕩[J].大地測(cè)量與地球動(dòng)力學(xué),2007,(6)：37 -44.

12 任佳,等.數(shù)字化水位儀觀測(cè)的蘇門答臘大地震激發(fā)的地球球型自由振蕩[J].地震研究,2009,32(4)：333-338.

13 Stein S and Okal E A.Speed and size of the Sumatra earthquake[J].Nature,2005a,434：581-582.

14 鄧娜,等.利用雙陽(yáng)臺(tái)水管觀測(cè)資料檢測(cè)地球自由振蕩[J].大地測(cè)量與地球動(dòng)力學(xué),2009,(增刊)：41-44,58.

15 Stein S and Okal E A.The 2004 Sumatra earthquake and India ocean tsunami：What happened and why[J].The Earth’s Scientist,2005b,(2)：6-11.

16 萬(wàn)永革,盛書中,周公威.中國(guó)數(shù)字地震臺(tái)網(wǎng)記錄的蘇門答臘—安達(dá)曼地震激發(fā)的地球球型自由振蕩的檢測(cè)[J].地震學(xué)報(bào),2007,29(4)：369-381.

17 雷湘鄂,等.蘇門達(dá)臘地震激發(fā)的地球自由振蕩及其譜線分裂的檢測(cè)與討論 [J].中國(guó)科學(xué) (D輯),2007,37 (4)：504-511.

18 邱澤華,等.鉆孔差應(yīng)變儀觀測(cè)的蘇門答臘大地震激發(fā)的地球環(huán)型自由振蕩[J].地球物理學(xué)報(bào),2007,50(3)：797-805.

19 萬(wàn)永革.數(shù)字信號(hào)處理的MATLAB實(shí)現(xiàn)[M].北京：科學(xué)出版社,2007.

20 Park J,et al.Earth’s free oscillations excited by the 26 December 2004 SumatraAndaman earthquake[J].Science, 2005,308：1 139-1 146.

21 Beroza G C and Jordan T H.Searching for slow and silent earthquakes using free oscillations[J].Geophys Res., 1990,95(B3)：2 485-2 510.

DETECTING FREE OSC ILLATION OF EARTH BY USE OF OBSERVATIONS OF EXTENSOM ETER

Yang Yuewen1),Ren Jia2),Yang Xing3),YangLingying3),He Yuquan1)and LiLihong1)

(1)L ijiang Seism ostation,Earthquake Adm inistration of Yunnan Province,Lijiang 674100 2)Zhangjiakou Central Seism ostation,Earthquake Adm inistration of Hebei Province,Zhangjiakou 075000 3)Earthquake Adm inistration of Yunnan Province,Kunm ing 650224)

W ith power spectrum density esti mation method,we extract the spheroidal free oscillation of fundamental-frequency(0S3～0S26)of the earth excited by theM9.0 Sumatra-Andaman great earthquake on 26 December 2004 from the data of digital extensometer at Yunlong seis mostation in Yunnan province and at Taian seis mostation in Shandong province without eliminating the earth’s tide.Comparing the period of the extracted oscillation with that of the theoreticaloscillation by using the Preliminary Reference EarthModel(PREM),it is found that they are accordantwith each other.

extensometer;M9.0 Sumatra earthquake;PREM model;free oscillation;spheroidal oscillation

1671-5942(2010)Supp.(Ⅰ)-0078-04

2010-07-21

云南省自然科學(xué)基金(2008CD220)

楊躍文,男,1969年生,工程師,主要從事地震監(jiān)測(cè)及研究.E-mail：72310@163.com

P315.5

A