摘要:本文对长江口细颗粒泥沙过程研究进行了总结, 并提出了今后重点研究方向:1)长江口细颗粒泥沙是非均匀沙,其运动机理究竟如何?2)大规模水利工程究竟如何影响长江口最大浑浊带? 3)在长江口细颗粒泥沙过程的数学模拟中, 如何考虑河口波、流相互作用(耦合)及其对近底细颗粒泥沙输移的影响? 4)整个长江口水域瞬时、连续的水深、含沙量、地形变化资料的获取技术和方法的改进,以提高长江口细颗粒泥沙过程数学模拟的精度。 5)长江口悬沙以拉格朗日模式输运, 而过去大多悬沙观测调查是在欧拉模式中进行,如何进行欧拉和拉格朗日模式对比研究? 6)风暴潮、台风等对长江口细颗粒泥沙运动的影响。
关键词:长江口 最大浑浊带 细颗粒泥沙
1 引言 长江口是长江注入东海的入海口,自徐六泾以下经过三次分汊,共形成四个入海通道。崇明岛将长江口分为南支和北支;长兴岛和横沙岛又将南支分为南港和北港;南港又进一步被九段沙分为南槽和北槽(图1)。根据大通水文站多年统计资料,长江多年平均流量29500m3s-1。长江口外潮汐为正规半日潮,口内为非正规半日浅海潮。长江口为中潮河口,根据中浚站多年统计资料,多年平均潮差为2.66m。长江口波浪主要为风浪以及风浪和涌浪的混合浪。长江口水域的沿岸流主要为苏北沿岸流。长江口实测最大年输沙量为6.78亿t,最小年输沙量3.41亿t,年平均输沙量大约有4.66亿t[1]。每年由上游携带来的泥沙中有50%左右在长江口水下三角洲地区沉积下来,成为形成长江口拦门沙的主要成份。 |
|
长江口细颗粒泥沙过程的研究主要起因于: 1) 海洋科学陆海相互作用中长江河口动力学的基础理论研究;2) 长江口深水航道的维持(整治、疏浚); 3) 长江河口水环境、污染物的处理的日益恶化。长江口细颗粒泥沙过程主要研究径流、潮流、波浪、科氏力及沿岸流与细颗粒泥沙、淤泥质底床相互作用。长江口细颗粒泥沙过程研究可以揭示长江口淤泥底床侵蚀、堆积过程的机制[2]和长江河口演变规律[3],深入理解长江口细颗粒泥沙过程, 是海岸带陆海相互作用(LOICZ)的重要基础研究内容, 具有重要的科学理论意义。而且对港口航道回淤有重要影响的浮泥形成运动机制、时空分布规律得到进一步认识,并从而得到工程措施所需的科学依据,特别对上海市的长江口深水航道建设有重要现实意义。 长江口细颗粒泥沙过程研究成果大量来自海洋学家(河口海岸学、物理海洋学)、土木工程师(海岸工程师)、流体力学家(环境流体力学)的文献、著作。本文的目的: 1) 是力图把这些文献(以正式发表的文献为准,不包括研究报告)汇集起来,对长江口细颗粒泥沙过程的研究进行总结; 2) 究竟我们对长江口细颗粒泥沙过程了解多少?3) 究竟长江口细颗粒泥沙过程还有哪些问题值得研究?
2 长江口细颗粒泥沙特征
河口水体中含有大量细颗粒泥沙,即粒径小于0.004mm的粘土和一定比例的粉沙(粒径在0.0044~0.062mm之间)。细颗粒悬沙一般并不是以单颗粒的形式存在, 往往同附近其它大量的悬沙颗粒结合在一起,并形成一定的,这种相邻颗粒在一定条件下结合成集合体的作用被称为“絮凝作用”, 通过絮凝作用结合而成的这种体称为“絮凝体”。长江口悬沙或底沙,粒径均以&<0.032mm的颗粒为主,在悬沙中占90%以上,底沙中占75%[4]。长江口河床质主要由粉沙和粘土所组成,是一种粘性细颗粒非均匀沙,&<0.050mm的颗粒占88%, 其中&<0.005mm的粘土约占13%[2]。长江口细颗粒悬沙粒径0.032~0.016mm以下各粒级占90.5% [5]。
3 长江口细颗粒泥沙过程
不少学者对长江口细颗粒泥沙输移做了定性研究[6~9]。Milliman et al.[10~12]就长江口悬沙通量进行了半定量研究, 北槽似乎是悬沙进入海洋的主要通道[10]。长江口门附近的悬沙浓度和净输沙量在时间上有明显的涨落潮、大小潮和洪枯季变化。北港是长江输水输沙的主要通道,122°30'E-123°E是悬沙向东扩散的一条重要界线,入海泥沙主要向东偏南方向扩散[13]。借助于环行水槽的实验装置,沈承烈、阮文杰[2]对长江口天然河床质的冲淤特性和机理进行了初步研究。在动水淤积过程中,长江口悬移质与床沙有着粗细颗粒泥沙交换, 而在冲刷过程中无[2]。长江口细颗粒泥沙过程与环流有关[6,14~16]。基于现场水文泥沙调查资料,潘定安[17]就长江口南港的水文泥沙环境及污染物输移进行了探讨。
4 长江口悬沙平面扩散、悬沙锋
遥感卫星资料被应用于长江口入海悬沙扩散问题研究, 进而定量分析了长江口不同季节的悬沙浓度、扩散方向、扩散范围[8,18~20]。依据悬沙浓度分布,分四类:①高浑浊舌状构成长江向海冲泻的主体;②浑水带(122°25′~122°30′E);③羽状体(122°30′~122°35′E);④冲淡水带(122°40′~122°45′E)[8,18,19]。这四个带的变化受径流、风向、风级影响,主要由长江径流、台湾暖流、苏北沿岸流控制, 它们的输移方向是冬季朝南、夏季朝东北东。由于卫星照片能够反映同一时刻广阔区域内的悬沙分布情况,这一研究工作取得了很大进步。新的卫星图像资料(1988~1993)进一步提供了长江口表层悬浮泥沙浓度场空间分布和动态变化规律[21,22],并给出了表层悬浮泥沙浓度场的定量分布及其与垂线平均含沙量的关系。此外, 利用现场的调查资料,谷国传、曹沛奎[23]、邵秘华等[24]分析了长江口及其毗邻区域悬沙成份和浓度平面分布。利用长江口水文泥沙现场观测资料, 曹沛奎、严肃庄[25]对长江口含沙量梯度较大的悬沙锋进行了研究, 发现了内、外两个锋面。悬沙锋对长江口的物质输移和冲淤变化起积极作用。它们的变化规律仍有待深入研究。
转贴于5 长江口最大浑浊带
自80年代,河口学家对长江口最大浑浊带细颗粒泥沙输移过程进行了深入的实验研究[26]。沈焕庭等[27]系统分析了长江口最大浑浊带形成的环境背景、时空变化规律、泥沙来源和絮凝作用对悬沙落淤的影响、浮泥的特性与分布以及悬沙的富集机制。时伟荣[28]、Shi and Li[29]利用长江口南槽底层含沙量资料,分析计算了不同潮时的泥沙垂向紊动扩散系数,研究了底床泥沙再悬浮与最大浑浊带形成的关系。时伟荣、李九发[30]运用机制分解的方法计算了长江口南北槽浑浊带区域内的余流和悬沙单宽输移量,并据此探讨了浑浊带的发育过程。通过对长江口不同河段的悬沙特性(大小、沉速、含沙量)和输移规律的分析对比, 李九发等[31]、Li and Zhang[32]探讨了这些因素对长江口最大浑浊带形成过程的重要性,认识到最大浑浊带潮流强劲,引起床沙再悬浮。沈健等[33]对长江口最大浑浊带的水沙输运机制进行了定量分析。贺松林、孙介民[34]研究了长江口最大浑浊带悬沙输移过程中的“潮泵效应”(在涨落潮中,悬沙颗粒在水体与底床之间周期性地做上下悬扬、沉降的现象), 并指出南、北槽之间的大尺度平面环向悬沙输移和南、北槽内的次级尺度的平面环向悬沙输移的重要性。Zhou and Wu[35]、周华君[36]利用通量分析方法研究了长江口最大浑浊带的形成机制。结果表明支持最大浑浊带的长周期泥沙源是陆向雷诺(Reynolds)输运和较低层的陆向余流(平流)输运。潘定安等[37]研究了长江口浑浊带的形成机理与特点,指出径流潮流相互作用形成潮汐浑浊带,而盐淡水交汇混合形成盐水浑浊带。是否这两种机制是可以统一起来?在这些研究中, 所采用的方法主要是现场实验观测, 悬沙浓度资料大多依赖于现场水样的采集、室内分析。利用现场水文、悬沙资料,并结合数学模型,时钟、陈伟民(2000)[38]探讨了长江口北槽最大浑浊带的泥沙过程。研究结果表明:在北槽口内,最大浑浊带形成的主要动力过程是潮汐的不对称性和河口重力环流。在北槽口外,最大浑浊带形成的主要动力过程则是河口底部泥沙的周期性再悬浮。在长江口北槽口内、口外最大浑浊带中,细颗粒泥沙的再悬浮过程也存在着一定的周期性。此外,由盐度、悬沙浓度层化引起的“层化抑制紊流”也是长江口北槽口内、口外最大浑浊带的成因机制之一。
从以上研究可见,长江口最大浑浊带细颗粒泥沙运动主要表现为:1)垂向交换、2)南北槽平面输移3)水平富集沉降。这几种运动方式分别发生于何种条件下呢?作者认为可能同时发生, 从已有的文献看, 这一问题有待进一步研究。
6 长江口浮泥特性
河口学家对长江口浮泥也进行了现场实验观测研究[39,40]。通过对提纯粘粒做X射线、电子显微镜、差热和脱水分析,张志忠等[39]定性地鉴定出长江口浮泥的粘土矿物成份,并探讨了浮泥的形成条件和分层沉淀现象。结果发现:① 长江口南槽底质中粘土的矿物成分主要是伊利石;② 长江口的浮泥可能主要是汛季上游带来大量粘土的细颗粒物质在长江口受咸水环境作用絮凝沉降所致;③ 长江口浮泥层是由浓度高、密度大、颗粒细的物质组成的分层的半流体的透镜体。长江口浮泥的特性:① 长江口浮泥面容重1.04g/cm3,底板容重1.25g/cm3;② 组成浮泥的细颗粒泥沙中值粒径6.44~6.97;③ 组成浮泥泥沙的粘土矿物为伊利石-绿泥石 高岭石 蒙脱石组合;④ 自然状态下一般浮泥厚度为0.1~0.4m[40]。
7 长江口细颗粒泥沙絮凝机理和絮凝体沉降速率
细颗粒泥沙絮凝是潮汐河口的普遍现象, 其机理具有多样性, 如:盐度、高含沙量条件下颗粒碰撞、有机物或生物作用引起的絮凝,湍流强度变化产生的变化。长江口细颗粒泥沙絮凝主要侧重于盐度絮凝, 絮凝现象的产生必须同时具备两个条件:①絮凝体的电化学反应和颗粒之间的碰撞;②絮凝体所受剪切力小于其抗剪强度[41,42]。实验表明,产生絮凝的最小悬沙粒径是0.030mm,当粒径处于0.01至0.03mm时, 絮凝作用是很微弱的。由于絮凝作用,颗粒大小可以在0.001mm到0.1mm范围内迅速变化。
长江口细颗粒泥沙絮凝机理的研究也得到了发展[5,41~54]。通过现场新鲜水样进行的显微摄影,关许为等[46]测得长江口不同区域絮凝体分布、形状、和尺寸,并分析了盐度分布和潮汐现象对泥沙絮凝现象的影响。絮凝体直径分布在0.01~0.5mm之间,大部分出现在0.01~0.2mm范围内。絮凝团主要有条状、团状和蜂窝状三种形式。盐淡水交汇锋面附近、拦门沙河段临底区域为出现大量絮凝体的区域。低流速情况下,特别是憩流期,絮凝较为显著;而高流速情况下,例如涨急与落急期,絮散比较明显。根据长江口细颗粒泥沙絮凝沉降试验和长江口一年的实测阳离子浓度资料,蒋国俊、张志忠[48]分析了阳离子浓度对细颗粒泥沙动力絮凝沉降的影响。结果发现:① 长江口细颗粒泥沙动力絮凝存在最佳离子浓度;② 在阳离子浓度相同的水体中,高价离子比低价离子更能促进细颗粒泥沙的动力絮凝沉降;③ 长江口阳离子浓度的时空变化,影响细颗粒泥沙絮凝沉降强度和沉积部位,对拦门沙的发育及冲淤变化有重要影响。
Han and Lu[43]得出如下长江口细颗粒泥沙絮凝沉速的综合关系式
ω/ω0=26.5(1/K·Cmf/C)-0.8(絮凝沉降加速阶段) |
(1) |
ω/ω0=26.5(1/K·Cmf/C)0.892(絮凝沉降减速阶段) |
(2) |
式中ω0为中值粒径泥沙颗粒沉降速度;K为经验系数;C为悬沙浓度;Cmf为最适絮凝悬沙浓度。
通过对长江口细颗粒泥沙静水絮凝沉速试验研究, 彭润泽等[44]得出如下长江口泥沙絮凝沉速的综合关系式
F=ωf50/ω50=(df50/d50)2=1+b0cb1sb2Hb3d50b4 |
(3) |
式中 F静水絮凝系数,等于ωf50/ω50, 沉距H, d50中值粒径,ω50中值沉速,ωf50絮凝级配的中值粒径沉速,df50表示wf50对应的当量粒径,bi,i=0,1,2,3,为回归系数。
依据悬沙扩散理论,李炎等[55]、Li等[56]提出运用分粒级的泥沙垂线分布推算长江口各粒级相对沉速、絮凝因素比和絮凝指标的方法。研究结果表明:① 长江近口段悬浮泥沙处于絮凝程度低,群体沉速较低的亚稳态;② 近口段至河口段间盐度急剧变化的区间,存在着一个絮凝程度和沉降特征的迅速变化带;③ 河口段最大浑浊带的絮凝程度高,群体沉速大,各粒级沉速接近,絮凝团形成已近极限状态;④ 口外海滨段盐度大于20‰时絮凝程度又有所降低。利用流速仪和声学悬浮泥沙观测系统,时钟等[57]获得了长江口北槽口外大潮水流、悬沙浓度垂线分布资料, 通过Rouse公式拟合法计算得到细颗粒悬沙沉降速度为3.0~4.0mms-1, 此值可能偏大, 更多组的计算仍需加强。
8 长江口近底边界层细颗粒泥沙运动
在淤泥质河口近底边界层,强劲的潮流、波流能引起淤泥底床沉积物的侵蚀或再悬浮,产生悬沙浓度的垂直分层现象,反过来又影响潮流、波流。占水体10~20%的近底高含沙层(浓度5~10kgm-3), 这作为河口细颗粒泥沙输移的主要形式, 普遍存在于许多高浑浊的河口内。深入理解河口底部边界层高含沙层的动力学机制具有重要的科学意义,这将加深我们对河口近底水-泥相互作用的认识,而这种作用控制着整个河口细颗粒泥沙的输移过程。
近年来,海洋水声学被成功应用于长江口近底边界层细颗粒泥沙输移的实验研究[58-64]。这应为深入理解河口近底边界层高含沙层的侵蚀、挟运、输移和堆积过程提供新的实验研究途径[65],对河口最大浑浊带成因的进一步认识应有重要的科学意义。利用声学悬浮泥沙观测系统, 在长江口深水航道北槽口内(横沙东)、口外(牛皮礁)获取到近底高浓度悬浮泥沙高时、空分辨率连续分布图象及运动变化过程;并得到了大潮汛时段高浓度含沙水体(2~5 kgm-3)、流动浮泥层(浓度5~10kgm-3)、固定浮泥层(&>10kgm-3)及粘性淤泥底床三个不连续面的判别指标, 这些成果对泥沙再悬浮、近底泥沙输移, 浮泥的形成与消散、底沙与悬沙交换等问题的深入研究开创了新的实验途径。具体地讲,有以下几个方面的进展:
8.1 连续的、高时空分辨率的河口细颗粒悬沙浓度垂线分布
在长江口深水航道北槽,利用“声学悬沙观测系统”观测到大潮典型的高时空分辨率细颗粒悬沙浓度声学垂向分布图和垂线分布曲线。研究结果表明: ①在涨潮时呈L形,悬沙浓度的垂向变化梯度小(浓度小于1.0kgm-3)。② 在涨急时呈射流形,射流顶的悬沙浓度达10 kgm-3,悬沙浓度的垂向变化梯度大。③ 在接近涨憩时,悬沙浓度的垂向梯度小。④ 在落潮时从水面到水底悬沙浓度按指数增加,代表恒定均匀流中悬沙处于平衡条件的分布,泥沙垂向扩散系数εs在数值上大约是泥沙颗粒沉降速度ωs的两倍。⑤ 在落急时不连续的,水体中部呈射流形,悬沙浓度的垂向变化梯度大。⑥ 在接近落憩时,悬沙浓度的垂向梯度适中。
8.2 河口近底细颗粒泥沙再悬浮、再挟运过程
利用声学悬浮泥沙观测仪获得的长江口粘性悬沙浓度垂向分布图表明,近底高浓度悬沙层,以一个浓度约为2kgm-3的天然分界点被分为上、下两层。声学图象还揭示了周期为数秒的近底高频率再悬浮现象,这一过程主要促进了近底高浓度悬浮液的形成。上、下泥跃层也分别显示了上层和下层高浓度悬浮泥沙的再挟运。近底高频率再悬浮可能是由湍流造成的,而再挟运可能归因于内波作用。
8.3 河口近底高含沙层(浮泥)和高、低频率界面波
利用一个悬浮泥沙声学观测仪(ASSM),在大潮时,对长江口外泥沙浓度进行了实验观测。结果发现底部泥沙存在分层:从上往下分别是低浓度悬沙层;高浓度悬沙层和浮泥层。声学图像显示有四种不同频率的界面波: ① 介于低浓度悬沙层和高浓度悬沙层之间的低频界面波。② 介于高浓度悬沙层和流动浮泥层之间的低频内波。③ 介于流动浮泥层和固定浮泥层之间的波,被谨慎地解释为低频界面波。④ 在固定浮泥层上产生的是高频界面波。这些界面波起到次级的动力作用,增强了局部粘性泥沙的再悬浮和再挟运。
在长江口细颗粒泥沙输移过程现场实验研究中, 由于絮凝作用,细颗粒泥沙絮凝体大小发生变化,这将会如何影响悬沙浓度的声学估算精度呢? 在长江口细颗粒泥沙输移过程现场实验中, 与泥沙浓度相匹配的连续的、高时空分辨率的水流流速、密度垂线分布资料仍然缺乏。国内外河口水沙数学模型大多利用六点法测得的流速、泥沙浓度资料验证, 精度低、误差大。如何利用声学悬沙浓度现场实验资料并结合水动力资料验证长江口水沙数学模型以及现有的河口细颗粒悬沙紊流扩散理论? 如何定量研究长江口近底高含沙层和上层水体之间的泥沙交换通量以及近底高含沙层和淤泥底床之间泥沙交换通量?
尽管在过去20多年中长江口细颗粒泥沙运动研究取得了很大的进展,但仍有深入研究的必要。这是因为长江口细颗粒泥沙运动本身是非常复杂的。河口悬沙以拉格朗日(Lagrangian)模式输运, 但是, 大多悬沙调查是在欧拉(Eulerian)模式中进行。这表明,悬沙浓度的变化可归因于泥沙平流以及底床泥沙侵蚀。由于细颗粒泥沙的絮凝是依赖于悬沙浓度和湍流。所以,将悬沙浓度变化的平流成分和湍流成分分开是十分重要的。这样, 能定量地表明这两种参数对絮凝过程的意义。从而,可以经验性地进行欧拉和拉格朗日模式对比研究。此外,还有利于对絮凝过程的了解,对影响泥沙分布、输移的水流的认识有利,从而加强对细颗粒泥沙运动过程的了解。
9 长江口细颗粒泥沙运动的数学模拟
从90年代开始, 河口学家、海岸工程师和环境流体力学家对长江口细颗粒泥沙运动及冲淤变化进行了数学模拟研究,这些包括垂向一维悬沙数学模型、垂向二维悬沙数学模型、平面二维悬沙数学模型等等。
9.1 垂向一维悬沙数学模型
在实测资料分析的基础上,认为长江口地区瞬时含沙量的变化主要取决于泥沙沿垂线方向的扩散和沉降,徐建益等[66]对非恒定二维泥沙扩散模型进行了简化,建立了长江口南支非均匀沙垂向分层的数学模型,通过求解垂向扩散方程给出垂线上分层的含沙量和悬沙粒径级配。利用垂向一维悬沙数学模型,时钟、陈伟民[38]对长江口北槽口内、口外的细颗粒泥沙过程进行了模拟。
9.2 垂向二维悬沙数学模型
利用一维水流、垂向二维泥沙数学模型,黄永健[67]进行了水流越槽的不平衡输沙计算,求出沿水流各断面垂向含沙量分布,从而推算出挖槽的回淤量。类似地,利用一个一维水流、垂向二维泥沙数学模型,李家春等[68]、王涛、李家春[69]、周济福等[70]研究了长江口北槽拦门沙地区不同径流、潮流组合作用下泥沙输运规律。李家春等[68]分析了河口泥沙输移受波流边界层、湍流泥沙相互作用、河口混合过程的影响,在此基础上建立了河口泥沙输运的模型。
9.3 平面二维悬沙数学模型
利用平面二维悬沙数学模型,赵棣华等[71]对长江口南支悬沙运动进行了模拟。在大范围平面二维水流和含沙量计算的基础上,将所关心的局部区域视为一个封闭的区域,通过分别计算其各条边界上的泥沙输移量,来预测该区域内泥沙冲淤分布及河床变形的趋势[72]。
9.4 数理统计模型
选用长江口南槽逐月水下地形图、相应的大通站的流量和含沙量,中浚、横沙和高桥等站的潮位资料,应用数理统计理论和方法,朱惠芳、周纪芗[73]对长江口南槽航道拦门沙的冲淤变化进行了数学模拟。
长江口受径流、潮流、沿岸流影响,是非定常的往复流动。同时,湍流对泥沙输运有重要影响。湍流边界层中的猝发现象与泥沙起动有关,层结现象又可抑制剪切产生的湍流脉动。在那里,由于淡水和盐水混合,又会导致絮凝、加速沉降等。因此,河口边界层中流动,泥沙与底床的相互作用是关键问题。因长江泥沙浓度较低,泥沙运动可视为被动标量,而且冲淤计算的时间尺度较长,所以流场、浓度、冲淤计算可依次进行。结果表明,流场分布与观测资料定量一致,浓度分布虽然定性一致,但在定量上尚有一些误差。河口泥沙研究必须建立在波流场、分层流、两相流、湍流研究成果的基础上[74]。今后,仍应加强长江口最大浑浊带的数学模拟的研究,并在模拟中考虑水流和悬沙的耦合。此外,最大浑浊带中悬沙会影响阻力系数和底床切应力的分布。近年来,人工神经网络也不断应用到水科学中,能否将其应用到长江口水动力、泥沙过程研究中?蔡树棠等[75]从两相流体运动的观点出发来讨论潮汐河口泥沙的运动和淤积问题。今后,能否考虑用两相流运动理论研究长江口泥沙运动呢?
10 长江口入海航道疏浚、回淤和治理
沈焕庭等[76]探讨了长江河口入海航道治理研究。顾伟浩等[77]对长江口铜沙老航槽自然回淤进行了初步分析。陈吉余[78]深入讨论了长江口治理。1984年2~3月以前,长江口入海航道原为南槽,之后,开挖北槽航道,当年4月25日正式通航。根据现场观测的水文、泥沙资料和数值模拟计算结果,盛升国等[79]对长江口北槽航道抛泥区疏浚泥沙的运动进行了研究。利用平面二维数学模型,盛升国等[79]对抛入不同抛泥区的疏浚土,在潮流作用下的扩散方向和范围做了研究。顾伟浩[80]分析了台风对长江口北槽航槽挖槽段回淤影响。钟修成、任苹[81]对长江口拦门沙航道(北槽)回淤进行了分析。根据1990年2~4月的现场实测水文泥沙资料及有关航槽检测图,关许为、顾伟浩[82]就寒潮对长江口北槽回淤影响进行了分析。径流、潮流、风浪和风吹流是长江口航槽冲淤四大主要动力因素,泥沙运动是这些动力因素综合作用的结果。他们认为只有7级以上风力的寒潮天,才会对长江口航槽淤浅产生显著影响。通过大量现场实测资料,张栋梁、孙建国[83]长江口北槽挖槽段泥沙淤积机理探讨。张栋梁、姚金元[84]就长江口北槽挖槽段泥沙淤积特性进行了分析研究。研究表明,北槽挖槽段泥沙淤积呈三峰型,南槽、北港水沙与北槽交换的影响使得北槽拦门沙位置上移。三峰型淤积机理各不相同,如洪汛时泥沙在北槽全槽淤积,汛后期泥沙在挖槽段内淤积,实质是洪汛北槽下段的淤积泥沙在时间和空间上的调整。
11 结论
本文对长江口细颗粒泥沙过程研究进行了回顾总结。在过去20年中, 长江口细颗粒泥沙特征、细颗粒泥沙过程、悬沙平面扩散、悬沙锋、最大浑浊带、浮泥特性、细颗粒泥沙絮凝机理和絮凝体沉降速率、近底边界层细颗粒泥沙运动和细颗粒泥沙运动的数学模拟等方面取得了一定成就, 为长江口入海航道疏浚、治理等提供了科学依据。 海洋水声学被成功应用于长江口近底边界层细颗粒泥沙输移的实验研究这应为深入理解河口近底边界层高含沙层的侵蚀、挟运、输移和堆积过程提供新的实验研究途径,对河口最大浑浊带成因的进一步认识应有重要的科学意义。数学模拟已成为长江口细颗粒泥沙过程研究及应用的重要手段。
就作者的认识和兴趣, 今后应在以下几个方面加强进一步研究:大规模工程势必对长江口悬沙场等产生影响,究竟对长江口最大浑浊带的影响程度如何?在长江口最大浑浊带细颗粒泥沙过程的数学模拟中, 如何考虑河口波、流相互作用(耦合)及其对近底细颗粒泥沙输移的影响?长江口地形变化大,整个水域瞬时的水深、地形变化资料很难获取;此外,整个水域瞬时、连续的含沙量资料缺少,这些都给长江泥沙过程数学模拟精度带来不少问题。整个长江口水域瞬时、连续的水深、含沙量、地形变化资料的获取技术和方法的改进,可以提高长江泥沙过程数学模拟精度。长江口悬沙以拉格朗日(Lagrangian)模式输运, 而过去大多悬沙观测调查是在欧拉(Eulerian)模式中进行,如何进行欧拉和拉格朗日模式对比或转换研究?极端气候条件下,如风暴潮、台风等对长江口细颗粒泥沙运动的影响。
参 考 文 献
[1] 林承坤。长江口泥沙的数量与输移。中国科学, 1988, 1: 104-112.
[2] 沈承烈,阮文杰。长江口河床质冲淤特性的试验研究。泥沙研究, 1986,(3) 62-72.
[3] 黄胜。长江河口演变特征。泥沙研究, 1986,(4),1-12.
[4] 吕全荣,王效京。长江口细颗粒沉积物的粘土矿物及地球化学特征。沉积学报, 1985, 3 (4): 141-151.
[5] 张志忠。长江口细颗粒泥沙基本特性。泥沙研究, 1996,(1),67-73.
[6] Shen, H.T., Zhu, H.F. and Mao, Z.C. Circulation of the Changjiang and its effect on the transport of suspended sediment. In: V.S. Kennedy (Ed.), Estuarine Comparisons, Academic Press, New York, 1982, 677- 691.
[7] Shen, H.T., Li, J.F., Zhu, H.F. and Zhou, F.G. The suspended sediment transport in the Changjiang Estuary. Proceedings of International Symposium on Sedimentation on the Continental Shelf (Hangzhou, China). China Ocean Press, Beijing, 1983, 359-369.
[8] Yun, C.X. and Wan, J.R. A study of diffusion of upper layer suspended sediments in discharges from the Changjiang estuary into the sea, based on satellite imagery. In: V.S. Kennedy, (Ed.), Estuarine Comparisons, Academic Press, New York, 1982, 693-704.
[9] 张志忠,徐志刚。长江口悬沙及其运移。海洋科学, 1983, 5: 6-11.
[10] Milliman, J.D., Yang, Z.S. and Meade, R.H. Flux of suspended sediment in the Changjiang Estuary. Proceedings of the International Symposium on Sedimentation on the Continental Shelf with Special Reference to the East China Sea. China Ocean Press, Beijing, 1983: 413-430.
[11] Milliman, J.D., Hsueh, Y., Hu, D.X., Pashinski, D.J., Shen, H.T., Yang, Z.S。and Hacker, P. Tidal phase control of sediment discharge from the Yangtze River. Estuarine, Coastal and Shelf Science , 1984, 19: 119-128.
[12] Milliman, J.D., Shen, H.T., Yang, Z.S. and Meade, R.H. Transport and deposition of river sediment in the Changjiang estuary and adjacent continental shelf. Continental Shelf Research , 1985, 4: 37-45.
[13] 沈焕庭,李九发,朱慧芳,周福根。长江河口悬沙输移特性。泥沙研究, 1986, (1), 1-13.
[14] 沈焕庭,朱慧芳,茅志昌。长江河口环流及其对悬沙输移的影响。海洋与湖沼, 1986, 17 (1): 26-35.
[15] Su, J.L. and Wang, K.S. The suspended sediment balance in Changjiang Estuary. Estuarine, Coastal && Shelf Science,1986, 23: 81-98.
[16] 王康NFDA2,苏纪兰。长江口南港环流及悬移质输运的计算分析。海洋学报, 1987, 9: 627-637.
[17] 潘定安。长江口南港的水文泥沙环境及污染物输移。地理研究, 1996, 15 (1): 39-45.
[18] 恽才兴,万嘉若。应用卫星象片分析长江入海水体表层悬浮泥沙的扩散。遥感文选, 北京,科学出版社, 1981: 175-179.
[19] 恽才兴,蔡孟裔,王宝全。利用卫星象片分析长江入海悬浮泥沙扩散问题。海洋与湖沼, 1981, 12 (5): 391-478.
[20] 陈鸣,李士鸿,刘小靖。长江口悬浮泥沙遥感信息处理和分析。学报, 1991,(5),47-51.
[21] He, Q. and Yun, C.X. Preliminary study on transport and distribution of sediment load in the Yangtze Estuary. International Journal of Sediment Research , 1998, 13 (4): 59-69.
[22] 何青,恽才兴,时伟荣。长江口表层水体悬沙浓度场遥感分析。自然科学进展, 1999, 9 (2): 160-165.
[23] 谷国传,曹沛奎。长江口外水域悬沙分布特征。东海海洋, 1986, 4 (1): 10-15.
[24] 邵秘华,李炎,王正方,张素香。长江口海域悬浮物的分布时空变化特征。海洋环境科学, 1996, 15 (3): 36-39.
[25] 曹沛奎,严素庄。长江口悬沙锋及其对物质输移的影响。华东师范大学学报, 1996, 1: 85-94.
[26] 沈焕庭,郭成涛,朱慧芳,徐海根,恽才兴,陈邦林。长江河口最大浑浊带的变化规律及其成因探讨。中国海洋湖沼学会河口海岸分会编, 海岸河口区动力、地貌、沉积过程论文集, 科学出版社, 1984: 76-89.
[27] 沈焕庭,贺松林,潘定安,李九发。长江河口最大浑浊带研究。地理学报, 1992, 47 (5): 472-479.
[28] 时伟荣。长江口浑浊带含沙量的潮流变化及其成因分析。地理学报, 1993, 48 (5): 412-419.
[29] Shi, W.R. and Li, J.F. Sediment resuspension in the Changjiang Estuary. International Journal of Sediment Research, 1995, 10 (1): 32-42.
[30] 时伟荣,李九发。长江河口南北槽输沙机制及浑浊带发育分析。海洋通报, 1993, 12 (4): 69-76.
[31] 李九发,时伟荣,沈焕庭。长江河口最大浑浊带的泥沙特性和输移规律。地理研究, 1994, 13 (1): 51-59.
[32] Li, J.F. and Zhang, C. Sediment resuspension and implications for turbidity maximum in the Changjiang Estuary. Marine Geology, 1998, 148: 117-124.[33] 沈健,沈焕庭,潘定安,肖成猷。长江河口最大混浊带水沙输运机制分析。地理学报, 1995, 50: 411-419.
[34] 贺松林,孙介民。长江口最大浑浊带的悬沙输移特征。海洋与湖沼, 1996, 27 (1): 60-65.
[35] Zhou, H.J. and Wu, S.R. Mechnisms of formation and maintenance of the turbidity maximum in the Changjiang Estuary. Proceedings of the Fourth Pacific/Asia Offshore Mechanics Symposium, ISOPE, Pusan, Korea, 1996: 99-104.
[36] 周华君。潮汐河口悬浮物通量分解方法探讨。周连第,邵维文、戴世强(主编), 第十二届全国水动力学研讨会文集,北京,海洋出版社, 1998:371-376.
[37] 潘定安,沈焕庭,茅志昌。长江口浑浊带的形成机理与特点。海洋学报, 1999, 21 (4): 62-69.
[38] 时钟,陈伟民。长江口北槽最大浑浊带泥沙过程。泥沙研究, 2000, (1), 28-39.
[39] 张志忠,黄文盛,杨晓顺,沈金山,龚尧圭。长江口浮泥若干特性的初步研究。杭州大学学报, 1977, 1: 83-96.
[40] 徐海根,徐海涛,李九发。长江口浮泥层“适航水深”初步研究。华东师范大学学报, 1994, 2: 91-97.
[41] 阮文杰。细颗粒泥沙动水絮凝的机理分析。海洋科学, 1991, 5: 46-49.
[42] 阮文杰。长江口天然水流中细颗粒泥沙的絮凝作用。海洋科学, 1991, 6: 39-43.[43] Han, N.B. and Lu, Z.Y. Settling properties of the sediments of the Changjiang Estuary in salt water. Proceedings of International Symposium on Sedimentation on the Continental Shelf (Hangzhou, China). China Ocean Press, Beijing, 1983: 483-493.
[44] 彭润泽,黄永健,蒋如琴,张振秋,颜燕。长江口泥沙静水絮凝沉速试验研究。 水电科学研究院科学研究论文集,北京,电力出版社,1990, 33: 67-73.
[45] 关许为,陈英祖。长江口泥沙絮凝静水沉降动力学模式的试验研究。海洋工程, 1995, 13 (1): 46-50.
[46] 关许为,陈英祖,林以安,李炎。长江口泥沙絮凝体的现场显微观测。泥沙研究, 1992, 3: 54-59.
[47] 关许为,陈英祖,杜心慧。长江口絮凝机理的试验研究。学报, 1996, 6: 70-74.
[48] 蒋国俊,张志忠。长江口阳离子浓度与细颗粒泥沙絮凝沉积。海洋学报, 1995, 17: 76-82.
[49] 张志忠。长江口细颗粒泥沙基本特性研究。全国泥沙基本理论研究学术讨论会论文集,北京, 中国建材工业出版社, 1992: 369-375.
[50] 张志忠,王允菊,徐志刚。长江口细颗粒泥沙絮凝若干特性探讨。第二届河流泥沙国际学术讨论会论文集, 北京:电力出版社, 1983: 274-285.
[51] 张志忠、阮文杰、蒋国俊。长江口动水絮凝沉降与拦门沙淤积的关系。第七届全国海岸工程学术讨论会论文集, 北京,海洋出版社, 1993: 378-383.
[52] 张志忠,阮文杰,蒋国俊。长江口动力絮凝沉降与拦门沙淤积的关系。海洋与湖沼, 1995, 26 (6): 632-642.
[53] 林以安,李炎,唐仁友。长江口絮凝聚沉特征与颗粒表面物理化学因素作用-I悬浮颗粒絮凝沉降特征。泥沙研究, 1997,(1),42-47.
[54] 林以安,唐仁友,李炎。长江口絮凝聚沉特征与颗粒表面物理化学因素作用-II 颗粒表面性质对聚沉的作用。泥沙研究, 1997,(4),76-84.
[55] 李炎,王亚琴,林以安。运用分粒级的泥沙垂线分布分析长江口泥沙絮凝沉降特征。海洋工程, 1994, 12 (1): 71-78.
[56] Li, Y., Wang, Y.Q. and Lin, Y.A. Investigation on settling characteristics of flocs in the Changjiang Estuary, China. China Ocean Engineering , 1995, 9(3): 325-334.
[57] 时钟,朱文蔚,周洪强。长江口北槽口外细颗粒悬沙沉降速度。上海交通大学学报, 2000, 34 (1): 18-23.
[58] 时钟。河口粘性泥沙再悬浮和再挟运的声散射观测。第八届全国海岸工程学术讨论会暨1997海峡两岸港口及海岸开发研讨会论文集, 北京: 海洋出版社, 1997: 67-74.
[59] 时钟,周洪强。长江口深水航道北槽口外悬沙浓度垂向分布。上海交通大学学报, 1997, 32 (11): 132-138.
[60] 时钟,凌鸿烈。长江口细颗粒悬沙浓度垂向分布。泥沙研究, 1999,(2),59-64
[61] 时钟,张叔英,Hamilton, L.J. 河口近底细颗粒悬沙运动的声散射观测。声学学报, 1998, 23 (3): 221-227.
[62] Shi, Z., Ren, L.F. and Lin, H.L. Vertical suspension profile in the Changjiang Estuary. Marine Geology, 1996, 130: 72-81
[63] Shi, Z., Ren, L.F., Chen, J.Y. and Zhang, S.Y. Acoustic imaging of cohesive sediment resuspension and re-entrainment in the Changjiang Estuary, East China Sea. Geo-Marine Letters, 1997, 17: 162-168.
[64] Shi, Z., Ren, L.F. and Hamilton, L.J. Acoustic profiling of fine suspension concentration in the Changjiang Estuary. Estuaries, 1999, 23 (3A): 648-656.[65] Hamilton, L.J., Shi, Z. and Zhang, S.Y. Acoustic backscatter measurements of estuarine suspended cohesive sediment concentration profiles. Journal of Coastal Research, 1998, 14 (4): 1213~1224.
[66] 徐建益,陶学为,方良田,楼越平,韩曾萃。长江口南支非均匀沙垂向分层的数学模型。泥沙研究, 1995, 2: 74-79.
[67] 黄永健。长江口挖槽自然回淤的计算。泥沙研究, 1997, (2),69-73.
[68] 李家春,周济福,王涛。河口泥沙输运研究。‘97海岸海洋资源与环境学术研讨会论文集,香港, 1997: 15-24.
[69] 王涛,李家春。长江口拦门沙地区泥沙输运。力学与实践, 1997, 19 (6):33-36.
[70] 周济福,王涛,李家春,刘青全。径流与潮流对长江口泥沙输运的影响。水动力学研究与进展, 1999, 14 (1): 90-100.
[71] 赵棣华,谭仁忠,谭维炎。长江口南支河段悬移质含沙量计算模型。泥沙研究, 1990, (2), 54-62.
[72] 徐建益,袁建中。长江河口局部区域内泥沙冲淤计算方法及应用。学报, 1992, (12), 63-69.
[73] 朱慧芳,周纪芗。长江河口航道拦门沙冲淤变化的数学模拟和预测。海洋与湖沼, 1993, 24 (3): 294-303.
[74] 李家春,吴承康。环境力学与可持续发展。庄逢甘 (主编), 现代力学与科技进步,清华大学出版社,北京:1997: 35-42.
[75] 蔡树棠,黄菊卿。潮汐对河口泥沙运动和淤积的影响。水动力学研究与进展, 1986, 1 (2): 1-9.
[76] 沈焕庭,徐海根,马相奇。长江河口入海航道治理研究。海洋科学, 1983, (4),5-9.
[77] 顾伟浩,姚金元,袁金林,曾守源。长江口铜沙老航槽自然回淤的初步分析。海洋科学, 1984,(5), 10-13.
[78] 陈吉余。论长江口治理。科技导报, 1993, 3: 56-58.
[79] 盛升国,施厚庆,于福生,黄永健。长江口北槽航道抛泥区泥沙扩散的试验研究。泥沙研究, 1986, (3), 37-45.
[80] 顾伟浩。台风对长江口北槽航槽挖槽段回淤影响的初步分析。水运工程, 1986, (9),32-50.
[81] 钟修成,任苹。长江口拦门沙航道(北槽)回淤分析。河海大学学报, 1988, 16 (6): 50-57.
[82] 关许为,顾伟浩。一九九0年寒潮对长江口北槽回淤影响的分析。泥沙研究, 1992,(1), 55-60.
[83] 张栋梁,孙建国。长江口北槽挖槽段泥沙淤积机理探讨。水运工程, 1991, (11), 27-35.
[84] 张栋梁,姚金元。长江口北槽挖槽段泥沙淤积特性研究。泥沙研究, 1993, (2),66-78.
转贴于