?
區(qū)域地磁場(chǎng)諧和樣條模型研究
馮麗麗
(中國(guó)地震局地球物理研究所,北京100081)
地震過程是巖石圈能量積累、聚集和釋放的過程。伴隨此過程,巖石圈磁場(chǎng)會(huì)出現(xiàn)異常變化。近年來,我國(guó)的地震預(yù)報(bào)研究實(shí)踐證明,巖石圈磁場(chǎng)異常變化與地震發(fā)生地點(diǎn)有很好的對(duì)應(yīng)關(guān)系,并取得了良好效果。
準(zhǔn)確提取地震引起的巖石圈磁場(chǎng)變化的主要途徑是建立高分辨率的巖石圈磁場(chǎng)模型。由于針對(duì)地震預(yù)報(bào)的觀測(cè)都是區(qū)域性的,因此,模型計(jì)算方法應(yīng)同時(shí)具備3個(gè)方面的優(yōu)點(diǎn): 一是適合區(qū)域分析; 二是能準(zhǔn)確剔除或剝離地磁主磁場(chǎng); 三是能精細(xì)刻畫巖石圈磁場(chǎng)細(xì)微變化??赡艿脑?,模型還應(yīng)能采用衛(wèi)星數(shù)據(jù)與地面數(shù)據(jù)聯(lián)合運(yùn)算,以彌補(bǔ)我國(guó)西部無人區(qū)地面觀測(cè)數(shù)據(jù)稀疏的不足及控制邊界畸變,即模型應(yīng)能徑向延拓。但在實(shí)際工作中,許多模型很難三者兼顧,更別說兼具徑向延拓。該論文從模型比較開始,精心選擇模型,經(jīng)過模型設(shè)計(jì)和程序?qū)崿F(xiàn)、模型調(diào)優(yōu)和模型驗(yàn)證等,初步獲得了能四者兼顧的中國(guó)地區(qū)區(qū)域地磁場(chǎng) “諧和樣條” 模型及其計(jì)算方法。論文最后將中國(guó)地區(qū)區(qū)域地磁場(chǎng) “諧和樣條” 模型與其他模型進(jìn)行了對(duì)比分析,發(fā)現(xiàn)該模型初步達(dá)到了設(shè)計(jì)目標(biāo)。
(1) 模型選擇: 論文調(diào)研和對(duì)比分析了球諧分析、矩諧分析、球冠諧模型、修正球冠模型、等效源、多項(xiàng)式擬合、樣條函數(shù)擬合、雙調(diào)和擬合、諧和樣條擬合。論文分析認(rèn)為,調(diào)和樣條擬合具備前述的4個(gè)方面優(yōu)點(diǎn)。
(2) 模型設(shè)計(jì)和程序?qū)崿F(xiàn): 論文完成了對(duì) “諧和樣條” 方法的數(shù)學(xué)公式推導(dǎo)、代數(shù)公式構(gòu)建以及matlab程序編寫等工作。
(3) 模型調(diào)優(yōu): 論文完成了 “諧和樣條” 模型定解過程中的 “本征值截?cái)唷?以及球諧階數(shù)問題分析討論,發(fā)現(xiàn) “諧和樣條” 方法定解線性方程組本征值可能對(duì)應(yīng) “諧和樣條” 模型所描述場(chǎng)分布的空間特征尺度。因此,通過本征值截?cái)嚯A數(shù)的選擇,能夠?qū)崿F(xiàn)對(duì)不同空間特征尺度磁場(chǎng)成分的剝離。
(4) 模型驗(yàn)證: 論文基于大量的計(jì)算實(shí)驗(yàn),選擇了球諧截?cái)嚯A數(shù)L=90階,并分別利用地面測(cè)量數(shù)據(jù)以及衛(wèi)星數(shù)據(jù)的組合分別定解,計(jì)算對(duì)比分析了 “中國(guó)及周邊區(qū)域諧和樣條模型” 。結(jié)果表明: 與球諧方法相比,“諧和樣條” 模型具有更高的 “效率” ,引入衛(wèi)星數(shù)據(jù)確能改善計(jì)算結(jié)果。此外,當(dāng) “諧和樣條” 方法用于全球模型和區(qū)域模型計(jì)算時(shí),發(fā)現(xiàn)球諧截?cái)嚯A數(shù)的確定是影響模型計(jì)算正確性的關(guān)鍵因素之一。
(5) 與其他模型對(duì)比分析: 論文分別以1~15階EMM2010模型和 “2010.0中國(guó)及周邊區(qū)域諧和樣條模型” 20階本征值模型作為中國(guó)區(qū)域地磁正常場(chǎng),從90階 “2010.0中國(guó)及周邊區(qū)域諧和樣條模型” 中分別剔除兩種算法得到的中國(guó)區(qū)域地磁正常場(chǎng),獲得2個(gè)中國(guó)區(qū)域地磁異常場(chǎng)模型,對(duì)比分析了其差異性和優(yōu)劣性,論文最后采用 “2010.0中國(guó)及周邊區(qū)域諧和樣條模型” 地磁異常場(chǎng)、EMM2010模型、MF6模型和CGRF模型的巖石圈磁場(chǎng)Z分量,與中國(guó)區(qū)域地質(zhì)構(gòu)造進(jìn)行了聯(lián)合對(duì)比分析,達(dá)到了論文初期目標(biāo)。
關(guān)鍵詞諧和樣條; 區(qū)域地磁場(chǎng); 截?cái)嚯A數(shù); 本征值截?cái)?/p>
(作者電子信箱,馮麗麗: ynufll@sina.com)
中圖分類號(hào):P318;
文獻(xiàn)標(biāo)識(shí)碼:A;
doi:10.3969/j.issn.0235-4975.2016.01.010