Search citation statements
Paper Sections
Citation Types
Year Published
Publication Types
Relationship
Authors
Journals
Viscous fingering in porous media is an instability which occurs when a low-viscosity injected fluid displaces a much more viscous resident fluid, under miscible or immiscible conditions. Immiscible viscous fingering is more complex and has been found to be difficult to simulate numerically and is the main focus of this paper. Many researchers have identified the source of the problem of simulating realistic immiscible fingering as being in the numerics of the process, and a large number of studies have appeared applying high-order numerical schemes to the problem with some limited success. We believe that this view is incorrect and that the solution to the problem of modelling immiscible viscous fingering lies in the physics and related mathematical formulation of the problem. At the heart of our approach is what we describe as the resolution of the “M-paradox”, where M is the mobility ratio, as explained below. In this paper, we present a new 4-stage approach to the modelling of realistic two-phase immiscible viscous fingering by (1) formulating the problem based on the experimentally observed fractional flows in the fingers, which we denote as $$ f_{\rm w}^{*} $$ f w ∗ , and which is the chosen simulation input; (2) from the infinite choice of relative permeability (RP) functions, $$ k_{\rm rw}^{*} $$ k rw ∗ and $$ k_{\rm ro}^{*} $$ k ro ∗ , which yield the same $$ f_{\rm w}^{*} $$ f w ∗ , we choose the set which maximises the total mobility function, $$ \lambda_{\text{T}}^{{}} $$ λ T (where $$ \lambda_{\text{T}}^{{}} = \lambda_{\text{o}}^{{}} + \lambda_{\text{w}}^{{}} $$ λ T = λ o + λ w ), i.e. minimises the pressure drop across the fingering system; (3) the permeability structure of the heterogeneous domain (the porous medium) is then chosen based on a random correlated field (RCF) in this case; and finally, (4) using a sufficiently fine numerical grid, but with simple transport numerics. Using our approach, realistic immiscible fingering can be simulated using elementary numerical methods (e.g. single-point upstreaming) for the solution of the two-phase fluid transport equations. The method is illustrated by simulating the type of immiscible viscous fingering observed in many experiments in 2D slabs of rock where water displaces very viscous oil where the oil/water viscosity ratio is $$ (\mu_{\text{o}} /\mu_{\text{w}} ) = 1600 $$ ( μ o / μ w ) = 1600 . Simulations are presented for two example cases, for different levels of water saturation in the main viscous finger (i.e. for 2 different underlying $$ f_{\rm w}^{*} $$ f w ∗ functions) produce very realistic fingering patterns which are qualitatively similar to observations in several respects, as discussed. Additional simulations of tertiary polymer flooding are also presented for which good experimental data are available for displacements in 2D rock slabs (Skauge et al., in: Presented at SPE Improved Oil Recovery Symposium, 14–18 April, Tulsa, Oklahoma, USA, SPE-154292-MS, 2012. 10.2118/154292-MS, EAGE 17th European Symposium on Improved Oil Recovery, St. Petersburg, Russia, 2013; Vik et al., in: Presented at SPE Europec featured at 80th EAGE Conference and Exhibition, Copenhagen, Denmark, SPE-190866-MS, 2018. 10.2118/190866-MS). The finger patterns for the polymer displacements and the magnitude and timing of the oil displacement response show excellent qualitative agreement with experiment, and indeed, they fully explain the observations in terms of an enhanced viscous crossflow mechanism (Sorbie and Skauge, in: Proceedings of the EAGE 20th Symposium on IOR, Pau, France, 2019). As a sensitivity, we also present some example results where the adjusted fractional flow ($$ f_{\rm w}^{*} $$ f w ∗ ) can give a chosen frontal shock saturation, $$ S_{\rm wf}^{*} $$ S wf ∗ , but at different frontal mobility ratios, $$ M(S_{\rm wf}^{*} ) $$ M ( S wf ∗ ) . Finally, two tests on the robustness of the method are presented on the effect of both rescaling the permeability field and on grid coarsening. It is demonstrated that our approach is very robust to both permeability field rescaling, i.e. where the (kmax/kmin) ratio in the RCF goes from 100 to 3, and also under numerical grid coarsening.
Viscous fingering in porous media is an instability which occurs when a low-viscosity injected fluid displaces a much more viscous resident fluid, under miscible or immiscible conditions. Immiscible viscous fingering is more complex and has been found to be difficult to simulate numerically and is the main focus of this paper. Many researchers have identified the source of the problem of simulating realistic immiscible fingering as being in the numerics of the process, and a large number of studies have appeared applying high-order numerical schemes to the problem with some limited success. We believe that this view is incorrect and that the solution to the problem of modelling immiscible viscous fingering lies in the physics and related mathematical formulation of the problem. At the heart of our approach is what we describe as the resolution of the “M-paradox”, where M is the mobility ratio, as explained below. In this paper, we present a new 4-stage approach to the modelling of realistic two-phase immiscible viscous fingering by (1) formulating the problem based on the experimentally observed fractional flows in the fingers, which we denote as $$ f_{\rm w}^{*} $$ f w ∗ , and which is the chosen simulation input; (2) from the infinite choice of relative permeability (RP) functions, $$ k_{\rm rw}^{*} $$ k rw ∗ and $$ k_{\rm ro}^{*} $$ k ro ∗ , which yield the same $$ f_{\rm w}^{*} $$ f w ∗ , we choose the set which maximises the total mobility function, $$ \lambda_{\text{T}}^{{}} $$ λ T (where $$ \lambda_{\text{T}}^{{}} = \lambda_{\text{o}}^{{}} + \lambda_{\text{w}}^{{}} $$ λ T = λ o + λ w ), i.e. minimises the pressure drop across the fingering system; (3) the permeability structure of the heterogeneous domain (the porous medium) is then chosen based on a random correlated field (RCF) in this case; and finally, (4) using a sufficiently fine numerical grid, but with simple transport numerics. Using our approach, realistic immiscible fingering can be simulated using elementary numerical methods (e.g. single-point upstreaming) for the solution of the two-phase fluid transport equations. The method is illustrated by simulating the type of immiscible viscous fingering observed in many experiments in 2D slabs of rock where water displaces very viscous oil where the oil/water viscosity ratio is $$ (\mu_{\text{o}} /\mu_{\text{w}} ) = 1600 $$ ( μ o / μ w ) = 1600 . Simulations are presented for two example cases, for different levels of water saturation in the main viscous finger (i.e. for 2 different underlying $$ f_{\rm w}^{*} $$ f w ∗ functions) produce very realistic fingering patterns which are qualitatively similar to observations in several respects, as discussed. Additional simulations of tertiary polymer flooding are also presented for which good experimental data are available for displacements in 2D rock slabs (Skauge et al., in: Presented at SPE Improved Oil Recovery Symposium, 14–18 April, Tulsa, Oklahoma, USA, SPE-154292-MS, 2012. 10.2118/154292-MS, EAGE 17th European Symposium on Improved Oil Recovery, St. Petersburg, Russia, 2013; Vik et al., in: Presented at SPE Europec featured at 80th EAGE Conference and Exhibition, Copenhagen, Denmark, SPE-190866-MS, 2018. 10.2118/190866-MS). The finger patterns for the polymer displacements and the magnitude and timing of the oil displacement response show excellent qualitative agreement with experiment, and indeed, they fully explain the observations in terms of an enhanced viscous crossflow mechanism (Sorbie and Skauge, in: Proceedings of the EAGE 20th Symposium on IOR, Pau, France, 2019). As a sensitivity, we also present some example results where the adjusted fractional flow ($$ f_{\rm w}^{*} $$ f w ∗ ) can give a chosen frontal shock saturation, $$ S_{\rm wf}^{*} $$ S wf ∗ , but at different frontal mobility ratios, $$ M(S_{\rm wf}^{*} ) $$ M ( S wf ∗ ) . Finally, two tests on the robustness of the method are presented on the effect of both rescaling the permeability field and on grid coarsening. It is demonstrated that our approach is very robust to both permeability field rescaling, i.e. where the (kmax/kmin) ratio in the RCF goes from 100 to 3, and also under numerical grid coarsening.
Viscous fingering in porous media occurs when the (miscible or immiscible) displacing fluid has a lower viscosity than the displaced fluid. For example, immiscible fingering is observed in experiments where water displaces a much more viscous oil. Modelling the observed fingering patterns in immiscible viscous fingering has proven to be very challenging, which has often been identified as being due to numerical issues. However, in a recent paper (Sorbie et al. in Transp. Porous Media 133:331–359, 2020) suggested that the modelling issues are more closely related to the physics and formulation of the problem. They proposed an approach based on the fractional flow curve, $${f}_{w}^{*}$$ f w ∗ , as the principal input, and then derived relative permeabilities which give the maximum total mobility. Sorbie et al. were then able to produce complex, well-resolves immiscible finger patterns using elementary numerical methods. In this paper, this new approach to modelling immiscible viscous fingering is tested by performing direct numerical simulations of previously published experimental water/oil displacements in 2D sandstone porous media. Experiments were modelled at adverse viscosity ratios ($${\mu }_{o}/{\mu }_{w}$$ μ o / μ w ), with oil viscosities ranging from μo = 412 to 7000 cP, i.e. for a viscosity ratio range, ($${\mu }_{o}/{\mu }_{w}$$ μ o / μ w ) $$\sim$$ ∼ 400–7000. These experiments have extensive production data as well as in situ 2D immiscible fingering images, measured by X-ray scanning. In all cases, very good quantitative agreement between experiment and modelling results is found, providing a strong validation of the new modelling approach. The underlying parameters used in the modelling of these unstable immiscible floods, the $${f}_{w}^{*}$$ f w ∗ functions, show very consistent and understandable variation with the viscosity ratio, ($${\mu }_{o}/{\mu }_{w}$$ μ o / μ w ).
Realistic immiscible viscous fingering, showing all of the complex finger structure observed in experiments, has proven to be very difficult to model using direct numerical simulation based on the two-phase flow equations in porous media. Recently, a method was proposed by the authors to solve the viscous-dominated immiscible fingering problem numerically. This method gave realistic complex immiscible fingering patterns and showed very good agreement with a set of viscous unstable 2D water → oil displacement experiments. In addition, the method also gave a very good prediction of the response of the system to tertiary polymer injection. In this paper, we extend our previous work by considering the effect of wettability/capillarity on immiscible viscous fingering, e.g. in a water → oil displacements where viscosity ratio $$\left( {\mu_{{\text{o}}} /\mu_{{\text{w}}} } \right) \gg 1$$ μ o / μ w ≫ 1 . We identify particular wetting states with the form of the corresponding capillary pressure used to simulate that system. It has long been known that the broad effect of capillarity is to act like a nonlinear diffusion term in the two-phase flow equations, denoted here as $$D(S_{w} )$$ D ( S w ) . Therefore, the addition of capillary pressure, $$P_{c} (S_{w} )$$ P c ( S w ) , into the equations acts as a damping or stabilisation term on viscous fingering, where it is the derivative of this quantity that is important, i.e. $$D(S_{w} )\sim\left( {dP_{c} (S_{w} )/dS_{w} } \right)$$ D ( S w ) ∼ d P c ( S w ) / d S w . If this capillary effect is sufficiently large, then we expect that the viscous fingering to be completely damped, and linear stability theory has supported this view. However, no convincing numerical simulations have been presented showing this effect clearly for systems of different wettability, due to the problem of simulating realistic immiscible fingering in the first place (i.e. for the viscous-dominated case where $$P_{c} = 0$$ P c = 0 ). Since we already have a good method for numerically generating complex realistic immiscible fingering for the $$P_{c} = 0$$ P c = 0 case, we are able for the first time to present a study examining both the viscous-dominated limit and the gradual change in the viscous/capillary force balance. This force balance also depends on the physical size of the system as well as on the length scale of the capillary damping. To address these issues, scaling theory is applied, using the classical approach of Rapport (1955), to study this scaling in a systematic manner. In this paper, we show that the effect of wettability/capillarity on immiscible viscous fingering is somewhat more complex and interesting than the (broadly correct) qualitative description above. From a “lab-scale” base case 2D water → oil displacement showing clear immiscible viscous fingering which we have already matched very well using our numerical method, we examine the effects of introducing either a water wet (WW) or an oil wet (OW) capillary pressure, of different “magnitudes”. The characteristics of these two cases (WW and OW) are important in how the value of corresponding $$D(S_{w} )$$ D ( S w ) functions, relate to the (Buckley–Leverett) shock front saturation, $$S_{wf}$$ S wf , of the viscous-dominated ($$P_{c} = 0$$ P c = 0 ) case. By analysing this, and carrying out some confirming calculations, we show clearly why we expect to see much clearer immiscible fingering at the lab scale in oil wet rather than in water wet systems. Indeed, we demonstrate why it is very difficult to see immiscible fingering in WW lab systems. From this finding, one might conclude that since no fingering is observed for the WW lab-scale case, then none would be expected at the larger “field” scale. However, by invoking scaling theory—specifically the viscous/capillary scaling group, $$C_{{{\text{VC1}}}}$$ C VC1 , (and a corresponding “shape group”, $$C_{{{\text{S}}1}}$$ C S 1 ), we demonstrate very clearly that, although the WW viscous fingers do not usually appear at the lab scale, they emerge very distinctly as we “inflate” the system in size in a systematic manner. In contrast, we demonstrate exactly why it is much more likely to observe viscous fingering for the OW (or weakly wetting) case at the lab scale. Finally, to confirm our analysis of the WW and OW immiscible fingering conclusions at the lab scale, we present two experiments in a lab-scale bead pack where $$\left( {\mu_{{\text{o}}} /\mu_{{\text{w}}} } \right) = 100$$ μ o / μ w = 100 ; no fingering is seen in the WW case, whereas clear developed immiscible fingering is observed in the OW case.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.