Gravity Segregation of Miscible Fluids in Linear Models
Authors: Gardner G.H.F. , Downie J. , Kendall H.A. 
Published: December 2007

Abstract
Some cases of the motion of two miscible fluids in uniform linear models are discussed. There is no bulk flow through the models, and the convection currents are caused solely by density gradients. Horizontal, vertical and inclined positions of the models are treated separately. Analytical formulas for the motion of a narrow transition zone between the fluids are obtained and compared with measurements, and the limitations are discussed. Observations of the ultimate stabilizing effect of diffusion on the motion are compared with the theory. It is shown that fluid motion in an inclined model with a constant rate of bulk flow can be directly deduced from the cases discussed, provided the viscosity ratio is small.

Introduction
In some reservoirs, convection currents caused by variations in the fluid densities may be significant. For practical purposes it may be desirable either to minimize the effect, particularly in miscible displacement, or to maximize the effect as in some schemes for the production of heavy crude oils by combustion. Experiments by Craig, et al, have shown the effectiveness of gravity in reducing the efficiency of a miscible displacement when the mixed zone between the displacing and displaced fluids is narrow compared with the height of the system. The presence of a board transition zone from one fluid to the other has been shown by Perkins, et al, to nullify considerably the influence of a difference in density. Indeed, the theoretical discussion given by Perrine has shown that when the transition zone exceeds a certain critical length the efficiency approaches 100 per cent. Many additional results have been published on miscible and immiscible displacements which show the importance of density difference and also mobility ratio. The aim of the present paper is to discuss some cases of gravity segregation of miscible fluids in linear models such as are often used in the laboratory to evaluate recovery processes, and to supply additional experimental evidence in support of the conclusions. An attempt is made to calculate simply the magnitude of the fluid motion caused by density gradients and mobility ratios, and to quantify the modifications in the fluid motion produced by the combined effect of molecular and convective dispersion. However, it should be pointed out that the uniform permeability of the models and the use of only a single pair of fluids (whereas several pairs of fluids may be employed in a recovery scheme) severely restricts direct conclusions about natural reservoirs. For conciseness, the motion generated by density gradients is emphasized for the case of no bulk flow through the models. This is not a drastic limitation because, as shown in Appendix A, a fluid motion with a constant bulk flow through any linear model is equivalent to a fluid motion with no bulk flow in the same model inclined at a different angle to the horizontal, provided only that the density is a linear function of the viscosity of the fluid mixture.

HORIZONTAL MODEL
In this case the model is assumed to have a rectangular cross-section with a horizontal breadth b, vertical height h and a great length L. Throughout the paper it will be assumed that x-axis points along the length, the y-axis along the width and the z-axis along the height. Thus, in the present case the x, y plane is horizontal. Initially, the transition zone between the fluids is assumed to be thin and parallel to the y, z plane. An example of this situation in a practical problem arises in a scheme for using air for cushion gas in the storage of fuel gas in aquifers. With the passage of time the denser fluid sinks and spreads along the bottom of the model while the lighter fluid rises. The interface extends from the top plane of the model to the bottom. Experiments have shown that for practical purposes this interface may be treated as a plane surface at all times. The surface is significantly curved, however, for viscosity ratios greater than about 10. Other complications arise with large viscosity ratios because of the nonlinear dependence of the viscosity on the composition of the mixture, so that the simple results given here are restricted to fairly low viscosity ratios.
SPEJ
P. 95^