張媛媛

（河北省滄州水文水資源勘測(cè)局，河北 滄州 061000）

到21世紀(jì)中葉，我國(guó)城市化水平將達(dá)到60%，城市人口將增加到9.6億[1]。隨著城市規(guī)模不斷擴(kuò)大，城市化水平的提高，城區(qū)不透水面積急劇增加，導(dǎo)致入滲、滯蓄量不斷減少，地表徑流系數(shù)增大，產(chǎn)流時(shí)間變短，徑流量變大，同時(shí)，城市熱島效應(yīng)、阻障效應(yīng)、凝結(jié)核效應(yīng)導(dǎo)致城區(qū)降雨模式和分布發(fā)生了顯著變化，使降雨量城區(qū)大于郊區(qū)。文獻(xiàn)資料顯示北京市增雨系數(shù)（城區(qū)／郊區(qū)）達(dá)到了1.04[2]。當(dāng)前，我國(guó)的城市排水管網(wǎng)建設(shè)與城市化發(fā)展不相適應(yīng)，設(shè)計(jì)標(biāo)準(zhǔn)低，建設(shè)滯后，往往不能滿足排水需要，因而造成了各地城市內(nèi)澇嚴(yán)重的現(xiàn)象。

滄州市位于河北平原的中東部，總面積10456km2，總?cè)丝?17萬(wàn)。市區(qū)面積64.56km2，人口50.3萬(wàn)。按照城市總體規(guī)劃，到2020年市區(qū)人口將達(dá)100萬(wàn)，面積95km2，市域城鎮(zhèn)化水平達(dá)到63%。滄州市屬于資源型缺水地區(qū)，自20世紀(jì)80年代以來(lái)，為了支撐社會(huì)經(jīng)濟(jì)持續(xù)發(fā)展，連年超采地下水尤其是深層地下水，形成了全國(guó)聞名的地下水漏斗，并引發(fā)了地面沉降等地質(zhì)災(zāi)害，2005年滄州市沉降區(qū)中心地面沉降量達(dá)到2236mm[3]，成為我國(guó)地面沉降最嚴(yán)重的城市之一。滄州市區(qū)地勢(shì)平坦，自然排泄能力較差，加之地面沉降影響，市區(qū)內(nèi)積水不得不依靠泵站進(jìn)行外排，排水容量達(dá)到29.8m3／s[4]。較低的排水能力，致使一旦發(fā)生城區(qū)瀝澇災(zāi)害，同等淹沒條件下造成的生命和財(cái)產(chǎn)損失也會(huì)更嚴(yán)重，因此，研究滄州市除澇問題具有較強(qiáng)的現(xiàn)實(shí)意義和深遠(yuǎn)的歷史意義。

1 研究區(qū)域概況

本次研究區(qū)域?yàn)闇嬷菔兄鞒菂^(qū)，主要指渤海路、長(zhǎng)蘆大道、海河路和迎賓大道包圍的范圍，面積為64.56km2。所選區(qū)域位于滄州市轄區(qū)的東偏北位置，地勢(shì)平坦，由西南向東北方向傾斜，坡度約為1／8500，地面海拔高程為7.2～9.5m。該區(qū)降水年際變化較大，最大年降水量1160.7mm（1964年），最小年降水量264.5mm（1968年），多年（1956～2011年）平均降水量563.4mm；降水年內(nèi)分布不均，主要分布在汛期，汛期降水占全年降水的80%左右。

研究資料主要包括降雨量、內(nèi)澇積水量、排水量等基礎(chǔ)數(shù)據(jù)，其中降雨量取用1971～2010年降雨系列，內(nèi)澇積水取用2009～2012年29處監(jiān)測(cè)點(diǎn)的監(jiān)測(cè)數(shù)據(jù)，排水資料取用與內(nèi)澇積水同步的13處泵站的實(shí)測(cè)數(shù)據(jù)。

2 SWMM模型構(gòu)建

2.1 SWMM模型簡(jiǎn)介

SWMM（Storm Water Management Model）是由Metcalf和Eddy有限公司、佛羅里達(dá)大學(xué)和美國(guó)水資源有限公司研制的城市暴雨洪水管理模型，最新版本為SWMM5.0。其主要用于模擬單個(gè)降雨事件或長(zhǎng)時(shí)間尺度（連續(xù)）降雨的徑流水量與水質(zhì)。本研究中主要應(yīng)用SWMM5.0中城區(qū)暴雨徑流、合流制管道、污水管道和其它排水系統(tǒng)的規(guī)劃、分析與設(shè)計(jì)。

2.2 滄州市城區(qū)暴雨內(nèi)澇SWMM模型簡(jiǎn)介

滄州市城區(qū)暴雨內(nèi)澇模型的構(gòu)建，基礎(chǔ)工作包括下墊面信息提取、管網(wǎng)概化、匯水區(qū)劃分和建立模型信息數(shù)據(jù)庫(kù)4部分。

（1）下墊面信息提取。參考研究區(qū)域地形地貌圖和遙感影像圖等地圖資料，利用ArcGIS 9.3在30m×30m的海河流域DEM圖中提取出研究范圍，獲得坡度、坡向數(shù)據(jù)、不透水面積比、曼寧糙率、下滲率等下墊面參數(shù)，存儲(chǔ)于地理信息系統(tǒng)（GIS）軟件中，以方便模型演算過(guò)程中提取。

（2）管網(wǎng)概化。依據(jù)研究區(qū)排水管網(wǎng)13個(gè)泵站、34個(gè)節(jié)點(diǎn)、49條管道、11個(gè)出水口的空間分布及屬性特征，建立管網(wǎng)空間位置、高程、管長(zhǎng)、管徑、流向、坡度等數(shù)據(jù)庫(kù)，并根據(jù)排水標(biāo)準(zhǔn)確定，排水區(qū)劃分，排水口選定，地面積水深度、時(shí)間及變化過(guò)程，雨水口的位置、數(shù)量，管網(wǎng)位置、管長(zhǎng)、管徑、管材、流速、流量，雨水排水泵站位置、數(shù)量、裝機(jī)容量等要素，對(duì)管網(wǎng)進(jìn)行概化。

（3）匯水區(qū)劃分。如圖1所示，依據(jù)研究區(qū)域下墊面特征，泵站分布、管網(wǎng)概化成果等資料，采用信息相似性、地理位置臨近性確定匯水區(qū)排水邊界，參考下水管道入口位置，劃分為81個(gè)集水區(qū)。

圖1 管網(wǎng)概化后研究區(qū)域的節(jié)點(diǎn)與管網(wǎng)圖

2.3 模型參數(shù)率定

SWMM模型參數(shù)較多，匯水區(qū)參數(shù)包括坡度、不滲透面積、不滲透面積粗糙系數(shù)、滲透面積粗糙系數(shù)及不滲透面積洼地蓄水量、滲透面積洼地蓄水量以及Horton公式中的最小下滲率、最大下滲率、下滲衰減系數(shù)等，節(jié)點(diǎn)／管段參數(shù)主要包括節(jié)點(diǎn)內(nèi)底、節(jié)點(diǎn)最大深度、節(jié)點(diǎn)積水面積、管渠長(zhǎng)度、管渠幾何尺寸、管渠粗糙系數(shù)、流量單位、管段偏移等。對(duì)于不確定的參數(shù)，研究過(guò)程中根據(jù)下墊面特征，按經(jīng)驗(yàn)參數(shù)進(jìn)行初始值預(yù)估，通過(guò)優(yōu)選獲得參數(shù)最優(yōu)值。

參數(shù)率定過(guò)程中，根據(jù)2009年8月16日、2011年7月1日以及2012年8月1日實(shí)測(cè)降雨量及實(shí)測(cè)積水點(diǎn)積水深度，以子流域漫流寬度、不透水面積比例等產(chǎn)流參數(shù)和糙率、管道埋深等匯流參數(shù)等主要參數(shù)調(diào)整為基礎(chǔ)，以積水點(diǎn)的積水深度效率函數(shù)Ens為模型參數(shù)最優(yōu)化調(diào)試的目標(biāo)函數(shù)。

積水點(diǎn)的積水深度效率函數(shù)Ens為

式中，h0為積水深度觀測(cè)值，cm；hs為積水深度模擬值，cm；hav為積水深度觀測(cè)值的平均值，cm。

經(jīng)過(guò)反復(fù)試算調(diào)整，最終確定內(nèi)澇模型的主要參數(shù)，見表1。

表1 內(nèi)澇模型主要參數(shù)率定結(jié)果

3 設(shè)計(jì)暴雨研究

暴雨強(qiáng)度是指一定頻率的暴雨在單位時(shí)間內(nèi)的平均降雨深度，是描述暴雨特征的重要指標(biāo)。城市短歷時(shí)暴雨強(qiáng)度公式是設(shè)計(jì)城市雨水排水設(shè)施的主要依據(jù)，目前普遍采用推理公式，主要是通過(guò)城市暴雨強(qiáng)度公式來(lái)確定雨水設(shè)計(jì)流量[5]。通過(guò)對(duì)暴雨資料分析計(jì)算，確定滄州市城區(qū)暴雨強(qiáng)度公式，借助典型暴雨，利用同頻率放大法推求設(shè)計(jì)暴雨過(guò)程。

3.1 暴雨資料分析

（1）基礎(chǔ)資料選擇。研究顯示，資料系列越長(zhǎng)，暴雨強(qiáng)度公式精度越高[6]。本次研究采用滄州市南運(yùn)河捷地站、南運(yùn)河滄州站、滄浪區(qū)南支小元站1971年～2010年40年暴雨資料，每場(chǎng)降雨選擇5min，10min，15min，20min，30min，45min，60min，90min，120min，180min，360min，720min，1440min共9個(gè)歷時(shí)進(jìn)行雨量統(tǒng)計(jì)[7]。

（2）暴雨資料選樣。暴雨選樣方法主要包括年最大值法和非年最大值法，非年最大值法又包括年超大值法、超定量法和年多個(gè)樣法，其中年最大值法主要被水利水電部門和氣象部門采用，而年多個(gè)樣法主要被市政排水部門所采用。本研究中，采用年多個(gè)樣法，即選擇每年各個(gè)歷時(shí)中的6個(gè)最大值。進(jìn)行降雨量降序排列，從上到下選擇資料年限4倍的數(shù)據(jù)，作為編制暴雨強(qiáng)度公式的基礎(chǔ)資料，因此本研究中每個(gè)歷時(shí)的樣本數(shù)為40×4＝160個(gè)。

（3）頻率計(jì)算。進(jìn)行樣本資料經(jīng)驗(yàn)頻率計(jì)算，采用P-III型分布曲線進(jìn)行理論頻率曲線擬合，再由暴雨強(qiáng)度頻率曲線反查不同重現(xiàn)期的設(shè)計(jì)暴雨強(qiáng)度。

3.2 暴雨強(qiáng)度公式參數(shù)求解

暴雨強(qiáng)度公式分為單一重現(xiàn)期公式和統(tǒng)一公式。

單一重現(xiàn)期公式

統(tǒng)一公式

式中：A1，C，b，n分別為待求參數(shù)，無(wú)量綱；i為暴雨強(qiáng)度，mm／min；t為降雨歷時(shí)，min；P 為重現(xiàn)期，a。

傳統(tǒng)暴雨強(qiáng)度公式參數(shù)求解較為繁瑣，研究中引入R軟件和MATLAB軟件進(jìn)行參數(shù)求解，大大減少計(jì)算量的同時(shí)，避免了計(jì)算錯(cuò)誤的發(fā)生。首先以重現(xiàn)期1年為例，確定單一重現(xiàn)期公式的值A(chǔ)，b，n。由 R 軟件計(jì)算參數(shù)的迭代初值A(chǔ)1（0），C（0），b（0），n（0），再通過(guò) MATLAB軟件最終確定各個(gè)參數(shù)。最終確定暴雨強(qiáng)度公式：

3.3 設(shè)計(jì)暴雨過(guò)程

滄州市區(qū)設(shè)計(jì)暴雨過(guò)程的推求可以不考慮雨量在空間分布的不均勻性，忽略點(diǎn)雨量與市區(qū)面平均雨量的差別。雖然一場(chǎng)實(shí)際暴雨的空間分布是不均勻的，雨量在空間分布的梯度可能相當(dāng)大，但是由于暴雨中心的位置是不確定的，使點(diǎn)雨量與面雨量之間的差別不那么明顯，所以采用市區(qū)附近雨量站設(shè)計(jì)雨量代替市區(qū)面平均設(shè)計(jì)雨量。研究中以2005年7月23日10時(shí)～24日10時(shí)典型暴雨，利用同頻率放大法推求設(shè)計(jì)暴雨過(guò)程。

計(jì)算結(jié)果顯示，P＝0.33a，P＝1a，P＝3a，P＝5a，P＝10a，P＝20a的各種重現(xiàn)期下，歷時(shí)24h的設(shè)計(jì)暴雨量分別為28.8mm，72mm，115.2mm，133.5mm，158.4mm，187.2mm，設(shè)計(jì)暴雨過(guò)程見圖2。

圖2 不同重現(xiàn)期設(shè)計(jì)暴雨過(guò)程

4 內(nèi)澇模擬

滄州暴雨內(nèi)澇模型，可以根據(jù)不同情景的降雨參數(shù)，對(duì)積水開始時(shí)間、積水深度、積水歷時(shí)、積水面積、峰值流量、總流量等進(jìn)行模擬。由重現(xiàn)期P＝0.33a，P＝1a，P＝3a，P＝5a的設(shè)計(jì)暴雨過(guò)程，模擬相應(yīng)的內(nèi)澇情況。以P＝1a為例，模擬結(jié)果見表2。將模擬成果在ArcGIS 9.3中進(jìn)行可視化顯示輸出，見圖3。

表2 1年一遇管段節(jié)點(diǎn)內(nèi)澇結(jié)果

5 結(jié)論及建議

圖3 1年一遇內(nèi)澇積水空間分布

對(duì)滄州市城區(qū)內(nèi)澇的研究，從城市排水管網(wǎng)系統(tǒng)模擬研究現(xiàn)狀出發(fā)，以排水管網(wǎng)的降雨模擬為研究對(duì)象，基于地理信息系統(tǒng)的數(shù)據(jù)處理和空間分析等功能，利用SWMM排水管網(wǎng)模擬系統(tǒng)的理論與方法，建立了暴雨內(nèi)澇模型，使城市內(nèi)澇模擬更加簡(jiǎn)化和直觀。

滄州市暴雨內(nèi)澇模型研究以雨量、積水、排水等為基礎(chǔ)資料，相比雨量資料，積水、排水等基礎(chǔ)資料主要來(lái)源于2009年滄州城市水文開展以來(lái)的實(shí)測(cè)數(shù)據(jù)，資料系列較短，對(duì)模型參數(shù)率定及模型模擬成果精度有一定影響，應(yīng)隨著積水、排水資料系列的延長(zhǎng)，不斷修訂模型參數(shù)，提高模型的模擬精度。

[1]錢正英，張興斗.中國(guó)可持續(xù)發(fā)展水資源戰(zhàn)略研究綜合報(bào)告——中國(guó)工程院“21世紀(jì)中國(guó)可持續(xù)發(fā)展水資源戰(zhàn)略研究”項(xiàng)目組[J].中國(guó)工程科學(xué)，2000（8）：1-17.

[2]周翠寧，任樹梅，楊培嶺，等.城市化對(duì)降雨特征影響研究[J].水利水電技術(shù)，2007，38（10）：62-65.

[3]劉振鐸，王春橋，吳興國(guó).滄州市地下水超采區(qū)核定及禁采區(qū)、限采區(qū)劃分[J].水利科技與經(jīng)濟(jì)，2011（10）：65-68.

[4]劉振鐸，朱艷飛，哈建強(qiáng).滄州積水成因及對(duì)策簡(jiǎn)析[J].甘肅水利水電技術(shù)，2011（06）：37-39.

[5]王寶泉.淺談我國(guó)城市雨水排水系統(tǒng)中的問題與對(duì)策[J].河北工程技術(shù)高等專科學(xué)校學(xué)報(bào)，2008（4）：16-19.

[6]邱兆富，曾曉嵐，郝以瓊，等.重慶主城暴雨強(qiáng)度公式推求[J].重慶建筑大學(xué)學(xué)報(bào)，1999（6）：16-19.

[7]上海市政工程設(shè)計(jì)研究總院.GB 50014-2006室外排水設(shè)計(jì)規(guī)范[S].北京：中國(guó)計(jì)劃出版社，2006.