# A Generalized Capillary Imbibition Model for Porous Media in Tight Reservoirs

Authors: Zhiyuan Wang; Zhengming Yang; Yunhong Ding; Wei Lin; Ying He; Xiaolang Duan
Journal: Advances in Civil Engineering
 (2018)

---

## Abstract

Capillary imbibition models have been widely studied in oil and gas development field over the past decades. However, the existing models applied to the tight reservoirs rarely take fluid flow resistance and apparent viscosity into account. To investigate the capillary imbibition characteristics of fluids in tight porous media, a generalized capillary imbibition model considering the flow resistance and apparent viscosity of fluids in tight porous media is derived. By comparing with the results of other capillary imbibition models and experimental data, the derived capillary imbibition model is verified. In addition, compared with the conventional capillary imbibition models, the derived capillary imbibition model is more consistent with the experimental results and has a wider applicability. The imbibition distance of water in tight reservoirs can also be obtained using the derived capillary imbibition model, which will facilitate the study on water injection development in tight oil and gas reservoirs.


---

## Body

## 1. Introduction

With the application of large-scale volume fracturing technology in the development of tight oil and gas, the imbibition is playing an increasingly important role in enhancing the recovery of tight reservoirs [1–5]. As the capillary imbibition of liquids in porous media has been utilized extensively in reservoir development, paper treatment, and antiseepage of the dam and other fields [6–11], the subject on imbibition had received wide attention and research from the last century. Bell and Cameron [12], Lucas [13], and Washburn [14] have put forward capillary imbibition equations, based on Hagen–Poseuille (H-P) flow in a cylindrical tube model, which laid the foundation for the development of imbibition model. Due to the complexity of natural porous media, the classical LW equation has been modified by considering the different structural characteristics of porous media [15–21].Dullien et al. [15] derived the equation of capillary rising rate in sandstone samples, based on three-dimensional straight pore structure network model which is composed of a repeat capillary element with a step change. Hammecker et al. [16] proposed a new cylindrical model to describe capillary imbibition in porous sedimentary rocks, consisting of a periodic succession of a single hollow spherical element. Benavente et al. [17] modified the LW equation by introducing the correction factors of the pore shape and tortuosity to describe the pore space of natural rocks. Fries and Dreyer [18] derived an analytic solution for the capillary rise of liquids in a cylindrical tube or a porous medium in terms of height h as a function of time t. Cai [19–21] introduced the fractal theory into the capillary imbibition model, and he derived a new generalized capillary imbibition model through the improved H-P equation and the modified Laplace–Young (L-Y) equation, which takes into account the size and shape of the pores, the tortuosity of random porous media, and the initial wetting-phase saturation [8]. However, these models cannot be well applied to the capillary imbibition of tight rocks, because tight rock is a kind of complex porous media, which has fine pore-throat structure and obvious fluid-solid effect.So far, it is generally believed that Darcy flow only occurs under certain conditions, and non-Darcy flow (nonlinear percolation) is more common [22]. For example, Darcy’s law cannot accurately describe the porous flow of high viscosity crude oil in the reservoir. And for the flow in low-permeability rocks, especially for that in tight rocks, the Darcy’s law is not applicable as well, due to the complex pore-throat structure and strong fluid-solid action of tight reservoir. In terms of non-Darcy flow, scholars have studied it by means of direct experiments, numerical derivation, or combined conceptual models. And a series of results have been obtained in the study of nonlinear percolation mechanism [22–28]. Miller and Low found the presence of a threshold gradient for water flow in clays and proved that no flow occurs when the hydraulic gradient is below the threshold gradient [23]. Huang et al. expounded the characteristics of nonlinear percolation in low-permeability porous media according to experimental data. In his opinion, the porous media have selectivity for the passage of fluids, and the apparent viscosity of the same fluid in different sizes of pores is different, which makes the seepage law of low-permeability reservoirs obey nonlinear law with a start-up pressure gradient. Therefore, the characteristics of the nonlinear seepage cannot be ignored in the study of imbibition of tight cores [22].In this paper, a new generalized capillary imbibition model is derived by considering the resistance of fluids in tight porous media and the change of apparent viscosity in tight pore space on the basis of previous studies.


## 2. Characteristics of Porous Flow in Tight Cores

Tight core is defined as the core collected from the tight reservoir, whose air permeability is less than 1 × 10−3 μm2, and the pore-throat radius of tight core is mainly distributed in the range from 10 nm to 1000 nm [29]. We assume the core as a capillary model with the same average pore radius.


### 2.1. Flow Resistance of Fluids in Tight Cores

A large number of experiments [30, 31] indicate that the fluid flow in tight cores follows the nonlinear law with a start-up pressure gradient (Ga), as shown in Figure 1. Figure 1 illustrates that the fluid does not flow when the pressure gradient is less than a certain value, and the pressure gradient (Ga) is the minimum resistance gradient that the fluid needs to overcome when flowing in the tight core. When the pressure gradient is greater than Ga, the curve in Figure 1 is a concave trend, because the fluid begins to overcome fluid-solid force to move. When the pressure gradient is greater than a certain value (Gb), the curve in Figure 1 becomes a straight line, this is because most of the fluid in pores is fully utilized and start flowing. The pressure gradient (Gb) is called the maximum resistance gradient. Hence, the fluid flow resistance in porous media should be fully considered when it comes to seepage of fluids in tight reservoirs.Figure 1
Classical curve and diagram of nonlinear porous flow in tight cores under different pressure gradients.


### 2.2. Apparent Viscosity of Fluids in Tight Cores

According to the study by Huang [22], the fluid is divided into two parts: boundary fluid and bulk fluid. Moreover, the boundary fluid takes up too large space of tight cores to ignore. The apparent viscosity of fluid in tight core is available under the expression(1)μc=Aμ1+1−Aμ=μ1μ−1A+1μ=ξμ,where μc is the apparent viscosity of fluid in tight core (MPa·s), A is the ratio of volume of boundary fluid and volume of total fluid in tight core, μ1 is the average viscosity of boundary fluid in tight core (MPa·s), μ is the viscosity of bulk fluid (MPa·s), and ξ is the apparent viscosity coefficient.


## 2.1. Flow Resistance of Fluids in Tight Cores

A large number of experiments [30, 31] indicate that the fluid flow in tight cores follows the nonlinear law with a start-up pressure gradient (Ga), as shown in Figure 1. Figure 1 illustrates that the fluid does not flow when the pressure gradient is less than a certain value, and the pressure gradient (Ga) is the minimum resistance gradient that the fluid needs to overcome when flowing in the tight core. When the pressure gradient is greater than Ga, the curve in Figure 1 is a concave trend, because the fluid begins to overcome fluid-solid force to move. When the pressure gradient is greater than a certain value (Gb), the curve in Figure 1 becomes a straight line, this is because most of the fluid in pores is fully utilized and start flowing. The pressure gradient (Gb) is called the maximum resistance gradient. Hence, the fluid flow resistance in porous media should be fully considered when it comes to seepage of fluids in tight reservoirs.Figure 1
Classical curve and diagram of nonlinear porous flow in tight cores under different pressure gradients.


## 2.2. Apparent Viscosity of Fluids in Tight Cores

According to the study by Huang [22], the fluid is divided into two parts: boundary fluid and bulk fluid. Moreover, the boundary fluid takes up too large space of tight cores to ignore. The apparent viscosity of fluid in tight core is available under the expression(1)μc=Aμ1+1−Aμ=μ1μ−1A+1μ=ξμ,where μc is the apparent viscosity of fluid in tight core (MPa·s), A is the ratio of volume of boundary fluid and volume of total fluid in tight core, μ1 is the average viscosity of boundary fluid in tight core (MPa·s), μ is the viscosity of bulk fluid (MPa·s), and ξ is the apparent viscosity coefficient.


## 3. Mathematical Model

### 3.1. Modified H-P and L-Y Equations

According to the study by Cai [8], we can get the modified H-P and L-Y equations, respectively, for tortuous capillaries with variably shaped apertures:(2)q=παrh4ΔP8μτL0=πr4ΔP8μτL0,Pc=2σcosθαrh=2σcosθr,where L0 is the length of the capillary tube, m, τ is the tortuosity, ΔP is the drop pressure along flow direction, MPa, rh = 2AC/C, AC is the capillary’s cross-sectional area, C is the capillary’s circumference of cross-sectional, r is the equivalent pore radius, μm, α is the dimensionless geometry correction factor, Pc is the capillary pressure, MPa, σ is the surface tension of the liquid-gas interface, N/m, and θ is the contact angle.We assume that a tight core hasn pores, and we can get a new modified H-P and L-Y equations for apparent viscosity of fluid in tight cores:(3)qa=∑nqin=παra4ΔP8μcτL0=πra4ΔPξ8μτL0,Pca=∑nPCn=2σcosθra,where ra is the average equivalent pore radius of the cores, Pca is the average capillary pressure, MPa, and qa is the average flow rate, m3/d.According to the simultaneous connection of H-P and Darcy equations, the average equivalent pore radius can be written as(4)ra=8Kϕ.


### 3.2. Capillary Imbibition Model for Tight Porous Media

We assume that the tight core is equivalent to a cylinder consisting ofn tortuous, irregular cross-sectional capillary tubes. When the wetting liquid imbibes in the pore of the tight core, it will be subjected to capillary forces (Pca), the hydrostatic pressure (Pha), and the resistance of porous flow (PGa). We can calculate the total pressure drop in the following equation:(5)ΔP=Pca+Pha+PGa,where Pha=ρgL0sinψ and PGa=GL0.Inserting Equation (5) into Equation (3):(6)qa=πra4ΔP8μcτL0tPca+Phat+PGat,vt=qπαra2=πra2ΔP8μcτL0t2σcosθra−ρg sinψL0t−GτL0t.According to equationvt=dLt/d, we can get the relation equation of the imbibition distance and imbibition time:(7)dLatdt=ra28μcτL0t2σcosθra−ρg sinψL0t−GτL0t,dL0tdt=raσcosθ4μcτ21L0t−ra2ρg sinψ+τG8μcτ2.Rearranging Equation (7) gives(8)dL0tdt=aL0t−b,where(9)a=raσcosθ4μcτ2,b=ra2ρg sinψ+τG8μcτ2.The relationship oft-h can be obtained through the integral of Equation (8):(10)t=−L0tb−ab2ln1−bL0ta.Due to the condition “1−bL0t/a>0,” we can know the maximum distance of imbibition is Lmax=a/b.We can assume that whenL0t=0.99Lmax, the total time of imbibition tmax can be calculated through substituting into Equation (10):(11)tmax=ab2ln−0.99−ln0.01≈3.62ab2.The LambertW is introduced to obtain an explicit analytical solution for Equation (7). The function W(x) is defined by an inverse exponential function [18]:(12)x=WxeWx,where(13)Wx≈2ex−10.7036+7.568592+2ex12.7036+5.135012+2ex−e−1≤x≤0.Therefore, the explicit analytical solution of Equation (8) can be obtained:(14)L0t=ab1+W−e−1−b2/at,where e is Euler’s number (2.718282…).The accumulated weight of a wetting liquid imbibing into a tight core can be given by(15)Mt=ρAϕSwf−SwiL0t=abρAϕSwf−Swi1+W−e−1−b2/at,where ρ is the density of fluid, kg/m3, A is the cross-sectional area of the cores, m2, ϕ is the porosity, Swf is the saturation of wetting liquid after imbibition, and Swi is the initial saturation of wetting liquid.


## 3.1. Modified H-P and L-Y Equations

According to the study by Cai [8], we can get the modified H-P and L-Y equations, respectively, for tortuous capillaries with variably shaped apertures:(2)q=παrh4ΔP8μτL0=πr4ΔP8μτL0,Pc=2σcosθαrh=2σcosθr,where L0 is the length of the capillary tube, m, τ is the tortuosity, ΔP is the drop pressure along flow direction, MPa, rh = 2AC/C, AC is the capillary’s cross-sectional area, C is the capillary’s circumference of cross-sectional, r is the equivalent pore radius, μm, α is the dimensionless geometry correction factor, Pc is the capillary pressure, MPa, σ is the surface tension of the liquid-gas interface, N/m, and θ is the contact angle.We assume that a tight core hasn pores, and we can get a new modified H-P and L-Y equations for apparent viscosity of fluid in tight cores:(3)qa=∑nqin=παra4ΔP8μcτL0=πra4ΔPξ8μτL0,Pca=∑nPCn=2σcosθra,where ra is the average equivalent pore radius of the cores, Pca is the average capillary pressure, MPa, and qa is the average flow rate, m3/d.According to the simultaneous connection of H-P and Darcy equations, the average equivalent pore radius can be written as(4)ra=8Kϕ.


## 3.2. Capillary Imbibition Model for Tight Porous Media

We assume that the tight core is equivalent to a cylinder consisting ofn tortuous, irregular cross-sectional capillary tubes. When the wetting liquid imbibes in the pore of the tight core, it will be subjected to capillary forces (Pca), the hydrostatic pressure (Pha), and the resistance of porous flow (PGa). We can calculate the total pressure drop in the following equation:(5)ΔP=Pca+Pha+PGa,where Pha=ρgL0sinψ and PGa=GL0.Inserting Equation (5) into Equation (3):(6)qa=πra4ΔP8μcτL0tPca+Phat+PGat,vt=qπαra2=πra2ΔP8μcτL0t2σcosθra−ρg sinψL0t−GτL0t.According to equationvt=dLt/d, we can get the relation equation of the imbibition distance and imbibition time:(7)dLatdt=ra28μcτL0t2σcosθra−ρg sinψL0t−GτL0t,dL0tdt=raσcosθ4μcτ21L0t−ra2ρg sinψ+τG8μcτ2.Rearranging Equation (7) gives(8)dL0tdt=aL0t−b,where(9)a=raσcosθ4μcτ2,b=ra2ρg sinψ+τG8μcτ2.The relationship oft-h can be obtained through the integral of Equation (8):(10)t=−L0tb−ab2ln1−bL0ta.Due to the condition “1−bL0t/a>0,” we can know the maximum distance of imbibition is Lmax=a/b.We can assume that whenL0t=0.99Lmax, the total time of imbibition tmax can be calculated through substituting into Equation (10):(11)tmax=ab2ln−0.99−ln0.01≈3.62ab2.The LambertW is introduced to obtain an explicit analytical solution for Equation (7). The function W(x) is defined by an inverse exponential function [18]:(12)x=WxeWx,where(13)Wx≈2ex−10.7036+7.568592+2ex12.7036+5.135012+2ex−e−1≤x≤0.Therefore, the explicit analytical solution of Equation (8) can be obtained:(14)L0t=ab1+W−e−1−b2/at,where e is Euler’s number (2.718282…).The accumulated weight of a wetting liquid imbibing into a tight core can be given by(15)Mt=ρAϕSwf−SwiL0t=abρAϕSwf−Swi1+W−e−1−b2/at,where ρ is the density of fluid, kg/m3, A is the cross-sectional area of the cores, m2, ϕ is the porosity, Swf is the saturation of wetting liquid after imbibition, and Swi is the initial saturation of wetting liquid.


## 4. Results and Discussion

### 4.1. Comparison and Analysis

Through research of the existing imbibition models, it is found that Cai’s model, which takes into account the size and shape of pores, the tortuosity of imbibition streamlines in random porous media, and the initial wetting-phase saturation, is generalized and can be simplified to other models, such as Fries–Dreyer model [8]. However, because the resistance of imbibition is little in conventional porous media with relatively large pores, Cai’s model omits the resistance term. As we all know, if there is no horizontal imbibition resistance in tight reservoirs, the distance of imbibition will be infinite, which is obviously unreasonable.Therefore, based on Cai’s model [8], our model considers the resistance of imbibition and the apparent viscosity of fluid in tight cores in order to apply in tight reservoirs. Fries–Dreyer model [18], Cai’s model, and our model can be formed into “dL/dt=a/L−b” mode. Table 1 shows the values of “a and b” in the mode.Table 1
Values of “a and b” in “dL/dt=a/L−b” mode.
Imbibition model
Direction
a
b
Washburn model
Vertical
r
σ
cos
θ
/
2
μ
0
Horizontal
r
σ
cos
θ
/
2
μ
0
Fries–Dreyer model
Vertical
r
σ
cos
θ
/
4
μ
r
2
ρ
g
/
8
μ
Horizontal
r
σ
cos
θ
/
4
μ
0
Cai’s model
Vertical
r
a
σ
cos
θ
/
4
μ
τ
2
r
a
2
ρ
g
/
8
μ
τ
2
Horizontal
r
a
σ
cos
θ
/
4
μ
τ
2
0
Our model
Vertical
r
a
σ
cos
θ
/
4
μ
c
τ
2
r
a
2
ρ
g
+
τ
G
/
8
μ
c
τ
2
Horizontal
r
a
σ
cos
θ
/
4
μ
c
τ
2
r
a
2
G
/
8
μ
c
τGenerally speaking, our model takes more geometrical parameters into account than other imbibition models in order to better match the real process of imbibition in tight cores. Accordingly, our model can derive other existing imbibition models under conditions ofG = 0, μc = μ as well.We choose tight sandstone data (ra = 1 × 10−7 m, σ  = 0.074 N·m, θ  = 80°, μc = 1.08 cp, τ  = 2, G = 0.43 MPa/s, ρ  = 1000 g/m3, and g = 9.8 m/s2) for calculation in these imbibition models at horizontal direction to quantitatively compare our model with the others.As shown in Figure2, since the Washburn Model and Fries–Dreyer Model do not take into account the tortuosity, the calculated value has a large gap with Cai’s and our model’s calculated value. The bigger the value of tortuosity is, the greater the difference is. Cai’s model fits well in the initial stage of our model. However, since Cai’s model does not consider the influence of resistance, there is no limit as the imbibition distance increases with imbibition time. As the imbibition continues, our model calculated curve tends to be flat, marking the end of the process of imbibition. Therefore, our model can be used for the calculation of the whole process of imbibition, especially in tight sandstone conditions.Figure 2
Quantitative comparison of our model to other models with tight sandstone data.


### 4.2. Model Validation

WhenG = 0 and μc = μ, our model is same as Cai’s model, which has been verified through porous materials with porosity of 0.045–0.7 and permeability of 5–1200 mD.The experimental data [32] of horizontal imbibition in tight cores (k < 1 mD) are collected to verify these models, and petrophysical parameters of tight cores for imbibition experiment are shown in Table 2.Table 2
Petrophysical parameters of tight cores for imbibition experiment.
Number
Diameter (cm)
Length (cm)
Porosity (%)
Permeability (mD)
Average radius (m)
Density of fluid (kg/m3)
Viscosity of fluid (MPa·s)
Exp1
2.53
3.118
9.48
0.14
0.118
1000
1.183
Exp2
2.52
2.794
5.85
0.215
0.294
1000
1.183
Exp3
2.51
3.000
6.89
0.0229
0.026
1000
1.183It can be seen from Figure3 that the result of Fries–Dreyer model has significant difference with the experimental data, and there is a good agreement between the result of Cai’s model and the experimental data in the early time when τ = 8, but Cai’s model cannot fit the experimental data at the later stage very well.Figure 3
Variation of imbibition distance with imbibition time in tight cores: simulation results of Fries–Dreyer Model and Cai’s Model with different tortuosity vs. experimental data of Exp1.In our model, the comprehensive coefficient (ζ) is defined to be tortuosity (τ) multiplied with apparent viscosity coefficient (ξ), which can characterize more influence parameters. Figure 4 shows that the result of our model (ζ = 7.5 and G = 0.4 MPa/m) can fit the experimental data better than other models. We also verify our model by different experimental data, as shown in Figure 5. It can be known that the resistance (G) is positively correlated with average radius, and the comprehensive coefficient (ζ) is negatively correlated with average radius.Figure 4
Variation of imbibition distance with imbibition time in tight cores: simulation results of our model with different composite factors (ζ) and fluid flow resistance (G) vs. experimental data of Exp1.Figure 5
Comparison between imbibition experimental data of different tight cores and calculated results of our model under different composite factors (ζ) and fluid flow resistance (G).As shown in Figure6, it can be concluded that the resistance (G) reduces the rate of imbibition, and the influence of resistance (G) increases with the increase of imbibition time. Therefore, our model can fit the result of experiment better. According to Equation (14) and Exp3, we can predict that the total time of imbibition is 56 days and the maximum distance of imbibition is 0.09 m in corresponding tight reservoirs.Figure 6
The variation curves of imbibition rate with imbibition time under different fluid flow resistances (G) in tight porous media.Due to a large amount of experimental data [33, 34], the empirical equation of imbibition resistance (G) and average pores radius (ra) can be written by(16)G=MeNra,where M and N are the empirical constants. According to experimental data of Exp1, 2, and 3 and their fitting data with our model, we can obtain the M = 3.2347 and N = −2 × 107.When the imbibition resistance (G) is not considered, the smaller the average pore radius (r) is, the slower the imbibition will be, and the imbibition distance is unlimited. When the imbibition resistance (G) is considered, as shown in Figure 7(a), in the early stage of imbibition, the smaller the average pore radius (r) is, the slower the imbibition will be. The early stage of imbibition is mainly controlled by capillary force. In the middle and late stages of imbibition, as the imbibition distance increases, the imbibition resistance enhances, and the rate of imbibition gradually slows down. According to Equation (3) and Equation (16), the smaller the average pore radius (r) is, the bigger the capillary force and the imbibition resistance will be. The stronger the capillary force relative to the imbibition resistance, the smaller the deceleration of imbibition rate will be. Figure 7(b) shows that when the average pore radius is more than 50 nm, the maximum imbibition distance increases as the average pore radius increases; when the average pore radius is less than 50 nm, the maximum imbibition distance increases as the average pore radius decreases. Therefore, the imbibition is common controlled by capillary force and imbibition resistance during the middle and late stages of imbibition.Figure 7
The variation curves of imbibition rate with imbibition time under different average pore radius (r) and the curve maximum distance of imbibition with average pore radius.
(a)
(b)


## 4.1. Comparison and Analysis

Through research of the existing imbibition models, it is found that Cai’s model, which takes into account the size and shape of pores, the tortuosity of imbibition streamlines in random porous media, and the initial wetting-phase saturation, is generalized and can be simplified to other models, such as Fries–Dreyer model [8]. However, because the resistance of imbibition is little in conventional porous media with relatively large pores, Cai’s model omits the resistance term. As we all know, if there is no horizontal imbibition resistance in tight reservoirs, the distance of imbibition will be infinite, which is obviously unreasonable.Therefore, based on Cai’s model [8], our model considers the resistance of imbibition and the apparent viscosity of fluid in tight cores in order to apply in tight reservoirs. Fries–Dreyer model [18], Cai’s model, and our model can be formed into “dL/dt=a/L−b” mode. Table 1 shows the values of “a and b” in the mode.Table 1
Values of “a and b” in “dL/dt=a/L−b” mode.
Imbibition model
Direction
a
b
Washburn model
Vertical
r
σ
cos
θ
/
2
μ
0
Horizontal
r
σ
cos
θ
/
2
μ
0
Fries–Dreyer model
Vertical
r
σ
cos
θ
/
4
μ
r
2
ρ
g
/
8
μ
Horizontal
r
σ
cos
θ
/
4
μ
0
Cai’s model
Vertical
r
a
σ
cos
θ
/
4
μ
τ
2
r
a
2
ρ
g
/
8
μ
τ
2
Horizontal
r
a
σ
cos
θ
/
4
μ
τ
2
0
Our model
Vertical
r
a
σ
cos
θ
/
4
μ
c
τ
2
r
a
2
ρ
g
+
τ
G
/
8
μ
c
τ
2
Horizontal
r
a
σ
cos
θ
/
4
μ
c
τ
2
r
a
2
G
/
8
μ
c
τGenerally speaking, our model takes more geometrical parameters into account than other imbibition models in order to better match the real process of imbibition in tight cores. Accordingly, our model can derive other existing imbibition models under conditions ofG = 0, μc = μ as well.We choose tight sandstone data (ra = 1 × 10−7 m, σ  = 0.074 N·m, θ  = 80°, μc = 1.08 cp, τ  = 2, G = 0.43 MPa/s, ρ  = 1000 g/m3, and g = 9.8 m/s2) for calculation in these imbibition models at horizontal direction to quantitatively compare our model with the others.As shown in Figure2, since the Washburn Model and Fries–Dreyer Model do not take into account the tortuosity, the calculated value has a large gap with Cai’s and our model’s calculated value. The bigger the value of tortuosity is, the greater the difference is. Cai’s model fits well in the initial stage of our model. However, since Cai’s model does not consider the influence of resistance, there is no limit as the imbibition distance increases with imbibition time. As the imbibition continues, our model calculated curve tends to be flat, marking the end of the process of imbibition. Therefore, our model can be used for the calculation of the whole process of imbibition, especially in tight sandstone conditions.Figure 2
Quantitative comparison of our model to other models with tight sandstone data.


## 4.2. Model Validation

WhenG = 0 and μc = μ, our model is same as Cai’s model, which has been verified through porous materials with porosity of 0.045–0.7 and permeability of 5–1200 mD.The experimental data [32] of horizontal imbibition in tight cores (k < 1 mD) are collected to verify these models, and petrophysical parameters of tight cores for imbibition experiment are shown in Table 2.Table 2
Petrophysical parameters of tight cores for imbibition experiment.
Number
Diameter (cm)
Length (cm)
Porosity (%)
Permeability (mD)
Average radius (m)
Density of fluid (kg/m3)
Viscosity of fluid (MPa·s)
Exp1
2.53
3.118
9.48
0.14
0.118
1000
1.183
Exp2
2.52
2.794
5.85
0.215
0.294
1000
1.183
Exp3
2.51
3.000
6.89
0.0229
0.026
1000
1.183It can be seen from Figure3 that the result of Fries–Dreyer model has significant difference with the experimental data, and there is a good agreement between the result of Cai’s model and the experimental data in the early time when τ = 8, but Cai’s model cannot fit the experimental data at the later stage very well.Figure 3
Variation of imbibition distance with imbibition time in tight cores: simulation results of Fries–Dreyer Model and Cai’s Model with different tortuosity vs. experimental data of Exp1.In our model, the comprehensive coefficient (ζ) is defined to be tortuosity (τ) multiplied with apparent viscosity coefficient (ξ), which can characterize more influence parameters. Figure 4 shows that the result of our model (ζ = 7.5 and G = 0.4 MPa/m) can fit the experimental data better than other models. We also verify our model by different experimental data, as shown in Figure 5. It can be known that the resistance (G) is positively correlated with average radius, and the comprehensive coefficient (ζ) is negatively correlated with average radius.Figure 4
Variation of imbibition distance with imbibition time in tight cores: simulation results of our model with different composite factors (ζ) and fluid flow resistance (G) vs. experimental data of Exp1.Figure 5
Comparison between imbibition experimental data of different tight cores and calculated results of our model under different composite factors (ζ) and fluid flow resistance (G).As shown in Figure6, it can be concluded that the resistance (G) reduces the rate of imbibition, and the influence of resistance (G) increases with the increase of imbibition time. Therefore, our model can fit the result of experiment better. According to Equation (14) and Exp3, we can predict that the total time of imbibition is 56 days and the maximum distance of imbibition is 0.09 m in corresponding tight reservoirs.Figure 6
The variation curves of imbibition rate with imbibition time under different fluid flow resistances (G) in tight porous media.Due to a large amount of experimental data [33, 34], the empirical equation of imbibition resistance (G) and average pores radius (ra) can be written by(16)G=MeNra,where M and N are the empirical constants. According to experimental data of Exp1, 2, and 3 and their fitting data with our model, we can obtain the M = 3.2347 and N = −2 × 107.When the imbibition resistance (G) is not considered, the smaller the average pore radius (r) is, the slower the imbibition will be, and the imbibition distance is unlimited. When the imbibition resistance (G) is considered, as shown in Figure 7(a), in the early stage of imbibition, the smaller the average pore radius (r) is, the slower the imbibition will be. The early stage of imbibition is mainly controlled by capillary force. In the middle and late stages of imbibition, as the imbibition distance increases, the imbibition resistance enhances, and the rate of imbibition gradually slows down. According to Equation (3) and Equation (16), the smaller the average pore radius (r) is, the bigger the capillary force and the imbibition resistance will be. The stronger the capillary force relative to the imbibition resistance, the smaller the deceleration of imbibition rate will be. Figure 7(b) shows that when the average pore radius is more than 50 nm, the maximum imbibition distance increases as the average pore radius increases; when the average pore radius is less than 50 nm, the maximum imbibition distance increases as the average pore radius decreases. Therefore, the imbibition is common controlled by capillary force and imbibition resistance during the middle and late stages of imbibition.Figure 7
The variation curves of imbibition rate with imbibition time under different average pore radius (r) and the curve maximum distance of imbibition with average pore radius.
(a)
(b)


## 5. Conclusions

In this paper, a more generalized capillary imbibition model is derived by considering the resistance and the apparent viscosity of fluids in porous media. This model takes more geometrical parameters into account than other imbibition models to better match the result of imbibition experiment in tight cores and can derive other existing imbibition models in special cases. Through this mode and experimental data, we can predict the maximum distance and total time of imbibition of water in tight reservoirs, which will facilitate the study on water injection development in tight oil and gas reservoirs.



---

*Source: 4148734-2018-10-08.xml*