一种基于元胞自动机的城市地表径流与管网汇流耦合方法技术

本发明专利技术公开了一种基于元胞自动机的城市地表径流与管网汇流耦合方法，通过构建地表元胞地理场景，概化城市地下排水管网，以雨水井口为耦合节点，进行地表径流与管网汇流的耦合计算。本发明专利技术首先进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件，再进行管网汇流计算获得当前时间步长管道流量和节点水头，根据降雨过程中的地表入流、节点溢流、退水回流三种不同的流量过程确定当前时刻节点净流量，并代入管道流量基本方程和排水管网节点控制方程，获得下一时刻管道流量和节点水头。本发明专利技术方法耦合了管道汇流和基于元胞自动机的地表径流过程，实现地表和管网水量交换的动态模拟，更加符合真实的城市降雨径流过程。

【技术实现步骤摘要】
一种基于元胞自动机的城市地表径流与管网汇流耦合方法
本专利技术涉及城市水文和管网排水领域，尤其涉及一种基于元胞自动机的城市地表径流与管网汇流耦合方法。
技术介绍
随着城市化的快速发展，城市不透水面不断扩大，再加上以前建设的城市排水系统排水能力不足，在暴雨事件发生时易形成城市内涝，影响社会生活和经济生产。城市雨洪模型能够模拟降雨事件下的城市管网排水情况，是进行城市管网规划建设重要决策工具之一。目前的城市雨洪模型主要为一维排水管网模型，难以模拟排水能力不足情境下的溢流水量在地表的二维运动，无法确定内涝的影响范围。元胞自动机是一种时间、空间、状态都离散的网格动力学模型，具有模拟复杂系统时空演化过程的能力，在城市复杂地形条件下的地表径流过程模拟中具有一定的优势。将地下管网排水模型与基于元胞自动机的城市地表径流模型耦合，实现地上地下水量的交换，模拟城市内涝过程，可为城市地下排水管网规划设计提供参考。因此，考虑降雨事件周期下的水流地上地下运动过程，构建地表元胞地理场景，概化地下城市排水管网，以雨水井口作为地表地下水量交换节点，实现基于元胞自动机的地表径流与地下管网汇流的耦合计算，模拟暴雨事件下真实城市降雨径流过程，确定城市内涝产生时长以及地表积水淹没范围，对于城市规划和管网设计具有重要的指导意义。
技术实现思路
专利技术目的：本专利技术旨在提供一种基于元胞自动机的城市地表径流与管网汇流耦合方法，该方法可以较好的描述城市降雨过程中地表入流、节点溢流以及回水退流过程，更加符合真实降雨径流情况。r>技术方案：本专利技术所述的一种基于元胞自动机的城市地表径流与管网汇流耦合方法，主要包括如下步骤：(1)构建地表元胞地理场景，根据研究区下垫面类型、DEM高程和随时间变化的降雨量，获得每个元胞的下垫面类型、高程和随时间变化的降雨量，定义模型地理场景的元胞大小、邻域关系以及模型边界；(2)概化城市地下排水管网，根据研究区排水管网建设情况，获得排水管网管道参数、节点参数以及排水出口位置，并将雨水井口作为地表和管网水流的交换节点；(3)以雨水井口为耦合节点，进行地表径流与管网汇流耦合计算，包括：进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件；进行管网汇流计算获得当前时间步长管道流量和节点水头，再根据降雨过程中的地表入流、节点溢流以及退水回流三种不同的过程确定当前时刻节点净流量；将当前时刻管道流量、节点水头与节点净流量代入管道流量基本方程和排水管网节点控制方程，获得下一时刻管道流量和节点水头，再确定节点净流量，重复计算，直至达到模拟总时长，结束计算。作为优选，所述步骤(3)中进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件包括：利用初损后损法进行产流计算获得每个元胞水深；比较中心元胞与邻居元胞的水位，排除高于平均水位的邻居元胞，确定可以进行水量转移的邻居元胞；在产流和流向确定的基础上，采用曼宁公式计算一个时间步长内中心元胞可以转移的水量，同时根据元胞间水位差得出中心元胞可以向邻居元胞转移的潜在最大水量，将二者相比，得出一个时间步长内中心元胞转移的真实水量；采用侧收缩堰流量计算公式计算各邻居元胞流入雨水井口节点的流量，并求和得到地表流入雨水井口节点的总流量。作为优选，所述步骤(3)中通过管道流量基本方程和排水管网节点控制方程计算管道流量和节点水头；其中管道流量基本方程和节点控制方程的有限差分形式表示为：其中，为t时刻的管道流量，为t时刻的节点净流量，为t时刻的节点水头，Δt为模型时间步长，为管道平均流速，A1、A2分别为节点上、下游管道断面面积，ΔA为一个时间步长内平均断面面积变化值，L为管道的长度，为管道最宽处的断面面积，g为重力加速度，H1、H2分别为管道上、下游节点水头，k为损失系数，为平均水力半径，Ask为节点面积。作为优选，所述步骤(3)中，处于降雨初期水流由地表流入管网的地表入流过程时，节点净流量计算方式为：其中，为t时刻节点地表流入流量，该流量为采用侧收缩堰流量计算公式计算得到的节点元胞的各邻居元胞流入流量的总和，和分别为t时刻节点上、下游管道流量。作为优选，所述步骤(3)中，处于降雨强度增大后的管网水流溢出地表的节点溢流过程时，节点净流量计算方式为：当节点元胞转移出去的水量Volttrans_out小于其邻居元胞转移进来的水量Volttrans_in时，其中，和分别为t时刻节点上、下游管道流量；当节点元胞转移出去的水量Volttrans_out大于其邻居元胞转移进来的水量Volttrans_in时，其中作为优选，所述步骤(3)中，处于降雨后期水流重新进入管道的退水回流过程时，节点流量计算方式为：退水回流过程中，节点水头高于地表高程时，按节点溢流过程计算节点净流量；地表积水退去，节点水头低于地表高程时，按地表入流过程计算节点净流量。有益效果：与现有技术相比，本专利技术公开的基于元胞自动机的城市地表径流与管网汇流耦合方法具有以下有益效果：1、本专利技术将元胞自动机理论应用于城市降雨径流模型研究中，用其“自下而上”的建模方式取代繁复的方程，相比于传统的集总式水文模型，更有利于对城市水文过程机理进行分析研究。2、本专利技术基于元胞自动机模型耦合城市地表径流和管网汇流模型，构建了可以动态模拟地表径流过程、溢流内涝过程以及退水回流过程的城市降雨径流模型。附图说明图1是本专利技术实施例的流程图。图2是本专利技术实施例中元胞边界条件示意图。图3是本专利技术实施例中地表径流模型和排水管网模型耦合节点图。图4是本专利技术实施例中地表径流与管网汇流在节点处的流量交换示意图。图5是本专利技术实施例中降雨径流过程模拟示意图。具体实施方式如图1所示，本专利技术实施例公开的一种基于元胞自动机的城市地表径流与管网汇流耦合方法，包括以下步骤：(1)构建地表元胞地理场景，根据研究区下垫面类型、DEM高程和随时间变化的降雨量，获得每个元胞的下垫面类型、高程和随时间变化的降雨量，定义地理场景模型元胞大小、邻域关系以及模型边界。该步骤具体包括：(1.1)获取研究区下垫面类型和数字高程数据(DEM)构建元胞基础地理场景，获得每个元胞的下垫面类型(例如路面、屋面、草地等)和高程，设置模型边界条件和模型出口，定义元胞大小、曼宁糙率系数、坐标、邻域关系、模型时间步长以及模拟时长。元胞边界设置如图2所示，本专利技术中元胞的邻域关系为Moore型八邻域，边界元胞设置为极大值从而保证水量无法进入。(1.2)获取研究区域随时间变化的降雨量，可以是真实情况下的降雨过程，也可以是根据防洪标准假设的不同情境下的降雨过程，进而获得每个元胞不同时刻的降雨量。(2)概化城市地下排水管网，根据研究区排水管网建设情况，获得排水管网的出口位置、管道长度、管径、管道数量、节点数量及其所在位置，并将雨水井口作为地表和管网水流的交换节点，将本文档来自技高网...

【技术保护点】
1.一种基于元胞自动机的城市地表径流与管网汇流耦合方法，其特征在于，包括如下步骤：/n(1)构建地表元胞地理场景，根据研究区下垫面类型、DEM高程和随时间变化的降雨量，获得每个元胞的下垫面类型、高程和随时间变化的降雨量，定义模型地理场景的元胞大小、邻域关系以及模型边界；/n(2)概化城市地下排水管网，根据研究区排水管网建设情况，获得排水管网管道参数、节点参数以及排水出口位置，并将雨水井口作为地表和管网水流的交换节点；/n(3)以雨水井口为耦合节点，进行地表径流与管网汇流耦合计算，包括：进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件；进行管网汇流计算获得当前时间步长管道流量和节点水头，再根据降雨过程中的地表入流、节点溢流以及退水回流三种不同的过程确定当前时刻节点净流量；将当前时刻管道流量、节点水头与节点净流量代入管道流量基本方程和排水管网节点控制方程，获得下一时刻管道流量和节点水头，再确定节点净流量，重复计算，直至达到模拟总时长，结束计算。/n

【技术特征摘要】
1.一种基于元胞自动机的城市地表径流与管网汇流耦合方法，其特征在于，包括如下步骤：
(1)构建地表元胞地理场景，根据研究区下垫面类型、DEM高程和随时间变化的降雨量，获得每个元胞的下垫面类型、高程和随时间变化的降雨量，定义模型地理场景的元胞大小、邻域关系以及模型边界；
(2)概化城市地下排水管网，根据研究区排水管网建设情况，获得排水管网管道参数、节点参数以及排水出口位置，并将雨水井口作为地表和管网水流的交换节点；
(3)以雨水井口为耦合节点，进行地表径流与管网汇流耦合计算，包括：进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件；进行管网汇流计算获得当前时间步长管道流量和节点水头，再根据降雨过程中的地表入流、节点溢流以及退水回流三种不同的过程确定当前时刻节点净流量；将当前时刻管道流量、节点水头与节点净流量代入管道流量基本方程和排水管网节点控制方程，获得下一时刻管道流量和节点水头，再确定节点净流量，重复计算，直至达到模拟总时长，结束计算。


2.根据权利要求1所述的基于元胞自动机的城市地表径流与管网汇流耦合方法，其特征在于，所述步骤(3)中进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件包括：
利用初损后损法进行产流计算获得每个元胞水深；
比较中心元胞与邻居元胞的水位，排除高于平均水位的邻居元胞，确定可以进行水量转移的邻居元胞；
在产流和流向确定的基础上，采用曼宁公式计算一个时间步长内中心元胞可以转移的水量，同时根据元胞间水位差得出中心元胞可以向邻居元胞转移的潜在最大水量，将二者相比，得出一个时间步长内中心元胞转移的真实水量；
采用侧收缩堰流量计算公式计算各邻居元胞流入雨水井口节点的流量，并求和得到地表流入雨水井口节点的总流量。


3.根据权利要求1所述的基于元胞自动机的城市地表径流与管网汇流耦合方法，其特征在于，所述...

【专利技术属性】
技术研发人员：陈求稳，陈诚，张建云，李港，何湖滨，何梦男，王智源，严晗璐，林育青，刘东升，黄玉，
申请(专利权)人：水利部交通运输部国家能源局南京水利科学研究院，
类型：发明
国别省市：江苏;32