1. 引言
悬浮物指悬浮在水中的固体物质,包括不溶于水中的无机物、有机物及泥砂、黏土、微生物等常常悬浮在水流之中,使水产生的浑浊现象,其含量与分布情况决定了水域环境的优劣 [1] 。因此研究悬浮物的输送过程对确定海区环境质量和环境治理有着重要意义。
目前,国内外学者已对悬浮物输送问题进行了一定的研究。一些学者通过对相关水域的海流和悬浮物进行实测得到悬浮物输送量结果,如Qiu等 [2] 分两个时间段,采用定点观测法对黄河口水流及颗粒浓度进行了实测,估算了悬浮物的下行通量,这一方法操作简单,准确度高,但样本数据的时间连续性和覆盖范围有限,适合于小范围水域;另一种则是采用卫星遥感与实测相结合的方法,尹小青等 [3] 通过建立MODIS卫星影像数据与现场实测数据的对应关系,利用海流速度反演得到悬浮物的输送量,这种方法避免了大面积采集样本,但受天气条件和卫星分辨率影响较大,无法获得任意时刻的悬浮物运动数据。
鉴于以上方法的优势和局限性,计算机数值模拟渐渐成为新的研究热点,Cai [4] [5] 通过引入风浪潮流和泥沙运动的MIKE21模型以及SWAN模型,建立了一种能够模拟风浪、潮流和泥沙的数值模型;Yang等 [6] 基于sedtrans05悬浮物输运模型和破波带内的泥沙输运公式,建立一个计算推移质输运的悬浮物输运模型。以上研究通过建立悬浮物运动模型,研究特定水域的悬浮物输送率,但模型计算结果未考虑初始条件、颗粒、水深、温度、海床高等因素的影响,适用范围有限。
针对以上问题,本文选取近岸河口悬浮物中悬浮泥沙为研究对象,考虑水的密度和粘度、沉降流速,以及传输方程,搭建了基于sedtran05垂直剖面输沙模型的计算方法;采用模型对黄海环境进行模拟,计算了不同泥沙粒径、水深、流速和温度时的悬浮物输送率,并对模型相关参数进行了修正;通过将仿真结果与部分文献数据进行对比,验证了方法的有效性并明确了模型的适用范围。
2. 垂直剖面输沙模型
Sedtrans模型由Nova Scotia于20世纪80年代提出 [7] ,能够用于模拟稳定流或在海浪、海流综合影响下,沙子或粘性悬浮物的输运问题 [8] 。
基于Sedtrans05建立的数值模型(垂直剖面输沙模型),在计算粘性悬浮物输运量时,对算法做出改进,计算中考虑波、流共同作用效果,增加了影响输送率的参数计算过程,如水的密度、粘度、温度、悬浮物颗粒大小、沉降流速、再悬浮临界速度、传输方程等,能够逼真地模拟洋流作用下的悬浮物输送过程,而不受模型计算效率和复杂性的限制,进而计算悬浮物输送量。
2.1. 水的密度和粘度
水的密度和动力学粘度,与海水的温度和盐度有关,如式(1)所示
(1)
式中计算误差在温度范围0到
28 摄氏度
时小于1%。
2.2. 沉降流速
采用Gibbs公式来计算球形颗粒的沉降速度 [9]
(2)
其中v是
的运动学粘度,D为中值粒径,
为无量纲。
2.3. 启动再悬浮的临界剪切速度
泥沙悬浮的临界剪切速度
,如式(3)所示
(3)
2.4. 传输方程
泥沙传输速率与粒子速度的速率,跃移高度和泥沙浓度有关,瞬时泥沙传输速率和剪切应力参数有关,泥沙在纯净状态下的传输速率为 [10] :
(4)
其中s是悬浮物和水的密度的比率,
是一个常数等于0.053。在这里,
是瞬间肤面摩擦流切应力,而
是启动泥沙运动的临界剪切应力。
合并流和波的情况下瞬时推移质输沙率有:
(5)
上式子中
,是个校准系数,
为显著波高度(m),H为水深(m),
是瞬时结合皮肤摩擦剪切应力,时间平均推移质输沙率由一个波的周期内的平均值得到。
3. 仿真实验
实验主要以黄海海洋环境为背景,研究粒径大小、水深、流速、和温度等海洋参数,对悬浮物输送率的影响。黄海北部东岸的悬浮物为黏土质粉砂,其密度为1.469 g/cm3。将参数分别设置为水深20 m,流速为0.42 m,流向为从北为正顺时针150˚,需要计算的位置为海床上方5 m,波高1.5 m,波的周期设为3.0 s,波的方向为以北顺时针145˚,泥沙粒径设为0.025 mm,盐度30,温度为25℃。图1为sedtrans05模型中五种不同方程所计算出的悬浮物输送率。
3.1. 泥沙粒径的影响(0.01~0.075 mm)
图2为悬浮物输送率与泥沙粒径的变化关系,由图中可知,悬浮物输送率在泥沙粒径为0.062 mm时出现突变,因此可以将图像分为两种情况进行分析,分别是粒径大于0.062 mm和小于0.062 mm。可以得到,悬浮物的输送率会随着粒径的增大而减小,且由图上显示为凹函数可知减小的趋势不断增大。
![](//html.hanspub.org/file/14-1541384x24_hanspub.png)
Figure 1. Result of conveying sand with different output equation in sedtran05
图1. sedtran05输出不同方程计算输沙结果
![](//html.hanspub.org/file/14-1541384x25_hanspub.png)
Figure 2. The relationship between the transport rate and sediment particle size
图2. 悬浮物输送率与泥沙粒径的变化关系
3.2. 水深的影响(5~30 m)
图3为悬浮物输送率与水深变化的关系,对于悬浮物输送率,由图中可知,在水深小于13 m时,输送率随深度增加而减小,而当水深超过13 m后,输送率受水深影响较小,这一临界值区别了浅水和深水的悬浮物输送情况。
3.3. 水流流速的影响(0.1~1 m/s)
图4为悬浮物输送率与流速的变化关系,有图中可知,悬浮物输送率会随着水流速的加快而增大,且为增大速率会不断变大,表明海水流速对悬浮物的输送起着重要作用,验证了波浪掀沙,海流输沙这水深一规律。
![](//html.hanspub.org/file/14-1541384x26_hanspub.png)
Figure 3. The relationship between the transport rate and depth of water
图3. 悬浮物输送率与水深变化的关系
![](//html.hanspub.org/file/14-1541384x27_hanspub.png)
Figure 4. The relationship between the transport rate and velocity
图4. 悬浮物输送率与流速的变化关系
3.4. 海床上方高度的影响(0.1~20 m)
图5为海床上方高度与输沙率的关系,有图中可知,当某位置越接近海床时,输送率越高,反之,越接近海面时,输送率越低。
![](//html.hanspub.org/file/14-1541384x28_hanspub.png)
Figure 5. The relationship between the transport rate and sea level
图5. 悬浮物输送率与海床上方高度变化的关系
3.5. 温度的影响(10℃~25℃)
图6为悬浮物输送率与温度的关系,由图中可知,当温度升高时会降低悬浮物的输送率,且图中曲线为凸函数,表明降低速率会随着温度的升高略微减小,即输送率整体减小,但减小的趋势逐渐缓慢。
![](//html.hanspub.org/file/14-1541384x29_hanspub.png)
Figure 6. The relationship between the transport rate and temperature
图6. 悬浮物输送率与温度的关系
由以上仿真分析可知,在不同的泥沙颗粒大小、水深、水流速度、海床高度等条件下,输沙量具有较大差别,且结果表明,Sedtran05模型所做的是推移质数值模拟,不适合于中国黄海等泥沙粒径较小地区,适合于泥沙粒径大,并且沿岸海水较深的地方,如中国南海,海州湾,东海陆架海地区等。
4. 结语
本文选取近岸河口悬浮物中悬浮泥沙为研究对象,首先,考虑水的密度和粘度、沉降流速,以及传输方程,搭建了垂直剖面输沙模型的计算方法;然后,采用sedtran05模型对黄海环境进行模拟,计算了不同泥沙粒径、水深、流速和温度时的悬浮物输送率,并对模型相关参数进行了修正;最后,通过将仿真结果与文献数据进行对比,验证了方法的有效性并明确了模型的适用范围。
*通讯作者。