A Method to Correct Steady-State Relative Permeability Measurements for Inhomogeneous Saturation Profiles in One-Dimensional Flow

Traditionally, steady-state relative permeability is calculated from measurements on small rock samples using Darcy’s law and assuming a homogenous saturation profile and constant capillary pressure. However, these assumptions are rarely correct as local inhomogeneities exist; furthermore, the wetting phase tends to be retained at the outlet–the so-called capillary end effect. We have introduced a new method that corrects the relative permeabilities, analytically, for an inhomogeneous saturation profile along the flow direction. The only data required are the measured pressure drops for different fractional flow values, an estimate of capillary pressure, and the saturation profiles. An optimization routine is applied to find the range of relative permeability values consistent with the uncertainty in the measured pressure. Assuming a homogenous saturation profile systematically underestimates the relative permeability and this effect is most marked for media where one of the phases is strongly wetting with a noticeable capillary end effect. Relative permeabilities from seven two-phase flow experiments in centimetre-scale samples with different wettability were corrected while reconciling some hitherto apparently contradictory results. We reproduce relative permeabilities of water-wet Bentheimer sandstone that are closer to other measurements in the literature on larger samples than the original analysis. Furthermore, we find that the water relative permeability during waterflooding a carbonate sample with a wide range of pore sizes can be high, due to good connectivity through the microporosity. For mixed-wet media with lower capillary pressure and less variable saturation profiles, the corrections are less significant.


Introduction
Many natural and industrial processes in fields such as hydrology, petroleum engineering, carbon sequestration, and fuel cell design depend fundamentally on multiphase flow in porous media. Relative permeability and capillary pressure are the two critical parameters needed to quantify multiphase flow (Adler and Brenner 1988;Bear 1988). In recent years, pore-scale X-ray imaging has been combined with special core analysis to determine capillary pressure simultaneously with relative permeability (Gao et al. 2017;Lin et al. 2018;Zou et al. 2018;Zou and Armstrong 2019;Jackson et al. 2020). Small centimetresized cylindrical rock samples (cores) are typically used for flow experiments to allow the use of X-ray imaging at a high resolution of a few microns. Both the pore geometry and fluid-fluid interfaces can be visualized and thus pore-scale multiphase flow mechanisms can be quantified. However, the samples are often smaller than those used in traditional core flooding studies without pore-scale imaging.
Some discrepancies have been observed in relative permeability measurements between samples of different size (Qu et al. 2022) as well as between experimental measurements and pore-scale modelling studies (Raeini et al. 2022). Typically, lower relative permeabilities in the pore-scale imaging experiments are observed. The reason for this may be threefold. First, both experiments and predictions may have a large degree of uncertainty (Raeini et al. 2022). Second, local heterogeneity exists, leading to inhomogeneous saturation and the development of capillary pressure gradients. Third, in smaller samples with a diameter of a few millimetres, but lengths of several centimetres, the flow is largely onedimensional: this means that any inhomogeneity that locally reduces the saturation may significantly restrict flow.
Traditionally, the steady-state relative permeability is calculated using Darcy's law with experimental measurements of pressures and average saturation. This assumes that the saturation profile is homogeneous with the same pressure drop in each phase (the capillary pressure is constant). However, a high pressure gradient will build at a region where the saturation of one phase is low as it will restrict the flow of that phase. Although the average saturation for the whole core may be higher, the overall pressure drop may be dominated by the low saturation regions. As a result, a high pressure drop and a low relative permeability are then assigned to a high average saturation. In this case, the relative permeability is underestimated if a homogeneous saturation profile is assumed.
Furthermore, experimental measurements of relative permeability are affected by the capillary end effect-the capillary pressure vanishes at the outlet, causing retention of the wetting phase (Leverett 1941;Hassker et al. 1944;Osoba et al. 1951). One common method for reducing this effect is to use a flow rate or a high pressure drop in experiments, or to perform the experiment in a larger core (Osoba et al. 1951;Chen and Wood 2001;Gupta and Maloney 2016). However, it is not feasible to perform experiments with high flow rates in many cases, particularly if the aim is to maintain capillary-controlled conditions at the pore scale. To account for the capillary end effect, relative permeabilities are found consistent with the experimental data using numerical techniques, but additional information is often required, such as an independent measurement of capillary pressure as a function of saturation (Virnovsky et al. 1995;Huang and Honarpour 1998;Qadeer et al. 1988;Gupta and Maloney 2016;Zou et al. 2020). Capillary heterogeneity in rocks also leads to inhomogeneous saturation profiles (Krevor et al. 2011). Jackson et al. (2020 developed three-dimensional numerical models to predict the variations in trapping and relative permeability considering the capillary heterogeneity. However, in pore-scale imaging experiments, a long and thin core is typically used and thus the flow is approximately one-dimensional making a three-dimensional analysis unnecessary. In this study, we propose an analytical method to correct steady-state relative permeabilities for inhomogeneous saturation profiles including the capillary end effect. This method does not rely on an independent measurement of capillary pressure if capillary pressure can be estimated from interfacial curvature in the images, and does not rely on the use of complex history-matching techniques (Lai and Brandt 1988;Li et al. 2012;Wang et al. 2010;Zhang et al. 2012). To illustrate the approach, seven oil-water two phase flow experiments with porescale imaging in the literature with measured saturation profiles and capillary pressures are reanalysed and discussed.

Relative Permeability Correction Method for Inhomogeneous Saturation Profiles
In the experiments we study, the fluid flow is upwards, and we assume that the flow is approximately one-dimensional in a long and thin core (see Table 1). The multiphase Darcy's law for one-dimensional vertical two-phase flow is where q is the Darcy velocity, k r is the relative permeability, μ is the viscosity, P is the pressure, ρ is the density, and g is the gravitational acceleration. The subscripts refer to the two phases, water (phase 1) and a less dense phase 2 (gas or oil). K is the absolute permeability.
(1) and (2) between 0 (the inlet) and the sample length L (the outlet) we find the pressure drop, ΔP, across the sample, which we assume is measured in phase 2 (this will be the case in the experiments we analyse later): where we account for pressure differences between the phases due to differences in capillary pressure and buoyancy. The capillary pressure at the inlet can be estimated from curvature measurements on the fluid-fluid interfaces or from independent measurements, whereas the capillary pressure at the outlet is assumed to be zero since here both phases are produced with no pressure difference between them. We assume there is a unique relative permeability as a function of saturation: k r (S 1 ), which can be defined locally in sliceaverages of the image.
where we replace the integrals in Eqs. (6) and (7) using the linear sums. Δx is increment of sample length in the images, which is the voxel size in imaging experiments, j is the slice number of the images along the measured saturation profile, and n is the total number of slices. Here i labels the fractional flow: in the experiments we study there are N = 7 or 8 distinct fractional flow values, f 1 , for a waterflood experiment (see Appendix 1)-the first value is a water fractional flow of 0 , f 1 1 = 0 , and the last value is 1, f N 1 = 1. We determine the initial saturation of phase 1 (water in all the cases we study), S 1i , and the residual saturation of phase 2, S 2r , from the lowest slice-averaged saturations of phases 1 and 2 respectively measured in these experiments at f 1 1 and f N 1 respectively: we define We divide the measured saturation range from S 1i to 1-S 2r into the same number of increments as the number of fractional flows in each experiment, N. Hence we describe the relative permeabilities of phases 1 and 2 using seven or eight unknowns at discrete points k r1,i and k r2,i . The saturations for each relative permeability value are the average saturations at each fractional flow-that is the same saturations as for the uncorrected data. We use linear interpolation to define relative permeability values at any saturation between these discrete points as a weighted linear sum of the unknowns.
To determine the relative permeability values, we use differential evolution which is a stochastic optimization technique to minimize the errors between the pressure gradient integrated along the measured saturation profile and the experimentally measured pressure drop. Differential evaluation (DE) is a stochastic population-based optimization algorithm introduced by Storn (1996). This metaheuristic method can search a large space of candidate solutions. This method makes a few assumptions about problem being optimized and does not use the gradient of the problem being optimized which means does not require the optimization problem to be differentiable. The DE algorithm starts with a population of candidate solutions (called agents). These agents are moved around in the search-space using simple mathematical formulae to combine the existing agents from the population. If the new agent is an improvement, then it is accepted and forms part of the population, otherwise it is simply discarded. The process is repeated until a satisfactory solution is found. The objective function is: where k r1 and k r2 are the adjusted parameters. The solution is constrained to ensure that the relative permeabilities are monotonic: The monotonicity constraints (Eqs. 9 and 10) are implemented using a penalty function. The objective function (Eq. 8) was minimized to have two characteristics. The optimized relative permeabilities should give less error than the error based on relative permeabilities obtained from an assumed uniform saturation profile. Also, the optimized relative permeabilities should predict pressure differences within the range of uncertainty of the measurements, which we assume to be 20%. This uncertainty in the experimental measurement leads to an uncertainty in the optimization and all the solutions with the above characteristics will be accepted. We find a range of relative permeability consistent with the uncertainty in the experimental measurements.

Results and Discussion
Relative permeabilities from seven two-phase flow experiments obtained using porescale imaging from our research group with the results available in the literature were corrected for inhomogeneous saturation profiles using the method developed here. These seven two-phase flow experiments are the only publicly available datasets with measured pressure drops, saturation profile and capillary pressure combined. The pressure drops, saturation profiles along the length of the samples, and capillary pressures for each fractional flow are available in these experiments, see Appendix 1. These seven datasets include water-wet and mixed-wet Bentheimer sandstones, water-wet and oil-wet Estaillades carbonates, a mixed-wet reservoir carbonate, and a reservoir sandstone with the same geometry under two wettability conditions (Table 1). In these experiments phase 1 is water (or brine) and phase 2 is oil.

Bentheimer Sandstone
Bentheimer sandstone comes from shallow marine formation deposited during the Lower Cretaceous, which forms a reservoir for oil found on the border between the Netherlands and Germany (Peksa et al. 2015). Bentheimer is considered as an ideal rock for reservoir studies because of its lateral continuity and homogeneous block-scale nature. Bentheimer has been used as a benchmark for multiphase flow experiments (Øren et al. 1998;Herring et al. 2013;Alizadeh and Piri 2014;Gao et al. 2017;Lin et al. 2018Lin et al. , 2019. Gao et al. (2017) and Lin et al. (2018) conducted steady-state two-phase flow experiments with coinjection of water and oil into water-wet Bentheimer samples. Lin et al. (2019) also conducted steady-state two-phase flow experiments on a mixed-wet Bentheimer sample. The measured saturation profiles and relative permeabilities before and after correction for the measured inhomogeneous saturation profiles can be found in Figs. 1 and 2.
The originally reported relative permeabilities are lower than those of the previously published for larger, full-core experiments under water-wet conditions (Øren et al. 1998;Alizadeh and Piri 2014). Lin et al. (2018) corrected the water relative permeabilities by considering a constant gradient in capillary pressure which did lead to an increase in the estimated water relative permeability, but it was still lower than in the other experiments. Qu et al. (2022) measured the rate of spontaneous imbibition on samples from the same block of stone used by Lin et al. (2018). Using a semi-analytical solution for the imbibition rate, they estimated a water relative permeability that was consistent with previous measurements, but higher than reported by Lin et al. (2018). This discrepancy in measurements on a homogeneous sandstone, although just within the reported experimental uncertainty, was concerning and raised the possibility that measurements of macroscopic properties on the small samples used for pore-scale imaging might not be as accurate as those on larger cores.
In Fig. 1 it is evident that in the water-wet case, not only there is a noticeable gradient in saturation but also the profiles themselves show large fluctuations. Where the water saturation is locally small, flow is restricted and the pressure gradient is large: associating a relative permeability, consistent with the measured pressure drop, to the average saturation will tend to underestimate the true value. For the mixed-wet sample, although small fluctuations are also present, they are less significant, while the capillary pressure is very low, and so little correction in relative permeability is needed. This is confirmed in Fig. 2, where the water (phase 1) relative permeabilities in water-wet Bentheimer systematically increase after our correction. The end point values are close to those of the other experiments in the literature: the apparent discrepancy is resolved through a proper accounting for variations in the saturation profile. Relative permeabilities for phase 1 (water/brine, blue) and phase 2 (oil or gas, red) with and without the inhomogeneous saturation profile correction shown on a linear axis (left) and on a semi-logarithmic axis (right): a and b water-wet Bentheimer sandstone (Lin et al. 2018); c and d water-wet Bentheimer sandstone (Lin et al. 2018) compared to the pore-scale imaging experimental data by Gao et al. (2017), and traditional core-flood data by Øren et al. (1998) and Alizadeh and Piri (2014), and from spontaneous imbibition measurements by Qu et al. (2022); e and f mixed-wet Bentheimer sandstone (Lin et al. 2019). The error bars indicate the uncertainty in the measurement of saturation and pressure drop, while the shaded areas represent the uncertainty of the corrected relative permeabilities after optimization For the mixed-wet Bentheimer, the relative permeabilities show little change after correction. This is reassuring, and adds confidence to the results, which indicate potentially favourable waterflood recovery in this case.

Estaillades and Reservoir Carbonates
Three carbonates, namely water-wet Estaillades, oil-wet Estaillades and a mixed-wet reservoir carbonate are selected in our study. Estaillades carbonate is a bioclastic limestone, which contains 99% calcite, quarried near the town of the Oppède, south of France. It has a complex sub-micron pore space with a large range of pore sizes, which is a benchmark microporous carbonate for relative permeability analysis. The reservoir sample comes from a producing oil field in the Middle East; wettability was restored by placing the rock in contact with crude oil from the same reservoir for several weeks at high temperature and pressure (Alhammadi et al. 2020). Figure 3 shows the measured saturation profiles in the experiments while Fig. 4 compares the relative permeability results with and without correction.
Both oil and water relative permeability increase significantly after correction as there is a clear saturation gradient in the water-wet Estaillades carbonate experiments with significant local fluctuations, Fig. 3a. Estaillades carbonate has a large fraction of  (Lin et al. 2021); c mixed-wet reservoir carbonate (Alhammadi et al. 2020) macropores that are connected through sub-resolution microporosity with a high initial water saturation at the beginning of waterflooding. Compared to the Bentheimer sandstone results shown in Fig. 2, although there is trapping of oil in macropores, subresolution microporosity provides good connectivity leading to a much higher end-point water relative permeability even in this water-wet sample. Relative permeabilities for water (brine, phase 1, blue) and oil (phase 2, red) with and without the inhomogeneous saturation profile correction shown on a linear axis (left) and on a semi-logarithmic axis (right): a and b water-wet Estaillades carbonate (Gao et al. 2019); c and d oil-wet Estaillades carbonate (Lin et al. 2021); e and f mixed-wet reservoir carbonate (Alhammadi et al. 2020). The error bars indicate the uncertainty of the measurements of saturation and pressure drop, while the shaded areas represent the uncertainty of the corrected oil and water relative permeability after optimization In oil-wet Estaillades carbonate, the saturation profiles are smoother with no obvious gradient: the oil relative permeabilities increase whereas there are no obvious changes in the water relative permeabilities after correction.
The reservoir sample shows a relatively smooth saturation profile (Fig. 3c) and the correction only makes small adjustments to the oil and water relative permeabilities. This confirms the validity of the original analysis (Alhammadi et al. 2020) and the suggestion that mixed-wet media can lead to high recoveries in waterflooding. Gao et al. (2020) simultaneously measured relative permeability and capillary pressure on the same reservoir sandstone sample under water-wet and mixed-wet conditions, before and after exposing the rock to a crude oil at high temperature. Figure 5 shows the measured saturation profiles for each fractional flow for both experiments, while Fig. 6 compares the relative permeability results before and after correction. Raeini et al. (2022) used a generalized pore network model to predict capillary pressure and relative permeability for this case. The modelling results showed a reasonable match with the measurements within the uncertainty of the experiments. However, the presented experimental water relative permeabilities were lower than the predictions. Here, we see an increase in the relative permeability values after correction, which are then closer to the modelling results (Raeini et al. 2022). Here the increase in relative permeability values is more obvious in the mixed-wet experiment, as more heterogeneous saturation profiles are observed compared to the water-wet case.

Conclusions
We have introduced a novel method to correct for inhomogeneous saturation profiles and the so-called capillary end effect in steady-state measurements of relative permeability. This method only requires pressure drops, saturation profiles and capillary pressures for each fractional flow. It assumes one-dimensional flow along the sample and that the relative permeability can be written as a single function of saturation. The method does not require additional measurements from independent techniques if the capillary pressure can be estimated from the interfacial curvature in the pore-space images. Where there are significant variations in the saturation profile, our methodology leads to larger estimates of relative permeability than would be assumed assigning a value to the average saturation. Relative permeabilities for oil (red) and water (blue) with and without the inhomogeneous saturation profile correction shown on a linear axis (left) and on a semi-logarithmic axis (right): a and b water-wet reservoir sandstone ; c and d mixed-wet reservoir sandstone . The error bars indicate the uncertainty of the measured saturation and pressure drop, while the shaded areas represent the uncertainty of the corrected oil and water relative permeability after optimization We applied this methodology to seven experimental datasets in the literature. We resolved an apparent discrepancy in results for water-wet Bentheimer sandstone between pore-scale imaging experiments and the results of traditional core floods, as well as measurements of the rate of spontaneous imbibition. For a carbonate sample with a wide range of pore size, we demonstrated that the water relative permeability can be high, even under water-wet conditions, due to good connectivity through microporosity. For two of the mixed-wet media studied, the measured saturation profile showed only small fluctuations and the capillary pressure was low, and hence the correction to relative permeability was insignificant. This strengthens the original conclusion that mixed-wet conditions may lead to favourable local displacement efficiency in waterflooding.
Overall, we suggest applying this correction to experiments where the flow is largely one-dimensional and where the saturation profile is measured.

Appendix 1: Saturation, pressure, and relative permeability with and without correction in oil-water two phase flow in seven experiments
See Tables 2, 3, 4, 5, 6, 7 and 8.  Author Contributions All authors contributed to the study conception and design, investigation, original draft writing, edit and review. The software development was performed by SF. The funding and supervision of the work were conducted by MJB and BB.

Funding
The authors have not disclosed any funding.

Code Availability
The code used for the analysis in this paper can be found at the link: https:// github. com/ Imper ialCo llege London/ relPe rmCor recti on

Declarations
Conflict of interest The author declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.