劉家宏，周晉軍，邵薇薇

(1.中國水利水電科學研究院流域水循環(huán)模擬與調(diào)控國家重點實驗室，北京 100038; 2.水利部水資源與水生態(tài)工程技術(shù)研究中心，北京 100038; 3.清華大學水利水電工程系，北京 100084)

隨著城市化進程的不斷推進，城市面積增大，人口增多，用水量也持續(xù)增加[1-2]。城市耗水一般包括城市面積的蒸散發(fā)、入滲補給[3]、工業(yè)耗水[4]和產(chǎn)品輸出所包含水量等，是城市水文循環(huán)的重要環(huán)節(jié)。城市耗水量的大小決定了城市區(qū)域水循環(huán)通量的大小，同時直接影響城市供需水的規(guī)劃管理，其研究緊迫而重要。

由于城市下墊面的復(fù)雜性和強人類活動干擾，城市的水循環(huán)過程相比流域自然水循環(huán)要復(fù)雜[5]。城市耗水是指發(fā)生在城市區(qū)域的各種形式的水的耗散，是重要的城市水循環(huán)過程[6]，既包括植被蒸騰蒸發(fā)，也包括城市建筑物內(nèi)部人類用水產(chǎn)生的水分耗散。城市耗水量的大小對空氣濕度影響較大。承德市監(jiān)測數(shù)據(jù)表明，過去50年的城市化過程中，市區(qū)空氣濕度明顯增加[7]。蒸散發(fā)是典型的自然側(cè)耗水形式，城市區(qū)域地表蒸散發(fā)量通?；谀芰科胶夥匠?，通過遙感反演土地利用進行計算[8-9]。城市化過程對城市區(qū)域的蒸散發(fā)也存在影響，其中人為熱源釋放對城市蒸散發(fā)影響較大[10]。另一方面，城市生活用水量的比例增大。2014年中國城市統(tǒng)計年鑒數(shù)據(jù)表明，全國城市生活用水量占城市用水總量的1/3以上。隨著城鎮(zhèn)化進程加快，生活用水量比例增大，由此產(chǎn)生的耗水量也增大，水文效應(yīng)持續(xù)增強[11]。有學者研究了居民生活用水效率[12]、生活用水指標[13]等方面內(nèi)容。除了居民生活用水，城市建筑物內(nèi)部用水也包括生產(chǎn)、室內(nèi)種植等其他用水，這些建筑物內(nèi)部用水同樣會產(chǎn)生水汽耗散，這部分耗水也是城市耗水的組成部分。在城市蒸散發(fā)計算方面，郭偉等[14]以佳木斯為例，將城市下墊面劃分為水面、不透水面和透水面進行蒸發(fā)計算。Ramamurthy等[15]研究城市不透水地面的貢獻率。在城市耗水方面，Qin等[16]采用用水定額、耗水系數(shù)和水量平衡法計算城市耗水量，周晉軍等[6]建立了基于不同耗水特征的分布式城市耗水計算模型。由于城市下墊面的高度異質(zhì)性和人類活動的干擾，城市區(qū)域的耗水過程呈現(xiàn)多樣性，耗水機理較為復(fù)雜。城市耗水作為城市水文循環(huán)的重要過程，耗水機理解析和耗水量計算，對城市水循環(huán)機理和水循環(huán)通量研究具有重要意義，因此，解析城市耗水機理，建立具有物理機制的城市耗水計算模型，既是城市水文研究的難點，也是的熱點。

1 城市高耗水現(xiàn)象

為描述城市高耗水現(xiàn)象，選取典型工廠、建筑物和小區(qū)3種不同尺度的城市典型下墊面開展用耗水及水循環(huán)研究，采用水平衡分析、實驗監(jiān)測、水量平衡等方法解析典型城市耗水單元的耗水類型，為探究城市耗水機理做基礎(chǔ)。

1.1 典型工廠的水平衡測試數(shù)據(jù)分析

選取藍星石化有限公司天津石油化工廠作為水平衡分析對象，通過解析其“供-用-耗-排-循環(huán)再生-回用”過程，對其各個生產(chǎn)環(huán)節(jié)的耗水量進行分析計算，從物理、化學機制上分析其水分的運移轉(zhuǎn)化規(guī)律，查清高耗水的原因和水分的去向。重點分析其產(chǎn)品耗水、廠區(qū)辦公/生活耗水、廠區(qū)綠化耗水等，計算工廠耗水總量。對產(chǎn)品攜帶走的“自由態(tài)”水、“結(jié)晶態(tài)”水和“化合態(tài)”水進行甄別計算，揭示工廠區(qū)超高強度耗水的內(nèi)在成因。分析結(jié)果顯示：藍星石化的新水取用量是3 843.55 m3/d，耗水量是2 242.84 m3/d，排水量1 600.77 m3/d。折合成年新水取用量為140.3萬m3/a，耗水量81.9萬 m3/a，排水量58.4萬 m3/a。按照工廠占地面積78.33 萬m2計算，藍星石化的人工取用水耗水強度為1 046 mm，疊加上降水造成的蒸發(fā)耗水量，藍星石化的單位面積的耗水量要遠大于當?shù)氐慕邓浚撬Y源的凈消耗區(qū)[17]。

1.2 典型城市建筑物耗水監(jiān)測研究

選取辦公樓和住宅公寓作為典型城市建筑物進行耗水監(jiān)測研究，表1給出了典型建筑物基本信息。通過用水量、排水量的監(jiān)測，結(jié)合各用水活動水量耗散的實驗監(jiān)測，解析典型住宅樓和辦公樓的耗水過程，計算耗水量。建筑物用水量通過安裝在自來水供水管的水表計量，排水量通過中國水利水電科學研究院自主研發(fā)的間歇量水設(shè)備(專利號碼：CN104614030A)進行監(jiān)測。需要說明的是，所選兩個建筑物的供水水源均為自來水。建筑物排水監(jiān)測設(shè)備安裝在建筑物內(nèi)部排水管末端，忽略建筑物內(nèi)部供水量檢測與排水監(jiān)測中的管道滲漏損失，認為用水量減去排水量等于建筑物內(nèi)部耗水量。

表1 典型城市建筑物基本信息

圖1 辦公樓排水系數(shù)測定

圖1為典型辦公樓用水-排水過程曲線。從圖1可以看出，在單位小時時間尺度上，用水量在0.05～0.94 m3之間波動，但是排水量始終低于用水量，并且排水系數(shù)(排水量與用水量之比)在0.76～0.95間波動。根據(jù)前面的實驗條件，可以得出建筑物對應(yīng)的耗水系數(shù)變化區(qū)間為0.05～0.24，說明建筑物內(nèi)部耗水現(xiàn)象存在，并且具有相對穩(wěn)定的變化特征。

分別對典型建筑物進行2個月的供水和排水監(jiān)測，計算得出辦公建筑物的平均排水系數(shù)為0.89，住宅樓的平均排水系數(shù)為0.83。測量得到實驗期內(nèi)典型辦公樓和住宅樓的日平均用水量分別是5.5 m3和10.6 m3。結(jié)合表1信息，知實驗建筑物單位建筑面積年耗水強度為辦公樓37 mm，住宅樓 183 mm，由于典型建筑物樓層均為3層，因此對應(yīng)的單位地面面積的耗水強度分別是111 mm和549 mm。由此說明城市辦公樓和住宅樓屬于高耗水下墊面，建筑物耗水量與建筑物樓層成正比。

1.3 典型小區(qū)水循環(huán)觀測試驗

選取典型小區(qū)作為區(qū)域典型代表，以2014年數(shù)據(jù)為例進行區(qū)域耗水現(xiàn)象分析。耗水計算方法選用水量平衡方程進行：

WD=WS+P-D-WSUR-ΔWG

(1)

式中：WD為區(qū)域年耗水總量，m3；WS為年供水總量，m3；P為年降雨總量，m3；D為區(qū)域年排水總量，m3；WSUR是地表徑流總量，m3；ΔWG是地下水變化量，m3。

所選區(qū)域面積是3.5 km2，居民60 000人，區(qū)域內(nèi)綠化率55%，氣象站區(qū)域內(nèi)測得2014年降水量為442 mm。區(qū)域具有相對獨立的供排水系統(tǒng)，由于地勢較高，周邊地區(qū)的降雨徑流不會流入該小區(qū)，因此地表徑流量數(shù)據(jù)與區(qū)域面積相對應(yīng)。監(jiān)測數(shù)據(jù)顯示該區(qū)域2014年總供水量是432萬m3，總降雨量是155萬m3，區(qū)域的污水排放量是254萬m3，地表徑流總量是55萬m3，地下水變化量近似認為是0，由此計算得出區(qū)域2014年單位面積耗水量是794 mm，高于遙感反演得到的蒸散發(fā)520 mm，說明區(qū)域耗水強度遠高于傳統(tǒng)蒸散發(fā)計算結(jié)果。

周琳[18]以北京市為例，將北京市劃分為山區(qū)和平原，其中平原又分為城區(qū)和郊區(qū)，利用水量平衡原理研究北京市城市蒸散發(fā)，結(jié)果表明，城市不同分區(qū)的同期蒸散發(fā)量差異很大，2002—2013年計算結(jié)果顯示城區(qū)平均蒸散發(fā)量最大，599 mm，山區(qū)最小，479 mm。城區(qū)蒸發(fā)量大的主要原因是區(qū)域內(nèi)高密度的建筑和人類活動。

2 城市高耗水機理分析

典型城市單元耗水現(xiàn)象研究表明，城市的耗水發(fā)生在城市的各個用水部位，除了城市綠地、水面等自然側(cè)蒸散發(fā)，還包括人類用水所產(chǎn)生的社會側(cè)蒸散發(fā)，特別是建筑體內(nèi)部人類活動產(chǎn)生的水量耗散，而這部分耗水量在以往城市蒸散發(fā)研究中是被忽略的，因此將城市的耗水劃分為建筑體內(nèi)部耗水和建筑體外部耗水(即露天耗水，包括硬化地面的耗水和綠地、水面等自然側(cè)的下墊面蒸散發(fā))。

2.1 建筑物內(nèi)部耗水

建筑物內(nèi)部的耗水主要是水由液態(tài)轉(zhuǎn)化為氣態(tài)形式，然后隨著空氣流通，經(jīng)由建筑體的窗戶、門等部位(類似于植被葉片的氣孔)進入室外的空氣中，參與水循環(huán)過程。根據(jù)建筑物分類標準，可以將建筑物分為生產(chǎn)性建筑物和非生產(chǎn)性建筑物。生產(chǎn)性建筑物包括工業(yè)建筑物和農(nóng)業(yè)建筑物。這些建筑物在城市建成區(qū)內(nèi)分布較少，大多分布在城郊，因此在本研究中不考慮生產(chǎn)性建筑物。非生產(chǎn)性建筑物也稱為民用建筑，相比生產(chǎn)性建筑物，民用建筑廣泛分布在城市建成區(qū)域內(nèi)，其耗水是有限的。不同功能的建筑物包含的耗水項不同，耗水特征也不同。相同功能的建筑物耗水量與人數(shù)、建筑面積相關(guān)。本文將城市民用建筑劃分為2級7種建筑大類，并與耗水類型建立關(guān)聯(lián)關(guān)系，構(gòu)建建筑物內(nèi)部耗水模型框架(圖2)。

圖2 建筑物內(nèi)部耗水結(jié)構(gòu)

基于仿生學原理，建筑物可以類比為城市的“建筑樹”組成的“混凝土森林”[2]，每一個建筑物相當于“森林”中的一棵“樹”，通過類似于樹木的葉片氣孔的門、窗、通風口等部位將建筑物的水汽釋放到空氣中，參與城市區(qū)域的大氣水循環(huán)。由于城市建筑物多數(shù)為多層建筑，建筑物的耗水在垂直空間上也呈均勻分層特征，對于上下均勻的多層建筑物，可以認為不同樓層的耗水特征相同，建筑物投影面積的耗水強度與樓層數(shù)成正比。

2.2 城市露天耗水

城市露天耗水可以根據(jù)二元水循環(huán)原理劃分為自然側(cè)耗水和社會側(cè)耗水。自然側(cè)耗水主要指綠地、土壤、植被、水面等自然形式的下墊面的蒸騰蒸發(fā)量以及入滲補給；社會側(cè)耗水主要指人工澆灌植被、綠地的耗水，市政道路灑水，硬化地面截留雨水蒸發(fā)的耗水等?；诤乃卣骺梢詫⒊鞘邢聣|面劃分為建筑物屋頂、不透水硬化地面、透水硬化地面、綠地、水面、裸土5類。

建筑物屋頂、不透水硬化地面的耗水形式主要是截留雨水蒸發(fā)，其中部分城市道路的耗水來源還有市政灑水，這部分水量也全部蒸發(fā)到空氣中。由于建筑屋頂和不透水硬化地面的止水能力有限，同時與下層結(jié)構(gòu)水分隔絕，因此蒸發(fā)呈現(xiàn)間歇性，在有降雨或者灑水發(fā)生時表面截留水量會有蒸發(fā)現(xiàn)象發(fā)生[15]。對于城市道路，汽車車輪“霧化效應(yīng)”增強了蒸發(fā)的動力項，提高了蒸發(fā)效率。

透水硬化地面主要包括部分鋪磚人行道和海綿城市建設(shè)中推行的透水鋪裝(透水磚、透水混凝土)[19]。與不透水地面相比，透水鋪裝增加了滲透的耗水方式，耗水量由蒸發(fā)和下滲共同組成。對于一場降雨，降雨量減去產(chǎn)流量即為耗水量。對于一次降雨，忽略期間的蒸發(fā)量。在月尺度、年尺度或者更長時間序列中，耗水量由蒸發(fā)量和深層滲漏水量組成。

綠地、水面和裸土的蒸散發(fā)過程，認為與自然狀態(tài)的過程基本一致，也稱之為城市自然側(cè)耗水。不同的是，城市區(qū)域的綠地耗水在水源上增加了人工灌溉的水量來源，在一定程度降低了水量限制。同時由于城市人為熱源的釋放以及城市熱島效應(yīng)的影響[12]，城市自然側(cè)的蒸散發(fā)的能量約束力減弱，一定程度上了增加了蒸散發(fā)的動力項。

3 城市耗水計算模型

圖3 城市建成區(qū)耗水計算方法

基于城市水循環(huán)的“供水-用水-耗水-排水-循環(huán)再生-回用”過程，耗水量及耗水機理是較為核心的內(nèi)容。筆者以耗水量的計算作為切入點，建立不同類型城市斑塊耗水計算的理論范式，采用“自上而下”和“自下而上”兩種思路進行耗水量計算(圖3)，并在過程中搜集、計算其他水循環(huán)過程數(shù)據(jù)，系統(tǒng)解析了典型城市單元水循環(huán)的結(jié)構(gòu)和路徑，最終完成了城市耗水計算理論框架的構(gòu)建。 “自上而下”，即將耗水量作為輸出水量的未知項，再分別計算區(qū)域水量的總輸入和輸出，通過水量平衡關(guān)系計算耗水量；“自下而上”，即將耗水量按下墊面類型進行分類計算，如計算水面蒸發(fā)、綠地蒸散發(fā)、硬化地面蒸發(fā)、建筑耗水等，最后通過不同面積上的耗水積分得到總的耗水量。

該耗水計算方法是適用于城市不同類型斑塊上的“雙源”(天然降水和人工供水)耗水定量計算模型。模型以模擬城市耗水過程為基礎(chǔ)，進行城市耗水單元的劃分，模型突出耗水發(fā)生的界面和產(chǎn)生的物理機制，主要包括硬化地面及建筑物截留蒸發(fā)模塊、馬路灑水蒸發(fā)模塊、各建筑體(住宅、辦公、商場、工廠等)內(nèi)用水的蒸發(fā)耗散計算模塊、城市綠地蒸發(fā)/蒸騰計算模塊、城市水面蒸發(fā)計算模塊等。

基于上述城市耗水機理和耗水計算方法，對模型進行概化，基于“自下而上”思路，建立城市區(qū)域耗水計算模型：

WD=0.365PNPDI+(AHEH+AGEG+AWEW+ASES)/1000

(2)

式中：PN為區(qū)域內(nèi)常住人口量，人；PDI為人均日耗水量，L/d；AH為區(qū)域硬化地面面積，m2；EH為硬化地面年蒸散發(fā)強度，mm；AG為區(qū)域綠地面積，m2；EG為綠地年蒸散發(fā)強度，mm；AW為區(qū)域水面面積，m2；EW為水面年蒸散發(fā)強度，mm；AS為區(qū)域裸土地面面積，m2；ES為裸土地面年蒸散發(fā)強度，mm。

同樣地，基于北京市的土地利用類型建立覆蓋北京市城區(qū)和郊區(qū)的耗水計算模型，根據(jù)2014年TM影像，30 m分辨率(云覆蓋較多的用2013年補充，其余選擇5—9月份遙感影像)解譯制作，主要土地利用類型包括耕地、林地、草地、水域、城鄉(xiāng)和工礦建設(shè)用地，以及未利用土地6種類型。其中城鄉(xiāng)和工礦建設(shè)用地耗水量通過式(2)計算得到，其他土地利用類型的蒸散發(fā)通過彭曼公式計算，2014年耗水計算結(jié)果見圖4中。需要說明的是，本研究沒有計算北京市工業(yè)的耗水量，在城鄉(xiāng)和工礦用地中沒有考慮人工灑水產(chǎn)生的蒸散發(fā)。耗水空間分布格局表明人類活動和用水量集中的城市區(qū)域耗水量高于郊區(qū)。圖4中紅色區(qū)域為北京市建成區(qū)主要區(qū)域，說明城市區(qū)域人類用水產(chǎn)生的耗水量不容忽視，也證實了城市是高耗水區(qū)域的判斷。

圖4 北京市2014年耗水空間格局

4 結(jié)論和展望

以不同尺度的城市單元為例，描述了城市高耗水現(xiàn)象，分析了城市耗水過程，解析了城市耗水機理，提出了城市耗水計算方法，得出了以下主要結(jié)論：①城市高耗水現(xiàn)象確實存在，其主要原因是受人類用水活動的影響，建筑物內(nèi)部耗水是其高耗水的主要增加項；②城市耗水是二元水循環(huán)的重要環(huán)節(jié)，耗水量的大小不僅影響城市供排水管網(wǎng)規(guī)劃設(shè)計，同時影響城市區(qū)域的水熱平衡關(guān)系；③城市的耗水可以分為建筑物內(nèi)部耗水和露天耗水，基于耗水特征差異，露天耗水可以劃分為建筑屋頂與不透水硬化地面、透水硬化地面、綠地、水面和裸土耗水；④考慮人類用水過程中產(chǎn)生的耗水，北京市城區(qū)的耗水量高于郊區(qū)，同時高于林地、草地、耕地以及未利用土地。

本研究基于水平衡數(shù)據(jù)分析，以及典型小區(qū)、典型耗水過程的觀測實驗，初步闡釋了城市耗水的機理，建立了相關(guān)的計算模型。但研究成果還是初步的，需要在未來的研究中進一步深化，同時城市耗水強度與城市局部的水熱平衡關(guān)系密切，與城市熱島效應(yīng)等的內(nèi)在聯(lián)系和相互作用機理也有待進一步研究揭示。

參考文獻：

[1] MCDONALD R I,DOUGLAS I,REVENGA C,et al.Global urban growth and the geography of water availability,quality,and delivery[J].Ambio,2011,40(5):437.

[2] 胡慶芳,王銀堂,李伶杰,等.水生態(tài)文明城市與海綿城市的初步比較[J].水資源保護,2017,33(5):13-18.(HU Qingfang,WANG Yintang,LI Lingjie,et al.Preliminary comparison between water-ecological civilization city and sponge city[J].Water Resources Protection,2017,33(5):13-18.(in Chinese))

[3] 索海生.文峪河灌區(qū)開柵高耗水區(qū)節(jié)水改造噴灌工程的分析[J].水資源保護,2005,21(1):69-71.(SUO Haisheng.Analysis on the water saving modified spray irrigation project in high water consumption area of the Wenyu River irrigation[J].Water Resources Protection,2005,21(1):69-71.(in Chinese))

[4] 何月峰,沈海萍,馮曉飛,等.基于壓力-狀態(tài)-響應(yīng)模型和"五水共治"決策的浙江省水環(huán)境安全評價[J].水資源保護,2016,32(6):104-109.(HE Yuefeng,SHEN Haiping,FENG Xiaofei,et al.Assessment of water environmental security in Zhejiang Province based on PSR model and policy of co-governance of fiver water categories[J].Water Resources Protection,2016,32(6):104-109.(in Chinese))

[5] 黃國如,王欣.基于城市雨洪模型的市政排水與水利排澇標準銜接研究[J].水資源保護,2017,33(2):1-5.(HUANG Guoru,WANG Xin.Assessment of water environmental security in Zhejiang Province based on PSR model and policy of co-governance of fiver water categories[J].Water Resources Protection,2017,33(2):1-5.(in Chinese))

[6] 周晉軍,劉家宏,董慶珊,等.城市耗水計算模型[J].水科學進展,2017,28(2):276-284.( ZHOU Jinjun,LIU Jiahong,DONG Qingshan,et al.Simulation model for urban water dissipation[J].Advances in Water Science,2017,28(2):276-284.(in Chinese))

[7] 馬鳳蓮,丁力,王宏.承德市干濕島效應(yīng)及其城市化影響分析[J].氣象與環(huán)境學報,2009,25(3):14-18.(MA Fenglian,DING Li,WANG Hong.Dry and moisture island effects and its responses to the urbanization in Chengde,Hebei province[J].Journal of Meteorology & Environment,2009,25(3):14-18.(in Chinese))

[8] 敬書珍,倪廣恒.基于遙感的城市蒸散發(fā)對土地利用/覆蓋的響應(yīng)研究[C]// 2008年全國城市水利學術(shù)研討會暨工作年會資料.北京：中國水利學會，2008.

[9] 吳炳方,邵建華.遙感估算蒸騰蒸發(fā)量的時空尺度推演方法及應(yīng)用[J].水利學報,2006,37(3):286-292.(WU Bingfang,SHAO Jianhua.Temporal and spatial extension of evaportranspiration estimated from remote sensing[J].Journal of Hydraulic Engineering,2006,37(3):286-292.(in Chinese))

[10] 潘文祥,高學睿,王玉寶,等.城市家庭生活單元耗用水過程和機理研究進展[J].南水北調(diào)與水利科技,2017,15(2):140-147.(PANG Wenxiang,GAO Xuerui,WANG Yubao,et al.Research progress on the process and mechanism of urban household water consumption[J].South-to-North Water Transfers and Water Science & Technology,2017,15(2):140-147.(in Chinese))

[11] 方詩標,賈仁甫,方詩彬,等.城市生活用水效率驅(qū)動因子及調(diào)控研究[J].人民黃河,2013,35(3):47-50.(FANG Shibiao,JIA Renpu,FANG Shibin,et al.Study on drive factors and regulation of urban domestic water use efficiency[J].Yellow River,2013,35(3):47-50.(in Chinese))

[12] CONG Z T,SHEN Q N,ZHOU L,et al.Evapotranspiration estimation considering anthropogenic heat based on remote sensing in urban area[J].Science China Earth Sciences,2017,60(4):659-671.

[13] 劉家宏,王建華,李海紅,等.城市生活用水指標計算模型[J].水利學報,2013,44(10):1158-1164.(LIU Jiahong,WANG Jiahua,LI Haihong,et al.A mathematic model for rational domestic water demand considering climate and economic development factors[J].Journal of Hydraulic Engineering,2013,44(10):1158-1164.(in Chinese))

[14] 郭偉,劉洪福,張文清.佳木斯市城市蒸散發(fā)量計算[J].黑龍江大學工程學報,1998(3):51-53.(GUO Wei,LIU Hongfu,ZHANG Wenqing.Evaporation calculation for Jiamusi city[J].Journal of Hei longjiang Hydraulic Engineering College,1998(3):51-53.(in Chinese))

[15] RAMAMURTHY P,BOU-ZEID E.Contribution of impervious surfaces to urban evaporation[J].Water Resources Research,2014,50(4):2889-2902.

[16] QIN D Y,LYU J Y,LIU J H,et al.Theories and calculation methods for regional objective ET[J].Chinese Science Bulletin,2009,54(1):150-157.

[17] 劉家宏,王浩,高學睿,等.城市水文學研究綜述[J].科學通報,2014(36):3581-3590.(LIU Jiahong,WANG Hao,GAO Xuerui,et al.Review on urban hydrology[J].Chinese Science Bulletin,2014(36):3581-3590.(in Chinese))

[18] 周琳.北京市城市蒸散發(fā)研究[D].北京：清華大學,2015.

[19] 夏軍,石衛(wèi),王強,等.海綿城市建設(shè)中若干水文學問題的研討[J].水資源保護,2017,33(1):1-8.(XIA Jun,SHI Wei,WANG Qiang,et al.Discussion of several hydrological issues regarding sponge city construction[J].Water Resources Protection,2017,33(1):1-8.(in Chinese))