# A Pressure Transient Model for Power-Law Fluids in Porous Media Embedded with a Tree-Shaped Fractal Network

Authors: Xiao-Hua Tan; Jian-Yi Liu; Jia-Hui Zhao; Xiao-Ping Li; Guang-Dong Zhang; Chuan Tang
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

This work studies the pressure transient of power-law fluids in porous media embedded with a tree-shaped fractal network. A pressure transient model was created based on the fractal properties of tree-shaped capillaries, generalized Darcy’s law and constitutive equation for power-law fluids. The dimensionless pressure model was developed using the Laplace transform and Stehfest numerical inversion method. According to the model’s solution, the bi-logarithmic type curves of power-law fluids in porous media embedded with a tree-shaped fractal network are illustrated. The influences of different fractal factors and Power-law fluids parameters on pressure transient responses are discussed.


---

## Body

## 1. Introduction

Power-law fluids flow in porous media has always been a subject of great interest owing to its fundamental and pragmatic significance. Consequently, over the years, a voluminous amount of research has been conducted to gain insight into phenomena that is related to fluids flow. Traditionally, numerical methods have played a significant role in the analysis of power-law fluids flow in porous media. Lopez et al. [1] determined the flow of power-law fluids in porous media using network models. Through the use of numerical methods, a network model is capable of predicting each of the parameters associated with a pore, such as pressure drop, velocity, and superficial viscosity. Network model predictions are superior to those that are based on simplified capillary bundle models. Kumari and Nath [2] used numerical methods to study non-Darcy mixed convection of power-law fluids in porous media. Orgéas et al. [3] created a flow model for power-law fluids in anisotropic porous media at the low-pore Reynolds number. Pascal [4] researched the transient flow of power-law fluids in a porous incline using the Kármán momentum integral method. Non-Newtonian constitutive equations, used as one of the numerical methods, have been incorporated into the conventional lattice Boltzmann method (LBM) to simulate flow of power-law fluids [5–9]. However, the LBM demands tedious calculations and necessitates accuracy of model. These constraints limit the LBM’s ability to describe flow of power-law fluids in porous media.Yun et al. [10] have proposed the starting pressure gradient model for the flow of Bingham fluids in fractal porous media, based on the fractal nature of pore size distribution in porous media. They also analyzed the effect of capillary pressure on the starting pressure gradient. Results showed that the effect of capillary pressure on the starting pressure gradient at a low porosity is significant and cannot be neglected. However, at high porosity, the effect of capillary pressure on the starting pressure gradient is found to be negligible. Li and Yu [11] proposed a starting pressure gradient model for Bingham fluids flow in fractal networks, which could be used to analyze the influence of network microstructural parameters on the starting pressure gradient and on network permeability. The model can also be used to study Bingham fluids flow in dual porosity media. The theory proposed by Yun et al. [10] and Li and Yu [11] was further developed by Wang and Yu [12]. They derived the starting pressure gradient model for Bingham fluids flow in tree-shaped fractal porous media by embedding the tree-shaped fractal network into porous media. Based on the tree-shaped fractal model, Wang et al. [13] proposed a starting pressure gradient model for Bingham fluids flow in porous media that is embedded with randomly distributed tree-shaped fractal networks.


## 2. Physical Modeling

Figure1 shows a coordinate system for radial flow toward the center point, O, from the external boundary. The center point, O, has an internal boundary, r
in. h is the system thickness. The coordinate system consists of two parts: the tree-shaped fractal network and the porous media. During a tree-shaped fractal network generation, the branches on each level must end up on the same circle with all circles having the same center of origin, O. N tubes, which start at O, make up the tree-shaped fractal network. The tube’s initial length and diameter are l
0 and d
0, respectively. The double branches (
b
=
2
), whose angles are θ
(
θ
<
π
/
2
) and total network branch levels are M, are applied in this network. Furthermore, two scale factors are used in this fractal network, length ratio, α, and diameter ratio, β.Figure 1
Porous media embedded with a tree-shaped fractal network.The physic model assumptions are as follows.(1)
The branch tube in tree-shaped fractal network is assumed to be smooth and tube wall thickness is ignored.
(2)
The porous media is divided intoM annular sections by a tree-shaped fractal network.
(3)
Each of the tree-shaped fractal network section’s properties is different, but each of the porous media section’s properties is identical.
(4)
Porous media permeability is much lower than tree-shaped fractal network permeability, and fluids only flow to the center point,O, through the tree-shaped fractal network.
(5)
Rock and liquid are considered slightly compressible, with each having a constant and small compressibility.
(6)
Isothermal and single direction flow is considered.
(7)
Capillary pressure and gravity effects are neglected.
(8)
Non-Newtonian fluids obey the power-law principle and are considered as pseudoplastic fluids(
0
<
n
<
1
).
(9)
The internal boundary is considered to have constant flow rate, while the external boundary is considered to be closed or at constant pressure.
(10)
At timet
=
0, pressure is uniformly distributed in all system and is equal to initial pressure (p
i).


## 3. Mathematical Modeling

### 3.1. Flow Parameters in the Tree-Shaped Fractal Network

The flow velocity through a single tube for power-law fluids is given by [14]
(1)
v
=
d
1
+
(
1
/
n
)
2
1
+
(
2
/
n
)
m
1
/
n
(
3
+
(
1
/
n
)
)
(
-
d
p
d
L
)
1
/
n
,
where τ is shear stress, m is the consistency index, γ
˙ is shear rate, and n is the power index.Ikoku and Ramey Jr [15] defined the flow velocity for power-law fluids as
(2)
v
=
(
-
K
μ
eff
d
p
d
L
)
1
/
n
,
where K is permeability and μ
eff is effective viscosity.Substituting (1) into (2), permeability for power-law fluids in a single tube can be expressed by effective viscosity:
(3)
K
=
μ
eff
d
n
+
1
2
n
+
2
m
(
3
+
(
1
/
n
)
)
n
.Whenn
=
1 and m
=
μ
eff
=
μ in (3), the permeability expression for power-law fluids reduces to that for Newtonian fluids
(4)
K
=
d
2
32
which agrees with the research result of Xu et al. [16].To characterize the branching structures, let the length and diameter of a typical branch at some intermediate levelk
(
k
=
0,1
,
…
,
M
) be d
k and l
k, respectively. Two scale factors length ratio, α, and diameter ratio, β, are defined as
(5)
α
=
l
k
+
1
l
k
,
β
=
d
k
+
1
d
k
.For thekth level tube, length and diameter are given by
(6)
l
k
=
α
k
l
0
,
(7)
d
k
=
β
k
d
0
.The interface radius is defined as the distance from the internal boundary to each section’s boundary [16]:
(8)
r
k
=
l
0
[
1
+
α
(
1
-
α
k
)
cos
⁡
θ
1
-
α
]
.The external flow radius is defined as the distance from the internal boundary to the external boundary:(9)
r
ex
=
r
M
=
l
0
[
1
+
α
(
1
-
α
M
)
cos
⁡
θ
1
-
α
]
.The single tube permeability inkth section is expressed by
(10)
K
f
s
k
=
μ
eff
d
k
n
+
1
2
n
+
2
m
(
3
+
(
1
/
n
)
)
n
.Consideringkth section in tree-shaped fractal network tortuosity T
k, the single tube permeability in kth section can be expressed as
(11)
K
f
s
k
=
μ
eff
d
k
n
+
1
2
n
+
2
m
T
k
(
3
+
(
1
/
n
)
)
n
,
where
(12)
T
k
=
l
k
r
k
-
r
k
-
1
=
{
1
k
=
0
1
cos
⁡
θ
k
>
0
.Substituting (7) into (11), the permeability in kth section of tree-shaped fractal network is
(13)
K
f
k
=
{
N
μ
eff
d
k
n
+
1
2
n
+
2
m
(
3
+
(
1
/
n
)
)
n
k
=
0
N
μ
eff
b
k
β
2
k
d
k
n
+
1
cos
⁡
θ
2
n
+
2
m
T
k
(
3
+
(
1
/
n
)
)
n
k
>
0
.The total system volume in thekth section of tree-shaped fractal network can be calculated by
(14)
V
k
=
{
π
h
r
0
2
k
=
0
π
h
(
r
k
2
-
r
k
-
1
2
)
k
>
0
.The volume in thekth section of tree-shaped fractal network can be calculated by
(15)
V
f
k
=
N
π
b
k
α
k
β
2
k
l
0
d
0
2
4
.The tree-shaped fractal network porosity in thekth section of tree-shaped fractal network can be obtained by dividing (15) by (14):
(16)
ϕ
f
k
=
V
f
k
V
k
=
{
N
d
0
2
4
h
l
0
k
=
0
N
b
k
α
k
β
2
k
l
0
d
0
2
4
h
(
r
k
2
-
r
k
-
1
2
)
k
>
0
.Under the condition of double branches(
b
=
2
), the permeability, K
f
k, and porosity, ϕ
f
k, in the kth section of tree-shaped fractal network are equal; that is,
(17)
K
f
k
=
K
f
(
k
+
1
)
ϕ
f
k
=
ϕ
f
(
k
+
1
)
.Substituting (13) and (16) into (17), we can obtain b
=
2, α
=
1, and β
=
0.707.


### 3.2. Mathematical Model

According to the physical model, the mathematical pressure transient model for power-law fluids in porous media embedded with a tree-shaped fractal network can be described as follows.Governing differential equations for tree-shaped fractal network [17]:
(18)
K
k
n
μ
eff
(
q
2
π
h
)
1
-
n
r
n
-
1
(
∂
2
p
f
k
∂
r
2
+
n
r
∂
p
f
k
∂
r
)
+
a
K
m
μ
eff
(
p
m
k
-
p
f
k
)
=
(
ϕ
C
t
)
f
k
∂
p
f
k
∂
t
r
k
-
1
≤
r
≤
r
k
.Governing differential equations for porous media [17]:
(19)
-
a
K
m
m
u
eff
(
p
m
k
-
p
f
k
)
=
(
ϕ
C
t
)
m
∂
p
m
k
∂
t
r
k
-
1
≤
r
≤
r
k
.Initial condition:(20)
p
f
k
(
r
,
0
)
=
p
m
k
(
r
,
0
)
=
p
i
(
k
=
0,1
,
…
,
M
)
.Interface connecting conditions of each section, pressure continuity [18]:
(21)
p
f
k
|
r
=
r
k
=
p
f
(
k
+
1
)
|
r
=
r
k
(
k
=
0,1
,
…
,
M
-
1
)
.Interface connecting conditions of each section, rate continuity [18]:
(22)
∂
p
f
k
∂
r
|
r
=
r
k
=
k
f
(
k
+
1
)
k
f
k
∂
p
f
(
k
+
1
)
∂
r
|
r
=
r
k
(
k
=
0,1
,
…
,
M
-
1
)
.Internal boundary condition [19]:
(23)
p
in
=
p
f
0
|
r
=
r
in
.External boundary condition (infinite):(24)
p
f
M
(
r
⟶
∞
,
t
)
=
p
i
.External boundary condition (constant pressure):(25)
p
f
M
(
r
=
r
M
,
t
)
=
p
i
.External boundary condition (closed):(26)
∂
p
f
M
∂
r
|
r
=
r
M
=
0
.


### 3.3. Solution to Mathematical Model

The dimensionless mathematical model for power-law fluids in porous media embedded with a tree-shaped fractal network is shown in the AppendixB. Introduce the Laplace transform based on t
D. The dimensionless mathematical model in Laplace space is as follows.Governing differential equations for tree-shaped fractal network:(27)
d
2
p
-
D
f
k
d
r
D
2
+
n
r
D
d
p
-
D
f
k
d
r
D
+
λ
k
r
D
1
-
n
(
K
f
k
K
m
p
-
D
m
k
-
p
-
D
f
k
)
=
ω
k
r
D
1
-
n
K
f
0
K
f
k
z
p
-
D
f
k
r
D
(
k
-
1
)
≤
r
D
≤
r
D
k
.Governing differential equations for porous media:(28)
-
λ
k
(
p
-
D
m
k
-
K
m
K
f
k
p
-
D
f
k
)
=
(
1
-
ω
k
)
K
f
0
K
f
k
z
p
-
D
m
k
r
D
(
k
-
1
)
≤
r
D
≤
r
D
k
.Initial condition:(29)
p
-
D
f
k
(
r
D
,
0
)
=
p
-
D
m
k
(
r
D
,
0
)
=
0
(
k
=
0,1
,
…
,
M
)
.Interface connecting conditions of each zone, pressure continuity:(30)
p
-
D
f
k
|
r
D
=
r
D
k
=
p
-
D
f
(
k
+
1
)
|
r
D
=
r
D
k
(
k
=
0,1
,
…
,
M
-
1
)
.Interface connecting conditions of each zone, rate continuity:(31)
d
p
-
D
f
k
d
r
D
e
|
r
D
=
r
D
k
=
K
f
(
k
+
1
)
K
f
k
d
p
-
D
f
(
k
+
1
)
d
r
D
|
r
D
=
r
D
k
(
k
=
0,1
,
…
,
M
-
1
)
.Internal boundary condition:(32)
d
p
-
D
in
d
r
D
|
r
D
=
1
=
-
1
.External boundary condition (infinite):(33)
p
-
D
f
M
(
r
D
⟶
∞
,
z
)
=
0
.External boundary condition (constant pressure):(34)
p
-
D
f
M
(
r
D
=
r
D
M
,
z
)
=
0
.External boundary condition (closed):(35)
d
p
-
D
f
M
d
r
D
|
r
D
=
r
D
M
=
0
.Substituting (28) into (29), we can obtain
(36)
d
2
p
-
D
f
k
d
r
D
2
+
n
r
D
k
d
p
-
D
f
k
d
r
D
-
r
D
k
1
-
n
z
S
(
z
)
p
-
D
f
k
=
0
(
k
=
0,1
,
…
,
M
)
.In (36), The dimensionless interface radius, r
D
k, and the function, S
k
(
z
), can be directly expressed by the tree-shaped fractal network factors.The dimensionless interface radius expression can be derived by substituting (8) into (A.3):
(37)
r
D
k
=
l
0
r
in
[
1
+
α
(
1
-
α
k
)
cos
⁡
θ
1
-
α
]
.The expression of the functionS
k
(
z
) is as follows:
(38)
S
k
(
z
)
=
(
K
f
0
/
K
f
k
)
λ
k
(
1
-
ω
k
)
(
K
f
0
/
K
f
k
)
(
1
-
ω
k
)
z
+
λ
k
+
K
f
0
K
f
k
ω
k
,
where
(39)
K
f
0
K
f
k
=
(
b
β
2
)
-
k
.The fluid capacitance coefficient,ω
k, can be obtained by substituting (16) into (A.5), and the interporosity flow coefficient, λ
k, can be obtained by substituting (13) into (A.6).The general solution ofp
-
D
f
k in (36) is calculated
(40)
p
-
D
f
k
=
A
k
f
I
(
r
D
k
,
S
k
(
z
)
)
+
B
k
f
K
(
r
D
k
,
S
k
(
z
)
)
(
k
=
0,1
,
…
,
M
)
,
where
(41)
f
I
(
r
D
k
,
S
k
(
z
)
)
=
r
D
k
(
1
-
n
)
/
2
I
v
(
r
D
k
(
3
-
n
)
/
2
2
S
k
(
z
)
3
-
n
)
,
f
K
(
r
D
k
,
S
k
(
z
)
)
=
r
D
k
(
1
-
n
)
/
2
K
v
(
r
D
k
(
3
-
n
)
/
2
2
S
k
(
z
)
3
-
n
)
,
v
=
(
1
-
n
)
(
3
-
n
)
.Then, the derivative ofp
-
D
f
k in (40) can be expressed as
(42)
d
p
-
D
f
k
d
r
D
=
A
k
f
I
(
r
D
k
,
S
k
(
z
)
)
d
r
D
+
B
k
f
K
(
r
D
k
,
S
k
(
z
)
)
d
r
D
(
k
=
0,1
,
…
,
M
)
.


## 3.1. Flow Parameters in the Tree-Shaped Fractal Network

The flow velocity through a single tube for power-law fluids is given by [14]
(1)
v
=
d
1
+
(
1
/
n
)
2
1
+
(
2
/
n
)
m
1
/
n
(
3
+
(
1
/
n
)
)
(
-
d
p
d
L
)
1
/
n
,
where τ is shear stress, m is the consistency index, γ
˙ is shear rate, and n is the power index.Ikoku and Ramey Jr [15] defined the flow velocity for power-law fluids as
(2)
v
=
(
-
K
μ
eff
d
p
d
L
)
1
/
n
,
where K is permeability and μ
eff is effective viscosity.Substituting (1) into (2), permeability for power-law fluids in a single tube can be expressed by effective viscosity:
(3)
K
=
μ
eff
d
n
+
1
2
n
+
2
m
(
3
+
(
1
/
n
)
)
n
.Whenn
=
1 and m
=
μ
eff
=
μ in (3), the permeability expression for power-law fluids reduces to that for Newtonian fluids
(4)
K
=
d
2
32
which agrees with the research result of Xu et al. [16].To characterize the branching structures, let the length and diameter of a typical branch at some intermediate levelk
(
k
=
0,1
,
…
,
M
) be d
k and l
k, respectively. Two scale factors length ratio, α, and diameter ratio, β, are defined as
(5)
α
=
l
k
+
1
l
k
,
β
=
d
k
+
1
d
k
.For thekth level tube, length and diameter are given by
(6)
l
k
=
α
k
l
0
,
(7)
d
k
=
β
k
d
0
.The interface radius is defined as the distance from the internal boundary to each section’s boundary [16]:
(8)
r
k
=
l
0
[
1
+
α
(
1
-
α
k
)
cos
⁡
θ
1
-
α
]
.The external flow radius is defined as the distance from the internal boundary to the external boundary:(9)
r
ex
=
r
M
=
l
0
[
1
+
α
(
1
-
α
M
)
cos
⁡
θ
1
-
α
]
.The single tube permeability inkth section is expressed by
(10)
K
f
s
k
=
μ
eff
d
k
n
+
1
2
n
+
2
m
(
3
+
(
1
/
n
)
)
n
.Consideringkth section in tree-shaped fractal network tortuosity T
k, the single tube permeability in kth section can be expressed as
(11)
K
f
s
k
=
μ
eff
d
k
n
+
1
2
n
+
2
m
T
k
(
3
+
(
1
/
n
)
)
n
,
where
(12)
T
k
=
l
k
r
k
-
r
k
-
1
=
{
1
k
=
0
1
cos
⁡
θ
k
>
0
.Substituting (7) into (11), the permeability in kth section of tree-shaped fractal network is
(13)
K
f
k
=
{
N
μ
eff
d
k
n
+
1
2
n
+
2
m
(
3
+
(
1
/
n
)
)
n
k
=
0
N
μ
eff
b
k
β
2
k
d
k
n
+
1
cos
⁡
θ
2
n
+
2
m
T
k
(
3
+
(
1
/
n
)
)
n
k
>
0
.The total system volume in thekth section of tree-shaped fractal network can be calculated by
(14)
V
k
=
{
π
h
r
0
2
k
=
0
π
h
(
r
k
2
-
r
k
-
1
2
)
k
>
0
.The volume in thekth section of tree-shaped fractal network can be calculated by
(15)
V
f
k
=
N
π
b
k
α
k
β
2
k
l
0
d
0
2
4
.The tree-shaped fractal network porosity in thekth section of tree-shaped fractal network can be obtained by dividing (15) by (14):
(16)
ϕ
f
k
=
V
f
k
V
k
=
{
N
d
0
2
4
h
l
0
k
=
0
N
b
k
α
k
β
2
k
l
0
d
0
2
4
h
(
r
k
2
-
r
k
-
1
2
)
k
>
0
.Under the condition of double branches(
b
=
2
), the permeability, K
f
k, and porosity, ϕ
f
k, in the kth section of tree-shaped fractal network are equal; that is,
(17)
K
f
k
=
K
f
(
k
+
1
)
ϕ
f
k
=
ϕ
f
(
k
+
1
)
.Substituting (13) and (16) into (17), we can obtain b
=
2, α
=
1, and β
=
0.707.


## 3.2. Mathematical Model

According to the physical model, the mathematical pressure transient model for power-law fluids in porous media embedded with a tree-shaped fractal network can be described as follows.Governing differential equations for tree-shaped fractal network [17]:
(18)
K
k
n
μ
eff
(
q
2
π
h
)
1
-
n
r
n
-
1
(
∂
2
p
f
k
∂
r
2
+
n
r
∂
p
f
k
∂
r
)
+
a
K
m
μ
eff
(
p
m
k
-
p
f
k
)
=
(
ϕ
C
t
)
f
k
∂
p
f
k
∂
t
r
k
-
1
≤
r
≤
r
k
.Governing differential equations for porous media [17]:
(19)
-
a
K
m
m
u
eff
(
p
m
k
-
p
f
k
)
=
(
ϕ
C
t
)
m
∂
p
m
k
∂
t
r
k
-
1
≤
r
≤
r
k
.Initial condition:(20)
p
f
k
(
r
,
0
)
=
p
m
k
(
r
,
0
)
=
p
i
(
k
=
0,1
,
…
,
M
)
.Interface connecting conditions of each section, pressure continuity [18]:
(21)
p
f
k
|
r
=
r
k
=
p
f
(
k
+
1
)
|
r
=
r
k
(
k
=
0,1
,
…
,
M
-
1
)
.Interface connecting conditions of each section, rate continuity [18]:
(22)
∂
p
f
k
∂
r
|
r
=
r
k
=
k
f
(
k
+
1
)
k
f
k
∂
p
f
(
k
+
1
)
∂
r
|
r
=
r
k
(
k
=
0,1
,
…
,
M
-
1
)
.Internal boundary condition [19]:
(23)
p
in
=
p
f
0
|
r
=
r
in
.External boundary condition (infinite):(24)
p
f
M
(
r
⟶
∞
,
t
)
=
p
i
.External boundary condition (constant pressure):(25)
p
f
M
(
r
=
r
M
,
t
)
=
p
i
.External boundary condition (closed):(26)
∂
p
f
M
∂
r
|
r
=
r
M
=
0
.


## 3.3. Solution to Mathematical Model

The dimensionless mathematical model for power-law fluids in porous media embedded with a tree-shaped fractal network is shown in the AppendixB. Introduce the Laplace transform based on t
D. The dimensionless mathematical model in Laplace space is as follows.Governing differential equations for tree-shaped fractal network:(27)
d
2
p
-
D
f
k
d
r
D
2
+
n
r
D
d
p
-
D
f
k
d
r
D
+
λ
k
r
D
1
-
n
(
K
f
k
K
m
p
-
D
m
k
-
p
-
D
f
k
)
=
ω
k
r
D
1
-
n
K
f
0
K
f
k
z
p
-
D
f
k
r
D
(
k
-
1
)
≤
r
D
≤
r
D
k
.Governing differential equations for porous media:(28)
-
λ
k
(
p
-
D
m
k
-
K
m
K
f
k
p
-
D
f
k
)
=
(
1
-
ω
k
)
K
f
0
K
f
k
z
p
-
D
m
k
r
D
(
k
-
1
)
≤
r
D
≤
r
D
k
.Initial condition:(29)
p
-
D
f
k
(
r
D
,
0
)
=
p
-
D
m
k
(
r
D
,
0
)
=
0
(
k
=
0,1
,
…
,
M
)
.Interface connecting conditions of each zone, pressure continuity:(30)
p
-
D
f
k
|
r
D
=
r
D
k
=
p
-
D
f
(
k
+
1
)
|
r
D
=
r
D
k
(
k
=
0,1
,
…
,
M
-
1
)
.Interface connecting conditions of each zone, rate continuity:(31)
d
p
-
D
f
k
d
r
D
e
|
r
D
=
r
D
k
=
K
f
(
k
+
1
)
K
f
k
d
p
-
D
f
(
k
+
1
)
d
r
D
|
r
D
=
r
D
k
(
k
=
0,1
,
…
,
M
-
1
)
.Internal boundary condition:(32)
d
p
-
D
in
d
r
D
|
r
D
=
1
=
-
1
.External boundary condition (infinite):(33)
p
-
D
f
M
(
r
D
⟶
∞
,
z
)
=
0
.External boundary condition (constant pressure):(34)
p
-
D
f
M
(
r
D
=
r
D
M
,
z
)
=
0
.External boundary condition (closed):(35)
d
p
-
D
f
M
d
r
D
|
r
D
=
r
D
M
=
0
.Substituting (28) into (29), we can obtain
(36)
d
2
p
-
D
f
k
d
r
D
2
+
n
r
D
k
d
p
-
D
f
k
d
r
D
-
r
D
k
1
-
n
z
S
(
z
)
p
-
D
f
k
=
0
(
k
=
0,1
,
…
,
M
)
.In (36), The dimensionless interface radius, r
D
k, and the function, S
k
(
z
), can be directly expressed by the tree-shaped fractal network factors.The dimensionless interface radius expression can be derived by substituting (8) into (A.3):
(37)
r
D
k
=
l
0
r
in
[
1
+
α
(
1
-
α
k
)
cos
⁡
θ
1
-
α
]
.The expression of the functionS
k
(
z
) is as follows:
(38)
S
k
(
z
)
=
(
K
f
0
/
K
f
k
)
λ
k
(
1
-
ω
k
)
(
K
f
0
/
K
f
k
)
(
1
-
ω
k
)
z
+
λ
k
+
K
f
0
K
f
k
ω
k
,
where
(39)
K
f
0
K
f
k
=
(
b
β
2
)
-
k
.The fluid capacitance coefficient,ω
k, can be obtained by substituting (16) into (A.5), and the interporosity flow coefficient, λ
k, can be obtained by substituting (13) into (A.6).The general solution ofp
-
D
f
k in (36) is calculated
(40)
p
-
D
f
k
=
A
k
f
I
(
r
D
k
,
S
k
(
z
)
)
+
B
k
f
K
(
r
D
k
,
S
k
(
z
)
)
(
k
=
0,1
,
…
,
M
)
,
where
(41)
f
I
(
r
D
k
,
S
k
(
z
)
)
=
r
D
k
(
1
-
n
)
/
2
I
v
(
r
D
k
(
3
-
n
)
/
2
2
S
k
(
z
)
3
-
n
)
,
f
K
(
r
D
k
,
S
k
(
z
)
)
=
r
D
k
(
1
-
n
)
/
2
K
v
(
r
D
k
(
3
-
n
)
/
2
2
S
k
(
z
)
3
-
n
)
,
v
=
(
1
-
n
)
(
3
-
n
)
.Then, the derivative ofp
-
D
f
k in (40) can be expressed as
(42)
d
p
-
D
f
k
d
r
D
=
A
k
f
I
(
r
D
k
,
S
k
(
z
)
)
d
r
D
+
B
k
f
K
(
r
D
k
,
S
k
(
z
)
)
d
r
D
(
k
=
0,1
,
…
,
M
)
.


## 4. Flow Behavior Characteristics

### 4.1. Flow Regime Recognition

p
-
D
f
k, A
k and B
k
(
k
=
0,1
,
…
,
M
) can be obtained by solving the general solution equations (40) and (42) and definite condition equations (32)–(35). We can obtain dimensionless internal boundary pressure in real space, p
D
in, by the Stehfest numerical inversion method. The bilogarithmic type curves for power-law fluids in porous media embedded with a tree-shaped fractal network can be illustrated. In a condition of closed external boundary, the transient flow process, which has four flow regimes, can be clearly shown (Figure 2). The full and dashed lines represent pressure and pressure derivative curves, respectively.Regime 1 is interporosity flow in tree-shaped fractal network regime. Power-law fluids start to flow in tree-shaped fractal network. Dimensionless pressure increases rapidly before slowing slightly, because it is hard to flow for power-law fluids (0
<
n
<
1) at the beginning.
Regime 2 is radial flow in tree-shaped fractal network regime. All fluids have flowed in tree-shaped fractal network, and pressure waves spread all over the tree-shaped fractal network. The pressure derivative curve converges to a horizontal line, which depicts the response of the pressure dynamic in tree-shaped fractal network.
Regime 3 is the interporosity flow regime of porous media to tree-shaped fractal network. Power-law fluids in porous media start to flow into tree-shaped fractal network. The pressure derivative curve is V-shaped, which depicts the response of interporosity flow between the tree-shaped fractal network and porous media. It is influenced by spread of the pressure wave through the porous media.
Regime 4 is the total system radial flow regime in porous media embedded with a tree-shaped fractal network. The pressure derivative curve converges to a horizontal line, which depicts the response of the pressure dynamic balance state in the whole system.Figure 2
Pressure type curves for power-law fluids in porous media embedded with a tree-shaped fractal network (n
=
0.9, α
=
1, β
=
0.707, θ
=
1, N
=
4, M
=
10, l
0
=
10
m, d
0
=
0.05
m, r
in
=
0.1
m, h
=
5
m, C
m
t
=
2.2
×
10
-
5
MP
a
-
1, and C
f
t
=
1
×
10
-
4
MP
a
-
1).


### 4.2. Parameter Influence

Figure3 shows the type curve characteristics affected by the power index of power-law fluids, n. As n decreases, it is harder to flow for power-law fluids (
0
<
n
<
1
) in tree-shaped fractal network and porous media. A larger n leads to a lower location of the dimensionless pressure curve and regime 4 occurs earlier, resulting in a deeper and wider V-shaped type curve. When n
=
1, the power-law fluids pressure transient model reduces to the Newtonian fluids pressure transient model. The horizontal line representing regime 4 equals 0.5.Figure 3
Effect of power index (n) on type curves (α
=
1, β
=
0.707, θ
=
1, N
=
4, M
=
10, l
0
=
10
m, d
0
=
0.05
m, r
in
=
0.1
m, h
=
5
m, C
m
t
=
2.2
×
10
-
5
MP
a
-
1, and C
f
t
=
1
×
10
-
4
MP
a
-
1).Figure4 exhibits the type curve characteristics affected by the length ratio, α. The porosity of tree-shaped fractal network, ϕ
f, decreases with an increase in α. A smaller ϕ
f leads to a lower supplying capacity in tree-shaped fractal network, which would result in a lower location of the dimensionless pressure curve in regimes 2, 3, and 4 and an earlier transition to a longer lasting regime 3. When α increases, regime 3 occurs earlier, which manifests in a deeper and wider V-shaped type curve.Figure 4
Effect of length ratio (α) on type curves (n
=
0.9, β
=
0.707, θ
=
1, N
=
4, M
=
10, l
0
=
10
m, d
0
=
0.05
m, r
in
=
0.1
m, h
=
5
m, C
m
t
=
2.2
×
10
-
5
MPa
-
1, and C
f
t
=
1
×
10
-
4
MPa
-
1).Figure5 exhibits the type curve characteristics affected by diameter ratio, β. The permeability of tree-shaped fractal network, K
f, and the porosity of tree-shaped fractal network, ϕ
f, increase with an increase in β. A larger K
f leads to larger flow capacity in tree-shaped fractal network and a later transition to regime 2. A larger ϕ
f leads to greater supplying capacity in tree-shaped fractal network with a later and shorter regime 3. A larger β leads to a lower location of the dimensionless pressure curve and regime 3 occurs later, resulting in a shallower and narrower V-shaped type curve. When β equals the critical value of 0.707 (see in (17)), the horizontal line representing regime 4 equals 0.5. When β is smaller than 0.707, K
f increases with an increase in radius, r, and the horizontal line is higher than 0.5. When β is greater than 0.707, K
f decreases with radius, r, and the horizontal line is higher than 0.5.Figure 5
Effect of diameter ratio (β) on type curves (n
=
0.9, α
=
1, θ
=
1, N
=
4, M
=
10, l
0
=
10
m, d
0
=
0.05
m, r
in
=
0.1
m, h
=
5
m, C
m
t
=
2.2
×
10
-
5
MP
a
-
1, and C
f
t
=
1
×
10
-
4
MP
a
-
1).Figure6 shows the type curve characteristics affected by branch angle, θ. Permeability of tree-shaped fractal network, K
f, decreases with an increase in θ. A smaller K
f leads to lower flow capacity in tree-shaped fractal network, with an earlier transition to regime 3. Porosity of tree-shaped fractal network, ϕ
f, increases with an increase in θ. A larger ϕ
f leads to greater supplying capacity in the tree-shaped fractal network, and regime 3 occurs later with a shorter duration.Figure 6
Effect of branch angle (θ) on type curves (n
=
0.9, α
=
1, β
=
0.707, N
=
4, M
=
10, l
0
=
10
m, d
0
=
0.05
m, r
in
=
0.1
m, h
=
5
m, C
m
t
=
2.2
×
10
-
5
MP
a
-
1, and C
f
t
=
1
×
10
-
4
MP
a
-
1).


## 4.1. Flow Regime Recognition

p
-
D
f
k, A
k and B
k
(
k
=
0,1
,
…
,
M
) can be obtained by solving the general solution equations (40) and (42) and definite condition equations (32)–(35). We can obtain dimensionless internal boundary pressure in real space, p
D
in, by the Stehfest numerical inversion method. The bilogarithmic type curves for power-law fluids in porous media embedded with a tree-shaped fractal network can be illustrated. In a condition of closed external boundary, the transient flow process, which has four flow regimes, can be clearly shown (Figure 2). The full and dashed lines represent pressure and pressure derivative curves, respectively.Regime 1 is interporosity flow in tree-shaped fractal network regime. Power-law fluids start to flow in tree-shaped fractal network. Dimensionless pressure increases rapidly before slowing slightly, because it is hard to flow for power-law fluids (0
<
n
<
1) at the beginning.
Regime 2 is radial flow in tree-shaped fractal network regime. All fluids have flowed in tree-shaped fractal network, and pressure waves spread all over the tree-shaped fractal network. The pressure derivative curve converges to a horizontal line, which depicts the response of the pressure dynamic in tree-shaped fractal network.
Regime 3 is the interporosity flow regime of porous media to tree-shaped fractal network. Power-law fluids in porous media start to flow into tree-shaped fractal network. The pressure derivative curve is V-shaped, which depicts the response of interporosity flow between the tree-shaped fractal network and porous media. It is influenced by spread of the pressure wave through the porous media.
Regime 4 is the total system radial flow regime in porous media embedded with a tree-shaped fractal network. The pressure derivative curve converges to a horizontal line, which depicts the response of the pressure dynamic balance state in the whole system.Figure 2
Pressure type curves for power-law fluids in porous media embedded with a tree-shaped fractal network (n
=
0.9, α
=
1, β
=
0.707, θ
=
1, N
=
4, M
=
10, l
0
=
10
m, d
0
=
0.05
m, r
in
=
0.1
m, h
=
5
m, C
m
t
=
2.2
×
10
-
5
MP
a
-
1, and C
f
t
=
1
×
10
-
4
MP
a
-
1).


## 4.2. Parameter Influence

Figure3 shows the type curve characteristics affected by the power index of power-law fluids, n. As n decreases, it is harder to flow for power-law fluids (
0
<
n
<
1
) in tree-shaped fractal network and porous media. A larger n leads to a lower location of the dimensionless pressure curve and regime 4 occurs earlier, resulting in a deeper and wider V-shaped type curve. When n
=
1, the power-law fluids pressure transient model reduces to the Newtonian fluids pressure transient model. The horizontal line representing regime 4 equals 0.5.Figure 3
Effect of power index (n) on type curves (α
=
1, β
=
0.707, θ
=
1, N
=
4, M
=
10, l
0
=
10
m, d
0
=
0.05
m, r
in
=
0.1
m, h
=
5
m, C
m
t
=
2.2
×
10
-
5
MP
a
-
1, and C
f
t
=
1
×
10
-
4
MP
a
-
1).Figure4 exhibits the type curve characteristics affected by the length ratio, α. The porosity of tree-shaped fractal network, ϕ
f, decreases with an increase in α. A smaller ϕ
f leads to a lower supplying capacity in tree-shaped fractal network, which would result in a lower location of the dimensionless pressure curve in regimes 2, 3, and 4 and an earlier transition to a longer lasting regime 3. When α increases, regime 3 occurs earlier, which manifests in a deeper and wider V-shaped type curve.Figure 4
Effect of length ratio (α) on type curves (n
=
0.9, β
=
0.707, θ
=
1, N
=
4, M
=
10, l
0
=
10
m, d
0
=
0.05
m, r
in
=
0.1
m, h
=
5
m, C
m
t
=
2.2
×
10
-
5
MPa
-
1, and C
f
t
=
1
×
10
-
4
MPa
-
1).Figure5 exhibits the type curve characteristics affected by diameter ratio, β. The permeability of tree-shaped fractal network, K
f, and the porosity of tree-shaped fractal network, ϕ
f, increase with an increase in β. A larger K
f leads to larger flow capacity in tree-shaped fractal network and a later transition to regime 2. A larger ϕ
f leads to greater supplying capacity in tree-shaped fractal network with a later and shorter regime 3. A larger β leads to a lower location of the dimensionless pressure curve and regime 3 occurs later, resulting in a shallower and narrower V-shaped type curve. When β equals the critical value of 0.707 (see in (17)), the horizontal line representing regime 4 equals 0.5. When β is smaller than 0.707, K
f increases with an increase in radius, r, and the horizontal line is higher than 0.5. When β is greater than 0.707, K
f decreases with radius, r, and the horizontal line is higher than 0.5.Figure 5
Effect of diameter ratio (β) on type curves (n
=
0.9, α
=
1, θ
=
1, N
=
4, M
=
10, l
0
=
10
m, d
0
=
0.05
m, r
in
=
0.1
m, h
=
5
m, C
m
t
=
2.2
×
10
-
5
MP
a
-
1, and C
f
t
=
1
×
10
-
4
MP
a
-
1).Figure6 shows the type curve characteristics affected by branch angle, θ. Permeability of tree-shaped fractal network, K
f, decreases with an increase in θ. A smaller K
f leads to lower flow capacity in tree-shaped fractal network, with an earlier transition to regime 3. Porosity of tree-shaped fractal network, ϕ
f, increases with an increase in θ. A larger ϕ
f leads to greater supplying capacity in the tree-shaped fractal network, and regime 3 occurs later with a shorter duration.Figure 6
Effect of branch angle (θ) on type curves (n
=
0.9, α
=
1, β
=
0.707, N
=
4, M
=
10, l
0
=
10
m, d
0
=
0.05
m, r
in
=
0.1
m, h
=
5
m, C
m
t
=
2.2
×
10
-
5
MP
a
-
1, and C
f
t
=
1
×
10
-
4
MP
a
-
1).


## 5. Conclusions

A pressure transient model for power-law fluids in porous media embedded with a tree-shaped fractal network has been developed and is expressed as a function of the power index of power-law fluids, tree-shaped branches number, diameter ratio, length ratio, and other parameters of porous media. This model is established and solved, type curves are illustrated, and dual fractal flow behavior characteristics are analyzed. Four flow regimes for pressure type curves can be established. Type curves are dominated by parameters of tree-shaped fractal network and porous media. These various parameters affect type curves differently. Power-law fluids flow in a porous media embedded with a tree-shaped fractal network is an interesting and challenging topic, and this work is in processing.



---

*Source: 405640-2014-02-24.xml*