深圳市重演鄭州“7?20”特大暴雨情景分析

摘要：深圳極易遭遇以“臺風、暴雨、海潮”三碰頭為代表的極端天氣事件，有必要分析其在極端天氣下內澇災害發生、發展和消退全過程。以深圳市羅湖中心區為例，基于分布式水文水動力耦合模型，采用暴雨移置方法，重演了鄭州“7·20”特大暴雨造成的內澇災害風險，分析了羅湖中心區現狀工程體系對極端暴雨的響應。結果表明：深圳市羅湖中心區內澇積水深度大于0.3 m的區域面積占比高達13.58%，積水范圍與小于2 a一遇設計重現期的雨水管渠分布范圍高度重合。現狀內澇成因主要為管渠輸水能力不足，澇水不能迅速通過管網輸送至泵站前池，從而造成遠端積澇。研究成果有助于提高深圳市極端暴雨災害認識和應對能力，可為超標準洪水應急預案編制和汛前演練等提供技術支撐。

關鍵詞：極端暴雨； 城市內澇； 分布式水文水動力耦合模型； 鄭州“7·20”特大暴雨； 深圳市

中圖法分類號：TV122

文獻標志碼：A

DOI：10.15974/j.cnki.slsdkb.2024.11.001

文章編號：1006-0081（2024）11-0010-07

0 引 言

根據國家統計局2022年末統計數據，中國城鎮化率已經達到65.22%。在經歷數十年的高速發展與城鎮化建設后，中國城市內澇成因及致災機理正在發生一系列深刻改變，加之全球氣候變化和人類活動的加劇，極端暴雨導致洪澇發生的可能性與不確定性增大［1-2］。根據住房和城鄉建設部2010年對351個城市的內澇調研顯示，213個城市發生過不同程度的積水內澇，內澇災害一年超過3次的城市就有137個，其中甚至囊括了干旱少雨的西安、沈陽等北部城市［3］。城鎮水安全問題正處于漸進累積到突變激增的狀態，內澇問題正在成為中國城市水安全保障的重大隱患。

2021年7月17～23日，在西太平洋副熱帶高壓、中國西部的青藏高壓、東海“煙花”臺風和南海“查帕卡”臺風共同作用下［4-5］，河南省遭遇歷史罕見特大暴雨。從過程累計面雨量來看，鶴壁最大（589 mm）、鄭州次之（534 mm）、新鄉第三（512 mm）。從小時最強點雨量來看，鄭州最大，發生在7月20日16∶00～17∶00（鄭州國家氣象站201.9 mm）。據核查評估，河南省共有150個縣（市、區）1 478.6萬人受災，因災死亡398人，直接經濟損失1 200.6億元［6］。

鄭州“7·20”特大暴雨災害為研究特大城市洪澇風險演變特征與趨向提供了實例驗證。王帆等［7］以山東日照市傅疃河流域為研究對象，基于分布式水文模型和二維水動力模型分析了其在鄭州“7·20”特大暴雨情景下的洪水形勢；陳升［8］分別以上海市徐匯區、黃浦區和浦東新區世紀大道為暴雨中心，重演了鄭州“7·20”特大暴雨在上海城區的內澇風險分布；李敏等［9］基于鄭州“7·20”特大暴雨，對武漢市中心城區內澇風險區域進行了定量分析，并提出了應對極端暴雨內澇的策略；張孝奎等［10］以福州濱海核心區遭遇鄭州“7·20”暴雨為例，分析了濱海新城核心區可能面臨的內澇風險；張鴻等［11］將鄭州“7·20”特大暴雨移置到洪汝河流域，分析小洪河流域在“7·20”特大暴雨下的洪水過程。

深圳市羅湖區是深圳市的中心城區之一，地處低緯度濱海臺風頻繁登陸地區，受鋒面雨、臺風雨影響，暴雨頻發，洪、潮災害常常發生，且在氣候變化條件下，近年來“臺風、暴雨、海潮”三碰頭概率增大。為分析羅湖中心區遭遇極端降雨的受災情況，本文以鄭州“7·20”特大暴雨作為輸入條件，構建了深圳市羅湖中心區城市分布式水文水動力耦合模型，分析了其在現狀治澇工程體系、閘站調度方案條件下的內澇災害。

1 研究區概況

1.1 區域基本情況

深圳市羅湖中心區（圖1）為羅湖區的核心地帶，南臨深圳河，與香港新界隔河相望，西至紅嶺路，東抵深圳水庫排洪河，北以泥崗路—布心路—沿河北路為界，分布有深圳火車站、羅湖口岸、文錦渡口岸、東門老商業步行街、國貿商業中心等，是深圳市最大的人流集散地、最早的商業中心和重要的區域性交通樞紐。羅湖中心區地勢總體呈北高南低趨勢，一旦遭遇強降雨，北部澇水部分就近排入布吉河和深圳水庫排洪河，部分迅速向區域南部地勢低洼區域集中。

1.2 治澇現狀

為解決深圳市羅湖中心區和平路、嘉賓路和建設路一帶低洼地段每遇暴雨便內澇成災的現象，深圳市政府于1996年投資建成羅湖中心區防洪治澇工程，將南起深圳河、北至解放路、東起春風路口、西至布吉河6.29 km2范圍的降雨“圍起來、抽出去”，以解決該片區的內澇問題。

羅湖中心區現狀內澇防治工程體系主要包括雨水管渠系統、雨水泵站和水閘。其中，雨水管渠總長90.92 km，集水口447處，外江自排口9處。目前片區基本實現了雨污分流，90%以上的污水管段已與雨水管網分離。雨水管渠中，暗涵總長20.8 km，所占比例高達22.8%。片區共有雨水泵站9座，分別為羅雨、船步街、建設路、蔡屋圍、文錦渡、東廣場、羅湖橋、火車站1號及火車站2號泵站，總設計流量63.7 m3/s，其中羅雨泵站設計流量48 m3/s。片區共有水閘3座，分別為羅雨、文錦渡、布吉河河口水閘，其中羅雨水閘和文錦渡水閘分別為羅雨泵站和文錦渡泵站的配套水閘。

羅湖中心區為雨污合流制，為防止污水直排深圳河、布吉河等，在東門渠、春風渠、嘉賓路箱涵、羅雨干渠等主干雨水管渠內設置了5處截污堰，用于攔截污水下河，分布位置詳見圖2。5處截污堰中堰高最高達0.6 m，占雨水管渠截面積比例最大達42.86%，在一定程度上造成雨水管渠排水能力的下降。

目前羅湖中心區現狀排澇標準整體偏低，內澇防治重現期僅10～20 a。管渠、泵站外排能力嚴重不足，加之受天文大潮和臺風增水疊加影響，深圳河、布吉河水位較高，泵站外排能力不能有效發揮，澇水難以排出；同時，區域澇水僅能通過唯一的羅雨干渠及其附屬管渠進行抽排，外排通道極為單一。

2 模型構建

一個完整的城市排水系統通常包括雨水管渠、內河排水溝渠、排澇泵站等組成的市政排水系統，以及城市內湖、外江等組成的河湖承泄區。單獨將上述某一環節剝離出來建模分析不能準確模擬城市雨洪的發生、發展和消退全過程。本文構建一種降雨產流-地表徑流-管網匯流-泵站外排-河道下泄的城市雨洪耦合模型，包括降雨產匯流模擬、一維河道洪水演進模擬、二維地表漫流模擬和一維管網匯流模擬等模塊。

2.1 降雨產匯流模擬

城市下墊面具有匯水面積較小、地表覆蓋情況復雜、匯流區域邊界不明顯等特點，采用時間-面積曲線法描述城市下墊面的匯流過程，假定匯水區內各點的雨水匯流速度相同，區域出水口流量由每個子匯水流面積上的徑流按匯流至出水口的時間依次線性疊加而成，其計算公式為

Q（ti）=αFtk∑tkj=1Iti-jΔt（1）

式中：Q（ti）為ti時刻雨水口的入匯流量；F為雨水口的集雨面積；Iti-j為ti-j時刻的降雨強度；α為單位換算系數；Δt為時間間隔。為計算匯水區的匯流過程，還需提供匯水區的匯水時間-匯水面積關系。對于任意形狀的匯水區，其匯水時間-匯水面積曲線可由矩形、發散型、匯集型等3種基本的曲線類型通過比例參數組合確定，具體見參考文獻［12］。

2.2 一維河道洪水演進模擬

在假定河道底坡較緩、斷面壓力沿垂線按靜水壓力分布、斷面流速均勻分布等條件下，一維河道水流可由圣維南方程組描述：

BZt+Qx=q（2）

Qt+xQ2A+gAZx+gAn2QQA2R4/3=0（3）

式中：Z和Q分別為明渠水位和流量；A和B分別為明渠斷面過水面積和寬度；q為旁側入流或出流的單寬流量；n為明渠河床糙率系數；R為水力半徑；g為重力加速度；x為沿明渠水流方向的河道里程；t為時間。

一維圣維南方程的數值解法已較成熟，本文采用四點Pressimann隱格式離散方法，詳見參考文獻［13］。

2.3 二維地表漫流模擬

雨水在城市地表的漫流過程可采用二維淺水方程組描述：

ht+huixi=0（4）

huit+huiujxj+xig2h2+ghZbxi=xjhυuixj-τbiρ（5）

式中：xi為沿i方向（x方向或y方向）的平面坐標；h為水深；ui為沿i方向的流速；Zb為河床高程；ρ，υ分別為水的密度和運動黏度；τbi為沿i方向的摩阻，可表示為τbi=ρCbuiujuj，其中Cb為阻力系數，重復下標滿足求和約定。

采用格子Boltzmann方法求解二維淺水方程組，其中離散速度模型采用D2Q9模型，推導及求解步驟詳見參考文獻［14-16］。

2.4 雨水管渠匯流模擬

一維管網水動力學模型控制方程如下：

Qx+At=0（6）

Qt+xαQ2A+gAyx=gAI0-If（7）

式中：Q為管道流量；A為管道過水面積；y為管道中的液面高度；I0為管道底坡；If為阻力坡度；g為重力加速度。一維圣維南方程組僅能模擬重力流，對于城市雨水管網，當降雨超過管道過水能力時（管道超載），其流態將由重力流轉變為壓力流。為了讓圣維南方程組仍然適用于壓力流，可采用Pressimann窄縫假設，在管道頂端引入一段假想的無限高狹槽，狹槽的寬度取決于管道橫斷面面積和波速，具體處理方法詳見參考文獻［17］。

2.5 管渠匯流模型與地表漫流模型的耦合

將一維管渠匯流模型與二維地表匯流模型進行耦合，模擬地表雨水通過雨水口進入管渠和管渠超載后雨水溢出地表的雙向過程，兩者間的水量交換采用孔口流量方程定量計算：

Qex=signHP-HGCDAin 2gHP-HG（8）

式中：Qex為管渠和地表的交換流量；HP為管渠水位；HG為地表水位；Ain為雨水口的過水面積；CD為流量系數。sign（x）為符號函數，當xgt;0時，sign（x）=1；當x=0時，sign（x）=0；當xlt;0時，sign（x）=-1。

3 模型參數率定

考慮到上述模型中參數眾多，逐一對其取值進行率定既無必要也不現實。為驗證采用本模型進行后續情景模擬的合理性，各參數取值除參考相關行業規范和以往研究成果外，本節采用2018年8月29日暴雨期間，深圳市羅湖中心區實測內澇積水情況進行率定，率定結果詳見表1。可以看出，2018年實測積澇點位置、范圍、水深、歷時等參數與模型計算結果大致吻合。需要說明的是，因為文錦渡口岸關卡限制勘測人員進入，因此未收集到沿河南路與深圳河干流右岸之間的管網資料，所以片區南部靠近深圳河干流的兩處實測積澇點在模擬結果中未能得到反映。

4 重演情景設定

4.1 降雨邊界條件設定

本次采用2021年鄭州“7·20”特大暴雨鄭州國家氣象站最大日降雨過程線，詳見圖3，其中最大小時降雨量高達201.9 mm。羅湖中心區所在深圳河流域200 a一遇1 h設計雨量為152 mm，24 h設計雨量為610 mm，可見鄭州“7·20”特大暴雨實際降雨過程在深圳市均超過了200 a一遇標準。

4.2 外海潮汐邊界條件設定

除雨水泵站外，羅湖中心區有大量雨水排口直排深圳河和布吉河，其外排流量受管渠上游來水和外江水位雙重影響，計算時需要通過一維洪水演進模型提供深圳河和布吉河水位作為其邊界條件。深圳河下游直排南海，屬感潮河流，采用《深圳市防洪（潮）排澇規劃（2021～2035）》擬定的洪潮組合，對于以雨為主的洪潮遭遇，推薦采用設計頻率降雨遭遇5 a一遇的高潮位過程，詳見圖4。

5 情景模擬結果分析

5.1 內澇積水發展過程

遭遇鄭州“7·20”特大暴雨最大24 h降雨情況下，深圳市羅湖中心區不同時刻積水分布詳見圖5。現狀條件下，羅湖中心區內澇防治標準為10～20 a一遇，最大可抵御約90 mm的小時降雨量。在7月19日16∶00之前，本場降雨過程最大小時降雨量未超過61 mm，區域未出現內澇災害。7月19日16∶00的小時降雨量達201.9 mm，開始大面積出現內澇災害，災害范圍與小于2 a一遇設計重現期的雨水管渠分布范圍高度重合。鄭州“7·20”特大暴雨衰減較快，加之片區面積較小、地表不透水率較高，退水時間較快，7月20日23∶00全域積水深度基本退至50 cm以下。深圳羅湖區東門老商業街和區政府大樓等典型位置積水過程線見圖6，東門老商業街最大積水深度達0.64 m，積水時間53 min；區政府大樓最大積水深度達1.96 m，積水時間長達90 min。

本場降雨過程中，積水深度大于1.5 m的區域面積為0.243 km2，占羅湖中心區總面積的3.86%；積水深度1.0～1.5 m的區域面積為0.341 km2，占羅湖中心區總面積的5.42%；積水深度0.5～1.0 m的區域面積為0.252 km2，占羅湖中心區總面積的4.01%。內澇積水深度大于0.3 m的區域面積占比高達13.58%，積水范圍與小于2 a一遇設計重現期的雨水管渠分布范圍高度重合。

5.2 治澇工程運行情況

作為羅湖中心區治澇的關鍵性工程，羅雨閘站承擔了片區80%以上的雨水外排任務，其中泵站共配備了6臺軸流泵，總設計流量48 m3/s；水閘共3孔，單孔尺寸6 m（寬）×3 m（高），設計流量150 m3/s。深圳火車站下沉廣場因地勢較低，是片區易澇點之一，其排澇任務由東廣場排澇泵站承擔，配備2臺潛水泵，總設計流量0.72 m3/s。本次降雨過程中，羅雨閘站與東廣場泵站外排過程詳見圖7。羅雨泵站最多僅開啟2臺機組，因羅雨泵站服務面積較大，雨水匯流時間較長，泵站開啟時間相對于雨峰位置有所滯后，首臺機組于雨峰后25 min開啟，開啟后未能將泵站前池水位維持在1.8 m的控制水位以下，第二臺機組于雨峰后31 min開啟。羅雨水閘最大泄流量為71.47 m3/s。因服務面積較小，東廣場泵站于雨峰后9 min開啟，并持續開啟了105 min。

即便遭遇鄭州“7·20”特大暴雨，區域最大的羅雨泵站僅開啟6臺泵站中的2臺即可將泵站前池水位降至控制水位以下，而區域依然嚴重積澇。可見片區形成內澇的重要原因是遠端管渠輸水能力不足，澇水不能迅速通過管渠輸送至泵站前池，進而造成遠端積澇。同時，鄭州“7·20”特大暴雨降雨強度遠超現狀排水管渠排水能力，也是片區形成內澇的原因之一。

為分析截污堰對片區排澇的影響，以東門渠1號截污堰為例，分析阻水截污堰對管渠排水能力的影響。東門渠1號截污堰位于東門中路西側雨水箱涵，該處雨水管渠現狀為1.6 m（寬）×1.6 m（高）矩形箱涵，截污堰尺寸為1.6 m（寬）×0.3 m（高），占雨水管渠過水面積的18.75%。鄭州“7·20”特大暴雨條件下，東門渠1號截污堰上、下游管段最高水頭線如圖8所示。可見，截污堰上、下水頭差約為0.76 m；在羅雨干渠接納東門渠位置，羅雨干渠處于超載狀態（水頭線高于管渠頂部高程），內部水流狀態為壓力流。在羅雨干渠高水頂托和截污堰阻水的雙重影響下，本段東門渠整體處于超載狀態。圖9為鄭州“7·20”特大暴雨情況下，東門渠1號截污堰上、下游水頭過程對比情況。拆除門渠1號截污堰后，堰上水頭最大降低約0.73 m，可顯著改善截污堰上游管渠排水條件。可見，截污堰雖然可攔截部分污水，對承泄區水質改善起到一定積極作用，但其占用了部分雨水管渠過水斷面，較大程度上降低了雨水管渠的排水能力，導致澇水不能及時通過排水管渠輸送至末端泵站。

6 結 論

本文在構建城市分布式水文水動力耦合模型基礎上，計算分析了深圳市羅湖中心區重演鄭州“7·20”特大暴雨的內澇災害，主要結論如下。

（1） 鄭州“7·20”特大暴雨條件下，深圳市羅湖中心區積水深度大于1.5 m的區域占總面積的3.86%；積水深度1.0～1.5 m的區域占總面積的5.42%。深圳市羅湖中心區內澇積水深度大于0.3 m的區域面積占比高達13.58%，積水范圍與小于2 a一遇設計重現期的雨水管渠分布范圍高度重合。

（2） 羅雨閘站作為片區治澇的骨干工程，該場降雨僅開啟6臺泵站中的2臺，說明形成內澇的主要原因為現狀管渠排水標準偏低，難以將雨水及時輸送至末端泵站。

（3） 片區現狀5處截污堰部分占用了雨水管渠斷面，降低了雨水管渠的排水能力。以東門渠1號截污堰為例，其占雨水管渠過水面積的18.75%，拆除后堰上水頭最大降低約0.73 m，可顯著改善截污堰上游管渠排水條件。

參考文獻：

［1］ 程曉陶，劉昌軍，李昌志，等.變化環境下洪澇風險演變特征與城市韌性提升策略［J］.水利學報，2022，53（7）：757-768.

［2］ 張利平，杜鴻，夏軍，等.氣候變化下極端水文事件的研究進展［J］.地理科學進展，2011，30（11）：1370-1379.

［3］ 張建云，王銀堂，賀瑞敏，等.中國城市洪澇問題及成因分析［J］.水科學進展，2016，27（4）：485-491.

［4］ 張入財，田金華，陳超輝，等.鄭州“7·20”特大暴雨極端性成因分析［J］.氣象與環境科學，2022，45（2）：52-64.

［5］ 彭瑋怡，薛宇峰，劉春雷.河南鄭州“7·20”極端暴雨成因分析［J］.水利水電快報，2022，43（12）：8-15.

［6］ 國務院災害調查組.河南鄭州“7·20”特大暴雨災害調查報告［R］.北京：國務院災害調查組，2022.

［7］ 王帆，王曉童，莊玲，等.基于鄭州“7·20”特大暴雨移置的山東省傅疃河流域洪水形勢分析［J］.中國防汛抗旱，2023，33（9）：19-24.

［8］ 陳升.上海模擬重演鄭州“7·20”特大暴雨情景分析［J］.中國防汛抗旱，2022，32（9）：56-59.

［9］ 李敏，謝珊，蔣佳鑫，等.武漢市極端暴雨城市內澇風險分析及應對［J］.中國防汛抗旱，2022，32（9）：46-50.

［10］ 張孝奎，馮立超，崔曉莉.鄭州“7.20”暴雨情景下城市韌性內澇防治體系構建研究——以福州濱海新城核心區為例［J］.災害學，2022，37（2）：79-83.

［11］ 張鴻，馬濤，石東凱，等.河南省西平縣小洪河鄭州“7·20”特大暴雨移置洪水響應分析［J］.水利水電快報，2023，44（11）：25-30.

［12］ 孫曉晨，熊麗華，周密，等.基于MIKE URBAN模型的武漢市典型城區內澇研究——I：模型構建與模擬驗證［J］.水資源研究，2021，10（5）：459-469.

［13］ 劉恒恒，齊鵬云，萬能勝，等.巢湖閘下河網洪水演進及閘門調度數值模擬［J］.人民長江，2022，53（5）：41-46.

［14］ 葉永，張芯菲，鄭保敬.基于格子玻爾茲曼法的三維潰壩流數值模擬與驗證［J］.人民長江，2018，49（12）：93-97.

［15］ 吳家陽.發展浸沒邊界-格子Boltzmann方法模擬水電工程流場［D］.武漢：武漢大學，2019.

［16］ 劉峰，胡非.淺水動力學的格子Boltzmann模擬［J］.海洋與湖沼，2004，35（1）：35-40.

［17］ 楊東，侯精明，李東來，等.基于擴散波方法的管網排水過程數值模擬［J］.中國給水排水，2020，36（15）：113-120.

（編輯：江 文）

Simulation analysis of Zhengzhou \"7·20\" extreme rainstorm scenario in Shenzhen City

WU Jiayang 1，LIU Jiaming1，LIU Xiaoxia2

（1.Changjiang Survey，Planning，Design and Research Co.，Ltd.，Wuhan 430010，China； 2.Water Resources Bureau of Shenzhen Municipality，Shenzhen 518035，China）

Abstract：

Shenzhen City is extremely vulnerable to extreme weather events such as typhoon，rainstorm and sea tide.It is necessary to analyze the entire process of the occurrence，development，and disappearance of city waterlogging disasters under extreme weather conditions.Based on urban distributed hydrological-hydrodynamic coupling model and rainstorm transplantation method，the waterlogging disaster caused by the Zhengzhou \"7·20\" extreme rainstorm in Shenzhen Luohu central area was simulated，and the response of waterlogging control system was analyzed.The results indicated that the proportion of area with waterlogging depth greater than 0.3 meter was 13.58%，and the disaster area was highly coincide with those of rainstorm pipes with return period of less than 2 years.The main cause of waterlogging is the insufficient pipe′s conveyance capacity，resulting in the inability to efficiently transport waterlogged water to the pump station′s front pool and remote flooding.This research results can help to improve the understanding and coping ability of extreme rainstorm disasters in Shenzhen City，and provide a reference for the preparation of over-standard flood emergency plans and pre-flood drills.

Key words：

extreme rainstorm； urban waterlogging； distributed hydrological-hydrodynamic coupling model； Zhengzhou \"7·20\" extreme rainstorm； Shenzhen City