Unidad 1
Unidad 1
Unidad 1
Villahermosa
Semestre: Séptimo.
Bibliografía. ......................................................................................................................... 30
1.1 Naturaleza del flujo en yacimiento.
Los depósitos de petróleo contienen gas y agua, o aceite y agua, con el petróleo
generalmente en las partes centrales de los poros y el agua en anillos pendulares alrededor de
los contactos de los granos. Esta agua es aparentemente inmóvil y la saturación de agua
irreducible suele ser del 20 al 40% del volumen de los poros.
Las presiones en el petróleo son más altas que las presiones hidrostáticas normales
del agua debido a su menor densidad de peso.
I. Ecuación de continuidad.
a. Conservación de masa.
b. Conservación de energía.
c. Conservación de momento.
II. Ecuación de transporte.
III. Ecuación de estado.
lim ∆𝑦 → 0
lim ∆𝑧 → 0
lim ∆𝑡 → 0
Considerando que un solo fluido de densidad “ρ“, pasa a través de cada una de las caras de
un volumen elemental de porosidad Φ a una velocidad “v“, se tiene un flujo másico
por unidad de superficie tal como lo muestra la ecuación:
𝑀 𝐿 𝑀
𝑝𝑣 = 𝑃 [ 3 ] 𝑣 [ ] = 𝑝𝑣 [ 2 ]
𝐿 𝑇 𝐿𝑇
Volumen elemental respectivo.
Conservación de masa.
De donde:
ρ es la densidad de un fluido.
1 𝜕(𝑟𝜌𝑣𝑟 ) 𝜕(𝜌∅)
=
𝑟 𝜕𝑟 𝜕𝑡
𝑘𝜌
𝑉= ∆Φ
𝜇
Donde:
Φ: potencial.
𝜌
𝑑𝜌
Φ= ∫ + 𝑔𝑧
𝜌
𝜌0
Según la definición anterior podemos deducir cual es la ecuación de Darcy para las
direcciones x, y, z, y pueden expresarse como sigue:
𝑘𝑥 𝑑𝑝
𝑣𝑥 = −
𝜇 𝑑𝑥
𝑘𝑦 𝑑𝑝
𝑣𝑦 =
𝜇 𝑑𝑦
𝑘𝑧 𝑑𝑝
𝑣𝑧 = [ + 𝜌𝑔]
𝜇 𝑑𝑧
Para el caso de flujo radial, despreciando los efectos gravitacionales, se tiene que la ley de
Darcy queda como:
𝑘𝑟 𝑑𝑝
𝑣𝑟 =
𝜇 𝑑𝑟
1 𝑑𝑣
𝑐=− ( )
𝑉 𝑑𝑝 𝑇
1 𝑑𝜌
𝑐= ( )
𝜌 𝑑𝑝 𝑇
𝜌=𝜌 𝑐(𝑝−𝑝𝑜 )
𝑜𝑒
Donde:
Los tipos de flujo que se generan en el yacimiento durante una prueba tienen un efecto
importante en el comportamiento de presión y generalmente se asocia una geometría de flujo
con un patrón de variación de la presión de fondo en el tiempo, sin embargo, puede existir
confusión en casos como el que se mencionó anteriormente, donde las líneas de flujo pueden
seguir varios patrones. Lo anterior se resuelve, considerando que la variación de la presión
en el pozo es afectada por la geometría de flujo de la zona que más aporta a la expansión que
genera el flujo. Esto es, si la zona que más se expande durante cierto periodo de la prueba
exhibe líneas de flujo que siguen rectas, entonces la presión en el pozo varía de acuerdo a las
ecuaciones de flujo lineal.
La zona que mayor expansión aporta se mueve a través del yacimiento y al inicio de la
producción se encuentra localizada en las vecindades del pozo, de tal manera, que se aleja y
cubre un mayor volumen a medida que transcurre el tiempo.
Geometrías de flujo.
Geometrías de flujo.
1.3 Ecuaciones de flujo y gráficas.
Ecuación de continuidad.
Si no hay pérdidas de fluido dentro de un tubo uniforme, la masa de fluido que entra es igual
a la que sale en un tubo en un tiempo. Por ejemplo, en la figura, la masa (∆m1) que entra y la
masa (∆m2) que salen en el tubo durante un tiempo corto (∆t) son:
donde A1 y A2 son las áreas transversales del tubo en la entrada y salida. Puesto que la masa
se conserva (∆m1 = ∆m2) se tiene:
𝜌1 𝐴1 𝑣1 ∆𝑡 = 𝑝1 𝐴2 𝑣2 ∆𝑡
En flujo permanente:
ṁ1 = ṁ2 = ṁ
𝜌1 ∙ 𝑄1 = 𝜌2 ∙ 𝑄2 = 𝜌 ∙ 𝑄
ṁ = 𝜌1 ∙ 𝑉1 ∙ 𝑆1 = 𝜌2 ∙ 𝑉2 ∙ 𝑆2
Para líquidos tenemos que:
(𝜌1 = 𝜌1 = 𝜌)
𝑄 = 𝑉1 ∙ 𝑆1 = 𝑉2 ∙ 𝑆2
𝐻1 = 𝐻2 + 𝐻𝑟12
𝑉11 𝑝1 𝑉22 𝑝2
+ + 𝑧1 = + + 𝐻𝑟12
2𝑔 𝛾 2𝑔 𝛾
𝑉 2⁄ = altura cinética (desnivel entre LE y LP).
2𝑔
𝑝
⁄𝛾 = altura de presión (desnivel entre LP y conducto).
𝑝
⁄𝛾 + 𝑧 = altura piezometrica.
Cuando a lo largo de un volumen de control, la velocidad del flujo varia, es por que actúan
fuerzas sobre el que lo aceleran:
⃑⃑⃑⃑⃑⃑⃑⃑⃑
∑𝐹 = ṁ ∙ 𝑎
̅̅̅̅̅
El impulso (∑ 𝐹 ∙ 𝑑𝑡) sobre la masa del volumen de control provocara una variación de su
⃑ )]:
cantidad de movimiento [𝑑(𝑚 ∙ 𝑉
̅̅̅̅̅̅̅
⃑ ) = 𝑑𝑝
∑ 𝐹 ∙ 𝑑𝑡 = 𝑑(𝑚 ∙ 𝑉
Esta variación 𝑑𝑝 del sistema es la que corresponde al instante (𝑡 + 𝑑𝑡), menos la que tenía
en t:
Flujo esférico.
Régimen de flujo que se produce cuando el patrón de flujo predominante del yacimiento es
hacia un punto. El flujo esférico se produce en terminaciones con penetración parcial y
entrada limitada. Este régimen de flujo se reconoce por una pendiente de -1/2 en la derivada
de la presión en el gráfico de diagnóstico logarítmico. Su presencia permite la determinación
de la permeabilidad esférica. Cuando el flujo esférico es seguido por el flujo radial, se pueden
cuantificar la permeabilidad horizontal y vertical.
Flujo bilineal.
El área de la sección transversal se representa con el subíndice r, para indicar que es una
función de radio de drene. De esta forma, el gradiente de presión es también una función del
radio y de forma similar se representa con un subíndice r. Para un radio r y un gasto q en
particular, la pendiente de la gráfica de presión contra el radio permanece constante todo el
tiempo. En consecuencia, durante el tiempo en que permanece constante el gasto, la
distribución de presión también permanece constante.
Esta idea puede ser extendida para aplicarse a fluidos compresibles, por ejemplo, un gas, si
el gasto se establece en unidades convenientes, por ejemplo, pies cúbicos estándar. De este
modo, el historial de presión y gasto del pozo pueden ser empleados para determinar si un
pozo se encuentra en régimen estacionario. Si el gasto es constante y la presión de fondo
permanece constante, no hay duda que el área de drene de este pozo se encuentra en régimen
de flujo estacionario.
Para que tal situación ocurra, el flujo a través del radio de drene, , deberá ser igual a el flujo
a través del radio del pozo, , y el mismo fluido deberá cruzar ambos radios. Esta condición
estrictamente no es conocida en un yacimiento. Sin embargo, un fuerte empuje por agua, por
medio del cual un gasto de agua es igual al gasto de producción proporciona un historial de
presión y gasto casi idéntico al descrito.
Las ecuaciones de flujo estacionario son además útiles en el análisis de las condiciones
cercanas al pozo. Al igual que en un sistema de flujo transitorio, el gasto cerca del pozo es
aproximadamente constante de tal forma que las condiciones alrededor del pozo son casi
constantes. Así, las ecuaciones de flujo estacionario pueden ser aplicadas a esta porción del
yacimiento sin que se presenten errores significativos. Es decir, las ecuaciones de flujo
estacionario pueden ser utilizadas para representar períodos cortos de tiempo para el flujo
alrededor del pozo.
1.3.3 Flujo pseudo-estacionario.
Después de un período inicial de producción con presión y gasto no constante, es decir, flujo
transitorio, las condiciones de frontera externa (no flujo y p = cte) comienzan a afectar la
producción en el pozo y el flujo se estabiliza. Cuando la estabilización en el yacimiento se
lleva a cabo, la condición de frontera externa de presión constante da origen al flujo
denominado como flujo pseudo-estacionario. La condición de frontera externa de presión
constante representa la frontera en la cual la presión del yacimiento se mantiene en su valor
inicial. La condición de frontera externa de presión constante es usualmente causada ya sea
por la entrada de agua de un acuífero asociado o por la inyección de agua o gas a través de
pozos inyectores, o bien, la combinación de los tres.
Para yacimientos que tengan diferente forma geométrica se puede utilizar la siguiente
ecuación:
Donde:
También se aplicará este principio para el caso de varios pozos produciendo dentro deun
mismo yacimiento. Incluso si algunos de estos produjeran a varias tasas.
Hay tener en cuenta que la aplicación del principio de superposición solo es posible para
sistemas lineales (en el sentido matemático); sin embargo, estos incluyen la mayoría de las
funciones estándares de respuestas de flujo usadas en el análisis de prueba de pozos, tales
como soluciones para el flujo radial a tasa constante, doble porosidad, pozo fracturado y
pozo con fronteras.
Note que, para la caída de presión de A en A, se incluye el factor de daño del pozo A. Esto
obedece a que estamos modelando también la caída de presión dentro de la zona de
permeabilidad cercana al pozo A. No así las caídas en A dados por los otros pozos distantes
a A. La existencia del factor de daño diferente a cero para los pozos B y C afecta solamente
a la presión dentro de su zona de permeabilidad alterada y no tiene influencia sobre la
presión en el pozo A.
Usando este método, podemos analizar simultáneamente algunos pozos fluyentes, a tasa
constante, en un yacimiento que actúa como infinito. Así, podemos modelar las llamadas
pruebas de interferencia, las cuales básicamente son diseñadas para determinar las
propiedades del yacimiento a través de la respuesta de presión obtenida en un pozo dado,
generalmente de observación (tal como el pozo A) a la producción que ocurre en uno o más
pozos (tal como los pozos B y C) dentro de un mismo yacimiento.
L L
q q
Método de imagen.
Falla
se concentra en un espacio infinitesimal de la pared del pozo, deahí el otro nombre: efecto
de piel.
Esta técnica puede ser usada también para modelar otros casos, como, por ejemplo:
- Distribución de presión para un pozo entre dos limites que se intersectan a 90º.
IMAGEN
IMAGEN
IMAGEN
Este último caso ha sido estudiado completamente por Matthews, Brons y Hazebroek y es
uno de los métodos más frecuentemente usados para estimar la presión promedio del área
de drenaje a partir de las pruebas de restauración de presión.
q
POZO A
q2
q1
q3
t1 t2
t
a)
q POZO A-1
q1
pozo 1
b)
t
c)
q POZO A-3
t2
pozo 3
(q3 – q2)
t
d)
Pozos produciendo con tasas variables.
La primera contribución a la caída de presión del pozo real es la del pozo produciendo a una
tasa q1 que empieza a fluir a un tiempo t = 0. Este pozo, en general, estará dentro de la zona
de permeabilidad alterada; así, su contribución a la caída de presión del yacimiento es:
Nótese que este primer pozo no sólo produce por tiempo t1 sino por todo el tiempo t >t2.
Empezando a un tiempo t1, la nueva tasa total real es q2. Introduciremos ficticiamenteun
pozo 2, produciendo a una tasa (q2 – q1) empezando a un tiempo t1, así que la tasa total
real después de t1 es la requerida q2 = q1 + (q2-q1). Note que el tiempo total que produce
este segundo pozo ficticio es (t – t1); note además que este pozo esta todavía dentro de la
zona de permeabilidad alterada.
Así, la contribución del pozo 2 ficticio a la caída de presión del yacimiento es (primer
cambio de tasa):
Así, la caída total para el pozo real con dos cambios en la tasa es:
Procediendo de manera similar, podemos modelar un pozo actual con docenas decambios
en su historia; podemos también modelar la historia de tasas para un pozo con tasa
continuamente variable (con una secuencia de períodos de tasa constante a una tasa
promedio durante el período) pero, muchos casos resultan en una larga ecuación, tediosa
para cálculo manual. Note, sin embargo, que tal procedimiento sólo es válido si la ecuación
4.4 ó 4.68 es válida para el tiempo total transcurrido desde que el pozo empezóa fluir en su
tasa inicial- es decir, que para el tiempo t > t2, ri debe ser menor ó igual a re de modo que
el yacimiento actúa como si fuera infinito.
Bibliografía.
https://www.sciencedirect.com/science/article/abs/pii/S037673610870092X
1?inline=1