长江口水动力学及其泥沙运输规律.docVIP

长江口水动力学及其泥沙运输规律 长江口概况： 长江河口地处我国东部沿海，受到来自流域径流、泥沙和外海潮流、成水入侵、风、波 浪及河口科氏力和复杂地形等绪多园了的影响，动力条件多变，泥沙输运复杂。从陆海相互 作用的角度看，长江河口至少存在几个水沙特性不同的典型河段，而每个典型河段又存在不 同性质的界面，如：大通河段(潮区界)、江阴河段(潮流界)、徐六径河段(盐水入侵界)、拦 门沙河段(涨落潮流优势转换界面)、口外海滨区(泥沙向海扩散的外边界和长江冲淡水扩散 的外边界)。每个典型河段及关键界而都涉及到物质和能量的传输；每个典型河段及关键界 面都有其固有性质，且相互影响，可以说河口过程在很大程度上被发生在每个典型河段的界 面上各种现象所制约。 水动力方程及验证 长江口水动力过程的研究进展（长江口水动力过程的研究进展） 在过去 20多年中 , 长江口水动力过程研究成果大量来自河口海岸学家、 物理海洋学家、 海岸工程师、环境流体力学家的文献、 著作。本文的目的是力图把这些文献(以正式发表的文献为准 ,不包括研究报告)汇集起来 ,对长江口潮流、 余流、 波浪、 盐水入侵的研究进行总结 , 究竟我们对长江口水动力过程了解多少 ?究竟长江口水动力过程还有哪些问题值得研究 ? 1.1 长江口余流、 环流、 水团、 长江冲淡水 基于现场实测资料 , 胡辉等 1985年对长江口外海滨余流的运动变化特性进行了一定的研究。 研究结果表明: 长江口外余流约为潮流的 1/ 2～1/ 5 , 上层余流以向东为主 , 中层余流多偏北 , 底层余流有偏西的趋势。径流是长江口外上层余流的重要组成部分 ,并以冲淡水的形式存在; 中、 下层余流则与台湾暖流的顶托和牵引有关。 王康　、 苏纪兰 1987年研究了长江口南港的横向环流、 垂直环流及其对悬移质输运的影响。 在前人基础上导出了长江口相对观测层次的物质断面传输公式 ,增加了反映环流及振荡切变的各种相互关系的有关项。 基于现场观测资料 ,Wang等 1990年研究了长江口水团、 长江冲淡水团等的基本特征。根据 1996年长江口南港水道枯、 丰水期大、 中、 小潮两次各 26 h的全潮水文实测资料 ,杨许侯等[1 ]统计分析了实测潮流的特征和潮流类型、 运动形式、 潮流垂直变化、 余流、 分潮对涨、 落潮流不等的影响。崔茂常1984年 , 张庆华等 1993年 , 朱建荣、沈焕庭 1997年对长江冲淡水进行了研究。 长江口潮流数学模拟 长江口水域开阔 , 口内多浅滩和沙岛 , 流场分布规律比较复杂。在长江口水域建立平面二维数学模型 , 有利于对长江口的水动力特性从宏观上加以研究。 同时 , 河口平面二维数学模型的差分求解方法已比较成熟 , 且运算量相对较小。 采用平面二维数学模型能够利用现有的计算设备 , 并在较短的运算时间内完成对长江口水域内流场平面分布特性的数学模拟。对平面二维数学模型的有限差分解法中 , ADI法由于其较好的稳定性且计算量相对较小 ,赵士清 1985年 ,刘上煊、 叶永 1987年 ,韩丕康、 黄国玲 1987年 ,刘世康、 徐建益 1987年 ,刘上煊、 陶学为 1987年 ,王船海、 程文辉 1991年 ,徐建益、 袁建中 1992年 ,许朋柱、毛锐 1993年 , 刘桦[2 ]等已将此法应用于长江口。此外 ,汪德　等 1987年用特征线法、 破开算子法等差分方法进行长江口平面二维潮流数学模型求解。 除有限差分法外 , 唐苓等 1992年将边界元法用于平面二维数学模型的求解;易家豪、 叶雪祥 1983年; 成安生等1987年采用有限单元法及限体积法等 , 赵士清 1985年采用三角形网格 , 并将局部的有限元法和有限差分法结合起来 , 在保证计算稳定性的同时 ,减少了计算量。参照国外学者的方法 , 刘上煊、 叶永 1987年将长江口平面二维潮流方程分别划分为: 动量平流项、水平动量扩散项以及波传播项。其中 ,动量平流项按特征曲线差分格式展开 , 水平动量扩散项按有限差分展开 , 波传播项采用隐式方向交替格式 (ADI ) 展开。 韩丕康、 黄国玲 1987年建立了长江口南支河段平面二维不稳定流计算的数学模型 , 着重长导堤整治工程情况下的数值解法。 从二维不稳定流的基本方程和流体连续方程出发 , 在加权剩余的 G alerkin 方法上 , 成安生等 1987年导出了长江口北支二维不稳定流的一种有限元数值计算模式。 长江口南槽各参量横向差异较小 , 潮流主要呈往复流形态。在该区域 , 横向积分的垂向二维模式显然是合理的。基于潮流有限差分方程 , 黄世昌等1994年选用 Prandtl 混合长度理论处理涡动扩散和黏滞系数 , 同时考虑由盐度差引起的密度异重流 , 采用Mu