Artigo Completo:

We present a locally conservative numerical methodology to simulate the two-phase flow (water and CO2) with mass absorption between the fluid phases and reaction between the CO2 phase and rock in a homogeneous reservoir. This problem is modeled by a system of partial differential equations, which basically consists in a parabolic subsystem for determining the velocity field and two non-linear hyperbolic equations for the transport of phases that flow (equations of saturation and concentration). From the numerical point of view, we use the operator splitting technique to properly treat the time scale of each physical phenomenon. We propose the application of a locally conservative finite element method for the total Darcy velocity and a high-order non-oscillatory central- scheme finite volume method for nonlinear hyperbolic equations that govern the saturation and concentration of phases. Furthermore, we treat numerically the mass flux between fluid phases, the dissolution of CO2 in the aqueous phase, using the flash methodology that treats numerically equilibrium reactions. The reaction of CO2 with rock (precipitation), which causes changes in porosity and permeability, was treated by applying principles of kinetic theory.