楊偉光,鄭有業(yè),2,劉 婷,王成松

(1.中國(guó)地質(zhì)大學(xué)(北京) 地球科學(xué)與資源學(xué)院,北京 100083;2.中國(guó)地質(zhì)大學(xué) 地質(zhì)過程與礦產(chǎn)資源國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100083;3.西藏自治區(qū)地質(zhì)礦產(chǎn)勘查開發(fā)局 第二地質(zhì)大隊(duì),西藏 拉薩 850000)

0 引 言

礦山遙感監(jiān)測(cè)是了解礦產(chǎn)資源開發(fā)現(xiàn)狀及礦山地質(zhì)環(huán)境的重要參考資料[1-2]。隨著遙感技術(shù)的快速發(fā)展，高空間分辨率遙感數(shù)據(jù)具有宏觀、快速和同步等優(yōu)點(diǎn)，在礦山遙感監(jiān)測(cè)方面的應(yīng)用越來越廣泛[3-9]。

鉻鐵礦是我國(guó)重要的戰(zhàn)略資源，近年來，隨著中國(guó)經(jīng)濟(jì)的快速發(fā)展，礦產(chǎn)資源需求強(qiáng)勁，作為我國(guó)緊缺礦種的鉻鐵礦對(duì)外依存度進(jìn)一步提高。西藏羅布莎是中國(guó)最大的鉻鐵礦基地，幾十年來的礦業(yè)開發(fā)促進(jìn)了當(dāng)?shù)亟?jīng)濟(jì)的快速發(fā)展。為了礦產(chǎn)經(jīng)濟(jì)及礦區(qū)地質(zhì)環(huán)境的可持續(xù)發(fā)展，針對(duì)當(dāng)前西藏羅布莎礦區(qū)的礦山遙感監(jiān)測(cè)尚未開展工作的情況下，利用最新的WorldView-2衛(wèi)星數(shù)據(jù)對(duì)羅布莎鉻鐵礦工作區(qū)進(jìn)行遙感調(diào)查工作，獲取地表及礦產(chǎn)開發(fā)的客觀基礎(chǔ)數(shù)據(jù)，為國(guó)土資源部門規(guī)劃資源礦產(chǎn)、維護(hù)和綜合整治礦區(qū)環(huán)境提供技術(shù)支撐及決策依據(jù)。

1 數(shù)據(jù)及研究區(qū)概況

WorldView-2衛(wèi)星由美國(guó)DigitalGlobe公司于2009年10月6日發(fā)射升空，是全球第一批使用了控制力矩陀螺(CMGs)、成像空間分辨率可達(dá)0.5 m的商業(yè)衛(wèi)星，因此在發(fā)射后很長(zhǎng)一段時(shí)間內(nèi)被認(rèn)為是全球響應(yīng)最敏捷、分辨率最高的商業(yè)成像衛(wèi)星[10]，詳細(xì)參數(shù)見表1。本文影像成像時(shí)間為2015年2月8日，含有1個(gè)全色波段和8個(gè)多光譜段，影像無云層覆蓋，地物細(xì)節(jié)清晰。

羅布莎地處西藏自治區(qū)山南地區(qū)，位于首府拉薩市東南方向約200 km，行政區(qū)劃屬西藏曲松縣境內(nèi)，交通十分便利[11]。礦產(chǎn)資源以鉻鐵礦石為主，是我國(guó)規(guī)模最大、礦石質(zhì)量最高的鉻鐵礦床。

表1 WorldView-2衛(wèi)星及波段參數(shù)Table 1 Satellite parameters of WorldView-2

2 影像預(yù)處理

2.1 正射校正

研究區(qū)地形起伏較大，高分辨率WorldView-2數(shù)據(jù)應(yīng)通過正射校正來消除地形起伏引起影像數(shù)據(jù)的畸變[12]。本文正射校正采用RPC有理多項(xiàng)式模型，利用DEM數(shù)據(jù)和地面控制點(diǎn)參與來提高校正精度，達(dá)到《礦山遙感調(diào)查與監(jiān)測(cè)項(xiàng)目技術(shù)要求》的相關(guān)要求。

2.2 大氣校正

大氣校正的目的是消除大氣等對(duì)影像成像的影響[13]，本文采用FLAASH大氣校正，它是基于MODTRAN5輻射傳輸模型，是ENVI平臺(tái)下精度較高的大氣校正模型。

2.3 影像融合

本文影像融合的目的是提高空間分辨率，增強(qiáng)目標(biāo)特征，提高分類精度[13]。融合方法采用PCI下超分辨率貝葉斯融合算法，即PanSharpening，這種方法使全色波段數(shù)據(jù)與多光譜波段數(shù)據(jù)自動(dòng)對(duì)齊，保留了光譜信息，豐富了地面信息，提高了地物分辨能力。

2.4 波段選擇

根據(jù)本文研究的目的，要最大化利用WorldView-2高空間分辨率的優(yōu)勢(shì)，采用最佳指數(shù)法(OIF)進(jìn)行波段選擇[14]，經(jīng)過數(shù)學(xué)統(tǒng)計(jì)分析，選擇近紅外2(Band 8)，紅色邊緣(Band 6)，黃色(Band 4)3個(gè)波段，分別置于R、G、B通道，作為信息提取基礎(chǔ)影像。同時(shí)結(jié)合紅、綠、藍(lán)真彩色合成影像，接近自然色彩，有利于礦山信息的識(shí)別與解譯[9]。

3 礦區(qū)遙感解譯

礦山遙感監(jiān)測(cè)涉及的解譯對(duì)象種類繁多，主要涉及的目標(biāo)地物有采場(chǎng)、礦山道路、開采斜井、尾礦、礦山建筑、附近的村莊、水電站、耕地等基本礦區(qū)環(huán)境。本文以WorldView-2最佳波段組合影像(Band 864)及野外驗(yàn)證來說明研究區(qū)內(nèi)目標(biāo)地物影像特征，準(zhǔn)確建立高分辨率遙感解譯標(biāo)志，對(duì)于監(jiān)測(cè)礦山開采和環(huán)境治理具有重要意義[9]。

3.1 采場(chǎng)及采區(qū)道路

羅布莎鉻鐵礦區(qū)多采取露天開采方式，規(guī)模巨大的采坑，色調(diào)與周圍地表環(huán)境明顯不同，人為活動(dòng)與地貌破壞明顯，表面光滑，無植被發(fā)育，采坑范圍易于判讀，部分礦坑在遙感影像表現(xiàn)為紋理較為粗糙，有不規(guī)則斑塊分布。運(yùn)輸?shù)缆非逦梢?，呈寬窄不一的線狀圖案(圖1)。

3.2 開采斜井

開采斜井指傾斜直達(dá)地下礦體的巷道。開采的硐口與地表呈近垂直，即使是在高空間分辨率影像上也難以直接解譯，因此通常利用間接解譯標(biāo)志進(jìn)行判讀[15]。例如為方便生產(chǎn)的礦山用房，或礦石滑道、礦車鐵軌及就近順山坡堆放的開采廢渣等。圖2展示的為江南礦業(yè)4266號(hào)斜井及西藏礦業(yè)4030號(hào)斜井的影像特征及野外查證。

圖1 羅布莎礦區(qū)露天采場(chǎng)及道路影像Fig.1 Open stopes and roads image(Band 864) of Luobusa mining area

3.3 排石(土)場(chǎng)

排石(土)場(chǎng)又稱廢石場(chǎng)，是指礦山采礦過程中產(chǎn)生的大量固體廢棄物[15]。一般靠近開采場(chǎng)附近堆積，呈近圓形、扇形圖斑，是一種巨型人工松散堆墊體[15]，較周圍自然地物，其粒徑小，影像上紋理細(xì)膩，表面平整，同時(shí)較周圍地物明亮，在羅布莎礦區(qū)內(nèi)十分常見，如圖3所示。

3.4 尾礦庫(kù)

羅布莎鉻鐵礦尾礦地處山谷，地物特征較為明顯，可見較直或呈弧度彎曲的壩堤攔截谷口，影像上具明顯的線條，距離礦區(qū)不遠(yuǎn)(圖4)，但要注意的是，羅布莎地區(qū)屬于高原溫帶季風(fēng)濕潤(rùn)氣候，尤其正值夏秋兩季礦山活動(dòng)時(shí)，降水較多，要防止水質(zhì)污染。

3.5 礦區(qū)建筑

在羅布莎礦區(qū)開展工作的單位主要包括西藏江南礦業(yè)股份有限公司、西藏礦業(yè)發(fā)展股份有限公司、西藏自治區(qū)地質(zhì)礦產(chǎn)勘查開發(fā)局第二地質(zhì)大隊(duì)等，因此礦區(qū)建筑相對(duì)分散復(fù)雜，但主要包括了礦部、礦區(qū)宿舍、工棚、炸藥庫(kù)等，一般具有規(guī)則的矩形特征，部分建筑立體效果顯著，頂棚材料反射率高，同時(shí)可見部分機(jī)械裝備(圖5)。

圖2 羅布莎礦區(qū)斜井或平硐遙感影像及野外查證Fig.2 Inclined shaft images(Band 864) and field photos of Luobusa mining area

圖3 羅布莎礦區(qū)開采坑、排石場(chǎng)遙感影像及野外查證Fig.3 Open stope,waste dump image(Band 864)and field photo of Luobusa mining area

3.6 地表地物

在高空間分辨率遙感影像上，羅布莎地區(qū)地表各地物輪廓清晰，易判讀[16]。主要有之字形道路穿越羅布莎村，民用房屋清晰可見，近南北向水系平行羅布莎村；沃卡水電站距鉻鐵礦勘查區(qū)北3 km左右，以水力發(fā)電為主，兼具小規(guī)模灌溉功能，最后匯流至雅魯藏布江；羅布莎地區(qū)農(nóng)耕土地多呈階梯狀分布于雅魯藏布江沿岸，農(nóng)作物以青稞、土豆、大豆為主，附近可見呈線性展布的高壓電塔(圖6)。

圖4 羅布莎礦區(qū)尾礦庫(kù)影像Fig.4 Mine tailing pond images(Band 864) of Luobusa mining area

圖5 羅布莎礦區(qū)建筑影像及野外照片F(xiàn)ig.5 Mine construction images(Band 864)and field photo of Luobusa mining area

圖6 羅布莎地表地物影像Fig.6 Surface feature images(Band 864) of Luobusa

3.7 生態(tài)環(huán)境

圖7 羅布莎地區(qū)地表生態(tài)環(huán)境影像Fig.7 Surface ecological environment images(Band 864) of Luobusa area

羅布莎礦山開采帶來的環(huán)境影響主要體現(xiàn)在破壞地貌景觀和地表植被，同時(shí)可能引發(fā)滑坡等非自然地質(zhì)災(zāi)害，如圖7所示。

4 數(shù)據(jù)分析與建議

4.1 礦山開發(fā)現(xiàn)狀統(tǒng)計(jì)

本文利用WorldView-2遙感數(shù)據(jù)對(duì)羅布莎鉻鐵礦區(qū)進(jìn)行了礦山開發(fā)遙感調(diào)查，調(diào)查面積4 500 hm2。

調(diào)查結(jié)果顯示，研究區(qū)范圍內(nèi)露天采坑較多，呈群呈帶狀分布，開采的礦種主要是鉻鐵礦石，采場(chǎng)占地總面積約152 hm2，礦山建筑面積約3.64 hm2，固體堆積物或尾礦庫(kù)占地面積約13 hm2。

4.2 礦區(qū)生態(tài)環(huán)境分析

上述占用的土地類型主要是草地和農(nóng)耕地。首先這些采場(chǎng)、排土場(chǎng)等礦業(yè)活動(dòng)占用耕種土地，破壞了地面植被，造成土地利用率低；其次西藏羅布莎地區(qū)降雨頻繁，在表生帶強(qiáng)氧化條件下，廢石中的重金屬元素等的遷移，可能對(duì)水質(zhì)與土壤產(chǎn)生影響；再者采礦活動(dòng)和廢石堆積可能加劇滑坡等地質(zhì)災(zāi)害發(fā)生的頻率和強(qiáng)度；同時(shí)目前羅布莎礦區(qū)鉻鐵礦的開采方式基本依靠地下開采，這可能加速礦區(qū)周圍地形地貌的變化。

4.3 建議

在野外實(shí)地勘察期間，當(dāng)?shù)鼐用癖磉_(dá)了對(duì)礦山環(huán)境問題的擔(dān)憂。因此本文首次開展羅布莎礦山遙感監(jiān)測(cè)是十分必要和有意義的，針對(duì)如何合理規(guī)劃礦山開采活動(dòng)及如何恢復(fù)當(dāng)?shù)赜刹傻V活動(dòng)引起的生態(tài)環(huán)境問題已迫在眉睫。

在羅布莎地表礦開采基本結(jié)束的情況下，應(yīng)積極對(duì)地表廢棄礦坑進(jìn)行復(fù)墾復(fù)綠，增加土地利用率，對(duì)排土場(chǎng)進(jìn)行人工改造，加強(qiáng)對(duì)廢石的綜合利用，同時(shí)要跟進(jìn)礦山開發(fā)環(huán)境影響評(píng)價(jià)工作，制定針對(duì)羅布莎礦山環(huán)境保護(hù)治理規(guī)劃與措施，并在下一步礦業(yè)活動(dòng)開始之前，對(duì)可能產(chǎn)生的環(huán)境影響進(jìn)行分析和評(píng)估。

5 結(jié) 論

(1)高空間分辨率遙感在礦區(qū)遙感監(jiān)測(cè)中應(yīng)用效果顯著，成為礦山開采狀況及環(huán)境調(diào)查中不可缺少的基礎(chǔ)數(shù)據(jù)。將高空間分辨率遙感地質(zhì)調(diào)查與礦山遙感監(jiān)測(cè)工作結(jié)合起來，可實(shí)現(xiàn)一種數(shù)據(jù)，多目標(biāo)調(diào)查，增加數(shù)據(jù)利用價(jià)值。

(2)開展前期野外踏勘及后期野外查證是相當(dāng)必要的工作程序，有利于提高遙感解譯精度。本文在結(jié)合野外工作的前提下，建立了研究區(qū)內(nèi)礦山開采過程中涉及到的地物及地表環(huán)境的高空間分辨率遙感解譯標(biāo)志。

(3)遙感調(diào)查結(jié)果顯示，羅布莎礦山環(huán)境的恢復(fù)治理工作任重道遠(yuǎn)，在關(guān)注礦產(chǎn)經(jīng)濟(jì)發(fā)展的同時(shí)，應(yīng)加大對(duì)礦山環(huán)境治理工作的投入，因此本次的研究成果可以服務(wù)于下一次羅布莎地區(qū)礦產(chǎn)資源環(huán)境規(guī)劃決策。

參考文獻(xiàn)：

[1] 楊金中，秦緒文，張志，等.礦山遙感監(jiān)測(cè)理論方法與實(shí)踐[M].北京: 測(cè)繪出版社,2011:1-200.

[2] 秦緒文，楊金中，康高峰，等.礦山遙感監(jiān)測(cè)技術(shù)方法研究[M].北京: 測(cè)繪出版社,2011:1-308.

[3] 王燕波，羅偉，李名勇，等.基于高分辨率遙感影像的礦山開發(fā)監(jiān)測(cè)研究[J].熱帶地理,2011,31(4): 377-382.

[4] 張焜，馬世斌，劉麗萍.基于SPOT5數(shù)據(jù)的鹽湖礦產(chǎn)開發(fā)及礦山環(huán)境遙感監(jiān)測(cè)[J].國(guó)土資源遙感,2012,24(3): 146-153.

[5] 路云閣，劉采，王姣.基于國(guó)產(chǎn)衛(wèi)星數(shù)據(jù)的礦山遙感監(jiān)測(cè)一體化解決方案——以西藏自治區(qū)為例[J].國(guó)土資源遙感,2014,26(4): 85-90.

[6] 楊顯華，黃潔，田立，等.基于高分辨率遙感數(shù)據(jù)的礦山環(huán)境綜合治理研究——以冕寧牦牛坪稀土礦為例[J].國(guó)土資源遙感,2015,27(4): 115-121.

[7] 代晶晶，王登紅，吳亞楠.基于高分遙感數(shù)據(jù)的稀有礦山監(jiān)測(cè)——以江西宜春414稀有礦山為例[J].國(guó)土資源遙感,2017,29(3): 104-110.

[8] 王耿明，朱俊鳳，武國(guó)忠，等.廣東大寶山多金屬礦開發(fā)環(huán)境遙感監(jiān)測(cè)與礦山地質(zhì)環(huán)境評(píng)價(jià)[J].物探化探計(jì)算技術(shù),2017,39(1): 122-128.

[9] 于博文，田淑芳，趙永超，等.高分一號(hào)衛(wèi)星在京津礦山遙感監(jiān)測(cè)中的應(yīng)用[J].現(xiàn)代地質(zhì),2017,31(4): 843-850.

[10] MULLER S J,NIEKERKA V.Identification of WorldView-2 spectral and spatial factors in detecting salt accumulation in cultivated fields[J].Geoderma,2016,273: 1-11.

[11] 梁鳳華，許志琴，巴登珠，等.西藏羅布莎—澤當(dāng)蛇綠巖體的構(gòu)造產(chǎn)出與侵位機(jī)制探討[J].巖石學(xué)報(bào),2011,27(11): 3255-3268.

[12] 楊金中，孫延貴，秦緒文，等.高分辨率遙感地質(zhì)調(diào)查[M].北京: 測(cè)繪出版社,2013：1-225.

[13] 趙英時(shí).遙感應(yīng)用分析原理與方法[M].北京: 科學(xué)出版社,2002：1-478.

[14] 焦?jié)櫝?，秦彥平，張淑云，?WorldView-2數(shù)據(jù)在沉積巖地區(qū)的遙感巖性增強(qiáng)方法初探[J].西北地質(zhì),2014,47(4): 278-283.

[15] 王蔚，傅濤.礦業(yè)開發(fā)活動(dòng)的高分辨率遙感影像解譯標(biāo)志[J].云南地質(zhì),2012,31(2): 242-244.

[16] 安志宏，聶洪峰，王昊，等.ZY-102C星數(shù)據(jù)在礦山遙感監(jiān)測(cè)中的應(yīng)用研究與分析[J].國(guó)土資源遙感,2015,27(2): 175-182.