基于MIKE FLOOD洪泽湖周边滞洪区洪水演进模拟.docVIP

基于MIKE FLOOD洪泽湖周边滞洪区洪水演进模拟 　　摘要：基于实测断面资料建立了研究区的一维水动力模型，基于高精度DEM以及1：10000地形图建立了研究区的二维水动力模型，并用MIKE FLOOD将一维模型和二维模型进行耦合，构建了洪泽湖周边滞洪区一、二维耦合的洪水演进数学模型。利用2003年历史洪水资料对模型参数进行了率定，并以2007年历史洪水资料进行了验证。以洪泽湖百年一遇设计洪水为模型上边界，二河闸、三河闸以及高良涧闸的现行调度方案的水位-流量关系为模型下边界，对洪泽湖百年一遇设计洪水方案进行模拟计算，当蒋坝水位达到14 33 m时，洪泽湖周边滞洪区开始滞洪，得到开始滞洪后不同时段研究区内各类洪水风险要素的动态分布情况以及最大淹没水深、淹没历时，验证了模型的合理性，可用于蓄滞洪区洪水演算分析。 　　关键词：蓄滞洪区；洪水演算；洪泽湖；风险要素；MIKE FLOOD 　　中?D分类号：TV 122 文献标识码：A 文章编号：1672 1683（2017）05-0056-07 　　蓄滞洪区是我国防洪减灾的重点对象，适时启用蓄滞洪区削减洪峰、分蓄洪水是保障下游和防洪重点保护区域安全的有效措施。为做好防洪减灾工作，我国高度重视防洪规划、防洪减灾工程建设与非工程措施部署。其中，洪水演进模拟是非工程措施的重要组成部分，是洪水风险分析、洪水风险图绘制的重要依据。洪水演进模拟的研究方法主要有差分法、有限元法、有限体积法等。对于洪水数值模拟，国内外已有不少专家利用有限差分法和矩形网格建立了二维洪水演进模型，但是独立的二维模型无法计算获得河道中某个断面、某个时刻的水位和流量等水文要素，尤其在研究区内水工建筑物众多、控制调度复杂的情况下，更是难以高精度地模拟水库调度、溃口等对河道洪水演进的影响。 　　为了更好更精确的反映决堤后洪水演进的情况，本文拟采用以一、二维耦合的水力学方法，依据水量守恒及水动量守恒原理建立洪水数值模型，模拟洪水的恒定及非恒定演进过程，分为一维水力学法和二维水力学法。其中一维水力学法是指采用一维非恒定流微分方程即圣维南方程组对河道非恒定流的水动力学进行模拟，该方法能够获得每个模拟时刻河道内的水面线、流量等情况，便于检查和验证，也便于过程控制。二维水力学方法采用浅水方程来模拟洪水在淹没区域内耕地、林地、山丘区、道路等演进的情况，与一维数学模型相比，二维数学模型能够提供更加详细的水情信息，如：淹没范围、淹没水深、淹没历时等情况。目前世界上最为广泛应用的丹麦水利研究所（Danish Hydraulic Institute）研发的MIKE软件，包含MIKE11、MIKE21、MIKEFLOOD等多个功能强大的模块，其中MIKEFLOOD模块可以模拟一维河网水动力学系统以及二维的洪泛区，实现一维和二维区域之间自由的水体交换，适用于宏观上的流域控制性工程规模论证分析和流域洪水调度研究等。本文以洪泽湖地区为例，利用MIKE FLOOD建立一、二维耦合的洪水演进数学模型进行模拟研究。 　　1研究区概况 　　洪泽湖是淮河流域最大的湖泊型水库，地处苏北平原中部偏西，位于淮河中下游结合部，承泄淮河上、中游15.8万km3的来水，总库容135亿m3。地貌上，湖西、湖南为低山岗阜，湖西岗陇和洼地宽窄不等、“三洼四岗”高低相间；湖南区蒋坝至盱眙县城是连绵的低山，湖东、湖北由河湖冲刷堆积而成的平原，地势低下。洪泽湖整体呈西北高东南低的形态。洪泽湖周边滞洪区位于洪泽湖大堤以西，废黄河以南，泗洪县西南高地以东，包含盱眙县的沿湖、沿淮地区，地面高程在16.83 m（1985国家高程，下同）以下，涉及淮安、宿迁两市的淮阴区、洪泽县、盱眙县和宿城区、泗洪县、泗阳县六县区部分地区，还有省属洪泽湖、三河农场，滞洪区面积2132 k m3。 　　洪泽湖入湖河流主要在湖西部，有淮河、怀洪新河、新老濉河、新老汴河和徐洪河、安东河等，在湖北侧入成子湖的河流有古山河、五河、肖河等，南侧主要入湖河道为维桥河、高桥河，淮河入量占流入总量的70%以上。区域内地势外高内低，唯一排水出路即洪泽湖湖区，洪泽湖洪水由三河闸、二河闸、高良涧闸分别经入江水道、入海水道、分淮入沂、灌溉总渠入江入海。区域防洪工程包括外围防洪、内部防洪工程，主要为迎湖挡洪堤和河道堤防。入湖控制建筑物主要有团结闸、高松河闸、古山河闸、西民便河闸、安东河闸、濉河闸、老汴河闸等。出湖控制建筑物主要是三河闸、二河闸、高良涧闸；区域内防洪工程主要有滞洪区间内部河道堤防和排涝涵闸、泵站。详见图1。 　　2模型构建 　　2.1模型控制方程 　　描述一维水流运动的Saint-Venant方程组是建立在质量和能量守恒的基础上的，以水位和流量为研究对象，其表达式为： 　　2.2河道一维模型 　　一