波流共同作用下的河口悬沙数值模拟

摘要

河口段的泥沙淤积阻碍通航，降低河道行洪能力，已经成为制约河口经济发展的主要因素之一。河口区域是淡水与海水的交汇处，在潮流与波浪的共同作用下河口区域的悬移质呈现季节性的冲淤变化，由气象条件及泥沙运移引起的河口淤积难以估计，需要建立一个稳定准确的泥沙运输数学模型。
　　 本文介绍了MIKE21模拟河口区域流场及粘性沙运移的基本原理。水动力模型基于水流连续性方程和不可压缩的雷诺时均Navier-Stokes方程；泥沙模型基于对流扩散方程，以波浪辐射应力来表示波浪对对水流和泥沙运移的影响，模拟波流共同作用下悬沙的运移，考虑了粘性沙的絮凝及干扰沉降，可进行多级粒径，多层底床冲淤模拟计算。
　　 通过搜集整理马颊河多年水文泥沙资料，应用MIKE21MT模块建立起马颊河河口波流共同作用下的平面二维泥沙数学模型，并以实测潮流及悬沙浓度资料进行模型率定。模拟了河口的流场分布及悬沙浓度分布变化，为满足河道防洪及上游供水要求，制定清淤、建闸、加高堤岸等河道治理方案，并根据模拟结果预测各方案的效果及治理后的河道淤积情况。
　　 由于对马颊河汇合口以下泥沙淤积的研究以前从未进行过，相关的历史淤积资料匮乏，以及受泥沙运移问题的复杂性及时间的限制，有些问题没有考虑，如河道突然开挖可能出现的骤淤问题，底床的固结，以及长期的冲淤变化情况，还有待深入研究。