PSEUDO PRESION
ECUACION DE DIFUSIVIDAD DE GASES
La produccin de gas es ms complejo que el petrleo, por el hecho que las propiedades del
gas son ms fuertes funciones de la presin. De aqu de la ecuacin de difusividad desarrollada
para petrleo no se puede lineal izar.
Los problemas principales, son:
- La alta dependencia de la viscosidad del gas con respeto de la presin
- La compresibilidad del gas
Por ende, para poder resolver este problema se introdujo el concepto de pseudo presin m(p)
o potencial real del gas, se define como:
Como se ve, esta pseuda presin considera la viscosidad y la compresibilidad z funciones de
presin. Donde Po es la presin arbitraria elegida inferior a las pseudo presiones durante el
ensayo.
As, la ecuacin de difusividad nos quedara como:
Para realizar esta ecuacin, en la prctica se consideran a estos parmetros independientes de
P, y se los calcula a P=Pi.
Considerando flujo estabilizado:
Anlisis de la ecuacin:
Para calcular m(p), es necesario evaluar su integral. Una forma de realizar esto, es aplicando la
regla de Simpson o del trapezoide:
Los datos de viscosidad y compresibilidad salen de ensayos PVT de laboratorio o de
correlaciones.
Ejemplo:
Otra forma de calcular m(p), es poder establecer una relacion con lapresion para cada pozo. Si
nosotros analizamos el grafico de m(p) vs p, vemos se verifica que a bajas presiones, m(p), se
ajusta a una expresion cuadratica y a altas presiones a una relacion lineal.
En la mayoria de los casos, se considera la relacion cuadratica, no por que sean mas exactas
sino porque esta ecuacion es de simplecomparacion con los metodos antiguos de analisis de
ensayos para gas.
Entonces generalizando el modelo para cualquier presion p del reservorio, que puede ser de
mucho menor que la presion inicial despues de aos de produccion:
CALCULO DE LA FUNCION
Este ejemplo se ilustra el clculo de la funcin por una simple integracin numrica.
S.G= 0.61 TPC= 357R
TY= 120F I = 0.0114 CP
PPC= 664 psia TPr= 1.625 R
COLUMNA TERMINOS UNIDADES FUENTE
D 1 P Psia Dato
A
T 2 Z Dato
O
S 3 Cp Dato
I N 4 Mpsia/cp 2*1/2/3
2*
N U
T M 5 Promedio 2* Mpsia/cp Promedio de las dos entradas
E E de la colum. 4
G R 6 P Psia Diferencia de dos sucessivas
R I entradas de 1
A C 7 MM 5*6/100
Promedio(2* )
C A psia2/cp
I 8
O ( ) MM 7+entrada previa de 8
N psia2/cp
9 Pr =1/Pc
10 Valor de la integral De la tabla 2-2
11 MM 2*Pc2(10)/i
( )
psia2/cp
Datos de la integracin numrica usando la table 2-
27
P Z (5)*(6)/100 *106 Pr Val
Integr
1 2 3 4 5 6 7 8 9 10 11
0.0 0.0 0.0 0.0 0.0 0.kk 0.0 0.0 0.0 0.0 0
400 0.955 0.0118 70.99 35.50 400 14.2 14.2 0.6 0.1824 14.1
800 0.914 0.0125 140.0 105.5 400 42.2 56.4 1.2 0.7253 56.1
1200 0.879 0.0134 203.8 171.9 400 68.8 125.2 1.8 1.6102 124.5
1600 0.853 0.0145 258.7 231.2 400 92.5 217.2 2.4 2.8054 217.0
2000 0.838 0.0156 306.0 282.4 400 112.9 330.7 3.0 4.2605 329.6
FUNCION DEL POTENCIAL DEL GAS REAL
EJM: Con los siguientes datos PVT de un pozo de gas del campo provenir se tiene a
continuacin:
P(psia) g(cp) Z 2P/ (psia2/cp)x106
g*Z(psia/cp)
0 0.01270 1.000 0 0
400 0.01286 0.937 66391.0 13.2
800 0.01390 0.882 130507.8 52.0
1200 0.01530 0.832 188537.0 113.1
1600 0.01680 0.794 239824.4 198.0
2000 0.01840 0.770 282326.4 304.0
2400 0.02010 0.763 312982.9 422.0
2800 0.02170 0.775 332986.5 542.4
3200 0.02340 0.797 343167.2 678.0
3600 0.02500 0.827 348246.7 816.0
4000 0.02600 0.860 357781.8 950.0
4400 0.02831 0.896 346924.4 1089.0
El pozo esta produciendo con una presin estabilizada de fondo de 3600 psi. Tiene un radio de
0.3 ft.
Los datos adicionales son las siguientes:
K= 65md h=15ft T= 600R
Pe= 4400psi re=1000ft Pwf= 3600psi
Calcular el caudal delgas natural en Mscf/dia.
Paso 1.- Calcular la funcin m(p) R=m(PR)=
Paso 2.- Graficar el termino () versus la presin
Paso 3.- Calcular numricamente el rea bajo la curva para algunos valores de p.
Estas reas corresponden al potencial del gas real para con una presin determinada.
Paso 4.- Calcular el caudal del gas aplicando la correspondiente ecuacin.
Pw= 816.0x106 Pe=1089x106
()()()
Qg= = , / wf=m(Pwf)=
()()( )
.
Presin *106
0 0 0
400 66931.0 13.20
800 130507.0 52.00
1200 188537.0 113.10
1600 239834.0 198.00
2000 282326.0 304.00
2400 312982.1 422.00
2800 332986.0 542.40
3200 343147.2 678.00
3600 348246.0 816.00
4000 357781.0 950.00
4400 348824.4 1089.00
ANALISI DE FLUJO LIT
La ecuacin tiene en cuenta los efectos de flujo laminar, inercial y turbulento y tiene en cuenta
las siguientes suposiciones:
Condicin isotrmica
Efectos gravitacionales insignificantes
Flujo monofsico
Medio homogneo e isotrpico de porosidad constante
Permeabilidad independiente de la presin
Viscosidad y compresibilidad constantes
Se aplica el modelo de flujo radial
Debido a estas limitaciones se utiliza la aproximacin con pseudo presiones
La pseudo presin est definido por:
Donde po es la presin de referencia: en cualquier problema que implique la utilizacin de
es recomendable una tabla de conversin o un grfico de conversin. Tambin de puede
realizar por integracin numrica.
El grafico de anlisis es el siguiente:
