ACS PRF | ACS | All e-Annual Reports

Reports: AC9

Back to Table of Contents

45169-AC9
Physical Foundations of the Vibratory Recovery of Residual Oil

Igor Beresnev, Iowa State University

            To better predict the mobilization of oil entrapped in porous media, the process of disintegration of oil into droplets should be thoroughly understood.  We have continued building such an understanding in two directions: (1) by refining the analytical models of the evolution of an oil/water interface to a complete snap-off and (2) by validating these models against computational fluid dynamics (CFD).  We are now at the stage at which the analytical and computational results have closely matched each other, which attests to the correct theoretical description achieved. 

            Technically, this is a hydrodynamic problem of the dynamics of an interface between a non-wetting core fluid (such as oil) surrounded by a wetting annulus (water) in a constricted capillary tube.  We are interested in describing the evolution of the interface toward its end configurations as the oil phase moves through pore constrictions.  As established during Year 1 of the project, the geometry of the constrictions plays pivotal role in whether the oil continues to flow as a stream or disintegrates into fluid particles; the resulting state of the oil phase has fundamental significance for the methods of its recovery, including the vibratory stimulation. 

            The theory reported as part of the project implementation during Year 1 has now been significantly improved and validated.  Whereas the earlier theory had to assume comparable thicknesses of the oil core and water annulus, as well as a much greater viscosity of the oil, the present theory is valid for arbitrary viscosities and thicknesses.  Another milestone is the verification of the theory in an accurate computational-fluid-dynamics experiment. 

            The key elements that led to a partial-derivative equation for the evolution of the oil/water interface are as follows.  In the surface-tension-driven flow, the fluids in the annulus and the core move in opposite directions, and the conservation of mass dictates that Qwater = –Qoil, where Q is the volume flow rate.  The former condition means that the tangential fluid velocity at the interface is zero.  With the small-slope assumption, the fluid flow in the core and the annulus can be approximated by Poiseuille flow.  Combining the conservation of mass, Laplace's law for capillary pressure, and the equations for Poiseuille flow in the core and the annulus, we arrive at the evolution equation for the radial position of the interface κ(x, τ) as a function of the dimensionless axial coordinate x and time τ:

,                      (1)

where μwater/μoil is the viscosity ratio, and λ(x) is the shape of the pore wall.  Equation (1) is a much more general, unrestricted formulation of the quantitative interface dynamics than the simplified version reported during Year 1.  Its solutions are expected to better describe the reality. 

            The validity of the numerical solutions of equation (1) has been verified using computational fluid dynamics.  The commercial code FLUENT was used to compute the evolution of the interfaces for the same initial conditions as in the numerical solution of (1), for a variety of different scenarios. 

            Figure 1 shows a typical comparison between the interface shape predicted by the theory and that obtained by the fully dynamic CFD simulations.  The scenario uses a sinusoidal shape of the wall; the initial oil/water interface is parallel to it.  The evolution from the initial interface configuration to the break-up, obtained through numerical integration of equation (1), is shown in Figure 1, upper panel.  The evolution is from the initial profile in the “back” toward the viewer; the end profile just before it snaps off is at the “front” of the plot.  The formation of a large satellite drop in the neck of the constriction is predicted.  The lower panel shows the same final interface shape obtained through computational fluid dynamics (water is yellow and oil is red).  It is very close to the shape predicted by the evolution equation, and a satellite drop forms as well.  

            The close match between the theory and the CFD simulations indicates that the theoretical model correctly captures all essential features of the evolution of the oil/water interface driven by capillary forces.  A correct prediction of the sizes of the oil drops forming in porous channels is an important prerequisite for the calculation of the possibility of the drops' mobilization by vibrations.  The theory of the break-up can therefore be used in conjunction with the theory of vibratory stimulation to provide a complete description of the mobilization process. 

Figure 1.  Oil/water interface profiles just before the break-up.  Upper: numerical integration of the evolution equation, lower: computational-fluid-dynamics simulation

Back to top