(江蘇省水文水資源勘測局連云港分局,江蘇連云港222004)
水資源是重要的生命資源,也是一種重要的戰(zhàn)略資源。如今,水資源短缺、水污染等問題已嚴重制約了城市的經濟發(fā)展和現(xiàn)代化進程。隨著中國水資源供需矛盾的日趨緊張和最嚴格水資源管理制度的實施,眾多學者已開始重視水資源對社會經濟的支撐能力的研究,即水資源承載力。承載力最開始是在物理學領域提出的,是指物體在不產生任何破壞時所能承受的最大負荷。后來,國外學者將這一概念引入生態(tài)學領域,水資源承載力是承載力概念在水資源領域的具體應用。在國內,最早由施雅風院士于20世紀90年代初明確提出水資源承載力的概念:水資源承載能力是指在某一歷史發(fā)展階段,以可以預見的技術、經濟和社會發(fā)展水平為依據,以可持續(xù)發(fā)展為原則,以維護生態(tài)環(huán)境良性發(fā)展為條件,經過合理的優(yōu)化配置,水資源對該區(qū)社會經濟發(fā)展的最大支撐能力[1]。水資源承載力是可持續(xù)發(fā)展與水資源安全戰(zhàn)略中的一個突出課題,也是水資源科學研究中的一個熱點問題[2]。城市的水資源承載能力受到“自然—社會”二元模式的影響,自然系統(tǒng)、社會系統(tǒng)以及水資源系統(tǒng)相互影響、相互作用。因此,對城市水資源承載力進行評價,必然要涉及到資源、環(huán)境和社會經濟發(fā)展水平等眾多因素的影響。目前,研究水資源承載力方法眾多,主要有指標體系評價法、經驗公式估算法和系統(tǒng)動力學法等。水資源承載力評價就是采用合適的水資源承載力評價指標,建立相關的數學模型,對一個地區(qū)的水資源開發(fā)利用程度進行評價,分析其水資源承載力狀況,并為水資源可持續(xù)利用提供科學依據和決策[3]。迄今為止,對水資源承載力的評價還未有一個統(tǒng)一的模型,各個方法都有各自的優(yōu)缺點。對城市水資源承載力進行研究,對于實行最嚴格水資源管理制度與“三條紅線”控制具有重要的意義。
江蘇省是個經濟大省,水資源對社會經濟的可持續(xù)發(fā)展具有重要作用,眾多學者對江蘇省的水資源承載力進行了研究。黃莉新[4]使用系統(tǒng)動力學方法計算了包含經濟增長率、供水水平和保證率3種參數8種典型方案的江蘇省水資源承載能力。宋曉猛[5]等人以承載人口規(guī)模為綜合指標按照溫飽型、過渡性、小康型以及初步富裕型4種標準對江蘇省水資源承載力進行了評價。馬學良[6]等人運用半升Γ型分布函數對城鎮(zhèn)化進程中的江蘇省水資源承載力進行了縱向分析與研究。許朗[7]等人則從時間與空間上運用主成分分析法對江蘇省水資源承載力進行了綜合評價。本文在總結已有研究成果的基礎上,綜合考慮主觀賦權和客觀賦權各自的優(yōu)缺點,選取不同的指標構建了江蘇省水資源承載力評價的數學模型,對江蘇省2004—2016年的水資源承載力進行了量化分析,為水資源承載力的科學評價提供了依據。
綜合指標體系法是對水資源承載力進行量化評價的常用方法,其基本方法是選擇單項或者多項指標,通過一定的數學模型對水資源承載力進行量化計算,從而判斷水資源承載力的相對大小。
按照水資源承載力的概念,水資源承載力是在一定條件下,區(qū)域水資源對人類社會經濟活動支撐能力的閾值,它是在一定條件約束下的目標值或者理想值[8],本文通過量化分析過程建立評價指標體系,最終計算出來的承載力值是一個介于0~1之間的無量綱值,通過這個數值反映區(qū)域水資源承載力的相對大小。
影響水資源承載力的各個單項指標與綜合承載力之間并不是簡單的線性關系,為保證綜合承載力評價結果可靠,單項指標承載力的計算需要合理進行。常規(guī)的做法是對各項指標進行歸一化和無量綱化,然后直接將指標代入模型計算。這實際上假定了單向指標與綜合承載力關系是線性的,這一假定并非對每一項指標都是合理的。在此引入單項指標影響函數的概念,影響函數不僅可以準確給出單項指標對上一層指標的影響變化趨勢,同時也完成了對單項指標的無量綱化和歸一化[9]。
在選取單項指標承載力的計算模型時,借鑒了水安全度的計算模型。對于各單項指標承載力值來說,正向指標承載力隨著指標值的增大而增大,當指標值增加到一定程度后,承載力的增加則越來越緩慢,并趨于飽和;負向指標承載力的值則隨著指標值的增加而減小,而當指標值增加到一定程度后承載力的減小越來越慢,并趨于零。即是x的減函數,常見函數中的對數函數、指數小于1的冪函數等符合這一變化規(guī)律[10]。對數函數關系簡單,適合單純指標的計算,還可以準確地反映二者之間的關系。因此,本文選用對數函數[11]作為單項指標承載力的計算模型,即:
y=a+blgx
(1)
式中y、x——因變量和自變量;a、b——模型參數。
為了確定承載力模型中的模型參數,需明確各評價指標的核算標準,根據國際公認的部分標準值、相關標準、區(qū)域實際情況以及現(xiàn)有專家的研究成果上,經過綜合考量,最終確定各指標的核算標準,并確定相應的評分,從而建立各指標的承載力模型。
綜合承載力計算公式為:
(2)
式中Ei——各單項指標所對應的承載力;wi——相應的權重。
計算出來的綜合承載力是一個處于0~1之間的無量綱值,該評價值反映水資源承載力的相對大小,越接近1,說明水資源承載力狀況越好,水資源對社會經濟的支撐能力越大;反之,越接近0,說明水資源承載力狀況越差,水資源對社會經濟的支撐能力越小。結合現(xiàn)有研究成果與有關專家意見,對水資源承載力的取值劃分為4個區(qū)間[12]:Ⅰ級為不可承載狀態(tài),判斷區(qū)間為(0~0.2);Ⅱ級為弱承載狀態(tài),判斷區(qū)間為(0.2~0.5);Ⅲ級為基本可承載狀態(tài),判斷區(qū)間為(0.5~0.8) ;Ⅳ級為良好可承載狀態(tài),判斷區(qū)間為(0.8~1)。
在進行水資源承載力評價時,需確定各評價指標的大小。關于權重的確定,一般分為主觀賦權法和客觀賦權法兩大類,主觀賦權法主要是專家根據經驗進行主觀判斷得到權重,容易受到主觀偏重的影響,包括層次分析法、專家調查法、二項系數法等;客觀賦權法是從指標的數值差異程度客觀地反映指標的重要程度,容易受到原始數值的影響,主要有熵權法、變異系數法、主成分分析法等。為了使指標權重更加合理,避免單一賦權的片面性,將主觀賦權法和客觀賦權法結合,確定組合賦權。
1.2.1主觀權重的確定
主觀權重的確定選用層次分析法(簡稱AHP)。層次分析法的原理是將與決策有關的元素分解成目標、準則、指標等層次,在同一層次對該層次元素間進行兩兩比較,建立判斷矩陣,然后對判斷矩陣進行一致性檢驗,得出指標層元素對準則層元素的權重。層次分析法的方法與步驟見圖1。
圖1 層次分析法確定權重的流程
1.2.2客觀權重的確定
客觀權重的確定方法采用變異系數法,是直接利用各項指標所包含的信息,通過計算得到指標的權重,其基本原理是:在評價指標體系中,指標取值差異越大的指標,也就是越難以實現(xiàn)的指標,這樣的指標更能反映被評價單位的差距。
由于評價指標體系中的各項指標的量綱不同,不宜直接比較其差別程度。為了消除各項評價指標的量綱不同影響,需要用各項指標的變異系數來衡量各項指標取值的差異程度。各項指標的變異系數公式如下:
(3)
各項指標的權重為:
(4)
1.2.3組合權重的確定
設用層次分析法得到的指標層權重為Wi=(Wi1,Wi2,…,Wik),由變異系數法得到的指標層權重為Wj=(Wj1,Wj2,…,Wjk),則指標層第k個指標組合權重為:
(5)
江蘇省地處長江、淮河下游,東瀕黃海,省內分布長江、太湖、淮河和沂沭泗四大水系,長江橫穿東西,大運河縱貫南北,經濟發(fā)達,城市密集,人口眾多。其氣候具有明顯的季風特征,處于亞熱帶向暖溫帶過渡地帶。全省年降水量為700~1 250 mm,年徑流深為150~400 mm,降水主要集中在6—9月汛期。
江蘇本地水資源不足,過境水豐沛,全省多年平均本地水資源量322億m3,多年平均過境水量為9 490億m3,是本地水資源量的30倍。豐富的過境水資源,特別是長江的過境水資源,為江蘇省的開發(fā)利用提供了得天獨厚的條件。然而長江、淮河匯集200萬km2國土面積的流域洪水穿境入海,縱橫交錯的流域行洪河道把全省分割成眾多區(qū)域,平原洼地普遍處于洪潮水位之下,面臨外洪內澇威脅。降水量在年際之間變化大,年內又集中在汛期,容易出現(xiàn)突發(fā)性、災害性的暴雨洪水,以及旱澇急轉、干旱缺水。特定的地理位置和氣候特點,導致江蘇省局部洪水頻繁、區(qū)域洪水常見、流域洪水易發(fā),決定了防洪任務的艱巨性;流域上中游來水主要集中在汛期,水資源時空分布不均,調蓄能力差,大部分排泄入江入海,形成棄水而不能利用,淮北、沿海、丘陵高地階段性水資源緊缺。
結合江蘇省實際情況,并根據指標體系選擇的原則,以水資源承載力為目標層,準則層選取水資源數量因素、社會經濟用水水平、水資源質量因素,選取人均水資源占有量、產水模數、折合降水總量、年入境水量、人均供水量、萬元GDP用水量、萬元工業(yè)增加值用水量、農田灌溉畝均用水量、污徑比、優(yōu)于Ⅲ類水控制河長比例、水功能區(qū)達標率、廢污水排放量等影響因素作為指標層,建立評價指標體系。為方便起見,用字母表示各層要素,見表1。
表1 水資源承載力評價指標體系
通過變異系數法確定各指標的權重以及單項指標承載力值的計算需要根據各指標的原始值。
由各指標層因素,查閱相關資料,根據《中國統(tǒng)計年鑒》《江蘇省統(tǒng)計年鑒》《江蘇省水資源公報》等相關資料,統(tǒng)計出各單項指標2004—2016年的原始數據。統(tǒng)計結果見表2。
表2 江蘇省2004—2016年水資源承載力各指標原始值
水資源單項指標承載力的計算用的是對數函數,根據《國家級生態(tài)縣、生態(tài)市、生態(tài)省建設指標》《全國水生態(tài)文明城市評價指標體系》和國際公認的標準以及相關專家的研究成果上,將各指標的核算標準分為三級并給出評分。最差值評分確定為0,及格值為0.6,最優(yōu)值確定為1。評價指標有正向指標與負向指標之分。正向指標表示指標越大承載力值越大,用“+”表示,負向指標表示指標越大承載力值越小,用“-”表示。
根據各指標的核算標準,通過對數函數確定模型參數,從而確定各單項指標承載力模型,見表3。
表3 單項指標評價標準及承載力模型
對評價指標兩兩比較,構造判斷矩陣,使用層次分析法確定相應的權重,再根據各指標的原始值利用變異系數法公式確定各指標權重,最后由式(5)得出各指標組合權重,結果見表4。
表4 水資源承載力各指標的權重
由各指標承載力的原始數據,按照確定的各單項指標承載力模型,計算歷年單項指標承載力。當正向指標小于最差值或負向指標大于最差值時,所得出的單項承載力值為負數,表示不可承載,記為零;而當正向指標大于最優(yōu)值或負向指標小于最優(yōu)值,計算出的承載力值大于1,表示完全可承載,記為1。計算結果見表5。
表5 水資源承載力各單項指標承載力值
由表5各指標單項承載力的結果和表4中的各指標的權重,將其代入式(2),得出綜合承載力,并根據水資源承載力劃分區(qū)間得出結果。
表6 水資源承載力評價結果
由圖2及表6可以發(fā)現(xiàn),2004—2016年江蘇省水資源承載力基本上呈現(xiàn)一個穩(wěn)中上升的過程,由2004年的0.304上升至2016年的0.614。2004—2014年,江蘇省水資源承載力處于Ⅱ級狀態(tài),2013年以后,承載力上升幅度變大,2014年以后上升為Ⅲ級。
分析圖2,江蘇省水資源承載力在2005—2006年、2007—2008年以及2011—2013年有小幅度下降的過程。分析各指標單項承載力,2006年水資源承載力由2005年的0.359下降至0.330主要是由水資源數量減小以及農田灌溉畝均用水量、污徑比、廢污水排放量增加引起的;2008年水資源承載力由2007年的0.373下降至0.335主要是由水資源數量減少以及農田灌溉畝均用水量增加、污徑比增加以及優(yōu)于Ⅲ類水控制河長比例減小引起的;而2011—2013年水資源承載力下降幅度最為明顯,由2011年的0.440下降至2013年的0.394,下降了10.5%,主要原因是水資源數量減小、污徑比增加以及廢污水排放量增加。2013年以后,隨著最嚴格水資源管理制度的實施,水資源承載力增加幅度明顯變大,2015、2016年已達到Ⅲ級。
圖2 2004—2016年江蘇省水資源承載力變化趨勢
本文基于組合權重結合對數函數量化分析了2004—2016年的江蘇省水資源承載力。為了使評價結果更具科學性,通過層次分析法和變異系數法相結合計算組合賦權,避免了單一賦權方法的片面性,并通過對數函數對各評價指標進行了歸一化和無量綱化,最后得出的承載力值是一個無量綱數值。計算結果表明,近年來,江蘇省水資源承載力整體上處于穩(wěn)中上升的狀態(tài),然而仍有較大的改善需求和空間。從各項評價指標來看,水資源數量以及水資源質量是水資源承載力的主要影響因素。隨著最嚴格水資源管理制度的實施和水資源“三條紅線”的逐漸落實,江蘇省水資源承載力必然會繼續(xù)增長。
許朗[7]等人通過主成分分析法得出2000—2009年江蘇省水資源承載力整體上呈現(xiàn)穩(wěn)步上升的趨勢;馬海良[6]使用半升Γ型分布函數得出2003—2013年江蘇省的水資源承載力經歷了小幅度波動增長的發(fā)展階段。雖然由于指標的選取與評價方法的不同,評價結果存在一定的誤差,但是總體趨勢基本一致,因此,本文的評價模型具有較好的可信度。為有效提高江蘇省水資源承載力,可以從以下2個方面著手。
a) 開源與節(jié)流并舉。可以開展河流的疏浚工作,因地制宜開發(fā)新水源,例如沿海地區(qū)可以考慮海水利用以及海水淡化,有條件的地區(qū)開展再生水利用。同時,要加強民眾的節(jié)水意識,大力開展節(jié)水型企業(yè)與節(jié)水型社會的建設,全面提高水資源的利用效率。
b) 嚴格落實最嚴格水資源管理制度,守住水資源“三條紅線”,改善水環(huán)境質量。加快轉變經濟發(fā)展模式,嚴格控制排污總量,加大對高能耗、高污染企業(yè)的管理,控制工業(yè)廢水的排放以及農業(yè)面源污染。