Conificacion PDF
Conificacion PDF
Conificacion PDF
1 INTRODUCCIN
Uno de los problemas que se presenta en la explotacin de un reservorio de hidrocarburo es la
produccin de agua asociada al mismo. Estudiaremos en este trabajo el caso de un reservorio
de petrleo. Una prctica comn es realizar los punzados en los pozos productores a una
distancia segura sobre el contacto agua-petrleo. Si la distancia no es muy grande se genera
un gradiente vertical del potencial de fluido, con un mximo en el eje del pozo, que produce
una deformacin en forma de cono de la superficie de contacto. Este efecto es conocido como
conificacin de agua. Un efecto similar se produce con el contacto gas-petrleo para el caso
de reservorios con casquete de gas.
Entre los efectos adversos de la presencia de conificacin de agua se pueden destacar los
siguientes:
Cuando en un pozo se requiere producir petrleo por encima del caudal crtico, es necesario
realizar una evaluacin de (1) el tiempo en el cual el cono alcanza el extremo inferior del
intervalo punzado y por lo tanto comienza la produccin de agua (tiempo de breakthrough) y
(2) la evolucin del corte de agua a partir del tiempo de breakthrough.
Existen dos trabajos publicados en relacin con la determinacin del tiempo de breakthrough,
uno de Sobocinski y Cornelius[1] y otro de Bournazel y Jeanson[2]. En ambos trabajos utilizan
las siguientes formas adimensionales para la altura del cono y el tiempo:
H D,v =
2 g h 02 ( w - o ) k r
h D,v
q o B of o
tD =
h D,v =
a
g ( w - o )k v 1+ M wo
t
2 o f ho
donde
hv
ho
(1)
(2)
ho
hv
rw,ro
kr,kv
qo
Bof
mo
mw
f
(t )
SC
D BT
(H )
D,v BT
(3)
(t D )BJBT =
(H )
3 - 0.7 (H )
D,v BT
(4)
D,v BT
Los tiempos estimados con la correlacin de Bournazel y Jeanson son menores que los
estimados por Sobocinsky y Cornelius; de 40 a 50% menores para valores de HD,v entre 0.1 y
3[5].
En cuanto al estudio del comportamiento del corte de agua en funcin del tiempo, Kuo y
DesBrisay[3] desarrollaron una correlacin a partir de una simulacin numrica. Partiendo de
un caso de referencia con valores prefijados de las variables, realizaron algunos clculos de
sensibilidad variando de a una algunas de las variables involucradas en el problema,
obteniendo la siguiente expresin para el corte de agua:
WC = 0
para
t
t BT
< 0.56
WC
t
+ 0.24
= 0.94 log
(WC)lim t
t BT
WC
= 1.0
(WC)lim t
para
t
t BT
para 0.56
t
t BT
5.7
(5)
> 5.7
r
petrleo
lp
ho
agua
P,T cte
P>Pbbl
WOC inicial
hw
z
Figura 1. Modelo utilizado para el estudio de conificacin de agua.
S - Siw
= w
1 - Siw
( ) (2 - S )
k ro = S*OF
*
OF
con S
*
OF
(6)
12
donde:
Siw = saturacin de agua irreducible
Sor = saturacin de petrleo residual
Podemos hacer la hiptesis adicional de que Swi=Sor ya que ambos valores no suelen diferir
demasiado (salvo cuando la mojabilidad de la roca es un factor dominante), por lo que la
nica variable independiente resultante es Siw.
Otras variables an necesarias para caracterizar el sistema fueron fijadas considerando valores
promedios en reservorios de petrleo[6] por lo que se mencion anteriormente:
Presin (150 bar). Esta presin asegura que en todo el modelo se est por encima de la
presin de burbuja .
Temperatura (60)
Espesor de la capa de petrleo (30m)
Espesor de la capa de agua (30m). Como era de esperar, no se encontr ninguna
variacin en los resultados al variar este espesor.
Factor de volumen del petrleo (Bo)
Viscosidad del agua: 0.5 cP
Viscosidad del petrleo: 7.2 cP
Densidad del agua: 1 g/cm3
6
Estos valores describen el tipo de fluidos que en promedio suelen encontrarse en un reservorio
de petrleo negroa pero no se hicieron estudios acerca de la sensibilidad del modelo a las
variaciones de los mismos. Es recomendable realizar este estudio en el futuro para estar
seguro de que no hay que incluir alguna de estas variables dentro del conjunto que s se vara.
Con estas caractersticas se construy un modelo de simulacin tipo black oil basado en el
software ECLIPSE 100. Este programa es un simulador de reservorios de petrleo (black
oil) con una amplia capacidad de configuracin, que permite simular la mayora de las
operaciones que se realizan en un yacimiento de petrleo. El movimiento de los fluidos se
determina basndose en la ecuacin de continuidad y la ley de Darcy.
La ecuacin de continuidad para cada fase k se escribe as[8]:
(r k fS k ) + (r k vr k ) = -q k
t
(7)
r
donde la velocidad superficial de la fase k ( v k ) se define como:
r
r
v k = fS k u k
(8)
r
siendo u k la velocidad poral de la fase k, rk la densidad, qk la fuente de masa por unidad de
volumen, Sk la fraccin de volumen que ocupa cada fase y f la porosidad de la roca.
La ecuacin de Darcy para flujo multifsico es la siguiente[8]:
kk
r
r
v f = - rf (p f + r f g )
mf
(9)
r
donde f corresponde al tipo de fluido (petrleo, gas o agua), v f es la velocidad del fluido, krf
es la permeabilidad relativa al fluido f (que es tomada como un escalar en el simulador), k es
r
la permeabilidad absoluta, pf es la presin y g es la aceleracin de la gravedad.
El programa resuelve las ecuaciones anteriores utilizando un mtodo implcito de diferencias
finitas[7].
Los tiempos de corrida son de aproximadamente 8 segundos por caso. Se realizaron diversos
estudios de sensibilidad destinados a optimizar la discretizacin espacial y temporal de forma
a
Estas hiptesis pueden no ser correctas para petrleos voltiles, gas-condensado y petrleos pesados
tal de obtener una buena convergencia en tiempos de CPU razonablemente bajos. Al aumentar
la discretizacin se obtiene una mejor convergencia pero con un aumento de los tiempos de
CPU.
Malla Espacial
En direccin r, se utiliz una malla ms fina en los alrededores del pozo productor de
petrleo central.
En direccin z se represent en forma ms detallada la zona cercana al contacto aguapetrleo, especialmente en la zona de petrleo.
Discretizacin temporal.
Se compararon diferentes discretizaciones en tiempo, de forma tal de obtener una
buena convergencia en tiempos de CPU razonables. Se utilizaron pasos temporales de
5 das en los primeros dos meses y de un mes en el lapso restante.
En la figura 2 se presenta el modelo utilizado. Los pozos reinyectores de petrleo y agua se
ubican en la ltima celda de la grilla en direccin radial.
La variable que utilizamos para caracterizar la produccin de agua en funcin del tiempo es el
corte de agua (WWCT), es decir, el porcentaje de agua producida respecto a la produccin
total de fluidos.
Las curvas de corte de agua en funcin del tiempo tienen la forma general que se muestra en
la figura 3.
Para realizar las simulaciones, se utiliz el programa Eclipse, instalado sobre un sistema
operativo Linux.
Para la ejecucin en forma automtica de la serie de 15625 casos necesarios para la
9
realizacin del presente estudio, se escribi un script de unix con la siguiente funcionalidad:
1. Se construye un archivo base de entrada al simulador, que contiene claves en las
posiciones donde deben ir los valores de los parmetros a modificar.
2. Para cada uno de los casos, el script sustituye los valores correspondientes de las
variables en el archivo base de entrada, generando un nuevo archivo de entrada
particular para el caso correspondiente.
3. Se ejecuta el simulador con el archivo de entrada generado.
4. Se procesa la salida del programa, extrayendo de la misma todos los valores de
inters y almacenndolos en un archivo de texto.
5. Se eliminan todos los archivos generados por el simulador que ya no son necesarios.
Este paso es indispensable para evitar la utilizacin de un enorme espacio de
almacenamiento.
6. Se repite el proceso desde el paso 2 para todos los casos.
El archivo de texto que contiene los resultados de todas las simulaciones es posteriormente
procesado utilizando el programa Matlab.
3 RESULTADOS DEL ESTUDIO
Se obtuvo el corte de agua en funcin del tiempo para cada una de las 15625 combinaciones
de valores de parmetros definidas. Se descartaron los casos en los que las condiciones del
reservorio no permiten producir petrleo durante 5 aos al caudal constante predeterminado y,
con los casos resultantes se gener una base de datos que contiene el tiempo de breakthrough
tBT, el corte de agua mximo wwct5 (5 aos) y el tiempo en alcanzar el 70% del corte de agua
mximo t70 en funcin de los parmetros operacionales y del reservorio.
En la figura 4 se muestra el corte de agua (normalizado al corte de agua mximo en cada caso,
wwct5) en funcin del tiempo (normalizado al tiempo de breakthrough, tBT), para cuatro casos
seleccionados.
Si definimos tmax como el tiempo en que se llega a un porcentaje fijo de la produccin
mxima de agua (por ejemplo 0.99*wwct5), se puede observar que en los casos presentados
tmax vara entre aproximadamente 10*tBT y 100*tBT.
Se puede observar claramente que, haciendo variar las diferentes variables simultneamente,
no es posible caracterizar el corte de agua en funcin del tiempo con una sola curva. En el
trabajo de Kuo y DesBrissay[3] logran una nica correlacin debido a que realizan el anlisis
partiendo de un caso de referencia y haciendo variar los parmetros de a uno por vez. Por otra
parte, al adimensionalizar el tiempo se lo divide por el tiempo de breakthrough calculado por
10
otra aproximacin, la de Bournazel y Jeanson[2], por lo que la curva resultante tiene un corte
de agua mayor que cero para un tiempo igual al de breakthrough.
Una metodologa para obtener una descripcin completa del comportamiento del corte de
agua en funcin del tiempo consistira en obtener correlaciones, en funcin de los parmetros
operacionales y del reservorio, para (1) el tiempo de breakthrough, (2) el tiempo en alcanzar
el corte de agua mximo y (3) el corte de agua mximo.
1
Figura 4. Corte de agua en funcin del tiempo para cuatros casos seleccionados.
Para la altura adimensional desde el contacto agua-petrleo hasta el extremo inferior del
intervalo punzado utilizamos la frmula (1), mientras que para el tiempo de breakthrough
adimensional utilizamos la siguiente frmula:
g ( w - o ) k v 1 + M awo
tD =
2o f ho
) 1 + S
11
b - S t BT
iw
iw
(10)
donde
a = 0.6
b = 0.64
g = 1.43
Notar que la definicin de la altura adimensional es la misma que fue utilizada en las
estimaciones de tiempo de breakthrough de Sobocinsky-Cornelius y Bournazel-Jeanson.
g
1 + Siw
y se elige utilizar a = 0.6
En cuanto al tiempo adimensional, se agrega el trmino
b
S
iw
12
Figura 5. Tiempo de breakthrough adimensional vs. altura adimensional desde el contacto agua-petrleo hasta la
base del intervalo punzado: todos los puntos y detalle para alturas adimensionales menores que 1.
13
Con la eleccin de los parmetros de ajuste a = 0.6, b = 0.64 y g = 1.43, resulta la siguiente
relacin entre el tiempo de breakthrough adimensional tD y la altura de punzado adimensional
HD,v:
t D = H D ,v
H D,v < 1
k r 100md
(11)
Como verificacin, y con el objeto de comparar nuestros resultados con los de BournazelJeanson, comprobamos que aquella correlacin da resultados similares a la nuestra para un
valor de Siw apenas superior a 0.1. Esto parece indicar que esa fue la saturacin irreducible en
el experimento que ellos utilizaron para encontrar su correlacin.
La correlacin que proponemos en la Ec.(11), entonces, es ms general que la de BournazelJeanson, ya que incluye la dependencia adicional con la saturacin de agua irreducible, a
1 + Siw
travs del agregado del coeficiente
b - Siw
en la adimensionalizacin temporal.
4 CONCLUSIONES
hidrocarburos Eclipse fue posible gracias a la cesin de una licencia acadmica por parte de
Schlumberger-Geoquest.
6 REFERENCIAS
[1] Sobocinski, D.P and Cornelius, A.J.: " A Correlation for Predicting Water Coning Time",
Journal of Petroleum Technology (May 1965) 594-600. Paper SPE 894.
[2] Bournazel, C. and Jeanson, B: "Fast Water-Coning Evaluation Method" (1971). Paper
SPE 3628.
[3] Kuo, M.C.T and DesBrisay, C.L.; "A Simplified Method for Water Coning Predictions".
Paper SPE 12067 presented at the 1983 SPE Annual Technical Conference and Exhibition,
San Francisco.
[4] "Principles of Petroleum Reservoir Engineering". Chierichi, G.L. (1994) Springer-Verlag
[5] "Applied Reservoir Engineering". Smith, C.R., Tracy, G.W. and Farrar, R.L. (1992) Oil &
Gas Consultants International Publications.
[6] "Petroleum Engineering Handbook". H.B. Bradley, SPE (1987).
[7] "Eclipse Technical Description", Manuals 2002A, 2002.
[8] "Simulacin de Reservorios Aplicada". - Notas del curso dictado en UADE del 13 al 17 de
mayo de 2002. P. Carrica, F. Nissero.
15