CN112270115B - 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法 - Google Patents

一种基于元胞自动机的复杂地形洪水淹没进程模拟方法 Download PDF

Info

Publication number
CN112270115B
CN112270115B CN202011336996.5A CN202011336996A CN112270115B CN 112270115 B CN112270115 B CN 112270115B CN 202011336996 A CN202011336996 A CN 202011336996A CN 112270115 B CN112270115 B CN 112270115B
Authority
CN
China
Prior art keywords
cell
water
time
downstream
cells
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.)
Active
Application number
CN202011336996.5A
Other languages
English (en)
Other versions
CN112270115A (zh
Inventor
刘曙光
顾家悦
周正正
钟桂辉
甄亿位
顾圣华
鞠玉梅
林荷娟
姜桂花
陈甜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN202011336996.5A priority Critical patent/CN112270115B/zh
Publication of CN112270115A publication Critical patent/CN112270115A/zh
Application granted granted Critical
Publication of CN112270115B publication Critical patent/CN112270115B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，包括：获取待模拟区域数据信息；基于元胞自动机模型，建立洪水淹没进程模型的基础地理场景；确定对每一个中心元胞邻域元胞中的下游元胞；计算t时刻元胞的流速；确定t时刻的下一个时间步长Δt内，对于每一个中心元胞i转移水量至其下游元胞的分配权重；计算t时刻的下一个时间步长Δt内，进行元胞的水量转移计算；确定模型的自适应时间步长Δt；计算t+Δt时刻元胞的水深和流速；返回步骤4，继续进行迭代，直至满足模拟时长要求；完成对洪水淹没进行的模拟。与现有技术相比，本发明具有实现实时模拟以及大量计算模拟需求、精度高、稳定性好等优点。

Description

一种基于元胞自动机的复杂地形洪水淹没进程模拟方法
技术领域
本发明涉及洪水过程模拟和防灾减灾技术领域，尤其是涉及一种基于元胞自动机的复杂地形洪水淹没进程模拟方法。
背景技术
洪水灾害是我国发生频率最高、影响范围最广、威胁人类生命财产最严重的自然灾害之一。随着计算机技术、地理信息技术的发展，洪水计算模拟技术也在不断进步。
目前，作为洪水淹没计算模拟方法，水面高程推算法和水动力模型法是理论和技术较为成熟的常用方法。水面高程推算法虽然能快速计算得到大范围洪水淹没情况，但是不能动态模拟洪水淹没进程。水动力模型法虽然能够得到详细洪水淹没进程，但是通过数值求解二维水动力方程，所需计算时间较长。
中国专利CN110442988B中公开了一种基于元胞自动机的城市地表径流流向计算方法及装置，该专利利用元胞自动机进行经验化的径流流向与流量计算，在模拟过程中使用了复杂的经验化参数用于模拟洪水过程中的动态水深和流速，经验参数的物理可解释性较差。
中国专利CN109190168A中公开了一种实时测量数据支持下的洪水演进模拟动态修正方法，该专利采用实时测量数据对洪水演进模拟结果进行修正，依赖于大量的洪水实时数据(如，水位等)，对于目前大部分地区还不具有适用性。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种实现实时模拟以及大量计算模拟需求、精度高、稳定性好的基于元胞自动机的复杂地形洪水淹没进程模拟方法。
本发明的目的可以通过以下技术方案来实现：
一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，所述的模拟方法包括以下步骤：
步骤1：获取待模拟区域数据信息；
步骤2：基于元胞自动机模型，建立洪水淹没进程模型的基础地理场景；
步骤3：确定对每一个中心元胞邻域元胞中的下游元胞；
步骤4：计算t时刻元胞的流速；
步骤5：确定t时刻的下一个时间步长Δt内，对于每一个中心元胞i转移水量至其下游元胞的分配权重；
步骤6：计算t时刻的下一个时间步长Δt内，进行元胞的水量转移计算；
步骤7：确定模型的自适应时间步长Δt；
步骤8：计算t+Δt时刻元胞的水深和流速；
步骤9：返回步骤4，继续进行迭代，直至满足模拟时长要求；
步骤10：完成对洪水淹没进行的模拟。
优选地，所述步骤1中待模拟区域数据信息包括：待模拟区域数字高程数据DEM、下垫面类型、曼宁糙率系数以及初始状态下待模拟区域的洪水水深数据。
更加优选地，所述步骤2包括：
将元胞自动机模型的元胞形状结构设置为正方形，元胞的领域关系采用vonNeumann型四邻域；
根据获得的数字高程数据DEM的精度选取元胞的面积；
采用笛卡尔坐标系定义每个元胞的位置；
定义边缘细胞的开边界或闭边界属性。
优选地，所述的步骤3具体为：
当元胞水深
Figure BDA0002797530560000022
大于等于干湿分界的最小水深ε1，即
Figure BDA0002797530560000023
时，元胞中水量可以流向其邻域元胞；
比较中间元胞与其邻域元胞的水位值，当中心元胞的水位高于邻域元胞，并且水位之差大于最小水位差ε2时，即
Figure BDA0002797530560000021
时，则此邻域元胞为中心元胞的下游元胞，水量仅从中心元胞转移到其邻域的下游元胞。
更加优选地，所述的步骤4具体为：
在下一个时间步长Δt内，从中心元胞流向其某个下游元胞的流速值为：
Figure BDA0002797530560000031
其中，
Figure BDA0002797530560000032
为t时刻中心元胞i与其下游元胞j之间的水力坡度，计算式为
Figure BDA0002797530560000033
Figure BDA0002797530560000034
为t时刻，中心元胞的水位值与两元胞中心元胞i与其下游元胞j中地表高程较高值的差。
更加优选地，所述的步骤5具体为：
对于每一个中心元胞i，根据t时刻中心元胞i流向下游元胞的流速，确定t时刻的下一个时间步长Δt内，其下游元胞的水量分配权重：
Figure BDA0002797530560000035
其中，M表示对于中心元胞i，其邻域有M个下游元胞。
更加优选地，所述的步骤6中进行元胞水量转移时，需要遵循从水位高处转移到水位低处的原则，因此中心元胞在转移水量后，不得低于任一个下游元胞；
所述的元胞水量转移包括两种情况，分别为中心元胞转出水量后的水位与下游元胞获得水量后的水位最高值相等以及中心元胞转出其t时刻所有水量后，高程仍然高于下游元胞获得水量后的水位最高值。
更加优选地，所述的元胞水量转移情况为中心元胞i转出水量
Figure BDA0002797530560000036
后的水位与下游元胞获得水量后的水位最高值相等时，具体计算方法为：
Figure BDA0002797530560000037
即
Figure BDA00027975305600000313
其中，A为元胞的面积；
所述的元胞水量转移情况为中心元胞i转出其t时刻所有水量后，高程仍然高于下游元胞获得水量后的水位最高值，此时中心元胞i转出的水量为
Figure BDA0002797530560000038
具体计算方法为：
Figure BDA0002797530560000039
t时刻的下一个时间步长Δt内，中心元胞i转出水量为
Figure BDA00027975305600000310
根据下式计算：
Figure BDA00027975305600000311
即
Figure BDA00027975305600000312
t时刻的下一个时间步长Δt内，中心元胞i的下游元胞获得的水量为
Figure BDA0002797530560000041
根据下式计算：
Figure BDA0002797530560000042
更加优选地，所述的步骤7具体为：
根据下式计算模型的自适应时间步长Δt：
当
Figure BDA0002797530560000043
即
Figure BDA0002797530560000044
时，那么
Figure BDA0002797530560000045
当
Figure BDA0002797530560000046
时，
Figure BDA0002797530560000047
当
Figure BDA0002797530560000048
时，
Figure BDA0002797530560000049
综上，
Figure BDA00027975305600000410
即
Figure BDA00027975305600000411
即Δt取待模拟区域内的最小值。
更加优选地，所述的步骤8具体为：
对于元胞i，其在t时刻的下一个时间步长Δt内，交换的水量是其作为中心元胞转移水量和作为其他元胞的下游元胞获得水量的总和，因此元胞i在t+Δt时刻的水位为：
Figure BDA00027975305600000412
水深为：
Figure BDA00027975305600000413
t时刻的下一个时间步长Δt内，中心元胞的流速
Figure BDA00027975305600000414
为中心元胞流向其下游元胞流速的矢量和，即
Figure BDA00027975305600000415
与现有技术相比，本发明具有以下有益效果：
一、精度高，稳定性好：本发明中的模拟方法通过多流向的权重分配以及由最大转移水量约束下的自适应时间步长计算式，能够更好地减少洪水淹没进程模拟中水位振荡现象，提高模拟精度及稳定性。
二、实现实时模拟以及大量计算模拟需求：本发明中的模拟方法基于元胞自动机，方法框架适用于并行计算，计算成本与用时相比于水动力模型方法大幅度减小，能够实现洪水淹没进程的实时模拟以及洪水风险中大量计算模拟的需求。
三、为洪水风险分析提供技术支持：本发明中的模拟方法避免了复杂地经验化参数选取，主要利用曼宁公式、水量平衡原理，模拟洪水淹没进程，计算洪水过程中研究区域内的动态水深和流速，为研究区域洪水风险分析提供高效的技术手段以及相关的数据支持。
附图说明
图1为本发明中模拟方法的流程示意图；
图2为本发明中von Neumann元胞邻域关系示意图；
图3为本发明中第一种极限情况下元胞之间水量转移的示意图；
图4为本发明中第二种极限情况下元胞之间水量转移的示意图；
图5为本发明中元胞水量流动算法的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都应属于本发明保护的范围。
一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其流程如图1所示，包括：
步骤1：获取待模拟区域数据信息，包括：待模拟区域数字高程数据DEM、下垫面类型(例如：草地、林地、房屋、道路等)、曼宁糙率系数以及初始状态下待模拟区域的洪水水深数据；
元胞自动机洪水模型参数如表1所示。
表1元胞自动机洪水模型参数
Figure BDA0002797530560000061
步骤2：基于元胞自动机模型，建立洪水淹没进程模型的基础地理场景；
元胞自动机模型具体为：
将元胞自动机模型的元胞形状结构设置为正方形，元胞的领域关系采用vonNeumann型四邻域，本实施例中von Neumann领域关系如图2所示；
根据获得的数字高程数据DEM的精度选取元胞的面积；
采用笛卡尔坐标系定义每个元胞的位置；
定义边缘细胞的开边界或闭边界属性；
步骤3：确定对每一个中心元胞邻域元胞中的下游元胞，具体为：
当元胞水深
Figure BDA0002797530560000062
大于等于干湿分界的最小水深ε1，即
Figure BDA0002797530560000063
时，元胞中水量可以流向其邻域元胞；
比较中间元胞与其邻域元胞的水位值，当中心元胞的水位高于邻域元胞，并且水位之差大于最小水位差ε2时，即
Figure BDA0002797530560000064
时，则此邻域元胞为中心元胞的下游元胞，水量仅从中心元胞转移到其邻域的下游元胞；
步骤4：计算t时刻元胞的流速，具体为：
在下一个时间步长Δt内，从中心元胞流向其某个下游元胞的流速值为：
Figure BDA0002797530560000071
其中，
Figure BDA0002797530560000079
为t时刻中心元胞i与其下游元胞j之间的水力坡度，计算式为
Figure BDA0002797530560000072
Figure BDA0002797530560000073
为t时刻，中心元胞的水位值与两元胞中心元胞i与其下游元胞j中地表高程较高值的差；
步骤5：确定t时刻的下一个时间步长Δt内，对于每一个中心元胞i转移水量至其下游元胞的分配权重，具体为：
对于每一个中心元胞i，根据t时刻中心元胞i流向下游元胞的流速，确定t时刻的下一个时间步长Δt内，其下游元胞的水量分配权重：
Figure BDA0002797530560000074
其中，M表示对于中心元胞i，其邻域有M个下游元胞；
步骤6：计算t时刻的下一个时间步长Δt内，进行元胞的水量转移计算；
进行元胞水量转移时，需要遵循从水位高处转移到水位低处的原则，因此中心元胞在转移水量后，不得低于任一个下游元胞；
所述的元胞水量转移包括两种情况，分别为中心元胞转出水量后的水位与下游元胞获得水量后的水位最高值相等以及中心元胞转出其t时刻所有水量后，高程仍然高于下游元胞获得水量后的水位最高值；
第一种极限情况：如图3所示，元胞水量转移情况为中心元胞i转出水量
Figure BDA0002797530560000075
后的水位与下游元胞获得水量后的水位最高值相等时，具体计算方法为：
Figure BDA0002797530560000076
即
Figure BDA0002797530560000077
其中，A为元胞的面积；
第二种极限情况：如图4所示，元胞水量转移情况为中心元胞i转出其t时刻所有水量后，高程仍然高于下游元胞获得水量后的水位最高值，此时中心元胞i转出的水量为
Figure BDA0002797530560000078
具体计算方法为：
Figure BDA0002797530560000081
t时刻的下一个时间步长Δt内，中心元胞i转出水量为
Figure BDA0002797530560000082
根据下式计算：
Figure BDA0002797530560000083
即
Figure BDA0002797530560000084
t时刻的下一个时间步长Δt内，中心元胞i的下游元胞获得的水量为
Figure BDA0002797530560000085
根据下式计算：
Figure BDA0002797530560000086
步骤7：确定模型的自适应时间步长Δt，具体为：
根据下式计算模型的自适应时间步长Δt：
当
Figure BDA0002797530560000087
即
Figure BDA0002797530560000088
时，那么
Figure BDA0002797530560000089
当
Figure BDA00027975305600000810
时，
Figure BDA00027975305600000811
当
Figure BDA00027975305600000812
时，
Figure BDA00027975305600000813
综上，
Figure BDA00027975305600000814
即
Figure BDA00027975305600000815
即Δt取待模拟区域内的最小值；
步骤8：计算t+Δt时刻元胞的水深和流速，具体为：
如图5所示，对于元胞i，其在t时刻的下一个时间步长Δt内，交换的水量是其作为中心元胞转移水量和作为其他元胞的下游元胞获得水量的总和，因此元胞i在t+Δt时刻的水位为：
Figure BDA00027975305600000816
水深为：
Figure BDA0002797530560000091
t时刻的下一个时间步长Δt内，中心元胞的流速
Figure BDA0002797530560000092
为中心元胞流向其下游元胞流速的矢量和，即
Figure BDA0002797530560000093
步骤9：返回步骤4，继续进行迭代，直至满足模拟时长要求；
步骤10：完成对洪水淹没进行的模拟。
本发明基于元胞自动机，通过其灵活的元胞状态转换规则以及适应于并行计算的方法框架，能够适用于复杂地形下多流向的洪水淹没进程快速动态模拟，为实时及动态洪水风险分析提供高效的技术方法。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。

Claims (4)

1.一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其特征在于，所述的模拟方法包括以下步骤：
步骤1：获取待模拟区域数据信息；
步骤2：基于元胞自动机模型，建立洪水淹没进程模型的基础地理场景；
步骤3：确定对每一个中心元胞邻域元胞中的下游元胞；
所述的步骤3具体为：
当元胞水深
Figure FDA0003515845970000011
大于等于干湿分界的最小水深ε1，即
Figure FDA0003515845970000012
时，元胞中水量可以流向其邻域元胞；
比较中心元胞与其邻域元胞的水位值，当中心元胞的水位高于邻域元胞，并且水位之差大于等于最小水位差ε2时，即
Figure FDA0003515845970000013
时，则此邻域元胞为中心元胞的下游元胞，水量仅从中心元胞转移到其邻域的下游元胞；
步骤4：计算t时刻元胞的流速；
步骤5：确定t时刻的下一个时间步长Δt内，对于每一个中心元胞i转移水量至其下游元胞的分配权重；
所述的步骤5具体为：
对于每一个中心元胞i，根据t时刻中心元胞i流向下游元胞的流速，确定t时刻的下一个时间步长Δt内，其下游元胞的水量分配权重：
Figure FDA0003515845970000014
其中，M表示对于中心元胞i，其邻域有M个下游元胞；
Figure FDA0003515845970000015
为t时刻中心元胞的水位值与两元胞中心元胞i与其下游元胞j中地表高程较高值的差；
步骤6：计算t时刻的下一个时间步长Δt内，进行元胞的水量转移计算；
所述的步骤6中进行元胞水量转移时，需要遵循从水位高处转移到水位低处的原则，因此中心元胞在转移水量后，不得低于任一个下游元胞；
所述的元胞水量转移包括两种情况，分别为中心元胞转出水量后的水位与下游元胞获得水量后的水位最高值相等以及中心元胞转出其t时刻所有水量后，高程仍然高于下游元胞获得水量后的水位最高值；
所述的元胞水量转移情况为中心元胞i转出水量Q1 t+Δt后的水位与下游元胞获得水量后的水位最高值相等时，具体计算方法为：
Figure FDA0003515845970000021
即
Figure FDA0003515845970000022
其中，A为元胞的面积；
所述的元胞水量转移情况为中心元胞i转出其t时刻所有水量后，高程仍然高于下游元胞获得水量后的水位最高值，此时中心元胞i转出的水量为Q2 t+Δt，具体计算方法为：
Figure FDA0003515845970000023
其中，
Figure FDA0003515845970000024
为元胞i在t时刻的水深；
t时刻的下一个时间步长Δt内，中心元胞i转出水量为
Figure FDA0003515845970000025
根据下式计算：
Figure FDA0003515845970000026
即
Figure FDA0003515845970000027
t时刻的下一个时间步长Δt内，中心元胞i的下游元胞获得的水量为
Figure FDA0003515845970000028
根据下式计算：
Figure FDA0003515845970000029
步骤7：确定模型的自适应时间步长Δt；
所述的步骤7具体为：
根据下式计算模型的自适应时间步长Δt：
当
Figure FDA00035158459700000210
即
Figure FDA00035158459700000211
时，那么
Figure FDA00035158459700000212
当
Figure FDA00035158459700000213
时，
Figure FDA00035158459700000214
当
Figure FDA00035158459700000215
时，
Figure FDA00035158459700000216
综上，
Figure FDA00035158459700000217
即
Figure FDA0003515845970000031
即Δt取待模拟区域内的最小值，其中，n为曼宁糙率系数；
步骤8：计算t+Δt时刻元胞的水深和流速；
所述的步骤8具体为：
对于元胞i，其在t时刻的下一个时间步长Δt内，交换的水量是其作为中心元胞转移水量和作为其他元胞的下游元胞获得水量的总和，因此元胞i在t+Δt时刻的水位为：
Figure FDA0003515845970000032
水深为：
Figure FDA0003515845970000033
t时刻的下一个时间步长Δt内，中心元胞的流速
Figure FDA0003515845970000034
为中心元胞流向其下游元胞流速的矢量和，即
Figure FDA0003515845970000035
步骤9：返回步骤4，继续进行迭代，直至满足模拟时长要求；
步骤10：完成对洪水淹没进行的模拟。
2.根据权利要求1所述的一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其特征在于，所述步骤1中待模拟区域数据信息包括：待模拟区域数字高程数据DEM、下垫面类型、曼宁糙率系数以及初始状态下待模拟区域的洪水水深数据。
3.根据权利要求2所述的一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其特征在于，所述步骤2包括：
将元胞自动机模型的元胞形状结构设置为正方形，元胞的领域关系采用vonNeumann型四邻域；
根据获得的数字高程数据DEM的精度选取元胞的面积；
采用笛卡尔坐标系定义每个元胞的位置；
定义边缘细胞的开边界或闭边界属性。
4.根据权利要求1所述的一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其特征在于，所述的步骤4具体为：
在下一个时间步长Δt内，从中心元胞流向其某个下游元胞的流速值为：
Figure FDA0003515845970000041
其中，
Figure FDA0003515845970000042
为t时刻中心元胞i与其下游元胞j之间的水力坡度，计算式为
Figure FDA0003515845970000043
d为元胞边长。
CN202011336996.5A 2020-11-25 2020-11-25 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法 Active CN112270115B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011336996.5A CN112270115B (zh) 2020-11-25 2020-11-25 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011336996.5A CN112270115B (zh) 2020-11-25 2020-11-25 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法

Publications (2)

Publication Number Publication Date
CN112270115A CN112270115A (zh) 2021-01-26
CN112270115B true CN112270115B (zh) 2022-06-14

Family

ID=74340320

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011336996.5A Active CN112270115B (zh) 2020-11-25 2020-11-25 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法

Country Status (1)

Country Link
CN (1) CN112270115B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113916746B (zh) * 2021-10-14 2023-08-04 深圳大学 一种基于元胞自动机的非饱和地层入渗分析方法
CN113743836A (zh) * 2021-11-08 2021-12-03 国网江西省电力有限公司电力科学研究院 一种基于实时数据的洪涝灾害风险分布图绘制方法及装置
CN114611188B (zh) * 2022-03-09 2022-09-02 生态环境部卫星环境应用中心 基于元胞自动机的尾矿库溃坝泄漏模拟分析方法及系统
CN114970388B (zh) * 2022-05-12 2023-07-25 河海大学 一种自适应时间步长的场地地表污染扩散方法
CN115081740B (zh) * 2022-07-20 2023-07-28 中山大学 基于元胞自动机的河口地区盐度预测方法、系统及设备

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103544390A (zh) * 2013-10-18 2014-01-29 西南交通大学 基于元胞自动机的溃决洪水快速演进模拟分析方法
CN109145396A (zh) * 2018-07-27 2019-01-04 同济大学 一种基于植被分布的河道糙率分区率定方法
CN109190168A (zh) * 2018-07-31 2019-01-11 中国科学院遥感与数字地球研究所 一种实时测量数据支持下的洪水演进模拟动态修正方法
CN109556823A (zh) * 2018-10-30 2019-04-02 同济大学 一种三维含植物河道水动力及水质模型装置
CN109583628A (zh) * 2018-11-05 2019-04-05 水利部交通运输部国家能源局南京水利科学研究院 一种基于元胞自动机的洪涝人员动态疏散路径分析方法
CN109583035A (zh) * 2018-11-05 2019-04-05 水利部交通运输部国家能源局南京水利科学研究院 基于元胞自动机的城市面源污染物累积冲刷过程计算方法
CN110442988A (zh) * 2019-08-12 2019-11-12 水利部交通运输部国家能源局南京水利科学研究院 一种基于元胞自动机的城市地表径流流向计算方法及装置
CN111310322A (zh) * 2020-02-10 2020-06-19 河海大学 基于元胞自动机的干旱区自然绿洲空间动态模拟方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10614255B2 (en) * 2010-10-29 2020-04-07 Bentley Systems, Incorporated Computer-implemented land planning system and method with GIS integration

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103544390A (zh) * 2013-10-18 2014-01-29 西南交通大学 基于元胞自动机的溃决洪水快速演进模拟分析方法
CN109145396A (zh) * 2018-07-27 2019-01-04 同济大学 一种基于植被分布的河道糙率分区率定方法
CN109190168A (zh) * 2018-07-31 2019-01-11 中国科学院遥感与数字地球研究所 一种实时测量数据支持下的洪水演进模拟动态修正方法
CN109556823A (zh) * 2018-10-30 2019-04-02 同济大学 一种三维含植物河道水动力及水质模型装置
CN109583628A (zh) * 2018-11-05 2019-04-05 水利部交通运输部国家能源局南京水利科学研究院 一种基于元胞自动机的洪涝人员动态疏散路径分析方法
CN109583035A (zh) * 2018-11-05 2019-04-05 水利部交通运输部国家能源局南京水利科学研究院 基于元胞自动机的城市面源污染物累积冲刷过程计算方法
CN110442988A (zh) * 2019-08-12 2019-11-12 水利部交通运输部国家能源局南京水利科学研究院 一种基于元胞自动机的城市地表径流流向计算方法及装置
CN111310322A (zh) * 2020-02-10 2020-06-19 河海大学 基于元胞自动机的干旱区自然绿洲空间动态模拟方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
A weighted cellular automata 2D inundation model for rapid flood analysis;Michele Guidolin;《Environmental Modelling & Software》;20160722;378-394 *
An approach for flood monitoring by the combined use of Landsat 8 optical imagery and COSMO-SkyMed radar imagery;Xiaohua Tong;《ISPRS Journal of Photogrammetry and Remote Sensing》;20180111;144-153 *
Statistical Analysis of Temporal Characteristics of Natural Runoff in Golmud River Basin;Yiwei Zhen;《IOP Conference Series: Materials Science and Engineering》;20200531;1-7 *
Storm Catalog‐Based Analysis of Rainfall Heterogeneity and Frequency in a Complex Terrain;Zhengzheng Zhou;《Water Resources Research》;20190304;1871-1889 *
Velocity and turbulence affected by submerged rigid vegetation under waves, currents and combined wave–current flows;Shuguang Liu;《Coastal Engineering》;20200516;1-13 *
城市化进程中的防洪排涝体系建设;刘曙光;《科学》;20200925;32-37 *
暴雨致涝及洪水风险分析虚拟仿真实验的设计与实现;钟桂辉;《实验技术与管理》;20190819;160-163,214 *

Also Published As

Publication number Publication date
CN112270115A (zh) 2021-01-26

Similar Documents

Publication Publication Date Title
CN112270115B (zh) 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法
Balzano Evaluation of methods for numerical simulation of wetting and drying in shallow water flow models
Begnudelli et al. Conservative wetting and drying methodology for quadrilateral grid finite-volume models
CN109918821B (zh) 一种迎风守恒型河道漫溢出流数值模拟方法
Feng et al. Development of an unstructured-grid wave-current coupled model and its application
CN109582996A (zh) 一种小尺度岸滩剖面与大尺度岸线变化的耦合模拟方法
CN114329950B (zh) 基于动态概化的斜坡式潜堤波浪水动力影响数值模拟方法
Wei et al. Simulation of nearshore wave processes by a depth-integrated non-hydrostatic finite element model
CN114580283B (zh) 山溪性强潮河口分汊段中长期动力地貌演变数值模拟方法
Sebastian et al. Storm surge hydrodynamics at a concave coast due to varying approach angles of cyclone
CN113807008A (zh) 一种基于深度学习的城市暴雨内涝模拟方法
CN104750972A (zh) 复杂水网地区楔形水体容积的计算方法
CN112257313A (zh) 一种基于gpu加速的污染物输移高分辨率数值模拟方法
CN107038308B (zh) 一种基于线性内插的规则格网地形建模方法
CN103774605B (zh) 一种提高环抱式港池水体交换能力的设计方法
CN108196951B (zh) Gpu流域径流模拟分布式调度系统及方法
CN110990926B (zh) 一种基于面积修正率的城市地表建筑水动力学仿真方法
CN105956327A (zh) 地面灌溉地表水流运动过程模拟方法
CN110765694B (zh) 一种基于简化型浅水方程组的城市地表水流数值模拟方法
Walters Design considerations for a finite element coastal ocean model
Chen et al. A 3-D numerical study of solitary wave interaction with vertical cylinders using a parallelised particle-in-cell solver
CN110543679B (zh) 一种多波况大尺寸固定式海工平台水平力计算方法
CN110555236B (zh) 一种大尺寸固定式海工平台波浪垂向力计算方法
CN117313325A (zh) 一种用于海洋环流模式的σ-Z-σ混合坐标系构建方法及系统
CN107748926B (zh) 一种南海深海盆中尺度涡预报方法

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