張涵俐，崔雯雯，王 越，謝祥財

(1.福建農(nóng)林大學園林學院，福建 福州 350002；2.國家水利風景區(qū)研究中心，福建 福州 350002)

在快速城市化背景下，水資源的供需量與城市社會經(jīng)濟發(fā)展之間存在著緊密復雜的協(xié)同博弈關(guān)系。隨著人類活動的需求逐步擴大，城市出現(xiàn)水源短缺、河湖聯(lián)通阻塞等一系列水網(wǎng)空間退化的現(xiàn)象，嚴重制約其可持續(xù)發(fā)展[1]。當前，水資源的需求量在不斷增加，水環(huán)境容量與水生態(tài)空間的約束日益凸顯，水資源安全保障的主要矛盾逐漸從供給不足向過度承載轉(zhuǎn)變[2]。根據(jù)聯(lián)合國在2018年巴西利亞第八屆世界水論壇上發(fā)布的報告指出，由于人口增長、經(jīng)濟發(fā)展，全球?qū)λY源的需求每年增加1%，且這一速度將會在未來大大加快[3]。在國內(nèi)，據(jù)2015年環(huán)境狀況報告顯示，在202個城市的地下水監(jiān)測點中，僅有9.1%的水質(zhì)良好，而61.3%被認為是差或更差[4]。如何應對城市社會與水資源負載能力之間的矛盾，是水環(huán)境現(xiàn)階段可持續(xù)發(fā)展亟待解決的問題。

城市水資源承載力是指地區(qū)在一定的社會發(fā)展水平下，水資源所能支持的最大工業(yè)量、人口數(shù)和城市規(guī)模度，其反映了水資源的利用程度[5-6]。近年來，隨著水資源問題的惡化，探究水資源與社會經(jīng)濟發(fā)展、生態(tài)環(huán)境保護間協(xié)調(diào)與利用程度的水資源承載力研究成為熱點。學者們選取多個子系統(tǒng)來統(tǒng)籌分析水資源承載力，對其進行不同維度的綜合分析，如Jia等[7]發(fā)現(xiàn)鄭州市水資源的驅(qū)動力子系統(tǒng)和壓力子系統(tǒng)的水資源承載力呈現(xiàn)波動性變化；何偉等[8]在子系統(tǒng)評價的基礎(chǔ)上，發(fā)現(xiàn)區(qū)域水資源承載力總體呈波浪式上升趨勢，其中生態(tài)子系統(tǒng)和社會經(jīng)濟子系統(tǒng)的貢獻力度最大；楊海燕等[9]發(fā)現(xiàn)煙臺市水資源承載力受社會、水環(huán)境及其他子系統(tǒng)的影響，在不同階段存在不同程度的波動。除上述學者外，就研究方法來看，目前有關(guān)水資源承載力評價的方法包括生態(tài)足跡法[10]、系統(tǒng)動力學[11-12]等。然而，傳統(tǒng)的指標模型并不能很好地綜合考慮水生態(tài)系統(tǒng)與多重因子間的關(guān)系。且在權(quán)重確定方面，多運用主觀的指標分析法來確定[10-12]，以此構(gòu)建的判斷矩陣一致性強，但在量化過程中容易受單一因子的臆斷影響。TOPSIS算法與樣本的數(shù)量并無較強聯(lián)系，且能充分表達數(shù)據(jù)的原有含義，同時PSR模型能有效反映水資源承載力與各子系統(tǒng)之間的關(guān)聯(lián)度，因此二者結(jié)合能作為水資源承載力量化的有效工具。本研究基于PSR模型框架下，運用組合賦權(quán)-TOPSIS法確定評估指標的加權(quán)因子，避免出現(xiàn)單個評價樣本等主觀因素的影響，使研究結(jié)果更具可靠性與科學性。

福建東瀕大海，近北回歸線，全省水資源總量優(yōu)于全國平均。但省內(nèi)水資源分布不協(xié)調(diào)，尤其是東部沿海經(jīng)濟發(fā)達地區(qū)水資源貧乏，研究發(fā)現(xiàn)晉江流域、木蘭流域開發(fā)利用率已經(jīng)接近國際公認水資源開發(fā)警戒線40%[13]，水資源嚴重短缺。目前，各學者對于長江水資源環(huán)境[14-15]的研究較為豐富，但是對福建省的相關(guān)研究較少[16]，特別是針對水資源承載力及其制約因素的研究。為此，為支持東南沿海城市水資源的空間均衡與可持續(xù)發(fā)展，本研究基于福州市水環(huán)境現(xiàn)狀來構(gòu)建評價體系，并運用綜合賦權(quán)法建立基于PSR體系的福州市水資源承載力綜合評價模型，以研究2010-2020年間水資源承載力狀況及其制約因素。

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

1.1 研究區(qū)概況

福建省的省會福州市，位于我國東部沿海地區(qū)，地處閩江中下游，客水資源豐富，為全省的政治、經(jīng)濟中心。福州市地域總面積11968km2，屬典型的亞熱帶季風氣候，氣溫適宜，年均降水量為900～2100mm，年均氣溫為16～20℃，年相對濕度約77%，7-8月為臺風高發(fā)期[17-18]。近年來，福州境內(nèi)的人均水資源占有量僅達726m3(不包括客水補給)，低于全省人均水資源量的50%。在供水方面，福州市域河流流量小、流速快，致使蓄水工程的調(diào)節(jié)能力相對較差。這種水資源與經(jīng)濟發(fā)展水平的不匹配，導致出現(xiàn)區(qū)域不透水面積迅速擴張、河網(wǎng)萎縮等現(xiàn)象[13]，使得福州市水環(huán)境問題日益突出。因此，科學研究水資源承載力，把握水資源現(xiàn)存問題是該區(qū)域亟待解決的問題。

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

福州市2010-2020年水資源承載力評價指標原始數(shù)據(jù)值來源于《福州市水資源公報》《福州市國民經(jīng)濟和社會發(fā)展統(tǒng)計公報》《福州市統(tǒng)計年鑒》《福建省統(tǒng)計年鑒》和福建省科學技術(shù)廳等。

1.3 水資源承載力評價體系構(gòu)建

綜合考慮福州市當?shù)厮Y源及社會經(jīng)濟現(xiàn)狀，本研究采用有關(guān)環(huán)境診斷的“壓力-狀態(tài)-響應”(PSR)模型來構(gòu)建水資源承載力評價指標體系[19-20]。

在借鑒相關(guān)水資源承載指標分級研究成果[7，10，21]的基礎(chǔ)上，對福州市實際情況進行綜合分析，共選取17個指標(見表1)。

表1 福州市水資源承載力評價指標體系

經(jīng)濟的發(fā)展不可避免地帶動人口增長，因此本研究采用與水資源相關(guān)的人類活動來評價壓力，分別為城鎮(zhèn)化率、人口密度、人均GDP以及萬元GDP用水量，同時選取生態(tài)環(huán)境用水量、工業(yè)用水量以及農(nóng)田灌溉畝均用水量來表征人類社會活動對水資源承載負荷所造成的影響；在選取狀態(tài)類指標時考慮到水環(huán)境當前呈現(xiàn)出的現(xiàn)狀水平，包括人均水資源量、產(chǎn)水模數(shù)、年降水量以及人均供水量，再者考慮到排水情況也能反映水環(huán)境狀況，故選取廢水排放量因子；響應類指標表征人類為推動水資源可持續(xù)發(fā)展所采取的措施，選取了R&D經(jīng)費支出占GDP比重和第一產(chǎn)業(yè)支出占GDP比重，同時為反映生態(tài)資源對福州市水資源承載力的影響，選取城市污水處理率、供水能力以及森林覆蓋率3項指標。

1.4 研究方法

1.4.1組合賦權(quán)法

層次分析法等主觀賦權(quán)法缺乏客觀性和穩(wěn)定性，熵權(quán)法等客觀賦權(quán)法易受樣本隨機性影響。因此本文采用組合賦權(quán)法來確定指標權(quán)重，在保留經(jīng)驗打分的同時保證精度[22-23]。結(jié)合AHP法與熵權(quán)法，建立組合權(quán)重計算公式：

(1)

1.4.2TOPSIS評價模型

TOPSIS(Technique for Order Preference by Similarity to an Ideal Solution)即優(yōu)劣解距離法，是一種基于評價對象與最優(yōu)解、最劣解之間相對距離的排序方法，通過不斷逼近理想解的方式進行評價[24]。

(1)構(gòu)建矩陣

Y=[yij]n·m

(2)

(2)確定理想解

(3)

(4)

(3)計算指標值與正負理想解之間的距離Dj+與Dj-以及貼近度

(5)

式中，Cj-歷年評價數(shù)據(jù)與最優(yōu)方案的貼近程度；貼近度Cj是介于[0，1]的數(shù)值。

(4)排序

根據(jù)貼近度結(jié)果Cj將水資源承載力水平進行分級。借鑒國內(nèi)外已有相關(guān)研究成果[24-25]，本研究將水資源承載力水平分為5個等級(見表2)。

表2 水資源承載力分級標準表

1.4.3障礙度模型

障礙度模型是以指標偏離度為基礎(chǔ)建立的數(shù)學診斷模型，通過探索阻礙事物發(fā)展的影響因素來診斷分析各項指標[23]。因此，引入因子貢獻度Fj、指標偏離度Ij以及障礙度Nj三個變量，其計算公式為：

(6)

式中，Nj-第j個指標對水資源承載力的障礙度。

2 結(jié)果與分析

2.1 水資源承載力原始數(shù)據(jù)及分析

根據(jù)查閱歷年數(shù)據(jù)[17-18]得到福州市水資源利用及評價體系數(shù)據(jù)各指標值(見表3)。在水資源利用方面，福州市的供水主要以地表水為主，地表供水比例穩(wěn)定在95%以上。在福州城鎮(zhèn)化進程中，用水結(jié)構(gòu)發(fā)生部分變化，工業(yè)用水常年穩(wěn)定在10億～15億m3上下，居民生活用水則穩(wěn)定在4.5億m3上下。但在2019年工業(yè)用水量降至最低值3.36億m3，居民生活用水降至3.9億m3。人均水資源量呈先上升后下降的趨勢，于2016年達到高峰值2435m3，2020年到達最低值726m3。

表3 福州水資源承載力各指標值

據(jù)表3可知，福州市降水量受氣候影響較大，在11年間，年徑流量屬偏枯水年至豐水年不等。城鎮(zhèn)化率、人均GDP呈逐年快速增長狀態(tài)，2020年人均GDP是2010年的2.7倍左右，持續(xù)領(lǐng)先全國平均。萬元GDP用水量從2010年的108m3下降到2020年的35m3，第一產(chǎn)業(yè)占GDP比重、農(nóng)田灌溉畝均用水量逐年下降，工業(yè)用水量趨于平穩(wěn)，但工業(yè)廢水排水量逐年上升。隨著城鎮(zhèn)化進程的加快，城鎮(zhèn)污水處理率不斷上漲，近十年間存在80%以上年份的城市污水處理率高于90%。R&D經(jīng)費支出占GDP比重有少許上升，提高了區(qū)域水環(huán)境安全指數(shù)。與此同時，森林覆蓋率呈上升趨勢，意味著福州市的綠地水土保持凈化能力逐步增強。

2.2 水資源承載力綜合評價分析

根據(jù)式(1)得到水資源承載力各指標組合權(quán)重(見表4)。其中，壓力類指標所占權(quán)重值最高，達0.3774，表明城鎮(zhèn)化及重工業(yè)的新興發(fā)展對水資源承載力有著重要影響[26]。狀態(tài)類指標所占權(quán)重最低，僅為0.2781，其內(nèi)部最高權(quán)重指標為年降水量，達0.0838。

表4 水資源承載力指標權(quán)重

根據(jù)綜合權(quán)重與TOPSIS模型，計算得到2010-2020年福州市水資源承載力綜合系統(tǒng)評價值(如圖1所示)。從時間維度上來看，2010-2020年福州市水資源承載能力水平呈先上升后下降再上升的趨勢，平均值為0.4473，整體波動較小，等級水平在Ⅳ-Ⅲ(短缺-合理)之間變動(如圖1所示)。在11a間，最低值于2013年出現(xiàn)，僅有0.3443，分析可知該年份人均水資源量、年降水量以及R&D經(jīng)費支出占GDP比重等正向指標較低，分別為1240m3·人-1、1393.9mm及0.3%(見表3)，從而影響福州市整體水資源承載力。11年間的最高值在2016年出現(xiàn)，該年份水資源條件優(yōu)渥、降水豐沛，人均水資源與降水量達2435m3·人-1與2316.5mm(見表3)。與此同時，2016年國家逐步加強對生態(tài)環(huán)境保護的重視度，水利廳財政廳等五廳局的《關(guān)于共同推進萬里生態(tài)水系建設(shè)的通知》[27]一系列環(huán)境保護政策的出臺，使得未來幾年福州市水資源承載力出現(xiàn)上升趨勢。

圖1 福州市水資源承載力綜合系統(tǒng)評價結(jié)果

2.3 水資源承載力各子系統(tǒng)評價分析

運用相同模型，對福州市水資源承載力子系統(tǒng)進行分析，計算得到2010-2020年各子系統(tǒng)水資源承載力評價值(如圖2所示)。此外，各子系統(tǒng)對水資源承載力綜合可持續(xù)利用水平的貢獻逐年增加(如圖3所示)。在壓力子系統(tǒng)中，2010-2019年較為穩(wěn)定，保持在0.40～0.52間，但2020年出現(xiàn)明顯下降。分析其原因是由于隨著城鎮(zhèn)化的加快，工業(yè)用水量、河道外生態(tài)用水量以及農(nóng)田灌溉用水量明顯增加，分別達到14.35億、1.49億、632m3(見表3)，由此導致水資源承載壓力也逐步上漲；在狀態(tài)子系統(tǒng)中，2010-2020年波動幅度較大，整體趨勢不平穩(wěn)，于2016年到達最高值0.8746，2020年到達最低值0.1914。2020年人均水資源量歷年最低，僅有726m3·人-1(見表3)，導致該年份水資源總量相對減少，供需矛盾突出；在響應子系統(tǒng)中，2010-2020年出現(xiàn)顯著上升趨勢，尤其是2018-2019年段，由0.6880增加至0.8592，增長率達24.8%。該年段內(nèi)，城市污水處理率增長至97%(見表3)，為歷年來的最大值，有效提高了福州市水資源響應子系統(tǒng)承載力。

圖2 福州市各子系統(tǒng)水資源承載力評價結(jié)果

圖3 福州市水資源承載力各子系統(tǒng)的貢獻率

2.4 水資源承載力障礙度因素診斷

根據(jù)式(6)計算出水資源承載力各指標障礙度，選擇比重居前5位的指標作為主要障礙因子并對其進行排序(如圖4所示)。2010-2016年，比重居于前二的為A17森林覆蓋率、A16供水能力；2017-2020年，森林覆蓋率不再是主要障礙因素，A9年降水量與A6工業(yè)用水量輪流領(lǐng)先，且A10產(chǎn)水模數(shù)、A8人均水資源量障礙度進一步提升。

圖4 福州市水資源承載力主要障礙因子(前五名)

按照障礙度Nj≥5.0%的標準進行篩選，將頻率從高到低排序，以確定障礙因子主次關(guān)系(見表5)。

表5 福州市水資源承載力主要障礙因子頻率表

近11年來，供水能力和年降水量是水資源承載力障礙度中出現(xiàn)頻率最高的兩個因子，出現(xiàn)頻率同為90.90%，其次是產(chǎn)水模數(shù)與森林覆蓋率，出現(xiàn)頻率同為81.81%。在出現(xiàn)頻次較高的因子中，供水能力與產(chǎn)水模數(shù)均為水資源供需量方面的因子。

針對表5中發(fā)生頻率大于50%的8項指標，繪制2010-2020年福州水資源承載力主要障礙因子熱點圖(如圖5所示)。A8人均水資源量、A9年降水量、A10產(chǎn)水模數(shù)的障礙度有逐年增大的趨勢，表明這3項指標對福州市的水資源承載力具有深遠影響，其中應注重年降水量，在2016年后增幅較大，最大增幅達23.9%，可能會成為福州市水資源承載力水平的限制因素。

圖5 福州市水資源承載力主要障礙因子熱點圖

3 討論

(1)人均水資源量與年降水量是阻礙水資源承載力狀態(tài)子系統(tǒng)的關(guān)鍵。人均水資源量可以體現(xiàn)區(qū)域自然地理條件與水資源的協(xié)調(diào)程度，年降水量則影響著地表水資源的豐富程度，對水文循環(huán)有直接的影響。程度良等人[28]發(fā)現(xiàn)影響鄂爾多斯市水資源承載力提升的主要障礙因素是人均水資源量和產(chǎn)水模數(shù)，水資源子系統(tǒng)則是阻礙鄂爾多斯市水資源承載力提升的首要子系統(tǒng)；徐政華等人[29]得出提升長春市水資源承載力的關(guān)鍵在于提升單位面積水資源量、人均水資源利用量等指標。王金星等人[30]研究1999-2020年甘肅省WRCC得到人均水資源量、人均用水量、年均降水量、產(chǎn)水模數(shù)和水功能區(qū)水質(zhì)達標率為該區(qū)域水資源承載力的五大障礙因子；金菊良等人[31]發(fā)現(xiàn)黑龍江省的水資源承載力總體狀況較好，而“降水量”、“第三產(chǎn)業(yè)占比”處于反勢，對于這些處于反勢的脆弱性指標應予以重點關(guān)注、調(diào)控。本研究中，在權(quán)重方面，人均水資源量與年降水量是狀態(tài)類指標權(quán)重占比最大的兩個因子，分別達到了0.0675、0.0838(見表4)；在障礙度方面，自2016年起人均水資源量與年降水量對水資源承載力的阻礙作用明顯增強，且11年間的出現(xiàn)頻率均高于70%(見表5)。在2020年，二者對福州市水資源承載力的阻礙度達到了43.16%，是障礙度中影響力占比最大的兩個因子(如圖4所示)因此，阻礙水資源承載力狀態(tài)子系統(tǒng)的關(guān)鍵就是人均水資源量與年降水量。

(2)響應子系統(tǒng)呈現(xiàn)持續(xù)增長趨勢。響應子系統(tǒng)表征著人類為促進水資源可持續(xù)發(fā)展所采取的對策與措施，反映人們積極響應策略的效應，為外部干預作用。葉正偉等人[32]在研究鹽城市水環(huán)境安全中發(fā)現(xiàn)，響應指數(shù)出現(xiàn)顯著的升高趨勢，在3個子系統(tǒng)中響應子系統(tǒng)增幅最大；Kaize Zhang等人[33]發(fā)現(xiàn)北京市水資源安全領(lǐng)域的整體績效水平呈波動增長態(tài)勢，同時壓力子系統(tǒng)與響應子系統(tǒng)的績效水平也不斷上升；Xiaoqiong Liu等人[34]發(fā)現(xiàn)響應子系統(tǒng)的年均增長率高達4.06%，表明社會各方面采取的響應措施在提高水環(huán)境治理效率方面發(fā)揮了重要作用。在本研究中，從評價結(jié)果來看，響應子系統(tǒng)保持持續(xù)上升趨勢，自2010年的0.2419提高到2020年的0.9809，增長率高達305%(如圖2所示)；從子系統(tǒng)貢獻率來看，2010-2020年響應子系統(tǒng)占總體比值由23%提升到67%(如圖3所示)。響應子系統(tǒng)是為促進社會-經(jīng)濟-生態(tài)環(huán)境協(xié)調(diào)所人為采取的各項措施行為，為主觀可控因素。因此，在意識到水環(huán)境需求量逐步擴張的現(xiàn)實后，政府采取一系列手段來進行調(diào)控，使得現(xiàn)階段響應子系統(tǒng)持續(xù)保持上漲態(tài)勢。

(3)提升水資源狀態(tài)子系統(tǒng)是提升水資源承載力的重點。水資源承載力狀態(tài)子系統(tǒng)是生態(tài)環(huán)境的現(xiàn)實狀態(tài)，與水資源子系統(tǒng)相互關(guān)聯(lián)，是指水資源子系統(tǒng)在壓力下顯示的狀態(tài)，同時是水環(huán)境治理評價目標的直接體現(xiàn)[34]。Xiaoqiong Liu等[34]發(fā)現(xiàn)狀態(tài)子系統(tǒng)是水環(huán)境治理效率的一個直觀表現(xiàn)；楊海燕等[9]在研究中提及，由于水資源子系統(tǒng)中水資源總量的不確定性，煙臺市水資源承載力易受洪澇或旱情影響發(fā)生短期內(nèi)的突升與突降；賈瓊等人[35]表明水資源子系統(tǒng)的承載力水平呈現(xiàn)出明顯的滯后性特征，并且隨著時間的推移呈現(xiàn)出波動下降趨勢。關(guān)中地區(qū)水資源量受降雨影響顯著，可能與降雨量變化有關(guān)。在本研究中，壓力子系統(tǒng)維持基本平穩(wěn)狀態(tài)，水資源承載力值在0.2882～0.5271間變動，響應子系統(tǒng)則保持持續(xù)上升的趨勢，由0.2419提升到0.9809(如圖2所示)，且二子系統(tǒng)各指標為人為調(diào)控因素居多，趨向良好穩(wěn)定的發(fā)展態(tài)勢；而對于狀態(tài)子系統(tǒng)，上下波動起伏較大，在2016年達到最高值0.8746，在2020年達到最低值0.1914，二者間相差0.6832(如圖2所示)。2016年福建省臺風接踵而至，致使全年四季降水偏多[13]，水資源狀態(tài)相對較優(yōu)，整體承載力評價達到最高。再者，狀態(tài)子系統(tǒng)表征研究區(qū)的自然稟賦條件，與洪澇或旱情等特殊天氣緊密聯(lián)系。至此，狀態(tài)子系統(tǒng)的提升對福州市水資源承載力的整體提升就尤為重要，在發(fā)展經(jīng)濟的同時，需在現(xiàn)有水資源量和降水量的基礎(chǔ)上維持對水資源的有效利用，對水資源總量的主觀調(diào)控則是提升水資源承載力的重點所在。

4 結(jié)論

(1)運用TOPSIS模型評價2010-2020年福州市水資源承載力的分析結(jié)果表明，福州市水資源承載能力水平整體較為合理，但其上升趨勢相對緩慢，未來有較大的提升空間；子系統(tǒng)方面，壓力子系統(tǒng)基本維穩(wěn)，響應子系統(tǒng)呈現(xiàn)持續(xù)上升趨勢，而狀態(tài)子系統(tǒng)則表現(xiàn)出較大波動，分析知人均水資源量與年降水量對其有著顯著影響。

(2)通過障礙度模型對福州市水資源承載力的阻礙分析結(jié)果表明，年降水量、工業(yè)用水量以及人均水資源量等因子已代替森林覆蓋率與供水能力成為影響福州市水資源承載力的主要障礙因素。就現(xiàn)階段整體趨勢來看，福州市水資源承載力的主要障礙因子集中在狀態(tài)子系統(tǒng)中，且障礙度日漸增強。

(3)為實現(xiàn)福州水系的生態(tài)保護和高質(zhì)量發(fā)展，建議應立足于現(xiàn)有水環(huán)境條件，認識到在人水關(guān)系中人需要主動適應水環(huán)境。地方政府應加大生態(tài)環(huán)保建設(shè)，引進客水資源進行生態(tài)補水，同時在重點行業(yè)實施用水總量和強度雙控制度，統(tǒng)籌協(xié)調(diào)福州水系水域管理，發(fā)揮宏觀調(diào)控作用；城鎮(zhèn)居民應適當減少用水消耗，因地制宜節(jié)約用水，從而有效提高水資源承載能力。

(4)本研究的研究區(qū)位于福建省省會福州市，其地域特色可借鑒到東南沿海城市區(qū)域，研究思路與方法可適用于其他區(qū)域乃至全國。本研究由于福州市各縣域部分指標不足，僅評估了福州市時間序列上的水資源承載力，尚未對福州市空間序列進行詳細探討，在今后的研究中可以具體到區(qū)縣的時空分析，使得評價結(jié)果更加完善。