3D image-based characterization of fluid displacement in a Berea core

Improved network flow models require the incorporation of increasingly accurate geometrical characterization of the microscale pore structure as well as greater information on fluid–fluid interaction (interfaces) at pore scales. We report on three dimensional (3D) pore scale medium characterization, absolute permeability computations for throat structures, and pore scale residual fluid distribution in a Berea core. X-ray computed microtomography combined with X-ray attenuating dopants is used to obtain 3D images of the pore network and to resolve phase distributions in the pore space.     

Two-dimensional power series solution for non-axisymmetrical transport in a radially convergent tracer test with scale-dependent dispersion

It has been known for many years that dispersivity increases with solute travel distance in a subsurface environment. The increase of dispersivity with solute travel distance results from the significant variation of hydraulic properties of heterogeneous media and was identified in the literature as scale-dependent dispersion. This study presents an analytical solution for describing two-dimensional non-axisymmetrical solute transport in a radially convergent flow tracer test with scale-dependent dispersion. The power series technique coupling with the Laplace and finite Fourier cosine transform has been applied to yield the analytical solution to the two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates with variable-dependent coefficients. Comparison between the breakthrough curves of the power series solution and the numerical solutions shows excellent agreement at different observation points and for various ranges of scale-related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at any observation point.     

Laboratory simulation of the salt weathering of schist: II. Fragmentation of fine schist particles

Recent developments in long term landform evolution modelling have created a new demand for quantitative salt weathering data, and in particular data describing the size distribution of the weathered rock fragments. To enable future development of rock breakdown models for use in landscape evolution and soil production models, laboratory work was undertaken to extend existing schist/salt weathering fragmentation studies to include an examination of the breakdown of sub‐millimetre quartz chlorite schist particles in a seasonally wet tropical climate. Laser particle sizing was used to assess the impact of different experimental procedures on the resulting particle size distribution. The results reveal that salt weathering under a range of realistic simulated tropical wet season conditions produces a significant degree of schist particle breakdown. The fragmentation of the schist is characterized by splitting of the larger fragments into mid‐sized product with finer material produced, possibly from the breakdown of mid‐sized fragments when weathering is more advanced. Salinity, the salt addition method and temperature were all found to affect weathering rates. Subtle differences in mineralogy also produce variations in weathering patterns and rates. It is also shown that an increase in drying temperature leads to accelerated weathering rates, however, the geometry of the fracture process is not significantly affected. Copyright © 2006 John Wiley & Sons, Ltd.     

Nonlinear determinism in river flow: prediction as a possible indicator

Whether or not river flow exhibits nonlinear determinism remains an unresolved question. While studies on the use of nonlinear deterministic methods for modeling and prediction of river flow series are on the rise and the outcomes are encouraging, suspicions and criticisms of such studies continue to exist as well. An important reason for this situation is that the correlation dimension method, used as a nonlinear determinism identification tool in most of those studies, may possess certain limitations when applied to real river flow series, which are always finite and often short and also contaminated with noise (e.g. measurement error). In view of this, the present study addresses the issue of nonlinear determinism in river flow series using prediction as a possible indicator. This is done by (1) reviewing studies that have employed nonlinear deterministic methods (coupling phase‐space reconstruction and local approximation techniques) for river flow predictions and (2) identifying nonlinear determinism (or linear stochasticity) based on the level of prediction accuracy in general, and on the prediction accuracy against the phase‐space reconstruction parameters in particular (termed as the ‘inverse approach’). The results not only provide possible indications to the presence of nonlinear determinism in the river flow series studied, but also support, both qualitatively and quantitatively, the low correlation dimensions reported for such. Therefore, nonlinear deterministic methods are a viable complement to linear stochastic ones for studying river flow dynamics, if sufficient caution is exercised in their applications and in interpreting the outcomes. Copyright © 2006 John Wiley & Sons, Ltd.     

Velocity and turbulence variations at the edge of saltmarshes

Saltmarsh vegetation significantly influences tidal currents and sediment deposition by decelerating the water velocity in the canopy. In order to complement previous field results, detailed profiles of velocity and turbulence were measured in a laboratory flume. Natural Spartina anglica plants were installed in a 3 m length test section in a straight, recirculating flume. Different vegetation densities, water depths and surface velocities were investigated. The logarithmic velocity profile, which existed in front of the vegetation, was altered gradually to a skimming-flow profile, typical for submerged saltmarsh vegetation. The flow reduction in the denser part of the canopy also induced an upward flow (the current was partially deflected by the canopy). The skimming flow was accompanied by a zone of high turbulence co-located with the strongest velocity gradient. This gradient moved upward and the turbulence increased with distance from the edge of the vegetation. Below the skimming flow, the velocity and the turbulence were low. The structure of the flow in the canopy was relatively stable 2 m into the vegetation. The roughness length (z₀) of the vegetation depends only on the vegetation characteristics, and is not sensitive to the current velocity or the water depth. Both the reduced turbulence in the dense canopy and the high turbulence at the top of the canopy should increase sediment deposition. On the other hand, the high turbulence zone just beyond the vegetation edge and the oblique upward flow may produce reduced sedimentation; a phenomenon that was observed near the vegetation edge in the field.     

Analytical models of contaminant transport in coastal aquifers

The Henry formulation, which couples subsurface flow and salt transport via a variable-density flow formulation, can be used to evaluate the extent of sea water intrusion into coastal aquifers. The coupling gives rise to nontrivial flow patterns that are very different from those observed in inland aquifers. We investigate the influence of these flow patterns on the transport of conservative contaminants in a coastal aquifer. The flow is characterized by two dimensionless parameters: the Péclet number, which compares the relative effects of advective and dispersive transport mechanisms, and a coupling parameter, which describes the importance of the salt water boundary on the flow. We focus our attention on two regimes – low and intermediate Péclet number flows. Two transport scenarios are solved analytically by means of a perturbation analysis. The first, a natural attenuation scenario, describes the flushing of a contaminant from a coastal aquifer by clean fresh water, while the second, a contaminant spill scenario, considers an isolated point source.     

Coarsening of three-dimensional structured and unstructured grids for subsurface flow

We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells.     

Multi-scale modeling of three-phase–three-component processes in heterogeneous porous media

Flow and transport processes in porous media occur on different spatial and temporal scales and may also be locally different. Additionally, the structure of the porous medium itself generally shows a high dependence on the spatial scale.     

Flow and deposition of pyroclastic granular flows: A type example from the 1975 Ngauruhoe eruption,New Zealand

Small-volume pyroclastic density currents (PDCs) are generated frequently during explosive eruptions with little warning. Assessing their hazard requires a physical understanding of their transport and sedimentation processes which is best achieved by the testing of experimental and numerical models of geophysical mass flows against natural flows and/or deposits. To this end we report on one of the most detailed sedimentological studies ever carried out on a series of pristine small-volume PDC deposits from the 1975 eruption of Ngauruhoe volcano, whose emplacement were also witnessed during eruption. Using high-resolution GPS surveys, a series of lateral excavations across the deposits, and bulk sedimentological analysis we constrained the geomorphology, internal structure and texture of the deposits with respect to laterally varying modes of deposition.     

Magmatic structures in the Krkonoše–Jizera Plutonic Complex,Bohemian Massif: evidence for localized multiphase flow and small-scale thermal–mechanical instabilities in a granitic magma chamber

The present paper examines magmatic structures in the Jizera and Liberec granites of the Krkonoše–Jizera Plutonic Complex, Bohemian Massif. The magmatic structures are here interpreted to preserve direct field evidence for highly localized magma flow and other processes in crystal-rich mushes, and to capture the evolution of physical processes in an ancient granitic magma chamber. We propose that after chamber-wide mixing and hybridization, as suggested by recent petrological studies, laminar magma flow became highly localized to weaker channel-like domains within the higher-strength crystal framework. Mafic schlieren formed at flow rims, and their formation presumably involved gravitational settling and velocity gradient flow sorting coupled with interstitial melt escape. Local thermal or compositional convection may have resulted in the formation of vertical schlieren tubes and ladder dikes whereas subhorizontal tubes or channels formed during flow driven by lateral gradients in magma pressure. After the cessation or deceleration of channel flow, gravity-driven processes (settling of crystals and enclaves, gravitational differentiation, development of downward dripping instabilities), accompanied by compaction, filter pressing and melt segregation, dominated in the crystal mush within the flow channels. Subsequently, magmatic folds developed in schlieren layers and the magma chamber recorded complex, late magmatic strains at high magma crystallinities. Late-stage magma pulsing into localized submagmatic cracks represents the latest events of magmatic history of the chamber prior to its final crystallization. We emphasize that the most favorable environments for the formation and preservation of magmatic structures, such as those hosted in the Jizera and Liberec granites, are slowly cooling crystal-rich mushes. Therefore, where preserved in plutons, these structures may lend strong support for a “mush model” of magmatic systems.