戴 苗 李軍輝 劉 堅 李 明 蔣玲霞
1)湖北省地震局,武漢 430071
2)安徽省地震局,合肥 230031
皖西南地區中強地震前后地磁資料變化特征*1
戴 苗1)李軍輝2)劉 堅1)李 明1)蔣玲霞1)
1)湖北省地震局,武漢 430071
2)安徽省地震局,合肥 230031
計算并分析了2011年1月19日安慶M4.8地震和2011年9月10日瑞昌-陽新M4.6地震前后震中周邊主要地磁臺站地磁轉換函數參數的變化特征,總結了安慶地震前華東片主要地磁臺站地磁資料垂直分量的變化特征。結果表明:在兩次中強地震前,武漢、涇縣、金寨三個地磁臺的轉換函數的部分參數出現了不同程度的異常變化;安慶地震前,中國大陸東部地磁低點時間出現位移,地磁日變形態亦發生了明顯的畸變。
轉換函數;低點位移;中強地震;地磁;皖西南地區
目前,從地磁場的變化中提取震磁異常常見的方法,通常是在掌握地磁場正常變化規律的前提下,根據地磁場時空變化規律,利用多臺站、多測點的比較分析,排除地磁場本身的正常背景及外空場的干擾等非震因素的影響,尋找偏離正常變化的成分,并通過物理和數學的分析,篩選出有用的震磁信號,如在地磁場長期變化規律中常用的提取異常的簡單差值法[1]、復雜差值法[2]、時空參考場法[3]、空間相關法和加權差分法[4];在變化的磁場中提取地磁前兆信息的幅相法[5]、地磁垂直分量日變形態[6]、地磁加卸載響應比法[7]、低點位移法[8,9]、譜分析方法;利用地磁場的短周期變化的地磁轉換函數法[10]和地磁諧波振幅比法[11]等。
2011年1月19日,安徽安慶發生了 M4.8地震,同年,9月10日,江西瑞昌與湖北陽新交界發生M4.6地震,兩個地震時間間隔不到8個月,震級相近,震中距相差190 km。本文利用此兩次地震震中周邊主要地磁臺站絕對地磁觀測資料,采用地磁轉換函數、地磁低點位移、垂直分量日變形態三種方法探討地震前后地震震中周邊地磁臺地磁場變化特征,探索地磁場變化與中強地震的孕育和發生的相關性。
地磁轉換函數法可了解和監視震源及鄰區的地下電性結構,如果在地震發生前捕捉到這種變化,就可進行地震預報[4]。早在上個世紀90年代,汪江田等[12]就有過利用轉換函數法對我國上海及鄰近地區中強地震進行預報的嘗試,并確認該方法對中強地震中期預報的有效性。2008年,李琪等[10]分析了2006年7月4日河北文安地震前后,距離震中分別為52 km和298 km的靜海地磁臺和泰安地磁臺不同周期的轉換函數參數隨時間的變化過程,發現靜海地磁臺轉換函數參數的數據異常點會隨著時間的向后推移從短周期向長周期遷移。近年來,地磁轉換函數在我國不同區域、不同省份區域地磁臺網均有不同程度的應用研究,積累了部分地磁轉換函數預報地震的寶貴經驗和實例[13,14]。
地磁低點位移是指地磁垂直分量日變極小值出現時間因受到外來因素的影響偏離正常時間的現象,研究結果表明,全國較大范圍的地磁低點位移異常與我國中強地震活動有一定的關系[1]。目前研究人員利用地磁日變低點位移的兩條以上分界線交匯來研究其與地震之間的關系,以此縮小預測范圍[8]。由于應力的積累,孕震介質的結構等性質發生變化而產生的各種前兆也呈多樣性,地磁低點位移異常在我國不同地區有著不同的表現形式,反應不同地區、不同強度的地震信息。
變化磁場內場強弱既與外場有關,又與地下物質電導率的分布、大小有關。強震前孕震體環境由于應力增強引起介質電導率等電磁性質的改變,可觀測到與地震有關的局部或區域性磁異常。國內多次震例顯示,中強地震發生前,區域地磁場垂直分量日變化會出現的相位和幅度變化。
地磁短周期變化的垂直分量ΔZ、水平分量ΔH和磁偏角ΔD存在如下關系:
式中的A和B為轉換函數。轉換函數A和B的時空變化特征為地下電性構造的函數,因而可用于研究局部和部分地區地殼和上地幔電性結構及其變化,轉換函數A和 B是地點(λ,φ)、深度(h)、頻率(ω)和時間(t)的函數。由于每一個地磁短周期變化都可以對應一個方程(如式(1)),于是可以利用多個相同類型的地磁短周期變化建立一個矩陣方程,用最小二乘法得出轉換函數A和B、方位角αp以及總方差σz等。
本方法計算選用湖北、安徽、湖南等震中相鄰省份7套FHD質子矢量磁力儀數字化分鐘值數據,資料選用時間段為2009年1月—2012年5月。數據計算前先采用粗差處理技術對數據進行高通濾波,去除分鐘值中數據突跳等不可用數據,然后,計算預處理后地磁各要素分鐘值虛實富氏譜,在此基礎之上,計算并選取不同周期的轉換函數參數,最后,采用滑動平均方法依次去除轉換函數中的月以及年變趨勢。安慶地震前后,地磁轉換函數出現異常的臺站主要有湖北九峰地磁臺、安徽涇縣地磁臺和安徽金寨地磁臺,研究區域內異常地磁臺站及顯著地震事件震中分布如圖1所示,三個地磁臺站涇縣臺、金寨臺、九峰臺三個臺站距離安徽安慶地震震中依次為120、160、150 km,距離湖北陽新地震震中依次為310、230、120 km。經臺站所屬臺網人員核實,異常出現時間段觀測儀器的運行率較高,數據資料穩定可靠。

圖1 區域內異常地磁臺站及顯著地震事件震中分布Fig.1 Explanation of the map of epicenter distribution and geomagnetic observatory in the studied area
在安徽安慶4.8級地震前,2010年3月初開始,安徽涇縣臺轉換函數A不同周期曲線同步出現異常低值,異常時間持續6月余,10月開始,不同周期異常曲線逐步恢復,12月底恢復到正常水平(圖2)。在此期間,湖北九峰臺,除了周期為10分鐘的時間序列曲線外,轉換函數A各周期的時間曲線均出現了高值異常,其中周期為50分鐘的時間序列曲線出現的異常持續時間段最長。安徽金寨臺轉換函數A在震前5~6個月不同周期也同步“成片”出現高值異常。隨后,2011年1月19日發生安慶M4.8地震。2011年3月中旬開始,安徽涇縣臺轉換函數不同周期曲線同步再次出現異常高值,高值超過兩倍均方差,高值持續兩月余,異?;謴驼K胶笕齻€月,2011年9月10日發生江西瑞昌和湖北陽新交界M4.6地震。江西瑞昌和湖北陽新交界地震發生前一個月內,安徽金寨臺部分短周期段轉換函數A出現高值異常。安徽安慶地震、江西瑞昌和湖北陽新交界地震前,安徽涇縣臺轉換函數A雖均出現異常但曲線異常形態的表現形式卻不一致(前者為低值異常,后者為高值異常),可能與這兩次地震發震構造差異性[13]有一定關聯。
“低點位移”是指一個大區域各臺的日變低點時間明顯地與另一個大區域各臺的低點時間不同,而每個大區域內部低點時間又相對一致,兩個區域之間突變分界線的低點時間相差2小時以上,這種少見的磁場突變分界線與強震發生有較好的相關性。根據震例統計,地震的發震日期一般為異常日后第27天(或41天)前后4天,即:
低點位移日期+27(或41)天=發震日期

圖2 九峰(JF)、涇縣(JX)、金寨(JZ)地磁臺不同周期(T)轉換函數A時間序列曲線Fig.2 Variation curves of different period transfer function A in Jiufeng,Jingxian and Jinzhai Station
通過對安慶Ms 4.8級地震前中國大陸東部地區42個臺站的地磁Z分量的幅度和相位分析,結果顯示安慶地震前,華東地區出現地磁Z分量日變化幅度的極小值時間在空間上出現地磁低點位移現象,地磁低點時間的分界線沿九峰、金寨、新沂、安丘、大山等臺站延伸,整體走向與郯廬斷裂帶走向基本一致,呈現NNE走向,分界線兩側低點時間相差超過2小時,安慶地震發震日期為異常發生后的29(27±4)天,震中位于低點位移分界線附近,如圖3所示。

圖3 2010年12月21日地磁低點位移分布Fig.3 Low point displacement image on Dec.21,2010
安慶M4.8地震前29天(2010年12月21日),沿郯廬斷裂帶分布及附近的信陽、金寨、蒙城、郯城、陵陽、泰安等六個臺站地磁垂直分量除了日變幅度和相位變化顯著外,地磁日變形態還發生了明顯的畸變異常,出現了明顯的地磁日變“雙低點”現象(圖4)。
分析2011年1月19日安慶M4.8地震和2011年9月10日瑞昌-陽新M4.6地震前后地震周邊主要地磁臺站地磁轉換函數時間序列曲線,以及安慶地震前華東片主要地磁臺站地磁資料垂直分量變化特征,可得結論如下:
1)兩次顯著地震事件前后部分地磁臺站地磁轉換函數參數存在正負異常現象,但有的臺站反應較為明顯,出現“成片”現象,即幾個相鄰的周期同時出現異常,這與過去的震例基本一致。有的臺站幾乎沒有反應或者反應不夠顯著。這主要跟臺站所處的地理位置以及地震發震構造差異性有關。
2)依據各臺站地磁轉換函數參數值變化曲線,轉換函數參數異常為漸變的趨勢性異常,持續時間較長,而且異常消失后往往需要持續一段時間后發震,間隔的時間不一,該異常具有中期異常特征,對臨震異常特征識別不好區分和把握。
3)安慶M4.8震前29天,中國大陸東部地磁低點時間出現了明顯的位移,地磁日變形態還發生了明顯的畸變異常,即出現了明顯的地磁日變“雙低點”現象。
湖北及周邊地區顯著地震少,地磁震例不多,僅從本文兩個震例還不足以得出確定性的結論,還需在以后的研究中繼續總結和探索。

圖4 2010年12月21日中國大陸東部地區地磁臺站垂直分量日變曲線Fig.4 Daily curves of Z component in East China on Dec.21,2010
1 丁鑒海,盧振業,黃雪香.地震地磁學[M].地震出版社,1994.
2 任熙憲,祁貴仲,詹志佳.唐山地震前后北京地區地磁場總強度的變化[J].地震學報,1984,6(3):271-286.
3 姚法章.地磁總強度資料在地震預報中的應用[J].西北地震學報,1984,6(3):44-48.
4 林云芳,等.地磁方法在地震預報中的應用[J].地震地磁觀測與研究,1999,20(6):35-44.
5 蔣延林,王建宇,趙衛紅.鹽城臺-杭州臺地磁幅相法臨震異常標志體系[J].地震學刊,2000,20(1):45-48.
6 馮志生,等.地磁垂直分量日變幅逐日比及其與地震關系的探討[J].華南地震,2001,21(2):20-27.
7 曾小蘋,等.地球磁場對太陽風的加卸載響應與地震[J].地震地磁觀測與研究,1996,17(1):49-52.
8 丁鑒海,余素榮,肖武軍.地磁“低點位移”現象與昆侖山口西8.1級地震[J].西北地震學報,2003,25(1):16-21.
9 吳小平.2003.利用地磁低點位移預測發震時間[J].四川地震,2003,1:33 -36.
10 李琪,王建軍,楊冬梅.文安地震前后的轉換函數參數變化分析[J].地震地磁觀測與研究,2008,29(3):35-40.
11 馮志生,等.中強震前地磁諧波振幅比的趨勢性變化特征步研究[J].西北地震學報,2004,26(1):50-56.
12 汪江田,顧學章,龔耀.從巖石的應變、破裂與電性的變化關系探索對地震的預報[J].地震學刊,2000,20(3):1-8.
13 鄭在壯,沈瑞童.地磁短周期轉換函數在地震預報中的應用[J].地震地磁觀測與研究,2010,31(3):13-17.
14 周瓊,曹志磊,夏仕安.安徽地區地磁轉換函數變化與地震的關系[J].地震地磁觀測與研究,2011,32(5):55-59.
CHARACTERISTIC ANALYSIS OF GEOMAGNETIC DATA BEFORE AND AFTER MODERATELY STRONG EARTHQUAKES IN SOUTHWEST AREA OF ANHUI PROVINCE
Dai Miao1),Li Junhui2),Liu Jian1),Li Ming1)and Jiang Lingxia1)
1)Earthquake Administration of Hubei Province,Wuhan 430071
2)Earthquake Administration of Anhui Province,Hefei230031
The variation characteristic of transfer function parameter of geomagnetic data is analyzed which recorded by observatories in the time before and after M4.8 Anqing earthquake in January 19,2011 and M4.6 Ruichang-Yangxin earthquake in September 10,2011.The variation of geomagnetic field in the Z component before Anqing M4.8 earthquake in east China is also summarized.The results show that some transfer function parameters of Wuhan observatory,Jingxian observatory and Jinzhai observatory appeared different degrees of abnormalities before and after the two moderately strong earthquakes.The low-point time of the geomagnetic vertical component in east China appears anomalous features before the Anqing earthquake.Daily variation of geomagnetic field arises obvious distortion.
transfer functions;low-point time of the geomagnetic vertical component;moderately strong earthquake;geomagnetism;southwest area of Anhui province
P318.5
A
1671-5942(2013)Supp.(Ⅰ)-0143-05
2012-09-03
武漢市科技計劃項目(2013060501010166)
戴苗,男,1982年生,工程師,理學碩士,主要從事地球地磁學理論、方法和應用研究.E-mail:theocean123@163.com