Vertical integration of three-phase flow equations for analysis of light hydrocarbon plume movement

A mathematical model is derived for areal flow of water and light hydrocarbon in the presence of gas at atmospheric pressure. Vertical integration of the governing three-dimensional, three-phase flow equations is performed under the assumption of local vertical equilibrium to reduce the dimensionality of the problem to two orthogonal horizontal directions. Independent variables in the coupled water and hydrocarbon areal flow equations are specified as the elevation of zero gauge hydrocarbon pressure (air-oil table) and the elevation of zero gauge water pressure (air-water table). Constitutive relations required in the areal flow model are vertically integrated fluid saturations and vertically integrated fluid conductivities as functions of air-oil and air-water table elevations. Closed-form expressions for the vertically integrated constitutive relations are derived based on a three-phase extension of the Brooks-Corey saturation-capillary pressure function. Closed-form Brooks-Corey relations are compared with numerically computed analogs based on the Van Genuchten retention function. Close agreement between the two constitutive models is observed except at low oil volumes when the Brooks-Corey model predicts lower oil volumes and transmissivities owing to the assumption of a distinct fluid entry pressure. Nonlinearity in the vertically integrated constitutive relations is much less severe than in the unintegrated relations. Reduction in dimensionality combined with diminished nonlinearity, makes the vertically integrated water and hydrocarbon model an efficient formulation for analyzing field-scale problems involving hydrocarbon spreading or recovery under conditions for which the vertical equilibrium assumption is expected to be a satisfactory approximation.