秦歡歡，黃碧賢，吳 昊，甘家宇，周宇卓，熊文躍，易振貴

（1.東華理工大學(xué)核資源與環(huán)境國(guó)家重點(diǎn)實(shí)驗(yàn)室，南昌330013；2.東華理工大學(xué)水資源與環(huán)境工程學(xué)院，南昌330013）

0 引言

一般來(lái)說(shuō)，需水量預(yù)測(cè)是區(qū)域水資源可持續(xù)利用的重要依據(jù)[1]，涉及到諸多復(fù)雜的水文、氣象、經(jīng)濟(jì)、社會(huì)、工程、技術(shù)等因素[2]。如何準(zhǔn)確預(yù)測(cè)某個(gè)區(qū)域未來(lái)的需水量，是水利管理部門(mén)必須要考慮的關(guān)鍵問(wèn)題，對(duì)于該區(qū)域未來(lái)可持續(xù)的水資源利用具有十分重要的意義[3,4]。用于需水量預(yù)測(cè)的方法有很多，如指標(biāo)法[3]、時(shí)間序列法[5]、神經(jīng)網(wǎng)絡(luò)法[6]、定額法[7]、灰色預(yù)測(cè)模型法[8]、混沌預(yù)測(cè)法[9]、系統(tǒng)動(dòng)力學(xué)法[2,10-13]等。然而，與系統(tǒng)動(dòng)力學(xué)法相比，其他方法的應(yīng)用過(guò)程對(duì)于影響需水量預(yù)測(cè)的諸多復(fù)雜因素的考慮不足，無(wú)法定量刻畫(huà)這些因素之間相互作用、相互反饋的關(guān)系，也就無(wú)法捕獲需水量預(yù)測(cè)過(guò)程中的動(dòng)態(tài)變化過(guò)程[14,15]。系統(tǒng)動(dòng)力學(xué)法較好地解決了這個(gè)問(wèn)題，被廣泛地應(yīng)用于區(qū)域需水量預(yù)測(cè)[12-19]。

情景分析法是假定某種現(xiàn)象或趨勢(shì)將持續(xù)到未來(lái)的前提下，對(duì)預(yù)測(cè)對(duì)象可能出現(xiàn)的情況或引起的后果做出預(yù)測(cè)的方法。情景分析法最大的優(yōu)勢(shì)是使管理者能發(fā)現(xiàn)未來(lái)變化的某些趨勢(shì)和避免兩個(gè)最常見(jiàn)的決策錯(cuò)誤：過(guò)高或過(guò)低估計(jì)未來(lái)的變化及其影響。自誕生以來(lái)，情景分析法在需水量預(yù)測(cè)、水資源利用等諸多方面得到了廣泛的應(yīng)用[20-22]。

作為我國(guó)的政治、經(jīng)濟(jì)和文化中心，華北平原屬于“資源型”缺水地區(qū)，是我國(guó)水資源供需矛盾最突出的區(qū)域之一。由于華北平原人口密集、經(jīng)濟(jì)發(fā)達(dá)，其人均水資源占有量只有全國(guó)平均值得23%[23]，加上區(qū)域內(nèi)降水量的時(shí)空分布極不均勻，容易發(fā)生旱澇災(zāi)害[24]，這進(jìn)一步加劇了華北平原的水資源短缺問(wèn)題，導(dǎo)致地下水的大量超采，以滿(mǎn)足工農(nóng)業(yè)及生活用水的需求。水資源短缺給華北平原帶來(lái)了諸多生態(tài)環(huán)境問(wèn)題，如地下水超采、地面沉降、塌陷等[25]，深刻地影響著華北平原社會(huì)經(jīng)濟(jì)的發(fā)展與資源環(huán)境的可持續(xù)管理。

為此，本文以華北平原為研究對(duì)象，通過(guò)構(gòu)建華北平原需水量系統(tǒng)動(dòng)力學(xué)模型，在模型中考慮影響需水量的諸多社會(huì)經(jīng)濟(jì)和工程技術(shù)因素，采用情景分析法，定量預(yù)測(cè)在不同的發(fā)展情景下華北平原未來(lái)30年的需水量和水資源供需平衡狀況。通過(guò)對(duì)不同情景的對(duì)比，挑選出既符合可持續(xù)發(fā)展內(nèi)涵又滿(mǎn)足華北平原實(shí)際的發(fā)展情景，作為水利管理部門(mén)決策的參考，為華北平原未來(lái)水資源的可持續(xù)利用和管理提供科學(xué)的依據(jù)和有效的建議。

1 材料與方法

1.1 研究方法簡(jiǎn)介

系統(tǒng)動(dòng)力學(xué)(System Dynamics， SD)是1956年創(chuàng)立的研究方法，距今已有60 多年的歷史。SD 是一門(mén)分析研究信息反饋系統(tǒng)的學(xué)科，也是一門(mén)認(rèn)識(shí)系統(tǒng)問(wèn)題和解決系統(tǒng)問(wèn)題的交叉綜合學(xué)科。從系統(tǒng)方法論來(lái)說(shuō)，SD 是結(jié)構(gòu)的方法、功能的方法和歷史的方法的統(tǒng)一。它基于系統(tǒng)論，吸收了控制論、信息論的精髓，是一門(mén)綜合自然科學(xué)和社會(huì)科學(xué)的橫向?qū)W科[26]。

SD 對(duì)問(wèn)題的理解，是基于系統(tǒng)行為與內(nèi)在機(jī)制間的相互緊密的依賴(lài)關(guān)系，并且透過(guò)數(shù)學(xué)模型的建立與操弄的過(guò)程而獲得的，逐步發(fā)掘出產(chǎn)生變化形態(tài)的因果關(guān)系。隨著相關(guān)理論、方法和技術(shù)的不斷發(fā)展與完善，SD 方法已經(jīng)在需水量預(yù)測(cè)、水資源管理、生態(tài)環(huán)境規(guī)劃等諸多領(lǐng)域得到了較廣泛的應(yīng)用[1-2,10,14,27]，在處理多重復(fù)雜反饋系統(tǒng)的問(wèn)題方面具有顯著優(yōu)勢(shì)。

1.2 研究區(qū)概況

總面積14 萬(wàn)km2的華北平原是我國(guó)三大平原（東北平原、華北平原和長(zhǎng)江中下游平原）之一，指的是黃河以北、燕山以南和太行山以東的沖積平原區(qū)，即36°～41°N、114°30′～118°30′E之間地帶，包括京津冀魯豫等省市的21個(gè)地級(jí)市210個(gè)縣（區(qū)），是我國(guó)的政治、經(jīng)濟(jì)、文化中心。由于半濕潤(rùn)半干旱的氣候原因、密集的人口壓力和發(fā)達(dá)的經(jīng)濟(jì)水平，華北平原幾乎每年都要遭受旱災(zāi)的威脅，是我國(guó)水資源壓力最大的地區(qū)。

華北平原冬春寒冷干燥，夏季炎熱多雨，多年平均降水量為554 mm，6-9月的汛期降水量占全年的75%以上，多年平均蒸散發(fā)為1 550 mm[28]。研究區(qū)內(nèi)有大小河流近60 條，包括潮白河、永定河、灤河、薊運(yùn)河等。近20年以來(lái)，區(qū)內(nèi)大部分河流常年干涸或僅在汛期短時(shí)過(guò)流，或成為城鎮(zhèn)生活及工業(yè)的排污河。

2 華北平原需水量預(yù)測(cè)SD模型

2.1 概念模型與系統(tǒng)流圖

華北平原水資源系統(tǒng)包括需水和供水兩部分（圖1）。需水系統(tǒng)由農(nóng)業(yè)、生活和工業(yè)需水三部分組成，其中農(nóng)業(yè)需水由灌溉與牲畜需水組成。農(nóng)業(yè)需水受灌溉面積、牲畜數(shù)量及用水定額影響；生活需水由總?cè)丝诩叭丝谟盟?xí)慣決定；工業(yè)需水取決于工業(yè)總產(chǎn)值和工業(yè)用水定額。供水系統(tǒng)由地表水、地下水、南水北調(diào)工程引水、灌溉回歸水及污廢水回用等組成。華北平原水資源現(xiàn)狀決定了地表水和地下水可供水量，南水北調(diào)工程增加了地表水供水量，灌溉方法與技術(shù)、污廢水處理與回用技術(shù)決定了灌溉回歸水量和污廢水回用量。供需水不平衡會(huì)造成缺水問(wèn)題，缺水程度在一定程度影響各部分需水量。通過(guò)缺水程度可以反應(yīng)華北平原水資源短缺嚴(yán)重程度，從而以定量方式衡量華北平原水資源利用情況。

系統(tǒng)流圖是SD 基本變量和符號(hào)的有機(jī)組合。根據(jù)系統(tǒng)內(nèi)部各因素之間因果關(guān)系設(shè)計(jì)系統(tǒng)流圖，可以將系統(tǒng)內(nèi)部各因果關(guān)系中未能反映出來(lái)的不同變量的性質(zhì)和特點(diǎn)反映出來(lái)。通過(guò)流圖中關(guān)系的量化（采用SD 模擬軟件VENSIM Professional 進(jìn)行模型的構(gòu)建與量化）就達(dá)到政策模擬的目的。華北平原水資源系統(tǒng)是一個(gè)復(fù)雜的系統(tǒng)，它與人口、經(jīng)濟(jì)、社會(huì)關(guān)系密切，因此把它分為人口、農(nóng)業(yè)、工業(yè)、水環(huán)境、水資源5個(gè)子系統(tǒng)，各個(gè)子系統(tǒng)相互聯(lián)系，相互影響（圖2）。

圖2 華北平原需水量預(yù)測(cè)系統(tǒng)流圖Fig.2 Flow chart of water demand prediction system in North China Plain

2.2 模型變量、參數(shù)與數(shù)據(jù)來(lái)源

在對(duì)系統(tǒng)中各個(gè)變量之間因果關(guān)系分析的基礎(chǔ)上，構(gòu)建了5個(gè)狀態(tài)方程、大量的速率方程和輔助方程以及表函數(shù)。通過(guò)這些方程，就把各變量之間的邏輯關(guān)系“翻譯”成為了數(shù)學(xué)語(yǔ)言，以顯示其定量關(guān)系。生活子系統(tǒng)、工業(yè)子系統(tǒng)和農(nóng)業(yè)子系統(tǒng)主要變量的計(jì)算方式見(jiàn)公式(1)~(2)。

式中：WT、WD、WI和WA分別表示總需水量、生活需水量、工業(yè)需水量和農(nóng)業(yè)需水量；WUD和WRD分別表示城鎮(zhèn)生活需水量和農(nóng)村生活需水量；WIrr、WLL和WSL分別表示灌溉需水量、大牲畜需水量和小牲畜需水量；PUD和PRD分別表示城鎮(zhèn)人口和農(nóng)村人口；OI表示工業(yè)產(chǎn)值；AIrr、PLL和PSL分別表示灌溉面積、大牲畜頭數(shù)和小牲畜頭數(shù)；QUD和QRD分別表示城鎮(zhèn)生活用水定額和農(nóng)村生活用水定額；QI表示萬(wàn)元工業(yè)產(chǎn)值用水量；QIrr、QLL和QSL分別表示灌溉定額、大牲畜用水定額和小牲畜用水定額；ST、SS、SG和SSNWDP分別表示總供水量、地表水供水量、地下水供水量和南水北調(diào)工程供水量；SRW和SWW分別表示灌溉回歸水和污廢水回用量；WD表示供需差額。

模擬時(shí)間為2000-2050年，基準(zhǔn)年是2000年，時(shí)間步長(zhǎng)是1 a，結(jié)果輸出間隔也是1 a。模擬時(shí)間分為兩階段：模型校準(zhǔn)與模型預(yù)測(cè)。模型校準(zhǔn)階段（2000-2018年）是根據(jù)歷史數(shù)據(jù)對(duì)模型進(jìn)行校準(zhǔn)，從而獲得和歷史數(shù)據(jù)擬合度比較好的模型參數(shù)，以供模型預(yù)測(cè)階段使用。模型預(yù)測(cè)階段（2019-2050年）是利用前一階段獲得的參數(shù)，對(duì)不同發(fā)展情景下華北平原未來(lái)的需水量和水資源供需平衡情況進(jìn)行模擬預(yù)測(cè)。

模型中參數(shù)有常數(shù)類(lèi)、表函數(shù)和初始值。常數(shù)是指在模型模擬過(guò)程中不發(fā)生改變的量；表函數(shù)是用插值方法（本文選的是Euler 插值法）來(lái)確定無(wú)法用簡(jiǎn)單函數(shù)表示出來(lái)的量；初始值是模擬最初始時(shí)刻流位變量的值。表1是模型中常數(shù)和初始值，表2是模型中的表函數(shù)。

表1 模型中的常數(shù)和初始值Tab.1 Constants and initial values in the model

表2 模型校準(zhǔn)階段的表函數(shù)Tab.2 Table function at model calibration stage

華北平原需水量預(yù)測(cè)SD 模型中用到的數(shù)據(jù)可以分為社會(huì)經(jīng)濟(jì)數(shù)據(jù)、水資源數(shù)據(jù)及污水回用數(shù)據(jù)等三大類(lèi)，具體的數(shù)據(jù)來(lái)源包括：①人口、工業(yè)、灌溉、牲畜、作物等社會(huì)經(jīng)濟(jì)數(shù)據(jù)來(lái)源于2001-2019年華北平原各省市（北京市、天津市、河北省、山東省、河南省，下同）的統(tǒng)計(jì)年鑒，分別搜集華北平原包含的21 個(gè)地級(jí)市的社會(huì)經(jīng)濟(jì)數(shù)據(jù)，經(jīng)加總得到華北平原的數(shù)據(jù)；②地表水供水量、地下水可開(kāi)采量、用水定額、灌溉定額等水資源數(shù)據(jù)來(lái)源于華北平原各省市2000-2018年的水資源公報(bào)和2000-2018年《海河流域水資源公報(bào)》[29]；③污水回用數(shù)據(jù)，來(lái)源于2001-2019年華北平原各省市的統(tǒng)計(jì)年鑒。

2.3 模型校準(zhǔn)

模型校準(zhǔn)是將歷史參數(shù)輸入到模型經(jīng)運(yùn)行后的仿真結(jié)果，與實(shí)際發(fā)生的行為數(shù)據(jù)進(jìn)行比較，驗(yàn)證其吻合程度，為模型行為模擬的可靠性和準(zhǔn)確性做出判斷。由于SD 模擬過(guò)程中采用插值的形式對(duì)表函數(shù)中間年份參數(shù)進(jìn)行取值，而短期歷史數(shù)據(jù)波動(dòng)較大，因此模型校準(zhǔn)著重于檢驗(yàn)?zāi)Ｐ湍M的趨勢(shì)與歷史數(shù)據(jù)表現(xiàn)的趨勢(shì)是否大致一致。本次校準(zhǔn)項(xiàng)目為人口子系統(tǒng)、經(jīng)濟(jì)子系統(tǒng)（工業(yè)、農(nóng)業(yè)）、環(huán)境子系統(tǒng)與資源子系統(tǒng)，被校準(zhǔn)變量有生活需水量、工業(yè)需水量、農(nóng)業(yè)需水量和總需水量。表3 列出了模型校準(zhǔn)結(jié)果。從表3 可以看出，模型模擬結(jié)果與歷史數(shù)據(jù)基本吻合，相對(duì)誤差基本都在5%以?xún)?nèi)，這說(shuō)明模擬結(jié)果與和歷史數(shù)據(jù)比較符合，參數(shù)可運(yùn)用于模型預(yù)測(cè)階段。

表3 模型校準(zhǔn)結(jié)果 億m3Tab.3 Results of model calibration

2.4 發(fā)展情景設(shè)計(jì)

需水量變化受到諸多因素影響，不同因素的組合可能導(dǎo)致不同的水資源供需情況。為此，本文從影響華北平原需水量主要因素著手，設(shè)計(jì)了趨勢(shì)保持型情景S1、經(jīng)濟(jì)發(fā)展型情景S2、節(jié)約用水型情景S3和可持續(xù)發(fā)展型情景S4共4種情景，用來(lái)預(yù)測(cè)華北平原未來(lái)在不同發(fā)展情景下的需水量和水資源供需平衡狀況，具體如下：

（1）發(fā)展情景S1：又稱(chēng)為趨勢(shì)保持型情景。此情景假設(shè)華北平原的發(fā)展政策及系統(tǒng)結(jié)構(gòu)不發(fā)生大的調(diào)整，城市人口、社會(huì)經(jīng)濟(jì)發(fā)展速度沒(méi)有太大的變化?？?cè)丝谠鲩L(zhǎng)率在未來(lái)設(shè)定為0.68%，工業(yè)總產(chǎn)值增長(zhǎng)率在未來(lái)設(shè)定為2%，城鎮(zhèn)化率在2030年、2040年和2050年分別為70%、75%和80%。

（2）發(fā)展情景S2：又稱(chēng)為經(jīng)濟(jì)發(fā)展型情景。此情景假設(shè)在當(dāng)前和今后相當(dāng)長(zhǎng)的一段時(shí)間內(nèi)，發(fā)展經(jīng)濟(jì)仍是華北平原的重中之重。為此，把2030年、2040年和2050年工業(yè)總產(chǎn)值增長(zhǎng)率分別提高到5%、7%和8%，人口增長(zhǎng)率在2030年、2040年和2050年分別提高為0.8%、0.9%和0.9%，其他參數(shù)與情景S1相同。

（3）發(fā)展情景S3：又稱(chēng)節(jié)約用水型情景。假設(shè)華北平原將水資源保護(hù)列為重點(diǎn)，采用多種節(jié)約用水的方式，將城鎮(zhèn)生活用水定額和農(nóng)村生活用水定額分別降低10%，萬(wàn)元工業(yè)產(chǎn)值用水定額、灌溉定額、大牲畜用水定額和小牲畜用水定額分別降低20%，其他參數(shù)與情景S1相同。

（4）發(fā)展情景S4：又稱(chēng)為可持續(xù)發(fā)展型情景。在此情景下，既強(qiáng)調(diào)經(jīng)濟(jì)的發(fā)展，同時(shí)也注意對(duì)水資源的保護(hù)。人口增長(zhǎng)率在2030年、2040年和2050年分別提高為0.8%、0.9%和0.9%，把2030年、2040年和2050年工業(yè)總產(chǎn)值增長(zhǎng)率分別提高到5%、7%和8%，將城鎮(zhèn)生活用水定額和農(nóng)村生活用水定額分別降低10%，萬(wàn)元工業(yè)產(chǎn)值用水定額、灌溉定額、大牲畜用水定額和小牲畜用水定額分別降低20%，其他參數(shù)與情景S1相同。

3 結(jié)果與分析

3.1 需水量結(jié)果

圖3 是4 種情景下華北平原總需水量預(yù)測(cè)結(jié)果。從總需水量角度來(lái)看，現(xiàn)有狀況（情景S1）下華北平原2050年總需水量為421.8 億m3，而經(jīng)濟(jì)發(fā)展型情景（情景S2）、節(jié)約用水型情景（情景S3）和可持續(xù)發(fā)展型情景（情景S4）2050年總需水量分別為590.1、353.4 和489.1 億m3。預(yù)測(cè)期內(nèi)4 種情景平均總需水量分別為383.56、430.79、321.73和359.75 億m3。預(yù)測(cè)期內(nèi)，全部情景的總需水量都隨時(shí)間遞增，情景S2 的總需水量最大，情景S1和情景S4次之，情景S3最小。說(shuō)明在未來(lái)年份里，無(wú)論采取什么樣的發(fā)展情景，要使得經(jīng)濟(jì)保持發(fā)展（適度發(fā)展或者快速發(fā)展），總需水量都會(huì)隨著經(jīng)濟(jì)的發(fā)展而不斷增加，這是必然的趨勢(shì)。

圖3 總需水量預(yù)測(cè)結(jié)果Fig.3 Total water demand prediction results

具體到各用水部門(mén)，圖4 列出了模擬期內(nèi)4 種情景下農(nóng)業(yè)、工業(yè)和生活平均需水量柱狀圖。對(duì)于農(nóng)業(yè)需水量，情景S3 和S4（172.01 億m3）小于情景S1 和S2（215.02 億m3）；對(duì)于工業(yè)需水量，情景S2（104.10 億m3）最大，情景S4（83.39億m3）和S1（58.07 億m3）次之，情景S3（46.46 億m3）最??；對(duì)于生活需水量，情景S2（73.21 億m3）最大，情景S1（72.00億m3）和S4（65.89億m3）次之，情景S3（64.80億m3）最小。4種情景下，農(nóng)業(yè)需水量均為最大的用水部門(mén)，占總需水量的平均比例分別為56.06%、49.91%、53.47%和47.81%；工業(yè)需水量占總需水量的平均比例分別為15.14%、24.16%、14.44%和23.18%；生活需水量占總需水量的平均比例分別為18.77%、16.99%、20.14%和18.32%。因此，控制農(nóng)業(yè)需水量的增長(zhǎng)速度，在很大程度上能控制華北平原水資源供需平衡的狀況。

圖4 4種情景下各部門(mén)平均需水量柱狀圖Fig.4 Histogram of average water demand of each section in four scenarios

3.2 水資源供需平衡

從水資源供需平衡的角度來(lái)看（圖5），4種情景下華北平原2050年供需差額分別為35.2 億m3、18.1 億m3、-30.27 億m3和87.35 億m3?？傮w來(lái)看（圖5），情景S2 的供需差額最大，情景S1和情景S4次之，情景S3最小。可持續(xù)型發(fā)展情景（情景S4）只在2042-2050年出現(xiàn)缺水的情況，其余年份都不存在缺水問(wèn)題。4種情景在模擬末期基本都會(huì)缺水，在現(xiàn)有條件下可持續(xù)型發(fā)展型情景由于既考慮了經(jīng)濟(jì)發(fā)展，又考慮了水資源保護(hù)，缺水量只在模擬后期5年比情景S1 多，這是可持續(xù)型發(fā)展情景的優(yōu)勢(shì)。同時(shí)，情景S3 雖然缺水量最少，但經(jīng)濟(jì)發(fā)展也受到限制。在華北平原大力發(fā)展經(jīng)濟(jì)的宏觀(guān)背景下，這樣的情景顯然不會(huì)被決策者采納。

圖5 水資源供需差額預(yù)測(cè)結(jié)果Fig.5 Prediction results of water resources supply-demand balance

3.3 水環(huán)境模擬結(jié)果

4 個(gè)發(fā)展情景的側(cè)重點(diǎn)不同，污廢水回用量也不同（圖6）。情景S2 由于過(guò)分強(qiáng)調(diào)經(jīng)濟(jì)發(fā)展，必然導(dǎo)致污廢水排放量增多，也即導(dǎo)致污廢水回用量在4 種情景中最多；情景S3 由于采用了節(jié)約用水措施，此情景下污廢水回用量是最少的，而情景S4 由于既考慮了經(jīng)濟(jì)發(fā)展，又考慮了節(jié)約用水措施，該情景下污廢水回用量?jī)H比情景S3 略多。污廢水回用量多少可以作為環(huán)境好壞的一個(gè)衡量標(biāo)準(zhǔn)，而環(huán)境好壞與水資源有著直接關(guān)系。由此可見(jiàn)，情景S4 既注重經(jīng)濟(jì)的發(fā)展，也注重保護(hù)環(huán)境，可以很好地協(xié)調(diào)經(jīng)濟(jì)發(fā)展與水資源及環(huán)境保護(hù)之間的關(guān)系。

圖6 污廢水回用預(yù)測(cè)結(jié)果Fig.6 Prediction results of wastewater reuse

3.4 發(fā)展情景比較

表4是4種發(fā)展情景的綜合比較。在用水方面，情景S2用水最多，情景S3用水最少，而情景S1和情景S4居中。但是由于情景S4 只在模擬期后6年需水量大于情景S1（或缺水量大于情景S1），所以情景S4 用水總體比情景S1 少。在環(huán)境指標(biāo)方面，情景S2的污廢水回用量最多，情景S3最少，而情景S1和情景S4 居中。綜合來(lái)看，情景S4 在環(huán)境指標(biāo)方面優(yōu)于情景S1；從經(jīng)濟(jì)發(fā)展方面看，情景S2 發(fā)展速度最快，情景S3 由于強(qiáng)調(diào)水資源保護(hù)，其經(jīng)濟(jì)發(fā)展速度最慢，而情景S4 發(fā)展速度比情景S1要快（表4）。

綜合來(lái)看，在4 種情景中，情景S4 是最優(yōu)的發(fā)展情景（表4）。經(jīng)濟(jì)發(fā)展規(guī)模、總需水量、缺水量及環(huán)境保護(hù)等方面，情景S4 都最符合可持續(xù)發(fā)展內(nèi)涵，情景S4 是華北平原未來(lái)發(fā)展的最優(yōu)情景。在情景S4 下發(fā)展，華北平原保持適中的經(jīng)濟(jì)發(fā)展速度，對(duì)水資源及環(huán)境的需求和壓力也適中，經(jīng)濟(jì)與水資源環(huán)境協(xié)調(diào)可持續(xù)地發(fā)展。通過(guò)本文模擬可以發(fā)現(xiàn)，華北平原未來(lái)的缺水問(wèn)題不能夠得到徹底根除，這是由華北平原的實(shí)際情況所決定的。但是，采取節(jié)水措施和調(diào)整產(chǎn)業(yè)結(jié)構(gòu)，可以使華北平原缺水狀況得到改善。

表4 華北平原4種情景的對(duì)比排序（從大到小）Tab.4 Comparative ranking of four scenarios in North China Plain

4 結(jié)論

本文利用系統(tǒng)動(dòng)力學(xué)方法原理和情景分析的方法，建立了華北平原需水量預(yù)測(cè)的系統(tǒng)動(dòng)力學(xué)模型，設(shè)定了4個(gè)不同的發(fā)展情景，對(duì)華北平原未來(lái)需水量和水資源供需狀況進(jìn)行了模擬預(yù)測(cè)，可以得到以下結(jié)論：

（1）模擬期內(nèi)（2019-2050年），4種情景下華北平原平均總需水量分別為383.56、430.79、321.73和359.75 億m3，平均供需差額分別為0.68、41.62、-58.85和-25.88億m3。4種情景下，農(nóng)業(yè)、工業(yè)和生活需水量占總需水量的平均比例分別為56.06%、 49.91%、 53.47% 和47.81%， 15.14%、 24.16%、14.44%和23.18%以及18.77%、16.99%、20.14%和18.32%。

（2）伴隨著社會(huì)經(jīng)濟(jì)發(fā)展，華北平原需水量會(huì)越來(lái)越大。華北平原在未來(lái)不能徹底解決缺水問(wèn)題。如果華北平原在社會(huì)經(jīng)濟(jì)的發(fā)展中選擇可持續(xù)發(fā)展情景，那么在未來(lái)可以緩解缺水問(wèn)題。

（3）可持續(xù)型發(fā)展情景（S4）是華北平原未來(lái)應(yīng)該選擇的發(fā)展情景。過(guò)分強(qiáng)調(diào)經(jīng)濟(jì)發(fā)展（S2）、過(guò)分強(qiáng)調(diào)水資源保護(hù)（S3）或者保持目前趨勢(shì)（S1），勢(shì)必在水資源使用和經(jīng)濟(jì)發(fā)展之間失衡。從華北平原實(shí)際情況來(lái)看，決策者既不會(huì)采納過(guò)分強(qiáng)調(diào)經(jīng)濟(jì)發(fā)展的情景，也不會(huì)采納過(guò)分強(qiáng)調(diào)水資源保護(hù)的情景。只有保持經(jīng)濟(jì)發(fā)展與水資源保護(hù)的平衡，社會(huì)才能可持續(xù)地發(fā)展下去。

（4）水資源的可持續(xù)發(fā)展，不僅需要政府的政策引導(dǎo)，更加需要社會(huì)公眾的參與?！伴_(kāi)源”與“節(jié)流”是解決缺水問(wèn)題的兩個(gè)有效途徑，是辯證統(tǒng)一的。政府在“開(kāi)源”方面的作用比較大，而公眾則更多的充當(dāng)“節(jié)流”的主力軍。只有政府與公眾的協(xié)調(diào)一致，才能保證水資源的可持續(xù)發(fā)展，進(jìn)而保證社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展。