US2016115768A1PendingUtilityA1
Method for analyzing multi phase and heat flow of fluids in reservoir and recording media therefor
Assignee: KOREA INST GEOSCIENCE & MINERAPriority: Oct 24, 2014Filed: Oct 6, 2015Published: Apr 28, 2016
Est. expiryOct 24, 2034(~8.3 yrs left)· nominal 20-yr term from priority
E21B 41/0092E21B 47/00E21B 43/00
33
PatentIndex Score
0
Cited by
0
References
0
Claims
Abstract
A method for analyzing multi phase and heat flow in a reservoir is disclosed. The method comprises a stage for specifying fracture surfaces in a reservoir as multiple three-dimensional coordinates, and specifying a three-dimensional coordinates of a intersecting line which intersects between the fracture surfaces to grasp connectivity between the fracture surfaces, a stage for generating a mathematical equation having main variables of pressure, saturation and temperature in the cell from at least two or more of fluids among oil, water, gas that are moving within the fracture network and each flow model of fluids with heat, and a stage for grasping the flow of fluids and heat within the fracture network in the reservoir by activating the flow model by specifying initial values of pressure, saturation and temperature in each cell, and specifying production condition or injection condition of a drilling well connected to the fracture network.
Claims
exact text as granted — not AI-modified1 . A method for analyzing multi phase and heat flow of fluids in a reservoir, the method comprising:
(a) specifying a fracture surface in a reservoir as multiple three-dimensional coordinates and specifying a three-dimensional coordinate of an intersecting line between the fracture surfaces to grasp a connectivity between the fractures surfaces; (b) forming a fracture network by partitioning the fracture surface into multiple cells; (c) generating a flow model of each of at least two fluids among oil, water, and gas that move in the fracture network, and heat in a mathematical equation having main variables of pressure, saturation and temperature in the cell; and (d) grasping the flow of the fluids and heat within the fracture network in the reservoir layer by specifying initial values of the pressure, saturation and temperature in each of the cells, and specifying production condition or injection condition of a drilling well connected to the fracture network to operate the flow model.
2 . The method of claim 1 , wherein, when the fracture surface is partitioned into the multiple cells, the intersecting line between the fracture surfaces becomes a boundary line between the cells.
3 . The method of claim 1 , wherein the flow model of each of the fluids in the fracture network is a non-linear partial differential equation which is represented by a sum of flow amounts of the respective fluids moving along a longitudinal direction and a width direction of the cell within the cell.
4 . The method of claim 1 ,
wherein the fluid comprises oil and water, wherein the flow model of oil is represented by Equation 1 below, and wherein the flow model of water is represented by Equation 2 below.
∂
∂
L
[
T
oL
∂
∂
L
(
p
o
-
γ
o
z
)
]
Δ
L
+
∂
∂
W
[
T
oW
∂
∂
W
(
p
o
-
γ
o
z
)
]
Δ
W
+
Q
o
+
Q
Io
=
(
φ
V
f
)
∂
∂
t
(
S
o
B
o
)
Equation
1
∂
∂
L
[
T
wL
∂
∂
L
(
p
w
-
γ
w
z
)
]
Δ
L
+
∂
∂
W
[
T
wW
∂
∂
W
(
p
w
-
γ
w
z
)
]
Δ
W
+
Q
w
+
Q
Iw
=
(
φ
V
f
)
∂
∂
t
(
S
w
B
w
)
,
Equation
2
where L is a cell length, W is a cell width, T L is a longitudinal conductivity of fluid in cell, T W is width directional conductivity of fluid in cell, p is a fluid pressure in cell, γ is a cell density, z is a cell depth, Q is an amount of fluid produced in or injected from the drilling well, Q I is a flow amount of fluid on intersecting line between fracture surfaces, φ is a cell porosity, V f is a volume of fracture surface, t is time, S is a saturation rate of fluid in cell, B is a volume coefficient (volume) of fluid in cell, subscript o represents oil, and subscript w represents water.
5 . The method of claim 1 , wherein the fluid comprises gas, and
wherein the flow model of gas is represented by summing up the flow amounts of gas in state of gas and in state of gas being dissolved in oil.
6 . The method in claim 5 , wherein
the flow model of the gas is represented by Equation 3 below:
∂
∂
L
[
T
gL
∂
∂
L
(
p
g
-
γ
g
z
)
+
R
so
T
oL
∂
∂
L
(
p
o
-
γ
o
z
)
]
Δ
L
+
∂
∂
W
[
T
gW
∂
∂
W
(
p
g
-
γ
g
z
)
+
R
so
T
oW
∂
∂
W
(
p
o
-
γ
o
z
)
]
Δ
W
+
Q
g
+
R
so
Q
o
+
Q
Ig
+
R
so
Q
Io
=
(
φ
V
f
)
∂
∂
t
(
S
g
B
g
+
R
so
S
o
B
o
)
,
Equation
3
where L is a cell length, W is a cell width, T L is a longitudinal conductivity of fluid in cell, T W is a width directional conductivity in cell, p is a pressure of fluid in cell, γ is a cell density, z is a cell depth, R s is an amount of gas dissolved in another fluid, Q is an amount of fluid produced in or injected from the drilling well, Q I is a flow amount on the intersecting line between the fractures surfaces, φ is a cell porosity, V f is a volume of the fracture surface, t is time, S is a saturation rate of fluid in cell, B is a volume coefficient of fluid (volume) in cell, subscript o represents oil, and subscript g represents gas.
7 . The method of claim 1 , wherein the flow model of heat in the fracture network is represented by a sum of temperature variation amounts in a longitudinal direction and a width direction in the cell, and a value obtained by multiplying enthalpy by the flow models of the fluids.
8 . The method in claim 7 , wherein the flow model of heat is represented by Equation 4 below:
∂
∂
L
[
H
o
T
oL
∂
∂
L
(
p
o
-
γ
o
z
)
]
Δ
L
+
∂
∂
W
[
H
w
T
oW
∂
∂
W
(
p
o
-
γ
o
z
)
]
Δ
W
+
∂
∂
L
[
H
g
T
gL
∂
∂
L
(
p
g
-
γ
g
z
)
+
H
g
T
oL
∂
∂
L
(
p
o
-
γ
o
z
)
]
Δ
L
++
∂
∂
W
[
H
g
T
gW
∂
∂
W
(
p
g
-
γ
g
z
)
+
H
o
R
so
T
oW
∂
∂
W
(
p
o
-
γ
o
z
)
]
Δ
W
+
∂
∂
L
[
H
w
T
wL
∂
∂
L
(
p
w
-
γ
w
z
)
]
Δ
L
+
∂
∂
W
[
H
w
T
wW
∂
∂
W
(
p
w
-
γ
w
z
)
]
Δ
W
+
H
o
(
Q
o
+
Q
Io
)
+
H
g
[
Q
g
+
Q
Ig
+
R
so
(
Q
o
+
Q
Io
)
]
+
H
w
(
Q
w
+
Q
Iw
)
+
∂
∂
L
(
T
h
∂
θ
∂
L
)
+
∂
∂
W
(
T
h
∂
θ
∂
W
)
=
φ
V
f
∂
∂
t
(
U
)
Equation
4
where L is a cell length, W is a cell width, H is enthalpy of fluid in cell, T L is conductivity of fluid in a longitudinal direction in cell, T W is conductivity of fluid in a width direction in cell, p is a fluid pressure in cell, γ is a cell density, z is a cell depth, R s is an amount of gas dissolved in another fluid, Q is an amount of fluid produced in or injected from the drilling well, Q I is a flow amount of fluid on intersecting line between fracture surfaces, φ is a cell porosity, V f is a volume of fracture surface, t is time, S is a saturation rate of fluid in cell, B is a volume coefficient (volume) in cell, θ is a temperature in cell, U is internal energy, subscript o represents oil, subscript w represents water, subscript g represents gas, and subscript h represents heat.
9 . A computer-readable recording media in which the method of any one of claims 1 - 8 is implemented by a program to be executed by a computer.Join the waitlist — get patent alerts
Track US2016115768A1 — get alerts on status changes and closely related new filings.
We store only your email — no account needed. See our privacy policy.