武汉市暴雨内涝数学模型的研究与应用 刘晓 （湖北工业大学，湖北，武汉，120330270） 摘要：暴雨内涝对城市的影响日益严重，为了城市能够更好的应对暴雨带来的冲击，本文以城市的街道路面与河道水流的运动为对象进行模拟，建立了武汉市暴雨内涝积水数学模型。模型以平面二维非恒定流基本方程和不规则网格划分技术为框架，采用简化分类处理的方法，将通道分为路面型、河道型以及特殊通道型，根据不同类型简化动量方程，求任一网格各个通道上的单宽流量。根据不规则网格的方法，按照武汉市的地形进行多边形计算网格的设计。介绍了数学模型在武汉市的应用和误差分析以及城市路面降雨量的计算。 关键词：城市暴雨 内涝灾害 数学模型 误差分析 武汉市 Research and Application of Wuhan Waterlogging Mathematical Model Liu Xiao (,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 computational 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 error analysis;Wuhan
1 引言 城市内涝是由于强降雨超过城市排水能力而产生的城市内积水的灾害。降雨强度大是造成城市内涝的原因客观原因，范围集中。降雨强度较大、时间较长有可能形成积水，降雨急的地方也可能形成积水。近年来，暴雨对城市内涝灾害的影响日益严峻，不仅加重了城市的防汛排涝的任务，同时对经济建设和城市居民的日常生活也造成了较大影响。因此，为了进一步提高城市排水的服务能力和水平，对城市内涝进行研究已然十分必要，从而减少暴雨对城市的影响。 1.1国内外的经验 巴黎的下水道始建于1854年，用时24年巴黎修建了600多公里的下水道。截止今日，巴黎的下水道总长度近2300公里，相当于巴黎到伊斯坦布尔的距离。巴黎的下水道位于地面以下50米。水道纵横交错，密如蛛网，下水道十分宽敞，中间是宽约3米的排水道，两侧是宽约1米的供检修人员通行的便道。巴黎的下水道按沟道大小可以分为小下水道、中下水道和排水渠三种。每天有120亿立方米的水经此管道排出，这些下水道犹如这座大城市的消化系统。因此，暴雨不会使巴黎的形成内涝。 青岛是目前被认为是管道设置最合理的城市，青岛率先采用雨水污水分流的排水模式。德式排污管道设计巧妙，管道为倒置的鸭蛋型，上大下小，下半部分呈“V”型，在流量相同的情况下，污水流速更快，即使污水流量比较小，污水依然能保持一定的速度，使得污水中所夹杂的污物被冲击的没有沉淀的机会，有效保证了管道的畅通。 1.2武汉市的现状 武汉湖泊的减少，大大减小了城市的调蓄能力，从而大大增加了城市内涝的发生率。1950年代，武汉市湖泊面积达到1581平方公里，城中湖有127个，80年代已经缩减到874平方公里。近30年来，又减少了228.9平方公里，城中湖也缩减为38个。与之相对应的却是，城市建成区总面积从220平方公里增加到475平方公里。天上的雨水没有减少，调蓄的湖泊在减少，城市建成区的面积却在增加，如果抽排能力没有加强，城市内涝就不可避免。 1.3城市暴雨内涝数学模型的基本原理及网格划分 以城市的河道水流与街道路面的运动为对象进行模拟，建立了城市暴雨内涝积水的数学模型。模型以平面二维非恒定流基本方程和不规则网格划分技术为主要方法。一般的数值模式计算，基本上采用规则矩形网格或三角形网格对计算区域进行划分，但是城市的局部地区往往地形复杂，河流、街道分布也各不相同，对于建筑物，尤其是防洪排涝工程在南、北方有着巨大的差异，对径流的形成影响很大。因此，根据城市暴雨径流的特征，模型选取无结构不规则网格。同时，根据城市各异的地貌特征可将网格设计为三边形、四边形或五边形，网格的各边则定义为通道。选择网格单元为研究体，取网格中心位置得水位为计算水位H，取网格周边的通道中点处得流量为计算流量Q。水位代表网格的平均值，流量也表示通道的平均值。此外，在时间上，流量和水位采取时间交错的计算方式。根据武汉市的地形设计计算网格。武汉市地形平缓，海河较宽河道设为河道型网格，较窄的明沟和河道为特殊通道。 2 基本方程 由高斯定理可将连续方程转化为
dtqLQAHHdtTiikkKkiiTidtTi2dt212 (1)
式中Hi——i单元的水深； iA——i单元的网格面积；
Qik——i单元的第k条通道上的流量； ikL——i单元的第k条通道长度；
qi——i单元的源； dt——T时刻时间的增量。 根据简化的动量方程的离散格式来计算Qik。对于不规则的网格，如果按动量方程严格计算任一方向通道的流量，则计算量太大，所以该模型中采取简化分类处理的方法。假设同一时段同一网格水位的变化不大，采取分类处理的方法去求任一网格各通道上的单宽流量。通道分为路面型、河道型、特殊通道型、其他型，由不同类型的简化动量方程，计算流量。 河道型通道对应的动量方程：
dtHQQngdtdLZZgHQidtTikdtTikikikTikTikidtTikdtTk37212i22Q (2)
式中 g——重力加速度； Zik1、Zik2——为i单元的第k条通道左右侧1、2断面的水位； dLik——i单元的第k条通道左右侧单元形心距离； nik——i单元的第k条通道上的糙率。 地面型通道对应的动量方程：
dtTdNkNjddtTdikdtTikdiTdiiijijLQLQAdtHq22/2Q1212dtTik (3)
对于有缺口堤或连续堤的通道，其流量采取宽顶堰流的计算公式，即 2/3ik2бQiikdtTikHgm (4)
式中Σik——i单元的第k条通道下游淹没系数； mik——i单元的第k条通道流量系数。 武汉市较窄的河道认为是特殊通道，它与其左右侧网格之间的流量，选取宽顶堰流公式(4)进行计算；对于河道的单宽流量Qs选择与河道型通道的动量方程式(2)相近的公式来计算；此外，特殊单元水深则按照以下的离散格式计算：
dtLQLQAdtHdtTdNkNjddtTdikdtTikdTddtTiijijiiiq22/2H1212
d (5)
式中Hdi——i特殊单元的水深； Adi——i特殊单元的面积； Qdij——i通道周边的单宽流量； Ldij——i通道周边的边长； qdi——i特殊单元上的降雨强度，一般选取其特殊节点处的降雨强度来作为这个单元的平均值。
3 武汉市的排水系统 武汉市的排水系统由水井、管道、泵站、明沟、闸门等组成，排水管道错综复杂，通过雨水井来承担和接收来自地面的水，通过泵、闸或淹没出流管道来实现管道水的排出。对于地下排水管网内的水流，一般选用一维非恒定流方程进行计算。地面单元与地下排水管网的水量交换，一般可以通过连续方程源汇项来进行控制，即可将源汇项分成单元排水强度和有效降雨强度两部分。单元的排水量可依据单元的最大排水能力、单元积水深度，以及降雨强度来确定：当单元最大排水能力大于降雨强度并且单元积水相对较浅时，降雨量将全部转化成排水量流入到地下管网，故单元排水强度则为降雨强度；反之，当单元最大排水能力小于降雨强度且单元的积水相对较深时，单元排水强度可达到最大的排水强度，单元最大排水强度与降雨强度之差则转化成地表径流。 对于计算管道内的流量选用一维非恒定流方程的离散形式：
dtRAQQngdtdLZZgAQpipidtTpidtTpiiTiTipidtTpi3/4212dtTpi22Q （6）
式中Qpi——计算管道的断面流量； dLi——相临的两网格的形心到通道中点的距离之和； Rpi——水力半径； 计算管道内的水深一般选用连续方程的离散形式: