平面二维泥沙输移模型及其应用

第 11 卷 2011 年
第8期 8月
中 国 水 运 Chi na W er Tr a ns por t at
V . 11 ol Augus t
N 8 o. 2011
平面二维泥沙输移模型及其应用
夏雪瑾，高程程
（上海市水务规划设计研究院，上海 200232） 摘 要：为研究海岸工程对周边水沙环境和河床冲淤的影响，运用 Mik e2 1 建立了平面二维水沙数值模型， 并将模型应用
工程实例。通过实测资料验证流场和含沙量场，结果表明模型能较好地模拟含沙量场。在此基础上，从工程实施后引起的 水流变化以及年冲淤变化强度这些指标来判断工程实施对海域的影响。研究成果可为今后类似的工程提供一定的参考。 关键词：二维悬沙模型；数值模拟；含沙量；泥沙回淤 中图分类号：TV 14 文献标识码：A 文章编号：1006- 7973（2011）08- 0082- 03 为零；开边界条件（水边界条件） s (x , y, t ) = s (x, y, t ) ： 二、应用实例 应用该 模型，建立舟山群岛海域二 维泥沙数学模型，计 算分析工 程给周边海域泥沙带来变化 ，为工程的决策提供技 术依据。 1．模型范围 考虑计 算水边界取在基本上不受本 工程影响的海域，同 时兼顾到 水文条件等相关资料获取的 方便，选取模型（右图 1 ）开边界的北边界为 3 0 °19 ′00 ″N，东边界 1 22 °4 0 ′0 0″E ，南边界 2 9°2 6′0 0″N，西边界 1 21 °37 ′0 0 ″E ； 模型闭边界为自然岸线， 包括象山港全域， 最西在 12 1 °2 5′2 4 ″E。模型范围约 1 20 km ×10 0 km ，采用三角形 网格， 共布有 46 43 3 个三角形单元， 域内最大水深达 1 00 m 以上，最小空间步长约为 35 m 。
在河口海岸地 区兴建工程时，特别是兴 建重要的航运或 围垦工程时， 往往要求评估工程实施对周 边水域水沙环境的 影响。河口海 岸泥沙运动规律研究比较复 杂，而在规划设计 阶段，通常要 求迅速给出比较概括性的解 答。因此，数学模 型研究泥沙运 动被越来越广泛地应用于解 决近岸工程实际问 题。水流泥沙 数学模型可对岸滩的演变、 海床冲淤等提供长 期预报，为工程的规划设计提供科学依据。 目前在悬沙、 底沙输移以及河床演变研 究中，二维泥沙 数值模型应用 最为广泛。国内外相继出现 了一批功能强大， 通用性好、 成熟的商业化综合数学模型， TRIM 2D 模型[1 ]、 如 美国的 Miss ips ipp i 大学水科学计算中心 CC HE 2D 模型 、 丹麦水利所的 M iKE2 1 模型等， 大多数情况下具有足够的精 度能满足工程要求，被广泛应用。本文主要利用 Mike 2 1 建 立二维平面水 沙数值模型研究工程对大范 围水域泥沙运动的 影响，为今后类似工程的规划设计提供一些借鉴。 一、潮流泥沙数学模型 1．基本方程 基本方程包括二维浅水方程（此处略）和二维悬沙方程。 悬沙基本方程：
S t + u S x + v S y = 1 h x h Dx S x + 1 h y hD y S y + Fs h
[2]
图 1 模型范围示意图 2．边界条件 本计 算域潮位开 边界由东中 国海潮波数学 模型提供[3 ]。 悬沙模块 由于只有同步的实时实测含 沙量资料，整个边界取 所有实测含沙量的平均值 0.2 5k g/ m 3 作为悬沙的边界条件。 3．计算参数的选取 （1 ）糙率 n ：由于模型计算域较大，根据模型研究区域 实际床面情况通过试算作适当调整，n 取值在 0 .02 5~ 0.0 4 之间。 （2 ）沉速：工程区域悬移质平均粒径为 0 .02 9m m ，易 受 海水 影响 发生 絮凝而 加速 沉降 ，因 此沉 速计 算公式 为： 。 ω = ω FD （F 为絮凝因子；D 为衰减系数） 0 （3 ）泥沙临界淤积切应力： 经过反复调试，选取临界淤 积切应力为 0 .1 3N/ m 2。 （4 ）泥沙临界冲刷切 应力τ 经过反复调试，选取泥 e ： 沙临界冲刷切应力在 0.1 5- 0.2 5N/ m 2 之间。
其中：S 为悬沙浓度，D x、D y 分别为 x、y 方向上的泥 沙扩散系数；F s 为泥沙冲淤函数。 2．泥沙冲淤函数 模型中泥沙冲 淤函数采用的是切应力法 ，由床面临界淤 积切应力和临界冲刷切应力确定源汇项： ω cb (1 ττ ) s d
Fs = 0 E (1 ττ)n e
τ< τ d τ < τ< τ d e τ≥ e τ
式中： 为瞬时底床剪切应力， d 为临界淤积切应力， e τ τ τ 为临界冲刷切应力， 为床面泥沙冲刷系数， E 由率定计算确定。 3．定解条件
初始条件：整 个计算域内每一个节点的 水位和流速、流
向由流场模型计算结果提供， 悬沙浓度初始值在开始时取零。 边界条件：闭 边界（即陆地边界）取含 沙量的法向梯度 收稿日期：2 01 1- 0 6 - 14 作者简介：夏雪瑾，上海市水务规划设计研究院。

第8期 4．模型验证
夏雪瑾等：平面二维泥沙输移模型及其应用
83
域的泥沙 运动规律，因此本文中采用 河口挟沙力公式。根据 实测资料相关性分析拟合出来的河口挟沙力公式为： 大潮：S = 64.845(
v2 gh )0.4388
2 v2 0 .80 62 ) ； 中潮：S = 87.191( v ) 0.85 97 ； 小潮： gh gh
模型采用了 2 0 07 年地形资料和工程区附近水域的大、 中、小潮水文 实测资料进行验证，潮位验 证取定海、岱山、 桃花岛、沈家门、镇海共 5 个潮位站的逐时实测潮位，流速 和含沙量验证取用 1 1 条水文测线的实测值。 （1 ）潮位和流速验证 根据潮流模型验证结果（由于篇幅所限，验证图略） ：潮 位、 流速流向的计算值与实测值吻合良好， 符合相关规程要求； 工程区测流断面的涨落潮量计算值与实测值较为一致，很好地 反映了工程区域进出潮量的实际情况； 模型能够较好地反映工 程区的水流特性，并为悬沙模型提供较为准确的水动力条件。 （2 ）含沙量验证 泥沙数学模型 的困难之处在于泥沙运动 本身的复杂性和 验证资料的匮 乏。本文从两方面来对含沙 量进行验证：①对 有实测资料的 1 1 个测站的含沙量过程线进行局部验证；② 整体含沙量分布趋势要能反映工程区海域泥沙运动的趋势。 1 ）含沙量过程线验证
含沙量 （ kg/m 3） 0 4 . 0 .3 5 0 3 . 0 .2 5 0 2 . 0 .1 5 0 1 . 0 .0 5 时间 0 1 1 0 :0 0 1 2 0 :0 0 1 3 0 :0 0 1 4 0 :0 0 1 5 0 :0 0 1 6 0 :0 0 17 0 00 : 18 0 00 : 19 0 00 : 20 0 00 : 2 1 0 0 : 0
S = 2.34(
式中：s 为全潮平均挟沙力，v 为全潮平均流速，h 为全 潮平均水深。 比较挟 沙力公式和悬沙模型计算的 工程海域的全潮平均 含沙量场分布图（图 4 ，以中潮为例） ，可见挟沙力公式与悬 沙模型计 算的全潮平均含沙量场分布 相当相似，但公式计算 的结果量 级上偏大一些。主要是因为 资料在测验期间有偏北 风，风浪 大于年平均情况，因此采用 此次实测数据来拟合挟 沙力公式 也隐含考虑了波浪影响因素 ，同时也表明波浪对当 地海域含 沙量实际分布有一定影响。 总的来说，群岛东南向 外部海域 的含沙量相当小，而靠近杭 州湾的海域含沙量则远 大于舟山 群岛内部水域及东南方向的 群岛外部水域，这种分 布情况与 实际情况是吻合的，各个海 域含沙量的量级也与实 测资料中含沙量的量级及分布相一致。
P2
计 算值 实 测值
图 2 p2 含沙 量验证曲线 结果（代表测点 p 2 验证见图 2）表明，模型泥沙参数选 取比较合理， 计算值基本可以反映含沙量 随潮变化过程，涨 落潮最大误差、平均误差基本控制在 50 %以内，符合有关规 范规程的要求，基本反映了该海域的悬沙输运特征。 2 ）泥沙运动趋势验证 由于本区缺乏 整个工程研究区域完整的 、长历时的泥沙 实测资料，而 且前人对该区研究成果也较 少，因此，本文中 从国内常用水 流挟沙能力公式（刘家驹公 式、窦国仁公式、 张瑞瑾公式和 河口挟沙力公式）中选用公 式来计算得出大中 小潮平均含沙 量场，然后与用模型模拟的 大中小潮平均含沙 量场作比较，以做校核。 因为缺乏波浪 资料，因此没有选用刘家 驹公式和窦国仁 公式，主要在 张瑞瑾公式和河口挟沙力公 式之间考虑。根据 实测大中小潮 数据并结合这两种公式对流 速与含沙量的相关 性进行拟合分析，相关分析拟合线（以小潮为例）见图 3 。
0
（a）拟合公式
（b）悬沙模型
图 4 中潮全潮平均含沙量 综上分 析，模型计算的含沙量场不 仅与实测资料吻合较 好，而且 能较好反映工程区泥沙运动 整体趋势，该模型可用 来进一步研究工程区海域泥沙淤积问题。 4．工程对附近海域的影响 在工程区海域平面二维水沙模型验证较好的基础上，应用 模型对工程（右图 5 ）后水流、泥沙进行计算，并采用涨落流 速的变化、 海床冲淤变化等指标分析工程实施后对海域的影响。
0
- 1
-1
- 2
-2 实测
- 3
实测 相关线 -3
相关 性
图5
对附近海域水流流场的影响
- 4
- 5 -6 -5 - 4 -3 -2 - 1 0
-4 -1 0 -8 -6 -4 -2 0
（1 ）对附近海域水流流场的影响 根据计 算，由于工程后截断了登步 岛与桃花岛间的潮通 道，使得 工程区近旁水域流速减小， 主要的流速减弱区为工 程区西北 和东南水域，此两水域为涨 、落潮流的流影区。而 工程两侧 进潮通道流速均有所增大， 其中朱家尖与登步岛间 的流速增大较多 ，最大增大 0 .1m / s ，这种状况将更有利于 航道水深 的维护。但是，螺头水道的 流速略有减少，主要是 由于总进潮量的减少导致，需加以注意。
图 3 小潮流速与含沙量相关性分析图（张瑞瑾公式左，河口 挟沙力公式右） 从相关系数 R 可以看出，张瑞瑾公式拟合相关系数太低 在 0 .3 以下， 而河口挟沙力公式拟合相关系数都在 0 .6 以 上， 表明全潮平均含沙量 S 与全潮平均流速 v 以及全潮平均水深
h 之间存 在相关的函数关系，公式能较好的拟合出工程区海