Fundamentos de Ingenieria de Yacimientos-Fredy Escobar
Fundamentos de Ingenieria de Yacimientos-Fredy Escobar
Fundamentos de Ingenieria de Yacimientos-Fredy Escobar
FUNDAMENTOS
DE INGENIERIA
DE YACIMIENTOS
Autor:
© de esta edición
Primera edición:
Xxxxx de 200?
ISBN xxx-yyyy-yy-z
Diseño de Portada:
xxxxxxxxxxxxxxxxx
Fotografía portada:
xxxxxxxxxxxxxxxxxxxxxxxxx
Diseño y diagramación:
María Constanza Cardoso Perdomo
Impresión y encuadernación
Editora Guadalupe Ltda..
INTRODUCCIÓN
El contenido del libro se ha dividido en ocho capítulos. El capítulo 1 presenta una breve
revisión de los conceptos geológicos que definen y clasifican los yacimientos de
hidrocarburos. Al igual que una clasificación actualizada de los yacimientos de acuerdo al
estado de los fluidos. Además, se estudiarán los regímenes de presión de un yacimiento
(diagramas profundidad-presión) con su importancia en la caracterización de yacimientos y
las escalas de medición que se usan en Ingeniería de Yacimientos. El capítulo dos se centra
en el entendimiento de las propiedades físicas del medio poroso, la ley de Darcy con sus
aplicaciones y limitaciones, las clasificaciones de la permeabilidad, los problemas
asociados con la geometría del sistema, el daño a la formación, la ecuación de difusividad,
propiedades inherentes a la roca y fluido: presión capilar, función J de Leverett,
permeabilidades relativas, etc. El capítulo 3 estudia las propiedades PVT de los fluidos del
yacimiento, y se presentan diferentes correlaciónes y ecuaciones para la determinación de
las propiedades PVT del agua, gas y petróleo, además de presentar los factores que afectan
el recobro de crudo. El capítulo cuatro introduce los primeros conceptos de evaluación de
reservas de gas y petróleo mediante el método volumétrico, construcción de mapas
manualmente y asistidos por computador y determinación de volúmenes. El capítulo cinco
se enfoca a los métodos de balance de materia. Este es con seguridad el capítulo más largo
y quizá más importante del curso de Ingeniería de Yacimientos. En él, se desarrollará la
ecuación de balance de materia y se estudiarán sus múltiples aplicaciones, entre las que se
cuenta la linealización de dicha ecuación y el estudio de yacimientos de condensados, y
yacimientos naturalmente fracturados. Este capítulo, también hace referencia a las técnicas
de predicción del comportamiento y recobro final de un yacimiento de petróleo mediante
los métodos de Schilthuis, Tarner, Muskat, Pirson y Tracy, las cuales se presentan en
detalle en los apéndices. El capítulo sexto complementa la discusión capítulo del quinto
pero se enfoca porincipalmente a los cálculos de intrusión de agua a un yacimiento
DEDICATORIAS
Con mi más profundo amor dedico este trabajo a mi esposa Matilde Montealegre y a mis
adorados hijos Jennifer Andrea y Freddy Alonso Escobar. A mis padres Sotero Escobar
(QEPD) y Delfina Macualo Vda. De Escobar (QEPD), a mis hermanos Sotero Alonso
(QEPD), Dayra Stella y Leonardo Fabio Escobar. A mis sobrinos Daniel Alfredo Escobar,
Samuel Alejandro Escobar , Juan David Betancourt y Gabriel Fernando Betancourt. A mis
ahijados Raul Alejandro Vageón López, Claudia Patricia Montealegre Molina y José
Gabriel Rancel López.
AGRADECIMIENTOS
De igual forma, deseo expresar mi gratitud al Dr. José Gregorio Osorio Gallego, Profesor
de la Universidad Nacional de Colombia sede Medellín por el tiempo dedicado a la
evaluación del presente libro y a la formulación de valiosas recomendaciones.
PROLOGO
TABLA DE CONTENIDO
INTRODUCCION.................................................................................................................. 3
TABLA DE CONTENIDO .................................................................................................... 8
CAPITULO 1 ....................................................................................................................... 12
CONCEPTOS FUNDAMENTALES, DEFINICION Y CLASIFICACION....................... 12
1. ASPECTOS GEOLOGICOS............................................................................................ 12
1.1. FUENTE DE HIDROCARBUROS .............................................................................. 12
1.2. MIGRACION DE HIDROCARBUROS....................................................................... 13
1.3. DEFINICIONES DE INGENIERIA DE YACIMIENTOS .......................................... 13
1.4. DEFINICION DE YACIMIENTO................................................................................ 13
1.4.1. Clasificación Geológica de los Yacimientos.............................................................. 13
1.4.2. Clasificación de Acuerdo al Punto de Burbuja........................................................... 15
1.4.3. Clasificación de Acuerdo al Estado de los Fluidos .................................................... 16
1.4.4. Clasificación de Acuerdo al Mecanismo de Producción............................................ 25
1.4.5. Clasificación de Acuerdo a Variaciones del Volumen Originalmente Disponible a
Hidrocarburos ....................................................................................................................... 27
1.6. REGIMENES DE PRESION DE FLUIDOS ................................................................ 28
1.7. APLICACIONES DE REPEAT FORMATION TESTERS, RFT ................................ 42
1.8. ESCALAS DEL YACIMIENTO .................................................................................. 44
CAPITULO 2 ....................................................................................................................... 48
PROPIEDADES FISICAS DEL MEDIO POROSO............................................................ 48
2.1. POROSIDAD ................................................................................................................ 48
2.1.1. Clasificación Ingenieril de la porosidad .................................................................... 48
2.1.1.1. Porosidad absoluta................................................................................................... 48
2.1.1.2. Porosidad efectiva ................................................................................................... 48
2.1.1.3. Porosidad no efectiva .............................................................................................. 48
2.1.2. Clasificación Geológica de la porosidad .................................................................... 48
2.1.2.1. Porosidad primaria o intergranular.......................................................................... 49
2.1.2.1.1. Porosidad intercristalina ....................................................................................... 49
2.1.2.1.3. Planos estratificados ............................................................................................. 49
2.1.2.1.4. Espacios Sedimentarios Misceláneos .................................................................. 49
2.1.2.2. Porosidad secundaria, inducida o vugular ............................................................... 49
2.1.2.2.1. Porosidad de disolución........................................................................................ 49
2.1.2.2.2. Dolomitización ..................................................................................................... 50
2.1.2.2.3. Porosidad de Fractura ........................................................................................... 50
2.1.2.2.4. Espacios secundarios misceláneos ....................................................................... 50
2.1.3. Factores que afectan la porosidad............................................................................... 50
2.1.3.1. Tipo de empaque ..................................................................................................... 50
2.1.3.2. Material cementante ................................................................................................ 53
2.1.3.3. Geometría y distribución de granos........................................................................ 54
2.1.3.4. Presión de las capas suprayacentes......................................................................... 54
2.1.3.5. Presencia de partículas finas................................................................................... 54
2.1.4. Promedio de la porosidad .......................................................................................... 54
2.1.4.1. Promedio aritmético ................................................................................................ 54
2.1.4.2. Promedio ponderado................................................................................................ 54
CAPITULO 1
CONCEPTOS FUNDAMENTALES, DEFINICIÓN Y CLASIFICACIÓN
DE LOS YACIMIENTOS
1. ASPECTOS GEOLOGICOS
Los escudos continentales en áreas de aguas tranquilas, como lagos, cuencas profundas y
pendientes continentales, poseen las condiciones favorables para la depositación de los
sedimentos ricos en materia orgánica. Las tres etapas principales para la evolución de la
materia orgánica son diagénesis, catagénesis y metagénesis.
Puesto que la mayoría de las rocas en el subsuelo se hallan saturadas con agua, el
movimiento de hidrocarburos tiene que ser debido a flujo activo de agua, o flujo
independiente de la fase acuosa, por desplazamiento o por difusión. Puesto que
normalmente los hidrocarburos son menos densos que el agua, su acumulación toma lugar
en la parte más alta de la trampa. La migración es detenida por rocas superiores
relativamente impermeables.
Se entiende por yacimiento una unidad geológica de volumen limitado, poroso y permeable
que contiene hidrocarburos en estado líquido y/o gaseoso. Los cinco ingredientes básicos
que deben estar presentes para tener un yacimiento de hidrocarburos son: (1) fuente, (2)
Camino migratorio, (3) Trampa, (4) Almacenaje/porosidad, (5) Transmisibilidad/
Permeabilidad.
ESTRATIGRAFICO
GAS
PETROLEO
ESTRUCTURAL
PETROLEO
AGUA
COMBINADO
PETROLEO
FALLA
AGUA
3. Combinados: Hace referencia a las posibles combinaciones que se presenten entre los
dos grupos anteriores. Ver Fig. 1.1.b.
Pto Crit.
Aceite
A E
Subsaturado
Pb
Pto
B
roc
io
F
Presion
Aceite + Gas
C Saturado
Temperatura
PB
PRESION
PC
PC
Petróleo
Gas
PB
PR
PR
TEMPERATURA
Fig. 1.2.b. Diagrama de fases de la capa de gas y la zona de crudo mostrando capa de gas
retrógrada2
PC
Gas
PRESION
PB
PB
Petróleo
PR
PC
PR
TEMPERATURA
Fig. 1.2.c. Diagrama de fases de la capa de gas y la zona de crudo mostrando capa de gas
no retrógrada2
1. Petróleo negro. Consiste de una amplia variedad de especies químicas que incluyen
moléculas grandes, pesadas y no volátiles. El punto crítico está localizado hacia la
pendiente de la curva. Las líneas (iso-volumétricas o de calidad) están uniformemente
espaciadas y tienen un rango de temperatura amplio. Los primeros crudos de este tipo
fueron de color negro, de allí su nombre. También se le llama crudo de bajo encogimiento
o crudo ordinario. Estos crudos tienen GOR ≤ 1000 pcs/STB, el cual se incrementa por
debajo del punto de burbuja. Bo ≤ 2 y API ≤ 45 y el contenido de C7+ mayor o igual a 30
%, ver Fig. 1.3.c. Las temperaturas del yacimiento son menores de 250 °F. La gravedad
decrece lentamente con el tiempo hasta bien avanzada la vida del yacimiento donde vuelve
a incrementarse ligeramente. Este crudo es normalmente negro (compuestos pesados)
aunque pude ser marrón o verduzco. (Figs. 1.3.a, 1.3.b, 1.3.c y 1.7.b)5.
PETROLEO NEGRO
2,500
Condiciones
iniciales
Presió
2,000 Punto Crítico n de
Rocío
1,500 90 80
ja 70
bu 60
1,000 Bur 50
e
nd 40
e sió
Pr
30
20
500
Separador
0
0 50 100 150 200 250 300 350 400 450
PETROLEO NEGRO
1800
rocío
80 %
2
1000 70 %
os de
% de líquido
60 %
Pu n t
800
50 %
600 40 %
30 %
400
20 %
200 10 %
Separador
0
60 110 160 210 260 310
Temperatura, F
100000
Gases secos
90000
80000
Producida Inicial, scf/STB
Relación Gas-Petróleo
70000
60000
50000
40000
30000
20000
10000
Aceites negros
0
0 10 20 30 40 50 60 70
PETROLEO VOLATIL
4,000
3,500 Presió
n de
Condiciones Pto Crítico Rocío
3,000 iniciales
2,500
70 60
2,000
50 40 % líquido
30
a
uj
rb
20
Bu
1,500
de
n
ió
10
es
Pr
1,000
5
500 ocío
Separador
ión de R
Pres
0
0 100 200 300 400 500 600
PETROLEO
3500 NEGRO Puntos d
e rocío
Punto crítico
3000
2 PETROLEO
VOLATIL
Presión, psia
2500 60 %
GAS
RETROGRADO
2000 50 %
Punto crítico
1500 20 % 40 %
30 %
15 %
1000
20 %
Separador 10 %
10 %
500 5%
cío
s de ro
Punto
0
0 100 200 300 400 500
Temperatura, F
20000 12500
Aceite STB/D GOR, scf/STB
16000 10000
8000 5000
4000 2500
0 0
0 5 10 15 20 25 30
Tiempo, meses
GAS RETROGRADO
3,000
Condiciones Presió
iniciales n de
Rocío
2,500
2,000
Pto Crítico
1,500
% líquido
a
uj 30
rb 40 20
Bu 15
1,000 de
n
ió
es 10
Pr
5
500
Separador
0
0 50 100 150 200 250 300 350 400
1600 60000
Aceite STB/D GOR, scf/STB
1200 45000
800 30000
400 15000
0 0
0 200 400 600 800 1000 1200 1400
Tiempo, días
6000
Puntos de Rocio Puntos de burbuja
5000
Producida Inicial, scf/STB
Relación Gas-Petróleo
4000
3000
2000
1000
0
7.5 10 12.5 15 17.5 20 22.5
transparente. GOR > 15000 pcs/STB y permanece constante durante toda la vida del
yacimiento. Se producen menos de 60 STB crudo por cada millón de pies cúbicos normales
de gas. Ver Fig. 1.6.a y 1.6.b.
GAS HUMEDO
3,500
Condiciones
iniciales
3,000
o
ocí
2,500
eR
nd
sió
2,000
Pre
1,500 % líquido
uja
Pto
rb
1,000 5
Bu
Crítico
de
1
n
50
sió
25
e
500
Pr
Separador
0
0 100 200 300 400 500 600 700
50000
40000
Producida Inicial, scf/STB
Relación Gas-Petróleo
30000
20000
10000
0
0 2 4 6 8 10 12
1.6.b. Variación en el porcentaje de C7+ de acuerdo a la relación gas petróleo inicial (gas
húmedo)
GAS SECO
4,500
Condiciones
4,000 iniciales
ío
Roc
3,500
e
ión d
3,000
Pres
2,500
% líquido
2,000
25
1,500 Pto 1
Crítico
1,000 50
500
Separador
0
0 100 200 300 400 500 600 700
5. Gas seco. Está formado principalmente por metano y algunos intermedios. El diagrama
de fases muestra una mezcla de hidrocarburos gaseosa tanto en superficie como en el
yacimiento. No hay presencia de líquidos ni en yacimiento ni superficie. Sin embargo, a
temperaturas criogénicas, menores de 50 °F, se puede obtener luidos de estos gases. La
EBM puede aplicarse tanto a gas como gases húmedos para determinar gas original in-situ
y predecir reservas de gas. Ver Fig. 1.7.a.
6. Asfalténicos. En estos yacimientos, las condiciones iniciales del yacimiento están muy
por encima y a la izquierda del punto crítico. El rango de temperatura es bastante amplio.
Estos no se vaporizan ni tiene punto crítico. Ver Fig. 1.7.b.
Cuando la presión del yacimiento localiza a éste en la zona de una sola fase, normalmente
la composición se mantiene constante. Sin embargo, cuando la presión localiza al
yacimiento por debajo de la envolvente, el gas puede producirse más fácilmente y la
relación gas-petróleo no se mantiene causando un consecuente cambio del punto de burbuja
y un desplazamiento del diagrama de fases. Caso típico ocurre en un yacimiento de
condensado retrógrado. Al pasar por el punto de rocío la condensación toma lugar y el
líquido queda atrapado dentro de los poros de la roca, luego el producto en superficie
tendrá menos contenido líquido y como consecuencia se incrementa el GOR. Una vez el
punto de rocío se alcanza, la composición del fluido cambia y la envolvente se desplaza a
la derecha lo que agrava la pérdida de líquido en los poros.
Las Figs. 1.7.b, 1.7.c, 1.8 y 1.9 también ilustran la carecterización de los fluidos del
yacimiento.
Pyto, Tyto
Pto Crit.
Pto Crit.
gas seco
Pto Crit.
gas humedo
Pto Crit.
Presion
Gas condensado
Petroleo volatil
Pto Crit.
Psep, Tsep
Temperatura
Condensado Volátil 2
5000
Presión, psia
4000
3000
2000
Gas
húmedo
Aceite
negro
1000
Gas
seco
0
-200 -100 0 100 200 300 400 500 600 700 800
Temperatura, F
30000
20000
húmedo
10000
Gas
0
0 5 10 15 20 25 30
Fig. 1.8. Variación en el porcentaje de C7+ de acuerdo a la relación gas petróleo inicial
(todos los fluidos)
GOR
GOR
GOR
No hay
GOR
líquido
No hay
API
API
API
API
API
líquido
Fig. 1.9. Comportamiento del GOR y la gravedad API a través del tiempo
100
1 Expansion Roca y fluido
2 Empuje por gas disuelto
3 Expansion capa de gas
80 4 Intrusion de agua
5 Drenaje gravitacional
60
4
3
40
2
5
1
20
0
0 10 20 30 40 50 60
Eficiencia de recobro, % OOIP
OP =FP + GP
d(FP)=-d(GP)
⎛ dP ⎞
Pw = ⎜ ⎟ D + 14.7
⎝ dD ⎠agua
Esta ecuación asume que hay continuidad desde la superficie al fondo luego la salinidad es
constante. El gradiente del agua dulce es 0.4335 psi/ft. Para un caso anormal de agua:
⎛ dP ⎞
Pw = ⎜ ⎟ D + 14.7 + C
⎝ dD ⎠ agua
Po = 0.35 D + C
Presion, psi
14.7
Profundidad, ft
Presion Presion
de fluido de grano
FP GP
Presion de
sobrecarga, OP
Sobrepresion
5000
5100
Gas
Po GOC
5200
sib
l e
gr
ad
ie
5300 nt
e
no
rm Oil
5400 al
WOC
5500
5600
Agua
En el GOC a 5200 ft, Pg = Po, luego Po = 0.35 D + 565 = 0.35(5200) + 565 = 2385 psi
Puesto que Pg = Po y a 5200 ft, 2385 psi, se tiene 2385 = 0.08(5200) + C de donde C =
1969, luego Pg = 0.08D + 1969 .
Pg = 0.08 D + 1969 = 0.08(5000) + 1969 = 2369 psi. Si se asume que se está perforando una
formación con agua dulce, a la profundidad de 5000 ft (tope) la presión será de:
Lo que equivale a que se tendría una patada de 2369 – 2266 = 103 psi al entrar a la zona de
petróleo.
EJEMPLO
Un pozo penetra la arena First Bromide a una prefundida de 7500 ft. El pozo cruza el
contacto gas-petróleo a 7510 ft y el contacto agua-petróleo a 7550 ft. La formación está
normalmente presionada. Qué presiones se esperan en los contactos y en el tope de la
formación. Asuma que los gradientes de agua, gas y crudo son 0.45, 0.35 y 0.08 psi/ft,
respectivamente.
SOLUCIÓN
Pw = 0.45D + 14.7
Po = 0.45D + C
3412.2 = 0.45(7550) + C de donde C = 769.7 psi. En el contacto gas petróleo (D=7510 ft):
WOC 7550'
Una milla
GOC 7510'
Po = 0.35D + 769.7
Pg = 0.08D + C
EJEMPLO
Un pozo perforado en la ten Mile Falt encontró la formación Red Fork probó positivo para
producción de crudo. El operador midió una presión de 3880 psia en el fondo del pozo a
una profundidad de 8600 ft. También sobre la formación Red Fork a una milla de distancia
se midió una presión de 4050 psia en un pozo seco a una profundidad de 9000 ft.
SOLUCIÓN
De acuerdo con la Fig. 1.14 a la profundidad de 9000 ft, la ecuación gobernante (agua) está
dada por:
⎛ dG ⎞
Pw = ⎜ w ⎟ D
⎝ dP ⎠
⎛ dG ⎞
4050 = ⎜ w ⎟ 9000
⎝ dP ⎠
De donde resulta un gradiente de 0.451 psi/ft por lo que se considera que la formación está
normalmente presionada.
Se sabe que en el WOC, las presiones de agua y petróleo (Po = Pw) tienen el mismo valor
para que halla equilibrio:
Po = 0.35D
Pw = 0.45 D
Se sabe, además, que la presión a 8600 ft es de 3880 psia. Luego en el contacto WOC, se
tiene:
EJEMPLO
De acuerdo con la Fig. 1.15, el pozo A está siendo perforado con un lodo base agua cuya
densidad es 9 ppg y en la actualidad se encuentra a una profundidad de 4000 pies con
respecto al nivel del terreno. Seis meses antes, se terminó de perforar el pozo B que resultó
seco. El operador prudentemente midió la presión de fondo de este pozo correspondiente a
3100 psia. Con base en la información suministrada en la Fig. 1.15 determine las presiones
en cada uno de los contactos y la densidad mínima para perforar la formación gasífera.
Asuma que el gradiente del crudo es 0.35 psi/ft y el del gas es de 0.08 psi/ft.
SOLUCIÓN
⎛ dG ⎞
P=⎜ ⎟D
⎝ dP ⎠
⎛ dG ⎞
3100 = ⎜ ⎟ 6000
⎝ dP ⎠
⎛ dG ⎞
⎜ ⎟ = 0.517 psi / ft
⎝ dP ⎠
Con este valor de gradiente se determina la presión en el contacto agua petróleo (D=5500
ft):
200 millas
2000 ft B A
Nivel del terreno
Nivel del
1000 ft
mar 4900 ft
Fo 5000 ft
r m
ac 5500 ft
ió GOC
n
WOC
6000 ft
A
Presión de 3100 psia fue
medida en el punto A
⎛ dG ⎞
P5500 = ⎜ ⎟D
⎝ dP ⎠
Entonces:
⎛ dG ⎞
Po = ⎜ o ⎟ D + C
⎝ dP ⎠
2843.3 = 0.35(5500) + C
De donde C = 918 psi. En el contacto gas petróleo (D=5000 ft), la presión será:
La ecuación en la zona de gas inmediatamente después de la zona del contacto está dada
por:
Pg = 0.08D + C
De donde:
2668 = 0.08(5000) + C
Puesto que:
Pg = 0.052 ρ h
El tema del gradiente de presión toma auge con la aparición de herramientas, tales como el
RFT y el aumento de herramientas similares, como el MDT. El primero ha estado en
operación desde mediados de 1970 y ha revolucionado la relación presión-profundidad en
columnas de hidrocarburos y acuíferos. La implicación de la ingeniería a de yacimientos
en este tema está en la ubicación de los fluidos en contacto con la formación que permite
calcular el volumen neto de de roca que aparece en la ecuaciones volumétricas (Capítulo
4). La situación representada en la Fig. 1.16a muestra una estructura masiva de
cabalgamiento donde todos los pozos penetran directamente el contacto agua-aceite (OWC)
el cual será detectado sobre corazones y registros. Las capas de los yacimientos con lutitas
impermeables (Fig. 1.16b), sin embargo, representa un desafío para la ingeniería en cuanto
al establecimiento de los contactos entre fluidos. En esta complicada, las presiones son
controladas por la presión común del acuífero, pero el contenido de hidrocarburos en
cualquiera de las capas individuales es dictado por trayectorias de migración, propiedades
de la roca del yacimiento, etc.
Vsc
ρ gr = ρ gsc = 0.0763γ g E (lb/ft3)
Vr
(b)
Pozo gas
(a) Pozo
gas cru
do
gas
a
ag gua
cru ua
OWC gas do
cru
do ag
ua
Falla
ag
ua
Fig. 1.16. (a) Estructura de cabalgamiento masiva. (b) sistema multi-capa
Presión Pozo 1
Gas 0.08 psi/ft
Pozo 2
Profundidad
GOC
GOC
Petróleo 0.33 psi/ft
dpg 0.0763γ g E
= (psi/ft)
dD 144
conservación de la masa usando densidades de superficie del crudo y del gas juntos con las
propiedades PVT3.
El gradiente de presión de crudo puede ser calculado como ρor/144 psi/ft. Con la
combinación de un solo punto de presión y el gradiente, es posible realizar las líneas de
presión-profundidad, las cuales se interceptan con las líneas de gas localizadas en el
contacto crudo-gas3.
Po − Pb
∆D = (ft)
dp / dD
Si el valor calculado de ∆D localiza el GOC dentro del yacimiento, entonces, puede haber
una capa de gas libre, pero esto no es seguro. En la ecuación dada anteriormente se asume
que el gradiente de presión en el crudo, dP/dD es constante, pero en algunos casos,
especialmente en yacimientos con suficiente espesor, las propiedades PVT y por tanto, el
gradiente varían con la profundidad, lo cual distorsiona los cálculos del GOC. La única
manera segura de encontrar la presencia de capa de gas es perforando un pozo en la cresta
del yacimiento. La Fig. 2.19b ilustra una incertidumbre similar asociada con la estimación
de campos de gas. El gas ha sido visto únicamente debajo el nivel bajo de gas (gas-down-
to, GDT), pero esto permite la posibilidad de una zona de crudo, como se verá en el
siguiente ejercicio.
Gas? GOC
Pozo
Profundidad
Petróleo Petróleo
OWC
Agu
a
(b) Presión
Pozo
Gas
Profundidad
GDT
OWC?
Posible Agu WOC más
GWC a
profundo
EJEMPLO
Los campos de gas costa afuera se evaluaron con pozos exploratorios, A1 y después con
dos pozos A2 y A3. Estos pozos penetraron dos yacimientos delgados que producen gas, X
e Y, como muestra la Fig. 1.20 y una secuencia de yacimientos más profundos de petróleo
(no mostrados). Los resultados de DSTs en los tres pozos se muestran en la tabla 1.2. En las
pruebas realizadas en areniscas productoras de petróleo, debajo de los yacimientos de gas,
se encontró que estaban sobrepresionados debido al incremento encontrado con la
profundidad, el cual, es atribuido a un extenso intervalo de lutitas entre cada yacimiento
que cambia el equilibrio hidrostático con el enterramiento. Han sido evaluados
adecuadamente estos yacimientos con los tres pozos para permitir el desarrollo del campo,
en caso contrario, cuál es la mejor estimación requerida3?
SOLUCIÓN
P 2797
E = 35.37 ( scf / ft 3 ) = 35.37 = 154.87
zT 0.91(242 + 460)
Luego:
A2 A3
A1
SP inducción SP inducción
SP inducción
5900`
5300` 6900`
DST 9` X
DST 8`
5100`
5500`
7100`
DST 7`
5300` Y
DST 4`
5700`
7300`
6000 A2/8
A2/7
GWC-6430´
6500
Si esto es solamente gas, el contacto gas-agua (GWC) puede estar a 6430 ft.ss, pero allí,
también es posible encontrar subyaciendo una acumulación de crudo. La máxima
extensión de éste se indica por la línea interrumpida en la Fig. 1.22. Esto se da por el hecho
que el nivel por encima del agua (water-up-to, WUT) en el yacimiento X está a 6928 ft.ss,
pero es posible, en el caso más optimista que halla crudo inmediatamente sobre esta
profundidad. El gradiente de crudo en yacimientos más profundos es 0.3 psi/ft y usando
esta cifra, como muestra la Fig. 1.22, se puede llegar a un posible GOC a 6150 ft.ss
implicando una máxima columna3 de crudo subyacente de 780 ft.
6000
GOC-6150´
7500 Columna
máxima de
Profundidad petróleo,
(ft.ss) 8000 ft.ss)
8000
En el yacimiento más profundo, Y, hay dos pruebas en la zona de gas. Podría resultar
incorrecto unir las presiones con una línea recta porque como puede verse en la tabla 1.2 el
gradiente de presión resultante puede ser (3112-3100)/(6212-5544) = 0.018 psi/ft el cual es
un valor físicamente irreal. Las dos pruebas de gas en este yacimiento fueron efectuadas
usando diferentes registradores de presión y por lo tanto no hay validez en la conexión de
los puntos. Considerando que el gradiente promedio para las dos pruebas calculado con el
PVT es 0.061 psi/ft,
Po − Pg
∆D = (ft)
dPo / dD − dPg / dD
En el cual la presión está evaluada a 8000 ft.ss (Fig. 1.22) para Po = 3512 psia (al posible
OWC) y Pg = 3240 psia (por extrapolación) y usando el gradiente de petróleo y gas de 0.3 y
0.061 psi/ft, entonces el ∆D = 1138 ft el cual localiza el GOC más profundo a 6862 ft.ss
3512 − 3240
∆D = = 1138 ft
0.3 − 0.061
El RTF fue introducido a mediados de los 70s. Su mayor ventaja sobre su antecesor el FIT
(prueba de intervalos de formación), fue que este puede medir un número ilimitado de
puntos de presión en un sólo viaje al pozo mientras que el FIT se restringió a uno.
Originalmente se consideró que la aplicación más importante del RFT era para muestreo de
fluidos, pero después se observó su eficacia para proporcionar valores de presión-
"B"
5800 ft.ss
5900
6000
6100
6200
OWC 6100 "A"
6300
Fig. 1.23. El mapa de contorno estructural muestra las locaciones de los pozos "A" y "B"3
Durante la etapa evaluativa de un campo, las lecturas del RTF proveen mayor calidad de
datos de presión y se corre continuamente para establecer el contacto entre fluidos. A
manera de ilustración del uso de RTF valoración de un campo, el mapa estructural de la
Fig. 1.23 muestra que la capa que contiene crudo tiene 50 ft de espesor y un acuífero
establecido por otros medios. El primer pozo perforado, A, fue realizado en el acuífero,
pero el operador había previsto medir la presión de agua sobre un intervalo de 160 ft para
hallar la presión de agua en esta nueva área. El segundo pozo, B, perforado varios
kilómetros al este fue más afortunado, se descubrió un yacimiento de 50 ft de espesor
productor de crudo con buena porosidad y permeabilidad. Se midieron seis presiones RTF a
través de este intervalo. Las presiones registradas en los dos pozos se muestran a
continuación.
"B" columna
5800 de crudo
28 psi/ft
6000
"A" Agua
Profundidad,
(ft.ss)
6400
Debe notarse que todas las profundidades deben ser convertidas a valores bajo el nivel del
mar para facilitar la comparación directa de presiones. Igualmente, las presiones deben
estar en psia. Las medidas anteriores de presión y de tensión registraban en psig, lo cual es
relativo a la presión atmosférica: psia = psig + 14.7. Los ingenieros deben especificar
totalmente la presión usada. No es común encontrar presiones únicamente en unidades de
psi, ya que, no es una unida de presión absoluta.
La gravedad específica del crudo con gas disuelto a condiciones iniciales del yacimiento,
determinado en análisis PVT, fue 0.644 relativo al agua. La densidad del agua pura es
62.43 lb/ft3, el gradiente de presión de crudo es (0.646 * 62.43)/144 = 0.28 psi/ ft, tal línea
de gradiente se ajusta a través de la medida de puntos de presión y extrapolando esto da un
OWC de 1635 ft.ss: 422 ft bajo el oil-down-to (ODT) de 5813 ft.ss.
Una manera de que tan bien entendemos el yacimiento puede obtenerse considerando la
fracción del yacimiento que está siendo muestreada mediante las diferentes técnicas. Por
ejemplo, supongamos que se desea hallar el tamaño del área muestreada desde un pozo que
tiene un radio de 6 pulgadas. Si se asume un área circular, el área se puede estimar como
MICRO
MACRO
MEGA
GIGA
Una señal de un registro eléctrico expande el área siendo muestreada. Suponga que un
registro pueda penetrar la formación unos 5 pies desde el pozo, lo cual es razonable. La
fracción del área siendo muestreada es 4 partes en 10000. El tamaño muestreado dentro del
área de drenaje (5 acres) es todavía una fracción de un porcentaje.
Los corazones y registros eléctricos dan una visión muy limitada del yacimiento. Una
sección sísmica expande la fracción del área muestreada, pero la interpretación de datos
sísmicos es menos precisa. La credibilidad de la sísmica se puede mejorar
Correlaciónándola con datos de análisis de corazones o perfiles eléctricos.
La Fig. 1.16 presenta la definición de escala de yacimiento. Note que éstas no son
universalmente aceptadas, pero ilustran la escala relativa asociada con la propiedad del
yacimiento medida. La escala Giga incluye información asociada con geofísica, tales como
arquitectura del yacimiento. Esta también incluye teorías de caracterización regional como
REFERENCIAS
1. Smith, C.R., and Tracy, G.W. “Applied Reservoir Engineering”. Oil & Gas
Consultants, Inc. Tulsa, Ok. 1987.
2. Craft, B.C. and M.F., Hawkins. “Applied Reservoir Engineering”. Prentice-Hall
International. New Jersey, 1991.
3. Dake, L.P. “Fundamental of Reservoir Engineering”. Elsevier Scientific Publishing
Co. 1978.
4. Dake, L.P. “The Practice of Reservoir Engineering - Revised edition”. Elsevier
Developments in Petroleum Science. Second impresión. Amstrerdam, Holanda,
2004.
5. McCain, W. “The Properties of the Petroleum Fluids”. Gulf Publishing Co. 1988.
6. Abdus S. and Ganesh T. “Integrated Petroleum Reservoir Management: A Team
Approach”. PennWell Books. Tulsa, Ok. 1994.
CAPITULO 2
PROPIEDADES FÍSICAS DEL MEDIO POROSO
2.1. POROSIDAD
Vp
φ=
Vt
Vp = volumen poroso
Vt = volumen total
2.1.1.3. Porosidad no efectiva. Es la diferencia que existe entre las porosidad absoluta y
efectiva.
A medida que los sedimentos se depositaron en los mares antiguos, el agua fue el primer
fluido que llenó el espacio poroso. Esta agua se le denomina agua connata. Un método
común de clasificación de la porosidad se basa en la condición si porosidad se formó
inicialmente o si fue producto de una diagénesis subsiguiente (dolomitización),
catagénesis, campo de esfuerzos o percolación de agua.
2.1.2.1.1. Porosidad intercristalina. Se refiere a los espacios existentes entre los planos
de un cristal o espacios vacíos entre cristales. Muchos de éstos poros son sub-capilares, v.g.
poros menores de 0.002 mm de diámetro. La porosidad que se encuentra entre cristales o
partículas tamaño lodo se llama comúnmente “microporosidad”.
2.1.2.1.2. Porosidad Integranular. Es función del espacio vacío entre granos, es decir, de
los espacios intersticiales de toda clase en todo tipo de roca. Esta porosidad comprende
tamaño sub-capilar a super-capilar. Generalmente, los espacios tienen un diámetro mayor
de 0.5 mm.
espacios vacíos causados por organismos vivientes pueden sufrir alargamiento debido a
dilución.
Algunas rocas carbonatas están constituidas solamente por calizas. Si el agua circulante a
través del espacio poroso contiene suficientes cantidades de magnesio disuelto, el calcio en
la roca puede intercambiarse por el magnesio en solución. Como el magnesio es
considerablemente más pequeño que el calcio, la resultante dolomita tendrá una porosidad
mayor, cuyo incremento oscila entre el 12-13 %.
90° 90°
90°
90°
90°
Empaquetamiento cúbico Empaquetamiento romboédrico
de esferas uniformes, φ = 47.64 % de esferas uniformes, φ = 25.96 %
2r
o
60
a
c
o
b 45
Vt = a ⋅ b ⋅ c
donde:
a = 4r cos 30o
b = 4r
c = 4r
Vt = 55.426r 3
⎛4 ⎞
Vgr = 8 ⎜ π r 3 ⎟
⎝3 ⎠
Vgr = 33.51r 3
Vt = a ⋅ b ⋅ c
a = 4r cos 45o
b = 4r
c = 4r
Vt = 45.25r 3
⎛4 ⎞
Vgr = 8 ⎜ π r 3 ⎟
⎝3 ⎠
Vgr = 33.51r 3
2.1.3.2. Grado de cementación o consolidación. Cemento que une los granos y que se
forma posterior a la depositación ya sea por dilución de los mismos granos o por transporte.
Son cementos: carbonato de calcio, carbonato de magnesio, carbonato de hierro, limonita,
hematita, dolomita de sulfato de calcio, arcillas, y combinación de estos. Las areniscas
altamente cementadas presentan bajas porosidades. Lo contrario ocurre con rocas poco
consolidadas. La cementación toma lugar tanto en el tiempo de mitificación como en el
proceso de alteración de la rocas causada por agua circulante.
2.1.3.4. Presión de las capas suprayacentes. Las capas suprayacentes pueden compactar
el yacimiento y reducir el espacio poroso. La compactación tiende a cerrar los espacios
vacíos, forzar el fluido a salir y permitir un mayor acercamiento de las partículas minerales,
especialmente en rocas sedimentarias de grano fino.
∑
n
φ
φ= i =1 i
∑
n
φx
φ= i =1 i i
∑
n
x
i =1 i
φ = n φ1φ2φ3 ...φn
Es imposible tener una descripción detallada del espacio poroso pero, con respecto a los
granos en un material granular, es posible contar con una descripción estadística de una
clase u otra. La determinación de la distribución del tamaño de poro es de importancia
especial en materiales consolidados donde no puede obtenerse una distribución del tamaño
de grano. Puesto que los medios porosos poseen una estructura porosa más o menos
aleatoria, no es sorprendente que las muestras pequeñas del mismo material no tengan la
misma porosidad o permeabilidad. Generalmente se observa que a mayor volumen, para
muestras individuales, hay más probabilidad de que los mismos valores de permeabilidad y
porosidad se observen. Estas características de los materiales porosos pueden entenderse
por el siguiente análisis. Considere un volumen bruto de material poroso e imagine que este
se divide en paralelepípedos muy pequeños. Estos elementos poseen una distribución de
porosidad debido a la estructura aleatoria del medio poroso. Denótese esta función de
distribución como F(φ), de modo que F(φ)dφ es la fracción de los elementos que tienen
porosidades entre φ y φ + dφ. La porosidad media de estos elementos, y la porosidad real
del material bruto es:
1
φ = ∫ φ F (φ )dφ
0
Suponga que se toman varias muestras del material poroso. La porosidad promedio de la
muestra puede definirse como:
1 n
φ = ∑ φi
n i +1
ln F(φ) B
(φ − φ ) 2
Fig. 2.2. Función de distribución de porosidad
(φ −φ )
2
−
1 2 (σ / n )
F (φ ) = e
⎛ σ ⎞
2π ⎜ ⎟
⎝ n⎠
30
25
20
F(φ)
15
10
0
0-5 5 - 10 10 - 15 15 - 20 20 - 25 25 - 30
Rango de porosidad, %
2
F (φ ) = Ae − B (φ −φ )
ln F (φ ) = ln A − B(φ − φ ) 2
La muestra de roca de la Fig. 2.3 presenta la distribución gausiana mostrada en la Fig. 2.4.
Es la relación que expresa la cantidad de fluido que satura el medio poroso. Conocida dicha
cantidad y la extensión del volumen poroso se puede volumétricamente determinar cuanto
fluido existe en una roca.
unidad de masa retirada está siendo reemplazada por una misma cantidad que se adiciona al
sistema. Este toma lugar en yacimientos con empuje de agua o capa de gas. Además, puede
aplicarse sin un margen de error significativo en las zonas aledañas al pozo. El flujo
inestable presenta variaciones de la presión con el tiempo y el flujo pseudoestable (falso
estable o semiestable) es un flujo inestable que se puede considerar temporalmente estable.
dP
≠0
dt dP
Presión
Presión
=0
dt
Tiempo Tiempo
Steady state
q at t=0 to t=inf
Caudal
Pressure
p at t=0 to t=inf
Radius
p a t0 Unsteady state
p a t1 p a t2
p a t3
Caudal
Presión
q a t2
Pwf q a t3
q a t1
r1
Radio
p a t0
p a t3
p a t1 p a t2 p a t4
p a t5
p a t6
Presión
Caudal
q a t2, t3, t4 y t5
Pwf q a t6
q a t1
r1
Radio
1 dV
c=−
V dP
luego,
dV = −cVdP
Para que haya expansión tiene que haber una caída de presión. Mientras que la presión no
afecte la frontera el comportamiento es infinito. Todos los yacimientos tienen ese
comportamiento. En el estado pseudoestable el caudal en el fondo del pozo se mantiene
constante lo que nindica que para que esto ocurra tienen que haber líneas de presión
paralelas como lo indica la Fig. 2.6.c. El principio es similar al estado inestable, pero
cuando la presión afecta las fronteras, en todo punto del yacimiento dP/dt es el mismo y se
obtienen líneas paralelas. Matemáticamente, dP/dt es constante y entonces la rata será
constante hasta que la presión de fondo fluyente no se pueda mantener.
Existen otras clasificaciones de los estados de flujo de acuerdo con la geometría que
presenten las líneas isobáricas (flujo radial, lineal, esférico, etc.). Las figs. 2.7.a y 2.7.b
describen las geometrías de flujo que se presentan más frecuentemente en yacimientos de
hidrocarburos.
Cuando el pozo se completa en una pequeña parte en el interior del estrato, la geometría de
una esfera puede tomar lugar. Este se le llama flujo esférico.
Flujo radial
Pozo
LENTE DE ARENA
VISTA DE PERFIL
LENTE DE ARENA
Flujo radial VISTA DE PLANTA
FLUJO HEMISFERICO
Pozo VISTA DE PERFIL
FLUJO HEMISFERICO
VISTA DE PLANTA
Pozo
Slider1 dice que si tratáramos de resolver todos los problemas de Ingeniería utilizando la
ley de Darcy es comparada como “ir a la luna en caballo”. Por lo tanto, la Ley de Darcy
puede aplicarse a condiciones muy particulares.
k
q = −cA ⎡ ∆P + ρ gh senθ ⎤⎦
µ⎣
Donde c = 1.127 para conversión. En unidades de campo esta ecuación se escribe como:
k ⎛ dp ⎞
v = −0.00127 ⎜ − 0.433 γ cos θ ⎟
µ ⎝ ds ⎠
k dp flujo
v = −0.00127
µ ds
∆Pflujo = ∆P ± 0.433γ∆D
FLUJO LINEAL
k ∆P
v=−
µ ∆x
. kA( P1 − P2 )
1127
qo =
µL
El flujo gaseoso difiere del flujo líquido ya que la rata de flujo volumétrico, q, varía con la
presión. Es decir, un pie3/día a 1000 psi no representa la misma masa de gas que el mismo
sistema a 100 psi. Para aliviar el anterior problema se introduce la ley de los gases reales
pV = znRT
Si se emplean psi, pies cúbicos, grados Ranking, la constante universal de los gases, R, es
10.73. Si un mol de gas a condiciones normales (14.7 psi y 60 °R) ocupa 379 ft3, se puede
tener la siguiente expresión1:
T
V / scf = z (1/ 379 ) (10.73)
P
Si se multiplica esta expresión por la rata de flujo expresada en miles de pie3/día, qg, y se
hace la conversión a barriles se tiene:
T ⎛ 1000 ⎞
V / scf = z (1/ 379 ) (10.73)qg ⎜ ⎟
P ⎝ 6.615 ⎠
5.04qg zT
q=
P
1.127kA ∆P
q =−
µ ∆x
se tiene:
5.04qg zT 1.127 kA ∆P
=−
P µ ∆x
Separando integrandos:
L P2
1.127 kA
5.04qg zT ∫ dx = − ∫ PdP
0
µ P1
1.127kA
5.04qg zTL = − ( P2 − P1 )
2µ
Resulta en:
FLUJO RADIAL
2π krh ∆P
q=
µ ∆r
Integrando:
Pe r
q µ e dr
∫ dP =
Pw
2π kh r∫w r
Integrando y despejando q;
7.08kh ( Pe − Pw )
q=
µ ln(re / rw )
7.08kh ( P1 − P2 )
qo =
µ ln(r1 / r2 )
En una forma similar al caso lineal gaseoso, igualando el flujo volumétrico de gas con el
flujo en geometría radial, se tiene:
Separando integrandos:
re Pe
7.08kh
5.04qg zT ∫ rdr = ∫ PdP
r
µ Pw
integrando:
re 7.08kh
5.04qg zT ln = ( Pe − Pw )
rw 2µ
De donde:
El término n es un parámetro que se introduce para responder por los efectos de flujo
rápido que toma lugar cada vez que el flujo se aproxima al pozo. Este término se evalúa
mediante pruebas de presión y para efectos prácticos podría asumirse igual a la unidad.
FLUJO HEMISFÉRICO
El flujo hemisférico, ver Fig. 2.7.b, para fluidos compresible e incompresibles puede
derivarse de la misma manera que se hizo para la geometría radial y lineal, con la salvedad
que para este caso la sección transversal de flujo es 2πr2 (media esfera). Con base en lo
anterior para flujo radial de crudo se tiene1:
1.127kA ∆P
qo =
µ ∆r
7.0811kr 2 ∆P
qo =
µ ∆r
7.08kh ( P1 − P2 )
qo =
µ (1/ r1 − 1/ r2 )
El radio interno r1 no es el radio del pozo. Es una función del radio del pozo y el intervalo
perforado. Una aproximación a éste se puede obtener fijando el área seccional de flujo en el
pozo igual a un área seccional hipotética:
2π rw h = 2π (r1 hemisferico ) 2
r1 hemisferico = rw h
En estas ecuaciones
n = constante de turbulencia
qo = bbl/día medidos a condiciones de yacimiento
Contacto agua-petróleo
lla
Fa
o
θ = 135
Pozo
Falla
lla
r e2
Fa
θ = 75 o
θ = 60
o
Pozo
Falla
re1
60 7.08kh ( Pwoc − Pw ) 75 7.08kh ( Pwoc − Pw )
q= +
360 µ ln(re1 / rw ) 360 µ ln(re 2 / rw )
π re2
Area de drene = 2re L +
2 Flujo radial Falla
agua-petróleo
Contacto
Flujo lineal re
qµ L 2q µ ln(re / rw )
Pwe − Pw = +
1.127 kA 7.08kh
Fig. 2.7.e. Aproximación de las ecuaciones de flujo a dos sistemas en serie
qg = Mft3/D
A = Área en pie²
k = Permeabilidad, Darcy
P = Presión en psia
T = Temperatura del yacimiento en °R
Existen muchos casos donde no se presenta una geometría de flujo completamente definida
sino parcial o combinada con dos muy diferentes. En estos casos se puede tener sistemas en
serie y sistemas en paralelo1.
El primer caso se describe en la Fig. 2.7.c donde al pozo no convergen líneas de flujo
alrededor de los 360°, sino que se cubre un área de 135°, por lo que el flujo efectivo hacia
el pozo resultaría siendo 135/360 = 0.375. Esto indica que el 37.5 % del círculo es
disponible al flujo por el ello se precisa modificar todas las ecuaciones de flujo mediante la
siguiente expresión:
θ 7.08kh ( P1 − P2 )
qo =
360 µ ln(r1 / r2 )
dar una aproximación del área total. El modelo resultante asume que las caídas de presión
en ambos subsistemas son aditivas. Al sumar la caída de presión desde el contacto de agua
hasta el radio de drene del pozo y del radio de drene del pozo hasta el pozo se puede tener
una expresión que permita obtener el caudal o la presión en algún punto.
Si el sistema se divide de forma tal que la combinación de las geometrías implica sistemas
en paralelo, entonces las ratas de flujo son aditivas. Ver Fig. 2.7.e. Esta geometría es
similar a la mostrada en la Fig. 2.7.d, sin embargo, las distancias a los contactos no son
uniformes. Por lo tanto una mejor aproximación es la combinar en dos sistemas radiales.
Cuando se perfora un pozo la invasión de los fluidos hacia la formación para controlar el
filtrado causa disminución en la permeabilidad en las zona aledañas al pozo. Entre otros
muchos factores, ésto hace que se genere una caída de presión adicional a la que se debería
obtener en condiciones normales. Además, la estimulación es un mecanismo usado
ampliamente en la industria petrolera para incrementar la productividad de un pozo. Lo
anterior ocasiona que la caída de presión esperada sea menor. Luego existe una caída de
presión adicional, pero favorable, en los alrededores del pozo. La Fig. 2.7.f clarifica este
fenómeno que toma lugar en una zona infinitesimal alrededor del pozo. En dicha figura la
zona de daño ha sido exagerada1,2,3.
La caída o ganancia de presión adicional causada por el daño o skin factor, s, está dada por:
141.2qµB
∆Ps = s
kh
Con base a lo anterior, la ley de Darcy incluyendo efectos de daño resulta ser:
7.08kh ( Pe − Pw − ∆Ps )
q=
µ ln(re / rw )
o también:
141.2q µ B
( Pe − Pw − s)
7.08kh kh
q=
µ ln(re / rw )
141.2qµ B re
Pe = Pw + ln
kh rw
Estimulación,
s<0
⎧
∆Ps ⎨
⎩
Daño,
⎧ s<0
∆Ps ⎨
⎩
Zona de
daño
Radio, ft
141.2qµ B re 141.2q µ B
Pe = Pw + ln + s
kh rw kh
Rearreglando:
141.2q µ B ⎛ re ⎞
Pe = Pw + ⎜ ln + s ⎟
kh ⎝ rw ⎠
7.08kh ( Pe − Pw )
q=
µ ⎛r ⎞
ln ⎜ e ⎟ + s
⎝ rw ⎠
Cuando un pozo está dañado, se tiene la impresión de que el radio del pozo tiene un valor
diferente, el cual es más grande que el real cuando el pozo está estimulado o mas pequeño
que el real cuando el pozo está dañado. Este radio aparente, rw’ se calcula mediante1,2:
rw ' = rwe− s
7.08kh ( Pe − Pw )
q=
µ ln(re / rw ')
⎛ k ⎞ r
s = ⎜⎜ − 1 ⎟⎟ ln s
⎝ ks ⎠ rw
El significado del las fracturas como medios para permitir el paso de fluidos se puede
evaluar considerando una fractura simple extendida cierta distancia dentro de la roca, ver
Fig. 2.8.a, usando la ecuación de hidrodinámica para flujo a través de placas paralelas6:
w3 h∆P
qf =
12 µ L
q w 2 ∆P
v= =
w ⋅ h 12 µ L
∆P
v = 9.869 × 10−9 k f
µL
P1 P2
wf
w
k f = 8.444 × 10−9 w2
v
vr =
φ f (1 − S wc )
Donde la porosidad de la fractura es la relación que hay entre su ancho por el ancho total de
la muestra. Bajo estas condiciones la permeabilidad de la fractura en Darcies es:
k f = 8.444 × 10−9 φ f (1 − S wc ) w2
Vo = Vom + Vof
Vof = Ahφ f (1 − S wf ) / Bo
VoR ≅ Vof FR f
⎛ n wh ⎞ ⎛ n f wh ⎞
kmf = ⎜ f ⎟ k f + ⎜1 − ⎟ km
⎝ A ⎠ ⎝ A ⎠
EJEMPLO
SOLUCIÓN
1.127kA( P1 − P2 ) 1.127(1)(1)(10)
qo = = = 0.0075 bbl / dia
µL 1.5(1)
P2 P1
canal q
rc
w
L
⎛ n f wh ⎞ ⎛ n f wh ⎞
kmf = ⎜ ⎟ k f + ⎜1 − ⎟
⎝ A ⎠ ⎝ A ⎠
π rc4 ∆P
q=
8 µL
∆P
q = 9.369 ×10−9 π rc2 kc
µL
Ar = φcπ rc2 (1 − S wi )
⎛ n π r2 ⎞ ⎛ n π r2 ⎞
kmc = ⎜ c c ⎟ kc + ⎜ 1 − c c ⎟ km
⎝ A ⎠ ⎝ A ⎠
EJEMPLO
SOLUCIÓN
kc = 12.6 ×106 φc rc2 (1 − S wi ) = 12.6 × 106 (0.03)(0.052 )(1 − 0.18) = 775 Darcy
⎛ n π r2 ⎞ ⎛ n π r2 ⎞
kmc = ⎜ c c ⎟ kc + ⎜ 1 − c c ⎟ km
⎝ A ⎠ ⎝ A ⎠
⎛ 5π (0.052 ) ⎞ ⎛ 5π (0.052 ) ⎞
kmc =⎜ 2 ⎟ 775000 + ⎜ 1 − ⎟ (1) = 36.5 mD
⎝ 929 ⎠ ⎝ 9292 ⎠
La ley de Darcy como tal considera que un solo fluido satura 100 % el medio poroso, por lo
tanto, el estado estable prevalece. Otra consideración hecha por Darcy es que el flujo es
homogéneo y laminar. La ecuación de Forchheimer tiene en cuenta los factores inerciales
que determinan que el flujo no es laminar o no Darcy6.
dP µ v
− = + βρ v 2
dL k
4.851x104
β=
φ 5.5 k
2.6x1010
β=
k 1.2
Estando k en md.
Aunque este fenómeno no se presenta a menudo en campos petroleros, puede ser común en
los laboratorios, donde a bajas presiones la molécula de gas puede tener el mismo tamaño
que el de los poros por lo que no se presenta un perfil de flujo o no existe flujo viscoso.
Realmente, para el caso de los gases, el fluido no se pega a la pared de los poros como
requiere la ley de Darcy y un fenómeno llamado deslizamiento toma lugar. Este
escurrimiento o desplazamiento del fluido a lo largo de los poros sugiere una aparente
dependencia de la permeabilidad con la presión. La relación propuesta por Klinkenber entre
presión y permeabilidad es1,2:
⎛ b⎞
k = k∞ ⎜ 1 + ⎟
⎝ P⎠
q = q1 + q2 + q3
k h + k 2 h2 + k3 h3
∑k h
j =1
j j
k= 1 1 =
h h
Para el promedio en serie, Ver Fig. 2.9, la caída de presión total corresponde a las
contribuciones individuales de caídas de presión en cada uno de los elementos
constituyentes. Matemáticamente:
qµ L q µL q µL q µL
= 1 1 + 2 2 + 3 3
1.127kA 1.127k1 A 1.127 k2 A 1.127 k3 A
q k1 , h 1
k2, h2
k 3 , h3
w
L
q
k1 k2 k3
L1 L2 L3
k1
q1, A1
q2, A2 k2
q3, A3 k3
L
q1, k1, h1
q2, k2, h2
q3, k3, h3
L L
k= =
L1 L2 L3 n
+ +
k1 k 2 k3 ∑L
j =1
j /kj
q = q1 + q2 + q3
re
ke
h
rs
ks
rw
k A + k 2 A2 + k3 A3
∑k j =1
j Aj
k= 1 1 =
At At
q = q1 + q2 + q3
k h + k 2 h2 + k3 h3
∑k h
j =1
j j
k= 1 1 =
h h
∆P = ∆Pw− s + ∆Ps −e
k s k e ln(re / rw )
k=
k s ln(re / rs ) + k ea ln(rs / rw )
Flujo
Agua
El fenómeno presente en la Fig. 2.13 se debe a que Pc1 (presión capilar) quiere mover el oil
hacia la izquierda y Pc2 hacia la derecha, puesto que el radio en la izquierda es menor, la
presión capilar, Pc, es mayor y se requiere cierto gradiente de presión (mayor que Pc) para
mover la burbuja.
1.0
krg
Fase no mojante
kr
kro
Fase mojante
Concava hacia
arriba
0
Swc+Sor 1.0
SL
kro
0.9
Fase no mojante
(Forma de S)
0.8
0.7
0.6
kr krw
0.5
0.4
Fase mojante
0.3
(Concava hacia
arriba)
0.2
0.1
0
0 10 20 30 40 50 60 70 80 90 100
Sw
1
Gas
0.9 Aceite
0.8
Permeabilidad relativa
0.7
0.6
0.5
0.4
0.3
0.2
0.1
Inicio Máximo
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Saturación de líquidos
1 4000
0.8 3200
Saturación de líquidos
0.7 2800
Presión, psi
0.6 2400
S critica condensados
0.5 2000
0.4 1600
0.3 1200
S wi
0.2 800
0.1 400
0 0
0.1 1 10 100 1000 10000
Radio, ft
• REGIÓN 3 (PR > PD – Presión de rocío) : Fluye gas condensado en una fase, igual a la
saturación inicial de gas (Sgi = 1- Swi).
• REGION 2 (PR < PD): Fluye gas condensado en una fase. Se acumula condensado
retrógrado porque SL < Scc (saturación crítica de condensados). Disminuye la krg.
• REGION 1 (PR < PD): Fluyen 2 fases a composiciones constantes. Hay pérdida de
productividad debido a las disminución de krg.
krg
= A(0.0435 + 0.4556 A)
kro
1 − S gc − S wi − So
A=
So − Sor
krw = (1 − S ) 2
kro = S 3
So
S=
1 − S wi
krw = S 3
kro = (1 − S )3
S w − S wi
S=
1 − S wi
krw = (1 − S ) 4
kro = S 3 (2 − S )
So
S=
1 − S wi
krw = S 4
kro = (1 − 2 S )1.5 ⎡⎣ 2 − 1 − 2S ⎤⎦
S w − S wi
S=
1 − S wi
kro = (1 − S ) 4
krg = S 3 (2 − S )
Sg
S=
1 − S wi
krw = (1 − S )3
kro = S 3
S w − S wi
S=
1 − S wi
krw = (1 − S ) 2 (1 − S 1.5 )
kro = S 3.5
krw = (1 − S ) 2 (1 − S 2 )
kro = S 4
Correlación de Torcaso y Willie1. Esta solo es válida para procesos de drenaje en areniscas
no consolidadas.
krg (1 − S * ) 2 ⎡⎣1 − ( S * ) 2 ⎤⎦
=
kro ( S * )4
So
S* =
1 − S wi
Correlación de Pirson. Para rocas con porosidad intergranular y para flujo de gas en
condiciones de drenaje, las ecuaciones son:
krw = S wn S w3
krg = (1 − S wn ) ⎡⎣ S w 1 − S wn ⎤⎦
S w − S wi
S wn =
1 − S wi
Para rocas con porosidad intergranular y para flujo simultáneo de agua y petróleo en
condiciones de imbibición, las ecuaciones son:
krw = S wn S w3
⎛ 1 − S w − S wi ⎞
kro = ⎜ ⎟
⎝ 1 − Sor − S wi ⎠
Correlación de Spivak.
kro = Son4
krw = S wn
2
⎛ S − Sor ⎞
Son = ⎜ o ⎟
⎝ 1 − Sor − S wi ⎠
⎛ S − S wi ⎞
S wn = ⎜ w ⎟
⎝ 1 − Sor − S wi ⎠
3
⎛ S − S wc ⎞
krw = ⎜ w ⎟
⎝ 1 − S wc ⎠
( S − S wc )
2.1
kro = w
(1 − Swc )
2
S g3
krg =
(1 − Swi )
3
So3
kro =
(1 − Swi )
3
( S w − S wi )3
krw =
(1 − Swi )
3
S g2 (1 − S wi ) 2 − ( S w + So − S wi ) 2
krg =
(1 − Swi )
4
So3 (2S w + So − 2S wi )
kro =
(1 − Swi )
4
( S w − S wi ) 4
kro =
(1 − Swi )
4
Las ecuaciones anteriores pueden aplicarse a sistemas mojados por petróleo siempre y
cuando la saturación de agua, Sw, se cambie por la saturación de petróleo, So, y viceversa.
Para unas saturaciones de agua y gas dadas, obtenga krw y krow de las curvas para el sistema
agua-petróleo usando la saturación de agua y obtenga (usando la saturación de gas) krg y
krog de las curvas para el sistema gas-petróleo. En caso que kro esté definida con respecto a
EJEMPLO
Un pozo produce de un yacimiento formado por una arena consolidada cuya saturación
inicial era de 20 % y en la actualidad es de 35 %. Ver Fig. 2.15. El yacimiento está rodeado
por un acuífero cuya presión es de 3000 psi y que se mantiene constante a través del
tiempo. La arena consta de dos estratos con diferente permeabilidades iniciales. La
permeabilidad de la arena superior es de 10 darcies y de la arena inferior es de 40 darcies.
Iniciada la producción se presentó producción de escamas en un radio de 200 pies desde el
pozo que causó una reducción de la permeabilidad en un 50 % en ambas capas.
Adicionalmente, debido a la depositación de parafinas la permeabilidad se reduce en ambos
estratos en un 30 % en un radio de 2000 pies medidos desde el pozo. Otros datos
adicionales:
Pwf=500 psia
5000 ft
5 ft
oil
15 ft
ft 5000 ft
5000 6000 ft
P = 3000 psi
SOLUCIÓN
a) Debido a los dos mecanismos de daño que ha sufrido la formación las permeabilidades
quedan de la siguiente manera:
Efectuando un promedio aritmético ponderado por espesor en cada una de las zonas se
tiene:
k s hs + ki hi (5000)(5) + (20000)(15)
k = = = 16250 md
hs + hi 5 + 15
Radio k , ponderado por espesor, md
0.5 - 200 ft 16250
200 - 2000 ft 22750
2000 - 5000 32500
ft
De donde k = 18490.51 md
0.00708kkro h∆Pt
qo =
µ ln(re / rw )
141.5 141.5
γ= = = 0.865
API + 131.5 32 + 131.5
S w − S wi 0.35 − 0.2
S= = = 0.1875
1 − S wi 1 − 0.2
0.83
⎛R ⎞ ⎛ 5 ⎞
0.83
Pb = 18 ⎜ sb ⎟ 10 = 18 ⎜
yg
⎟ 10−0.2362 = 47.83 psi
⎜γ ⎟ ⎝ 0.8 ⎠
⎝ g ⎠
EJEMPLO
SOLUCIÓN
Asuma un valor de rw = 0.5 ft y re = 1000 ft. Se sabe que para flujo radial en serie:
De donde:
k = 1.17 k1
k = 1.363k1
q k ⎡∂ P ⎤
u= =− ⎢ + ρ g sin θ ⎥ (2.1)
A µ ⎣ ∂l ⎦
q k ∆P
u= = (2.2)
A µ ∆L
Discusión. Las propiedades de los materiales porosos usados en las pruebas de flujo y los
datos de diferencial de presión y ratas de flujo se presentan en las tablas 2.1 y 2.2,
respectivamente. La Fig. 2.16 muestra el gráfico de los datos de q/A vs. ∆p/∆L obtenidos
experimentalmente usando una salmuera saturada que fluía a través de arenisca arcillosa,
una arenisca marrón y empaquetamientos. Observe que estos dan una línea recta que está
algo desplazada del origen. Prada y Civan sugieren que este desplazamiento indica que la
Ec. 2.2 debe corregirse, así:
q k ⎡ ∆P ⎛ ∆P ⎞ ⎤
u= = ⎢ −⎜ ⎟ ⎥
A µ ⎣ ∆L ⎝ ∆L ⎠cr ⎦
∆P ⎛ ∆P ⎞
cuando >⎜ ⎟ (2.3a)
∆L ⎝ ∆L ⎠cr
U = q/A = 0 (2.3b)
Donde (∆P/∆L)cr representa el umbral del gradiente de presión por debajo del cual el fluido
no puede fluir porque la fuerza de flujo no es suficiente para contrarrestar los efectos de
fricción.
La tabla 3 muestra que el umbral del gradiente de presión obtenido de la Fig. 2.16 decrece a
medida que la movilidad del fluido, k/µ, se incrementa. Típicamente, el valor del umbral
del gradiente de presión es despreciable para gases pero puede ser significante para
líquidos. Un gráfico de estos datos se reporta en la Fig. 2.17 en un gráfico log-log sobre un
amplio rango de movilidades, 57 < (k/µ) < 2.8x104 md/cp. El mejor ajuste por mínimos
cuadrados de estos datos:
−0.8
⎛ ∆P ⎞ ⎛ 1000k ⎞
⎜ ⎟ = 16 ⎜ ⎟ (2.4)
⎝ ∆L ⎠cr ⎝ µ ⎠
Es importante clarificar que la permeabilidad en la Ec. 2.4 está dada en mD. El coeficiente
de correlación es 0.96, indicando una fuerte correlación e los datos medidos con la Ec. 2.4.
Basado en esta ecuación, es aparente que el umbral del gradiente de presión se desprecia
solamente cuando existen altas movilidades de los fluidos. Por lo tanto, para bajas
movilidades, la ley de Darcy debe corregirse como sigue:
q k ⎡∂ P ⎛ ∂ P ⎞ ⎤
u= =− ⎢ −⎜ ⎟ + ρ g sin θ ⎥
A µ ⎢⎣ ∂ l ⎝ ∂ l ⎠cr ⎥⎦
Tabla 2.1. Propiedades de los materiales porosos usados en los experimentos de flujo8
Sand packs
Sandpack #1 Sandpack #2 Sandpack #3
∆p (psi) q (cc/min) ∆p (psi) q (cc/min) ∆p (psi) q (cc/min)
0.50 3.81 0.60 2.47 0.50 4.57
0.80 13.34 1.00 12.34 1.00 11.59
1.20 25.25 1.20 15.12 1.40 21.91
1.80 34.30 1.50 20.34 1.80 29.92
Shaly sandstone
∆p (psi) q (cc/min)
9.89 2
14.95 4
19.52 6
Tabla 2.3. Umbral del gradiente de presión vs. la movilidad del fluido8
Muestra (∆P/∆L)cr k/µ (md/cp)
(psi/cm)
Brown sandstone #1 8.46E-03 6.52E+03
Brown sandstone #2 l.40E-02 4.81E+03
Brown sandstone #3 2.03E-02 4.73E+03
Brown sandstone #4 l.70E-02 3.13E+03
Sandpack #1 4.60E-03 2.79E+04
Sandpack #2 7.89E-03 2.35E+04
Sandpack #3 5.84E-03 2.38E+04
Shaly sandstone 8.13E-0l 5.66E+0l
0.10
Brownsanstone # 4
Brownsandstone # 3
Brownsandstone # 2
0.08
Brownstone # 1
Sandpack # 1
(∆p/∆L), psi/cm
Sandpack # 3
0.06 Sandpack # 2
0.04
0.02
0.00
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
q/A, cm/min
⎛∂P ⎞ ⎛∂P⎞
cuando − ⎜ + ρ g sin θ ⎟ > − ⎜ ⎟ (2.5a)
⎝ ∂l ⎠ ⎝ ∂ l ⎠cr
q
u= =0
A
⎛∂P ⎞ ⎛∂P⎞
cuando − ⎜ + ρ g sin θ ⎟ ≤ − ⎜ ⎟ (2.5b)
⎝ ∂l ⎠ ⎝ ∂ l ⎠cr
1.E+00
−0.8
⎛ ∆p ⎞ ⎛k⎞
⎜ ⎟ = 16 ⎜ ⎟
⎝ ∆L ⎠cr ⎝µ⎠
2
R = 0.96
, psi/cm
1.E-01
cr
(∆p/∆L)
1.E-02
1.E-03
1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05
K/µ, md/cp
Fig. 2.17. Gráfica del umbral del gradiente de presión vs. Movilidad del fluido8
−n
⎛∂P⎞ ⎛ 1000k ⎞
−⎜ ⎟ = m⎜ ⎟ (2.6)
⎝ ∂ l ⎠cr ⎝ 1 ⎠
k ⎪⎧ ⎡ ∂ P ⎡1000k ⎤ ⎪⎫
−n
q ⎤
u = = − ⎨⎢ + ρ g sin θ ⎥ − m ⎢ ⎥ ⎬
A µ ⎩⎪ ⎣ ∂ l ⎦ ⎣ µ ⎦ ⎭⎪
−n
⎛∂P ⎞ ⎛ 1000k ⎞
cuando − ⎜ + ρ g sin θ ⎟ > m ⎜ ⎟ (2.7a)
⎝ ∂l ⎠ ⎝ µ ⎠
−n
q ⎛∂P ⎞ ⎛ 1000k ⎞
u = = 0, −⎜ + ρ g sin θ ⎟ ≤ m ⎜ ⎟ (2.7b)
A ⎝ ∂l ⎠ ⎝ µ ⎠
Como puede observarse en las Ecs. 2.7a y 2.7b, la corrección es aditiva. El siguiente
ejercicio indica que los fluidos fluyen a través de medios porosos solo si la fuerza sobre el
fluido es suficiente para contrarrestar el umbral del gradiente de presión y, por lo tanto, la
ley de Darcy debe corregirse para tener en cuenta este efecto. La corrección en la forma de
la Ec. 2.6 ha sido obtenida mediante regresión de un número limitado de medidas y usando
solo tres diferentes medios porosos y salmuera. Por lo tanto, para otros tipos de fluido y
rocas, se requiere mejorar las Correlaciónes del umbral del gradiente de presión. En
términos generales, el valor del umbral del gradiente de presión es una medida del grado de
la mojabilidad de la roca, de la tensión interfacial petróleo-agua, y del tamaño de los poros
más grandes en el exterior de la muestra de roca. Una alta presión de desplazamiento indica
ya sea un fuerte grado de mojabilidad o poros pequeños, o ambos efectos.
Entonces 0.1 > 0.0637. Por ende, hay suficiente presión para que exista flujo. La velocidad
de flujo es calculada usando la Ec. 2.7a;
k ⎪⎧ ⎡ ∂ P ⎡1000k ⎤ ⎫⎪
−n
1 ⎪⎧⎛ 3 − 1 ⎞ ⎛ 1000(1) ⎞ ⎪⎫
−0.8
⎤
u = − ⎨⎢ + ρ g sin θ ⎥ − m ⎢ ⎥ ⎬ = − 1 ⎨⎜ 30 ⎟ − 16 ⎜ 1 ⎟ ⎬ = −0.0694 cm / s
µ ⎪⎩ ⎣ ∂ l ⎦ ⎣ µ ⎦ ⎪⎭ ⎪⎩⎝ ⎠ ⎝ ⎠ ⎪⎭
q k ∆P 1 ⎛ 3 − 1 ⎞
u= =− = = −0.1 cm / s
A µ ∆L 1 ⎜⎝ 30 ⎟⎠
Como puede verse, hay una marcada diferencia entre las dos velocidades. La mínima
presión para tener flujo es:
−0.8
⎛∂P⎞ ⎛ 1000k ⎞
−⎜ ⎟ = 16 ⎜ ⎟
⎝ ∂l ⎠ ⎝ µ ⎠
⎛∂P⎞
⎜ ⎟ = 0.0637 atm / cm
⎝ ∂l ⎠
Para la mayoría de los fluidos hidrocarburos, el esfuerzo de corte y la rata de corte pueden
describirse mediante la ley de fricción de Newton la cual combinada con la ecuación de
movimiento resulta en la bien conocida ecuación de Navier-Stokes. La solución de dicha
ecuación para las condiciones de frontera apropiadas da lugar a la distribución de velocidad
del problema dado. Sin embargo, la geometría de los poros, no permite la formulación
adecuada de las condiciones de frontera a través del medio poroso. Luego, una
aproximación diferente se debe tomar. Darcy descubrió una relación simple entre el
gradiente de presión y el vector velocidad para una sola fase4.
Pero,
1 dV
c=−
V dP
Luego;
dV = −cVdP
dV = −c(2π rhdr )φ dP
∂V
Si dq = entonces:
∂t
∂P
dq = −cφ (2π rhdr )
∂t
ó;
∂q ∂P
= −cφ (2π rh) (2.9)
∂r ∂t
k ∂P
q = − (2π rh) (2.10)
µ ∂r
r+dr
P+dP P
r
Well
h
∂q k ⎡ ∂P ∂2P ⎤
= − (2π h) ⎢ + r 2 ⎥ (2.11)
∂r µ ⎣ ∂r ∂r ⎦
∂P k ⎡ ∂P ∂2P ⎤
− cφ (2π rh) = − (2π h) ⎢ + r 2 ⎥
∂t µ ⎣ ∂r ∂r ⎦
ó;
∂P k ⎡ ∂P ∂2P ⎤
cφ r = ⎢ +r 2 ⎥
∂t µ ⎣ ∂r ∂r ⎦
Rearreglando,
∂ 2 P 1 ∂P cφ µ ∂P
+ = (2.12)
∂r 2 r ∂r k ∂t
∂ 2 P 1 ∂P 1 ∂ ⎛ ∂P ⎞
+ = ⎜r ⎟=0
∂r 2 r ∂r r ∂r ⎝ ∂r ⎠
Integrando:
∂P
r = c1 (2.13)
∂r
Imponiendo una condición de frontera en la cara del pozo (de la ley de Darcy):
k ∂P
q = (2π rh)
µ ∂r
∂P qµ 1
= (2.14)
2πk h rw
rw
∂r
⎛ qµ 1 ⎞
rw ⎜ ⎟ = c1
⎝ 2π k h rw ⎠
qµ
c1 =
2π k h
Reemplazando en (2.14):
∂P qµ
r =
∂r 2πk h
Separando variables:
qµ
re Pe
dr
∫ = ∫ dp
2πk h rw r Pw
Integrando:
7.08kh ( Pe − Pw)
q=
µ ln(re / rw )
re
2
P= 2
re ∫ Prdr
rw
(2.15)
si,
7.08kh ( P − Pw)
q=
µ ln ( r / rw )
2
re
⎛ qµ r ⎞
P=
re2 ∫ ⎜⎜⎝ P r + 2πkh r ln r
rw
w
⎟⎟dr
w ⎠
x2 x2
∫ x ln xdx = − 4
+ ln x
2
Simplificando:
⎛ qµ ⎡ 1 re ⎤ ⎞
P = ⎜⎜ Pw + ⎢− + ln ⎥ ⎟⎟
⎝ 2πkh ⎣ 2 rw ⎦ ⎠
de donde:
7.08kh( P − Pw )
q=
µ B [ ln(re / rw ) − 0.5]
dP − q
=
dt cV p
1 ∂ ⎛ ∂P ⎞ cφ µ q
⎜r ⎟=−
r ∂r ⎝ ∂r ⎠ k cV p
1 ∂ ⎛ ∂P ⎞ qµ
⎜r ⎟=−
r ∂r ⎝ ∂r ⎠ πkre2 h
Separando:
⎛ ∂P ⎞ qµ
∫ ⎜⎝ r ∂r ⎟⎠ = − πkr h ∫ rdr
e
2
Después de integrar:
⎛ ∂P ⎞ qµ r 2
⎜r ⎟=− + c1
⎝ ∂r ⎠ 2πkre2 h
Aplicando la condición de frontera cuando r = re, dP/dr =0 porque el sistema debe ser
cerrado. Entonces6:
q µ re2
c1 =
2π kre2 h
Luego:
⎛ ∂P ⎞ qµ r qµ
⎜ ⎟=− +
⎝ ∂r ⎠ 2πkre h 2πkhr
2
Separando;
qµ qµ
Pe re re
dr
∫Pw dP = − 2π kre2 h rw∫ rdr + 2π kh rw∫ r
Integrando:
q µ ⎛ re2 − rw2 ⎞ q µ r
Pe − Pw = − ⎜ ⎟+ ln e
2π kh ⎝ 2re ⎠ 2π kh rw
2
Puesto que re2 es mucho mayor que rw2, el termino rw2 se puede ignorar.
q µ ⎛ re ⎞
Pe − Pw = ⎜ ln − 0.5 ⎟
2π kh ⎝ rw ⎠
Si e-0.5 = 0.606
q µ ⎛ 0.606re ⎞
Pe − Pw = ⎜ ln ⎟
2π kh ⎝ rw ⎠
7.08kh( P − Pw ) 7.08kh( P − Pw )
q= =
⎛ r ⎞ ⎛ 0.472re ⎞
µB⎜⎜ ln e − 0.75 ⎟⎟ µB⎜⎜ ln ⎟
⎝ rw ⎠ ⎝ rw ⎟⎠
q µ ⎛ r 2 − rw2 ⎞ qµ r
Pe − Pw = − ⎜ ⎟+ ln
2π kh ⎝ 2re2 ⎠ 2π kh rw
qµ ⎛ r 2 ⎞ qµ r
Pe − Pw = − ⎜ 2 ⎟+ ln
2π kh ⎝ 2re ⎠ 2π kh rw
Factorizando:
qµ ⎛ r r2 ⎞
P − Pw = ⎜ ln − ⎟ (2.16)
2π kh ⎝ rw 2re2 ⎠
re
P=
∫ PdV rw
re
∫ dV rw
Puesto que:
dV re
rw = 2π rhφ dr re
rw entonces:
re
2π hφ ∫ Prdr
rw
P= re
2π hφ ∫ rdr
rw
re
2 ∫ Prdr
rw
P=
(re2 − rw2 )
2 re
P=
re2 ∫ rw
Prdr
2 ⎧⎪ qµ ⎛ r
re r2 ⎞ ⎫⎪
P=
re2 ∫rw ⎪ 2π kh ⎝ rw 2re ⎠ ⎬⎪ rdr
⎨ ⎜ ln − 2 ⎟ + Pw
⎩ ⎭
qµ ⎛ r
re r2 ⎞ 2 re
P=
π khre2 ∫rw ⎜⎝ ln rw − 2re2 ⎟⎠ rdr + re2 ∫ rw
Pw rdr
qµ ⎛
re r r3 ⎞ 2 re2 − rw2
P=
π khre2 ∫rw ⎜⎝ r ln rw − 2re2 ⎟⎠ dr + re2 Pw 2
qµ ⎡ re r ⎛ re2 ⎞ ⎤
π khre2 ⎣ ∫rw
P= ⎢ r ln dr − ⎜ ⎟ ⎥ + Pw (2.17)
rw ⎝ 8 ⎠⎦
re r
Pero, ∫ rw
r ln
rw
dr debe integrarse por partes:
∫ udv = uv − ∫ vdu
r 1
u = ln entonces: du =
rw (r − rw )rw
r2
dv = rdr entonces: v =
2
re r r2 r re r2 1
∫ dr = ln −∫
re
r ln rw dr
rw rw 2 rw rw 2 r
re r r2 r r2 r ⎡ r2 r2 ⎤
∫rw
r ln
rw
dr = e ln e − w ln w − ⎢ e − w ⎥
2 rw 2 rw ⎣ 4 4 ⎦
re r re2 re re2
∫rw r ln rw dr = 2 ln rw − 4
Reemplazando la anterior respuesta en la Ec. 2.17, se tiene:
q µ ⎡ 1 re 1 1 ⎤
P − Pw = ln − −
π kh ⎢⎣ 2 rw 4 8 ⎥⎦
q µ ⎡ re 1 1 ⎤
P − Pw = ⎢ln − − ⎥
2π kh ⎣ rw 2 4 ⎦
q µ ⎡ re 3 ⎤
P − Pw = ⎢ln − ⎥
2π kh ⎣ rw 4 ⎦
2.12. MOVILIDAD
kf
λf =
µf
En un proceso de flujo multifásico, existe relación entre las movilidades de los fluidos, a
esto se le conoce como relación de movilidad, M, normalmente se expresa como la relación
entre el fluido desplazante sobre el desplazado (algunos autores consideran la definición
contraria). Si el fluido desplazante es agua2:
λw
M=
λo
Si M < 1, significa que el crudo se mueve más fácilmente que el agua, si M = 1 significa
que ambos fluidos tienen igual movilidad y si M > 1, significa que el agua es muy móvil
con respecto al crudo.
La interfase que separa a dos fases es una región con solubilidad limitada, que a lo sumo
tiene el espesor de unas pocas moléculas. Se puede visualizar como una barrera que se
forma debido a que las fuerzas atractivas entre las moléculas de la misma fase son mucho
mayores que aquellas que existen en dos fases diferentes. La tensión superficial es una
propiedad termodinámica fundamental de la interfase. Se define como la energía
disponible para incrementar el área de la interfase en una unidad. Cuando dos fluidos
están en contacto, las moléculas cerca a la interfase se atraen desigualmente por sus vecinas
porque unas son más grandes que las otras, esto origina una superficie de energía
libre/unidad de área que se llama tensión interfacial. En otras palabras, es la unidad de
P = 40 + 2.38 * PM liq
2.14. MOJABILIDAD
Sistema repulsivo
θ > 90 grados
Aire
Mercurio
El mercurio trata
de extenderse al
VIDRIO mínimo
Oil o aire
Agua
De acuerdo con la Fig. 2.22, para que exista el proceso de IMBIBICIÓN se requiere que θ
< 90°, P2 > P1, Pc = P2 - P1 > 0. Otras particularidades de la imbibición son a) Pe > Ps
Imbibición forzada, b) Pe = Ps Imbibición espontánea, y c) Pe < Ps Imbibición inhibida. Si θ
< 90° existe succión. Aquí se requiere contrapresión (igual a Pc) para que no avance el
menisco.
En un sistema con DRENAJE, ver Fig. 2.23, θ > 90°, P1 > P2, Pc = P2 - P1 < 0 Pe - Ps > Pc
existe drenaje forzado.
Water-wet Neutrally-wet
Oil Oil
θ
Agua Agua
Oil-wet
Oil
θ
Agua
V
θ
Pe P1 P2 Ps
oil
agua
+ mojante - mojante
P P2
Pe
P1
P2 Ps
Pe Ps
P1
P2
Pe Ps
P1
x
Fig. 2.22. Principio de la imbibición
θ
Pe P1 P2 Ps
oil
agua
- mojante + mojante
P
Pe P1
P2
Ps
. x
σ ow
OIL
WATER
ANGULO DE CONTACTO
σ os σ ws
SOLID
∑ F = −σ OS + σ OW Cosθ + σ WS
⎡ σ OS − σ WS ⎤
θ = arccos⎢ ⎥
⎣ σ OW ⎦
Si las tensiones de adhesión son iguales, θ = 90°. De aquí también se tiene que:
θ es medido con microscopio. θ cambia con el tiempo, pero se estabiliza después de unos
pocos días. θ es función de la fuerza de adhesión líquido-sólido. En medios porosos, se
recurre a una medida indirecta de la mojabilidad, la cual se conoce como índice de
mojabilidad, IW, o índice USBM que va desde 1 para altamente humectable por agua hasta
-1.5 para fuertemente humectable por petróleo. Ver Figs. 2.25 y 2.26. Luego esto incide en
las curvas de permeabilidad relativa. Fig. 2.27.
⎛A ⎞
IW = log ⎜ 1 ⎟
⎝ A2 ⎠
10
1
2
A1
Pc
A2
4
3
-10
Sw
oil oil
oil oil
Permeabilidad relativa
0.6
Aceite
0.4
0.2
Agua
0
0 20 40 60 80 100
Sw
Oil Oil
Agua Agua
θ θ
Una propiedad fundamental de los líquidos es la tendencia a contraerse y dar la menor área
superficial posible produciendo una forma esférica en gotas pequeñas. La explicación a este
comportamiento puede expresarse como un desbalance de fuerzas moleculares atractivas en
la superficie del líquido. Considere un líquido coexistiendo con gas. Las moléculas en el
interior están rodeadas por otras en todos los lados, sujetándolas a una atracción molecular
uniforme en todas las direcciones. En la superficie externa, sin embargo, las moléculas son
atraídas hacia adentro en todos lados, pero no existe fuerza atracción hacia fuera para
balancear la fuerza hacia adentro. Este desbalance de fuerzas causa que la superficie se
contraiga a la menor área posible y produce la tensión superficial, σ. Se debe efectuar
trabajo para extender la superficie en oposición a la tensión superficial forzando las
moléculas desde el interior a la superficie. Esto indica que existe una energía libre asociada
con la superficie que tiene las mismas dimensiones que la tensión superficial.
Considere un segmento de une superficie interfacial que separa a dos fluidos con presiones
diferenciales a través de la interfase, produciendo un rectángulo curvilíneo como se ilustra
en la Fig. 2.29.a. Ambos centros de curvatura están al mismo lado, por lo tanto, R1 y R2 son
positivos. El trabajo realizado para expandir la superficie mediante el incremento de la
presión en el lado convexo, es el trabajo contra la tensión superficial. Las longitudes de los
arcos en los lados se incrementa de L1 a L1 + (L1/R1)dz y de L2 a L2 + (L2/R2)dz. El área de
la superficie original ABCD se expande a la superficie A´B´C´D´ donde6:
D'
A'
dz
D
A B'
C'
B
α
β C
R1 R2
O'
O'
Area ABCD = L1 × L2
⎛ L ⎞ ⎛ L ⎞
Area A ' B ' C ' D ' = ⎜ L1 + 1 dz ⎟ × ⎜ L2 + 2 dz ⎟
⎝ R1 ⎠ ⎝ R2 ⎠
⎛ dz dz dz 2 ⎞
Area A ' B ' C ' D ' = L1 L2 ⎜1 + + + ⎟
⎝ R1 R2 R1 R2 ⎠
⎛ 1 1 ⎞
A ' B ' C ' D '− ABCD = L1 L2 dz ⎜ + ⎟
⎝ R1 R2 ⎠
El trabajo isotérmico contra la tensión superficial requerido para expandir el área es:
⎛1 1 ⎞
Trabajo(1) = σ L1 L2 dz ⎜ + ⎟
⎝ R1 R2 ⎠
Trabajo(2) = PL1 L2 dz
Aceite Agua
x
σnws σws
Línea de contacto
Pnw Pw
σnws σws
Fig. 2.29.c. Diagrama de fuerzas
Igualando los dos trabajos y cancelando términos comunes da la presión capilar en función
de la tensión interfacial y los radios de curvatura:
⎛ 1 1 ⎞
Pc = σ ⎜ + ⎟
⎝ R1 R2 ⎠
Cuando se habla de medios porosos, R2 es negativo (ver Fig. 2.31), por lo tanto, una forma
más general de la ecuación anterior es:
⎛ 1 1 ⎞
Pc = σ ⎜ ± ⎟
⎝ R1 R2 ⎠
Considere el sistema de la Fig. 2.29.c donde existe equilibrio estático. El agua moja
fuertemente la superficie con un ángulo de contacto. Se puede observar que la presión de
petróleo es mayor que la de agua, sin importar la longitud del tubo. El agua puede
desplazarse mediante inyección de aceite. El aceite se desplazará espontáneamente si la
presión de la fase de aceite se reduce, aunque la presión en la fase de agua es menor que la
de aceite. Este fenómeno puede verse en el análisis de fuerzas dado a continuación. En
equilibrio, ΣFx = 0, luego:
2(σ nws − σ ws )
Pnw − Pw =
r
σ os − σ ws = σ owCosθ
σ nws − σ ws = σ nwwCosθ
Luego:
2σ nwwCosθ
Pnw − Pw =
r
Por convención la presión capilar Po - Pw es negativa para sistemas mojados por aceite. En
términos generales, la presión capilar se define como la diferencia de presión entre la
presión de la fase mojante y no mojante y siempre se considera positiva. Existen
formaciones menos mojables, intermediamente mojables y fuertemente mojable. La presión
capilar tiene aplicaciones en simulación de yacimientos y en ingeniería de yacimientos para
calcular, principalmente, la altura de la zona de transición y la saturación de agua
irreducible.
Pc = Pnw − Pw
Po
h
Pw θ
Oil
Water
Gravedad
r2
r1
2π rσ cos θ
Pc =
π r2
2 σ cos θ
Pc =
r
⎛1 1⎞
Pc = Pnw − Pw = Po − Pw = σ ⎜ + ⎟
⎝ r1 r2 ⎠
Po
Pc
Pw
h
Zona de
transicion
Sw
donde r1 se mide sobre el petróleo. En términos generales r1 y r2 son los radios de curvatura
en el mismo lado de la interfase, para el caso del tubo capilar ambos tienen valores
positivos. De lo contrario el radio más pequeño es positivo y el más grande es negativo. r1
y r2 son iguales a r/cosθ para un capilar uniforme de radio r. En la Fig. 2.30, a la altura h,
cuando existe equilibrio hidrostático capilaridad-gravedad, las presiones deben ser iguales,
para el caso del aceite y agua, respectivamente:
P = Po + ρ o gh
P = Pw + ρ w gh
Mediante substracción:
Pc = Po − Pw = ∆ρ gh
⎛1 1⎞
Pc = σ ⎜ − ⎟ = ∆ ρ gh
⎝R R⎠
⎛2⎞
Pc = σ ⎜ ⎟ = ∆ ρ gh
⎝R⎠
2 σ cos θ
Pc = Po − Pw = = ∆ ρ gh
r
Existe un grupo adimensional que representa la relación de las fuerzas viscosas y las
fuerzas capilares. Este grupo se conoció como número capilar, Nca:
Fuerzas Viscosas v µw
N ca = =
Fuerzas Capilares σ ow cos θ
u µw
N ca =
σ ow
EJEMPLO
SOLUCIÓN
La Tensión superficial es la deferencia de presiones entre dos fases, en este caso el líquido
y su vapor.
Pc r
σw = (a)
2Cosθ
Pc = ρwgh (c)
ρ w ghr
σw = (d)
2Cosθ
EJEMPLO
Calcule el gradiente de presión para el flujo de agua con viscosidad de 1 cp que fluye a una
rata de un yacimiento de 1 ft/día a través de un capilar con diámetro de 0.04 pulgadas.
SOLUCIÓN
∆P µV
= − 6.22*10−8 2
L r gc
r t = 0.04 in
r gota = 0.104 in
∆P (1cp ) (1 ft / día )
= − 6.22*10−8 = −0.0007 psi / ft
lbm ft
L
(1.66*10 )
−3 2
(32.2
lbf s 2
)
EJEMPLO
Para determinar la tensión interfacial entre agua y crudo, el agua está siendo desplazada a
través de un tubo en un baño de crudo. El agua forma una gota esférica la cual se rompe o
separa del agua en el tubo de inyección cuando alcanza cierto tamaño y cae al fondo del
recipiente. Ver Fig. 2.32.a. Otros datos son: ρw = 1.02 gr/cc, ρo = 0.8 gr/cc12.
SOLUCIÓN
( ρ w − ρo )Vgota g − IFT ⋅ L = 0
El hecho que las curvas de presión capilar de casi todos los materiales naturalmente
porosos tienen muchas características en común y ha conducido a introducir algunas
ecuaciones generales desde el punto de vista adimensional. Basado en el hecho que la
presión depende de la porosidad, la fuerza interfacial y la geometría del poro, Leverett
definió su función adimensional de saturación la cual la llamó la función J. La Función J
de Leverett es una herramienta poderosa para el desarrollo de técnicas nuevas y veloces de
mejoramiento de procesos que intentan obtener un punto de vista más exacto de la
distribución de los fluidos del yacimiento. Algunos de estas aplicaciones son: recobro de
inyección de agua, inicialización de modelos de simulación, que mejoran la distribución
inicial de la saturación de agua para una representación tridimensional, distribución de
presiones y saturaciones en yacimientos dinámicos y otros. Esta función tiene una gran
importancia y una relación fuerte con el mejor entendimiento del comportamiento del
yacimiento9,10. Referirse a las Figs. 2.33.b a la 2.36.
Pc k
J (Sw ) =
σ φ
Esta función adimensional sirve muy bien en muchos casos para remover discrepancias en
las curvas de presión capilar y las reduce a una curva común. Para derivar J, se tiene en
cuenta que la permeabilidad tiene unidades L2, luego el radio en la Ec. de Pc se substituye
por (k/φ)0.5. Las curvas de permeabilidad relativa también sufren de histéresis que está
relacionada con la histéresis de las curvas de presión capilar. Yacimientos que son
candidatos a ser inyectados con agua usualmente tienen saturaciones de petróleo inferior a
la inicial. Luego, la condición inicial del yacimiento estará sobre uno de los caminos de
desplazamiento. Si el proceso de recobro moviliza aceite, creando un banco de petróleo con
saturación incrementada de petróleo, el camino de saturación cambia de imbibición a
drenaje o viceversa. La histéresis del ángulo de contacto y la rugosidad de las superficies
son posibles contribuyentes a la histéresis de las permeabilidades relativas. La pendiente de
la curva de presión capilar durante drenaje, es una buena medida cualitativa del rango de
distribución del tamaño del poro. Entre más cerca de la horizontalidad o entre más plana
sea la curva de presión capilar es más uniforme el tamaño de los poros dentro de la roca.
Dadas las curvas de presión capilar medidas sobre 75 corazones tomadas al azar frente a
intervalos de yacimientos de arena. Cinco de estas curvas son mostradas en la Fig. 2.41.
Otros datos requeridos son los siguientes:
Transicion
capilar
Pc
Drenaje
Imbibicion Histeresis
Sw
Pc
Agua
imbibe
Umbral de
presión
Sw
Fig. 2.34. Curva de presión capilar en un sistema mojado por agua
El aceite drena
Sw
Fig. 2.35. Curva de presión capilar en un sistema mojado por aceite
30
20
Presion capilar, cm de Hg
Drenaje
10
Imbibicion
espontanea
0
-10
Imbibicion
forzada
-20
-30
0 20 40 60 80 100 120
Sw
Fig. 2.36. Curva de presión capilar para un sistema con mojabilidad intermedia – Se
presenta en el 47 % de los casos
2
100
2 100
No mojante
Mojante
80 80 No mojante
Permeabilidad relativa, %
Drenaje
60 60
primario
40 40
Imbibicion
20 20 secundaria
Mojante
1 3
0 0
0 20 40 60 80 100 0 20 40 60 80 100
0.8
kro
0.6
Permeabilidad relativa
0.4
Inyección
de agua
0.2
Inyección
de Petróleo
krw
0
0.3 0.4 0.5 0.6 0.7
Sw
Fig. 2.38. Histéresis de las curvas de permeabilidad relativa durante inyección de agua y de
petróleo, respectivamente
100
kro Saturación residual de
la salmuera = 13.51 %
80
Permeabilidad relativa
60 Dri-Film
40 krw
20
0
0 0.2 0.4 0.6 0.8 1
Sw
Fig. 2.39. Curvas de permeabilidades relativas para una arenisca Berea antes y después de
un tratamiento Dri-film
100
kro
krw
Incrementa mojabilidad
Permeabilidad relativa
Angulo de contacto
0
47
90
1 138
180
0.1
0 20 40 60 80 100
Sw
Fig. 2.40. Permeabilidades relativas para un rango determinado de mojabilidades (dado por
el ángulo de contacto) de la arenisca Torpedo
kav, Permeabilidad absoluta promedio del yacimiento obtenida de 1000 muestras analizadas
155 md.
Determinar:
1. Correlaciónes entre permeabilidad y saturación de agua.
2. Curva de presión capilar promedio para el yacimiento.
3. Factor de conversión para relacionar la escala de presión capilar con la elevación desde
el nivel de agua libre.
4. Saturación mínima de agua intersticial.
90
10
25
80
50
300
70 950
Presión Capilar, psi
60
50
40
30
20
10
0
0 10 20 30 40 50 60 70 80 90 100
Fig. 2.41. Curva típica de presión capilar para un yacimiento de arena (aire-agua)10
1000
Permeabilidad, md
10
0 20 40 60
Se traza una línea promedio a través de los datos para dar una relación entre permeabilidad
y saturación de agua a una presión capilar constante (ésto se hizo para presión capilar de
75 psi en la Fig. 2.42). La Fig. 2.43 muestra curvas de cinco presiones capilares; ésto fue
derivado individualmente, en la misma forma que la curva de 75 psi.
2. Las Figs. 2.41, 2.42 y 2.43 muestran un método aceptable para promediar los datos de
presión capilar. Según todos los datos de análisis de corazones se encontró una
permeabilidad promedio del yacimiento y ésta se usa para obtener la curva promedio de
presión capilar del yacimiento de datos promediados de la Fig. 2.43. Los datos de la tabla
2.5 fueron leídos de la Fig. 2.44 para una permeabilidad promedio de yacimiento de 155
md y se presentan en la Fig. 2.44. Las condiciones promedio de la Fig. 2.44 se muestran
como la curva típica de permeabilidad relativa para agua y aceite de la Fig. 2.45. La última
gráfica se construyó para mostrar la relación entre las curvas de presión capilar y
permeabilidad relativa.
1000
75 psi
50 psi
900 25 psi
10 psi
800 5 psi
Permeabilidad, md 700
600
500
400
300
200
100
0
0 10 20 30 40 50 60 70 80 90 100
Tabla 2.5. Datos de presión capilar promedio de yacimiento, (kav = 155 md)10
3. Se encontró una correlación entre presión capilar y distancia desde el nivel de agua libre
de los datos dada y la ecuación es:
σw
Pcl = ( ρ − ρ )h
σ ow w o
ó h = 3.3 Pc1
En este caso 3.3 es el factor para relacionar los datos de presión capilar promedia medida
en el laboratorio, Pcl, con la correspondiente distancia vertical desde el nivel de agua libre
(extremo inferior al final de la zona de transición o elevación más alta a la cual se
encuentra el 100% de la saturación de agua).
70 210
intersticial de agua
40 120
Zona de transición
30 90
20 60
10 Nivel de
30
agua libre
0
0 20 40 60 80 100
Fig. 2.44. Datos leídos de la Fig. 2.43 para una permeabilidad promedio de 155 md y
usados para graficar curva de presión capilar promedio de yacimiento10
0.9
Permeabilidad relativa al petróleo, k ro
0.7
0.6
0.5
0.4
0.3
0.2 Mínima
Saturación de
agua intersticial Presión
0.1 Capilar= 1 psi
0
0 10 20 30 40 50 60 70 80 90 100
DISCUSIÓN: La presión capilar esta definida como la presión diferencial a través de una
interfase entre dos fases de fluido. Esto es dado en la ecuación:
Pc = σ (1 / R1 + 1 / R2 )
Pc = ∆ρ gh
Los datos de presión capilar han sido ampliamente usados para encontrar saturación de
agua intersticial. Es necesario analizar bastantes corazones para asegurar que los datos
sean representativos. En este caso (yacimiento grande) se analizaron 75 muestras. La
Fig. 2.41 muestra la curva de presión capilar para cinco de estas muestras. Note que la
presión capilar incrementa cuando la saturación de agua disminuye. Se alcanza un punto
cuando cambios significativos en presión no cambian apreciablemente la saturación. Esta
región de la curva da la llamada mínima saturación de agua intersticial. Otra característica
de algunas de estas curvas es que una cantidad considerable de desaturación puede tomar
lugar con un pequeño incremento en la presión capilar. Esto es verdad en la curva de la
muestra de 950 md. Dicho comportamiento indica una ligera uniformidad en la muestra.
La presión capilar a la cual la fase no mojante (aceite en un sistema humectado por agua),
puede ser forzada primero dentro de los poros (por ejemplo la presión capilar mas baja a la
cual existe cierta saturación de aceite) es llamada presión de desplazamiento. En los 950
md del ejemplo de la Fig. 2.41 esta presión es 3 psi. Los conceptos teóricos indican que la
permeabilidad y la presión capilar están relacionadas. Esto es, para una roca particular y
una saturación de agua dada, la presión capilar debería incrementarse más ampliamente en
muestras ordenadas con la permeabilidad decreciente. Este fenómeno se observa en la Fig.
2.41. La distribución de la Fig. 2.42 muestra la influencia de otros factores. La
distribución de puntos de la Fig. 2.43 es similar a la Fig. 2.42.
La Fig. 2.44 muestra la curva promedio de presión capilar del yacimiento obtenida de la
Fig. 2.43 usando la permeabilidad promedio absoluta del yacimiento. Observe que esta
gráfica tiene dos escalas de presión capilar, una a condiciones de laboratorio y otra a
condiciones de yacimiento. Con esta gráfica es posible obtener un valor de saturación de
agua (para la permeabilidad promedio de la roca) en algún punto por encima del nivel de
agua libre. Por ejemplo a 90 pies la saturación de agua libre es el 40 %. Si desea, las
curvas de la Fig. 2.43 pueden ser usadas para determinar las curvas de presión capilar
medias para segmentos del yacimiento (compartimiento). A partir de datos de corazones
para cada segmento respectivo, las permeabilidades promedio son calculadas y usadas en la
Fig. 2.43.
Muchos ingenieros y geólogos definen el contacto agua-aceite como la elevación más baja
a la cual puede obtenerse la producción de aceite libre de agua. La máxima saturación de
agua por debajo de la cual solamente aceite libre será producido es definida por las
características de permeabilidad relativa petróleo-agua de la roca del yacimiento con sus
fluidos. La Fig. 2.45 muestra que aun a una saturación del 40 % muy poca agua fluye con
el aceite producido. Cuando los datos de permeabilidad relativa no son disponibles,
usualmente se asume la saturación de agua en la zona que produce aceite libre, se
representa por aquella correspondiente presión capilar de 45 a 75 psi. Esta es una
consideración razonable para arenas limpias de buena y alta porosidad y permeabilidad
porque la zona de transición posee solo unos pocos pies de altura y el mínimo valor de agua
intersticial puede ser aproximado al orden de 10 pies por encima del nivel de agua libre.
Además, hay una pequeña diferencia en la saturación de agua a este nivel y a niveles altos
con sus correspondientes presiones capilares más altas.
EJEMPLO
Cuatro juegos de presiones capilares petróleo-agua se dan en la tabla 2.6, para cuatro
corazones del mismo yacimiento. La tensión interfacial del sistema aceite-agua es 72
dinas/cm y θ = 45º. Calcule y grafique la función J. Los factores de conversión son
6894.757 Pa = 6894.57 dinas/cm2 = 1 psi, lo cual da 1 dina/cm2 = 1.450377x10-5 y 1 darcy
= 0.989633 µm2.
Pc, psi
Sw % k = 1 md k = 15 md k = 100 md k = 500 md
φ=8% φ = 11 % φ = 15 % φ = 22 %
15 93.59 28.34 12.82 6.94
16 66.48 20.13 9.10 4.93
17 46.54 14.09 6.37 3.45
18 39.89 12.08 5.46 2.96
19 33.24 10.06 4.55 2.47
20 31.82 9.63 4.36 2.36
25 22.46 6.80 3.08 1.67
30 16.85 5.10 2.31 1.25
35 13.1 3.97 1.79 0.97
40 11.23 3.40 1.54 0.83
50 7.49 2.27 1.03 0.56
60 6.55 1.98 0.90 0.49
70 6.08 1.84 0.83 0.45
80 5.62 1.70 0.77 0.42
90 5.62 1.70 0.77 0.42
100 5.62 1.70 0.77 0.42
SOLUCIÓN
Se debe tener especial cuidado en las unidades cuando se resuelve este problema. Haciendo
los cálculos para Sw = 20 %, k = 500 y φ = 0.22 se tiene9:
Pc k
J (Sw ) =
σ cos θ φ
100
90
80
70
60
J(Sw )
50
40
30
20
10
0
0 20 40 60 80 100
Sw , %
J(Sw)
Sw % k = 1 md k = 15 md k = 100 md k = 500 md
φ=8% φ = 11 % φ = 15 % φ = 22 %
15 1.4097 1.4097 1.4097 1.4097
16 1.0014 1.0014 1.0014 1.0014
17 0.7010 0.7010 0.7010 0.7010
18 0.6008 0.6008 0.6008 0.6008
19 0.5007 0.5007 0.5007 0.5007
20 0.4793 0.4793 0.4793 0.4793
25 0.3383 0.3383 0.3383 0.3383
30 0.2537 0.2537 0.2537 0.2537
35 0.1974 0.1974 0.1974 0.1974
40 0.1692 0.1692 0.1692 0.1692
50 0.1128 0.1128 0.1128 0.1128
60 0.0987 0.0987 0.0987 0.0987
70 0.0916 0.0916 0.0916 0.0916
80 0.0846 0.0846 0.0846 0.0846
90 0.0846 0.0846 0.0846 0.0846
100 0.0846 0.0846 0.0846 0.0846
En la Fig. 2.47, todos los tubos tienen el mismo radio, r. La tortuosidad, τ expresa la
relación entre el largo del corazón y el largo de un camino6.
2
⎛L ⎞
τ =⎜ h ⎟
⎝ L⎠
⎛ nπ r 4 ⎞ ∆P
q=⎜ ⎟ (2.18)
⎝ 8µ ⎠ L
⎛ kA ⎞ ∆P
q=⎜ ⎟ (2.19)
⎝ µ ⎠ L
nπ r 4
k= (2.20)
8A
Si;
Vpnπ r 2 L
φ= = (2.21)
Vt AL
Lh
nπ r 2
A= (2.22)
φ
⎛ r2 ⎞
k = ⎜ ⎟φ (2.23)
⎝ 8⎠
k
rh = 2
φ
Defina Spv = Area superficial por unidad de volumen poroso, 1/cm. En otras palabras, es la
superficie total de los poros por unidad de volumen de poros de la roca.
As n(2π rL) 2 4
S pv = = = = (2.24)
Vp nπ r L
2
rp d p
rp y dp son el radio y el diámetro del capilar, respectivamente. Defina Sgv = Area superficial
por unidad de volumen de grano.
As n 2π rL nπ r 2 2 1
S gv = = = (2.25)
AL(1 − φ ) AL(1 − φ ) A r (1 − φ )
En la Ec. 2.25:
nπ r 2
φ= (2.26)
A
⎡ φ ⎤
S gv = S pv ⎢ ⎥ (2.27)
⎣1 − φ ⎦
4
Puesto que Spv = 2/rp, luego rp2 = 2
, entonces la Ec. (2.23) se convierte en:
S pv
⎛ r2 ⎞ ⎛ 4 / Spv 2 ⎞
k = ⎜⎜ ⎟⎟φ = ⎜⎜ ⎟⎟φ
⎝ 8 ⎠ ⎝ 8 ⎠
φ
k= (2.28)
2 S pv 2
1 φ3
k= (2.29)
2S gv 2 (1 − φ ) 2
2
⎛L ⎞ Camino real
τ =⎜ h ⎟ = (2.30)
⎝ L⎠ Longitud corazon
de donde;
Lh = L τ
⎛ nπ r 4 ⎞ ∆p
q=⎜ ⎟ (2.31)
⎝ 8µ ⎠ L τ
9/dg
6/dg
6/dg
5/dg
4/dg
Aplicando un procedimiento similar. En otras palabras combinando la Ec. (2.31) con la ley
de Darcy, se llega a:
⎛ r2 ⎞
k =⎜ p
⎜ 8τ ⎟⎟ φ (2.32)
⎝ ⎠
y;
φ
k= (2.33)
2τ S pv 2
1 φ3
k= (2.34)
2τ S gv 2 (1 − φ ) 2
Usando granos esféricos, Carman halló que τ puede ser aproximado a 2.5, luego las Ecs.
(2.33) y (2.34) se convierten en:
φ
k= (2.35)
5 S pv 2
1 φ3
k= (2.36)
5S gv 2 (1 − φ ) 2
φ 4.4
k = 0.135
S wi 2
Limitaciones:
2
⎛ 1 − S wi ⎞
k = (10φe ) ⎜ 4
⎟
⎝ S wi ⎠
REFERENCIAS
CAPITULO 3
PROPIEDADES DE LOS FLUIDOS DEL YACIMIENTO
Un análisis PVT es costoso y muchas veces se trata de yacimientos viejos que no poseen ésta
información o muy nuevos que todavía no han sido evaluados. Por éstas razones se han
desarrollado una serie de ecuaciones o Correlaciónes empíricas que permitan determinar las
propiedades de los fluidos del yacimiento. A continuación se presentan las principales
características del gas y el petróleo y las Correlaciónes más usadas en su determinación. En
general, el PVT se refiere al conjunto de propiedades de volumen medidas a una presión y
temperatura determinada. Estas propiedades son factor volumétrico de formación del
petróleo, que es función de la presión del yacimiento, del factor volumétrico y de la
compresibilidad del crudo, factores volumétricos del gas y el agua., gas disuelto en crudo
que es función de la gravedad del crudo, temperatura, presión y gravedad del gas. La
viscosidad del crudo es función de la Temperatura, presión y gas disuelto. El factor
volumétrico bifásico, β t = β o + ( Rsi − Rs ) β g . La compresibilidad del crudo es función de
P, API, T y γg. La compresibilidad del agua es función del gas disuelto en agua y la
concentración de cloruro de sodio dado en parte por millón, ppm.
NOTA: En los gráficos que se presentan a continuación, los puntos de arranque de cada
curva son muy diferentes y no se ven reflejados allí.
Volumen en yacimiento
βg =
Volumen en superficie
znRT
V =
p
PARAMETRO
µo Bo
Rs
Pb
PRESION
PARAMETRO
Co
Bt
Pb
PRESION
PARAMETRO
Bg
µg
Pb
PRESION
Bw
PARAMETRO
Rsw
Pb
PRESION
znRTyto p yto
βg =
nRTsup psup
⎛T ⎞
β g = cte ⋅ z ⎜ ⎟
⎝ p ⎠ yto
⎛T ⎞ ft 3
β g = 0.00282 ⋅ z ⎜ ⎟ ;
⎝ p ⎠ yto pcn
⎛T ⎞ bbl
β g = 0.00504 ⋅ z ⎜ ⎟ ;
⎝ p ⎠ yto pcn
⎛T ⎞ PCN
β g = 35.35 ⋅ z ⎜ ⎟ ; 3
⎝ p ⎠ yto ft
⎛T ⎞ pcn
β g = 198.4 ⋅ z ⎜ ⎟ ;
⎝ p ⎠ yto bbl
1 dV 1 d ρ
c=− =
V dP ρ dP
La compresibilidad del petróleo juega un papel importante entre Pb y P > Pb, de modo que:
p V dV
c ∫ dP = − ∫
Pb Vb V
Integrando resulta:
V
c( Pb − P ) = − ln
Vb
V
ec ( Pb − P ) =
Vb
Puesto que existe una relación de volúmenes, se pueden reemplazar éstos por los factores
volumétricos, Por lo tanto:
β o = β ob ec ( P − P ) b
Psr = P /Pcr
Tsr = T/Tcr
Con estos valores se va a la gráfica de Standing y Katz y se lee el valor de Z. Para una
mezcla de gases, las propiedades críticas del gas y su peso molecular promedio se obtienen
por ponderación con la composición molar, de modo que:
n
PM = ∑ PM i ⋅ xi
i =1
n
Pcr = ∑ Pcri ⋅ xi
i =1
n
Tcr = ∑ Tcr i ⋅ xi
i =1
pV = znRT
nPM
ρ=
V
nPM PM ⋅ p
ρ= =
znRT / p zRT
Puesto que la gravedad específica del gas es la relación entre la densidad del gas y la
densidad del aire, se tiene:
ct = co So + cw S w + cg S g + c f
ct
ce =
1 − Sw
1 dV 1 N − Nr 1 Np
ce = − = =
V dp N ∆p N ∆p
Es un proceso similar al anterior, la diferencia radica en que a medida que miden los
volúmenes, el gas es retirado de la muestra en cada decremento de presión.
Se denota como Pb. Es la presión a la cual la primer burbuja de gas comienza a liberarse del
petróleo. También es llamada presión de saturación. Cada yacimiento tiene su presión de
burbuja particular. La presión del punto de burbuja se determina en función de la temperatura,
la gravedad específica del gas, γg, la gravedad específica del petróleo, γo, y la cantidad de gas
disuelto en el crudo, Rs. La correlación de Standing1,2,3 fue obtenida en 105 experimentos,
utilizando petróleos de California. El promedio de error de la correlación fue aplicado a los
datos usados para desarrollar el método y resultó ser 4,8 % a una presión de 106 psi. El rango
de datos utilizados para desarrollar el método se da en la tabla 3.1. Los gases involucrados en
el desarrollo del método no contenían Hidrógeno ni Sulfuro de Hidrógeno. Alguno de los
gases contenían Dióxido de Carbono en cantidades insignificantes. La correlación podrá
aplicarse a otros petróleos teniendo en cuenta que las características de los petróleos y el gas
sean similares a las usadas para desarrollar el método. Fuera de éste rango se incurre en un
margen de error. La correlación de Standing, matemáticamente se expresa como:
0.83
⎛R ⎞
Pb = 18 ⎜ sb ⎟ 10 yg
⎜γ ⎟
⎝ g ⎠
Esta ecuación graficada en papel doblemente logarítmico es una línea recta. El rango de
aplicación de la correlación de Standing es:
Parámetro Rango
Pb 130 a 7000 psi
TF 100 a 258 °F
Rsb 20 a 1.425 pcn/BF
API 16.5 a 63.8
γg 0.59 a 0.95
Existen otras dos Correlaciónes para calcular la presión del punto de burbuja, a saber:
1/ b
⎡ GOR ⎤
Pb = ⎢ c API / (T + 460) ⎥
⎢⎣ a γ g e ⎥⎦
b 1.937 1.870
c 25.240 23.310
Parámetro Rango
Pb 50 a 5250 psi
TF 70 a 295 °F
GOR 20 a 2070 pcn/BF
API 16 a 50 °API
γg 0.56 a 1.18
Base de datos: 600 análisis PVT de crudos de todo el mundo. Un total de 6004 datos
medidos.
Correlación de Lasater3
Rb
A=
379.3
A
Yg =
A + 350γ o / Mo
Para API ≤ 40
Mo = 630 − 10 API
Si Yg ≤ 0.6
Si Yg > 0.6
Parámetro Rango
Pb 48 a 5780 psi
TF 82 a 272 °F
API 17.9 a 51.1 °API
γg 0.574 a 1.223
Rsb 3 a 2905 pcn/BF
Base de datos: 137 sistemas independientes de crudos del Canadá, Estados Unidos y Sur
América. Un total de 158 presiones en el punto de burbuja medidas experimentalmente.
Error promedio 3.8%
Correlación de Glaso5
Pb = 10
0.816
⎛R ⎞ ⎛ TF 0.172 ⎞
F = ⎜ sb ⎟⎟ ⎜ 0.989 ⎟
⎜ γg ⎝ API ⎠
⎝ ⎠
0.816
⎛R ⎞ ⎛ TF 0.130 ⎞
Pb = ⎜ ⎟⎟ ⎜ 0.989 ⎟
⎜γg ⎝ API ⎠
⎝ ⎠
Parámetro Rango
Pb 165 – 7142 psi
TF 80 – 280 °F
API 22.3 - 48.1 °API
γg 0.650 - 1.276
Rsb 90 - 263 pcn/BF
Base de datos: Datos de 45 muestras de aceite de crudos del Mar del Norte, Medio Oriente,
Algeria y Estados Unidos. Error promedio 1.28 %.
b
⎛R ⎞ ⎛ 10C *TF ⎞
Pb = a ⎜ sb ⎟⎟ ⎜⎜ d * API ⎟⎟
⎜ γg
⎝ ⎠ ⎝γg ⎠
°API a b c d
≤10 12.8470 0.9636 0.000993 0.034170
10-35 25.2755 0.7617 0.000835 0.011292
35-45 216.4711 0.6922 -0.000427 0.023140
Base de datos: 289 análisis PVT del Area Mayor de Oficina, en el Oriente de Venezuela.
Correlación de Al-Marhoun6
Parámetro Rango
Pb 130 - 3513 psi
TF 74 – 240 °F
API ---------
γg 0.752 - 1.367
Rsb 20 - 1602 pcn/BF
Correlación de Kartoatmodjo6,13
Para ºAPI ≤ 30
0.9086
⎛ Rsf ⎞
Pb = ⎜ ⎟⎟
⎜ 0.05958γ g100 ×1013.1405 API /(TF + 460)
0.7972
⎝ ⎠
0.9143
⎛ Rsf ⎞
Pb = ⎜ ⎟⎟
⎜ 0.03150 γ +
g 100 × 10
0.7587 11.2895 API /( TF 460)
⎝ ⎠
Parámetro Rango
Pb 14.7 - 6054.7 psi
TF 75 - 320 °F
API 14.4 - 58.9 °API
γg100 0.752 - 1.367
Rsf 0.0 - 2890 pcn/BF
Base de datos: Reportes PVT del Sur de Asia, Norte América, Medio Oriente y
Latinoamérica, en total 5392 datos de 740 crudos diferentes. Error promedio relativo de
3.34 %.
Parámetro Rango
Pb 31.7 - 4137.7 psi
TF 95 -260 °F
Temp. Tanque 71-108 °F
API 14.4 - 58.9 °API
γg 0.657-1.731
Rs 6-1645 pcn/BF
Bob 1.060 – 2.064 rb/STB
Parámetro Rango
Pb 549.7-3789.7 psi
TF 95-256 °F
Temp. Tanque 68-108 °F
API 14.3-45.1 °API
γg en separador 0.573-1.337
Rs 6-1645 pcn/BF
Bob 1.045 – 1.995rb/STB
Presión de separación 34.7-514.7 psi
Rs en el separador 66-1571 pcn/BF
⎣ ⎦
Parámetro Rango
Pb 3190 - 4058.7 psi
TF 162 - 260 °F
API 40.3-46.5 °API
γg 0.657-1.731
Rs total 961.3 - 1400.1 pcn/BF
Bob 1.569-1.905 rb/STB
Glaso5 reportó factores correctivos para la Presión en el punto de burbuja, cuando están
presentes el dióxido de carbono, el nitrógeno y el sulfuro de hidrógeno.
( )
2
+ ⎡⎣1.954 × 10−11 API 4.699TF + (0.027 API − 0.36666) ⎤⎦ YN 2
Se denota como Rs. También se le denomina solubilidad del gas en petróleo, razón gas
disuelto y relación gas petróleo, RGP (en inglés GOR). Se define como la cantidad de gas
medido a condiciones de superficie, que se disuelven en un barril de petróleo, también
medido a condiciones de superficie. Los factores que afectan la solubilidad del gas en el
petróleo, Rs, son:
1.204
⎛ P ⎞
Rrb = γ g ⎜ Yg ⎟
⎝ 18 ×10 ⎠
Parámetro Rango
Pb 130 - 7000 psi
TF 100 - 258 °F
API 16.5 - 63.8 °API
γg 0.590 - 0.95
Rsb 20 - 1425 pcn/BF
⎡ ⎡ Ps ⎤ ⎤
γ gc = γ g ⎢1 + 5.91× 10−5 API × Ts log ⎢ ⎥⎥
⎣ ⎣114.7 ⎦ ⎦
Base de datos: 600 análisis PVT de crudos de campos de todo el mundo. Un total de 6004
datos medidos. Error promedio de 0.7 %.
Parámetro Rango
Pb 50-5250 psi
TF 70 - 295 °F
API 16 - 58 °API
γg 0.56 - 1.18
Rsb 20-2070 pcn/BF
Correlación de Lasater3
γg
A= P
TF
Para A ≥ 3.29
Yg = (0.121A − 0.236)0.281
γ gγ o
Rs = 132.755
Mo(1 − Yg )
Mo, se obtiene de la misma forma que para el cálculo de la presión del punto de burbuja
usando la correlación de Lasater. Base de datos: 137 sistemas independientes del Canáda,
Estados Unidos y Sur América.
Parámetro Rango
Pb 48 - 5780 psi
TF 82 - 275 °F
API 17.9 – 51.1 °API
γg 0.574 - 1.123
Rsb 3 – 2090 pcn/BF
Correlación de Glaso5
1.2255
⎡ API 0.989 ⎤
RSb = γ g * ⎢ F * 0.172 ⎥
⎣ T ⎦
Parámetro Rango
Pb 165-7142 psi
TF 82 - 275 °F
API 22.3-48.1 °API
γg 0.650-1.276
Rsb 90-2637 pcn/BF
Base de datos: 45 muestras de aceite del Mar del Norte, Medio Oriente, Algeria y Estados
Unidos. Error promedio 1.28 %.
d
⎛ 10( a* API −b*T ⎞
Rsb = γ g ⎜ Pb ⎟
⎝ c ⎠
°API a b c d
≤10 0.30405 0.000000 12.2651 0.9699
10-35 0.015200 0.4484x10-3 15.0057 1.0950
35-45 0.024800 -0.001469 112.9250 1.1290
Base de datos: 289 análisis PVT de Area Mayor de Oficina, en el Oriente de Venezuela.
Correlación de Al-Marhoun6
e
Rs = ⎡⎣ aγ g bγ o cTF d * P ⎤⎦
Parámetro Rango
TF 74 - 240 °F
Rs 26 - pc 1602 pcn/BF
Base de datos: 69 análisis PVT de muestras de fondo de pozo de crudos del Medio Oriente.
Correlación deKartoatmodjo5
Para ºAPI ≤ 30
0.7972
Rsf = 0.05958 γ g100 P 1.0014 × 1013.1405 API /(TF + 460)
Donde:
⎛ Ps ⎞
γ g100 = γ gsep × ⎜ 1 + 0.1595 API 0.4078 Ts−0.2466 log
⎝ 114.7 ⎟⎠
Parámetro Rango
TF 75 - 320 °F
API 14.4 - 58.9 °API
γg100 0.752 - 1.367
γgsep 0.379 - 1.709
Rsf 0.0 - 2890 pcn/BF
Ts 65 – 186 °F
Ps 14.7 - 514.7 psia
Base de datos: Reportes PVT del Sudeste de Asia, Norte América, Medio Oriente y
Latinoamérica. Fueron usados 5392 datos representados en 740 diferentes crudos. Error
promedio -4.6839%.
Rs = 0.01936 Pb1.1574γ gc
0.73495
× 100.000337TF + 0.01771 API
Se denota como co. Cuando la presión es mayor que la presión del punto de burbuja, el
petróleo en el yacimiento tiene todo el gas en solución. Cuando se aplica presión en exceso a
éste sistema, el líquido sufre una disminución no lineal en su volumen que depende de la
temperatura y composición del petróleo. Esa pequeña variación en el volumen es lo que se
conoce como factor de compresibilidad del petróleo, que es muy significativa en cálculos de
Ingeniería de Yacimientos aunque en las demás ciencias se desprecia tal factor debido a la
poca compresibilidad (en algunos casos ninguna) de los líquidos. La correlación de Vásquez
y Beggs es un medio excelente para obtener valor de la compresibilidad del petróleo. Vásquez
y Beggs usaron aproximadamente 2.000 medidas experimentales de compresibilidad del
petróleos en más de 600 crudos diferentes para desarrollar una correlación en función de Rsb,
TF, γg, API y presión. La correlación obtenida por Vásquez y Beggs3 fue:
Correlación de Ahmed5
1
co = e −0.00018473 P
24841.0822 + 14.07428745 Rs
Correlación de Kartoatmodjo5,13
10−6 0.5002
co = 6.8257 Rsf API 0.3613TF 0.76606γ g0.35505
100
P
Parámetro Rango
TF 75 - 320 °F
API 14.4 - 58.9 °API
Rsf 0.0 - 2890 pcn/BF
Base de datos: Reportes PVT del Sureste de Asia, Norte América, Medio Oriente y
Latinoamérica. Fueron usados 5392 datos representados en 740 diferentes crudos. Error
promedio: 0.30078 %.
Se denota por Bo o ßo. Se define como el volumen de petróleo (más su gas en solución) en el
yacimiento, requerido para producir un barril de petróleo medido a condiciones de superficie.
Por ejemplo, ßo = 1.5 B/BF significa que para tener un barril de petróleo en superficie (Barril
Fiscal, BF) se requiere 1.5 barriles de petróleo en el yacimiento. Lógicamente, el valor de ßo
será mayor de la unidad debido al gas que entra en solución. De otro modo, al pasar el
petróleo de yacimiento a superficie sufre disminución en presión y temperatura y ocurre
liberación de gas presente en el líquido (petróleo). Este proceso conduce a una merma del
volumen de petróleo del yacimiento al pasar a superficie. El ßo se calcula en función de Rs,
API, γg y temperatura. La correlación de Standing3 para determinar ßo en petróleos
saturados, matemáticamente se expresa como:
⎧⎪ γg ⎫⎪
β o = 0.972 + 0.000147 ⎨ Rs + 1.25T ⎬
⎩⎪ γo ⎭⎪
Correlación de Arps4
⎡ API ⎤ ⎡ API ⎤
Bob = 1 + aRsb + b (TF − 60 ) ⎢ ⎥ + cRsb (T − 60 ) ⎢ ⎥
⎢⎣ γ gc ⎥⎦ ⎢⎣ γ gc ⎥⎦
Parámetro Rango
Bob 1.042 - 1.545 rb/STB
TF 70 - 295 °F
API 16 - 58 °API
γgc 0.590 - 0.950
Rsb 20-2070 pcn/BF
Correlación de Glaso5
Bob = 1 + 10
Parámetro Rango
Bob 1.025 - 2.588 rb/STB
TF 80 - 280 °F
API 16 - 58 °API
γo 0.788 - 0.920
γg 0.650 - 1.276
Rsb 90 - 2637 pcn/BF
Base de datos: Datos de 45 muestras de aceites de crudos del Mar del Norte, Medio
Oriente, Algeria y Estados Unidos. Error promedio -0.43%
API API
Bob = A + B × Rsb + C (TF − 60) + D(T − 60) Rsb
γg γg
Base de datos: 289 análisis PVT del Area Mayor de oficina, en el Oriente de Venezuela.
Correlación de Al-Marhoun6
F = Rsb0.742390γ 0.323294
g
γ o−1.202040
Parámetro Rango
Bob 1.032 – 1.997 rb/STB
TF 74 - 240 °F
API 16 - 58 °API
γo 0.788 - 0.920
γg 0.752 - 1.367
Rsb 26 - 1602 pcn/BF
Base de datos: 69 análisis PVT de muestras de fondo de 69 yacimientos de crudo del Medio
Oriente. Error promedio de 0.03 %.
Donde:
F = Rsf0.755γ g0.25
100
γ o−1.5 + 0.45TF
⎛ ⎛ Ps ⎞ ⎞
γ g100 = γ gsep ⎜ 1 + 0.1595 API 0.4078Ts−0.2466 log ⎜ ⎟⎟
⎝ ⎝ 114.7 ⎠ ⎠
Parámetro Rango
Bof 1.007 - 2.747 rb/STB
TF 75 - 320 °F
API 14.4 - 58.9 °API
γg sep 0.379 - 1.709
Rsf 0 - 2890 pcn/BF
Ts 65 – 186 °F
Ps 14.7 - 514.7 psia
Bof es el aceite requerido a una presión de yacimiento menor que la Pb para producir 1 STB
de aceite a 60 ºF, cuando se somete a liberación flash a través del separador a condiciones
de de tanque, Rsf es el gas en el separador y en el tanque a una presión menor que la Pb en 1
STB de aceite.
Base de datos: Reportes PVT del Sudeste de Asia, Norte América, Medio Oriente y
Latinoamérica. Fueron usados 5392 datos representados en 740 diferentes crudos. Error
promedio: -0.1040%
Bo = 1 + 10 (
−2.6541 + 0.5576log( A ) + 0.3331*[ log( A )])
Donde:
A = Rs0.5956γ gc γo
0.2369 −1.3282
+ 0.0976TF
Bo = 1 + 3 ×10−5 Rsc1.3267 + A
Donde:
β o = β ob ec ( P − P )
o b
Correlación de Vazquez3
P
− A ln
Bo = Bob e Pb
Parámetro Rango
Bob 1.025-2.588 rb/STB
TF 70 - 295 °F
API 16 - 58 °API
γg 0.560-1.180
Rs 20-2070 pcn/BF
Pb 50-550 psi
Base de datos: 600 análisis PVT de crudos de todo el mundo. Un total de 6004 datos
medidos.
Correlación de Calhoun8
Esta correlación se puede usar para calcular la densidad del aceite a presiones P < Pb,
haciendo Bob igual Bo y Rsb igual a Rs.
1
ρo = ρob
1 − Mo(1 − ( P / Pb ))
M o = 10a Rsep
b
PscTF d 10− e (TF )
Base de datos: Se utilizaron un total de 177 PVT disponibles de crudos colombianos. Error
absoluto maximo: 3.15%
ρo = ρob ec ( P − P )
o b
P
A ln
ρo = ρob e Pb
Base de datos: 600 análisis PVT de crudos de campos de todo el mundo. Un total de 6004
datos medidos.
Correlación de Ahmed3
(
B e-0.00018473P − e-0.00018473BP )
ρo = ρob e
B = −(4.588893 + 0.025999 Rs ) −1
β t = β o + ( Rsi − Rs )β g
Correlación de Glaso5
2
Bt = 100.080135+ 0.47257log( F ) + 0.17351log( F )
⎡ TF ⎤ ( −0.00027*Rs )
F = Rs ⎢ 0.3 ⎥ γ o2.9×10 P −1.1089
⎢⎣ γ g ⎥⎦
Parámetro Rango
TF 80 -280 °F
γg 0.650-1.276
Rs 90 – 263 pcn/BF
Base de datos: 45 muestras de crudos del Mar del Norte, Medio Oriente, Algeria y Estados
Unidos. Error promedio de 6.54 %.
Correlación de Al-Marhoun6
Parámetro Rango
TF 70 - 240 °F
γg 0.752 - 1.367
γo 0.804 - 0938
Rs 20 - 1602 pcn/BF
Base de datos: 69 análisis PVT de muestras de fondo de 69 yacimientos de crudo del Medio
Oriente. Error absoluto promedio de 4.11 %.
Se denota como µo. Se define como la medida de la resistencia del petróleo al flujo. Es
usualmente medida en centipoises, cp, (gr/cm*seg). La resistencia al flujo es causada por
fricción interna generada cuando las moléculas del fluido tratan de desplazarse unas sobre
otras. Los valores de µo se requieren a diferentes presiones, tanto en Ingeniería de
Yacimientos como en Ingeniería de Producción. Si se dispone de un análisis PVT las medidas
de la viscosidad se reportan a la presión y temperatura del yacimiento y a diferentes
presiones. No obstante, a medida que el fluido fluye a superficie su temperatura disminuye, lo
que involucra la necesidad de corregir la viscosidad para cambios de temperatura, mediante
Correlaciónes empíricas. Los principales factores de interés en Ingeniería de Petróleos que
afectan µo son: La composición del petróleo, la temperatura, el gas disuelto y la presión. La µo
aumenta cuando disminuye la API y también aumenta con un decremento en la temperatura.
El efecto del gas disuelto es alivianar el petróleo y por tanto disminuir su viscosidad. Mientras
exista un incremento en la presión sobre un petróleo subsaturado, su viscosidad se
incrementará. El método más común de obtener la viscosidad del petróleo, µo, consiste en
evaluar la propiedad para petróleo muerto (petróleo sin gas disuelto) para luego corregirla por
efectos del gas disuelto. Para sistemas saturados se utilizan la correlación de Beggs y
Robinson1,2, desarrollada con más de 2000 mediciones de viscosidad, usando 600 muestras
distintas de petróleo. El rango de variables de los datos usados se da a continuación:
Variable Rango
Rs 20 a 2070 pcn/BF
API 16 a 58
P 0 a 5250 psi
TF 70 a 245 °F
La correlación de Beggs y Robinson1 para determinar la viscosidad del petróleo muerto, µod,
es:
µod = 10 X − 1
X = T −1.163e6.9824−0.04658 API
La ecuación para corregir la viscosidad por efecto del gas disuelto es:
µob = A µod B
B = 5.44( Rs + 150)−0.338
µo = µob (P/Pb)m
µob = A [ µod ]
B
Variable Rango
µod 0.377-50 cp
Rsb 51-3544 pcn/BF
API 16 a 58
Pb 132 a 564 psi
TF 72 a 297 °F
Base de datos: 457 muestras de aceite, de crudos de Canadá, Estados Unidos y Sur
América.
Correlación de Khan
0.09 γ g
µob = 4.5
⎛ TF ⎞
(1 − γ o )
3
Rs0.333 ⎜ ⎟
⎝ 460 ⎠
Variable Rango
Rs 24 - 1901 pcn/BF
API 14.3 – 44.6
P 100 – 4315 psi
TF 75 – 240 °F
µob = A [ µod ]
B
−7
Rsb − 7.4×10−4
A = 10 Rsb (2.2×10
Correlación de Khan9
−4
µo = µob ( P / Pb ) −0.14 e−2.5×10 *( P − Pb )
Base de datos: 1503 datos experimentales de crudos de Arabia Saudita. Error promedio del
2 %.
Correlación de Kartoatmodjo5
Donde:
−0.001127*Rsf
f = (0.2001 + 0.8428 × 10 ) * µod (0.43+0.5165 y )
−0.00081Rsf
y = 10
Variable Rango
µo 0.1 – 586 cp
µod 0.5 – 682 cp
Rsf 0 – 289 rb/STB
Base de datos: Se usaron 5392 reportes PVT del Sureste de Asia, Norte de América, Medio
Oriente y Latinoamérica. Error promedio de 0.08034%
µo = A( µod ) x
−4
A = 0.1651 + 0.6165 ×10−6.0866×10 Rs
−3
B = 0.5131 + 0.5109 ×10−1.183×10 Rs
Base de datos: Un total de 126 análisis PVT de el Golfo de México fueron utilizados para
el desarrollo de las Correlaciónes. Error relativo promedio de -3.12 %.
Correlación de Beal1
Variable Rango
µo 0.160 - 315 cp
µob 0.142 - 127 cp
P 1515 – 5515 psi
Correlación de Khan9
−5
µo = µob e9.6×10 *( P − Pb )
Correlación de Kartoatmodjo5
Variable Rango
µo 0.2 - 517 cp
µob 0.2 – 184.8 cp
Pb 14.7 – 6054.7 psi
Base de datos: Reportes PVT del Sureste de Asia, Norte de América, Medio Oriente y
Latinoamérica. Fueron usados 5392 datos representados en 740 diferentes crudos. Error
promedio -4.2870%
A = −1.0146 + 1.3322 log( µob ) − 0.4876 [ log( µob ) ] − 1.15036 ⎡⎣log ( µob ) ⎤⎦ )
2 3
Base de datos: Un total de 126 datos experimentales del Golfo de México fueron utilizados,
para el desarrollo de las Correlaciónes. Error relativo promedio: -0.19%
Se denota como γg. La gravedad específica de un gas, se define como la relación entre su
densidad y la densidad del aire. En el campo ésta propiedad se mide fácilmente en la cabeza
del pozo, por lo tanto, es usada como una indicación de la composición del gas. La gravedad
específica del gas es proporcional a su peso molecular, PM, cuando éste se mide a bajas
presiones donde el comportamiento del gas se aproxima a la idealidad. Una vez se obtiene el
peso molecular, la gravedad del gas se determina dividiendo su peso molecular entre 28.966
(peso molecular del aire). La gravedad específica también se usa para Correlaciónar otras
propiedades físicas de los gases como las propiedades críticas. En algunas ocasiones cuando
existe la presión y temperatura en el separador, es necesario corregir la gravedad específica
del gas para tener unos datos de PVT óptimos. La correlación de Vázquez y Beggs1,3 permite
efectuar éste proceso:
⎛ Ps ⎞
γ gs = γ g ⎜1 + 2.215 × 10−5 API Ts ln ⎟
⎝ 114.7 ⎠
Correlacion de Katz11
Base de datos: 76 análisis PVT del Area Mayor de oficina, en el Oriente de Venezuela.
Error promedio: El 85% de los 76 valores de γg calculados con la correlación, presentan un
error menor del 20%
Para condensados1
γ g1 + γ g 2
γg =
2
Las propiedades críticas están sujetas a variaciones por presencia de contaminantes (Dióxido
carbónico, CO2 y Sulfuro de Hidrógeno, H2S). La correlación de Wichert y Aziz1 es utilizada
para efectuar éstas correcciones:
[
e = 120 YCO 2 + YH 2 S − (YCO 2 + YH 2 S )
1.6
] + 15( YH 2 S − YH 2 S
4
)
Método de Carr-Kobayashi-Burrows1
P
Psr =
(677 + 15γ g − 37.7γ g2 )
T
Tsr =
(168 + 325γ g − 12.5γ g2 )
Generalmente, las propiedades críticas y la gravedad específica del gas se desconocen, pero,
casi siempre se da la composición de la mezcla de gases, es decir, se lista cada uno de los
componentes del gas con su respectivo porcentaje o fracción volumétrica de la cantidad
existente en la mezcla. A continuación se da un ejemplo de un reporte típico de una mezcla de
gases.
Para evaluar las propiedades críticas de la mezcla de gases se toma la presión crítica, la
temperatura crítica tabulados a continuación, cada una de las cuales se multiplica por su
fracción volumétrica y la sumatoria constituye la presión crítica, Pcr, y la temperatura crítica,
Tcr, respectivamente.
n
Pcr = ∑ Pcr i Yi
i =1
n
Tcr = ∑ Tcr i Yi
i =1
Para evaluar la gravedad específica de la mezcla gaseosa se lee de la tabla anterior los pesos
moleculares de cada compuesto presente en la mezcla, cada uno de ellos se multiplica por la
fracción volumétrica, Yi, y su sumatoria se divide entre 28,966.
n
1
γg= ∑ PM i * Yi
28,966 i =1
Los Heptanos y compuestos más pesados (C7+) no fueron tabulados anteriormente debido a
que no tiene una composición ni peso molecular fijo. Por tal razón, las propiedades críticas se
determinan por medio de la correlación gráfica en función del peso molecular y la gravedad
específica de éste grupo de compuestos. La correlación gráfica mostrada a continuación
permite obtener estos parámetros a partir de su peso molecular y su gravedad específica.
1400
0,95
1300
1200
0,85
1100
0,80
0,75
1000 0,70
500
PRESION CRITICA, lpc
0.85
0.80
400 0.95
0.70
0.75
300 0.90
200
100 120 140 160 180 200 220 240
PESO MOLECULAR
Las Correlaciónes de Sutton and Whitson para las propiedades críticas de los C7+ son:
Tcr = 341.7 + 811γ C7+ + (0.4244 + 0.1174γ C7+ )TB + (0.4669 − 3.2623γ C7+ )105 / TB
( )
3
TB = 4.5579 PM C7+ 0.15178γ C7+ 0.15427
En la realidad no existen gases ideales o perfectos; sin embargo, muchos gases cerca de la
temperatura y presión atmosféricas se aproximan a la idealidad. El gas ideal puede definirse
como el gas cuyo volumen se reduce a la mitad al duplicarse la presión y cuya presión se
dobla si se duplica la temperatura manteniendo el volumen constante. Esto no es más que los
enunciados de la leyes de Boyle y Charles Mariotte. En muchos gases en particular los gases
naturales de interés para los Ingenieros de Petróleos, se ha observado que si el volumen del
gas se comprime a la mitad, la presión resulta ser menor del doble de la presión inicial; es
decir, el gas es más compresible que el gas ideal.
Psr6
D = 0.32
109(Tsr −1)
E=B+C+D
Z = A + (1 − A)e − E + F * PsrG
f (ρ r )
ρ r (i +1) = ρ r (i ) −
g (ρ r )
2
g (ρ r ) = 6 A ρ r3 + 3B ρ r2 + 2 C ρ r + D + E ρ r2 * (3 + F ρ r (3 − 2F ρ r2 ))e ( − F ρ r )
A = 0.06424
D = Tsr
E = 0.6816/Tsr²
F = 0.6845
G = 0.27 Psr
Correlación de Beggs y Brill4. Este método radica en un ajuste efectuado sobre una de las
curvas de la gráfica de Standing y Katz.
D = 0.0657/(Tsr - 0.86)
E = 0.32−19.53(Tsr −1)
F = 0.122−11.53(Tsr −1)
Debe tenerse en cuenta, que si los exponentes de e son menores de 100, todo este valor se
anula.
(1 + X + X 2 + X 3 ) + A X + B X C
Z=
(1 − X ) 3
Xi = b ρr /4
Método de Burnett4
⎡ 90 Psr ⎤
Z = 1 + ( Z ´−1) ⎢sin ⎥
⎣ Psr ' ⎦
Donde:
T
Tsr =
(315.1G + 170.4)
P
Psr =
(697.2 − 43G )
⎡ P ⎤ e Psr / Tsr
G = ⎢1.1 + 0.26Tsr + (1.04 − 1.42Tsr ) sr ⎥
⎣ Psr ' ⎦ Tsr
Método de Hakinson-Thomas-Phillips4
⎡ A ⎤ ⎡ Ppr ⎤ ⎡ Psr2 ⎤
Z − 1 + ⎢ A 4 *Tsr − A2 − 62 ⎥ * ⎢ 2 ⎥ + ( A3 × Tsr − A1 ) * ⎢ 3 3⎥
+
⎣ Tsr ⎦ ⎣ Z × Tsr2 ⎦ ⎣ Z × Tr ⎦
− A P2
⎡ A1 × A5 × A7 × Psr5 ⎤ ⎛ A8 × Psr2 ⎞ Z 2 ×Tsr2
1 sr
⎢ ⎥⎜ 1 + 2 ⎟
e =0
⎣ Z 6
× T 6
sr ⎦⎝ Z 2
× Tsr ⎠
Esta correlación requiere la densidad reducida obtenida iterativamente para la solución del
factor de compresibilidad del gas. Mediante la correlación de Dranchuk et al, se obtiene la
compresibilidad del gas.
Cg = CR / Pcr
2
dZ 5 A ρ r5 + 2 B ρ r2 + C ρ r + 2 E ρ r2 (1 + F ρ r2 − F 2 ρ r4 )e( − F ρr )
=
d ρr ρ r Tsr
1 1
CR =
Psr (1 + ρ r ) dZ
Z d ρr
Se denota como Bg o ßg. Es un factor que relaciona el número de barriles de gas que se
necesitan en el yacimiento para producir un pie cúbico de gas a condiciones normales (pcn).
Las ecuaciones para determinar ßg se presentaron anteriormente.
Se denota como µg. Es una propiedad importante para determinar la resistencia al flujo que
presenta el gas durante su producción y transporte. Generalmente, la viscosidad del gas
aumenta con los incrementos de presión. A presiones bajas la viscosidad del gas (al contrario
que los líquidos) se incrementa con la temperatura. Sus unidades son los centipoises, cp. la
viscosidad del gas también se determina por medio de Correlaciónes empíricas con un
margen de error inferior al 2 %. La correlación de Carr et al permite obtener la viscosidad del
gas a la presión atmosférica, µ1. Este valor, es llevado a cualquier presión, mediante la
correlación de Dempsey.
Correlación de Dempsey1,11
⎛ µg ⎞
⎟ = A0 + A1 Psr + A2 Psr + Tsr ( A4 + A5 Psr + A6 Psr + A7 Psr ) + ξ
2 2 3
ln ⎜ Tsr
⎝ µ1 ⎠
ξ = Tsr 2 ( A8 + A9 Psr + A10 Psr 2 + A11 Psr 3 ) + Tsr 3 ( A12 + A13 Psr + A14 Psr 2 + A15 Psr 3 )
ßw = (A + B P + C P²) SALCOR
SALCOR = 5.1×10−8 P + (5.47 ×10−6 −1.95×10−10 P)(T − 60) + (−3.23×10−8 + 0.0085 P)(T − 60)Nacl +1
247.8
−140
µ w = PRESCOR( SALCOR)0.02414 × 10 TK
SALCOR = (1 − 0.00187 Nacl + 0.000218 Nacl 2.5 + T − 0.0135T )(0.0827 Nacl − 0.00034 Nacl 1.5 )
Cuando se desconoce la gravedad específica del agua, pero se conoce el contenido de sólidos
disueltos en mg/lt, SDC, la gravedad específica del agua puede calcularse con:
γw = 1 + 0.695x10-6 SDC
FCOR = 0.99984x10-3.229268012P
FCOR es el factor de corrección por efecto del gas disuelto. Se obtiene por el ajuste
estadístico de la correlación gráfica de Baker y Swerdloff.
σ280 = 0.0063 P + 53
1200 2500
Pi=500 psi
Pi=1000 psi
Pi=1500 psi
P Pi=2000 psi
1000 Pi=2500 psi
2000
Presión, psi
1500
600
1000
3
400
500
200
GOR
0 0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Fig. 3.3. Efecto de presión sobre el recobro, Presión de abandono de 14.7 psi12
1200 2500
Pi=500 psi
Pi=1000 psi
Pi=1500 psi
P Pi=2000 psi
1000 Pi=2500 psi
2000
1500
600
1000
3
400
500
200
GOR
0 0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Fig. 3.4. Efecto de presión sobre el recobro, Presión de abandono de 100 psi12
26
µ = 0.69 cp
12
µ = 1.38 cp
µ = 2.76 cp
P µ = 5.52 cp 21
16
8
Presión, psi
11
6
6
4
3
1
2
GOR
0 -4
0 2 4 6 8 10 12 14 16 18 20
3) Gas en solución. A menor gas en solución mayor será el recobro de crudo. El crudo del
yacimiento que contiene menor gas disuelto requiere una mayor contracción para tener
una saturación de gas crítica además existe una mayor caída de presión en las primeras
etapas de la vida del yacimiento. El petróleo con mayor solubilidad se crea un vacío
para permitir la formación de la saturación de gas crítica.
5) Saturación de agua connata. Cuando hay agua connata hay más recuperación ya que
las curvas de permeabilidades relativas se desvían hacia la región de saturación baja de
petróleo. En otras palabras, sin agua connata, el agua fluye más rápido. Ver Figs. 3.7 y
3.8.
15
25 Solub. Gas = 267 cu. ft/bbl
Solub. Gas = 534 cu. ft/bbl
13
P
15
10 5
3
3
GOR 1
0 -1
0 2 4 6 8 10 12 14 16 18 20
6
25 Sgc = 10
Sgc = 0, Swc = 0.3
20
Presión x 1000, psi
15
10
3
5
1
GOR
0 0
0 2 4 6 8 10 12 14
16
25 Swc = 0
Swc = 0.3
14
P
10
15
10 6
3
4
GOR
0 0
0 2 4 6 8 10 12 14
m=0.5
1500
Presión, psi
m=0.1
m=0
1000
500
0
0 5 10 15 20 25 30 35 40 45 50
1500
m=1.0
1000 `
m=0.5
m=0.1
500
m=0
0
0 5 10 15 20 25 30 35 40 45 50
1500
100 %
Presión, psi
90 %
1000
70 %
500 50 %
0%
0
0 10 20 30 40 50 60
7000
6000
100 %
Relación Gas-Petróleo, cu. ft/SCF
5000
4000
3000
2000
90 %
1000
70 %
0% 50 %
0
0 10 20 30 40 50 60
2.0 120
Equilibrio
Desequilibrio
1.8
100
1.6
1.4
80
1.2
1.0 60
0.8
40
0.6
0.4
20
0.2
0.0 0
0 500 1000 1500 2000 2500
Presión, psi
1200 3000
Equilibrio
Desequilibrio
1000 2500
600 1500
400 1000
200 500
GOR
0 0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
10) Ratas de producción de los fluidos y caída de presión. Un alto caudal trae como
consecuencia una alta presión diferencial lo que impide el establecimiento de equilibrio
en la vaporización quedando el petróleo sobresaturado de gas lo que causa histéresis de
vaporización y la consecuente pérdida de recobro.
siendo la saturación de líquido mayor en la vecindad del pozo productor. Si este gradiente
es muy alto perjudica el recobro de petróleo.
11) Volatilidad del crudo en el yacimiento. Cuando el crudo es de alta volatilidad, las
proporciones relativas de líquido y vapor no son las mismas a condiciones de superficie
y yacimiento. A mayor volatilidad menor recobro por la formación e bancos de gas.
REFERENCIAS
CAPITULO 4
El cálculo volumétrico de petróleo y/o gas es de una de las herramientas para la estimación
de reservas. Los métodos para cuantificar reservas son:
a) Método volumétrico
b) Ecuación de balance de materia
c) Curvas de declinación
d) Simulación numérica y/o matemática de yacimientos
7758 Ah φ Soi
N=
β oi
G = 43560 Ah φ S gi / Bgi
G está dado en pcn, ßgi está dado en bbl/pcn y h está dado en pies (intervalo de gas)
Yacimientos de condensados
0.2
⎛R ⎞
G p * = −2229.4 + 148.43 ⎜ pi ⎟ + 124130 / T + 21.831API
⎝ 100 ⎠
Las unidades de Np* están en STB/rb volumen poroso de hidrocarburos (HCPV). Válido
para presiones de abandono de 500 psia.
FR es el factor de recobro
k en permeabilidad absoluta en milidarcies
Pa es la presión de abandono del yacimiento en psia
T temperatura del yacimiento, °R
Rpi relación gas producido/petróleo producido a condiciones iniciales
Para hallar volúmenes de roca a partir de mapas geológicos se utiliza el planímetro. Con
este dispositivo se miden las áreas de cada contorno. Los volúmenes se proyectan, ver Figs.
4.1.a y 4.1.b, de modo que se tengan áreas para leer.
247 acre /1 km
A = lectura planimetro *
lectura inicial
h
∆V = ( A n + An+1 )
2
h
∆V = ( A 0 +2 A1 + 2 A2 "" + 2 An−1 + An ) + hn ⋅ An
2
∆V =
h
3
(
A 1 + A2 + A1 A2 )
Con esta regla piramidal, el volumen total se da mediante1,4:
∆V =
h
3
(
A 1 +2 A2 + 2 A3 + ... + 2 An −1 + 2 An + A1 A2 + A2 A3 + ... + An −1 An )
Cono
A
Hemiesfera
Perspectiva
Sección
B A Transversal
B
C
C
D
D
A
B
A
B C
C
D
Mapa de contorno
Mapa de contorno
Pirámide Elipsoide
A Sección A
B
B Transversal
C
C D
Perspectiva
A D
A B
B
C
D Mapa de contorno
Mapa de contorno E
40
30
Altura
20
10
0 500 1000
Area
60'
40'
20'
0'
60'
40'
20'
0'
Utilizando la regla de integración numérica de Simpson del área bajo la curva extendida a
un volumen for n figuras pares:
h
∆V = ( A1 +4 A2 + 2 A3 + 4 A4 ... + 2 An −2 + 4 An−1 + An )
3
h, ft A, Ac
A0 0 1000
A1 10 720
A2 30 640
A3 30 30
EJEMPLO
El mapa isópaco para la arena petrolífera Gusher se presenta en la Fig. 4.3. La porosidad
promedio es de 23 %, la saturación de agua connata es 32 % y el factor volumétrico inicial
del crudo fue de 1.32 rb/STB. Este yacimiento ha producido a la fecha 3.53 MM STB. Las
áreas dentro de los contornos son: A0 = 640 Ac, A20= 480 Ac, A40=320 Ac y A60 = 160 Ac.
h
HCPV = [ Ao + 2 A20 + 2 A40 + A80 ]
2
20
HCVP = [ 640 + 2 × 480 + 2 × 320 + 160] × 43650 = 104544000 ft 3 = 1.862 x108 bbl
2
Np 3.53 × 106
FR = = = 16 %
N 2.206 × 107
Este es el método más recomendado. Los mapas mostrados en las Figs. 4.4 y 4.5 fueron
generados por el programa Surfer con los datos de la tabla 4.1 usando el método de
interpolación Krigging. Algunos problemas que no se presentan en el trazado manual son
obvios en estos mapas. Primero, la estructura de los contornos tiende a agrandarse o
abultarse hacia el sur, suroeste y noreste, áreas donde no existen puntos. Otros métodos de
interpolación son: curvatura mínima, distancia inversa, método de Shepherd, triangulación,
regresión polinómica y funciones de bases radiales.
10000
8000
6000
4000
2000
-2000
-4000
10000
8000
6000
4000
2000
-2000
-4000
10000
8000
6000
4000
2000
-2000
-4000
10000
8000
6000
4000
2000
-2000
-4000
10000
8000
6000
4000
2000
-2000
-4000
10000
8000
6000
4000
2000
-2000
-4000
Pozo x, ft y, ft D, ft h, ft φ Sw
1 0 -110 -4900 25 0.2 0.5
2 4500 4970 -4920 33 0.205 0.55
3 2640 2640 -4800 44 0.2 0.2
4 660 2100 -4900 36 0.19 0.22
5 2000 770 -4850 22 0.023 0.27
6 2640 4500 -4850 45 0.242 0.23
7 4000 3000 -4880 23 0.235 0.22
8 660 3400 -5100 34 0.23 1
9 5500 1320 -5050 24 0.19 1
10 2200 6500 -5010 11 0.21 1
11 1000 -2000 -5000 24 0.18 1
12 3000 5000 -4990 28 0.12 1
Prom. 31.73 0.2023
Este método tiende a producir curvas más suavizadas que satisfacen los datos en una
manera muy cercana. Sin embargo, lejos de los datos, este método tiende a reversar la
pendiente y producir puntos de cabalgamiento. La Fig. 4.6 fue generada con los mismos
datos de la tabla 4.1. Otro problema adicional, es que se puede observar que el contacto
agua petróleo está a -4900 ft cuando en realidad está a -4960 ft como lo indicará el próximo
ejemplo.
Este algoritmo crea triángulos mediante el trazado de líneas entre puntos. Los triángulos
resultantes se usan para interpolar entre puntos produciendo resultados aceptables en las
zonas donde hay datos. El método trabaja mejor cuando hay cientos o miles de puntos
espaciados uniformemente a lo largo y ancho de la malla, pero no produce ningún contorno
donde no hay datos4. Es muy bueno para reproducir fallas cuando hay suficientes datos.
Ver Fig. 4.8.
EJEMPLO
Un campo petrolero contiene siete y cinco pozos secos cuyos datos se dan en la tabla 4.1.
B = 1.435 rb/STB y OWC = -4960 ft sobre el nivel del mar. Use la información
suministrada en la tabla 4.1 para determinar petróleo in-situ por el método volumétrico y
con la ayuda del Surfer.
SOLUCIÓN
La tabla 4.2 presenta las áreas leídas con el planímetro del mapa isosaturación de la Fig.
4.4. Los contornos se tomaron de 0 a 0.9 con intervalos de 0.1 (h = 0.1). La regla piramidal
se usa para calcular el volumen de crudo:
∆V =
h
3
(A 1 +2 A2 + 2 A3 + ...2 A9 + A10 + A1 A2 + A2 A3 + .. + A9 A10 )
∆V = 3402030 ft2 = 78.1 Ac (incluye la saturación de crudo)
7758 Ah φ 7758(79.1)(31.73)(0.2023)
N= = = 10658020 STB
β oi 1.435
Otra forma más exacta es generando el valor de cada una de las propiedades (espesor,
porosidad y saturación de crudo) para cada una de las diferentes celdas en que puede
dividirse el yacimiento, se debe tener cuidado con anular aquellos valores donde el espesor
o la saturación son cero y el programa puede colocarle algún valor. Surfer permite exportar
los valores de cada celda. Para este ejemplo se tomaron 6500 celdas (100*65) teniendo
cada celda un área de 38600.86 ft2 de acuerdo con los límites del yacimiento de la Fig. 4.4.
Valores de saturación menores de 0.01 (incluyendo valores negativos fueron asumidos
como cero). Usando esta forma el petróleo in-situ resulto ser 19724505 STB.
VOLUMETRICO
INICIAL FINAL
NO-VOLUMETRICO
INICIAL FINAL
Se considera acuífero infinito aquellos cuya extensión relativa al yacimiento guardan una
relación 10:1. Si el acuífero es infinito habrá dos posibilidades a) presión constante, o b)
presión variable (declinable), Qgas > We. Cuando el yacimiento es finito existe empuje
parcial, poco dinámico y la presión declina.
We = 0, Sg y Sw son constantes
Gp = gas producido
Gi = gas inicial
Gr = gas remanente = G - Gp
Luego;
Gi − Gr β gi − β gr
FR = = ; si β gi en pcn / pie3
Gi β gi
1/ β gi − 1/ β gr
FR = ; si β gi en pcn / pie3
1/ β gi
Si P es variable:
G − Gr S gi ⋅ β gi − S gr ⋅ β gr
FR = =
G S gi ⋅ β gi
Esto porque hubo una caída de presión. Si la presión se mantiene constante, entonces Bg se
mantiene constante:
G − Gr S gi − S gr
FR = =
G S gi
EJEMPLO
Calcular el factor de recobro considerando los tres casos si Sgi es 65 %, Sgr es 15 %, βgi =
1215 pcn/pie3 y βgr = 635 pcn/pie3.
SOLUCIÓN
Para volumétrico 47.8 %, para no volumétrico con empuje parcial 87.9 % y con empuje
total 76.9 %.
β oi
FR = 1 −
βo
⎛ 1 − S w − S gr ⎞ β oi
FR = 1 − ⎜ ⎟⎟
⎜ 1 − S w − S gi
⎝ ⎠βo
⎛ 1 − S w − Sor ⎞
FR = 1 − ⎜ ⎟
⎝ 1 − S wi ⎠
⎛ 1 − Sw ⎞ β oi
FR = 1 − ⎜ ⎟
⎝ 1 − S wi ⎠ βo
REFERENCIAS
CAPITULO 5
BALANCE DE MATERIA
El balance de materia expresa que la masa de hidrocarburos producidos es igual a los moles
iniciales menos los remanentes, matemáticamente1;
x p = xi − xr
Vi = Vr
El volumen remanente considera el volumen inicial menos el agua que intruye más el agua
que se produce, es decir:
Vr = Vi − (We − β wW p )
Vi = Vr ; G β gi = (G − G p ) β g , simplificando;
G β gi = G β g − G p β g
G ( β gi − β g ) = −G p β g
β g − β gi G p
= (5.1)
βg G
G G − Gp
=
β gi βg
Rearreglando;
G β g = G β gi − G p β gi
Dividiendo por G;
Gp
β g = β gi − β gi (5.2)
G
pT
β g = cte
z
Dicha constante incluye la temperatura del yacimiento que para un ejemplo dado se
mantiene invariable luego, la Ec. (5.2) se convierte en:
p p p G
cte = cte i − cte i p
z zi zi G
Reorganizando términos:
p p p
= − i Gp + i
z zi G zi
P
Z Caso 1
Caso 2
Qaband G
Gp
Los resultados al graficar dichos datos se ajustarán a una línea recta si el yacimiento es
volumétrico, sin embargo, dos casos pueden ocurrir, como se muestra en la Fig. 5.1.a.
CASO 1:
1. Datos erróneos
2. Condensación retrógrada
3. Drenaje o escape de fluidos del yacimiento
4. Se está drenando más de la cuenta
5. Puede existir un yacimiento con presiones anormales (sobrepresionado)
6. Existe caída de presión a medida que hay reducción de permeabilidad y se pierde
comunicación
Caso 2:
1. Datos erróneos
2. Subsidencia
3. Hay comunicaciones o escapes hacia el yacimiento a lo largo de fallas o problemas
operacionales
pi
m=
zi G
EJEMPLO
Un yacimiento de gas seco puede representarse por un prisma rectangular de 5000 pies de
largo por 3000 pies de ancho y un espesor promedio de 50 pies. La porosidad promedia es
de 12.5 % y la saturación de agua connata es 20 %. La presión inicial fue de 4200 psia (Zi =
0.96) y la temperatura del yacimiento fue de 160 °F.
SOLUCIÓN
V φ (1 − S w ) 93750000 (1 − 0.2 )
G= = = 1.89012 × 1011 scf
Bgi 3.968 ×10−4
T (160 + 460) ft 3
Bgi = 0.0028 Z i = 0.0028(0.96) = 3.968 x10−4
Pi 4200 scf
p P P
= − i Gp + i
Z ZiG Zi
c) Se cree que este yacimiento volumétrico puede producir a una rata constante de 5 MM
scf/día hasta que la presión promedia del yacimiento llega a los 2000 psia. Cuántos años
puede el yacimiento proporcionar esta rata de flujo?
Luego:
Gp 7.53 × 1010
t= = = 41.26 anos
q 5MM × 365
EJEMPLO
SOLUCIÓN
Dividiendo el peso molecular del gas por el del aire resulta una gravedad específica de
0.6845. Calculando las presiones y temperaturas pseudoreducidas es posible obtener el
valor de Z, a saber:
5000
4500
4000
3500
3000
y = -6.3711x + 4558.5
2
P/Z, psi R = 0.9927
2500
2000
1500
1000
500
0
0 50 100 150 200 250
P Z P/Z Gp
4000 0.8882 4503.648 0
3825 0.8758 4367.226 34
3545 0.8579 4131.952 75
3388 0.8491 3990.198 91.3
3192 0.8394 3802.714 120
2997 0.8315 3604.222 145
2738 0.8242 3322.006 192
Al graficar P/Z se obtiene cuando P/Z=0 que G = 4558.5 MMM scf de gas. Ver Fig. 5.1.b.
n f = ni − n p
Aplicando la ley de los gases reales y asumiendo que los moles producidos son referidos a
condiciones iniciales:
PV f PV PG
= i i
− i p
ZRT Z i RT Z i RT
De donde:
PV f PV PG
= i i
− i p
Z Zi Zi
P(Vi − We ) PV PG
= i i− i p
Z Zi Zi
P(Vi − We ) Pi
= (G − Gp )
Z Zi
P(Vi − We ) Pi
= (1 − FR )
Z Zi
De donde:
P Pi (1 − FR )
=
Z Z i (Vi − We )
P Pi (1 − FR )
=
Z Z i (1 − We / Vi )
EJEMPLO
P Pi ⎛ G p ⎞ ⎛ We ⎞
= ⎜1 − ⎟⎟ / ⎜ 1 − ⎟
Z Z i ⎜⎝ G ⎠ ⎝ Vi ⎠
Una suposición razonable para un yacimiento con empuje de agua de fondo es que la
diferencia de presión puede transmitirse a través del acuífero. La cantidad de agua intruida
puede determinarse como la cantidad de agua expandida. Si la compresibilidad total del
acuífero es de 5x10-6 psi-1 y el volumen de acuífero es 20 veces el volumen poroso
disponible a hidrocarburos de este yacimiento gasífero, cual es el factor de recobro
adicional a 2500 psia comparado con un yacimiento volumétrico?
SOLUCIÓN
P Pi W
= (1 − FR ) /1 − e
Z Zi Ui
Vacui = 20 ( V φ )
De la definición de compresibilidad:
1 dV 1 dWe
C=− luego C = −
V dP Vi ∆P
P Pi ⎛ G p ⎞ ⎛ We ⎞
= ⎜1 − ⎟ / ⎜1 − ⎟
Z Zi ⎝ G ⎠ ⎝ Vi ⎠
0.561797 =
[1 − FR ]
0.75
P Pi
= (1 − FR )
Z Zi
4000 5000 ⎛ 42 ⎞
= ⎜ 1 − MMM ⎟
0.96 1 ⎝ G ⎠
G p 2500 = G × FR
P Pi ⎛ G p ⎞
= ⎜1 − ⎟
Z Zi ⎝ G ⎠
Gp⎛ P Zi ⎞
− =⎜ ⎟ −1
G ⎝ Z Pi ⎠
⎛ 2500 1 ⎞
FR = 1 − ⎜ ⎟ − 43.82 %
⎝ 0.89 5000 ⎠
G p β gi = G ⎡⎣ β gi − β g ⎤⎦ + We − W p β w
F: extracciones
E: expansiones
Entonces;
F = G p β gi + W p β w
E = β gi − β g , luego:
F = GEg + We
De donde;
F W
=G+ e ⇔ y = mx + b
Eg Eg
F
Eg
WeD
Eg
Para solucionar y conocer las incógnitas se utiliza el método de Hurst & Van Everdingen
quienes resolvieron la ecuación de difusividad para dos casos:
Con el método de Hurst y van Everdingen1,2,4-7, se tiene una intrusión de agua adimensional
dada por;
WeD = β ∑ ∆p j ⋅ QtD
F W
= G + β eD
Eg Eg
Una vez conocido We se construye una gráfica similar a la de la Fig. 5.1.c, para estimar el
valor de G.
EJEMPLO:
r = 1500 ft, φ = 25 %, ct = 5x10-5 /psia, Pi = 3274 psia, viscosidad = 30 cp y permeabilidad
0.01 md. Los resultados de entrada y salida son:
θ
β = 1.119φ ce ryto
2
h
360
Siendo θ el ángulo subtendido por la circunferencia del yacimiento semicircular, los demás
parámetros se dan en unidades de campo.
co So + cw S w + cg S g + c f
ce =
1 − Sw
Puesto que se trata del acuífero, la ecuación para la compresibilidad efectiva queda:
ce = cw + c f
La Fig. 5.1.d muestra un tanque (yacimiento) con capa de gas a condiciones iniciales. Allí,
m es el tamaño de la capa de gas, ésta se obtiene a partir de registros, datos de corazones,
datos de completamiento, presiones de fondo y mapas estructurales. También se ha
presentado una disminución en el volumen poroso disponible a hidrocarburos causados por
expansión de roca y fluido.
GAS
mNBoi
mNBoi(Bg/Bgi)
Np, Gp, Wp
OIL
NBoi
(N-Np) Bo
Gas libre
(We-Wp)Bw
WOC
Np = Petróleo producido, BF
N = Petróleo original in-situ, BF
G = Gas inicial en el yacimiento
m = Tamaño inicial de la capa de gas o volumen inicial de la capa de gas/volumen de la
zona de petróleo (N)
Np = Petróleo producido acumulado, BF
Gp = Gas producido acumulado, pcn
Wp = Agua producido acumulado, BF
Rp = Relación gas-petróleo acumulada, Gp/Np, pcn /BF
Rs = Relación gas-petróleo, pcn/BF
ßo, ßw = Factor volumétrico de formación del petróleo y del agua, bbl/BF
ßg = Factor volumétrico de formación del gas, bbl/pcn
We = Intrusión acumulada de agua, BF
Sw = Saturación de agua, fracción
cw, co, cg = Compresibilidad del agua, del petróleo y de gas, 1/psi
cf = Compresibilidad del volumen poroso, 1/psi
P = Presión estática del yacimiento, psia
∆P = Pi - P
i = inicial
De modo que:
βg
Capa de gas original expandida = mN β oi , bbl
β gi
Petróleo remanente = (N - Np) ßo, bbl
Gas neto liberado = Gas liberado - gas libre producido
gas liberado = N (Rsi - Rs)ßg, bbl
gas libre producido = Np (Rp - Rs) ßg, bbl
Gas en solución producido = Np Rs, bbl
Gas neto liberado = N (Rsi - Rs)ßg - Np (Rp - Rs) ßg, bbl
Volumen hidrocarburos remanentes =
(N - Np)ßo + mNßoi (ßg/ßgi) + N (Rsi - Rs) ßg – Np (Rp - Rs) ßg, bbl
Intrusión neta de agua = (We - Wp) ßw
1 dV
ce = −
V dP
⎛ cw S w + c f ⎞
− d ( HCVP) = (1 + m) N β oi ⎜ ⎟ ∆P
⎝ 1 − Sw ⎠
⎛ βg ⎞
(1 + m) N β oi = ( N − Np ) β o + mN β oi ⎜ + N ( Rsi − Rs ) β g − Np ( R p − Rs )
⎜ β ⎟⎟
⎝ gi ⎠
⎛ cw S w + c f ⎞
+ (W − W p ) β w + (1 + m) N β oi ⎜ ⎟ ∆P
⎝ 1 − S w ⎠
Agrupando los términos con N al lado izquierdo y con Np al lado derecho, se tiene:
⎡ ⎛β ⎞ ⎛ c S + cf ⎞ ⎤
N ⎢ β o − β oi + mβ oi ⎜ g − 1⎟ + β g ( Rsi − Rs ) + β oi (1 + m) ⎜ w w ⎟ ∆P ⎥ =
⎜β ⎟
⎢⎣ ⎝ gi ⎠ ⎝ 1 − Sw ⎠ ⎥⎦
N p ⎡⎣ β o + β g ( R p − Rs ) ⎤⎦ − (We − W p ) β w
Despejando N;
N p ⎡⎣ β o + β g ( R p − Rs ) ⎤⎦ − (We − W p ) β w
N=
⎛ βg ⎞ ⎛ c S + cf ⎞
β o − β oi + mβ oi ⎜⎜ − 1⎟ + β g ( Rsi − Rs ) + β oi (1 + m) ⎜ w w ⎟ ∆P
⎟
⎝ β gi ⎠ ⎝ 1 − Sw ⎠
DISCUSIÓN
1) Desarrollar una relación entre los fluidos producidos con la presión para predecir el
comportamiento del yacimiento
2) Verificar la existencia de la capa de gas
3) Comparar con el petróleo calculado volumétricamente. Esto no significa que los valores
sean cercanos. Balance de materia siente la presencia de una falla (2 yacimientos) o
puede haber una extensión en el yacimiento que el método volumétrico no detecta.
4) No es bueno tratar de hallar N con la Ecuación de balance de materia en yacimientos
con alta intrusión de agua porque Pi - P es pequeño.
Ejercicio:
Calcule el petróleo in-situ para un yacimiento de petróleo cuya presión inicial es de 5000
psia y su presión de burbuja es de 2750 psia6. Otros datos son:
compresibilidad del petróleo. Para el segundo caso parta la EBM y no. Desprecie los
efectos de cw y cf.
Solución6: Primero que todo es necesario calcular cf, ya que el dato dado fue cr, la
compresibilidad de la matriz.
φ = 1 − Vr (5.3)
1 dVr
cr = − (5.4)
Vr dP
1 dφ
cf = (5.5)
φ dP
Derivando la Ec. (5.3) con respecto a P y sustituyendo este resultado en la Ec. (5.5)
1 ⎛ dVr ⎞
cf = − (5.6.a)
φ ⎜⎝ dP ⎟⎠
Despejando dVr/dP de la Ec. (5.4) y reemplazando ese resultado en la Ec. (5.6.a) se tiene:
1
cf = ( cr ⋅Vr ) . De la Ec. (5.3) Vr = 1 - φ, luego:
φ
1
cf = cr (1 − φ ) = 2.7 × 10−6 / psia
φ
N p βo
N= (5.6.b)
⎛c S +c ⎞
β o − β oi + β oi ⎜ w w f ⎟ ∆P
⎝ 1 − Sw ⎠
N p βo 1510000(1.338)
N= = = 61.2 MMBF
β o − β oi 1.338 − 1.305
b) Considerando cw y cf
N p βo
N=
⎛c S +c ⎞
β o − β oi + β oi ⎜ w w f ⎟ ∆P
⎝ 1 − Sw ⎠
1510000(1.338)
N=
⎛ 3.5 ×10−6 (0.216) + 2.7 × 10−6 ⎞
1.338 − 1.305 + 1.305 ⎜ ⎟1650
⎝ 1 − 0.216 ⎠
1.51MM (1.338)
N= = 47.5 MMBF
0.033 + 1.305 ( 4.408 ×10−6 )1650
∆ = 29 %
⎣ ⎦
Los términos con exponentes mayores que uno tienden a cero por lo tanto:
β o = β oi ⎡⎣1 + co ( Pi − P ) ⎤⎦ , luego
β o = β oi + β oi co ( Pi − P ) , de donde:
N p βo
N=
β oi co ∆P + β oi ce ∆P
N p βo 1.51MM (1.338)
N= = = 47.6 MMBF
β oi ce ' ∆P 1.305 (15.3 + 4.408 ) × 10−6 (1650)
TZ (460 + 20)0.9
β g = 0.00503 = 0.0050) = 0.00211 bbl /
p 1500
N p βo
N=
β oi coe ∆P
De donde;
Eo
F = N p ⎡⎣ β o + ( R p − Rs ) β g ⎤⎦ + W p β w
Eo = ( β o − β oi ) + ( Rsi − Rs ) β g
⎛β ⎞
Eg = β oi ⎜ g − 1⎟
⎜ β gi ⎟
⎝ ⎠
⎡ cw S w + c f ⎤
E f , w = β oi (1 + m ) ⎢ ⎥ ∆P
⎣ 1 − Sw ⎦
F = N ⎡⎣ Eo + mEg + E f , w ⎤⎦ + We β w
1) CASO: m = 0, We = 0, Cw,f = 0
F = N Eo
Al graficar F vs. Eo debería dar una recta que pasa por el origen. Fig. 5.3. Este es el
tratamiento más indicado para arrancar un análisis. Si el gráfico no es lineal se está
definiendo la existencia de otro mecanismo de producción.
m asumido
muy pequeño
F
N
m asumido
muy grande
Eo + mEg
F
Eo mN
Eg
Eo
Pi
Relación gas-petróleo
Presión
producida (Rp)
Parámetro
Corte de agua
Rsi
Tiempo
Fig. 5.6. Comportamiento de la presión, gas en solución y corte de agua (caso 3)1
F = N ( Eo + mEg )
El caso normal es que m no se conozca. Se suponen varios m’s. El m correcto da una línea
recta. Si hay incertidumbre tanto en el valor de m como en el de N se sugiere rescribir la
ecuación de Havlena y Odeh de la siguiente manera:
F = NEo + mNEg
de donde:
F E
= N + mN g
Eo Eo
F = N ( Eo + We Bw )
F WB
=N+ e w
Eo Eo
We muy pequeño
We correcto
F We muy grande
Eo + mEg
45°
We Bw
Eo + mEg
Tiempo
Fig. 5.8. Comportamiento de la presión, gas en solución y corte de agua (caso 4)1
Este caso se presenta cuando el acuífero no es muy activo, pero se asume que tanto m como
N se conocen. De modo que:
F = N ⎡⎣ Eo + mEg ⎤⎦ + We β w
F We β w
=N+
Eo + mEg Eo + mEg
EJEMPLO
SOLUCIÓN
⎛β ⎞ ⎛ 0.000857 ⎞
Eg = β oi ⎜ g − 1⎟ = 1.5 ⎜ − 1⎟ = 0.127215
⎜ β gi ⎟ ⎝ 0.00079 ⎠
⎝ ⎠
3.5E+08
3.0E+08
2.5E+08
2.0E+08
F/Eo
1.5E+08
1.0E+08
5.0E+07
0.0E+00
0 0.5 1 1.5 2 2.5
Eg/Eo
1) Al graficar estos datos, ver Fig. 5.9, se tiene que la pendiente (mN) es 9x107 y el
intercepto es 1x108, el cual corresponde a N. Luego el tamaño de la capa de gas es
9x107/1x108= 0.9. Luego la cantidad de gas inicial está formada por el gas en solución
(RsixN) = 1x1011 scf más el gas en la capa de gas (mN*5.615) que corresponde a
505.35x106 scf, para un total de 1.00505x1011 scf.
3) El gas remanente en solución (SCF) a 1500 psi se calcula estimando la cantidad de crudo
remanente a 1500 psi, lo cual es 1x108 – 39.5 MM = 6.05x107 STB y multiplicado por el
gas en solución (380 scf/STB) da 2.3x1010 scf.
4) El gas libre en el yacimiento (SCF) a 1500 psi se estima con los 6.05x107 STB
remanentes multiplicado por (Rsi-Rs) = (1000-380) = 3.751x1010 scf.
⎛ 39.5MM ⎞
(1 − Sw ) ⎛⎜1 − ⎞
NP
⎟ βo (1 − 0.23) ⎜1 − ⎟1.272
SL = Sw + ⎝ N ⎠ = 0.23 + ⎝ 100 MM ⎠ = 0.625
β ob 1.5
Otra forma es estimando la cantidad de gas libre en bbl es 3.751x1010 scf x 0.00216 bbl/scf =
81021600 bbl. El volumen total del yacimiento es 150399326.5 bbl/(1-0.23) = 195323800.6
bbl. Luego 81021600/195323800.6 es igual a 41.5 %. Como el yacimiento es volumétrico, la
cantidad de crudo es 0.625-0.23 = 0.395, el remanente es gas, o sea 0.375. La diferencia es
que la ecuación es válida para sistemas sin capa de gas.
Es el yacimiento más simple que se considera para efectos de cálculo. Se considera que la
permeabilidad vertical es limitada, no hay capa de gas, la intrusión de agua es despreciable.
Por encima del punto de burbuja hay una sola fase (aunque puede existir una fase sólida de
asfaltenos precipitados). La ecuación de balance de materia para éste tipo de yacimiento,
teniendo en cuenta las consideraciones mencionadas anteriormente es1:
N p βo
N=
⎡c S + c ⎤
β oi co ∆P + β oi ⎢ w w f ⎥ ∆P
⎣ 1 − Sw ⎦
cw S w + c f
coe = co +
1 − Sw
Luego;
N p βo
N=
β oi coe ∆P
Las estimaciones del petróleo original basado en la historia de producción sobre el punto de
burbuja a menudo no son confiables por la sensibilidad al gradiente de presión.
EJEMPLO
SOLUCIÓN
N p ⎡⎣ β t + β g ( R p − Rsi ) ⎤⎦
N= , despejando Rp:
β
βt − βti + m ti ( β g − β gi )
β gi
N ⎡ β ti ⎤
⎢ β t − β ti + m ( β g − β gi ) ⎥ − β t
Np ⎢⎣ β gi ⎥⎦
Rp = + Rsi
βg
56 ⎡ 1.498 ⎤
⎢1.523 − 1.498 + (0.5) (0.001155 − 0.001048) ⎥ − 1.523
3.741 ⎣ 0.001048 ⎦
Rp = + 721
0.001155
Como puede verse en la tabla 5.4 y en la Fig. 5.10.a a partir de un a presión de 1500 psi se
observa diferencias entre los valores de Rp indicando un aumento en la producción de gas,
por lo que se concluye que en este punto inicia la inyección de gas.
2) Para estimar la cantidad de gas inyectada se obtiene la diferencia de los Rp’s a 500 psi, lo
cual corresponde a 1802 scf por cada barril de crudo producido, luego multiplicando por Np
(19.59 MM STB) se tiene una inyección total de 35.3 MMM scf de gas.
6000
Rp Medido
Rp Calc.
5000
4000
Rp, scf/STB
3000
2000
1000
0
0 500 1000 1500 2000 2500
P, psi
A medida que la presión cae por debajo del punto de burbuja se libera gas. Al principio son
burbujas de gas muy pequeñas y no continuas. A medida que cae la presión las burbujas se
revientan y se forma una capa de gas continua y el gas empieza a fluir.
N p ⎡⎣ β o + β g ( R p − Rsb ) ⎤⎦
Nb =
⎛ cw S w + c f ⎞
β o − β ob + β g ( Rsb − Rs ) + β oi ⎜ ⎟ ∆P
⎝ 1 − Sw ⎠
N p ⎡⎣ β o + β g ( R p − Rsb ) ⎤⎦
Nb =
β o − β ob + β g ( Rsb − Rs )
CONSIDERACIONES:
1) Yacimiento volumétrico
2) Proceso isotérmico
3) Existe equilibrio de presiones. Las propiedades del fluido son las mismas en todas
partes.
4) Los datos PVT son confiables
5) m = 0
6) El recobro es independiente de la rata de flujo
7) Datos de producción confiables
8) Las condiciones iniciales son el punto de burbuja
9) We = 0
Aunque algunas reglas son violadas, la anterior ecuación es una herramienta poderosa ya
que a ella es que se aplican los diferentes métodos de ensayo y error para predecir el
comportamiento y recobro final de un yacimiento de petróleo.
βt = β o + ( Rsi − Rs ) β g
N p ⎡⎣ β t + β g ( R p − Rsi ) ⎤⎦ − (We − W p ) β w
N=
⎛ cw S w + c f ⎞ βti
βt − βti + ⎜ ⎟ ∆P β ti + m ( β − β gi )
⎝ 1 − Sw ⎠ β gi g
β ti ⎛ cw S w + c f ⎞
N p ⎡⎣ βt + β g ( R p − Rsi ) ⎤⎦ = N ( βt − β ti ) + Nm ( β g − β gi ) + (We − W p ) + ⎜ ⎟ N ∆P β ti
β gi ⎝ 1 − Sw ⎠
βti
Nm ( β − β gi )
N ( β t − βti ) β gi g (We − W p )
1= + + +
N p ⎡⎣ β t + β g ( R p − Rsi ) ⎤⎦ Np ⎡⎣ βt + β g ( R p − Rsi ) ⎤⎦ N p ⎡⎣ β t + β g ( R p − Rsi ) ⎤⎦
⎛ cw S w + c f ⎞
⎜ ⎟ N ∆P β ti
⎝ 1 − Sw ⎠
N p ⎡⎣ βt + β g ( R p − Rsi ) ⎤⎦
Cinco métodos son estudiados para tal fin, estos se presentan con detalle en los anexos del
presente texto, a saber:
Método de Schilthuis9,10
Método de Tarner9,10
Método de Tracy9,10
Método de Muskat10
Método de Pirson4,10
⎪⎧ ⎡ ⎛ Bg ⎞ ⎛ c S + c m ⎞ ⎤ ⎪⎫
N p ⎣⎡ Bo + ( R p − R s ) B g ⎦⎤ = N 1 ⎨ ( Bt − Bti ) + Boi ⎢ m ⎜ − 1 ⎟ + (1 + m )⎜ w wmi ⎟ ∆p ⎥ ⎬ +
⎩⎪ ⎢⎣ ⎜⎝ B gi ⎟
⎠ ⎝ 1 − S wmi ⎠ ⎥⎦ ⎪⎭
⎧⎪ ⎡ ⎛B ⎞ ⎛ c S + cf ⎞ ⎤ ⎫⎪
N 2 ⎨ ( Bt − Bti ) + Boi ⎢ m ⎜ g − 1 ⎟ + (1 + m )⎜ w wfi ⎟⎟ ∆ p ⎥ ⎬ (5.7)
⎪⎩ ⎢⎣ ⎜⎝ B gi ⎟
⎠
⎜ 1 − S wfi
⎝ ⎠ ⎥⎦ ⎪⎭
en donde el primer termino representa el volumen total (gas y petróleo) que ocupa el
espacio poroso en todo el sistema; el segundo término se refiere a el volumen de
transferencia (crudo y gas) de matriz a fractura y el tercer término expresa la transferencia
de fluido (crudo y gas) de las fracturas a superficie. Además:
Bt = Bo + (Rsi − Rs ) B g
G B gi
m=
N Boi
CASO 0
Este considera el caso convencional de Havlena y Odeh (sin capa de gas) que consiste en
reorganizar la ecuación para obtener grupos de variables de manera que se comporten como
una línea recta al ser graficados, donde la EBM general es:
⎪⎧ ⎛c S +c ⎞ ⎪⎫
N p ⎡⎣ Bo + ( R p − Rs ) Bg ⎤⎦ = N1 ⎨ Bo − Boi + ( Rsi − Rs ) Bg + ⎜ w wi m ⎟ ∆P Boi ⎬ +
⎩⎪ ⎝ 1 − S wi ⎠ ⎭⎪
⎧⎪ ⎛ cw S wfi + c f ⎞ ⎫⎪
N 2 ⎨ Bo − Boi + ( Rsi − Rs ) Bg + ⎜ ⎟⎟ ∆P Boi ⎬
⎜ 1 − S wfi
⎩⎪ ⎝ ⎠ ⎭⎪
[
F = N p Bo + (R p − Rs ) Bg ]
⎛c S +c ⎞
Eo1 = Bo − Boi + ( Rsi − Rs ) Bg + ⎜ w wi m ⎟ ∆P Boi
⎝ 1 − S wi ⎠
⎛ cw S wfi + c f ⎞
Eo 2 = Bo − Boi + ( Rsi − Rs ) Bg + ⎜ ⎟⎟ ∆P Boi
⎜ 1 − S wfi
⎝ ⎠
F = N1 Eo1 + N 2 Eo 2
donde Eo1 representa la expansión neta del petróleo original en el sistema matriz y Eo2 es la
expansión neta del petróleo original en el sistema de fracturas. La anterior expresión puede
escribirse como:
F E
= N1 + N 2 o 2
Eo1 Eo1
De modo que una gráfica de F/Eo1 vs. Eo2/Eo1 da una recta de pendiente N2 y corte N1 como
ilustra la Fig. 5.10.b.
CASO 1
Conociendo el tamaño de la capa de gas, desconociendo el volumen de petróleo original en
cada uno de los medios porosos, dividiendo toda la ecuación (5.7) entre el factor
volumétrico del petróleo a la presión inicial y utilizando las siguientes definiciones:
Np
F=
Boi
[B + (R
o p − Rs ) B g ]
N2
F/E o1
N1
E o2/E o1
3
Fig. 5.10.b. Gráfico de F/Eo1 vs. Eo2/Eo1
N1
F/x 2
N2
x1/x2
⎛ c w S wfi + c f ⎞
cef =⎜ ⎟
⎜ 1− S ⎟
⎝ wfi ⎠
⎛ c S + cm ⎞
cem =⎜⎜ w wmi ⎟⎟
⎝ 1 − S wmi ⎠
⎪⎧ ⎛ B ⎞ ⎛ B ⎞ ⎪⎫ ⎪⎧ ⎛ B ⎞ ⎛B ⎞ ⎪⎫
F = N 1 ⎨ ⎜ t − 1 ⎟ + m ⎜ g − 1 ⎟ + (1 + m ) c em ∆ P ⎬ + N 2 ⎨ ⎜ t − 1 ⎟ + m ⎜ g − 1 ⎟ + (1 + m ) cef ∆ P ⎬
⎜ ⎟ ⎜ ⎟
⎩⎪ ⎝ Bti ⎠ ⎝ B gi ⎠ ⎭⎪ ⎪⎩ ⎝ Bti ⎠ ⎝ B gi ⎠ ⎪⎭
(5.8)
Como se conoce m, se tiene una expresión lineal de la forma y = a1 x1 + a2 x2 con:
⎛ B ⎞ ⎛ Bg ⎞
x1 = ⎜ t − 1⎟ + m ⎜ − 1⎟ + (1 + m ) cem ∆P
⎜ Bgi ⎟ ⎜ Bgi ⎟
⎝ ⎠ ⎝ ⎠
⎛ B ⎞ ⎛ Bg ⎞
x2 = ⎜ t − 1⎟ + m ⎜ − 1⎟ + (1 + m ) cef ∆P
⎜ Bgi ⎟ ⎜ Bgi ⎟
⎝ ⎠ ⎝ ⎠
que puede solucionarse por regresión lineal para N1 y N2. Una alternativa de este caso es
realizar un gráfico que describa una línea recta en donde y/x2 es la variable dependiente de
x1/x2, cuya pendiente es el petróleo original en la matriz y cuya ordenada es el petróleo
original en las fracturas. Ver Fig. 5.10.c. Al tomar la ecuación lineal y hacer el cociente con
la variable x2:
y x
= a1 1 + a2
x2 x2
CASO 2
N1cem + N 2 cef
cem− f = (5.9)
N1 + N 2
N2
γ=
N1 + N 2
⎛B ⎞ ⎛ Bg ⎞
F = N ⎜⎜ t − 1⎟⎟ + mN ⎜ − 1⎟ + N (1 + m ) cem− f ∆p
⎜B ⎟
⎝ Bti ⎠ ⎝ gi ⎠
a1 = N
a2 = mN
a3 = N (1 + m )cem− f
⎛B ⎞
x1 = ⎜⎜ t − 1 ⎟⎟
⎝ B ti ⎠
⎛ Bg ⎞
x2 = ⎜ − 1⎟
⎜B ⎟
⎝ gi ⎠
x3 = ∆P
Una vez conocidos y ,x1 , x2 y x3, se hallan los coeficientes a1, a2 y a3, a partir de los cuales
se calculan las incógnitas de la siguiente forma:
N = a1
a2
m=
a1
a3
c em − f =
N (1 + m )
Para la solución por regresión no lineal, la ecuación (5.8) puede ser escrita de manera que
aparezcan todas las incógnitas posibles, de la siguiente forma:
2.40E+07
2.35E+07
2.30E+07
2.20E+07
2.15E+07
2.10E+07
2.00E+07
0.64 0.65 0.66 0.67 0.68 0.69 0.7 0.71 0.72 0.73 0.74
x1/x2
⎧⎪⎛ B ⎞ ⎛ Bg ⎞ ⎫⎪ ⎧⎪⎛ B ⎞ ⎛ Bg ⎞ ⎫⎪
F = N (1 − γ ) ⎨⎜ t − 1⎟ + m ⎜ − 1⎟ + (1 + m ) cem ∆p ⎬ + N γ ⎨⎜⎜
t
− 1⎟ + m ⎜ − 1⎟ + (1 + m ) cef ∆p ⎬
⎜ ⎟ ⎜ Bgi ⎟ ⎟ ⎜ Bgi ⎟
⎩⎪⎝ Bgi ⎠ ⎝ ⎠ ⎭⎪ ⎩⎪⎝ Bgi ⎠ ⎝ ⎠ ⎭⎪
Esta ecuación puede ser solucionada mediante la utilización de regresión no lineal entre
cualquier combinación de las variables N, γ, m y cef.
EJEMPLO
SOLUCION
La tabla 5.5 presenta los resultados de los cálculos de F, x1 y x2 con los cuales se construyó
la Fig. 5.10.d. Obsérvese que inicialmente los datos no se ajustan a una línea recta
probablemente porque las suposiciones del modelo no se ajustan con la realidad. Se obtuvo
que la cantidad de crudo en la matriz es de 30’421.540 STB y en el sistemas fracturado es
de 1’416.798 STB.
Ver Fig. 5.11. Para un sistema de separación de dos etapas, la gravedad específica, γw, del
fluido total del pozo es:
141.5
γo =
API + 131.5
Et Et
B. Sin capa de gas, m=0. Nota: El gráfico de P/Z vs. F usualmente se prefiere
en vez de este
N
e=
ent
F i
nd
Pe
Et
= mN
1
i ente te
=
ien
F/Eg
nd
F/E o Pe Pe
nd
N G
E g / Eo W e / Eg
Et Et Et Et
We correcto We correcto
We muy Pendiente=1 We muy Pendiente=1
pequeño pequeño
F/Et
F/Et
Geometria Geometria
incorrecta incorrecta
N We muy G We muy
grande grande
W e / Et o Σ∆Pj Q tD / E t W e / Et o Σ∆Pj Q tD / E t
4. Método de Campbell 4. Método de Cole
Empuje de agua. Grafique F/Et vs. F. El intercepto Empuje de agua. Grafique F/Et vs. F. El intercepto
es N es G
n
ió n rus
ió
us int
i n tr n
n
Co agua Co agua
F/Et
F/Et
d e de
Sin intrusión
Sin intrusión
de agua
de agua
G
N
F
F
P
Np Gp
1. Método P/Z
Cualquier yacimiento. Grafique P/Z vs. Gp.
El intercepto con el eje X da G. Se obtiene línea
recta si We=Ef,w=0.
P ⎛ Gp ⎞ Pi
= ⎜1 − ⎟
Z ⎝ G ⎠ Zi
EMPUJE
DE
DE AGUA
PLE
P/Z
CI
O
N
PRESIONES
ANORMALES
G
Gp
Cada método se aplica a un yacimiento específico. Los siguientes métodos de línea recta deberían usarse en
los siguientes casos:
1. Yacimientos de gas en solución. Incógnita N. Use el método 1.
2. Yacimiento con capa de gas. Incógnitas N y m. Use el método 2.
3. Yacimientos con empuje de agua. Incógnitas We y N. Use método 3, 4.
Todas las tres incógnitas pueden calibrarse durante ajuste histórico. Para ésto se requiere la siguiente
información:
⎛ co So + cw S w + c f ⎞
N p β o = N β oi ⎜ ⎟ ∆P
⎝ 1 − Sw ⎠
N p β o = N β oi coe ∆P
Subsaturado F = NEo N p ⎣⎡ β o + β g ( R p − Rs ) ⎦⎤ = N ⎣⎡ β o − β oi + β g ( Rsi − Rs ) ⎦⎤ So<Soi F-vs-Eo, pend=N, b=0. Si existe desviación hacia
volumétrico, Sw=Swi arriba indica la presencia de otro mecanismo
Pi < Pb Sg=Sgi
N p ⎡⎣ β o + β g ( R p − Rs ) ⎤⎦ = N [ β t − β ti ]
βo
S g = 1 − (1 − FR )(1 − S wc )
β oi
Tabla 5.5.b. Resumen factores de recobro por método volumétrico y balance de materia
GP (Superficie)
Gsep. primario Gtanque almac.
Fluido Tanque de
del pozo Separador almacena-
primario miento
GP (Superficie)
Gsep. primario Gsep. secundario Gtanque almac.
Fluido Tanque de
del pozo Separador Separador
almacena-
primario secundario
miento
5954 42.43γ o
M wo = =
API − 8.811 1.008 − γ o
donde, La cantidad inicial de gas in-situ total por ac-ft de yacimiento es:
pV
G ' = 35.35
zT
ng R1 + R3
fg = =
ng + no 132790γ o
R1 + R3 +
M wo
G = fgG '
G
N=
R1 + R3
Gas diario
∆G p =
fg
⎛ ZT ⎞
∆V = ∆G p ⎜ ⎟
⎝ 35.35 p ⎠
EJEMPLO
Calcula la cantidad de gas y petróleo inicial por ac-ft para un yacimiento de gas
condensado1. Dados:
Petróleo diario en tanque, STB = 242 Pi, psia = 2740 φ, % = 25
TR, °F = 215 Swc, % = 30 API = 48°
Gas diario en separador, Mcf = 3100 γg en separador = 0.65
Gas diario en tanque, Mcf = 120 γg en tanque = 1.2
SOLUCIÓN
141.5
γo = = 0.788
48 + 131.5
5954
M wo = = 151.9
48 − 8.811
Con la gravedad del gas se hallan las condiciones críticas, Tcr = 423 °R y Pcr = 637 psia.
Luego z = 0.825.
pV (2740)(43560*0.25*(1 − 0.3))
G ' = 35.35 = 35.35 = 1326 Mft 3 / ac − ft
zT (0.825)(675)
12810 + 496
fg = = 0.951
132790(0.788)
12810 + 496 +
151.9
1261000
N= = 94.8 STB / ac − ft
12810 + 496
⎛ 0.825(675) ⎞
∆V = 3386000 ⎜ ⎟ = 19450 ft / dia
3
⎝ 35.35(2740) ⎠
G ( Bg − Bgi ) + We = G p Bg + BwW p
Esta ecuación se usa para halla ya sea G o We. Para presiones por debajo del punto de rocío,
el factor volumétrico debe ser bifásico, el cual está contenido en el término Z.
REFERENCIAS
CAPITULO 6
INTRUSIÓN DE AGUA
Existen tres métodos básicos para predecir y/o estimar la cantidad de agua que
intruye al yacimiento. Ellos son: El método de Schilthuis, el método de Hurst y Van
Everdingen y el método de Fetckovick.
Es el método más simple que se utiliza para fines de estimación de la intrusión de agua. Se
prefiere usar al principio. Se asume acuífero gigante y altamente permeable (Mayor de 50
md) para que la presión nunca caiga y esquematizado en la Fig. 6.1. Este método integra la
ley de Darcy, la cual es estado estable. Schilthuis arranca con la ley de Darcy:
Qw = Cs ( Pi − p )
( e )n s ∑ ⎡⎢ Pi − ( Pj −1 + Pj )⎤⎥∆t j
n
1
W = C
j =1 ⎣ 2 ⎦
Siendo;
Qw = flujo de agua a través del WOC
Cs = Constante del acuífero. Contiene las constantes de Darcy: k, µ, A y sus unidades son
bbl/tiempo/psia.
Pi = Presión del acuífero
P = Presión estática del yacimiento.
Pj = presión estática a tj
Pj-1 = Presión estática entre tj-1 y tj.
Altura de agua
Pi
constante
Pi
P1
P2
P3
t1 t2 t3
t
Fig. 6.2. Método de Schilthuis
N p ⎡⎣ β o + β g ( R p − Rs ) ⎤⎦ + W p β w − We
N=
β g ( Rsi − Rs ) − ( β oi − β o )
N p ⎡⎣ β o + β g ( R p − Rs ) ⎤⎦ + W p β w
N=
D
Si no hay intrusión de agua el Na permanece constante para todos los periodos. Pero si
existe intrusión de agua el Na se incrementa en cada periodo. Ver Fig. 6.3. El método de
Schilthuis se usa con la ecuación de balance de materia para:
Si Na = N + We/D, entonces:
We ≠ 0
Na
We = 0
Np
Na
m = Cs
N = OOIP
∆P∆t
∑ D
⎛ ∆p∆t ⎞
N a = N + Cs ⎜ ∑ ⎟
⎝ D ⎠
D ' = β g − β gi
G p β g + Wp β w
Ga =
β g − β gi
Si el yacimiento es de gas o petróleo hay dos formas para obtener la “mejor” línea recta:
2) Mínimos cuadrados
Pi
P
variable
We = 2 π α h φ ce ryto
2
[ ∆P QtD ]
α = 0 - 1. Representa la extensión a la cual el acuífero rodea el yacimiento
We = Cantidad de agua intruida en cm3
h = Espesor del estrato en cm
ce = Compresibilidad efectiva del acuífero, 1/atm.
ryto = Radio del yacimiento de petróleo o gas, cm
∆P = Caída de presión constante a través del acuífero, atm
QtD = Función acumulativa de entrada de agua
tD = Tiempo adimensional
P0 P1
P2
P3
P4
P5
Presión
t0 t1 t2 t3 t4 t5
Tiempo
0.00633kt
tD = ; t en dias
φ µ wce ryto
2
0.0002637 kt
tD = ; t en horas
φ µ wce ryto
2
We = 1.119 α h φ ce ryto
2
[ ∆p QtD ]
Este método se basa en la solución de la ecuación de difusividad. Aplica el principio de
superposición que indica que las soluciones son aditivas. El principio de superposición al
contacto agua petróleo original es mostrado en la siguiente figura. Para usar superposición,
la curva se aproxima a una serie de pequeños incrementos de presión de modo que los
pasos de tiempo sean pequeños para aproximar la curva con una recta. Los cambios de
presión de un intervalo al otro se calculan con el valor promedio del comienzo y el fin de
los intervalos de presión. Ver. Fig. 6.6.
Generalizando;
1
∆Pj =
2
( Pj −2 − Pj )
Para j = 1;
1
∆P1 = ( P0 − P1 )
2
Cuando la ecuación de Hurst and Van Everdingen se arregla para permitir variaciones de
presión queda:
n
We = Cv ∑ ∆p j Q ⎡⎣ A ( tn − t j −1 ) ⎤⎦
j =1
0.00633k
A=
φ µ wce ryto
2
n
We = Cv ∑ ∆p j Q ⎡⎣ A ( n − j + 1) ∆t ⎤⎦
j =1
n
Cv ∑ ∆p j Q ⎡⎣ A ( n − j + 1) ∆t ⎤⎦
j =1
Na = N +
D
Esta ecuación tiene tres incógnitas a saber: N, Cv y A. Estas dos últimas son las constantes
de Hurst and Van Everdingen. Los valores de la entrada de agua adimensional se
encuentran tabulados.
Radiales Infinitos5
Los resultados de este modelo se aproximan bastante a los del método de Hurst y Van
Everdingen para acuíferos finitos. Sin embargo, la teoría de Fetkovich es más simple y su
aplicación más fácil. Se utiliza un procedimiento de ensayo y error para evaluar las
constantes de intrusión de agua con datos de producción. El índice de productividad es
adecuado para describir la intrusión de agua desde un acuífero finito a un yacimiento de
hidrocarburos; la rata de intrusión de agua es directamente proporcional a la caída de
presión entre la presión promedia del acuífero y la presión en el contacto agua-aceite1,3,6.
Tabla 6.2. Funciones QtD a presión constante de Hurst y van Everdingen para Yacimientos
Radiales Finitos5
rD=1.5 rD=2 rD=2.5 rD=3
tD QtD tD QtD tD QtD tD QtD
5.0E -2 0.276 5E-2 0.278 0.1 0.408 0.3 0.755
6.0E -2 0.304 7.5E-2 0.345 0.15 0.509 0.4 0.895
7.0E -2 0.33 1E-1 0.404 0.2 0.599 0.5 1.023
8.0E -2 0.54 1.25E-1 0.458 0.25 0.681 0.6 1.143
9.0E -2 0.375 1.5E-1 0.507 0.3 0.758 0.7 1.256
0.1 0.395 1.75E-1 0.553 0.35 0.829 0.8 1.363
0.11 0.414 2E-1 0.597 0.4 0.897 0.9 1.465
0.12 0.431 2.25E-1 0.638 0.45 0.962 1 1.563
0.13 0.446 2.5E-1 0.678 0.5 1.024 1.25 1.791
0.14 0.461 2.75E-1 0.715 0.55 1.083 1.5 1.997
0.15 0.474 3E-1 0.751 0.6 1.14 1.75 2.184
0.16 0.486 3.25E-1 0.785 0.65 1.195 2 2.353
0.17 0.497 3.50E-01 0.817 0.7 1.248 2.25 2.507
0.18 0.507 3.75E-1 0.848 0.75 1.299 2.5 2.646
0.19 0.517 4E-1 0.877 0.8 1.348 2.75 2.772
0.2 0.525 0.425 0.905 0.85 1.395 3 2.886
0.21 0.533 0.45 0.932 0.9 1.44 3.25 2.99
0.22 0.541 0.475 0.958 0.95 1.484 3.5 3.084
0.23 0.548 0.5 0.983 1 1.526 3.75 3.17
0.24 0.554 0.55 1.028 1.1 1.605 4 3.247
0.25 0.559 0.6 1.07 1.2 1.679 4.25 3.317
0.26 0.569 0.65 1.108 1.3 1.747 4.5 3.381
0.28 0.574 0.7 1.143 1.4 1.811 4.75 3.439
0.3 0.582 0.75 1.174 1.5 1.87 5 3.491
0.32 0.588 0.8 1.203 1.6 1.924 5.5 3.581
0.34 0.594 0.9 1.253 1.7 1.975 6 3.656
0.36 0.599 1 1.295 1.8 2.022 6.5 3.717
0.38 0.603 1.1 1.33 2 2.106 7 3.767
0.4 0.606 1.2 1.358 2.2 2.178 7.5 3.809
0.45 0.613 1.3 1.382 2.4 2.241 8 3.843
0.5 0.617 1.4 1.402 2.6 2.294 9 3.894
0.6 0.621 1.6 1.432 2.8 2.34 10 3.928
0.7 0.623 1.7 1.444 3 2.38 11 3.951
0.8 0.624 1.8 1.453 3.4 2.444 12 3.967
. 2 1.468 3.8 2.491 14 3.985
2.5 1.487 4.2 2.525 16 3.993
3 1.495 4.6 2.551 18 3.997
4 1.499 5 2.57 20 3.999
5 1.5 6 2.599 22 3.999
7 2.613 24 4
8 2.619
9 2.622
10 2.624
A muy pequeña
Na
m = Cv
A muy grande
n ∆Pj A ⎡⎣tn − t j −1 ⎤⎦
∑j =1 D
∂We
qw = J ( Pa − Pf ) = (6.1)
∂t
donde:
We = Va ce ( Pi − Pa ) (6.2)
∂Pa
−Va ce = J ( Pa − Pf )
∂t
⎛ Jt ⎞
−⎜ ⎟
⎝ Va ce ⎠
Pa = Pf + ( Pi − Pf )e (6.3)
⎛ Jt ⎞
1 ∂We −⎜
Vc
⎟
Pf + = Pf + ( Pi − Pf )e ⎝ a e ⎠
J ∂t
⎡ ⎛ Jt ⎞
−⎜ ⎟⎤
We = ( Pi − Pf )Va ce 1 − e ⎝ a e ⎠ ⎥
⎢ Vc
⎢ ⎥
⎣ ⎦
⎡ ⎟⎤
⎛ J ∆tn ⎞
−⎜
= ( Pa ,n −1 − Pf ,n )Va ce 1 − e ⎝ a e ⎠ ⎥
⎢ Vc
∆We, n (6.4)
⎢ ⎥
⎣ ⎦
donde:
⎡ ⎟⎤
⎛ J ∆tn ⎞
−⎜ ⎡ ⎟⎤
⎛ J ∆t ⎞
−⎜
⎢
Va ce 1 − e ⎝ Va ce ⎠ ⎥
= Va ce 1 − e ⎝ a e ⎠ ⎥ = cte
⎢ Vc
⎢ ⎥ ⎢ ⎥
⎣ ⎦ ⎣ ⎦
∆We, n = A( Pa , n −1 − Pf ,n ) (6.5)
siendo A:
⎡ ⎟⎤
⎛ J ∆t ⎞
−⎜
A = Va ce ⎢1 − e a e ⎠ ⎥
⎝ Vc
⎢ ⎥
⎣ ⎦
⎡ n −1 ⎤
⎢ ∑ ∆We, j ⎥
Pa , n −1 = Pi ⎢1 − j =1 ⎥ (6.6)
⎢ Va ce Pi ⎥
⎢ ⎥
⎣ ⎦
Pf ,n = ( Pf , n −1 + Pf ,n ) / 2
La Ec. 6.5 proporciona la presión promedia del acuífero al final del paso de tiempo n-1.
Escribiendo una ecuación similar para la presión promedia del acuífero al final del paso de
tiempo n:
⎡ n
⎤
⎢ ∑ ∆We, j ⎥
Pa , n = Pi ⎢1 − j =1 ⎥ (6.7)
⎢ Va ce Pi ⎥
⎢ ⎥
⎣ ⎦
n
1
Pa , n = Pi −
Va ce
∑ ∆W
j =1
e, j
1 n
Pa , n = Pi − ∑ A( Pa, j −1 − Pf , j )
Va ce j =1
Tabla 6.3. Índices de productividad para acuíferos radiales y lineales1
⎛ n n
⎞
Pa , n = Pi − B ⎜ ∑ Pa , j −1 − ∑ Pf , j ⎟
⎝ j =1 j =1 ⎠
donde:
A
B=
Va ce
⎛ J ∆t ⎞
−⎜ ⎟
⎝ Va ce ⎠
B = 1− e
We ,n = Va ce ( Pi − Pa , n )
∆We,n = We,n − We ,n −1
Los índices de productividad se pueden estimar con las ecuaciones de la tabla 6.3.
EJEMPLO
Verifique los resultados del método de Fetkovich con el de Van Everdingen and Hurst para
un acuífero finito6. Dados los siguientes datos (ver además tabla 6.4):
Tabla 6.4. Datos de presión en el contacto agua-petróleo6
Cv = 1.119 α h φ ce ryto
2
= 1.119 (1)20)(0.2)(7 ×10−6 )(50002 ) = 783.3
0.00633k 0.00633(500)
A= = = 0.113 dias −1
φ µ wce ryto (0.2)(0.8)(7 ×10 )(5000 )
2 −6 2
Como la presión se da en términos anuales, el paso de nivel de tiempo es 365 días, luego:
π (rac2 − ryto
2
)hφ ce π (400002 − 50002 )(20)(0.2)(7 ×10−6 )
Va ce = = = 24673.95 bbl
5.615 5.615
Usando estos valores, se determina la intrusión de agua por el método de Hurst y van
Everdingen. Ver Tabla 6.5. Por medio del método de Fetkovich, el parámetro “B” se
calcula mediante ensayo y error usando el siguiente procedimiento:
Tabla 6.5. Volumen de agua intruida obtenida por el método de Hurst y van Everdingen
n
1 n
∑P
j =1
f,j = ∑ ( Pf , j −1 + Pf , j )
2 j =1
donde Pf,0 = Pi
∑P
j =1
a , j −1
De la Ec. 6.3:
⎛ Jt ⎞
−⎜ ⎟
⎝ Va ce ⎠
Pa , j −1 = Pf , j −1 + ( Pi − Pf , j −1 )e
Para el primer valor. Los demás valores se obtienen consecutivamente. Ver ejemplo. Note
que J y Vace pueden ser calculados. Esta sumatoria se evalúa si existen valores previos de
tiempo. Para n = 1, Pa ,0 = Pi
⎛ n n
⎞
Pa , n = Pi − B ⎜ ∑ Pa , j −1 − ∑ Pf , j ⎟
⎝ j =1 j =1 ⎠
5. Calcule Vace
(We ) H &VE
(Va ce ) n =
Pi − Pa , n
6. Determine el valor promedio de (Vace)n. Este valor variará debido a (1) incertidumbre en
la teoría, y (2) inexactitudes en los datos de campo. Por lo tanto, las diferencias en este
valor se deben primariamente a que la teoría no es confiable. Con valores correcto de B las
desviaciones son mínimas.
Los cálculos se realizan de la siguiente manera. Para el nivel de tiempo número 1 estima
Pa,n de la siguiente ecuación:
⎛ Jt ⎞ ⎛ J ∆t ⎞
−⎜ ⎟ −⎜ ⎟
⎝ Va ce ⎠ ⎝ Va ce ⎠
Pa , n = Pf ,n + ( Pi − Pf ,n )e = Pf , n + ( Pi − Pf , n )e
Acumule el valor de Pa , j −1 , para el primer caso será 2976.5 + 3000 = 5976.6 y obtenga Pa,n
de la siguiente expresión:
⎛ n n
⎞
Pa , n = Pi − B ⎜ ∑ Pa , j −1 − ∑ Pf , j ⎟ = 3000 − 0.6075(5976.6 − 5863) = 2930.988 ≈ 2931 psi
⎝ j =1 j =1 ⎠
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11)
n Pf,n Pf ,n
n n (5)-(4) B*(6) Pa , n We(H&VE) (Vace)n We(fetk)
∑P
j =1
f ,n
∑P a , j −1
j =1
0 3000 0.0
1 2923 2961.5 2961.5 3000 38.5 23.4 2976.6 581.4 24.8581 576.739
2 2880 2901.5 5863.0 5976.6 113.6 69.0 2931.0 1705.2 24.7087 1701.75
3 2848 2864.0 8727.0 8907.6 180.6 109.7 2890.3 2699.9 24.6084 2705.43
4 2821 2834.5 11561.5 11797.9 236.4 143.6 2856.4 3532.3 24.596 3541.32
5 2800 2810.5 14372.0 14654.3 282.3 171.5 2828.5 4219.7 24.6051 4228.92
6 2780 2790.0 17162.0 17482.8 320.8 194.9 2805.1 4798.0 24.6195 4805.65
7 2762 2771.0 19933.0 20287.9 354.9 215.6 2784.4 5310.6 24.6315 5316.48
8 2744 2753.0 22686.0 23072.3 386.3 234.7 2765.3 5782.1 24.6385 5786.86
9 2730 2737.0 25423.0 25837.6 414.6 251.9 2748.1 6207.9 24.6473 6210.8
10 2713 2721.5 28144.5 28585.7 441.2 268.0 2732.0 6608.2 24.6548 6609.27
11 2698 2705.5 30850.0 31317.7 467.7 284.1 2715.9 7005.3 24.6555 7006.25
12 2680 2689.0 33539.0 34033.6 494.6 300.5 2699.5 7408.6 24.6567 7409.22
13 2665 2672.5 36211.5 36733.1 521.6 316.9 2683.1 7814.1 24.6601 7813.68
14 2650 2657.5 38869.0 39416.2 547.2 332.4 2667.6 8197.9 24.661 8197.18
15 2633 2641.5 41510.5 42083.8 573.3 348.3 2651.7 8588.7 24.6603 8588.16
16 2620 2626.5 44137.0 44735.5 598.5 363.6 2636.4 8966.8 24.6619 8965.66
17 2607 2613.5 46750.5 47371.9 621.4 377.5 2622.5 9309.9 24.662 9308.71
18 2593 2600.0 49350.5 49994.4 643.9 391.2 2608.8 9647.3 24.6627 9645.76
19 2580 2586.5 51937.0 52603.2 666.2 404.7 2595.3 9982.2 24.6647 9979.82
20 2568 2574.0 54511.0 55198.5 687.5 417.7 2582.3 10301.0 24.6638 10298.9
24.6588
⎛ n n ⎞
Pa , n = Pi − B ⎜ ∑ Pa , j −1 − ∑ Pf , j ⎟ = 3000 − 0.6075(8907.6 − 8727) = 2890.25 psi
⎝ j =1 j =1 ⎠
⎛ n n ⎞
Pa , n = Pi − B ⎜ ∑ Pa , j −1 − ∑ Pf , j ⎟ = 3000 − 0.6075(11797.9 − 11561.5) = 2856.4 psi
⎝ j =1 j =1 ⎠
7. Como una medida de variación entre los dos métodos, utilice la siguiente función de
desviación (Esta función depende del valor asumido de B):
1 n ⎪⎧ ⎡ (We ) Fetk ⎤ ⎪⎫
Desv( B) = ∑ ⎨abs ⎢1 −
n j =1 ⎩⎪
⎥⎬
⎣ (We ) H &VE ⎦ ⎭⎪ j
⎛ P ⎞
P = − ⎜ i ⎟ We + Pi
⎝ Wei ⎠
Siendo P la presión promedia del acuífero después de haber extraído We barriles de agua,
Pi es la presión inicial del acuífero y Wei es la cantidad de agua in-situ inicial o
potencialmente intruible a la presión inicial.
qw Bw = J ( P − PR ) ma
Wei ⎛ − i ⎞
JPt
We = ( Pi − PR ) ⎜1 − e Wei
⎟
Pi ⎜ ⎟
⎝ ⎠
donde:
⎛ θ ⎞
⎟ π (rac − ryto )hφ Pi
2 2
ct ⎜
Wei = ⎝
360 ⎠
5.615
Wei ⎛ JP ∆t
− i n ⎞
∆We,n = ( Pn −1 − PR ,n ) ⎜1 − e Wei ⎟
Pi ⎜ ⎟
⎝ ⎠
⎛ W ⎞
Pn −1 = Pi ⎜1 − e ⎟
⎝ Wei ⎠
PR ,n = ( PR ,n −1 + PR , n ) / 2
EJEMPLO1
Calcule la intrusión de agua al tercer y cuarto año de producción para un yacimiento que
fallado (θ = 180), φ = 20.9 %, k = 275 md (para acuífero y yacimiento), µ = 0.25 cp, ct =
6x10-6 /psi, h = 19.2 ft, Ayto = 1216 ac, Aac = 250000 ac.
SOLUCIÓN
1216(43560)
ryto = = 5807 ft
0.5π
250000(43560)
ryto = = 83263 ft
0.5π
⎛ θ ⎞ ⎛ 180 ⎞
0.00708kh ⎜ ⎟ 0.00708(275)(19.2) ⎜ ⎟
J= ⎝ 360 ⎠ = ⎝ 360 ⎠ = 39.08
µ ⎡⎣ln(rac / ryto ) − 0.75⎤⎦ 0.25 [ ln(83263 / 5807) − 0.75]
⎛ θ ⎞
⎟ π (rac − ryto )hφ Pi
2 2
ct ⎜
Wei = ⎝
360 ⎠
5.615
⎛ 180 ⎞
6 × 10−6 ⎜ ⎟ π (83263 − 5807 )(19.2)(0.209)3793
2 2
Wei ⎛ JP ∆t
− i n ⎞
∆We,n = ( Pn −1 − PR ,n ) ⎜1 − e Wei ⎟
Pi ⎜ ⎟
⎝ ⎠
176.3MM ⎛ −
(39.08)(3793)(91.3)
⎞
∆We, n = ( Pn −1 − PR ,n ) ⎜1 − e 176.3 MM
⎟ = 3435( Pn −1 − PR ,n ) (A)
3793 ⎝ ⎠
⎛ W
Pn −1 = Pi ⎜1 − e
⎞ ⎛
⎟ = 3793 ⎜⎜1 −
∑ ∆We ⎞
⎟⎟ (B)
⎝ Wei ⎠ ⎝ 176.3MM ⎠
Resolviendo las Ecs. (A) y (B), se obtienen los resultados dados en la tabla 8. El
procedimiento se resume de la siguiente forma para el nivel 1:
2) Estime la diferencia entre Pn −1 − PR ,n =3793-3790.5 = 2.5 (columna 5). Estime ∆We con la
Ec. A.
∆We, n = 3435( Pn −1 − PR , n ) = 3435(3793 − 3790.5) = 8587.5 ≅ 8600 bbl
⎛
Pn = 3793 ⎜⎜1 −
∑ ∆We ⎞ ⎛
⎟⎟ = 3793 ⎜ 1 −
8600 ⎞
⎟ = 3792.82 psi
⎝ 176.3MM ⎠ ⎝ 176.3MM ⎠
⎛
Pn = 3793 ⎜⎜1 −
∑ ∆We ⎞ ⎛
⎟⎟ = 3793 ⎜ 1 −
49100 ⎞
⎟ = 3791.94 psi
⎝ 176.3MM ⎠ ⎝ 176.3MM ⎠
REFERENCIAS
CAPITULO 7
CURVAS DE DECLINACIÓN
7.1. INTRODUCCIÓN
El análisis de curvas de declinación podría ser una de las técnicas de ingeniería que más
están en desuso y al mismo tiempo parece ser una de las técnicas que menos atención
ofrece ya que ellas se aplican siempre y cuando las condiciones mecánicas del pozo y el
área de drene del yacimiento permanecen constantes. Sin embargo, el uso de curvas tipo
incluye soluciones que alivian los problemas en mención. Sin embargo, para hacer
predicciones del yacimiento debería emplearse dichos análisis. El típico análisis consiste en
graficar datos de producción contra tiempo en papel semilog e intentar ajustar estos datos
con una recta la cual se extrapola hacia el futuro. Las reservas se calculan con base en una
rata de producción promedia anual1-5.
Por muchos años, un gráfico de q vs. t para muchos pozos puede extrapolarse, lo cual se
convirtió en un arte. Es una de las técnicas menos usadas. Las reservas se calculan con base
en una producción promedia anual para las ratas de producción extrapoladas. La
declinación hiperbólica da mejores resultados. Sin embargo, puesto que es más difícil se
prefiere la armónica. Además, la excusa, es que la diferencia entre una y otra curva, con el
tiempo, no es muy significativa.
Rata de producción, q
q1
∆q
∆t
t1
Tiempo, t
⎛ ∆q / q ⎞
a = −⎜ ⎟ (7.1)
⎝ ∆t ⎠
D = ( qt =0 − qt =1 year ) (7.2)
a = − ln (1 − D ) (7.3)
Este tipo de curva de declinación parece ser la más usada por los ingenieros de
yacimientos, por su facilidad, e incluso cuando se es consciente que la declinación
hiperbólica describe mejor las características de la mayoría de los pozos. Es definida por
una función exponencial. Arreglando la ecuación (7.1)4:
⎛ ∆q ⎞
a∆t = − ⎜ ⎟
⎝ q ⎠
t q
∆q
a ∑ ∆t = −∑
0 qi q
Integrando;
qi q
at = ln = 2.303 log i
q q
q = qi e− at
t2
∆N p = ∑ q∆t
t1
q2
∆q
∆N p = −∑
q1 a
q1 − q2
∆N p =
a
Esta considera que la rata de declinación varía con el tiempo. Es buena para yacimientos
que producen por gas en solución. Esta técnica es muy consumidora de tiempo. La rata de
declinación varía así4:
n
a ⎛q⎞
=⎜ ⎟
ai ⎝ qi ⎠
⎛ ∆q ⎞
Si a = − ⎜ ⎟ entonces:
⎝ q∆t ⎠
− [ ∆q q∆t ]
n
⎛q⎞
=⎜ ⎟
ai ⎝ qi ⎠
Separando variables:
t q
ai ∫ dt = −qin ∫ q − ( n +1) dq
0 qi
⎛ q − n qin ⎞
ai t = − qin ⎜ − ⎟
⎝ −n −n ⎠
nai t = qin q − n − 1
qin
nai t = −1
qn
1
q = qi [ nai t + 1] n
De igual forma:
t2
∆N p = ∫ qdt
t1
t2
∆q
∆N p = − ∫
t1
a
n
⎛q⎞
Si a = ai ⎜ ⎟
⎝ qi ⎠
Entonces;
q2
⎛ qin ⎞ − n
∆N p = − ∫ ⎜ − ⎟ q dq
q1 ⎝
ai ⎠
q2
qn
∫q
−n
∆N p = i dq
ai q1
qin
H=
ai (1 − n )
∆N p = H ( q11− n − q12− n )
q = qi / [ nai t + 1]
q2
∆q
∆N p = −∑
q1 a
como se manifestó:
q
a= ai
qi
Entonces:
q2
∆q
∆N p = −∑
q
q1
ai
qi
qi q1
∆N p = ln
ai q2
No existen curvas tipo para declinación armónica debido a que ésta ocurre muy
esporádicamente.
Una forma más práctica en usar la declinación hiperbólica es comparar los datos reales de
declinación con curvas tipo, las cuales viene para varios valores de n y ai. Estas curvas son
deferentes que aquellas que se usan para análisis de presiones de fondo. Una vez, se ha
determinado cual curva es la que mejor se ajusta los datos de declinación, se han
determinado los valores de n, ai y qi.
EJEMPLO
La tabla 7.1 presenta los datos de producción para un pozo de crudo. Cuál será la rata de
producción a los 5 años? Cuál es la vida del pozo a Jun-82 si el límite económico es 1 BPD.
100
n=0.3
ai = 0.0025/mes
ai = 0.0020/mes
ai = 0.0015/mes
ai = 0.0010/mes
ai = 0.0005/mes
10
Rata ó q/q i
0.1
Dic -120
-110
-100
-90
-80
-70
-60
-50
-40
-30
-20
-10
10
20
30
40
50
60
70
80
Jun 90
100
110
Dic 120
Jun
Dic
Dic
Jun
Jun
Dic
Tiempo, meses
100
n=0.5
ai = 0.0025/mes
ai = 0.0020/mes
ai = 0.0015/mes
ai = 0.0010/mes
ai = 0.0005/mes
10
Rata ó q/q i
0.1
Dic -120
-110
-100
-90
-80
-70
-60
-50
-40
-30
-20
-10
10
20
30
40
50
60
70
80
Jun 90
100
110
120
Jun
Dic
Dic
Dic
Jun
Jun
Dic
Tiempo, meses
100
n=0.7
ai = 0.0025/mes
ai = 0.0020/mes
ai = 0.0015/mes
ai = 0.0010/mes
ai = 0.0005/mes
10
Rata ó q/q i
0.1
Dic -120
-110
-100
-90
-80
-70
-60
-50
-40
-30
-20
-10
10
20
30
40
50
60
70
80
90
100
110
120
Jun
Jun
Dic
Dic
Dic
Jun
Jun
Dic
Tiempo, meses
100
b/d
10
i
Rata ó q/q
1 Años
0.1
Dic -120
-110
-100
-90
-80
-70
-60
-50
-40
-30
-20
-10
10
20
30
40
50
60
70
80
Jun 90
100
110
120
Jun
Dic
Dic
Dic
Jun
Dic
Jun
Tiempo, meses
10
n = 0.0
n = 0.1
n = 0.2
n = 0.3
n = 0.4
n = 0.5
n = 0.6
n = 0.7
n = 0.8
n = 0.9
q i/q
n = 1.0
1
0.1
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
∆Np / qi t
Fig. 7.6. Curva tipo que relaciona rata de producción con producción acumulada4
n = 0.0
0.1 n = 0.1
n = 0.2
n = 0.3
n = 0.4
i
q/q
n = 0.5
n = 0.6
n = 0.7
n = 0.8
n = 0.9
0.01
n = 1.0
q 1
= p ara n > 0
(1 + n a i t )
1/ n
qi
Exponen
Ar
m
ón
q 1
ica
= ai t para n = 0
qi e cial
0.001
0.01 0.1 1 10 100 1000
ai t
Fig. 7.7. Curvas tipo de Arps4
1000
100
q, STB/D
10
1
Jul-72 Dec-73 Apr-75 Aug-76 Jan-78 May-79 Oct-80 Feb-82 Jul-83
Tiempo, anual
SOLUCIÓN
Para resolver este problema refiérase a la Fig. 7.8. A Jul-83 el caudal (extrapolado) es de 10
a Dic-76, q = 42 BPD. Entre estos dos valores hay un lapso de 80 meses. La constante de
declinación es:
La vida del pozo para un límite económico de 1 BPD (30.4 bbl/mes) es:
EJEMPLO
Asuma que el pozo del ejemplo anterior se fracturó en Jun-82 y la rata pasó a 52 BPD.
Asuma, además, que el pozo entró a declinación constante, cuánto es el petróleo recuperado
después de Jun-82? Cuál es el incremento en el recobro y el cambio en la vida del pozo
para el mismo límite económico4.
SOLUCIÓN
q1 − q2
∆N p =
a
(13 − 1)30.4
∆N p antes = = 20267 bbl
0.018
(52 − 1)30.4
∆N p despues = = 21533 bbl
0.072
Las reservas recuperadas se incrementaron tan solo (21533-20667) en 866 bbl pero fue más
rápido.
EJEMPLO
Usando las gráficas para n =0.3, 0.5 y 0.7 halle la vida remanente del pozo y las reservas si
el límite económico es de 10 bbl/mes.
100
n=0.5
ai = 0.0025/mes
ai = 0.0020/mes
ai = 0.0015/mes
ai = 0.0010/mes
1000 ai = 0.0005/mes
10
Rata ó q/q i
100
10
0.1
Dic -120
-110
-100
-90
-80
-70
-60
-50
-40
-30
-20
-10
10
20
30
40
50
60
70
80
Jun 90
100
110
120
Jun
Dic
Dic
Dic
Jun
Dic
Jun
Tiempo, meses
1
Mar-71 Jul-72 Dec-73 Apr-75 Aug-76 Jan-78 May-79 Oct-80 Feb-82 Jul-83
El mejor ajuste se obtuvo para la gráfica de n = 0.5. Ver Fig. 7.9. Del ajuste, se tiene que ai
= 0.0015/mes. La vida útil del pozo será:
1
q = qi [ nai t + 1] n
1
10(30.4) = 510(30.4) [ 0.5(0.0015)t + 1] 0.5
REFERENCIAS
CAPITULO 8
CONIFICACIÓN Y DIGITACIÓN
8.1. GENERALIDADES
Por otra parte, si la producción de agua ocurre en la vida temprana de un yacimiento sin
buzamiento donde el contacto se vuelve inestable y el agua fluye más rápido que el crudo
debido a su baja viscosidad y a la permeabilidad horizontal, este fenómeno se llama
digitación o lengüeteo.
Esta sección está dedicada a presentar estos dos procesos que conducen a la ruptura
prematura de agua. Estos fenómenos son muy importantes porque causan un decremento en
la rentabilidad del proyecto en diversas formas. Primero, la productividad de crudo se
reduce debido a efectos de permeabilidad relativa, segundo, los costos de levantamiento
aumentan como resultado de un fluido promedio más pesado y la inyección y/o
acondicionamiento de agua es sustancialmente cara. Tercero, la eficiencia de recobro se
reduce porque el límite económico del corte de agua se alcanza con petróleo remanente
producible en el área de drene del pozo2-3.
Esto se debe a que la caída de presión que causa el flujo o producción de agua es mayor
que la fuerza gravitacional:
La ecuación (8.1) no habla de tiempo. Sin embargo, expresa que a cualquier tiempo esta
inigualdad se satisface ocurriendo instantáneamente la conificación. También cuando la
inigualdad no se satisface no habrá conificación puesto que la presión de fondo fluyente es
controlable mediante regulación de la rata de flujo. Luego, la conificación es controlable.
Aceite y agua
Consideraciones:
1) Flujo radial
2) El agua fluye verticalmente desde el contacto agua petróleo al fondo del pozo
3) El gradiente de presión controlante de la producción se restringe al gradiente
gravitacional debido a la diferencia de densidades entre los fluidos y la distancia entre
el contacto agua petróleo y las perforaciones más bajas. Con base en estas
consideraciones se halla la siguiente ecuación:
0.001535 ( ρ w − ρo ) k (h 2 − D 2 )
qc =
⎛ re ⎞
µo Bo ln ⎜ ⎟
⎝ rw ⎠
qc = Caudal crítico de petróleo o máxima rata de petróleo para precluir producción de agua,
BF/D
ρ = Densidad de los fluidos, gr/cc
k = Permeabilidad de la formación
h = Espesor de la formación, pie
D = Intervalo perforado, pie
EJEMPLO
Halle el qc con los siguientes datos: k = 100 md, h = 50 pies, D = 10 ft, densidad del agua =
1.05, densidad del aceite = 0.8 gr/cc, Viscosidad del petróleo = 1, factor volumétrico de
formación del aceite = 1.2 bbl/BF, radio del yacimiento 745 pies, radio del pozo = 0.25
pulgadas.
SOLUCIÓN
Aunque la rata calculada es una guía, es evidente que estas ratas son muy bajas para ser
económicamente rentable. Por lo que este método es muy restrictivo.
Aunque producir un pozo a ratas por debajo de la crítica evita la conificación, muchos
creen que hacerlo causa baja rentabilidad. Sobocinski y Cornelious estudiaron, mediante un
prototipo de plexiglás en forma de pastel empacado con arena, ver Fig. 8.2, el problema de
definir el tiempo requerido después de que el pozo esté puesto en producción para resultar
la conificación o para irrumpir cuando el pozo se produce a mayor caudal que el crítico. A
medida que ellos cambiaban el color del agua y del crudo de suministro pudieron
monitorear la posición del contacto agua-petróleo. Los resultados se presentan en la Fig.
8.3.
0.00307 ( ρ w − ρo ) k h hc
Z= (8.2)
µo Bo qo
MODELO DE PLEXIGLAS
Arena empacada
Producción Suministro
CRUDO de crudo
Irupción de agua
Cono de agua
AGUA Suministro
de agua
4.0
3.5
Z, Altura del cono, Adimensional
3.0
Curva de
ruptura
2.5
Curvas de
2.0 partida
1.5
1.0
0.0
0 1 2 3 4 5 6 7 8
Tiempo adimensional, t D
0.000137 ( ρ w − ρo ) k H (1 + M α ) t
tD = (8.3)
µo φ h FK
kH = Permeabilidad horizontal
FK = Relación entre kH y kV
h = Espesor zona de la aceite
hc = Altura del cono. Distancia inicial entre el contacto agua petróleo y las perforaciones
más bajas
t = Tiempo, días
α = Exponente experimental que depende de la relación de movilidades, M, α = 0.5 para M
< 1 y 0.6 para M ≥ 1.
t Dbt =Z
(16 + 7 Z − 3Z )2
(8.4)
28 − 8Z
La Fig. 8.3 presenta tres diferentes tipos de curvas: ruptura, restauración y partida. Para
determinar ruptura solamente se requiere la curva de ruptura. Las otras dos normalmente no
se usan, pero permiten estudiar la forma como el cono se construye. Después que se
determina Z, se obtiene un punto sobre la curva de ruptura. A partir de este punto, trace una
curva de partida (paralela a las mostradas) hacia la curva básica. Luego considere la
primera parte de la historia de construcción del cono, la curva básica se usa de tiempo
adimensional cero hasta alcanzar la curva de ruptura. La formación del cono es muy lenta
hasta alcanzar un punto cerca al pozo donde la velocidad del cono se acelera. El cono crece
cada vez más rápido hasta que ocurre ruptura.
Inspeccionando la Ec. 8.4 indica que para Z=3.5, el valor de tD se vuelve infinito. Este valor
puede usarse para hallar la rata crítica con la Ec. 8.2. El valor de Z que corresponde a un
tiempo de ruptura infinito corresponde a la rata crítica:
0.000887 ( ρ w − ρo ) k h hc
qc = (8.5)
µo Bo
EJEMPLO
Determine el tiempo de producción hasta que la ruptura toma lugar para un pozo, de
acuerdo con la siguiente información: Permeabilidad horizontal es 50 md, h es 40 pies, la
viscosidad del petróleo es 1 cp, el factor volumétrico es 1.2 bbl/BF, la diferencia de
densidades es 0.2 gr/cc, la permeabilidad vertical es 5 md, D es 10 pies, la relación de
movilidades es 1.5, la porosidad es 0.2 y el caudal de petróleo es 50 BF/D.
SOLUCIÓN
Como M > 1, α = 0.6. La relación de permeabilidades Fk = 50/5 = 10. La altura del cono
hc=40-10 = 30 ft. Z se halla por medio de la Ec. 8.2:
0.00307 ( ρ w − ρ o ) k h hc 0.00307(0.2)(50)(40)(30)
Z= = = 0.61
µo Bo qo (1)(1.2)(50)
Con Z = 0.61, se usa la Fig. 8.3 para hallar tD = 0.5. Despejando el tiempo de la Ec. 8.3 se
tiene:
µo φ h FK t D (1)(0.2)(40)(10)(0.5)
t = = = 1283 dias = 3.5 años
0.000137 ( ρ w − ρo ) k H (1 + M ) 0.000137(0.2)(50) (1 + 1.50.6 )
α
EJEMPLO
Determine la rata crítica con la Ec. 8.5 para el problema dado en la sección 8.2.1.
SOLUCIÓN
0.000887 ( ρ w − ρo ) k h hc 0.000887(0.25)(100)(50)(40)
qc = = = 36.5 STB / D
µo Bo (1)(1.2)
El problema de la sección 8.2.1 dio un valor de 9.6 STB/D. Tracy sugiere que se use un
valor entre la rata crítica de Chaney y otros (sección 8.2.3) y la rata crítica de Sobocinski y
Cornelious.
Los resultados se presentaron gráficamente en las Figs. 8.4.a a 8.4.e. Estas figuras son para
diferentes espesores. Las figuras presentan el caudal crítico en función de la distancia desde
el tope de las perforaciones hasta el tope de la formación o al contacto gas-aceite. Estas
curvas se prepararon con una permeabilidad de 1000 md, una viscosidad de crudo de 1 cp y
un diferencial de gravedad entre crudo y agua de 0.3 gr/cm. La ecuación de la rata crítica
es:
100
1.25 ft
Conificación
2.5 ft
de gas
3.75 ft
5 ft
6.25
1.25 ft
10
2.5 ft
3.75 ft
q curva
5 ft
6.25
Conificación
de agua
0.1
0 2 4 6 8 10 12
1000
2.5 ftft
1.25
Conificación
2.5 ft
5.0
de gas
3.75ftft
7.5
5 ftft
10
100 6.25 ft
12.5
2.5 ftft
1.25
2.5 ft
5.0
3.75ftft
7.5
q curva
5 ftft
10
10 6.25 ft
12.5
Conificación
Conificación
de agua
de agua
0.1
0 5 10 15 20 25
1000
Conificación
de gas
5 ft
10 ft
100 15 ft
20 ft
q curva
25 ft
5 ft
10 ft
15 ft
10 20 ft
25 ft
Conificación
de agua
1
0 5 10 15 20 25 30 35 40 45 50
10000
7.5 ft
15 ft
22.5 ft Conificación
30 ft de gas
37.5 ft
1000
q curva
7.5 ft
100 15 ft
22.5 ft
30 ft
37.5 ft
10
Conificación
de agua
1
0 10 20 30 40 50 60 70 80
10000
10 ft
20 ft
Conificación
30 ft
de gas
40 ft
50 ft
1000
q curva
10 ft
100 20 ft
30 ft
40 ft
50 ft
10
Conificación
de agua
1
0 10 20 30 40 50 60 70 80 90 100
0.00333qcurva k ( ρ w − ρo )
qc = (8.6)
µo Bo
Un método alternativo es usar la Ec. 8.7 que tiene un error dentro del 5 %. Para h = 25 ft el
error es ligeramente mayor:
2 2
⎛ 57 − h ⎞ ⎛D ⎞⎛ h − 10 ⎞
qcurve = 0.1313(h − D ) + 34 ⎜
2 2
⎟ − 250 ⎜ − 0.3 ⎟⎜ ⎟ − 40 (8.7)
⎝ 44 ⎠ ⎝h ⎠⎝ 90 ⎠
EJEMPLO
Use los datos del problema de la sección 8.2.1 para determinar el caudal crítico.
SOLUCIÓN
Puesto que hay 50 ft desde el tope de la zona de petróleo hasta el contacto gas-petróleo, se
usa la Fig. 8.4.c. Se asume que el intervalo completado está al tope de la zona de petróleo,
por lo tanto el valor de entrada en el eje x es cero. Cortando con la curva para un intervalo
perforado de 10 ft. De allí qcurva = 280 bbl/D. Insertando este valor en la Ec. 8.6 se tiene:
0.00333qcurva k ( ρ w − ρo ) 0.00333(280)(100)(0.25)
qc = = = 19.4 STB / D
µo Bo (1)(1.2)
2 2
⎛ 57 − 50 ⎞ ⎛ 10 ⎞⎛ 50 − 10 ⎞
qcurve = 0.1313(50 − 10 ) + 34 ⎜
2 2
⎟ − 250 ⎜ − 0.3 ⎟⎜ ⎟ − 40 = 280.9 B / D
⎝ 44 ⎠ ⎝ 50 ⎠⎝ 90 ⎠
Por lo tanto:
Ellos le hicieron dos cambios al método de Sobocinski y Cornelious: (1) Se desarrolló una
ecuación para reemplazar la altura del cono vs. el tiempo adimensional, (2) Encontraron
que α = 0.7 para 0.14 ≤ M ≤ 7.3. Las nuevas ecuaciones son:
0.00307 ( ρ w − ρo ) k H h hc
Z= (8.2)
µo Bo qo
Z
t Dbt = (8.8)
3 − 0.7 Z
µo φ h FK tDbt
t = (8.9)
0.000137 ( ρ w − ρ o ) k H (1 + M α )
Los tiempos de ruptura de este método son menores que los predichos por Sobocinski y
Cornelious. Para valores de Z entre 0.1 y 3, los tiempos de ruptura de este método están
entre el 40 y 50 % del método de Sobocinski y Cornelious.
Inspeccionando la Ec. 8.8 puede verse que el tiempo adimensional se hace infinito cuando
Z = 4.28. Luego la rata crítica será:
0.000717 ( ρ w − ρo ) k H h hc
qc = (8.10)
µo Bo
EJEMPLO
0.00307 ( ρ w − ρ o ) k h hc 0.00307(0.2)(50)(40)(30)
Z= = = 0.61
µo Bo qo (1)(1.2)(50)
Z 0.61
t Dbt = = = 0.237
3 − 0.7 Z 3 − 0.7.(0.61)
µo φ h FK t D (1)(0.2)(40)(10)(0.237)
t = = = 594 dias = 1.6 años
0.000137 ( ρ w − ρo ) k H (1 + M ) 0.000137(0.2)(50) (1 + 1.50.7 )
α
0.000717 ( ρ w − ρ o ) k H h hc 0.000717(0.25)(100)(50)(40)
qc = = = 29.9 STB / D
µo Bo (1)(1.2)
A=
( ρ w − ρo ) k (h 2 − D 2 ) (8.12)
2049µo Bo
π
B = 0.432 + (8.13)
⎛r ⎞
ln ⎜ e ⎟
⎝ rw ⎠
0.14
⎛h⎞
C =⎜ ⎟ (8.14)
⎝ re ⎠
EJEMPLO
Usando los datos del ejemplo de la sección 8.2.1 determine el caudal crítico por el método
de Kuo y Desbrisay.
SOLUCIÓN
A=
( ρ w − ρo ) k (h2 − D 2 ) = (0.25)(100) (502 − 102 ) = 24.402
2049µo Bo 2049(1)(1.2)
π π
B = 0.432 + = 0.432 + = 0.825
⎛r ⎞ ⎛ 745 ⎞
ln ⎜ e ⎟ ln ⎜ ⎟
⎝r ⎠ ⎝ 0.25 ⎠
0.14 0.14
⎛h⎞ ⎛ 50 ⎞
C =⎜ ⎟ =⎜ ⎟ 0.685
⎝ re ⎠ ⎝ 745 ⎠
Kuo y Desbrisay3 desarrollaron un método simple para desarrollar una correlación para
predecir el comportamiento del corte de agua en un yacimiento con empuje de agua de
fondo. Ellos usan el método de Bournazel y Jeanson como se describió en la sección 8.2.4.
Basado en sus resultados, ellos tomaron la producción de agua a comenzarse a la mitad del
tiempo predicho por Bournazel y Jeanson. Ellos definieron dos variables adimensionales:
t
tD = (8.15)
t Dbt
WC
(WC ) D = (8.16)
WClim
Mhw
WClim = (8.17)
Mhw + ho
Para evaluar WClim se hicieron otras consideraciones. El flujo vertical de agua a presión
constante se asumió tal que la producción del pozo causa movimiento equivalente de agua.
También se asume que el área seccional es constante. Con estas consideraciones,
presentaron la siguiente ecuación de balance da materia que se resuelve por ensayo y error.
⎡ ⎛N ⎞ 1 − S wc ⎤
ho = H o ⎢1 − ⎜ p ⎟ ⎥ (8.18)
⎣ ⎝ N ⎠ 1 − S wc − Sor ⎦
⎡ ⎛N ⎞ 1 − S wc ⎤
hw = H w + H o ⎢1 − ⎜ p ⎟ ⎥ (8.19)
⎣ ⎝ N ⎠ 1 − S wc − Sor ⎦
La Ec. 8.16 relaciona el corte real de agua adimensional. El corte de agua adimensional y el
tiempo se relacionan de la siguiente forma:
Slider presentó un método para analizar la inclinación del contacto agua petróleo como
resultado de la producción de pozos estructuralmente altos.
tan β =
( G − [ M − 1]) tan α (8.20)
G
C
WO
β
α WOCi
A B
0.488 ( γ w − γ o ) k A Sen α
qc = (8.22)
µo µw
−
kro krw
EJEMPLO
Determine el máximo caudal para crear un contacto agua petróleo estable con el tiempo. k
= 0.5 D, A = 50000 pie², diferencia de densidades = 0.2 gr/cc, ángulo α = 10, viscosidad del
petróleo 2 cp, kro = 1, krw = 0.25, y viscosidad del agua 0.4 cp.
SOLUCIÓN
Asumiendo que los pozos que están en la parte superior de la estructura producen solo
crudo, entonces:
qc 1061.4
qsup = = = 786.2 STB / D
Bo 1.35
REFERENCIAS
APENDICE A
MÉTODO DE MUSKAT3,6,10,13
Este método fue presentado por Muskat en 1945. El método es aplicable a yacimientos de
petróleo que producen por gas en solución con la posibilidad de incluir capa de gas e
inyección de gas. Este método hace uso de la ecuación de balance de materiales en forma
diferencial y se aplica a volúmenes pequeños del yacimiento con gradientes de presión
despreciables.
Muskat consideró el yacimiento como un medio poroso homogéneo a lo largo del cual la
presión es uniforme. Lo comparó con un tanque con válvulas de salida distribuidas
continua y uniformemente utilizadas para drenar fluidos. Cada elemento de volumen del
yacimiento produce a través de su propia salida y no existe intercambio de fluido entre los
elementos de volumen. El comportamiento del yacimiento total se determina a partir del
comportamiento de cualquiera de los elementos de volumen que forman el yacimiento.
Considere una unidad de volumen poroso (un barril) en un yacimiento de volumétrico sin
capa de gas. Considere además que el petróleo remanente después del período de
producción está uniformemente distribuido. Si So representa la saturación de petróleo a
cualquier tiempo durante la producción del yacimiento, el volumen de petróleo a
condiciones del yacimiento, por barril de espacio poroso, será So barriles y a condiciones
normales será So/ßo. El volumen de gas se da en pies cúbicos normales por barril de a
condiciones de yacimiento y está determinado por el gas en solución más el gas libre, que
matemáticamente se expresa:
So 1
Rs + (1 − S o − S w ) (A.1)
βo βg
Los caudales de aceite y gas a condiciones normales se expresan como el cambio del
volumen de petróleo y gas existente en el espacio poroso unitario con el tiempo. Si Qo es
el caudal de petróleo y Qg el de gas, a condiciones normales, se tiene:
d ⎛ So ⎞
Qo = ⎜ ⎟ (A.2)
dt ⎜⎝ β o ⎟⎠
d ⎡ So 1 ⎤
Qg = ⎢ Rs + (1 − S o − S w ) ⎥ (A.3)
dt ⎢⎣ β o β g ⎥⎦
La razón gas-petróleo instantánea, Ri, resulta de dividir el caudal de gas por el caudal de
petróleo, así:
d ⎡ So 1 ⎤
⎢ Rs + (1 − S o − S w ) ⎥
Qg dt ⎢⎣ β o β g ⎥⎦
Ri = = (A.4)
Qo d ⎛ So ⎞
⎜ ⎟
dt ⎜⎝ β o ⎟⎠
En esta ecuación los diferenciales están dados en función del tiempo, pero pueden
expresarse en función de presión ya que ésta variable independiente. Por lo tanto, puede
escribirse:
d d dp
= (A.5)
dt dp dt
d ⎡ So 1 ⎤
⎢ Rs + (1 − S o − S w ) ⎥
dp ⎣⎢ β o β g ⎦⎥
Ri = (A.6)
d ⎛ So ⎞
⎜ ⎟
dp ⎜⎝ β o ⎟⎠
S o dRs d ⎛ 1 ⎞
+ (1 − S o − S w ) ⎜ ⎟
β o dp dp ⎜⎝ β g ⎟
⎠
Ri = Rs + (A.7)
1 ⎛ dS o dβ ⎞
⎜β
2 ⎜ o
− S o o ⎟⎟
β o ⎝ dp dp ⎠
k g µo β o
Ri = Rs + (A.8)
ko µ g β g
S o dRs d ⎛ 1 ⎞
+ (1 − S o − S w ) ⎜ ⎟
k g µo β o β o dp dp ⎜⎝ β g ⎟
⎠
= Rs + (A.9)
ko µ g β g 1 ⎛ dS o dβ ⎞
⎜β
2 ⎜ o
− S o o ⎟⎟
β o ⎝ dp dp ⎠
dS o
Despejando ;
dp
S o β g dRs d ⎛ 1 ⎞
⎟ β g + S o µ o g dβ o
k
+ (1 − S o − S w ) ⎜
dS o β o dp dp ⎜⎝ β g ⎟
⎠ β o µ g ko dp
= (A.10)
dp k g µo
1+
ko µ g
Defina:
β g dRs
λ= (A.11)
β o dp
d ⎛ 1⎞
σ = βg ⎜ ⎟ (A.12)
dp ⎜⎝ β g ⎟⎠
1 µ o dβ o
η= (A.13)
β o µ g dp
kg
ψ= (A.14)
ko
dS o
y reemplace estos términos en la ecuación de
dp
dS o S o λ + (1 − S o − S w )σ + S oψ η
= (A.15)
dp µ
1 +ψ o
µg
⎡ ⎤
⎢ S λ + (1 − S − S )σ + S ψ η ⎥
∆S o = ∆p ⎢ o o w o ⎥ (A.16)
⎢ µo ⎥
⎢ 1 +ψ ⎥
⎣ µg ⎦
Para calcular el valor de ∆So, correspondiente a un valor dado de ∆P = P1 - P2, los valores
de λ ,σ ,η y µo/µg deben evaluarse a la presión promedia del intervalo, o sea, a (P1+P2)/2
ya que para intervalos pequeños de ∆P, puede asumirse una variación lineal entre P1 y P2
para un decremento de saturación ∆So = ( So1 - So2). La determinación de kg/ko y (1-So-
Sw) requiere del conocimiento de la saturación promedio correspondiente al decremento de
presión dado o mediante el uso de Correlaciónes empíricas como se efectúa en el presente
trabajo. Este valor también podría determinarse por tanteo pero es demasiado arriesgado y
poco recomendable. Si los decrementos de presión son pequeños puede, tomarse el valor de
saturación de petróleo al comienzo del intervalo, es decir, el correspondiente a la presión
P1, obteniéndose buenos resultados. Esto se hace normalmente en la práctica. Sin embargo
es bueno tener presente que existe determinado error. Además los errores por este concepto
son acumulativos. Si se quiere obviar este error, Muskat recomienda usar el método de
Runge - Kutta para tener más precisión en la solución numérica de la ecuación diferencial.
A.2. PROCEDIMIENTO
El cálculo de predicción por el método de Muskat puede resumirse en los siguientes pasos:
Np2 So2 β oi
= 1− (A.17)
N (1 − Sw ) β o 2
⎡ ⎛ Ri I ⎞ ⎤
⎢ S o λ + (1 − S o − S w )σ + S oη ⎜ψ − ⎟⎥
⎢ ⎝ α ⎠⎥
∆S o = ∆p (A.18)
⎢ µ ⎛ Ri I ⎞ ⎥
⎢ 1 + o ⎜ψ − ⎟ ⎥
⎣ µg ⎝ α ⎠ ⎦
βo µo
α= (A.19)
βg µg
⎛ N p ⎞ βo
S o = (1 − S w )⎜⎜1 − ⎟ (A.20)
⎝ N ⎟⎠ β oi
⎡ ⎛ Ri I ⎞ ⎤
⎢ S o λ + [(1 + m )(1 − S w ) − S o ]σ + S oη ⎜ψ − ⎟⎥
⎢ ⎝ α ⎠⎥
∆S o = ∆p (A.21)
⎢ µ ⎛ Ri I ⎞ ⎥
⎢ 1 + o ⎜ψ − ⎟ ⎥
⎣ µg ⎝ α ⎠ ⎦
APENDICE B
METODO DE PIRSON3,6,10,13
B.1. GENERALIDADES
La ecuación de balance de materia para un yacimiento que produce únicamente por gas en
solución es:
Np (β o − Rs β g ) + β g Gp
N= (B.1)
β g (Rsi − Rs ) − (β oi − β o )
Considere que la producción es una fracción del petróleo original en el yacimiento (factor
de recobro) y que la cantidad de petróleo original, N, es 1. Dividiendo la Ec. (B.1) por ßg,
se tiene:
⎛β ⎞
N p ⎜ o − Rs ⎟ +Gp
⎜β ⎟
1= ⎝ g ⎠ (B.2)
⎛ βo ⎞ ⎛β ⎞
⎜ − Rs ⎟ − ⎜ oi − Rsi ⎟
⎜β ⎟ ⎜β ⎟
⎝ g ⎠ ⎝ g ⎠
Defina pj y pj+1 dos presiones sucesivas del yacimiento (pj > pj+1) para un periodo de
producción fraccional ∆ jj+ 1Np . Desde el inicio de la producción han ocurrido decrementos
Ri 0 + Ri1
Ri1 = (B.3)
2
Ri1 + Ri 2
Ri 2 = (B.4)
2
Ri j −1 + Ri j
Ri j = (B.5)
2
j
Gp j = ∑ ∆Np j Ri j (B.6)
j =1
⎛ βo ⎞ ⎛ β oi ⎞ ⎛ βo ⎞ j
⎜⎜ − Rs⎟⎟ − ⎜⎜ − Rsi⎟⎟ = Np j ⎜⎜ − Rs⎟⎟ + ∑ ∆Np j Ri j (B.7)
⎝βg ⎠ j ⎝ βg ⎠j ⎝βg ⎠ j j =1
⎛ βo ⎞ ⎛β ⎞ ⎛β ⎞ j +1
⎜⎜ − Rs⎟⎟ − ⎜⎜ oi − Rsi ⎟⎟ = Np j ⎜⎜ o − Rs⎟⎟ + ∑ ∆Np j +1 Ri j +1 (B.8)
⎝βg ⎠ j +1 ⎝ β g ⎠ j +1 ⎝βg ⎠ j + 1 j =1
Puesto que:
Np j +1 = Np j + ∆ jj+1 Np (B.9)
j j +1
Ri jj +1 ∆ jj+1 Np = ∑ Np j Ri j − ∑ Np j +1 Ri j +1 (B.10)
j =1 j =1
Defina:
⎛β ⎞ ⎛β ⎞ ⎛β ⎞
∆ jj+1 ⎜⎜ o − Rs⎟⎟ = ⎜⎜ o − Rs⎟⎟ − ⎜⎜ o − Rs⎟⎟ (B.11)
⎝βg ⎠ ⎝βg ⎠ j +1 ⎝ β g ⎠j
⎛ 1 ⎞ ⎛ 1 ⎞ ⎛ 1 ⎞
∆ jj+1 ⎜⎜ ⎟⎟ = ⎜⎜ ⎟⎟ − ⎜⎜ ⎟⎟ (B.12)
⎝ β g ⎠ ⎝ β g ⎠ j +1 ⎝ β g ⎠ j
Restando la Ec. (B.7) de la Ec. (B.8), reemplazando las ecuaciones anteriores y despejando
∆ jj+1 Np , se tiene;
⎛β ⎞ ⎛ 1 ⎞
(1 − Np j )∆ jj+1 ⎜ o − Rs ⎟ − β oi ∆ jj+1 ⎜ ⎟
⎜β ⎟ ⎜β ⎟
∆ jj+1 Np = ⎝ g ⎠ ⎝ g ⎠ (B.13)
⎛ βo ⎞
− Rs ⎟ + (R i ) j
j +1
⎜
⎜β ⎟
⎝ g ⎠ j +1
⎛ Ri j + Ri j +1 ⎞
( Ri )
j +1
=⎜ ⎟ (B.14)
j
⎝ 2 ⎠
B.2. PROCEDIMIENTO
Los pasos para efectuar pronósticos mediante le método de Pirson se listan a continuación:
Al inicio de la producción para presiones mayores a la presión del punto de burbuja, Pb, la
relación kg/ko = 0 y el cálculo de ∆ jj+1 Np es directo usando la Ec. (B.13) ya que ( Ri ) j
j +1
es
únicamente función de la presión. Una vez la relación de permeabilidades kg/ko tenga un
valor mayor de cero, se considera un valor de presión igual a pj+1 y se calculan los
parámetros función de presión. Todos los valores en el punto j han sido calculados en la
etapa anterior. Se asume luego un valor de ∆ jj+1 Np .
Se calcula Npj+1 utilizando la Ec. (B.9) y luego la saturación correspondiente con la Ec.
(A.20).
Con el valor de la saturación obtenido en el paso 2, se obtiene la relación de
permeabilidades kg/ko ya sea por medio de datos de laboratorio o por Correlaciónes
empíricas como se hizo en este trabajo. Luego se calcula ( Ri ) j
j +1
por medio de la Ec.
(B.14).
Se calcula ∆ jj+1 Np utilizando la Ec. (B.13). Si este valor calculado es igual al asumido
(dentro de un límite de error preestablecido) en el paso 1, se continúa con la siguiente paso.
En caso contrario, se toma el valor recientemente calculado como el nuevo valor asumido y
se repiten los cálculos hasta cuando exista convergencia. Cuide de restar el valor de
producción acumulada de la Ec. (B.9).
Se considera el siguiente valor de pj+1 y se repiten los cálculos desde el paso B.
APENDICE C
METODO DE TARNER3,6,10,13
La ecuación de balance de materiales, escrita para la forma de producción de gas, Ec. (B.1),
para el caso de un yacimiento que produce únicamente por gas en solución, es:
∆Gp ⎛ Gp 2 − Gp1 ⎞
=⎜ ⎟=
N ⎝ N ⎠
⎡ ⎛ β − β o1 ⎞ ⎛ β oi − β o 2 ⎞ ⎤ Np1 ⎛ β o1 ⎞ Np 2 ⎛ β o 2 ⎞
⎢( Rs1 − Rs2 ) + ⎜⎜ oi ⎟⎟ − ⎜⎜ ⎟⎟ ⎥ + ⎜⎜ − Rs1 ⎟⎟ − ⎜⎜ − Rs2 ⎟⎟
⎢⎣ ⎝ β g1 ⎠ ⎝ β g 2 ⎠ ⎥⎦ N ⎝ β g1 ⎠ N ⎝ β g2 ⎠
(C.1)
Las ecuaciones de saturación y relación gas-petróleo son respectivamente Ec. (A.20) y Ec.
(A.8). La Ec. (A.8), puede escribirse como:
k g µo β o
Ri = Rs + (C.2)
ko µ g β g
kg
Ri = Rs + F (C.3)
ko
µo βo
Donde F = . La producción de gas en un intervalo entre dos producciones
µg βg
acumuladas de petróleo fiscal, Np1-Np2, Ec. (C.1), también puede calcularse a partir de la
razón gas-petróleo instantánea y la producción de aceite durante el intervalo, puesto que el
gas producido resulta de:
Integrando la Ec. (C.4) entre los límites del intervalo, considerando una relación gas-
petróleo promedia para el mismo intervalo, se tiene:
⎡ Ri + Ri 2 ⎤
∆Gp = (Gp 2 − Gp1 ) = Ri( Np 2 − Np1 ) = ⎢ 1 ⎥( Np2 − Np1 ) (C.5)
⎣ 2 ⎦
6. Para cada valor obtenido en el paso anterior, se calcula la cantidad de gas producida
durante el intervalo (p1-p2) utilizando la Ec. (C.7), basada en la razón gas-petróleo
instantáneas. Por lo tanto, se obtienen tres valores de ∆Gp/N denominados A’1, A’2 y
A’3, respectivamente.
7. Los valores de ∆Gp/N del paso 2; A1, A2 y A3 se grafican contra los valores de Np2/N
asumidos. En el mismo papel se grafican los valores A’1, A’2 y A’C. En el punto de
cruce de estas curvas, se lee el valor verdadero de Np2/N correspondiente a p2; también
se lee el valor de ∆Gp/N verdadero correspondiente al intervalo en estudio. Cuando los
valores seleccionados son bien adecuados y espaciados resulta una línea recta. Al
seleccionar los valores de Np2/N debe hacerse de modo tal que el punto de corte de
líneas ocurra en un punto intermedio entre los valores de Np2/N asumidos para evitar
extrapolaciones y con ello crecimiento del error.
8. Una vez obtenido el valor correcto de Np2/N, se reemplaza en la Ec. (1.20) obteniendo
así el valor correcto de la saturación. Con este valor se obtiene el valor verdadero de la
relación de permeabilidades y luego el valor de Ri2. Los valores correctos de Np2/N y
Ri2 serán los valores de Np1/N y Ri1 para la siguiente etapa.
Gp 2 − Gp1
N ∆Gp
Ri = = (C.8)
Np 2 − Np1 ∆Np
N
APENDICE D
METODO DE TRACY3,6,10,13
• Temperatura de yacimiento.
• Presión original de yacimiento.
• Petróleo inicial in-situ.
• Saturación de agua intersticial.
• Saturación de petróleo.
• Saturación de gas.
• Factor volumétrico original de formación de petróleo.
• Viscosidad del gas a condiciones de yacimiento, para presiones entre el punto de
burbuja y la presión de abandono.
• Factores de expansibilidad promedios (cf, co y cw).
• Compresibilidad del petróleo.
1 = Npn φ on + Gp n φ gn (D.1)
Donde:
⎡ βo ⎤
⎢ − Rs ⎥
⎢ βg ⎥⎦
φo = ⎣ (D.3)
X
1
φg = (D.4)
X
⎡β ⎤ ⎡β ⎤
X = ⎢ o − Rs ⎥ − ⎢ oi − Rsi ⎥ (D.5)
⎢⎣ β g ⎥⎦ ⎢⎣ β g ⎥⎦
⎡1 − N p ⎤
S L = S w + (1 − S w ) ⎢ ⎥βo (D.6)
⎣ β oi ⎦
k g µo β o
Ri = Rs + (D.7)
ko µ g β g
cw S w cf
ce = co + + (D.8)
1 − Sw 1 − Sw
En las ecuaciones (D.1), (D.2) y (D.6), N es mantenido igual a la unidad y todos los
términos de producción son expresados como fracción del N actual. Normalmente la
solución se inicia a la presión en el punto de burbuja. Pero cualquier presión más baja en la
cual los recobros acumulativos de gas y petróleo son conocidos puede ser usada.
Usualmente se utilizan decrementos de presión de 100 a 300 psi a partir de la presión
inicial.
La ecuación (D.2) comprende dos incógnitas ∆Np y Ri. Una vez es solucionada esta
ecuación, puede hacerse lo mismo con las ecuaciones (D.6), (D.7) y (D.1). Para cualquier
presión seleccionada las ecuaciones (3) y (4) pueden ser resueltas ya que comprenden
solamente propiedades del fluido que son función de la presión.
Los factores de presión son muy sensitivos a condiciones de presión cercanas al punto de
burbuja. El cambio rápido característico de estos factores cerca a la presión de saturación y
el hecho de que la presión promedia del yacimiento no pueda determinarse de forma
precisa en estos puntos, la ecuación de balance de materia no es confiable para ser usada en
presiones cercanas al punto de burbuja. Sin embargo, después que la presión del yacimiento
disminuye a un valor tal que las pendientes de los factores con presión sean menos
pronunciadas, el petróleo original en el yacimiento puede determinarse con mejor
exactitud.
Los valores de los factores F siempre disminuyen con la presión. El factor Fn es negativo
para determinado intervalo de presión, pasa por un mínimo y luego aumenta poco a poco
hasta llegar a tener de nuevo un valor positivo a presión atmosférica.
Los datos de entrada para el método, son datos que pueden ser obtenidos en campo, tales
como presión inicial, Rsi, Temperatura del yacimiento, saturación de agua, gravedad del
petróleo, saturación de gas crítica, saturación residual de petróleo y la composición del gas.
Con estos datos es posible obtener el comportamiento de propiedades del yacimiento tales
como: Petróleo producido acumulado, Np, relación gas producido-petróleo producido, Rp,
gas producido, Gp, y propiedades de los fluidos con respecto a la presión en el punto de
burbuja durante la vida productiva del yacimiento.
D.2. CONDICIONES
• Yacimiento volumétrico.
• Solo es posible calcular el comportamiento a presiones por debajo del punto de burbuja
.
• La solución del método se inicia en la presión del punto de burbuja, pero puede
utilizarse una presión más baja a la cual loe recobros acumulativos de gas y petróleo
puedan ser usados.
APENDICE E
3,6,10,13
METODO DE SCHILTHUIS
E.1. INTRODUCCION
La ecuación de balance de materia se usa para evaluar las cantidades de fluido presentes en el
yacimiento a cualquier tiempo. De un modo especial, se usa para estimar la cantidad de
hidrocarburos iniciales en el yacimiento y predecir el comportamiento futuro y recuperación
total de los mismos.
N P ⎡⎣ βt + ( RP − Rsi ) β g ⎤⎦ − A
N= (E.1)
B − ( β ti − βt )
A = We − W p
⎡⎛ β ⎞ ⎤
B = mβti ⎢⎜ g ⎟⎟ − 1⎥ (E.2)
⎢⎣⎜⎝ β gi ⎠ ⎥⎦
N P ⎡⎣ βt + ( RP − Rsi ) β g ⎤⎦
N= (E.3)
β t − βti
En resumen, la ecuación E.3 no contempla ni empuje hidráulico ni empuje por capa de gas.
Schilthuis parte del principio que la presión inicial es la presión en el punto de burbuja. Con
base en la anterior asumpción, la cantidad de gas disuelto inicial, Rsi, es la misma cantidad de
gas disuelto, Rsb, a la presión del punto de burbuja. Luego Rsi = Rsb. El factor volumétrico
total, β t = β o + ( Rsi − Rs ) β g (visto en la unidad 3), se expresaría así:
β t = β o + ( RP − RS ) β g (E.4)
Utilizando la ecuación E.4 se calcula el factor volumétrico total inicial en el punto de burbuja:
β ti = β ob + ( Rsb − Rs ) β g
N P ⎡⎣ β t + ( RP − Rsb ) β g ⎤⎦
N= (E.5)
β t − β ob
NP
⎡⎣ β t + ( RP − Rsb ) β g ⎤⎦
1= N (E.6)
β t − β ob
Para dar solución a la ecuación E.6 se requieren dos ecuaciones adicionales a saber:
K g µo β o
Ri = Rs +
Ko µg β g
(1 − Sw ) ⎛⎜1 − ⎞
NP
⎟ βo
⎝ N ⎠
SL = Sw +
β ob
SCHILTHUIS
• Datos de propiedades de los fluidos para (ßo, ßg, Rs, µo y µg) para cada presión de
trabajo.
• Presión inicial y temperatura del yacimiento.
• Yacimiento saturado y volumétrico.
• Cantidad de petróleo inicial en el yacimiento, N, expresado en barriles a condiciones
normales.
• Saturación de agua, Sw.
Np n ⎛ Np ⎞
= ∑⎜ ⎟
N i =1
⎝ N ⎠i
4. Con el Np/N, determine la saturación de líquidos, SL, utilizando la ecuación E.8, para la
presión de interés, Pn. Con el valor de saturación de líquidos, SL, determine el valor de la
relación de permeabilidades, Kg/Ko. Para efectos de programación determine Kg/Ko
utilizando la correlación de Torcaso y Willie.
So = S L − SW
So
S* =
(1 − S w )
Kg
=
(1 − S ) (1 − ( S )
* 2 * 2
)
Ko (S ) * 4
K g µo β o
Ri = RS +
Ko µg β g
Gp n ⎛ Gp ⎞
= ∑⎜ ⎟
N i =1
⎝ N ⎠i
GP
GP N
RP = =
NP NP
N
10. Determine N − N P =
( N − NP ) = 1− NP
N N
NP
= valor obtenido en el paso 3. Despeje Np.
N
de donde:
12. Pase al siguiente valor de presión Pn+1 e inicie desde el paso 2. El ensayo termina una vez
se hallan evacuado todas las presiones de trabajo.
Presión ßo Rs Bg(x10-3) µo µg µt
INFORMACION: Los datos de propiedades de los fluidos (ßo, ßg, ßt, Rs, µo, µg) para 13
presiones diferentes (se trabajan las tres primeras presiones) de un yacimiento saturado y
volumétrico, dados en la TABLA siguiente. El yacimiento se encuentra en el punto de
burbuja a la presión inicial de 2.500 psi y a la temperatura de 180 °F.
Presión ßo Rs Bg(x10-3) µo µg ßt
7. Gp/N = 0
8. Rp = 0/0 = indeterminado
9. Indeterminado
10. 1 - 0 = 1
11. Np/N = 0, luego Np = 0
12. Para presión de 2.300 psi ir al paso 2.
Primer ensayo P = 2.300 psi.
2. Se asume ∆Np/N = 0,015
3. Np/N = 0 + 0,015
4. SL = 0,2+ 0,8*(1 - 0,015) * 1,463/1,498 = 0,969
5. So = 0,969 - 0,2 = 0,769
S* = 0,769/(1 - 0,2) = 0,961
Kg/Ko = (1-0,961)²*(1-0,961²)/0,9614 =0,00013638
Ri=669+0,00013638*(0,539/0,0166)*(1,463/0,001155)
Ri=674,6
6. ∆Gp/N = 0,015 * (721 + 674,6)/2 = 10,467
7. Gp/N = 0 + 10,467 = 10,467
8. Rp = 10,467/0,015 = 697,8
9. 0,015 [1,523+0,001155(679,8-721)]/(1,523-1,498)
= 0,897. Ensayo erróneo, asume ∆Np/N mayor.
Np/N = 0
Gp/N = 0
Ir al paso 2.