Next Article in Journal
Invert Drilling Fluids with High Internal Phase Content
Next Article in Special Issue
Research on Fresh and Hardened Sealing Slurries with the Addition of Magnesium Regarding Thermal Conductivity for Energy Piles and Borehole Heat Exchangers
Previous Article in Journal
Method for Sizing of a PV System for Family Home Using Economic Indicators
Previous Article in Special Issue
Single Well Productivity Prediction Model for Fracture-Vuggy Reservoir Based on Selected Seismic Attributes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on the Correction Method of the Capillary End Effect of the Relative Permeability Curve of the Steady State

1
College of Energy Resources, China University of Geosciences, Beijing 100083, China
2
State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Southwest Petroleum University), Chengdu 610500, China
3
Reservoir Department, Tianjin Dagang Oilfield Group Engineering Construction Co., Ltd., Binhai 300280, China
4
Oil and Gas Field Development Department, Shengli Oil Field Exploration and Development Research Institute, Dongying 257001, China
5
Reservoir Department, Shandong Province Zhong Di Yi Cai Petroleum Technology Co., Ltd., Dongying 257001, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(15), 4528; https://doi.org/10.3390/en14154528
Submission received: 17 May 2021 / Revised: 1 July 2021 / Accepted: 13 July 2021 / Published: 27 July 2021
(This article belongs to the Special Issue Advancements in Thermal and Energy Geotechnics)

Abstract

:
Relative permeability curve is a key factor in describing the characteristics of multiphase flow in porous media. The steady-state method is an effective method to measure the relative permeability curve of oil and water. The capillary discontinuity at the end of the samples will cause the capillary end effect. The capillary end effect (CEE) affects the flow and retention of the fluid. If the experimental design and data interpretation fail to eliminate the impact of capillary end effects, the relative permeability curve may be wrong. This paper proposes a new stability factor method, which can quickly and accurately correct the relative permeability measured by the steady-state method. This method requires two steady-state experiments at the same proportion of injected liquid (wetting phase and non-wetting phase), and two groups of flow rates and pressure drop data are obtained. The pressure drop is corrected according to the new relationship between the pressure drop and the core length. This new relationship is summarized as a stability factor. Then the true relative permeability curve that is not affected by the capillary end effect can be obtained. The validity of the proposed method is verified against a wide range of experimental results. The results emphasize that the proposed method is effective, reliable, and accurate. The operation steps of the proposed method are simple and easy to apply.

Graphical Abstract

1. Introduction

With the vigorous progress of tight oil and gas exploration and development, experts and scholars have continuously improved their understanding of tight reservoirs [1,2,3,4]. It is generally believed that tight reservoirs have small pores and throats, low permeability, and strong capillary forces [5,6]. The study of the mechanism of the fluids in porous media and influencing factors of tight reservoirs has important guiding significance for improving the recovery of this type of reservoir [7,8]. The relative permeability curve is the basic parameter for studying two-phase seepage in porous media and is an indispensable important data for tight oilfield development dynamic analysis, numerical simulation, and production forecasting [9,10,11,12,13,14,15].
The relative permeability data is usually obtained from the core displacement experiment of the tight reservoir in the laboratory. In the experiment, the end of the core is in contact with the outside which could be viewed as a medium with zero capillary pressure. The capillary pressure inside the core is greater than zero. The capillary force at the outlet of the core disappeared. This discontinuity of capillary force is called the capillary end effect. Figure 1 is the physical schematic figure of the capillary end effect. In order to achieve capillary balance, the wetting phase will gather at the outlet end. This phenomenon continues until the saturation of the wetting phase reaches the saturation when the capillary force is zero. It will affect the saturation distribution of the wetting and non-wetting phases of the entire core, and will also affect the pressure drop across the entire core. Therefore, correcting the capillary end effect is very important for low-permeability reservoirs with strong capillary forces. Leverett [16] pointed out the existence of the capillary end effect as early as 1941. He showed that the capillary end effect makes the flow of oil and water more difficult. Osoba [17] pointed out the relative permeability curve will not be affected by the flow rate if there is no capillary end effect. If the flow of fluid in the core is affected by the capillary end effect, the relative permeability curve will change with the flow rate. Chen and Wood [18] confirmed the above viewpoint and pointed out that the relative permeability curve of mixed wet and water-wet samples has nothing to do with the flow rate. Later, Richardson [19] pointed out the capillary end effect will lead to the deviation of the relative permeability value and saturation value. The error of the saturation of the wetting phase is 15%. The error of relative permeability value of the wetting phase reaches 12%, and the error of relative permeability value of the non-wetting phase reaches 28%. Therefore, it is necessary to correct the capillary end effect.
Many scholars have made some progress in correcting the capillary end effect problem. In 1957, Richardson [20] pointed out that increasing the length of the core sample or increasing the fluid flow rate can reduce the influence of the capillary end effect. Kyte and Rapoport [21] also found in 1958 that the influence of the capillary end effect will decrease with the increase of core length and fluid flow rate. Among them, the method of increasing the length of the core is not convenient for practical operation. The method of increasing the fluid flow rate is obviously not suitable for studying the characteristics of oil-water two-phase seepage under low-speed conditions. After theoretical research, Handley [22] introduced two dimensionless parameters to quantify the influence of the capillary end effect under steady-state conditions and obtained the core pressure profile representation equation and the non-wetting phase saturation profile representation equation. Huang and Honarpour [23] proposed a new method in 1998 to correct the capillary end effect using a capillary-saturation curve. This method only accurately corrects one point on the relative permeability curve. The other data depends on the empirical formula of relative permeability so that the error is large. Moreover, the calculation process of this method is relatively complicated. In 1988, Qadeer [24] established a mathematical model for estimating the relative permeability curve of oil and water in 1988 using an optimization algorithm. He provided a formula for correcting the capillary end effect. Lei [25] further improved this formula and corrected the endpoints of the relative permeability curve. This correction formula is mainly for high permeability cores (k > 500 mD) and is not applicable to low-permeability cores.
Romanenko et al. [26] used the nuclear magnetic resonance (NMR) method to quantitatively describe the capillary end effect in 2013. However the NMR measurement has certain errors, and this method is only suitable for high permeability rocks, not for strongly heterogeneous or fractured rocks. Su [27] pointed out in a 2018 invention patent that numerical simulation method can be used to eliminate the capillary end effect and make a corresponding correction chart [28]. This method is more dependent on the selection of the relative permeability formula. There are many types of relative permeability formulas. This makes the correction result inaccurate. Ashrafi and Helalizadeh [29] proposed a method to correct the capillary end effect of carbonate rock using a genetic algorithm in 2014. In this method, the selection of the empirical formulas for relative permeability and capillary may make the correction result not unique. Gupta and Maloney [30] proposed an intercept method to correct the capillary end effect by conducting steady-state experiments under different flow rates in 2016. However, this method assumes that the proportion of the length of the end effect is the same as the proportion of the pressure intercept, although it did not confirm it. Whether the hypothesis holds is open to question. Moghaddam [31] proposed that by conducting core relative permeability experiments at four different flow rates, the relative fluid permeability and fluid saturation can be obtained. This method requires four steady-state experiments to calculate the core length affected by the capillary end effect, relative fluid permeability, and fluid saturation. This method is more accurate than the previous method, but it requires more steady-state experiments, which is difficult to popularize on a large scale.
The existing methods for correcting the capillary end effect can be divided into two categories (Table 1). The first category is represented by the correction method proposed by Huang and Honarpour [23]. This method is based on the existing empirical formula of relative permeability to adjust the relative permeability curve. This type of correction method is more dependent on the selection of the empirical formula. The correction results of this type of method are diversified and accurate correction data cannot be obtained because the empirical formula of the correlation is more diversified. The second category is represented by the correction method proposed by Rasoul Nazari Moghaddam [31]. This method corrects the capillary end effect by constructing the pressure relationship between the region affected by the capillary end effect and the region not affected by the capillary end effect inside the core. Compared with the first type of correction method, this type of method does not rely on the selection of the relative permeability formula, so the correction result is more accurate. Moghaddam’s method of measuring relative permeability experiments generally requires different proportions of injection at the same total flow rate. However, this method requires four experiments at each proportion of injection. As a result, this kind of method requires a large number of steady-state relative permeability experiments.
Based on previous work, this paper summarizes the law between the pressure drop and the length of the core affected by the capillary end effect. This paper introduces a stability factor to express the relationship between pressure drop and core length. In this paper, the relationship among core pressure drop, flow rate, and stability factor is derived. By this method, the relative permeability of the non-wetting phase that is not affected by the capillary end effect is obtained. This method simplifies the experimental process of correcting the end effect. It reduces the number of experiments without affecting the accuracy of the data. It provides a simple, effective, and feasible method for correcting the capillary end effect of the steady-state method for measuring relative permeability.

2. Materials and Methods

In the steady-state relative permeability measurement experiment, it is assumed that the two fluids are incompressible and immiscible. The following equation can be used to characterize the one-dimensional two-phase percolation law (taking the non-wetting phase as an example):
q n w = K k rnw A μ nw Δ P nw Δ L
As shown in Figure 2, the Δ P nw , CEE is the pressure drop caused by the capillary end effect. At this time, the measured pressure drop of the entire core is divided into the pressure drop not affected by the capillary end effect and the pressure drop affected by the capillary end effect:
Δ P nw , exp = Δ P nw , uneff + Δ P nw , CEE
According to Darcy’s law, we can get the relationship between the permeability and the flow rate in the region that is not affected by the capillary end effect.
q nw = K k rnw , uneff A μ nw Δ P nw , uneff Δ L uneff
Figure 2 shows the relationship between the length of each region in the core:
Δ L uneff = Δ L Δ L CEE
In Equation (4), the Δ L CEE is the length of the region affected by the capillary end effect. Substituting the Δ P nw , CEE in Equation (3) and the Δ L uneff in Equation (4) into Equation (2), the corresponding relationship between the measured pressure drop of the core and the relative permeability of the unaffected region can be obtained:
Δ P nw , exp = q nw μ nw Δ L Δ L CEE K k rnw , uneff A + Δ P nw , CEE
Since the relationship between Δ P nw , CEE and Δ L CEE in Equation (5) is not clear, we cannot solve k rnw , uneff directly. Experiments are now designed to study the relationship between the pressure drop caused by the capillary end effect and the length of the core region affected by the capillary end effect. In order to ensure the accuracy of the experimental results, the authors of this paper design three sets of steady-state experiments for measuring the oil-water relative permeability. The instrument used in the experiment is a new type of constant pressure and speed duplex pump manufactured by Jiangsu Shili Petroleum Scientific Research Instrument Co., Ltd. This instrument can inject two-phase fluid into the core at the same time, and record the inlet flow rate and pressure at any time. The oil used in this experiment is crude oil which viscosity is 20 cp at 20 °C. The water used in the experiment was pure water. The water viscosity is 1 cp at 20 °C. The core used in the experiment is the Changqing outcrop core. The core is mainly fine sand with a small amount of medium sand. As for Table 2, Table 3 and Figure 3, the composition of the core is mainly composed of quartz and calcite. Since the content of clay minerals in the core is small, the swelling effect of clay minerals has little influence on the experimental results. Figure 4 is the experimental flowchart.
The experimental steps are as follows:
  • First, the rock sample is cleaned. The original wettability of the reservoir is water wet. In the experiment, benzene and alcohol were used to clean the rock sample;
  • Then the core is evacuated and dried. Weigh the dry core and record the dry weight. Then the core is saturated with water, and the water-saturated rock sample is weighed. The wet weight of the core can be obtained. The effective pore volume and porosity are obtained by using the dry weight and wet weight of the core;
  • The oil is injected into the core which is saturated with water. Then carry out the core displacement experiment and make the core reach the state of irreducible water;
  • Record the amount of water and oil at the outlet; it prepares for calculating the relative permeability of oil and water under a certain oil-water ratio and different total flow rates. According to Table 4, oil and water are injected into the rock sample at a specific ratio. When the flow is stable (the pressure difference between the two ends of the rock sample is stable), we record the inlet and outlet pressures. Then use a density meter to get the water cut of the liquid at the outlet. The water saturation of the core is obtained according to the principle of material balance;
  • Under the condition of the same proportion of injected liquid (wetting phase and non-wetting phase), change the total flow rate. Systematically varying the flow rate, repeat the previous step until the end of the experimental run with that imposed ratio between phases.
Table 4 is a summary table of flow rates under three different injections rates. A total of 12 displacement experiments were carried out, for four total flow rates, 0.3, 0.6, 1.1, and 2.2 mL/min, and each of them for three oil-water ratios, 5:1, 1:5, and 1:10. After each displacement run with fixed phase ratios, experiment 3 is repeated in order to ensure the accuracy of the procedure; corresponding tests are labeled 5, 10, and 15 in Table 5. It is observed that each time the results do not have significant change and the experiments are reproducible. This can lead to two conclusions. First, the core pore structure has not changed; second, the experimental data is accurate. We make sure that the internal structure of the core has not changed by comparing the changes in flow rate, pressure difference, and water saturation. This method can ensure the accuracy of subsequent experimental results. In this paper, 12 displacement experiments and 3 repeated experiments were carried out.

3. Results

In the experiment, the pressure and total flow rate are recorded; then by using the intercept method [30], Δ P nw , CEE and Δ L nw , CEE are evaluated, in order to explore the correlation between them. Draw the relationship between Δ P nw , CEE Δ P nw , exp Δ P nw , CEE and Δ L CEE Δ L Δ L CEE according to the corrected result. The intercept method is shown in Appendix A. It can be found from Figure 5a that under the PI (the same proportion of injected liquid (the wetting phase and the non-wetting phase), the same ratio between Δ P nw , CEE Δ P nw , exp Δ P nw , CEE and Δ L CEE Δ L Δ L CEE can be obtained. There is a stable proportional relationship between Δ P nw , CEE Δ P nw , exp Δ P nw , CEE and Δ L CEE Δ L Δ L CEE . We have made a linear fit between Δ P nw , CEE Δ P nw , exp Δ P nw , CEE and Δ L CEE Δ L Δ L CEE . The fitting coefficients are all greater than 0.95, and the fitting effect is very satisfactory. Now the ratio of Δ P nw , CEE Δ P nw , exp Δ P nw , CEE to Δ L CEE Δ L Δ L CEE is defined as the stability factor, as shown in Equation (6).
Δ P nw , CEE Δ P nw , exp - Δ P nw , CEE ÷ Δ L CEE Δ L Δ L CEE = β
The capillary end effect factor is now defined to represent the ratio of the length affected by the capillary end effect to the total length of the core, as like α = Δ L C E E / Δ L . The introduction of the concept of stability factor (β) in this paper can reduce the number of experiments for the subsequent correction of capillary end effects. We can obtain the relationship between the capillary end effects factor and the stability factor:
α = Δ P nw , CEE Δ P nw , CEE Δ P nw , CEE β + Δ P nw , CEE
Then Equation (7) is incorporated into Equation (5) to obtain the corresponding relationship between the measured pressure drop of the core and the stability factor:
Δ P nw , exp = Δ P nw , CEE + q nw μ nw Δ L K A × Δ P nw , exp Δ P nw , CEE β Δ P nw , exp Δ P nw , CEE β + Δ P nw , CEE k rnw , uneff
After finishing Equation (8), we can get:
Δ P nw , exp = q nw μ nw Δ L K A 1 k rnw , uneff + Δ P nw , CEE Δ P nw , CEE β
Moghaddam [31] pointed out that the value of relative permeability is related to the proportion of injected liquid and does not depend on the flow rate. According to Equation (9) and two steady-state displacement experiments of oil and water, the values of k rnw , uneff and Δ P nw , CEE Δ P nw , CEE β can be obtained.
Δ P nw , exp 1 = q nw 1 μ nw Δ L K A × 1 K rnw , uneff + Δ P nw , CEE Δ P nw , CEE β Δ P nw , exp 2 = q nw 2 μ nw Δ L K A × 1 K rnw , uneff + Δ P nw , CEE Δ P nw , CEE β
According to Equation (10), take q nw i μ nw Δ L K A as the abscissa and Δ P nw , e x p i as the ordinate to establish a rectangular coordinate system. The reciprocal of the slope is k rnw , uneff , and the intercept is Δ P nw , CEE Δ P nw , CEE β . Therefore, the relative permeability of the water phase that is not affected by the capillary end effect can be obtained by two steady-state experiments. This method reduces the workload by half compared with the previous experimental design [31].
The stability factor is used to calculate the relative permeability of the wetting phase and the non-wetting phase in this paper. The repeated process is as follows. First of all, we need to test the relative permeability of oil and water with two steady-state experiments of the core. It should be noted that the two experiments need to be carried out at the same oil-water proportion. After the experiment, we can get Δ P nw , e x p 1 , Δ P nw , e x p 2 , q nw 1 and q nw 2 . Second, we can get k rnw , uneff and Δ P nw , CEE Δ P nw , CEE β . Since this method is calculated based on the stability factor, this method is named the stability factor method. And the correction results are compared with the intercept method proposed by Gupta and Maloney [30], as shown in Figure 5b. Figure 5b shows the comparison of the relative permeability of the wetting phase between the stability factor method and the intercept method under different total flow rates. Under the same oil-water ratio, the relative permeability should not change with the total flow rate. The intercept method obviously does not conform to this law. This method believes that the seepage capacity of the A and B (Figure 2) is the same. This view is obviously wrong. It can be found from Figure 5b that whether it is for the wetting phase or the non-wetting phase, the stability factor method is free from these drawbacks and can reasonably be considered as a more accurate one. In order to ensure the accuracy of the experimental results, this paper uses the correction results of Moghaddam [31] to compare the results of the method in this paper, as shown in Figure 6. As mentioned earlier, the correction result of Moghaddam [31] is accurate, but the number of experiments is more than other methods. This paper will use this method as a standard to measure the accuracy of the correction result of the method in this paper. It can be seen from Figure 6 that the correction result error between the stability factor method and the standard method is small, which further demonstrates the accuracy of the stability factor method.
Figure 7 shows the comparison among the original relative permeability curve, the relative permeability curve corrected by the intercept method, and the stability factor method. The OR/IC/SF in Figure 7 are the original experimental data, the repeated result of the intercept method, and the repeated result of the stability factor method. It is observed that the relative permeability value of the corrected wetting phase and non-wetting phase is significantly enhanced. The end effect hinders the flow of the non-wetting phase. Affected by the capillary end effect, the residual oil saturation of the non-wetting phase increased by 30%. If the error of capillary end effect is not corrected, accurate relative permeability data cannot be provided for reservoir production.

4. Discussion

To ensure the accuracy of the conclusions obtained, we use the experimental data of Moghaddam [31] to verify the feasibility of the stability factor method. The experimental data of Moghaddam [31] is now used to study the relationship between the pressure drop caused by the capillary end effect and the length of the core region affected by the capillary end effect.
The basic data of the core is shown in Table 6.
The pressure drop and flow data at both ends of the core measured in the experiment are shown in Table 7.
Table 7 is a summary table of pressure and flow rate under two different injection ratios of wetting phase and non-wetting phase. In each injection ratio, steady-state methods to measure the relative permeability of oil and water experiments were carried out at four different flow rates. The flow rate and pressure drop of each experiment were recorded separately. Table 8 shows the corrected pressure drop and core length in the region where the capillary end effect occurs in the experimental data of Moghaddam [31]. Each of the same proportion of injected liquid of wetting phase and non-wetting phase conditions has the same value of Δ P nw , CEE . According to the corrected results, the relationship between Δ P nw , CEE Δ P nw , exp Δ P nw , CEE and Δ L CEE Δ L Δ L CEE is shown in Figure 8a.
It can be found from Table 8 that under the same injection ratio, there is a certain proportional relationship between Δ P nw , CEE and Δ L CEE . With the difference in the total flow rate, the stability factor is constant. It is only related to the proportion of injected liquid of the wetting phase to the non-wetting phase in the steady-state method of measuring the relative permeability of oil and water and has nothing to do with the total flow rate.
It is necessary to verify whether the reduction of the number of experiments has an impact on the accuracy of the experimental results. In this paper, the corrected result calculated by the stability factor is compared with the original result. Table 9 shows the specific numerical values of x i , y i , k rnw , eff and Δ P nw , CEE Δ P nw , CEE β . In the table, the unit of x i and y i is atm.
In Figure 8b, q nw i μ nw Δ L K A is the abscissa and Δ P nw , exp i is the ordinate. According to the stability factor, it can be found that q nw i μ nw Δ L K A and Δ P nw , exp i have an obvious linear relationship. The reciprocal of the slope in the figure is the relative permeability of the non-wetting phase. The intercept is the result of the interaction of Δ P nw , exp i and q nw i μ nw Δ L K A . Under the same proportion of injected liquid of wetting and non-wetting phases, they have the same slope and intercept. That is to say, there is no need for the four-time steady-state method to measure relative permeability experimental data. Only two sets of experiments are needed to obtain the relative permeability of the non-wetting phase that is not affected by the capillary end effect. The linear relationship between Δ P nw , exp i and q nw i μ nw Δ L K A under the four flow rates can be found to be consistent with the two groups when calculated separately. These four points are exactly on the same straight line. The R2 of each fitted straight line is above 1. It shows that the change of the proportion of injected liquid of the wetting phase and the non-wetting phase does not affect the linear law.
From Figure 8b, it can be concluded that how to change the proportion of injected liquid of the wetting phase and the non-wetting phase. As long as the proportion of injected liquid is the same, there is no need for the previous four steady-state experiments to measure relative permeability. In the method shown in this paper, only two sets of experimental data are needed to obtain the relative permeability of the non-wetting phase that is not affected by the capillary end effect. This method reduces the number of experiments and greatly saves time and cost.
As shown in Table 10, k rnw is the relative permeability value of the non-wetting phase obtained by the stability factor method. The k rnw * is the relative permeability value obtained by Moghaddam’s method. The k rnw * is the absolute value between the two errors. The magnitude of k rnw * k rnw is much smaller than k rnw and k rnw * . Therefore, the error caused by the stability factor method is within the allowable range of error. It is proved that this method not only reduces the number of experiments but also calculates the relative permeability of the non-wetting phase very accurately. The problem of correcting the capillary end effect of relative permeability under low-velocity seepage conditions can be dealt with more efficiently.

5. Conclusions

The stability factor method summarizes the law between the pressure drop caused by the end effect and the core region. It can reduce the number of correction experiments by half, thereby reducing the manpower, material resources, and time required. Research shows that under the same injection ratio, the core has the same stability factor in the experiment of measuring the relative permeability of oil and water by the steady-state method. The stability factor does not change with the total flow rate. In addition, under the same injection ratio, the stability factor method only requires two steady-state methods to measure the relative permeability experiment to obtain the corrected relative permeability value of the non-wetting phase. It not only guarantees the accuracy of the results (the error is less than 0.01) but also simplifies the experimental process.

Author Contributions

Conceptualization, Y.L. and S.W.; methodology, Y.L. and S.W.; validation, Y.L.; formal analysis, Y.L.; investigation, Y.L. and S.W.; resources, Q.Y., X.X. and C.Y.; data curation, X.Z.; writing—original draft preparation, Y.L.; writing—review and editing, Z.K., Q.Y., X.X. and C.Y.; visualization, Y.L.; supervision, S.W.; project administration, S.W.; funding acquisition, Y.L. and S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Open Foundation of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Southwest Petroleum University).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

List of Symbols

q n w Non-wetting phase flow rates, mL/min
Δ L The length of the core, cm
K Absolute permeability, µm2
A Area, cm2
β Stability factor
i Number of experiments
Δ L CEE Core length affected by capillary end effect, cm
Δ P nw , exp Pressure drop of the non-wetting phase measured in the experiment, atm
Δ P nw , uneff Pressure drop unaffected by the capillary end effect, atm
Δ P nw , CEE The pressure drop affected by the capillary end effect, atm
k rnw , uneff Relative permeability of non-wetting phase unaffected by capillary end effect
k rnw Relative permeability of non-wetting phase calculated by the stability factor method
k rnw * The relative permeability of the non-wetting phase calculated by Moghaddam’s method

Appendix A

Gupta and Maloney [30] proposed an intercept method to correct the capillary end effect by conducting steady-state experiments under different flow rates in 2016. This method requires running a steady-state relative permeability test with several different flow rates for each fractional flow. First, take the non-wetting phase as an example to establish Darcy’s formula as follows:
q t ( 1 F w ) = K k rnw , uneff A μ nw Δ L Δ P nw , exp Δ P nw , CEE
Rearranging Equation (A1) gives:
Δ P nw , exp = q t μ nw Δ L K k rnw , uneff A ( 1 F w ) + Δ P nw , CEE
The Δ P nw , exp can be expressed as the pressure drop of the non-wetting phase measured in the experiment. The Δ P nw , CEE can be expressed as the pressure drop affected by the capillary end effect. The Δ P nw , uneff is the pressure drop unaffected by the capillary end effect. The q t is the total flow rate and the F w is the wetting phase fractional flow. Using the above concept, Gupta and Maloney proposed to obtain Δ P nw , CEE from the intercept of the plot of Δ P nw , exp and the injected total flow rate ( q t ).
To correct the saturation data, they used the following overall saturation balance equation for a given fractional flow condition:
Δ L S nw , exp = Δ L Δ L CEE S n w , u n e f f + S nw , CEE Δ L CEE
The S nw , exp can be expressed as the saturation of the non-wetting phase measured in the experiment. The S nw , CEE can be expressed as the saturation affected by the capillary end effect. The S n w , u n e f f is the saturation unaffected by the capillary end effect.
And they proposed a capillary end effect length factor β, and hypothesized:
β = Δ L CEE Δ L = Δ P nw , CEE Δ P nw , exp Δ P nw , CEE
Rearranging Equation (A3) and using Equation (A4) gives the expression:
1 1 β S nw , exp = S nw , CEE β 1 β + S n w , u n e f f
Based on Equation (A5), S n w , u n e f f is the intercept of the plot of 1 1 β S nw , exp and β 1 β .

References

  1. Zhu, X.M.; Pan, R.; Sun, S.F. Research progress and core issues in tight reservoir exploration. Earth Sci. Front. 2018, 25, 141–146. [Google Scholar]
  2. Weng, D.; Yan, X.; Lu, Y.; Peng, X. Optimization and realization of stress interference in tight oil and gas reservoir. J. China Univ. Min. Technol. 2014, 43, 639–645. [Google Scholar]
  3. Fan, D.; Yao, J.; Sun, H.; Zeng, H. A composite model of hydraulic fractured horizontal well with stimulated reservoir volume in tight oil & gas reservoir. J. Nat. Gas Sci. Eng. 2015, 24, 115–123. [Google Scholar]
  4. Liang, T.; Guo, J.C.; Zhao, Z.H.; Yin, Q.W. Refracturing candidate selection for MFHWs in tight oil and gas reservoirs using hybrid method with data analysis techniques and fuzzy clustering. J. Cent. South Univ. 2020, 27, 277–287. [Google Scholar]
  5. Hu, H.Y. Accumulation mechanism of ultro-low permeability sandstone reservoir, Huaqing oilfield, Ordas Basin. Adv. Mater. Res. 2014, 868, 66–69. [Google Scholar] [CrossRef]
  6. Zou, C. Tight Oil and Gas. In Unconventional Petroleum Geology, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2017; pp. 239–273. [Google Scholar]
  7. Ren, X.; Li, A.; Fu, S.; Tian, W. Influence of micro-pore structure in tight sandstone reservoir on the seepage and water-drive producing mechanism—A case study from Chang 6 reservoir in Huaqing area of Ordos basin. Energy Sci. Eng. 2019, 7, 741–753. [Google Scholar] [CrossRef] [Green Version]
  8. Han, B.; Bian, X.Q. A hybrid PSO-SVM-based model for determination of oil recovery factor in the low-permeability reservoir. Petroleum 2017, 4, 43–49. [Google Scholar] [CrossRef]
  9. Bryant, S.; Blunt, M. Prediction of relative permeability in simple porous media. Phys. Rev. A 2015, 46, 2004–2011. [Google Scholar] [CrossRef] [PubMed]
  10. Li, B.; Pu, W.F.; Li, K.X.; Jia, H.; Wang, K.Y.; Yang, Z.G. The characteristics and impacts factors of relative permeability curves in high temperature and low-permeability limestone reservoirs. Adv. Mater. Res. 2014, 1010–1012, 1676–1683. [Google Scholar] [CrossRef]
  11. Berg, S.; Unsal, E.; Dijk, H. Sensitivity and uncertainty analysis for parameterization of multiphase flow models. Transp. Porous Media 2021, 1–31. [Google Scholar]
  12. Feng, Y.L.; Zhang, Y.; Ji, B.Y.; Mu, W.Z. Micro-acting force in boundary layer in low-permeability porous media. Chin. Phys. Lett. 2011, 28, 024703. [Google Scholar] [CrossRef]
  13. Zhang, D.; Tan, J.; Yang, D.; Mu, S.; Peng, Q. The residual potential of bottom water reservoir based upon genetic algorithm for the relative permeability inversion. J. Geosci. Environ. Prot. 2019, 07, 192–201. [Google Scholar] [CrossRef] [Green Version]
  14. Varela, O.J.; Torres-Verdín, C.; Lake, L.W. On the value of 3D seismic amplitude data to reduce uncertainty in the forecast of reservoir production. J. Pet. Sci. Eng. 2006, 50, 269–284. [Google Scholar] [CrossRef]
  15. Hamoud, M.; Chalaturnyk, R.; Leung, J. Experimental studies of shear-induced changes in two-phase relative permeability. Arab. J. Geosci. 2020, 13, 1–22. [Google Scholar] [CrossRef]
  16. Leverett, M.C. Capillary Behavior in porous solids. Transactions 1941, 142, 152–169. [Google Scholar] [CrossRef]
  17. Osoba, J.S.; Richardson, J.G.; Kerver, J.K.; Hafford, J.; Blair, P. Laboratory measurements of relative permeability. J. Pet. Technol. 1951, 3, 47–56. [Google Scholar] [CrossRef]
  18. Chen, A.L.; Wood, A.C. Rate Effects on water-oil relative permeability. In Proceedings of the International Symposium of the Society of Core Analysts, Edinburgh, Scotland, 17 September 2001; pp. 17–19. [Google Scholar]
  19. Richardson, J.G.; Kerver, J.K.; Hafford, J.A.; Osoba, J. Laboratory determination of relative permeability. J. Pet. Technol. 1952, 4, 187–196. [Google Scholar] [CrossRef]
  20. Richardson, J.G. The calculation of waterflood recovery from steady-state relative permeability data. J. Pet. Technol. 1957, 9, 64–66. [Google Scholar] [CrossRef]
  21. Kyte, J.R.; Rapoport, L.A. Linear waterflood behavior and end effects in water-wet porous media. J. Pet. Technol. 1958, 10, 47–50. [Google Scholar] [CrossRef]
  22. Hadley, G.F.; Handy, L.L. A Theoretical and experimental study of the steady state capillary end effect. J. Pet. Technol. 1956, 707-G. [Google Scholar]
  23. Huang, D.D.; Honarpour, M.M. Capillary end effects in coreflood calculations. J. Pet. Sci. Eng. 1998, 19, 103–117. [Google Scholar] [CrossRef]
  24. Qadeer, S.; Dehghani, K.; Ogbe, D.O.; Ostermann, R.D. Correcting oil/water relative permeability data for capillary end effect in displacement experiments. In Proceedings of the SPE California Regional Meeting, San Jose, CA, USA, 1 January 1988. [Google Scholar]
  25. Lei, X.; Zha, Y.Q.; Jiang, P.; Zhang, Q.L.; Li, F.Y. An improved correction method for end-effect in oil-water relative permeability experiment. China Offshore Oil Gas 2016, 28, 49–53. [Google Scholar]
  26. Romanenko, K.; Balcom, B.J. An assessment of non-wetting phase relative permeability in water-wet sandstones based on quantitative MRI of capillary end effects. J. Pet. Sci. Eng. 2013, 110, 225–231. [Google Scholar] [CrossRef]
  27. Li, R.; Wang, J. Study on Numerical simulation of oil-water phase non-darcy flow in low-permeability reservoirs. Open Access Libr. J. 2016, 3, 1–8. [Google Scholar] [CrossRef]
  28. Su, Y.L.; Cheng, A.Q.; Zhan, S.Y.; Li, D.; Sheng, G. The calculating method for correcting end effect and relative permeability calculation error in tight core. Sci. Technol. Eng. 2018, 18, 180–185. [Google Scholar]
  29. Ashrafi, M.; Helalizadeh, A. Genetic algorithm for estimating relative permeability and capillary pressure from unsteady-state displacement experiments including capillary end-effect. Energy Sources Part A Recovery Util. Environ. Eff. 2014, 36, 2443–2448. [Google Scholar] [CrossRef]
  30. Gupta, R.; Maloney, D.R. Intercept method-a novel technique to correct steady-state relative permeability data for capillary end effects. SPE Reserv. Eval. Eng. 2014, 19, 316–330. [Google Scholar] [CrossRef]
  31. Moghaddam, R.N.; Jamiolahmady, M. Steady-state relative permeability measurements of tight and shale rocks considering capillary end effect. Transport. Porous Media 2019, 128, 1–22. [Google Scholar] [CrossRef]
Figure 1. Physical schematic figure of the capillary end effect. In order to achieve capillary balance, the wetting phase (blue) will gather at the outlet end. This phenomenon continues until the saturation of the wetting phase reaches the saturation when the capillary force is zero. However, the non-wetting phase (red) will not accumulate at the outlet end.
Figure 1. Physical schematic figure of the capillary end effect. In order to achieve capillary balance, the wetting phase (blue) will gather at the outlet end. This phenomenon continues until the saturation of the wetting phase reaches the saturation when the capillary force is zero. However, the non-wetting phase (red) will not accumulate at the outlet end.
Energies 14 04528 g001
Figure 2. Diagram of the internal conditions of the core (according to whether it is affected by the capillary end effect, the core is divided into two regions: A and B. The fluid in the core is divided into a wetting phase (blue) and a non-wetting phase (red)).
Figure 2. Diagram of the internal conditions of the core (according to whether it is affected by the capillary end effect, the core is divided into two regions: A and B. The fluid in the core is divided into a wetting phase (blue) and a non-wetting phase (red)).
Energies 14 04528 g002
Figure 3. Diagram of the internal conditions of the core.
Figure 3. Diagram of the internal conditions of the core.
Energies 14 04528 g003
Figure 4. Experimental flowchart (the wetting phase and the non-wetting phase are injected into the core through a pump. The flow rates at the inlet and outlet ends of the core are measured and transmit the measurement data to the control system).
Figure 4. Experimental flowchart (the wetting phase and the non-wetting phase are injected into the core through a pump. The flow rates at the inlet and outlet ends of the core are measured and transmit the measurement data to the control system).
Energies 14 04528 g004
Figure 5. PI = 0.1 (orange); PI = 0.2 (green); PI = 5 (blue). The calculation and application of the stability factor: (a) Calculation of stability factor; (b) comparison of wetting phase correction results between the stability factor method (SF) and the intercept method (IC).
Figure 5. PI = 0.1 (orange); PI = 0.2 (green); PI = 5 (blue). The calculation and application of the stability factor: (a) Calculation of stability factor; (b) comparison of wetting phase correction results between the stability factor method (SF) and the intercept method (IC).
Energies 14 04528 g005
Figure 6. Comparison of the correction results of the wetting phase and the non-wetting phase of the stability factor method (SF) and the method proposed by Moghaddam (SS) when the PI = 0.1, PI = 0.2, and PI = 5.
Figure 6. Comparison of the correction results of the wetting phase and the non-wetting phase of the stability factor method (SF) and the method proposed by Moghaddam (SS) when the PI = 0.1, PI = 0.2, and PI = 5.
Energies 14 04528 g006
Figure 7. Comparison among the original relative permeability curve (OR), the relative permeability curve corrected by the intercept method (IC), and the stability factor method (SF).
Figure 7. Comparison among the original relative permeability curve (OR), the relative permeability curve corrected by the intercept method (IC), and the stability factor method (SF).
Energies 14 04528 g007
Figure 8. Application examples of stability factor: (a) Calculation of stability factor; (b) application of stability factor.
Figure 8. Application examples of stability factor: (a) Calculation of stability factor; (b) application of stability factor.
Energies 14 04528 g008
Table 1. Category of correction methods for capillary end effects.
Table 1. Category of correction methods for capillary end effects.
CategoryThe First CategoryThe Second Category
CharacteristicThis type of correction method is more dependent on the selection of the empirical formula.This method corrects the capillary end effect by constructing the pressure relationship between the region affected by the capillary end
effect and the
region not affected by the capillary end effect inside the core.
AdvantageThis kind of method requires once unsteady-state relative permeability experiments.This type of method does not rely on the selection of the relative permeability formula, so the correction
result is more accurate.
DisadvantageThe correction results of this type of method are diversified and accurate correction data cannot be obtained because the empirical formula of the correlation is more diversified.This kind of method requires a large
number of steady-state relative permeability experiments.
Reference number[23,24,25,26,27,28][30,31]
Table 2. Basic data of the core.
Table 2. Basic data of the core.
Core NumberLength/cmDiameter/cmArea/cm2
60-17.9902.5401.613
Porosity/%Absolute permeability/DDry weight/gWet weight/g
19.7780.063184.73993.027
Table 3. X-ray diffraction test results table.
Table 3. X-ray diffraction test results table.
Core NumberComponent Content (%)
QuartzPotash FeldsparPlagioclaseCalciteDolomite
60-16815251
PyriteBariteAnhydriteAnalciteHematite
00000
Table 4. Experimental design table.
Table 4. Experimental design table.
Experiment NumberProportion of Injected LiquidTotal Flow Rates/(mL/min) q nw / ( mL / min ) q w / ( mL / mim )
11:100.30.0270.273
21:100.60.0550.545
31:101.10.1001.000
41:102.20.2002.000
51:101.10.1001.000
65:010.30.250
0.500
0.050
0.100
75:010.6
85:011.10.917
1.833
0.183
0.367
95:012.2
101:101.10.1001.000
111:050.30.050
0.100
0.250
0.500
121:050.6
131:051.10.183
0.367
0.917
1.833
141:052.2
151:101.10.1001.000
Table 5. Data from repeated experiments.
Table 5. Data from repeated experiments.
Experiment NumberProportion of Injected LiquidPressure Drop/MPaWater Saturation
31:107.4200.420
51:107.9240.412
101:107.9910.437
151:107.2100.457
Table 6. Basic core data [31].
Table 6. Basic core data [31].
Core Number Δ L / cm D / cm A / cm 2 Φ / % K / μ D
-4.613.7711.1613.59.9
Table 7. Experimental data of core affected by capillary end effect under different flow rates.
Table 7. Experimental data of core affected by capillary end effect under different flow rates.
Experiment Number q w / q nw q nw / mL min 1 Δ P nw , exp / atm
10.080.073.95
20.084.82
30.2513.51
40.5026.55
50.020.186.46
60.258.42
70.4213.31
80.8325.53
Table 8. Stability factor calculation table.
Table 8. Stability factor calculation table.
Number Δ P nw , CEE / atm Δ L CEE / cm β
12.873.181.20
22.872.541.20
32.870.851.20
42.870.421.20
56.044.251.22
66.043.111.22
76.041.871.22
86.040.931.22
Table 9. Calculation results of relative permeability.
Table 9. Calculation results of relative permeability.
Number y i x i k rnw , uneff Δ P nw , CEE Δ P nw , CEE β
13.950.00540.15100.048
24.820.0067
313.510.02020.15100.048
426.550.0404
56.460.01480.26850.110
68.420.0202
713.310.03370.26850.110
825.530.0674
Table 10. The experimental data of Moghaddam is corrected for the end effect. Comparison of Moghaddam’s method and the relative permeability correction result of the stability factor method.
Table 10. The experimental data of Moghaddam is corrected for the end effect. Comparison of Moghaddam’s method and the relative permeability correction result of the stability factor method.
Number k rnw k rnw * k rnw * k rnw
10.1550.1553.53 × 10−6
2
30.1556.8 × 10−5
4
50.2760.1559.4 × 10−5
6
70.2762.5 × 10−5
8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Y.; Wang, S.; Kang, Z.; Yuan, Q.; Xue, X.; Yu, C.; Zhang, X. Research on the Correction Method of the Capillary End Effect of the Relative Permeability Curve of the Steady State. Energies 2021, 14, 4528. https://doi.org/10.3390/en14154528

AMA Style

Li Y, Wang S, Kang Z, Yuan Q, Xue X, Yu C, Zhang X. Research on the Correction Method of the Capillary End Effect of the Relative Permeability Curve of the Steady State. Energies. 2021; 14(15):4528. https://doi.org/10.3390/en14154528

Chicago/Turabian Style

Li, Yanyan, Shuoliang Wang, Zhihong Kang, Qinghong Yuan, Xiaoqiang Xue, Chunlei Yu, and Xiaodong Zhang. 2021. "Research on the Correction Method of the Capillary End Effect of the Relative Permeability Curve of the Steady State" Energies 14, no. 15: 4528. https://doi.org/10.3390/en14154528

APA Style

Li, Y., Wang, S., Kang, Z., Yuan, Q., Xue, X., Yu, C., & Zhang, X. (2021). Research on the Correction Method of the Capillary End Effect of the Relative Permeability Curve of the Steady State. Energies, 14(15), 4528. https://doi.org/10.3390/en14154528

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop