唐 淋,祁國亮,張欏訢
(1. 四川省地震局,四川 成都 610041;2. 成都市應急管理局,四川 成都 610042)
中國是世界上地震災害最嚴重的國家之一,地震災害嚴重威脅我國經濟、社會可持續(xù)發(fā)展。當破壞性地震發(fā)生后,及時評估地震災情,是高效、有序、快速地進行抗震救災應急指揮決策的關鍵[1]。地震烈度速報系統(tǒng)可通過觀測儀器直接計算地震烈度,無需到現(xiàn)場調查即可快速獲知地震影響程度,為人員傷亡估計、經濟損失評估、應急救援指揮和工程搶險修復決策等提供依據,具有顯著的減災實效[2,3]。日本、美國、臺灣等國家和地區(qū)以及成都、蘭州、唐山、西安等國內大中型城市均建設了地震烈度速報系統(tǒng)或城市地震烈度速報臺網[4-13]。由中國地震局實施、覆蓋全國的國家地震烈度速報與預警工程項目也已開始正式建設,預計建成后將大幅提高我國地震災害應急處置能力[14]。
觀測數據是地震烈度速報臺網發(fā)揮地震烈度速報功能的基礎,數據質量優(yōu)劣將直接決定產出成果的可靠性和準確性,影響人員傷亡估計、經濟損失評估、地震應急救援指揮和工程搶險修復決策的科學性,因此,保障觀測數據質量是地震烈度速報臺網運行維護的核心。
但地震烈度速報系統(tǒng)基于強震動觀測系統(tǒng),其大規(guī)模建設晚于測震和前兆等其它地震監(jiān)測系統(tǒng),臺站數量和臺網規(guī)模長期發(fā)展緩慢,對觀測資料質量的關注程度也嚴重滯后。當前,強震觀測資料質量評價側重于記錄獲取率、記錄質量、記錄應用價值、記錄報送等一般指標[15],沒有充分考慮地震烈度速報對記錄獲取實時性、強震動記錄的可用性和臺網整體運行率的特殊要求,不適用于地震烈度速報臺網觀測資料的質量評價。成都市地震烈度速報臺網經過多年運行,近年來逐漸出現(xiàn)臺站運行率降低、數據傳輸實時性不足、資料質量下降等問題,導致個別臺站的速報地震動參數、部分地震事件的地震烈度圖與現(xiàn)場調查結果和實際情況有較大出入,烈度速報系統(tǒng)的穩(wěn)定性、準確性和可靠性有所降低。
因此,明確地震烈度速報臺網觀測資料質量內涵,分析成都市地震烈度速報臺網觀測資料存在的質量問題,提出合理可行的改進措施,不僅有利于提高成都市地震烈度速報系統(tǒng)觀測資料質量和使用價值,也可為國家地震烈度速報與預警工程建成后的高效運維提供借鑒。
地震烈度速報臺網觀測資料既包括原始強震動加速度記錄,也包括對原始記錄分析、處理后生成的地震動參數和地震烈度圖等成果資料。因此,地震烈度速報臺網觀測資料質量由原始記錄質量和成果資料質量兩部分構成。其中,原始記錄質量包括記錄的完整性、可靠性和可用性,成果資料質量則主要指準確性。原始記錄的質量是地震烈度速報臺網觀測資料質量的基礎,成果資料的質量則是核心。
1.1.1 原始記錄的完整性
原始記錄的完整性是包括地震烈度速報臺網在內的各類地震監(jiān)測網數據質量評價最為基礎的一項指標。它體現(xiàn)了地震烈度速報臺網實際產出的數據總量在應產出數據總量中的比重,主要反映數據采集和網絡通訊設備的運行狀態(tài)。將完整率作為衡量原始記錄完整性的主要參數,定義為:
完整率=(總時間-總丟數時間)/總時間×100%
1.1.2 原始記錄的可靠性
原始記錄的可靠性用于評價觀測數據是否真實客觀反映了臺站所在地的實際振動行為。加速度計、數據采集器、電源避雷器等設備的機械故障、電路故障和一些干擾都可能造成臺站輸出波形異常,使得原始記錄不能真實反映臺站所在地的實際振動行為,降低原始記錄的可靠性。以事件波形正常率作為衡量原始記錄可靠性的主要指標,定義某一次地震事件波形正常率為:
事件波形正常率=正常波形臺站個數/臺站總數×100%
1.1.3 原始記錄的可用性
原始記錄可靠并不意味著可用,當加速度計被安裝在較高的樓層、臺站周邊存在的大型擾動源時,容易造成單個臺站與周邊其它臺站的地震動參數差距懸殊,當守時精度不高和數據傳輸延時過大時,容易導致單個臺站的地震動參數計算結果有誤,此時的原始記錄雖然真實客觀反映了臺站所在地的實際振動行為,是可靠的,但無法使用,不能參與地震動參數和烈度平面圖的繪制。原始記錄的可用性決定觀測數據能否用于地震動參數計算和成果圖件繪制等成果資料的產出,臺址條件、臺基環(huán)境、守時狀態(tài)和數據傳輸效率是其主要影響因素。以事件波形可用率作為衡量原始記錄可用性的主要參數,定義某一次地震事件波形可用率為:
事件波形可用率=可用波形臺站個數/臺站總數×100%
1.1.4 成果資料的準確性
原始資料具備較好的完整性、可靠性和可用性,是保證成果資料較高準確性的基礎。但同時,地震動參數計算方法、臺站密度和布局方式以及離散點的網格化方法等因素對地震動參數和地震烈度圖等成果資料準確性有重要影響。將成果資料的準確性定義為評價地震烈度速報系統(tǒng)產出成果的準確度和真實性、衡量成果資料標識的地震破壞范圍和程度是否與現(xiàn)場調查結果一致的重要指標。
決定地震烈度速報臺網觀測資料質量的因素復雜多樣,總結起來主要有以下幾點。
1.2.1 設備運行狀態(tài)
地震烈度速報臺站包含加速度傳感器、數據采集器、傳輸設備和電源系統(tǒng)等設備,分別擔負數據采集、傳輸、供電等任務,其中任何一項設備出現(xiàn)故障,都會導致臺網中心接收不到數據或接收到的數據異常,影響原始記錄的完整性、可靠性和可用性,進而對成果資料的準確性造成影響。臺站設備運行狀態(tài)是決定地震烈度速報臺網觀測資料質量的首要因素。
1.2.2 數據傳輸效率
不同于傳統(tǒng)的強震動臺網,地震烈度速報臺網對成果的產出時限有較高要求,一般震后5~10 分鐘內即要產出自動地震烈度速報圖,臺站監(jiān)測數據必須要實現(xiàn)實時或準實時在線傳輸,通信速率、信道誤碼率和系統(tǒng)時延等技術指標均要滿足相應的要求,否則將直接影響地震烈度速報成果的產出時效和質量,故數據傳輸效率對地震烈度速報臺網觀測資料質量有重要影響。
1.2.3 臺站場地條件
場地條件直接影響地震動記錄質量,理想的場地條件是保證強震動記錄真實反映自由地表地面運動的基本前提。根據相關行業(yè)標準,遠離大型建筑、局部地形起伏變化不大和軟硬均勻的自由場地是較為理想的強震動觀測臺址,它可以避免周圍環(huán)境建筑和結構振動、地形、土層效應的影響,獲得真實的強震動原始記錄[16]。而條狀突出的山嘴、高聳孤立的山丘、陡坡、河岸、盆地、河谷等復雜地形、附近的高大建筑物以及軟硬不均勻的場地則極可能引起地震動放大、與附近地面運動不一致的情況。
1.2.4 環(huán)境噪聲水平
環(huán)境振動是由隨機振源激發(fā)并經場地不同性質的巖土層界面多次反射和折射后傳播到場地地面的振動。雖然它對強震動臺站觀測記錄的影響不如測震臺站那么明顯,但幅值較大的一些環(huán)境振動,因震動源強度、傳播介質特性等因素影響,對強震動記錄質量仍有較大的影響,尤其是大型的馬達、泵站、發(fā)電機、塔柱狀結構、重型車輛通路、大型管道、機場和鐵路設施等干擾源,會顯著提高臺站的環(huán)境地噪聲水平,降低觀測記錄質量[16]。
1.2.5 儀器烈度的計算方法
由于地震烈度定義的模糊性、地震破壞及地震動記錄本身的復雜性,宏觀地震烈度與實際儀器記錄的地震動參數間沒有明確、直接的函數關系,但研究發(fā)現(xiàn),地震烈度與地震動參數之間仍存在某種統(tǒng)計對應關系[17-19],根據強震動記錄計算的儀器地震烈度仍可較為直觀的反映地震造成的潛在破壞和危險程度,選擇合適的儀器地震烈度計算方法對提高地震烈度速報結果的準確性具有重要作用。
1.2.6 臺站密度及其布局方式
不論是完全基于儀器觀測值的地震烈度速報,還是在實測地震烈度基礎上根據地震動衰減模型約束插值后的地震烈度速報,臺站密度及其布局方式都直接影響地震烈度圖的產出時效和質量。
在觀測臺站密集且分布相對均勻的地區(qū),將實測記錄進行場地響應校正后,進行常規(guī)的空間網格化處理,即可較好的勾畫地震能量的釋放過程,準確得到地震動影響場的實際情況,精確繪制出地震烈度等值線圖。而在臺站分布稀疏、布局不均勻的地區(qū),由于中強地震的地震動場在空間上高度離散,較少的臺站不能如實地反映地震動場的復雜性,需要根據地震動參數衰減關系由鄰近臺站記錄數據通過插值方法來獲得空白區(qū)的地震動參數值,從而彌補實測數據的欠缺[20-21]。
1.2.7 分析軟件使用的地震烈度衰減模型和公式
當臺站稀疏且布局不均勻時,地震動參數衰減模型、公式的準確性和適用性在很大程度上直接影響地震烈度圖等成果資料的可靠性。為此,要根據地震震級和破裂范圍大小,從點源模型、橢圓模型和線源模型中選擇合適的模型作為地震動參數衰減模型,在此基礎上,根據區(qū)域地震活動特點和歷史地震資料,選擇適應當地實際的衰減公式,從而為成果資料準確產出奠定基礎。
1.2.8 離散點的網格化方法
由各個臺站觀測地震動參數計算得到的地震烈度離散數據,必須經過網格化處理后,才能繪制出地震烈度等值線圖。根據不同的標準,可以對網格化方法進行分類,而不同的網格化方法有不同的效果,針對目標地區(qū)的實際情況,通過分析和研究不同的網格化方法的適應范圍及參數設置,選取一種最適合本地區(qū)的模型,對于客觀、正確地描述地震烈度情況有重要的意義。
成都市于2009—2010 年建設了由1 個臺網中心、67 個子臺組成的城市地震烈度速報臺網。臺站設備包括力平衡式加速度傳感器、GPS 授時系統(tǒng)、24 位數據采集器、不間斷電源等,采用CDMA 無線傳輸方式向臺網中心發(fā)送觀測數據。臺網中心配置有數字強震數據處理系統(tǒng),可對監(jiān)測數據實現(xiàn)實時自動化處理。當成都市及周邊發(fā)生強烈地震時,該臺網能在震后5 分鐘左右生成全市地震動參數和地震烈度分布圖。
2.2.1 原始記錄不完整
(1)事件記錄部分缺失。主要表現(xiàn)為強震動記錄的丟頭(圖1a)、斷尾和中間缺失等現(xiàn)象,其主要原因多為數據采集器或傳輸故障。這種情況將對后續(xù)加速度峰值、持續(xù)時間和反應譜計算產生嚴重影響,尤其是加速度和速度峰值計算錯誤會進一步導致儀器地震烈度計算錯誤,從而降低地震烈度圖的可靠性。
(2)事件記錄完全缺失。主要表現(xiàn)是強震動記錄完全缺失地震事件發(fā)生期間的數據(圖1b)。主要原因可能是信號傳輸故障或強震計、數據采集器等儀器故障。出現(xiàn)事件記錄完全缺失的臺站無法計算地震動參數和地震烈度,不能參與地震烈度圖繪制,降低了地震烈度圖的精度。
2.2.2 部分臺站的原始記錄不可靠
(1)記錄波形有明顯奇異。如個別分量數據異常、階梯狀基線漂移、尖刺(圖1c)等,這些異常主要是由強震計故障造成。記錄波形明顯異常容易造成地震動參數和地震烈度計算錯誤,導致地震烈度圖可靠性降低。
(2)臺站沒有記錄到地震事件。地震發(fā)生期間,附近其它臺站均對當次地震事件有清晰的記錄,獲得較大的峰值加速度值,但個別臺站卻沒有記錄到相應的事件(圖1d)。出現(xiàn)此類現(xiàn)象的原因主要有兩個方面,一是數據采集器故障,記錄顯示的是數采的自噪聲;二是臺站所在地噪聲水平較高,強震動記錄信噪比較低,信號被噪聲所淹沒。同樣,這些異常將直接導致地震動參數計算錯誤,由此獲得的儀器地震烈度偏低,影響地震烈度圖的可靠性。
2.2.3 個別臺站的原始記錄不可用
部分臺站的強震計因安裝樓層較高(3 層以上)或沒有穩(wěn)定的固定在擺墩上等場地原因,與附近其它臺站相比,會出現(xiàn)地震烈度明顯偏高或偏低的現(xiàn)象。如2013 年4 月20 日8 時02分蘆山M7.0 級地震,新津金華臺與雙流黃龍溪臺相距僅4km,峰值加速度分別為181.69gal 和66.08gal,儀器烈度相差近2 度(圖1e、圖1f)。當這些臺站參與地震烈度圖繪制時,會出現(xiàn)一些以臺站為中心的孤立小圈,即常說的“牛眼”現(xiàn)象(圖2)。


圖1 臺站原始記錄圖Fig.1 Original record of the seismic stations

圖2 2013 年4 月20 日8 時02 分蘆山M7.0 級地震速報烈度圖Fig.2 Earthquake intensity rapid map of the M7.0 Lushan earthquake at 8:02 am on April 20,2013
2.2.4 烈度圖精度有待提高
成都市地震烈度速報臺網使用的烈度速報處理軟件能對原始加速度記錄和積分后的速度記錄進行自動濾波處理,使用地震動速度峰值(PGV)計算地震烈度,具備等距離法、克里格法等地震動參數網格化方法選項,可利用震源、斷層和烈度衰減模型等約束插值過程,還可根據需要分別生成地震動加速度峰值(PGA)分布圖、地震動速度峰值(PGV)分布圖和地震烈度圖等各類圖件,總體上功能強大、計算高效、圖件美觀,由其建立的地震烈度影響場基本可以作為震災應急救援指揮決策的依據。
但該軟件的地震烈度計算方法與即將頒布的《儀器地震烈度計算》 行業(yè)標準(以下簡稱行業(yè)標準)存在一定的差異,據其繪制的地震烈度圖精度低于行業(yè)標準地震烈度圖精度。以三方向地震動參數為例,現(xiàn)軟件計算儀器地震烈度的公式為:

行業(yè)標準計算儀器地震烈度的公式為:

其中:

從式中可以看出,現(xiàn)軟件僅使用PGV 計算儀器地震烈度,而行業(yè)標準在Ⅵ度及以上區(qū)域使用PGV,Ⅵ度以下區(qū)域則聯(lián)合使用PGV 和PGA。由于PGV 對小地震的中高頻成分較為敏感,在低烈度區(qū)聯(lián)合使用PGV 和PGA 可以有效避免因高頻成分影響而出現(xiàn)不合實際的大烈度異常。圖3 為蘆山M7.0 級地震行業(yè)標準儀器烈度與兩種儀器烈度差的交會圖,可以看出,現(xiàn)有的儀器烈度值較行業(yè)標準儀器烈度值整體偏大,且烈度值越低,差異越明顯。

圖3 行業(yè)標準儀器烈度與兩種儀器烈度差的交會圖Fig.3 The crossplot of the industry standard instrument intensity and the difference between the two instruments
圖4 與圖2(b)使用了相同的成圖參數和方法,但儀器地震烈度計算方法不同,分別由式(1)和式(2)計算得到。兩圖對比發(fā)現(xiàn),圖4 中的Ⅵ度、Ⅶ度區(qū)面積較圖2(b)大,Ⅴ度區(qū)較圖2(b)小,且缺少Ⅳ度區(qū)。與中國地震局發(fā)布的人工調查地震烈度圖(圖5)相比,圖2(b)與之更接近,尤其是成都境內Ⅶ度區(qū)和Ⅵ度區(qū)的面積、分布范圍兩者基本一致,說明根據行業(yè)標準計算的儀器地震烈度與現(xiàn)場調查烈度更為接近,成都市地震烈度速報臺網使用的儀器地震烈度計算公式有待改進。

圖4 2013 年4 月20 日8 時02 分蘆山M7.0 級地震速報烈度圖Fig.4 Earthquake intensity rapid map of the M7.0 Lushan earthquake at 8:02 am on April 20,2013

圖5 2013 年4 月20 日8 時02 分蘆山M7.0 級人工調查地震烈度圖(據中國地震局)Fig.5 Manual earthquake intensity map of the M7.0 Lushan earthquake at 8:02 am on April 20,2013
通常,各級地震部門只公布M5.0 級以上地震的人工調查地震烈度圖,且標注的最低烈度為Ⅵ度。成都市地震烈度速報臺網建成運行以來,四川省內發(fā)生的M5.0 級以上地震中,對成都影響較大、烈度達到Ⅵ度以上的主要是2013年4 月20 日蘆山M7.0 級地震,因上一節(jié)已對該次地震烈度速報成果資料的準確性進行了詳細分析,故以下僅對兩次典型地震事件的原始觀測記錄質量進行分析(圖6)。

圖6 典型地震事件震中位置圖Fig.6 Epicenter distribution map of typical seismic events
2.3.1 2013 年4 月20 日蘆山M7.0 級地震
(1)原始記錄完整率。接入臺網的66 個臺站中,彭州銀廠溝臺等5 個臺站數據整體缺失,大邑西嶺臺等4 個臺站數據部分缺失,其中大邑西嶺臺數據缺失較多,臺網原始記錄完整率為90.19%。
(2)事件波形正常率。在原始記錄基本完整的60 個臺站中,金堂淮口臺沒有真實記錄到此次地震,該臺站原始記錄不可靠,整個臺網事件波形正常率為89.40%。
(3)事件波形可用率。在事件波形可靠的59 個臺站中,成華區(qū)建設局臺等3 個臺站因強震計安裝于較高樓層或沒有與地面很好的固定,導致其與附近其它臺站的烈度值相比明顯偏高,無法直接用于繪制地震烈度圖,整個臺網的事件波形可用率為84.85%。
2.3.2 2017 年11 月10 日綿竹M4.2 級地震
(1)原始記錄完整率。接入臺網的66 個臺站中,植物園臺等13 個臺站數據整體缺失,大邑霧山臺等2 個臺站數據部分缺失,其中邛崍城區(qū)臺數據缺失較多,臺網原始記錄完整率為78.16%。
(2)事件波形正常率。在原始記錄基本完整的52 個臺站中,金堂趙鎮(zhèn)臺等2 個臺站的各分量均沒有真實記錄到此次地震,大邑霧山臺等2 個臺站的部分分量沒有真實記錄到此次地震,這些臺站原始記錄不可靠,整個臺網事件波形正常率為72.73%。
(3)事件波形可用率。在事件波形可靠的48 個臺站中,新津金華臺因強震計安裝于較高樓層,導致其與附近其它臺站的峰值加速度值相比明顯偏高,整個臺網的事件波形可用率為71.21%。
兩次典型地震事件原始記錄資料分析表明,成都市地震烈度速報臺網建成初期原始記錄質量較高,但隨著臺網的持續(xù)運行,儀器設備故障逐漸增多,原始記錄質量呈現(xiàn)下降趨勢(圖7)。

圖7 典型地震事件原始記錄質量柱狀圖Fig.7 Histogram of the original recording quality of typical seismic events
針對成都市地震烈度速報臺網原始記錄和成果資料存在的質量問題,需要從決定地震烈度速報臺網觀測資料質量的各種因素出發(fā),采取有針對性的措施,切實提高觀測資料質量。
3.1.1 全面完善規(guī)章制度
建立地震烈度速報臺網運維管理辦法,完善值班、巡查、應急、數據分析處理等各項制度,進一步界定維護管理部門、人員的職能職責,健全臺網(站)設備故障維修機制,明確軟(硬)件運行狀態(tài)檢查、儀器維修、強震計標定檢測的流程和周期,規(guī)范臺站巡檢、工作日志填寫、數據整理分析、檔案資料整理的程序,確保地震烈度速報臺網安全、穩(wěn)定和可靠運行。
3.1.2 定期檢查設備狀態(tài)
借助現(xiàn)場巡檢、遠程抽檢、定期標定等手段,重點對儀器參數設置、加速度傳感器零位電壓、定時設備狀態(tài)、電池電壓、軟件功能等內容進行核驗,全面及時掌握加速度傳感器、數據采集器、供電系統(tǒng)、避雷系統(tǒng)、傳輸設備等臺站硬件設備及監(jiān)控軟件、分析處理軟件和網絡設備等臺網中心軟硬件的運行狀態(tài)。
3.1.3 及時維修故障設備
日常做好車輛、人員、工具等方面的準備,加強技術、通訊線路和備件的儲備,當臺站儀器出現(xiàn)故障時,能及時派人前往現(xiàn)場維修,保證較高的修復率,不能現(xiàn)場修復時,也能通過更換備件及時恢復臺站運行;當臺網中心服務器、軟件和通訊線路故障時,能立即切換備份服務器,啟用備用通訊線路,從而保證臺網的連續(xù)觀測。
3.2.1 果斷停測或搬遷個別臺站
對于個別場地條件不理想、環(huán)境噪聲水平較高、周邊干擾較大的臺站,無論如何精細管理和維修維護,其強震動記錄都無法真實反映自由地表地面運動和用于地震烈度圖繪制,任其長期無效運行意義不大,應該果斷停測或勘選新的臺址,進行異地搬遷,使其發(fā)揮應有的作用。
3.2.2 充分利用其它行業(yè)的資源和服務優(yōu)勢
與地震行業(yè)類似,氣象、通信、油氣開發(fā)等行業(yè)在全國各地分布有大量的氣象觀測站、通信基站、廢棄井場等各類規(guī)模不等的站點,具有優(yōu)越的站址、通訊、電力和人力資源服務優(yōu)勢。應該摒棄以往獨立建設和運維的定式思維,通過與這些行業(yè)創(chuàng)新合作,依托現(xiàn)有資源,實現(xiàn)地震烈度速報臺網的高效建設和運維,進而提高觀測數據的綜合質量。
3.2.3 著力提高監(jiān)控智能化水平
通過安裝臺站環(huán)境監(jiān)控單元、網絡視頻攝像頭、軟件運行狀態(tài)檢測程序和數據質量監(jiān)控軟件,對臺站設備、臺網中心軟件的運行狀態(tài)和臺網數據質量進行實時監(jiān)控,一旦發(fā)現(xiàn)異常可及時報警,提醒臺網運行人員處理故障,進一步減小軟硬件故障對臺網數據質量的影響。
3.3.1 改進儀器地震烈度計算方法
即將頒布的儀器地震烈度計算行業(yè)標準所提出的指標建立在充分的強震動數據、震后震害數據和各類科學研究成果基礎之上,它繼承、拓展并銜接《中國地震烈度表》 (GB/T 17742-2017)、《中國地震動參數區(qū)劃圖》(GB/T 18306-2015)等相關標準,以易于計算和理解的加速度峰值和速度峰值作為表征儀器地震烈度的地震動參數,具有較強的科學性、繼承性、一致性和可操作性,是以后儀器地震烈度計算的推薦方法,因此,成都市地震烈度速報臺網也應將其作為儀器地震烈度計算的首選方法。
3.3.2 優(yōu)化地震動衰減模型和公式
近年來,我國許多學者對西部地區(qū),尤其是四川地區(qū)的地震烈度衰減模型做過深入研究,并給出了一些適用于各研究區(qū)的衰減公式[22-25]。與現(xiàn)軟件中已有的東北經驗衰減公式和ShakeMap公式相比,這些衰減模型更適用于成都市地震烈度速報臺網,能更好的反映出真實的地震動衰減特征。因此,在繪制地震烈度圖使用地震動衰減模型時,應采用上述更加符合實情的地震動衰減關系替換系統(tǒng)中原有的地震動衰減。