南極冰下地形測繪路線規劃及制圖研究

王連仲, 孫 波, 吳文會, 郭景學, 崔祥斌, 鄭福海

(1.黑龍江測繪地理信息局極地測繪工程中心，黑龍江 哈爾濱 150086；2.中國極地研究中心，上海 200000；3.黑龍江地理信息工程院，黑龍江 哈爾濱 150086)

南極冰下地形測繪路線規劃及制圖研究

王連仲1, 孫 波2, 吳文會1, 郭景學2, 崔祥斌2, 鄭福海3

(1.黑龍江測繪地理信息局極地測繪工程中心，黑龍江 哈爾濱 150086；2.中國極地研究中心，上海 200000；3.黑龍江地理信息工程院，黑龍江 哈爾濱 150086)

南極冰下地形測繪在國際上還沒有全面開展，獲取南極冰下地形為研究南極冰蓋運動機理，了解南極冰下全貌具有重要意義。本文通過對南極Dome A核心區域冰下地形測繪實踐，并對冰下地形圖成圖方法進行闡述，得出適合現階段南極冰下地形測繪的最優方法——深冰雷達探測法。通過嚴密的規劃測線，可以勘測不同比例尺和精度的冰下地形圖，直觀地獲取冰下形態，為研究南極冰蓋運動機理提供重要依據進行。

南極；冰下地形；測繪；制圖

南極大陸總面積約1 400萬km2，98%的陸地常年被冰雪覆蓋，冰雪平均厚度約為2 000 m，最大厚度約為4 800 m，全球 90%的冰雪儲存在這里。其中Dome A區域是南極海拔最高區域，我國的昆侖站建在這個區域(南緯80°25′01″，東經77° 06′58″，橢球高4 087 m)，自從我國建立昆侖站以來，已經成為南極冰蓋科學考察的重要地區。

南極冰蓋下地形地貌，一直以來就是世界各國科學家想要了解的重要領域，但是，由于南極冰下陸地地形被冰覆蓋，了解冰蓋下的地形地貌狀況是十分困難的。科學家們根據現有的技術條件，研究和利用各種方法和手段來測繪冰下地形，但都限于南極的特殊環境影響勘測結果不夠理想。本文總結出一種測繪冰下地形圖的路線方法，經過實地查驗，在南極昆侖站20 km×20 km的區域對冰下甘比采夫(Gamburtsev)山脈成功實施，測繪出世界上第一幅1∶500 00比例尺精度的冰下地形圖和5 m精度的DEM，為今后南極冰下地形圖測繪，研究南極冰蓋運動機理提供技術支持。

1 冰下地形測繪方法

確定冰下地形測繪采用機載雷達法、引力傳感器測深法、地震勘探法和地面冰雷達測量法。

1)機載雷達和引力傳感器法。2009年2月，來自澳大利亞、英國、加拿大、中國、德國、日本和美國的專家組利用測繪飛機使用機載雷達和引力傳感器對甘比采夫冰下山脈首次進行了勘測，發現了一處驚人的古山脈，位于南極冰下4 km深處，形狀和大小都類似于阿爾卑斯山脈，并繪制出其首批地形圖，如圖1所示。

此方法測繪面積大，但費用高，同時受飛機的不穩定的影響，測量精度不高，對于局部地區詳細冰下地形狀況無法細致的了解。

圖1 使用機載雷達配合引力傳感器測繪的甘比采夫區域地形圖

2)地震勘探法。尋找地下資源可采用的地震勘探法，其原理是利用線纜連接多個傳感器，采用炸藥或空氣槍作為震源，通過人工激發地震波，根據地下不同地質構造，傳感器接收波的傳播信號不同而區分地層介質[2]，如圖2所示。此方法可用于南極冰下地形測繪，但輪式工程車無法進入內陸，需要改裝到雪地車上。

圖2 地震勘探原理

3)冰雷達測深法。2008年1月，在我國第24次南極期間，我國極地研究中心有關專家首次完成了南極Dome A區域200 km×30 km冰下地形測量，受現場條件限制，線路從最窄處幾十米到最寬處十幾公里不等，以此繪制了冰下三維地形圖，見圖3，其成果在Nature雜志上發表[1]。

2010年1月，武漢大學使用簡易冰雷達在南極格羅夫山成功測繪了深度約1 200 m的冰下地形[3]。

綜合分析，采用地面冰雷達測量冰下地形的方法是比較科學和適用的，為此，根據此方法測繪了Dome A昆侖站區域的冰下地形圖。

2 數據獲取及處理

1)深冰雷達。測量用冰深雷達是由中國極地研究中心自主研發的、穿透能力超過3 500 m、擁有自主知識產權的深冰雷達探測系統。它工作于調頻脈沖壓縮體制，工作頻率為500 MHz，發射機最大發射功率為100 W，天線采用對數周期天線，系統的主要技術指標如表1所示。

圖3 冰雷達冰下地形探測圖

表1 冰蓋深部冰雷達系統主要技術指標

2012年12月10—16日完成了所需雷達設備的安裝與調試工作，首先使用的是深部雷達與SIR雷達兩套系統。深部雷達的對數周期天線共分接收與發射天線兩支。首先進行雷達天線自身的組合安裝，每個天線有一個饋電用的集合線和13對陣子，每對陣子有兩支，從集合線上的饋電端口看去，振子的長度逐漸減小；將13對振子通過專用螺絲安裝于集合線上的13對安裝孔。組裝好的天線長為3.35 m，根據系統特性，將兩支天線通過橫梁支架(橫梁支架固定于雷達艙頂部)分別架設于雷達工作艙的左右兩側，天線底端距離冰面約2.3 m，為了天線在行駛過程中的穩固，在橫梁支架的兩端裝配上下斜撐各兩根，天線的連接電纜經雷達工作艙兩側的通線孔進入到艙內。用于淺層探測的SIR雷達天線置于雷達艙后面，根據SIR雷達的天線特性，將天線與冰面進行偶合，將天線固定的特制玻璃鋼雪撬內，并將該小雪撬拖掛在大雪撬的雪撬板正后方，使得玻璃鋼雪撬在行進過程中一直位于大雪撬的撬轍內，以保證淺層雷達天線緊貼冰面平穩的向前滑行。淺層SIR雷達天線的連接電纜通過雷達艙后面的通線也進行艙內，見圖4。

艙內部的深部雷達系統由接收機、發射機、上位機(筆記本)三部分組成。它工作于調頻脈沖壓縮體制，工作頻率是100～200 MHz ，冰體厚度、冰下地形探測深度：>3 500 m，冰體內部結構分辨率：2 m。

2)數據獲取。為了詳細了解昆侖站區域冰下地形情況，2013年1月我國第29次南極內陸考察期間，重新對以昆侖站為中心20 km×20 km范圍冰下地形測量進行嚴密規劃，在已有的1∶500 00比例地形圖上，按相同比例尺對冰下地形測繪進行規劃設計，根據現場情況分4個測量區域，規劃每條航線的起終點坐標，同時，區域間要保證有效連接，考慮到雪地車每天的行走距離，規劃每個區域線路總長不大于100 km，線路間距1.5 km，保證地形圖每個網格的高程點數不少于9個(地形圖測繪基本要求)。

圖4 雷達工作艙冰面測量

利用深冰雷達經過4 d的測量，共采集冰深點3萬余個，覆蓋了基本以昆侖站為中心的附近20 km×20 km的范圍，實際冰深點采集路線如圖5所示。

圖5 實際冰深數據采集路線

對3萬余個冰深數據進行了分析處理，首先，按照國家地形圖規范要求，在每個成圖單元格內(10 cm×10 cm)選取9個以上最深和最淺冰深點，剔除高低連續的中間點，并通過與臨格比較，獲得了有效冰深點數據，編制了此區域的1∶500 00比例尺冰深圖和5 mDEM,同時獲得了南極冰蓋最厚處位置為距離昆侖站東0.7 km處。

為了得到冰下地形點數據，根據同區域的冰面地形圖，內插獲得冰下同平面坐標的冰上高程，計算出冰下地形點高程：

式中：H底為冰下地形點高程；H面為冰面地形點高程；h為此處的冰蓋厚度，即冰層深度。

其中，冰面點高程從1∶500 00比例尺DOME A冰面地形圖獲得。

經過計算，繪制出昆侖站區域1∶500 00比例尺冰下地形圖和DEM為5 m精度的冰下三維地形圖，見圖6。

圖6 1∶5萬Dome A冰下甘比采夫山脈昆侖站區域三維效果圖

3 冰下地形圖精度

通過與文獻1中同區域測量成果對比，可以看出，此區域冰下甘比采夫山脈地形兩次測量在走向和山體形狀上基本一致，個別區域有差別，主要原因是兩次測量點密度不同造成的。通過同區域成圖對比，采用地面冰雷達測量獲取冰下地形圖的方法是比較先進的，如果進行精密規劃，其成圖精度可以滿足1∶1萬～1∶5萬比例尺精度要求，其前提條件是，對測量區域按不同比例尺的規劃要有合理的路線設計，保證獲取數據有效。

4 結 論

采用地面冰雷達測量方法測繪冰下地形圖是現階段最實用、最方便和精度最高的測量方法。在內業成圖處理上，需要計算冰下地形點高程，同時需要對大量的數據進行比較篩選，獲得比較真實的冰下高程，從而繪制出冰下地形。

由于在冰下測量時是按照規劃路線進行，受規劃線路疏密及無法看見冰下地形情況的影響，對冰下地形最高點和最低點可能漏測，不能完全符合冰下地貌，但基本走向和地貌形態不會因此改變。同時，本項工作的開展在研究南極冰蓋運動機理及全球環境變化奠定基礎。

[1] SUN Bo, etc.《The Gamburtsev mountains and the origin and rarly evolution of the Antarctic Ice Sheet》,Nature,Vol.459,4June 2009.

[2] 王澤民，譚智，艾松濤，等.南極格羅夫山核心區冰下地形測繪[J].極地研究，2014，26(4)：399-404.

[3] 王連仲,毛偉娣,陳華.南極地圖符號的標準化探討及在GIS中的實現[J].測繪與空間地理信息，2009，32(5)：68-70.

[4] 韓惠軍,吳文會.南極測量控制點新標志的設計與實現[J].測繪與空間地理信息，2013，36(4)：192-194.

[5] 徐敏,初啟鳳,陽俊,等.南極1∶5000地形圖測繪關鍵技術及相關問題處理——以南極維多利亞地站為例[J].測繪工程，2016，24(1):73-76.

[6] 馮守珍，雷寧.南極中山站海灣地形特征研究[J].海洋測繪，2005，25(5)：28-30.

[7] 吳文會，殷福忠，吳迪.基于直升機平臺及非量測數碼相機的南極航空攝影測量技術研究[J].極地研究，2010，22(2)：190-198.

[8] 王連仲，孔繁宇，朱虹.WalkField軟件在南極測繪中的應用[J].測繪工程，2007，15(3):1-5.

[9] 吳學峰，魏寶安，喬輝.高精度GPS測量在南極內陸測繪中的應用研究[J].測繪與空間地理信息，2010，33(4)：89-91.

[10] 潘家訓.南極麥克斯威爾灣的測繪及海底地形分析[J].海洋測繪，1995，15(2)：30-37.

[11] 溫家洪，王清華，吳健平.地理信息系統在南極考察和研究中的應用[J].極地研究，2001，13(3)：217-228.

[責任編輯：李銘娜]

Design and research of subglacial topography surveying and mapping in Antarctic

WANG Lianzhong1, SUN bo2, WU Wenhui1, GUO Jingxue2, CUI Xiangbin2，Zheng Fuhai3

(1.Heilongjiang Bureau of Surveying and Mapping,Harbin 150086，China;2.China′s Institute of Polar Research,Shanghai 200000,China;3.Heilongjiang Institute of Geomatics Engineering, Harbin 150086,China)

Internationally,subglacial topography surveying and mapping has not been comprehensively carried out in Antarctic,which is the key of research on the ice sheet motion to obtain the whole subglacial topography.This paper gives the optimal method at the present stage through summarizing subglacial topography surveying and mapping in Dome A area,Antarctic.At the same time,a simple method of subglacial topography mapping is provided.

Antarctic;subglacial topography;surveying;mapping

引用著錄：王連仲,孫波,吳文會，等.南極冰下地形測繪路線規劃及制圖研究[J].測繪工程，2017，26(5)：15-19.

10.19349/j.cnki.issn1006-7949.2017.05.004

2016-09-27

王連仲(1963-)，男，高級工程師.

P715

A

1006-7949(2017)05-0015-05