黄河宁蒙段河道水沙输移过程的仿真分析

通过建立河道水沙仿真系统模型,系统的模拟了天然状态下宁蒙河道泥沙的输移过程,并与现状条件下该河段泥沙输移过程进行比较.结果表明:龙羊峡和刘家峡水库建成运行后,由于水库拦蓄作用,宁蒙河段入口下河沿站1969-2003年累计减少输沙量25.9966×108t,说明上游梯级水库的蓄水拦沙作用显著地减轻了宁蒙河段的淤积过程.宁...     

小浪底水库“腾库迎洪”期异重流形成分析

小浪底水库异重流的形成与运行规律是关系黄河调水调沙方案拟定的重要指标。基于现场调研、理论探讨与预测分析结合等手段,明晰了2018年“腾库迎洪”期小浪底水库来水来沙过程,根据水库调度与库区淤积形态等边界条件,观测了异重流的形成与运行,采用基于流速分布的水库异重流潜入点预测公式进行了计算分析。结果表明:2018年7月库区产生的剧烈冲刷,在潜入点处发生大量淤积,水深大幅减小,引起潜入点下移;潜入点位置随入库水沙增大、坝前水位降低及库区淤积的推进向坝前移动,计算与实测结果反映了该变化过程。     

小浪底水库“腾库迎洪”期异重流形成分析

小浪底水库异重流的形成与运行规律是关系黄河调水调沙方案拟定的重要指标。基于现场调研、理论探讨与预测分析结合等手段,明晰了2018年"腾库迎洪"期小浪底水库来水来沙过程,根据水库调度与库区淤积形态等边界条件,观测了异重流的形成与运行,采用基于流速分布的水库异重流潜入点预测公式进行了计算分析。结果表明:2018年7月库区产生的剧烈冲刷,在潜入点处发生大量淤积,水深大幅减小,引起潜入点下移;潜入点位置随入库水沙增大、坝前水位降低及库区淤积的推进向坝前移动,计算与实测结果反映了该变化过程。     

黄河下游河道悬移质泥沙与床沙交换计算研究

在前人研究成果的基础上,根据挟沙水流任一粒径组泥沙在输移过程中质量守恒原理,建立了一维非恒定挟沙水流悬移质泥沙和床沙交换基本方程;通过引入平衡冲淤物粒径的概念,建立了河床处于淤积与冲刷时冲淤物粒径的计算公式,并提出了一套完整的一维非恒定挟沙水流悬移质泥沙和床沙交换计算方法。将该成果引入黄河下游一维扩展泥沙数学模型中,采用黄河下游1977年高含沙洪水与1999年汛后至2002年汛前冲刷系列进行了验证。结果表明,该方法能较好地模拟悬沙与床沙的交换过程,克服分组挟沙力方法的缺陷,使得非均匀沙计算理论上更加完善,应用上更加方便。     

水库浑水异重流潜入点判别条件

泥沙淤积是影响多沙河流水库寿命的一大难题,而异重流排沙是减少库区淤积的重要措施之一。异重流的潜入现象是异重流开始形成的直观标志,研究异重流潜入条件的判别方法有助于掌握异重流在库区内的演进规律。总结了水库异重流潜入条件的定性描述及定量计算方法,指出已有的潜入点判别公式的优缺点及适用范围,改进了描述异重流运动的动量方程,同时分析了异重流流速与含沙量沿垂线不均匀分布对动量传递的影响;在此基础上提出新的异重流潜入条件判别式,并用多组室内及野外实测资料对该判别条件进行率定与验证。分析结果表明,新的计算公式可用于判别小浪底库区异重流的潜入条件。     

浊流和泥石流的异重流初期潜入点的实验研究

潜入点的水流泥沙条件是异重流的发生条件和持续条件,受到了国内外学者的广泛关注。初期潜入点Fr代表异重流的发生条件,而稳定潜入点Fr则代表异重流的持续条件。通过一系列的低浓度浊流和高浓度浊流及泥石流的异重流潜入点的实验研究,分析对比低浓度浊流和高浓度浊流及泥石流的异重流的关系,得到在均匀顺直水槽中的异重流初期的潜入点Fr规律。提出在一定的水槽宽度、泥沙浓度条件下,异重流的初期潜入点Fr与头部流速成正比。与其它文献不同研究条件的实验结果对比有很好的一致性,说明该研究结论有很好的可靠性。该研究结论的适用范围为初期潜入点,不适用于稳定潜入点和初期潜入点与稳定潜入点之间的过渡阶段。但高浓度浊流和泥石流的异重流初期潜入点Fr与其稳定潜入点Fr很接近,粘性泥石流的异重流更接近。     

淤泥质河口水沙运动数学模型相关问题

为研究淤泥质河口的水沙运动规律,建立了用于模拟淤泥质河口水沙运动的二维数学模型。该模型采用基于无结构三角网格下的有限体积法对方程组进行离散,结合Roe-MUSCL方法及时间方向的预测-校正格式,使模型在时空方向具有二阶计算精度。模型中分别采用不同方法计算粘性和非粘性泥沙的输移源项,并引入粘性泥沙的起动流速和冲刷率计算公式。采用已有的概化水槽试验数据对模型进行了初步验证。然后模拟了1995年10月小潮及大潮期间海河口的潮流运动与泥沙输移过程,计算得到的潮位、潮流速及含沙量过程与实测过程符合较好,结果表明模型能够用来模拟淤泥质河口粘性和非粘性泥沙的不平衡输移过程。同时还比较了泥沙输移源项的不同处理方式对计算结果的影响,计算表明在淤泥质河口水沙运动数学模型中必须同时考虑粘性和非粘性泥沙的输移。     

斜坡上异重流的三维数值模拟

针对异重流的流动特征,建立了适用于具有各向异性浮力紊动特征的三维异重流运动的数学模型,并模拟了异重流在15°斜坡底面上的潜行过程。计算结果准确地模拟了异重流的运动特征和形态,其前锋的潜行速度与实验结果相当吻合。该模型采用非结构同位网格上的SIMPLEC算法能适应复杂边界和地形,可应用于自然界实际环境中异重流的演进计算。     

黄河下游高含沙洪水过程一维水沙耦合数学模型

采用浑水控制方程,建立了基于耦合解法的一维非恒定非均匀沙数学模型,用于模拟高含沙洪水演进时的河床冲淤过程.然后采用黄河下游游荡段1977年7—8月实测高含沙洪水资料对该模型进行率定,基于水沙耦合解法的各水文断面流量、总含沙量及分组含沙量的计算过程与实测过程符合更好,计算的沿程最高水位及累计河段冲淤量与实测值也较为符合.最后还采用2004年8月高含沙洪水资料对该模型进行了验证.模型率定及验证计算结果表明,采用一维水沙耦合模型计算高含沙洪水过程,能取得较高的精度.     

黄河三门峡水库水沙电耦合模型研究及应用

为实现多沙河流水库短期兴利效益和长期减淤效益的平衡,采用数学模型研究水库减淤与发电综合效益具有重要意义。充分考虑水沙输移和河床变形的耦合作用,通过结合一维水沙计算模块、水库调度模块以及发电计算模块,建立了考虑水库调度的水沙电耦合计算模型。将该模型应用到三门峡水库,采用2019水文年和2020水文年实测水沙资料分别对模型进行率定及验证。并应用该模型研究了不同水沙条件和调度方式对水库冲淤和发电的影响。计算结果表明：各测站流量、水位以及含沙量的计算过程与实测过程符合较好,发电量计算值误差均在15%以内,该模型能较为准确地模拟库区冲淤及电站发电过程;较丰水枯沙过程,其余3种典型水沙条件下冲淤量和发电量变幅分别达129%～360%及18%～52%,受入库水沙条件影响大;降低非汛期起调水位有利于控制非汛期库区淤积,提高汛期敞泄的临界入库流量有利于提高减淤及发电综合效益。     

Turbidite depositional patterns and flow characteristics, Navy Submarine Fan, California Borderland

The late Pleistocene and Holocene stratigraphy of Navy Fan is mapped in detail from more than 100 cores. Thirteen ¹⁴C dates of plant detritus and of organic-rich mud beds show that a marked change in sediment supply from sandy to muddy turbidites occurred between 9000 and 12,000 years ago. They also confirm the correlation of several individual depositional units. The sediment dispersal pattern is primarily controlled by basin configuration and fan morphology, particularly the geometry of distributary channels, which show abrupt 60° bends related to the Pleistocene history of lobe progradation. The Holocene turbidity currents are depositing on, and modifying only slightly, a relict Pleistocene morphology. The uppermost turbidite is a thin sand to mud bed on the upper-fan valley levées and on parts of the mid-fan. Most of its sediment volume is in a mud bed on the lower fan and basin plain downslope from a sharp bend in the mid-fan distributary system. Little sediment occurs farther downstream within this distributary system. It appears that most of the turbidity current overtopped the levée at the channel bend, a process referred to as flow stripping. The muddy upper part of the flow continued straight down to the basin plain. The residual more sandy base of the flow in the distributary channel was not thick enough to maintain itself as gradient decreased and the channel opened out on to the mid-fan lobe. Flow stripping may occur in any turbidity current that is thick relative to channel depth and that flows in a channel with sharp bends. Where thick sandy currents are stripped, levée and mid-fan erosion may occur, but the residual current in the channel will lose much of its power and deposit rapidly. In thick muddy currents, progressive overflow of mud will cause less declaration of the residual channelised current. Thus both size and sand-to-mud ratio of turbidity currents feeding a fan are important factors controlling morphologic features and depositional areas on fans. The size-frequency variation for different types of turbidity currents is estimated from the literature and related to the evolution of fan morphology.     

Patterns of deposition from experimental turbidity currents with reversing buoyancy

Turbidity currents are turbulent, sediment‐laden gravity currents which can be generated in relatively shallow shelf settings and travel downslope before spreading out across deep‐water abyssal plains. Because of the natural stratification of the oceans and/or fresh water river inputs to the source area, the interstitial fluid within which the particles are suspended will often be less dense than the deep‐water ambient fluid. Consequently, a turbidity current may initially be denser than the ambient sea water and propagate as a ground‐hugging flow, but later reverse in buoyancy as its bulk density decreases through sedimentation to become lower than that of the ambient sea water. When this occurs, all or part of the turbidity current lofts to form a buoyant sediment‐laden cloud from which further deposition occurs. Deposition from such lofting turbidity currents, containing a mixture of fine and coarse sediment suspended in light interstitial fluid, is explored through analogue laboratory experiments complemented by theoretical analysis using a ‘box and cloud’ model. Particular attention is paid to the overall deposit geometry and to the distributions of fine and coarse material within the deposit. A range of beds can be deposited by bimodal lofting turbidity currents. Lofting may encourage the formation of tabular beds with a rapid pinch‐out rather than the gradually tapering beds more typical of waning turbidity currents. Lofting may also decouple the fates of the finer and coarser sediment: depending on the initial flow composition, the coarse fraction can be deposited prior to or during buoyancy reversal, while the fine fraction can be swept upwards and away by the lofting cloud. An important feature of the results is the non‐uniqueness of the deposit architecture: different initial current compositions can generate deposits with very similar bed profiles and grading characteristics, highlighting the difficulty of reconstructing the nature of the parent flow from field data. It is proposed that deposit emplacement by lofting turbidity currents is common in the geological record and may explain a range of features observed in deep‐water massive sands, thinly bedded turbidite sequences and linked debrites, depending on the parent flow and its subsequent development. For example, a lofting flow may lead to a well sorted, largely ungraded or weakly graded bed if the fines are transported away by the cloud. However, a poorly sorted, largely ungraded region may form if, during buoyancy reversal, high local concentrations and associated hindered settling effects develop at the base of the cloud.     

Bedload transport and bed resistance associated with density and turbidity currents

Turbidity currents in the ocean are driven by suspended sediment. Yet results from surveys of the modern sea floor and turbidite outcrops indicate that they are capable of transporting as bedload and depositing particles as coarse as cobble sizes. While bedload cannot drive turbidity currents, it can strongly influence the nature of the deposits they emplace. This paper reports on the first set of experiments which focus on bedload transport of granular material by density underflows. These underflows include saline density flows, hybrid saline/turbidity currents and a pure turbidity current. The use of dissolved salt is a surrogate for suspended mud which is so fine that it does not settle out readily. Thus, all the currents can be considered to be model turbidity currents. The data cover four bed conditions: plane bed, dunes, upstream‐migrating antidunes and downstream‐migrating antidunes. The bedload transport relation obtained from the data is very similar to those obtained for open‐channel flows and, in fact, is fitted well by an existing relation determined for open‐channel flows. In the case of dunes and downstream‐migrating antidunes, for which flow separation on the lee sides was observed, form drag falls in a range that is similar to that due to dunes in sand‐bed rivers. This form drag can be removed from the total bed shear stress using an existing relation developed for rivers. Once this form drag is subtracted, the bedload data for these cases collapse to follow the same relation as for plane beds and upstream‐migrating antidunes, for which no flow separation was observed. A relation for flow resistance developed for open‐channel flows agrees well with the data when adapted to density underflows. Comparison of the data with a regime diagram for field‐scale sand‐bed rivers at bankfull flow and field‐scale measurements of turbidity currents at Monterey Submarine Canyon, together with Shields number and densimetric Froude number similarity analyses, provide strong evidence that the experimental relations apply at field scale as well.     

Thresholds of intrabed flow and other interactions of turbidity currents with soft muddy substrates

Controlled laboratory experiments reveal that the lower part of turbidity currents has the ability to enter fluid mud substrates, if the bed shear stress is higher than the yield stress of the fluid mud and the density of the turbidity current is higher than the density of the substrate. Upon entering the substrate, the turbidity current either induces mixing between flow‐derived sediment and substrate sediment, or it forms a stable horizontal flow front inside the fluid mud. Such ‘intrabed’ flow is surrounded by plastically deformed mud; otherwise it resembles the front of a ‘bottom‐hugging’ turbidity current. The ‘suprabed’ portion of the turbidity current, i.e. the upper part of the flow that does not enter the substrate, is typically separated from the intrabed flow by a long horizontal layer of mud which originates from the mud that is swept over the top of the intrabed flow and then incorporated into the flow. The intrabed flow and the mixing mechanism are specific types of interaction between turbidity currents and muddy substrates that are part of a larger group of interactions, which also include bypass, deposition, erosion and soft sediment deformation. A classification scheme for these types of interactions is proposed, based on an excess bed shear stress parameter, which includes the difference in the bed shear stress imposed by the flow and the yield stress of the substrate and an excess density parameter, which relies on the density difference between the flow and the substrate. Based on this classification scheme, as well as on the sedimentological properties of the laboratory deposits, an existing facies model for intrabed turbidites is extended to the other types of interaction involving soft muddy substrates. The physical threshold of flow‐substrate mixing versus stable intrabed flow is defined using the gradient Richardson number, and this method is validated successfully with the laboratory data. The gradient Richardson number is also used to verify that stable intrabed flow is possible in natural turbidity currents, and to determine under which conditions intrabed flow is likely to be unstable. It appears that intrabed flow is likely only in natural turbidity currents with flow velocities well below ca 3·5 m s^?1, although a wider range of flows is capable of entering fluid muds. Below this threshold velocity, intrabed flow is stable only at high‐density gradients and low‐velocity gradients across the upper boundary of the turbidity current. Finally, the gradient Richardson number is used as a scaling parameter to set the flow velocity limits of a natural turbidity current that formed an inferred intrabed turbidite in the deep‐marine Aberystwyth Grits Group, West Wales, United Kingdom.     

Linking submarine channel–levee facies and architecture to flow structure of turbidity currents: insights from flume tank experiments

Submarine leveed channels are sculpted by turbidity currents that are commonly highly stratified. Both the concentration and the grain size decrease upward in the flow, and this is a fundamental factor that affects the location and grain size of deposits around a channel. This study presents laboratory experiments that link the morphological evolution of a progressively developing leveed channel to the suspended sediment structure of the turbidity currents. Previously, it was difficult to link turbidity current structure to channel–levee development because observations from natural systems were limited to the depositional products while experiments did not show realistic morphodynamics due to scaling issues related to the sediment transport. This study uses a novel experimental approach to overcome scaling issues, which results in channel inception and evolution on an initially featureless slope. Depth of the channel increased continuously as a result of levee aggradation combined with varying rates of channel floor aggradation and degradation. The resulting levees are fining upward and the grain‐size trend in the levee matches the upward decrease in grain size in the flow. It is shown that such deposit trends can result from internal channel dynamics and do not have to reflect upstream forcing. The suspended sediment structure can also be linked to the lateral transition from sediment bypass in the channel thalweg to sediment deposition on the levees. The transition occurs because the sediment concentration is below the flow capacity in the channel thalweg, while higher up on the channel walls the concentration exceeds capacity resulting in deposition of the inner levee. Thus, a framework is provided to predict the growth pattern and facies of a levee from the suspended sediment structure in a turbidity current.     

Turbidite stacking patterns in salt‐controlled minibasins: Insights from integrated analogue models and numerical fluid flow simulations

The sea floor of intraslope minibasins on passive continental margins plays a significant role in controlling turbidity current pathways and the resulting sediment distribution. To address this, laboratory analogue modelling of intraslope minibasin formation is combined with numerical flow simulations of multi‐event turbidity currents. This approach permits an improved understanding of evolving flow–bathymetry–deposit interactions and the resulting internal stacking patterns of the infills of such minibasins. The bathymetry includes a shelf to slope channel followed by an upper minibasin, which are separated by a confining ridge from two lower minibasins that compares well with analogous bathymetries reported from natural settings. From a wider range of numerical flow experiments, a series of 100 consecutive flows is reported in detail. The turbidity currents are released into the channel and upon reaching the upper minibasin follow a series of stages from short initial ponding, ‘filling and spilling’ and an extended transition to long retrogradational ponding. Upon reaching the upper minibasin floor, the currents undergo a hydraulic jump and therefore much sediment is deposited in the central part of the minibasin and the counterslope. This modifies the bathymetry such that in the fill and spill stage, flow stripping and grain‐size partitioning cause some finer sediment to be transported across the confining ridge into the lower minibasins. Throughout the basin infill process, the sequences retrograde upstream, accompanied by lateral switching into locally formed depressions in the upper minibasin. After the fill and spill stage, significant deposition occurs in the channel where retrograding cyclic steps with wavelengths of 1 to 2 km develop as a function of pulsating flow criticality. These results are at variance with conventional schemes that emphasize sequential downstream minibasin filling through ponding dominated by vertical aggradation. Comparison of these results with published field and experimental examples provides support for the main conclusions.     

Deposits of depletive high-density turbidity currents: a flume analogue of bed geometry, structure and texture

Flume experiments were performed to study the flow properties and depositional characteristics of high‐density turbidity currents that were depletive and quasi‐steady to waning for periods of several tens of seconds. Such currents may serve as an analogue for rapidly expanding flows at the mouth of submarine channels. The turbidity currents carried up to 35 vol.% of fine‐grained natural sand, very fine sand‐sized glass beads or coarse silt‐sized glass beads. Data analysis focused on: (1) depositional processes related to flow expansion; (2) geometry of sediment bodies generated by the depletive flows; (3) vertical and horizontal sequences of sedimentary structures within the sediment bodies; and (4) spatial trends in grain‐size distribution within the deposits. The experimental turbidity currents formed distinct fan‐shaped sediment bodies within a wide basin. Most fans consisted of a proximal channel‐levee system connected in the downstream direction to a lobe. This basic geometry was independent of flow density, flow velocity, flow volume and sediment type, in spite of the fact that the turbidity currents of relatively high density were different from those of relatively low density in that they exhibited two‐layer flow, with a low‐density turbulent layer moving on top of a dense layer with visibly suppressed large‐scale turbulence. Yet, the geometry of individual morphological elements appeared to relate closely to initial flow conditions and grain size of suspended sediment. Notably, the fans changed from circular to elongate, and lobe and levee thickness increased with increasing grain size and flow velocity. Erosion was confined to the proximal part of the leveed channel. Erosive capacity increased with increasing flow velocity, but appeared to be constant for turbidity currents of different grain size and similar density. Structureless sediment filled the channel during the waning stages of the turbidity currents laden with fine sand. The adjacent levee sands were laminated. The massive character of the channel fills is attributed to rapid settling of suspension load and associated suppression of tractional transport. Sediment bypassing prevailed in fan channels composed of very fine sand and coarse silt, because channel floors remained fully exposed until the end of the experiments. Lobe deposits, formed by the fine sand‐laden, high‐density turbidity currents, contained massive sand in the central part grading to plane parallel‐laminated sand towards the fringes. The depletive flows produced a radial decrease in mean grain size in the lobe deposits of all fans. Vertical trends in grain size comprised inverse‐to‐normal grading in the levees and in the thickest part of the lobes, and normal grading in the channel and fringes of the fine sandy fans. The inverse grading is attributed to a process involving headward‐directed transport of relatively fine‐grained and low‐concentrated fluid at the level of the velocity maximum of the turbidity current. The normal grading is inferred to denote the waning stage of turbidity‐current transport.     

Modelling of turbidity currents on Navy Submarine Fan, California Continental Borderland

Several Holocene turbidites can be correlated across much of Navy Fan through more than 100 sediment core localities. The uppermost muddy turbidite unit is mapped throughout the northern half of the fan; its volume, grain-size distribution and the maximum height of deposition on the basin slopes are known. These parameters can be related to the precise channel morphology and mesotopography revealed by deep-tow surveys. Thus there is sufficient information to estimate detailed flow characteristics for this turbidity current as it moved from fan valley to distal basin plain. On the upper fan, the gradient and the increasing downstream width of the channel and only limited flow overspill suggest that the flow had a Froude number close to 1.0. The sediment associated with the channel indicates friction velocities of about 0.06 m s^?1 and flow velocities of about 0.75 m s^?1. Using this flow velocity and channel dimensions, sediment concentration (～2×10^?3) and discharge are estimated, and from a knowledge of the total volume of sediment deposited, the flow duration is estimated to be from 2 to 9 days. It is shown that the estimates of Froude number, drag coefficient, and sediment concentration are not likely to vary by more than a factor of 2. On the mid-fan, the flow was much thicker than the height of the surface relief of the fan and it spread rapidly. The cross-flow slope, determined from the horizontal extent of turbidite sediment, is used to estimate flow velocity, which is confirmed by consideration of both sediment grain size and rate of deposition. This again allows sediment concentration and discharge to be estimated. The requirements of flow continuity, entrainment of water during flow expansion, and observed sediment deposition provide checks on all these estimates, and provide an integrated picture of the evolution of the flow. The flow characteristics of this muddy turbidity current are well constrained compared to those for more sand-rich late Pleistocene and early Holocene turbidity currents on the fan.     

渤海湾盆地东营凹陷碎屑流主控型深水体系沉积过程及模式

【研究目的】碎屑流是深水环境沉积物搬运和分散的重要机制,其相关的砂岩储层是含油气盆地重要的勘探目标,然而,与经典浊流及浊积系统相比,对碎屑流主控型深水体系的发育规律目前仍知之甚少。【研究方法】本文基于岩心、测井及全三维地震资料,通过系统的岩心观察描述、测井及地震资料解释,对渤海湾盆地东营凹陷始新统沙三中亚段深水体系沉积过程及模式开展研究。【研究结果】结果表明,沙三中深水体系发育九种异地搬运岩相,可概括为四大成因类型,反映了块体及流体两种搬运过程。岩相定量统计表明,该深水体系主要由碎屑流沉积构成,浊流沉积很少,碎屑流中又以砂质碎屑流为主。重力流在搬运过程中经历了滑动、滑塌、砂质碎屑流、泥质碎屑流及浊流等5个阶段演变,发育5类主要的深水沉积单元,包括滑动体、滑塌体、碎屑流水道、碎屑流朵体及浊积薄层砂。从发育规模及储层物性上,砂质碎屑流水道、朵体及砂质滑动体构成了本区最重要的深水储层类型。【结论】认为沙三中时期充足的物源供给、三角洲前缘高沉积速率、断陷期频繁的断层活动以及较短的搬运距离是碎屑流主控型深水体系形成及演化的主控因素,最终基于沉积过程、沉积样式及盆地地貌特征综合建立了碎屑流主控型深水体系沉积模式。本研究将进一步丰富深水沉积理论,为陆相深水储层预测提供借鉴。