新冠疫情防控前后湖北恩施地震臺(tái)背景噪聲變化分析

吳麗慧, 廖武林, 黎金玲, 趙學(xué)彧(. 湖北第二師范學(xué)院, 湖北 武漢 43005; . 湖北省地震局, 湖北 武漢 43007)

0 引言

隨著社會(huì)經(jīng)濟(jì)的發(fā)展和人口的增長(zhǎng),人類(lèi)活動(dòng)強(qiáng)度急劇增加,活動(dòng)范圍快速擴(kuò)張,地震計(jì)觀測(cè)到的人為噪聲(cultural noise)受到越來(lái)越多學(xué)者的關(guān)注。人為噪聲即人類(lèi)活動(dòng)引起的環(huán)境噪聲,主要集中于短周期頻段和高頻頻段,尤其是高頻頻段[1]。人類(lèi)活動(dòng)引起的噪聲源大多是分散或移動(dòng)的,來(lái)自各個(gè)方向的波疊加形成相當(dāng)復(fù)雜且近似穩(wěn)定的隨機(jī)噪聲場(chǎng),以高頻面波的方式傳播,隨距離和深度迅速衰減[2]。地震計(jì)將城市中人類(lèi)日?；顒?dòng)記錄為連續(xù)的噪聲信號(hào),且具有高頻特征[3]。在人口密集地區(qū),人類(lèi)活動(dòng)產(chǎn)生的噪聲通常在晚上和周末較弱,白天和工作日較強(qiáng)[4],各地法定節(jié)日期間呈現(xiàn)出較安靜狀態(tài),比如國(guó)慶節(jié)、新年和中國(guó)農(nóng)歷新年[5]。Wu等[6]在高校建立了一個(gè)小型校園微震臺(tái)網(wǎng),獲得了為期一個(gè)月的連續(xù)地震波形記錄數(shù)據(jù),發(fā)現(xiàn)地震高頻噪聲強(qiáng)度曲線(xiàn)清晰地展示了學(xué)生日?；顒?dòng)的強(qiáng)度時(shí)序規(guī)律。新冠疫情暴發(fā)以來(lái),地震背景噪聲再次受到地震學(xué)界廣泛關(guān)注,大量研究結(jié)果發(fā)現(xiàn)由于全球各地執(zhí)行疫情防控封鎖措施(lockdown),各地地震背景噪聲出現(xiàn)了不同比例的下降[5,7-8]。Lecocq等[5]研究發(fā)現(xiàn)在新冠疫情防控封鎖期間,人類(lèi)活動(dòng)的減少導(dǎo)致地震高頻噪聲在2020年3—5月期間平均下降了50%,認(rèn)為這段安靜期很可能是由防控封鎖措施、服務(wù)業(yè)和工業(yè)關(guān)閉、旅游和旅行人數(shù)下降等總的影響造成的,并發(fā)現(xiàn)這段時(shí)期是有地震噪聲記錄以來(lái)最長(zhǎng)、最明顯的地震背景噪聲安靜時(shí)期。

自20世紀(jì)50年代,地震儀器中記錄的地震噪聲就引起了地震學(xué)家的廣泛興趣[9]。深入認(rèn)識(shí)地震儀器記錄的地震噪聲時(shí)變特征,可獲得對(duì)地震臺(tái)站勘選和降噪、地震觀測(cè)技術(shù)與儀器性能及人類(lèi)活動(dòng)的深入認(rèn)識(shí)。起初,地震噪聲模型主要聚焦于地球本底背景噪聲。近年來(lái),隨著大數(shù)據(jù)技術(shù)的快速發(fā)展和地震數(shù)據(jù)的逐漸公開(kāi)共享,地震觀測(cè)記錄的人類(lèi)活動(dòng)噪聲受到關(guān)注,并被用于城市規(guī)劃、防災(zāi)減災(zāi)、公共安全管理等領(lǐng)域。

系統(tǒng)研究地震噪聲并量化其強(qiáng)度的一個(gè)標(biāo)志性成果是高/低噪聲模型的提出。Peterson[10]基于SRO和ASRO地震觀測(cè)數(shù)據(jù),通過(guò)人工挑取較安靜時(shí)段的觀測(cè)數(shù)據(jù),首先定量計(jì)算了各個(gè)臺(tái)站在0.1～100 s頻帶的噪聲強(qiáng)度,提出了低噪聲模型(LNM)。隨后,通過(guò)融合其他數(shù)據(jù),這個(gè)噪聲模型的適用頻帶被展寬。并在此基礎(chǔ)上得出了高噪聲模型(HNM)。

基于全球75個(gè)地震臺(tái)站噪聲觀測(cè)數(shù)據(jù),Peterson[11]更新了全球噪聲模型,通過(guò)計(jì)算選定地震觀測(cè)數(shù)據(jù)的功率譜密度(Power Spectral Density,PSD),得出新的全球噪聲模型,新高噪聲模型(NHNM)和新低噪聲模型(NLNM)。通過(guò)全球比對(duì)試驗(yàn)和觀測(cè)實(shí)踐,此模型已成為衡量臺(tái)站噪聲水平的重要依據(jù)和客觀標(biāo)準(zhǔn)。在傳統(tǒng)方法中,通常選擇相對(duì)平靜的一段地震記錄,通過(guò)計(jì)算其功率譜密度評(píng)估臺(tái)站噪聲水平[11]。該方法雖然在臺(tái)站背景噪聲水平的評(píng)估中獲得了廣泛的應(yīng)用,但在選擇噪聲記錄方面通常受一定主觀因素的影響,難以客觀反映臺(tái)站噪聲水平的整體特征。在此基礎(chǔ)上,McNamara等4]提出了基于概率密度的地震噪聲功率譜評(píng)估方法,即地震噪聲概率密度函數(shù)(power density function,PDF)。該方法使用連續(xù)地震觀測(cè)數(shù)據(jù),無(wú)需剔除諸如地震、人為活動(dòng)等時(shí)窗段,通過(guò)概率形式獲得對(duì)臺(tái)站噪聲變化特征的準(zhǔn)確把握。

2020年1月底,新冠疫情在以武漢為首的湖北省內(nèi)蔓延。根據(jù)疫情防控需要,2020年1月23日武漢封城,此時(shí)正值中國(guó)春運(yùn)高峰期并臨近中國(guó)農(nóng)歷春節(jié)(1月25日),人員跨區(qū)域流動(dòng)量大。1月24日恩施州封城,同日恩施州內(nèi)新增11名新冠肺炎感染者。圖1為疫情期間恩施州新冠肺炎每日新增及累計(jì)感染人數(shù)。恩施針對(duì)新冠疫情采取了嚴(yán)格的疫情防控措施,具體措施如下:

(1) 1月24日0時(shí)起,暫停恩施州城區(qū)公交、室內(nèi)公共文化娛樂(lè)場(chǎng)所等營(yíng)運(yùn)營(yíng)業(yè)活動(dòng)。

圖1 恩施州新冠肺炎每日新增及累計(jì)感染人數(shù)(數(shù)據(jù)來(lái)源:湖北省衛(wèi)生健康委員會(huì)官網(wǎng)統(tǒng)計(jì))Fig.1 Daily new and cumulative infected persons of COVID-19 in Enshi (Data source: Health Commission of Hubei Province)

(2) 1月26日0時(shí)起,城鄉(xiāng)道路客運(yùn)、鄉(xiāng)村公交、跨市客運(yùn)、旅游及包車(chē)客運(yùn)暫停營(yíng)運(yùn)。

(3) 2月19日0時(shí)起,進(jìn)一步強(qiáng)化新冠肺炎疫情防控措施,強(qiáng)化車(chē)輛通行管控,強(qiáng)化城鄉(xiāng)社區(qū)(村)封閉管理,所有城鄉(xiāng)居民小區(qū)實(shí)行24小時(shí)封閉管理。

(4) 2月底,根據(jù)湖北省新型冠狀病毒感染肺炎疫情防控指揮部發(fā)布的疫情管控文件,“省內(nèi)各類(lèi)企業(yè)先按不早于3月10日24時(shí)前復(fù)工”,引導(dǎo)涉及民生保障的相關(guān)工業(yè)企業(yè)分批分期復(fù)工復(fù)產(chǎn)。

(5) 3月15日起,恩施實(shí)行安全有序“放”政策,實(shí)行分區(qū)分級(jí)差異化防控策略,撤除鄉(xiāng)(鎮(zhèn)、辦、處)和村、組交通卡口,有序恢復(fù)城區(qū)出租車(chē)、市鄉(xiāng)客運(yùn)班車(chē)、村村通客車(chē)營(yíng)運(yùn),有序復(fù)工復(fù)產(chǎn)。

(6) 3月25日0時(shí)起,恩施州解除管控,恩施機(jī)場(chǎng)復(fù)航。

湖北恩施地震監(jiān)測(cè)中心站位于湖北省恩施市土橋壩五峰山,距離恩施城約3 km,始建于1976年,是全球數(shù)字地震臺(tái)網(wǎng)臺(tái)站、國(guó)家測(cè)震基本臺(tái)和國(guó)家地磁基本臺(tái)。該臺(tái)觀測(cè)歷史悠久,臺(tái)基條件好(圖2)。 本文選取恩施地震臺(tái)站2018年1月至2020年3月期間記錄到的連續(xù)波形數(shù)據(jù),通過(guò)計(jì)算功率譜密度PSD和地震噪聲概率密度函數(shù)PDF,并將2018年、2019年和2020年垂直向地震噪聲的加速度幅值進(jìn)行對(duì)比,分析新冠疫情防控前后恩施州地震背景噪聲水平的變化特征。

1 數(shù)據(jù)處理和方法

本文選取2018年1月至2020年3月恩施地震臺(tái)站記錄的連續(xù)波形數(shù)據(jù),參照 Peterson[11]、McNamara等[4]、McNamara等[12]和王芳等[13]的方法,不對(duì)干擾數(shù)據(jù)進(jìn)行篩選,將各臺(tái)站垂直分量的連續(xù)地震波形記錄去傾和去平均等預(yù)處理后,將其劃分為長(zhǎng)度是1 h(3 600 s)的數(shù)據(jù)段,相鄰數(shù)據(jù)段50%重合,估算功率譜密度和概率密度函數(shù)。

圖2 湖北恩施地震監(jiān)測(cè)中心站Fig.2 Enshi seismic station in Hubei Province

1.1 功率譜密度的估算

噪聲功率譜密度PSD是衡量地震臺(tái)站噪聲水平的常規(guī)方法,最為常用的PSD值估算方法是通過(guò)離散傅里葉變換,即對(duì)觀測(cè)數(shù)據(jù)序列在有限時(shí)間范圍內(nèi)通過(guò)快速傅里葉變換進(jìn)行計(jì)算[14-15]。

周期時(shí)間序列y(t)的有限范圍傅里葉變換可表示為:

大美新疆（柳京） .........................................................................................................................................11-56

(1)

式中:tr為時(shí)間序列段的長(zhǎng)度;f為頻率。

對(duì)離散頻率值fk,傅里葉變換可定義為:

(2)

式中:fk=k/(NΔt),其中,k=1,2,3,…,N-1;Δt為采樣間隔;N為截取時(shí)段的采樣點(diǎn)數(shù)。

功率譜密度PSD定義為:

(3)

將式(3)中速度功率譜密度PSD轉(zhuǎn)換成加速度PSD,計(jì)算公式為:

Pα,k(f)=(2πf)2Pk(f)

(4)

將每小時(shí)計(jì)算樣本劃分為13個(gè)時(shí)間段,為減少PSD的方差,每段數(shù)據(jù)之間在時(shí)間上50%重疊。通過(guò)上式(4)計(jì)算每時(shí)間段的功率譜密度值,每小時(shí)的PSD值由13個(gè)時(shí)間段的PSD值平均得到。

1.2 概率密度函數(shù)的計(jì)算

(5)

式中:NPTc是功率譜密度值落在某PSD取值范圍內(nèi)的記錄段個(gè)數(shù);NTc是中心周期Tc在整個(gè)功率譜密度值范圍內(nèi)估計(jì)值的總數(shù)。概率密度函數(shù)PDF反映的是在所選取的地震觀測(cè)數(shù)據(jù)數(shù)列中,對(duì)某個(gè)特定中心周期,某個(gè)特定功率值出現(xiàn)的概率。

通過(guò)此方法對(duì)周期-功率值進(jìn)行頻數(shù)統(tǒng)計(jì),可以繪制地震噪聲概率密度函數(shù)PDF分布圖,具體計(jì)算由IRIS提供的程序來(lái)實(shí)現(xiàn)。該圖可與Peterson[11]的全球高低噪聲模型進(jìn)行對(duì)比。

2 地震背景噪聲強(qiáng)度變化

本文基于湖北省恩施地震臺(tái)站記錄到的連續(xù)觀測(cè)數(shù)據(jù),系統(tǒng)計(jì)算了2018年1月至2020年3月間垂直向分量記錄的連續(xù)地震數(shù)據(jù)的功率譜密度PSD和概率密度函數(shù)PDF,具體結(jié)果如下。

圖3展示的是2018至2019年以及2020年防控期間的概率密度函數(shù)圖。從圖中可以看出,封城前后的地震噪聲功率譜在1 s及更長(zhǎng)周期上并沒(méi)有出現(xiàn)明顯變化。而封城后短于1 s周期地震噪聲功率譜出現(xiàn)了明顯降低,如在0.1 s左右,功率譜出現(xiàn)了10～15 dB下降。這些高頻地震噪聲的顯著變化,對(duì)應(yīng)了封城期間抑制的人類(lèi)活動(dòng)和相關(guān)工業(yè)設(shè)施的運(yùn)轉(zhuǎn)。

圖3 疫情防控前后恩施地震臺(tái)站垂直分量的PDF分布圖Fig.3 Vertical PDFs of Enshi seismic station before and after the prevention and control of COVID-19 pandemic

圖4是恩施地震臺(tái)垂直向地震噪聲功率譜在2018至2019年間及2020年恩施封城期間的日變規(guī)律圖?？梢钥闯?2018—2019年間恩施地區(qū)高頻(～0.1 s)地震噪聲在早晨6時(shí)為-125 dB,隨后持續(xù)升高,并在中午12時(shí)左右到達(dá)峰值-115 dB;隨后(因午休？)降低數(shù)dB,在16時(shí)再次達(dá)到峰值;在晚上20時(shí)左右地震噪聲功率逐漸降低。此變化規(guī)律與恩施作息時(shí)間一致。而在疫情防控期間,日間幅值顯著下降,高頻噪聲功率譜相比之前同期降低～10 dB。需要指出的是,夜間噪聲低值略有下降,1 s以上的長(zhǎng)周期地震噪聲并沒(méi)有出現(xiàn)明顯變化,這清晰地表明恩施地區(qū)最主要的高頻地震噪聲來(lái)源于人類(lèi)活動(dòng)。

圖4 恩施地震臺(tái)垂直向地震噪聲功率譜日變規(guī)律圖Fig.4 Daily variation law of vertical seismic noise power spectrum at Enshi seismic station

需要注意的是,2020年1月24日恩施開(kāi)始封城,日期與中國(guó)農(nóng)歷春節(jié)重疊。春節(jié)期間,地震高頻噪聲常因人類(lèi)活動(dòng)減弱而減低。為了比對(duì)和消除往年春節(jié)期間人類(lèi)活動(dòng)減弱的影響,本文將2018年、2019年和2020年垂直向地震噪聲的加速度幅值進(jìn)行對(duì)比(圖5)。在2018—2020年間春節(jié)前幾日的加速度幅值變化曲線(xiàn)高度一致,顯示了臺(tái)站所在區(qū)域高頻噪聲信號(hào)和地震儀器非常穩(wěn)定。在2020年1月24日(圖5的X軸上第754儒略天)0時(shí)起,恩施封城,地震高頻噪聲迅速降低。

圖5 2018年、2019年和2020年恩施地震臺(tái)記錄的垂直向加速度波形分鐘中值(1～10 Hz濾波)Fig.5 Minute mid-value of vertical acceleration waveform recorded by Enshi seismic station in 2018, 2019, and 2020(filtered at 1～10 Hz)

圖5為恩施地震臺(tái)2018年(黑點(diǎn))、2019年(藍(lán)點(diǎn))和2020年(紅點(diǎn))垂直向記錄的加速度波形(1～10 Hz濾波)的分鐘中值。圖中X軸上的第754 儒略天(Julian天)為農(nóng)歷除夕。將三年的觀測(cè)數(shù)據(jù)按中國(guó)農(nóng)歷除夕進(jìn)行了對(duì)齊。

2018年和2019年春節(jié)期間,地震高頻噪聲也出現(xiàn)了一定程度地下降,但其僅下降約20%,并且在春節(jié)7天假期結(jié)束后迅速恢復(fù)至節(jié)前水平,隨后逐漸變強(qiáng)趨勢(shì)(圖5)。然而,2020年新冠疫情防控期間,高頻噪聲加速度值一直維持在低幅值。

3 討論和總結(jié)

基于恩施地區(qū)2018年1月至2020年3月期間地震臺(tái)網(wǎng)監(jiān)測(cè)數(shù)據(jù),通過(guò)計(jì)算噪聲功率譜密度PSD和概率密度函數(shù)PDF,對(duì)恩施地區(qū)疫情前后地球背景噪聲進(jìn)行定量化分析,得出恩施地區(qū)地震高頻噪聲與人類(lèi)活動(dòng)密切相關(guān)。

通過(guò)觀察和比對(duì)恩施州疫情前后的概率密度函數(shù)PDF,發(fā)現(xiàn)疫情封控期間(2020年1月24日至3月25日),恩施地區(qū)采取嚴(yán)格的疫情管控措施,該地區(qū)高頻噪聲(～10 Hz)明顯下降,且整體水平下降了10～20 dB。然而,在周期大于1 s 的長(zhǎng)周期頻帶上,背景噪聲功率沒(méi)有出現(xiàn)明顯變化。在周期小于1 s的短周期頻段或高頻頻段,測(cè)振噪聲主要來(lái)源于人文噪聲,即人類(lèi)活動(dòng)引起的環(huán)境噪聲(如公路、鐵路、工廠等引起的噪聲)[2,16]。圖6顯示的是2019至2020年春節(jié)前后恩施州內(nèi)出行強(qiáng)度、遷出及遷入規(guī)模趨勢(shì)圖。由圖可以得出,與2019年相比,2020年恩施封城后州內(nèi)出行強(qiáng)度急劇下降,遷出和遷入恩施州的趨勢(shì)降幅更加顯著,出行強(qiáng)度與圖5噪聲數(shù)據(jù)趨勢(shì)基本一致。所以,疫情封鎖措施使得由人類(lèi)活動(dòng)振動(dòng)產(chǎn)生的高頻地震噪聲明顯下降,該結(jié)果與Lecocq等[5]、Xiao等[17]得出的結(jié)果一致。

圖6 2019及2020年春節(jié)前后恩施州內(nèi)出行強(qiáng)度、遷出及遷入規(guī)模趨勢(shì)(數(shù)據(jù)來(lái)源: 百度遷徙https://qianxi.baidu.com/)Fig.6 Travel intensity, emigration, and immigration in Enshi Prefecture during the Spring Festival in 2019 and 2020 (Data source: Baidu migration https://qianxi.baidu.com/)

基于恩施地震臺(tái)垂直向地震噪聲功率譜在2018—2019年間及2020年恩施封城期間的日變規(guī)律發(fā)現(xiàn),疫情前,在周期小于1 s頻帶內(nèi),背景噪聲具有明顯的日變化特征,從早晨6時(shí)開(kāi)始慢慢升高,中午12時(shí)達(dá)到峰值,隨后中午12—14時(shí)下降,而后逐漸上升在16時(shí)再次達(dá)到峰值,在晚上20時(shí)左右地震噪聲功率逐漸降低,該頻段噪聲日變化幅值變化與恩施生活作息時(shí)間一致;在疫情封控期間,高頻噪聲功率譜比之前同期降低～10 dB,且沒(méi)有之前的明顯的時(shí)段變化規(guī)律。周期大于1 s頻帶內(nèi),疫情前后沒(méi)有明顯變化。由此也表明,由于人類(lèi)活動(dòng)的減少,高頻段噪聲下降。

地震臺(tái)站記錄到的高頻噪聲變化與人類(lèi)活強(qiáng)度相關(guān),所以可以通過(guò)檢測(cè)震動(dòng)噪聲的強(qiáng)度變動(dòng)來(lái)監(jiān)測(cè)社會(huì)公共事件以及各級(jí)機(jī)構(gòu)對(duì)相關(guān)重大政策的響應(yīng)情況。新冠疫情防控為研究人為噪聲與人類(lèi)活動(dòng)關(guān)系提供了客觀依據(jù),由此從地震學(xué)視角為人類(lèi)活動(dòng)監(jiān)測(cè)、城市規(guī)劃、防災(zāi)減災(zāi),公共管理等提供參考。

致謝:本PDF計(jì)算所采用的IRIS程序來(lái)自IRIS網(wǎng)站,中國(guó)地質(zhì)大學(xué)(武漢)陳世驕研究生協(xié)助繪制了部分圖件,在此一并表示感謝。