ORIGINAL RESEARCH article

Front. Earth Sci., 11 February 2025

Sec. Solid Earth Geophysics

Volume 13 - 2025 | https://doi.org/10.3389/feart.2025.1508283

The relationship between water production rate and resistivity in sand-based porous materials under two-phase fluid flow and its application

  • 1. School of Petroleum, China University of Petroleum (Beijing) at Karamay, Karamay, China

  • 2. State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum (Beijing), Beijing, China

  • 3. China Petroleum Logging Co., Ltd., Xian, Shanxi, China

Abstract

Accurately predicting the water production rate of multiple-phase fluid flows through porous rock is important for many engineering and geological applications. Taking into account irreducible water and capillary tortuosity, the equivalent rock element model from previous studies has been improved. Based on the improved capillary equivalent rock element model, this study proposes a relationship model between the water production rate in two-phase fluid systems and the resistivity index. The relationship model is verified using rock samples, and the result shows that the water production rate calculated by the new model closely matches experimental values, improving parameter calculation accuracy. In addition, applications of the new model to well logging data show that the calculated water production rate is consistent with the actual production situation. This approach is more efficient and accurate in reservoir evaluation. The new model provides a new idea for studying the seepage characteristics of rocks with multiple-phase fluids.

1 Introduction

The seepage characteristics of porous media with multiple-phase fluids are a research hot spot (; ). Understanding the seepage characteristics of rocks in oil–water two-phase flow and oil–gas–water three-phase flow is fundamental for reservoir development planning, production forecasting, and reservoir evaluation (; ; ).

The water production rate, defined as the water volume ratio during reservoir exploitation, is an important seepage characteristic parameter and serves as a dynamic indicator for evaluating the reservoir’s oil-bearing properties. The theoretical results show that there is a clear functional relationship between water production rate and relative permeability (). In order to accurately calculate the water production rate, several researchers have developed methods for relative permeability. Honarpour et al. derived a relationship model between water production rate and water saturation based on the empirical relative permeability model proposed by , , and . This model is one of the most widely used models, but it only considers the general relationship between water production rate and water saturation, and the accuracy is low in low porosity and low permeability reservoirs (; ). analyzed the relationship between water production rate and relative permeability during natural gas exploitation in a fractured reservoir. Clarkson described the variation in the water production rate during coal-bed gas exploitation using relative permeability analysis (). A theoretical model of gas–water relative permeability in gas hydrate reservoirs is established by . These models can accurately predict the water production rate in fractured formations, coal-bed gas formations, and gas hydrate reservoirs, but their application range is limited. Schlachter et al. examined the water production rate using Modular Dynamics Tester (MDT) testing and logging (; ; ). Wang et al. used machine learning and deep learning to predict relative permeability and the oil content ratio during reservoir exploitation, achieving good results (; ).

Previous studies have shown that relative permeability is a key parameter in determining the water production rate, and relative permeability itself is a function of water saturation (; ). Therefore, the relationship between water production rate and water saturation plays a crucial role in determining the water production rate (). Water saturation is a static parameter describing the oil-bearing properties of a reservoir. A common method for determining water saturation is using rock electrical properties, with widely used existing models such as the Archie model and the three-water model (; ; ; ). The equivalent rock element model (EREM), which incorporates the impact of the pore structure on rock resistivity, is also extensively utilized in the investigation of rock seepage (; ; ). Consequently, using rock electrical properties to determine the water production rate proves to be a valuable approach. However, research on this topic is limited.

Research has shown that the seepage characteristics and electrical properties of rocks are closely related to the pore structure (; ; ). proposed a relationship model between relative permeability and rock electrical properties, considering the effect of the pore structure. Xu P et al. developed a relative permeability model based on the fractal theory (; ). The pore structure is an essential factor in establishing the relationship model between water production rate and rock electrical properties. Taking into account the effect of the pore structure, this paper has improved the EREM, and a relationship model between water production rate and resistivity index is proposed. Seven rock samples with different pore structure types were used to validate the new model, and the average relative error between the water production rate calculated by the new model and the experimental data is 5%. The results of both verification and application show that the new model can directly calculate the water production rate from the resistivity index and accurately predict the water production rate of rocks with different pore structures. The novel model incorporates pore structure characterization parameters to mitigate the impact of the pore structure on the water production rate, thereby enhancing the precision of water production rates. In addition, the new model provides a new idea for studying the seepage characteristics of rocks and enriches the data foundation for guiding reservoir development planning, production forecasting, and reservoir evaluation.

2 Relationship between water production rate and electrical resistivity in two-phase fluid flows

2.1 Improved capillary equivalent rock element model

According to the EREM (; ; ), we obtainwhere in Equation 1, (v/v) is the total porosity, (v/v) is the macropore porosity, (v/v) is the throat-pore porosity, and (dimensionless) is the ratio of to . Obviously, is closely related to the pore structure. Figure 1 shows the equivalent rock element model of the porous media. Figure 1A shows the pore distribution of the porous media. Figure 1B shows the equivalent rock element model of the porous media.

FIGURE 1

In Figure 1, is the total porosity, is the macropore porosity, is the throat-pore porosity, and is the matrix volume.

Generally, the flow characteristics are mainly featured by throat pores. Figure 2 shows the flow model of the pore fluids. The equivalent capillary porosity can be expressed aswhere in Equation 2, (v/v) is the equivalent capillary porosity.

FIGURE 2

In Figure 2, is the equivalent capillary porosity and is the matrix volume.

Actually, water-wetted rocks contain irreducible water. The effective flow porosity (v/v) is expressed aswhere (v/v) is the equivalent capillary irreducible water saturation (the ratio of irreducible water porosity to equivalent capillary porosity). (v/v) is the bound water volume. The equivalent capillary irreducible water saturation can be calculated as follows:where (v/v) is the total irreducible saturation (the ratio of total irreducible water to total porosity).

By substituting Equation 4 into Equation 3, the effective porosity of flow fluids is expressed as

The equivalent rock element model assumes that the effective flow pores are straight capillaries. However, studies have shown that the flow of pore fluids in porous media is a complex process (). In order to simulate the actual flow of pore fluids, the effective flow pores are assumed to be curved capillaries. Figure 3 shows the improved capillary equivalent rock element model (ICEREM), where is the tortuosity. The ICEREM introduces the volume of bound water, which is more consistent with the microscopic pore characteristics of the rock, and its pore structure characterization accuracy is further improved than that of the EREM model.

FIGURE 3

In Figure 3, is the effective flow porosity, is the equivalent tortuosity per unit volume of the rock, is the bound water volume, and is the matrix volume.

2.2 Relationship between water production rate and electrical resistivity in two-phase fluid flows

According to Darcy’s law and Poiseuille equation (; ), the water flow through the rock section can be expressed aswhere (m3) is the water flow, (md) is the permeability, (Pa) is the pressure difference in the rock section, (s) is the water viscosity, and (dimensionless) is the equivalent tortuosity per unit volume of the rock.

From Equation 6, it is assumed that fluids flow only in the effective flow pores, as shown in Figure 3.

As the non-wetting phase fluid entered the pores, it was mainly distributed in the central part of the capillary pores. Figure 4 shows the distribution of the two-phase flow fluids in the capillary pores.

FIGURE 4

In Figure 4, is the effective flow porosity, is the equivalent tortuosity per unit volume of the rock, is the bound water volume, is the oil porosity in the effective flow porosity, is the water porosity in the effective flow porosity, is the oil porosity in the macropore porosity, is the water porosity in the macropore porosity, and is the matrix volume.

The water flow and non-wetting phase fluid flow across the rock section can be expressed aswhere (m3) is the water flow, (m3) is the non-wetting phase fluid flow, (s) is the non-wetting phase fluid viscosity, (v/v) is the water saturation of the effective flow porosity, is the equivalent tortuosity of the water flow in the effective flow capillary pores, and (dimensionless) is the equivalent tortuosity of the non-wetting phase fluid flow in the effective flow capillary pores.

According to Equation 7, the water production rate of a two-phase fluid flow can be calculated as follows:

Rock conductivity is the parallel conductance between equivalent capillary porosity and the invalid flow part (invalid flow capillary porosity and the rock matrix). The conductivity of the invalid flow part is assumed to be 0, and the conductivities of water-saturated and partially saturated rock can be obtained as follows:

In Equation 9, (S/M) is the water-saturated rock conductivity, (S/M) is the partially saturated rock conductivity, (S/M) is the water conductivity, (S/M) is the non-wetting phase fluid conductivity, (dimensionless) is the equivalent tortuosity of the irreducible water, and (S/M) is the conductivity of the rock matrix.

The three-water model considers rock conductivity as a parallel conductance among shaly irreducible water, micropore irreducible water, and free water, and each conductance component follows the Archie model (; ). According to the three-water model, the conductivity of water-saturated and partially saturated rock can be further revised as follows:

In Equation 10, both shaly irreducible water and micropore irreducible water are treated as irreducible water and are assumed to have the same conductivity. and (dimensionless) are the cementation exponent and saturation exponent of the equivalent capillary pores, respectively. and represent the cementation exponent and saturation exponent of the irreducible water pores, respectively.

By comparing Equation 9 with Equation 10, we can derive the equivalent tortuosity for the water flow in effective flow capillary pores (dimensionless) and non-wetting phase fluid flow in the same effective flow capillary pores (dimensionless):

In addition, when the non-wetting phase fluid conductivity equals the pore water conductivity, according to Equations 9, 10, can be written as follows:However, the difference between actual non-wetting phase fluid conductivity and wetting phase fluid conductivity may cause a variation in . To eliminate the difference, according to Equations 11, 12, is defined as

By substituting Equations 4, 5, 11, 13 into Equation 8, the water production rate (v/v) of the two-phase fluid flow can be simplified as

According to Equations 9, 10, the relationship between and the resistivity index is as follows:where (dimensionless) is the resistivity index and . Generally, the conductivity of the non-wetting phase fluid (usually oil, gas, or kerosene) is much lower than that of pore water. So, we can assume that . Then, Equation 15 can be simplified aswhere , , and (dimensionless) are the exponents (functions of , , , and ).

By comparing Equation 16 with Equation 14, the relationship between the resistivity index and water production rate can be obtained.

In Equation 17, , , and are introduced to account for the difference between the equivalent tortuosity of pore fluids and that of the electrical current. Generally, , , , , , and are regional coefficients, which can be calculated using the least squares method.

In Equations 16, 17, pore structure characterization parameters and are innovatively introduced, which can accurately characterize the seepage law of the two-phase fluid in rocks with different pore structure characteristics.

3 Verification

To verify the model in this paper, the selected rock samples have the same sedimentary environment and lithology, ensuring that the resistivity and seepage capacity of the samples are mainly affected by the pore structure while minimizing the influence of other factors, such as clay content. Seven rock samples selected from the Dongying Formation, Nanpu Sag, are used to verify the relationship model between the and water production rate. The porosity ranges from 0.19 to 0.255 (v/v), the permeability ranges from 6.09 to 292 (md), and the irreducible water content ranges from 0.518 to 0.861 (v/v). The rock samples are tested using nuclear magnetic resonance (NMR) experiments and combined tests of resistivity and water–oil two-phase seepage. The seven rock samples, with low clay content, are selected to minimize the influence of shale on resistivity. Table 1 presents the characteristic parameters of the seven rock samples.

TABLE 1

NO.Shale content v/vPorosity v/vPermeability mdFormation factorIrreducible water saturation v/vp
10.0170.24146.81112.1310.5440.239
20.0440.25024.4009.0720.5660.075
30.0250.24255.47310.5120.6420.140
40.0290.19445.40015.2740.5180.039
50.0220.255292.0008.0150.5330.095
60.0370.1996.09010.6410.8610.033
70.0380.216102.0006.9710.5360.005
Range
average value
0.017–0.044
0.03
0.19–0.255
0.23
6.09–292
81.7
6.97–15.3
10.37
0.518–0.861
0.6
0.005–0.24
0.089

Characteristic parameters of rock samples.

Figure 5 shows the T2 spectrum of rock samples; the variation in the distribution of the T2 spectra and the position of T2 peak values indicate differences in pore structures. For example, sample No. 1 shows a double-peak right-skew type with the highest peak position of 100 ms; sample No. 2 shows a double-peak balanced type with the highest peak position of 10 ms; sample No. 3 shows a double-peak left-skew type with the highest peak position of 3 ms; sample No. 4 shows a double-peak right-skew type with the highest peak position of 30 ms; sample No. 5 shows a medium peak with the peak position of 10 ms; sample No. 6 shows a single-peak left-skew type with the peak position of 2 ms; and sample No. 7 shows a single-peak left-handed type with the peak position of 3 ms.

FIGURE 5

Figures 68 show the relationship among (resistivity index), (water production rate), and (water saturation). Figure 6 shows the relationship between resistivity and water saturation for the seven rock samples, which exhibit non-Archie phenomenon, and the Archie saturation index n ranges from 1.14–3.1.

FIGURE 6

FIGURE 7

FIGURE 8

It can be observed that data and data for the seven rock samples in Figures 7, 8 are highly dispersed, which is primarily caused by differences in pore structures. This means that significant errors in the water production rate can occur without considering the pore structure (; ; ; ).

In Equation 16, parameter is related to the pore structure and is used to calculate the water production rate. In this section, mercury injection capillary pressure (MICP) and NMR experimental data are used to determine macropores and throat pores.

The Swanson parameter (inflexion point) of the mercury injection curve is used to distinguish large and small pores (throat pores) in rock samples (; ; ). The integral curve of the NMR T2 spectrum corresponds well with the mercury injection curve. According to the definition of the Swanson parameter, T2 time at the inflexion point between macropore and small pores is determined, and then parameter is determined. In this paper, T2 time at the inflexion point is 110 ms, determined by averaging data from 25 rock samples. Figure 9 shows the classification of macropores and throat pores.

FIGURE 9

In Figure 9, the NMR T2 spectra of the core sample (red line) and the pore size distribution of the core sample derived from the MICP experimental data (blue dash-dotted line) are presented. T2 (ms) is the transverse relaxation time. S (T2) (%) is the T2 pore component. ra (μm) is the pore radius. S (ra) (%) is the pore component derived from MICP experimental data. Demarcation T2 time (ms) is the inflexion point between macropores and throat pores, which is 110 ms in this part.

Regional coefficients , , , , , and are calculated using porosity, irreducible water saturation, p, and water production rate through the least squares method: , , , , , and .

In order to compare with previous research results, the water production rate is expressed aswhere the relative permeability of oil and water phases in the Equation 18 can be calculated according to the model proposed by .where in Equation 19, is water saturation, is bound water saturation, is residual oil saturation, is water phase relative permeability, and is oil phase relative permeability. and are regional experience indexes, and and in this region. The calculation results are shown in Figure 10.

FIGURE 10

The water production rate can be calculated using Equation 17 and Wang and Cardenas model. Figure 10 shows the comparison between the calculated water production rate using Equation 17 and experimental data (the red point in Figure 10), where the goodness of fit is 0.97945, and the average relative error is 5%. Figure 10 shows the comparison between the calculated water production rate by the Wang and Cardenas model and experimental data, where the goodness of fit is 0.7995 and the average relative error is 23%. The result indicates a strong agreement between the calculated water production rate and experimental data.

4 Application in well logging data

The new method is applied to Well NPX49 in Nanpu Sag, Huanghua depression, China. The NMR logging measurement mode of Well NPX49 is D9TWE3. Specific parameter areas are as follows: for Group A, TWl = 12.988 s, TEs = 0.9 ms, and the echo number is 500; for Group B, TWs = 1 s, TEs = 0.9 ms, and the echo number is 500; for Group C, TWc = 0.02s, TEc = 0.6 ms, and the echo number is 20; for Group D, TWl = 12.988 s, TEl = 3.6 ms, and the echo number is 125; and for Group E, TWs = 1 s, TEl = 3.6 ms, and the echo number is 125. Drilling mud resistivity is 2.26 OHMM/18°C.

Figure 11 shows the evaluation result of Well NPX49 (layers 1–3) using the new method. The reservoir parameters and water production rate of the produced fluid from layers 1–3 are shown in Table 2.

FIGURE 11

TABLE 2

No.Rt ohm.mPOR %PERM mdSIRR %Sw_Archie %Calculated Fw %Production Fw %
Layer 13921753840.35.50
Layer 22220553665.887.195.6
Layer 32718523361.767.560.4

Reservoir parameters.

In Figure 11, tracks from left to right include tracks 1–4: natural gamma-ray logging (GR: GAPI)/spontaneous potential logging (SP: MV), depth (meters), apparent resistivity logs (RLLD/RLLS: OHMM), acoustic-wave slowness logs (AC: us/m), bulk density (DEN: g/cm3), and neutron porosity (CNL:%). Track 5: NMR logging T2 spectra measured with parameters TE = 0.9 ms and TW=12988 ms (NMR.TA:ms). Track 6: NMR logging T2 spectra measured with parameters TE = 3.6 ms and TW=12988 ms (NMR.TB: ms). Track 7: clay-bound water porosity (PWC: V/V)/ capillary-bound water porosity (PWI: V/V)/ movable water porosity (POR: V/V) computed from NMR logging. Track 8: permeability calculated from NMR logging (NMR.PERM: MD). Track 9: water saturation by the Archie model (ARCHIE.SW: V/V)/ irreducible water saturation from NMR (NMR.SIRR: V/V). Track 9: the new method calculated the water production rate of percolating fluid (FW: V/V). Track 10: the number of the layer.

In layers 1–3, the low natural gamma-ray logging values present the characteristics of sandstone with low clay content. The resistivity of layer 1 is high, with a value of 39 Ω·m. The distribution of NMR T2 spectra for both short-TE and long-TE is broad. The calculated water saturation is 40.3%, and the bound water saturation calculated based on NMR is 38%, demonstrating excellent oil-bearing characteristics. The daily oil production of layer 1 is 14 m³ with no water (the water production rate of layer 1 is 0, as shown in Table 2; represents the volume ratio of water in the produced liquid volume). This indicates that pore water in layer 1 is in a bound state and cannot flow. The resistivity of layers 2 and 3 is low, with values of 22 and 27 Ω·m, respectively. The distribution of NMR T2 spectra for short-TE and long-TE is also broad. The calculated water saturation is 65.8% and 61.7% for layers 2 and 3, respectively, indicating the oil-bearing characteristics of layers 2 and 3.

The distribution of long-TE T2 spectra and short-TE T2 spectra narrows from layer 1 to layer 3. The oil saturations of layers 1–3 are 59.7%, 34.2%, and 38.3%, respectively. The calculated water saturation reduced from layer 1 to layer 3, indicating that layer 1 is more oil-rich than layers 2 and 3. The calculated results of layers 2 and 3 show some differences. However, the daily oil production of layer 2 is 0.2 m3, with a water production rate of 95.6%. The daily oil production of layer 3 is 3.2 m3, with a water production rate of 60.4%. The main reason for the difference between layers 2 and 3 is the difference in the pore structure; for instance, the NMR T2 spectrum in layer 2 shows obvious bimodal right-skew characteristics, while in layer 3, it shows single-modal left-skew characteristics. The application results show that oil saturation does not reflect actual oil production. The water production rate calculated by the new model is consistent with the actual production situation.

5 Discussion and future work

  • (1) Accurately estimating the pore structure is a crucial factor in improving the accuracy of water production rate prediction models. The new model can directly calculate the water production rate from the resistivity index, avoiding errors introduced by water saturation calculation in previous relationship models. On one hand, the new model incorporates pore structure characterization methods, and its accuracy has been verified through experimental data and practical engineering applications. On the other hand, compared with the Wang and Cardenas model, which requires sequential calculations of water saturation, bound water saturation, and residual oil saturation before determining the water production rate, the new model leverages pore structure models and the similarity between current conduction and seepage flow patterns, directly estimating the water production rate based on rock resistivity. This approach is straightforward and minimizes intermediate error propagation.

  • (2) The new model can accurately predict the water production rate of rocks with different pore structures with the introduction of a pore structure-related parameter . The pore structure-related parameter should be obtained from the NMR T2 spectrum. Therefore, resistivity and nuclear magnetic resonance logging data need to be obtained and applied to actual logging data.

  • (3) The error in the water production rate calculated by the new model may be caused by the following factors: ①pore structure differences between the actual model and the ICEREM, ② flow path differences between the electrical and fluid flows, ③measurement errors in experimental data, and ④ the presence of additional conductive minerals in the rock skeleton.

  • (4) Since the ICEREM model mainly considers irreducible water and capillary tortuosity, other pore structure parameters, such as wettability and shaly content, are assumed to be uniform; the accuracy of the new model can be further improved by classifying pore structures and redefining the empirical parameters in Equations 16, 17 of each pore structure type.

  • (5) The wetting phase saturation value in a reservoir does not accurately reflect the wetting phase ratio of the output fluid volume in two-phase fluid seepage, so the calculated water saturation/oil saturation is not consistent with the actual water/oil ratio, as shown in Figure 11. However, the water production rate calculated by the new model is consistent with the actual production situation, indicating that the model can be used to guide reservoir production or the application of two-fluid seepage in a porous material.

Further research directions concerning the electrical properties of sand-based porous materials and the conductivity equation of rocks may include the following aspects:

  • a. A new equation should be established for rock resistivity based on sandstone pore structure characteristics. However, in carbonate, volcanic rocks, and other types of rock, the pore characteristics differ from those of sandstone, such as fractures and caves in carbonate and high- and low-angle fractures in volcanic rocks (). The new model’s adaptability needs further validation, and model parameters must be determined by sample experiments.

  • b. The model in this paper does not consider skeleton conductivity, but some mineral skeletons in rock samples are conductive, such as wet clay and pyrite. Considering the influence of conductive minerals on rock conductivity is the next step to improve the accuracy of the model.

  • c. The conductivity models of rock are widely used in well logging research, reservoir evaluation, reservoir research, and geological research (). One of the next steps is to extend the new model in this paper and apply it to these research areas.

  • d. The conductivity theory of sandstone in this paper can be extended and applied to cement, coal, and other types of porous materials, but it needs to be verified and analyzed by conducting experiments on different types of porous media.

6 Conclusion

  • (1) Pore structure significantly influences rock electrical properties and seepage characteristics. The EREM has been improved by considering irreducible water and capillary tortuosity.

  • (2) Based on the ICEREM, the relationship model between the water production rate and resistivity index has been proposed.

  • (3) Compared to the previous water production rate models ( relationship models), the new model, on one hand, avoids errors caused by water saturation calculation; on the other hand, it introduces a pore structure-related parameter to address the effect of the pore structure on rock seepage characteristics and electrical properties.

  • (4) The results of verification and application have shown that the water production rate calculated by the new model is in accordance with experimental values and the actual water production rate.

  • (5) The new model can predict the water production rate based on electrical properties, dynamically enhance reservoir evaluation, and strengthen the database for reservoir development planning, production forecasting, and reservoir evaluation.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

WX: conceptualization, methodology, writing–original draft, and writing–review and editing. QY: conceptualization, methodology, and writing–review and editing. XD: writing–review and editing. YF: data curation, validation, and writing–review and editing. PZ: writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This article was supported by the Research Foundation of Karamay, China (2024hjcxrc0075), the Research Foundation of the China University of Petroleum-Beijing at Karamay (No. XQZX20230012), “Tianchi Talent” Introduction Plan Foundation of Xinjiang, and the Natural Science Foundation of Xinjiang (No. 2021D01E22).

Conflict of interest

Author YF was employed by China Petroleum Logging Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    ArchieG. E. (1942). The electrical resistivity log as an aid in determining some reservoir characteristics. Trans. AIME146 (01), 5462. 10.2118/942054-g

  • 2

    BearJ. (2013). Dynamics of fluids in porous media. Cour. Corp.

  • 3

    BehrangA.MohammadmoradiP.TaheriS.KantzasA. (2016). A theoretical study on the permeability of tight media; effects of slippage and condensation. Fuel181, 610617. 10.1016/j.fuel.2016.05.048

  • 4

    BurdineN. T. (1953). Relative permeability calculations from pore size distribution data. Petr.trans.am Inst.mining Metall.eng.98 (3), 7178.

  • 5

    CaiJ.LinD.SinghH.ZhouS.MengQ.ZhangQ. (2019). A simple permeability model for shale gas and key insights on relative importance of various transport mechanisms. Fuel252, 210219. 10.1016/j.fuel.2019.04.054

  • 6

    ChenX.YaoG. (2017). An improved model for permeability estimation in low permeable porous media based on fractal geometry and modified Hagen-Poiseuille flow. Fuel210, 748757. 10.1016/j.fuel.2017.08.101

  • 7

    ChimáA.ChávezE. A.CalderónZ. (2010). “An equation to predict two-phase relative permeability curves in fractures,” in SPE Latin American and Caribbean petroleum engineering conference (Society of Petroleum Engineers). 10.2118/138282-MS

  • 8

    ClarksonC. R.RahmanianM.KantzasA.MoradK. (2011). Relative permeability of CBM reservoirs: controls on curve shape. Int. J. Coal Geol.88 (4), 204217. 10.1016/j.coal.2011.10.003

  • 9

    CoreyA. T. (1954). The interrelation between gas and oil relative permeabilities. Prod. Mon.19 (1), 3841.

  • 10

    DarcyH. (1856). Les Fontaines Publiques de la Ville de Dijon. Paris: Dalmont.

  • 11

    HonarpourM. M. (2018). Relative permeability of petroleum reservoirs. CRC Press.

  • 12

    IscanA. G. (2021). Water saturation calculation using fractional flow and production logging data in a Caspian region sandstone petroleum reservoir. J. Petroleum Sci. Eng.200, 108355. 10.1016/j.petrol.2021.108355

  • 13

    JinY.LiS.YangD. (2020). Experimental and theoretical quantification of the relationship between electrical resistivity and hydrate saturation in porous media. Fuel269, 117378378. 10.1016/j.fuel.2020.117378

  • 14

    KashaA.MyersM.HathonL.Sakhaee-PourA.SadooniF.NasserM.et al (2023). Integrated approach for closure correction of mercury injection capillary pressure measurements. Geoenergy Sci. Eng.230, 212245. 10.1016/j.geoen.2023.212245

  • 15

    LeiG.LiaoQ.LinQ.ZhangL.XueL.ChenW. (2020). Stress dependent gas-water relative permeability in gas hydrates: a theoretical model. Adv. Geo-Energy Res.4 (3), 326338. 10.46690/ager.2020.03.10

  • 16

    McKennaS. A.AkhrievA.CiaurriD. E.ZhukS. (2020). Efficient uncertainty quantification of reservoir properties for parameter estimation and production forecasting. Math. Geosci.52 (2), 233251. 10.1007/s11004-019-09810-y

  • 17

    MoX. W.HeD. H.LiZ.BoWenX. C.LiG. R. (2001). The application of three-water conduction model in the interpretation of low-resistivity reservoir. J. Changchun Univ. ence Technol (01), 9295. 10.13278/j.cnki.jjuese.2001.01.019

  • 18

    PengXuQiuS.YuB.JiangZ. (2013). Prediction of relative permeability in unsaturated porous media with a fractal approach. Int. J. Heat Mass Transf.64, 829837. 10.1016/j.ijheatmasstransfer.2013.05.003

  • 19

    PiusT. O.OlamigokeO. (2020). “Investigating the correlation between water saturation obtained from cased-hole saturation tool measurements and produced water cut in strong water drive reservoirs,” in SPE Nigeria annual international conference and exhibition, D013S018R002.

  • 20

    PurcellW. R. (1949). Capillary pressures - their measurement using mercury and the calculation of permeability therefromJ. Pet. Technol.1 (02), 3948.

  • 21

    RangelovM.NassiriS. (2018). Empirical time-dependent tortuosity relations for hydrating mortar mixtures based on modified archie's law. Constr. Build. Mater.171, 825838. 10.1016/j.conbuildmat.2018.03.173

  • 22

    Reza-E-RabbiS.AhmmedS. F.ArifuzzamanS. M.SarkarT.KhanM. S. (2020). Computational modelling of multiphase fluid flow behaviour over a stretching sheet in the presence of nanoparticles. Eng. Sci. Technol. Int. J.23 (3), 605617. 10.1016/j.jestch.2019.07.006

  • 23

    RutqvistJ.WuY. S.TsangC. F.BodvarssonG. (2002). A modeling approach for analysis of coupled multiphase fluid flow, heat transfer, and deformation in fractured porous rock. Int. J. Rock Mech. Min. Sci.39 (4), 429442. 10.1016/s1365-1609(02)00022-9

  • 24

    SchlachterG. (2007). “Using wireline formation evaluation tools to characterize coalbed methane formations,” in Eastern regional meeting (Society of Petroleum Engineers). 10.2118/111213-ms

  • 25

    SchönJ. (2011). Physical properties of rocks: a workbook, 8. Elsevier.

  • 26

    ShangB. Z.HammanJ. G.CaldwellD. H. (2003). “A physical model to explain the first Archie relationship and beyond,” in SPE annual technical conference and exhibition (Society of Petroleum Engineers). 10.2118/84300-MS

  • 27

    ShangB. Z.HammanJ. G.CaldwellD. H. (2004). “Water saturation estimation using equivalent rock element model,” in SPE annual technical conference and exhibition (Society of Petroleum Engineers). 10.2118/90143-MS

  • 28

    ShangB. Z.HammanJ. G.CaldwellD. H. (2008). “Improved water saturation estimation using equivalent rock element model and applications to different rock types,” in Europec/EAGE conference and exhibition (Society of Petroleum Engineers). 10.2118/113342-MS

  • 29

    SmirnovN. N.NikitinV. F.SkrylevaE. I. (2019). Microgravity investigation of seepage flows in porous media. Microgravity Sci. Technol.31 (1), 629639. 10.1007/s12217-019-09733-7

  • 30

    SongR.WangY.LiuJ.CuiM.LeiY. (2019). Comparative analysis on pore‐scale permeability prediction on micro‐ct images of rock using numerical and empirical approaches. Energy ence and Eng.7 (2), 28422854. 10.1002/ese3.465

  • 31

    SwansonB. F. (1981). A simple correlation between permeabilities and mercury capillary pressures. J. Petroleum Technol.33 (12), 24982504. 10.2118/8234-pa

  • 32

    TiabD.DonaldsonE. C. (2024). Petrophysics: theory and practice of measuring reservoir rock and fluid transport properties. Elsevier.

  • 33

    WangF.JiaoL.ZhaoJ.CaiJ. (2019). A more generalized model for relative permeability prediction in unsaturated fractal porous media. J. Nat. Gas Sci. Eng.67, 8292. 10.1016/j.jngse.2019.04.019

  • 34

    WangL.CardenasM. B. (2018). Connecting pressure-saturation and relative permeability models to fracture properties: the case of capillary-dominated flow of supercritical CO2 and brine. Water Resour. Res.54 (9), 69656982. 10.1029/2018wr023526

  • 35

    WangY.LiH.XuJ.LiuS.WangX. (2022). Machine learning assisted relative permeability upscaling for uncertainty quantification. Energy245, 123284. 10.1016/j.energy.2022.123284

  • 36

    WhittleT. M.LeeJ.GringartenA. C. (2003). “Will wireline formation tests replace well tests?,” in SPE annual technical conference and exhibition (Society of Petroleum Engineers). 10.2118/84086-MS

  • 37

    WuY.ZhangH.YangD.ZhangH.YuanW. (2018). Oil-gas breakthrough time and water content predict or mdt pump in low permeability reservoir. Well logging Technol.42 (5), 557561. 10.16489/j.issn.1004-1338.2018.05.012

  • 38

    XiaoL.LiuD.WangH.LiJ.LuJ.ZouC. (2017). The applicability analysis of models for permeability prediction using mercury injection capillary pressure (MICP) data. J. Petroleum Sci. Eng.156, 589593. 10.1016/j.petrol.2017.06.042

  • 39

    XieW.YinQ.GuanW.WangG.LaiJ. (2020). Estimating the relative permeability from the electrical parameters of sandstone with a complex pore structure. Energy&Fuels34, 1412414131. 10.1021/acs.energyfuels.0c03218

  • 40

    XieW.YinQ.WuL.YangF.ZhaoJ.WangG. (2024b). A new nuclear magnetic resonance-based permeability model based on two pore structure characterization methods for complex pore structure rocks: permeability assessment in Nanpu Sag, China. Geophysics89 (1), MR43MR51. 10.1190/geo2023-0026.1

  • 41

    XieW.YinQ.ZengJ.WangG.FengC.ZhangP. (2023). Fractal-based approaches to pore structure investigation and water saturation prediction from NMR measurements: a case study of the gas-bearing tight sandstone reservoir in Nanpu sag. Fractal Fract.7, 273. 10.3390/fractalfract7030273

  • 42

    XieW.YinQ.ZengJ.YangF.ZhangP.YanB. (2024a). An improved rock resistivity model based on multi-fractal characterization method for sandstone micro-pore structure using capillary pressure. Fractal Fract.8, 118. 10.3390/fractalfract8020118

  • 43

    XuF.MuL.WuX.SunT.DingY.TianX.et al (2014). New expression of oil/water relative permeability ratio vs. water saturation and its application in water flooding curve. Energy Explor. and exploitation32 (5), 817830. 10.1260/0144-5987.32.5.817

  • 44

    YinX.JiangS.YanLuL.WeiG.JungangL.PengW.et al (2020). Impact of pore structure and clay content on the water-gas relative permeability curve within tight sandstones: a case study from the ls block, eastern ordos basin, China. J. Nat. Gas ence Eng.81, 103418. 10.1016/j.jngse.2020.103418

  • 45

    Yushu Wu (2015). Multiphase fluid flow in porous and fractured reservoirs. Houston: Gulf professional publishing.

  • 46

    ZhangL.PanB.LiZ.MoX.XiaZ.XuW. (2010). New three-water conduction model and its application in evaluation of low porosity and low permeability reservoir. Oil Geophys. Prospect.45 (3), 431435. 10.13810/j.cnki.issn.1000-7210.2010.03.025

  • 47

    ZhangZ.WellerA. (2021). A comparative study of permeability prediction for Eocene sandstones—Part 1: application of modified Swanson models to mercury injection capillary pressure and nuclear magnetic resonance data. Geophysics86 (6), M233M243. 10.1190/geo2021-0194.1

  • 48

    ZhongZ.SunA. Y.RenB.WangY. (2021). A deep-learning-based approach for reservoir production forecast under uncertainty. SPE J.26 (03), 13141340. 10.2118/205000-pa

Summary

Keywords

pore structure, resistivity index, water production rate, two-phase fluid, sand-based porous material

Citation

Xie W, Yin Q, Dai X, Fan Y and Zhang P (2025) The relationship between water production rate and resistivity in sand-based porous materials under two-phase fluid flow and its application. Front. Earth Sci. 13:1508283. doi: 10.3389/feart.2025.1508283

Received

09 October 2024

Accepted

09 January 2025

Published

11 February 2025

Volume

13 - 2025

Edited by

Weichao Yan, Ocean University of China, China

Reviewed by

Suhaib Umer Ilyas, Jeddah University, Saudi Arabia

Li Dan, CNOOC Research Institute Ltd, China

Wang Song, China University of Petroleum, China

Updates

Copyright

*Correspondence: Weibiao Xie, ; Qiuli Yin,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics