Low salinity water flooding: Evaluating the effect of salinity on oil and water relative permeability curves using coupling of DLVO and geochemical reactions

Abstract

Enhancing the oil recovery of carbonate reservoirs with acidic oil depends on injected brine composition and salinity throughout the low salinity water flooding. The typical description of this phenomenon on oil and water relative permeability curves is being put forward by the process-dependent interpolation function, which barely addresses the underlying fundamentals governing the direct impact of wettability alteration and geochemical reactions on flow saturation functions. The main objective of this work is to perform an integrated investigation of the effect of salinity on oil and water relative permeability curves by dynamic coupling of the Derjaguin, Landau, Verwey, and Overbeek (DVLO) theory and geochemical reactions. We studied both the single and two-phase systems using the basis of the geochemical reactions to track the variations of the wetting and non-wetting phase relative permeability curves. In the single-phase scenario, the results show the superiority of the reactive flow modeling technique with kinetic-controlled reactions, in the range between [2–4] of the Damkohler number, to yield a better result for the experimental effluent ions' concentrations than the equilibrium approach. In addition, the single-phase results under dynamic conditions propose Civan correlation as the best predictor of residual resistance factor (RRF) by numerically adjusting the correlation of coefficients based on the physical variables - salinity, dimensionless consumed ion concentration, and Peclet number. Accordingly, the results indicate a novel linkage between the single-phase and two-phase results in the sense that RRF under single-phase mode can be used to estimate the change in the wetting phase relative permeability under the two-phase mode. The single-phase methodology provides the capability for an improvement – up to 15% - in the prediction of the pressure drop curve. For two-phase cases, our procedure involved coupling the characterized geochemical reactions and the DLVO theory to consider the effect of injected water composition on contact angle as the static two-phase parameter. The phase-field method is coupled with DLVO theory to calculate the relative permeability curves for the non-wetting phase, which corresponds to the contact angle at each simulation time step. The two-phase results indicate an improvement – up to 30% - in the prediction of the recovery factor curve, on a history-matching basis, using the DLVO theory with contact angle-dependent correlation of non-wetting phase relative permeability curve. The main novelties of this study are using the single-phase dynamic permeability impairment results in calculating the variation of the wetting-phase relative permeability curve. Moreover, dynamic coupling of DLVO and non-linear contact angle-dependent correlation is used to calculate the non-wetting phase relative permeability curve. According to these novelties, more accurate prediction of the laboratory production profiles can be obtained.