天津某環(huán)境背景噪聲分析

李曄 郭巍 許可 曹付陽 馬超群 徐小遠

[摘要] ?連續(xù)波形檢查中發(fā)現(xiàn)，天津爾王莊（EWZ）等臺明顯受到某未知噪聲的干擾，且EWZ臺受影響程度最大，此噪聲常年如此，不隨晝夜時間發(fā)生變化，且已存在多年。長期觀測還發(fā)現(xiàn)，此噪聲每年會有幾天觀測不到，但時間并不確定。與通常的環(huán)境背景噪聲不同，此噪聲具有一些固定特征。本文對其從本身特征、影響范圍、頻譜特性、對地震的影響等幾個方面進行了研究。研究表明，噪聲不受天氣、季節(jié)等自然因素影響，具有隨距EWZ臺距離增大而振幅變小的變化特征，其主要頻率范圍在1～2 Hz之間。噪聲的影響范圍較廣泛，可被周圍12個臺站記錄到，最遠到達YGZ臺，距離EWZ臺約58 km。通過偏振分析，發(fā)現(xiàn)噪聲的傳播具有特定的方向性，初步判斷此噪聲具有固定源噪聲的特點。此噪聲會對地震記錄造成影響，使震相拾取率降低，甚至完全淹沒在噪聲中，通過帶阻頻率濾波會有所改善。

[關(guān)鍵詞] 噪聲源; 時頻域特性; 影響范圍; 偏振分析; 濾波

[DOI] 10.19987/j.dzkxjz.2023-165

0 ?引言

地震監(jiān)測中，地震臺的環(huán)境條件影響著臺址的背景噪聲水平，較好的噪聲水平，能提升信噪比，獲得優(yōu)質(zhì)的觀測記錄^[1]；而不理想的背景噪聲水平會影響地震臺網(wǎng)的監(jiān)測，降低監(jiān)測能力。地震噪聲能反應(yīng)人類活動及大氣活動，同時能反應(yīng)地球嗡鳴、地球脈動等信息。地震噪聲有很多種，按照地震噪聲的來源，可將地震噪聲分為海洋、風(fēng)雨等自然環(huán)境噪聲，車輛行駛、工業(yè)生產(chǎn)等人為噪聲，儀器自身噪聲等；按照噪聲的頻率大小可分為短周期噪聲（如激流等）和高頻噪聲（如機械運作等）；按照地震噪聲的相干性，可分為有一定頻率和視速度的相干噪聲和隨機噪聲^[2-3]。

天津地處首都圈中心地區(qū)，是集各種職能于一身的大型城市，隨著各行各業(yè)發(fā)展進一步加深，地震觀測條件也面臨著前所未有的挑戰(zhàn)^[4-6]。地震臺站干擾因素增多、背景噪聲水平降低已成為必然。特殊噪聲的出現(xiàn)，會降低地震記錄信噪比^[7-8]。近年來，越來越多的學(xué)者將地震環(huán)境背景噪聲作為研究對象^[9]，他們通過背景噪聲互相關(guān)的方法研究地層速度結(jié)構(gòu)^[10]，在研究中發(fā)現(xiàn)了固定干擾源的干擾，發(fā)現(xiàn)固定噪聲源干擾會影響背景噪聲互相關(guān)的結(jié)果^[11-13]，同時還通過一定的方法找到了固定噪聲源的位置^[14-16]。

地震編目過程中，通過連續(xù)波形檢查發(fā)現(xiàn)，天津爾王莊（EWZ）等臺，明顯受到某未知噪聲的干擾，此噪聲常年如此，且已存在較長時間。長期觀測還發(fā)現(xiàn)，此噪聲每年會有幾天觀測不到，但時間并不確定。本研究對此噪聲進行了分析，從噪聲對時間域記錄的影響，噪聲的頻率特性，噪聲的影響范圍，噪聲源的存在性，以及濾波對地震中震相拾取率是否有所提升幾個方面進行研究，加深對此噪聲的了解。

1 ?噪聲基本情況

天津地震臺網(wǎng)在地震觀測過程中，發(fā)現(xiàn)爾王莊（EWZ）、糙甸（CAD）等臺站的連續(xù)波形記錄中存在干擾。多年觀察表明，此干擾已存在較長時間，多數(shù)情況下，不隨晝夜時間而發(fā)生明顯變化，振幅在夜間會略大于白天，周期較均勻（圖1）。同時，此干擾不會隨著季節(jié)更替及天氣變化而變，且連續(xù)不斷，一直出現(xiàn)在環(huán)境背景噪聲中。噪聲對三分向背景噪聲數(shù)據(jù)的影響并不一致，對垂直向的影響最大，對兩個水平分向的影響略小，且兩水平分向幅值比較一致。調(diào)查發(fā)現(xiàn)，爾王莊臺站記錄到的噪聲幅值最大，故本研究將此噪聲命名為EWZ噪聲。編目過程中發(fā)現(xiàn)，此噪聲會影響到地震中震相的拾取，進一步會降低本地區(qū)的地震監(jiān)測能力。

爾王莊地震臺位于天津市寶坻區(qū)，地處渤海灣盆地的北部邊界，與燕山造山帶相鄰^[17-18]，布設(shè)儀器為短周期地震計，在井下250 m左右處^[19]。為了排查噪聲的屬性，對EWZ臺站周圍100 km范圍內(nèi)環(huán)境進行調(diào)查，發(fā)現(xiàn)有水庫、河流、鐵路、自然保護區(qū)、農(nóng)場、石材城等場所的分布。

2 ?噪聲對地震記錄的影響

為了研究噪聲對地震時間域記錄的影響，選取本區(qū)域內(nèi)不同震級大小地震10條（表1）。通過對這些地震的觀察，發(fā)現(xiàn)對于附近100 km范圍內(nèi)的1級以下0.5級以上的較小地震，當震中距較小時，受影響程度小，隨著震中距變大，受影響程度有所增加，同時EWZ臺受影響程度較CAD臺大。在此范圍內(nèi)，隨著震級的增大、震中距的減小，兩臺受到噪聲干擾的情況會逐漸減輕。從圖2可以看出，爾王莊噪聲不僅會影響初動方向的拾取，還會影響地震到時類震相的拾取，使到時類震相無法辨認，甚至完全淹沒在噪聲中，使初至被混雜在噪聲中，無法判別；同時也影響振幅量取的準確性，進而降低定位精度，影響地震監(jiān)測能力。

為了進一步確定噪聲對地震的影響，這里將不同地震中EWZ臺和CAD臺的受影響情況進行了統(tǒng)計，本研究將地震記錄受到噪聲影響的程度分為10級，級別越高受影響程度越高，級別越低受影響程度越低（表1）。發(fā)現(xiàn)震級較小且震中距較大時，兩臺站的記錄情況不佳，受到噪聲影響的程度較大；當震級一定時，震中距越小記錄會越好；當震中距一定時，較大震級的地震會有較好的記錄。由于統(tǒng)計的地震多數(shù)在寶坻地區(qū)，震中距離CAD臺的距離要小于EWZ臺，同時噪聲源距離EWZ臺距離較近，使得多數(shù)地震中EWZ臺的受影響程度大于CAD臺。

3 ?噪聲頻譜分析

截取EWZ、CAD兩臺站夜間2小時的背景噪聲記錄，計算功率譜密度及振幅譜（圖3），從圖中可以看出，在1～2 Hz之間兩臺站都受到較明顯的噪聲干擾。噪聲的頻率并非為一固定值，而是在一定范圍內(nèi)。同時，EWZ臺受到的干擾程度較CAD臺大。通過對周圍其他臺站進行功率譜密度計算，發(fā)現(xiàn)ANK、ZTZ等十幾個臺都受到此頻段噪聲的影響。

圖3中采用了兩個臺站同時間段的背景噪聲數(shù)據(jù)，圖3b為三分向的振幅譜圖。從振幅譜圖中可以看出，背景噪聲中一組明顯的干擾，頻率在1～2 Hz之間，兩個臺站記錄到的噪聲的振幅有所不同，EWZ臺的振幅要大于CAD臺，可以判斷此噪聲源距離EWZ臺較CAD臺近。

4 ?爾王莊噪聲影響范圍

通過對EWZ臺周圍的臺站進行頻譜分析，發(fā)現(xiàn)除EWZ臺站外，此頻段噪聲還可被CAD、BAD、ANK等12個臺站記錄到，臺站分布情況如圖4所示，其中ANK臺已于2020年停止觀測，最遠距爾王莊臺的距離可達58 km，為YGZ臺。同時，為了更好的研究，在爾王莊臺與寶坻臺之間架設(shè)了兩個流動觀測臺，流動臺布設(shè)在地表，經(jīng)過了2個月的觀測，新架設(shè)流動臺也能記錄到此頻段信號，且距離EWZ臺近的流動臺振幅較大，進一步證明此噪聲在EWZ臺附近。

5 ?噪聲源位置固定性分析

觀察發(fā)現(xiàn)，爾王莊噪聲的振幅具有隨距離變化的特征，結(jié)合各臺站記錄到的噪聲的頻率特征及頻率一致性，排除了風(fēng)、海浪等非固定位置自然環(huán)境的干擾，初步判斷噪聲為固定源噪聲。為了進一步了解驗證，選取能記錄到此噪聲的臺站，對背景噪聲記錄作去均值、去傾斜、濾波等數(shù)據(jù)預(yù)處理，進行質(zhì)點振動分析，部分臺站計算結(jié)果如圖5所示。從圖中結(jié)果可以看出，此頻段內(nèi)噪聲的傳播，具有明顯的方向性，可進一步判斷噪聲為固定位置噪聲源。

6 ?地震數(shù)據(jù)的濾波處理

本研究針對表1中的地震，對受影響的臺站記錄進行1～2 Hz帶阻濾波處理，EWZ、CAD兩臺濾波前后對比圖如圖6所示。將濾波后的地震記錄的清晰情況也分為10級，級別越高記錄越清晰，震相越容易識別，級別越低震相越不易識別，見表1。對于發(fā)生在EWZ臺100 km范圍內(nèi)的地震，通過濾波，會使Pg震相、Sg震相有所凸顯，Sg震相會更加明顯。EWZ臺的震相拾取率會提升50%以上，CAD臺的震相拾取率會提升80%以上。除此之外，本研究還對其他受影響臺站進行了分析，發(fā)現(xiàn)通過對1～2 Hz頻帶噪聲進行濾波，會很大程度上提升震相拾取率，同時也會有助于初動方向的拾取。從結(jié)果可以看出，此濾波處理會在一定程度上提升震相的識別率，進一步提升地震的定位精度。

對于震級較大地震，采用普通頻率濾波，會影響地震信號原有振幅大小，如需較好的將此噪聲信號去除，還需進一步使用一些特殊的濾波方法。

7 ?結(jié)論與討論

本研究對EWZ臺附近的噪聲研究表明，此噪聲不同于通常意義上的環(huán)境背景噪聲。通常意義上的噪聲一般具有時間與位置上的隨機性，且不會造成太大規(guī)模的影響；而本研究中的噪聲明顯不是隨機的，它具有固定范圍的頻率特性和會隨著與EWZ臺距離的增大、影響逐漸變小的特點。布設(shè)在地表的流動臺，同樣能記錄到相同的噪聲信號，說明了此噪聲不是井下原因的影響。同時其地震記錄的粒子振動還具有明顯的方向性。此外，在我們不斷觀測的過程中發(fā)現(xiàn)，此噪聲會在每年中有不定期的幾天，是觀測不到的，這也充分說明了此噪聲的特殊性。

在工業(yè)及城市發(fā)展日益增強的今天，地震監(jiān)測受到了前所未有的挑戰(zhàn)，人類活動、工業(yè)生產(chǎn)越來越密集，尤其是對于像天津這樣較大型城市，如何在這種條件下保持或是提升地震監(jiān)測水平，已成為一個難題。同時，對于較大的干擾是否可以得到更多的利用，將其作為有用震源，也成為了我們需要研究的課題。

參考文獻

?

[1]中國地震局監(jiān)測預(yù)報司. 測震學(xué)原理與方法[M]. 北京：地震出版社，2017 ???Monitoring and Forecasting Department of China Earthquake Administration. Principles and methods of seismology[M]. Beijing：Seismological Press，2017

[2]彼得·鮑曼. 新地震觀測實踐手冊[M]. 中國地震局監(jiān)測預(yù)報司，譯. 北京：地震出版社，2006 ???Bormann P. New manual of seismological observatory practice[M]. Monitoring and Forecasting Department of China Earthquake Administration，trans. Beijing：Seismological Press，2006

[3]Kanai K. On microtremors Ⅷ[J]. Bulletin of Earthquakes Research Institute，1961，39：97-114

[4]孫路強，劉磊，卞真付，等. 利用時間序列算法實現(xiàn)測震數(shù)據(jù)存儲服務(wù)器磁盤容量預(yù)測[J]. 地震工程學(xué)報，2018，40（增刊1）：170-173，193 ???Sun L Q，Liu L，Bian Z F，et al. Disk capacity prediction for the seismic data storage server using the time series algorithm[J]. China Earthquake Engineering Journal，2018，40（S1）：170-173，193

[5]宋程，馬婷，畢金孟，等. 基于背景噪聲的環(huán)渤海地區(qū)波速變化初探[J]. 地震地磁觀測與研究，2022，43（增刊1）：18-20 ???Song C，Ma T，Bi J M，et al. Preliminary study on seismic wave velocity variation in the circum-Bohai region based on background seismic noise[J]. Seismological and Geomagnetic Observation and Research，2022，43（S1）：18-20

[6]王熠熙，邵永新，李悅，等. 基于多種方法的寶坻新井水位異常分析[J]. 地震，2020，40（1）：172-183 ???Wang Y X，Shao Y X，Li Y，et al. A study on Baodi well water level anomalies using multiple methods[J]. Earthquake，2020，40（1）：172-183

[7]張小艷，趙星，熊峰，等. 內(nèi)蒙古地傾斜觀測主要干擾信號的功率譜密度特征[J]. 華北地震科學(xué)，2022，40（3）：62-68 ???Zhang X Y，Zhao X，Xiong F，et al. Power spectral density characteristics of main interference signals of tilt observation in Inner Mongolia[J]. North China Earthquake Sciences，2022，40（3）：62-68

[8]楊亞運，汪建，傅卓，等. 臺基觀測方式對地震臺站背景噪聲影響分析[J]. 地震科學(xué)進展，2023，53（6）：241-250 ???Yang Y Y，Wang J，F(xiàn)u Z，et al. Influence analysis of observation mode for seismic station on background noise[J]. Progress in Earthquake Sciences，2023，53（6）：241-250

[9]李奇，張智，侯爵，等. 背景噪聲提取體波方法研究進展[J]. 地震科學(xué)進展，2021，51（10）：433-451 ???Li Q，Zhang Z，Hou J，et al. Research progress of the extraction body waves from ambient noise[J]. Progress in Earthquake Sciences，2021，51（10）：433-451

[10]王芳. 利用背景噪聲研究新疆呼圖壁儲氣庫周邊淺層介質(zhì)結(jié)構(gòu)及其變化[J]. 地震科學(xué)進展，2018，48（1）：45-46 ???Wang F. Shallow structure and temporal variation monitoring around Hutubi underground gas storage in Xinjiang based on seismic ambient noise[J]. Progress in Earthquake Sciences，2018，48（1）：45-46

[11]Shapiro N M，Campillo M，Stehly L，et al. High-resolution surface-wave tomography from ambient seismic noise[J]. Science，2005，307（5715）：1615-1618

[12]Gu Y J，Dublanko C，Lerner-Lam A，et al. Probing the sources of ambient seismic noise near the coasts of southern Italy[J]. Geophysical Research Letters，2007，34（22）：L22315

[13]孫和平，陳曉東，危自根，等. 基于地震和重力觀測的深地超寬頻帶背景噪聲初探[J]. 地球物理學(xué)報，2022，65（12）：4543-4554 ???Sun H P，Chen X D，Wei Z G，et al. A preliminary study on the ultra-wide band ambient noise of the deep underground based on observations of the seismometer and gravimeter[J]. Chinese Journal of Geophysics，2022，65（12）：4543-4554

[14]王爍帆，倪四道，王偉濤，等. 基于背景噪聲經(jīng)驗格林函數(shù)的地震準確定位精度分析?以2008年甘肅武都地震為例[J]. 地球物理學(xué)報，2022，65（8）：2904-2916 ???Wang S F，Ni S D，Wang W T，et al. Accuracy analysis of locating earthquakes with ambient noise empirical Greens functions：A case study of the 2008 Wudu earthquake （Gansu，China）[J]. Chinese Journal of Geophysics，2022，65（8）：2904-2916

[15]王偉濤，倪四道，王寶善，等. 中國中東部地震臺站噪聲互相關(guān)函數(shù)中面波前驅(qū)信號的分析研究[J]. 地球物理學(xué)報，2012，55（2）：503-512 ???Wang W T，Ni S D，Wang B S，et al. Studies on a persistent localized microseism source that produces precursors on noise correlation function observed using stations in central-eastern China[J]. Chinese Journal of Geophysics，2012，55（2）：503-512

[16]曲浩鑫，王鋒，周金玲，等. 海拉爾地震臺陣噪聲源調(diào)查[J]. 地震地磁觀測與研究，2020，41（6）：82-89 ???Qu H X，Wang F，Zhou J L，et al. Research on noise source in Hailar Seismic Array[J]. Seismological and Geomagnetic Observation and Research，2020，41（6）：82-89

[17]河北省地質(zhì)礦產(chǎn)局. 河北省 北京市 天津市區(qū)域地質(zhì)志[M]. 北京：地質(zhì)出版社，1989 ???Hebei Bureau of Geology and Mineral Resources. Regional geology of Hebei Province，Beijing Municipality and Tianjin Municipality[M]. Beijing：Geology Press，1989

[18]閆麗莉，高文晶，王建國，等. 震前衛(wèi)星紅外亮溫異常時空演化特征分析[J]. 地震工程學(xué)報，2021，43（1）：90-98 ???Yan L L，Gao W J，Wang J G，et al. Spatio-temporal evolution characteristics of satellite infrared brightness temperature anomalies before earthquakes[J]. China Earthquake Engineering Journal，2021，43（1）：90-98

[19]李曄，郭巍，高也，等. 3·21響水爆炸的地震學(xué)分析及譜比法研究[J]. 地震科學(xué)進展，2021，51（5）：215-222 ???Li Y，Guo W，Gao Y，et al. Seismic analysis and amplitude-style spectrum ratio study of the 21 March blast at Xiangshui[J]. Progress in Earthquake Sciences，2021，51（5）：215-222

Analysis of one ambient seismic noise in Tianjin

Li Ye^*， Guo Wei， Xu Ke， Cao Fuyang， Ma Chaoqun， Xu Xiaoyuan

Tianjin Earthquake Agency， Tianjin 300202， China

[Abstract] ?During checking continuous waveform record， it was found that some stations of Tianjin such as EWZ were significantly affected by one unknown noise， and EWZ station was the most seriously affected. Long term observations had also found that this noise cannot be observed for uncertainly several days each year. This noise was all year round and did not change with day and night time， and it had been existed for many years. This noise was different from typical environmental background noise and had some fixed characteristics. This article studied the noise from several aspects such as its own characteristic， its coverage of influence， spectral characteristics， and impact on earthquake records. The research has shown that the noise was not affected by natural factors such as weather and season， and had a characteristic of amplitude decreasing as the distance increasing from EWZ station. The main frequency of this noise ranged 1～2 Hz. The impact of noise was relatively wide and it can be recorded by 12 surrounding stations， with YGZ station being the farthest and approximately 58 km away from EWZ station. Through polarization analysis， it was found that the propagation of noise had a specific directionality， so it is preliminary judged that this noise had the characteristic of fixed source noise. This noise will have an impact on seismic records， reducing the seismic phase pickup rate， and even completely submerging it in the noise. Frequency filtering can improve it.

[Keywords] noise source; time frequency domain characteristics; influence scope; polarization analysis; filtering