8 OMEKEH, EVJE, AND FRIISfor a dynamic change from an initial high salinity state towards a low salinity state controlled bythe degree of desorption of calcium and magnesium from the core.3.3. Modeling of transition from high salinity to low salinity conditions. In this work, wefollow along the lines of [14, 15], and model wettability alteration by defining relative permeabilitycurves through an interpolation between the high salinity and low salinity curves given by (14),combined with the data specified in Appendix B as given by (97) and (98).More precisely, motivated by the proposed hypothesis that transition from a high salinity wettingstate towards low salinity conditions depends on the desorption of calcium and magnesiumcaused by the cation exchange process, the following interpolation is proposed:(19) k(s, β ca , β mg ) = H(β ca , β mg )k HS (s) + [1 − H(β ca , β mg )]k LS (s),where H(β ca , β mg ) is defined by (18). Hence, when no desorption of divalent cations from the claysurface takes place it follows that H(β ca , β mg ) = 1, implying that k(β ca , β mg ) = k HS (s). Thisreflects the initial high salinity wetting state. Then, as desorption of divalent cations takes place,it follows that m(β ca , β mg ) increases. In particular, if the desorption effects becomes large enough,it follows from the above discussion (see equation (18)) that H(β ca , β mg ) ≈ 0, which means thatk(s, β ca , β mg ) ≈ k LS (s), reflecting that a wettability alteration has taken place which results in alow salinity wetting state.4. The coupled model for water-oil flow and multiple ion exchangeWe now want to take into account convective and diffusive forces associated with the brine aswell as the oil phase. In order to include such effects we must consider the following equations forthe total concentrations ρ o , ρ l , ρ ca , ρ so , ρ mg , ρ na , ρ cl (mol per liter core):(20)∂ t ρ o + ∇ · (ρ o v o ) = 0,∂ t ρ l + ∇ · (ρ l v l ) = 0,∂ t ρ na + ∂ t (M c β na ) + ∇ · (ρ na v g ) = 0,∂ t ρ cl + ∇ · (ρ cl v g ) = 0,∂ t ρ ca + ∂ t (M c β ca ) + ∇ · (ρ ca v g ) = 0,∂ t ρ so + ∇ · (ρ so v g ) = 0,∂ t ρ mg + ∂ t (M c β mg ) + ∇ · (ρ mg v g ) = 0,(oil flowing through the pore space)(water flowing through the pore space)(Na + -ions in water)(Cl − -ions in water)(Ca 2+ -ions in water)(SO 2−4 -ions in water)(Mg 2+ -ions in water).Here v o , v l and v g are, respectively, the oil, water and species ”fluid” velocities, whereas M crepresents the mass of the clay. The subsequent derivation of the model closely follows the work[14, 15], however, the water-rock chemistry in our current model is given in terms of a multipleion exchange process (MIE), not dissolution/precipitation as in [14, 15].Let s o denote the oil saturation, i.e. the fraction of volume of the pore space represented byporosity ϕ that is occupied by the oil phase, and s the corresponding water saturation. Thetwo saturations are related by the basic relation s o + s = 1. Furthermore, we define the porousconcentration C o associated with the oil component as the concentration taken with respect to thevolume of the pore space occupied by oil and represented by ϕs o . Hence, C o and ρ o are related by(21) ρ o = ϕs o C o .Similarly, the porous concentrations of the various components in the water phase are defined as theconcentration taken with respect to the volume of the pores occupied by water ϕs. Consequently,the porous concentrations C l , C na , C cl , C ca , C mg , and C so are related to the total concentrationsbyρ l = ϕsC l , ρ na = ϕsC na , ρ cl = ϕsC cl , ρ ca = ϕsC ca , ρ mg = ϕsC mg , ρ so = ϕsC so .Following, for example, [1, 2], we argue that since oil, water, and the ions in water Na + , Cl − ,Ca 2+ , Mg 2+ , and SO 2−4 flow only through the pores of the calcite specimen, the ”interstitial”velocity v o associated with the oil, v l associated with the water, and v g associated with theions, have to be defined with respect to the concentrations inside the pores, and differ from the
<strong>MODEL<strong>IN</strong>G</strong> <strong>OF</strong> <strong>LOW</strong> <strong>SAL<strong>IN</strong>ITY</strong> <strong>EFFECTS</strong> <strong>IN</strong> <strong>SANDSTONE</strong> <strong>OIL</strong> ROCKS 9respective seepage velocities V o , V l and V g . The velocities are related by the Dupuit-Forchheimerrelations, see [1] and references therein,(22) V o = ϕs o v o , V l = ϕsv l , V g = ϕsv g .Consequently, the balance equations (20) can be written in the form(23)∂ t (ϕs o C o ) + ∇ · (C o V o ) = 0,∂ t (ϕsC l ) + ∇ · (C l V l ) = 0,∂ t (ϕsC na ) + ∂ t (M c β na ) + ∇ · (C na V g ) = 0,∂ t (ϕsC cl ) + ∇ · (C cl V g ) = 0,∂ t (ϕsC ca ) + ∂ t (M c β ca ) + ∇ · (C ca V g ) = 0,∂ t (ϕsC so ) + ∇ · (C so V g ) = 0,∂ t (ϕsC mg ) + ∂ t (M c β mg ) + ∇ · (C mg V g ) = 0.In order to close the system we must determine the seepage velocities V o , V l and V g . For thatpurpose we consider the concentration of the water phase (brine) C that occupies the pore spaceas a mixture of water C l and the various species Na + , Cl − , Ca 2+ , Mg 2+ , and SO 2−4 representedby C g . In other words,(24) C g = C na + C cl + C ca + C mg + C so , C = C g + C l .Then, we define the seepage velocity V associated with C by(25) CV := C g V g + C l V l .Now we are in a position to rewrite the model in terms of V and the diffusive velocity U g givenby(26) U g = V g − V.Then the model (23) takes the form(27)∂ t (ϕs o C o ) + ∇ · (C o V o ) = 0,∂ t (ϕsC l ) + ∇ · (C l V l ) = 0,∂ t (ϕsC na ) + ∂ t (M c β na ) + ∇ · (C na U g ) = −∇ · (C na V),∂ t (ϕsC cl ) + ∇ · (C cl U g ) = −∇ · (C cl V),∂ t (ϕsC ca ) + ∂ t (M c β ca ) + ∇ · (C ca U g ) = −∇ · (C ca V),∂ t (ϕsC so ) + ∇ · (C so U g ) = −∇ · (C so V),∂ t (ϕsC mg ) + ∂ t (M c β mg ) + ∇ · (C mg U g ) = −∇ · (C mg V).Furthermore, we can assume that the seepage velocity V associated with the water phase representedby C, is given by Darcy’s law [1, 7, 26](28) V = −κλ(∇p − ρg∇d), λ = k µ ,where κ is absolute permeability, k is water relative permeability, and µ is viscosity, and p pressurein water phase. Similarly, for the oil phase(29) V o = −κλ o (∇p o − ρ o g∇d), λ o = k oµ o.The diffusive velocity U g is expressed by Fick’s law by(30) C i U g = −D∇C i , i = na, cl, ca, so, mg,where D is the diffusion coefficient. In view of (24) and (30), it follows that(31) C g U g = −D∇C g .Note that we assume that the diffusion coefficient D is the same for all species i = na, cl, ca, so, mg.This is a reasonable assumption as long as the concentration is not too high. D is the diffusion