14.11.2014 Views

simulación en hysys de los procesos productivos de enap magallanes

simulación en hysys de los procesos productivos de enap magallanes

simulación en hysys de los procesos productivos de enap magallanes

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

UNIVERSIDAD DE MAGALLANES<br />

FACULTAD DE INGENIERÄA<br />

DEPARTAMENTO DE QUÄMICA<br />

ÅSIMULACIÇN EN HYSYS É DE LOS PROCESOS<br />

PRODUCTIVOS DE ENAP MAGALLANESÑ<br />

LUIS EMILIO VILLEGAS VIVAR<br />

2007


UNIVERSIDAD DE MAGALLANES<br />

FACULTAD DE INGENIERÄA<br />

DEPARTAMENTO DE QUÄMICA<br />

ÅSIMULACIÇN EN HYSYS É DE LOS PROCESOS<br />

PRODUCTIVOS DE ENAP MAGALLANESÑ<br />

Trabajo <strong>de</strong> titulaciÄn pres<strong>en</strong>tado<br />

<strong>en</strong> conformidad a <strong>los</strong> requisitos<br />

para obt<strong>en</strong>er el tÅtulo <strong>de</strong><br />

Ing<strong>en</strong>iero Civil QuÅmico<br />

Profesor GuÅa: Sr. Lor<strong>en</strong>zo Lazaneo Cerda<br />

Supervisor ENAP: Sr. Nelson M<strong>en</strong>Çn<strong>de</strong>z DÅaz<br />

LUIS EMILIO VILLEGAS VIVAR<br />

2007


i<br />

RESUMEN<br />

En el pres<strong>en</strong>te trabajo <strong>de</strong> titulaciÄn, d<strong>en</strong>ominado ÅSimulaciÄn <strong>en</strong> HYSYS Ç <strong>de</strong> <strong>los</strong><br />

<strong>procesos</strong> <strong>productivos</strong> <strong>de</strong> ENAP MagallanesÉ, se han caracterizado <strong>los</strong> <strong>procesos</strong> tÑpicos llevados a<br />

cabo <strong>en</strong> las instalaciones <strong>de</strong> ENAP Magallanes, mediante simulaciones computacionales<br />

realizadas <strong>en</strong> el software Asp<strong>en</strong> HYSYSÇ.<br />

El proceso <strong>de</strong> creaciÄn <strong>de</strong> estas simulaciones se <strong>en</strong>cu<strong>en</strong>tra explicado completam<strong>en</strong>te<br />

<strong>en</strong> el ÅManual para simular <strong>procesos</strong> <strong>productivos</strong> <strong>de</strong> ENAP Magallanes <strong>en</strong> HYSYSÉ, que se<br />

adjunta. Dicho manual constituye a<strong>de</strong>mÖs una guÑa gradual <strong>de</strong> <strong>en</strong>seÜanza, que mediante<br />

ejemp<strong>los</strong>, casos guiados y consejos prÖcticos, permite obt<strong>en</strong>er las habilida<strong>de</strong>s necesarias para<br />

recrear, <strong>en</strong> el software ya m<strong>en</strong>cionado, difer<strong>en</strong>tes situaciones relacionadas con <strong>los</strong> equipos<br />

tÑpicos que utiliza ENAP Magallanes <strong>en</strong> su proceso productivo.<br />

Se ha realizado, <strong>en</strong>tre otras, la simulaciÄn completa <strong>de</strong> Planta Cull<strong>en</strong> y Planta PosesiÄn,<br />

g<strong>en</strong>erando programas que permit<strong>en</strong> estimar el comportami<strong>en</strong>to <strong>de</strong> las variables <strong>de</strong> proceso<br />

fr<strong>en</strong>te a cambios <strong>en</strong> las condiciones <strong>de</strong> operaciÄn <strong>de</strong> las plantas.<br />

A<strong>de</strong>mÖs se han comparado <strong>los</strong> valores normales <strong>de</strong> operaciÄn <strong>de</strong> cada planta con el<br />

resultado que <strong>en</strong>trega HYSYS, y se ha concluido que las similitu<strong>de</strong>s son a<strong>de</strong>cuadas d<strong>en</strong>tro <strong>de</strong>l<br />

rango <strong>de</strong> precisiÄn requerido, y por lo tanto es importante que el personal <strong>de</strong> ENAP Magallanes<br />

pueda t<strong>en</strong>er acceso a las simulaciones realizadas, y conozca la forma <strong>de</strong> <strong>de</strong>sarrollar ástas para<br />

realizar <strong>los</strong> cambios <strong>en</strong> valores <strong>de</strong> variables o estructura <strong>de</strong>l mo<strong>de</strong>lo. Esto se lograrÖ mediante la<br />

difusiÄn <strong>de</strong>l ÅManual para simular <strong>procesos</strong> <strong>productivos</strong> <strong>de</strong> ENAP Magallanes <strong>en</strong> HYSYSÉ <strong>en</strong> las<br />

difer<strong>en</strong>tes instalaciones <strong>de</strong> la empresa.<br />

Debido a la imposibilidad <strong>de</strong> activar el modo dinÖmico <strong>en</strong> HYSYS, por limitaciones <strong>de</strong> la<br />

lic<strong>en</strong>cia usada <strong>en</strong> la empresa, las simulaciones se han limitado exclusivam<strong>en</strong>te a casos<br />

estacionarios.


ii<br />

ÄNDICE<br />

RESUMEN................................................................................................................................................ I<br />

ÄNDICE.................................................................................................................................................... II<br />

ILUSTRACIONES Y TABLAS .....................................................................................................................IV<br />

INTRODUCCIÅN ....................................................................................................................................VI<br />

I - ANTECEDENTES GENERALES............................................................................................................... 1<br />

1.1 ENAP MAGALLANES........................................................................................................................1<br />

1.2 PLANTEAMIENTO DEL PROBLEMA..................................................................................................3<br />

1.3 OBJETIVOS ......................................................................................................................................5<br />

1.3.1 OBJETIVO GENERAL ..................................................................................................................5<br />

1.3.2 OBJETIVOS ESPECÄFICOS............................................................................................................5<br />

II - ANTECEDENTES TEÅRICOS ................................................................................................................ 7<br />

2.1 USO DE SIMULADORES EN LA INDUSTRIA QUàMICA.......................................................................7<br />

2.2 ETAPAS GENERALES DE UNA SIMULACIâN...................................................................................10<br />

2.3 EL SIMULADOR HYSYS Ç .................................................................................................................12<br />

2.4 CONSIDERACIONES TEâRICAS UTILIZADAS POR HYSYS................................................................15<br />

2.4.1 COMPRESORES RECÄPROCOS ..................................................................................................16<br />

2.4.2 INTERCAMBIADORES DE CALOR DE TUBO Y CARCASA...........................................................19<br />

2.4.3 REACCIONES QUÄMICAS Y REACTORES EN HYSYS...................................................................24<br />

2.4.4 COLUMNAS DE DESTILACIÅN..................................................................................................27<br />

2.4.5 AJUSTADORES .........................................................................................................................33<br />

2.4.6 RECICLOS.................................................................................................................................35<br />

III - PARTE EXPERIMENTAL Y RESULTADOS........................................................................................... 38<br />

3.1 ESQUEMA DE DESARROLLO..........................................................................................................38<br />

3.2 SIMULACIONES REALIZADAS PARA ENAP MAGALLANES..............................................................41<br />

3.3 CONSIDERACIONES ESPECàFICAS DE CADA SIMULACIâN.............................................................43<br />

3.3.1 COMPRESOR HRA-5 DE PLANTA CULLEN................................................................................43<br />

3.3.2 COMPRESOR HRA-2 DE ESTACIÅN COMPRESORA CALAFATE ................................................44


iii<br />

3.3.3 COMPRESOR TLA-3 DE PLANTA POSESIÅN.............................................................................44<br />

3.3.4 COMPRESOR TC-10 DE PLANTA CULLEN.................................................................................46<br />

3.3.5 TURBINA DEL TREN SULZER ....................................................................................................46<br />

3.3.6 ESQUEMA DE PRODUCCIÅN DE PLANTAS POSESIÅN Y CABO NEGRO ...................................48<br />

3.3.7 ELEMPLO DE CARACTERIZACIÅN DE PETRÅLEOS EN HYSYS...................................................49<br />

3.3.8 EJEMPLO DE LOOP DE GASODUCTOS .....................................................................................49<br />

3.3.9 PLANTA CULLEN......................................................................................................................49<br />

3.3.10 PLANTA POSESIÅN ................................................................................................................52<br />

3.4 DATOS Y RESULTADOS..................................................................................................................56<br />

3.4.1 COMPRESOR HRA-5 DE PLANTA CULLEN................................................................................58<br />

3.4.2 COMPRESOR HRA-5 DE ESTACIÅN COMPRESORA CALAFATE ................................................59<br />

3.4.3 COMPRESOR TLA-3 DE PLANTA POSESIÅN.............................................................................60<br />

3.4.4 COMPRESOR TC-10 DE PLANTA CULLEN.................................................................................61<br />

3.4.5 TURBINA DEL TREN SULZER ....................................................................................................62<br />

3.4.6 PLANTA CULLEN......................................................................................................................63<br />

3.4.7 PLANTA POSESIÅN ..................................................................................................................64<br />

3.5 VALIDACIâN DE LAS SIMULACIONES OBTENIDAS ........................................................................66<br />

3.5.1 VALIDACIÅN DE PLANTA CULLEN ...........................................................................................67<br />

3.5.2 VALIDACIÅN DE PLANTA POSESIÅN........................................................................................71<br />

3.6 ANäLISIS Y DISCUSIâN DE RESULTADOS.......................................................................................76<br />

IV - CONCLUSIONES.............................................................................................................................. 80<br />

V - BIBLIOGRAFÄA ................................................................................................................................. 83<br />

5.1 FUENTES CONSULTADAS ..............................................................................................................83<br />

5.2 INTERNET ã SITIOS CONSULTADOS...............................................................................................83<br />

ANEXOS ............................................................................................................................................... 84<br />

ASPECTO GRÇFICO DE LAS SIMULACIONES........................................................................................... 85<br />

INFORMACIÅN DE COMPRESORES ....................................................................................................... 93


iv<br />

ILUSTRACIONES Y TABLAS<br />

ILUSTRACIâN 1.- ASPECTO GENERAL DE PLANTA POSESIâN...........................................................................................................2<br />

ILUSTRACIâN 2.- ASPECTO GENERAL DE PLANTA CULLEN..............................................................................................................2<br />

ILUSTRACIâN 3.- SIMULATION BASIS MANAGER DE HYSYS ........................................................................................................13<br />

ILUSTRACIâN 4.- VENTANA DE PROPIEDADES PARA UN INTERCAMBIADOR DE TUBO Y CARCASA ..........................................................13<br />

ILUSTRACIâN 5.- ASPECTO DE UNA SIMULACIâN TàPICA EN HYSYS...............................................................................................14<br />

ILUSTRACIâN 6.- ESQUEMA DEL CILINDRO DE UN COMPRESOR RECàPROCO.....................................................................................16<br />

ILUSTRACIâN 7.- ESPECIFICACIâN EN MODO WEIGHTED.............................................................................................................22<br />

ILUSTRACIâN 8.- INGRESO DE PARäMETROS FàSICOS DEL INTERCAMBIADOR...................................................................................22<br />

ILUSTRACIâN 9.- ALGUNAS REACCIONES PREDEFINIDAS EN HYSYS...............................................................................................24<br />

ILUSTRACIâN 10.- TIPOS DE REACTORES GENERALES..................................................................................................................26<br />

ILUSTRACIâN 11.- ASPECTO DE UNA COLUMNA DE DESTILACIâN EN HYSYS...................................................................................28<br />

ILUSTRACIâN 12.- ESQUEMA DE UN PROCESO DE SEPARACIâN POR ETAPAS DE EQUILIBRIO ...............................................................29<br />

ILUSTRACIâN 13.- ESQUEMA TEâRICO DE UNA COLUMNA DE DESTILACIâN ....................................................................................31<br />

ILUSTRACIâN 14.- ASPECTO DEL AJUSTADOR EN HYSYS.............................................................................................................35<br />

ILUSTRACIâN 15.- ASPECTO DEL MâDULO DE RECICLO EN HYSYS................................................................................................36<br />

ILUSTRACIâN 16.- MODELO INICIAL PARA EL COMPRESOR DE DOS ETAPAS .....................................................................................45<br />

ILUSTRACIâN 17.- MODELO AJUSTADO PARA CUMPLIR EL BALANCE DE MASA.................................................................................45<br />

ILUSTRACIâN 18.- ESQUEMA DE LA UNIDAD REGENERADORA DE GLICOL DE PLANTA CULLEN ...........................................................51<br />

ILUSTRACIâN 19.- ASPECTO DE LA UNIDAD REGENERADORA DE GLICOL SIMULADA EN HYSYS..........................................................52<br />

ILUSTRACIâN 20.- ESQUEMA DE UN REHERVIDOR TIPO TERMOSIFâN ............................................................................................53<br />

TABLA 1.- MåTODOS DE RESOLUCIâN DE COLUMNAS DE DESTILACIâN EN HYSYS ...........................................................................33<br />

TABLA 2.- COMPOSICIONES INGRESADAS PARA LAS SIMULACIONES DE COMPRESORES......................................................................57<br />

TABLA 3.- COMPOSICIONES INGRESADAS A LAS SIMULACIONES DE TREN SULZER, CULLEN Y POSESIâN ..............................................57<br />

TABLA 4.- DEFINICIONES DE CORRIENTES EN LA SIMULACIâN DEL COMPRESOR HRA-5.....................................................................58<br />

TABLA 5.- DATOS DEL COMPRESOR HRA-5 .............................................................................................................................58<br />

TABLA 6.- RESULTADOS OBTENIDOS HRA-5 ............................................................................................................................58<br />

TABLA 7.- DEFINICIONES DE CORRIENTES EN LA SIMULACIâN DEL COMPRESOR HRA-2.....................................................................59<br />

TABLA 8.- DATOS DEL COMPRESOR HRA-2 DIVIDIDO EN DOS COMPRESORES TEâRICOS ...................................................................59<br />

TABLA 9.- RESULTADOS OBTENIDOS PARA EL COMPRESOR HRA-2 ...............................................................................................59


v<br />

TABLA 10.- DEFINICIONES DE CORRIENTES EN LA SIMULACIâN DEL COMPRESOR TLA-3....................................................................60<br />

TABLA 11.- DATOS DEL COMPRESOR TLA-3 DIVIDIDO EN DOS ETAPAS. .........................................................................................60<br />

TABLA 12.- RESULTADOS OBTENIDOS TLA-3 PARA PRESIâN INTERETAPA DE 1,100 PSIG..................................................................60<br />

TABLA 13.- ESPECIFICACIONES REALIZADAS AL COMPRESOR TC-10 ..............................................................................................61<br />

TABLA 14.- RESULTADOS OBTENIDOS AL SIMULAR EL COMPRESOR TC-10 .....................................................................................61<br />

TABLA 15.- DATOS INGRESADOS A LA SIMULACIâN DE LA TURBINA DEL TREN SULZER.....................................................................62<br />

TABLA 16.- RESULTADOS OBTENIDOS PARA EL TREN SULZER .....................................................................................................62<br />

TABLA 17.- COMPOSICIONES DEL GAS DE ENTRADA A TURBINA MEDIANTE LOS DOS MåTODOS...........................................................62<br />

TABLA 18.- DATOS INGRESADOS PARA LA SIMULACIâN DE PLANTA CULLEN ...................................................................................63<br />

TABLA 19.- CROMATOGRAFàA DEL GAS RESIDUAL OBTENIDA EN LA SIMULACIâN .............................................................................63<br />

TABLA 20.- RESULTADOS DE LA SIMULACIâN DE PLANTA CULLEN.................................................................................................63<br />

TABLA 21.- LISTADO DE ESPECIFICACIONES UTILIZADAS PARA SIMULAR PLANTA POSESIâN................................................................64<br />

TABLA 22.-TEMPERATURAS Y PRESIONES OBTENIDAS DE LA SIMULACIâN DE PLANTA POSESIâN.........................................................64<br />

TABLA 23.- RESULTADOS OBTENIDOS AL SIMULAR LA COLUMNA V-5............................................................................................65<br />

TABLA 24.- COMPOSICIâN DE LAS CORRIENTES OBTENIDAS EN LA SIMULACIâN...............................................................................65<br />

TABLA 25.- VALORES NORMALES DE OPERACIâN PLANTA CULLEN Y PLANTA POSESIâN....................................................................66<br />

TABLA 26.- COMPOSICIONES DEL GAS DE ENTRADA PARA LAS SITUACIONES ANALIZADAS ..................................................................68<br />

TABLA 27.- SIMULACIâN DE PLANTA CULLEN DEL DàA 24/07/07................................................................................................68<br />

TABLA 28.- SIMULACIâN DE PLANTA CULLEN DEL DàA 25/07/07................................................................................................69<br />

TABLA 29.- SIMULACIâN DE PLANTA CULLEN DEL DàA 26/07/07................................................................................................69<br />

TABLA 30.- SIMULACIâN DE PLANTA CULLEN DEL DàA 27/07/07................................................................................................70<br />

TABLA 31.- SIMULACIâN DE PLANTA CULLEN DEL DàA 28/07/07................................................................................................70<br />

TABLA 32.- VALORES ESTABLES DE OPERACIâN DE PLANTA POSESIâN...........................................................................................71<br />

TABLA 33.- INFORMACIâN DE LABORATORIO PARA LA SIMULACIâN DE PLANTA POSESIâN................................................................71<br />

TABLA 34.- SIMULACIâN DE PLANTA POSESIâN DàA 31/07/07 ..................................................................................................72<br />

TABLA 35.- CROMATOGRAFàA OBTENIDA PARA LA SIMULACIâN DE PLANTA POSESIâN DàA 31/07/07................................................72<br />

TABLA 36.- SIMULACIâN DE PLANTA POSESIâN DàA 01/08/07 ..................................................................................................73<br />

TABLA 37.- CROMATOGRAFàA OBTENIDA PARA LA SIMULACIâN DE PLANTA POSESIâN DàA 01/08/07................................................73<br />

TABLA 38.- SIMULACIâN DE PLANTA POSESIâN DàA 04/08/07 ..................................................................................................74<br />

TABLA 39.- CROMATOGRAFàA OBTENIDA PARA LA SIMULACIâN DE PLANTA POSESIâN DàA 04/08/07................................................74<br />

TABLA 40.- SIMULACIâN DE PLANTA POSESIâN DàA 05/08/07 ..................................................................................................75<br />

TABLA 41.- CROMATOGRAFàA OBTENIDA PARA LA SIMULACIâN DE PLANTA POSESIâN DàA 05/08/07................................................75


vi<br />

INTRODUCCIÇN<br />

Las simulaciones computacionales <strong>de</strong> <strong>procesos</strong> <strong>productivos</strong> <strong>en</strong> la industria quÑmica<br />

constituy<strong>en</strong> una herrami<strong>en</strong>ta muy çtil para el <strong>de</strong>sarrollo <strong>de</strong> <strong>los</strong> <strong>procesos</strong> <strong>de</strong> ENAP Magallanes,<br />

pues aprovechan la rapi<strong>de</strong>z <strong>de</strong> cÖlculo <strong>de</strong> un ord<strong>en</strong>ador para <strong>en</strong>tregar resultados <strong>de</strong> forma<br />

instantÖnea. En una industria <strong>en</strong> constante cambio, como lo es la industria <strong>de</strong> <strong>los</strong> combustibles,<br />

la habilidad para pre<strong>de</strong>cir rÖpidam<strong>en</strong>te quá pasarÖ si se modifican una o mÖs variables <strong>de</strong><br />

proceso es fundam<strong>en</strong>tal a la hora <strong>de</strong> controlar a<strong>de</strong>cuadam<strong>en</strong>te <strong>los</strong> <strong>procesos</strong>.<br />

Para lograr esto, ENAP Magallanes posee una lic<strong>en</strong>cia para utilizar el software Asp<strong>en</strong><br />

HYSYS.Plant Ç (comçnm<strong>en</strong>te llamado HYSYS). Dicho software permite recrear <strong>los</strong> equipos e<br />

instalaciones que llevan a cabo <strong>los</strong> <strong>procesos</strong> <strong>productivos</strong>, asignar condiciones <strong>de</strong> trabajo y<br />

obt<strong>en</strong>er resultados estimativos <strong>de</strong> variables <strong>de</strong> operaciÄn.<br />

El trabajo que a continuaciÄn se pres<strong>en</strong>ta <strong>en</strong>trega <strong>los</strong> resultados obt<strong>en</strong>idos al realizar<br />

simulaciones, <strong>en</strong> el software ya m<strong>en</strong>cionado, <strong>de</strong> <strong>los</strong> <strong>procesos</strong> relacionados con la producciÄn <strong>de</strong><br />

gas natural y petrÄleo que <strong>de</strong>sarrolla ENAP Magallanes <strong>en</strong> la actualidad.<br />

La investigaciÄn se ha realizado dividi<strong>en</strong>do el proceso productivo <strong>de</strong> ENAP Magallanes<br />

<strong>en</strong> unida<strong>de</strong>s bÖsicas repres<strong>en</strong>tativas <strong>de</strong> las operaciones diarias (compresores, expansores,<br />

gasoductos, separadores, intercambiadores <strong>de</strong> calor, etc.), se han <strong>de</strong>sarrollado aplicaciones<br />

simples para ejemplificar el uso <strong>de</strong> estas unida<strong>de</strong>s bÖsicas y a continuaciÄn se han integrado <strong>en</strong><br />

simulaciones mÖs complejas que son capaces <strong>de</strong> repres<strong>en</strong>tar sectores mÖs amplios <strong>de</strong> la red <strong>de</strong><br />

sistema productivo.<br />

Desarrollar simulaciones computacionales <strong>de</strong> <strong>procesos</strong> <strong>de</strong> petrÄleo o gas implica<br />

conocer a fondo el sistema que se <strong>de</strong>sea simular, tanto <strong>en</strong> el rango permitido <strong>de</strong> operaciÄn <strong>de</strong><br />

ciertos equipos como <strong>en</strong> valores tÑpicos <strong>de</strong> variables <strong>de</strong> proceso (presiÄn, temperatura, flujos o<br />

composiciones, etc.).


vii<br />

Una vez realizada la simulaciÄn, que se basa <strong>en</strong> el flowsheet <strong>de</strong> planta, se <strong>de</strong>be<br />

comparar con <strong>los</strong> resultados medidos <strong>en</strong> terr<strong>en</strong>o, y para esto se ha utilizado diversa<br />

informaciÄn; informes <strong>de</strong> laboratorio, hojas <strong>de</strong> diseÜo <strong>en</strong>tregadas por <strong>los</strong> fabricantes, informes<br />

<strong>de</strong> mant<strong>en</strong>imi<strong>en</strong>to predictivo y, principalm<strong>en</strong>te, mediante mediciÄn y observaciÄn directa <strong>de</strong> las<br />

variables o <strong>procesos</strong> <strong>de</strong> interás.<br />

Se han reunido todas las simulaciones <strong>en</strong> el ÅManual para simular <strong>procesos</strong> <strong>productivos</strong><br />

<strong>de</strong> ENAP Magallanes <strong>en</strong> HYSYSÉ, docum<strong>en</strong>to que sirve como manual <strong>de</strong> uso <strong>de</strong>l software, y que<br />

<strong>de</strong>sarrolla <strong>de</strong>s<strong>de</strong> <strong>los</strong> aspectos mÖs bÖsicos necesarios para <strong>en</strong>t<strong>en</strong><strong>de</strong>r las simulaciones, hasta<br />

situaciones <strong>de</strong> carÖcter avanzado que requier<strong>en</strong> un conocimi<strong>en</strong>to global <strong>de</strong> lo <strong>de</strong>tallado <strong>en</strong> el<br />

manual.<br />

Este manual servirÖ como instrum<strong>en</strong>to <strong>de</strong> <strong>en</strong>seÜanza y apoyo para <strong>los</strong> ing<strong>en</strong>ieros <strong>de</strong><br />

<strong>procesos</strong> <strong>de</strong> ENAP Magallanes y operadores <strong>de</strong> terr<strong>en</strong>o, pues <strong>de</strong> una forma metÄdica y <strong>de</strong>tallada<br />

<strong>en</strong>trega todos <strong>los</strong> conocimi<strong>en</strong>tos necesarios para po<strong>de</strong>r simular las situaciones tÑpicas <strong>de</strong> <strong>los</strong><br />

<strong>procesos</strong> <strong>productivos</strong> llevados a cabo <strong>en</strong> ENAP Magallanes. Todas las simulaciones analizadas <strong>en</strong><br />

este trabajo se <strong>en</strong>cu<strong>en</strong>tran <strong>en</strong> un disco compacto, que se adjunta al manual. A<strong>de</strong>mÖs, mediante<br />

el <strong>de</strong>sarrollo <strong>de</strong> <strong>los</strong> ejemp<strong>los</strong> y ejercicios <strong>de</strong>tallados <strong>en</strong> el manual, una vez finalizado el<br />

<strong>de</strong>sarrollo se lograrÖ un conocimi<strong>en</strong>to completo <strong>de</strong> las opciones que HYSYS ofrece a <strong>los</strong><br />

profesionales vinculados con <strong>los</strong> <strong>procesos</strong> <strong>de</strong> ENAP Magallanes.


CAPÄTULO I<br />

ANTECEDENTES GENERALES


1<br />

I - ANTECEDENTES GENERALES<br />

1.1 ENAP MAGALLANES<br />

La Empresa Nacional <strong>de</strong>l PetrÄleo ENAP, mediante sus instalaciones <strong>en</strong> la regiÄn <strong>de</strong><br />

Magallanes, se <strong>de</strong>sarrolla <strong>en</strong> el negocio <strong>de</strong> explotaciÄn <strong>de</strong> combustibles fÄsiles.<br />

Con sus plantas <strong>de</strong> PosesiÄn, Cabo Negro, Gregorio y Cull<strong>en</strong>, <strong>de</strong>sarrolla mçltiples<br />

<strong>procesos</strong> <strong>productivos</strong>, como por ejemplo:<br />

Ä Transporte <strong>de</strong> gas <strong>de</strong> productores ubicados <strong>en</strong> Arg<strong>en</strong>tina, para abastecer la planta<br />

productora <strong>de</strong> metanol Methanex Chile Ltd.<br />

Ä SeparaciÄn y comercializaciÄn <strong>de</strong>l propano, butano y gasolinas a partir <strong>de</strong>l gas<br />

obt<strong>en</strong>ido <strong>de</strong> yacimi<strong>en</strong>tos ubicados <strong>en</strong> Tierra <strong>de</strong>l Fuego y <strong>en</strong> contin<strong>en</strong>te.<br />

Ä ProducciÄn <strong>de</strong> combustibles mediante refinaciÄn <strong>de</strong>l petrÄleo prov<strong>en</strong>i<strong>en</strong>te <strong>de</strong><br />

estos mismos yacimi<strong>en</strong>tos y <strong>de</strong> contratos <strong>de</strong> compra con difer<strong>en</strong>tes<br />

productores.<br />

ENAP es la principal empresa <strong>de</strong> este Ömbito <strong>en</strong> la regiÄn. Los <strong>procesos</strong> llevados a cabo<br />

por ENAP <strong>en</strong>globan una amplia red <strong>de</strong> gasoductos, poliductos, estaciones compresoras, plantas<br />

<strong>de</strong>shidratadoras, pozos <strong>de</strong> producciÄn, etc.<br />

Para lo que a este trabajo se refiere, se consi<strong>de</strong>raron dos <strong>de</strong> sus plantas:<br />

Ä Planta PosesiÄn: Ubicada <strong>en</strong> el contin<strong>en</strong>te, a aproximadam<strong>en</strong>te 220 km. De Punta<br />

Ar<strong>en</strong>as, realiza separaciÄn <strong>de</strong> propano y mÖs pesados a partir <strong>de</strong>l gas que<br />

ingresa a planta. G<strong>en</strong>era por lo tanto dos productos, Raw Product y Gas<br />

Residual.


2<br />

Ä Planta Cull<strong>en</strong>: En la actualidad Planta Cull<strong>en</strong> sirve como planta <strong>de</strong>shidratadora <strong>de</strong>l<br />

gas que provi<strong>en</strong>e <strong>de</strong> yacimi<strong>en</strong>tos <strong>de</strong> la isla <strong>de</strong> Tierra <strong>de</strong>l Fuego. Separa a<strong>de</strong>mÖs<br />

<strong>los</strong> cond<strong>en</strong>sables para obt<strong>en</strong>er una corri<strong>en</strong>te <strong>de</strong> Gasolinas y un Gas Residual.<br />

Este çltimo ti<strong>en</strong>e varias finalida<strong>de</strong>s, gas <strong>de</strong> reinyecciÄn, suministro <strong>de</strong> gas<br />

combustible para <strong>los</strong> equipos asociados, suministro <strong>de</strong> gas para las plantas <strong>de</strong><br />

Methanex Chile Ltd., etc.<br />

IlustraciÄn 1.- Aspecto g<strong>en</strong>eral <strong>de</strong> Planta PosesiÄn<br />

IlustraciÄn 2.- Aspecto g<strong>en</strong>eral <strong>de</strong> Planta Cull<strong>en</strong>


3<br />

1.2 PLANTEAMIENTO DEL PROBLEMA<br />

ENAP Magallanes estÖ interesado <strong>en</strong> <strong>de</strong>sarrollar un proyecto que permita capturar, a<br />

travás <strong>de</strong> un trabajo <strong>de</strong> docum<strong>en</strong>taciÄn y anÖlisis, el conocimi<strong>en</strong>to adquirido por sus ing<strong>en</strong>ieros<br />

al simular <strong>en</strong> el software Asp<strong>en</strong> HYSYS <strong>los</strong> <strong>procesos</strong> relacionados con la producciÄn <strong>de</strong> gas<br />

natural y petrÄleo. La utilizaciÄn <strong>de</strong> este simulador por parte <strong>de</strong> <strong>los</strong> ing<strong>en</strong>ieros <strong>de</strong> proceso ha<br />

resultado satisfactoria y precisa, para <strong>de</strong>terminar condiciones <strong>de</strong> operaciÄn, posibles mejoras <strong>en</strong><br />

el sistema investigado y respuestas <strong>de</strong> operaciÄn fr<strong>en</strong>te a cambios <strong>en</strong> la materia prima.<br />

Se <strong>de</strong>sea <strong>de</strong>sarrollar simulaciones <strong>en</strong> HYSYS, para <strong>los</strong> casos mÖs tÑpicos <strong>de</strong> <strong>procesos</strong><br />

reales que actualm<strong>en</strong>te estÖn <strong>en</strong> operaciÄn <strong>en</strong> el sistema <strong>de</strong> producciÄn <strong>de</strong> ENAP Magallanes.<br />

El trabajo se realizarÖ dividi<strong>en</strong>do el proceso productivo <strong>de</strong> ENAP Magallanes <strong>en</strong><br />

unida<strong>de</strong>s bÖsicas repres<strong>en</strong>tativas <strong>de</strong> las operaciones diarias (compresores, expansores,<br />

gasoductos, oleoductos, unida<strong>de</strong>s <strong>de</strong>shidratadoras, etc) para posteriorm<strong>en</strong>te integrarlas <strong>en</strong><br />

simulaciones mÖs complejas que sean capaces <strong>de</strong> repres<strong>en</strong>tar sectores mÖs amplios <strong>de</strong> la red <strong>de</strong><br />

sistema productivo.<br />

El objetivo <strong>de</strong> una simulaciÄn no es siempre recrear la planta <strong>en</strong> su totalidad, para<br />

ENAP es mÖs çtil po<strong>de</strong>r analizar casos pequeÜos y especÑficos (un compresor, una turbina, etc.),<br />

pero para dar mayor alcance al trabajo, se simularÖn las plantas <strong>de</strong> PosesiÄn y Cull<strong>en</strong> <strong>en</strong> su<br />

totalidad.<br />

Basado <strong>en</strong> estas premisas, se <strong>de</strong>be <strong>de</strong>sarrollar el ÅManual para simular <strong>procesos</strong><br />

<strong>productivos</strong> <strong>de</strong> ENAP Magallanes <strong>en</strong> HYSYSÉ, que explique paso a paso cÄmo simular <strong>los</strong><br />

<strong>procesos</strong> <strong>de</strong> ENAP Magallanes para, llegado el caso, a<strong>de</strong>lantar una posible soluciÄn fr<strong>en</strong>te a un<br />

problema <strong>de</strong>terminado.


4<br />

El usuario que <strong>de</strong>sarrolle por completo el manual, habrÖ adquirido <strong>los</strong> conocimi<strong>en</strong>tos<br />

necesarios para crear simulaciones <strong>de</strong> cualquier situaciÄn que se pres<strong>en</strong>ta <strong>en</strong> ENAP; re<strong>de</strong>s <strong>de</strong><br />

gasoductos, tr<strong>en</strong> <strong>de</strong> intercambio <strong>de</strong> calor, torres <strong>de</strong> <strong>de</strong>stilaciÄn, <strong>procesos</strong> <strong>de</strong> expansiÄn y<br />

compresiÄn, separaciÄn lÑquido-vapor, etc. Todos estos casos estÖn explicados paso a paso <strong>en</strong> el<br />

manual que se adjunta como anexo a este trabajo.<br />

La exactitud <strong>de</strong> las simulaciones <strong>de</strong>sarrolladas serÖ cotejada con datos obt<strong>en</strong>idos<br />

directam<strong>en</strong>te <strong>en</strong> planta, informes <strong>de</strong> control <strong>de</strong> calidad <strong>de</strong> <strong>los</strong> laboratorios <strong>de</strong> ENAP Magallanes<br />

e Informes <strong>de</strong> Mant<strong>en</strong>imi<strong>en</strong>to Predictivo asociado a <strong>los</strong> difer<strong>en</strong>tes equipos simulados.


5<br />

1.3 OBJETIVOS<br />

1.3.1 OBJETIVO GENERAL<br />

Implem<strong>en</strong>tar simulaciones, <strong>de</strong>sarrolladas <strong>en</strong> el simulador Asp<strong>en</strong> HYSYS Ç , <strong>de</strong> <strong>los</strong><br />

<strong>procesos</strong> diarios tÑpicos <strong>de</strong> ENAP Magallanes, <strong>de</strong>bidam<strong>en</strong>te calibrados con <strong>los</strong> datos <strong>de</strong> terr<strong>en</strong>o.<br />

1.3.2 OBJETIVOS ESPECÄFICOS<br />

Ä Simular <strong>en</strong> Asp<strong>en</strong> HYSYS Ç las unida<strong>de</strong>s bÖsicas <strong>de</strong> <strong>los</strong> <strong>procesos</strong> <strong>productivos</strong> <strong>de</strong><br />

ENAP Magallanes.<br />

Ä Realizar s<strong>en</strong>sibilizaciÄn <strong>de</strong> las simulaciones a difer<strong>en</strong>tes condiciones <strong>de</strong> operaciÄn<br />

contrastando <strong>los</strong> resultados <strong>de</strong> Asp<strong>en</strong> HYSYSÇ con <strong>los</strong> datos reales <strong>de</strong> operaciÄn.<br />

Ä Recopilar la informaciÄn y experi<strong>en</strong>cia adquirida para crear un ÅManual para<br />

simular <strong>procesos</strong> <strong>productivos</strong> <strong>de</strong> ENAP Magallanes <strong>en</strong> HYSYSÉ, <strong>de</strong>stinado a <strong>los</strong><br />

ing<strong>en</strong>ieros <strong>de</strong> <strong>procesos</strong> y operadores <strong>de</strong> ENAP Magallanes y a <strong>los</strong> alumnos <strong>de</strong><br />

Ing<strong>en</strong>ierÑa QuÑmica <strong>de</strong> la Universidad <strong>de</strong> Magallanes.


CAPÄTULO II<br />

ANTECEDENTES TEÇRICOS


7<br />

II - ANTECEDENTES TEÇRICOS<br />

2.1 USO DE SIMULADORES EN LA INDUSTRIA QUÄMICA<br />

Todos <strong>los</strong> profesionales relacionados al <strong>de</strong>sarrollo diario <strong>de</strong> una planta quÑmica ti<strong>en</strong><strong>en</strong><br />

la responsabilidad <strong>de</strong> <strong>de</strong>cidir quá hacer fr<strong>en</strong>te a difer<strong>en</strong>tes situaciones <strong>de</strong> operaciÄn. Si el flujo<br />

<strong>de</strong> <strong>en</strong>trada a planta aum<strong>en</strong>ta, ése <strong>de</strong>be aÜadir mÖs glicol <strong>de</strong> <strong>de</strong>shidrataciÄn? éCuÖnto? Si la<br />

temperatura <strong>de</strong>l separador ha bajado, équá vÖlvula se <strong>de</strong>be abrir? Casos como estos, y quizÖ<br />

mÖs complicados, se pres<strong>en</strong>tan a diario, y es importante que la persona a cargo sepa quá acciÄn<br />

resulta ser la mÖs a<strong>de</strong>cuada para solucionar la problemÖtica pres<strong>en</strong>te.<br />

Lam<strong>en</strong>tablem<strong>en</strong>te las condiciones <strong>de</strong> trabajo <strong>en</strong> la industria no permit<strong>en</strong> el proceso <strong>de</strong><br />

Åprueba y errorÉ, pues implem<strong>en</strong>tar una <strong>de</strong>cisiÄn equivocada pue<strong>de</strong> llegar a t<strong>en</strong>er<br />

consecu<strong>en</strong>cias in<strong>de</strong>seadas <strong>en</strong> las condiciones <strong>de</strong> especificaciÄn <strong>de</strong>l producto.<br />

Para obt<strong>en</strong>er una refer<strong>en</strong>cia <strong>de</strong>l efecto que produce cierta acciÄn sobre un sistema, se<br />

recurre a un mo<strong>de</strong>lo. Un mo<strong>de</strong>lo es una repres<strong>en</strong>taciÄn simplificada <strong>de</strong>l sistema estudiado,<br />

basado <strong>en</strong> <strong>los</strong> esquemas teÄricos subyac<strong>en</strong>tes <strong>en</strong> el sistema real. Al proceso <strong>de</strong> experim<strong>en</strong>tar<br />

con el mo<strong>de</strong>lo, esto es, cambiar parÖmetros, modificar esquemas, etc., se le d<strong>en</strong>omina simular.<br />

Un dispositivo que permita simular se d<strong>en</strong>omina simulador.<br />

Un simulador no es un recurso exclusivo <strong>de</strong> la industria quÑmica. Por ejemplo, para la<br />

construcciÄn <strong>de</strong> un pu<strong>en</strong>te a m<strong>en</strong>udo se recurre a maquetas a escala, <strong>en</strong> lo que se podrÑa<br />

d<strong>en</strong>ominar ÅsimuladorÉ <strong>de</strong> un pu<strong>en</strong>te.<br />

En primera instancia, para validar <strong>los</strong> resultados <strong>de</strong> un esquema especÑfico <strong>de</strong><br />

producciÄn, se recurre a una planta piloto, una repres<strong>en</strong>taciÄn a escala <strong>de</strong>l proceso. Pero<br />

cuando <strong>los</strong> cambios son tan rÖpidos que no se pue<strong>de</strong> ÅarmarÉ una planta a escala, o<br />

econÄmicam<strong>en</strong>te no es conv<strong>en</strong>i<strong>en</strong>te, se pue<strong>de</strong> utilizar un simulador computacional, para recrear


8<br />

el sistema real e implem<strong>en</strong>tar y comprobar <strong>en</strong> áste <strong>los</strong> cambios que el operador consi<strong>de</strong>ra<br />

a<strong>de</strong>cuados y <strong>los</strong> efectos que t<strong>en</strong>drÑa dicho cambio <strong>en</strong> la instalaciÄn.<br />

Sin embargo utilizar un simulador no es siempre la çnica o mejor alternativa. Una<br />

bu<strong>en</strong>a forma <strong>de</strong> <strong>de</strong>cidir si un simulador es a<strong>de</strong>cuado, es consi<strong>de</strong>rar las sigui<strong>en</strong>tes situaciones:<br />

Ä Cuando no existe una formulaciÄn matemÖtica.<br />

Ä Cuando existe la formulaciÄn matemÖtica pero es difÑcil obt<strong>en</strong>er una soluciÄn<br />

analÑtica: La resoluciÄn <strong>de</strong> muchos esquemas teÄricos resulta ser tan ext<strong>en</strong>sa<br />

que a m<strong>en</strong>udo se recurre a simplificaciones que merman la exactitud <strong>de</strong>l<br />

mo<strong>de</strong>lo.<br />

Ä Si no existe el sistema real<br />

Ä Si es imposible experim<strong>en</strong>tar directam<strong>en</strong>te con el sistema real: Debido a que <strong>en</strong><br />

planta se <strong>de</strong>be mant<strong>en</strong>er un estricto control <strong>de</strong> calidad, no se pue<strong>de</strong> interv<strong>en</strong>ir<br />

el sistema las veces que sea necesario hasta <strong>en</strong>contrar la soluciÄn Äptima.<br />

Ä El sistema evoluciona muy l<strong>en</strong>tam<strong>en</strong>te, o muy rÖpidam<strong>en</strong>te: La posibilidad que<br />

ofrec<strong>en</strong> <strong>los</strong> simuladores actuales <strong>de</strong> modificar el tiempo virtual <strong>de</strong> la simulaciÄn,<br />

es un recurso invaluable para estudiar problemas <strong>de</strong> este estilo. Si por ejemplo<br />

se <strong>de</strong>be esperar una hora para observar el cambio <strong>de</strong>seado, <strong>en</strong> el simulador se<br />

pue<strong>de</strong> obt<strong>en</strong>er el resultado <strong>en</strong> unos pocos segundos.<br />

Como <strong>de</strong>sv<strong>en</strong>tajas se pued<strong>en</strong> citar:<br />

Ä El <strong>de</strong>sarrollo <strong>de</strong>l mo<strong>de</strong>lo pue<strong>de</strong> ser costoso, laborioso y l<strong>en</strong>to. Es preciso un<br />

amplio conocimi<strong>en</strong>to <strong>de</strong>l sistema a mo<strong>de</strong>lar, incluy<strong>en</strong>do consi<strong>de</strong>raciones<br />

especÑficas tomadas <strong>de</strong> plantas o diseÜos similares.<br />

Ä Existe la posibilidad <strong>de</strong> cometer errores. Nunca se <strong>de</strong>be olvidar que la simulaciÄn<br />

se lleva a cabo sobre un mo<strong>de</strong>lo, y no sobre la situaciÄn real; <strong>en</strong>tonces, si el<br />

mo<strong>de</strong>lo estÖ mal formulado, o se com<strong>en</strong>t<strong>en</strong> errores <strong>en</strong> su manejo, <strong>los</strong><br />

resultados serÖn incorrectos<br />

Ä No se pue<strong>de</strong> conocer el grado <strong>de</strong> imprecisiÄn <strong>de</strong> <strong>los</strong> resultados. Por lo g<strong>en</strong>eral, el<br />

mo<strong>de</strong>lo se utiliza para experim<strong>en</strong>tar situaciones nunca planteadas <strong>en</strong> el sistema


9<br />

real, por lo que no existe informaciÄn previa para estimar el grado <strong>de</strong><br />

correspond<strong>en</strong>cia <strong>en</strong>tre la respuesta <strong>de</strong>l mo<strong>de</strong>lo y la <strong>de</strong>l sistema real.<br />

G<strong>en</strong>eralizando, las conv<strong>en</strong>i<strong>en</strong>cias <strong>de</strong> utilizar el simulador computacional son, <strong>en</strong>tre<br />

otras, que evita complicados cÖlcu<strong>los</strong> y <strong>en</strong>trega resultados rÖpidam<strong>en</strong>te, si hay cambios que<br />

requier<strong>en</strong> <strong>de</strong> un tiempo consi<strong>de</strong>rable para observar <strong>los</strong> resultados <strong>en</strong> terr<strong>en</strong>o, <strong>en</strong> el simulador se<br />

pued<strong>en</strong> obt<strong>en</strong>er <strong>de</strong> inmediato. Y fundam<strong>en</strong>talm<strong>en</strong>te, al trabajar con sistemas abstractos, no<br />

existe el riesgo <strong>de</strong> inutilizar un equipo ni la necesidad <strong>de</strong> <strong>de</strong>t<strong>en</strong>er la producciÄn por <strong>de</strong>cisiones<br />

equivocadas.


10<br />

2.2 ETAPAS GENERALES DE UNA SIMULACIÅN<br />

FormulaciÉn <strong>de</strong>l problema: En este paso <strong>de</strong>be quedar perfectam<strong>en</strong>te establecido el<br />

objetivo <strong>de</strong> la simulaciÄn. Se <strong>de</strong>b<strong>en</strong> especificar lo mÖs <strong>de</strong>talladam<strong>en</strong>te posible <strong>los</strong> sigui<strong>en</strong>tes<br />

factores: <strong>los</strong> resultados que se esperan <strong>de</strong>l simulador, el plan <strong>de</strong> experim<strong>en</strong>taciÄn, el tiempo<br />

disponible, las variables <strong>de</strong> interás, el tipo <strong>de</strong> perturbaciones a estudiar, el tratami<strong>en</strong>to<br />

estadÑstico <strong>de</strong> <strong>los</strong> resultados, la complejidad <strong>de</strong> la interfaz <strong>de</strong>l simulador, etc. Se <strong>de</strong>be establecer<br />

si el simulador serÖ operado por el usuario o si el usuario sÄlo recibirÖ <strong>los</strong> resultados.<br />

DefiniciÉn <strong>de</strong>l sistema: El sistema a simular <strong>de</strong>be estar perfectam<strong>en</strong>te <strong>de</strong>finido, tanto<br />

<strong>en</strong> variables a <strong>de</strong>finir como <strong>en</strong> <strong>los</strong> resultados que se espera obt<strong>en</strong>er.<br />

FormulaciÉn <strong>de</strong>l mo<strong>de</strong>lo: Se comi<strong>en</strong>za con un mo<strong>de</strong>lo simple que captura <strong>los</strong> aspectos<br />

relevantes <strong>de</strong>l sistema real. Este mo<strong>de</strong>lo se irÖ <strong>en</strong>riqueci<strong>en</strong>do como resultado <strong>de</strong> varias<br />

consi<strong>de</strong>raciones adicionales que aportarÖn precisiÄn.<br />

ColecciÉn <strong>de</strong> datos: La naturaleza y cantidad <strong>de</strong> datos necesarios se <strong>de</strong>terminan<br />

directam<strong>en</strong>te por la formulaciÄn <strong>de</strong>l problema y <strong>de</strong>l mo<strong>de</strong>lo. Se pue<strong>de</strong> consi<strong>de</strong>rar como fu<strong>en</strong>te<br />

<strong>de</strong> datos tanto registros histÄricos y mediciones <strong>de</strong> laboratorio como observaciones realizadas<br />

<strong>en</strong> el sistema real. Estos datos <strong>de</strong>berÖn ser procesados a<strong>de</strong>cuadam<strong>en</strong>te para darles el formato<br />

exigido <strong>en</strong> el mo<strong>de</strong>lo.<br />

Implem<strong>en</strong>taciÉn <strong>de</strong>l mo<strong>de</strong>lo <strong>en</strong> el ord<strong>en</strong>ador: Para el caso <strong>de</strong>l simulador HYSYS, el<br />

mismo programa cu<strong>en</strong>ta con mo<strong>de</strong><strong>los</strong> pre<strong>de</strong>terminados que sirv<strong>en</strong> <strong>de</strong> base para realizar las<br />

modificaciones necesarias.<br />

VerificaciÉn: En esta etapa se comprueba que no se haya cometido errores durante la<br />

implem<strong>en</strong>taciÄn <strong>de</strong>l mo<strong>de</strong>lo. Se <strong>de</strong>be revisar cada cÖlculo, estructura <strong>de</strong> programaciÄn o<br />

mátodo utilizado.


11<br />

ValidaciÉn: En esta etapa se comprueba la exactitud <strong>de</strong>l mo<strong>de</strong>lo <strong>de</strong>sarrollado. Esto se<br />

lleva a cabo comparando las predicciones <strong>de</strong>l mo<strong>de</strong>lo con mediciones realizadas <strong>en</strong> el sistema<br />

real, datos histÄricos o datos <strong>de</strong> sistemas similares. Como resultado <strong>de</strong> esta etapa pue<strong>de</strong> surgir<br />

la necesidad <strong>de</strong> modificar el mo<strong>de</strong>lo o recolectar datos adicionales.<br />

DiseÑo <strong>de</strong> experim<strong>en</strong>tos: En esta etapa se <strong>de</strong>cid<strong>en</strong> las caracterÑsticas <strong>de</strong> <strong>los</strong><br />

experim<strong>en</strong>tos a realizar; el tiempo <strong>de</strong> arranque, el tiempo <strong>de</strong> simulaciÄn y el nçmero <strong>de</strong><br />

simulaciones necesarias.<br />

Experim<strong>en</strong>taciÉn: Se realizan las simulaciones <strong>de</strong> acuerdo con el diseÜo previo. Los<br />

resultados son <strong>de</strong>bidam<strong>en</strong>te recolectados y procesados.<br />

InterpretaciÉn: Se analiza la s<strong>en</strong>sibilidad <strong>de</strong>l mo<strong>de</strong>lo con respecto a <strong>los</strong> parÖmetros que<br />

ti<strong>en</strong><strong>en</strong> asociados una mayor incertidumbre. El mo<strong>de</strong>lo serÖ s<strong>en</strong>sible a <strong>de</strong>terminados parÖmetros<br />

si ante pequeÜos cambios <strong>en</strong> <strong>los</strong> valores <strong>de</strong> <strong>los</strong> mismos, las respuestas varÑan notablem<strong>en</strong>te.<br />

Implem<strong>en</strong>taciÉn: Se difun<strong>de</strong> la simulaciÄn y el mo<strong>de</strong>lo obt<strong>en</strong>idos. El responsable <strong>de</strong> la<br />

simulaciÄn <strong>de</strong>be guiar <strong>en</strong> esta etapa, para evitar que <strong>los</strong> resultados se utilic<strong>en</strong> mÖs allÖ <strong>de</strong>l rango<br />

<strong>de</strong> aplicaciÄn consi<strong>de</strong>rado <strong>en</strong> el estudio.<br />

Docum<strong>en</strong>taciÉn: Elaborar la docum<strong>en</strong>taciÄn tácnica, con una <strong>de</strong>scripciÄn <strong>de</strong>tallada <strong>de</strong>l<br />

mo<strong>de</strong>lo y <strong>de</strong> <strong>los</strong> datos, y manuales <strong>de</strong> uso con las consi<strong>de</strong>raciones particulares <strong>de</strong> cada sistema.


12<br />

2.3 EL SIMULADOR HYSYS Ç<br />

HYSYS Ç es un software <strong>de</strong>sarrollado por la empresa Asp<strong>en</strong>Tech, <strong>de</strong>l cual tanto ENAP<br />

Magallanes como el Departam<strong>en</strong>to <strong>de</strong> QuÑmica <strong>de</strong> la Universidad <strong>de</strong> Magallanes pose<strong>en</strong> lic<strong>en</strong>cia<br />

<strong>de</strong> uso para su versiÄn 3.2, <strong>en</strong> estado estacionario.<br />

HYSYS se ha mant<strong>en</strong>ido durante 20 aÜos como la alternativa lÑ<strong>de</strong>r <strong>en</strong> simulaciones<br />

relacionadas con la industria petroquÑmica, y <strong>en</strong> ENAP se utilizan a diario simulaciones creadas<br />

por <strong>los</strong> ing<strong>en</strong>ieros <strong>de</strong> <strong>procesos</strong>, para estimar soluciones y planificar futuros proyectos.<br />

HYSYS es un programa <strong>de</strong> simulaciÄn <strong>de</strong> <strong>procesos</strong> quÑmicos ori<strong>en</strong>tado a objetos. Posee<br />

un aspecto visual caracterÑstico y trabaja sobre sistemas operativos Windows. Permite mo<strong>de</strong>lar<br />

sistemas complejos mediante una avanzada interfaz grÖfica, que ofrece al usuario la posibilidad<br />

<strong>de</strong> ÅarmarÉ el flowsheet <strong>de</strong> planta <strong>en</strong> una v<strong>en</strong>tana llamada PFD, Process Flowsheet Diagram.<br />

Posteriorm<strong>en</strong>te ingresando <strong>los</strong> datos necesarios (presiones, temperaturas, composiciones, flujos<br />

y especificaciones tácnicas <strong>de</strong> equipos, bÖsicam<strong>en</strong>te), se pue<strong>de</strong> simular el sistema real y obt<strong>en</strong>er<br />

<strong>los</strong> resultados o estimaciones a<strong>de</strong>cuadas.<br />

Todas las consi<strong>de</strong>raciones <strong>de</strong> manejo y creaciÄn <strong>de</strong> programas se <strong>en</strong>cu<strong>en</strong>tran<br />

totalm<strong>en</strong>te <strong>de</strong>talladas <strong>en</strong> el ÅManual para simular <strong>procesos</strong> <strong>productivos</strong> <strong>de</strong> ENAP Magallanes <strong>en</strong><br />

HYSYSÉ, que se adjunta. RazÄn por la cual sÄlo se esbozarÖ brevem<strong>en</strong>te <strong>los</strong> aspectos bÖsicos <strong>de</strong><br />

uso <strong>de</strong> HYSYS.<br />

El primer paso necesario es <strong>de</strong>finir <strong>los</strong> compon<strong>en</strong>tes que estÖn pres<strong>en</strong>tes <strong>en</strong> la<br />

simulaciÄn, y el paquete termodinÖmico asociado. Esto se realiza <strong>en</strong> una secciÄn llamada<br />

Simulation Basis Manager. HYSYS provee <strong>de</strong> informaciÄn <strong>de</strong>tallada para bastantes mo<strong>de</strong><strong>los</strong><br />

termodinÖmicos, e inclusive avisa al usuario cuando se escog<strong>en</strong> compon<strong>en</strong>tes que no pued<strong>en</strong><br />

ser satisfactoriam<strong>en</strong>te mo<strong>de</strong>lados mediante el sistema termodinÖmico seleccionado. A<strong>de</strong>mÖs, si<br />

se <strong>de</strong>sea cambiar <strong>de</strong> termodinÖmica para ciertos equipos o secciones <strong>de</strong> la simulaciÄn, se da la<br />

posibilidad <strong>de</strong> elegir diversos paquetes termodinÖmicos y asociar<strong>los</strong> a difer<strong>en</strong>tes listados <strong>de</strong><br />

compon<strong>en</strong>tes.


13<br />

IlustraciÄn 3.- Simulation Basis Manager <strong>de</strong> HYSYS<br />

Posteriorm<strong>en</strong>te se ingresa al Simulation Envirom<strong>en</strong>t, que conjuga las particularida<strong>de</strong>s<br />

<strong>de</strong> tres tipos <strong>de</strong> v<strong>en</strong>tanas (PFD, Workbook y v<strong>en</strong>tanas <strong>de</strong> propieda<strong>de</strong>s <strong>de</strong> cada objeto) y <strong>en</strong> las<br />

que es posible ingresar informaciÄn <strong>de</strong> la manera que mÖs acomo<strong>de</strong> al usuario.<br />

Una vez seleccionado el mo<strong>de</strong>lo asociado con el sistema real, el posterior paso es<br />

crearlo <strong>en</strong> el PFD, tal como si se diseÜara un plano, asignando nombres a equipos y corri<strong>en</strong>tes<br />

<strong>de</strong> proceso, creando las conexiones y asignando <strong>los</strong> datos que se han recopilado a las<br />

respectivas casillas que HYSYS pres<strong>en</strong>ta para tal fin.<br />

IlustraciÄn 4.- V<strong>en</strong>tana <strong>de</strong> propieda<strong>de</strong>s para un intercambiador <strong>de</strong> tubo y carcasa


14<br />

Una vez <strong>de</strong>finidas todas las corri<strong>en</strong>tes <strong>de</strong> proceso y caracterizadas todas las condiciones<br />

y variables relativas a <strong>los</strong> difer<strong>en</strong>tes equipos que toman parte <strong>en</strong> la simulaciÄn, se pued<strong>en</strong><br />

observar directam<strong>en</strong>te <strong>los</strong> resultados, ya que HYSYS posee un sistema <strong>de</strong> resoluciÄn que por<br />

<strong>de</strong>fecto se <strong>en</strong>cu<strong>en</strong>tra siempre activado, lo que significa que el programa calcula todas las<br />

propieda<strong>de</strong>s y resultados factibles <strong>de</strong> calcular <strong>en</strong> todo mom<strong>en</strong>to, y una vez <strong>de</strong>finida por<br />

completo la simulaciÄn, <strong>en</strong> el instante <strong>en</strong> que el usuario <strong>de</strong>fina la çltima variable, el sistema<br />

automÖticam<strong>en</strong>te calcularÖ y <strong>en</strong>tregarÖ <strong>los</strong> resultados. Obviam<strong>en</strong>te esta Årespuesta automÖticaÉ<br />

se podrÑa ev<strong>en</strong>tualm<strong>en</strong>te <strong>de</strong>sactivar, pero dada la actual velocidad <strong>de</strong> <strong>los</strong> ord<strong>en</strong>adores, se vuelve<br />

innecesario.<br />

Cuando el sistema estÖ completam<strong>en</strong>te <strong>de</strong>finido, se <strong>de</strong>be asignar nuevos valores a las<br />

variables especificadas, y conoci<strong>en</strong>do las respuestas reales <strong>de</strong>l sistema ante dichos cambios, se<br />

pue<strong>de</strong> monitorear el grado <strong>de</strong> similitud que pres<strong>en</strong>tan ambos sistemas, el real con el simulado.<br />

Si la similitud obt<strong>en</strong>ida es satisfactoria, se pue<strong>de</strong> concluir que la simulaciÄn HYSYS sirve para<br />

estimar condiciones futuras <strong>de</strong> proceso, y <strong>en</strong>tregarÖ estimaciones a<strong>de</strong>cuadas para analizar cÄmo<br />

se comportarÖ el sistema real ante <strong>los</strong> cambios ingresados.<br />

IlustraciÄn 5.- Aspecto <strong>de</strong> una simulaciÄn tÑpica <strong>en</strong> HYSYS


15<br />

2.4 CONSIDERACIONES TEÅRICAS UTILIZADAS POR HYSYS<br />

Consi<strong>de</strong>rando la totalidad <strong>de</strong> simulaciones realizadas <strong>en</strong> este trabajo <strong>de</strong> titulaciÄn, se<br />

hace uso <strong>de</strong> muchos equipos y propieda<strong>de</strong>s disponibles <strong>en</strong> HYSYS. Algunos <strong>de</strong> <strong>los</strong> objetos<br />

utilizados simulan equipos efectivam<strong>en</strong>te instalados <strong>en</strong> terr<strong>en</strong>o, como compresores o<br />

intercambiadores <strong>de</strong> calor, sin embargo otras propieda<strong>de</strong>s son <strong>de</strong> uso exclusivo <strong>de</strong>l simulador,<br />

como ajustadores, planillas <strong>de</strong> cÖlculo o recic<strong>los</strong>.<br />

En la imposibilidad <strong>de</strong> <strong>de</strong>tallar <strong>en</strong> el pres<strong>en</strong>te informe cada uno <strong>de</strong> <strong>los</strong> equipos<br />

utilizados, y las respectivas consi<strong>de</strong>raciones teÄricas subyac<strong>en</strong>tes, se han elegido tres<br />

operaciones que son repres<strong>en</strong>tativas <strong>de</strong> cÄmo HYSYS interpreta <strong>los</strong> mo<strong>de</strong><strong>los</strong> teÄricos<br />

disponibles y <strong>los</strong> adapta a su propio l<strong>en</strong>guaje <strong>de</strong> simulaciÄn. Las operaciones correspond<strong>en</strong> a:<br />

Ä Compresores recÑprocos.<br />

Ä Intercambiadores <strong>de</strong> calor <strong>de</strong> tubo y carcasa.<br />

Ä Columnas <strong>de</strong> <strong>de</strong>stilaciÄn.<br />

Y como parÖmetros propios <strong>de</strong> HYSYS, se <strong>de</strong>sarrollarÖn <strong>los</strong> aspectos teÄricos <strong>de</strong> dos<br />

herrami<strong>en</strong>tas fundam<strong>en</strong>tales a la hora <strong>de</strong> armar simulaciones:<br />

Ä Adjust.<br />

Ä Recycle.<br />

Se <strong>de</strong>be consi<strong>de</strong>rar a<strong>de</strong>mÖs que HYSYS es un producto comercial cuyo cÄdigo fu<strong>en</strong>te es<br />

confid<strong>en</strong>cial, razÄn por la cual no se pued<strong>en</strong> <strong>de</strong>tallar <strong>en</strong> profundidad las consi<strong>de</strong>raciones<br />

teÄricas o algoritmos numáricos <strong>de</strong> soluciÄn que utiliza el programa. Como çnica fu<strong>en</strong>te <strong>de</strong><br />

informaciÄn se ha consi<strong>de</strong>rado la pres<strong>en</strong>tada <strong>en</strong> <strong>los</strong> manuales que acompaÜan al programa.


16<br />

2.4.1 COMPRESORES RECÄPROCOS<br />

Un compresor recÑproco es un equipo que mediante la adiciÄn <strong>de</strong> <strong>en</strong>ergÑa,<br />

normalm<strong>en</strong>te cedida por un motor acoplado, eleva la presiÄn <strong>de</strong>l fluido que ingresa a áste.<br />

El compresor consta <strong>de</strong> un cilindro por el cual avanza un pistÄn, que comprime el<br />

fluido. Si <strong>en</strong> el movimi<strong>en</strong>to <strong>de</strong> retroceso el pistÄn a<strong>de</strong>mÖs comprime el fluido que ingresa por la<br />

parte trasera (cÖrter), <strong>en</strong>tonces se le d<strong>en</strong>omina pistÄn <strong>de</strong> doble efecto. Este tipo <strong>de</strong><br />

compresores es el mÖs utilizado <strong>en</strong> las instalaciones <strong>de</strong> ENAP Magallanes.<br />

IlustraciÄn 6.- Esquema <strong>de</strong>l cilindro <strong>de</strong> un compresor recÑproco<br />

recÑproco son:<br />

Las variables que son <strong>de</strong> interás para estimar el comportami<strong>en</strong>to <strong>de</strong> un compresor<br />

Ef<br />

adiab<br />

Pot<strong>en</strong>cia adiabÄtica requerida<br />

(%) Å<br />

Ä100%<br />

Pot<strong>en</strong>cia real requerida<br />

[Ec. 1]<br />

La pot<strong>en</strong>cia adiabÖtica requerida correspon<strong>de</strong> al trabajo mecÖnicam<strong>en</strong>te reversible W:<br />

W<br />

P2<br />

Å ÇVdP<br />

P1<br />

[Ec. 2]


17<br />

Para este caso HYSYS calcula el resultado utilizando directam<strong>en</strong>te las ecuaciones <strong>de</strong>l<br />

paquete termodinÖmico que se haya seleccionado. La pot<strong>en</strong>cia real requerida correspon<strong>de</strong> a la<br />

difer<strong>en</strong>cia <strong>de</strong> <strong>en</strong>talpÑas que, si se d<strong>en</strong>omina como H, significa:<br />

Pot<strong>en</strong>cia real requerida Å H É H<br />

Salida<br />

Entrada<br />

[Ec. 3]<br />

En el caso <strong>en</strong> que se conozcan las presiones <strong>de</strong> <strong>en</strong>trada y salida, la temperatura <strong>de</strong><br />

<strong>en</strong>trada y la efici<strong>en</strong>cia, la pot<strong>en</strong>cia real requerida se calcula como la divisiÄn <strong>en</strong>tre la pot<strong>en</strong>cia<br />

adiabÖtica calculada y la efici<strong>en</strong>cia especificada. Posteriorm<strong>en</strong>te, utilizando el mátodo<br />

termodinÖmico seleccionado por el usuario, se calcula la temperatura <strong>de</strong> salida que satisface la<br />

ecuaciÄn [Ec. 3].<br />

Para el cÖlculo <strong>de</strong> la efici<strong>en</strong>cia politrÄpica, ásta se obti<strong>en</strong>e <strong>de</strong> la sigui<strong>en</strong>te ecuaciÄn:<br />

Ef<br />

Ñ n Ö<br />

ä Ü á ã<br />

n 1<br />

P<br />

à É â<br />

åÑ Ö<br />

sal<br />

äÑ n Ö Ñ k É1Öã<br />

É1 ç<br />

Ä Ä<br />

åÜ á<br />

P ç åÜ á Ü á<br />

<strong>en</strong>t<br />

n 1 k<br />

ç<br />

éà É â à âè<br />

å<br />

à â<br />

ç<br />

Å<br />

é<br />

è<br />

ÄEf<br />

pol. Ñ k É1<br />

Ö<br />

adiab.<br />

ä Ü á ã<br />

Ñ<br />

k<br />

P Öà<br />

â<br />

sal<br />

å<br />

É1<br />

ç<br />

åÜ<br />

á<br />

P ç<br />

à <strong>en</strong>t<br />

å<br />

â<br />

é çè [Ec. 4]<br />

Don<strong>de</strong><br />

log( Psal<br />

/ P<strong>en</strong>t<br />

)<br />

n Å<br />

actual<br />

log( Ä / Ä )<br />

sal<br />

<strong>en</strong>t<br />

y<br />

log( Psal<br />

/ P<strong>en</strong>t<br />

)<br />

k Å<br />

i<strong>de</strong>al<br />

log( Ä / Ä )<br />

sal<br />

<strong>en</strong>t<br />

Las ecuaciones ya pres<strong>en</strong>tadas son comunes a cualquier compresor, ya sea c<strong>en</strong>trÑfugo o<br />

recÑproco. Adicionalm<strong>en</strong>te HYSYS, para el cÖlculo <strong>de</strong> un compresor recÑproco, evalça tres<br />

parÖmetros mÖs: Volum<strong>en</strong> <strong>de</strong>splazado por cada cilindro (PD), Clearance <strong>de</strong> <strong>los</strong> cilindros (Cl) y<br />

efici<strong>en</strong>cia volumátrica (VE).


18<br />

ê Cl cada cilindro<br />

Cl Å<br />

PD [Ec. 5]<br />

PD se calcula como el producto <strong>en</strong>tre el Örea transversal neta <strong>de</strong> compresiÄn <strong>de</strong>l<br />

cilindro y la longitud o carrera (Stroke) que se <strong>de</strong>splaza el pistÄn.<br />

1<br />

ä ä<br />

ãã<br />

k<br />

Z Ñ<br />

s<br />

P Ö<br />

d<br />

VE Å<br />

å<br />

ë1 Lí<br />

C å 1ç<br />

ç<br />

å<br />

É É É<br />

å Ü á Z<br />

d P ç<br />

à<br />

ç<br />

s<br />

å<br />

â<br />

é åé<br />

çè<br />

çè [Ec. 6]<br />

Pd<br />

Å PresiÅn <strong>de</strong> <strong>de</strong>scarga<br />

P Å PresiÅn <strong>de</strong> succiÅn<br />

s<br />

L Å Efectos <strong>de</strong> friccion, pÇrdidas <strong>de</strong> carga <strong>en</strong> vÉlvulas, filtraciones<br />

k Å Cp / Cv<br />

Zd Å factor <strong>de</strong> compresibilidad <strong>en</strong>trada<br />

Zs Å factor <strong>de</strong> compresibilidad salida<br />

C Å Volum<strong>en</strong> <strong>de</strong> clearance<br />

El flujo comprimido (F) es funciÄn directa <strong>de</strong> la velocidad <strong>de</strong> rotaciÄn <strong>de</strong>l compresor, se<br />

calcula como flujo molar y estÖ <strong>de</strong>finido por:<br />

1<br />

ä<br />

ä ãã<br />

ä N ã<br />

k Ä PD ÄÄ<br />

Ñ L Ö Zs<br />

Ñ Pd<br />

Ö å<br />

F 1 C 1 Ä 60 ç<br />

Å<br />

å<br />

å<br />

çç<br />

åÜ<br />

É á É Ü á É å ç<br />

à 100 â å Z<br />

d P çç<br />

à s â å PM<br />

å<br />

å<br />

çç<br />

ç<br />

é é<br />

èè é è [Ec. 7]<br />

N Å Velocidad <strong>de</strong> rotaciÅn (rpm)<br />

Ä=D<strong>en</strong>sidad <strong>de</strong>l gas<br />

PM = Peso molecular <strong>de</strong>l gas


19<br />

En el caso especÑfico <strong>en</strong> que la velocidad <strong>de</strong> rotaciÄn <strong>de</strong>l compresor es exactam<strong>en</strong>te<br />

cero, el flujo a travás <strong>de</strong>l compresor se calcula mediante una fÄrmula <strong>de</strong> presiÄn-flujo, don<strong>de</strong><br />

especificando la resist<strong>en</strong>cia a velocidad cero, k vel.cero , se calcula como:<br />

F Å kvel. ceroÄ Ä<br />

ìP<br />

fricciÅn<br />

[Ec. 8]<br />

ìP fricciÅn<br />

Å PÇrdida <strong>de</strong> carga por efectos <strong>de</strong> fricciÅn<br />

La presiÄn mÖxima <strong>de</strong> <strong>de</strong>scarga que se pue<strong>de</strong> obt<strong>en</strong>er <strong>en</strong> un compresor recÑproco es:<br />

mÄx<br />

ä Z<br />

ã<br />

d<br />

Pd<br />

Å Ps<br />

Äå<br />

Ä(1 É L ÉVE î C)<br />

ç<br />

é ZsÄC<br />

è [Ec. 9]<br />

k<br />

En g<strong>en</strong>eral HYSYS resuelve presiÄn y flujo <strong>de</strong> forma in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te. Consi<strong>de</strong>rando la<br />

corri<strong>en</strong>te <strong>de</strong> <strong>en</strong>trada completam<strong>en</strong>te <strong>de</strong>finida, si se especifica la presiÄn <strong>de</strong> <strong>de</strong>scarga y la<br />

efici<strong>en</strong>cia, se obti<strong>en</strong>e la <strong>en</strong>ergÑa necesaria, temperatura <strong>de</strong> <strong>de</strong>scarga y la velocidad <strong>de</strong> rotaciÄn.<br />

Si se especifica la velocidad <strong>de</strong> rotaciÄn y la <strong>en</strong>ergÑa ingresada, se calcula la presiÄn y<br />

temperatura <strong>de</strong> <strong>de</strong>scarga y el flujo. Pero no se pue<strong>de</strong> especificar el flujo y la velocidad, ambos a<br />

la vez.<br />

2.4.2 INTERCAMBIADORES DE CALOR DE TUBO Y CARCASA<br />

Entre todos <strong>los</strong> tipos <strong>de</strong> intercambiadores <strong>de</strong> calor posibles <strong>de</strong> utilizar <strong>en</strong> HYSYS, el mÖs<br />

utilizado para las simulaciones <strong>de</strong> ENAP Magallanes es el <strong>de</strong> tubo y carcasa.<br />

Este tipo <strong>de</strong> intercambiador <strong>de</strong> calor realiza balances <strong>de</strong> materia y <strong>en</strong>ergÑa para ambas<br />

corri<strong>en</strong>tes<br />

simultÖneam<strong>en</strong>te. Como ya se ha especificado, el trabajo se <strong>de</strong>sarrollÄ


20<br />

exclusivam<strong>en</strong>te <strong>en</strong> estado estacionario (Steady-State), por lo cual las ecuaciones que sigu<strong>en</strong><br />

tratarÖn exclusivam<strong>en</strong>te dicho aspecto.<br />

El balance g<strong>en</strong>eral aplicado por HYSYS al intercambiador correspon<strong>de</strong> a:<br />

ë í Äë í<br />

ä M frÇoÄ<br />

H sal<br />

É H <strong>en</strong>t<br />

É Q ã<br />

frÇo filtrado<br />

É ä M H H Q Balance<br />

cali<strong>en</strong>te <strong>en</strong>t<br />

É ã<br />

sal<br />

É<br />

cali<strong>en</strong>te perdido<br />

Å<br />

é è é<br />

è<br />

error<br />

[Ec. 10]<br />

M Å Flujo mÉsico <strong>de</strong>l fluido<br />

H Å EntalpÑa<br />

El parÖmetro Balance error es una especificaciÄn directa a HYSYS, que ti<strong>en</strong>e unida<strong>de</strong>s <strong>de</strong><br />

flujo <strong>de</strong> calor (<strong>en</strong>ergÑa por tiempo) y por lo g<strong>en</strong>eral se asume igual a cero, pero se estÖ <strong>en</strong><br />

libertad <strong>de</strong> especificar un valor difer<strong>en</strong>te.<br />

El calor total transferido <strong>en</strong>tre <strong>los</strong> tubos y la carcasa, se pue<strong>de</strong> <strong>de</strong>finir <strong>en</strong> tárminos <strong>de</strong>l<br />

coefici<strong>en</strong>te global <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> calor, el Örea <strong>de</strong> transfer<strong>en</strong>cia y la temperatura media<br />

logarÑtmica <strong>de</strong> acuerdo a la sigui<strong>en</strong>te ecuaciÄn:<br />

Q Å UÄAÄìT LM<br />

ÄFt<br />

[Ec. 11]<br />

U Å Coefici<strong>en</strong>te global <strong>de</strong> transfer<strong>en</strong>cia <strong>de</strong> calor<br />

A Å Örea total <strong>de</strong> transfer<strong>en</strong>cia<br />

ìT<br />

t<br />

LM<br />

Å Difer<strong>en</strong>cia <strong>de</strong> temperaturas media logarÑtmica (LMTD)<br />

F Å Factor <strong>de</strong> correciÅn <strong>de</strong> LMTD<br />

La LMTD se calcula como:<br />

ìT<br />

LM<br />

ìT1 É ìT2<br />

Å<br />

ln( ìT1 / ìT2<br />

)<br />

[Ec. 12]


21<br />

Don<strong>de</strong><br />

ìT Å T ÉT<br />

1<br />

2<br />

salida<br />

cali<strong>en</strong>te<br />

ìT Å T ÉT<br />

<strong>en</strong>trada<br />

cali<strong>en</strong>te<br />

<strong>en</strong>trada<br />

frÇo<br />

salida<br />

frÇo<br />

Se pue<strong>de</strong> elegir si calcular el factor F t al especificar la configuraciÄn especÑfica <strong>de</strong>l<br />

intercambiador, o asumirlo igual a 1. Se suele agrupar <strong>los</strong> factores U y A y <strong>de</strong>finir directam<strong>en</strong>te<br />

el valor para UA <strong>en</strong> conjunto.<br />

La párdida <strong>de</strong> carga se pue<strong>de</strong> calcular <strong>de</strong> tres maneras:<br />

Ä EspecificaciÄn directa.<br />

Ä CÖlculo automÖtico <strong>de</strong> HYSYS <strong>de</strong> acuerdo a la configuraciÄn y geometrÑa <strong>de</strong>l<br />

intercambiador.<br />

Ä Defini<strong>en</strong>do una relaciÄn presiÄn-flujo y especificando un factor k.<br />

Para la çltima opciÄn, la fÄrmula es exactam<strong>en</strong>te la [Ec. 8].<br />

Los mátodos <strong>de</strong> resoluciÄn, y una breve <strong>de</strong>scripciÄn <strong>de</strong> ástos, son:<br />

Ä End-Point: Realiza el balance consi<strong>de</strong>rando sÄlo puntos <strong>de</strong> <strong>en</strong>trada y salida.<br />

Consi<strong>de</strong>ra un factor UA y valores <strong>de</strong> Cp constantes para todo el intercambiador.<br />

Para esta opciÄn, HYSYS pue<strong>de</strong> calcular el factor F t <strong>en</strong> funciÄn <strong>de</strong> <strong>los</strong> parÖmetros<br />

fÑsicos <strong>de</strong>l intercambiador (nçmero <strong>de</strong> pasos por <strong>los</strong> tubos, contracorri<strong>en</strong>te o<br />

paralelo, etc.).<br />

Ä Weighted: Divi<strong>de</strong> el intercambiador <strong>en</strong> interva<strong>los</strong> pequeÜos, <strong>en</strong> cada intervalo se<br />

calculan <strong>los</strong> difer<strong>en</strong>tes parÖmetros <strong>de</strong> la ecuaciÄn [Ec. 11] (LMTD, UA, etc).<br />

Sumam<strong>en</strong>te çtil para intercambiadores con alto rango <strong>de</strong> temperaturas,<br />

cond<strong>en</strong>saciÄn, etc. SÄlo calcula el factor F t para intercambiadores <strong>en</strong><br />

contracorri<strong>en</strong>te.<br />

Ä Steady-State Rating: Para ser usado <strong>en</strong> estado estacionario. Realiza las mismas<br />

consi<strong>de</strong>raciones que el modo End Point, pero a<strong>de</strong>mÖs permite estudiar el<br />

diseÜo <strong>de</strong>l equipo (modo Rating).


22<br />

Ä Dynamic: Estudia el comportami<strong>en</strong>to dinÖmico <strong>de</strong>l intercambiador. No<br />

consi<strong>de</strong>rado <strong>en</strong> este trabajo.<br />

IlustraciÄn 7.- EspecificaciÄn <strong>en</strong> modo Weighted<br />

Se pue<strong>de</strong> <strong>de</strong>cir que <strong>en</strong> g<strong>en</strong>eral, para <strong>los</strong> objetivos que se <strong>de</strong>sarrollan <strong>en</strong> este trabajo, no<br />

ha sido necesario modificar el modo que por <strong>de</strong>fecto <strong>en</strong>trega HYSYS (End-Point), ya que <strong>los</strong><br />

intercambiadores no juegan un papel clave <strong>en</strong> <strong>los</strong> <strong>procesos</strong> investigados, y las difer<strong>en</strong>cias <strong>en</strong> <strong>los</strong><br />

valores <strong>de</strong> temperatura obt<strong>en</strong>idos prÖcticam<strong>en</strong>te no son <strong>de</strong>cisivas.<br />

IlustraciÄn 8.- Ingreso <strong>de</strong> parÖmetros fÑsicos <strong>de</strong>l intercambiador


23<br />

Los parÖmetros fÑsicos posibles <strong>de</strong> ingresar cubr<strong>en</strong> la totalidad <strong>de</strong> variables utilizadas <strong>en</strong><br />

<strong>los</strong> cÖlcu<strong>los</strong> acostumbrados <strong>de</strong> intercambiadores <strong>de</strong> calor, como configuraciÄn <strong>de</strong> <strong>los</strong> tubos,<br />

nçmero <strong>de</strong> pasos por la carcasa, pitch, espaciado <strong>de</strong> <strong>los</strong> baffles, etc.<br />

En modo Steady-State Rating, a<strong>de</strong>mÖs, se calculan <strong>los</strong> coefici<strong>en</strong>tes locales <strong>de</strong><br />

transfer<strong>en</strong>cia <strong>de</strong> calor, para el fluido <strong>de</strong> <strong>los</strong> tubos y carcasa, <strong>de</strong> acuerdo a la correlaciÄn:<br />

h<br />

i<br />

0.8 1/3 0.14<br />

0.027k m<br />

Ñ D i<br />

G i<br />

Ö Ñ Cp i<br />

Å i<br />

Ö Ñ Å i<br />

Ö<br />

w<br />

Di Åi km Åi<br />

Å Ü á Ü á Ü á<br />

à â à â à â [Ec. 13]<br />

G Å Velocidad mÉsica <strong>de</strong>l fluido <strong>en</strong> <strong>los</strong> tubos (velocidadÄd<strong>en</strong>sidad)<br />

i<br />

Å Å Viscosidad <strong>de</strong>l fluido <strong>de</strong> <strong>los</strong> tubos<br />

i<br />

Å Å Viscosidad <strong>de</strong>l fluido <strong>en</strong> las pare<strong>de</strong>s <strong>de</strong>l tubo<br />

w<br />

i<br />

Cp Å Capacidad calorÑfica <strong>de</strong>l fluido <strong>en</strong> <strong>los</strong> tubos<br />

i<br />

Y la relaciÄn que utiliza HYSYS <strong>en</strong>tre estos coefici<strong>en</strong>tes locales y el factor U es:<br />

U<br />

1<br />

Å<br />

D Ñ<br />

o<br />

1 Ö<br />

h0<br />

î ro î rw î Ü ri<br />

î á<br />

Di<br />

à hi<br />

â [Ec. 14]<br />

h0 Å Coefici<strong>en</strong>te local <strong>de</strong> transf. <strong>de</strong> calor por la carcasa<br />

r Å Factor <strong>de</strong> <strong>en</strong>suciami<strong>en</strong>to por la carcasa<br />

r<br />

o<br />

w<br />

D<br />

Å Resist<strong>en</strong>cia <strong>de</strong> las pare<strong>de</strong>s <strong>de</strong> <strong>los</strong> tubos<br />

o<br />

i<br />

Å DiÉmetro exterior <strong>de</strong> <strong>los</strong> tubos<br />

D Å DiÉmetro interior <strong>de</strong> <strong>los</strong> tubos<br />

r Å Factor <strong>de</strong> <strong>en</strong>suciami<strong>en</strong>to por d<strong>en</strong>tro <strong>de</strong> <strong>los</strong> tubos<br />

i<br />

h Å Coefici<strong>en</strong>te <strong>de</strong> transf. <strong>de</strong> calor por d<strong>en</strong>tro <strong>de</strong> <strong>los</strong> tubos<br />

i<br />

Exist<strong>en</strong> muchas mÖs opciones para especificar un intercambiador <strong>de</strong> tubo y carcasa, y<br />

Asp<strong>en</strong>Tech ofrece mçltiples programas especÑficos para intercambiadores, que pued<strong>en</strong> ser<br />

aÜadidos a HYSYS para dar mayor exactitud a la simulaciÄn, y serÑa poco prÖctico seguir


24<br />

<strong>de</strong>tallando cada uno <strong>de</strong> <strong>los</strong> botones o parÖmetros. Por otra parte, <strong>los</strong> resultados obt<strong>en</strong>idos <strong>en</strong><br />

modo End-Point pres<strong>en</strong>tan difer<strong>en</strong>cias poco apreciables con respecto a utilizar todas las<br />

alternativas <strong>de</strong> especificaciÄn que se puedan ingresar <strong>en</strong> Steady-State Rating, y no justifican,<br />

para <strong>los</strong> objetivos que ENAP Magallanes requiere, el gasto <strong>de</strong> tiempo y esfuerzo necesario.<br />

2.4.3 REACCIONES QUÄMICAS Y REACTORES EN HYSYS<br />

Si bi<strong>en</strong> <strong>en</strong> las instalaciones <strong>de</strong> ENAP Magallanes no se llevan a cabo reacciones quÑmicas<br />

como parte <strong>de</strong> <strong>los</strong> <strong>procesos</strong> <strong>productivos</strong>, una <strong>de</strong> las simulaciones <strong>de</strong>sarrolladas correspon<strong>de</strong> a la<br />

turbina <strong>de</strong>l tr<strong>en</strong> SULZER, <strong>en</strong> la que se ha simulado la combustiÄn <strong>de</strong>l gas natural que g<strong>en</strong>era la<br />

mezcla que ingresa a turbina mediante dos opciones; <strong>de</strong>fini<strong>en</strong>do las reacciones <strong>de</strong> combustiÄn y<br />

utilizando un reactor pre<strong>de</strong>finido <strong>de</strong> HYSYS d<strong>en</strong>ominado Gibbs Reactor. A continuaciÄn se<br />

especifican ambos mátodos.<br />

En g<strong>en</strong>eral, para <strong>de</strong>finir reacciones quÑmicas, HYSYS ofrece un mÄdulo especÑfico, con<br />

algunas reacciones pre<strong>de</strong>finidas, tÑpicas <strong>de</strong> la industria quÑmica, como pue<strong>de</strong> observar a<br />

continuaciÄn:<br />

IlustraciÄn 9.- Algunas reacciones pre<strong>de</strong>finidas <strong>en</strong> HYSYS<br />

HYSYS divi<strong>de</strong> las reacciones <strong>en</strong> cinco grupos: Conversion, Equilibrium, Heterog<strong>en</strong>eous<br />

Catalytic, Kinetic y Simple Rate. Una breve <strong>de</strong>scripciÄn es:


25<br />

Ä Conversion: Requiere la estequiometrÑa <strong>de</strong> las reacciones y la conversiÄn <strong>en</strong><br />

funciÄn <strong>de</strong> un compon<strong>en</strong>te <strong>de</strong> la reacciÄn.<br />

Ä Equilibrium: Requiere la estequiometrÑa y ord<strong>en</strong> <strong>de</strong> las reacciones, y pres<strong>en</strong>ta<br />

difer<strong>en</strong>tes formas <strong>de</strong> calcular o especificar la constante <strong>de</strong> equilibrio ln(K).<br />

Ä Heterog<strong>en</strong>eous Catalytic: Requiere estequiometrÑa y parÖmetros cináticos <strong>de</strong> la<br />

reacciÄn, como la EnergÑa <strong>de</strong> ActivaciÄn, Factor <strong>de</strong> Frecu<strong>en</strong>cia etc.<br />

Ä Kinetic: Requiere <strong>los</strong> parÖmetros cináticos <strong>de</strong> la ecuaciÄn <strong>de</strong> Arrh<strong>en</strong>ius y la<br />

estequiometrÑa. Se pue<strong>de</strong> <strong>de</strong>finir el ord<strong>en</strong> <strong>de</strong> cada reacciÄn.<br />

Al <strong>de</strong>finir <strong>los</strong> coefici<strong>en</strong>tes estequiomátricos, HYSYS calcula el balance <strong>de</strong> moles e indica<br />

<strong>en</strong> el caso que <strong>los</strong> coefici<strong>en</strong>tes ingresados están mal balanceados.<br />

Las reacciones <strong>de</strong> combustiÄn <strong>de</strong> gas natural que se consi<strong>de</strong>raron son:<br />

CH î 2O ï 2H O î CO<br />

4 2 2 2<br />

2C H î 7O ï 6H O î 4CO<br />

2 6 2 2 2<br />

C H î 5O ï 4H O î 3CO<br />

3 8 2 2 2<br />

2C4H10 î 13O2 ï 10H 2O î 8CO2<br />

[Ec. 15]<br />

Se han consi<strong>de</strong>rado metano, etano, propano y butano, puesto que el gas que se quema<br />

<strong>en</strong> la turbina correspon<strong>de</strong> al gas residual <strong>de</strong> Planta PosesiÄn, y <strong>de</strong> hecho la reacciÄn <strong>de</strong> <strong>los</strong><br />

butanos no se lleva a cabo pues es casi seguro que nunca habrÖ butano <strong>en</strong> el gas residual.<br />

La forma <strong>de</strong> <strong>de</strong>finir estas reacciones se <strong>en</strong>cu<strong>en</strong>tra completam<strong>en</strong>te explicada <strong>en</strong> el<br />

manual que se adjunta.<br />

Los tipos <strong>de</strong> reactores que se pued<strong>en</strong> <strong>de</strong>finir <strong>en</strong> HYSYS correspond<strong>en</strong> a:<br />

Ä CSTR: Continuous Stirred Tank Reactor, reactor totalm<strong>en</strong>te agitado <strong>de</strong> flujo<br />

continuo.<br />

Ä PFR: Plug Flow Reactor, reactor <strong>de</strong> flujo pistÄn.


26<br />

Ä Conversion Reactor: Reactor <strong>de</strong> conversiÄn.<br />

Ä Equilibrium Reactor: Reactor <strong>de</strong> equilibrio.<br />

Ä Gibbs Reactor: Reactor <strong>de</strong> Gibbs.<br />

Para el caso que concierne a este trabajo, se utilizaron dos tipos <strong>de</strong> reactores; <strong>de</strong><br />

conversiÄn y <strong>de</strong> Gibbs.<br />

IlustraciÄn 10.- Tipos <strong>de</strong> reactores g<strong>en</strong>erales<br />

El reactor <strong>de</strong> conversiÄn <strong>de</strong>sarrolla la transformaciÄn <strong>de</strong> reactantes <strong>en</strong> productos <strong>de</strong><br />

acuerdo a la informaciÄn estequiomátrica ingresada, consi<strong>de</strong>rando el reactivo limitante y la<br />

fÄrmula <strong>de</strong> conversiÄn <strong>en</strong> funciÄn <strong>de</strong> la temperatura <strong>de</strong>l reactor, <strong>de</strong> acuerdo a:<br />

% Ä Ä<br />

2<br />

conv<br />

Å Co î C1 T î C2<br />

T<br />

[Ec. 16]<br />

Si se <strong>de</strong>sea una conversiÄn <strong>de</strong> 100% simplem<strong>en</strong>te se asigna 100 a la variable Co. Para<br />

mo<strong>de</strong><strong>los</strong> mÖs <strong>de</strong>tallados se <strong>de</strong>be t<strong>en</strong>er informaciÄn <strong>de</strong>l porc<strong>en</strong>taje <strong>de</strong> conversiÄn a<br />

<strong>de</strong>terminadas temperaturas y ajustar <strong>los</strong> datos experim<strong>en</strong>tales a una correlaciÄn polinomial<br />

cuadrÖtica para obt<strong>en</strong>er <strong>los</strong> parÖmetros C 0 , C 1 y C 2 .<br />

El reactor <strong>de</strong> Gibbs, por su parte, no necesita informaciÄn <strong>de</strong> reacciones ni coefici<strong>en</strong>tes<br />

estequiomátricos, ya que <strong>los</strong> resultados se obti<strong>en</strong><strong>en</strong> al especificar que la corri<strong>en</strong>te <strong>de</strong> salida<br />

<strong>de</strong>be poseer la mÑnima cantidad posible <strong>de</strong> <strong>en</strong>ergÑa libre <strong>de</strong> Gibbs, cuyo valor se calcula<br />

mediante el mo<strong>de</strong>lo termodinÖmico elegido.


27<br />

Ambos tipos <strong>de</strong> reactores <strong>en</strong>tregan la opciÄn <strong>de</strong> funcionar como un separador lÑquido<br />

vapor, pero a<strong>de</strong>mÖs, el reactor <strong>de</strong> Gibbs sirve como un reactor <strong>de</strong> equilibrio si se especifica un<br />

set <strong>de</strong> reacciones <strong>de</strong> equilibrio.<br />

Al funcionar <strong>en</strong> estado estacionario, las dim<strong>en</strong>siones fÑsicas que se pued<strong>en</strong> especificar<br />

para el reactor (altura, diÖmetro, nivel <strong>de</strong> lÑquido, etc.) no ti<strong>en</strong><strong>en</strong> importancia.<br />

2.4.4 COLUMNAS DE DESTILACIÅN<br />

Sin lugar a dudas esta operaciÄn unitaria es para la que HYSYS <strong>en</strong>trega la mayor<br />

cantidad <strong>de</strong> opciones, pues es un equipo fundam<strong>en</strong>tal <strong>en</strong> <strong>los</strong> <strong>procesos</strong> <strong>de</strong> la industria <strong>de</strong>l<br />

petrÄleo y gas natural. HYSYS incluye un ambi<strong>en</strong>te especÑfico, un sub-flowsheet in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te<br />

<strong>de</strong> la simulaciÄn global, <strong>en</strong> el que se pue<strong>de</strong> modificar la columna <strong>de</strong> <strong>de</strong>stilaciÄn ingresada.<br />

EL sub-flowsheet <strong>de</strong> la columna conti<strong>en</strong>e <strong>los</strong> equipos y corri<strong>en</strong>tes asociados, e<br />

intercambia informaciÄn con el flowsheet global, mediante las conexiones <strong>de</strong> corri<strong>en</strong>tes <strong>de</strong><br />

<strong>en</strong>trada y salida. Des<strong>de</strong> el punto <strong>de</strong> vista <strong>de</strong> la simulaciÄn global, la columna se muestra como<br />

una operaciÄn con mçltiples <strong>en</strong>tradas y salidas, corri<strong>en</strong>tes <strong>en</strong>ergáticas, etc.<br />

Haci<strong>en</strong>do doble clic <strong>en</strong> dicho Ñcono, se pue<strong>de</strong> ingresar al sub-flowsheet especÑfico <strong>de</strong> la<br />

columna, y al hacer dicho cambio la simulaciÄn global <strong>en</strong>tra a modo Holding, o sea, pausa el<br />

motor <strong>de</strong> cÖlculo hasta que se hayan realizado las modificaciones internas <strong>de</strong> la columna.<br />

Normalm<strong>en</strong>te, si la columna consta <strong>de</strong> un cond<strong>en</strong>sador <strong>de</strong> tope, una torre <strong>de</strong> platos y<br />

un rehervidor <strong>de</strong> fondo, no es imperativo ingresar al sub-flowsheet <strong>de</strong> la columna, y se pue<strong>de</strong><br />

especificar <strong>de</strong>s<strong>de</strong> la simulaciÄn global. Pero si el sistema repres<strong>en</strong>tado incluye extraciones<br />

laterales, recirculaciones, rectificadores, etc., es necesario ingresar a la columna y especificar<strong>los</strong><br />

<strong>en</strong> la misma metodologÑa que la simulaciÄn global; mediante Ñconos que repres<strong>en</strong>tan equipos y<br />

corri<strong>en</strong>tes, y lÑneas <strong>de</strong> conexiÄn.


28<br />

se pue<strong>de</strong> citar:<br />

Entre las principales v<strong>en</strong>tajas <strong>de</strong> consi<strong>de</strong>rar un flowsheet especÑfico para las columnas<br />

Ä In<strong>de</strong>p<strong>en</strong>dizar el mátodo <strong>de</strong> resoluciÄn<br />

Ä Uso opcional <strong>de</strong> una termodinÖmica difer<strong>en</strong>te a la global<br />

Ä ConstrucciÄn <strong>de</strong> columnas estÖndar (Templates) para ser ocupadas <strong>en</strong> otras<br />

simulaciones, como Åcaja negraÉ.<br />

Ä Se pued<strong>en</strong> resolver mçltiples columnas simultÖneam<strong>en</strong>te.<br />

La principal difer<strong>en</strong>cia <strong>en</strong>tre el flowsheet global y el <strong>de</strong> la columna se pue<strong>de</strong> notar a<br />

continuaciÄn, para la misma columna el aspecto es:<br />

IlustraciÄn 11.- Aspecto <strong>de</strong> una columna <strong>de</strong> <strong>de</strong>stilaciÄn <strong>en</strong> HYSYS<br />

Entre las unida<strong>de</strong>s mÖs complejas que HYSYS pue<strong>de</strong> simular se <strong>en</strong>cu<strong>en</strong>tran torres <strong>de</strong><br />

fraccionami<strong>en</strong>to, <strong>de</strong>stilaciÄn <strong>en</strong> vacÑo, columnas <strong>de</strong>etanizadoras, absorbedores y <strong>de</strong>stilaciÄn<br />

extractiva. Todos estos equipos consi<strong>de</strong>ran una serie <strong>de</strong> etapas <strong>de</strong> equilibrio <strong>en</strong>tre una corri<strong>en</strong>te<br />

<strong>de</strong> vapor que asci<strong>en</strong><strong>de</strong> y un lÑquido que <strong>de</strong>sci<strong>en</strong><strong>de</strong>, a<strong>de</strong>mÖs <strong>de</strong> consi<strong>de</strong>rar mçltiples<br />

extracciones, alim<strong>en</strong>taciones, recirculaciones, etc.


29<br />

IlustraciÄn 12.- Esquema <strong>de</strong> un proceso <strong>de</strong> separaciÄn por etapas <strong>de</strong> equilibrio<br />

2.4.4.1 ECUACIONES DE UNA COLUMNA DE DESTILACIÅN<br />

HYSYS aplica balances <strong>de</strong> masa y <strong>en</strong>ergÑa a cada plato <strong>de</strong> la columna <strong>en</strong> cuestiÄn.<br />

Consi<strong>de</strong>rando la nom<strong>en</strong>clatura <strong>de</strong> la IlustraciÄn 12, las ecuaciones son:<br />

Balance <strong>de</strong> masa:<br />

F î L îV Å L îV î R îVSD î LSD<br />

[Ec. 17]<br />

Global:<br />

j jÉ1 jî1<br />

j j j j j<br />

Compon<strong>en</strong>te: FjÄz j<br />

î LjÉ 1Äx jÉ1 îVjî 1Ä y<br />

jî1<br />

Å ( Lj î LSDj )Ä x<br />

j<br />

î ( V<br />

j<br />

îVSDj )Äy j<br />

î RjÄz<br />

j<br />

[Ec. 18]<br />

Balance <strong>de</strong> <strong>en</strong>ergÑa, con H <strong>en</strong>talpÑa <strong>de</strong>l vapor y h <strong>en</strong>talpÑa <strong>de</strong>l lÑquido:<br />

FjÄH Fj<br />

î LjÉ 1ÄhjÉ1 îVjî 1Ä H<br />

jî1<br />

î Qj Å ( Lj î LSDj )Ä hj î ( Vj îVSDj )ÄH j<br />

î RjÄH<br />

j<br />

[Ec. 19]<br />

Y la forma <strong>en</strong> que se relacionan las composiciones <strong>de</strong>l lÑquido y <strong>de</strong>l vapor que sal<strong>en</strong> <strong>de</strong><br />

cada plato estÖ dada por la sigui<strong>en</strong>te ecuaciÄn <strong>de</strong> equilibrio:


30<br />

yi Å kiÄ xi<br />

Para cada compon<strong>en</strong>te " i"<br />

[Ec. 20]<br />

En don<strong>de</strong> el parÖmetro k i se calcula mediante el paquete termodinÖmico seleccionado.<br />

Estas ecuaciones explican el comportami<strong>en</strong>to <strong>de</strong> <strong>los</strong> fluidos <strong>en</strong> cada plato <strong>de</strong> la<br />

columna. Si bi<strong>en</strong> se pued<strong>en</strong> especificar las condiciones especÑficas <strong>de</strong> un plato <strong>de</strong>terminado, <strong>en</strong><br />

HYSYS interesa mÖs <strong>de</strong>finir las variables globales (flujos <strong>de</strong> salida, recirculaciones, composiciones<br />

<strong>de</strong> productos, etc). Las ecuaciones globales se aplican a la columna vista Åpor fueraÉ y <strong>de</strong><br />

acuerdo a la IlustraciÄn 13, se pued<strong>en</strong> <strong>de</strong>ducir las sigui<strong>en</strong>tes expresiones:<br />

Balance <strong>de</strong> masa aplicado a toda la columna:<br />

Global: F Å D î W [Ec. 21]<br />

Compon<strong>en</strong>te: FÄZ Å DÄZ îW ÄX<br />

F D W<br />

[Ec. 22]<br />

Balance <strong>de</strong> <strong>en</strong>ergÑa, omiti<strong>en</strong>do <strong>los</strong> valores <strong>de</strong> párdida Q Ln y Q Lm :<br />

FÄH î Q Å Q î DÄH îWÄH<br />

F W C F W<br />

[Ec. 23]<br />

Se <strong>de</strong>fine la razÄn <strong>de</strong> reflujo (R Ratio ) <strong>en</strong> el tope como el coci<strong>en</strong>te <strong>en</strong>tre el reflujo <strong>de</strong> tope<br />

y el flujo <strong>de</strong> <strong>de</strong>stilado:<br />

R<br />

Ratio<br />

Å<br />

L0<br />

D [Ec. 24]


31<br />

IlustraciÄn 13.- Esquema teÄrico <strong>de</strong> una columna <strong>de</strong> <strong>de</strong>stilaciÄn<br />

No se <strong>de</strong>be confundir la razÄn <strong>de</strong> reflujo (Reflux Ratio <strong>en</strong> HYSYS) con el flujo <strong>de</strong><br />

recirculaciÄn (Reflux Rate, <strong>en</strong> HYSYS). El segundo correspon<strong>de</strong> al flujo (mÖsico, molar o<br />

volumátrico) <strong>de</strong> la corri<strong>en</strong>te L 0 .


32<br />

2.4.4.2 ESPECIFICACIONES DE UNA COLUMNA DE DESTILACIÅN EN HYSYS<br />

Al aÜadir una columna <strong>de</strong> <strong>de</strong>stilaciÄn <strong>en</strong> HYSYS, luego <strong>de</strong> <strong>de</strong>finir las corri<strong>en</strong>tes<br />

<strong>en</strong>ergáticas y materiales asociadas, <strong>de</strong> no realizar modificaciones internas, las variables que se<br />

pued<strong>en</strong> especificar por <strong>de</strong>fecto son:<br />

Ä Overhead Vapour Flowrate: Correspon<strong>de</strong> al flujo <strong>de</strong> vapores <strong>de</strong> tope <strong>de</strong>l<br />

cond<strong>en</strong>sador parcial.<br />

Ä Distillate Flowrate: Flujo <strong>de</strong> <strong>de</strong>stilado que se <strong>de</strong>sea obt<strong>en</strong>er.<br />

Ä Bottoms Flowrate: Flujo <strong>de</strong> producto <strong>de</strong> fondo.<br />

Ä Reflux Ratio: RazÄn <strong>de</strong> reflujo <strong>de</strong> tope.<br />

Ä Reflux Rate: Flujo (molar, mÖsico o volumátrico) <strong>de</strong>l reflujo <strong>de</strong> tope.<br />

Sin embargo se pued<strong>en</strong> omitir estas especificaciones y agregar las que el usuario estime<br />

conv<strong>en</strong>i<strong>en</strong>te, ya sea temperatura <strong>en</strong> un plato especÑfico, fracciÄn molar <strong>de</strong> un compon<strong>en</strong>te <strong>en</strong><br />

un plato o corri<strong>en</strong>te, flujo <strong>de</strong> calor <strong>en</strong> el rehervidor o cond<strong>en</strong>sador, etc.<br />

HYSYS indica <strong>los</strong> grados <strong>de</strong> libertad que posee el sistema diseÜado, y se <strong>de</strong>b<strong>en</strong><br />

especificar tantas variables como grados <strong>de</strong> libertad, para po<strong>de</strong>r empezar el proceso iterativo<br />

que da lugar a la soluciÄn.<br />

Es importante escoger el mátodo iterativo que ocuparÖ HYSYS <strong>en</strong> la resoluciÄn <strong>de</strong> la<br />

columna. Se pue<strong>de</strong> escoger <strong>en</strong>tre seis mátodos, <strong>los</strong> cuales se <strong>de</strong>tallan a continuaciÄn:


33<br />

MÖtodo<br />

DescripciÉn<br />

HYSIM Insi<strong>de</strong>-Out Mátodo g<strong>en</strong>eral, çtil para la mayorÑa <strong>de</strong> <strong>los</strong> problemas. Es el que vi<strong>en</strong>e activado<br />

por <strong>de</strong>fecto al agregar una columna <strong>de</strong> <strong>de</strong>stilaciÄn.<br />

Modified HYSIM Mátodo g<strong>en</strong>eral, que permite agregar mezcladores, divisores, intercambiadores<br />

Insi<strong>de</strong>-Out<br />

<strong>de</strong> calor o separadores al sub-flowsheet <strong>de</strong> la columna.<br />

Newton Raphson Permite <strong>de</strong>finir reacciones cináticas <strong>en</strong> la fase lÑquida.<br />

Insi<strong>de</strong>-Out<br />

Sparse Continuation Soporta dos fases lÑquidas <strong>en</strong> la columna, se utiliza principalm<strong>en</strong>te para resolver<br />

Solver<br />

sistemas quÑmicos altam<strong>en</strong>te no i<strong>de</strong>ales, y <strong>de</strong>stilaciÄn reactiva.<br />

Simultaneous Similar al mátodo Sparse, pero no permite la adiciÄn <strong>de</strong> rectificadores laterales o<br />

Correction<br />

pump around.<br />

OLI Solver<br />

Se <strong>de</strong>be usar solam<strong>en</strong>te para cÖlcu<strong>los</strong> <strong>en</strong> sistemas electrolÑticos.<br />

Tabla 1.- Mátodos <strong>de</strong> resoluciÄn <strong>de</strong> columnas <strong>de</strong> <strong>de</strong>stilaciÄn <strong>en</strong> HYSYS<br />

Para <strong>los</strong> casos que se estudiarÖn, sirve el primer mátodo, y si se <strong>de</strong>b<strong>en</strong> realizar<br />

modificaciones <strong>en</strong> el esquema que <strong>en</strong>trega HYSYS inicialm<strong>en</strong>te, se <strong>de</strong>be escoger Modified<br />

HYSIM Insi<strong>de</strong>-Out.<br />

Lam<strong>en</strong>tablem<strong>en</strong>te, al ser HYSYS un software comercial cuyo cÄdigo fu<strong>en</strong>te se manti<strong>en</strong>e<br />

<strong>en</strong> estricto secreto, no se ti<strong>en</strong>e acceso al algoritmo resolutivo especÑfico <strong>de</strong> cada mátodo; la<br />

çnica informaciÄn que se ha podido obt<strong>en</strong>er es la Tabla 1, por otra parte <strong>en</strong> la pÖgina web <strong>de</strong>l<br />

producto se indica que el mátodo HYSIM Insi<strong>de</strong>-Out es el especificado por Russell (Ver<br />

bibliografÑa). Pero las modificaciones que Asp<strong>en</strong> Tech realizÄ sobre este mátodo, para obt<strong>en</strong>er<br />

<strong>los</strong> restantes, no estÖn disponibles al pçblico.<br />

2.4.5 AJUSTADORES<br />

Un ajustador <strong>en</strong> HYSYS es una operaciÄn que ajusta el valor <strong>de</strong> una variable especÑfica<br />

(variable in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te) a fin <strong>de</strong> obt<strong>en</strong>er un valor especÑfico <strong>en</strong> otra variable o operaciÄn<br />

(variable <strong>de</strong>p<strong>en</strong>di<strong>en</strong>te). Es una herrami<strong>en</strong>ta fundam<strong>en</strong>tal ya que automatiza el proceso <strong>de</strong><br />

prueba y error para obt<strong>en</strong>er cierto valor requerido.<br />

BÖsicam<strong>en</strong>te el ajustador pue<strong>de</strong> cumplir dos funciones:


34<br />

Ä Ajustar la variable in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te hasta lograr que la variable <strong>de</strong>p<strong>en</strong>di<strong>en</strong>te llegue<br />

a un valor especificado.<br />

Ä Ajustar la variable in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te hasta que la variable <strong>de</strong>p<strong>en</strong>di<strong>en</strong>te se iguale al<br />

valor <strong>de</strong> otro objeto <strong>de</strong>terminado.<br />

Aparte <strong>de</strong> <strong>de</strong>finir las variables in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te y <strong>de</strong>p<strong>en</strong>di<strong>en</strong>te, y el valor objetivo, se<br />

<strong>de</strong>b<strong>en</strong> especificar ciertos parÖmetros que son fundam<strong>en</strong>tales a la hora <strong>de</strong> realizar el<br />

procedimi<strong>en</strong>to automatizado <strong>de</strong> prueba y error, y son:<br />

Ä Method: Se ofrec<strong>en</strong> dos mátodos <strong>de</strong> resoluciÄn; Secant (no tan rÖpido<br />

comparativam<strong>en</strong>te, pero seguro) y Broyd<strong>en</strong> (rÖpido pero no tan estable)<br />

Ä Tolerance: La tolerancia indica quá error se pue<strong>de</strong> aceptar <strong>en</strong> la variable objetivo.<br />

Este error se calcula como absoluto, no es relativo o porc<strong>en</strong>tual. La expresiÄn<br />

es:<br />

Error Å Valor obt<strong>en</strong>ido ÉValor objetivo<br />

Ä Step Size: Correspon<strong>de</strong> a la mÖxima difer<strong>en</strong>cia que la variable in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te<br />

podrÑa variar <strong>en</strong>tre un int<strong>en</strong>to y otro. Este valor se utiliza hasta que la soluciÄn<br />

haya sido id<strong>en</strong>tificada <strong>en</strong> un rango, y <strong>en</strong>tonces se ocupa un algoritmo <strong>de</strong><br />

converg<strong>en</strong>cia especificado. Un valor positivo inicialm<strong>en</strong>te increm<strong>en</strong>ta el valor<br />

inicial <strong>de</strong> la variable in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te, mi<strong>en</strong>tras que un valor negativo la<br />

disminuye. Si el valor se aleja <strong>de</strong>l rango <strong>de</strong> soluciÄn, la direcciÄn <strong>de</strong> las<br />

iteraciones se revierte automÖticam<strong>en</strong>te.<br />

Ä MÜximum/Minimum: Se pue<strong>de</strong> acotar la variable in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te a un rango, si se<br />

sabe que el resultado <strong>de</strong>be estar <strong>en</strong> tal rango, con lo que se acelera el proceso<br />

<strong>de</strong> converg<strong>en</strong>cia. Con esto tambián se pue<strong>de</strong> evitar que el ajustador <strong>en</strong>tregue<br />

respuestas incoher<strong>en</strong>tes, como flujos negativos.<br />

Ä Maximum Iterations: Por <strong>de</strong>fecto el ajustador realiza 10 iteraciones, pero se<br />

pue<strong>de</strong> especificar cualquier valor.


35<br />

IlustraciÄn 14.- Aspecto <strong>de</strong>l ajustador <strong>en</strong> HYSYS<br />

Tambián se ofrece la opciÄn <strong>de</strong> utilizar mçltiples ajustadores y resolver<br />

simultÖneam<strong>en</strong>te el sistema, y <strong>en</strong> tal caso se <strong>de</strong>sactiva la opciÄn Method, y se utiliza solam<strong>en</strong>te<br />

el algoritmo <strong>de</strong> Lev<strong>en</strong>berg y Marquardt, modificado por <strong>los</strong> fabricantes <strong>de</strong>l programa.<br />

2.4.6 RECICLOS<br />

La capacidad <strong>de</strong> cualquier simulador <strong>de</strong> <strong>procesos</strong> quÑmicos <strong>de</strong> repres<strong>en</strong>tar sistemas con<br />

recic<strong>los</strong> resulta <strong>de</strong> un valor fundam<strong>en</strong>tal. HYSYS pres<strong>en</strong>ta un esquema especÑfico para resolver<br />

recic<strong>los</strong> que estÖ evaluado como uno <strong>de</strong> <strong>los</strong> mejores <strong>de</strong>l Ömbito.<br />

Utilizar un reciclo implica instalar un mÄdulo teÄrico <strong>en</strong>tre las corri<strong>en</strong>tes <strong>de</strong> proceso,<br />

que realizarÖ un procedimi<strong>en</strong>to iterativo especÑfico. La gran versatilidad <strong>de</strong> dicho mÄdulo radica<br />

<strong>en</strong> que las condiciones <strong>de</strong> proceso pued<strong>en</strong> ser transferidas hacia atrÖs o hacia <strong>de</strong>lante (<strong>en</strong> el<br />

s<strong>en</strong>tido <strong>de</strong>l flujo) <strong>en</strong>tre la <strong>en</strong>trada y salida <strong>de</strong>l reciclo. En tárminos <strong>de</strong> la soluciÄn que se busca,<br />

hay valores estimados y valores calculados, para cada una <strong>de</strong> las variables <strong>en</strong> las corri<strong>en</strong>tes <strong>de</strong><br />

<strong>en</strong>trada y salida. Dep<strong>en</strong>di<strong>en</strong>do <strong>de</strong> la direcciÄn <strong>en</strong> la que se transfier<strong>en</strong> <strong>los</strong> datos, <strong>los</strong> valores<br />

estimados pued<strong>en</strong> existir tanto <strong>en</strong> la <strong>en</strong>trada como <strong>en</strong> la salida. Por ejemplo, si el usuario escoge<br />

como direcciÄn <strong>de</strong> transfer<strong>en</strong>cia Backwards (hacia atrÖs) para la temperatura, el valor estimado


36<br />

es la temperatura <strong>de</strong> la corri<strong>en</strong>te <strong>de</strong> <strong>en</strong>trada, y el calculado es la temperatura <strong>de</strong> la corri<strong>en</strong>te <strong>de</strong><br />

salida.<br />

Durante el proceso <strong>de</strong> converg<strong>en</strong>cia se suced<strong>en</strong> las sigui<strong>en</strong>tes etapas <strong>de</strong> cÖlculo,<br />

aproximadam<strong>en</strong>te:<br />

Ä HYSYS utiliza <strong>los</strong> valores estimados y resuelve el flowsheet.<br />

Ä Luego HYSYS compara <strong>los</strong> valores estimados con <strong>los</strong> que se obti<strong>en</strong><strong>en</strong> por<br />

resoluciÄn <strong>de</strong>l flowsheet <strong>en</strong> la otra corri<strong>en</strong>te relacionada <strong>en</strong> el reciclo<br />

Ä BasÖndose <strong>en</strong> la difer<strong>en</strong>cia <strong>en</strong>tre <strong>los</strong> valores calculados y estimados, HYSYS<br />

g<strong>en</strong>era nuevos valores que se toman como nuevas estimaciones.<br />

Ä El proceso se repite hasta que <strong>los</strong> valores obt<strong>en</strong>idos <strong>en</strong> las dos corri<strong>en</strong>tes<br />

relacionadas <strong>en</strong> el reciclo difier<strong>en</strong> <strong>en</strong> un valor especificado lo sufici<strong>en</strong>tem<strong>en</strong>te<br />

bajo para asegurar la converg<strong>en</strong>cia.<br />

IlustraciÄn 15.- Aspecto <strong>de</strong>l mÄdulo <strong>de</strong> reciclo <strong>en</strong> HYSYS<br />

De todas las opciones, bÖsicam<strong>en</strong>te numáricas, la que mÖs interesa es el mo<strong>de</strong>lo <strong>de</strong><br />

resoluciÄn, que pue<strong>de</strong> elegirse <strong>en</strong>tre Nested y Simultaneous. Como su nombre lo indica, <strong>en</strong><br />

caso que la simulaciÄn pres<strong>en</strong>te difer<strong>en</strong>tes recic<strong>los</strong> interconectados o relacionados <strong>de</strong> alguna<br />

manera, es necesario utilizar la opciÄn Simultaneous. Para simulaciones con un solo reciclo, o<br />

varios recic<strong>los</strong> cuyos resultados son in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>tes, se pue<strong>de</strong> utilizar la opciÄn Nested.<br />

Nuevam<strong>en</strong>te, el <strong>de</strong>sarrollo interno <strong>de</strong> estos algoritmos es informaciÄn que se manti<strong>en</strong>e<br />

<strong>en</strong> estricto secreto, razÄn por la cual no se pue<strong>de</strong> dar una mayor exposiciÄn <strong>de</strong> cÄmo HYSYS<br />

realiza cÖlcu<strong>los</strong> <strong>de</strong> recic<strong>los</strong>.


CAPÄTULO III<br />

PARTE EXPERIMENTAL Y RESULTADOS


38<br />

III - PARTE EXPERIMENTAL Y RESULTADOS<br />

En el pres<strong>en</strong>te capÑtulo se <strong>de</strong>sarrollan las consi<strong>de</strong>raciones y metodologÑa utilizada <strong>en</strong> el<br />

<strong>de</strong>sarrollo <strong>de</strong> las simulaciones HYSYS <strong>de</strong> <strong>los</strong> <strong>procesos</strong> <strong>productivos</strong> seleccionados como mÖs<br />

caracterÑsticos d<strong>en</strong>tro <strong>de</strong> la actividad <strong>de</strong> ENAP Magallanes, se <strong>en</strong>umeran y explican las<br />

simulaciones creadas, se analizan <strong>los</strong> resultados obt<strong>en</strong>idos y se contrastan (<strong>en</strong> <strong>los</strong> casos<br />

posibles) con la situaciÄn real.<br />

3.1 ESQUEMA DE DESARROLLO<br />

El trabajo se ha <strong>de</strong>sarrollado parti<strong>en</strong>do <strong>de</strong> la premisa que todos <strong>los</strong> pasos necesarios<br />

para crear las simulaciones puedan ser recreados por <strong>los</strong> operadores e ing<strong>en</strong>ieros <strong>de</strong> <strong>procesos</strong><br />

<strong>de</strong> ENAP Magallanes, por lo tanto <strong>en</strong> el ÅManual para simular <strong>procesos</strong> <strong>productivos</strong> <strong>de</strong> ENAP<br />

Magallanes <strong>en</strong> HYSYSÉ se <strong>de</strong>b<strong>en</strong> agregar <strong>de</strong>s<strong>de</strong> <strong>los</strong> conocimi<strong>en</strong>tos bÖsicos <strong>de</strong> uso <strong>de</strong> HYSYS,<br />

completam<strong>en</strong>te <strong>de</strong>tallados, hasta las simulaciones mÖs ext<strong>en</strong>sas, <strong>de</strong> planta PosesiÄn y Cull<strong>en</strong>.<br />

Por lo tanto se ha dividido el trabajo <strong>en</strong> las sigui<strong>en</strong>tes etapas <strong>de</strong> <strong>de</strong>sarrollo:<br />

Ä Dar a conocer <strong>los</strong> aspectos bÖsicos y g<strong>en</strong>eralida<strong>de</strong>s <strong>de</strong> las simulaciones <strong>de</strong> HYSYS.<br />

Ä Desarrollar aplicaciones rÖpidas y g<strong>en</strong>erales que utilic<strong>en</strong> las unida<strong>de</strong>s bÖsicas <strong>de</strong><br />

proceso <strong>de</strong> ENAP Magallanes.<br />

Ä Profundizar <strong>en</strong> las herrami<strong>en</strong>tas es<strong>en</strong>ciales <strong>de</strong> HYSYS que no ti<strong>en</strong><strong>en</strong> equival<strong>en</strong>cia<br />

con <strong>los</strong> equipos instalados <strong>en</strong> terr<strong>en</strong>o (ajustadores, planillas <strong>de</strong> cÖlculo, etc.)<br />

Ä Desarrollar casos reales <strong>de</strong> ENAP Magallanes.<br />

Ä Analizar y ajustar las simulaciones para lograr coher<strong>en</strong>cia <strong>en</strong>tre el mo<strong>de</strong>lo y la<br />

situaciÄn real.<br />

Se ha d<strong>en</strong>ominado unida<strong>de</strong>s bÖsicas <strong>de</strong> proceso a todos <strong>los</strong> equipos que sirv<strong>en</strong> para<br />

llevar a cabo <strong>los</strong> <strong>procesos</strong> <strong>de</strong> ENAP Magallanes, que estÖn pres<strong>en</strong>tes <strong>en</strong> todas las instalaciones<br />

productivas <strong>de</strong> la empresa y que <strong>en</strong> conjunto conforman lo que se d<strong>en</strong>omina ÅplantaÉ. Son


39<br />

operaciones unitarias que realizan una acciÄn <strong>de</strong>terminada y por lo g<strong>en</strong>eral consi<strong>de</strong>ran una<br />

corri<strong>en</strong>te <strong>de</strong> <strong>en</strong>trada, una corri<strong>en</strong>te <strong>de</strong> salida y la <strong>en</strong>ergÑa asociada.<br />

En HYSYS estas unida<strong>de</strong>s bÖsicas correspond<strong>en</strong> a <strong>los</strong> objetos que se <strong>de</strong>b<strong>en</strong> aÜadir a las<br />

simulaciones, para con las interconexiones a<strong>de</strong>cuadas, obt<strong>en</strong>idas <strong>de</strong>l flowsheet <strong>de</strong> planta<br />

obt<strong>en</strong>er la simulaciÄn final <strong>de</strong>l proceso investigado.<br />

Las difer<strong>en</strong>tes unida<strong>de</strong>s bÖsicas <strong>de</strong> proceso que se han consi<strong>de</strong>rado como necesarias<br />

para el <strong>de</strong>sarrollo <strong>de</strong> las difer<strong>en</strong>tes simulaciones son:<br />

Ä Separadores <strong>de</strong> lÑquido y vapor.<br />

Ä Intercambiadores <strong>de</strong> calor:<br />

o Coolers y Heaters.<br />

o Intercambiadores <strong>de</strong> tubo y carcasa.<br />

Ä Mezcladores y divisores <strong>de</strong> flujo.<br />

Ä VÖlvulas.<br />

Ä Bombas.<br />

Ä Compresores.<br />

o Compresores c<strong>en</strong>trÑfugos.<br />

o Compresores recÑprocos.<br />

Estas unida<strong>de</strong>s <strong>en</strong> HYSYS pres<strong>en</strong>tan difer<strong>en</strong>tes opciones <strong>de</strong> especificaciÄn; <strong>en</strong> primera<br />

instancia se han <strong>de</strong>sarrollado aplicaciones que sÄlo necesitan una <strong>de</strong>finiciÄn rÖpida y<br />

aproximada, para posteriorm<strong>en</strong>te <strong>de</strong>sarrollar simulaciones <strong>de</strong> casos reales cuyo principal<br />

objetivo es repres<strong>en</strong>tar <strong>de</strong> la forma mÖs exacta posible estas unida<strong>de</strong>s bÖsicas, y esto<br />

necesariam<strong>en</strong>te implica utilizar las opciones avanzadas <strong>de</strong> las unida<strong>de</strong>s investigadas.<br />

En una segunda etapa se <strong>de</strong>sarrollan ejemp<strong>los</strong> para lograr un domino <strong>de</strong> las<br />

herrami<strong>en</strong>tas especÑficas <strong>de</strong> HYSYS, y para esto se utilizan las sigui<strong>en</strong>tes aplicaciones:


40<br />

Ä Ajustadores: Se <strong>de</strong>sarrolla la simulaciÄn <strong>de</strong> la unidad reg<strong>en</strong>eradora <strong>de</strong> etil<strong>en</strong>glicol<br />

<strong>de</strong> Planta Cull<strong>en</strong>.<br />

Ä Set: Se <strong>de</strong>sarrolla un esquema <strong>de</strong> expansor y compresor acoplados, simulando el<br />

sistema utilizado <strong>en</strong> Planta PosesiÄn<br />

Ä Planilla <strong>de</strong> cÜlculo: Se <strong>de</strong>sarrolla el proceso <strong>de</strong> flasheo <strong>de</strong> gasolinas <strong>de</strong> Planta<br />

Cull<strong>en</strong> y mediante el Spreadsheet (planilla <strong>de</strong> cÖlculo) se cuantifica la ganancia<br />

obt<strong>en</strong>ida por v<strong>en</strong>ta <strong>de</strong> Raw Product.


41<br />

3.2 SIMULACIONES REALIZADAS PARA ENAP MAGALLANES<br />

Las simulaciones <strong>de</strong> casos especÑficos <strong>de</strong> ENAP Magallanes que se han realizado<br />

durante el pres<strong>en</strong>te trabajo <strong>de</strong> titulaciÄn son:<br />

Ä Compresor recÑproco HRA-5, <strong>de</strong> Planta Cull<strong>en</strong>.<br />

Ä Compresor recÑproco HRA-2, <strong>de</strong> EstaciÄn Compresora Calafate.<br />

Ä Compresor recÑproco TLA-3, <strong>de</strong> Planta PosesiÄn.<br />

Ä Compresor c<strong>en</strong>trÑfugo TC-10, <strong>de</strong> Planta Cull<strong>en</strong>.<br />

Ä SimulaciÄn completa <strong>de</strong> Planta Cull<strong>en</strong>, que consta <strong>de</strong>:<br />

o Circuito <strong>de</strong> refrigeraciÄn por propano.<br />

o Sistema <strong>de</strong> acondicionami<strong>en</strong>to y <strong>de</strong>shidrataciÄn <strong>de</strong> gas.<br />

o Flasheo <strong>de</strong> gasolinas.<br />

o Unidad reg<strong>en</strong>eradora <strong>de</strong> etil<strong>en</strong>glicol.<br />

Ä Turbina <strong>de</strong>l tr<strong>en</strong> SULZER, <strong>de</strong> Planta PosesiÄn.<br />

Ä Esquema <strong>de</strong> producciÄn <strong>de</strong> Plantas PosesiÄn y Cabo Negro.<br />

Ä Ejemplo <strong>de</strong> caracterizaciÄn <strong>de</strong> petrÄleos <strong>en</strong> HYSYS.<br />

Ä Ejemplo <strong>de</strong> loop <strong>de</strong> gasoductos.<br />

Ä SimulaciÄn completa <strong>de</strong> Planta PosesiÄn.<br />

A<strong>de</strong>mÖs se han <strong>de</strong>sarrollado dos simulaciones basadas <strong>en</strong> materiales <strong>de</strong> <strong>en</strong>tr<strong>en</strong>ami<strong>en</strong>to<br />

que Asp<strong>en</strong> ofrece a <strong>los</strong> usuarios registrados, <strong>en</strong> la secciÄn <strong>de</strong> ÅAdvanced Process Mo<strong>de</strong>ling Using<br />

Asp<strong>en</strong> HYSYSÉ <strong>de</strong> su pÖgina web, las cuales son:<br />

Ä Gas Gathering; repres<strong>en</strong>taciÄn <strong>de</strong> una red <strong>de</strong> gasoductos.<br />

Ä Mo<strong>de</strong>ling Real Separators in Asp<strong>en</strong> HYSYS; esquema <strong>de</strong> un separador lÑquido<br />

vapor que permite simular el efecto <strong>de</strong> arrastre <strong>de</strong> lÑquido <strong>en</strong> la salida <strong>de</strong> vapor.<br />

Tambián es <strong>de</strong>l interás <strong>de</strong> la empresa explicar el funcionami<strong>en</strong>to <strong>de</strong> un software<br />

adicional asociado a HYSYS, llamado HYSYS Browser, disponible para su <strong>de</strong>scarga <strong>de</strong>s<strong>de</strong> el sitio<br />

web <strong>de</strong>l software, para usuarios registrados.


42<br />

HYSYS Browser es un complem<strong>en</strong>to a Microsoft Excel, y permite la intercomunicaciÄn<br />

<strong>en</strong>tre Microsoft Excel y HYSYS. Esto significa que <strong>los</strong> datos y especificaciones que se <strong>de</strong>see<br />

ingresar a HYSYS pued<strong>en</strong> prov<strong>en</strong>ir <strong>de</strong> una planilla <strong>de</strong> cÖlculo <strong>de</strong> Microsoft Excel, y viceversa, que<br />

<strong>los</strong> resultados obt<strong>en</strong>idos <strong>en</strong> la simulaciÄn <strong>de</strong> HYSYS se puedan exportar a Microsoft Excel; con lo<br />

cual, relacionando a<strong>de</strong>cuadam<strong>en</strong>te estos dos programas mediante el uso <strong>de</strong> las opciones que<br />

<strong>en</strong>trega HYSYS Browser, ev<strong>en</strong>tualm<strong>en</strong>te se lograrÑa trabajar con HYSYS manejando cualquier<br />

aspecto <strong>en</strong> el ambi<strong>en</strong>te grÖfico <strong>de</strong> Microsoft Excel.<br />

Se ha <strong>de</strong>sarrollado un ejemplo que muestra todas las posibilida<strong>de</strong>s que ofrece HYSYS<br />

Browser, basado <strong>en</strong> un caso ficticio <strong>de</strong> un compresor <strong>de</strong> dos etapas.<br />

Todas estas simulaciones se pued<strong>en</strong> recrear por completo al <strong>de</strong>sarrollar el ÅManual<br />

para simular <strong>procesos</strong> <strong>productivos</strong> <strong>de</strong> ENAP Magallanes <strong>en</strong> HYSYSÉ, que se anexa a este trabajo.


43<br />

3.3 CONSIDERACIONES ESPECÄFICAS DE CADA SIMULACIÅN<br />

Es importante recordar que una simulaciÄn <strong>de</strong> HYSYS no repres<strong>en</strong>ta el ci<strong>en</strong> por ci<strong>en</strong>to<br />

<strong>de</strong> <strong>los</strong> equipos, conexiones o corri<strong>en</strong>tes que <strong>en</strong> planta se pued<strong>en</strong> observar. Debido a esto, se<br />

<strong>de</strong>b<strong>en</strong> realizar ciertas simplificaciones, que conllevan a trabajar con un flowsheet aproximado al<br />

esquema real, pero que <strong>de</strong>p<strong>en</strong>di<strong>en</strong>do <strong>de</strong>l criterio <strong>de</strong> qui<strong>en</strong> diseÜe las simulaciones podrÖ ser tan<br />

exacto como se <strong>de</strong>see.<br />

Para id<strong>en</strong>tificar la informaciÄn disponible relativa al caso estudiado, y po<strong>de</strong>r juzgar<br />

a<strong>de</strong>cuadam<strong>en</strong>te el mo<strong>de</strong>lo diseÜado, se pres<strong>en</strong>ta a continuaciÄn una breve <strong>de</strong>scripciÄn <strong>de</strong> las<br />

simulaciones, las limitaciones intrÑnsecas al trabajo <strong>de</strong> mo<strong>de</strong>lar un proceso y <strong>los</strong> datos que se<br />

han consi<strong>de</strong>rado para la posterior validaciÄn:<br />

3.3.1 COMPRESOR HRA-5 DE PLANTA CULLEN<br />

El compresor HRA-5 estÖ instalado <strong>en</strong> la sala <strong>de</strong> compresores <strong>de</strong> Planta Cull<strong>en</strong>. Es un<br />

compresor recÑproco con 4 cilindros compresores, todos <strong>de</strong> igual diÖmetro <strong>de</strong> pistÄn y vÖstago.<br />

La simulaciÄn se ha realizado ingresando <strong>los</strong> parÖmetros fÑsicos <strong>en</strong>tregados por el<br />

Informe <strong>de</strong> Mant<strong>en</strong>imi<strong>en</strong>to Predictivo. SÄlo se han especificado las variables necesarias para<br />

resolver el sistema, y <strong>los</strong> resultados han sido contrastados con el resto <strong>de</strong> informaciÄn<br />

prov<strong>en</strong>i<strong>en</strong>te <strong>de</strong>l Informe <strong>de</strong> Mant<strong>en</strong>imi<strong>en</strong>to.<br />

Se ha consi<strong>de</strong>rado que es necesario esquematizar el separador <strong>de</strong> lÑquido y vapor a la<br />

<strong>en</strong>trada <strong>de</strong>l compresor, pues aunque al ingresar las propieda<strong>de</strong>s <strong>de</strong> la corri<strong>en</strong>te <strong>de</strong> <strong>en</strong>trada se<br />

ha observado que dicho flujo no pres<strong>en</strong>ta lÑquidos (fracciÄn <strong>de</strong> vapor igual a 1), aÜadir el<br />

separador es un proceso rÖpido, y da mayor amplitud <strong>de</strong> estudio al sistema, preparando la<br />

simulaciÄn para un ev<strong>en</strong>tual caso <strong>en</strong> que la corri<strong>en</strong>te <strong>de</strong> <strong>en</strong>trada cont<strong>en</strong>ga lÑquidos.<br />

Se han ingresado a<strong>de</strong>mÖs <strong>los</strong> parÖmetros fÑsicos <strong>de</strong> cada cilindro (diÖmetro <strong>de</strong> pistones,<br />

carrera, clearances, etc.), lo cual da mayor similitud al mo<strong>de</strong>lo.


44<br />

3.3.2 COMPRESOR HRA-2 DE ESTACIÅN COMPRESORA CALAFATE<br />

Este compresor es similar al HRA-5 ya m<strong>en</strong>cionado, con la çnica salvedad que pres<strong>en</strong>ta<br />

difer<strong>en</strong>tes caracterÑsticas para <strong>los</strong> cuatro cilindros compresores, y dichas difer<strong>en</strong>cias se <strong>de</strong>b<strong>en</strong><br />

repres<strong>en</strong>tar <strong>en</strong> la simulaciÄn como dos compresores individuales, cada cual con caracterÑsticas<br />

<strong>de</strong>finidas.<br />

Para simular este sistema se <strong>de</strong>be agregar un divisor, y HYSYS calcularÖ<br />

automÖticam<strong>en</strong>te <strong>los</strong> porc<strong>en</strong>tajes <strong>de</strong> divisiÄn <strong>de</strong>bido a que el flujo que pasa por cada compresor<br />

es funciÄn <strong>de</strong> la velocidad <strong>de</strong> rotaciÄn <strong>de</strong>l motor asociado al compresor, y dicha variable es una<br />

especificaciÄn.<br />

Posteriorm<strong>en</strong>te las dos corri<strong>en</strong>tes se juntan mediante un mixer, para conformar la<br />

corri<strong>en</strong>te <strong>de</strong> salida.<br />

correspondi<strong>en</strong>te.<br />

Los datos y resultados se han contrastado con el Informe <strong>de</strong> Mant<strong>en</strong>imi<strong>en</strong>to Predictivo<br />

3.3.3 COMPRESOR TLA-3 DE PLANTA POSESIÅN<br />

Este compresor ha sido configurado para trabajar <strong>en</strong> dos etapas <strong>de</strong> compresiÄn, y la<br />

forma <strong>de</strong> simular esta situaciÄn <strong>en</strong> HYSYS correspon<strong>de</strong> a <strong>de</strong>finir cada etapa como un compresor<br />

in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te, <strong>en</strong> el cual la salida <strong>de</strong> la primera etapa se <strong>en</strong>frÑa y pasa a la <strong>en</strong>trada <strong>de</strong> la<br />

segunda etapa.<br />

Sin embargo hay que consi<strong>de</strong>rar un factor bastante importante, el hecho que cada<br />

compresor calcula el flujo <strong>de</strong> acuerdo a las condiciones <strong>de</strong> <strong>en</strong>trada y salida <strong>de</strong> manera<br />

in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>te para cada etapa, y pue<strong>de</strong> suce<strong>de</strong>r que <strong>los</strong> flujos requeridos por cada compresor<br />

sean difer<strong>en</strong>tes. Consi<strong>de</strong>rando la alternativa inicialm<strong>en</strong>te esbozada (una sola lÑnea <strong>de</strong><br />

corri<strong>en</strong>tes), <strong>de</strong> un simple balance <strong>de</strong> masa se obti<strong>en</strong>e que el flujo es constante a lo largo <strong>de</strong> las<br />

dos etapas <strong>de</strong>l compresor, y este esquema conllevarÑa a un error <strong>de</strong> especificaciones por parte


45<br />

<strong>de</strong> HYSYS. Para subsanar este efecto se <strong>de</strong>be dividir la corri<strong>en</strong>te que sale <strong>de</strong> la primera etapa <strong>en</strong><br />

dos flujos difer<strong>en</strong>tes; una corri<strong>en</strong>te que continçe hacia la segunda etapa y otra corri<strong>en</strong>te<br />

d<strong>en</strong>ominada ÅantorchaÉ, que se agrega para asumir el exceso o falta <strong>de</strong> flujo que requiere la<br />

segunda etapa, y asÑ cumplir con el balance <strong>de</strong> masas.<br />

Este cambio <strong>en</strong> el mo<strong>de</strong>lo se pue<strong>de</strong> observar <strong>en</strong> las sigui<strong>en</strong>tes ilustraciones:<br />

IlustraciÄn 16.- Mo<strong>de</strong>lo inicial para el compresor <strong>de</strong> dos etapas<br />

IlustraciÄn 17.- Mo<strong>de</strong>lo ajustado para cumplir el balance <strong>de</strong> masa<br />

Y <strong>en</strong> una posterior etapa se pue<strong>de</strong> buscar la presiÄn Äptima <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> la primera<br />

etapa, que minimiza el flujo <strong>de</strong> la corri<strong>en</strong>te ÅAntorchaÉ evitando párdidas o recirculaciones.<br />

respectivo.<br />

Los datos, nuevam<strong>en</strong>te, se han obt<strong>en</strong>ido <strong>de</strong>l Informe <strong>de</strong> Mant<strong>en</strong>imi<strong>en</strong>to Predictivo


46<br />

3.3.4 COMPRESOR TC-10 DE PLANTA CULLEN<br />

Este compresor es c<strong>en</strong>trÑfugo, y trabaja <strong>en</strong> dos etapas <strong>de</strong> compresiÄn, por lo cual el<br />

esquema a mo<strong>de</strong>lar es exactam<strong>en</strong>te el mismo al <strong>de</strong>sarrolllado para el compresor TLA-3. La<br />

difer<strong>en</strong>cia mÖs importante radica <strong>en</strong> el uso <strong>de</strong> curvas <strong>de</strong> diseÜo para especificar el compresor.<br />

Se <strong>de</strong>be aclarar que, durante el <strong>de</strong>sarrollo <strong>de</strong>l pres<strong>en</strong>te trabajo <strong>de</strong> titulaciÄn, el equipo<br />

<strong>en</strong> cuestiÄn estaba <strong>en</strong> proceso <strong>de</strong> instalaciÄn <strong>en</strong> Planta Cull<strong>en</strong>, razÄn por la cual no se ti<strong>en</strong><strong>en</strong><br />

datos reales <strong>de</strong> funcionami<strong>en</strong>to <strong>de</strong>l compresor TC-10.<br />

Como <strong>en</strong> el proceso <strong>de</strong> adquisiciÄn <strong>de</strong> dicho compresor ENAP Magallanes <strong>en</strong>viÄ al<br />

fabricante las especificaciones <strong>de</strong>l gas a comprimir (cromatografÑa promedio, variables <strong>de</strong><br />

proceso, etc.) y <strong>los</strong> requerimi<strong>en</strong>tos <strong>de</strong> salida, y el fabricante <strong>en</strong>tregÄ una hoja con resultados, se<br />

ha optado por contrastar dichos valores con <strong>los</strong> que garantiza el fabricante.<br />

Por estas razones, <strong>los</strong> datos para la simulaciÄn correspond<strong>en</strong> a las especificaciones<br />

<strong>en</strong>viadas por ENAP Magallanes y <strong>los</strong> resultados <strong>en</strong>tregados por el fabricante.<br />

3.3.5 TURBINA DEL TREN SULZER<br />

HYSYS no posee un mÄdulo especÑfico para <strong>de</strong>finir turbinas, pero se pued<strong>en</strong> estimar las<br />

propieda<strong>de</strong>s globales <strong>de</strong> funcionami<strong>en</strong>to <strong>de</strong> la turbina utilizando un expansor.<br />

En la pÖgina web <strong>de</strong>l programa se recomi<strong>en</strong>da un mo<strong>de</strong>lo para simular el proceso <strong>de</strong><br />

combustiÄn <strong>de</strong>l gas residual con que se alim<strong>en</strong>ta a la turbina, y se ha utilizado dicho mo<strong>de</strong>lo<br />

conceptual para simular la turbina <strong>de</strong>l tr<strong>en</strong> SULZER <strong>de</strong> Planta PosesiÄn.<br />

Se ha investigado a<strong>de</strong>mÖs las difer<strong>en</strong>cias para <strong>los</strong> sigui<strong>en</strong>tes esquemas <strong>de</strong> mo<strong>de</strong>laciÄn<br />

<strong>de</strong>l proceso <strong>de</strong> combustiÄn:


47<br />

Ä Defini<strong>en</strong>do las reacciones <strong>de</strong> combustiÄn <strong>de</strong>l gas natural [Ec. 15] y utilizando un<br />

reactor <strong>de</strong>l tipo Conversion.<br />

Ä Omiti<strong>en</strong>do las reacciones <strong>de</strong> combustiÄn y utilizando un reactor <strong>de</strong> tipo Gibbs.<br />

Se aÜa<strong>de</strong> a<strong>de</strong>mÖs el cÖlculo <strong>de</strong> ciertos parÖmetros que son <strong>de</strong> interás a la hora <strong>de</strong><br />

evaluar el correcto funcionami<strong>en</strong>to <strong>de</strong>l sistema:<br />

Ä Pot<strong>en</strong>cia disponible total: Consi<strong>de</strong>rando que la <strong>en</strong>ergÑa involucrada <strong>en</strong> la<br />

compresiÄn <strong>de</strong>l aire se obti<strong>en</strong>e <strong>de</strong> la misma turbina, la pot<strong>en</strong>cia disponible para<br />

comprimir el gas residual y para alim<strong>en</strong>tar al g<strong>en</strong>erador se pue<strong>de</strong> obt<strong>en</strong>er <strong>de</strong><br />

acuerdo a:<br />

HPTOTAL Å HPAire î HPDisponible<br />

[Ec. 25]<br />

Ä Uso <strong>de</strong> la pot<strong>en</strong>cia disponible total: Se pue<strong>de</strong> calcular quá porc<strong>en</strong>taje se ocupa<br />

tanto <strong>en</strong> el g<strong>en</strong>erador como <strong>en</strong> el compresor C-3, asumi<strong>en</strong>do que la pot<strong>en</strong>cia<br />

consumida <strong>en</strong> el g<strong>en</strong>erador (HP G<strong>en</strong>erador ) es <strong>de</strong> 3,200 HP.<br />

HP<br />

G<strong>en</strong>erador Å<br />

HP<br />

%<br />

G<strong>en</strong>erador<br />

HP<br />

% C É 3 Å<br />

Disponible<br />

HP<br />

Disponible<br />

É HP<br />

Disponible<br />

G<strong>en</strong>erador<br />

[Ec. 26]<br />

[Ec. 27]<br />

Ä Efici<strong>en</strong>cia <strong>de</strong> uso <strong>de</strong> gas combustible: De la <strong>en</strong>ergÑa total que ti<strong>en</strong>e el gas <strong>de</strong><br />

<strong>en</strong>trada a turbina, se calcula el porc<strong>en</strong>taje que <strong>en</strong>trega la turbina.<br />

% Ef<br />

GC<br />

HP<br />

Å<br />

TOTAL Turbina<br />

He [Ec. 28]


48<br />

Ä<br />

Efici<strong>en</strong>cia <strong>de</strong> la turbina: Se compara la <strong>en</strong>ergÑa <strong>de</strong> salida <strong>de</strong> <strong>los</strong> gases con el valor<br />

teÄrico calculado si <strong>los</strong> gases salieran a 0 èC. Si se d<strong>en</strong>ota He como la <strong>en</strong>talpÑa<br />

<strong>de</strong>l gas <strong>de</strong> <strong>en</strong>trada y Hs <strong>en</strong>talpÑa <strong>de</strong>l gas <strong>de</strong> salida, se obti<strong>en</strong>e:<br />

Real<br />

ë<br />

H Å Hs É He<br />

I<strong>de</strong>al<br />

ë<br />

T Ü Actual<br />

H Å Hs É He<br />

0ÜC<br />

í<br />

í<br />

Ef<br />

Turbina<br />

Å<br />

H<br />

H<br />

Real<br />

I<strong>de</strong>al<br />

[Ec. 29]<br />

3.3.6 ESQUEMA DE PRODUCCIÅN DE PLANTAS POSESIÅN Y CABO<br />

NEGRO<br />

En esta simulaciÄn el objetivo es realizar un mo<strong>de</strong>lo rÖpido e intuitivo <strong>de</strong>l proceso <strong>de</strong><br />

fraccionami<strong>en</strong>to <strong>de</strong>l gas que provi<strong>en</strong>e <strong>de</strong> <strong>los</strong> yacimi<strong>en</strong>tos y que ingresa a Planta PosesiÄn para<br />

g<strong>en</strong>erar Gas Residual y Raw Product. En Cabo Negro el Raw Product se separa <strong>en</strong> sus<br />

compon<strong>en</strong>tes comerciales propano, butano y gasolinas.<br />

Este esquema <strong>de</strong> separaciÄn, cuando no se requiere el grado <strong>de</strong> exactitud <strong>en</strong>tregado<br />

por la simulaciÄn especÑfica <strong>de</strong> Planta PosesiÄn, se <strong>de</strong>sarrolla utilizando la herrami<strong>en</strong>ta Splitter<br />

<strong>de</strong> HYSYS, que permite separar una corri<strong>en</strong>te <strong>de</strong> acuerdo a porc<strong>en</strong>tajes <strong>de</strong> separaciÄn.<br />

El ánfasis se hace <strong>en</strong> po<strong>de</strong>r obt<strong>en</strong>er una estimaciÄn <strong>de</strong> las ganancias econÄmicas<br />

g<strong>en</strong>eradas por la separaciÄn <strong>de</strong>l Raw Product <strong>de</strong>l Gas Rico, y una vez obt<strong>en</strong>ida la simulaciÄn se<br />

pue<strong>de</strong> compara esta situaciÄn con la ev<strong>en</strong>tualidad <strong>de</strong> no realizar <strong>procesos</strong> <strong>de</strong> separaciÄn y<br />

v<strong>en</strong><strong>de</strong>r todo el flujo <strong>de</strong> Gas Rico a precio <strong>de</strong> Gas Residual.<br />

Tambián se pue<strong>de</strong> investigar quá suce<strong>de</strong> al variar la especificaciÄn <strong>de</strong> calidad <strong>de</strong>l Raw<br />

Product (permiti<strong>en</strong>do mÖs o m<strong>en</strong>os etano) o <strong>los</strong> efectos que un alza <strong>en</strong> el precio <strong>de</strong> v<strong>en</strong>ta <strong>de</strong>l<br />

gas residual ti<strong>en</strong>e <strong>en</strong> las ganancias econÄmicas.


49<br />

3.3.7 ELEMPLO DE CARACTERIZACIÅN DE PETRÅLEOS EN HYSYS<br />

Para esta simulaciÄn se ha consi<strong>de</strong>rado un AnÖlisis <strong>de</strong> Laboratorio efectuado a la<br />

Gasolina Natural <strong>en</strong> el Terminal Gregorio, y se muestra el procedimi<strong>en</strong>to necesario para<br />

repres<strong>en</strong>tar la informaciÄn <strong>de</strong> dicho AnÖlisis <strong>en</strong> una simulaciÄn <strong>de</strong> HYSYS.<br />

Este caso no correspon<strong>de</strong> a ninguna instalaciÄn o equipo <strong>de</strong> terr<strong>en</strong>o, y se agrega sÄlo<br />

para <strong>de</strong>mostrar el procedimi<strong>en</strong>to correcto <strong>de</strong> especificaciÄn.<br />

3.3.8 EJEMPLO DE LOOP DE GASODUCTOS<br />

De igual forma que para la caracterizaciÄn <strong>de</strong> petrÄleos, esta simulaciÄn ejemplifica un<br />

loop <strong>de</strong> gasoductos g<strong>en</strong>árico; no ti<strong>en</strong>e sÑmil con situaciones reales <strong>de</strong> terr<strong>en</strong>o.<br />

Repres<strong>en</strong>ta la simulaciÄn <strong>de</strong> un caso <strong>en</strong> que se ti<strong>en</strong><strong>en</strong> dos gasoductos por <strong>los</strong> cuales se<br />

ha <strong>de</strong> transportar una corri<strong>en</strong>te <strong>de</strong>terminada. HYSYS <strong>en</strong>trega informaciÄn valiosa a la hora <strong>de</strong><br />

<strong>de</strong>cidir cuÖnto porc<strong>en</strong>taje <strong>de</strong>be ir <strong>en</strong> cada gasoducto, <strong>de</strong> manera que <strong>en</strong> el punto <strong>de</strong> uniÄn <strong>de</strong><br />

ambos gasoductos se obt<strong>en</strong>ga igual presiÄn <strong>de</strong> llegada.<br />

Se aÜa<strong>de</strong> a<strong>de</strong>mÖs <strong>en</strong> el punto <strong>de</strong> recepciÄn la corri<strong>en</strong>te <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> un compresor<br />

g<strong>en</strong>árico, y se muestra una aplicaciÄn <strong>de</strong> las herrami<strong>en</strong>tas Adjust y Set que es fundam<strong>en</strong>tal para<br />

resolver la presiÄn <strong>de</strong> salida <strong>de</strong>l gasoducto y <strong>de</strong> <strong>de</strong>scarga <strong>de</strong>l compresor.<br />

3.3.9 PLANTA CULLEN<br />

El proceso efectuado actualm<strong>en</strong>te <strong>en</strong> Planta Cull<strong>en</strong> consiste <strong>en</strong> <strong>de</strong>shidratar el gas,<br />

acondicionÖndolo para su transporte hacia contin<strong>en</strong>te, y separar las gasolinas mediante<br />

<strong>en</strong>friami<strong>en</strong>to.<br />

La simulaciÄn que mo<strong>de</strong>larÖ esta planta contempla las sigui<strong>en</strong>tes etapas <strong>de</strong> <strong>de</strong>sarrollo:


50<br />

Ä Circuito <strong>de</strong> refrigeraciÄn por propano<br />

Correspon<strong>de</strong> a un circuito cerrado <strong>de</strong> compresiÄn, cond<strong>en</strong>saciÄn, expansiÄn y<br />

cal<strong>en</strong>tami<strong>en</strong>to; la etapa <strong>de</strong> cal<strong>en</strong>tami<strong>en</strong>to se logra <strong>en</strong> un intercambiador <strong>de</strong> tubo y<br />

carcasa <strong>en</strong> que el propano se evapora y por <strong>los</strong> tubos circula el fluido <strong>de</strong> proceso que se<br />

<strong>de</strong>sea <strong>en</strong>friar. El compresor es c<strong>en</strong>trÑfugo, la cond<strong>en</strong>saciÄn se logra mediante una serie<br />

<strong>de</strong> aero<strong>en</strong>friadores y la expansiÄn mediante una vÖlvula <strong>de</strong> Joule-Thompson.<br />

Ä Tr<strong>en</strong> <strong>de</strong> <strong>en</strong>friami<strong>en</strong>to <strong>de</strong>l gas <strong>de</strong> <strong>en</strong>trada a planta<br />

El gas <strong>de</strong> proceso se pre<strong>en</strong>frÑa <strong>en</strong> un intercambiador <strong>de</strong> tubo y carcasa por interacciÄn<br />

con el gas prov<strong>en</strong>i<strong>en</strong>te <strong>de</strong> un separador <strong>de</strong> tres fases, para posteriorm<strong>en</strong>te, gracias a un<br />

circuito <strong>de</strong> propano, <strong>en</strong>friar açn mÖs la corri<strong>en</strong>te <strong>de</strong> proceso. Este esquema se simula<br />

mediante dos intercambiadores <strong>de</strong> calor <strong>de</strong> tubo y carcasa.<br />

Ä Flasheo <strong>de</strong> gasolinas<br />

Las gasolinas g<strong>en</strong>eradas por el <strong>en</strong>friami<strong>en</strong>to <strong>de</strong>l gas <strong>de</strong> proceso pasan por dos etapas<br />

<strong>de</strong> flasheo que logran disminuir la presiÄn <strong>de</strong> la corri<strong>en</strong>te hasta 15 kg/cm2<br />

aproximadam<strong>en</strong>te. Esto se realiza <strong>en</strong> HYSYS mediante una serie <strong>de</strong> separadores y<br />

vÖlvulas<br />

Ä Unidad reg<strong>en</strong>eradora <strong>de</strong> glicol<br />

A la corri<strong>en</strong>te <strong>de</strong> <strong>en</strong>trada a planta se le aÜa<strong>de</strong> una mezcla <strong>de</strong> etil<strong>en</strong>glicol y agua. En el<br />

proceso <strong>de</strong> <strong>en</strong>friami<strong>en</strong>to se logra separar cierta cantidad <strong>de</strong> agua, que es atrapada por<br />

el etil<strong>en</strong>glicol. El resultado (una mezcla <strong>de</strong> etil<strong>en</strong>glicol y agua con mayor cont<strong>en</strong>ido <strong>de</strong><br />

agua) se <strong>de</strong>be purificar para po<strong>de</strong>r recircularlo al punto <strong>de</strong> <strong>en</strong>trada a planta. Para esto<br />

se dispone <strong>de</strong> una serie <strong>de</strong> intercambiadores <strong>de</strong> calor que evaporan cierta cantidad <strong>de</strong><br />

agua hasta obt<strong>en</strong>er una pureza <strong>de</strong> glicol sufici<strong>en</strong>te para ser recirculada. Este esquema<br />

se ha simulado <strong>en</strong> HYSYS con la mayor similitud posible, pues para muchas variables


51<br />

que eran necesarias para la simulaciÄn (principalm<strong>en</strong>te temperaturas) no hay<br />

medidores instalados. Sin embargo, como se verÖ con posterioridad, <strong>los</strong> resultados<br />

obt<strong>en</strong>idos <strong>en</strong> HYSYS al simular la mezcla glicol agua no resultan satisfactorios <strong>de</strong>bido a<br />

errores internos <strong>de</strong>l programa.<br />

a continuaciÄn.<br />

El esquema <strong>de</strong> <strong>los</strong> equipos que conforman la unidad reg<strong>en</strong>eradora <strong>de</strong> glicol se pres<strong>en</strong>ta<br />

IlustraciÄn 18.- Esquema <strong>de</strong> la Unidad Reg<strong>en</strong>eradora <strong>de</strong> Glicol <strong>de</strong> Planta Cull<strong>en</strong><br />

Se pue<strong>de</strong> notar que, por lo m<strong>en</strong>os <strong>en</strong> aspecto grÖfico, pres<strong>en</strong>ta similitu<strong>de</strong>s bastante<br />

razonables con la simulaciÄn <strong>de</strong>sarrollada, que se pres<strong>en</strong>ta a continuaciÄn:


52<br />

IlustraciÄn 19.- Aspecto <strong>de</strong> la Unidad Reg<strong>en</strong>eradora <strong>de</strong> Glicol simulada <strong>en</strong> HYSYS<br />

Los datos que se han utilizado para realizar y validar la simulaciÄn provi<strong>en</strong><strong>en</strong> <strong>de</strong><br />

diversas fu<strong>en</strong>tes; Informes <strong>de</strong> Rutina <strong>de</strong> Laboratorio, datos obt<strong>en</strong>idos <strong>en</strong> las Hojas <strong>de</strong> Estado<br />

Diario medidas por el operador <strong>de</strong> <strong>procesos</strong>, datos <strong>en</strong>tregados por <strong>los</strong> fabricantes <strong>de</strong> <strong>los</strong> equipos<br />

simulados y principalm<strong>en</strong>te, por mediciÄn directa <strong>de</strong> qui<strong>en</strong> redacta este trabajo durante el<br />

tiempo <strong>de</strong> estadÑa <strong>en</strong> planta.<br />

3.3.10 PLANTA POSESIÅN<br />

El principal objetivo <strong>de</strong> Planta PosesiÄn es separar el gas <strong>de</strong> <strong>en</strong>trada a planta <strong>en</strong> dos<br />

corri<strong>en</strong>tes; Raw Product, constituido <strong>de</strong> propano, compon<strong>en</strong>tes mÖs pesados y trazas <strong>de</strong> etano,<br />

y una corri<strong>en</strong>te d<strong>en</strong>ominada Gas Residual, que es bÖsicam<strong>en</strong>te metano, etano y una pequeÜa<br />

cantidad <strong>de</strong> propano.


53<br />

Para efectuar este proceso <strong>de</strong> separaciÄn se llevan a cabo etapas sucesivas <strong>de</strong><br />

<strong>en</strong>friami<strong>en</strong>to al gas que <strong>en</strong>tra a planta, llegando a trabajar temperaturas <strong>de</strong> hasta -90 èC.<br />

En una primera etapa se efectça integraciÄn <strong>de</strong> calor con corri<strong>en</strong>tes frÑas <strong>de</strong> otros<br />

puntos <strong>de</strong> la planta, logrando disminuir la temperatura <strong>de</strong> la corri<strong>en</strong>te <strong>de</strong> <strong>en</strong>trada hasta<br />

temperaturas cercanas a -45 èC. Producto <strong>de</strong> esta caÑda <strong>de</strong> temperatura se obti<strong>en</strong>e<br />

cond<strong>en</strong>saciÄn, y se separa el lÑquido y el vapor <strong>en</strong> un separador.<br />

El <strong>en</strong>friami<strong>en</strong>to se ve favorecido por dos expansiones sucesivas <strong>de</strong>l gas <strong>de</strong> <strong>en</strong>trada, que<br />

aproximadam<strong>en</strong>te ingresa a planta a 70 kg/cm 2 . En el primer expansor se baja la presiÄn <strong>de</strong>l gas<br />

a 42 kg/cm 2 , y como resultado se obti<strong>en</strong>e cond<strong>en</strong>saciÄn <strong>de</strong> algunos compon<strong>en</strong>tes pesados, que<br />

pasan por un separador. En una segunda etapa se baja la presiÄn a 26 kg/cm 2 , nuevam<strong>en</strong>te hay<br />

cond<strong>en</strong>saciÄn y se separan <strong>los</strong> constituy<strong>en</strong>tes. Todas las corri<strong>en</strong>tes <strong>de</strong> lÑquido que se han<br />

g<strong>en</strong>erado hasta el mom<strong>en</strong>to se un<strong>en</strong> e ingresan a una columna <strong>de</strong> <strong>de</strong>stilaciÄn, d<strong>en</strong>ominada V-5.<br />

Esta columna <strong>de</strong> <strong>de</strong>stilaciÄn posee un cond<strong>en</strong>sador total. Los vapores g<strong>en</strong>erados <strong>en</strong><br />

este cond<strong>en</strong>sador constituy<strong>en</strong> el gas residual, que se precali<strong>en</strong>ta por intercambio con otras<br />

corri<strong>en</strong>tes <strong>de</strong> proceso hasta llegar a dos compresores c<strong>en</strong>trÑfugos acoplados a <strong>los</strong> dos<br />

expansores ya m<strong>en</strong>cionados. Por el fondo <strong>de</strong> la columna V-5 el rehervidor pres<strong>en</strong>ta una<br />

configuraciÄn <strong>de</strong> termosifÄn, recirculando completam<strong>en</strong>te el fluido que pasa por áste.<br />

IlustraciÄn 20.- Esquema <strong>de</strong> un rehervidor tipo termosifÄn


54<br />

Consi<strong>de</strong>re que el rehervidor real correspon<strong>de</strong> a un horno alim<strong>en</strong>tado por gas residual.<br />

Mediante una extracciÄn lateral <strong>de</strong> lÑquido ubicada <strong>en</strong> el plato nè 1 (la nom<strong>en</strong>clatura<br />

utilizada <strong>en</strong>umera <strong>los</strong> platos <strong>de</strong>s<strong>de</strong> abajo hacia arriba) se obti<strong>en</strong>e la corri<strong>en</strong>te <strong>de</strong> Raw Product,<br />

que es el principal producto <strong>de</strong> Planta PosesiÄn.<br />

La columna V-5 posee 30 platos, y trabaja a una presiÄn <strong>de</strong> diseÜo <strong>de</strong> 30 kg/cm 2 <strong>en</strong> un<br />

amplio rango <strong>de</strong> temperaturas, <strong>de</strong>s<strong>de</strong> -50 èC <strong>en</strong> el tope hasta 115 èC <strong>en</strong> el fondo. Las<br />

caracterÑsticas propias <strong>de</strong> este sistema, <strong>en</strong> particular el rango <strong>de</strong> temperaturas <strong>de</strong> la columna,<br />

hac<strong>en</strong> que simular este equipo <strong>en</strong> HYSYS sea un procedimi<strong>en</strong>to <strong>de</strong>licado y se <strong>de</strong>be efectuar una<br />

bu<strong>en</strong>a elecciÄn <strong>de</strong> las variables a especificar.<br />

Consi<strong>de</strong>rando la columna estÖndar que ofrece HYSYS, se han realizado las<br />

modificaciones necesarias para mo<strong>de</strong>lar <strong>de</strong> forma mÖs acertada el <strong>en</strong>torno <strong>de</strong> la columna V-5,<br />

utilizando la opciÄn <strong>de</strong> resoluciÄn <strong>de</strong> columnas Modified HYSIM Insi<strong>de</strong>-Out.<br />

La columna final posee 3 grados <strong>de</strong> libertad, y el proceso <strong>de</strong> converg<strong>en</strong>cia necesario<br />

para resolver la columna consi<strong>de</strong>ra por lo tanto la especificaciÄn <strong>de</strong> 3 variables.<br />

Se han escogido las sigui<strong>en</strong>tes especificaciones:<br />

Ä Reflujo Tope: Correspon<strong>de</strong> al flujo volumátrico (Std I<strong>de</strong>al Vol) <strong>de</strong> recirculaciÄn <strong>de</strong><br />

tope.<br />

Ä Etano <strong>en</strong> Raw Product: La fracciÄn molar <strong>de</strong> etano <strong>en</strong> el plato 1_Main TS.<br />

Ä ProducciÉn <strong>de</strong> Raw Product: El flujo volumátrico estÖndar <strong>de</strong> Raw Product que se<br />

obti<strong>en</strong>e por la extracciÄn <strong>en</strong> el plato 1_Main TS <strong>de</strong> la columna V-5.<br />

Ä Flujo P-3: Correspon<strong>de</strong> al flujo volumátrico <strong>de</strong> producto <strong>de</strong> fondo que ingresa al<br />

horno H-1 y recircula completam<strong>en</strong>te hacia la torre V-5.


55<br />

Se han <strong>de</strong>finido cuatro variables, <strong>de</strong>bido a que <strong>en</strong> el proceso <strong>de</strong> converg<strong>en</strong>cia <strong>de</strong> la<br />

columna se <strong>de</strong>b<strong>en</strong> realizar sustituciones <strong>de</strong> especificaciÄn escogi<strong>en</strong>do tres <strong>de</strong> ástas hasta<br />

obt<strong>en</strong>er valores normales <strong>de</strong> operaciÄn.<br />

Entre las limitaciones <strong>de</strong>l mo<strong>de</strong>lo realizado se <strong>de</strong>be com<strong>en</strong>tar la aus<strong>en</strong>cia <strong>de</strong> las curvas<br />

<strong>de</strong> diseÜo correspondi<strong>en</strong>tes a <strong>los</strong> expansores y compresores <strong>de</strong> Planta PosesiÄn. Debido a la<br />

antigêedad <strong>de</strong> estos equipos no se ti<strong>en</strong>e registro <strong>de</strong> esta informaciÄn.<br />

Se pue<strong>de</strong> com<strong>en</strong>tar que, <strong>de</strong> todas las simulaciones <strong>de</strong>sarrolladas, la <strong>de</strong> Planta PosesiÄn<br />

es la que consi<strong>de</strong>ra la mayor cantidad <strong>de</strong> equipos e interacciÄn <strong>en</strong>tre corri<strong>en</strong>tes. Es necesario<br />

agregar 3 mÄdu<strong>los</strong> <strong>de</strong> reciclo para lograr una a<strong>de</strong>cuada simulaciÄn.<br />

Los datos se han obt<strong>en</strong>ido <strong>de</strong> informaciÄn tácnica obt<strong>en</strong>ida <strong>en</strong> planta, las variables <strong>de</strong><br />

proceso y valores <strong>de</strong> cromatografÑa han sido recopilados <strong>de</strong> las Hojas <strong>de</strong> Estado Diario obt<strong>en</strong>idas<br />

por el operador <strong>de</strong> <strong>procesos</strong>, <strong>de</strong> Informes <strong>de</strong> Rutina <strong>de</strong> Laboratorio y por mediciÄn directa <strong>de</strong> las<br />

variables <strong>de</strong> interás por parte <strong>de</strong> qui<strong>en</strong> redacta el trabajo. Los resultados han sido validados <strong>en</strong><br />

base a las mismas fu<strong>en</strong>tes.


56<br />

3.4 DATOS Y RESULTADOS<br />

En la pres<strong>en</strong>te secciÄn se <strong>de</strong>tallan <strong>los</strong> datos y resultados que se han especificado para<br />

las simulaciones realizadas para ENAP Magallanes, y se pres<strong>en</strong>ta el valor real medido <strong>en</strong> planta<br />

u obt<strong>en</strong>ido <strong>de</strong> las diversas fu<strong>en</strong>tes <strong>de</strong> informaciÄn ya nombradas. A<strong>de</strong>mÖs se calcula la<br />

difer<strong>en</strong>cia porc<strong>en</strong>tual <strong>en</strong>tre el valor obt<strong>en</strong>ido por HYSYS y el valor real.<br />

En el anexo se pres<strong>en</strong>tan las capturas <strong>de</strong> pantalla para cada simulaciÄn realizada, <strong>en</strong><br />

don<strong>de</strong> se pue<strong>de</strong> comprobar el esquema utilizado. El ÅManual para simular <strong>procesos</strong> <strong>productivos</strong><br />

<strong>de</strong> ENAP Magallanes <strong>en</strong> HYSYSÉ conti<strong>en</strong>e un disco compacto con <strong>los</strong> archivos <strong>de</strong> las simulaciones<br />

respectivas.<br />

Una forma rÖpida <strong>de</strong> obt<strong>en</strong>er la difer<strong>en</strong>cia <strong>en</strong>tre el mo<strong>de</strong>lo y el sistema real es <strong>de</strong>finir el<br />

error relativo (E R ) <strong>en</strong>tre ambas medidas; <strong>de</strong> la sigui<strong>en</strong>te forma:<br />

E<br />

R<br />

ä<br />

Abs<br />

é<br />

ë<br />

Valor<br />

ÉValor<br />

HYSYS<br />

Real<br />

Å å ç<br />

ValorReal<br />

í<br />

ã<br />

Ä100% [Ec. 30]<br />

è<br />

Para analizar la similitud <strong>en</strong> <strong>los</strong> valores obt<strong>en</strong>idos <strong>de</strong> composiciÄn para <strong>de</strong>terminadas<br />

corri<strong>en</strong>tes, como se trabaja con magnitu<strong>de</strong>s <strong>en</strong> un rango muy amplio (por ejemplo se contrastan<br />

valores <strong>de</strong> metano a 92% y valores <strong>de</strong> CO2 a 0.2%), se ha calculado el error individual para cada<br />

mediciÄn, pero a<strong>de</strong>mÖs se obti<strong>en</strong>e un error global, d<strong>en</strong>ominado E G , que consi<strong>de</strong>ra la relevancia<br />

<strong>de</strong> cada compuesto <strong>en</strong> la cromatografÑa total, <strong>de</strong> la sigui<strong>en</strong>te forma:<br />

E Å ê x ÄE<br />

[Ec. 31]<br />

G i,Re al R,<br />

i<br />

Don<strong>de</strong><br />

X<br />

i,Real<br />

correspon<strong>de</strong> a la fracciÄn reportada por laboratorio para cada<br />

compon<strong>en</strong>te i, y el error relativo para la medida <strong>de</strong> dicho compon<strong>en</strong>te es E<br />

R,<br />

i<br />

.


57<br />

Con esto se logra quitar relevancia al 50% <strong>de</strong> error relativo que se obt<strong>en</strong>drÑa si, por<br />

ejemplificar, el valor real <strong>de</strong> la fracciÄn <strong>de</strong> nitrÄg<strong>en</strong>o es 0.02 y HYSYS <strong>en</strong>trega 0.03.<br />

A continuaciÄn se muestra <strong>en</strong> tablas <strong>los</strong> valores <strong>de</strong> composiciÄn que se ha ingresado <strong>en</strong><br />

cada simulaciÄn pres<strong>en</strong>tada <strong>en</strong> este capÑtulo.<br />

Compon<strong>en</strong>te HRA-5 Cull<strong>en</strong> HRA-2 Calafate TLA-3 PosesiÄn TC-10 Cull<strong>en</strong><br />

(%mol)<br />

Metano 92.12 92.12 91.50 92.15<br />

Etano 3.83 3.83 3.92 3.83<br />

Propano 1.20 1.20 1.53 1.22<br />

i-Butano 0.25 0.25 0.41 0.25<br />

n-Butano 0.30 0.30 0.47 0.35<br />

i-P<strong>en</strong>tano 0.09 0.09 0.15 0.14<br />

n-P<strong>en</strong>tano 0.06 0.06 0.12 0.10<br />

n-Hexano 0.01 0.01 0.12 0.15<br />

n-Heptano 0.00 0.00 0.02 0.13<br />

NitrÄg<strong>en</strong>o 2.04 2.04 1.69 1.48<br />

DiÄxido <strong>de</strong> carbono 0.10 0.10 0.06 0.20<br />

Agua 0.00 0.00 0.01 0.00<br />

Tabla 2.- Composiciones ingresadas para las simulaciones <strong>de</strong> compresores<br />

Compon<strong>en</strong>te Turbina SULZER Planta Cull<strong>en</strong> Planta PosesiÄn<br />

(%mol)<br />

Metano 92.35 91.16 91.04<br />

Etano 5.18 4.47 5.37<br />

Propano 1.05 1.48 1.12<br />

i-Butano 0.00 0.30 0.26<br />

n-Butano 0.00 0.42 0.34<br />

i-P<strong>en</strong>tano 0.00 0.14 0.13<br />

n-P<strong>en</strong>tano 0.00 0.12 0.09<br />

n-Hexano 0.00 0.17 0.21<br />

NitrÄg<strong>en</strong>o 1.04 1.58 1.24<br />

DiÄxido <strong>de</strong> carbono 0.38 0.16 0.20<br />

Agua 0.00 0.00 0.00<br />

Tabla 3.- Composiciones ingresadas a las simulaciones <strong>de</strong> Tr<strong>en</strong> SULZER, Cull<strong>en</strong> y PosesiÄn<br />

Para el caso <strong>de</strong> Planta PosesiÄn y Planta Cull<strong>en</strong> se <strong>de</strong>sg<strong>los</strong>an <strong>los</strong> resultados obt<strong>en</strong>idos<br />

<strong>de</strong> la simulaciÄn <strong>de</strong>sarrollada ingresando <strong>los</strong> valores normales promedio <strong>de</strong> operaciÄn, para<br />

posteriorm<strong>en</strong>te validar el mo<strong>de</strong>lo al ingresar valores <strong>de</strong> informes <strong>de</strong> laboratorio<br />

correspondi<strong>en</strong>tes a difer<strong>en</strong>tes dÑas <strong>de</strong> producciÄn.


58<br />

3.4.1 COMPRESOR HRA-5 DE PLANTA CULLEN<br />

Datos ingresados a la simulaciÄn:<br />

Corri<strong>en</strong>te Variable Valor especificado<br />

Temperatura [èF] 75<br />

Pozos <strong>de</strong> baja<br />

PresiÄn [psig] 700<br />

(Entrada al separador)<br />

ComposiciÄn Ver Tabla 2<br />

Descarga<br />

PresiÄn [psig] 1,123<br />

(Descarga <strong>de</strong>l compresor)<br />

Tabla 4.- Definiciones <strong>de</strong> corri<strong>en</strong>tes <strong>en</strong> la simulaciÄn <strong>de</strong>l compresor HRA-5<br />

Equipo Variable Valor especificado<br />

Number of cylin<strong>de</strong>rs 4<br />

Cyl. Type<br />

Double Acting, No-Tail Rod Type<br />

Bore [in] 8.5<br />

Compresor HRA-5<br />

Stroke [in] 14<br />

Piston Rod Diameter [in] 2.5<br />

Speed [RPM] 270<br />

Adiabatic Effici<strong>en</strong>cy [%] 85<br />

Vol. Clearances [in3]<br />

Ver Manual<br />

Tabla 5.- Datos <strong>de</strong>l compresor HRA-5<br />

Equipo o corri<strong>en</strong>te Variable Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

Compresor HRA-5 Pot<strong>en</strong>cia [HP] 1,053 1,045 0.8<br />

Descarga Temperatura [èF] 145.8 143 2.0<br />

(Descarga <strong>de</strong>l Flujo comprimido 1,217 1,269 4.1<br />

compresor) [MMMCSD]<br />

Tabla 6.- Resultados obt<strong>en</strong>idos HRA-5


59<br />

3.4.2 COMPRESOR HRA-5 DE ESTACIÅN COMPRESORA CALAFATE<br />

Datos ingresados y resultados obt<strong>en</strong>idos <strong>en</strong> la simulaciÄn:<br />

Corri<strong>en</strong>te Variable Valor especificado<br />

Temperatura [èF] 82<br />

Calafate<br />

PresiÄn [psig] 988<br />

(Entrada al separador)<br />

ComposiciÄn Ver Tabla 2<br />

Descarga<br />

PresiÄn [psig] 1,259<br />

(Descarga <strong>de</strong>l compresor)<br />

Tabla 7.- Definiciones <strong>de</strong> corri<strong>en</strong>tes <strong>en</strong> la simulaciÄn <strong>de</strong>l compresor HRA-2<br />

Equipo Variable Valor especificado<br />

Number of cylin<strong>de</strong>rs 2<br />

Cyl. Type<br />

Double Acting, No-Tail Rod Type<br />

Bore [in] 5.5<br />

Cil. 1y2<br />

Stroke [in] 14<br />

Piston Rod Diameter [in] 3<br />

Speed [RPM] 329<br />

Adiabatic Effici<strong>en</strong>cy [%] 85<br />

Vol. Clearances [in3]<br />

Ver Manual<br />

Number of cylin<strong>de</strong>rs 2<br />

Cyl. Type<br />

Double Acting, No-Tail Rod Type<br />

Bore [in] 7.25<br />

Cil. 3y4<br />

Stroke [in] 14<br />

Piston Rod Diameter [in] 3<br />

Adiabatic Effici<strong>en</strong>cy [%] 85<br />

Vol. Clearances [in3]<br />

Ver Manual<br />

Tabla 8.- Datos <strong>de</strong>l compresor HRA-2 dividido <strong>en</strong> dos compresores teÄricos<br />

Equipo o corri<strong>en</strong>te Variable Valor HYSYS Valor real Difer<strong>en</strong>cia<br />

[%]<br />

Compresor HRA-2 Pot<strong>en</strong>cia [HP] 881 891 1.1<br />

Temperatura [èF] 123.4 133.4 7.5<br />

Descarga1<br />

Flujo comprimido 0.618 0.594 4.0<br />

(<strong>de</strong>l compresor Cyl. 1y2)<br />

[MMMCSD]<br />

Descarga1 Temperatura [èF] 123.1 124 0.7<br />

(Descarga <strong>de</strong>l compresor Flujo comprimido 1.16 1.12 3.6<br />

Cyl. 1y2)<br />

[MMMCSD]<br />

Descarga total Temperatura [èF] 123.1 128 3.8<br />

(Suma <strong>de</strong>scargas <strong>de</strong> ambos Flujo comprimido 1.78 1.68 6.0<br />

compresores) [MMMCSD]<br />

Tabla 9.- Resultados obt<strong>en</strong>idos para el compresor HRA-2


60<br />

3.4.3 COMPRESOR TLA-3 DE PLANTA POSESIÅN<br />

Datos ingresados a la simulaciÄn:<br />

Corri<strong>en</strong>te Variable Valor especificado<br />

Temperatura [èF] 66<br />

Entrada<br />

PresiÄn [psig] 425<br />

(Entrada al separador)<br />

ComposiciÄn Ver Tabla 2<br />

Entrada2<br />

Temperatura [èF] 91<br />

(Entrada a la segunda etapa)<br />

Descarga<br />

PresiÄn [psig] 1,700<br />

(Descarga <strong>de</strong>l compresor)<br />

Tabla 10.- Definiciones <strong>de</strong> corri<strong>en</strong>tes <strong>en</strong> la simulaciÄn <strong>de</strong>l compresor TLA-3<br />

Equipo Variable Valor especificado<br />

Number of cylin<strong>de</strong>rs 2<br />

Cyl. Type<br />

Double Acting, No-Tail Rod Type<br />

Bore [in] 9.75<br />

Compresor Etapa1<br />

Stroke [in] 19<br />

Piston Rod Diameter [in] 4<br />

Speed [RPM] 263<br />

Adiabatic Effici<strong>en</strong>cy [%] 95<br />

Vol. Clearances [in3]<br />

Ver Manual<br />

Number of cylin<strong>de</strong>rs 2<br />

Cyl. Type<br />

Double Acting, No-Tail Rod Type<br />

Bore [in] 5.625<br />

Compresor Etapa2 Stroke [in] 19<br />

Piston Rod Diameter [in] 4<br />

Adiabatic Effici<strong>en</strong>cy [%] 90<br />

Vol. Clearances [in3]<br />

Ver Manual<br />

Tabla 11.- Datos <strong>de</strong>l compresor TLA-3 dividido <strong>en</strong> dos etapas.<br />

Equipo o corri<strong>en</strong>te Variable Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

Compresor Etapa1 Pot<strong>en</strong>cia [HP] 1,242 1,238 0.3<br />

Compresor Etapa2 Pot<strong>en</strong>cia [HP] 481 496.8 3.2<br />

Temperatura [èF] 196.4 200.5 2.0<br />

Descarga Etapa1 Flujo comprimido 0.780 0.720 8.3<br />

[MMMCSD]<br />

Descarga Etapa2<br />

Temperatura [èF] 155.5 148 5.1<br />

Flujo comprimido 0.657 0.650 1.1<br />

[MMMCSD]<br />

Tabla 12.- Resultados obt<strong>en</strong>idos TLA-3 para presiÄn interetapa <strong>de</strong> 1,100 psig.


61<br />

3.4.4 COMPRESOR TC-10 DE PLANTA CULLEN<br />

A<strong>de</strong>mÖs <strong>de</strong> ingresar las curvas <strong>de</strong> diseÜo para ambas etapas <strong>de</strong>l compresor c<strong>en</strong>trÑfugo,<br />

que se pres<strong>en</strong>tan <strong>en</strong> el manual anexo, se han realizado las sigui<strong>en</strong>tes especificaciones <strong>en</strong> la<br />

simulaciÄn:<br />

Corri<strong>en</strong>te o equipo Variable Valor especificado<br />

Temperatura [èC] 20<br />

Entrada<br />

PresiÄn [kg/cm 2 ] 6<br />

(Entrada al separador)<br />

ComposiciÄn Ver Tabla 2<br />

Descarga Etapa1 PresiÄn [kg/cm 2 ] 23.28<br />

SucciÄn2<br />

Temperatura [èC] 30<br />

(Entrada Etapa2)<br />

Solar Etapa1 Speed [RPM] 21,216<br />

Solar Etapa2 Speed [RPM] 21,216<br />

Tabla 13.- Especificaciones realizadas al compresor TC-10<br />

Equipo o corri<strong>en</strong>te Variable Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

Pot<strong>en</strong>cia [HP] 2,548 2,585 1.4<br />

Etapa1 Efici<strong>en</strong>cia [%] 69 69.3 0.4<br />

Head [ftëlbm/lbf] 71,000 70,945 0.1<br />

Pot<strong>en</strong>cia [HP] 1,833 1,833 0.0<br />

Etapa2 Efici<strong>en</strong>cia [%] 70 70.6 0.8<br />

Head [ftëlbm/lbf] 51,797 51,265 1.0<br />

Temperatura [èC] 155.1 154.6 0.3<br />

Flujo a condiciones 3,020 3,017 0.1<br />

Descarga Etapa1<br />

actuales<br />

[Actual ft 3 /min]<br />

Flujo comprimido 677 675 0.3<br />

[Nm 3 /d]<br />

Temperatura [èC] 132.5 130.0 1.9<br />

Descarga Etapa2 Flujo a condiciones 801.9 813.6 1.4<br />

actuales<br />

[Actual ft 3 /min]<br />

Tabla 14.- Resultados obt<strong>en</strong>idos al simular el compresor TC-10<br />

Se <strong>de</strong>be notar que se ha agregado a las curvas <strong>de</strong> diseÜo el punto <strong>de</strong> operaciÄn real, a<br />

21,216 RPM, pues como el rango <strong>de</strong> efici<strong>en</strong>cias que pres<strong>en</strong>tan las curvas es estrecho, se<br />

observaban errores <strong>en</strong> la simulaciÄn, y <strong>los</strong> resultados no eran satisfactorios.


62<br />

3.4.5 TURBINA DEL TREN SULZER<br />

Los datos especificados, y resultados obt<strong>en</strong>idos <strong>de</strong> la simulaciÄn <strong>de</strong>l proceso <strong>de</strong><br />

combustiÄn y turbina <strong>de</strong>l tr<strong>en</strong> SULZER son <strong>los</strong> sigui<strong>en</strong>tes:<br />

Corri<strong>en</strong>te o equipo Variable Valor especificado<br />

Temperatura [èC] 12<br />

Gas Combustible PresiÄn [psia] 100<br />

ComposiciÄn Ver Tabla 3<br />

Aire a compresiÄn<br />

Temperatura [èC] 12<br />

PresiÄn [psia] 14<br />

Gases <strong>de</strong> <strong>en</strong>trada a la Temperatura [èC] 750<br />

turbina<br />

Compresor <strong>de</strong> aire Efici<strong>en</strong>cia [%] 90<br />

Turbina Efici<strong>en</strong>cia 75<br />

Tabla 15.- Datos ingresados a la simulaciÄn <strong>de</strong> la turbina <strong>de</strong>l tr<strong>en</strong> SULZER<br />

Equipo o corri<strong>en</strong>te Variable Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

Compresor <strong>de</strong> aire Pot<strong>en</strong>cia [HP] 35,093 35,000 0.3<br />

Pot<strong>en</strong>cia [HP] 49,870 50,000 0.3<br />

Temperatura <strong>de</strong> 459.3 450 2.1<br />

salida <strong>de</strong> gases [èC]<br />

Efici<strong>en</strong>cia <strong>de</strong> uso <strong>de</strong> 51.93 No disponible -<br />

gas combustible [%]<br />

Turbina Pot<strong>en</strong>cia disponible 14,778 15,000 1.5<br />

total [HP]<br />

Porc<strong>en</strong>taje utilizado 21.65 No disponible -<br />

por g<strong>en</strong>erador [%]<br />

Efici<strong>en</strong>cia <strong>de</strong><br />

38 No disponible -<br />

combustiÄn [%]<br />

Tabla 16.- Resultados obt<strong>en</strong>idos para el tr<strong>en</strong> SULZER<br />

ComposiciÄn obt<strong>en</strong>ida Compon<strong>en</strong>te Reacciones <strong>de</strong> Reactor <strong>de</strong> Difer<strong>en</strong>cia [%]<br />

(%mol)<br />

combustiÄn Gibbs<br />

OxÑg<strong>en</strong>o 16.41 16.41 0<br />

Gases a turbina NitrÄg<strong>en</strong>o 77.38 77.38 0<br />

Agua 04.09 04.09 0<br />

DiÄxido <strong>de</strong> carbono 02.12 02.12 0<br />

Tabla 17.- Composiciones <strong>de</strong>l gas <strong>de</strong> <strong>en</strong>trada a turbina mediante <strong>los</strong> dos mátodos


63<br />

3.4.6 PLANTA CULLEN<br />

Datos ingresados y resultados obt<strong>en</strong>idos:<br />

Corri<strong>en</strong>te o equipo Variable Valor especificado<br />

Temperatura [èC] 30<br />

PresiÄn [kg/cm 2 ] 55<br />

Gas <strong>de</strong> <strong>en</strong>trada Flujo estÖndar <strong>de</strong> gas [MMMCSD] 2.0<br />

ComposiciÄn Ver Tabla 3<br />

Water Dew Point [èC] -14<br />

Glicol <strong>de</strong> recirculaciÄn<br />

Flujo [USGPM] 1.575<br />

ComposiciÄn EGlicol [%masa] 73.1<br />

Propano <strong>de</strong> <strong>en</strong>friami<strong>en</strong>to Flujo volumátrico [MCSD] 149,300<br />

Tabla 18.- Datos ingresados para la simulaciÄn <strong>de</strong> Planta Cull<strong>en</strong><br />

ComposiciÄn <strong>de</strong> gas Compon<strong>en</strong>te Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

residual obt<strong>en</strong>ida (%mol)<br />

Metano 91.42 91.55 0.1<br />

Etano 4.45 4.44 0.2<br />

Propano 1.44 1.43 0.7<br />

i-Butano 0.28 0.28 0.0<br />

n-Butano 0.38 0.36 5.6<br />

Gas a compresores i-P<strong>en</strong>tano 0.11 0.10 10<br />

n-P<strong>en</strong>tano 0.09 0.07 28<br />

n-Hexano 0.08 0.02 300<br />

NitrÄg<strong>en</strong>o 1.59 1.61 1.2<br />

DiÄxido <strong>de</strong> carbono 0.16 0.16 0.0<br />

ERROR GLOBAL RESULTANTE [%] 0.28<br />

Tabla 19.- CromatografÑa <strong>de</strong>l gas residual obt<strong>en</strong>ida <strong>en</strong> la simulaciÄn<br />

Corri<strong>en</strong>te o equipo Variable Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

Glicol <strong>de</strong> separador<br />

ComposiciÄn EGlicol 72.7 70.7 2.8<br />

[%masa]<br />

Unidad <strong>de</strong> Temperatura <strong>de</strong> 241.1 242 0.4<br />

reg<strong>en</strong>eraciÄn <strong>de</strong> glicol reg<strong>en</strong>eraciÄn [èF]<br />

Gas a compresores Water Dew Point [èC] -25.3 -24 5.4<br />

Gasolinas<br />

Flujo volumátrico 36 35 ã 45 Valor<br />

estÖndar [m 3 /d]<br />

Aceptable<br />

Separador <strong>de</strong> Temperatura [èC] -19.3 -17 13<br />

gasolinas<br />

Tabla 20.- Resultados <strong>de</strong> la simulaciÄn <strong>de</strong> Planta Cull<strong>en</strong>


64<br />

3.4.7 PLANTA POSESIÅN<br />

Datos ingresados y resultados obt<strong>en</strong>idos:<br />

Corri<strong>en</strong>te o equipo Variable Valor especificado<br />

Temperatura [èC] 23<br />

Gas <strong>de</strong> <strong>en</strong>trada<br />

PresiÄn [kg/cm 2 ] 72<br />

Flujo estÖndar <strong>de</strong> gas [MMMCSD] 7.5<br />

ComposiciÄn Ver Tabla 3<br />

Expansor 1<br />

PresiÄn <strong>de</strong> <strong>de</strong>scarga [kg/cm 2 ] 42<br />

Flujo que pasa por vÖlvula JT [%] 10<br />

Expansor 2<br />

PresiÄn <strong>de</strong> <strong>de</strong>scarga [kg/cm 2 ] 27<br />

Flujo que pasa por vÖlvula JT [%] 1<br />

Separador V-21 Temperatura [èC] -74<br />

Intercambiador E-6<br />

Temperatura <strong>de</strong> corri<strong>en</strong>te <strong>de</strong> salida -66<br />

por <strong>los</strong> tubos [èC]<br />

Reflujo <strong>de</strong> tope [m 3 /d] 1,000<br />

Columna V-5 Etano <strong>en</strong> Raw Product [fracciÄn molar] 0.016<br />

Flujo P-3 (A horno) [m 3 /d] 1,280<br />

Tabla 21.- Listado <strong>de</strong> especificaciones utilizadas para simular Planta PosesiÄn<br />

Corri<strong>en</strong>te o equipo Variable Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

Separador V-1 PresiÄn [kg/cm 2 ] 65.88 69.5 5.2<br />

Separador V-21 PresiÄn [kg/cm 2 ] 42 41.1 2.2<br />

PresiÄn Entrada Ex-1 65.88 63.1 4.4<br />

[kg/cm 2 ]<br />

Expansor Ex-1<br />

PresiÄn Entrada Gas 30.6 27.0 13.3<br />

<strong>en</strong> C-2 [kg/cm 2 ]<br />

PresiÄn Salida Gas <strong>en</strong> 37.6 30.7 22.5<br />

C-2 [kg/cm 2 ]<br />

PresiÄn Entrada Ex-2 42 40 5.0<br />

[kg/cm 2 ]<br />

Expansor Ex-2<br />

PresiÄn Entrada Gas 26.5 23 15.2<br />

<strong>en</strong> C-1 [kg/cm 2 ]<br />

PresiÄn Salida Gas <strong>en</strong> 30 26 15.4<br />

C-1 [kg/cm 2 ]<br />

Tabla 22.-Temperaturas y presiones obt<strong>en</strong>idas <strong>de</strong> la simulaciÄn <strong>de</strong> Planta PosesiÄn<br />

Con respecto a <strong>los</strong> resultados relacionados con la columna <strong>de</strong> <strong>de</strong>stilaciÄn V-5 y <strong>los</strong><br />

equipos asociados, <strong>los</strong> resultados obt<strong>en</strong>idos son:


65<br />

Columna V-5 Variable Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

Corri<strong>en</strong>te <strong>de</strong> <strong>en</strong>trada<br />

Temperatura [èC] -45 -47 4.3<br />

PresiÄn [kg/cm 2 ] 32 30 6.7<br />

Fondo <strong>de</strong> la torre Temperatura [èC] 106.7 110 3.0<br />

Reflujo <strong>de</strong> fondo Temperatura [èC] 111 120 7.5<br />

Vapores <strong>de</strong> tope Temperatura [èC] -40 -38 5.3<br />

Reflujo <strong>de</strong> tope Temperatura [èC] -50 -47 6.4<br />

Gas residual Flujo [MMMCSD] 1.28 1.15 11.3<br />

Raw Product Flujo [m 3 /d] 668.5 750 10<br />

Tabla 23.- Resultados obt<strong>en</strong>idos al simular la columna V-5<br />

La simulaciÄn <strong>en</strong>trega <strong>los</strong> sigui<strong>en</strong>tes valores para la cromatografÑa <strong>de</strong> <strong>los</strong> gases<br />

reportados <strong>en</strong> el Informe <strong>de</strong> Laboratorio:<br />

ComposiciÄn<br />

Compon<strong>en</strong>te Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

(%mol)<br />

Metano 78.07 81.69 4.4<br />

Etano 21.56 17.35 24.3<br />

Gases <strong>de</strong> V-4<br />

Propano 0.01 0.01 0.0<br />

NitrÄg<strong>en</strong>o 0.32 0.27 18.5<br />

DiÄxido <strong>de</strong> carbono 0.05 0.69 92.8<br />

ERROR GLOBAL RESULTANTE [%] 8.5<br />

Metano 96.29 96.29 0.0<br />

Etano 2.05 2.30 10.9<br />

Gases <strong>de</strong> V-3<br />

Propano 0.04 0.08 50.0<br />

NitrÄg<strong>en</strong>o 1.47 1.04 41.3<br />

DiÄxido <strong>de</strong> carbono 0.15 0.29 48.3<br />

ERROR GLOBAL RESULTANTE [%] 0.9<br />

Etano 1.60 1.35 18.5<br />

Propano 50.52 50.10 0.8<br />

i-Butano 12.07 12.62 4.4<br />

Raw Product<br />

n-Butano 15.81 16.94 6.7<br />

i-P<strong>en</strong>tano 6.05 5.74 5.4<br />

n-P<strong>en</strong>tano 4.19 4.15 1.0<br />

n-Hexano 9.77 9.11 7.2<br />

ERROR GLOBAL RESULTANTE [%] 3.3<br />

Metano 93.00 93.28 0.3<br />

Etano 5.49 5.34 2.8<br />

Gas Residual<br />

Propano 0.03 0.06 50.0<br />

NitrÄg<strong>en</strong>o 1.26 0.92 37.0<br />

DiÄxido <strong>de</strong> carbono 0.20 0.40 50.0<br />

ERROR GLOBAL RESULTANTE [%] 1.0<br />

Tabla 24.- ComposiciÄn <strong>de</strong> las corri<strong>en</strong>tes obt<strong>en</strong>idas <strong>en</strong> la simulaciÄn


66<br />

3.5 VALIDACIÅN DE LAS SIMULACIONES OBTENIDAS<br />

Consi<strong>de</strong>rando que las simulaciones mÖs importantes correspond<strong>en</strong> a las <strong>de</strong> Planta<br />

PosesiÄn y Planta Cull<strong>en</strong>, se ha investigado el comportami<strong>en</strong>to <strong>de</strong> <strong>los</strong> mo<strong>de</strong><strong>los</strong> propuestos <strong>en</strong><br />

HYSYS para estas plantas.<br />

Se ha contrastado el resultado <strong>en</strong>tregado por las simulaciones creadas <strong>en</strong> HYSYS con<br />

valores reales obt<strong>en</strong>idos durante cinco dÑas <strong>de</strong> operaciÄn <strong>de</strong> planta (para el caso <strong>de</strong> Planta<br />

PosesiÄn se han consi<strong>de</strong>rado cuatro dÑas, <strong>de</strong>bido a la estabilidad <strong>de</strong> <strong>los</strong> valores observados),<br />

ingresando la informaciÄn correspondi<strong>en</strong>te a <strong>los</strong> Informes <strong>de</strong> Rutina <strong>de</strong> Laboratorio y Hojas <strong>de</strong><br />

Estado Diario registradas por el operador <strong>de</strong> <strong>procesos</strong> <strong>de</strong> planta. Las fechas correspond<strong>en</strong> a:<br />

Ä Planta Cull<strong>en</strong>: Entre el 24 y el 28 <strong>de</strong> Julio <strong>de</strong> 2007.<br />

Ä Planta PosesiÄn: Entre el 31 <strong>de</strong> Julio y el 5 <strong>de</strong> Agosto <strong>de</strong> 2007.<br />

Las condiciones <strong>de</strong> presiÄn, temperatura y flujo <strong>de</strong> <strong>en</strong>trada a planta se han consi<strong>de</strong>rado<br />

estables para el lapso <strong>de</strong> tiempo estudiado, y correspond<strong>en</strong> a:<br />

Condiciones <strong>de</strong> Variable Planta Cull<strong>en</strong> Planta PosesiÄn<br />

proceso<br />

Temperatura [èC] 30 23<br />

Entrada a Planta<br />

PresiÄn [kg/cm 2 ] 55 70<br />

Flujo volumátrico 2.0 7.0<br />

estÖndar [MMMCSD]<br />

Tabla 25.- Valores normales <strong>de</strong> operaciÄn Planta Cull<strong>en</strong> y Planta PosesiÄn


67<br />

3.5.1 VALIDACIÅN DE PLANTA CULLEN<br />

Las variables que se <strong>de</strong>b<strong>en</strong> especificar, para obt<strong>en</strong>er <strong>los</strong> resultados y analizar la<br />

simulaciÄn pres<strong>en</strong>tada <strong>de</strong> Planta Cull<strong>en</strong> son:<br />

Ä CromatografÑa realizada al gas <strong>de</strong> <strong>en</strong>trada a planta: Laboratorio realiza dos<br />

medidas diarias <strong>de</strong> cromatografÑa al gas <strong>de</strong> <strong>en</strong>trada. Se ingresa a la simulaciÄn el<br />

valor promedio <strong>de</strong> dichos resultados.<br />

Ä Porc<strong>en</strong>taje <strong>de</strong> pureza <strong>de</strong> etil<strong>en</strong>glicol <strong>de</strong> recirculaciÄn; tambián reportado por<br />

laboratorio.<br />

Ä Punto <strong>de</strong> rocÑo <strong>de</strong>l gas <strong>de</strong> <strong>en</strong>trada: medido por laboratorio dos veces al dÑa. Se<br />

ingresa a la simulaciÄn el promedio <strong>de</strong> dichos valores.<br />

La variable que es fundam<strong>en</strong>tal a la hora <strong>de</strong> obt<strong>en</strong>er un bu<strong>en</strong> resultado es la<br />

temperatura a la cual se lleva a cabo la separaciÄn; que se regula mediante el flujo <strong>de</strong> propano<br />

<strong>de</strong> <strong>en</strong>friami<strong>en</strong>to. Esta variable pres<strong>en</strong>ta constantes variaciones, razÄn por la cual no se pue<strong>de</strong><br />

dar un dato exacto, pero se ha observado que un valor razonable estÖ <strong>en</strong>tre 140,000 y 160,000<br />

MCSD. Se ajusta dicha variable <strong>en</strong> el mo<strong>de</strong>lo <strong>de</strong> HYSYS hasta obt<strong>en</strong>er resultados satisfactorios.<br />

y 45 m 3 /dÑa.<br />

Bajo las mismas consi<strong>de</strong>raciones, se consi<strong>de</strong>ra aceptable un flujo <strong>de</strong> gasolinas <strong>en</strong>tre 35<br />

Una vez realizadas estas modificaciones se pued<strong>en</strong> recopilar <strong>los</strong> resultados y<br />

contrastar<strong>los</strong> con el valor reportado por Laboratorio, principalm<strong>en</strong>te cromatografÑa y punto <strong>de</strong><br />

rocÑo <strong>de</strong>l gas residual.<br />

<strong>en</strong>trada son:<br />

Los valores reportados por Laboratorio, para cromatografÑa y punto <strong>de</strong> rocÑo <strong>de</strong>l gas <strong>de</strong>


68<br />

Fecha 24-jul 25-jul 26-jul 27-jul 28-jul<br />

Metano 91.68 91.69 91.78 91.67 91.75<br />

Etano 4.33 4.36 4.29 4.24 4.27<br />

Propano 1.50 1.51 1.46 1.43 1.48<br />

i-Butano 0.34 0.34 0.33 0.34 0.34<br />

n-Butano 0.46 0.46 0.46 0.45 0.45<br />

i-P<strong>en</strong>tano 0.16 0.16 0.16 0.16 0.16<br />

n-P<strong>en</strong>tano 0.13 0.13 0.13 0.28 0.13<br />

n-Hexano 0.28 0.25 0.25 0.23 0.28<br />

NitrÄg<strong>en</strong>o 1.02 1.00 1.03 1.12 1.04<br />

CO2 0.12 0.12 0.12 0.11 0.12<br />

Agua 0.00 0.00 0.00 0.00 0.00<br />

Dew Point [èC] -11.00 -12.00 -12.00 -13.00 -11.00<br />

Tabla 26.- Composiciones <strong>de</strong>l gas <strong>de</strong> <strong>en</strong>trada para las situaciones analizadas<br />

A continuaciÄn se pres<strong>en</strong>tan <strong>los</strong> resultados obt<strong>en</strong>idos. Se d<strong>en</strong>omina URG a la Unidad<br />

Reg<strong>en</strong>eradora <strong>de</strong> Glicol. Se <strong>en</strong>listan <strong>los</strong> datos <strong>de</strong> cromatografÑa y punto <strong>de</strong> rocÑo <strong>de</strong>l gas residual,<br />

temperatura <strong>de</strong> reg<strong>en</strong>eraciÄn <strong>de</strong> glicol y porc<strong>en</strong>taje <strong>de</strong> pureza <strong>de</strong>l glicol que ingresa al sistema<br />

<strong>de</strong> reg<strong>en</strong>eraciÄn. Para el flujo <strong>de</strong> propano necesario y la producciÄn <strong>de</strong> gasolina, se especifica si<br />

el valor obt<strong>en</strong>ido estÖ d<strong>en</strong>tro <strong>de</strong>l rango aceptable ya <strong>de</strong>finido.<br />

24/07/07 Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

Metano 92.11 92.00 0.1<br />

Etano 4.30 4.45 3.4<br />

Propano 1.44 1.43 0.7<br />

i-Butano 0.31 0.30 3.3<br />

n-Butano 0.40 0.39 2.6<br />

i-P<strong>en</strong>tano 0.12 0.11 9.1<br />

n-P<strong>en</strong>tano 0.09 0.08 12<br />

n-Hexano 0.11 0.07 57<br />

NitrÄg<strong>en</strong>o 1.03 1.05 1.9<br />

CO2 0.10 0.13 23<br />

ERROR GLOBAL [%] 0.4<br />

Dew Point Salida[èC] -21.24 -26.0 18.3<br />

Temperatura URG [èF] 237.6 242 1.8<br />

Glicol a URG[%] 69.7 68.4 1.9<br />

Gasolinas [m 3 /d] 43.27 Aceptable<br />

Flujo <strong>de</strong> propano [MCSD] 148 Aceptable<br />

Tabla 27.- SimulaciÄn <strong>de</strong> Planta Cull<strong>en</strong> <strong>de</strong>l dÑa 24/07/07


69<br />

25/07/07 Valor HYSYS Valor real<br />

Difer<strong>en</strong>cia<br />

[%]<br />

Metano 92.12 91.98 0.2<br />

Etano 4.33 4.45 2.7<br />

Propano 1.45 1.44 0.7<br />

i-Butano 0.31 0.30 3.3<br />

n-Butano 0.40 0.40 0.0<br />

i-P<strong>en</strong>tano 0.12 0.12 0.0<br />

n-P<strong>en</strong>tano 0.09 0.09 0.0<br />

n-Hexano 0.10 0.07 43<br />

NitrÄg<strong>en</strong>o 1.01 1.05 3.8<br />

CO2 0.10 0.13 23<br />

ERROR GLOBAL [%] 0.38<br />

Dew Point Salida[èC] -23.1 -24.5 1.6<br />

Temperatura URG [èF] 241.0 242 0.4<br />

Glicol a URG[%] 72.5 68.8 4.4<br />

Gasolinas [m 3 /d] 41.12 Aceptable<br />

Flujo <strong>de</strong> propano [MCSD] 149.3 Aceptable<br />

Tabla 28.- SimulaciÄn <strong>de</strong> Planta Cull<strong>en</strong> <strong>de</strong>l dÑa 25/07/07<br />

26/07/07 Valor HYSYS Valor real<br />

Difer<strong>en</strong>cia<br />

[%]<br />

Metano 92.17 92.12 0<br />

Etano 4.26 4.39 2.9<br />

Propano 1.40 1.38 2.2<br />

i-Butano 0.30 0.29 4.4<br />

n-Butano 0.40 0.38 7.1<br />

i-P<strong>en</strong>tano 0.12 0.10 22<br />

n-P<strong>en</strong>tano 0.09 0.07 30<br />

n-Hexano 0.10 0.06 84<br />

NitrÄg<strong>en</strong>o 1.04 1.09 5<br />

CO2 0.12 0.13 7.7<br />

ERROR GLOBAL [%] 0.35<br />

Dew Point Salida[èC] -22.31 -23.50 5.1<br />

Temperatura URG [èF] 238.20 242 1.6<br />

Glicol a URG[%] 70.4 68.90 5<br />

Gasolinas [m 3 /d] 39.27 Aceptable<br />

Flujo <strong>de</strong> propano [MCSD] 148 Aceptable<br />

Tabla 29.- SimulaciÄn <strong>de</strong> Planta Cull<strong>en</strong> <strong>de</strong>l dÑa 26/07/07


70<br />

27/07/07 Valor HYSYS Valor real<br />

Difer<strong>en</strong>cia<br />

[%]<br />

Metano 92.10 92.03 0.1<br />

Etano 4.21 4.40 4.3<br />

Propano 1.37 1.36 0.7<br />

i-Butano 0.31 0.30 3.3<br />

n-Butano 0.39 0.38 2.6<br />

i-P<strong>en</strong>tano 0.12 0.11 9.1<br />

n-P<strong>en</strong>tano 0.18 0.08 125<br />

n-Hexano 0.09 0.11 18<br />

NitrÄg<strong>en</strong>o 1.13 1.13 0.0<br />

CO2 0.11 0.12 8.3<br />

ERROR GLOBAL [%] 0.4<br />

Dew Point Salida[èC] -20.89 -27.5 24<br />

Temperatura URG [èF] 239.7 242 1.0<br />

Glicol a URG[%] 71.66 68.80 4.2<br />

Gasolinas [m 3 /d] 45.27 Aceptable<br />

Flujo <strong>de</strong> propano [MCSD] 146.7 Aceptable<br />

Tabla 30.- SimulaciÄn <strong>de</strong> Planta Cull<strong>en</strong> <strong>de</strong>l dÑa 27/07/07<br />

28/07/07 Valor HYSYS Valor real<br />

Difer<strong>en</strong>cia<br />

[%]<br />

Metano 92.14 92.03 0.1<br />

Etano 4.24 4.40 3.6<br />

Propano 1.42 1.40 1.4<br />

i-Butano 0.31 0.30 3.3<br />

n-Butano 0.39 0.39 0.0<br />

i-P<strong>en</strong>tano 0.12 0.12 0.0<br />

n-P<strong>en</strong>tano 0.09 0.08 12<br />

n-Hexano 0.12 0.09 33<br />

NitrÄg<strong>en</strong>o 1.05 1.09 3.7<br />

CO2 0.12 0.13 7.7<br />

ERROR GLOBAL [%] 0.4<br />

Dew Point Salida[èC] -21.1 -26.0 18.8<br />

Temperatura URG [èF] 233.2 242 3.6<br />

Glicol a URG[%] 65.6 66.00 0.6<br />

Gasolinas [m 3 /d] 40.9 Aceptable<br />

Flujo <strong>de</strong> propano [MCSD] 146.7 Aceptable<br />

Tabla 31.- SimulaciÄn <strong>de</strong> Planta Cull<strong>en</strong> <strong>de</strong>l dÑa 28/07/07


71<br />

3.5.2 VALIDACIÅN DE PLANTA POSESIÅN<br />

Durante el proceso <strong>de</strong> recolecciÄn <strong>de</strong> datos las condiciones operativas <strong>de</strong> planta se<br />

mantuvieron sumam<strong>en</strong>te estables, para las variables que interesan <strong>en</strong> la simulaciÄn. Esto<br />

significa que las sigui<strong>en</strong>tes variables, para todos <strong>los</strong> casos estudiados, se especificarÖn al valor<br />

que se pres<strong>en</strong>ta a continuaciÄn:<br />

Corri<strong>en</strong>te o equipo Variable Valor especificado<br />

Temperatura [èC] 13<br />

Gas <strong>de</strong> <strong>en</strong>trada PresiÄn [kg/cm 2 ] 70<br />

Flujo estÖndar <strong>de</strong> gas [MMMCSD] 7.0<br />

Expansor 1<br />

PresiÄn <strong>de</strong> <strong>de</strong>scarga [kg/cm 2 ] 41.5<br />

Flujo que pasa por vÖlvula JT [%] 10<br />

Expansor 2<br />

PresiÄn <strong>de</strong> <strong>de</strong>scarga [kg/cm 2 ] 26.5<br />

Flujo que pasa por vÖlvula JT [%] 1<br />

Separador V-21 Temperatura [èC] -73<br />

Columna V-5<br />

Reflujo <strong>de</strong> tope [m 3 /d] 1,000<br />

Flujo P-3 (A horno) [m 3 /d] 1,280<br />

Tabla 32.- Valores estables <strong>de</strong> operaciÄn <strong>de</strong> Planta PosesiÄn<br />

La çnica variable que se modificarÖ <strong>en</strong> cada validaciÄn es la cromatografÑa <strong>de</strong>l gas <strong>de</strong><br />

<strong>en</strong>trada a planta, y se ingresarÖ la fracciÄn <strong>de</strong> etano <strong>en</strong> Raw Product reportada por Laboratorio<br />

como especificaciÄn para la columna V-5. Los valores reportados por Laboratorio son:<br />

Fecha 31-ago 01-sep 04-sep 05-sep<br />

Metano 91.23 91.2 91.44 90.58<br />

Etano 5.33 5.24 5.04 6.00<br />

Propano 1.07 1.15 1.12 1.06<br />

i-Butano 0.24 0.26 0.24 0.25<br />

n-Butano 0.31 0.35 0.31 0.33<br />

i-P<strong>en</strong>tano 0.12 0.12 0.12 0.12<br />

n-P<strong>en</strong>tano 0.09 0.09 0.09 0.09<br />

n-Hexano 0.20 0.22 0.17 0.18<br />

NitrÄg<strong>en</strong>o 1.21 1.16 1.25 1.20<br />

CO2 0.20 0.21 0.22 0.19<br />

Tabla 33.- InformaciÄn <strong>de</strong> laboratorio para la simulaciÄn <strong>de</strong> Planta PosesiÄn<br />

Y <strong>los</strong> resultados obt<strong>en</strong>idos son:


72<br />

Corri<strong>en</strong>te o equipo Variable Valor HYSYS Valor real Difer<strong>en</strong>cia<br />

[%]<br />

Separador V-1 PresiÄn [kg/cm 2 ] 63.9 68.7 7.0<br />

PresiÄn Entrada [kg/cm 2 ] 63.9 62.5 2.2<br />

Expansor Ex-1 PresiÄn Entrada C-2 [kg/cm 2 ] 29.7 26.7 11.3<br />

PresiÄn Salida C-2 [kg/cm 2 ] 34.1 30.3 12.6<br />

Expansor Ex-2<br />

PresiÄn Entrada [kg/cm 2 ] 41.5 40.0 3.8<br />

PresiÄn Entrada C-1 [kg/cm 2 ] 26.0 23.4 11.1<br />

Entrada V-5 Temperatura [èC] -55.0 -55.0 0.0<br />

Fondo <strong>de</strong> la torre Temperatura [èC] 106.9 108.7 1.7<br />

Reflujo <strong>de</strong> fondo Temperatura [èC] 111.6 115.8 3.6<br />

Vapores <strong>de</strong> tope Temperatura [èC] -36.9 -42.0 12.1<br />

Reflujo <strong>de</strong> tope Temperatura [èC] -48.2 -52.3 7.9<br />

Gas residual Flujo [MMMCSD] 1.0 1.0 0.0<br />

Raw Product Flujo [m 3 /d] 585.8 580.0 1.0<br />

Tabla 34.- SimulaciÄn <strong>de</strong> Planta PosesiÄn dÑa 31/07/07<br />

ComposiciÄn<br />

Compon<strong>en</strong>te Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

(%mol)<br />

Metano 75.91 81.28 6.6<br />

Etano 23.30 17.82 30.8<br />

Gases <strong>de</strong> V-4 NitrÄg<strong>en</strong>o 0.29 0.23 25.1<br />

DiÄxido <strong>de</strong> carbono 0.50 0.68 27.2<br />

ERROR GLOBAL RESULTANTE [%] 11.1<br />

Metano 96.14 96.40 0.3<br />

Etano 2.26 2.27 0.6<br />

Gases <strong>de</strong> V-3<br />

Propano 0.05 0.08 39.8<br />

NitrÄg<strong>en</strong>o 1.40 0.99 41.6<br />

DiÄxido <strong>de</strong> carbono 0.15 0.28 45.4<br />

ERROR GLOBAL RESULTANTE [%] 0.8<br />

Etano 1.42 1.42 0.0<br />

Propano 51.05 49.68 2.8<br />

i-Butano 11.86 12.48 5.0<br />

Raw Product<br />

n-Butano 15.35 16.85 8.9<br />

i-P<strong>en</strong>tano 5.95 6.13 3.0<br />

n-P<strong>en</strong>tano 4.46 4.21 6.0<br />

n-Hexano 9.91 9.23 7.4<br />

ERROR GLOBAL RESULTANTE [%] 4.6<br />

Metano 93.03 93.41 0.4<br />

Etano 5.45 5.26 3.7<br />

Gas Residual<br />

Propano 0.04 0.06 30.6<br />

NitrÄg<strong>en</strong>o 1.27 0.89 42.4<br />

DiÄxido <strong>de</strong> carbono 0.20 0.38 46.2<br />

ERROR GLOBAL RESULTANTE [%] 1.1<br />

Tabla 35.- CromatografÑa obt<strong>en</strong>ida para la simulaciÄn <strong>de</strong> Planta PosesiÄn dÑa 31/07/07


73<br />

Corri<strong>en</strong>te o equipo Variable Valor HYSYS Valor real Difer<strong>en</strong>cia<br />

[%]<br />

Separador V-1 PresiÄn [kg/cm 2 ] 63.9 69.9 8.6<br />

PresiÄn Entrada [kg/cm 2 ] 63.9 59.0 8.3<br />

Expansor Ex-1 PresiÄn Entrada C-2 [kg/cm 2 ] 29.7 26.9 10.5<br />

PresiÄn Salida C-2 [kg/cm 2 ] 34.1 31.3 9.0<br />

Expansor Ex-2<br />

PresiÄn Entrada [kg/cm 2 ] 41.5 40.0 3.8<br />

PresiÄn Entrada C-1 [kg/cm 2 ] 26.0 23.3 11.5<br />

Entrada V-5 Temperatura [èC] -53.0 -53.5 0.9<br />

Fondo <strong>de</strong> la torre Temperatura [èC] 107.0 109.8 2.5<br />

Reflujo <strong>de</strong> fondo Temperatura [èC] 111.7 115.7 3.4<br />

Vapores <strong>de</strong> tope Temperatura [èC] -37.2 -38.0 2.1<br />

Reflujo <strong>de</strong> tope Temperatura [èC] -48.6 -50.2 3.2<br />

Gas residual Flujo [MMMCSD] 1.0 1.0 0<br />

Raw Product Flujo [m 3 /d] 631.4 573.0 10.2<br />

Tabla 36.- SimulaciÄn <strong>de</strong> Planta PosesiÄn dÑa 01/08/07<br />

ComposiciÄn<br />

Compon<strong>en</strong>te Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

(%mol)<br />

Metano 79.45 79.45 79.45<br />

Etano 19.50 19.50 19.50<br />

Gases <strong>de</strong> V-4 NitrÄg<strong>en</strong>o 0.41 0.41 0.41<br />

DiÄxido <strong>de</strong> carbono 0.65 0.65 0.65<br />

ERROR GLOBAL RESULTANTE [%] 7.0<br />

Metano 96.24 96.22 0.0<br />

Etano 2.21 2.42 8.9<br />

Gases <strong>de</strong> V-3<br />

Propano 0.05 0.09 43.7<br />

NitrÄg<strong>en</strong>o 1.35 0.99 36.0<br />

DiÄxido <strong>de</strong> carbono 0.16 0.28 42.8<br />

ERROR GLOBAL RESULTANTE [%] 0.7<br />

Etano 1.30 1.30 0.0<br />

Propano 50.94 51.24 0.6<br />

i-Butano 11.92 12.46 4.4<br />

Raw Product<br />

n-Butano 16.08 17.22 6.6<br />

i-P<strong>en</strong>tano 5.52 4.41 25.1<br />

n-P<strong>en</strong>tano 4.14 4.20 1.5<br />

n-Hexano 10.12 9.17 10.3<br />

ERROR GLOBAL RESULTANTE [%] 4.1<br />

Metano 93.24 92.98 0.3<br />

Etano 5.32 5.69 6.5<br />

Gas Residual<br />

Propano 0.04 0.07 40.9<br />

NitrÄg<strong>en</strong>o 1.19 0.88 34.7<br />

DiÄxido <strong>de</strong> carbono 0.21 0.38 43.6<br />

ERROR GLOBAL RESULTANTE [%] 1.1<br />

Tabla 37.- CromatografÑa obt<strong>en</strong>ida para la simulaciÄn <strong>de</strong> Planta PosesiÄn dÑa 01/08/07


74<br />

Corri<strong>en</strong>te o equipo Variable Valor HYSYS Valor real Difer<strong>en</strong>cia<br />

[%]<br />

Separador V-1 PresiÄn [kg/cm 2 ] 63.9 70.2 9.0<br />

PresiÄn Entrada [kg/cm 2 ] 63.9 63.7 0.3<br />

Expansor Ex-1 PresiÄn Entrada C-2 [kg/cm 2 ] 29.7 27.0 10.1<br />

PresiÄn Salida C-2 [kg/cm 2 ] 34.0 31.0 9.6<br />

Expansor Ex-2<br />

PresiÄn Entrada [kg/cm 2 ] 41.5 39.9 4.0<br />

PresiÄn Entrada C-1 [kg/cm 2 ] 26.0 23.9 8.7<br />

Entrada V-5 Temperatura [èC] -48.0 -48.0 0.0<br />

Fondo <strong>de</strong> la torre Temperatura [èC] 105.1 111.5 5.7<br />

Reflujo <strong>de</strong> fondo Temperatura [èC] 109.0 116.8 6.7<br />

Vapores <strong>de</strong> tope Temperatura [èC] -37.3 -37.0 0.8<br />

Reflujo <strong>de</strong> tope Temperatura [èC] -49.1 -49.5 0.9<br />

Gas residual Flujo [MMMCSD] 1.0 1.1 9.1<br />

Raw Product Flujo [m 3 /d] 585.3 575.0 1.8<br />

Tabla 38.- SimulaciÄn <strong>de</strong> Planta PosesiÄn dÑa 04/08/07<br />

ComposiciÄn<br />

Compon<strong>en</strong>te Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

(%mol)<br />

Metano 76.59 81.65 6.2<br />

Etano 22.56 17.25 30.8<br />

Gases <strong>de</strong> V-4 NitrÄg<strong>en</strong>o 0.29 0.25 17.9<br />

DiÄxido <strong>de</strong> carbono 0.56 0.85 34.7<br />

ERROR GLOBAL RESULTANTE [%] 10.7<br />

Metano 96.17 96.28 0.1<br />

Etano 2.16 2.28 5.2<br />

Gases <strong>de</strong> V-3<br />

Propano 0.05 0.75 92.9<br />

NitrÄg<strong>en</strong>o 1.44 1.03 40.0<br />

DiÄxido <strong>de</strong> carbono 0.17 0.33 48.9<br />

ERROR GLOBAL RESULTANTE [%] 1.5<br />

Etano 1.28 1.28 0.0<br />

Propano 52.97 52.13 1.6<br />

i-Butano 11.78 12.15 3.1<br />

Raw Product<br />

n-Butano 15.25 16.36 6.8<br />

i-P<strong>en</strong>tano 5.91 5.60 5.5<br />

n-P<strong>en</strong>tano 4.43 3.97 11.7<br />

n-Hexano 8.37 8.51 1.6<br />

ERROR GLOBAL RESULTANTE [%] 3.2<br />

Metano 93.35 93.28 0.1<br />

Etano 5.11 5.28 3.3<br />

Gas Residual<br />

Propano 0.05 0.06 24.0<br />

NitrÄg<strong>en</strong>o 1.28 0.92 38.7<br />

DiÄxido <strong>de</strong> carbono 0.22 0.46 51.2<br />

ERROR GLOBAL RESULTANTE [%] 0.8<br />

Tabla 39.- CromatografÑa obt<strong>en</strong>ida para la simulaciÄn <strong>de</strong> Planta PosesiÄn dÑa 04/08/07


75<br />

Corri<strong>en</strong>te o equipo Variable Valor HYSYS Valor real Difer<strong>en</strong>cia<br />

[%]<br />

Separador V-1 PresiÄn [kg/cm 2 ] 63.9 69.8 8.5<br />

PresiÄn Entrada [kg/cm 2 ] 63.9 63.3 0.9<br />

Expansor Ex-1 PresiÄn Entrada C-2 [kg/cm 2 ] 29.8 27.1 9.9<br />

PresiÄn Salida C-2 [kg/cm 2 ] 33.9 31.0 9.3<br />

Expansor Ex-2<br />

PresiÄn Entrada [kg/cm 2 ] 41.5 39.2 5.9<br />

PresiÄn Entrada C-1 [kg/cm 2 ] 26.0 24.1 7.8<br />

Entrada V-5 Temperatura [èC] -48.0 -48.0 0.0<br />

Fondo <strong>de</strong> la torre Temperatura [èC] 108.0 110.2 2.0<br />

Reflujo <strong>de</strong> fondo Temperatura [èC] 112.0 97.0 15.5<br />

Vapores <strong>de</strong> tope Temperatura [èC] -37.2 -36.2 2.7<br />

Reflujo <strong>de</strong> tope Temperatura [èC] -47.2 -49.5 4.7<br />

Gas residual Flujo [MMMCSD] 1.0 1.3 30.0<br />

Raw Product Flujo [m 3 /d] 581.3 581.0 0.0<br />

Tabla 40.- SimulaciÄn <strong>de</strong> Planta PosesiÄn dÑa 05/08/07<br />

ComposiciÄn<br />

Compon<strong>en</strong>te Valor HYSYS Valor real Difer<strong>en</strong>cia [%]<br />

(%mol)<br />

Metano 75.10 83.69 10.3<br />

Etano 24.16 15.43 56.5<br />

Gases <strong>de</strong> V-4 NitrÄg<strong>en</strong>o 0.29 0.25 16.5<br />

DiÄxido <strong>de</strong> carbono 0.45 0.65 31.0<br />

ERROR GLOBAL RESULTANTE [%] 17<br />

Metano 96.05 96.51 0.5<br />

Etano 2.35 2.10 11.8<br />

Gases <strong>de</strong> V-3<br />

Propano 0.04 0.07 40.4<br />

NitrÄg<strong>en</strong>o 1.42 1.06 33.9<br />

DiÄxido <strong>de</strong> carbono 0.14 0.27 47.9<br />

ERROR GLOBAL RESULTANTE [%] 1.2<br />

Etano 0.48 0.48 0.0<br />

Propano 51.18 50.64 1.1<br />

i-Butano 12.44 12.69 2.0<br />

Raw Product<br />

n-Butano 16.45 17.16 4.2<br />

i-P<strong>en</strong>tano 5.99 5.85 2.3<br />

n-P<strong>en</strong>tano 4.49 4.19 7.2<br />

n-Hexano 8.98 8.99 0.1<br />

ERROR GLOBAL RESULTANTE [%] 1.9<br />

Metano 92.44 93.68 1.3<br />

Etano 6.11 4.98 22.7<br />

Gas Residual<br />

Propano 0.03 0.05 30.7<br />

NitrÄg<strong>en</strong>o 1.23 0.91 34.7<br />

DiÄxido <strong>de</strong> carbono 0.19 0.38 49.1<br />

ERROR GLOBAL RESULTANTE [%] 2.9<br />

Tabla 41.- CromatografÑa obt<strong>en</strong>ida para la simulaciÄn <strong>de</strong> Planta PosesiÄn dÑa 05/08/07


76<br />

3.6 ANÉLISIS Y DISCUSIÅN DE RESULTADOS<br />

Consi<strong>de</strong>rando <strong>los</strong> resultados expuestos <strong>en</strong> la secciÄn anterior se pued<strong>en</strong> analizar <strong>los</strong><br />

resultados que <strong>en</strong>tregan las difer<strong>en</strong>tes simulaciones HYSYS realizadas para <strong>los</strong> <strong>procesos</strong> <strong>de</strong> ENAP<br />

Magallanes.<br />

Con respecto a las simulaciones referidas a <strong>los</strong> compresores recÑprocos HRA-2, HRA-5 y<br />

TLA-3 se pue<strong>de</strong> comprobar que <strong>los</strong> resultados <strong>en</strong>tregados por las simulaciones resultan<br />

satisfactorios para la precisiÄn que se requiere <strong>de</strong> una estimaciÄn computacional, las difer<strong>en</strong>cias<br />

<strong>en</strong> la mayorÑa <strong>de</strong> las variables estudiadas no superan el 5% <strong>de</strong>l valor real, lo cual se consi<strong>de</strong>ra<br />

aceptable.<br />

Para la simulaciÄn <strong>de</strong>l compresor c<strong>en</strong>trÑfugo TC-10, se han ingresado las curvas <strong>de</strong><br />

diseÜo y para la primera etapa <strong>de</strong> compresiÄn se ha aÜadido el punto <strong>de</strong> operaciÄn <strong>de</strong>finido por<br />

las condiciones que garantiza el fabricante, <strong>de</strong>bido a que durante el proceso <strong>de</strong> creaciÄn <strong>de</strong> la<br />

simulaciÄn se observaban oscilaciones apreciables y no se observaban <strong>los</strong> mismos resultados al<br />

variar las condiciones <strong>de</strong> <strong>en</strong>trada y volver al caso inicial, lo cual se <strong>de</strong>be al poco rango <strong>de</strong><br />

efici<strong>en</strong>cias <strong>en</strong> que se pres<strong>en</strong>tan las curvas. Sin embargo la segunda etapa <strong>de</strong> compresiÄn sÄlo<br />

consi<strong>de</strong>ra las curvas <strong>de</strong> diseÜo, y repres<strong>en</strong>ta eficazm<strong>en</strong>te la efici<strong>en</strong>cia <strong>de</strong> operaciÄn que<br />

garantiza el fabricante. Este factor, sumado a que el resto <strong>de</strong> variables consi<strong>de</strong>radas no difier<strong>en</strong><br />

<strong>de</strong>l valor real <strong>en</strong> mÖs <strong>de</strong> un 2%, hace concluir que la simulaciÄn repres<strong>en</strong>ta perfectam<strong>en</strong>te la<br />

situaciÄn garantizada, y que <strong>de</strong> instalarse dicho equipo e ingresar datos reales a la simulaciÄn, se<br />

obt<strong>en</strong>drÖn similitu<strong>de</strong>s satisfactorias.<br />

Para la simulaciÄn <strong>de</strong> Planta Cull<strong>en</strong>, consi<strong>de</strong>rando el proceso <strong>de</strong> validaciÄn realizado<br />

para cinco dÑas <strong>de</strong> operaciÄn <strong>de</strong> planta, se pue<strong>de</strong> concluir que el valor obt<strong>en</strong>ido para las<br />

gasolinas producidas y el flujo <strong>de</strong> propano necesario para cumplir <strong>los</strong> requerimi<strong>en</strong>tos se<br />

<strong>en</strong>cu<strong>en</strong>tran d<strong>en</strong>tro <strong>de</strong> <strong>los</strong> valores esperados, si bi<strong>en</strong> es necesario un flujo <strong>de</strong> propano<br />

relativam<strong>en</strong>te mÖs bajo que el valor real.


77<br />

El valor para el punto <strong>de</strong> rocÑo es mÖs alto que el reportado por laboratorio. Esto<br />

suce<strong>de</strong> <strong>de</strong>bido a la imposibilidad <strong>de</strong> simular la interacciÄn <strong>de</strong> glicol y agua con las<br />

termodinÖmicas disponibles <strong>en</strong> el software.<br />

Con respecto a la cromatografÑa <strong>de</strong>l gas <strong>de</strong> salida, el error global no supera el 1%, y<br />

consi<strong>de</strong>rando la variabilidad <strong>de</strong> <strong>los</strong> valores reales se pue<strong>de</strong> concluir que <strong>los</strong> resultados son<br />

repres<strong>en</strong>tativos d<strong>en</strong>tro <strong>de</strong>l rango <strong>de</strong> precisiÄn necesario. Por otra parte sÄlo se realizan dos<br />

muestreos diarios, y las variaciones <strong>en</strong>tre dichas medidas se equiparan a la variaciÄn <strong>en</strong>tre el<br />

dato real y la estimaciÄn que <strong>en</strong>trega HYSYS. De haber algçn sistema <strong>de</strong> cromatografÑa <strong>en</strong> lÑnea<br />

se podrÑa haber <strong>de</strong>purado mÖs <strong>los</strong> resultados obt<strong>en</strong>idos <strong>en</strong> este punto, pero <strong>los</strong> valores<br />

refer<strong>en</strong>ciales ya pres<strong>en</strong>tan discrepancias y por lo tanto se consi<strong>de</strong>ra cumplido el objetivo.<br />

Para la secciÄn <strong>de</strong> reg<strong>en</strong>eraciÄn <strong>de</strong> etil<strong>en</strong>glicol es don<strong>de</strong> se observan las mayores<br />

discrepancias. Interfiere directam<strong>en</strong>te <strong>en</strong> la predicciÄn <strong>de</strong>l punto <strong>de</strong> rocÑo <strong>de</strong> salida, <strong>en</strong>tregando<br />

un valor mÖs alto <strong>de</strong>l real. Esto suce<strong>de</strong> fundam<strong>en</strong>talm<strong>en</strong>te ya que se utilizÄ el paquete<br />

termodinÖmico <strong>de</strong> P<strong>en</strong>g-Robinson, y serÑa preferible un mo<strong>de</strong>lo <strong>de</strong> actividad para simular la<br />

interacciÄn <strong>en</strong>tre el glicol y el agua. Sin embargo la importancia <strong>de</strong> simular la cond<strong>en</strong>saciÄn <strong>de</strong><br />

<strong>los</strong> hidrocarburos impi<strong>de</strong> cambiar <strong>de</strong> termodinÖmica, pues se ha podido comprobar que P<strong>en</strong>g-<br />

Robinson <strong>en</strong>trega muy bu<strong>en</strong>os resultados para simular gas natural. Los fabricantes <strong>de</strong> HYSYS<br />

estÖn <strong>en</strong> conocimi<strong>en</strong>to <strong>de</strong> las fal<strong>en</strong>cias <strong>de</strong> simular glicol, agua y gas natural, y han pres<strong>en</strong>tado un<br />

paquete <strong>de</strong> actualizaciÄn especÑfico para simular dichas situaciones. Lam<strong>en</strong>tablem<strong>en</strong>te dicha<br />

actualizaciÄn se pue<strong>de</strong> implem<strong>en</strong>tar <strong>en</strong> versiones HYSYS 2004 o superiores, y ya se ha dicho que<br />

ENAP posee lic<strong>en</strong>cia sÄlo para la versiÄn 3.2 <strong>en</strong> estado estacionario, que es anterior a la versiÄn<br />

2004. En la imposibilidad <strong>de</strong> mejorar este aspecto <strong>de</strong> la simulaciÄn, se <strong>de</strong>be consi<strong>de</strong>rar esta<br />

fal<strong>en</strong>cia como un error propio <strong>de</strong>l programa y no <strong>de</strong>l mo<strong>de</strong>lo propuesto.<br />

Para la simulaciÄn <strong>de</strong> Planta PosesiÄn <strong>los</strong> resultados obt<strong>en</strong>idos <strong>en</strong> la simulaciÄn son<br />

aceptables. Las discrepancias se <strong>en</strong>cu<strong>en</strong>tran bÖsicam<strong>en</strong>te <strong>en</strong> las presiones <strong>de</strong> <strong>los</strong> compresores<br />

C-1 y C-2, pero <strong>de</strong>bido a que no se posee informaciÄn <strong>de</strong> diseÜo refer<strong>en</strong>te a estos compresores<br />

(son necesarias las curvas, principalm<strong>en</strong>te), el çnico parÖmetro disponible <strong>de</strong> modificar es la<br />

efici<strong>en</strong>cia adiabÖtica, pero esta variable <strong>en</strong> un compresor g<strong>en</strong>árico sÄlo actça sobre la


78<br />

temperatura <strong>de</strong> <strong>de</strong>scarga, y por lo tanto no se pued<strong>en</strong> ajustar mÖs <strong>los</strong> valores <strong>de</strong> presiÄn<br />

obt<strong>en</strong>idos.<br />

Si bi<strong>en</strong> el error global al simular la corri<strong>en</strong>te <strong>de</strong> tope <strong>de</strong> la columna (Gases <strong>de</strong> V-4, <strong>en</strong> las<br />

tablas) es aproximadam<strong>en</strong>te un 10%, la relevancia <strong>de</strong> dicha corri<strong>en</strong>te no es fundam<strong>en</strong>tal a la<br />

hora <strong>de</strong> juzgar el comportami<strong>en</strong>to <strong>de</strong> la simulaciÄn, pues interactça con pocos equipos. Con<br />

respecto a flujos y temperatura, <strong>en</strong> g<strong>en</strong>eral el comportami<strong>en</strong>to es satisfactorio. Consi<strong>de</strong>rando<br />

que se han simulado casos <strong>en</strong> que la composiciÄn <strong>de</strong> etano <strong>de</strong> la corri<strong>en</strong>te <strong>de</strong> Raw Product varÑa<br />

<strong>en</strong>tre 0.4% y 1.8%, y no se ha <strong>de</strong>bido modificar las especificaciones <strong>de</strong> la columna ni la<br />

estructura <strong>de</strong>l mo<strong>de</strong>lo para obt<strong>en</strong>er resultados a<strong>de</strong>cuados, se concluye que la simulaciÄn <strong>de</strong> la<br />

columna V-5 y equipos asociados resulta satisfactoria.<br />

La simulaciÄn g<strong>en</strong>eral <strong>de</strong> Planta PosesiÄn es coher<strong>en</strong>te con lo observado <strong>en</strong> terr<strong>en</strong>o, y<br />

resulta satisfactoria d<strong>en</strong>tro <strong>de</strong>l rango <strong>de</strong> precisiÄn esperado.<br />

En resum<strong>en</strong>, las simulaciones <strong>de</strong> <strong>los</strong> <strong>procesos</strong> <strong>de</strong> ENAP Magallanes <strong>en</strong> HYSYS que se han<br />

<strong>de</strong>sarrollado <strong>en</strong> el pres<strong>en</strong>te trabajo <strong>de</strong> titulaciÄn cumpl<strong>en</strong> el objetivo <strong>de</strong> dar una bu<strong>en</strong>a<br />

aproximaciÄn a la situaciÄn real, y permitirÖn al operador obt<strong>en</strong>er una estimaciÄn <strong>de</strong>l efecto que<br />

producirÑa realizar algçn cambio <strong>en</strong> las condiciones <strong>de</strong> proceso.


CAPÄTULO IV<br />

CONCLUSIONES


80<br />

IV - CONCLUSIONES<br />

situaciones:<br />

Consi<strong>de</strong>rando lo expuesto <strong>en</strong> el pres<strong>en</strong>te trabajo se pued<strong>en</strong> concluir las sigui<strong>en</strong>tes<br />

Se ha logrado implem<strong>en</strong>tar simulaciones HYSYS <strong>de</strong> <strong>los</strong> <strong>procesos</strong> <strong>productivos</strong> diarios <strong>de</strong><br />

ENAP Magallanes. Los resultados que <strong>en</strong>tregan estas simulaciones se han contrastado con <strong>los</strong><br />

valores reales <strong>de</strong> operaciÄn y se ha concluido que dichas simulaciones son efectivam<strong>en</strong>te<br />

repres<strong>en</strong>tativas <strong>de</strong> las situaciones que se <strong>de</strong>sea mo<strong>de</strong>lar, y por lo tanto servirÖn para obt<strong>en</strong>er<br />

una estimaciÄn inicial bastante certera <strong>de</strong> <strong>los</strong> efectos que suce<strong>de</strong>rÑan si se modifica algçn set<br />

point o cambian las condiciones <strong>de</strong> operaciÄn <strong>de</strong> planta.<br />

En lo referido a las simulaciones <strong>de</strong> <strong>los</strong> compresores, se concluye que HYSYS provee<br />

muy bu<strong>en</strong>os resultados para las variables mÖs importantes <strong>en</strong> estos equipos (flujos<br />

comprimidos, velocida<strong>de</strong>s <strong>de</strong> rotaciÄn, temperaturas <strong>de</strong> <strong>de</strong>scarga, etc.) y se pue<strong>de</strong> adaptar a<br />

diversas configuraciones <strong>de</strong> compresiÄn (<strong>de</strong> una o mÖs etapas, recÑprocos, c<strong>en</strong>trÑfugos, etc.), lo<br />

que permite simular difer<strong>en</strong>tes situaciones mediante una misma metodologÑa.<br />

Para la simulaciÄn <strong>de</strong> Planta Cull<strong>en</strong> se obti<strong>en</strong><strong>en</strong> resultados coher<strong>en</strong>tes y precisos para<br />

el proceso <strong>de</strong>sarrollado. La predicciÄn <strong>de</strong> cromatografÑa <strong>de</strong>l gas residual es a<strong>de</strong>cuada para<br />

obt<strong>en</strong>er estimaciones rÖpidas <strong>de</strong> la situaciÄn actual.<br />

Con respecto a la imposibilidad <strong>de</strong> repres<strong>en</strong>tar <strong>de</strong> manera mÖs fiel el proceso <strong>de</strong><br />

reg<strong>en</strong>eraciÄn <strong>de</strong> etil<strong>en</strong>glicol, se plantea la opciÄn <strong>de</strong> actualizar el software a la versiÄn HYSYS<br />

2004, que pres<strong>en</strong>ta una actualizaciÄn especÑfica para dichas situaciones. De todas formas, dicha<br />

parte <strong>de</strong>l proceso global <strong>de</strong> Planta Cull<strong>en</strong> es sumam<strong>en</strong>te estable <strong>en</strong> planta, no pres<strong>en</strong>ta mayores<br />

complejida<strong>de</strong>s <strong>de</strong> operaciÄn y se <strong>de</strong>berÑa evaluar apropiadam<strong>en</strong>te quá b<strong>en</strong>eficios reportarÑa el<br />

po<strong>de</strong>r simular a<strong>de</strong>cuadam<strong>en</strong>te dicho sistema.<br />

Para Planta PosesiÄn se pres<strong>en</strong>ta una simulaciÄn que <strong>en</strong>trega resultados confiables, y<br />

las especificaciones elegidas para la columna V-5 g<strong>en</strong>eran que sea estable fr<strong>en</strong>te a cambios <strong>en</strong>


81<br />

las condiciones <strong>de</strong> <strong>en</strong>trada a planta. De todas formas las temperaturas que <strong>en</strong>trega la<br />

simulaciÄn pres<strong>en</strong>tan ciertas difer<strong>en</strong>cias, esto <strong>de</strong>bido a que el rango <strong>de</strong> trabajo <strong>de</strong> la columna V-<br />

5 es bastante amplio (<strong>en</strong>tre -50 y 120 èC) y teÄricam<strong>en</strong>te difÑcil <strong>de</strong> simular.<br />

Los compresores y expansores <strong>de</strong> Planta PosesiÄn, <strong>de</strong>bido a la aus<strong>en</strong>cia <strong>de</strong> curvas <strong>de</strong><br />

diseÜo, se han simulado como equipos g<strong>en</strong>áricos, y esto provoca una difer<strong>en</strong>cia <strong>en</strong> las presiones<br />

que <strong>en</strong>trega HYSYS <strong>de</strong> hasta un 10% comparado con el valor real. Sin embargo para <strong>los</strong> efectos<br />

que se requiere analizar (se hace ánfasis <strong>en</strong> la predicciÄn <strong>de</strong> la corri<strong>en</strong>te <strong>de</strong> Raw Product), la<br />

presiÄn <strong>de</strong>l gas <strong>en</strong> <strong>los</strong> compresores no es un factor relevante. Una mejor precisiÄn (asociada a<br />

ingresar las curvas <strong>de</strong> operaciÄn) no afectarÖ mayorm<strong>en</strong>te la simulaciÄn global. En caso <strong>de</strong><br />

querer investigar especÑficam<strong>en</strong>te la secciÄn <strong>de</strong> expansores y compresores <strong>de</strong> Planta PosesiÄn sÑ<br />

es fundam<strong>en</strong>tal el ingreso <strong>de</strong> las curvas <strong>de</strong> operaciÄn y, como se ha visto por la simulaciÄn <strong>de</strong>l<br />

compresor c<strong>en</strong>trÑfugo TC-10, <strong>los</strong> resultados son a<strong>de</strong>cuados y sumam<strong>en</strong>te certeros.<br />

En el proceso <strong>de</strong> <strong>de</strong>sarrollo <strong>de</strong> este trabajo se ha confeccionado el ÅManual para<br />

simular <strong>procesos</strong> <strong>productivos</strong> <strong>de</strong> ENAP Magallanes <strong>en</strong> HYSYSÉ, docum<strong>en</strong>to que <strong>en</strong>trega la<br />

informaciÄn completa y <strong>de</strong>tallada <strong>de</strong>l proceso <strong>de</strong> creaciÄn <strong>de</strong> las simulaciones que se han<br />

analizado <strong>en</strong> el pres<strong>en</strong>te informe y otras. A<strong>de</strong>mÖs, <strong>de</strong> una forma clara y metodolÄgica, el manual<br />

<strong>en</strong>trega <strong>los</strong> conocimi<strong>en</strong>tos necesarios para que el usuario pueda recrear estas simulaciones y<br />

analizar <strong>los</strong> resultados, o g<strong>en</strong>erar el caso especÑfico que sea <strong>de</strong> su interás.<br />

Este manual permitirÖ difundir el uso <strong>de</strong> HYSYS <strong>en</strong>tre <strong>los</strong> profesionales <strong>de</strong> ENAP<br />

Magallanes, y servirÖ tambián <strong>de</strong> material <strong>de</strong> estudio para <strong>los</strong> cursos relacionados con<br />

simulaciones <strong>de</strong>l Departam<strong>en</strong>to <strong>de</strong> QuÑmica <strong>de</strong> la UMAG.<br />

Se pue<strong>de</strong> concluir que se han cumplido a cabalidad <strong>los</strong> objetivos planteados <strong>en</strong> el<br />

trabajo <strong>de</strong> titulaciÄn.


CAPÄTULO V<br />

BIBLIOGRAFÄA


83<br />

V - BIBLIOGRAFÄA<br />

5.1 FUENTES CONSULTADAS<br />

[1.] N. Clem<strong>en</strong>t, P. Smith, ÅHYSYS 3.2 USER GUIDEÉ, Hyprotech, a subsidiary of<br />

Asp<strong>en</strong> Technology Inc., 2003.<br />

[2.] R. Russell, ÅA FLEXIBLE AND RELIABLE METHOD SOLVES SINGLE-TOWER AND<br />

CRUDE-DISTILLATION-COLUMN PROBLEMSÉ, CHEMICAL ENGINEERING, Octubre<br />

1983, p. 53-59.<br />

[3.] E. Tarifa, ÅSIMULACIâN DE SISTEMAS QUàMICOSÉ, VersiÄn <strong>en</strong> PDF obt<strong>en</strong>ido<br />

<strong>de</strong>s<strong>de</strong> el sitio web http://www.mo<strong>de</strong>lado<strong>en</strong>ing<strong>en</strong>ieria.edu.ar.<br />

[4.] R. Treyball, OPERACIONES DE TRANSFERENCIA DE MASA, 2è Ed., Editorial Mc-<br />

Graw Hill, Mexico, 1980.<br />

5.2 INTERNET Ñ SITIOS CONSULTADOS<br />

[1.] http://support.asp<strong>en</strong>tech.com<br />

[2.] http://www.hyprotech.com<br />

[3.] http://www.<strong>en</strong>ap.cl<br />

[4.] http://www.mo<strong>de</strong>lado<strong>en</strong>ing<strong>en</strong>ieria.edu.ar/in<strong>de</strong>x.htm


ANEXOS


85<br />

ASPECTO GRÖFICO DE LAS SIMULACIONES<br />

En las sigui<strong>en</strong>tes pÖginas se pres<strong>en</strong>tan las capturas <strong>de</strong> pantalla <strong>de</strong> las simulaciones<br />

analizadas <strong>en</strong> el pres<strong>en</strong>te trabajo.<br />

Correspond<strong>en</strong> a:<br />

Ä SimulaciÄn Compresor recÑproco HRA-5.<br />

Ä SimulaciÄn Compresor recÑproco HRA-2.<br />

Ä SimulaciÄn Compresor recÑproco TLA-3.<br />

Ä SimulaciÄn Compresor c<strong>en</strong>trÑfugo TC-10.<br />

Ä SimulaciÄn Turbina <strong>de</strong>l tr<strong>en</strong> SULZER.<br />

Ä SimulaciÄn <strong>de</strong> Planta Cull<strong>en</strong>.<br />

Ä SimulaciÄn <strong>de</strong> Planta PosesiÄn.


SimulaciÉn Compresor recáproco HRA-5<br />

86


SimulaciÉn Compresor recáproco HRA-2<br />

87


SimulaciÉn Compresor recáproco TLA-3<br />

88


SimulaciÉn Compresor c<strong>en</strong>tráfugo TC-10<br />

89


SimulaciÉn Turbina Tr<strong>en</strong> Sulzer<br />

90


SimulaciÉn Planta Cull<strong>en</strong><br />

91


SimulaciÉn Planta PosesiÉn<br />

92


93<br />

INFORMACIÇN DE COMPRESORES<br />

Se pres<strong>en</strong>tan <strong>los</strong> informes <strong>de</strong> mant<strong>en</strong>imi<strong>en</strong>to predictivo para <strong>los</strong> compresores<br />

recÑprocos analizados, y las curvas <strong>de</strong> diseÜo <strong>de</strong>l compresor TC-10.<br />

Compresor HRA-2


Compresor HRA-5<br />

94


Compresor TLA-3<br />

95


Curvas compresor TC-10 - Etapa 1<br />

96


Curvas compresor TC-10 - Etapa 2<br />

97


Curvas compresor TC-10 - Pot<strong>en</strong>cia<br />

98

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!