_ K

p(. m ) ii es un factor de agotamiento y me indica el estado de presion del yacimiento y

J-1" Bo
(I - R2) es un factor de caida de presion que me indica el diferencial de presion que se esta
aelicando al yacimiento para tener la tasa de produccion q; en este factor R es la relacion Pwf
/ P
Si la ecuacion (1 .102) se aplica ala tasa de flujo de dos intervalos consecutivos y se tomara la
relacion entre elias se tendria

q ,:: , ~P" [jl~~J , (I - R' ), '

(1.103)
q, P, (1-J?2)~-
[ Km )
J-1()B" J

EI lade derecho de la ecuacion (1 .103) se puede ir evaluando y como se conoce la tasa de flujo
del intervalo anterior al que se esta calculando 0 sea qj-1 se puede calcular qj y suponiendo que
el intervalo de presion LlP i la tasa de flujo se mantuvo constante en q y se calculo LlNp)i ' el
tiempo requerido para producir este volumen , para que la presion cayera de Pi-1 a Pi ' se
obtiene de

I:!N,,) ,
~/) ,
q,

y el tiempo total transcurrido hasta que la presion haya caido a Pi es

,
I , = L8.I,
/= 1

y de esta forma se puede ir transformando la escala de Np en una escala de tiempo .

1.6 - Prediccion de Yacimientos de Condensado 0 Petr61eos Volatiles.

En los metodos tradicionales de prediccion de yacimientos (Pirson , Muskat, etc) se hace un

balance de materiales de tipo volumetrico y se supone que el gas que se va liberando del petroleo
permanecera como gas dentro 0 fuera del yacimiento . Esto es cierto siempre y cuando el gas que
se libera no sufra comportamiento retrogrado en el yacimiento, 10 cual ocurre normalmente con
yacimientos de petr61eos medianos 0 pesados 0 yacimientos de gas humedo 0 seco.

Cuando se trata de yacimientos de petroleo volatil, lIamados tambien de petr61eo claro 0 cercanos
al punto critico, y yacimientos de condensado , se presentara condensacion retrograda en el
yacimiento y parte de 10 que era gas a una determinada presion pasara a Ifquido a una presion
menor.

La ocurrencia de condensacion retrograda en el caso de un yacimiento de petroleo volatil puede

implicar que la saturacion de liquido despues de un decremento de presion aumente 0 se
mantenga constante a pesar de haberse producido un cierto volumen de petroleo, yen el caso de
un yacimiento de condensado significa que se presenta formacion de saturacion de hidroearburos
liquidos y esto desafortunadamente implica perdida de los componentes mas pesados y mas
valiosos del gas, pues no habra flujo de Ifquido hasta que no se tenga la saturaci6n critiea .

56
Si no se tiene en cuenta la condensaci6n retrograda en un yacimiento de petr61eo se tendra como
resultado de la predicci6n un factor de recobro mas bajo de 10 que se puede tener, pues durante el
proceso de condensacion retrograda se esta suponiendo un valor de saturacion menor que el
real , 10 cual implica mayores valores de Kg/Ko y final mente menor capacidad de produccion; y en
el caso de un yacimiento de condensado, si se permite la ocurrencia de la condensacion
retr6grada, se presentara la perdida de los componentes mas valiosos del gas y una reduccion en
el rendimiento de /iquidos.

Para evitar la condensaci6n retr6gada en un yacimiento de condensado se acostumbra hacerle al

yacimiento mantenimiento de presi6n para evitar que la presi6n del yacimiento lIegue a la presion
de r~cio , esto se puede hacer con inyeccion de agua 0 de gas; en el caso de inyeccion de gas , el
gas que se inyecta es gas residual proveniente del mismo yacimiento que al lIegar a superficie es
procesado para removerle sus componentes pesados y el residuo, que es basicamente meta no,
se inyecta luego al yacimiento y de esta manera se previene la condensacion retrograda por dos
razones pues por una parte se presuriza el yacimiento y por otra el gas residual va modificando la
composici6n del fluido del yacimiento haciendolo cada vez mas liviano y por tanto reduciendo la
zona bifasica del diagrama P-T del fluido del yacimiento desplazando la linea de puntos de rocio
hacia la izquierda ; puede ocurrir que el desplazamiento de la linea de puntos de rocio hacia la
izquierda 0 la reduccion de la zona bifasica del diagrama P-T sea tal que la isoterma del
yacimiento ya no entre a la zona bifasica y p~r tanto no habra condensaci6n retr6grada , en este
caso se puede suspender la inyecci6n de gas y continuar la despresurizacion del yacimiento sin
peligro de que se presente la condensacion de los elementos mas valiosos del fluido en el
yacimiento .

Para tener en cuenta la condensacion retrograda en la predicci6n del yacimiento se debe hacer un
balance de materiales de tipo composicional y suponer que en el yacimiento se esta presentando
una despresurizacion a volumen constante (CVO)

En general el procedimiento se puede plantear de la siguiente manera:

i) Se toma una mol del sistema a las condiciones iniciales del yacimiento y al que se Ie conoce la
composici6n . Se calcula el volumen de esta mol y este sera el volumen que se va a considerar
constante; es el volumen poroso donde esta almacenado el hidrocarburo. Para lIevar este
volumen a condiciones normales se debe ademas calcular Bo, Rs, y Bg a las condiciones
iniciales del yacimiento.

Ii) Se considera una carda de presion la cual debe ser baja ( unas 50 LPC maximo) y a la nueva
presi6n se hace un calculo de fases tomando 1 mol y la composicion inicial del sistema . Este
calculo de fases permite conocer las moles de gas y vapor, su composicion y sus volUmenes.

EI volumen total del sistema ahara es

(1 .104)

y la diferencia entre este volumen y el volumen inicial del sistema es el volumen que tendra que
salir del yacimiento, ya que el volumen de este permanece constante.

Se debe tam bien calcular a esta presi6n los factores volumetricos Bo, Rs , Bg , !J<" 119

EI volumen drenado sera entonces:

Vo=Vs-V (1 .105)

iii) Cuando se trata de yacimientos de condensado el volumen drenado sera todo gas, mientras el
condensado no alcance la saturaci6n critica . En el caso de yacimientos' de petr61eo volatil 0
condensado con saturaci6n de /iquido ya por encima de la saturaci6n critica, el volumen
drenado estara compuesto de fase vapor y fase liquida y es necesario conocer cuanto sale de
cada fase. Para ello se procede de la siguiente manera .

57
 Se supone un cierto valor VOd para el liquido drenado, 0 sea que la saturacion de liquido sera :

SI VI - V"" = So (1 .106)
V
donde VI es el volumen de liquido obtenido en el calculo de fases.

Con la saturacion de petroleo se puede obtener Kg/Ko .

Con los valores Kg/Ko se obtiene la relacion gas petroleo instantanea de

K e * fie!. * !!:.~ + R.\. (1 .107)

R= K B . 11K
" i'

La relacion gas - petroleo instantanea. R, tambien debe ser igual a :

V,,!! + Vno * R,. (V, - Vf)() ) + (v/)()2 * Rs

R=
B. ~
V
s" S J.: S ()
(1 .108)
lJo V /)()

So So

De la ecuacion (1 .108) se puede despejar Voo Y compararlo con el valor supuesto. Si son
iglJales se puede conocer los volumenes de cada fase que han salido del yacimiento y los que
permanecen , al igual que sus composiciones, y por 10 tanto la composicion del sistema de
hidrocarburos que permanece en el yacimiento. Si no son iguales se asume un nuevo valor
pC:lra Voo , que puede ser el que se acaba de calcular 0 el promedio de el supuesto y el
caiculado , y se repiten los cinco pasos anteriores. Se continua este procedimiento hasta que en
u" a iteracion el valor supuesto sea igual al calculado.

iv) Lma vez encontrados el volumen de petroleo y gas removidos durante el decremento de
presion se procede a encontrar la composicion del sistema remanente de la siguiente manera:

~ 10les de gas remanentes (n gr):

ngr=ng-ngO (1109)

ng=n*V (1 .110)

p * V /:/) p * (V, - V()n )

n ed = Z / ' RT (1 .111)
' Z v RT

donde ng y ngO son las moles de gas obtenidas en el calculo de fases, y las moles de gas
removidas respectivamente, V es la fraccion molar del sistema que es gas de acuerdo con el
calculo de fases y Zv es el factor de compresibilidad para la fase vapor.

Moles de Uquido Remanentes (N Lr ):

nLR=nL-nLD (1 .112)

nL=n*L (1 .113)
P * Von
l1 /.d= Z RT (1 .114)
I.

58
donde nL YnLD son las moles de IIquido del calculo de fases y las drenadas respectivamente, n
son las moles iniciales del sistema sobre las que se realizo el calculo de fases, L es la fraccion
molar del sistema que es liquido de acuerdo con el calculo de fases y ZI es el factor de
compresibilidad para la fase liquida

Conociendo las moles remanentes de cada fase se puede obtener la composicion del sistema
remanente asi:

(1.115)
n~r + n,.r
v) Se asume un nuevo decremento de presion y se repiten los pasos ii) - iv) teniendo en cuenta
que el calculo de fases se realiza sobre las moles remanentes obtenidas en el decremento
anterior y cuya composicion se ha obtenido con la ecuacion (1.108).

vi) EI procedimiento continua hasta que se Iiegue a la presion de abandono preestablecida hasta
la cual se quiere realizar la prediccion.

vii) La informacion obtenida se puede presentar finalmente de muchas formas; por ejemplo en el
caso de la prediccion de un yacimiento, a cada presion se puede obtener Np/N, RGP
acumulada y R. En el caso de un condensado se presenta la cantidad de gas producida y el
rendimiento de liquido ademas del contenido de fraccion pesada en el sistema remanente.

La Figura 9 muestra el diagrama de bloques para este tipo de aplicacion del calculo de equilibrio
liquido - vapor.

Otra forma de hacer la prediccion del comportamiento de un yacimiento de condensado es usando

los resultados de una prueba de laboratorio conocida como prueba CVD (Constant Volume
Depletion) en la cual se toma en una celda PVT una muestra de volumen V y composicion Zi una
presion Po igual a la presion inicial del yacimiento. Cuando la muestra se despresuriza hasta una
presion P, < Poa esta presion ocupara un volumen igual a V+tN, =V, y habra condensacion de un
volumen L, de liquido, tal como 10 muestra la Figura 10; este incremento en volumen es 10 que se
expande el fluido del yacimiento al ocurrir la disminucion de presion y es el mismo volumen que
tendra que salir del yacimiento pues su capacidad de almacenamiento es V. EI liquido L, no saldra
del yacimiento y por tanto la cantidad tN, que sale del yacimiento es todo gas. EI gas que se retira
de la celda PVT se hace manteniendo la presion P, y cuando se haya removido el gas la celda
tendra nuevamente un volumen V de hidrocarburos, compuesto por (V, - t; V, - L,) unidades de
volumen de gas y L, unidades de volumen de condensado, y estara lista para la despresurizacion
desde P, hasta P2 . AI gas que se retira en cada despresurizacion se determina su composicion.

AI final de la prueba se puede elaborar una tabla como la tabla 1.

La columna t; Vi (CY) da el volumen de gas removido de la celda en cada despresurizacion,

medido a condiciones de la celda (CY) y la columna t; Vi (CN) da el mismo volumen pero medido a
condiciones normales; la columna Zi es la relacion entre el volumen de gas removido en cada
despresurizacion medido a condiciones de la celda, que sera el volumen real, y ese mismo
volumen de gas a condiciones de la celda calculado usando la ecuacion de estado del gas ideal,
para calcular este ultimo volumen se requiere conocer el numero de moleS de gas removidas en
cada despresurizacion y esto se puede hacer dividiendo la columna t;Vi (CN) por 379 ..

59
Tabla 1-. Resultados de una Prueba CVD.

Presion Composicion del Gas tNi Li tNi Zi

C1 C2 C3 C4 Cn (CY) (CY) (CN)
Po ---­ --­ ---­ ---­ --­ 0 0 0 Zo
P1 --­ --­ --­ --­ --- tN i L1 --­
P2 --­ --­ --­ --­ --­ t:,V2 L2 --­ Z2

Pn --­ --­ --­ --­ --­ t:,Vn Ln --- Zn

La Figura 10 muestra las despresurizaciones desde Po a P 1 y desde P1a P2.

Con la informacion de la prueba CVD dada en la tabla anterior se puede proceder a hacer la
predicci6n del yacimiento tomando una unidad de volumen de yacimiento, por ejemplo un acre.pie
3
que es igual a 43560 pies , y elaborando una tabla como la que se muestra en la tabla 2.

Tabla 2-. Calculos para la Prediccion de un Yacimiento de Condensado Usando los

Resultados de un Prueba CVD.

(1 ) (2) -(3) (41 (5) (6) (7) (8) (9) (10) ( 11)

P t:, V, ,,£t:,V i t:,V ir ,,£t:,V ir Vu "£ Vu Ri (~,) (~'J

LV! / .1

(PCN) (PCN) (PCN) (PCN) (Bls.) (Bls.) (PCN/BN) N

(%) (%)
(%L
-

Para elaborar la tabla anterior se hace 10 siguiente:

EI volumen de gas contenido en un acre-pie de yacimiento es

V If =43560*';'*(I-S
If' w)

La relacion entre el volumen de gas contenido en un acre-pie de yacimiento y el volumen de la

celda donde se realize la prueba CVD es

43560*¢ * (I- S w )
r= (1.116)
V ('c1J a

Por 10 tanto la columna (2) de la tabla anterior es la columna de (t:,VilcN de la tabla de

resultados de la prueba CVD multiplicada por r. La columna (3) es la sumatoria de la columna
(2) y es el acumulado de gas humedo en PCN que se han recuperado del yacimiento a la
presion dada .

EI gas residual es 10 que queda despues de retirarle en superficie al gas producido sus
componentes pesados, 0 sea despues de removerle al gas humedo producido su GPM . Por 10
tanto la columna (4) se obtiene de :

( ~V) 7' F *Y
~VIr = I CN *(1- L -'- ' ) (1 .117)
1=1

60
donde F; es la fracci6n del componente i que se puede retirar del gas humedo.

La columna (5) es la sumatoria de la columna (4) .

Suponiendo que el liquido condensado en el yacimiento permanece en el yacimiento , el

volumen de liquido que aparece en la columna (6) es el removido del gas que aparece en la
columna (2) 0 sea

(L.lAVI ) = (~ Vi tv * L.
f F,Y,MW, (
Bl )
s. ( 1. 118)
. , 379*5.6\5 ,=) p"

La columna (7) es la sumatoria de la columna (6) .

La columna (8) se obtiene de dividir la columna (4) por la columna (6) .

Entrar Datos
~,Tc" PCi , Wi: i=1,n
p, T, Pa , 6P

1.
Reallzar Calcu lo de Fases a
PyT
En con trar el vol umen que
- -. - - .. - . : _ l _ _ • _ ... _
-'

~
C alcu lar B o v R _a P v

G~ L -_ _ _. - -_ _- - - '

Real izar calculo de lases a P

Enc ontrar VL , Vg Y Vs

Calc ular So y Calcul ar

KglK"

Figura 9. Diagrama de Bloques para los Calculos de una prediccion de Yacimientos.

61
 Caleular R de

*~*~*R s
R= K,
Ko JJ, B, ~ <~

C aleul ar Voo de
Vs
VOll = --oBo-----''---­
-.!.(R-Rs)+l
B.

NO
VO O)C3IC= Voo)sup

SI

Obtener Moles Remanentes Liquido

y Va por 0)

Impn m ir P Y cantidades Produei das de

Iiquido y Vapor

p >p. ~-
? SI
.10
~)

NO

Figura 9. Continua cion.

Para obtener las columnas (9) y (10) se debe obtener G , y esto se puede hacer con

G = r *V,ddaI S )!, ( p eN) ( 1.119)

Una vez conocido G, la columna (9) es la columna (3) dividido por G y la (10) es la columna (5)
dividida por G .

Para calcular la columna (11) se requiere Ny este se puede obtener de

O
N = G * 7· F, yi MW
379 * 5.6 15 ~ , (Sis.) (1.120)
,=.. P"

62
donde Y,O es la fracci6n molar del gas inicialmente en el yacimiento .

Calculado N la columna (11) se obtiene de dividir la columna (7) por N.

PI r­
P=P i ~

Ir=Z=

~ ;r'I l

t
vI VLI

i
~

T'
P=P 1 Vgd2
P2

1 If V,
. _. _ .- -- . _-- -- --­

\T +
VL VL2

Figura 10-. Prueba CVD

63
1.7 -.Comportamiento de Acuiferos

Se conoce como acuifero el yacimiento con saturacion 100% de agua que' rodea un yacimiento de
hidrocarburos, EI acuifero puede tener contacto con el yacimiento a traves del limite inferior de este
y en este caso se habla de acuifero de fondo, 0 a traves de su limite externo en la direccion radialy
en este caso se habla de acuifero lateral 0 de curia ,

AI disminuir la presion del yacimiento esta caida de presion IIegara eventualmente al contacto de
este con eJ acuifero y empezara a pasar agua de este ultimo al primero esto es 10 que se conoce
como empuje hidraulico y es un mecanismo de conservacion de la presion del yacimiento pues los
fluidos que salen de este son reemplazados por el agua procedente del acuifero,

EI agua normal mente tiene mucha mas movilidad que el petroleo y trata de irrumpir a traves del
yacimienlo ocasionando en el caso de los acuiferos de fondo 10 que se conoce como conificacion
del agua y en el caso de los laterales el fenomeno conocido como digitacion , En ambos casos se
tiene reduccion de la capacidad de flujo del petroleo pero el problema mas critico se puede
presentar con la conificacion pues rapidamente se puede tener presencia de agua en el pozo
bloqueando la salida del aceite y adem as fluyendo hacia el pozo

Los acuiferos, tambien se pueden clasificar de acuerdo al flujo que se presenta en ellos y en este
caso se habla de acuiferos lineales y radiales ; los primeros seran los de fondo y los laterales seran
radiales cuando contactan el yacimiento a traves de los 360 0 0 una fraccion apreciable de la
circunferencia , cuando esta fraccion es pequeria su comportam iento se aproxima mas al caso
lineal. En los dos casos el acuifero puede ser finito 0 infinito , Un acuifero se considera infinito
cuando su tamario es mucho mayor que el tamario del yacimiento y por tanto una perturbacion de
presion que viaje a traves del yacimiento y Ilegue al contacto de este con el acuifero, conocido
como limite interior del acuifero requerira demasiado tiempo para lIegar allimite exterior del mismo;
de todos modos cuando ya la perturbacion de presion IIegue al limite exterior del acuifero este ya
se consldera finito pero el tiempo para que esto ocurra, dependiendo del tamario del acuifero puede
ser mayor incluso que la vida del yacimiento; un acuifero con una relacion radio exterior del
acuiferol radio exterior del yacimiento mayor de 10 generalmente se puede considerar infinito Un
acuifero se considera finito cuando su tamario no es mucho mayor que el tamario del yacimiento y
por tanto una perturbacion de presion que ha viajado a traves del yacimiento y IIegue al contacto de
este con el acuifero , conocido como limite interior de este, requerira poco tiempo para IIegar al
limite exterior del acuifero; de todas maneras aun cuando el acuifero es pequeno mientras la
perturbacion de presion no IIegue a su limite exterior se considera infinito; acuiferos de relaciones
radio exterior del acuiferol radio exterior del yacimiento menores de 10 se consideran normal mente
finitos ,

Cuando el acuifero empieza actuar es necesano conocer su comportamiento porque afecta el del
yacimiento, especial mente en 10 que tiene que ver con la presion del mismo,

EI comportamiento del acuifero se puede describir con expresiones que se obtienen dependiendo
de la geometria y tamario del mismo y de las condiciones iniciales y de frontera que se establezcan
en cada caso, Las soluciones mas comunes son las de presion y tasa terminal constante tanto para
acuiferos lineales como radiales y en cada caso de tamano infinito 0 finito , Sin embargo las
soluciones mas usadas son las de presion terminal constante para acuiferos infinitos y para los
finitos cerrados , 0 sea los que no tienen comun icacion en ' su limite exterior, Ademas se puede
hablar de una expresion conocida como solucion trivial, la cual es aplicable cuando se tienen
acuiferos finitos y de tamano pequerio y con buena permeabilidad de tal forma que se pueda
pensar que una perturbacion de presion que IIega al acuifero se pueda considerar que se propaga
inmediatamente a todo el sistema ,

1.7.1 -.Descripcion del Comportamiento de Acuiferos Radiales.

64
Es necesario tener una expresion que me relacione la cantidad de agua , 0 la tasa de agua, que
pasa del acuifero al yacimiento con la caida de p,esion y con el tiempo para poder, por eJemplo,
calcular la cantidad de agua que ha entrado al yacimiento y aplicar a este la ecuacion de Balance
de Materiales.

Para todas los casos de soluciones de presion terminal con stante, que ya se dijo es la solucion que
mas se usa en el caso de acuiferos, la expresion que describe el comportamlento del acuffero se
expresa por

W" = ewe * Y(t , p) (1 ,121)

donde ewe se conoce como la constante de intrusion del acuifero y su valor depende de las
propiedades ffsicas y la geometria del acuifero, y Y(t,P) es una ex presion que es funcion del tiempo
y la presion y su forma depende del tamano del acuifero, su geometria y el modelo y regimen de
flujo que se presenta a traves de el.

Existen diferentes modelos para describir el comportamiento de un acuffero depend iendo de su

tamano y el regimen de flujo que se presenta para que el agua pase hacia el yacimiento, Se
analiza ran a continuacion algunos.

1.7.1.1-. Modelo para acuiferos pequenos , Conocida tambien como modelo de acuiferos"Pot" se
basa en la solucion trivial dada por la ecuacion (1 .22):

(1.22)

es una solucion que se puede aplicar en las condiciones mencionadas antes y facilita el analisis del
comportamiento del yacimiento.

Para acomodar la ecuacion (1 .22) a la ecuacion (1 .121) se hace

- (r 2 ­ r"2)* h * qJ * f (e)
ewe -Jr "(I
(1122)
y

Y(t , p)=f')p(t) = (~ -p(t)) (1 .123)

1.7.1.2-. Modelo para Flujo Continuo(8)

Supone que la tasa de intrusion de agua se puede describir como

dW .
e =qw =C g * (Pi - P) =~ (1.124)

don de P es la presion en el contacto acuifero yacimiento y Pi es la presion inicial del sistema

acu ifero - Yacimiento.

como se ve este metodo no requiere conocer caracterfsticas del yacimiento pues C g , que se conoce
como la constante de Intrusion de Schilthuis , no depende de elias y Pi - P son las presiones del
acuifero (inicial) y del CAP . Esta solucion es aplicable en el caso de acuiferos abiertos, presion
constante en ellimite exterior del acuifero.

65
Como Cs es una constante se procede a obtener W e de la siguiente manera :

Supongamos que la historia de comportamiento de la presion del yacimiento es la que se muestra

en la fig ura 11 De la ecuaci6n (1 .124) se tiene:

We = Cs * i(f~ - p) * dl (1 .125)

y la integral se obtiene de la historia de presion asi

P,

l' P 2 fo ( Pi - P) * dt

Pi I
i
k,

Pr.

tl ~ ti t"

-
Figura 11 -.Comportamiento de la Presion en el Contacto Yacimiento - Acuifero

Si la historia de presion es como se muestra en la figura 11 la integral es el area sobre la curva , la

cual para cada tiempo se puede calcular como las areas de los diferentes trapecios desde t = 0
hasta t; 0 sea

f~ (Pi p) * dt =Al + A2 + A3 . An

= ~
2 (p' _ p.1 1112 * [(p _ p ) + (p _ p )] + 111J * [(p _ p ) + (p _ p )]
) +2 ,1, 2 ,2, J
2

+ 11~" [(Pi -Pn-I )+(Pi - Pn)]

l
r (p' _ p) * dl = 11/1
2
* (p' _ p.1
) + 11/ 2 * (2P
2 '
_ P. _ P ) + I1l 1 (2P _ P _ P )
12 , 21
2
+ l1~n * (2Pj - P II
_
1
- Pn)

66
 (1.126)

teniendo en cuenta que Po = P" Si adem as los ~tj se consideran constantes, entonces

( 1.126a)

Llevando la ecuacion (1 .126a) ala ecuacion (1 .125) se tiene

w" = C.\ * t11 ~ (2 P -

2 L.... I
P
,- I
- P )
.I
(1 .127)
)=1

Para lIevar la ecuacion (1 .127) a la forma de la ecuacion (1121) se hace

(1 .128)

e H', "
=c "
(1 .129)

1.7.1.3-. Modelo General de Comportamiento de un Acuifero(8).

Se basa en la solucion de la ecuacion de difusividad para el acuifero baja unas determinadas

condiciones iniciales y de limite aplicables al caso cie acuiferos .

Como se dijo, para el acuifero al igual que para el yacimiento se pueden plantear dos tipos de
expresiones generales dependiendo de si se considera constante la presion en el contacto agua
petroleo , en este caso se puede tener una expresion que me permite calcular la intrusion de agua
We que ha entrado al yacimiento y la expresion se conoce como de presion terminal constante ; 0 si
es constante la tasa de intrusion de agua del aculfero al yacimiento y en este caso la expresion me
permite calcular la caida de presion que se presentara en un intervalo de tiempo dado y la
expresion se conoce como solucion de tasa terminal con stante. La expresion mas usada es la de
presion terminal constante.

• Acuiferos Radiales.

Su forma general se obtiene de acuerdo con Van Everdingen and Hurst (VEH) . Ellos definen una
tasa de intrusion adimensional dada por

(1 .130)

67
Tomando integrales a ambos lad os de la ecuacion (1 .130) sobre el tiempo, es posible obtener una
relacion para W e Y W D donde W e es la intrusion acumulada y W D es la funcion de intrusion
adimensional acumulada .

r CD) * ~
IJ

J) q dl
* dl I> = 21fkhM
j.1
j)
* Jr) qdl (1 .131)

Para el caso de acuiferos tD esta dado por

_ kl it ¢j.1Cro2
I f) - 2 ~ - = -'--'--~ (1 .132 )
¢j.1Cro it I> k

donde ro es el radio del contacto agua - petroleo. Llevando la expresion anterior a la ecuacion
(1131) se tiene

21fkh6.P
j.1
__ * fq dl = "'"Cr02
'f'r
K
*
f
' q I> dl I>

y lIamando W. = fq dl Y WeI> = fq I> dl n se tiene finalmente

W e = 27I h ~C r~ • t-,p · W eD

0
La ecuacion anterior supone un acuifero radial y que actua sobre los 360 del contacto agua ­
0
petroleo; cuando no actua sobre los 360 entonces se incluye un factor f dado por

e e radianes
Ice) = 360 0
21f

y la ecuacion se convierte en

27Ih~C r; . (~)
We=
360 t-, P ·WeD

1
Y si se usan unidades practicas: W e en barriles , h y r en pies , C en Lpc.- y Pen Lpc. se tiene

W e= 1.119 h~C r'; · f( e)· t-, p . W eD(t D) (1 .133)

Para presentar la ecuacion (1 .133) en la forma de la ecuacion (1 .121) se hace

ewc=1 . 119h~C r; (3:0) (1 .134 )

Y(/ , p) = 6.P * W.n(I 1» (1.135)

La ecuacion (1.133) no es general porque supone que la presion en el contacto acuifero-yacimiento
es constante 10 cual no es cierto y por tanto para tener la forma general de Y(t, P) se debe recurrir al
principio de superposicion que se vera mas adelante .

La expresion para WeD depende del tipo de acuifero . En general para el caso de soluciones de
presion terminal constante se habla de tres posibilidades :

• Acuifero infinito
• Acuifero finito cerrado, donde no hay comunicacion hidraulica del acuifero con los alrededores
del limite exterior
• Acuifero finito abierto, donde hay comunicacion del acuifero con los alrededores del limite
exterior y la presion del acuifero en el limite exterior se mantiene constante y por tanto en todo
el acuifero porque el agua que sale hacia el yacimiento es reemplazada par la fuente externa
del acuifero; en este caso en el acuifero se presentara una situacion de flujo continuo 0 estable,
cuyo de modelo de comportamiento ya se presento antes.

Cuando el acuifero es infinito la expresion para WeD presentada por VEH es:

(1.136)

donde J o Y Yo son las funciones Bessel de primera y segunda clase y orden 0 respectivamente .

Klins , Bouchard y Cable (7) muestran que la ~cuacion (1.136) toma las siguientes formas,
dependiendo del valor de to

= 1.1284* , I f) (1 .137)

Si 0.01 < to < 200

W e/) (I /) =
1.12955 ,.jl f) + 1.16044 * 1IJ + 0.2642821 i\12 + 0.0113181
- c - - .... - .
;t 912 . (1 .138 )
) + 0.5900113 ) 1I) + 0.045897 * I/)
12
Si 200<to <2.0*10

Wef) (I /) = 10 t (4.39890 + 0.43693 * ]n 1 f) - 4.16078 * (In 1IJ ) O 09) (1139)

Cuando el acuifero es finito cerrado Weo(t o) depende de reO Y de to y la expresion general dada por
VEH es

(1 .140)

donde Jo Y J 1 son las funciones Bessel de primera clase y orden 0 y 1 respectivamente y los an son
las raices de la ecuacion

69
 J 1 ( all r ell ) Y o( a ll ) - Y 1 ( a n r eI! ) J o( an) =0 (1 .141 )

Donde Y1 Y Yo son las funciones Bessel de segunda clase y arden 0 y 1.

Las funciones Bessel de primera clase estan definidas par las siguientes expresiones dependiendo
del orden y del valor de x:

• Funciones Bessel de Primera clase y Orden 0

O:s;x<3.0

Jo(x)=h o +hl(X/3)2 +h2eX/3) 4 +hl,(xl3)6 +h4(X/3)8 +bs (x/3)' O+hf,(x/3)1 2 (1 .142)

3.0:S;x<oo

.1 o ex) =(xt" 2 F,, (x)[cos(Boex))] (1 .143)

Fo(x)=h o +hle3 / x)+b2e3 / x) 2 +hl, e3 / x) ' +h4(3/X)4 +hs(3 /x )s +h6e3/x)6 (1 .144)
Boex)=x+ho +h,(3 / x)+h2(3 / x) 2 +b J e3 / x) 1, +h4e3/x)4 +hse3 / x) 5 +hc,(3/x)C, (1 .145)

• Funciones Bessel .de Primera Clase y Orden 2:

O:s;x<3.0
(, 2"

(xt'.1 1 (x)= Ib, (x / 3) ( 1.146)

,=0

3.0:s;x<oo

J1ex) =(Xt' 12 F, (x)[cos(B, (x))] (1.147)

F,(x)=h o +b,(3 / x)+b2e3 / x) 2 +hl, e3 / x)' +h4(3 / x)4 +hs(3 / x) 5 +h6(3/x)6 (1.148)

B,ex)=x+ho +hle3 / x)+h2e3 / x) 2 +bl, (3/x)' +b 4(3/x)4 +hse3/x) s +h6e3/x)6 (1.149)

Los b; de las ecuaciones (1 .142), (1.44)-(1 .145), (1 .146) Y (1.148)-(1 .149) se dan en la tabla 3.

Tabla 3-.Constantes para las Ecuaciones(1.142), (1.144-45), (1.146)y (1.148-49) (7)

Jo J1
O:s;x<3 3.0:s;x<oo O:s;x<3 3.0:S; x<oo
Fo 90 F1 91
bo 1 0.79788456 -0.78539816 0.5 0.79788456 -2 .3561944
b1 -2.2499997 -0.00000077 -0 .04166397 -0.5624998 0.00000156 0.12499612
b2 1.2656208 -000552740 0.00003954 0.21093573 0.01659667 0.00005650
b3 -.3163866 -0 .00009512 0.00262573 -0.0395429 0.00017065 -0 .0063788
b4 0.0444479 0.00137237 -0 .00054125 0.00443319 -0.0024951 0.00074348
bs 0.0039444 -0.00072805 -0 .00029333 -0 .0003176 0.00113653 0.00079824
b6 _ 0.0002100 0.00014476 0.00013558 0.00001109 -00002003 -00002917

70