尹建平

（中煤平朔集團(tuán)有限公司，山西 朔州 036006）

隨著社會(huì)經(jīng)濟(jì)及基本物質(zhì)條件的快速發(fā)展，以及自然資源、空間環(huán)境的不斷改變，導(dǎo)致土壤退化、水土流失、生物多樣性減少等一系列環(huán)境問(wèn)題，給生態(tài)文明建設(shè)和社會(huì)可持續(xù)發(fā)展帶來(lái)嚴(yán)重挑戰(zhàn)。如何準(zhǔn)確快速地評(píng)估區(qū)域內(nèi)生態(tài)環(huán)境狀況，成為當(dāng)今迫切需要解決的問(wèn)題。傳統(tǒng)的生態(tài)環(huán)境質(zhì)量評(píng)價(jià)方法包括模糊評(píng)價(jià)法[1]、層次分析法[2]、PSR 模型[3]等，這些方法受主觀因素影響較大，評(píng)價(jià)結(jié)果往往偏離區(qū)域內(nèi)實(shí)際狀況；主成分分析法[4]、人工神經(jīng)網(wǎng)絡(luò)法[5]的出現(xiàn)解決了主觀權(quán)重賦值的問(wèn)題，使評(píng)估更具科學(xué)性。2006 年環(huán)境保護(hù)部發(fā)布的《生態(tài)環(huán)境狀況評(píng)價(jià)技術(shù)規(guī)范（試行）》為生態(tài)環(huán)境狀況評(píng)價(jià)提供了標(biāo)準(zhǔn)，生態(tài)環(huán)境狀況指數(shù)EI 得到廣泛應(yīng)用，但其指標(biāo)權(quán)重需要人為設(shè)定，且存在評(píng)價(jià)結(jié)果無(wú)法空間可視化的問(wèn)題。近年來(lái)，眾多專家學(xué)者利用遙感與地理信息技術(shù)對(duì)城市、沙漠、黃土流失區(qū)和自然保護(hù)區(qū)等不同區(qū)域進(jìn)行生態(tài)環(huán)境評(píng)價(jià)[6-7]。

平朔目前已經(jīng)成為我國(guó)規(guī)模最大，現(xiàn)代化程度最高的煤炭基地之一。隨著煤炭的開采，對(duì)環(huán)境的破壞也是日益嚴(yán)重。1985 年平朔礦區(qū)開始進(jìn)行初級(jí)的復(fù)墾工作，隨著國(guó)家對(duì)生態(tài)文明建設(shè)工作的重視和人們環(huán)保意識(shí)的不斷增強(qiáng)，平朔礦區(qū)開始實(shí)行更高級(jí)的生態(tài)重建工作，實(shí)行“采，運(yùn)，排，復(fù)”一體化工藝，“地貌重塑、土壤重構(gòu)，植被重建，水土保持，植被配置模式，堆狀排棄工藝”等生態(tài)環(huán)境重構(gòu)的綜合技術(shù)，對(duì)我國(guó)黃土高原區(qū)生態(tài)恢復(fù)治理工作起到了積極的示范作用。

目前，對(duì)于平朔露天礦的研究主要集中在植被恢復(fù)與土壤因子的關(guān)系、地表覆被變化及土地利用結(jié)構(gòu)變化和景觀格局等領(lǐng)域，鮮見針對(duì)礦區(qū)進(jìn)行生態(tài)環(huán)境狀況評(píng)估的研究。為此，基于遙感生態(tài)指數(shù)RSEI，集成綠度、濕度、熱度和干度4 項(xiàng)指標(biāo)，利用主成分分析法建立評(píng)估模型，綜合分析1989—2018年礦區(qū)涉及生態(tài)環(huán)境狀況的指標(biāo)變化，著重分析近30 年的生態(tài)環(huán)境變化程度，評(píng)估平朔露天礦的生態(tài)環(huán)境狀況，揭示平朔露天礦開始復(fù)墾工作后其生態(tài)環(huán)境恢復(fù)情況，以期為礦區(qū)生態(tài)環(huán)境治理提供科學(xué)依據(jù)和理論支撐。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況

礦區(qū)位于山西省朔州市平魯區(qū)寧武煤田北端（N39°23′，E112°30′），礦區(qū)南北長(zhǎng)23 km，東西寬22 km，總面積為380 km2，已探明地質(zhì)儲(chǔ)量127.5 億t?，F(xiàn)包括安太堡露天礦、安家?guī)X露天礦、東露天礦、井工一礦、井工二礦（已停產(chǎn)）、井工三礦6 個(gè)煤礦及其配套選煤廠。該區(qū)屬黃土丘陵地貌，溫帶半干旱大陸性氣候，海拔在1 250～1 350 m 水平，年均氣溫4.8～7.5℃，年均降水量428.2～449 mm，礦區(qū)土壤主要為栗褐土，其土壤結(jié)構(gòu)差，水土嚴(yán)重流失。

1.2 數(shù)據(jù)源及研究方法

在中國(guó)地理空間數(shù)據(jù)云網(wǎng)站（http：//www.gscloud.cn/）獲取了1989 年6 月20 日、1998 年7 月11 日、2007 年6 月18 日的Landsat5 TM 數(shù)據(jù)，2018年7 月1 日的Landsat 8 OLI 數(shù)據(jù)影像資料及矢量邊界數(shù)據(jù)和Google Earth 影像數(shù)據(jù)。云層覆蓋面積小于10%，空間分辨率15 m/30 m，研究主要使用ENVI Classic 軟件來(lái)處理圖像并且結(jié)合使用ENVI5.3和ArcGIS 10.2 空間處理軟件。

選取易獲取的歸一化植被指數(shù)、濕度分量、干度指數(shù)和地表溫度代表綠度NDVI、濕度WET、干度NDSI 和熱度LST，對(duì)各項(xiàng)指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理，利用ArcGIS 軟件的主成分分析工具將遙感生態(tài)指數(shù)RSEI 表示為4 項(xiàng)指標(biāo)合成的函數(shù)。

1.2.1 評(píng)價(jià)因子標(biāo)準(zhǔn)化處理

RSEI 各分量指標(biāo)綠度、濕度、干度以及熱度，它們的量綱不同，沒法進(jìn)行相互比較，故需要進(jìn)行標(biāo)準(zhǔn)化處理。

式中：X、Y 分別為指數(shù)原值和標(biāo)準(zhǔn)化值；Xmax為指數(shù)的最大值；Xmin為指數(shù)的最小值。

1.2.2 遙感生態(tài)指數(shù)計(jì)算

首先將經(jīng)過(guò)標(biāo)準(zhǔn)化處理的4 個(gè)評(píng)價(jià)因子利用ENVI 5.3 軟件Layer Stacking 工具進(jìn)行波段合成，然后再利用軟件中主成分分析模塊進(jìn)行主成分分析，將主成分1（PCA1）作為初始遙感生態(tài)指數(shù)RSEI0：

RSEI0=PCA[（fNDVI，WET，NDBSI，LSTI）]

為了方便各時(shí)期遙感生態(tài)指數(shù)之間進(jìn)行比較，同樣對(duì)RSEI0進(jìn)行標(biāo)準(zhǔn)化處理，按照如下公式進(jìn)行：

式中：RSEI為對(duì)應(yīng)各期經(jīng)過(guò)標(biāo)準(zhǔn)化的生態(tài)遙感指數(shù)，值域[0，1]；RSEI0－min為原始遙感生態(tài)指數(shù)中最小值；RSEI0-max為原始遙感生態(tài)指數(shù)中最大值。

RSEI 值越大表示生態(tài)環(huán)境狀況越好，反之則代表生態(tài)環(huán)境狀況越差。

2 結(jié)果與分析

2.1 平朔露天礦區(qū)生態(tài)環(huán)境狀況分析

1989—2018 年4 個(gè)指標(biāo)的主成分分析結(jié)果見表1（已篩選第1 主成分）。

表1 1989—2018 年4 個(gè)指標(biāo)的主成分分析結(jié)果（CP1 值）

由計(jì)算結(jié)果可知1989—2018 年綠度、濕度、干度和熱度的值，其中綠度和濕度與生態(tài)環(huán)境狀況呈正相關(guān)關(guān)系，共同對(duì)礦區(qū)的生態(tài)環(huán)境狀況起正面的促進(jìn)作用。干度和熱度與生態(tài)環(huán)境狀況呈負(fù)相關(guān)關(guān)系，兩者共同起負(fù)面作用。由表1 可知各年度4 項(xiàng)指標(biāo)第1 主成分的貢獻(xiàn)率均大于65%，貢獻(xiàn)較好，表明第1 主成分已經(jīng)集成4 個(gè)指標(biāo)的大部分特征。由各指標(biāo)的第1 主成分載荷可知，與該區(qū)域生態(tài)環(huán)境狀況呈正相關(guān)關(guān)系的濕度、綠度中，兩者的貢獻(xiàn)率較低甚至有的年份為負(fù)值，表明植被生長(zhǎng)及分布情況是礦區(qū)生態(tài)環(huán)境狀況的直接影響因素。

2.2 平朔礦區(qū)生態(tài)環(huán)境狀況變化

1989—2018 年平朔礦區(qū)遙感生態(tài)指數(shù)RSEI 見表2。

表2 1989—2018 年平朔礦區(qū)遙感生態(tài)指數(shù)RSEI

由表2 及空間分布（圖略）來(lái)看，1989 年生態(tài)環(huán)境狀況為優(yōu)的區(qū)域零性散落在整個(gè)礦區(qū)內(nèi)，1994 年生態(tài)狀況為優(yōu)的區(qū)域主要位于礦區(qū)南部，整體看來(lái)生態(tài)環(huán)境質(zhì)量比1989 年有所好轉(zhuǎn)。東部地區(qū)生態(tài)環(huán)境質(zhì)量下降明顯，這是因?yàn)椴傻V初期東部地區(qū)是采礦的重心位置，建設(shè)用地面積增多，植被稀疏，覆蓋度低，生態(tài)破壞較為嚴(yán)重，而北部地區(qū)由于采礦活動(dòng)較少，對(duì)植被、耕地等保護(hù)措施較好，植被茂盛，也與該時(shí)期林地面積增加有關(guān)。但是在1994—2007年，由于采礦活動(dòng)重心位置的北移，礦區(qū)生態(tài)環(huán)境仍處于下降階段，在該段時(shí)期內(nèi)，礦區(qū)草地、耕地持續(xù)減少，采礦用地、裸地持續(xù)增加是導(dǎo)致生態(tài)環(huán)境質(zhì)量下降的主要原因，在圖上我們可以明顯的看到礦區(qū)北部及南部區(qū)域生態(tài)質(zhì)量明顯惡化。

2007 年之后礦區(qū)逐步開展復(fù)墾工作，復(fù)墾范圍主要集中在礦區(qū)南部安家?guī)X、安太堡露天礦周圍，由圖2 實(shí)例可以看出10 年間生態(tài)環(huán)境質(zhì)量明顯變好。2018 年RSEI 上升0.018，隨著復(fù)墾工作的進(jìn)行，礦區(qū)生態(tài)質(zhì)量持續(xù)變好，林地、耕地面積增加，建設(shè)用地、裸地面積減少。復(fù)墾工作進(jìn)行至2018 年RSEI 值已基本與1989 年建礦初期時(shí)的RSEI 值接近，說(shuō)明礦區(qū)在擴(kuò)大采礦的同時(shí)兼顧到了礦區(qū)生態(tài)環(huán)境的維護(hù)治理，近10 年來(lái)的復(fù)墾工作頗有成效。

3 結(jié)論

1）平朔露天礦1989、1994、1998、2003、2007、2013、2018 年的遙感生態(tài)指數(shù)值分別為0.564、0.592、0.547、0.527、0.512、0.558、0.576，數(shù)據(jù)表明RSEI 值呈先升高后降低再升高的變化過(guò)程，結(jié)合空間分布圖整體來(lái)看，礦區(qū)生態(tài)環(huán)境狀況為良性發(fā)展趨勢(shì)。

2）綠度和濕度對(duì)礦區(qū)生態(tài)環(huán)境狀況起顯著的正向作用，主要是因?yàn)檠芯繀^(qū)植被覆蓋度隨復(fù)墾年限的增加而增加，而干度和熱度指標(biāo)對(duì)礦區(qū)生態(tài)環(huán)境狀況起負(fù)向作用。

3）通過(guò)分析遙感數(shù)據(jù)及空間分布情況可實(shí)現(xiàn)快速地對(duì)平朔露天礦大面積區(qū)域進(jìn)行生態(tài)環(huán)境狀況評(píng)估，同時(shí)運(yùn)用主成分分析法可有效避免評(píng)價(jià)模型主觀因素的影響，提高了評(píng)估結(jié)果的客觀性，對(duì)于其他礦區(qū)進(jìn)行相關(guān)方面的研究有一定的指導(dǎo)意義。