蔣延林 袁桂平 李鴻宇 唐廷梅 趙衛(wèi)紅

1）江蘇省地震局高郵地震臺，江蘇省高郵市南郊地震臺路2號 225601

2）江蘇省地震局，南京 210014

3）山東省地震局，濟南 250014

0 引言

在地震孕育過程中，隨著孕震區(qū)應(yīng)力的增加，孕震區(qū)介質(zhì)電阻率一般會發(fā)生下降變化。馮志生等（2009）指出，基于電磁感應(yīng)理論的各種計算方法可用于探測孕震區(qū)深部電阻率變化信息，但對于提取與地震孕育過程有關(guān)的異常變化信息這一目標(biāo)而言，在目前對地震孕育過程中深部電阻率異常變化的機理并不十分清楚、深部電阻率異常變化所引起的變化磁場的異常很微弱的情況下，所采用的計算方法對實際地球介質(zhì)的假定和約束應(yīng)越少越好，而地磁諧波振幅比則是能滿足該條件的較好的方法。馮志生等（2004）初步研究表明，地磁諧波振幅比異常變化特征與（直流）地電阻率的變化特征類似，主要表現(xiàn)為下降-轉(zhuǎn)折-恢復(fù)上升的變化過程，地震基本發(fā)生在趨勢性下降后的轉(zhuǎn)折期、恢復(fù)上升期或恢復(fù)后的初期階段。異常持續(xù)時間與震級間的關(guān)系符合地電阻率的異常持續(xù)時間與震級間的關(guān)系，即時間越長，震級越高，但同震級地震的地磁諧波振幅比異常持續(xù)時間比地電阻率異常持續(xù)時間長。馮志生等（2013）進一步研究發(fā)現(xiàn)，地震一般發(fā)生在異?；謴?fù)上升期出現(xiàn)SN向與EW向不同步或長短周期不同步現(xiàn)象的臺站附近300km范圍內(nèi)以及僅有短周期出現(xiàn)異常且異常持續(xù)時間僅1～2年的臺站附近300km范圍內(nèi)，但此情況下震級一般僅為5級左右。

在以往一些研究中，研究區(qū)域內(nèi)、時段內(nèi)地震較多，地震與異常的一一對應(yīng)關(guān)系不明顯，此情形對提取異常特征和判據(jù)可能會產(chǎn)生干擾，而與此不同的是，孤立地震事件則有利于提取地震前后的異常特征。2012年7月20日江蘇省高郵市與寶應(yīng)市交界發(fā)生4.9級地震，該地震周邊200km范圍內(nèi)地磁臺站分布較密，其中高郵臺、鹽城臺和淮安臺的震中距分別僅為37、57、67km，同時該地震發(fā)生前后附近地區(qū)也未發(fā)生其他中等以上地震，此情形為研究地震地磁異常時空特征提供了便利。本文初步研究了該地震周圍地磁諧波振幅比異常特征并對其進行了定性解釋。

1 基本原理、資料處理與結(jié)果

地磁諧波振幅比 YZHx（NS）和 YZHy（EW）的定義為

其中，Z（ω）、Hx（ω）、Hy（ω）分別為地磁場的垂直分量、水平分量 SN向和 EW向譜；ω為圓頻率。對于隨時間周期變化的不均勻場源，在地球介質(zhì)為均勻各向同性的平面導(dǎo)體的條件下，則有

式中，μ為磁導(dǎo)率；σ為電導(dǎo)率；ω為圓頻率；λ為常數(shù)，是變化磁場的波數(shù)。在這種情況下，YZHx（NS）和 YZHy（EW）與地下介質(zhì)的電阻率 ρ呈正比，當(dāng) YZHx（NS）和 YZHy（EW）下降時，可能對應(yīng)了地下介質(zhì)的電阻率ρ的下降。

FHD質(zhì)子矢量磁力儀采用質(zhì)子磁力儀與線圈組合觀測原理（錢家棟等，2010），每分鐘觀測1組地磁總場F、水平分量（水平分量全矢量）H和磁偏角D，其中F和H為絕對值，D為相對值，F(xiàn)、H、D的觀測精度分別為0.3nT、0.6nT、0.1′。本次計算分析采用的是地震周圍臺站FHD質(zhì)子矢量磁力儀分鐘值三分量觀測資料，地震震中與臺站分布見圖1，各臺站FHD磁力儀開始觀測時間見表1，其中無錫臺資料質(zhì)量一直不佳，海安臺2009年資料不穩(wěn)定，凌陽臺2010年資料不穩(wěn)定，崇明臺一直受地鐵干擾，故以上4個臺的資料均未參與計算分析。

資料處理及計算過程采用馮志生編制的軟件①馮志生，2012，地磁諧波振幅比地磁轉(zhuǎn)換函數(shù)專用計算程序，具體計算步驟為：①將每日F、H和D資料轉(zhuǎn)換為Z、Hx和Hy；②去噪聲、線性去傾、加漢寧窗；③按周期等間隔計算5～65m in的每分鐘的三分量富氏譜 Z（ω）、Hx（ω）和 Hy（ω）；④按每10m in一個頻帶計算譜的頻帶均值；⑤計算各頻帶（10、20、30、40、50、60m in）的地磁諧波振幅比；⑥計算年滑動平均，消除年變化。各臺站計算結(jié)果如圖2所示。

2 異常特征

圖2（a）高郵臺震中距37km，各周期的地磁諧波振幅比約自2010年初開始出現(xiàn)下降，轉(zhuǎn)折時間大部分出現(xiàn)在2010年底，但除YZHx（NS）的10、20、30min出現(xiàn)完整的下降-轉(zhuǎn)折-恢復(fù)回升變化，并于震后的2013年初恢復(fù)到原有水平（圖2（a）中陰影部分），其它基本沒有出現(xiàn)恢復(fù)回升變化，即恢復(fù)階段不同步。另外，YZHx（NS）50m in未出現(xiàn)下降-轉(zhuǎn)折-恢復(fù)回升變化。

圖1 高郵-寶應(yīng)4.9級地震震中與地磁臺站分布

表1 地磁臺站觀測開始時間

圖2（b）鹽城臺震中距57km，各周期的地磁諧波振幅比下降出現(xiàn)時間前后跨度大，約從2009年底～2010年底，轉(zhuǎn)折時間大部分出現(xiàn)在2011年底，比高郵臺晚約1年，恢復(fù)時間大部分出現(xiàn)在震后的2012年底～2013年初，與高郵臺恢復(fù)時間大體一致，但恢復(fù)水平或不完全或過頭（圖2（b）中陰影部分）。另，YZHx（NS）10min一直到2011年中才出現(xiàn)小幅下降變化，而YZHy（EW）20min無下降變化。

圖2（c）淮安臺震中距 67km，YZHy（EW）10、20、30、40min未出現(xiàn)下降-轉(zhuǎn)折-恢復(fù)回升變化或變化不明顯，其他周期從2010年初開始出現(xiàn)下降變化，2011年先后發(fā)生轉(zhuǎn)折，恢復(fù)時間大部分出現(xiàn)在震后的2012年底～2013年初，除 YZHx（NS）10、20min外，其他周期都未能恢復(fù)到原有水平（圖2（c）中陰影部分）。

圖2（d）連云港臺震中距178km，YZHx（NS）各周期都出現(xiàn)了下降-轉(zhuǎn)折-恢復(fù)回升變化且形態(tài)完整，下降時間基本出現(xiàn)在2009年下半年，2012年先后發(fā)生轉(zhuǎn)折，恢復(fù)時間出現(xiàn)在震后的2012年底～2013年。YZHy（EW）各周期都出現(xiàn)了下降和轉(zhuǎn)折，且與 YZHx（NS）同步，但未出現(xiàn)恢復(fù)回升變化或恢復(fù)回升變化不明顯。

圖2 地磁諧波振幅比變化與地震

圖2（e）新沂臺震中距187km，各周期地磁諧波振幅比都出現(xiàn)了下降-轉(zhuǎn)折-恢復(fù)回升變化且形態(tài)完整，下降時間基本出現(xiàn)在2009年下半年，2011年先后發(fā)生轉(zhuǎn)折，恢復(fù)時間出現(xiàn)在震后的2012年底～2013年。

圖2（f）泰安臺震中距419km，各周期地磁諧波振幅比都出現(xiàn)了下降-轉(zhuǎn)折-恢復(fù)回升變化且形態(tài)完整，恢復(fù)回升同步性較好。

3 定性解釋

上述地磁諧波振幅比異常變化特征可概括為：基本變化形態(tài)為下降-轉(zhuǎn)折-恢復(fù)上升，其中，下降時間為2009～2010年，轉(zhuǎn)折時間為2011年，恢復(fù)上升時間為2012～2013年，高郵-寶應(yīng)4.9級地震發(fā)生在恢復(fù)上升期。距該地震震中較近的高郵、鹽城和淮安等臺站異?；謴?fù)期不同步現(xiàn)象明顯，距震中較遠的新沂臺和泰安臺異常形態(tài)完整，恢復(fù)期沒有出現(xiàn)不同步現(xiàn)象。

龔紹京等（2001）應(yīng)用三維有限差分法對長方形高導(dǎo)體的地磁轉(zhuǎn)換函數(shù)分布特征進行了數(shù)值模擬計算，根據(jù)計算結(jié)果的空間分布立體圖，地磁轉(zhuǎn)換函數(shù)A和B在高導(dǎo)體中央的變化小，極值主要分布在高導(dǎo)體4個犄角和邊界附近，且A與B在高導(dǎo)體4個犄角和邊界附近的變化分布明顯不同。由于地磁轉(zhuǎn)換函數(shù)A和B分別與地磁諧波振幅比YZHx和YZHy對應(yīng)，因此該模擬計算結(jié)果很好地解釋了地震附近臺站地磁諧波振幅比YZHx和YZHy的不同步現(xiàn)象，即當(dāng)孕震區(qū)出現(xiàn)高導(dǎo)異常后，位于高導(dǎo)異常體邊界附近臺站的地磁諧波振幅比YZHx和YZHy會出現(xiàn)不同步現(xiàn)象，故利用地磁諧波振幅比YZHx和YZHy出現(xiàn)不同步現(xiàn)象可以初步界定高導(dǎo)異常的邊界位置。

4 結(jié)語

本文通過孤立地震事件研究，進一步驗證了地磁諧波振幅比的異常時空特征，即地磁諧波振幅比異常特征表現(xiàn)出與地電阻率類似的下降-轉(zhuǎn)折-恢復(fù)上升的異常變化過程，地震一般發(fā)生在異常轉(zhuǎn)折或恢復(fù)上升的過程中，距震中較近的臺站會出現(xiàn)不同步現(xiàn)象——SN向與EW向的不同步以及長短周期的不同步。