波浪及波流边界层泥沙起动规律

基于波浪边界层理论及单向流泥沙起动Shields曲线,推证出波浪泥沙起动Shields曲线;基于波流边界层理论,提出表述波流边界层动力特征的波流比因子X及非线性作用因子Y,并建立了Y与X的相关关系;在此基础之上,结合单向流及波浪泥沙起动Shields曲线,推证出波流共同作用下泥沙起动Shields曲线。结果表明:波浪泥沙起动Shields曲线在层流区与单向流光滑紊流区曲线保持一致,粗糙紊流区与单向流粗糙紊流区曲线保持一致,过渡区线型为折线,由层流区及粗糙紊流区曲线延长交汇获得;X及Y能够合理地表征波流边界层动力对比特征及非线性作用特征;波流泥沙起动Shields曲线介于波浪及单向流泥沙起动Shields曲线之间,随着波流比因子X的不同,依据非线性作用因子Y,自动在波浪及单向流泥沙起动Shields曲线之间非线性过渡。建立的波流泥沙起动Shields曲线与试验结果吻合较好,且能够概括单向流、波浪及波流等不同动力及细沙、粗沙等不同粒径的泥沙起动条件。     

单向流边界层泥沙起动规律

Shields曲线常用于表示泥沙起动的临界条件,基于边界层理论,对Shields曲线各个流区的线型进行了推证;考虑粘结力的作用,对Shields参数及Shields曲线进行了修正,并给出修正Shields曲线表达式;在此基础之上,从边界层角度重新阐述了Shields曲线。结果表明:Shields曲线在光滑紊流及层流区呈直线分布,在过渡区与阻力系数线型保持一致,在粗糙紊流区呈水平直线分布;修正后Shields曲线与原始Shields曲线在形式上保持一致,修正Shields曲线表达式与实测数据吻合较好,适用于粗、细颗粒泥沙起动条件的计算;Shields曲线事实上代表了Shields参数与沙粒周围绕流流态的关系,同一颗粒处于不同流区起动时,其起动切应力不同。     

波流边界层水沙运动数值模拟——II.泥沙运动模拟

为解析波流边界层内泥沙运动,建立了基于水动力-泥沙-床面互馈过程的波流边界层1DV泥沙数学模型,可用于模拟不同床面形态下粉沙-沙的含沙量过程。床面形态模块提供床面形态类型和相应参数;给出了平底和沙波床面粗糙高度和泥沙扩散系数的确定方法;采用了适宜粉沙及沙的制约沉速、底部参考浓度和起动剪切应力等公式;引入含沙量层化效应和制约沉降反映水动力与泥沙之间的相互影响。水槽试验资料验证表明,建立的模型较好地模拟了不同床面不同波流组合条件下的含沙量剖面。在此基础上,讨论了不同床面含沙量剖面模拟方法的差异,指出床面形态是决定含沙量变化的重要因素之一,仅通过改变床面粗糙高度不足以反映漩涡沙波床面的含沙量剖面特征。该模型可为研究波流边界层内泥沙运动和物质输运提供工具。     

层流区Shields曲线分布规律试验

为探明层流条件下Shields曲线的分布规律,运用电荷耦合组件(Charge Coupled Device,CCD)成像技术结合激光扫描及计算机图像处理技术,在水和甘油混合液中对无黏性均匀玻璃颗粒的起动拖曳力进行测量,并根据层流绕流理论建立了无黏性均匀颗粒的滚动起动模型,提出层流条件下无黏性均匀颗粒的量纲一起动拖曳力公式.试验观测发现在层流范围内希尔兹数呈规则的带状分布,流体作用引起颗粒床表面粗化并导致床面颗粒突起减小,使颗粒起动拖曳力增大1倍以上.结果表明颗粒床表面结构性状直接影响着颗粒起动拖曳力的大小,在层流区Shields曲线具有带状分布特性.     

波流边界层水沙运动数值模拟——I:水动力模拟

波流边界层水动力模拟对研究波流相互作用和泥沙运动具有重要的理论意义和实践价值。开发了波流边界层1DV垂向一维水动力数值模型,可用于模拟漩涡沙波床面和平底床面水动力特征。模型的构建基于边界层控制方程,平底床面采用k-ε模型,沙波床面采用双层模型,提出了漩涡层和紊动扩散层交界面紊动动能和紊动耗散率表达式。试验资料验证表明,模型较好地模拟了波浪-水流-床面共同作用下的边界层水动力特征,包括波周期内不同相位流速分布、紊动动能、剪切应力等以及波致时均流速分布和波流相互作用下的时均流速分布等。根据所建模型,讨论了不同床面和波流组合条件下的水动力特征。该模型可为研究波流边界层内水动力特征提供工具。     

单向明渠流与波浪作用下植被对水沙运动影响研究综述

作为河、湖以及滨海湿地生态系统中必不可少的组成部分,水生植被具有重要的生态服务价值,且许多生态服务价值是通过改变水体动力条件实现的。含植被水流研究不仅可用于科学阐明水生植被的生态环境效应,还能指导河湖生态系统修复及污染治理的工程实践。本文考虑单向明渠流与波浪2种水动力环境,对国内外有关水生植被对水流结构以及泥沙运动影响研究的主要成果进行梳理。单向明渠流条件下,植被对水动力的影响研究主要集中于植被对水流阻力的影响以及冠层内水体的紊动结构与紊动尺度特征;波浪条件下,植被对波高与波浪流速的减弱作用以及冠层内水体的时均与紊动结构特征是研究重点。受水动力条件控制,植被影响下的泥沙运动特征也受到广泛关注,且研究焦点为单向明渠流条件下水生植被对泥沙起动与输移的影响以及波浪条件下植被对床底泥沙再悬浮的影响。     

波浪作用下的床面切应力研究进展

波浪作用下的床面切应力是估算近岸泥沙起动、输移的重要基本参数之一,其形成的底部摩阻效应也会对近岸水动力环境产生影响。由于现场观测较为困难,对波浪边界层与波浪作用下床面切应力的认识主要建立在室内试验观测基础上。回顾国内外相关理论模型和试验测量研究,梳理各类研究方法与测量技术的优缺点及适用条件;整理20世纪70年代至今的大量试验资料,对已有研究成果进行归纳分析,包括波浪非线性、波浪破碎等因素对床面切应力的影响;总结现有研究存在的局限性,提出今后的研究重点。超大型水槽是未来突破波浪边界层理论研究瓶颈的重要设施;高适应性水下二维切应力传感器的发展是复杂动力条件下床面切应力研究取得突破的关键。     

考虑渗流力的海床临界冲刷机理及计算方法

从波浪引起的海床内部渗流与床面泥沙运动耦合的角度出发,研究波浪作用下海床临界冲刷机理及计算方法。研究表明,波浪作用下床面泥沙起动及冲刷是一个不断向下发展的过程,最终达到一个临界冲刷深度。波浪引起的渗流力能够显著降低泥沙的临界起动切应力,促进泥沙起动,是影响海床冲刷的一个重要因素。将渗流力引入到传统泥沙起动公式中,推导并给出了波浪作用下海床临界冲刷深度的计算方法。结合室内和现场两个算例,很好地解释了波浪水槽试验中海床"流化"现象和黄河水下三角洲粉砂流冲沟等灾害地貌特征及成因,初步验证了该方法用于评价和计算海床冲刷的有效性。     

波流边界层水沙运动数值模拟——Ⅲ.底部高含沙层模拟

为了解析底部高含沙层特征,研究了物理影响机制和时均含沙量剖面。采用波流边界层1DV泥沙数学模型,结合试验资料,通过设置不同的计算工况进行了影响因素敏感性分析;在此基础上,提出波浪相关的时均泥沙扩散系数分布,考虑主要影响机制,推导了基于边界层物理过程的波浪作用下时均含沙量剖面表达式。结果表明,底部高含沙层与波浪边界层密切相关,是受水动力和床面形态综合影响的结果。仅建立高含沙层与水动力或床面形态的单一关系是有局限性的。含沙量层化效应和制约沉速对底部高含沙层具有重要影响。提出的平底床面含沙量剖面表达式为幂函数-Rouse-指数分布,漩涡沙波床面为指数-幂函数-Rouse分布。预期可应用于二维和三维泥沙数值模拟。     

排海污水挟沙射流研究进展

排海污水中往往含有一定数量的悬浮物, 这些悬浮物在排放口附近中呈现出复杂的扩散与沉降特性, 具有挟沙射流的典型运动特征。本文系统回顾了静水、恒定流、潮流和波浪等环境下排海污水挟沙射流的研究进展, 发现在不同动力环境下挟沙射流泥沙沉积形态存在显著差异性。如何精细刻画波浪-水流-射流多动力耦合过程及其对泥沙输运和沉降过程的影响是排海污水挟沙射流研究的难点所在, 对此提出了基于PIV/PLIF固相分离物理试验、粒子追踪两相流数值模型和拉格朗日积分模型, 系统开展波流共同作用下挟沙射流运动规律研究的设想, 以期为排海污水工程的优化设计提供重要参考。     

Sedimentary processes in the Selvage sediment-wave field, NE Atlantic: new insights into the formation of sediment waves by turbidity currents

An integrated geophysical and sedimentological investigation of the Selvage sediment-wave field has revealed that the sediment waves are formed beneath unconfined turbidity currents. The sediment waves occur on the lower continental rise and display wavelengths of up to 1 km and wave heights of up to 6 m. Wave sediments consist of interbedded turbidites and pelagic/hemipelagic marls and oozes. Nannofossil-based dating of the sediments indicates a bulk sedimentation rate of 2·4 cm 1000 years^–1, and the waves are migrating upslope at a rate of 0·28 m 1000 years^–1. Sediment provenance studies reveal that the turbidity currents maintaining the waves are largely sourced from volcanic islands to the south. Investigation of existing models for sediment-wave formation leads to the conclusion that the Selvage sediment waves form as giant antidunes. Simple numerical modelling reveals that turbidity currents crossing the wave field have internal Froude numbers of 0·5–1·9, which is very close to the antidune existence limits. Depositional flow velocities range from <6 to 125 cm^–1. There is a rapid increase in wavelength and flow thickness in the upper 10 km of the wave field, which is unexpected, as the slope angle remains relatively constant. This anomaly is possibly linked to a topographic obstacle just upslope of the sediment waves. Flows passing over the obstacle may undergo a hydraulic jump at its boundary, leading to an increase in flow thickness. In the lower 15 km of the wave field, flow thickness decreases downslope by 60%, which is comparable with results obtained for other unconfined turbidity currents undergoing flow expansion.     

深水牵引流形成的床形单元组合

海底上发育深水牵引流形成的各种床形单元,包括等深岩丘及大型沉积物波在内的不同床形单元有规律地组合在一起,对于此类床形组合的确切形成机理,目前仍在探讨之中。本文经详细研究得出以下几点进展：① 运用内波理论可对海底上大型沉积物波各组成单元的成因作出较为合理的解释,向深海方向传播的内波可形成向上坡方向迁移的大型沉积物波;② 在等深流与上覆低密度水体之间的界面上具备产生大规模界面内波的条件;③ 某些底流成因的床形单元组合属于等深流与等深流所引发内波的联合作用的产物,首次提出了等深流－内波沉积组合的概念;④ 在现代海底上及古代地层记录中均发现了等深流－内波沉积组合的实例;⑤ 建立了一个古代地层记录中的深水牵引流沉积组合综合模式。     

Threshold studies of gravel size particles under the co-linear combined action of waves and currents

Gravel size sediment beds are tested under the combined influence of simulated wave action and co-linear currents in a laboratory flume. Critical current speed, at threshold, increases with increasing size. Superimposed wave energy causes a small reduction in the unidirectional current energy. For low values of wave-induced near bed current velocities, the resistance to erosion increases when the wave period decreases from 10 to 6 s. Finally, combined critical shear stresses are found to be lower than those predicted using the Shields curve, as modified for oscillatory flow. Grain protrusion is suggested as a mechanism to explain this divergence.     

Sediment suspension dynamics and a new criterion for the maintenance of turbulent suspensions

The vertical component of the turbulent flow acceleration term, , is used to determine the net positive vertical force that may support a suspended sediment load. A dimensionless criterion, Λ, is proposed for the maintenance of suspension, defined as the ratio of the maximum vertical turbulent stress to immersed weight of the suspended load above a unit bed area. In order that a suspension be maintained: where v ′ is instantaneous vertical turbulent velocity, σ and ρ are solid and fluid densities, respectively and m is the suspended load dry mass. The Λ criterion is dynamic, being a ratio of stresses and is analogous in this respect to Shields dimensionless stress criterion, θ, for the initiation of bedload motion. The new criterion is successful in predicting the maintenance of steady-state suspended sediment transport in open channel shear flow and deposition from non-uniform particulate density flows of wall jet type.     

Review of the Seabed Sediment Resuspension by Internal Solitary Wave

Internal Solitary Waves (ISWs) are nonlinear, large amplitude motions of the interface between two fluids with different densities in the stratified ocean. Because of their strong vertical and horizontal current velocity, and the vortex, turbulent mixing caused by breaking, they affect marine environment, seabed sediment and man-made structures in the ocean. In the paper, we systematically analyzed and summarized the ISW-induced shoaling break mechanisms, models of suspension, and seabed dynamical response. Then, we discussed the ISW-induced sediment resuspension criteria, forming bottom and intermediate nepheloid layer and the capacity to suspend sediments in the seabed, and further put forward the unsolved problems based on the conducted work and related achievements. In shallow seas with complex terrain variations, shoaling can cause ISWs to deform, break, and split. Studies on the propagation of ISWs of depression over sloping topography have shown that an adverse pressure gradient causes the rotation of the flow separation, which produces vortices, and this results in global instability of the boundary layer and ISW burst. The separation vortices increase the bottom shear stress, vertical velocity, and near-bottom Reynolds stress, which leads to sediment resuspension and transport in the flow and vortex core. Although episodic, ISW-induced resuspension is hypothesized to be important enough to shape the topography. Shoaling ISWs may erode, resuspend and transport mud-like sediments, first towards shore by boluses, and subsequently offshore through the generation of intermediate nepheloid layers. Shoaling ISWs might be an important mechanism of muddy sediment dispersal along continental shelves. Furthermore, recent hypotheses suggest that sediment mobilization and transport caused by internal waves in general, and ISWs in particular, may be at the origin of some sedimentary structures found in the sedimentary rock record and also the hummocky-cross stratification. Observed on-shelf propagating frontal ISW most likely interacts with the sand waves, sediment waves or sand dunes. ISWs contribute to their generation, as they are trailed by considerable shear-induced turbulence and high-frequency internal waves close to the buoyancy frequency. This work is of great value for further understanding the process of ISW-induced sediment resuspension, transportation, and the capacity to suspend sediments in the seabed. It helps further study of the dynamic process of marine ecological environment dynamic process by ISW and the deep sea sedimentation process.