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

[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]
三、为洪水风险分析提供技术支持：本发明中的模拟方法避免了复杂地经验化参数选取，主要利用曼宁公式、水量平衡原理，模拟洪水淹没进程，计算洪水过程中研究区域内的动态水深和流速，为研究区域洪水风险分析提供高效的技术手段以及相关的数据支持。
附图说明
[0067]
图1为本发明中模拟方法的流程示意图；
[0068]
图2为本发明中von neumann元胞邻域关系示意图；
[0069]
图3为本发明中第一种极限情况下元胞之间水量转移的示意图；
[0070]
图4为本发明中第二种极限情况下元胞之间水量转移的示意图；
[0071]
图5为本发明中元胞水量流动算法的流程示意图。
具体实施方式
[0072]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都应属于本发明保护的范围。
[0073]
一种基于元胞自动机的复杂地形洪水淹没进程模拟方法，其流程如图1所示，包括：
[0074]
步骤1：获取待模拟区域数据信息，包括：待模拟区域数字高程数据dem、下垫面类型(例如：草地、林地、房屋、道路等)、曼宁糙率系数以及初始状态下待模拟区域的洪水水深数据；
[0075]
元胞自动机洪水模型参数如表1所示。
[0076]
表1元胞自动机洪水模型参数
[0077][0078]
步骤2：基于元胞自动机模型，建立洪水淹没进程模型的基础地理场景；
[0079]
元胞自动机模型具体为：
[0080]
将元胞自动机模型的元胞形状结构设置为正方形，元胞的领域关系采用von neumann型四邻域，本实施例中von neumann领域关系如图2所示；
[0081]
根据获得的数字高程数据dem的精度选取元胞的面积；
[0082]
采用笛卡尔坐标系定义每个元胞的位置；
[0083]
定义边缘细胞的开边界或闭边界属性；
[0084]
步骤3：确定对每一个中心元胞邻域元胞中的下游元胞，具体为：
[0085]
当元胞水深大于等于干湿分界的最小水深ε1，即时，元胞中水量可以流向其邻域元胞；
[0086]
比较中间元胞与其邻域元胞的水位值，当中心元胞的水位高于邻域元胞，并且水位之差大于最小水位差ε2时，即时，则此邻域元胞为中心元胞的下游元胞，水量仅从中心元胞转移到其邻域的下游元胞；
[0087]
步骤4：计算t时刻元胞的流速，具体为：
[0088]
在下一个时间步长δt内，从中心元胞流向其某个下游元胞的流速值为：
[0089][0090]
其中，为t时刻中心元胞i与其下游元胞j之间的水力坡度，计算式为为t时刻中心元胞i与其下游元胞j之间的水力坡度，计算式为为t时刻，中心元胞的水位值与两元胞中心元胞i与其下游元胞j中地表高程较高值的差；
[0091]
步骤5：确定t时刻的下一个时间步长δt内，对于每一个中心元胞i转移水量至其下游元胞的分配权重，具体为：
[0092]
对于每一个中心元胞i，根据t时刻中心元胞i流向下游元胞的流速，确定t时刻的下一个时间步长δt内，其下游元胞的水量分配权重：
[0093][0094]
其中，m表示对于中心元胞i，其邻域有m个下游元胞；
[0095]
步骤6：计算t时刻的下一个时间步长δt内，进行元胞的水量转移计算；
[0096]
进行元胞水量转移时，需要遵循从水位高处转移到水位低处的原则，因此中心元胞在转移水量后，不得低于任一个下游元胞；
[0097]
所述的元胞水量转移包括两种情况，分别为中心元胞转出水量后的水位与下游元胞获得水量后的水位最高值相等以及中心元胞转出其t时刻所有水量后，高程仍然高于下游元胞获得水量后的水位最高值；
[0098]
第一种极限情况：如图3所示，元胞水量转移情况为中心元胞i转出水量后的水位与下游元胞获得水量后的水位最高值相等时，具体计算方法为：
[0099][0100]
即
[0101]
其中，a为元胞的面积；
[0102]
第二种极限情况：如图4所示，元胞水量转移情况为中心元胞i转出其t时刻所有水量后，高程仍然高于下游元胞获得水量后的水位最高值，此时中心元胞i转出的水量为具体计算方法为：
[0103][0104]
t时刻的下一个时间步长δt内，中心元胞i转出水量为根据下式计算：
[0105][0106]
即
[0107]
t时刻的下一个时间步长δt内，中心元胞i的下游元胞获得的水量为根据下式计算：
[0108][0109]
步骤7：确定模型的自适应时间步长δt，具体为：
[0110]
根据下式计算模型的自适应时间步长δt：
[0111]
当即时，那么
[0112]
当时，
[0113]
当时，
[0114]
综上，即
[0115]
即δt取待模拟区域内的最小值；
[0116]
步骤8：计算t+δt时刻元胞的水深和流速，具体为：
[0117]
如图5所示，对于元胞i，其在t时刻的下一个时间步长δt内，交换的水量是其作为中心元胞转移水量和作为其他元胞的下游元胞获得水量的总和，因此元胞i在t+δt时刻的水位为：
[0118][0119]
水深为：
[0120][0121]
t时刻的下一个时间步长δt内，中心元胞的流速为中心元胞流向其下游元胞流速的矢量和，即
[0122][0123]
步骤9：返回步骤4，继续进行迭代，直至满足模拟时长要求；
[0124]
步骤10：完成对洪水淹没进行的模拟。
[0125]
本发明基于元胞自动机，通过其灵活的元胞状态转换规则以及适应于并行计算的方法框架，能够适用于复杂地形下多流向的洪水淹没进程快速动态模拟，为实时及动态洪水风险分析提供高效的技术方法。
[0126]
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。