1887

Abstract

Summary

Water-alternate-gas injections for improved recovery, slug-wise CO2 injections for permanent storage, and seasonal storage of natural gas and hydrogen are all processes that lead to cyclic flows of drainage and imbibition in the reservoir. Reservoir simulations of these processes need flow functions with hysteresis, which is typically implemented using correlations and hysteresis loop logic that could be inaccurate, especially for the description of higher-order scanning curves. Typically, measuring enough hysteresis-loop data from core-scale experiments is not feasible, and pore-scale simulation of these relations directly on micro-CT images is computationally demanding.

The discrete-domain model (DDM) represents an efficient, physics-based, method to describe hysteresis. DDM divides the porous rock into a set of compartments where each compartment is described by a Helmholtz energy and a local saturation. These energy functions describe a rugged energy landscape where hysteresis occurs due to irreversible jumps across energy barriers separating the local energy minima. For pressure-controlled displacement, the evolution equations of the local saturations are uncoupled, while for saturation-controlled (or, rate-controlled) displacement they are coupled by global saturation constraints where capillary pressure is a Lagrange multiplier. The global coupling leads to fluid redistribution among compartments (cooperative behavior) and pressure fluctuations.

Thus far, the DDM has only employed simple phenomenological energy functions. The objective of this work is to explore the applicability of the DDM on realistic data from rock samples. For that purpose, we simulate primary drainage, imbibition, and subsequent scanning loops, on segmented micro-CT images of Castlegate sandstone, using a level set model (LSM) for capillary-controlled displacement. The generated data is used to calculate the energy functions for different compartment architectures in the DDM. From the data we also explore how the interfacial area leads to differences in the energy between drainage and imbibition.

We find that the DDM reproduces the capillary pressure curves from both the pressure- and saturation-controlled LSM simulations using the energy functions from the saturation-controlled case. A finer compartment division of the rock sample leads to more energy minima and smoother results. Using the same energy landscape for either drainage or imbibition on all processes (including scanning curves) leads to a slight deviation from the LSM results, whereas the case with different energy landscapes for drainage and imbibition shows excellent agreement. Hence, the DDM emerges as a suitable pore-to-core upscaling approach for hysteresis as its compartmental description is based on extensive properties.

© EAGE Publications BV
Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.202531055
2025-04-02
2025-12-10

  • From This Site

    /content/papers/10.3997/2214-4609.202531055
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. Andrew, M., Bijeljic, B. and Blunt, M.J. [2014] Pore-by-pore capillary pressure measurements using X-ray microtomography at reservoir conditions: Curvature, snap-off, and remobilization of residual CO2. Water Resour. Res., 50, 8760–8774.
     [Google Scholar]
  2. Andrew, M., Bijeljic, B. and Blunt, M.J. [2015] The imaging of dynamic multiphase fluid flow using synchrotron-based x-ray microtomography at reservoir conditions. Transp. Porous Med., 110, 1–24.
     [Google Scholar]
  3. Berg, S., Ott, H., Klopp, S.A., Schwing, A., Neiteler, R., Brussee, N., Makurat, A., Leu, L., Enzmann, F, Schwartz, J.O., Kersten, M., Irvine, S. and Stampanoni, M. [2013] Real-time 3D imaging of Haines jumps in porous media flow. Proc. Natl. Acad. Sci. U. S. A., 110(10), 3755–3759. Doi: 10.1073/pnas.l2213731.0.
     https://doi.org/10.1073/pnas.l2213731.0 [Google Scholar]
  4. Blunt, M. [2000] An empirical model for three-phase relative permeability. SPE Journal, 5, 435–445.
     [Google Scholar]
  5. Brooks, R.H. and Corey, A.T. [1964] Hydraulic properties of porous media. Hydraulic paper no. 3, Colorado State University.
     [Google Scholar]
  6. Carlson, FM. [1981] Simulation of relative permeability hysteresis to the nonwetting phase, presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, October 1981, SPE-10157-MS.
     [Google Scholar]
  7. Cueto-Felgueroso, L. and Juanes, R. [2016] A discrete-domain description of multiphase flow in porous media: Rugged energy landscapes and the origin of hysteresis. Geophys. Res. Lett., 43, 1615–1622.
     [Google Scholar]
  8. Ebadi, M., McClure, J., Mostaghimi, P. and Armstrong, R.T [2024] Investigating an extended multiphase flow model that includes specific interfacial area. Computer Methods in Applied Mechanics and Engineering, 418, 116594.
     [Google Scholar]
  9. Gao, Y., Lin, Q., Bijeljic, B. and Blunt, M.J. [2017] X-ray Microtomography of Intermittency in Multiphase Flow at Steady State Using a Differential Imaging Method. Water Resources Research, 53(12), 10274–10292.
     [Google Scholar]
  10. van Genuchten, M.T [1980] A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J., 44, 892–898.
     [Google Scholar]
  11. Hassanizadeh, S.M. and Gray, W.G. [1993] Thermodynamic basis of capillary pressure in porous media. Water Resources Research, 29(10), 3389–3405.
     [Google Scholar]
  12. Helland, J.O., Friis, H.A., Jettestuen, E. and Skjaeveland, S.M. [2017] Footprints of spontaneous fluid redistribution on capillary pressure in porous rock. Geophys. Res. Lett., 44, 4933–4943.
     [Google Scholar]
  13. Helland, J.O., Jettestuen, E. and Friis, H.A. [2021] A discrete-domain approach to three-phase hysteresis. Water Resources Research, 57, e2021WR029560.
     [Google Scholar]
  14. Hilfer, R. and Øren, PE. [1996] Dimensional analysis of pore scale and field scale immiscible displacement. Transp. Porous Med., 22, 53–72.
     [Google Scholar]
  15. Jerauld, G. and Salter, S. [1990] The effect of pore-structure on hysteresis in relative permeability and capillary pressure: pore-level modelling. Transport in Porous Media, 5, 103–151.
     [Google Scholar]
  16. Jettestuen, E., Aursjø, O., Helland, J.O., Vinningland, J.L. and Hiorth, A. [2024] Workflow for Direct Pore-Scale Simulation of Relative Permeability and Capillary Pressure Curves with Hysteresis at Low Capillary Numbers, presented at the SPE Norway Subsurface Conference, Bergen, Norway, April 2024, SPE-218427-MS.
     [Google Scholar]
  17. Jettestuen, E., Friis, H.A. and Helland, J.O. [2021] A locally conservative multiphase level set method for capillary-controlled displacements in porous media. J. Comput. Phys., 428, 109965.
     [Google Scholar]
  18. Jettestuen, E., Helland, J.O. and Prodanovic, M. [2013] A level set method for simulating capillary-controlled displacements at the pore scale with nonzero contact angles. Water Resour. Res., 49(8), 4645–4661. Doi: 10.1002/wrcr.203.4.
     https://doi.org/10.1002/wrcr.203.4 [Google Scholar]
  19. Juanes, R., Spiteri, E.J., OrrJr., F.M. and Blunt, M.J. [2006] Impact of relative permeability hysteresis on geological CO2 storage. Water Resources Research, 42(12), W12418.
     [Google Scholar]
  20. Killough, J.E. [1976] Reservoir simulation with history-dependent saturation functions. SPE Journal, 16(01), 37–48.
     [Google Scholar]
  21. Land, C.S. [1968] Calculation of imbibition relative permeability for two- and three-phase flow from rock properties. SPE Journal, 8, 149–156.
     [Google Scholar]
  22. Lomeland, F. and Ebeltoft, E. [2013] Versatile three-phase correlations for relative permeability and capillary pressure. SCA2013-034. paper presented at the International Symposium of the Society of Core Analysts held in Napa Valley, California, USA, 16–19 September, 14.
     [Google Scholar]
  23. Long, J., Hyder, M.N., Huang, R.Y.M. and Chen, P. [2005] Thermodynamic modeling of contact angles on rough heterogeneous surfaces. Adv. Int. Sci, 118, 173–190.
     [Google Scholar]
  24. Lysyy, M., Føyen, T, Birkeland Johannesen, E., Fernø, M. and Ersland, G. [2022] Hydrogen Relative Permeability Hysteresis in Underground Storage. Geophysical Research Letters, 49(17), e2022GL100364.
     [Google Scholar]
  25. Masalmeh, S. [2003] The effect of wettability heterogeneity on capillary pressure and relative permeability. J. Petrol. Sci. Eng., 39(3–4), 399–408.
     [Google Scholar]
  26. McClure, J.E., Armstrong, R.T., Berrill, M.A., Schluter, S., Berg, S., Gray, W.G. and Miller, C.T. [2018] Geometric state function for two-fluid flow in porous media. Phys Ref Fluids, 3, 084306.
     [Google Scholar]
  27. Morrow, N.R. [1975] The effects of surface roughness on contact angle with special reference to petroleum recovery. J. Can. Pet. Tech., 14(4), 42–53.
     [Google Scholar]
  28. Rücker, M., Berg, S., Armstrong, R.T., Georgiadis, A., Ott, H., Schwing, A., Neiteler, R., Brussee, N., Makurat, A., Leu, L., Wolf, M., Khan, E, Enzmann, F. and Kersten, M. [2015] From connected pathway flow to ganglion dynamics. Geophys. Res. Lett., 42, 3888–3894.
     [Google Scholar]
  29. Sheppard, A. and Prodanovic, M. [2015] Network Generation Comparison Forum, https://www.digitalrocksportal.org/projects/16.
     [Google Scholar]
  30. Skjæveland, S., Siqveland, L., Kjosavik, A., Hammervold Thomas, W. and Virnovsky, G. [2000] Capillary pressure correlation for mixed-wet reservoirs. SPE Reservoir Evaluation & Engineering, 3(1), 60–67.
     [Google Scholar]
  31. Spiteri, E.J. and Juanes, R. [2006] Impact of relative permeability hysteresis on the numerical simulation of WAG injection. J. Pet. Sci. Eng., 50, 115–139.
     [Google Scholar]
  32. Spiteri, E.J., Juanes, R., Blunt, M.J. and OrrJr., F.M. [2008] A new model of trapping and relative permeability hysteresis for all wettability characteristics. SPE Journal, 13, 277–288.
     [Google Scholar]
  33. The MathWorks Inc. [2022] MATLAB version: 9.13.0 (R2022b).
     [Google Scholar]
/content/papers/10.3997/2214-4609.202531055
Loading
Discrete-Domain Model for Hysteresis and its Application on Porous Rocks
Conference Proceedings 2025, 1 (2025); https://doi.org/10.3997/2214-4609.202531055
/content/papers/10.3997/2214-4609.202531055
/content/papers/10.3997/2214-4609.202531055
Loading

Data & Media loading...

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error