洪水与风暴潮共同作用下溃堤洪水一维二维耦合模型及应用.docVIP

洪水与风暴潮共同作用下溃堤洪水一维二维耦合模型及应用 　　摘要：为了解决复杂条件下河道洪水漫堤、溃堤和潮水倒灌的问题，建立了在洪水和风暴潮共同作用下的天然河道漫溃堤洪水在防洪保护区的一二维水动力耦合模型。在兼顾计算效率和精度的基础上，分区剖分网格，设置合理的网格面积；根据植被和地物的不同，分区设置糙率；通过宽顶堰的方式将河道一维模型和防洪保护区二维模型进行侧向耦合；采用基于侵蚀的渐变溃破坏方式模拟河道溃堤过程，实现了河道上游发生洪水和下游遭遇外海风暴潮的复杂情况下精细化模拟。实例研究表明，所建模型可灵活处理复杂多变的水力条件，模拟结果合理可靠，可为防汛部门制定决策提供有力的科学依据。 　　关键词：一二维耦合模型；溃、漫堤洪水；风暴潮；渐变溃；太阳河 　　中图分类号：TV122.4 文献标识码：A 文章编号：1672-1683（2017）06-0043-07 　　近几年，我国注重综合采用工程措施和非工程措施应对洪涝灾害，在全国各地推行洪水风险图编制项目，使得我国的水动力数学模型有了很大的进展。在我国，黄金池等利用模型代替了多个单独使用的水动力学模型，模拟了土石坝溃坝后的洪水演进过程。宋霁云等描述了模型参数和模型效率间的关系，为模型参数的合理设置提供了理论依据。李大鸣等建立二维水动力学模型研究了桥涵概化尺度对洪水演进的影响。陈文龙等通过构造并求解Riemann问题实现一维二维模型耦合。苑希民等建立了堤洪水的二维水动力学模型，模拟了黄河内蒙段南岸灌区溃堤洪水流速、流态及水深的变化情况。Lai等提出了适用于大尺度水动力模拟的一维二维耦合方法，为模型的耦合提供了理论基础。孙秀丽等采用标准连接的方式，建立了一维、二维耦合模型模拟溃堤洪水。而在国外，Paul D.Bates等学者基于浅水动力学模型的改进理论，研究出一种简单高效的二维水动力学模型。Maarten Breckpot等利用水文模型解决洪水问题，推动了洪水模型的多元发展。Morales-Heril ndez等模拟并分析了浅水流的演进情况，系统地研究了水动力模型的机理，Dewals等将二维模型和集总模型结合模拟水库的溃堤。Dushmanta等运用有限差分法模拟了河流洪水漫堤过程。以上研究中，主要针对水库溃坝、河道溃堤洪水的模拟，并没有考虑风暴潮倒灌加剧河道洪水的影响。而我国沿海地区往往河流上游洪水伴随着区间降雨，入海口由于风暴潮的作用致使潮位顶托河道洪水，使洪涝灾害更加严重。而且在传统建模过程中，对地形、河道以及某些参数进行过度概化，造成了模拟结果的不准确性，如何准确地模拟洪水的发生过程成为当前亟需解决的问题。基于此，本文建立上游洪水叠加外海风暴潮的一二维耦合精细化模型模拟复杂水力条件下洪水的漫堤、溃堤以及在防洪保护区的演进过程，为沿海地区制定正确防灾减灾措施提供科学依据。 　　1一维、二维水动力耦合模型原理 　　1.1一维水动力模型 　　单一河道非恒定流控制方程采用圣维南方程组，圣维南方程是反映有关物理定律的微分方程，包括连续方程和动量方程。 　　一二维耦合模型计算时，为提高计算效率，二维模型应选择基于CFL（Courant Friedrichs Lewy）条件的自适应时间步长模式。为有效解决一维模型和二维模型时间步长不一致问题，以一维模型的固定时间步长为基准，二维模型时间步长进行自适应匹配。 　　2研究实例 　　2.1研究区概况 　　本次以海南省太阳河为研究对象。太阳河流域系海南岛暴雨区，年平均来水量6.6亿m3，是万宁市最大的河流。万宁市位于海南省东南部，东濒南海，西毗琼中，南邻陵水，北与琼海接壤。万宁水库将太阳河分为上、下游两段，本研究区上游白万宁水库起，下游到太阳河入海口，流域面积120 km2。图 　　2为太阳河水系图。 　　2.2精细化模型的构建 　　2.2.1一维河道模型 　　河道断面是对一维模型精确度影响最大的因素。由于河道洪水传播的不稳定性，因此使用尽可能多的横断面以尽可能准确地描述河道地形是非常必要的，特别是在横断面变化迅速的地方更应布置更多的断面。在一维模型创建过程中，增加河道的断面数量，减少两断面间的?g隔里程，将会极大地提高模型的精确性。本研究区域河长16.2 km，共设置70个断面，平均断面间距230 m。图3为太阳河下游段河道断面布置图。 　　2.2.2二维模型的构建 　　网格面积的大小是影响二维模型精确度的重要因素。通常情况下，网格面积越小，模型的计算结果更精确，但这将极大地延长模型的计算时间；而减少二维网格量，则不能真实反映地物地貌，计算误差较大。考虑模拟结果的精确性又兼顾模型的计算效率，本模型采用了分区剖分网格的方式。在靠近河道和堤防的位置，水流交换频繁，所以进行网格细化处理；而远离河道的二维区域，适当地增大网