[发明专利]基于元胞自动机的溃决洪水快速演进模拟分析方法

摘要

本发明属于地理信息系统的虚拟地理环境研究领域，特别涉及到溃决洪水演进模拟分析技术。本发明提供一种基于元胞自动机的溃决洪水快速演进模拟分析方法，该方法能够解决目前由元胞尺度划分的随意性而导致的最终数值分析结果与实际水情之间误差较大的情况，同时，合理的元胞邻域范围的设置，结合GIS的空间分析功能和元胞自动机模拟的时空建模能力，能对溃坝洪水演进及影响进行快速准确的综合模拟与分析，包括构建二维元胞自动机模拟模型，对元胞空间数据结构进行优化，基于虚拟地理环境的可视化模拟分析。

主权项

基于元胞自动机的溃决洪水快速演进模拟分析方法，其特征在于，包括如下步骤：S1、构建二维元胞自动机模拟模型，包括：确定元胞状态；建立元胞的邻域关系；建立元胞状态的转换规则；所述元胞状态包括元胞所在处的河床高度，元胞所在处的河床糙率，元胞所在处的水深，元胞所在处的x方向的单宽流量和元胞所在处的方向的单宽流量，根据所述元胞状态，对一个确定的元胞空间可以设置10种元胞状态属性，即高程、t时刻的水深、t+1时刻的水深、x方向t时刻的单宽流量、y方向t时刻的单宽流量、x方向t+1时刻的单宽流量、y方向t+1时刻的单宽流量、x方向的流速、y方向的流速、糙率；所述元胞邻域关系的建立为：建立冯·诺依曼型元胞邻域管理，即中心元胞有上、下、左、右4个邻元单元；所述建立元胞状态的转换规则为：在保证模型精度和计算结果合理为前提下对非线性微分方程进行简化，满足符合圣维南方程求解原理的条件，对所述简化后的非线性微分方程采用离散计算方式来推导元胞状态转换规则；所述元胞状态的转换规则具体如下：S11、设除源头地区的其他地区为陆地地区，所述陆地地区元胞状态初始值中x，y方向单宽流量和水面深度均设为0，取源头地区洪水量和源头地区水面高度量测值作为初始状态数据，S12、根据流速和元胞空间格网大小确定元胞状态变化的时间步长，令表示为第i行j列的元胞在t时刻的状态，S13、所述元胞状态转换规则在满足圣维南方程求解原理的条件下进行二维转换，即，中心元胞t+1时刻的水深和单宽流量由中心元胞相邻的4个元胞t时刻的水深和单宽流量决定，所述相邻的4个元胞分别为中心元胞的上、下、左、右邻域的元胞，利用t时刻的元胞状态，计算t+1时刻元胞的水深和单宽流量，具体方法如下：由t时刻的相邻的元胞的水深计算t+1时刻的中心元胞的单宽流量，计算方法如下：Mi,jk+1=Mi,jk-gΔt(hi+1,jk+hi,jk)(zi+1,jk-zi,jk)Δx-gn2i,ju‾i,jΔt(ui,jk)2+(vi,jk)2[(hi+1,jk+hi,jk)/2]1/3Ni,jk+1=Ni,jk-gΔt(hi,j+1k+hi,jk)(zi,j+1k-zi,jk)Δy-gn2i,jv‾i,jΔt(ui,jk)2+(vi,jk)2[(hi,j+1k+hi,jk)/2]1/3,]]>由t时刻的相邻的元胞的单宽流量计算t+1时刻的中心元胞的水深，计算方法如下：hi,jk+1=ki,jk-Δt(Mi+1,jk+1-Mi,jk+1)Δx-Δt(Ni,j+1k+1-Ni,jk+1)Δy,]]>S14、当元胞所处区域水深＞0时，设该元胞所在区域为河流区域，其中，元胞尺度的确定根据不同比例尺的地形图，遥感影像，实地考察资料，历史水温资料对研究区域进行解译，提取数据整合得到溃决演进元胞自动机数值模型单一元胞的不同空间尺寸，其中，元胞的演化由洪水演进的计算方法和已有的观测数据，用源头即模拟河段的最上游的流速和水位变化决定，其中，M为x方向的单宽流量、N为y方向的单宽流量，i为元胞所在的行、j为元胞所在的列，k为元胞模拟变化的次数，h为柱体平均水深，u和v分别为x、y方向的平均流速,z为水面高程，t为时间，n为糙率，g为重力加速度；S2、对元胞空间数据结构进行优化，包括：S21、元胞尺度效应分析，具体包括如下步骤：S211、收集不同比例尺地形图数据，生成对应空间分辨率的栅格数据，依次对应为溃决洪水元胞自动机模型中的单一尺度元胞的空间尺寸；S212、根据遥感影像、实地考察、历史水文资料，提取对应空间分辨率数据的河长、河流各高程面的面积值、各河段的河长值以及初始河段上端面各高程处的河宽值；S213、根据S212所述的数据计算各河段不同高程处的河宽，以及各河段不同高程处的地形糙率，利用计算结果分别进行模拟实验，记录不同元胞尺寸对溃决洪水演进模型运算效果的影响，分析不同元胞尺度下的效应，比较分析找到适合的溃决模拟模型的元胞尺度和邻域范围并设置参数；S22、基本流域骨干线提取，具体步骤如下：S221、根据S21所得尺度效应分析结果，针对地形的起伏变化状态分别计算最优元胞尺度和邻域范围设置参数，进而提取基本流域骨干线；S222、设定流域范围内细节层次变化要素，进行元胞空间数据结构多细节层次优化，为了简化元胞状态传递变化计算，设定相邻间元胞细节层次变化不超过1层，从而获得要设定元胞尺度变化范围的阈值，以及对应的地形糙率，溃决元胞状态主要包括水深、速度分量状态值，元胞自动机模型运行的关键在于不同时刻相邻元胞状态值的传递；S23、元胞分裂与合并，具体步骤如下：在元胞分裂的情况下，每个元胞t+1时刻的水深由t时刻的水面高减去分裂后元胞高程得到，所述元胞t+1时刻的x方向单宽流量u和y方向单宽流量v的当前时刻分量由分裂元胞的坡度和坡向与分裂前的夹角来计算，而在元胞合并的情况下，元胞t+1时刻的水深由合并前元胞的水面高均值减去合并后的元胞高程得到，所述元胞t+1时刻的x方向单宽流量u和y方向单宽流量v由合并前元胞的坡度坡向值与合并后的夹角关系进行矢量求积得到；S3、基于虚拟地理环境的可视化模拟分析，包括：S31、基于GIS获取用于溃坝数值模型计算的基础地理信息，实现数字高程模型和遥感影像在全球尺度下的分层分块，纹理映射和叠加，进而生成具有真实感的三维地形场景模型，实现对洪水事件对象模型和洪水基地对象模型的一体化显示；S32、根据S31所述的三维地形场景模型，根据溃坝洪水演进模型计算结果，进行动态水流模型的构建，并集成到虚拟地理环境中，通过设计基于视点简化的虚拟场景“实时”绘制算法，可在普通PC机上较流畅实现溃决洪水演进模拟的可视化表达与分析。