洞庭湖区典型蓄洪垸分洪洪水演进模拟与可视化——以澧南垸为例教学 讲义.pptVIP

澧南垸分洪洪水三维淹没程序介绍 程序包括程序操作的基本功能、图像数据的处理部分、分洪洪水给定水位的淹没模拟、分洪洪水淹没分析四个部分，程序界面如下图所示。 澧南垸分洪洪水三维淹没程序的应用界面 基本功能 地图的加载功能。加载已创建的澧南垸图层。 地图的视图浏览、图层控制功能。 洞庭湖区典型蓄洪垸分洪洪水演进模拟与可视化研究 ——以澧南垸为例 毛德华 湖南师范大学资源与环境科学学院 报告提纲 研究区域概况 澧南垸分洪洪水二维演进模拟与可视化 澧南垸分洪洪水三维淹没程序开发与可视化 成果的应用 一、 研究区域概况 1 洞庭湖区概况 二、澧南垸分洪洪水二维演进模拟与可视化 本项目采用的是DHI水力研究所开发的平面二维数学模型 MIKE21 水动力学(HD)模型，以澧南垸分洪闸流量过程为边界条件，采用矩形结构网格剖分计算，模拟澧南垸分洪洪水演进过程。 模型特征 HD 模块是 MIKE21 软件中的基本模块，为泥沙传输和环境水文学提供了水动力学的计算基础。HD 模块利用 ADI 二阶精度的有限差分法对动态流的连续方程和动量守衡方程求解，采用交错矩形网格离散。模拟湖泊、河口和海岸地区的水位变化、和由于各种力的作用而产生的水流变化。为模型提供了地形、底部糙率、风场和水动力学边界条件等输入数据后，可以计算出相应位置的水位和水流变化。HD 模型的主要特征可以概括为： 具有强大的前后处理功能，可以根据地形资料进行网格划分。 采用干湿动边界处理技术，可以更加准确地模拟滩地水流。 能定义流量、流速或水位等多种类型的水边界条件。 采用 ADI 计算方法，计算稳定性好，并且速度快 。 控制方程 澧南垸分洪洪水的演进模拟采用水力学方法进行模拟计算，这种方法基于平面二维非恒定流数值模型，包含二阶紊动项，其基本方程为： 数值求解方法 差分采用非迭代的交替方向的隐式(ADI)算法分步计算，所差分方程空间离散采用矩形交错网格。水面高程等标量定义于网格结点，流量分量定义在相应方向的网格线的中点。空间导数的有限差分近似值采用中心格式，时间上的积分采用时间中心隐式格式。 水流方程的离散和数值求解应按下列步骤进行： （1）将一个时间步长 △ t分为两个半步长，在每半个时间步长内，计算水位和x 向流速、y 向流速； （2）在 内，将水流方程离散为下三式 联立前二式，采用追赶法按 x 方向逐行求解u和h，然后利用后式逐点求出v。 计算方案 计算范围通过澧南垸的栅格图像来进行提取，包括堤垸内保证水位（46.07 m）以下的所有地区，在此范围内进行网格剖分。采用边长为 80m的矩形网格，网格总计13776个网格点。 计算时间步长取Δt = 30s。 初始条件：初始的高程由Arcmap软件生成，从DEM图中将每个网格点的地面高程提取出来作为初始的高程。 计算边界条件：入口边界条件为分洪闸流量过程，闭边界取无滑移边界条件；动边界采用干湿网络法进行处理。 特征点选取 在数值模拟产生的洪水数据中，提取16个特征点的水位、流量及流速等信息，来概述整个垸内的分洪过程、洪水特征。 特征线及特征点的位置 模拟结果 开闸分洪半小时后淹没范围与流场 开闸分洪1小时后淹没范围与流场 开闸分洪2小时后淹没范围与流场 特征点洪水特征 L1上4个特征点（L1a曹家巷、L1b丰坪村、L1c吴家弄子、L1d刘家屋场）的水位及流速变化： L2上8个特征点（L2a向家屋场、L2b前进村、L2c雷家岗、L2d罗家巷、L2e徐家岗、L2f乌鸦山、L2g陈家河、L2h廖家坪）的水位及流速变化： L3上4个特征点（L3a回龙村、L3b毛堰地、L3c何家溶、L3d魏家台）的水位及流速变化： 三、澧南垸分洪洪水三维淹没程序的开发与可视化 以ArcGis Engine的二次开发为基础，采用C#作为编程语言，在开发平台Visual Studio 2008里来开发程序。 基于Visual Studio2008的Arcgis Engine的二次开发界面 图像数据的前期处理 （1）澧南垸Tin图的创建。 根据区域的有限个点集将区域划分为相等的三角面网络，数字高程由连续的三角面组成，三角面的形状和大小取决于不规则分布的测点的密度和位置。 通过创建澧南垸的Tin图来作为澧南垸蓄洪区域的三维地形图。 （2）栅格图层的生成。 根据所获取的澧南垸Tin图，选取了从TIN中创建栅格表面的方法来创建栅格图层。 （3）澧南垸遥感图像数据的处理。 选取澧南垸遥感图像是为了和TIN数据叠加显示，以获取更好的虚拟现实及三维显示效果。遥感图像数据的处理大致包括投影坐标