張永庭，陳建平，劉自增，閆 軍，黃 濤，嚴 瑾，樊新剛

(1.中國地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院，北京 100083；2.寧夏遙感測繪勘查院，寧夏 銀川 750021；3.寧夏大學(xué)社會服務(wù)處，寧夏 銀川 750021)

寧夏礦山環(huán)境質(zhì)量空間分布特征及其分區(qū)

張永庭1，2，陳建平1，劉自增2，閆 軍2，黃 濤2，嚴 瑾2，樊新剛3

(1.中國地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院，北京 100083；2.寧夏遙感測繪勘查院，寧夏 銀川 750021；3.寧夏大學(xué)社會服務(wù)處，寧夏 銀川 750021)

為明確寧夏礦山環(huán)境質(zhì)量空間分布特征和問題區(qū)域，提供修復(fù)和管理的依據(jù)，以全區(qū)22個縣區(qū)為研究單元，構(gòu)建了土地占壓、地質(zhì)災(zāi)害、環(huán)境污染、水土流失四個維度的準(zhǔn)則及對應(yīng)的共16項指標(biāo)的礦山環(huán)境綜合評價指標(biāo)體系，利用遙感解譯與礦山環(huán)境統(tǒng)計資料獲取數(shù)據(jù)，基于熵權(quán)模糊綜合方法評價，利用Spss22.0進行模糊聚類分組，結(jié)合實地踏勘確定分級含義后分區(qū)。結(jié)果顯示：寧夏礦山環(huán)境質(zhì)量空間差異較大，中北部礦山環(huán)境質(zhì)量總體優(yōu)于南部山區(qū)、但北部局部縣區(qū)問題突出，主要環(huán)境問題在北部老舊礦區(qū)集中的惠農(nóng)區(qū)和大武口區(qū)、新興礦區(qū)靈武市(縣級市)，以及中部的沙坡頭區(qū)。全區(qū)劃分了“嚴重區(qū)、較為嚴重區(qū)、一般區(qū)、良好區(qū)”四類礦山環(huán)境區(qū)，存在“面上連續(xù)、點上鑲嵌”的特征。同時，應(yīng)關(guān)注新興礦山靈武市的開采進度，持續(xù)監(jiān)測和預(yù)防，以最大限度降低環(huán)境危害。利用聚類分析得出的分級閾值存在不連續(xù)性，推廣應(yīng)用的前提是應(yīng)加大案例樣本以找出閾值分布的規(guī)律，從而實現(xiàn)連續(xù)分級。

礦山環(huán)境質(zhì)量；空間分區(qū)；熵權(quán)；模糊綜合評價；寧夏

寧夏是我國西北礦產(chǎn)資源比較豐富的省區(qū)之一，已發(fā)現(xiàn)各類礦產(chǎn)46種，煤炭、石膏、冶金用砂巖、水泥配料用板巖、瓦用黏土等能源礦產(chǎn)和非金屬礦產(chǎn)儲量居于全國前列[1]。目前，自治區(qū)內(nèi)以能源(煤)、非金屬(石灰?guī)r、硅石及各類建材礦產(chǎn))等礦山為主，形成了以石嘴山市、靈武市(縣級市)等資源型城市為代表的礦山集中區(qū)，以及煤炭為支柱產(chǎn)業(yè)的經(jīng)濟發(fā)展格局。由于經(jīng)濟欠發(fā)達背景下“礦山開采的經(jīng)濟依賴、礦產(chǎn)資源深加工能力不足，生態(tài)環(huán)境和地質(zhì)環(huán)境治理投入不足”等原因，長期開發(fā)積累了較為嚴重的礦山環(huán)境和地質(zhì)災(zāi)害問題。截至2015年底，全區(qū)22個縣區(qū)共開采礦山849處，全年發(fā)生礦山地質(zhì)災(zāi)害達136處，其中崩塌11處，滑坡15處，泥石流11處，地面塌陷及地裂縫84處[2]。因此，調(diào)查評價寧夏礦山環(huán)境質(zhì)量，分析空間分布特征，明確主要問題區(qū)，對指導(dǎo)國土與環(huán)保部門的科學(xué)治理具有重要的現(xiàn)實意義。同時，礦山環(huán)境質(zhì)量作為“綠色礦山”建設(shè)的主要指標(biāo)和區(qū)域可持續(xù)發(fā)展的重要領(lǐng)域，日益得到學(xué)界關(guān)注[3-4]，但目前我國礦山環(huán)境評價多數(shù)圍繞單個礦山或礦區(qū)開展，區(qū)域礦山綜合環(huán)境質(zhì)量評價正在探索之中，評價指標(biāo)體系多樣、龐雜，影響了結(jié)果的合理性和應(yīng)用指導(dǎo)性[5-7]，本研究將礦山地質(zhì)環(huán)境和生態(tài)環(huán)境統(tǒng)一納入評價指標(biāo)體系，可為礦山綜合環(huán)境評價方法研究提供借鑒案例。

1 寧夏礦山開發(fā)及環(huán)境問題空間分布

據(jù)我國大地構(gòu)造綜合地層分區(qū)，寧夏屬華北地層大區(qū)，境內(nèi)成礦由北向南分為：鄂爾多斯西緣Fe-Pb-Zn-磷-石膏-芒硝成礦帶、鄂爾多斯(盆地)U-石油-天然氣-煤-鹽類成礦區(qū)、河西走廊 Fe-Mn-螢石-鹽類-石油成礦帶、北祁連 Cu-Pb-Zn-Fe-Cr-Au-Ag-硫鐵礦-石棉成礦帶[8]。目前，礦山開發(fā)空間格局為(圖1)：無煙煤、煉焦用煤、硅石在賀蘭山北段；非煉焦用煤、石油(氣)在靈武市、鹽池縣；石膏在海原縣東北部、鹽池縣中部及中衛(wèi)市西部；石灰?guī)r大礦集中在賀蘭山中段、牛首山、天景山、青龍山等地；白云巖在青龍山和賀蘭山中段；石英砂在青銅峽市和原州區(qū)；金屬礦產(chǎn)主要產(chǎn)出于賀蘭山、衛(wèi)寧北山和南-西華山地區(qū)[9]。寧夏地處干旱與半干旱地區(qū)，降水量小、植被覆蓋低，境內(nèi)多數(shù)礦山巖層裸露，主要環(huán)境問題為地面塌陷和裂縫、崩塌滑坡、泥石流等地質(zhì)災(zāi)害，礦山占地破壞植被、土壤和水污染，以及煤矸石自燃的SO2環(huán)境空氣污染等生態(tài)環(huán)境問題。境內(nèi)能源礦產(chǎn)(煤)、非金屬礦產(chǎn)(石灰?guī)r、硅石及各類建材礦產(chǎn))等礦山企業(yè)較多，產(chǎn)生礦山環(huán)境地質(zhì)問題最多的礦種是煤礦，其次是灰?guī)r礦、硅石礦、建筑用砂和建筑用石料[2]。

圖1 寧夏礦山開發(fā)現(xiàn)狀分布

2 礦山環(huán)境質(zhì)量評價及空間分區(qū)方法

2.1 礦山環(huán)境質(zhì)量評價方法

2.1.1 評價指標(biāo)體系構(gòu)建

評價指標(biāo)體系構(gòu)建是礦山環(huán)境質(zhì)量評價的關(guān)鍵環(huán)節(jié)。在借鑒甘肅[10]、貴州[11]、安徽[6]等指標(biāo)體系基礎(chǔ)上，綜合陳橋等[12]、徐友寧等[13]對礦山綜合環(huán)境評價指標(biāo)體系構(gòu)建的專項成果，以地質(zhì)災(zāi)害為核心，地表直觀反映為途徑，基于系統(tǒng)性、代表性、可量化性、直觀性原則，設(shè)計了土地占壓(B1)、地質(zhì)災(zāi)害(B2)、環(huán)境污染(B3)、水土流失(B4)四個維度的準(zhǔn)則及對應(yīng)的16項指標(biāo)(表1)，以反映礦山環(huán)境質(zhì)量水平。其中，土地占壓利用各縣區(qū)“礦山、固廢堆積、開采和中轉(zhuǎn)場地”的面積4項指標(biāo)表達，反映區(qū)域礦山開采對土地和植被資源的占用規(guī)模及內(nèi)部結(jié)構(gòu)水平。地質(zhì)災(zāi)害方面，利用各縣區(qū)的塌陷區(qū)數(shù)量和面積表達面上問題的嚴重程度，點上設(shè)計了滑坡、崩塌、泥石流3個指標(biāo)表達縣區(qū)重點環(huán)境問題。環(huán)境污染選擇在土壤和水體污染面積指標(biāo)的基礎(chǔ)上，由于寧夏部分煤礦區(qū)如石嘴山市、靈武市煤矸石含硫量高，其自燃和露天堆放對于大氣環(huán)境影響較大，選擇環(huán)境空氣中SO2含量表達。寧夏屬干旱半干旱區(qū)域，地表植被稀疏，表土和巖層裸露，因此，選擇了降水量、最大坡度(礦區(qū)范圍內(nèi))、侵蝕力、植被覆蓋度4項指標(biāo)表達地表破壞下的水土流失程度。研究表明，觀察指標(biāo)的增加本來是為了使認識目標(biāo)趨于完整，但反過來，為使結(jié)果清晰而一味增加觀察指標(biāo)又讓評價陷入混亂[14]。因此，控制指標(biāo)規(guī)模，選擇的16項指標(biāo)全部可量化，力求直接和簡化，避免過度關(guān)注系統(tǒng)性和全面性而破壞可操作性和代表性等問題。

表1 寧夏礦山環(huán)境質(zhì)量評價指標(biāo)體系

2.1.2 熵權(quán)模糊綜合評價方法

對于復(fù)雜系統(tǒng)的綜合評價，涉及因素眾多、層次結(jié)構(gòu)復(fù)雜且包含著各種不確定性、模糊性和隨機性，因此對多層次因素的綜合以及對各類不確定性的處理成為復(fù)雜系統(tǒng)綜合評價的關(guān)鍵所在[15]。熵權(quán)法是把各評價指標(biāo)的信息進行量化與綜合，計算各指標(biāo)反映的信息熵，通過信息熵來確定權(quán)重的客觀賦權(quán)方法[16]。熵權(quán)法有效地避免了人為因素的干擾，使評價結(jié)果更符合實際，從而給出客觀可靠的評價結(jié)果。與主觀賦權(quán)法相比，通過客觀賦權(quán)法得出的權(quán)重是以實際數(shù)據(jù)為依據(jù)，對數(shù)據(jù)本身進行挖掘，從而找出事物間的內(nèi)在規(guī)律，更具客觀性[17]。一般來說，若某個指標(biāo)的信息熵越小，表明指標(biāo)值得變異程度越大，提供的信息量越多，在綜合評價中所能起到的作用也越大，其權(quán)重也就越大。相反，某個指標(biāo)的信息熵越大，表明指標(biāo)值得變異程度越小，提供的信息量也越少，在綜合評價中所起到的作用也越小，其權(quán)重也就越小[18]?；陟貦?quán)模糊綜合評價是利用熵權(quán)替代隸屬度進行的模糊綜合評價，利用熵權(quán)權(quán)重的客觀性與模糊評價的主觀性相整合，實現(xiàn)了二者的統(tǒng)一，更為科學(xué)地進行目標(biāo)評價[19]。熵權(quán)模糊綜合評價步驟如下所述。

1)原始數(shù)據(jù)矩陣及其標(biāo)準(zhǔn)化。設(shè)指標(biāo)體系有m個評價指標(biāo)、n個評價對象，得到的原始數(shù)據(jù)矩陣見式(1)。

(1)

之后，在對m個評價指標(biāo)性質(zhì)判斷下，分為正向指標(biāo)和負向指標(biāo)，分別采用rij=(xij-xmin)/(xmax-xmin)、rij=(xmax-xij)/(xmax-xmin)對該矩陣標(biāo)準(zhǔn)化，形成無量綱矩陣R(式(2))。正向指標(biāo)值越大，則環(huán)境質(zhì)量越好；負向指標(biāo)值越大，則環(huán)境質(zhì)量越差。

(2)

式中：Rij為第j個評價對象第i個指標(biāo)的標(biāo)準(zhǔn)化值；xij為第j個評價對象第i個指標(biāo)的原始值；xmax為相同指標(biāo)下不同評價對象中的最大值；xmin為相同指標(biāo)下不同評價對象中的最小值。

2)定義熵。在m個指標(biāo)，n個被評價對象的評價問題中，第i 個指標(biāo)的熵值定義見式(3)。

(3)

3)定義熵權(quán)。定義了第i 個指標(biāo)的熵值之后，可得到第i 個指標(biāo)的熵權(quán)值，見式(4)。

(4)

4)計算模糊綜合評價結(jié)果。對于權(quán)重A={ω1，ω2，…ωn}與模糊關(guān)系矩陣R，計算模糊評價結(jié)果，見式(5)。

(5)

式中：“°”為模糊算子；Z為綜合評價值。

2.2 環(huán)境質(zhì)量空間分區(qū)方法

目前研究礦山環(huán)境質(zhì)量分級的方法主要分為“以劃分評價指標(biāo)目標(biāo)值等級、專家意見分級”兩類，存在指標(biāo)值確定合理性和分級閾值主觀性等問題，對研究區(qū)礦山環(huán)境的實際狀況的結(jié)合程度不夠，實踐指導(dǎo)性不強[20-21]。事實上，由于多指標(biāo)體系方法是通過指標(biāo)間相互重要性或?qū)δ繕?biāo)貢獻程度逐步匯總得出的值，且必須通過原始數(shù)據(jù)的無量綱化，因此，最終評價值屬相對量，即最終確定的分級等次是以多目標(biāo)中的最大量和最小量為參照基準(zhǔn)的，屬相對結(jié)果，基于熵權(quán)模糊綜合評價方法也不例外，必須進行驗證。鑒于此，基于2010年、2013年、2014年的礦山環(huán)境質(zhì)量指數(shù)值，利用Spss22.0軟件對寧夏22個縣區(qū)進行系統(tǒng)聚類分組，進而確定分段閾值。同時，通過對閾值所在縣區(qū)礦山環(huán)境實地踏勘，驗證劃分的合理性后，確定分級名稱，最終利用ArcGIS10.0制圖形成分區(qū)結(jié)果。

2.3 數(shù)據(jù)來源及處理

土地占壓、地質(zhì)災(zāi)害、環(huán)境污染，以及最大坡度和植被覆蓋度的數(shù)據(jù)利用2010年、2013年、2014年三個時段的GF2、GF1、02C、ZY3、SJ9、YG14和P1等衛(wèi)星圖像提取，總面積5.18萬km2，空間分辨率為0.5～2m，經(jīng)平面坐標(biāo)體系和高程基準(zhǔn)轉(zhuǎn)換，幾何校正，圖像鑲嵌和增強，在ArcGIS10.0平臺進行信息提取、屬性錄入。2015年7～9月，進行了外業(yè)勘查，驗證信息提取的可靠性，共完成中轉(zhuǎn)場地、固體廢棄物、尾礦庫、礦山建筑等占地斑塊野外驗證點410個，礦山地質(zhì)災(zāi)害和環(huán)境污染斑塊野外驗證點33個。降水量數(shù)據(jù)主要來自《寧夏回族自治區(qū)資源環(huán)境地圖集》[22]。侵蝕力數(shù)據(jù)(侵蝕力是由降雨引起土壤侵蝕的潛在能力，侵蝕力具有明顯的空間分異特征和年際變化。)基于“中國多年平均降雨侵蝕力等值線圖”[23]，利用ArcGIS10.0軟件空間內(nèi)插提取。各縣區(qū)環(huán)境空氣SO2含量數(shù)據(jù)來源于銀川市、石嘴山市、吳忠市、固原市、中衛(wèi)市的環(huán)境質(zhì)量年報。

3 結(jié)果及分析

3.1 環(huán)境質(zhì)量現(xiàn)狀分布

圖2 寧夏22個縣區(qū)礦山環(huán)境質(zhì)量綜合指數(shù)值

寧夏22個縣區(qū)2010年、2013年、2014年的平均礦山環(huán)境質(zhì)量綜合指數(shù)值差異較大(圖2)，呈“中北部整體高位且高低波峰共存、南部偏低且平緩”特征，即中北部礦山環(huán)境質(zhì)量總體優(yōu)于南部山區(qū)、但北部局部縣區(qū)問題突出。首先，北部的惠農(nóng)區(qū)、大武口區(qū)及靈武市、中部的沙坡頭區(qū)礦山環(huán)境質(zhì)量最差，分別為0.4729、0.5003、0.5320、0.4893，是寧夏礦山治理的重點區(qū)域。經(jīng)實地踏勘，大武口區(qū)屬寧夏開發(fā)歷史最久(超過60年)、規(guī)模最大的礦區(qū)，區(qū)內(nèi)沿賀蘭山分布有石炭井、汝芨溝等老舊礦區(qū)，煤矸石等廢棄物堆積達5697.01hm2，占礦區(qū)總面積的82.02%，地表破壞嚴重，由于坡度較大，水土流失風(fēng)險大，導(dǎo)致了以上結(jié)果?；蒉r(nóng)區(qū)與大武口區(qū)同屬煤炭資源型城市石嘴山市，開發(fā)歷史久，有石嘴山礦務(wù)局一二三礦和沙巴臺煤礦等大中型礦山，雖然進行了塌陷區(qū)治理(治理的圖斑面積910hm2)和植被恢復(fù)[2]，但仍存在固廢占地(1605.33hm2)和中轉(zhuǎn)場地占地(1256.09hm2)面積大(分別占礦區(qū)總面積的46.78%、36.60%)、環(huán)境空氣中SO2含量高等問題，礦山環(huán)境質(zhì)量低。銀川市下轄的靈武市屬新興礦區(qū)，以井工開采為主，已經(jīng)形成了4處塌陷區(qū)(面積698.1hm2)，中轉(zhuǎn)場地和采場面積大(占礦區(qū)總面積的45.82%)，存在506.86hm2的土壤和水體污染，總體環(huán)境質(zhì)量較差。沙坡頭區(qū)地處騰格里沙漠邊緣，植被覆蓋度不足30%，礦山泥石流和崩塌災(zāi)害達6處，固體廢棄物占地(874.69hm2)為礦區(qū)總面積的52.07%，是造成環(huán)境質(zhì)量值低的主要原因。其次，介于0.6～0.7的為寧南山區(qū)的海原縣、原州區(qū)、西吉縣、彭陽縣、隆德縣、涇源縣和中部吳忠市下轄的青銅峽縣級市，整體屬于礦山環(huán)境相對較差區(qū)域。寧南山區(qū)6縣區(qū)礦山開采總體規(guī)模較小(總面積2273.13hm2)，環(huán)境壓力較小，但由于地處六盤山和黃土高原丘陵區(qū)，坡度較大(15～30°)，降雨量大(367.4～619.9mm)，特別是侵蝕力大(>1000MJ·mm·hm-2·h-1· a-1)，使其礦山環(huán)境質(zhì)量總體上偏低，惡化風(fēng)險較大。其中，彭陽縣分布著寧南山區(qū)最大規(guī)模的王洼煤礦和長慶油田采區(qū)，已經(jīng)出現(xiàn)了32.8hm2的塌陷區(qū)、4處滑坡，使其成為7縣市中環(huán)境質(zhì)量最低的縣。青銅峽市由于礦山地表占地破壞相對較小，廢棄物堆積828.32hm2，占礦區(qū)總面積的28.61%，但土壤和水污染土地面積達525.25hm2，占礦區(qū)總面積的18.15%，礦山衍生環(huán)境破壞壓力較為突出。最后，其余11個縣區(qū)質(zhì)量值均高于0.7，其中，北部石嘴山市平羅縣、銀川市興慶區(qū)和金鳳區(qū)、吳忠市利通區(qū)礦山環(huán)境質(zhì)量指數(shù)分別達到了0.9150、0.9505、0.9537、0.8669，礦山環(huán)境質(zhì)量最優(yōu)。

3.2 環(huán)境質(zhì)量空間分區(qū)

圖3 寧夏22個縣區(qū)礦山環(huán)境質(zhì)量系統(tǒng)聚類

利用寧夏22個縣區(qū)2010年、2013年、2014年的礦山環(huán)境質(zhì)量綜合指數(shù)值，基于離差平方和法聚類(圖3)，礦山環(huán)境質(zhì)量可分為“寧南山區(qū)6縣區(qū)+靈武市、惠農(nóng)區(qū)+大武口區(qū)+沙坡頭區(qū)、興慶區(qū)+金鳳區(qū)+利通區(qū)+平羅縣、永寧縣+中寧縣+紅寺堡區(qū)+鹽池縣+同心縣+賀蘭縣+西夏區(qū)+青銅峽市”4組，結(jié)合前文分析，靈武市礦山環(huán)境問題嚴重且有加重趨勢，將其與寧南6縣區(qū)置于同一分組與實際存在差距，因此將其放入“惠農(nóng)區(qū)+大武口區(qū)+沙坡頭區(qū)”。青銅峽市2010年、2013年礦山環(huán)境質(zhì)量綜合指數(shù)值分別為0.7349、0.7043，因此得到了以上聚類分析結(jié)果。據(jù)此，按照2010年至2014年礦山環(huán)境綜合指數(shù)平均值確定閾值范圍：≤0.5320(惠農(nóng)區(qū)+大武口區(qū)+沙坡頭區(qū)+靈武市)、0.6477～0.6859(寧南山區(qū)6縣區(qū))、0.7239～0.8110(永寧縣+中寧縣+紅寺堡區(qū)+鹽池縣+同心縣+賀蘭縣+西夏區(qū)+青銅峽市)、≥0.8669(興慶區(qū)+金鳳區(qū)+利通區(qū)+平羅縣)?？梢?，以上閾值存在不連續(xù)性，但作為指導(dǎo)礦山環(huán)境修復(fù)管理的目的，不影響其結(jié)果。楊金中等對全國各省區(qū)礦山環(huán)境遙感評價分級中，礦山數(shù)量規(guī)模小的省區(qū)內(nèi)，區(qū)域間礦山環(huán)境質(zhì)量差異較大[24]，印證了這一結(jié)果。

圖4 寧夏礦山環(huán)境質(zhì)量空間分區(qū)

通過實地踏勘，依據(jù)以上閾值，各組分別按順序定義為環(huán)境質(zhì)量“嚴重區(qū)、較為嚴重區(qū)、一般區(qū)、良好區(qū)”。其中，嚴重區(qū)按照評價指標(biāo)體系準(zhǔn)則層含義，表達為土地占壓面積大且廢棄物占地面積大，土壤和水及環(huán)境空氣污染嚴重，地面塌陷、崩塌和泥石流等地質(zhì)災(zāi)害嚴重，水土流失壓力大。剩余級別含義依據(jù)名稱以此類推。2014年，寧夏礦山開發(fā)占地總面積30433.06hm2(約占寧夏總面積5.18萬km2的0.59%)。其中，嚴重區(qū)面積17756.9hm2，占礦山開發(fā)總面積的58.35%；較為嚴重區(qū)面積2273.13hm2，占7.47%；一般區(qū)面積8880.29hm2，占29.18%；良好區(qū)面積1522.74hm2，占5.04%。結(jié)合空間分布可見(圖4)，嚴重區(qū)、一般區(qū)和良好區(qū)的礦山全部分布在寧夏平原經(jīng)濟發(fā)達的中北部區(qū)域，表明臨近城市產(chǎn)業(yè)區(qū)的礦山開發(fā)規(guī)模較大，產(chǎn)生了局部環(huán)境嚴重破壞區(qū)域，未來應(yīng)重點監(jiān)控嚴重區(qū)的礦山開采活動。文獻[2]的研究也表明“寧夏以煤炭產(chǎn)業(yè)為主的惠農(nóng)區(qū)、大武口區(qū)、靈武市，以及中部的中衛(wèi)市沙坡頭區(qū)，由于開發(fā)規(guī)模大，地面塌陷、地裂縫、泥石流、地表植被破壞和土水污染問題突出?！保c本研究結(jié)果一致。從全區(qū)看，在“北密南疏”基礎(chǔ)上，礦山環(huán)境質(zhì)量空間上呈“面上從北至南連續(xù)分布、點上四類分區(qū)鑲嵌分布”的特征，尤其是良好區(qū)與嚴重區(qū)夾雜于其他分區(qū)之間。

4 討論及結(jié)論

1)寧夏礦山環(huán)境質(zhì)量總體特征表現(xiàn)為22個縣區(qū)間差異較大，北部優(yōu)于南部山區(qū)、但北部局部縣區(qū)問題突出。礦山環(huán)境問題集中在北部老舊礦區(qū)集中的石嘴山市的惠農(nóng)區(qū)和大武口區(qū)、新興礦區(qū)銀川市的靈武市(縣級市)，以及中部的中衛(wèi)市沙坡頭區(qū)。

2)寧夏全區(qū)可分為“嚴重區(qū)、較為嚴重區(qū)、一般區(qū)、良好區(qū)”四類礦山環(huán)境區(qū)，在空間上呈“面上連續(xù)、點上鑲嵌”的特征，嚴重區(qū)在北部和中部均有分布，寧南山區(qū)6縣區(qū)屬較為嚴重區(qū)，均存在一定礦山環(huán)境風(fēng)險，是礦山環(huán)境管理的重點區(qū)域。

3)新興礦區(qū)靈武市目前其開采量位居寧夏首位，已經(jīng)形成了較大面積的塌陷區(qū)和環(huán)境惡化趨勢。由于其主要為井工開采，早期地表破壞具有隱蔽性，往往在生產(chǎn)一定年限后由量變轉(zhuǎn)化為質(zhì)變，出現(xiàn)更大面積的地面塌陷、水土流失、環(huán)境污染問題。因此，未來應(yīng)關(guān)注與此類似的新興礦山開采進度，持續(xù)監(jiān)測，提前預(yù)防地質(zhì)災(zāi)害和環(huán)境污染，以最大限度降低危害。

4)利用聚類分析得出的寧夏礦山環(huán)境分級閾值存在不連續(xù)性，與以往各類分級方法結(jié)果存在差異，推廣應(yīng)用的前提是應(yīng)加大此種方法的案例樣本，從而找出熵權(quán)模糊綜合評價值分級閾值的分布規(guī)律，實現(xiàn)連續(xù)分級和推廣應(yīng)用。

[1] 宋新華.寧夏礦產(chǎn)資源產(chǎn)業(yè)布局研究[D].北京：中國地質(zhì)大學(xué)(北京)，2013.

[2] 閆軍，劉自增，張永庭，等.寧夏自治區(qū)礦產(chǎn)資源開發(fā)環(huán)境遙感監(jiān)測成果報告[R].寧夏地質(zhì)調(diào)查院，2015:12-15.

[3] 鄒長新，沈渭壽，劉發(fā)民.礦山生態(tài)環(huán)境質(zhì)量評價指標(biāo)體系初探[J].中國礦業(yè)，2011，20(8)：56-59，68.

[4] 亢亢，趙曉林，高建國.可拓聚類預(yù)測在礦山地質(zhì)環(huán)境質(zhì)量預(yù)測評價中的應(yīng)用[J].安全與環(huán)境學(xué)報，2013(2)：111-116.

[5] 任幼蓉，吳慶璞.礦山地質(zhì)災(zāi)害危險性指數(shù)探討[J].地質(zhì)災(zāi)害與環(huán)境保護，2004，15(1)：1-4.

[6] 江松林，孫世群，王輝.安徽省礦山環(huán)境質(zhì)量綜合評價研究[J].合肥工業(yè)大學(xué)學(xué)報：自然科學(xué)版，2008，31(1)：112-115.

[7] 趙曉亮，齊慶杰，趙國智，等.礦山環(huán)境質(zhì)量集對分析定量評價模型的建立與應(yīng)用研究[J].干旱區(qū)資源與環(huán)境，2014，28(11)：72-77.

[8] 梁峰.寧夏回族自治區(qū)鐵礦資源潛力評價報告[R].寧夏地質(zhì)調(diào)查院，2010:59-60.

[9] 陸彥俊，倪萬魁，尚慧，等.寧夏礦山環(huán)境評價整治與監(jiān)測技術(shù)方法研究[R].寧夏國土資源調(diào)查監(jiān)測院，2012.

[10] 王念秦，王永鋒，王得楷.甘肅礦山生態(tài)地質(zhì)環(huán)境現(xiàn)狀綜合評價分區(qū)研究[J].水土保持研究，2009，16(5)：225-228，232.

[11] 羅炳佳，沈誠.貴州礦山地質(zhì)環(huán)境影響評估[J].水文地質(zhì)工程地質(zhì)，2013，40(1)：134-138.

[12] 陳橋，胡克，雒昆利，等.基于AHP法的礦山生態(tài)環(huán)境綜合評價模式研究[J].中國礦業(yè)大學(xué)學(xué)報，2006，35(3)：377-383.

[13] 徐友寧，袁漢春，何芳，等.礦山環(huán)境地質(zhì)問題綜合評價指標(biāo)體系[J].地質(zhì)通報，2003(10)：829-832.

[14] 樊新剛.寧夏限制開發(fā)生態(tài)區(qū)生態(tài)經(jīng)濟系統(tǒng)協(xié)調(diào)發(fā)展研究[D].銀川：寧夏大學(xué)，2015.

[15] 韓英強，吳曉平，王甲生.一種基于熵權(quán)和經(jīng)驗因子的模糊綜合評價方法[J].計算機與數(shù)字工程，2012，40(11)：105-107.

[16] 李曉雪，盧新衛(wèi).基于熵權(quán)的西安市大氣環(huán)境質(zhì)量模糊綜合評價[J].環(huán)境科學(xué)與技術(shù)，2012，35(4)：195-199.

[17] 楊曉艷，魯紅英.基于模糊綜合評判的城市環(huán)境空氣質(zhì)量評價[J].中國人口·資源與環(huán)境，2014(S2)：143-146.

[18] 王一任.綜合評價方法若干問題研究及其醫(yī)學(xué)應(yīng)用[D].長沙：中南大學(xué)，2012.

[19] Liu R，Qiang F U.Entropy weight coefficient model and its application in evaluation of groundwater vulnerability of the Sanjiang Plain[J].Journal of Northeast Agricultural University，2007，14(4)：368-373.

[20] 陸秋琴，潘長波，黃光球.基于未確知屬性層次法的礦山資源環(huán)境安全評價模型[J].安全與環(huán)境學(xué)報，2013(1)：234-240.

[21] 何芳，徐友寧，袁漢春.礦山環(huán)境地質(zhì)問題綜合評價客觀權(quán)值確定方法探討[J].中國地質(zhì)，2008，35(2)：337-343.

[22] 孫志文，屈彥.寧夏回族自治區(qū)資源環(huán)境地圖集[M].北京：中國地圖出版社，2006：30-31.

[23] 章文波，謝云，劉寶元.中國降雨侵蝕力空間變化特征[J].山地學(xué)報，2003，21(1)：33-40.

[24] 楊金中，秦緒文，聶洪峰，等.中國礦山遙感監(jiān)測[M].北京：測繪出版社，2014：111-112.

Environmental quality characteristics and spatial division of Ningxia mines

ZHANG Yongting1，2，CHEN Jianping1，LIU Zizeng2，YAN Jun2，HUANG Tao2，YAN Jin2，F(xiàn)AN Xingang3

(1.School of Earth Sciences and Resources，China University of Geosciences(Beijing)，Beijing 100083, China； 2.Ningxia Institute of Remote Sensing Surveying and Mapping，Yinchuan 750021，China； 3.Deparment of Social Service，Ningxia University，Yinchuan 750021，China)

Mines environment quality spatial distribution is important to government，which can provide the basis to repair and manage the mines environment of 22 counties in Ningxia province of China.Mines environment comprehensive evaluation index system was been built，include land use，geologic damage，environment pollution，water loss and soil erosion，and 16 indexs.The data was gotten by remote sensing data and mines environment statistics.Entropy weight fuzzy comprehensive method was used to evaluate the mines environmental quality，and then fuzzy clustering method by spss22.0 was used to divide the spatial partitions of mines environmental quality in Ningxia，and combined with field reconnaissance to determine result.The results show that the mines environmental quality in 22 counties is very differences.The northern part of Ningxia is generally better than the southern mountain area，but some of northern counties in Ningxia are worst.The mines environmental quality variation of 22 counties are relatively stable，and the main environmental problems are in the northern part of the old mining area of Huinong and Dawukou county，emerging mining area of Lingwu city and middle part of Ningxia i.e.Shapotou county.Then，Ningxia mines can be divided into four categories：“serious area，more serious area，general area，good area”.The distribution of mines environmental quality partitions is continuous from the province scale，but it concentrates from the counties scale.It should pay attention to the mining progress of the emerging mines，and continuously monitor on the mines environment，so as to early discovery and early management.The thresholds of categories by means of cluster analysis are discontinuity.It should increase the sample of cases，to find the threshold distribution of entropy weight fuzzy comprehensive evaluation value，so as to realize the continuous grading.

mine environmental quality；spatial division；entropy weight；fuzzy comprehensive evaluation；Ningxia

2016-07-14

全國礦產(chǎn)資源開發(fā)環(huán)境遙感監(jiān)測子項目“寧夏自治區(qū)礦產(chǎn)資源開發(fā)環(huán)境遙感監(jiān)測”資助(編號：1210115061301)

張永庭(1975-)，男，寧夏平羅縣人，博士研究生，教授級高級工程師，主要從事資源環(huán)境遙感研究，E-mail：fanxg@nxu.edu.cn。

簡介：陳建平(1959-)，男，福建閩候人，博士后，教授，博士生導(dǎo)師，主要從事礦產(chǎn)資源預(yù)測評價、遙感與地理信息系統(tǒng)應(yīng)用開發(fā)研究，E-mail：516948563@qq.com。

X321

A

1004-4051(2017)01-0049-06