Two-Phase Flow of Volatile Hydrocarbons
Authors: Kniazeff V.J. , Naville S.A. 
Published: December 2007

The problem of unsteady-state condensate-gas flow through porous media leads to a set of second-order non-linear partial differential equations. Such a set of equations is numerically solved in the case of radial two-phase flow around a well, taking into consideration both the thermodynamical properties of the fluid and the mechanical properties of the reservoir. The fluid properties, reflecting the PVT relationship of the gaseous and liquid phases, are expressed by using the partial specific masses of the two main separator products in these phases. The flow properties of the reservoir rock are expressed by the generalized Darcy's law for the liquid phase and by a quadratic relationship between the rate of flow and the pressure gradient for the gaseous phase. The numerical solution of the equations for pressure and saturation vs radius and time is worked out through programs written for a computer. The evolution of bottom-hole pressures, well productivities or deliverabilities and effluent compositions with the depletion of the reservoir is easily derived. The application to the Saharian gas-condensate field Hassi R'Mel led to a better understanding of the drainage mechanism. A zone of fairly high liquid saturation develops around the wells, reducing the effective permeability, and represents a loss of condensible products in addition to the PVT-like retrograde condensation. Inside this zone, near the well, the deviation from Darcy's law for the flow of the gaseous phase governs the well deliverability. A back-pressure test has been computed and correlates with the field results.

Introduction

Two-phase flow of volatile hydrocarbons, like condensate gas or light crude oil, may be treated as the flow of a binary mixture by an arbitrary division of the chemical components into two groups. This is translated into two equations of mass continuity, which constitute a set of relationships for the pressure and the saturations vs the space coordinates and the time. The equations contain the laws governing the composition and the motion of the phases. The problem so defined is solved with the assumption that the compositions of the phases at any pressure are respectively the same as those observed in a PVT measuring cell under differential liberation. In a first series of computations, it was assumed that the flow obeys the generalized Darcy's law. A satisfactory representation of the retrograde condensation around the well was thus obtained. In addition, the trend toward decreasing effective permeabilities was obtained, and the computed composition of the effluent checked the laboratory values. However, it has not been possible within this basic assumption to represent the non-linear relationship between the production rate and the bottom-hole pressure drawdown as observed for gas wells in the field. Following the advice of M. A. Houpeurt it was decided to consider the relative permeability to gas as a function of the velocity of the gas phase. The necessary physical determinations were made by E. Costaseque using a method devised by Messieurs A. Houpeurt and R. Iffy. As numerical processing of the equations progressed, several difficulties were encountered which were overcome through collaboration with the computer manufacturer. This mathematical model of two-phase flow in porous media had been primarily intended for and extensively applied to the case of the Hassi R'Mel gas-condensate field, operated by SN Repal for SEHR, a joint subsidiary of SN Repal and CFP (A). The programs have also found their applications to forecast the behavior of several fields in the Sahara area containing light and volatile hydrocarbons.