Profa.: Ing. Varinia Marcano D.

INTRODUCCIÓN

Una gran cantidad de yacimientos de petróleo y gas tienen un acuífero

asociado que representa una fuente importante de energía.

Esta energía provee un mecanismo de empuje para la producción de

fluidos cuando los yacimientos son sometidos a producción.

Se cree que el gran número de yacimientos con empuje de agua está

relacionado con el origen marino de muchos yacimientos.

En los casos en que el volumen del acuífero es menos de 10 veces el volumen

del yacimiento, el mecanismo de empuje por agua es considerado pequeño. Si
el tamaño del acuífero es significativamente mayor (> 10x), el mecanismo de
empuje por agua puede ser la principal fuente de energía de yacimiento.
 MODELOS DE ACUÍFEROS

REGÍMENES DE FLUJO
Hay básicamente tres tipos de regímenes de flujo que debe ser reconocidos con el
fin de describir el comportamiento del flujo de fluidos y la distribución de la presión
en función del tiempo. Estos tres regímenes de flujo son los siguientes:

 Flujo en estado estacionario o estable;

 Flujo en estado inestable;
 Flujo en estado pseudoestable

FLUJO EN ESTADO ESTACIONARIO

El régimen de flujo se identifica como un flujo en estado estable, si la presión en cada
lugar en el yacimiento se mantiene constante, es decir, no cambia con el tiempo.
Matemáticamente, esta condición es expresa como:

𝜕𝑝
= 0
𝜕𝑡 i
 MODELOS DE ACUÍFEROS

Esta ecuación establece que la tasa de cambio de la presión, con respecto al

tiempo, en cualquier lugar es igual a cero.

FLUJO EN ESTADO INESTABLE

La tasa de cambio de presión con respecto al tiempo en cualquier posición del
yacimiento no es cero o constante.

𝜕𝑝
= 𝑓(𝑖, 𝑡)
𝜕𝑡

FLUJO EN ESTADO PSEUDOESTABLE

Establece que la tasa de cambio de presión con respecto al tiempo en cada
posición es constante.
𝜕𝑝
= constante
𝜕𝑡 i
MODELOS DE ACUÍFEROS

Fuente: Ribón, H. (2019)

 MODELOS DE ACUÍFEROS

Los modelos matemáticos de influjo de agua comúnmente utilizados en la industria

petrolera son:

(a) Pot
1. Estado Estable (b) Schithuis (1936)
(c) Hurst (1943)

(a) Van Everdingen-Hurst (1949)

2. Estado inestable (b) Carter-Tracy (1960)
(c) Fetkovich (1971)
 MODELOS DE ACUÍFEROS

Modelo Pot
Es el modelo más simple, que puede ser utilizado para estimar el influjo de agua
en un yacimiento.
- Esta basado en la definición básica de compresibilidad.

- Una caída de presión en el yacimiento debido a la producción de fluidos, causa

que el agua del acuífero se expanda y fluya hacia el yacimiento.

- Usualmente se utiliza para acuíferos pequeños, del mismo tamaño del yacimiento.

Aplicando la definición de compresibilidad al acuífero se tiene:

𝑊𝑒 = (𝑐𝑤 +𝑐𝑓 ) ∗ 𝑊𝑖 ∗ (𝑝𝑖 − 𝑝) (1)

 MODELOS DE ACUIEROS

donde:
We: influjo de agua acumulado [MMbbl]
cw: compresibilidad del agua [psi−1]
cf : compresibilidad de la roca [psi−1]
Wi: volumen de agua inicial en el acuıfero [bbl]
pi: presión inicial del yacimiento [psi]
p: presión actual del yacimiento (en el OWC) [psi]

El volumen de agua inicial en un acuífero radial es:

𝜋 𝑟𝑎2 − 𝑟02 ∗ ℎ ∗ ∅
𝑊𝑖 = (2)
5.615
 MODELOS DE ACUÍFEROS

donde:
ra: radio del acuífero [ft]
ro: radio del yacimiento [ft]
h: espesor del acuífero [ft]
θ: porosidad en el acuífero

En el caso que la influencia del acuífero no sea completamente radial, se define un

factor de forma:

𝑊𝑒 = (𝑐𝑤 +𝑐𝑓 ) ∗ 𝑊𝑖 ∗ 𝑓 ∗ 𝑝𝑖 − 𝑝 (3)

donde:
f= θ/360°
MODELOS DE ACUÍFEROS

Figura 1.- Modelo de acuífero radial

 MODELOS DE ACUÍFEROS

MODELO SCHILTHUIS

- El comportamiento de flujo esta descrito por la Ley de Darcy.

- Régimen de flujo en estado estable

La tasa de influjo de agua se puede describir aplicando la Ley de Darcy:

(5) (6)

Donde:
k: permeabilidad del acuífero [md] t: tiempo [d]
h: espesor del acuífero [ft] C: constante de influjo de agua [bbl/d/psi]
ra: radio del acuífero [ft] ro: radio del yacimiento [ft]
 MODELOS DE ACUIFEROS

Integrando obtenemos:

(7)

Utilizando un método de integración numérico obtenemos:

(8)

(9)
 MODELOS DE ACUÍFEROS

donde:
j: paso de tiempo
k: número de intervalos de tiempo

MODELO HURST

- El radio “aparente” del acuífero ra se incrementa con el tiempo.

- La relación adimensional ra/ro se reemplaza por una función que depende del
tiempo.

ra/ro = at
Sustituyendo en la Ec. 5 obtenemos:
 MODELOS DE ACUÍFEROS

(10) (11)

Integrando obtenemos:

(12)

Utilizando un método de integración numérico obtenemos:

(13)
 MODELOS DE ACUÍFEROS

El modelo de acuífero de estado estable de Hurst contienen dos parámetros

desconocidos: a y C. Estos parámetros se pueden determinar a partir del
comportamiento de presión e historia de influjo de agua. Utilizando la Ec. 6 se tiene:

(14)

La Ec. 14 indica que un grafico de pi−p/dWe /dt en función de ln t debe ser una línea
recta con pendiente 1/ C y cuando t = 1 se obtiene 1/ C ln a
MODELOS DE ACUIFEROS

Figura 3. Método pi−p/dWe/dt vs. ln t

 MODELOS DE ACUÍFEROS

MODELO VAN EVERDINGEN Y HURST

(14)

Esto conduce a la determinación del influjo de agua como función de una caída de
presión dada en el borde interno del sistema yacimiento - acuífero.

Figura 4.- Influjo de agua a un yacimiento cilíndrico

 MODELOS DE ACUÍFEROS

Van Everdingen-Hurst propusieron una solución a la ecuación adimensional de

difusividad que utiliza la condición de presión constante y las siguientes condiciones
iníciales y de borde:

- Condición inicial: p = pi, ∆t

- Condición de borde interno: p = (pi − p), r = ro, ∆t
- Condición de borde exterior:
- Acuífero infinito: p = pi, r → ∞
- Acuífero finito:

Adicionalmente, Van Everdingen-Hurst asumieron que el acuífero estaba caracterizado

por:

- Espesor uniforme
- Permeabilidad constante
- Porosidad constante
- Compresibilidad de roca y agua constante
 MODELOS DE ACUÍFEROS

La solución de la ecuación de difusividad para un sistema yacimiento - acuífero,

considerando las condiciones de borde descritas, permite calcular el influjo de agua en
forma de un parámetro adimensional, denominado influjo de agua adimensional (WeD), el
cual es función del tiempo adimensional (tD ) y el radio adimensional (rD).

WeD = f (tD, rD)

WeD se encuentra en forma tabular para diversas geometrías de sistema yacimiento-

acuífero. El influjo acumulado de agua se calcula de la siguiente expresión:

We = B∆pWeD

donde:
We: influjo de agua acumulado [bbl]
B: constante de influjo de agua (depende del modelo geométrico) [bbl/psi]
p = pi − p
WeD: influjo de agua adimensional
 MODELOS DE ACUÍFEROS

El valor de tD y B se muestran a continuación:

𝑟𝑎
𝑟𝐷 =
𝑟𝑒

Donde:
k: permeabilidad del acuífero [md]
t: tiempo [años]
h: espesor del yacimiento [ft]
μw: viscosidad del agua en el acuífero [cp]
ra: radio del acuífero [ft], re: radio del yacimiento [ft]
cw: compresibilidad del agua [psi−1], cr: compresibilidad de la roca [psi−1]
ct = cw + cr: compresibilidad total [psi−1]
L: longitud del acuífero [ft]
W: ancho del yacimiento [ft]
Va: volumen del acuífero [ft3],
 MODELOS DE ACUÍFEROS

PRINCIPIO DE SUPERPOSICIÓN. VAN EVERDINGEN-HURST

- Existe una caída de presión en el contacto agua-petróleo debido a la producción de

fluidos en un yacimiento asociado a un acuífero.

- El agua se expande y la caída de presión se propaga dentro del acuífero hacia el

borde exterior.

- Debido a que las caídas de presión ocurren en forma independiente, el agua se

expande a consecuencia de sucesivas caídas de presión.
 MODELOS DE ACUÍFEROS

Figura 5. Presión en el contacto agua-petróleo

 MODELOS DE ACUÍFEROS

La presión promedio es:

(17)

La caída de presión es:

(18)

Para calcular el influjo acumulado de agua a un tiempo arbitrario, el cual corresponde

al paso de tiempo, se requiere la superposición de las soluciones de la ecuación 19.
 MODELOS DE ACUIFEROS

(19)

Sumando obtenemos:

(20)
 MODELOS DE ACUÍFEROS

Figura 6. Ilustración del principio de superposición

 MODELOS DE ACUÍFEROS

La constante del acuífero B puede ser determinado mediante la solución del método
gráfico de balance de materiales. Para ello se tiene:

(21)

Por lo que:

(22)

La solución de la ecuación lineal de balance de materiales mediante el método gráfico

puede ser utilizada para determinar el valor de un parámetro desconocido del acuífero,
cuando el resto de los parámetros son conocidos.
MODELOS DE ACUÍFEROS
MODELOS DE ACUÍFEROS

2009
 MODELOS DE ACUÍFEROS

Influjo de agua para el Modelo CARTER TRACY

 Up t Dj 1   We t Dj 1 P t Dj  
We t Dj   We t Dj 1    t Dj  t Dj 1  (23)

 P t Dj   t Dj 1 P t DJ  


Donde:

W  influjo acumulado agua

e
U  constante acuífero
p  caída de presión, psi
t  tiempo adimensional
D
 MODELOS DE ACUÍFEROS

Pt D   ao  a1t D  a 2 ln t D  a3 ln t D 

2

Pt D   solución adimens. ec. difusivida d

Pt D   dPt D  / dt D

 kt 
t D  0.00634  
2 
días
 cro  (24)
 kt 
t D  2.309  
2 
años
 cro 
 MODELOS DE ACUÍFEROS

- Acuífero de geometría radial

U = 1.119 fØhcr2 bbl/psi

- Acuífero de geometría lineal

U = 0.1781wLhØc bbl/psi
 MODELOS DE ACUÍFEROS

MODELO FETKOVITCH

 J  pa  p 
dWe
qw  (25)
dt

Donde:

q w  influjo agua
J  indice productividad acuífero
p  presión yacimiento- contactos
p  presión promedio acuífero
 MODELOS DE ACUIFEROS

We 
Wei
 pa  pn 1  e  Jp t
n 1
i n / Wei
 (26)
n
pi
pa pres. prom. acuíf. al tiempo (n  1) th
n

pn p. prom. en límite yac. - acuif. a n th

 n 1

  We 
pa  pi 1  j 1 
j

n 1
 Wei 
 
  (27)
pn 1  pn
p
2
 MODELOS DE ACUÍFEROS

ACUIFERO RADIAL FINITO

7.08 103 f  k  h
J (28)
 re 3 
u w  ln  
 ro 4 