稀土礦山地質環境調查中的三維遙感技術研究

劉鳳梅，曾 敏

(中國地質調查局武漢地質調查中心，武漢 430205)

稀土礦山地質環境調查中的三維遙感技術研究

劉鳳梅，曾 敏

(中國地質調查局武漢地質調查中心，武漢 430205)

以遙感和GIS為主要技術手段，對CBERS－02B數據進行幾何糾正、鑲嵌融合和正射糾正處理后，在ArcGIS平臺上，利用1∶1萬地形圖建立數字高程模型(DEM)，通過疊加遙感影像對地形地貌進行三維模擬，制作了江西省安遠縣稀土礦開采區三維遙感影像圖，并對礦區地質環境進行了三維分析實驗。實踐證明，三維遙感技術對野外地質調查困難區域的野外調查路線選擇、典型區域勘察以及指導礦山生產等都有重要的現實意義;且能快速而準確地揭示礦山地質環境的空間分布規律，計算各類地質環境要素的面積、體積;還可以疊加地質環境要素等矢量數據，通過緩沖區分析對地質災害的危害范圍及發展趨勢進行分析和預測，為減災排險和恢復治理提供決策依據。

CBERS－02B;礦山地質環境;三維分析;緩沖區分析

0 引言

遙感圖像的三維技術基于DEM數據和遙感影像，可實時顯示、繪制動態的三維地形地貌景觀。三維地形的數字化顯示不單是提供真實觀感的問題，也是對空間數據及所處環境的一種非常好的描述，是二維現實地形地貌和水系分布延伸到三維地形世界的真實再現，以直接、可視化、易于理解的形式表達二維地表信息，可應用于所有三維實時顯示的場合(特別適合大范圍、不規則地形的三維可視化表達)。對三維遙感技術的研究和應用，不僅為本文進行礦山地質環境調查研究提供了有利工具，也是真實三維影像、視景仿真、虛擬現實、虛擬地理信息系統等領域的前沿和熱點研究技術［1－2］。

本文采用中巴地球資源衛星02B星(CBERS－02B)遙感數據，經CCD與HR數據融合后(空間分辨率達2.36 m)，與利用1∶1萬地形圖建立的DEM相結合，對稀土礦開采區的地形地貌進行三維模擬，快速而準確地查清礦山地質環境要素的基礎信息和空間分布規律。經野外驗證結果可信，并取得為減災排險和恢復治理提供決策依據的社會效應。

1 研究區與數據源

1.1 研究區概況

安遠縣車頭稀土礦區位于江西省贛州市安遠縣城以西 6 km，地理坐標為 E115°11＇57″～115°14＇06″，N25°01＇52″～25°03＇52″，治理礦區面積 16.3 km2。礦區屬亞熱帶溫暖濕潤季風氣候區，降雨強度大、暴雨多，降雨量年際與月際變化大，旱澇災害頻發，尤其突發的山洪常常使江河水位猛漲，沖毀農田，嚴重剝蝕表土層與植被蓋層，造成巨大的生命財產損失;開采過度的礦山極易發生泥石流災害并形成淤積區，毀壞礦區周邊農田。

1.2 數據預處理

1.2.1 數據源

本次研究采用CBERS－02B數據，采取遙感解譯與野外實地調查相結合的手段，進行安遠縣車頭稀土礦山地質環境調查。

采用的CBERS－02B數據為2008年12月13日獲取的CCD數據與2009年12月19日獲取的HR數據。通過融合處理，生成模擬真彩色圖像。

1.2.2 遙感數據預處理

CBERS－02B圖像數據均為二級產品，只經過輻射校正和幾何粗糾正處理，所以對CBERS－02B衛星圖像又進行了配準、幾何精糾正與圖像增強、多波段圖像合成及裁剪等處理。

采用主成分變換法對CBERS－02B衛星CCD傳感器的1～4波段與HR傳感器的高分辨率波段圖像進行融合，生成含上述4個波段信息的融合圖像。根據雪氏熵值法統計分析和地物波譜特征，選取CBERS－02B 融合圖像的 Band 4(R)、Band 3(G)、Band 2(B)組合制作的假彩色合成圖像為最佳解譯波段組合的圖像;以 Band 3(R)、Band 2(G)、Band 1(B)組合制作的模擬真彩色圖像適用于專題制圖。

采用1∶1萬地形圖對融合后的CBERS－02B數據進行幾何精糾正。糾正過程中盡量選取圖像中明顯、清晰的永久性定位識別標志(如道路交叉點、建筑物邊界、橋梁等)作為控制點，使糾正后圖像的幾何誤差控制在1個像元之內;從糾正后的圖像中裁剪出研究區范圍進行信息提取。

2 礦山地質環境遙感解譯

2.1 遙感解譯方法

遙感解譯基于地學原理對基礎圖像上的地物進行判讀和識別，并對其進行定性、定量、時間和空間分析，獲取各類地物基本屬性特征和影像特征［3］。

根據安遠縣稀土礦區的開采特征，本次解譯的礦山地質環境要素主要包括采剝區、堆積區、尾砂淤積區、農田淤積區和采礦后已恢復區，根據前期調查結果建立了礦山地質環境要素解譯標志(表1)。

表1 解譯標志一覽表Tab.1 List of interpretation keys

根據遙感圖像的空間分辨率情況，本次遙感解譯在ENVI和ArcGIS軟件平臺上采用人機交互的方式進行，提取礦山地質環境要素，構建拓撲關系，并按照實際規模大小在1∶1萬地形圖上進行標注。遙感解譯是一個初步解譯—野外驗證—詳細解譯的綜合、反復過程，具體步驟如圖1所示。

圖1 遙感解譯技術流程Fig.1 Flow chart of RS interpretation

2.2 野外調查驗證

礦山地質環境野外調查的目的主要是驗證遙感解譯圖斑的可靠性和準確性，特別是對采剝區與堆積區界線的劃分、堆積區與尾礦淤積區等問題進行重點調查，對室內解譯結果進行再修正，并對室內解譯與野外驗證結果不一致的原因進行分析。

采取GPS實地測量的方法進行野外調查，并填寫野外調查記錄。對于采剝區、尾砂堆積區及淤積區分界線的地段，使用GPS布置代表性路線進行追蹤驗證［4］。在野外共驗證了23個礦區的50個圖斑，驗證結果表明遙感解譯的判對率高達94.7%以上(表2)。

表2 野外驗證結果與圖斑屬性精度評價Tab.2 Precision evaluation of field verification results and spot attributes

2.3 解譯結果分析

對第一次室內解譯結果進行野外驗證后，共解譯圖斑108個(包括采剝區55個、尾砂堆積區31個、尾砂淤積區17個、農田淤積區3個和采礦后已恢復區2個)，詳見表3。

表3 解譯結果Tab.3 Interpretation results

3 礦區地質環境三維分析

為了更清楚、更直觀地表現稀土礦開采所引起的一系列礦山地質環境問題，利用高分辨率圖像數據和數字高程模型(DEM)制作了研究區的三維顯示效果圖。本文選用ENVI 4.4和ArcGIS 9.1生成有關礦區的三維模型，礦區三維模型效果如圖2所示。

圖2 研究區三維模型效果圖Fig.2 3D model of study area

3.1 DEM模型的建立

建立DEM模型的方法主要有數據內插和三角測量。內插方法包括反距離權重、樣條函數和Kriging插值法［5］。

DEM內插是根據參考點的高程求出其他待定點的高程，原理上雖屬于數學插值問題，但比數學插值更復雜(因為DEM內插必須充分考慮地形的連續性和光滑性，或者說鄰近地面高程點之間有很大的相關性)。因此，不是所有的插值方法都適用于DEM數據的內插。本文從地形的連續性和光滑性這一特點出發，采用Kriging插值法對DEM數據進行插值;使用ArcMap的三維分析模塊，以1∶1萬地形圖數據為基礎，建立區域DEM模型，并生成1∶1萬DEM 暈渲圖(圖3)。

圖3 研究區1∶1萬DEM暈渲圖Fig.3 DEM hill－ shading map of study area(1 ∶10000)

3.2 遙感影像三維模型的建立

以DEM數據為基礎，在ArcGIS中疊加遙感影像及各類矢量數據，構建真實地表三維顯示圖像。可以通過不同的瀏覽器從多個角度瀏覽圖像上的不同場景，有助于識別目標地物的影像特征，分析影像與地形的相關性［6］。

在ArcGIS的ArcScene模塊中，添加DEM、遙感影像和解譯的地質環境矢量數據，場景初始狀態的基準高程值為0。影像中的地面紋理是地形的主要信息來源，將遙感影像、地理要素、解譯矢量數據和文字符號標注等多種數據疊加到DEM模型上;以DEM模型高程為基準，在Scene Properties對話框設置場景中所有圖層的屬性，生成三維地形影像［7］。由于為遙感影像和地質環境要素矢量數據添加了高程信息，地面紋理、各類地質環境要素和地形地貌之間的關系就會顯示出來，建立起地形地貌的三維顯示場景(圖4)。

圖4 礦區地質環境三維解譯與野外調查對比Fig.4 Comparison between 3D interpretation and field survey of mining area

3.3 礦山地質環境三維可視化分析

建立場景逼真的三維地形的主要目的是對研究區的地理、地形、地貌有更清楚的了解和掌握，不僅有助于了解各類礦山地質環境要素的平面分布特征，還有助于進一步查明各地質環境點的地形地貌分布特征，從而有利于從宏觀上發現地質環境問題，便于及時決策及解決問題。結合研究區的具體地質環境問題進行的三維分析主要包括以下內容:

(1)礦山地質環境要素的空間計算。從DEM模型中提取坡度、坡向、地面剖面曲率、地面平面曲率等相關的地形因子，利用三維分析擴展模塊，結合解譯的地質環境要素矢量界線所對應的DEM高程值，量算各類地質環境要素的面積和體積。

(2)三維礦山地質環境緩沖區分析。利用遙感影像結合地面采集的各類地質災害數據，能夠迅速、準確地監測和確定各類地質災害的地理位置、面積、規模等內容。先根據地形地貌特征，在三維可視化顯示下模擬地質災害發生的情況，對所解譯地質環境要素的矢量數據進行分段緩沖區分析(緩沖區半徑根據各地質環境要素的高程計算而得);再結合專業數學模型，對地質災害的危害范圍及發展趨勢做出分析和預測。

4結論

(1)在實際應用方面，利用三維遙感可視化技術取得的影像色彩真實、立體感顯著增強，對礦山地質環境野外調查路線選擇和野外驗證具有真實觀感;形象反映了研究區的地形地貌，擺脫了以往二維平面遙感影像解譯的空間局限性，對解譯標志的建立及解譯準確率的提高具有重要作用。

(2)通過三維空間建模，采用不同的空間尺度和觀察視角，可準確、高效地提取出稀土礦開采區、尾礦堆積區、尾砂淤積區和農田淤積區的范圍、面積、空間分布等基礎信息。

(3)通過遙感影像與DEM相結合進行三維空間分析，能快速了解地質災害的影響范圍，從遙感影像中實時分析出受地質災害影響的地物，為減災排險和恢復治理提供決策依據。

［1］徐友寧.礦山地質環境調查研究現狀及展望［J］.地質通報，2008，27(8):1235 －1244.

［2］付碧宏，二宮芳樹，董彥芳，等.三維衛星遙感圖像生成技術及其在第四紀構造地貌研究中華的應用［J］.第四紀研究，2008，28(2):189－196.

［3］況順達，趙震海.遙感技術在貴州礦山地質環境調查中的應用［J］.中國礦業，2006，15(11):49 －52.

［4］孫立新，白喜慶.“3S”技術在礦山環境地質調查中的應用［J］.中國煤田地質，2004，16(S1):70 －71，83.

［5］湯國安，劉學軍，閭國年，等.數字高程模型及地學分析的原理與方法［M］.北京:科學出版社，2005.

［6］張心彬，魯孟勝，黃春慧.礦區遙感解譯與三維景觀的制作方法及應用［J］.煤炭工程，2006(3):88－90.

［7］楊武年，廖崇高，濮國梁，等.數字區調新技術新方法——遙感圖像地質解譯三維可視化及影像動態分析［J］.地質通報，2003，22(1):60 －64.

A Study of 3D Remote Sensing Technology in Geological Environment Investigation of the REE Mine

LIU Feng－mei，ZENG Min
(Wuhan Center of Geological Survey of China Geological Survey，Wuhan 430205，China)

According to the methods of RS and GIS and after geometric correction，image fusion，image mosaicking and ortho－rectification for CBERS－02B image，the remote sensing image can be overlapped with DEM generated by TIN from the relief map.On the basis of 3D remote sensing images of a REE ore exploration area in Anyuan County，three － dimensional simulation can be applied to the three－ dimensional analysis so as to study the geological environment of the mining area.Practices show that 3D remote sensing technology has great practical significance for selection of geological investigation route in areas where it is difficult to conduct field work as well as typical regional reconnaissance and mining production.In addition，it can find out the distribution pattern of the geological environment with great speed and accuracy，and can also be superimposed with vector data for the buffer analysis.Finally，the method can be used to analyze and predict the scale and development tendency of the geological hazards so as to provide the basis for decision－making on eliminating and tackling the hazards.

CBERS－02B;Geological environment in mining area;Three－dimensional analysis;Buffer analysis

TP 79

A

1001－070X(2011)04－0136－04

2011－03－22;

2011－04－22

劉鳳梅(1984－)，女，碩士，工程師，主要從事遙感、地理信息系統等地學應用方面的研究。

(責任編輯:劉心季)