Aziz Rimsha, Joekar-Niasar Vahid, Martínez-Ferrer Pedro J, Godinez-Brizuela Omar E, Theodoropoulos Constantinos, Mahani Hassan
University of Manchester, School of Chemical Engineering and Analytical science, Sackville St, Manchester, M139PL, United Kingdom.
Manchester Metropolitan University, Centre for Mathematical Modelling and Flow Analysis, Chester Street, Manchester, M1 5GD, United Kingdom.
Sci Rep. 2019 Jun 25;9(1):9257. doi: 10.1038/s41598-019-45434-2.
Low salinity waterflooding has proven to accelerate oil production at core and field scales. Wettability alteration from a more oil-wetting to a more water-wetting condition has been established as one of the most notable effects of low salinity waterflooding. To induce the wettability alteration, low salinity water should be transported to come in contact with the oil-water interfaces. Transport under two-phase flow conditions can be highly influenced by fluids topology that creates connected pathways as well as dead-end regions. It is known that under two-phase flow conditions, the pore space filled by a fluid can be split into flowing (connected pathways) and stagnant (deadend) regions due to fluids topology. Transport in flowing regions is advection controlled and transport in stagnant regions is predominantly diffusion controlled. To understand the full picture of wettability alteration of a rock by injection of low salinity water, it is important to know i) how the injected low salinity water displaces and mixes with the high salinity water, ii) how continuous wettability alteration impacts the redistribution of two immiscible fluids and (ii) role of hydrodynamic transport and mixing between the low salinity water and the formation brine (high salinity water) in wettability alteration. To address these two issues, computational fluid dynamic simulations of coupled dynamic two-phase flow, hydrodynamic transport and wettability alteration in a 2D domain were carried out using the volume of fluid method. The numerical simulations show that when low salinity water was injected, the formation brine (high salinity water) was swept out from the flowing regions by advection. However, the formation brine residing in stagnant regions was diffused very slowly to the low salinity water. The presence of formation brine in stagnant regions created heterogeneous wettability conditions at the pore scale, which led to remarkable two-phase flow dynamics and internal redistribution of oil, which is referred to as the "pull-push" behaviour and has not been addressed before in the literature. Our simulation results imply that the presence of stagnant regions in the tertiary oil recovery impedes the potential of wettability alteration for additional oil recovery. Hence, it would be favorable to inject low salinity water from the beginning of waterflooding to avoid stagnant saturation. We also observed that oil ganglia size was reduced under tertiary mode of low salinity waterflooding compared to the high salinity waterflooding.
低盐度水驱已被证明能在岩心和油藏尺度上加速原油生产。从更亲油状态转变为更亲水状态的润湿性改变已被确认为低盐度水驱最显著的效果之一。为引发润湿性改变,低盐度水应被输送至与油水界面接触。在两相流条件下的输送会受到流体拓扑结构的强烈影响,流体拓扑结构会形成连通通道以及死端区域。众所周知,在两相流条件下,由于流体拓扑结构,被一种流体填充的孔隙空间可被划分为流动区域(连通通道)和停滞区域(死端)。流动区域的输送受平流控制,停滞区域的输送主要受扩散控制。为全面了解通过注入低盐度水实现岩石润湿性改变的情况,了解以下几点很重要:i)注入的低盐度水如何驱替并与高盐度水混合;ii)持续的润湿性改变如何影响两种不混溶流体的重新分布;以及iii)低盐度水与地层盐水(高盐度水)之间的流体动力输送和混合在润湿性改变中的作用。为解决这两个问题,使用流体体积法对二维区域内耦合的动态两相流、流体动力输送和润湿性改变进行了计算流体动力学模拟。数值模拟表明,注入低盐度水时,地层盐水(高盐度水)通过平流从流动区域被驱替出来。然而,滞留在停滞区域的地层盐水向低盐度水的扩散非常缓慢。停滞区域中地层盐水的存在在孔隙尺度上产生了非均一的润湿性条件,这导致了显著的两相流动动力学和油的内部重新分布,这被称为“推拉”行为,此前文献中尚未涉及。我们的模拟结果表明,三次采油中停滞区域的存在阻碍了通过润湿性改变实现额外原油采收的潜力。因此,从水驱开始就注入低盐度水以避免停滞饱和度是有利的。我们还观察到,与高盐度水驱相比,在低盐度水驱的三次采油模式下油滴尺寸减小。