ALBERT

Description: We consider the instantaneous release of a finite volume of fluid in a porous medium saturated with a second, immiscible fluid of different density. The resulting two-phase gravity current exhibits a rich array of behaviours due to both the residual trapping of fluid as the current recedes and the differing effects of surface tension between advancing and receding regions of the current. We develop a framework for the evolution of such a current with particular focus on the large-scale implications of the form of the constitutive relation between residual trapping and initial saturation. Pore-scale hysteresis within the current is represented by families of scanning curves relating capillary pressure and relative permeability to saturation. In the resulting vertically integrated model, all capillary effects are incorporated within specially defined saturation and flux functions specific to each region. In the long-time limit, when the height of the current and the saturations within it are low, the saturation and flux functions can be approximated by mathematically convenient power laws. If the trapping model is approximately linear at low saturations, the equations admit a similarity solution for the propagation rate and height profile of the late-time gravity current. We also solve the governing partial differential equation numerically for the nonlinear Land's trapping model, which is commonly used in studies of two-phase flows. Our investigation suggests that for trapping relations for which the proportion of trapped to initial fluid saturation increases and tends to unity as the initial saturation decreases, both of which are properties of Land's model, a gravity current slows and eventually stops. This trapping behaviour has important applications, for example to the ultimate distance contaminants or stored carbon dioxide may travel in the subsurface. © 2017 Cambridge University Press.

Description: We develop a model describing the buoyancy-driven propagation of two-phase gravity currents, motivated by problems in groundwater hydrology and geological storage of carbon dioxide (CO2). In these settings, fluid invades a porous medium saturated with an immiscible second fluid of different density and viscosity. The action of capillary forces in the porous medium results in spatial variations of the saturation of the two fluids. Here, we consider the propagation of fluid in a semi-infinite porous medium across a horizontal, impermeable boundary. In such systems, once the aspect ratio is large, fluid flow is mainly horizontal and the local saturation is determined by the vertical balance between capillary and gravitational forces. Gradients in the hydrostatic pressure along the current drive fluid flow in proportion to the saturation-dependent relative permeabilities, thus determining the shape and dynamics of two-phase currents. The resulting two-phase gravity current model is attractive because the formalism captures the essential macroscopic physics of multiphase flow in porous media. Residual trapping of CO2 by capillary forces is one of the key mechanisms that can permanently immobilize CO2 in the societally important example of geological CO2 sequestration. The magnitude of residual trapping is set by the areal extent and saturation distribution within the current, both of which are predicted by the two-phase gravity current model. Hence the magnitude of residual trapping during the post-injection buoyant rise of CO2 can be estimated quantitatively. We show that residual trapping increases in the presence of a capillary fringe, despite the decrease in average saturation. © 2011 Cambridge University Press.

Description: The effect of confining boundaries on gravity currents in porous media is investigated theoretically and experimentally. Similarity solutions are derived for currents when the volume increases as tα in horizontal channels of uniform cross-section with boundary height b satisfying b ∼ a|y/a|n, where y is the cross-channel coordinate and a is a length scale of the channel width. Experiments were carried out in V-shaped and semicircular channels for the case of gravity currents with constant volume (α=0) and constant flux (α=1). These showed generally good agreement with the theory. Typically, we find that the propagation of the current is well described by L ∼ tc for some scalar c. We study the dependence of c on the time exponent of the volume of fluid in the current,α, and the geometry of the channel, parameterized by n. For all channel shapes, there exists a critical value of α,αc = 1/2, above which increasing n causes an increase in c and below which increasing n causes a decrease in c, where increasing n corresponds to opening up the channel boundary to the horizontal. The current height increases or decreases with respect to time depending on whether is greater or less than αc. It is this fact, along with global mass conservation, which explains why varying the channel shape n affects the propagation rate c in different ways depending on α. We also consider channels inclined at an angle Øto the horizontal. When the slope of the channel is much greater than the slope of the free surface of the current, the component of gravity parallel to the slope dominates, causing the current to move with a constant velocity, Vf say, regardless of channel shape n and flux parameter α, in agreement with results for a two-dimensional gravity current obtained by Huppert & Woods (1995) and some initially axisymmetric gravity currents presented by Vella & Huppert (2006). If the effect of the component of gravity perpendicular to the channel may not be neglected, i.e. if the slopes of the channel and free surface of the current are comparable, we find that, in a frame moving with speed Vf, the form of the governing equation for the height of a current in an equivalent horizontal channel is recovered. We calculate that the height of a constant flux gravity current down an inclined channel will tend to a fixed depth, which is determined by the channel shape, n, and the physical properties of the fluid and rock. Experimental and numerical results for inclined V-shaped channels agree very well with this theory. © 2010 Cambridge University Press.