楊振華,周秋文,郭 躍,蘇維詞,4*,張鳳太

?

基于SPA-MC模型的巖溶地區(qū)水資源安全動態(tài)評價——以貴陽市為例

楊振華1,周秋文2,郭 躍3,蘇維詞3,4*,張鳳太5

(1.貴州師范大學(xué)喀斯特研究院,貴州貴陽 550001；2.貴州師范大學(xué)地理與環(huán)境學(xué)院,貴州貴陽550001；3.重慶師范大學(xué)地理與旅游學(xué)院,重慶 400047；4.貴州科學(xué)院山地資源研究所,貴州貴陽550001；5.貴州師范學(xué)院資源環(huán)境與災(zāi)害研究所,貴州貴陽 550003)

為合理評價巖溶地區(qū)水資源安全態(tài)勢,基于PESBR概念模型構(gòu)建涵蓋巖溶地區(qū)工程性缺水特性的評價指標(biāo)體系,并采用SPA-MC(集對分析--馬爾科夫鏈)耦合模型對2002~2014年貴陽市水資源安全狀況進(jìn)行動態(tài)評價.研究表明:貴陽市水資源安全除2008、2012、2014年屬“較安全”等級外,其余年份為“臨界安全”等級,受“安全”等級隸屬度1和累計等級聯(lián)系度2逐漸上升的影響,貴陽市水資源安全等級逐漸向“較安全”的趨勢轉(zhuǎn)變;依據(jù)貴陽市2015~2050年水資源安全預(yù)測結(jié)果可知,貴陽市未來水資源安全將長期處于“較安全”狀態(tài),其中,“安全”等級隸屬度1在2015~2050年期間從0.365持續(xù)上升至0.435,而“不安全”等級5則由0.194下降至0.138,年均下降率達(dá)1.93%;實例證明SPA-MC耦合模型對水資源動態(tài)評價結(jié)果合理、客觀.

SPA-MC；巖溶；水資源安全；動態(tài)評價

巖溶地區(qū)是一個由雙重含水介質(zhì)組成的“二元三維”空間結(jié)構(gòu)系統(tǒng),具有特殊的地貌-水文結(jié)構(gòu)與動態(tài)過程,地表儲水能力差,可利用、便于利用的水資源相對有限[1-3].同時,由于城市人口不斷增長,工農(nóng)業(yè)迅速發(fā)展,大量生活、生產(chǎn)污水的排放,嚴(yán)重破壞水環(huán)境質(zhì)量,加劇區(qū)域用水困難,導(dǎo)致水資源成為許多巖溶地區(qū)社會、經(jīng)濟發(fā)展的重要制約因素,水資源安全問題也逐漸成為巖溶水資源研究領(lǐng)域的熱點[4-6].因此,構(gòu)建科學(xué)評價指標(biāo)體系,運用合適的巖溶地區(qū)水資源安全評價模型,開展巖溶地區(qū)水資源安全評價研究,對制定巖溶地區(qū)城市水資源安全利用規(guī)劃、保障城市水資源安全尤為重要.

目前,關(guān)于水資源安全評價的研究成果主要集中在評價指標(biāo)體系構(gòu)建、評價模型設(shè)計和選取等方面,如從RSR、DPSIR、DPSIRM概念模型的角度構(gòu)建指標(biāo)體系[7-9],選取WEAP模型[10]、CN-AM模型[11]、投影尋蹤模型[12]、系統(tǒng)動力學(xué)模型[13]、SPA模型[14]、支持向量機[15]、聚類分析法[16]等眾多評價模型、方法應(yīng)用于水資源安全評價中,而這些模型、方法在巖溶地區(qū)水資源安全評價中也取得了一定的研究進(jìn)展,如張鳳太等構(gòu)建熵權(quán)集對耦合評價模型對貴州省水資源安全進(jìn)行整體評價[17];鄒勝章等提出巖溶地區(qū)水資源安全評價應(yīng)涵蓋水量、水質(zhì)、水資源利用以及水災(zāi)害防治等內(nèi)容[18];段琪彩等結(jié)合巖溶地形地貌,水循環(huán)時空特征和人口經(jīng)濟狀況分析出水資源安全的影響因素[19];不難看出,已有的巖溶地區(qū)水資源安全評價研究多是直接借鑒了傳統(tǒng)的評價指標(biāo)體系,使所選取的指標(biāo)體系不能體現(xiàn)巖溶地區(qū)工程性缺水特性,更不能突顯巖溶水資源利用特征,導(dǎo)致其評價結(jié)果的客觀性也有待商榷.另外,受評價模型的限制,其評價過程多注重歷年水資源安全狀態(tài)而忽略各影響因素的動態(tài)變化過程,對歷年數(shù)據(jù)的線性趨勢分析不能突出水資源安全的變化特征,對水資源安全的預(yù)測研究也比較薄弱.

模糊集對與馬爾科夫鏈的耦合評價模型(SPA-MC)是利用集對聯(lián)系度衡量水資源安全系統(tǒng)內(nèi)的確定性與不確定性,然后,對聯(lián)系度的穩(wěn)態(tài)轉(zhuǎn)化規(guī)律來預(yù)測各序列狀態(tài)的變化趨勢,從而實現(xiàn)對水資源安全的動態(tài)評價與預(yù)測.其中,模糊集對分析(SPA)將水資源安全系統(tǒng)因素的不確定性和確定性作為一個系統(tǒng)進(jìn)行綜合考察[20-21].馬爾科夫鏈(MC)主要通過研究各水資源安全狀態(tài)的初始概率及各狀態(tài)之間的轉(zhuǎn)移概率,來預(yù)測各序列狀態(tài)的變化趨勢,具有對于短期水資源安全水平預(yù)測較為準(zhǔn)確的優(yōu)勢[22].

基于此,本文根據(jù)巖溶地區(qū)“三水”資源轉(zhuǎn)化迅速, “工程性”缺水顯著、工農(nóng)業(yè)用水效率低下等水資源系統(tǒng)特征,依照PSR模型原理,構(gòu)建基于PESBR概念模型的評價指標(biāo)體系,采用SPA-MC模型動態(tài)評價貴陽市2002~2014年水資源安全狀況,并預(yù)測其2015~2050年變化趨勢,為協(xié)調(diào)巖溶地區(qū)水資源可持續(xù)利用與區(qū)域經(jīng)濟發(fā)展提供科學(xué)支撐和決策依據(jù).

1 水資源安全評價模型

本文依據(jù)水資源安全的內(nèi)涵,巖溶地區(qū)水資源供需利用特征,設(shè)計出水資源安全PESBR概念模型(圖1).其中, “工程性”缺水(E)主要表現(xiàn)在水資源獲取基礎(chǔ)性條件,保障城鎮(zhèn)用水的主要水庫蓄水率和滿足農(nóng)田用水的有效灌溉面積比等3方面(表1).相比于非巖溶地區(qū),巖溶地區(qū)地表產(chǎn)流系數(shù)較低,地表、地下水漏失嚴(yán)重,地形起伏度大,導(dǎo)致水資源獲取條件差,其獲取方式以高差提水和長距離引水為主,供水成本高;巖溶地質(zhì)、地貌復(fù)雜,地層不穩(wěn)定性高,水利設(shè)施建設(shè)難度大,工程滲漏問題嚴(yán)重,蓄水工程的蓄水率低,城市化進(jìn)程中水源保障能力有限;巖溶地區(qū)耕地保水能力弱,農(nóng)田有效灌溉絕大部分需農(nóng)業(yè)水利工程的支持,故用有效灌溉面積代表代表農(nóng)田水利工程的完善程度,可突顯“工程性”缺水對農(nóng)業(yè)發(fā)展的影響.

1.1 指標(biāo)體系構(gòu)建

根據(jù)PESBR概念模型內(nèi)涵,建立涵蓋需水壓力(P)、工程性缺水(E)、承載狀態(tài)(S)、生態(tài)基礎(chǔ)(B)、人為響應(yīng)(R)等5個評價準(zhǔn)則27個評價指標(biāo)的巖溶地區(qū)水資源安全評價指標(biāo)體系.根據(jù)指標(biāo)屬性分成正向型與負(fù)向型指標(biāo)(表1),對不同屬性指標(biāo)分別進(jìn)行賦權(quán)評價.

表1 巖溶地區(qū)水資源安全評價指標(biāo)體系及釋義

注 :X＋/－(=1,2,…,27)中的“＋/－”代表該項指標(biāo)屬性為正向性/負(fù)向性,“(-)”代表該指標(biāo)無量綱.

1.2 指標(biāo)賦權(quán)

在構(gòu)建評價指標(biāo)體系的基礎(chǔ)上,融合AHP法與變異系數(shù)法[26]對各指標(biāo)進(jìn)行組合賦權(quán),實現(xiàn)主觀偏好與客觀相關(guān)的結(jié)合權(quán)重[27],其值見表1.

1.3 指標(biāo)分級

考慮到巖溶地區(qū)生態(tài)系統(tǒng)脆弱性和“三水”轉(zhuǎn)化規(guī)律性,參照已有水資源安全標(biāo)準(zhǔn)的相關(guān)文獻(xiàn)[7-9,13-16],結(jié)合有關(guān)水資源安全指標(biāo)臨界值以及國內(nèi)外政府頒布的標(biāo)準(zhǔn)和規(guī)劃目標(biāo),并根據(jù)專家意見將安全標(biāo)準(zhǔn)做出適當(dāng)調(diào)整,將水資源安全評價等級劃分為“安全、較安全、臨界安全、較不安全、不安全”5個等級,并確定評價指標(biāo)分級標(biāo)準(zhǔn)(表2).

表2 水資源安全評價指標(biāo)等級劃分標(biāo)準(zhǔn)

1.4 評價模型

1.4.1 模糊集對分析 SPA理論將評價指標(biāo)體系與安全等級當(dāng)作一個集對,根據(jù)指標(biāo)體系X和評價等級B構(gòu)造集對,得到=5的多元聯(lián)系度[28],即:

令

則公式(1)變成:

(2)

對于負(fù)向性指標(biāo),指標(biāo)值x與該指標(biāo)級評價標(biāo)準(zhǔn)的聯(lián)系度如下[28]:

(3)

對于正向性指標(biāo),指標(biāo)值x與該指標(biāo)級評價標(biāo)準(zhǔn)的聯(lián)系度(x,B)如下:

根據(jù)參考文獻(xiàn)[28],本文采用置信度準(zhǔn)則來確定樣本的等級.即:

(5)

式中:為置信度,其取值范圍一般建議為[0.5,0.7];h為樣本的元聯(lián)系度.

1.4.2 馬爾科夫鏈 設(shè)是一概率空間, {(),30}為定義在概率空間內(nèi)的整數(shù)隨機序列,若對于任意的31,當(dāng)1,2,…,t(其中1, <2<,…,<t)時刻,對應(yīng)()的觀測(1), C(2),…,C(t)滿足條件[29-30]:

((t)|(t-1),(t-2),…,(1))=((t)|(t-1)) (6)

則稱{(),30}為馬爾科夫鏈(Markov Chain).馬爾科夫鏈性質(zhì)表明:t時刻的{(),30}只與t-1時刻的數(shù)值有關(guān),有其他時刻的數(shù)值無關(guān),((t)|(t-1))即為t-1與t時段的轉(zhuǎn)移概率.設(shè)時間狀態(tài)集合=1,2,3,…,S,則將=((t)|(t-1))稱為轉(zhuǎn)移概率矩陣.

在[t,t+D]時段內(nèi),評估指標(biāo)的安全等級通常會因D不同而不同.由此,假設(shè)在t時刻有S個評價指標(biāo)的安全等級為非常安全,到t+D時有S1個評價指標(biāo)為非常安全等級,有S2、S3、S4、S5個評價指標(biāo)等級分別轉(zhuǎn)變成安全、臨界安全、不安全和極不安全,且=1+2+3+4+5.故S個安全等級評價指標(biāo)在[t,t+D]周期內(nèi)的狀態(tài)轉(zhuǎn)移概率向量為:

1=(11,12,13,14,15)=

式中:11+12+13+14+15=1,為[t,t+D]周期內(nèi)的各轉(zhuǎn)移方向的權(quán)重之和.同理,可以得到2、3、4、5對應(yīng)的轉(zhuǎn)移概率向量,進(jìn)而可得出評價指標(biāo)體系在[t,t+D]期間的狀態(tài)轉(zhuǎn)移概率矩陣為:

(8)

通過上述SPA確定的安全等級聯(lián)系度和狀態(tài)轉(zhuǎn)移概率矩陣,采用隸屬度的原理建立水資源安全評價模型,即[t,t+D]時期內(nèi),水資源安全評價值為:

~B(t+D)=(t+D)+(t+D)1+(t+D)2+(t+D)3+(t+D)=[(t),(t),(t),g(t),(t)]×

PD×(1,1,2,3,)(9)

式(9)中,(t)、(t)、(t)、(t)、(t)分別代表t時刻原始同一度、對立度分量集和差異度;PDt為D時段后的概率轉(zhuǎn)移矩陣;(1,1,2,3,)為同一度、對立度和差異度系數(shù).

由MC的遍歷性可知,轉(zhuǎn)移概率矩陣PDt符合A.N.柯爾莫哥洛夫方程[31],即隨著變化周期的遞增,PDt將趨于穩(wěn)定.因此,時刻的安全狀態(tài)評估值在經(jīng)過多個變化周期的轉(zhuǎn)移矩陣后,最終達(dá)到穩(wěn)定態(tài)勢.同時,考慮到聯(lián)系度歸一化的性質(zhì),則可由方程組(10)得到水資源安全評價穩(wěn)態(tài)值.

2 實例分析

2.1 研究區(qū)概況

貴陽市地處云貴高原的東部,屬我國東部丘陵向西部高原過渡地帶.地形、地貌走勢大致呈東西向延展,南北高,中部低,地形起伏度大.巖溶地貌發(fā)育廣泛,高原山地、丘陵、盆地、河谷等地貌分布廣泛,海拔介于872~1659m之間.貴陽市下轄六區(qū)一市三縣,總面積8034km2,市區(qū)面積2403km2.貴陽市常年受西風(fēng)帶控制,屬于亞熱帶濕潤溫和型氣候,并有明顯的高原季風(fēng)氣候特征,極端最高最低溫分別為32.5℃和-4.1℃,年平均氣溫為15.3℃,年平均相對濕度為78%,年平均降水量為1129.5mm,年徑流量45.2×108m3,徑流密度達(dá)12.4km/km2.

2.2 數(shù)據(jù)來源

研究數(shù)據(jù)主要來源于《貴陽市水資源公報》(2002~2014),《貴陽市環(huán)境狀況公報》(2002~2014)以及《貴陽市統(tǒng)計年鑒》、《貴陽年鑒》、《國民經(jīng)濟和社會發(fā)展統(tǒng)計公報》(2002~2014)等資料整理.部分指標(biāo)涉及其他數(shù)據(jù)來源簡述如下:降水季節(jié)變異性采用貴陽市國家氣象站點逐日數(shù)據(jù)中計算得到;地表徑流密度(13)來源于貴陽市水系圖(1:50000)提取;巖溶石漠化率(20)一律采用國土部門提供的數(shù)據(jù),缺值年份通過該年度石漠化治理面積得到;另外,工業(yè)廢水達(dá)標(biāo)率(24)和生活污水處理率(26)最近兩年的數(shù)據(jù)根據(jù)以往年份的平均增長率插值得到.

2.3 評價結(jié)果與分析

2.3.1 水資源安全評價 根據(jù)式(3~4)構(gòu)造水資源安全評價因子與評價等級集對分析集合,并作同一性(1)、差異性(2-4)、對立性(5)等聯(lián)系度分析,得到貴陽市2002~2014年水資源安全評價指標(biāo)的聯(lián)系度,然后,將各指標(biāo)的聯(lián)系度與權(quán)重相結(jié)合,得到貴陽市2002~2014年歷年水資源安全水平與評價等級的聯(lián)系度,取=0.5為界,3=1+2+3=0.55>,由置信度準(zhǔn)則[32]可以判斷2002年貴陽市水資源安全為臨界安全等級,同理,可以計算出2002~2014年貴陽市水資源安全等級(表3).

由表3可知,貴陽市2002~2014年水資源安全等級除了2005、2012、2014年為較安全外,其余年份均為臨界安全水平,歷年水資源安全水平變化較小,且呈現(xiàn)出波動式上升的趨勢,累計聯(lián)系度H和H與置信度的距離也有逐漸變小的趨勢.這表明貴陽市在供水能力和用水效率上有較大提升,這也是貴陽市生態(tài)文明城市建設(shè)進(jìn)程中,水利設(shè)施不斷完善,產(chǎn)業(yè)結(jié)構(gòu)不斷優(yōu)化,生態(tài)環(huán)境趨于良性循環(huán)的結(jié)果.結(jié)合圖2所示,2002~2008年、2008~2012年和2012~2014期間,分別以2004年、2010年和2013年為分界點,H皆呈現(xiàn)出現(xiàn)出先上升后下降的過程,這與貴陽市供水能力不足,人口聚集,經(jīng)濟迅速發(fā)展,需水量大且密集的水資源利用特征相關(guān),如2009~2010年受西南大旱影響,水資源缺乏,導(dǎo)致2010年左右全市水資源安全降為臨界安全狀態(tài).

表3 貴陽市2002~2014年水資源安全等級聯(lián)系度

2.3.2 水資源安全動態(tài)預(yù)測 水資源利用系統(tǒng)的不確定性和社會經(jīng)濟因素的年際差異性,導(dǎo)致水資源安全評價指標(biāo)的歷年變異系數(shù)值較高(Cv=0.235),屬于離散時間序列,符合馬爾科夫過程.文章在各指標(biāo)等級聯(lián)系度的基礎(chǔ)上,利用馬爾科夫鏈(MC)對貴陽市2015~2020年水資源安全趨勢進(jìn)行預(yù)測.將貴陽市2015~2020年水資源安全評價指標(biāo)等級與權(quán)重相結(jié)合,計算出2002~ 2003、2003~2004、…、2013~2014等各時段的轉(zhuǎn)移向量,代入式(7)中得到2002~2014年水資源安全轉(zhuǎn)移概率矩陣的穩(wěn)態(tài)值(圖4).

由于Dt滿足C-K方程,即經(jīng)過多個D周期后,Dt趨于穩(wěn)定.根據(jù)馬爾科夫過程的遍歷性,代入式(10)可得到未來年份的貴陽市水資源安全等級聯(lián)系度的預(yù)測值(表4)以及貴陽市水資源安全的穩(wěn)態(tài)聯(lián)系度:(,)=[0.2006, 0.2001, 0.1999, 0.1998, 0.1995],處于臨界安全狀態(tài),與歷年貴陽市的水資源安全系統(tǒng)的整體態(tài)勢吻合.另外,文章將2002~2010年作為水資源安全轉(zhuǎn)移概率矩陣模擬年份,2011~2014年為驗證年份,從驗證結(jié)果(表4)可以看出,多元聯(lián)系度的預(yù)測整體精度較高,平均相對誤差僅為0.069,故達(dá)到MC對貴陽市未來年份的水資源安全狀態(tài)預(yù)測精度要求.

從表5可以看出,2015~2050年貴陽市水資源安全等級將長期處于較安全狀態(tài),呈臨界安全向較安全-安全的過渡趨勢.123的等級聯(lián)系度整體上逐年上升,其中1從2015年的0.365上升至2050年的0.435,逐漸向置信度0.5靠近,說明未來將有更多的水資源安全單一指標(biāo)向安全等級變化.將45的等級聯(lián)系度則持續(xù)下降,特別是5的等級聯(lián)系度下降0.056,下降率達(dá)28.86%,表明在人為因素的影響下,水資源安全單一指標(biāo)等級上升明顯.

表4 2011~2014貴陽市水資源安全水平預(yù)測精度驗證

表5 2015~2030年貴陽市水資源安全預(yù)測值(MC)

2.3.3 驅(qū)動因素分析 根據(jù)各個評價準(zhǔn)則的加權(quán)等級聯(lián)系度,可識別出貴陽市2002~2014年水資源安全等級主要驅(qū)動因素類別.結(jié)合表3和圖5可知,貴陽市2002~2014年水資源安全等級主要由P、E、S、B、R等因素類對123的加權(quán)等級聯(lián)系度決定,對高等級聯(lián)系度越高,水資源安全等級也越高,反之則越低.同時,從2002年以來,貴陽市不安全、較不安全要素類越來越少,安全和較安全的要素類越來越多,特別是P、E、S、B、R等3個要素類逐漸從不安全、較不安全向較安全、臨界安全轉(zhuǎn)變.

通過分析貴陽市2002~2014年水資源資源安全評價指標(biāo)體系中各個指標(biāo)的歷年變化特征,及其對水資源安全要素類的貢獻(xiàn)率,對各類水資源安全驅(qū)動因素的主要因子簡述如下:

(1)需水壓力(P)

從需水壓力(P)(圖5a)來看,2002~2003年、2008年、2010年、2013~2014年為較不安全等級,其余年份為臨界安全等級,呈波動下降趨勢,說明在2002~2005年以前經(jīng)濟技術(shù)相對落后,萬元GDP耗水量平均值達(dá)845m3/萬元(圖6),而隨著產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整,新興低耗產(chǎn)業(yè)比重增加,萬元GDP耗水量降至2014年的176m3/萬元,但2011以來人口密度持續(xù)性增加至2014年的567人/km2,導(dǎo)致人均水資源量也由2002年的415m3逐步下降至2014年的232m3,生活用水比重以年均20.81%的增長率上升,促使需水壓力(P)逐漸由臨界安全向較不安全轉(zhuǎn)化.從整體來說,貴陽市面臨的需水壓力將會不斷加大,特別是在城市化水平進(jìn)程加快,居民生活用水量增加的背景下,城市供水壓力會持續(xù)上升.

(2)工程性缺水(E)

從工程性缺水(E)(圖5b)來看,2002~2006年、2011年為較不安全等級,2007~2010年、2012~2013年為臨界安全等級,2014年為較安全等級,整體由較不安全向臨界安全轉(zhuǎn)變的態(tài)勢.貴陽市雖地處巖溶丘陵地帶,地表水資源調(diào)節(jié)能力弱,在“十二五”期間黔中水利骨干水源工程、引提灌溉工程和地下水開發(fā)利用工程建設(shè)以來,城市供水缺口2010年將缺水18萬t/d降至2013年的9~10萬t/d.全市有效灌溉面積由2002年的46%提升至2014年的75%(圖7),農(nóng)田設(shè)施對農(nóng)業(yè)發(fā)展的滿足率不斷提高,使工程性缺水狀況得到一定的改善,但大中型水庫蓄水量受降水影響較大,如2011年降水量僅783mm,距平值達(dá)246mm,導(dǎo)致大中型水庫蓄水量僅有23.68億m3,年末蓄水率為50.27%,較歷年平均值少9億m3,該年份水資源安全等級降低,表明貴陽市仍然屬于工程性缺水城市.提引水比重受地表蓄水量和地形特征的影響較大,一定程度表明貴陽市用水的不便性,其平均值為33.1%,且變化特征不明顯.

(3)承載狀態(tài)(S)

從水資源承載狀態(tài)(S)(圖5c)來看,除2010~ 2011年受“西南大旱”的影響為臨界安全等級外,其余年份均為安全或較安全等級.水資源承載狀態(tài)指標(biāo)可分成3類:

1)降水特征:貴陽市雖地處亞熱帶季風(fēng)氣候區(qū),降水總量豐富,但降水的時空變異性也相當(dāng)顯著,特別是降水季節(jié)變異性平均指數(shù)值0.918,相比空間異質(zhì)性高出0.804,其出入境水量也與產(chǎn)水模數(shù)保持顯著正相關(guān)性(圖8b),表明降雨直接決定區(qū)域出入境水量差.

2)水量特征:貴陽市地形起伏度大,巖溶地貌廣布,導(dǎo)致地表徑流密度大,但地表產(chǎn)流系數(shù)(圖8.a)介于0.37~0.6之間,表明地表徑流量較小,較大一部分的地表水直接通過巖溶裂隙、落水洞、豎井等流入地下,導(dǎo)致地下水比重達(dá)24.91%,接近全國平均水平,地表便于利用的河川徑流量較小.

3)水質(zhì)特征:不管是飲用水還是城市地表水,水質(zhì)皆相對良好,歷年水質(zhì)達(dá)標(biāo)率平均值達(dá)96.43%,其主要原因為貴陽市重工業(yè)、高污染企業(yè)比重小,水體污染物中的COD和氨氮的主要來源于生活污水排放,特別是貴陽市人口密集的中心城區(qū).

(4)生態(tài)基礎(chǔ)(B)

由于生態(tài)基礎(chǔ)(B)(圖5d)的年際變化較小,故其等級聯(lián)系度均為安全等級,對水資源安全系統(tǒng)的影響也較小.貴陽市近10年以來,持續(xù)推進(jìn)石漠化治理工程和封山育林的相關(guān)工作,森林覆蓋率基本呈平穩(wěn)上升的趨勢,使全市石漠化率持續(xù)降低,由2002年的29.87%降至2014年的16.55%(圖9).而對于輸沙模數(shù)來說,其變化幅度較大,主要是受監(jiān)測站點單一的影響.

(5)人為響應(yīng)(R)

從人為響應(yīng)(R)(圖5e)來看,貴陽市工程性缺水難題一直是困擾全市水資源安全水平提升的關(guān)鍵問題,因此,通過人為響應(yīng)加強水利工程建設(shè),提升水資源利用效率,合理開發(fā)與保護(hù)水源成為保障貴陽市供水的重要途徑.貴陽市先后啟動實施了水利建設(shè)“三大會戰(zhàn)”、滋黔水利工程、黔中水利工程、“小康水”行動計劃、水利建設(shè)“三年行動計劃”、水利八大改革、山區(qū)現(xiàn)代水利建設(shè)示范等一系列水利工作,使得病險水庫治理率接近100%,水資源開發(fā)利用率也降低至18.21%.在生態(tài)文明建設(shè)的驅(qū)動下,全市在節(jié)能減排成效集中于萬元GDP污染物減排率和生活污水處理率等方面,其中,GDP污染物減排率平均值達(dá)15.68%,減排成效顯著.生活污水處理率以年均提升1.63%,降低水體污染物對水源的影響(圖10).

綜合上述PESBR概念模型的各驅(qū)動因素分析結(jié)果可知,2002~2014年貴陽市水資源安全等級變化不明顯,但各個因素類年際變化顯著,整體向安全和較安全等級轉(zhuǎn)變,尤其是工程性缺水(E)和人為響應(yīng)(R),表明人為活動對貴陽市水資源安全具有重要意義.

水資源安全評價研究的重難點之一在于如何構(gòu)建評價指標(biāo)體系和分級標(biāo)準(zhǔn).本文根據(jù)巖溶地區(qū)“工程性”缺水特征,結(jié)合國內(nèi)外研究現(xiàn)狀和國家水資源安全利用相關(guān)規(guī)劃,確定各指標(biāo)分級標(biāo)準(zhǔn),彌補指標(biāo)體系未考慮工程性缺水狀況的缺陷,為科學(xué)評價巖溶地區(qū)水資源安全提供新的指標(biāo)體系構(gòu)建依據(jù).但囿于數(shù)據(jù)有限,未能對貴陽市轄區(qū)內(nèi)各區(qū)縣市水資源安全利用特征進(jìn)行研究.另外,對于全球變化環(huán)境下巖溶地區(qū)水資源的安全響應(yīng)機理及效應(yīng)仍有待探討.

3 結(jié)論

3.1 在構(gòu)建貴陽市PESBR評價指標(biāo)體系的基礎(chǔ)上,采用SPA-MC耦合模型對其水資源安全安全動態(tài)評價,由評價結(jié)果可知,貴陽市2002~2014年水資源安全等級為“臨界安全”,歷年評價等級逐漸向“較安全”轉(zhuǎn)變,歷年“安全”等級隸屬度1和累計等級聯(lián)系度2逐漸上升.未來2015~2050年貴陽市水資源安全等級將逐漸由“較安全”向“安全”等級轉(zhuǎn)變,其“安全”、“較安全”和“臨界安全”的等級隸屬度123整體上逐年上升,水資源安全水平逐年提升.

3.2 基于PESBR概念模型的驅(qū)動因子分析表明,解決貴陽市“工程性”缺水主導(dǎo)下的水資源安全問題,不僅需要加大水利設(shè)施建設(shè),提升水利工程蓄水效益,充分利用有效水資源量,還需控制城市人口密度和國土空間開發(fā)強度,降低城市生活用水與廢水排放,提高水資源重復(fù)利用率.

3.3 與常規(guī)的DPSIR、DPSIRM概念模型相比,PESBR概念模型將巖溶地區(qū)“工程性”缺水指標(biāo)納入評價體系中,突顯“工程性”缺水對巖溶地區(qū)水資源安全的影響,特別是在驅(qū)動力分析方面,常規(guī)模型只強調(diào)供、用水量的直接影響因素,而忽略了巖溶地區(qū)水資源獲取的難度和有限性,誤判貴陽市水資源安全處于“危機”狀態(tài),其評價結(jié)果也難免有失偏頗.其次,SPA-MC耦合模型動態(tài)評價水資源的變化特征,與單純評價水資源歷年狀況更具指導(dǎo)性,符合貴陽市長期工程性缺水明顯的現(xiàn)狀.

[1] 楊明德,譚 明,梁 虹.喀斯特流域水文地貌系統(tǒng)[M]. 北京:地質(zhì)出版社, 1998.

[2] 孫 萍,肖飛鵬,黎志鍵.廣西大石山區(qū)干旱災(zāi)害識別與特大干旱成因分析[J]. 中國農(nóng)村水利水電, 2012,(1):91-94.

[3] 李靜文,楊旭宏,李青文.巖溶含水介質(zhì)各向異性綜合探究 [J]. 湖北農(nóng)業(yè)科學(xué), 2014,53(8):1737-1841.

[4] Proclamation of the Minister-class Conference on the 2nd World Water of Man of the United Nations.[R]. Proclamation of the Minister-class conference in Hague, 2000.

[5] David Grey, Claudia W Sadoff. Sink or swim ? Water security for growth and development [J]., 2007,9:545-571.

[6] 夏 軍,朱一中.水資源安全度量:水資源承載力的研究與挑戰(zhàn)[J]. 自然資源學(xué)報, 2002,17(3).

[7] 蘇 印,官冬杰,蘇維詞.基于SPA的喀斯特地區(qū)水安全評價 [J]. 中國巖溶, 2015,34(6):560-569.

[8] 曹 琦,陳興鵬,師滿江.基于DPSIR概念的城市水資源安全評價及調(diào)控 [J]. 資源科學(xué), 2012,34(8):591-1599.

[9] 張鳳太,王臘春,蘇維詞.基于DPSIRM概念框架模型的巖溶區(qū)水資源安全評價 [J]. 中國環(huán)境科學(xué), 2015,35(11):3511-3520.

[10] 王鵬全,方三朝,張永明,等.基于WEAP模型的石羊河水資源安全分析[J]. 中國農(nóng)村水利水電, 2014,(9):92-99.

[11] 金菊良,吳開亞,魏一鳴.基于聯(lián)系數(shù)的流域水安全評價模型 [J]. 水利學(xué)報, 2008,39(4):401-409.

[12] 吳開亞,金菊良.基于變權(quán)重和信息熵的區(qū)域水資源安全投影尋蹤評價模型 [J]. 長江流域資源與環(huán)境, 2011,20(9):1085-1090.

[13] 代 穩(wěn).基于SD模型的水資源安全模擬研究--以貴州省為例 [J]. 水科學(xué)與工程技術(shù), 2010,(4):4-8.

[14] 王宏偉,張 鑫,邱俊楠.模糊集對分析法在水資源安全評價中的應(yīng)用 [J]. 西北農(nóng)林科技大學(xué)(自然科學(xué)版), 2011,39(10): 209-214.

[15] 暢明琦,劉俊萍,馬 惟.基于支持向量機的水資源安全評價 [J]. 自然災(zāi)害學(xué)報, 2011,20(6):167-171.

[16] 高媛媛,王紅瑞,許新宜,等.水資源安全評價模型構(gòu)建與應(yīng)用--以福建省泉州市為例[J]. 自然資源學(xué)報, 2012,27(2):204-214.

[17] 張鳳太,王臘春,蘇維詞,等.基于熵權(quán)集對耦合模型的表層巖溶帶“二元”水資源安全評價 [J]. 水力發(fā)電學(xué)報, 2012,31(6): 70-76.

[18] 段琪彩,黃 英,王 杰.云南省水循環(huán)特性及水安全剖析[J]. 人民長江, 2012,43(13):7-20.

[19] 鄒勝章,朱明秋,唐建生,等.西南喀斯特區(qū)水資源安全與對策 [J]. 地質(zhì)學(xué)報, 2006,80(10):1637-1642.

[20] 趙克勤.集對分析對不確定性的描述和處理 [J]. 信息與控制, 1995,24(3):162-165.

[21] 秦 鵬,孫國政,秦植海,等.基于多元聯(lián)系數(shù)的水環(huán)境安全評價模型 [J]. 數(shù)學(xué)的實踐與認(rèn)識, 2016,46(2):165-171.

[22] 石朝陽.灰色馬爾科夫鏈的改進(jìn)及其應(yīng)用 [D]. 蘭州:蘭州大學(xué), 2014.

[23] 楊 宇.多指標(biāo)綜合評價中賦權(quán)方法評析 [J]. 統(tǒng)計與決策, 2006,(7):17-19.

[24] 暢明琦,劉俊萍,黃 強.水資源安全Vague集多目標(biāo)評價及預(yù)警 [J]. 水力發(fā)電學(xué)報, 2008,27(3):81-87.

[25] 王文圣,李躍清,金菊良,等.水文水資源集對分析區(qū) [M]. 北京:科學(xué)出版社, 2010.

[26] 張一瑫,吳詩輝,劉曉東,等.基于集對分析和馬爾科夫鏈的航空維修安全動態(tài)評估 [J]. 中國安全科學(xué)學(xué)報, 2016,26(1):122- 128.

[27] 許智慧.馬爾科夫狀態(tài)轉(zhuǎn)移概率矩陣的求解方法研究[D]. 哈爾濱:東北農(nóng)業(yè)大學(xué), 2013.

[28] 亢 永.城市燃?xì)饴竦毓艿老到y(tǒng)風(fēng)險研究 [D]. 沈陽:東北大學(xué), 2013.

[29] 王文圣,李躍清,金菊良,等水文水資源集對分析區(qū) [M]. 北京:科學(xué)出版社, 2010.

Dynamic assessment of water resources security in karst area based on SPA-MC model--a case study of Guiyang city.

YANG Zhen-hua1, ZHOU Qiu-wen2, GUO-Yue3, SU Wei-ci3,4*, ZHANG Feng-tai5

(1.Institute of Karst, Guizhou Normal University, Guiyang 550001, China；2.School of Geography and Environment, Guizhou Normal University, Guiyang 550001, China；3.School of Geography and Tourism, Chongqing Normal University, Chongqing 400047, China；4.Institute of Mountain Resources, Guizhou Academy of Sciences, Guiyang 550001, China；5.Institute of Environment, Resources and Disaster, Guizhou Normal College, Guiyang 550018, China).

To reasonable evaluate water resources security situation in karst area, this paper based on the PESBR conceptual model, build evaluation index system of water resources safety which covered engineering water shortage characteristics, then the SPA-MC (Set pair analysis-Markov chain) coupling model was adopted to dynamically evaluated water resources security situation in Guiyang city during 2002~2014. According to the forecast results of the safety of water resources in Guiyang city 2015~2050, Guiyang city in the future water resources security will remain relatively safe state, the "safe" level of membership in the12015~2050 period lasted from 0.365 rose to 0.435, while the "unsafe" level of5decreased from 0.194 to 0.138, with an average annual decline rate up to 1.93%; the example proved that the SPA-MC coupling model to evaluate the water resources dynamic results are reasonable and objective.

SPA-MC；karst；water resources security；dynamic assessment

X824

A

1000-6923(2017)04-1589-12

2016-08-01

貴州省重大科技專項([2012]6015號);國家科技支撐技術(shù)項目(2014BAB03B01);貴州省創(chuàng)新人才團(tuán)隊([2014]4014號)

楊振華(1991-),男,江西贛州人,貴州師范大學(xué)碩士研究生,主要研究方向為水資源開發(fā)與評價.發(fā)表論文10余篇.

* 責(zé)任作者, 研究員, suweici@sina.com

, 2017,37(4)：1589~1600