基于“9·7”特大暴雨的深圳河流域防洪能力分析

關鍵詞：特大暴雨；暴雨洪水；防洪能力；深圳河流域

中圖分類號：TV87 文獻標識碼：A 文章編號：1001-9235（2024）11-0034-09

洪水是全球危害最大的自然災害之一，具有發生頻率高、破壞力強及影響范圍大等特點［1-2］。據統計，2000—2018年，全球共有2. 7億人遭受到了洪水的影響，經濟損失高達6 400億美元，共造成約10萬人死亡（全球緊急災害數據庫）［3］。作為全球受洪水影響最大的國家之一，中國約有三分之二的國土面積、人口和工農業生產值均遭受洪水威脅［4］，2000—2018年的洪災共造成22293人死亡和33562. 22億元的直接經濟損失［5］。同時，隨著全球氣候變暖，海平面持續上升，極端暴雨及洪水的頻率和強度均顯著增強［6-7］，且隨著城市化的迅速發展，不少城市防洪工程的防洪能力降低或已不能滿足現有需求，洪災問題日益嚴重［8-9］。如2021年鄭州發生的“7·20”特大暴雨水災，導致150個縣、1 478 萬人受災［10］。因此，對現有防洪工程的防洪能力進行再評估顯得尤為必要。

深圳河作為深圳和香港的界河，是深圳市最重要的排洪入海河流，具有重要的區位優勢。1985年以來，經過1—4期的河道治理工程建設，深圳河流域的防洪標準已提升至50年一遇［11］。但在全球氣候變化及人類活動的影響下，極端暴雨天氣頻發（如2018年“8·29”暴雨和2023年“9·7”暴雨），深圳河下游潮位極值上抬，城市化導致土地硬化加速了暴雨的產匯流速度，增加了流域內的洪水風險。同時，河道清淤后發生快速回淤，導致河道淤積并高于泄洪設計底高程，減小了深圳河的泄洪能力［12］。2023年9月7—8日，“海葵”臺風殘余環流、季風和弱冷空氣的共同作用下，深圳河流域發生極端特大暴雨，強度超強、持續時間超長、范圍超廣，打破了1952年有氣象記錄以來的最大2、3、6、12、24、48、72 h共7項歷史雨量極值，導致深圳河干流防洪壓力劇增［13-14］。因此，亟需對深圳河流域現有防洪能力進行分析。本文通過對深圳河流域（圖1）2023年“9·7”特大暴雨下的洪水模擬，對該流域的防洪能力進行復核分析，發掘深圳河流域防洪所面臨的問題與挑戰，從流域層面提出一套防洪能力提升的總體方案，以期為深圳河流域的防洪治理提供切實可行的參考方案。

1研究區域及水情概況

1. 1研究區域概況

深圳河流域（22°27'～22°39'N，114°～ 114°13'E）位于珠江口的東側，發源于梧桐山牛尾嶺，自東北向西南流入深圳灣，全長37 km，共有6條支流，其中干流為沙灣河與蓮塘河匯合后的三岔河口至河口處，長18 km［15］，整個流域面積為312. 5 km²。深圳河流域多年平均降雨量為2 023 mm，主要集中在4—9月，約占全年降雨量的85%。深圳河干流主要以潮流動力為主，大潮潮差在2. 5 m以上，徑流動力較弱，河口處多年平均徑流量為4. 73億m³。深圳河水系呈扇形分布，三岔河口以上的上游區域位于丘陵山地，坡度較陡，河床比降在2‰～4‰，中下游區域位于平坦的沖積平原，坡度較緩，河床比降僅在0. 10‰～0. 17‰。因此，支流河流短、比降陡，暴雨后約1 h洪峰即可到達干流，位于流域中心位置的鹿丹村干流河段的洪水風險增大。同時，由于城市化進程加劇引起的流域降雨產匯流速度進一步加快、氣候變化導致的極端災害頻發導致的河道洪潮歷史水位不斷新高，深圳河流域面臨著越來越顯著的洪災風險，對深圳及香港區域內的人民生命安全和社會經濟造成巨大威脅［16］。

1. 2水情概況

深圳河流域“ 9·7”特大暴雨發生于2023年9月7日17時至9月8日15時30分，歷時近23 h，主要集中在9月7日17時至9月8日7時，共14 h。流域降水總量為1. 4億m³，全流域雨量均較大，最大24 h降雨暴雨中心位于小梧桐一帶、布吉河筍崗滯洪區至香港梧桐河一帶，1、6、24 h最大點雨強均超100年一遇洪水標準（表1）。

由于深圳河流域內城市化程度高，不透水表面面積占比大，徑流系數大，且“蘇拉”和“海葵”雙臺風間隔短，土壤含水飽和度高，同時降雨中心自流域上游向下游移動，洪峰與雨峰疊加，多重因素的共同作用導致本次暴雨具有產匯流時間短、產流量大、洪峰出現時間提前等特征。其中雨峰和洪峰間隔最近的為正坑水大望站附近，僅25min［ 13-14］。

“9·7”暴雨雖然發生在小潮的高潮位，受到的頂托作用較弱，但深圳河干流洪水位均超歷史最大值，河口以上歷史最高水位提升了0. 64～1. 43 m（表2）。深圳河的主要河流全部超警戒水位，其中鹿丹村斷面由于主要支流集中匯入，最大流量達到了1 840 m³/s，超過了 50年一遇設計洪峰流量，為有實測資料以來的最大洪水。暴雨期間，深圳河沿程水位未超過防洪墻，但皇崗至羅湖口岸部分河段水位超出堤面，平均超堤面高程約0. 35m，深圳側流域內水深超過0. 27m的淹沒面積達13. 1 km²，超過1m的淹沒面積達3. 7km²。

2結果與討論

2. 1洪水模擬分析

為了對“9·7”暴雨期間深圳河流域的防洪能力進行分析，本研究基于Mike模型構建深圳河流域一維模型，模型采用2022年最新地形，上邊界包含上游支流蓮塘河（四期段）、梧桐河、布吉河和水庫排洪河，下游邊界為深圳灣赤灣斷面（圖2）。模擬時段為2023年9月7日18時至8日15時，共22 h。將流域劃分圖3a所示的匯水區，根據閘泵出流、單位線統計和徑流系數等方法計算整個流域內的產匯流流量，其中閘泵出流為實測水庫、閘泵出流流量數據，單位線統計為Guo等［17-19］提出的廣義單位線匯流計算方法，計算支流流域面積大于10 km²的區域，徑流系數為綜合徑流系數法［20］，計算支流流域面積小于10 km²的區域。計算得到該時段內進入深圳河道的總水量達9 850 萬m³，平均徑流系數為0. 69，各支流的流量過程見圖3b。

模擬前根據清淤策略研究，將深圳河的糙率設置為0. 025。由于蓮塘河香港側存在大量樹木，導致其糙率要高于設計糙率。因此根據模型率定及驗證后，將該河段糙率定為0. 05，蓮塘河的模擬水位與實測水位基本吻合。沿程模擬水位結果見圖4，其中橫坐標為距離赤灣斷面距離，模擬水位與實測水位的平均誤差在0. 10m以內，說明模擬結果可信。深圳河流域沿程各斷面逐時水位變化的模擬值與實測值見圖5。

模擬結果表明，“9·7”暴雨期間深圳河流量量級達到了50年一遇等級，其中蓮塘河的蓮麻坑以上流量接近200年一遇洪水等級，部分河段發生了漫堤，模擬水位比四期工程設計水面線高程平均升高0. 80 m，最高超高處為蓮麻坑河段，超過1. 40m（圖6）。深圳河干流水位未達到50年一遇等級，其中皇崗河口以下、布吉河口以上河段水位未達20年一遇等級，皇崗河口至布吉河口河段的平均水位比50年一遇低0. 32 m，比20年一遇高0. 18 m（圖7）。

2. 2防洪能力分析

“9·7”暴雨期間，通過工程調度方法最大限度地發揮了深圳水庫的調蓄滯洪錯峰作用，以及筍崗滯洪區對布吉河洪峰流量的削減作用，最大程度地保證了深圳河流域的安全［12-13］。但此次洪水仍有暴露出防洪問題，如全線超警戒水位、部分河段出現漫堤，因此有必要對現深圳河防洪現狀能力進行復核分析。深圳河流域的洪澇風險主要是由強降雨導致，但河道自身的地形變化及糙率變化在洪澇災害的防御中不可忽視，且這兩個因素更易受到人為影響。因此，本研究主要從深圳河的地形及糙率變化對防洪能力進行分析，以期對后續深圳河流域的防洪整治提供借鑒。

2. 2. 1地形變化對防洪能力的影響

2016—2022年，深圳河流域進行了兩次大規模的清淤工程。2017年6月至2017年11月，實施了深圳河2016年度清淤工程，對深港兩側共計41個涵口及重點河段進行了清淤，累計清淤量約15. 68萬m3。2018年7—12月，實施了深圳河梧桐河口至上步立交河段清淤工程，累計清淤量約49. 76 萬m3。通過對比深圳河2016年和2022年的實測地形，發現相較于2016年，2022年深圳河淤積量減小了32. 7萬m³，且以平原河口為界，河道下游以沖刷為主，上游以淤積為主，地形平均加深了0. 22 m。

基于2016、2022年的地形，分別模擬了不同地形下深圳河在50、200年一遇洪水情景下的沿程水位變化，見圖8。結果表明，在兩場洪水情景下，2022年地形下的深圳河在皇崗河口至平原河口的水位均低于2016年地形下的水位，其中福田河口至布吉河口降低了0. 04～0. 13 m，布吉河口至羅湖口岸降低了0. 11～0. 25m，羅湖口岸至沙灣河口降低了0. 22～0. 25 m。但在深圳河口至皇崗河口、平原河口至蓮塘水文站則呈現相反的表現，即2022年地形下的深圳河水位高于2016年地形下的水位。其中兩場洪水情景下2022年地形下的蓮塘河口水位均比2016年地形下的水位高了0. 35 m左右，在50年一遇洪水情景下，2022年地形下的蓮塘水文站水位比2016年地形下的水位高了0. 11m，而在200年一遇洪水情景下高了0. 21 m。造成該結果的原因考慮為主要與河道地形相關，在下游沖刷處水位以降低為主，在上游淤積處水位以抬高為主。

根據《海堤工程設計規范》計算，深圳河流域防洪標準需考慮防洪墻0. 4 m 超高。通過對深圳河2016年和2022年不同地形下的水位進行模擬后發現，基于2016年的地形，福田河口至沙灣河口總體上達不到50年一遇洪水標準，其他河段均能滿足50年一遇的防洪標準，其中深圳河口至皇崗河口、四期蓄滯洪區的防洪能力達到了200年一遇洪水標準（圖9）。而基于2022年的地形模擬發現，深圳河干流的整體防洪能力基本都能達到50年一遇洪水標準，但除了深圳河口至皇崗河口外，上游河段均未達到200年一遇洪水標準（圖10）。

2. 2. 2糙率變化對防洪能力的影響

除了河床地形變化對防洪有影響外，河道糙率變化的作用也極為重要。本研究基于2022年的地形，對蓮塘河河段設計糙率（0. 025）和推算糙率（0. 05）下50年一遇洪水及200年一遇洪水情景下的沿程水位變化進行了模擬（圖11）。結果表明，糙率的增加會促使蓮塘河上游水位壅高，壅高值在0. 4～0. 9m，其中在50年一遇洪水情景下，蓮塘口岸和蓮塘水文站的水位分別抬高了0. 48、0. 77m，在200年一遇洪水情景下，則分別抬高了0. 56、0. 87m。糙率的增加減緩了洪水的流速，增大了洪水的風險。在考慮防洪墻0. 4m超高的條件下，蓮塘河上游河段已不能滿足50年一遇的防洪標準，蓮塘口岸至蓮塘水文站河段的平均欠高為0. 24m，最大欠高達到了0. 58m。

3防洪挑戰及治理建議

3. 1面臨的防洪挑戰

3. 1. 1不利的下墊面條件

a）扇形集雨地貌疊加城市化發展，導致暴雨產匯速度加快。從流域地貌來看，深圳河流域呈扇形，三面環山，支流坡降大，導致暴雨匯流速度快。同時疊加城市化進程加速的不利影響，下墊面硬化，暴雨匯水速度快且集中于深圳河中游。深圳河上游主要支流梧桐河、布吉河、沙灣河和蓮塘河的集雨面積共占流域面積的75%，各支流的徑流主要集中在布吉河口至梧桐河口約1 km的河段內，再通過下游約9 km的河道進行排洪。由于深圳河干流的河床比降小，且受海潮頂托作用及各支流的相互頂托作用，布吉河口河段極易形成高水位，成為防洪敏感區。此外，近幾十年來，深圳市區城市化發展迅速，流域內的大量農田、植被等自然土地已被修建為房屋、混凝土/瀝青道路等不透水面，加劇了暴雨的產匯流［13］。“9·7”暴雨期間，深圳河流域布吉站的雨峰與布吉河草埔斷面的洪峰只相隔1 h 左右，產匯流時間短。在流域中下游的老城區內發生嚴重內澇，加劇受災風險。

b）河道清淤后回淤明顯，現有過洪斷面面積遠小于設計過洪斷面。深圳河河道長期呈淤積趨勢，雖經過幾次大規模或局部清淤，但由于河道干流徑流作用小，潮流作用明顯，且深圳河口存在的大范圍淺灘為深圳河提供了源源不斷的淤積泥沙來源，導致河道清淤后泥沙回淤顯著。通過對深圳河一二期河段多年淤積量分析發現，在2006—2007年大規模清淤和2017—2018年局部清淤工程結束后，深圳河河道均發生了快速回淤，并在長期保持緩慢淤積態勢。深圳河干流多年平均徑流量僅為15 m³/s，泥沙的回淤使得河道無法維持安全泄洪斷面，無法保障近2 500 m³/s流量的安全宣泄。現有斷面面積不足設計斷面的1/2，計算得到的泥沙沖淤平衡斷面面積僅能維持原設計斷面的1/3，行洪能力與原設計目標差距大（表3）。

c）河口紅樹林不斷擴張，嚴重擠占河口泄洪斷面。近幾十年來，深圳河口的紅樹林不斷向下游擴張。2000年，香港側紅樹林前沿位于分汊口以上900 m，目前已擴張至下游600 m。同時淺灘紅樹林的生長擠壓河道，導致深圳河河口內平均河寬為200 m，而最窄處僅150 m，嚴重擠占河口泄洪斷面。

3. 1. 2極端災害天氣頻發

極端災害天氣頻發導致深圳河流域面臨更加嚴峻的防洪挑戰，主要體現在超強臺風、暴雨的頻次及強度顯著增加。2017年以來，連續發生天鴿、山竹臺風以及“8·29”和“9·7”暴雨。極端天氣迫使深圳河水位發生抬升，其中2場臺風使得赤灣潮位抬升0. 48 m；“8·29”洪水期間深圳河水位超過了治理工程設計水位；相較于“8·29”洪水，“9·7”洪水導致深圳河水位又抬升0. 99～1. 46m。此外，海平面的持續上升將導致設計潮位的抬升。根據香港《暴雨排放系統手冊》所示，2018年尖鼻咀200年一遇設計潮位較1995年抬升了0. 82 m，并預測在21世紀中期將抬升0. 23m，在21世紀末將抬升0. 49m。

3. 1. 3內澇風險嚴重

河道高水位的頂托作用是導致城市發生內澇的一個重要原因。由于周邊區域地勢低洼，深圳河干、支流的高水位頂托作用將導致城市內部難以排水。“9·7”暴雨期間，洪水將深圳河河口上游的最高水位提高了0. 99～1. 42 m。雖然深圳河干流未發生漫頂現象，但高水位頂托作用導致深圳河支流排洪和流域內的陸域排水不暢，導致兩岸發生了嚴重內澇。此外，本次暴雨期間僅遭遇小潮，若遭遇大潮高潮位、更大量級暴雨或海平面升高變化，深圳市將面臨更強的洪水頂托作用，城市內澇排水問題更加嚴峻。

3. 2治理思路

在以往對深圳河防洪治理過程中，實施了不少對流域防洪能力提升的措施，包括工程和非工程措施。目前已有的工程措施包括加高堤防和清淤清障。但加高堤防會導致束水歸槽、抬高水位、加重內澇、影響景觀等問題，而清淤清障則面臨疏浚量大和快速回淤的問題。已有的非工程措施包括深圳水庫挖潛，但目前起調水位已降低至26. 3m附近，調蓄庫容仍有限［21-22］。因此，在海平面上升和極端天氣頻發的氣候變化背景下，需要改變以往加高堤防、擋洪為主的治理對策，從流域層面對防洪能力的提升進行總體方案布局，可以采取蓄、滯、分、泄、擋等各項措施相結合的方法，進一步削減洪峰流量和增加河道行洪能力［23-25］。

4結論

本文通過對深圳河流域2023年“9·7”特大暴雨洪水的模擬，復核分析了深圳河流域的防洪能力，并探討了深圳河流域防洪未來可能面臨的挑戰，提出了一套綜合治理方案建議，主要結論如下。

a）2016年深圳河福田河口至沙灣河口河段總體上達不到50年一遇洪水標準，但經過清淤工程的開展，2022年深圳河干流的整體防洪能力基本都能達到50年一遇洪水標準。

b）由于河床淤積抬高和糙率的增加，深圳河支流蓮塘河上游河段會發生雍高，已不能滿足50年一遇的防洪標準，蓮塘口岸至蓮塘水文站河段的平均欠高為0. 24m，最大欠高達到了0. 58 m。

c）深圳河的防洪面臨著匯流快、清淤后回淤明顯、紅樹林擠壓泄洪斷面等不利下墊面條件，以及極端災害天氣頻發、內澇風險嚴重等挑戰。在海平面上升和極端天氣頻發的氣候變化背景下，需要改變以往加高堤防、擋洪為主的治理對策，從流域層面對防洪能力的提升進行總體方案布局，采取蓄、滯、分、泄、擋等各項措施相結合的方法，進一步削減洪峰流量和增加河道行洪能力。