INTEGRANTES:

BECERRA COCA ANA MELIZZA

CAIGUARA LARA MARISELA
PEREZ QUIROZ MARCOS DANIEL
PINTO FLORES JUAN ARMANDO
 CONTENIDO

1. OBJETIVO GENERAL ..............................................................................................................1

2. OBJETIVOS ESPECIFICOS ....................................................................................................1
3. INTRODUCCION .......................................................................................................................2
4. DESARROLLO ...........................................................................................................................3
4.1. ECUACIONES QUE DESCRIBEN EL FLUJO DE FLUIDOS .....................................3
4.2. ECUACION DE CONTINUIDAD ......................................................................................3
4.3. LEY DE DARCY .................................................................................................................6
4.4. ECUACIONES DE DIFUSIVIDAD .................................................................................17
4.5. ECUACION DE DIFUSIVIDAD PARA FLUJO RADIAL .............................................17
4.6. ECUACION DE DIFUSIVIDAD PARA FLUJO LINEAL ..............................................22
4.6.1. Ecuación de difusividad para flujo lineal fluido incompresible ...........................27
4.6.2. Ecuacion de difusividad para flujo lineal fluido ligeramente compresible ........30
4.6.3. Ecuacion de difusividad para flujo lineal: fluido compresible (gases) ...............33
4.7. CONDICIONES NECESARIAS PARA LA SOLUCION ..............................................40
4.8. SOLUCIÓN DE LA ECUACIÓN DE DIFUSIVIDAD RADIAL ....................................41
4.8.1. SOLUCIÓN PARA FLUJO ESTABLE ....................................................................43
4.8.2. SOLUCIÓN PARA FLUJO PSEUDO ESTABLE ..................................................43
5. CONCLUSION Y RECOMENDACIÓN .................................................................................45
6. BIBLIOGRAFIA ........................................................................................................................45
1. OBJETIVO GENERAL

Enunciar las ecuaciones que describen el flujo de fluidos en los yacimientos para
una correcta evaluación del yacimiento.

2. OBJETIVOS ESPECIFICOS

• Describir la ecuación de continuidad

• Describir la ecuación de la ley de Darcy para flujo lineal y radial
• Desarrollar la ecuación de difusividad mediante la aplicación de la ley de
continuidad y la ley de Darcy.
• Desarrollar las ecuacion de difusividad para flujo lineal y radial para los
diferentes tipos de fluidos.

1
3. INTRODUCCION

La evaluación de yacimientos tiene el objetivo de presentar los fundamentos teóricos

y los principios físicos sobre los cuales se basa el análisis de las pruebas de presión,
proporciona las bases teóricas que permiten el entendimiento de las relaciones
matemáticas que se plantean bajo diferentes condiciones de flujo de fluidos en un
medio poroso y permeable.

El área de pruebas de presión requiere del conocimiento y comprensión de términos

específicos para poder interpretar los fenómenos que se analizan; y así, poder dar un
correcto uso de la información que se obtiene de una prueba de presión. Las pruebas
de presión permiten determinar cualitativamente muchas características de un
yacimiento; por ejemplo: estimar el petróleo original in situ y el petróleo que puede
ser extraído, lo cual refleja el estado de agotamiento del yacimiento.

Una explotación racional de los hidrocarburos, asegura una recuperación aceptable

de los mismos, la optimización de la producción de un yacimiento depende en gran
parte del conocimiento de los parámetros físicos del sistema, los cuales se obtienen
con el análisis de datos de pruebas de presión.

El estudio de flujo de fluidos en los yacimientos está basado en dos conceptos

fundamentales: ley de Darcy y el Balance de materia y sus variaciones dependen de
las características del yacimiento, entre las cuales deben considerarse: la
configuración geométrica del sistema, la compresibilidad de los fluidos, la
invariabilidad de las tasas de flujo y de la presión con el tiempo, o bien el flujo de una
o más fases simultáneas.

El principal objetivo de este trabajo es el de presentar las diferentes relaciones

matemáticas que han sido diseñadas para describir el comportamiento de flujo de
fluidos del yacimiento y señalar como el ingeniero puede reconocer las características
de cada uno de ellos.

2
 4. DESARROLLO

4.1. ECUACIONES QUE DESCRIBEN EL FLUJO DE FLUIDOS

Todas las ecuaciones prácticas que describen el flujo de fluidos en medios porosos se
basan en dos conceptos básicos: la ley de Darcy y la ley de conservación de masa.
Los conceptos más simples de ingeniería de yacimientos se basan en la primera de estas
leyes, pero los complejos y generalmente los más útiles se basan en ambas.

4.2. ECUACION DE CONTINUIDAD

La ecuación de continuidad es una expresión matemática del principio de conservación

de masa. En los problemas de transferencia (fluidos, calor, etc.), uno de los postulados
más útiles es el Principio de Conservación, por medio del cual podemos establecer
una continuidad en alguno de los parámetros dentro de un yacimiento.

En el caso de flujo de fluidos en el medio poroso, el parámetro a controlar es la masa. El

principio de conservación de masa establece que:

En el proceso general de la derivación de la ecuación se selecciona un volumen

representativo, de acuerdo a la geometría de flujo del problema en interés.

Considérese el elemento de flujo que se muestra en la siguiente figura tiene un ancho dr

y está localizado a una distancia r del centro del pozo.

El elemento poroso tiene un volumen diferencial dV. Aplicando el concepto de balance

de materiales, la tasa másica de flujo que entra al elemento menos la tasa másica de
flujo que sale de el durante un tiempo diferencial Δt, debe ser igual a la acumulación
durante este intervalo de tiempo, o sea:

3
 Ecuación 1

A continuación se expresan individualmente los términos de la ecuación 55.

• Primer término: Masa que entra al elemento de volumen durante un intervalo de

tiempo Δt
Ecuación 2

donde v es la velocidad del fluido fluyendo en pie/dia; ρ la densidad del fluido a (r + dr)
en lb/pie3; A el área a r+dr en pie2; y ∆t el intervalo de tiempo en días.

El área del elemento de entrada es:

Ecuación 3

Combinando la ecuación 2 con la ecuación 3 resulta:

Ecuación 4

• Segundo término: Masa que sale del elemento de volumen durante un intervalo de
tiempo Δt.

Adoptando el mismo criterio para esta masa se tiene:

4
 Ecuación 5

• Tercer término: Masa total acumulada en el elemento de volumen durante un

intervalo de tiempo Δt.

El volumen total de un elemento de radio r viene dado por: V = π r 2 h

dV
Diferenciando la ecuacion anterior con respecto a r, resulta: = 2πrh
dr
O también:

Ecuación 6

Luego la masa total acumulada en el tiempo

Sustituyendo por dV:

Masa total acumulada=

Reemplazando los términos calculados en la ecuación 1 se obtiene:

Dividiendo la ecuación anterior por (2𝜋𝑟ℎ)𝑑𝑟 y simplificando, resulta:

O también:

Ecuación 7

Donde Ø es la porosidad, ρ la densidad en lb/pie3, v la velocidad del fluido en pie/dia.

La ecuación 7 se denomina ecuación de continuidad y muestra el principio de

conservación de masa en coordenadas radiales.

5
 4.3. LEY DE DARCY

La ecuación de Darcy es la más aplicada para el problema de flujo de fluidos en medios

porosos.

Esta ecuación es el resultado de estudios experimentales de flujo de agua a través de

filtros de arena no consolidada, Henry Darcy dedujo la fórmula que lleva su nombre. La
ley se ha extendido, con ciertas limitaciones, al movimiento de otros fluidos no miscibles,
en rocas consolidadas y otros medios porosos. La ley de Darcy enuncia que la velocidad
de un fluido homogéneo en un medio poroso es proporcional al gradiente de presión e
inversamente proporcional a la viscosidad del fluido.

Ecuación 8

𝑣 es la velocidad aparente en centímetros por segundo e igual a q/A donde q es la rata

volumétrica del flujo, en centímetros cúbicos por segundo y A es el área de la sección
transversal total aparente de la roca, en centímetros cuadrados. En otras palabras, A
incluye tanto el área del material solido (esqueleto mineral) de la roca como también el
área de los canales porosos. La viscosidad del fluido µ se expresa en centipoises y el
gradiente de presión dp/ds tomado en la misma dirección que q y 𝑣 en atmósfera por
centímetro. La constante de proporcionalidad, k es la permeabilidad puede denominarse
como la conductividad del fluido y es similar al termino de conductividad eléctrica. El
signo negativo índica que, si se toma el flujo positivo en la dirección positiva de s, la
presión disminuye en esa dirección y la pendiente dp/ds es negativa.

La ley de Darcy se aplica solo en la región de flujo laminar. En flujo turbulento que ocurre
a altas velocidades el gradiente de presión aumenta a una rata mayor que la del flujo.
Afortunadamente excepto en casos de muy altas ratas de inyección o de producción en
la vecindad del pozo, el flujo en el yacimiento y en la mayoría de las pruebas hechas en
el laboratorio es laminar y la ley de Darcy se cumple.

La ley de Darcy no se aplica a flujos en canales porosos individuales sino a partes de la

roca de dimensiones razonablemente grandes comparadas con el tamaño de los canales
porosos; en otras palabras, es una ley estadista que promedia el comportamiento de
muchos canales porosos. Por esta razón con muestras de formación de areniscas
uniformes de uno o dos centímetros se obtiene resultados satisfactorios de mediciones
de permeabilidad es cambio rocas de tipo fracturado o drusoso exigen muestras mucho
más grandes.

Debido a la porosidad de la roca a la tortuosidad de las líneas de flujo y a la ausencia de

flujo en algunos de los espacios porosos (incomunicados) la velocidad real del fluido varia

6
de lugar a lugar dentro de la roca, y mantiene un promedio mucho más alto que la
velocidad aparente. Como las velocidades reales no son medibles por lo general, y para
mantener la porosidad y permeabilidad independientes, las velocidades aparentes
constituyen la base de la ley de Darcy. Es decir, la velocidad real promedia de avance
de un fluido es la velocidad aparente dividida por la porosidad, cuando el fluido satura
por completo la roca.

La unidad de permeabilidad es el Darcy. Se dice que una roca tiene la permeabilidad de

una Darcy cuando un fluido con una viscosidad de un centipoise avanza a una velocidad
de un centímetro por segundo bajo un gradiente de presión de una atmosfera por
centímetro. Como es una unidad bastante alta para la mayoría de las rocas productoras
la permeabilidad generalmente se expresa en milésimas de Darcy es decir milidarcys
0,001 darcy. Las permeabilidades de formaciones de gas y petróleo comercialmente
productoras varias desde pocos milidarcys a varios de miles. Las permeabilidades de
calizas intergranulares pueden ser solo una fracción de un milidarcy y aun tener
producción comercial, siempre y cuando la roca contenga fracturas u otro tipo de
aberturas adicionales naturales o artificiales. Rocas con fracturas y tipo drusoso pueden
tener permeabilidades muy altas y algunas calizas cavernosas se aproximan al
equivalente de tanques subterráneos.

El gradiente dp/ds es la fuerza de empuje y se debe a los gradientes de presión del fluido
y en parte o totalmente a los gradientes hidráulicos (gravitacionales), superpuestos e
intercambiables. En muchos casos de interés practico los gradientes hidráulicos, aunque
siempre presentes son pequeños comparados con los gradientes de presión del fluido y
son frecuentemente despreciables. En otros casos los gradientes hidráulicos son
importantes y deben considerarse en particular en operaciones de producción por
bombeo en yacimientos con presiones agotadas o en yacimientos con capa de gas en
expansión con buenas características de drenaje gravitacional.

En caso de que ambos actúen simultáneamente la ecuación puede modificarse:

𝑘 𝑑𝑝
𝑣 = − µ ( 𝑑𝑠 − 9,67𝑥10−4 𝜌 𝑐𝑜𝑠 𝜃) Ecuación 9

𝜌 es la densidad del fluido en gramos por centímetro cubico y 𝜃 es el ángulo entre la

dirección positiva de s y la línea vertical en dirección hacia abajo 9,67 𝑥10−4 convierte el
gradiente de presión en gramos por centímetro cuadrado por centímetro a atmosfera por
centímetro.

En unidades de campo seria:

𝑘 𝑑𝑝
𝑣 = −1,127 ∗ ( 𝑑𝑠 − 0,433 𝜌 𝑐𝑜𝑠𝜃) Ecuación 10
µ

7
0,433 𝜌 𝑐𝑜𝑠𝜃 es el gradiente hidráulico en lpc por pie, 𝜌 la gravedad especifica relativa al
agua.

Los gradientes hidráulicos en los yacimientos varias de un máximo alrededor de 0,5000

lpc/p para salmueras a 0,4333 lpc/p par agua dulce a 60°F y varían de acuerdo con la
presión, temperatura y salinidad del agua.

Los gradientes de petróleos de gas a alta presión y de condensado de gas varían entre
0,10 y 0,30 lpc/p y varían de acuerdo con la presión, temperatura y composición del
fluido.

Los gases a presiones bajas tienen gradientes muy bajos, alrededor de 0,002 lpc/p para
gas natural a 100 lpca.

La ley de Darcy es aplicable solo si se cumple las siguientes condiciones:

1. Flujo laminar: si el número de Reynolds es mayor que uno, no es aplicable.

2. Flujo continuo
3. Fluidos incompresibles
4. Formación homogénea: no existe reacción entre el fluido y la roca.

Para flujo turbulento, que se produce a altas velocidades el gradiente de presión aumenta
en mayor porcentaje que la tasa de flujo y por lo tanto la aplicación de la ley de darcy
puede dar lugar a errores graves, haciéndose necesario realizar una modificación en la
ecuación de Darcy.

4.3.1. Ley de Darcy para flujo radial

Para un sistema radial el gradiente de presión es positivo como se muestra en la figura

y la ecuación de darcy puede expresarse en la siguiente forma generalizada:

8
 Ecuación 11

Donde q es la tasa de flujo volumétrica en el radio r; Ar el área de la sección transversal

de flujo en el radio r; (dp/dt)r el gradiente de presión al radio r y v la velocidad aparente
del fluido al radio r.

El área de la sección transversal al flujo a un radio r es el área de un cilindro y para un

pozo totalmente penetrado con un espesor neto h viene dada por: Ar=2πrh

4.3.1.1. Flujo Radial de fluidos incompresibles

En un sistema radial, los fluidos se mueven hacia el pozo productor en todas las
direcciones. Antes de que el flujo tenga lugar, existe una diferencia de presión y, por lo
tanto, si un pozo es productor de petróleo, el flujo de los fluidos ocurre de la formación
hacia la boca del pozo, donde la presión debe ser menor que la existente en un punto de
la formación situado a cierta distancia del pozo.

La presión de la formación en el hoyo del pozo productor se conoce como presión de

flujo de fondo, Pwf.

Considérese un pozo tal como se muestra en la

Figura, situado en una formación cilíndrica horizontal
de radio exterior re y espesor h. El pozo tiene un radio
rw y las presiones en el pozo y en el radio exterior
son Pwf y Pe , respectivamente.

Para la aplicación de la ecuación que rige este flujo

se parte de la siguiente ecuación

Ecuación 12

Donde V es la velocidad aparente del fluido en BY

/día-𝑝𝑖𝑒 2 .

Haciendo en esta ecuación 𝑥 = 𝑟 𝑦 𝐴 = 2𝜋𝑟ℎ, y tomando en cuenta que, para un sistema

radial, r aumenta en la misma dirección que la presión p, lo que implica que 𝛿𝑝/𝛿𝑟 es
positivo y, por lo tanto, no se requiere el signo menos, resulta:

𝑞 𝑞 𝑘 𝑑𝑝
𝑉= = = 0,001127
𝐴𝑟 2𝜋𝑟ℎ 𝜇 𝑑𝑟

9
La tasa de flujo para un sistema de petróleo crudo se expresa generalmente en unidades
de superficie, esto es, en BN. Usando el símbolo Qo para representar el flujo de petróleo
en BN /día, tenemos:

𝑞 = 𝐵𝑜×𝑄𝑜

Donde Bo es el factor volumétrico del petróleo en la formación en BY /BN. Luego la tasa

de flujo en la ecuación: de Darcy se puede expresar de la siguiente manera:

𝐵𝑜𝑄𝑜 𝑘 𝑑𝑝
= 0,001127
2𝜋𝑟ℎ 𝜇 𝑑𝑟

Separando variables y considerando k, h, 𝜇𝑜 constantes, se puede integrar la ecuación

anterior entre los radios r1 y r2 cuando las presiones son p1 y p2, que resulta en:

Ecuación 13

Para un sistema incompresible en una formación homogénea, la ecuación 13 se

simplifica para dar:
r2 P2
Qo dr 0,001127K
∫ = ∫ dp
2πh r μoBo
r1 P1

Resolviendo la integral:

0,00708Kh(P2 − P1)
Qo =
μoBo ln(r2⁄r1)

Frecuentemente, los dos radios de interés son el radio del pozo rw y el radio de drenaje
re. Así se tiene:

Ecuación 14

En este caso Qo es la tasa de flujo de petróleo en BY /día; Pc, la presión en el radio de

drenaje en lpc; Pwf, la presión de fondo fluyente en lpc; k, la permeabilidad absoluta en
md; μo, la viscosidad del petróleo en cp; h, el espesor en pies; Bo, el factor volumétrico
del petróleo en la formación en BY /BN; rw y re son los radios del pozo y de drenaje en
pies, respectivamente.

10
El radio de drenaje o radio exterior re, se determina usualmente a partir del área de
espaciamiento entre los pozos. Si se considera, por ejemplo, que esta es un circulo, se
tiene:

πre2 = 43560A, o sea:

Ecuación 15

Donde A es el espaciamiento entre pozos en acres.

En la práctica, ni el radio externo ni el pozo se conocen con precisión, afortunadamente,

estos términos entran en la ecuación como logaritmos, de manera que el error en ella
será menor que los errores en los radios.

La ecuación 14 se puede arreglar para calcular la presión p a cualquier radio r; así, se

tiene:

Ecuación 16

La presión en el radio de drenaje, Pe, usada en la ecuación 14 no puede ser medida

fácilmente, pero en caso de estar presente un acuífero fuerte y activo, se mantiene
cercana a la presión inicial del yacimiento. Craft y Hawkins demostraron que la presión
promedio del yacimiento está localizada a una distancia cercana a 61% del radio de
drenaje, si se mantienen condiciones de flujo continuo. Así, en la ecuación 16 se sustituye
r por 0,61re se obtiene:

𝑄𝑜𝐵𝑜𝜇𝑜 0,61𝑟𝑒
𝑝(𝑟 = 0,61𝑟𝑒 ) = ̅̅̅
𝑝𝑟 = 𝑝𝑤𝑓 + [ ] ln ( )
0,00708𝐾ℎ 𝑟𝑤

O en términos de la tasa de flujo:

Ecuación 17

0,61𝑟𝑒 𝑟
Pero, como ln ( ) = ln (𝑟𝑒 ) − 0,5; entonces:
𝑟𝑤 𝑤

11
 Ecuación 18

4.3.2. Flujo radial de fluidos ligeramente compresibles

La ecuación 19 describe la relación que existe entre la presión el volumen para los fluidos
ligeramente compresibles, esto es:
Ecuación 19

Esta ecuación puede modificarse para expresarla en términos de la tasa de flujo, así:

Ecuación 20

Craft et al. Utilizaron la ecuación 20 para expresar la dependencia de la tasa de flujo con
presión para fluidos ligeramente compresibles. Si esta ecuación se sustituye en la forma
radial de la ley de Darcy, entonces:

q q ref [1 + c(pref − p)] k dp

= = 0,001127
Ar 2πrh μ dr

Donde q ref es la tasa de flujo a una presión de referencia pref .

Separando variables e integrando, se obtiene la siguiente expresión:

re pe
q ref μ dr dp
∫ = 0,001127 ∫
2πrh r 1 + c(pref − p)
rw pwf

Despejando q ref , resulta:

0,00708Kh 1 + c(pe − pref )

q ref = [ re ] ln [ ]
μoBoCo ln( ⁄rw) 1 + c(pwf − pref )

Donde q ref es la tasa de flujo a la presión de referencia pref . Seleccionando la presión de

fondo fluyente pwf, como la presión de referencia y expresando la tasa de flujo en BN
/día, se obtiene finalmente:

12
 Ecuación 21

Donde Qo , es la tasa de flujo de petróleo en BN /día, co , el coeficiente isotérmico de

compresibilidad del petróleo en lpc −1y k, la permeabilidad absoluta en md.

4.3.3. Flujo radial de gases compresibles

La ecuación básica de la ley de Darcy en su forma diferencial para flujo horizontal laminar
es válida para describir el flujo de ambos sistemas líquidos y gases. Para flujo radial de
un gas, la forma de la ecuación de Darcy es:

Ecuación 22

Donde q g es una tasa de flujo de gas a un radio r en BY /día; r, la distancia radial en pies;
h, el espesor de la zona en pies; μg , la viscosidad del gas en cp, p, la presión en lpc;
0,001127, una constante de conversión para llevar de unidades Darcy a unidades de
campo.

La tasa de flujo de gas se expresa generalmente en PCN/día. Designando la tasa de flujo

de gas en condiciones normales como Qg , la tasa de flujo q g a una determinada presión
y temperatura se puede convertir a las condiciones normales aplicando la ecuación de
estado para un gas real para ambas condiciones, Entonces:

5,615q g p Qg psc
=
zRT zsc RTsc

Despejando q g resulta:

Ecuación 23

Donde psc y Tsc son la presión en lpca y la temperatura en °R en condiciones normales,

respectivamente; Qg , la tasa de flujo de gas en PCN/día, q g , la tasa de flujo de gas en el
radio r en BY /día; r, la distancia radial en pies; p, la presión a un radio r en lpca, T, la
temperatura del yacimiento en °R; z, el factor de compresibilidad del gas a p y T; y zsc ,
el factor de compresibilidad del gas en condiciones normales ≅ 1.

13
Combinando las ecuaciones 22 y 23, se obtiene:

zT psc 0,001127(2πrh)k dp
Qg ( ) ( )=
p 5,615Tsc μg dr

Suponiendo que T = 520 °R y psc = 14.7 lpca

Ecuación 24

Integrando la ecuación 24 desde las condiciones en el pozo (rw, pwf) hasta cualquier
punto en el yacimiento (r, p), resulta:

Ecuación 25

Si la ecuación 25 se le imponen las condiciones de la ley de Darcy, esto es:

• Flujo continuo, lo cual requiere que Qg sea constante para cualquier radio, y
• Formación homogénea, que implica que k y h son constantes, se tiene:

Combinando las dos expresiones anteriores, resulta:

Ecuación 26

14
 Ecuación 27

La ecuación de seudopresión fue introducida por primera vez por Al-Hussainy et al. Con
el fin de linearizar la ecuacion que describe el flujo de un gas real.

La ecuación 26 escrita en términos del potencial de un gas real es:

O bien:

Ecuación 28

La ecuación 28 indica que un gráfico de da una línea recta de pendiente

e intercepto

Como se observa en la figura siguiente

La tasa de flujo viene dado por:

15
 Ecuación 29

En el caso particular en que r=re, entonces:

Ecuación 30

Donde: Qg es la tasa de flujo de gas en PCN/dia

el potencial de una gas real evaluado desde 0 hasta Pe en lpc 2 /cp el

potencial de una gas real evaluado desde 0 hasta Pwf en lpc 2 /cp

k la permeabilidad absoluta en md

re el radio de drenaje en pies

rw el radio del pozo en pies
h el espesor en pies

Como la tasa de flujo de gas se expresa comúnmente en MPCN/dia se tiene:

Ecuación 31

Donde Qg es la tasa de flujo de gas en MPCN/dia.

La ecuación 31 puede expresarse en términos de la presión promedio del yacimiento 𝑃̅𝑟

en lugar de la presión inicial del yacimiento Pe:

Ecuación 32

Usando las expresiones de la ley de Darcy, la ecuacion 25 y 32 pueden aproximarse

2
sacando el término (µ𝑔 𝑍) fuera de la integral como una constante. Es importante señalar

16
que (µg Z) se considera constante solo para presiones menores a los 2000 lpc. Luego,
la ecuación 31 puede reescribirse así:

Sacando el término e integrando, resulta:

Ecuación 33

Donde Qq es la tasa de flujo de gas en MPCN/dia y k la permeabilidad en md.

El término (μ̅ 𝑔 𝑧̅) se evalua a una presión promedio 𝑃̅ definida por:

Pwf 2 − Pc 2
̅=√
P
2

El método de aproximación anterior se conoce como método de las presiones al

cuadrado y esta limitado a cálculos de flujo cuando la presión del yacimiento es menor
a 2000 lpc.

4.4. ECUACIONES DE DIFUSIVIDAD

Es la combinación de las principales ecuaciones que describen el proceso físico del

movimiento de fluido dentro del reservorio, combina la ecuación de continuidad (que es
el principio de la conservación de la masa, y de aquí obtenemos el balance de materia),
la ecuación de flujo (ecuación de Darcy) y la ecuación de estado (compresibilidad).

4.5. ECUACION DE DIFUSIVIDAD PARA FLUJO RADIAL

Para relacionar la velocidad de flujo con el gradiente de presión dentro de una volumen
DV, debe introducirse la ecuación de transporte en la ecuación de continuidad. La
ecuación de Darcy es esencialmente la ecuación básica de movimiento, la cual establece
𝜕𝑝
que la velocidad v es proporcional al gradiente de presión (𝜕𝑟 ).

De la ecuación 12:

17
 Ecuación 34

Donde k es la permeabilidad en md y v la velocidad del fluido en pie/dia.

Combinando la ecuacion 34 con la ecuación 7 resulta:

Ecuación 35

Expandiendo la parte derecha de la ecuación 35, se obtiene:

Ecuación 36

Como la porosidad se relaciona con la compresibilidad de la formación por medio de la

ecuación:

Ecuación 37

𝜕Ø
Al aplicar la regla de la cadena de la diferenciación a resulta:
𝜕𝑡

Sustituyendo la ecuación 37 en la expresión anterior, queda:

Finalmente, Sustituyendo esta expresión en la ecuación 36 y el resultado en la ecuación

35, se obtiene:

Ecuación 38

La ecuación 38 es la ecuación general en derivadas parciales que se utiliza para describir

el flujo de cualquier fluido que fluye en dirección radial en el medio poroso.

Adicionalmente a las suposiciones iniciales, en la deducción se tomó en cuenta la

ecuación de Darcy, lo cual implica que el flujo es laminar. Es decir, la ecuación no está

18
restringida para ningún tipo de fluido y es igualmente válida tanto para liquidos como
para gases. Sin embargo, los fluidos compresibles y los ligeramente compresibles deben
tratarse separadamente para poder desarrollar las ecuaciones prácticas aplicables a
cada uno de ellos. A continuación se estudiaran cada uno de estos casos.

4.5.1. Flujo radial de fluidos ligeramente compresibles

Para simplificar la ecuación 36 supongamos que la permeabilidad y la viscosidad son

constantes en el rango de presión, tiempo y distancia. Esto lleva a:

Ecuación 39

Expandiendo la ecuación anterior:

Usando la regla de la cadena en las expresiones anteriores:

Dividiendo la expresión anterior entre la densidad del fluido ρ se obtiene:

Ahora, bien como la compresibilidad de cualquier fluido viene dada por la siguiente
1 𝜕𝜌
ecuacion: 𝑐 = 𝜌 𝜕𝑝 al sustituirla en la expresión anterior queda:

𝜕𝑝
El término 𝑐(𝜕𝑟 )2 se considera muy pequeño y puede ser ignorado, así se tiene:

Ecuación 40

19
Si la compresibilidad total está definida por:

Ecuación 41

Al combinar la ecuación 41 con la 40 y Rearreglando, resulta finalmente:

Ecuación 42

Donde el tiempo t esta presado en días.

La ecuación 42 se denomina ecuacion de difusividad y es particularmente usada en el

análisis de datos de pruebas de pozos donde el tiempo t se expresa comúnmente en
horas. Rearreglando la ecuación resulta:

Ecuación 43

Donde k es la permeabilidad en md; r la distancia radial en pies; p la presión en lpca;

𝑐𝑡 la compresibilidad total en lpc-1; t el tiempo en horas, ∅ la porosidad en fracción; y µ la
viscosidad en cp.

Cuando el yacimiento contiene más de un flujo, la comprensibilidad total puede

calcularse por:

Ecuación 44

Donde Co, Cw, Cg son las comprensibilidades del petróleo, agua y gas, respectivamente,
mientras que So, Sw y Sg son saturaciones de estos fluidos, nótese que la introducción
de 𝑐𝑡 en la ecuación 42, no hace que la misma sea aplicable a flujo multifásico, el uso
de 𝑐𝑡 , definida por la ecuación 44 simplemente toma en cuenta la compresibilidad del
cualquier fluido inmóvil que puede estar en el yacimiento con el fluido que está fluyendo.

0,000264𝑘
El término de la ecuación 43, se denomina constante de difusividad y se
∅𝜇𝑐𝑡
denota por el símbolo η, es decir:

Ecuación 45

20
Así la ecuación de difusividad puede escribirse en una forma más conveniente de la
siguiente manera:

Ecuación 46

4.5.2. Flujo radial de fluidos compresibles

La viscosidad y densidad de un gas varían significativamente con la presión y por lo tanto

las suposiciones de la ecuación 46 no satisfacen los sistemas de gas, esto es, de los
fluidos compresibles. Para poder desarrollar una función matemática que describa el flujo
de fluidos compresibles en el yacimiento, deben considerarse las siguientes ecuaciones
para gases reales:

• Ecuación de densidad

• Ecuación de compresibilidad del gas

Combinando estas dos ecuaciones básicas con la ecuación 38 resulta:

Ecuación 47

Donde t es el tiempo en horas; k la permeabilidad en md; Ct la compresiblidad isotérmica

en 𝑙𝑝𝑐 −1; y Ø la porosidad.
Al-Hussainy, Ramey y Crawford linealizaron la ecuación anterior introduciendo la
definición de potencial para un gas real m(p). Así se tiene:

Ecuación 48

Diferenciando la ecuación anterior con respecto a p, resulta:

21
 Ecuación 49

De donde se obtienen las siguientes relaciones aplicando la regla de la cadena:

Ecuación 50

Ecuación 51

Sustituyendo la ecuación 44 en las ecuaciones 49 y 50 resulta:

Ecuación 52

Ecuación 53

Combinando las ecuaciones 52 y 53 con la ecuación 47 se obtiene, finalmente:

Ecuación 54

La ecuación 54 es la ecuación de difusividad para fluidos compresibles, la cual relaciona

el potencial de un gas real con el tiempo t y el radio r.

4.6. ECUACION DE DIFUSIVIDAD PARA FLUJO LINEAL

La ecuación de difusividad es una ecuación de movimiento que permite analizar el
comportamiento de presión por posición y tiempo en un medio poroso y como se dijo
anteriormente se requiere de 3 ecuaciones: ecuación de continuidad, ecuación de flujo y
la ecuacion de estado.

Para expresar la ecuación de continuidad debe suponerse un medio poroso como se

muestra en la figura; en este medio se presenta un flujo en las direcciones x, y, z.
22
 Fig. Medio poroso con flujo lineal

La masa que está entrando al medio en un intervalo de tiempo Δt es:

𝐮𝐱 𝛒(∆𝐲∆𝐳)∆𝐭 + 𝐮𝐲 𝛒(∆𝐱∆𝐳)∆𝐭 + 𝐮𝐳 𝛒(∆𝐲 ∆𝐱)∆𝐭

Y la masa que sale del mismo medio es.

[𝐮𝐱 𝛒 + ∆(𝐮𝐱 𝛒)](∆𝐲 ∆𝐳)∆𝐭 + (𝐮𝐲 𝛒 + ∆(𝐮𝐲 𝛒)) (∆𝐱 ∆𝐳)∆𝐭 + [𝐮𝐳 𝛒 + ∆(𝐮𝐳 𝛒)](∆𝐲 ∆𝐱)∆𝐭

Durante el intervalo Δt, la masa que se acumula en el medio Δx Δy Δz es:

[masa acumulada] = [masa que entra] – [masa que sale]

𝐮𝐱 𝛒 (∆𝐳 ∆𝐲)∆𝐭 − [𝐮𝐱 𝛒 + ∆(𝐮𝐱 𝛒)] (∆𝐲 ∆𝐳)∆𝐭 + 𝐮𝐲 𝛒(∆𝐱 ∆𝐳)∆𝐭 − [𝐮𝐲 𝛒 + ∆(𝐮𝐱 𝛒)](∆𝐱∆𝐳)∆𝐭 +
𝐮𝐳 𝛒(∆𝐱 ∆𝐲)∆𝐭 − [𝛍𝐳 𝛒 + ∆(𝐮𝐳 𝛒)](∆𝐱 ∆𝐲) ∆𝐭 = −[∆𝐮𝐳 𝛒(∆𝐱 ∆𝐲) + ∆𝐮𝐲 𝛒(∆𝐱 ∆𝐳) +
∆𝐮𝐱 𝛒(∆𝐲 ∆𝐳)]∆𝐭

∆(𝒖𝒙 𝝆) ∆(𝒖𝒚 𝝆) ∆(𝒖𝒛 𝝆)

= ∆𝒙∆𝒚∆𝐳 ( + + ) ∆𝒕 Ecuación 55
∆𝒙 ∆𝒚 ∆𝒛

A la expresión anterior hay que agregarle la masa que entra o sale del medio a través de
fuentes o sumideros (el caso de un pozo productor o inyector).

Así que la expresión para la acumulación de masa, teniendo en cuenta lo anterior, se

convertirá en.

∆(𝐮𝐱 𝛒) ∆(𝐮𝐲 𝛒) ∆(𝐮𝐳 𝛒)

= ∆𝐱∆𝐲∆𝐳 ( + + ) ∆𝐭 − 𝐪𝐦 ∆𝐱∆𝐲∆𝐳∆𝐭 Ecuación 56
∆𝐱 ∆𝐲 ∆𝐳

23
Además, la acumulación de masa Δx Δy Δz puede obtenerse también de

((∆𝒙∆𝒚∆𝒛)𝚽𝝆)𝒕+∆𝒕 − (∆𝒙∆𝒚∆𝒛 𝚽𝛒)𝒕 Ecuación 57

E igualando las ecuaciones (56) y (57) y dividiendo por ΔxΔyΔzΔt se tiene.

∆(𝐮𝐱 𝛒) ∆(𝐮𝐲 𝛒) ∆(𝐮𝐳 𝛒) (𝚽𝛒)𝐭+∆𝐭 −(𝚽𝛒)𝐭

−( + + ) − 𝐪𝐦 = Ecuación 58
∆𝐱 ∆𝐲 ∆𝐳 ∆𝐭

Y cuando Δx Δy Δz y Δt tienden a cero, se tiene:

𝛛(𝐮 𝛒) 𝛛(𝐮 𝐲 𝛒) 𝛛(𝐮 𝛒) 𝛛(𝚽𝛒)
− 𝐱 − − 𝐳 − 𝐪𝐦 = Ecuación 59
𝛛𝐱 𝛛𝐲 𝛛𝐳 𝛛𝐭

Si se tiene en cuenta que la velocidad del fluido es una cantidad vectorial esta puede
expresarse como.
Ecuación 60

Donde ux, uy y uz son las componentes del vector de la velocidad del fluido en las
direcciones x, y, z.

Usando la ecuación (59), la ecuación (60) puede escribirse como

Ecuación 61

Si se introduce ahora la ecuación de Darcy para flujo lineal en la ecuación (59), se

tiene.

𝛛 𝐤 𝐱 𝛛𝐏 𝛛 𝐤 𝐲 𝛛𝐏 𝛛 𝐤 𝐳 𝛛𝐏 𝛛(𝛒𝚽)
( 𝛒) + ( 𝛒) + ( 𝛒) − 𝐪𝐦 = Ecuación 62
𝛛𝐱 𝛍 𝛛𝐱 𝛛𝐲 𝛍 𝛛𝐲 𝛛𝐳 𝛍 𝛛𝐳 𝛛𝐭

La ecuación (62) también puede escribirse como.

Ecuación 63

Si el flujo va solo en una dirección, los componentes en las otras direcciones no

aparecerán en la ecuación (62). En este caso la ecuación de continuidad se plantea de
la siguiente forma, suponiendo el volumen del medio que se muestra en la figura:

24
 Fig. Volumen para flujo lineal en una dimensión

Dividiendo ahora por Δx Δt, la ecuación anterior se convierte en.

∆(𝛒𝐮𝐱 𝐀(𝐗))
− − 𝐪𝐦 𝐀(𝐗) = [(𝐀(𝐗)𝛒𝚽)𝐭+𝚫𝐭 − (𝐀(𝐗)𝛒𝚽)𝐭 ]/𝚫𝐭 Ecuación 64
∆𝐱

Y expresando la ecuación (64) en forma diferencial se convierte en.

𝛅 𝛅
−
𝛅𝐱
(𝛒𝐮𝐱 𝐀(𝐗)) − 𝐪𝐦 𝐀(𝐗) = 𝛅𝐭 (𝐀(𝐗)𝛒𝚽) Ecuación 65

Si se introduce la ecuación de Darcy en la ecuación (65), suponiendo que se pueda

aplicar la ecuación (8), se tiene

𝛅 𝐤 𝛅𝐏 𝛅
𝛅𝐱
(𝛒 𝛍 𝐀(𝐗) 𝛅𝐱 ) − 𝐪𝐦 𝐀(𝐗) = 𝛅𝐭 (𝐀(𝐗)𝛒𝚽) Ecuación 65

La ecuación (65) puede acomodarse a la ecuación (63), aunque en este caso el tensor
de permeabilidades solo tendrá la componente en la dirección x, k x. La ecuación (63)
podrá escribirse en forma general como.

Ecuación 66

25
Donde α=1 cuando el flujo es en tres dimensiones y α=A(x) cuando el flujo es en una
sola dimensión.

Cuando el flujo es en dos dimensiones x, y, el volumen elemental se representa como

en la figura, y la ecuación de continuidad para cada una de las direcciones de flujo se
plantea así: en la dirección x

Fig. Volumen para flujo lineal en dos dimensiones

Si la ecuación anterior se divide por ΔxΔyΔt, se tiene.

∆(𝛒𝐮𝐱 𝐇(𝐗)) ∆(𝛒𝐮𝐲 𝐇(𝐘)) (𝐇(𝐱,𝐲)𝛒𝚽)𝐭+𝚫𝐭 −(𝐇(𝐱,𝐲)𝛒𝚽)𝐭

− − − 𝐪𝐦 𝐇(𝐱, 𝐲) =
∆𝐱 ∆𝐲 𝚫𝐭

26
Y si la ecuación anterior se expresa en forma diferencial y, además, se tiene en cuenta
que H no es función del tiempo, se tiene.

𝛅(𝛒𝐮𝐱 𝐇(𝐱)) 𝛅(𝛒𝐮𝐲 𝐇(𝐲)) (𝛒𝚽)

− − − 𝐪𝐦 𝐇(𝐱, 𝐲) Ecuación 67
𝛅𝐱 𝛅𝐲 𝛅𝐭

Si se introduce la ecuación de Darcy, la ecuación (8), en la ecuación 67 se tiene.

𝐤 𝛅𝐏 𝐤𝐲 𝛅𝐏
𝛅(𝛒 𝐱 𝐇(𝐱) ) 𝛅(𝛒 𝐇(𝐲) ) 𝛅(𝛒𝚽)
𝛍 𝛅𝐱
+
𝛍 𝛅𝐲
− 𝐪𝐦 𝐇(𝐱, 𝐲) = 𝐇(𝐱, 𝐲) Ecuación 68
𝛅𝐱 𝛅𝐲 𝛅𝐭

La ecuación (68) también puede ajustarse a la ecuación (66) teniendo en cuenta que, en
este caso, el tensor de permeabilidades solo tiene componentes k x en la dirección x, y,
ky, en la dirección y; también que α vale H(x).

La ecuación (66) es entonces la forma general de la ecuación de difusividad para flujo

lineal; y es así mismo la base para obtener las diferentes formas de tal ecuación,
dependiendo del fluido que corre a través del medio poroso cuando se tiene flujo lineal.

4.6.1. Ecuación de difusividad para flujo lineal fluido incompresible

Cuando las propiedades petrofísicas no dependen de la presión, puede plantearse la
relación.

𝐤 𝐤𝐲 𝐤
𝛅( 𝛍𝐱 ) 𝛅( 𝛍 ) 𝛅( 𝛍𝐳 ) 𝛅(𝛒𝛟)
= = = =𝟎
𝛅𝐱 𝛅𝐲 𝛅𝐳 𝛅𝐭

Y, por tanto, al realizar los operadores indicados en la ecuación (66) y al tener en

cuenta la expresión anterior, se tiene para flujo lineal en tres dimensiones

𝐤 𝛅 𝛅𝐏 𝐤𝐲 𝛅 𝛅𝐏 𝐤 𝛅 𝛅𝐏
𝛒 ( 𝐱 ) ( ) + 𝛒 ( ) ( ) + 𝛒 ( 𝐳 ) ( ) − 𝐪𝐦 = 𝟎 Ecuación 69
𝛍 𝛅𝐱 𝛅𝐱 𝛍 𝛅𝐲 𝛅𝐲 𝛍 𝛅𝐳 𝛅𝐳

Y suponiendo que kx = ky = kz se tiene:

𝛅𝟐 𝐏 𝛅𝟐 𝐏 𝛅𝟐 𝐏 𝐪𝐯 𝛍
𝟐 + 𝟐 + − =𝟎 Ecuación 70
𝛅𝐱 𝛅𝐲 𝛅𝐳 𝟐 𝐤

Y si, finalmente, no se tiene en cuenta el factor de fuentes o sumideros

𝛅𝟐 𝐏 𝛅𝟐 𝐏 𝛅𝟐 𝐏 Ecuación 71
𝛅𝐱 𝟐 + 𝛅𝐲 𝟐 + 𝛅𝐳 𝟐
=𝟎

Donde q v = q m /ρ
27
Para el flujo en dos dimensiones:

𝐤 𝛅 𝛅𝐏 𝐤𝐲 𝛅 𝛅𝐏
𝛒 ( 𝐱 ) 𝐡 ( ) + 𝛒 ( ) (𝐡 ) − 𝐪𝐦 𝐡 = 𝟎 Ecuación 72
𝛍 𝛅𝐱 𝛅𝐱 𝛍 𝛅𝐲 𝛅𝐲

Y suponiendo que kx = ky = kz se tiene:

𝛅 𝛅𝐏 𝛅 𝛅𝐏 𝛍
𝛅𝐱
(𝐡 𝛅𝐱 ) + 𝛅𝐲 (𝐡 𝛅𝐲 ) − 𝐪𝐯 𝐤 𝐡 = 𝟎 Ecuación 73

Y finalmente, si no se tiene en cuenta el factor de fuentes y sumideros y h se mantiene

constante

𝛅 𝛅𝐏 𝛅 𝛅𝐏 𝛅𝟐 𝐏 𝛅𝟐 𝐏
( ) + 𝛅𝐲 ( 𝛅𝐲 ) = 𝟎 = 𝛅𝐱𝟐 + 𝛅𝐲𝟐
𝛅𝐱 𝛅𝐱
Ecuación 74

De igual manera, para el flujo en una direcci6n, se tendrá al suponer los mismos casos
de dos y tres dimensiones.

𝐤 𝛅 𝛅𝐏
𝛒 ( 𝐱 ) (𝐀 ) − 𝐪𝐦 𝐀 = 𝟎 Ecuación 75
𝛍 𝛅𝐱 𝛅𝐱

𝛅 𝛅𝐏 𝛍
𝛅𝐱
(𝐀 𝛅𝐱 ) − 𝐪𝐦 𝐤 𝐀 = 𝟎 Ecuación 76

𝛅 𝛅𝐏 𝛅𝟐 𝐏
( ) = 𝟎 = 𝛅𝐱𝟐
𝛅𝐱 𝛅𝐱
Ecuación 77

Si las propiedades petrofísicas dependen de la presión: partiendo de la ecuación (66)

se tiene para tres dimensiones

𝛛 𝐤 𝐱 𝛛𝐏 𝛛 𝐤 𝐲 𝛛𝐏 𝛛 𝐤 𝐳 𝛛𝐏 𝛛(𝛟)
[ ]+ [ ]+ [ ] − 𝐪𝐯 = Ecuación 78
𝝏𝐱 𝛍 𝝏𝐱 𝝏𝐲 𝛍 𝝏𝐲 𝝏𝐳 𝛍 𝝏𝐳 𝝏𝐭

Y desarrollando las derivadas de La expresión anterior

𝛛 𝐱 𝐤 𝛅𝐏 𝐤 𝐱 𝛅𝟐 𝐏 𝛛 𝐤 𝐲 𝛅𝐏 𝐤 𝐲 𝛅𝟐 𝐏 𝛛 𝐳𝐤 𝛅𝐏 𝐤 𝐳 𝛅𝟐 𝐏 𝛛(𝛟)
( ) +
𝝏𝐱 𝛍 𝛅𝐱 𝛍 𝛅𝐱 𝟐
+ ( ) +
𝝏𝐲 𝛍 𝛅𝐲 𝛍 𝛅𝐲 𝟐
+ ( ) +
𝝏𝐳 𝛍 𝛅𝐳 𝛍 𝛅𝐳 𝟐
− 𝐪𝐯 =
𝝏𝐭

Recordando ahora que

𝛛 𝐤𝐱 𝛛 𝐤 𝐱 𝛅𝐏 𝛛 𝐤 𝐲 𝛛 𝐤 𝐲 𝛅𝐏 𝛛 𝐤 𝐳 𝛛 𝐤 𝐳 𝛅𝐏 𝛅𝛟 𝛅𝛟 𝛅𝐏
( ) = ( ) ; ( ) = ( ) ; ( ) = ( ) ; = ( )
𝝏𝐱 𝛍 𝛅𝐏 𝛍 𝛅𝐱 𝝏𝐲 𝛍 𝛅𝐏 𝛍 𝛅𝐲 𝝏𝐳 𝛍 𝛅𝐏 𝛍 𝛅𝐳 𝛅𝐭 𝛅𝐏 𝛅𝐭

Ecuación 79
La expresión anterior, para la ecuación de difusividad, queda así

28
𝛛 𝐱 𝐤 𝛅𝐏 𝟐 𝐤 𝐱 𝛅𝟐 𝐏 𝛛 𝐤𝐲 𝛅𝐏 𝟐 𝐤 𝐲 𝛅𝟐 𝐏 𝛛 𝐤
𝐳 𝛅𝐏 𝟐 𝐤 𝐳 𝛅𝟐 𝐏
(
𝛅𝐏 𝛍
) ( 𝛅𝐱 ) + 𝛍 𝛅𝐱 𝟐
+ ( ) ( 𝛅𝐲 ) +
𝛅𝐏 𝛍 𝛍 𝛅𝐲 𝟐
+ (
𝛅𝐏 𝛍
) ( 𝛅𝐳 ) + 𝛍 𝛅𝐳 𝟐
− 𝐪𝐯 =
𝛛(𝛟) 𝛅𝐏
𝝏𝐏 𝛅𝐭

Si se considera que

𝛅𝐏 𝟐 𝛅𝐏 𝟐 𝛅𝐏 𝟐 Ecuación 80
( ) =( ) =( ) =𝟎
𝛅𝐱 𝛅𝐲 𝛅𝐳

Lo cual es válido, porque cuando se tiene flujo estable el gradiente de presión es bajo, y
porque de acuerdo con la definición de compresibilidad de poros puede escribirse

𝛅𝛟 Ecuación 81
= 𝐂𝐏 𝛟
𝛅𝐏

Donde Cp es la compresibilidad de poros de la formación, y si se supone kx = ky = kz se

tiene

𝛅𝟐 𝐏 𝛅𝟐 𝐏 𝛅𝟐 𝐏 𝛍 𝛟𝛍𝐂𝐏 𝛅𝐏
𝛅𝐱 𝟐 + 𝛅𝐲 𝟐 + 𝛅𝐳 𝟐
− 𝐪𝐯 =
𝐤 𝐤 𝛅𝐭 Ecuación 82

y cuando no se considera el factor de fuentes y sumideros

𝛅𝟐 𝐏 𝛅𝟐 𝐏 𝛅𝟐 𝐏 𝛟𝛍𝐂𝐏 𝛅𝐏
𝟐 + 𝟐 + = Ecuación 83
𝛅𝐱 𝛅𝐲 𝛅𝐳 𝟐 𝐤 𝛅𝐭

Para dos dimensiones y aplicando el mismo procedimiento anterior se tendrá.

𝛅 𝛅𝐏 𝛅 𝛅𝐏 𝛍 𝛟𝛍𝐂𝐏 𝐡 𝛅𝐏
(𝐡 )+ (𝐡 ) − 𝐪𝐯 𝐡 = Ecuación 84
𝛅𝐱 𝛅𝐱 𝛅𝐲 𝛅𝐲 𝐤 𝐤 𝛅𝐭

Y si se considera que el espesor es constante y se desprecia el efecto de fuentes y

sumideros

𝛅 𝛅𝐏 𝛅 𝛅𝐏 𝛟𝛍𝐂𝐏 𝛅𝐏 𝛅𝟐 𝐏 𝛅𝟐 𝐏
( ) + 𝛅𝐲 ( 𝛅𝐲 ) = = + Ecuación 85
𝛅𝐱 𝛅𝐱 𝐤 𝛅𝐭 𝛅𝐱 𝟐 𝛅𝐲 𝟐

Para flujo en una dirección se tendrá

𝛅 𝛅𝐏 𝛍 𝛟𝛍𝐂𝐏 𝛅𝐏
𝛅𝐱
(𝐀 𝛅𝐱 ) − 𝐪𝐯 𝐀 𝐤 = 𝐀 𝐤 𝛅𝐭
Ecuación 86

29
Y si se considera A constante y se desprecia el efecto de fuentes y sumideros

𝛅𝟐 𝐏 𝛟𝛍𝐂𝐏 𝛅𝐏
= Ecuación 87
𝛅𝐱 𝟐 𝐤 𝛅𝐭

Con respecto a las ecuaciones de (82) a (86) es importante establecer dos aclaraciones.
La primera, el término Cp es la compresibilidad de poro; pero es normal que en lugar de
aparecer Cp aparezca Ct 0 C que representa la compresibilidad total, que está dada por:

𝐂𝐭 = 𝐂𝐩 + 𝐂𝐟 Ecuación 88

Donde Cf es la compresibilidad de fluido. Y como en el caso de los fluidos incompresibles

esta compresibilidad es cero, entonces C se convierte en Cp. Y la segunda, el termino
ϕµC / k se conoce como el inverso del coeficiente de difusividad ƞ, el cual está definido
por

𝐤
𝛈= Ecuación 89
𝛟𝛍𝐂

EI coeficiente de difusividad es una medida de la velocidad de propagación de una

perturbación de presión en un medio poroso. Se expresa en términos de área barrida
por unidad de tiempo, pues un análisis dimensional de este coeficiente muestra que
tiene unidades de área sobre tiempo

𝐤 𝐋𝟐 𝐋𝟐 Ecuación 90
[ ]= (𝐌/𝐋𝐭)(𝐋𝐭 𝟐 /𝐌)
=
𝛟𝛍𝐂 𝐭

EI coeficiente de difusividad depende de las propiedades de la roca y del fluido

almacenado en ella, y para un yacimiento dado depende del tipo de fluido contenido en
sus poros. Por ejemplo, y como se verá más adelante, el valor de ƞ será mayor cuando
en los poros hay petróleo que cuando hay gas. Por tanto, la perturbación de presión viaja
más rápido en un yacimiento cuando hay gas que cuando hay petróleo; a su vez, esto
quiere decir que el periodo transiente es más corto en un yacimiento de gas que en uno
petróleo. Finalmente, y de acuerdo con su definición, el valor de ɳ disminuye con la
presión.

4.6.2. Ecuacion de difusividad para flujo lineal fluido ligeramente compresible

Para fluido ligeramente compresible, recordando la ecuación (91), puede escribirse

𝛛𝐏 𝟏 𝛛𝛒 𝛛𝐏 𝟏 𝛛𝛒 𝛛𝐏 𝟏 𝛛𝛒 𝛛𝐏 𝟏 𝛛𝛒
= ; = ; = = Ecuación 91
𝝏𝐱 𝐂𝛒 𝝏𝐱 𝝏𝐲 𝐂𝛒 𝝏𝐲 𝝏𝐳 𝐂𝛒 𝝏𝐳 𝝏𝐭 𝐂𝛒 𝝏𝐭

30
Realizando las derivadas indicadas en la ecuación (66), suponiendo flujo en tres
dimensiones y aplicando la ecuación, (91), se tiene

𝛛 𝐤 𝛅𝐏 𝐤𝐱 𝛛𝐏 𝟐 𝛛 𝐤 𝐲 𝛅𝐏 𝐤𝐲 𝛛𝐏 𝟐 𝛛 𝐤 𝛅𝐏
( 𝐱 )+
𝝏𝐱 𝛍 𝛅𝐱 𝛍
𝐂𝐟 𝛒 ( ) + 𝛒 (
𝝏𝐱 𝝏𝐲 𝛍 𝛅𝐲
)+ 𝛍
𝐂𝐟 𝛒 ( ) + 𝛒 ( 𝐳 ) +
𝝏𝐲 𝝏𝐳 𝛍 𝛅𝐳
𝐤 𝛛𝐏 𝟐 𝛛 𝛅𝐏 𝛛𝛟
( 𝛍𝐳 ) 𝐂𝐟 𝛒 ( 𝝏𝐳 ) − 𝐪𝐦 = 𝝏𝐭 (𝛒𝛟) = 𝛟𝛒𝐂𝐟 𝛅𝐭 + 𝛒 𝝏𝐭

Recordando lo planteado para justificar la ecuación (80) y luego de dividir por la

densidad, se tiene

𝛛 𝐱 𝐤 𝛅𝐏 𝛛 𝐤 𝐲 𝛅𝐏 𝛛 𝐤 𝛅𝐏 𝛛 𝛛𝛟 𝛅𝐏
( ) + 𝝏𝐲 ( 𝛍 ) + ( 𝐳 ) − 𝐪𝐯 = (𝛒𝛟) = (𝛟𝐂𝐟 + )
𝝏𝐱 𝛍 𝛅𝐱 𝛅𝐲 𝝏𝐳 𝛍 𝛅𝐳 𝝏𝐭 𝝏𝐏 𝛅𝐭

Ecuación 92

La ecuación (92) es la forma general de la ecuación de difusividad en tres dimensiones

y de ella pueden obtenerse casos particulares, es decir, casos, cuando las propiedades
petrofísicas no dependen de la presión y cuando sí.

Cuando las propiedades petrofísicas no dependen de la presión, la ecuación (92) se

transforma en

𝐤 𝐱 𝛅𝟐 𝐏 𝐤 𝐲 𝛅𝟐 𝐏 𝐤 𝐳 𝛅𝟐 𝐏 𝛛 𝛅𝐏
𝛍
(𝛅𝐱𝟐 ) + 𝛍
(𝛅𝐲𝟐 ) + 𝛍
( 𝛅𝐳𝟐 ) −𝐪𝐯 = 𝝏𝐭 (𝛒𝛟) = (𝛟𝐂𝐟 ) 𝛅𝐭 Ecuación 93

y suponiendo un medio isotrópico, se tiene finalmente

δ2 P δ2 P δ2 P µ ϕµ δP Ecuación 94
( δx2 ) + (δy2 ) + ( δz2 ) − qv k = (Cf )
k δt

Para llegar a la ecuación (94) se ha supuesto una constante en la viscosidad del fluido,
lo cual no es correcto porque se trata de un fluido ligeramente compresible. Por tanto,
debe especificarse el nivel de presión en el que se halla esta variable, y para considerarla
constante podrá calcularse la presión promedio del intervalo de presión o calcular
también una viscosidad promedio en el mismo intervalo.

Cuando las propiedades petrofísicas no pueden considerarse constantes, se tiene a

partir de la ecuación (93).

𝛛𝐏 𝛛 𝐤𝐱 𝐤 𝐱 𝛛𝟐 𝐏 𝛛𝐏 𝛛 𝐤𝐲 𝐤 𝐲 𝛛𝟐 𝐏 𝛛𝐏 𝛛 𝐤𝐳 𝐤 𝐳 𝛛𝟐 𝐏
( )+( ) + ( )+ + ( )+( ) −
𝝏𝐱 𝝏𝐱 𝛍 𝛍 𝛛𝐱 𝟐 𝝏𝐲 𝝏𝐲 𝛍 𝛍 𝛛𝐲 𝟐 𝝏𝐳 𝝏𝐳 𝛍 𝛍 𝛛𝐳 𝟐
𝛛𝛟 𝛅𝐏
𝐪𝐯 = (𝛟𝐂𝐟 + )
𝝏𝐏 𝛅𝐭

31
Y aplicando ahora las ecuaciones (79) y (81) se tiene

𝐤 𝐱 𝛛𝟐 𝐏 𝐤 𝐲 𝛛𝟐 𝐏 𝐤 𝐳 𝛛𝟐 𝐏 𝛅𝐏
( ) 𝟐 + 𝟐 +( ) 𝟐 − 𝐪𝐯 = 𝛟(𝐂𝐟 + 𝐂𝐩 )
𝛍 𝛛𝐱 𝛍 𝛛𝐲 𝛍 𝛛𝐳 𝛅𝐭

y suponiendo que la permeabilidad es igual en todas las direcciones

𝛛𝟐 𝐏 𝛛𝟐 𝐏 𝛛𝟐 𝐏 𝛍 𝛟𝛍 𝛅𝐏
𝛛𝐱 𝟐 + 𝛛𝐲 𝟐 + 𝛛𝐳 𝟐 − 𝐪𝐯 =
𝐤 𝐤
(𝐂𝐟 + 𝐂𝐩 )
𝛅𝐭 Ecuación 95

Las ecuaciones (94) y (95) son similares y solo difieren en que, en la primera, se tiene
al lado derecho Cf y, en la segunda, se tiene (Cf+Cp).

Por tanto, ambas ecuaciones pueden escribirse en forma general como

𝛛𝟐 𝐏 𝛛𝟐 𝐏 𝛛𝟐 𝐏 𝛍 𝛟𝛍𝐂 𝛅𝐏
𝛛𝐱 𝟐 + 𝛛𝐲 𝟐 + 𝛛𝐳 𝟐 − 𝐪𝐯 =
𝐤 𝐤 𝛅𝐭
Ecuación 96

Donde C es la compresibilidad total del medio, que está dada por la ecuación (88).

Para un fluido ligeramente compresible, cuando las propiedades petrofísicas no

dependen de la presión, la compresibilidad de poro es cero y la compresibilidad total es
la del fluido. Cuando las propiedades petrofísicas si dependen de la presión, la
compresibilidad total es la suma de las compresibilidades de poro y de fluido.

Cuando se tiene flujo en dos dimensiones α, en la ecuación (66), es igual a h, y, por

tanto, siguiendo el mismo procedimiento para obtener la ecuación (92), se llega a

𝛅 𝐤 𝐱 𝛅𝐏 𝛅 𝐤 𝐲 𝛅𝐏 𝛅𝛟 𝛅𝐏
(𝐡 )+ (𝐡 ) − 𝐪𝐯 𝐡 = 𝐡 (𝐂𝐟 𝛟 + ) Ecuación 97
𝛅𝐱 𝛍 𝛅𝐱 𝛅𝐲 𝛍 𝛅𝐲 𝛅𝐏 𝛅𝐭

Así (97) es la ecuación general de difusividad para flujo lineal en dos dimensiones de un
fluido ligeramente compresible; a partir de ella pueden obtenerse las siguientes
ecuaciones, útiles cuando las propiedades petrofísicas dependan o no de la presión.
Siguiendo el mismo procedimiento para obtener las ecuaciones (94) y (95)

𝛅 𝛅𝐏 𝛅 𝛅𝐏 𝛍 𝛟𝛍 𝛅𝐏
(𝐡 )+ (𝐡 ) − 𝐪𝐯 𝐡 = 𝐡 (𝐂𝐟 ) Ecuación 98
𝛅𝐱 𝛅𝐱 𝛅𝐲 𝛅𝐲 𝐤 𝐤 𝛅𝐭

𝛅 𝛅𝐏 𝛅 𝛅𝐏 𝛍 𝛟𝛍 𝛅𝐏
(𝐡 )+ (𝐡 ) − 𝐪𝐯 𝐡 = 𝐡 (𝐂𝐟 + 𝐂𝐩 ) Ecuación 99
𝛅𝐱 𝛅𝐱 𝛅𝐲 𝛅𝐲 𝐤 𝐤 𝛅𝐭

Y si se tiene en cuenta la ecuación (88), las ecuaciones (98 y 99) pueden escribirse en
forma general

32
𝛅 𝛅𝐏 𝛅 𝛅𝐏 𝛍 𝛟𝛍𝐂 𝛅𝐏
(𝐡 )+ (𝐡 ) − 𝐪𝐯 𝐡 = 𝐡 Ecuación 100
𝛅𝐱 𝛅𝐱 𝛅𝐲 𝛅𝐲 𝐤 𝐤 𝛅𝐭

Finalmente, si se considera h constante y, además, no se tiene en cuenta el término de

fuentes y sumideros, se tendrá en el caso de la ecuación (100).

𝛅𝟐 𝐏 𝛅𝟐 𝐏 𝛟𝛍𝐂 𝛅𝐏 Ecuación 101

𝛅𝐱 𝟐 + 𝛅𝐲 𝟐
=
𝐤 𝛅𝐭

De manera similar se obtendrá las siguientes ecuaciones para flujo lineal, en una
dimensión, partiendo de la ecuación (66), y recordando que para este caso α =A.

𝛅 𝐤 𝐱 𝛅𝐏 𝛅𝛟 𝛅𝐏
(𝐀 ) − 𝐪𝐯 𝐀 = 𝐀 (𝐂𝐟 𝛟 + ) Ecuación 102
𝛅𝐱 𝛍 𝛅𝐱 𝛅𝐏 𝛅𝐭

𝛅 𝛅𝐏 𝛍 𝛍𝛟 𝛅𝐏
(𝐀 ) − 𝐪𝐯 𝐀 = 𝐀 (𝐂𝐟 ) Ecuación 103
𝛅𝐱 𝛅𝐱 𝐤 𝐤 𝛅𝐭

𝛅 𝛅𝐏 𝛍 𝛍𝛟 𝛅𝐏
(𝐀 ) − 𝐪𝐯 𝐀 = 𝐀 (𝐂𝐟 + 𝐂𝐩 ) Ecuación 104
𝛅𝐱 𝛅𝐱 𝐤 𝐤 𝛅𝐭

𝛅 𝛅𝐏 𝛍 𝛍𝛟𝐂 𝛅𝐏 Ecuación 105

(𝐀 ) − 𝐪𝐯 𝐀 = 𝐀
𝛅𝐱 𝛅𝐱 𝐤 𝐤 𝛅𝐭

La ecuación (103) es para cuando las propiedades petrofísicas no dependan de la

presión, la ecuación (104) para cuando tales propiedades si dependan y la (105) es la
forma general de ambas ecuaciones (103 y 104).

Nuevamente, si se considera A constante y se deprecian las fuentes y sumideros, se

tendrá

𝛅𝟐 𝐏 𝛍𝛟𝐂 𝛅𝐏
= Ecuación 106
𝛅𝐱 𝟐 𝐤 𝛅𝐭

4.6.3. Ecuacion de difusividad para flujo lineal: fluido compresible (gases)

Para el caso de flujo de gas se parte de la ecuación (66) y se usa, para la densidad del
fluido la ecuación dada por la ecuación de estado de los gases,

ρ = PM/ZRT Ecuación 107

Ahora bien, si se trata de gases ideales, la ecuación de estado puede presentarse

como

33
ρ = PM/RT

y retornando a la ecuación (66) para aplicarla al flujo lineal en tres dimensiones, se

tiene al introducir en ella la expresión.
𝛛 𝐏𝐌 𝛛𝐏 𝐊𝐱 𝛛 𝐏𝐌 𝛛𝐏 𝐊𝐲 𝛛 𝐏𝐌 𝛛𝐏 𝐊𝐳 𝛅 𝐏𝐌
[ ]+ [ ]+ [ ] − 𝐪𝐦 = ( ∅)
𝛛𝐱 𝐑𝐓 𝛛𝐱 𝛍 𝛛𝐲 𝐑𝐓 𝛛𝐱 𝛍 𝛛𝐳 𝐑𝐓 𝛛𝐱 𝛍 𝛅𝐭 𝐑𝐓

𝛛 𝛛𝐏 𝐊𝐱 𝛛 𝛛𝐏 𝐊𝐲 𝛛 𝛛𝐏 𝐊𝐳 𝛅
[𝐏 ]+ [𝐏 ]+ [𝐏 ] − 𝐪𝐦 = (𝐏∅)
𝛛𝐱 𝛛𝐱 𝛍 𝛛𝐲 𝛛𝐱 𝛍 𝛛𝐳 𝛛𝐱 𝛍 𝛅𝐭

Que puede llevarse a

𝛛 𝛛𝐏𝟐 𝐊 𝐱 𝛛 𝛛𝐏𝟐 𝐊 𝐲 𝛛 𝛛𝐏𝟐 𝐊 𝐳 𝐑𝐓 ∅ 𝛅∅ 𝛅𝐏𝟐

[ ]+ [ ]+ [ ] − 𝟐𝐪𝐦 =( + )
𝛛𝐱 𝛛𝐱 𝛍 𝛛𝐲 𝛛𝐲 𝛍 𝛛𝐳 𝛛𝐳 𝛍 𝐌 𝐏 𝛅𝐏 𝛅𝐏

Ecuación 108

A partir de la ecuación (108) pueden obtenerse expresiones para el flujo de gas ideal,
pero esto depende de si las propiedades petrofísicas son independientes o no de la
presión. En el primer caso, en un medio isotrópico y calculando la viscosidad del gas bajo
una presión inicial promedio y constante, la ecuación será

𝝏𝟐 𝑷𝟐 𝝏𝟐 𝑷𝟐 𝝏𝟐 𝑷𝟐 𝑹𝑻 𝝁 ∅𝝁 𝟏 𝜹𝑷𝟐
[ + + − 𝟐𝒒𝒎 = ( )
𝝏𝒙𝟐 𝝏𝒚𝟐 𝝏𝒛𝟐 𝑴 𝐊 𝑲 𝑷 𝜹𝒕

Y como 1/P es la compresibilidad del gas ideal, la ecuación anterior es también

𝝏𝟐 𝑷𝟐 𝝏𝟐 𝑷𝟐 𝝏𝟐 𝑷𝟐 𝑹𝑻 𝝁 ∅𝝁𝑪𝒇 𝜹𝑷𝟐
[ + + 𝟐 ] − 𝟐𝒒𝒎 = Ecuación 109
𝝏𝒙𝟐 𝝏𝒚𝟐 𝝏𝒛 𝑴 𝐊 𝑲 𝜹𝒕

Si k y ϕ no pueden considerarse constantes con la presión, después de aplicar los

operadores indicados a la ecuación (108) y luego las relaciones entre las ecuaciones (80
y 81), se tiene.

𝐊 𝐱 𝛛𝟐 𝐊 𝐲 𝛛𝟐 𝐊 𝐳 𝛛𝟐 𝐑𝐓 𝟏 𝛅𝐏 𝟐
[ 𝟐
+ 𝟐
+ 𝟐
− 𝟐𝐪𝐦 = ∅ ( + 𝐂𝐩)
𝛍 𝛛𝐱 𝛍 𝛛𝐲 𝛍 𝛛𝐳 𝐌 𝐏 𝛅𝐭

Y al suponer un medio isotrópico y, a su vez, aplicar la relación de la ecuación (88) y

calcular la viscosidad bajo una presión inicial o promedio como constante, se tiene

𝜹 𝟐 𝑷𝟐 𝝏𝟐 𝑷𝟐 𝝏𝟐 𝑷𝟐 𝑹𝑻 𝝁 ∅𝝁𝑪 𝜹𝑷𝟐
[ + + ] − 𝟐𝒒𝒎 = Ecuación 110
𝜹𝒙𝟐 𝝏𝒚𝟐 𝝏𝒛𝟐 𝑴 𝐊 𝑲 𝜹𝒕

34
Y si igualmente no se tiene en cuenta el efecto de fuentes 0 sumideros

𝛅𝟐 𝐏 𝟐 𝛛𝟐 𝐏𝟐 𝛛𝟐 𝐏𝟐 ∅𝛍𝐂 𝛅𝐏𝟐
[ + + ]= ) Ecuación 111
𝛅𝐱 𝟐 𝛛𝐲 𝟐 𝛛𝐳 𝟐 𝐊 𝛅𝐭

También, para el flujo en dos dimensiones, a partir de la ecuación (66) se tiene

𝛅 𝐤 𝐱 𝐡 𝛅𝐏𝟐 𝛅 𝐤 𝐲 𝐡 𝛅𝐏𝟐 𝐑𝐓 ∅ 𝛅∅ 𝛅𝐏𝟐

( )+ ( ) − 𝟐𝐪𝐦 𝐡 = 𝐡( + ) Ecuación 112
𝛅𝐱 𝛍 𝛅𝐱 𝛅𝐲 𝛍 𝛅𝐲 𝐌 𝐏 𝛅𝐏 𝛅𝐭

A partir de la ecuación (112) pueden obtenerse las siguientes expresiones, dependiendo

de que las propiedades petrofísicas sean respectivamente independientes o no de la
presión.

𝛅 𝛅𝐏𝟐 𝛅 𝛅𝐏 𝟐 𝐑𝐓 ∅𝛍 𝟏 𝛅𝐏𝟐
(𝐡 )+ (𝐡 ) − 𝟐𝐪𝐦 𝐡 =𝐡 ( ) Ecuación 113
𝛅𝐱 𝛅𝐱 𝛅𝐲 𝛅𝐲 𝐌 𝐤 𝐏 𝛅𝐭

𝛅 𝛅𝐏𝟐 𝛅 𝛅𝐏 𝟐 𝐑𝐓 ∅𝛍 𝟏 𝛅𝐏𝟐
(𝐡 )+ (𝐡 ) − 𝟐𝐪𝐦 𝐡 =𝐡 ( + 𝐂𝐩) Ecuación 114
𝛅𝐱 𝛅𝐱 𝛅𝐲 𝛅𝐲 𝐌 𝐤 𝐏 𝛅𝐭

Para llegar a las ecuaciones (113 y 114) se aplicaron las mismas suposiciones que para
llegar a las ecuaciones (109 y 110); además, teniendo en cuenta (88), la ecuación general
para flujo lineal de gas ideal en dos dimensiones es.

𝛅 𝛅𝐏𝟐 𝛅 𝛅𝐏 𝟐 𝐑𝐓 ∅𝛍𝐂 𝛅𝐏𝟐

(𝐡 )+ (𝐡 ) − 𝟐𝐪𝐦 𝐡 =𝐡 Ecuación 115
𝛅𝐱 𝛅𝐱 𝛅𝐲 𝛅𝐲 𝐌 𝐤 𝛅𝐭

Si se desprecia el efecto de fuentes o sumideros y, además, se considera h constante

se tiene

𝛅𝟐 𝐏 𝟐 𝛅𝟐 𝐏 𝟐 ∅𝛍𝐂 𝛅𝐏 𝟐
+ = Ecuación 116
𝛅𝐱 𝟐 𝛅𝐲 𝟐 𝐤 𝛅𝐭

Finalmente, si el flujo es en una dimensión, entonces a partir de la ecuación (66) y

recordando que en este caso α=A(x), se tiene

𝛅 𝐤 𝐱 𝐀(𝐱) 𝛅𝐏𝟐 𝐑𝐓 ∅ 𝛅∅ 𝛅𝐏𝟐

[ ] − 𝟐𝐪𝐦 𝐀 = 𝐀( + ) Ecuación 117
𝛅𝐱 𝛍 𝛅𝐱 𝐌 𝐏 𝛅𝐏 𝛅𝐭

A partir de la ecuación (117); y planteando las mismas suposiciones de los casos

anteriores, se logran las siguientes ecuaciones para casos particulares de flujo de gas
ideal en un medio poroso. Primero, para cuando las propiedades petrofísicas no
dependan de la presión

35
𝛅 𝛅𝐏𝟐 𝛍 𝐑𝐓 𝐀∅𝛍 𝟏 𝛅𝐏𝟐 𝐀∅𝛍𝐂𝐟 𝛅𝐏𝟐 Ecuación 118
[𝐀(𝐱) ] − 𝟐𝐪𝐦 𝐀 = ( ) =
𝛅𝐱 𝛅𝐱 𝐤 𝐌 𝐤 𝐏 𝛅𝐭 𝐤 𝛅𝐭

Y para cuando las propiedades petrofísicas dependan de ella

𝛅 𝛅𝐏𝟐 𝛍 𝐑𝐓 𝐀∅𝛍 𝟏 𝛅𝐏𝟐 𝐀∅𝛍𝐂𝐟 𝛅𝐏𝟐

[𝐀(𝐱) ] − 𝟐𝐪𝐦 𝐀 = ( + 𝐂𝐩 ) =
𝛅𝐱 𝛅𝐱 𝐤 𝐌 𝐤 𝐏 𝛅𝐭 𝐤 𝛅𝐭

Ecuación 119

Cuando el área puede considerarse constante y se desprecia el efecto de fuentes o

sumideros, se tiene

𝛅𝟐 𝐏 𝟐 ∅𝛍 𝟏 𝛅𝐏𝟐 ∅𝛍𝐂 𝛅𝐏𝟐

= ( + 𝐂𝐩) = Ecuación 120
𝛅𝐱 𝟐 𝐤 𝐏 𝛅𝐭 𝐤 𝛅𝐭

En el caso contrario, cuando el gas es real, la ecuación para la densidad es p =PM/ZRT

y al reemplazar la densidad por esta expresión en la ecuación (66) se tiene, para el caso
de flujo en tres dimensiones

𝛛 𝐏𝐌 𝐊𝐱 𝛛𝐏 𝛛 𝐏𝐌 𝐊𝐲 𝛛𝐏 𝛛 𝐏𝐌 𝐊𝐳 𝛛𝐏 𝛛 𝐏𝐌
[ ]+ [ ]+ [ ]= ( ∅) Ecuación 121
𝛛𝐱 𝐑𝐓 𝛍 𝛛𝐱 𝛛𝐲 𝐑𝐓 𝛍 𝛛𝐱 𝛛𝐳 𝐑𝐓 𝛍 𝛛𝐱 𝛛𝐭 𝐙𝐑𝐓

𝛛 𝐏𝐤 𝐱 𝛛𝐏 𝛛 𝐏𝐤 𝐲 𝛛𝐏 𝛛 𝐏𝐤 𝐳 𝛛𝐏 𝛛 𝐏∅
[ ]+ [ ]+ [ ]= ( ) Ecuación 122
𝛛𝐱 𝛍𝐙 𝛛𝐱 𝛛𝐲 𝛍𝐙 𝛛𝐲 𝛛𝐳 𝛍𝐙 𝛛𝐱 𝛛𝐭 𝐙

La presencia del término Z en la ecuación (122) la torna más compleja para su

manipulación y obtención de una ecuación diferencial parcial de segundo orden, similar
a las obtenidas hasta ahora. Sin embargo, si se plantean algunas suposiciones, es
igualmente posible obtener ecuaciones tan simples como las anteriores.

De este modo; si se supone que el factor (µZ) es constante, en la ecuación (122);

después de algunas simplificaciones ella se convierte en.

𝛛 𝛛𝐏 𝟐 𝛛 𝛛𝐏 𝟐 𝛛 𝛛𝐏 𝟐 𝐑𝐓 𝛅∅ 𝟏 𝟏 𝛅𝐙 𝛅𝐏 𝟐
𝛛𝐱
(𝐤 𝐱 𝛛𝐱
) + 𝛛𝐲 (𝐤 𝐲 𝛛𝐲
) + 𝛛𝐳 (𝐤 𝐳 𝛛𝐳
) − 𝟐𝐪𝐦 𝐌
= 𝛍(
𝛅𝐏
− ∅( −
𝐏 𝐙 𝛅𝐏
)) 𝛅𝐭

Ecuación 123
La ecuación (123) puede ser válida en un intervalo de presión en donde, al aumentar
esta, Z disminuya y la viscosidad aumente. Este comportamiento se da bajo presiones
bajas y hasta un valor del orden de 2.000 Lpc. Por tanto, esta suposición puede aceptarse
para presiones menores de unas 2.000 Lpc.

36
La ecuación (123) es similar a la ecuación (108), solo que en lugar del término (1/P) de
la ecuación (108), que es la compresibilidad del gas ideal, se tiene la expresión 1/P -
(1/Z) (δZ/δP), que es la compresibilidad de un gas real. Por tanto, a partir de ella pueden
obtenerse ecuaciones similares a las ecuaciones (109, 110 y 111) mediante las mismas
suposiciones.

En conclusión, la ecuación de difusividad para un gas real, cuando puede realizarse la

suposición de que (µZ) es constante (0 sea, para presiones menores de unas 2.000 Lpc),
es idéntica a la ecuación de difusividad para un gas ideal.

Solo que, al calcular la compresibilidad del gas, debe tenerse en la cuenta que, para un
gas real, la compresibilidad se calcula mediante la ecuación (124)

𝟏 𝟏𝛛𝐙
𝐂𝐠 = [ − ] Ecuación 124
𝐏 𝐙𝛛𝐏

De igual manera, pueden conseguirse las ecuaciones para una y dos dimensiones
similares a las ecuaciones (112) y (117).

Si en la ecuación (122) se considera constante el termino (P/µZ), esta se convierte en

𝛅 𝛅𝐏 𝟐 𝛅 𝛅𝐏 𝟐 𝛅 𝛅𝐏 𝟐 𝐑𝐓 𝛍𝐙 𝛅∅ 𝟏 𝟏 𝛅𝐙 𝛅𝐏
𝛅𝐱
(𝐤 𝐱 𝛅𝐱
) + 𝛅𝐲 (𝐤 𝐲 𝛅𝐲
) + 𝛅𝐳 (𝐤 𝐳 𝛅𝐳
) − 𝐪𝐦 𝐌 𝐏
= 𝛍(
𝛅𝐏
− ∅( −
𝐏 𝐙 𝛅𝐏
)) 𝛅𝐭

Ecuación 125

La ecuación (125) es similar a la ecuación para un fluido ligeramente compresible, como

se señala en (92), solo que en lugar del término Cf, que es la compresibilidad constante
del fluido, se tiene 1/P-(1/Z) (δZ/δP) que es la compresibilidad del gas real. Por tanto, a
partir de la ecuación (125) pueden obtenerse, por medio de las mismas suposiciones,
ecuaciones similares a las (94, 95 y 96). De igual manera se pueden obtener expresiones
para flujo en una y dos dimensiones.

En resumen, la ecuación de difusividad para flujo de gases reales cuando puede

suponerse de (P/µZ) constante, es idéntica a la ecuación de difusividad para un fluido
ligeramente compresible, solo que para el gas la compresibilidad se debe calcular con la
ecuación (124).

La suposición de (P/µZ) constante puede hacerse cuando al aumentar P tanto la

viscosidad como Z aumentan esto ocurre bajo presiones mayores a unas 3.000 Lpc, 0
sea, cuando el gas empieza a tener comportamientos que se aproximan al de los
líquidos. De allí que la ecuación de difusividad resulte similar a la de estos.

37
Es importante insistir en las suposiciones para obtener algunas formas de la ecuación de
difusividad para gases. Primero, la suposición de que (δP/δs)2 ≈ 0 es menos valida que
en el caso de los líquidos, pues para los gases el gradiente de presión en la dirección de
flujo no es tan pequeño como en el otro caso. Segundo, la viscosidad del gas no puede
considerarse constante y, por tanto, el valor de viscosidad que aparece en el coeficiente
de difusividad de las ecuaciones de difusividad para gases debe medirse normalmente
a la presión inicial o a la presión promedio del intervalo en el que se está trabajando.
Tercero, la ecuación de difusividad para gases es aún menos lineal que la ecuación para
líquidos, pues el termino (µC) es mucho más dependiente de Ia presión en el caso de un
gas que en el otro.

Las ecuaciones de difusividad obtenida para gases reales, como indican las ecuaciones
(123) y (125), implican suposiciones y, además, son altamente no lineales. Para no hacer
las suposiciones de (µZ) 0 (P/µZ) constantes y para tratar de obtener una ecuación
menos no lineal que las anteriores, se introduce en la deducción de la ecuación de
difusividad para gases el concepto de pseudopresión. La función pseudopresión está
definida por
𝐏 𝟐𝐏
𝐦(𝐏) = ∫𝐏𝐛 𝐝𝐏 Ecuación 126
𝛍𝐳

donde m (P) se conoce como la función pseudopresión y Pb es una presión de referencia

0 base; esta puede ser la presión correspondiente a la presión normal 0 estándar.

De acuerdo con la ecuación (126) pueden tenerse las relaciones

𝟐𝐏
𝐝𝐦(𝐏) = 𝐝𝐩 Ecuación 127
𝛍𝐳

𝐝𝐏 𝛍𝐙 𝐝(𝐦)
= Ecuación 128
𝐝𝐢 𝟐𝐏 𝐝𝐢

Donde i puede ser cualquier variable que afecte la presión.

Ahora, llevando a la ecuación (66) la ecuación (128) y la expresi6n de densidad para un

gas real, se tiene

𝛛 𝐏𝐌 𝐤 𝐱 𝛍𝐙 𝛅𝐦(𝐏) 𝛛 𝐏𝐌 𝐤 𝐲 𝛍𝐙 𝛅𝐦(𝐏) 𝛛 𝐏𝐌 𝐤 𝐳 𝛍𝐙 𝛅𝐦(𝐏)

( )+ ( ) + ( ) − 𝐪𝐦 =
𝛛𝐱 𝐙𝐑𝐓 𝛍 𝟐𝐏 𝛅𝐱 𝛛𝐲 𝐙𝐑𝐓 𝛍 𝟐𝐏 𝛅𝐲 𝛛𝐳 𝐙𝐑𝐓 𝛍 𝟐𝐏 𝛅𝐳
𝛅 𝐏𝐌 𝐌 𝛅 𝐏∅ 𝛍𝐙 𝛅𝐦(𝐏)
( ∅) = ( )
𝛅𝐭 𝐙𝐑𝐓 𝐑𝐓 𝛅𝐏 𝐙 𝟐𝐏 𝛅𝐭

Y después de las simplificaciones del caso

38
𝛛 𝛅𝐦(𝐏) 𝛛 𝛅𝐦(𝐏) 𝛛 𝛅𝐦(𝐏) 𝐑𝐓 𝛅∅
(𝐤 𝐱 )+ (𝐤 𝐲 )+ (𝐤 𝐳 ) − 𝟐𝐪𝐦 = 𝛍( +
𝛛𝐱 𝛅𝐱 𝛛𝐲 𝛅𝐲 𝛛𝐳 𝛅𝐳 𝐌 𝛅𝐏
𝟏 𝟏 𝛅𝐙 𝛅𝐦(𝐏)
∅( − ))
𝐏 𝐙 𝛅𝐏 𝛅𝐭

Y si recordamos la ecuación (124), la expresión anterior se escribe como

𝛛 𝛅𝐦(𝐏) 𝛛 𝛅𝐦(𝐏) 𝛛 𝛅𝐦(𝐏) 𝐑𝐓 𝜹∅ 𝜹𝒎(𝑷)

𝛛𝐱
(𝐤 𝐱 𝛅𝐱
) + 𝛛𝐲 (𝐤 𝐲 𝛅𝐲
) + 𝛛𝐳 (𝐤 𝐳 𝛅𝐳
) − 𝟐𝐪𝐦 𝐌
= 𝝁(
𝜹𝑷
+ ∅(𝑪𝒇))
𝜹𝒕

Ecuación 129

La ecuación (129) es similar a la ecuación (92), solo que esta última se expresa en
términos de la presión y aquella en términos de la pseudopresión; así que una solución
para la ecuación (92) en términos de la presión también lo será para la ecuación (129)
en términos de la pseudopresión. La ecuación (129) es también similar a las ecuaciones
(123 y 125), pero para llegar a estas hubo que realizar respectivamente las suposiciones
de (µZ) y (P/µZ) como constantes; en cambio, para llegar a la ecuación (129) no hubo
necesidad de ninguna suposición y por tanto es más general que las anteriores.

De la ecuación (129) se obtendrán ecuaciones similares a las (94, 95 y 96), dependiendo

de si se considera que la presión afecta o no las propiedades petrofísicas. Además,
cuando las propiedades petrofísicas dependen de la presión, debe hacerse la suposición
de que (δm(P)/δs)2 ≈ 0 donde s es cualquier dirección de flujo.

De igual manera, si se considera el flujo en una 0 dos dimensiones se llegan

respectivamente a ecuaciones similares a (97 y 102).

Cuando se trabaja con m (P) debe haber la posibilidad de convertir presión a m (P), o lo
contrario, m (P) a presión. Para ello, debe contarse con un gráfico de m (P) por presión
mediante el siguiente procedimiento.

1. Se toma un intervalo amplio de presión, dependiendo de la presión a Ia que se

encuentre el yacimiento que se está analizando, normalmente puede ser desde 14,7 Psi
hasta la presión del yacimiento Pi.

Este intervalo se subdivide en otros de unas 50 Psi.

2. Para la presión inicial y final de cada uno de los intervalos, en que se subdividió el
intervalo de presión, se calcula (2P/µZ) y luego se grafica

(2P/µZ) vs P.

39
3. EI valor de m (P) a Ia presión final Pn del intervalo n es el área bajo la curva de (2P/µZ)
vs P entre 14,7 Lpc y Pn y puede obtenerse aplicando el método trapezoidal cuya fórmula
general es

𝛁𝐏 𝟐𝐏 𝟐𝐏 𝟐𝐏
𝐦(𝐏) = [( ) + 𝟐 ∑𝐧−𝟏
𝐢=𝟏 ( ) + ( ) ] Ecuación 130
𝟐 𝛍𝐳 𝐏 𝛍𝐳 𝐏 𝛍𝐳 𝐏
𝐨 𝐢 𝐧

Dónde: (2P/µZ) po es el valor de (2P/µZ) evaluado a 14,7 Lpc y (2P/µ) pi es el valor de

(2P/µZ) evaluado bajo la presión final de cada uno de los n intervalos de amplitud P
comprendidos en el intervalo 14,7 Psi.

4. Se repite el procedimiento anterior hasta que se haga el recorrido de todos los

intervalos en que se dividió el intervalo 14,7 Psi y luego puede dibujarse un gráfico de m
(P) vs P.

5. Con el grafico obtenido puede lograrse el valor de m (P) correspondiente a una presión
dada o la presión correspondiente a un m (P) dado. Pero es más práctico, especialmente
si se desea programar el procedimiento de convertir P a m (P), o lo contrario, obtener
mediante regresión una relación entre las dos variables.

4.7. CONDICIONES NECESARIAS PARA LA SOLUCION

Las soluciones de la ecuación de difusividad para el periodo pseudoestable y para un

tiempo determinado, se conocen como ecuaciones de flujo o afluencia. Estas se pueden
obtener a partir de alguna de las formas de la ecuación de difusividad para flujo radial.
Debe tenerse en cuenta que el periodo pseudoestable se caracteriza por δP/δt = 0 o
δP/δt = constante dependiendo esto si el flujo es continuo o pseudoestable
respectivamente.

𝛛 𝛛𝐏 ∅𝛍𝐂 𝛛𝐏
[𝐫 ] =
𝛛𝐫 𝛛𝐫 𝐤 𝛛𝐭
Las suposiciones hechas para obtener la ecuación anterior son:
• Flujo en una sola fase
• Fluido ligeramente compresible
• Ley de Darcy válida
• Viscosidad constante
• Medio homogéneo e isótropo
• Efectos de gravedad despreciables
• Gradientes de presión pequeños en el yacimiento
• Flujo isotérmico

40
La Ecuación anterior es una ecuación diferencial en derivadas parciales y requiere para
su solución establecer condiciones iníciales y de frontera.
La condición inicial es la siguiente:
𝐩(𝐫, 𝐭 = 𝟎) = 𝐩𝐢
La condición de frontera son las siguientes:

a) Gasto especifico
b) Presión especificada

Para obtener la solución de la ecuación de difusión aplicable para un caso particular es

necesario definir bajo qué condiciones actúan las fronteras; es decir, si la frontera es
impermeable o se mantiene constante o si hay producción a través de la frontera.

4.8. SOLUCIÓN DE LA ECUACIÓN DE DIFUSIVIDAD RADIAL

Es posible encontrar una solución general que dependerá de la geometría del medio
poroso, del tipo de fluido y de las condiciones iníciales y de límite que se apliquen. En
estos casos, y con el fin de tener una ecuación diferencial más sencilla, se trabaja con la
ecuación de difusividad en variables adimensionales las cuales tienen el objetivo de
reducir el número de parámetros para resolver dicha ecuación. Se hace uso de las
variables adimensionales para caracterizar las soluciones que representan el flujo de
fluidos en un medio poroso.

Existen tres parámetros adimensionales importantes para flujo radial: 𝒑𝑫 ,𝒕𝑫 ,𝒓𝑫 .

lineal Radial Esférico

Presión 𝒑 𝒑 𝒑
Adimensional 𝒌𝑨(𝑷𝒊− 𝑷) 𝟐𝝅𝒌𝒉(𝑷𝒊− 𝑷) 𝟒𝝅𝒌𝒓𝒘 (𝑷𝒊− 𝑷)
𝑫= 𝑫= 𝑫=
𝒒𝝁𝑳 𝒒𝝁 𝒒𝝁

Tiempo 𝒕 𝒕
Adimensional
𝒕 𝒌𝒕 𝒌𝒕 𝒌𝒕
𝑫= 𝑫= 𝑫=
∅𝝁𝑪𝒕 𝑳𝟐 ∅𝝁𝑪𝒕 𝒓𝟐𝒘 ∅𝝁𝑪𝒕 𝒓𝟐𝒘
Posición
Adimensional
𝒙𝑫 = 𝑿⁄𝑳 𝒓𝑫 = 𝒓⁄𝒓𝒘 𝒓𝑫 = 𝒓⁄𝒓𝒘

41
 lineal Radial Esférico
Presión 𝒌𝑨(𝑷𝒊− 𝑷) 𝑷𝑫 𝑷
Adimensional 𝑷𝑫 = 𝟐𝝅𝒌𝒉(𝑷𝒊− 𝑷) 𝟒𝝅𝒌𝒓𝒘 (𝑷𝒊− 𝑷)
𝟖𝟖𝟕. 𝟐𝒒𝑩𝝁𝑳 = 𝑫=
𝟏𝟒𝟏.𝟐𝒒𝑩𝝁 𝟕𝟎.𝟔𝒒𝑩𝝁

Tiempo 𝒕𝑫
Adimensional 𝟐.𝟔𝟑𝟕∗𝟏𝟎−𝟒 𝒌𝒕 𝒕𝑫 𝟐.𝟔𝟑𝟕∗𝟏𝟎−𝟒 𝒌𝒕 𝒕 𝟐.𝟔𝟑𝟕∗𝟏𝟎−𝟒 𝒌𝒕
= = 𝑫=
∅𝝁𝑪𝒕 𝑳𝟐 ∅𝝁𝑪𝒕 𝒓𝟐𝒘 ∅𝝁𝑪𝒕 𝒓𝟐𝒘
Posición
Adimensional
𝒙𝑫 = 𝑿⁄𝑳 𝒓𝑫 = 𝒓⁄𝒓𝒘 𝒓𝑫 = 𝒓⁄𝒓𝒘

Para el caso de un yacimiento circular con flujo radial, yacimiento homogéneo e infinito
y produciendo a gasto constante.

Se presenta la solución para las consideraciones mencionadas anteriormente partiendo

de:

𝟏 𝛛 𝛛𝐩 𝛟𝛍𝐂𝐭 𝛛𝐩
(𝐫 ) =
𝐫 𝛛𝐫 𝛛𝐫 𝐤 𝛛𝐭
La ecuación anterior representa la ecuación de difusión en coordenadas cilíndricas para
la solución supondremos las siguientes condiciones de frontera.

𝐏(𝐫, 𝟎) = 𝐩𝐢 Condición inicial.

𝛛𝐩 𝐪𝛍
(𝐫 ) =− Condición de frontera que considera gasto constante en el
𝛛𝐫 𝐫=𝐫𝐰 𝟐𝛑𝐤𝐡
pozo.

𝐥𝐢𝐦 𝒑(𝒓, 𝒕) = 𝒑𝒊 Condición de frontera infinita.

𝒓→∞

Considerando las condiciones de frontera anteriores la solución de la ecuación de

difusión es la siguiente:

𝐪𝛍 𝛟𝛍𝐂𝐭𝐫 𝟐
𝐩(𝐫, 𝐭) = 𝐩𝐢 + 𝐄 (− )
𝟒𝛑𝐤𝐡 𝐈 𝟒𝐤𝐭

Donde 𝐸1 (x) es la función Integral Exponencial y se define como:

42
 ∞
𝐞−𝐮
−𝐄𝐢 (−𝐱) = ∫ 𝐝𝐮
𝐮
𝐱

La ecuación anterior representa el comportamiento de presión en cualquier punto del

yacimiento y a cualquier tiempo, esta solución de la ecuación de difusión es de las más
usadas en el análisis de pruebas de presión.

Expresada en variables adimensionales:

𝟏 𝐫𝐃 𝟐
𝐩𝐃 (𝐫𝐃 , 𝐭 𝐃 ) = − 𝐄𝟏 ( )
𝟐 𝟒𝐭 𝐃

4.8.1. SOLUCIÓN PARA FLUJO ESTABLE

La ecuación para flujo estable es δP/δt=0

Y cuando se evalúa en r = re se tiene:

𝐪𝛍 𝐫𝐞
𝐏𝐞 − 𝐏𝐰𝐟 = 𝐥𝐧
𝟐𝛑𝐤𝐡 𝐫𝐰
Que es la ecuación para flujo estable, idéntica a la ecuación de Darcy para flujo radial.

Al igual que en la ecuación el único termino difícil de conocer es Pₑ por tanto, es más
conveniente hablar de P* que puede obtenerse de pruebas de presión, y es en si la
presión promedio de todo el yacimiento. Siguiendo un procedimiento análogo para
presentar la ecuación de flujo pseudoestable en términos de la presión promedio del
yacimiento, como consta en la ecuación, se tiene para el caso de flujo estable:
𝐪𝛍 𝐫𝐞 𝟏
𝐏 ∗ − 𝐏𝐰𝐟 = [𝐥𝐧 − ]
𝟐𝛑𝐤𝐡 𝐫𝐰 𝟐

Cuando se tienen unidades prácticas del sistema Ingles se transforma en:

𝐪𝛍𝐁 𝐫
𝐏 ∗ − 𝐏𝐰𝐟 = 𝟏𝟒𝟏, 𝟐 [𝐥𝐧 𝐞 ]
𝐤𝐡 𝐫𝐰

4.8.2. SOLUCIÓN PARA FLUJO PSEUDO ESTABLE

Las soluciones de la ecuación de difusividad para el periodo pseudoestable y para un

tiempo determinado se conocen como ecuaciones de flujo o afluencia. Estas se pueden
obtener a partir de alguna de las formas de la ecuación de difusividad para flujo radial
planteadas anteriormente. Debe tenerse en cuenta que el periodo pseudoestable se
caracteriza por δP/δt= 0 o δP/δt = Cte. dependiendo esto de si el flujo es continuo o

43
pseudoestable respectivamente. Se usará la ecuación de la cual se tomará el termino
(Cp + Cf) como C, se despreciarán los efectos de fuentes o sumideros y se considerara
h constante

𝟏 𝛛 𝛛𝐏 ∅𝛍𝐂 𝛛𝐏
[𝐫 ]=
𝐫 𝛛𝐫 𝛛𝐫 𝐤 𝛛𝐭

La expresión para definir el flujo pseudoestable es δP/δt=Constante

Y las de limite en el exterior son de flujo cero; pero en el pozo, limite interior, puede darse
que la tasa de flujo sea constante o que la presión sea constante.

Se tomará la tasa de flujo constante. Entonces, en el exterior

𝛅𝐏 𝛅𝐏 𝐪𝛍
( ) = 𝟎 𝐲 𝐞𝐧 𝐞𝐥 𝐢𝐧𝐭𝐞𝐫𝐢𝐨𝐫 (𝐫 ) = = 𝐜𝐨𝐧𝐬𝐭𝐚𝐧𝐭𝐞
𝛅𝐫 𝐫=𝐫𝐨 𝛅𝐫 𝐫=𝐫𝐰 𝟐𝛑𝐤𝐡

Cuando se tiene flujo pseudoestable, gráficamente el comportamiento de la presión es

como se muestra en la figura y como δP/δt=Constante el valor de esta constante puede
obtenerse cuando ha transcurrido un tiempo t de producción también constante, y del
yacimiento ha salido entonces una cantidad de fluido dV que, de acuerdo con la ecuación
de compresibilidad, puede aproximarse a:

𝐝𝐕 = −𝐂𝐕𝐝𝐏 = −𝐂𝛑𝐫𝐞𝟐 ∅𝐝𝐏

Comportamiento de la presión en el periodo pseudoestable

44
 𝐪𝛍 𝐫𝐞 𝟏
𝐏𝐞 − 𝐏𝐰𝐟 = [𝐥𝐧 − ]
𝟐𝛑𝐤𝐡 𝐫𝐰 𝟐

Que es la ecuación general para flujo al pozo en estado pseudoestable, y que es similar
a la ecuación para flujo en estado estable (ecuación de Darcy para flujo radial).

Que finalmente se tendrá

𝐪𝛍 𝐫𝐞 𝟑
𝐏 ∗ − 𝐏𝐰𝐟 = [𝐥𝐧 − ]
𝟐𝛑𝐤𝐡 𝐫𝐰 𝟒

5. CONCLUSION Y RECOMENDACION

• Fue posible modelar el comportamiento del flujo de fluidos bajo una serie de
ecuaciones que me permiten determinar el movimiento de estos fluidos, llegando
a la conclusión que son los mismos principios fundamentales que regulan el
movimiento de fluidos en tuberías.

• Los principios fundamentales que permiten representar el movimiento de los

fluidos en un yacimiento son la Conservación de la Masa, Momento y Energía.
Como el fenómeno es enfocado al flujo de fluidos a través de un medio poroso, el
principio de la Conservación del Momento es reemplazado por una ecuación más
experimental como lo es la Ley de Darcy. Adicionalmente a estas relaciones que
se han establecido hasta ahora, hay que tener muy en cuenta las propiedades
físicas de los fluidos del sistema, pues deberían estar representados como
funciones de las variables independientes.

6. BIBLIOGRAFIA

• Ingeniería Aplicada de Yacimientos Petrolíferos, B.C. Craft y M.F. Hawkins.

• Fundamentos de Ingeniería de Yacimientos, Magdalena Paris de Ferrer
• Flujo de fluidos en medios porosos, Ing. Gabriel J. Colmont.
• Evaluación de Yacimientos de Hidrocarburos, Abel Naranjo Agudelo.
• Trabajo de investigación de la escuela politécnica nacional ingeniería de
petróleos-ecuaciones que describen el flujo de fluidos en medios porosos.

45