Multiscale and Multiphysical Approaches to Fluids Flow in Unconventional Reservoirs
View this Special IssueResearch Article  Open Access
Huiyuan Bian, Kewen Li, Binchi Hou, Xiaorong Luo, "A New Model to Calculate OilWater Relative Permeability of Shaly Sandstone", Geofluids, vol. 2020, Article ID 8842276, 11 pages, 2020. https://doi.org/10.1155/2020/8842276
A New Model to Calculate OilWater Relative Permeability of Shaly Sandstone
Abstract
Oilwater relative permeability curves are the basis of oil field development. In recent years, the calculation of oilwater relative permeability in sandstone reservoirs by resistivity logging data has received much attention from researchers. This article first analyzed the existing mathematical models of the relationship between relative permeability and resistivity and found that most of them are based on Archie formula, which assumes the reservoir is clean sandstone. However, in view of the fact that sandstone reservoir is commonly mixed with shale contents, this research, based on the dual water conductivity model, Poiseuille’s equation, Darcy’s law, and capillary bundle model, derived a mathematical model (DW relative permeability model) for shaly sandstone reservoir, which calculates the oilwater relative permeability with resistivity. To test and verify the DW relative permeability model, we designed and assembled a multifunctional core displacement apparatus. The experiment of core oilwater relative permeability and resistivity was designed to prove the effectiveness of the DW relative permeability model in shaly sandstone reservoirs. The results show that the modified Li model can well express the transformational relation between resistivity and relative permeability in sandstone reservoir with low clay content. Compared with the modified Li model and the Pairoys model, the DW relative permeability model is more helpful to collect better results of relative permeability in shaly sand. These findings will play a significant role in the calculation of oilwater relative permeability in reservoirs based on resistivity logging data and will provide important data and theory support to the shaly sandstone reservoir characterized oil field development.
1. Introduction
The evaluation of tight reservoir has always been an important part of petroleum geology research, while oilwater relative permeability, which is vital to the evaluation of fluid flow in porous media, is used in all aspects of the reservoir engineering [1, 2]. Traditionally, relative permeability is obtained in laboratory. However, in many cases, especially in low permeability reservoirs, or when phase transformation or mass transfer happens with the change of pressure, oilwater relative permeability experiments are difficult, expensive, and time consuming simultaneously [3, 4]. Alternately, it is difficult to maintain the samples the same as in the reservoirs; moreover, relative permeability is almost impossible to obtain in real time. Despite all these difficulties, experiment serves as the main method to calculate relative permeability curves for oil fields.
Conventional resistivity logging data, which is the basic information of oil and gas well standard logging, is in large amount and available. In recent years, more and more scholars indicated that there is a relationship between relative permeability and resistivity [5–9]. Cai et al. [10] presented a review of the electrical conductivity models using fractal, percolation, and effective medium theories. In another article, Cai et al. [11] proposed a combined model including porethroat ratio, tortuosity, and connectivity, exactly estimating the influence of complex pore structure on the transport behavior associated with electrical parameters. Li [12] based on Archie formula put forward a mathematical model that uses resistivity to calculate the relative permeability of gaswater, and Li verified the model with experimental data. Li et al. [13, 14] conducted a lot of research on this field. Li together with Horne and Williams worked out methods to calculate twophase relative permeability with resistivity logging data in uniform medium. Mohammed and Birol [15] modified Li model by taking the fluid viscosity and the average water saturation at the time of water breakthrough into consideration. Alex et al. [16] proposed a method using resistivity to calculate relative permeability in dual porosity model, but the model has not been verified experimentally. Pairoys et al. [17] verified the Li model and BrooksCorey model [6] with gaswater relative permeability experiments and found that the Li model works better than BrooksCorey model in that situation. Then, the Li model is modified by replacing the pore size distribution index to index saturation exponent Li model [18]. Pairoys [19] analyzed the change of resistivity under different frequency in the process of unsteady twophase flow displacement, based on which the Li model was verified again with gaswater relative permeability experimental data and oil flooding data by Bian and Li [20]. The abovemodified models based on the Li model were established under the condition of homogeneous clean sandstone reservoir. However, most real sandstone reservoirs contain shale contents, which influence rock resistivity and relative permeability significantly.
This study, based on the dual water conductivity model, Poiseuille’s equation, Darcy’s law, and capillary bundle model, proposed and verified a mathematical model (DW relative permeability model) to calculate relative permeability using resistivity. To improve the Li model, a new model named “dual water relative permeability model (DW model)” was proposed in the consideration of better expressing the transformational relations between resistivity and relative permeability in the shaly sand reservoir. According to experiments, the DW model achieved the goal of reflecting the relation between resistivity and relative permeability in a better way than the modified Li model and Pairoys model, which is helpful in both the calculation of oilwater relatively permeability in shaly sand reservoir based on resistivity and the oil field development.
2. Mathematical Background
2.1. Relationship between Water Saturation and Relative Permeability
There exist many relationship models between water saturation and relative permeability, among which the most common one is as shown below [20]: where and are relative permeabilities of the wetting and nonwetting phase, and are the saturation and the irreducible saturation of wetting phase, is the maximum when , is the maximum when , is the residual saturation of nonwetting phase, and is the normalized saturation of wetting phase.
2.2. Li Model for the Relationship between Relative Permeability and Resistivity
Fluid flow in porous media is similar to current flow in conductive media [13]. According to the Li model, gas/water relative permeability is calculated using resistivity. where is the resistivity index, and is the wettingphase normalized relative permeability.
2.3. Pairoys Model
Many modified models grew out of the Li model that is suitable for gaswater twophase flow and oil flooding, but not for water flooding. Pairoys worked out the following model after analyzing water flooding situation [19]. where is the saturation of the wetting phase, is the irreducible saturation of the wetting phase, is the residual saturation of the nonwetting phase, is the normalized saturation of the wetting phase, is the resistivity index, when , and and are normalized relative permeabilities of the wetting and nonwetting phases is the resistivity index.
2.4. Modified Li Model
Based on the Li model, Bian and Li proposed a model for the relationship between resistivity and oilwater relative permeability of water wet sandstone reservoirs with low shaly contents [20]. where is the formation resistivity when and and are normalized relative permeabilities of the water and oil phases.
3. Relationship between Resistivity and Relative Permeability of Shaly Sand Reservoir
The relationship models of resistivity and relative permeability mentioned above all assume that the reservoirs are homogeneous and pure sandstone. However, in reality, most sandstone reservoirs contain shale contents. Therefore, in order to calculate the relative permeability in shaly sand reservoirs accurately, a new model suitable for shaly sandstone should be established.
The crosssectional area, length, and volume of the water wet shaly sandstone are , , and , respectively (as shown in Figure 1(a)). The effective pore space of rock is considered to be composed of large bore capillary columns with equal crosssectional area and small bore capillary columns with equal crosssectional area. The large columns are filled with movable water and oil, while the small columns are filled with immovable water (irreducible water) and residual oil. The crosssectional area, crosssection radius, length, and volume of the large bore capillary columns are , , , and , while those of the small columns are , , , and , respectively (as shown in Figure 1(b)).
(a) Shaly sand core
(b) Capillary columns pack bundle (water wetting)
(c) The cross section of capillary columns pack bundle (water wetting)
When the water saturation of the rock is , in the () large capillary columns, the oil crosssectional area and oil crosssection radius are and , while the crosssectional area, length, and volume of the movable water are , , and , respectively. In the () small capillary columns, the oil crosssectional area and oil crosssection radius are and , while the crosssectional area, length, and volume of the immovable water are , , and , respectively. Due to the existence of shale contents, it is assumed that the immobile water in the small capillary columns contains clay water. The crosssectional area, length, and volume of the clay water are , , and , respectively (as shown in Figure 1(c)).
When the core sample is saturated with water, according to Poiseuille flow formula, the liquid flow in the () large capillary column is calculated as follows. where is the pressure difference, and is the fluid viscosity.
The total flow of water in the rock is
According to Darcy’s formula,
Suppose that the length of the large capillary column is equal to that of the small capillary column. The permeability is obtained as follows.
Similarly, when the water saturation is , the small capillary columns are filled with bound water and residual oil, while the large capillary columns are filled with movable water and movable oil.
According to the Poiseuille flow formula, the total flow of movable water in the rock is
According to Darcy’s formula,
So, the permeability is
And the water relative permeability is obtained as follows.
The electrical conductivity of the model is analyzed below. The resistivity of free water is , while the resistivity of clay water is . In the () small capillary columns, the crosssectional area, length, and volume of the bound free water are , , and , respectively. When the core sample is saturated with water, the resistivity is . When the saturation is and the resistivity is , there is the following equation.
Assuming that all the capillary columns have the same length,
The resistivity index is
So, the water relative permeability can be obtained.
Define the resistivity .
Define can be calculated as follows.
Therefore, the normalized water relative permeability in the dual water relative permeability model can be expressed as follows:
When , . When , and , which satisfies the boundary condition.
When the clay water content is 0 (), the model is simplified to a clean sandstone model, and the normalized water relative permeability is as follows.
3.1. Determination of the Parameters in the DW Relative Permeability Model
3.1.1. Calculation of the Resistivity of Clay Water
Diffusion factor of Na^{+} ion diffusion layer is calculated as follows [21]: where is the salinity of formation water, is the salinity of formation water when , and is the thickness of Na^{+} ion diffusion layer (10^{8} cm).
The pore volume occupied by clay water when is calculated as follows:
The equivalent conductivity of compensation Na^{+} ion in clay water (S/m) (mmol/L) is calculated as follows:
The clay water resistivity is calculated as follows:
It can be seen from the above formula that the resistivity of formation water is affected by both and temperature . When the water salinity is high, . Therefore, the resistivity of clay water is independent from the equilibrium cation concentration and clay types.
3.1.2. Calculation of the Clay Water Saturation
Substitute equation (33) into equation (37),
From the above formula, it can be seen that the clay water saturation increases with the increase of . As the temperature increases, decreases. With the increase of , increases.
3.1.3. Calculation of the Resistivity
3.1.4. Calculation of the Parameter
3.1.5. Calculation of the Normalized Water Saturation
3.1.6. Calculation of the Normalized Water Relative Permeability
In conclusion, the DW relative permeability model can be expressed as follows.
3.2. Sensitivity Analysis of Parameters of DW Model of Shaly Sandstone
3.2.1. The Effect of Irreducible Water Saturation on
Suppose that , , , , , and , oilwater relative permeability curves under different irradiated water saturation are shown in Figure 2. It indicates that the normalized water relative permeability decreases while increases. The reason is that the movable water saturation decreases under the same normalized water saturation as increases, thus the normalized relative permeability of water phase decreases.
3.2.2. The Effect of Residual Oil Saturation on
Suppose that , , , , , and , the oilwater relative permeability curves under different residual oil saturation are shown in Figure 3. It illustrates that the normalized water relative permeability decreases with the increase of because the movable water saturation decreases under the same normalized water saturation when increases, which leads to the decrease of normalized relative permeability of water phase .
3.2.3. The Effect of Cation Exchange Capacity on
Suppose that there are a set of cores with the same parameters as follows. , , , , , , and . Normalized oilwater relative permeability curves under different are shown in Figure 4. It illustrates that the normalized water relative permeability decreases with the increase of .
Studies [22–24] show that clay mineral content is one of the main factors affecting the shape of oilwater relative permeability curve of rock. When water is injected into the core sample, it first enters into larger pores, where the relative permeability of the water phase increases rapidly. Soon after the injection, water gradually enters into small pores, where the flow resistance increases. At the same time, the oil in large pore paths is separated into small oil droplets by the water. If the oil droplets migrate to the vicinity of the pore throat, the socalled “liquid resistance effect” will emerge when the diameter of the oil droplets is similar to that of the pore throat. In this case, the capillary force of the orifice throat must be overcome if the oil droplets want to move [25, 26]. The hydrophilic particles in the pores will move to the pore throat and cause blockage. With the increase of water saturation , the amount of plugging particles will increase, and the relative permeability of water phase will decrease accordingly.
3.2.4. The Effect of Total Salinity on
Suppose that , , , and . Figure 5 shows the normalized oilwater relative permeability curves under different . It illustrates that the normalized water relative permeability increases with increased .
In shaly sand reservoir, with the decrease of salt content in free water, the salinity in free water decreases, clay mineral crystal layer expands, and the formation permeability continues to decline, which is called the saltsensitive phenomenon (Meng, 2012). There are a large number of clay minerals in shaly sand reservoir. Therefore, with the decrease of free water salinity, the salt content in free water decreases, the saltsensitive phenomenon gets worse, and the relative permeability of water phase goes down.
4. Experimental Verification of DW Relative Permeability Model
4.1. Experiments
In order to verify the relationship between resistivity and relative permeability in DW Model, a multifunctional core displacement experiment device was designed (Figure 6). The resistivity and relative permeability of core samples with different water saturations were measured.
(a) Schematic of the apparatus used for the simultaneous measurements
(b) Experimental apparatus
Core samples from Wells A and B are tested to explore the relationship between resistivity and relative permeability. Core samples are divided into two groups. Group I contains less clay and smaller than group II. The basic physical parameters of cores are shown in Table 1. The salinity of formation water in wells A and B are 7500 ppm and 8000 ppm, respectively. The density of brine used is 1.02 g/cm^{3}. The oil viscosity in both wells is 8.8 mPa s at 20°C, and its density is 0.845 g/cm^{3}.

4.2. Verification of Relationship Models between Resistivity and Relative Permeability of Sandstone
The Li model, modified Li model, and Pairoys model are verified with the experimental data of core samples in group I. Figure 7 shows the oilwater relative permeability curves of group I cores. In Figure 7(a), the solid blue triangular and pink dots are the normalized water and oil relative permeability obtained from the water displacing oil experiment, respectively; the brown and sky blue chain dotted lines are the normalized water and oil relative permeability calculated with resistivity based on the Li model; the green and purple broken line are the normalized water and oil relative permeability calculated with resistivity based on Pairoys model; and the blue and red solid line are the ones calculated with resistivity based on the modified Li model. Figure 7(a) indicates that the Li model does not work well in the data process of water flooding. The problems of Pairoys model is that the normalization oil relative permeability it calculates is nonnegligibly larger than the experimental value. However, the normalization water and oil relative permeabilities calculated by the modified Li model are in good agreement with experimental data.
Figure 7(b) shows the oilwater relative permeability curve of cores in group I. The water relative permeabilities calculated by the Li model and Pairoys model are smaller than the experimental data, while the oil relative permeabilities are larger. The oil and water relative permeabilities calculated by the modified Li model fit well with the experimental data in sandstone reservoir with less shale contents.
4.3. Verification of Relationship Models between Resistivity and Relative Permeability of Shaly Sandstone
Experiment is designed to measure the resistivity and relative permeability of shaly sandstone samples in group II. Figure 8 shows the comparison of the experimental results and the model calculated results.
Figure 8(a) shows the normalized relative permeability curves, and Figure 8(b) shows the relative permeability curves. The filled dots are the unsteady oilwater relative permeability experiment data, the broken lines are the relative permeability curves calculated with resistivity by the modified Li model, and the solid lines are calculated with resistivity by the DW model. As shown in Figure 8, the Li model, the water relative permeability calculated by the modified Li model fits well with the experimental data, but the calculated oil relative permeability curve is smaller than the experimental data. Meanwhile, the relative permeability curves of oil and water calculated by the DW relative permeability model better fit the experimental data.
5. Conclusions
This study established the relationship model between the resistivity and oilwater relative permeability of the shaly sandstone reservoir based on the rock physics experiment and the logging response of the shaly sandstone reservoir. According to existing research results, the following main conclusions can be drawn. (1)In view of the influence of shale, the DW relative permeability model, suitable for shaly sandstone reservoir, was derived to calculate oilwater relative permeability using resistivity based on the dual water conductivity model, Poiseuille’s equation, and Darcy’s law.(2)According to the sensitivity analysis, with other conditions being the same, the relative permeability of water phase will decrease as the irreducible water saturation increases, residual oil saturation increases, cation exchange capacity of rock increases, or free water salinity decreases.(3)With the core water flooding experimental device, the resistivity and oilwater relative permeability of two groups of sandstone samples with different shale contents were tested. The experimental results show that the modified Li model is suitable for clean sandstone reservoirs, and the DW relative permeability model is suitable for shaly sandstone reservoirs.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declared that they have no conflicts of interest to this work.
Acknowledgments
The study has been supported by the National Major Science and Technology Projects of China (No. 2017ZX05030002), by the Natural Science Basic Research Plan in Shaanxi Province of China (Grant no. 2020JQ747), and by the Scientific Research Plan Projects of Shaanxi Education Department (Grant no. 18JK0517).
References
 M. T. Szabo, “New methods for measuring imbibition capillary pressure and electrical resistivity curves by centrifuge,” Society of Petroleum Engineers Journal, vol. 14, no. 3, pp. 243–252, 2013. View at: Google Scholar
 G. I. Barenblatt and A. A. Gil'Man, “Nonequilibrium counterflow capillary impregnation,” Journal of Engineering Physics, vol. 52, no. 3, pp. 335–339, 1987. View at: Publisher Site  Google Scholar
 J. Toth, T. Bodi, P. Szucs, and F. Civan, “Direct determination of relative permeability from nonsteadystate constant pressure and rate displacements”. View at: Google Scholar
 E. C. Kumbur, K. V. Sharp, and M. M. Mench, “Validated Leverett approach for multiphase flow in PEFC diffusion media,” Journal of The Electrochemical Society, vol. 154, no. 12, pp. B1295–B1304, 2007. View at: Publisher Site  Google Scholar
 A. T. Corey, “The interrelation between gas and oil relative permeabilities,” Producers Monthly, vol. 19, no. 1, pp. 38–41, 1954. View at: Google Scholar
 R. H. Brooks and A. T. Corey, “Properties of porous media affecting fluid flow,” Journal of the Irrigation and Drainage Division, vol. 92, no. 2, pp. 61–88, 1966. View at: Google Scholar
 J. F. App and J. E. Burger, “Experimental determination of relative permeabilities for a rich gas/condensate system using live fluid,” SPE Reservoir Evaluation & Engineering, vol. 12, no. 2, pp. 263–269, 2013. View at: Google Scholar
 W. Wei, J. C. Cai, X. Y. Hu, and Q. Han, “An electrical conductivity model for fractal porous media,” Geophysical Research Letters, vol. 42, no. 12, pp. 4833–4840, 2015. View at: Publisher Site  Google Scholar
 K. Li, B. Hou, H. Bian, H. Liu, C. Wang, and R. Xie, “Verification of model for calculating capillary pressure from resistivity using experimental data,” Fuel, vol. 252, pp. 281–294, 2019. View at: Publisher Site  Google Scholar
 J. C. Cai, W. Wei, X. Y. Hu, and D. A. Wood, “Electrical conductivity models in saturated porous media: a review,” EarthScience Reviews, vol. 171, pp. 419–433, 2017. View at: Publisher Site  Google Scholar
 J. Cai, Z. Zhang, W. Wei, D. Guo, S. Li, and P. Zhao, “The critical factors for permeabilityformation factor relation in reservoir rocks: porethroat ratio, tortuosity and connectivity,” Energy, vol. 188, p. 116051, 2019. View at: Publisher Site  Google Scholar
 G. E. Archie, “The electrical resistivity log as an aid in determining some reservoir characteristics,” Transaction of American Institute of Mining, Metallurgical, and Petroleum Engineers, vol. 146, no. 2, pp. 54–62, 1942. View at: Google Scholar
 K. Li and R. N. Horne, “Comparison of methods to calculate relative permeability from capillary pressure in consolidated waterwet porous media,” Water resources research, vol. 42, no. 6, pp. W06405.1–W06405.9, 2006. View at: Google Scholar
 K. Li and W. Williams, “Determination of capillary pressure function from resistivity data,” Transport in Porous Media, vol. 67, no. 1, pp. 1–15, 2007. View at: Publisher Site  Google Scholar
 I. A. Mohammed and D. Birol, “An improved relative permeability model to match displacement experiments,” International Journal of Applied Science and Technology, vol. 1, no. 2, pp. 41–58, 2011. View at: Google Scholar
 A. ValdesPerez, H. Pulido, H. CincoLey, and G. GaliciaMuñoz, “Discretization of the resistivity, capillary pressure and relative permeability for naturally fractured reservoirs. Proceedings”. View at: Google Scholar
 F. Pairoys, A. AlZoukani, and A. Keskin, “Interrelationship between resistivity and relative permeability of a carbonate rock during drainage and imbibition experiments”. View at: Google Scholar
 K. Li and R. N. Horne, “A semianalytical method to calculate relative permeability from resistivity well logs”. View at: Google Scholar
 F. Pairoys, “Twophase semidynamic method with electrical measurements: determination of relative permeability and capillary pressure from a single experiment”. View at: Google Scholar
 H. Y. Bian and K. Li, “A modified method and experimental verification for estimating relative permeability from resistivity logging data”. View at: Google Scholar
 H. Q. Xia, X. H. Zhang, and X. J. Liu, “A new approach to calculate clay bound water volume and cation exchange capacity,” Journal of Southwest Petroleum Institute, vol. 22, no. 1, pp. 55–58, 2000. View at: Google Scholar
 R. Zhang, “An approach to the origin of anomalous shape of oilwater relative permeability curves in sandstone and conglomerate reservoirs,” Petroleum Exploration and Development, vol. 23, no. 2, pp. 79–83, 1996. View at: Google Scholar
 H. M. Gao, H. Q. Jiang, M. F. Chen, and Y. F. He, “Microscopic simulation of the effect of microscopic reservoir parameters on oilwater relative permeability,” Petroleum Exploration and Development, vol. 33, no. 6, pp. 734–737, 2007. View at: Google Scholar
 G. M. Linpeng, Z. S. Guohua, and L. F. Lixia, “Feature controls and influence for petroleum development of relative permeability curve in low permeable sandstone reservoirs,” Acta Sedimentologica Sinica, vol. 26, no. 3, pp. 445–451, 2008. View at: Google Scholar
 H. Z. Cui, G. Q. Yao, and F. D. Zhou, “Type and the variety characteristics of relative permeability curve in low permeable glutenite oil reservoirs,” Geological Science and Technology Information, vol. 22, no. 1, pp. 88–91, 2003. View at: Google Scholar
 J. M. He, “Influencing factors of abnormal oil/water relative permeability curve,” Petroleum Geology and Recovery Efficiency, vol. 16, no. 2, pp. 74–76, 2009. View at: Google Scholar
Copyright
Copyright © 2020 Huiyuan Bian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.