張向宇， 關(guān)永賢， 崔秀云

(國(guó)土資源部 海底礦產(chǎn)資源重點(diǎn)實(shí)驗(yàn)室，中國(guó)地質(zhì)調(diào)查局廣州海洋地質(zhì)調(diào)查局，廣州 510075)

利用回歸分析法補(bǔ)償海底日變數(shù)據(jù)

張向宇， 關(guān)永賢， 崔秀云

(國(guó)土資源部 海底礦產(chǎn)資源重點(diǎn)實(shí)驗(yàn)室，中國(guó)地質(zhì)調(diào)查局廣州海洋地質(zhì)調(diào)查局，廣州 510075)

海洋磁力測(cè)量中常會(huì)投放海底日變站，但有時(shí)會(huì)因投放深度較深導(dǎo)致回收的海底日變站數(shù)據(jù)出現(xiàn)失真的情況，從而影響日變改正后磁測(cè)數(shù)據(jù)的精度。通過比較不同緯度處地磁日變數(shù)據(jù)的相關(guān)分析結(jié)果，發(fā)現(xiàn)可近似認(rèn)為地磁日變數(shù)據(jù)依緯度變化具有一定的線性相關(guān)性，進(jìn)而引入線性回歸分析方法推算某一緯度處日變數(shù)據(jù)。通過南海多個(gè)工區(qū)實(shí)測(cè)數(shù)據(jù)的對(duì)比分析，發(fā)現(xiàn)回歸分析法推算的日變數(shù)據(jù)可一定程度上補(bǔ)償海底日變數(shù)據(jù)，從而改善磁測(cè)數(shù)據(jù)的精度。

海洋磁測(cè); 海底日變站; 回歸分析; 日變改正

0 引言

在海洋磁測(cè)中，地磁日變站的控制范圍是有限的，在海洋調(diào)查規(guī)范中，要求地磁日變站有效控制半徑為300 km～500 km[1]，而南海周邊可用陸地地磁臺(tái)站有限，部分區(qū)域離陸地臺(tái)站距離較遠(yuǎn)，超過了臺(tái)站控制范圍，針對(duì)這種情況，野外測(cè)量時(shí)會(huì)根據(jù)情況選擇投放海底日變站，但受多種因素制約，使得回收的數(shù)據(jù)會(huì)出現(xiàn)因投放深度及海水層過濾作用等原因?qū)е率д?，尤其是?dāng)日變站投放深度超過1 000 m時(shí)，這種情況會(huì)更明顯。另外回收的數(shù)據(jù)經(jīng)常會(huì)出現(xiàn)缺失，這些都將影響日變改正效果，進(jìn)而影響改正后磁測(cè)數(shù)據(jù)的精度，因此，找尋合適的方法補(bǔ)償海底日變數(shù)據(jù)是很有必要的[2-5]。

這里將采用推算日變數(shù)據(jù)的方法來達(dá)到補(bǔ)償?shù)哪康摹＝陙?，?guó)內(nèi)學(xué)者提出了多種推算方法：?jiǎn)稳陜€等[6]提出了幾種距離加權(quán)推算法，包括二維多項(xiàng)式最小二乘擬合法、時(shí)空擬合法和線性內(nèi)插法對(duì)局部地磁日變數(shù)據(jù)進(jìn)行擬合；邊剛等[7]采用多站同步實(shí)例分析了加權(quán)平均法和函數(shù)擬合法在海洋磁力測(cè)量地磁日變改正中的應(yīng)用；卞光浪等[8]提出了脫離距離加權(quán)法的基于緯差加權(quán)法的海洋磁力測(cè)量多站地磁日變改正值計(jì)算。通過對(duì)沿岸幾個(gè)日變站的數(shù)據(jù)進(jìn)行計(jì)算證實(shí)了該方法效果優(yōu)于距離加權(quán)法。以上各種推算方法都各具特點(diǎn)，應(yīng)用效果也不盡相同，直接用來補(bǔ)償海底日變數(shù)據(jù)都不甚理想，因此，需要在已有推算方法的基礎(chǔ)上做一定地改進(jìn)。

筆者從緯差加權(quán)法推算公式入手，并根據(jù)不同臺(tái)站地磁日變數(shù)據(jù)相關(guān)性分析的結(jié)果，發(fā)現(xiàn)不同緯度處的日變數(shù)據(jù)可近似認(rèn)為是呈線性相關(guān)的，因此可通過統(tǒng)計(jì)學(xué)中常用的線性回歸分析方法擬合出日變數(shù)據(jù)推算公式，進(jìn)而得到待推算位置處的地磁日變數(shù)據(jù)。選取南海海域投放過海底日變站的幾個(gè)工區(qū)的實(shí)測(cè)數(shù)據(jù)，由與海底日變站數(shù)據(jù)相關(guān)性較好的兩個(gè)陸地臺(tái)站數(shù)據(jù)，通過線性回歸分析推算出投放位置處的地磁日變數(shù)據(jù)，對(duì)比幾組日變數(shù)據(jù)及改正后數(shù)據(jù)的交點(diǎn)差，發(fā)現(xiàn)線性回歸分析法得到的數(shù)據(jù)對(duì)海底日變數(shù)據(jù)進(jìn)行了一定程度的補(bǔ)償，提升了改正后磁測(cè)數(shù)據(jù)精度，是對(duì)多臺(tái)站地磁日變數(shù)據(jù)推算方法地有效改進(jìn)。

1 回歸分析法原理

緯差加權(quán)法[8]假設(shè)測(cè)點(diǎn)P的地磁日變改正值是由局部鄰域內(nèi)的多個(gè)同步地磁日變站檢測(cè)信息加權(quán)平均得到。設(shè)地磁日變站坐標(biāo)為(φi,λj)，在t時(shí)刻各站的地磁日變改正值為ΔTi(φi,λj,t)，記各站權(quán)函數(shù)為Wi(φ,λ)，則測(cè)點(diǎn)的地磁日變改正值可表示為式(1)。

(1)

設(shè)權(quán)函數(shù)Wi(φ,λ)與緯差的μ(μ≥0)次方的倒數(shù)成正比，則權(quán)函數(shù)表達(dá)式為式(2)。

在文獻(xiàn)[8]中，應(yīng)用某三個(gè)臺(tái)站的實(shí)際地磁日變數(shù)據(jù)對(duì)式(1)中μ和k分別取0、1、2時(shí)緯差加權(quán)法的計(jì)算精度進(jìn)行統(tǒng)計(jì)，發(fā)現(xiàn)當(dāng)μ=1、k=1時(shí)計(jì)算出的中誤差最小，這時(shí)測(cè)點(diǎn)P的地磁日變改正值為：

ΔT(φ,λ,t)=

(3)

其中：ΔT1(φ1,λ1,t)為A臺(tái)站地磁日變改正值； ΔT2(φ2,λ2,t)為B臺(tái)站地磁日變改正值；φA、φB分別為A臺(tái)站和B臺(tái)站的緯度；ΔφAP、ΔφBP分別為測(cè)點(diǎn)P與A站和B站的緯差。

在式(3)中，三個(gè)日變站的日變數(shù)據(jù)間是呈現(xiàn)線性關(guān)系的，通過對(duì)多組實(shí)測(cè)地磁日變數(shù)據(jù)地分析發(fā)現(xiàn)，不同臺(tái)站的日變數(shù)據(jù)具有一定的相關(guān)性，緯差越小的臺(tái)站數(shù)據(jù)相關(guān)性越好，因此通過多方面地分析，可近似認(rèn)為地磁臺(tái)站間的日變數(shù)據(jù)是線性相關(guān)的[9]，即

ΔT=aΔT1+bΔT2

(4)

式中：a、b為常數(shù)；ΔT1、ΔT2分別為兩個(gè)基臺(tái)站的地磁場(chǎng)值；ΔT為待推算位置處地磁場(chǎng)值，公式的關(guān)鍵在于求取常數(shù)a、b的值，可以利用統(tǒng)計(jì)學(xué)中的線性回歸分析方法。

因此，在知曉兩個(gè)基臺(tái)站(ΔT1和ΔT2)及待推算位置處(ΔT)相同時(shí)段的日變數(shù)據(jù)后，可選取適當(dāng)時(shí)段數(shù)據(jù)作為樣本，通過線性回歸分析方法擬合出公式(4)，即利用最小二乘法對(duì)參數(shù)a、b進(jìn)行無偏估計(jì)，進(jìn)而得到最終的推算公式，這樣便可求取任一時(shí)段待推算位置處的日變數(shù)據(jù)。在做線性回歸分析前應(yīng)先對(duì)樣本數(shù)據(jù)進(jìn)行相關(guān)性分析，當(dāng)相關(guān)系數(shù)較大，一般在0.9以上時(shí)才適宜做線性回歸分析。

從不同位置處的地磁日變曲線形態(tài)可看到，地磁日變化受緯度影響較大，不同緯度處的日變曲線幅值相差較大，而經(jīng)度的不同使日變曲線的相位產(chǎn)生移動(dòng)，且國(guó)、內(nèi)外多位學(xué)者研究成果均認(rèn)為通常情況下，緯度效應(yīng)引起的地磁日變化較經(jīng)度效應(yīng)大很多，因此這里我們只關(guān)注緯度的影響。由文獻(xiàn)[9]可知應(yīng)用緯差加權(quán)法時(shí)，兩基臺(tái)站與待推算位置的緯度越接近推算結(jié)果越準(zhǔn)確，兩基臺(tái)站緯度差應(yīng)盡量控制在10°內(nèi)，外推時(shí)兩基臺(tái)站中要有一個(gè)距離待推算位置處的緯差不超過3°；回歸分析法只能在待推算位置處有作業(yè)期間相鄰時(shí)間段內(nèi)的日變數(shù)據(jù)來做樣本數(shù)據(jù)才可使用，在內(nèi)插時(shí)其偏差較緯差加權(quán)法小，外推時(shí)當(dāng)待推算位置距離最近基臺(tái)站緯度差超過3°時(shí)，偏差明顯小于緯差加權(quán)法。另外使用回歸分析法時(shí)，樣本數(shù)據(jù)應(yīng)盡量保證選擇連續(xù)磁平靜日的數(shù)據(jù)，且樣本數(shù)據(jù)量越大擬合結(jié)果越準(zhǔn)確[10]。

2 實(shí)例分析

選取四個(gè)南海海域工區(qū)實(shí)測(cè)數(shù)據(jù)，這些工區(qū)在野外作業(yè)期間均投放了海底日變站，且投放深度各不相同。各工區(qū)計(jì)算時(shí)需用到南海周邊幾個(gè)地磁臺(tái)站的數(shù)據(jù)，這些臺(tái)站數(shù)據(jù)均來源于國(guó)際地磁臺(tái)網(wǎng)中心網(wǎng)站，各臺(tái)站信息見表1，相對(duì)位置見圖1。

表1 地磁臺(tái)站信息表

圖1 地磁臺(tái)位置示意圖Fig.1 Diagram of the geomagnetic stations

2.1 工區(qū)一計(jì)算結(jié)果

該工區(qū)位于南海中西部，調(diào)查測(cè)線在112°E～114°E、12°N～16°N范圍內(nèi)布設(shè)，野外采集時(shí)間為5月11日至5月29日，所用磁力儀為Seaspy海洋磁力儀，觀測(cè)精度為0.25 nT。野外測(cè)量時(shí)在工區(qū)附近投放了Sentinel地磁日變觀測(cè)站，5月11日在112.45°E，16.85°N，水深917 m處投放，5月27日回收，缺失5月28日、29日及部分27日的數(shù)據(jù)。

將該工區(qū)投放的海底日變站(代號(hào)HD-1)與陸地臺(tái)站數(shù)據(jù)做相關(guān)性分析，結(jié)果見表2。

表2 工區(qū)一日變數(shù)據(jù)相關(guān)系數(shù)統(tǒng)計(jì)表

從表2可以看出，與海底日變數(shù)據(jù)相關(guān)性最好的為越南DLT地磁臺(tái)和海南瓊中地磁臺(tái)數(shù)據(jù)，其相關(guān)系數(shù)均在0.97以上，因此選擇這兩個(gè)地磁臺(tái)推算海底日變站緯度(16.85°N)處的地磁日變數(shù)據(jù)。

選取5月12日至5月27日期間海底日變數(shù)據(jù)、越南DLT地磁臺(tái)數(shù)據(jù)和海南瓊中地磁臺(tái)數(shù)據(jù)作為樣本數(shù)據(jù)，推算出16.85°N處的日變數(shù)據(jù)，選取5月19日的數(shù)據(jù)成圖，同時(shí)為了有所對(duì)比，由式(3)計(jì)算16.85°N處的緯差加權(quán)法推算結(jié)果，成圖前將三組日變數(shù)據(jù)基值進(jìn)行統(tǒng)一(圖2)。

圖2 5月19日推算日變數(shù)據(jù)對(duì)比圖Fig.2 Comparison for speculated diurnal data in May 19th

從圖2可以看出，回歸分析法推算的數(shù)據(jù)與緯差加權(quán)法推算的數(shù)據(jù)較接近，海底日變數(shù)據(jù)較推算數(shù)據(jù)并未損失掉較多的細(xì)節(jié)信息，但幅值略小。

對(duì)野外采集到的實(shí)測(cè)數(shù)據(jù)依照海洋地質(zhì)調(diào)查規(guī)范要求做各項(xiàng)改正后，將各交點(diǎn)處的交點(diǎn)差展布于測(cè)線分布圖上，發(fā)現(xiàn)交點(diǎn)差有所減小，統(tǒng)計(jì)24個(gè)交點(diǎn)處不同數(shù)據(jù)日變改正后的測(cè)線交點(diǎn)差，結(jié)果見表3。

從表3可以看出，回歸分析法推算的日變數(shù)據(jù)改正后交點(diǎn)差小于海底日變改正后結(jié)果，且明顯優(yōu)于緯差加權(quán)法改正后結(jié)果，同時(shí)以回歸分析法得到的處理后數(shù)據(jù)精度最高(處理后數(shù)據(jù)交點(diǎn)差均方根差值，即磁異常ΔT精度)。

表3 工區(qū)一數(shù)據(jù)交點(diǎn)差統(tǒng)計(jì)表

2.2 工區(qū)二計(jì)算結(jié)果

該工區(qū)位于南海中部，調(diào)查測(cè)線在112°E～118°E、12°N～16°N范圍內(nèi)布設(shè)，野外采集日期為5月16日至7月8日，所用磁力儀為G-882SX海洋磁力儀，觀測(cè)精度小于0.3 nT。野外測(cè)量時(shí)在工區(qū)附近投放了Sentinel地磁日變觀測(cè)站，5月16日在113.05°E，14.6°N，水深2 810 m處投放，7月8日收回觀測(cè)站?；厥盏娜兆償?shù)據(jù)在6月13至14日及6月17至21日期間出現(xiàn)缺失。

對(duì)表1中所列四個(gè)地磁臺(tái)數(shù)據(jù)連同該工區(qū)投放的海底日變站(代號(hào)HD-2)數(shù)據(jù)做相關(guān)性分析，樣本數(shù)據(jù)選擇與海底日變回收數(shù)據(jù)同時(shí)間段的各地磁臺(tái)站數(shù)據(jù)，結(jié)果見表4。

從表4可以看出，與海底日變站數(shù)據(jù)相關(guān)系數(shù)最高的為越南DLT地磁臺(tái)和海南瓊中地磁臺(tái)數(shù)據(jù)，其相關(guān)系數(shù)均在0.9以上，因此，本次選擇由這兩個(gè)地磁臺(tái)通過線性回歸分析方法，推算海底日變站緯度(14.6°N)處的地磁日變數(shù)據(jù)。

表4 工區(qū)二日變數(shù)據(jù)相關(guān)系數(shù)統(tǒng)計(jì)表

選取5月16日至7月8日期間海底日變數(shù)據(jù)、越南DLT地磁臺(tái)數(shù)據(jù)和海南瓊中地磁臺(tái)數(shù)據(jù)作為樣本數(shù)據(jù)，進(jìn)行線性回歸分析，得到14.6°N的日變數(shù)據(jù)擬合公式，同樣由式(3)計(jì)算14.6°N處的緯差加權(quán)法得到的推算結(jié)果。選取6月3日的數(shù)據(jù)成圖，成圖前統(tǒng)一各組日變數(shù)據(jù)基值，結(jié)果如圖3所示。

從圖3可以看出，回歸分析法推算數(shù)據(jù)在數(shù)值上更接近海底日變數(shù)據(jù)，而海底日變數(shù)據(jù)因受到深度及海水過濾等因素的影響，與兩種推算數(shù)據(jù)相比幅值稍小，且損失了一些細(xì)節(jié)信息，導(dǎo)致了一定程度的數(shù)據(jù)失真，而推算數(shù)據(jù)也從一定程度上復(fù)原了這些損失掉的信息。

對(duì)野外采集到的實(shí)測(cè)數(shù)據(jù)做各項(xiàng)改正后，分別用推算得到的兩組日變數(shù)據(jù)和海底日變數(shù)據(jù)對(duì)實(shí)測(cè)數(shù)據(jù)進(jìn)行日變改正，將各交點(diǎn)處的交點(diǎn)差展布于測(cè)線分布圖上，比較后發(fā)現(xiàn)交點(diǎn)差有所減小，統(tǒng)計(jì)30個(gè)交點(diǎn)處(海底日變有部分?jǐn)?shù)據(jù)缺失，改正后只有22個(gè)交點(diǎn))不同數(shù)據(jù)日變改正后的測(cè)線交點(diǎn)差(表5)。

圖3 6月3日推算日變數(shù)據(jù)對(duì)比圖Fig.3 Comparison for speculated diurnal data in June 3rd

交點(diǎn)差最大值/nT交點(diǎn)差均方根差值/nT處理后數(shù)據(jù)交點(diǎn)差均方根差值/nT海底日變24．227．814．25回歸分析法14．244．423．74緯差加權(quán)法15．425．223．92

從表5可以看出，回歸分析法和緯差加權(quán)法推算的日變數(shù)據(jù)改正后交點(diǎn)差均明顯小于海底日變改正后結(jié)果，且回歸分析法改正后效果最好。在測(cè)線分布圖上逐個(gè)交點(diǎn)比較三組數(shù)據(jù)改正后的交點(diǎn)差后，同樣發(fā)現(xiàn)回歸分析法改正后的交點(diǎn)差更小些。

2.3 工區(qū)三計(jì)算結(jié)果

該工區(qū)位于南海北部，調(diào)查測(cè)線在114.5°E～111.0°E、20.5°N～22.5°N范圍內(nèi)布設(shè)，野外采集時(shí)間為3月24日至4月6日，所用磁力儀為Seaspy海洋磁力儀，觀測(cè)精度為0.25 nT。野外測(cè)量時(shí)在工區(qū)附近投放了Sentinel地磁日變觀測(cè)站，3月23日在115.40° E，18.65° N，水深3 798 m處投放，5月6日收回觀測(cè)站，回收的日變數(shù)據(jù)覆蓋整個(gè)測(cè)量日期。

因越南DLT地磁臺(tái)在野外作業(yè)日期內(nèi)無可用數(shù)據(jù)，故將該工區(qū)投放的海底日變站(代號(hào)HD-3)與另三個(gè)陸地臺(tái)站數(shù)據(jù)做相關(guān)性分析，結(jié)果見表6。

表6 工區(qū)三日變數(shù)據(jù)相關(guān)系數(shù)統(tǒng)計(jì)表

從表6可以看出，三個(gè)陸地地磁臺(tái)站與海底日變數(shù)據(jù)相關(guān)性均在0.9以上，在此基礎(chǔ)上考慮到兩個(gè)基臺(tái)站經(jīng)度差越小越好的原因，選擇廣州地磁臺(tái)和海南瓊中地磁臺(tái)這來推算海底日變站緯度(18.65°N)處的地磁日變數(shù)據(jù)。

選取3月24日至4月6日期間海底日變數(shù)據(jù)、廣州地磁臺(tái)數(shù)據(jù)和海南瓊中地磁臺(tái)數(shù)據(jù)作為樣本數(shù)據(jù)，推算出18.65°N處的日變數(shù)據(jù)，選取4月1日的數(shù)據(jù)成圖，因海底日變站投放位置緯度離海南瓊中地磁臺(tái)較近，為了更好地進(jìn)行對(duì)比，將瓊中地磁臺(tái)數(shù)據(jù)一起成圖，成圖前將各組日變數(shù)據(jù)基值進(jìn)行統(tǒng)一，結(jié)果見圖4。

從圖4可以看出，兩組推算數(shù)據(jù)與海南瓊中地磁臺(tái)數(shù)據(jù)較為接近，海底日變數(shù)據(jù)在細(xì)節(jié)信息上的損失較實(shí)例一更多，其幅值較小，與其他三組數(shù)據(jù)的幅值差距較大，這也與海底日變站投放深度較深有關(guān)，深度越深，磁場(chǎng)值的幅值越小，受海水層影響，其細(xì)節(jié)信息損失也越大。

圖4 4月1日推算日變數(shù)據(jù)對(duì)比圖Fig.4 Comparison for speculated diurnal data in April 1st

對(duì)野外采集到的實(shí)測(cè)數(shù)據(jù)做各項(xiàng)改正后，將各交點(diǎn)處的交點(diǎn)差展布于測(cè)線分布圖上，逐點(diǎn)比較交點(diǎn)差，大部分交點(diǎn)差有一定幅度的減小，統(tǒng)計(jì)32個(gè)交點(diǎn)處不同數(shù)據(jù)日變改正后的測(cè)線交點(diǎn)差，結(jié)果見表7。

從圖4和表7中可以看出，推算數(shù)據(jù)改正后結(jié)果，優(yōu)于海底日變和瓊中地磁臺(tái)日變數(shù)據(jù)改正后結(jié)果，其中以回歸分析法改正后數(shù)據(jù)精度為最優(yōu)。相較瓊中地磁臺(tái)的日變數(shù)據(jù)，由海底日變數(shù)據(jù)回歸分析得到的推算數(shù)據(jù)雖然在可靠性上略差，但其實(shí)際上是對(duì)海底日變數(shù)據(jù)補(bǔ)償后的結(jié)果，位置上更接近工區(qū)，因此改正效果要優(yōu)于陸地臺(tái)站。

表7 工區(qū)三數(shù)據(jù)交點(diǎn)差統(tǒng)計(jì)表

3 結(jié)論

通過對(duì)多臺(tái)站地磁日變數(shù)據(jù)做相關(guān)分析，發(fā)現(xiàn)不同緯度處的日變數(shù)據(jù)是具有一定的相關(guān)關(guān)系，緯度差距越小，相關(guān)程度越高。因此可以利用某一處海底日變站數(shù)據(jù)聯(lián)合相近緯度處的兩個(gè)地磁臺(tái)站數(shù)據(jù)，通過統(tǒng)計(jì)學(xué)中的線性回歸分析方法求得海底日變站位置處的日變數(shù)據(jù)，以達(dá)到對(duì)海底日變數(shù)據(jù)進(jìn)行補(bǔ)償?shù)哪康摹?/p>

文中所涉的幾個(gè)工區(qū)投放的海底日變站深度各不相同，導(dǎo)致回收的數(shù)據(jù)與相近緯度處的陸地臺(tái)站數(shù)據(jù)相比存在不同程度的差距，主要體現(xiàn)在細(xì)節(jié)信息丟失和磁場(chǎng)值幅值減小，投放深度越深，這種差距越大，從而影響到日變改正后數(shù)據(jù)的精度。通過分析改正后數(shù)據(jù)交點(diǎn)差的統(tǒng)計(jì)結(jié)果可知，由回歸分析法及緯差加權(quán)法推算得到的數(shù)據(jù)進(jìn)行日變改正后的數(shù)據(jù)精度，均優(yōu)于海底日變數(shù)據(jù)改正后的結(jié)果，且以前者得到的數(shù)據(jù)精度為最優(yōu)。另外通過日變曲線的對(duì)比也可以看出，推算的數(shù)據(jù)一定程度上復(fù)原了因水深和海水層影響等因素造成的海底日變數(shù)據(jù)的細(xì)節(jié)信息丟失和幅值減小，因此可以說明回歸分析法推算日變數(shù)據(jù)可以補(bǔ)償海底日變數(shù)據(jù)，在一定程度上解決了海底日變數(shù)據(jù)日變改正效果不佳的問題，是一種改進(jìn)的海洋磁測(cè)數(shù)據(jù)處理方法。

[1] GB/T13909-92.海洋調(diào)查規(guī)范 海洋地質(zhì)地球物理調(diào)查[S].北京：中國(guó)標(biāo)準(zhǔn)出版社，1993. GB/T13909-92.Specification for oceanographic survey marine geology and geophysics investigation[S]. Beijing: National technology supervised department,1993.(In Chinese)

[2] 管志寧.地磁場(chǎng)與磁力勘探[M].北京：地質(zhì)出版社，2005. GUAN Z N. Geomagnetic field and magnetic exploration[M]. Beijing:Geology Press,2005.(In Chinese)

[3] 劉天佑，陳國(guó)新.海磁日變資料的處理[J].石油地球物理勘探，1987，22(4)：454-460. LIU T Y, CHEN G X. The processing of marine magnetic diurnal variation data[J].Oil Geophysical Explorition ,1987,22(4):454-460. (In Chinese)

[4] 卞光浪，翟國(guó)君，劉雁春，等.海洋磁力測(cè)量中地磁日變站有效控制范圍確定[J].地球物理學(xué)進(jìn)展，2010,25(3)：817-822. BIAN G L, ZHAI G J, LIU Y C, et al. Effective operation range of base stations in marine geomagnetic survey[J]. Progress in Geophys,2010,25(3):817-822. (In Chinese)

[5] 徐行，廖開訓(xùn)，陳邦彥，等.多臺(tái)站地磁日變觀測(cè)數(shù)據(jù)對(duì)遠(yuǎn)海磁測(cè)經(jīng)度的影響分析[J].海洋測(cè)繪，2007,27(1)：38-40. XU X, LIAO K X, CHEN B Y, et al. The effect analysis of the observed diural magnetic variation from multi-stations on the accuracy of marine magnetic survey[J]. Hydrographic surveying and charting,2007,27(1):38-40. (In Chinese)

[6] 單汝儉，金國(guó)，曾志成.局部地區(qū)地磁日變及擬合方法研究[J].長(zhǎng)春地質(zhì)學(xué)院學(xué)報(bào)，1990,(3)：315-322. SHAN R J, JIN G, ZENG Z C. Study on geomagnetic diurnal variation and its fitting methods in local area[J]. Journal of Changchun University of earth science,1990,(3):315-322. (In Chinese)

[7] 邊剛，劉雁春，卞光浪，等.海洋磁力測(cè)量中多站地磁日變改正值計(jì)算方法研究[J].地球物理學(xué)報(bào)，2009,52(10)：2613-2618. BIAN G,LIU Y C,BIAN G L,et al.Reserch on computation method of multi-station diurnal variation corrpction in marine magnetic surveys[J].Chinese Geophys,2009,52(10):1613-2618.(In Chinese)

[8] 卞光浪，劉雁春，翟國(guó)君，等.基于緯差加權(quán)法的海洋磁力測(cè)量多站地磁日變改正值計(jì)算[J].測(cè)繪科學(xué)，2010,35(3):118-120. BIAN G L, LIU Y C, ZHAI G J, et al. Based on the difference of latitude weighting method of Marine magnetic survey multistation geomagnetic divrnal variation correction value is calcalateol[J].Surveying and mapping scienec,2010,355(3):118-120.(In Chinese)

[9] 張向宇，關(guān)永賢，張錫林.回歸分析法在日變數(shù)據(jù)推算中的應(yīng)用[J].物探與化探，2016,40(3)：603-608. ZHANG X Y, GUAN Y X, ZHANG X L. The application os regression to estimating geomagnetic data[J]. Geophusical and Geochemical Exploration, 2016,40(3):603-608. (In Chinese)

[10]時(shí)立文.SPSS 19.0統(tǒng)計(jì)分析從入門到精通[M].北京：清華大學(xué)出版社，2012. SHI L W. SPSS 19.0 statistics analysis from threshold to mastery[M]. Beijing: Qing hua University Press,2012. (In Chinese)

Compensating the data of geomagnetism observation on the seafloor by regression method

ZHANG Xiangyu, GUAN Yongxian, CUI Xiuyun

(Key Laboratory of Marine Mineral Resources of Ministry of Land and Resources, Guangzhou Marine Geological Survey, Guangzhou 510075, China)

We always release geomagnetism observation on the seafloor when the marine magnetic survey is carrying on, but sometimes the recycling data may be fuzzy because of the filter by seawater, and then it will affect the precision. We found that diurnal data from different observation have linear dependence according to the latitude by correlation analysis. We then import regression method to speculate geomagnetic data at certain latitude. We knew that regression method can compensate the data of geomagnetism observation on the seafloor in a certain extent, and improve the precision finally.

marine magnetic survey; geomagnetism observation; regression; diurnal correction

2016-06-21 改回日期：2016-08-06

“127”國(guó)家專項(xiàng)(DD20160227)

張向宇(1987-)，女，碩士，工程師，主要從事海洋重磁數(shù)據(jù)處理技術(shù)研究，E-mail：zhang5821421@163.com 。

1001-1749(2017)03-0327-06

P 631.2

A

10.3969/j.issn.1001-1749.2017.03.05