叢東來, 于少鵬, 陳 曦, 劉 獻

(1.哈爾濱學院 地理與旅游學院, 黑龍江 哈爾濱 150086; 2.黑龍江省寒區(qū)濕地生態(tài)與環(huán)境研究重點實驗室,黑龍江 哈爾濱 150086; 3.哈爾濱濕地研究院, 黑龍江 哈爾濱 150086; 4.哈爾濱學院 食品工程學院, 黑龍江 哈爾濱 150086)

城鎮(zhèn)化是以人口和經(jīng)濟活動向特定地域集中為基本特征的發(fā)展過程，是現(xiàn)代化的必由之路。按國家行政建制設立的市、鎮(zhèn)是城鎮(zhèn)化發(fā)展的重要空間載體，而可供開發(fā)利用的各類資源則是城鎮(zhèn)化發(fā)展不可或缺的物質(zhì)保障。當前中國城鎮(zhèn)化率已接近60%，未來25 a中國城鎮(zhèn)化水平每提高1%所需水量、建設用地和造成的生態(tài)超載分別是過去25 a的1.88倍，3.45，2.42倍[1]，表明未來取水難度、用地難度和生態(tài)重建難度會越來越大。水是自然系統(tǒng)中最為活躍的因子，具有原生的生態(tài)屬性和次生的“生態(tài)—社會”關(guān)聯(lián)屬性[2]，水作為基礎性的自然資源和戰(zhàn)略性的經(jīng)濟資源[3]，已成為中國建設“生態(tài)文明”和進行“綠色發(fā)展”的重要支撐條件。新世紀以來，伴隨著中國城鎮(zhèn)數(shù)量的增加和城鎮(zhèn)面積的不斷擴大，城鎮(zhèn)化發(fā)展與水資源、水環(huán)境和水生態(tài)的相互作用日趨復雜。因此，探討水資源開發(fā)利用與城鎮(zhèn)化發(fā)展的復雜關(guān)系具有現(xiàn)實價值和長遠的戰(zhàn)略意義。

由于發(fā)達國家城鎮(zhèn)化發(fā)展較早，并且目前多數(shù)發(fā)達國家已進入城鎮(zhèn)化進程后期，因此國外對城鎮(zhèn)化進程與水資源開發(fā)利用相互關(guān)系的研究成果較多。主要涉及水資源管理、政策與城鎮(zhèn)化環(huán)境的關(guān)系[4]、水污染對城鎮(zhèn)化的影響[5]、城鎮(zhèn)化增長與水資源危機[6]、城鎮(zhèn)化發(fā)展階段與水資源開發(fā)利用階段的對應關(guān)系[7]等。國內(nèi)學者對城鎮(zhèn)化與水資源開發(fā)利用關(guān)系的研究，從內(nèi)容上看主要包括城鎮(zhèn)化與水資源系統(tǒng)的動態(tài)耦合及其效應分析[8]，城鎮(zhèn)化與水足跡、水污染的關(guān)系[9]，城鎮(zhèn)化進程中水資源承載力的測度[10]，城鎮(zhèn)化與水系統(tǒng)的演變[11]，城鎮(zhèn)化與水資源利用的互動機理和調(diào)控模式[12]，水資源對城鎮(zhèn)化發(fā)展的約束作用[13]等；從研究區(qū)域看，多集中在西北干旱、半干旱地區(qū)、華北地區(qū)和兩湖地區(qū)[14-16]；從研究的空間尺度看，以全國范圍、流域范圍和規(guī)模較大的城市群地區(qū)為主[17-19]。整體而言，現(xiàn)有成果以宏觀、中觀尺度的城鎮(zhèn)化發(fā)展為背景，探討城鎮(zhèn)化進程中水資源的承載力、約束力問題較多，而針對人水關(guān)系、地水關(guān)系相對緊張的城鎮(zhèn)化核心城市的研究尚有不足?；诖?，本文選取北方地區(qū)典型的缺水型中心城市——哈爾濱市為研究對象，評價其水資源開發(fā)利用與城鎮(zhèn)化發(fā)展狀態(tài)，并建立二者的響應關(guān)系模型，以期揭示城鎮(zhèn)化進程與水資源開發(fā)利用的作用機制。

1 研究區(qū)域與研究方法

1.1 研究區(qū)概況

哈爾濱市是黑龍江省省會，地處東北亞中心位置，是中國東北北部的區(qū)域性中心城市、哈長城市群的核心城市之一和第一批國家新型城鎮(zhèn)化綜合試點城市。同時，哈爾濱也是中國重要的制造業(yè)基地和商品糧基地。作為東北地區(qū)“T”型經(jīng)濟網(wǎng)絡的關(guān)鍵節(jié)點和哈大齊工業(yè)走廊的起點，哈爾濱是地區(qū)經(jīng)濟的增長極和人口的聚集地。截止到2017年末，哈爾濱市轄9區(qū)9縣(市)，全市常住人口1092.9萬人，其中戶籍總?cè)丝?55萬人，常住人口城鎮(zhèn)化水平64.5%；市域總面積53 076 km2，其中市轄區(qū)面積10 193 km2，建成區(qū)面積434 km2；國內(nèi)生產(chǎn)總值6.355億元，占全省國內(nèi)生產(chǎn)總值39.2%。

哈爾濱市為北溫帶大陸性季風氣候，年降水變率較大。水資源空間分布呈東—中—西依次遞減的格局，具有東豐西欠、貧富不均的特征。市域多年平均降水量為611 mm，6月至9月汛期徑流量占多年平均徑流量的60%～75%，多年平均水資源量為1.14×1010m3，地下水資源總量約為3.00×109m3。松花江干流中段在哈爾濱穿城而過，城市過境水資源相對豐富。2017年，全市人均綜合用水量708.3 m3，萬元工業(yè)增加值用水量22.9 m3。按照國際通行的水資源緊缺指標判斷[20]，哈爾濱市屬于中、重度缺水地區(qū)。2008年到2013年，哈爾濱市先后被列為第三批全面節(jié)水型社會建設試點地區(qū)，國家級水生態(tài)系統(tǒng)保護與修復試點城市，第一批全國水生態(tài)文明城市建設試點。但受水資源時空分布不均，水質(zhì)污染和用水供需矛盾突出等因素影響[21]，哈爾濱市社會經(jīng)濟發(fā)展與水資源開發(fā)利用的形勢依然嚴峻。

1.2 指標體系和研究方法

系統(tǒng)論的視角有利于揭示“人—地—水”復合系統(tǒng)內(nèi)部的相互作用及其整體性特征，按照“要素—子系統(tǒng)—系統(tǒng)”的構(gòu)成關(guān)系[22]，本研究通過選取要素指標構(gòu)建城鎮(zhèn)化系統(tǒng)和水資源系統(tǒng)，評價其發(fā)展狀況，并運用響應度模型分析二者的響應關(guān)系。

(1) 要素指標。通過理論分析和頻度統(tǒng)計法初步選取城鎮(zhèn)化系統(tǒng)和水資源系統(tǒng)的基本要素指標。運用SPSS19.0軟件篩選出相關(guān)系數(shù)較大的指標并予以去除，再利用因子分析法，提取累計貢獻率超過85%的指標作為終選指標[16]。由于水資源系統(tǒng)要素指標的統(tǒng)計數(shù)量及統(tǒng)計年份都有一定的局限，因此綜合考慮數(shù)據(jù)的可得性，對少部分相關(guān)系數(shù)較高但對系統(tǒng)評價重要的要素指標予以保留，最終確定的要素指標詳見表1。

表1 哈爾濱市城鎮(zhèn)化系統(tǒng)與水資源系統(tǒng)評價體系及指標權(quán)重

(2) 系統(tǒng)、子系統(tǒng)綜合評價指數(shù)。采用層次分析法，結(jié)合熵權(quán)法對初始指標進行修正[23]，并進行標準化處理。運用逐層加權(quán)法分別計算系統(tǒng)、子系統(tǒng)綜合評價指數(shù)。

(3) 水資源系統(tǒng)與城鎮(zhèn)化系統(tǒng)的響應指數(shù)和響應度模型。響應指數(shù)通常用響應度來表達，其內(nèi)涵來自于經(jīng)濟學點彈性的概念，表示因變量對自變量變化的敏感程度。響應度模型[24]可以揭示水資源開發(fā)利用與城鎮(zhèn)化發(fā)展的響應關(guān)系及作用強度。

(1)

(4) 水資源子系統(tǒng)和城鎮(zhèn)化子系統(tǒng)的相互作用。采用SPSS19.0軟件，對水資源和城鎮(zhèn)化的子系統(tǒng) 相關(guān)系數(shù)矩陣進行運算，得出相關(guān)系數(shù)及其p值，據(jù)此判斷二者之間的關(guān)聯(lián)程度和作用機制。

1.3 數(shù)據(jù)來源

本研究采用的2006—2017年城鎮(zhèn)化發(fā)展和水資源開發(fā)利用的數(shù)據(jù)，來自于相應的《哈爾濱統(tǒng)計年鑒》《哈爾濱市國民經(jīng)濟和社會發(fā)展統(tǒng)計公報》《中國城市統(tǒng)計年鑒》《哈爾濱市水資源公報》，以及哈爾濱市水務局年度工作總結(jié)和部門決算(網(wǎng)址：http:∥www.hrbwrb.gov.cn/n_zwgk/n_bmysxx/)。對個別年份數(shù)據(jù)缺失的百戶城市居民汽車擁有量、生態(tài)環(huán)境用水量、城鎮(zhèn)生活用水量、灌溉水有效利用系數(shù)指標，進行了相關(guān)部門補充調(diào)查或相鄰年份插值法處理。

2 結(jié)果與分析

2.1 城鎮(zhèn)化系統(tǒng)綜合評價

2.1.1 城鎮(zhèn)化系統(tǒng)發(fā)展的主要影響因素 哈爾濱市城鎮(zhèn)化系統(tǒng)的主要影響因素可通過其子系統(tǒng)和要素指標的權(quán)重來體現(xiàn)(見表1)。4個子系統(tǒng)的影響力從大到小依次為經(jīng)濟城鎮(zhèn)化(0.307)、人口城鎮(zhèn)化(0.278)、社會城鎮(zhèn)化(0.242)和空間城鎮(zhèn)化(0.173)。前三者權(quán)重合計82.7%，對哈爾濱城鎮(zhèn)化系統(tǒng)的影響較大，經(jīng)濟城鎮(zhèn)化子系統(tǒng)起主導作用。在評價哈爾濱城鎮(zhèn)化系統(tǒng)的18個要素指標中，二三產(chǎn)業(yè)經(jīng)濟密度(0.178)、城鎮(zhèn)人口比重(0.136)、人均GDP(0.091)等9個指標的權(quán)重合計約80%，是哈爾濱城鎮(zhèn)化系統(tǒng)發(fā)展的主要影響因素，主要影響因素指標的權(quán)重分布與城鎮(zhèn)化子系統(tǒng)權(quán)重值相互對應，可以看出影響哈爾濱城鎮(zhèn)化系統(tǒng)發(fā)展的主要影響因素集中在二三產(chǎn)業(yè)經(jīng)濟結(jié)構(gòu)、城鎮(zhèn)人口的規(guī)模分布和生產(chǎn)消費方面。

2.1.2 城鎮(zhèn)化系統(tǒng)的時序演化特征 由圖1可知，2006—2017年哈爾濱城鎮(zhèn)化系統(tǒng)及經(jīng)濟、社會、空間子系統(tǒng)均呈同步上升、增速放緩的態(tài)勢。人口城鎮(zhèn)化的發(fā)展變化較為劇烈和復雜，呈前期陡增、中期波動、后期趨穩(wěn)態(tài)勢。2006—2014年，經(jīng)濟城鎮(zhèn)化和社會城鎮(zhèn)化滯后于人口城鎮(zhèn)化和空間城鎮(zhèn)化，說明這期間哈爾濱城鎮(zhèn)化發(fā)展主要受人口城鎮(zhèn)化驅(qū)動，其次是空間城鎮(zhèn)化驅(qū)動。具體表現(xiàn)為2006—2014年城鎮(zhèn)人口數(shù)量增長了1.19倍，建成區(qū)面積增加了1.33倍，城鄉(xiāng)人口轉(zhuǎn)變加快，外延型城市土地擴張明顯，是典型的粗放型城鎮(zhèn)化發(fā)展；同期經(jīng)濟城鎮(zhèn)化和社會城鎮(zhèn)化雖然持續(xù)增長，但整體效率不高，城鎮(zhèn)化質(zhì)量有待提升。2014以后空間城鎮(zhèn)化、經(jīng)濟城鎮(zhèn)化和社會城鎮(zhèn)化穩(wěn)步上升，人口城鎮(zhèn)化的發(fā)展趨于平穩(wěn)。原因在于哈爾濱在城市發(fā)展戰(zhàn)略上注重人口—經(jīng)濟—空間—社會發(fā)展的協(xié)調(diào)性，在空間布局上延續(xù)“北躍—南拓—中興—強縣”的空間均衡思想，人口集聚更加趨于理性；而空間城鎮(zhèn)化仍有較快發(fā)展原因在于哈爾濱市雙城等地撤市劃區(qū)的行政區(qū)劃變動、二三產(chǎn)業(yè)的布局調(diào)整，促進了哈爾濱都市圈內(nèi)部土地利用的優(yōu)化。

圖1 哈爾濱市城鎮(zhèn)化系統(tǒng)及其子系統(tǒng)綜合評價指數(shù)

根據(jù)哈爾濱城鎮(zhèn)化系統(tǒng)綜合評價指數(shù)及其年增長率的變化，可將2006—2017年哈爾濱城鎮(zhèn)化系統(tǒng)的發(fā)展分為兩個階段： ①2006—2010年城鎮(zhèn)化快速增長階段。此階段是東北振興戰(zhàn)略啟動以后的紅利釋放期，哈爾濱城鎮(zhèn)化發(fā)展迎來較大機遇。②2011—2017年城鎮(zhèn)化階梯式增長階段。這一階段中國先后提出了新型城鎮(zhèn)化內(nèi)涵式發(fā)展和經(jīng)濟新常態(tài)，東北經(jīng)濟處在轉(zhuǎn)型升級的進程中，哈爾濱市通過出臺一系列以“緊湊型城市”和“精明增長”為目標的政策法規(guī)，有效控制了城市建設用地的無序蔓延，逐步提升城鎮(zhèn)化效率。

2.2 水資源系統(tǒng)綜合評價

2.2.1 水資源系統(tǒng)開發(fā)利用的影響因素 哈爾濱市水資源系統(tǒng)的4個子系統(tǒng)影響作用從大到小依次為水資源利用(0.383)、水資源管理(0.270)、水資源稟賦(0.179)、水資源荷載(0.168)。在構(gòu)成水資源系統(tǒng)的18個要素指標中，水資源開發(fā)利用率(0.143)、萬元工業(yè)增加值廢水排放量(0.105)、萬元GDP用水量(0.099)等9個指標的權(quán)重和為80%，其中4個與工業(yè)用水密切相關(guān)的指標約占總權(quán)重的1/3以上，即哈爾濱水資源系統(tǒng)開發(fā)利用的主要影響因素是工業(yè)經(jīng)濟的用水需求和用水效率，其次是水利行政管理和農(nóng)業(yè)用水因素。

2.2.2 水資源系統(tǒng)的時序演化特征 由圖2可知，2006—2017年哈爾濱市水資源系統(tǒng)整體呈上升發(fā)展態(tài)勢，研究期內(nèi)系統(tǒng)波動變化明顯，波動幅度逐漸減小。各子系統(tǒng)中： ①水資源稟賦波動最為劇烈，主要原因是受季風區(qū)自然條件影響，哈爾濱市年降水變率較大，研究期內(nèi)降雨量最大值為2013年的4.10×1010m3，最小值為2007年的2.15×1010m3，分別比多年平均降雨量浮動+26.4%和-33.6%； ②水資源荷載在研究期的首末時間點上變動不大，水資源荷載略有增加，但整體發(fā)展過程中的變動十分明顯。2006—2008年荷載的指標值明顯高于水資源其他3個子系統(tǒng)，表明哈爾濱市在快速的城鎮(zhèn)化進程中，人口和社會經(jīng)濟的發(fā)展壓力造成水資源負荷過重；2008—2013年水資源負荷持續(xù)降低，與推行嚴格的節(jié)水制度、建設節(jié)水型城市以及產(chǎn)業(yè)用水調(diào)整有很大關(guān)聯(lián)；2013—2017年用水負荷再次持續(xù)增加，主要是受氣候條件影響，農(nóng)業(yè)用水增加較多； ③水資源管理整體上不斷提升進步，反映在水利行政經(jīng)費投入占比逐年增加、城鎮(zhèn)污水處理率、水質(zhì)監(jiān)測達標率持續(xù)上升等方面，管理系統(tǒng)的穩(wěn)步發(fā)展得益于哈爾濱水資源管理制度的高效落實和水處理技術(shù)的應用推廣； ④水資源利用與水資源系統(tǒng)的變動趨勢大體保持同步。萬元GDP用水量連續(xù)下降，下降幅度達132%，工業(yè)產(chǎn)業(yè)的水依賴逐漸減??；工業(yè)用水重復率和灌溉水有效利用系數(shù)持續(xù)增大，得益于水處理技術(shù)的提升和水利工程設施的增加。但是，在技術(shù)進步和城鎮(zhèn)化內(nèi)涵發(fā)展的背景下，哈爾濱市人均用水量和水資源開發(fā)利用率(用水量與水資源總量的比值)兩項指標仍有明顯增加，具體表現(xiàn)為工業(yè)用水量減小的同時，農(nóng)業(yè)用水量不斷增加，并且哈爾濱市農(nóng)田灌溉用水占全市各產(chǎn)業(yè)用水比重常年保持在80%以上。

綜合考慮水資源系統(tǒng)及其子系統(tǒng)的變化，可以將2006—2017年哈爾濱市水資源系統(tǒng)的發(fā)展劃分為兩個階段。①2006—2011年波動發(fā)展階段，這一階段前期哈爾濱市水資源系統(tǒng)發(fā)展退化的壓力較大，后期因國家生態(tài)、環(huán)保政策調(diào)整和區(qū)域增長方式轉(zhuǎn)變使得水資源系統(tǒng)快速恢復和提升； ②2012—2017年階梯式發(fā)展階段，此階段在生態(tài)水文明背景下哈爾濱市進行了水生態(tài)保護與修復試點城市建設，以技術(shù)進步和產(chǎn)業(yè)升級為基礎水資源系統(tǒng)得以穩(wěn)步發(fā)展。水資源系統(tǒng)與城鎮(zhèn)化系統(tǒng)的階段劃分在時間節(jié)點上較為接近，表明二者的發(fā)展在系統(tǒng)內(nèi)外部環(huán)境上存在一定關(guān)聯(lián)。由于城鎮(zhèn)化系統(tǒng)和水資源系統(tǒng)在自身各階段的發(fā)展速度和內(nèi)容并不一致，必然導致兩個系統(tǒng)間產(chǎn)生差異性的響應。

圖2 哈爾濱市水資源系統(tǒng)及其子系統(tǒng)綜合評價指數(shù)

2.3 水資源系統(tǒng)與城鎮(zhèn)化系統(tǒng)的響應關(guān)系分析

2.3.1 水資源系統(tǒng)與城鎮(zhèn)化系統(tǒng)響應關(guān)系的函數(shù)表達 一般而言，兩系統(tǒng)的響應關(guān)系是由A→B和B→A共同組成。A對B的響應(B→A)，是指以B為自變量，以A為因變量，考察B變動時A的變動幅度，反之亦然。因此，表征水資源系統(tǒng)與城鎮(zhèn)化系統(tǒng)響應關(guān)系的函數(shù)如表2和表3所示。使用SPSS19.0軟件，進行線性模型、指數(shù)模型、冪函數(shù)模型、二次項模型、三次項模型、邏輯斯蒂模型、生長曲線模型等多種曲線估計分析，對城鎮(zhèn)化和水資源的系統(tǒng)及子系統(tǒng)的響應指數(shù)分別進行曲線擬合，通過比較擬合優(yōu)度建立最佳響應函數(shù)(表2和表3)。從響應函數(shù)的類型看，除水資源稟賦對城鎮(zhèn)化的響應方程為指數(shù)函數(shù)、經(jīng)濟城鎮(zhèn)化和社會城鎮(zhèn)化對水資源系統(tǒng)的響應方程為二次函數(shù)以外，其他響應函數(shù)均為三次函數(shù)。對擬合的響應函數(shù)進行檢驗，當F>F0.1(1,10)=3.285，p<0.1時，響應函數(shù)具有顯著性意義。測定系數(shù)R2值越大表明擬合優(yōu)度越高。

表2 2006-2017年哈爾濱市水資源系統(tǒng)對城鎮(zhèn)化系統(tǒng)的響應函數(shù)及其檢驗

注：表中Y,Y1,Y2,Y3,Y4表示水資源系統(tǒng)和水資源稟賦、荷載、利用、管理子系統(tǒng)；X表示城鎮(zhèn)化系統(tǒng)。

表3 2006-2017年哈爾濱市城鎮(zhèn)化系統(tǒng)對水資源系統(tǒng)的響應函數(shù)及其檢驗

注：X，X1，X2，X3，X4表示城鎮(zhèn)化系統(tǒng)和城鎮(zhèn)化人口、經(jīng)濟、社會、空間子系統(tǒng)；Y表示水資源系統(tǒng)。

2.3.2 水資源系統(tǒng)與城鎮(zhèn)化系統(tǒng)的響應關(guān)系 從水資源對城鎮(zhèn)化的響應來看(見表2)： ①水資源管理與城鎮(zhèn)化擬合程度最高，其次是水資源荷載和水資源利用，水資源稟賦與城鎮(zhèn)化系統(tǒng)擬合度較低。說明城鎮(zhèn)化發(fā)展通過促進技術(shù)進步和提高社會組織效能，進而對水資源管理能力提升產(chǎn)生顯著影響。而水資源稟賦受自然條件影響較大，受人類活動影響相對較?。虎谒Y源荷載和水資源管理對城鎮(zhèn)化系統(tǒng)的響應度較大，表明哈爾濱市的用水和水資源管理對城鎮(zhèn)化發(fā)展有更強的敏感性。一方面2006—2017年水資源荷載的發(fā)展呈劇烈波動的雙U結(jié)構(gòu)(見圖2)，與城鎮(zhèn)化的發(fā)展軌跡并不具有對應性；另一方面在某些年份，如2006，2013，2017年水資源荷載與水資源稟賦的變化呈反向?qū)P(guān)系，反映出水資源荷載的響應強度不僅與城鎮(zhèn)化自身演進有關(guān)，也受區(qū)域降水等自然因素影響； ③水資源系統(tǒng)綜合評價指數(shù)位于擬合曲線上方(見圖3)，反映出這些年份哈爾濱市水資源系統(tǒng)具有較大的承載空間，能夠?qū)Τ擎?zhèn)化發(fā)展起支撐作用，保障城鎮(zhèn)化健康發(fā)展；其余年份水資源系統(tǒng)綜合評價指數(shù)位于擬合曲線下方，表明城鎮(zhèn)化對水資源系統(tǒng)的消耗較大，水資源對城鎮(zhèn)化發(fā)展產(chǎn)生了約束作用。

從城鎮(zhèn)化系統(tǒng)對水資源系統(tǒng)的響應來看(表3)： ①空間城鎮(zhèn)化與水資源系統(tǒng)擬合度最高，其次是經(jīng)濟城鎮(zhèn)化、社會城鎮(zhèn)化和人口城鎮(zhèn)化。表明哈爾濱市水資源系統(tǒng)與社會—經(jīng)濟—空間的城鎮(zhèn)化關(guān)聯(lián)性較強，而與人口城鎮(zhèn)化的關(guān)聯(lián)相對較弱； ②城鎮(zhèn)化系統(tǒng)及其子系統(tǒng)對水資源系統(tǒng)的響應度數(shù)值較為接近，極值差僅為0.288，說明響應的敏感性差異較小，反映出生態(tài)文明背景下哈爾濱的城鎮(zhèn)化發(fā)展追求生產(chǎn)、生活、生態(tài)的全面效益，進而優(yōu)化了水資源系統(tǒng)結(jié)構(gòu)； ③從城鎮(zhèn)化和水資源各類指數(shù)的響應度數(shù)值來看(表2—3)，前者明顯大于后者，表明哈爾濱市城鎮(zhèn)化發(fā)展對水資源系統(tǒng)的變動更具敏感性，水資源變動對城鎮(zhèn)化發(fā)展的影響有放大效應。尤其是哈爾濱城鎮(zhèn)化粗放式快速增長時期，大量新增的用水需求對水資源系統(tǒng)形成脅迫作用，增大了水資源系統(tǒng)衰退的風險。

圖3 水資源系統(tǒng)與城鎮(zhèn)化系統(tǒng)發(fā)展的關(guān)系

2.3.3 水資源子系統(tǒng)和城鎮(zhèn)化子系統(tǒng)的相關(guān)性分析 對哈爾濱水資源子系統(tǒng)評價指數(shù)和城鎮(zhèn)化子系統(tǒng)評價指數(shù)進行相關(guān)性分析，形成的4×4關(guān)系矩陣可反映出子系統(tǒng)間的關(guān)系(見表4)。①水資源荷載與經(jīng)濟城鎮(zhèn)化呈“顯著”的負相關(guān)，與空間城鎮(zhèn)化呈“非常顯著”的負相關(guān)，表明哈爾濱市水資源的負荷壓力與經(jīng)濟活動和空間擴張關(guān)聯(lián)較大。哈爾濱新區(qū)建設和城鎮(zhèn)空間拓展，導致水需求加大，加劇了水資源負荷；②水資源利用、水資源管理與經(jīng)濟、社會和空間城鎮(zhèn)化都存在“非常顯著”的正相關(guān)，說明城鎮(zhèn)化的發(fā)展加強了水資源開發(fā)利用程度、提升了水資源管理能力，從而提高水資源利用效率； ③水資源稟賦與城鎮(zhèn)化各子系統(tǒng)均不相關(guān)，體現(xiàn)出水資源稟賦的自然屬性與城鎮(zhèn)化發(fā)展的社會經(jīng)濟屬性具有一定的獨立性。

表4 哈爾濱市水資源子系統(tǒng)與城鎮(zhèn)化子系統(tǒng)的相關(guān)關(guān)系

注：**表示p<0.01，差異非常顯著； *表示p<0.05,差異顯著。

3 結(jié)論與建議

3.1 結(jié) 論

(1) 從發(fā)展過程看，2006—2017年哈爾濱市水資源開發(fā)利用與城鎮(zhèn)化發(fā)展均呈階段性增長趨勢。水資源開發(fā)利用過程變化劇烈，城鎮(zhèn)化發(fā)展過程相對平穩(wěn)。二者均以2011年前后為時間節(jié)點，表現(xiàn)出前期波動、后期趨穩(wěn)的階段性特征。

(2) 從系統(tǒng)構(gòu)成看：①哈爾濱市城鎮(zhèn)化系統(tǒng)中經(jīng)濟、空間和社會3個子系統(tǒng)的上升發(fā)展較為同步。哈爾濱市城鎮(zhèn)化由2014年以前單一的人口城鎮(zhèn)化驅(qū)動轉(zhuǎn)變?yōu)?014年以后的經(jīng)濟、空間和社會城鎮(zhèn)化共同驅(qū)動。空間城鎮(zhèn)化的發(fā)展表明，在人地協(xié)調(diào)和均衡開發(fā)的前提下，哈爾濱市的外圍縣(市)、鎮(zhèn)尚有較大的城鎮(zhèn)化提升空間。人口城鎮(zhèn)化在研究期內(nèi)起伏較大，雖然發(fā)展?jié)u趨平穩(wěn)，但是需要警惕持續(xù)人口外流帶來系統(tǒng)衰退的風險；②水資源系統(tǒng)中，受哈爾濱市自然條件影響水資源稟賦表現(xiàn)出明顯的不穩(wěn)定性。盡管水資源利用和水資源管理整體上升發(fā)展，但水資源負荷仍有少量增加，主要原因在于哈爾濱市農(nóng)業(yè)耗水居高不下抵消了工業(yè)節(jié)水和其他節(jié)水的效果。

(3) 從響應關(guān)系來看：①城鎮(zhèn)化系統(tǒng)對水資源系統(tǒng)的響應強度明顯高于水資源系統(tǒng)對城鎮(zhèn)化系統(tǒng)的響應程度。表明在哈爾濱市這類缺水城市中，城鎮(zhèn)化發(fā)展對水資源的變動具有更強的敏感性，水資源的量和質(zhì)是制約城鎮(zhèn)化發(fā)展的重要因素；②空間、經(jīng)濟和社會城鎮(zhèn)化對水資源系統(tǒng)響應較強，說明哈爾濱市的經(jīng)濟集聚和空間擴張產(chǎn)生了明顯的資源消耗效應。人口城鎮(zhèn)化對水資源系統(tǒng)的響應較弱，原因在于受人口外流影響，哈爾濱市人口城鎮(zhèn)化的發(fā)展階段與水資源系統(tǒng)改善并不具有對應關(guān)系；③水資源荷載與經(jīng)濟城鎮(zhèn)化和空間城鎮(zhèn)化呈明顯的負相關(guān)，表明哈爾濱市水資源的負荷壓力與經(jīng)濟活動和空間擴張關(guān)聯(lián)較大，水資源開發(fā)利用受到城鎮(zhèn)化發(fā)展的脅迫，城鎮(zhèn)化對資源環(huán)境的剝奪效應明顯。水資源利用、水資源管理與經(jīng)濟、社會和空間城鎮(zhèn)化都存在“非常顯著”的正相關(guān)，說明城鎮(zhèn)化的發(fā)展能夠加強水資源開發(fā)利用程度、提升水資源管理能力，促進水資源系統(tǒng)的優(yōu)化。

3.2 建 議

(1) 運用系統(tǒng)觀和整體性思維解決“人—地—水”的復雜矛盾問題。①整合管理資源，協(xié)調(diào)水利、城建、土地和環(huán)保等多個部門探索新的管理機制，處理好哈爾濱市水利規(guī)劃、國土規(guī)劃和城市總體規(guī)劃的關(guān)系；②以水資源時空分布為依據(jù)，適時調(diào)整產(chǎn)業(yè)結(jié)構(gòu)和布局，構(gòu)建節(jié)水型產(chǎn)業(yè)體系，形成“以供定需”和“分質(zhì)供水”的水資源利用格局；③保持城鎮(zhèn)化適度發(fā)展，推進主城區(qū)海綿城市建設，加速水生態(tài)修復和城市修補的“雙修工程”進度，實現(xiàn)水生態(tài)文明和新型城鎮(zhèn)化建設的同步發(fā)展。

(2) 加快城鄉(xiāng)基礎設施信息化建設，增強水資源時空調(diào)控能力和區(qū)域供水保障能力。與“數(shù)字龍江”和“智慧城市”云平臺對接，提升水資源管理的集成化和智能化水平，實現(xiàn)大數(shù)據(jù)驅(qū)動的水資源優(yōu)化配置，實時評估和預測水資源與城鎮(zhèn)化的匹配程度，減少水資源對城鎮(zhèn)化發(fā)展的硬約束。

與社會經(jīng)濟發(fā)展數(shù)據(jù)相比，水資源指標統(tǒng)計類型、統(tǒng)計年限都較為有限，在一定程度上制約了水資源與城鎮(zhèn)化響應關(guān)系的規(guī)律性總結(jié)。并且由于樣本數(shù)的限制，基于“點彈性”的響應模型在變量擬合的精度上受到一定限制。此外，哈爾濱市建成區(qū)與外圍縣(市)鄉(xiāng)在經(jīng)濟發(fā)展和資源環(huán)境方面差異較大，因此市域內(nèi)部差異帶來的復雜影響、人口城鎮(zhèn)化與水資源開發(fā)利用的作用機制等問題的解決還有待于統(tǒng)計數(shù)據(jù)的積累和測度模型的進一步完善。