城市暴雨内涝数学模型地地地研究与的应用

X晓

〔某某工业大学，某某，某某，120330270〕

摘要：暴雨内涝对城市的影响日益严重，为了城市能够更好的应对暴雨带来的冲击，本文以城市的街道路面与河道水流的运动为对象进展模拟，建立了某某市暴雨内涝积水数学模型。模型以平面二维非恒定流根本方程和不规如此网格划分技术为框架，采用简化分类处理的方法，将通道分为路面型、河道型以与特殊通道型，根据不同类型简化动量方程，求任一网格各个通道上的单宽流量。根据不规如此网格的方法，按照某某市的地形进展多边形计算网格的设计。介绍了数学模型在某某市的应用和误差分析以与城市路面降雨量的计算。 关键词：城市暴雨 内涝灾害 数学模型 误差分析 某某市

Research and Application ofWuhan Waterlogging Mathematical Model

LiuXiao

(,Hubei University of Technology, Hubei,Wuhan,120330270)

Abstract:Waterlogging increasingly serious impact on the city, in order to respond to storm the city the impact of urban road surface better and the main river flow motion simulation object, the mathematical model of urban storm water waterlogging.The basic equation model for unsteady flow and irregular unstructured meshing technology as the backbone, the use of simplified classification method,the channel into the river type, road type,special channel type, depending on the type of simplified momentum equation,seeking grid unit discharge any individual channel.According unstructured irregular grid design ideas, according to the terrain features are designed in Wuhan polygon putational grid.Describes analysis methods and mathematical models to calculate surface rainfall in the city of Wuhan and application errors.

Keywords: urban storm; waterlogging disasters; mathematical model;model erroranalysis;Wuhan 1 引言

城市内涝是由于强降雨超过城市排水能力而产生的城市内积水的灾害。降雨强度大是造

1 / 8

word 成城市内涝的原因客观原因，X围集中。降雨强度较大、时间较长有可能形成积水，降雨急的地方也可能形成积水。近年来，暴雨对城市内涝灾害的影响日益严峻，不仅加重了城市的防汛排涝的任务，同时对经济建设和城市居民的日常生活也造成了较大影响。因此，为了进一步提高城市排水的服务能力和水平，对城市内涝进展研究已然十分必要，从而减少暴雨对城市的影响。

巴黎的下水道始建于1854年，用时24年巴黎修建了600多公里的下水道。截止今日，巴黎的下水道总长度近2300公里，相当于巴黎到伊斯坦布尔的距离。巴黎的下水道位于地面以下50米。水道纵横交织，密如蛛网，下水道十分宽敞，中间是宽约3米的排水道，两侧是宽约1米的供检修人员通行的便道。巴黎的下水道按沟道大小可以分为小下水道、中下水道和排水渠三种。每天有120亿立方米的水经此管道排出，这些下水道犹如这座大城市的消化系统。因此，暴雨不会使巴黎的形成内涝。

某某是目前被认为是管道设置最合理的城市，某某率先采用雨水污水分流的排水模式。德式排污管道设计巧妙，管道为倒置的鸭蛋型，上大下小，下半局部呈“V〞型，在流量一样的情况下，污水流速更快，即使污水流量比拟小，污水依然能保持一定的速度，使得污水中所夹杂的污物被冲击的没有沉淀的机会，有效保证了管道的畅通。

某某湖泊的减少，大大减小了城市的调蓄能力，从而大大增加了城市内涝的发生率。1950年代，某某市湖泊面积达到1581平方公里，城中湖有127个，80年代已经缩减到874平方公里。近30年来，又减少了228.9平方公里，城中湖也缩减为38个。与之相对应的却是，城市建成区总面积从220平方公里增加到475平方公里。天上的雨水没有减少，调蓄的湖泊在减少，城市建成区的面积却在增加，如果抽排能力没有加强，城市内涝就不可防止。

以城市的河道水流与街道路面的运动为对象进展模拟，建立了城市暴雨内涝积水的数学模型。模型以平面二维非恒定流根本方程和不规如此网格划分技术为主要方法。一般的数值模式计算，根本上采用规如此矩形网格或三角形网格对计算区域进展划分，但是城市的局部地区往往地形复杂，河流、街道分布也各不一样，对于建筑物，尤其是防洪排涝工程在南、北方有着巨大的差异，对径流的形成影响很大。因此，根据城市暴雨径流的特征，模型选取无结构不规如此网格。同时，根据城市各异的地貌特征可将网格设计为三边形、四边形或五边形，网格的各边如此定义为通道。选择网格单元为研究体，取网格中心位置得水位为计算

2 / 8

word 水位H，取网格周边的通道中点处得流量为计算流量Q。水位代表网格的平均值，流量也表示通道的平均值。此外，在时间上，流量和水位采取时间交织的计算方式。根据某某市的地形设计计算网格。某某市地形平缓，海河较宽河道设为河道型网格，较窄的明沟和河道为特殊通道。

2 根本方程

由高斯定理可将连续方程转化为

HT2dti2dtKHQikLik2qiTdtdt (1) Aik1Ti式中Hi——i单元的水深；

Ai——i单元的网格面积；

Qik——i单元的第k条通道上的流量；

Lik——i单元的第k条通道长度；

qi——i单元的源； dt——T时刻时间的增量。

根据简化的动量方程的离散格式来计算Qik。对于不规如此的网格，如果按动量方程严格计算任一方向通道的流量，如此计算量太大，所以该模型中采取简化分类处理的方法。假设同一时段同一网格水位的变化不大，采取分类处理的方法去求任一网格各通道上的单宽流量。通道分为路面型、河道型、特殊通道型、其他型，由不同类型的简化动量方程，计算流量。

河道型通道对应的动量方程：

2TdtTdtTTnQQZikZikikikik12gHi2dt2gdt (2) 7dLikHi3QTdtikQTdtik式中 g——重力加速度；

Zik1、Zik2——为i单元的第k条通道左右侧1、2断面的水位； dLik——i单元的第k条通道左右侧单元形心距离； nik——i单元的第k条通道上的糙率。

地面型通道对应的动量方程：

3 / 8

word QT2dtik2N2dtNTdtTdtTdt (3) HQLQL/22qikikdijdijdiAdik1j1Tdi对于有缺口堤或连续堤的通道，其流量采取宽顶堰流的计算公式，即

dtQTmikбik2gHi3/2 (4) ik式中Σik——i单元的第k条通道下游淹没系数； mik——i单元的第k条通道流量系数。

某某市较窄的河道认为是特殊通道，它与其左右侧网格之间的流量，选取宽顶堰流公式(4)进展计算；对于河道的单宽流量Qs选择与河道型通道的动量方程式(2)相近的公式来计算；此外，特殊单元水深如此按照以下的离散格式计算：

HT2dtdi2N2dtNTdtTdtTdtQikLikQdHL/22qddidt (5) ijijAdik1j1Tdi式中Hdi——i特殊单元的水深； Adi——i特殊单元的面积； Qdij——i通道周边的单宽流量；

Ldij——i通道周边的边长；

qdi——i特殊单元上的降雨强度，一般选取其特殊节点处的降雨强度来作为这个单元的平均值。

3 某某市的排水系统

某某市的排水系统由水井、管道、泵站、明沟、闸门等组成，排水管道错综复杂，通过雨水井来承当和接收来自地面的水，通过泵、闸或淹没出流管道来实现管道水的排出。对于地下排水管网内的水流，一般选用一维非恒定流方程进展计算。地面单元与地下排水管网的水量交换，一般可以通过连续方程源汇项来进展控制，即可将源汇项分成单元排水强度和有效降雨强度两局部。单元的排水量可依据单元的最大排水能力、单元积水深度，以与降雨强度来确定：当单元最大排水能力大于降雨强度并且单元积水相对较浅时，降雨量将全部转化成排水量流入到地下管网，故单元排水强度如此为降雨强度；反之，当单元最大排水能力小于降雨强度且单元的积水相对较深时，单元排水强度可达到最大的排水强度，单元最大排水强度与降雨强度之差如此转化成地表径流。

对于计算管道内的流量选用一维非恒定流方程的离散形式：

4 / 8

word dtdtQTQTpipi2TdtTdtnQQZiT2ZiTpipi12gApidt2gdt 〔6〕 4/3dLiApiRpi式中Qpi——计算管道的断面流量；

dLi——相临的两网格的形心到通道中点的距离之和； Rpi——水力半径；

计算管道内的水深一般选用连续方程的离散形式:

HT2dtpi2dtNTdtTdtAH Q2qdt (7) 1,pk ,ApiAk1piTpi式中Hpi——管道水深； q1——管道排水强度；

N——单元中有排水管道通过的通道数；

A,pi——计算单元内管道的等效底面积。

降雨的过程当中，地面积水会向管道聚集，然后管道沿着管网系统集合直至出口，再在出口处经泵站或闸门向河道泄水。抽排水的管道内水量控制着泵站的排水量，当管道内水量小于排水能力时，泵的排水量可以自动调节。闸门的排水量同样按自流口门的设计流量和管道内水量为条件进展调节。河道与管道的水位差调节着淹没出流管道的排水量，当河道水位低于管道水位时，管道水将流向河道。反之，将会出现倒流。

4 某某市面雨量计算

一般必须采取有效的插值法来获得降雨量，进而转换为无结构且不规如此网格的面雨量，此处选用二次曲面法来计算面雨量。

假设有N个雨量观测站。对计算X围内任意一位置(x，y)，可选取与其最接近的n〔n≤N〕个观测点进展观测，获得这n个观测点实测的降雨量，记为p1，p2，…，pn。通过最小二乘法拟合成一个二次曲面，由此二次曲面可确定点〔x，y〕的降雨量。

假设二次曲面的方程为

px,ya1x2a2y2a3xya4xa5ya6(8)

n点误差之和为

R2 (9) ppii1n5 / 8

word 根据最小二乘法，要使n点误差之和R达到最小，如此有

R0 i=1，2，…，n (10) ai由(10)式可求出ai(i=1，2，…，n)，并代入(8)式中，即可求出(x，y)点的插值降雨量p。

5 数学模型的验证

一般情况下，暴雨所造成的内涝灾害用最大积水深度来衡量。因此，验证此数学模型的的方法是选择暴雨个例，代入该数学模型中来计算，比拟该模型的实测结果和计算结果，并分析两者的绝对误差和误差的分布特征。

该数学模型的计算X围选取为某某市三环线以内，由某某市的地形等特征，将长江等划分为河道型网络，城区东湖、南湖等水体划分为湖泊型网格，其余的如此划分为面积各异的不规如此的陆地型网格。其中，对于某某市的中心地区和易积水区选取较为稠密的网格，其他区域如此选取较为稀疏的网格。综上，某某市共划分2336个不规如此网格，河道、湖泊共占有218个。模型选择2004年7月17~19日的特大暴雨个例，对模型进展调试，调试结果见表1。在调试的根底上选择2006年10月15~17日暴雨资料进展验证〔表2〕。模型调试、验证结果趋势正确，大局部精度符合要求。

表1 最大水深的调试结果

积水位置 竹叶山涵洞 姑嫂树铁路桥 琴台大道体育馆 琴台大道红太阳建材市场 晒湖涵洞 解放大道下延线铁路桥 长虹桥 21号公路石化南门 徐东路汪家墩车站 书城路金地格林小城 雄楚大街石牌岭 三环线常青立交—古田立交 东某某路上岭小区 实测值/m 计算值/m 绝对误差/m 相对误差/% 6 / 8

word 前进一路 开展大道银墩路 鹦鹉大道江鹦路 九支沟涵洞 五支沟涵洞 四支沟南路

表2 最大深度的验证结果

积水位置 开展大道东方酒店至江汉分局 开展大道银墩路 团山路白宫酒店 四支沟南路 常青花园武商量贩 友谊大道工业路口 21号公路石化南门 长虹桥 姑嫂树兴业路口附近 姑嫂树红旗渠路口 姑嫂树常青三路口 实测值/m 计算值/m 绝对误差/m 相对误差/%

6 误差产生原因

〔1〕某某市降雨X围分布不均。在短时间内，暴雨往往在空间分布上存在差异较大。某某市的自动雨量站也仅有17个且分布不均，降雨的空间分布难以准确地反映出来。 〔2〕数学模型网格的误差。此数学模型网格的大小主要涉与到市区地形、建筑物和街道，但由于此模型存在计算时间和计算稳定性的要求，网格的尺度不能太小。网格尺度较大时，通常取地势的平均值为该区域的地形高程，且计算水深值往往不能反映局部低洼点处的实际的积水深度。

〔3〕排水管网的概化误差。某某市的排水管网是非常复杂的树枝状系统构成的，由于此模型中每个网格都含有假如干管道，为了便于计算，管道的参数取其概化值，所以会影响其计算精度。

〔4〕积水的观测误差。实测积水时往往在最深处测量水深，而局部积水深度往往不能

7 / 8

word 代表整个网格的平均积水的深度。 7 结论

本文给出了城市暴雨内涝数学模型的建设思路，并应用该模型验证了某某市的暴雨过程。结果明确：某某市的暴雨内涝数学模型根本可信，并且模型具有较好的适用性。主要结论如下：

〔1〕以明渠河道与市内地表的水流运动为主要研究对象，以平面的二维非恒定流根本方程为框架，结合一维明渠非恒定流方程的计算方法，专门对于小于离散网格尺度的河道或排水渠涌，设计此模型。

〔2〕通过城市暴雨内涝的数学模型了验证了某某市的暴雨过程，其结果明确：城市暴雨内涝数学模型所得的结果是可信、可靠的，故而说明了此模型具有较好的适用性。此数学模型在某某市经历多场暴雨的考验，从侧面更证明了此模型的适用性。

参考文献：

[1] 解以扬，李大鸣，李培彦等.城市暴雨内涝数学模型的研究与应用[J].气象科学，2005,16(3)：333-335.

[2] 谭维炎.计算浅水动力学[M].：清华大学，1998.62-79. [3] 汪德.计算水力学理论与应用[M].某某：河海大学，1989.4-123.

[4] 徐向阳，X俊，郝庆庆等.城市暴雨积水过程模拟[J].水科学进展，2003,14(1)：193-196. [5] 仇劲卫，李娜，程晓陶等.某某市城区暴雨沥涝仿真模拟系统[J].水利学报，2000(11)：34-42.

[6] 蔡亦钢.流体传输管道动力学[M].某某:某某大学，1990.15-120.

[7] 某某大学水力学与水文学教研室.水力学(下册)[M].：人民教育，1980.45-109. [8] X俊.城市雨洪模型研究[J].河海大学学报，1997,25(6)：21-24.

8 / 8