Tesis PDF

DISEO E IMPLEMENTACIN DE UN SISTEMA DE
FISCALIZACIN DE GAS PARA LA EMPRESA ATLAS

(ENGINEERING & AUTOMATION)
NESTOR FERNEY SANABRIA VELANDIA
JUAN DAVID RODRGUEZ FONSECA
Monografa para optar el ttulo de Ingeniero en Control
Tutor:
EDUARDO ALBERTO DELGADILLO GOMEZ
INGENIERO EN ELCTRNICA
Universidad Distrital Francisco Jos de Caldas
Facultad Tecnolgica
Tecnologa en Electrnica
Bogot, 2013
ndice general
Agradecimientos 9
Introduccin 11
Resumen 15
1. Captulo 1 17
1.1. Planteamiento del problema . . . . . . . . . . . . . . . . . . . . . . . 17
1.1.1. Descripcin . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.1.2. Identicacin . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.1.3. Formulacin . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.2. Objetivos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.2.1. General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.2.2. Especcos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.3. Justicacin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2. Captulo 2 21
2.1. Antecedentes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.1.1. Actualizacin Del Mdulo leo Del Sistema Centinela Apli-
cando La Norma Api Para El Clculo De Cantidades De Pe-
trleo, En Los Tanques De Almacenamiento Presentes En Los
Patios De Tanques Y Terminales De Embarque Del Oriente
De Venezuela . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.1.2. Disminucin Del ndice De Prdidas En Diferentes reas De
La Gerencia Complejo Barrancabermeja (GCB) De ECOPE-
TROL A Partir Del Mejoramiento De La Medicin Esttica
Y Dinmica. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.1.3. Diseo De Un Sistema De Descarga De Crudo Pesado Desde
Carrotanques Para Campo Rubiales, Meta En La Empresa
Ingeniera STRYCO, LTDA . . . . . . . . . . . . . . . . . . . 22
2.1.4. Implementacin De Un Sistema De Medicin Y Fiscalizacin
De Custodia De Crudo En La Estacin De Transferencia Moriche 23
2.1.5. Diseo De Un Sistema De Medicin Y Control De Uso De
Combustibles Para Las Centrales De Generacin Termoelc-
trica Del S. N. I . . . . . . . . . . . . . . . . . . . . . . . . . . 24
i
ndice general ndice general
2.1.6. Diseo De Un Sistema De Medicin Automtico En Lnea De
Produccin En Las Principales Estaciones De Flujo Del La
Unidad De Explotacin Tierra Oeste (PDVSA) . . . . . . . . 24
2.2. Marco terico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.2.1. HMI . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.2.2. Medidor de Coriolis . . . . . . . . . . . . . . . . . . . . . . . . 26
2.2.3. PLC (Programing Logic Control) . . . . . . . . . . . . . . . . 27
2.2.4. Sistema de Descarga de Hidrocarburos. . . . . . . . . . . . . . 28
2.2.5. GLP Gas Licuado de Petrleo . . . . . . . . . . . . . . . . . . 29
2.2.6. Elementos Principales Del Proceso . . . . . . . . . . . . . . . 30
2.2.7. NEMA 4X . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.2.8. RS-485 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
3. Captulo 3 39
3.1. Metodologa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
3.1.1. Seal de pulsos . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3.1.2. Seal 4-20mA . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3.1.3. PLC KOYO DL06 . . . . . . . . . . . . . . . . . . . . . . . . 41
3.1.4. HMI . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3.1.5. Impresora EPSON TM-T90 . . . . . . . . . . . . . . . . . . . 44
3.2. Desarrollo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.2.1. PRUEBAS INICIALES . . . . . . . . . . . . . . . . . . . . . . 44
3.2.2. TOTALIZADOR . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.2.3. FUNCIONES DEL PROGRAMA DEL HMI . . . . . . . . . . 50
3.2.4. MONTAJE DE COMPONENTES EN UN TABLERO ELC-
TRICO NEMA 4X . . . . . . . . . . . . . . . . . . . . . . . . 65
3.2.5. INSTALACIN EN CAMPO . . . . . . . . . . . . . . . . . . 71
4. Captulo 4 75
4.1. CORRECCIN DE VOLUMEN PARA GLP . . . . . . . . . . . . . 75
4.1.1. API M.P.M.S 11.2.4 . . . . . . . . . . . . . . . . . . . . . . . 75
4.1.2. Eleccin de la Tabla apropiada . . . . . . . . . . . . . . . . . 76
4.1.3. TABLA 54E . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
4.1.4. Ejemplo de Aplicacin . . . . . . . . . . . . . . . . . . . . . . 85
4.1.5. Criterios para una implementacin del algoritmo en PLC . . . 85
4.2. AGA 8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
4.2.1. Factor de compresibilidad Z . . . . . . . . . . . . . . . . . . . 87
4.2.2. MTODO DE LA CARACTERIZACIN DEL DETALLE . . 87
4.2.3. ECUACIONES GENERALES . . . . . . . . . . . . . . . . . . 88
4.2.4. ECUACIONES PARA EL FACTOR DE COMPRESIBILIDAD 89
4.2.5. ALGORITMO GENERAL AGA REPORTE 8 EN MATLAB 92
ii
ndice general
5. Captulo 5 95
5.1. Resultados . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
5.1.1. Comparacin entre LectroCount LCR y el sistema Totalizador
desarrollado. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
5.1.2. Instalacin en Campo . . . . . . . . . . . . . . . . . . . . . . 97
5.1.3. Estudios AGA reporte 8 . . . . . . . . . . . . . . . . . . . . . 102
6. CONCLUSIONES 117
Bibliografa 119
Nomenclatura 121
A. PLANO P&ID 129
B. PLANO UNIFILAR 131
C. PLANO ELCTRICO 133
D. CONEXIONES DE TRANSMISION DE DATOS 135
E. COMPARACIN DE FUNCIONES ENTRE TOTALIZADORES 137
F. TABLAS Y ALGORITMO DEL AGA REPORTE 8 139
G. TABLA 54E API MPMS 11.2.4 155
iii
ndice de guras
1.1. Sistema a cambiar de la empresa PEGASUS BLENDING . . . . . . . 17
1.2. Sensor coriolis RHM 20 instalado en PEGASUS BLENDING . . . . . 18
2.1. Aplicacin HMI. [1] . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.2. Vibracin del tubo del sensor de Coriolis. [2] . . . . . . . . . . . . . . 27
2.3. Torsin de los tubos del sensor de coriolis. [2] . . . . . . . . . . . . . 27
2.4. PLC. [3] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.5. PLC KOYO DL06. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.6. MAPLE HMI5043N. . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.7. Impresora EPSON TM-T90. . . . . . . . . . . . . . . . . . . . . . . . 32
3.1. Diagrama general para la solucin del sistema . . . . . . . . . . . . . 39
3.2. Conguracin de salida de pulsos del transmisor[4]. . . . . . . . . . . 41
3.3. PLC Koyo DL06[5]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
3.4. Conguracin de entradas y salidas del PLC[5]. . . . . . . . . . . . . 42
3.5. Conguracin de entradas y salidas drenadoras/surtidoras[5] . . . . . 43
3.6. Mdulo analgico F0-04AD-1[6]. . . . . . . . . . . . . . . . . . . . . . 44
3.7. Montaje del sistema en el laboratorio ATLAS . . . . . . . . . . . . . 45
3.8. Plano P&ID Laboratorio ATLAS . . . . . . . . . . . . . . . . . . . . 45
3.9. Gabinete del totalizador inicial . . . . . . . . . . . . . . . . . . . . . 46
3.10. Vista interna del gabinete inicial . . . . . . . . . . . . . . . . . . . . . 46
3.11. Seleccin de Pulsos/Kg, mediante tabla . . . . . . . . . . . . . . . . . 48
3.12. Ajuste del K-factor . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.13. Eleccin de totalizar masa/volumen . . . . . . . . . . . . . . . . . . . 49
3.14. Unidades de medida para el totalizador. . . . . . . . . . . . . . . . . 50
3.15. Pantalla Principal del programa . . . . . . . . . . . . . . . . . . . . . 51
3.16. pantalla PRODUCTO . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.17. Datos de medida . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.18. Formato de Ticket numero 1 . . . . . . . . . . . . . . . . . . . . . . 53
3.19. Ticket tipo 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.20. Pantalla del submenu REGISTRO . . . . . . . . . . . . . . . . . . . . 54
3.21. Opcin IMPRIMIR PANTALLA . . . . . . . . . . . . . . . . . . . . 55
3.22. Impresin del total del da . . . . . . . . . . . . . . . . . . . . . . . . 55
3.23. Mensaje: Elemento protegido por contrasea . . . . . . . . . . . . . . 56
3.24. submenu CONFIG . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.25. Password . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
1
ndice de guras ndice de guras
3.26. Unidad de medida . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.27. Canales de 4-20mA . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.28. Unidades para ujo msico . . . . . . . . . . . . . . . . . . . . . . . . 59
3.29. Error en el Canal de 4-20mA . . . . . . . . . . . . . . . . . . . . . . . 60
3.30. Conguracin de fecha y hora . . . . . . . . . . . . . . . . . . . . . . 60
3.31. Conguracin de datos del recibo . . . . . . . . . . . . . . . . . . . . 61
3.32. submenu CALIBRAR . . . . . . . . . . . . . . . . . . . . . . . . . 62
3.33. Ventana de Calibracin-Delta . . . . . . . . . . . . . . . . . . . . . . 63
3.34. Calibracin-numero de bombas . . . . . . . . . . . . . . . . . . . . . 64
3.35. ABOUT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
3.36. Placa caracterstica del gabinete . . . . . . . . . . . . . . . . . . . . . 65
3.37. Conexin entre la fuente de 24VDC y el HMI . . . . . . . . . . . . . 66
3.38. Distribucin del tablero . . . . . . . . . . . . . . . . . . . . . . . . . . 66
3.39. Elementos de alimentacin y proteccin . . . . . . . . . . . . . . . . . 67
3.40. Entadas y salidas del sistema . . . . . . . . . . . . . . . . . . . . . . 68
3.41. SSRs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
3.42. Cableado interno y montaje de componentes . . . . . . . . . . . . . . 69
3.43. Vista frontal del gabinete . . . . . . . . . . . . . . . . . . . . . . . . . 69
3.44. Vista interior de la puerta del gabinete . . . . . . . . . . . . . . . . . 70
3.45. Entradas y salidas del gabinete . . . . . . . . . . . . . . . . . . . . . 71
3.46. Instalacin de los elementos . . . . . . . . . . . . . . . . . . . . . . . 72
3.47. Alimentacin disponible. . . . . . . . . . . . . . . . . . . . . . . . . . 72
3.48. Punto donde se desmonto el sensor . . . . . . . . . . . . . . . . . . . 73
4.1. Ejemplo de la tabla 54E. [7] . . . . . . . . . . . . . . . . . . . . . . . 86
5.1. LectroCount LCR 600. [8] . . . . . . . . . . . . . . . . . . . . . . . . 95
5.2. Totalizador ATLAS E&A . . . . . . . . . . . . . . . . . . . . . . . . . 96
5.3. Zona de instalacin . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
5.4. Instalacin Coriolis . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
5.5. Instalacin totalizador y transmisor RHE08 . . . . . . . . . . . . . . 98
5.6. Soporte para el totalizador y transmisor . . . . . . . . . . . . . . . . 99
5.7. Bombas dentro del proceso . . . . . . . . . . . . . . . . . . . . . . . . 100
5.8. Totalizador y transmisor RHE12 . . . . . . . . . . . . . . . . . . . . . 100
5.9. Sistema en funcionamiento . . . . . . . . . . . . . . . . . . . . . . . . 101
5.10. Visualizacin del sistema en funcionamiento . . . . . . . . . . . . . . 101
5.11. Funcionamiento totalizador y RHE08 . . . . . . . . . . . . . . . . . . 102
5.12. Funcionamiento totalizador y RHE12 . . . . . . . . . . . . . . . . . . 102
5.13. Metano . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
5.14. Nitrgeno . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
5.15. Dixido de carbn . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
5.16. Etano . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
5.17. Propano . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
5.18. 6. Agua (Vapor) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
2
ndice de guras
5.19. 7. Sulfato de hidrgeno . . . . . . . . . . . . . . . . . . . . . . . . . 113
5.20. 8. Hidrgeno . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
5.21. 10 Oxgeno . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
5.22. 20. Helio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
5.23. 21. Argn . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
5.24. 11. i- Butano . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
5.25. 13. i-pentano . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
F.1. Diagrama de ujo para el algortmo del clculo del factor de compre-
sibilidad Z . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
G.1. Prodecimiendo de implementacin para la Tabla 54E[7] . . . . . . . . 156
3
ndice de cuadros
2.1. Clasicacin reas restringidas . . . . . . . . . . . . . . . . . . . . . . 36
3.1. Eleccin de cantidad de pulsos sobre unidad[4]. . . . . . . . . . . . . 40
4.1. Parmetros de los uidos de referencia[7] . . . . . . . . . . . . . . . . 77
5.1. GULF COAST . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
5.2. AMARILLO . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
5.3. EKOFISK . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
5.4. HIGH N2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
5.5. HIGH CO2-N2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
F.1. Tabla 4 del AGA reporte 8 . . . . . . . . . . . . . . . . . . . . . . . . 140
F.2. Continuacin Tabla 4 del AGA reporte 8 . . . . . . . . . . . . . . . . 141
F.3. Tabla 5 AGA reporte 8 . . . . . . . . . . . . . . . . . . . . . . . . . . 142
F.4. Tabla 6 AGA reporte 8 . . . . . . . . . . . . . . . . . . . . . . . . . . 143
F.5. Continuacin Tabla 6 AGA reporte 8 . . . . . . . . . . . . . . . . . . 144
F.6. Tabla A.5-1 AGA reporte 8 . . . . . . . . . . . . . . . . . . . . . . . 145
F.7. Tabla A. 5-2, AGA reporte 8 . . . . . . . . . . . . . . . . . . . . . . . 146
5
ndice de algoritmos
F.1. AGA reporte 8 en MATLAB . . . . . . . . . . . . . . . . . . . . . . 148
G.1. Implementacin de la tabla 54E de API MPMS 11.2.4 . . . . . . . . . 157
7
Agradecimientos
Primero y antes que nada, dar gracias a Dios, por estar con nosotros en cada paso,
por darnos fortaleza e iluminar nuestras mentes y por haber puesto en el camino
a aquellas personas que han sido soporte y compaa durante todo el periodo de
estudio.
Agradecer hoy y siempre a nuestras familias por el esfuerzo realizado por ellos. El
apoyo en los estudios, de ser as no hubiese sido posible. A nuestros padres y dems
familiares ya que siempre nos brindan el apoyo, la alegra y nos dan la fortaleza
necesaria para seguir adelante.
Del mismo modo agradecemos a la Universidad Distrital Francisco Jos de Caldas
y al estado de Colombia por habernos dado la oportunidad de tener un estudio y
una formacin que nos permita salir adelante y poder competir por los altos cargos
en nuestra rea de estudio.
Un agradecimiento especial al ingeniero Alberto Delgadillo y a la empresa ATLAS
E&A, por la enseanza, la colaboracin, el apoyo en cada fase del proceso de elabo-
racin de este proyecto y sobre todo por la oportunidad brindada.
9
Introduccin
El presente trabajo de grado hace referencia al desarrollo e implementacin de un
sistema de scalizacin para la empresa ATLAS E&A para un proceso de medicin
en una transferencia de custodia de Gas licuado de petrleo (GLP) , hidrocarburo
derivado del petrleo, con el n de obtener mayor control y facilitar la labor de
operacin sobre este proceso.
El sistema se compone de varios equipos entre los que se encuentran el sensor y
transmisor de ujo, un totalizador y una impresora. Los dos primeros (sensor y
transmisor) se encargan de la medida directa del proceso y transmisin de datos al
totalizador, el totalizador est compuesto por un PLC y pantalla HMI y conjun-
tamente son los encargados de realizar todas la operaciones para la supervisin y
control del proceso, y por ltimo la impresora es la encargada de entregar un tiquete
como constancia de la descarga realizada.
El sistema cuenta con las funciones bsicas y avanzadas que se listan a continuacin:
FUNCIONES BSICAS
- Totalizar por masa.
- Seleccin PRESET mltiple.
- Control vlvula.
- Control bombas.
- Registro de transferencias.
- Impresin de tiquete.
- Visualizacin de mltiples datos.
- Parada de emergencia.
- Opciones de rearme.
FUNCIONES AVANZADAS
- Conguracin Fecha/Hora.
- Seleccin del tipo de tiquete.
11
Introduccin
- Seleccin unidades de medida.
- Seleccin modo manual.
- Proteccin por contrasea.
- Lectura de 2 canales 4-20mA
- Calibracin K-Factor
- Calibracin Delta para cierre de vlvula y bombas.
Adicional a esto el sistema fue elaborado de tal manera que cumpla con las normas
mnimas para operacin en campo y por tanto el sistema es de tipo NEMA 4X
INTEMPERIE.
Dentro del desarrollo fuerte de este trabajo se encuentra el desarrollo de un algoritmo
para correccin de la medida por medio del factor de compresibilidad y su clculo
basado en la norma AGA reporte 8, esto con el n de generar un estudio para la
correccin de la medida de volumen de distintos tipos de mezclas de gases.
Adicionalmente se incluye el clculo de compensacin por temperatura para el GLP
desarrollado en base a las normas API, esto como un estudio y un fortalecimiento
del conocimiento para la elaboracin de sistemas de totalizacin, igualmente visto
para futuros trabajos o mejoras para el sistema.
El sistema totalizador implementado y entregado para la empresa ATLAS E&A no
cuenta con la implementacin de estos algoritmos de correccin debido a que no
eran necesarios puesto que la totalizacin del GLP se hace por masa principalmente
y no requiere de correcciones por la gran exactitud y precisin del instrumento de
medida, por lo tanto las normas AGA 8 y API fueron desarrolladas y trabajadas en
MATLAB ms como un objeto de estudio fuerte para la parte de ingeniera y como
cumplimiento de uno de los objetivos de este proyecto mas no para la implementacin
como tal del sistema.
El propsito de desarrollar este trabajo es el dar solucin a una necesidad para un
proceso en especco pero que al mismo tiempo pueda ser fcilmente expandible
y aplicable a otros procesos del mismo tipo, usando tecnologa de calidad pero de
menor costo comparado con instrumentos que pudieran realizar la misma funcin y
con esto no solo se ahorran costos sino que se obtiene un sistema al gusto del usuario
y se contribuye al desarrollo de la ingeniera colombiana.
En el marco metodolgico de este proyecto se realiz un previo estudio y recopilacin
de informacin referente a la necesidad que se tena, con esto se conform un estado
del arte y se hizo un anlisis de los requerimientos dados por la empresa para el
desarrollo del sistema. Posteriormente se compilan las posibles soluciones concluidas,
se desarrollan los modelos del sistema y conforme se elabora se va vericando el
funcionamiento para inmediatamente realizar los ajustes que se requieran.
Puesto que el proceso es delicado y no se puede detener para probar el sistema en
campo, es necesario probar el sistema antes en un montaje que simula el proceso y
que permita que en campo solo se realice la instalacin y puesta en funcionamiento.
12
Introduccin
Finalmente en este trabajo se presenta un marco terico con conceptos bsicos para
el desarrollo y entendimiento del sistema implementado, los antecedentes compuestos
por trabajos y proyectos que permitan la comparacin de similitudes y diferencias
con nuestro proyecto, y as mismo de un punto de partida para el desarrollo de este
proyecto, y nalmente se elabora todo el desarrollo de este proyecto paso a paso,
adjuntando material de apoyo como planos y tablas.
13
Resumen
Desde el descubrimiento del petrleo y sus mltiples usos, el mundo se ha movido
a costa de este recurso el cual no es renovable, esto a tal punto que nos hemos
hecho bastante dependientes del mismo y de agotarse no hay una fuente alternativa
para suplir el recurso. Mientras se encuentran recursos alternativos es necesario que
la medicin de este recurso sea correcta durante su comercializacin para evitar
desperdicios y prdidas.
La medicin correcta de hidrocarburos tiene una suprema importancia en la indus-
tria de petrleo, puesto que se relaciona de manera directa con las negociaciones que
se realizan con los clientes. La empresa ATLAS deseaba implementar un sistema de
medicin automtico para GAS tipo GLP (Gas Licuado de Petrleo), un hidrocar-
buro derivado del petrleo, de tal manera que les permita a ellos tener la lectura de
la masa de este tipo de hidrocarburo y de igual forma poder realizar una supervisin
al proceso de descarga. Por esto se dise e implement un sistema de scalizacin,
el cual se compone de un medidor de ujo de efecto coriolis, transmisor de ujo, un
sistema totalizador de masa y una impresora de tiquetes. Con el sistema de sca-
lizacin se desea que los operarios tengan lectura en todo momento de la medicin
del GLP durante el proceso de descarga mediante el uso de distintas funciones del
totalizador, y al nal una impresin de un tiquete de dicha descarga.
Puesto que inicialmente se tena la medicin para cualquier tipo de GAS se plante
un objetivo adicional para corregir posibles errores en la medida que el instrumento
pueda tener y se desarrolla un algoritmo de correccin del volumen por medio del
factor de compresibilidad por presin y temperatura basados en el reporte AGA 8.
Puesto que nalmente el GAS medido es de tipo GLP se adiciona un estudio basado
en la norma API para la correcin de medida de volumen para este tipo de GAS,
sin embargo en el sistema nal implementado no fueron necesarias estas normas de
correccin y estas quedan solo como un estudio y aporte para mejoras del sistema y
futuras implementaciones que si requieran este tipo de correcciones.
15
1. Captulo 1
1.1. Planteamiento del problema
Actualmente el sector de la industria de hidrocarburos tiene la necesidad de im-
plementar sistemas de medida que le permitan obtener valores exactos para tener
un mayor control y evitar prdidas por malgaste o robo de hidrocarburos y de esta
manera tener ciertas ganancias econmicas para la empresa, por tal motivo se hacen
tan importantes los sistemas de scalizacin los cuales pretenden dar solucin a es-
te problema de medicin, combinando equipos y herramientas automticas para la
supervisin de los hidrocarburos y haciendo uso de medidores de ujo incorporados
al sistema.
La empresa ATLAS E&A se dispuso a disear implementar un sistema de totali-
zacin y registro de masa para la empresa PEGASUS BLENDING, en la cual se
encuentraba instalado un sistema totalizador de esta variable ver Figura 1.1 para el
control del proceso scal de GLP (Gas Licuado de Petrleo).
Figura 1.1.: Sistema a cambiar de la empresa PEGASUS BLENDING
El sistema mostrado en Figura 1.1 cumple la funcin de realizar un proceso de
17
Captulo 1 Captulo 1
scalizacin en descarga del producto en carro-tanques y est conformado por un
sensor de ujo tipo coriolis de referencia RHM 20 ver Figura 1.2 el cual es el sensor
primario y mide directamente el ujo msico de la sustancia que pasa atra vez de
el. Transmisor de ujo RHE 08 el cual totaliza la masa que pasa a traves del sensor
y entrega este valor por medio de una seal de pulsos. Contador Lectro Count LCR
II el cual muestra en pantalla la masa totalizada y envia una orden de impresin a
la impresora EPSON TMU295 la cual imprime un recibo de la descarga realizada.
Figura 1.2.: Sensor coriolis RHM 20 instalado en PEGASUS BLENDING
Este sistema no cuenta con una base de datos para llevar el control de las operaciones
realizadas por la empresa. El proceso de apertura, cierre y circulacin del uido se
realizaba por medio de una valvula manual, y el encedido de las bombas las cuales
hacian circular el uido eran activadas de forma manual.
1.1.1. Descripcin
En el sector petrolfero se cuenta con gran cantidad de equipos para desarrollar un
sistema de scalizacin, como computadores de ujo, transmisores y totalizadores,
pero estos equipos son de alto costo y para necesidades generales. Por esta razn se
requiere de un sistema de menor costo y que se adapte a necesidades ms especcas
de la empresa.
1.1.2. Identicacin
En la empresa ATLAS se implement un sistema de scalizacin para la medida del
ujo msico de GLP, el cual cuenta con una lectura de la medida con y haciendo
uso de equipos de bajo costo comparados a los equipos del mercado.
18
1.2 Objetivos
Un sistema de scalizacin hace uso de elementos de ltima tecnologa tales co-
mo medidores de ujo, transmisores y herramientas de comunicacin que se pueden
incorporar a un sistema HMI SCADA y de esta manera tener una completa su-
pervisin del sistema en el rea de medicin de GAS, dndole mayor control a los
trabajadores involucrados, facilitndoles sus tareas y mejorando la conabilidad en
el sistema.
1.1.3. Formulacin
El propsito del proyecto fue desarrollar un sistema de scalizacin que haga uso de
un medidor digital de ujo de gran exactitud y un transmisor de ujo incorporados
a un sistema totalizador para conformar un sistema de scalizacin y darle mayor
control a la empresa sobre su proceso de medicin realizando una base de datos de
las descargas que se lleven acabo y manejando de forma autompatica los actuadores
referentes al proceso.
Adicionalmente generar un estudio de la norma AGA por medio del desarrollo de
un algortimo de correcin por el factor de compresibilidad tratado en el reporte 8
se hara un acercamiento de el proceso de transferencia de custodia en gases.
1.2. Objetivos
De acuerdo a el problema planteado se plantearon los siguientes objetivos para dar
solucin a la propuesta descrita.
1.2.1. General
Disear e implementar un sistema de scalizacin para la empresa ATLAS para la
medicin de ujo msico de GAS en la zona de descarga mediante el uso de un
medidor de coriolis y un totalizador por medio del uso de un PLC, pantalla HMI e
impresora de tiquetes.
1.2.2. Especcos
Realizar un totalizador utilizando un transmisor de ujo, PLC y una pantalla
HMI en donde se muestre la masa totalizada y se guarden los datos registrados.
Disear e implementar un men para el totalizador dentro del HMI el cual ten-
dr como principales funciones: impresin de recibo de la descarga realizada,
visualizacin del estado de la descarga del GAS y de los datos guardados.
19
Captulo 1 Captulo 1
Desarrollar un algoritmo el cual simule la compensacin de gas descrita en
AGA 8 (American Gas Association report 8).
Integrar el PLC y HMI en un tablero elctrico para la instalacin y operacin
en campo.
1.3. Justicacin
Actualmente el sector de la industria de hidrocarburos tiene la necesidad de imple-
mentar sistemas de medida que le permitan obtener valores exactos para tener un
mayor control y evitar prdidas por malgaste o robo de hidrocarburos en el proceso
de transferencia de custodia y de esta manera obtener ciertas ganancias econmicas
para la empresa, por tal motivo se hacen tan importantes los sistemas de scalizacin
los cuales pretenden dar solucin a este problema de medicin, combinando equipos
y herramientas automticas para la supervisin de los hidrocarburos y haciendo uso
de medidores de ujo incorporados al sistema.
El propsito de este proyecto es desarrollar un sistema de scalizacin que haga uso
de un medidor digital de ujo de gran exactitud y un transmisor de ujo incorporados
a un sistema totalizador para conformar un sistema de scalizacin y darle mayor
control a la empresa sobre su proceso de medicin al momento de la descarga o
transferencia de custodia del producto.
Un aspecto muy importante en cuanto a estos sistemas, es que hay gran variedad
de instrumentos pero la gran mayora estn diseados para aplicaciones generales y
son de un muy alto costo. Por este motivo se desea elaborar un totalizador para el
sistema de scalizacin, de tal manera que sea bastante parecido a los del mercado
pero que se adapte en una forma ms especca a nuestra aplicacin, pero a un
menor costo, todo esto genera que se tenga desempeo y funcionamiento al gusto
del cliente y al tiempo se est contribuyendo a desarrollar ingeniera colombiana con
este tipo de mecanismos.
Finalmente estos sistemas deben estar soportados sobre normas para que tengan
valides, por lo cual se asegura que los equipos cumplan las normas mnimas para
operar en campo bajo la norma NEMA 4X
Por otra parte se desea desarrollar un algoritmo de correccin de compresibilidad
por temperatura y presin basados en el reporte AGA 8 el cual hace referencia a
manipulacin de GAS, esto como solamente como objeto de estudio para el fortale-
cimiento del conocimiento de este tipo de correciones y mejoras para otros sistemas
totalizadores o iimplementaciones futuras que puedan requerir de esta norma.
20
2. Captulo 2
2.1. Antecedentes
A continuacin se realiza un compilado de trabajos realizados en el sector de los hi-
drocarburos, ms especcamente en cuanto a su medicin, supervisin y transporte,
con el n de sentar las bases de comparacin con el proyecto de grado a desarrollado
Diseo E Implementacin De Un Sistema De Fiscalizacin De Hidrocarburos Para
La Empresa ATLAS (ENGINEERING & ATUMATION) y de esta manera apo-
yarse en estos trabajos para tener un mejor punto de partida y conocimiento para
el buen desarrollo del proyecto a implementar.
2.1.1. Actualizacin Del Mdulo leo Del Sistema Centinela
Aplicando La Norma Api Para El Clculo De Cantidades
De Petrleo, En Los Tanques De Almacenamiento
Presentes En Los Patios De Tanques Y Terminales De
Embarque Del Oriente De Venezuela
Esta es una tesis desarrollada en la Universidad de Oriente de Venezuela, cuyo n
es controlar las actividades petroleras en Venezuela, pertinentes a la produccin,
manejo y almacenamiento de crudo en el pas, se utiliza la herramienta compu-
tacional Sistema CENTINELA, que permite a los ingenieros y trabajadores de la
Industria Petrolera: consultar, actualizar e incluir los datos, referidos al negocio de
hidrocarburos.
Entre el Sistema CENTINELA, se encuentra el Mdulo leo, que se encarga del re-
gistro y procesamiento de las cantidades de crudo en los tanques de almacenamiento,
para su utilizacin en todas las locaciones donde sea necesaria esa informacin. Esta
unidad presentaba un desfase en cuanto a la medicin de crudos en tanques, con
respecto las normas publicadas por el Instituto Americano de Petrleo (API) en su
reglamento Manual de Estndares de Medicin de Hidrocarburos, especcamente
en la seccin de clculo de cantidades de petrleo. Las diferencias que se presenta-
ban causaban varios problemas en la exportacin y comercio del crudo por lo cual
el objetivo era ajustar el sistema CENTINELA para que se ajustara a las normas
API correspondientes.
21
Captulo 2 Captulo 2
Para poder realizar esto se hizo un estudio previo de todas las locaciones de medicin
del combustible y una vez determinadas las variables del proceso y el algoritmo
implementado a la medicin se cambia dicho algoritmo para ajustarlo a las normas
API basado en los estudios previamente realizados.
La principal diferencia entre este trabajo realizado y el que se quiere realizar, se basa
en que en esta tesis ya se contaba con un sistema de medicin especco que requeria
ajustes para que cumplieran con las normas API, mientras que en el proyecto a
realizar no se cuenta con ningn sistema de medicin, por lo cual se implementa
uno completamente diferente pero de igual manera basados en la norma API entre
otros. [9]
2.1.2. Disminucin Del ndice De Prdidas En Diferentes reas
De La Gerencia Complejo Barrancabermeja (GCB) De
ECOPETROL A Partir Del Mejoramiento De La
Medicin Esttica Y Dinmica.
Esta es una tesis desarrollada en la Universidad Industrial de Santander con el
n de realizar actividades que ayudarn a mejorar las prcticas de medicin a los
trabajadores y la vericacin y control de las mediciones manuales y automticas.
El objetivo principal de este trabajo fue el de mejorar el porcentaje de prdidas
no identicadas para mejorar el balance msico de la renera. En este se realiza el
estudio de las herramientas y tecnologas existentes, siguiendo con la actualizacin y
validacin de las tablas en los sistemas de informacin, el estudio de la medicin del
combustleo, la realizacin de procedimientos y recomendaciones que mejoran los
procesos en las plantas. Para esto se hace uso de la telemetra o medicin automtica
de nivel y temperatura.
En cuanto a nuestro proyecto, este es diferente en cuanto a que ya se tiene un sistema
de medicin que aplican varias normas dentro de una empresa de gran prestigio como
ECOPETROL y lo que se quiere es realizar un estudio para mejorar este sistema,
mientras nuestro proyecto se desea implementar un sistema de medicin partiendo
de 0, es decir, uno completamente nuevo, pero en general los estudios y modo de
implementacin es similar puesto que se trabajan con combustibles. [10]
2.1.3. Diseo De Un Sistema De Descarga De Crudo Pesado
Desde Carrotanques Para Campo Rubiales, Meta En La
Empresa Ingeniera STRYCO, LTDA
Tesis elaborada en la Universidad Industrial de Santander, el cual se desarroll
para mejorar las condiciones de descarga de crudo desde Carrotanques en el campo
petrolero Campo Rubiales, Meta, con el n de hacer esta operacin ms rpida
22
2.1 Antecedentes
y eciente, para aumentar el nmero de carrotanques descargados y aumentar el
volumen de crudo pesado enviado diariamente al rea de procesos.
Para la realizacin de este proyecto se realiz un estudio previo de las tuberas y
bombas del sistema que cumplieran con los requisitos para disminuir el tiempo de
descarga en los carrotanques y dar mayor exibilidad de operacin para el aprove-
chamiento de las facilidades. Una vez realizado el estudio y obtenidos la tubera,
bombas y dems elementos necesarios se realiz la implementacin y se obtuvo un
tiempo de descarga de 40 minutos para cada carrotanque.
La diferencia en cuanto al proyecto que deseamos implementar es el uido que se
maneja, puesto que en este proyecto se trabaja con crudo pesado mientras en nuestra
implementacin se trabaja con combustibles ligeros tales como gas y diesel, adems
en nuestro caso la tubera ya est implementada y nos concentramos en la medicin
y el trabajo en un sistema de HMI y SCADA. [11]
2.1.4. Implementacin De Un Sistema De Medicin Y
Fiscalizacin De Custodia De Crudo En La Estacin De
Transferencia Moriche
Proyecto desarrollado e implementado en la Universidad Industrial de Santander,
en el cual se hizo la implementacin de la unidad de medicin dinmica de crudo en
el Campo Moriche, Puerto Boyac.
Para el proyecto se realiz un estudio para hacer la diferenciacin entre medicin
esttica y dinmica de hidrocarburos lquidos, de igual manera la medicin y ca-
libracin de instrumentos, teniendo en cuenta normas y procedimientos de funcio-
namiento as como el estudio econmico para la implementacin de la unidad de
medicin.
Una vez obtenidos los datos de los estudios se hace la eleccin de los mejores equipos
y herramientas para el sistema de medicin entre estos el ms importante el medidor
de desplazamiento o de turbina dependiendo el tipo de hidrocarburo lquido que se
maneje. De igual manera se tienen en cuenta las normas API para cada momento
de la implementacin del sistema de medicin.
Este proyecto es uno de los ms similares al que se est desarrollando y sobre todo
los estudios realizados en este sirven de base para poder realizar la seleccin de los
mejores instrumentos en nuestra implementacin para el sistema de medicin. [12]
23
Captulo 2 Captulo 2
2.1.5. Diseo De Un Sistema De Medicin Y Control De Uso
De Combustibles Para Las Centrales De Generacin
Termoelctrica Del S. N. I
Este es un proyecto realizado en la Escuela Politcnica Nacional en Ecuador, con
el objeto de disear un sistema de medicin y control de uso de los combustibles
aplicable a todas las centrales de generacin termoelctrica que actualmente operan
en el pas y forman parte del sistema nacional interconectado.
En este proyecto se comienza realizando un estudio de la generacin termoelctrica
en el pas y de acuerdo a esto se determinan los tipos de medidores de ujos nece-
sarios para la implementacin. Posterior a esto se hace la descripcin del sistema
SCADA local y de las interfaces hombre-mquina. Ya obtenidos todos los elementos
necesarios y diseos terminados se presenta la normativa necesaria para los proce-
dimientos tcnicos elctricos y mecnicos necesarios para la implementacin, junto
con una gua para la determinacin de las reas clasicadas dentro de cada planta.
Todo esto es con el n de obtener un mayor control sobre el uso de combustibles en
el sector elctrico, evitando as el robo y malgaste de los recursos hidrocarburferos
y as alcanzando un ahorro econmico notable.
Este proyecto es similar puesto que se encarga de realizar un sistema de supervisin
y medicin de los hidrocarburos de tal manera que se tenga un mayor control y se
logre un ahorro evitando prdidas dentro del sistema, sin embargo, los sistemas a
implementar son distintos puesto que los requerimientos no son iguales debido al
uso nal de los combustibles y de igual manera el sistema SCADA tambin tendr
otras caractersticas diferentes de acuerdo a lo pedido por el cliente. [13]
2.1.6. Diseo De Un Sistema De Medicin Automtico En
Lnea De Produccin En Las Principales Estaciones De
Flujo Del La Unidad De Explotacin Tierra Oeste
(PDVSA)
Proyecto realizado en la Universidad de los Andes en Venezuela con el objetivo de
realizar una investigacin compuesta por un conjunto de instrumentos y equipos
de medicin, las cuales capturan el estado de las variables que intervienen en el
proceso de ujo del crudo desde el sitio de extraccin (Pozos) hasta los sitios de
recoleccin (estaciones de ujo), y por un sistema de supervisin que interpreta las
seales provenientes de estos dispositivos.
Esto se desarroll puesto que no exista un sistema de medicin automtica en la
Unidad de Explotacin Tierra Oeste de la empresa de Petrleos de Venezuela que
permitan conocer la cantidad y calidad del crudo recolectado en las principales
estaciones de ujo.
24
2.2 Marco terico
Para ello se disea un sistema que se ajusta a las normas tcnicas de scalizacin de
hidrocarburos lquidos, el cual permita obtener valores reales de los volmenes diarios
producidos ya que estos representan el pago de regalas de la nacin. Dicho diseo
est compuesto por un toma muestras automtico para determinar el porcentaje de
agua de la mezcla, un medidor de ujo, transmisores de presin y temperatura entre
otros elementos.
Bsicamente la diferencia en cuanto al proyecto que se trabaja es el tipo de uido,
ya que aqu se trabaja con crudo pesado, mientras que con el que se desea trabajar
son combustibles ligeros. [14]
2.2. Marco terico
2.2.1. HMI
En la interaccin del hombre con mquinas a travs de medios computacionales, hay
que tener en cuenta el trmino HMI. Trmino que hacer referencia a la supervisin
de procesos en su mayora industriales, mediante el uso de herramientas HMI
SCADA.[15]
HMI signica Human Machine Interface, es decir es el dispositivo o sistema que
permite el interfaz entre la persona y la mquina. Tradicionalmente estos sistemas
consistan en paneles compuestos por indicadores y comandos, tales como luces pilo-
tos, indicadores digitales y anlogos, registradores, pulsadores, selectores y otros que
se interconectaban con la mquina o proceso. En la actualidad, dado que las mqui-
nas y procesos en general estn implementadas con controladores y otros dispositivos
electrnicos que dejan disponibles puertas de comunicacin, es posible contar con sis-
temas de HMI bastantes ms poderosos y ecaces, adems de permitir una conexin
ms sencilla y econmica con el proceso o mquinas. [1]
- Tipos de HMI:
Descontando el mtodo tradicional, podemos distinguir bsicamente dos tipos de
HMIs:
- Terminal de Operador, consistente en un dispositivo, generalmente construido para
ser instalado en ambientes agresivos, donde pueden ser solamente de despliegues
numricos, o alfanumricos o grcos. Pueden ser adems con pantalla sensible al
tacto (touch screen).
- PC + Software, esto constituye otra alternativa basada en un PC en donde se
carga un software apropiado para la aplicacin. Como PC se puede utilizar cualquiera
segn lo exija el proyecto, en donde existen los llamados Industriales (para ambientes
agresivos), los de panel (Panel PC) que se instalan en gabinetes dando una apariencia
de terminal de operador, y en general veremos muchas formas de hacer un PC,
25
Captulo 2 Captulo 2
pasando por el tradicional PC de escritorio. Respecto a los softwares a instalar en
el PC de modo de cumplir la funcin de HMI.
Software HMI:
Estos softwares permiten entre otras cosas las siguientes funciones: Interface grca
de modo de poder ver el proceso e interactuar con l, registro en tiempo real e
histrico de datos, manejo de alarmas. Si bien es cierto slo con la primera funcin
enunciada es la propiamente HMI, casi todos los proveedores incluyen las otras
dos ya sea en el mismo paquete o bien como opcionales. Tambin es normal que
dispongan de muchas ms herramientas. Al igual que en los terminales de operador,
se requiere de una herramienta de diseo o desarrollo, la cual se usa para congurar la
aplicacin deseada, y luego debe quedar corriendo en el PC un software de ejecucin
(Run Time). Por otro lado, este software puede comunicarse directamente con los
dispositivos externos (proceso) o bien hacerlo a travs de un software especializado
en la comunicacin, cual es la tendencia actual.[1]
Figura 2.1.: Aplicacin HMI. [1]
Con los programas de diseo, se puede utilizar un editor de pantallas para desarrollar
toda la visualizacin de datos del proceso. Estos moldes son guardados y almacenan
la forma como sern visualizados los datos en las pantallas. Por lo tanto la interfaz
hombre-mquina se encarga de refrescar las variables de la base de datos en la
pantalla, adems de generar una base de datos de estas variables. [16]
2.2.2. Medidor de Coriolis
El principio de Coriolis es una ciencia exacta. La fuerza de Coriolis es la ms famosa
para los efectos de rotacin en la tierra, por ejemplo, causa vientos desde un rea
26
2.2 Marco terico
de alta presin hasta un sensor de vibracin de tubos espirales hacia el exterior en
una direccin hacia la derecha en el hemisferio norte.
Caudalmetros de Coriolis operan bajo los principios de la fuerza de Coriolis, en
donde se debe tener en cuenta la fuerza de inercia para describir los cuerpos en
movimiento en un plato giratorio. Si la masa M (en velocidad V) es guiado por la
placa giratoria (en la velocidad angular O). La fuerza de Coriolis (Fc) se ejerce sobre
la placa como se muestra en la Ecuacin 2.1:
Fc = 2 M O V (2.1)
Los medidores de ujo Coriolis tienen uno o ms tubos vibrantes como se muestra
en Figura 2.2. Cuando el lquido a ser medido pasa a travs de los tubos vibrantes,
se produce un movimiento de torsin en los tubos. El ujo de masa es proporcional
a los niveles de un movimiento de torsin (es decir, el movimiento angular). [2]
Figura 2.2.: Vibracin del tubo del sensor de Coriolis. [2]
Figura 2.3.: Torsin de los tubos del sensor de coriolis. [2]
2.2.3. PLC (Programing Logic Control)
Un controlador lgico programable (Programmable Logic ControllerPLC) Figura 2.4
es un dispositivo operado digitalmente, que usa una memoria para el almacenamiento
27
Captulo 2 Captulo 2
interno de instrucciones con el n de implementar funciones especcas, tales como
lgica, secuenciacin, registro y control de tiempos, conteo y operaciones aritmticas,
para controlar a travs de entradas/salidas digitales o analgicas, varios tipos de
mquinas o procesos. Los PLCs operan de manera secuencial y cclica, es decir,
una vez nalizado el recorrido completo de un programa, comienza a ejecutar su
primera instruccin.[17]
Figura 2.4.: PLC. [3]
Los elementos que contiene un PLC son:
Unidad Central de proceso Mdulos de entrada Mdulos de salida Fuente
de Alimentacin Dispositivos perifricos Interfaces
Un PLC o Autmata Programable posee las herramientas necesarias, tanto de soft-
ware como de hardware, para controlar dispositivos externos, recibir seales de sen-
sores y tomar decisiones de acuerdo a un programa que el usuario elabore segn el
esquema del proceso a controlar.
2.2.4. Sistema de Descarga de Hidrocarburos.
Se entiende como Produccin Fiscalizada de Hidrocarburos a los hidrocarburos
provenientes de determinada rea, producidos y medidos bajo trminos y condicio-
nes acordados en cada contrato. Esta medicin resulta importante para tres nes
principales:
- Registro histrico de las cifras ociales de la produccin scalizada por lote.
- Valorizacin de la produccin scalizada para efectos del pago/cobro de la rega-
la/retribucin por parte de los contratistas petroleros, lo que a su vez tiene efectos
sobre la recaudacin del estado por concepto de la renta petrolera.
- Facturacin por la compra-venta de hidrocarburos entre productores y usuarios de
los hidrocarburos.
28
2.2 Marco terico
LACT Lease: Automatic Custody Transfer Transferencia Automtica de Hidro-
carburos, con nes de compra/venta. Esta instalacin est conectada en la lnea de
descarga, para medir hidrocarburos en movimiento, en nuestro caso procedentes de
un tanque preparado para ser scalizado automticamente en calidad y cantidad
del hidrocarburo transferido; dispone de una unidad de rechazo para asegurar que
la calidad del producto est dentro de los rangos aceptables. Desde antes de su apa-
ricin y mientras no se instalen estn en mantenimiento estas unidades LACT, el
petrleo se scaliza por la medicin manual en volmenes y calidad en un tanque en
reposo, luego de abre la vlvula de descarga y se inicia la transferencia, cuando se
completa el volumen a entregarse, se cierra la vlvula y se espera un tiempo razona-
ble de reposo para proceder a la medida nal. Por diferencia se calcula el volumen
entregado, todo este proceso es manual.
La unidad LACT es una importante herramienta en la evolucin de la automatiza-
cin de transferencia de hidrocarburos, la cual se va mejorando con el tiempo y es
ampliamente aceptada en la Industria del Petrleo ya que esta normada por el API.
Su instalacin se justica para reemplazar un trabajo repetitivo hecho por uno o ms
medidores, que redunda en un ahorro de tiempo y eliminacin del error humano,
aunque cabe notar que la medicin automtica tampoco es exacta mientras que est
afectada por el error mecnico, que a diferencia del error humano, aquel puede ser
mejor controlado por una vericacin peridica en el mantenimiento y calibracin
del instrumento.
Para que estos instrumentos tengan la conabilidad exigida, tanto el comprador
como el vendedor de crudo, deben acreditar su exactitud con un certicado expedido
por el fabricante, basado en pruebas hechas en el laboratorio y campo, de acuerdo
a procedimientos estndar aprobados por el API.
2.2.5. GLP Gas Licuado de Petrleo
El GLP es un hidrocarburo que a condicin normal de presin y temperatura se en-
cuentra en estado gaseoso, sin embargo este en condiciones normales de temperatura
y moderadamente alta presin se encuentra en estado licuable.
El GLP puede ser propano o butano solamente, pero de igual menera puede ser
una mezcla de ambos hidrocarburos C3 y C4. La composicin del GLP puede va-
riar dependiendo la regin, por ejemplo en zonas donde la temperatura promedio
es alta la composicin se encuentra con un mayor porcentaje de C4 (abutanado),
mientras que en las zonas con menor temperatura promedio la composicin tiene
mayor porcentaje de C3 (apropanado). [18]
Las principales diferencias entre el propano y el butano son:
- El propano se vaporiza a temperaturas mayores de -44F (-42C) a presin atmos-
frica.
29
Captulo 2 Captulo 2
- A diferencia del propano, el butano no vaporiza adecuadamente a temperaturas
menores de 32F (0C).
- El propao tiene una mayor presin que el butano a una temperatura determinada.
- Cuando se produce una combustin, el propano produce menos calor en compara-
cin con la misma cantidad de gas butano.
- Un galn de gas propano pesa menos que un galn de gas butano.
Caractersticas del GLP:
- El GLP es un combustible limpio.
- No es txico, pero puede provocar asxia.
- Puede ocasionar irritaciones en contacto con la piel y con los ojos.
- Es altamente inamable, su combustin es muy rpida generando altas tempera-
turas.
- El GLP est compuesto, mayoritariamente, por propano y butano.
- El GLP se licua a bajas presiones entre 60 y 120 psi aproximadamente, dependiendo
de la mezcla propano butano.
- Posee una gran capacidad de expansin, de estado lquido a gaseoso aumenta su
volumen 270 veces aproximadamente.
- El GLP en estado gaseoso, es ms pesado que el aire, por ello, en caso de fugas
tiende a ubicarse o depositarse en lugares bajos.
- En estado lquido el GLP es ms liviano que el agua.
- El GLP es un combustible que en determinados porcentajes con el aire forma una
mezcla explosiva, presentando un Lmite de Inamabilidad para el propano entre
2.15 y 9.60 % de gas en aire, y para el butano, entre 1.55 y 8.60 % de gas en aire.
- El GLP producido de los lquidos de gas natural o de gases de renera es incoloro e
inodoro, por lo que para percibir su presencia en el ambiente se le aade un qumico
especial agente odorante denominado mercaptano.
2.2.6. Elementos Principales Del Proceso
2.2.6.1. PLC KOYO DL06
Este es tipo de PLC usado para el desarrollo del sistema de scalizacin y pertenece
a la familia de PLCs DL06 es una lnea de productos muy verstil que combina
caractersticas poderosas en un tamao compacto. Este PLC ofrece mdulos de
expansin de entradas y salidas discretas y anlogas, contadores de alta velocidad,
matemtica del punto otante, controladores PID, programacin de secuenciador de
tambor, varias opciones de comunicaciones seriales y un visor LCD.
30
2.2 Marco terico
Este PLC ofrece un conjunto de instrucciones muy similares a la poderosa CPU D2-
260 de la familia DirectLogic incluyendo instrucciones MODBUS y ASCII. Todos
los PLCs DL06 tienen dos puertos de comunicacin seriales que pueden ser usados
para programacin, para interfaces de operador, para redes, etc.
Las unidades con entradas de corriente continua tienen caractersticas de alta velo-
cidad en cuatro puntos de entrada. Las unidades con salidas de corriente continua
tienen pulsos de salidas de alta velocidad en diferentes formatos en el primero y
segundo punto de salida.
Figura 2.5.: PLC KOYO DL06.
2.2.6.2. Pantalla HMI MAPLE HMI5043N
Esta es la pantalla usada para desarrollar la interfaz hombre - mquina, es brillante
y de 4.3" TFT soporta hasta 65 mil colores y tiene una resolucin ms alta (480 x 272
pxeles) que la mayora de 6" pulgadas. Este modelo est diseado para aplicaciones
que requieren una interfaz de operador compacto resistente con baja potencia. El
HMI5043N opera en 24 V CC y tiene un modo de ahorro de energa que reduce la
corriente a la pantalla para ahorrar energa. Este terminal industrial requiere slo
1,5" de profundidad cuando se monta a un panel.
2.2.6.3. Impresora EPSON TM-T90
Esta impresora es la nica que se encuentra actualmente en el mercado y trabaja
con la pantalla HMI de MAPLE, motivo por el que fue seleccionada. Posee una lata
velocidad de impresin, auto corte de papel, cuchilla para corte manual de papel,
botn de alimentacin de papel.
31
Captulo 2 Captulo 2
Figura 2.6.: MAPLE HMI5043N.
Figura 2.7.: Impresora EPSON TM-T90.
2.2.7. NEMA 4X
NEMA es el acrnimo de Asociacin Nacional de fabricantes de electrnica (por sus
siglas en ingls: National Electrical Manufacturers Association). Esta organizacin
ofrece gabinetes de componentes elctricos en calicaciones, como NEMA 4X, basa-
das en sus cualidades protectoras. NEMA especica si un gabinete de componentes
elctricos es seguro para su uso en exteriores, y qu tipo de condiciones ambientales
pueden soportar.
El gabinete fabricado para el proceso es de tipo NEMA 4X y este dicta que son
cajas fabricadas para el uso de interior o al aire libre para proporcionar un grado
proteccin al personal contra el acceso a partes peligrosas, para proporcionar un
grado de proteccin del equipo en el interior de la envolvente contra la penetracin
de cuerpos slidos extraos (polvo arrastrado por el viento), para proporcionar un
grado de proteccin con respecto a los efectos nocivos en el equipo debido a la entrada
de agua (lluvia, aguanieve, nieve, agua salpica y manguera de agua dirigido), que
proporciona un nivel adicional de proteccin contra la corrosin y no se daan por
32
2.2 Marco terico
la formacin externa de hielo en el gabinete. [19]
Los otros tipos de cajas o gabinetes que manejan la NEMA son los siguientes:
Tipo 1: Cajas fabricadas para uso en interiores para proporcionar un grado de protec-
cin al personal contra el acceso a partes peligrosas y para proporcionar un grado de
proteccin del equipo en el interior de la envolvente contra la penetracin de cuerpos
slidos extraos (suciedad que cae).
Tipo 2: Cajas fabricadas para uso en interiores para proporcionar un grado de pro-
teccin al personal contra el acceso a partes peligrosas, para proporcionar un grado
de proteccin del equipo en el interior de la envolvente contra la penetracin de
cuerpos slidos extraos (suciedad que cae), y para proporcionar un grado de pro-
teccin con res pecto a los efectos nocivos en el equipo debido a la entrada de agua
(goteo y salpicaduras ligeras).
Tipo 3: Cajas fabricadas para el uso de interior o al aire libre para proporcionar un
grado de proteccin al personal contra el acceso a partes peligrosas, para propor-
cionar un grado de proteccin del equipo en el interior de la envolvente contra la
penetracin de cuerpos slidos extraos (cada de suciedad y el polvo arrastrado por
el viento), a disposicin de un grado de proteccin con respecto a los efectos nocivos
en el equipo debido a la entrada de agua (lluvia, aguanieve, nieve); y no se daan
por la formacin externa de hielo en el gabinete.
Tipo 3R: Cajas fabricadas para el uso de interior o al aire libre para proporcionar
un grado de proteccin al personal contra el acceso a partes peligrosas, para pro-
porcionar un grado de proteccin del equipo en el interior de la envolvente contra
la penetracin de cuerpos slidos extraos (suciedad que cae), para proporcionar
un grado de proteccin con respecto a los efectos nocivos en el equipo debido a la
entrada de agua (lluvia, aguanieve, nieve), y no se daan por la formacin externa
de hielo en el gabinete.
Tipo 3S: Cajas fabricadas para el uso de interior o al aire libre para proporcionar
la penetracin de cuerpos slidos extraos (cada de suciedad y el polvo arrastrado
por el viento), a disposicin de un grado de proteccin con respecto a los efectos
nocivos en el equipo debido a la entrada de agua (lluvia, aguanieve, nieve), y para
los cuales el mecanismo externo (s) permanezca operable cuando el hielo cargado.
Tipo 3X: Cajas fabricadas para el uso de interior o al aire libre para proporcionar
por el viento), para proporcionar un grado de proteccin con respecto a los efectos
nocivos en el equipo debido a la entrada de agua (lluvia, aguanieve, nieve), que
proporciona un nivel adicional de proteccin contra la corrosin y no se daan por
la formacin externa de hielo en el gabinete.
33
Captulo 2 Captulo 2
Tipo 3RX: Cajas fabricadas para el uso de interior o al aire libre para proporcionar
entrada de agua (lluvia, aguanieve, nieve), que se daan por la formacin externa
de hielo en el envolvente que proporciona un nivel adicional de proteccin contra la
corrosin, y que habr daos por la formacin externa de hielo en el gabinete.
Tipo 3SX: Cajas fabricadas para el uso de interior o al aire libre para proporcionar
nocivos en el equipo debido a la entrada de agua (lluvia, aguanieve, nieve), que
proporciona un nivel adicional de proteccin contra la corrosin, y para los cuales
el mecanismo externo (s) permanezca operable cuando el hielo cargado .
Tipo 4: Cajas fabricadas para el uso de interior o al aire libre para proporcionar
nocivos en el equipo debido a la entrada de agua (lluvia, aguanieve, nieve, agua,
salpicaduras, y manguera de agua dirigido), y no se daan por la formacin externa
de hielo en el gabinete.
Tipo 5: Cajas fabricadas para uso en interiores para proporcionar un grado de pro-
teccin al personal contra el acceso a partes peligrosas, para proporcionar un grado
de proteccin del equipo en el interior de la envolvente contra la penetracin de
cuerpos slidos extraos (suciedad que cae y liquidacin de polvo en el aire, pe-
lusa, bras y yings), y para proporcionar un grado de proteccin con respecto a
los efectos nocivos en el equipo debido a la entrada de agua (goteo y salpicaduras
ligeras).
Tipo 6: Cajas fabricadas para el uso de interior o al aire libre para proporcionar
entrada de agua (manguera dirigida de agua y la entrada de agua de vez en cuando
durante la inmersin temporal a una profundidad limitada), y que sern daos por
la formacin externa de hielo en el gabinete.
Tipo 6P: Cajas fabricadas para el uso de interior o al aire libre para proporcionar
34
2.2 Marco terico
entrada de agua (manguera de agua dirigida y la entrada de agua durante la inmer-
sin prolongada a una profundidad limitada), que proporciona un nivel adicional de
proteccin contra la corrosin y no se daan por el exterior formacin de hielo en el
gabinete.
Tipo 12: Cajas fabricadas (sin nocauts) para su uso en interiores para proporcionar
un grado de proteccin al personal contra un l acceso a partes peligrosas, para
proporcionar un grado de proteccin del equipo en el interior de la envolvente contra
la penetracin de cuerpos slidos extraos (suciedad que cae y la circulacin de polvo,
pelusas , bras y partculas voltiles), y para proporcionar un grado de proteccin
con respecto a los efectos nocivos en el equipo debido a la entrada de agua (goteo y
salpicaduras ligeras).
Tipo 12K: Cajas fabricadas (con nocauts) para su uso en interiores para propor-
cionar un grado de proteccin al personal contra el acceso a partes peligrosas, para
proporcionar un grado de proteccin del equipo en el interior de la envolvente con-
tra la penetracin de cuerpos slidos extraos (suciedad que cae y la circulacin
de polvo, pelusas , bras y partculas voltiles), y para proporcionar un grado de
proteccin con respecto a los efectos nocivos en el equipo debido a la entrada de
agua (goteo y salpicaduras ligeras).
Tipo 13: Cajas fabricadas para uso en interiores para proporcionar un grado de
proteccin al personal contra el acceso a partes peligrosas, para proporcionar un
grado de proteccin del equipo en el interior de la envolvente contra la penetracin
de cuerpos slidos extraos (suciedad que cae y la circulacin de polvo, pelusa, bras
y yings ), para proporcionar un grado de proteccin con respecto a los efectos
nocivos sobre el equipo debido a la entrada de agua (goteo y salpicaduras ligeras),
y para proporcionar un grado de proteccin contra las fumigaciones, salpicaduras,
y la ltracin de petrleo y no corrosivo refrigerantes.
En reas peligrosas, cuando se instala y se mantiene completa y adecuadamente,
Tipo 7 y 10 cajas estn diseadas para contener una explosin interna sin causar
un peligro externo. Tipo 8 cajas estn diseadas para evitar la combustin a travs
del uso de equipos en bao de aceite. Tipo 9 recintos estn diseados para evitar la
ignicin de polvo combustible.
Tipo 7: Cajas fabricadas para uso en interiores en lugares peligrosos (clasicados)
clasicados como Clase I, Divisin 1, Grupos A, B, C, o D, segn se dene en la
norma NFPA 70.
Tipo 8: Cajas fabricadas para el uso de interior o al aire libre en lugares peligrosos
(clasicados) clasicados como Clase I, Divisin 1, Grupos A, B, C, y D como se
dene en la norma NFPA 70.
Tipo 9: Cajas fabricadas para uso en interiores en lugares peligrosos (clasicados)
clasicados como Clase II, Divisin 1, Grupos E, F o G como se dene en la norma
NFPA 70.
35
Captulo 2 Captulo 2
Tipo 10: Cajas fabricadas para cumplir con los requisitos de la Administracin de
Seguridad y Salud, 30 CFR, Parte 18.
La Tabla 2.1 provee informacin sobre el tipo de reas clasicadas y los grupos
principales de los cuales se basa la NEMA para la fabricacin de las cajas.
Cuadro 2.1.: Clasicacin reas restringidas
2.2.8. RS-485
RS-485 o tambin conocido como EIA-485, que lleva el nombre del comit que lo
convirti en estndar en 1983. Es un estndar de comunicaciones en bus de la capa
fsica del Modelo OSI.
Est denido como un sistema en bus de transmisin multipunto diferencial, es
ideal para transmitir a altas velocidades sobre largas distancias (35 Mbit/s hasta 10
metros y 100 kbit/s en 1200 metros) y a travs de canales ruidosos, ya que reduce
los ruidos que aparecen en los voltajes producidos en la lnea de transmisin. El
medio fsico de transmisin es un par entrelazado que admite hasta 32 estaciones en
1 solo hilo, con una longitud mxima de 1200 metros operando entre 300 y 19200
bit/s y la comunicacin half-duplex (semiduplex). Soporta 32 transmisiones y 32
receptores. La transmisin diferencial permite mltiples drivers dando la posibilidad
de una conguracin multipunto. Al tratarse de un estndar bastante abierto permite
muchas y muy diferentes conguraciones y utilizaciones.
Desde 2003 est siendo administrado por la Telecommunications Industry Associa-
tion (TIA) y titulado como TIA-485-A.222. [20]
Especicaciones
Interfaz diferencial Conexin multipunto Alimentacin nica de +5V Hasta 32 es-
taciones (ya existen interfaces que permiten conectar 256 estaciones) Velocidad m-
36
2.2 Marco terico
xima de 10 Mbit/s (a 12 metros) Longitud mxima de alcance de 1200 metros (a
100 kbit/s) Rango de bus de -7V a +12V
Aplicaciones
SCSI -2 y SCSI-3 usan esta especicacin para ejecutar la capa fsica. RS-485 se usa
con frecuencia en las UARTs para comunicaciones de datos de poca velocidad en las
cabinas de los aviones.
Algunas unidades de control del pasajero lo utilizan, equipos de monitoreo de siste-
mas fotovoltaicos.
Requiere el cableado mnimo, y puede compartir el cableado entre varios asientos.
Por lo tanto reduce el peso del sistema.
RS-485 se utiliza en sistemas grandes de sonido, como los conciertos de msica y
las producciones de teatro, se usa software especial para controlar remotamente
el equipo de sonido de una computadora, es utilizado ms generalmente para los
micrfonos.
RS-485 tambin se utiliza en la automatizacin de los edicios pues el cableado
simple del bus y la longitud de cable es larga por lo que son ideales para ensamblar
los dispositivos que se encuentran alejados.
RS-485 Tiene la mayor parte de su aplicacin en las plantas de produccin automa-
tizadas.
37
3. Captulo 3
Este captulo aborda el desarrollo del sistema totalizador donde se empieza descri-
biendo la metodologa empleada seguida de una descripcin de los instrumentos,
funcionamiento e instalacin del sistema. De igual manera se muestra el sistema ini-
cial desarrollado y el montaje realizado dentro del laboratorio de la empresa ATLAS
para las pruebas y demostraciones para la empresa.
3.1. Metodologa
La metodologa de solucin fue desarrollada de acuerdo a requerimientos que deba
tener el sistema junto a las caractersticas que poseen el sensor y transmisor de ujo.
Con las caractersticas del sistema se planteo la siguiente solucin:
Figura 3.1.: Diagrama general para la solucin del sistema
39
Captulo 3 Captulo 3
El sistema que exista anteriormente utilizaba el sensor RHM 20 y el transmisor RHE
08. Estos mismos elementos se usaron para la realizacin del sistema totalizador, por
lo que se tuvo en cuenta las caractersticas del transmisor para poder seleccionar
un PLC apropiado. EL transmisor RHE 08 entrega 1 salida de pulsos y 2 seales
analgicas 4-20mA las cuales deban ser ledas apropiadamente por el PLC, las
caractersticas de estas seales se muestran a continuacin.
3.1.1. Seal de pulsos
EL transmisor de ujo RHE 08 totaliza la masa o volumen que pasa por el sensor
de coriolis, y entrega el valor totalizado en una seal de pulsos la cual puede ser
congurada para que tenga una frecuencia mxima de 10kHz y entrega una cier-
ta cantidad de pulsos sobre unidad. La cantidad de pulsos puede ser congurada
mediante un jumper el cual segn su posicin realiza una divisin por decada, con
esta conguracin y segn el ujo de trabajo se elije una opcin de pulsos, como lo
muestra la Tabla 3.1.
Cuadro 3.1.: Eleccin de cantidad de pulsos sobre unidad[4].
La salida del transmisor para la seal de pulsos es de colector abierto, como lo
muestra Figura 3.2, esta conguracin necesita de una fuente exterior de 24 VDC y
una resistencia para regular el voltaje de la seal.
3.1.2. Seal 4-20mA
El transmisor posee ademas 2 salidas analgicas 4-20mA las cuales pueden ser con-
guradas para transmitir las variables: ujo msico, densidad o temperatura. El
transmisor posee 3 modos de operacin para las salidas 4-20mA[4]:
1. 4-20mA: El rango de la seal esta entre 4 y 20mA. Valor analgico de estado
de error igual a 2mA.
40
3.1 Metodologa
Figura 3.2.: Conguracin de salida de pulsos del transmisor[4].
2. 3.7-20mA: El rango de la seal esta entre 3.7 y 20mA. Valor analgico de
estado de error es menor a 3.7mA.
3. 4-22mA: El rango de la seal esta entre 4 y 22mA. Valor analgico de estado
de error mayor a 22mA.
3.1.3. PLC KOYO DL06
De acuerdo a las seales de el transmisor se desidi por utilizar el PLC de referencia
D0-06DD2 el cual posee 20 entradas sinking y source, 16 salidas de tipo Source,
ademas cuenta con 2 entradas rpidas HSIO las cuales se usaron para leer los pulsos
que entregaba el transmisor de ujo. El PLC no posee entradas analgicas pero
se puede expandir con diversos mdulos Figura 3.3. Para las seales de 4-20mA
se adquiri un mdulo analgico de referencia F0-04AD-1 el cual posee 4 canales
analgicos. Ademas cuanta con 2 puertos de comunicaciones los cuales manejan los
protocolos RS 232, RS-485, RS-422, protocolo K y ASCII. Este PLC entrega un
voltaje de 24 VDC el cual pede ser utilizado para alimentar dispositivos externos al
PLC y tambin para poder proporcionar un voltaje necesario para la conguracin
de entradas y salidas como lo muestra Figura 3.4
41
Captulo 3 Captulo 3
Figura 3.3.: PLC Koyo DL06[5].
Figura 3.4.: Conguracin de entradas y salidas del PLC[5].
42
3.1 Metodologa
3.1.3.1. Entradas y salidas sinking/Sourcing (Drenadora/surtidora)
Las entradas y salidas drenadoras hacen referencia a la ruta para suministrar el
negativo (-), las surtidoras hacen referencia al positivo (+). La terminologa de
drenadoras y surtidoras se aplica slo a circuitos de entradas y salidas de corriente
continua (CC). Los puntos de entradas y las salidas que son drenadoras o surtidoras
pueden conducir corriente en una direccin solamente. Esto signica que es posible
conectar el aparato externo de suministro y el campo al punto de entradas y salidas
de tal forma que la corriente uya en la direccin errada y el circuito no operar.[5].
Las conguraciones de este tipo de entradas y salidas se muestra en Figura 3.5.
Figura 3.5.: Conguracin de entradas y salidas drenadoras/surtidoras[5]
3.1.3.2. Mdulo analgico
El mdulo analgico F0-04AD-1 posee 4 canales analgicos de corriente para seales
de 0/4-20mA con una resolucin de 12 bit. Para realizar la conexin a dispositivos
de campo es necesario una fuente de alimentacin externa, pero la conguracin del
transmisor RHE 08 hace que no sea necesario utilizar una fuente adicioneal, pues
este entrega la seal ya lista para ser utilizada.
3.1.4. HMI
La pantalla HMI utilizada es la Mapple 5043N la cual maneja el protocolo K-secuence
el cual es soportado por el PLC y cuenta con las opciones de guardar los datos en
una base de datos, impresin de pantalla en varias impresoras listadas en el manual,
posee ademas proteccin Nema 4X. El HMI posee 2 puertos en los que se puede
utilizar los protocolos de comunicacin RS-232 y RS-485 de 2 y 4 hilos
43
Captulo 3 Captulo 3
Figura 3.6.: Mdulo analgico F0-04AD-1[6].
3.1.5. Impresora EPSON TM-T90
En la lista de las opciones de impresoras compatibles con el HMI esta la EPSON
TM-T90 la cual se adquiri, la EPSON TM-T90 es una impresora para tickets, posee
la funcin de auto-cortado de papel y maneja el protocolo RS-232
3.2. Desarrollo
En este apartado se mostrara el desarrollo de el sistema y cumplimiento de los pri-
meros objetivos propuestos, ademas de otras funciones que se adicionaron al sistema.
3.2.1. PRUEBAS INICIALES
Al principio, se realizaron pruebas con un gabinte inicial ver Figura 3.9 y se realiz
una demostracin a la empresa ATLAS E&A para que ellos realizan las repectivas
observaciones y se hicieran las correcciones necesarias, el montaje para la prueba
el cual se observa en Figura 3.7 se hizo en el laboratorio de la misma empresa.
Este montaje contaba con el sensor coriolis instalado en un circuito cerrado para la
circulacin de agua la cual era movida por una bomba que se controlaba por medio
de una fuente variable.
44
3.2 Desarrollo
Figura 3.7.: Montaje del sistema en el laboratorio ATLAS
El siguiente plano muestra las conexones y tipos de seales que se manejaban en el
montaje realizado en el laboratorio de ATLAS
Figura 3.8.: Plano P&ID Laboratorio ATLAS
Para esa prueba se hicieron las respectivas conexiones entre el sensor de coriolis, el
transmisor y el sistema totalizador. El sistema totalizador estaba elaborado dentro
45
Captulo 3 Captulo 3
de un gabinete sencillo el cual se muestra en Figura 3.9 el cual era para aplicaciones
sencillas y solo soportaba las condiciones ambientales ms bsicas.
Figura 3.9.: Gabinete del totalizador inicial
Figura 3.10.: Vista interna del gabinete inicial
En la prueba y demostracin del funcionamiento del sistema totalizador los resulta-
dos fueron favorables, en donde el sistema realiz la medida de masa comparando la
46
3.2 Desarrollo
medida del totalizador con la medida mostrada por el transmisor y se observo una
completa igualdad sin error. De igual manera las funciones instaladas al totaliza-
dor como las conguraciones y calibraciones funcionaron sin problema alguno, sin
embargo se realizaron las siguientes observaciones para mejorar el sistema:
- La ms importante fue el cambio necesario del gabinete puesto que el rea de
instalacin es clasicada y requeria como mnimo un gabinete con norma NEMA 4X
- Adicionar funciones de calibracin y conguracin (DELTA, Nmero de bombas,
Totalizador por volumen, Unidades de medida).
- Mejorar el registro de datos del sistema.
- Adicionar datos a los tiquetes impresos (Nmero de placa, nombre de la empresa,
tipo de producto)
-Agregar funciones al modo de totalizacin como el de un PRESET mltiple y
visualizacin de otras variables del proceso
3.2.2. TOTALIZADOR
Como se menciono anteriormente el transmisor RHE 08 totaliza por masa o por vo-
lumen segn la eleccin que se desee. En la empresa Pegasus Blending se desarrolla
el proceso scal mediante masa, por lo que principalmente se desarrollo un sistema
totalizador de masa, aunque el sistema se desarrollo para poder hacer un proceso s-
cal tambin por volumen
1
, El totalizador desarrollado en el PLC y HMI se realizo de
forma que las conguraciones que se tuvieran en el transmisor respecto ala variable
totalizada, el sistema las obtuviera tambin, principalmente se trabajo primero con
la entrada rpida del PLC, las entradas HSIO, en donde se desarrollo la conexin
mostrada en Figura 3.2, este valor se llevo a la entrada de un contador del PLC el
cual tiene la caracterstica de contar hasta 99.999.999 pulsos. EL PLC tiene puede
leer entrada de pulsos los cuales posean una frecuencia hasta 7KHz, para probar
si el conteo realizado era el correcto se utilizo un circuito el cual entregaba 10mil
pulsos a una frecuencia de 5KHz. Con los pulsos ledos correctamente se prosedia a
realizar la divisin de pulsos (Pulsos/ Kilogramo) en donde simplemente se divida
el nmero de pulsos ledo por el PLC con el valor seleccionado.
Pulsos ledos
Pulsos/Kg
= Masa(Kg) (3.1)
Se desarrollaron 3 formas para poder introducir este valor de pulsos/Unidad:
Mediante una tabla: Se desarrollo una interfaz en el HMI la cual tuviera la
misma forma que la Tabla 3.1 la cual aparece en la hoja de datos de el trans-
misor. En la pantalla al presionar sobre el valor deseado, el HMI envia este
valor al PLC para poder hacer la divisin.
47
Captulo 3 Captulo 3
Figura 3.11.: Seleccin de Pulsos/Kg, mediante tabla
K-factor introducido: se desarrollo una interfaz como lo muestra la Figura 3.12
en donde existen 2 formas de ingresar el K-factor el cual es el valor de divisin
de pulsos, esta forma se realizo por si al utilizar otro transmisor de ujo este no
tuviera las mismas opciones que en la tabla anterior. en este modo el operario
ingresa el valor del K-factor el cual puede ser un nmero real con punto otante.
Figura 3.12.: Ajuste del K-factor
K-factor ajustado: en la misma pantalla que muestra la Figura 3.12 se puede
observar que esta ventana muestra la conguracin actual de pulsos/Kg y
debajo de esta hay un recuadro en donde se muestra la cantidad actual de masa
calculada. Si se posee un patrn o para un ajuste rpido se puede ingresar el
valor de masa que debera marcar el totalizador y as el PLC se encargara de
calcular el factor apropiador y ajustara el valor al que debera mostrar.
1
La medicin de volumen entregada es la misma medicin que el valor totalizado en el transmisor
48
3.2 Desarrollo
3.2.2.1. Unidades de medida
Si se desea totalizar por volumen, con el transmisor congurado para que totalice
por volumen se puede hacer esta opcin en el men del HMI como lo muestra la
Figura 3.13.
Figura 3.13.: Eleccin de totalizar masa/volumen
Adems de seleccionar como se va a totalizar tambin se creo un interfaz en donde
el operario puede elegir en que unidades de medida desea totalizar, estas unidades
son: para masa; g (gramo), kg (Kilogramo), t (toneladas) y lb (libras). Volumen; L
(litros), m
3
(metros cbicos), ml (mililitros), gal (galn americano).
Para la conversin de unidades se utilizaron las equivalencias:
1kg=1000g
1000kg=1t
1 lb =0.4536 kg
1 L= 0.001m
3
1 L= 1000ml
1 galn americano (gal) es igual a 3.7854L
Las conversiones se hacen desde el Kilogramo o desde el Litro, donde estas variables
son las que se obtienen directamente desde el transmisor despues de hacer la divisin
de pulsos.
Al cambiar las unidades en que se totaliza, la funcin de divisin por pulsos no
cambia, esta tiene una unidad base la cual es kg (kilogramo) o L (litros)
Cada ves que el totalizador realiza un nuevo proceso de carga en el carro tanque,
guarda los datos del proceso como son la cantidad descargada, placa del vehculo,
49
Captulo 3 Captulo 3
Figura 3.14.: Unidades de medida para el totalizador.
y la empresa a la cual se descargo el producto. esto se describe detalladamente mas
adelante.
3.2.3. FUNCIONES DEL PROGRAMA DEL HMI
El programa del sistema totalizador se muestra a continuacin en donde se muestra
el programa desarrollado en el HMI el cual fue realizado en el software Easy Builder
5000.
EL programa del HMI es el programa principal del sistema, pues este es el que
va a utilizar el operario y es el que enva las ordenes al PLC para poder realizar
diversas acciones, esta conformado por diversas funciones las cuales son mostradas
en el men de la pantalla que se puede ver en la Figura 3.15. El men esta dividido
en varias partes como son: PRODUCTO (donde se realiza la descarga del producto),
REGISTRO (la base de datos del sistema), CONFIGURACIN (conguraciones que
se pueden realizar para que el sistema funcione de forma ptima), CALIBRACIN
(calibracin de la entradas del sistema), ABOUT (informacin del sistema). Tambin
se puede observar los logos de las empresas PEGASUS BLENDING international y
ATLAS E&A S.A.S.
3.2.3.1. PRODUCTO
Al ingresar en la primera opcin del men principal observamos la pantalla mostra-
da en la Figura 3.16. Esta pantalla a su ves posee diversas funciones como son la
visualizacin de la cantidad actual, visualizacin de estado de los procesos, visuali-
zacin de la fecha y hora, botones de control, ingreso de datos del ticket, ingreso del
valor del Preset.
50
3.2 Desarrollo
Figura 3.15.: Pantalla Principal del programa
En la visualizacin del producto actual se observa la cantidad totalizada junto a las
unidades en que se totaliza. Si se toca la pantalla sobre el recuadro de cantidad,
aparecern otras variables del proceso, estas variables pueden ser vistas en la panta-
lla mostrada en la Figura 3.17 en donde se muestran la Masa, Volumen y Densidad,
tambin se puede mostrar otra variable la cual puede ser temperatura o ujo, depen-
diendo de la conguracin realizada tanto en el programa como en el Transmisor.
Se observa tambin los estados de 3 actuadores: Vlvula, Bomba1 y Bomba 2. Estos
actuadores son los encargados de hacer circular el producto por el sistema para po-
der realizar la descarga. segn el modo de totalizacin seleccionado (masa/volumen)
el sistema calcula la otra variable por medio de la lectura de la densidad la cual se
lee en kg/l y la ecuacin = m/V esta variable es transmitida por el transmisor
mediante un canal 4-20mA. El valor del volumen calculado no esta compensado,
osea no esta llevado a una temperatura base como se aria en un proceso scal, esto
se muestra mas adelante en la correccin de volumen mediante la tabla 54E de API.
Las unidades de masa y volumen son las mismas que se seleccionan en el submenu
CONFIGURACIN.
Al volver a la pantalla anterior, la pantalla PRODUCTO observamos unos indi-
cadores :
VLVULA: muestra el estado de la vlvula y bombas instaladas.
MANUAL: Si esta prendido indica que el modo de operacin del sistema esta
en manual, si esta apagado esta en modo automtico. En modo automtico
cada vez que se realiza un proceso de descarga el operario deber ingresar un
valor de PRESET al sistema, este valor de preset es el valor que se entregara en
la descarga, el sistema prender y apagara los diversos actuadores que hacen
circular el producto. Cuando el sistema alcance el valor deseado apagara los
elementos y esperara a la orden de impresin de recibo o se puede realizar otra
descarga por si el el carro-tanque no se alcanzo a cargar a la medida deseada.
51
Captulo 3 Captulo 3
Figura 3.16.: pantalla PRODUCTO
Cuando se da la orden de impresin al sistema, este imprime un recibo de
la descarga y guarda el dato en la base de datos. Si el sistema esta en modo
manual, el operario deber prender y apagar los actuadores de control mediante
el botn del HMI de INICIO, tambin deber reiniciar la cuenta manualmente
y guardar el dato mediante el botn de guardado.
GUARDADO: cuanto este indicador se enciende signica que el sistema ha
guardado un dato.
Se tienen tambin los botones de control:
INICIO : Es el botn de color verde, en modo automtico inicia el proceso
de descargue. En modo manual prende y paga los actuadores de control.
RESET: Botn de color azul, este botn reinicia el contador, solo funciona
si el sistema se encuentra en modo manual.
REINICIAR PRESET: Botn en forma de X, cuando se pulsa vuelve el
valor del preset a 0.
GUARDAR: En modo manual el sistema guarda el dato actual. En modo
automtico no realiza ninguna accin.
IMPRIMIR: Cuando es pulsado, el sistema imprime un recibo el cual posee
los datos actuales del recibo como son: la cantidad entregada, la hora, la fecha,
el nmero del ticket e informacin adicional del ticket la cual puede ser selec-
cionada en donde hay 2 formatos de ticket. Uno de los 2 formato se observa
en la Figura 3.18.
Datos para ingresar:
52
3.2 Desarrollo
Figura 3.17.: Datos de medida
Figura 3.18.: Formato de Ticket numero 1
PRESET: Este dato funciona en modo automtico, y es el valor deseado de la
descarga. Este dato puede ser un dato real con punto otante.
Placa: cuando la opcin 1 del ticket esta habilitada es necesario ingresar este
dato. Este dato como su nombre lo indica es el nmero de la placa del vehculo
al cual se le ha realizado la descarga.
Producto: Cuando la opcin 1 del ticket o recibo es seleccionada se activa esta
casilla en la cual se ingresa el tipo de producto a Fiscalizar.
Empresa: con la opcin 1 de ticket, se habilita la casilla empresa, en esta casilla
se introduce el nombre de la empresa a la cual se le entrega el producto.
Precio: Cuando la opcin 2 del ticket es seleccionada esta opcin se habilita.
El valor que se ingresa en este dato es el valor de precio/Unidad.
IVA: es el valor del impuesto que se cobra sobre el producto. Este Impuesto
se aplica al total del producto entregado y es ingresado en porcentaje.
53
Captulo 3 Captulo 3
Figura 3.19.: Ticket tipo 2
3.2.3.2. REGISTRO
Al entrar el submenu REGISTRO se observa la pantalla mostrada en la Figura 3.20
en donde se muestran los registros de las descargas cada historial de la base de datos
es mostrada por das, as si se desea visualizar un da diferente al mostrado se debe
oprimir en la parte superior izquierda donde muestra un 0 y escribir el nmero de
das contando desde el da actual al cual se quiera consultar. El da actual es mostra-
do por el nmero 0, si se desea consultar el historial del da anterior, se debe escribir
el nmero 1 y si se desea ver el historial de hace 20 das se escribir el nmero 20.
en el caso de no encontrar la fecha deseada se puede ver en el historial a que fecha
corresponde cada nmero y as encontrar ms rpido el historial deseado. Este base
de datos est congurada para que guarde registros con un tiempo igual a 365 das
(1 ao) pero si en ese tiempo se excede la capacidad de memoria del HMI los datos
se irn sobrescribiendo y se irn borrando los datos ms antiguos.
Figura 3.20.: Pantalla del submenu REGISTRO
54
3.2 Desarrollo
Los datos guardados en la base de datos son los que aparecen en el recibo; hora, fecha,
cantidad, unidad de medida, placa, empresa. Para poder llevar un control de las
descargas realizadas en el da, se incorporaron 2 botones en la pantalla REGISTRO
los cuales son de impresin. Si se oprime el primer botn, el que se encuentra en
la parte inferior izquierda, el sistema imprimir lo que se encuentra en la pantalla
actual (se realiza una impresin de la pantalla) como lo muestra la Figura 3.21 en
donde se observan 2 registros que se encuentran en el historial.
Figura 3.21.: Opcin IMPRIMIR PANTALLA
La otra opcin de impresin se llama IMPRIMIR TOTAL en donde al oprimirse, el
sistema imprimir el total que se ha descargado el da actual, se imprime tambin el
nmero de descargas que se llevan en el da, tambin se imprimen los datos de fecha
y hora. cada ves que se inicia un nuevo da, el sistema reinicia la cuenta de total y
tambin el del nmero de descargas. En la Figura 3.22 se muestra el recibo asociado
a el historial mostrado en la Figura 3.21, en donde se llevan en el da 2 descargas
realizadas y un total de 0.87Kg.
Figura 3.22.: Impresin del total del da
El recibo posee los datos de la masa totalizada y el volumen totalizado. este dato
depende de como este congurado el sistema totalizador y el transmisor. En el recibo
de la Figura 3.22 se muestra que aparece 0 gal, esto es debido a que durante el da
no se ha totalizado por volumen.
55
Captulo 3 Captulo 3
3.2.3.3. CONFIGURACIN
En el men conguracin abreviado como CONFIG se encuentran diversas opcio-
nes, ver Figura 3.24, como son: Manual (eleccin del sistema en modo manual o
automtico), Totalizador (masa o volumen), Unidades de medida, Fecha y hora, da-
tos de recibo (tipo de recibo), ingresar password. Para poder ingresar a todas estas
opciones el operario necesita ingresar un password, si no se ingresa la contrasea
y se oprime alguna opcin de conguracin se obtendr el mensaje mostrado en la
Figura 3.23.
Figura 3.23.: Mensaje: Elemento protegido por contrasea
Figura 3.24.: submenu CONFIG
Password
Si se desea ingresar a alguna de la opciones que el submenu CONFIG posee es ne-
cesario ingresar un password, pues solo personas autorizadas pueden cambiar estas
opciones de conguracin. para poder ingresar la contrasea es necesario presionar
el icono en donde aparece la pantalla mostrada en la Figura 3.25 en donde el ope-
rario debe ingresar la contrasea, si la contrasea es incorrecta la pantalla indicara
que se ha introducido una clave incorrecta. una vez introducida esta contrasea, el
operario ademas de poder ingresar a todas las conguraciones que posee el sistema,
tambin puede cambiar la contrasea introducida por una nueva.
56
3.2 Desarrollo
Figura 3.25.: Password
Al pasar 15 minutos de introducida la contrasea y si el operario no realiza ninguna
accin, el sistema se volver a bloquear para no quedar abierto a cualquier personal.
El sistema tambin posee un botn de LOG OUT para poder bloquear el ingreso a
las conguraciones del sistema.
Manual
Con esta opcin se puede seleccionar si el sistema operara de modo manual o auto-
mtico. La diferencia de estos modos radca que en el modo automtico el operario
debe ingresar un valor para el PRESET, despus dar inicio y el sistema prende-
r y apagara los elementos de control (vlvulas y bombas). En el modo manual el
operario prender y apagara los actuadores en el momento que el desee. En el mo-
do automtico, despus de realizar una impresin el sistema guarda un registro en
el historial, mientras que en el modo manual el operario debe guardar el registro
mediante el botn de guardado.
Totalizador
En este men el operario elije la forma en que se va a totalizar; por masa o por
volumen. As como lo muestra la Figura 3.13.
Unidad de medida
Al seleccionar esta opcin el HMI muestra en pantalla la ventana mostrada en la
Figura 3.26. En esta pantalla se muestran las variables que son medidas por el
sistema: Masa, Volumen, densidad y otra variable la cual puede ser temperatura o
ujo msico. si se presiona el botn UNIDADES aparece en pantalla la ventana
mostrada en la Figura 3.14 de donde se seleccionan las unidades en que se va a
totalizar la masa y el volumen.
57
Captulo 3 Captulo 3
Figura 3.26.: Unidad de medida
Si se oprime el botn 4-20mA en pantalla se muestra la ventana de la Figura 3.27
Figura 3.27.: Canales de 4-20mA
En esta ventana se muestra las conguraciones para los canales analgicos. En el
canal 1 se tiene la variable densidad, esta variables es usada para poder calcular ya
sea la masa o volumen dependiendo de la conguracin que se tenga en el transmi-
sor. L signica el nivel bajo de la seal, osea el valor que corresponde a 4mA, H
es el valor alto de la seal el cual corresponde al valor de corriente de 20mA, D es
el valor actual de la variable densidad la cual es medida por el sensor y llevada por
el transmisor al mdulo analgico del PLC. I es el valor de corriente que entrega
el transmisor. EL canal 2 puede ser congurado para leer la variable temperatura
en
o
C o ujo msico el cual es medido segn la conguracin que se tenga para esta
58
3.2 Desarrollo
variable. las unidades para ujo msico pueden ser seleccionadas como se muestran
en la Figura 3.28. En esta ventana se pueden tener hasta 9 posibilidades de congu-
racin para la lectura de unidades, como son kg/h, kg/m, kg/s, g/h, g/m, g/s, t/h,
t/m, t/s.
Figura 3.28.: Unidades para ujo msico
Para convertir el valor que entrega el mdulo analgico a unidades se aplica la
Ecuacin 3.2 y la corriente se calcula la Ecuacin 3.3:
H= Limite alto de la unidad de la variable
L=limite bajo de la variable
A=valor analgico que entrega el mdulo (0-4095)
Unidades = A
H L
4095
+L (3.2)
Corriente =
16A
4095
+ 4 (3.3)
Si algn canal presenta alguna falla se mostrara la ventana de la Figura 3.29. Lo
que puede quiere decir que no existe comunicacin entre el Transmisor y el PLC.
59
Captulo 3 Captulo 3
Figura 3.29.: Error en el Canal de 4-20mA
Fecha y hora
Esta opcin sirve para congurar los datos de fecha y hora del HMI ver Figura 3.30
, puesto que la pantalla HMI no posee un RTC para poder llevar un registro de
la hora y fecha. Esta conguracin se puede realizar de dos maneras, la primera es
seguir los pasos que aparecen en pantalla y la segunda es directamente de la hora
del PLC y simplemente es tocar cada uno de los parmetros y luego dar ENTER.
Para esta segunda opcin el ao si debe ser ingresado por completo es decir, dar
clic en el parmetro del ao y posteriormente ingresar el ao completo, esto debido
a que el PLC lee el ao solamente con los dos ltimos dgitos, es decir si es el ao
2013, el PLC slo muestra el nmero 13, mientras el HMI si toma el ao completo
por lo cual se debe ingresar en esta casilla el ao completo, es decir, 2013.
Figura 3.30.: Conguracin de fecha y hora
60
3.2 Desarrollo
Datos de recibo
En esta opcin se pueden congurar los datos que tendr el recibo o ticket el cual se
imprime con cada descarga. Esta conguracin se puede observar en la Figura 3.31.
Entre las opciones de conguracin se encuentran:
Figura 3.31.: Conguracin de datos del recibo
NIT: es el NIT de la empresa la cual realiza los proceso de scalizacin, en
este caso es el NIT de la empresa Pegasus Blending.
# tiquete inicial: es el nmero de tiquete inicial en donde al congurar este
parmetro, se modica el conteo inicial de los recibos. Si se pone a 0 la cuenta,
el siguiente recibo saldr con el nmero 1.
Opcin de recibo: Existen 2 formatos para el recibo, estos son los que se
muestran en la Figura 3.18 y en la Figura 3.19.
Precio Unidad: Este es el precio que se pone por unidad. Este valor se puede
cambiar en la pantalla PRODUCTO si la opcin 2 del formato de recibo est
seleccionada.
IVA: el impuesto en porcentaje el cual se le aplica a la entrega del producto.
Placa: placa del vehculo al cual se carga con el producto.
Empresa: Empresa ala cual se le ha entregado el producto
Producto: en este caso es GLP, pero se podran utilizar otros productos para
realizar un proceso scal.
3.2.3.4. CALIBRACIN.
En el submenu calibracin se encuentran diferentes opciones para congurar el sis-
tema de una forma optima ver . Dentro de estas conguraciones se encuentran:
61
Captulo 3 Captulo 3
Pulsos/Kilogramo (la conguracin de divisin de pulsos), K-factor (conguracin
de divisin por pulsos), Delta (conguracin de apertura y cierre de actuadores de
control), # de bombas (Conguracin de la activacin de las bombas del sistema).
Figura 3.32.: submenu CALIBRAR
Para poder ingresar a estas conguraciones es necesario ingresar la contrasea de
acceso la cual es usada solo por personal autorizado.
Pulsos/kg
Esta conguracin es usada para la conguracin de la conversin de pulsos a uni-
dades de masa o volumen. Esta ventana se puede observar en la Figura 3.11. En
donde se selecciona de una tabla la cantidad de pulsos que enva el transmisor por
kilogramo.
K-factor
En esta conguracin se pueden congurar la divisin de pulsos ya sea introduciendo
la cantidad de pulsos sobre kilogramo o introduciendo la cantidad deseada de masa
ver Figura 3.12 para que el sistema se autoajuste a la conguracin deseada
2
.
Delta
Esta funcin es activada cuando el sistema se encuentra en modo automtico y se
reere a una compensacin. Cuando se da la orden de cierre de la vlvula y apagado
de las bombas, puede ocurrir segn la instalacin y distancias de estos dispositivos
2
Las conguraciones de la divisin de pulsos puede ser vista en Subseccin 3.2.2 TOTALIZADOR
62
3.2 Desarrollo
al sensor puede pasar mas o menos producto que el deseado, para ello se incorporo
esta funcin en la cual el operario despus de la primer descarga hacia un carro-
tanque se puede dar cuenta de a falta o sobra de producto entregado, entonces el
operario debe ingresar de cuanto se paso o cuanto falto para alcanzar la referencia
introducida en el valor de PRESET, as cuando se realice una nueva descarga, el
sistema se acercara ms al valor introducido en el PRESET. El ingreso de este valor
se puede observar en la Figura 3.33, en donde se observa que aparte de el valor de
delta ingresado, se puede seleccionar si se va a sumar o a restar el valor introducido.
Figura 3.33.: Ventana de Calibracin-Delta
N
o
de bombas
En esta ventana el operario puede congurar como van a trabajar las bombas en-
cargadas de la circulacin del producto. Dado que en la empresa Pegasus Blending
en el proceso de entrega del GLP poseen 2 bombas para la circulacin del produc-
to en donde pueden elegir si trabajan 1 o las 2, esta conguracin se introdujo en
el programa del totalizador en donde por medio de la eleccin que se realice ver
Figura 3.34 se puede prender cualquiera de las 2 bombas o las 2 al tiempo. En la
activacin de una sola bomba se puede seleccionar con cual de las 2 bombas se ira a
trabajar, esto resulta til al momento que se desee realizar un mantenimiento sobre
alguna de la bombas sin necesidad de inhabilitar el funcionamiento del sistema.
63
Captulo 3 Captulo 3
Figura 3.34.: Calibracin-numero de bombas
3.2.3.5. About
En este submenu se muestra informacin general del sistema ver Figura 3.35en donde
se puede observar los elementos generales que componen el sistema: PLC Koyo DL06,
HMI 5043N MAPPLE (con sus respectivas medidas), Impresora EPSON TM-T90P
y por ltimo se muestra la versin del sistema, v1.00, y la empresa a la cual pertenece
el sistema desarrollado: ATLAS E&A.
Figura 3.35.: ABOUT
64
3.2 Desarrollo
3.2.4. MONTAJE DE COMPONENTES EN UN TABLERO
ELCTRICO NEMA 4X
Para el montaje de los elementos de control en un tablero elctrico, se tuvo en cuen-
ta el lugar donde deba ser instalado el gabinete la recomendacin de la empresa
Pegasus Blending de que grado de proteccin debera poseer el gabinete elctrico.
El gabinete elctrico seleccionado posee una proteccin NEMA 4X como se pue-
de observar en la Figura 3.36 la cual muestra la placa caracterstica del gabinete
adquirido.
Figura 3.36.: Placa caracterstica del gabinete
El cableado interno del tablero y montaje se realizo en base a los manuales de
instalacin de los equipos.
3.2.4.1. Consideracin de instalacin para el PLC en gabinete
El PLC koyo DL06 debe ser instalado horizontalmente y no se debe instalar
verticalmente, esto se hace para que el PLC tenga una ventilacin apropiada.
Se debe dejar un espacio mnimo de 39mm entre el PLC y todos los lados del
gabinete.
Debe haber tambin por lo menos 78mm de espacio libre entre la unidad y
cualquier conducto de alambrado que corra en paralelo con las terminales.
El terminal de tierra debe ser conectado a un solo punto de tierra. Se debe
usar alambre trenzado de cobre para conseguir una baja impedancia.
EL PLC no posee un fusible interno por esto se instalo un fusible en el gabinete de
forma que el fusible instalado sea puesto en la fase de alimentacin del PLC y no es
necesario un fusible en el neutro, este fusible debe ser de por lo menos 1A (rpido).
65
Captulo 3 Captulo 3
3.2.4.2. Consideracin para la instalacin del HMI
Para el correcto funcionamiento del HMI segn las recomendaciones del fabricante
se uso una fuente de 24VDC dedicada. y se utilizo cable de calibre 18AWG de 2
conductores con shield. La correcta conexin de la fuente de alimentacin se muestra
en la Figura 3.37. El HMI posee un fusible interno para proteccin de sobrecorriente.
Figura 3.37.: Conexin entre la fuente de 24VDC y el HMI
Figura 3.38.: Distribucin del tablero
Tanto en el HMI como en el PLC se recomienda no tener en la misma canaleta
cables de 110VAC y cables de 24VDC y de activacin o seales de entrada y salida
del PLC, por esto se tuvo en cuenta la distribucin del cableado interno. si llegado
el caso los cables de poder y 24VDC se deben cruzar, estos deben hacerlo de forma
perpendicular y no paralela. Esta distribucin puede ser vista en la Figura 3.38 en
donde se observa la separacin de la linea de alimentacin del sistema con las lneas
66
3.2 Desarrollo
de datos y 24VDC. En la izquierda se coloco el cableado referente a la alimentacin
de 110VAC y a la derecha se coloco la entrada y salida de datos. Segn la norma IEC
60445 que dene el marcado e identicacin delos bornes de equipos, las regletas de
un gabinete se nombran Xn donde n es el nmero de la regleta. Para la identicacin
de borneras se numeran de acuerdo a grupos formados por las regletas, en donde la
numeracin es creciente y se hace de izquierda a derecha.
En la regleta X1 se instalaron elementos de proteccin y alimentacin de los dife-
rentes elementos que componen el sistema, estos elementos son: 1 breaker para la
proteccin de alimentacin principal, 1 breaker para el PLC, un fusible de protec-
cin para el PLC, una caja de distribucin para realizar las conexiones de fase (L),
neutro (N) y tierra (GND), tambin en esta regleta se encuentra la fuente de 24VDC
la cual alimenta el HMI. esto puede ser visto en la Figura 3.39.
En la regleta X2 se instalaron el PLC, y las borneras correspondientes para la
alimentacin del HMI. En la regleta X3 se instalaron las borneras de entrada, salida
y alimentacin de 24 VDC para equipos que requieran este voltaje ver Figura 3.40.
Este voltaje disponible en la regleta X3 es proporcionado por el PLC, y es usado para
la conexin de entradas, salidas y elementos que requieran de una fuente externa
como el transmisor. Las borneras de entradas y salidas se nombraron de la misma
forma que se encuentran en el PLC: X0-X7 son entradas discretas del sistema, X0
es la entrada de pulsos, Y14-Y17 son salidas del sistema en donde se encuentran la
activacin de SSRs y piloto de emergencia. En esta regleta tambin se encuentran
las entradas analgicas de corriente del PLC nombradas como CH1 y CH2, de donde
se obtienen las seales (+) y (-) de cada canal.
Figura 3.39.: Elementos de alimentacin y proteccin
67
Captulo 3 Captulo 3
Figura 3.40.: Entadas y salidas del sistema
El montaje de el sistema en el doble fondo del gabinete se puede observar en la
gura Figura 3.42 donde se pueden ver lo elementos del sistema con sus perspectivas
marquillas. En la Figura 3.41 se pueden observar los rels de estado solido que
activan las bombas del sistema, las caractersticas de estos rels de estado slido
son: voltaje de alimentacin 240VAC, corriente 40A esta corriente se selecciono para
que fuera de 5 a 10 veces la corriente nominal en cargas inductivas [21]. Para el
correcto funcionamiento de los SSRs se debe tener en cuenta factores contra los que
hay que proteger un SSR, las cuales son:
Sobretemperatura
Sobreintensidad
Variaciones bruscas en la Tensin.
Para la proteccin de estos factores se utilizaron: un disipador para evitar efectos
de sobretensin. El disipador usado es de montaje en riel y al cual se le aplico pasta
disipadora para que la transferencia de calor entre el SSR y el disipador sea uniforme.
Los SSRs usados son de referencia HFS15 D-240A25Z-L. El cual posee una red RC
snubber la cual limita las variaciones bruscas de tensin.
Figura 3.41.: SSRs
68
3.2 Desarrollo
Figura 3.42.: Cableado interno y montaje de componentes
Figura 3.43.: Vista frontal del gabinete
69
Captulo 3 Captulo 3
Puerta del gabinete
En la puerta del gabinete se instalaron los botones de control junto a el HMI, tam-
bin se instalaron los pilotos de encendido del sistema y estado de emergencia ver
Figura 3.43. Para la instalacin de los pulsadores e indicadores se perforo el tablero
con crculos de dimetro 22mm. Para el HMI se perforo un rectngulo de dimen-
siones 119x93mm. En el centro de la puerta se encuentra instalado el HMI, encima
de este esta el selector de 2 posiciones el cual hace la funcin de encendido y apa-
gado de el sistema. Al lado de el selector se encuentra un indicador luminoso de
color verde el cual al estar encendido indica que todo el sistema esta alimentado y
encendido. debajo de el HMI se encuentran 3 pulsador los cuales son: Pulsador de
emergencia (pulsador en forma de hongo color rojo), pulsador de rearme (pulsador
azul), pulsador de reset (pulsador negro) e indicador de estado de emergencia (in-
dicador luminoso rojo). Tanto los pulsadores como los pilotos fueron seleccionados
segn lo que establece la norma IEC 60204-1 la cual establece el cdigo de colores
para visualizadores y pilotos.
Piloto rojo: Emergencia-condicin peligrosa que requiere una accin inmediata
Pulsador rojo: Emergencia- accin en caso de peligro.
Los pulsadores, indicadores e interruptor On/OFF instados son de la casa SCHNEI-
DER ELECTRIC los cuales poseen grados de proteccin NEMA 1, 2, 3, 4, 4X, 12, y
13, donde al ser instalados en el gabinete de forma correcta se garantiza el grado de
proteccin NEMA 4X. La parte frontal de la pantalla HMI tambin posee un grado
de proteccin NEMA 4 y 4X y al ser instalada en un gabinete de el mismo grado
de proteccin se garantiza la proteccin que posea dicho gabinete, en este caso La
proteccin general del gabinete es NEMA 4X. Tanto el selector de 2 posiciones y el
indicador de encendido del sistema funcionan a 110VAC, por este motivo se instalo
de una forma alejada de los otros pulsadores los cuales trabajan a 24VDC.
Figura 3.44.: Vista interior de la puerta del gabinete
70
3.2 Desarrollo
Entradas y salidas
El gabinete en su cara inferior posee 8 posibles conexiones para entradas y salidas
tal como lo muestra la Figura 3.45 en estas conexiones a prensacables se encuentran:
Conexin a la red de 110VAC, Entradas de pulsos y 4-20mA del Transmisor, salidas
para la activacin de las bombas que hacen circular el GLP, una conexin para una
vlvula de apertura y cierre y por ultimo se tiene una salida la cual es el cable de
comunicaciones de la impresora la cual fue instalada de forma remota. Se dejaron
conexiones de sobre para posibles conexiones adicionales que se necesiten por si en
un futuro se planea ampliar el sistema, con ms entradas y salidas.
Figura 3.45.: Entradas y salidas del gabinete
3.2.5. INSTALACIN EN CAMPO
Una vez nalizado el montaje de el sistema se desarrollaron pruebas de funciona-
mianto y se entrego el equipo a ATLAS E&A en donde se dio una capacitacin
al personal de la empresa del funcionamianto del equipo y se entrego un manual
requerido por la empresa
3
El personal de ATLAS E&A fu el encargado de la instalacin del dispositivo en la
empresa PEGASUS BLENDING. Para la instalacin de los dispositivos en el lugar
donde se realiza el proceso de descarga se tuvo en cuenta los grados de proteccin de
los elementos, pues si alguno de los elementos que componen el sistema no cumplen
con el suciente grado de proteccin, puede ocurrir un accidente. Tanto el transmisor,
sensor, y elementos de activacin como la bombas son elementos que se encontraban
instalados en el sistema totalizador anterior. El transmisor de ujo RHE 08 esta
aprobado para ser instalado en exteriores para reas de clase 1, divisin 2 y posee
un grado de proteccin IP 65/ Nema 4X. El sensor de ujo coriolis es aceptado para
trabajar en reas de clase I, divisin 1 y 2 y posee una proteccin IP 65 Nema 4X. El
gabinete del sistema totalizador posee un grado de proteccin Nema 4X, mientras
que la impresora de tickets EPSON TM-T90 no cumple con lo grados de proteccin
nombrados anteriormente. Por su manejo no poda ser instalada en el gabinete, por
esto se planteo una solucin para la instalacin del sistema, esta se muestra en la
Figura 3.46. Donde se pens en instalar la impresora de forma remota en una ocina
3
Este manual no se anexo a el documento por permisos de la empresa.
71
Captulo 3 Captulo 3
ubicada a unos 50m de distancia. Como la comunicacin entre los elementos y la
impresora era RS-232, solo se poda instalar esta a una distancia mxima de 15m
como mximo, por esto se decidi utilizar el estndar de comunicaciones RS-485 el
cual permite una longitud mxima de 1200m (a 100kbit/s).
Figura 3.46.: Instalacin de los elementos
En el sitio de instalacin se encontraban las conexiones realizadas del sistema ante-
rior, las cuales fueron tenidas en cuenta para poner en marcha el sistema desarro-
llado. El sistema de alimentacin que proporciona el voltaje tanto AC como DC ver
Figura 3.47 fue utilizado para alimentar el gabinete, este necesitaba una alimenta-
cin de 110 VAC, el sistema de alimentacin tambin proporciona una alimnetacin
de 24 VDC la cual era utilizada para alimentar el transmisor del totalizador anterior.
Figura 3.47.: Alimentacin disponible.
El sensor RHM 20 fue desmontado de la lnea de en la cual circula el GLP ver jjj
para que la empresa ATLAS E&A realizara la respectiva calibracin se este, con el
72
3.2 Desarrollo
fn de tener un a medida adeacuada. La conexin es de 11/2 ANSI 150, adems
cuando se intale de nuevo el sensor se deben cambiar los empaques y los tornillos.
Figura 3.48.: Punto donde se desmonto el sensor
3.2.5.1. Impresora
La impresora EPSON TM-T90 es una impresora para tickets la cual es soportada
por el HMI 5043N de Mapple. Esta impresora aunque posee una funcin de auto-
cortado de papel, El HMI no enva ese comando para realizar esta funcin. Para
realizar la funcin de auto-cortado de papel despus de una impresin se utilizo
el PLC para que enviara el comando necesario. La impresora utiliza los comandos
ESC/POS[22] para realizar diferentes acciones, entre ellas se encuentra la funcin
de cortado de papel. El comando correspondiente a la funcin de cortado de papel
es:
ASCII: GS V m
Hex: 1D 56 m
Si m=0, la impresora corta el papel totalmente, mientras si m=1 el papel sera
cortado parcialmente. Tanto el PLC como el HMI pueden enviar datos a travs del
protocolo RS-485, pero la impresora no puede recibir este tipo de datos ya que solo
puede recibir datos mediante RS-232. Para esto se uso un conversor RS-485/RS-232
de la marca Hexin technology. Los voltajes de alimentacin del conversor pueden ser
entregados por los pines RTS y DTR de el conector DB25[23] de la impresora as
que no es necesario una forma adicional de alimentacin. El tipo de conexin RS-485
la cual es de 2 hilos. El diagrama de conexin de el PLC, el HMI y la impresora
puede ser visto en Apndice D.
73
4. Captulo 4
En este captulo Se hace referencia a el estudio para la implementacin de algorit-
mos de compensacin y correccin de volumen para GLP (tabla API 54E ) y gas
AGA reporte 8. Estos algoritmos fueron desarrollados en Matlab en donde se estu-
dio los requerimientos necesarios para la implementacin en sistemas que necesiten
compensaciones para GLP y gas los cuales pueden en un futuro ser implementados
en campo. El algoritmo AGA 8 fue desarrollado para el cumplimiento de uno de los
objetivos que se haban propuesto en el anteproyecto realizado. En donde se tuvo
una confusin por parte de la empresa en cuanto a el producto que se hiba a to-
talizar. Al principio se tenia pensado para gas, pero despus se cambio a GLP. En
el captulo se desarrolla el algoritmo para correccin de volumen para GLP el cual
no fue implementado por requerimientos de la empresa PEGASUS BLENDING en
donde se totaliza por masa y si ellos desean conocer el valor del volumen, utilizan
un cromatgrafo para saber el valor de la densidad y despus hacen la conversin.
4.1. CORRECCIN DE VOLUMEN PARA GLP
En esta ltima seccin del capitulo se presenta un mtodo para la correccin de
medicin de volumen el cual puede ser implementado en una posible ampliacin
del sistema Totalizador, para poder realizar un proceso de transferencia de custodia
adecuado. Tambin se presenta lo necesario para la aplicacin de este mtodo en el
sistema actual instalado en la empresa Pegasus Blending.
Para propsitos de transferencia de custodia, gas natural licuado (GNL) y gas li-
cuado de petrleo (GLP), la medicin de volumen es corregida a una temperatura
base y presin de saturacin. En la mayora de transferencias de volumen, la tem-
peratura y presiones cambian de las condiciones estndar. Estos volmenes deben
ser ajustados a condiciones estndar a travs de factores de correccin[7].
4.1.1. API M.P.M.S 11.2.4
En el libro MPMS, captulo 11, seccin 2, parte 4, establece mtodos para calcular los
factores de correccin para GLP, en el cual aparecen procedimientos de implemen-
tacin, ejemplos de tablas, diagramas de ujo, y ejemplos especcos. Este estndar
cubre rangos para densidades relativas a 60
o
F desde 0.3500 hasta 0.6880 kg/L los
75
Captulo 4 Captulo 4
cuales son normalmente iguales a una densidad a 15
o
C de 351.5 hasta 687.8 kg/m
3
.
Los rangos de temperatura de este estndar estn en el rango de -50.8 hasta 199.4
o
F
(-46 hasta 93
o
C).
El procedimiento de implementacin describe:
1. EL clculo de C
TL
(factor de correccin de temperatura) dando un factor de
densidad apropiado a temperatura base y a una temperatura observada.
2. El clculo apropiado del factor de densidad a una temperatura base dando una
densidad relativa y una temperatura observada.
Las tablas 24E, 54E y 60E calculan el C
TL
,las tablas 23E, 53E y 59E calculan la
densidad relativa a la temperatura base.
4.1.2. Eleccin de la Tabla apropiada
Para saber con cual tabla se va a trabajar, hay que conocer primero que variables
entrega el sistema y lo que se desea calcular. En nuestro caso necesitamos calcular el
factor de correccin C
TL
para corregir el volumen que pueda ser ledo por el sensor
coriolis u otro sensor el cual pueda entregar una lectura del volumen actual. Para el
calculo del C
TL
se disponen de las tablas 24E, 54E y 60E. En las 3 tablas se ingresan
los valores de temperatura y densidad. La diferencia es el dato de la densidad base
utilizada. Para el calculo de C
TL
la tabla 24E necesita los valores de la densidad
relativa a 60
o
F y la temperatura observada en
o
F. La tabla 54E necesita los datos
de la densidad a 15
o
C y la temperatura observada en
o
C. La tabla 60E necesita los
datos de densidad a 20
o
C y la temperatura observada en
o
C. El conjunto sensor
y transmisor de temperatura entregan el valor de la temperatura del producto el
cual esta siendo medido. El sistema totalizador instalado muestra el dato de esta
temperatura en
o
C. En cuanto a la densidad, la empresa en donde se realice el
proceso de transferencia de custodia debe tener el dato de la densidad ya sea a 15
o 20
o
C. En pegasus Blending poseen el dato de la densidad a 15
o
C, lo que quiere
decir que se debe trabajar con la tabla 54E.
4.1.3. TABLA 54E
En [7] se describe el procedimiento computacional que se debe realizar para calcular
el factor C
TL
el cual es usado para calcular volmenes de uidos a una temperatura
base desde volmenes a una temperatura conocida. La tabla 54E es usada para
uidos caracterizados con la especicacin de densidad a temperatura base, 15
o
C.
Las entradas y salidas de el algoritmo a desarrollar son:
Entradas: Densidad a 15
o
C,
15
(kg/m
3
)
Temperatura Observada, T
F
(
o
C)
Salida: Factor de correccin de temperatura, C
TL
(desde T
F
hasta T
B
)
76
4.1 CORRECCIN DE VOLUMEN PARA GLP
Para el clculo del factor se usan 2 uidos de referencia a 60
o
F los cuales son selec-
cionados para que uno de ellos posea un valor de densidad un poco mas alto que el
de referencia y el otro que sea un poco mas bajo. Las densidades de los uidos de
referencia tienen que ser escalados a la temperatura observada reducida (reducida
por la temperatura crtica de el uido de inters). en la Tabla 4.1 se observan los
uidos de referencia.
Cuadro 4.1.: Parmetros de los uidos de referencia[7]
60
: Densidad relativa a 60
o
F
T
c
: Temperatura crtica en Kelvin
Z
c
: Factor de compresibilidad crtico
c
: Densidad crtica en gramos-moles por litro
k
1
, k
2
, k
3
, k
4
: Parmetros de ajuste de saturacin de densidad
EE(68/32) Denota una mezcla de 68 mole % ethane + 32 % ethylene
EP(65/35) Denota una mezcla de 65 mole % ethane + 35 % propane
EP(35/65) Denota una mezcla de 35 mole % ethane + 65 % propane
El procedimiento de implementacin esta descritos en 11 pasos en donde realizan
diversas operaciones para llevar la densidad a una densidad relativa a 60
o
F y despus
hacer clculos de densidades y factores de correccin para los uidos de referencia.
Los clculos son realizados con decimales de aproximadamente 16 dgitos . En [7]
se presenta tambin un diagrama de ujo de el proceso para el clculo del factor de
correccin ver Apndice G.
4.1.3.1. PROCEDIMIENTO DE IMPLEMENTACIN
El procedimiento descrito en [7] para la implementacin de la tabla 54E consta de
11 pasos, los cuales necesita primero los valores de temperatura actual T
F
(
o
C) y el
valor de la densidad base
15
(kg/m
3
).
77
Captulo 4 Captulo 4
Paso 1
Redondear el valor de densidad a el valor mas cercano 0.1 y la temperatura observada
a 0.05
o
C.
Paso 2
Convertir el valor redondeado de temperatura a unidades Kelvin, T
x
.
T
x
= T
F
+ 273.15 (4.1)
Paso 3
La temperatura en Kelvin, T
x
y la densidad deben estar dentro de un rango apropia-
do
1
.
La temperatura debe estar dentro de 227.15 y 336.15K (-46 y 93
o
C).
Densidad entre 351.7 y 687.8kg/m
3
Paso 4
Convertir la densidad a 15
o
C a una densidad relativa, relativa a la densidad de el
agua a 60
o
F.
TB
=

15
999.016
(4.2)
Paso 5
Usar el procedimiento descrito en la tabla 23 para calcular una densidad relativa a
60
o
F desde la densidad relativa conocida a 15
o
C. Las entradas de implementacin
de la tabla 23 son los valores de T
BK
y
TB
, donde T
BK
es la temperatura base
15
o
C en Kelvin (288.15
o
K ) y
TB
es la densidad a la temperatura base 15
o
C. El
procedimiento a implementar de la tabla 23 es desde el paso 4 hasta el paso 11.
1
Si los valores de temperatura y densidad no estn dentro del rango, entonces el estndar no
aplica
78
Paso 5.1
Los uidos de referencia se escogen para el desarrollo del clculo de densidades. Los
uidos de referencia tienen que ser calculados ala temperatura observada, T
x
. Se
debe usar los parmetros de los uidos de referencia descritos en Tabla 4.1. Se usa
la temperatura crtica T
c,ref
de cada uido para calcular su temperatura observada
reducida, T
r,x
:
T
r,x
=
T
x
T
c,ref
(4.3)
Si T
r,x
1, se debe calcular la densidad de saturacin para cada uido de referencia a
su temperatura reducida T
r,x
. Este procedimiento est descrito en el paso nmero 10
de la Tabla 24. Se debe repetir este paso para 60
o
F usando la temperatura reducida
T
r,60
:
T
r,60
=
519.67
1.8T
c,ref
(4.4)
Si por el contrario T
r,x
> 1 EL uido de referencia no sera un lquido a su tempera-
tura observada y el valor
x,ref
no puede ser calculado.
Paso 5.1.1 Clculo de densidad de saturacin.
Se calcula la densidad de saturacin para ambos uidos uidos de referencia a 60
o
F
usando la temperatura reducida, T
r,60
.Para cada uido las ecuaciones que calculan
la densidad de saturacin de cualquier temperatura reducida son:
= 1 T
r
(4.5)
sat
=
c
_
1 +
(k
1

0.35
) + (k
3

2
) + (k
4

3
)
1 + (k
2

0.65
)
_
(4.6)
Donde k
1
, k
2
, k
3
y k
4
son parmetros para cada uido de referencia.
Finalmente se calcula la densidad relativa a la temperatura observada,
x,ref
=
60,ref
_

sat
x,ref
sat
60,ref
_
Donde
60,ref
es la densidad de relativa a 60
o
F de cada uido de referencia.
79
Captulo 4 Captulo 4
Paso 5.2
Determinar 2 uidos adyacentes de referencia para se usados en los clculos. Se debe
elegir el uido de referencia con la densidad mas baja que el valor mas grande que
x
y referir este uido con el nmero 2. Tambin se usa el el siguiente uido de
referencia con la densidad mas baja y se reere esta al nmero 1.
Paso 5.3
Se inician los lmites de la iteracin para la densidad relativa de los uidos obser-
vados. Para la mayora de los casos, la densidad relativa de los uidos observados
debera estar entre los 2 uidos de referencia 1 y 2,
60,1
y
60,2
.
Se inicia el lmite superior para la densidad relativa del uido observado,
60,high
,
como:
60,high
=
60,2
(4.7)
La correspondiente densidad relativa a la temperatura observada,
x,high
, como:
x,high
=
x,2
(4.8)
El limite inferior como:
60,low
=
60,1
(4.9)
La correspondiente densidad relativa a la temperatura observada,
x,high
, como:
x,low
=
x,1
(4.10)
Paso 5.3
Se calcula un valor de densidad relativa 60
o
F intermedia,
60,mid
.
=

x
x,low
x,high
x,low
(4.11)
80
si el valor de es menor que 0.001 entonces poner este valor igual a 0.001; si el valor
es mas grande que 0.999, entonces poner este valor igual a 0.999.
60,mid
=

60,high
+
60,low
2
(4.12)
Calcular el factor de correccin C
TL
, usando el valor de
60,mid
y T
x
, utilizando el
procedimiento descrito en los pasos 5 al 13 de la tabla 24. La densidad relativa
x,mid
a la temperatura observada T
x
, como:
x,mid
= CTL
60,mid
(4.13)
Paso 5.3.1 Clculo del CTL
Usando la tabla de el Tabla 4.1, las densidades relativas a 60
o
F, calcular la variable
de interpolacin, :
=

60
60,1
60,2
60,1
(4.14)
De la tabla, con los valores de temperaturas crticas, calcular la temperatura crtica
del uido, T
c
:
T
c
= T
c,1
+(T
c,2
T
c,1
) (4.15)
Despus se calcula la temperatura observada reducida, T
r,x
:
T
r,x
=
T
x
T
c
(4.16)
La temperatura reducida a 60
o
F, T
r,60
:
T
r,60
=
519.67
1.8T
c
(4.17)
De la tabla de referencia y con los factores de compresibilidad Z
c
y las densidades
crticas
c
, se calcula el factor escalar h
2
.
h
2
=
Z
c,1

c,1
Z
c,2

c,2
(4.18)
81
Captulo 4 Captulo 4
Luego se calcula la densidad de saturacin para ambos uidos utilizando la Ecuacin 4.5
y Ecuacin 4.6.
Calcular el factor de interpolacin X:
X =

sat
60,1
1 +
__

sat
60,1
h
2
sat
60,2
_
1
_
(4.19)
Por ltimo se calcula el factor de correccin.
C
TL
=

sat
x,1
X
_
1 +
_

sat
x,1
h
2
sat
x,2
1
__
(4.20)
Paso 5.4
Se chequea si el valor converge, para que exista convergencia debe existir
x
debe estar
x,low
y
x,mid
y la diferencia entre
60,low
y
60,mid
es menor que
0.00000001 (10
8
)
x
debe estar
x,high
y
x,mid
y la diferencia entre
60,high
y
60,mid
es menor
que 0.00000001 (10
8
)
Si el valor converge, entonces,
60
=
60,mid
Paso 5.5
Existen 3 pares de puntos de los valores de densidades relativas: (
x,low
,
60,low
) , (
x,mid
,
60,mid
) y (
x,high
,
60,high
).
Una ecuacin cuadrtica puede ser arreglada a travs de esos 3 puntos. Esta ecua-
cin cuadrtica debera ser una buena aproximacin a la relacin actual entre
x
y
60
. Usando el valor de la densidad relativa observada en la ecuacin cuadrtica,
debera dar una muy buena estimacin a
60
.
Primero hay que calcular los parmetros para la ecuacin cuadrtica:
= (
60,high
60,low
)
=
_
2
x,high
2
x,low
_
82
=

x,high
x,low
x,mid
x,low
A =
(
60,mid
60,low
)

_
2
x,mid
2
x,low
_
B =
AB
x,high
x,low
El valor asociado a
60,trial
se calcula usando estos valores de A B y C.
60,trial
= A
2
x
+B
x
+C (4.21)
Finalmente se calcula el factor de correccin de temperatura, usando el valor
60,trial
para calcular la densidad relativa a la temperatura observada
x,trial
:
x,trial
= C
TL

60,trial
(4.22)
Paso 5.6
Se revisa la convergencia de la densidad relativa. El calculo converge cuando el valor
absoluto de la diferencia entre
x,trial
y
x
es menor que 0.00000001 (10
8
). Si el
valor converge, entonces:
60
=
60,trial
Paso 6
El resultado
60
debe estar en el rango de 0.3500 y 0.6880. Si es as se puede continuas
y si por el contrario est por debajo de este rango, entonce el estndar no aplica.
83
Captulo 4 Captulo 4
Paso 7
Se debe usar el procedimiento dela Tabla 24 para calcular el factor de correccin
de temperatura C
TL1
de 60
o
F a la temperatura observada T
x
, Este paso provee un
factor usado para reducir un volumen observado a T
x
a un volumen a 60
o
F, cuando
el valor de densidad relativa 60
o
F,
60
es conocida.
Por denicin:
C
TL1
=
V
60
V
Tx
=

Tx
60
(4.23)
Paso 8
Usar el procedimiento de la Tabla 24 para calcular el factor de correccin de tempe-
ratura C
TL2
de 60
o
F a la nueva temperatura base 15
o
C, Este paso provee un factor
usado para reducir un volumen observado a 15
o
C a un volumen a 60
o
F, cuando el
valor de densidad relativa 60
o
F,
60
es conocida.
Por denicin:
C
TL2
=
V
60
V
15
=

TB
60
(4.24)
Paso 9
En este paso se calcula el factor C
TL
deseado para reducir un volumen desde una
temperatura observada, T
F
, a la condicin base de 15
o
C. Las frmulas denidas
muestran que el calculo es hecho por la razn C
TL1
/C
TL2
.
C
TL1
C
TL2
=
_
V
60
V
Tx
_
_
V
60
V
15
_
=
V
15
V
Tx
C
TL
=
V
15
V
Tx
=

Tx
15
(4.25)
Paso 10
Comprobar que el valor de C
TL
sea positivo, si es menor o igual a 0 se debe generar
un error.
84
Paso 11
Por ltimo se redondea el factor de correccin a el valor mas cercano a 0.00001.
4.1.4. Ejemplo de Aplicacin
Aplicando el algoritmo desarrollado en MATLAB ver Apndice G se desarrollo el
ejemplo 54/3 el cual utiliza EP(65/35) y EP (35/65) como uidos de referencia. Este
ejemplo tiene como entradas una densidad base de 451.09kg/m
3
y una temperatura
observada de 30.774
o
C.
La simulacin del algoritmo genera como resultado:
CTL =
0.930274971426775
CTLrounded =
0.930270000000000
Al comparar con el ejemplo mostrado en la Figura 4.1 se observa que se obtuvo el
mismo resultado, este resultado es C
TL
= 0.93027.
Al aplicar el factor de correccin a la lectura de volumen observada se deber aplicar
la ecuacin.
V
15
= C
TL
V
Tx
4.1.5. Criterios para una implementacin del algoritmo en PLC
Para la implementacin del algoritmo en el PLC, primero se desarrollo un algoritmo
en MATLAB y se comparo los resultados con uno de los ejemplos
2
de la tabla 54E.
Cuando se obtuvo el resultado deseado, el algoritmo desarrollado resulto ser mas
largo de lo esperado. Despus se vio la posibilidad de realizar el mismo algoritmo
en el PLC del sistema totalizador
3
, pero existen limitaciones para las funciones que
este tipo de PLC ofrece, este tipo de PLC solo puede ser programado bajo lenguaje
ladder
4
y no soporta el calculo de funciones como exponenciales, potencias y lo-
garitmos. Aunque en la implementacin del algoritmo se realizan clculos usando
las funciones bsicas como multiplicaciones, divisiones, sumas y restas y procesos
repetitivos los cuales son fciles de aplicar, pero en una parte del algoritmo para el
2
En el apartado del libro API M.P.M.S 11.2.4 aparecen varios ejemplos para cada tabla de im-
plementacin
3
PLC koyo DL06, remitirse a Subseccin 3.2.2
4
Mas conocido como programacin en escalera, es uno de los 5 tipos de lenguajes mencionados
en el IEC61131-3 para la programacin de PLCs
85
Captulo 4 Captulo 4
Figura 4.1.: Ejemplo de la tabla 54E. [7]
clculo de la densidad de saturacin ver Ecuacin 4.6 para los uidos de referencia
usando la temperatura reducida, hay que realizar clculos de potencias cuyo expo-
nente es un numero racional, en este caso
0.35
y
0.65
. Los cuales resultan imposibles
de calcular para el PLC y para poder calcular estos valores, habra que valerse de
mtodos como el CORDIC
5
, lo que conllevara a un algoritmo Ladder ms extenso
y que ocupara una cantidad de memoria la cual es limitada.
sat
=
c
_
1 +
(k
1

0.35
) + (k
3

2
) + (k
4

3
)
1 + (k
2

0.65
)
_
Para la implementacin de este algoritmo sera aconsejable el uso de PLCs los cuales
puedan ser programados con lenguaje estructurado en donde se pueda realizar el
calculo de potencias racionales, y teniendo el algoritmo desarrollado en lenguajes
estructurado, se puede utilizar cualquier PLC que soporte este tipo de lenguaje.
5
Cordinate Rotation Digital Computer. Algoritmo para calcular funciones trigonomtricas me-
diante rotacin de vectores.
86
4.2 AGA 8
4.2. AGA 8
En esta seccin se hace referencia al factor de compresibilidad Z para con este realizar
la correccin de la medida del instrumento dependiendo de la temperatura y presin
presentes en el proceso. Para esto nos basamos en la norma AGA reporte nmero 8
el cual trata a fondo el clculo de este factor de compresibilidad y brinda las bases
para el desarrollo de un algoritmo que pueda ser implementado en un instrumento
totalizador.
De acuerdo a lo anterior esta seccin dene lo que es el factor de compresibilidad,
muestra las bases y el paso a paso de lo necesario para realizar su clculo y se
desarrolla un algoritmo genrico en MATLAB el cual pueda ser fcilmente llevado
a cualquier tipo de PLC o instrumento totalizador en el que se quiera implementar
este tipo de correccin, de igual manera los resultados obtenidos con este algoritmo
diseado son comparados con los resultados experimentales y de referencia dados en
la norma AGA reporte 8.
4.2.1. Factor de compresibilidad Z
El Factor de compresibilidad (Z) se dene como la razn entre el volumen molar de
un gas real (Vreal) y el correspondiente volumen de un gas ideal (Videal).
Se utiliza para comparar el comportamiento de un gas real respecto al establecido
por la ecuacin de los Gases Ideales. Partiendo de esta denicin Z representa un
factor de correccin para la ecuacin de los gases ideales.
Con base en esto se encuentra tres tipos de comportamiento distintos:
Z = 1, comportamiento de Gas Ideal. (Altas temperaturas y bajas presiones).
Z > 1, gases como el Hidrgeno y Nen, difcilmente compresibles (altas tem-
peraturas y presiones).
Z < 1, gases como el O2, Argn y CH4, fcilmente compresibles (bajas tem-
peraturas y altas presiones).
4.2.2. MTODO DE LA CARACTERIZACIN DEL DETALLE
El mtodo de la caracterizacin del detalle fue desarrollado para describir exacta-
mente el comportamiento de la presin - densidad de la fase de gas de las mezclas
del gas natural sobre una amplia gama de condiciones. Puede tambin describir
exactamente el comportamiento de la densidad con base en la presin - temperatura
de la fase de gas de los componentes de metano puro, etano, bixido de carbono,
nitrgeno e hidrgeno y mezclas de estos componentes. Adems, una correlacin
de la baja densidad fue desarrollada para el propano y mezclas ms pesadas del
hidrocarburo y binarias de estos componentes con bixido de metano, de etano, de
87
Captulo 4 Captulo 4
nitrgeno y de carbono. Este mtodo reduce la incertidumbre de los clculos del fac-
tor de compresibilidad y de la densidad para los gases naturales de los separadores
de la produccin, que pueden contener porcentajes de moles de hexanos ms hidro-
carburos ms pesados el mayores que el 1 %. Las correlaciones del comportamiento
de la densidad del sulfuro puro del hidrgeno y de las mezclas binarias del sulfuro
de hidrgeno con bixido de metano, de etano, de nitrgeno y de carbono fueron
desarrolladas para reducir la incertidumbre del clculo para los gases naturales que
contenan el sulfuro de hidrgeno (gas amargo).[24]
Finalmente, las segundas correlaciones viriales fueron desarrolladas para el agua
y las mezclas binarias del agua con bixido de metano, de etano, de nitrgeno y
de carbono para reducir la incertidumbre del clculo para los gases naturales que
contenan el vapor de agua (gas mojado).
4.2.3. ECUACIONES GENERALES
Esta seccin contiene las ecuaciones generales que implican el factor de compresibi-
lidad, la densidad molar y la densidad de masa de los gases del gas natural y otros
hidrocarburos.
El factor Z de la compresibilidad es denido por la ecuacin:
Z =
PV
RT
(4.26)
El mtodo de caracterizacin de detalle expresa el factor Z de la compresibilidad en
trminos de la densidad molar d:
d =
n
V
(4.27)
La masa molar M
r
de la mezcla del gas, se calcula de la composicin usando la
relacin:
M
r
=
N
i=1
x
i
M
ri
(4.28)
La suma en la Ecuacin 4.28 es sobre todos los componentes en la mezcla del gas.
La densidad total p es relacionada con la densidad molar d por la relacin:
p = M
r
d (4.29)
88
4.2 AGA 8
Entonces, usando la Ecuacin 4.27 y Ecuacin 4.29, en la Ecuacin 4.26, las ecua-
ciones siguientes para densidad molar d y la densidad de masa p en trminos del
resultado del factor de compresibilidad:
d =
P
ZRT
(4.30)
p =
M
r
P
ZRT
(4.31)
4.2.3.1. Densidad total en las condiciones de referencia
La densidad total p
b
en la condicin de referencia (condicin base) T
b
P
b
puede ser
calculada con la siguiente relacin:
p
b
=
M
r
P
b
Z
b
RT
b
(4.32)
4.2.3.2. Condicin de referencia del factor de Supercompresibilidad
Las tabulaciones del factor F de supercompresibilidad presentados por el A.G.A.
son denidos por la relacin siguiente:
F
2
pv
=
Z(60F.14, 73psia)
Z(T.P)
(4.33)
En esta ecuacin, la condicin 60F, 14.73 psia es una condicin de referencia espe-
cca. Observe que la condicin de referencia iguala la condicin de referencia para
el factor del supercompresibilidad solamente cuando la condicin de referencia es
60F, 14.73 psia.
4.2.4. ECUACIONES PARA EL FACTOR DE
COMPRESIBILIDAD
Las siguientes ecuaciones son para hallar el factor de compresibilidad Z por medio
del mtodo de caracterizacin detallado
89
Captulo 4 Captulo 4
4.2.4.1. Ecuacin de estado del mtodo de la caracterizacin de detalle
para el factor de compresibilidad
Las ecuaciones, las constantes, y los parmetros necesitaron calcular los factores de
la compresibilidad para las mezclas del gas natural usando el DETALLE que EL
MTODO de la CARACTERIZACIN se da aqu.
La ecuacin de estado usada aqu es una formulacin hbrida. El mtodo se basa
en el trabajo de Starling y otros. Combina caractersticas de la ecuacin virial de
estado (una serie de energa en densidad) para las condiciones de baja densidad y las
funciones exponenciales para los usos en las condiciones de alta densidad (ecuacin
extendida de Benedict-Webb-Rubin). [24]
Esta formulacin proporciona alta exactitud, la amplia gama del uso de la temperatura-
presin-composicin, y capacidades termodinmicas derivadas de la caracterstica.
La ecuacin de estado para el factor z de compresibilidad para el mtodo de la
caracterizacin de detalle es dada en su forma condensada por la ecuacin siguiente:
Z = 1+
DB
K
3
D
18
n=13
C
n
T
u
n
+
58
n=13
C
n
T
u
n
(b
n
c
n
k
n
D
k
n
)D
k
n
exp(c
n
D
k
n
) (4.34)
La densidad D de la reduccin es relacionada con la densidad molar d por la ecuacin:
D = K
3
d (4.35)
Se calcula el parmetro K del tamao de la mezcla usando la ecuacin:
K
5
=
_
N
i=1
x
i
k
5/2
i
_
2
+ 2
N1
i=1
N
j=i+1
x
i
x
j
(K
5
ij
1)(K
i
K
j
)
5/2
(4.36)
El subscrip i reere al i esimo componente en la mezcla del gas y el subndice j
reere al j esimo componente en la mezcla. En la suma simple, i se encuentra en los
valore enteros a partir del 1 a N. por ejemplo, para una mezcla de 12 componentes,
de N = 12 y all seran 12 trminos en la suma simple. En la suma doble, i se
extiende a partir de 1 a N1 y, para cada valor de i, j se extiende de i +1 al N. Por
ejemplo, para la mezcla de 12 componentes, habra 66 trminos en la suma doble si
todos los valores de K
ij
dieren de uno a excepcin de los valores en la tabla 6.
El segundo coeciente virial B es dado por las ecuaciones siguientes, donde est el
nmero N de componentes en la mezcla del gas y de los valores de i y de j ambos
a partir de 1 a N:
B =
18
n=1
a
n
T
U
n
N
i=1
N
j=1
x
i
x
j
E
U
n
ij
(K
i
K
j
)
3/2
B
*
nij
(4.37)
90
4.2 AGA 8
B
*
nij
= (G
ij
+1g
n
)
g
n
(Q
i
Q
j
+1q
n
)
q
n
(F
1/2
i
F
1/2
j
+1f
n
)
f
n
(S
i
S
j
+1s
n
)
s
n
(W
i
W
j
+1w
n
)
w
n
(4.38)
Debe ser observado que W
i
es cero para todos los componentes excepto el agua y
esa F
j
es cero para todos los componentes excepto el hidrgeno.
Los parmetros binarios E
ij
y G
ij
se calculan usando las ecuaciones siguientes:
E
ij
= E
ij
(E
i
E
j
)
1/2
(4.39)
G
ij
=
G
ij
(G
i
+G
j
)
2
(4.40)
Se observa que todos los valores de los parmetros binarios de la iteracin son uno
a excepcin de los valores en la tabla 6.
Los coecientes (n=13 a 58) estn dados por la ecuacin:
C
n
= a
n
(G + 1 g
n
)
g
n
(Q
2
+ 1 q
n
)
q
n
(F + 1 f
n
)
f
n
U
u
n
(4.41)
Se calculan los parmetros U, G, Q, y F de la mezcla usando las ecuaciones siguien-
tes, donde en las sumas dobles, i se extiende a partir de 1 a N 1 y, para cada valor
de i, j se extiende dei + 1 a N:
U
5
=
_
N
i=1
x
i
E
5/2
i
_
2
+ 2
N1
i=1
N
j=i+1
x
i
x
j
(U
5
ij
1)(E
i
E
j
)
5/2
(4.42)
Debe ser observado que todos los valores de los parmetros binarios de la interaccin
K
ij
, E
ij
, G
ij
, y U
ij
son uno a excepcin de los valores mostrados en la tabla 6.
Tambin observe que F
i
es cero para todos los componentes excepto el hidrgeno.
En Apndice F se encuentran las tablas 4, 5 y 6 del AGA reporte 8
Con lo anterior se tiene lo necesario para realizar el clculo del factor de compresi-
bilidad Z, para esto desarrollamos un algoritmo general el cual pueda ser fcilmente
trasladado a otros lenguajes y ser implementado en casi cualquier tipo de PLC o
dispositivo involucrado en la medicin de GAS.
91
Captulo 4 Captulo 4
4.2.5. ALGORITMO GENERAL AGA REPORTE 8 EN
MATLAB
El algoritmo general es desarrollado en MATLAB el cual trabaja con un lenguaje
sencillo y fcilmente entendible, adems este cuenta con las mejores herramien-
tas para trabajar cualquier clculo matemtico. Este algoritmo puede ser visto en
Apndice F.
El algoritmo desarrollado en MATLAB utiliza todas las ecuaciones mencionadas en
el AGA reporte 8 para la solucin de la ecuacin de estado del mtodo de caracteriza-
cin detallado, ahora bien para comprobar que este algoritmo ha sido correctamente
implementado, el AGA reporte 8 brinda dos tablas, la tabla A.5-1 y la tabla A.5-2.
La tabla A.5-1 muestra 5 tipos de mezclas de gases y sus porcentajes de compo-
nentes de mol, es decir X
i
, y que estos datos sirvan como entrada para el algoritmo
desarrollado. Esta tabla se encuentra en Apndice F
Por su parte la tabla A.5-2 muestra los valores del factor de compresibilidad obte-
nidos de forma experimental para cada una de la mezclas de gases de la anterior
tabla.
Estos valores de Z se dan para 4 valores diferentes de temperatura y 8 valores di-
ferentes de presin, valores que son usados tambin como entradas en el algoritmo
desarrollado para comparar los valores obtenidos por medio del mtodo de caracte-
rizacin detallada.
El Tabla F.6 y Tabla F.7 para cada tipo de mezcla de gas fueron preparadas por el
AGA para que sirvan como referencia y con propsitos de vericacin de resultados
de los algoritmos desarrollados.
Puesto que los resultados obtenidos experimentalmente (tabla A.5-2) no siempre
son exactos comparados con los valores obtenidos por medio del mtodo de carac-
terizacin detallado, es necesario generar estadsticas entre los valores calculados y
los experimentales para obtener de esta manera una incertidumbre de los valores
calculados.
Para ello el AGA tambin provee las siguientes frmulas para generar las estadsticas
necesarias de vericacin.
Z
diff
=
Z
calc
Z
exp
Z
exp
100 % (4.43)
BIAS =
1
N
N
i=1
Z
diff,i
(4.44)
92
4.2 AGA 8
AAD =
1
N
N
i=1
_
(Z
diff,i
)
2
_
1/2
(4.45)
Std.Dev. =
_
1
N 1
N
i=1
(Z
diff,i
BIAS)
2
_
1/2
(4.46)
RMS =
_
1
N
N
i=1
(Z
diff,i
)
2
_
(4.47)
Max.Dev = (Z
diff
)
max
(4.48)
Z
diff
Porcentaje relativo de diferencia entre los factores de compresibilidad calcula-
dos y experimentales
Z
diff,i
Zdi para el i-simo dato
(Z
diff
)
max
Valor mximo de
Z
calc
Factor de compresibilidad calculado
Z
exp
Factor de compresibilidad experimental
AAD Promedio de desviacin absoluta. BIAS =bias
Std.Dev Desviacin estndar
RMS Desviacin de la raz media cuadrtica
Max.Dev. Valor mximo de Zdi
93
5. Captulo 5
Luego del desarrollo de la metodologa, se obtuvo una serie de resultados para ca-
da uno de los objetivos propuestos, entre otros resultados. Dichos resultados son
presentados y analizados a continuacin.
5.1. Resultados
5.1.1. Comparacin entre LectroCount LCR y el sistema
Totalizador desarrollado.
En la construccin y elaboracin del totalizador se bas en otros equipos y sistemas
totalizadores para adaptarle las funciones necesarias o mnimas. En el anterior proce-
so de descarga se manejaba un totalizador con gran variedad de funciones por lo cual
se hizo nfasis en este totalizador para realizar la comparacin con el desarrollado y
as poder comprobar los resultados obtenidos con el nuevo sistema implementado.
El totalizador usado anteriormente era un LectroCount LCR II de LIQUID CON-
TROLS GROUP que se observa en Figura 5.1
Figura 5.1.: LectroCount LCR 600. [8]
95
Captulo 5 Captulo 5
Este totalizador recibe las seales del transmisor para poder realizar diversos clculos
y poder mostrar las medidas en pantalla por lo cual funciona de manera muy similar
al totalizador desarrollado para este proyecto el cual se muestra en Figura 5.2
Figura 5.2.: Totalizador ATLAS E&A
Comparando ambos totalizadores se obtuvo como resultado la tabla mostrada en
Apndice E.
De la anterior comparacin se puede decir que el equipo desarrollado cumple con
las funciones bsicas que debe tener un totalizado y las funciones adicionales se
realizaron de acuerdo a los requerimientos hechos por la empresa para el proceso de
descarga.
Debido a que el totalizador desarrollado era ms personalizado de acuerdo a los
requerimientos hechos, no se implementaron ciertas funciones que pueden tener el
LCR II u otros totalizadores y en cambio se adicionaron funciones especcas, de
igual manera que cada funcin de desarrollo de manera personalizada.
A pesar de ser el totalizador personalizado, se tuvieron en cuenta funciones bsicas y
adicionales de tal manera que se pueda adaptar el sistema fcilmente a otros procesos
de descarga o aplicaciones similares.
Una de las desventajas fue que al querer realizar el totalizador de tal manera que
fuera de menor costo se usaron equipos igualmente de bajo costo que pudieran
realizar operaciones bsicas, esto sin embargo ocasion que al querer implementar
funciones ms complejas los equipos usados no estaban en la capacidad de hacerlo
o requera de procesos y tareas ms complejas de implementar. Un ejemplo claro
96
5.1 Resultados
de esto fue el querer realizar la implementacin de la compensacin de volumen por
norma API para lquidos y de la norma AGA reporte 8, puesto que para implementar
estas correcciones era necesario utilizar funciones matemticas las cuales no eran
soportadas por este tipo de PLC y para lograr implementarlas se necesitaba aplicar
otros mtodos de clculo como el uso del CORDIC, el cual es complejo y requiri de
estudio adicional por lo cual no se alcanz a implementar para la fecha de entrega.
Tanto el API como el AGA 8 son funciones que se pueden adicionar al sistema
como una mejora o una implementacin futura mas no estrictamente necesarias
puesto que estas hacen correccin del volumen medido y para el proceso solo se iba
a utilizar la totalizacin por masa. Pero de igual manera se gener un estudio a
fondo de este tipo de compensaciones y se desarrollaron algoritmos que pueden ser
fcilmente traducidos a otros lenguajes para PLC y dispositivos totalizadores para
su implementacin.
Entre las funciones adicionales que se implementaron y no eran requeridas por la
empresa estn: La totalizacin por volumen, visualizacin y conguracin de fecha y
hora, mltiples opciones de visualizacin de datos, calibracin del delta para cierre
de vlvula, impresin de tiquete de precio y lectura de canales analgicos.
5.1.2. Instalacin en Campo
La instalacin del sistema totalizador fue realizada en la zona de descarga y medicin
que se observa en la Figura 5.3, en donde los equipos como el sensor de coriolis y el
transmisor que previamente haban sido desmontados para ser llevados a calibracin
y para realizar pruebas con el equipo totalizador desarrollado, fueron instalados
nuevamente. Detras de la estacin donde se instalo el sistema se pueden observar
los tanques que almacenan el GLP.
Figura 5.3.: Zona de instalacin
97
Captulo 5 Captulo 5
El sensor de coriolis como se observa en la Figura 5.4 despues de ser calibrado se
instalo dentro de la lnea por la que circula el GLP en el punto donde antes se
encontraba instalado, mientras que el transmisor fue instalado de tal manera que
quedara junto al equipo totalizador como se observa en la Figura 5.5. El transmisor
al poder instalarse en interperie, cumpliendo con los grados de proteccin para su
funcionamiento en campo, se intalo sobre una placa la cual esta sujeta al sistema
totalizador desarrollado.
Figura 5.4.: Instalacin Coriolis
Figura 5.5.: Instalacin totalizador y transmisor RHE08
El transmisor inicialmente instalado fue el RHE08 de RHEONIK el cual entrega las
seales requeridas por el sistema totalizador para poder visualizar las varibles de el
98
5.1 Resultados
proceso: densidad, temperatura o ujo, adems de realizar el calculo del volumen
mediante la densidad entregada por este transmisor. Al instalar el transmisor en el
mismo lado que el sistema totalizador se da la opcin de manipular el sistema y
congurar los 2 elemtos sin tener que desplazarse a otro lado, esto con el n de que
sea mas comodo para el operador al momento de congurar el sistema.
Para el montaje del equipo totalizador y el transmisor fue necesario armar un soporte
como se observa en la Figura 5.6, esto fue realizado en el rea clasicada o de
proceso pero en un lugar con el suciente espacio y sin obstculos que permitiera la
manipulacin del equipo sin ningn problema para prever posibles accidentes.
Figura 5.6.: Soporte para el totalizador y transmisor
Luego una vez montados y cableados los equipos, se realizaron los respectivos ajustes
y calibraciones para su puesta en funcionamiento, dentro de estas pruebas realizadas
se comprueba que las bombas a controlar, las cuales se observan en la Figura 5.7,
tengan un correcto funcionamiento a momento de su accionamiento. Esta fue una
de las principales funciones que el sistema anterior no poseia, el activar los actua-
dores y desactivarlos al momento de realizr una descarga del producto. En la gura
Figura 5.7 se observan las bombas 1 y 2 las cuales son los actuadores principales
del proceso y hacen circular el GLP para realizar la decarga correspondiente a los
carrotanques.
99
Captulo 5 Captulo 5
Figura 5.7.: Bombas dentro del proceso
Una vez hechas las pruebas y que todo quedara funcionando correctamente, el sis-
tema se mostro a los operarios y se les enseo su funcionamiento para que puedan
empezar a realizar las descargas. De igual manera para comprobar la versatilidad del
equipo este es probado con el otro transmisor disponible el cual es el RHE12 como
se observa en la Figura 5.8, aunque este no entrega seales como la de densidad al
totalizador para poder calcular el volumen, sin embargo puede entregar el valor de
la temperatura o ujo del sistema los cuales pueden ser tiles para el operario.
Figura 5.8.: Totalizador y transmisor RHE12
Este transmisor al igual que el RHE 08 puede entregar la masa o volumen totalizado
por medio de una seal de pulsos.
100
5.1 Resultados
En la Figura 5.9 se observa el sistema totalizador en funcionamiento ya en un pro-
ceso de transferencia de custodia para el carrotanque el cual se puede observar en
la parte posterior de el sistema, de igual forma se ve ms claramente del posiciona-
miento del totalizador y del transmisor, posicin en la que es fcil su manipulacin
y se encuentra cerca a carrotanque que se desea cargar. Como se haba planeado ini-
cialmente la impresin y entrega del tiquete se realizo en una de las casetas aledaas
al lugar y por tanto la impresora no se instalo en el rea de descarga.
Figura 5.9.: Sistema en funcionamiento
El funcionamiento del sistem puede ser visto en la Figura 5.10 donde se el sistema se
encuentra realizando una descarga, la cual fue puesta para que se entregaran 6703.2
kg. El sistema al no haver llegado al valor deseado, los actuadores se mantienen
encedidos.
Figura 5.10.: Visualizacin del sistema en funcionamiento
Finalmente en la Figura 5.11 y Figura 5.12 se observan las comparaciones en la
medida del totalizador con el transmisor, en la Figura 5.11 se encuentra el totalizados
101
Captulo 5 Captulo 5
con el transmisor RHE08 y en la Figura 5.12 con el transmisor RHE12, en ambas
se ve la medida entregada tanto por totalizador y transmisor es la misma, en las
imgenes hay un pequeo desfase debido al tiempo de toma de las fotos puesto que
no se contaba sino con una sola cmara.
Figura 5.11.: Funcionamiento totalizador y RHE08
Figura 5.12.: Funcionamiento totalizador y RHE12
5.1.3. Estudios AGA reporte 8
Desarrollando el algoritmo con los valores de las tablas A.5-1 y A.5-2 e implemen-
tando las ecuaciones descritas en Subseccin 4.2.5 en Excel con base en los datos
entregados por MATLAB a partir de entradas simuladas de presin y temperatura
se obtienen las siguientes tablas y resultados para cada una de las mezclas de gases.
Estas mezclas de gases son las trabajadas en la norma AGA 8, la cual provee re-
sultados experimentales para determinadas presiones y temperaturas para ser com-
parados con los resultados generados por los algoritmos desarrollados, de esta com-
102
5.1 Resultados
paracin se genera igualmente anlisis de acuerdo a clculo de desviacin estndar,
erros de BIAS, RMS, entre otros.
Las siguientes tablas 5 tablas hacen referencia a 5 mezclas de gases diferentes (GULF
COAST, AMARILLO, EKOFISK, HIGH N2 Y HIGH CO2-N2) en donde se mues-
tran los diferentes valores de presin y temperatura para los cuales se realizaron los
clculos del factor de supercompresibilidad (Fpv), Factor de compresibilidad calcu-
lado por el algoritmo desarrollado en MATLAB (Zcalc), Factor de compresibilidad
experimental entregado por la norma AGA 8 para comparacin de resultados (Zexp)
y la diferencia (Zdi) entre el valor clculado (Zcalc) y el valor experimental de re-
ferencia (Zexp).
103
Captulo 5 Captulo 5
Cuadro 5.1.: GULF COAST
BIAS = 0,869575
AAD = 0,2483
Std.Dev = 1,1207
RMS = 1,4046
Max.Dev.= 4.4013
104
5.1 Resultados
Cuadro 5.2.: AMARILLO
BIAS = 0,883983
AAD = 0,2647
Std.Dev = 1,2277
RMS = 1,4972
Max.Dev.= 4.7202
105
Captulo 5 Captulo 5
Cuadro 5.3.: EKOFISK
BIAS = 1,281059
AAD = 0,4667
Std.Dev = 2,3453
RMS = 2,6400
Max.Dev.= 8.2121
106
5.1 Resultados
Cuadro 5.4.: HIGH N2
BIAS = 0,531931
AAD = 0,1580
Std.Dev = 0,7298
RMS = 0,8938
Max.Dev.= 2.8301
107
Captulo 5 Captulo 5
Cuadro 5.5.: HIGH CO2-N2
BIAS = 0,826496
AAD = 0,2556
Std.Dev = 1,2051
RMS = 1,4457
Max.Dev.= 4.7728
108
5.1 Resultados
De los resultados anteriores podemos observar que el error de los valores del factor
de compresibilidad calculados con respecto a los valores experimentales es muy bajo
y que el algoritmo desarrollado tiene muy buena precisin.
El error que se puede generar en los valores calculados puede ser principalmente por
el uso de decimales puesto que en el clculo se implementan algunas conversiones de
unidades, constantes con decimales y adicionalmente a eso el mtodo de clculo de
MATLAB tiende a redondear los decimales perdiendo un poco la presin y causando
quizs estas pequeas diferencias con respecto a los valores experimentales.
Se observa de todos modos un error muy pequeo lo que concluye que el algoritmo
implementado y las ecuaciones usadas para el clculo del factor de compresibilidad
son muy precisos y conables por lo cual pueden ser perfectamente implementados
en dispositivos totalizadores o PLCs.
Lo anterior fue para comprobar el algoritmo implementado en ciertas mezclas de
gases, ahora se muestra a continuacin el valor del factor de compresibilidad para
cada uno de los gases que se muestran en la tabla 5 del AGA reporte 8, esto con
el n de generar una mayor profundizacin en el uso del mtodo de caracterizacin
detallado y observar el comportamiento de cada uno de estos gases a determinada
presin y temperatura.
Para los valores de temperatura y humedad se seleccionaron los ya usados para los
clculos anteriores, es decir, para la temperatura se tienen los cuatro valores (32F,
50F, 100F y 130F), mientras para la presin se tiene un rango de 14.73 psia a
1200 psia.
Los valores del factor de compresibilidad se muestran en varias grcas cada una
para cada gas donde el eje X es la presin y el eje Y es el valor del factor de
compresibilidad, las cuatro respuestas generadas son para cada valor de temperatura
(32F, 50F, 100F y 130F).
Los gases del 11 al 18 de la tabla 5 del AGA reporte 8 corresponden especialmente
a fracciones de gases que se encuentran presentes en algunas mezclas de gases por
lo cual no son evaluadas de forma individual para el clculo del factor Z.
De acuerdo a lo anterior las grcas a continuacin muestran el comportamiento
de la compresibilidad del gas de acuerdo a las temperaturas y presiones dadas pa-
ra los gases de la tabla 5 del AGA reporte 8 a excepcin de los ya mencionados
anteriormente.
109
Captulo 5 Captulo 5
Figura 5.13.: Metano
La grca muestra que el metano desde las temperaturas y presiones bajas, hasta las
temperaturas y presiones altas el gas en fcilmente compresible debido a que el valor
de Z siempre es menor a 1 y a medida que se aumente la presin o la temperatura
el gas es ms susceptible a la compresin.
Figura 5.14.: Nitrgeno
El nitrgeno a temperatura y presin baja se comporta ms como un gas ideal sin
embargo su comportamiento vara mucho en presiones y temperaturas altas. Como
se observa el nitrgeno a presiones altas pero con baja temperatura el gas se vuelve
fcilmente compresible (Z<1), sin embargo si en estas presiones altas la temperatura
tambin es alta el gas se vuelve difcilmente compresible (Z>1).
110
5.1 Resultados
Figura 5.15.: Dixido de carbn
El dixido de carbn a temperatura y presin baja es un gas con comportamiento
ideal y al ir aumentando la temperatura o presin el gas es fcilmente compresible
(Z<1).
Figura 5.16.: Etano
Al igual que el dixido de carbn o el metano, el etano es ms fcilmente compresible
a medida que la temperatura y la presin aumentan (Z<1).
temperatura (Z>1).
111
Captulo 5 Captulo 5
Figura 5.17.: Propano
El propano es un gas que como se observa tiene un comportamiento particular, en
donde hasta cierto punto de presin (alrededor de 3 MPa) tiene un comportamiento
muy similar al de un gas ideal (Z=1), pero si despus de este punto de presin la
temperatura es baja el gas es fcilmente compresible (z<1), pero si se tiene altas
temperaturas el gas mantiene un comportamiento de gas ideal.
Figura 5.18.: 6. Agua (Vapor)
Inicialmente a simple vista se podra decir que el vapor de agua tiene un compor-
tamiento similar al del propano, sin embargo al observar con ms detalle se logra
observar que los valores de Z se encuentran el orden de 104 es decir valores que no
112
5.1 Resultados
corresponde a un valor de factor de compresibilidad, esto es debido a que en los
valores de temperatura evaluados el agua an no es vapor por lo cual las formulas
aplicadas y el algoritmo por tanto no aplican y debera cambiarse este rango de
temperatura para evaluar el vapor de agua.
Figura 5.19.: 7. Sulfato de hidrgeno
Comportamiento similar a la mayora de gases en donde a mayor temperatura o
presin el gas es fcilmente compresible (Z<1).
Figura 5.20.: 8. Hidrgeno
Este tiene gas tiene comportamiento opuesto a la mayora de los gases hasta ahora
113
Captulo 5 Captulo 5
donde a mayor presin o temperatura el gas es ms difcilmente compresible (Z>1).
9. Monxido de carbono
El monxido de carbono tiene un comportamiento similar al ya visto en el nitrgeno
en donde a presiones altas y temperaturas bajas el gas es fcilmente compresible
(Z<1), pero en presiones altas y temperaturas altas el gas es difcilmente compresible
(Z>1).
Figura 5.21.: 10 Oxgeno
El oxgeno posee comportamiento de la mayora de gases vistos, en donde a mayor
temperatura o mayor presin el gas es fcilmente compresible (Z>1).
Figura 5.22.: 20. Helio
114
5.1 Resultados
Este gas se comporta tambin de manera inversa a la mayora de gases y se vuelve
difcilmente compresible a medida que aumenta la presin o la temperatura (Z>1).
Figura 5.23.: 21. Argn
Este tipo de gas tiene comportamiento tpico y se vuelve menos compresible a medida
que la presin o la temperatura aumentan (Z<1).
Figura 5.24.: 11. i- Butano
115
Captulo 5 Captulo 5
Figura 5.25.: 13. i-pentano
Como se puede observar en las dos anteriores grcas tanto el i-butano como el
i-pentano poseen con el clculo por medio del mtodo de caracterizacin detallada
unos valores de Z que se encuentran fuera de rango, esto es debido a que estos tipos
de gases al igual que los gases de la tabla 5 del AGA reporte 8 del 11 al 18 hacen
referencia a fracciones de gas dentro de una mezcla de gas y por tal motivo su factor
de compresibilidad individual no puede ser hallado por este mtodo y por tanto el
algoritmo implementado arroja estos datos fuera de rango para el valor de Z.
116
6. CONCLUSIONES
El sistema totalizador desarrollado puede ser aplicado para procesos de scalizacin
de masa, donde al utilizar un sensor de tipo coriolis no se hace necesario realizar
alguna correccin sobre esta medida. El rango de productos Oil and Gas en que puede
ser aplicado este sistema est determinado por el tipo de sensor coriolis el cual debe
estar diseado para procesos de transferencia de custdia. El sistema ademas de
mostrar en un pantalla HMI el valor totalizado por el transmisor de ujo, realiza un
proceso de registro de las operaciones realizadas, en donde por medio de la impresin
de recibos y el desarrollo de una base de datos se logro llevar un control de el proceso.
El sistema puede ser extendido para la realizacin de transferencia de custdia por
volumen, por medio de la utilizacin de una correccin de volumen descrita en el
API MPMS 11.2.4 la cual fu estudiada para su aplicacin. El sistema totalizador
desarrollado puede mostrar en la pantalla HMI el volumen de la descarga por medio
de un clculo simple utilizando el valor de densidad leda por el equipo. Este valor
de volumen para transferencia de custdia debe ser corregido a condiciones base.
Para la implementacin de la correccin de volumen en el sistema actual se puede
desarrollar mediante la utilizacin de mtodos como el CORDIC para el clculo de
funciones las cuales son necesarias para la aplicacin del algoritmo de correccin.
Aunque es implemetable la correccin de volumen, esta solo est descrita para l-
quidos, en este caso para GLP . Para lo correccin de volumen en gases es necesario
aplicar el reporte 8 de AGA el cual fue tambin estudiado en este trabajo. Para
la implementacin del AGA 8 y API MPMS 11.2.4 se necesita de un PLC el cual
soporte el ccculo de funciones como la exponencial, logaritmos y potencias racio-
nales. Para un mejor desarrollo de este tipo de algoritmos es recomendable el uso
de PLCs los cuales puedan ser programados mediante lenguaje estructurado, en el
cual se pueden realizar plantillas de los algoritmos descritos. Para la correccin de
volumen de otro tipo de lquido ser necesario referirse a otro tipo de tabla descrita
en el API MPMS
En procesos de scalizacin y transferencia de custdia normalmente se usan regis-
tadores y computadores de ujo los cuales son desarrollados en el exterior. Con este
trabajo se logra hacer una aproximacin al desarrollo de sistemas que sirvan para
la realizacin de procesos scales en Colombia, donde estos sistemas pueden poseer
un menor costo y pueden ser adaptados a necesidades de las empresas, como es el
caso el manejo de actuadores para controlar el sistema de forma automomtica.
El proceso de descarga y transferencia de custodia mejora con la utilizacin del
equipo puesto que se tiene mayor control sobre la medida del uido y se eliminan
117
Captulo 6 CONCLUSIONES
procesos que se realizaban de manera manual y los cuales no eran convenientes para
este tipo de proceso debido a que generaban mucha incertidumbre en la descarga
que podran generar prdidas.
El clculo del factor de coompresibilidad basados en el AGA reporte 8 a pesar de
requerir de gran cantidad de ecuaciones y procedimientos, es un proceso que se puede
trabajar de tal manera que el algoritmo desarrollado pueda ser traducido a otros
lenguajes de programacin y pueda ser implementado en gran variedad de PLCs
o dispositivos que se quierean utilizar para realizar una totalizacin para gases, sin
embargo se deben tener en cuenta la composicin que tiene la mezcla de gas con la
que se trabaja, puesto que su abilidad depende en gran medida de los porcentajes
exactos de cada gas que componen la mezcla, por lo tanto en una aplicacin real
sera necesario un anlisis de la mezcla por medio de un cromatgrafo para obtener
estos datos antes de la implementacin del AGA 8.
El algoritmo implementado en MATLAB gener buenos resultados de acuerdo a
las comparaciones realizadas con los resultados experimentales suministrados por la
norma AGA 8, de tal manera que se logr observar la gran precisin y abilidad
que este tiene comparado con los resultados de los datos experimentales obtenidos
por el AGA, adems su codigo aunque extenso es de fcil traduccin a lenguajes
de PLC para que pueda ser implementada la correccin en una mejora al sistema o
aplicaciones futuras de otros totalizadores.
118
Bibliografa
[1] R. Cobo, El abc de la automatizacin, tech. rep., Asociaciin De La Industria
Elctrica Electrnica, Chile.
[2] Y. Mori and M. Shitara, The resonance of coriolis mass owmeters in a pul-
sating ow, in SICE 2003 Annual Conference, vol. 1, pp. 1014 Vol.1, 2003.
[3] Rocatek, Ejemplo plc. disponible en: http://www..rocatek.com/images- /fo-
rum/plc_desc.jpg., [Accedido: 30-Abril-2013].
[4] RHEONIK., Operation manual rheonik mass owmeter rhe 07, 08, 11. rhm..nt,
etx, ht, 2004.
[5] AutomationDirect., Dl05/06 plc - dl0506select.pdf. disponible en:
http://www.automationdirect.com/static/specs/dl0506select.pdf., [Acce-
dido: 01-may-2013].
[6] AutomationDirect., Dl05/06 option modules user manual, Julio de 2011.
[7] American Petroleum nstitute, Manual of petroleum measurement standards
chapter 11-physical properties data. section2, part 4-temperature correction for
the volume od NGL and LPG. tables 23E, 24E, 53E 54E, 59E, and 60E, 2007.
[8] Liquid Control Groups, LectroCount LCR-ii electronic register, set-up and
operation manual, publicacn no. 500302, Feb. 2006.
[9] R. A. L. M., Actualizacin Del Mdulo Oleo Del Sistema Centinela Aplicando
La Norma Api Para El Clculo De Cantidades De Petrleo, En Los Tanques
De Almacenamiento Presentes En Los Patios De Tanques Y Terminales De
Embarque Del Oriente De Venezuela. Ingeniera, Universidad De Oriente Ncleo
De Anzotegui, Escuela De Ingeniera Y Ciencias Aplicadas, Departamento De
Petrleo., Barcelona, Espaa, Nov. 2011.
[10] C. F. Wilson, Disminucin Del ndice De Prdidas En Diferentes reas De La
Gerencia Complejo Barrancabermeja (GCB) De ECOPETROL A Partir Del
Mejoramiento De La Medicin Esttica Y Dinmica. Ingeniera, Universidad
Industrial de Santander, Facultad de Ingenieras Fsico-Mecnicas, Escuela de
Ingenieras Elctrica, Electrnica y de Telecomunicaciones, Bucaramanga, Co-
lombia, 2006.
[11] N. R. B. Mrquez, Diseo De Un Sistema De Descarga De Crudo Pesado Desde
Carrotanques Para Campo Rubiales, Meta En La Empresa Ingeniera STRY-
CO, LTDA. Ingeniera, Universidad Industrial de Santander, Facultad de In-
119
Bibliografa
genieras Fsico-Mecnicas, Escuela de Ingeniera Mecnica, Bucaramanga, Co-
lombia, 2010.
[12] D. F. B. Patarroyo, Implementacin De Un Sistema De Medicin Y Fiscaliza-
cin De Custodia De Crudo En La Estacin De Transferencia Moriche. Inge-
niera, Universidad Industrial de Santander, Colombia. Facultad de Ingenieras
Fsico-Qumicas. Ingeniera de Petrleos, Bucaramanga, Colombia.
[13] D. P. B. Salazar, Diseo De Un Sistema De Medicin Y Control De Uso De
Combustibles Para Las Centrales De Generacin Termoelctrica Del S. N. I.
Ingeniera, Escuela Politcnica Nacional. Escuela de Ingeniera, Quito, Ecuador,
Nov. 2009.
[14] K. C. Ruiz, Diseo De Un Sistema De Medicin Automtico En Lnea De Pro-
duccin En Las Principales Estaciones De Flujo De La Unidad De Explotacin
Tierra Oeste (PDVSA). Ingeniera, Universidad de los Andes. Facultad de In-
geniera. Escuela de Sistema. Departamento de Sistemas de Control, Repblica
Bolivariana de Venezuela, July 2006.
[15] K. C. Ruiz, Interfaz hombre-mquina, tech. rep., Unicauca, Colombia, Sept.
2009.
[16] Ingeniera en Automatizacin y Control, Introcuccin a HMI (interfaz hombre
mquina), laboratorio 2, Universidad Nacional de Quilmes, Argentina.
[17] A. C. Maza Cabello, Dise o de un sistema experto para el enderezado de chasis
en fro. Computacin, electrnica y mecatrnica, Universidad de las Amricas
Puebla, Escuela de Ingeniera, Mxico, 2009.
[18] Osinergim, Propiedades y caractersticas del GLP, Sept. 2011.
[19] NEMA, NEMA enclosure types, 2013.
[20] I. I. en Microcontroladores, Protocolo RS-485, 2013.
[21] O. E. S.A., Guia rapida reles de estado solido ssrs, Accedido en Junio de
2013.
[22] J. Seiko Epson Corporation, Nagano, Faq for esc/pos, 2002.
[23] Soarland and H. t. Hexin, Hxsp-485 rs-232/rs-485 converter user guide, Ac-
cedido en Mayo de 2013.
[24] J. S. K.e. Starling, Compressibility Factors of Natural Gas and Other Related
Hydrocarbon Gases. PhD thesis, Transmision Measurement Committee Report
No.8, American Petroleum Institute MPMS Chapter 14.3, GasbResearch Insti-
tute, Nov. 2003.
120
Nomenclatura
, , Parmetros de la ecuacin cuadrtica
Interpolacin variable
Nmero de topos del gas
60,high
Lmite superior para la densidad relativa observada del uido a 60
o
F
60,low
Lmite inferior para la densidad relativa observada del uido a 60
o
F
60,mid
Lmite intermedio para la densidad relativa observada del uido a 60
o
F
60,trial
Densidad relativa trial a 60
o
F
60
Densidad relativa a la temperatura base de 60
o
F
TB
Densidad relativa a la temperatura base, T
B
Tx
Densidad relativa ala temperatura observada, T
x
x,high
Densidad relativa a una temperatura observada correspondiente a el lmite
superior de 60
o
F densidad relativa
x,low
Densidad relativa a una temperatura observada correspondiente a el lmite
inferior de 60
o
F densidad relativa
x,mid
Densidad relativa a una temperatura observada correspondiente a el interme-
dio de 60
o
F densidad relativa
x,trial
Densidad relativa trial a la temperatura observada
x
Densidad relativa a una temperatura observada
sat
Densidad molar de saturacin (grammole/L)
15
Densidad a una temperatura base de 15
o
C (kg/m
3
)
20
o
C (kg/m
3
)
60
o
F (kg/m
3
)
121
Nomenclatura
sat
60
Densidad molar de saturacin a 60
o
F (grammole/L)
T
60
Densidad molar de saturacin a la temperatura observada (grammole/L)
c
Densidad molar crtica (grammole/L)
Parmetro en la ecuacin de densidad de saturacin
A, B, C Parmetros de la ecuacin cuadratica
a
n
Constante en la tabla 4 del AGA reporte 8
API American Gas Association
B Coeciente en segundo lugar virial
B
*
nij
coeciente binario de la caracterizacin
b
n
c
n
C
n
Coecientes que son funciones de la composicin
C
TL
Factor de correccin de temperatura
D Densidad del gas reducida
d densidad molar (moles por unidad de volumen)
DTR Input Data Terminal Ready
E
ij
Parmetro binario de la energa del coeciente en segundo lugar virial
E
ij
Parmetro binario de la interaccin de la energa en segundo lugar virial del
coeciente (tabla 6)
E
i
Parmetro caracterstico de la energa para el componente i-simo (tabla 5)
E
j
Parmetro caracterstico de la energa para el componente j-simo (tabla 5)
F Parmetro de la temperatura alta de la mezcla
F
i
Parmetro de alta temperatura para el componente del i-simo (tabla 5)
F
j
Parmetro de alta temperatura para el componente del j-simo (tabla 5)
f
n
G Parmetro de la orientacin
122
Nomenclatura
G
ij
Parmetro binario de la orientacin
G
ij
Parmetro binario de la interaccin para la orientacin (tabla 6)
G
i
Parmetro de la orientacin para el componente del i-simo (tabla 5)
G
j
Parmetro de la orientacin para el componente del j-simo (tabla 5)
g
n
GLP Gas Licuado de Petrleo
GNL Gas Natural Licuado
h
2
Factor de escalacin
HMI Human Machine Interface
IEC International Electro-technical Commission
IP International Protection
K Parmetro del tamao
k
1
, k
2
, k
3
, k
4
Parmetros en la ecuacin de densidad saturacin
K
ij
Parmetro binario de la interaccin para el tamao (tabla 6)
K
ij
Parmetro binario de la interaccin para la orientacin (tabla 6)
K
i
Parmetro del tamao para el componente del i-simo (tabla 5)
K
j
Parmetro del tamao para el componente del j-simo (tabla 5)
k
n
M
r
Masa molar (peso molecular)
M
r
i Masa molar del i-simo componente
MPMS Measurement Petroleum Manual Standard
N Nmero de componentes en mezcla del gas
n Nmero de moles del gas
NEMA National Electrical Manufacturers Association
P Presin esttica absoluta del gas
123
Nomenclatura
p Densidad total (masa por unidad de volumen)
P
b
Presin absoluta de la referencia del contrato
p
b
Densidad total en la condicin de referencia TbPb
PLC Programmable Logic Controler
Q Parmetro quadrupole
Q
i
Parmetro quadrupole para el componente del i-simo
Q
j
Parmetro quadrupole para el componente del j-simo
q
n
Constante (n = 1.2, etc.) en la tabla 4 del AGA reporte 8
R Constante de los gases.
RTS Input Request To Send
S
i
Parmetro del dipolo para el componente del i-simo (tabla 5)
S
j
Parmetro del dipolo para el componente del j-simo (tabla 5)
s
n
Constante (n = 1.2, etc.) en la tabla 4 del AGA reporte 8
SSR Solid State Rele
T Temperatura absoluta del gas
T
BK
Temperatura base (288.15K o 293.15K)
T
B
Temperatura base (60
o
F, 15
o
C o 20
o
C)
T
b
Temperatura absoluta de la referencia
T
c,ref
Temperatura de referencia crtica del uido (K)
T
c
Temperatura crtica del uido (K)
T
F
Medicin de la temperatura observada (
o
F o
o
C)
T
r,x
Temperatura observada reducida
T
x
Temperatura observada (K)
U Parmetro de la energa de la mezcla
U
ij
Parmetro binario de la interaccin para la energa conformada (tabla 6)
124
Nomenclatura
V Volumen del gas
V
60
/V
Tx
Razn de volumen a 60
o
F
W
i
Parmetro de la asociacin para el componente del i-simo (tabla 5)
W
j
Parmetro de la asociacin para el componente del j-simo (tabla 5)
w
n
X Factor de interpolacin
x
i
Fraccin de mol del i-simo componente en la mezcla del gas
x
i
Fraccin de moles del componente i en mezcla del gas
x
j
Fraccin de mol del j-simo componente en la mezcla del gas
Z Factor de compresibilidad del gas
Z
b
Factor de compresibilidad en condiciones de referencia del Tb Pb
Z
c
Factor de conpresibilidad crtico
125
Anexos
127
A. PLANO P&ID
129
Captulo A PLANO P&ID
130
B. PLANO UNIFILAR
131
Captulo B PLANO UNIFILAR
132
C. PLANO ELCTRICO
133
Captulo C PLANO ELCTRICO
134
D. CONEXIONES DE
TRANSMISION DE DATOS
135
Captulo D CONEXIONES DE TRANSMISION DE DATOS
136
E. COMPARACIN DE
FUNCIONES ENTRE
TOTALIZADORES
137
Captulo E COMPARACIN DE FUNCIONES ENTRE TOTALIZADORES
138
F. TABLAS Y ALGORITMO DEL
AGA REPORTE 8
139
Captulo F TABLAS Y ALGORITMO DEL AGA REPORTE 8
Cuadro F.1.: Tabla 4 del AGA reporte 8
140
TABLAS Y ALGORITMO DEL AGA REPORTE 8
Cuadro F.2.: Continuacin Tabla 4 del AGA reporte 8
141
Cuadro F.3.: Tabla 5 AGA reporte 8
142
Cuadro F.4.: Tabla 6 AGA reporte 8
143
Cuadro F.5.: Continuacin Tabla 6 AGA reporte 8
144
Cuadro F.6.: Tabla A.5-1 AGA reporte 8
145
Cuadro F.7.: Tabla A. 5-2, AGA reporte 8
146
Figura F.1.: Diagrama de ujo para el algortmo del clculo del factor de compre-
sibilidad Z
147
Algoritmo F.1 AGA reporte 8 en MATLAB
% AGA8.m
% P = Presin absoluta
% T = Temperatura absoluta
% N = Nmero de elementos
% R = Constante de los gases J/mol-K
% Entradas principales del sistema
%---------------------------------------------------------------------
N=input(Nmero de gases en la mezcla:);
P=[14.73 100 200 400 600 800 1000 1200]; % Valores de presin
T=[0 10 37.778 54.444]; % Valores de temperatura
V=1000; % V en dm^3
R=0.00831451; % R en dm^3*MPa/mol-K
T=T+273.15; % T en K (Kelvin)
for k=1:8
P(k)=P(k)*0.006894759086775369; % P en MPa (Megapascales)
end
%---------------------------------------------------------------------
% Tablas 4, 5 y 6
%---------------------------------------------------------------------
%Constantes del sistema
% Tabla 4 del AGA reporte 8
an = [0.1538326 1.341953 -2.998583 -0.04831228 0.3757965 -1.589575
-0.05358847 0.88659463 -0.71023704 -1.471722 1.32185035 -0.78665925
2.29129e-9 0.1576724 -0.436386 -0.04408159 -0.003433888 0.03205905
0.02487355 0.07332279 -0.001600573 0.6424706 -0.4162601 -0.06689957
0.2791795 -0.6966051 -0.002860589 -0.008098836 3.150547 0.007224479
-0.7057529 0.5349792 -0.07931491 -1.418465 -5.99905e-17 0.1058402
0.03431729 -0.007022847 0.02495587 0.04296818 0.7465453 -0.2919613
7.294616 -9.936757 -0.005399808 -0.2432567 0.04987016 0.003733797
1.874951 0.002168144 -0.6587164 0.000205518 0.009776195 -0.02048708
0.01557322 0.006862415 -0.001226752 0.002850908];
bn = [1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 3 3 3 3
3 3 3 3 3 3 4 4 4 4 4 4 4 5 5 5 5 5 6 6 7 7 8 8 8 9 9];
cn = [0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 0 1 1 1
1 1 1 1 1 1 0 0 1 1 1 1 1 0 1 1 1 1 0 1 0 1 1 1 1 1 1];
kn = [0 0 0 0 0 0 0 0 0 0 0 0 3 2 2 2 4 4 0 0 2 2 2 4 4 4 4 0 1 1 2
2 3 3 4 4 4 0 0 2 2 2 4 4 0 2 2 4 4 0 2 0 2 1 2 2 2 2];
un = [0 0.5 1 3.5 -0.5 4.5 0.5 7.5 9.5 6 12 12.5 -6 2 3 2 2 11 -0.5
0.5 0 4 6 21 23 22 -1 -0.5 7 -1 6 4 1 9 -13 21 8 -0.5 0 2 7 9 22 23
1 9 3 8 23 1.5 5 -0.5 4 7 3 0 1 0];
gn = [0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 0 0
1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 1 0 0];
qn = [0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0
0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 1 0 0 1 0 0 0 0 0 1];
fn = [0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0
0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
sn = [0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
148
wn = [0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
%---------------------------------------------------------------------
%Caracterizacin de parmetros
Mrn = [16.043 28.0135 44.01 30.07 44.097 18.0153 34.082 2.0159 28.01
31.9988 58.123 58.123 72.15 72.15 86.177 100.204 114.231 128.258
142.285 4.0026 39.948];
En = [151.3183 99.73778 241.9606 244.1667 298.1183 514.0156 296.355
26.95794 105.5348 122.7667 324.0689 337.6389 365.5999 370.6823 402.636293
427.72263 450.325022 470.840891 489.558373 2.610111 119.6299];
Kn = [0.4619255 0.4479153 0.4557489 0.5279209 0.583749 0.3825868
0.4618263 0.3514916 0.4533894 0.4186954 0.6406937 0.6341423 0.6738577
0.6798307 0.7175118 0.7525189 0.784955 0.8152731 0.8437826 0.3589888
0.4216551];
Gn = [0 0.027815 0.189065 0.0793 0.141239 0.3325 0.0885 0.034369
0.038953 0.021 0.256692 0.281835 0.332267 0.366911 0.289731 0.337542
0.383381 0.423754 0.469659 0 0];
Qn = [0 0 0.69 0 0 1.06775 0.633276 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
Fn = [0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0];
Sn = [0 0 0 0 0 1.5822 0.39 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
Wn = [0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
%---------------------------------------------------------------------
%Parmetro de valores de interaccin binaria
load TEij.txt;
load TUij.txt;
load TKij.txt;
load TGij.txt;
TE = TEij(1:21,1:21);
TU = TUij(1:21,1:21);
TK = TKij(1:21,1:21);
TG = TGij(1:21,1:21);
%---------------------------------------------------------------------
% Nmero de referencia del componente
% 1-Metano 2-Nitrogeno 3-Dixido de carbn 4-Etano 5-Propano 6-Agua
% 7-Sulfato de hidrgeno 8-Hidrgeno 9- Monxido de carbono 10-xigeno
% 11-i-Butano 12-n-Butano 13-i-Pentano 14-n-Pentano 15-n-Hexano
% 16-n-Heptano 17-n-Octano 18-n-Nonano 19-n-Decano 20-Helio 21-Argn
%---------------------------------------------------------------------
switch N
case 0,
disp(N debe ser mayor que 0)
otherwise
% Inicializacin de variables
F=0; Q=0; G1=0; G2=0; G=0;
U4=0; U3=0; U2=0; U1=0; U=0; K4=0; K3=0; K2=0; K1=0; K=0;
n=0; count=0; xi=0; xj=0; CID=0; Mr=0; i=0; j=0; k=0; l=0;
Gij=0; Eij=0; Gb=0; Qb=0; Fb=0; Sb=0; Wb=0;
Bnij=0; BK=0; B=0; D=0; d=0; Bb=0; Db=0; db=0;
B1=0; B2=0; B3=0; B4=0; B5=0; B6=0; B7=0; B8=0; B9=0;
B10=0; B11=0; B12=0; B13=0; B14=0; B15=0; B16=0; B17=0; B18=0;
149
Bb1=0; Bb2=0; Bb3=0; Bb4=0; Bb5=0; Bb6=0; Bb7=0; Bb8=0; Bb9=0;
Bb10=0; Bb11=0; Bb12=0; Bb13=0; Bb14=0; Bb15=0; Bb16=0;
Bb17=0; Bb18=0;
Z=0; Z1=0; Z2=0; Z3=0; Z4=0; Z5=0;
Zb=0; Z1b=0; Z2b=0; Z3b=0; Z4b=0; Z5b=0;
% Introduccin de cada uno de los elementos y fracciones de
% estos elementos que componen el gas
for count=1:N
CID(count)=input(Nmero del tipo de elemento:);
xi(count)=input(Fraccin molar del elemento:);
end
xt=sum(xi);
while(xt>1)
for count=1:N
CID(count)=input(Nmero del tipo de elemento:);
xi(count)=input(Fraccin molar del elemento:);
end
xt=sum(xi);
end
% Normalizacin de las fracciones de mol del gas
for n=1:N
xi(n)=(xi(n))/(xt);
xj(n)=xi(n);
end
% Clculo del factor de compresibilidad para cada valor de T y P
for k=1:8
for l=1:4
% Reiniciar valores
F=0; Q=0; G1=0; G2=0; G=0;
U4=0; U3=0; U2=0; U1=0; U=0; K4=0; K3=0; K2=0; K1=0; K=0;
n=0; Mr=0; i=0; j=0;
Gij=0; Eij=0; Gb=0; Qb=0; Fb=0; Sb=0; Wb=0;
Bnij=0; BK=0; B=0; D=0; d=0; Bb=0; Db=0; db=0; p=0;
B1=0; B2=0; B3=0; B4=0; B5=0; B6=0; B7=0; B8=0; B9=0;
B10=0; B11=0; B12=0; B13=0; B14=0; B15=0; B16=0;
B17=0; B18=0;
Bb1=0; Bb2=0; Bb3=0; Bb4=0; Bb5=0; Bb6=0; Bb7=0;
Bb8=0; Bb9=0;
Bb10=0; Bb11=0; Bb12=0; Bb13=0; Bb14=0; Bb15=0;
Bb16=0; Bb17=0; Bb18=0;
Z1=0; Z2=0; Z3=0; Z4=0; Z5=0;
Z1b=0; Z2b=0; Z3b=0; Z4b=0; Z5b=0;
% Clculo de cada uno de los parmetros
for n=1:N
F=(F)+((xi(n))^(2)*Fn(CID(n)));
Q=(Q)+((xi(n))*(Qn(CID(n))));
G1=(G1)+((xi(n))*(Gn(CID(n))));
U3=(U3)+((xi(n))*(En(CID(n)))^(5/2));
K3=(K3)+((xi(n))*(Kn(CID(n)))^(5/2));
end
for i=1:N-1
for j=(i+1):N
150
G2=(G2)+((xi(i))*(xj(j))*(((TG(CID(j),CID(i)))-1)*(Gn(CID(i))+
Gn(CID(j)))));
U2=(U2)+((xi(i))*(xj(j))*((((TU(CID(j),CID(i)))^(5))-1)*(En(CID(i))*
En(CID(j)))^(5/2)));
K2=(K2)+((xi(i))*(xj(j))*((((TK(CID(j),CID(i)))^(5))-1)*(Kn(CID(i))*
Kn(CID(j)))^(5/2)));
end
end
n=0; i=0; j=0;
U4=(U3)^(2);
U1=U4+(2*U2);
U=(U1)^(1/5);
K4=(K3)^(2);
K1=K4+(2*K2);
K=(K1)^(1/5);
G=G1+G2;
% Clculo del segundo coeficiente virial B
for i=1:N
for j=1:N
Eij(i,j)=(TE(CID(j),CID(i)))*(((En(CID(i)))*(En(CID(j))))^(1/2));
Gij(i,j)=((TG(CID(j),CID(i)))*((Gn(CID(i)))*(Gn(CID(j)))))/2;
BK(i,j)=(xi(i))*(xj(j))*(((Kn(i))*(Kn(j)))^(3/2));
B1=(B1)+(BK(i,j));
B2=(B2)+((BK(i,j))*((Eij(i,j))^(1/2)));
B3=(B3)+((BK(i,j))*(Eij(i,j)));
B4=(B4)+((BK(i,j))*((Eij(i,j))^(un(4))));
B5=(B5)+((BK(i,j))*((Eij(i,j))^(un(5)))*(Gij(i,j)));
B6=(B6)+((BK(i,j))*((Eij(i,j))^(un(6)))*(Gij(i,j)));
B7=(B7)+((BK(i,j))*((Eij(i,j))^(un(7)))*((Qn(i))*(Qn(j))));
B8=(B8)+((BK(i,j))*((Eij(i,j))^(un(8)))*((Sn(i))*(Sn(j))));
B9=(B9)+((BK(i,j))*((Eij(i,j))^(un(9)))*((Sn(i))*(Sn(j))));
B10=(B10)+((BK(i,j))*((Eij(i,j))^(un(10)))*((Wn(i))*(Wn(j))));
B13=(B13)+((BK(i,j))*((Eij(i,j))^(un(13)))*((Fn(i))*(Fn(j))));
B14=(B14)+((BK(i,j))*((Eij(i,j))^(un(14))));
B15=(B15)+((BK(i,j))*((Eij(i,j))^(un(15))));
B16=(B16)+((BK(i,j))*((Eij(i,j))^(un(16)))*((Qn(i))*(Qn(j))));
B17=(B17)+((BK(i,j))*((Eij(i,j))^(un(17))));
B18=(B18)+((BK(i,j))*((Eij(i,j))^(un(18))));
end
j=0;
end
% Coeficientes a temperatura base
Bb1=(B1)*((an(1))*((288.71)^(-un(1))));
Bb2=(B2)*((an(2))*((288.71)^(-un(2))));
Bb3=(B3)*((an(3))*((288.71)^(-un(3))));
Bb4=(B4)*((an(4))*((288.71)^(-un(4))));
Bb5=(B5)*((an(5))*((288.71)^(-un(5))));
Bb6=(B6)*((an(6))*((288.71)^(-un(6))));
Bb7=(B7)*((an(7))*((288.71)^(-un(7))));
Bb8=(B8)*((an(8))*((288.71)^(-un(8))));
151
Bb9=(B9)*((an(9))*((288.71)^(-un(9))));
Bb10=(B10)*((an(10))*((288.71)^(-un(10))));
Bb11=(B11)*((an(11))*((288.71)^(-un(11))));
Bb12=(B12)*((an(12))*((288.71)^(-un(12))));
Bb13=(B13)*((an(13))*((288.71)^(-un(13))));
Bb14=(B14)*((an(14))*((288.71)^(-un(14))));
Bb15=(B15)*((an(15))*((288.71)^(-un(15))));
Bb16=(B16)*((an(16))*((288.71)^(-un(16))));
Bb17=(B17)*((an(17))*((288.71)^(-un(17))));
Bb18=(B18)*((an(18))*((288.71)^(-un(18))));
% Coeficientes para el factor de compresibilidad
B1=(B1)*((an(1))*((T(l))^(-un(1))));
B2=(B2)*((an(2))*((T(l))^(-un(2))));
B3=(B3)*((an(3))*((T(l))^(-un(3))));
B4=(B4)*((an(4))*((T(l))^(-un(4))));
B5=(B5)*((an(5))*((T(l))^(-un(5))));
B6=(B6)*((an(6))*((T(l))^(-un(6))));
B7=(B7)*((an(7))*((T(l))^(-un(7))));
B8=(B8)*((an(8))*((T(l))^(-un(8))));
B9=(B9)*((an(9))*((T(l))^(-un(9))));
B10=(B10)*((an(10))*((T(l))^(-un(10))));
B11=(B11)*((an(11))*((T(l))^(-un(11))));
B12=(B12)*((an(12))*((T(l))^(-un(12))));
B13=(B13)*((an(13))*((T(l))^(-un(13))));
B14=(B14)*((an(14))*((T(l))^(-un(14))));
B15=(B15)*((an(15))*((T(l))^(-un(15))));
B16=(B16)*((an(16))*((T(l))^(-un(16))));
B17=(B17)*((an(17))*((T(l))^(-un(17))));
B18=(B18)*((an(18))*((T(l))^(-un(18))));
% B > > Segundo coeficiente virial
B=B1+B2+B3+B4+B5+B6+B7+B8+B9+B10+B11+B12+B13+B14+B15+B16+B17+B18;
% B > > Segundo coeficiente virial a temperatura base
Bb=Bb1+Bb2+Bb3+Bb4+Bb5+Bb6+Bb7+Bb8+Bb9+Bb10+Bb11+Bb12+Bb13+Bb14+Bb15+
Bb16+Bb17+Bb18;
% Clculo de la densidad molar d y la densidad
reducida D
d=(P(k))/(R*T(l));
D=((K)^(3))*(d);
% Clculo de la densidad molar d y la densidad
reducida D a
% condiciones base
db=(0.10155977)/(R*288.71);
Db=((K)^(3))*(db);
% Clculo del coeficiente Cn
for n=13:18
Cn(n)=(an(n))*((G+1-gn(n))^(gn(n)))*(((Q)^(2)+
1-qn(n))^(qn(n)))*((F+1-fn(n))^(fn(n)))*
((U)^(un(n)));
Z1=(Z1)+((Cn(n))*((T(l))^(-un(n))));
end
% Clculo del coeficiente Cn en condiciones base
for n=13:18
152
Cnb(n)=(an(n))*((G+1-gn(n))^(gn(n)))*(((Q)^(2)+
1-qn(n))^(qn(n)))*((F+1-fn(n))^(fn(n)))*
((U)^(un(n)));
Z1b=(Z1b)+((Cnb(n))*((288.71)^(-un(n))));
end
% Coeficiente en condiciones dadas
n=0; Cn=0;
for n=13:58
Cn(n)=(an(n))*((G+1-gn(n))^(gn(n)))*(((Q)^(2)+
1-qn(n))^(qn(n)))*((F+1-fn(n))^(fn(n)))*
((U)^(un(n)));
Z3(n)=((Cn(n))*((T(l))^(-un(n))));
Z4(n)=((bn(n))-((cn(n))*(kn(n))*((D)^(kn(n)))))
*((D)^(bn(n)));
Z5(n)=exp(-((cn(n))*((D)^(kn(n)))));
Z2=(Z2)+((Z3(n))*(Z4(n))*(Z5(n)));
end
% Paso anterior para condiciones base
n=0; Cnb=0;
for n=13:58
Cnb(n)=(an(n))*((G+1-gn(n))^(gn(n)))*(((Q)^(2)+
1-qn(n))^(qn(n)))*((F+1-fn(n))^(fn(n)))*((U)^(un(n)));
Z3b(n)=((Cnb(n))*((288.71)^(-un(n))));
Z4b(n)=((bn(n))-((cn(n))*(kn(n))*((D)^(kn(n)))))*((D)^(bn(n)));
Z5b(n)=exp(-((cn(n))*((Db)^(kn(n)))));
Z2b=(Z2b)+((Z3b(n))*(Z4b(n))*(Z5b(n)));
end
%Factor de compresibilidad
Z(k,l)=1+(d*B)-(D*Z1)+(Z2);
% Factor de compresibilidad en condiciones base
Zb(k,l)=1+(db*Bb)-(Db*Z1b)+(Z2b);
% Masa molar
n=0;
for n=1:N
Mr=(Mr)+((xi(n))*(Mrn(n)));
end
% Densidad msica
p=(P(k)*Mr)/(Z(k,l)*R*T(l));
% Densidad msica en condiciones base
pb=(0.10155977*Mr)/(Z(k,l)*R*288.71);
% Factor de supercompresibilidad
Fpv(k,l)=(Zb(k,l)/Z(k,l))^(1/2);
% Volumen corregido
Vc(k,l)=V*Z(k,l);
end
l=0;
end
disp(Elementos seleccionados:);
disp(CID);
disp(Fraccin de cada elemento:);
disp(xi);
disp(Factor de Compresibilidad Z:);
153
disp(Z);
disp(Factor de Supercompresibilidad Fpv:);
disp(Fpv);
disp(Volumen:);
disp(V);
disp(Volumen Corregido:);
disp(Vc);
end
154
G. TABLA 54E API MPMS 11.2.4
155
Captulo G TABLA 54E API MPMS 11.2.4
Figura G.1.: Prodecimiendo de implementacin para la Tabla 54E[7]
156
TABLA 54E API MPMS 11.2.4
Algoritmo G.1 Implementacin de la tabla 54E de API MPMS 11.2.4
% TABLA API 54 E
% Entradas: Densidad a 15C (Kg/m^3), Temperatura observada C
%Utilize EP (65/35) and EP (35/65)
%Den15=451.09;
%Tf=30.774;
% % 1, PASO redondeo de datos: Den15 redondeado a 0.1, Tf a 0.05
Den15=451.1; Tf=30.75;
% % 2, PASO Convertir temperatura a Kelvin
Tx=Tf+273.15
% % 3, PASO esta en el rango?
% % 4 PASO, convertir den15 a densidad relativa del agua a 60F
rTB=Den15/999.016 %
% 5 PASO, temperatura reducida y densidades de referencia
denx=rTB;
Tx=288.15;
% Constantes del fluido 1
r60ref1=0.429277;
Tcref1=333.67;
roc1=5.615;
Z1=0.28060;
k1=2.20970078464;
k2=-0.294253708172;
k3=-0.405754420098;
k4=0.319443433421;
%------------------------------------------------
Trx=Tx/Tcref1;
% debe ser menor a 1 para continuar
Tr60=519.67/(1.8*Tcref1);
tao=1-Trx; tao60=1-Tr60;
rosatxref=roc1*(1+(k1*tao^0.35+k3*tao^2+k4*tao^3)/(1+k2*tao^0.65));
rosat60ref=roc1*(1+(k1*tao60^0.35+k3*tao60^2+k4*tao60^3)/(1+k2*tao60^0.65));
% Densidad relativa del fluido 1 a 60F
rxref1=r60ref1*(rosatxref/rosat60ref);
% Constantes del fluido 2
r60ref2=0.470381;
Tcref2=352.46;
roc2=5.110;
Z2=0.27930;
k21=2.25341981320;
k22=-0.266542138024;
k23=-0.372756711655;
k24=0.384734185665;
%------------------------------------------------
Trx=Tx/Tcref2;
% debe ser menor a 1 para continuar
Tr60=519.67/(1.8*Tcref2);
tao=1-Trx;
tao60=1-Tr60;
rosatxref=roc2*(1+(k21*tao^0.35+k23*tao^2+k24*tao^3)/(1+k22*tao^0.65));
rosat60ref=roc2*(1+(k21*tao60^0.35+k23*tao60^2+k24*tao60^3)/(1+k22*tao60^0.65));
% Densidad relativa del fluido 2 a 60
F rxref2=r60ref2*(rosatxref/rosat60ref);
%--------Eleccin de fluidos de referencia, limites
r60high=r60ref2;
rxhigh=rxref2;
r60low=r60ref1;
rxlow=rxref1;
%-----------23/7
Delta=(denx-rxlow)/(rxhigh-rxlow)
%Si ? es menor que 0.001 entonces dejar igual a 0.001; Si ? es mayor que 0.999
%Entonces dejar igual a 0.999
r60mid=r60low+Delta*(r60high-r60low);
%------------------------- calculo de CTL-----------------------------------------
Delta2=(r60mid-r60low)/(r60high-r60low);
Tc2=Tcref1+Delta2*(Tcref2-Tcref1);
Trx2=Tx/Tc2;
157
Captulo G TABLA 54E API MPMS 11.2.4
Tr602=519.67/(1.8*Tc2);
h2=(Z1*roc1)/(Z2*roc2);
tao2=1-Trx2;
tao602=1-Tr602;
%Densidad de saturacion fluid 1
rosatx21=roc1*(1+(k1*tao2^0.35+k3*tao2^2+k4*tao2^3)/(1+k2*tao2^0.65));
rosat6021=roc1*(1+(k1*tao602^0.35+k3*tao602^2+k4*tao602^3)/(1+k2*tao602^0.65));
%fluid2 rosatx22=roc2*(1+(k21*tao2^0.35+k23*tao2^2+k24*tao2^3)/(1+k22*tao2^0.65));
%interpolating Factor X
X=rosat6021/(1+Delta2*((rosat6021/(h2*rosat6022))-1));
CTL=rosatx21/(X*(1+Delta2*((rosatx21/(h2*rosatx22))-1)));
%--------------------------------------------------------------------
rxmid=CTL*r60mid;
%Chequear convergencia Si rx esta entre rx,low y rx,mid y la diferencia
%entre r60,low y r60,mid es menor
%que 0.00000001 (108).
% Si rx esta entre rx,high y rx,mid y la diferencia entre r60,high y
% r60,mid es menor %que 0.00000001 (108).
denx rxlow rxmid
r60low-r60mid
r60high-r60mid
%T23/9
alpha=r60high-r60low;
beta=rxhigh^2-rxlow^2;
phi=(rxhigh-rxlow)/(rxmid-rxlow);
A=(alpha-phi*(r60mid-r60low))/(beta-phi*(rxmid^2-rxlow^2));
B=(alpha-A*beta)/(rxhigh-rxlow);
C=r60low-B*rxlow-A*rxlow^2;
r60trial=A*denx^2+B*denx+C;
% Este valor r60,trial puede res ajustado a r60,low o r60,high. con unas ecuaciones.
% Volver a calcular otro CTL para el r60trial
%pasos desde T24/4 to 13 en el T24/4 hay hay que elejir fluidos
%------------------------- calculo de CTL Trial-----------------------------------
Delta2=(r60trial-r60low)/(r60high-r60low);
Trx2=Tx/Tc2;
Tr602=519.67/(1.8*Tc2);
h2=(Z1*roc1)/(Z2*roc2);
tao2=1-Trx2;
tao602=1-Tr602;
%Densidad de saturacion
%fluid 1
%fluid2
%interpolating Factor X
X=rosat6021/(1+Delta2*((rosat6021/(h2*rosat6022))-1));
CTLtrial=rosatx21/(X*(1+Delta2*((rosatx21/(h2*rosatx22))-1)));
%--------------------------------------------------------------------
rxtrial=CTLtrial*r60trial;
%Chequearla convergencia de 60F desidad relativa.
%El calculo converge si la diferencia absoluta entre rx,trial and rx
%es menor que 0.00000001 (108).
rxtrial-denx
r60=r60trial
% % 6 PASO, La densidad relativa debe estar dentro de 0.34995 y 0.68805
% % 7 PASO, CTL1---factor de correccin de temperatura, reduce la
% % temperatura observada a un volumen a 60F, Se utiliza T24/4-13
%T24/4 eleccin de fluidos, Para este caso dieron los mismos
%t24/5:
Delta3=(r60-r60low)/(r60high-r60low);
%T24/6
158
TABLA 54E API MPMS 11.2.4
%T24/7:
Tx=Tf+273.15 Trx3=Tx/Tc3;
%Si la temperatura reducida Tr,x es mayor que 1.0, el fluido es supecritico
% y no existe como un lquido
%T24/8
Tr603=519.67/(1.8*Tc3);
%T24/9
h3=(Z1*roc1)/(Z2*roc2);
%T24/10 desnidad de saturacin
tao3=1-Trx3; tao603=1-Tr603;
%fluid 1
%fluid2
%T24/11 interpolating Factor X
X3=rosat6031/(1+Delta3*((rosat6031/(h3*rosat6032))-1));
%T24/13
CTL1=rosatx31/(X3*(1+Delta3*((rosatx31/(h3*rosatx32))-1)))
% % 8 PASO, CTL2 ---FACTOR CORRECCIN DE TEMPERATURA DESDE 60F a % % temperatura base 15C
%T24/4 eleccin de fluidos, Para este caso dieron los mismos %t24/5:
Delta3=(r60-r60low)/(r60high-r60low);
%T24/6 Tc3=Tcref1+Delta3*(Tcref2-Tcref1);
%T24/7: Tbk=288.15 Trx3=Tbk/Tc3;
%Si la temperatura reducida Tr,x es mayor que 1.0, entonces
%el fluido tiene condiciones supercriticas y no puede existir como un liquido.
%T24/8
Tr603=519.67/(1.8*Tc3);
%T24/9
h3=(Z1*roc1)/(Z2*roc2);
%T24/10 desnidad de saturacin
tao3=1-Trx3;
tao603=1-Tr603;
%fluid 1
%fluid2
%T24/11 interpolating Factor X
X3=rosat6031/(1+Delta3*((rosat6031/(h3*rosat6032))-1));
%T24/13
CTL2=rosatx31/(X3*(1+Delta3*((rosatx31/(h3*rosatx32))-1)))
% % 9 PASO, CTL REDUCE EL VOLUMEN DE UNA TEMPERATURA OBSERVADA
%A CONDICION BASE 15C
CTL=CTL1/CTL2
% % 10 PASO, Ctl debe ser positivo o mayor a 0
% % 11 PASO, redondear Correction Factor CTL a 0.00001.
CTLrounded=roundn(CTL,-5)
159

Tesis PDF

Cargado por

Copyright:

Formatos disponibles

Tesis PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tesis PDF

Cargado por

Copyright:

Formatos disponibles

DISEO E IMPLEMENTACIN DE UN SISTEMA DE

FISCALIZACIN DE GAS PARA LA EMPRESA ATLAS

También podría gustarte