Programación AvanzadaAplicaciones de Ingeniería Petrolera
Semestre 2019-1
Ing. Juan Carlos Sabido AlcántaraIngeniero Petrolero
Facultad de Ingeniería UNAM
Aplicaciones de Ingeniería Petrolera
• Dentro de la ingeniería petrolera en sus diferentes ramas esnecesario realizar una gran cantidad de cálculos que permitentomar decisiones sobre la forma en la que se explotará unyacimiento por medio de predicciones a futuro, el diseño de laperforación de los pozos que permitirán dicha explotación,conocer el tipo de fluidos producidos para así diseñar lasinstalaciones de producción y almacenaje ubicadas en lasuperficie, entre otros, todos estos cálculos pueden requerir demucho tiempo para obtener resultados, sin embargo con el usode una computadora esta inversión de tiempo puede reducirseconsiderablemente.
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Algunos cálculos al ser programados, requieren de una soluciónnumérica como las tratadas en esta asignatura, esto implica quese necesita de conocimientos avanzados de programación pararealizar un software que permita resolver este tipo deproblemas, sin embargo hay otros que no necesitan de estetipo de solución, pues los resultados se obtienen al aplicar losdatos de entrada indicados por el usuario a una ecuación queproporciona un valor que después se requerirá en otra, asíhasta que se obtenga el resultado final, éstos sólo requieren deconocimientos básicos de programación
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Newton – Raphson
• El método de Newton-Raphson es un método abierto,esto se dice porque no garantiza la convergencia global.Para obtener la convergencia del método se debe deseleccionar un valor inicial lo suficientemente cercano ala raíz que se busca. Lo anterior permite que la primeraiteración se haga con un valor muy cercano a cero,comúnmente se denomina a este como punto dearranque o valor supuesto.
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Newton – Raphson
• El punto inicial y su cercanía a la raíz va a depender de lanaturaleza de la función; si la misma presenta múltiplespuntos de inflexión o pendientes grandes cercanas a laraíz, la posibilidad de que el método diverja aumenta.Una vez que se ha seleccionado el punto inicial y seaplica el método, la función se linealiza por larecta tangente en ese valor supuesto. La abscisa alorigen de dicha recta es, de acuerdo al método, unamejor aproximación de la raíz que el valor anterior.
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Newton – Raphson
• Se realizan nuevas iteraciones tomando como nuevo punto elcalculado en la iteración anterior, hasta que el método converjahasta una cierta tolerancia.
• Tarea: Demostrar y obtener la formula de Newton – Raphson.Ing. Juan Carlos Sabido Alcántara
𝑥𝑖+1 = 𝑥𝑖 −𝑓 𝑥𝑖𝑓′ 𝑥𝑖
Aplicaciones de Ingeniería Petrolera
• Equilibrio líquido – vapor.
• Conocer la cantidad de líquido y vapor que existe en unamezcla de hidrocarburos permite calcular algunaspropiedades importantes, tales como: Bo, Bg, Bt, RGA,viscosidad, que serán usados en la evaluación de reservas.
• Algoritmo.9
• La suma de la fase líquida más la fase gas de una mezcladebe ser igual a 1, por lo tanto:
Ing. Juan Carlos Sabido Alcántara
)1...(0.1VL
Aplicaciones de Ingeniería Petrolera
• Equilibrio líquido – vapor.
• Ahora, considerando la fracción molar de cada componente de lamezcla:
• Donde n es el número de componentes de la mezcla, además setiene que:
• La constante de equilibrio 𝑘𝑖 definida como:
•
Ing. Juan Carlos Sabido Alcántara
VyLxz iii )2...(,,2,1 ni
)3...(1
1
1
1 1
vfxy
yx
n
i
ii
n
i
n
i
ii
i
i
ix
yk )4...(,,2,1 ni
Aplicaciones de Ingeniería Petrolera
• Equilibrio líquido – vapor.
• Combinando las anteriores ecuaciones, se tiene:
• Sustituyendo 5 y 6 en 3, se tiene lo siguiente:
Ing. Juan Carlos Sabido Alcántara
...(5) )0.1(0.1 vk
zx
i
ii
)6...(
0.10.1 vk
zky
i
iii
)7...(011
1
1
n
i i
ii
vk
kzvf
Aplicaciones de Ingeniería Petrolera
• Equilibrio líquido – vapor.• Comúnmente para resolver esta ecuación se suponen valores de v,
comenzando en 0.5 que implica 50% de líquido y 50% vapor, y sesustituye en la ecuación 5 para obtener las fracciones de líquido 𝑥𝑖 decada componente, entonces si la suma es igual a 1, se tendrá lafracción de vapor representada por v y la fracción de líquido estará dadapor:
• En caso de que el valor de v supuesto no cumpla con se debesuponer un nuevo valor hasta tener el que más se aproxime, para calcularlas fracciones de vapor 𝑦𝑖 de cada componente se utiliza la ecuación:
Ing. Juan Carlos Sabido Alcántara
n
i
ix1
)8...(1 vl
11
n
i
ix
)9...(* iii xky
Aplicaciones de Ingeniería Petrolera
• Equilibrio líquido – vapor.
• De igual forma que con las fracciones de líquido, la suma esigual a 1.
• La expresión 7 representa una ecuación no lineal y se puederesolver utilizando el método de Newton – Raphson, para esto esnecesario obtener la primera derivada de 7, esta es:
• Entonces el algoritmo de Newton - Raphson para resolver estaecuación queda expresado como:
Ing. Juan Carlos Sabido Alcántara
n
i
iy1
)10...(11
1)('
12
2
n
i i
ii
vk
kzvf
)11...('
1
K
KkK
Vf
VfVV
Aplicaciones de Ingeniería Petrolera
• Equilibrio líquido – vapor.• Si es menor que cierta tolerancia, el algoritmo converge, lo
que nos da el valor de v que permite tener el equilibrio líquido - vaporexacto de la mezcla de hidrocarburos.
Ing. Juan Carlos Sabido Alcántara
)10...(11
1)('
12
2
n
i i
ii
vk
kzvf
)11...('
1
K
KkK
Vf
VfVV
KK VVabs 1
)7...(011
1
1
n
i i
ii
vk
kzvf
Aplicaciones de Ingeniería Petrolera
• Equilibrio líquido – vapor.• Ejercicio: Calcular el equilibrio liquido - vapor de la siguiente mezcla aplicando
el método de Newton – Raphson, Los valores de k se obtienen considerando
Ing. Juan Carlos Sabido Alcántara
)10...(11
1)('
12
2
n
i i
ii
vk
kzvf
)11...('
1
K
KkK
Vf
VfVV
)7...(011
1
1
n
i i
ii
vk
kzvf0.0157 1.18
0.0214 4
0.0037 9.3
0.4921 5.2
0.1038 1.5
0.0594 0.64
0.012 0.32
0.0283 0.25
0.0121 0.132
0.017 0.11
0.0246 0.05
0.2099 0.023
ZiComponente
2H S
2CO
2N
1C
2C
3C
4i C
4n C
5i C
5n C
6C
7C
Suma
iksupuestaV 0.5
iy
ik ix ix ix ix
3000p psi
Aplicaciones de Ingeniería Petrolera
• Equilibrio líquido – vapor.• Ejercicio: Calcular el equilibrio liquido - vapor de la siguiente mezcla aplicando
el método de Newton – Raphson, Los valores de k se obtienen considerando
Ing. Juan Carlos Sabido Alcántara
3000p psi
Vsupuesta
ki f(v) f'(v) Vk+1 xi
H2S 0.01570 1.180 0.0025927 0.0004281 0.591973209 0.01418818
CO2 0.02140 4.000 0.0256800 0.0308160 0.00770916
N2 0.00370 9.300 0.0059631 0.0096104 0.00062570
C1 0.49210 5.200 0.6667161 0.9032928 0.14115302
C2 0.10380 1.500 0.0415200 0.0166080 0.08009342
C3 0.05940 0.640 -0.0260780 0.0114489 0.07548708
i-C4 0.01200 0.320 -0.0123636 0.0127383 0.02008509
n-C4 0.02830 0.250 -0.0339600 0.0407520 0.05089744
i-C5 0.01210 0.132 -0.0185562 0.0284572 0.02488855
n-C5 0.01700 0.110 -0.0272613 0.0437163 0.03592988
C6 0.02460 0.050 -0.0445143 0.0805497 0.05621245
C7+ 0.20990 0.023 -0.4009234 0.7657911 0.49781548
Suma 1.00000 0.17881512 -1.94420876 1.00508544
0.5ZiComponente
Vsupuesta
ki f(v) f'(v) Vk+1 xi
H2S 0.01570 1.180 0.0025539 0.0004154 0.588038761 0.01419726
CO2 0.02140 4.000 0.0231275 0.0249944 0.00774208
N2 0.00370 9.300 0.0051933 0.0072893 0.00062917
C1 0.49210 5.200 0.5928427 0.7142094 0.14182526
C2 0.10380 1.500 0.0400467 0.0154503 0.08021518
C3 0.05940 0.640 -0.0271753 0.0124327 0.07535145
i-C4 0.01200 0.320 -0.0136579 0.0155448 0.01999555
n-C4 0.02830 0.250 -0.0381731 0.0514906 0.05062875
i-C5 0.01210 0.132 -0.0216033 0.0385703 0.02471494
n-C5 0.01700 0.110 -0.0319776 0.0601510 0.03566592
C6 0.02460 0.050 -0.0534018 0.1159250 0.05573641
C7+ 0.20990 0.023 -0.4863657 1.1269729 0.49331808
Suma 1.00000 -0.00859065 -2.18344603 1.00002005
0.591973209Componente Zi
Aplicaciones de Ingeniería Petrolera
• Equilibrio líquido – vapor.• Ejercicio: Calcular el equilibrio liquido - vapor de la siguiente mezcla aplicando
el método de Newton – Raphson, Los valores de k se obtienen considerando
Ing. Juan Carlos Sabido Alcántara
3000p psi
...(5) )0.1(0.1 vk
zx
i
ii
)6...(
0.10.1 vk
zky
i
iii
Vsupuesta
ki f(v) f'(v) Vk+1 xi yi
H2S 0.0157 1.18 0.00255551 0.00041597 0.588023023 0.014197263 0.01675277
CO2 0.0214 4 0.02322663 0.02520917 0.007742077 0.03096831
N2 0.0037 9.3 0.00522226 0.00737082 0.000629174 0.00585132
C1 0.4921 5.2 0.59567743 0.72105588 0.141825257 0.73749134
C2 0.1038 1.5 0.04010783 0.01549748 0.080215182 0.12032277
C3 0.0594 0.64 -0.02712633 0.01238784 0.075351445 0.04822493
i-C4 0.012 0.32 -0.01359673 0.01540592 0.019995546 0.00639857
n-C4 0.0283 0.25 -0.03797076 0.05094624 0.050628751 0.01265719
i-C5 0.0121 0.132 -0.02145197 0.03803199 0.024714943 0.00326237
n-C5 0.017 0.11 -0.03174174 0.05926693 0.035665919 0.00392325
C6 0.0246 0.05 -0.0529478 0.11396217 0.055736412 0.00278682
C7+ 0.2099 0.023 -0.48195435 1.10662216 0.493318078 0.01134632 V 0.5880
Suma -5.299E-10 -2.16617256 1.000020048 0.99998596 L 0.4120
Componente Zi0.588023023
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite.
• El proceso de separación del gas y el aceite es realizado por unaserie de etapas de separación, una etapa de separación se definecomo la condición a la cual el aceite y el gas alcanzan el equilibrio ala presión y temperatura existentes en el separador, en la figura semuestra el esquema de una batería de separación gas – aceite entres etapas.
Ing. Juan Carlos Sabido Alcántara
1ª Etapa 2ª Etapa3ª Etapa
F1 F2,L1 F3,L2 F4,L3
V1,ng1 V3,ng3V2,ng2
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite.
• Para determinar las cantidades de gas y líquidos separados se llevana cabo cálculos de balance de materia, dado que es factible conocerla composición de alimentación al sistema y las condiciones deoperación, es posible determinar la constante de equilibrio físico,esta se determina por medio de correlaciones, gráficas,tabulaciones y ecuaciones de estado, esto lleva al equilibrio líquido– vapor visto anteriormente, obteniendo las fracciones líquido yvapor para cada etapa de separación, de estos balances se puedenobtener factores como: RGA, ºAPI del aceite residual, etc.
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite. Algoritmo.
• Se analiza un sistema en tres etapas como el mostrado en la figura.
• Si son los moles alimentados al sistema, los moles alimentados a la segunda etapa son:
• Y a la tercera
Ing. Juan Carlos Sabido Alcántara
1F
)12...(1112 FlLF
)13...(1121223 FllLlLF
1ª Etapa 2ª Etapa3ª Etapa
F1 F2,L1 F3,L2 F4,L3
V1,ng1 V3,ng3V2,ng2
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite. Algoritmo.
• Los moles retenidos en el tanque de almacenamiento son:
𝐿3 = 𝑙3𝐿2 = 𝑙1𝑙2𝑙3𝐹1…(14)
• En el caso en el que se tenga un número n de etapas se tiene que:
… (15)
Ing. Juan Carlos Sabido Alcántara
1ª Etapa 2ª Etapa3ª Etapa
F1 F2,L1 F3,L2 F4,L3
V1,ng1 V3,ng3V2,ng2
1 2 1 1 1
1
n
m m m i
i
L l l l l F F l
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite. Algoritmo.
• Si los moles en el tanque de almacenamiento por mol dealimentación al sistema son:
…(16)
Ing. Juan Carlos Sabido Alcántara
1ª Etapa 2ª Etapa3ª Etapa
F1 F2,L1 F3,L2 F4,L3
V1,ng1 V3,ng3V2,ng2
1 1.0F
1
n
t i
i
n l
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite. Algoritmo.
• De manera similar el número de moles de gas separado pueden ser calculados para primera, segunda y tercera etapa.
…(17)
…(18)
…(19)
Ing. Juan Carlos Sabido Alcántara
1ª Etapa 2ª Etapa3ª Etapa
F1 F2,L1 F3,L2 F4,L3
V1,ng1 V3,ng3V2,ng2
111 FVV
112222 FlVFVV
1123333 FllVFVV
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite. Algoritmo.
• El volumen de aceite en el tanque de almacenamiento puedeobtenerse a partir de la densidad y el peso molecular del aceite:
…(20)
• El volumen de gas separado por etapa es calculado mediante losmoles de gas liberado y el volumen molar del gas a condicionesestándar es:
…(21)
Ing. Juan Carlos Sabido Alcántara
ot
tt
ot
PMnV
gii VMVVG
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite. Algoritmo.
• La relación gas - aceite definida como volumen de gas a volumen delíquido, a condiciones estándar, puede ser determinada por mediode la expresión:
…(22)
• El peso molecular del aceite en el tanque de almacenamiento escalculado mediante la composición de la fase líquida y el pesomolecular de cada componente:
…(23)Ing. Juan Carlos Sabido Alcántara
tt
otgI
Ot
tt
GI
Ot
I
IPMn
VMVG
PMn
VMVG
V
VGRGA
I
IIt PMXPM
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite. Algoritmo.
• La relación gas - aceite total es la suma de las relaciones gas - aceite por etapa.
…(24)
• La densidad del aceite a condiciones del tanque, se determina apartir de la composición de la fase líquida en el tanque dealmacenamiento y la densidad de cada componente medido acondiciones estándar.
…(25)
Ing. Juan Carlos Sabido Alcántara
M
I
IRGARGAT1
I
II
I
IsiI
otPMX
PMX
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite.
• Tarea 4: A partir de la siguiente información calcular el valor de la para unsistema de separación en tres etapas usando Excel:
Ing. Juan Carlos Sabido Alcántara
totalRGA
1 1
2 2
3 3
654.7 , 140 º
100 , 135 º
14.7 , 104 º
P psia T F
P psia T F
P psia T F
3@ . . 0.8436
203
ot c s
t
gcm
lbPMlb mole
1a etapa 2a etapa 3a etapa
0.0157 1.18 8 32.86
0.0214 4 20 86
0.0037 9.3 138 560
0.4921 5.2 30 180
0.1038 1.5 7 37
0.0594 0.64 2.6 12
0.012 0.32 1.15 4.8
0.0283 0.25 0.9 3.4
0.0121 0.132 0.4 1.45
0.017 0.11 0.33 1.1
0.0246 0.05 0.12 0.37
0.2099 0.023 0.048 0.13
Componente Zi
2H S
2CO
2N
1C
2C
3C
4i C
4n C
5i C
5n C
6C
7C
ik ik ik
Aplicaciones de Ingeniería Petrolera
• Balance de materia en un proceso de separación gas –aceite.
• Tarea 4: El primer paso es calcular el equilibrio líquido – vapor paracada una de las etapas de separación, para esto se utiliza el métodode Newton – Raphson tal y como se explicó anteriormente.
• Los valores de se obtienen considerando . La fracción delíquido obtenida después de la primera etapa de separación es laque se considera como en la segunda etapa.
• Considerar para la ecuación 21 el valor de 𝑉𝑀𝑔 = 379𝑝𝑖𝑒3
𝑙𝑏−𝑚𝑜𝑙que
es el volumen que ocupa 1 𝑙𝑏 − 𝑚𝑜𝑙 de cualquier gas acondiciones estándar. Ing. Juan Carlos Sabido Alcántara
ik 3000p psi
iz
Aplicaciones de Ingeniería Petrolera
• Programa #4: Balance de Materia.• Realizar el programa correspondiente al balance de materia en una batería de
separación de tres etapas.
• El programa deberá de ser capaz de leer los datos de un archivo de texto externopermitiendo al usuario indicar de manera manual la ruta en dónde se ubica dichoarchivo, y presentarlos en pantalla. Deberán enviar el archivo en un bloc de notas conlas características necesarias para que su programa funcione.
• Los resultados deberán de imprimirse en pantalla y en archivo de texto externopermitiendo al usuario indicar de manera manual la ruta en dónde se ubicara dichoarchivo.
• Enviar el código, algoritmo y prueba de escritorio en un archivo de Word al correo de lamateria a más tardar a las 19:00 horas del día anterior a la siguiente clase.
• El nombre del archivo debe seguir la siguiente estructura:
• balance_primerapellido_segundoapellido_nombreIng. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Pruebas de Presión.
• Caracterización estática de un yacimiento: se realiza utilizandoherramientas como registros geofísicos en pozos, datos geofísicos ygeológicos, núcleos, sísmica, entre otros. Todos estos recursosbuscan conocer y definir el tipo de roca que forma el yacimiento, elsello, estructura, es decir el sistema petrolero.
• Caracterización dinámica de un yacimiento: se realiza paraidentificar y evaluar los factores que afectan la explotación de unyacimiento. Esto se consigue mediante el análisis de diferentesvariables que definen el comportamiento del sistema, estas son lapresión, temperatura, tipo de flujo, composición de los fluidos,propiedades petrofísicas, etc. Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Pruebas de Presión.
• Caracterización dinámica de un yacimiento: los elementos queafectan el comportamiento del flujo en los yacimientos son fallasgeológicas, doble porosidad, doble permeabilidad, fracturasinducidas, entre otros.
• Entre las herramientas utilizadas en la caracterización dinámica setienen:– Datos históricos de producción (agua, aceite,gas)
– Pruebas de variación de presión
– Registros de presión de fondo fluyendo y cerrado
– Datos del comportamiento del yacimiento
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Pruebas de Presión.
• Pruebas de Decremento: Se trata de una serie de mediciones depresión en el fondo del pozo durante un periodo de flujo a gastoconstante, generalmente el pozo se encuentra cerrado durante unintervalo de tiempo suficiente para que existan condiciones depresión estables.
• La prueba se realiza colocando un registrador de presión en el fondoen el punto medio del intervalo disparado o frente a la formaciónproductora en el caso de un agujero descubierto.
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Pruebas de Presión.
• Pruebas de Decremento: El objetivo de este tipo de pruebas esobtener parámetros como la permeabilidad (k), el área de drene delpozo, estimar el daño a la formación o de lo contrario laestimulación de la misma, también es posible determinar elvolumen poroso drenado, identificar heterogeneidades y fronterasdel yacimiento.– La principal ventaja técnica de las pruebas de decremento es la posibilidad de
estimar las fronteras físicas del yacimiento, detectar fallas, acuñamientos,heterogeneidades y el volumen poroso drenado.
– La principal desventaja es la dificultad de mantener constante la produccióndurante al prueba.
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Pruebas de Presión.
• Pruebas de Decremento: Aunque las pruebas de decremento noestán limitadas al inicio de la producción de un pozo, este es eltiempo ideal para realizar estas pruebas.
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Pruebas de Presión.
• Pruebas de Incremento: Para llevar a cabo una prueba deincremento de presión es necesario cerrar el pozo productor y queéste haya producido a gasto constante, ya que la mayoría de losmodelos utilizados en las ecuaciones de interpretación estánapoyadas en el principio de superposición y consideran gastoconstante, también se han desarrollado métodos que consideranpresión constante antes del cierre.
• Las pruebas de incremento se diseñan secuencialmente con laspruebas de decremento, con lo que se logran perturbaciones depresión importantes en el medio poroso.
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Pruebas de Presión.
• Pruebas de Incremento: Una prueba de incremento es la medicióncontinua de presión de cierre de un pozo después de un periodo deflujo.
Ing. Juan Carlos Sabido Alcántara
Aplicaciones de Ingeniería Petrolera
• Pruebas de Presión.
• Pruebas de Incremento: Su objetivo es:– Estimar parámetros del yacimiento
– Estimar el factor de daño del pozo
– Determinar la presión media del área de drene
• Las ventajas que se tiene son:– Mediciones suaves de presión
– Gasto constante (q=0)
• Y las desventajas que tienen:– Se tiene que cerrar el pozo (se difiere la producción)
– Dificultad en mantener el gasto constante antes del cierre
Ing. Juan Carlos Sabido Alcántara