(1.成都理工大學(xué) 國(guó)土資源部地學(xué)空間信息技術(shù)重點(diǎn)實(shí)驗(yàn)室，四川 成都 610059; 2.中國(guó)地質(zhì)調(diào)查局 成都地質(zhì)調(diào)查中心，四川 成都 610059)

城市生態(tài)環(huán)境問(wèn)題已然成為社會(huì)發(fā)展的限制性因素，而對(duì)于城市生態(tài)環(huán)境評(píng)價(jià)不僅可以幫助公眾了解自身生活環(huán)境的狀況，也能為城市規(guī)劃的設(shè)計(jì)者及生態(tài)環(huán)境監(jiān)測(cè)、治理人員提供幫助，因此如何準(zhǔn)確、有效地對(duì)城市生態(tài)環(huán)境進(jìn)行評(píng)價(jià)顯得尤為重要[1-2]。

目前，國(guó)內(nèi)外有許多針對(duì)城市生態(tài)環(huán)境評(píng)價(jià)的研究，并取得較大的進(jìn)步，研究策略也開(kāi)始從單一因素模型向多因素的模型評(píng)價(jià)模式轉(zhuǎn)變，并且逐漸從定性向定量過(guò)渡[2-5]。

雖然國(guó)內(nèi)已經(jīng)建立了許多城市環(huán)境評(píng)價(jià)體系，但是這些體系大都為某單一模型設(shè)計(jì)，其區(qū)域擴(kuò)展性較差，并且需要評(píng)價(jià)人員進(jìn)行大量繁雜的計(jì)算。因此，本文從多模型多指標(biāo)的角度出發(fā)，建立了一套通用的城市生態(tài)環(huán)境評(píng)價(jià)體系，并以GIS空間分析和可視化技術(shù)對(duì)評(píng)價(jià)模型的結(jié)果進(jìn)行顯示與分析，以此來(lái)解決體系擴(kuò)展性較差的問(wèn)題。而針對(duì)評(píng)價(jià)體系復(fù)雜多變，計(jì)算繁雜的問(wèn)題，本文利用ArcGIS Engine二次開(kāi)發(fā)出了一個(gè)評(píng)估系統(tǒng)，以提高人員的操作性和評(píng)價(jià)的準(zhǔn)確率。

1 城市生態(tài)環(huán)境評(píng)價(jià)體系的構(gòu)建

1.1 研究區(qū)及數(shù)據(jù)

本文選取天府新區(qū)成都直管區(qū)為研究區(qū)，研究區(qū)下轄有15個(gè)鄉(xiāng)鎮(zhèn)街道，常住人口54.19萬(wàn)人，氣候春夏日照足，秋冬云霧多，四季分明。該區(qū)地處成都平原南部邊緣地帶，所處大地構(gòu)造位置為新華夏系四川沉降帶成都斷陷的東南邊緣地帶，其地質(zhì)構(gòu)造主要表現(xiàn)為褶皺與斷裂。

遙感影像數(shù)據(jù)采用高分2號(hào)影像，過(guò)境時(shí)間段為2016年6月；在室內(nèi)遙感解譯之后，研究團(tuán)隊(duì)分別于2017年7～9月及2018年8月對(duì)研究區(qū)域進(jìn)行野外實(shí)地驗(yàn)證考察；并結(jié)合《成都市統(tǒng)計(jì)年鑒——2018》的數(shù)據(jù)整理出城市生態(tài)環(huán)境評(píng)價(jià)所需數(shù)據(jù)種類(lèi)，主要包括以下4種：① 行政區(qū)劃等基礎(chǔ)地理數(shù)據(jù)；② 地形地貌、道路、地質(zhì)災(zāi)害等城市生態(tài)環(huán)境評(píng)價(jià)因子指標(biāo)數(shù)據(jù)；③ 生態(tài)環(huán)境評(píng)價(jià)指標(biāo)量化分級(jí)數(shù)據(jù)；④ 評(píng)價(jià)指標(biāo)權(quán)重?cái)?shù)據(jù)。其中①、②為城市生態(tài)環(huán)境評(píng)價(jià)所需的矢量柵格數(shù)據(jù)，由于空間數(shù)據(jù)的復(fù)雜性，本文采用ArcGIS的Personal Geodatabase管理空間數(shù)據(jù)集；③為評(píng)價(jià)指標(biāo)量化分級(jí)提供數(shù)據(jù)支撐，有唯一賦值和區(qū)間賦值兩種。

1.2 評(píng)價(jià)指標(biāo)確定

由于城市生態(tài)環(huán)境影響指標(biāo)眾多，且各影響因素重要程度各不相同，因此要想建立一個(gè)包含完整影響因素的評(píng)價(jià)模型是不現(xiàn)實(shí)的。本文在兩次野外實(shí)際考察驗(yàn)證，以及充分考慮城市生態(tài)環(huán)境復(fù)雜性和多樣性基礎(chǔ)上，查閱相關(guān)文獻(xiàn)[6-9]和征詢(xún)專(zhuān)家意見(jiàn)，制定了適合于城市生態(tài)環(huán)境的評(píng)價(jià)的體系。

該評(píng)價(jià)體系一共分為3層，第1層為目標(biāo)層，第2層為要素層，第3層為指標(biāo)層。其中要素層又分為：自然要素層、社會(huì)要素層和生態(tài)環(huán)境層。自然要素層主要包括：植被覆蓋程度、地形地貌、水系、坡度、斷層、海拔、降水量、日照時(shí)間共8種自然環(huán)境指標(biāo)因子。社會(huì)要素層包括：公路交通、土地利用類(lèi)型、居民點(diǎn)分布、軌道交通、地區(qū)GDP共5種社會(huì)環(huán)境指標(biāo)因子。生態(tài)要素層包括：地質(zhì)災(zāi)害點(diǎn)、土壤重金屬污染、水體重金屬污染、噪聲污染、空氣污染共5種生態(tài)環(huán)境指標(biāo)因子，見(jiàn)表1。

表1 城市生態(tài)環(huán)境評(píng)價(jià)指標(biāo)體系Tab.1 Urban ecological environment assessment index system

通過(guò)分析、比較將目標(biāo)地區(qū)實(shí)際存在的影響環(huán)境因子作為評(píng)價(jià)指標(biāo)，并以此來(lái)建立適合目標(biāo)區(qū)域的評(píng)價(jià)體系。本文針對(duì)研究區(qū)實(shí)際情況，選?。旱匦蔚孛?N1)、河流(N2)、坡度(N3)、斷層(N4)、海拔(N5)，國(guó)道(S1)、鄉(xiāng)鎮(zhèn)道路(S2)、高速(S3)、地鐵(S4)、鐵路(S5)、居民地分布(S6)，土地利用類(lèi)型(E1)、地質(zhì)災(zāi)害點(diǎn)(E2)、土壤重金屬污染(E3)、水體重金屬污染(E4)共15個(gè)指標(biāo)作為評(píng)價(jià)因子，并對(duì)15個(gè)評(píng)價(jià)指標(biāo)按評(píng)價(jià)標(biāo)度賦值，見(jiàn)表2。

1.3 評(píng)價(jià)指標(biāo)權(quán)重確定

評(píng)價(jià)指標(biāo)的權(quán)重可以在多屬性評(píng)價(jià)中反映各個(gè)評(píng)價(jià)指標(biāo)對(duì)同一評(píng)價(jià)單元的相對(duì)重要程度，并直接影響評(píng)價(jià)結(jié)果的可信度。目前，權(quán)重確定的方法分為主觀(guān)賦權(quán)和客觀(guān)賦權(quán)兩類(lèi)。主觀(guān)賦權(quán)主要有：層次分析法、最小平方法、專(zhuān)家打分法等，該類(lèi)方法的優(yōu)點(diǎn)是考慮到人為主觀(guān)因素根據(jù)實(shí)際問(wèn)題進(jìn)行排序，但未考慮到指標(biāo)數(shù)據(jù)本身的依賴(lài)關(guān)系；客觀(guān)賦權(quán)主要有：變異系數(shù)法、熵值法、主成分分析法等，該類(lèi)方法的優(yōu)點(diǎn)是可以避免人為主觀(guān)因素，但是未考慮評(píng)價(jià)指標(biāo)之間的相對(duì)重要性[10-12]。

為了準(zhǔn)確地刻畫(huà)評(píng)價(jià)單元的重要程度，本文根據(jù)兩類(lèi)方法的優(yōu)缺點(diǎn)，擬采用主客觀(guān)綜合賦權(quán)法來(lái)最大程度地優(yōu)化權(quán)重取值。并通過(guò)查閱文獻(xiàn)和多次實(shí)驗(yàn)，最終采用層次分析法計(jì)算主觀(guān)權(quán)重，變異系數(shù)法計(jì)算客觀(guān)權(quán)重，再通過(guò)組合賦權(quán)法得到可信度較高的權(quán)重分配。

1.3.1基于層次分析法的權(quán)重賦值

層次分析法是一種將復(fù)雜目標(biāo)分解為多層次、多指標(biāo)的主流權(quán)重計(jì)算方法。它的核心是通過(guò)把目標(biāo)問(wèn)題分解為相關(guān)聯(lián)的元素，各關(guān)聯(lián)元素再通過(guò)屬性及關(guān)系形成若干層次，上一層元素作為準(zhǔn)則對(duì)下一層次有決定作用，這些層次分為3層：目標(biāo)層、中間層、決策層[13]。層次分析方法最關(guān)鍵部分是根據(jù)判斷矩陣計(jì)算中間層各個(gè)評(píng)價(jià)指標(biāo)的權(quán)重，從而在最終方案中衡量評(píng)價(jià)指標(biāo)對(duì)環(huán)境的影響程度，計(jì)算權(quán)重的主要步驟如下。

(1) 計(jì)算各個(gè)指標(biāo)因子的權(quán)重

(1)

(2) 計(jì)算判斷矩陣的最大特征值λmax及對(duì)判斷矩陣的一致性檢驗(yàn)。

表2 天府新區(qū)成都直管區(qū)指標(biāo)標(biāo)度等級(jí)Tab.2 Indicator scale of Chengdu direct management area of Tianfu new District

(2)

(3)

(4)

由于人們對(duì)客觀(guān)事物的認(rèn)知有限，在選取評(píng)價(jià)指標(biāo)的相對(duì)重要性會(huì)存在偏差，因此一般會(huì)借用式(3)式和(4)對(duì)判斷矩陣進(jìn)行一致性檢驗(yàn)。若C.R<0.1，則判斷矩陣滿(mǎn)足要求，否則需反復(fù)調(diào)整判斷矩陣直至滿(mǎn)足要求為止。計(jì)算得出的研究區(qū)評(píng)價(jià)指標(biāo)的主觀(guān)權(quán)重值見(jiàn)表3。

表3 評(píng)價(jià)指標(biāo)權(quán)重賦值Tab.3 Evaluation index weight assignment table

1.3.2基于變異系數(shù)法的權(quán)重賦值

變異系數(shù)法是通過(guò)衡量統(tǒng)計(jì)數(shù)據(jù)中特征值變異程度，來(lái)刻畫(huà)評(píng)價(jià)指標(biāo)的影響程度。在多指標(biāo)評(píng)價(jià)系統(tǒng)的指標(biāo)中，各觀(guān)測(cè)值的差異程度越大，表明此指標(biāo)重要性越高，也即此指標(biāo)權(quán)重越大[12,14-15]。計(jì)算步驟如下。

(1) 由于評(píng)價(jià)指標(biāo)的量綱不同，在進(jìn)行變異系數(shù)權(quán)重計(jì)算之前要先進(jìn)行無(wú)量綱化處理。

(5)

(6)

其中，Di為對(duì)i個(gè)指標(biāo)標(biāo)準(zhǔn)化之后的值，Xi為第i類(lèi)指標(biāo)的原始數(shù)據(jù)，max(Xi)為第i類(lèi)指標(biāo)數(shù)據(jù)的最大值，min(Xi)為第i類(lèi)指標(biāo)數(shù)據(jù)的最小值。式(6)針對(duì)正向指標(biāo)數(shù)據(jù)，式(7)針對(duì)負(fù)向指標(biāo)數(shù)據(jù)。

(2) 計(jì)算變異系數(shù)。

(7)

(3) 計(jì)算權(quán)重：

(8)

2 評(píng)價(jià)模型的建立與二次開(kāi)發(fā)系統(tǒng)

2.1 評(píng)價(jià)指標(biāo)量化分級(jí)

按構(gòu)建的評(píng)估體系標(biāo)準(zhǔn)對(duì)城市環(huán)境評(píng)價(jià)指標(biāo)進(jìn)行量化分級(jí)，其核心思想是對(duì)影像像元灰度值進(jìn)行二次賦值。量化分級(jí)數(shù)據(jù)的標(biāo)準(zhǔn)根據(jù)評(píng)價(jià)指標(biāo)數(shù)據(jù)類(lèi)型的不同而不同，一般分為兩種：唯一賦值和區(qū)間賦值，其中區(qū)間賦值是最常用的模式，使用者可自行指定起始、終止、新值的區(qū)間范圍 。對(duì)研究區(qū)15個(gè)指標(biāo)數(shù)據(jù)進(jìn)行量化分級(jí)，其中S1，N1，N4，E4評(píng)價(jià)指標(biāo)量化分級(jí)結(jié)果如圖1所示。

圖1 部分評(píng)價(jià)指標(biāo)量化分級(jí)Fig.1 Quantification and grading of some evaluation indicators

2.2 評(píng)價(jià)模型的建立

評(píng)價(jià)模型建立的基本思想是利用地圖代數(shù)理論將評(píng)價(jià)體系和量化分級(jí)數(shù)據(jù)組合。再通過(guò)層次分析法和變異系數(shù)法計(jì)算指標(biāo)權(quán)重之后，對(duì)單個(gè)評(píng)價(jià)單元引用線(xiàn)性組合賦權(quán)的模型計(jì)算指標(biāo)權(quán)重，并組合量化分級(jí)數(shù)據(jù)，模型如下：

(9)

其中，Ni為每個(gè)像元新的灰度值，Mi為第i類(lèi)指標(biāo)層次分析法計(jì)算權(quán)重，Ai為第i類(lèi)指標(biāo)變異系數(shù)法計(jì)算的權(quán)重，γ為喜好系數(shù)(0<<γ<<1)，PIXi為第i類(lèi)指標(biāo)的影像量化分級(jí)影像，為了后期分析展示，需要對(duì)數(shù)據(jù)進(jìn)行歸一化處理。

2.3 二次開(kāi)發(fā)的輔助系統(tǒng)

在實(shí)際過(guò)程中，層次分析法和變異系數(shù)法計(jì)算繁雜，對(duì)操作人員專(zhuān)業(yè)知識(shí)要求較高，為了便于后續(xù)評(píng)價(jià)結(jié)果的可視化與分析，本文另行設(shè)計(jì)了一個(gè)基于ArcGIS Engine的二次開(kāi)發(fā)系統(tǒng)。該評(píng)估系統(tǒng)屬于地理信息系統(tǒng)，基于地理信息系統(tǒng)基本理論，在計(jì)算機(jī)硬件與C#基礎(chǔ)應(yīng)用框架下，調(diào)用ArcGIS Engine中的GeoAnalyst類(lèi)庫(kù)及RasterMapAlegbraOp單元作為核心技術(shù)進(jìn)行二次開(kāi)發(fā)。由于系統(tǒng)開(kāi)發(fā)不屬于本文討論的主要內(nèi)容，因此僅在此作簡(jiǎn)要說(shuō)明。

3 評(píng)價(jià)模型的應(yīng)用

《成都市統(tǒng)計(jì)年鑒——2018》的數(shù)據(jù)顯示，天府新區(qū)成都直管區(qū)在規(guī)劃之前屬于待開(kāi)發(fā)區(qū)域，沒(méi)有大型污染企業(yè)，后期經(jīng)濟(jì)建設(shè)嚴(yán)格按照國(guó)家標(biāo)準(zhǔn)執(zhí)行，因此基本不存在由經(jīng)濟(jì)建設(shè)造成的嚴(yán)重生態(tài)污染；并且該區(qū)發(fā)展規(guī)劃也較為合理，即社會(huì)因素對(duì)城市生態(tài)環(huán)境影響相對(duì)較小，所以生態(tài)環(huán)境的主要影響因素為自然因素。為了選取合適的評(píng)價(jià)指標(biāo)以及對(duì)評(píng)價(jià)結(jié)果進(jìn)行驗(yàn)證，研究團(tuán)隊(duì)在2017年7～9月與2018年8月兩次對(duì)研究區(qū)進(jìn)行實(shí)地考察。考察發(fā)現(xiàn)：研究區(qū)東部沿龍泉山方向有較多連續(xù)磚瓦廠(chǎng)，開(kāi)采規(guī)模較大；龍泉山中發(fā)育許多零散分布的地質(zhì)災(zāi)害點(diǎn)，其中大部分已經(jīng)得到治理恢復(fù)；在西部主城區(qū)域內(nèi)部及城鎮(zhèn)周邊含有大量待建或正在建設(shè)的區(qū)域，地面大片裸露，部分水體污染也較為嚴(yán)重，生態(tài)環(huán)境有較大破壞；除此之外的大部分研究區(qū)域均為農(nóng)村未開(kāi)發(fā)區(qū)域，除有少量生活、生產(chǎn)垃圾，整體生態(tài)環(huán)境較好。

因?yàn)楦鱾€(gè)評(píng)價(jià)指標(biāo)具有不同的量綱，為了便于后期的統(tǒng)計(jì)分析，對(duì)各個(gè)評(píng)價(jià)指標(biāo)進(jìn)行無(wú)量綱化處理。在利用計(jì)算主觀(guān)權(quán)重賦值和客觀(guān)權(quán)重之后，為體現(xiàn)主客觀(guān)權(quán)重對(duì)評(píng)價(jià)結(jié)果的影響程度，根據(jù)多次實(shí)驗(yàn)調(diào)整優(yōu)化，分別取研究區(qū)的喜好系數(shù)γ為0.25,0.50,0.75進(jìn)行計(jì)算比較,其第i類(lèi)指標(biāo)對(duì)應(yīng)的綜合權(quán)重見(jiàn)表3。對(duì)每一個(gè)評(píng)價(jià)指標(biāo)進(jìn)行量化分級(jí)處理后利用設(shè)計(jì)的二次開(kāi)發(fā)系統(tǒng)中的評(píng)價(jià)模型編輯模塊將各個(gè)指標(biāo)因子分級(jí)結(jié)果與指標(biāo)因子權(quán)重進(jìn)行地圖代數(shù)處理，如圖2所示。然后對(duì)取不同偏好系數(shù)情況評(píng)價(jià)結(jié)果分析，并進(jìn)行彩色渲染展示，如圖3所示。為充分體現(xiàn)研究區(qū)域內(nèi)生態(tài)環(huán)境等級(jí)分布及波動(dòng)情況，在充分考慮3種喜好系數(shù)以及兩次野外調(diào)查結(jié)果之后，最終將本研究區(qū)的生態(tài)環(huán)境區(qū)間設(shè)定為：差、較差、一般、好、較好5個(gè)等級(jí)，對(duì)應(yīng)圖3的圖例等級(jí)1～5。

圖2 二次開(kāi)發(fā)系統(tǒng)的模型計(jì)算編輯界面Fig.2 Model calculation editing interface of the secondary development system

圖3 天府新區(qū)成都直管區(qū)城市生態(tài)等級(jí)空間分布Fig.3 Spatial distribution of urban ecological hierarchy in the Chengdu direct management area of Tianfu new District

將3種不同喜好系數(shù)的評(píng)價(jià)結(jié)果與兩次野外考察結(jié)論進(jìn)行對(duì)比后發(fā)現(xiàn)：γ=0.25時(shí),對(duì)研究區(qū)域的生態(tài)環(huán)境較差區(qū)域進(jìn)行準(zhǔn)確的劃分，評(píng)價(jià)結(jié)果與兩次實(shí)際考察結(jié)論最為接近；γ=0. 50時(shí),對(duì)西部城市區(qū)域的評(píng)價(jià)效果不明顯；γ=0.75時(shí),對(duì)西部城市區(qū)域的評(píng)價(jià)效果不佳，且對(duì)中部未開(kāi)發(fā)區(qū)域評(píng)價(jià)也存在過(guò)度劃分的情況。因此，γ=0.25時(shí)整體效果最優(yōu)，γ=0.50時(shí)次之，γ=0.75時(shí)效果最差。

將圖3(a)與兩次野外考察驗(yàn)證結(jié)果對(duì)比分析，并結(jié)合《成都市統(tǒng)計(jì)年鑒——2018》可以得出天府新區(qū)成都直管區(qū)生態(tài)等級(jí)較低區(qū)域主要分布在西北部和東部邊界區(qū)域，中間大部分區(qū)域評(píng)價(jià)等級(jí)較高。西部區(qū)域主要沿城市主干道分布，集中于人類(lèi)活動(dòng)較為強(qiáng)烈的低坡度低高程地區(qū)，植被覆蓋程度較低，生態(tài)環(huán)境等級(jí)以較差及差為主，其對(duì)應(yīng)面積占比為28.93%，人類(lèi)活動(dòng)是導(dǎo)致生態(tài)環(huán)境變差的重要因素。東部邊界呈條狀分布主要是由于該區(qū)域位于東北-西南走向的龍泉山斷裂帶上，該區(qū)域海拔高、坡度大以及包含一些地質(zhì)災(zāi)害點(diǎn)從而影響評(píng)價(jià)等級(jí)。在此之間的中部及南部區(qū)域大都為未開(kāi)發(fā)地區(qū)，植被覆蓋程度高，區(qū)域內(nèi)分布零星居住地及農(nóng)莊，但人類(lèi)活動(dòng)強(qiáng)度不大對(duì)整體評(píng)價(jià)等級(jí)影響較小，生態(tài)環(huán)境評(píng)價(jià)等級(jí)為一般、好、較好的面積占總研究區(qū)的比例為71.07%，整體生態(tài)環(huán)境良好。

4 結(jié) 語(yǔ)

基于主客觀(guān)賦權(quán)的城市生態(tài)環(huán)境評(píng)價(jià)方法，可以準(zhǔn)確地建立城市生態(tài)環(huán)境中各層次、其指標(biāo)與城市生態(tài)質(zhì)量等級(jí)之間的映射關(guān)系，而且對(duì)于城市環(huán)境評(píng)價(jià)精度較高，評(píng)價(jià)體系及模型簡(jiǎn)潔實(shí)用。在實(shí)際運(yùn)用中結(jié)合ArcGIS Engine對(duì)評(píng)價(jià)系統(tǒng)進(jìn)行二次開(kāi)發(fā)，進(jìn)一步簡(jiǎn)化了流程，操作人員無(wú)需進(jìn)行繁復(fù)的運(yùn)算和擁有很強(qiáng)的專(zhuān)業(yè)知識(shí)即可完成，為大范圍的城市生態(tài)環(huán)境評(píng)價(jià)提供了新的思路。