CN112052545A - 一种基于元胞自动机的城市地表径流与管网汇流耦合方法 - Google Patents
一种基于元胞自动机的城市地表径流与管网汇流耦合方法 Download PDFInfo
- Publication number
- CN112052545A CN112052545A CN202010860835.XA CN202010860835A CN112052545A CN 112052545 A CN112052545 A CN 112052545A CN 202010860835 A CN202010860835 A CN 202010860835A CN 112052545 A CN112052545 A CN 112052545A
- Authority
- CN
- China
- Prior art keywords
- node
- flow
- water
- pipe network
- pipeline
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000001413 cellular effect Effects 0.000 title claims abstract description 38
- 238000010168 coupling process Methods 0.000 title claims abstract description 23
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 143
- 238000000034 method Methods 0.000 claims abstract description 74
- 230000008569 process Effects 0.000 claims abstract description 53
- 238000004364 calculation method Methods 0.000 claims abstract description 45
- 230000008859 change Effects 0.000 claims abstract description 11
- 238000004519 manufacturing process Methods 0.000 claims abstract description 7
- 230000008878 coupling Effects 0.000 claims description 16
- 238000005859 coupling reaction Methods 0.000 claims description 16
- 238000011144 upstream manufacturing Methods 0.000 claims description 13
- 238000011160 research Methods 0.000 claims description 11
- 230000008602 contraction Effects 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 239000002352 surface water Substances 0.000 claims description 2
- 230000036962 time dependent Effects 0.000 claims description 2
- 238000005094 computer simulation Methods 0.000 abstract 1
- 210000004027 cell Anatomy 0.000 description 55
- 238000010586 diagram Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 2
- 230000008595 infiltration Effects 0.000 description 2
- 238000001764 infiltration Methods 0.000 description 2
- 238000010992 reflux Methods 0.000 description 2
- 102000020897 Formins Human genes 0.000 description 1
- 108091022623 Formins Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000003888 boundary cell Anatomy 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Sewage (AREA)
Abstract
本发明公开了一种基于元胞自动机的城市地表径流与管网汇流耦合方法,通过构建地表元胞地理场景,概化城市地下排水管网,以雨水井口为耦合节点,进行地表径流与管网汇流的耦合计算。本发明首先进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件,再进行管网汇流计算获得当前时间步长管道流量和节点水头,根据降雨过程中的地表入流、节点溢流、退水回流三种不同的流量过程确定当前时刻节点净流量,并代入管道流量基本方程和排水管网节点控制方程,获得下一时刻管道流量和节点水头。本发明方法耦合了管道汇流和基于元胞自动机的地表径流过程,实现地表和管网水量交换的动态模拟,更加符合真实的城市降雨径流过程。
Description
技术领域
本发明涉及城市水文和管网排水领域,尤其涉及一种基于元胞自动机的城市地表径流与管网汇流耦合方法。
背景技术
随着城市化的快速发展,城市不透水面不断扩大,再加上以前建设的城市排水***排水能力不足,在暴雨事件发生时易形成城市内涝,影响社会生活和经济生产。城市雨洪模型能够模拟降雨事件下的城市管网排水情况,是进行城市管网规划建设重要决策工具之一。目前的城市雨洪模型主要为一维排水管网模型,难以模拟排水能力不足情境下的溢流水量在地表的二维运动,无法确定内涝的影响范围。
元胞自动机是一种时间、空间、状态都离散的网格动力学模型,具有模拟复杂***时空演化过程的能力,在城市复杂地形条件下的地表径流过程模拟中具有一定的优势。将地下管网排水模型与基于元胞自动机的城市地表径流模型耦合,实现地上地下水量的交换,模拟城市内涝过程,可为城市地下排水管网规划设计提供参考。
因此,考虑降雨事件周期下的水流地上地下运动过程,构建地表元胞地理场景,概化地下城市排水管网,以雨水井口作为地表地下水量交换节点,实现基于元胞自动机的地表径流与地下管网汇流的耦合计算,模拟暴雨事件下真实城市降雨径流过程,确定城市内涝产生时长以及地表积水淹没范围,对于城市规划和管网设计具有重要的指导意义。
发明内容
发明目的:本发明旨在提供一种基于元胞自动机的城市地表径流与管网汇流耦合方法,该方法可以较好的描述城市降雨过程中地表入流、节点溢流以及回水退流过程,更加符合真实降雨径流情况。
技术方案:本发明所述的一种基于元胞自动机的城市地表径流与管网汇流耦合方法,主要包括如下步骤:
(1)构建地表元胞地理场景,根据研究区下垫面类型、DEM高程和随时间变化的降雨量,获得每个元胞的下垫面类型、高程和随时间变化的降雨量,定义模型地理场景的元胞大小、邻域关系以及模型边界;
(2)概化城市地下排水管网,根据研究区排水管网建设情况,获得排水管网管道参数、节点参数以及排水出口位置,并将雨水井口作为地表和管网水流的交换节点;
(3)以雨水井口为耦合节点,进行地表径流与管网汇流耦合计算,包括:进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件;进行管网汇流计算获得当前时间步长管道流量和节点水头,再根据降雨过程中的地表入流、节点溢流以及退水回流三种不同的过程确定当前时刻节点净流量;将当前时刻管道流量、节点水头与节点净流量代入管道流量基本方程和排水管网节点控制方程,获得下一时刻管道流量和节点水头,再确定节点净流量,重复计算,直至达到模拟总时长,结束计算。
作为优选,所述步骤(3)中进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件包括:
利用初损后损法进行产流计算获得每个元胞水深;
比较中心元胞与邻居元胞的水位,排除高于平均水位的邻居元胞,确定可以进行水量转移的邻居元胞;
在产流和流向确定的基础上,采用曼宁公式计算一个时间步长内中心元胞可以转移的水量,同时根据元胞间水位差得出中心元胞可以向邻居元胞转移的潜在最大水量,将二者相比,得出一个时间步长内中心元胞转移的真实水量;
采用侧收缩堰流量计算公式计算各邻居元胞流入雨水井口节点的流量,并求和得到地表流入雨水井口节点的总流量。
作为优选,所述步骤(3)中通过管道流量基本方程和排水管网节点控制方程计算管道流量和节点水头;其中管道流量基本方程和节点控制方程的有限差分形式表示为:
其中,为t时刻的管道流量,为t时刻的节点净流量,为t时刻的节点水头,Δt为模型时间步长,为管道平均流速,A1、A2分别为节点上、下游管道断面面积,ΔA为一个时间步长内平均断面面积变化值,L为管道的长度,为管道最宽处的断面面积,g为重力加速度,H1、H2分别为管道上、下游节点水头,k为损失系数,为平均水力半径,Ask为节点面积。
作为优选,所述步骤(3)中,处于降雨初期水流由地表流入管网的地表入流过程时,节点净流量计算方式为:其中,为t时刻节点地表流入流量,该流量为采用侧收缩堰流量计算公式计算得到的节点元胞的各邻居元胞流入流量的总和,和分别为t时刻节点上、下游管道流量。
作为优选,所述步骤(3)中,处于降雨强度增大后的管网水流溢出地表的节点溢流过程时,节点净流量计算方式为:
作为优选,所述步骤(3)中,处于降雨后期水流重新进入管道的退水回流过程时,节点流量计算方式为:退水回流过程中,节点水头高于地表高程时,按节点溢流过程计算节点净流量;地表积水退去,节点水头低于地表高程时,按地表入流过程计算节点净流量。
有益效果:与现有技术相比,本发明公开的基于元胞自动机的城市地表径流与管网汇流耦合方法具有以下有益效果:1、本发明将元胞自动机理论应用于城市降雨径流模型研究中,用其“自下而上”的建模方式取代繁复的方程,相比于传统的集总式水文模型,更有利于对城市水文过程机理进行分析研究。2、本发明基于元胞自动机模型耦合城市地表径流和管网汇流模型,构建了可以动态模拟地表径流过程、溢流内涝过程以及退水回流过程的城市降雨径流模型。
附图说明
图1是本发明实施例的流程图。
图2是本发明实施例中元胞边界条件示意图。
图3是本发明实施例中地表径流模型和排水管网模型耦合节点图。
图4是本发明实施例中地表径流与管网汇流在节点处的流量交换示意图。
图5是本发明实施例中降雨径流过程模拟示意图。
具体实施方式
如图1所示,本发明实施例公开的一种基于元胞自动机的城市地表径流与管网汇流耦合方法,包括以下步骤:
(1)构建地表元胞地理场景,根据研究区下垫面类型、DEM高程和随时间变化的降雨量,获得每个元胞的下垫面类型、高程和随时间变化的降雨量,定义地理场景模型元胞大小、邻域关系以及模型边界。该步骤具体包括:
(1.1)获取研究区下垫面类型和数字高程数据(DEM)构建元胞基础地理场景,获得每个元胞的下垫面类型(例如路面、屋面、草地等)和高程,设置模型边界条件和模型出口,定义元胞大小、曼宁糙率系数、坐标、邻域关系、模型时间步长以及模拟时长。元胞边界设置如图2所示,本发明中元胞的邻域关系为Moore型八邻域,边界元胞设置为极大值从而保证水量无法进入。
(1.2)获取研究区域随时间变化的降雨量,可以是真实情况下的降雨过程,也可以是根据防洪标准假设的不同情境下的降雨过程,进而获得每个元胞不同时刻的降雨量。
(2)概化城市地下排水管网,根据研究区排水管网建设情况,获得排水管网的出口位置、管道长度、管径、管道数量、节点数量及其所在位置,并将雨水井口作为地表和管网水流的交换节点,将雨水井口对应的地表元胞作为节点元胞。该步骤具体包括:
(2.1)获取研究区排水管网资料,取得排水管网管道长度、管道内径、管道坡度、管道糙率、节点数量及其位置以及排水出口位置等参数。
(2.2)将雨水井口所处的元胞作为地表和管网水流交换的节点元胞,通过雨水井口耦合地表径流模型和管网汇流模型进行地表和管网水量交换的计算。管网与地面在节点处耦合示意图如图3所示。
(3)地表径流与管网汇流耦合计算,进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件,进行管网汇流计算获得当前时间步长管道流量和节点水头,根据降雨过程中的地表入流、节点溢流以及退水回流三种不同的过程确定当前时刻节点净流量,联立管道流量方程和节点控制方程求解下一时刻的管道流量、节点水头以及节点净流量,重复计算直至降雨径流过程结束。该步骤具体包括:
(3.1)进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件,具体为:
(3.1.1)利用初损后损法进行产流计算获得每个元胞水深:初损后损法将下渗简化为初损和后损两个部分,降雨开始到出现产流时的水量损失为初损量I0;产流后的水量损失为后损,用后损平均入渗率f表示,mm/min;每个元胞的净雨深以及净雨量计算公式如下:
h=X-I0-ftc-Xn
Vol=0.001×h×A
式中,X为总降雨量,mm;h为X形成的每个元胞的水深,mm;tc为超渗产流历时,min;Xn为后损阶段非超渗历时的降雨量,mm;A为每个元胞的面积,m2;Vol为每个元胞的水量,m3。
(3.1.2)判断水流流向:比较中心元胞与邻居元胞的水位,排除高于平均水位的邻居元胞,确定可以进行水量转移的邻居元胞。
(3.1.3)进行地表汇流计算:在产流和流向确定的基础上,采用曼宁公式计算一个时间步长内中心元胞可以转移的水量,同时根据元胞间水位差得出中心元胞可以向邻居元胞转移的潜在最大水量,将二者相比,得出一个时间步长内中心元胞转移的真实水量;计算公式如下:
式中,Ti为中心元胞水深转移至第i个邻居元胞所需时间,min;h为中心元胞水深,mm;si为中心元胞与第i个邻居元胞间的水面坡度;n为曼宁糙率系数,s/m1/3;Δt为模型时间步长,s;D为中心元胞与邻居元胞的距离,m;表示中心元胞向第i个邻居元胞转移的水量,m3;hi′表示转移后第i个邻居元胞的水深,mm。
(3.1.4)计算地表流入雨水井口的流量:市道路排水***中,边沟流量一般不大于0.4m3/s,纵坡不大于0.08,因此边沟水流流入雨水口的流态为堰流,边沟即对应于雨水井口的各个邻居元胞,水流由节点元胞邻居元胞流入,水深为h′的元胞水流流入雨水口的流量(q,m3/s)采用侧收缩堰流量计算公式:
式中,ε为收缩系数,采用经验公式计算,m为堰流系数,按排水规范中堰流计算取m=0.375,B′为堰顶过水断面宽度,m;g为重力加速度,m2/s;H0为水流总水头高度,H0=0.001×h′+v2/2g,m;v为邻居元胞水流流速,m/s。
根据地表汇流计算获得t时刻节点元胞第i个邻居元胞水深h′i,t时刻第i个邻居元胞流入雨水井口节点的流量采用侧收缩堰流公式计算得到qi,地表流入雨水井口的总流量则为各个邻居元胞流入雨水井口的流量之和qsum=∑qi。
(3.2)进行管网汇流计算获得管道流量以及节点水头随时间的变化过程,具体为:
(3.2.1)获得管道流量基本方程和排水管网节点控制方程联立求解各个时刻管道水流流量和节点水头:
管道内流量基本方程:
排水管网节点控制方程:
其中,Qpipe为管道流量,m3/s;A为管道过水断面面积,m2;Hpipe为管道静压水头,m;t为时间,s;x为距离,m;k为损失系数,k=gn2;R为水力半径,m;V为管道内水流流速,m/s,使用绝对值可以保证摩阻与水流方向相反;Hnode为节点水头,m;Qnode为节点净流量,m3/s;Ask为节点面积,m2。
将管道流量基本方程和节点控制方程用有限差分形式表示:
其中H1、H2为管道上下游节点水头,m;A1、A2为节点上下游管道断面面积,m2;ΔA为一个时间步长内平均断面面积变化值;L为管道的长度,m;为管道过水断面平均面积,m2;当管道水流全部充满管道时,取为管道最宽处的断面面积;为平均水力半径,m;为管道平均流速,m3/s;为t时刻的节点净流量,m3/s;为t时刻的节点水头,m;为t时刻的管道流量,m3/s。
(3.2.2)方程组求解步骤:根据已知的初始的节点水头节点净流量管道流量和初始管道水头一般情况下t=0时管道流量和水头以及节点流量和水头均为0,也可以根据实际情况赋予初始节点水头和管道流量,但均为已知量,将初始条件代入步骤(3.2.1)中管道流量基本方程和节点控制方程,求解下一时间步长初的各个管道流量以及各个节点水头再结合步骤(3.3)中不同降雨时期节点处的流量过程求得节点净流量重复上述步骤,经过t/Δt个时间步长,即可求得t时刻节点水头和管道流量Qt pipe,重复计算直至模拟结束。
(3.3)节点净流量计算:根据降雨初期水流由地表流入管网的地表入流过程、降雨强度增大后的管网水流溢出地表的节点溢流过程以及降雨后期水流重新进入管道的退水回流三个过程,处于不同过程时节点的净流量计算方式不同,具体为:
地表入流过程:降雨初期,排水管网排水能力大于排水量,节点未溢流,处于地表入流状态,水流由地表流入雨水井口进入地下管网,节点水头低于节点地表高程(图4中的(a)),则t时刻根据步骤(3.1)求得节点地表流入流量根据步骤(3.2)求得的t时刻各个管道流量确定节点上下游管道流量和进而获得节点净流量:
节点溢流过程:当降雨量逐渐增大,雨水井口节点水头逐渐增加,节点溢流,当雨水井口处节点水头高于地表高程时,此时节点元胞地表高程以上水量根据步骤(3.1.1)-(3.1.3)中基于元胞自动机的地表径流计算过程计算节点元胞处的水量转移量,当节点元胞转移出去的水量Volt trans_out小于其邻居元胞转移进来的水量Volt trans_in时,表现为节点进流(图4中的(b)),t时刻节点地表进流流量 再根据步骤(3.2)求得节点上下游管道流量分别为和获得此时节点流量:
当节点元胞转移出去的水量Volt trans_out大于其邻居元胞转移进来的流量Volt trans _in时,表现为节点出流(图4中的(c)),节点出流流量再根据步骤(3.2)求得节点上下游管道流量分别为和获得此时节点净流量:
退水回流过程:当降雨量经过峰值后逐渐减小,排水量逐渐减小,地表积水逐渐进入地下管网,进入退水回流阶段,退水回流过程中,节点水头仍高于地表高程时,按节点溢流过程计算节点净流量,地表积水退去,节点水头低于地表高程时,按地表入流过程计算节点净流量,将节点净流量和节点水头代入管道流量基本方程和排水管网节点控制方程获得下一时刻的节点水头和管道流量,直至降雨径流过程结束。
根据设置的研究区域的地表与管网情况以及降雨过程,得到的降雨径流过程模拟情况如图5所示。
Claims (6)
1.一种基于元胞自动机的城市地表径流与管网汇流耦合方法,其特征在于,包括如下步骤:
(1)构建地表元胞地理场景,根据研究区下垫面类型、DEM高程和随时间变化的降雨量,获得每个元胞的下垫面类型、高程和随时间变化的降雨量,定义模型地理场景的元胞大小、邻域关系以及模型边界;
(2)概化城市地下排水管网,根据研究区排水管网建设情况,获得排水管网管道参数、节点参数以及排水出口位置,并将雨水井口作为地表和管网水流的交换节点;
(3)以雨水井口为耦合节点,进行地表径流与管网汇流耦合计算,包括:进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件;进行管网汇流计算获得当前时间步长管道流量和节点水头,再根据降雨过程中的地表入流、节点溢流以及退水回流三种不同的过程确定当前时刻节点净流量;将当前时刻管道流量、节点水头与节点净流量代入管道流量基本方程和排水管网节点控制方程,获得下一时刻管道流量和节点水头,再确定节点净流量,重复计算,直至达到模拟总时长,结束计算。
2.根据权利要求1所述的基于元胞自动机的城市地表径流与管网汇流耦合方法,其特征在于,所述步骤(3)中进行地表产汇流计算获得雨水井口入流流量随时间的变化过程作为管网汇流计算的边界条件包括:
利用初损后损法进行产流计算获得每个元胞水深;
比较中心元胞与邻居元胞的水位,排除高于平均水位的邻居元胞,确定可以进行水量转移的邻居元胞;
在产流和流向确定的基础上,采用曼宁公式计算一个时间步长内中心元胞可以转移的水量,同时根据元胞间水位差得出中心元胞可以向邻居元胞转移的潜在最大水量,将二者相比,得出一个时间步长内中心元胞转移的真实水量;
采用侧收缩堰流量计算公式计算各邻居元胞流入雨水井口节点的流量,并求和得到地表流入雨水井口节点的总流量。
6.根据权利要求1所述的基于元胞自动机的城市地表径流与管网汇流耦合方法,其特征在于,所述步骤(3)中,处于降雨后期水流重新进入管道的退水回流过程时,节点流量计算方式为:退水回流过程中,节点水头高于地表高程时,按节点溢流过程计算节点净流量;地表积水退去,节点水头低于地表高程时,按地表入流过程计算节点净流量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010860835.XA CN112052545B (zh) | 2020-08-25 | 2020-08-25 | 一种基于元胞自动机的城市地表径流与管网汇流耦合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010860835.XA CN112052545B (zh) | 2020-08-25 | 2020-08-25 | 一种基于元胞自动机的城市地表径流与管网汇流耦合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112052545A true CN112052545A (zh) | 2020-12-08 |
CN112052545B CN112052545B (zh) | 2021-10-08 |
Family
ID=73600157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010860835.XA Active CN112052545B (zh) | 2020-08-25 | 2020-08-25 | 一种基于元胞自动机的城市地表径流与管网汇流耦合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112052545B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112711825A (zh) * | 2021-03-26 | 2021-04-27 | 中国水利水电科学研究院 | 一种地表与管网分布式的直接双向耦合方法 |
CN113627817A (zh) * | 2021-08-19 | 2021-11-09 | 东莞理工学院 | 一种基于城市内涝条件下道路连通性分析方法 |
CN117252132A (zh) * | 2023-11-20 | 2023-12-19 | 中山大学 | 一种流域产汇流模拟方法、装置、设备和存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070288168A1 (en) * | 2006-03-24 | 2007-12-13 | Tristan Salles | Use of the Stationary State of Energy-Exchanging Cellular Automata for Modelling Sedimentary Architectures |
US20150261893A1 (en) * | 2014-04-22 | 2015-09-17 | China University Of Petroleum - Beijing | Method and apparatus for determining pipeline flow status parameter of natural gas pipeline network |
CN108388714A (zh) * | 2018-02-07 | 2018-08-10 | 杭州师范大学 | 流域水系和城市管网耦合的平原河网城市洪水模拟方法 |
CN109255185A (zh) * | 2018-09-11 | 2019-01-22 | 中国水利水电科学研究院 | 一种基于城市地表地下管网的一二维水动力学耦合分析方法 |
CN109284531A (zh) * | 2018-08-03 | 2019-01-29 | 中国水利水电科学研究院 | 一种基于空间拓扑的一二维水动力学耦合方法 |
CN109537537A (zh) * | 2018-11-30 | 2019-03-29 | 河北海王水工机械有限公司 | 一种量测水方法及其装置 |
CN109583035A (zh) * | 2018-11-05 | 2019-04-05 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于元胞自动机的城市面源污染物累积冲刷过程计算方法 |
CN110442988A (zh) * | 2019-08-12 | 2019-11-12 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于元胞自动机的城市地表径流流向计算方法及装置 |
-
2020
- 2020-08-25 CN CN202010860835.XA patent/CN112052545B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070288168A1 (en) * | 2006-03-24 | 2007-12-13 | Tristan Salles | Use of the Stationary State of Energy-Exchanging Cellular Automata for Modelling Sedimentary Architectures |
US20150261893A1 (en) * | 2014-04-22 | 2015-09-17 | China University Of Petroleum - Beijing | Method and apparatus for determining pipeline flow status parameter of natural gas pipeline network |
CN108388714A (zh) * | 2018-02-07 | 2018-08-10 | 杭州师范大学 | 流域水系和城市管网耦合的平原河网城市洪水模拟方法 |
CN109284531A (zh) * | 2018-08-03 | 2019-01-29 | 中国水利水电科学研究院 | 一种基于空间拓扑的一二维水动力学耦合方法 |
CN109255185A (zh) * | 2018-09-11 | 2019-01-22 | 中国水利水电科学研究院 | 一种基于城市地表地下管网的一二维水动力学耦合分析方法 |
CN109583035A (zh) * | 2018-11-05 | 2019-04-05 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于元胞自动机的城市面源污染物累积冲刷过程计算方法 |
CN109537537A (zh) * | 2018-11-30 | 2019-03-29 | 河北海王水工机械有限公司 | 一种量测水方法及其装置 |
CN110442988A (zh) * | 2019-08-12 | 2019-11-12 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于元胞自动机的城市地表径流流向计算方法及装置 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112711825A (zh) * | 2021-03-26 | 2021-04-27 | 中国水利水电科学研究院 | 一种地表与管网分布式的直接双向耦合方法 |
CN112711825B (zh) * | 2021-03-26 | 2021-07-30 | 中国水利水电科学研究院 | 一种地表与管网分布式的直接双向耦合方法 |
CN113627817A (zh) * | 2021-08-19 | 2021-11-09 | 东莞理工学院 | 一种基于城市内涝条件下道路连通性分析方法 |
CN113627817B (zh) * | 2021-08-19 | 2024-05-28 | 东莞理工学院 | 一种基于城市内涝条件下道路连通性分析方法 |
CN117252132A (zh) * | 2023-11-20 | 2023-12-19 | 中山大学 | 一种流域产汇流模拟方法、装置、设备和存储介质 |
CN117252132B (zh) * | 2023-11-20 | 2024-04-02 | 中山大学 | 一种流域产汇流模拟方法、装置、设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112052545B (zh) | 2021-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112052545B (zh) | 一种基于元胞自动机的城市地表径流与管网汇流耦合方法 | |
CN107832931B (zh) | 一种平原水网地区内涝风险的模块化分析方法 | |
CN110298076B (zh) | 一种基于gis和swmm的城市内涝智能建模及分析方法 | |
CN106884405B (zh) | 一种无资料地区溃堤型山洪灾害分析评价方法 | |
CN113610264B (zh) | 一种精细化电网台风洪涝灾害预测*** | |
CN109543275B (zh) | 一种城区地表径流二维数值模拟方法 | |
CN113723024B (zh) | 一种适用于滨海地区的“溪流”-“河道”-“河口”分布式洪水过程模拟方法 | |
CN114647881B (zh) | 顾及建筑物微观水文过程的城市内涝建模方法 | |
CN107239657A (zh) | 一种面向对象的水动力学建模要素管理方法 | |
CN111046551B (zh) | 一种城市群排水过程模拟方法 | |
CN112084671B (zh) | 城市时变增益降雨-径流过程模拟计算方法 | |
CN108269199B (zh) | 一种面向对象的小水库群时空分布式出流计算方法 | |
CN111062125A (zh) | 海绵型综合管廊水文效应评估方法 | |
CN114942596B (zh) | 一种城市防洪排水智能控制*** | |
CN112528563A (zh) | 一种基于svm算法的城市内涝预警方法 | |
CN114936505B (zh) | 一种城市雨水井多点位水深快速预报的方法 | |
CN114676473B (zh) | 一种基于人工智能算法的绿色基础设施空间布局优化方法 | |
CN115510771A (zh) | 一种基于水文水动力耦合模型的社区雨涝数值模拟方法 | |
Han et al. | A coupled 1-D and 2-D channel network mathematical model used for flow calculations in the middle reaches of the Yangtze River | |
CN112926786B (zh) | 一种基于关联规则模型和数值模拟的浅水湖泊目标水位逆向预测方法及*** | |
CN113887053A (zh) | 一种面向管网水流计算的市政排水资料质量评估方法及*** | |
Liaw et al. | ASSESSING FLOOD MITIGATION ALTERNATIVES IN SHIJR AREA IN METROPOLITAN TAIPEI 1 | |
CN111985082A (zh) | 一种高守恒性水文水动力耦合模拟方法 | |
Ghimire et al. | Application of cellular automata approach for fast flood simulation | |
Mendoza-Lugo et al. | Reliability analysis of the ancient Nezahualcoyotl’s dike: Investigating failure due to overflow using an improved hydrological model. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |