6 OMEKEH, EVJE, AND FRIISBeta Cabeta Cax 10 −34.543.532.521.510.50.10.20.30.40.5Mg2+ (mole/liter)0.60.70.20.410.80.6Ca2+ (mole/liter)1.21.4Figure 1. Plot showing β ca as a function of varying C ca and C mg concentrations,and with C na = 0.15 (moles per liter). Other parameters in (12) like K cana ,K mgna , and CEC, are listed in Appendix A.(12) β ca (C na , C ca , C mg ) =and(13) β mg (C na , C ca , C mg ) =K cana√γca C ca CEC2K cana√γca C ca + 2K mgna√γmg C mg + γ na C na,K mgna√γmg C mg CEC2K cana√γca C ca + 2K mgna√γmg C mg + γ na C na.We note that the equations (11), (12) and (13) are equivalent to a Langmuir-type adsorptionisotherm. Fig. 1 illustrates how the β ca function depends on the concentrations C na , C ca , C mg .At high magnesium concentration, the amount of calcium ion adsorbed on the rock, β ca , becomesquite low.3. Coupling of wettability alteration to changes on the rock surfaceThis section discusses aspects concerning the flow functions. Since we only consider flows withoutcapillary pressure in Section 6, we limit the discussion to the relative permeability functions.The following ideas are also employed in [14, 15], however, in the context of spontaneous imbibitionon chalk cores where capillary forces are the driving forces in oil recovery. It also partially followsideas employed in previous works in [22, 41, 45].3.1. Relative permeability functions. As a basic model for relative permeability functionsthe well-known Corey type correlations are used [12]. They are given in the form (dimensionlessfunctions)(14)k(s) = k ∗( s − s)wrNk,swr ≤ s ≤ 1 − s or ,1 − s or − s wr( 1 −k o (s) = ko∗ sor − s) Nko,swr ≤ s ≤ 1 − s or ,1 − s or − s wrwhere s wr and s or represent critical saturation values and Nk and Nk o are the Corey exponentsthat must be specified. In addition, k ∗ and k ∗ o are the end point relative permeability values thatalso must be given. Now, we define two extreme relative permeability functions corresponding tothe wetting state of the rock for high salinity and low salinity conditions.
<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 7Relative permeability0.90.80.70.60.50.40.30.2krw−HSkro−HSkrw−LSkro−LSfractional flow10.80.60.40.2HSLS0.100 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1Water Saturation00 0.2 0.4 0.6 0.8 1Water SaturationFigure 2. Left: Example of relative permeability curves corresponding to highsalinity(HS) and low-salinity (LS) conditions. Right: Corresponding fractionalflow functions.3.1.1. High salinity conditions. This is assumed to be the initial state of the core. The valuesfor generating the functions are listed in Appendix B. Applying the values of the high salinitycondition in (14) the following two relative permeability curves are obtained(15) k HS (s; s HSwr , s HSor ),koHS (s; s HSwr , s HSor ),We refer to Fig. 2 for a plot of these curves (red line).for s HSwr≤ s ≤ 1 − s HSor .3.1.2. Low salinity conditions. This wetting condition is assumed to be attained only whenthere is complete desorption of Ca 2+ and Mg 2+ ions from the rock surface. The values for generatingthe functions are listed in Appendix B. Applying the values of the low salinity condition in(14) gives corresponding relative permeability curves(16) k LS (s; s LSwr , s LSor ),koLS (s; s LSwr , s LSor ),for s LSwr ≤ s ≤ 1 − s LSor .We refer to Fig. 2 for a plot of these curves (blue line). The motivation for the choice of the valueswas to match the form of the relative permeability measured for a variety of high salinity and lowsalinity brines in [43].3.2. Cation exchange as a mechanism for wettability alteration. We let β ca0 and β mg0be the amount of calcium and magnesium, respectively, initially bounded to the clay surface. Wethen define the quantity(17) m(β ca , β mg ) := max(β ca0 − β ca , 0) + max(β mg0 − β mg , 0),as a measure for the desorption of cations from the clay. Moreover, we define(18) H(β ca , β mg ) :=11 + rm(β ca , β mg ) ,where r > 0 is a specified constant. Note that the choice of r determines the extent in which thedivalent ion desorbed leads to a certain change of the wetting state.The function H(β ca , β mg ) is a weighting function, and works such that H = 1, when there is nodesorption of calcium and magnesium from the rock, whereas 0 < H < 1 in case of desorption ofat least one of these cations. How fast H(β ca , β mg ) is approaching 0 as m(β ca , β mg ) is increasing,depends on the choice of r. Now, the weighting function H(β ca , β mg ) can be used to represent thewetting state in the core plug; H(β ca , β mg ) = 1 corresponds to the initial oil-wet state, whereasH(β ca , β mg ) ≈ 0 represents the water-wet state. By defining relative permeability curves by meansof the weighting function H(β ca , β mg ) as described in the next subsection, the model can account