Click to view article in PDF format.
Sedimentation from Jets: A Depositional Model for Clastic Deposits of all Scales and Environments*
By
D.C.J.D. Hoyal, J.C. Van Wagoner, N.L. Adair, M. Deffenbaugh, D. Li, T. Sun, C. Huh, and D.E. Giffin
Search and Discovery Article #40082 (2003)
ExxonMobil Upstream Research Company, Houston, Texas, P.O. Box 2189, Houston, Texas, 77252-2189, Tel: 713-431-4416, Fax: 713-431-6336 ([email protected])
*Adapted from "extended abstract" of poster session presentation at AAPG Annual Meeting, May 14, 2003, Salt Lake City, Utah.
NOTE: This is the second of six presentations on the general subject of the shapes of siliciclastic sedimentary bodies presented by this group of ExxonMobil researchers under the leadership of John C. Van Wagoner. Click to view a list of all these articles.
Sedimentary
geologists typically employ depositional models in order to predict subsurface
sand
body
morphology and distribution. One assumption of these models is that
sedimentary processes are environment and scale specific (Walker and James,
1992). However, because the physics of sediment and fluid transport operate
independently of environment and scale (e.g., Fischer et al., 1979; Middleton
and Southard, 1978), a more logical approach may be to develop depositional
models that are based on the properties of the decelerating flows responsible
for most clastic deposits (Bates, 1953).
Applying a
flow-based method significantly reduces the number of unique depositional
models, potentially simplifying and improving sand
body
prediction. It also
provides the context for a transition from the current qualitative depositional
models to more quantitative models that are based on well-established laws of
fluid and sediment transport dynamics. To be useful, these models must predict
the properties of 3-D sedimentary bodies over a wide range of complexity, from
simple deposits of a single flow to highly self-organized, avulsive complexes.
Real progress in this field can only be made with a coordinated program of
sediment transport and fluid mechanics research based on laboratory experiments,
mathematical models and fieldwork. Sedimentary process models used in the
forward sense can provide quantitative estimates of grain-size, layering, and
connectivity at various scales over complex 3D-topography and potentially
include syndepositional tectonics and changing sea level. These models can also
be used in the inverse sense to investigate the strong interdependence of
deposit properties, for example, grain-size distribution, relative thickness
(Sadler, 1982), facies (Bouma, 1962), and erosion patterns inherited from the
flow field and associated grain transport mechanism. Results from these
investigations can be used to constrain predictions of micro-scale properties
like grain-size distribution from some observable macro-scale property like
facies or relative thickness.
|
uFeatures of jet-plume pair deposits
uFeatures of jet-plume pair deposits
uFeatures of jet-plume pair deposits
uFeatures of jet-plume pair deposits
uFeatures of jet-plume pair deposits
uFeatures of jet-plume pair deposits
uFeatures of jet-plume pair deposits
uFeatures of jet-plume pair deposits
uFeatures of jet-plume pair deposits
|
We propose that in all environments and at all scales there is a dominant flow type that creates most clastic deposits. Consequently these deposits, including deltas, submarine fans, crevasse splays, overwash fans, tidal bars, lower shoreface deposits and bedforms, are constructed from fundamentally similar elements. They are point-sourced, down-flow expanding deposits with well-defined horizontal and vertical distributions of thickness, grain-size, facies, erosion and bedding type. Even deposits in channels, for example braid and point bars in rivers and submarine channels, exhibit features similar both to the above listed types and to one another. Allen (1982) proposed that the dominant mechanism controlling the geometry of most clastic deposits is deceleration associated with flow separation and the growth of turbulent shear layers. We suggest that the observed commonality in the shape and property distribution results from a subset of these separated turbulent flows, those that expand and decelerate from a submerged, focused slot or orifice into a lower velocity region. These flows are commonly called turbulent jets. While traditionally applied to coarse-grained parts of deltas (Bonham-Carter and Sutherland, 1968; Syvitski et al., 1998), a jet model may also explain the delta-like properties of other clastic deposits, for example submarine fans (Beaubouef et al., I, II, 2003) and river bars (Dunn et al., 2003). Other point-source expanding flows including: pure plumes or modified jets (e.g., bi-directional oscillating jets, jets bent by a cross-current), may produce deposits with some similarities to unidirectional jet-plume pair deposits, because these flows also expand and decelerate from a submerged orifice. However, these deposits may be more heterogeneous and geometrically complex. Fluid entrainment in turbulent jets is the most efficient mechanism for flow deceleration. Deceleration immediately downstream of a topographic discontinuity or submerged orifice requires the transfer of momentum from the moving fluid to the slower or stationary ambient fluid. In dilute Newtonian flows (<10% v/v), the major momentum transfer mechanism is turbulence (Fischer et al., 1979). If the flow boundary downstream from the orifice expands more than a few degrees the main flow and the surrounding fluid will separate causing a distinct velocity discontinuity. Strong turbulence will develop at this discontinuity due to a Kelvin-Helmholtz instability creating a shear layer, which will grow downstream as the momentum diffuses through a cascade of turbulent eddies. Viscous and gravitational forces inhibit turbulence development at the shear layer. If the viscous or gravitational forces significantly exceed the inertial force, turbulence will not develop (Allen, 1982; Turner, 1973). Consequently the rate of deceleration from an orifice is controlled by two variables: (1) Reynolds number: the ratio of inertial force to viscous force (Re = DrDUL/m); and (2) Froude number: the ratio of the inertial force to the gravity force (Fr = DU/[Drgh]0.5). Here L and h are characteristic length and height scales, g is gravitational acceleration and m is the dynamic viscosity. The velocity (DU) and density (Dr) gradients are calculated between the region of separated flow and the surrounding fluid environment and not necessarily integrated across the entire flow. At the orifice most natural sedimenting flows have a high Reynolds number and can be considered fully turbulent (Re > 5000). However, the Froude number may vary widely, producing either an inertial dominated supercritical jet (Fr > 1) or gravity dominated subcritical plume (Fr < 1) (Fischer et al., 1979). A pure, homopycnal (Dr = 0), completely inertial jet, similar to the high Fr jet in Figure 1a, represents an asymptotic end-member with the smallest possible expansion angle. As the Froude number decreases the jet becomes wider as gravitational collapse becomes more important. Flow separation in natural flows is commonly focused at a 3-D point-source orifice or a 2-D line-source orifice. For a channel exiting from a canyon, expansion is in 3D from a point-source. In flow over 2-D bedforms, expansion is in 2-D from a half orifice (step) producing a half jet (Jopling, 1965). In the classic literature this type of separation has been called both a jet (Jopling, 1965) and a wake (Raudkivi, 1963). Wakes and jets are geometrically equivalent as they are both products of relative fluid motion and differ only in the coordinate systems that define them (Allen, 1982). Therefore, the jet is not a single flow but a distinctive family of flows belonging to the larger class of separated flows that includes wakes (Allen, 1982). A wake behind a flow-blocking object creates two shear layers, which can also be considered as half jets expanding into the separated/blocked flow region. Flow between two objects creates two shear layers; this is often called a jet but can also be considered as two half wakes, one behind each object. In essence all jets need a region of slower moving fluid, or separated flow blocked behind some obstruction.
Examination of the
Froude number equation shows that jet flow is favored when one or more of the
factors dominates: (1) velocity of the flow is significantly higher that the
surrounding flow (high
DU = large inertial force) (2) the density contrast
between the flow and its surroundings is weak (small
Dr = small gravity force),
or (3) if the flow is very thin (small h= small gravity force). Consequently,
all turbulent sediment-laden flows are either jets themselves or, if they are
plumes, they will contain separated flow regions that are jets (Allen, 1982).
Because stratigraphically significant sedimentation is associated with As a point-source jet expands from the orifice, velocity decreases downstream, and the density (gravity force) of the flow begins to dominate the inertia (motion force). Since all sediment-laden flows have a different density than the surrounding fluid, they will transition at a submerged hydraulic jump, into a gravity-dominated plume that expands through gravitational collapse. We call this fundamental decelerating depositing flow the jet-plume pair. In an expanding jet-plume pair flowing over a horizontal or downward sloping surface the submerged hydraulic jump is usually weak and diffuse; so deposition is predominantly in the jet. Submerged hydraulic jump deposits may become important if the flow is one-dimensional or if the topography slopes upstream. The coarsest sediment is deposited in the jet region, producing better hydrocarbon reservoir deposits than the plume deposit. An experimental wall-bounded jet-plume pair is presented in Figure 1a. Jet-plume pair properties include: (1) characteristic spreading angle controlled by vertical confinement, bottom friction and Froude number, (2) strongly turbulent in the jet region, (3) maximum turbulence intensity occurring some distance downstream from orifice [8 orifice diameters for a free jet (Fischer et al., 1979)], and (4) power law downstream and gaussian across-stream velocity distributions (Fischer et al., 1979).
Although boundary
and stratification conditions vary between jets in different depositional
environments, the functional form of the velocity decay and dependence on
controlling variables is the same due to the universality of the jet
deceleration mechanism. Consequently, the only Features of Jet-Plume Pair Deposits All of the features of the deposit of a point-source jet-plume pair can be explained by convolving the sediment flux process at the bed (erosion g bypass g suspension deposition) with the spatially decaying velocity field. The jet flow field is more systematic and predictable than a plume because jet behavior is largely independent of slope and hence topography (Imran et al., 2002), while a plume is very slope sensitive. Consequently jet deposits have very characteristic bed patterns and vertical successions. These characteristics, illustrated in Figure 1b, include: (1) point-source planform expanding outwards at characteristic expansion angle controlled by the jet velocity field and grain-size, (3) near orifice erosion / incipient channel formation surrounded by levees (Imran et al., 2002), (4) progression of facies with decreasing velocity: erosion region g bedform region g pure suspension deposition region, and (5) thickness decreasing exponential-linear downstream and gaussian-like across-stream. In a larger scale deposit with paleo-flow indicators, vectors will be spreading from point-source orifice.
Strong turbulence
in the jet region has great potential for eroding the substrate creating a
flute-like erosional scour. A jet experiment that illustrates a
Grain-size decay
in a jet-plume pair deposit is also highly systematic, reflecting systematic
velocity and concentration decay. The results of an experimental high Froude
number flow (Fr = 14.6) with over 300 grain-size samples collected on a grid is
presented in Figure 3. Although downstream velocity decay in the jet region
follows a power law, the grain-size decay is exponential-linear in the
downstream direction due to convolution of the jet velocity field and
depositional flux process (Figure 3c). In the cross-stream direction,
convolution of a gaussian transverse velocity decay and deposition flux gives a
peaked gaussian-like distribution (Figure 3c). Because the thickness and
grain-size distributions are both ultimately controlled by the flow field, there
is remarkable similarity in the shape of thickness and grain-size contours
(Figure 3b), creating a strong correlation between statistics of the grain-size
distribution and thickness as observed in Figure 3d. The exact form or the
correlation depends on the position in the deposit, suspension, bypass or bedload region. As discussed earlier, for a particular boundary configuration,
the planform of the jet-plume pair flow and hence that of the deposit are a
function of Froude number. In Figure 4 we present experimental deposit shapes
and axial grain-size decay for multiple experiments on a flat horizontal plate
with different Froude numbers. The orientations of bedform crests illustrated in
Figure 4a generally mimic the deposit thickness contours which vary in shape
from elongate for a jet to almost circular for the low Froude number
(plume-like) case (Figure 4b). A cross-plot of plan shape against Froude number
is presented in Figure 4c where the plan shape is calculated as the ratio of
maximum length to maximum width of a specific thickness contour (0.5 relative
thickness). This plot (Figure 4c) shows that at high Froude numbers (>5) the deposit
approaches an elongate shape about 4 times longer than wide, a feature which
seems reasonable, because as discussed earlier, the pure jet flow field is the
narrowest and asymptotic in all properties. This relationship between
The exponential
coefficient of the down-axis grain-size decay (P[50]) is also a function of
Froude number but varies inversely to grain-size (Figure 4c). Spatial grain-size
decay is faster for low orifice Froude number flows and asymptotically reaches a
slower decay at high Froude numbers. This high Froude number limit may represent
the uniform suspension deposition rate of Martin and Nokes (1988), and faster
deposition rates at lower Froude numbers may result from collapse of the
concentration profile or flow thinning. Because both
So far the
discussion has been restricted implicitly to steady flows. However, all
sedimentary flows are fundamentally unsteady due to interaction of the flow with
the accreting deposits and variation in water and sediment discharge at the
orifice. Real deposits can be considered as a composite or sandwich of steady
deposits or layers, laid down sequentially in small time steps. Each layer of
the sandwich has a certain properties that are spatially correlated through the
flow velocity and concentration fields at every x, y position. These properties
include relative thickness, grain-size, and facies. Consequently a complex
We suggest that
environment of deposition and scale may not be first-order controls on the
properties of clastic sand bodies because the physics of turbulent flow
deceleration and sediment transport transcend many clastic depositional
environments and scales. There are a limited number of ways to decelerate a
flow; of these spatial expansion is the most common. Expansion in turbulent
flows creates flow separation, shear layers, and turbulent jets. Consequently,
turbulent jet deposits may be the most fundamental and ubiquitous identifiable
bodies in the clastic sedimentary record. Apart from expansion angle, jet-plume
pair deposits exhibit a large degree of similarity in all environments where
jets dominate over ambient flow and where re-working by other processes has not
obscured them. The properties of these jet-plume pair deposits are strongly
correlated to their shape through the turbulent flow field and depositional
process. Physics-based jet-plume pair depositional models can be used to
understand the link between 3D-sediment Albertson, M. L., Y. B. Dai, R. A. Jensen, and H. Rouse, 1950, Diffusion of submerged jets: ASCE Transactions, p. 639-697. Allen, J. R. L., 1982, Chapter 4. Sedimentation from jets and separated flows, Sedimentary Structures Their Character and Physical Basis, v. II: Amsterdam, Oxford, New York, Elsevier Scientific Publishing Company, p. 133-171. Bates, C. C., 1953, Rational theory of delta formation: Bulletin of the American Association of Petroleum Geologists, v. 37, p. 2119-2162. Beaubouef, R.T., V. Abreu, and N.L. Adair, 2003, Ultra-high resolution 3-D characterization of deep-water deposits- I: A new approach to understanding the stratigraphic evolution of intra-slope depositional systems, Search and Discovery Article #40083. Beaubouef, R.T., Van Wagoner, J.C., and N.L. Adair, 2003, Ultra-high resolution 3-D characterization of deep-water deposits- II: Insights into the evolution of a submarine fan and comparisons with river deltas: Search and Discovery Article #40084. Bouma, A. H., 1962, Sedimentology of Some Flysch Deposits: Amsterdam New York, Elsevier Publishing Company, 168 p. Bonham-Carter, G., and A. J. Sutherland, 1968, Mathematical model and Fortran IV program for computer simulation of deltaic sedimentation: Kansas Geological Survey Computer Contribution, v. 24, p. 56.
Dunn, Paul A., J.C. Van Wagoner, and
M. Deffenbaugh, 2003, Hierarchical, Self-Affine Fluvial Sand Fischer, H. B., J. E. List, R. C. Y. Koh, J. Imberger, and N. H. Brooks, 1979, Mixing in Inland and Coastal Waters: San Diego, Academic Press. Imran, J., G. Parker, and P. Harff, 2002, Experiments on incipient channelization of submarine fans: Journal Hydraulic Research, v. 40, p. 21-32. Jopling, A. V., 1965, Hydraulic factors controlling the shape of laminae in laboratory deltas: Journal of Sedimentary Petrology, v. 35, p. 777-791. Martin, D., and R. Nokes, 1988, Crystal settling in a vigorously convecting magma chamber: Nature, v. 332, p. 534-536. Middleton, G. V., and J. B. Southard, 1978, Mechanics of Sediment Movement: SEPM Short Course Number 3, March 29-30, 1977: Binghampton. Pilkey, O. H., S. D. Locker, and W. J. Cleary, 1980, Comparison of sand-layer geometry on flat floors of 10 modern depositional basins: AAPG Bulletin, v. 64, p. 841-856. Raudkivi, A. J., 1963, Study of sediment ripple formation: Journal of the Hydraulics Division, ASCE, v. 89, p. 15-33. Sadler, P. M., 1982, Bed-thickness and grain size of turbidites: Sedimentology, v. 29, p. 37-51. Sidorchuk, A., ed., 1996, The Structure of River Bed Relief: Coherent Flow Structures, in Open Channels, Ashworth, P.J.., Bennett, S.J., Best, J.L., McLelland, S.J., eds, Wiley, NY, p. 397-421. Syvitski, J. P. M., K. I. Skene, M. K. Nicholson, and M. D. Morehead, 1998, Plume 1.1: deposition of sediment from a fluvial plume: Computers & Geosciences, v. 24, p. 159-171. Turner, J. S., 1973, Bouyancy Effects in Fluids: Cambridge Monographs on Mechanics and Applied Mathematics, Cambridge University Press, 367p. Van Wagoner, J.C., D.C.J.D. Hoyal, N.L. Adair, T. Sun, R.T. Beaubouef, M. Deffenbaugh, P.A. Dunn, C. Huh, and D. Li, 2003, Energy Dissipation and the Fundamental Shape of Siliciclastic Sedimentary Bodies: Search and Discovery Article #40080. Walker, R. G., and N. P. James, 1992, Facies Models: Response to Sea Level Change, Geological Association of Canada, 454 p. Wright, L. D., 1977, Sediment transport and deposition at river mouths: a synthesis: Geological Society of America Bulletin, v. 88, p. 857-868. |
