李惠玲,張淑亮,李 麗,高云峰

(1.山西省地震局，山西 太原 030021；2.太原大陸裂谷動力學國家野外科學觀測研究站，山西 太原 030025)

0 引言

2010年河津M4.8地震是山西地球物理臺網(wǎng)數(shù)字化改造后發(fā)生在山西境內(nèi)的一次顯著地震。震前山西及鄰區(qū)多種地球物理場觀測記錄到豐富的異常信息。河津地震發(fā)生后，一些學者對這次地震前地球物理場異常特征進行研究并取得一些成果。張淑亮等通過對河津地震臨汾臺VP寬頻帶傾斜儀震前異常的研究，討論了野外實際觀測中的亞失穩(wěn)現(xiàn)象[1]；楊俊芳等通過對韓城臺金屬擺傾斜資料的分析，總結(jié)了汾渭地震帶山西及陜西地區(qū)1989年-2010年間的幾次顯著地震及汶川8.0級地震前的異常變化特征[2-3]；李艷、張紅秀等通過對臨汾地傾斜、地應變震前異常特征研究，討論了河津地震等一系列中強地震前捕捉到的顯著短臨異常信息[4-8]。上述研究對于了解河津地震地球物理場異常特征具有重要作用，但這些研究多是基于單臺站或單測項進行，對于多臺多測項異常的時空特征涉及較少。

該文擬在河津地震前震源及鄰近地區(qū)形變、流體、電磁學科大量觀測數(shù)據(jù)分析基礎上，系統(tǒng)梳理異常起始時間、異常形態(tài)、異常幅度等，分析震前異常的時、空演化特征，并應用相關數(shù)據(jù)處理方法分析多種地球物理場觀測數(shù)據(jù)的時頻特征。研究結(jié)果有可能為地震綜合預測與研究提供震例支撐。

1 地球物理場異常概況

根據(jù)2010年-2011年山西省年度地震趨勢研究報告以及張淑亮等人的研究成果(1)山西省地震局.山西省2010年度地震趨勢研究報告,2009.(2)山西省地震局.山西省2011年度地震趨勢研究報告,2010.[1-11]，按照中國震例編寫規(guī)范中研究范圍的要求(5級地震研究區(qū)為200 km)[12]，同時為避免遺漏敏感點[13]的異常，分析時將研究區(qū)擴大至距震中250 km的范圍。河津4.8級地震前，研究區(qū)內(nèi)共有22個地球物理場觀測臺站，17類地球物理場定點觀測項目(見第2頁圖1)，包括金屬擺、水管儀、伸縮儀、水平擺、VS垂直擺、VP垂直擺、體應變、水準、地電阻率、大地電場、地磁、水位、水溫、水氡、氣氡、水汞、氣汞等60個觀測手段。震前有19個測項存在異常，異常表現(xiàn)特征如第2頁表1及第3頁圖2至第4頁圖10所示。

可以看出，河津地震前地球物理場異常數(shù)量與張淑亮等人對山西地區(qū)中強地震研究結(jié)果[14](5級以上地震異常臺項數(shù)一般為15～17項，4.5～4.9級地震異常臺項數(shù)一般為10項左右)相比，異常數(shù)量是山西地區(qū)同等級別地震的2倍左右；異常在空間上主要分布在距離震中較近的山西南部至晉陜交界區(qū)，在空間上表現(xiàn)出不均勻性；在時間上以短臨異常為主，且具有明顯的階段性特征。

圖1 河津4.8級地震250 km范圍內(nèi)定點地球物理場觀測臺站分布圖Fig.1 Distribution map of fixed-point geophysical observation stations within 250 km of Hejin M4.8 earthquake

2 異常時空特征

由第4頁表2可知，2010年山西河津4.8級地震前異常持續(xù)時間大于6個月的中期異常3項，占異?？倲?shù)的15.79%；短期異常和臨震異常16項，占異?？倲?shù)的84.21%。異常出現(xiàn)時間最早的是夏縣水氡反年變高值異常，出現(xiàn)在震前8～9個月；其次是臨汾臺短水準、臨汾VP垂直擺，分別出現(xiàn)在震前7個月和5～6個月；16項短臨異常，出現(xiàn)的時間較為集中，主要在震前13天～4個月和2～6天兩個時段內(nèi)。異常在時間進程上具有階段性和準同步性特征，距離發(fā)震時間愈近，異常數(shù)量越多，震前20天就出現(xiàn)了11項異常表現(xiàn)出叢集性特征(見第4頁圖11)。

從河津4.8級地震多項地球物理測項異常空間演化圖11看出，震前的19項異常分布在距震中40～248 km范圍內(nèi)。距離最近的異常為韓城金屬擺和侯馬水管傾斜，分別距震中40和58 km；最遠為乾陵垂直擺，距震中248 km，異常分布范圍較廣。其中，距離震中100 km范圍內(nèi)的有11項，異常占比52%；101～200 km范圍內(nèi)異常有1項，異常占比4%；在201～250 km范圍內(nèi)有7項，異常占比41%(見表2)?？梢钥闯?，異常主要集中在0～100 km和200～250 km內(nèi)，空間上表現(xiàn)出不均勻性。其原因可能與觀測臺點空間分布和震中區(qū)構造的復雜性有一定的關系。由圖1看出，河津4.8級地震在250 km范圍內(nèi)，觀測臺點多數(shù)分布在與盆地走向一致的NE方向，在震中的東部與西部幾乎無觀測點，這種空間分布格局的不均勻性是造成異常空間不均勻性的主要原因；其次，由于震中及附近地區(qū)構造的復雜性，NE韓城斷裂、NEE向的峨嵋臺地北緣斷裂和孤山斷裂、NW向的西辛封斷裂和NWW向羅云山斷裂等在震中周圍密集交錯分布，在區(qū)域應力的作用下，因斷裂的不均一性在某些部位介質(zhì)結(jié)構發(fā)生了改變。河津4.8級地震波速比異常的分布特征也證實了這一推斷。震中以北的臨汾和安澤地震臺，在震前具有持續(xù)2～3 年的中期波速比低值異常，臨汾地震臺西南方位和隰縣地震臺正南方位的波速比時間曲線出現(xiàn)短期異常變化，震中距最近的夏縣臺在震前無明顯的波速異常[15]。

表1 河津4.8級地震前定點地球物理場觀測映震效能分析統(tǒng)計Table 1 Analysis and statistics of earthquake reflecting efficiency of fixed point geophysical field observation before Hejin M4.8 earthquake

圖2 河津4.8級地震夏縣臺水氡觀測曲線Fig.2 Water radon observation curve at Xiaxian station before Hejin M4.8 earthquake

圖3 河津4.8級地震前臨汾臺VP寬頻帶傾斜儀觀測曲線及局部放大圖Fig.3 VP broadband tiltmeter observation curveand partial enlarged figure at Linfen station before Hejin M4.8 earthquake

圖4 河津4.8級地震前韓城臺金屬擺NS分量日均值觀測曲線Fig.4 Daily mean observation curve of NS component of metal pendulum at Hancheng station before Hejin M4.8 earthquake

圖5 河津4.8級地震前西安臺VS垂直擺傾斜觀測曲線Fig.5 Tilt observation curve of VS vertical pendulum at Xi'an station before Hejin M4.8 earthquake

圖6 河津4.8級地震前西安臺水管傾斜觀測曲線Fig.6 Observation curve of water pipe inclination at Xi'an station before Hejin M4.8 earthquake

圖7 河津4.8級地震前西安臺洞體應變觀測曲線Fig.7 Observation curve of body strain in Xi'an station cave before Hejin M4.8 earthquake

圖8 河津4.8級地震前夏縣臺氣汞觀測曲線Fig.8 Gas mercury observation curve at Xiaxian station before Hejin M4.8 earthquake

綜上所述，河津4.8級地震異常在時間上具有階段性、準同步性和叢集性，在空間上表現(xiàn)出范圍廣和不均勻性，且異常圍繞震源區(qū)出現(xiàn)聚集或擴散現(xiàn)象。中期階段異常從震源區(qū)向外圍地區(qū)擴展，進入短期階段后，異常有一個快速的擴展過程，即從外圍地區(qū)向震源區(qū)快速收縮推進，之后進入以高、低頻波動為主的臨震異常階段(見圖11)。

圖9 河津4.8級地震前西安臺鉆孔水位觀測曲線Fig.9 Water level observation curve of boreholes at Xi'an station before Hejin M4.8 earthquake

圖10 河津4.8級地震前乾陵臺VS垂直擺觀測曲線Fig.10 Observation curve of VS vertical pendulum at Qianling station before Hejin M4.8 earthquake

表2 河津4.8級地震前定點地球物理場測項及異常數(shù)量統(tǒng)計Table 2 Statistics of fixed-point geophysical field measurement items and anomaly number before Hejin M4.8 earthquake

圖11 河津4.8級地震前地球物理場異常時空演化及進程圖Fig.11 The spatio-temporal evolution and process diagram of geophysical field anomalies before Hejin M4.8 earthquake

3 異常時頻特征

異常的時頻特征研究是定量識別異常的一種有效方法，前述的異常識別基本以定性判別為主，為了快速準確地判別異常，采用小波多分辨率時頻分析方法[16]對震前異常進行處理。這是由于這種信號分解的能力可將各種交織在一起的不同頻率組成的混合信號分解成不同頻帶的子信號，在研究非穩(wěn)態(tài)信號、識別異常信息方面具有優(yōu)勢[17-18]。

從河津地震前各異常測項小波時頻分析結(jié)果看出，震前3個多月始，有17個測項存在不同程度的異常(限于篇幅，小波分解僅給出異常最先出現(xiàn)的細節(jié)層)。如圖12a所示,震中距100 km范圍內(nèi)的異常主要集中在小波分解細節(jié)部分第5層(32～64 h)或第6層(64～128 h)；距震中100～200 km范圍內(nèi)的異常集中在第7層，周期為128～256 h(見第6頁圖12b)；距震中200～250 km范圍內(nèi)的異常主要集中在第6層(64～128 h)或第7層(128～256 h),如圖12c所示。

圖13中，橫軸為出現(xiàn)異常的細節(jié)層，縱軸為處理結(jié)果在此細節(jié)層上有異常的測項數(shù)與總異常項的百分比。可以看出，距震中100 km范圍內(nèi)，第4層出現(xiàn)異常的比率為0.08，第5層為0.69，第6層為0.23，第7層無異常出現(xiàn)(見圖13a)；距震中101～200 km范圍內(nèi)出現(xiàn)的異常僅1項(見圖13b),出現(xiàn)在第7層；距震中201～250km范圍內(nèi)異常資料處理結(jié)果，第4層和第5層無異常出現(xiàn)，第6層和第7層出現(xiàn)異常的比率分別為0.64和0.36(見圖13c)。結(jié)果表明，隨著臺站距震中的由近到遠，其頻率越來越低，周期亦越來越長。各測項異常信息特征如第7頁表3所示。

圖12 小波分解結(jié)果Fig.12 Wavelet decomposition results

圖13 小波分解異常頻段統(tǒng)計圖Fig.13 Statistical chart of abnormal frequency band by wavelet decomposition

表3 河津4.8級地震前定點地球物理場觀測小波時頻分析結(jié)果統(tǒng)計Table 3 Statistics of wavelet time-frequency analysis results of fixed-point geophysical field observations before Hejin M4.8 earthquake

4 結(jié)論

通過對山西河津4.8級地震前250 km范圍內(nèi)地球物理場測項異常信息分析，得出到如下認識：

(1) 在時間上，以短臨異常為主，主要集中在震前13天～4個月及2～6天內(nèi)，異常在時間上具有階段性和準同步性特征，并且距離發(fā)震時間愈近異常數(shù)量越多，表現(xiàn)出叢集性特征。

(2) 在空間上，異常分布范圍較廣(40～248 km)，主要集中在距離震中200 km范圍內(nèi)，中期異常由震中內(nèi)向外圍擴展，短臨階段表現(xiàn)為震中與外圍交替出現(xiàn)。異常在空間的分布隨時間而變化，具有不均勻性。

(3) 震前異常測項時頻分析結(jié)果顯示，異常包含的周期成分較復雜，隨著臺站距離震中的遠近，異常頻段表現(xiàn)出一定的規(guī)律性。離震中越近的測項，異常高頻信息相對豐富；反之，異常低頻信息相對豐富。這一規(guī)律對未來地震震中判定具有一定的參考意義。

致謝：國家地震臺網(wǎng)中心和陜西省地震局為該研究提供了數(shù)據(jù)支持，在此深表感謝！