EJEMPLO DE CALCULO DE INYECCION DE AGUA

El yacimiento que se toma como ejemplo tiene las propiedades indicadas en la Tabla 1.
Se harán los cálculos a los datos del yacimiento que se toma como ejemplo, con el
objeto de ilustrar el uso de los diferentes métodos de predicción de comportamiento de
la inyección de agua. Este yacimiento a sido sujeto de depleción de presión por empuje
de gas, y luego por inyección de gas por un período breve.

Tabla 1 Propiedades Del Yacimiento

PROPIEDADES DEL YACIMIENTO DATOS UNIDADES

Espaciamiento entre pozos 20 Acres

Espesor 50 Pies

Permeabilidad Promedio 10 md.

Porosidad 20 %

Saturación de agua congénita 10 % VP

Saturación Actual, promedio de gas 15 % VP

Viscosidad del petróleo, a la presión 1.0 Cp

actual de yacimiento

Viscosidad del agua 0,5 Cp

Presión del yacimiento 1,000 Lb./pulg2.

Distribución de la permeabilidad Ilustrada en la

Figura E.3

Características de permeabilidad Ilustradas en las

relativa Agua - Petróleo Figuras 2.17ª y
2.17b

Recuperación actual de petróleo 10.4 % del petróleo

inicial I.S.

Factor de volumen del petróleo a la 1.29

presión original de saturación

Presión Actual 1.20

Arreglo de Inyección Cinco pozos,

usando los
 existentes.

Área del Arreglo * 40 Acres

Radio del pozo 1 Pies

*El arreglo completo de cinco pozos contiene dos pozos en total un

inyector y un productor.

E.1 Cálculo de la curva de flujo fraccional y del comportamiento del

Desplazamiento.
E.2 Calculo de la Relación de Movilidad
E.3 Cálculo de la Recuperación Total por Inyección de Agua
E.4 Comportamiento Combinado WOR-Recuperación
E.5 Combinación de Ritmos de Producción e Inyección, WOR, y Recuperación
contra tiempo.
E.6 Comportamiento Individual de Pozos.

E.1 Cálculo de la curva de flujo fraccional y del comportamiento del

Desplazamiento.
Para determinar la curva de flujo para agua desplazando aceite, se usa La siguiente
ecuacion:
 1
𝑓𝑤 = 𝜇𝑤 𝑘𝑟𝑜 E. 1
1+
𝑘𝑟𝑤 𝜇0

El yacimiento que se toma como ejemplo, considere que se aplica las caracteristicas de
permeabilidad realtiva agua-aceite se muestra en la figura 1.
Figura 1 Caracteristicas de permeabilidad relativa agua-aceite, con unaroca

marcada mojabilidad por aceite y agua.

Agua aceite

Considere tambien que la viscosidad del aceite y del agua a las condiciones del
yacimiento son respectivamente 1,0 y 0,5 cp y que la satuacion inicial de agua es de
10% VP.

En la tabla se puede ver la relacion saturacion-flujo fraccional.

 Tabla 2 Curva de flujo Fraccional, problema ejemplo

La curva de flujo fraccional calculada, que aparece en la figura 2 es identica a las de las

Figura 2 Curva de Flujo Fraccional, problema ejemplo

Figuras 3.2. y 3.8.

Se traza una tangente a la curva del flujo fraccional desde la saturacion inicial de agua
10% VP. La saturacion de agua en la tangencia es de 0,469 , o 46.9% VP.
𝑆𝑊𝑆𝑍 = Saturacion de agua en el extremo corriente arriba de la zona estabilizada. Si
prolongamos la tangente hasta que intersecte la linea horizontal que corresponde a

𝑓𝑤 = 1.0, la interseccion acurrira a una saturacion de agua de 0,563, o 56.3% VP. Esta
es la saturacion promedio de agua detrás del frente de invasion a y antes de la surgencia
de agua, 𝑆𝑤𝑏𝑡 .

A partir de la curva del flujo fraccional, para saturaciones, de agua iguales y mayores a
la saturacion de agua en el extremo corriente arriba de la zona estabilizada,𝑆𝑤𝑠𝑧 (46,9%
VP), se determina la pendiente de la curva de flujo fraccional.

El valor de 𝑄𝑖 , los volumenes de poros de agua inyectada acumulativa requerida para

alcanzar la saturacion 𝑆𝑤2 en el extremo productor del sistema de roca, a y despues de la
surgencia, se determina mediante la ecuacion:

1
𝑄1 = 𝑑𝑓𝑤 E. 2
(𝑑𝑠 )𝑆𝑊2
𝑤

Usando la ecuacion de welge

𝑆𝑤 = 𝑆𝑤2 + 𝑄1 𝑓𝑛2 E. 3

E.2 Calculo de la Relación de Movilidad

La movilidad es la relacion entre la permeabiliad efectiva y la viscosidad de fase.

El valor de la relacion de movilidad agua-aceite se calcula de la siguiente ecuacion:

𝑘1𝑤 𝜇𝑜
𝑀= E. 4
𝜇𝑤 𝑘1𝑜

La permeabilidad relativa al agua a la saturacion d agua promedio (56.3%PV) se lee en

la figura 1 y resulta a 0,4. La permeabilidad relativa al aceite adelante del frente de
barrido es igual a 1,0.

Reemplazando datos:
0.4 1.0
𝑀 = 0.5 1.0 =M=0.80

E.3 Cálculo de la Recuperación Total por Inyección de Agua

En el yacimiento que se considera como ejemplo se tiene la siguiente orden de valores
de la permeabilidad desde el maximo hasta el minimo, ver la figura 2.

Figura 3 Distribucion de Permeabilidad, problema ejemplo

Esta distribucion de permeabilidades tiene unavariacion de permeabilidad de 0.5.

ademas el yacimiento ha sido parcialmente agotado por el empuje de gas disuelto y la
recuperacion a la fecha ha sido del 10,4% del aceite original in situ.

Para calcular la recuperacion total por inyeccion de agua hasta unafraccion de agua
producida del 98%, primero ir a la tabla 3 y determinamos que esta fraccion de agua
corresponde a una saturacion d agua promedio del 66.6% VP.
E.4 Comportamiento Combinado WOR-Recuperación

El metodo recomendado para calcular el comportamiento combinado WOR-

Recuperacion, es el propuesto por Dystra and Parsons.

Utilizando las siguientes figuras podemos calcular la recuperacion fraccional de aceite

Figura 4 Variacion de la permeabilidad graficada contra la relacion de movilidad,

mostrado lineas de constante 𝑬𝑹 (𝟏 − 𝑺𝑾 ), para un WOR de produccion de 1.

para WOR de produccion de 1,5,25 y 100.

Figura 5 Variacion de permeabilidad graficada contra la relacion de movilidad,

 mostrado lineas constante 𝑬𝑹 (𝟏 − 𝟎, 𝟕𝟐 𝑺𝑾 ), para un WOR de produccion de 5.

Figura 6 Variacion de permebilidad graficada contra la relacion de movilidad,

mostrado lineas constante𝑬𝑹 (𝟏 − 𝟎, 𝟓𝟐 𝑺𝑾 ) , para un WOR de produccion de 25.

Figura 7 Variacion de la permeabilidad graficada contra la relacion de movilidad,

mostrado lineas de constante 𝑬𝑹 (𝟏 − 𝟎, 𝟒𝟎 𝑺𝑾 ) para un WOR de produccion de 100.

Se entra en a figura 5 con la relacion de mojabilidad de 0,8 y una variacion de

permeabilidad de 0,5. Para un WOR de produccion igual a 1, tenemos:

𝐄𝐑 (𝟏 − 𝐒𝐖 ) = 𝟎, 𝟐𝟎𝟎

𝟎. 𝟐𝟎 𝟎, 𝟐𝟎
𝐄𝐑 = = = 𝟎, 𝟐𝟐𝟐
𝟏 − 𝑺𝒘 𝟎, 𝟗𝟎

Para un WOR de produccion de 5,(figura 6) :

𝐄𝐑 (𝟏 − 𝟎, 𝟕𝟐𝐒𝐖 ) = 𝟎, 𝟐𝟗

𝟎. 𝟐𝟗 𝟎, 𝟐𝟗
𝐄𝐑 = = = 𝟎, 𝟑𝟏𝟐
𝟏 − 𝟎, 𝟕𝟐𝑺𝒘 𝟎, 𝟗𝟐𝟖

E.5 Combinación de Ritmos de Producción e Inyección, WOR, y Recuperación

contra tiempo.

Este metodo es recomendado para obtener el comportamiento de la inyeccion de agua

con el tiempo para un arreglo de cinco pozos es (1) la solucion analitica de Craig,
Geffen y Morse para relacionar la recuperacion del aceite y la WOR de produccion con
el agua inyectada acumulativa, combinado esto con (2) la correlacion de Caudle y
Witte para calcular los gastos de inyeccion de agua en um arreglo de cinco pozos.

Como medio para incluir los efectos de la heterogeneidad del yacimiento, puede
calcularse el comportamiento de cinco pozos estraficados. Para determinar el minimo
numero de capas recomendado, ver la tabla 7.1. esta tabla indica que para una relacion
de movilidad de 0,8 y una variacion de permeabilidad de 0.5, 10 capas de igual espesor
daran el mismo comportamiento calculado por inyeccion de agua, que un numero muy
superior de capas, para una WOR de 2,5.

La tabla 4 presenta las permeabilidades y espesores para 10 capas, con una variacion de
permeabilidad de 0,5. Estas permeabilidades estan tomadas de la figura que representa
las permeabilidades promedio para cada incremento del 10% de la muestra acumulativa.
Estos valores son las permeabilidades selecionadas para 5, 10, 15…..95% de la muestra
acumulativa.
 Tabla 3 Propiedades de 10 capas de igual espesor con una variacion de
permeabilidad de 0,5 para el yacimiento ejemplo

CAPA PERMEABILIDAD (mD) ESPESOR(pies)

1 31.5 5.0
2 20.5 5.0
3 16.0 5.0
4 13.1 5.0
5 10.9 5.0
6 8.2 5.0
7 7.7 5.0
8 6.3 5.0
9 4.9 5.0
10 3.2 5.0

El yacimiento estratificado de este ejemplo esta compuesto por capas indeticas en todas
sus propiedades, con excepcion de la permeabilidad. Por lo tanto, para calcular el
comportamiento de una capa (la llamada capa base) y obtener despues, el
comportamiento para cada una de las capas restantes ajustandolo de acuerdo con el
contraste de permeabilidad.

El metodo de este ajuste se expondra mas adelante en este apendice. Para yacimientos
estratificados con capas de diferentes caracteristicas de permeabilidad relativa agua-
aceite, es necesario calcular individualmente el comportamiento de cada una de las
capas. En ambos casos, el comportamiento combinado es la suma de los
comportamientos individuales de cada una de las capas.

Estos calculos suponen una eficiencia vertical de barrido de 100% en cada capa es decir,
los fluidos no segrregan por gravedad. El comportamiento de la inyeccion de agua
puede dividirse en cuatro etapas.
  Etapa uno es un periodo de flujo radial a partir de los pozos inyectores, desde el
inicio de la inyeccion hasta que los bancos de aceite formados alrededor de los
pozos inyectores, se juntan. El encuentro de bancos adyacentes de aceite se
llama interferencia.
 Etapa dos es el periodo desde la interferencia hasta el llenado del espacio de gas
pre-existente. El llenado es el inicio de la repuesta en la produccion de aceite.
 Etapa tres es el periodo desde el llenado hasta la sugerencia de agua en los pozos
productores. La surgencia senala el comienzo de la produccion de agua.
 Etapa cuatro es el periodo desde la surgencia hasta la inundacion total.

En este apendice discutiremos el calculo del comportamiento de la inyeccion de agua,

durante cada una de estas etapas.

Para comenzar, se requieren algunos calculos iniciales, estos aparecen en la Tablae5

Estos calculos incluyen la determinacion de el volumen de poros; el aceite in situ en el
momento de la inyeccion de agua y la eficiencia areal de barrido a la sugerencia de
agua. Los otros calculos que aparecen en la tablae5 se discutiran mas adelante.

Comportamiento Antes De La Inteferencia

Durante esta etapa de la inyeccion, prevalece el flujo radial. El agua inyectada a la

inteferencia, 𝑊𝑖𝑖 , es igual a la saturacion de gas libre en la porcion cilindrica del
yacimiento, con un radio 𝑟𝑒𝑖 es la mitad de la distancia entre pozos inyectores
adyacentes. De la tabla, encontramos que el volumen de agua inyectada a la
interferencia es de 36,572 Bls.

el gasto de inyeccion de agua antes de la interferencia es:

7,07∗10−3 ℎ𝑘∆𝑝
𝑖𝑤 = 𝜇 𝑟 𝜇 𝑟 E. 5
( 𝑤 𝑙𝑛 + 𝑜 𝑙𝑛 𝑒 )
𝑘𝑟𝑤 𝑟𝑤 𝑘𝑟𝑜 𝑟

Donde:

𝑟𝑒 = radio exterior del banco de aceite

r = radio esterior del frente de barrido de agua.

La tabla muestra el calculo de la variacion del gasto de inyeccion desde el inicio de la
inyeccion deagua hasta la interferencia, la tabla E.6 utiliza algunas de las ecuaciones
predentadas en la tabla e5. Este calculo indica que la inteferencia ocurrira a los 75,3dias
despues del inicio de la inyeccion de agua.

1. Volumen de poros (𝑉𝑝 ) = 7,758 𝐴ℎ∅ = 7.758 ∗ 40 ∗ 5,0 ∗ 0,20 =

310,320 𝐵𝑙𝑠.
𝑉𝑝 ∗𝑆𝑜𝑖
2. Aceite in situ a condiciones de tanque al comenzar la inyeccion de agua= =
𝐵𝑜𝑖
310320∗0,75
= 193950 𝐵𝑙𝑠.
1,20

3. Eficiencia de area barrida a la surgencia de agua, 𝐸𝑎𝑏𝑡 = 0,717 de la figura e6 a

M=0,8
4. Agua inyectada a la surgencia, 𝑊𝑖𝑏𝑡 = 𝑉𝑝 ∗ 𝐸𝑎𝑏𝑡 ∗ (𝑆𝑤𝑏𝑡 − 𝑆𝑤𝑐 )
= 310320 ∗ 0,717 ∗ (0,563 − 0,10) = 103020 𝐵𝑙𝑠.
a. Valor maximo de saturacion inicial de gas, para la cual el metodo de prediccion

es exacto,
b. El valor de C se encuentra en la Figura e.7 para:

𝑀 = 0,8. 𝑆𝑔𝑖 ∗= 1,18(0,75 − 0,437) = 0,369 𝑜 36,9%

c. Puesto que la saturacion inicial de gas, 15% PV, es menos que el valor de 𝑆𝑔𝑖 ∗,
el metodo de prediccion es aplicable.
𝑥ℎ∅ 𝑆𝑔𝑖 (𝑟𝑒𝑖 )2
5. Agua inyectada a la interferencia, 𝑊𝑖𝑖 = , donde 𝑟𝑒𝑖 es la
5,61
mitad de la distancia entre pozos inyectores adyacentes, 0 660 pies,
𝑥 ∗ 5,0 ∗ 0,20 ∗ 0,15 ∗ (660)2
𝑊𝑖𝑖 = = 36,572 𝐵𝑙𝑠.
5,61
6. Agua inyectada al llenado, 𝑊𝑖𝑓 = 𝑉𝑝 ∗ 𝑆𝑔𝑖 = 310320 ∗ 0,15 = 46550 𝐵𝑙𝑠.
7. Eficiencia de area barrida, 𝐸𝐴, para un volumen de agua inyectado, 𝑊𝑖, antes de
la surgencia de agua.
 𝑊𝑖 𝑊𝑖 𝑊𝑖
𝐸𝐴 = = =
𝑉𝑝 (𝑆𝑤𝑏𝑡 − 𝑆𝑤𝑐 ) 310320(0,563 − 0,1) 143678

8. Radio exterior del banco de aceite, 𝑟𝑒 , antes de la interferencia:

5,61 𝑊𝑖 1⁄ 5,61 ∗ 𝑊𝑖 1 1
𝑟𝑒 = ( ) 2 =( ) ⁄2 = (11915𝑊𝑖 ) ⁄2
𝑥ℎ∅ 𝑆𝑔𝑖 𝑥5,0 ∗ 0,20 ∗ 0,15

9. Radio exterior del frente de barrido de agua, r, antes de la interferencia:

𝑆𝑔𝑖 1 0,15 1
𝑟 = 𝑟𝑒 ( ) ⁄2 = 𝑟𝑒 ( ) ⁄2 = 0,5692𝑟𝑒
𝑆𝑤𝑏𝑡 − 𝑆𝑤𝑐 0,563 − 0,10

10. Gasto de inyeccion de agua, 𝑖𝑤 , hasta la interferencia para una capa de 31,5 md
o 0,0315 darcies y ∆𝑝 = 3000𝑝𝑠𝑖:

𝜇𝑤 𝑟 𝜇𝑜 𝑟𝑒 −1
𝑖𝑤 = 7,07 ∗ 10−3 ℎ𝑘∆𝑝( 𝑙𝑛 + 𝑙𝑛 )
𝑘𝑟𝑤 𝑟𝑤 𝑘𝑟𝑜 𝑟

𝜇𝑤 𝑟 𝜇𝑜 𝑟𝑒 −1
𝑖𝑤 = 7,07 ∗ 10−3 ∗ 5,0 ∗ 31,5 ∗ 3000( 𝑙𝑛 + 𝑙𝑛 )
𝑘𝑟𝑤 𝑟𝑤 𝑘𝑟𝑜 𝑟

𝜇𝑤 𝑟 𝜇𝑜 𝑟𝑒 −1
𝑖𝑤 = 3340( 𝑙𝑛 + 𝑙𝑛 )
𝑘𝑟𝑤 𝑟𝑤 𝑘𝑟𝑜 𝑟

11. Gasto basico de inyeccion de agua:

3,541 ∗ 10−3 ℎ 𝑘 𝑘𝑟𝑜 ∆𝑝

𝑖𝑏𝑎𝑠𝑒 =
𝜇𝑜 (𝑙𝑛 𝑑⁄𝑟𝑤 − 0,619)

3,541 ∗ 10−3 ∗ 5,0 ∗ 31,5 ∗ 1,0 ∗ 3000

𝑖𝑏𝑎𝑠𝑒 =
1,0(𝑙𝑛 932⁄1 − 0,619)

𝑖𝑏𝑎𝑠𝑒 = 269,1 𝐵𝑙𝑠.⁄𝑑𝑖𝑎

Tabla 6. Continuacion
La figura e9 presenta el comportamiento calculado WOR-recuperacion para unasola
capa. El aumento brusco de WOR producida para aproximadamente 49% de
recuperacion, ocurre en el punto en el cual la eficiencia de area barrida llega al 100%.

Comportamiento de las capas restantes

Para calcular el comportamiento de otras capas que pueden diferir en espesor, porosidad
y permeabilidad, se sigue la logica siguiente: Si el indice 1 indica la primera capa y el
indice n indica la capa de orden n:

Tabla 4 Valores de para varios valores de eficiencia de area barrida a la surgencia

Figura 8 Eficiencia de area barrida a la surgencia de agua, arreglo de cinco pozos

Figura 9 Relacion entre el Coeficiente C y la relacion de movilidad

Figura 10 Efecto del volumen de fluido inyectado sobre el incremento de area barrida, despues de
la surgencia de agua

Figura 11 Comportamiento calculado WOR - recuperacion,

inyeccion de agua en una sola capa, cinco pozos
 Figura 12 Comportamiento calculado WOR-recuperacion
Figura 13 Comportamiento calculado de la inyeccion
,yacimiento estratificado, cinco pozos