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

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

【技术实现步骤摘要】
一种基于元胞自动机的复杂地形洪水淹没进程模拟方法


[0001]本专利技术涉及洪水过程模拟和防灾减灾
，尤其是涉及一种基于元胞自动机的复杂地形洪水淹没进程模拟方法。

技术介绍

[0002]洪水灾害是我国发生频率最高、影响范围最广、威胁人类生命财产最严重的自然灾害之一。随着计算机技术、地理信息技术的发展，洪水计算模拟技术也在不断进步。
[0003]目前，作为洪水淹没计算模拟方法，水面高程推算法和水动力模型法是理论和技术较为成熟的常用方法。水面高程推算法虽然能快速计算得到大范围洪水淹没情况，但是不能动态模拟洪水淹没进程。水动力模型法虽然能够得到详细洪水淹没进程，但是通过数值求解二维水动力方程，所需计算时间较长。
[0004]中国专利CN110442988B中公开了一种基于元胞自动机的城市地表径流流向计算方法及装置，该专利利用元胞自动机进行经验化的径流流向与流量计算，在模拟过程中使用了复杂的经验化参数用于模拟洪水过程中的动态水深和流速，经验参数的物理可解释性较差。
[0005]中国专利CN109190168A中公开了一种实时测量数据支持下的洪水演进模拟动态修正方法，该专利采用实时测量数据对洪水演进模拟结果进行修正，依赖于大量的洪水实时数据(如，水位等)，对于目前大部分地区还不具有适用性。

技术实现思路

[0006]本专利技术的目的就是为了克服上述现有技术存在的缺陷而提供一种实现实时模拟以及大量计算模拟需求、精度高、稳定性好的基于元胞自动机的复杂地形洪水淹没进程模拟方法。
[0007]本专利技术的目的可以通过以下技术方案来实现：
[0008]一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，所述的模拟方法包括以下步骤：
[0009]步骤1：获取待模拟区域数据信息；
[0010]步骤2：基于元胞自动机模型，建立洪水淹没进程模型的基础地理场景；
[0011]步骤3：确定对每一个中心元胞邻域元胞中的下游元胞；
[0012]步骤4：计算t时刻元胞的流速；
[0013]步骤5：确定t时刻的下一个时间步长Δt内，对于每一个中心元胞i转移水量至其下游元胞的分配权重；
[0014]步骤6：计算t时刻的下一个时间步长Δt内，进行元胞的水量转移计算；
[0015]步骤7：确定模型的自适应时间步长Δt；
[0016]步骤8：计算t+Δt时刻元胞的水深和流速；
[0017]步骤9：返回步骤4，继续进行迭代，直至满足模拟时长要求；
[0018]步骤10：完成对洪水淹没进行的模拟。
[0019]优选地，所述步骤1中待模拟区域数据信息包括：待模拟区域数字高程数据DEM、下垫面类型、曼宁糙率系数以及初始状态下待模拟区域的洪水水深数据。
[0020]更加优选地，所述步骤2包括：
[0021]将元胞自动机模型的元胞形状结构设置为正方形，元胞的领域关系采用von Neumann型四邻域；
[0022]根据获得的数字高程数据DEM的精度选取元胞的面积；
[0023]采用笛卡尔坐标系定义每个元胞的位置；
[0024]定义边缘细胞的开边界或闭边界属性。
[0025]优选地，所述的步骤3具体为：
[0026]当元胞水深大于等于干湿分界的最小水深ε1，即时，元胞中水量可以流向其邻域元胞；
[0027]比较中间元胞与其邻域元胞的水位值，当中心元胞的水位高于邻域元胞，并且水位之差大于最小水位差ε2时，即时，则此邻域元胞为中心元胞的下游元胞，水量仅从中心元胞转移到其邻域的下游元胞。
[0028]更加优选地，所述的步骤4具体为：
[0029]在下一个时间步长Δt内，从中心元胞流向其某个下游元胞的流速值为：
[0030][0031]其中，为t时刻中心元胞i与其下游元胞j之间的水力坡度，计算式为为t时刻中心元胞i与其下游元胞j之间的水力坡度，计算式为为t时刻，中心元胞的水位值与两元胞中心元胞i与其下游元胞j中地表高程较高值的差。
[0032]更加优选地，所述的步骤5具体为：
[0033]对于每一个中心元胞i，根据t时刻中心元胞i流向下游元胞的流速，确定t时刻的下一个时间步长Δt内，其下游元胞的水量分配权重：
[0034][0035]其中，M表示对于中心元胞i，其邻域有M个下游元胞。
[0036]更加优选地，所述的步骤6中进行元胞水量转移时，需要遵循从水位高处转移到水位低处的原则，因此中心元胞在转移水量后，不得低于任一个下游元胞；
[0037]所述的元胞水量转移包括两种情况，分别为中心元胞转出水量后的水位与下游元胞获得水量后的水位最高值相等以及中心元胞转出其t时刻所有水量后，高程仍然高于下游元胞获得水量后的水位最高值。
[0038]更加优选地，所述的元胞水量转移情况为中心元胞i转出水量后的水位与下游元胞获得水量后的水位最高值相等时，具体计算方法为：
[0039][0040]即
[0041]其中，A为元胞的面积；
[0042]所述的元胞水量转移情况为中心元胞i转出其t时刻所有水量后，高程仍然高于下游元胞获得水量后的水位最高值，此时中心元胞i转出的水量为具体计算方法为：
[0043][0044]t时刻的下一个时间步长Δt内，中心元胞i转出水量为根据下式计算：
[0045][0046]即
[0047]t时刻的下一个时间步长Δt内，中心元胞i的下游元胞获得的水量为根据下式计算：
[0048][0049]更加优选地，所述的步骤7具体为：
[0050]根据下式计算模型的自适应时间步长Δt：
[0051]当即时，那么
[0052]当时，
[0053]当时，
[0054]综上，即
[0055]即Δt取待模拟区域内的最小值。
[0056]更加优选地，所述的步骤8具体为：
[0057]对于元胞i，其在t时刻的下一个时间步长Δt内，交换的水量是其作为中心元胞转移水量和作为其他元胞的下游元胞获得水量的总和，因此元胞i在t+Δt时刻的水位为：
[0058][0059]水深为：
[0060][0061]t时刻的下一个时间步长Δt内，中心元胞的流速为中心元胞流向其下游元胞流速的矢量和，即
[0062][0063]与现有技术相比，本专利技术具有以下有益效果：
[0064]一、精度高，稳定性好：本专利技术中的模拟方法通过多流向的权重分配以及由最大转移水量约束下的自适应时间步长计算式，能够更好地减少洪水淹没进程模拟中水位振荡现象，提高模拟精度及稳定性。
[0065]二、实现实时模拟以及大量计算模拟需求：本专利技术中的模拟方法基于元胞自动机，方法框架适用于并行计算，计算成本与用时相比于水动力模型方法大幅度减小，能够实现洪水淹没进程的实时模拟以及洪水风险中大量计算模拟的需求。
[0066]三、为洪水风险分析提供技术支持：本专利技术中的模拟方本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其特征在于，所述的模拟方法包括以下步骤：步骤1：获取待模拟区域数据信息；步骤2：基于元胞自动机模型，建立洪水淹没进程模型的基础地理场景；步骤3：确定对每一个中心元胞邻域元胞中的下游元胞；步骤4：计算t时刻元胞的流速；步骤5：确定t时刻的下一个时间步长Δt内，对于每一个中心元胞i转移水量至其下游元胞的分配权重；步骤6：计算t时刻的下一个时间步长Δt内，进行元胞的水量转移计算；步骤7：确定模型的自适应时间步长Δt；步骤8：计算t+Δt时刻元胞的水深和流速；步骤9：返回步骤4，继续进行迭代，直至满足模拟时长要求；步骤10：完成对洪水淹没进行的模拟。2.根据权利要求1所述的一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其特征在于，所述步骤1中待模拟区域数据信息包括：待模拟区域数字高程数据DEM、下垫面类型、曼宁糙率系数以及初始状态下待模拟区域的洪水水深数据。3.根据权利要求2所述的一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其特征在于，所述步骤2包括：将元胞自动机模型的元胞形状结构设置为正方形，元胞的领域关系采用von Neumann型四邻域；根据获得的数字高程数据DEM的精度选取元胞的面积；采用笛卡尔坐标系定义每个元胞的位置；定义边缘细胞的开边界或闭边界属性。4.根据权利要求1所述的一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其特征在于，所述的步骤3具体为：当元胞水深大于等于干湿分界的最小水深ε1，即时，元胞中水量可以流向其邻域元胞；比较中间元胞与其邻域元胞的水位值，当中心元胞的水位高于邻域元胞，并且水位之差大于最小水位差ε2时，即时，则此邻域元胞为中心元胞的下游元胞，水量仅从中心元胞转移到其邻域的下游元胞。5.根据权利要求4所述的一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其特征在于，所述的步骤4具体为：在下一个时间步长Δt内，从中心元胞流向其某个下游元胞的流速值为：其中，为t时刻中心元胞i与其下游元胞j之间的水力坡度，计算式为为t时刻中心元胞i与其下游元胞j之间的水力坡度，计算式为为t时刻，中心元胞的水位值与两元胞中心元胞i与其下游元胞j中地表高程较高值...

【专利技术属性】
技术研发人员：刘曙光，顾家悦，周正正，钟桂辉，甄亿位，顾圣华，鞠玉梅，林荷娟，姜桂花，陈甜，
申请(专利权)人：同济大学，
类型：发明
国别省市：