馬 亮

1 甘肅省地震局，蘭州市東崗西路450號(hào)，730000

地磁低點(diǎn)位移法是一種使用地磁日變數(shù)據(jù)預(yù)報(bào)地震的方法。陳紹明[1]在研究1966-03邢臺(tái)6.8級(jí)、7.2級(jí)地震余震的前兆時(shí)發(fā)現(xiàn)地磁異常現(xiàn)象，并于1972年提出“日變低點(diǎn)位移法”。丁鑒海等[2-3]在此基礎(chǔ)上通過(guò)研究全國(guó)地磁臺(tái)站日變低點(diǎn)位移的空間分布特點(diǎn)，于1975年提出“低點(diǎn)位移平面圖”方法，并以此開(kāi)展地震預(yù)報(bào)實(shí)踐工作。1986年陳紹明[4]在研究華北地區(qū)地磁日變曲線的反相位現(xiàn)象時(shí)提出“日變反向”的概念。丁鑒海[3]研究發(fā)現(xiàn)，低點(diǎn)時(shí)間均出現(xiàn)在臺(tái)站當(dāng)?shù)貢r(shí)間12時(shí)左右，于1988年提出低點(diǎn)時(shí)間的“經(jīng)度效應(yīng)”，并在1994年用數(shù)學(xué)歸納法推理低點(diǎn)位移法的判斷標(biāo)準(zhǔn)，從而解決該方法有經(jīng)驗(yàn)、無(wú)指標(biāo)的缺點(diǎn)。2009年王亞麗等[5]通過(guò)建立地磁低點(diǎn)時(shí)間的期望值與地理經(jīng)度的線性回歸方程指出，低點(diǎn)時(shí)間標(biāo)準(zhǔn)差與緯度間具有正相關(guān)關(guān)系。2017年許康生等[6]在研究岷縣-漳縣M6.6地震前地磁低點(diǎn)位移現(xiàn)象時(shí)發(fā)現(xiàn)，單個(gè)臺(tái)站的異常幅度與該臺(tái)站的震源距有關(guān)。2018年郭增建等[7]提出地磁低點(diǎn)位移以磁暴倍九法作為補(bǔ)充預(yù)測(cè)的觀點(diǎn)。賈昕曄等[8]使用“理論低點(diǎn)時(shí)間”與“實(shí)測(cè)低點(diǎn)時(shí)間”來(lái)區(qū)分地方時(shí)與世界時(shí)，并繪制低點(diǎn)時(shí)間垂直分量等值線圖。2019馬亮[9]提出隨經(jīng)度變化的地磁臺(tái)站當(dāng)?shù)貢r(shí)的計(jì)算公式。戴勇等[10]使用克里格插值法求出岷縣-漳縣M6.6地震前的低點(diǎn)位移事件的突變界線，這是首次脫離手繪方式求取突變界線。本文將使用數(shù)學(xué)語(yǔ)言及MATLAB語(yǔ)言，對(duì)地磁低點(diǎn)位移法與其搭載軟件進(jìn)行改進(jìn)，通過(guò)建立突變界線與臺(tái)網(wǎng)地磁低點(diǎn)時(shí)間之間的幾何解析關(guān)系來(lái)實(shí)現(xiàn)自動(dòng)成圖。

1 數(shù)據(jù)選取與問(wèn)題提出

中國(guó)地震臺(tái)網(wǎng)中心調(diào)用全國(guó)124個(gè)地磁臺(tái)站的數(shù)據(jù)來(lái)監(jiān)測(cè)地磁低點(diǎn)位移現(xiàn)象，圖1為臺(tái)站空間分布情況。東五區(qū)和東六區(qū)臺(tái)站密度小，經(jīng)向間距大；臺(tái)灣省、上海市與港澳地區(qū)未建設(shè)地磁臺(tái)站。該臺(tái)網(wǎng)中觀測(cè)儀包含F(xiàn)HD-2質(zhì)子矢量磁力儀、FHDZ-M15自動(dòng)化臺(tái)站系統(tǒng)、GM4磁通門(mén)磁力儀、CTM-DI磁通門(mén)磁力儀、MINGEO磁通門(mén)磁力儀、Mag-01H磁通門(mén)地磁經(jīng)緯儀、G856質(zhì)子(旋進(jìn))磁力儀和GSM-19FD overhauser磁力儀，本文選用的數(shù)據(jù)均為該臺(tái)網(wǎng)的地磁低點(diǎn)時(shí)間的垂直分量。由于地磁變化常受到空間磁暴的影響，本文對(duì)中國(guó)科學(xué)院國(guó)家空間科學(xué)中心空間環(huán)境預(yù)報(bào)中心(http://www.sepc.ac.cn/Kp3HPred_chn.php)提供的磁暴信息進(jìn)行篩除，以排除磁暴干擾。

采用等距離圓錐投影以確保經(jīng)向距離不變形圖1 地磁臺(tái)站空間分布Fig.1 Spatial distribution of geomagnetic stations

王亞麗等[5]在研究低點(diǎn)時(shí)間的經(jīng)度效應(yīng)時(shí)發(fā)現(xiàn)，垂直分量的低點(diǎn)時(shí)間期望值與經(jīng)度之間具有很強(qiáng)的負(fù)相關(guān)關(guān)系，這反映出低點(diǎn)時(shí)間具有顯著的地方時(shí)依賴(lài)性。由于每個(gè)測(cè)點(diǎn)的低點(diǎn)時(shí)間約為測(cè)點(diǎn)地方時(shí)12:00，經(jīng)度相差15°的兩個(gè)地磁臺(tái)站的世界時(shí)剛好相差1 h(圖1)。如果兩個(gè)臺(tái)站的經(jīng)向間距超過(guò)一個(gè)時(shí)區(qū)，突變界線的走向會(huì)出現(xiàn)較大誤差；如果兩個(gè)臺(tái)站的經(jīng)向間距超過(guò)兩個(gè)時(shí)區(qū)，在未出現(xiàn)地磁異常的情況下也可能存在地磁低點(diǎn)位移現(xiàn)象。為了消除經(jīng)度效應(yīng)帶來(lái)的誤差，需采用標(biāo)準(zhǔn)公式將世界時(shí)統(tǒng)一轉(zhuǎn)化為地方時(shí)，再劃分突變界線。

中國(guó)地震臺(tái)網(wǎng)中心一般通過(guò)地震前兆信息處理與分析軟件EIS2000調(diào)用地磁數(shù)據(jù)并繪制低點(diǎn)時(shí)間分布圖。EIS2000中地磁低點(diǎn)時(shí)間采用格林威治時(shí)間，手繪突變界線并進(jìn)行平滑處理，但該方法會(huì)產(chǎn)生兩種誤差。第一種誤差由經(jīng)度效應(yīng)引起，同一時(shí)區(qū)內(nèi)最西側(cè)與最東側(cè)的格林威治時(shí)間會(huì)自動(dòng)產(chǎn)生1 h偏差(圖1)，每個(gè)臺(tái)站格林威治時(shí)間的時(shí)差隨經(jīng)度而不是隨時(shí)區(qū)序數(shù)而變化，尤其是突變界線沿南北向延展時(shí)，格林威治時(shí)間對(duì)突變界限的走向會(huì)產(chǎn)生較大影響。第二種誤差為手繪突變界線產(chǎn)生的隨機(jī)誤差，同時(shí)在平滑處理過(guò)程中又會(huì)產(chǎn)生二次誤差。第一種誤差較為簡(jiǎn)單，可用臺(tái)站當(dāng)?shù)貢r(shí)間代替格林威治時(shí)間予以消除；第二種誤差處理難度較大，本文將采用新思路進(jìn)行處理。

2 Voronoi剖分的應(yīng)用

從集合的觀點(diǎn)來(lái)看，突變界線上任意一點(diǎn)到相鄰兩個(gè)臺(tái)站的距離相等，因此突變界線一定沿這兩個(gè)臺(tái)站的垂直平分線延伸。為獲取垂直平分線的坐標(biāo)，本文引入Voronoi剖分技術(shù)，用Voronoi剖分線代替突變界線的可能走向線。Voronoi圖(以下簡(jiǎn)稱(chēng)V圖)與Delaunay三角剖分圖互為偶圖。以某一散點(diǎn)為核心的V圖內(nèi)部部分被稱(chēng)為該點(diǎn)的值守范圍，V圖中相鄰兩個(gè)散點(diǎn)之間的垂直平分線稱(chēng)為Dirichlet鑲嵌線[11]，圖2為甘肅-青海地磁臺(tái)網(wǎng)的Delaunay三角剖分網(wǎng)和Dirichlet鑲嵌線。按照地磁臺(tái)站低點(diǎn)時(shí)間的大小對(duì)每個(gè)臺(tái)站的值守范圍進(jìn)行涂色，如果Dirichlet鑲嵌線兩側(cè)的色差超過(guò)2 h，地磁低點(diǎn)位移的突變界線便會(huì)沿該條Dirichlet鑲嵌線延展。

圖2 甘肅-青海地磁臺(tái)網(wǎng)的Delaunay三角網(wǎng)和Voronoi圖Fig.2 Delaunay triangulation and Voronoi diagram of the Gansu-Qinghai geomagnetic station network

3 新時(shí)間格式與原時(shí)間格式的計(jì)算結(jié)果對(duì)比

從地磁日變曲線可以直觀地看出，地磁低點(diǎn)時(shí)間大多位于當(dāng)?shù)貢r(shí)間12:00前后，且服從正態(tài)分布[9]。為消除采用格林威治時(shí)間對(duì)突變界線的走向造成的誤差，本節(jié)地磁低點(diǎn)時(shí)間統(tǒng)一采用當(dāng)?shù)貢r(shí)間，與格林威治時(shí)間的轉(zhuǎn)化關(guān)系為：

Tl=T0+N

(1)

式中，Tl、T0分別為臺(tái)站當(dāng)?shù)貢r(shí)間與格林威治時(shí)間，N為臺(tái)站所在時(shí)區(qū)的序數(shù)。若將Tl表示為T(mén)0與臺(tái)站經(jīng)度的函數(shù)，則式(1)可寫(xiě)為[9]：

(2)

式中，lon為地磁臺(tái)站的經(jīng)度。為避免低點(diǎn)時(shí)間隨經(jīng)度變化而對(duì)突變界線的走向產(chǎn)生影響，本節(jié)的地磁低點(diǎn)時(shí)間統(tǒng)一采用式(2)進(jìn)行求解，并非采用當(dāng)?shù)貢r(shí)區(qū)時(shí)間。由于T0已知，因此只需知道臺(tái)站經(jīng)度便可計(jì)算地磁低點(diǎn)時(shí)間的當(dāng)?shù)貢r(shí)間，進(jìn)而用集合思想計(jì)算出突變界線的坐標(biāo)數(shù)組。

為評(píng)價(jià)使用式(2)的新算法相對(duì)于原算法的優(yōu)勢(shì)，本文以2017年九寨溝7.0級(jí)地震發(fā)生前07-20出現(xiàn)的地磁低點(diǎn)位移現(xiàn)象為例，分析時(shí)間格式的改進(jìn)對(duì)地磁低點(diǎn)位移突變界線走向的影響。從圖1中可以看出，全國(guó)地磁臺(tái)站的空間分布跨越5個(gè)時(shí)區(qū)，如果采用格林威治時(shí)間計(jì)算突變界線，東五區(qū)和東七區(qū)臺(tái)站的地磁低點(diǎn)時(shí)間在正常情況下將會(huì)存在2 h偏差，為了直觀地表現(xiàn)臺(tái)站間的偏差，通過(guò)最近鄰插值法將全國(guó)124個(gè)采用格林威治時(shí)間的地磁臺(tái)的低點(diǎn)時(shí)間繪制成空間等值線分布圖(圖3)。

由于世界時(shí)會(huì)隨經(jīng)度而變化，在地磁數(shù)據(jù)無(wú)異常時(shí)，臺(tái)站經(jīng)向跨度較大的地區(qū)低點(diǎn)時(shí)間可能會(huì)相差2 h以上，從而產(chǎn)生“假異常”現(xiàn)象。圖3中新疆5個(gè)地磁臺(tái)外邊緣的突變界線就屬于“假異?！蓖蛔兘缇€，“假異?！笔怯墒澜鐣r(shí)的固有性質(zhì)而產(chǎn)生，與場(chǎng)地條件無(wú)關(guān)，采用本文方法中式(2)將世界時(shí)轉(zhuǎn)化為當(dāng)?shù)貢r(shí)后即可消除。

圖3 采用格林威治時(shí)間的突變界線Fig.3 Abrupt boundaries based on GMT

各臺(tái)站地質(zhì)構(gòu)造、地下電性結(jié)構(gòu)存在差異，臺(tái)站資料差異對(duì)突變界線的影響已體現(xiàn)在原始數(shù)據(jù)中。因此只要使用相同的數(shù)據(jù)，無(wú)論采用原始手繪方法還是本文方法，在去除手繪的隨機(jī)誤差后，突變界線走向一致。手繪的原始突變界線為折線，但經(jīng)過(guò)后期平滑處理后可轉(zhuǎn)化為平滑曲線。

新疆、西藏、青海、甘肅地區(qū)地磁臺(tái)站稀疏，各臺(tái)站經(jīng)向跨度大，按照格林威治時(shí)間計(jì)算的突變界線的誤差會(huì)被放大，為體現(xiàn)與式(2)地方時(shí)的差異，本節(jié)按式(2)將圖3中中國(guó)西部臺(tái)站的原有低點(diǎn)時(shí)間轉(zhuǎn)化為地方時(shí)，并求出兩者的偏差值(表1)。

表1 西部地磁臺(tái)網(wǎng)的格林威治時(shí)間與地方時(shí)的低點(diǎn)時(shí)間對(duì)比Tab.1 Comparison of low point-time between GMT and local time of western geomagnetic network

從表1可以看出，獅泉河臺(tái)與喀什臺(tái)的低點(diǎn)時(shí)間與背景值之差超過(guò)60 min，各臺(tái)站格林威治地磁低點(diǎn)時(shí)間差異非常大，且經(jīng)度跨度越大，地磁低點(diǎn)時(shí)間的差異越大。將格林威治時(shí)間轉(zhuǎn)化為式(2)地方時(shí)后，各臺(tái)站的低點(diǎn)時(shí)間約為12:00，與文獻(xiàn)[3]研究結(jié)果一致，表明新方法可有效消除“經(jīng)度效應(yīng)”。但文獻(xiàn)[3]未定義地方時(shí)的計(jì)算公式，該研究中“地方時(shí)”概念可能源于臺(tái)站的時(shí)區(qū)時(shí)間，因此本文建議將各臺(tái)站按式(2)計(jì)算的地方時(shí)的地磁低點(diǎn)時(shí)間長(zhǎng)期觀測(cè)的期望定義為地磁低點(diǎn)時(shí)間的背景值，以消除原有方法的“經(jīng)度效應(yīng)”。

4 突變界線的定位

現(xiàn)虛擬一個(gè)由8個(gè)地磁臺(tái)站組成的地磁臺(tái)網(wǎng)，臺(tái)站緯、經(jīng)度分別用矩陣lat和lon表示，低點(diǎn)時(shí)間用矩陣data表示，突變界線在MATLAB中的坐標(biāo)計(jì)算代碼如下：

lat=[y1,y2,y3,y4,y5,y6,y7,y8];

lon=[x1,x2,x3,x4,x5,x6,x7,x8];

data=[a1,a2,a3,a4,a5,a6,a7,a8];

[X,Y]=meshgrid(linspace(y1,y8,n),

linspace(x1,x8,n));

Z=griddata(lat,lon,data,X,Y, ‘nearest’);

[C,h,CF]=contourf(Y,X,Z, 2.5);

其中，命令linspace(y1,y8,n)用來(lái)產(chǎn)生y1到y(tǒng)8之間等距離的n個(gè)點(diǎn)。[X,Y]=meshgrid(X,Y)(X=[x1,x2,x3],Y=[y1,y2,y3])函數(shù)用來(lái)產(chǎn)生網(wǎng)格矩陣，X為網(wǎng)格矩陣的橫坐標(biāo)范圍，Y是網(wǎng)格矩陣的縱坐標(biāo)范圍；返回網(wǎng)格矩陣的坐標(biāo)[X,Y]。Z=griddata(lat,lon,data,X,Y, ‘n’)為散點(diǎn)插值函數(shù)，用于在散點(diǎn)[lat,lon]之間進(jìn)行插值[11]，data為該點(diǎn)的色值，[X,Y]表示插值的網(wǎng)格(密度遠(yuǎn)大于[lat,lon])，n為插值類(lèi)型；返回插值結(jié)果矩陣Z。函數(shù)[C,h,CF]=contourf(Y,X,Z,m)可用來(lái)繪制二維等值線圖，Z為平面內(nèi)網(wǎng)格點(diǎn)[Y,X]處的色值，m為梯度線的個(gè)數(shù)，添加色彩后會(huì)產(chǎn)生m+1種顏色。同時(shí)返回與命令contourc中相同的等高線矩陣C，C也可被命令clabel使用；返回包含patch圖形對(duì)象的句柄向量h；返回填充色值的矩陣CF。

突變界線可通過(guò)設(shè)置參數(shù)m，將等值線階數(shù)降為1階來(lái)實(shí)現(xiàn)。將contourf命令中的m值設(shè)為[11]：

m=n-1

(3)

(4)

式中，m為等值線條數(shù)，n為梯度數(shù)，d為總梯度差，L為步長(zhǎng)。例如在某次地震事件中，低點(diǎn)時(shí)間等值線條數(shù)為m，所涂顏色共有n種，網(wǎng)內(nèi)地磁低點(diǎn)時(shí)間的最早值與最晚值之差為d，兩個(gè)區(qū)域之間的低點(diǎn)時(shí)間相差步長(zhǎng)為L(zhǎng)，且L≡2。例如網(wǎng)內(nèi)低點(diǎn)時(shí)間的最早值與最晚值分別為08:00和15:00，此時(shí)d=7，將d和L≡2代入式(3)和式(4)中求出m值為2.5。執(zhí)行contourf命令可輸出以m為參數(shù)的地磁低點(diǎn)時(shí)間的等值線，突變界線為其中一條或一段。

式(3)中m并不是突變界線的條數(shù)，而是確定等值線的參數(shù)。在地磁臺(tái)網(wǎng)中，所有臺(tái)站的低點(diǎn)時(shí)間的相互差異程度決定m值的范圍，當(dāng)0

在2015年皮山6.5級(jí)地震前06-21地磁低點(diǎn)位移事件中，m值為6，等值線有5條，突變界線有1條。在2016年門(mén)源6.4級(jí)地震前2015-12-23甘肅、青海地區(qū)出現(xiàn)的地磁低點(diǎn)位移現(xiàn)象中，m值為2，等值線有2條，突變界線有1條。在2016年甘肅金塔4.7級(jí)地震前02-28地磁低點(diǎn)位移事件中，m值為3，等值線有3條，突變界線有2條，本節(jié)案例參見(jiàn)文獻(xiàn)[11]。

5 數(shù)據(jù)分析

中國(guó)地震臺(tái)網(wǎng)中心于2017-07-20在華北、長(zhǎng)江中下游地區(qū)出現(xiàn)的地磁低點(diǎn)位移現(xiàn)象中發(fā)現(xiàn)，低點(diǎn)時(shí)間被分為內(nèi)外兩個(gè)區(qū)域，按原數(shù)據(jù)的格林威治時(shí)間統(tǒng)計(jì)的內(nèi)區(qū)域平均低點(diǎn)時(shí)間為6.73時(shí)，外區(qū)域平均低點(diǎn)時(shí)間為3.84時(shí)，兩者相差2.89 h；將格林威治時(shí)間轉(zhuǎn)化為式(2)地方時(shí)后，內(nèi)區(qū)域的平均低點(diǎn)時(shí)間為14.39時(shí)(表2)。外區(qū)域的平均低點(diǎn)時(shí)間為11.01時(shí)，與內(nèi)區(qū)域相差3.38 h。桃源臺(tái)儀器出現(xiàn)故障，因此用邵陽(yáng)臺(tái)地磁數(shù)據(jù)代替。在地磁低點(diǎn)位移現(xiàn)象發(fā)生后的第19天，即北京時(shí)間2017-08-08 21:19，四川省阿壩州九寨溝縣發(fā)生MS7.0地震，比文獻(xiàn)[3]中統(tǒng)計(jì)的極可能發(fā)震區(qū)間的最小值小4 d，震中距突變界線約185 km。2017-08-08的K指數(shù)為2，未受到外空間磁暴干擾，地磁低點(diǎn)數(shù)據(jù)真實(shí)可靠，因此認(rèn)為，07-20出現(xiàn)的地磁低點(diǎn)位移現(xiàn)象與此次地震具有較好的相關(guān)性，圖4為中國(guó)地震臺(tái)網(wǎng)中心采用人工繪制的突變界線。

表2 內(nèi)區(qū)域全部臺(tái)站的低點(diǎn)時(shí)間Tab.2 Low-point time of all stations in inner region

圖4 2017-07-20出現(xiàn)的地磁低點(diǎn)位移現(xiàn)象(采用當(dāng)?shù)貢r(shí)間)Fig.4 The geomagnetic low-point displacement phenomenon on July 20, 2017(use local time)

突變界線由乳山、濟(jì)南、紅山、定襄、太原、臨汾、萬(wàn)州、石柱、涪陵、重慶、桃源、九峰、高郵、海安臺(tái)共同組成，外區(qū)域西端存在榆林、漢中、成都、美姑、貴陽(yáng)、邵陽(yáng)6個(gè)控制性臺(tái)站，對(duì)突變界線的走向起決定性作用。為了直觀地表現(xiàn)兩種形式的低點(diǎn)時(shí)間在空間分布上的差異以及§4中自動(dòng)定位技術(shù)提取與人工畫(huà)線的差異，采用式(2)地方時(shí)結(jié)合§5方法求出全國(guó)124個(gè)臺(tái)站的低點(diǎn)時(shí)間等值線圖(圖5)。

圖5 采用當(dāng)?shù)貢r(shí)間繪制的突變界線Fig.5 Abrupt boundaries based on local time

由圖5可見(jiàn)，最近鄰插值結(jié)果中共有4條梯度線，第1條與人工繪制的突變界線形狀相似，第2條在襄陽(yáng)市、南陽(yáng)市周?chē)?，?條在山東省無(wú)棣大山地震臺(tái)周?chē)?條在新疆維吾爾自治區(qū)喀什市周?chē)?。由于單個(gè)地磁臺(tái)站低點(diǎn)時(shí)間的偏移未形成大區(qū)域，因此認(rèn)為，圖5中第2～4條梯度線不是地磁低點(diǎn)位移突變界線，只有第1條梯度線符合突變界線的定義。

圖5為使用式(2)地方時(shí)的結(jié)果，喀什臺(tái)與都蘭臺(tái)的相對(duì)偏差減小為185 min，烏什臺(tái)、紅淺臺(tái)、烏魯木齊臺(tái)與都蘭臺(tái)的低點(diǎn)時(shí)間色值相同。圖3為使用格林威治時(shí)間的結(jié)果，喀什臺(tái)和都蘭臺(tái)低點(diǎn)時(shí)間色值不同，相對(duì)偏差為274 min，烏什臺(tái)、紅淺臺(tái)、烏魯木齊臺(tái)和都蘭臺(tái)低點(diǎn)時(shí)間色值不同。上述分析表明，新方法可減小具有經(jīng)向跨度的臺(tái)站之間低點(diǎn)時(shí)間的偏差，可過(guò)濾掉新疆地區(qū)的“假異?！蓖蛔兘缇€。

對(duì)比圖4與圖5可知，前者是平滑曲線，后者為折線。圖4中東北地區(qū)的第2條突變界線在人工繪制時(shí)未被發(fā)現(xiàn)，而圖5中已繪制東北地區(qū)的突變界線。在華北、華中地區(qū)，兩個(gè)圖的突變界線走向大體一致，在河北北部略有差異。圖5中的突變界線在內(nèi)區(qū)域存在兩個(gè)孤島，由于單個(gè)地磁臺(tái)站低點(diǎn)時(shí)間的偏移未形成大區(qū)域，因此本文未把孤島外輪廓作為地磁低點(diǎn)位移突變界線。

6 結(jié) 語(yǔ)

本文提出的方法在不改變地磁低點(diǎn)位移法基本原理的情況下可有效排除干擾，減少誤差。九寨溝地震前的地磁低點(diǎn)位移現(xiàn)象分析表明，本文提出的突變界線定位法與原有算法相比具有以下優(yōu)點(diǎn)：

1)原有方法會(huì)產(chǎn)生“經(jīng)度效應(yīng)”， 新方法可很好地消除。

2)原有方法的突變界線可以反映突變界線的大致走向，但無(wú)法獲取突變界線內(nèi)區(qū)域的孤島；而新方法可以精確地繪制突變界線內(nèi)區(qū)域中的孤島。

3)原有地磁低點(diǎn)時(shí)間的數(shù)據(jù)格式無(wú)法統(tǒng)計(jì)與地磁低點(diǎn)背景時(shí)間的偏差，并且使用概率統(tǒng)計(jì)法識(shí)別異常會(huì)受來(lái)自外空?qǐng)龅母蓴_所影響，從而較難識(shí)別真正的異常；而新方法中的低點(diǎn)時(shí)間采用式(2)地方時(shí)，可以更精確地識(shí)別異常和來(lái)自外空?qǐng)龅母蓴_，過(guò)濾 “假異?！蓖蛔兘缇€。

4)如果將線條看作點(diǎn)的集合，且假定原方法手繪的突變界線足夠精確，則原方法的突變界線為新方法的突變界線的非空真子集，后者包含前者，人工劃線時(shí)容易遺漏一些隱藏的突變界線，而新方法不會(huì)遺漏隱藏的突變界線。