許冬亮
(遼寧省自然資源事務(wù)服務(wù)中心—遼寧省基礎(chǔ)測繪院,遼寧 錦州 121003)
“國家主席習(xí)近平2014年11月10日在APEC歡迎宴會致辭時表示,希望中國藍(lán)天常在、青山常在、綠水常在。習(xí)近平總書記在黨的十九大報告中要求:“改革生態(tài)環(huán)境監(jiān)管體制”。在推進(jìn)全面完成國家自然資源綜合治理方面,要嚴(yán)厲打擊偷采盜運(yùn)等違法行為,加快采煤沉陷區(qū)、破損礦山地質(zhì)環(huán)境整治,推進(jìn)礦山生態(tài)修復(fù)。因此需要利用遙感解譯技術(shù),監(jiān)測礦山及周邊環(huán)境變化,從而提高自然資源利用,為自然資源合法、合理管理提供技術(shù)支撐。
收集整理監(jiān)測區(qū)域范圍內(nèi)多源遙感影像、基礎(chǔ)測繪成果、地理國情監(jiān)測成果等資料,經(jīng)過坐標(biāo)轉(zhuǎn)換,統(tǒng)一轉(zhuǎn)換為2000國家大地坐標(biāo)系,對遙感影像進(jìn)行解譯試驗(yàn),得到基本解譯尺度,套合多種數(shù)據(jù)成果,首先進(jìn)行自動提取,然后進(jìn)行人工判讀,最后進(jìn)行野外核查,形成監(jiān)測數(shù)據(jù)成果。對監(jiān)測成果進(jìn)行質(zhì)量檢查,建立監(jiān)測數(shù)據(jù)評價體系,制作專題地圖,為土地利用、礦山演變、環(huán)境質(zhì)量定價提供參考數(shù)據(jù)。利用遙感解譯方法進(jìn)行礦山環(huán)境監(jiān)測的技術(shù)路線(如圖1所示):
圖1 礦山環(huán)境監(jiān)測總體技術(shù)路線
通過對已有資料分析,按照礦山監(jiān)測總體技術(shù)路線要求,可以將基于遙感解譯的礦山監(jiān)測分類為土地利用和礦山監(jiān)測兩大類,其中土地利用可以監(jiān)督農(nóng)田、林地等7個分類指標(biāo),礦山監(jiān)測需要監(jiān)督礦山的破壞區(qū)域和治理區(qū)域的對比。礦山監(jiān)測分類(如圖2所示):
圖2 礦山監(jiān)測分類
2.2.1 環(huán)境評價
環(huán)境評價是在一個特定的評價范圍內(nèi),對某一既定目標(biāo)制定相關(guān)標(biāo)準(zhǔn),例如以環(huán)境作為評價對象,在有限程度上預(yù)測其發(fā)展趨勢,綜合評估人類活動對環(huán)境產(chǎn)生的影響,最終得出評價結(jié)論。我國從二十世紀(jì)七十年代末開始引入環(huán)境評價的相關(guān)機(jī)制,經(jīng)過四十余年的不斷革新,從最初的一個分類、一個標(biāo)準(zhǔn)、單一體系、普適分析,發(fā)展到如今的多個分類、多個標(biāo)準(zhǔn)、層疊體系、廣義分析,并且能夠做到動態(tài)預(yù)測、實(shí)時播報,為環(huán)境評價方法的革新作出了重要貢獻(xiàn)。隨著計(jì)算機(jī)、軟件、通信、物聯(lián)網(wǎng)等高新技術(shù)的不斷發(fā)展,信息技術(shù)在進(jìn)行環(huán)境評價時的應(yīng)用較之前更多、更廣泛,特別是3S技術(shù)的持續(xù)發(fā)展,給環(huán)境評價工作帶來了新的方向。
環(huán)境評價一般包含數(shù)量評價、標(biāo)準(zhǔn)評價和方法評價三個部分。經(jīng)過多年驗(yàn)證,這樣的分類和評價方法是科學(xué)的,環(huán)境評價的三個部分具體要求如下:
(1)數(shù)量評價:定義在有效的評價方法之上,并與評價指標(biāo)相關(guān),具備客觀性、可測性和可比性。
(2)標(biāo)準(zhǔn)評價:評價指標(biāo)應(yīng)該是客觀的,必須實(shí)現(xiàn)評價方法和環(huán)境相統(tǒng)一,并且具有一定的包含或交叉領(lǐng)域。
(3)方法評價:旨在找到評價對象的內(nèi)在聯(lián)系,并且參照相關(guān)學(xué)科的基本理論,不能脫離自然和社會,盡量做到簡單易行,減小誤差。
2.2.2 礦山環(huán)境評價
礦產(chǎn)資源的不斷開發(fā),造成非常嚴(yán)重的問題,給人類的生存環(huán)境帶來了非常大的破壞。因此,科學(xué)地評價礦山周邊環(huán)境是一個值得研究的方向。礦山周邊環(huán)境的相關(guān)評價是建立在采集、分類、走訪、計(jì)算等步驟之上,將收集到的可靠數(shù)據(jù)導(dǎo)入評價模型進(jìn)行模擬分析,在有限范圍內(nèi)預(yù)測未來一段時間內(nèi)礦山周邊環(huán)境可能出現(xiàn)的問題,做到為決策支持提供可靠數(shù)據(jù)。而且在礦山周邊環(huán)境評價開始初期,先對其存在的問題進(jìn)行研究,其次根據(jù)調(diào)查結(jié)果和基礎(chǔ)資料的精度,選擇不同的方法進(jìn)行評價。最后給出解決其存在的問題、保護(hù)和修復(fù)其治理的對策。
評價方法是在礦山環(huán)境評價時一種不可或缺的方式。通過對各種環(huán)境指標(biāo)進(jìn)行監(jiān)測或分析,找出影響環(huán)境質(zhì)量的主要指標(biāo),從而計(jì)算單指標(biāo)環(huán)境質(zhì)量評價值和綜合環(huán)境質(zhì)量評價值。評價模型對比表(如表1所示),礦山周邊環(huán)境評價方法(如表2所示):
表1 評價模型對比表
表2 礦山周邊環(huán)境評價方法
監(jiān)測區(qū)環(huán)境質(zhì)量評價的目的是通過分析、研究礦山開采活動對監(jiān)測區(qū)環(huán)境質(zhì)量的影響,從而發(fā)現(xiàn)存在的問題,然后提出環(huán)境治理措施及建議。本文選用層次分析法及粗糙集理論相結(jié)合的方法進(jìn)行礦山環(huán)境監(jiān)測研究。
2.3.1 層次分析法
層次分析法的相關(guān)計(jì)算步驟:首先是建立層次結(jié)構(gòu),其次是代入判斷矩陣。根據(jù)兩個評價指數(shù)之間的重要性關(guān)系,構(gòu)建判斷矩陣如式(1)所示:
式中表示元素i與j的相對重要性,得出的結(jié)果可以用標(biāo)準(zhǔn)度方法來表示(如表3所示):
表3 矩陣標(biāo)準(zhǔn)度及其含義
2.3.2 計(jì)算權(quán)向量
計(jì)算每行矩陣元素的乘積mi如式(2)所示:
其中i,j=(1,2,......,n)。
計(jì)算mi的n次方根如式(3)所示:
對特征向量W=[W1,W2,……,Wn]T進(jìn)行歸一化如式(4)所示:
求矩陣P的最大特征值如式(5)所示:
2.3.3 求一致性指標(biāo)C.I.如式(6)所示:
2.3.4 判斷矩陣合理性檢驗(yàn)如式(7)所示:
其中,R.I.為一致性隨機(jī)指數(shù),其值(如表4所示):
當(dāng)C.R.<0.10時,認(rèn)為判斷矩陣具有合理性,當(dāng)C.R.≥0.10時,說明判斷矩陣不滿足要求,需要進(jìn)行修改。
表4 一致性隨機(jī)指數(shù)
2.3.5 粗糙集相關(guān)理論
粗糙集相關(guān)理論提供了一種數(shù)據(jù)挖掘方法,這種方法構(gòu)建在自然語言基礎(chǔ)上,一般的分析對象基于對數(shù)據(jù)的分析和處理。它不僅是人工智能領(lǐng)域的一個熱點(diǎn),而且為信息科學(xué)研究提供了一種新的思路。粗糙集理論定義:
定義1:S=(U,A,V,f)是一個決策系統(tǒng)。其中U≠Φ為論域;A=C∪D,C和D分別為條件和決策屬性集,且C∪D≠Φ;V為屬性的值域集,Va為屬性a的值域;f:U×A→V為一個信息函數(shù),對?x∈U,a∈A,存在f(x,a)∈Va。每一個屬性子集P?A存在一個二元不可區(qū)分關(guān)系如式(8)所示:
關(guān)系IND(P),P∈A表示U的一個等價劃分,記作U/IND(P),U//ND(P)中的任意元素[X]p稱為等價類。
定義2:設(shè)X∈U,R為U上的等價關(guān)系,U/R為R對U的劃分,則如式(9)和式(10)所示:
分別稱為X關(guān)于R的下近似和上近似。
定義3:令P和Q為U中的等價關(guān)系,Q的P正域記為POSP(Q),Q和P之間的依賴程度記為其中0≤γ(Q)≤1,屬性的差異帶來的直接影響就是導(dǎo)致決策結(jié)論不一致,因此需要逐一排查特定的屬性,經(jīng)過迭代后找到變化趨勢最大的屬性,確定該屬性為對決策結(jié)論影響最大的屬性。屬性變化的幅度直接影響到?jīng)Q策結(jié)論變化的幅度。
定義4:在定義1中的決策支持關(guān)系中,屬性子集C'∈C關(guān)于D的重要定義如式(11)所示:
若C={a1,a2,......,an},則屬性ai的權(quán)重如式(12)所示:
以遼東半島南部作為研究區(qū)域,面積約157300.34公頃。該區(qū)域?yàn)闇貛Ъ撅L(fēng)氣候,年平均降雨量約為610-720毫米。研究區(qū)域影像(如圖3所示):
圖3 研究區(qū)域影像圖
本論文選擇了高分一號、高分二號及資源三號衛(wèi)星遙感影像,基礎(chǔ)測繪DLG成果,2018年地理國情監(jiān)測影像分類成果,另外還收集了地區(qū)行政區(qū)劃、研究區(qū)域內(nèi)礦權(quán)文本數(shù)據(jù)等資料。
數(shù)據(jù)處理一般包括影像融合和幾何校正。其中影像融合通過融合算法,對影像進(jìn)行處理,避免畸變,融合后地物更加清晰,提高人工判讀準(zhǔn)確性。幾何校正的目的是糾正影像的幾何變化,便于影像配準(zhǔn),幾何校正是必須的,這樣做可以消除地物變形問題。
首先將融合、幾何校正后的影像進(jìn)行多尺度分割,通過試驗(yàn)對比,得出適合試驗(yàn)區(qū)域的分割尺度,采集特征點(diǎn),自動解譯土地利用中農(nóng)田、林地、草地、水庫等類別,通過基礎(chǔ)測繪數(shù)據(jù)自動提取建設(shè)用地、道路和礦區(qū),通過人工判讀進(jìn)行修正,并人工解譯歷史影像中礦山破壞區(qū)域和礦山治理區(qū)域。礦區(qū)恢復(fù)治理前后對比影像(如圖4所示):
圖4 礦區(qū)恢復(fù)治理前后對比影像
通過監(jiān)測結(jié)果分析可知,較為嚴(yán)重影響的礦區(qū)面積約為6338.17公頃,占試驗(yàn)區(qū)域面積的4.03%,次一級嚴(yán)重影響區(qū)域面積約為21961.28公頃,占試驗(yàn)區(qū)面積的13.96%,較輕影響區(qū)域面積約為129000.89公頃,占試驗(yàn)區(qū)面積82.01%。研究試驗(yàn)區(qū)域?qū)儆诘V區(qū)治理輕度影響區(qū)域,與研究區(qū)域?qū)嶋H情況相符,說明試驗(yàn)區(qū)域影像解譯尺度合理,特征點(diǎn)采集準(zhǔn)確,可以作為研究區(qū)域礦山演變、環(huán)境保護(hù)等自然資源管理的決策依據(jù)。
通過對監(jiān)測區(qū)域范圍內(nèi)多源遙感影像進(jìn)行解譯試驗(yàn),我們不難發(fā)現(xiàn),監(jiān)測試驗(yàn)區(qū)域開發(fā)進(jìn)程加快,采礦用地在逐年增大。礦山恢復(fù)治理速度跟不上開發(fā)進(jìn)程。通過遙感解釋,直觀體現(xiàn)出礦山恢復(fù)治理率的大幅下降。監(jiān)測區(qū)域內(nèi)有關(guān)部門應(yīng)進(jìn)一步加強(qiáng)對礦山監(jiān)管和治理。
遙感解譯監(jiān)測在礦山環(huán)境監(jiān)測中,是獲得礦山環(huán)境監(jiān)測數(shù)據(jù)較為直觀且精準(zhǔn)的技術(shù)路線,可以作為大范圍推廣的礦山環(huán)境監(jiān)測手段。為生態(tài)環(huán)境監(jiān)管、治理提供了重要的數(shù)據(jù)支持及治理依據(jù),對我國環(huán)境生態(tài)保護(hù)具有重要的意義。