温州浅滩围涂促淤工程泥沙数值模拟的研究.pdfVIP

泥 沙 研 究 年 月 第 期 温州浅滩围涂促淤工程泥沙数值模拟研究 李孟国 交通部天津水运工程科学研究所 天津 摘要 建立了波浪潮流泥沙数学模型 对瓯江口温州浅滩围涂促淤工程效果进行了数值模拟研究 通过多个 促淤方案的淤积效果模拟计算 推荐了最佳促淤方案 推荐的最佳促淤方案与定床浑水淤积物理模型试验推 荐的方案一致 关键词 促淤 泥沙 瓯江口 温州浅滩 数值模拟 中图分类号 文献标识码 文章编号 前言 温州浅滩是瓯江口外规模最大 发育最完善的拦门沙滩 位于灵昆岛与霓屿岛之间 图 温州浅 滩围涂工程就是沿浅滩南北两侧建两条连接灵昆岛和霓屿的海上大堤 北堤顺浅滩 等深线 南堤为 灵昆岛 小霓屿 霓屿堤 将浅滩人工围垦变为陆地 图 目前北堤建设已经接近完成南堤尚未建设 图 瓯江口附近地形示意图 温州浅滩围涂造陆面积达 万亩 南 北堤建设完成后 要使中间围涂区域成为陆地将需要大量的 收稿日期 基金项目 中国博士后科学基金 编号 资助 作者简介 李孟国 男 天津市人 研究员 博士 主要从事海岸河口水动力与泥沙的数值模拟研究 土石方量 如果能结合南堤的建设 使泥沙尽可能大量地在围涂区域内自然落淤将会节省大量的人力 物力和财力 因此有必要对泥沙促淤方案形式进行研究 更好地借助自然之力完成部分围涂 本文拟建立平面二维波浪潮流泥沙数学模型 采用数值模拟研究手段 对温州浅滩最佳促淤方案形 式进行研究 在河口海岸地区 引起泥沙运动的主要动力因素是潮流和波浪 波浪掀沙 潮流输沙 被公认为是 泥沙运动机制 目前考虑波浪潮流共同作用的泥沙数学模型被广泛采用 在目前的考虑波浪潮流 共同作用的泥沙数学模型中 通常的做法是通过波浪辐射应力 波流底切应力 波流挟沙力来考虑 波浪 掀沙 潮流输沙 因此 该类泥沙模型涉及 四场 的计算 即波浪场 潮流场 悬沙场和底床冲淤场 本 文建立了综合考虑多种因素的多向不规则波波浪场数学模型 考虑波浪作用的二维潮流场数学模型 二 维悬沙场数学模型 悬沙底沙造成的底床冲淤场数学模型 波浪数学模型采用矩形网格 潮流 悬沙及冲 淤数学模型采用三角形网格 网格节点处的波浪要素采用插值法从波浪模型的矩形网格得到 基于现 场资料得出了波流挟沙力公式 在充分验证的基础上 对瓯江口海区温州浅滩围涂促淤方案进行了数值 模拟研究 数学模型的建立 波浪场数学模型 基于缓坡方程和线性波动叠加原理导出的多向不规则波联合折绕射基本方程为 式中 表征频率间隔 及方向间隔 内组成波的能量 分别为组成波的位相函数梯度 波数 波向 波速及群速 采用矩形网格有限差分方法对式 进行数值求解 关于差分方程及对底摩擦波能损失 波浪破碎及岛屿影响的考虑参见文献 波浪计算域的西边界至 经度线 东边界至 经度线 北边界到乐清湾顶 南边界至 纬度线 东西距离 南北距离 用正方形网格剖分计算域 网格间距 网格数 个 以南麂岛波要素作为模型推算的验证条件 考虑波浪作用的潮流场数学模型 基本方程 连续方程 动量方程 其中 为直角坐标系坐标 分别为 方向的流速分量 为重力加速度 为海水密度 分别为 方向的水平涡动粘性系数 为时间 为柯氏参数 为相对 平面的水面起伏 为静 水深 为辐射应力张量的四个分量 分别为波流底摩擦应力在 方向的分量 当不考虑波浪作用时 方