王宵君,楊 歡,李楨萍,嚴(yán)文鳳,藍四嬌,杜 超
(贛南師范大學(xué)地理與環(huán)境工程學(xué)院,江西 贛州 341000)
水資源承載力(Water Resources Carrying Capacity)是以可預(yù)見的人口、環(huán)境、經(jīng)濟發(fā)展水平為依據(jù),以尋求環(huán)境經(jīng)濟平衡發(fā)展為前提,經(jīng)過資源的合理優(yōu)化配置,對該地區(qū)社會經(jīng)濟與人口規(guī)模發(fā)展所能提供的最大支撐能力預(yù)測[1-2]。因此,水資源承載力研究涉及一個復(fù)雜巨系統(tǒng),包括社會經(jīng)濟系統(tǒng)、水資源系統(tǒng)、水環(huán)境系統(tǒng)[3]。水資源承載力的可持續(xù)發(fā)展不僅涉及生態(tài)安全問題,也關(guān)系到一個區(qū)域的經(jīng)濟發(fā)展與社會穩(wěn)定[2]。目前,水資源承載力與區(qū)域社會經(jīng)濟發(fā)展的耦合研究已成為當(dāng)今地理科學(xué)與區(qū)域經(jīng)濟學(xué)的一個熱點話題。
目前,關(guān)于水資源承載力的研究,已有較多報道:①在研究內(nèi)容上,多集中于水資源承載力定義、水資源承載力影響因素、水資源承載力評價指標(biāo)體系確立、水資源承載力模型構(gòu)建等方面;②研究尺度上,目前的研究多集中在流域、區(qū)域、城市3個不同的空間尺度[3],而且這些區(qū)域多是缺水或者生態(tài)環(huán)境脆弱地區(qū),對于水資源相對豐富的區(qū)域研究較少;③研究方法上,代表性研究方法主要有常規(guī)趨勢法[4-5]、模糊綜合評價法[6]、多目標(biāo)分析法[7]、主成分分析法[8-9]及系統(tǒng)動力學(xué)法等[10-11]。其中,系統(tǒng)動力學(xué)(System Dynamic,SD)方法是以計算機仿真為框架,研究復(fù)雜系統(tǒng)內(nèi)部的動態(tài)演化行為的工具[12],它在復(fù)雜的非線性系統(tǒng)研究中具有不可替代的優(yōu)勢,目前已經(jīng)廣泛應(yīng)用于工業(yè)、農(nóng)業(yè)、生態(tài)、環(huán)境等諸多領(lǐng)域,在水資源系統(tǒng)中的應(yīng)用也日益凸顯其優(yōu)勢[13-15]。利用系統(tǒng)動力學(xué)模型能夠綜合影響水資源承載力的人口、經(jīng)濟、社會、生態(tài)等因素,充分刻畫影響水資源承載力的演化機理,突出經(jīng)濟發(fā)展、資源與人口的耦合關(guān)系,凸顯了水資源承載力評估管理的重要性[2]。贛州市作為江西省唯一的省域副中心城市,擁有九百多萬人口,地處江西省南部,屬亞熱帶季風(fēng)性濕潤氣候,基本屬于富水區(qū),但近年來,隨著贛州社會經(jīng)濟迅速發(fā)展,水資源消耗日趨增加,且因工業(yè)發(fā)展、稀土開采等造成水資源污染日趨嚴(yán)重。另一方面,贛州作為贛江、東江和北江三江源頭所在地區(qū),保持它的水資源、水環(huán)境安全對于長江經(jīng)濟帶和粵港澳大灣區(qū)建設(shè)意義重大。因此,科學(xué)評估贛州市水資源承載力具有重要的研究意義。目前學(xué)者們對贛州市水資源承載力的研究主要是用模糊綜合評價法[16]、因子分析和熵權(quán)法[17],而采用系統(tǒng)動力學(xué)法研究其水資源承載力可持續(xù)發(fā)展目前還沒有相關(guān)報道。因此,本文采用系統(tǒng)動力學(xué)方法,通過構(gòu)建贛州市水資源承載力可持續(xù)發(fā)展SD模型,對不同情景下水資源承載力進行模擬仿真,探索贛州市水資源承載力可持續(xù)發(fā)展的最佳模式,為緩解贛州市水資源矛盾和建立統(tǒng)一的水資源管理體制提供理論參考與現(xiàn)實依據(jù),以促進贛州市的社會發(fā)展。
贛州市(24°29′~27°09′N、113°54′~116°38′E)位于江西省南部(圖1),贛江上游,簡稱贛南,是江西省最大的行政區(qū)。贛州市屬亞熱帶季風(fēng)性濕潤氣候,基本屬于富水區(qū)。贛州是江西的母親河——贛江的源頭,也是廣東、香港同胞飲用水水源東江、北江的源頭,承擔(dān)著“一江清水給下游”的重任。
近年來,贛州市的社會經(jīng)濟得到飛速發(fā)展。根據(jù)贛州市國民經(jīng)濟和社會發(fā)展統(tǒng)計公報,2018年,贛州市地區(qū)生產(chǎn)總值(GDP)2 807.24億元,比2017年增長9.3%。其中,第一產(chǎn)業(yè)增加值340.30億元,增長3.7%;第二產(chǎn)業(yè)增加值1 194.24億元,增長8.9%;第三產(chǎn)業(yè)增加值1 272.70億元,增長11.5%。但是,隨著經(jīng)濟和工業(yè)的發(fā)展,以及水資源自身條件的變化,加上現(xiàn)有水利設(shè)施和水利管理體制的局限性,贛州市在水資源、水環(huán)境等方面都存在一些問題[18]。首先,近年來隨著社會經(jīng)濟發(fā)展速度提升,贛州市廢水排放量增多,水環(huán)境問題日益凸顯。贛州市水資源公報顯示,2018年,贛州廢污水排放量為54 300萬t,其中城鎮(zhèn)居民生活廢水排放量達18 900萬t,占總排放量34.8%;第二產(chǎn)業(yè)廢水29 300萬t,占總排放量54.0%;第三產(chǎn)業(yè)廢水6 100萬t,占總排放量11.2%。其次,水資源問題也更加突出。根據(jù)贛州市水資源公報顯示,2018年贛州市地表水資源量222.92億m3,比2017年減少19.1%,比多年平均減少33.8%;地下水資源量68.92億m3,比2017年減少14.8%,比多年平均減少22.4%。此外,近年來臍橙種植業(yè)的發(fā)展也使得一些地區(qū)水土流失嚴(yán)重;農(nóng)業(yè)灌溉和工業(yè)生活用水方面的水資源利用率也較低,而區(qū)域水資源的自然稟賦和調(diào)控狀況將直接影響贛州市的社會經(jīng)濟發(fā)展和科學(xué)技術(shù)水平[19]。因此,為了確保贛江、東江和北江源頭及其流域的水資源、水環(huán)境安全,促進贛州市社會經(jīng)濟的進一步發(fā)展,有必要對其水資源承載力進行研究分析。
圖1 贛州市區(qū)位
所建立的贛州市水資源承載力SD模型數(shù)據(jù)來源于以下數(shù)據(jù)集:《江西省統(tǒng)計年鑒》(2011—2016)、《江西省水資源公報》(2011—2016)、《江西省環(huán)境統(tǒng)計年報》(2011—2016)、DB 36-T 419—2017《江西省生活用水定額》《贛州市“十三五”節(jié)能減排綜合工作方案》《贛州市國民經(jīng)濟和社會發(fā)展第十四個五年規(guī)劃和二〇三五年遠景目標(biāo)綱要》《江西省節(jié)水行動實施方案》《贛州市城市總體規(guī)劃》(2017—2035年)等。
水資源承載力研究涉及一個復(fù)雜巨系統(tǒng),包括社會經(jīng)濟系統(tǒng)、水資源系統(tǒng)、水環(huán)境系統(tǒng),重在水資源承載力與區(qū)域社會經(jīng)濟發(fā)展各要素的耦合。利用系統(tǒng)動力學(xué)方法,以贛州市為研究區(qū),在Vensim? PLE軟件支持下,建立水資源承載力因果關(guān)系圖,并構(gòu)建水資源承載力SD仿真模擬模型,通過進行不同情景下的模擬分析和預(yù)測,分析水資源能夠承載的最佳可持續(xù)發(fā)展模式。具體方法流程見圖2。
圖2 研究路線
2.1.1模型的構(gòu)建
SD模型是一種從系統(tǒng)內(nèi)部關(guān)系入手的系統(tǒng)與綜合的研究方法,適合研究多變量、非線性、高階次、多回路的復(fù)雜系統(tǒng),通過分析系統(tǒng)作用機制,可全面反映事物發(fā)展過程及動因[20],是模擬預(yù)測水資源承載力動態(tài)發(fā)展的有效方法。根據(jù)贛州市的水資源水環(huán)境特征,構(gòu)建模型結(jié)構(gòu)和反饋關(guān)系(圖3)。由于研究區(qū)域?qū)儆诟凰畢^(qū)域,因此,在建立系統(tǒng)流程圖時,以反映贛州市主要因果關(guān)系,選擇影響贛州市水資源承載力問題的核心因子,最終確定參數(shù)。
SD系統(tǒng)模擬使用Vensim?PLE軟件建模,確定贛州市水資源承載力系統(tǒng)的邊界,空間邊界為贛州市,模型運行時間為2011—2035年,模擬時間步長為1年,其中,2020年為預(yù)測基準(zhǔn)年。系統(tǒng)由社會經(jīng)濟子系統(tǒng)、水資源子系統(tǒng)和水環(huán)境子系統(tǒng)組成,研究內(nèi)容主要集中于贛州市水資源供需水量和水環(huán)境污染2個方面。確定贛州市水資源承載力系統(tǒng)的主要水平變量(贛州市人口、GDP等)、速率變量(人口增長率、GDP增長率)及輔助變量(工業(yè)用水量、生活用水量、COD排放量等)。SD模型的系統(tǒng)流程見圖4,并確定模型的主要方程和反饋關(guān)系。
圖3 贛州市水資源承載力系統(tǒng)基本反饋回路
圖4 贛州市水資源承載力系統(tǒng)動力學(xué)模型
2.1.2模型的檢驗
為驗證建立的系統(tǒng)動力學(xué)模型的有效性,可通過直接將模擬值與歷史值進行誤差分析來評估模型結(jié)構(gòu)的準(zhǔn)確性。與2012—2016年的實際值進行對比,系統(tǒng)模擬結(jié)果的整體變化趨勢與歷史結(jié)果相吻合,與贛州市實際情況基本相符,變量的相對誤差基本在±10%以內(nèi),說明該模型有效。表1所示為3個子系統(tǒng)中總?cè)丝?、GDP總量、生態(tài)用水量、總供水量、COD排放總量5個關(guān)鍵性的指標(biāo)的歷史檢驗數(shù)據(jù)。對比相關(guān)研究成果[2-3,23],模型可以用于模擬預(yù)測贛州市2020—2035年水資源承載力系統(tǒng)的可持續(xù)發(fā)展情況。
表1 歷史性檢驗數(shù)據(jù) %
2.2.1情景設(shè)計
本研究中水資源承載能力最佳模式的確定將是在考慮贛州市未來發(fā)展趨勢和特征的基礎(chǔ)上,通過比較不同發(fā)展方案情景下SD模型的模擬結(jié)果,提出符合研究區(qū)水資源承載力可持續(xù)發(fā)展的最佳模式。因此,確定合適的系統(tǒng)模擬方案,并根據(jù)需要選擇模擬參數(shù)是本文的重要任務(wù)之一。
為了更好地實現(xiàn)贛州市水生態(tài)的可持續(xù)發(fā)展,參照前人研究基礎(chǔ),確定的4種模擬方案分別為正常發(fā)展型、經(jīng)濟優(yōu)先型、環(huán)境優(yōu)先型、綜合協(xié)調(diào)型。情景1:正常發(fā)展型,即現(xiàn)狀延續(xù)型。該模式根據(jù)研究區(qū)發(fā)展現(xiàn)狀和未來發(fā)展規(guī)劃設(shè)定指標(biāo)數(shù)據(jù),為延續(xù)現(xiàn)狀的發(fā)展模式。情景2:經(jīng)濟優(yōu)先型,即重點發(fā)展經(jīng)濟。突出社會經(jīng)濟發(fā)展的目標(biāo),將更多的資源投入到社會經(jīng)濟發(fā)展中,在現(xiàn)狀發(fā)展的基礎(chǔ)上,上調(diào)社會經(jīng)濟指標(biāo)數(shù)據(jù)。情景3:環(huán)境優(yōu)先型,即重點抓生態(tài)環(huán)境保護,突出環(huán)境保護目標(biāo),政策制定優(yōu)先考慮環(huán)境因素,降低經(jīng)濟發(fā)展速度和強度,污水排放和處理率降低。情景4:綜合協(xié)調(diào)型,既重視社會經(jīng)濟發(fā)展,又關(guān)注生態(tài)環(huán)境保護。具體方案設(shè)計見表2。
表2 SD模型情景設(shè)計
2.2.2參數(shù)的設(shè)定
在系統(tǒng)動力學(xué)模型中,通過調(diào)整相關(guān)的參數(shù),可以模擬承載力的提升以及其他年份的預(yù)測。在前人[21-22]研究基礎(chǔ)上,參考《江西省節(jié)水行動實施方案》《贛州市國民經(jīng)濟和社會發(fā)展第十四個五年規(guī)劃和二〇三五年遠景目標(biāo)綱要》《贛州市2020年國民經(jīng)濟和社會發(fā)展計劃》《贛州市水資源公報》《贛州市水污染防治工作計劃》等重要文件,設(shè)置各參考方案的參數(shù)。
a)情景1:正常發(fā)展型。該方案為現(xiàn)狀延續(xù)型,以近年來贛州市資源、環(huán)境、社會和經(jīng)濟等各方面的發(fā)展現(xiàn)狀和近遠期發(fā)展規(guī)劃為基礎(chǔ),將人口增長率定為4.8‰,城鎮(zhèn)化率定為0.53%,GDP增長率定為7.5%,萬元工業(yè)增加值用水量設(shè)置為41萬m3,工業(yè)污水排放系數(shù)設(shè)為0.27,生活污水排放系數(shù)設(shè)為0.49。
b)情景2:經(jīng)濟優(yōu)先型。以經(jīng)濟發(fā)展為重點,即根據(jù)相關(guān)規(guī)劃中較高水平調(diào)控參數(shù)設(shè)置為較高數(shù)值,均在情景1的基礎(chǔ)上,上調(diào)20%。
c)情景3:環(huán)境優(yōu)先型。該方案以環(huán)境保護為主,即調(diào)控參數(shù)設(shè)置為較低數(shù)值,在情景1的基礎(chǔ)上除城鎮(zhèn)化率響應(yīng)國家新型城鎮(zhèn)化戰(zhàn)略保持不變,其他參數(shù)均下調(diào)20%。
d)情景4:綜合協(xié)調(diào)型。調(diào)控參數(shù)中社會經(jīng)濟數(shù)據(jù)、水環(huán)境、水資源等參數(shù)均取最佳值,其中社會經(jīng)濟參數(shù)在情景1基礎(chǔ)上上調(diào)10%,水環(huán)境、水資源參數(shù)與環(huán)境優(yōu)先型保持一致。具體模擬方案主要參數(shù)設(shè)定見表3。
表3 方案及調(diào)控參數(shù)設(shè)定
將SD模型中的4種方案中的相關(guān)參數(shù)按照上述步驟進行設(shè)置,重新運行系統(tǒng)動力學(xué)模型,得到不同方案下的模擬數(shù)據(jù)。本文所建立的系統(tǒng)動力學(xué)模型中,GDP、總?cè)丝跀?shù)量可作為社會經(jīng)濟評價數(shù)據(jù);工業(yè)污水COD排放量、COD排放總量作為水環(huán)境評價數(shù)據(jù);工業(yè)用水量和生產(chǎn)用水量作為水資源評價數(shù)據(jù)。以上這些數(shù)據(jù)均可作為評價方案優(yōu)越性的主要參考數(shù)據(jù),2025、2035年主要變量仿真結(jié)果見表4。
表4 主要變量仿真結(jié)果
4種情景下,2020—2035年社會經(jīng)濟發(fā)展情況為:經(jīng)濟優(yōu)先型>綜合協(xié)調(diào)型>正常發(fā)展型>環(huán)境優(yōu)先型。在國內(nèi)生產(chǎn)總值方面,其中經(jīng)濟優(yōu)先型GDP增幅最高,2025、2035年,GDP總量分別達5 608.59億、13 277.60億元;環(huán)境優(yōu)先型的GDP增幅最低,將分別達到4 878.10億、8 735.93億元;綜合協(xié)調(diào)型增幅高于正常發(fā)展型,低于經(jīng)濟優(yōu)先型,2025、2035年,綜合協(xié)調(diào)型的GDP總量分別達5 418.27億、11 971.30億元。4種方案GDP動態(tài)仿真結(jié)果見圖5。在人口增長方面,經(jīng)濟優(yōu)先型呈現(xiàn)人口高速增長,至2025、2035年,總?cè)丝趯⒎謩e達到923.13萬、977.71萬人,環(huán)境優(yōu)先型呈現(xiàn)人口低速增長,至2025、2035年,總?cè)丝趯⒎謩e達到914.36萬、950.08萬人。而綜合協(xié)調(diào)型人口增幅雖低于經(jīng)濟優(yōu)先型,但高于正常發(fā)展型,至2025、2035年,總?cè)丝趯⒎謩e達到920.93萬、970.73萬人。4種方案總?cè)丝跀?shù)量動態(tài)仿真結(jié)果見圖6。
圖5 GDP動態(tài)仿真結(jié)果
圖6 總?cè)丝跀?shù)量動態(tài)仿真結(jié)果
4種情景下,2020—2035年水資源情況為:環(huán)境優(yōu)先型>綜合協(xié)調(diào)型>正常發(fā)展型>經(jīng)濟優(yōu)先型。在用水量方面,其中經(jīng)濟優(yōu)先型用水量增幅最高,2025、2035年,其工業(yè)用水量分別達9.71億、19.11億m3,生產(chǎn)用水量分別達到34.31億、44.48億m3;環(huán)境優(yōu)先型的用水量增幅最低,2025、2035年,其工業(yè)用水量將分別達到5.63億、8.38億m3,生產(chǎn)用水量分別達到30.23億、33.76億m3;綜合協(xié)調(diào)型用水量增幅略高于環(huán)境優(yōu)先型,但低于正常發(fā)展型,2025和2035年,綜合協(xié)調(diào)型的工業(yè)用水量分別達6.25億、11.49億m3,生產(chǎn)用水量分別達到30.86億、36.86億m3。4種方案工業(yè)用水量動態(tài)仿真結(jié)果見圖7,生產(chǎn)用水量動態(tài)仿真結(jié)果見圖8。
圖7 工業(yè)用水量動態(tài)仿真結(jié)果
圖8 生產(chǎn)用水量動態(tài)仿真結(jié)果
4種情景下,2020—2035年水環(huán)境情況為:環(huán)境優(yōu)先型>綜合協(xié)調(diào)型>正常發(fā)展型>經(jīng)濟優(yōu)先型。在COD排放總量方面,其中經(jīng)濟優(yōu)先型COD排放總量增幅最高,2025、2035年,其COD排放量分別達17.97萬、25.17萬t;環(huán)境優(yōu)先型的COD排放總量增幅最低,2025、2035年,其COD排放總量將分別達到9.62萬、11.21萬t;綜合協(xié)調(diào)型COD排放總量增幅略高于環(huán)境優(yōu)先型,但低于正常發(fā)展型,2025、2035年,綜合協(xié)調(diào)型的COD排放總量分別達10.39萬、13.29萬t。4種方案COD排放總量動態(tài)仿真結(jié)果見圖9。
圖9 COD排放總量動態(tài)仿真結(jié)果
綜上所述,在4種模擬方案中,從2020—2035年,社會經(jīng)濟發(fā)展趨勢來看,經(jīng)濟優(yōu)先型>綜合協(xié)調(diào)型>正常發(fā)展型>環(huán)境優(yōu)先型;從水資源、水環(huán)境發(fā)展情況趨勢來看,環(huán)境優(yōu)先型>綜合協(xié)調(diào)型>正常發(fā)展型>經(jīng)濟優(yōu)先型。這就表明,在延續(xù)現(xiàn)狀發(fā)展的正常發(fā)展模式和經(jīng)濟優(yōu)先發(fā)展模式下,贛州市社會經(jīng)濟高速增長時若不采取嚴(yán)格的節(jié)水、治理污水措施,工業(yè)用水量、污水排放量將大幅增加,違背我國資源節(jié)約、環(huán)境友好發(fā)展戰(zhàn)略。而在其環(huán)境優(yōu)先型發(fā)展模式下,經(jīng)濟規(guī)模增加趨勢較小,無法達到社會經(jīng)濟發(fā)展要求。根據(jù)系統(tǒng)預(yù)測結(jié)果,在綜合協(xié)調(diào)型模式下,至2035年,COD排放量保持較低水平,為13.29萬t;而經(jīng)濟發(fā)展處于相對較高水平,其中GDP達到11 971.3億元。即只有在綜合協(xié)調(diào)發(fā)展模式下,既保證社會經(jīng)濟維持中高速增長態(tài)勢,符合社會經(jīng)濟發(fā)展要求,又使得水資源水生態(tài)得到最大程度的保護,方可確保贛州市能夠?qū)崿F(xiàn)可持續(xù)發(fā)展。
通過構(gòu)建贛州市水資源承載力系統(tǒng)動力學(xué)模型,對贛州市的水資源承載力發(fā)展?fàn)顩r進行了情景模擬,在研究過程中,通過對社會經(jīng)濟、水資源和水環(huán)境3個子系統(tǒng)中的主要參數(shù)進行調(diào)整,對比分析了2020—2035年正常發(fā)展型、經(jīng)濟優(yōu)先型、環(huán)境優(yōu)先型和綜合協(xié)調(diào)型4種方案下的發(fā)展?fàn)顩r。研究表明:①在4種模擬方案中,從2020—2035年,社會經(jīng)濟發(fā)展趨勢來看,經(jīng)濟優(yōu)先型>綜合協(xié)調(diào)型>正常發(fā)展型>環(huán)境優(yōu)先型;②從水資源、水環(huán)境發(fā)展情況趨勢來看,環(huán)境優(yōu)先型>綜合協(xié)調(diào)型>正常發(fā)展型>經(jīng)濟優(yōu)先型。不顧水資源的承受能力單純追求經(jīng)濟與人口的高速增長,或是為了減少水資源消耗而犧牲經(jīng)濟與人口的增長,都無法保證水資源與經(jīng)濟社會的可持續(xù)發(fā)展。正常發(fā)展型和經(jīng)濟優(yōu)先型發(fā)展模式下,供需比例在2035年將分別下降至0.89、0.78,而且在未來還將持續(xù)下降,使贛州市水資源承載力壓力越來越大;環(huán)境優(yōu)先型對水體污染有緩解作用,但將在一定程度上影響社會經(jīng)濟的發(fā)展;而綜合協(xié)調(diào)型則同時考慮了經(jīng)濟和環(huán)境等各方面的要素,模擬結(jié)果較好,符合贛州市對實現(xiàn)社會經(jīng)濟與水資源的協(xié)調(diào)發(fā)展的要求,能夠更好地保障贛州市的水資源承載力得以持續(xù)健康發(fā)展。
因此,根據(jù)情景模擬的結(jié)果,在今后的水資源開發(fā)利用中,有關(guān)部門應(yīng)多管齊下,從多個方面綜合管理水資源,不僅要提高生產(chǎn)生活中水的利用效率,還要提高各種污水的處理能力,促進污水的再利用,從而降低需水總量和生產(chǎn)用水量的漲幅,提高供需比例。同時,要對企業(yè)進行嚴(yán)格的管控,促使企業(yè)根據(jù)法律法規(guī)集中合理排放廢水,對惡意排放廢水的企業(yè)必須采取措施立即停止生產(chǎn),令行禁止;加大城鄉(xiāng)水污染集中處理力度,加快城市及農(nóng)村生活污水集中處理設(shè)施建設(shè)與改造,實現(xiàn)城市、重點鄉(xiāng)鎮(zhèn)以及沿湖農(nóng)村污水處理設(shè)施建設(shè)全覆蓋,從而在保證經(jīng)濟發(fā)展的同時降低COD排放總量的增速,實現(xiàn)贛州市社會經(jīng)濟可持續(xù)發(fā)展。