陳燕飛 ，鄒志科，王 娜，張烈濤

(1.長(zhǎng)江大學(xué)地球環(huán)境與水資源學(xué)院油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，武漢 430100；2.河海大學(xué)水文水資源學(xué)院，南京 210098； 3.襄陽(yáng)水文水資源勘測(cè)局，湖北 襄陽(yáng) 441000；4.湖北省咸寧市水文局，湖北 咸寧 437100)

南水北調(diào)工程是緩解中國(guó)部分地區(qū)水資源短缺的特大型跨流域引水工程。中線工程從長(zhǎng)江支流漢江上的丹江口水庫(kù)引水，沿伏牛山和太行山山前平原開(kāi)渠輸水，終點(diǎn)北京。它對(duì)緩解京津及華北地區(qū)水資源短缺，改善受水區(qū)生態(tài)環(huán)境，促進(jìn)華北地區(qū)經(jīng)濟(jì)和社會(huì)的持續(xù)穩(wěn)定發(fā)展具有巨大作用。但是，南水北調(diào)中線工程實(shí)施后，改變了漢江中下游的水資源分配，使?jié)h江中下游水環(huán)境產(chǎn)生重大影響[1-3]。改變了流域間水資源的自然地理分布，導(dǎo)致江漢丹江口水庫(kù)下泄量減少，漢江中下游干流流量減少，水位降低，流速變緩[4]，也有學(xué)者[5]研究了南水北調(diào)對(duì)漢江中下游生態(tài)環(huán)境的綜合影響。隨著經(jīng)濟(jì)發(fā)展，該地區(qū)工業(yè)用水、生活用水以及生態(tài)環(huán)境用水日益增大,排入漢江的污染也在增加，在南水北調(diào)中線工程實(shí)施的條件下，如何提高中下游的水資源循環(huán)利用率，及減少工業(yè)用水、農(nóng)業(yè)用水、生活用水的同時(shí)，工農(nóng)業(yè)產(chǎn)值能夠穩(wěn)定增長(zhǎng)，是該地區(qū)經(jīng)濟(jì)發(fā)展的迫切需要。

系統(tǒng)動(dòng)力學(xué)能處理高階次、非線性、多重反饋復(fù)雜時(shí)變系統(tǒng)的問(wèn)題；對(duì)參數(shù)要求不高，著重研究系統(tǒng)的結(jié)構(gòu)和動(dòng)態(tài)行為，不需要特別精確的數(shù)字；可對(duì)系統(tǒng)的動(dòng)態(tài)發(fā)展及趨勢(shì)進(jìn)行考察，可作長(zhǎng)期的動(dòng)態(tài)預(yù)測(cè)[6]；國(guó)內(nèi)最早研究有陳成鮮等[7]對(duì)我國(guó)水資源可持續(xù)發(fā)展系統(tǒng)動(dòng)力學(xué)的研究，王武科等[8]基于系統(tǒng)動(dòng)力學(xué)模型對(duì)渭河流域關(guān)中地區(qū)水資源調(diào)度系統(tǒng)優(yōu)化。本文通過(guò)運(yùn)用系統(tǒng)動(dòng)力學(xué)法對(duì)漢江中下游城市供水區(qū)進(jìn)行建模、模擬，構(gòu)建不同方案運(yùn)行不同結(jié)果進(jìn)行比較，得出有效的方案緩南水北調(diào)中線工程對(duì)解中下游水資源的供需關(guān)系造成的影響。

1 資 料

漢江中下游主要流經(jīng)的城市有襄樊，天門、潛江、孝感、仙桃及武漢等6個(gè)城市。本文數(shù)據(jù)主要來(lái)源于《湖北省統(tǒng)計(jì)局》及2006-2012 年的《湖北省水資源公報(bào)》，統(tǒng)計(jì)各城市的地表供水量、地下供水量、工業(yè)耗水量、農(nóng)業(yè)耗水量、農(nóng)業(yè)總產(chǎn)值及工業(yè)總產(chǎn)值。對(duì)2010、2020年水資源供需情況預(yù)測(cè)。漢江中下游各城市分布圖見(jiàn)圖1。

圖1 漢江中下游主要城市分布圖

2 漢江中下游水資源系統(tǒng)動(dòng)力學(xué)模型構(gòu)建

2.1 系統(tǒng)總體結(jié)構(gòu)

首先確定系統(tǒng)的邊界。本研究將系統(tǒng)的邊界確定為整個(gè)漢江中下游干流流經(jīng)的整個(gè)流域。南水北調(diào)中線工程調(diào)水后，漢江中下游干流供水區(qū)水資源的利用方式更加關(guān)系到整個(gè)區(qū)域的可持續(xù)發(fā)展。根據(jù)本文研究的需要，確定水資源總量、各產(chǎn)業(yè)水資源需求量、水資源的利用率、利用方式；各產(chǎn)業(yè)產(chǎn)值、萬(wàn)元產(chǎn)值耗水量、人口總量、污水排放量等對(duì)整個(gè)系統(tǒng)有較大影響，故全部劃在邊界之內(nèi)。

水資源系統(tǒng)動(dòng)力學(xué)模型是以水資源為核心，由農(nóng)業(yè)子模塊、工業(yè)子模塊、人口子模塊、污水子模塊共同組成，各模塊因素之間互為因果關(guān)系。污水子模塊主要研究供水量的變化所導(dǎo)致的污水排放量的變化，以及由于污水排放量的變化引起人口數(shù)量的變化等。人口子模塊屬于發(fā)展子系統(tǒng)，通過(guò)人口數(shù)量的增加和人口素質(zhì)的提高影響著生活需水量以及水技術(shù)進(jìn)步、工農(nóng)業(yè)耗水量；工農(nóng)業(yè)的發(fā)展增加了對(duì)水資源的需求量，從而加劇了水資源的短缺程度，水資源的短缺又會(huì)促進(jìn)工農(nóng)業(yè)節(jié)水，提高水的重復(fù)利用率，促進(jìn)水資源的開(kāi)發(fā)利用，從而又為經(jīng)濟(jì)的發(fā)展提供了源泉和動(dòng)力。構(gòu)建的系統(tǒng)因果關(guān)系圖如圖2所示。采用系統(tǒng)動(dòng)力學(xué)進(jìn)行模擬仿真，需要確定狀態(tài)變量、目標(biāo)變量、控制變量及輔助變量等，根據(jù)對(duì)漢江中下游干流水資源系統(tǒng)結(jié)構(gòu)的分析，確定了變量集，見(jiàn)表1。系統(tǒng)所涉及的變量名稱、單位及其SD方程式見(jiàn)表1。將狀態(tài)變量和表函數(shù)作為決策變量，在模型中輸入已經(jīng)確定的參數(shù)進(jìn)行模擬。

圖2 水資源系統(tǒng)總體因果關(guān)系圖

變量名變量類型方程式可供水資源量/億m3A 可供水資源量=地下水供水量+地表供水量+0．8×中水回用量水資源短缺量/億m3A 水資源短缺量=水資源需求量-可供水資源量水資源需求量/億m3A 水資源需求量=工業(yè)需水量+農(nóng)業(yè)需水量+生活需水量地下水供水量/億m3T 地下水供水量=Table(time)地表水供水量/億m3T 地表水供水量=Table(time)工業(yè)總產(chǎn)值/億元L 工業(yè)總產(chǎn)值=INTER(工業(yè)產(chǎn)值增長(zhǎng)速度，初始工業(yè)總產(chǎn)值)工業(yè)產(chǎn)值增長(zhǎng)速度/億元R工業(yè)產(chǎn)值增長(zhǎng)速度=工業(yè)總產(chǎn)值×工業(yè)產(chǎn)值增長(zhǎng)率×(1-3×水資源短缺量/可供水資源量×0．01)

續(xù)表1 變量名稱及SD方程

注：L為狀態(tài)變量；R為速率變量；A為輔助變量；T為表函數(shù)。

2.2 系統(tǒng)流圖

在圖2系統(tǒng)總體因果關(guān)系的基礎(chǔ)上，建立水資源可持續(xù)發(fā)展的系統(tǒng)動(dòng)力學(xué)流圖，如圖3所示。

圖3 漢江中下游水資源系統(tǒng)動(dòng)力學(xué)流程圖

3 模型檢驗(yàn)

3.1 參數(shù)確定

利用歷史統(tǒng)計(jì)數(shù)據(jù)求取漢江流域各城市的算術(shù)平均值。這些數(shù)據(jù)包含平均人口出生率(0.8%)、平均人口死亡率(0.56%)、平均工業(yè)用水利用率(0.7～0.85)、工業(yè)和生活污水排放系數(shù)(0.7～0.8)等。

3.2 模型的檢驗(yàn)

為保證模型的正確性和精確度，必須進(jìn)行模型檢驗(yàn)。運(yùn)用Vensim PIE軟件進(jìn)行模擬。模型檢驗(yàn)有4種方法：直觀檢驗(yàn)、運(yùn)行檢驗(yàn)、歷史檢驗(yàn)、靈敏度。由于數(shù)據(jù)獲取難度較大，僅選取武漢市、仙桃市工業(yè)總產(chǎn)值、農(nóng)業(yè)總產(chǎn)值指標(biāo)為代表進(jìn)行歷史驗(yàn)證。以2006年的數(shù)據(jù)為起始，模擬輸出2006-2012 年的數(shù)據(jù)，輸出結(jié)果與歷史統(tǒng)計(jì)數(shù)據(jù)相比，如果二者擬合程度較好，則驗(yàn)證了模型的有效性和精確性。仙桃的模型檢驗(yàn)效果見(jiàn)表2，武漢的模型檢驗(yàn)效果見(jiàn)表3。

由表2和表3可知，模型檢驗(yàn)最大誤差不超過(guò) 2.0%， 可以認(rèn)為模型模擬系統(tǒng)與實(shí)際系統(tǒng)具有較好的一致性。設(shè)定時(shí)間邊界為2006-2020年，時(shí)間步長(zhǎng)為1 a，進(jìn)行近期( 2010年) 從丹江口水庫(kù)調(diào)水95億m3條件下的系統(tǒng)模擬。運(yùn)用此系統(tǒng)對(duì)未來(lái)2006-2020年水資源需求量、水資源可供量預(yù)測(cè)，結(jié)果見(jiàn)圖4。

表2 2006-2012年仙桃工農(nóng)業(yè)總產(chǎn)值模擬檢驗(yàn)結(jié)果

表3 2006-2012年武漢工農(nóng)業(yè)總產(chǎn)值模擬檢驗(yàn)結(jié)果

圖4 仙桃市和武漢市水資源供求圖

由圖4可知， 2006-2020年，漢江中下游流域水資源供不應(yīng)求的趨勢(shì)不斷上升，其中武漢市水資源需求量增長(zhǎng)幅度大于仙桃市，武漢市的水資源供需缺口也比仙桃市大。究其原因，隨著社會(huì)經(jīng)濟(jì)的發(fā)展，城鎮(zhèn)化進(jìn)程的加快，近年來(lái)，武漢市的經(jīng)濟(jì)快速發(fā)展，武漢市的人口膨脹較大。同時(shí)南水北調(diào)中線工程的通水，又使得漢江中下游流域的可供水資源量減少，加劇了水資源的供需不平衡。因此，必須采取相應(yīng)措施，使水資源、社會(huì)經(jīng)濟(jì)可持續(xù)穩(wěn)定發(fā)展。

4 方案設(shè)計(jì)與模擬結(jié)果分析

4.1 方案設(shè)計(jì)

方案1：現(xiàn)狀延續(xù)。保持現(xiàn)在狀況和發(fā)展趨勢(shì)不變，不采取任何措施。

方案2：節(jié)水減排及污水再回生利用。提高人們的節(jié)水意識(shí)，降低生活用水定額，減少生活污水排放量，即降低生活污水排放系數(shù)，同時(shí)增大污水處理率和中水回用率。

方案3：企業(yè)工藝改造。制定優(yōu)惠政策，引導(dǎo)企業(yè)改進(jìn)設(shè)備，提高水資源的利用效率并減少工業(yè)廢水排放系數(shù)，降低工業(yè)萬(wàn)元產(chǎn)值耗水量(降低5%)。

方案4：農(nóng)業(yè)生產(chǎn)現(xiàn)代化。主要是將農(nóng)業(yè)的生產(chǎn)方式逐步現(xiàn)代化，減少農(nóng)業(yè)萬(wàn)元生產(chǎn)耗水量(降低5%)。

方案5：綜合方案。將方案1、2、3、4這4個(gè)方案結(jié)合起來(lái)，綜合治理。

4.2 結(jié)果分析

通過(guò)改變方案中的決策變量值(常數(shù)變量和表函數(shù))，進(jìn)行模擬，以水資源短缺量在不同方案下的模擬結(jié)果為例對(duì)各個(gè)方案對(duì)比分析。具體結(jié)果見(jiàn)圖5。

圖5 不同方案下仙桃市和武漢市水資源短缺量變化趨勢(shì)

從圖5中曲線變化趨勢(shì)可以看出，首先，空間比較，各種方案下，仙桃市的水資源短缺量較武漢市波動(dòng)較大。仙桃市在2006-2009年期間水資源短缺量呈降低的趨勢(shì)，2009年達(dá)到最低，之后快速增長(zhǎng)，到2015年達(dá)到最大，之后緩慢下降，2017年開(kāi)始緩慢增加。這主要與仙桃市人口增長(zhǎng)速度有關(guān)。武漢市人口，經(jīng)濟(jì)發(fā)展比較平穩(wěn)，故水資源短缺量也呈緩慢增長(zhǎng)趨勢(shì)。其次，方案比較，與方案1相比，其他4種方案都能有效減緩武漢市和仙桃市水資源短缺量，而方案5水資源短缺量最小，方案2、方案3和方案4次之。方案5不僅更能緩解水資源短缺量，并且同時(shí)更能促進(jìn)工農(nóng)業(yè)增長(zhǎng)。方案3和方案4比較，對(duì)于以工業(yè)區(qū)為主的城市，工業(yè)工藝改造后降低的萬(wàn)元耗水量和農(nóng)業(yè)生產(chǎn)現(xiàn)代化后降低萬(wàn)元耗水量的百分比相同的情況下，減少工業(yè)用水比減少農(nóng)業(yè)用水更能減緩武漢市和仙桃市水資源短缺量。方案2和方案3、4比較發(fā)現(xiàn)，日常生活用水所占比重很大，可以減緩武漢市和仙桃市水資源短缺量。因此，對(duì)于漢江中下游以工業(yè)區(qū)為主的城市采取相應(yīng)措施減少生活用水和工業(yè)用水額度是緩解水資源短缺的關(guān)鍵。

5 結(jié) 語(yǔ)

(1)采用系統(tǒng)動(dòng)力學(xué)法，建立了漢江中下游水資源短缺系統(tǒng)，運(yùn)用Vensim PIE建模、仿真。研究發(fā)現(xiàn)，2006-2020年，漢江中下游流域水資源供不應(yīng)求的趨勢(shì)不斷上升，其中武漢市水資源需求量增長(zhǎng)幅度大于仙桃市，武漢市的水資源供需缺口也比仙桃市大。特別是在2014年南水北調(diào)通水之后，水資源短缺速率不斷增高，水資源短缺越來(lái)越嚴(yán)重。

(2)空間比較，各種方案下，仙桃市的水資源短缺量較武漢市波動(dòng)較大。方案比較，與方案1相比，其他4種方案都能有效減緩武漢市和仙桃市水資源短缺量，而方案5水資源短缺量最小，方案2、方案3和方案4次之。對(duì)于漢江中下游以工業(yè)區(qū)為主的城市采取相應(yīng)措施減少生活用水和工業(yè)用水額度是緩解水資源短缺的關(guān)鍵。

(3)建議積極鼓勵(lì)和引導(dǎo)企業(yè)通改造落后的制造工藝，降低萬(wàn)元工業(yè)產(chǎn)值年耗水量；引進(jìn)國(guó)內(nèi)外先進(jìn)的節(jié)水灌溉技術(shù)，降低萬(wàn)元農(nóng)業(yè)產(chǎn)值年耗水量，確定城鎮(zhèn)居民生活合理用水量，推廣節(jié)水型生活用水器具。

(4)建議把漢江中下游流域納入漢江生態(tài)保護(hù)體系中，積極爭(zhēng)取國(guó)家的財(cái)政投入，加大對(duì)流域內(nèi)的生態(tài)環(huán)境治理，切實(shí)保護(hù)漢江中下游各城市的水資源安全。

□

[1] 劉丙軍，邵東國(guó),許明祥，等. 南水北調(diào)中線與漢江中下游地區(qū)的水資源利用關(guān)系研究[J].南水北調(diào)與水利科技,2003,1(6):6-9.

[2] 胡志芳，張利平.南水北調(diào)對(duì)江漢中下游水環(huán)境影響分析[J].水利與建筑工程學(xué)報(bào),2008,6(3):16-18.

[3] 陳姍姍,劉俊新,王國(guó)義，等.南水北調(diào)對(duì)江漢中下游水環(huán)境的影響與可持續(xù)發(fā)展研究[J]. 西北大學(xué)學(xué)報(bào),2011,(6) :78-80.

[4] 肖 嬋, 謝 平, 唐 濤,等.南水北調(diào)中線工程對(duì)漢江中下游的水文情勢(shì)影響分析[J].水文，2009,29(1):26-29.

[5] 高永年，高俊峰.南水北調(diào)中線工程對(duì)漢江中下游流域生態(tài)環(huán)境影響的綜合評(píng)價(jià)[J].地理科學(xué)進(jìn)展，2010,1(29):59-64.

[6] 蔡 林.系統(tǒng)動(dòng)力學(xué)在可持續(xù)發(fā)展研究中的應(yīng)用[M].北京：中國(guó)環(huán)境科學(xué)出版社，2008.

[7] 陳成鮮，嚴(yán)廣樂(lè).我國(guó)水資源可持續(xù)發(fā)展系統(tǒng)動(dòng)力學(xué)模型研究[J].上海理工大學(xué)學(xué)報(bào),2000,22(2):154-159.

[8] 王武科，李同升，徐冬平，等.基于SD模型的渭河流域關(guān)中地區(qū)水資源調(diào)度系統(tǒng)優(yōu)化[J].資源科學(xué)，2008,30(7):984-989．