系統(tǒng)動力學模型在成都市水生態(tài)承載力評估方面的應用

馬涵玉+黃川友+殷彤+朱國宇

摘要：基于成都市水生態(tài)現狀，為提高其水生態(tài)承載力，促進經濟發(fā)展和提高人民生活水平，運用系統(tǒng)動力學（SD）的方法，建立了成都市水生態(tài)-經濟-人口-水資源-水環(huán)境的耦合系統(tǒng)，模擬了現狀延續(xù)型、節(jié)約用水型、污染防治型和綜合協(xié)調型四種情景模式。模擬結果顯示：在模擬年限內（2014年-2020年），現狀延續(xù)型和污染防治型未能有效降低水生態(tài)承載限制系數，水生態(tài)問題將進一步加??；節(jié)約用水型和綜合防治型都可以降低水生態(tài)承載限制系數，但節(jié)約用水型不能顯著減低該系數，只有通過節(jié)約用水和污染防治相結合的綜合防治型，才可以更加有效的減低該系數，該情景模式是提高水生態(tài)承載力的最佳模式；到2020年，該情景模式下成都市的水生態(tài)承載限制系數下降為0.297，與2010年的相比，下降了59.4%。研究結果可為成都市水生態(tài)保護提供技術依據。

關鍵詞：SD模型；水生態(tài)承載力；決策變量；情景模擬

中圖分類號：TV213 文獻標識碼：A 文章編號：1672-1683（2017）04-0101-10

Abstract：In order to enhance the water ecological carrying capacity of Chengdu city whilst promoting the development of economy and improving the living standards of the people，System Dynamics （SD）method was applied to build a coupled system of water ecology-economy-population-water resources-water environment of the city based on the status quo of water ecology.Four scenario patterns including status-continuing，water-saving，pollution-controlling，and integrated-coordination were simulated by the SD model.Simulation results showed that：in the analog period （2014 to 2020），the status-continuing and pollution-controlling scenario patterns can not effectively reduce the limit coefficient of water ecological carrying capacity，resulting in further deterioration of water ecological problems.However，the water-saving and integrated-coordination scenario patterns can reduce the limit coefficient of water ecological carrying capacity.It should be noted that the water-saving scenario pattern can not significantly reduce the coefficient.The integrated-coordination pattern combining water-saving and pollution-controlling measures can reduce the coefficient more effectively.It is the best pattern to improve the water ecological carrying capacity.With the aforementioned recommended pattern，the limit coefficient of water ecological carrying capacity can be reduced to 0.297 for Chengdu City by 2020，down by 59.4% compared with that of 2010.The research results can provide technical basis for water ecological protection in Chengdu city.

Key words：SD model；water ecological carrying capacity；decision variable；scenario simulation

系統(tǒng)動力學最初是由美國麻省理工學院Forrester教授所創(chuàng)立，20世紀70年代末引入我國。經過30多年來的發(fā)展，它幾乎被應用于自然科學和人類社會的全部領域。當前，系統(tǒng)動力學在可持續(xù)發(fā)展和生態(tài)系統(tǒng)變化的預測研究中應用相當廣泛[1]，它作為系統(tǒng)科學理論與計算機仿真緊密結合、研究系統(tǒng)反饋與行為的一門科學，是系統(tǒng)科學與管理科學的一個重要分支[2]。我國的水生態(tài)承載力研究主要是沿著水資源承載力→水環(huán)境承載力→水生態(tài)承載力和環(huán)境承載力→生態(tài)承載力→水生態(tài)承載力這兩條主線逐步發(fā)展起來的[3]。目前對水生態(tài)承載力定義的敘述，不同的學者對此表述不一。李靖[4]等將水生態(tài)承載力表述為在一定的歷史階段，某一流域的水生態(tài)系統(tǒng)在滿足自身健康發(fā)展的前提下，在一定的環(huán)境背景下所能持續(xù)支持人類社會經濟發(fā)展規(guī)模的閾值。彭文啟[5]給出了水生態(tài)承載力的廣義定義：流域水生態(tài)承載力指維持良好狀態(tài)的流域水資源及水環(huán)境系統(tǒng)所能承受的一定水平的人口的最大數量和一定技術水平下的最大經濟規(guī)模，從水資源、水環(huán)境等水體理化特性及水生態(tài)特征等方面考慮了水生態(tài)承載力的定義。

成都市是我國西部地區(qū)重要的中心城市。隨著城市化建設和經濟的發(fā)展，成都市對水資源的需求，特別是工業(yè)水量的需求不斷增加，供需矛盾日益突出，加之水資源利用率低，回用率更低，入河污染物排放過多，城市區(qū)域水環(huán)境污染嚴重，致使成都市在建設西部經濟高地不斷前進的道路上，面臨的水生態(tài)問題也越來越突出。針對以上問題，許多學者用不同的理論與方法進行了分析與研究。薛小妮等[6]以可承載人口數量作為指標，利用三層次分析方法和納污模型分別進行了成都市的水資源和水環(huán)境承載力預測分析，但存在指標單一的問題。陳琳[7]等運用PSR模型建立水安全評價指標體系，采用層次分析法計算指標權重，最后用水安全評價綜合指數對成都市的水安全狀況進行評價。歐陽銘[8]等利用生態(tài)足跡理論計算了2000年-2009年的水資源承載力，并運用主成分分析法對影響承載力的因子進行了分析，而未預測未來水資源承載力的發(fā)展變化。蘇敏[9]采用主成分分析法對主成分進行特征值和貢獻率的計算，得到了影響成都市水資源承載力的主成分因子，而未對承載力進行量化研究。因此，本文運用Vensim軟件建立成都市水生態(tài)承載力SD模型來定性與定量地預測分析水生態(tài)承載力的變化情況，旨在分析影響成都市水生態(tài)承載力的兩個方面——水資源與水環(huán)境，預測水生態(tài)承載力的變化趨勢，為成都市節(jié)約用水和污染防治政策制定提供科學依據，提高水生態(tài)承載力。

1 區(qū)域概況

成都市位于東經102°54′-104°53′，北緯30°05′-31°26′，地處四川省中部，四川盆地西部，屬長江水系岷江及沱江流域。其東部為龍泉山低山丘陵區(qū)和金堂縣沱江以東丘陵區(qū)，中部為成都平原，西部為龍門山和邛崍山山區(qū)。全市東西最大橫距192 km，南北最大縱距166 km，總面積為12 121.3 km2 [10]。成都市屬于亞熱帶季風氣候區(qū)，其氣候特點為：夏無酷暑，冬無冰雪，氣候溫和，夏長東短，無霜期長，秋雨和夜雨較多，風速小，濕度大，云霧多。全市多年平均降水量為900～1300 mm，其雨量主要集中在7月-8月，約占全年降水量的一半，年平均氣溫在15.2 ℃ ～16.6 ℃左右。

成都市是我國西南部經濟最為發(fā)達的城市之一， 2015年年末全市常住人口1 465.8萬人，實現地區(qū)生產總值10 801.2億元，占全四川省生產總值的35.9%[11]。成都市水資源量較為豐富，根據2004年-2013年統(tǒng)計[12]，年均水資源總量為81.72億m3；總用水量為54.15億m3，水資源開發(fā)利用率高達66.3%，多來年農業(yè)用水、工業(yè)用水、生活用水比例為3∶[KG-*4]1∶[KG-*4]1，農業(yè)用水占據了總用水量的60%。在水環(huán)境方面，多年來COD、氨氮的排放量分別為19.73萬t、2.31萬t，水質污染情況嚴重[13]；城市水體富營養(yǎng)化指數為59.85[14]，為輕度富營養(yǎng)；2015年成都市重要江河湖泊水功能區(qū)水質達標率為35%。這一系列亟待解決的問題，促使我們必須加快實施最嚴格的水資源管理制度和推進成都市水生態(tài)文明建設，以此為指導來改善成都市的水生態(tài)承載能力。

2 材料與方法

2.1 水生態(tài)承載力系統(tǒng)耦合關系分析

本文所研究的水生態(tài)承載力系統(tǒng)是人類活動影響下的系統(tǒng)，它是由人類社會系統(tǒng)、經濟系統(tǒng)、自然生態(tài)系統(tǒng)組成的耦合系統(tǒng)，可以將其細分為5個子系統(tǒng)——人口子系統(tǒng)、經濟子系統(tǒng)、水資源子系統(tǒng)、水環(huán)境子系統(tǒng)和水生態(tài)子系統(tǒng)。人口子系統(tǒng)以人口數量為特征，經濟子系統(tǒng)以GDP為主要指標，水資源子系統(tǒng)發(fā)揮水量供給的作用，水環(huán)境子系統(tǒng)具有自凈納污的重要功能，水生態(tài)子系統(tǒng)是水資源子系統(tǒng)和水環(huán)境子系統(tǒng)兩者的結合，通過水量供給和自凈納污來發(fā)揮作用。各個子系統(tǒng)之間相互作用，相互影響，共同對水生態(tài)承載力系統(tǒng)產生影響，決定水生態(tài)承載力的發(fā)展走向。成都市水生態(tài)承載力系統(tǒng)耦合關系見圖1。

2.2 SD模型原理及方法

SD模型具有分析速度快、構造簡單、可運用非線性方程等優(yōu)點[15]，是進行情景分析和系統(tǒng)模擬的一個很好工具。社會-經濟-生態(tài)組成的水生態(tài)系統(tǒng)屬于非線性的高階次復雜時變系統(tǒng)，經常表現出千差萬別、反直觀的動態(tài)特性，Vensim軟件作為系統(tǒng)動力學的專業(yè)建模軟件，其建立的SD模型可以作為此類系統(tǒng)的“實驗室”，通過情景分析來定性與定量的研究此類系統(tǒng)。本文選取影響水生態(tài)承載力的五個指標——總需水量、COD入河總量、水資源承載限制系數、水環(huán)境承載限制系數和水生態(tài)承載限制系數，進行以需水狀況和水環(huán)境納污水平為主的水生態(tài)承載力預測研究，其中水資源分析是建立在水資源開發(fā)利用率不超過40%的基礎上，水環(huán)境分析是以COD水環(huán)境容量[16]為衡量標準進行具體的研究分析。運用Vensim軟件建立成都市水生態(tài)承載力SD模型，并通過情景分析來定性與定量地分析水生態(tài)承載力的變化情況。

2.3 數據來源

本文中所使用的人口、工業(yè)GDP、糧食作物種植面積、禽畜養(yǎng)殖數量等資料來源于文獻[10，17-18]；水資源總量和工業(yè)萬元GDP用水量來自于文獻[12]；生活用水和農業(yè)用水定額數據根據文獻[19]確定；居民生活和農業(yè)的COD排放系數來自于文獻[20]；工業(yè)污水和污水處理廠排放的COD濃度數據取自于文獻[21]； COD水環(huán)境容量數據是按一維模型分別計算成都市每個水功能區(qū)的納污能力，匯總得到 [6]。并補充參考全國、四川及其他地方的相關數據和資料[22-24]，作為參考資料。

3 成都市水生態(tài)承載力SD模型的構建

3.1 模型邊界及基準年的確定

本模型將對水生態(tài)有直接、重要影響的因素劃在邊界范圍之內，同時考慮政策的可實施性和行政邊界的完整性，將模型邊界定義為成都市的行政區(qū)范圍。系統(tǒng)模型的模擬年限為2010年-2020年，其中2010年-2013年為歷史檢驗年份，2014年-2020年為預測年份，規(guī)劃基準年為2010年，模擬時間步長為年。模擬過程：首先以2010年的各項指標值為初始數據模擬預測2010年-2013年的各項指標值，然后檢驗模型的有效性，在此基礎上對模型進行測試；在各項模擬指標誤差都在某一范圍內后，進一步預測2014年-2020年成都市在不同情景方案下的水環(huán)境納污量和總需水量。

3.2 模型的結構劃分及變量方程式

3.2.1 模型的結構劃分及反饋回路

水生態(tài)承載力系統(tǒng)屬于高階非線性的復雜時變系統(tǒng)，里面主要包含納污和用水兩大子系統(tǒng)，經濟快速發(fā)展、城市常住人口增長和城鎮(zhèn)化率提高是城市排污量和用水量增長的主要驅動因素[25]。模型組成結構包括人口模塊、農業(yè)模塊、工業(yè)模塊、污水處理及回用模塊[26]、水資源模塊、水環(huán)境模塊、水生態(tài)模塊。

（1）人口模塊。

人口模塊是根據區(qū)域人口現狀，進行未來年份人口變化量預測，分為城鎮(zhèn)人口和農村人口預測，用來預測人口數量變化帶來的生活需水量和生活排污量的變化。其中城鎮(zhèn)居民生活污水進入污水處理廠處理后排放，作為點源污染進入自然水體；農村生活污水排放不經任何處理直接排放，作為面源污染進入自然水體。

（2）農業(yè)模塊。

農業(yè)主要包括種植業(yè)和禽畜養(yǎng)殖業(yè)，農業(yè)模塊主要根據禽畜數量、農作物種植面積和有效灌溉面積的變化量進行需水量和排污量的預測。其中種植業(yè)由于不同農作物用水量差異[27]比較大，將其分為水稻、其他農作物、油菜、蔬菜和水果五大類，分別進行需水量的計算；禽畜養(yǎng)殖業(yè)的禽畜數量通過豬與其他禽畜的數量換算關系統(tǒng)一折合成豬的數量；禽畜養(yǎng)殖業(yè)和種植業(yè)污染物作為面源污染進入自然水體。

（3）工業(yè)模塊。

工業(yè)模塊主要通過工業(yè)GDP和萬元工業(yè)增加值用水量來進行需水量和排污量的預測。其中工業(yè)污水的排放主要分兩部分：未經處理直排的量和經污水處理廠處理后排放的量，兩者都作為點源污染最終進入自然水體。

（4）污水處理及回用模塊。

污水處理及回用模塊主要包括城鎮(zhèn)居民生活污水處理模塊和工業(yè)污水處理模塊。污水處理廠的存在會大大降低污染物進入自然水體的風險，且經處理的污水可以實現一定程度的重新利用，增加水資源供給，本模型也設置了污水回用量。

（5）水資源模塊。

水資源模塊預測的總需水量主要包括農業(yè)需水量、工業(yè)需水量、居民生活需水量和生態(tài)需水量四部分。根據成都市的實際情況，將農業(yè)需水量分為灌溉需水量和禽畜養(yǎng)殖需水量；居民生活需水量分為城鎮(zhèn)居民生活需水量和農村居民生活需水量。這樣的分類便于數據統(tǒng)計也符合當地的實際情況，具有可操作性。

（6）水環(huán)境模塊。

水環(huán)境模塊中通過模型模擬成都市COD入河總量，與計算得到的水功能區(qū)COD水環(huán)境容量做比，分析水環(huán)境的納污狀況。其中城鎮(zhèn)居民生活和工業(yè)排放的COD作為點源入河，農村居民生活、農田徑流和禽畜養(yǎng)殖排放的COD作為非點源入河。[HJ1.95mm]

（7）水生態(tài)模塊。

水生態(tài)模塊通過構造水環(huán)境承載限制系數和水資源承載限制系數，最后取兩者的最大值來作為水生態(tài)承載限制系數。其中水環(huán)境承載限制系數由COD入河總量和COD水環(huán)境容量決定；水資源承載限制系數由總需水量、水資源總量、水資源開發(fā)利用率及污水回用量決定。

3.2.2 模型主要變量參數及方程式

根據前人對水生態(tài)承載力的定義和兩條研究主線，本模型主要通過三個系數：水環(huán)境承載限制系數、水資源承載限制系數和水生態(tài)承載限制系數來量化水環(huán)境、水資源和水生態(tài)的超載狀態(tài)。其中，水環(huán)境承載限制系數作為一個比值變量，首先計算COD入河總量與COD水環(huán)境容量的差值，然后用差值與COD水環(huán)境容量的比值來描述水環(huán)境的承載能力，當其比值大于0時，說明水環(huán)境已處于超載狀態(tài)；水資源承載限制系數也為比值變量，首先計算總需水量與可供水量的差值，然后用差值與可供水量的比值來描述水資源的承載能力，當其值大于0時，說明水資源已處于超載狀態(tài)；最后綜合前兩個系數的最大值，得水生態(tài)承載限制系數來描述水生態(tài)承載力超載狀態(tài)。承載限制系數的值越大，其反映的水環(huán)境、水資源和水生態(tài)的超載狀態(tài)越嚴重。主要變量參數及方程式見表1。

3.3 系統(tǒng)流圖的構建

通過分析各子系統(tǒng)之間的聯(lián)系和系統(tǒng)中各因素之間的反饋關系，本文運用系統(tǒng)動力學專用建模軟件Vensim建立了以成都市水生態(tài)承載力系統(tǒng)耦合關系圖（圖1）為基礎的成都市水生態(tài)承載力SD模型。模型系統(tǒng)總流圖見圖2。

3.4 模型參數的選取

系統(tǒng)動力學模型的行為模式與結果主要取決于模型結構，即反饋模型的行為對參數的變化是不敏感的，模型對參數的準確度滿足模型要求即可[29]。本文選取的主要參數變量的數值均取自“材料與方法—數據來源”中的數據。其初始值見表1。

4 成都市水生態(tài)承載力的有效性檢驗

本文對成都市水生態(tài)承載力SD模型進行有效性檢驗采用的是歷史性檢驗[29]，將2010年-2013年的相關數據輸入模型，將模擬的仿真結果與歷史數據對比，驗證其吻合度。

由于模型結構較為復雜，參數較多，本文僅對2010年-2013年的人口數、總需水量、工業(yè)需水量、COD入河總量的預測結果進行歷史驗證，結果見表3。

由誤差值分析可知，除2011年總需水量和2010年COD入河總量模擬值誤差略大于5%，其余模擬值誤差都低于5%，模擬結果與歷史數據基本吻合，因此認為模型模擬結果與實際數據擬合較好，模型結構合理，可信度較高。

5 成都市水生態(tài)承載力模型情景設計與分析

情景分析法[30-31]是一種以某種趨勢和現象將一直持續(xù)為前提，通過定性和定量的分析，對事件可能出現的情況作出預測的方法。在研究成都市水生態(tài)承載力的過程中，運用情景分析法提出不同的情景方案，在不同的情景模式下，對水生態(tài)承載力進行預測分析，選擇出最為有效緩解承載力的情景方案。

5.1 決策變量的選取

參考相關研究[32-35]，選取萬元工業(yè)增加值用水量、工業(yè)用水重復利用率、污水處理率、污水回用率、灌溉水利用系數、禽畜集約化養(yǎng)殖比例、城鎮(zhèn)居民生活用水定額、農村居民生活用水定額、工業(yè)污水直排COD平均濃度、污水處理廠COD平均出水濃度、城鎮(zhèn)居民生活人均COD產污系數、農村居民生活人均COD產污系數、農田單位面積COD排放量作為決策變量。

為綜合反映成都市水生態(tài)承載力狀況，選取總需水量、COD入河總量、水環(huán)境承載限制系數、水資源承載限制系數、水生態(tài)承載限制系數作為評價模型的主要指標。

5.2 情景設計

（1）現狀延續(xù)型。假定從2013年至2020年，工業(yè)用水重復利用率、污水處理率、污水回用率、禽畜集約化養(yǎng)殖比例、灌溉水利用系數、城市居民生活用水定額和農村居民生活用水定額按照現在的增長趨勢線性增加，即每年分別增加1%、0.5%、2%、2%、0.005、12 L/（人·d）、3 L/（人·d），則到2020年底，工業(yè)用水重復利用率達到90%、污水處理率達到85%、污水回用率達到20%、禽畜集約化養(yǎng)殖比例達到60%、灌溉水利用系數提高到0.52、城市居民生活用水定額提高至450 L/（人·d）、農村居民生活用水定額提高至180 L/（人·d）。根據相關工作方案和最嚴水資源管理制度[36-37]，到2020年，工業(yè)用水方面逐漸降低萬元工業(yè)增加值用水量降至2013年的70%。

（2）節(jié)約水量型。為保證水資源的可持續(xù)利用，大力推行節(jié)約用水、循環(huán)用水和污水回用。根據相關規(guī)劃、工作方案及最嚴水資源管理制度[36-39]，到2017年，全市農業(yè)灌溉水有效利用系數達到0.535；到2020年，農田灌溉水有效利用系數達到0.56以上的要求來制定節(jié)約用水方案。即從2014年開始到2020年底，居民生活用水定額方面逐步消減至2013年的80%，工業(yè)用水方面逐漸降低萬元工業(yè)增加值用水量至2013年的50%，工業(yè)用水重復利用率逐漸提高至95%，灌溉水利用系數提高至0.56，污水回用率逐步提高至40%。

（3）污染防治型。為保證水環(huán)境的可持續(xù)發(fā)展，提高水環(huán)境對外界干擾的抵抗力和水體的自凈能力，通過降低產污系數和提高污水處理率來提高水環(huán)境質量。根據相關工作方案及國家標準[21，36]，2020年底前，新建城鎮(zhèn)污水處理設施執(zhí)行一級A排放標準，中心城區(qū)、郊區(qū)（市）縣城市、鄉(xiāng)鎮(zhèn)污水集中處理率分別達到98%、88%、75%以上，推進農業(yè)農村污染防治和控制農業(yè)面源污染的基礎上來制定詳細的污染防治方案。即工業(yè)污水直排COD平均濃度從2014年起由150 mg/L降為100 mg/L，污水處理廠COD平均出水濃度由60 mg/L降為50 mg/L。從2014年開始到2020年底，居民生活COD人均產污系數和農村居民COD人均產污系數逐年消減至基準年的80%，農田單位面積COD排放量逐年消減至基準年的70%，污水處理率逐年提高到90%，禽畜集約化養(yǎng)殖比例逐年提高到80%。

（4）綜合協(xié)調型。綜合協(xié)調型即在同時考慮節(jié)約用水和污染防治的基礎上進行情景模擬。具體方案為：工業(yè)污水直排COD平均濃度從2014年起由150 mg/L降為100 mg/L，污水處理廠COD平均出水濃度由60 mg/L降為50 mg/L。從2014年開始到2020年底，居民生活用水定額方面逐步消減至2013年的80%，工業(yè)用水方面逐漸降低萬元工業(yè)增加值用水量至2013年的50%，工業(yè)用水重復利用率逐漸提高至95%，灌溉水利用系數提高至0.56，污水回用率逐步提高至40%。居民生活COD人均產污系數和農村居民COD人均產污系數逐年消減至基準年的80%，農田單位面積COD排放量逐年消減至基準年的70%，污水處理率逐年提高到90%，禽畜集約化養(yǎng)殖比例逐年提高到80%。

5.3 結果分析與討論

按照情景設計進行參數調整，運行成都市水生態(tài)承載力SD模型，模擬結果見圖3-圖7。

（1）總需水量模擬結果（圖3）中，現狀延續(xù)型和污染防治型總需水量逐年上升，綜合協(xié)調型和節(jié)約用水型總需水量逐年下降，且現狀延續(xù)型與污染防治型上升速度一致，綜合協(xié)調型與節(jié)約用水型下降速度一致。節(jié)約用水型的情景模擬，2020年總需水量降低到48.39 億m3，與2013年相比下降了8%，說明本情景中對居民生活用水、工業(yè)用水重復率、污水回用率、灌溉水有效利用系數的調整對總需水量的影響較大。

（2）COD入河總量模擬結果（圖4）中，綜合協(xié)調型<污染防治型<節(jié)約用水型<現狀延續(xù)型，在情景模擬中，污染防治型COD入河總量與2013年相比下降了11.71%，節(jié)約用水型下降了4.12%，而將兩種情景綜合協(xié)調后，COD入河總量下降了29.04%，遠遠大于兩者的加和。以上三種情景說明，污染防治型比節(jié)約用水型對減少COD入河總量效果更好，但是，如果將兩者綜合，對減少COD的入河總量的控制將會達到最佳。

（4）水環(huán)境承載限制系數模擬結果（圖6）中，綜合協(xié)調型<污染防治型<節(jié)約用水型<現狀延續(xù)型，其中除了現狀延續(xù)型，其他三種情景的水環(huán)境承載系數自2013年以來逐年變小，說明成都市的水環(huán)境現狀在這三種情境下會逐年變好。通過曲線可以看出，節(jié)約用水型情境下，曲線的下降速率最?。痪C合協(xié)調型情境下，曲線的下降速率最大；污染防治型情景下，曲線的下降速率居中。到2020年，節(jié)約用水型情境下的水環(huán)境承載限制系數為0.53，污染防治型情境下的水環(huán)境承載限制系數為0.41，綜合協(xié)調型情境下的水環(huán)境承載限制系數為0.13，說明了只有通過節(jié)約用水和污染防治的綜合協(xié)調，才能使水環(huán)境的質量得到更好的改善。

（5）水生態(tài)承載限制系數模擬結果曲線（圖7）是水資源承載限制系數和水環(huán)境承載限制系數綜合后的一條曲線，該曲線反映出了水生態(tài)的兩個方面水資源與水環(huán)境。在該曲線中，綜合協(xié)調型<節(jié)約用水型<污染防治型<現狀延續(xù)型，同時在污染防治型和現狀延續(xù)型情景下，水生態(tài)承載限制系數逐年變大，說明節(jié)約用水才是影響水生態(tài)承載力的主要因素。比較節(jié)約用水型和綜合協(xié)調型曲線的下降速率，發(fā)現節(jié)約用水型下降速率遠遠小于綜合協(xié)調型，更進一步說明了1+1>2這種現象在科研中和生活中的普遍性，只有通過節(jié)約用水和污染防治的綜合協(xié)調，才能使水生態(tài)的質量得到更快的改善，滿足水生態(tài)可持續(xù)發(fā)展的目標。

由模擬結果可知，從2014年-2020年，成都市的水生態(tài)承載力在四種情景下由大到小的順序為：綜合協(xié)調型，節(jié)約用水型，污染防治型，現狀延續(xù)型。若成都市按當前的現狀發(fā)展，其水生態(tài)承載力將越來越弱，最終可能造成生態(tài)嚴重破壞和經濟巨大損失；當單獨考慮污染防治或者節(jié)約用水時，得到的水生態(tài)承載限制年變化曲線減小速率較慢；若同時考慮污染防治和節(jié)約用水，就可以得到減小速率較快的水生態(tài)承載限制系數年變化曲線，到2020年的水生態(tài)承載限制系數與節(jié)約用水型情景相比，將降至節(jié)約用水型的一半。因此，污染防治與節(jié)約用水相結合，是提高成都市水生態(tài)承載力的關鍵。

6 結論

根據系統(tǒng)動力學方法建立了成都市水生態(tài)承載力模型，以圖形的形式直觀的反映出4種不同情景下的總需水量、COD入河總量、水環(huán)境承載限制系數、水資源承載限制系數、水生態(tài)承載限制系數隨時間的變化趨勢。模擬結果可為成都市的社會經濟和資源環(huán)境協(xié)調發(fā)展提供較為有價值的技術依據。

與前人在水生態(tài)承載力方面的研究相比，本文給出了影響水生態(tài)承載力主要因素的計算結果，同時提出的水生態(tài)承載系數可更加綜合、精確地反映承載力的大小，使不同地區(qū)的水生態(tài)承載力具有了可比性，有利于模型的推廣和使用。建議在以后的研究中考慮水環(huán)境容量的歷年變化情況，以期獲得更為準確、客觀的成果。

參考文獻（References）：

[1] 陳國衛(wèi)，金家善，耿俊豹.系統(tǒng)動力學應用研究綜述[J].控制工程，2012（6）：921-928.（CHEN Guo-wei，JIN Jia-shan，GENG Jun-bao.Application research overview of system dynamics [J].Control Engineering of China，2012，（6）：921-928.（in Chinese））DOI：10.3969/j.issn.1671-7848.2012.06.001

[2] 鐘永光，賈曉箐，錢穎，等.系統(tǒng)動力學[M].北京：科學出版社，2013.（ZHONG Yong-guang，JIA Xiao-qing，QIAN Ying，et al.System dynamics [M].Beijing：Science Press，2013.（in Chinese））

[3] 柴淼瑞.基于SD模型的流域水生態(tài)承載力研究—以鐵嶺控制單元為例[D].西安：西安建筑科技大學，2014.（CHAI Miao-rui.Research of carrying capacity of water ecological based on the sd model：illustrated by the tieling control unit [D].Xi′an：Xi′an University of Architecture and Technology，2014.（in Chinese））

[4] 李靖，周孝德.葉爾羌河流域水生態(tài)承載力研究[J].西安理工大學學報，2009（3）：249-255.（LI Jing，ZHOU Xiao-de.Research on water ecological carrying capacity of Yarkant river basin [J].Journal of Xi′an University of Technology，2009 （3）：249-255.（in Chinese））DOI：10.3969/j.issn.1006-4710.2009.03.001

[5] 彭文啟.流域水生態(tài)承載力理論與優(yōu)化調控模型方法[J].中國工程科學，2013（3）：33-43.（PENG Wen-qi.Watershed water ecological carrying capacity and optimizing regulation model [J].Engineering Science in China，2013（3）：33-43.（in Chinese））DOI：10.3969/j.issn.1009-1742.2013.03.007

[6] 薛小妮，甘泓，游進軍，等.成都市水資源及水環(huán)境承載能力分析[J].水利水電技術，2012（4）：14-18.（XUE Xiao-ni，GAN Hong，YOU Jin-jun，et al.Study on water resources and water environmental carrying capacities of Chengdu [J].Water Resources and Hydropower Engineering，2012 （4）：14-18.（in Chinese））DOI：10.3969/j.issn.1000-0860.2012.04.004

[7] 陳琳，鄒添丞，石杰，等.基于層次分析法的成都市水安全評價[J].南水北調與水利科技，2013，11（4）：41-45.（CHEN Lin，ZOU Tian-cheng，SHI Jie，et al.Water security assessment of Chengdu City based on analytic hierarchy process [J].South-to-North Water Transfers and Water Science & Technology，2013，11 （4）：41-45.（in Chinese））DOI：10.3724/SP.J.1201.2013.04041

[8] 歐洋銘，白若男，朱國宇.成都市水資源承載力動態(tài)變化及成因[J].南水北調與水利科技，2012，10（6）：110-114.（OU Yang-ming，BAI Ruo-nan，ZHU Guo-yu.Dynamic variation of water resources carrying capacity and its causes in Chengdu [J].South-to-North Water Transfers and Water Science & Technology，2012，10（6）：110-114.（in Chinese））DOI：10.3724/SP.J.1201.2012.06110

[9] 蘇敏.基于主成分分析法的水資源承載力影響因素分析——以成都市為例[J].科技與生活，2012（19）：218-219.（SU Min.Analysis of water resources carrying capacity based on principal component analysis of factors：a case study of Chengdu City [J].Technology and Living，2012（19）：218-219.（in Chinese））

[10] 成都市地方志編纂委員會.成都年鑒[Z].成都：成都年鑒社，2006-2013.（Chengdu Local Chronicles Compilation Committee.Chengdu yearbook [Z].Chengdu：Chengdu Yearbook Press，2006-2013.（in Chinese））

[11] 成都市統(tǒng)計局，國家統(tǒng)計局成都調查隊.2015年成都市國民經濟和社會發(fā)展統(tǒng)計公報[Z].成都：成都市統(tǒng)計局，2016.（Statistics Bureau of Chengdu，NBS Survey Office in Chengdu.Statistical bulletin of the national economic and social development of Chengdu city in 2015 [Z].Chengdu：Statistics Bureau of Chengdu，2016.（in Chinese））

[12] 成都市水務局.成都市水資源公報[Z].成都：成都市水務局，2004-2013.（Chengdu Water Authority.Water resources bulletin of Chengdu [R].Chengdu：Chengdu Water Authority，2004-2013.（in Chinese））

[13] 成都市環(huán)保局.2015年成都市總量減排工作計劃[Z].成都：成都市環(huán)保局，2015.（Environment Protection Bureau of Chengdu.Total emission reduction work plan of Chengdu City in 2015[Z].Chengdu：Environment Protection Bureau of Chengdu，2015.（in Chinese））

[14] 蔣偉，柴夏，顏飛.成都地區(qū)水質調查與富營養(yǎng)化評價[J].環(huán)境科學與技術，2012，35（6I）：370-372.（JIANG Wei，CHAI Xia，YAN Fei.Water quality surveys in Chengdu and eutrophication [J].Environmental Science & Technology，2012，35（6I）：370-372.（in Chinese））

[15] 崔寧，梁冬梅，蘇偉.水環(huán)境承載力量化方法述評[J].吉林水利，2009（11）：45-47.（CUI Ning，LIANG Dong-mei，SU Wei.Review on the methods of water environment carrying capacity [J].Jilin Water Resources，2009（11）：45-47.（in Chinese））

[16] 趙元慧.鐵嶺市水環(huán)境承載力研究[D].沈陽：沈陽理工大學，2012.（ZHAO Yuan-hui.Research on water environment carrying capacity of Tieling [D].Shenyang：Shenyang Ligong University，2012.（in Chinese））

[17] 四川省統(tǒng)計局，國家統(tǒng)計局四川調查總隊.四川統(tǒng)計年鑒 [Z].北京：中國統(tǒng)計出版社，2004-2014.（Statistics Bureau of Sichuan，NBS Survey Office in Sichuan.Sichuan statistical yearbook [Z].Beijing：China Statistics Press，2004-2014.（in Chinese））

[18] 成都市統(tǒng)計局，國家統(tǒng)計局成都調查隊，成都市統(tǒng)計學會.成都統(tǒng)計年鑒[Z].北京：中國統(tǒng)計出版社，2003-2013.（Statistics Bureau of Chengdu，NBS Survey Office in Chengdu，Statistical Society of Chengdu.Chengdu statistical yearbook [Z].Beijing：China Statistics Press，2003-2013.（in Chinese））

[19] 四川省水利科學研究院.四川省用水定額（修訂稿）[Z].成都：四川省水利廳，2010.（Sichuan Academy of Water Conservancy.Water quota of Sichuan Province （revised draft）[Z].Chengdu：Sichuan Provincial Water Resources Department，2010.（in Chinese））

[20] 環(huán)境保護部環(huán)境規(guī)劃院.全國水環(huán)境容量核定技術指南[Z].北京：環(huán)境保護部環(huán)境規(guī)劃院，2003.（Chinese Academy for Environmental Planning.National water environmental capacity verification manual [Z].Beijing：Chinese Academy for Environmental Planning，2003.（in Chinese））

[21] GB8978-1996，污水綜合排放標準[S].（GB8978-1996，Integrated standard for waste water discharge [S].（in Chinese））

[22] 錢俊，楊姍姍，肖杰，等.四川省生活污水產污系數研究[J].四川環(huán)境，2009（1）：27-32.（QIAN Jun，YANG Shan-shan，XIAO Jie，et al.Study on pollutant rate in Sichuan Province [J].Sichuan Environment，2009 （1）：27-32.（in Chinese））DOI：10.3969/j.issn.1001-3644.2009.01.009

[23] 楊立夢，付永勝，高紅濤.四川省茫溪河流域非點源污染負荷研究[J].重慶理工大學學報：自然科學，2014（11）：57-63.（YANG Li-meng，FU Yong-sheng，GAO Hong-tao.Study on non-point source pollution of Mangxi River Basin in Sichuan Province [J].Journal of Chongqing University of Technology：Natural Science，2014（11）：57-63.（in Chinese））DOI：10.3969/j.issn.1674-8425（z）.2014.11.012

[24] 馮猛，邱忠平，白磊磊，等.四川省規(guī)?；菪箴B(yǎng)殖污染防治現狀[J].畜牧與飼料科學，2010（3）：107-108.（FENG Meng，QIU Zhong-ping，BAI Lei-lei，et al.Pollution prevention and control of large-scale livestock and poultry production in current Sichuan Province [J].Animal Husbandry and Feed Science，2010 （3）：107-108.（in Chinese））DOI：10.3969/j.issn.1672-5190.2010.03.047

[25] 方國華，于鳳存，曹永瀟.中國水環(huán)境容量研究概述[J].安徽農業(yè)科學，2007（27）：8601-8602.（FANG Guo-hua，YU Feng-cun，CAO Yong-xiao.Review of water environment capacity in China [J].Journal of Anhui Agricultural Sciences，2007 （27）：8601-8602.（in Chinese））DOI：10.3969/j.issn.0517-6611.2007.27.090

[26] 劉曉臣，李小平，陳小華，等.SD模型在洱海流域營養(yǎng)物消減策略研究中的應用[J].湖泊科學，2013（5）：655-664.（LIU Xiao-chen，LI Xiao-ping，CHEN Xiao-hua，et al.Reduction strategies for TN and TP based on system dynamics model （SD）in Lake Erhai [J].Journal of Lake Sciences，2013（5）：655-664.（in Chinese））DOI：10.18307/2013.0506

[27] 鄭富玲.銀川市再生水回用SD模型的構建[D].銀川：寧夏大學，2015.（ZHENG Fu-ling.The construction of SD model for reclaimed water reuse in Yinchuan City [D].Yinchuan：Ningxia University，2015.（in Chinese））

[28] 曾肇京，石海峰.中國水資源利用發(fā)展趨勢合理性分析[J].中國水利，2000（8）：45-48.（ZENG Zhao-jing，SHI Hai-feng.Rationality analysis of development and utilization trend of water resources in China [J].China Water Resources，2000（8）：45-48.（in Chinese））

[29] 肖悅.基于SD模型的區(qū)域水資源承載力研究—以黃石港工業(yè)園區(qū)為例[D].武漢：武漢科技大學，2012.（XIAO Yue.Regional water resources carrying capacity based on SD model of Huangshi Port Industrial Park [D].Wuhan：Wuhan University of Science and Technology.（in Chinese））

[30] 劉俏.情景分析法在城市規(guī)劃區(qū)域污染物排放總量中的預測研究-以安慶市為例[D].合肥：合肥工業(yè)大學，2013.（LIU Qiao.Research on the prediction of the total discharge of pollutants of urban planning based on scenario analysis：the case study of Anqing City [D].Hefei：Hefei University of Technology，2013.（in Chinese））

[31] 李瑋，秦大庸，褚俊英，等.基于情景分析法的污染物排放趨勢研究[J].水電能源科學，2010（5）：36-39.（LI Wei，QIN Da-yong，ZHU Jun-ying，et al.Study on pollutant emission trend based on scenario analysis [J].Water Resources and Power，2010（5）：36-39.（in Chinese））DOI：10.3969/j.issn.1000-7709.2010.05.011

[32] 羅宇，姚幫松.基于SD模型的長沙市水資源承載力研究[J].中國農村水利水電，2015（1）：42-46.（LUO Yu，YAO Bang-song.The water resources carrying capacity based on SD model in Changsha [J].China Rural Water and Hydropower，2015（1）：42-46.（in Chinese））DOI：10.3969/j.issn.1007-2284.2015.01.012

[33] 馮海燕，張昕，李光永，等.北京市水資源承載力系統(tǒng)動力學模擬[J].中國農業(yè)大學學報，2006（6）：106-110.（FENG Hai-yan，ZHANG Xin，LI Guang-yong，et al.A system dynamic model and simulation for water resources carrying capacity in Beijing [J].Journal of China Agricultural University，2006（6）：106-110.（in Chinese））DOI：10.3321/j.issn：1007-4333.2006.06.021

[34] 陳雪.SD模型在阿什河流域水污染控制規(guī)劃中的應用[D].哈爾濱：東北林業(yè)大學，2010.（CHEN Xue.The appliance of SD model in the pollution control planning of Ashi River Basin [D].Harbin：Northeast Forestry University，2010.（in Chinese））

[35] 劉曉玉.水生態(tài)承載力政策設計與政策效率的系統(tǒng)仿真—以贛江流域水生態(tài)承載力為例[D].南昌：南昌大學，2011.（LIU Xiao-yu.System simulation of the policy design and efficiency for the water ecology carrying capacity：taking the water ecology carrying capacity of Gan River Basin for example [D].Nanchang：Nanchang University，2011.（in Chinese））

[36] 成都市人民政府.成都市水污染防治工作方案[Z].成都：成都市人民政府，2016.（People′s Government of Chengdu.Scheme of water pollution prevention and control in Chengdu City [Z].Chengdu：People′s Government of Chengdu，2016.（in Chinese））

[37] 成都市人民政府.成都市人民政府辦公廳關于實行最嚴格水資源管理制度考核方法的通知[Z].成都：成都市人民政府，2014.（People′s Government of Chengdu.Notice to implement the most stringent water management system assessment in Chengdu City [Z].Chengdu：People′s Government of Chengdu，2014.（in Chinese））

[38] 成都市人民政府.成都市水務發(fā)展“十二五”規(guī)劃[Z].成都：成都市人民政府，2012.（People′s Government of Chengdu.The 12th five- year plan of water development in Chengdu City [Z].Chengdu：People′s Government of Chengdu，2012.（in Chinese））

[39] 成都市人民政府.成都市循環(huán)經濟發(fā)展“十二五”規(guī)劃[Z].成都：成都市人民政府，2012.（People′s Government of Chengdu.The 12th five- year plan of circular economy development in Chengdu City [Z].Chengdu：People′s Government of Chengdu，2012.（in Chinese））