Discriminating between pore‐filling load and bed‐structure load: a new porosity‐based method,exemplified for the river Rhine

Sediments contained in the river bed do not necessarily contribute to morphological change. The finest part of the sediment mixture often fills the pores between the larger grains and can be removed without causing a drop in bed level. The discrimination between pore‐filling load and bed‐structure load, therefore, is of practical importance for morphological predictions. In this study, a new method is proposed to estimate the cut‐off grain size that forms the boundary between pore‐filling load and bed‐structure load. The method evaluates the pore structure of the river bed geometrically. Only detailed grain‐size distributions of the river bed are required as input to the method. A preliminary validation shows that the calculated porosity and cut‐off size values agree well with experimental data. Application of the new cut‐off size method to the river Rhine demonstrates that the estimated cut‐off size decreases in a downstream direction from about 2 to 0·05 mm, covariant with the downstream fining of bed sediments. Grain size fractions that are pore‐filling load in the upstream part of the river thus gradually become bed‐structure load in the downstream part. The estimated (mass) percentage of pore‐filling load in the river bed ranges from 0% in areas with a unimodal river bed, to about 22% in reaches with a bimodal sand‐gravel bed. The estimated bed porosity varies between 0·15 and 0·35, which is considerably less than the often‐used standard value of 0·40. The predicted cut‐off size between pore‐filling load and bed‐structure load (D_c,p) is fundamentally different from the cut‐off size between wash‐load and bed‐material load (D_c,w), irrespective of the method used to determine D_c,p or D_c,w. D_c,w values are in the order of 10^?1 mm and mainly dependent on the flow characteristics, whereas D_c,p values are generally much larger (about 10⁰ mm in gravel‐bed rivers) and dependent on the bed composition. Knowledge of D_c,w is important for the prediction of the total sediment transport in a river (including suspended fines that do not interact with the bed), whereas knowledge of D_c,p helps to improve morphological predictions, especially if spatial variations in D_c,p are taken into account. An alternative to using a spatially variable value of D_c,p in morphological models is to use a spatially variable bed porosity, which can also be predicted with the new method. In addition to the morphological benefits, the new method also has sedimentological applications. The possibility to determine quickly whether a sediment mixture is clast‐supported or matrix‐supported may help to better understand downstream fining trends, sediment entrainment thresholds and variations in hydraulic conductivity.     

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.     

Mechanical trapping of fine particles in a medium of mono‐sized randomly packed spheres

Mechanical trapping (or straining) of fine particles is a key mechanism in many filtration systems. For example, the performance of rapid sand filters depends in part on mechanical trapping of larger fine particles, while relying on adsorptive processes to trap very small fine particles and microbes. The ability to trap these particles is directly related to the construction of the packed bed used for filtration in this system. Thus, the ability to model the effect of the inner structure of the packed bed can lead to more efficient design for improved filtration. Because of its significant efficiency, gravitational sphere packing is employed in this work to simulate a bed of mono‐sized randomly packed spheres. The simulated bed provides a way to visualize the pore network and estimate the pore size distribution associated with the void space between particles. Furthermore, by subsequently introducing fine particles into the bed, we evaluate the mass‐rate of fine particles passing through and possibly saturating the packed bed. Results show that fine particles between 15% and 25% of the coarse particle size can be physically strained within the randomly packed bed. These results differ significantly from the results obtained assuming a periodically spaced bed. The technique therefore provides an efficient yet accurate alternative for understanding how fine particles pass through a coarse particulate medium. Copyright © 2014 John Wiley & Sons, Ltd.     

走滑环境中陆相盆地充填层序特征以云南先锋盆地为例

走滑陆相盆地在构造背景、沉积范围、物源供应方式、基底地形特征、沉积物特点和沉积层序特征方面均与大陆边缘海相盆地有重要区别，通常具有多物源、多沉积体系、相带窄、相变快的特点。本文在对云南先锋盆地充填序列、沉积格架、构造格架和区域构造背景分析基础上，讨论了在走滑构造背景下陆相盆地的一系列特征以及走滑型陆相盆地的层序地层特点.     

Hydrodynamic Modeling of the Roman Harbor of Portus in the Tiber Delta: The Impact of the North‐Eastern Channel on Current and Sediment Dynamics

This study examines the hydrodynamics of the Roman harbor of Portus during the Trajanic period (second century A.D.). We have evaluated the impact of the north‐eastern channel on water circulation and sediment resuspension within the harbor in relation to the problem of sediment infilling. We used a 2D numerical model to compute the distribution of the depth‐averaged current velocity, bed shear stress, and kinetic energy induced by each of the four prevailing local winds. First, the results confirm the persistent conditions for sediment infilling over the entire western harbor basin, even when the north‐eastern channel is operational; these conditions being present for 61.3% (40.3%) of the summer period and up to 68% (44.2%) of the winter period in the south‐western harbor basin. Second, the results show that favorable conditions for navigation occur in the central basin, leading to the landing installations, under the S 180°, SE 135°, and NE 22° winds when the north‐eastern channel is operational; these conditions being present for 21% of the summer period and 51.3% of the winter period. Inversely, the access to the landing installations with both channels operating remains affected by sediment infilling, for 40.3% of the summer period and 16.7% of the winter period.     

Hydro-mechanical coupled analysis of near-wellbore fines migration from unconsolidated reservoirs

Oil or gas production from unconsolidated reservoirs could be hampered by sand migration near the wellbore. This paper presents a numerical investigation of production-induced migration of fine sands towards a wellbore drilled in a gap-graded sediment. The solid–fluid interaction is simulated by coupling the discrete element method and the dynamic fluid mesh. With the merit of DEM and a dynamic mesh, the model is capable of naturally capturing particle movements and spatiotemporal variations of hydraulic properties of the sediment at the pore scale. The results show that fine particles are mobilized by radial flow under an imposed hydraulic gradient, and the increase in the hydraulic gradient causes an increase in the fines production. The microscopic pattern of sand migration is clearly visualized through the simulation. The presence of fine particles affects the process of fines migration through two competing mechanisms. Under a low fine content, fine sands mainly serve as the fines production source, and thus, fines production is enhanced as the fine content increases up to a critical value, beyond which fines production is weakened with a further increase in the fine content since the blocking effect gradually dominates. A barrier layer is likely formed during sand migration due to settling and jamming of fine sands at the throats of pores, as fine sands migrate with the radial flow towards the wellbore. This layer is helpful to slow down sand migration, while it could impede production due to reduced permeability in the affected reservoir.

     

Two methods for pore network of porous media

Two approaches of generating pore networks of porous media are presented to capture the pore fabric. The first methodology extracted pore structure from a computer simulated packing of spheres. The modified Delaunay tessellation was used to describe the porous media, and modified Nelder–Mead method in conjunction with three pore‐merging algorithms was used to generate the pore size and coordination number distributions of the randomly packed spheres. The Biconical Abscissa Asymmetric CONcentric bond was used to describe the connection between two adjacent voids. This algorithm was validated by predicting pore structure of a cubic array of spheres of equal radius with known pore sizes, throat sizes and coordination number distributions. The predicted distributions of pore structure agreed well with the measured. Then, the algorithm was used to predict pore structure and permeability of randomly packed spherical particles, and predicted permeability values were compared with published experimental data. The results showed that the predicted permeability values were in good agreement with those measured, confirming the proposed algorithm can capture the main flow paths of packed beds. The second methodology generated an equivalent pore network of porous media, of which the centers of voids were located in a regular lattice with constant pore center distance. However, this network allowed for matching both main geometrical and topological characteristics of the porous media. A comparison of the two approaches suggested that the second approach can also be used as a predictive tool to quantitatively study the microscopic properties of flow through porous media. Copyright © 2012 John Wiley & Sons, Ltd.     

漏斗式全沙排沙技术及其应用

主要介绍了一种新型泥沙处理技术,即“漏斗式全沙排沙术”.该项技术既可以排除粗颗粒的推移质泥沙,也可以排除细颗粒的悬移质泥沙.对0.5mm以上的粗颗粒泥沙的截沙率为100%,对0.5mm以下细颗粒泥沙的截沙率为90%,排沙耗水量仅占渠道引水量的2%～5%.     

Effects of bed longitudinal inflexion and sediment porosity on basal entrainment mechanism: insights from laboratory debris flows

Channel morphology and bed sediment erodibility are two crucial factors that significantly affect debris flow entrainment processes. Current debris flow entrainment models mostly hypothesize the erodible beds are infinite with uniform slopes. In this study, a series of small-scale flume experiments were conducted to investigate the effects of bed longitudinal inflexion and sediment porosity on basal entrainment characteristics. Experimental observations revealed that sediment entrainment is negligible at early stages and accelerates rapidly as several erosion points appear. Continual evolution of flow-bed interfaces changes interactions between debris flows and bed sediments, rendering the interfacial shear action involved into a mixed shear and frontal collisional action. Lower bed sediment porosity will change the spatial arrangement and orientation of particle mixture, strengthen the interlocking and anti-slide forces of adjacent sediment particles, and promote the formation of particle clusters, all of which will increase bed sediment resistance to erosion. By examining the post-experimental bed morphology, the slope-cutting amounts and topographic reliefs are determined to positively correlate with longitudinal transition angles. These high topographic reliefs may indicate the propensity of triangular slab erosion, rather than strip-shaped slab erosion, in non-uniform channels with relatively steep erodible beds. Empirical formulas are obtained that denote the relationships among bed sediment strength, channel curvature radius, and sediment porosity through a multi-parameter regression analysis. This study may aid in clarifying the complex coupling effects of spatial variations in debris flow dynamics as well as sediment erodibility and bed morphology in non-uniform channels with abundant seismic loose material.

     

Sensitivity of incipient particle motion to fluid flow penetration depth within a packed bed

A discrete element method is applied to a three‐dimensional analysis related to sediment entrainment on a micro‐scale. Sediment entrainment is the process by which a fluid medium accelerates particles from rest and advects them upward until they are either transported as bedload or suspended by the flow. Modelling of the entrainment process is a critically important aspect for studies of erosion, pollutant resuspension and transport, and formation of bedforms in environmental flows. Previous discrete element method studies of sediment entrainment have assumed the flow within the particle bed to be negligible and have only allowed for the motion of the topmost particles. At the same time, micro‐scale experimental studies indicate that there is a small slip of the fluid flow at the top of the bed, indicating the presence of non‐vanishing fluid velocity within the topmost bed layers. The current study demonstrates that the onset of particle incipient motion, which immediately precedes particle entrainment, is highly sensitive to this small fluid flow within the topmost bed layers. Using an exponential decay profile for the inner‐bed fluid flow, the discrete element method calculations are repeated with different fluid penetration depths within the bed for several small particle Reynolds numbers. For cases with slip velocity corresponding to that observed in previous experiments with natural sediment, the predicted particle velocity is found to be a few percent of the fluid velocity at the top of the viscous wall layer, which is a reasonable range of velocities for observation of incipient particle motion. This method for prescribing the fluid flow within the particle bed allows for the current discrete element method to be extended in future studies to the analysis of sediment entrainment under the influence of events such as turbulent bursting. Additionally, predictions for the slip velocities and fluid flow profile within the bed suggest the need for further experimental studies to provide the data necessary for additional improvement of the discrete element method models.     

A nearshore model to investigate the effects of seagrass bed geometry on wave attenuation and suspended sediment transport

The effects of seagrass bed geometry on wave attenuation and suspended sediment transport were investigated using a modified Nearshore Community Model (NearCoM). The model was enhanced to account for cohesive sediment erosion and deposition, sediment transport, combined wave and current shear stresses, and seagrass effects on drag. Expressions for seagrass drag as a function of seagrass shoot density and canopy height were derived from published flume studies of model vegetation. The predicted reduction of volume flux for steady flow through a bed agreed reasonably well with a separate flume study. Predicted wave attenuation qualitatively captured seasonal patterns observed in the field: wave attenuation peaked during the flowering season and decreased as shoot density and canopy height decreased. Model scenarios with idealized bathymetries demonstrated that, when wave orbital velocities and the seagrass canopy interact, increasing seagrass bed width in the direction of wave propagation results in higher wave attenuation, and increasing incoming wave height results in higher relative wave attenuation. The model also predicted lower skin friction, reduced erosion rates, and higher bottom sediment accumulation within and behind the bed. Reduced erosion rates within seagrass beds have been reported, but reductions in stress behind the bed require further studies for verification. Model results suggest that the mechanism of sediment trapping by seagrass beds is more complex than reduced erosion rates alone; it also requires suspended sediment sources outside of the bed and horizontal transport into the bed.     

Slurry-flow deposits in the Britannia Formation (Lower Cretaceous), North Sea: a new perspective on the turbidity current and debris flow problem

The Lower Cretaceous Britannia Formation (North Sea) includes an assemblage of sandstone beds interpreted here to be the deposits of turbidity currents, debris flows and a spectrum of intermediate flow types termed slurry flows. The term ‘slurry flow’ is used here to refer to watery flows transitional between turbidity currents, in which particles are supported primarily by flow turbulence, and debris flows, in which particles are supported by flow strength. Thick, clean, dish‐structured sandstones and associated thin‐bedded sandstones showing Bouma T_b–e divisions were deposited by high‐ and low‐density turbidity currents respectively. Debris flow deposits are marked by deformed, intraformational mudstone and sandstone masses suspended within a sand‐rich mudstone matrix. Most Britannia slurry‐flow deposits contain 10–35% detrital mud matrix and are grain supported. Individual beds vary in thickness from a few centimetres to over 30 m. Seven sedimentary structure division types are recognized in slurry‐flow beds: (M₁) current structured and massive divisions; (M₂) banded units; (M₃) wispy laminated sandstone; (M₄) dish‐structured divisions; (M₅) fine‐grained, microbanded to flat‐laminated units; (M₆) foundered and mixed layers that were originally laminated to microbanded; and (M₇) vertically water‐escape structured divisions. Water‐escape structures are abundant in slurry‐flow deposits, including a variety of vertical to subvertical pipe‐ and sheet‐like fluid‐escape conduits, dish structures and load structures. Structuring of Britannia slurry‐flow beds suggests that most flows began deposition as turbidity currents: fully turbulent flows characterized by turbulent grain suspension and, commonly, bed‐load transport and deposition (M₁). Mud was apparently transported largely as hydrodynamically silt‐ to sand‐sized grains. As the flows waned, both mud and mineral grains settled, increasing near‐bed grain concentration and flow density. Low‐density mud grains settling into the denser near‐bed layers were trapped because of their reduced settling velocities, whereas denser quartz and feldspar continued settling to the bed. The result of this kinetic sieving was an increasing mud content and particle concentration in the near‐bed layers. Disaggregation of mud grains in the near‐bed zone as a result of intense shear and abrasion against rigid mineral grains caused a rapid increase in effective clay surface area and, hence, near‐bed cohesion, shear resistance and viscosity. Eventually, turbulence was suppressed in a layer immediately adjacent to the bed, which was transformed into a cohesion‐dominated viscous sublayer. The banding and lamination in M₂ are thought to reflect the formation, evolution and deposition of such cohesion‐dominated sublayers. More rapid fallout from suspension in less muddy flows resulted in the development of thin, short‐lived viscous sublayers to form wispy laminated divisions (M₃) and, in the least muddy flows with the highest suspended‐load fallout rates, direct suspension sedimentation formed dish‐structured M₄ divisions. Markov chain analysis indicates that these divisions are stacked to form a range of bed types: (I) dish‐structured beds; (II) dish‐structured and wispy laminated beds; (III) banded, wispy laminated and/or dish‐structured beds; (IV) predominantly banded beds; and (V) thickly banded and mixed slurried beds. These different bed types form mainly in response to the varying mud contents of the depositing flows and the influence of mud on suspended‐load fallout rates. The Britannia sandstones provide a remarkable and perhaps unique window on the mechanics of sediment‐gravity flows transitional between turbidity currents and debris flows and the textures and structuring of their deposits.     

Shaking table tests to investigate the influence of various factors on the liquefaction resistance of sands

This paper presents the shaking table studies to investigate the factors that influence the liquefaction resistance of sand. A uniaxial shaking table with a perspex model container was used for the model tests, and saturated sand beds were prepared using wet pluviation method. The models were subjected to horizontal base shaking, and the variation of pore water pressure was measured. Three series of tests varying the acceleration and frequency of base shaking and density of the soil were carried out on sand beds simulating free field condition. Liquefaction was visualized in some model tests, which was also established through pore water pressure ratios. Effective stress was calculated at the point of pore water pressure measurement, and the number of cycles required to liquefy the sand bed were estimated and matched with visual observations. It was observed that there was a gradual variation in pore water pressure with change in base acceleration at a given frequency of shaking. The variation in pore water pressure is not significant for the range of frequency used in the tests. The frequency of base shaking at which the sand starts to liquefy when the sand bed is subjected to any specific base acceleration depends on the density of sand, and it was observed that the sand does not liquefy at any other frequency less than this. A substantial improvement in liquefaction resistance of the sand was observed with the increase in soil density, inferring that soil densification is a simple technique that can be applied to increase the liquefaction resistance.     

Combined sewer overflows, sediment accumulation and element patterns of river bed sediments: a quantitative study based on mixing models of composite fingerprints

Stormwater runoff from urban areas and subsequent stormwater-induced effluents from sewer systems may deteriorate the quality of the receiving water sediments by emitting particulate matter and associated pollutants. However, the relevance of stormwater and combined sewer effluents for the pollution of bed sediments was not yet quantified. Therefore, we applied a multivariate mixing model of composite fingerprints to investigate how much stormwater effluent may contribute to the accumulation of fines and associated pollutants in the bed sediment of the Bode River, Germany. In our study, stormwater and combined sewer effluents contribute about 10 % of the fines accumulated in the bed sediment. As stormwater overflow fines are a major carrier of C, N, P, Cu and Zn, up to 40 % of these pollutants in the bed sediment originate from stormwater effluents. Especially N and Zn have to be seen critically because high nutrient concentrations trigger excessive macrophyte growth within the studied river stretch and Zn contents exceed German sediment quality standards (LAWA 1998) in the bed sediment.     

Cave stratigraphy in western Norway; multiple Weichselian glaciations and interstadial vertebrate fauna

Skjonghclleren is a marine-cut cave with 15–20. m thick pre-Holocene sediments. Corings and excavations reveal three beds of extremely fine-grained, laminated sediments alternating with blocky sediments. The laminated beds are interpreted as glaciolacustrine sediments deposited subglacially at times when ice sheets covered the area, suggesting at least three glaciations after the cave was formed. The blocky/diamictic sediments were formed by frost-shattered blocks from the roof of the cave during ice-free periods, and mixing with the fines through slow mass movements along the floor of the cave. In the diamictic sediment beneath the uppermost laminated bed, almost 7,000 bone and teeth fragments of birds, mammals and fish were found. Birds dominated, with little auk and brunnich's guillemot as the most frequently occurring species. Arctic fox was the dominating mammal. During climatic optimum of the interstadial, conditions seem to have been similar to present-day coastal Finnmark, with North Atlantic warm water entering the Norwegian Sea. Two radiocarbon dates on bones and three Uranium series dates on speleothems from this bed all cluster around 30,000 B.P., i.e., the end of the Ålesund interstadial. Above the uppermost laminated bed, bone fragments of birds, fish and mammals, deposited between c . 12,000 and c . 10,000 B.P., were found. Little auk dominate. The occurrence of squirrel is worth noting since it is limited mainly to areas with coniferous forest today. The beds below the 30,000B.P. bed are poorly dated or undated, but it is tentatively concluded that the entire sediment sequence was deposited during the Weichselian stage. It seems that the cave was formed at a high relative sea-level stand sometime during the Early Weichsclian. Two recorded palaeomagnetic excursions seem to correlate with the Laschamp/Olby and the Lake Mungo events, respectively.     

New insight into the evolution of large-volume turbidity currents: comparison of turbidite shape and previous modelling results

The Marnoso Arenacea Formation provides the most extensive correlation of individual flow deposits (beds) yet documented in an ancient turbidite system. These correlations provide unusually detailed constraints on bed shape, which is used to deduce flow evolution and assess the validity of numerical and laboratory models. Bed volumes have an approximately log‐normal frequency distribution; a small number of flows dominated sediment supply to this non‐channelized basin plain. Turbidite sandstone within small‐volume (<0·7 km³) beds thins downflow in an approximately exponential fashion. This shape is a property of spatially depletive flows, and has been reproduced by previous mathematical models and laboratory experiments. Sandstone intervals in larger‐volume (0·7–7 km³) beds have a broad thickness maximum in their proximal part. Grain‐size trends within this broad thickness maximum indicate spatially near‐uniform flow for distances of ∼30 km, although the flow was temporally unsteady. Previous mathematical models and laboratory experiments have not reproduced this type of deposit shape. This may be because models fail to simulate the way in which near bed sediment concentration tends towards a constant value (saturates) in powerful flows. Alternatively, the discrepancy may be the result of relatively high ratios of flow thickness and sediment settling velocity in submarine flows, together with very gradual changes in sea‐floor gradient. Intra‐bed erosion, temporally varying discharge, and reworking of suspension fallout as bedload could also help to explain the discrepancy in deposit shape. Most large‐volume beds contain an internal erosion surface underlain by inversely graded sandstone, recording waxing and waning flow. It has been inferred previously that these characteristics are diagnostic of turbidites generated by hyperpycnal flood discharge. These turbidites are too voluminous to have been formed by hyperpycnal flows, unless such flows are capable of eroding cubic kilometres of sea‐floor sediment. It is more likely that these flows originated from submarine slope failure. Two beds comprise multiple sandstone intervals separated only by turbidite mudstone. These features suggest that the submarine slope failures occurred as either a waxing and waning event, or in a number of stages.