彭 博

(遼寧省撫順?biāo)木?遼寧 撫順 113015)

供給和需求是影響水資源是否短缺的兩個(gè)主要因素,若供給無法滿足用水需求就會(huì)產(chǎn)生水資源短缺風(fēng)險(xiǎn),而供水和需水受徑流、降雨及其它隨機(jī)因素影響也存在許多不確定因素[1]。因此,我國很多地區(qū)都存在不同程度的水資源短缺風(fēng)險(xiǎn),隨著環(huán)境的變化和經(jīng)濟(jì)的快速發(fā)展,水資源短缺問題逐漸成為制約城市可持續(xù)發(fā)展的關(guān)鍵因素。為實(shí)現(xiàn)經(jīng)濟(jì)、社會(huì)與資源的協(xié)調(diào)發(fā)展,必須全面分析水資源短缺程度,在此基礎(chǔ)上制定科學(xué)的調(diào)控措施與管理方案,以保證水資源安全。

水資源短缺風(fēng)險(xiǎn)評價(jià)是科學(xué)管理水資源的重要依據(jù),許多學(xué)者利用不同方法開展了廣泛研究,常用的評價(jià)方法有多目標(biāo)決策法、主成分法、投影尋蹤法、模糊綜合法、聚類分析法、蟻群算法、遺傳算法和集對分析法等[2-6]。由于水資源短缺風(fēng)險(xiǎn)的評價(jià)因子具有不確定性和模糊性,不同因素之間存在相互影響、互相作用的復(fù)雜關(guān)系,很難鎖定風(fēng)險(xiǎn)系統(tǒng),加之風(fēng)險(xiǎn)的形成是一個(gè)復(fù)雜的變化過程,也很難精確量化影響風(fēng)險(xiǎn)的各類模糊因素[7]。另外,受限于人類社會(huì)發(fā)展水平尚未厘清風(fēng)險(xiǎn)系統(tǒng)內(nèi)部的因果機(jī)理,由于觀測到的風(fēng)險(xiǎn)系統(tǒng)信息不全可能存在未預(yù)見或未知的風(fēng)險(xiǎn)因素,水資源短缺風(fēng)險(xiǎn)信息的不完備性也是風(fēng)險(xiǎn)精準(zhǔn)評價(jià)的瓶頸。因此,文章結(jié)合風(fēng)險(xiǎn)評價(jià)相關(guān)成果和方法,采用組合賦權(quán)法計(jì)算參評因子權(quán)重,通過灰關(guān)聯(lián)法和模糊綜合法評價(jià)沈陽市2008～2021年水資源短缺風(fēng)險(xiǎn),旨在為沈陽市水資源管理規(guī)劃和合理開發(fā)利用提供決策依據(jù)。

1 評價(jià)方法

1.1 組合權(quán)重計(jì)算

參評因子權(quán)重是決定水資源短缺風(fēng)險(xiǎn)評價(jià)科學(xué)合理性的關(guān)鍵環(huán)節(jié),其中主觀賦權(quán)法主要是考慮專家知識、經(jīng)驗(yàn)等計(jì)算參評因子權(quán)重,實(shí)際計(jì)算時(shí)往往未考慮參評因子的信息量值;客觀賦權(quán)法是參評因子所攜帶的數(shù)據(jù)信息為依據(jù),通過計(jì)算離散差異度確定因子權(quán)重,該方法具有較強(qiáng)的客觀性,往往會(huì)引起因子權(quán)重與重要性偏差較高的現(xiàn)象[8-9]。因此,文章綜合利用主、客觀賦權(quán)法計(jì)算參評因子組合權(quán)重,組合公式如下:

(1)

式中:ωm、θm為參評因子m的客觀與主觀權(quán)重。

目前,比較常用的客觀賦權(quán)法是熵權(quán)法。在信息論中熵值能夠反映信息量的大小,衡量信息的無序化程度,參評因子包含的信息越多則其發(fā)揮的決策作用越大,相應(yīng)的熵值就越小。參照熵的定義,利用公式(2)～(4)計(jì)算參評因子m的熵Em、權(quán)重ωm,即:

(2)

(3)

(4)

對于無法完全量化影響因素模糊賦值的情況,比較常用的主觀賦權(quán)法是G1法,經(jīng)專家決策將參評因子集R={u1,u2,…,um}按照重要程度進(jìn)行排序,從而構(gòu)造唯一序關(guān)系:u1>u2>…>uM。然后考慮相鄰指標(biāo)的重要程度計(jì)算參評因子權(quán)重δM,計(jì)算公式為:

δm=θm-1θm-1;m=M,M-1,…,2

(5)

(6)

式中:θk-1=δkθk。

1.2 灰關(guān)聯(lián)分析法

灰關(guān)聯(lián)法是以決策樣本的具體情況為基礎(chǔ),通過計(jì)算分析若干個(gè)比較序列與最優(yōu)序列的差異程度判定是否存在緊密聯(lián)系[10]。灰關(guān)聯(lián)分析相比于其它多因素分析法,評定結(jié)果較為客觀,計(jì)算量小且對原始數(shù)據(jù)要求低,其主要流程如下:

1)步驟1:根據(jù)待決策樣本的實(shí)際情況,確定影響系統(tǒng)行為的比較序列以及反映系統(tǒng)行為特征的參考序列。為準(zhǔn)確評定區(qū)域水資源短缺風(fēng)險(xiǎn),定義水資源需求量與供給量的比值Y為參考序列Y={y(n)|y=1,2,…,N},各影響因子值為比較序列,由n個(gè)決策樣本和m個(gè)因子構(gòu)成矩陣U:U=(umn)M×N。

2)步驟2:采用歸一化公式對極大型和極小型因子進(jìn)行預(yù)處理,將變量序列U=(umn)M×N轉(zhuǎn)變成標(biāo)準(zhǔn)矩陣V=(vmn)M×N,處理公式為:

(7)

式中:ummin、ummax為n個(gè)決策樣本中第m項(xiàng)參評因子的最小和最大值。

根據(jù)最優(yōu)序列X0計(jì)算關(guān)聯(lián)系數(shù)ξmn,其中△mn=|X0-vmn|為兩序列絕對差,分辨系數(shù)ρ取0.5,計(jì)算公式為:

(8)

3)步驟3:結(jié)合上述計(jì)算數(shù)據(jù)和權(quán)重計(jì)算結(jié)果ωm,采用下式計(jì)算水資源短缺風(fēng)險(xiǎn)指數(shù)Rn,并判定相應(yīng)的風(fēng)險(xiǎn)等級,即:

(9)

1.3 模糊綜合法

模糊綜合法是一種利用模糊數(shù)學(xué)隸屬度原理評價(jià)受多因素制約對象的常用方法,可以將定性分析轉(zhuǎn)化成定量計(jì)算[11-12],其計(jì)算流程如下:

1)步驟1:采用函數(shù)式R={u1,u2,…,um}表示水資源短缺風(fēng)險(xiǎn),其中u1、u2、…、um代表參評因子。考慮區(qū)域水資源特點(diǎn)和開發(fā)利用現(xiàn)狀,將短缺風(fēng)險(xiǎn)劃分成低(可忽略)、較低(可接受)、臨界、較高(較嚴(yán)重)和高(無法承受)5個(gè)等級,并建立等級集R′={1,2,3,4,5}。

2)步驟2:采用公式(2)～(3)計(jì)算各參評因子隸屬于不同風(fēng)險(xiǎn)等級的程度rm(j-1),根據(jù)隸屬度函數(shù)(1)確定極大型和極小性因子隸屬度函數(shù),即:

(10)

(11)

(12)

式中:m、n為參評因子和對象;umn為參評對象n的第m個(gè)因子值;j為風(fēng)險(xiǎn)等級j;βmj、βm(j-1)為參評因子m隸屬于風(fēng)險(xiǎn)等級j和(j-1)的臨界值。根據(jù)各因子隸屬度構(gòu)造模糊判斷矩陣B,然后將模糊矩陣B與權(quán)向量W進(jìn)行乘積運(yùn)算獲取評價(jià)水資源承載力風(fēng)險(xiǎn)評價(jià)向量,并按照最大隸屬度原則判定風(fēng)險(xiǎn)等級,具體如下:

(13)

(14)

式中:ω1,ω2,…,ωM為參評因子權(quán)重。

2 水資源短缺風(fēng)險(xiǎn)評價(jià)

2.1 組合權(quán)重值

充分考慮沈陽市水資源特點(diǎn)和開發(fā)現(xiàn)狀,從工程技術(shù)、社會(huì)和自然3個(gè)方面選擇降雨量、污水處理率、萬元GDP耗水量、環(huán)境用水量、農(nóng)業(yè)用水量等8個(gè)因子評價(jià)水資源短缺風(fēng)險(xiǎn),各因子數(shù)據(jù)來源于沈陽市水資源公報(bào)、統(tǒng)計(jì)年鑒和國民經(jīng)濟(jì)統(tǒng)計(jì)資料等[13-15]。將各因子劃分成低、較低、臨界、較高和高5個(gè)風(fēng)險(xiǎn)等級,然后利用G1法和熵權(quán)法計(jì)算參評因子權(quán)重如表1,綜合利用灰關(guān)聯(lián)法和模糊綜合法評價(jià)2008—2021年沈陽市水資源短缺風(fēng)險(xiǎn)。

表1 水資源短缺風(fēng)險(xiǎn)因子權(quán)重

2.2 模糊綜合評價(jià)分析

參考序列為需水量與供水量之比,采用模糊綜合法評價(jià)沈陽市2008—2021年水資源短缺風(fēng)險(xiǎn),并按照最大隸屬度原則確定風(fēng)險(xiǎn)等級及其變化趨勢如圖1所示。結(jié)果表明,沈陽市水資源短缺風(fēng)險(xiǎn)總體表現(xiàn)出下降趨勢,2008年短缺風(fēng)險(xiǎn)達(dá)到5級(無法承受)水平,2009年和2010年短缺風(fēng)險(xiǎn)達(dá)到4級(較嚴(yán)重)水平,除2016年達(dá)到3級臨界狀態(tài)外,從2013之后整體處于1級(可忽略)和2級(可接受)水平。

圖1 模糊綜合法評價(jià)的風(fēng)險(xiǎn)等級

2.3 灰關(guān)聯(lián)分析

采用灰關(guān)聯(lián)法評價(jià)分析沈陽市2008—2021年水資源短缺風(fēng)險(xiǎn),并進(jìn)一步揭示其變化趨勢如圖2所示。結(jié)果表明,沈陽市水資源系統(tǒng)面臨著一定的短缺風(fēng)險(xiǎn),2008—2021年達(dá)到臨界風(fēng)險(xiǎn)等級以上的年份所占比例為71.4%,并以2008年的風(fēng)險(xiǎn)值最高,風(fēng)險(xiǎn)等級接近于高級(無法承受)水平;隨后沈陽市水資源短缺風(fēng)險(xiǎn)呈波動(dòng)下降趨勢,波動(dòng)下降到臨界或較高等級后逐漸穩(wěn)定,2019年之后短缺風(fēng)險(xiǎn)有所上升,于2021年又達(dá)到臨界等級以上水平。

圖2 灰關(guān)聯(lián)法評價(jià)的風(fēng)險(xiǎn)等級

2.4 結(jié)果對比

通過對比灰關(guān)聯(lián)法和模糊綜合法評價(jià)結(jié)果可知,2008—2016年期間的評價(jià)結(jié)果基本相同,兩種方法評價(jià)結(jié)果均表明沈陽市水資源短缺風(fēng)險(xiǎn)呈下降趨勢,特別是2013—2014年風(fēng)險(xiǎn)水平下降到較低等級以下,這與沈陽市水資源系統(tǒng)狀況高度相符,該期間沈陽市降水量較為充沛,外調(diào)水量也有所增加,在一定程度上緩解了城市供水壓力。2015—2019年,沈陽市水資源短缺風(fēng)險(xiǎn)基本穩(wěn)定,雖然人口規(guī)模有所增大,但人均水資源占有量和水資源利用總量也不斷增大,沈陽市水資源供給壓力逐年減小,這表明城市水資源管理逐漸發(fā)揮成效?；谊P(guān)聯(lián)計(jì)算結(jié)果顯示,2020—2021年沈陽市水資源短缺風(fēng)險(xiǎn)呈上升趨勢,從較低逐漸上升到臨界風(fēng)險(xiǎn)等級,而模糊綜合法表明該時(shí)段為低風(fēng)險(xiǎn)水平。通過分析源數(shù)據(jù)可知,2020年和2021年沈陽市生活用量明顯增加,降水量相對減小,所以灰關(guān)聯(lián)分析結(jié)果更加符合實(shí)際情況,能夠更客觀真實(shí)地反映城市水資源變化情況。

從水資源總量上,當(dāng)前沈陽市人均水資源量還低于全國平均水平,按照國際標(biāo)準(zhǔn)屬嚴(yán)重缺水城市。隨著人口的增加和經(jīng)濟(jì)的不斷發(fā)展,在未來一段時(shí)間內(nèi)沈陽市將依然存在水資源短缺風(fēng)險(xiǎn),并且外調(diào)水通道單一使得沈陽市供用水面臨著較大的風(fēng)險(xiǎn)。因此,為保障沈陽市供水安全,還要從鼓勵(lì)全民節(jié)水、減少水環(huán)境污染和提高水資源利用率等方面采取措施,以保證水資源的高效持續(xù)利用和水系統(tǒng)良性循環(huán)。

3 結(jié) 論

文章利用灰關(guān)聯(lián)法和模糊綜合法評價(jià)2008—2021年沈陽市水資源短缺風(fēng)險(xiǎn),結(jié)果表明2008—2016年期間兩者的評價(jià)結(jié)果基本相同,均表明沈陽市水資源短缺風(fēng)險(xiǎn)呈下降趨勢。但個(gè)別年份兩種方法評價(jià)結(jié)果存在一定差異,特別是2019年后灰關(guān)聯(lián)法認(rèn)為沈陽市水資源短缺風(fēng)險(xiǎn)呈上升趨勢,從較低逐漸上升到臨界風(fēng)險(xiǎn)等級,而模糊綜合法表明該時(shí)段為低風(fēng)險(xiǎn)水平,結(jié)合源數(shù)據(jù)可知灰關(guān)聯(lián)分析結(jié)果更加符合實(shí)際情況。文章綜合利用主、客觀賦權(quán)法有效避免了單一賦權(quán)法存在的不足,保證了各參評因子權(quán)重計(jì)算和評價(jià)結(jié)果的科學(xué)合理性,但不同方法對參評因子權(quán)重仍有一定影響。為更加準(zhǔn)確地揭示水資源短缺風(fēng)險(xiǎn)變化特征,未來應(yīng)進(jìn)一步研究隸屬度函數(shù)的合理構(gòu)造和風(fēng)險(xiǎn)等級的科學(xué)分類等內(nèi)容。