长江口水文、泥沙过程与圆桩冲刷的数值模拟

本文主要利用ECOMSED模式对长江口及邻近海域的水文、泥沙过程进行三维数值模拟，并结合实测资料分析其水动力、泥沙输运、底床冲淤等特征；然后利用欧拉二相流模型模拟小尺度条件下长江口底床上圆桩周围的水流和泥沙冲刷、输运规律。 通过资料分析和ECOMSED数值模拟结果比较，我们得出：长江口口门内为非正规半日潮流区，潮流运动形式多为往复流，落潮流占优。落潮流速大于涨潮流速，流速垂向分布从表层到底层递减。悬沙浓度与流速关系密切，一般来说，流速越大，悬沙浓度越高；一个潮周期过程中会出现两次、三次或四次泥沙再悬浮，分别是涨急、落急、涨转落、落转涨时刻；盐水楔结构对粘性与非粘性悬沙浓度的分布起决定性作用，转流时泥沙再悬浮主要是由于这时会出现盐水楔，并形成垂向环流，使床面大量未被固结的泥沙再悬浮，形成峰值。悬浮泥沙垂向分布可分为垂线型，斜线型，抛物线型和L型。流场和底床冲淤变化与水深关系密切：深水区，流速较大，底床冲淤变化也较大。其中，受径流影响区表现为淤积，受潮流影响区表现为冲刷；浅水区基本表现为淤积。 从模式运行结果和实测资料比较可以看出，该模型可以较好的模拟长江口水流、悬浮泥沙分布与变化；能够再现在径流入海口处，盐水楔结构及其诱生的垂向环流从形成到发展，又到消失的完整过程；也能够展示底床的冲淤变化。对于我们模拟长江口背景流场，了解该区域内水动力变化、悬浮泥沙输运、底床冲淤等有重要意义。 在欧拉二相流模型对长江口底床上圆桩周围的水流和局部冲刷数值模拟过程中，我们不仅考虑水流和泥沙之间的作用，还引入泥沙颗粒之间的相互影响。模拟结果较合理的展示了圆桩周围的流场类型和底床冲刷变化：在圆桩前方，流速减小并形成垂向涡旋，从而产生局部冲刷；在圆桩两侧，水流加速，挟带上游泥沙向下输送，并在内侧堆积；而在圆桩后面，存在流速分离区。在该分离区内流速很小，并且当流速较大时，会产生回流，形成两个对称的漩涡。流速越大，圆桩前由垂向涡旋引起的局部冲刷就越明显；而当底床泥沙粒径变小时，泥沙临界起动流速变小，底床也更容易被冲刷。悬浮泥沙浓度分布受流场的影响，并且当粒径小而流速大时，能悬浮到更高的深度。