Rev. Int. Mt. Num. Clc. Dis. Ing.

e
a
Vol. 20, 1, 7795 (2004)

Revista Internacional de
Mtodos Numricos para
e
e
Clculo y Diseo en Ingenier
a
n
a

Anlisis de estabilidad de un modelo

a
semi-impl
cito de un pozo petrol
fero

Gabriela B. Savioli, M. de los Angeles Morelli, Sebastin P. Regis y M. Susana Bidner

a
Laboratorio de Ingenier de Reservorios, Universidad de Buenos Aires
a
Pabelln de Industrias, Ciudad Universitaria
o
1428 Buenos Aires, Argentina
Tel.: 54-11-47 84 80 85; Fax: 54-11-47 80 01 45
e-mail: gsavioli@[Link]
e-mail: sbidner@[Link]

Resumen
Se desarrolla un modelo semi-impl
cito del ujo radial bifsico petrleo-gas hacia un pozo de hidrocarburos.
a
o
Para ello se modica en forma creativa el modelo black oil basado en el mtodo IMPES. En los simuladores
e
black oil tradicionales el pozo se modela como trmino fuente en las ecuaciones diferenciales aplicando la
e
hiptesis de ujo monofsico. Nuestro modelo, en cambio, representa las condiciones de produccin del pozo
o
a
o
como condiciones de contorno, considerando el ujo bifsico. Se estudian dos condiciones de produccin:
a
o
caudal total (petrleo y gas) constante y caudal de petrleo constante. Se comparan los perles de presiones
o
o
y saturaciones para ambas condiciones.
Asimismo, se analizan las inestabilidades numricas generadas por IMPES, obtenindose las restricciones
e
e
del incremento temporal. Para caudal total constante, el mximo incremento temporal que suministra
a
nimo, pero a partir de all aumenta con el tiempo de produccin.

o
estabilidad, tmax, disminuye hasta un m
o
Para caudal de petrleo constante, tmax disminuye a medida que el tiempo de produccin aumenta.
o

Palabras clave: ujo bifsico, IMPES, medios porosos, diferencias nitas, anlisis de estaa
a
bilidad, simulacin de reservorios.
o
STABILITY ANALYSIS OF A SEMI-IMPLICIT OIL WELL MODEL

Summary
A semi-implicit model of the radial ow of two phases (oil-gas) towards a well is developed. For that aim,
the black oil model, based on the IMPES method is creatively modied. Traditional black oil simulators
consider the well as a source term in the dierential equations, with the hypothesis of single-phase ow. Our
model, in contrast, considers well production conditions as boundary conditions, which take into account the
two-phase ow. Two production conditions are studied: constant total (oil and gas) ow rate and constant
oil ow rate. The pressure and saturation distributions for both conditions are compared.
Also, the numerical instabilities generated by IMPES are analyzed. The time increment restrictions to
maintain stability are obtained. For the constant total ow rate, the maximum time increment that provides
stability, tmax, decreases until a minimum, but from there it increases with production time. For the
constant oil ow rate, tmax diminishes as the production time increases.

Keywords: two-phase ow, IMPES, porous media, nite dierences, stability analysis, reservoir simulation.

c Universitat Polit`cnica de Catalunya (Espa a).

e
n

ISSN: 02131315

Recibido: Junio 2003

Aceptado: Julio 2003

78

G.B. Savioli, M.A. Morelli, S.P. Regis y M.S. Bidner

INTRODUCCION
o
e
Se adapta el modelo black oil15 a la simulacin numrica del ujo radial hacia un
pozo en recuperacin primaria. El mecanismo de recuperacin se debe exclusivamente a
o
o
la energ proporcionada por las fracciones livianas de hidrocarburos, inicialmente disueltas
a
en el l
quido, que se liberan y expanden al disminuir la presin durante la produccin. La
o
o
fase acuosa permanece inmvil, no participando del proceso productivo. Este mecanismo se
o
denomina solution gas drive en la jerga petrolera6 .
Se desarrolla un simulador bifsico petrleo-gas, cuyo modelo matemtico surge de
a
o
a
combinar las ecuaciones de Darcy para cada fase, las ecuaciones de continuidad para cada
componente y ecuaciones de estado emp
ricas. Estas ultimas se basan en la medicin en

o
el laboratorio de los parmetros PVT: factores de vol menes de las fases y solubilidad del
a
u
gas en el petrleo. Resulta un sistema de dos ecuaciones diferenciales parciales no lineales,
o
cuyas variables dependientes son la presin y la saturacin de la fase oleosa. El sistema
o
o
de ecuaciones se resuelve aplicando diferencias nitas y el mtodo IMPES15 (Impl
e
cito en
Presiones y Expl
cito en Saturaciones). Este modelo es un caso particular del ujo radial
trifsico (gas-petrleo-agua) descrito en Regis et al.7 , que resulta de mantener constante la
a
o
saturacin de agua y de anular su movilidad.
o
La descripcin del modelo hecha hasta aqu coincide con el modelo black-oil utilizado en
o

los simuladores comerciales. Se alaremos en el prximo prrafo las diferencias.

n
o
a
En los simuladores black-oil comerciales, cuyo objetivo es simular todo un reservorio, los
pozos se consideran como trminos fuente o sumidero puntuales en las ecuaciones diferene
ciales. Para modelarlos se recurre a ecuaciones monofsicas con propiedades constantes,
a
basadas en fuertes hiptesis simplicatorias. A n de obtener resultados ms cercanos a la
o
a
realidad, a dichos modelos de pozo simplicados se les agregan parmetros emp
a
ricos que
debe proveer el usuario y que se obtienen ajustando mediciones de produccin.
o
Nuestro inters, en cambio, radica en simular solamente un pozo. Justamente, cuando
e
se simula el ujo hacia un pozo, conviene trabajar en coordenadas cil
ndricas. En este
caso se ha supuesto simetr axial plantendose un modelo radial. Las condiciones de
a
a
operacin (ujo en el pozo productor) se modelan como condiciones de contorno, teniendo en
o
cuenta el ujo bifsico sin necesidad de recurrir a hiptesis simplicatorias. Adems, en este
a
o
a
simulador, la malla radial tiene una distribucin logar
o
tmica de los incrementos espaciales.
Esto permite visualizar mejor los cambios abruptos de presiones y saturaciones que ocurren
en los alrededores del pozo.
Se consideran dos modos de operacin del pozo representados por sendas condiciones de
o
contorno:
1. Caudal total (petrleo + gas) constante. Se emplea para hacer pronsticos de produccin
o
o
o
de petrleo y gas de un pozo. La condicin de contorno de caudal total constante se
o
o
corresponde a la condicin de extraccin con modernas bombas multifsicas8 .
o
o
a
2. Caudal de petrleo constante. Se utiliza en la interpretacin de ensayos transitorios
o
o
de presin en pozos (test draw-down). El objetivo principal de los ensayos es estimar
o
parmetros de la roca reservorio (permeabilidad, porosidad, factor de da o, etc.), ajusa
n
tando las mediciones de presin a los resultados de un modelo terico. Es decir, se aplica
o
o
un mtodo inverso. Los modelos tericos tradicionales son monofsicos9,10,11 . Nuestro
e
o
a
simulador permite remover las restricciones de dichos modelos tericos tradicionales, pues
o
contempla el ujo bifsico.
a
La eleccin del mtodo IMPES se basa en su rapidez de clculo, ya que desacopla y
o
e
a
linealiza el sistema de ecuaciones evaluando los coecientes en las presiones y saturaciones
correspondientes al tiempo anterior. El problema principal fue la necesidad de aplicarlo
utilizando muy peque os incrementos temporales para evitar inestabilidades.
n
La estabilidad del mtodo IMPES ha sido analizada para modelos multifsicos considee
a

79

Anlisis de estabilidad de un modelo semi-impl

a
cito de un pozo petrol
fero

rando uidos incompresibles de hasta tres dimensiones. Este anlisis se realiza sobre una
a
formulacin de las ecuaciones de tipo hiperblico4,12 .
o
o
Aqu se presenta un anlisis de estabilidad, teniendo en cuenta la compresibilidad de

a
los uidos e incluyendo las condiciones de contorno. Se encuentran las restricciones de
estabilidad para el ujo bifsico compresible y radial (de malla logar
a
tmica) y se determina
que la cota del incremento temporal se comporta en forma diferente seg n sea la condicin
u
o
de contorno elegida. No hemos visto que estas cotas hayan sido presentadas y analizadas
anteriormente. Ahora bien, la condicin general encontrada tiene como caso particular la
o
correspondiente a ujo bifsico incompresible y lineal (de malla regular), presentada por
a
Coats12 y otros autores4 .
TEOR
IA
Modelo matemtico
a
Para simular el ujo bifsico petrleo-gas hacia un pozo productor se adopta el modelo
a
o
black-oil. Se consideran dos fases: oleosa (o) y gaseosa (g) en presencia de una tercera
fase inmvil: el agua connata. En el sistema hidrocarbonado petrleo-gas existen dos
o
o
componentes: el componente petrleo (O), que se dene como el l
o
quido residual a presin
o
atmosfrica despus de una vaporizacin diferencial, y el componente gas (G), uido liberado
e
e
o
en esa vaporizacin6 . El componente gas se solubiliza en la fase oleosa, pero el componente
o
petrleo no se vaporiza en la fase gaseosa. La proporcin de componente gas en la fase
o
o
oleosa se estima usando los propiedades f
sicas de los uidos (parmetros PVT): los factores
a
de volumen Bo , Bg y la solubilidad del gas en el petrleo, Rso .
o
Se estudia el ujo bifsico hacia un pozo productor. Por lo tanto, como modelo geolgicoa
o
geomtrico, se considera un reservorio cil
e
ndrico con un unico pozo situado en el eje del

mismo. Todo lo cual se muestra en la Figura 1, donde h representa el espesor de la capa

productiva, rw el radio del pozo y re el radio mximo que el pozo es capaz de drenar. Si
a
se supone que todas las propiedades del reservorio son axisimtricas e independientes de la
e
profundidad y los efectos gravitatorios son despreciables, la unica coordenada relevante es

la radial.
r

q
h

Figura 1. Flujo radial hacia un pozo de hidrocarburos

Bajo estas hiptesis, las ecuaciones diferenciales del modelo black-oil se obtienen combio
nando la ecuacin de conservacin de masa para cada componente (G, O) con la ecuacin
o
o
o
,
emp
rica de movimiento de Darcy para cada fase (g, o)2,4 . As las ecuaciones resultan:
para el componente petrleo
o
1
r r

rk

kro po
Bo o r

So
Bo

(1)

para el componente gas

1
r r

rk

Rso kro po
krg pg
+
Bg g r
Bo o r

Rso So
Sg
+
Bg
Bo

(2)

80

G.B. Savioli, M.A. Morelli, S.P. Regis y M.S. Bidner

Las variables dependientes son las presiones pj y las saturaciones Sj . Los parmetros (datos)
a
son la porosidad, k la permeabilidad absoluta en la direccin radial, krj las permeabilidades
o
relativas y j las viscosidades. El sub
ndice j indica las dos fases (j = g, o).
La suma de las saturaciones es 1. Entonces, Sg = 1 So Swc (Swc saturacin de agua
o
connata). Adems, considerando nula la presin capilar gas-petrleo, pg = po . As resulta
a
o
o

un sistema de dos ecuaciones diferenciales (1)(2) con dos incgnitas po y So. Este modelo
o
tambin se puede obtener como caso particular del ujo radial trifsico (gas-petrleo-agua)7,
e
a
o
anulando la permeabilidad relativa al agua y manteniendo constante la saturacin de agua
o
en el valor Swc .
Las condiciones iniciales establecen presiones y saturaciones constantes al comienzo de
la explotacin
o
po (r, t = 0) = poIN

(3)

So (r, t = 0) = SoIN

(4)

Es interesante notar que este caso de ujo radial hacia un pozo tiene condiciones de contorno
particulares, ya que en el borde interno (r = rw ) se encuentra el pozo productor y en el
borde externo (r = re ) est el radio de drenaje. As los trminos fuente o sumidero de
a
,
e
las ecuaciones de conservacin se consideran nulos. Las condiciones de contorno en el pozo
o
productor (r = rw ) consideran el ujo bifsico con las dos posibilidades de extraccin ya
a
o
mencionadas:
1. caudal total de extraccin conocido en condiciones de reservorio (fondo del pozo)
o
Qt = Qo + Qg = 2rhk

kro krg
+
o
g

po
r

(5)
r=rw

En la ecuacin (5) se introdujeron las ecuaciones de Darcy para evaluar los caudales de
o
petrleo y de gas.
o
2. caudal de extraccin de petrleo conocido en condiciones de reservorio
o
o
Qo = 2rhk

kro po
o r

(6)
r=rw

La condicin de contorno en el borde externo (r = re ) impone ujo nulo para cada una
o
de las fases
Qj = 2rhk

krj pj
j r

=0

j = g, o

(7)

r=re

RESOLUCION NUMERICA APLICANDO EL METODO IMPES

El sistema de ecuaciones diferenciales (1)(2), con condiciones iniciales (3)(4) y condiciones de contorno (5) y (7) o (6) y (7), se resuelve numricamente con una tcnica en

e
e
e
diferencias nitas y aplicando el mtodo IMPES4 . La idea principal de este mtodo es
e
reemplazar la ecuacin (2) por una combinacin de las ecuaciones (1) y (2) en la que slo
o
o
o
aparezcan derivadas de la presin2,7 .
o
Para obtener dicha ecuacin en presin, el primer paso es desarrollar el miembro dereo
o
cho de las ecuaciones diferenciales, tratando de lograr una expresin en la que aparezcan
o
solamente las derivadas temporales de la presin y de las saturaciones. Para ello se debe
o

81

Anlisis de estabilidad de un modelo semi-impl

a
cito de un pozo petrol
fero

tener en cuenta que los parmetros PVT (Bg , Bo , Rso ) y la porosidad son funciones
a
conocidas de la presin (mediciones de laboratorio). As el miembro derecho de la ecuacin
o
,
o
(1) resulta

So
Bo

So d
So dBo
2
Bo dpo Bo dpo

So
+
Bo t

po
t

(8)

y el de la ecuacin (2)
o

Sg Rso So
+
Bg
Bo
+ Rso

Sg
+
Bg t

So
+
Bo t

Sg d
Sg dBg
2
Bg dpo Bg dpo

So d
So dBo
2
Bo dpo Bo dpo

po
t

po
+
t

So dRso po
+
Bo dpo t

(9)

En las ecuaciones (8) y (9) slo aparecen las derivadas temporales de la presin de la
o
o
fase oleosa y de las saturaciones.
El segundo paso del mtodo IMPES es deshacerse de las derivadas temporales de las satue
raciones. Para lograrlo se combinan las ecuaciones diferenciales multiplicando la ecuacin
o
del componente petrleo (1) por el trmino (Bo Rso Bg ) y la del componente gas (2) por
o
e
Bg . Sumando, reordenando y cancelando los trminos correspondientes, el miembro derecho
e
de la ecuacin combinada resulta
o
(Bo Rso Bg )

So
Bo

+ Bg

Rso So
Sg
+
Bg
Bo

po
po
= ct
= [cr (1 Swc ) + co So + cg Sg ]
t
t

(10)

En la ecuacin anterior se introduce la compresibilidad total, ct ,

o
ct = cr (1 Swc ) + co So + cg Sg

(11)

donde cr , co y cg son las compresibilidades de la formacin, del petrleo y del gas, denidas
o
o
como
1 d
dpo

(12)

Bg dRso
1 dBo

Bo dpo
Bo dpo

(13)

cr =
co =

cg =

1 dBg

Bg dpo

(14)

Con este procedimiento se obtiene la ecuacin en presin de la fase oleosa que se usa en
o
o
el mtodo IMPES:
e
(Bo Rso Bg )

1
r r

ro

po
r

+ Bg

po
1
po
r(g + Rso o )
= ct
r r
r
t

(15)

En la ecuacin (15), j = (kkrj /Bj j ) son las movilidades, j = g, o.

o
En las proximidades del pozo productor se requiere un estudio ms detallado porque
a
las velocidades de fase se van incrementando rpidamente. Por eso se adimensionaliza la
a

82

G.B. Savioli, M.A. Morelli, S.P. Regis y M.S. Bidner

coordenada radial con una transformacin logar

o
tmica. As se obtiene mayor densidad de

puntos de discretizacin en las proximidades del pozo productor usando una malla regular
o
x = ln

r
rw

= r = rw ex

(16)

De este modo, las derivadas radiales en las ecuaciones diferenciales (1) y (15) resultan
1
r r

rj

pj
r

1 2x
e
2
rw
x

pj
x

j = g, o

(17)

Discretizando la ecuacin (15) con una malla equiespaciada de bloques centrados

o
xi = (i 0, 5)x, i = 1, . . . , N , se obtiene el siguiente sistema de ecuaciones algebraicas
(Bo,i Rso,i Bg,i ) o,i+ 1 (pn+1 pn+1 ) o,i 1 (pn+1 pn+1 ) +
o,i+1
o,i
o,i
o,i1
2
2
+ Bg,i (g,i+ 1 + Rso,i+ 1 o,i+ 1 )(pn+1 pn+1 ) (g,i 1 + Rso,i 1 o,i 1 )(pn+1 pn+1 ) =
o,i+1
o,i
o,i
o,i1
2
2
2
2
2
2
= i ct,i (pn+1 pn )
o,i
o,i

2
x2 rw 2xi
e
t

i = 2, . . . , N 1

(18)

En los puntos de la malla i = 1 e i = N se agrega la discretizacin de las condiciones de

o
contorno7 .
Para calcular los coecientes de las ecuaciones (18) se deben estimar los parmetros
a
PVT, las viscosidades y las permeabilidades relativas en los puntos intermedios i 1/2.
Las permeabilidades relativas, que dependen de las saturaciones de las fases, se representan
mediante los siguientes modelos potenciales6
So Sor
1 Sor Sgc Swc

no

kro = kro

1 So Swc Sgc
1 Sor Sgc Swc

ng

krg = krg

(19)
(20)

o
o
o
tica
donde Sor es la saturacin de petrleo residual y Sgc la saturacin de gas cr 6 . Las
permeabilidades relativas se eval an en un punto upstream. En este caso resulta krj,i+ 1 =
u
2
krj (So,i+1 ), j = g, o. En cambio, se aplica un promedio aritmtico para aproximar las
e
propiedades que dependen exclusivamente de la presin: parmetros PVT y viscosidades.
o
a
El mtodo IMPES linealiza el sistema resultante de la discretizacin, ecuacin (18), evae
o
o
luando las movilidades, parmetros PVT, compresibilidades y porosidad en las saturaciones
a
y presiones correspondientes al tiempo anterior tn .
Con las nuevas presiones pn+1 obtenidas resolviendo el sistema (18) se calculan las
o
n+1
saturaciones de petrleo So
o
en forma expl 7 de la discretizacin de la ecuacin (1)
cita
o
o
y tambin evaluando las propiedades en el tiempo tn .
e
La ventaja del mtodo IMPES es que requiere menor tiempo de cmputo y capacidad
e
o
de almacenamiento que los mtodos totalmente impl
e
citos. Su principal inconveniente es
necesitar muy peque os incrementos en el paso del tiempo para evitar inestabilidades.
n
Esto se debe al clculo expl
a
cito de las saturaciones de petrleo. Para determinar la
o
eleccin adecuada del incremento temporal, en el siguiente apartado se realiza un anlisis
o
a
de estabilidad.

83

Anlisis de estabilidad de un modelo semi-impl

a
cito de un pozo petrol
fero

Anlisis terico de estabilidad

a
o
Para simplicar el anlisis de estabilidad, la ecuacin (1) se modica4 introduciendo el
a
o
caudal de petrleo Qo de acuerdo a la ecuacin (6). Luego se aplica el ujo fraccional fo ,
o
o
denido por
fo =

Qo
Qt

(21)

De tal modo, la ecuacin (1) se transforma en

o
1
r r

Qt
fo
2hBo

So
Bo

(22)

Aplicando las ecuaciones de Darcy y despreciando la presin capilar gas-petrleo, el ujo

o
o
fraccional se puede escribir en funcin de las permeabilidades relativas
o
fo =

Qo
kro /o
=
Qt
kro /o + krg /g

(23)

Adems, considerando constantes las viscosidades, el ujo fraccional var slo con la
a
a o
saturacin de petrleo, fo = fo (So).
o
o
La discretizacin de la ecuacin (22), aplicando el cambio de variables x = ln(r/rw ) y
o
o
tomando constante la porosidad, resulta

n+1
n
So,i = So,i

n+1
n+1
Bo,i
Bo,i t
1
+
n
2
Bo,i
x 2hrw e2xi

Qt fo
Bo

n
i+ 1
2

Qt fo
Bo

i = 1, . . . , N (24)
i 1
2

Para cada tiempo, se considera que la variacin espacial del caudal total es despreciable.
o
Esta hiptesis se cumple en forma aproximada, conrmndose por los resultados de la simuo
a
lacin. Con esta suposicin y evaluando las permeabilidades relativas y, en consecuencia, el
o
o
ujo fraccional de petrleo en un punto upstream fo,i+ 1 = fo (So,i+1 ), se llega a
o
2
n+1
n
So,i = So,i

n+1
n+1
Bo,i t
Bo,i
Qn
t
+
n
2
Bo,i
x hrw e2xi

n
n
fo (So,i+1 )
fo (So,i )

n
n
n
n
(Bo,i+1 + Bo,i ) (Bo,i + Bo,i1 )

(25)

En el Apndice se aplica el mtodo de las perturbaciones a la ecuacin (25) y se obtiene

e
e
o
la siguiente condicin de estabilidad sobre el incremento temporal t para cada tiempo tn
o

e2xi hr 2 (B n + B n )
t n
w
o,i
o,i1
(t ) < min
(26)
i
dfo
x
n
n
B
(S ) Qn
o,i

dSo

o,i

n
n
deniendo Bo,0 : = Bo,1
Esta condicin est calculada para ujo radial bifsico gas-petrleo, considerando uidos
o
a
a
o
compresibles y despreciando la variacin espacial del caudal total. La ecuacin (26) tiene
o
o
como caso particular la condicin de estabilidad del IMPES para ujo lineal, obtenida por
o
otros autores despreciando la transferencia de masa y la compresibilidad de los uidos4,12,13 .

84

G.B. Savioli, M.A. Morelli, S.P. Regis y M.S. Bidner

RESULTADOS NUMERICOS
Datos
El simulador desarrollado se aplica para analizar el ujo bifsico petrleo-gas hacia
a
o
un pozo productor para un conjunto de datos t
picos que se muestra en la Tabla I. Los
parmetros PVT (Bo , Bg , Rso ) y las viscosidades (0 , g ) del petrleo y del gas en funcin
a
o
o
de la presin pueden verse en las Figuras 2ae. Los datos fueron tomados de Odeh14 . Las
o
permeabilidades relativas al gas y al petrleo fueron calculadas con las correlaciones de
o
Fari a y Quintela15 , descritas en el libro de Bidner6 . Se les ajust un modelo potencial,
n
o

resultando kro = 0, 92, krg = 0, 272, no = 1, 654 y ng = 1, 960. Las permeabilidades relativas
al gas y al petrleo se representan grcamente en la Figura 2f en funcin de la saturacin
o
a
o
o
de petrleo. La compresibilidad de la formacin se considera nula.
o
o
1.8

0 .00 11

o [ Pas ]

RC

1.4

0 .00 08
0 .00 07
0 .00 06

1.2

0 .00 05

)a

B o [ mp etr le o

0 .00 09

)b

/ mp etr le o

SC

0 .0 01
1.6

0 .00 04
0

1000

2000

3000

40 00

5000

6000

10 00

Presin [ psi ]

2 00 0

30 0 0

40 00

5 00 0

6 00 0

Presin [ psi ]

0.8

g [ Pas ]

0.000 03

0.6

0.4

0.000 02

0.000 01

0.2

)c

)d

B g [ m gas

RC

/ m gas

SC

0.000 04

1000

2000

30 00

4000

5000

6000

1000

Pr esin [ psi ]

2000

3000

4000

5000

6000

P resin [ psi ]

250

1.0

200

0 .8

0.8

150

0 .6

0.6

100

0 .4

0.4

50

0 .2

0.2

k rg

SC

)f

R so [ m g as

k ro

/ mp etr leo

SC

1 .0

)e
0

0 .0

1000

20 00

3000

4000

Presin [ psi ]

5000

600 0

0.0
0.0

0.2

0 .4

0 .6

0.8

1.0

So

Figura 2. Propiedades de los uidos en funcin de la presin y permeabilidades relativas del

o
o
gas y del petrleo en funcin de la saturacin de petrleo (1 psi = 68,94 hPa)
o
o
o
o

Se comparan los perles de presin y saturacin en el pozo con las dos condiciones de
o
o
contorno. Para ello se consideran los siguientes caudales de extraccin total y de petrleo
o
o
en condiciones de fondo, Qt = Qo = 300 barriles/d = 0, 00055 m3 /s.
a

85

Anlisis de estabilidad de un modelo semi-impl

a
cito de un pozo petrol
fero

Espesor de la capa productiva

h = 30 pies = 9,144 m

Radio externo de drenaje

re = 767, 29 pies = 233,87 m

Radio del pozo

rw = 0, 328 pies = 0,1 m

Saturacin inicial de petrleo

o
o

SoIN = 0, 88

Saturacin residual de petrleo

o
o

Sor = 0, 4

Saturacin de agua connata

o

Swc = 0, 12

Saturacin de gas cr
o
tica

Sgc = 0, 01

Presin inicial del reservorio

o

poIN = 4500 psia = 31,02 MPa

Porosidad

= 0, 3

Permeabilidad absoluta

k = 100 mD = 9, 8 1014 m2

Presin de burbuja
o

pb = 4014, 7 psia = 27,67 MPa

Caudal de extraccin
o

Qt Qo = 300 bbl/d = 0, 00055 m3 /s

o

Tabla I. Propiedades de los uidos y caracter

sticas del reservorio

Anlisis del comportamiento del ujo

a
Al comenzar la explotacin del reservorio, la presin inicial es mayor que la presin de
o
o
o
burbuja, por lo que no hay gas libre en el reservorio. En consecuencia, el espacio poral est
a
ocupado por petrleo y por agua connata, cuyas saturaciones iniciales son: SoIN = 0, 88 y
o
Swc = 0, 12.
En la Figura 3 se representan dos perles de saturacin de petrleo en el pozo en
o
o
funcin del tiempo, cuando se extrae a caudal total y a caudal de petrleo constante de
o
o
300 barriles/d respectivamente. En ambos casos se observa que, mientras la presin en
a,
o
el pozo est por encima de la presin de burbuja, las dos curvas coinciden. La saturacin
e
o
o
de petrleo quedar constante e igual a la saturacin de petrleo inicial (So = 0, 88). En el
o
a
o
o
momento en que se alcanza la presin de burbuja, el gas disuelto en el petrleo comienza a
o
o
liberarse y se produce un descenso brusco de la saturacin de petrleo.
o
o

0,9
Qt=300 barriles/da
Qo=300 barriles/da

So

0,8

0,7

0,6

0,5
0

200

400

600

800

1000

1200

1400

tiempo [das]
Figura 3. Saturacin de petrleo en la celda que contiene al pozo en funcin del tiempo para
o
o
o
ambos reg
menes de extraccin
o

86

G.B. Savioli, M.A. Morelli, S.P. Regis y M.S. Bidner

Cuando el gas comienza a ser mvil (es decir, la saturacin de gas supera al valor cr
o
o
tico
Sgc ), aparece un cambio de pendiente en ambas curvas y So decrece ms suavemente. A
a
partir de ese momento, la condicin de produccin Qo = 300 barriles/d es de mucho
o
o
a
mayor exigencia que la de Qt = 300 barriles/d
a. En este ultimo caso, a medida que

el gas se libera y se mueve hacia el pozo, el caudal de gas producido aumenta y el de

petrleo disminuye. En consecuencia, las curvas se separan, declinando muy rpidamente
o
a
la saturacin correspondiente al caso de caudal de extraccin de petrleo constante.
o
o
o
En la Figura 4 se representa la presin de petrleo en el pozo en funcin del tiempo
o
o
o
cuando se extrae a caudal total y a caudal de petrleo de 300 barriles/d En ambos casos
o
a.
se observa un cambio brusco en la pendiente en el momento en que el gas comienza a ser
mvil.
o

Qt=300 barriles/da
4000

Qo=300 barriles/da

po [psi]

3000

2000

1000

0
0

200

400

600

800

1000

1200

1400

tiempo [das]

Figura 4. Presin de petrleo en la celda que contiene al pozo en funcin del tiempo para
o
o
o
ambos reg
menes de extraccin (1 psi = 68,94 hPa)
o

Del mismo modo que en la curva de saturacin (Figura 3), las presiones coinciden,
o
mientras el gas no es mvil, y se van separando a medida que aumenta la movilidad del
o
gas. En el caso de caudal de extraccin de petrleo constante, se llega a la presin de
o
o
o
abandono muy rpido, alrededor de los 1 400 d
a
as. En el caso de caudal total de extraccin
o
constante la presin declina ms suavemente.
o
a
Se eligi representar las distribuciones de saturaciones y de presiones en el pozo donde
o
se producen los cambios ms bruscos. Las distribuciones son semejantes, pero ms suaves,
a
a
a medida que nos alejamos del pozo. Adems, para un tiempo jo, tanto la presin como la
a
o
saturacin de petrleo aumentan con el radio7,8 .
o
o
A n de lograr distribuciones estables se utiliz x = 0, 3694 y t = 20 s.
o
Anlisis numrico de estabilidad
a
e
Se eval a el comportamiento del simulador desarrollado utilizando distintos valores del
u
incremento temporal t. Se analizan en forma independiente cada una de las condiciones
de contorno, pues se obtienen diferentes resultados.

87

Anlisis de estabilidad de un modelo semi-impl

a
cito de un pozo petrol
fero

Extraccin a caudal total constante

o
Se observ que al aumentar el valor de t, para x jo, la saturacin se inestabiliza
o
o
durante cierto tiempo, para luego volver a un comportamiento normal. Esto se ilustra en
la Figura 5, donde se representa de modo grco la saturacin de petrleo en la celda que
a
o
o
contiene al pozo, para t [s] = 100, 200 y 300, y se compara su comportamiento con un
caso estable correspondiente a t = 20 s. Para t = 100 s, la saturacin oscila entre los
o
350 y los 1 600 d para t = 200 s, entre los 200 y los 2 640 d y para t = 300 s,
as,
as;
entre los 140 y 3 360 d
as. Por su parte, las presiones, que se calculan en forma impl
cita,
no presentan ning n tipo de oscilacin, mantenindose siempre estables.
u
o
e

2
1,6

)a(

So

1,2

t = 100 s
t = 20 s

0,8
0,4
0
2

)b(

1,6

t = 200 s
t = 20 s

So

1,2
0,8
0,4
0
2

)c(

1,6

t = 300 s
t = 20 s

So

1,2
0,8
0,4
0
0

1000

2000
3000
tiempo [das]

4000

5000

Figura 5. Saturacin de petrleo en la celda que contiene al pozo cuando la extraccin es a

o
o
o
caudal total constante. Comparacin del perl estable, t = 20 s, con los perles
o
inestables para t [s] = 100 (a), 200 (b) y 300 (c)

88

G.B. Savioli, M.A. Morelli, S.P. Regis y M.S. Bidner

De esta observacin emp

o
rica se desprende que el incremento temporal debe ir disminuyendo durante un per
odo, para luego volver a aumentar. Estas observaciones se verican
analizando la cota de estabilidad, ecuacin (26). Debido a que el caudal total de extraccin
o
o
es constante, se estima Qn = Qt .
t
La ecuacin (26) puede compactarse considerando los valores del factor de volumen del
o
petrleo (Bo = volumen de petrleo en condiciones de reservorio/volumen de petrleo en
o
o
o

condiciones estndar). Este es mximo en el punto de burbuja Bo = Bob . Por debajo de

a
a
la presin de burbuja, su valor disminuye al disminuir la presin (Figura 2). Su m
o
o
nimo se
producir cuando hipotticamente la presin de abandono fuera la atmosfrica, all Bo = 1.
a
e
o
e

Adems, como e2x1 e2xi , es suciente que t(tn ) satisfaga

a
t(tn ) <

2
e2x1 hrw 2x
Bob Qt
max
i

1
dfo
n
(So,i )
dSo

cte

(27)

dfo
n
(So,i )
dSo

max
i

En la Figura 6a se muestra la derivada del ujo fraccional respecto de la saturacin

o
de petrleo. Si bien se ha representado para este conjunto particular de datos, su como
portamiento cualitativo es representativo para la mayor de los casos. Se puede observar
a

que dfo /dSo presenta un mximo en So 0, 77. Teniendo en cuenta que la saturacin de
a
o
petrleo va disminuyendo mientras avanza el tiempo de simulacin, se inere que la cota de
o
o
estabilidad de t cambiar as
a :
df
n

n
1. Inicialmente, cuando So,i So para todo i, el max dSo (So,i ) ir aumentando al
a
o
i
n
a
disminuir So,i . Por tanto, la cota de t ir disminuyendo.
n

2. En la etapa intermedia, en que existe alg n i tal que So,i > So y otro i tal que So,i < So ,
u
df
df
n

se puede tomar como max dSo (So,i ) al valor dSo (So ) y la cota de t es m
nima.
o
o
i
n

a
a
3. En la etapa nal, cuando So,i < So para todo i, el mximo ir decreciendo al disminuir
n
So,i y la cota de t vuelve a agrandarse.

Este es el esquema ptimo de variacin del paso temporal que debe incorporarse en el
o
o
simulador. El comportamiento de los valores de la cota de t, calculada a partir de la
ecuacin (26), se ilustra en la Figura 7.
o
1,2

6
o
Sd
fd
o

)b(

o
Sd
fd
o

f
5

1,0

)a(

0,8

0,6

0,4

0,2

0,0
0,6

0,7
So

0,8

0,5

o
*

0,4

0,9

0,4

roS

0,5

0,6

0,7

0,8

0,9

So

roS

Figura 6. a) Derivada del ujo fraccional en funcin de la saturacin de petrleo; b) Flujo

o
o
o
fraccional dividido la derivada del ujo fraccional en funcin de la saturacin de
o
o
petrleo
o

89

Anlisis de estabilidad de un modelo semi-impl

a
cito de un pozo petrol
fero

Cota de

t [seg]

400

300

200

100

0
0

1000

2000

3000

4000

5000

tiempo [das]

Figura 7. Cota del incremento temporal en funcin del tiempo de simulacin, cuando se
o
o
produce a caudal total constante de 300 barriles/d
a

Extraccin a caudal de petrleo constante

o
o
En este caso, al aumentar el valor de t, con x jo, la saturacin se torna ms
o
a
inestable a medida que el tiempo de simulacin aumenta. Este problema se ilustra en
o
la Figura 8, donde se representa la saturacin de petrleo en la celda que contiene al
o
o
pozo, para t [s] = 50, 100, 200 y 300, y se compara su comportamiento con un caso
estable correspondiente a t = 20 s. Para t = 50 s, la saturacin comienza a oscilar
o
aproximadamente a partir de los 560 d para t = 100 s, a partir de los 320 d para
as;
as;
t = 200 s a los 190 d
as; y para t = 300 s, a partir de los 120 d
as. Nuevamente,
las presiones, que se calculan en forma impl
cita, se mantienen estables hasta que valores
absurdos de saturacin provocan una detencin del simulador por errores de clculo.
o
o
a

)a

So

So

1,5

1,5

0,5

0,5
0

200

400
600
tiempo [das]

800

1000

200

400
tiempo [das]

600

2,5

)d

)c

2
So

t = 200 s
t = 20 s

1,5

t = 300 s
t = 20 s

2,5

So

t = 100 s
t = 20 s

)b

2,5

t = 50 s
t = 20 s

2,5

1,5

0,5

0,5
0

100

200
tiempo [das]

300

400

100
200
tiempo [das]

300

Figura 8. Saturacin de petrleo en la celda que contiene el pozo cuando la extraccin es a

o
o
o
caudal de petrleo constante. Comparacin del perl estable, t = 20 s, con los
o
o
perles inestables para t [s] = 50 (a), 100 (b), 200 (c) y 300 (d)

90

G.B. Savioli, M.A. Morelli, S.P. Regis y M.S. Bidner

As para simular un per

,
odo largo de tiempo se necesita un t cada vez ms peque o.
a
n
Esto se verica analizando la condicin de estabilidad, ecuacin (26), en funcin del tiempo.
o
o
o
n
En este caso, se puede aproximar Qn = Qo /fo(So,1 ). Entonces, t(tn ) debe vericar que
t

e2xi hr 2 (B n + B n )f (S n )x e2x2 hr 2 2x f (S n )
o
w
o,i
o,i1 o
o,1
o,2
w
t(tn ) < min

=
dfo
n
2iN
dfo

Qo
(So,2 )
Bn
(S n ) Qo
dSo
o,i

dSo

o,i

n
fo (So,2 )
= cte dfo n
(So,2 )
dSo

(28)
En la Figura 6b se muestra grcamente fo /(dfo /dSo ). Se puede observar que
a

a
fo /(dfo /dSo ) So Sor 0. Entonces, el incremento temporal debe tomarse cada vez ms
n
a
peque o, pues, a medida que avanza la explotacin, So,2 decrece, acercndose a la satun
o
racin residual de petrleo, Sor . El decrecimiento de los valores de la cota de t, calculada
o
o
a partir de la ecuacin (26), se ilustra en la Figura 9.
o

Cota de t [seg]

140
120
100
80
60
40
20
0
0

500

1000

1500

tiempo [das]

Figura 9. Cota de incremento temporal en funcin del tiempo de simulacin, cuando se

o
o
produce a caudal de petrleo de 300 barriles/d
o
a

F
sicamente, las inestabilidades se explican por la presencia de la fase gaseosa, mucho
ms mvil y compresible que el petrleo. La fase gaseosa aumenta con el tiempo en ambas
a
o
o
condiciones de contorno. Pero este incremento es mucho mayor cuando se exige producir un
caudal jo de petrleo que cuando el caudal total de ambas fases est controlado.
o
a
La cota de estabilidad, ecuacin (26), aplicada a la condicin de contorno de caudal de
o
o
petrleo constante, establece valores cada vez ms peque os de t al aumentar el tiempo de
o
a
n
simulacin. Afortunadamente, la principal aplicacin de este simulador con dicha condicin
o
o
o
de contorno es la interpretacin automtica de ensayos de pozos, que duran unos pocos d
o
a
as.
Esto permite incluir el simulador dentro de un mtodo inverso9 que estime los parmetros
e
a
ajustando los resultados del simulador a las mediciones de presin realizadas durante dichos
o
ensayos de pozo.

Anlisis de estabilidad de un modelo semi-impl

a
cito de un pozo petrol
fero

91

CONCLUSIONES Y APLICACIONES
Se describe y aplica un simulador black-oil al ujo radial bifsico compresible petrleoa
o
gas hacia un pozo, a n de analizar su comportamiento f
sico y su estabilidad numrica. Se
e
consideran dos condiciones de extraccin: caudal total y caudal de petrleo constante. Las
o
o
conclusiones son:
1. Con las dos condiciones de extraccin, los perles de presiones y saturaciones cambian
o
abruptamente cuando el gas comienza a ser mvil. Los perles coinciden para ambas
o
condiciones de extraccin cuando las presiones son mayores que la presin de burbuja.
o
o
Tambin coinciden a presiones menores que la de burbuja mientras el gas no sea mvil.
e
o
Una vez que la saturacin de gas super el valor cr
o
o
tico (Sg > Sgc ), las curvas comienzan
a separarse. La separacin aumenta a medida que disminuye la presin (ya que aumenta
o
o
la movilidad del gas) o, en otras palabras, a medida que avanza el tiempo de produccin.
o
Esto representa la realidad f
sica del proceso productivo cuyo mecanismo de drenaje se
debe a la expansin del gas.
o
2. Los perles de saturaciones son estables para incrementos temporales muy peque os. Las
n
inestabilidades se deben al clculo expl
a
cito de las saturaciones en el mtodo IMPES. Se
e
obtienen las siguientes restricciones de estabilidad en el paso de tiempo:
en el caso de caudal total de extraccin constante, la cota de t, para x jo, presenta
o
un m
nimo y luego crece. El algoritmo ptimo, que minimiza el tiempo de clculo,
o
a
utiliza este criterio para ir variando el incremento temporal, es decir, disminuyndolo
e
en la zona inestable y aumentndolo despus.
a
e
en el caso de caudal de petrleo constante, la cota de t para x, jo, tiende a 0 a
o
medida que el tiempo aumenta. El incremento t, debe ir disminuyendo a medida
que avanza la simulacin.
o
3. Las restricciones de estabilidad encontradas para ujo bifsico radial compresible (de
a
malla logar
tmica) incluyendo las condiciones de contorno son novedosas y tienen como
caso particular a la restriccin publicada para ujo bifsico lineal incompresible (de malla
o
a
regular).
4. La eleccin del mtodo IMPES se basa en su rapidez de clculo, ya que desacopla y
o
e
a
linealiza el sistema de ecuaciones evaluando los coecientes en las presiones y saturaciones
correspondientes al tiempo anterior. Esto es particularmente util cuando el simulador

se incorpora en la funcin objetivo de un mtodo inverso, ya que se requieren varias

o
e
ejecuciones del simulador para alcanzar la solucin optima. Justamente, es el caso de la
o
interpretacin de ensayos de pozo a caudal de petrleo constante. Con esta condicin
o
o
o
el incremento temporal disminuye con el tiempo. Sin embargo, dado que un ensayo
dura pocos d
as, el valor de t resulta sucientemente alto para que cada ejecucin
o
del simulador sea rpida. Por lo tanto, este simulador bifsico con la tcnica IMPES
a
a
e
aventaja a otros mtodos fuertemente impl
e
citos para ser incorporado en una rutina de
interpretacin automtica de ensayos transitorios de presin en pozos. Adems, supera
o
a
o
a
la restriccin de ujo monofsico de la interpretacin tradicional9,10,11 .
o
a
o
AGRADECIMIENTOS
Este trabajo fue realizado con subsidios de la Universidad de Buenos Aires y de la
Agencia Nacional de Promocin Cient
o
ca y Tecnolgica de la Argentina, PECOM Oil Co.
o
y la Universidad de Cuyo. M.S. Bidner, G.B. Savioli y M.A. Morelli son miembros del
Consejo Nacional de Investigaciones Cient
cas y Tcnicas de la Argentina (CONICET).
e

92

G.B. Savioli, M.A. Morelli, S.P. Regis y M.S. Bidner

NOMENCLATURA
S
mbolo

Denicin
o

B
c
cr
ct
fo
G
h
k
kr

kr
n
p
Qo
Qt
r
re
rw
Rso
S
Sgc
Sor
Swc
t
x

factor de volumen
compresibilidad
compresibilidad de la formacin
o
compresibilidad total
ujo fraccional
matriz de amplicacin
o
espesor de la formacin
o
permeabilidad absoluta
permeabilidad relativa
punto mximo de la curva de permeabilidad relativa
a
exponente de la curva de permeabilidad relativa
presin
o
caudal de petrleo en condiciones de reservorio
o
caudal total en condiciones de reservorio
radio
radio mximo de drenaje
a
radio del pozo
relacin gas-petrleo disuelta
o
o
saturacin
o
saturacin de gas cr
o
tica
saturacin residual de petrleo
o
o
saturacin de agua connata
o
tiempo
variable adimensional espacial
perturbacin
o
porosidad
autovalor de la matriz de amplicacin
o
movilidad
viscosidad

Sub
ndices
b
i
IN
g
o

relativo a presin de burbuja

o
i-sima malla espacial
e
inicial
gas
petrleo
o

Supra
ndices
n

n-simo paso del tiempo

e

Dimensin
o
Lt2 /M
Lt2 /M
Lt2 /M
L
L2

M/Lt2
L3 /t
L3 /t
L
L
L

L3 t/M
M/Lt

Anlisis de estabilidad de un modelo semi-impl

a
cito de un pozo petrol
fero

93

REFERENCIAS
1 J.R. Fanchi, Principles of applied reservoir simulation, Gulf Publishing Company, Houston,
Texas, [Link]., (1997).
2 J.R. Fanchi, K.J. Harpole y S.W. Bujnowski, BOAST- A three-dimensional, three-phase black
oil applied simulation tool, U.S. Department of Energy Report DOE/BC/10033-3, 2 volmenes,
u
(1982).
3 K. Coats, L. Thomas y R. Pierson, Compositional and black oil reservoir simulation, Society of
Petroleum Engineers Reservoir Evaluation and Engineering, Vol. 1, No 4, pp. 372379, (1998).
4 K. Aziz y A. Setttari, Petroleum reservoir simulation, Elsevier Applied Science Publishers,
Gran Bretaa, (1985).
n
5 M.S. Bidner y S.C. Gabbanelli, Simulacin numrica de procesos convencionales de recuperacin
o
e
o
de petrleo. Una revisin, Latin American Journal of Heat and Mass Transfer , Vol. 10, pp.
o
o
131154, (1986).
6 M.S. Bidner, Propiedades de la roca y los uidos en reservorios de petrleo, Eudeba (Editorial
o
de la Universidad de Buenos Aires), Argentina, (2001).
7 S. Regis, G. Savioli, M. Morelli y S. Bidner, Simulacin numrica del ujo radial trifsico
o
e
a
petrleo-gas-agua hacia un pozo, Mecnica Computacional , Vol. XIX, pp. 3136, (2000).
o
a
8 S.P. Regis, Desarrollo y aplicaciones de un simulador del ujo trifsico hacia un pozo de
a
petrleo, Tesis de Grado de Ingenier Qu
o
a
mica, Facultad de Ingenier Universidad de Buenos
a,
Aires, (1999).
9 G.B. Savioli y M.S. Bidner, Aplicacin del mtodo inverso al anlisis de ensayos de pozos
o
e
a
petrol
feros, Revista Internacional de Mtodos Numricos para Clculo y Dise o en Ingenier
e
e
a
n
a,
Vol. 10, No 1, pp. 321, (1994).
10 G.B. Savioli, P.M. Jacovkis y M.S. Bidner, Stability analysis and numerical simulation of 1-D
and 2-D radial ow towards an oil well, Computers & Mathematics with Applications, Vol. 33,
No 3, pp. 121135, (1997).
11 R.N. Horne, Modern well test analysis, Petroway, Inc., Palo Alto, CA, USA, (1990).
12 K.H. Coats, IMPES stability: the CFL limit, Society of Petroleum Engineers Paper 66345,
SPE Reservoir Simulation Symposium, Houston, Texas, 1114 February, (2001).
13 K.H. Coats, A Note on IMPES and some IMPES-based Simulation Models, Society of
Petroleum Engineers Journal , Vol. 5, No 3, pp. 245251, (2000).
14 A.S. Odeh, Comparison of solutions to a three-dimensional black-oil reservoir simulation problem, Journal of Petroleum Technology, Vol. 33, pp. 1325, (1981).
15 G.B. Savioli, CORPERM: una herramienta para la obtencin de curvas de permeabilidades
o
relativas de sistemas bifsicos, Programa y manual del usuario, publicados por el Laboratorio
a
de Ingenier de Reservorios, Universidad de Buenos Aires, (1997).
a
16 G. Marshall, Solucin numrica de ecuaciones diferenciales, tomo I y II, Editorial Revert,
o
e
e
Buenos Aires, Argentina, (1985).
17 M.R. Todd, P.M. ODell y G.J. Hirasaki, Methods for increased accuracy in numerical reservoir
simulators, Society of Petroleum Engineers Journal, pp. 515530, (1972).

94

G.B. Savioli, M.A. Morelli, S.P. Regis y M.S. Bidner

APENDICE
n+1
cita, ecuacin
o
El mtodo IMPES calcula las saturaciones de petrleo So en forma expl
e
o
(25), por eso el anlisis de estabilidad se efect a sobre estas ecuaciones. La ecuacin (25) es
a
u
o
estrictamente vlida para 2 i N 1. En los puntos i = 1 e i = N se deben incluir las
a
condiciones de contorno. Las condiciones de operacin en el pozo productor se consideran
o
en i = 1, donde se distinguen dos casos:

1. cuando el caudal total de extraccin es conocido, ecuacin (5), la ecuacin en saturacin

o
o
o
o
resulta
n+1
n
So,1 = So,1

n+1
n+1
Bo,1 t
Bo,1
Qn
t
+
n
2
Bo,1
x hrw e2x1

n
n
fo (So,2 )
fo (So,1 )

n
n
n
(Bo,2 + Bo,1 )
2Bo,1

(A.1)

2. cuando el caudal de extraccin de petrleo es conocido, ecuacin (6), la ecuacin en

o
o
o
o
saturacin resulta
o

n+1
n
So,1 = So,1

n+1
n+1
n+1
n
fo (So,2 )
Bo,1 t
Bo,1 t
Bo,1
Qo
Qn
1
t

+
n
2 e2x1 (B n + B n )
2 e2x1 B n
Bo,1
x hrw
x 2hrw
o,2
o,1
o,1

(A.2)

En el punto i = N , donde se incluye la condicin de contorno de ujo nulo en el borde

o
externo ecuacin (7), se obtiene
o
n+1
n
So,N = So,N

n+1
n+1
n
Bo,N
fo (So,N )
Bo,N t
Qn
t

n
n
n
2
Bo,N
x hrw e2xN (Bo,N + Bo,N 1 )

(A.3)

Entonces, se realiza un anlisis de estabilidad sobre las ecuaciones (25), (A.1), (A.2) y
a
(A.3) para determinar la eleccin adecuada del incremento temporal t. Observando dichas
o
ecuaciones, se deduce que la saturacin de petrleo en xi y en el tiempo tn+1 es funcin de
o
o
o
las saturaciones en xi y en xi+1 del tiempo anterior tn (excepto en i = N ), es decir
n+1
n
n
So,i = gi (So,i , So,i+1 )

1 i N 1

n+1
n
So,N = gN (So,N )

(A.4)

Para el estudio de la estabilidad se aplicar el mtodo de las perturbaciones16 para

a
e
una ecuacin en diferencias no lineal, teniendo en cuenta que la saturacin de petrleo
o
o
o
n+1
n
es una solucin de decaimiento (So,i So,i i). Este mtodo consiste en introducir una
o
e
o
perturbacin de la variable dependiente en un tiempo tn y calcular cmo se propaga dicha
o
n
perturbacin en tn+1 , linealizando las ecuaciones. Sean n , n las perturbaciones de So,i ,
o
i
i+1
n
So,i+1 , respectivamente, entonces
n+1
n
n
So,i + n+1 = gi (So,i + n , So,i+1 + n )
i
i
i+1

(A.5)

Desarrollando por Taylor y aproximando por el polinomio lineal, se obtiene la relacin

o
entre los errores en dos tiempos sucesivos
n+1 =
i

gi
gi
n
n
(So,i , So,i+1 )n +
(S n , S n )n
i
n
n
So,i
So,i+1 o,i o,i+1 i+1

(A.6)

As el vector de perturbaciones evaluado en tn+1 , n+1 = (n+1 , . . . , n+1 , . . . , n+1 ) se

,
1
i
N
relaciona con el del tiempo tn mediante la matriz de amplicacin G
o
n+1 = G n

(A.7)

95

Anlisis de estabilidad de un modelo semi-impl

a
cito de un pozo petrol
fero

La condicin de estabilidad en el caso de decaimiento es que la perturbacin en el tiempo

o
o
tn+1 sea menor que en el tiempo anterior. Una forma de asegurar esta condicin es que los
o
autovalores de G sean, en mdulo, menor que la unidad. Pero en este caso, el autovalor
o
dominante puede ser negativo llevando a una estabilidad oscilatoria. Por ello, Todd et al.17
requieren una condicin ms restrictiva, que llaman estabilidad no-oscilatoria, exigiendo
o
a
que los autovalores sean positivos. Coats13 recomienda para el IMPES utilizar la condicin
o
de estabilidad no-oscilatoria, que es la que aplicamos en este trabajo. De la ecuacin (A.6)
o
resulta que los autovalores de G son
i =

gi
(S n , S n )
n
So,i o,i o,i+1

(A.8)

Entonces, para el caso 1-caudal total de extraccin conocido en condiciones de reservorio

o
(ecuaciones (25), (A.1) y (A.3)), se tiene
dfo
n+1 n
n
n+1
Bo,i
gi
t Bo,i Qt dSo (So,i )
n
n
(So,i , So,i+1 ) = n
i =
n
n
n
2
So,i
Bo,i
x hrw e2xi (Bo,i + Bo,i1 )

1iN

(A.9)

n
n
deniendo Bo,0 : = Bo,1 n. Para el caso 2- caudal de extraccin de petrleo conocido en
o
o
condiciones de reservorio (ecuaciones (25), (A.2) y (A.3)), la frmula (A.9) es vlida para
o
a
2 i N y para i = 1 se cumple
n+1
Bo,1
1 = n
Bo,1

(A.10)

Para obtener condiciones de estabilidad no-oscilatoria, los autovalores deben cumplir

que 0 < i < 1. La restriccin i < 1 aplicada a la ecuacin (A.9) se cumple siempre,
o
o
n+1
n
o
pues dfo /dSo > 0 y Bo,i < Bo,i (por debajo de la presin de burbuja, donde se producen
los problemas de estabilidad por la aparicin del gas, el factor de volumen de petrleo
o
o
disminuye al disminuir la presin). Por la misma razn, 1 < 1 se cumple siempre en la
o
o
ecuacin (A.10). La restriccin i > 0 tambin se satisface siempre al aplicarla a la ecuacin
o
o
e
o
(A.10). As de estas condiciones no se obtiene ninguna restriccin sobre t. Por lo tanto, la
,
o
cota de estabilidad presentada en la ecuacin (26) del texto, surge directamente de aplicar
o
la restriccin i > 0 a la ecuacin (A.9).
o
o