1、华北水利水电大学毕业论文I目录摘 要 IABSTRACT： .II第一章 绪论 .11.1 问题的提出 .11.2 生产堤的介绍 .11.2.1 生产堤的作用 11.2.2 生产堤的由来及现状 11.3 生产堤的积极作用及危害 .21.3.1 生产堤的积极作用 21.3.2 生产堤的危害 31.4 研究的意义 .51.5 研究思路 .6第二章 河段的基本情况 .7第三章 模型的建立与率定 .93.1 MIKE 21 软件的介绍 .93.2 模型结构 103.2.1 平面二维数学模型的建立 .103.2.2 控制方程 .103.2.3 模型结构 .143.3 MIKE21 模型的建立 .153.
2、3.1 地形文件的生成 .153.3.2 区域边界的生成及调整 .153.3.3 模拟区域的三角划分 .163.3.4 初始糙率的设定 .193.3.5 夹河滩高村河段平面二维数学模型的建立 .20华北水利水电大学毕业论文II3.4 模型率定 26第四章 计算方案及结果分析 314.1 计算工况 314.2 模拟结果分析思路 314.3 成果分析 324.3.1 来洪时主槽水面高程分析 .324.3.2 滩区淹没面积分析 .324.4 认识 33第五章 结 语 .34参考文献 .35致 谢 .36附 录 .37华北水利水电大学毕业论文III黄河下游生产堤对滩区分滞洪的影响摘 要黄河下游滩区河道
3、主槽两侧修建生产堤缩窄了洪水过流断面,对防洪造成了一定影响。留存还是废除,再一次成为黄河下游防洪与减灾决策的焦点。本文以夹河滩高村河段为研究区域,利用 DHI 公司开发的 MIKE21 软件，建立平面二维数值模型,对现状生产堤和废除生产堤两种工况条件进行典型频率洪水演进的模拟。全文主要的研究内容：1）针对研究课题，收集黄河下游河段的基本资料。经整理分析得到该课题题目的来源及研究的意义。2）整理夹河滩高村河段上现存的生产堤的资料，得到生产堤的来由、建造标准及其存在的利弊。3）学习 MIKE 软件，建立夹河滩高村河段的平面二维数值模型。模拟计算出现有生产堤和废除生产堤两种不同工况下“96.8”洪水
4、的演进过程、河道流速、河道洪水水位。4）比较分析两种工况下，禅房、东坝头、石头庄、青庄和榆林五个测站的水面高程及洪峰流量下的淹没面积。进而综合分析生产堤存在的利弊，考虑生产生活的需要以及生产堤对防洪的影响等方面，对它进行统一规划。为了维持黄河河槽活力，协调人河关系，抑制二级悬河发展，黄河下游河道生产堤的治理方案是比较适宜的。 关键词：生产堤；防洪；洪水演进；数值模拟；MIKE21 模型中图分类号：TV871华北水利水电大学毕业论文IVINFLUENCE OF THE YELLOW RIVER DOWNSTREAM PRODUCTION OFBEACH DISTINGUISH FLOOD EMB
5、ANKMENTABSTRACT：The production dike which is built on the both sides of the lower reaches of the Yellow River main channel has some impact on the flood control. Retention or abolition, becomes the focal point of flood control and decision of the disaster reduction once again. In this paper, choosing
6、 Jiahetan to high village river as study area, a two-dimensional numerical model was established to simulate the evolution of the typical flood frequency on two conditions of production dykes and demolishing, using the DHI MIKE21 software. The main research contents:1) According to the research,coll
7、ecting the basic data of the Yellow River then getting the source and the significance of the research.2) organizing the related material of production dykes between Jiahetan and high village, getting the production levee reason, standards of construction and its existing stakes.3) learning MIKE sof
8、tware, establishing a two-dimensional numerical model of Jiahetan to high village river. Simulating and calculating the “96.8“ flood evolution process, river flow, river flood water level under the two different kinds of working conditions of the existing production dike and demolishing. 4) analysis
9、 the water surface elevation, the flow field and the flooded area of the peak flow on Buddhist temple, Dongba head, the stone village, Qing Zhuang and Yulin .Comprehensive analysis of the advantages and disadvantages of the production dykes, considering the need of production and life, besides the e
10、ffect of the production of embankment on the flood control, unified planning for it .In order to maintain the vitality of the Yellow River channel, coordinator, inhibition of two levels of elevated development, the control plan of the production dykes in the lower reaches of the Yellow River is feas
11、ible.Key Words: Production dykes; flood protection; advance of freshet; numerical simulation; Mike21 model华北水利水电大学毕业论文1第一章 绪论1.1 问题的提出沿黄群众早期为了保护农田和滩区，逼近黄河自修自守了土堤，一般顶宽为 23m,高度为 24m，它们一般为单个村或镇所建，一般为独立的、封闭的村堤、村堰，没有连成长堤。随着生产的发展和需要，他们在黄河两岸大量修建了生产堤。从而束窄河道并妨碍行洪，制约了槽滩的水沙交换，导致“二级悬河”加剧。生产堤的存在一方面保护了耕地、村庄和滩区，另一
12、方面也在特定的来水来沙条件下，提高水流挟沙能力。但是它也提高了河道水面高程，增加了黄河下游的防洪难度。所以本文就对生产堤存在的利弊，研究它对防洪的影响，在确保群众利益的同时，减小防洪难度。1.2 生产堤的介绍1.2.1 生产堤的作用黄河的生产堤是为了减少洪水淹没滩地的一种措施，滩地群众在黄河滩内顺着河槽修筑了大量的生产堤。用来保护农田和居住地，生产堤的修建的确给滩区的群众在生产及生活方面带来了一定的益处，使原本每年都要遭受淹没之苦，变为只有遇到多年一遇的较大洪水时才能漫滩。1.2.2 生产堤的由来及现状黄河下游滩区涉及河南、山东两省 15 个市、43 个县(区)，居住着 180 万人口，有25
13、0000hm2( 375 万亩)耕地 1。20 世纪 8090 年代，黄河下游洪水漫滩流量一般在5000m3/s 左右，而 90 年代以后，漫滩流量逐年变小。2002 年开始调水调沙后，河道有所冲刷，但部分河段漫滩流量仍不足 3500m3/s。沿黄两岸保护滩地的生产堤，早些年一般华北水利水电大学毕业论文2顶宽为 23m,高度为 24m。近年来受经济利益的驱动，滩区生产堤不断向主河槽外延，其区域扩大、防御标准提高，大部分高达 2.55m，甚至有些地方修筑的生产堤可以防御5000m3/s 以上的洪水。黄河滩内最早的生产堤是沿黄群众为保护耕地、村庄，自修自守的土堤，它一般逼近黄河，为单个村或镇所建，
14、一般为独立的、封闭的村堤、村堰，没有连成长堤。1958 年在“大跃进”形势下，基于错误估计了防洪形势，为保护两岸良田、村庄，在黄河两岸大修生产堤，从而束窄河道并妨碍行洪，制约了槽滩的水沙交换，导致“二级悬河”加剧。1974 年国务院制定了“从全局和长远考虑，黄河滩区应迅速废除生产堤，修筑避水台，实行一水一麦 2，一季留足群众全年口粮”的政策，在黄河主管部门和各级政府的共同努力下，生产堤破开了进水口门，到 1993 年，破口门长度已达生产堤长度的半数。1993 年后，由于黄河来水减少，滩区内群众防汛意识淡化，受短期经济利益的驱使，沿黄生产堤又重新修筑起来。1.3 生产堤的积极作用及危害1.3.1
15、 生产堤的积极作用1）保护了滩区群众生产和生活黄河下游滩区是群众赖以生产和生活的重要场所，也是群众最直接最根本的利益所在。生产堤对稳定滩区群众生产、生活有一定积极作用 3。(1)减轻了中小洪水漫滩对群众生活的影响。生产堤能够防止中小洪水漫滩，减少较大洪水漫滩几率，避免了群众的田地、房屋、道路、渠道等生产生活设施被淹，不但减轻了群众的灾害损失，而且避免了洪水围困村庄对群众的吃、住、行及就医、孩子上学等造成的不利影响，稳定了群众的生产生活，有利于滩区经济的发展。(2)增加了大水来临前群众外迁的时间。若没有生产堤，大洪水到来前，水流将从低洼处串水，洪水进入堤沟河后，阻断撤退道路，影响群众搬迁。有了生
16、产堤，将延缓洪水进入滩区的时间，增加群众外迁的时间。2）在特定的来水来沙条件下对冲刷河道、输送泥沙的积极作用对含沙量较低的中小洪水，若没有生产堤，可能造成漫滩。减小河道流量，相应降低了水流挟沙能力，不利于冲刷河道与输沙入海。而且由于洪水含沙量低，漫滩后达不到华北水利水电大学毕业论文3淤滩的目的。有了生产堤，能够防止中小洪水漫滩，有效地冲刷河道，即束水攻沙。对于高含沙中常洪水情况较为复杂，一般情况下不漫滩虽有利于输沙入海，但易加重主槽淤积。近几年调水调沙清水下泄，防止了漫滩分水，取得了良好的冲刷效果。1.3.2 生产堤的危害1)生产堤是“二级悬河”加剧的一个主要因素 据了解，黄河“二级悬河”最早
17、出现于上个世纪 70 年代。上个世纪 90 年代以来，由于气候干旱、来水减少及人类对黄河水资源的过度利用，黄河河道输沙用水被大量挤占，进入黄河下游的水量急剧减少，加剧了黄河水沙的不平衡；同时，黄河下游滩区群众在河道内普遍修筑生产堤，人为缩窄了行洪河道，影响了滩槽水沙交换，致使泥沙大部分淤积在生产堤内的主河槽里，加剧了“二级悬河”的发展。黄河水利委员会调查表明，上个世纪 50 年代，黄河下游河道中的泥沙 77淤在滩里，23淤在主河槽中，而目前，80的泥沙淤在主河槽中。在“二级悬河”形态下，河槽外高于滩面 4 米6 米，内呈“浅碟状”，即使发生中小洪水，主槽也难以容纳，必然造成重大河势变化，易出现
18、横河、斜河，增大了“冲决”和“溃决”的危险。2003 年 9 月，黄河河南兰考段出现流量 2400 多立方米每秒的“小水”，却酿成重大漫滩灾情。河南兰考县、山东东明县 114 个村庄、近 12 万人被洪水围困，而兰考段，恰恰是“二级悬河”最为严重的河段之一。2)“二级悬河”加剧了滩地横比降的发展生产堤的修建，加剧了滩地横比降的发展，但不同漫滩洪水情况下的影响程度差异很大。天然条件下的中小漫滩洪水，入滩水流少、含沙量低、滩槽水沙交换较弱，生产堤存在与否对河道冲淤及滩地横比降的影响不大。即使没有生产堤，因滩区分流比例较少，也难以造成滩区尤其是远离河槽的堤河低洼地带的大量淤积，影响范围主要集中在近生
19、产堤附近的滩面上，对滩地横比降总体影响较小。相反，生产堤的存在对中小洪水具有一定的约束作用，水流集中，对河道的排洪输沙具有一定的有利影响：主槽（深槽）冲刷作用可能会有所增强。但对于大漫滩洪水的影响则截然不同，控导工程和生产堤的存在对滩槽水流交换具有较大的影响。滩区淤积量的明显减少，加剧了滩地横比降的发展。使得滩唇高仰、堤根低洼的局面进一步加剧。大大增大了生产堤决口后，水流沿滩区最大比降（纵向和横向比降合成后的比降）方向、或沿着滩区串沟、汊河方向直冲大堤、华北水利水电大学毕业论文4造成冲决的可能性。如果没有生产堤，也存在堤防冲决的可能，但由于滩区漫水早、滩区水位随大河水位的上涨而近乎同步抬升，对
20、大堤的冲击强度相对较弱。由于河道横比降远大于纵比降，如若中、小洪水生产堤一旦决口就会形成一条小河直冲大堤。今后中、小洪水甚至大洪水时的横河、斜河、滚河的情况均有可能产生。生产堤的修建使中、小洪水行洪的宽度大大缩小。而生产堤使行洪河槽缩窄，泥沙的淤积范围缩小，仅限于生产堤之间，加速主河槽淤积抬高的速率，增大滩唇与大堤根之间的断面横比降，这与黄河下游“宽河固堤”方针是相违背的。3)生产堤的修建加大了黄河防洪的难度近年来由于黄河多次加高加固黄河堤防，兴建河道整治工程，特别是 1998 年长江大水之后，国家加大了黄河治理力度，按照 2000 年设防标准，加高、加固两岸堤防，加修、改建了许多河道整治工程
21、，近期又进行了标准化堤防建设，将黄河堤防建成“防洪保障线、抢险交通线、生态景观线”，防洪工程的抗洪能力大大增强。但是如果任由“二级悬河”的发展，任由生产堤的加固，槽滩高差的加大，一旦生产堤决口，那么形成的串沟、滚河等险情将极大的影响防洪的安全。4)河槽萎缩，排洪能力降低由于黄河上游来沙去向主要是输送入海和入海口造陆、淤积在河道内，还有一小部分人工引放至黄河灌区及固堤，由于多年来黄河没有发生较大洪水，泥沙集中淤积在主河槽内，河道主槽淤积比例由 80 年代初的 30增加到 70，河槽严重萎缩，河道排洪能力降低 4。目前，开封以下河段横比降达 45，而河道纵比降为 1.71.8，一旦发生较大洪水，由
22、于河道横比降远大于纵比降，滩区过流比增大，极易发生横河、斜河，甚至发生“滚河”，产生顺堤行洪，危及堤防安全。由于河槽过洪能力明显降低，局部河段 2000m3/s 左右流量就可漫滩，小水就能造成大灾。5)严重影响洪水传播和洪水的准确预测黄河下游河道，在目前的边界条件下，洪水演进受生产堤阻隔，断面缩窄，水位雍高。滩区漫滩进水时间不仅受洪水量的影响，更受生产堤强度及防守效果的影响。人为因素较多地改变了自然状态下的洪水演进规律，对生产堤破口、决口时间、地点，漫滩水量及漫滩水如何演进、衰减等都难以预测，很难准确地进行洪水预测、预报，从而带来较大的测报误差，进而影响到防汛调度和决策。如“96.8”洪水在花
23、园口利津河段传播时间长达 369 h，约为多年均值的 2 倍，为有水文记录以来最长。如孙口断面“82.8”洪水漫滩，8 月 6 日 20 时至 7 日 11 时洪峰前后，洪峰流量报汛误差最大达 2 000m3/s，相华北水利水电大学毕业论文5对误差 20%，致使本应为 10100m3/s 的孙口洪峰流量，错报为 8100m3/s，严重影响了东平湖分洪的决策和防汛指挥调度 5。6)影响河道整治工程的安全，增加了河道整治的难度黄河下游河道整治工程是按照规划的河道治导线进行的，而生产堤的修建是一种民事行为，随意性较大，平面布置主要是依据护滩、保村的原则修建的，很多河段布置违背了河道洪水演进的规律，有
24、的生产堤甚至修建在主河槽内，严重干扰了洪水的正常演进，造成了河道整治工程条件恶化，易出现较大险情。同时，生产堤加剧了二级悬河的发展，使平滩流量减小，增大了河道整治的难度。7)对滩区群众生产生活也存在不利影响生产堤虽然保护了滩区群众免于中小洪水频繁漫滩受灾，但对群众的生产生活也有许多不利影响。主要表现为:漫滩水加重了对滩区群众安全的威胁。生产堤决口后水流集中，冲刷力强，对村庄、桥梁、渠道等群众生活、生产设施造成较大危害，加重漫滩损失。堤沟河得不到及时淤平，一旦漫滩或遇暴雨，滩内积水不易排出，不仅秋季作物不保甚至连冬麦也无法播种。不利于滩区种植和土壤改良。影响滩区群众及时迁移安置。群众易依赖生产堤
25、而疏于迁安，一旦生产堤溃决，洪水直冲地势较低的堤根，顺堤成河，此时开始迁移，道路已断，仅靠少量渡船外迁只能勉强救命。生产堤口门处粗沙沉淀不利于耕种。1.4 研究的意义现状生产堤严重影响防洪和滩槽淤积分布，但废除生产堤影响滩区 180 余万人的生存与发展，地方政府阻力很大。参考其它河流如淮河、长江的滩区治理模式，可以考虑废除小滩区和严重影响行洪安全的生产堤，设置 1.54.0km 的 8000m3/s 及以下洪水通道，构筑大滩区（20km 2以上）的滞洪区，有计划分洪淤滩，找出黄河下游滩区治理中经济发展与分洪安全的结合点，解决滩区 85%人口生存与滩区淤积分布问题。本项目采用平面二维数学模型，依
26、据 “96.8”典型洪水情况，按现状生产堤和废除生产堤两种条件进行了典型洪水的模拟，模拟计算给出了夹河滩高村河段不同边界洪水演进、漫滩的过程及滩区淹没程度。通过对模拟结论的分析，从而进一步验证黄河下游生产堤对滩区分滞洪的影响以及给出今后对生产堤的管理建议。华北水利水电大学毕业论文61.5 研究思路应用平面二维数学模型，按现状生产堤、废除生产堤两种条件进行了典型洪水 “96.8”的模拟，模拟计算给出了夹河滩高村河段不同边界洪水演进、漫滩的过程及滩区淹没程度。以平面二维数学模型模拟计算成果为基础，从而对比不同模式、不同工况下，河道水面高程、滩区淹没面积和洪水演进过程等，提出更加合理的生产堤的运用原
27、则、要求和方案。以平面二维数学模型模拟计算成果分析滩区淹没状况及相关迁安救护问题。研究工作流程见图 2图 2 研究工作流程图华北水利水电大学毕业论文7第二章 河段的基本情况20 世纪 50 年代黄河下游河道基本属于天然情况，60 年代以后滩区农民自发性地逐渐修建生产堤，对滩地农田有一定保护作用，故屡禁屡建。19871999 年期间，黄河下游河道经常出现小水带大沙的异常来水来沙条件。由于生产堤阻碍了滩槽水流泥沙的横向交换，泥沙淤积主要集中在主槽和嫩滩上，广大滩区淤积很少。滩槽淤积分布的不均匀性形成了滩唇高仰、堤根低洼，临河滩面明显低于滩唇，背河地面又明显低于临河滩面的的不利局面，形成所谓“二级悬
28、河”。河道典型断面状况如图 3 所示。 图 3： 黄河下游典型大断面形态示意黄河在孟津县白鹤由山区进入平原，河道上宽下窄，比降上陡下缓，花园口至高村河段，河宽水散，冲淤幅度大，主流摆动频繁，为典型的游荡性河段。两岸大堤堤距一般为 510km，最宽处达 20 多 km。由于下游河道水少沙多，河道淤积严重，长期以来“善淤、善徙、善决”而著称于世。人民治黄以来，通过大规模的整治，有效控制了河势摆动范围，减小了横河、斜河对地方冲决的威胁。陶城铺以下河道的河势已得到控制；高村至陶城铺河段，主流已趋于稳定，河势得到了基本控制；高村以上游荡性河段也已修建了一批布点工程，缩小了游荡摆动范围。经过多年洪水的考验
29、，河道整治工程在控导主流、护滩保堤、争取防洪主动方面发挥了重要作用。各河段具体的基本情况 6：孟津铁谢至京广铁路桥河段。该河段为禹王故道，长 93km。由于大河初出峡谷，骤入平原，比降较陡，流速较大。此河段过去主流摆动频繁，坍塌剧烈，滩地迁徙无常。70 年代以后，修建了部分控导护滩工程，滩地相对稳定。目前此河段滩地面积约519km2，其中耕地 2.72 万 hm2，有村庄 58 个，6.8 万人。这一河段滩槽高差较大，一般华北水利水电大学毕业论文8当地洪峰流量 8000m3/s 可漫滩，滩地主要集中在左岸的孟津市、温县、武陟县境内，面积广大，习惯上称为“温梦滩”。京广铁路桥至东坝头河段。该河段
30、河道长 124km，为明清河道，两岸均有黄河大堤，河道内有 1855 年铜瓦厢决口后河床下形成的“高滩”。一百多年来，由于主流摆动，部分高滩坍塌。现有滩地面积 844km2，其中高滩 452km2，低滩 392km2。高滩土地生产相对比较稳定，粮棉单产较高，村庄稠密，有 426 个村庄，35.7 万人，耕地 5.17 万 hm2。该河段自铜瓦厢决口以来，河槽淤高 2m3m，目前花园口水文站 15000m3/s 的洪水将大部漫滩，高滩已不高。东坝头至陶城铺河段。该河段长 227m，是铜瓦厢决口改道后形成的河段。一百多年来，此河段游荡多变，处于强烈的堆积状态，一般淤高 4m5m，特别是东坝头至高村
31、河段，滩槽高差小，主槽位置常有摆动，俗称“豆腐腰”河段。该河段堤距大，滩面广，滩地面积 1738km2，其中耕地 10.97hm2，是黄河下游的主要削峰区。有村庄 1140 个，84.66 万人。长垣县、东明县的滩区约占该河段滩区总面积的 42%。现在，4000m 3/s 的洪水有的串沟可能过流，局部漫滩受淹，8000m 3/s 的洪水可大部漫滩，10000m 3/s 以上的洪水滩地全部淹没，迁安任务十分艰巨。该河段漫滩机遇较多，生产发展缓慢。堤河、串沟比较多，在防洪方面，还存在着大洪水时夺流滚河，顺堤行洪的危险。陶城铺至西河口河段。该河段长 345km，是铜瓦厢改道后夺大清河演变形成的。此河
32、段已治理成弯曲性河道，河势流路比较稳定，滩槽高差较大。除长清、平阴两县的滩区为连片的大滩区外，其余全部是小片滩地。滩地面积 855.7km2，其中耕地 6.11 万 hm2，有 544 个村庄，41.6 万人。该河段的滩区，不仅伏秋大洪水漫滩几率高，而且还受凌汛漫滩的威胁，生产不稳定。西河口以下河段。长 50 多公里，属河口地区，有部分油田设施。华北水利水电大学毕业论文9第三章 模型的建立与率定3.1 MIKE 21 软件的介绍MIKE21 数值计算与分析软件是国际上比较成熟的 DHI 软件系列中关于水动力、波浪和泥沙输运等模型进行潮流场、代表波要素的波浪场和泥沙输运的数值模拟工具。MIKE2
33、1 是一个专业的工程软件包，用于模拟河流、湖泊、河口、海湾、海岸及海洋的水流、波浪、泥沙及环境。目前该软件在国内的应用发展很快，并在一些大型工程中广泛应用，如:长江口综合治理工程、杭州湾数值模拟、南水北调工程、重庆市城市排污评价、太湖富营养模型、香港新机场工程建设、台湾桃园工业港兴建工程等。它具有以下优点：（1）用户界面友好，属于集成的 Windows 图形界面。（2）具有强大的前、后处理功能。在前处理方面，能根据地形资料进行计算网格的划分;在后处理方面具有强大的分析功能，如流场动态演示及动画制作、计算断面流量、实测与计算过程的验证、不同方案的比较等。（3）可以进行热启动，当用户因各种原因需暂
34、时中断 MIKE21 模型时，只要在上次计算时设置了热启动文件，再次开始计算时将热启动文件调入便可继续计算，极大地方便了计算时间有限制的用户。（4）能进行干、湿节点和干、湿单元的设置，能较方便地进行滩地水流的模拟。（5）具有功能强大的卡片设置功能，可以进行多种控制性结构的设置，如桥墩、堰、闸、涵洞等。（6）可以定义多种类型的水边界条件，如流量、水位或流速等。（7）可广泛地应用于二维水力学现象的研究，潮汐、水流、风暴潮、传热、盐流、水质、波浪紊动、湖震、防浪堤布置、船运、泥沙侵蚀、输移和沉积等，被推荐为河流、湖泊、河口和海岸水流的二维仿真模拟工具。水动力(HD)模块是 MIKE21 软件包中的基
35、本模块，它为泥沙传输和环境水文学提供了水动力学的计算基础。HD 模块模拟湖泊、河口和海岸地区的水位变化、和由于各种力的作用而产生的水流变化。在为模型提供了地形、底部糙率、风场和水动力学边界条件等数据后，模型会计算出每个网格的水位和水流变化。华北水利水电大学毕业论文101）MIKE 21 模型的子模块主要有以下 4 个领域： 海岸水文学和海洋学 环境水文学 泥沙传输过程 波浪2）MIKE 21 主要应用 河口海岸结构物设计数据的评价 港口布局和海岸保护措施的优化 冷却水、海水淡化及再循环分析 河口海岸及海洋结构物的环境影响评价 海上安全操作和航行海情预报 沿海洪水和风暴潮预警 内陆洪水及坡面流模
36、拟3.2 模型结构3.2.1 平面二维数学模型的建立该数学模型已在一些复杂工程中得到成功应用，能较准确地模拟、预测在修建工程前后计算域内洪水演进、建筑物过流状况及河床变形状况。3.2.2 控制方程1)控制方程特点对于河口、海岸、湖泊、宽浅河道与坡面流等广阔的宽浅型水域，水平尺度远大于垂向尺度，水力要素（流速、水深等）在垂直方向的变化要远小于水平方向的变化，沿水深分布比较均匀，其流态、水力要素可用沿水深平均值来表示。采用水深平均的平面二维水沙运动方程可较好地反映这类流场中水流泥沙运动特征，较好地解决一些宽浅水域的河流模拟计算问题。平面二维水流运动所遵循的基本方程是由三维时均雷诺方程沿水深进行莱布
37、尼兹积分得到，在运动方程中以混长紊流模型求解紊动切应力。本次采用华北水利水电大学毕业论文11的平面二维数学模型综合考虑洪水演进、河床变形等河流演变过程，可以全面模拟计算域内水流及河床变形过程。2)控制方程MIKE 21 水动力学模型的主要控制方程为雷 诺 方 程 的 垂 向 平 均 如下：质 量 守 恒 方 程 ：（ 3-2-S=yq +xpt1）X 方 向 动 量 方 程 ：（ 3-2-ixaxxyxP SphfVqhCqgpt 1 222）Y 方 向 动 量 方 程 ：（ 3-2-iyayxyyq SphfVqhCqgpt 1 223）方程中 为时间 ； 、 为右手 Cartesian 坐
38、标系； 为水面相对于未扰动水面的高度tSxy即通常所说的水位； 为静止水深； 为重力加速度 ; ， 分 别 为hg2/smtyxp,tq,方 向 和 方 向 的 通 量 ， 、 分别为流速在 、 方向上的xyms/3vhqup,其 中 uxy分量； 为当地大气压； 为水密度； 为 Coriol 系数= ， 为地球自转apyxsin角速度， 为计算点所处的纬度，一般取 ；C 为 柯 西 阻 力 系 数140729.s（ m1/2/s） ； = 为风 摩 擦 因 素 函 数 ， 为风应力系数， 为空气密度，V 为 风Vfa2 aa速 ， 、 为 方 向 和 方 向 的 风 速 ； 为 源 项 ，
39、为 源 项 在 方 向 的 分 量 ，xyysm/Sixx为 源 项 在 方 向 的 分 量 ； 、 、 为 各 方 向 上 的 粘滞应力项。iySxy3）方程组求解华北水利水电大学毕业论文12在空间上，由于网格类型为非结构网格，采用中心有限体积法对原方程进行离散，把整体的计算区域细分为非重叠的单元；时间上采用显示积分。有 限 体 积 法 （ FVM） 又 称 为 控 制 体 积 法 。 其 基 本 思 路 是 ： 将 计 算 区 域 划 分 为 一 系 列不 重 复 的 控 制 体 积 ， 并 使 每 个 网 格 点 周 围 有 一 个 控 制 体 积 ； 将 待 解 的 微 分 方 程 对
40、 每 一个 控 制 体 积 积 分 ， 便 得 出 一 组 离 散 方 程 。 其 中 的 未 知 数 是 网 格 点 上 的 因 变 量 的 数 值 。为 了 求 出 控 制 体 积 的 积 分 ， 必 须 假 定 值 在 网 格 点 之 间 的 变 化 规 律 ， 即 假 设 值 分 段 的分 布 剖 面 。从 积 分 区 域 的 选 取 方 法 看 来 ， 有 限 体 积 法 属 于 加 权 剩 余 法 中 的 子 区 域 法 ； 从 未 知解 的 近 似 方 法 看 来 ， 有 限 体 积 法 属 于 采 用 局 部 近 似 的 离 散 方 法 。 简 言 之 ， 子 区 域 法 属于
41、 有 限 体 积 发 的 基 本 方 法 。 就 离 散 方 法 而 言 ， 有 限 体 积 法 可 视 作 有 限 单 元 法 和 有 限 差分 法 的 中 间 物 。 有 限 体 积 法 只 寻 求 的 结 点 值 ， 这 与 有 限 差 分 法 相 类 似 ； 但 有 限 体 积法 在 寻 求 控 制 体 积 的 积 分 时 ， 必 须 假 定 值 在 网 格 点 之 间 的 分 布 ， 这 又 与 有 限 单 元 法相 类 似 。 在 有 限 体 积 法 中 ， 插 值 函 数 只 用 于 计 算 控 制 体 积 的 积 分 ， 得 出 离 散 方 程 之 后 ，便 可 忘 掉 插 值
42、 函 数 ； 如 果 需 要 的 话 ， 可 以 对 微 分 方 程 中 不 同 的 项 采 取 不 同 的 插 值 函 数 。有 限 体 积 法 的 基 本 思 路 易 于 理 解 ， 并 能 得 出 直 接 的 物 理 解 释 。 离 散 方 程 的 物 理 意义 ， 就 是 因 变 量 在 有 限 大 小 的 控 制 体 积 中 的 守 恒 原 理 ， 如 同 微 分 方 程 表 示 因 变 量 在无 限 小 的 控 制 体 积 中 的 守 恒 原 理 一 样 。 有 限 体 积 法 得 出 的 离 散 方 程 ， 要 求 因 变 量 的 积分 守 恒 对 任 意 一 组 控 制 体 积
43、 都 得 到 满 足 ， 对 整 个 计 算 区 域 ， 自 然 也 得 到 满 足 。 这 是 有限 体 积 法 吸 引 人 的 优 点 。 有 限 体 积 法 即 使 在 粗 网 格 情 况 下 ， 也 显 示 出 准 确 的 积 分 守 恒 。在进行计算之前，首先要将计算区域离散化，区域离散化的实质是对空间上连续的计算区域进行划分，把它划分成许多个子区域，并确定每个区域中的节点，从而生成网格。然后，将控制方程在网格上离散，即将偏微分格式的控制方程转化为各个节点上的代数方程组，然后在计算机上求解离散方程组，得到节点上的解。节点之间的近似解一般可以认为光滑变化，原则上可以应用插值方法确定，从
44、而得到定解问题在整个计算区域上的近似解。因此当网格节点很密时，离散方程的解将趋近于相应微分方程的精确解。此外，对于瞬态问题，还需要涉及时间域离散。应用控制容积法导出离散方程的主要步骤： （1）将守恒型的控制方程在任一控制容积及时间间隔内对空间和时间做积分； 华北水利水电大学毕业论文13（2）选定未知函数及其导数对时间和空间的局部分布曲线，即“型线”（profile）。也就是如何从相邻节点的函数值来确定控制容积界面上的被求函数值的插值方式； （3）对方程各项按选定的型线做出积分并整理成节点上未知值的代数方程。图 4 中使用的是三角形控制体积，三角形的质心是计算节点，如图中的红点所示。由于使用了显
45、式迎风格式，时间步长要求严格满足 CFL1。CFL 的计算公式为：ytvghxtughCFL（ 3-2-4）针对所有的计算网格，在控制体积 P 及时间段 （时间从 到 ）上对控制方程积ttt分有：（ 3-2-5）VtVt tt Sdtdgraivuiv对于上式中的瞬态项，假定物理量 在整个控制体积 上均具有节点处的值 ，同PP时假定密度在时间段上的变化量极小，则式(3 5)中的瞬态项变为：（ 3-2-6）VdVtdVt Ptt 0在上式中，上标 表示物理量在时刻 的值，而在 时刻的物理量没有用上标来0tt标记，下标 表示物理量在控制体积 的节点 处取值。PP对于源项：（ 3-2-7）dtVSV
46、dtSVdtSdt PCtPCttVt 对于对流项：根据 Gauss 散度定理，将体积分转变为面积分后，有：（ 3-2-8） dtAvvAutdVtiv snwet stVt 图 4华北水利水电大学毕业论文14式中，A 是控制体积界面的面积。扩散项：同样根据 Gauss 散度定理，将体积分转变为面积分后，有：（ 3-2-dtxAxxAx tdVtgraiv sSPsnPNnwWPePEet setVt 9）在得到（ 3-2-5） 各项的单独表达式后，还要做两方面的工作：1）在对流项中需要引入特定的离散格式（ 3-2-8） 中界面物理量 、 、 和 用节ewns点物理量来表示，可用一阶迎风格式。
47、 2）在对流项，扩散项和源项中引入全隐式的时间积分方案，如 。这tdPt样，方程（ 3-2-5） 变为：（ 3-2-10）baaaNSEWP这就是在全隐式时间积分方案下得到的二维瞬态对流-扩散问题的离散方程。式中系数 、 、 和 取决于在对流项中引入的特定离散格式。使用一阶迎风格式有：WaENS（ 3-2-tVaSb VSaFFaDaPPC PsnweNSEWnnNssSeeEwwW 00 0,mx,0a11）3.2.3 模型结构1）计算方法 华北水利水电大学毕业论文15由于网格类型为非结构网格，采用中心有限体积法对原方程进行离散，把整体的计算区域细分为非重叠的单元。该方法的优点为计算速度较快
48、，非结构网格可以拟合复杂地形。2）初始条件及边界条件（1）初始条件对于给定的计算区域，在时间 t=0 时，令 t=0= 0（ x， y）， u t=0=u0（ x， y）， v t=0=v0（ x， y）（2）边界条件a）进口边界条件:给出进口开边界处的流量过程（包括上游大河进口条件与侧向入汇进口条件）：Q（ t）= Qopb（ t）上式中 Qopb 分别为开边界上已知的水位或流速分量以及流量，一般根据计算区域以上的产汇流模型计算或由实测水文资料确定。b）出口边界条件：出口开边界有两类：一类是自然开边界，主要是经下边界或侧边界出流的河流，可按实测水文资料（水位流量关系）确定。如无实测资料，则按附近河道纵坡，以均匀出流考虑。c）陆地边界:根据流体固壁不可穿越的原理，在不考虑渗透的情况下，可以认为陆地边界上法向速度为零；根据水流无滑动原理，水体在陆地边界上的切向流速也应为零。3.3 MIKE21 模型的建立3.3.1 地形文件的生成地形处理是建模时的第一步，它是需要在建模前首先完成的工作，地形处理的好坏非常关键，关系到模型搭建的成败。此次地形处理包括夹河滩高村河段的主槽和滩地的地形处理。创建计算网格需要河段的地形数据资料。夹高河段的主
