Scribd
Upload
6 views

Lec7

The document discusses reservoir modeling and simulation using finite difference approximations. It presents an example involving a point-centered grid system in a horizontal reservoir, detailing the application of both implicit and explicit finite difference methods to solve a partial differential equation governing the flow system. The document concludes with calculations of pressure at specific points and prompts a discussion on the acceptability of the results and potential improvements.

Uploaded by

rya336437
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
6 views

Lec7

The document discusses reservoir modeling and simulation using finite difference approximations. It presents an example involving a point-centered grid system in a horizontal reservoir, detailing the application of both implicit and explicit finite difference methods to solve a partial differential equation governing the flow system. The document concludes with calculations of pressure at specific points and prompts a discussion on the acceptability of the results and potential improvements.

Uploaded by

rya336437
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 7

RESERVOIR MODELLING & SIMULATION

Application of Finite Difference Approximation


(Examples)

1
Example 1
A point-centered grid system in a horizontal reservoir is shown in Figure
1. The initial reservoir pressure is 3500 psia. The pressure at the
boundaries BC and CD drops to 3000 psia, while the pressure at the
boundaries AB and DA is maintained at the initial reservoir pressure.
The flow system is governed by the following partial differential equation
(PDE): 2 2
 P  P 1 P
2
 2 
x y  t

a) Base on the PDE, identify type of the flow system in this reservoir.
b) Apply the implicit finite difference approximation to this PDE.
c) Apply the explicit finite difference to this PDE.
d) Calculate the pressure at points 1, 2, 3, and 4 using the explicit
finite difference approximation, then discuss the results. Given:
x  y  1000 ft ,  5  10 5 ft 2 / Day , and t  5 Days

2
Example 1

3
Solution
a) Base on the PDE, identify type of the flow system in this reservoir.
 2 P  2 P 1 P
2
 2 
x y  t
1. 2-Dimensional flow
2. Homogeneous reservoir
3. Uniform thickness reservoir
4. Single-phase flow
5. Slightly compressible fluid

b) Apply the implicit finite difference approximation to this PDE.

Pi n1,1j  2 Pi ,nj1  Pi n1,1j Pi ,nj11  2 Pi ,nj1  Pi ,nj11 P n 1


1 i, j  P n
i, j
2
  ( )
( x ) ( y ) 2  t
c) Apply the explicit finite difference to this PDE.
n 1
Pi n1, j  2 Pi ,nj  Pi n1, j Pi ,nj 1  2 Pi ,nj  Pi ,nj 1 n
1 Pi , j  Pi , j
2
  ( )
( x ) ( y ) 2  t
4
Solution
a) Base on the PDE, identify type of the flow system in this reservoir.
 2 P  2 P 1 P
2
 2 
x y  t
1. 2-Dimensional flow
2. Homogeneous reservoir
3. Uniform thickness reservoir
4. Single-phase flow
5. Slightly compressible fluid

b) Apply the implicit finite difference approximation to this PDE.

Pi n1,1j  2 Pi ,nj1  Pi n1,1j Pi ,nj11  2 Pi ,nj1  Pi ,nj11 P n 1


1 i, j  P n
i, j
2
  ( )
( x ) ( y ) 2  t
c) Apply the explicit finite difference to this PDE.
n 1
Pi n1, j  2 Pi ,nj  Pi n1, j Pi ,nj 1  2 Pi ,nj  Pi ,nj 1 n
1 Pi , j  Pi , j
2
  ( )
( x ) ( y ) 2  t
5
d) Calculate the pressure at points 1, 2, 3, and 4 using the explicit finite
difference approximation, then discuss the results. Given:
x  y  1000 ft ,  5  10 6 ft 2 / Day , and t  5 Days
n 1
Pi n1, j  2 Pi ,nj  Pi n1, j Pi ,nj 1  2 Pi ,nj  Pi ,nj 1 n
1 Pi , j  Pi , j
2
  ( )
( x ) ( y ) 2  t
Multiply both sides of the finite difference equation by η ∆t, and
since ∆x = ∆y, thus the equation can re-arranged as follows:

 t
( x ) 2
P n
i 1, j  P n
i 1, j  P n
i , j 1  P n
i , j 1  4 P n
i, j   P n 1
i, j  P n
i, j

 t
P n 1
i, j P  n
i, j
( x ) 2
P n
i 1, j  P n
i 1, j  P n
i , j 1  P n
i , j 1  4 Pi, j 
n

6
 t
Pi ,nj1  Pi ,nj 
( x ) 2
Pi
n
1, j  Pi
n
1, j  Pi
n
, j 1  Pi
n
, j 1  4 Pi
n
, j

 t 5  10 5  5
  2 .5
x 2 (1000 ) 2

For P1 (i  1, j  1)
P1n 1  3500  2.5 3500  3500  3500  3500  4  3500   3500 psia
For P2 (i  2, j  1)
P2n 1  3500  2.5 3000  3500  3500  3500  4  3500   2250 psia
For P3 (i  1, j  2)
P3n 1  3500  2.5 3500  3500  3500  3000  4  3500   2250 psia
For P4 (i  2, j  2)
P4n 1  3500  2.5 3000  3500  3000  3500  4  3500   1000 psia
Quiz:
Are these results acceptable? If not, why? How to improve the results
to be acceptable?
7

You might also like