Network models for two-phase flow in porous media Part 1. Immiscible microdisplacement of non-wetting fluids

• Type: Article in International Scientific Journal
• Date: 1986

Authors

Madalena M. Dias (Author) FEUP; Payatakes, AC (Author) Other

Journal

Title: Journal of Fluid Mechanics

Vol. 164

Pages: 305-336

ISSN: 0022-1120

Publisher: Cambridge University Press

Other information

Authenticus ID: P-00R-KJ6

DOI: 10.1017/s0022112086002574

Abstract (EN): <jats:p>A theoretical simulator of immiscible displacement of a non-wetting fluid by a wetting one in a random porous medium is developed. The porous medium is modelled as a network of randomly sized unit cells of the constricted-tube type. Under creeping-flow conditions the problem is reduced to a system of linear equations, the solution of which gives the instantaneous pressures at the nodes and the corresponding flowrates through the unit cells. The pattern and rate of the displacement are obtained by assuming quasi-static flow and taking small time increments. The porous medium adopted for the simulations is a sandpack with porosity 0.395 and grain sizes in the range from 74 to 148 ¿rn. The effects of the capillary number, <jats:italic>Ca</jats:italic>, and the viscosity ratio, ¿ = ¿<jats:sub>o</jats:sub>/¿<jats:sub>w</jats:sub>, are studied. The results confirm the importance of the capillary number for displacement, but they also show that for moderate and high Ca values the role of ¿ is pivotal. When the viscosity ratio is favourable (¿ < 1), the microdisplacement efficiency begins to increase rapidly with increasing capillary number for <jats:italic>Ca</jats:italic> > 10<jats:sup>¿5</jats:sup>, and becomes excellent as <jats:italic>Ca</jats:italic> ¿ 10<jats:sup>¿3</jats:sup>. On the other hand, when the viscosity ratio is unfavourable (¿ > 1), the microdisplacement efficiency begins to improve only for Ca values larger than, say, 5 × 10<jats:sup>¿4</jats:sup>, and is substantially inferior to that achieved with ¿ < 1 and the same <jats:italic>Ca</jats:italic> value. In addition to the residual saturation of the non-wetting fluid, the simulator predicts the time required for the displacement, the pattern of the transition zone, the size distribution of the entrapped ganglia, and the acceptance fraction as functions of <jats:italic>Ca</jats:italic>, ¿, and the porous-medium geometry.</jats:p>

Language: English

Type (Professor's evaluation): Scientific