陳理庭，蔡海生,2，張 婷，張學玲，曾 珩，陳 藝

(1.江西農業(yè)大學江西省鄱陽湖流域農業(yè)資源與生態(tài)重點實驗室/江西農業(yè)大學富硒農業(yè)產(chǎn)業(yè)發(fā)展研究中心，南昌 330045；2.江西旅游商貿職業(yè)學院，南昌 330100；3.江西省土地開發(fā)整理中心，南昌 330002)

【研究意義】國土空間作為人類生產(chǎn)生活的重要載體，理應受到嚴格保護。隨著社會經(jīng)濟的飛速發(fā)展，人類活動日漸頻繁、城鎮(zhèn)建設空間擴張無序等問題對國土空間造成了巨大壓力，生態(tài)問題成為制約社會發(fā)展的重要因素之一。與此同時，隨著國土空間規(guī)劃工作的全面鋪開，國土空間生態(tài)修復提升至國家戰(zhàn)略工程，生態(tài)文明思想深入人心。進行生態(tài)脆弱性評價可以有效掌握國土空間的本底特征[1]，為國土空間開發(fā)利用以及生態(tài)修復提供重要依據(jù)，具有重要現(xiàn)實意義。【前人研究進展】當今學者大都認為：生態(tài)環(huán)境脆弱性是指在自然要素和人為因素的共同影響下，生態(tài)系統(tǒng)對于干擾所表現(xiàn)的敏感性和恢復力，只要干擾程度超出了生態(tài)系統(tǒng)的適應能力就可能會出現(xiàn)相應的脆弱表征[2]。近年來，隨著對生態(tài)脆弱性研究的不斷深入，其評價體系正朝著多樣化和科學化邁進。目前，生態(tài)脆弱性評價方法多采用多因子綜合的定量研究方法，如綜合指數(shù)法[3]、主成分分析法[4-5]等。指標模型多為PSR(壓力—狀態(tài)—響應)模型[6-7]、SRP(Sensitivity-Recovery-Pressure)模型[8-10]以及基于“暴露—敏感—適應能力”的VSD模型[11-12]等,其研究尺度普遍選擇的是流域、省域以及市域。但這些現(xiàn)有研究多應用一組指標權重進行線性組合的方式如AHP-PCA熵權法[13]等進行動態(tài)評價，缺乏對不同決策風險的考慮。有序加權平均(OWA)法是一種通過調整風險系數(shù)以獲得不同評價結果的多準則決策法，可以有效解決多屬性決策的問題，其在城市減災[14]、洪水災害風險評價[15]、建設用地開發(fā)適宜性[16]、土地生態(tài)安全評價[17]、生態(tài)安全格局[18]、土壤適宜性評價[19]、生態(tài)風險評價[20]等領域得到廣泛應用，而在生態(tài)環(huán)境脆弱性研究中涉及較少?！颈狙芯壳腥朦c】基于SRP概念模型，從生態(tài)敏感性、生態(tài)恢復力、生態(tài)壓力度3個層面選取11個具體指標，運用綜合指數(shù)法對萬年縣2005—2018年生態(tài)脆弱性進行綜合分析，揭示其時空演變特征，同時結合空間自相關分析，探討其在空間上的集聚格局，以期為國土空間生態(tài)修復、分區(qū)治理提供參考。此外，引入OWA多準則決策方法對萬年縣2018年生態(tài)脆弱性進行情景模擬，分析其在不同決策風險下的評價結果，輔助決策者進行更為科學全面的判斷，以期為區(qū)域長遠發(fā)展提供依據(jù)?！緮M解決的關鍵問題】此研究為縣域生態(tài)脆弱性評價補充了多重視角，為當前生態(tài)脆弱性研究多情景模擬分析的不足提供了可行方案，為萬年縣進一步踐行生態(tài)文明建設、生態(tài)保護修復等提供了科學依據(jù)，對促進其社會—經(jīng)濟—生態(tài)協(xié)調發(fā)展具有重要意義。

1 材料與方法

1.1 研究區(qū)概況

萬年縣(116°46′41″～117°15′16″E、28°30′00″～28°54′08 ″N)處于江西省東南部、上饒市中部(圖1)。年均降水量約1808 mm，年均氣溫約17.4 ℃，縣內植被良好，雨量充足。作為稻米之鄉(xiāng)，全縣土地利用以耕地為主，居民點和工礦用地次之。境內以丘陵為主，地勢從東南向西北降低，中部丘陵連綿，西部湖泊眾多，水系發(fā)達。近幾年，萬年縣建設開發(fā)，人類活動頻率上升，對其生態(tài)環(huán)境造成了一定壓力。萬年縣域內礦產(chǎn)資源豐富，但礦產(chǎn)資源開發(fā)手段不完善對其生態(tài)造成了一定破壞。而其夏季容易受到副熱帶高壓影響，氣溫偏高且雨水較多，且全縣山地面積占比較大，容易產(chǎn)生生態(tài)安全隱患。

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

1.2.1 數(shù)據(jù)來源 主要包括DEM高程數(shù)據(jù)、氣象數(shù)據(jù)、土地利用數(shù)據(jù)以及遙感數(shù)據(jù)和人口數(shù)據(jù)。其中2005、2010和2018年3期土地利用數(shù)據(jù)于萬年縣自然資源局與江西省土地整理中心獲取。氣象數(shù)據(jù)在江西省氣象局處獲取，數(shù)據(jù)日期為2005、2010和2018年。DEM影像數(shù)據(jù)由地理空間數(shù)據(jù)云(http://www.gscloud.cn/)下載，影像分辨率為30 m。人口數(shù)據(jù)來自2005—2018年中國縣域統(tǒng)計年鑒及萬年縣各年份統(tǒng)計公報。

1.2.2 構建指標體系 SRP模型從生態(tài)敏感性、生態(tài)恢復力以及生態(tài)壓力度3個層面綜合考慮影響生態(tài)脆弱性的因素。其中，地形、地表和氣象是生態(tài)敏感性層面的主要考慮因子。生態(tài)恢復力層面主要包含植被因子。人類活動是造成生態(tài)壓力的重要來源，因此生態(tài)壓力度層面選擇的是人類活動因子。SRP模型囊括了生態(tài)脆弱性的兩個形成要素，即內部脆弱性和外部脆弱性，選擇SRP評價模型可以較為全面地分析生態(tài)狀況，并為動態(tài)分析打下基礎。本文選取11個具體指標，并進一步劃分為正向和逆向指標，對萬年縣2005—2018年生態(tài)脆弱狀況進行綜合評價，見表1。

圖1 研究區(qū)位Fig.1 Research location

1.2.3 指標計算及標準化 評價基于格網(wǎng)尺度，所有涉及的指標均轉化為格網(wǎng)尺度。研究通過借鑒前人利用格網(wǎng)GIS法對生態(tài)環(huán)境脆弱性的相關研究[21]，并結合萬年縣的實際情況，從萬年縣景觀類型分布以及土地利用進行綜合考慮，選擇1 km×1 km對萬年縣進行格網(wǎng)劃分。將提及的指標數(shù)據(jù)轉化成分辨率為30 m的柵格數(shù)據(jù)，并將各指標值通過ArcGIS 10.5軟件平臺進行分區(qū)統(tǒng)計，得到各格網(wǎng)的指標值。

各個指標統(tǒng)一量綱，以便指標間的統(tǒng)計和對比。采納極差標準化法，見下式：

(1)

表1 指標體系

表2 2005—2018年各指標組合權重

(2)

式中，Yij表示標準化后各指標的值，Xij代表第i個指標在第j個格網(wǎng)的原始數(shù)據(jù)值，Xjmax則是指標在格網(wǎng)中的最大值，同理，Xjmin代表指標在格網(wǎng)中的最小值。

1.2.4 確定指標權重 層次分析法運用范圍廣，但受主觀影響較大，而熵權法則是以各指標中所能提供的信息熵來確定各權重，指標所含信息熵越小，其變異程度越大則賦予的權重越大，反之亦然[22]。以熵權法對層次分析法的結果進行修正，得出各指標的綜合權重，公式如下：

(3)

式中，ari為層次分析法計算的權重，Sri為熵權法計算的權重，Ri為組合權重。計算2005、2010、2018年萬年縣生態(tài)脆弱性，各指標權重見表2。

1.2.5 綜合指數(shù)法 采用綜合指數(shù)法求得生態(tài)脆弱性指數(shù)EVI，計算公式如下：

(4)

式中，EVI表示生態(tài)脆弱性指數(shù)，Wi為各個指標的標準化值，而Ri是各個指標所占的權重。同時，對生態(tài)脆弱性指數(shù)進行統(tǒng)一標準化的計算，便于其數(shù)值比照，公式如下：

(5)

式中，SEVI代表標準化后的EVI值，EVI為全部格網(wǎng)單元的生態(tài)脆弱性指數(shù)值，EVImax為其中的最大值，EVImin代表其中的最小值。在得到標準化后的生態(tài)脆弱性指數(shù)后，使用自然斷點法[3]，生態(tài)脆弱性最后劃分為微度脆弱、輕度脆弱、中度脆弱、重度脆弱以及極度脆弱5個層級。

1.2.6 空間自相關分析 以Moran’ I指數(shù)和聚類圖表征萬年縣2005—2018年生態(tài)脆弱性的空間關聯(lián)關系，研究通過Geoda 1.12軟件繪制萬年縣局部空間自相關聚類圖，以分析萬年縣生態(tài)脆弱性高值區(qū)、低值區(qū)的空間分異規(guī)律。

1.2.7 OWA方法 OWA方法由美國數(shù)學家Yager率先提出，其核心在于對指標按照屬性的重要性重新排序，對不同位序的指標賦予不同的次序權重。關于OWA權重的確定方法眾多，其中單調規(guī)則遞增(RIM)法運用廣泛，其計算公式如下：

(6)

式中，vj為次序權重，vj∈[0,1]，且v1+v2+…+vn=1。a為決策風險系數(shù)，取值0到∞之間；Wk為指標重要性等級，其計算公式如下：

(7)

式中，n表示指標個數(shù)，rk根據(jù)指標數(shù)值大小對其進行重要性取值，最大取1，最小取n。

利用OWA方法計算生態(tài)脆弱性空間多準則評價的關鍵在于將準則權重和次序權重進行線性組合從而得出評價結果[1]，其公式如下：

(8)

式中，zij為第i個像元中第j項指標的屬性值；ui為準則權重，ui∈[0,1]，且u1+u2+…+un=1,研究中使用AHP方法計算準則權重；vj為次序權重。

本文基于單調規(guī)則遞增(RIM)法計算7種不同決策風險系數(shù)下的次序權重，結果見表3，在決策者持樂觀態(tài)度時(a<1),賦予重要性較高的指標更高的次序權重；在決策者持悲觀態(tài)度時(a>1)，重要性較低的指標被優(yōu)先考慮賦予更高的次序權重，重要性較高的反而不被重視；當決策者無偏好時(a=1)，只用準則權重進行綜合評價。

2 結果與分析

2.1 萬年縣生態(tài)脆弱性空間分布特征

由圖2可知，萬年縣生態(tài)脆弱性在2005—2010年間大體呈現(xiàn)西北低、東南高的現(xiàn)象，而2018年其生態(tài)脆弱性呈現(xiàn)中部低、北部及東南高的狀態(tài)，且2005—2018年間極度脆弱區(qū)集中在陳營鎮(zhèn)。

由表4可知，2005—2018年均以輕度脆弱為主，其中2005年萬年縣輕度脆弱面積為528.56 km2，占據(jù)其總面積45.96%，主要分布在齊埠鄉(xiāng)、石鎮(zhèn)鎮(zhèn)、汪家鄉(xiāng)、青云鎮(zhèn)、蘇橋鄉(xiāng)、珠田鄉(xiāng)及梓埠鎮(zhèn)南部區(qū)域。微度脆弱區(qū)主要分布于珠田鄉(xiāng)和石鎮(zhèn)鎮(zhèn)的中部區(qū)域，面積占比為28.12%，而萬年縣中度、重度脆弱區(qū)主要分布在陳營鎮(zhèn)周邊以及大源鎮(zhèn)南部、裴梅鎮(zhèn)東南部地區(qū)，其面積分別為203.19、88.98 km2，占比分別為17.67%、7.74%。而極度脆弱在萬年縣分布集中但面積不大，僅占0.52%。

與2005年相比，萬年縣2010年輕度脆弱面積有所下降，面積為473.77 km2，占據(jù)縣域面積41.19%，其中度脆弱和重度脆弱、極度脆弱分布更為集中，主要在萬年縣東南部各鄉(xiāng)鎮(zhèn)，且面積也有所增加，分別為213.52、90.93、8 km2，各自占縣域面積18.56%、7.91%和0.70%。而2018年，萬年縣中度、重度脆弱區(qū)相比2005、2010年有向北部聚集的趨勢，其中中度脆弱區(qū)主要集中在陳營鎮(zhèn)周邊以及萬年縣北部的梓埠鎮(zhèn)、湖云鄉(xiāng)等地區(qū)，重度脆弱區(qū)主要集中在大源鎮(zhèn)和裴梅鎮(zhèn)，且在北部占據(jù)零星地塊。中度脆弱區(qū)面積為237.57 km2，占比為20.66%，相比2005和2010年呈增加趨勢。另外，微度和極度脆弱變化不大，極度脆弱區(qū)仍集中于陳營鎮(zhèn)中心，亦為萬年縣的城區(qū)所在地，面積占比為0.80%。

2005、2010和2018年的生態(tài)脆弱性指數(shù)的均值分別為0.2633、0.2473、0.1914，呈逐年遞減趨勢，表征萬年縣生態(tài)脆弱性在研究期間內逐漸轉好。2018年中度脆弱區(qū)集中于北部耕地地區(qū)，而隨著社會經(jīng)濟的發(fā)展、人口的不斷增加，對糧食的需求也在不斷增加，人類活動增強使得生態(tài)壓力增強，造成輕度脆弱逐漸演變成中度脆弱，并出現(xiàn)零星重度脆弱區(qū)。而2005—2010年重度脆弱區(qū)集中于大源、裴梅等礦區(qū)開發(fā)地，2010年之后，萬年縣對廢棄礦區(qū)開展治理等工作，至2018年，重度脆弱區(qū)面積下降為86.81 km2，占7.55%。如持續(xù)進行生態(tài)修復等工作，其生態(tài)狀況可進一步改善。

2.2 萬年縣鄉(xiāng)鎮(zhèn)尺度生態(tài)脆弱性變化特征

由表5可知，2005—2018年，齊埠鄉(xiāng)、青云鎮(zhèn)、上坊鄉(xiāng)、石鎮(zhèn)鎮(zhèn)、蘇橋鄉(xiāng)、汪家鄉(xiāng)和珠田鄉(xiāng)均以微度、輕度脆弱區(qū)為主，均位于萬年縣丘陵區(qū)，植被覆蓋豐富，生態(tài)恢復力較強。2005—2010年間，大源鎮(zhèn)鎮(zhèn)域生態(tài)脆弱性明顯加重，從中度脆弱轉變?yōu)橹囟却嗳酰?010—2018年間，又逐步轉為中度脆弱，表明其社會經(jīng)濟活動已經(jīng)對生態(tài)環(huán)境造成了巨大威脅，對此，政府大力推進生態(tài)文明建設，其成果初步顯現(xiàn)，環(huán)境由差轉好，生態(tài)保護修復得到高度重視。隨著經(jīng)濟、人口的快速增長，耕地資源利用加強，造成湖云鄉(xiāng)、梓埠鎮(zhèn)這兩個以耕地為主的鄉(xiāng)鎮(zhèn)由輕度脆弱轉變?yōu)榱酥卸却嗳?。萬年縣主城區(qū)由于土地利用程度較高、人口壓力等對生態(tài)壞境造成一定破壞，始終為主要的極度脆弱區(qū)。而裴梅鎮(zhèn)雖然在研究期間主要以輕度脆弱為主，但其山地眾多，存在明顯的自然本底脆弱性，生態(tài)問題不容忽視。

表3 生態(tài)脆弱性評價指標次序權重

圖2 萬年縣2005—2018生態(tài)脆弱性分級Fig.2 Classification of ecological vulnerability in Wannian county from 2005 to 2018

表4 生態(tài)脆弱性面積占比

表5 各鄉(xiāng)鎮(zhèn)生態(tài)脆弱性面積占比

圖3 萬年縣生態(tài)脆弱性指數(shù)2005—2018年Moran’I散點 Fig.3 Moran scatter of ecological vulnerability index in Wannian county from 2005 to 2018

2.3 萬年縣生態(tài)脆弱性空間集聚特征

經(jīng)過空間自相關的探討，計算出全局Moran’I指數(shù)分別為0.504 848、0.680 825和0.742 960，表明2005、2010和2018年萬年縣生態(tài)脆弱性具有明顯的空間關聯(lián)關系(圖3)，2005—2018年生態(tài)脆弱性的空間自相關性逐漸增加，呈現(xiàn)高度空間正相關關系。

圖4顯示，2005年萬年縣生態(tài)脆弱性高高集聚區(qū)集中于東南部。鄉(xiāng)鎮(zhèn)尺度上，萬年縣生態(tài)脆弱性的高高集聚區(qū)主要分布在中心城區(qū)、裴梅鎮(zhèn)、大源鎮(zhèn)，基本分布在極度脆弱區(qū)和重度脆弱區(qū)。而低低聚集區(qū)主要集中在微度脆弱區(qū)，少部分零星分布在輕度脆弱區(qū)。而2010、2018年高高集聚區(qū)相比2005年更為集中，并逐漸向中心城區(qū)、縣域東南部靠攏，且2018年的高高集聚區(qū)有向北部拓展的趨勢。大源鎮(zhèn)、裴梅鎮(zhèn)山地較多，且礦產(chǎn)資源開發(fā)、廢棄資源處理方式的不完善，形成了生態(tài)脆弱性較高的區(qū)域以這兩個鎮(zhèn)為主向四周輻射的狀態(tài)，是高高聚集的重要原因。而低低集聚區(qū)在研究期間內變化較大，其分布狀況由零散逐漸變得集中，至2018年，基本分布在微度脆弱區(qū)?？傮w而言，以目前生態(tài)脆弱性的空間集聚效應來看，需要進一步加強對高高集聚區(qū)的管控，尤其要防范其向周邊地區(qū)輻射擴展效應。

2.4 基于OWA的萬年縣生態(tài)脆弱性情景模擬

以2018年萬年縣生態(tài)脆弱性評價結果為基礎，結合IDRISI軟件中的OWA模塊，將11個指標的標準化柵格圖、準則權重、次序權重輸入其中，得到不同決策風險下的萬年縣生態(tài)脆弱性評價結果。結合實地調查萬年縣社會經(jīng)濟發(fā)展水平以及生態(tài)文明建設情況，研究設置了3種情景模擬萬年縣生態(tài)脆弱性：a=1時，為維持現(xiàn)狀型；a=0.8時，為生態(tài)可持續(xù)發(fā)展的樂觀決策情景；a=3時，為生態(tài)不可持續(xù)的悲觀決策情景。

基于圖5分析3種決策情境下萬年縣的生態(tài)脆弱性，當a=0.8時，萬年縣全縣處于微度脆弱和輕度脆弱，極少數(shù)處于中度脆弱，其中微度脆弱占比5.62%、輕度脆弱占比93.40%。a=1時，萬年縣生態(tài)脆弱性狀況接近使用AHP-熵權法時的評價結果，全縣以輕度脆弱和中度脆弱為主，面積占比分別為63.29%、33.03%，表征OWA方法在評價生態(tài)脆弱性具有一定適用性。a=3時，萬年縣以中度脆弱為主，占比59.37%，且重度脆弱占比1.64%，除縣域西北部生態(tài)脆弱性較低，其余地區(qū)生態(tài)脆弱程度均較高。此時，高程、起伏度、均溫、降水等指標次序權重較高，萬年縣東南部原有自然因素成為造成生態(tài)環(huán)境脆弱的重要原因，生態(tài)環(huán)境問題不容忽視。

圖4 萬年縣生態(tài)脆弱性指數(shù)2005—2018聚類Fig.4 Clustering of ecological vulnerability index in Wannian county from 2005 to 2018

圖5 不同政策情境下萬年縣生態(tài)脆弱性狀況Fig.5 Ecological vulnerability in Wannian county under different policy scenarios

生態(tài)持續(xù)型情景下，全縣輕度脆弱面積超過90%，生態(tài)環(huán)境保護得到極大重視，此時全縣生態(tài)被破壞的區(qū)域幾乎沒有，決策策略可以適當考慮加大社會發(fā)展的投入，以提高未來生態(tài)修復的經(jīng)濟保障。維持現(xiàn)狀型情景下，萬年縣生態(tài)脆弱性評價結果趨近于反映區(qū)域現(xiàn)實脆弱性，此時縣域中心城區(qū)、西北耕作區(qū)脆弱程度較高，決策策略可以考慮在這些區(qū)域開展生態(tài)修復試點工作。生態(tài)不可持續(xù)型情景下，全縣脆弱程度較高，中度脆弱占比接近六成，且集中于東南方，此時決策策略需要重點恢復縣域生態(tài)環(huán)境，由于東南山區(qū)地勢較高、地形起伏較大，決策者可以考慮從水源涵養(yǎng)、生物多樣性維護、水土流失治理等方面制定修復措施，通過自然修復與人為干預的方式進行綜合治理。

3 討 論

運用SRP模型從生態(tài)敏感性、生態(tài)恢復力、生態(tài)壓力度3個層次遴選指標，較為綜合地分析了2005—2018年間萬年縣生態(tài)脆弱性時空演變特征，同時對其進行空間自相關分析，探討生態(tài)脆弱性在空間上的集聚狀況，期望為生態(tài)治理修復、分區(qū)管控提供參考。但由于研究尺度較小，部分數(shù)據(jù)獲取受到限制，且生態(tài)脆弱性分級方式、權重確定等多樣化，并未有一個統(tǒng)一標準，因此生態(tài)脆弱性評價體系值得進一步完善。

通過OWA方法設置不同決策風險，對生態(tài)脆弱性進行不同情景的模擬，有助于決策者更全面、科學地進行判斷。在對萬年縣進行生態(tài)脆弱性評價時，不同決策風險系數(shù)下的結果也各不相同，決策者通過得到多情景組合情況下的評價結果進而統(tǒng)籌社會經(jīng)濟發(fā)展和風險管控之間的關系。值得注意的是，OWA方法的評價結果只能為處在動態(tài)變化之中的區(qū)域發(fā)展提供決策思路，并不代表方案唯一優(yōu)解。本文基于對生態(tài)脆弱性的認識，嘗試探索不同決策情景下的脆弱性結果，期望為決策者提供不同的決策策略，同時彌補多準則決策在生態(tài)脆弱性研究中較少的情況。

4 結 論

采用綜合指數(shù)法對萬年縣2005—2018年的生態(tài)脆弱性進行評價，并基于OWA方法探討2018年不同情景方案下的國土空間脆弱性。

(1)2005—2018年，萬年縣生態(tài)脆弱性呈明顯下降趨勢，生態(tài)脆弱性指數(shù)均值由2005年的0.2633下降至2018年的0.1914，累積下降27.31%。區(qū)域生態(tài)脆弱性總體呈現(xiàn)“西北高—中部、西部低—東南高”的空間格局。

(2)2005—2018年，萬年縣生態(tài)脆弱性Moran’ I指數(shù)分別為0.504 848、0.680 825、0.74 296，表征萬年縣在研究期間生態(tài)脆弱性呈現(xiàn)高度正相關關系，空間集聚進一步增強。

(3)基于OWA方法模擬不同決策風險下萬年縣生態(tài)脆弱性，結果表明隨著風險系數(shù)的增加，萬年縣生態(tài)脆弱性逐漸升高，決策者態(tài)度由樂觀變?yōu)楸^。萬年縣處于“生態(tài)持續(xù)型”、“生態(tài)不可持續(xù)型”情景下的生態(tài)脆弱性分別處于輕度脆弱(93.40%)、中度脆弱(59.37%)水平。

(4)“維持現(xiàn)狀型”時，全縣以輕度脆弱和中度脆弱為主，面積占比分別為63.29%、33.03%，接近使用AHP-熵權法時的評價結果，表明OWA方法在生態(tài)脆弱性研究中具有一定適用性。