Dynamic Capillary Pressure Curves From Pore-Scale Modeling in Mixed-Wet-Rock Images
- Yingfang Zhou (International Research Institute of Stavanger) | Johan Helland (International Research Institute of Stavanger) | Espen Jettestuen (International Research Institute of Stavanger)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- March 2013
- Document Type
- Journal Paper
- 634 - 645
- 2013. Society of Petroleum Engineers
- 5.3.2 Multiphase Flow, 1.6.9 Coring, Fishing, 4.3.4 Scale
- 9 in the last 30 days
- 509 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
In reservoir multiphase-flow processes with high flow rates, both viscous and capillary forces determine the pore-scale fluid configurations, and significant dynamic effects could appear in the capillary pressure/saturation relation. We simulate dynamic and quasistatic capillary pressure curves for drainage and imbibition directly in scanning-electron-microscope (SEM) images of Bentheim sandstone at mixed-wet conditions by treating the identified pore spaces as tube cross sections. The phase pressures vary with length positions along the tube length but remain unique in each cross section, which leads to a nonlinear system of equations that are solved for interface positions as a function of time. The cross-sectional fluid configurations are computed accurately at any capillary pressure and wetting condition by combining free-energy minimization with a menisci-determining procedure that identifies the intersections of two circles moving in opposite directions along the pore boundary. Circle rotation at pinned contact lines accounts for mixed-wet conditions. Dynamic capillary pressure is calculated with volume-averaged phase pressures, and dynamic capillary coefficients are obtained by computing the time derivative of saturation and the difference between the dynamic and static capillary pressure. Consistent with previously reported measurements, our results demonstrate that, for a given water saturation, simulated dynamic capillary pressure curves are at a higher capillary level than the static capillary pressure during drainage, but at a lower level during imbibition, regardless of the wetting state of the porous medium. The difference between dynamic and static capillary pressure becomes larger as the pressure step applied in the simulations is increased. The model predicts that the dynamic capillary coefficient is a function of saturation and is independent of the incremental pressure step, which is consistent with results reported in previous studies. The dynamic capillary coefficient increases with decreasing water saturation at water-wet conditions, whereas, for mixed- to oil-wet conditions, it increases with increasing water saturation. The imbibition simulations performed at mixed- to oil-wet conditions also indicate that the dynamic capillary coefficient increases with decreasing initial water saturation.
The proposed modeling procedure provides insights into the extent of dynamic effects in capillary pressure curves for realistic mixed-wet pore spaces, which could contribute to the improved interpretation of core-scale experiments. The simulated capillary pressure curves obtained with the pore-scale model could also be applied in reservoir-simulation models to assess dynamic pore-scale effects on the Darcy scale.
|File Size||1 MB||Number of Pages||12|
Barenblatt, G.I. and Vinnitchenko, A.P. 1980. Nonequilibrium Seepage ofImmiscible Fluids. Adv. Mech. 3: 35-50.
Bottero, S., Hassanizadeh, S.M., Kleingeld, P.J. et al. 2006. ExperimentalStudy of Dynamic Capillary Pressure Effect in Two-Phase Flow on Porous Media.Proceedings of the XVI International Conference on ComputationalMethods in Water Resources, Copenhagen, Denmark.
Camps-Roach, G., O'Carroll, D.M., Newson, T.A. et al. 2010. ExperimentalInvestigation of Dynamic Effects in Capillary Pressure: Grain Size Dependencyand Upscaling. Water Resour. Res. 46: W08544. http://dx.doi.org/10.1029/2009WR008881.
Chen, S. and Doolen, G.D. 1998. Lattice Boltzmann Methods for Fluid Flows.Annu. Rev. Fluid Mech. 30: 329-364.
Dahle, H.K. and Celia, M.A. 1999. A Dynamic Network Model for Two-PhaseImmiscible Flow. Computat. Geosci. 3 (1): 1-22. http://dx.doi.org/10.1023/A:1011522808132.
Dahle, H.K., Celia, M.A., and Hassanizadeh S.M. 2005. Bundle-of-Tubes Modelfor Calculating Dynamic Effects in the Capillary-Pressure/SaturationRelationship. Transp. Porous Media 58 (1): 5-22. http://dx.doi.org/10.1007/s11242-004-5466-4.
Dong, M., Dullien, F.A.L., Dai, L. et al. 2005. Immiscible Displacement inthe Interacting Capillary Bundle Model Part I. Development of InteractingCapillary Bundle Model. Transp. Porous Media 59 (1): 1-18.http://dx.doi.org/10.1007/s11242-004-0763-5.
Frette, O.I. and Helland, J.O. 2010. A Semi-Analytical Model for Computationof Capillary Entry Pressures and Fluid Configurations in Uniformly Wet PoreSpaces From 2D Rock Images. Adv. Water Resour. 33 (8):846-866. http://dx.doi.org/10.16/j.advwatres.2010.05.002.
Hassanizadeh, S.M., Celia, M.A. and Dahle, H.K. 2002. Dynamic Effect in theCapillary Pressure/Saturation Relationship and Its Impacts on Unsaturated Flow.Vadose Zone J. 1: 38-57. http://dx.doi.org/10.2113/1.1.38.
Hassanizadeh, S.M. and Gray, W.G. 1993. Thermodynamic Basis of CapillaryPressure in Porous Media. Water Resour. Res. 29 (10):3389-3405. http://dx.doi.org/10.1029/93WR01495.
He, X., Doolen, G.D., and Clark, T. 2002. Comparison of the LatticeBoltzmann Method and the Artificial Compressibility Method for Navier-StokesEquations. J. Comput. Phys. 179: 439-451.
Helland, J.O. and Frette, O.I. 2010. Computation of Fluid Configurations andCapillary Pressures in Mixed-Wet 2D Pore Spaces From Rock Images.Proceedings of the XVIII International Conference on WaterResources, Barcelona, Spain.
Helland, J.O. and Skjæveland, S.M. 2006. Physically Based Capillary PressureCorrelation for Mixed-Wet Reservoirs From a Bundle-of-Tubes Model. SPEJ. 11 (2): 171-180. http://dx.doi.org/10.2118/89428-PA.
Huang, H., Thorn Jr., D.T., Schaap, M.G. et al. 2007. Proposed Approximationfor Contact Angles in Shan-and-Chen-Type Multicomponent Multiphase LatticeBoltzmann Models. Phys. Rev. E 7 6 (6): 066701. http://dx.doi.org/10.1103/PhysRevE.76.066701.
Hughes, R.G. and Blunt, M.J. 2000. Pore-Scale Modeling of Rate Effects inImbibition. Transp. Porous Media 40 (3): 295-322. http://dx.doi.org/10.1023/A:1006629019153.
Jadhunandan, P.P. and Morrow, N.R. 1995. Effect of Wettability on WaterfloodRecovery for Crude-Oil/Brine/Rock Systems. SPE Res Eng 10(1): 40-46. http://dx.doi.org/10.2118/22597-PA.
Joekar-Niasar, V. and Hassanizadeh, S.M. 2011. Effect of Fluids Propertieson Nonequilibrium Capillarity Effects: Dynamic Pore-Network Modeling. Int.J. Multiphase Flow 37 (2): 198-214.
Kalaydjian, F.J-M. 1992. Dynamic Capillary Pressure Curve for Water/OilDisplacement in Porous Media: Theory Vs. Experiment. Paper SPE 24813 presentedat the SPE Annual Technical Conference and Exhibition, Washington, DC, 4-7October. http://dx.doi.org/10.2118/24813-MS.
Korteland, S., Bottero, S., Hassanizadeh, S.M. et al. 2010. What Is theCorrect Definition of Average Pressure? Transp. Porous Media 84: 153-175. http://dx.doi.org/10.1007/s11242-009-9490-2.
Lago, M. and Araujo, M. 2001. Threshold Pressure in Capillaries WithPolygonal Cross Section. J. Coll. Interface Sci. 243 (1):219-226. http://dx.doi.org/10.1006/jcis.2001.7872.
Ma, S., Mason, G., and Morrow, N.R. 1996. Effect of Contact Angle onDrainage and Imbibition in Regular Polygonal Tubes. Coll. Surf. A Phys. Eng.Asp. 117 (3): 273-291. http://dx.doi.org/10.1016/0927-7757(96)03702-8.
Mirzaei, M. and Das, D.B. 2007. Dynamic Effects in CapillaryPressure-Saturations Relationships for Two-Phase Flow in 3D Porous Media:Implications of Micro-Heterogeneities. Chem. Eng. Sci. 62(7): 1927-1947. http://dx.doi.org/10.1016/j.ces.2006.12.039.
Nourgaliev, R.R., Dinh, T.N., Theofanous, T.G. et al. 2003. The LatticeBoltzmann Equation Method: Theoretical Interpretation, Numerics, andImplications. Intl. J. Multiphase Flows 29 (1): 117-169. http://dx.doi.org/10.1016/S0301-9322(02)00108-8.
O'Carroll, D.M., Mumford, K.G., Abriola, L.M. et al. 2010. Influence ofWettability Variations on Dynamic Effects in Capillary Pressure. WaterResour. Res. 46 (8): W08505. http://dx.doi.org/10.1029/2009WR008712.
Ruth, D. and Bartley, J. 2002. A Perfect Cross-Flow Model for Two-Phase Flowin Porous Media. Paper SCA2002-05 presented at the International Symposium ofSociety of Core Analysts, Monterey, California.
Schembre, J.M. and Kovscek, A.R. 2006. Estimation of Dynamic RelativePermeability and Capillary Pressure From Countercurrent Imbibition Experiments.Transp. Porous Media 65 (1): 31-51. http://dx.doi.org/10.1007/s11242-005-6092-5.
Shan, X. and Chen, H. 1993. Lattice Boltzmann Model for Simulating FlowsWith Multiple Phases and Components. Phys. Rev. E. 47 (3):1815-1819. http://dx.doi.org/10.1103/PhysRevE.47.1815.
Silin, D. and Patzek, T. 2002. On Barenblatt's Model of SpontaneousCountercurrent Imbibition. Transp. Porous Media 54 (3):297-322. http://dx.doi.org/10.1023/B:TIPM.0000003678.85526.b1.
Succi, S. 2001. The Lattice Boltzmann Equation for Fluid Dynamics andBeyond. Oxford, England: Oxford University Press.
Unsal, E., Mason, G., Morrow, N.R. et al. 2007a. Co-Current andCounter-Current Imbibition in Independent Tubes of NonAxisymmetric Geometry.J. Coll. Interface Sci. 306 (1): 105-117. http://dx.doi.org/10.16/j.jcis.2006.10.042.
Unsal, E., Mason, G., Ruth, D.W. et al. 2007b. Co- and Counter-CurrentSpontaneous Imbibition Into Groups of Capillary Tubes With Lateral ConnectionsPermitting Cross-Flow. J. Coll. Interface Sci. 315 (1):200-209. http://dx.doi.org/j.jcis.2007.05.070.
Wang, J. and Dong, M. 2011. Trapping of the NonWetting Phase in anInteracting Triangular Tube-Bundle Model. Chem. Eng. Sci. 66 (3): 250-259. http://dx.doi.org/10.16/j.ces.2010.10.009.
Zhou, Y., Helland, J.O., and Hatzignatiou, D.G. 2011. A Model for Imbibitionin Pore Spaces From 2D Rock Images. Proceedings of the InternationalConference on Flows and Mechanics in Natural Porous Media from Pore to FieldScale. Rueil-Malmaison, France: IFP Energies Nouvelles, 16-18 November.
Zhou, Y., Helland, J.O., and Hatzignatiou, D.G. 2012a. Pore-Scale Modellingof Water Flooding in Mixed-Wet Rock Images: Effects of Initial Saturation andWettability. Paper SPE 154284 presented at the SPE Improved Oil RecoverySymposium, Tulsa, Oklahoma, 14-18 April. http://dx.doi.org/10.2118/154284-MS.
Zhou, Y., Helland, J., and Hatzignatiou, D.G. 2012b. A DimensionlessCapillary Pressure Function for Imbibition Derived From Pore-Scale Modeling inMixed-Wet-Rock Images. SPE J. SPE-154129-PA (in press; posted 28December 2012). http://dx.doi.org/10.2118/154129-PA.