Numerical simulation of turbidity current flow and sedimentation: II. Results and geological applications

A process-based, forward computer model of turbidity current flow and sedimentation, termed the TCFS model, has been developed to trace the downslope evolution of individual turbidity flows. Details of the model itself have been presented in a preceding paper. We here outline a series of tests of the TGFS model. The sensitivity tests of the TCFS model to general geological controls reveal the quantitative relationship between these controls and the behaviour of turbidity flows and the geometry and textural features of the resulting turbidites. Experimental turbidity currents on relatively steep slopes accelerate more rapidly and reach higher velocities than those on gentle slopes. Flows with larger initial volumes have higher initial velocities, travel further downslope, and form beds of greater thickness and downslope extent than smaller flows. Experimental high-concentration flows with suspended-sediment concentrations of 25% accelerate more rapidly and reach higher downslope velocities than dilute flows with 5% suspended sediment. The higher velocities and enhanced hindered-settling effects of the high-concentration flows lead to much greater transport distances and reduced vertical and lateral sediment size grading in the resulting turbidites. Beds formed by experimental high-concentration flows are massive or show coarse-tail grading whereas beds formed by low-concentration flows show distribution-grading. Experimental flows fed by coarse sediment sources tend to deposit the bulk of their suspended sediment loads on the proximal slope, resulting in more rapid flow deceleration and sedimentation than flows fed by silt-rich, fine-grained sediment sources. Turbidites formed by coarse-sediment flows tend to have a wedge-shaped geometry, with low downslope extent and high surface relief, whereas turbidites formed by fine-sediment flows tend to have a tabular geometry, with greater downslope extent and lower surface relief. A specific geological test of the TCFS model is based on studies of modern turbidity currents in Bute Inlet, British Columbia, Canada. With the input initial and boundary conditions estimated from Bute Inlet, the model predicts the downslope velocity evolution of turbidity currents comparable to those of modern and ancient turbidity flows measured in Bute Inlet. Model-calculated vertical and downslope grain-size properties of turbidites are similar to those exhibited by surface and cored Bute Inlet turbidites. Model flows tend to decelerate more rapidly than some stronger turbidity currents in the Bute Inlet system, and model beds tend to decrease in grain-size downslope more rapidly than observed bottom sediments. This is probably because the TCFS model flows lacked clay, which is abundant in Bute Inlet; they do not fully simulate turbulent mixing of suspended sediments; and they better represent the unsteady, depositional stage of turbidity-currents than the preceding stage of more-or-less steady-flow conditions. These tests demonstrate that the TCFS model provides a semi-quantitative method to study the growth patterns of submarine turbidite systems. It can serve as a predictive tool for analysing the facies architecture of ancient turbidite systems through simulating multi-depositional events by improving its erosion function, and the compatibility between its numerical components.     

A numerical model for sedimentation from highly-concentrated multi-sized suspensions

The principal numerical approach to describing sedimentation from multi-component suspensions (Mirza and Richardson, 1979) has been applied only to systems containing two or three particle size populations, but is theoretically applicable to suspensions containing a wider range of particle sizes. In order to adapt this model to the simulation of sedimentation from natural high-density sediment-laden flows, we have conducted computational tests on the sedimentation of suspensions with up to ten particle size populations. The tests run smoothly with binary systems and sometimes ternary systems, but fail with systems containing more than three particle size populations. One cause of the discrepancy between theoretical predictions and computational tests arises from the method used to calculate the changing concentration of particles in each zone of a stratified settling suspension. In quantifying the changing sediment budget for each sedimentation zone, the equation set fails to include the apparent outflow of finer particle populations across the upper boundary of the zone. In the present study, we express the sediment budget of each particle population within each sedimentation zone as the net sediment flux, the algebraic sum of the apparent particle inflow from the zone's lower boundary and the apparent particle outflow across the zone's upper boundary. This revised model successfully predicts the evolution of multi-component suspensions containing up to ten particle size populations over the length of time required for complete sedimentation, up to 8×10⁵ seconds. It can be used to predict the sediment sorting and the vertical textural variation of beds formed by the simulated sedimentation of a multi-component suspension. The model provides the basis for future computer simulation of sedimentation from highly concentrated sediment flows and for the prediction of downslope textural and structural variations of turbidites.     

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.     

Response of unconfined turbidity current to deep‐water fold and thrust belt topography: Orthogonal incidence on solitary and segmented folds

Sea‐floor topography of deep‐water folds is widely considered to have a major impact on turbidity currents and their depositional systems, but understanding the flow response to such features was limited mainly to conceptual notions inspired by small‐scale laboratory experiments. High‐resolution three‐dimensional numerical experiments can compensate for the lack of natural‐scale flow observations. The present study combines numerical modelling of thrusts with fault‐propagation folds by Trishear3D software with computational fluid dynamics simulations of a natural‐scale unconfined turbidity current by MassFlow‐3D? software. The study reveals the hydraulic and depositional responses of a turbidity current (ca 50 m thick) to typical topographic features that it might encounter in an orthogonal incidence on a sea‐floor deep‐water fold and thrust belt. The supercritical current (ca 10 m sec^?1) decelerated and thickened due to the hydraulic jump on the fold backlimb counter‐slope, where a reverse overflow formed through current self‐reflection and a reverse underflow was issued by backward squeezing of a dense near‐bed sediment load. The reverse flows were re‐feeding sediment to the parental current, reducing its waning rate and extending its runout. The low‐efficiency current, carrying sand and silt, outran a downslope distance of >17 km with only modest deposition (<0·2 m) beyond the fold. Most of the flow volume diverted sideways along the backlimb to surround the fold and spread further downslope, with some overspill across the fold and another hydraulic jump at the forelimb toe. In the case of a segmented fold, a large part of the flow went downslope through the segment boundary. Preferential deposition (0·2 to 1·8 m) occurred on the fold backlimb and directly upslope, and on the forelimb slope in the case of a smaller fold. The spatial patterns of sand entrapment revealed by the study may serve as guidelines for assessing the influence of substrate folds on turbiditic sedimentation in a basin.     

海底浊流在坡道转换处的流动及沉积的数值模拟

根据一经多项试验数据验证的基于三维不可压缩流体Navier-Stokes方程和湍流 k-ε 模型的重力流数值计算的数学模型,模拟并分析了单粒径沉积物的海底浊流沿不同斜坡流至近似平坦坡的流动及沉积特征。模拟结果显示了有关海底浊流的一些重要特征:连续入流的浊流在斜坡上的流速随着斜坡的增大而增大,同时浊流厚度由于对环境水体的夹带而渐渐增厚,坡度越大,增厚越快;流至近水平坡时,流速均有明显的降低,但大斜坡入流依然保持相对较高的流速。在沉积方面,初步的模拟结果显示对给定的沉积物来说存在一相对应的临界坡度:当坡度小时,坡上沉积多,坡下少,这样整体的坡度有逐渐增大之势;当坡度大时,坡上沉积少或为侵蚀,而坡下沉积相对较多,坡度有整体减小之势。了解了不同坡度转换的浊流沉积的上述特点,对于我们根据实测的浊流沉积的剖面特征推测其形成的环境,进而推测相关油气储层的分布状况会有一定的参考作用。     

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.     

水库异重流一维水沙耦合模型

建立了基于库区不规则断面的一维非恒定异重流数学模型,并采用明流与异重流水沙输移模型交替运算的两步模式,即用潜入条件动态判别异重流计算的上游边界位置,将潜入点上游的明流浑水段与下游异重流段计算连接起来。水流运动、泥沙输移与河床变形过程完全耦合,采用TVD(Total Variation Diminishing)形式的MUSCL-Hancock格式进行数值求解。将该模型应用于恒定流量与释放定量悬沙两种条件下的异重流水槽实验模拟,比较了有无水面梯度项对模拟精度的影响,计算结果表明该模型能较为准确地预测异重流的厚度、含沙量分布及传播过程。     

Inferring the mass fraction of floc-deposited mud: application to fine-grained turbidites

Fine sediment deposition in the ocean is complicated by the cohesive nature of muds and their tendency to flocculate. The result is disaggregated inorganic grain size (DIGS) distributions of bottom sediment that are influenced by single‐grain and floc deposition. This study outlines a parametric model that characterizes bottom sediment DIGS distributions. Modelled parameters are then used to infer depositional conditions that account for the regional variation in the grain sizes deposited by turbidity currents on the Laurentian Fan–Sohm Abyssal Plain, offshore south‐eastern Canada. Results indicate that, on the channellized Laurentian Fan, the mass fraction of floc‐deposited mud increases only slightly downslope. The small evolution in this fraction arises because sediment concentration and turbulent energy are associated in turbidity currents. On the Sohm Abyssal Plain, however, the mass fraction of floc‐deposited mud decreases, probably as a result of lower sediment concentration at this source‐distal site. Estimates of the mass fraction of mud deposited as flocs suggest that floc deposition is the dominant mode by which sediment is lost from suspension, although single‐grain deposition contributes more to the depositional flux in proximal areas where high energy breaks flocs and in distal areas where low sediment concentration limits floc formation. It is concluded that, throughout the dispersal system, changes in the fraction of flocculated mud deposited from turbidity currents reflect changes in sediment concentration and energy downslope.     

扇三角洲砂体几何形态沉积数值模拟方法研究

扇三角洲是陆相碎屑岩重要的沉积类型之一.利用砂体沉积过程数值模拟方法可以预测扇三角洲砂体的几何形态.在利用泥沙冲淤动力学模式建立扇三角洲沉积过程的数学方程后,根据现代三角洲沉积特征,设计了其沉积过程的模拟条件,计算域长100 km、宽50 km,流量按50年一遇洪水设计,模拟过程到2 000年时,扇三角洲沉积过程基本达到平衡状态,此时前缘复合砂体最大厚度约39 m.通过对扇三角洲前缘砂体几何形态的模拟研究,总结出扇三角洲沉积体储层建筑结构的模拟与预测方法.计算过程中可识别出4种沉积砂体,包括水下分流河道(砂体的平均长宽比为2.68,平均宽厚比为79.8)、河口砂坝(平均长宽比为2.02,平均宽厚比为68.2)、远砂坝(平均长宽比为1.65,平均宽厚比为58.3)、水下溢岸沉积.主要沉积单元砂体几何参数之间的相关关系较好.实验结果表明,利用泥沙冲淤动力模式可以较好地揭示扇三角洲发育过程,进而可以预测扇三角洲砂体的几何形态.     

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.     

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.     

Flow structure of turbidity currents

A two‐dimensional numerical model is used to describe the flow structure of turbidity currents in a vertical plane. To test the accuracy of the model, it is applied to historical flows in Bute Inlet and the Grand Banks flow. The two‐dimensional spatial and temporal distributions of velocity and sediment concentration and non‐dimensionalized vertical profiles of velocity, turbulent kinetic energy and sediment concentration are discussed for several simple computational currents. The flows show a clear interaction between velocity, turbulence and sediment distribution. The results of the numerical tests show that flows with fine‐grained sediment have low vertical and high horizontal gradients of velocity and sediment concentration, show little increase in flow thickness and decelerate slowly. Steadiness and uniformity in these flows are comparable for velocity and concentration. In contrast, flows with coarse‐grained sediment have high vertical and low horizontal velocity gradients and high horizontal concentration gradients. These flows grow considerably in thickness and decelerate rapidly. Steadiness and uniformity in flows with coarse‐grained sediment are different for velocity and concentration. The results show the influence of spatial and temporal flow structure on flow duration and sediment transport.     

The motion of sediment-water mixtures during intense bedload transport: computer simulations

A new model, which couples fluid and particle dynamics, has been developed to study the motion of the sediment-water mixture during intense bedload transport, including the velocity profiles of both sediment and water, the roughness length of an upper plane bed and the thickness of moving sediment layers. Standard mixing length theory is used to model the motion of water above the boundary between the overlying water and the sediment-water mixture. The turbulent flow within the moving sediment layers is described by a shear stress model, in which the effective viscosity of the flowing water is proportional to the velocity difference between the fluid and the sediment. The particle dynamics method, in which the equations of motion of each of many particles are solved directly, is applied to model the movement of sediment particles. The particle-fluid interaction is expressed by a velocity-squared fluid drag force exerted on each sediment particle. Both computer simulation results and theoretical analysis have shown that the velocities of both sediment and fluid during intense sediment transport decrease exponentially with depth in the top layers of a fast-moving sediment—water mixture. The thickness of the moving sediment layers, obtained from the computer simulation results, is proportional to the shear stress, which agrees with previous experimental observations.     

Flow properties of turbidity currents in Bute Inlet, British Columbia

Bute Inlet, a fiord along the southwestern coast of British Columbia, Canada, includes a sea-floor sedimentation system 70 km in length which resembles those developed on some large submarine fans. Turbidity currents originate at the head of the flord on the submerged delta fronts of the Homathko and Southgate rivers. They move downslope for about 30 km within a single large incised channel, spill onto a depositional area termed the channel lobe complex, and finally spread out over a low-relief distal splay area that passes 55 km downslope into a flat basin floor. During the present study, turbidity currents in Bute Inlet were studied using sea-floor morphology, bottom sediment distribution, and in-situ instrument packages. The mean velocities of the most recent flows, estimated from surface sediment grain size, has varied between 100–120 cm s^–1 in the incised channel, 20–50 cms^–1 in the channel lobe complex, and < 5 cm s^–1 on the basin floor. Velocities based on channel morphology are poorly constrained but are in the range of 160-425 cm s^–1 in the upper part of the incised channel and 66 cm s^–1 in the lower channel. Calculated flow densities range from 1.049 to 1.028g cm^–3. Turbidity flows monitored in 1986 using submerged instrument packages exceeded 32 m in thickness in the upper part of the incised channel, where the maximum measured velocity was 330 cm s^–1. At the head of the channel lobe complex the maximum velocity had declined to 75 cm s^–1. The density of the monitored flows is estimated at 1.025-1.03g cm^–3. The cored sediments and channel morphology yield estimates of mean flow velocities that are generally greater than those measured by the in-situ instrument packages and estimated from modern surface sediments. The former suggest past flow velocities up to 500 cm s^–1 in the incised channel, about 20 cm s^–1 in spillover deposits along the lower part of the incised channel, and 100-140 cm s^–1 in the distal splay. The contrast between the velocities of modern and past flows suggests that past flows may have been considerably larger and more energetic than those presently occurring in Bute Inlet. The size properties of sediments in the monitored turbidity flows suggest a strong vertical size gradient in the suspended load during transport. The surface and cored sediments fine downslope from the channel lobe complex to distal splay area. Distinctive sedimentary sequences are recognized in cores from the spillover lobes, channel lobe complex, distal splay, and basin floor depositional areas. Many individual turbidites grade downslope from massive T_a divisions in the channel lobe complex and probably in the incised channel to T_a divisions overlain by slurried divisions on the distal splay and largely slurried beds on the basin floor. These facies suggest that individual currents commonly evolve from largely cohesionless suspensions in the incised channel and channel lobe complex to dilute cohesive slurries downslope on the distal splay and basin floor. Many flows in Bute Inlet fail to develop a traction state of sedimentation and the resulting turbidites lack well-developed T_b. T_c, and T_d divisions.     

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

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

河口最大浑浊带数学模拟研究的进展

分三个方面介绍了河口最大浑浊带研究中数学模拟研究的进展：（１）利用通量分析计算河口盐度、污染物和泥沙通量，分析带内泥沙富采机制；（２）由物质平衡原理建立一维和较简单的机理模型，讨论最大浑浊带的成因；（３）据水动力方程和物质平衡方程建立二维或三维数值模型，计算最大浑浊带的水流结构和物质浓度分布，模拟并探讨不同条件下最大浑浊带的成因和演化机制。     

泥沙异重流与环境物质交换经验式对比

采用异重流层平均水沙耦合数学模型,模拟开闸式和恒定入流式泥沙异重流水槽实验,对比分析异重流与环境之间物质交换经验式的适应性和不确定性。考虑4个水卷吸经验式:e_w59、e_w86、e_w87和e_w01,5个泥沙侵蚀经验式:E_s77、E_s86、E_s87、E_s93和E_s04。数值研究表明:水卷吸对于水槽异重流影响较小,应用综合考虑底床摩擦和剪切不稳定的e_w经验式时模拟结果较好;开闸式异重流对床面侵蚀能力有限。对于恒定入流式异重流,应用E_s87和E_s93侵蚀经验式计算所得淤积厚度与实测值吻合较好,可能是率定时综合考虑了异重流实验数据。     

Two-dimensional horizontal modeling of fine-sediment transport at the South Channel–North Passage of the partially mixed Changjiang River estuary,China

Tidal flow and fine-sediment transport at the South Channel–North Passage of the partially-mixed Changjiang River estuary were studied using a two-dimensional horizontal (2DH) numerical model. This 2DH model was achieved by depth-integrating the momentum and convection–diffusion equations. The Alternating Direction Implicit scheme was used to solve the governing equations. The iterative method was adopted for the calculation of convection and diffusion terms of momentum equation. Comparisons between calculated and measured results (tidal elevations and depth-averaged velocities) have shown reasonable agreement. Horizontal distributions of tidal current velocity and suspended sediment concentration were qualitatively consistent with observations. Those modeled results were analyzed to elucidate the mechanisms for the formation of the turbidity maximum and intratidal variations in fine-sediment transport processes.     

Trapping of sustained turbidity currents by intraslope minibasins

Depositional turbidity currents have filled many intraslope minibasins with sediment creating targets for petroleum exploration. The dynamics of sustained turbidity currents and their depositional characteristics are investigated in a scaled physical model of a minibasin. Each turbidity current deposited a downstream thinning wedge of sediment near the inlet. Farther downstream the turbidity current was ponded by a barrier. The ponded part of the turbidity current was separated from the sediment‐free water above by a relatively sharp, horizontal settling interface indicating highly Froude‐subcritical flow. The very slow moving flow within the ponded zone created conditions for the passive rainout of suspended sediment onto the bed. In the lower part of the ponded zone, the concentration and mean grain‐size of the sediment in suspension tended to be relatively uniform in both the vertical and streamwise directions. As a result, the deposit emplaced in the ponded zone showed only a weak tendency toward downstream fining and was passively draped over the bed in such a way that irregularities in the inerodible bed were accurately reflected. The discharge of suspended sediment overflowing the downstream end of the minibasin was significantly less than the inflow discharge, resulting in basin sediment trapping efficiencies >95%. A simple model is developed to predict the trapping of sediment within the basin based on the relative magnitudes of the input discharge of turbid water and the detrainment discharge of water across the settling interface. This model shows a limiting case in which an intraslope basin captures 100% of the sediment from a ponded turbidity current, even through a succession of sustained flow events, until sediment deposition raises the settling interface above the downstream lip of the minibasin. This same process defines one of the mechanisms for minibasin filling in nature, and, when this mechanism is operative, the trap efficiency of sediment can be expected to be high until the minibasin is substantially filled with sediment.     

Changes in fine sediment size distributions due to interactions with streambed sediments

In the classical view of fine sediment transport and deposition in streams, particles are expected to be removed from flowing water simply by direct sedimentation onto the streambed. However, recent research has demonstrated that fine sediments can propagate into pore spaces in the streambed due to hyporheic exchange and be removed by a combination of physical and chemical processes. This behaviour can significantly alter fine sediment size distributions during in-stream sediment transport because the physical transport of fine particles and their attachment to bed sediment grains are both a function of the particle size. Herein, we present model simulations for deposition of suspended sediments with a bimodal size distribution. We also applied this approach to analyse the results of laboratory flume observations of suspended sediment deposition. Results from model simulations and flume experiments clearly show that the rate of particle deposition increases with increasing particle size. Thus, the larger particles are preferentially removed from mixtures and there is a fining of the mixed suspensions over time. Both particle deposition mechanisms, i.e. particle sedimentation and filtration, contribute to the fining of the mixed fine particle suspensions over time, and their effects are clearly demonstrated using the fundamental process-based model. These results clearly demonstrate the effects of stream-subsurface exchange on the temporal evolution of the suspended fine sediment size distribution in downstream transport.