基于系統(tǒng)動力學(xué)的河北省水資源供需平衡預(yù)測

摘要：針對河北省長期存在水資源供需矛盾的問題，采用系統(tǒng)動力學(xué)法建立河北省水資源供需模型，依據(jù)研究區(qū)的經(jīng)濟發(fā)展水平、水資源利用狀況等因素，設(shè)計5種不同情景（現(xiàn)狀延續(xù)型情景W1、經(jīng)濟發(fā)展型情景W2、節(jié)約用水型情景W3、污染治理型情景W4、協(xié)調(diào)發(fā)展型情景W5），預(yù)測河北省2022—2050年的需水量。結(jié)果表明：①5種情景下總需水量均呈遞增趨勢，平均值分別為239.47億、258.51億、215.53億、239.47億、232.95億m3；②5種情景下均處于缺水狀態(tài)，平均缺水程度分別為0.243 8、0.302 6、0.137 5、0.146 7、0.129 0；③5種情景下農(nóng)業(yè)需水量均值占總需水量均值比重最大，分別為46.57%、43.14%、45.72%、46.57%和43.85%。因此，提高農(nóng)業(yè)用水效率對減少水資源消耗起到舉足輕重的作用。

關(guān)鍵詞：系統(tǒng)動力學(xué)；水資源供需平衡；缺水程度；情景分析；河北省

中圖分類號：TV213文獻標(biāo)識碼：A文章編號：1001-9235（2024）08-0092-12

Supply-demand Balance Prediction of Water Resources in Hebei ProvinceBased on System Dynamics

WANG Jianquan1，QIN Huanhuan1，2*，HUANG Lixiang1，ZHU Shiyi1

（1.School of Water Resources and Environmental Engineering，East China University of Technology，Nanchang 330013，China;2.StateKey Laboratory of Nuclear Resources and Environment，East China University of Technology，Nanchang 330013，China）

Abstract：To ease the long-term contradiction between the supply and demand of water resources in Hebei Province，the system dynamics method was used to construct a supply-demand model of water resources in Hebei Province.Based on factors in the studied area including economic development level and water resources utilization，five different development scenarios were designed（status quo continuation scenario W1，economic development scenario W2，water conservation scenario W3，pollution control scenario W4，and coordinated development scenario W5）.Water demand in Hebei Province from 2022 to 2050 was thus predicted.According to the results，①the total water demand shows an increasing trend under the five scenarios，with average values of 23.947，25.851，21.553，23.947 and 23.295 billion m3 respectively;②all five scenarios are in a state of water shortage，with an average shortage degree of 0.243 8，0.302 6，0.137 5，0.146 7 and 0.129 0 respectively;③the average value of agricultural water demand，under the five scenarios，accounts for the largest proportion of the average value of total water demand，with 46.57%，43.14%，45.72%，46.57%，and 43.85%respectively.Therefore，agricultural water use efficiency improvement plays an important role in reducing water consumption.

Keywords：systemdynamics;supply-demand balance of water resources;water shortage degree;scenarioanalysis;Hebei Province

經(jīng)濟發(fā)展和人口增長導(dǎo)致水資源短缺被認(rèn)為是人類社會最重要的威脅之一，也是社會可持續(xù)發(fā)展的制約因素［1］，尤其在人口經(jīng)濟集中的城市，城市空間的擴張和經(jīng)濟規(guī)模的增長均在持續(xù)驅(qū)動水資源需求的增加［2］。中國的水資源總量長期低于世界平均水平，人均水資源量遠(yuǎn)低于世界平均水平［3］。隨著城市化的發(fā)展，中國對水資源的需求不斷增加；由于人類活動不斷進行，全球范圍內(nèi)產(chǎn)生的廢水超過80%未經(jīng)處理就排入水環(huán)境當(dāng)中［4］，導(dǎo)致水資源受到不同程度的污染；許多地區(qū)的降雨呈現(xiàn)不穩(wěn)定變化，持續(xù)干旱天氣的范圍大、頻率高，眾多因素加劇了水資源的短缺［5］。河北省作為“首都圈”的省份之一，是國家區(qū)域發(fā)展規(guī)劃的重要一環(huán)，具有十分重要的樞紐作用，但全省多年人均水資源量占有量僅307 m3［6］，是全國平均值的1/7［7］，遠(yuǎn)低于國際公認(rèn)的人均500 m3的“極度缺水標(biāo)準(zhǔn)”，是典型的資源型缺水省份［8-9］。隨著城市化和經(jīng)濟不斷發(fā)展，水資源供需矛盾將更加突出，嚴(yán)重制約當(dāng)?shù)厣鐣?jīng)濟和水資源的協(xié)調(diào)發(fā)展。因此，科學(xué)預(yù)測河北省未來水資源供需平衡，對水資源管理部門制定相關(guān)的水資源管理政策具有重要的指示意義。

需水量預(yù)測是水資源規(guī)劃與管理的重要基礎(chǔ)和依據(jù)，準(zhǔn)確地預(yù)測區(qū)域需水量對于水資源可持續(xù)發(fā)展起到重要的指導(dǎo)作用。目前，常用的需水量預(yù)測方法有指標(biāo)法［10］、定額法［11］、時間序列法［12］、灰色預(yù)測法［13］、BP-神經(jīng)網(wǎng)絡(luò)法［14］、多元線性回歸法［15］、主成分分析法［16］、系統(tǒng)動力學(xué)法（System Dynamics，SD）［17］。影響需水量的因素諸多，通常涉及到水文、社會、工程、經(jīng)濟和科技等方面，并且因素間的關(guān)系較為復(fù)雜，傳統(tǒng)的預(yù)測方法不能系統(tǒng)、準(zhǔn)確地刻畫因素間的關(guān)系。然后SD法可以簡明地反映各因素間的本質(zhì)聯(lián)系，真實有效地組織和揭示復(fù)雜水資源系統(tǒng)內(nèi)各因素之間和因素內(nèi)部的相互反饋機制［18］。相比之下，SD法能定量考慮影響需水量的各種水文、社會經(jīng)濟、科技等諸多因素，進而較準(zhǔn)確地預(yù)測需水量。

系統(tǒng)動力學(xué)是美國的福瑞斯特教授于1958年為分析生產(chǎn)管理及庫存管理等企業(yè)問題而提出的系統(tǒng)仿真方法，被應(yīng)用于解決各種復(fù)雜系統(tǒng)性問題，尤其在處理高度非線性、高階次、多變量、多重反饋問題方面具有顯著優(yōu)勢［11，17-26，29-30］。系統(tǒng)動力學(xué)是以定性與定量相結(jié)合來研究如何有效解決復(fù)雜系統(tǒng)的方法。當(dāng)前，系統(tǒng)動力學(xué)在很多領(lǐng)域得到廣泛的應(yīng)用，如需水量預(yù)測、水資源規(guī)劃與管理、環(huán)境承載力、裝配式建筑、交通運輸管理等方面［11，19-25］。

河北省缺水問題引起了專家和學(xué)者的關(guān)注，因此做了一些研究工作。李騰等［27］采用水足跡法對河北省2004—2018年的水資源可持續(xù)性進行了研究，結(jié)果顯示水資源處于虧損狀態(tài)；段一明等［28］采用指標(biāo)體系法對河北省水資源脆弱性進行分析，結(jié)果顯示水資源處于強、中等脆弱性狀況。然而針對河北省缺水問題而展開的工作大多數(shù)是對現(xiàn)狀的研究。因此，河北省未來需水量亟待預(yù)測。影響需水量預(yù)測的因素諸多，并且各因素的關(guān)系較為復(fù)雜。SD法建立的模型可以捕抓系統(tǒng)內(nèi)每一時刻各因素的狀態(tài)，進而掌握其供需之間的關(guān)系。本文采用系統(tǒng)動力學(xué)法建立河北省水資源供需SD模型，對多因素耦合的水資源復(fù)雜系統(tǒng)進行需水量預(yù)測，為河北省水資源可持續(xù)發(fā)展提供科學(xué)依據(jù)。

1材料與研究方法

1.1研究區(qū)概況

河北?。?6°05′～42°40′N，113°27′～119°50′E）是位于中國華北平原的省級行政區(qū)，省會石家莊，東臨渤海，內(nèi)環(huán)京津，西為太行山地，北為燕山山地，是中國唯一兼有高原、山地、丘陵、平原、湖泊和海濱的省份，總面積為18.88萬km2（平原區(qū)81 459 km2、山地90 280 km2、高原15 954 km2），屬溫帶大陸性季風(fēng)氣候。從圖1可知，2006—2021年的平均降水量、水資源量、總用水量分別為521.59 mm、164.01億m3和190.17億m3，多年水資源量小于總用水量（2012、2016、2021年除外），絕大多數(shù)年份處于缺水狀態(tài)；2021年末，河北省常住總?cè)丝? 448萬人，其中，城鎮(zhèn)常住人口4 554萬人，農(nóng)村常住人口2 894萬人；社會生產(chǎn)總值穩(wěn)步增長，到2021年末，生產(chǎn)總值為40 391.2億元，其中，第一、二、三產(chǎn)業(yè)生產(chǎn)值分別為4 030.3億、16 364.2億、19 996.7億元。

1.2研究方法簡介

系統(tǒng)動力學(xué)由麻省理工學(xué)院的J.W.Forrester教授提出，是一種以定性分析為先導(dǎo)、定量分析為支持的系統(tǒng)仿真方法，其以反饋控制理論為基礎(chǔ)，借助數(shù)學(xué)計算機仿真技術(shù)，研究復(fù)雜系統(tǒng)行為。當(dāng)前，系統(tǒng)動力學(xué)在國內(nèi)的研究已十分成熟，主要被應(yīng)用于解決各種復(fù)雜的系統(tǒng)性問題，特別在需水量預(yù)測方面取得一系列的科研成果［11，17-21，29-30］。

通常，利用SD法建立水資源供需模型，可以劃分為生活、工業(yè)、農(nóng)業(yè)、環(huán)境以及水資源5個子系統(tǒng)，從而達到簡化模型的效果，在此基礎(chǔ)上再對5個子系統(tǒng)作進一步的描述，以使模型逐漸充實。通常，變量在系統(tǒng)動力學(xué)中可以用狀態(tài)、流率、輔助變量以及常量表達，應(yīng)具體考慮變量的含義而作進一步確定；明確變量之間的關(guān)系并以關(guān)系方程式進行定量化，搜集模型需要的所有數(shù)據(jù)，輸入模型后運行，將模擬數(shù)據(jù)與歷史數(shù)據(jù)進行擬合后獲得模型參數(shù)后，方可利用模型預(yù)測需水量。

SD模型的實質(zhì)是使用微分方程組對各個系統(tǒng)及其變量間的關(guān)系進行定量描述［29］。假設(shè)系統(tǒng)內(nèi)每個要素變量si（i=1，2，3，…n）在其定義域內(nèi)均為可微的函數(shù)，并且滿足式（1），則這n個要素變量構(gòu)成一個系統(tǒng)。若要考慮到依賴于時間和空間，則聯(lián)立方程組就須用偏微分方程表示。式（1）表示系統(tǒng)動力學(xué)的基本原理。

2河北省系統(tǒng)動力學(xué)模型構(gòu)建

2.1系統(tǒng)動力學(xué)流程圖

河北省水資源供需系統(tǒng)由生活需水量、工業(yè)需水量、環(huán)境需水量、農(nóng)業(yè)需水量和總供水量子系統(tǒng)組成（圖2）。本文利用VENSIM Professional軟件并以缺水程度為核心建立河北省水資源供需SD模型，將與水資源有直接影響的變量考慮進系統(tǒng)邊界內(nèi)，用歐文法數(shù)值積分形式定量刻畫變量之間的關(guān)系。在流圖中，狀態(tài)變量（表示累積效應(yīng)的變量）用來表示、流率變量（表示積累效應(yīng)變化快慢的變化率變量）用來表示、輔助變量（從積累效用變量到變化率變量及變化率之間的中間變量）、常量（某一時間區(qū)間內(nèi)不隨時間變化而變化的量）、表函數(shù)（用來表示無法用顯性方程來表示變量間關(guān)系的一種方程）用來表示，表示源匯項，→表示變量之間的因果關(guān)系。

2.2模型變量、方程與數(shù)據(jù)來源

河北省水資源供需SD模型包括7對狀態(tài)變量方程、大量的速率方程、輔助變量、常量以及表函數(shù)。通過這些方程和變量，把各變量間的邏輯關(guān)系“翻譯”成為數(shù)學(xué)語言，以顯示其定量關(guān)系。由于文章篇幅有限，本文列出部分主要方程，見表1。

模型中所用數(shù)據(jù)主要包括水資源數(shù)據(jù)、社會經(jīng)濟數(shù)據(jù)和環(huán)境污染數(shù)據(jù)。其中，地下水供水量、地表水供水量、各類用水定額來自2006—2021年河北省水資源公報；人口、區(qū)域生產(chǎn)總值、牲畜數(shù)量和有效灌溉面積來源于2007—2022年河北省統(tǒng)計年鑒；污水排放系數(shù)、污水處理率和污水回用率參考2007—2022年河北省統(tǒng)計年鑒、《第一次全國污染源普查城鎮(zhèn)生活源產(chǎn)排污系數(shù)手冊》。

模擬時間為2006—2050年，基準(zhǔn)年為2006年，時間步長為1 a。其中，2006—2021年為校準(zhǔn)年份（借助表函數(shù)將增長率、用水定額、地下和地表供水量、廢污水排放系數(shù)、污水處理率和回用率等均作為模型的輸入?yún)?shù)，輸入狀態(tài)變量的初始值、常量和參數(shù)后運行模型，直至模擬結(jié)束后，輸出生活、工業(yè)、農(nóng)業(yè)和總需水量的模擬數(shù)據(jù)，將其與歷史數(shù)據(jù)進行誤差分析，從而獲得和歷史數(shù)據(jù)擬合度較好的模型參數(shù)，以供模型的預(yù)測階段使用）；2022—2050年為預(yù)測年份（利用前一階段獲得的參數(shù)，模擬出5種情景下各行業(yè)未來需水量、總需水量、污水回用量、水資源供需缺口和缺水程度）。

2.3模型的有效性檢驗與敏感性檢測

有效性檢驗是判斷模型具備反映客觀事實能力的重要依據(jù)。為此，選取4個指標(biāo)（生活需水量、工業(yè)需水量、農(nóng)業(yè)需水量和總需水量）對模型進行有效性檢驗。將歷史數(shù)據(jù)輸入模型中，檢查變量單位、模擬時長以及變量間的定量關(guān)系，確定無誤后運行模型，直至模擬結(jié)束后，輸出4個待檢驗指標(biāo)的模擬數(shù)據(jù)，將與歷史數(shù)據(jù)進行誤差分析。相對誤差值均在±5%以內(nèi)（表2），兩者的吻合度較好，說明河北省SD模型是有效的［29］，能夠較準(zhǔn)確地刻畫、模擬研究區(qū)的系統(tǒng)行為。

敏感性檢測可以驗證模型的可靠性與準(zhǔn)確性。一般而言，具有較好強壯性的模型表現(xiàn)在目標(biāo)變量對大部分外生變量是不敏感的。由于文章篇幅有限，因此，本研究展示模型中僅有2個常量（綠化用水定額、林地用水定額）的敏感性檢測結(jié)果。假設(shè)綠化用水定額、林地用水定額分別在0.001～0.700m3/m2（圖3a、3c）和800～1 400 m3/hm2（圖3b、3d）區(qū)間內(nèi)服從隨機分布，分別觀察生態(tài)環(huán)境需水量、水資源供需缺口、林牧業(yè)需水量的靈敏度，見圖3。秦劍［30］構(gòu)建了北京市水資源供需平衡系統(tǒng)動力學(xué)模型，對模型進行了敏感性檢測，結(jié)果顯示模型具有可靠性?；谖墨I［30］的判別法則，從圖3可知，模型對外生變量不敏感。結(jié)合有效性檢驗結(jié)果，進一步說明模型是可靠的，具有準(zhǔn)確性。

2.4情景設(shè)計

尋找適合河北省未來發(fā)展情景對解決水資源供需矛盾具有重要的實際意義。因此，依據(jù)河北省經(jīng)濟發(fā)展現(xiàn)狀和水資源利用情況，充分考慮了河北省未來社會經(jīng)濟發(fā)展過程中影響水資源可持續(xù)利用的因素，綜合考慮而設(shè)計了5種情景（現(xiàn)狀延續(xù)型情景W1、經(jīng)濟發(fā)展型情景W2、節(jié)約用水型情景W3、污染治理型情景W4和協(xié)調(diào)發(fā)展型情景W5），各情景下決策變量參數(shù)值見表3。決策變量預(yù)測值是依據(jù)河北省人口增長情況、經(jīng)濟發(fā)展水平、水資源利用情況和工農(nóng)業(yè)節(jié)水技術(shù)水平等因素，綜合評價并考慮各情景強調(diào)的側(cè)重點不同而確定的參數(shù)值。

3結(jié)果與分析

3.1需水量模擬結(jié)果

通過圖4、5和表4，分析河北省未來的工業(yè)、農(nóng)業(yè)、生活、生態(tài)環(huán)境和總需水量。

a）從工業(yè)需水量角度分析，5種情景的工業(yè)需水量均呈遞增趨勢。其中，情景W2（均值為75.69億m3）的增長速度最快，從20.81億m3升至189.30億m3，上升了8.1倍，其占總需水量比重在11.23%～46.88%，比重呈逐年遞增趨勢，出現(xiàn)此現(xiàn)象的主要原因是情景W2更強調(diào)未來工業(yè)發(fā)展，從而導(dǎo)致工業(yè)需水量增加；情景W1（均值為57.49億m3）、情景W4（均值為57.49億m3）和情景W5（均值為62.28億m3）的增長速度次之。情景W1、W4均從20.81億m3升至106.2億m3，上升了4.1倍，其占總需水量比重分別為11.23%～33.36%、11.23%～33.38%。情景W5從20.65億m3升至140.6億m3，上升了5.8倍，其占總需水量比重為11.20%～42.68%；情景W3（均值為49.25億m3）的增長速度最慢，從20.58億m3升至85.1億m3，上升了3.14倍，是情景W2的0.39倍，其占總需水量比重為11.20%～31.75%。

b）從農(nóng)業(yè)需水量角度分析，5種情景的農(nóng)業(yè)需水量均值占總需水量均值比例最大（表4），均值依次為111.51億、111.51億、98.54億、111.51億、102.15億m3，比重依次為46.57%、43.14%、45.72%、46.57%、43.85%，均超過43.0%，表明河北省在未來仍以農(nóng)業(yè)經(jīng)濟為主。農(nóng)業(yè)需水量均呈增長趨勢，其中，情景W1、W2和W4均從97.53億m3升至137.90億m3，漲幅為41.39%；情景W5從96.95億m3升至117.5億m3，漲幅為21.20%，漲幅是前三者的0.51倍；情景W3呈先降后升趨勢，總體上從96.53億m3升至112.20億m3，漲幅為16.23%，占總需水量比重依次為52.63%～43.32%、52.63%～34.15%、52.52%～41.87%、52.63%～43.33%和52.58%～35.67%，比重均呈遞減趨勢，可能與提高農(nóng)業(yè)用水效率、改變種植結(jié)構(gòu)和耕地面積減少等因素有關(guān)聯(lián)。

c）從生活需水量角度分析，相比工農(nóng)業(yè)需水量，生活需水量的漲幅處于穩(wěn)定狀態(tài)。5種情景的漲幅依次為25.6%、34.3%、14.3%、25.6%和16.1%，依次從27.75億m3漲至34.85億m3、從27.75億m3漲至37.26億m3、從27.44億m3漲至31.37億m3、從27.75億m3漲至34.85億m3、從27.56億m3漲至32.00億m3，占總需水量比重均呈遞減趨勢，依次從14.98%降至10.95%、從14.98%降至9.23%、從14.93%降至11.71%、從14.98%降至10.94%、從14.95%降至9.71%，比重逐年下降可能與推廣節(jié)水技術(shù)、提高節(jié)水意識有關(guān)。

d）從生態(tài)需水量角度分析，5種情景的生態(tài)環(huán)境需水量的變化幅度較小，均值依次為39.26億、39.27億、39.26億、39.26億、39.27億m3，其占總需水量比重呈遞減趨勢。河北省在未來需重點加強對這一方面的管理力度，確保生態(tài)環(huán)境與經(jīng)濟協(xié)調(diào)發(fā)展。

e）5種情景的總需水量均呈遞增趨勢。其中，情景W2（均值為258.51億m3）的增長速度最快，情景W1（均值為239.47億m3）、情景W4（均值為239.47億m3）、情景W5（均值為232.95億m3）的增長速度次之，情景W3（均值為215.53億m3）的增長速度最慢。3.2水資源供需平衡和污水回用量模擬結(jié)果

5種情景的水資源供需缺口均呈遞增趨勢（圖6），供需不平衡程度加劇。其中，情景W2的增長速度最快，從3.35億m3升至177.9億m3，上升了174.55億m3，表明強調(diào)經(jīng)濟發(fā)展將導(dǎo)致需水量增加，進而加劇水資源供需不平衡；情景W1的增長速度次之，從3.35億m3升至113.5億m3，上升了110.15億m3，是情景W2上升量的0.631倍；情景W3、W4和W5供需缺口增加量分別為67.32億、66.97億、80.60億m3，在模擬前期，情景W4供需缺口量均大于情景W3和W5，但在模擬后期，情景W5供需缺口增長速度變快，超過情景W3和W4。供需缺口均值依次為47.77億、62.12億、26.50億、31.86億、28.09億m3。

5種情景的污水回用量均呈遞增趨勢（圖7）。情景W2、W4和W5的增長速度比情景W1和W3快，出現(xiàn)此現(xiàn)象的原因是情景W4較W1強調(diào)污水處理率、回用率，2種情景的平均總需水量均為239.47億m3，但缺水程度分別為0.146 7、0.243 8（表5），表明提高污水處理技術(shù)可以緩解缺水壓力；情景W2較W1強調(diào)未來的經(jīng)濟發(fā)展，在污水處理率、回用率相等情況下，情景W2（26.31億m3）的平均污水回用量比W1（21.63億m3）多4.68億m3，但情景W2的缺水程度比W1大0.059，說明污水回用量對緩解缺水壓力起到較小的作用；情景W5較W1更加全面考慮了人口、社會經(jīng)濟、節(jié)水技術(shù)等因素，情景W5的污水回用量均值比W1多5.22億m3，缺水程度均值比情景W1少0.114 8，表明情景W5對緩解缺水壓力取得較好的效果。

3.3缺水程度模擬結(jié)果

缺水程度在模型中以IF THEN ELSE（水資源供需缺口大于0，水資源供需缺口/總供水量，0）的方式進行定量表達。5種情景均處于缺水狀態(tài)（圖8），且缺水程度不斷加劇。其中，情景W2（均值為0.302 6）的缺水速度最快，情景W1（均值為0.243 8）的缺水速度次之，情景W3（均值為0.137 5）、W4（均值為0.146 7）、W5（均值為0.129 0）的缺水速度較緩慢。在未來發(fā)展過程中，若河北省注重經(jīng)濟發(fā)展（情景W2），則缺水程度較情景W1將進一步加?。蝗艉颖笔〖訌娦麄鞴?jié)水力度，降低用水定額（情景W3），缺水矛盾較情景W1得到緩解，但缺水問題仍存在，需從多方面采取行動解決缺水問題；若河北省提高污水處理技術(shù)同時引進水資源量（情景W4），缺水程度較情景W1得到緩解，且在模擬期末，其缺水程度值最??；若河北省綜合考慮社會經(jīng)濟發(fā)展、污水處理技術(shù)、引水工程（情景W5），那么缺水矛盾將得到有效的緩解。

總體而言，情景W2在未來經(jīng)濟發(fā)展速度最快，但缺水問題較其他情景更為嚴(yán)重，不利于河北省水資源長期安全利用；情景W3和W4在未來的節(jié)水力度最大和污水處理效率最高，缺水矛盾較情景W1得到緩解，但大力支撐河北省經(jīng)濟發(fā)展背景下很難被采納；情景W5在未來適度強調(diào)經(jīng)濟發(fā)展、提高污水處理技術(shù)，同時考慮水資源利用情況，從模擬結(jié)果看來，此情景較符合河北省未來發(fā)展情況。

3.4情景對比結(jié)果

從表5可知，在預(yù)測期內(nèi)，情景W5的平均需水量僅大于情景W3的215.53億m3，而小于情景W1、W2和W4；情景W5的平均缺水程度是所有情景中最低，均小于情景W1、W2、W3和W4；情景W5的平均供需缺口在所有情景中處于倒數(shù)第二，僅大于情景W3的26.50億m3；情景W5的平均污水回用量僅小于情景W4（29.66億m3），大于情景W1、W2和W3；情景W5的平均GDP增長率大于情景W1、W3和W4（均為5%），僅小于情景W2的8%。綜上所述，較其余4種情景，情景W5既能以適當(dāng)?shù)乃俣缺ＷC未來經(jīng)濟發(fā)展，又能緩解水資源利用壓力，是最符合河北省未來的發(fā)展情景。

3.5政策建議

從情景對比結(jié)果可知，情景W5在模擬期內(nèi)平均需水量為232.95億m3，大于情景W3（215.53億m3），小于情景W1和W4（239.47億m3）、情景W2（258.51億m3）；情景W5的平均缺水程度為0.129 0，小于其余4種情景（W1為0.243 8，W2為0.302 6，W3為0.137 5，W4為0.146 7）。故而，在5種情景中，情景W5綜合考慮人口、社會經(jīng)濟、節(jié)水技術(shù)、居民用水習(xí)慣等因素，是符合綠色發(fā)展理念，也是研究區(qū)在未來應(yīng)該采納的發(fā)展方案。研究區(qū)的政府決策人員在制定社會經(jīng)濟發(fā)展和水資源保護等政策時，應(yīng)該考慮情景W5所提出的社會經(jīng)濟發(fā)展的方案。

河北省的經(jīng)濟發(fā)展對國家區(qū)域發(fā)展起到至關(guān)重要的作用。然而，河北省缺水問題將極大地制約當(dāng)?shù)厣鐣?jīng)濟發(fā)展。在現(xiàn)有的水資源情況下，一味地發(fā)展社會經(jīng)濟而不考慮水資源承載力，顯然對河北省的長遠(yuǎn)發(fā)展是極為不利。因此，河北省在未來的發(fā)展過程中，發(fā)展社會經(jīng)濟的同時應(yīng)兼顧對水資源的保護。情景W5既保證經(jīng)濟發(fā)展有一個適當(dāng)?shù)脑鲩L速度，又采取一定的措施保護水資源，有效地緩解缺水壓力，適合研究區(qū)的未來發(fā)展，可以給政府有關(guān)部門制定政策提供參考。

4結(jié)論

a）5種情景的總需水量均呈遞增趨勢。其中，情景W2（均值為258.51億m3）的增長速度最快，情景W1（均值為239.47億m3）、情景W4（均值為239.47億m3）、情景W5（均值為232.95億m3）的增長速度次之，情景W3（均值為215.53億m3）的增長速度慢；污水回用量均值依次為21.63億、26.31億、18.95億、29.66億、26.85億m3。

b）5種情景下水資源供需缺口均值依次為47.77億、62.12億、26.50億、31.86億、28.09億m3，缺水程度均值依次為0.243 8、0.302 6、0.137 5、0.146 7、0.129 0，表明河北省在未來仍存在缺水問題，這是由河北省是典型的資源型缺水省份的實際情況所決定，現(xiàn)有水資源不足以支撐未來發(fā)展，但通過采取積極的措施可以緩解未來水資源利用壓力。

c）尋找河北省水資源供需矛盾問題的解決方案，情景W5是最應(yīng)該被采納的方案，既能保證河北省在未來以適度的速度發(fā)展經(jīng)濟，又能緩解水資源利用壓力。在未來水資源利用過程中，需加強水資源規(guī)劃管理、推廣節(jié)水技術(shù)、加強污水處理技術(shù)、發(fā)展水資源多元化利用，才能實現(xiàn)水資源可持續(xù)利用。

參考文獻：

［1］GOHARI A，ESLAMIAN S，MIRCHI A，et al.Water transfer as a solution to water shortage：a fix that can backfire［J］.Journal of Hydrology，2013，491：23-39.

［2］邱國玉，張曉楠.21世紀(jì)中國的城市化特點及其生態(tài)環(huán)境挑戰(zhàn)［J］.地球科學(xué)進展，2019，34（6）：640-649.

［3］JIANG Y.China′s water scarcity［J］.Journal of Environmental Management，2009，90（11）：3185-3196.

［4］馬丁，李碩.中國地表水水質(zhì)變化趨勢及治理政策應(yīng)對［J］.中國人口·資源與環(huán)境，2023，33（5）：27-39.

［5］CHENG H F，HU Y A.Improving China’s water resources management for better adaptation to climate change［J］.Climatic Change，2012，112（2）：253-282.

［6］曹曉磊，楊健，梁坤，等.河北省水資源現(xiàn)狀及保護對策［J］.水科學(xué)與工程技術(shù)，2014（4）：45-48.

［7］康愛榮，于法穩(wěn).河北省農(nóng)業(yè)水資源可持續(xù)利用的對策研究［J］.唐山學(xué)院學(xué)報，2007（4）：38-42，52.

［8］余嶺，畢星，付祥偉.河北省水資源短缺預(yù)測研究：基于離散二階差分方法分析［J］.河南科學(xué)，2015，33（4）：614-620.

［9］刁子乘，趙晶，韓宇平，等.河北省城鎮(zhèn)居民生活層次化用水的影響因素分析［J］.中國農(nóng)村水利水電，2023（6）：215-221.

［10］吳彎.城鎮(zhèn)需水量預(yù)測方法研究［D］.廣州：華南理工大學(xué)，2015.

［11］秦歡歡，黃碧賢，吳昊，等.基于SD模型和情景分析法的華北平原需水量預(yù)測［J］.節(jié)水灌溉，2022（2）：59-65.

［12］ZHAI Y Z，WANG J S，TENG Y G，et al.Water demand forecasting of Beijing using the time series forecasting method［J］.Journal of Geographical Sciences，2012，22（5）：919-932.

［13］吳凡，陳伏龍，張志君，等.基于改進灰色模型的石河子市需水量預(yù)測［J］.水資源與水工程學(xué)報，2020，31（3）：76-82.

［14］賈振亞.基于BP-神經(jīng)網(wǎng)絡(luò)的晉中市需水預(yù)測及水資源優(yōu)化配置［D］.太原：太原理工大學(xué)，2021.

［15］張惠芬.城市生態(tài)環(huán)境需水量研究［D］.昆明：昆明理工大學(xué)，2012.

［16］孔祥仟，陳園，劉博懿，等.基于主成分分析的神經(jīng)網(wǎng)絡(luò)在需水預(yù)測中的應(yīng)用［J］.水電能源科學(xué)，2018，36（4）：26-28.

［17］張騰，張震，徐艷.基于SD模型的海淀區(qū)水資源供需平衡模擬與仿真研究［J］.中國農(nóng)業(yè)資源與區(qū)劃，2016，37（2）：29-36.

［18］SUN Y H，LIU N N，SHANG J X，et al.Sustainable utilization of water resources in China：A system dynamics model［J］.Journal of Cleaner Production，2017，142：613-625.

［19］秦歡歡，賴冬蓉，萬衛(wèi)，等.基于系統(tǒng)動力學(xué)的北京市需水量預(yù)測及缺水分析［J］.科學(xué)技術(shù)與工程，2018，18（21）：175-182.

［20］曹琦，陳興鵬，師滿江.基于SD和DPSIRM模型的水資源管理模擬模型：以黑河流域甘州區(qū)為例［J］.經(jīng)濟地理，2013，33（3）：36-41.

［21］胡棋智.基于系統(tǒng)動力學(xué)的武漢市水資源管理策略研究［D］.武漢：華中科技大學(xué)，2023.

［22］張百川，張揚，李政，等.湖泊流域水環(huán)境–經(jīng)濟系統(tǒng)耦合模擬方法研究［J］.北京大學(xué)學(xué)報（自然科學(xué)版），2022，58（5）：937-948.

［23］陶夢婷.基于系統(tǒng)動力學(xué)的裝配式建筑施工安全風(fēng)險評價與控制研究［D］.西安：西安理工大學(xué)，2022.

［24］殷悅.基于利益相關(guān)者和SD的裝配式建筑供應(yīng)鏈績效評價與優(yōu)化研究［D］.鎮(zhèn)江：江蘇科技大學(xué)，2022.

［25］程志友，羅鋒，吳蕊，等.基于系統(tǒng)動力學(xué)模型的港口?；愤\輸管理演化博弈［J］.上海海事大學(xué)學(xué)報，2023，44（3）：71-78，99.

［26］楊子江，韓偉超，楊恩秀.昆明市水資源承載力系統(tǒng)動力學(xué)模擬［J］.長江流域資源與環(huán)境，2019，28（3）：594-602.

［27］李騰，劉晴，隋鵬.河北省水資源變化及可持續(xù)性評估［J］.中國農(nóng)業(yè)大學(xué)學(xué)報，2021，26（6）：11-20.

［28］段一明，張戈，于大濤.河北省水資源脆弱性及其調(diào)控分析［J］.國土與自然資源研究，2021（1）：84-85.

［29］黃麗想，王健泉，邸齊夢，等.基于生態(tài)足跡和SD模型的山東省水資源可持續(xù)利用評價［J］.人民珠江，2024，45（1）：1-12.

［30］秦劍.水環(huán)境危機下北京市水資源供需平衡系統(tǒng)動力學(xué)仿真研究［J］.系統(tǒng)工程理論與實踐，2015，35（3）：671-676.

（責(zé)任編輯：向飛）