Computing Relative Permeability and Capillary Pressure of Heterogeneous Rocks Using Realistic Boundary Conditions

Relative permeability and capillary pressure are key parameters in multiphase flow modelling. In heterogeneous porous media, flow direction- and flow-rate dependence result from non-uniform saturation distributions that vary with the balance between viscous, gravitational, and capillary forces. Typically, relative permeability is measured using constant inlet fractional-flow—constant outlet fluid pressure conditions on samples mounted between permeable porous plates to avoid capillary end-effects. This setup is replicated in numeric experiments but ignores the extended geologic context beyond the sample size, impacting the saturation distribution and, consequently, the upscaled parameters. Here, we introduce a new workflow for measuring effective relative permeability and capillary pressure at the bedform scale while considering heterogeneities at the lamina scale. We harness the flexibility of numeric modelling to simulate continuum-REV-scale saturation distributions in heterogeneous rocks eliminating boundary artefacts. Periodic fluid flux boundary conditions are applied in combination with arbitrarily oriented, variable-strength pressure gradient fields. The approach is illustrated on a periodic model of cross-bedded sandstone. Stepping saturation while applying variable-strength pressure-gradient fields with different orientations, we cover the capillary-viscous force balance spectrum of interest. The obtained relative permeability and capillary pressure curves differ from ones obtained with traditional approaches highlighting that the definition of force balances needs consideration of flow direction as an additional degree of freedom. In addition, we discuss when the common viscous and the capillary limits are applicable and how they vary with flow direction in the presence of capillary interfaces.