Pore-Scale Simulation of the Interplay between Wettability, Capillary Number, and Salt Dispersion on the Efficiency of Oil Mobilization by Low-Salinity Waterflooding

Summary Nonuniform mixing during low-salinity waterflooding (LSWF) is a function of the pore geometry and flow patterns within the porous system. Salinity-dependent wettability alteration (WA) changes the entry capillary pressure, which may mobilize the trapped oil depending on the flow regime and salt dispersion pattern. The complex interplay between the wettability, capillary number (NCa), and salt dispersion caused by pore-scale heterogeneity on the efficiency of LSWF is not well understood. In this paper, direct numerical simulations in a pore-doublet model (PDM) were carried out with OpenFOAM® (OpenCFD, Berkshire, UK) using the volume-of-fluid (VOF) method. Oil trapping and remobilization were studied at relevant NCa as low as 10−6 under different initial wettability states. Depending on the effective salinity ranges (ESRs) for the low-salinity effect (LSE), three WA models were implemented, and the effects of WA degree and salinity distribution on LSWF flow dynamics were investigated. The slow process of WA by means of thin film phenomena was captured by considering a diffuse interface at the three-phase contact line. Because of the pore structure of the pore doublet, only in nonwater-wet cases, oil is trapped in the narrower side channel (NSC) after high-salinity waterflooding (HSWF) and may be remobilized by LSWF. In strongly oil-wet cases, oil is recovered gradually by LSWF by means of a film-flow mechanism near the outlet. In moderately oil-wet cases, however, the entire trapped oil ganglion can be mobilized, provided that the entry capillary pressure is sufficiently reduced. The degree of WA, ESR, kinetics of WA, and the wettability of pore surface at the outlet are determining factors in the drainage of the trapped oil. The salt dispersion pattern in the flowing region [i.e., wider side channel (WSC)] controls the wettability distribution and the rate and magnitude of oil recovery from the stagnant region (i.e., NSC). The difference between the WA models is more apparent near the outlet, where the salinity profile is more dispersed. The ESR in which WA occurs determines the speed of the entry capillary pressure reduction and, thus, the recovery factor. In cases where WA occurs at a salinity threshold (ST), the highest recovery is obtained, whereas with the full-salinity-range WA model, the oil recovery performance is lowest. From the capillary desaturation perspective, it is found that the LSE becomes more pronounced when NCa is less than 10−5, and the dispersion regime is in the power-law interval. Because the adverse effect of salt dispersion in the flowing region is delayed, the LSE is intensified. For the simulations to be representative of the actual conditions in the porous medium, much lower NCa than currently used in many research works must be studied. Otherwise, the simulations may lead to over- or underestimation of the LSE. The synergetic or antagonistic effects caused by the interplay between viscous and capillary forces and dispersion may lead to total recovery or entrapment of oil, regardless of WA. Based on the pore geometry, initial wettability state, and balance of forces, the mobilized oil may flow past the conjunction (favorable) or in the backward direction (unfavorable) to the WSC and get retrapped. Successful drainage of oil from the pore system after WA is essential for observing incremental oil recovery by LSWF.