肖 攀， 宋國(guó)華， 李露露， 吳 輝， 宋 馳， 張偉峰

(安徽省地震局蚌埠地震臺(tái)， 安徽 蚌埠 233000）

0 引言

蚌埠地震臺(tái)位于江淮地區(qū)， 淮河南岸， 蚌埠市西郊小黃山（高新技術(shù)開發(fā)區(qū)內(nèi)）， 郯廬斷裂帶西側(cè)， 臺(tái)基屬花崗巖。 1976年9月建臺(tái)， 1977年6月投入觀測(cè)， 1980年4月擴(kuò)建。距市區(qū)3 km， 海拔40 m， 屬安徽省12個(gè)區(qū)域有人值守?cái)?shù)字地震臺(tái)之一。 測(cè)震儀器現(xiàn)使用KS-2000 型寬頻帶地震計(jì)。 數(shù)字化測(cè)震儀較以往使用的模擬測(cè)震儀， 具有頻帶寬、 動(dòng)態(tài)范圍大等特性， 因此記錄到的地震數(shù)目比模擬記錄的多。 但往往在分析地震時(shí)， 地震波被一些低頻信號(hào)或者是高頻信號(hào)相互疊加干擾， 使得震相無(wú)法精確識(shí)別， 這些干擾嚴(yán)重影響地震臺(tái)測(cè)震記錄波形質(zhì)量。 為了能更好地提高觀測(cè)資料質(zhì)量的精度， 本文分析了這些干擾源的特性， 并制作出相應(yīng)的濾波器， 以便于解決實(shí)際工作中存在的問(wèn)題。

1 干擾源介紹

1.1 高頻干擾介紹

隨著蚌埠市社會(huì)經(jīng)濟(jì)的迅速發(fā)展， 城市化進(jìn)程的快速推進(jìn)， 蚌埠地震臺(tái)所處的地理位置已屬于規(guī)劃中高新技術(shù)開發(fā)區(qū)范圍內(nèi)。 依2006年高新技術(shù)開發(fā)區(qū)規(guī)劃局設(shè)計(jì)， 2012年初， 開發(fā)區(qū)開始施工修建45 m 六車道大慶路主干道， 離蚌埠地震臺(tái)僅180 m。 目前快車道已經(jīng)通車。

作為城市外環(huán)的主干道， 每天都有大量的載重汽車通過(guò)， 由于蚌埠地震臺(tái)測(cè)震儀器直接放置在基巖上， 無(wú)山洞保護(hù)， 因此載重汽車的通行對(duì)我臺(tái)的地震觀測(cè)環(huán)境產(chǎn)生了較大的影響（見圖1）。

圖1 載重汽車通過(guò)時(shí)的干擾波形Fig.1 Interference signals as the truck through

西山黃山村窯廠（距臺(tái)站直線距離不到1000 m）采用地表爆破進(jìn)行挖土， 對(duì)蚌埠地震臺(tái)觀測(cè)環(huán)境造成了很大的威脅。 由于窯廠工作的干擾， 測(cè)震記錄波形背景干擾增大， 小震記錄不清， 震相難以識(shí)別， 嚴(yán)重影響了測(cè)震的觀測(cè)質(zhì)量。 另外， 大量爆破頻頻發(fā)生， 易造成數(shù)字觀測(cè)系統(tǒng)不穩(wěn)定。 如果爆破與地震事件重合， 其波形還會(huì)疊加干擾， 給地震分析工作帶來(lái)一定的困難（見圖2）。

圖2 臺(tái)站記錄大量爆破波形數(shù)據(jù)Fig.2 Blasting waveforms recorded by Bengbu seismic station

1.2 低頻干擾介紹

地震臺(tái)站觀測(cè)的數(shù)據(jù)中除了地震震源激發(fā)的地震波外， 還包含了如地脈動(dòng)、 海浪干擾等低頻數(shù)據(jù)， 為地震波數(shù)據(jù)分析帶來(lái)麻煩[1]。 以我臺(tái)記錄的ML1.5 級(jí)近震的UD 分向記錄為例說(shuō)明（見圖3）。

2 干擾源分析

場(chǎng)地干擾背景水平是影響臺(tái)站地震監(jiān)測(cè)能力的主要原因， 因此， 有必要對(duì)地脈動(dòng)噪聲進(jìn)行測(cè)試和分析。 本文選取了蚌埠地震臺(tái)凌晨0 點(diǎn)和上午10 點(diǎn)這2個(gè)時(shí)段KS-2000 寬頻帶地震計(jì)連續(xù)波形作平均值噪聲功率譜分析（見圖4、 5）。 可以看出， 蚌埠地震臺(tái)的功率譜密度曲線介于地球的高噪聲模型和低噪聲模型曲線之間， 其RMS 值如表1 所示。

圖3 淹沒在地脈動(dòng)中的地震波Fig.3 Seismic wave covered by microtremor

圖4 00 點(diǎn)鐘噪聲功率譜Fig.4 Nosie power spectrum at 0 o'clock

圖5 10 點(diǎn)鐘噪聲功率譜Fig.5 Nosie power spectrum at at 10 o'clock

表1 蚌埠地震臺(tái)RMS 值Table 1 RMS of Bengbu seismic station

根據(jù)GB/T19531.1-2004 地震臺(tái)站觀測(cè)環(huán)境技術(shù)要求， 在干擾高峰期和平靜期， 蚌埠地震臺(tái)臺(tái)址均達(dá)到了國(guó)家C 類地區(qū)環(huán)境噪聲水平Ⅱ級(jí)標(biāo)準(zhǔn)。 但是從表1 中可以看出， 在干擾高峰期的0.008～0.1 Hz、 1～20 Hz 其干擾明顯大于平靜期。 根本原因在于低頻部分有地脈動(dòng)干擾， 高頻部分有爆破和載重汽車過(guò)境干擾。

3 干擾的去除

3.1 高頻干擾的去除

地震臺(tái)站觀測(cè)的數(shù)據(jù)中除了地震震源激發(fā)的地震波外， 還包含了如地方爆破等高頻干擾數(shù)據(jù)。 這里以蚌埠臺(tái)記錄的2011年3月12日23 時(shí)53 分日本大地震余震的UD 分向記錄為例說(shuō)明高頻干擾的去除。 爆破頻率相對(duì)地震波的頻率明顯較高， 去掉該爆破就需要設(shè)計(jì)低通濾波器， 把高頻爆破濾掉。 由于遠(yuǎn)震波形數(shù)據(jù)的地震波頻成分較低。 我們選擇低于1.5 Hz 的頻率成分即可滿足要求（濾波結(jié)果見圖6）。 圖6 為濾波前后的波形圖， 在270 s處有一個(gè)高頻爆破， 濾波后完全去除了高頻爆破的干擾， 使得震相更加清晰。

圖6 運(yùn)用低通濾波器濾波前后波形的比較Fig.6 Comparison of original and filtered signals by low-pass filter

蚌埠地震臺(tái)180 m 處， 目前正在進(jìn)行道路建設(shè)工程， 在施工過(guò)程中， 蚌埠地震臺(tái)的地震監(jiān)測(cè)受到了較大的干擾。 圖7 可見， 長(zhǎng)周期地震信號(hào)中疊加了高頻的載重汽車干擾， 嚴(yán)重影響到數(shù)據(jù)的分析和使用。 在頻譜圖中可以看出， 高頻干擾集中在25 Hz 左右。 因此我們可以用低通濾波器將高頻干擾濾除， 分析結(jié)果見圖8。

圖7 蚌埠臺(tái)受干擾的波形圖（a) 和Fourier 振幅圖(b)Fig.7 Seismic wave disturbed by road construction 180 meters away（a) and Fourier amplitude (b)

圖8 濾波前和濾波后的波形圖Fig.8 Comparison of original and filtered signals

通過(guò)輸出結(jié)果， 我們可以直觀看出， 完全去除了高頻干擾。

3.2 低頻干擾的去除

我們以蚌埠臺(tái)記錄的ML1.5 級(jí)近震的UD 分向記錄為例說(shuō)明。 在原始波形圖上， 有一個(gè)低頻的波浪式背景干擾與這個(gè)高頻地震波形疊加， 故需要設(shè)計(jì)高通濾波器， 把低頻的波浪式背景干擾波濾掉。 我們選擇高于0.2 Hz 的頻率成分即可， 濾波結(jié)果見圖9。

通過(guò)圖9 的對(duì)比， 我們可以很清晰的看到， 在濾波前， 地震信號(hào)完全被干擾信號(hào)所覆蓋， 根本分不清地震信號(hào)的到時(shí)， 濾波后可清楚分辨Pg 出現(xiàn)的時(shí)間應(yīng)該是47.2 s 處， 起到了預(yù)期的效果。

圖9 濾波前和濾波后的波形圖Fig.9 Comparison of original and filtered signals

4 結(jié)語(yǔ)和討論

數(shù)字濾波器的應(yīng)用十分廣泛， 運(yùn)用MATLAB 語(yǔ)言能夠便捷地設(shè)計(jì)出滿足要求的濾波器。 在排除數(shù)字地震記錄中的干擾波時(shí)， 希望盡量濾除干擾成分而保留正常的頻率成分。但需要注意的是， 在設(shè)計(jì)濾波器參數(shù)時(shí),一定要對(duì)干擾信號(hào)進(jìn)行分析,找準(zhǔn)干擾波頻率的邊界范圍,盡可能做到真實(shí)信號(hào)不丟失。

數(shù)字地震記錄的干擾源還有很多， 若能針對(duì)臺(tái)站的實(shí)際情況對(duì)各種干擾波進(jìn)行分析，找出其干擾頻率和特點(diǎn)， 做為臺(tái)站的干擾背景資料， 將對(duì)今后臺(tái)站的資料處理和應(yīng)用起到很大的幫助。

致謝： 本文中所設(shè)計(jì)的濾波器原型來(lái)自于萬(wàn)永革教授所著 《數(shù)字信號(hào)處理的MATLAB實(shí)現(xiàn)》 一書， 特此致謝。

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