陳 娟,王 森*,楊建濤,何興昌

(1.珠江水利委員會(huì)珠江水利科學(xué)研究院,廣東 廣州 510611;2 水利部珠江河口治理與保護(hù)重點(diǎn)實(shí)驗(yàn)室,廣東 廣州 510611;3 河海大學(xué)水文水資源學(xué)院,江蘇 南京 210098)

水是生命之源,水資源是城市發(fā)展的剛性約束條件。在氣候變化及人類活動(dòng)的影響下,社會(huì)經(jīng)濟(jì)發(fā)展水平與水資源條件不相匹配,水資源供需問題突出,成為制約城市發(fā)展的瓶頸[1]。因此,處理好社會(huì)經(jīng)濟(jì)與水資源間協(xié)調(diào)關(guān)系,對(duì)推動(dòng)水資源與社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展有著重要意義。大灣區(qū)地處中國沿海開放前沿,是中國開發(fā)程度最高、經(jīng)濟(jì)活力最強(qiáng)的區(qū)域之一。2019年《粵港澳大灣區(qū)發(fā)展規(guī)劃綱要》印發(fā),大灣區(qū)發(fā)展戰(zhàn)略的深入實(shí)施,社會(huì)經(jīng)濟(jì)發(fā)展對(duì)大灣區(qū)水資源的有效利用提出了更高要求。

近年來,針對(duì)水資源及相關(guān)系統(tǒng)間耦合協(xié)調(diào)關(guān)系研究已成為研究熱點(diǎn)。安莉娜等[2]以遵義為例建立耦合協(xié)調(diào)數(shù)學(xué)模型分析水資源承載力水平與城鎮(zhèn)化質(zhì)量水平的時(shí)空協(xié)調(diào)關(guān)系。張?zhí)褡说萚3]利用水資源-經(jīng)濟(jì)社會(huì)-生態(tài)環(huán)境耦合模型,分析了甘肅省水資源-經(jīng)濟(jì)社會(huì)-生態(tài)環(huán)境耦合協(xié)調(diào)的驅(qū)動(dòng)因素。常烴等[4]構(gòu)建了京津冀水資源-經(jīng)濟(jì)-社會(huì)系統(tǒng)耦合協(xié)調(diào)度模型,分析了各子系統(tǒng)和綜合系統(tǒng)的發(fā)展水平以及系統(tǒng)的協(xié)調(diào)程度。李波等[5]建立水-能源-環(huán)境系統(tǒng)的耦合協(xié)調(diào)評(píng)價(jià)指標(biāo)體系,定量分析黃河流域水-能源-環(huán)境系統(tǒng)與子系統(tǒng)之間的耦合協(xié)調(diào)發(fā)展水平。

綜上研究發(fā)現(xiàn),國內(nèi)學(xué)者研究水資源-經(jīng)濟(jì)社會(huì)-生態(tài)環(huán)境、水資源-經(jīng)濟(jì)-社會(huì)、水-能源-環(huán)境的關(guān)系較多,很少單獨(dú)研究水資源與社會(huì)經(jīng)濟(jì)的協(xié)調(diào)關(guān)系。因此,開展大灣區(qū)水資源-社會(huì)經(jīng)濟(jì)耦合協(xié)調(diào)關(guān)系評(píng)價(jià),揭示系統(tǒng)間時(shí)空演變特征,分析其主要障礙因子,對(duì)提高水資源利用率,促進(jìn)大灣區(qū)可持續(xù)發(fā)展具有一定的理論價(jià)值。

1 研究區(qū)概況與數(shù)據(jù)來源

由于資料限制,本文研究區(qū)域僅包含粵港澳大灣區(qū)中的廣東省9市,并簡稱大灣區(qū)。大灣區(qū)位于廣東省東南部,地處珠江流域尾閭,由廣州、深圳、珠海、佛山、惠州、東莞、中山、江門、肇慶9個(gè)珠三角城市組成,總面積約5.6萬km2,在國家發(fā)展大局中具有重要戰(zhàn)略地位。2021年大灣區(qū)內(nèi)常住人口7 860.60萬人,國內(nèi)生產(chǎn)總值為10.06億元,占全國生產(chǎn)總值的8.75%。

大灣區(qū)多年平均水資源總量為568億m3,入境水資源量約為2 880億m3(高要+石角+河源)。2021年大灣區(qū)水資源總量為391.3億m3,用水總量為221.0億m3,其中工業(yè)用水占29.85%、農(nóng)業(yè)用水占30.58%、生活用水占37.09%、人工生態(tài)環(huán)境補(bǔ)水占2.48%。水資源年內(nèi)、年際分布極不均勻,年徑流主要集中在汛期 4—9 月,占全年的71.3%～85.0%,而枯水季 10 月—翌年 3 月僅占全年的15.0%～28.7%。隨著大灣區(qū)社會(huì)和經(jīng)濟(jì)的發(fā)展,對(duì)水資源的需求越來越大,水資源的有限性及由此帶來的污染問題和用水量不斷增長之間的矛盾越來越突出。

數(shù)據(jù)來源于2011—2021年的《廣東省水資源公報(bào)》《廣東統(tǒng)計(jì)年鑒》。

2 研究方法與數(shù)據(jù)來源

2.1 指標(biāo)體系構(gòu)建

本文從水資源、社會(huì)經(jīng)濟(jì)2個(gè)系統(tǒng)出發(fā),水資源從水資源條件及水資源利用2方面考慮,社會(huì)經(jīng)濟(jì)從經(jīng)濟(jì)水平、社會(huì)結(jié)構(gòu)和城市建設(shè)3方面考慮。水資源條件從水資源特征出發(fā),選取了產(chǎn)水系數(shù)、產(chǎn)水模數(shù)和人均水資源量;水資源利用從生活、工業(yè)、農(nóng)業(yè)、及地方生產(chǎn)總值出發(fā),選取了人均綜合用水量、萬元工業(yè)增加值用水量、耕地實(shí)際灌溉畝均用水量和萬元地區(qū)生產(chǎn)總值用水量。經(jīng)濟(jì)水平主要考慮了人均GDP和人均可支配收入;社會(huì)結(jié)構(gòu)從人口、城市化及產(chǎn)業(yè)結(jié)構(gòu)出發(fā)選取了人口密度、城鎮(zhèn)化和第三產(chǎn)業(yè)生產(chǎn)總值占比;城市建設(shè)主要從生態(tài)環(huán)境考慮,選取了城市污水處理率和城市人均公園綠地面積。社會(huì)經(jīng)濟(jì)涵蓋了社會(huì)、經(jīng)濟(jì)、環(huán)境的相關(guān)指標(biāo)。構(gòu)建系統(tǒng)-準(zhǔn)則-指標(biāo)3個(gè)層次共計(jì)14個(gè)評(píng)價(jià)指標(biāo)的大灣區(qū)水資源-社會(huì)經(jīng)濟(jì)評(píng)價(jià)指標(biāo)體系,見表1。

表1 大灣區(qū)水資源-社會(huì)經(jīng)濟(jì)評(píng)價(jià)指標(biāo)體系

表2 耦合協(xié)調(diào)狀態(tài)判別標(biāo)準(zhǔn)

2.2 研究方法

2.2.1指標(biāo)標(biāo)準(zhǔn)化

對(duì)指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理,保證各指標(biāo)具有可比性,見式(1)、(2)。

正向指標(biāo):yij=(xij-minxj)/(maxxj-minxj)

(1)

負(fù)向指標(biāo):yij=(maxxj-xij)/(maxxj-minxj)

(2)

式中yij——指標(biāo)標(biāo)準(zhǔn)化后的值;maxxj——指標(biāo)的最大值;minxj——指標(biāo)的最小值;i——第i個(gè)樣本;j——第j個(gè)指標(biāo)。

2.2.2指標(biāo)權(quán)重

采用組合權(quán)重w,其中主觀權(quán)重w1采用層次分析法,具體計(jì)算見文獻(xiàn)[6],客觀權(quán)重w2采用熵權(quán)法,具體計(jì)算見文獻(xiàn)[7]。利用最小信息熵原理構(gòu)造模型[8],見式(3)、(4)。

(3)

(4)

用拉格朗日乘數(shù)法求解,見式(5)。

(5)

式中w(i)——第i個(gè)指標(biāo)組合權(quán)重;w1(i)——第i個(gè)指標(biāo)主觀權(quán)重;w2(i)——第i個(gè)指標(biāo)客觀權(quán)重;n——指標(biāo)數(shù)。

2.2.3綜合值

水資源系統(tǒng)、社會(huì)經(jīng)濟(jì)系統(tǒng)的綜合值利用線性加權(quán)法進(jìn)行求解,見式(6)。

(6)

式中U1——水資源系統(tǒng)綜合值;U2——社會(huì)經(jīng)濟(jì)系統(tǒng)綜合值;ωj——第j個(gè)指標(biāo)組合權(quán)重;m——水資源系統(tǒng)指標(biāo)數(shù);n——社會(huì)經(jīng)濟(jì)系統(tǒng)指標(biāo)數(shù)。

2.2.4耦合協(xié)調(diào)度模型

為了判定系統(tǒng)間協(xié)調(diào)發(fā)展關(guān)系,采用耦合協(xié)調(diào)度模型,通過計(jì)算系統(tǒng)間的相對(duì)離差系數(shù)量化系統(tǒng)間協(xié)調(diào)程度[9],見式(7)、(8)。

(7)

(8)

式中α——水資源系統(tǒng)權(quán)重,β——社會(huì)經(jīng)濟(jì)系統(tǒng)權(quán)重,考慮兩系統(tǒng)同等重要,α=β=0.5;C——耦合度;D——耦合協(xié)調(diào)度,值在0～1。

根據(jù)D值大小劃分為 5 級(jí)[10-11],結(jié)合水資源與社會(huì)經(jīng)濟(jì)相對(duì)關(guān)系,對(duì)于失調(diào)和過渡階段進(jìn)一步劃分為社會(huì)經(jīng)濟(jì)過快發(fā)展和社會(huì)經(jīng)濟(jì)過慢發(fā)展,具體見表 2。

2.2.5障礙因子診斷模型

為了量化指標(biāo)對(duì)指標(biāo)體系的影響程度,引入障礙度函數(shù),判別主要障礙因子。根據(jù)評(píng)價(jià)指標(biāo)在體系中權(quán)重和指標(biāo)偏離度來推求障礙度,見式(9)、(10)[12-13]。

Iij=1-yij

(9)

(10)

式中Iij——第j個(gè)樣本第i個(gè)指標(biāo)偏離度;Hij——第j個(gè)樣本第i個(gè)指標(biāo)障礙度;n——指標(biāo)數(shù)。

指標(biāo)障礙度越小,說明體系受該項(xiàng)指標(biāo)的影響越弱;指標(biāo)障礙度越大,說明體系受該項(xiàng)指標(biāo)的影響越強(qiáng)。

3 結(jié)果與分析

3.1 水資源-社會(huì)經(jīng)濟(jì)耦合協(xié)調(diào)時(shí)空演變

3.1.1時(shí)間演變

為探索大灣區(qū)水資源-社會(huì)經(jīng)濟(jì)耦合協(xié)調(diào)發(fā)展的時(shí)間演變,采用2011—2021年指標(biāo)數(shù)據(jù)計(jì)算耦合協(xié)調(diào)度,繪制趨勢(shì),見圖1。

2011—2021年水資源系統(tǒng)綜合值變化不大,在0.18～0.26,其中2011、2021年值最低,為0.18。2011、2021年降水量均低于常年,屬偏枯水年,水資源條件相對(duì)較差。2011—2021年社會(huì)經(jīng)濟(jì)系統(tǒng)綜合值呈緩慢上升趨勢(shì),2019年達(dá)到0.30;受疫情影響,2020、2021年社會(huì)經(jīng)濟(jì)略低于2019年,均為0.29。

2011—2021年大灣區(qū)水資源-社會(huì)經(jīng)濟(jì)協(xié)調(diào)耦合度在0.43～0.51,處于過渡階段低度協(xié)調(diào)狀態(tài);相對(duì)于水資源系統(tǒng),社會(huì)經(jīng)濟(jì)系統(tǒng)發(fā)展過快,且與水資源系統(tǒng)間差距有增大趨勢(shì)。

3.1.2空間特征

為揭示大灣區(qū)水資源-社會(huì)經(jīng)濟(jì)耦合協(xié)調(diào)發(fā)展的空間特征,選取2011、2016、2021年3個(gè)典型年份繪制空間分布,見圖2。

a)2011年

2011年,大灣區(qū)各市耦合協(xié)調(diào)度在0.35～0.47。處于輕度失調(diào)狀態(tài)的有佛山,其他均處于低度協(xié)調(diào)狀態(tài),其中耦合協(xié)調(diào)度最高的為珠海。廣州、佛山和肇慶3個(gè)市屬于社會(huì)經(jīng)濟(jì)過快發(fā)展,深圳、珠海、惠州、東莞、中山和江門6個(gè)市屬于社會(huì)經(jīng)濟(jì)過慢發(fā)展。

2016年,大灣區(qū)各市耦合協(xié)調(diào)度在0.44～0.56,均處于低度協(xié)調(diào)狀態(tài),其中耦合協(xié)調(diào)度最高的為珠海,最低的為佛山。廣州、佛山、東莞和中山4個(gè)市屬于社會(huì)經(jīng)濟(jì)過快發(fā)展,深圳、珠海、惠州、江門和肇慶5個(gè)市屬于社會(huì)經(jīng)濟(jì)過慢發(fā)展。

2021年,大灣區(qū)各市耦合協(xié)調(diào)度在0.42～0.53,均處于低度協(xié)調(diào)狀態(tài),其中耦合協(xié)調(diào)度最高的為珠海,最低的為惠州。9個(gè)市均屬于社會(huì)經(jīng)濟(jì)過快發(fā)展。隨著大灣區(qū)社會(huì)經(jīng)濟(jì)發(fā)展,水資源滯后于社會(huì)經(jīng)濟(jì)的城市也有增多的趨勢(shì)。

從區(qū)域看,珠海水資源條件較好,珠江流域八大口門中的崖門、磨刀門、雞啼門和虎跳門都流經(jīng)珠海,其中縱貫珠海的磨刀門水道多年平均徑流量約923億m3,為八大口門之首,過境水資源量豐富;珠海屬亞熱帶海洋性氣候,海洋對(duì)本地氣候調(diào)節(jié)作用明顯,常年降雨充沛。珠海社會(huì)經(jīng)濟(jì)也比較突出,人均GDP、人均可支配收入均高于大灣區(qū)平均值,而人口密度僅為大灣區(qū)平均值的1/2。佛山地處珠江三角洲,本市用水主要依賴于過境水量,人均水資源量接近廣東省平均值的1/3,大灣區(qū)平均值的1/2。佛山人均GDP、人均可支配收入接近大灣區(qū)平均值,城鎮(zhèn)化率明顯高于大灣區(qū)平均值,第三產(chǎn)業(yè)生產(chǎn)總值占比明顯低于大灣區(qū)平均值。

3.2 障礙因子分析

利用障礙因子診斷模型計(jì)算大灣區(qū)2011—2021年14個(gè)指標(biāo)的障礙度(%),見圖3。

圖3 水資源-社會(huì)經(jīng)濟(jì)耦合協(xié)調(diào)障礙因子及障礙度

人均水資源量是大灣區(qū)水資源-社會(huì)經(jīng)濟(jì)耦合協(xié)調(diào)的主要障礙因子,障礙度為14.5%。第二為人均GDP,障礙度為13.4%。其后3個(gè)指標(biāo)依次為第三產(chǎn)業(yè)生產(chǎn)總值占比、城鎮(zhèn)化率和耕地實(shí)際灌溉畝均用水量,其障礙度分別為10.9%、9.8%和9.5%。在水資源系統(tǒng)中,人均水資源量受降水影響較大,不易管控,可通過調(diào)水工程、非常規(guī)水資源利用來部分緩解。正在實(shí)施的珠江三角洲水資源配置工程從西江干流水道向珠江三角洲東部地區(qū)引水,解決廣州、東莞、深圳城市缺水問題。加強(qiáng)非常規(guī)水資源利用,推廣再生水、雨水、淡化海水等利用方式。未來以增加人均GDP、增大第三產(chǎn)業(yè)生產(chǎn)總值占比,適當(dāng)控制城市發(fā)展速度以及優(yōu)化農(nóng)業(yè)灌溉模式減小灌溉畝均用水量成為推動(dòng)大灣區(qū)水資源-社會(huì)經(jīng)濟(jì)耦合協(xié)調(diào)發(fā)展的重點(diǎn)。加快經(jīng)濟(jì)結(jié)構(gòu)轉(zhuǎn)型升級(jí),大力發(fā)展新興產(chǎn)業(yè),提高經(jīng)濟(jì)活力,促進(jìn)用水方式轉(zhuǎn)變。

4 討論

大灣區(qū)各市嚴(yán)格落實(shí)水資源總量控制制度,積極推進(jìn)節(jié)水型社會(huì)、水生態(tài)文明城市和海綿城市建設(shè),強(qiáng)化節(jié)水減排、再生水利用等。大灣區(qū)用水效率明顯提高,萬元地區(qū)生產(chǎn)總值用水量和萬元工業(yè)增加值用水量均顯著下降,人均綜合用水量逐漸下降,耕地實(shí)際灌溉畝均用水量略有下降。隨著大灣區(qū)經(jīng)濟(jì)快速發(fā)展,社會(huì)經(jīng)濟(jì)系統(tǒng)發(fā)展過快,與水資源系統(tǒng)間差距有增大趨勢(shì),說明水資源形勢(shì)仍比較嚴(yán)峻,影響水資源社會(huì)經(jīng)濟(jì)協(xié)調(diào)發(fā)展的主要障礙因子有人均水資源量、第三產(chǎn)業(yè)生產(chǎn)總值占比、城鎮(zhèn)化率和耕地實(shí)際灌溉畝均用水量。因此,在落實(shí)水資源節(jié)約與保護(hù)現(xiàn)有制度和措施的前提下,進(jìn)一步優(yōu)化水資源配置,大力推進(jìn)雨洪資源利用等節(jié)約水、涵養(yǎng)水的工程建設(shè),構(gòu)建綠色生態(tài)水網(wǎng),推進(jìn)水資源全面節(jié)約和循環(huán)利用。強(qiáng)化城市建設(shè)管理,集中力量著力調(diào)整用水結(jié)構(gòu),提高用水效率,促進(jìn)經(jīng)濟(jì)發(fā)展方式加快轉(zhuǎn)變,促進(jìn)社會(huì)經(jīng)濟(jì)與水資源協(xié)調(diào)發(fā)展。

本研究中不確定性主要源于評(píng)價(jià)指標(biāo)體系的構(gòu)建,在指標(biāo)體系上主要考慮了水資源與社會(huì)經(jīng)濟(jì)耦合協(xié)調(diào)關(guān)系,其中社會(huì)經(jīng)濟(jì)涵蓋了社會(huì)、經(jīng)濟(jì)、環(huán)境的相關(guān)指標(biāo)。由于數(shù)據(jù)的可獲取性及體系構(gòu)建的科學(xué)性等因素,指標(biāo)體系并不能全面反映系統(tǒng)間的耦合協(xié)調(diào)關(guān)系,今后的研究仍需繼續(xù)完善。

因資料數(shù)據(jù)限制,尚未將香港、澳門納入分析,而香港、澳門水資源與大灣區(qū)水資源緊密聯(lián)系,大灣區(qū)對(duì)香港供水滿足了香港80%左右的用水需求,對(duì)澳門供水滿足了澳門98%以上的用水需求。后續(xù)應(yīng)加強(qiáng)粵港澳大灣區(qū)整體水資源研究。

5 結(jié)論

通過構(gòu)建指標(biāo)體系,運(yùn)用耦合協(xié)調(diào)度模型和障礙因子診斷模型對(duì)大灣區(qū)2011—2021年水資源-社會(huì)經(jīng)濟(jì)耦合協(xié)調(diào)時(shí)空特征和主要障礙因子進(jìn)行分析,主要結(jié)論如下。

a)2011—2021年,大灣區(qū)水資源與社會(huì)經(jīng)濟(jì)協(xié)調(diào)水平較低,社會(huì)經(jīng)濟(jì)發(fā)展過快,且差距增大;九市基本處于低度協(xié)調(diào)狀態(tài),耦合協(xié)調(diào)度珠海較高,佛山相對(duì)較低。

b)水資源與社會(huì)經(jīng)濟(jì)均存在影響協(xié)調(diào)發(fā)展的障礙因子。水資源主要障礙因子有人均水資源量、耕地實(shí)際灌溉畝均用水量。社會(huì)經(jīng)濟(jì)主要障礙因子有人均GDP、第三產(chǎn)業(yè)生產(chǎn)總值占比、城鎮(zhèn)化率。大灣區(qū)各市主要障礙因子空間差異性不大。