Simulador Petroleo Agua en Fortran

1
Modelamiento de Flujo lineal 3D de Petrleo Agua en sistema

Mtodo IMPES para plantear el sistema de ecuaciones
Mtodo GAUSS SEIDEL para resolver el sistema de ecuaciones

Elaborado Por
ANDRES DAVID HERRERA PALACIO
LUIS ALEJANDRO ROBLEDO RUIZ
ARLEN ZAPTA BENITEZ

GRUPO 5

Presentado a
I.P ABEL NARANJO

UNIVERSIDAD NACIONAL DE COLOMBIA
SEDE MEDELLIN FACULTAD DE MINAS
SIMULACION DE YACIMIENTOS
FEBRERO 2012
2

TABLA DE CONTENIDO
Resumen
1. MODELO FSICO......4

2. MODELO MATEMTICO.5
2.1 Ecuaciones de difusividad para aceite y agua. ..5
2.2 Ecuaciones adicionales. ...5

3. MODELO NUMRICO.6
3.1 Anlisis para la ecuacin de difusividad del aceite. .6
3.1.1 Expansin del primer trmino en X.7
3.1.2 Expansin del segundo trmino en Y7
3.1.3 Expansin del segundo trmino en Z.8
3.2 Anlisis para la ecuacin de difusividad del agua. ..8
3.3 Calculo de transmisibilidades. ...8
3.3.1 Para el caso del aceite...8
3.3.2 Para la fase del agua.9

4. INTRODUCCIN AL MTODO DE SOLUCIN SECUENCIAL IMPIS (IMPLICIT PRESSURE
IMPLICIT SATURATION)11

5. SISTEMAS DE ECUACIONES PARA TODA LA MAYA. ..13
5.1 Cara superior Z=nz14
5.2 Cara del fondo Z=117
5.3 Caras intermedias 2< Z> nz.21

6. SOLUCIN DEL SISTEMA DE ECUACIONES. METODO LSOR..25

7. DIAGRAMA DE BLOQUES..26

8. MODELO COMPUTACIONAL. ..27
8.1 Algoritmo principal.27
8.2 Subrutinas Y Mdulos32
8.2.1 Modulo General32
8.2.2 Dimensionamiento...32
8.2.3 caracterizacin.33
8.2.4 Bloques Fantasmas..34
8.2.5 Clculo De Permeabilidades Relativas...35
8.2.6 Transmisibilidades36
3

8.2.7 Stenciles38

ANEXOS40
Anexo 1. Curvas de permeabilidad relativa dato de entrada.40
Anexo 2. Curva de presin capilar aceite-agua41
Anexo 3. Caracterizacin del sistema. (Fortran)..43

4

Informe de Simulacin de Yacimientos Flujo Petrleo Agua
Mtodo IMPIS para plantear el sistema de ecuaciones
Mtodo LSOR para resolver el sistema de ecuaciones

RESUMEN
Se procedi a realizar un simulador para un yacimiento bifsico, con una inclinacin de 5 en donde
tenemos constantes los caudales y las presiones son desconocidas, para nuestro caso flujo petrleo-
agua lineal en 3D con malla irregular y para las tres direcciones varan las permeabilidades, las
ecuaciones fueron planteadas por el mtodo IMPES y se resolvi por medio de mtodo GAUSS SEIDEL.

Las direcciones de los ejes sera:

Z

X

Y

El objetivo es desarrollar tericamente los pasos para el desarrollo del simulador y luego realizar el
algoritmo lgico de programacin en fortran, para encontrar del comportamiento de las presiones de
la fase de petrleo y las saturaciones en cualquier punto de un yacimiento para un tiempo de
determinado. Por los mtodos propuestos.

1. MODELO FSICO.

Para el sistema de flujo aceite-agua se asume que los componentes aceite y agua son incompresibles
e inmiscibles, de all que no haya transferencia de masa entre el agua y el aceite. Adicionalmente se
asume que el flujo es isotrmico y las fases estn en equilibrio termodinmico. Tambin se tiene
anisotropa en el sistema y el flujo es lineal en 3 dimensiones. Se eligi malla de bloque centrado.

Otras consideraciones son:
- El total de bloques en direccin x es de 4, en direccin y son 3, y en direccin z son 1.
- El tamao de paso en x es x= 100fts, en direccin y es y= 100fts, y en direccin z es z= 100fts.
- Hay inclinacin en la direccin X de 5.
- Asumimos como datos: ,
donde estas propiedades seran

constantes para todo los puntos del yacimiento.
- El comportamiento de la permeabilidad tendr un corpotamiento lineal con la distancia.
5

Figura 1. Sistema de bloque tridiagonal
2. MODELO MATEMTICO

La formulacin del modelo matemtico para flujo multifsico en yacimientos de petrleo consiste
bsicamente de las ecuaciones de flujo para todos los fluidos que componen el sistema, ecuaciones
adicionales que complementen la descripcin de los fluidos y las condiciones iniciales y de frontera.
Esto se realizar involucrando las ecuaciones de conservacin de masa, ecuaciones de estado y la ley
de Darcy. Las ecuaciones adicionales incluyen la saturacin de las fases y las presiones capilares en
funcin de la saturacin de las fases. Las condiciones iniciales y de frontera son necesarias para
resolver el modelo matemtico.

Expresin general de la ecuacin de difusividad

(

Para flujo lineal 3D el valor de es igual a 1, por tanto se tiene que las eecuaciones de difusividad
para el aceite y el agua:

2.1 Ecuaciones de difusividad para aceite y agua.

))
))
))
))
))
))

Donde
o
u
y
w
u
son los potenciales de flujo de las fases petrleo y agua respectivamente dados por:
6

w w w
P z u =
o o o
P z u =

En las ecuaciones de potencial
w
y
o
son los gradientes hidrostticos del agua y el petrleo y z es la
altura, la cual se considera positiva hacia abajo.

El uso de las ecuaciones de difusividad en funcin del potencial de flujo permite tener en cuenta el
flujo inclinado.

2.2 Ecuaciones adicionales.

Ecuaciones de saturacin:
1.0
o w
S S + =

Ecuacin de presin capilar:

Ecuacin de tasas:
o w t
q q q + =

( ) 1
w t w t t
f q f q q + =

/
/ /
rw w
w
ro o rw w
k
f
k k

=
+

Condiciones de lmite: Tipo Von Neumman, yacimiento cerrado.

Estas ecuaciones estn dadas para cada una de las fases, en nuestro caso para el petrleo y para el
agua.

7

3. MODELO NUMRICO

En esta etapa del trabajo se expanden las ecuaciones difusividad en diferencias finitas.

3.1 Anlisis para la ecuacin de difusividad del aceite.

))
))
))

3.1.1 Expansin del primer trmino en X:

))

Aplicando expansin central en Xi:

[(
]

Expandiendo las derivadas internas que dependen de i +1/2 y de i -

*(
) (
)+

[(
[(
]

3.1.2 Expansin del segundo trmino en Y:

))

Aplicando expansin central:
8

[(
]

Expandiendo las derivadas internas que dependen de j +1/2 y de j -

*(
) (
)+

[(
[(
]

3.1.3 Expansin del segundo trmino en Z:

))

Aplicando expansin central:

[(
]

Expandiendo las derivadas internas que dependen de k +1/2 y de k -

[(
) (
)]

*(
*(
+ (3)

3.2 Anlisis para la ecuacin de difusividad del agua.

Procediendo de manera similar al aceite para el agua se obtiene:

9

[(
[(
]

[(
[(
]

[(
[(
]

3.3 Calculo de transmisibilidades.

3.3.1 Para el caso del aceite:

Partiendo de la Definicin de transmisibilidad.

Usando promedios armnicos

(
((
((
) (
))

((
))
Invirtiendo

10

Similarmente:

(8)

(9)

(10)

(11)

(12)

3.3.2 Para la fase del agua

(13)

(14)

(15)

(16)

(17)

(18)

Habiendo definido las transmisibilidades podemos pasar a diferencias finitas:
- Para el agua
11

( ) ( )
p
w w w t w
V
T q S
t
A Au = A
A
(19)
- Para el aceite

( ) ( )
p
o o o t o
V
T q S
t
A Au = A
A
(20)

Las ecuaciones a resolver son la (19) y (20) y las incgnitas son los potenciales de flujo y las
saturaciones de agua y petrleo.

En cuanto a las tasas de produccin en un bloque dado debemos tener en cuenta lo siguiente:

En pozos de produccin:

Se da
, i j t
q q =
, , , , , i j w i j o i j
q q q = + . Ambos valores se desconocen pero deben cumplir con la ecuacin de tasas.

4. INTRODUCCIN AL MTODO DE SOLUCIN SECUENCIAL IMPIS (IMPLICIT PRESSURE IMPLICIT
SATURATION)

Las ecuaciones de flujo en forma de diferencias finitas, para las fases agua y petrleo
respectivamente, quedan como:

( ) ( )
p
w w w t w
V
T q S
t
A Au = A
A

( ) ( )
p
o o o t o
V
T q S
t
A Au = A
A

Las ecuaciones anteriores al tener en cuenta conceptos como potencial, presin capilar y ecuaciones
de saturacin se pueden transformar en las siguientes:
( ) ( )
p
w o c w w t w
V
T P P z q S
t
A A = A (

A
(21)

12

( ) ( )
p
o o o o t w
V
T P z q S
t
( A A = A

A
(22)

Sumando las dos expresiones anteriores se tiene:

( ) ( ) ( ) ( ) ( ) 0
w o o o w c w w o o
T P T P T P T z T z q A A +A A A A A A A A = (23)

En la ecuacin (23) hay trminos que dependen de la saturacin, las transmisividades y la presin
capilar, y si estos trminos se calculan a un valor ya conocido de saturacin, saturacin al tiempo n,
cuando se quiere calcular la presin de la fase petrleo al tiempo n+1 esta ecuacin se puede expresar
como

| | | | | | | |
o w c w w o o
T P q T P T z T z A A = +A A +A A +A A (24)
Donde,
w o
T T T = +

El sistema de ecuaciones tambin se puede dar en trminos de la presin para la fase agua o en
trminos de los potenciales de flujo de alguna de las fases.
Por ejemplo tomando como variable
o
u :

Combinando la definicin de potencial de flujo y la ecuacin de presin capilar, se puede establecer
para el potencial de flujo de la fase agua

w w w o c w
P z P P z u = =

Y restando y sumando
o
z
( ) ( )
w o o c w o o c
P z P z P z u = =u A (25)

Reemplazando
w
u en la ecuacin (19) se tiene
( ) ( )
p
w o c w t w
V
T P z q S
t
( A A u A = A

A
(26)
Sumando las ecuaciones (20) y (26) y teniendo en cuenta que
t w t o
S S A = A , se tiene

| | | | ( ) * 0
o w c w
T T P T z q A Au A A A A A ( =

(27)

Y la ecuacin a resolver es entonces

13

( )
1
o c
n n n
w w
T q T P T z
+
( ( A Au = + A + A A A (

(28)

Donde nuevamente T est dado como se mencion antes y
w o
A = .

En este caso para resolver el sistema de ecuaciones dado por las expresiones (24), (28) o la ecuacin
de transmisibilidades al tiempo (n+1) se requiere la saturacin tambin al tiempo (n+1) y puesto que
no se conoce se debe suponer. Con el valor supuesto de la saturacin se pueden calcular las
transmisividades y la presin capilar y, por lo tanto, los coeficientes y el trmino independiente del
sistema de ecuaciones. Se resuelve el sistema de ecuaciones y luego se calcula la saturacin. Luego se
compara la saturacin calculada con la supuesta y si no coinciden, dentro de la tolerancia establecida,
con los valores calculados se repite el procedimiento. Una vez se consiga coincidencia entre la
saturacin calculada y la supuesta se ha encontrado la saturacin al tiempo (n+1) y se procede luego a
obtener presin capilar y la presin de la otra fase al tiempo (n+1), para luego pasar al siguiente
tiempo y aplicar el mismo procedimiento.
En resumen los pasos del procedimiento IMPIS son:
1. Se expresan todas las ecuaciones de diferencias finitas en funcin de la presin (potencial) de una
de las fases.
2. Se eliminan los trminos de la saturacin aplicando la ecuacin de saturacin, de acuerdo con la
cual dS
w
= -dS
o
. De esta manera se obtiene el sistema de ecuaciones dado por la ecuacin (24) o
(28).
3. Se supone la saturacin al tiempo (n+1) y se calculan los coeficientes y el trmino independiente
del sistema de ecuaciones.
4. Se resuelve el sistema de ecuaciones obtenido en el paso anterior y se obtiene la presin
(potencial) de la fase en trminos de la cual se obtuvo el sistema de ecuaciones.
5. Con la presin de la fase conocida se aplica la ecuacin de diferencias finitas para dicha fase y se
tiene una ecuacin de saturacin que se resuelve en forma directa para obtener la saturacin de
dicha fase.
6. Se compara la saturacin calculada en el paso 5 con la supuesta en el paso 3 y si no coinciden se
toma el valor calculado como el nuevo valor supuesto y se regresa al paso 3 a repetir el
procedimiento.
7. Cuando se tenga coincidencia entre la saturacin supuesta y la calculada se conoce la saturacin
al tiempo (n+1).
8. Con la saturacin al tiempo n+1 se puede obtener la presin capilar al tiempo (n+1).
9. Con la presin capilar al tiempo (n+1) y la presin de una de las fases, tambin al tiempo (n+1), se
obtiene la presin de la otra fase.
10. Se repite el procedimiento del paso 3 en adelante para otros tiempos.

14

5. SISTEMAS DE ECUACIONES PARA TODA LA MALLA.

Las condiciones de lmite empleadas para la descripcin del sistema de ecuaciones es Von Neuman-
yacimiento cerrado. El sistema de ecuaciones se establece teniendo en cuenta el mtodo LSOR, el
cual recorre la malla por filas.

Partiendo de la ecuacin general para aceite en diferencias finitas:

En trminos de estncil:

m
k j oi k j i
m
k j oi k j i
m
k j oi k j i
m
k j oi k j i k j i k j oi k j k j oi k j i k j oi k j i
TC N BC S F E C W
1 , , , , , 1 , , ,
1
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , , + +
+
+
+
u u u u = u + u + u

Donde:
( )
n
i x o w k j
T T E
2 / 1 , , , +
+ =
( )
n
i x o w k j i
T T W
2 / 1 , , ,
+ =

( )
n
j y o w k j i
T T N
2 / 1 , , , +
+ =
( )
n
j y o w k j i
T T S
2 / 1 , , ,
+ =

( )
n
k z o w k j i
T T BC
2 / 1 , , ,
+ =

( )
n
k z o w k j i
T T TC
2 / 1 , , , +
+ =

C y F dependen de la ubicacin del bloque en la maya.

A continuacin se describe el sistema de ecuaciones para cada cara del sistema de bloques
tridimensional.

5.1 Cara superior z=nz

Sabiendo que los siguientes estncils son siempre igual a:
( )
n
i x o w k j
T T E
2 / 1 , , , +
+ =
( )
n
i x o w k j i
T T W
2 / 1 , , ,
+ =

( )
n
j y o w k j i
T T N
2 / 1 , , , +
+ =
( )
n
j y o w k j i
T T S
2 / 1 , , ,
+ =

( )
n
k z o w k j i
T T BC
2 / 1 , , ,
+ =

( )
n
k z o w k j i
T T TC
2 / 1 , , , +
+ =

15

Para (I,j,k)= (1,1,Nz)

Se anula:
0
. . , , , , 1 , , , 1 , , , 1 , ,
1
1 , ,
1
, 1 ,
*
, , 1
*
, ,
= = = = = = u = u = u = u
+
+
+
+
k j i k j i k j i
n
k j ci
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
k j oi k j oi k j oi
TC S W P P P P

La ecuacin correspondiente es:
m
k j oi k j i
m
k j oi k j i k j i k j oi k j k j oi k j i
N BC F E C
, 1 , , ,
1
1 , , , , , ,
*
, , 1 , ,
*
, , , , +
+
+
u u = u + u

Donde:
( ) ( ) ( ) | |
n
j y o w
n
i x o w
n
k z o w k j i
T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + =
| |
| | ( )
k j i
n
j y w k j i
n
i x w k j i
n
j y w
n
i x w
n
k z w k j i
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
j y w
n
i x w
n
k z w
n
k j ci
n
k z w k j i
z T z T z T T T z T
P T P T P T T T P T q F
, 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 ,
, 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , ,
+ + + + + +
+ + + + + +
+ + + + A +
+ + + + + =

Para (I,j,k)= (Nx,1,Nz)

Se anula:
0
. . , , , , 1 , , , 1 , , , 1 , ,
1
1 , ,
1
, 1 ,
*
, , 1
*
, ,
= = = = = = u = u = u = u
+ +
+
+
+
+ k j i k j i k j i
n
k j ci
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
TC E S P P P P

m
k j oi k j i
m
k j oi k j i k j i k j oi k j i k j oi k j i
N BC F C W
, 1 , , ,
1
1 , , , , , ,
*
, , , ,
*
, , 1 , , +
+

u u = u + u

Donde:
( ) ( ) ( ) | |
n
j y o w
n
k z o w
n
i x o w k j i
T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , , , +
+ + + + + =
| |
| | ( )
k j i
n
j y w k j i
n
j y w
n
k z w
n
i x w k j i
n
k z w k j i
n
i x w
n
k j ci
n
j y w
n
k j ci
n
j y w
n
k z w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
i x w k j i
z T z T T T z T z T
P T P T T T P T P T q F
, 1 , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , , 1 2 / 1 ,
, 1 , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , , 1 2 / 1 , , ,
+ + +
+ + +
+ + + + A +
+ + + + + =

Para (I,j,k)= (1,Ny,Nz)

Se anula:
0
. . , , , , 1 , , , 1 , , , 1 , ,
1
1 , , , 1 ,
*
, , 1
*
, ,
= = = = = = u = u = u = u
+ +
+
+ + k j i k j i k j i
n
k j ci
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
TC N W P P P P

1
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , ,
+
+
+
u u = u + u
m
k j oi k j i
m
BC S F E C

Donde:
16

( ) ( ) ( ) | |
n
i x o w
n
k z o w
n
j y o w k j i
T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , , , +
+ + + + + =

| |
| | ( )
k j i
n
i x w k j i
n
j y w
n
k z w
n
i x w k j i
n
k z w k j i
n
j y w
n
k j ci
n
i x w
n
k j ci
n
i x w
n
k z w
n
j y w
n
k j ci
n
k z w
n
k j ci
n
j y w k j i
z T z T T T z T z T
, , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 ,
, , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , ,
+ + + +
+ + +
+ + + + A +
+ + + + + =

Para (I,j,k)= (Nx,Ny,Nz)
Se anula:
0
. . , , , , 1 , , , 1 , , , 1 , ,
1
1 , , , 1 ,
*
, , 1
*
, ,
= = = = = = u = u = u = u
+ + +
+
+ + + k j i k j i k j i
n
k j ci
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
TC N E P P P P

Ecuacin correspondiente:
1
1 , , , ,
1
, 1 , , , , ,
*
, , , ,
*
, , 1 , ,
+
+

u u = u + u
m
k j oi k j i
m
BC S F C W

Donde:

( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
i x o w k j i
T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , , ,
+ + + + + =
| |
| | ( )
k j i
n
k z w
n
j y w
n
i x w k j i
n
k z w k j i
n
j y w k j i
n
i x w
n
k j ci
n
k z w
n
j y w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i
z T T T z T z T z T
P T T T P T P T P T q F
, , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,

+ + + + A +
+ + + + + =

Para (I,j,k)= (1<i<Nx,1,Nz)
Se anula:
0
. . , , 1 , , , 1 , , ,
1
1 , ,
1
, 1 ,
*
, ,
= = = = u = u = u
+
+
+
+
k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
k j oi k j oi
TC S P P P
m
k j oi k j i
m
N BC F E C W
, 1 , , ,
1
1 , , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , , +
+
+
u u = u + u + u

Donde:

( ) ( ) ( ) ( ) | |
n
j y o w
n
i x o w
n
k z o w
n
i x o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + + + =
| |
| | ( )
k j i
n
j y w k j i
n
i x w k j i
n
j y w
n
i x w
n
k z w
n
i x w k j i
n
k z w k j i
n
i x w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
j y w
n
i x w
n
k z w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
i x w k j i
z T z T z T T T T z T z T
P T P T P T T T T P T P T q F
, 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , , 1 2 / 1 ,
, 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , , 1 2 / 1 , , ,
+ + + + + +
+ + + + + +
+ + + + + + A +
+ + + + + + + =

Para (I,j,k)= (1<i<Nx,Ny,Nz)

Se anula:
0
. . , , 1 , , , 1 , , ,
1
1 , , , 1 ,
*
, ,
= = = = u = u = u
+ +
+
+ + k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
k j oi k j oi
TC N P P P

17

1
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , ,
+
+
+
u u = u + u + u
m
k j oi k j i
m
BC S F E C W

Donde:
( ) ( ) ( ) ( ) | |
n
i x o w
n
k z o w
n
j y o w
n
i x o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , +
+ + + + + + + =
| |
| | ( )
k j i
n
i x w k j i
n
i x w
n
k z w
n
j y w
n
i x w k j i
n
k z w k j i
n
j y w k j i
n
i x w
n
k j ci
n
i x w
n
k j ci
n
i x w
n
k z w
n
j y w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i
z T z T T T T z T z T z T
P T P T T T T P T P T P T q F
, , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
, , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + +
+ + +
+ + + + + + A +
+ + + + + + + =

Para (I,j,k)= (1,1<j<Ny,Nz)

Se anula:
0
. . , , 1 , , , , 1 , ,
1
1 , ,
*
, , 1
*
, ,
= = = = u = u = u
+
+
+ k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
TC W P P P

m
k j oi k j i
m
k j oi k j i
m
N BC S F E C
, 1 , , ,
1
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , , +
+
+
+
u u u = u + u

Donde:

( ) ( ) ( ) ( ) | |
n
j y o w
n
i x o w
n
k z o w
n
j y o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + + + =
| |
| | ( )
k j i
n
j y w k j i
n
i x w k j i
n
j y w
n
i x w
n
k z w
n
j y w k j i
n
k z w
n
k j i
n
j y w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
j y w
n
i x w
n
k z w
n
j y w
n
k j ci
n
k z w
n
k j ci
n
j y w k j i
, 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 ,
, 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , ,
+ + + + + +
+ + + + + +
+ + + + + + A +
+ + + + + + + =

Para (I,j,k)= (Nx,1<j<Ny,Nz)

Se anula:
0
. . , , 1 , , , , 1 , ,
1
1 , ,
*
, , 1
*
, ,
= = = = u = u = u
+ +
+
+ + k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
TC E P P P

m
k j oi k j i
m
k j oi k j i
m
N BC S F C W
, 1 , , ,
1
1 , , , ,
1
, 1 , , , , ,
*
, , , ,
*
, , 1 , , +
+
+

u u u = u + u

Donde:
( ) ( ) ( ) ( ) | |
n
j y o w
n
k z o w
n
j y o w
n
i x o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , +
+ + + + + + + =
| |
| | ( )
k j i
n
j y w k j i
n
j y w
n
k z w
n
j y w
n
i x w k j i
n
k z w
n
k j i
n
j y w k j i
n
i x w
n
k j ci
n
j y w
n
k j ci
n
j y w
n
k z w
n
j y w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i
, 1 , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
, 1 , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + +
+ + +
+ + + + + + A +
+ + + + + + + =

18

Para (I,j,k)= (1<i<Nx, 1<j<Ny,Nz)

Se anula:
0
. . 1 , , , ,
1
1 , ,
*
, ,
= = u = u
+
+
+ k j i
n
k j ci
n
k j ci
m
k j oi k j oi
TC P P

m
k j oi k j i
m
k j oi k j i
m
N BC S F E C W
, 1 , , ,
1
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , , +
+
+
+
u u u = u + u + u

Donde:
( ) ( ) ( ) ( ) ( ) | |
n
j y o w
n
i x o w
n
k z o w
n
i y o w
n
i x o w k j i
T T T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + + + + + =
| |
(
| | )
k j i
n
j y w k j i
n
i x w k j i
n
j y w
n
i x w
n
k z w
n
j y w
n
i x w
k j i
n
k z w
n
k j i
n
j y w k j i
n
i x w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
j y w
n
i x w
n
k z w
n
j y w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i
z T z T z T T T T T
z T z T z T P T P T
P T T T T T P T P T P T q F
, 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + + + + +
+ + + +
+ +
+ + + + + +
+ + A + + +
+ + + + + + + =

5.2 CARA DEL FONDO z= 1

Para (I,j,k)= (1,1,1):

Se anula:
0
. . , , , , 1 , , , 1 , , , 1 , ,
1
1 , ,
1
, 1 ,
*
, , 1
*
, ,
= = = = = = u = u = u = u

+
+
k j i k j i k j i
n
k j ci
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
BC S W P P P P

Educacin correspondiente:
m
k j oi k j i
m
TC N F E C
1 , , , , , 1 , , , , ,
*
, , 1 , ,
*
, , , , + + +
u u = u + u

Donde:
( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
i x o w k j i
T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , , , + + +
+ + + + + =
| |
| | ( )
k j i
n
k z w
n
j y w
n
i x w k j i
n
k z w k j i
n
j y w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
j y w
n
i x w k j i
z T T T z T z T z T
P T P T P T P T T T q F
, , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , , ,
+ + + + + + + + +
+ + + + + + + + +
+ + + + A +
+ + + + + =

Para (I,j,k)= (Nx,1,1)

Se anula:
0
, , . . , , 1 , , , 1 , , , 1 , ,
1
1 , ,
1
, 1 ,
*
, , 1
*
, ,
= = = = = = u = u = u = u
+
+
+
+ k j i k j i k j i
n
k j ci
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
E BC S P P P P

19

m
k j oi k j i
m
TC N F C W
1 , , , , , 1 , , , , ,
*
, , , ,
*
, , 1 , , + +
u u = u + u

Donde:

( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
i x o w k j i
T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + =
| |
| | ( )
1 , , 2 / 1 , , 1 , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + + + + +
+ + + + + +
+ + + + A +
+ + + + + =
k j i
n
k z w k j i
n
j y w k j i
n
k z w
n
j y w
n
i x w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
k z w
n
j y w
n
i x w
n
k j ci
n
i x w k j i
z T z T z T T T z T

Para (I,j,k)= (1,Ny,1)

Se anula:
0
. . , , , , 1 , , , 1 , , , 1 , ,
1
1 , , , 1 ,
*
, , 1
*
, ,
= = = = = = u = u = u = u
+
+
+ k j i k j i k j i
n
k j ci
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
BC N W P P P P

m
k j oi k j i
m
TC S F E C
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , , +
+
+
u u = u + u

Donde:
( ) ( ) ( ) | |
n
k z o w
n
i x o w
n
j y o w k j i
T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + =
| |
| | ( )
1 , , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 ,
1 , , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 , , ,
+ + + + + + +
+ + + + + +
+ + + + A +
+ + + + + =
k j i
n
k z w k j i
n
i x w k j i
n
k z w
n
j y w
n
i x w k j i
n
j y w
n
k j ci
n
k z w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
i x w
n
j y w
n
k j ci
n
j y w k j i
z T z T z T T T z T

Para (I,j,k)= (Nx,Ny,1)

Se anula:

0
, , , , . . 1 , , , 1 , , , 1 , ,
1
1 , , , 1 ,
*
, , 1
*
, ,
= = = = = = u = u = u = u
+ +
+
+ + k j i k j i k j i
n
k j ci
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
N E BC P P P P

m
k j oi k j i
m
TC S F C W
1 , , , ,
1
, 1 , , , , ,
*
, , , ,
*
, , 1 , , +
+

u u = u + u

Donde:

( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
i x o w k j i
T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , , , +
+ + + + + =
| |
| | ( )
1 , , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + +
+ + +
+ + + + A +
+ + + + + =
k j i
n
k z w k j i
n
k z w
n
j y w
n
i x w k j i
n
j y w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
k z w
n
j y w
n
i x w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i
z T z T T T z T z T

20

Para (I,j,k)= (1<i<Nx,1,1)

Se anula:
0
. . , , 1 , , , 1 , , ,
1
1 , ,
1
, 1 ,
*
, ,
= = = = u = u = u

+
+
k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
k j oi k j oi
BC S P P P

m
k j oi k j i
m
TC N F E C W
1 , , , , , 1 , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , , + + +
u u = u + u + u

Donde:
( ) ( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
i x o w
n
i x o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + + +
+ + + + + + + =
| |
| | ( )
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + + + + + + + +
+ + + + + + + + +
+ + + + + + A +
+ + + + + + + =
k j i
n
k z w k j i
n
j y w k j i
n
i x w k j i
n
k z w
n
j y w
n
i x w
n
i x w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
j y w
n
i x w
n
i x w
n
k j ci
n
i x w k j i
z T z T z T z T T T T z T
P T P T P T P T T T T P T q F

Para (I,j,k)= (1<i<Nx,Ny,1)

Se anula:
0
, , . . 1 , , , 1 , , ,
1
1 , , , 1 ,
*
, ,
= = = = u = u = u
+
+
+ k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
k j oi
m
k j oi k j oi
N BC P P P

m
k j oi k j i
m
TC S F E C W
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , , +
+
+
u u = u + u + u

Donde:
( ) ( ) ( ) ( ) | |
n
k z o w
n
i x o w
n
j y o w
n
i x o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + + + =
| |
| | ( )
1 , , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + + + + +
+ + + + + +
+ + + + + + A +
+ + + + + + + =
k j i
n
k z w k j i
n
i x w k j i
n
k z w
n
i x w
n
j y w
n
i x w k j i
n
j y w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
i x w
n
j y w
n
i x w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i

Para (I,j,k)= (1,1<j<Ny,1)

Se anula:
0
. . , , 1 , , , , 1 , ,
1
1 , ,
*
, , 1
*
, ,
= = = = u = u = u

+
k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
BC W P P P

m
k j oi k j i
m
k j oi k j i
m
TC N S F E C
1 , , , , , 1 , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , , + +
+
+
u u u = u + u

21

Donde:

( ) ( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
i x o w
n
j y o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + + +
+ + + + + + + =
| |
| | ( )
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 , , ,
+ + + + + + + + +
+ + + + + + + + +
+ + + + + + A +
+ + + + + + + =
k j i
n
k z w k j i
n
j y w k j i
n
i x w k j i
n
k z w
n
j y w
n
i x w
n
j y w
n
k j i
n
j y w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
j y w
n
i x w
n
j y w
n
k j ci
n
j y w k j i

Para (I,j,k)= (Nx,1<j<Ny,1)

Se anula:
0
, , . . 1 , , , , 1 , ,
1
1 , ,
*
, , 1
*
, ,
= = = = u = u = u
+
+
+ k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
E BC P P P

m
k j oi k j i
m
k j oi k j i
m
TC N S F C W
1 , , , , , 1 , , ,
1
, 1 , , , , ,
*
, , , ,
*
, , 1 , , + +
+

u u u = u + u

Donde:
( ) ( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
j y o w
n
i x o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + + + =
| |
| | ( )
1 , , 2 / 1 , , 1 , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + + + + +
+ + + + + +
+ + + + + + A +
+ + + + + + + =
k j i
n
k z w k j i
n
j y w k j i
n
k z w
n
j y w
n
j y w
n
i x w
n
k j i
n
j y w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
k z w
n
j y w
n
j y w
n
i x w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i

Para (I,j,k)= (1<i<Nx, 1<j<Ny ,1)

Se anula
0
. . 1 , , , ,
1
1 , ,
*
, ,
= = u = u

+
k j i
n
k j ci
n
k j ci
m
k j oi k j oi
BC P P

m
k j oi k j i
m
k j oi k j i
m
TC N S F E C W
1 , , , , , 1 , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , , + +
+
+
u u u = u + u + u

Donde:

( ) ( ) ( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
i x o w
n
j y o w
n
i x o w k j i
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + + +
+ + + + + + + + + =
| |
(
| | )
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 ,
, 1 , 2 / 1 , , , 1 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + + + + + + + +
+ + + + + +
+ + +
+ + + + + + +
+ A + + + +
+ + + + + + =
k j i
n
k z w k j i
n
j y w k j i
n
i x w k j i
n
k z w
n
j y w
n
i x w
n
j y w
n
i x w
n
k j i
n
j y w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
j y w
n
i x w
n
j y w
n
i x w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i
z T z T z T z T T T T T
z T z T P T P T P T
P T T T T T P T P T q F

22

5.3 CARAS INTERMEDIAS 1<z>nz

Para (I,j,k)= Para (1,1,1<k<Nz)

Se anula:
0
. . , , , 1 , , , 1 , ,
1
, 1 ,
*
, , 1
*
, ,
= = = = u = u = u

+
k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
S W P P P

m
k j oi k j i
m
k j oi k j i
m
TC N S F E C W
1 , , , , , 1 , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , , + +
+
+
u u u = u + u + u

Donde:
( ) ( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
i x o w
n
k z o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + + +
+ + + + + + + =
| |
| | ( )
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , ,
+ + + + + + + + +
+ + + + + + + + +
+ + + + + + A +
+ + + + + + + =
k j i
n
k z w k j i
n
j y w k j i
n
i x w k j i
n
k z w
n
j y w
n
i x w
n
k z w
n
k j i
n
k z w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
j y w
n
i x w
n
k z w
n
k j ci
n
k z w k j i

Para (I,j,k)= (Nx,1,1<k<Nz)

Se anula:
0
, , . . , 1 , , , 1 , ,
1
, 1 ,
*
, , 1
*
, ,
= = = = u = u = u
+
+
+ k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
E S P P P

m
k j oi k j i
m
k j oi k j i
m
TC N BC F C W
1 , , , , , 1 , , ,
1
1 , , , , , ,
*
, , , ,
*
, , 1 , , + +
+

u u u = u + u

Donde:
( ) ( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
k z o w
n
i x o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + + + =
| |
| | ( )
1 , , 2 / 1 , , 1 , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , , 1 2 / 1 , , ,
+ + + + + +
+ + + + + +
+ + + + + + A +
+ + + + + + + =
k j i
n
k z w k j i
n
j y w k j i
n
k z w
n
j y w
n
k z w
n
i x w
n
k j i
n
k z w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
k z w
n
j y w
n
k z w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
i x w k j i

Para (I,j,k)= (1,Ny, 1<k<Nz )

Se anula:
0
. . , , , 1 , , , 1 , , , 1 ,
*
, , 1
*
, ,
= = = = u = u = u
+ + k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
N W P P P

m
k j oi k j i
m
k j oi k j i
m
TC BC S F E C
1 , , , ,
1
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , , +
+
+
+
u u u = u + u
23

Donde:
( ) ( ) ( ) ( ) | |
n
k z o w
n
i x o w
n
k z o w
n
j y o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + + + =
| |
| | ( )
1 , , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 ,
1 , , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , ,
+ + + + + +
+ + + + + +
+ + + + + + A +
+ + + + + + + =
k j i
n
k z w k j i
n
i x w k j i
n
k z w
n
i x w
n
k z w
n
j y w
n
k j i
n
k z w k j i
n
j y w
n
k j ci
n
k z w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
i x w
n
k z w
n
j y w
n
k j ci
n
k z w
n
k j ci
n
j y w k j i

Para (I,j,k)= (Nx,Ny,1<k<Nz)

Se anula:
0
. . , , , 1 , , , 1 , , , 1 ,
*
, , 1
*
, ,
= = = = u = u = u
+ + + + k j i k j i
n
k j ci
n
k j ci
n
k j ci
m
N E P P P

m
k j oi k j i
m
k j oi k j i
m
TC BC S F C W
1 , , , ,
1
1 , , , ,
1
, 1 , , , , ,
*
, , , ,
*
, , 1 , , +
+
+

u u u = u + u

Donde:
( ) ( ) ( ) ( ) | |
n
k z o w
n
k z o w
n
j y o w
n
i x o w k j i
T T T T T T T T C
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , +
+ + + + + + + =
| |
| | ( )
1 , , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + +
+ + +
+ + + + + + A +
+ + + + + + + =
k j i
n
k z w k j i
n
k z w
n
k z w
n
j y w
n
i x w
n
k j i
n
k z w k j i
n
j y w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
k z w
n
k z w
n
j y w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i

Para (I,j,k)= (1<i<Nx,1,1<k<Nz)

Se anula:
0
. . , 1 , , ,
1
, 1 ,
*
, ,
= = u = u

+
k j i
n
k j ci
n
k j ci
m
k j oi k j oi
S P P

m
k j oi k j i
m
k j oi k j i
m
TC N BC F E C W
1 , , , , , 1 , , ,
1
1 , , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , , + +
+
+
u u u = u + u + u

Donde:
( ) ( ) ( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
i x o w
n
k z o w
n
i x o w k j i
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + + +
+ + + + + + + + + =
| |
| | (
)
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , , 1 2 / 1 , , ,
+ + + + + +
+ + +
+ + + + + +
+ + +
+ + +
+ + + + + A +
+ + +
+ + + + + + =
k j i
n
k z w k j i
n
j y w k j i
n
i x w
k j i
n
k z w
n
j y w
n
i x w
n
k z w
n
i x w
n
k j i
n
k z w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
j y w
n
i x w
n
k z w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
i x w k j i
z T z T z T
z T T T T T z T z T
P T P T P T
P T T T T T P T P T q F

24

Para (I,j,k)= (1<i<Nx,Ny,1<k<Nz)

Se anula:
0
. . , 1 , , ,
1
, 1 ,
*
, ,
= = u = u
+
+
+ k j i
n
k j ci
n
k j ci
m
k j oi k j oi
N P P

m
k j oi k j i
m
k j oi k j i
m
TC BC S F E C W
1 , , , ,
1
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , , +
+
+
+
u u u = u + u + u

Donde:
( ) ( ) ( ) ( ) ( ) | |
n
k z o w
n
i x o w
n
k z o w
n
j y o w
n
i x o w k j i
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + + + + + =
| |
| | (
)
1 , , 2 / 1 , , , 1 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + + +
+ +
+ + + + + +

+ +
+ + + + + + A +
+ + + + + +
+ + + =
k j i
n
k z w k j i
n
i x w
k j i
n
k z w
n
i x w
n
k z w
n
j y w
n
i x w
n
k j i
n
k z w k j i
n
j y w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
i x w
n
k z w
n
j y w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i
z T z T
z T T T T T z T z T z T
P T P T P T T T T T
P T P T P T q F

Para (I,j,k)= (1,1<j<Ny,1<k<Nz)

Se anula:
0
, , , , 1 , ,
*
, , 1
*
, ,
= = u = u
k j i
n
k j ci
n
k j ci k j oi k j oi
W P P

m
k j oi k j i
m
k j oi k j i
m
k j oi k j i
m
TC N BC S F E C
1 , , , , , 1 , , ,
1
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , , + +
+
+
+
u u u u = u + u

Donde:
( ) ( ) ( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
i x o w
n
k z o w
n
j y o w k j i
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + + +
+ + + + + + + + + =
| |
| | (
)
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 ,
, , 1 2 / 1 , , , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , ,
+ + + + + +
+ + + + + + +
+ + + + +
+ + +
+ + + + + A + + +
+ + + + + + + =
k j i
n
k z w k j i
n
j y w k j i
n
i x w
k j i
n
k z w
n
j y w
n
i x w
n
k z w
n
j y w
n
k j i
n
k z w k j i
n
j y w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
j y w
n
i x w
n
k z w
n
j y w
n
k j ci
n
k z w
n
k j ci
n
j y w k j i
z T z T z T
z T T T T T z T z T P T P T
P T P T T T T T P T P T q F

Para (I,j,k)= (Nx,1<j<Ny,1<k<Nz)

Se anula:
0
, , , , 1 , ,
*
, , 1
*
, ,
= = u = u
+ + k j i
n
k j ci
n
k j ci k j oi k j oi
E P P
25

m
k j oi k j i
m
k j oi k j i
m
k j oi k j i
m
TC N BC S F C W
1 , , , , , 1 , , ,
1
1 , , , ,
1
, 1 , , , , ,
*
, , , ,
*
, , 1 , , + +
+
+

u u u u = u + u

Donde:
( ) ( ) ( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
k z o w
n
j y o w
n
i x o w k j i
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + + + + + =
| |
| | (
)
1 , , 2 / 1 , , 1 , 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + + +
+ +
+ + + +
+ +
+ +
+ + + + + + A +
+ +
+ + + + + + + =
k j i
n
k z w k j i
n
j y w
k j i
n
k z w
n
i x w
n
k z w
n
j y w
n
i x w
n
k j i
n
k z w k j i
n
j y w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
k z w
n
j y w
n
k z w
n
j y w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i
z T z T
z T T T T T z T z T z T
P T P T
P T T T T T P T P T P T q F

Para (I,j,k)= (1<i<Nx,1<j<Ny,1<k<Nz)

m
k j oi k j i
m
k j oi k j i
m
k j oi k j i
m
TC N BC S F E C W
1 , , , , , 1 , , ,
1
1 , , , ,
1
, 1 , , , , ,
*
, , 1 , ,
*
, , , ,
*
, , 1 , , + +
+
+
+
u u u u = u + u + u

Donde:
( ) ( ) ( ) ( ) ( ) | |
n
k z o w
n
j y o w
n
k z o w
n
j y o w
n
i x o w k j i
2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , , , + +
+ + + + + + + + + =
| |
| | (
)
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 ,
, , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 2 / 1 , 1 , , 2 / 1 , , 1 , 2 / 1 , , , 1 2 / 1 , , ,
+ + + + + +
+ + +
+ + + + + +
+ + +
+ + +
+ + + + + + + A +
+ + +
+ + + + + + + + =
k j i
n
k z w k j i
n
j y w k j i
n
i x w
k j i
n
k z w
n
j y w
n
i x w
n
k z w
n
j y w
n
i x w
n
k j i
n
k z w k j i
n
j y w k j i
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w
n
k j ci
n
k z w
n
j y w
n
i x w
n
k z w
n
j y w
n
i x w
n
k j ci
n
k z w
n
k j ci
n
j y w
n
k j ci
n
i x w k j i
z T z T z T
z T T T T T T z T z T z T
P T P T P T
P T T T T T T P T P T P T q F

6. SOLUCIN DEL SISTEMA DE ECUACIONES. METODO LSOR

El sistema de ecuaciones heptadiagonal generado por el recorrido de los bloques por filas ser
resuelto usando el mtodo LSOR, descrito por la siguiente ecuacin:

Sistema de ecuaciones:

Las (
*
P
) se hallaran por Thomas (eliminacin de Gauss).
El procedimiento para resolver el problema por LSOR es el siguiente:
26

1. Se ingresan los datos de entrada, los cuales son: el coeficiente de relajacin para el mtodo
LSOR, el tiempo a simular (en meses), el tamao del paso
t A
(das), y propiedades del fluido
como la API.
2. Se hallan los valores de las presiones intermedias (
*
P
) haciendo el recorrido por filas, de la
ecuacin anterior resulta una matriz tridiagonal la cual se soluciona mediante el algoritmo de
Thomas con la siguiente secuencia:
3. Ahora, con las presiones intermedias (
*
P
) halladas en el paso anterior, se da inicio al proceso
iterativo dado por la siguiente ecuacin:
m
k j i k j i
m
k j i
wP P w P
, ,
*
, ,
1
, ,
) 1 ( + =
+
donde se busca
encontrar las presiones en cada bloque de la malla durante el tiempo a simular.
Para el proceso iterativo se tiene el siguiente procedimiento:
1. Suponer
( ) M
k j i
P
, ,
,
j i ,
en este caso se asumen igual a la presin inicial.

2. Calcular
( ) 1
,
+ M
j i
P de la ecuacin general de LSOR.
m
k j i k j i
m
k j i
wP P w P
, ,
*
, ,
1
, ,
) 1 ( + =
+

3. Si
1
, ,
+
~
M
j i
M
j i
P P
,
j i ,
es decir
tolerancia P P
M
j i
M
j i
s
+
,
1
,
, se detiene el proceso iterativo.
4. Si
1
, ,
+
=
M
j i
M
j i
P P
,
j i ,
, entonces:
1
, ,
+
=
M
j i
M
j i
P P
,
j i ,

Repetir proceso para los dems tiempos, hasta llegar al tiempo deseado.

7. DIAGRAMA DE BLOQUES
A continuacin se describen pasos secuenciales para dar solucin al sistema numrico propuesto.
27

8. MODELO COMPUTACIONAL ALGORITMO EN FORTRAN.
8.1 Algoritmo principal
Diagrama de Bloques
Llamar los archivos Call
archivos
Suponer Sw
Calcular Kr con las curvas y
ecuaciones correspondientes
Calcular transmisibilidades del
aceite y del agua
Definicin de Stencils para
presin implcita y saturacin
implcita
Resolver por Thomas la matriz
pentadiagonal de presiones P*
Aplicar LSOR calular P
M+1
y
comparar con la P supuesta
Se calcula saturacin para el
agua con ecuacin proveniente
de diferencias finitas
Comparar Sw calculada con Sw
supuesta
No converge
Hacer Swcalculada =
Sw supuesta
Repetir desde el paso
2 hasta lograr
convergencia.
Converge
Sw=Sw
n+1

Calcular Pc
n+1
Po
n+1
Pw
n+1

28

29

30

31

32

33

8.2 Subrutinas Y Mdulos
8.2.1 Modulo General

8.2.2 Dimensionamiento
34

8.2.3 caracterizacin

35

8.2.4 Bloques Fantasmas

36

8.2.5 Clculo De Permeabilidades Relativas

37

8.2.6 Transmisibilidades

38

39

8.2.7 Stenciles

40

41

ANEXOS
Anexo 1. Curvas de permeabilidad relativa dato de entrada

Sw Krw Kro
0,18 0 1
0,21 0 0,92
0,24 0,00002 0,85
0,27 0,00014 0,79
0,3 0,00045 0,71
0,33 0,00111 0,64
0,36 0,00232 0,57
0,39 0,0043 0,51
0,42 0,00733 0,45
0,45 0,01175 0,40
0,48 0,01791 0,34
0,51 0,02623 0,29
0,54 0,03714 0,25
0,57 0,0516 0,21
0,6 0,06882 0,17
0,63 0,09069 0,14
0,66 0,117441 0,11
0,69 0,14963 0,087
0,72 0,18807 0,066
0,75 0,23347 0,048
0,78 0,28664 0,033
Krw = 2,225sw
4
- 1,622sw
3
+ 0,452sw
2
- 0,058sw + 0,003
R = 1
Kro = -2,522sw
4
+ 5,906sw
3
- 2,915sw
2
- 1,874sw + 1,400
R = 1
0
0.2
0.4
0.6
0.8
1
1.2
0 0.2 0.4 0.6 0.8 1
Krw
Kro
Poly. (Krw)
Poly. (Kro)
42

0,81 0,34842 0,021
0,84 0,41968 0,013
0,87 0,50135 0,007
0,9 0,59439 0,003

Anexo 2. Curva de presin capilar aceite-agua

Sw Pcow
0,2 8
0,25 4,3
0,3 3
0,4 1,78
0,5 1,21
0,6 0,79
0,7 0,43
0,8 0,1
0,9 0

Pcow= -853,3Sw
5
+ 2625,sw
4
- 3137,Sw
3
+ 1820,Sw
2
-
517,6sw + 59,78
R = 0,997
0
1
2
3
4
5
6
7
8
9
0 0.2 0.4 0.6 0.8 1
Curva de presin capilar o/w
Pcow Poly. (Pcow)
P( psi)
Saturacin.
43

Anexo 3. Caracterizacin del sistema. (Fortran)

44

Anexo 4. Dimensionamiento del modelo matemtico.

Simulador Petroleo Agua en Fortran

Cargado por

Copyright:

Formatos disponibles

Simulador Petroleo Agua en Fortran

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Simulador Petroleo Agua en Fortran

Cargado por

Copyright:

Formatos disponibles

1

Modelamiento de Flujo lineal 3D de Petrleo Agua en sistema

donde estas propiedades seran

en este caso se asumen igual a la presin inicial.

También podría gustarte

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.