李 蕊, 陳國清, 李瑋賢, 蒙 榮, 王明玖, 郭洋楠

(1.內(nèi)蒙古農(nóng)業(yè)大學(xué) 沙漠治理學(xué)院, 內(nèi)蒙古 呼和浩特 010018; 2.內(nèi)蒙古農(nóng)業(yè)大學(xué) 草原與資源環(huán)境學(xué)院, 內(nèi)蒙古 呼和浩特 010018; 3.神華神東煤炭集團(tuán)有限公司技術(shù)研究所, 陜西 榆林 719315)

隨著經(jīng)濟(jì)全球化的發(fā)展，生態(tài)環(huán)境問題儼然成為全球性關(guān)注熱點(diǎn)。礦業(yè)是中國經(jīng)濟(jì)支柱產(chǎn)業(yè)，其發(fā)展在一定程度上以破壞環(huán)境為代價(jià)，礦業(yè)在取得巨大成就的同時(shí)，也對(duì)生態(tài)環(huán)境造成了嚴(yán)重破壞。神東礦區(qū)是中國重要的區(qū)域生態(tài)脆弱帶以及國家級(jí)水土流失重點(diǎn)監(jiān)督區(qū)和重點(diǎn)治理區(qū)，生態(tài)環(huán)境脆弱，干旱少雨；原生植被種類單調(diào)，平均植被覆蓋率僅3%～11%，風(fēng)蝕區(qū)面積占總面積的70%，長(zhǎng)此以往，會(huì)使本來脆弱的生態(tài)環(huán)境變得徹底無法修復(fù)，面臨生態(tài)環(huán)境惡化的嚴(yán)重制約，不僅經(jīng)濟(jì)難以持續(xù)發(fā)展，還會(huì)影響晉陜蒙接壤地區(qū)的生態(tài)安全。所以深入研究礦區(qū)生態(tài)環(huán)境質(zhì)量時(shí)空變化規(guī)律，對(duì)于恢復(fù)礦區(qū)生態(tài)功能，提高礦區(qū)生態(tài)環(huán)境質(zhì)量等具有十分重要的理論價(jià)值和現(xiàn)實(shí)意義[1-2]。

2015年環(huán)保部發(fā)布了《生態(tài)環(huán)境狀況評(píng)價(jià)技術(shù)規(guī)范》(HJ192-2015)，提出生態(tài)環(huán)境狀況評(píng)價(jià)指數(shù)(ecological environment index，EI)[3-4]，而構(gòu)成EI指數(shù)的各個(gè)指標(biāo)在數(shù)據(jù)可獲取性、權(quán)重合理性和評(píng)價(jià)結(jié)果時(shí)空分析等方面仍存有弊端[5]。徐涵秋基于遙感技術(shù)提出了遙感生態(tài)指數(shù)(remote sensing ecology index，RSEI)，由綠度、濕度、熱度和干度4個(gè)遙感生態(tài)指標(biāo)構(gòu)建[6-7]，該指數(shù)不僅數(shù)據(jù)易獲得、權(quán)重非人為設(shè)定，而且具有較強(qiáng)的客觀性、評(píng)價(jià)結(jié)果可視化效果好等優(yōu)點(diǎn)，因此與EI指數(shù)相比該指數(shù)更適用于區(qū)域生態(tài)環(huán)境質(zhì)量評(píng)價(jià)。由于礦區(qū)生態(tài)環(huán)境受損的直接體現(xiàn)就是礦區(qū)土壤理化性質(zhì)、土地利用方式、植被的時(shí)空變化等，所以對(duì)礦區(qū)生態(tài)環(huán)境質(zhì)量的研究主要集中在礦區(qū)生態(tài)效應(yīng)[8]、生態(tài)安全[9]、承載力等[10]方面。但是這些研究和規(guī)范的目標(biāo)是揭示礦區(qū)生態(tài)環(huán)境質(zhì)量的現(xiàn)狀及變化，大部分文獻(xiàn)評(píng)價(jià)了生態(tài)環(huán)境質(zhì)量包括水環(huán)境[11]、植被環(huán)境和土壤環(huán)境[12-13]，部分涉及到土地利用景觀格局[14]，在神東礦區(qū)生態(tài)環(huán)境質(zhì)量時(shí)空演變規(guī)律的研究比較缺乏。本文在RSEI指數(shù)的基礎(chǔ)上，利用多時(shí)相遙感數(shù)據(jù)綜合研究神東礦區(qū)生態(tài)環(huán)境質(zhì)量變化及發(fā)展趨勢(shì)，同時(shí)利用重心遷移法來探究生態(tài)環(huán)境質(zhì)量空間變化狀況，以期為該研究區(qū)的空間變化情況及生態(tài)環(huán)境建設(shè)提供一定的科學(xué)依據(jù)，也為神東礦區(qū)完善生態(tài)治理制度體系提供參考依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

本研究選取神東礦區(qū)(為神府—東勝礦區(qū)的簡(jiǎn)稱)作為研究區(qū)，位于內(nèi)蒙古自治區(qū)鄂爾多斯東南部和陜西北部交界處，地理坐標(biāo)在北緯38°52′—39°41′，東經(jīng)109°51′—110°46′。實(shí)際研究區(qū)坐標(biāo)為38°56′—39°49′N，108°58′—110°25′E。地處西北地區(qū)毛烏素沙漠東南邊緣與黃土高原北端的過渡地帶，西北高，東南低，屬典型的干旱、半干旱荒漠高原氣候。礦區(qū)北部東勝梁呈東西向展布，形成南北向河流的分水嶺。東部及東北部為黃土丘陵山區(qū)，溝壑縱橫，形成梁峁、溝壑和土原3種地貌。

1.2 數(shù)據(jù)來源及預(yù)處理

本文所用的數(shù)據(jù)為神東礦區(qū)(行列號(hào)127，33)1995—2020年的Landsat 5 TM和Landsat 8 OLI遙感影像，數(shù)據(jù)通過“美國地質(zhì)調(diào)查局USGS”(http:∥glovis.usgs.gov/)進(jìn)行免費(fèi)數(shù)據(jù)的獲取。其中包括1995，2000，2005，2010年的Landsat 5 TM與2015，2020年的Landsat 8 OLI共6期的遙感影像數(shù)據(jù)，采用的遙感影像的質(zhì)量較好。由于綠度指標(biāo)(NDVI)選用的數(shù)據(jù)為一期影像而非多期求均值，為避免綠度指標(biāo)過高影響評(píng)價(jià)指數(shù)準(zhǔn)確性，且受到數(shù)據(jù)獲取限制，所以選用的時(shí)間基本集中在9—11月份影像晴空無云。

運(yùn)用遙感分析軟件ENVI 5.3對(duì)6個(gè)不同時(shí)期的影像進(jìn)行輻射定標(biāo)、大氣校正等，為滿足研究精度的需要，將均方根誤差控制在0.5個(gè)像元內(nèi)，最后使用神東礦區(qū)的矢量邊界對(duì)圖形進(jìn)行裁剪。

1.3 研究方法

1.3.1 生態(tài)環(huán)境質(zhì)量評(píng)價(jià) RSEI指數(shù)是用于綜合研究區(qū)域生態(tài)環(huán)境現(xiàn)狀的新型遙感生態(tài)指數(shù)，由濕度指標(biāo)、綠度指標(biāo)、熱度指標(biāo)、干度指標(biāo)組成的，利用主成分分析構(gòu)建神東礦區(qū)遙感生態(tài)指數(shù)，即：

RSEI0=1-{PC1〔f(NDVI，WET，LST，NDSI)〕}

(1)

(2)

式中：RSEI為遙感生態(tài)指數(shù); NDVI為綠度指標(biāo); WET為濕度指標(biāo); LST為熱度指標(biāo); NDSI為干度指標(biāo); PC1為主成分分析后得到的第一主成分; RSEI0為遙感生態(tài)指數(shù)初始值;其中RSEI0-min，RSEI0-max分別代表RSEI0的最小值、最大值。

(1) 濕度。濕度指標(biāo)選用遙感纓帽變換的濕度分量，濕度分量對(duì)神東礦區(qū)的土壤和植被的濕度有較好的反映[15-16]。Landsat 5 TM和Landsat 8 OLI的濕度計(jì)算公式分別為：

WETTM=0.031 5ρblue+0.202 1ρgreen+

0.310 2ρred+0.159 4ρNIR+

0.680 6ρSWIR1+0.610 9ρSWIR2

(3)

WETOLI=0.151 1ρblue+0.197 2ρgreen+

0.328 3ρred+0.340 7ρNIR+

0.711 7ρSWIR1+0.455 9ρSWIR2

(4)

式中：ρblue為藍(lán)波段反射率;ρgreen為綠波段反射率;ρSWIR1為短波紅外1波段反射率;ρSWIR2為短波紅外2波段反射率。

(2) 綠度。綠度指標(biāo)選用歸一化植被指數(shù)(NDVI)代替綠度，它主要反映植物生物量、葉面積指數(shù)以及植被覆蓋度，是使用最廣泛的植被指數(shù)[17]。具體表達(dá)式為：

NDVI=(ρNIR-ρRed)/(ρNIR+ρRed)

(5)

式中：ρNIR，ρRed分別為L(zhǎng)andsat 5 TM,Landsat 8 OLI的近紅外波段和紅波段的反射率。

(3) 熱度。熱度指標(biāo)選用地表溫度代表熱度，Landsat 5 TM影像主要利用熱紅外6波段輻射定標(biāo)參數(shù)，將像元灰度值轉(zhuǎn)換為傳感器處的輻射亮度值，通過Planck輻射函數(shù)求出包含了大氣影響的像元溫度亮度，進(jìn)而通過比輻射率轉(zhuǎn)換為地表溫度[5,18]，提取公式如下：

L6=gain·DN+bias

(6)

(7)

TS=Tb/〔1+(λTb/α)Inε6〕

(8)

式中：gain和bias分別是Landsat 5 TM第6波段的增益和偏置值;K1和K2為定標(biāo)參數(shù):K1=607.76 W/(m2·sr·μm),K2=1 260.56 W/(m2·sr·μm);中心波長(zhǎng)λ=11.48 μm,α取1.438×10-2mK。DN為像元灰度值。

Landsat 8 OLI影像主要利用熱紅外10波段進(jìn)行輻射定標(biāo)，得到熱紅外波段的輻射亮度值，最后通過普朗克定律反函數(shù)求得Ts，其表達(dá)式為：

(9)

(10)

TS=

(11)

(4) 干度。裸土和建筑用地均會(huì)造成地表“干化”，因此選用建筑指數(shù)(IBI)和裸土指數(shù)(SI)構(gòu)建干度指標(biāo)，記為NDSI[19-20]，其計(jì)算公式為：

NDSI=(SI+IBI)/2

(12)

式中：

(13)

(14)

式中：ρBlue,ρGreen,ρRed,ρNIR,ρSWIR1分別為L(zhǎng)andsat 5 TM和Landsat 8 OLI藍(lán)、綠、紅、近紅外、短波紅外1波段的反射率。

1.3.2 RSEIⅠ—Ⅴ級(jí)重心遷移 重心遷移模型是運(yùn)用ArcGIS平臺(tái)繪制出不同時(shí)期、不同等級(jí)的生態(tài)環(huán)境質(zhì)量重心坐標(biāo)，根據(jù)重心坐標(biāo)能夠更加直觀地反映出每個(gè)時(shí)期RSEI不同等級(jí)的遷移方向、遷移的距離等空間動(dòng)態(tài)變化規(guī)律。為了分析1995—2020年神東礦區(qū)不同等級(jí)的RSEI空間變化情況，運(yùn)用ArcGIS 10.4中的工具將5個(gè)等級(jí)進(jìn)行處理，再分別計(jì)算6期的5個(gè)等級(jí)的重心，最后繪制重心移動(dòng)路徑圖。重心坐標(biāo)的運(yùn)算方式如下：

(15)

(16)

2 結(jié)果與分析

2.1 神東礦區(qū)生態(tài)環(huán)境質(zhì)量指標(biāo)主成分分析

利用上述RSEI模型和4個(gè)分量模型的計(jì)算平均值(圖1)，可分析神東礦區(qū)的生態(tài)環(huán)境質(zhì)量變化過程：神東礦區(qū)1995—2020年的生態(tài)環(huán)境質(zhì)量呈現(xiàn)出下降—上升—下降—上升的變化趨勢(shì)；其中，2000年生態(tài)環(huán)境質(zhì)量最低，較1995年生態(tài)環(huán)境質(zhì)量下降11%，2010年生態(tài)環(huán)境質(zhì)量最高，較1995年生態(tài)環(huán)境質(zhì)量提高了5.4%。各分量指標(biāo)(歸一化后)的平均值變化情況: 1995—2020年，綠度指標(biāo)NDVI與RSEI的趨勢(shì)成正比，總體上呈上升趨勢(shì)；干度指標(biāo)NDSI與RSEI的趨勢(shì)成反比，總體上呈下降趨勢(shì)；濕度指標(biāo)WET呈現(xiàn)下降—上升—下降趨勢(shì)，2000年WET出現(xiàn)低值為0.46，總體上呈下降趨勢(shì)；熱度指標(biāo)LST呈現(xiàn)上升—下降—上升—下降—上升的波動(dòng)趨勢(shì)，2005年LST最低值為0.44，總體上呈上升趨勢(shì)。在各指標(biāo)變化中，WET與LST，NDVI與NDSI變化趨勢(shì)成反比?？偟膩碚f，神東礦區(qū)的生態(tài)環(huán)境質(zhì)量主要與NDVI，NDSI指標(biāo)相關(guān)，RSEI變化雖有起伏，但波動(dòng)范圍仍保持在一般狀態(tài)下，生態(tài)環(huán)境質(zhì)量有所提升。

注：NDVI為綠度指標(biāo)；RSEI為遙感生態(tài)指數(shù)；WET為濕度指標(biāo)；LST為熱度指標(biāo)；NDSI為干度指標(biāo)。下同。

利用ArcGIS 10.4在研究區(qū)內(nèi)隨機(jī)生成了180個(gè)點(diǎn)，提取各指標(biāo)及RSEI的值。通過分析1995—2020年NDVI，WET，NDSI和LST與RSEI擬合(圖2)可以發(fā)現(xiàn)：擬合結(jié)果在0.05水平下，均顯著不同于0。NDVI，NDSI，LST與RSEI的R2值分別為0.21，0.29，0.24，WET與RSEI的R2值為0.02，說明NDVI，NDSI，LST是影響RSEI的主要指標(biāo)。通過對(duì)NDVI，WET，NDSI和LST的主成分分析，得到如表1所示結(jié)果。由表1可知，第1主成分的貢獻(xiàn)率均達(dá)到50%以上，表明PC1已集中了4個(gè)指標(biāo)的大部分特征。在1995—2020年，PC1的NDVI均為正值，WET多為正值，說明NDVI，WET與RSEI變化呈正相關(guān)關(guān)系，LST均為負(fù)值，NDSI多為負(fù)值，說明 LST，NDSI與RSEI變化呈負(fù)相關(guān)關(guān)系。相較4個(gè)指標(biāo)對(duì)神東礦區(qū)生態(tài)環(huán)境質(zhì)量貢獻(xiàn)率，NDVI，NDSI對(duì)神東礦區(qū)生態(tài)環(huán)境質(zhì)量影響更大，表明植被生長(zhǎng)分布情況以及人類活動(dòng)是研究區(qū)生態(tài)環(huán)境質(zhì)量的直接影響因素。

圖2 各指標(biāo)環(huán)境生態(tài)質(zhì)量與RSEI擬合

表1 1995-2020年遙感生態(tài)指數(shù)的主成分分析

2.2 神東礦區(qū)生態(tài)環(huán)境質(zhì)量分級(jí)

經(jīng)過歸一化后的神東礦區(qū)RSEI數(shù)值范圍在[0,1]之間，RSEI值越接近于1表示神東礦區(qū)生態(tài)環(huán)境質(zhì)量越好，反之神東礦區(qū)生態(tài)環(huán)境質(zhì)量則越差。根據(jù)神東礦區(qū)RSEI值的范圍，選用前人分類標(biāo)準(zhǔn)，進(jìn)一步將各年份的RSEI以0.2為間隔分成5級(jí)(圖3)，分別代表差、較差、中等、良和優(yōu)5個(gè)等級(jí)，對(duì)應(yīng)的等級(jí)編號(hào)分別為Ⅰ(0～0.2)，Ⅱ(0.2～0.4)，Ⅲ(0.4～0.6)，Ⅳ(0.6～0.8)，Ⅴ(0.8～1.0)。

圖3 1995-2020年神東礦區(qū)遙感生態(tài)指數(shù)(RSEI)等級(jí)分布

通過RSEI分類圖及統(tǒng)計(jì)不同等級(jí)面積所占礦區(qū)總面積的比例可知(圖4)，1995—2020年神東礦區(qū)的生態(tài)環(huán)境質(zhì)量整體呈上升趨勢(shì)，2010年生態(tài)環(huán)境質(zhì)量最好，隨后由于經(jīng)濟(jì)發(fā)展及人類活動(dòng)生態(tài)環(huán)境質(zhì)量有所下降，但整體上呈變好趨勢(shì)。神東礦區(qū)生態(tài)環(huán)境質(zhì)量以“較差”和“中等”為主，合計(jì)平均占比為57%；“優(yōu)”“良”級(jí)合計(jì)平均占比為31.80%，“差”級(jí)平均占比約為11.19%。

圖4 各級(jí)別土地面積占神東礦區(qū)總面積的比例

1995—2010年神東礦區(qū)的生態(tài)環(huán)境質(zhì)量總體趨于穩(wěn)定，并保持繼續(xù)向生態(tài)環(huán)境質(zhì)量好的方向發(fā)展，其中2010年生態(tài)環(huán)境狀況最好，“中等”級(jí)面積變化幅度較小，相較1995年的“良”、“優(yōu)”等級(jí)之和2010年提升了33.77%。2010—2015年生態(tài)環(huán)境質(zhì)量有所下降，2015—2020年生態(tài)環(huán)境質(zhì)量有所提升。1995年生態(tài)環(huán)境質(zhì)量為“差”和“較差”的面積分別占整個(gè)研究區(qū)域的14.45%和28.51%，而2020年分別為9.34%和21.60%，“差”、“較差”級(jí)面積明顯下降，“中等”、“良”、“優(yōu)”等級(jí)面積均有提升，相較于1995年提升了12.01%。

通過研究1995—2020年RSEIⅠ—Ⅴ等級(jí)的土地利用類型及轉(zhuǎn)移矩陣(表2)發(fā)現(xiàn)：1995—2020年神東礦區(qū)生態(tài)環(huán)境質(zhì)量各等級(jí)主要土地利用類型為草地、耕地。神東礦區(qū)生態(tài)環(huán)境質(zhì)量整體上有上升趨勢(shì)，且耕地、草地面積整體上也在增加。說明神東礦區(qū)生態(tài)環(huán)境質(zhì)量受耕地、草地面積變化影響較大。

表2 1995-2020年神東礦區(qū)土地利用類型轉(zhuǎn)移矩陣 %

在1995—2020年，西南方向土地類型為內(nèi)陸湖泊，湖心有島嶼。按照模型因子，該地區(qū)主要評(píng)價(jià)指標(biāo)為濕度指標(biāo)(WET)和干度指標(biāo)(NDSI)。通過主成分分析的結(jié)果可以看到NDSI在2010年貢獻(xiàn)程度最大，且為正值，而在其他年份NDSI貢獻(xiàn)程度較小，且為負(fù)值，故在2010年內(nèi)陸湖泊質(zhì)量相對(duì)最差。

通過1995年與2020年生態(tài)環(huán)境質(zhì)量等級(jí)轉(zhuǎn)移變化分析，發(fā)現(xiàn)神東礦區(qū)生態(tài)環(huán)境質(zhì)量變化的區(qū)域較為分散，生態(tài)環(huán)境質(zhì)量變差的區(qū)域占總面積的23.09%，生態(tài)環(huán)境質(zhì)量變好的區(qū)域占總面積的42.22%(圖5)，整體呈變好趨勢(shì)。

圖5 1995-2020年神東生態(tài)質(zhì)量變化特征

2.3 不同等級(jí)RSEI重心轉(zhuǎn)移

通過分析1995—2020年神東礦區(qū)生態(tài)環(huán)境質(zhì)量各等級(jí)重心遷移情況(圖6)可知：神東礦區(qū)生態(tài)環(huán)境質(zhì)量Ⅰ—Ⅴ級(jí)重心在伊金霍洛旗和神木縣范圍內(nèi)進(jìn)行遷移。在1995—2020年“差”、“較差”、“優(yōu)”等級(jí)的重心轉(zhuǎn)移年幅度較大，其他等級(jí)的重心移動(dòng)幅度均較小。其中“差”等級(jí)的重心轉(zhuǎn)移在2000—2015年遷移均值達(dá)到16.32 km；“較差”等級(jí)的重心轉(zhuǎn)移在2005—2015年遷移均值達(dá)到11.13 km；“優(yōu)”等級(jí)的重心轉(zhuǎn)移在2005—2020年遷移均值達(dá)到14.02 km；其他等級(jí)的重心遷移均值為3.11 km。說明1995年—2020年，神東礦區(qū)不同生態(tài)環(huán)境質(zhì)量等級(jí)的區(qū)域變化明顯，波動(dòng)范圍較大。從5個(gè)等級(jí)的重心移動(dòng)路徑可以看出，“差”等級(jí)向西南方向偏移，“較差”、“中等”、“良”和“優(yōu)”等級(jí)向西北方向偏移。

圖6 1995-2020神東礦區(qū)RSEI Ⅰ-Ⅴ級(jí)移動(dòng)路徑

3 討 論

神東礦區(qū)作為中國重要的區(qū)域生態(tài)脆弱帶和重點(diǎn)治理區(qū)，其生態(tài)環(huán)境質(zhì)量的影響不言而喻。當(dāng)前，對(duì)礦區(qū)生態(tài)環(huán)境質(zhì)量評(píng)價(jià)多以綜合指數(shù)法、層次分析等方法進(jìn)行研究[21-24]，但就本研究區(qū)域而言，神東礦區(qū)是自然環(huán)境較多和人居較少的研究區(qū)，RSEI指數(shù)適用性比較強(qiáng)。而運(yùn)用RSEI指數(shù)研究市縣、流域、自然遺產(chǎn)、濕地等研究區(qū)較多[25-27]，對(duì)礦區(qū)研究較少。

已有文獻(xiàn)[28-29]運(yùn)用RSEI指數(shù)對(duì)礦區(qū)生態(tài)環(huán)境質(zhì)量變化進(jìn)行研究，但是沒考慮生態(tài)環(huán)境質(zhì)量變化是具有空間區(qū)域性的，本文增加了生態(tài)環(huán)境質(zhì)量在空間方面變化的研究；并在土地利用類型與生態(tài)環(huán)境質(zhì)量變化上做了進(jìn)一步的分析，對(duì)神東礦區(qū)生態(tài)環(huán)境質(zhì)量變化做了更深的研究。岳輝等[29]對(duì)神東礦區(qū)生態(tài)環(huán)境質(zhì)量變化的研究，得出1989—2015年神東礦區(qū)的生態(tài)環(huán)境質(zhì)量整體呈上升趨勢(shì)。這與其研究結(jié)果綜合可以得出，1989—2020年神東礦區(qū)整體生態(tài)環(huán)境質(zhì)量總體上呈向好發(fā)展的趨勢(shì)。礦區(qū)生態(tài)環(huán)境質(zhì)量與人類活動(dòng)強(qiáng)度、綠色植被的關(guān)系比較密切，綠色植被生長(zhǎng)情況主要決定了綠度指標(biāo)的高低，干度指標(biāo)很大程度受人為開發(fā)活動(dòng)影響，在二者綜合影響下決定了區(qū)域生態(tài)環(huán)境質(zhì)量的總體水平。在礦區(qū)生態(tài)環(huán)境質(zhì)量分級(jí)分析中，1995—2020年各級(jí)土地利用類型均以草地、耕地為主，且與礦區(qū)生態(tài)環(huán)境質(zhì)量總體變化趨勢(shì)一致。Ⅱ級(jí)、Ⅲ級(jí)占礦區(qū)50%以上的面積，且遷移幅度較小，說明在生態(tài)環(huán)境質(zhì)量治理過程中，大部分的區(qū)域生態(tài)環(huán)境質(zhì)量得到了一定的維持，加大生態(tài)環(huán)境保護(hù)力度。Ⅰ級(jí)、Ⅳ級(jí)和Ⅴ級(jí)的生態(tài)環(huán)境質(zhì)量遷移幅度較大，說明在較差的區(qū)域加大了生態(tài)環(huán)境的治理力度，導(dǎo)致較差較好的區(qū)域變化幅度較大。

在以后的研究中，可以豐富生態(tài)環(huán)境質(zhì)量評(píng)價(jià)指標(biāo)，對(duì)神東礦區(qū)生態(tài)環(huán)境質(zhì)量變化的驅(qū)動(dòng)因子進(jìn)行細(xì)致的分析，以更細(xì)致的討論神東礦區(qū)生態(tài)環(huán)境質(zhì)量變化原因，為對(duì)神東礦區(qū)的生態(tài)環(huán)境的保護(hù)、治理提供依據(jù)。

4 結(jié) 論

(1) 1995—2020年神東礦區(qū)的生態(tài)環(huán)境質(zhì)量總體呈變好趨勢(shì)。綠度指標(biāo)與神東礦區(qū)生態(tài)質(zhì)量變化一致，干度指標(biāo)與神東礦區(qū)生態(tài)質(zhì)量變化相反，綠度指標(biāo)是影響生態(tài)環(huán)境質(zhì)量的主要原因，干度指標(biāo)是影響生態(tài)環(huán)境質(zhì)量的次要原因。

(2) 神東礦區(qū)生態(tài)環(huán)境質(zhì)量變化較為分散，以Ⅱ級(jí)、Ⅲ級(jí)為主，其面積接近于研究區(qū)總面積的57%。對(duì)比1995年礦區(qū)生態(tài)環(huán)境質(zhì)量，2018年變好的區(qū)域是研究區(qū)總面積的42.22%，神東礦區(qū)生態(tài)環(huán)境質(zhì)量有所提升。

(3) 在1995—2020年神東礦區(qū)生態(tài)環(huán)境質(zhì)量各級(jí)土地利用類型均以草地、耕地為主，礦區(qū)生態(tài)環(huán)境質(zhì)量與耕地、草地面積均呈現(xiàn)總體上升趨勢(shì)。

(4) 神東礦區(qū)RSEIⅠ—Ⅴ級(jí)的重心轉(zhuǎn)移幅度較大等級(jí)有：Ⅰ級(jí)、Ⅳ級(jí)和Ⅴ級(jí)，遷移距離在10～15 km之間，說明處于該區(qū)域的神東礦區(qū)生態(tài)質(zhì)量空間遷移幅度較大。重心轉(zhuǎn)移幅度較小的等級(jí)有Ⅱ級(jí)、Ⅲ級(jí)，遷移距離在3 km左右，說明處于該區(qū)域的神東礦區(qū)生態(tài)質(zhì)量空間遷移幅度相對(duì)穩(wěn)定。