高偉亮，呂 睿，梁 艷，丁大業(yè)，鄭樹平

（1.山西省地震局，山西 太原 030021；2.太原大陸裂谷動(dòng)力學(xué)國(guó)家野外科學(xué)觀測(cè)研究站，山西 太原 030025）

0 引言

山西測(cè)震臺(tái)網(wǎng)已運(yùn)行7年，目前共有57個(gè)測(cè)震臺(tái)站。隨著山西經(jīng)濟(jì)建設(shè)的飛速發(fā)展，礦山開采、大型工廠的建設(shè)、高速公路建設(shè)、鐵路建設(shè)等，都可能對(duì)地震觀測(cè)的場(chǎng)地環(huán)境造成影響。因此，有必要對(duì)全省測(cè)震臺(tái)站的臺(tái)基背景噪聲進(jìn)行一次全面的分析研究，對(duì)存在的問題提出科學(xué)、合理的解決方案，為山西測(cè)震臺(tái)網(wǎng)的臺(tái)站優(yōu)化工作提供一定的參考依據(jù)。

1 山西數(shù)字測(cè)震臺(tái)網(wǎng)簡(jiǎn)介

山西數(shù)字測(cè)震臺(tái)網(wǎng)由3個(gè)國(guó)家數(shù)字測(cè)震臺(tái)、38個(gè)區(qū)域數(shù)字測(cè)震臺(tái)、16個(gè)市縣數(shù)字測(cè)震臺(tái)，共57個(gè)臺(tái)站組成。全部為基巖臺(tái)站，其中有44個(gè)臺(tái)站為地表觀測(cè)室，13個(gè)臺(tái)站為觀測(cè)山洞。儀器設(shè)備主要由港震公司EDAS－24IP數(shù)據(jù)采集器和EDAS－24GN數(shù)據(jù)采集器、BBVS－60／120寬頻帶地震計(jì)和英國(guó) Guralp公 司CMG－DM24數(shù)據(jù)采集器、CMG－3ESPC寬頻帶地震計(jì)組成。

2 資料選取

“十五”山西數(shù)字測(cè)震臺(tái)網(wǎng)自2007年9月運(yùn)行以來(lái)，運(yùn)行穩(wěn)定。通常臺(tái)基的干擾背景與時(shí)間有關(guān)，白天噪聲大，夜間較平靜。為達(dá)到對(duì)全省測(cè)震臺(tái)站進(jìn)行綜合評(píng)價(jià)的目的，在資料選取時(shí)，白天和夜間分別選取樣品記錄。隨機(jī)選擇57個(gè)臺(tái)站2012年至2014年所記錄的寬頻帶噪聲數(shù)據(jù)，白天選取11時(shí)至13時(shí)，夜間23時(shí)至次日01時(shí)，并確保所選的數(shù)據(jù)是無(wú)震、無(wú)明顯干擾情況下的背景噪聲記錄，記錄長(zhǎng)度足夠長(zhǎng)（1小時(shí)）。

3 數(shù)據(jù)預(yù)處理

3.1 數(shù)據(jù)獲取

臺(tái)網(wǎng)的臺(tái)站記錄速度值采用SEED格式。通過MSDP模塊從臺(tái)網(wǎng)中心JOPENS系統(tǒng)的流服務(wù)器上讀取原始波形數(shù)據(jù)，本文以各個(gè)臺(tái)站單通道每小時(shí)的實(shí)時(shí)記錄波形為計(jì)算樣本，首先將1小時(shí)的波形數(shù)據(jù)進(jìn)行分道并轉(zhuǎn)換為ASCⅡ碼格式文本。

3.2 去直流偏移

在臺(tái)站觀測(cè)記錄中會(huì)出現(xiàn)直流偏移，并不代表真實(shí)的地脈動(dòng)，在時(shí)間運(yùn)算過程中必須剔除。把整個(gè)記錄長(zhǎng)度求和除以記錄的總數(shù)，即把記錄的平均值作為該記錄的直流偏離量，然后，逐點(diǎn)減掉直流偏離量作為校正過程第一步，用數(shù)學(xué)公式可以寫成：

式中：fi為每個(gè)采樣點(diǎn)的值；N為記錄長(zhǎng)度中的采樣點(diǎn)個(gè)數(shù)［1］。

3.3 速度PSD值估算

估算地震噪聲PSD值最為常用的方法被稱為直接傅里葉變換法。該方法通過計(jì)算有限長(zhǎng)度數(shù)據(jù)序列的FFT變換來(lái)計(jì)算PSD值。為了得到以頻率為自變量的FFT變換，y（t）的FFT變換可以表示為：

式中：f＝k／NΔt（k＝0，1，…，N－1）；N 為采樣點(diǎn)個(gè)數(shù)。文中的功率譜密度PSD表示為總功率譜密度的形式，即：

式中：f＝k／NΔt（k＝0，1，…，N／2），k的取值從0至N／2，僅包含Y（f）的正頻率成分［2］。

3.4 扣除儀器響應(yīng)

為得到地動(dòng)噪聲的物理量值，需對(duì)公式（3）計(jì)算的Pv（f）扣除儀器響應(yīng)，即：

式中：H（f）為儀器的傳遞函數(shù)［1］。

4 臺(tái)網(wǎng)各子臺(tái)的噪聲觀測(cè)資料計(jì)算結(jié)果

4.1 臺(tái)站地動(dòng)噪聲均方根值和有效測(cè)量動(dòng)態(tài)范圍計(jì)算

數(shù)字地震儀記錄的數(shù)據(jù)是電壓數(shù)字值，實(shí)際計(jì)算中要先把數(shù)字轉(zhuǎn)換成地動(dòng)速度，計(jì)算公式為：

式中：v為實(shí)測(cè)地脈動(dòng)速度值；N為實(shí)際記錄背景噪聲的數(shù)字?jǐn)?shù)；V0為模擬輸入峰值電壓；R為儀器分辨率（counts）；S1為地震計(jì)的工作靈敏度；Gm為數(shù)據(jù)采集器實(shí)際工作時(shí)的增益。

利用地動(dòng)噪聲的均方根值（RMS）可以衡量臺(tái)站的背景噪聲水平，用該值衡量噪聲水平的優(yōu)點(diǎn)在于對(duì)來(lái)自不同噪聲源的噪聲可按照一個(gè)相同的尺度進(jìn)行比較。地動(dòng)噪聲均方根值按下式進(jìn)行計(jì)算：

式中：n為測(cè)量次數(shù)；vi為某一點(diǎn)實(shí)測(cè)地動(dòng)速度值；ˉv為實(shí)測(cè)地動(dòng)速度均值。

根據(jù)數(shù)字地震觀測(cè)技術(shù)規(guī)范要求，臺(tái)基的干擾水平在1～20Hz內(nèi)，地動(dòng)速度的噪聲水平應(yīng)低于1×10－7m／s［2］。

4.2 環(huán)境地噪聲水平（Enl）計(jì)算

按照中國(guó)地震局測(cè)震臺(tái)站觀測(cè)環(huán)境技術(shù)要求（GB／T19531.1－2004）的規(guī)定，采用1／3倍頻程濾波器在1～20Hz頻帶范圍內(nèi)用PSD計(jì)算RMS值，如果噪聲記錄是地面運(yùn)動(dòng)速度進(jìn)入，則計(jì)算結(jié)果得到的就是環(huán)境地噪聲水平Enl。

由速度功率譜密度PSD計(jì)算RMS的公式為：

式中：fc為分度倍頻程中心頻率；RBW 為相對(duì)帶寬，計(jì)算公式為：

式中：fu為分度倍頻程上限頻率；fl為分度倍頻程下限頻率［3］。

臺(tái)站有效測(cè)量的動(dòng)態(tài)范圍反映了儀器本身性能指標(biāo)確定后，臺(tái)站環(huán)境干擾水平對(duì)記錄地震信號(hào)最大動(dòng)態(tài)范圍的限制。由于臺(tái)站環(huán)境噪聲與地震信號(hào)的疊加，使得儀器達(dá)不到設(shè)計(jì)的動(dòng)態(tài)范圍，只能達(dá)到有效測(cè)量的動(dòng)態(tài)范圍。有效測(cè)量的動(dòng)態(tài)范圍計(jì)算公式為：

式中：V0為模擬輸入峰值電壓；Gm為數(shù)據(jù)采集器實(shí)際工作時(shí)的增益；S1為地震計(jì)的工作靈敏度；RMS為地動(dòng)噪聲的均方根值［2］。

4.3 環(huán)境地噪聲的劃分方法

根據(jù)中國(guó)大陸背景地噪聲區(qū)域劃分，山西省數(shù)字測(cè)震臺(tái)網(wǎng)臺(tái)站所屬的背景地噪聲區(qū)域?yàn)锽類地區(qū)，安放寬頻帶數(shù)字地震計(jì)臺(tái)站的環(huán)境地噪聲水平應(yīng)不大于Ⅱ級(jí)環(huán)境地噪聲水平，安放甚寬頻帶數(shù)字地震計(jì)臺(tái)站的環(huán)境地噪聲水平應(yīng)不大于Ⅰ級(jí)環(huán)境地噪聲水平。按照中國(guó)地震局測(cè)震臺(tái)站觀測(cè)環(huán)境技術(shù)要求（GB／T19531.1－2004）規(guī)定，臺(tái)基背景噪聲在1～20Hz頻帶范圍的速度RMS值作為評(píng)估標(biāo)準(zhǔn)，RMS＜3.16×10－8m／s的臺(tái)站劃為Ⅰ類臺(tái)，3.16×10－8m／s≤RMS＜1.00×10－7m／s的臺(tái)站劃為Ⅱ類臺(tái)，1.00×10－7m／s≤RMS＜3.16×10－7m／s的臺(tái)站劃為Ⅲ類臺(tái)，3.16×10－7m／s≤RMS＜1.00×10－6m／s的臺(tái)站劃為Ⅳ類臺(tái)，1.00×10－6m／s≤RMS＜3.16×10－6m／s的臺(tái)站劃為Ⅴ類臺(tái)［3］。

4.4 山西測(cè)震臺(tái)網(wǎng)測(cè)震臺(tái)站環(huán)境地噪聲水平計(jì)算

對(duì)57個(gè)數(shù)字測(cè)震臺(tái)站進(jìn)行地動(dòng)噪聲功率譜密度計(jì)算，由于臺(tái)基情況比較復(fù)雜，在短周期頻段1～20Hz范圍內(nèi)，地動(dòng)噪聲水平較低；在10～1 000s頻段，地動(dòng)噪聲水平較高。因此，把臺(tái)基按地噪聲高低分類，根據(jù)觀測(cè)頻段劃分比較合適，按全頻段劃分比較困難。計(jì)算得到57個(gè)臺(tái)站白天和夜間時(shí)段RMS、有效測(cè)量的動(dòng)態(tài)范圍和地噪聲環(huán)境水平如表1所示。

表1 山西測(cè)震臺(tái)網(wǎng)各臺(tái)站RMS、有效測(cè)量動(dòng)態(tài)范圍和環(huán)境地噪聲水平Table 1 The root－mean－square，the effective dynamic range and background noise level of each station in Shanxi seismic network

表1 （續(xù)）

表1 （續(xù)）

5 噪聲功率譜分析

5.1 速度PSD值估算

令時(shí)間函數(shù)f（t）表示地震記錄中的位移，對(duì)其進(jìn)行傅里葉變換：

將其變換到頻率域中得F（ω），實(shí)函數(shù)｜F（w）｜2通常稱為功率譜或能量譜，更確切地叫功率譜密度或能量譜密度，即在頻率的尺度上每單位間隔的功率或能量。計(jì)算功率譜的方法通常有2種：快速傅里葉變換法（FFT）和相關(guān)變換法。一種是首先計(jì)算給定信號(hào)的快速傅里葉變換，然后取其絕對(duì)值的平方；另一種是首先計(jì)算自相關(guān)，然后計(jì)算其傅里葉變換。本文采用第一種方法，利用快速傅里葉變換法（FFT）計(jì)算功率譜。公式如下：

式中：N 的取值為2k，即2，4，8，16，…。

5.2 山西測(cè)震臺(tái)網(wǎng)各臺(tái)基噪聲功率譜密度曲線分析

對(duì)所選的資料進(jìn)行計(jì)算處理后，得出山西數(shù)字測(cè)震臺(tái)網(wǎng)各子臺(tái)的地動(dòng)噪聲有效值和功率譜。計(jì)算功率譜密度和地動(dòng)噪聲有效值的程序采用中國(guó)地震局港震機(jī)電技術(shù)有限公司童汪練老師研制開發(fā)的軟件。下面將山西測(cè)震臺(tái)網(wǎng)各臺(tái)基噪聲功率譜密度曲線形態(tài)分幾種情況分析。

（1）臺(tái)基噪聲正常，一致性非常好的臺(tái)站。

從右玉、上皇莊臺(tái)臺(tái)基背景噪聲功率譜密度曲線形態(tài)上可以看出，垂直、東西和南北向的功率譜曲線從長(zhǎng)周期到高頻段的一致性非常好，而且3條曲線在20s～20Hz頻段內(nèi)非常貼近地球地噪聲模型（NLNM），屬于正常臺(tái)基的地動(dòng)噪聲功率譜形態(tài)（見圖1）。

（2）臺(tái)基噪聲大致正常，低頻端噪聲峰有疑問，高頻端儀器噪聲明顯的臺(tái)站。

從長(zhǎng)治、恒山臺(tái)臺(tái)基背景噪聲功率譜密度曲線形態(tài)上可以看出，垂直、東西和南北向的功率譜曲線從長(zhǎng)周期到高頻段的一致性要差一些，而且3條曲線在低頻段內(nèi)有一個(gè)噪聲峰，高頻段內(nèi)噪聲毛刺多，初步判定為儀器噪聲（見圖2）。

（3）臺(tái)基噪聲大致正常，低頻段東西向和南北向噪聲偏高的臺(tái)站。

從安澤、永濟(jì)臺(tái)臺(tái)基背景噪聲功率譜密度曲線形態(tài)上可以看出，低頻段東西向和南北向噪聲偏高。分析認(rèn)為，臺(tái)站附近可能存在長(zhǎng)周期干擾源（見圖3）。

（4）臺(tái)基噪聲大致正常，低頻段垂直向噪聲偏高。

從婁煩、偏關(guān)臺(tái)臺(tái)基背景噪聲功率譜密度曲線形態(tài)上可以看出，2個(gè)臺(tái)大致屬于正常臺(tái)基噪聲功率譜形態(tài)，但低頻段垂直向噪聲比水平向偏高，而且垂直向噪聲遠(yuǎn)離水平向噪聲功率譜曲線（見圖4）。

圖1 右玉、上皇莊臺(tái)臺(tái)基背景噪聲功率譜密度曲線Fig.1 Power spectral density of background noise at Youyu and Shanghuangzhuang Station

圖2 長(zhǎng)治、恒山臺(tái)臺(tái)基背景噪聲功率譜密度曲線Fig.2 Power spectral density of background noise at Changzhi and Hengshan Station

圖3 安澤、永濟(jì)臺(tái)臺(tái)基背景噪聲功率譜密度曲線Fig.3 Power spectral density of background noise at Anze and Yongji Station

圖4 婁煩、偏關(guān)臺(tái)臺(tái)基背景噪聲功率譜密度曲線Fig.4 Power spectral density of background noise at Loufan and Pianguan Station

（5）選取靈丘臺(tái)為例，分析臺(tái)基背景噪聲功率譜密度曲線年變特征。

從靈丘臺(tái)臺(tái)基背景噪聲功率譜密度曲線可以看到，2012年3月6日山西信息安全服務(wù)工程升級(jí)改造安裝儀器后，5月的臺(tái)基背景噪聲功率譜曲線高頻段毛刺多，而且三分向一致性較差，表明儀器剛開始運(yùn)行1個(gè)多月來(lái)不穩(wěn)定，可以視為儀器噪聲，直至2013年11月曲線形態(tài)的一致性才得到改善。從2014年11月臺(tái)基背景噪聲功率譜曲線看到三分向的臺(tái)基噪聲功率譜曲線形態(tài)一致性非常好，說明儀器運(yùn)行三年多來(lái)，趨于穩(wěn)定（見圖5）。

圖5 靈丘臺(tái)2012—2014年臺(tái)基背景噪聲功率譜密度曲線變化圖Fig.5 Variation of power spectral density of background noise at Lingqiu Station from 2012to 2014

（6）以保德臺(tái)為例，分析當(dāng)臺(tái)站環(huán)境遭到破壞時(shí)的年變特征。

從第23頁(yè)保德臺(tái)臺(tái)基背景噪聲功率譜密度曲線圖6看到，2012年5月，臺(tái)基噪聲功率譜密度曲線低頻段水平向接近于地球高噪聲模型（NHNM）范圍，到11月，低頻段水平向已遠(yuǎn)遠(yuǎn)超過地球高噪聲模型（NHNM）范圍，到2014年11月，低頻段三分向均遠(yuǎn)遠(yuǎn)超過地球高噪聲模型（NHNM）范圍，說明臺(tái)站周圍環(huán)境隨著時(shí)間的推移在逐步惡化。現(xiàn)場(chǎng)勘查發(fā)現(xiàn)，保德測(cè)震臺(tái)距黃河不到1 000m，雨季黃河水量變化時(shí)會(huì)有長(zhǎng)周期變化的干擾；離臺(tái)站100m左右的地方打了一口抽水井，水泵工作時(shí)會(huì)有一定程度的高頻干擾；忻保高速公路距臺(tái)站不到2 000m，在施工過程中，臺(tái)站附近山梁的一些活動(dòng)和大型機(jī)械打鉆開采煤層氣對(duì)臺(tái)站的觀測(cè)數(shù)據(jù)都有一些干擾。臺(tái)站的水平向與垂直向在低頻段的各項(xiàng)數(shù)據(jù)都相差甚大，一致性較差，分析認(rèn)為是由于臺(tái)站氣流、氣壓、溫度等變化大的原因。臺(tái)站的臺(tái)基噪聲功率譜密度曲線三分向在高頻段毛刺多，而且接近于地球高噪聲模型（NHNM）范圍，認(rèn)為與抽水井水泵工作和大型機(jī)械打鉆開采煤層氣有關(guān)。

圖6 保德臺(tái)2012—2014年臺(tái)基背景噪聲功率譜密度曲線變化圖Fig.6 Variation of power spectral density of background noise at Baode Station from 2012to 2014

（7）臺(tái)基噪聲不正常的臺(tái)站。

從襄垣、霍州臺(tái)臺(tái)基背景噪聲功率譜密度曲線圖7看出，2個(gè)臺(tái)的低頻端水平向噪聲高過地球高噪聲模型（NHNM）范圍，說明臺(tái)站附近可能存在固定干擾源，也可能跟臺(tái)站當(dāng)?shù)氐牡刭|(zhì)構(gòu)造或固定場(chǎng)地環(huán)境有關(guān)，在今后做進(jìn)一步的研究。

圖7 襄垣、霍州臺(tái)臺(tái)基背景噪聲功率譜密度曲線Fig.7 Power spectral density of background noise at Xiangyuan and Huozhou Station

5.3 山西數(shù)字測(cè)震臺(tái)網(wǎng)噪聲功率譜密度分析

將山西數(shù)字測(cè)震臺(tái)網(wǎng)57個(gè)臺(tái)站垂直向的噪聲功率譜密度曲線繪制到一張圖上進(jìn)行比較（見圖8），得出：

圖8 山西省數(shù)字測(cè)震臺(tái)網(wǎng)子臺(tái)地動(dòng)噪聲譜曲線圖Fig.8 Ground pulsation noise spectrum of substations of Shanxi Digital Seismic Network

（1）長(zhǎng)周期和高頻段噪聲功率譜曲線差異大。分析地噪聲源主要是傾斜、氣流、氣壓、溫度等因素，表現(xiàn)在噪聲模型20s以上。對(duì)于風(fēng)擾問題，建議增強(qiáng)測(cè)震臺(tái)站密封性，減少測(cè)震臺(tái)站及周圍因大風(fēng)而產(chǎn)生震動(dòng)的因素。

（2）10～5s頻段的噪聲變化差異較小。

（3）20Hz以上頻段的噪聲功率譜呈線性上升趨勢(shì)的可視為儀器噪聲。

（4）低頻端噪聲高的原因主要有兩個(gè)，一是環(huán)境溫度變化大，二是儀器噪聲大。對(duì)于溫差、濕度大的問題，地表型測(cè)震臺(tái)站一般日溫差都較大，要增強(qiáng)測(cè)震臺(tái)站觀測(cè)室的密封性，未安裝擺罩的，盡快安裝。對(duì)于濕度較大的臺(tái)站，適時(shí)開門通風(fēng)干燥；山洞臺(tái)站應(yīng)在山洞內(nèi)采取一定的干燥措施，所有臺(tái)站均應(yīng)在擺罩內(nèi)安放干燥劑。

6 結(jié)論與討論

通過對(duì)山西測(cè)震臺(tái)網(wǎng)各臺(tái)基背景噪聲數(shù)字化記錄進(jìn)行計(jì)算和分析，得到臺(tái)基背景噪聲地動(dòng)速度均方根值（RMS）、有效測(cè)量動(dòng)態(tài)范圍、噪聲信號(hào)功率譜結(jié)果，表明：

（1）除婁煩臺(tái)、左權(quán)臺(tái)、臨汾臺(tái)、霍州臺(tái)、陽(yáng)城臺(tái)外，其余臺(tái)基地動(dòng)速度RMS值均達(dá)到區(qū)域數(shù)字測(cè)震臺(tái)網(wǎng)布設(shè)原則，即臺(tái)基在1～20Hz觀測(cè)頻帶內(nèi)，干擾地動(dòng)速度的平均水平（RMS值）應(yīng)低于1.00×10－7m／s。

（2）白天人為等干擾比夜間較強(qiáng)，地脈動(dòng)噪聲RMS值白天高于夜晚。

（3）13個(gè)山洞測(cè)震臺(tái)相對(duì)于地表基巖臺(tái)，噪聲RMS值相對(duì)小，受人為因素干擾小。

（4）山西測(cè)震臺(tái)網(wǎng)各臺(tái)基噪聲背景優(yōu)勢(shì)頻率各有特征，既有低頻噪聲干擾，又有高頻噪聲干擾。

（5）從表1看出，大同、山自皂、右玉、鎮(zhèn)川、岢嵐、保德、太谷、侯馬、古交、陽(yáng)泉、應(yīng)縣、神池、壽陽(yáng)、臨縣和中陽(yáng)的相對(duì)噪聲差值都超過了100%，5Hz頻點(diǎn)觀測(cè)動(dòng)態(tài)范圍也明顯變小。說明這些臺(tái)站白天噪聲相比晚上的差異太大，分析原因?yàn)槿藶榛顒?dòng)干擾。

（6）把各個(gè)臺(tái)站的環(huán)境地噪聲水平分析作為日常工作的一部分。當(dāng)臺(tái)站環(huán)境地噪聲水平發(fā)生明顯變化而周圍觀測(cè)環(huán)境無(wú)大的變化時(shí)，可能是儀器出現(xiàn)故障，要及時(shí)維修，可以作為檢測(cè)儀器狀態(tài)的一種手段。

（7）開展臺(tái)基分析工作，及時(shí)掌握臺(tái)站觀測(cè)環(huán)境的變化，為臺(tái)站的觀測(cè)環(huán)境保護(hù)提供充分依據(jù)，可以有針對(duì)性地對(duì)部分臺(tái)站進(jìn)行環(huán)境改造。

由于各子臺(tái)所在地的經(jīng)濟(jì)發(fā)展、城市規(guī)模擴(kuò)大、礦山的開采、大型工程的建設(shè)以及高速公路、鐵路的建設(shè)等，使部分臺(tái)基觀測(cè)條件受到影響，建設(shè)高質(zhì)量的測(cè)震臺(tái)站面臨著越來(lái)越大的阻礙和困難。為了減少噪聲的影響，建設(shè)山洞臺(tái)和深井臺(tái)成為未來(lái)建臺(tái)的趨勢(shì)。通過對(duì)山西數(shù)字測(cè)震臺(tái)網(wǎng)各種臺(tái)基類型的臺(tái)站噪聲和有效測(cè)量動(dòng)態(tài)范圍、功率譜的計(jì)算與分析，可以對(duì)各臺(tái)基類型臺(tái)站和設(shè)備作出客觀的評(píng)價(jià)，對(duì)山西省地震局“十二五”期間新建臺(tái)站的各項(xiàng)指標(biāo)設(shè)計(jì)（如臺(tái)基選擇、點(diǎn)位選擇、儀器選型等）提供理論依據(jù)，為山西省主動(dòng)震源技術(shù)發(fā)展和應(yīng)用提供背景資料。因此，對(duì)數(shù)字測(cè)震臺(tái)的臺(tái)基噪聲進(jìn)行分析是一項(xiàng)有意義和重要的工作。

［1］徐嘉雋.福建測(cè)震臺(tái)網(wǎng)測(cè)震臺(tái)站環(huán)境地噪聲水平用［J］.地震地磁觀測(cè)與研究，2013，34（1）：101－105.

［2］裴曉，尹繼堯，楊庭春.上海遙測(cè)臺(tái)網(wǎng)各類型臺(tái)基背景噪聲分析［J］.地球物理學(xué)進(jìn)展，2012，27（5）：1 897－1 903.

［3］中國(guó)地震局，云南省地震局.GB／T19531.1－2004 地震臺(tái)站觀測(cè)環(huán)境技術(shù)要求［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社，2004.