Análisis de Pozos de Gas y Presión
Análisis de Pozos de Gas y Presión
Un reservorio de gas seco cuenta con dos pozos de producción donde el pozo x1 de tipo (J) cuenta
con un radio de 1385 ft ya que se registra un desplazamiento horizontal de 1226 ft en la
profundidad donde el Angulo comienza a mantenerse constante el cual es 6310 ft. La trayectoria
del pozo x2 cuenta con un Kop de 2845 pies seguida de una sección donde se genera el ángulo y
se mantiene constante hasta una profundidad de 7283 pies en el cual comienza a disminuir a un
ritmo de 2.3º por cada cien pies hasta generar una sección vertical de 1186 pies con el cual se
obtiene un desplazamiento horizontal hasta el objetivo de 4183 pies ya que toda la longitud
medida desde superficie hasta el objetivo es de 11577 pies.
El pozo x2 produce gas con un caudal de 10.5 MMPCD con una gravedad especifica del gas de
0.66 con un API de 50º ya que está conformado con una tubería de producción de 2.44” donde
se tiene una rugosidad de 0.0022” Y también se registra una presión de cabeza de 2320 psi a una
temperatura de 68ºF en cabeza y 237,3 ºF en el fondo del pozo y un índice de productividad de
1,83PC/Psi^2
El pozo x1 producía a un caudal de 11.5 MMPCD con una relación de gas condensado de 89876PCS
por barril a una temperatura de 68 ºF en cabeza y un gradiente geotérmico de 16,8ºF/1000 ft ya
que se llegó a determinar un índice de productividad de 2,175 PC/Psi^2 donde más después se
procedió a cerrar el pozo con una presión de 1345 psi en cabeza donde la gravedad especifica del
gas es de 0.65 con un API de 48º en el cual el diámetro es de 3” con una rugosidad de 0.0028”
𝜃 𝐷3 − 𝐷2 𝜃
𝑀𝐷𝑋2 = 𝐷1 + + + + 1186 → 𝑫𝟐 = 𝑫𝟏 + 𝒓𝟏 ∗ 𝐬𝐢𝐧 𝜽 = 𝟐𝟖𝟒𝟓 + 𝒓𝟏 ∗ 𝐬𝐢𝐧 𝜽
𝑞1 cos 𝜃 𝑞2
𝜽 𝟕𝟐𝟖𝟑 − 𝑫𝟐 𝜽
𝟏𝟏𝟓𝟕𝟕 = 𝟐𝟖𝟒𝟓 + + + + 𝟏𝟏𝟖𝟔
𝒒𝟏 𝐜𝐨𝐬 𝜽 𝟐. 𝟑
𝟏𝟎𝟎 𝟏𝟎𝟎
Sistema de ecuaciones 4 ecuaciones con 4 incógnitas
𝟓, 𝟓𝟐𝟒º
𝒓𝟏 = 𝟏𝟎𝟑𝟕 𝒇𝒕 𝑫𝟐 = 𝟑𝟓𝟑𝟎, 𝟐𝟏 𝒇𝒕 𝒒𝟏 = 𝜽 = 𝟒𝟏, 𝟑𝟓º
𝟏𝟎𝟎𝒇𝒕
𝐷4 = 𝐷3 + 𝑟2 ∗ sin 𝜃 𝑑𝑜𝑛𝑑𝑒 𝑻𝑽𝑫𝒙𝟐 = 𝑫𝟒 + 𝟏𝟏𝟖𝟔𝟎 = 𝟏𝟎𝟏𝟏𝟒, 𝟖 = 𝟏𝟎𝟏𝟏𝟓 𝒇𝒕
Calculamos la presión de fondo con el método de presión y temperatura promedio del pozo X2 en estado
monofasico
Dens Pwh
Pwf* Z (gr/cm3) u (cp) Nre f S calculado
0,17648201 2701588, 0,01930492 0,489916 3447,42340
4000 0,83452779 5 0,02102588 6 6 2 4
3447,42340 0,82653074 0,16260990 0,01992229 2851241, 0,01929676 0,494656 3450,12263
4 8 6 7 4 8 3 3
3450,12263 0,82655983 0,16268028 0,01992767 2850472, 0,01929680 0,494638 3450,11277
3 1 6 1 5 8 9 9
3450,11277 0,82655972 0,16268002 0,01992765 2850475, 0,01929680 0,494639 3450,11281
9 5 9 1 3 8 0 5
3450,11281 0,82655972 0,01992765 2850475, 0,01929680 0,494639 3450,11281
5 5 0,16268003 1 3 8 0 5
Pwf x2=3450 psi
b)
Aplicamos cullender estático para el pozo X1 y calcular la presión.
16,8º𝐹
𝑇𝑤𝑓𝑥1 𝑎 7856 𝑓𝑡 = 68 + ( )∗ 7856 = 199,98º𝐹 = 200º𝐹
1000 𝑓𝑡
Pm* Th Z Α Im Pm calcu
2000 593,67 0,816986785 95,745 0,242510272 1519,569622
1519,569622 593,67 0,843225372 95,745 0,329433808 1495,687772
1495,687772 593,67 0,844969913 95,745 0,335386367 1494,289171
1494,289171 593,67 0,84507331 95,745 0,335741356 1494,206583
1494,206583 593,67 0,84507942 95,745 0,33576234 1494,201704
Pm=1494, 2 psi
C)
Calcular la longitud medida del pozo para la profundidad de 8567 ft ya que sabemos que esa profundidad
D3<8567 ft<D4
𝜃 𝐷3 − 𝐷2 ∝ 8567 − 𝐷3
𝑀𝐷𝑋2 𝑎 8567 𝑓𝑡 = 𝐷1 + + + 𝑑𝑜𝑛𝑑𝑒 𝑠𝑒𝑛 ∝=
𝑞1 cos 𝜃 𝑞2 𝑟2
8567 − 𝐷3 8567 − 7283
∝= sin−1 ( ) = sin−1 ( ) = 31,026º
𝑟2 2491,12
40,35 7283−3530,21 31,026
𝑀𝐷𝑋2 𝑎 8567 𝑓𝑡 = 2845 + + + 2,3 =
5,524/100𝑓𝑡 cos 40,35 𝑓𝑡
100
9848,7 𝑓𝑡
El pozo x1 produce un caudal de 10.5 MMPCD con una relación de gas condensado de 169876PCS
por barril donde se registra una presión de 2345 psi a la profundidad del Kop con una
temperatura de 68ºF en cabeza y 295ºF de fondo donde la gravedad especifica del gas es de 0.66
con un API de 50º en el cual el diámetro de la tubería de producción de 3” con una rugosidad de
0.0012” ya que se registra un índice de productividad de 0,824 PC/psi^2
Datos que se tienen del pozo X1 son (q1=2.85º/100 ft) (x2=325ft) (D2=5842) (MD=11567ft)
180 1 180 1
𝑟1 = ∗ = ∗ = 2010 𝑓𝑡
𝜋 𝑞1 𝜋 2.85
100 𝑓𝑡
𝑟−𝑋2 2010−325
𝑋2 = 𝑟1 ∗ (1 − cos 𝜃) 𝑑𝑒𝑙 𝑐𝑢𝑎𝑙 𝑠𝑒 𝑜𝑏𝑡𝑖𝑒𝑛𝑒 𝜃 = cos −1 ( ) = cos −1 ( ) = 33º
𝑟 2010
Datos que se tienen del pozo X2 son (D1=2530ft) (q1=1.75º/100 ft) (D3=7283ft) (q2=2.3ª/100ft)
(x4=3428ft)
180 1
𝑟1 = ∗ = 3274 𝑓𝑡 𝑑𝑜𝑛𝑑𝑒 𝑋2 = 𝑟1 ∗ (1 − cos 𝜃) = 4092.55 ∗ (1 − cos 𝜃)
𝜋 1.75º
100 𝑓𝑡
180 1
𝑟2 = ∗ = 2491 𝑓𝑡 𝑑𝑜𝑛𝑑𝑒 𝑋3 = 𝑋2 + tan 𝜃 ∗ (𝐷3 − 𝐷2)
𝜋 2.3º
100 𝑓𝑡
𝟐𝟒𝟖𝟐 = 𝟑𝟐𝟕𝟒 ∗ (𝟏 − 𝐜𝐨𝐬 𝜽) + 𝐭𝐚𝐧 𝜽 ∗ (𝟕𝟐𝟖𝟑 − (𝟐𝟓𝟑𝟎 + 𝟑𝟐𝟕𝟒 ∗ 𝐬𝐢𝐧 𝜽)) + 𝟐𝟒𝟗𝟏 ∗ (𝟏 − 𝐜𝐨𝐬 𝜽)
𝑇𝑉𝐷𝑋2 = 𝟏𝟎𝟑𝟓𝟓 𝒇𝒕
𝜃 𝐷3 − 𝐷2 𝜃
𝑀𝑉𝐷𝑋2 = 𝐷1 + + + + 1864 → 𝐷2 = 𝐷1 + 𝑟1 ∗ sin 𝜃 = 2530 + 3274 ∗ sin 29 = 4117 𝑓𝑡
𝑞1 cos 𝜃 𝑞2
𝑴𝑽𝑫𝑿𝟐 = 𝟏𝟎𝟗𝟑𝟐 𝒇𝒕
𝜃 𝐷3−𝐷2
Calculamos el MD desde el Kop hasta el fondo para el pozo x1 𝑀𝐷 = + = 6820 𝑓𝑡
𝑞 𝑐𝑜𝑠𝜃
De la misma forma ya que se tiene la Twh, Twf del pozo x1 con el cual calculamos el gradiente geotérmico
mediante el TVD del pozo X1 con la ecuación de 𝑇𝑤𝑓 = 𝑇𝑤ℎ + 𝐺º ∗ 𝑇𝑉𝐷
𝑇𝑤𝑓 − 𝑇𝑤ℎ 295 − 68
º𝐺 = = = 0,021433º𝐹/𝑓𝑡
𝑇𝑉𝐷𝑥1 10591 𝑓𝑡
𝑇𝑘𝑜𝑝 = 𝑇𝑤ℎ + 𝐺º ∗ 𝐷1 = 169.74º𝐹 = 629,4º𝑅
a) Para el análisis de caudales críticos comenzamos calculando la presión de cabeza del pozo x1 con el
método de cullender en el cual se considerará el Sg de mezcla y estado bifásico para el cálculo de Z
4591 ∗ 𝛾
𝑆𝑔 + 𝑅𝐺𝐶 𝑜 141.5
𝑆𝑔𝑚 = 𝑑𝑜𝑛𝑑𝑒 𝛾𝑜 =
1123 131.5 + 𝐴𝑃𝐼
1 + 𝑅𝐺𝐶
4591 ∗ 𝛾𝑜
0,66 + 141.5
𝑆𝑔𝑚 = 169876 = 0,6766 𝑑𝑜𝑛𝑑𝑒 𝛾𝑜 = = 0,77961
1123 131.5 + 50
1 + 169876
Para los caudales críticos calcular la presión de cabeza partiendo del Kop como presión de
fondo y profundidad del Kop en estado monofásico ya que es un reservorio de gas seco.
Nº
iteraciones Pm* Z Α F Im Pm calcu
1 2000 0,78392035 60221,62875 0,064300072 187,0419985 2183,227896
2 2183,227896 0,778460993 60221,62875 0,064300072 175,4590096 2178,032685
3 2178,032685 0,77857371 60221,62875 0,064300072 175,7696687 2178,176373
4 2178,176373 0,77857056 60221,62875 0,064300072 175,7610624 2178,172395
5 2178,172395 0,778570648 60221,62875 0,064300072 175,7613006 2178,172505
Nº
iteraciones Pwh* Z α F Iwh Pwh calcu
1 3300 0,750841516 60221,62875 0,064300072 113,3698129 1969,887552
2 1969,887552 0,705861281 60221,62875 0,064300072 164,8055778 2001,344781
3 2001,344781 0,704367256 60221,62875 0,064300072 162,611517 2000,198203
4 2000,198203 0,704419418 60221,62875 0,064300072 162,6904567 2000,239714
5 2000,239714 0,704417526 60221,62875 0,064300072 162,6875975 2000,23821
Pwh del pozo X1=2000 Psi con el cual Realizamos el análisis de caudales criticos
ρo = 62.4 ∗ 𝛾𝑜 = 62.4 ∗ 0.77961 = 48,6476 lpc
Calculamos Z con Pwh y Twh donde Z se encuentra en estado bifásico ya que Pd>Pwh Z=0.827558
Nº
iteraciones Pm* Z α F Im Pm calcu
1 2000 0,895054377 86520,225 0,064300072 247,079726 2534,227983
2 2534,227983 0,887620493 86520,225 0,064300072 220,3921254 2545,957538
3 2545,957538 0,887587461 86520,225 0,064300072 219,8184922 2546,225643
4 2546,225643 0,887586772 86520,225 0,064300072 219,8053996 2546,23177
5 2546,23177 0,887586756 86520,225 0,064300072 219,8051004 2546,23191
Nº
iteraciones Pwf* Z α F Iwf Pwf calcu
1 3300 0,94270392 86520,225 0,064300072 205,504923 2749,660376
2 2749,660376 0,934920182 86520,225 0,064300072 227,2551866 2739,763219
3 2739,763219 0,934853826 86520,225 0,064300072 227,6704968 2739,583599
4 2739,583599 0,934852646 86520,225 0,064300072 227,678039 2739,58034
5 2739,58034 0,934852625 86520,225 0,064300072 227,6781758 2739,580281
b)
Pfw del pozo X1=2739,58=2740 psi
𝑸𝒈
𝑰𝒑 =
𝑷𝒓𝟐 − 𝑷𝒘𝒇𝟐
Calculamos la presión estática de reservorio mediante el índice de productividad del pozo x1 donde la
presión estática de reservorio se puede considerar el mismo valor para los dos pozos X1 Y X2:
𝑄𝑔 10,5 𝑀𝑀𝑃𝐶𝐷
𝑃𝑟 = √ + 𝑃𝑓𝑤 2 = √ + 27402 = 𝟒𝟓𝟎𝟎 𝑷𝒔𝒊
𝐼𝑝 8,24E − 7 MMPCD/Psi2
c)
En este caso para calcular el caudal optimo del pozo X2 para una presión de cabeza de 2000
psi realizamos la gráfica del IPR mediante el método de análisis de reservorio para datos de
prueba y graficamos la curva VLP mediante el método de caída de presión en estado
monofásico para gas seco calculando la pfw para diferentes caudales con el dato de la
presión de cabeza.
IPR METODO DE JONES BLOUND GLAZE
800000
600000
400000
y = 122,83x + 1E+06
200000
0
0 2 4 6 8 10 12 14
Qg MMPCD
A=1116655 y B=122,83
−𝐶 + √𝐶 2 + 4𝐷(𝑃𝑟 2 − 𝑃𝑤𝑓 2 )
𝑄𝑔 =
2𝐷
Pfw Qg MMPCD
4500 0
4050 3,444253491
3600 6,523744977
3150 9,239214308
2700 11,59131222
2250 13,58060111
1800 15,20755572
1350 16,47256373
900 17,37592616
450 17,91785778
0 18,09848735
IPR
5000
4500
4000
3500
3000
Pfw psi
2500
2000
1500
1000
500
0
0 2 4 6 8 10 12 14 16 18 20
Qg MMPCD
Seguidamente calculamos para cada caudal de gas la Pwf con la Pwh=1500 psi
con el método de P y T promedio:
Ejemplo para el caudal de 6,52 MMPCD nos suponemos una Pwf*=2000 psi
T Pwf
promedio P promedio Z u (cp) Nre f S calculado
639,5 2000 0,85355187 0,016729158 2611249,183 0,017619793 0,476634755 3072,445565
639,5 2536,222783 0,844431045 0,018498767 2361454,683 0,017641345 0,481782957 3074,578863
639,5 2537,289431 0,844429283 0,018502496 2360978,796 0,01764139 0,481783963 3074,580412
639,5 2537,290206 0,844429282 0,018502498 2360978,45 0,01764139 0,481783963 3074,580413
639,5 2537,290206 0,844429282 0,018502498 2360978,45 0,01764139 0,481783963 3074,580413
Pwf=3074,58 psi
De la misma forma realizamos para cada caudal del IPR el cálculo de la Pwf para una
presión de cabeza de 200 psi con el método de presión y temperatura promedio.
IPR VLP
6000
5000
4000
Pfw psi
3000
2000
1000
0
0 2 4 6 8 10 12 14 16 18 20
Qg MMPCD
El caudal óptimo para el pozo X2 y para una presión de cabeza de 2000 psi Qgx2=8,3 MMPCD que es la
intersección de la gráfica IPR y VLP
MISMO EJERCICIO RESUELTO CON EL SIMULADOR PIPESIM
Qg =8,0381 MMPCD
EJEMPLO
Un reservorio de gas el cual consta de tres pozos productores donde el primer pozo (X1) convencional se
encuentra a una profundidad de 167 pies mayor que el pozo (X2) ya que el pozo X2 se encuentra a 1456
ft sobre el nivel del mar donde la trayectoria del pozo consiste en una sección de incremento el cual comienza
a una profundidad de 7495 pies hasta una profundidad de 8821 pies seguida de una sección donde el ángulo
se mantiene constante hasta una profundidad de 13268 pies en el cual comienza a disminuir a un ritmo
determinado donde generar una sección vertical de 935 pies con el cual se obtiene un desplazamiento
horizontal hasta el objetivo de 4124 pies y un MD de 16479 ft ya que a una profundidad de 11567 ft del
pozo x2 se registra una presión de 3602 psi y el pozo X3 de tipo Sland cuenta con un Kop de 2865 ft donde
el ritmo de incremento es de 1,6 º/100 ft y un MD de 17589 ft dese superficie hasta el objetivo.
El pozo X3 cuenta con una tubería de 3 ¼” y un radio del pozo de 4 ½” con un caudal de 12,7 MMPCD con un
índice de productividad de 3,178 PC/psi^2 ya que cuenta con una rugosidad de 0,0023” donde la gravedad
especifica del gas es de 0,67 con un API de 51º y la RGC de 205678 PC/Bbl ya que la temperatura en cabeza
es de 65ºF donde se considera un gradiente geotérmico de 16,5ºF/1000 ft
El pozo x2 produce gas con una gravedad especifica de 0.66 y un API de 50º ya que está conformado con una
tubería de producción de 2 ½ ” con un radio de pozo de 3 ¼” donde se tiene una rugosidad de 0.0028 in Y
también se registra una presión del separador de 1368 psi ya que la temperatura en cabeza del pozo es de
h) Realizar el análisis de caudales críticos y el caudal optimo en el pozo X2 para una presión de cabeza
de 2500 psi.
i) Calcular la caída de presión presión desde el reservorio hasta el separador
j) Calcular el AOF del pozo X3
Considerar para el cálculo de caída de presión el método de cullender para el pozo X1 y X3 Método
de P y T promedio para el pozo X2 y jhones Blount and Glaze para datos de prueba y Darcy en caso
de datos de reservorio
SOLUCION
𝜃 𝐷3 − 𝐷2 𝜃
𝑀𝑉𝐷𝑋2 = 𝐷1 + + + + 935 → 𝐷2 = 𝐷1 + 𝑟1 ∗ sin 𝜃
𝑞1 cos 𝜃 𝑞2
𝜃 13268 − 8821 𝜃
16479 = 7495 + + + + 935
𝑞1 cos 𝜃 𝑞2
𝟑, 𝟓 𝟐, 𝟔
𝒒𝟐 = 𝒇𝒕 𝒒𝟏 = 𝒇𝒕 𝜽 = 𝟑𝟕º 𝒓𝟏 = 𝟐𝟐𝟎𝟑, 𝟕 𝒓𝟐 = 𝟏𝟔𝟑𝟕
𝟏𝟎𝟎 𝟏𝟎𝟎
Calculamos el TVD del pozo X1 y X2
Comenzamos por el pozo X1 a calcular la Pfw en estado dinámico donde mediante el Índice de productividad de
obtendrá le presión de reservorio con el cual se aplicará el método Jones and Glaze para trabajar en datos de prueba
0.10796∗𝑄𝑔 2 0.10796∗13 2
𝐹2 = ( ) =( ) = 0.006338 Para este caso como se parte desde cabeza hasta el fondo ya que
𝐷 2.612 32.612
el pozo es vertical TVD=MD donde TVDx1=H=15355 ft y donde h=H/2=7677,5 ft para una Twh=66ºF=525,67ºR y una
Twf=759ºR y para una Pwh=2970 psi
Con la presión de reservorio de 4500 psi que es válida para los dos pozos ya que atraviesan al mismo reservorio
realizamos el IPR para el pozo X2
DP/Qg
1200000
1000000
y = 121005x - 536776
DP/Qg Psi^2/MMPCD
800000
600000
400000
200000
0
0 2 4 6 8 10 12 14
Qg MMPCD
Pfw Qg MMPCD
4500 0
IPR tranciente
4000 8,545869936
5000
3000 12,11197556
4000
2000 14,01678555
Pwf psi
3000
1000 15,02439272
2000
0 15,34308334 1000
0
0 5 10 15 20
Qg MMPCD
IPR estabilizado
Pfw Qg MMPCD
5000
4500 0 4000
Del estado estabilizado calculamos realizamos la VLP para una presión de cabeza de 2500 psi en el pozo x2
determinar el caudal óptimo para el pozo X2
IPR vs VLP
5000
4500
4000
3500
3000
Pwf psi
2500
2000
1500
1000
500
0
0 2 4 6 8 10 12 14 16 18
Qg MMPCD
caudal óptimo de 8,3 MMPCD
Con la presión de fondo del pozo X2 calculamos la Pwh para caudales críticos donde determinamos la fase
mediante el cálculo de la presión de Roció
16º𝐹
𝑇𝑤𝑓𝑥2 = 65 + ( ) ∗ 15188 𝑓𝑡 = 308º𝐹 = 767,678º𝑅
1000
141.5
𝛾𝑜 = = 0,77961
131.5 + 50
4591 ∗ 𝛾𝑜 4591 ∗ 0,77961
𝛿𝑔 + 0,66 +
𝛿𝑚 = 𝑅𝐺𝐶 = 104131 = 0,68696
1123 1123
1+ 1+
𝑅𝐺𝐶 104131
Pwh del pozo X2 Pwh=2500 psi
Realizamos el análisis de caudales críticos para el caudal de 8,5 del
pozo X2
ρo = 62.4 ∗ 𝛾𝑜 = 62.4 ∗ 0.77961 = 48,6476 lpc
= 𝟑. 𝟖𝟕𝟔𝟗𝟗 𝒇𝒕/𝒔
Calculamos Z con Pwh y Twh donde Z se encuentra en estado bifásico ya que Pd>Pwh Z=0.705296
Ya que el caudal del pozo x2 es mayor a los caudales críticos la producción es óptima 8,5
MMPCD>3,95MMPCD>2,72MMPCD
Aplicamos nuevamente la ecuación para el pozo X3 con el valor obtenido de K*h ya que se está hablando del mismo
reservorio ya que estamos en estado Monofasico la ecuación es no acumulativa.
𝐾 ∗ ℎ(𝑃𝑟 2 − 𝑃𝑤𝑓𝑥32 )
𝐴𝑂𝐹 = 𝑄𝑔𝑥3 = 𝑟𝑒
1424 ∗ 𝑇𝑟 ∗ 𝜇𝑔 ∗ 𝑧 (ln (𝑟𝑤𝑥3) − 0,75 + 𝑠)
𝐾 ∗ ℎ(𝑃𝑟 2 )
𝐴𝑂𝐹 = 𝑄𝑔𝑥3 = 𝑟𝑒
1424 ∗ 𝑇𝑟 ∗ 𝜇𝑔 ∗ 𝑧 (ln (𝑟𝑤𝑥3) − 0,75 + 𝑠)
EJEMPLO
Un reservorio que costa de dos formaciones productoras que producen gas mediante dos pozos donde cada uno
atraviesa una formación . Ya que se lleva a cabo una prueba de flujo tras flujo en la formación de gas ICLA que
produce gas mediante el pozo X1 con el cual se cuenta con los siguientes datos:
La formacion de gas Santa Rosa el cual atraviesa el pozo X2 con un diámetro de 4 in ya que cuenta con un factor
de daño de 6 donde se registra una porosidad de 32% a una profundidad de 7940 pies con una temperatura de
256 F y una presión estática de 5760 psi donde la permeabilidad absoluta es de 14 md y la permeabilidad efectiva
del agua es igual a la del condensado que equivale a 6 md en el cual se determinara un espaciamiento de 78 acres
del reservorio con un espesor neto de 36 ft de la formación productora del cual se obtuvo una gravedad especifica
del gas de 0.66 y un gravedad especifica del condensado de 0.7796 en el cual indican que la relación gas
condensado es de 116789 pc por barril. Realizar el IPR compuesto por el método de LIT
Datos recompilados
Calculo de RGC
𝐐𝐠
RGC =
𝐐 𝐜𝐨𝐧𝐝𝐞𝐧𝐬𝐚𝐝𝐨
Q condensado
Qg(MMPCD)
(BBP) RGC
4.465 64 69765.625
5.78 73 79178.0822
6.51 78 83461.5385
7.26 81 89629.6296
6.45 71 90845.0704
promedio 82575.99
Calculo de SGo
141.5
SGo = = 0.759
131.5 + 55
Calculo de SGmix
4591 ∗ 0.759
0.65 +
SGmix = 82575.99 = 0.683
1123
1+
82575.99
Calculo la precion de roció
182575.99 −0.8207
%C7+= ( ) = 0.88
70680
82575.990.09741 0.28126 −0.00068∗0.881.90633 )
Pd = 346.77648 ∗ [ −0.29478 ∗ 8.41763 ∗ 55(−0.04783∗245 ] = 3433.9
0.88
Tpc 371.291 °R
Tpc 373.388 °R
Calculo de Z
Z monofasico
1.383
μ gas = 10−4 ∗ 143.8997 ∗ e5.0869∗0.2507 = 0.0305
Pwf ρ
Fase Ppr B Z μ gas(cp)
[Psi] gas(gr/cm3)
Po P1
𝑀𝑃 = ( + ) + (𝑃𝑜 − 𝑃1)
Uo ∗ Zo U1 ∗ Z1
7500 6250
𝑀𝑃 = ( + ) + (7500 − 6250) = 5.15𝐸8
1.193 ∗ 0.0305 0.0279 ∗ 1.090
Calculo de MP acumulado
Pwf MP
Fase P/U*Z MP
[Psi] acumulado
1000,00
500,00
0,00
0 1000 2000 3000 4000 5000 6000 7000 8000
a=3E-5
b=0.1192
MP=3E-5*7500^2+0.1192*7500=2581.5
𝑀𝑃 2581.5 − 2059.6
= = 116.88
𝑄 4.465
Pwf
Fase Qg(MMPCD) MP acumulado
(psia) MP/Q
inicio 7500 0 2581.5 -
1 6534 4.465 2059.6 116.88
2 5213 5.78 1436.7 198.07
3 4756 6.51 1245.5 205.22
4 3720 7.26 858.6 237.32
extendido 4295 6.45 1065.4 235.06
Estado transiente
MP/Qg VS Qg
300
250
y = 41.827x - 61.745
R² = 0.944
200
150
100
50
0
0 1 2 3 4 5 6 7 8
A=-61.745
B=41.827
B=41.827
7000
6000
5000
4000
3000
2000
1000
0
0.000 0.002 0.004 0.006 0.008 0.010
Qg transiente Qg estabilizado
(MMPCD) (MMPCD)
Datos recompilados
PM=28.97*0.66=19.1202
Calculo de la permeabilidad del gas
𝐾 = 𝐾𝑔 + 𝐾𝑜 + 𝐾𝑤
14 = 𝐾𝑔 + 6 + 6
Krg=2/14=0.1428
Kg=2md
43560 ∗ 78
𝑟𝑒 = √ = 1039.96 𝑓𝑡
𝜋
Rw=4 in=0.33 ft
Calculo de API
141.5
𝐴𝑃𝐼 = − 131.5 = 50
0.7796
Calculo de SGmix
4591 ∗ 0.7796
0.66 +
SGmix = 116789 = 0.6841
1123
1+
116789
116789 −0.8207
%C7+= ( ) = 0.662
70680
1167890.09741 0.28126 −0.00068∗0.6621.90633 )
Pd = 346.77648 ∗ [ −0.29478 ∗ 8.41763 ∗ 50(−0.04783∗256 ] = 3303.65
0.662
Propiedades Pseudocriticas:
Tpc 373.655 °R
Tpc 373.661 °R
Calculo de Z
Pwf ρ
Fase Ppr B Z μ gas(cp)
[Psi] gas(gr/cm3)
Monofásic
Pwf Mp M(P)
Fase P/U*Z MP
[Psi] acumulado [Mmpsi/cp]
5760 202530.3 3.85E+08 1.63E+09
Monofási
1616.6
4800 198057.7 3.70E+08 1.24E+09
co
1230.9
3840 186933.7 3.34E+08 8.75E+08 861.84
2880 160742.2 2.76E+08 5.40E+08 530.7
Bífasico
monofasico
Pwf Qg Qg Qg
Fase
[Psi] [MSCFD] [MMSCFD] acumulado
5760
Monofásic
0.0 0 0.000
4800 2049.4 2.0494 2.049
o
8000
7000 Qg estabilizado Qg Q*
(MMPCD) acumulado
6000
5000
4000
3000
2000
1000
0
0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 0.018
Realizar el análisis de caudales críticos para la Pwf=5500 psi considerando una caída de presión de 1200 psi
Pwf=5500
Twh=65 F=525 R
Δ𝑃 ∗ = 𝑃𝑤𝑓 − 𝑃𝑤ℎ
Pwh=4300 psi
Calculo de Z
Tpc 371.291 °R
Ppr 6.440
Ppr=4300/667.705=6.440
Tpr=525/371.291=1.414
Z monofasico
Z = 0.8685
3.06 ∗ 4.119 ∗ 0.2673 ∗ 4300
𝑄𝑔 − 𝑎𝑔𝑢𝑎 = = 31.77 𝑀𝑀𝑃𝐶𝐷
525 ∗ 0.8685
3.06 ∗ 2.869 ∗ 0.2673 ∗ 4300
𝑄𝑔 − 𝑎𝑔𝑢𝑎 = = 22.13 𝑀𝑀𝑃𝐶𝐷
525 ∗ 0.8685
Calculo para el caudal a Pwf=5500
MP=3E-5*5500^2+0.1192*5500=1563.1
B=41.827
A=-34.72415 B=41.827
Un reservorio de gas el cual cuenta con dos formaciones productoras consta de dos pozos de producción en el cual
cada uno atraviesa una formación ya que se lleva a cabo una prueba isocronal de producción en la formación B el
cual cuenta con el pozo convencional X2 en el cual se obtienen los siguientes datos:
FORMACION ( A ) Y POZO X1
Pwh=4658 psi Pr=5663 psi Pws =5517 psi Twh=68ºF
API=43 º e=53 Acre RGC=78000PC/Bbl Cg=1.4752E-6 psi-1
ºG=18ºF/1000 ft MD=5678 m TVD=4356 m rw=3”
Swc =25% Cf = 3.1121E-6 psi-1 Kw=8.3 md Ko=6 md
COMPOSICION
Componentes %másica
C1 71%
C2 16%
C3 8%
CO2 5%
FORMACION ( B ) Y POZO X2
Sg=0.63 Pr=6275 psi Twh=68ºF Tr=262ºF
API=42 º e=67 Acre TVD=15397 ft Cg=1.8932E-6 psi-1
a) Realizar el IPR y calcular el AOF de la formación A mediante el método de lit 40%
b) Calcular el AOF de la formación B mediante el método simplificado
en estado tranciente y estabilizado 25%
c) Calcular la Presión de Rocío en la formación A 15%
d) Realizar el IPR compuesto del reservorio de gas 20 %
SOLUCION FORMACION A
Convertir la composición másica a molar
COMPOSICION
Componentes %molar
C1 84.26%
C2 10.13%
C3 3.45 %
CO2 2.16%
%𝐶7 +
𝐶7+=
100
𝑃𝑑 = 𝐸𝑋𝑃 {𝐴[0.2 ∗ %𝑁2 + %𝐶𝑂2 + %𝐻2𝑆 + 0.4 ∗ %𝑀𝑒𝑡ℎ + %𝐸𝑡ℎ + 2(%𝑃𝑟𝑜𝑝 + %𝐼𝐵𝑢𝑡 − %𝑁𝐵𝑢𝑡) + %𝐼𝑃𝑒𝑛
%𝑀𝑒𝑡ℎ
+ %𝑁𝑃𝑒𝑛 + %𝑁𝐻𝑒𝑥] + 𝐵 ∗ 𝑫𝑬𝑵 𝑪𝟕 + +𝐶 [ ] + 𝐷 ∗ 𝑇𝑟 + 𝐸 ∗ 𝐿 + 𝐹 ∗ 𝐿2 + 𝐺 ∗ 𝐿3
%𝑪𝟕 + +0.2
+ 𝐻 ∗ 𝑀 + 𝐼 ∗ 𝑀2 + 𝐽 ∗ 𝑀3 + 𝐾}
Donde:
𝐴 = −2.0623054 ∗ 10−2
𝐵 = 6.6259728
𝐶 = −4.4670559 ∗ 10−3
𝐷 = 1.0448346 ∗ 10−4
𝐸 = 3.2673714 ∗ 10−2
𝐹 = −3.6453277 ∗ 10−3
𝐺 = 7.4299951 ∗ 10−5
𝐻 = −0.11381195
𝐼 = 6.2476497 ∗ 10−4
𝐽 = −1.0716866 ∗ 10−6
𝐾 = 10.746622
Simplificando en función a los componentes
𝐿=0 𝑀=0
%𝑀𝑒𝑡ℎ
𝑃𝑑 = 𝐸𝑋𝑃 {𝐴[%𝐶𝑂2 + 0.4 ∗ %𝑀𝑒𝑡ℎ + %𝐸𝑡ℎ + 2(%𝑃𝑟𝑜𝑝)] + 𝐶 [ ] + 𝐷 ∗ 𝑇𝑟 + 𝐾}
0.2
84.26
𝑃𝑑 = 𝐸𝑋𝑃 {𝐴[2.16 + 0.4 ∗ 84.26 + 10.13 + 2(3.45)] + 𝐶 [ ] + 𝐷 ∗ 257.23 + 𝐾}
0.2
Pd=2502 psi
Calculamos Sgm de mezcla
𝟏𝟒𝟏. 𝟓 𝟏𝟒𝟏. 𝟓
𝜸𝒐 = 𝒐
= = 𝟎. 𝟖𝟏𝟎𝟖
𝟏𝟑𝟏. 𝟓 + 𝑨𝑷𝑰 𝟏𝟑𝟏. 𝟓 + 𝟒𝟑
𝟕𝟎𝟑∗𝟏𝟎−𝟔 ∗𝑲∗𝒉∗(𝒎(𝒑𝒓)−𝒎(𝒑𝒇𝒘))
𝑸𝒈 = 𝒓𝒆 (MPCD)
𝑻∗(𝐥𝐧(𝒓𝒘)−𝟎.𝟕𝟓+𝑺)
Estado Bifásico:
𝟕𝟎𝟑∗𝟏𝟎−𝟔 ∗𝑲∗𝑲𝒓𝒈∗𝒉∗(𝒎(𝒑𝒓)−𝒎(𝒑𝒇𝒘))
𝑸𝒈 = 𝒓𝒆 (MPCD)
𝑻∗(𝐥𝐧(𝒓𝒘)−𝟎.𝟕𝟓+𝑺)
Pfw (psi) Qg MMPCD Qg acumulad MMPCD
5663 0 0
4854 8.84695 8.84695
4045 17.38365 26.23060
3263 25.08846 51.31906
2427 8.35859 59.67765
1618 9.79943 69.47708
809 10.78348 80.26056
0 11.14156 91.40212
SOLUCION FORMACION B
SIMPLIFICADO
1000
100
PFW2 MMPSI2
10
1
1 10 100 1000
QG MMPCD
𝐼2 ∗ 𝑃𝑤𝑓22 + 𝐼1 ∗ 𝑃𝑤𝑓12
𝑃𝑤𝑓 ∗ = √
𝐼2 + 𝐼1
Se desea realizar un estudio al reservorio que produce de una arena no consolidada para implementar un sistema
de levantamiento artificial a fututo considerando que se generó una caída de presión del 27% donde actualmente
se registra un índice de productividad de 0,732 PC/psi^2 a una presión de fondo de 3456 psi y se cuentan con los
siguientes datos:
COMPOSICION
Componentes %molar
CH4 70%
C2H6 10%
C3H8 7%
n-C4H10 3%
CO2 6%
N2 4%
Calcular:
A) Ya que la temperatura del reservorio es mayor a 320ºF la correlación para el cálculo de la presión
de rocío no es aplicable debido a que es máximo para una temperatura de 320ºF con el cual se
procede aplicar la correlación de Calgary para calcular la presión de rocío.
𝐺𝐶𝑅𝐾2 𝐾5 −𝐾6∗%𝐶7𝐾7 )
𝑃𝑑 = 𝐾1 [ ∗ 𝐾8 ∗ 𝐴𝑃𝐼 (𝐾4∗𝑇𝑟 ] Tpromedio=182 (ºF) GCR (Pc/bbl) %𝐶7+=
%𝐶7𝐾3
𝐺𝐶𝑅 −0.8207
( )
70680
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
𝐾3 = −0.29478 𝐾7 = 1.90633
𝐾4 = −0.04783 𝐾8 = 8.41763
𝑷𝒅 = 𝟐𝟖𝟐𝟔. 𝟕𝟒 𝑷𝒔𝒊
B)
ENCONTRANDO LA COMPRESIBILIDAD DE FORMACION
𝑪𝑻 = 𝑪𝒐 ∗ 𝑺𝑶 + 𝑪𝒈 ∗ 𝑺𝒈 + 𝑪𝒘 ∗ 𝑺𝒘 + 𝑪𝑭 EC (1)
DESPEJANDO DE LA EC 1 LA “𝑪𝑭”
𝑪𝑭 = 𝑪𝑻 − 𝑪𝒐 ∗ 𝑺𝑶 − 𝑪𝒈 ∗ 𝑺𝒈 − 𝑪𝒘 ∗ 𝑺𝒘 EC (A)
𝑺𝒐 + 𝑺𝒈 + 𝑺𝒘 = 𝟏 EC (2)
DESPEJANDO EL “𝑺𝒈 ” DE LA EC (2) Y SUSTITUYENDO
𝐶𝐹 = 3.668𝐸 − 6𝑃𝑠𝑖 −1
∅ = 0.1923 = 19,23%
𝟐
∅𝟑
𝑲 = 𝟐. 𝟓 ∗ (𝑺 ) ∗ 𝟏𝟎𝟎𝟎 EC (B) 𝑲𝒈 = 𝑲 − 𝑲𝒐 − 𝑲𝒘
𝑾𝑪
2
0.19233
𝐾 = 2.5 ∗ ( ) ∗ 1000 = 10.44 𝑚𝑑
0.11
AOF=102 MMPCD
IPR
8000
6000
Pwf psi
4000
2000
0
0 20 40 60 80 100 120
Q g MMPCD
EJEMPLO
Un reservorio que costa de dos formaciones productoras que producen gas mediante dos pozos donde cada uno
atraviesa una formación. Ya que se lleva a cabo una prueba de flujo tras flujo en la formación de gas ICLA que
produce gas mediante el pozo X1 con el cual se cuenta con los siguientes datos:
La formación de gas Santa Rosa el cual atraviesa el pozo X2 con un diámetro de 4 in ya que cuenta con un factor
de daño de 6 donde se registra una porosidad de 32% a una profundidad de 7940 pies con una temperatura de
256 F y una presión estática de 5760 psi donde la permeabilidad absoluta es de 14 md y la permeabilidad efectiva
del agua es igual a la del condensado que equivale a 6 md en el cual se determinara un espaciamiento de 78 acres
del reservorio con un espesor neto de 36 ft de la formación productora del cual se obtuvo una gravedad especifica
del gas de 0.66 y un gravedad especifica del condensado de 0.7796 en el cual indican que la relación gas
condensado es de 116789 pc por barril. Realizar el IPR compuesto por el método de LIT
Datos recompilados
Calculo de RGC
𝐐𝐠
RGC =
𝐐 𝐜𝐨𝐧𝐝𝐞𝐧𝐬𝐚𝐝𝐨
Q condensado
Qg(MMPCD)
(BBP) RGC
4.465 64 69765.625
5.78 73 79178.0822
6.51 78 83461.5385
7.26 81 89629.6296
6.45 71 90845.0704
promedio 82575.99
Calculo de SGo
141.5
SGo = = 0.759
131.5 + 55
Calculo de SGmix
4591 ∗ 0.759
0.65 +
SGmix = 82575.99 = 0.683
1123
1+
82575.99
Calculo la precion de roció
182575.99 −0.8207
%C7+= ( ) = 0.88
70680
82575.990.09741 0.28126 −0.00068∗0.881.90633 )
Pd = 346.77648 ∗ [ ∗ 8.41763 ∗ 55(−0.04783∗245 ] = 3433.9
0.88−0.29478
Tpc 371.291 °R
Tpc 373.388 °R
Calculo de Z para una Pfw=3750 psi (bifásico) y 7500 psi (monofásico) y Tr=705ºR
Z monofásico correlacion de bexanbrill
1 − 0.5906
Z = 0.5906 + + 0.043 ∗ 11.2301.091 = 1.193
𝑒 5.412
Para el primer calculo 3750
−3.56539 1.53428 11.230
Z = 2.24353 + −0.0375281 ∗ 11.230 + + 0.00082922 ∗ 11.232 + 2
+ 0.131987 ∗ . = 0.994
1.888 1.888 1.888
Calculo de densidad del gas
0.65 ∗ 7500 ∗ 28.97
ρ gas = 0.016018463374 ∗ = 0.2507
10.732 ∗ 1.193 ∗ 705
Calculo de viscosidad del gas
K 143.8997 X 5.0869 Y 1.383
1.383
μ gas = 10−4 ∗ 143.8997 ∗ e5.0869∗0.2507 = 0.0305
Pwf ρ
Fase Ppr B Z μ gas(cp)
[Psi] gas(gr/cm3)
Po P1
𝑀𝑃 = ( + ) + (𝑃𝑜 − 𝑃1)
Uo ∗ Zo U1 ∗ Z1
7500 6250
𝑀𝑃 = ( + ) + (7500 − 6250) = 5.15𝐸8
1.193 ∗ 0.0305 0.0279 ∗ 1.090
Calculo de MP acumulado
Pwf MP
Fase P/U*Z MP
[Psi] acumulado
Monofásic
6250 205720.0 5.08E+08 1816.76
o
5000 201070.3 4.78E+08 1308.27
3750 180956.0 4.13E+08 830.73
Bífasico
Pwf Vs MP acumulado
2500,00
1000,00
500,00
0,00
0 1000 2000 3000 4000 5000 6000 7000 8000
a=3E-5
b=0.1192
MP=3E-5*pfw^2+0.1192*pfw
𝑀𝑃 2581.5 − 2059.6
= = 116.88
𝑄 4.465
Pwf
Fase Qg(MMPCD) MP acumulado
(psia) MP/Q
inicio 7500 0 2581.5 -
1 6534 4.465 2059.6 116.88
2 5213 5.78 1436.7 198.07
3 4756 6.51 1245.5 205.22
4 3720 7.26 858.6 237.32
extendido 4295 6.45 1065.4 235.06
Estado transiente
MP/Qg VS Qg
300
250
y = 41.827x - 61.745
R² = 0.944
200
150
100
50
0
0 1 2 3 4 5 6 7 8
A=-61.745
B=41.827
A=-34.72415
B=41.827
7000
6000
5000
4000
3000
2000
1000
0
0.000 0.002 0.004 0.006 0.008 0.010
Qg transiente Qg estabilizado
(MMPCD) (MMPCD)
Datos recompilados
PM=28.97*0.66=19.1202
Calculo de la permeabilidad del gas
𝐾 = 𝐾𝑔 + 𝐾𝑜 + 𝐾𝑤
14 = 𝐾𝑔 + 6 + 6
Krg=2/14=0.1428
Kg=2md
43560 ∗ 78
𝑟𝑒 = √ = 1039.96 𝑓𝑡
𝜋
Rw=4 in=0.33 ft
Calculo de API
141.5
𝐴𝑃𝐼 = − 131.5 = 50
0.7796
Calculo de SGmix
4591 ∗ 0.7796
0.66 +
SGmix = 116789 = 0.6841
1123
1+
116789
Calculo de la presio de rocio
116789 −0.8207
%C7+= ( ) = 0.662
70680
1167890.09741 0.28126 −0.00068∗0.6621.90633 )
Pd = 346.77648 ∗ [ −0.29478 ∗ 8.41763 ∗ 50(−0.04783∗256 ] = 3303.65
0.662
Propiedades Pseudocriticas:
Monofásico Ppc 667.043 Psi
Tpc 373.655 °R
Tpc 373.661 °R
Calculo de Z
Pwf ρ
Fase Ppr B Z μ gas(cp)
[Psi] gas(gr/cm3)
5760 8.635 3.4493 1.055 0.218 0.0268
Monofásic
Pwf Mp M(P)
Fase P/U*Z MP
[Psi] acumulado [Mmpsi/cp]
5760 202530.3 3.85E+08 1.63E+09
Monofási
1616.6
4800 198057.7 3.70E+08 1.24E+09
co
1230.9
3840 186933.7 3.34E+08 8.75E+08 861.84
2880 160742.2 2.76E+08 5.40E+08 530.7
Bífasico
monofasico
Pwf Qg Qg Qg
Fase
[Psi] [MSCFD] [MMSCFD] acumulado
5760
Monofásic
0.0 0 0.000
4800 2049.4 2.0494 2.049
o
3840 4010.2 4.0102 4.010
2880 823.9 0.82388 4.834
Bífasico
7000 Qg estabilizado Qg Q*
(MMPCD) acumulado
6000
5000
4000
3000
2000
1000
0
0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 0.018
Twh=65 F=525 R
Δ𝑃 ∗ = 𝑃𝑤𝑓 − 𝑃𝑤ℎ
Pwh=4300 psi
Calculo de Z
Tpc 371.291 °R
Ppr 6.440
Ppr=4300/667.705=6.440
Tpr=525/371.291=1.414
Z monofasico
Z = 0.8685
3.06 ∗ 4.119 ∗ 0.2673 ∗ 4300
𝑄𝑔 − 𝑎𝑔𝑢𝑎 = = 31.77 𝑀𝑀𝑃𝐶𝐷
525 ∗ 0.8685
3.06 ∗ 2.869 ∗ 0.2673 ∗ 4300
𝑄𝑔 − 𝑎𝑔𝑢𝑎 = = 22.13 𝑀𝑀𝑃𝐶𝐷
525 ∗ 0.8685
Calculo para el caudal a Pwf=5500
MP=3E-5*5500^2+0.1192*5500=1563.1
B=41.827
A=-34.72415 B=41.827
Hay acumulación de liquido en el fondo del pozo por lo cual la producción noes optima
EJEMPLO
Un reservorio de gas donde se tiene dos formaciones productoras cuenta con dos pozos productores que
atraviesan la misma formación productora donde el pozo X1 es un pozo direccional y el pozo X2 es un pozo
horizontal ya que en el pozo X1 se realiza las pruebas de producción periódica con el objetivo de analizar y
determinar las propiedades de la formación productora (A) en el cual se cuentan con los siguientes datos:
DATOS DEL RESERVORIO O LA FORMACION PRODUCTORA
ºG=27ºF/1000 ft e=73 Acre Cg=1.4752E-6 psi-1 Pr=7845 psi Cw=2.56E-6 psi^-1
Swc =31% Cf = 2.7121E-6 psi-1 Kw=38.75 md Kro=0.37 S=7
Swi =37% h=64 ft Ct = 4.4571E-6 psi-1 So=31% API=55º
POZO X1
Pwh=3452 psi Twh=70ºF MD=8105 ft TVD=14678 ft rw=6 ¼ “
COMPOSICION
Componentes %másica
C1 74.11
C2 13.81
C3 4.91
POZO X2
CO2 4.69
El pozo cuenta con un Kop=3643 ft, q1=1.8º/100ft el cual alcanza hasta una N2 2.48
profundidad de 6574 ft, q2=2.4º/100ft inicia a una profundidad de 7861 ft ya que cuenta con un
desplazamiento horizontal hasta el objetivo de 6158 ft
reH =1153 ft KHg=18.15 md rw=2 ½” Twh=68ºF RGC=784952 PC/Bbl SoH=0.24
POXO X1
Calculamos con el TVD la Twfx1 donde RGPpromedio=814065
𝑻𝒘𝒙𝟏 = 𝑻𝒘𝒉 + 𝑮º ∗ 𝑻𝑽𝑫 = 𝟒𝟔𝟔, 𝟐𝟓º𝑭 = 𝟗𝟐𝟔º𝑹
Calculamos la presión de rocío por el método de Calgari ya que la correlación de
cromatografía cuenta con un límite máximo de aplicación para una temperatura de
320ºF y la temperatura de reservorio es de 466ºF
1ro calcular %C7+ (Heptano plus)
𝐺𝐶𝑅 −0.8207
%𝐶7+= ( ) 𝑜 %𝐶7+= 10260 ∗ [(𝐺𝐶𝑅 ∗ 𝛿𝑔)−0.8499 ]
70680
2do Calcular la Presión de rocío
𝐺𝐶𝑅𝐾2 𝐾5 −𝐾6∗%𝐶7𝐾7 )
𝑃𝑑 = 𝐾1 [ ∗ 𝐾8 ∗ 𝐴𝑃𝐼 (𝐾4∗𝑇𝑟 ] Tr (ºF) GCR (Pc/bbl)
%𝐶7𝐾3
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
Pd=2067
Realizamos el cálculo de los pseudopotenciales corrigiendo el valor de Z por impurezas
Qg vs DMP/Qg
68
66
64
DMP/Qg
62
60
58
y = -0,8912x + 71,048
56
54
0 2 4 6 8 10 12 14 16 18
Qg MMPCD
A tranciente=71,08
B tranciente (Sele toma valor positivo) =0,8912
−𝐴 + √𝐴 2 + 4𝐵(𝑚(𝑃𝑟) − 𝑚(𝑃𝑤𝑓))
𝑄𝑔 =
2𝐵
CALCULAMOS ( L)
𝐷2 − 𝐷1 180 1 180 1 180 1 180 1
sin 𝜃 = 𝑑𝑜𝑛𝑑𝑒 𝑟1 = ∗ = ∗ = 3183 𝑓𝑡 𝑦 𝑟2 = ∗ = ∗
𝑟1 𝜋 𝑞 𝜋 1.8 𝜋 𝑞 𝜋 2.4
100 𝑓𝑡 100 𝑓𝑡
= 2387 𝑓𝑡
6574 − 3643
𝜃 = sin−1 ( ) = 67º 𝑑𝑜𝑛𝑑𝑒 𝜃 + 𝛽 = 90º 𝑝𝑜𝑟 𝑙𝑜 𝑡𝑎𝑛𝑡𝑜 𝛽 = 23º
3183
𝑋 = 𝑟1 ∗ (1 − 𝑐𝑜𝑠𝜃) + tan(𝜃) ∗ (𝐷3 − 𝐷2) + 𝑟2 ∗ sin 𝛽 + 𝐿
Calculamos TVD
𝑇𝑉𝐷 = 𝐷3 + 𝑟2 ∗ (1 − 𝑐𝑜𝑠𝛽) = 7861 + 2387(1 − 𝑐𝑜𝑠23º) = 8051 𝑓𝑡
Calcular la porosidad del Cf = 2.7121E-6 psi-1 para próximamente calcular K absoluta y después Kg
𝟏. 𝟕𝟖𝟐
𝑪𝒇 = ( 𝟎,𝟒𝟑𝟖 ) ∗ 𝟏𝟎−𝟔 𝑷𝒔𝒊−𝟏 𝒅𝒐𝒏𝒅𝒆 𝒔𝒆 𝒅𝒆𝒔𝒑𝒆𝒋𝒂 ∅ = 𝟎. 𝟑𝟖𝟑𝟑 = 𝟑𝟖. 𝟑𝟑%
∅
Calcular la permeabilidad absoluta con la porosidad y Swc=25%
𝟐 𝟐
∅𝟑 𝟎. 𝟑𝟖𝟑𝟑𝟑
𝑲 = 𝟐, 𝟓 ∗ ( ) ∗ 𝟏𝟎𝟎𝟎 ( 𝒎𝒅) = 𝑲 = 𝟐, 𝟓 ∗ ( ) ∗ 𝟏𝟎𝟎𝟎 = 𝟖𝟐. 𝟓 𝒎𝒅
𝑺𝑾𝑪 𝟎. 𝟑𝟏
𝑲𝒐
Calcular Ko donde 𝑲𝒓𝒐 = 𝒚 𝑲𝒐 = 𝑲 ∗ 𝑲𝒓𝒐 = 𝟑𝟎. 𝟓𝟐𝟓 𝒎𝒅
𝑲
𝐾𝐻 ∗ℎ(𝑃𝑟 2 −𝑃𝑤𝑓2 )
𝑄𝑔 = ( MPCD Monofasico)
𝐿 2
𝑎+√𝑎2 −(2)
𝐼𝑎𝑛𝑖 ∗ℎ 𝐼 ∗ℎ
1424∗𝜇𝑔∗𝑧∗𝑇𝑟 ln 𝐿 + [ln( 𝑎𝑛𝑖 )−0.75+𝑆]
𝐿 𝑟𝑤(𝐼𝑎𝑛𝑖 +1)
2
{ [ ] }
Donde:
𝐾𝐻 18.15
𝐼𝑎𝑛𝑖 = √ =√ = 3.1623
𝐾𝑣 1.815
0.5
4 0.5
254 1153 𝐿
𝑎= 0.5 + [0.25 + ( ) ] = 1156.5 𝑝𝑎𝑟𝑎 < 0.9𝑟𝑒𝐻
2 254 2
{ 2 }
Calculamos el caudal con la ecuación del pozo horizontal y para una presión de reservorio de 7845 psi
𝟕𝟐𝟕𝟏𝟔𝟎𝟎𝟎 𝒑𝒄
𝑰𝒑 = = 𝟐, 𝟖𝟒𝟔𝟕 𝑷𝒄/𝒑𝒔𝒊^𝟐
𝟕𝟖𝟒𝟓^𝟐 − 𝟔𝟎𝟎𝟎^𝟐
EJEMPLO DE FLUJO DINAMICO y ESTATICO
Un reservorio de gas el cual consta de dos pozos productores donde el primer pozo (X1) se encuentra a
una profundidad de 760 pies mayor que el segundo pozo (X2). El primer pozo cuenta con un MD de 20284
pies y un ángulo de 44º ya que la trayectoria del segundo pozo consiste en una sección de incremento el
cual comienza a una profundidad de 7460 pies a un ritmo de 1.4º por cada cien pies seguida de una sección
donde el ángulo se mantiene constante hasta una profundidad de 13640 pies en el cual comienza a
disminuir a un ritmo de 2.3º por cada cien pies hasta generar una sección vertical de 1235 pies con el cual
se obtiene un desplazamiento horizontal de 6145 pies.
El caudal de gas del pozo X1 es de 4.5 MMPCD con una gravedad especifica de 0.67, una relación de gas
condensado de 43678 PC/Bbl y un API de 36º ya que está conformado con una tubería de producción de
2.44” donde se tiene una rugosidad de 0.0017 Y también se registra una presión de fondo fluyente de 5782
psi a una temperatura de fondo fluyente registrada en los dos pozos de 216ºF y 68ºF en cabeza.
El pozo x2 produce un caudal de 5.7 MMPCD con un índice de productividad de 1.4E-7 MMPC/Psi2 a una
presión de cabeza de 1457 Psi y una gravedad especifica del gas de 0.64 con un API de 39º en el cual el
diámetro es de 3.2” con una rugosidad de 0.0026.
a) Determinar la presión de cabeza del pozo x1.
b) Determinar la presión de fondo del pozo x2 en estado estático si se procede a cerrar el pozo con
una presión de cabeza de 1800 psi. Considerar método de cullender.
SOLUCION
(X2) (X1)
Pwh Pwh
D1
D1
𝜃 r1 D2
TVD D3
D3 MD
TVD
D4 r2
935 ft
X4 Pfw Pfw
RESERVORIO
Pr
180 1
𝑟1 = ∗ = 4092.55 𝑓𝑡 𝑑𝑜𝑛𝑑𝑒 𝑋2 = 𝑟1 ∗ (1 − cos 𝜃) = 4092.55 ∗ (1 − cos 𝜃)
𝜋 1.4º
100 𝑓𝑡
180 1
𝑟2 = ∗ = 2491.12 𝑓𝑡 𝑑𝑜𝑛𝑑𝑒 𝑋3 = 𝑋2 + tan 𝜃 ∗ (𝐷3 − 𝐷2)
𝜋 2.3º
100 𝑓𝑡
(𝐷4 − 𝐷3)
𝑋4 = 𝑋3 + (𝑟2 − ) 𝑑𝑜𝑛𝑑𝑒 𝐷2 = 𝐷1 + 𝑟1 ∗ sin 𝜃 𝑦 𝐷4 = 𝐷3 + 𝑟2 ∗ sin 𝜃
tan 𝜃
𝟔𝟏𝟒𝟓 = 𝟒𝟎𝟗𝟐. 𝟓𝟓 ∗ (𝟏 − 𝐜𝐨𝐬 𝜽) + 𝐭𝐚𝐧 𝜽 ∗ (𝟏𝟑𝟔𝟒𝟎 − (𝟕𝟒𝟔𝟎 + 𝟒𝟎𝟗𝟐. 𝟓𝟓 ∗ 𝐬𝐢𝐧 𝜽)) + 𝟐𝟒𝟗𝟏. 𝟏𝟐 ∗ (𝟏 − 𝐜𝐨𝐬 𝜽)
Una ecuación con una incógnita del cual obtenemos 𝜽 = 𝟓𝟏. 𝟎𝟎º
𝑇𝑉𝐷 = 𝐷4 + 1235 = 𝐷3 + 𝑟2 ∗ sin 𝜃 + 1235 = 13640 + 2491.12 ∗ sin 51º + 1235
𝜃 𝐷3 − 𝐷2 𝜃
𝑀𝐷𝑋2 = 𝐷1 + + + + 1235 → 𝐷2 = 𝐷1 + 𝑟1 ∗ sin 𝜃 = 7460 + 4092.55 ∗ sin 51 = 10640.5 𝑓𝑡
𝑞1 cos 𝜃 𝑞2
51 13640 − 10640.5 51
𝑀𝐷𝑋2 = 7460 + + + + 1235 = 𝟏𝟗𝟑𝟐𝟏. 𝟓 𝒇𝒕
1.4 cos 51 2.3
100 100
a)
INICIAMOS CALCULANDO LA PWH EN EL POZO x1
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
𝐾3 = −0.29478 𝐾7 = 1.90633
𝐾4 = −0.04783 𝐾8 = 8.41763
𝑃𝑑 = 4629.70 𝑃𝑠𝑖
TABLA ITERATIVA EN EL EXEL
b)
Iniciamos calculando la Pfw en estado estático para el pozo X2
Tabla de iteración y resultado en el Excel
EJEMPLO
Un reservorio de gas cuenta con un pozo direccional donde la trayectoria del pozo consiste en una sección
de incremento el cual se da a un ritmo de 1.2º por cada cien pies seguida de una sección donde el ángulo
se mantiene constante generando una sección vertical de 3716 ft que alcanza hasta la profundidad de
11950 pies en el cual se registra un desplazamiento horizontal en ese punto de 5250 pies ya que a partir
de la profundidad de 11950 ft comienza a disminuir a un ritmo de 1.45º por cada cien pies hasta generar
un sección vertical de 456 pies que llega al objetivo.
El pozo produce un caudal de gas es de 24 MMPCD donde el gas cuenta con una gravedad especifica de
0.65 con un API de 46º del condensado ya que está conformado con una tubería de producción de 3.44”
donde se tiene una rugosidad de 0.0029”. El reservorio cuenta con una presión estática de 4789 psi donde
la temperatura de reservorio es de 321ºF en el cual se registra una temperatura de 72ºF en superficie con
una presión del separador de 956 psi. Considerando una presión de 3456 psi a la profundidad del Kop
calcular:
a) La presión en cabeza y fondo si se procede a cerrar el pozo
b) La presión a una profundidad de 7850 ft.
SOLUCION
Twh=78ºF
Pwh=?
D1
Pkop=? 𝜃 r1
MD=?
D3=11950
X3=5250ft r2
𝜃 q2=1.55
1134 ft
Pfw
Tr=287ºF
RESERVORIO
𝜽 = 𝟒𝟓. 𝟎𝟎º
𝑇𝑉𝐷 = 𝐷3 + 𝑟2 ∗ sin 𝜃 + 1134 = 11950 + 3696.5 ∗ sin 45º + 1134 = 15697.8 𝑓𝑡
b) calculamos la Pkop
∆𝑍 = 𝐷3 − 𝐷1 = 11950 − 4550.9 = 7399.1
𝜃 𝐷3 − 𝐷2
𝐿= + = 9346.1 𝑓𝑡
𝑞1 cos 𝜃
RESULTADOS
PROPIEDADES GAS Condensado Agua Velocidades
Densidad 13.28392 37.20010 59.69234 Vg 24.92072
Viscocidad 0.025732104 0.25517 0.249604 Vw 0.24628
factor bolumetrico 0.003971 1.73171 1.04536 Vo 0.79504
CONSTANTES DEL HL Vm 25.96204
N1 844004.1767 Nre 36762.80678
N2 106946.8584 f 0.023333
Calculamos la Pfw
∆𝑍 = 𝑇𝑉𝐷 − 𝐷3 = 3747.8 𝑓𝑡
𝜃
𝐿= + 1134 = 4037.22 𝑓𝑡
𝑞2
RESULTADOS
PROPIEDADES GAS Condensado Agua Velocidades
Densidad 13.28392 37.20010 59.69234 Vg 24.92072
Viscocidad 0.025732104 0.25517 0.249604 Vw 0.24628
factor bolumetrico 0.003971 1.73171 1.04536 Vo 0.79504
CONSTANTES DEL HL Vm 25.96204
N1 844004.1767 Nre 36762.80678
N2 106946.8584 f 0.023333
𝑇𝑤𝑓 − 𝑇𝑤ℎ
= 𝐺º = 0.01331º𝐹/𝑓𝑡
𝑇𝑉𝐷
𝑇𝑘𝑜𝑝 = 78 + 11950 ∗ 0.01331 = 138.57º𝐹 = 598.2 º𝑅
∆𝑍 = 𝐿 = 𝐷1 = 4550.9
RESULTADOS
PROPIEDADES GAS Condensado Agua Velocidades
Densidad 14.56525 38.90426 61.79853 Vg 22.72840
Viscocidad 0.026124498 0.39146 0.515404 Vw 0.23789
factor bolumetrico 0.003621 1.65481 1.00973 Vo 0.75973
CONSTANTES DEL HL Vm 23.72602
N1 663293.866 Nre 29263.80444
N2 104304.189 f 0.024580
SOLUCION
Twh=78ºF
Pwh=?
D1
Pkop=? 𝜃 r1
MD=?
D3=11950
X3=5250ft r2
𝜃 q2=1.55
1134 ft
Pfw
Tr=287ºF
RESERVORIO
𝜽 = 𝟒𝟓. 𝟎𝟎º
𝑇𝑉𝐷 = 𝐷3 + 𝑟2 ∗ sin 𝜃 + 1134 = 11950 + 3696.5 ∗ sin 45º + 1134 = 15697.8 𝑓𝑡
b) calculamos la Pkop
∆𝑍 = 𝐷3 − 𝐷1 = 11950 − 4550.9 = 7399.1
𝜃 𝐷3 − 𝐷2
𝐿= + = 9346.1 𝑓𝑡
𝑞1 cos 𝜃
RESULTADOS
PROPIEDADES GAS Condensado Agua Velocidades
Densidad 13.28392 37.20010 59.69234 Vg 24.92072
Viscocidad 0.025732104 0.25517 0.249604 Vw 0.24628
factor bolumetrico 0.003971 1.73171 1.04536 Vo 0.79504
CONSTANTES DEL HL Vm 25.96204
N1 844004.1767 Nre 36762.80678
N2 106946.8584 f 0.023333
Calculamos la Pfw
∆𝑍 = 𝑇𝑉𝐷 − 𝐷3 = 3747.8 𝑓𝑡
𝜃
𝐿= + 1134 = 4037.22 𝑓𝑡
𝑞2
RESULTADOS
PROPIEDADES GAS Condensado Agua Velocidades
Densidad 13.28392 37.20010 59.69234 Vg 24.92072
Viscocidad 0.025732104 0.25517 0.249604 Vw 0.24628
factor bolumetrico 0.003971 1.73171 1.04536 Vo 0.79504
CONSTANTES DEL HL Vm 25.96204
N1 844004.1767 Nre 36762.80678
N2 106946.8584 f 0.023333
RESULTADOS
PROPIEDADES GAS Condensado Agua Velocidades
Densidad 14.56525 38.90426 61.79853 Vg 22.72840
Viscocidad 0.026124498 0.39146 0.515404 Vw 0.23789
factor bolumetrico 0.003621 1.65481 1.00973 Vo 0.75973
CONSTANTES DEL HL Vm 23.72602
N1 663293.866 Nre 29263.80444
N2 104304.189 f 0.024580
FORMACION ( A ) Y POZO X1
Pwh=4658 psi Pr=5663 psi Pws =5517 psi Twh=68ºF
API=43 º e=53 Acre RGC=78000PC/Bbl Cg=1.4752E-6 psi-1
ºG=18ºF/1000 ft MD=5678 m TVD=4356 m rw=3”
Swc =25% Cf = 3.1121E-6 psi-1 Kw=8.3 md Ko=6 md
Soc =20% hp=20 ft Kp =39 md S=6
Sw=32% h=74 ft Cw = 4.4571E-6 psi-1
COMPOSICION
Componentes %másica
C1 71%
C2 16%
C3 8%
CO2 5%
FORMACION ( B ) Y POZO X2
Sg=0.63 Pr=6275 psi Twh=68ºF Tr=262ºF
API=42 º e=67 Acre TVD=15397 ft Cg=1.8932E-6 psi-1
COMPOSICION
Componentes %molar
C1 84.26%
C2 10.13%
C3 3.45 %
CO2 2.16%
𝑃𝑀 𝐶7 +
𝐿 = (𝐶7 +) ∗ (𝑃𝑀𝐶7 +) 𝑀=
𝐷𝐸𝑁 𝐶7 + +0.0001
%𝐶7 +
𝐶7+=
100
𝑃𝑑 = 𝐸𝑋𝑃 {𝐴[0.2 ∗ %𝑁2 + %𝐶𝑂2 + %𝐻2𝑆 + 0.4 ∗ %𝑀𝑒𝑡ℎ + %𝐸𝑡ℎ + 2(%𝑃𝑟𝑜𝑝 + %𝐼𝐵𝑢𝑡 − %𝑁𝐵𝑢𝑡) + %𝐼𝑃𝑒𝑛
%𝑀𝑒𝑡ℎ
+ %𝑁𝑃𝑒𝑛 + %𝑁𝐻𝑒𝑥] + 𝐵 ∗ 𝑫𝑬𝑵 𝑪𝟕 + +𝐶 [ ] + 𝐷 ∗ 𝑇𝑟 + 𝐸 ∗ 𝐿 + 𝐹 ∗ 𝐿2 + 𝐺 ∗ 𝐿3
%𝑪𝟕 + +0.2
+ 𝐻 ∗ 𝑀 + 𝐼 ∗ 𝑀2 + 𝐽 ∗ 𝑀3 + 𝐾}
Donde:
𝐴 = −2.0623054 ∗ 10−2
𝐵 = 6.6259728
𝐶 = −4.4670559 ∗ 10−3
𝐷 = 1.0448346 ∗ 10−4
𝐸 = 3.2673714 ∗ 10−2
𝐹 = −3.6453277 ∗ 10−3
𝐺 = 7.4299951 ∗ 10−5
𝐻 = −0.11381195
𝐼 = 6.2476497 ∗ 10−4
𝐽 = −1.0716866 ∗ 10−6
𝐾 = 10.746622
%𝑀𝑒𝑡ℎ
𝑃𝑑 = 𝐸𝑋𝑃 {𝐴[%𝐶𝑂2 + 0.4 ∗ %𝑀𝑒𝑡ℎ + %𝐸𝑡ℎ + 2(%𝑃𝑟𝑜𝑝)] + 𝐶 [ ] + 𝐷 ∗ 𝑇𝑟 + 𝐾}
0.2
84.26
𝑃𝑑 = 𝐸𝑋𝑃 {𝐴[2.16 + 0.4 ∗ 84.26 + 10.13 + 2(3.45)] + 𝐶 [ ] + 𝐷 ∗ 257.23 + 𝐾}
0.2
Pd=2502 psi
Calculamos Sgm de mezcla
𝟏𝟒𝟏. 𝟓 𝟏𝟒𝟏. 𝟓
𝜸𝒐 = 𝒐
= = 𝟎. 𝟖𝟏𝟎𝟖
𝟏𝟑𝟏. 𝟓 + 𝑨𝑷𝑰 𝟏𝟑𝟏. 𝟓 + 𝟒𝟑
𝟕𝟎𝟑∗𝟏𝟎−𝟔 ∗𝑲∗𝒉∗(𝒎(𝒑𝒓)−𝒎(𝒑𝒇𝒘))
𝑸𝒈 = 𝒓𝒆 (MPCD)
𝑻∗(𝐥𝐧(𝒓𝒘)−𝟎.𝟕𝟓+𝑺)
Estado Bifásico:
𝟕𝟎𝟑∗𝟏𝟎−𝟔 ∗𝑲∗𝑲𝒓𝒈∗𝒉∗(𝒎(𝒑𝒓)−𝒎(𝒑𝒇𝒘))
𝑸𝒈 = 𝒓𝒆 (MPCD)
𝑻∗(𝐥𝐧(𝒓𝒘)−𝟎.𝟕𝟓+𝑺)
Pfw (psi) Qg MMPCD Qg acumulad MMPCD
5663 0 0
4854 8.84695 8.84695
4045 17.38365 26.23060
3263 25.08846 51.31906
2427 8.35859 59.67765
1618 9.79943 69.47708
809 10.78348 80.26056
0 11.14156 91.40212
SOLUCION FORMACION B
SIMPLIFICADO
1000
100
PFW2 MMPSI2
10
1
1 10 100 1000
QG MMPCD
AOF estabilizado=29 MMPCD
AOF tranciente = 37 MMPCD
𝐺𝐶𝑅𝐾2 𝐾5 −𝐾6∗%𝐶7𝐾7 )
𝑃𝑑 = 𝐾1 [ ∗ 𝐾8 ∗ 𝐴𝑃𝐼 (𝐾4∗𝑇𝑟 ] Tr (ºF) GCR (Pc/bbl)
%𝐶7𝐾3
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
𝐾3 = −0.29478 𝐾7 = 1.90633
𝐾4 = −0.04783 𝐾8 = 8.41763
%𝐶7+= 10260 ∗ [(𝑅𝐺𝐶 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 ∗ 𝛿𝑔)−0.8499 ]
72514.6𝐾2 𝐾5 −𝐾6∗%𝐶7𝐾7 )
𝑃𝑑 = 𝐾1 [ ∗ 𝐾8 ∗ 58(𝐾4∗220 ] = 4169 Psi Tr (ºF) GCR (Pc/bbl)
%𝐶7𝐾3
𝑃1 𝑃𝑜
+ 0+0
𝜇1𝑧1 𝜇𝑜𝑧𝑜
𝑚(𝑃11) = 2 [ ] ∗ (𝑃1 − 𝑃0) + 0 = 2 [ ]∗ (0 − 0) + 0 = 0 (MMPsi2/cp)
2 2
𝑃2 𝑃1 69258.7+0
+ 2[ ]∗(784−0)
𝜇2𝑧2 𝜇1𝑧1 2
𝑚(𝑃10) = 2 [ ] ∗ (𝑃2 − 𝑃1) + 𝑚(𝑃11) = + 0 = 54.299 (MMPsi2/cp)
2 1000000
𝑃1 𝑃𝑜 116464.2+69258.7
+
𝜇1𝑧1 𝜇𝑜𝑧𝑜 2[ ]∗(1568−784)
2
𝑚(𝑃9) = 2 [ ] ∗ (𝑃1 − 𝑃0) + 𝑚(𝑃10) = + 54.299 = 199.906 (MMPsi2/cp)
2 1000000
𝑃2 𝑃1 146240.3−116464.2
+ 2[ ]∗(2352−1568)
𝜇2𝑧2 𝜇1𝑧1 2
𝑚(𝑃8) = 2 [ ] ∗ (𝑃2 − 𝑃1) + 𝑚(𝑃8) = + 199.906 = 405.866 (MMPsi2/cp)
2 1000000
REALIZAMOS EL AJUSTE DE GRAFICA
Pfw vs MP
3000,000
2000,000
MP MMpsi/cp
1500,000
1000,000
500,000
0,000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000
Pfw psi
4 ∗ 𝑎0 + 𝑎1 ∗ ∑ 𝑋 + 𝑎2 ∗ ∑ 𝑋2 = ∑ 𝑌 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 1
𝑎0 ∗ ∑ 𝑋 + 𝑎1 ∗ ∑ 𝑋2 + 𝑎2 ∗ ∑ 𝑋3 = ∑ 𝑋 ∗ 𝑌 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 2
𝑎0 ∗ ∑ 𝑋2 + 𝑎1 ∗ ∑ 𝑋3 + 𝑎2 ∗ ∑ 𝑋4 = ∑ 𝑋2 ∗ 𝑌 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 3
𝑎1 ∗ ∑ 𝑋2 + 𝑎2 ∗ ∑ 𝑋 3 = ∑ 𝑋 ∗ 𝑌 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 2
𝑎1 ∗ ∑ 𝑋3 + 𝑎2 ∗ ∑ 𝑋 4 = ∑ 𝑋 2 ∗ 𝑌 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 3
Pfw vs MP
3000,000
2500,000 y = 2E-05x2 + 0.1378x
MP MMpsi/cp
2000,000
1500,000
1000,000
500,000
0,000
0 2000 4000 6000 8000 10000
Pfw psi
MP recalculado
2520.534912
2136.57317
1781.924372
1456.58852
1160.565612
893.855648
656.4586291
448.3745549
269.6034253
120.1452403
0
FLUJO Pfw (psi) Qg (MMPCD) MP Mmpsi2/cp Delt MP DELTMP/Q
1 6720 3.2 1980.991488 539.543424 168.60732
2 5845 4.56 1601.097616 919.4372959 201.6309859
3 4651 5.98 1141.611647 1378.923265 230.5891748
4 3263 7.1 692.9249208 1827.609991 257.4098579
Extendido 5926 4 1634.731367 885.8035448 221.4508862
DELTMP/Q vs Qg
300
250
y = 22.497x + 97.349
DELTMP/Qg
200
150
100
50
0
0 1 2 3 4 5 6 7 8
Qg
A tranciente=97.349
B =22.479
Para el estado estabilizado
DELTMP/Q vs Qg
300
250
DELTMP/Qg
200
150
100
50
0
0 1 2 3 4 5 6 7 8
Qg
A estabilizado=131.53
B =22.479
Calculamos los caudales de gas Qg
−𝐴 + √𝐴 2 + 4𝐵 (𝑀𝑃𝑝𝑟 − 𝑀𝑃𝑝𝑓𝑤)
𝑄𝑔 =
2𝐵
Pfw MP recalculado Qg MMPC tranc Qg MMPC estabilizado
7840 2520.534912 0 0
7056 2136.57317 2.500455257 2.137946593
6272 1781.924372 3.962184981 3.509930405
5488 1456.58852 5.047111651 4.550272404
4704 1160.565612 5.90858947 5.384473539
3920 893.855648 6.612655115 6.070065576
3136 656.4586291 7.194897329 6.639058431
2352 448.3745549 7.6769601 7.111301536
1568 269.6034253 8.073003796 7.499950368
784 120.1452403 8.392712864 7.814084014
0 0 8.642863419 8.060096497
IPR tranciente
9000
8000
7000
6000
5000
Pfw
4000
3000
2000
1000
0
0 1 2 3 4 5 6 7 8 9 10
Qg
IPR ESTABILIZADO
9000
8000
7000
6000
5000
Pfw
4000
3000
2000
1000
0
0 1 2 3 4 5 6 7 8 9
Qg
EJEMPLO
Un reservorio de gas donde se tiene dos formaciones productoras, cuenta con dos pozos productores que
atraviesan la misma formación productora donde el pozo X1 es un pozo direccional y el pozo X2 es un pozo
horizontal ya que en el pozo X1 se realiza las pruebas de producción periódica con el objetivo de analizar y
determinar las propiedades de la formación productora (A) en el cual se cuentan con los siguientes datos:
POZO X1
POZO X2
El pozo cuenta con un Kop=3643 ft, q1=1.8º/100ft el cual alcanza hasta una profundidad de 6574 ft,
q2=2.4º/100ft inicia a una profundidad de 7861 ft ya que cuenta con un desplazamiento horizontal hasta
el objetivo de 6158 ft
reH =1153 ft KHg=18.15 md rw=2 ½” Twh=68ºF RGC=784952 PC/Bbl SoH=0.24
CALCULAMOS ( L)
𝐷2 − 𝐷1 180 1 180 1 180 1 180 1
sin 𝜃 = 𝑑𝑜𝑛𝑑𝑒 𝑟1 = ∗ = ∗ = 3183 𝑓𝑡 𝑦 𝑟2 = ∗ = ∗ = 2387 𝑓𝑡
𝑟1 𝜋 𝑞 𝜋 1.8 𝜋 𝑞 𝜋 2.4
100 𝑓𝑡 100 𝑓𝑡
6574 − 3643
𝜃 = sin−1 ( ) = 67º 𝑑𝑜𝑛𝑑𝑒 𝜃 + 𝛽 = 90º 𝑝𝑜𝑟 𝑙𝑜 𝑡𝑎𝑛𝑡𝑜 𝛽 = 23º
3183
𝑋 = 𝑟1 ∗ (1 − 𝑐𝑜𝑠𝜃) + tan(𝜃) ∗ (𝐷3 − 𝐷2) + 𝑟2 ∗ sin 𝛽 + 𝐿
𝐿 = 6158 − 3183 ∗ (1 − 𝑐𝑜𝑠67º) − tan(67º) ∗ (7861 − 6574) − 2387 ∗ sin 23º = 254 𝑓𝑡
Calculamos TVD
𝑇𝑉𝐷 = 𝐷3 + 𝑟2 ∗ (1 − 𝑐𝑜𝑠𝛽) = 7861 + 2387(1 − 𝑐𝑜𝑠23º) = 8051 𝑓𝑡
Calcular la porosidad del Cf = 2.7121E-6 psi-1 para próximamente calcular K absoluta y después Kg
𝟏. 𝟕𝟖𝟐
𝑪𝒇 = ( 𝟎,𝟒𝟑𝟖 ) ∗ 𝟏𝟎−𝟔 𝑷𝒔𝒊−𝟏 𝒅𝒐𝒏𝒅𝒆 𝒔𝒆 𝒅𝒆𝒔𝒑𝒆𝒋𝒂 ∅ = 𝟎. 𝟑𝟖𝟑𝟑 = 𝟑𝟖. 𝟑𝟑%
∅
Calcular la permeabilidad absoluta con la porosidad y Swc=25%
𝟐 𝟐
∅𝟑 𝟎. 𝟑𝟖𝟑𝟑𝟑
𝑲 = 𝟐, 𝟓 ∗ ( ) ∗ 𝟏𝟎𝟎𝟎 ( 𝒎𝒅) = 𝑲 = 𝟐, 𝟓 ∗ ( ) ∗ 𝟏𝟎𝟎𝟎 = 𝟖𝟐. 𝟓 𝒎𝒅
𝑺𝑾𝑪 𝟎. 𝟑𝟏
𝑲𝒐
Calcular Ko donde 𝑲𝒓𝒐 = 𝒚 𝑲𝒐 = 𝑲 ∗ 𝑲𝒓𝒐 = 𝟑𝟎. 𝟓𝟐𝟓 𝒎𝒅
𝑲
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
𝐾3 = −0.29478 𝐾7 = 1.90633
𝐾4 = −0.04783 𝐾8 = 8.41763
𝑃𝑑 = 2331 𝑃𝑠𝑖
0.5
4 0.5
𝐿 𝑟𝑒𝐻 𝐿
𝑎= 0.5 + [0.25 + ( ) ] 𝑝𝑎𝑟𝑎 < 0.9𝑟𝑒𝐻
2 𝐿 2
{ 2 }
0.5
4 0.5
254 1153 75 𝑚𝑑
𝑎= {0.5 + [0.25 + ( 254 ) ] } = 1156.5 𝑓𝑡 𝑦 𝐼𝑎𝑛𝑖 = √ = 3.1622
2 7.5 𝑚𝑑
2
Qg MMPCD Qg MMPCD (acumulado)
0 0
91.24279534 91.24279534
201.203193 201.203193
331.0532046 331.0532046
466.4869399 466.4869399
132.3529942 598.8399341
158.4307541 757.2706882
IPR
8000
7000
6000
5000
4000
3000
2000
1000
0
0 200 400 600 800
POXO X1
𝟐𝟕 ∗ 𝟖𝟏𝟎𝟓
𝑻𝒓 = º𝐆 ∗ 𝐓𝐕𝐃 = + 𝟕𝟎 = 𝟐𝟖𝟖. 𝟖𝟑º𝑭 = 𝟕𝟒𝟖. 𝟓 º𝑹
𝟏𝟎𝟎𝟎
𝑹𝑮𝑪𝒑𝒓𝒐𝒎𝒆𝒅𝒊𝒐 = 𝟖𝟏𝟐𝟎𝟑𝟐 𝒑𝒄/𝑩𝒃𝒍
Calculamos la presión de roció CON Tr, RGC promedio y el API=65º
Calcular la Presión de rocío con temperatura promedio
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
𝐾3 = −0.29478 𝐾7 = 1.90633
𝐾4 = −0.04783 𝐾8 = 8.41763
𝑃𝑑 = 2278.3 𝑃𝑠𝑖
Pfw (psi) ESTADO Ppc Tpc Z Dens (gr/cm3) u (cp) P/UZ MP
(MMpsi2/Cp)
7245 MONOFASICO 669.496 386.29 1.172898 0.246220 0.03048 202654.296 2170.06
6037.5 MONOFASICO 669.496 386.29 1.077049 0.223443 0.02786 201193.507 1682.41
4830 MONOFASICO 669.496 386.29 0.996012 0.193298 0.02488 194945.584 1204.08
3622.5 MONOFASICO 669.496 386.29 0.941331 0.153395 0.02164 177855.883 753.924
2415 MONOFASICO 669.496 386.29 0.92461 0.104112 0.01857 140664.424 369.310
1207.5 BIFASICO 664.409 375.645 0.895543 0.053910 0.01633 82591.5726 99.7293
0 BIFASICO 664.409 375.645 0.840624 0 0.015 0 0
Presión (Pwh)
Periodo de flujo Presión (Pfw) psi Q gas (MMPCD) Mp MMpsi2/cp Delta Mp DELTA MP/Qg
psi
Flujo 1 7056 3875 6 2231.67168 100.85607 16.809345
Flujo 2 6678 3737 9 2036.38932 296.13843 32.90427
Flujo 3 6123 3372 14 1765.19967 567.32808 40.52343429
Flujo 4 5756 2973 17 1596.02368 736.50407 43.32376882
Flujo estabilizado 6222 3538 13.5 1812.21972 520.30803 38.5413324
Pfw vs Mp
2500
y = 3E-05x2 + 0,1046x
2000
MP (MMPSI2)
1500
1000
500
0
0 1000 2000 3000 4000 5000 6000 7000 8000
Pfw (psi)
Delta MP vs Qg
50
45 y = 2,2586x + 7,4164
40
35
30
Delta MP
25
20
15
10
5
0
0 2 4 6 8 10 12 14 16 18
Qg MMPCD
(B) 2.2586
(A) tranciente 7.4164
(A) estabilizado 8.050235556
IPR tranciente
8000
6000
Pfw
4000
2000
0
0 5 10 15 20 25 30 35
Qg MMPCD
Para una presión de fondo fluyente de 3400 psi calculamos el caudal del gas mediante la ecuación de
segundo grado o mediante la gráfica IPR en el que es de 24.75 MMPCD en el estado tranciente y 24.61
MMPCD en el estado estabilizado con el cual se recurre hacer los cálculos para determinar los caudales
críticos y hacer la comparación
141.5
𝛾𝑜 = = 0.72010
131.5 + 65
= 𝟑. 𝟒𝟗𝟎𝟓 𝒇𝒕/𝒔
Calculamos Z con Pwh y Twh donde Z se encuentra en estado bifásico ya que Pd>Pwh Z=0.710745
Lo cual señala que se está generando acumulación de líquido en el fondo del pozo en el
cual el caudal no es lo suficiente para desplazar el agua del fondo pero es suficiente para
desplazar el condenado lo cual señala que la producción no está de manera óptima.
EJEMPLO
Un reservorio de gas cuenta con un pozo direccional donde la trayectoria del pozo consiste en una sección
de incremento el cual se da a un ritmo de 1.1º por cada cien pies seguida de una sección donde el ángulo
se mantiene constante generando una sección vertical de 132,6m que alcanza hasta la profundidad de
3642,36m en el cual se registra un desplazamiento horizontal en ese punto de 1600.2 metros ya que a
partir de la profundidad de 3642,36 comienza a disminuir a un ritmo de 1.55º por cada cien pies hasta
formar para llegar al objetivo un tramo vertical de 345,6 m ya que cuando se realizaba la perforación
antes de llegar al objetivo a una profundidad de 4495.5 m se registraba una longitud medida de 5176,7
metros.
se realiza una prueba isocronal con tiempo variable donde la presión estática de reservorio es de 7840 psi
con un gradiente geotérmico de 0.0349ºF/m en el cual se registra una gravedad especifica del gas de 0.67
con un API de 46º. Considerando que la prueba de producción se realizó cuando inicialmente se registraba
una presión de 6183 psi y una temperatura de 71ºF en cabeza donde el pozo cuenta con una tubería de
producción de diámetro de 3” con una rugosidad de 0,0026”
a) Determinar el AOF en estado tranciente y estabilizado
b) Calcular el caudal óptimo de producción para una presión de cabeza de 3500 psi
c) Realizar el análisis de caudales críticos para una presión de cabeza de 3500 psi y una densidad del
agua
DATOS DE LA PRUEBA ISOCRONAL
de
Periodo de Tiempo CHOKE Presión Presión Q gas Qcond Qagua
flujo (hr) (CK/36") (Pfw) psi (Pwh) psi (MPCD) (BPD) (BPD)
Flujo 1 14 24 6720 3183 3200 144 16
Flujo 2 23 28 5845 2789 4560 163 26
Flujo 3 28 32 4651 2234 5980 198 38
Flujo 4 37 36 3263 1789 7100 216 40
Flujo
estabilizado 62 32 5926 2876 4103 165 27
1,08925 gramos por centímetro cubico.
SOLUCION
𝐷3 − 𝐷2 = 3716 𝑓𝑡 𝑑𝑜𝑛𝑑𝑒 𝑠𝑒 𝑜𝑏𝑡𝑖𝑒𝑛𝑒 𝑞𝑢𝑒 𝐷2 = 𝐷3 − 3716 = 11950 − 3716 = 8234 𝑓𝑡
180 1
𝑟1 = ∗ = 5208.7 𝑓𝑡 𝑑𝑜𝑛𝑑𝑒 𝑋2 = 𝑟1 ∗ (1 − cos 𝜃) = 4092.55 ∗ (1 − cos 𝜃)
𝜋 1.1º
100 𝑓𝑡
180 1
𝑟2 = ∗ = 3696.5 𝑓𝑡 𝑑𝑜𝑛𝑑𝑒 𝑋3 = 𝑋2 + tan 𝜃 ∗ (𝐷3 − 𝐷2)
𝜋 1.55º
100 𝑓𝑡
𝑋3 = 𝑟1 ∗ (1 − cos 𝜃) + tan 𝜃 ∗ (𝐷3 − 𝐷2) = 5208.7 ∗ (1 − cos 𝜃) + tan 𝜃 ∗ (11950 − 8234)
𝜽 = 𝟒𝟓. 𝟎𝟎º
𝑇𝑉𝐷 = 𝐷3 + 𝑟2 ∗ sin 𝜃 + 1134 = 11950 + 3696.5 ∗ sin 45º + 1134 = 15697.8 𝑓𝑡
𝐺𝐶𝑅𝐾2 𝐾5 −𝐾6∗%𝐶7𝐾7 )
𝑃𝑑 = 𝐾1 [ ∗ 𝐾8 ∗ 𝐴𝑃𝐼 (𝐾4∗𝑇𝑟 ] Tr (ºF) GCR (Pc/bbl)
%𝐶7𝐾3
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
𝐾3 = −0.29478 𝐾7 = 1.90633
𝐾4 = −0.04783 𝐾8 = 8.41763
%𝐶7+= 10260 ∗ [(𝑅𝐺𝐶 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 ∗ 𝛿𝑔)−0.8499 ]
72514.6𝐾2 𝐾5 −𝐾6∗%𝐶7𝐾7 )
𝑃𝑑 = 𝐾1 [ ∗ 𝐾8 ∗ 58(𝐾4∗220 ] = 4169 Psi Tr (ºF) GCR (Pc/bbl)
%𝐶7𝐾3
Calculamos la Ppc y Tpc para el estado monofásico
𝑃1 𝑃𝑜
+ 0+0
𝜇1𝑧1 𝜇𝑜𝑧𝑜
𝑚(𝑃11) = 2 [ ] ∗ (𝑃1 − 𝑃0) + 0 = 2 [ ]∗ (0 − 0) + 0 = 0 (MMPsi2/cp)
2 2
𝑃2 𝑃1 69258.7+0
+ 2[ ]∗(784−0)
𝜇2𝑧2 𝜇1𝑧1 2
𝑚(𝑃10) = 2 [ ] ∗ (𝑃2 − 𝑃1) + 𝑚(𝑃11) = + 0 = 54.299 (MMPsi2/cp)
2 1000000
𝑃1 𝑃𝑜 116464.2+69258.7
+
𝜇1𝑧1 𝜇𝑜𝑧𝑜 2[ ]∗(1568−784)
2
𝑚(𝑃9) = 2 [ ] ∗ (𝑃1 − 𝑃0) + 𝑚(𝑃10) = + 54.299 = 199.906 (MMPsi2/cp)
2 1000000
𝑃2 𝑃1 146240.3−116464.2
+ 2[ ]∗(2352−1568)
𝜇2𝑧2 𝜇1𝑧1 2
𝑚(𝑃8) = 2 [ ] ∗ (𝑃2 − 𝑃1) + 𝑚(𝑃8) = + 199.906 = 405.866 (MMPsi2/cp)
2 1000000
REALIZAMOS EL AJUSTE DE GRAFICA
Pfw vs MP
3000,000
2000,000
MP MMpsi/cp
1500,000
1000,000
500,000
0,000
0 1000 2000 3000 4000 5000 6000 7000 8000 9000
Pfw psi
4 ∗ 𝑎0 + 𝑎1 ∗ ∑ 𝑋 + 𝑎2 ∗ ∑ 𝑋2 = ∑ 𝑌 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 1
𝑎0 ∗ ∑ 𝑋 + 𝑎1 ∗ ∑ 𝑋2 + 𝑎2 ∗ ∑ 𝑋3 = ∑ 𝑋 ∗ 𝑌 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 2
𝑎0 ∗ ∑ 𝑋2 + 𝑎1 ∗ ∑ 𝑋3 + 𝑎2 ∗ ∑ 𝑋4 = ∑ 𝑋 2 ∗ 𝑌 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 3
𝑎1 ∗ ∑ 𝑋2 + 𝑎2 ∗ ∑ 𝑋3 = ∑ 𝑋 ∗ 𝑌 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 2
𝑎1 ∗ ∑ 𝑋3 + 𝑎2 ∗ ∑ 𝑋4 = ∑ 𝑋2 ∗ 𝑌 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 3
𝒂𝟏 = 𝟎. 𝟏𝟑𝟒𝟓𝟓𝟐 𝑦 𝒂𝟐 = 𝟐. 𝟑𝟖𝟒𝟓𝑬 − 𝟓
Pfw vs MP
3000,000
2500,000 y = 2E-05x2 + 0.1378x
MP MMpsi/cp
2000,000
1500,000
1000,000
500,000
0,000
0 2000 4000 6000 8000 10000
Pfw psi
MP recalculado
2520.534912
2136.57317
1781.924372
1456.58852
1160.565612
893.855648
656.4586291
448.3745549
269.6034253
120.1452403
0
FLUJO Pfw (psi) Qg (MMPCD) MP Mmpsi2/cp Delt MP DELTMP/Q
1 6720 3.2 1980.991488 539.543424 168.60732
2 5845 4.56 1601.097616 919.4372959 201.6309859
3 4651 5.98 1141.611647 1378.923265 230.5891748
4 3263 7.1 692.9249208 1827.609991 257.4098579
Extendido 5926 4 1634.731367 885.8035448 221.4508862
DELTMP/Q vs Qg
300
250
y = 22.497x + 97.349
DELTMP/Qg
200
150
100
50
0
0 1 2 3 4 5 6 7 8
Qg
A tranciente=97.349
B =22.479
Para el estado estabilizado
DELTMP/Q vs Qg
300
250
DELTMP/Qg
200
150
100
50
0
0 1 2 3 4 5 6 7 8
Qg
A estabilizado=131.53
B =22.479
Calculamos los caudales de gas Qg
−𝐴 + √𝐴 2 + 4𝐵(𝑀𝑃𝑝𝑟 − 𝑀𝑃𝑝𝑓𝑤)
𝑄𝑔 =
2𝐵
Pfw MP recalculado Qg MMPC tranc Qg MMPC estabilizado
7840 2520.534912 0 0
7056 2136.57317 2.500455257 2.137946593
6272 1781.924372 3.962184981 3.509930405
5488 1456.58852 5.047111651 4.550272404
4704 1160.565612 5.90858947 5.384473539
3920 893.855648 6.612655115 6.070065576
3136 656.4586291 7.194897329 6.639058431
2352 448.3745549 7.6769601 7.111301536
1568 269.6034253 8.073003796 7.499950368
784 120.1452403 8.392712864 7.814084014
0 0 8.642863419 8.060096497
IPR tranciente
9000
8000
7000
6000
5000
Pfw
4000
3000
2000
1000
0
0 1 2 3 4 5 6 7 8 9 10
Qg
IPR ESTABILIZADO
10000
8000
6000
Pfw
4000
2000
0
0 1 2 3 4 5 6 7 8 9
Qg
IPR VLP
10000
8000
Pfw psi
6000
4000
2000
0
0 1 2 3 4 5 6 7 8 9
Qg MMPCD
Series1 Series2
Realizamos el análisis de caudales críticos para el caudal de 13
MMPCD del pozo X1
Densidad del agua 1,08925g/cm3=68 lPC
ρo = 68 ∗ 𝛾𝑜 = 68 ∗ 0.79718 = 54,208 lpc
Calculamos Z con Pwh y Twh donde Z se encuentra en estado bifásico ya que Pd>Pwh Z=0.798196
El caudal del gas no es suficiente para remover el liquido del fondo del pozo hasta superficie con lo cual se
considera que se requiere de un levantamiento artificial.
EJEMPLO DE FLUJO ESTATICO
Un reservorio de gas el cual cuenta con un pozo direccional de tipo slant donde se tiene un desplazamiento
horizontal hasta el objetivo de 4780 ft donde se tiene un ritmo de incremento angular de 2º/100ft ya que
se determina una longitud medida de 8579 ft hasta el punto donde el ángulo comienza a mantenerse
constante el cual se encuentra a una profundidad de 7980 ft.
El pozo produce a un caudal de 8.3 MMPCD de gas con una gravedad específica de 0.66, también produce
396 BPD de condensado que tiene un API de 40º y 228 BPD de agua ya que el pozo se encuentra con una
presión de fondo fluyente de 3285 psi donde la temperatura de reservorio es de 275 ºF :
a) Calcular la presión de cabeza a una temperatura de 63ºF.
b) Optimizar la producción del pozo mediante caudales críticos para un diámetro de 4”.
c) Calcular la presión de fondo fluyente considerando una presión de cabeza de 5000 psi y una
temperatura de 63ºF. SOLUCION
Pwh=? Twh=63ºF
D1
D2=7980ft
Kop r 𝜃 q=2º/100ft
D3 r
X2 𝜃
SOLUCION:
180 1
𝑟= ∗ = 2864.79 𝑓𝑡 𝑑𝑜𝑛𝑑𝑒 𝑋2 = 𝑟 ∗ (1 − cos 𝜃) = 2864.79 ∗ (1 − cos 𝜃)
𝜋 2º
100 𝑓𝑡
𝜃
𝐷1 + = 8579 𝑓𝑡 𝑑𝑜𝑛𝑑𝑒 𝐷2 = 𝐷1 + 𝑟 ∗ sin 𝜃 → 7980 = 𝐷1 + 2864.79 ∗ sin 𝜃
2º
100𝑓𝑡
𝜃
𝐶𝑜𝑚𝑏𝑖𝑛𝑎𝑛𝑑𝑜 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛𝑒𝑠 𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟𝑒𝑠 7980 − 2864.79 ∗ sin 𝜃 + = 8579 𝑓𝑡 → 𝜃 = 63º
2º
100𝑓𝑡
𝑋3 = 𝑋2 + tan 𝜃 ∗ (𝐷3 − 𝐷2) → 4780 = 2864.79 ∗ (1 − cos 63º) + tan 63º ∗ (𝐷3 − 7980) → 𝑫𝟑 = 𝟗𝟔𝟏𝟖. 𝟓𝟐 𝒇𝒕
(𝐷3 − 𝐷2) (9618.52 − 7980)𝑓𝑡
𝑀𝐷 = 8579 𝑓𝑡 + → 𝑀𝐷 = 8579 𝑓𝑡 + = 𝟏𝟐𝟏𝟖𝟖. 𝟐 𝒇𝒕
cos 𝜃 cos 63º
Cálculos adicionales
141.5 141.5 8.3𝑀𝑀𝑃𝐶𝐷
𝛾𝑜 = = = 0.825072 𝑑𝑜𝑛𝑑𝑒 𝑅𝐺𝐶 = = 20959.6 𝑝𝑐/𝐵𝑏𝑙
131.5 + °𝐴𝑃𝐼 131.5 + 40 396 𝐵𝑃𝐷
4591 ∗ 𝛾𝑜 4591 ∗ 0.825072
𝛿𝑔 + 0.66 +
𝛿𝑚 = 𝐺𝐶𝑅 = 20959.6 = 0.797969
1123 1123
1+ 1+
𝐺𝐶𝑅 20959.6
Calcular la Presión de rocío
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
𝐾3 = −0.29478 𝐾7 = 1.90633
𝐾4 = −0.04783 𝐾8 = 8.41763
𝑃𝑑 = 3282.5 𝑃𝑠𝑖
a) Método de presión y temperatura promedio:
𝑃𝑤𝑓 3285 𝑝𝑠𝑖
𝑃𝑤ℎ∗ = = = 2648.2 𝑝𝑠𝑖
(1 + 2.5𝐸 − 5 ∗ 𝑇𝑉𝐷) (1 + 2.5𝐸 − 5 ∗ 9618.52 𝑓𝑡)
𝑃𝑤𝑠 + 𝑷𝒘𝒉 ∗ 3285 + 2648.2 275 + 63
𝑃̅ = = = 2966.6 𝑝𝑠𝑖 𝑇̅ = = 169º𝐹 = 629º𝑅
2 2 2
Calculamos Z con Pwh y Twh donde Z se encuentra en estado bifásico ya que Pd>Pwh Z=0.680884
El reservorio produce gas a un caudal de 6,3MMPCD con una presión de fondo fluyente de 2250 psi en el cual la
gravedad especifica del gas es de 0,65 con un API de 45 o donde la relación gas condensado es de 123000PC/Bbl en el
cual el diámetro es de 3.2” con una rugosidad de 0,0026” ya que se registra una presión en el separador de 460 psi
Calcular:
Solución:
180 1
Datos 𝑟1 = ∗ 1,2 = 4774,65 𝑓𝑡
𝜋
100
q=1,2o/100ft 𝜃 = 51,6
51,6
Lbc=1865 ft 𝐿𝑎𝑏 = 1,2 = 4300 𝑓𝑡
100
D2=5764 ft
𝐷1 = 5764 − 1865 ∗ cos(51,6) = 2022,138 𝑓𝑡
a)
Datos :
TVD=6922,44 ft
MD=8187,14 ft
Go=39oc/1000m
API= 45
D=3,2 “
E=0,0026”
SG=0,65
METODO P Y T
ESTADO ESTATICO
Algoritmo
1 suponer Pwh#
5 calcular Pwh:
0,01875∗𝑇𝑉𝐷∗𝑆𝐺
( )
𝑃𝑤ℎ = 𝑃𝑤ℎ ∗ 𝑒 𝑇∗𝑍̅
1
2250
𝑃𝑤ℎ = = 1918 𝑝𝑠𝑖
(1 + 2,5 ∗ 10−5 ∗ 6922,44)
2
2250 + 1918
𝑝̅ = = 2084 𝑝𝑠𝑖
2
3
220,11 + 72
𝑇̅ = = 146,1( 𝐹) = 605,8(𝑅)
2
4 Calculo de Z
Z BIFASICO
Sg 0,65 RGC 123000
T®R 605,8
Api 45
cálculos
Ppc 666,902 Sgo 0,80169972 Tpr 1,63248969
Tpc 371,090 Sgm 0,67377201
Pwf Ppr Z
2084 3,12489888 0,77869434
5
2250
𝑃𝑤ℎ = 0,01875∗0,67377∗6922,44
= 1869,3 𝑝𝑠𝑖
[ ]
𝑒 605,8∗0,778694
1869,3 ≠ 1918
Volver al paso uno
𝑍̅ Pwh
Pwh# 𝑃̅ calculado
1918 2084 0,77869 1869,3
1869,3 2059,7 0,77692 1868,7
1868,7 2059,35 0,7769 1868,5
1868,5 2059,25 0,7768954 1868,5
ESTADO DINAMICO
Este método también tiene un algoritmo de cálculo:
Algoritmo
25 ∗ 𝑄𝑔2 ∗ 𝑇̅ ∗ 𝑍̅ ∗ 𝑓 ∗ 𝑈𝑔 ∗ 𝑠𝑔 ∗ (𝑒 𝑠 − 1)
𝐺𝑟𝑢𝑝𝑜(1) =
𝑠 ∗ 𝐷5
6 calcular Pwh
𝑃𝑤𝑓 2 − 𝑔𝑟𝑢𝑝𝑜(1)
𝑃𝑤ℎ = √
𝑒𝑠
1 calcular
2250
𝑃𝑤ℎ# = = 1918 𝑝𝑠𝑖
(1 + 2,5 ∗ 10−5 ∗ 6922,44)
2250 + 1918
𝑝̅ = = 2084 𝑝𝑠𝑖
2
220,11 + 72
𝑇̅ = = 146,1( 𝐹) = 605,8(𝑅)
2
𝑍̅ = 0,7787
2
0,0375 ∗ 0,674 ∗ 6922,44
𝑆= = 0,37073
605,8 ∗ 0,778694
3
20 ∗ 6,3 ∗ 1000 ∗ 0,674
𝑁𝑅𝐸 = = 1526818
0,01737 ∗ 3,2
4
0,0026 21,25
𝑓 = 4 ∗ [2,28 − 4 ∗ log ( + )]−2 = 0,01897
3,2 15268180,5
5
22502 − 176686,45
𝑃𝑤ℎ = √ = 1836,42 𝑝𝑠𝑖
𝑒 0,37073
Ecuaciones:
1
5,3 ∗ (𝜌𝐻2𝑂 − 0,00279 ∗ 𝑃ℎ𝑤)ˆ ( )
𝑉𝑔𝑎𝑠 − 𝑎𝑔𝑢𝑎 = 4
1
(0,00279 ∗ 𝑃𝑤ℎ)2
1
4,03 ∗ (𝜌(𝑐𝑜𝑛𝑑𝑒𝑠𝑎𝑑𝑜) − 0,00279 ∗ 𝑃ℎ𝑤)ˆ ( )
𝑔𝑎𝑠 𝑐𝑜𝑛𝑑𝑒𝑛𝑠𝑎𝑑𝑜 = 4
1
(0,00279 ∗ 𝑃𝑤ℎ)2
3,06 ∗ 𝑉𝑔 ∗ 𝐴 ∗ 𝑃𝑤ℎ
𝑄=
𝑇∗𝑍
Resolución:
Estado dinámico
Datos:
Pwh=1835,1306 psi
P H2O=64,5
Z= 0,77568396
D=3.2”
T=605,8 oR
API=45
1
5,3 ∗ (64,5 − 0,00279 ∗ 1835,1306)ˆ ( )
𝑉𝑔𝑎𝑠 − 𝑎𝑔𝑢𝑎 = 4 = 6,502
1
(0,00279 ∗ 1835,1306)2
141,5
𝛾𝑜 = = 0,8017
131,5 +
𝜌𝑜 = 𝑦𝑜 ∗ 64,5 = 51,71
1
4,03 ∗ (𝜌𝑜 − 0,00279 ∗ 1835,1306)ˆ ( )
𝑔𝑎𝑠 𝑐𝑜𝑛𝑑𝑒𝑛𝑠𝑎𝑑𝑜 = 4 = 4,653
1
(0,00279 ∗ 1835,1306)2
𝜋
𝐴= ∗ 3,22 = 8,0424 = 0,05585 𝑓𝑡^2
4
3,06 ∗ 𝑉𝑔 ∗ 𝐴 ∗ 1835,1306
𝑄 𝑔𝑎𝑠 − 𝑎𝑔𝑢𝑎 = = 4,339 𝑀𝑀𝑃𝐶𝐷
605,8 ∗ 0,77568396
3,06 ∗ 𝑉𝑔 ∗ 𝐴 ∗ 1835,1306
𝑄 𝑔𝑎𝑠 𝑐𝑜𝑛𝑑𝑒𝑛𝑠𝑎𝑑𝑜 = = 3,105 𝑀𝑀𝑃𝐶𝐷
605,8 ∗ 0,77568396
Comparando el Qg>Qgas-agua>Qgas-condensado con lo cual el pozo produce de manera optima
EJEMPLO
Un reservorio que costa de dos formaciones productoras que producen gas mediante dos pozos donde cada uno
atraviesa una formación. Ya que se lleva a cabo una prueba de flujo tras flujo en la formación de gas ICLA que
produce gas mediante el pozo X1 con el cual se cuenta con los siguientes datos:
La formación de gas Santa Rosa el cual atraviesa el pozo X2 con un diámetro de 4 in ya que cuenta con un factor de
daño de 6 donde se registra una porosidad de 32% a una profundidad de 7940 pies con una temperatura de 256 F y
una presión estática de 5760 psi donde la permeabilidad absoluta es de 14 md y la permeabilidad efectiva del agua
es igual a la del condensado que equivale a 6 md en el cual se determinara un espaciamiento de 78 acres del
reservorio con un espesor neto de 36 ft de la formación productora del cual se obtuvo una gravedad especifica del
gas de 0.66 y un gravedad especifica del condensado de 0.7796 en el cual indican que la relación gas condensado es
de 116789 pc por barril. Realizar el IPR compuesto por el método de LIT
Datos recompilados
Calculo de RGC
𝐐𝐠
RGC =
𝐐 𝐜𝐨𝐧𝐝𝐞𝐧𝐬𝐚𝐝𝐨
Q condensado
Qg(MMPCD)
(BBP) RGC
4.465 64 69765.625
5.78 73 79178.0822
6.51 78 83461.5385
7.26 81 89629.6296
6.45 71 90845.0704
promedio 82575.99
Calculo de SGo
141.5
SGo = = 0.759
131.5 + 55
Calculo de SGmix
4591 ∗ 0.759
0.65 +
SGmix = 82575.99 = 0.683
1123
1+
82575.99
Calculo la precion de roció
182575.99 −0.8207
%C7+= ( ) = 0.88
70680
82575.990.09741 0.28126 −0.00068∗0.881.90633 )
Pd = 346.77648 ∗ [ ∗ 8.41763 ∗ 55(−0.04783∗245 ] = 3433.9
0.88−0.29478
Calculo de Propiedades Pseudocriticas:
Tpc 371.291 °R
Tpc 373.388 °R
Calculo de Z para una Pfw=3750 psi (bifásico) y 7500 psi (monofásico) y Tr=705ºR
Z monofásico correlacion de bexanbrill
1.383
μ gas = 10−4 ∗ 143.8997 ∗ e5.0869∗0.2507 = 0.0305
Pwf ρ
Fase Ppr B Z μ gas(cp)
[Psi] gas(gr/cm3)
Po P1
𝑀𝑃 = ( + ) + (𝑃𝑜 − 𝑃1)
Uo ∗ Zo U1 ∗ Z1
7500 6250
𝑀𝑃 = ( + ) + (7500 − 6250) = 5.15𝐸8
1.193 ∗ 0.0305 0.0279 ∗ 1.090
Calculo de MP acumulado
Pwf MP
Fase P/U*Z MP
[Psi] acumulado
Pwf Vs MP acumulado
2500,00
1000,00
500,00
0,00
0 1000 2000 3000 4000 5000 6000 7000 8000
a=3E-5
b=0.1192
MP=3E-5*pfw^2+0.1192*pfw
𝑀𝑃 2581.5 − 2059.6
= = 116.88
𝑄 4.465
Pwf
Fase Qg(MMPCD) MP acumulado
(psia) MP/Q
inicio 7500 0 2581.5 -
1 6534 4.465 2059.6 116.88
2 5213 5.78 1436.7 198.07
3 4756 6.51 1245.5 205.22
4 3720 7.26 858.6 237.32
extendido 4295 6.45 1065.4 235.06
Estado transiente
MP/Qg VS Qg
300
250 y = 41.827x - 61.745
200 R² = 0.944
150
100
50
0
0 1 2 3 4 5 6 7 8
A=-61.745
B=41.827
A=-34.72415
B=41.827
co
Qg transiente Qg estabilizado
(MMPCD) (MMPCD)
Datos recompilados
PM=28.97*0.66=19.1202
Calculo de la permeabilidad del gas
𝐾 = 𝐾𝑔 + 𝐾𝑜 + 𝐾𝑤
14 = 𝐾𝑔 + 6 + 6
Krg=2/14=0.1428
Kg=2md
43560 ∗ 78
𝑟𝑒 = √ = 1039.96 𝑓𝑡
𝜋
Rw=4 in=0.33 ft
Calculo de API
141.5
𝐴𝑃𝐼 = − 131.5 = 50
0.7796
Calculo de SGmix
4591 ∗ 0.7796
0.66 +
SGmix = 116789 = 0.6841
1123
1+
116789
Calculo de la presio de rocio
116789 −0.8207
%C7+= ( ) = 0.662
70680
1167890.09741 0.28126 −0.00068∗0.6621.90633 )
Pd = 346.77648 ∗ [ −0.29478 ∗ 8.41763 ∗ 50(−0.04783∗256 ] = 3303.65
0.662
Propiedades Pseudocriticas:
Tpc 373.655 °R
Tpc 373.661 °R
Calculo de Z
Pwf ρ
Fase Ppr B Z μ gas(cp)
[Psi] gas(gr/cm3)
5760 8.635 3.4493 1.055 0.218 0.0268
Monofásic
Pwf Mp M(P)
Fase P/U*Z MP
[Psi] acumulado [Mmpsi/cp]
5760 202530.3 3.85E+08 1.63E+09
Monofási
1616.6
4800 198057.7 3.70E+08 1.24E+09
co
1230.9
3840 186933.7 3.34E+08 8.75E+08 861.84
2880 160742.2 2.76E+08 5.40E+08 530.7
Bífasico
monofasico
8000
7000 Qg estabilizado… Qg… Q*
6000
5000
4000
3000
2000
1000
0
0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 0.018
Twh=65 F=525 R
Δ𝑃 ∗ = 𝑃𝑤𝑓 − 𝑃𝑤ℎ
Pwh=4300 psi
Calculo de Z
Tpc 371.291 °R
Ppr 6.440
Ppr=4300/667.705=6.440
Tpr=525/371.291=1.414
Z monofasico
Z = 0.8685
3.06 ∗ 4.119 ∗ 0.2673 ∗ 4300
𝑄𝑔 − 𝑎𝑔𝑢𝑎 = = 31.77 𝑀𝑀𝑃𝐶𝐷
525 ∗ 0.8685
3.06 ∗ 2.869 ∗ 0.2673 ∗ 4300
𝑄𝑔 − 𝑎𝑔𝑢𝑎 = = 22.13 𝑀𝑀𝑃𝐶𝐷
525 ∗ 0.8685
Calculo para el caudal a Pwf=5500
MP=3E-5*5500^2+0.1192*5500=1563.1
B=41.827
A=-34.72415 B=41.827
Hay acumulación de líquido en el fondo del pozo por lo cual la producción noes optima
EJEMPLO
Un reservorio de gas el cual cuenta con un pozo direccional de tipo slant donde se tiene una profundidad
de 7984 pies en el cual el Kop se encuentra a 2950 pies donde se registra una presión de 1290 psi ya que
la profundidad medida es de 9345 pies.
El pozo produce gas con una gravedad específica de 0.66 donde la relación de gas condensado es de
14500PC/Bbl donde el diámetro es de 2.6 in con una rugosidad de 0.003 in en el que se registra una
permeabilidad de gas de 8 md con un espesor neto de 48 pies considerando una presión estática de 2640
psi con una temperatura de reservorio de 250 ºF y una temperatura de 68ºF en superficie, calcular:
Tramo2
D1=2950 pies
Tramo 1
Qg=?
Pwf= ? Psi
Twf=250ºF
RESERVORIO Pr= 2641 Psi S=9 h=48 ft
Kg= 8 md
Re=1560 ft
Considerando el método de calender dinámico con el cual realizamos el próximo procedimiento de cálculo para
trabajar en el que lo denominamos tramo 1 el cual llega a ser desde el punto del KOP hasta el objetivo donde se
encuentra la Pfw con el cual elaboramos el siguiente algoritmo para determinar el caudal del gas.
Calculamos H para el tramo 1 el cual llega a ser 𝑇𝑉𝐷´ = 𝐻 = 7984 − 2950 = 5034
Como se puede observar la relación se gas condensado es mayor a los 120000 PC/Bbl con el cual se puede
trabajar para el cálculo de Z netamente en estado monofásico.
Algoritmo matemático a seguir
𝛼
9) Calcular 𝑃𝑚 = 𝑃𝑤ℎ + 𝐼 con el cual comparamos con el supuesto en el paso 5
𝑡𝑠 +𝐼𝑚
𝛼
13) Calcular Pfw para comparar con el Pfw* 𝑃𝑤𝑓 = 𝑃𝑚 + 𝐼
𝑤𝑓 +𝐼𝑚
703∗10−6∗𝐾∗ℎ(𝑃𝑟 2−𝑃𝑤𝑓2 )
14) Calcular mediante la ecuación de Darcy el Qg y comparar con el supuesto 𝑄𝑔 = 𝑟𝑒 (MPCD)
𝑇𝑟∗𝜇𝑔∗𝑧[𝑙𝑛( )−0.75+𝑆]
𝑟𝑤
𝑇𝑓𝑤−𝑇𝑤ℎ 250−68
Paso 3 calculamos 𝑇𝐾𝑂𝑃 = 𝑇𝑤ℎ + ( )∗ ℎ = 68 + ( ) ∗ 2950 = 135.246º𝐹 = 595º𝑅
𝐻 7984
1290
595∗0.85
Paso 4 Con Pkop y Tkop calcular Z para determinar la integral en el Kop 𝐼𝑘𝑜𝑝 = 1290 2 5035 = 501.4020
0.001∗( ) ∗ +0.00388
595∗0.85 6395
1365.73
652.5∗0.8986
Paso 8 Calcular 𝐼𝑚 = 1365.73 2 5035 = 543.52
0.001∗( ) ∗ +0.00388
652.5∗0.8986 6395
𝛼 79138.125
Paso 9 Calcular 𝑃𝑚 = 𝑃𝑤ℎ + 𝐼 = 1290 + 501+543.52 = 1365.7353 con el cual comparamos con el supuesto en el paso 5
𝑡𝑠 +𝐼𝑚
1436.043
710∗0.93069
Paso 12 Calcular la integral a condiciones de fondo 𝐼𝑤𝑓 = 1436.043 2 5035
= 582.2019
0.001∗( ) ∗ +0.00388
710∗0.93069 6395
𝛼 79138.125
Paso 13 Calcular Pfw para comparar con el Pfw* 𝑃𝑤𝑓 = 𝑃𝑚 + 𝐼 = 1365.73 + 582.20+543.52 = 1436.03
𝑤𝑓 +𝐼𝑚
Calculamos Rw =0.10833 ft
Paso 14 Calcular mediante la ecuación de Darcy el Qg y comparar con el supuesto
703∗10−6 ∗8∗48(26412 −14362 )
𝑄𝑔 𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 1560 = 7000.5 (MPCD)=7.0005 MMPCD
710∗0.0160839∗0.93069[𝑙𝑛(0.10833)−0.75+9]
Mediante el algoritmo respondemos al inciso a y C donde el caudal del gas es de 7 MMPCD Y la Pfw=1436 Psi
B) Para la presión de cabeza calculamos trabajando en el tramo dos que es desde el kop hasta la superficie donde se
puede considerar H=2950 ft y como Tfw=Tkop donde la Twh=528ºR donde TVD=MD
𝛼
6) Calcular 𝑃𝑚 = 𝑃𝐾𝑜𝑝 − con el cual comparamos con el supuesto en el paso 2
𝐼𝐾𝑂𝑃 +𝐼𝑚
7) Suponer Pwh*
8) Calcular Z del gas con Pfw* y Tr
𝑃𝑤ℎ∗
𝑇𝑟 ∗𝑍
9) Calcular la integral a condiciones de fondo 𝐼𝑤ℎ = 𝑃𝑤ℎ∗ 2 𝑇𝑉𝐷
0.001∗( ) ∗ +𝐹2
𝑇𝑟 ∗𝑍 𝑀𝐷
𝛼
10) Calcular Pwh para comparar con el Pwh* 𝑃𝑤ℎ = 𝑃𝑚 − 𝐼
𝑤ℎ +𝐼𝑚
Por lo tanto se puede ver que la Pwh es 1192 .17 Psi con el cual en función a los resultados
obtenidos podemos ver que se genera mayor perdida de presion en el medio poroso con el
cual se pierden mas de 1000 p
EJEMPLO
Un reservorio de gas el cual consta de dos pozos productores donde el primer pozo (X1) convencional se
encuentra a una profundidad de 167 pies mayor que el pozo (X2) ya que el pozo X2 se encuentra a 1456
ft sobre el nivel del mar donde la trayectoria del pozo consiste en una sección de incremento el cual comienza
a una profundidad de 7495 pies hasta una profundidad de 8821 pies seguida de una sección donde el ángulo
se mantiene constante hasta una profundidad de 13268 pies en el cual comienza a disminuir a un ritmo
determinado donde generar una sección vertical de 935 pies con el cual se obtiene un desplazamiento
horizontal hasta el objetivo de 4124 pies y un MD de 16479 ft ya que a una profundidad de 11567 ft del
pozo x2 se registra una presión de 3602 psi.
El pozo x2 produce gas con una gravedad especifica de 0.66 y un API de 50º ya que está conformado con una
tubería de producción de 2 ½ ” donde se tiene una rugosidad de 0.0028 in Y también se registra una presión
del separador de 1368 psi ya que la temperatura en cabeza del pozo es de 65ºF donde se considera un
gradiente geotérmico de 16ºF/1000 ft donde la densidad del agua es de 68 Lpc en ambos pozos ya que el
pozo X2 se planifico con cerrarlo si el caudal de gas de 8,5 MMPCD se encontrara muy próximo al del valor
critico donde a ese caudal se registra una presión de fondo de 3938 psi.
El pozo x1 el cual estaba cerrado durante una etapa de intervención de pozo cuando producía un caudal de
9MMPCD mediante el cual tras la intervención se llegó a incrementar el caudal del gas en un 44,44% donde
actualmente se registra un índice de productividad de 4,6803 PC/psi^2 a una RGC de 89456 PC/Bbl con una
gravedad especifica de 0.65 y un API de 48 se registra una presión de 2970 psi ya que el pozo cuenta con una
tubería de 3” de diámetro con una rugosidad de 0.0030 in donde el gradiente geotérmico es de 15,2ºF/1000
ft y a la temperatura de cabeza de 66 ºF
DATOS DEL RESERVORIO
a) Verificar mediante caudales críticos la decisión que se debe tomar para el pozo X2
b) Determinar la presión del pozo x2 a una profundidad de 13800 pies
c) Determinar la caída de presión desde el reservorio hasta la cabeza del pozo x1
d) Determinar cual es el pozo que cuenta con la producción más óptima y porque?
Considerar para el cálculo de caída de presión el método de cullender para el pozo X1 y Método de P y T
promedio para el pozo X2
SOLUCION
𝜃 𝐷3 − 𝐷2 𝜃
𝑀𝑉𝐷𝑋2 = 𝐷1 + + + + 935 → 𝐷2 = 𝐷1 + 𝑟1 ∗ sin 𝜃
𝑞1 cos 𝜃 𝑞2
𝜃 13268 − 8821 𝜃
16479 = 7495 + + + + 935
𝑞1 cos 𝜃 𝑞2
𝟑, 𝟓 𝟐, 𝟔
𝒒𝟐 = 𝒇𝒕 𝒒𝟏 = 𝒇𝒕 𝜽 = 𝟑𝟕º 𝒓𝟏 = 𝟐𝟐𝟎𝟑, 𝟕 𝒓𝟐 = 𝟏𝟔𝟑𝟕
𝟏𝟎𝟎 𝟏𝟎𝟎
Calculamos el TVD del pozo X1 y X2
b)
Comenzamos por el pozo X1 a calcular la Pfw en estado dinámico donde mediante el Índice de productividad de
obtendrá le presión de reservorio con el cual se trabajara con el caudal de gas actual que hace referencia a que se
incrementó un 44,44% donde el caudal actual es
0.10796∗𝑄𝑔 2 0.10796∗13 2
𝐹2 = ( ) =( ) = 0.006338 Para este caso como se parte desde cabeza hasta el fondo ya que
𝐷 2.612 32.612
el pozo es vertical TVD=MD donde TVDx1=H=15355 ft y donde h=H/2=7677,5 ft para una Twh=66ºF=525,67ºR y una
Twf=759ºR y para una Pwh=2970 psi
b) Caída de Presión desde el reservorio hasta la cabeza del pozo X1=4500-2970=1530 psi
a)
Del estado estabilizado calculamos analíticamente la Pwf para el pozo X2 a un caudal de 8,5 MMPCD
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
𝐾3 = −0.29478 𝐾7 = 1.90633
𝐾4 = −0.04783 𝐾8 = 8.41763
𝑃𝑑 = 3202,22 𝑃𝑠𝑖
T P Pwh
Pwh* promedio promedio ESTADO Z u (cp) Nre f S calculado
3100 646,17 3519 MONOFASICO 0,876128971 0,021731867 2065169,997 0,020361245 0,663989114 2507,757579
2507,757579 646,17 3222,878789 MONOFASICO 0,865429629 0,020698898 2168231,346 0,020353052 0,672198039 2500,209112
2500,209112 646,17 3219,104556 MONOFASICO 0,865316247 0,020685737 2169610,877 0,020352947 0,672286117 2500,128085
2500,128085 646,17 3219,064043 MONOFASICO 0,865315033 0,020685596 2169625,694 0,020352946 0,67228706 2500,127218
2500,127218 646,17 3219,063609 MONOFASICO 0,86531502 0,020685594 2169625,852 0,020352946 0,67228707 2500,127209
Pwh del pozo X2 Pwh=2500 psi
Realizamos el análisis de caudales críticos para el caudal de 8,5 del pozo X2
ρo = 62.4 ∗ 𝛾𝑜 = 62.4 ∗ 0.77961 = 48,6476 lpc
Calculamos Z con Pwh y Twh donde Z se encuentra en estado bifásico ya que Pd>Pwh Z=0.705296
Ya que el caudal del pozo x2 es mayor a los caudales críticos la producción es óptima 8,5 MMPCD>3,95MMPCD>2,72MMPCD por lo tanto no se procede a
cerrarlo el pozo x2
d)Realizamos el análisis de caudales críticos para el caudal de 13 MMPCD del pozo X1
ρo = 62.4 ∗ 𝛾𝑜 = 62.4 ∗ 0.7883 = 49,1899 lpc
5.3 ∗ (𝜌𝐿 − 0.00279 ∗ 𝑃𝑤ℎ )1/4 5.3 ∗ (68 − 0.00279 ∗ 2970)1/4
𝒗𝒈𝒂𝒔−𝒂𝒈𝒖𝒂 = = = 𝟓. 𝟏𝟏𝟖𝟏𝟔𝟐 𝒇𝒕/𝒔
(0.00279 ∗ 𝑃𝑤ℎ )1/2 (0.00279 ∗ 2970)1/2
Calculamos Z con Pwh y Twh donde Z se encuentra en estado bifásico ya que Pd>Pwh Z=0.752825
Ya que el caudal del pozo x1 es mayor a los caudales críticos la producción es óptima 13 MMPCD>6,322MMPCD>3,991MMPCD
Ya que el caudal del pozo x2 es mayor a los caudales críticos la producción es óptima 8,5 MMPCD>3,95MMPCD>2,72MMPCD
Analizando los resultados de los caudales críticos se puede considerar el pozo X1 como el pozo con la producción mas optima ya que el caudal de
producción se encuentra a un valor más alejado del caudal critico comparado con el del pozo X2
Considerar para el cálculo de caída de presión el método de cullender para el pozo X1 y Método de P y T
promedio para el pozo X2
SOLUCION
𝜃 𝐷3 − 𝐷2 𝜃
𝑀𝑉𝐷𝑋2 = 𝐷1 + + + + 935 → 𝐷2 = 𝐷1 + 𝑟1 ∗ sin 𝜃
𝑞1 cos 𝜃 𝑞2
𝜃 13268 − 8821 𝜃
16479 = 7495 + + + + 935
𝑞1 cos 𝜃 𝑞2
𝟑, 𝟓 𝟐, 𝟔
𝒒𝟐 = 𝒇𝒕 𝒒𝟏 = 𝒇𝒕 𝜽 = 𝟑𝟕º 𝒓𝟏 = 𝟐𝟐𝟎𝟑, 𝟕 𝒓𝟐 = 𝟏𝟔𝟑𝟕
𝟏𝟎𝟎 𝟏𝟎𝟎
Calculamos el TVD del pozo X1 y X2
Comenzamos por el pozo X1 a calcular la Pfw en estado dinámico donde mediante el Índice de productividad de
obtendrá le presión de reservorio con el cual se aplicará el método Jones and Glaze para trabajar en datos de prueba
0.10796∗𝑄𝑔 2 0.10796∗13 2
𝐹2 = ( ) =( ) = 0.006338 Para este caso como se parte desde cabeza hasta el fondo ya que
𝐷 2.612 32.612
el pozo es vertical TVD=MD donde TVDx1=H=15355 ft y donde h=H/2=7677,5 ft para una Twh=66ºF=525,67ºR y una
Twf=759ºR y para una Pwh=2970 psi
DP/Qg
1200000
1000000
y = 121005x - 536776
DP/Qg Psi^2/MMPCD
800000
600000
400000
200000
0
0 2 4 6 8 10 12 14
Qg MMPCD
IPR tranciente
Pfw Qg MMPCD
5000
4000
4500 0
Pwf psi
3000
4000 8,545869936
2000
3000 12,11197556
1000
2000 14,01678555
0
1000 15,02439272 0 5 10 15 20
0 15,34308334 Qg MMPCD
4500 0 4000
Pwf psi
3000 11,78226353 2000
2000 13,69639061
1000
1000 14,7077989
0
0 15,02756964
0 5 10 15 20
Qg MMPCD
Del estado estabilizado calculamos analíticamente la Pwf para el pozo X2 a un caudal de 8,5 MMPCD
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
𝐾3 = −0.29478 𝐾7 = 1.90633
𝐾4 = −0.04783 𝐾8 = 8.41763
𝑃𝑑 = 3202,22 𝑃𝑠𝑖
T P Pwh
Pwh* promedio promedio ESTADO Z u (cp) Nre f S calculado
3100 646,17 3519 MONOFASICO 0,876128971 0,021731867 2065169,997 0,020361245 0,663989114 2507,757579
2507,757579 646,17 3222,878789 MONOFASICO 0,865429629 0,020698898 2168231,346 0,020353052 0,672198039 2500,209112
2500,209112 646,17 3219,104556 MONOFASICO 0,865316247 0,020685737 2169610,877 0,020352947 0,672286117 2500,128085
2500,128085 646,17 3219,064043 MONOFASICO 0,865315033 0,020685596 2169625,694 0,020352946 0,67228706 2500,127218
2500,127218 646,17 3219,063609 MONOFASICO 0,86531502 0,020685594 2169625,852 0,020352946 0,67228707 2500,127209
Pwh del pozo X2 Pwh=2500 psi
Realizamos el análisis de caudales críticos para el caudal de 8,5 del pozo X2
ρo = 62.4 ∗ 𝛾𝑜 = 62.4 ∗ 0.77961 = 48,6476 lpc
Calculamos Z con Pwh y Twh donde Z se encuentra en estado bifásico ya que Pd>Pwh Z=0.705296
Ya que el caudal del pozo x2 es mayor a los caudales críticos la producción es óptima 8,5 MMPCD>3,95MMPCD>2,72MMPCD
Calculamos Z con Pwh y Twh donde Z se encuentra en estado bifásico ya que Pd>Pwh Z=0.752825
Ya que el caudal del pozo x1 es mayor a los caudales críticos la producción es óptima 13 MMPCD>6,322MMPCD>3,991MMPCD
Ya que el caudal del pozo x2 es mayor a los caudales críticos la producción es óptima 8,5 MMPCD>3,95MMPCD>2,72MMPCD
Analizando los resultados de los caudales críticos se puede considerar el pozo X1 como el pozo con la producción mas optima ya que
el caudal de producción se encuentra a un valor más alejado del caudal critico comparado con el del pozo X2
SOLUCION DEL 1ª PARCIAL DE EXPLOTACION GESTION 2017
COMPONETE % PORCENTAJE MOLAR PESO MOLECULAR
C1 85 16.01
C2 9 30.1
C3 4 44.1
n-C4 2 58.1
Donde el 𝑃𝑀𝑎𝑝𝑎𝑟𝑒𝑛𝑡𝑒 = 0.85 ∗ 16.01 + 0.09 ∗ 30.1 + 0.04 ∗ 44.1 + 0.02 ∗ 58.1 = 19.2435
Calculamos la RGP promedio de los 4 flujos incluido el extendido el cual se hará uso para
determinar la gravedad específica de la mezcla (Sgm) y la presión de Rocío (Pd).
62500 + 63636 + 70408 + 76141 + 89888
𝑅𝐺𝑃𝑃𝑟𝑚𝑒𝑑𝑖𝑜 = = 72514.6 𝑃𝐶/𝐵𝑏𝑙
5
Determinamos la gravedad específica del condensado con el API
141.5 141.5
𝛾𝑜 = = = 0.7467
131.5 + °𝐴𝑃𝐼 131.5 + 58
4591 ∗ 𝛾𝑜 4591 ∗ 0.7467
𝛿𝑔 +
𝐺𝐶𝑅𝑃𝑜𝑟𝑚𝑒𝑑𝑖𝑜 0.664256 + 72514.6
𝛿𝑚 = = = 0.7007
1123 1123
1+ 1+
𝐺𝐶𝑅𝑃𝑟𝑜𝑚𝑒𝑑𝑖𝑜 72514.6
𝐺𝐶𝑅𝐾2 𝐾5 −𝐾6∗%𝐶7𝐾7 )
𝑃𝑑 = 𝐾1 [ ∗ 𝐾8 ∗ 𝐴𝑃𝐼 (𝐾4∗𝑇𝑟 ] Tr (ºF) GCR (Pc/bbl)
%𝐶7𝐾3
Donde:
𝐾1 = 346.77648 𝐾5 = 0.28126
𝐾2 = 0.09741 𝐾6 = 0.00068
𝐾3 = −0.29478 𝐾7 = 1.90633
𝐾4 = −0.04783 𝐾8 = 8.41763
%𝐶7+= 10260 ∗ [(𝐺𝐶𝑅 ∗ 𝛿𝑔)−0.8499 ]
72514.6𝐾2 𝐾5 −𝐾6∗%𝐶7𝐾7 )
𝑃𝑑 = 𝐾1 [ ∗ 𝐾8 ∗ 58(𝐾4∗220 ] = 3647.9 Psi Tr (ºF) GCR (Pc/bbl)
%𝐶7𝐾3
Calculando el %C7 de la manera más correcta con Sg.
%𝐶7+= 10260 ∗ [(72514.6 ∗ 0.664256)−0.8499 ] = 1.0746
Aplicando de la segunda forma donde se genera más error la Pd=3551 Psi
𝑃2 𝑃1
+
𝜇2𝑧2 𝜇1𝑧1
𝑚(𝑃2) = 2 [ ] ∗ (𝑃2 − 𝑃1) + 𝑚(𝑃1) (Psi/cp)
2
MP vs Pfw
3000
2000
1500
1000
500
0
0 1000 2000 3000 4000 5000 6000 7000 8000 9000
-500
y = 2.2214E-05x2 + 0.145579x
Recalculamos los MP para cada flujo con la ecuación del hp y = 2.2214E-05x 2 + 0.145579x
Pfw Qg MP
FLUJO (psi) (MMPCD) Mmpsi/cp Delt MP DELTMP/Q
1 7784 4.5 2474.66406 80.99721872 17.99938194
2 6680 5.6 1960.407656 595.2536226 106.2952897
3 4546 6.9 1119.349942 1436.311336 208.1610632
4 3355 8.3 737.6259385 1818.03534 219.0404024
Extendido 6233 8 1767.539305 788.1219731 98.51524664
DELTMP/Q vs Qg
250
200
y = 37,712x - 121,16
150
100
50
0
0 1 2 3 4 5 6 7 8 9
A B
-121.16 37.712
Ya que existe otras formas de calcular A y B aplicamos de la siguiente forma:
∆𝑚(𝑃)
∑ ∑ 𝑞 2 − ∑ ∆𝑚(𝑃) ∑ 𝑞
𝑞
𝐴𝑡𝑟𝑎𝑛𝑐𝑖𝑒𝑛𝑡𝑒 = = −121.16
𝑁 ∑ 𝑞2 − ∑ 𝑞 ∑ 𝑞
∆𝑚(𝑃)
𝑁 ∑ ∆𝑚(𝑃) − ∑ ∑𝑞
𝑞
𝐵= = 37.712
𝑁 ∑ 𝑞2 − ∑ 𝑞 ∑ 𝑞
A estabilizado
-203.1807534
−𝐴 + √𝐴 2 + 4𝐵(𝑚(𝑃𝑟) − 𝑚(𝑃𝑤𝑓))
𝑄𝑔 =
2𝐵
Donde recalculamos los MP para cada presión de fondo fluyente con la ecuación:
(MMPCD)Qg (MMPCD) Qg
Pfw (MMPSI^2/Cp)MP tarncinet estabilizado
7948 2555.661279 0 0
6000 1670.514 6.710472835 8.237136321
5000 1281.395 7.637130199 9.100580328
4000 936.556 8.352784077 9.778354401
3000 635.997 8.919643591 10.32012369
2000 379.718 9.370368501 10.75336218
1000 167.719 9.724326012 11.09487739
0 0 9.99377617 11.35552559
8000
6000
4000
2000
0
0 2 4 6 8 10 12