董 蕾，楊立明

（中國(guó)地震局蘭州地震研究所，甘肅蘭州 730000）

汶川大震前數(shù)字地震儀位移地脈動(dòng)低頻異?，F(xiàn)象研究①

董 蕾，楊立明

（中國(guó)地震局蘭州地震研究所，甘肅蘭州 730000）

選取了2008年4月25日至5月12日間的四川臺(tái)網(wǎng)14個(gè)臺(tái)寬頻帶地震儀記錄到的地脈動(dòng)波形資料，把速度波形做適當(dāng)?shù)臄?shù)據(jù)校正和濾波處理后，經(jīng)積分轉(zhuǎn)換成位移波形。頻譜分析結(jié)果表明：地震發(fā)生前存在低頻異常顫動(dòng)，地震前2～3天開始0．1～0．3Hz頻段傅里葉譜幅值逐漸增大；各個(gè)臺(tái)的位移振幅值逐漸增強(qiáng)。其中靠近龍門山斷裂帶的臺(tái)站記錄到的地脈動(dòng)變化較明顯。

汶川地震；低頻波；地脈動(dòng)；位移波形；傅里葉譜

Abstract：Through selection the broadband digital seismic data from April 25，2008to May 12，2008，recorded by fourteen seismic stations in Sichuan province，the microtremor wave form records are transformed to displacement wave forms after the data correction and filtering．The results of Fourier analysis indicate that before the earthquake occurred the low－frequency abnormal quiver is found；Fourier spectrum value at 0．1～0．3Hz band is increased in two to three days．Amplitude values of displacement waves are gradually increased at every seismic stations．The microtremor records change greater at stations near the Longmenshan fault zone．

Key words：Whenchuan earthquake；Low－frequency wave；Microtremor；Displacement wave；Fourier spectrum

0 引言

地震儀記錄到的不僅是地震事件，還記錄到了沒有地震時(shí)的微小地面運(yùn)動(dòng)，它是由風(fēng)、城市等的高頻噪聲，海洋波和周期長(zhǎng)些的天氣過(guò)程效應(yīng)產(chǎn)生的波所造成的［1］。

Longuet－Higgins［2］對(duì)長(zhǎng)周期地脈動(dòng)進(jìn)行了研究，所研究的長(zhǎng)周期地脈動(dòng)又分為原生地脈動(dòng)和次生地脈動(dòng)。Oliver和Page［3］首次注意到了原生地脈動(dòng)和次生地脈動(dòng)與同樣的大氣干擾的關(guān)系，觀測(cè)到了原生地脈動(dòng)的優(yōu)勢(shì)周期是次生地脈動(dòng)兩倍；在陸地上觀測(cè)到的原生地脈動(dòng)頻率在40～80mHz；次生或雙倍頻率的地脈動(dòng)是經(jīng)常觀測(cè)到的，優(yōu)勢(shì)峰值頻率在100～160mHz，或接近2倍的海浪峰值頻率。之后人們?cè)陂L(zhǎng)周期事件中發(fā)現(xiàn)了慢地震事件［4－5］。并在一些大震發(fā)生前幾天，發(fā)現(xiàn)了長(zhǎng)周期形變波［5］。Astiz［1］比較了全球臺(tái)網(wǎng)的100個(gè)不同臺(tái)站的典型噪聲水平，發(fā)現(xiàn)脈動(dòng)峰值在5～8s的周期出現(xiàn)。楊立明［6－7］利用連續(xù)數(shù)字儀速度地脈動(dòng)記錄的資料，研究發(fā)現(xiàn)震前臨震階段0．1～0．3Hz頻率的幅值出現(xiàn)快速、持續(xù)的增大現(xiàn)象，該現(xiàn)象具有重現(xiàn)性。根據(jù)以上的研究結(jié)果，本文將地震臺(tái)的速度脈動(dòng)記錄轉(zhuǎn)成位移記錄，并嘗試基于位移地脈動(dòng)波形記錄新的分析方法，探討大震前低頻段的地脈動(dòng)變化，觀察汶川地震前地脈動(dòng)的走向在時(shí)空上的變化情況，以及地脈動(dòng)與大震發(fā)生的關(guān)系。

1 資料和方法

1．1 資料的來(lái)源

本文資料來(lái)源四川臺(tái)網(wǎng)中心。四川地震臺(tái)網(wǎng)由52個(gè)數(shù)字地震臺(tái)站構(gòu)成，使用CMG－3ESPC、CTS－1E等寬頻帶地震儀，觀測(cè)結(jié)果穩(wěn)定。采樣間隔為0．01s，Nyquist頻率是50Hz。本文選取了14個(gè)臺(tái)站，都圍繞震中相對(duì)對(duì)稱分布，大致分布分在龍門山斷裂帶以西，以東及斷裂帶上。臺(tái)站的位置如圖1所示，震中距見表1。

表1 臺(tái)站與震中的距離

圖1 臺(tái)站和震中分布圖Fig．1 Distribution of the selected seismic stations and the epicenter of Wenchuan earthquake．

1．2 位移積分

寬頻帶數(shù)字地震儀具有記錄頻帶寬、分辨率高、動(dòng)態(tài)范圍大的特點(diǎn)。數(shù)字化地震儀基本上采用速度平坦型設(shè)計(jì)。數(shù)字地震學(xué)中，使用位移型、速度型、加速度型的說(shuō)法主要是指儀器的頻率特性，地面運(yùn)動(dòng)可以是位移、速度和加速度。數(shù)字地震學(xué)中提到數(shù)字記錄只是一種數(shù)字，單位是“數(shù)字?jǐn)?shù)”［8］。而從數(shù)字記錄恢復(fù)到位移波形一直是地震學(xué)中的一個(gè)基本問題，位移波形是確定震源深度、破裂面和地震矩等震源參數(shù)以及反演震源破裂過(guò)程的基礎(chǔ)［9－10］。

本文使用寬頻帶速度地震儀的數(shù)字地震記錄，采用積分轉(zhuǎn)換計(jì)算來(lái)獲得地面運(yùn)動(dòng)的位移波形。時(shí)間域的一次積分運(yùn)算相當(dāng)于頻率域除以圓周率，則原始記錄低頻成分會(huì)放大，產(chǎn)生低頻漂移，考慮到這個(gè)問題在積分前必須做濾波處理。

利用Matlab中復(fù)合梯形公式函數(shù)cumtrapz進(jìn)行速度與位移的轉(zhuǎn)換計(jì)算，得到位移波形記錄。并利用快速傅里葉變換（FFT）求取頻率特征。

1．3 地脈動(dòng)速度的計(jì)算

初始速度記錄的數(shù)據(jù)是evt格式的文件，由二進(jìn)制數(shù)組成，經(jīng)由earthquake data transform＿v3．0軟件轉(zhuǎn)換成文本格式的文件，所得數(shù)值的單位是“數(shù)字?jǐn)?shù)”。除以系統(tǒng)靈敏度后可得近似的真實(shí)的地動(dòng)速度（沒有除去儀器響應(yīng)和場(chǎng)地響應(yīng)）。

地動(dòng)速度（μm／s）＝數(shù)字?jǐn)?shù)（counts）／

系統(tǒng)靈敏度（counts／μm／s

1．4 記錄的基線校正與濾波

由于地震記錄中含有直流分量及長(zhǎng)周期漂移，直流分量在很大程度上影響FFT變換和積分結(jié)果，所以在進(jìn)行FFT變換和積分前，要進(jìn)行水平基線校正。水平基線校正簡(jiǎn)單地是將所有記錄點(diǎn)減去它們的平均值。由于0．005Hz的低頻限可以濾掉由積分造成的長(zhǎng)周期漂移的低頻波［10］，并且考慮到本文所涉及的頻段，用巴特沃斯4階高通濾波器濾掉小于0．05Hz的低頻波。

2 結(jié)果分析

2．1 汶川臺(tái)連續(xù)記錄頻率特征

本文以長(zhǎng)周期地脈動(dòng)為對(duì)象。長(zhǎng)周期的脈動(dòng)源是如雨、海浪、地質(zhì)內(nèi)力作用等的自然因素，它的優(yōu)勢(shì)周期在2～40s之間?？紤]到儀器的幅頻特性和積分等問題，本文所選的頻段是0．05～0．5Hz之間。圖2為汶川臺(tái)東西向位移波形的一些整點(diǎn)傅立葉譜，可以看出：（1）地脈動(dòng)的位移波形的譜峰值有3個(gè)，且對(duì)于不同的時(shí)間點(diǎn)，它所對(duì)應(yīng)的頻率段不一樣，但大致分別是0．05～0．1Hz，0．1～0．16Hz，0．16～0．3Hz。（2）從5月11日開始，各頻段傅氏譜的振幅值增大明顯。

對(duì)其他臺(tái)站的分析也一樣。分析發(fā)現(xiàn)，遠(yuǎn)震及小震對(duì)脈動(dòng)的0．05～0．1Hz頻率段影響較大，這段頻率段的幅值不穩(wěn)定。本文以下選取0．1～0．3Hz頻率段進(jìn)行分析。

圖2 汶川8．0地震前汶川臺(tái)EW道一些整點(diǎn)脈動(dòng)位移記錄0．05～0．5Hz頻段的頻譜圖Fig．2 Fourier spectrums in 0．05～0．5Hz of microtremor displacement records in EW channel at Whenchuan seismic station before Whenchuan MS8．0earthquake．

汶川地震前4月25日到5月12日間全球共發(fā)生了5次大于5級(jí)的地震，其中2008年5月2日9時(shí)33分發(fā)生在安德烈亞諾夫群島的7級(jí)地震和2008年5月8日0時(shí)45分發(fā)生在日本本州東海岸近海的7．1級(jí)地震對(duì)本文所選的臺(tái)站記錄到的波形有很大影響。除此之外，還有一些小震和儀器故障等因素造成脈動(dòng)頻譜的突變，所以扣除了5月2日 9時(shí)，10時(shí)，5月8日0時(shí)，1時(shí)這些時(shí)間段的記錄。

為了更清楚展示地震前4月25日到5月12日間脈動(dòng)傅里葉譜值變化特征，求取每小時(shí)0．1～0．3 Hz頻段內(nèi)幅度值包絡(luò)線的極大值，并擴(kuò)大1 000倍作為該時(shí)段的頻率幅值，最后得到汶川臺(tái)整點(diǎn)傅里葉譜峰值的變化進(jìn)程，如圖3所示?？梢钥闯雒}動(dòng)序列傅里葉譜峰值在10日后開始增大。

圖3 汶川8．0地震前汶川臺(tái)三通道整點(diǎn)脈動(dòng)序列傅里葉譜峰值0．1～0．3Hz頻段變化進(jìn)程Fig．3 Change of Fourier spectrum peak values in 0．1～0．3Hz of microtremor displacement in three channels records at Whenchuan station before Whenchuan MS8．0earthquake．

2．2 震前各臺(tái)連續(xù)脈動(dòng)記錄頻譜特征的時(shí)間進(jìn)程

經(jīng)過(guò)處理，選取的寬頻帶地震儀脈動(dòng)記錄的UD方向震前整點(diǎn)記錄傅里葉譜0．1～0．3Hz頻段峰均值隨時(shí)間變化曲線如圖4，其中，圖4（a）為龍門山斷裂帶西邊的臺(tái)站，圖4（b）為龍門山斷裂帶上的臺(tái)站，圖4（c）為龍門山斷裂帶東邊的臺(tái)站。可以看出，汶川地震前從5月10日開始，每個(gè)臺(tái)的脈動(dòng)的幅值都有不同程度的持續(xù)增大，而斷裂帶附近的臺(tái)站所記錄到的地脈動(dòng)幅值增大明顯。

圖4 震前分區(qū)各臺(tái)連續(xù)脈動(dòng)記錄頻譜特征的時(shí)間進(jìn)程Fig．4 The time histories of spectrum characteristic of continuous microtremor displacement records from stations in three areas before the Wenchuan earthquake．

2．3 脈動(dòng)位移波形的連續(xù)記錄

為了直觀的看出地震前脈動(dòng)的變化，本文選取了成都臺(tái)5月7日到5月12日震前的波形記錄，對(duì)其進(jìn)行積分得到位移波形（圖5）?？梢钥闯?，從10日開始脈動(dòng)幅值增大。

圖5 汶川地震前5月7日到12日成都臺(tái)記錄的地脈動(dòng)位移波形三分量Fig．5 The microtremor displacement wave forms in three channels recorded at Chengdu station from May 7to May 12，2008，before the Wenchuan earthquake．

類似的分析，得到其他臺(tái)的10日到12日震前脈動(dòng)位移波形，如圖6所示。

可以看出，每個(gè)臺(tái)的脈動(dòng)幅值都有不同程度的增加，斷裂帶以西的臺(tái)站增加幅度相對(duì)較弱。

3 結(jié)論與討論

（1）本文在分析了一些臺(tái)站的地脈動(dòng)位移波形之后，發(fā)現(xiàn)震前有低頻前兆信息，臨近地震的發(fā)生，脈動(dòng)在0．1～0．3Hz頻段的傅里葉幅值越來(lái)越大。震前2到3天內(nèi)脈動(dòng)位移波形振幅值持續(xù)增加?？梢钥闯鑫灰撇ㄐ蔚恼穹蹬c0．1～0．3Hz頻段的傅里葉幅值是對(duì)應(yīng)增加的。

圖6 一些臺(tái)站5月10日到震前的NS道地脈動(dòng)位移波形變化Fig．6 The change of the microtremor displacement wave forms in NS channel from stations in different areas．

（2）關(guān)于脈動(dòng)異常已有大量的研究。楊又陵［5］研究發(fā)現(xiàn)2001年11月14日昆侖山口西8．1級(jí)地震前發(fā)生了緩慢運(yùn)動(dòng)事件，并認(rèn)為該慢地震事件是8．1級(jí)大震前中昆侖斷裂臨震預(yù)滑動(dòng)引起的。許昭永等［11］通過(guò)巖石高溫高壓破裂實(shí)驗(yàn)與有機(jī)玻璃破裂試驗(yàn)認(rèn)為慢地震是臨界或亞臨界破裂或預(yù)先滑動(dòng)所致，是低頻應(yīng)變波動(dòng)；同時(shí)也發(fā)現(xiàn)并非所有巖石主破裂前都有明顯的臨界破裂或預(yù)滑動(dòng)現(xiàn)象，所以并非每個(gè)地震前都有慢地震。Linde等［12－13］認(rèn)為慢地震是地震斷裂過(guò)程的一個(gè)組成部分，在地震成核作用過(guò)程中可能起著重要的作用——是導(dǎo)致地震發(fā)生的觸發(fā)因素。張晁軍等［12］認(rèn)為慢地震可能消除一個(gè)正趨于斷裂的地震威脅，也可能轉(zhuǎn)移應(yīng)力而觸發(fā)一個(gè)正常地震。

然而在2008年5月8日至12日期間，西太平洋上產(chǎn)生了強(qiáng)臺(tái)風(fēng)Rammasun（威馬遜），在較早的長(zhǎng)周期地脈動(dòng)信號(hào)研究中認(rèn)為脈動(dòng)信號(hào)的產(chǎn)生與海洋中的暴風(fēng)系統(tǒng)有關(guān)。對(duì)于地脈動(dòng)波形的產(chǎn)生有2種不同的產(chǎn)生機(jī)制：（1）在wiechert理論（Wiechert曾提出脈動(dòng)是在激烈的碎浪對(duì)陡峭的海岸作用下產(chǎn)生的）中的海浪（1905，1907）；（2）根據(jù)Longuet－Higgins的理論（海浪波動(dòng)非線性理論）提出的海洋駐波模型［2，4］。胡小剛［15］在研究后認(rèn)為路徑和強(qiáng)度類似于Rammasun的臺(tái)風(fēng)很難在中國(guó)內(nèi)陸引起明顯的異常擾動(dòng)，因此，對(duì)于汶川地震震前四川的臺(tái)站記錄到的低頻脈動(dòng)異常所產(chǎn)生的原因需做進(jìn)一步的探討。

（3）本文首次將脈動(dòng)速度記錄轉(zhuǎn)換為位移波形記錄，并對(duì)此進(jìn)行傅里葉分析，結(jié)果與楊立明［7］利用速度記錄發(fā)現(xiàn)汶川地震前甘肅臺(tái)網(wǎng)存在地脈動(dòng)特殊頻率波動(dòng)的異?，F(xiàn)象是一致的。

本文是在前人工作的基礎(chǔ)上進(jìn)行研究的，對(duì)于臨震前的低頻波動(dòng)異?，F(xiàn)象的研究和應(yīng)用還需要大量進(jìn)一步的工作，本文的研究?jī)H是初步的，對(duì)此問題的深入研究有助于提高臨震預(yù)報(bào)的認(rèn)識(shí)水平。

致謝：感謝四川臺(tái)網(wǎng)中心提供寬頻帶地震儀記錄到的數(shù)據(jù)。

［1］ Peter M Shearer，著．陳章立，譯．地震學(xué)引論［M］．北京：地震出版社，2008：168－169．

［2］ Longuet－Higgins M S．A theory of the origin of microseisms［J］．Phil．Trans．R．Soc．London，1950，A 243：1－35．

［3］ Robert K Cessaro．Sources of Primary and Secondary Microseisms［J］．Bulletin of the Seismological Society of America，1994，84（1）：142－148．

［4］ 王迪晉，李正媛，呂品姬．慢地震研究綜述［J］．大地測(cè)量與地球動(dòng)力學(xué)，2007，27：21－25．

［5］ 楊又陵，趙根模，高國(guó)英，等．2001年11月14日昆侖山口西M8．1地震前的緩慢地震事［J］．國(guó)際地震動(dòng)態(tài)，2003，（9）：1－4．

［6］ 楊立明．汶川地震臨震地脈動(dòng)記錄特殊頻率波動(dòng)現(xiàn)象及其重現(xiàn)性初究［J］．國(guó)際地震動(dòng)態(tài)，2009，（1）：14－19．

［7］ 楊立明，王建軍，馮建剛，等．汶川地震前地脈動(dòng)低頻波動(dòng)現(xiàn)象及其應(yīng)用的初步研究［J］．中國(guó)地震，2009，25（4）：356－366．

［8］ 陳運(yùn)泰，吳忠良，王培德，等．?dāng)?shù)字地震學(xué)［M］．北京：地震出版社，2000：16－18，52－53，57－59．

［9］ 李鴻吉．用FFT和現(xiàn)代控制論方法恢復(fù)地面位移［J］．地球物理學(xué)報(bào)，1992，35（1）：37－43．

［10］ 周云好，陳章立，繆發(fā)軍．2001年11月14日昆侖山口西MS8．1地震震源破裂過(guò)程研究［J］．地震學(xué)報(bào)，2004，26（7）：9－20．

［11］ 許昭永，楊潤(rùn)海，胡毅力．慢地震慢機(jī)制的研究［J］．地震，2003，23（2）：12－20．

［12］ 張晁軍，石耀霖，馬麗．慢地震研究中的一些問題［J］．中國(guó)科學(xué)院研究生學(xué)報(bào)，2005，22（3）：258－269．

［13］ Linde A T，Gladwin M T，Johnston M J S．A slow earthquake sequence on the San Andreas fault［J］．Nature，1996，383：65－68．

［14］ HasselmannK A．A statistical analysis of the generation of microseisms［J］．Rev．Geophys．，1963，1：177－209．

［15］ 胡小剛，赫小光．強(qiáng)臺(tái)風(fēng)對(duì)汶川大地震和昆侖山大地震“震前擾動(dòng)”影響的分析［J］．地球物理學(xué)報(bào)，2009，5（7）：1363－1375．

Research on Low－frequency Anomalies of Microtremor Displacement Recorded by Digital Seismigraph before the Whenchuan Earthquake

DONG Lei，YANG Li－ming

（Lanzhou Institute of Seismology，China Earthquake Administration，Lanzhou 730000，China）

P315．31

A

1000－0844（2012）03－0268－06

10．3969／j．issn．1000－0844．2012．03．0268

2011－05－10

中國(guó)地震局蘭州地震研究所論著編號(hào)：LC2012022

董 蕾（1986－），女（漢族），碩士，研究方向?yàn)榈卣鹁C合預(yù)測(cè)理論．