小清河口海域悬沙输运格局数值研究.pdfVIP

第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 小清河口海域悬沙输运格局的数值研究 陈斌¨’2’3 黄海军4 严立文4 梅冰5 国人民解放军理工大学，南京211101．) 摘要：本研究构建了小清河附近海域三维悬浮泥沙输运的数值模型，并通过网格嵌套技术为 研究海域提供边界条件，其中水动力模型中主要考虑了潮汐潮流的作用，输运模型中考虑了 黏性泥沙的絮凝、再悬浮等过程。并与观测水位、潮流资料及悬浮泥沙观测资料进行对比。 从模拟结果来看悬沙含量总体趋势为由岸向海逐渐增大，并呈现出北高南低的分布特点，在 小清河入海口东北方向海域有一高值区。小清河入海口北侧水域活动性较强，且水流输沙能 力大于南侧海域。河口北侧海域有一较大侵蚀区，淤积区主要分布在河槽以南的滩涂和近岸 区域。 关键词：小清河口；泥沙输运；泥沙浓度 1研究背景 泥沙运动对人类从事近岸河口地区的开发和利用有着重要的影响，如港口建设、开挖航 道、修建防波堤、滩涂保护、河口治理、修建水坝等都必须先要了解周围海域的水动力条件 和泥沙运动情况。海岸河口地区主要受潮流、径流和波浪等水动力要素的影响，由此引起近 岸泥沙的冲刷和淤积，造成港口航道淤积等，倍受人们关注。随着计算技术的发展，数值模 拟可为解答这些问题提供了有力的手段。 小清河河口呈喇叭状，潮间浅滩宽阔，坡度平缓．近口门附近有牡蛎滩和拦门沙两段碍 航浅滩，其闻局部深槽逐年向东南摆动。长期以来，由于受拦门沙及河道摆动的困扰．羊口 港发展缓慢，与本地区经济发展的需要和优越的地理位置形成鲜明的反差。本研究根据大量 的实测资料和相关的数值结果，对小清河口附近海域泥沙运移规律、冲淤变化等问题进行了 深人的分析论证，为生活生产、港口扩建等提供了重要依据。 近几十年来诸多专家学者对河口的泥沙输运展开研究，如李本臣等通过大量实测资料的 比较，分析了河口海域泥沙来源、冲淤演变规律…；沈焕庭等乜1、李九发等【31基于观测资料 对长江口海域悬沙分布特征和输运规律做过研究，王家生等H1应用水槽试验研究了粘性泥沙的 运动特性及水流条件和泥沙浓度对粘性泥沙运动的影响；陆永军等瞄1针对渤海湾曹妃甸海域 第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 波浪、潮流、泥沙及海床演变特点，应用了波流共同作用下二维泥沙数学模型研究港区开发 方案。朱建荣等旧1应用改进的ECOM—SI模式，耦合泥沙输运模型，通过对理想河口数值试验， 探讨最大浑浊带形成的动力机制，并指出潮流掀沙是形成最大浑浊带泥沙的重要因素。本研 究构建三维水动力一悬沙输运模式盯川1，主要考虑潮汐潮流作用，模拟小清河口附近海域悬浮 泥沙浓度的分布情况，并与实测资料进行了比较，获得令人满意的结果。 2模型描述及配置 2．1模型的一般特征 sigIila坐标，以便于很好地拟合实际海域的地形；水平方向使用Arakawa C网格，并可采用正 交曲线网格以利于更好地贴合实际计算域的岸形。垂向湍黏性系数由湍封闭模型阳1计算，水平 湖泊等浅水海域的水动力模拟。SED模块是泥沙输运模块，可描述黏性泥沙和非黏性泥沙的沉 积、输运及再悬浮过程，其中对于黏性泥沙考虑了絮凝作用，对于非黏性泥沙输运、沉积及 再悬浮过程采用的是VanRijn方法¨0。。有关ECoMSED模式可参考文献[11]。 2．2悬浮泥沙输运扩散方程 ∞ J 氆 a)c a2 瓠、n瓠。 aZ、naz 等+警+等+塑誓墨=昙c4等，+品c4等，+昙c畅警，卸 a、}、it卸。 式(1)中(U巧彩)为三维流速；o为悬浮泥沙浓度：％为悬浮泥沙沉降速度；么何：为 水平扩散系数；翰：为垂向涡黏系数。 泥沙边界条件： ％孥：吣专，7 oZ K，车：EK—DK’z斗一H 韶 其中，风为再悬浮通量：风为沉积通量。 黏性泥沙再悬