Ing. Franz Gustavo Vargas M.

1. Conceptos generales: para obtener las ecuaciones

se tiene 3 principios básicos:
• Ley de conservación de masas: se considera que alguna cantidad
física es considerada, donde no se crea ni se destruye.
• La ley de Darcy: expresa el hecho de que la rata volumétrica o flujo
por unidad de área de sección transversal o cualquier punto en un
medio poroso es proporcional al gradiente en potencial en la
dirección de flujo a ese punto.
• La ecuación de estado: pueden ser obtenidas de las ecuaciones de
flujo. Una ecuación de estado, especifica la dependencia de la
densidad del fluido sobre la presión de dicho fluido y la
temperatura.
1.1. ley de conservación de masas.-
Masa de entrada MASA Masa de salida

Aportes adicionales

Movimiento en Cantidad de Cantidad de Aportes

el contenido de = masa entrante
- masa saliente
+ adicionales
masa
Propiedades petrofísicas de la ROCA: Porosidad – permeabilidad - saturación
Propiedades petrofísicas del FLUIDO: viscosidad – densidad - °API (instituto
americano del petroleo)
𝑃1 𝑃2

𝑓𝑙𝑢𝑗𝑜 (𝑙𝑎𝑚𝑖𝑛𝑎𝑟)
𝜇
𝑞
𝑃1 > 𝑃2

> °API (liviano) = Bolivia °API = 55-65

< °API (pesado) = Venezuela °API = 10 -20
 1.2. para flujo laminar, se considera (Nre) bajos
donde se puede expresar en forma vectorial, con:
Nre < 2000 = flujo laminar
Nre 2000 – 4000 = comportamiento de transición
Nre > 4000 = flujo turbulento

𝑘
𝑣 =− ∙𝜌∙𝛻∙𝛟
𝜇
𝑣 = velocidad o rata de volumen 𝜌 = 𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 − = 𝒅𝒊𝒓𝒆𝒄𝒄𝒊𝒐𝒏 𝒅𝒆𝒍 𝒇𝒍𝒖𝒋𝒐
𝑘 = permeabilidad 𝛻 = 𝑑𝑖𝑟𝑒𝑐𝑐𝑖𝑜𝑛 𝑑𝑒𝑙 𝑓𝑙𝑢𝑗𝑜
(𝑝𝑎𝑟𝑡𝑒 𝑓𝑖𝑠𝑖𝑐𝑎 𝑑𝑒𝑙 𝑟𝑒𝑠𝑒𝑟𝑣𝑜𝑟𝑖𝑜)
𝜇 = viscosidad
𝛻 ∙ 𝛟 = gradiente potencial en dirección del flujo 𝛟 = potencial del flujo
 Vector : dirección – modulo - sentido
z 𝜌 𝑑𝛟
𝑣𝑧 = − ∙ 𝑘𝑧 ∙
𝜇 𝑑𝑧
𝒌
𝒗=− ∙𝝆∙𝜵∙𝝓
𝝁

𝑥 𝜌 𝑑𝛟
𝑣𝑥 = − ∙ 𝑘𝑥 ∙
𝜇 𝑑𝑥
La ley de Darcy puede expresarse como:
𝝆 𝒅𝝓
𝒗𝒚 = − ∙ 𝒌𝒚 ∙ 𝑘𝑥 𝑑𝑝
𝝁 𝒅𝒚 𝜌 𝒅𝝓 𝜌 1 𝑑𝑝
𝑣𝑥 = − ∙
𝑦 𝑣𝑥 = − ∙ 𝑘𝑥 ∙ 𝑣𝑥 = − ∙ 𝑘𝑥 ∙ ( ∙ ) 𝜇 𝑑𝑥
𝜇 𝒅𝒙 𝜇 𝜌 𝑑𝑥
Pero: 𝒅𝝓 1 𝑑𝑝
= ∙ 𝝆 𝒅𝝓 𝝆 𝟏 𝒅𝒑 𝒌𝒚 𝒅𝒑
𝒅𝒙 𝜌 𝑑𝑥 𝒗𝒚 = − ∙ 𝒌𝒚 ∙ 𝒗𝒚 = − ∙ 𝒌𝒚 ∙ ( ∙ ) 𝒗𝒚 = − ∙
𝝁 𝒅𝒚 𝝁 𝝆 𝒅𝒚 𝝁 𝒅𝒚
𝒅𝝓 𝟏 𝒅𝒑
= ∙
𝒅𝒚 𝝆 𝒅𝒚 𝜌 𝑑𝛟 𝜌 1 𝑑𝑝 𝒌𝒛 𝒅𝒑
𝑣𝑧 = − ∙ 𝑘𝑧 ∙ 𝑣𝑧 = − ∙ 𝑘𝑧 ∙ ( ∙ + 𝑔𝑍) 𝒗𝒛 = − ∙ (∙ + 𝒈𝒁)
𝒅𝝓 𝟏 𝒅𝒑 𝜇 𝜌 𝑑𝑧 𝝁 𝒅𝒛
𝜇 𝑑𝑧
= ∙ + 𝒈𝒁
𝒅𝒛 𝝆 𝒅𝒛 Condiciones: X=i, Z=1
Base para flujo Radial y Caudales Bajos:
- bajas velocidades
- la gravedad puede ser despreciable (g=0)
Flujo Radial 𝑘𝑥 𝑑𝑝
𝑣𝑥 = − ∙
𝜇 𝑑𝑥

Pozo

𝑘𝑟 𝑑𝑝 𝜽
r 𝑣𝑥 = − ∙ r
𝜇 𝑑𝑟
Base para flujo Radial a altas Velocidades: la ecuación de Darcy no es
totalmente valida, se añade una constante de condensación de
velocidad cuadrática.
Donde:
𝐷1 = constante en función (estructura final)
𝑘𝜎 𝑑𝑝
− ∙ = 𝑣 + 𝐷1 ∙ 𝑣 2 𝑠𝑖𝑔𝑚𝑎 𝜎 = 𝑑𝑖𝑟𝑒𝑐𝑐𝑖𝑜𝑛 𝑑𝑒𝑙 𝑓𝑙𝑢𝑗𝑜
𝜇 𝑑𝜎
 1.2.1. la ecuación de continuidad para flujo
𝜌𝒗 + ∆ 𝜌𝑣
lineal.- 𝜌𝑣
𝒛 𝑧
𝑦
z
∆ 𝜌𝑣𝑧 = aportes adicionales

𝜌𝑣𝑥 𝜌𝒗𝒙 + ∆ 𝜌𝑣𝑥

𝒗𝒛 , 𝒗𝒙 , 𝒗𝒚 = corresponde al flujo en las direcciones x, y, z.

𝜌𝒗𝒚 + ∆ 𝜌𝑣𝑦
𝑦

𝜌𝑣𝑧
Fig. 1.1.- elemento de volumen para dar conocimiento a la ecuación de CONTINUIDAD , flujo lineal
 MASA DE FLUJO DE ENTRADA EN LA DIRECCION x :

𝜌𝒗𝒛 + ∆ 𝜌𝑣𝑧 𝜌𝑣𝑥 ∆𝑦 ∆𝑧

𝜌𝑣𝑦
z
MASA DE FLUJO DE SALIDA EN LA DIRECCION x :

𝜌𝒗𝒙 + ∆ 𝜌𝑣𝑥 ∆𝑦 ∆𝑧

𝜌𝑣𝑥 𝜌𝒗𝒙 + ∆ 𝜌𝑣𝑥

𝑥
FLUJO NETO ¨x¨ = Cantidad que entra – Cantidad que Sale

=𝜌𝑣𝑥 ∆𝑦 ∆𝑧 − 𝜌𝒗𝒙 + ∆ 𝜌𝑣𝑥 ∆𝑦 ∆𝑧

𝜌𝒗𝒚 + ∆ 𝜌𝑣𝑦
𝑦 =𝜌𝑣𝑥 ∆𝑦 ∆𝑧 − 𝜌𝒗𝒙 ∆𝑦 ∆𝑧 − ∆ 𝜌𝑣𝑥 ∆𝑦 ∆𝑧

𝜌𝑣𝑧 FLUJO NETO ¨x¨ = −∆ 𝝆𝒗𝒙 ∆𝒚 ∆𝒛 …….(a)

 FLUJO NETO ¨x¨ = −∆ 𝝆𝒗𝒙 ∆𝒚 ∆𝒛 …….(a)

• En forma similar para y, z

FLUJO NETO ¨y¨ = −∆ 𝝆𝒗𝒚 ∆𝒙 ∆𝒛 …….(b)

FLUJO NETO ¨z¨ = −∆ 𝝆𝒗𝒛 ∆𝒙 ∆𝒚 …….(c)

Además por conservación de masa para un incremento, ∆t

Flujo Neto Total:

(∆𝒙 ∆𝒚 ∆𝒛 ) ∙ 𝝓𝝆 − ∆𝒙 ∆𝒚 ∆𝒛 ∙ 𝝓𝝆 … . . (𝒅)
𝒕+∆t 𝒕

Igualando la ecuación (d) con la sumatoria de (a) + (b)+(c) para un ∆t :

− ∆t ∆ 𝜌𝑣𝑥 ∆𝑦 ∆𝑧 + ∆ 𝝆𝒗𝒚 ∆𝒙∆𝒛 + ∆ 𝝆𝒗𝒛 ∆𝒙 ∆𝒚 = (∆𝒙 ∆𝒚 ∆𝒛 ) ∙ 𝝓𝝆 − ∆𝒙 ∆𝒚 ∆𝒛 ∙ 𝝓𝝆

𝒕+∆t 𝒕
− ∆t ∆ 𝜌𝑣𝑥 ∆𝑦 ∆𝑧 + ∆ 𝝆𝒗𝒚 ∆𝒙∆𝒛 + ∆ 𝝆𝒗𝒛 ∆𝒙 ∆𝒚 = (∆𝒙 ∆𝒚 ∆𝒛 ) ∙ 𝝓𝝆 − ∆𝒙 ∆𝒚 ∆𝒛 ∙ 𝝓𝝆
𝒕+∆t 𝒕

Dividiendo entre: ∆𝒙 ∆𝒚 ∆𝒛 ∆t

(∆𝒙 ∆𝒚 ∆𝒛 )∙𝝓𝝆
− ∆t
∆ 𝜌𝑣𝑥 ∆𝑦 ∆𝑧 +∆ 𝝆𝒗𝒚 ∆𝒙 ∆𝒛 +∆ 𝝆𝒗𝒛 ∆𝒙 ∆𝒚
= 𝒕+∆ t− ∆𝒙 ∆𝒚 ∆𝒛 ∙𝝓𝝆
𝒕
∆𝒙 ∆𝒚 ∆𝒛 ∆t ∆ 𝒙 ∆𝒚 ∆ 𝒛 ∆ t

∆ 𝜌𝑣𝑥 ∆ 𝝆𝒗𝒚 ∆ 𝝆𝒗𝒛 𝝓𝝆 𝝓𝝆

− + + = −
∆𝒙 ∆𝒚 ∆𝒛 ∆t 𝒕+∆t
∆t 𝒕

∆ 𝜌𝑣𝑥 ∆ 𝝆𝒗𝒚 ∆ 𝝆𝒗𝒛 𝝓𝝆 𝒕+∆t − 𝝓𝝆 𝒕

− + + =
∆𝒙 ∆𝒚 ∆𝒛 ∆t
En el limite : ∆𝒙 , ∆𝒚 , ∆𝒛 0

𝑑 𝑑 𝑑 𝑑
𝜌𝑣𝑥 + 𝜌𝑣𝑦 + 𝜌𝑣𝑧 = - 𝝓𝝆
𝑑𝑥 𝑑𝑦 𝑑𝑧 𝑑𝑡
 𝑑 𝑑 𝑑 𝑑
𝜌𝑣𝑥 + 𝑑𝑦 𝜌𝑣𝑦 +𝑑𝑧 𝜌𝑣𝑧 = - 𝑑𝑡 𝝓𝝆
𝑑𝑥

En forma vectorial:
𝑑 𝑑 𝑑
(𝑑𝑥 + 𝑑𝑦 + 𝑑𝑧
) =𝛻

𝑑
𝛻 = 𝜌𝑣 =- 𝝓𝝆
𝑑𝑡
 1.3 LA ECUACIÓN DE CONTINUIDAD PARA FLUJO
RADIAL:
Flujo neto = Flujo de Entrada – Flujo de salida

𝜌𝑣𝑟
→ Flujo de Entrada = 𝜃 𝑟 + ∆𝑟 ℎ(𝜌𝑣𝑟 )

→ Flujo de salida = 𝜃𝑟ℎ 𝜌𝑣𝑟 + ∆ 𝜌𝑣𝑟

h Donde:
𝜌𝑣𝑟 + ∆ 𝜌𝑣𝑟
𝜌𝑣𝑟 = Caudal por unidad de área
𝑣𝑟 = Flujo Radial

r + Δr

r
Figura. 1.2 Elemento de Volumen para dar Conocimiento a la
Ecuación de Continuidad para Flujo Radial

Flujo Neto a un Tiempo (Δt):

−∆𝑡 𝜽 𝒓 + ∆𝒓 𝒉 𝝆𝒗𝒓 − 𝜽𝒓𝒉 𝝆𝒗𝒓 + ∆ 𝝆𝒗𝒓 = ɸ𝜌ℎ𝜃𝑟∆𝑟ȁ ⋮𝑡+∆𝑡 −ɸ𝜌ℎ𝜃𝑟∆𝑟ห⋮𝑡

 Dividiendo entre el Angulo: 𝜃𝑟∆𝑟 ℎ∆𝑡

+∆𝑡 𝜃𝑟ℎ 𝜌𝑣𝑟 + 𝜃∆𝑟ℎ 𝜌𝑣𝑟 − 𝜃𝑟ℎ𝜌𝑣𝑟 − 𝜃𝑟ℎ∆(𝜌𝑣𝑟 ) = ɸ𝜌ℎ𝜃𝑟∆𝑟ȁ ⋮𝑡+∆𝑡 −ɸ𝜌ℎ𝜃𝑟∆𝑟ห⋮𝑡

+∆𝑡𝜃𝑟ℎ 𝜌𝑣𝑟 + ∆𝑡𝜃∆𝑟ℎ 𝜌𝑣𝑟 − ∆𝑡𝜃𝑟ℎ𝜌𝑣𝑟 − ∆𝑡𝜃𝑟ℎ∆ 𝜌𝑣𝑟 = ɸ𝜌ℎ𝜃𝑟∆𝑟ȁ ⋮𝑡+∆𝑡 − ɸ𝜌ℎ𝜃𝑟∆𝑟ห⋮𝑡

∆𝑡𝜃𝑟ℎ 𝜌𝑣𝑟 ∆𝑡𝜃∆𝑟ℎ 𝜌𝑣𝑟 ∆𝑡𝜃𝑟ℎ𝜌𝑣𝑟 𝜃𝑟ℎ∆ 𝜌𝑣𝑟 ∆𝑡 ɸ𝜌ℎ𝜃𝑟∆𝑟 ɸ𝜌ℎ𝜃𝑟∆𝑟

+ − − = ቤ⋮ − ቤ⋮
𝜃𝑟ℎ ∆𝑟 ∆𝑡 𝜃𝑟ℎ ∆𝑟 ∆𝑡 𝜃𝑟ℎ ∆𝑟 ∆𝑡 𝜃𝑟ℎ ∆𝑟 ∆𝑡 𝜃𝑟ℎ∆𝑟∆𝑡 𝑡+∆𝑡 𝜃𝑟ℎ∆𝑟∆𝑡 𝑡

𝜌𝑣𝑟 𝜌𝑣𝑟 𝜌𝑣𝑟 ∆ 𝜌𝑣𝑟 ɸ𝜌 ɸ𝜌

+ − − = ቤ⋮ − ቤ⋮
∆𝑟 𝑟 ∆𝑟 ∆𝑟 ∆𝑡 𝑡+∆𝑡 ∆𝑡 𝑡

𝜌𝑣𝑟 ∆ 𝜌𝑣𝑟 ɸ𝜌 ɸ𝜌
− = ቤ ⋮𝑡+∆𝑡 − ቤ ⋮𝑡
𝑟 ∆𝑟 ∆𝑡 ∆𝑡
1 ∆(ɸ𝜌)
𝜌𝑣𝑟 ∆𝑟 − 𝑟∆ 𝜌𝑣𝑟 =−
∆ 𝜌𝑣𝑟 𝜕 𝜌𝑣𝑟 𝑟∆𝑟 ∆𝑡
𝑠𝑖: → −
∆𝑟 𝜕𝑟 entonces:
𝑠𝑖 𝑒𝑛 𝑒𝑙 𝑙𝑖𝑚𝑖𝑡𝑒 ∶ 𝟏 𝛛 𝛛
∆𝑟 → 0 × 𝒓𝝆𝒗𝒓 = − ɸ𝝆
𝒓 𝛛𝒓 𝛛𝒕
∆𝑡 → 0 “La ecuación de continuidad para flujo radial”
1.7 Soluciones para Flujo Radial para Fluidos de Compresibilidad
Pequeña Constante
→ RECORDANDO : Las Ecuaciones: (1.10) , ( 1.11)

𝜕 2 𝜌 𝜕 2 𝜌 𝜕 2 𝜌 ∅𝜇𝑐 𝜕𝜌
+ + = ∙ (1.10)
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2 𝑘 𝜕𝑡

𝜕 2 𝜌 1 𝜕𝜌 ∅𝜇𝑐 𝜕𝜌
+ = ∙ (1.11)
𝜕𝑟 2 𝑟 𝜕𝑟 𝑘 𝜕𝑡

Consideraciones para Resolver la Ecuación (1.10)

 Flujo radial en todo el reservorio
 Medio poroso e isotropico
 Espesor uniforme
 ø y k constantes e independientes de la presión
 Compresibilidad pequeña y ctte
 Efectos de gravedad despreciables (g=0)
CASOS:
1) Reservorio infinito 𝑟𝑒 > 𝑟𝑤

2) Reservorio limitado , pozo en el centro de un reservorio cilíndrico ,sin flujo hacia el limite exterior

3) Presión ctte fuera del reservorio.

Esquemáticamente:
1) caso de reservorio infinito
𝑝 → 𝑝𝑖

. 𝑟→∞
2) caso de reservorio circular limitado
𝜕𝑝
⋮ =0
𝜕𝑟 𝑟𝑒

3) caso de presión constante fuera del reservorio

𝑝 → 𝑝𝑖
𝑟 → 𝑟𝑒
𝜕𝑝 𝑞𝑢 1
⋮𝑟𝑤 = ∙
𝜕𝑟 2𝜋𝑘ℎ 𝑟𝑤
si 𝒓 = 𝒓𝒆 ; la ecuación se puede escribir de la siguiente manera:

2𝜋𝑘ℎ 𝑑𝑝
𝑞= 𝑟
𝜇 𝑑𝑟 𝑟𝑤

Para un flujo constante, el gradiente de presión será:

𝑑𝑝 𝑞𝑢 1
= ∙
𝑑𝑟 𝑟𝑤
2𝜋𝑘ℎ 𝑟𝑤

Si no Existe flujo en el límite externo

Velocidad de flujo = 0
𝑟 = 𝑟𝑒
Gradiente de presión = 0

𝑑𝑝
=0
𝑑𝑟 𝑟𝑒
 SOLUCIONES MATEMATICA PARA CADA UNO DE LOS CASOS:
1) CASO DE UN RESERVORIO INFINITO

𝑞𝑢 1 ∅ 𝜇 𝑐 𝑟2
𝑃 𝑟,𝑡 = 𝑃𝑖 − − ∙ 𝐸𝑖 −
2𝜋𝑘ℎ 2 4𝑘𝑡 T.P # 1
2) CASO : RESERVORIO CIRCULAR
Donde: (LIMITADO)
∞ −𝑢
𝑒 3) CASO: RESERVORIO DE PRESIÓN
𝐸𝑖 −𝑥 = න 𝑑𝑢 Ei=integral Exponencial CONSTANTE EXTERNA
𝑥 𝑢

𝑥 𝑥2 𝑥3
𝐸𝑖 −𝑥 = 𝐼𝑛 𝑥 − + − +⋯
log(𝑥) 1! 2 ∙ 2! 3 ∙ 3!

Cuando: x < 0.01

1
𝐸𝑙 𝑣𝑎𝑙𝑜𝑟: −𝐸𝑖 −𝑥 ≅ − 𝑙𝑛 ℘𝑥 = 𝑙𝑛 − 0.5772
𝑥

𝐸𝑖 −𝑥 𝐷𝑜𝑛𝑑𝑒: ℘𝑒𝑠 𝑙𝑎 𝑐𝑡𝑡𝑒 𝑑𝑒 𝐸𝑢𝑙𝑒𝑟`𝑠 = 1.78

Cuando: x > 0.01
 4𝑘𝑡
𝑝𝑎𝑟𝑎 ∶ > 100
𝜙𝜇𝐶𝑟 2

𝑞𝜇 ℘𝜙µ𝐶𝑟 2 𝑞𝜇 𝑘𝑡
𝑃 𝑟,𝑡 = 𝑃𝑖 + ∙ 𝑙𝑛 ó 𝑃 𝑟,𝑡 = 𝑃𝑖 +
4𝜋𝑘ℎ
∗ 𝑙𝑛
𝜙𝜇𝐶𝑟 2
+ 0.80907
4𝜋𝑘ℎ 4𝑘𝑡

𝑝𝑎𝑟𝑎: 𝑟 = 𝑟𝑤 (hablamos de una presión de pozo)

Entonces la presión de fondo fluyente será:

𝑞𝜇 ℘𝜙µ𝐶𝑟𝑤2 𝑞𝜇 𝑘𝑡
𝑃𝑤𝑓 = 𝑃𝑖 + 𝑙𝑛 𝑃𝑤𝑓 = 𝑃𝑖 − 𝑙𝑛 + 0.80907
4𝜋𝑘ℎ 4𝑘𝑡 4𝜋𝑘ℎ 𝜙µ𝐶𝑟𝑤2

Donde: C = Compresibilidad de la roca (Psi-1) caída de presión:

𝑃𝑤𝑓 =Presión de fondo Fluyente (Psia) ΔP = 𝑃𝑖 – 𝑃𝑤𝑓

h = Espesor de la arena productiva
𝑃𝑖 = Presión Inicial (Psia)
Ejemplo 1.1 se tiene un reservorio con las siguientes características

µo = 0.72 cp h = 15 pies Ø = 23.4 % k = 500 md

Bo = 1.475 bbl/bls Cf = 15x10-6 Psi-1 𝑃𝑖 = 3000 psia

Después de que el pozo produzca 200 bls/dia durante 10 días, Calcular la presión de fondo fluyente para un radio de
1000 pies y la caída de presión.

Solución: 𝑞𝜇 𝑘𝑡
𝑃𝑤𝑓 = 𝑃𝑖 − 𝑙𝑛 + 0.80907
q = 200bls/día 4𝜋𝑘ℎ 𝜙µ𝐶𝑟𝑤2

t = 10 días Remplazamos:
200𝐵𝑙𝑠
r= 1000 pies 0,72 𝐶𝑝 0,5 𝑑𝑎𝑟𝑐𝑦 10 𝑑𝑖𝑎𝑠
𝑑í𝑎
𝑃𝑤𝑓 = 3000 𝑝𝑠𝑖𝑎 − ∙ 𝑙𝑛 + 0.80907
incognitas: 4𝜋 0,5 𝑑𝑎𝑟𝑐𝑦 15 𝑝𝑖𝑒𝑠 0,234 0,72 𝐶𝑝 15𝑥10−6 𝑃𝑠𝑖 −1 1000 𝑝𝑖𝑒𝑠 2

𝑃𝑤𝑓 = ???
1,5278 1,8968
ΔP = ??? 𝑃𝑤𝑓 = 3000 𝑝𝑠𝑖𝑎 − (1,5278 ∙ 1,8968) = 2997.1 psia.

ΔP = 𝑃𝑖 – 𝑃𝑤𝑓 ΔP = (3000 – 2997.1) psia = 2.9 psia

1.8 Conceptos de Comportamiento de flujo Transiente,
Estable y Semiestable.
Flujo Transiente

Flujo estable (POST-TRANSIENTE)

Pwf

Flujo Semi-estable(Pseudoestable)

t
Declinación de las presiones reservorio Circular, Flujo constante
1.9 El principio de superposición. -
toma en cuenta el caudal variable o rata de flujo con respecto al
tiempo. El principio de superposición es un modelo matemático que
facilita la readecuación de las ecuaciones para caudal múltiple
partiendo de las correspondientes ecuación o caudal constante.
a) Caso de Comportamiento de Varios Pozos
Pozo de observación o monitor

𝒓𝟑
𝒓𝟏
𝒓𝟐
Pozo 3
Pozo1 ∆𝑷𝟑
Pozo2 𝒒𝟑
∆𝑷𝟏
∆𝑷𝟐
𝒒𝟏 𝒒𝟐

Reservorio con 3 pozos fluyendo

Aplicando el principio de superposición será:

∆𝑃𝑡 = ∆𝑃1 + ∆𝑃2 + ∆𝑃3

Y como calcularemos cada caída de presión considerando un

reservorio infinito en unidades de campo.

70.6 𝑞 𝜇 𝐵𝑜 Ø𝑢𝐶𝑟²
∆𝑃 = 𝑃𝑖 – 𝑃(𝑟, 𝑡) = ∙ −𝐸𝑖 ( − )
𝐾ℎ 0,00105 𝑘𝑡

Para “n” pozos fluyendo:

𝑛

∆𝑃𝑡 = ෍ ∆𝑃𝑖
𝑖=1
 Ejemplo 1.2 se tiene 4 pozos fluyendo durante 10 días, calcular la caída de presión con respecto a un pozo
de observación según la Figura Posteriormente descrita. 𝑞₁ = 260𝑏𝑙𝑠/𝑑𝑖𝑎 𝑞₄ = 250𝑏𝑙𝑠/𝑑𝑖𝑎
(𝑋₁, 𝑌₁) 𝑋₄, 𝑌₄
µo = 0.35 cp h = 37 pies Ø = 14 % 1 4
Bo = 1.35 bbl/bls Cf = 12x10-6 Psi-1 k = 40 md UTM

SOLUCIÓN 5 𝑋₅, 𝑌₅
24 ℎ𝑜𝑟𝑎𝑠
10 𝑑𝑖𝑎𝑠 ∗ = 240 ℎ𝑜𝑟𝑎𝑠
1 𝑑𝑖𝑎
2 3
70.6 𝑞 𝜇 𝐵𝑜 Ø𝑢𝐶𝑟² 𝑋₃, 𝑌₃
∆𝑃 = 𝑃𝑖 – 𝑃(𝑟, 𝑡) = ∙ −𝐸𝑖 ( − ) (𝑋₂, 𝑌₂)
𝐾ℎ 0,00105 𝑘𝑡 𝑞₃ = 290𝑏𝑙 𝑠Τ𝑑 𝑖𝑎
𝑞₃ = 280𝑏𝑙 𝑠Τ𝑑 𝑖𝑎
𝑏𝑙𝑠
70.6 𝑞𝑢 𝐵𝑜 70.6 260 ∗ 0.35 𝑐𝑝 ∗ (1,35𝑏𝑏𝑙/𝑏𝑙𝑠൰
𝑑𝑖𝑎 5.86
= =
𝑘ℎ 40 𝑚𝑑) ∗ (37 𝑝𝑖𝑒𝑠

∅𝑢 𝑐𝑟 2 0.14 ∗ 0,35 𝑐𝑝 ∗ 12𝑥10−6 𝑃𝑠𝑖 −1 ∗ 1500 𝑝𝑖𝑒𝑠 2

− =− = −0, 13125 x > 0.01
0.00105𝑘𝑡 0.00105 ∗ 40𝑚𝑑 ∗ 240 ℎ𝑜𝑟𝑎𝑠
 0.13125 − 0.1
𝑦= ∗ 1.22265 − 1.82292 + 1.82292
0.2 − 0.1

𝑦 = 1.63533 = 𝐸𝑖(−𝑥 )

70.6 𝑞 𝜇 𝐵𝑜 Ø𝑢𝐶𝑟²
Δ𝑃1 = ∙ −𝐸𝑖 ( − )
𝐾ℎ 0,00105 𝑘𝑡

Cuando: x > 0.01

Δ𝑃1 = 5.86 ∗ (1.63533) = 9.58 𝑝𝑠𝑖
x = −0, 13125
Δ𝑃2 =

Δ𝑃3 =

Δ𝑃4 =

Δ𝑃total =
 b) Caso del Comportamiento de un Pozos.- ∆𝑃 = ∆𝑃1 + ∆𝑃2 + ∆𝑃3 +….
𝑡𝑜𝑡𝑎𝑙
q
𝑏𝑙𝑠
𝑞3  Para 𝟎 < 𝒕 < 𝒕𝟏
( )
𝑑𝑖𝑎
(𝑞3 − 𝑞2 ) 𝑞1 𝜇𝐵0
𝑞2 ∆𝑃 𝑡 = −𝐸𝑖 −𝑥
2𝜋𝑘ℎ
(𝑞2 − 𝑞1 ) Para (t)
𝑞1 Cambio de choke
 Para 𝒕𝟏 < 𝒕 < 𝒕𝟐

t 𝑞1 𝜇𝐵0 𝑞2 −𝑞1 𝜇𝐵0

𝑡1 𝑡2 𝑡3 𝑡4 ∆𝑃(𝑡) = [−𝐸𝑖 −𝑥 ] + [−𝐸𝑖 −𝑥 ]
2𝜋𝑘ℎ 2𝜋𝑘ℎ

𝑃𝑖 Para (𝑡 − 𝑡1 )
Δ𝑃1
 Para 𝒕𝟐 < 𝒕 < 𝒕𝟑
𝑃𝑤𝑓
Δ𝑃2
𝑞3 − 𝑞2 𝜇𝐵0
∆𝑃 𝑡 = 𝑖𝑑 + 𝑖𝑑 + + [−𝐸𝑖 −𝑥 ቉
Δ𝑃3
2𝜋𝑘ℎ

Para (𝑡 − 𝑡2 )

t
𝑡1 𝑡2 𝑡3 𝑡4
HISTORIAL DE PRODUCCIÓN Y PRESIÓN PARA UN POZO
Generalizado:

𝑛
𝑞1 𝜇𝐵0 ൫𝑞𝑖 − 𝑞𝑖−1 )𝜇𝐵0
∆𝑃 𝑡 = −𝐸𝑖 −𝑥 +෍ −𝐸𝑖 −𝑥
2𝜋𝑘ℎ 2𝜋𝑘ℎ
𝑖=2

ECUACIÓN GENERAL DE SUPERPOSICIÓN PARA UN POZO CON “Producción Variable”

Para un reservorio infinito de unidades de campo:

70.6 𝑞1 𝜇 𝐵𝑜 Ø𝑢𝐶𝑟 2 70.6 𝑞2 − 𝑞1 𝜇 𝐵𝑜 Ø𝑢𝐶𝑟 2

∆𝑃 = ∙ −𝐸𝑖 − + ∙ −𝐸𝑖 − +
𝑘ℎ 0,00105 𝑘𝑡 𝑘ℎ 0,00105 𝑘 (𝑡 − 𝑡1 )

70.6 𝑞3 − 𝑞2 𝜇 𝐵𝑜 Ø𝑢𝐶𝑟 2
+ ∙ −𝐸𝑖 −
𝑘ℎ 0,00105 𝑘(𝑡 − 𝑡2 )
Ejemplo 1.3.- Un pozo petrolífero a estado fluyendo durante 15 días a un caudal de 300
bls/día, entonces se incrementa la producción a 400 bls/día y se mantiene este caudal por
20 días para luego disminuir 280 bls/día con un tiempo de 15 días.
Calcular la caída de presión en el pozo al final de la 3er etapa, considerar ala reservorio como
𝑞 (𝑏𝑙𝑠/𝑑𝑖𝑎)
infinito.
𝑞2
400
DATOS:
300 𝑞1 𝑞3
rw = 0.25 pies 280

200
µ = 0.50 cp.
Φ = 17% 𝑡 (𝑑𝑖𝑎𝑠)
0 15 35 50
15 15 𝑑𝑖𝑎𝑠
Bo = 1.51 bbl/bls
𝑃𝑤𝑓
H = 40 pies
𝑃𝑖
Ct = 17.2 x10-6 psi-1
K = 35 md

𝑡 (𝑑𝑖𝑎𝑠)
0 15 35 50
15 𝑑𝑖𝑎𝑠 20 𝑑𝑖𝑎𝑠 15 𝑑𝑖𝑎𝑠
 24 ℎ𝑜𝑟𝑎𝑠
50 𝑑𝑖𝑎𝑠 ∙ = 1200 ℎ𝑜𝑟𝑎𝑠
1 𝑑𝑖𝑎
24 ℎ𝑜𝑟𝑎𝑠
Solución: 35 𝑑𝑖𝑎𝑠 ∙ 1 𝑑𝑖𝑎
= 400 ℎ𝑜𝑟𝑎𝑠

Tiempo Total = (15 + 20 + 15) días = 50 días 15 𝑑𝑖𝑎𝑠 ∙

24 ℎ𝑜𝑟𝑎𝑠
= 360 ℎ𝑜𝑟𝑎𝑠
1 𝑑𝑖𝑎

70.6 𝑞1 𝜇 𝐵𝑜 ∅ 𝜇 𝑐 𝑟𝑤 2 70.6 𝑞2 − 𝑞1 𝜇𝐵𝑜 ∅ 𝜇 𝑐 𝑟𝑤 2

∆P = −𝐸𝑖 − + + −𝐸𝑖 −
𝑘ℎ 0.00105 𝑘 𝑡 𝐾ℎ 0.00105 𝑘 𝑡 − 𝑡1
70.6 𝑞3 − 𝑞2 𝜇𝐵𝑜 ∅ 𝜇 𝑐 𝑟𝑤 2
+ −𝐸𝑖 −
𝐾ℎ 0.00105 𝑘 𝑡 − 𝑡2

70.6 𝑞1 𝜇 𝐵𝑜 ∅ 𝜇 𝑐 𝑟𝑤 2 70.6 𝑞2 − 𝑞1 𝜇𝐵𝑜 ∅ 𝜇 𝑐 𝑟𝑤 2

∆P = −𝐸𝑖 − + + −𝐸𝑖 −
𝑘ℎ 0.00105 𝑘 𝑡 𝐾ℎ 0.00105 𝑘 𝑡 − 𝑡1
70.6 𝑞3 − 𝑞2 𝜇𝐵𝑜 ∅ 𝜇 𝑐 𝑟𝑤 2
+ −𝐸𝑖 −
𝐾ℎ 0.00105 𝑘 𝑡 − 𝑡2

𝟕𝟎. 𝟔 𝝁 𝑩𝒐 70.6 0.50𝑐𝑝 1.51 𝑏𝑏𝑙/𝑏𝑙𝑠

= = 0.03807
𝒌𝒉 35 𝑚𝑑 40 𝑝𝑖𝑒𝑠
 ∅ 𝝁 𝑪 𝒓𝟐 0.17 0.50 𝑐𝑝 17.2 ∗ 10−6 𝑝𝑠𝑖 −1 0.25𝑝𝑖𝑒𝑠 2
= = 2.48 ∗ 10−6
𝟎. 𝟎𝟎𝟏𝟎𝟓 𝒌 0.00105 35 𝑚𝑑

𝑏𝑏𝑙 2.48 ∗ 10−6 𝑏𝑏𝑙 𝑏𝑏𝑙 2.48 ∗ 10−6

∆𝑃 = 0.03807 300 −𝐸𝑖 − + 400 − 300 −𝐸𝑖 −
𝑑𝑖𝑎 1200 ℎ𝑜𝑟𝑎𝑠 𝑑𝑖𝑎 𝑑𝑖𝑎 1200 − 360 ℎ𝑜𝑟𝑎𝑠

𝑏𝑏𝑙 𝑏𝑏𝑙 2.48 ∗ 10−6

+ 280 − 400 −𝐸𝑖 −
𝑑𝑖𝑎 𝑑𝑖𝑎 1200 − 840 ℎ𝑜𝑟𝑎𝑠

∆𝑃 = 0.03807 ∙ 300 −𝐸𝑖 − 2.06 x10−9 + 100 −𝐸𝑖 − 2.95 ∗ 10−9 + −200 −𝐸𝑖 − 6.88 ∗ 10−9

Como x < 0.01, usaremos:

1
−𝐸𝑖 −𝑥 ≅ − ln 𝛿𝑥 = ln − 0.5775
𝑥
 1
−𝐸𝑖 −𝑥 = ln − 0.5775 = 19.42
2.06 ∗ 10−9

1
−𝐸𝑖 −𝑥 = ln − 0.5775 = 19.06
2.95 ∗ 10−9

1
−𝐸𝑖 −𝑥 = ln − 0.5775 = 18.21
6.88 ∗ 10−9

∆𝑃 = 0.03807 ∙ 300 (19.42) + 100 19.06 − 200 (18.21)

∆𝑃 = 0.03807 5546.8

∆𝑃 = 211.2 𝑝𝑠𝑖 ∙ 𝑏𝑏𝑙 𝑐r. 𝑎 𝑐𝑜𝑛𝑑𝑖𝑐𝑖𝑜𝑛𝑒𝑠 𝑑𝑒 reservorio

1𝑏𝑙𝑠
∆𝑃 = 211.2 𝑝𝑠𝑖 ∗ = 139.84 𝑝𝑠𝑖 ∙ 𝑏𝑙𝑠 𝑐𝑠. 𝑎 𝑐𝑜𝑛𝑑𝑖𝑐𝑖𝑜𝑛𝑒𝑠 𝑑𝑒 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒
1.51𝑏𝑏𝑙
1.10 SISTEMA DE UNIDADES DE DARCY Y DE CAMPO

PARAMETRO O UNIDADES UNIDADES DE

VARIABLE DARCY CAMPO
Ct Vol/vol/atm Vol/vol/psi
∅ Fracción Fracción
h cm pies
k Darcies md
𝜇 cp cp
P atm psi
q cm3/seg bls/dia
r cm pies
t Seg horas
 Termino Termino
Unidades Darcy unidades practicas
𝑘𝑡 𝑘𝑡
𝑡𝐷 = 𝑡𝐷 = 0,000264
Ø 𝜇 𝐶 𝑟² Ø 𝜇 𝐶 𝑟𝑤2
𝑞𝜇 70.6 𝑞 𝜇 𝛽𝑂
4𝜋−𝑘ℎ 𝑘ℎ
162.6𝑞𝜇𝐵𝑜 0.000264 𝑘 𝑡 162.6 𝑞𝜇𝐵𝑜 𝑘𝑡
𝑃𝑖 − 𝑃𝑤𝑓 = {𝑙𝑜𝑔 + 0.351} 𝑃𝑖 − 𝑃𝑤𝑓 = {𝑙𝑜𝑔 – 3.23}
𝐾ℎ Ø 𝜇 𝑐 𝑟𝑤2 𝑘ℎ Ø 𝜇 𝑐 𝑟𝑤2
1.11 SUPERPOSICIÓN EN RESERVORIOS LIMITADOS O PARCIALMENTE LIMITADOS.

Falla

Pozo Pozo Imagen

La presión en el
pozo es
influenciada por
la falla ,cuando el
pozo fluye envía
un pulso u onda
𝜕𝑃
=0 de presión ,este
𝜕𝑟
pulso es reflejada
No existe flujo en el limite por la falla .
La onda de presión reflejada puede ser visualizada como si emanara de una
imagen del pozo situada a una distancia “r” , igual detrás de la falla entre el
pozo actual y el pozo imagen existiera una distancia de 2r .
Pozo Imagen 1
(∆𝑃)1

𝑟1

𝑟1 El concepto de pozo
Pozo Actual
imagen puede ser
𝑟𝑤 = 𝑟𝑎𝑑𝑖𝑜 𝑑𝑒 𝑝𝑜𝑧𝑜 𝑟2 extendido para
calcular el
𝑟2 comportamiento de
Pozo Imagen 2 presión de un pozo
(∆𝑃)2
localizado dentro de
una serie de fallas.
• La caída de presión total será: la caída de presión de ida ala propia
producción del pozo mas las caídas adicionales de presión causadas
por el o los pozos imagen.
𝑛

∆𝑃 𝑇𝑂𝑇𝐴𝐿 = 𝑃𝑖 − 𝑃𝑤𝑓 = ∆𝑃 𝑝𝑜𝑧𝑜 𝑎𝑐𝑡𝑢𝑎𝑙 + ෍ ቀ∆𝑃)𝑝𝑜𝑧𝑜 𝑖𝑚𝑎𝑔𝑒𝑛 (𝑖)

𝑖=1

162.6 𝑞𝑜 𝐵𝑜 𝜇𝑜 𝑘𝑡 70.6 𝑞𝑜 𝐵𝑜 𝜇𝑜 Ø 𝜇 𝑐𝑡 2 𝑟1 2
∆𝑃 𝑇𝑂𝑇𝐴𝐿 = × 𝑙𝑜𝑔 ( ) – 3,23 + −𝐸𝑖 −
𝑘ℎ Ø𝜇𝑜 𝑐 𝑟𝑤2 𝑘ℎ 0,00105 𝑘𝑡

∆𝑃𝑤=Caída de presión de pozo actual ∆𝑃i= Caída de presión imagen (i) Donde:
q: bls/día
Considerando solo el efecto del pozo sin tomar en cuenta los pozos T: horas
imagen. Se tendrá: K: md
∆𝑃𝑤 = 𝑃𝑖 − 𝑃𝑤 h: pies
r: pies
PRACTICA N°1: Un pozo a estado fluyendo por 12 días con 250 bls / día,
este pozo está localizado a 550 pies y 650 pies de distancia de 2 fallas
sellantes considerando los siguientes datos:

ɸ= 0,19 pi = 4000 psia

Ct = 1,80E-05 psi-1 μ= 1,2 cp
k= 70 md
rw = 0,28 pies
h= 30 pies
Bo= 1.33 bbl/bls

a) Calcular el PERFIL del pozo como función del tiempo según formato adjunto.
b) Graficar Pwf vs t en escala semi – log mostrando las pendientes para cada una de las fallas, calcular las
pendientes m1 y m2 en psi/ciclo.
c) Analizar los resultados
• Variaciones de ∆𝑃𝑤 y ∆𝑃𝑡𝑜𝑡𝑎𝑙 en función al tiempo 𝑓(𝑡)
• Resultados de ∆𝑃1 y ∆𝑃2 en funciones al tiempo
• ∆𝑃𝑤𝑓 y 𝑃𝑤 en funciones al tiempo
• Incidencias de las fallas
 a) Calcular el PERFIL del pozo como función del tiempo según formato adjunto.

Tiempo Δpw -x1 -Ei(-x1) Δp1 -x2 -Ei(-x) Δp2 Δp total Pwf Pw
∆𝑃𝑇𝑂𝑇𝐴𝐿 = ∆𝑃𝑤 + ∆𝑃1 + ∆𝑃2
(horas) (psi) (psi) (psi) (psi) (psi) (psi)

2 167,09 - 33,781 0.00000 0.00000 -47,182 0.00000 0.00000 167,09 3832.91 3832.91

∆𝑃𝑇𝑂𝑇𝐴𝐿 = 𝑃𝑖 − 𝑃𝑤𝑓

𝑃𝑤𝑓 = 𝑃𝑖 − ∆𝑃𝑇𝑂𝑇𝐴𝐿

∆𝑃𝑤 = 𝑃𝑖 - 𝑃𝑤

𝑃𝑤 = 𝑃𝑖 − ∆𝑃𝑤
Solución.-

Caída de presión de pozo actual Varian respecto a la FALLA Varian respecto a la FALLA
AP imagen 1 AP imagen 2

Tiempo = 2 horas
Caída de presión de pozo actual ɸ= 0,19
pi = 4000 psia
Ct = 1,80E-05 psi-1
μ= 1,2 cp
k= 70 md
30,894 8,639 - 3.23 = 167,09343 rw = 0,28 pies
h= 30 pies
Bo= 1.33 bbl/bls
r1= 550 pies
r2= 650 pies
qsc= 250 bls/d
Varian respecto a la FALLA
AP imagen 1

= 0.00000

>0,01 POR TABLA

13,414 -33,781224 <0,001 ECUACION EULER

Varian respecto a la FALLA = 0.00000

AP imagen 2

= 0.00000

>0,01 POR TABLA ɸ= 0,19

13,414 -47,182041 <0,001 ECUACION EULER pi = 4000 psia
Ct = 1,80E-05 psi-1
μ= 1,2 cp
k= 70 md
= 0.00000 rw = 0,28 pies
h= 30 pies
Bo= 1.33 bbl/bls
r1= 550 pies
r2= 650 pies
qsc= 250 bls/d
Tiempo Δpw -x1 -Ei(-x1) Δp1 -x2 -Ei(-x) Δp2 Δp total Pwf Pw
(horas) (psi) (psi) (psi) (psi) (psi) (psi)

0 ∆𝑃𝑇𝑂𝑇𝐴𝐿 = ∆𝑃𝑤 + ∆𝑃1 + ∆𝑃2

0,1
0,2
0,3
0,5
0,7
∆𝑃𝑇𝑂𝑇𝐴𝐿 = 𝑃𝑖 − 𝑃𝑤𝑓
0,9
𝑃𝑤𝑓 = 𝑃𝑖 − ∆𝑃𝑇𝑂𝑇𝐴𝐿
1
2 167,09 - 0.00000 0.00000 - 0.00000 0.00000 167,09 3832.91 3832.91
33,781 47,182

3
4 ∆𝑃𝑤 = 𝑃𝑖 - 𝑃𝑤
8
10
32
𝑃𝑤 = 𝑃𝑖 − ∆𝑃𝑤
56
80
104
128
152
200
224
248
272
296
b.- Graficar Pwf vs t en escala semi – log mostrando las pendientes para cada una de las fallas, calcular
las pendientes m1 y m2 en psi/ciclo.

m1

m2

ciclo ciclo ciclo ciclo