衛(wèi)星測高反演重力場在西太平洋海域的應用?

孫 楊， 姜效典， 王 燚， 李朝陽(中國海洋大學 .海洋地球科學學院；.計算機基礎部，山東 青島 6600)

衛(wèi)星測高反演重力場在西太平洋海域的應用?

孫 楊1， 姜效典1， 王 燚2， 李朝陽1
(中國海洋大學 1.海洋地球科學學院；2.計算機基礎部，山東 青島 266100)

利用高精度衛(wèi)星測高數(shù)據(jù)反演重力是獲取海洋重力異常的重要方法。本文介紹了自2008年以來新一代高精度測高衛(wèi)星的相關信息。在此基礎上，聯(lián)合多顆新型衛(wèi)星，進行交叉點平差，通過中心差分代替常規(guī)Sandwell方法計算垂線偏差，利用改進的逆Vening-Meinesz法和Stokes數(shù)值反解法對西太平洋某區(qū)域進行了重力場反演。反演結果與NGDC提供的船測重力值比較表明，逆Vening-Meinesz法計算誤差小于Stokes數(shù)值反解法。通過添加隨機噪聲數(shù)值實驗，進一步證明逆Vening-Meinesz法對于測高數(shù)據(jù)誤差的抗差能力更強，在研究區(qū)的重力異常計算結果可信度較高。

測高衛(wèi)星；海洋重力；逆Vening-Meinesz；垂線偏差

衛(wèi)星測高技術作為大地測量的重要手段，具有全天候，高重復率，高精度，低成本的特點。近年來隨著測高數(shù)據(jù)精度的不斷提高，其在地球重力場推算，尤其是快速推求大區(qū)域重力數(shù)據(jù)過程中的突出優(yōu)勢越來越引起人們的重視。與船測重力異常相比，衛(wèi)星反演重力數(shù)據(jù)穩(wěn)定性好，所有數(shù)據(jù)都基于同一基準面觀測，基本不受參考框架統(tǒng)一性問題的限制，避免了不同船測區(qū)域間需要數(shù)據(jù)融合的困擾。國內(nèi)外眾多學者對聯(lián)合多顆衛(wèi)星測高數(shù)據(jù)進行重力場反演開展了大量工作，萬劍華利用Cryosat-2衛(wèi)星數(shù)據(jù)在南海進行了重力異常反演，分辨率達到了2’×2’，精度達到了4.5×10-5m/s2[1]。彭富清，李洋等分析了衛(wèi)星測高誤差對海洋重力場反演精度的影響，指出反演誤差與測高隨機誤差呈正比變化，與空間分辨率成反比變化[2-3]。Sandwell聯(lián)合Cryosat-2，Envisat和Jason-1衛(wèi)星在墨西哥灣進行了重力反演，精度達到了2×10-5m/s2[4]。在此基礎上，建立了新一代的海洋重力模型，分析了目前衛(wèi)星重力反演存在的主要問題，指出了未來的研究方向[5]。此外，Hwang[6]，Keating[7]和Garcia[8]等也先后開展了針對淺水海域衛(wèi)星重力反演的研究，并取得了卓有成效的結果。鑒于針對新型測高衛(wèi)星數(shù)據(jù)聯(lián)合反演的系統(tǒng)研究較少，本文以近年來多顆新測高衛(wèi)星為研究對象，對西太平洋某海域進行了重力場的計算，并與NGDC船測數(shù)據(jù)對比分析誤差結果。

1 數(shù)據(jù)來源與研究區(qū)概況

利用衛(wèi)星測高數(shù)據(jù)反演得到的重力場精度主要依賴于四個因素，即軌道的空間密度，測高計高程精度，異軌數(shù)據(jù)的方向(方位角)，以及近岸涌浪模型的精度。近幾年，多顆新型測高衛(wèi)星進入太空，同時有已服役衛(wèi)星進行了軌道調(diào)整，大量新觀測數(shù)據(jù)的獲得,為進行高精度，高分辨率的衛(wèi)星重力反演提供了可能[3]。下表列出了這些衛(wèi)星的主要參數(shù)。

其中，Envisat New和Jason-1 Geodetic分別在執(zhí)行了8和11年任務之后實施了變軌操作，而Jason-2，Cryosat-2，Altika和HY-2則發(fā)射于2008—2014的6年間。

研究區(qū)位于西太平洋帕里西維拉海盆和馬里亞納弧盆系統(tǒng)的接合部(140.5°E～143.5°E，15°N～21.5°N)，西部為帕里西維拉海盆，水深較深，東部為西馬里亞納島弧，水深較淺。區(qū)中有2006年實測的重力觀測數(shù)據(jù)，由NGDC公布，編號BD0601，共35條主測線，3條聯(lián)絡測線，214 189個采樣點。船測重力異常最小值為-13.29×10-5m/s2，最大為148.55×10-5m/s2，平均值為22.8×10-5m/s2。

多星聯(lián)合進行重力計算是提高分辨率的重要手段。研究區(qū)不同時期衛(wèi)星地面軌跡如圖2所示。從圖中可以看出，不同時期衛(wèi)星軌跡密度有較大差別，在一定程度上會影響最終重力反演的結果。

2 重力異常計算

不同衛(wèi)星之間存在參考框架，測高精度以及各種校正模型的系統(tǒng)偏差，在進行重力反演之前需要先進行交點平差計算，減小系統(tǒng)誤差。之后利用測高數(shù)據(jù)和海面地形數(shù)據(jù)計算大地水準面高，從而可以得到垂線偏差數(shù)據(jù)。在此基礎上，利用逆Vening-Meinesz公式和“移去-恢復”技術即可反演出研究區(qū)的自由空間重力異常。下面介紹具體的計算過程。

2.1 交叉點平差

交點平差可以減弱衛(wèi)星軌道徑向誤差造成的影響[14-16]。若將交叉點海面高的不符值作為觀測量，采用線性模型模擬軌道誤差，可建立平差觀測方程

(1)

(2)

2.2 大地水準面確定

將得到的平差后的海面高數(shù)據(jù)利用Kriging法進行網(wǎng)格化，然后減去由AVISO提供的平均海平面MSS_CNES_CLS11模型校正后的動態(tài)海面地形網(wǎng)格數(shù)據(jù)，即可得到大地水準面高。下圖給出了西太平洋地區(qū)10′×10′的大地水準面分布。

2.3 計算垂線偏差

Sandwell和Olgiati等先后提出了不同的垂線偏差計算方法[17-19]。兩種方法得到的數(shù)據(jù)不規(guī)則地分布在觀測點和交叉點上，須額外將其處理成網(wǎng)格點的平均值。Sandwell方法利用了星下點速度參量，Olgiati方法則需要內(nèi)插垂直于軌跡的垂線偏差，這些步驟都引入了不必要的數(shù)值誤差。此外，褚坤曾嘗試利用大地水準面的網(wǎng)格數(shù)據(jù)直接計算垂線偏差[8]。

為了建立快速，準確的計算方法，本文從直觀物理意義出發(fā)，利用大地水準面空間差分形式計算垂線偏差。網(wǎng)格點(i,j)處的垂線偏差子午和卯酉分量分別為

(3)

(4)

式中：ξ為垂線偏差子午分量；η為垂線偏差卯酉分量；N為大地水準面高；R為地球平均半徑；φ為大地緯度；λ為大地經(jīng)度。在邊界處，可以相應的用前向差分和后向差分代替中心差分。通過這種方法，可以快速的進行垂線偏差的計算。圖4為利用差分法計算得到的10′×10′垂線偏差子午和卯酉分量。

2.4 計算重力異常

根據(jù)衛(wèi)星測高數(shù)據(jù)計算重力異常目前常用的有逆Vening-Meinesz方法和Stokes數(shù)值反解法。本文主要利用改進的逆Vening-Meinesz方法計算重力異常。其二維球面譜計算式為[22]

(5)

其中：H′(φ，λ)為積分核函數(shù)H(φ，λ)對球面角距φ的導數(shù)；α為積分點到計算點之間的方位角。當Ψ=0時，核函數(shù)H′(Ψ)將發(fā)生奇異。應當單獨考慮計算點所在方塊數(shù)據(jù)的影響，并對積分核函數(shù)做出相應的改進。其內(nèi)環(huán)帶的影響量可以表示為

(6)

其中：S0為內(nèi)環(huán)帶的球冠半徑；ξx和ηy分別表示ξ和η在x和y方向上的變化率。

與逆Vening-Meinesz法不同，Stokes反解法直接將測高數(shù)據(jù)作為輸入量，其二維譜計算式為

(7)

其中：N為大地水準面高；S(φ，λ)為Stokes函數(shù)。

2.5 移去-恢復技術

在計算過程中，還需要采用“移去-恢復”技術，消除長波長數(shù)據(jù)造成的影響[4]。通過EGM2008模型球諧系數(shù)計算模型垂線偏差[20-24]。從2.3步得到的垂線偏差中減去模型垂線偏差，即為殘余垂線偏差，進而計算得出殘余重力異常。殘余重力異常與EGM2008模型重力異常疊加，即為最終重力異常[17]。

3 實際資料處理

采用逆Vening-Meinesz法，聯(lián)合圖2右圖所示的Jason-2，Cryosat-2，Altika和HY-2四顆衛(wèi)星對研究區(qū)進行自由空間重力異常的反演。結果如下圖所示。

從圖中可以看出，研究區(qū)以正異常為主，正異常高值集中在東側(cè)西馬里亞納島弧處，最大值為130.5×10-5m/s2。西馬里亞納島弧西側(cè)異常變化劇烈，沿島弧分布有南北向的異常梯度帶，顯示此處地殼結構變化劇烈。最小值位于研究區(qū)中部， 為-11.6×10-5m/s2。進一步選取圖5中的AA’和BB’兩條測線，并利用Stokes數(shù)值反解法計算測線上的重力異常，最后將船測重力異常與對應點兩種方法的反演結果進行對比。如下圖所示，圖中縱坐標為自由空間重力異常，橫坐標分別為測線上點距A點或B點的水平距離。

從圖中可以看出，兩種方法得到的衛(wèi)星重力異常與船測重力異常波動趨勢一致，但都較為平滑。測高衛(wèi)星數(shù)據(jù)的主要信息集中在中長波段，對于波長小于10km的短波特征的表現(xiàn)較弱。整體來看，逆Vening-Meinesz法結果與實測重力更接近。而且對于細節(jié)的反映要好于Stokes數(shù)值反解法。具體到測線上，AA’測線上兩種方法的反演結果要優(yōu)于BB’測線。在BB’測線上，Stokes數(shù)值反解的結果已經(jīng)與實測重力有較大的誤差。對研究區(qū)和兩條測線的數(shù)據(jù)分別統(tǒng)計誤差信息如下表所示。

誤差統(tǒng)計表明，對于整個研究區(qū)，逆Vening-Meinesz法精度在5.77×10-5m/s2左右。兩種方法在AA’測線的結果要好于BB’測線的結果。這主要由兩個原因造成的，一是兩種方法在重力數(shù)據(jù)變化劇烈的區(qū)域反解精度要低于數(shù)值變化平緩的區(qū)域。另一方面，BB’測線附近衛(wèi)星測高數(shù)據(jù)精度較差。由于此處水深較淺，衛(wèi)星測高數(shù)據(jù)受回波影響精度大大降低，同時一些校正量如大氣濕對流層校正、海洋潮汐校正以及大氣高頻因數(shù)校正在淺海不準確，也會對結果造成影響。AA’測線處水深在4 500m左右，衛(wèi)星數(shù)據(jù)質(zhì)量高，逆Vening-Meinesz法誤差約為1.45×10-5m/s2?？紤]到船測重力也存在(1～3)×10-5m/s2左右的偏差，利用衛(wèi)星測高數(shù)據(jù)計算得到的重力異常在該區(qū)域已經(jīng)接近于普通科考船的數(shù)據(jù)精度。

通過兩種方法間誤差對比可以看出，逆Vening-Meinesz法的精度要高于Stokes數(shù)值反解法。在AA’測線處逆Vening-Meinesz法精度要高出2.76×10-5m/s2，在BB’測線處精度更是高出4.26×10-5m/s2，說明其對于測高數(shù)據(jù)的抗差能力更強。為進一步論證該結論，將BB’測線的測高數(shù)據(jù)添加平均值為0，均方差為0.05m的隨機誤差。分別計算兩種方法受隨機誤差的影響。結果如圖7所示。

從圖7中可以看出，逆Vening-Meinesz法和Stokes數(shù)值反解法都一定程度地受到噪聲的影響。噪聲數(shù)據(jù)通過逆Vening-Meinesz法產(chǎn)生誤差的標準差為8.263 4×10-5m/s2，Stokes數(shù)值反解法的誤差則為17.330 6×10-5m/s2?？梢?，逆Vening-Meinesz法的穩(wěn)定性要優(yōu)于Stokes數(shù)值反解法。

研究區(qū)內(nèi)既存在重力場劇烈變化的區(qū)域，也包含了淺水區(qū)。這些不利影響對利用衛(wèi)星測高數(shù)據(jù)反解重力場造成了很大的困難。逆Vening-Meinesz法可以有效地減弱測高誤差對于結果的影響，同時劇烈變化的重力場對其造成的干擾也相對較小。因此，利用逆Vening-Meinesz法得到的研究區(qū)重力異常數(shù)據(jù)具有較高的精度，可以很好地反映區(qū)域內(nèi)重力異常的空間分布。

4 結語

本文在介紹近年來新發(fā)射的多顆測高衛(wèi)星的基礎上，以改進的逆Vening-Meinesz方法為基礎，引入一階中間差分算法代替Sandwell垂線偏差法，計算了西太平洋某海區(qū)重力異常。同時利用Stokes數(shù)值反解法計算的測線結果和NGDC提供的船測實測數(shù)據(jù)與其進行了對比。結果顯示逆Vening-Meinesz法受重力場變化劇烈程度影響較小，計算結果與實測重力異常的符合度更高，很好地反映了研究區(qū)內(nèi)中長波重力異常分布。

數(shù)值實驗進一步證明，相對于Stokes數(shù)值反解對測高數(shù)據(jù)的抗差能力更強。由于目標區(qū)域內(nèi)包含淺水海域，衛(wèi)星測高數(shù)據(jù)誤差增大，故利用逆Vening-Meinesz法得到的重力場分布更加準確。

衛(wèi)星測高數(shù)據(jù)反演在獲取大區(qū)域海洋重力場分布中正起到越來越重要的作用。針對重力場變化劇烈和淺水海域等特殊研究區(qū)，逆Vening-Meinesz法可以提供更好的重力異常數(shù)據(jù)。

致謝：感謝法國國家空間研究中心(CNES)提供衛(wèi)星數(shù)據(jù)。

[1] 萬劍華, 李家軍, 劉善偉, 等. 基于Cryosat-2 數(shù)據(jù)的南海 2'× 2'重力異常計算與分析[J]. 中國石油大學學報, 2015, 39(3): 70-75.WanJH,LiJJ,LiuSW,etal.Calculationandanalysisof2′×2′gravityanomaliesovertheSouthChinaSeabasedonCryosat-2satellitealtimeterdata[J].JournalofChinaUniversityofPetroleum, 2015, 39(3): 70-75.

[2] 彭富清, 陳雙軍, 金群峰. 衛(wèi)星測高誤差對海洋重力場反演的影響[J]. 測繪學報, 2014, 43(4): 337-340.PengFQ,ChenSJ,JinQF.Influenceofaltimetryerrorsonmarinegeopotentialrecovery[J].ActaGeodaeticaetCartographicaSinica, 2014, 43(4): 337-340.

[3] 李洋, 張潤寧. 高度計測距精度對沿軌跡重力異常反演的影響[J]. 測繪學報, 2015, 44(4):G363-G369.LiY,ZhangRN.Inferenceofaltimeteraccuracyonalong-trackgravityanomalyrecovery[J].ActaGeodaeticaetCartographicaSinica, 2015, 44(4):G363-G369.

[4]SandwellD,GarciaE,SoofiK,etal.Toward1-mGalaccuracyinglobalmarinegravityfromCryoSat-2,Envisat,andJason-1[J].TheLeadingEdge, 2013, 32(8): 892-899.

[5]SandwellDT,MullerRD,SmithWHF,etal.NewglobalmarinegravitymodelfromCryoSat-2andJason-1revealsburiedtectonicstructure[J].Science, 2014, 346(6205): 65-67.

[6]HwangC,HsuHY.Shallow-watergravityanomaliesfromsatellitealtimetry:CasestudiesintheeastchinaseaandTaiwanstrait[J].JournaloftheChineseInstituteofEngineers, 2008, 31(5): 841-851.

[7]KeatingP,PinetN.Comparisonofsurfaceandshipbornegravitydatawithsatellite-altimetergravitydatainHudsonBay[J].TheLeadingEdge, 2013, 32(4): 450-458.

[8]GarciaES,SandwellDT,SmithWHF.RetrackingCryoSat-2,EnvisatandJason-1radaraltimetrywaveformsforimprovedgravityfieldrecovery[J].GeophysicalJournalInternational, 2014:ggt469.

[9]JiangX,LinM,LiuJ,etal.TheHY-2satelliteanditspreliminaryassessment[J].InternationalJournalofDigitalEarth, 2012, 5(3): 266-281.

[10] 彭海龍, 林明森, 穆博, 等.HY-2A衛(wèi)星雷達高度計數(shù)據(jù)的全球統(tǒng)計評價及質(zhì)量分析[J]. 海洋學報, 2015(7): 54-66.PengHL,LinMS,MuB,etal.HY-2AGlobalstatisticalevaluationandperformanceanalysisofHY-2Asatelliteradaraltimeterdata[J].ActaOceanologicaSinica, 2015(7): 54-66.

[11] 楊磊, 周興華, 彭海龍, 等. 基于Jason-2 的Saral/AltiKa高度計全球統(tǒng)計評估與交叉定標[J]. 海洋科學進展, 2014, 32(4): 482-490.YangL,ZhouXH,PengHL,etal.Globalassessmentandcross-calibrationofSaral/AltikabasedonJason-2Altimeter[J].AdvancesinMarineScience, 2014, 32(4): 482-490.

[12]WinghamDJ,FrancisCR,BakerS,etal.CryoSat:AmissiontodeterminethefluctuationsinEarth’slandandmarineicefields[J].AdvancesinSpaceResearch, 2006, 37(4): 841-871.

[13]ChaudharyA,BasuS,KumarR,etal.RetrievingsignificantwaveheightintheIndianOceannearVisakhapatnamusingJason-2altimeterdata[J].RemoteSensingLetters, 2015, 6(4): 286-294.

[14] 褚坤.Geosat高度計數(shù)據(jù)處理與南海重力異常反演精度評價[D]. 青島: 中國石油大學, 2010.ChuK.GeosatAltimeterDataProcessingandAccuracyAssessmentontheInversionofGravityAnomaliesintheSouthSea[D].Qingdao:ChinaUniversityofPetroleum, 2010.

[15] 王虎彪, 王勇, 陸洋, 等. 用衛(wèi)星測高和船測重力資料聯(lián)合反演海洋重力異常[J]. 大地測量與地球動力學, 2005, 25(1): 81-85.WangHB,WangY,LuY,etal.Inversionofmarinegravityanomaliesbycombinatingmulti-altimeterdataandshipbornegravimetricdata[J].JournalofGeodesyandGeodynamics, 2005, 25(1): 81-85.

[16] 黃謨濤, 翟國君, 歐陽永忠, 等. 利用多代衛(wèi)星測高數(shù)據(jù)反演海洋重力場[J]. 測繪科學, 2006, 31(6): 37-39.HuangMT,ZhaiGJ,OuyangYZ,etal.Recoveryofmarinegravityfieldusingintegrateddatafrommulti-satellitemissions[J].ScienceofSurveyingandMapping, 2006, 31(6): 37-39.

[17] 王莉娟. 多源衛(wèi)星高度計中國近海及其鄰域重力異常反演與精度檢驗[D]. 青島: 中國石油大學, 2012.WangLJ.TheInversionoftheGravityAnomaliesandtheAssessmentofItsAccuracyBasedonMulti-satelliteAltimeterDatainChinaSeaandItsVicinity[D].Qingdao:ChinaUniversityofPetroleum, 2012.

[18]SandwellDT.AdetailedviewoftheSouthPacificgeoidfromsatellitealtimetry[J].JournalofGeophysicalResearch:SolidEarth, 1984, 89(B2): 1089-1104.

[19]OlgiatiA,BalminoG,SarrailhM,etal.Gravityanomaliesfromsatellitealtimetry:comparisonbetweencomputationviageoidheightsandviadeflectionsofthevertical[J].Bulletingéodésique, 1995, 69(4): 252-260.

[20] 黃謨濤. 海洋重力場測定及其應用 [M]. 北京： 測繪出版社, 2005.HuangMT.DeterminationandApplicationofMarineGravityField[M].Beijing：SurveyingandMappingPress, 2005.

[21] 胡淑梅, 文漢江, 李洪超, 等. 利用多種測高數(shù)據(jù)反演中國南海海域重力異常[J]. 大地測量與地球動力學, 2011, 31(4): 56-59.HuSM,WenHJ,LiHC,etal.InversionofgravityanomaliesoverSouthChinaSeabyuseofcombinationofmulti-satellitealtimeterdata[J].JournalofGeodesyandGeodynamics, 2011, 31(4): 56-59.

[22] 孫鳳華, 郭玉良, 陳春旺, 等. 我國沿海及海域 1'× 1'垂線偏差模型的研究與建立[J]. 海洋測繪, 2002, 22(5): 12-15.SunFH,GuoYL,ChenCW,etal.Researchandconstructionof1' × 1'deflectionoftheverticalmodelforChinaseaarea[J].HydrographicSurveyingandCharting, 2002, 22(5): 12-15.

[23] 黃謨濤, 翟國君, 管錚, 等. 利用衛(wèi)星測高數(shù)據(jù)反演海洋重力異常研究[J]. 測繪學報, 2001, 30(2): 179-184.HuangMT,ZhaiGJ,GuanZ,etal.Ontherecoveryofgravityanomaliesfromaltimeterdata[J].ActaGeodaeticaetCartographicaSinica, 2001, 30(2): 179-184.

[24] 楊樂, 林明森, 張有廣,等. 中國近岸海域高度計JASON-1測量數(shù)據(jù)的波形重構算法研究[J]. 海洋學報, 2010, 32(6):91-101.YangL,LinMS,ZhangYG,etal.ImprovingthequalityofJASON-1altimetrydatabywaveformretrackingincoastalwatersoffChina[J].ActaOceanologicaSinica, 2010, 32(6):91-101.

[25]SandwellDT,SmithWHF.RetrackingERS-1altimeterwaveformsforoptimalgravityfieldrecovery[J].GeophysicalJournalInternational, 2005, 163(1): 79-89.

責任編輯 徐 環(huán)

Application of Gravity Inversion From Satellite Altimeter Data in the West Pacific

SUN Yang1，JIANG Xiao-Dian1，WANG Yi2，LI Chao-Yang1

(Ocean University of China 1. College of Marine Geosciences; 2. Department of Comupter Foudation，Qingdao 266100， China)

The application of gravity anomaly inversion from high quality satellite altimeter data is an important method to acquire marine gravity. The main characteristics of newly available satellite altimeter data since 2008 are reviewed in this paper. On the basis of multi-satellite data, we apply crossover adjustment and use central difference method to calculate deflection of the vertical instead of conventional Sandwell method.After that, improved inverse Vening-Meinesz method and numerical inverse-Stokes approach are adopted to calculate the gravity filed of certain area in the West Pacific.Comparision with the shipboard-based gravity provided by NGDC shows that inverse Vening-Meinesz method has smaller errors than numerical inverse-Stokes approach.With further numerical expericment ofadding stochastic noise, we show that inverse Vening-Meinesz method is more robust so that it can provide reliable gravity estimation for the target area.

Satellinte altimeter; marine gravity; inverse Vening-Meinesz; deflection of the vertical

海洋公益性行業(yè)科研專項(201305029-02)資助 Supported by the Public Science and Technology Research Funds Projects of Ocean (201305029-02)

2015-12-07;

2016-03-07

孫楊(1987-)，男，博士。E-mail：sunyang_ouc@163.com

P312.1

A

1672-5174(2017)06-087-08

10.16441/j.cnki.hdxb.20150411

孫楊，姜效典，王燚,等.衛(wèi)星測高反演重力場在西太平洋海域的應用[J].中國海洋大學學報(自然科學版),2017, 47(6): 87-94.

SUN Yang，JIANG Xiao-Dian，WANG Yi，et al.Application of gravity inversion from satellite altimeter data in the West Pacific[J].Periodical of Ocean University of China, 2017, 47(6): 87-94.