水利工程论文-长江口细颗粒悬沙浓度垂向分布.doc

水利工程论文-长江口细颗粒悬沙浓度垂向分布摘要：在长江口利用“声学悬沙观测系统”观测到大潮四种典型的高时空分辨率细颗粒悬沙浓度声学垂向分布图和垂线分布曲线。悬沙浓度垂线分布曲线表明：(1)在涨潮时呈L形，悬沙浓度的垂向变化梯度小(浓度小于1.0g/l)。(2)在涨急时呈射流形，射流顶的悬沙浓度达10g/l，悬沙浓度的垂向变化梯度大。(3)在落潮时从水面到水底悬沙浓度按指数增加，可能代表恒定均匀流中悬沙处于平衡条件的分布，泥沙垂向扩散系数S在数值上大约是泥沙颗粒沉降速度的两倍。(4)在落急时是不连续，并且在水体中部呈射流形，悬沙浓度的垂向变化梯度大。关键词：长江口细颗粒悬沙浓度高时空分辨率声散射1前言Kirby和Parker1首次用LUTOCLINE(泥跃层)一词描述潮汐河口粘性细颗粒悬沙浓度垂向分布的不连续性，即垂向悬沙浓度出现明显的变化。此后，泥跃层和近底高含沙层被发现普遍存在于许多高浑浊的河口海湾25。Ross和Mehta6、Smith和Kirby7以及Wolanski等59详细地研究了泥跃层和近底高含沙层的动力特性，并建立了相应的数学模型。这些研究加深了我们对潮汐河口近底水泥相互作用的认识，这种相互作用控制整个潮汐河口细颗粒泥沙的输移过程。在海岸工程应用上，如河口航道、港口、码头和系船池，需要对细颗粒泥沙运动定量化。由于河口悬沙浓度随深度迅速变化，因此，悬沙浓度的精确观测是一个关键。至今，为了获得近瞬时的悬沙浓度的垂向分布，已发展了许多新的技术。特别是声学技术被广泛应用于海洋中悬沙浓度的测量1013，使我们能得到高时空分辨率悬沙浓度垂向分布；但是声学技术尚未广泛应用于河口粘性细颗粒悬沙浓度的连续观测。国内，在河口细颗粒泥沙输移现场观测研究中，悬沙浓度资料大多依赖于野外水样的采集1416。Li等17利用光学测沙系统研究了浙江椒江河口细颗粒泥沙输移。本文作为“声学悬沙观测系统”在长江口北槽口外悬沙运动研究应用的系列成果1820，其目的是为了解河口细颗粒悬沙浓度近瞬时垂向分布提供一种较好的途径。2研究的区域和方法长江口是一个中潮河口，平均潮差是2.0m，悬沙年平均达4.9108t。许多科研工作者对长江口悬沙输移问题进行了研究15，16，21，悬沙浓度资料大多依赖于野外水样的采集。在长江口的陆架地区水面悬沙浓度超过0.02g/l，在水底超过0.1g/l22。我们的现场观测是于1994年5月2829日大潮时在长江口9405站(东经122274，北纬310515)(图1)进行了28h连续观测(图2)。观测数据的采集工作是在东海舰队东调221军舰上完成的。用于测量悬沙浓度垂向分布的ASSM型声学悬沙观测系统是由中国科学院东海研究站研制的，与Libicki等12叙述的观测系统相类似。整个系统(图3)在PC计算机控制下采集数据，并进行预处理，再储存。ASSM型系统发射声脉冲频率为0.5MHz，脉宽为40s。从悬沙背向声散射测得的垂向分布其垂直(深度)分辨率为10cm，横向(时间)分辨率为1.5s。每小时记录12min数据，其采样频率为75KHz。每幅图是由480条海底和换能器之间悬沙的背向声散射曲线所组成。由泥沙背向声散射信号的强度转化为悬沙浓度的定标工作是在有机玻璃沉降圆筒内完成的。基于实验室定标，绝对悬沙浓度测量的误差约为16%。ASSM系统换能器工作时通常放置在水下2m处。同时，利用流速仪测量水流的流速和流向。每小时分别在0.0D、0.2D、0.4D、0.6D、0.8D、1.0D六个点采水样(D为水深)，此法通称“六点法”。由采水器采集到的含悬沙水样，通过63微米滤纸，分离出粗的悬浮颗粒。通过过滤的水样留下来的悬浮物质用蒸馏水漂洗，并在40恒温箱内烘干，过滤器烘到60，再称重量。悬沙浓度是从最后沉淀物重量(包括粗颗粒)与被滤去总的水体积而换算出来的。图1长江口9405站位Locationofstation9405intheChangjiangEstuary图2潮位和垂线平均悬沙浓度过程图图3实验布设示意图Hydrographsoftidallevelanddepthaveragedsuspensionconcentrationhourlysampledatstation9405Sketchofexperiment3结果在上述观察区域内，盛行半日大潮。垂线平均悬沙浓度随潮水位的变化而变化(图2)。平均悬沙浓度从0.45g/l变化到2.41g/l。图4(ad)是每15min的悬沙浓度垂向分布的声学灰度图。图4(a)表明整个水体中均有悬沙分布；图4(b)中、下部有一高含沙层；图4(c)中70%的水体有较低的悬沙浓度；在图4(d)中的水体中部有较高的悬沙浓度。为进一步了解悬沙浓度垂线分布特征，将图4(ad)中每一个悬沙浓度声学图象等分成7段，悬沙浓度数值是按10cm水深间隔绘图，由此可得到7条瞬时悬沙浓度的垂线分布曲线(图5)。图5(ad)揭示了大潮四种典型的悬沙浓度的垂线分布曲线。第类是在涨潮时呈L形(图5(a)，近底悬沙浓度大于2.5g/l，大部分水体中悬沙浓度小于1.0g/l。平均悬沙浓度随水深从0.43g/l变化到0.84g/l。第类是在涨急时呈射流形(图5(b)，近底悬沙浓度大于2.0g/l，悬沙浓度的垂向梯度大。第类是在落潮时从水面到水底悬沙浓度按指数增加(图5(c)。近底悬沙浓度大于6.0g/l。平均悬沙浓度随水深从0.86g/l变化到0.96g/l。第类是在落急时，悬沙浓度在水体中部呈射流形(图5(d)。近底悬沙浓度大于5.0g/l，悬沙浓度的垂向变化梯度大，平均悬沙浓度从0.77g/l变化到1.24g/l。事实上，图5(ad)分别代表在涨潮、涨急、落潮和落急时细颗粒泥沙再悬浮(resuspension)过程的图像。图4声学灰度图Grayscaleplotsat(a)floodtide(1000h)；(b)max.floodtide(1200h)；(c)ebbtide(1500h)；(d)max.ebbtide(1800h)；28May1994.图5悬沙浓度垂线分布Verticalprofilesofburstaveragedsuspensionconcentrationover12min.periodat(a)floodtide(1000h)；(b)max.floodtide(1200h)；(c)ebbtide(1500h)；(d)max.ebbtide(1800h).28May1994.图6表示在涨潮(1000h)时最大流速为1.17m/s，流速和盐度垂向变化梯度均较大。在涨急时流速和盐度垂向变化梯度小。落潮(1500h)时最大流速为1.03m/s，盐度垂线分布表明底部有一盐水楔，水底盐度为26.0。在落急(1800h)时，流速和盐度的垂向变化梯度渐缓。4讨论4.1悬沙浓度垂向分层结构图4和图5中最明显的特征是悬沙浓度垂向呈分层结构，形成这种特征主要是由于紊流能量引起强烈而均匀混合，从而产生陡的密度梯度23。图6流速、盐度的垂线分布Verticalprofilesofvelocityandsalinityat(a)floodtide(1000h)；(b)max.floodtide(1200h)；(c)ebbtide(1500h)；(d)max.ebbtide(1800h),28May1994.一个稳定的悬沙浓度分层结构意味着当河口保持盐度和温度充分混合时悬沙有强的分层现象3。图4(a)和图5(a)表示在涨潮时，由于物质扩散率较低形成稳定的泥跃层和在上面的水域有较低的泥沙浓度。图4(b)和图5(b)表示在涨急时形成一个浮泥层和二个稳定的泥跃层。图4(c)和图5(c)表示在落潮时形成一个浮泥层(高含沙层)，一个持续的泥跃层以及在上面水域很低的泥沙浓度。图4(d)和图5(d)表示在落急时迅速形成一个流动的浮泥层(高含沙层)和二个起伏不定的泥沙跃层(分别在浮泥层上、下)。图5(d)的瞬时悬沙浓度垂线分布表明有二个泥跃层。Mehta4认为下泥跃层主要是由于泥沙沉降受阻引起的，上泥跃层形成与向上的紊流扩散的非线性有很大关系。上泥跃层和下泥跃层二者之间是最小的垂向混合区域4。Wolanski等8详细地研究了泥跃层混合的问题，他们发现泥沙混合受浮力效应和在高浑浊底层里泥沙引起的紊流破裂二者控制，悬沙受到的浮力阻碍泥沙垂向扰动混合，扰动能够缓和悬沙的沉降速度，并且在扰动方式上有一反馈效应9。4.2悬沙平衡在河口最大浑浊带，如不考虑平流影响，悬沙扩散方程可简单表示为(1)式中C是平均悬沙浓度；t是时间；S是泥沙颗粒沉降速度；z是距底床的高度；S是紊动扩散系数。事实上，图4中悬沙浓度的垂向分布表示了悬沙向下沉降和向上紊动扩散二者之间的平衡。在涨急和落急时，悬沙强烈地输移，泥沙沉降少。在1500h，水体上部、中部悬沙浓度均低，近底部悬沙浓度高(图4(c)和图5(c)；而且图5(c)反映了悬沙浓度随水深按指数增加。这可能是由于底床切应力使泥沙再悬浮和泥沙颗粒沉降保持平衡所致，类似的特征在文献里有过报道24。这种悬沙浓度垂向分布代表恒定均匀流平衡条件，即垂线上泥沙上下交换质量守恒，并可用以下方程表示C=Cexp(z)(2)式中是与半对数图的梯度成正比。公式(1)中对距离z及时间t的偏微分项均等于零，恒定均匀平衡条件下泥沙紊动扩散方程可表达为SC+S(dC/dz)=0(3)式中S为垂向泥沙扩散系数，它是随水深恒定的，并具有值S=S/(4)根据公式(2)计算出的值平均为0.58m-1，因此，由公式(4)可知，泥沙垂向扩散系数S在数值上大约是泥沙颗粒沉降速度S的两倍。图6中流速垂线分布也表明水体上、中部流速变化梯度较小，可视为恒定均匀流；而近底高含沙层对应着盐水楔异重流(图6)。4.3再悬浮过程图5(ad)分别代表在涨潮、涨急、落潮和落急时细颗粒泥沙再悬浮(resuspension)过程的图像。这种泥沙的再悬浮过程是受近底紊流大小和速度密度分布结构的控制23。在长江口，由风激起的表面波，潮流、内波破裂和江苏沿岸流(向南)是造成泥沙再悬浮过程的原因。近底观测到大的流速梯度，从而可以认为底部存在大的切应力(图6)，这些海底的切应力是泥沙再悬浮的动力。波浪诱导的河口淤泥再悬浮问题有人做了水槽试验研究25以及现场观测研究26。在实验室水槽研究基础上，Maa和Mehta25认为波引起的底床侵蚀是产生浮泥的原因。由波侵蚀的泥沙保留在近底处，占整个水深的20%，水流仅在水平方向迁移。在长江口，波浪和潮流的特性还需要进一步研究，由此可以估价它们在泥沙再悬浮中的作用和重要性。Graber等27在中国东海研究了风暴产生的表面波和泥沙的再悬浮问题。可以相信由风暴产生的波是使泥沙再悬浮的重要动力。5结论在长江口北槽，我们利用“声学悬沙观测系统”，对河口细颗粒悬沙浓度垂向分布进行高时空分辨率观测。悬沙浓度的垂向分布有明显的变化，高含沙层存在于近底或中部。高含沙层是起因于与波、流相关的淤泥底床泥沙的再悬浮和水体中泥沙的沉降。悬沙浓度的垂向分布反映了向下沉降和引起再悬浮的向上紊流扩散