李長(zhǎng)泓
(河南理工大學(xué)礦山空間信息技術(shù)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室,河南焦作454003)
長(zhǎng)期以來(lái),我國(guó)的經(jīng)濟(jì)增長(zhǎng)是一種以消耗大量礦產(chǎn)資源為代價(jià)的粗放式增長(zhǎng)。部分礦區(qū),只顧眼前經(jīng)濟(jì)效益,對(duì)開采引起的環(huán)境問(wèn)題重視不足,廢水廢渣任意排放堆積,污染地下水、土壤和生態(tài)環(huán)境,嚴(yán)重影響了人們的生活質(zhì)量,制約了區(qū)域經(jīng)濟(jì)健康發(fā)展。因此,加強(qiáng)礦區(qū)環(huán)境質(zhì)量檢測(cè)與評(píng)價(jià)已經(jīng)成為保護(hù)礦區(qū)環(huán)境的重要內(nèi)容。目前,正在開展的地理國(guó)情普查及相應(yīng)的專題檢測(cè)工作可為礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)提供有效的數(shù)據(jù)支持。本文首先論述了地理國(guó)情普查對(duì)礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)的支持作用,然后將地理國(guó)情數(shù)據(jù)與礦區(qū)專題數(shù)據(jù)進(jìn)行融合與集成,運(yùn)用模糊數(shù)學(xué)方法建立礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)模型,對(duì)礦區(qū)環(huán)境質(zhì)量做出綜合評(píng)定,為礦區(qū)環(huán)境保護(hù)和治理提供科學(xué)依據(jù)。
環(huán)境質(zhì)量評(píng)價(jià)[1]是按照一定的標(biāo)準(zhǔn)和方法對(duì)環(huán)境質(zhì)量給予定性和定量的說(shuō)明與描述。其中,礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)屬于環(huán)境質(zhì)量評(píng)價(jià)按照區(qū)域分類中的一個(gè)組成部分。礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)可以定義為通過(guò)研究礦區(qū)質(zhì)量環(huán)境質(zhì)量變化情況,按照環(huán)境質(zhì)量評(píng)價(jià)方法,對(duì)礦區(qū)環(huán)境質(zhì)量的優(yōu)劣做出定性或者定量的評(píng)定。
礦區(qū)活動(dòng)勢(shì)必會(huì)對(duì)礦區(qū)周圍的自然和地理環(huán)境產(chǎn)生擾動(dòng),影響整個(gè)礦區(qū)的環(huán)境質(zhì)量,而礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)正是對(duì)礦區(qū)環(huán)境質(zhì)量變化產(chǎn)生的影響,做出符合當(dāng)前實(shí)際情況的評(píng)估,從而能夠讓管理和決策有據(jù)可依,對(duì)保護(hù)和改善礦區(qū)環(huán)境質(zhì)量將更加有利有效。
針對(duì)礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)前人已經(jīng)做了很多工作,主要表現(xiàn)為兩個(gè)方面:一是在數(shù)學(xué)模型模擬下,用簡(jiǎn)單疊加影響因子的方式進(jìn)行環(huán)境質(zhì)量評(píng)價(jià);二是在局部區(qū)域內(nèi),用各自實(shí)測(cè)數(shù)據(jù)為依據(jù)進(jìn)行環(huán)境質(zhì)量評(píng)價(jià)。顯然,這兩種模式已經(jīng)不能滿足現(xiàn)在和今后礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)的需要。當(dāng)前在全國(guó)范圍內(nèi)進(jìn)展的地理國(guó)情普查工作,可為礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)的研究者提供統(tǒng)一有效的數(shù)據(jù),使得研究成果更具有可信度和說(shuō)服力。
礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)數(shù)據(jù)主要分為地理國(guó)情普查數(shù)據(jù)的獲取及相應(yīng)專題數(shù)據(jù)的收集。地理國(guó)情普查是為了獲取國(guó)土自然地理信息、人文地理信息等重要國(guó)情信息數(shù)據(jù),堅(jiān)持“應(yīng)采盡采”的原則[2]進(jìn)行的全國(guó)性、多方位的數(shù)據(jù)收集和調(diào)查研究。地理國(guó)情普查數(shù)據(jù)獲取主要包括地理國(guó)情數(shù)據(jù)的來(lái)源和地理國(guó)情數(shù)據(jù)的處理兩個(gè)方面。
2.1.1 地理國(guó)情數(shù)據(jù)的來(lái)源
地理國(guó)情普查數(shù)據(jù)主要包括內(nèi)業(yè)數(shù)據(jù)和外業(yè)數(shù)據(jù)。內(nèi)業(yè)數(shù)據(jù)主要來(lái)源于四部分[3]:(1)正射影像。以景為單位的全色及多波段衛(wèi)星影像和標(biāo)準(zhǔn)分幅DOM數(shù)據(jù);(2)區(qū)域地圖。如1∶10000和1∶50000DLG數(shù)據(jù);(3)專題數(shù)據(jù)資料。各職能部門提供的相關(guān)領(lǐng)域數(shù)據(jù)。如國(guó)土資源部門可以提供礦區(qū)的地理信息數(shù)據(jù);(4)其它途徑或者網(wǎng)絡(luò)公開的與之相關(guān)的數(shù)據(jù)信息。外業(yè)數(shù)據(jù)包括對(duì)道路、水系、構(gòu)筑物和地理單元及界線的實(shí)地量測(cè)。
2.1.2 地理國(guó)情數(shù)據(jù)的處理
對(duì)已有的基礎(chǔ)地理信息數(shù)據(jù)進(jìn)行提取、分析、整理和重構(gòu)。通過(guò)坐標(biāo)變換、圖幅變形校正、空間數(shù)據(jù)的壓縮處理、數(shù)據(jù)類型的轉(zhuǎn)換、數(shù)據(jù)內(nèi)插、圖幅匹配等過(guò)程,依據(jù)《地理國(guó)情普查內(nèi)容與指標(biāo)(試行稿)》的要求,提取并形成地理國(guó)情普查基礎(chǔ)數(shù)據(jù)。礦區(qū)專題數(shù)據(jù)資料對(duì)以形成的地理國(guó)情普查基礎(chǔ)數(shù)據(jù)進(jìn)行補(bǔ)充,形成地理國(guó)情普查參考數(shù)據(jù)。通過(guò)外業(yè)數(shù)據(jù)的核查,對(duì)地理國(guó)情普查參考數(shù)據(jù)進(jìn)行確認(rèn)和補(bǔ)調(diào)。補(bǔ)調(diào)是通過(guò)補(bǔ)測(cè)的形式進(jìn)行的,最終確定所獲得地理國(guó)情數(shù)據(jù)與真實(shí)地理信息準(zhǔn)確無(wú)誤。
用于礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)的地理國(guó)情普查數(shù)據(jù)和礦區(qū)專題數(shù)據(jù)都包含著礦區(qū)大量的地理信息,可以把兩種數(shù)據(jù)進(jìn)行融合,構(gòu)建礦區(qū)地理環(huán)境數(shù)據(jù)庫(kù),從而更加全面清楚地對(duì)整個(gè)礦區(qū)環(huán)境與生態(tài)要素進(jìn)行了解和掌握。要想實(shí)現(xiàn)礦區(qū)地理環(huán)境數(shù)據(jù)庫(kù)的構(gòu)建,就要把兩者的數(shù)據(jù)融合,這就必須對(duì)兩種數(shù)據(jù)進(jìn)行處理。首先,對(duì)于紙質(zhì)數(shù)據(jù)資料,用掃描儀或者數(shù)字化軟件將地理國(guó)情普查數(shù)據(jù)和礦區(qū)專題數(shù)據(jù)數(shù)字化,進(jìn)而矢量化,建立幾何數(shù)據(jù)庫(kù);然后,對(duì)兩幅數(shù)據(jù)進(jìn)行幾何校正、配準(zhǔn)、統(tǒng)一坐標(biāo)系;接著,將礦區(qū)專題數(shù)據(jù)融合到地理國(guó)情普查數(shù)據(jù)中,形成新的幾何數(shù)據(jù)庫(kù);最后,對(duì)融合數(shù)據(jù)進(jìn)行屬性編輯,添加各個(gè)地理實(shí)體的屬性,建立屬性數(shù)據(jù)庫(kù)。
不同數(shù)據(jù)有著其各自的比例尺,融合的關(guān)鍵在于對(duì)數(shù)據(jù)地圖的配準(zhǔn)[4]。宜選擇采取適當(dāng)個(gè)數(shù)的高精度控制點(diǎn)作為配準(zhǔn)控制點(diǎn)和雙線性內(nèi)插法進(jìn)行地理配準(zhǔn)。
除了建立礦區(qū)地理環(huán)境數(shù)據(jù)庫(kù),還要建立環(huán)境質(zhì)量評(píng)價(jià)因子數(shù)據(jù)庫(kù)?;贏rcGIS軟件平臺(tái),環(huán)境質(zhì)量評(píng)價(jià)因子數(shù)據(jù)庫(kù)可以有三部分構(gòu)成,即環(huán)境質(zhì)量評(píng)價(jià)因子實(shí)測(cè)數(shù)據(jù)庫(kù)、環(huán)境質(zhì)量評(píng)價(jià)因子評(píng)價(jià)標(biāo)準(zhǔn)數(shù)據(jù)庫(kù)和數(shù)學(xué)模型分析預(yù)測(cè)數(shù)據(jù)庫(kù)。借助GIS的強(qiáng)大分析和模塊擴(kuò)展功能,將礦區(qū)地理環(huán)境數(shù)據(jù)庫(kù)和環(huán)境質(zhì)量評(píng)價(jià)因子數(shù)據(jù)庫(kù)進(jìn)行集成,從而實(shí)現(xiàn)查詢、可視、預(yù)測(cè)、決策等功能。
環(huán)境質(zhì)量評(píng)價(jià)的過(guò)程中,評(píng)價(jià)的主體、評(píng)價(jià)方法和需要評(píng)價(jià)的對(duì)象都存在不確定因素,這就使得評(píng)價(jià)的結(jié)果也具有不確定性,而模糊數(shù)學(xué)為確定與不確定之間的溝通尋找了一條新的途徑。
礦山環(huán)境質(zhì)量的好壞不是某個(gè)單因素因子可以決定的,而是由多因素因子共同決定的,同時(shí)許多因子本身也具有模糊性,但是模糊綜合評(píng)價(jià)法在解決模糊性的定性和定量方面表現(xiàn)出特別的優(yōu)越性,也使得其結(jié)果更加客觀可信。因此,將模糊數(shù)學(xué)引入到礦山環(huán)境質(zhì)量評(píng)價(jià),能夠有效地處理評(píng)價(jià)因子對(duì)評(píng)價(jià)結(jié)果的不確定性問(wèn)題。
模糊綜合評(píng)價(jià)法數(shù)學(xué)模型[5-6]的建立步驟為:
設(shè)有環(huán)境質(zhì)量評(píng)價(jià)因子集
式中Ui是進(jìn)行評(píng)價(jià)體系的評(píng)價(jià)因子。
環(huán)境質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)集
式中:Vm表示評(píng)價(jià)因子根據(jù)評(píng)價(jià)對(duì)象的環(huán)境類質(zhì)量標(biāo)準(zhǔn)劃分的等級(jí);m根據(jù)劃分等級(jí)級(jí)別多少來(lái)確定。
在U、V確定以后,建立環(huán)境質(zhì)量評(píng)價(jià)因子和環(huán)境質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)之間的邊界模糊關(guān)系矩陣R,可表示為:
式中rnm表示第n個(gè)評(píng)價(jià)因子被劃分為等級(jí)Vm的可能性,即評(píng)級(jí)因子Un對(duì)評(píng)價(jià)等級(jí)Vm的隸屬度。
由于評(píng)價(jià)因子對(duì)評(píng)價(jià)主體的環(huán)境質(zhì)量影響存在正向作用和負(fù)向作用,因此評(píng)價(jià)因子對(duì)評(píng)價(jià)主體的環(huán)境質(zhì)量影響程度會(huì)偏向某一個(gè)評(píng)價(jià)等級(jí)。為了說(shuō)明這種偏向性,引入偏離度,運(yùn)用偏離度來(lái)描述隸屬度。偏離度是評(píng)價(jià)因子實(shí)測(cè)值對(duì)評(píng)價(jià)等級(jí)標(biāo)準(zhǔn)值的偏離程度。隸屬度是評(píng)價(jià)因子實(shí)測(cè)值對(duì)評(píng)價(jià)等級(jí)標(biāo)準(zhǔn)值的靠攏程度。偏離度與隸屬度恰好是兩個(gè)相反的過(guò),即位于相鄰等級(jí)標(biāo)準(zhǔn)值之間的評(píng)價(jià)因子,對(duì)前一等級(jí)標(biāo)準(zhǔn)的偏離度與對(duì)后一等級(jí)標(biāo)準(zhǔn)的隸屬度相同。因此,隸屬度可以依據(jù)隸屬度函數(shù)或者與偏離度的關(guān)系求得,從而得到隸屬度矩陣R。本文采用隸屬度與偏離度的關(guān)系確定各評(píng)價(jià)因子的隸屬函數(shù)。偏離度可以表示為:
式中:Ui表示評(píng)價(jià)因子實(shí)測(cè)值;Vt、Vt+1表示第t、t+1個(gè)評(píng)價(jià)等級(jí)標(biāo)準(zhǔn)值。
隸屬度的確定:
①Ui在V1的左側(cè)時(shí),Ui對(duì)V1的隸屬度為 1,對(duì){V2,V3,…,Vj},(j=2,3,…,m)的隸屬度為 0;Ui在Vj的右側(cè)時(shí),Ui對(duì)Vj的隸屬度為 1,對(duì){V1,V2,V3,…,Vj-1},(j=1,2,3,…,m)的隸屬度為 0。
②Ui在V1和Vj的之間時(shí),Ui∈(vt,vt+1)∈(v1,vj),Ui對(duì)Vt的偏離度為pit=n;j=1,2,3,…,m),Ui對(duì)Vt+1的偏離度為
因此,Ui∈(vt,vt+1)∈(v1,vj)時(shí),Ui對(duì)Vt的隸屬度為rit=Pit+1,Ui對(duì)Vt+1的隸屬度為rit+1=Pit。
根據(jù)模糊數(shù)學(xué)原理,對(duì)環(huán)境評(píng)判目標(biāo)集A,有以下模糊變換:
式中“·”表示模糊矩陣合成算子。
式(5)即為模糊綜合評(píng)價(jià)模型。依照隸屬度最大選擇原則,環(huán)境評(píng)判目標(biāo)集A中最大的值A(chǔ)max最為環(huán)境評(píng)判目標(biāo)的評(píng)判標(biāo)準(zhǔn)。
本文構(gòu)建了礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)模式。首先,選定評(píng)價(jià)礦區(qū),確定礦區(qū)主要污染源與污染物;其次,篩選出環(huán)境監(jiān)測(cè)因子,收集監(jiān)測(cè)數(shù)據(jù);最后,根據(jù)環(huán)境質(zhì)量評(píng)價(jià)方法,綜合評(píng)定礦區(qū)環(huán)境質(zhì)量的程度。
根據(jù)當(dāng)前的礦山開采現(xiàn)狀和環(huán)境質(zhì)量評(píng)價(jià)指標(biāo)[7],應(yīng)用等標(biāo)負(fù)荷法評(píng)判,可以確定主要影響礦區(qū)的環(huán)境質(zhì)量評(píng)價(jià)因素主要有水環(huán)境質(zhì)量、大氣質(zhì)量、礦井開發(fā)利用和開發(fā)利用后土地狀況。按照層次分析法原理和地理國(guó)情分類體系[8-9],構(gòu)建礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)指標(biāo)結(jié)構(gòu)[10-12](見(jiàn)圖1)。層次分析法是一種將總目標(biāo)問(wèn)題逐級(jí)分解成子目標(biāo)來(lái)研究,具有系統(tǒng)性的權(quán)重決策方法。它能有效解決無(wú)法準(zhǔn)確計(jì)量的模糊目標(biāo)定性或定量問(wèn)題。
在不同的層級(jí)上,不同的評(píng)價(jià)因子對(duì)上一層級(jí)的貢獻(xiàn)能力不一樣。為了確定各個(gè)評(píng)價(jià)因子在該層級(jí)上對(duì)上一級(jí)的貢獻(xiàn)能力大小,可以引入德爾菲法,以各因子所占權(quán)重大小作為衡量指標(biāo)。在圖1中,可以建立各層級(jí)的權(quán)系數(shù)矩陣,即亞目標(biāo)層(B)對(duì)目標(biāo)層(A)的權(quán)重 u=(0.24,0.24,0.39,0.13),指標(biāo)層(C)對(duì)亞目標(biāo)層(B)的權(quán)重 b1=(0.65,0.35)、b2=(0.49,0.31,0.20)、b3=(0.55,0.45),b4=(0.40,0.27,0.33)。
根據(jù)礦區(qū)的特點(diǎn)和礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)的目的,本文采用國(guó)內(nèi)外目前正在使用的環(huán)境質(zhì)量綜合評(píng)價(jià)方法:一種是兩步法,即首先進(jìn)行單要素評(píng)價(jià),然后歸納評(píng)價(jià)結(jié)果,得出礦區(qū)環(huán)境質(zhì)量綜合評(píng)價(jià);另一種是一步法,即根據(jù)評(píng)價(jià)目的,篩選出某些最能反映礦區(qū)環(huán)境質(zhì)量的環(huán)境要素和具有代表性的參數(shù),直接推出礦區(qū)環(huán)境質(zhì)量綜合評(píng)價(jià)值[13]。本文基于模糊數(shù)學(xué)的方法,構(gòu)建模糊綜合評(píng)價(jià)模型,對(duì)礦山環(huán)境質(zhì)量進(jìn)行評(píng)價(jià),符合兩步法(見(jiàn)圖2)的評(píng)價(jià)規(guī)則。
圖1 礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)指標(biāo)結(jié)構(gòu)
圖2 兩步法評(píng)價(jià)程序
以某礦區(qū)為例,結(jié)合相關(guān)評(píng)價(jià)因子以及礦區(qū)地理國(guó)情普查數(shù)據(jù)和礦區(qū)專題數(shù)據(jù)為數(shù)據(jù)源(見(jiàn)表1、表2),評(píng)價(jià)該礦區(qū)環(huán)境質(zhì)量的級(jí)別。
表1 某礦區(qū)2012年評(píng)價(jià)指標(biāo)數(shù)據(jù)(1)
表2 某礦區(qū)2012年評(píng)價(jià)指標(biāo)數(shù)據(jù)(2)
結(jié)合環(huán)境質(zhì)量評(píng)價(jià)原則和礦區(qū)管理實(shí)際需要,將礦區(qū)環(huán)境質(zhì)量劃分為四個(gè)等級(jí),即無(wú)污染、輕污染、中污染和重污染。與之相應(yīng)的各因子評(píng)價(jià)標(biāo)準(zhǔn)集V為:
根據(jù)實(shí)測(cè)數(shù)據(jù),計(jì)算各個(gè)評(píng)價(jià)因子的隸屬度矩陣:
水質(zhì)量(B1)的隸屬矩陣為R1,有
大氣質(zhì)量(B2)的隸屬矩陣為R2,有
開發(fā)利用(B3)的隸屬矩陣為R3,有
土地狀況(B4)的隸屬矩陣為R4,有
進(jìn)一步由式(5)計(jì)算出亞目標(biāo)層級(jí)對(duì)目標(biāo)層級(jí)的綜合權(quán)重系數(shù),即
其中矩陣A中各值代表各環(huán)境質(zhì)量評(píng)價(jià)因子對(duì)環(huán)境質(zhì)量評(píng)價(jià)等級(jí)的隸屬度,也就是各環(huán)境影響因子對(duì)環(huán)境質(zhì)量貢獻(xiàn)率的表征。取矩陣A中最大值A(chǔ)max=0.39作為相應(yīng)的評(píng)價(jià)結(jié)果,該數(shù)據(jù)顯示出該礦區(qū)的環(huán)境質(zhì)量處于輕污染等級(jí)級(jí)別。
地理國(guó)情普查是一項(xiàng)重要的國(guó)情監(jiān)測(cè)任務(wù),可以獲取到大量的地理信息數(shù)據(jù),這為礦山環(huán)境質(zhì)量的監(jiān)測(cè)和評(píng)價(jià)提供了準(zhǔn)確的空間屬性數(shù)據(jù),為建立礦山空間屬性數(shù)據(jù)倉(cāng)庫(kù)做了良好的準(zhǔn)備。結(jié)合模糊數(shù)學(xué)的方法,構(gòu)建環(huán)境質(zhì)量評(píng)價(jià)的模糊綜合評(píng)價(jià)模型,基于GIS平臺(tái),對(duì)地理國(guó)情數(shù)據(jù)和礦山專題數(shù)據(jù)的融合與集成,使得礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)智能化,為智慧礦山的建立提供有效的參考。模糊綜合評(píng)價(jià)法在礦山環(huán)境質(zhì)量評(píng)價(jià)中,通過(guò)隸屬度將模糊的定性環(huán)境質(zhì)量因子轉(zhuǎn)向定量化因子,實(shí)現(xiàn)了對(duì)礦區(qū)環(huán)境質(zhì)量的評(píng)價(jià)更加準(zhǔn)確、客觀和可信。
[1] 張從.環(huán)境評(píng)價(jià)教程[M].北京:中國(guó)環(huán)境科學(xué)出版社,2002.
[2] 虞欣,楊伯鋼,晁春浩,等.北京市地理國(guó)情普查試點(diǎn)的實(shí)踐[J].測(cè)繪通報(bào),2014(6):102-104.
[3] 董冬,龔偉.淺談地理國(guó)情普查基本要素內(nèi)容[J].測(cè)繪與空間地理信息,2013,36(8):199-201,206.
[4] 徐靜海,韓燁.基于多比例尺技術(shù)的南京市地震應(yīng)急數(shù)據(jù)庫(kù)建設(shè)[J].世界地震工程,2011,27(2):51-56.
[5] 高明美,孫濤,張坤.基于超標(biāo)倍數(shù)賦權(quán)法的濟(jì)南市大氣質(zhì)量模糊動(dòng)態(tài)評(píng)價(jià)[J].干旱區(qū)資源與環(huán)境,2014,28(9):150-154.
[6] 張敬花.基于模糊數(shù)學(xué)的華亭礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)研究[J].中央民族大學(xué)學(xué)報(bào):自然科學(xué)版,2011,20(3):41-43.
[7] 李政大,袁曉玲,楊萬(wàn)平.環(huán)境質(zhì)量評(píng)價(jià)研究現(xiàn)狀、困惑和展望[J].資源科學(xué),2014,36(1):175-181.
[8] 冷海芹.淺析地理國(guó)情普查內(nèi)容分類體系[J].測(cè)繪與空間地理信息,2014,37(6):118-121.
[9] 宋尚萍,陳世培,李井春,等.地理國(guó)情普查內(nèi)容與指標(biāo)體系構(gòu)建方法研究[J].測(cè)繪與空間地理信息,2014,37(6):60-62.
[10] 李堂軍,孫承愛(ài).礦區(qū)開發(fā)對(duì)環(huán)境的影響與評(píng)價(jià)方法[J].煤礦環(huán)境保護(hù),2000,14(4):56-59.
[11] 劉景凡,石忠旭.礦區(qū)環(huán)境質(zhì)量模糊綜合評(píng)價(jià)級(jí)別法及其應(yīng)用[J].青島建筑工程學(xué)院學(xué)報(bào),2002,23(2):31-34.
[12] 李堂軍,孫承愛(ài),魏志恒,等.兗州礦區(qū)環(huán)境質(zhì)量評(píng)價(jià)與適應(yīng)性對(duì)策[J].礦業(yè)研究與開發(fā),2002,22(3):20-22.
[13] 張從.環(huán)境評(píng)價(jià)教程[M].北京:中國(guó)環(huán)境科學(xué)出版社,2002.