2016年太原兩次有感地震前太原臺(tái)大地電場異常信號(hào)分析

李 艷，陳?？。瑥埩炼?/p>

(1.山西省地震局，山西 太原 030021；2.山西省地震局太原基準(zhǔn)地震臺(tái)，山西 太原 030025；3.太原大陸裂谷動(dòng)力學(xué)國家野外科學(xué)觀測研究站，山西 太原 030025)

0 引言

地電場是重要的地球物理場之一，根據(jù)不同場源，可將地電場劃分為大地電場和自然電場。大地電場是由外部的各種場源(如電離層、磁層、對(duì)流層等中的各種電流體系)在地球表面感應(yīng)產(chǎn)生的分布或較大地區(qū)的變化電場，具有廣域性；自然電場場源是指地球內(nèi)部各種因素形成的場，具有局部性。隨著“九五”“十五”項(xiàng)目建設(shè)，大地電場觀測得到普及，研究者也發(fā)現(xiàn)了地震前出現(xiàn)的異常電信號(hào)。統(tǒng)計(jì)發(fā)現(xiàn)，異常電信號(hào)一般出現(xiàn)在發(fā)震前數(shù)天至一個(gè)月左右[1-4]，最長為震前三個(gè)月以內(nèi)[4]，屬典型的短期異?；蚨膛R異常。

太原基準(zhǔn)地震臺(tái)大地電場(以下簡稱太原臺(tái)電場)建設(shè)于2006年10月，2008年1月投入正式運(yùn)行，已積累11 a年的觀測數(shù)據(jù)。2016年11月18日，距臺(tái)站北東方向5 km的小店區(qū)發(fā)生ML3.9地震，市區(qū)普遍有感;一個(gè)月后的12月18日，在臺(tái)站正南方向11 km的清徐發(fā)生ML4.6中等地震，市區(qū)震感明顯。根據(jù)多極距去噪音觀測原理，對(duì)2016年7月以后的觀測資料分析發(fā)現(xiàn)，10月至12月底共提取到5組異常。其中，10月有2組：1組為低頻突跳、1組呈低頻突跳并夾雜高頻擾動(dòng)；11月18日ML3.9地震前5 d又提取到2組，其特征與10月的2組一樣；12月18日ML4.6地震前10 d提取到1組高頻擾動(dòng)異常信號(hào)。5組異常距2次地震最短5 d、最長76 d。

1 觀測概況

太原臺(tái)電場屬于武家寨電磁臺(tái)觀測項(xiàng)目，其觀測場地位于太原市晉源區(qū)晉源鎮(zhèn)武家寨村，太原市西南20 km，電磁臺(tái)占地面積約3 000 m2，電場測區(qū)地形平坦、開闊。臺(tái)站地形地質(zhì)構(gòu)造上處于太原盆地中北部，交城大斷裂(走向NE)東側(cè)5 km，地勢(shì)開闊平坦，第四紀(jì)沉積厚度約250 m，其下部巖層巖性為沙質(zhì)泥巖、泥質(zhì)沙巖、砂礫巖，成巖較好。

太原臺(tái)電場觀測儀器型號(hào)為ZD9A-II，采用Pb-PbC12固體不極化電極,埋設(shè)在地表凍土層以下的潮濕土壤中，深度3 m，“雙L型”多極距布極方式(見第27頁圖1)。其中，O1A1為EW向長極距240 m，O1B1為NS向長極距240 m，A1B1為NW向長極距339 m，O2A2為EW向短極距160 m，O2B2為NS向短極距 160 m，A2B2為NW向短極距226 m。在無干擾、無地電暴影響等正常背景下，3個(gè)測向呈雙峰、雙谷弱日變化形態(tài)，兩峰幅度一強(qiáng)一弱，兩谷谷值接近，其中NS和NW向日變形態(tài)特征一致(見第27頁圖2)。

2 異常信號(hào)分析

該文異常信號(hào)識(shí)別以馬欽忠老師提出的多極距去噪音方法為依據(jù)[1-4]，該方法原理為：地震異常電信號(hào)一般來自測區(qū)數(shù)公里以外遠(yuǎn)源場，該場是均勻的或準(zhǔn)均勻的；工業(yè)游散電流、降雨、灌溉等噪音源一般具有近源特點(diǎn)，即噪音源位于測點(diǎn)不遠(yuǎn)的位置上,其在長短極距上產(chǎn)生的電場分量在測區(qū)局部不均勻性明顯，即近源場是非均勻的。異常電信號(hào)至少同時(shí)滿足三個(gè)條件：未被所有臺(tái)站記錄到，不是廣域性的大地電磁效應(yīng)；同時(shí)出現(xiàn)在同測向的長極距和短極距上，變化形態(tài)相同、幅度相近；長極距上記錄到的電位極性必須與平行于它的短極距上記錄到的電位極性相同。通過上述方法篩選確定與地震有關(guān)的異常電信號(hào)。

圖1 太原臺(tái)大地電場觀測布極示意圖Fig.1 Schematic diagram of geoelectric field observation distribution at Taiyuan station

圖2 太原臺(tái)電場2015年6月18日分鐘值曲線圖Fig.2 Minute value curve of Taiyuan station electric field on June 18, 2015

2.1 干擾信號(hào)

圖3是噪音干擾類不均勻性的例子，如2015年10月23日記錄到3組擾動(dòng)信號(hào)(虛線框標(biāo)識(shí))。其中，EW向長短極距在5:00-6:00、17:00-19:00、23:00-24:00三組時(shí)間段內(nèi)擾動(dòng)形態(tài)相同，變化幅度差別大；NW向在5:00-6:00第1組的長短極距記錄到信號(hào)形態(tài)相反(電位極性相反)，18:00-19:00第2組的短極距上記錄到幅度為11 mV/km的負(fù)極性信號(hào)(向下)，長極距則未記錄到類似信號(hào)；23:00-24:00第3組的長極距記錄到最大幅度為18 mV/km的負(fù)極性信號(hào)，短極距則無相似信號(hào)出現(xiàn)。根據(jù)前述判據(jù)，3組擾動(dòng)信號(hào)呈非均勻性特征，屬于近源場的噪音干擾，噪音源應(yīng)位于A1極附近。經(jīng)現(xiàn)場巡查，發(fā)現(xiàn)在A1極東南約400 m處有市政在開挖供熱管道和污水管道工程施工，期間的電焊施工可能是產(chǎn)生上述干擾的原因。

圖3 太原臺(tái)電場2015年10月23日分鐘值曲線圖Fig.3 Minute value curve of Taiyuan station electric field on October 23, 2015

2.2 震前信號(hào)

2016年11月18日，距電場測區(qū)北東方向5 km的太原小店區(qū)發(fā)生ML3.9地震；12月18日，距電場測區(qū)正南方向11 km的清徐發(fā)生ML4.6地震，如此近距離的2次地震事件為檢驗(yàn)太原臺(tái)電場監(jiān)測效能提供了條件。據(jù)文獻(xiàn)[4]所述，地震異常電信號(hào)最長出現(xiàn)在震前三個(gè)月內(nèi)，基于此，對(duì)2016年7月至12月太原臺(tái)電場6個(gè)測道的觀測數(shù)據(jù)進(jìn)行逐日分析。通過查詢觀測日志，排除各種已知干擾源(如漏電、降雨、農(nóng)田灌溉、地電暴等)，再依據(jù)馬欽忠老師提出的多極距去噪音觀測原理的三項(xiàng)必要條件，提取到5組異常信號(hào)。其中有兩組連續(xù)兩天出現(xiàn)，小店和清徐地震均發(fā)生在此異常信號(hào)之后。

第28頁圖4顯示2016年10月2日太原臺(tái)電場記錄到第1組異常，根據(jù)多極距去噪音原理方法，確定為異常電信號(hào)。異常形態(tài)為顯著的低頻突跳，各測道上的突跳形態(tài)相似、極性相同(見線框標(biāo)識(shí))，持續(xù)約2 h。其中，NS向擾動(dòng)幅度最大為15 mV/km，EW向和NW向約5m V/km,第1組異常距11月18日太原ML3.9地震46 d，距離12月18日清徐ML4.6地震76 d。

第28頁圖5顯示10月14日太原臺(tái)電場記錄到第2組異常，根據(jù)多極距去噪音原理方法，確定為異常電信號(hào)。異常特征以低頻突跳為主，同時(shí)夾雜有少許高頻擾動(dòng)信息，各方向擾動(dòng)形態(tài)相似、極性相同(見線框標(biāo)識(shí))，持續(xù)近一個(gè)半小時(shí)。NS向擾動(dòng)幅度最大約14 mV/km，NW向約8 mV/km,EW向約3 mV/km,第2組異常距11月18日太原ML3.9地震34 d，距12月18日清徐ML4.6地震64 d。

第28頁圖6顯示11月12日至13日，太原臺(tái)電場連續(xù)兩天記錄到3組異常信號(hào)，根據(jù)多極距去噪音原理方法研判，有大地電磁效應(yīng)(虛線框)，也有異常電信號(hào)(實(shí)線框)。其中，12日的異常為低頻突跳，六測道擾動(dòng)形態(tài)相似、極性相同，持續(xù)約2 h。突跳幅度特征與第1組類似，NS向突跳幅度最大為11 mV/km，EW向和NW向擾動(dòng)幅度為6 mV/km，此為第3組異常；第4組異常發(fā)生在13日，異常特征以低頻突跳為主，同時(shí)夾雜有高頻擾動(dòng)信息，各方向長短測道擾動(dòng)形態(tài)相似、極性相同，持續(xù)6小時(shí)左右。擾動(dòng)幅度特征與第2組類似，NS向擾動(dòng)幅度最大為17 mV/km，NW向約10 mV/km,EW向約6 mV/km,此組異常距11月18日太原ML3.9地震5 d，距12月18日清徐ML4.6地震35 d。

圖4 太原臺(tái)電場2016年10月2日分鐘值曲線圖Fig.4 Minute value curve of Taiyuan station electric field on October 2, 2016

圖5 太原臺(tái)電場2016年10月14日分鐘值曲線圖Fig.5 Minute value curve of Taiyuan station electric field on October 14, 2016

上述4組異常信號(hào)與隨后2次地震事件間隔在5 d至兩個(gè)半月不等，符合文獻(xiàn)[4]中最長3個(gè)月的異常統(tǒng)計(jì)。因2次地震事件震中位置不同，無法準(zhǔn)確判斷某組異常歸因于某個(gè)地震，僅從異常信號(hào)形態(tài)和幅度來看，1、3組有一定的相似性，即變化形態(tài)相近，異常信號(hào)呈負(fù)極性變化(向下)，NS向異常幅度最大，EW和NW兩測向變化幅度近似相等，前者是后者變化的2～3倍；2、4組異常存在相似性，即形態(tài)特征相近，各方向長短測道極性相同，異常幅度NS向最大，NW向次之，EW向最??；4組異常在三測向的擾動(dòng)幅度明顯大于日變形態(tài)幅度。

圖6 太原臺(tái)電場2016年11月12日至13日分鐘值曲線圖Fig.6 Minute value curve of Taiyuan station electric field from November 12 to 13, 2016

11月18日距電場北東方向5 km發(fā)生ML3.9地震，20天后(即12月8日)出現(xiàn)連續(xù)的高頻擾動(dòng)信號(hào)，持續(xù)約4 h(見圖7線框標(biāo)識(shí))。與前4組異常不同的是此次異常幅度不大，NS向擾動(dòng)幅度最大約7 mV/km，EW向和NW向分別約4 mV/km，均未超當(dāng)日最大變化幅度；擾動(dòng)信號(hào)以高頻形式出現(xiàn)，其中NS向和EW向高頻信號(hào)特征比NW向明顯。根據(jù)多極距去噪音觀測原理方法排除了噪音干擾的可能，確定為異常電信號(hào)(第5組)。圖中顯示，各測道長短極距異常形態(tài)和幅度呈高度一致性，10 d后的12月18日，距電場觀測場地正南11 km的清徐發(fā)生ML4.6地震。

圖7 太原臺(tái)電場2016年12月8日分鐘值曲線圖Fig.7 Minute value curve of Taiyuan station electric field on December 8, 2016

3 結(jié)論與討論

通過以上分析可知，太原臺(tái)電場在2016年11月18日ML3.9地震、12月18日ML4.6的2次地震前提取到5組顯著異常信號(hào)。其中，第1組異常時(shí)間距第一次地震約一個(gè)半月，距第二次地震兩個(gè)半月；第5組異常信號(hào)在第一次地震發(fā)生后20 d，在ML4.6地震前10 d。異常信號(hào)特征如下：

(1) 異常信號(hào)存在方向指示性。5組異常信號(hào)中，NS向異常擾動(dòng)幅度大于EW和NW向，說明來自NS方向的地震異常電信號(hào)強(qiáng)度優(yōu)于其他兩測向，2次震中的位置分別位于臺(tái)站的北東和正南，說明異常信號(hào)似乎有一定的方向指示性。

(2) 異常信號(hào)距震中非常近。太原臺(tái)電場測區(qū)距2次地震震中分別為5 km和11 km，查閱有關(guān)電場異常信號(hào)的震例研究文獻(xiàn)，發(fā)現(xiàn)電場信號(hào)源一般距震中幾十公里到幾百公里不等，位于震源區(qū)附近的異常震例文獻(xiàn)更是少見，震源區(qū)附近出現(xiàn)的異常電信號(hào)，可能在異常機(jī)理上需要更明確的物理解釋。

(3) 不同組異常信號(hào)圖像存在“相似性”特征。10月2日和11月12日出現(xiàn)的第1、3組異常為低頻信號(hào)，其變化幅度相近、變化極性相同，在圖像特征方面存在“相似性”；10月14日和11月13日記錄到的2、4組異常電信號(hào)(低頻+高頻擾動(dòng))，在幅度、極性等方面也存在“相似性”，兩者可能反映了不同地震事件的電場信號(hào)屬性；12月8日記錄的第5組異常(高頻擾動(dòng))，可能是18日清徐ML4.6地震前的異常反應(yīng)。

(4) 5組異常信號(hào)持續(xù)時(shí)間不等。第1、3組異常持續(xù)約2 h；第2組異常(11月12日出現(xiàn))持續(xù)時(shí)間最短，約一個(gè)半小時(shí)；第4組異常(11月13日)持續(xù)時(shí)間最長，約6 h；第5組持續(xù)約4 h。