顾及建筑物微观水文过程的城市内涝建模方法

1.本发明涉及地理信息领域，尤其涉及顾及建筑物微观水文过程的城市内涝建模方法。


背景技术：

2.近年来，我国城市内涝灾害频发，严重威胁到人们的生命及财产安全，给城市发展带来诸多困扰。为有效应对洪涝灾害及其不利影响，越来越多的工程与非工程措施在洪涝灾害管理中实施。城市内涝模拟作为一项有效的非工程性措施，通过城市降雨径流过程模拟获取一定降雨条件下地表与地下的径流情况，能有效辅助灾情预报及分析，在国内外被广泛应用。
3.建筑物是城市的重要组成部分，随着地表建筑物密度增加，进而影响城市原有的水文循环，且改变了城市原有的雨水汇流过程，对城市内涝灾害的产生具有重要影响。建筑物微观水文过程是由建筑物雨水口、建筑物排水管、地表雨水篦及地下雨水管网等形成的从屋顶到地表再到地下的立体化汇水网络，对地表渗透率、地表积水和径流形成具有较强的影响，且对城市内涝模拟具有重要作用。
4.然而，传统的城市内涝模拟缺乏对建筑物微观水文过程的深入剖析，忽略雨水从建筑物屋顶到地表再到地下的立体化汇水过程，缺乏对建筑物汇水单元的准确表达，未考虑建筑物屋顶特征对地表高程修正的影响，难以满足城市复杂环境下内涝灾害时空模拟及分析的需要，如何有效顾及建筑物微观水文过程对内涝产生的影响已经成为城市背景下内涝模拟的重要内容。
5.综上所述，目前缺少一种顾及建筑物微观水文过程的城市内涝建模方法，来提高城市内涝模拟的准确性与精度，为内涝治理措施制定提供科学依据。


技术实现要素：

6.本发明所要解决的技术问题是针对背景技术的不足提供了顾及建筑物微观水文过程的城市内涝建模方法，本发明针对建筑物微观水文过程对内涝形成过程的影响，创新地构建了从屋顶到地表再到地下的立体化汇水网络，实现雨水从建筑物屋顶到地表的排水过程模拟，改进了传统城市内涝模拟缺乏建筑物汇水过程的不足之处，进一步满足城市水文单元复杂环境下内涝灾害时空模拟的需要。
7.本发明为解决上述技术问题采用以下技术方案：
8.顾及建筑物微观水文过程的城市内涝建模方法，包括如下步骤：
9.步骤s1：首先对建筑物屋顶水文过程深入剖析，针对建筑物屋顶雨水口及排水管网进行处理，构建基于建筑物排水设施的立体化汇水网络；
10.步骤s2：分析建筑物立体空间结构特征，在升高建筑物所在的地面高程基础上，顾及屋顶矮墙进一步修正地表高程；
11.步骤s3：结合地形、道路、建筑物及各类人工排水设施对城市汇流过程的影响，遵
循从大到小、逐层细化的原则对城市地表雨水汇水区进行划分；
12.步骤s4：基于地表与地下管网水流过程分析，提出地表与地下管网水流交换耦合方法，并对地表径流过程与地下管网水流输送过程建模，构建城市内涝模型。
13.进一步的，步骤s1中包括以下步骤：
14.步骤s11：将建筑物屋顶的雨水口概化成真实的雨水篦，进行属性赋值并确定雨水口对应的雨水点坐标；
15.步骤s12：对概化而成的雨水点的雨水管点点号、管点埋深、管点高程值等属性进行赋值；将管点埋深设为0，其管点高程值设为雨水口对应地面高程值与建筑物高度之和；
16.步骤s13：将建筑物屋顶雨水口概化而成的雨水点与建筑物周边最近雨水点相连接，建立建筑物汇水网络与城市雨水管网连通关系；
17.步骤s14：通过对排水管网进行实际调查与采集，对雨水管网的管段名称、管段起点、管段终点、管段起点埋深、管段终点埋深以及管径等属性进行赋值，形成与现实相符合的建筑物立体化汇水网络。
18.进一步的，步骤s2中还包括以下步骤：
19.s21：将建筑物边缘与栅格单元边界对齐，以便屋顶标高可以正确描述建筑物的位置；
20.s22：将建筑物面要素转成同样分辨率的栅格数据，建筑物高度概化为栅格高度；
21.s23：将建筑物对应栅格高度叠加到原始dem上，从而完成dem与建筑物信息的融合；
22.s24：由于建筑物屋顶周围通常建设有矮墙，通过获取建筑物边界对应的栅格高程值进行修正，使得建筑物边界高程值满足实际建筑物屋顶边界情况。
23.进一步的，步骤s3中包括以下步骤：
24.步骤s31：对待划分区域的dem数据进行流向提取、伪洼地填充、汇流累积量计算、天然水系提取及汇水区生成操作，得到基于地形划分的汇水区；
25.步骤s32：自动提取道路主干道中心线和建筑物轮廓线，对已有汇水区进一步划分；
26.步骤s33：基于上述子汇水区，筛选出含有管点且管点数大于1的子汇水面要素，基于汇水区中的管点进行泰森多边形的划分；
27.步骤s34：根据s1中构建的建筑物立体化汇水网络，将建筑物屋顶的雨水口作为雨水点，基于雨水口概化的雨水点对建筑物子汇水单元进行划分。
28.进一步的，步骤s41：通过提取dem栅格的行列号，判断并计算出对应栅格的坐标值(x1，y1)；
29.步骤s42：提取地下管点的坐标值(x2，y2)，根据该坐标遍历研究区域的dem栅格，找到存在雨水管点的栅格，确定对应关系；
30.步骤s43：使用一维地下管网模型计算每个管网节点对应的水头，水头减去管点高程即为雨水管点对应水深h
1d
；再通过雨水管点坐标找到地表dem对应的栅格，利用二维地表漫流模型计算出对应dem栅格上的水深h
2d
；
31.步骤s44：判断是否满足条件z
2d
≤h
1d
≤h
2d
，z
2d
为地表高程；若满足，则选用孔口出流公式计算交换流量，否则进入步骤s45；其中孔口出流公式如下：
[0032][0033]
其中，qv表示流量，cq表示流量系数，a表示孔口断面面积，g表示重力加速度，h1-h2表示水位差；
[0034]
步骤s45：判断管点所在dem单元栅格高程值是否小于所有周围dem栅格单元对应高程值，即当前雨水完全覆盖雨水井，若满足，则选用孔口出流公式计算交换流量，否则进入步骤s46；
[0035]
步骤s46：选用堰流公式进行计算交换流量，堰流公式如下：
[0036][0037]
其中，q表示流量，m表示流量系数，b表示溢流宽度，g表示重力加速度，h0表示堰上总水头；
[0038]
步骤s47：选用horton模型进行下渗量计算，非线性水库模型对地表径流进行模拟，通过输入子汇水区的面积、宽度、坡度以及地表曼宁系数、滞蓄量等参数，并联立曼宁公式和连续方程进行求解，完成地表径流过程建模；其中horton模型经验公式如下：
[0039]fp
＝f
∞
(f
0-f
∞
)e-αt
[0040]
其中，f
p
为土壤的下渗率(m2/h)，f
∞
为稳定下渗率(m2/h)，f0为初始下渗率(m2/h)，t为降雨历时，α为衰减指数；
[0041]
其中连续方程如下：
[0042][0043]
其中，d为水深，v＝a
×
d为地表积水量(m3)，a为排水区表面积(m2)，i*为净雨，q为流出量(m2/s)；
[0044]
其中，曼宁方程如下：
[0045][0046]
其中，s为子流域的宽度，n为曼宁系数，w为子流域漫流宽度(m)，d
p
为地表滞蓄水深(m)；
[0047]
步骤s48：基于上述地表汇流模拟预测得到的指定地点的排放水流出口，结合动力波水流计算方法，求解一维圣维南方程，算得管道中的流速和水深，并利用栅格上的存储单元概念来模拟地下管网水流输送过程；其中圣维南方程由一维全浅水方程的动量方程和连续性方程组成：
[0048][0049][0050]
其中，q
x
是x笛卡尔方向的体积流量，a为水流的横截面积，h为水深，z为高层，g为重力，n为曼宁的摩擦系数，r为水力半径，t为时间，x为x笛卡尔方向的距离。
[0051]
本发明采用以上技术方案与现有技术相比，具有以下技术效果：
[0052]
1、通过对建筑物微观水文过程的深入剖析，充分考虑雨水从建筑物屋顶到地表再
到地下的立体化汇水过程，实现了对建筑物汇水单元的准确表达，改善了传统城市内涝模拟对建筑物微观水文过程研究的不足，进一步满足城市复杂环境下内涝灾害时空模拟及分析的需要。
[0053]
2、在此基础上，结合雨水汇水区划分模块、地下管网模型与地表漫流模型耦合模块，构建了完整的城市内涝模型，不仅丰富了城市内涝时空模拟的内涵，还可以促进对微观水文过程的理解，同时能够为社区尺度内涝灾害管理提供科学指导。
附图说明
[0054]
图1为本实施例的流程图。
[0055]
图2为本实施例的具体应用框图。
[0056]
图3为本实施例的基于建筑物排水设施构建的立体化汇水网络示意图。
[0057]
图4中a为建筑物屋顶雨水排水示意图，b为建筑物汇水网络与城市雨水管网连通示意图。
[0058]
图5为本实施例的地表与地下管网水流交换耦合方法算法流程图。
具体实施方式
[0059]
下面结合附图对本发明的技术方案做进一步的详细说明：
[0060]
本发明的描述中，需要理解的是，术语“左侧”、“右侧”、“上部”、“下部”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，“第一”、“第二”等并不表示零部件的重要程度，因此不能理解为对本发明的限制。本实施例中采用的具体尺寸只是为了举例说明技术方案，并不限制本发明的保护范围。
[0061]
如图1所示，本发明公开一种基于顾及建筑物微观水文过程的城市内涝建模方法，包括如下步骤：
[0062]
步骤s1：首先对建筑物屋顶水文过程深入剖析，针对建筑物屋顶雨水口及排水管网进行处理，构建基于建筑物排水设施的立体化汇水网络；
[0063]
步骤s2：分析建筑物立体空间结构特征，在升高建筑物所在的地面高程基础上，顾及屋顶矮墙进一步修正地表高程；
[0064]
步骤s3：结合地形、道路、建筑物及各类人工排水设施对城市汇流过程的影响，遵循从大到小、逐层细化的原则对城市地表雨水汇水区进行划分；
[0065]
步骤s4：基于地表与地下管网水流过程分析，提出地表与地下管网水流交换耦合方法，并对地表径流过程与地下管网水流输送过程建模，构建城市内涝模型。
[0066]
本实施例中，使用某典型内涝易发区域(面积约为447031m2)的雨水管网数据、降雨与积水数据、建筑物数据、数字高程数据(dem)、路网数据等进行建模。区域为城市核心城区，建筑面积较大，地表硬化率相对较高，地表建筑物微观排水过程对内涝过程具有重要影响。
[0067]
本实施例的具体实施步骤如图2所示：
[0068]
步骤s1：首先对建筑物屋顶水文过程深入剖析，针对建筑物屋顶雨水口及排水管网进行处理，构建基于建筑物排水设施的立体化汇水网络，其具体步骤如下：
[0069]
步骤s11：将建筑物屋顶的雨水口概化成真实的雨水篦，进行属性赋值并确定雨水口对应的雨水点坐标；
[0070]
步骤s12：对概化而成的雨水点的雨水管点点号、管点埋深、管点高程值等属性进行赋值。将管点埋深设为0，其管点高程值设为雨水口对应地面高程值与建筑物高度之和；
[0071]
步骤s13：将建筑物屋顶雨水口概化而成的雨水点与建筑物周边最近雨水点相连接，建立建筑物汇水网络与城市雨水管网连通关系；
[0072]
步骤s14：通过对排水管网进行实际调查与采集，对雨水管网的管段名称、管段起点、管段终点、管段起点埋深、管段终点埋深以及管径等属性进行赋值，形成与现实相符合的建筑物立体化汇水网络，如图3、图4所示。
[0073]
步骤s2：分析建筑物立体空间结构特征，在升高建筑物所在的地面高程基础上，顾及屋顶矮墙进一步修正地表高程，其具体步骤如下：
[0074]
步骤s21：将建筑物边缘与栅格单元边界对齐，以便屋顶标高可以正确描述建筑物的位置；
[0075]
步骤s22：将建筑物面要素转成同样分辨率的栅格数据，建筑物高度概化为栅格高度；
[0076]
步骤s23：将建筑物对应栅格高度叠加到原始dem上，从而完成dem与建筑物信息的融合；
[0077]
步骤s24：由于建筑物屋顶周围通常建设有矮墙，通过获取建筑物边界对应的栅格高程值进行修正，使得建筑物边界高程值满足实际建筑物屋顶边界情况。
[0078]
步骤s3：结合地形、道路、建筑物及各类人工排水设施对城市汇流过程的影响，遵循从大到小、逐层细化的原则对城市地表雨水汇水区进行划分，其具体步骤如下：
[0079]
步骤s31：对待划分区域的dem数据进行流向提取、伪洼地填充、汇流累积量计算、天然水系提取及汇水区生成操作，得到基于地形划分的汇水区；
[0080]
步骤s32：自动提取道路主干道中心线和建筑物轮廓线，对已有汇水区进一步划分；
[0081]
步骤s33：基于上述子汇水区，筛选出含有管点且管点数大于1的子汇水面要素，基于汇水区中的管点进行泰森多边形的划分；
[0082]
步骤s34：根据s1中构建的建筑物立体化汇水网络，将建筑物屋顶的雨水口作为雨水点，基于雨水口概化的雨水点对建筑物子汇水单元进行划分。
[0083]
步骤s4：基于地表与地下管网水流过程分析，提出地表与地下管网水流交换耦合方法，并对地表径流过程与地下管网水流输送过程建模，构建城市内涝模型，其具体实施步骤如下如图5所示：
[0084]
步骤s41：通过提取dem栅格的行列号，判断并计算出对应栅格的坐标值(x1，y1)；
[0085]
步骤s42：提取地下管点的坐标值(x2，y2)，根据该坐标遍历研究区域的dem栅格，找到存在雨水管点的栅格，确定对应关系；
[0086]
步骤s43：使用一维地下管网模型计算每个管网节点对应的水头，水头减去管点高程即为雨水管点对应水深h
1d
。再通过雨水管点坐标找到地表dem对应的栅格，利用二维地表漫流模型计算出对应dem栅格上的水深h
2d
；
[0087]
步骤s44：判断是否满足条件z
2d
≤h
1d
≤h
2d
(z
2d
为地表高程)，若满足，则选用孔口出
流公式计算交换流量，否则进入步骤s45；其中孔口出流公式如下：
[0088][0089]
其中，qv表示流量，cq表示流量系数，a表示孔口断面面积，g表示重力加速度，h1-h2表示水位差。
[0090]
步骤s45：判断管点所在dem单元栅格高程值是否小于所有周围dem栅格单元对应高程值(即当前雨水完全覆盖雨水井)，若满足，则选用孔口出流公式计算交换流量，否则进入步骤s46；
[0091]
步骤s46：选用堰流公式进行计算交换流量，堰流公式如下：
[0092][0093]
其中，q表示流量，m表示流量系数，b表示溢流宽度，g表示重力加速度，h0表示堰上总水头。
[0094]
步骤s47：选用horton模型进行下渗量计算，非线性水库模型对地表径流进行模拟，通过输入子汇水区的面积、宽度、坡度以及地表曼宁系数、滞蓄量等参数，并联立曼宁公式和连续方程进行求解，完成地表径流过程建模；其中horton模型经验公式如下：
[0095]fp
＝f
∞
(f
0-f
∞
)e-αt
[0096]
其中，f
p
为土壤的下渗率(m2/h)，f
∞
为稳定下渗率(m2/h)，f0为初始下渗率(m2/h)，t为降雨历时，α为衰减指数。
[0097]
其中连续方程如下：
[0098][0099]
其中，d为水深，v＝a
×
d为地表积水量(m3)，a为排水区表面积(m2)，i*为净雨，q为流出量(m2/s)。
[0100]
其中，曼宁方程如下：
[0101][0102]
其中，s为子流域的宽度，n为曼宁系数，w为子流域漫流宽度(m)，d
p
为地表滞蓄水深(m)。
[0103]
步骤s48：基于上述地表汇流模拟预测得到的指定地点的排放水流出口，结合动力波水流计算方法，求解一维圣维南方程，算得管道中的流速和水深，并利用栅格上的存储单元概念来模拟地下管网水流输送过程。其中圣维南方程由一维全浅水方程的动量方程和连续性方程组成：
[0104][0105][0106]
其中，q
x
是x笛卡尔方向的体积流量，a为水流的横截面积，h为水深，z为高层，g为重力，n为曼宁的摩擦系数，r为水力半径，t为时间，x为x笛卡尔方向的距离。
[0107]
本技术领域技术人员可以理解的是，除非另外定义，这里使用的所有术语(包括技
术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是，诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义，并且除非像这里一样定义，不会用理想化或过于正式的含义来解释。
[0108]
以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。上面对本发明的实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以再不脱离本发明宗旨的前提下做出各种变化。