漳州地震臺陣背景噪聲的特性研究

鄭重 郝春月 黃昭 張爽

1）中國地震局地球物理研究所，北京市海淀區民族大學南路5號 100081

2）福建省地震局地震災害防御中心，福州 350003

0 引言

建設一個新的臺陣／臺站時，首先應研究這個臺陣／臺站的噪聲場特性，因為噪聲狀況與該臺陣／臺站今后的信噪比密切相關。臺陣的噪聲場特性表現為臺陣勘址地點的噪聲日變化、噪聲水平、測點對間噪聲的相關性信息等，而研究目的是為了明確臺陣／臺站的基底背景噪聲，使地震監測、地震定位、地震數據分析以及更深入的地震研究工作得以順利地進行。早在20世紀70年代，人們便開始進行地震臺陣的噪聲結構（Hilmar et al，1971）以及地震臺臺址與噪聲源之間最小推薦距離的研究（Willmore，1979）。在最佳臺網布局的研究中，都要以信噪比為依據（Kijko，1977；Rabinowitz et al，1990；Steinberg et al，1995）。此類的研究還有Bormann等（1997）等對德國局域臺網（GRSN）中臺站的地震背景噪聲的研究，Harjes（1990）對中歐區域臺陣的噪聲研究等。近年來，Voulgaris等（2004）對希臘小孔徑臺陣、Radulian等（2005）對羅馬尼亞BURAR臺陣、Thomas等（2008）對意大利中部的臨時3分向小孔徑臺陣、郝春月等（2003、2006）對蘭州和那曲臺陣等分別進行了此方面的研究工作。由于漳州臺陣濱臨海域，為研究它與大陸內部地震臺陣所不同的噪聲場特性，本文主要對福建省漳州臺陣在2007年10～12月間的勘址測點進行噪聲場特性分析，并根據噪聲場特性進行噪聲源研究。應用背景噪聲均方根、噪聲功率譜、相干函數分析等研究方法力圖解決漳州臺陣噪聲場評估的諸方面問題。

1 數據資料

漳州臺陣測點布設在水庫周圍，附近是畬族居住區。在漳州臺陣的勘址過程中，布設了18個測點（圖1），臺間距為198～3031m。勘址過程中記錄了大量的原始波形數據，其中大部分數據是噪聲，這就為勘址地點噪聲場特性的研究提供了便利。

圖1 漳州臺陣測點分布及編號

2 噪聲場特性

2.1 各勘址測點的功率譜估計

噪聲功率譜是將白噪聲信號的平均功率分成各個頻率分量所占有的成分，再按照頻率大小依次畫出各頻率分量所占的功率。功率譜估計是信號處理的一項基本內容，它表明了被測信號的功率在各頻域隨頻率的分布情況，因此，又稱為功率譜密度。

2.1.1 Welch平均周期圖經典譜估計

Welch平均周期圖法是對直接法的改進，即把一長度為N的數據xN（n）分成L段（在分段時可允許每一段的數據有部分的重疊），每一段的長度為M。分別求每一段的功率譜，然后加以平均。每一段的功率譜可由下式表示（胡廣書，1997）

式中

U為歸一化因子，使用它是為了保證所得到的譜是漸近無偏估計。d（n）是漢寧窗（0.5±0.5cos（2pin／N））。這樣，幾段的平均功率譜可由下式決定

2.1.2 噪聲功率譜的計算與結果

本文對每個測點采用隨機選出的8個白天和8個晚間的數據樣品段進行計算，以便求出每個測點噪聲的平均功率譜。樣品數據段為30000個點，采用8192點長度的漢寧窗和50%的重疊進行計算。原始噪聲數據經過去均值、去傾等預處理后，進行分段、加窗、傅立葉變換等，最后去除儀器響應，獲得平均噪聲功率譜值（圖2）。圖2中，上、下2條虛線表示Peterson高低噪聲模型，該模型是Peterson（1993）根據世界范圍內75個數字地震臺在噪聲平靜期和噪聲高潮期測定的地動加速度功率譜密度而建立的，被稱為新低噪聲模型（NLNM）和新高噪聲模型（NHNM），是現今普遍被認同的最高和最低噪聲界線的標準。而舊模型也是Peterson于1963年建立的。圖2中的所測點都是垂直向地震記錄的噪聲功率譜結果。

由圖2可見，所有的測點顯示了相似的特征，即頻率小于0.2Hz時，噪聲功率譜值更接近高噪聲模型（NHNM）。我們知道，0.6～1.5Hz頻段主要由中周期的海洋風暴組成，而0.07～0.60Hz頻段主要為海洋風暴次生的脈動，由于漳州臺陣距海域較近，所以，小于0.2Hz的高噪聲主要是由海浪、風波等引起的次生脈動造成的。所有測點在1.5、5.0Hz左右分別有一個高的脈沖值，表示漳州臺陣周圍在這2個頻率段具有噪聲源。

總體來說，該陣址的低頻背景噪聲幅值較高，主要是因為選址地點緊靠大海，海浪對陸地拍擊的影響強烈，而高頻背景噪聲幅值尚可，在5.0Hz處具有噪聲源。

2.2 各測點的噪聲均方根值

均方根（rootmean square即RMS）是一段數據大小的量度，它給出了數據的量級。在許多科學領域，均方根是計算數據規模的行之有效的方法。在這里，為了找出各測點的各種噪聲源和噪聲的普遍變化規律，對各測點的均方根值作了計算。關于噪聲，從大的方面可分為短周期的人為噪聲和海洋、氣象因素引起的大于2s的較長周期噪聲（張誠，1986）。

2.2.1 均方根定義

均方根的定義是一段信號的振幅值取平方后，再求其在該時段內的平均，然后求其平方根。對于離散信號，均方根計算公式為

其中，n是樣本數，xi是第i個樣本的幅值。

2.2.2 各測點噪聲的均方根值

計算各測點噪聲的均方根值時采集了連續4天內的噪聲數據，對于每10m in的噪聲段，先對其進行濾波，然后根據式（4）計算其幅值的均方根值。濾波器是截止頻率為3.0Hz的高通濾波器和帶通濾波器（0.05～1.50Hz），它們分別對應人為噪聲干擾和海洋、氣象噪聲的主要頻率范圍（彼得·鮑曼，2006）。由于功率譜計算結果顯示在1.5、5.0Hz左右具有噪聲源，所以均方根值的計算還增加了1.0～2.0Hz、4.0～6.0Hz頻段的濾波結果。圖3給出了部分測點在2007年12月21～24日間的噪聲均方根值日變化。

圖2 各測點的平均噪聲功率譜

由圖3可見，各測點在0.05～1.50Hz之間的噪聲水平略微高于頻率大于3.0Hz時的噪聲水平，并且大部分測點在12：00或 12：00～18：00之間出現強振幅值。1.0～2.0Hz頻段的噪聲均方根值日變化曲線與0.05～1.50Hz頻段日變化曲線形狀相似（也就是在相同的時段出現強振幅值），但幅值偏低。經過3.0Hz高通濾波的噪聲均方根值，呈現出明顯的日變化，即白天時段的噪聲均方根值要比晚間的高，高出的幅值從10多個到上百個COUNT值不等。經過4.0～6.0Hz濾波的噪聲均方根值曲線形狀與經過3.0Hz高通濾波的噪聲均方根值曲線相似，也呈現了明顯的日變化，但幅值偏低。在12：00處也出現強峰值。

2.3 各測點對的噪聲相干函數

場地勘址時首先應該對當地的噪聲進行相關性分析，這是確定臺陣場地布局的關鍵點。本文主要研究信號在頻域中的相關性曲線，即相干函數。相干函數是2個信號在一個固定距離內的相關曲線，通過它可得出最小的噪聲相關距離，從而取得臺陣布局的依據。

圖3 2007年12月21～24日兩測點噪聲均方根值的日變化

2.3.1 相干函數原理

在頻域中，經常用到以功率譜表示的相干函數 rxy（f），相干函數rxy（f）（又稱凝聚函數）給出了2個隨機信號在頻率域的相似性，其中f表示2個隨機信號的頻率。一般來說，相干函數的值在0～1之間變化。若兩信號完全不相關，則相干函數值為零。對于2個相同的信號，其相干函數為 1。相干函數定義（Kulhánek，1973）為

其中，Gxx（f）和 Gyy（f）分別為 2個隨機信號 x（t）和 y（t）的自功率譜密度；Gxy（f）為它們的互譜密度。

2.3.2 勘址測點對間的噪聲相干函數

為了確定各測點之間的噪聲在頻域的相關性，并根據此相關性對理想的臺陣內外環半徑進行初步評估，對各測點的噪聲相干函數進行了計算與分析。

噪聲相干函數的計算采用了隨機選取的8個晚間和8個白天的10m in噪聲樣品段。利用4096點的漢寧窗（40.96s）和50%的重疊。勘址期間各子臺兩兩配對組成的臺間距共有30個，間距值為198～3031m。圖4為24個不同間距的噪聲相干函數圖，由圖4可見，隨著距離的增加，噪聲相干函數值在各頻率上逐漸減小。也就是說，測點對距離越近，噪聲相干性越好；距離越遠，噪聲相干性越差。距離為198～797m時，在0.1～1.0Hz之間的噪聲高度相干（接近于1），這個距離之后，0.1～1.0Hz頻段的噪聲相干性開始下降，到1700m左右，0.6Hz以上頻率的相干函數基本降到0.5以下。到1760m左右0.8Hz處的噪聲就基本不相干了，510m左右2.0Hz處的噪聲基本不相干。而遠震的優勢頻率段一般在0.8～2.5Hz，近震的優勢頻率段一般在2.0～10.0Hz。在不對信號進行相關性處理的情況下，可利用1760、510m這2個距離作為同心圓臺陣半徑的參考值。根據以上結果，內環半徑參考值定為510m，而外環半徑參考值定為1700m左右。

3 噪聲源分析

地震噪聲的季節變化主要是由自然界變化引起的，且噪聲周期明顯大于2s。海洋風暴產生的脈動頻譜峰值在周期T接近7s處，變化可達20dB。高頻噪聲主要來自人工源（如交通工具運行、機械施工等），通常有一個明顯的10～20dB的日變化（彼得·鮑曼，2006）。信噪比是評定臺址質量的主要參數，噪聲源在很大程度上影響同一頻段的地震信號，所以對于噪聲源的分析是非常重要的工作。

漳州臺陣陣址的噪聲數據處理結果表明，該臺陣存在明顯的噪聲源，主要的2個噪聲源分布在1.5、5.0Hz左右。下面分別利用3種方法對噪聲源進行分析。

3.1 噪聲功率譜方法

漳州臺陣各測點的噪聲功率譜在1.5、5.0Hz處的落差值如表1所示。功率譜是噪聲信號的功率隨頻率的分布情況，也就是能量分布情況，所以能量越強，功率譜值越高。表1表明，在2203測點處，1.5Hz處的噪聲源產生了比背景噪聲多10.5dB的能量，在所有測點里能量最強。在2101、2201、2202、2203測點組成的區域內和2301～2305組成的環內，噪聲源產生的能量均較大，大部分均比背景噪聲多8.5dB以上，預示著此區域就是1.5Hz噪聲源的發生地（而來自交通干道或其他人為噪聲和地質條件引起的噪聲干擾，皆會隨距離增加而減弱）。5.0Hz處的噪聲源在2402、2403、2404、和2405測點處產生了大小相等的強能量值，均比背景噪聲強17dB以上，其中，2402測點處產生了高于背景噪聲19.9dB的能量。這4個測點的連線方向為NW-SE方向，說明該噪聲源來自EN-WS方向，而在4個測點西南向的所有測點在5.0Hz處記錄到的信號功率譜值均較小，表明在5.0Hz處產生的噪聲源來自EN方向。

圖5為1.5、5.0Hz處噪聲最大范圍示意圖。由圖5可見，圓內（包括圓周）所有測點都在1.5Hz出現大振幅的噪聲干擾，根據這些測點的地理位置，估計是庫區大風引起的水浪拍打岸邊所致。外圍幾個測點從上到下分別為2404、2403、2402和2405，表示的是在5.0Hz單頻處噪聲分貝最大的臺站的連線，而其他測點記錄的5.0Hz噪聲功率譜均較小，預示著噪聲源來自EN方向。根據實地考察，認為5.0Hz噪聲源來自EN方向距臺陣中心點6.4km的發電廠。

圖4 測點對間的相干函數

表1 部分測點在2個頻率噪聲源處產生的功率譜差值

圖5 2個頻率處噪聲最大范圍示意圖

3.2 噪聲均方根值方法

根據噪聲均方根值的分析結果，在1.0～2.0Hz濾波的噪聲的日變化平穩，在12：00～18：00時段具有高強度脈沖。據勘址報告稱，勘址期間的下午都在刮大風，所以推斷1.0～2.0Hz的噪聲源有可能是大風引起的風浪拍打大壩所造成的。經過4.0～6.0Hz濾波的噪聲具有明顯的日變化，白天噪聲明顯升高，晚間噪聲明顯降低，呈現出清晰的日變化特征，表明5.0Hz噪聲源是人工源。

3.3 噪聲相干函數方法

根據噪聲相干函數分析的結果可知，最外圍測點（臺點編號第2位以4開頭命名的測點）之間組成的測點對在1.5Hz處的相干函數偏低，而包含其他測點（臺點編號第2位以1、2或3開頭命名的測點）的測點對在1.5Hz處的相干函數偏高，表明了噪聲源來自以1～3開頭命名的測點內部，如圖5所示。根據5.0Hz處的噪聲相干函數可知，與部分高分貝噪聲測點2402、2404、2406組成的測點對組合，相干函數明顯偏高，一般大于0.5，最大可達到0.8，而與低分貝測點2409組成的測點對組合，相干函數偏低，小于0.5。該結果也應與噪聲功率譜的結果相符合（圖5）。

4 結論

綜上所述，漳州臺陣陣址的低頻背景噪聲幅值較高，主要是因為勘址地點緊鄰海區，海浪對陸地拍擊的影響強烈，而高頻背景噪聲幅值為平均水平。大于3.0Hz的背景噪聲均方根值出現明顯的日變化特征。

以信噪比達到最大為原則，使用相干函數方法得出該臺陣內圓半徑參考值為510m，而外環半徑參考值為1700m左右。

利用場地噪聲進行噪聲源分析非常重要，這為今后的地震檢測、地震定位、地震數據分析及更深入的地震研究工作得以順利地進行提供了重要的基礎資料信息。

根據以上結果可判定，1.5Hz的噪聲源在3字頭測點區域內，主要是由大風引起的風浪拍打庫區大壩造成的。根據實地考察認為，引起單頻5.0Hz左右的噪聲源來自EN方向距臺陣中心點6.4km的發電廠。

致謝：福建地震災害預防中心和福建省地震局監測中心提供了勘址觀測數據，謹致謝忱！