Carbonated water flooding (CWI) increases oil production due to favorable dissolution effects and viscosity reduction. Accurate modeling of CWI performance requires a simulator with the ability to capture the true physics of such process. In this study, compositional modeling coupled with surface complexation modeling (SCM) are done, allowing a unified study of the influence in oil recovery of reduction of salt concentration in water. The compositional model consists of the conservation equations of total carbon, hydrogen, oxygen, chloride and decane. The coefficients of such equations are obtained from the equilibrium partition of chemical species that are soluble both in oleic and the aqueous phases. SCM is done by using the PHREEQC program, which determines concentration of the master species. Estimation of the wettability as a function of the Total Bound Product (TBP) that takes into account the concentration of the complexes in the aqueous, oleic phases and in the rock walls is performed. We solve analytically and numerically these equations in 1 - D in order to elucidate the effects of the injection of low salinity carbonated water into a reservoir containing oil equilibrated with high salinity carbonated water.

This paper is concerned with the study of the main wave interactions in a system of conservation laws in geochemical modeling. We study the modeling of the chemical complexes on the rock surface. The presence of stable surface complexes affects the relative permeability. We add terms representing surface complexes to the accumulation function in the model presented in [14]. This addition allows to take into account the interaction of ions with the rock surface in the modeling of the oil recovery by the injection of carbonated water. Compatibility hypotheses with the modeling are made on the coefficients of the system to obtain meaningful solutions. We developed a Riemann solver taking into account the complexity of the interactions and bifurcations of nonlinear waves. Such bifurcations occur at the inflection and resonance surfaces. We present the solution of a generalized eigenvalue problem in a (n+1)-dimensional case, which allows the construction of rarefaction curves. A method to find the discontinuous solutions is also presented. We find the solution path for some examples.

We study systems of partial differential-algebraic equations (PDAEs) of first order. Classical solutions of the theory of hyperbolic partial differential equation such as discontinuities (shock and contact discontinuities), rarefactions and diffusive traveling waves are extended for variables living on a surface S, which is defined as solution of a set of algebraic equations. We propose here an alternative formulation to study numerically and theoretically the PDAEs by changing the algebraic conditions into partial differential equations with relaxation source terms (PDREs). The solution of such relaxed systems is proved to tend to the surface S, i.e., to satisfy the algebraic equations for long times. We formulate a unified numerical scheme for systems of PDAEs and PDREs. This scheme is naturally parallelizable and has faster convergence. We do not perform a rigorous analysis about the convergence or accuracy for the method, the evidence of its effectiveness is presented by means of simulations for physical and synthetical problems.