Calculo de Reservas BM Gas Seco

Reservas y
Comportamientos de
Yacimientos de Gas Seco
Ing. William Navarro
Ecuaciones de Balance de Materiales
Ley de conservacin de la masa

Masa de fludos originalmente en sitio =
Fludos producidos + Fludos remanentes en
sitio
Ecuaciones de Balance de Materiales

SUPOSICIONES DE LA EBM
El espacio poroso se encuentra inicialmente ocupado por gas y

agua connata
La composicin del gas permanece constante
Se considera Rsw = 0
Yacimiento Isotrmico
Volumen poroso homogeneo
Uniforme distribucin de la presin
No hay direccin para el flujo de fludos
Clculo de Balance de Materiales

Caractersticas
Simulan el yacimiento como un tanque de volumen constante o
variable
Es un simulador cero dimensional, OD
Da buenos resultados cuando:
No hay comunicacin de las arenas con otras
Los gradientes de presin son pequeos
Se tiene data confiable de presin y produccin
Rsw es despreciable
Clasificacin de los Yacimientos de Gas

Yacimientos Volumtricos:
Son yacimientos cerrados sin entrada de agua. El Volumen
poroso ocupado por el hidrocarburo (VPH), slo vara con el
cambio de presin y la compresibilidad Cw y Cf.
Yacimientos No Volumtricos:
Son yacimientos con empuje hidrulico. El volumen poroso
ocupado por el gas disminuye con el tiempo.
Calculo de la Presin promedia del Reservorio

Definir un plano de referencia a donde se llevaran todas las presiones para
posteriormente promediar las mismas.
Es conveniente que el plano de referencia sea el baricentro del reservorio.
Esto significa que debera dividirlo en dos regiones que tuviesen similares
volmenes de gas.
Para llevar la presin de un punto cualquiera al plano del baricentro,
podemos aplicar la ley general de la hidrosttica:
P: presin en el plano de referencia
Pm: presin en el punto de medicin
D: densidad del fluido cuya presin se esta midiendo, a condiciones de fondo.
g: aceleracin de la gravedad
Delta h: diferencia de altura, entre el punto de medicin y el plano de referencia.
Confeccionar un plano isobarico. Hay que tener presente que todos los
valores de presin del plano isobarico deben estar referidos a un mismo
plano el del referencia y deben tener la misma fecha. Si la campaa de
presiones se extendi mucho, habr que llevar todas las mediciones a un
horizonte de tiempo comn.
B.M PARA RESERVORIOS DE GAS SECO

Aplicando la ley de la conservacin de la masa a un reservorio de gas:
Moles producidos = Moles iniciales in place - Moles remanentes
n producidos niniciales nremanentes
Sabemos que :
pv
pv ZnRT n
ZRT
Aplicando la ley de los gases reales:
p Vi We W p
piVi
Z sc RTsc Z i RTr
ZRTr
psc G p
Tr: Temperatura del reservorio

Vi: volumen inicial de gas a condiciones de reservorio
Gp: Gas producido a condiciones standar, scf
Pi: presin inicial del reservorio, y
P: presin del reservorio despus de producir Gp, scf
We: entrada de agua acumulada
Wp: produccin acumulada de agua
METODO DE P/Z
El volumen de gas del reservorio puede ser puesto en unidades de pie 3
estandar, utilizando el Bg.
Vi GBgi
G, es el gas orinal in-place a condiciones standar.
Combinando las dos ecuaciones anteriores y despejando p/z tenemos que:
psc G p
Tsc Z sc
psc G pTr
Tsc GBgi
pi GBgi
Tr Z i
pGBgi
(Considerando que no existe entrada de agua)
Tr Z
pi p

Zi Z
p T
p pi
sc r G p
Z Z i Tsc Bgi G
Ecuacin conocida como el mtodo de p/z para calcular el gas in place. Cuando se grafica p/z vs
Gp se tendr una lnea recta de pendiente (Trpsc/TscBgiG) e intercepto pi/zi a Gp=0.
Variacin de P/Z con GP para un Yacimiento

Volumtrico de Gas
EJEMPLO DE APLICACIN DEL METODO P/Z

La siguiente historia de produccin fue obtenida de un reservorio volumtrico de
gas.
a) Encontrar el gas in place y el Pi/Zi.
b) Estime la recuperacin de gas cuando la presin del reservorio caiga a 300
psia.
Tr: 200F,G.Eg: 0.9.
p
psia
Gp
Bcf
p/z
1885
1620
1205
888
645
0.767
0.787
0.828
0.866
0.9
6.873
14.002
23.687
31.009
36.207
2458
2058
1455
1025
717
VISCOSIDAD DEL GAS
VISCOSIDAD Vs Tpr
Flujo Radial- Darcy
2f
pw
2e
c
s
q
p
(
h
k
3
0
7
.
0
La ley de Darcy enuncia que la velocidad aparente, de un fluido

homogeneo en un medio poroso es proporcional al gradiente de
presin e inversamente proporcional a la viscosidad del fluido.
Tz ln( re / rw)
Z = Factor de desviacin del gas, evaluado a Tres, p avg
Presin Pseudocrtica
647.9
Temperatura Pseudocrtica
439.9
Psr
Tsr
(3000+1500/2)/647.9= 3.5
1.5
0.75
T = Temperatura, R
qsc = Caudal de gas a 14.7 psia, 60 F, scf/day
UgA @ 1atm y 200F

U/Uga
U
g = Viscosity del gas a T, pavg, cp

kg = permeabilidad al gas, md
h = espesor neto, pies
re = radio de drenaje, pies
De la Figura
rw = radio del pozo, pies
Diametro de la broca
Pe = Presin del reservorio, psi

Pwf= Presin fluyente, psi
0.0116
1.9
0.02204
VISCOSIDAD DEL GAS

USANDO EL METODO DE LEE
Lee, presento una expresin analtica para la viscosidad de hidrocarburos
Gaseosos en 1966.
4
y
K
10
exp(
X
La ecuacin es:
g
g)
Donde:
9.4 0.02M T 1.5

K
209 19M T
986
X 3.5
0.01M
T
y 2 .4 0 .2 X
En estas ecuaciones:
T en R
g en centipoise
M es el peso molecular
g en gr/cm3
VISCOSIDAD - Z DEL GAS

P g
PM 28.97 P g
g
2.7
ZRT 10.73ZT
ZT
g 0.04327
P g
ZT
Lb/ft3
gr/cc
Beggs& Brill
(1 A)
D
CP
pr
eB
A 1.39(Tpr 0.92) 0.5 0.36Tpr 0.101
Z A
0.066
2
0.32
B (0.62 0.23Tpr ) Ppr
0.037 Ppr 9
Ppr6
10 (Tpr 1)
(Tpr 0.86)
C (0.132 0.32 log Tpr )

D 10
2
( 0.3106 0.49T pr 0.1824T pr
)
Z
(1 A)
D
CP
pr
eB
A 1.39(Tpr 0.92) 0.5 0.36Tpr 0.101
Z A
0.066
2
0.32
B (0.62 0.23Tpr ) Ppr
0.037 Ppr 9
Ppr6
10 (Tpr 1)
(Tpr 0.86)
C (0.132 0.32 log Tpr )

D 10
2
( 0.3106 0.49T pr 0.1824T pr
)
BM General
Sabemos que:
Vi GBgi
Reemplazando esta relacin en la ecuacin, se tiene
p Vi We W p Bw
piVi
Z sc RTsc Z i RTr
ZRTr
psc G p
piVi p Vi We W p Bw
Z scTsc Z iTr
ZTr
Z scTsc piVi p Vi We W p Bw
Gp
psc Z iTr
ZTr
psc G p
Z scTsc
Gp
psc
pi Gi Bgi p Gi Bgi We W p Bw
Z
T
ZT
r
i r
1/Bgi
Ecuacin General de Materia

Gi Bgi Gi Bgi We W p Bw
Gp
Bg
Bgi
Gi Bgi We W p Bw
G p Gi
Bg
Bg
Gp
GB
Gi Bgi
Bg
W p Bw
Bg
G p Bg Gi Bg Bgi We W p Bw
G p Bg W p Bw Gi Bg Bgi We
Ecuacin de Balance de Materiales

Ecuaciones Generales
GP
G
Bgi
Bg
Bgi CwSwi Cf
Bg
1Swi
We WpBw
GBg
Recobro por intrusin de agua

Recobro por expansin de agua connata y
reduccin de porosidad
Recobro por expansin de Gas
Recobro total de Gas
METODO DE P/Z
Sabemos que:
Bgi
Psc Z iTr
Pi Z scTsc
Combinando las dos ecuaciones anteriores y despejando p/z tenemos que:
p pi

Z Zi
pscTr
psc Z iTr
piTsc
Tsc
p pi pi G p

Z Zi Zi G
Gp
p pi G p
1
Z Zi
G
La ecuacin puede ser expresada en trminos del factor de recuperacin.
METODO DE SAMANIEGO
La ecuacin la podemos escribir como:
Gp
pZ i
p pi G p

1
1
Z Zi
G
pi Z
G
Gp
pZ
1 i
G
pi Z
Tomando logaritmos
Gp
pZ i
log
log1
pi Z
G
pZ i
log G p log G
log1
pi Z
Como podemos ver la ecuacin anterior es la ecuacin de una lnea recta de pendiente 1, es decir
tiene un ngulo de 45 cuando se grafica en la escala log-log. La ecuacin anterior se cumple
cuando We=Wp=0
METODO DE SAMANIEGO
Extrapolando la recta para p=0 se obtiene el gas original in-place.
We=0
We0
pZ
log1 i
pi Z
) 45
Com. alta y variable
Log Gp
Gi
Entrada de agua: el mantenimiento de la presin hace que la funcin:
pZ i
log G p log G
log1
pi Z
Se ubique debajo de la recta de 45 grados.

Compresibilidad: en aquellas formaciones con alta compresibilidad ( arenas no
consolidadas), ante las primeras extracciones de fluidos cambian su
empaquetamiento, disminuyendo el volumen poral, el efecto es un
mantenimiento de presin inicial, y al disminuir la compresibilidad tiende a
buscar la linearidad.
Este grafico tambin se conoce como plots de energa. Son muy tiles en
detectar entrada de agua a inicios de la vida de un reservorio.
METODO DE HAVLENA Y ODEH

Havlena y Odeh parten de la siguiente ecuacin:
C w S wi C p
Gi B g B gi We Gi
pi p G p B g W p
1 S wi
Expansion del gas

(Eg)
Entrada de agua
Expansin de poros y
agua intersticial
(Epw)
La ecuacin anterior la escriben en la forma:
Gi E g We Gi E pw F
Siendo:
F G p Bg W p
E g Bg Bgi
E pw
Cw S wi C p
pi p
1 S wi
Fluidos
producidos en
condiciones de
fondo (F).
METODO DE HAVLENA Y ODEH

Despreciando Epw:
F Gi E g We
Acuifero
potente
F
Eg
Acuifero
debil
F es funcin de la produccin y presin

Eg es funcin de la presin.
Al graficar F/Eg vs Acumulado de gas:
- Si no hay entrada de agua se obtiene
una horizontal
- Si hay acuifero hay una desviacin
kte
We=0
Este mtodo es mas sensible para

analizar la entrada de agua.
Gas original in-place
Gp
EJEMPLO N2 DETERMINACION DE LA ENTRADA DE AGUA

Problema:
Para el reservorio cuyos datos figuran a continuacin, determinar si ha habido entrada
de agua utilizando las siguientes tcnicas:
A) P/Z, b) Havlena y Odeh, c) Samaniego
Nota:
El GOIS volumtrico dio 823 x109 scf.
La temperatura de fondo es de 210 F
METODO DE P/Z vs GAS ACUMULADO
GOIS:
1280 BCF
METODO DE HAVELENA Y ODHE
METODO DE HAVELENA Y ODHE
GOIS:
823 BCF
GOIS:
1280 BCF
GOIS:
823 BCF
METODO DE SAMANIEGO
METODO DE SAMANIEGO
Respuestas tpicas
P/Z
P/Z
Gp
P/Z vs Gp para un yacimiento de Gas

con Empuje Hidrulico
P/Z
Gp
con Presin Anormal
P/Z
Gp

con Condensacin Retrgrada
Gp
P/Z vs Gp para un yacimiento
con variacin de K
ndice de Empuje
CwSwi Cf P We WpBw
GpBg GBg Bgi 1GBgi
Swi
G BgBgi
GpBg
IE de Expansin
de Gas
G BgiCwSwi Cf P
1Swi GpBg
IE de Expansin de
Agua Connata y
reduccin del Vp
We WpBw
GpBg
IE de Intrusin de
Agua
Validacin del OGIP Volumtrico por BM
F
We
G
Eg Efw
Eg Efw
We Pequeo
We correcto
F
Eg Efw
We grande
OOGIP
We
Eg Efw
Yacimientos de Gas con Presiones Anormales
Wallace:
Son aquellos que tienen una presin que excede a la calculada
con un gradiente de 0.456 lpc/pie
Caractersticas:
Formados por secuencias de arena / arcilla
Cf Alta 28 x 10-6 Lpc-1
Cf Normal 6 x 10-6 Lpc-1
Presin de Relajacin:
Presin a la cual el yacimiento geopresurizado alcanza
compresibilidad normal
Reservas de Yacimientos Geopresurizados

Ecuacin a resolver:
P
Pi PiGp
Cpa
Z
Zi ZiG
CwSwi Cf Pi P
Cpa 1
1 Swi
Cpa = Correccin por Presin Anormal
Ecuacin a resolver:
Graficar P Cpa vs Gp
Z
Extrapolar la lnea recta para hallar el GOES a
P
Cpa 0
Z
De la recta hallar las reservas, Gpa a unas condiciones de abandono

P
Cpa
Z
Reduccin del volumen poroso por compactacin del

yacimiento N S2 B. Campo North Ossun, LA.
Determinacin del OGIP de un Yacimiento de Gas con

Presiones Anormal
Yacimientos de Gas con Empuje Hidrulico
Estos yacimientos adems de la energa que tienen como

resultado de la compresibilidad del gas, tienen una entrada
adicional de energa con el agua de intrusin que mantiene la
presin de los mismos.
Determinacin de We:
Mtodo de BM
Ecuaciones de Van Everdingen y Hurts
Acufero Pote
Mtodo de BM para calcular We

Ecuacin:
We GpBg GBg Bgi WpBw

Se dispone del:
Histrico de Presin
Histrico de Produccin de gas y agua
GOES
La Ecuacin se aplica a varias presiones durante el agotamiento
Mtodo de clculo de Intrusin de Agua, We

A) Van Everdingen Hurts
Fetkovitch
Acuferos grandes activos

Influjo de agua transitorio
B) Schilthuis (Acufero Pote)
Acuferos pequeos
rD 2
Influo de agua contnuo

Influjo Transitorio:
Tasa de influjo de agua menor que la

tasa de drenaje de gas
Influjo Continuo:
Tasa de influjo de agua igual a la tasa de

drenaje de gas
Ecuaciones de Van Everdingen - Hurts

Para una sola cada de presin DP durante un tiempo t
We C * P * Q td
2
C 1.119 * * Ce * rR * h * f
Q td f t d * rD
t d 6.323 10
rA
rD
rR
Ce Cw Cf
K*t
* w * Ce * rR
Ejemplos ilustrativos de estimacin de f y rD
Principio de Superposicin en el Tiempo

El influjo total de agua durante un tiempo tn de explotacin de
yacimiento, es igual a la suma de los influjos parciales de agua
producidos por DPo, , DPj, , DPn-1 cadas de presin.
We = C DPo Qtdn (We producido por DP0 durante tn)
+ C DP1 Q(tdn-td1) (We producido por DP1 durante tn t1)
+ C DPj Q(tdn-tdj) (We producido por DPj durante tn tj)
+ C DPn-1 Q(tdn-tdn-1) (We producido por DPn-1 durante tn tn-1)
Po Pi P1
Pi P1 Pi P1
P
P1 Pi P2
Pi P1 P1 P2 Pi P2
Pj Pj Pj1
2
Pj1 Pj
2
2
2
Pj Pj1 Pj1 Pj1
2
2
Ajuste de la distribucin continua de Geometra del Yacimiento de Gas con

Empuje Hidrulico
presin en el contacto Agua-Gas a
travs de una Distribucin
Escalonada de Presin
Influjo de Agua durante un tiempo Tn

Yacimientos Radiales
Agua
n 1
We C * Pj * Q t dn t dj
Pj
j 1
j0
Pj1
2
t dj 6.323 10 3
Yacimientos Lineales
Gas
Kt j
w CerR
n 1
KCe
We 3.26 10 3 hw
P j t n t j
w j0
Agua
Gas
Determinacin de Qtd
A) Tablas de Van Everdingen y Hurts (trans. AIME 1949,
Libros de Craft y Hawkins, IKOKU)
B) Ajuste Matemtico (Bird y Cols, Exxon)
Para rD < 100
1.6179
Q td Y1
1 Y2
Y1 1.07054 t d
td
Y2
1
0.5
1.21257
Y2
3
1 Y2
0.500712t d
1 0.532526rD 1
2.33849
2 2.72055 0.00401826rD
3 rD 1 0.5
0.90484
EBM como una Lnea Recta Mtodo de Havlena y

Odeh
Se resuelve simultneamente la EBM y
y ax b
la Ecuacin de Van Everdingen y Hurts
WpBw
y GpBg
BgBgi
ac
para determinar GOES y We:
n 1
PjQ tdn tdj

j0
bG
BgBgi
Para un Yacimiento de Gas con Empuje Hidrulico, y vs x es una

lnea recta
Obtener:
Del Intercepto
GOES
De la Pendiente
Procedimiento para hallar GOES y We
1. Suponer la geometra y tamao del acufero.

2. Determinar Gp, Bg, Bgi, Wp, Bw y S DPQtd a diferentes
presiones
3. Graficar Y vs X
4. La variacin de Y vs X es lineal?
SI:
Determinar GOES y B
NO:
El acufero no est bien caracterizado y se
debe variar el tamao y/o la geometra del
mismo
P Q
Comportamiento
Lineal
td
muy pequeo
Cambiar la
geometra del
acufero
P Q
td
muy grande
EBM como una Lnea Recta para Yacimientos de Gas

con Empuje Hidrulico
Modelo de Acufero POTE

Determinacin Simultnea de G y W
La EBM para Yacimientos de Gas puede escribirse en la forma en
la forma siguiente:
F
Eg
con
G PiEg P GBgi
SwiCw Cf
1 Swi
We W Cw Cf Pi P
Cw Cf W
vs PiEgP
Lnea Recta:
F
Eg
Intercepto:
bG
Pendiente:
Cf
Cw Cf W
m GBgi SwiCw
1Swi
Cf GBgi
m SwiCw
1Swi
W
Cw Cf
Determinacin de G y W.
Modelo del Acufero POTE

FACTOR DE RECOBRO
FORMA CONCEPTUAL DE UN YACIMIENTO DE GAS INVADIDO POR
AGUA

FACTOR DE RECOBRO
43560.Vb. .Sgr.Ev 43560.Vb. .Sgi 1 Ev

Gp G
, PCN
Bg
Bg
Pi P
1 Ev
Gp Zi Z
Sgi.Pi / Zi Sgr.P / Z Ev
, fraccin
Pi
G
Sgi.Pi / Zi
Zi
Recobro por expansin
del gas de la zona no
invadida por agua
% Recobro = (Gp / G) * 100
Recobro de gas por

desplazamiento de agua y
por expansin del gas en la
zona invadida

SATURACIN RESIDUAL DE GAS EN ZONAS INVADIDAS
POR AGUA
Debido a que el desplazamiento de gas por agua es un proceso de
imbibicin, Sgr es alta.
Pruebas de Laboratorio : Sgr = 15 - 50 %
Correlaciones:
- Geffen y Cols: Sgr = 35 % (promedio)
- Naar y Henderson: Sgr = Sgi / 2
- Legatasky y Cols: Sgr = 0.625 1.3125*Porosidad
- Land: (1 / Sgr) = (1 / Sgi) + [C / (1 Swi)]
C = 1 para arena no consolidada y C = 4 para arena consolidada

FACTORES QUE AFECTAN EL RECOBRO DE GAS
POR DESPLAZAMIENTO CON AGUA
TASA DE PRODUCCIN DE GAS
Se obtienen mayores recobros a altas tasas de produccin de gas
que permitan agotar la presin del yacimiento antes que los pozos
sean invadidos por agua.
SATURACIN RESIDUAL DE GAS
A medida que disminuye Sgr aumenta el recobro de gas.
A mayor consolidacin y contenido de arcilla de las areniscas,
mayor Sgr.
EFICIENCIA VOLUMTRICA DE BARRIDO
A mayor Ev, mayor es el recobro de gas.
Entre ms heterogneo sea el yacimiento, menor la Ev.

Calculo de Reservas BM Gas Seco

Cargado por

Copyright:

Formatos disponibles

Calculo de Reservas BM Gas Seco

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Calculo de Reservas BM Gas Seco

Cargado por

Copyright:

Formatos disponibles

Reservas y

Ing. William Navarro

Ecuaciones de Balance de Materiales

Ley de conservacin de la masa

Ecuaciones de Balance de Materiales

El espacio poroso se encuentra inicialmente ocupado por gas y

Clculo de Balance de Materiales

Clasificacin de los Yacimientos de Gas

Calculo de la Presin promedia del Reservorio

B.M PARA RESERVORIOS DE GAS SECO

n producidos niniciales nremanentes

Aplicando la ley de los gases reales:

Tr: Temperatura del reservorio

G, es el gas orinal in-place a condiciones standar.

Combinando las dos ecuaciones anteriores y despejando p/z tenemos que:

(Considerando que no existe entrada de agua)

Variacin de P/Z con GP para un Yacimiento

EJEMPLO DE APLICACIN DEL METODO P/Z

VISCOSIDAD DEL GAS

Flujo Radial- Darcy

La ley de Darcy enuncia que la velocidad aparente, de un fluido

Z = Factor de desviacin del gas, evaluado a Tres, p avg

UgA @ 1atm y 200F

g = Viscosity del gas a T, pavg, cp

rw = radio del pozo, pies

Pe = Presin del reservorio, psi

VISCOSIDAD DEL GAS

9.4 0.02M T 1.5

VISCOSIDAD - Z DEL GAS

C (0.132 0.32 log Tpr )

C (0.132 0.32 log Tpr )

Reemplazando esta relacin en la ecuacin, se tiene

Ecuacin General de Materia

Ecuacin de Balance de Materiales

Recobro por intrusin de agua

Combinando las dos ecuaciones anteriores y despejando p/z tenemos que:

Com. alta y variable

Entrada de agua: el mantenimiento de la presin hace que la funcin:

Se ubique debajo de la recta de 45 grados.

METODO DE HAVLENA Y ODEH

Expansion del gas

La ecuacin anterior la escriben en la forma:

METODO DE HAVLENA Y ODEH

F es funcin de la produccin y presin

Este mtodo es mas sensible para

Gas original in-place

EJEMPLO N2 DETERMINACION DE LA ENTRADA DE AGUA

METODO DE P/Z vs GAS ACUMULADO

METODO DE P/Z vs GAS ACUMULADO

METODO DE HAVELENA Y ODHE

METODO DE HAVELENA Y ODHE

METODO DE P/Z vs GAS ACUMULADO

P/Z vs Gp para un yacimiento de Gas

P/Z vs Gp para un yacimiento de Gas

Validacin del OGIP Volumtrico por BM

Yacimientos de Gas con Presiones Anormales

Reservas de Yacimientos Geopresurizados

Extrapolar la lnea recta para hallar el GOES a

De la recta hallar las reservas, Gpa a unas condiciones de abandono