流動地磁監(jiān)測數(shù)據(jù)處理流程

陳斌+袁潔浩+王粲+徐如剛+王振東

摘要：闡述了中國地震局流動地磁團隊現(xiàn)在采用的數(shù)據(jù)處理流程，包括最近臺站替代法的日變通化和自然正交分量模型的長期變化改正。并采用該流程對野外地磁測量數(shù)據(jù)進行了日變通化和長期變化改正，得到巖石圈磁場年度變化數(shù)值，殘差分析表明誤差不大于10 nT，一般在5 nT之內(nèi)。

關(guān)鍵詞：流動地磁監(jiān)測；數(shù)據(jù)處理流程；日變通化；長期變化改正；誤差分析

中圖分類號：P318文獻標識碼：A文章編號：1000-0666（2017）03-0335-05

0引言

地震的孕育及發(fā)生過程中，地下能量的積累與釋放過程必然會引起巖石圈介質(zhì)性質(zhì)，如溫度、應(yīng)力等狀態(tài)的改變。這一性質(zhì)的改變將通過各種震磁效應(yīng)引起巖石磁化率的改變（黃平章等，1990；郝錦綺等，1992，1993；侯登錄等，2003；李海燕，張世紅，2005），從而在地面產(chǎn)生能觀測到的巖石圈磁場異常變化。因此通過監(jiān)測地面的巖石圈磁場異常變化研究地震活動性是地磁學科主要研究方向之一。我國的地震學者通過監(jiān)測地磁臺站的連續(xù)觀測數(shù)據(jù)，發(fā)現(xiàn)了某些可能的地震地磁異?，F(xiàn)象，取得了如“低點位移”“轉(zhuǎn)換函數(shù)”等理論（馮志生等，2000，2004；劉哲，2002；史勇軍等，2002；王亶文，2003；陳伯舫，2003；丁鑒海等，2004；王世昌等，2004；張淑蘭，江在森，2004；吳小平等，2004；李樹華，2004；李琪等，2005；張繼紅，喬惠珍，2005；蔣延林等，2005；楊濤等，2005），但由于地磁臺站的空間間隔較大，這一類研究對空間的辨識率不高。流動地磁研究者也使用空間高分辨率的流動地磁數(shù)據(jù)進行了地震地磁相關(guān)性的分析（解用明等，2001，2002，2005；詹景立，黃健，2004），但僅重點分析以震中為中心的某一段時間-空間區(qū)域的變化，缺乏全時間-空間尺度的辨析，且不同地區(qū)、不同時間的地震活動對應(yīng)的地磁變化形態(tài)并不一致，這可歸因于其中的某些變化是磁場的正常變化。自2010年初以來，中國地震局流動地磁技術(shù)團隊在大華北、南北地震帶和南北天山3個區(qū)域持續(xù)開展了年尺度的矢量地磁復(fù)測，發(fā)現(xiàn)了一些與地震相關(guān)地磁異?，F(xiàn)象，如圍繞地震震中出現(xiàn)高梯度分布或者磁偏角零變線等現(xiàn)象（顧春雷等，2010，2012；倪喆等，2014a，b，c；蘇樹鵬等，2016）。本文主要闡述了流動地磁技術(shù)團隊采用的的數(shù)據(jù)處理過程，并分析數(shù)據(jù)處理過程中誤差的可能數(shù)值范圍。

1監(jiān)測目標

按照來源的不同，地球磁場可分為內(nèi)源磁場和外源磁場。內(nèi)源磁場可細分為源于地球液體外核磁流體動力學過程的主磁場和源于地殼及上地幔物質(zhì)磁性的巖石圈磁場。外源場可細分為源于固體地球之外的各種電流體系產(chǎn)生的磁場和固體地球內(nèi)部產(chǎn)生的感應(yīng)磁場。

由于巖石圈磁場包含了地殼和上地幔物質(zhì)狀態(tài)和結(jié)構(gòu)的豐富信息，因此通過對巖石圈磁場的研究可以推進對地球演化、地質(zhì)構(gòu)造演變及成礦研究等相關(guān)領(lǐng)域的基礎(chǔ)研究。從地震監(jiān)測預(yù)報目的來看，發(fā)生在地表60 km范圍內(nèi)的淺源地震往往比同級別的深源地震會帶來更大破壞（傅淑芳等，1980；許忠淮，2014），所以地磁場的地震監(jiān)測目標重點應(yīng)為對應(yīng)淺源地震深度的巖石圈磁場。巖石磁學的實驗室研究表明：巖石的狀態(tài)，如溫度、應(yīng)力狀態(tài)、巖石顆粒度及鐵磁性物質(zhì)含量等發(fā)生改變時，其磁性會發(fā)生相應(yīng)的改變，即產(chǎn)生如壓磁效應(yīng)、膨脹磁效應(yīng)、感應(yīng)磁效應(yīng)、熱磁效應(yīng)等震磁效應(yīng)。地震的孕育及發(fā)生過程中，地下能量的積累與釋放過程必然會引起巖石圈介質(zhì)性質(zhì)狀態(tài)的改變，并產(chǎn)生能觀測到的磁場異常變化。因此通過監(jiān)測巖石圈磁場的空間分布及時間變化的異常特征來研究地震活動性是我們的目標。

合理、準確地描述巖石圈磁場的空間分布及時間變化需要盡可能精準地從觀測數(shù)據(jù)中消除各種外源磁場以及地球主磁場的影響。因此，中國地震局制定了觀測技術(shù)規(guī)范中國地震局監(jiān)測預(yù)報司.2015.關(guān)于印發(fā)《流動地磁測量基本技術(shù)要求（試行）》的通知.以及日變通化、長期變化改正的數(shù)據(jù)處理流程，以便得到盡可能精確的巖石圈磁場數(shù)據(jù)。

2日變通化

本文闡述了中國地震局流動地磁團隊現(xiàn)有的數(shù)據(jù)處理流程，該流程依據(jù)地磁學的相關(guān)知識，利用臺站連續(xù)觀測數(shù)據(jù)，建立了獨立的地磁學日通化和長期變化模型，形成了穩(wěn)定可靠的數(shù)據(jù)處理流程，得到了高精度的巖石圈磁場年度變化數(shù)據(jù)。

選用測點最近地磁臺日變通化和自然正交模型長期變化改正的數(shù)據(jù)處理方法過程中需要有諸多特別需要注意的地方。該處理流程帶來的整體誤差小于10 nT，大部分數(shù)據(jù)誤差保持在5 nT之內(nèi)。但該流程仍然有諸多可以改進的地方，如日變通化處理過程中運用的三點近似處理方法和長期變化改正過程中對數(shù)據(jù)的強烈依賴性都是今后值得仔細思考和改進之處。

本文是中國地震局流動地磁團隊集體工作成果的闡述，所用流動地磁測量數(shù)據(jù)為全團隊的工作結(jié)果，在此向全體團隊成員表示感謝。日變通化和長期變化改正所使用的地磁臺站數(shù)據(jù)來自國家地磁臺網(wǎng)中心，感謝臺網(wǎng)中心和全國地磁臺站工作人員對流動地磁團隊的長期大力支持。中國地震局流動地磁技術(shù)團隊希望各位審者和讀者給我們提出寶貴的意見和建議，以促進我們團隊的健康進步，促進流動地磁監(jiān)測乃至地磁學科的科學研究水平的提高，促進我國的防震減災(zāi)事業(yè)長遠穩(wěn)固發(fā)展，在此再次感謝。

參考文獻：

陳伯舫.2003.日本鹿屋臺地磁轉(zhuǎn)換函數(shù)的變化[J].華南地震，23（1）：8-12.

丁鑒海，余素榮，肖武軍.2004.地磁異?，F(xiàn)象及其在地震預(yù)報中的應(yīng)用[J].石油地球物理勘探，39（增刊1）：27-30.

馮志生，梅衛(wèi)萍，張秀霞，等.2004.中強震前地磁諧波振幅比的趨勢性變化特征初步研究[J].西北地震學報，26（1）：50-56.

馮志生，王建宇，梅衛(wèi)萍，等.2000.江蘇地區(qū)地磁Z21測值空間相關(guān)異常及其標志體系[J].地震地磁觀測與研究，21（1）：44-49.endprint

傅淑芳，劉寶成，李文藝.1980.地震學教程[M].北京：地震出版社.

顧春雷，張毅，顧左文，等.2012.華北地震區(qū)巖石圈磁異常場零值線與中強震震中分布關(guān)系[J].西北地震學報，25（2）：174-179.

顧春雷，張毅，徐如剛，等.2010.地震前后巖石圈磁場變化特征分析[J].地球物理學進展，25（2）：472-477.

郝錦綺，黃平章，周建國.1992.巖石流變磁效應(yīng)——震磁效應(yīng)的一種可能機理[J].中國地震，8（2）：53-60.

郝錦綺，黃平章，周建國.1993.微破裂對巖石剩磁的影響——對地震預(yù)報的意義[J].地球物理學報，36（2）：203-211.

侯登錄，徐靜，唐貴德，等.2003.SrRuO3的磁化率隨溫度的變化關(guān)系[J].稀有金屬，27（5）：652-654.

黃平章，郝錦綺，周建國，等.1990.溧陽震區(qū)巖石磁性的溫度和應(yīng)力效應(yīng)[J].中國地震，6（1）：56-64.

蔣延林，趙衛(wèi)紅，趙永紅，等.2005.高郵地震臺地磁異常與江蘇及鄰區(qū)地震的關(guān)系研究[J].華北地震科學，23（1）：32-36.

解用明，鞠永，韓和平，等.2005.震前地磁幅值比異常研究[J].華北地震科學，23（1）：37-40.

解用明，喬子云，陳建國.2002.流動磁測在地震預(yù)報中的作用[J].內(nèi)陸地震，16（2）：170-174.

解用明，喬子云.2001.河北省及鄰區(qū)中等地震活動前的地磁異常[J].華北地震科學，19（4）：54-58.

李海燕，張世紅.2005.黃鐵礦加熱過程中的礦相變化研究——基于磁化率隨溫度變化特征分析[J].地球物理學報，48（6）：1384-1391.

李琪，楊冬梅，張素琴，等.2005.地磁轉(zhuǎn)換函數(shù)方法之典型震例及數(shù)模對比[J].地震，25（3）：85-90.

李樹華.2004.云南地磁場變化與異常特征分析[J].四川地震，（1）：31-36.

劉哲.2002.用地磁場總強度空間線性相關(guān)方法對遼寧省內(nèi)及周圍地區(qū)幾次中強地震的預(yù)報總結(jié)[J].華北地震研究，18（2）：46-50.

倪喆.2014a.洱源5.5級地震前后地磁場變化異常特征分析[J].地震研究，37（3）：426-432.

倪喆，陳雙貴，袁潔浩，等.2014b.蘆山7.0級地震前后巖石圈磁場異常變化研究[J].地震研究，37（1）：61-65.

倪喆，袁潔浩，王粲，等.2014c.2014年云南魯?shù)?.5級、永善5.0級地震前巖石圈磁場局部異常特征分析[J].地震研究，37（4）：537-541.

史勇軍，孫燕萍，朱燕，等.2002.伽師強震群活動過程中喀什臺地磁加卸載響應(yīng)比變化特征研究[J].內(nèi)陸地震，16（1）：84-88.

蘇樹鵬，李博，宋志剛，等.2016.晉冀蒙地磁總強度加密區(qū)觀測資料分析[J].地震地磁觀測與研究，37（2）：84-89.

王亶文.2003.地磁日變規(guī)律的破壞與地震的關(guān)系研究[J].地震地磁觀測與研究，24（1）：1-6.

王世昌，郅紅魁，宋建鎖.2004.地震孕育過程中地磁日變規(guī)律的分析[J].地震地磁觀測與研究，25（5）：18-23.

吳小平，辛華，張永久.2004.地磁“低點位移”與“磁暴”相結(jié)合預(yù)測發(fā)震時間[J].四川地震，（2）：21-26.

徐文耀，李衛(wèi)東.1994.地磁場Sq的經(jīng)度效應(yīng)和UT變化[J].地球物理學報，37（2）：157-166.

徐文耀.2003.地磁學[M].北京：地震出版社.

許忠淮.2014.地震學百科知識（九）中國淺源地震震源機制圖[J].國際地震動態(tài)，（3）：39-40.

楊濤，劉慶生，付媛媛，等.2005.隴西及周邊地區(qū)地震活動與震磁效應(yīng)的高階統(tǒng)計量研究[J].華北地震科學，23（1）：52-56.

詹景立，黃健.2004.鮮水河斷裂帶流動地磁變化特征與地震[J].四川地震，（2）：33-35.

張繼紅，喬惠珍.2005.地磁前兆異常的動態(tài)從屬函數(shù)[J].地球物理學進展，20（1）：142-146.

張淑蘭，江在森.2004.我國西部地磁絕對測量大臺距差值發(fā)對中強地震的反應(yīng)[J].地震，24（增刊1）：98-102.endprint