黄河口潮流与泥沙输移过程的数值研究 Tidal current and its sediment transport capacity in Yellow River Estuary.pdfVIP

水 利 学 报 2008年10月 SHUILIXUEBAO 第39卷第10期 文章编号：0559．9350(2008)10-1256—07 黄河口潮流与泥沙输移过程的数值研究 王崇浩，曹文洪，张世奇 (中国水利水电科学研究院泥沙研究所，北京100044) 摘要：黄河口的海洋动力输沙能力受径流、潮流和波浪等动力条件的影响，是个十分复杂的问题。本文利用MIKE21 模型建立了渤海水域大模型和黄河口水域的小模型，大模型为小模型提供外海边界条件。利用渤海水域验潮站的分 潮资料对大模型进行了验证，计算值与实测值符合良好。黄河口水域小模型计算研究了多个人海流路口外水域的潮 流过程和潮流输沙能力，计算结果得到黄河口水域新、老口门附近存在三个高流速区，最大潮流速在1．7m／s左右，自 渤海湾向莱洲湾的强大沿岸流输沙是黄河口潮流输沙的主要形式，分析计算El外泥沙净通量得到刁口河、原河道和 现行清8流路、马新河、北汉流路和十八户流路口外水域向口门两侧lOkm外和外海的输沙量占利津沙量百分数分别 t,0．72亿t和0．40亿t。计算得到的清水沟原河道的潮流输沙能力与多年资料分析结果基本一致。 关键词-黄河口；入海流路；数值模拟；潮流过程；输沙能力 中图分类号：TVl48 文献标识码：A 1研究背景 黄河以多沙闻名于世，大量泥沙被输送至河口地区，造成尾间河道的淤积抬高和1：3门的延伸n一】， 严重影响了黄河下游防洪和经济发展。泥沙入海后，在潮流和波浪等动力的作用下，向外海和沿岸两侧 输移、扩散，部分被带往深海。然而，黄河口海洋动力输沙是一个复杂的过程，受诸多因素的共同影响， 不同流路变化大、同一流路不同年代变化也很大。陆俭益、洪尚池等H1分析实测资料得到神仙沟流路年 均输往外海的沙量为9．62亿t、刁口河流路年均输往外海泥沙约4．07亿t；曹文洪b1通过地形变化对比 分析认为黄河口清水沟流路时期的海洋动力输送泥沙的作用是相对稳定的，年输沙量基本维持在2。0 亿一2．4亿t。 随着黄河人海水沙的变化和河口地形的演变，黄河121的潮流场和海洋动力输沙能力也必然会有新 的调整，需要有更深入的研究和探讨，为黄河口的综合治理提供科学依据。 目前，黄河口潮流过程和泥沙输移过程的认识主要来自实测资料分析。由于黄河口的实测资料非 常有限，且口f-j演变迅速，资料分析的结果往往很难全面了解黄河口的水流泥沙运动特征，难以满足生 产实践的需要。近年来，应用数学模型计算研究黄河12I的水沙运动的工作逐渐开展起来【6。】。然而，这 些研究大多只考虑M2分潮的作用，研究的对象仅限于清水沟单一流路。本文通过渤海大模型向黄河 口模型提供外海潮位边界，研究黄河口的水流运动特征，并分别计算研究了现行清8、清水沟原河道、北 汉流路、十八户流路、刁口河和马新河流路的泥沙通量和潮流输沙能力。图l为黄河口主要流路位置示 意图。 收稿日期：2007．11．05 基金项目：国家自然科学基金重点资助项目 作者简介：王崇浩(1968一)，男，浙江奉化人，博士。高级工程师，主要从事河流及河口海岸泥沙运动数模计算研究。 E—mail：ehwang@iwhr．c∞ 一1256～ 万方数据 图1黄河口主要流路 2计算原理 2．1基本方程采用基于MIKE21模型的数值模拟方法研究黄河口的潮流过程和海洋动力输沙能力。 MIKE21模型用于模拟垂向均匀的二维水流泥沙运动，已被广泛应用。以垂向积分后的物质守恒和 Navier．Stokes动量方程描述水位和水流的变化。 石方向动量方程： (1) Y方向动量方程： (2) 垂线积分后的水流连续方程： 蔷+羲+蔷2瓦薯+誊+碧=警