Esta obra está bajo una Licencia de Creative Commons. Atribución - No Comercial - Compartir Igual
DOI 10.15517/ri.v29i2.32611
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica.
Análisis del control dinámico de un tren de destilación
multicomponente del área de recuperación de solvente de una planta
de polietileno lineal
Dynamic control analysis of a multicomponent distillation train of the
solvent recovery area of a linear polyethylene plant
Gustavo Enrique Ramos Montero.
Laboratorio de Ingeniería Química, Universidad del Zulia. Maracaibo, Venezuela.
Correo: geramos@ng.luz.edu.ve
Charles Albert Gutiérrez Mendoza.
Laboratorio de Ingeniería Química, Universidad del Zulia. Maracaibo, Venezuela.
Correo: chgutierrez@ng.luz.edu.ve
Beatriz María Marrufo Pachano.
Laboratorio de Ingeniería Química, Universidad del Zulia. Maracaibo, Venezuela.
Correo: bmarrufo@ng.luz.edu.ve
Recibido: 9 de marzo 2019 Aceptado: 17 de julio 2019
_________________________________________________________
Resumen
Esta investigación aporta el análisis dinámico de un tren de destilación, que en lugar del tradicional estudio
de una mezcla binaria, este se realizó para una mezcla multicomponente del área de recuperación de solvente
de una planta de polietileno lineal. Esto se llevó a cabo a través de una simulación dinámica en Aspen Plus®
con la nalidad de predecir el comportamiento de las variables operacionales (temperatura, presión, ujo, nivel
o composición) al momento de ocurrir una interrupción del ujo de alimentación al área de destilación para
modelar una parada del área de reacción, y denir los procedimientos que compensen las desviaciones en las
variables controladas. Se evaluaron varios métodos termodinámicos: Soave-Redlich-Kwong, Peng-Robinson,
Lee-Kesler-Plock, Soave-Redlich-Kwong-Boston-Mathias y Peng-Robinson-Boston-Mathias; se eligió el
método SRKBM, ya que presenta menores porcentajes de error en la mayoría de las propiedades estimadas.
Esta comparación rigurosa le añade mayor conabilidad a los valores reportados en las simulaciones. Los
resultados mostraron una desestabilización de los controladores de nivel de las columnas de destilación A
y B, además de un aumento de 14 °C aproximadamente en el perl de temperaturas de la columna A. La
columna C no presentó desviaciones de sus variables operacionales. Finalmente, para estabilizar el sistema
fue necesario realizar tres pasos que involucran el cambio a modo manual y del porcentaje de apertura de la
válvula de los controladores de ujo de vapor a la columna A (50%), de presión del separador D-101 (50%)
y de ujo de tope de la columna B (43%).
Palabras clave:
Simulación, modelaje, control de procesos, dinámica, lazos de control
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
60
Abstract
This research provides the dynamic analysis of a distillation train, which instead of the traditional
study of a binary mixture, this was carried out for a multicomponent mixture of the solvent recovery
area of a linear polyethylene plant. This was done through a dynamic simulation in Aspen Plus® with
the purpose of predicting the behavior of the operational variables (temperature, pressure, ow, level or
composition) when an interruption of the feed ow occurs in the distillation area to model a reaction area
shutdown, and then establish procedures that compensate for deviations in the controlled variables. Several
thermodynamic methods were evaluated: Soave-Redlich-Kwong, Peng-Robinson, Lee-Kesler-Plock, Peng-
Robinson-Boston-Mathias y Soave-Redlich-Kwong-Boston-Mathias; the SRKBM method was chosen,
since it presents lower error percentages in most of the estimated properties. This rigorous comparison
adds greater reliability to the values reported in the simulations. The results showed a destabilization of
the level controllers of distillation columns A and B, in addition to an increase of approximately 14 °C in
the temperature prole of column A. Column C did not show deviations from its operational variables.
Finally, to stabilize the system, it was necessary to carry out three steps involving the switch to manual
mode and changing the valve opening percentage of the column A steam ow controller (50%), the
D-101separator pressure (50%) and the column B top ow controller (43%).
Keywords:
Simulation, modeling, process control, dynamics, control loops
61
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica. DOI 10.15517/ri.v29i2.32611
1. INTRODUCCIÓN
El polietileno lineal de baja densidad (PELBD) ha sido ampliamente utilizado en la vida cotidiana,
como material de embalaje o piezas de autoviles (Kusolsongtawee y Bumroongsri, 2017) debido a
que tiene una mejor resistencia mecánica y procesabilidad en comparación con el polietileno de alta
densidad, causado por la presencia de ramicaciones cortas (Sabetzadeh, Bagheri y Masoomi, 2016).
Típicamente, el PELBD se obtiene por copolimerización de etileno con αolenas como comonómeros
con catalizadores de Ziegler-Natta (Gagieva et al., 2015). La densidad de este polímero varía entre 915
y 930 kg/m
3
mediante el uso de distintos comonómeros como lo son: 1-buteno, 1-hexeno o 1-octeno.
Este proceso en solución opera entre 160-220 ºC y 3.4-34 MPa; bajo estas condiciones el polímero
es disuelto en el solvente, comúnmente ciclohexano (Malpass, 2010).
La planta de polietileno lineal en estudio en esta investigación tiene como nalidad producir
resinas lineales de baja densidad. Para esto, el proceso cuenta con las áreas de reacción, puricación,
acabado y servicios; asimismo, utiliza una tecnología en donde el reactor es alimentado con etileno
como monómero y el comonómero (ya sea 1-buteno, 1-octeno o una mezcla de estos) con ciclohexano
como solvente.
Como en la mayoría de las plantas químicas, el área de reacción de una planta de polietileno lineal
representa el núcleo del proceso productivo, por lo que su correcto funcionamiento es fundamental.
No obstante, es necesario realizar frecuentemente paradas de sus distintos procesos. Esto trae como
consecuencia que las áreas subsecuentes que dependen del área de reacción, deban adaptarse para
seguir funcionando a pesar de suspender la operacn de alguno de sus procesos (Whiteley, 2012).
Estas secciones incluyen el área de recuperación de solvente que se encarga de separar y puricar
cada uno de los componentes una vez terminada la polimerización.
Durante este periodo, es necesario responder rápidamente a los requerimientos del proceso
mediante la manipulación de las variables, ya que esto contribuye a la eficiencia y seguridad de la
planta para poder volver a operar normalmente en el menor tiempo posible una vez restablecido
el servicio en el área de reaccn (Fabro, Arruda y Neves, 2005). Estos procedimientos de
parada pueden ser mejorados con la ayuda de la automatización y control de procesos. Las
ventajas de incluir procedimientos automáticos residen en reducir el tiempo de respuesta,
incrementa la precisión y repetitividad ya que minimiza el error humano (Leitao, P., Colombo
W. y Karnouskos S., 2015).
Actualmente, la sección de síntesis del polímero de la planta de polietileno lineal de baja
densidad cuenta con una secuencia de parada del área de reacción, este procedimiento involucra
la interrupción del suministro de 1-buteno e hidrógeno, adecuación de intercambiadores de calor,
ajuste de presión en los separadores, entre otros pasos.
Sin embargo, el área de recuperación de solvente, que incluye tres columnas de destilación, no
cuenta con un procedimiento automatizado que permita una respuesta inmediata al momento de
detener la reacción. Esto, teniendo en cuenta que las columnas de destilación son unidades complejas
que exhiben dinámicas no lineales debido a sus relaciones de equilibrio no ideales y conguraciones
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
62
de operación (Abdullah, Ahmad y Aziz, 2009); así como también, a las interacciones entre las
unidades causadas por la integración del calor y recirculación de corrientes (Xu y Bao, 2017), y
la necesidad de alcanzar alta pureza de sus productos; hace necesario la implementación de una
secuencia automática de control que permita monitorear y regular variables como temperatura,
presn, ujo, nivel y composicn en las columnas de destilación de manera segura y rápida con
la nalidad de alcanzar condiciones estables y requerimientos en las especicaciones con el menor
consumo de energía e insumos al momento de ocurrir una parada de la sección de reacción. El
arranque y paro programado de la sección de recuperación de solvente requieren de perfecta
sincronización de eventos en secuencia que ameritan el uso de computadoras e instrumentos de
medición en línea para evitar accidentes y daños en los equipos (Batres, 2013).
La simulación dinámica es ampliamente utilizada para evaluar sistemas de control y automatización,
para regular la operación de un proceso cercano a su estado estable. Además, permite predecir el
comportamiento de las variables al ocurrir alguna emergencia o falla en los equipos. Y así poder
establecer distintos procedimientos para compensar los cambios de las condiciones normales de
operacn (Luyben, 2012).
Recientes investigaciones de modelaje y simulacn dinámica de procesos de destilacn se enfocan
en mezclas binarias (Taqvi, Tufa y Muhadizir, 2016), mezclas multicomponentes (Niculescu, et al., 2017)
o trenes de destilacn (Hosgor, et al., 2014). Sin embargo, pocas investigaciones involucran trenes de
destilación de sistemas multicomponentes en estado dinámico, de allí la relevancia de esta investigación.
Debido a lo anterior, se estudia mediante el uso de Aspen Plus® y Aspen Dynamic, el área
de recuperación de solvente, el cual incluye tres columnas de destilación, cuya función es separar el
etileno, buteno, ciclohexano y oligómeros formados durante la reacción. Con el objetivo de proponer un
procedimiento que permita una respuesta inmediata al momento de presentarse una parada del área de
reacción, y sea capaz de monitorear y regular variables: temperatura, presn, ujo, nivel y composición
en las columnas de destilación de modo seguro y rápido de tal manera de alcanzar condiciones estables.
2. METODOLOGÍA
2.1Denicióndeloscomponentesdelsistema
Para iniciar el proceso de simulación en estado estacionario se tomaron como referencia los
componentes descritos en los balances de masa por diseño de plantas de producción de polietileno
lineal; estos son: etileno (C
2
H
4
) como monómero, ciclohexano (C
6
H
12
) como solvente, 1-buteno
(1-C
4
H
8
) como comonómero, 2-buteno (2-C
4
H
8
) producto de la isomerización del 1-buteno durante
la reacción y undecano (C
11
H
24
) para modelar los oligómeros.
2.2 Selección del modelo termodinámico
Una vez denidos los componentes del sistema, se tomó en cuenta la naturaleza de los mismos
y las condiciones de temperatura y presión del sistema para realizar un ltrado preliminar de los
métodos termodinámicos a seleccionar para la reproducción del comportamiento PVT del sistema
63
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica. DOI 10.15517/ri.v29i2.32611
(Carlson, 1996). Posteriormente, se llevó a cabo un estudio comparativo que permitió seleccionar
el método termodinámico que se adapte más ecientemente a las propiedades de las corrientes
y reproduzca el comportamiento PVT del sistema. En este sentido, se compararon los métodos
Soave-Redlich-Kwong (SRK), Peng-Robinson (PR), Lee-Kesler-Plock (LKP), Peng-Robinson-
Boston-Mathias (PRBM) y Soave-Redlich-Kwong-Boston-Mathias (SRKBM). Las propiedades
comparadas fueron el ujo de las fases, densidad, calor especíco y viscosidad disponibles dentro
de la documentación de la planta. Se seleccionó la corriente 5 de alimentación al D-101, ya que
esta es una corriente de mezcla vapor-líquido, Figura 1.
2.3 Descripción del proceso
La alimentación del proceso consiste en un ujo de mezcla de C
6
H
12
y C
11
H
24
(corriente 1) y
un ujo de C
2
H
4
, 1-C
4
H
8
y 2-C
4
H
8
(corriente 2) que se combinan para luego ser enviados tanto a la
columna A como al intercambiador de calor H-101 y a su vez al tanque de retención D-101. Además,
este tanque es alimentado por un ujo de isómeros de C
4
H
8
y C
6
H
12
(corriente 4) proveniente de la
zona de reacción. La primera columna (A) separa principalmente el C
2
H
4
, 1-C
4
H
8
y 2-C
4
H
8
por el
tope y C
6
H
12
y C
11
H
24
por el fondo, como alimentación lateral esta columna cuenta con un ujo de
isómeros de C
4
H
8
(corriente 9), la segunda columna (B) permite separar el C
6
H
12
del C
11
H
24
, el fondo
de esta columna se envía a la última columna (C) que sirve como extensión de la segunda columna
para terminar de separar el C
11
H
24
y recircular el C
6
H
12
, como se muestra en la Figura 1.
Figura 1. Diagrama de la simulación del área de recuperación de solvente
2.4 Desarrollo del modelo mediante Aspen Plus
Luego de denir los pasos del proceso, se procedió a realizar el modelaje de este. Las columnas
fueron formadas mediante el uso del módulo Radfrac (Luyben, 2012 y 2013), el tanque se modeló
como un separador Flash. Por su parte, para el intercambiador de calor se utilizó el modelo Heater,
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
64
Tabla 1. Especicaciones de las columnas de destilación y equipos adicionales
Columnas A B C
Etapas teóricas 40 28 30
Etapa de alimentación
Vapor: 15 y 16
Líquido: 20 y 21
Corriente 9: 1
A: 8
C: 30
B: 3
Geometría de los platos
Diámetro (mm) 2500 2300 1000
Espaciado (mm) 600 500 550
Altura del vertedero (mm) 40 55 40
Tambores de reujo
Longitud (mm) 5000 7000 3000
Diámetro (mm) 2000 3000 1000
Condensador Total
Rehervidores
Tipo Horizontal/Elíptico
Longitud (mm) 4500 5800 3600
Diámetro (mm) 1300 1600 1000
Tanque de alimentación D-101
Tipo Horizontal/Elíptico
Longitud (mm) 7000
Diámetro (mm) 3000
Tabla 2 Especicaciones de diseño de las corrientes de alimentación al área de destilación
Corriente 1 2 4 9
T (ºC) 243 243 155 37
P (MPa) 2.65 2.65 2.6 2.98
Composición kg/h %p/p kg/h %p/p kg/h %p/p kg/h %p/p
C
6
H
12
78189.9 99.62 - - 24730.8 92.39 - -
1-C
4
H
8
- - 8300.6 76.52 625.2 2.36 2086.5 98.99
2-C
4
H
8
- - 1460.9 13.47 1340.4 5.14 21.1 1.01
C
2
H
4
- - 1085.9 10.01 - - - -
C
11
H
24
298.8 0.38 - - - - - -
Total 78488.7 100 10847.4 100 26696.4 100 2107.6 100
65
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica. DOI 10.15517/ri.v29i2.32611
las bombas y válvulas fueron hechas con los módulos Pump y Valve, respectivamente. En la
Tabla 1 se presentan los datos de dimensiones y características de estos equipos. Por otro,
en la Tabla 2 se visualizan los datos de las especicaciones de las corrientes de alimenta-
ción al modelo de simulación estudiado; y en la Tabla 3 las condiciones de operación. Cada
uno de estos datos fueron obtenidos de la revisión de una planta de polietileno lineal de
baja densidad.
Tabla 3. Especicaciones de diseño de las corrientes de alimentación al área de destilación
Columna A B C
Relación de reujo (másica) 3.1 0.3 0.6
Relación D/F (másica) 0.1350 0.9319 0.9474
Separadores S-1 S-2 S-3
Fracción de separación 0.80 0.54 0.69
Intercambiador de calor H-101
Temperatura (ºC) 205
Presión (MPa) 2.5
Tanque de alimentación D-101
Temperatura (ºC) 211
Presión (MPa) 2.5
Bombas P-101 P-102
Presión de descarga (MPa) 2.76 1.2
2.5 Desarrollo del modelo de simulación en estado dinámico
Una vez denido el modelo de simulación con datos reales del sistema, se procedió a gene-
rar el dinámico del proceso. Primero, se debe especicar la geometría de los equipos (Tabla 1),
así como también se deben denir las principales válvulas de control asociadas al área de desti-
lación (Figura 1).
2.6 Migración del modelo en estado estacionario a estado dinámico con el simulador
Aspen Plus Dynamics
Al denir las válvulas necesarias para el sistema de control se realizó la migración al estado
dinámico. Se denió una exportación de la simulación controlada por presión (P Driven Dyn
Simulation) bajo las siguientes premisas: las presiones de alimentación y productos quedaron
especicadas y los ujos fueron obtenidos a través de diferencias de presión. Cada uno de los
lazos de control fue congurado dentro de la simulación dinámica.
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
66
2.7Conguracióndeloslazosdecontroldeláreaderecuperacióndesolvente
Luego, al iniciar el modelaje en estado dinámico se describieron los lazos de control del
área de recuperación de solvente, de acuerdo con su función y estructura, tomando en cuenta el
equipo al cual están asociados y su relevancia en análisis dinámico.
En la Figura 1 se muestran la ubicación de los lazos de control utilizados de acuerdo con la
distribución actual del sistema automático de control del área de recuperación de solvente de la
planta de polietileno lineal de baja densidad.
A su vez, en la Tabla 4 fueron los parámetros de los controladores y la función de cada lazo;
se utilizó el método de Ziegler-Nichols para entonar los controladores (Ziegler y Nichols, 1993),
basado en la respuesta de la variable controlada del controlador en lazo abierto luego de realizar
una perturbación de su variable manipulada.
2.8 Análisis dinámico de las variables de proceso del área de recuperación de solvente
Una vez establecida la simulación dinámica y garantizada su operación estable, se procedió a
interrumpir el suministro de C
2
H
4
e isómeros de C
4
H
8
de alimentación a la sección de recuperación de
solvente. Lo anterior se llevo a cabo cerrando completamente la válvula del controlador FIC-102, con la
nalidad de modelar una parada de la sección de reacción en donde el ujo de monómeros y comonómeros
se interrumpe y analizar el comportamiento de las variables de cada uno de los equipos. Finalmente,
cabe mencionar que se tomaron en cuenta las siguientes variables: presión, nivel, temperatura y ujo,
considerando su desviación de los valores de operación estable.
2.9 Secuencia de pasos para la adecuación de la sección de recuperación de solvente
Luego de realizar el análisis dinámico de las variables del proceso al interrumpir el suminis-
tro de C
2
H
4
y C
4
H
8
, se estudiaron los pasos que deban llevarse a cabo para poder estabilizar las
variables principales del área de recuperación de solvente. Se realizaron cambios en las variables
manipuladas (apertura de válvulas o ujo) hasta estabilizar los lazos de control, comenzando por
aquellos que se encuentren más alejados de su valor deseado.
3. DISCUSIÓN DE RESULTADOS
3.1 Selección del método termodinámico
Las corrientes en estudio están compuestas principalmente de C
6
H
12
, C
2
H
4
y C
4
H
8
, por lo que
se consideran como una mezcla no polar de compuestos reales. Según los diagramas de selección
propuestos por Aspen Technology (Carlson, 1996), para mezclas reales de compuestos no polares, se
recomienda usar los siguientes métodos termodinámicos: Soave-Redlich-Kwong (SRK), Peng-Ro-
binson (PR), Lee-Kesler-Plock (LKP), Peng-Robinson-Boston-Mathias (PRBM) y Soave-Red-
lich-Kwong-Boston-Mathias (SRKBM).
67
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica. DOI 10.15517/ri.v29i2.32611
En la Tabla 5 se presentan los resultados de la comparación de los métodos en el cálculo de las
propiedades de la corriente 5 (Figura 1). Las propiedades comparadas fueron los ujos de cada fase,
densidad, calor especíco y viscosidad dinámica.
Tabla 4 Descripción y parámetros de los controladores del modelo de simulación
Lazo Variable Variable Manipulada K
C
t
I
t
D
Acción
Controlada Variable Apertura de la válvula 1.25 2 0 Inversa
Manipulada KC tI tD
Ac-
ción
0 Inversa
FIC-101 Flujo de ciclohexano Apertura de la válvula 1.25 2 0 Inversa
FIC-102 Flujo de etileno y buteno Apertura de la válvula 1 2 0 Inversa
FIC-103 Flujo de vapor a la columna A Apertura de la válvula 3.06 2 0 Inversa
FIC-104 Flujo de buteno al D-101 Apertura de la válvula 0.88 2 0 Inversa
LIC-101 Nivel del D-101 FIC-104 (Cascada) 11.75 6.58 1.65 Inversa
PIC-101 Presión del D-101 Apertura de la válvula 55.52 1.5 0 Directa
FIC-105 Flujo de líquido a la columna A Apertura de la válvula 15.22 1.55 1.14 Inversa
LIC-103 Nivel del rehervidor columna A FIC-105 (Cascada) 19.41 1.25 2 Inversa
PIC-102 Presión de tope de la columna A Flujo de agua de enfriamiento 48.02 1.95 0 Directa
LIC-102 Nivel del tanque de reujo columna A Apertura de la válvula 3.86 8.79 2.20 Inversa
FIC-106 Flujo de fondo de la columna A Apertura de la válvula 0.94 1.98 0 Inversa
LIC-104 Nivel del tanque de reujo de la columna B FIC-106 (Cascada) 3.72 9.02 2.25 Inversa
FFIC-101 Relación de reujo, columna B Flujo de reujo 1 2 0 Inversa
PIC-103 Presión de tope de la columna B Flujo de agua de enfriamiento 8.99 2.33 0 Directa
LIC-105 Nivel del rehervidor de la columna B Flujo de vapor al rehervidor 0.29 10.61 2.65 Directa
FIC-107 Flujo de tope de la columna B Apertura de la válvula 1.36 2.30 0 Inversa
FIC-108 Flujo de fondo de la columna B Apertura de la válvula 0.91 1.98 0 Inversa
FFIC-102 Relación de reujo, columna C Flujo de reujo 0.90 1.99 0 Inversa
FIC-109 Flujo de tope de la columna C Apertura de la válvula 4.55 1.99 0 Inversa
PIC-104 Presión de tope de la columna C Flujo de agua de enfriamiento 5.64 3.70 0 Directa
LIC-106 Nivel del tanque de reujo de la columna C FIC-109 (Cascada) 39.73 18.26 4.56 Directa
LIC-107 Nivel del rehervidor de la columna C Flujo de vapor al rehervidor 6.10 7.32 1.83 Directa
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
68
Se puede observar que, de manera general, se obtuvo un ajuste satisfactorio en las propiedades
estimadas. Sin embargo, usando los modelos PR, SRK y LKP se obtuvieron errores elevados en el
cálculo de la viscosidad de la fase líquida, por lo que se descartó su uso en esta investigación. Los
métodos SRKBM y PRBM predicen ecientemente las propiedades del sistema y su uso es recomendado
para aplicaciones de renería, procesamiento de gas y procesos petroquímicos (Al Ghafri, Maitland
y Trusler, 2014). Se eligió el método SRKBM, ya que este presenta los menores porcentajes de error
en la mayoría de las propiedades estimadas, exceptuando la densidad de la fase vapor.
Tabla 5. Comparación de los modelos termodinámicos
PR SRK LKP SRKBM PRBM
Diseño Calc. %Error Calc. %Error Calc. %Error Calc. %Error Calc. %Error
Fase vapor
Flujo (kg/h)
23911.5 21877.5 8.51 23000 3.81 23421.04 2.05 24406.08 2.07 25669.65 7.35
Densidad (kg/cm
3
)
59.84 62.98 5.25 61.47 2.73 61.33 2.49 63.336 5.84 61.82 3.31
Viscosidad (cP)
0.014 0.013 3.63 0.014 0.53 0.013 3.64 0.014 0.40 0.014 0.41
Cp (kcal/kgºC)
0.723 0.605 16.32 0.608 15.99 0.626 13.42 0.605 16.32 0.607 16.04
Fase líquida
Flujo (kg/h)
67521.76 69555.9 3.01 68433.4 1.35 68012.35 0.73 67027.32 0.73 65763.75 2.60
Densidad (kg/cm
3
)
544.41 528.28 2.96 528.83 2.86 529.81 2.68 522.33 4.06 462.64 15.02
Viscosidad (cP)
0.089 0.149 67.29 0.149 66.90 0.149 67.95 0.098 10.46 0.098 10.08
Cp (kcal/kgºC)
0.744 0.769 3.37 0.779 4.75 0.707 4.95 0.767 3.04 0.777 4.43
69
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica. DOI 10.15517/ri.v29i2.32611
3.2 Reproducción del balance de masa a partir de los datos de diseño en estado
estacionario
En la Tabla 6 se muestra la comparación de los datos de diseño, con los calculados con el
simulador para las corrientes de tope y fondo de la columna A. Se observa que la temperatura y
presión calculadas para cada corriente se ajustan adecuadamente a los datos de diseño con errores
relativos menores al 2%; de igual manera para los ujos de cada componente los cuales arrojaron
errores menores al 10%.
De manera similar, la comparación entre los ujos de las corrientes de tope y fondo calculados y
de diseño de la columna B (Tabla 7), arrojó errores relativos menores al 10% para todas las variables
indicadas, concordando aceptablemente con los datos de diseño.
Comparando las corrientes de tope y fondo de la columna C (Tabla 8), se observa que el
modelo reproduce satisfactoriamente la temperatura y presión, así como también para el C
6
H
12
en
el tope y C
11
H
24
en el fondo, con errores menores al 4%. Sin embargo, el ujo de C
11
H
24
en el tope
de la columna C es mayor en los datos de diseño, que el reportado por el simulador, posiblemente
debido a las diferencias del informe publicado en el balance por diseño la cantidad de C
11
H
24
a la
salida de la columna C.
Tabla 6. Comparación de los ujos de diseño con los calculados con el simulador,
para las corrientes de tope y fondo de la columna A
Corriente Tope A Fondo A
Diseño Simulación %Error Diseño Simulación %Error
Temperatura (ºC) 53.4 53.4 0 237 237.32 0.13
Presión (MPa) 2.25 2.85 1.51 23.54 23.44 0.42
Componentes Flujo (kg/h) Flujo (kg/h)
C
6
H
12
1000.0 1010.32 1.03 104688,4 103869.00 0.78
1-C
4
H
8
11051.8 10977.70 0.67 23 20.88 9.20
2-C
4
H
8
2940.8 2927.47 0.45 29.5 28.84 2.24
C
2
H
4
1078.4 1075.10 0.31 0.0 0.0 0.0
C
11
H
24
0.0 0.0 0.0 298.8 295.80 1.00
Total 16086.11 16027.87 0.36 105163.43 104319.08 0.80
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
70
Tabla 7 Comparación de los ujos de diseño con los calculados con el simulador,
para las corrientes de tope y fondo de la columna B
Corriente Tope B Fondo B
Diseño Simulación %Error Diseño Simulación %Error
Temperatura (ºC) 181 182.38 0.76 185.9 187.9 1.07
Presión (MPa) 1.00 1.01 1.40 1.01 1.07 5.92
Componentes Flujo (kg/h) Flujo (kg/h)
C
6
H
12
104599.3 103842.0 0.72 7134.4 7119.70 0.21
1-C
4
H
8
23.0 20.88 9.20 0.0 0.0 0.0
2-C
4
H
8
29.5 28.84 2.24 0.0 0.0 0.0
C
2
H
4
0.0 0.0 0.0 0.0 0.0 0.0
C
11
H
24
0.0 0.0 0.0 304.8 297.88 2.27
Total 104694.33 103918.98 0.74 7575.10 7216.18 4.74
Tabla 8 Comparación de los ujos de diseño con los calculados con el simulador,
para las corrientes de tope y fondo de la columna C
Corriente Tope C Fondo C
Diseño Simulación %Error Diseño Simulación %Error
Temperatura (ºC) 144.4 144.4 0 247 237.22 3.96
Presión (MPa) 0.49 0.5 2.04 0.52 0.54 3.05
Componentes Flujo (kg/h) Flujo (kg/h)
C
6
H
12
7116.0 7173.12 0.80 18.4 26.57 44.41
1-C
4
H
8
0.0 0.0 0.0 0.0 0.0 0.0
2-C
4
H
8
0.0 0.0 0.0 0.0 0.0 0.0
C
2
H
4
0.0 0.0 0.0 0.0 0.0 0.0
C
11
H
24
6.0 2.30 61.61 298.8 295.57 1.08
Total 7176.60 7194.60 0.25 398.40 399.45 0.26
71
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica. DOI 10.15517/ri.v29i2.32611
3.3 Análisis de las variables de proceso del área de recuperación de solvente
Una vez realizada la simulación en estado estacionario y congurados los controladores en
el módulo dinámico, se procedió a analizar el comportamiento de los controladores y variables
del proceso, luego de interrumpir el suministro de C
2
H
4
y C
4
H
8
en la sección de recuperación de
solvente. En la Figura 2 se muestra un descenso rápido del ujo de vapor de alimentación de la
columna A, a pesar de que la válvula abre al 100% no se alcanza el valor deseado, y el ujo se
estabiliza en 14100 kg/h.
De manera similiar ocurre una caída de presión de 0.07 MPa en el D-101, por lo que la válvula
se cierra completamente para intentar aumentar la presión. Se nota que el nivel del tanque D-101
(Figura 2.c) aumenta hasta 70%, lo que causa una disminución del ujo que entra al tambor (Figura
2.d) y poder estabilizar el lazo de nivel. Este aumento se debe a que la válvula de salida del vapor
del tanque se cierra completamente.
24000
22000
20000
18000
16000
50
45
40
35
80
70
60
50
40
30
22000
20000
18000
16000
2.30
2.28
2.26
2.24
2.22
2.20
40
30
20
10
0
35000
30000
25000
20000
15000
10000
100
90
80
70
60
50
0 2 4 6 8 10 0 2 4 6 8 10 0 2 4 6 8 10 0 2 4 6 8 10
Figura 2. Respuesta de: (a) Flujo de vapor a la columna A, (b) Presión, (c) Nivel y (d) Flujo de entrada
del tanque D-101. Variable del proceso: –––, Valor deseado: ---
Por su parte en la Figura 3 se ve que el nivel del rehervidor de la columna A decae completamente,
a pesar de abrir la válvula y aumentar el ujo hasta 79500 kg/h (Figura 3.b), no fue suciente para
aumentar el nivel en el rehervidor. En cambio, la presión en el tope de la columna A (Figura 3.c)
logra estabilizarse disminuyendo el ujo de agua de enfriamiento de 5100 a 3000 kg/h, debido a la
caída de presión. Posteriormente, ocurre un aumento del ujo hasta llegar a 4000 kg/h a causa de
un aumento de temperatura, como consecuencia de la disminución del nivel del tanque de reujo
de la columna, como se muestra en la Figura 3.d.
Respecto a la Figura 4.a, se presenta el perl de temperatura de la columna A. Se puede observar
un descenso en la temperatura debido a la disminución del ujo de vapor, la cual es más evidente
en las primeras etapas de la torre. Luego ocurre un aumento de la temperatura causado por el ajuste
de la presión de la torre (Figura 3.b). Seguidamente, se aprecia un segundo aumento de temperatura
aproximadamente a las 3 h, como consecuencia del descenso del nivel del tanque de reujo (Figura
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
72
3.d), hasta estabilizarse en valores desplazados 14 ºC de las condiciones de operación normal para
las primeras 11 etapas de la torre.
Figura 3. Respuesta de: (a) Nivel del rehervidor, (b) Flujo de líquido, (c) Presión y (d) Nivel del
tanque de reujo de la columna A. Variable del $proceso: –––, Valor deseado: ---
Figura 4. Perles de temperatura y presión de la columna A
Análogamente, se presenta en la Figura 4.b el perl de presión de la columna A. Después de la
interrupción del ujo de C
2
H
4
y C
4
H
8
, dicha presión alcanza un perl adecuado rápidamente, con un
desplazamiento de aproximadamente 0.005 MPa de las condiciones normales de operación de la torre.
En relación a la columna B, se observa en la Figura 5.a que el nivel del tanque de reujo dis-
minuye lentamente hasta alcanzar el nivel mínimo. El controlador intenta contrarrestar este evento
aumentando el ujo de entrada de la columna (Figura 5.b), pero no fue posible alcanzar los valo-
res deseados en ninguno de los controladores. En la Figura 5.c se observa que se logra mantener
la relación entre ujo de reujo y el ujo de alimentación.
73
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica. DOI 10.15517/ri.v29i2.32611
Por otro lado, la presión de la columna B (Figura 5.d) logra estabilizarse en el valor deseado
al bajar el ujo de agua de enfriamiento hasta 14500 kg/h. para compensar la caída de presión. El
ujo del tope de la columna B (Figura 6.a), se estabiliza en el valor requerido hasta aproximada-
mente 4 h, tiempo en el cual el nivel del tanque de reujo (Figura 5.a) desciende por completo y el
ujo llega a bajar a 91500 kg/h, a pesar de abrir completamente la válvula. El nivel en el rehervi-
dor (Figura 6.b) presenta uctuaciones a 1h, 2.5h y 4.5 h causadas por la interrupción del ujo de
C
2
H
4
y C
4
H
8
, por la disminución del nivel del tanque de reujo de la columna A (Figura 3.d) y del
nivel del tanque de reujo de la columna B (Figura 5.a), respectivamente.
Figura 5. Respuesta de: (a) Nivel del tanque de reujo, (b) Flujo de entrada, (c) Flujo de reujo y
(d) Presión de la columna B. Variable del proceso: –––, Valor deseado: ---
Figura 6. Respuesta de: (a) Flujo de tope y (b) Nivel del rehervidor de la columna B, (c) Flujo de entrada
y (d) Nivel del tanque de reujo de la columna C. Variable del proceso: –––, Valor deseado: ---
Ahora bien, en la Figura 7.a se muestra el perl de temperaturas de la columna B luego de
interrumpir el ujo de C
2
H
4
y C
4
H
8,
se observa que este se estrecha ligeramente, con variaciones de
1ºC aproximadamente. El perl de presiones se mantiene constante con una ligera disminución de
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
74
la presión en el fondo de la columna B, como se muestra en la Figura 7.b. Tanto el perl de tem-
peratura como el de presión se consideran adecuados para operar luego del corte de C
2
H
4
y C
4
H
8
.
Figura 7. Perles de temperatura y presión de la columna B
Figura 8. Respuesta de: (a) Flujo de tope, (b) Presión (c) Temperatura de fondo y (d) Nivel del rehervidor de la
columna C. Variable del proceso:–––, Valor deseado: ---
El ujo de entrada a la columna C permanece en el valor deseado con leves cambios en la
apertura de la válvula (Figura 6.c). El nivel del tanque de reujo se mantiene igual durante todo el
periodo de evaluación (Figura 6.d). Asimismo, el ujo de tope de la columna presenta ligeras per-
turbaciones mediante la manipulación de la apertura de la válvula, gura 8.a.
Se puede observar que tanto la presión del tope (Figura 8.b) como la temperatura del fondo de
la columna C (Figura 8.c) se mantienen en el valor esperado luego del corte de C
2
H
4
y C
4
H
8
. De
igual manera, pasa lo mismo con el nivel del rehervidor de la columna C (Figura 8.d). Los perles
de temperatura y presión de la columna C se mantienen constante asntes y después de interrumpir
el ujo de C
2
H
4
y C
4
H
8
.
75
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica. DOI 10.15517/ri.v29i2.32611
3.4 Secuencia de pasos para la adecuación de la sección de recuperación de solvente
Del análisis de las variables del área de recuperación de solvente, se propusieron distin-
tos pasos para estabilizar los lazos de control y adecuar las condiciones de la sección, con la
nalidad de alcanzar rápidamente las condiciones estables una vez se logre restablecer el ujo
de C
2
H
4
y C
4
H
8
.
3.4.1 Ajuste de la presión de tanque de alimentación D-101 (PIC-101)
Este controlador en modo automático cierra la válvula completamente para aumentar la presión
en el tanque, esto trae como consecuencia que el ujo de vapor de alimentación a la columna A
disminuya (Figura 2.b). Al cambiar a modo manual con 50% de apertura de la válvula, permite que
el ujo de vapor llegue a la columna A, a pesar de la caída de la presión del tanque.
3.4.2 Ajuste del ujo de vapor de alimentación a la columna A (FIC-103)
Como se mostró en la Figura 2.a, este controlador en modo automático realiza la acción de
abrir completamente la válvula para contrarrestar la disminución del ujo de vapor causado por la
interrupción del ujo de C
2
H
4
y C
4
H
8
. En este paso se propone cambiar el modo del controlador a
manual con una apertura de la válvula del 50%, ya que con la apertura de la válvula al 100% no se
puede regular el ujo de vapor que se alimenta a la columna A.
3.4.3 Ajuste de los lazos de control de nivel
El tanque de reujo y rehervidor de la columna A, así como también, el tanque de reujo de la
columna B se vacían completamente a causa de la interrupción del ujo de C
2
H
4
y C
4
H
8
. En primer
lugar se buscó estabilizar el tanque de reujo de la columna A, manipulando el ujo líquido de ali-
mentación (FIC-105). Para esto, se cambió el modo a manual y se disminuyó la apertura de la válvula
hasta que el nivel en el tanque aumente, esto se logró con 47% de apertura de válvula del controlador
FIC-105, valores menores a este porcentaje causan la inundación de tanque de reujo de la columna A.
Una vez alcanzado el valor deseado, el controlador del nivel del tanque de reujo regula la entrada
del ujo de la corriente 9 (LIC-102) para estabilizar el nivel en 50%, como se observa en la Figura 9.a.
Esto se debe a que al disminuir el ujo de líquido que sale del tanque D-101 aumenta el ujo de vapor
que se alimenta a la columna a 19200 kg/h además de aumentar la presión del tanque a 2.26 MPa.
Similarmente, se propuso disminuir el ujo del producto de fondo de la columna A (FIC-106)
para poder estabilizar el nivel del rehervidor, esto se logró al alcanzar el 43% de apertura de la vál-
vula. Como se muestra en la Figura 10.a el nivel del rehervidor se estabiliza rápidamente en el valor
deseado, por lo que el controlador del ujo de entrada de líquido (FIC-105) es capaz de regular el nivel.
Además, se logró asegurar el nivel del tanque de reujo debido a que el controlador del ujo
de alimentación líquida a la columna A se estabiliza al 47% de apertura de la válvula como en el
caso anterior. De igual manera, se propuso disminuir el ujo de salida del tope de la columna B
(FIC-107) para ajustar el nivel del tanque de reujo de la columna B. De este estudio se obtuvo que
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
76
es necesario cerrar la válvula hasta 43% para ayuda a estabilizar el nivel (Figura 11.a), el lazo de
nivel del rehervidor de la columna B permanece estable durante este procedimiento (Figura 11.b).
Figura 9. Respuesta de: (a) Nivel del tanque de reujo de la columna A, (b) Flujo de
vapor a la columna A y (c) Presión del D-101 con 47% de apertura de la
válvula de líquido a la columna A
Figura 10. Respuesta del Nivel de (a) Rehervidor y (b)
Tanque de reujo de la columna A con 43% de apertura de la
válvula del producto de fondo de la columna A
Este paso permitió estabilizar el nivel del rehervidor y tanque de reujo de la columna A, como
se explicó anteriormente, esto se debe a que al disminuir el ujo de salida de la columna B se esta-
blece un balance de masa general del proceso que permite la operación de la columna con ujos
menores. El ajuste de los lazos de control mediante la manipulación del ujo de producto de tope
de la columna B trae como consecuencia un desplazamiento de las temperaturas de los platos a
temperaturas mayores (Figura 12.a). El perl de presión no se vio afectado debido a que los lazos
de nivel logran estabilizarse.
77
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica. DOI 10.15517/ri.v29i2.32611
Figura 11. Respuesta del Nivel de (a) Tanque de reujo y (b) Rehervidor de
la columna B con 43% de apertura de la válvula de producto de tope de la
columna B
Figura 12. Perl de temperatura de la columna A con (a) 43% apertura de la
válvula de producto de tope de la columna B y (b) 40000 kg/h de ujo de vapor
al rehervidor de la columna A
3.4.4. Ajuste del perl de temperatura de la columna A
En este paso se propuso manipular el ujo de vapor al rehervidor de la columna A luego de que
se ajustaran los lazos de nivel. Al disminuir este ujo de 52000 a 40000 kg/h, se lograron estabilizar
las temperaturas de la columna a los valores iniciales que corresponden a las temperaturas antes
de interrumpir el ujo de C
2
H
4
y C
4
H
8
. Sin embargo, esta acción trajo como consecuencia que el
nivel del tanque de reujo de la columna A se disminuya completamente. Intentar restablecer este
nivel, variando las entradas o salidas de la columna, desplazan el perl de temperatura nuevamente
a valores mayores.
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
78
3.4.5. Reajuste del nivel tanque de reujo de la columna A
El controlador de presión de la columna A (PIC-102) varía el ujo de agua de enfriamiento
del condensador para ajustar la presión del tope. Por lo tanto, al cambiar el modo de operación
a manual se obtuvo que es necesario aumentar el ujo hasta 5500 kg/h para lograr estabilizar el
lazo del nivel del tanque de reujo luego de más de 20 h de iniciada la prueba como se muestra en
la Figura 13. A pesar de estabilizar el nivel, permanece vacío durante dos horas hasta nalmente
alcanzar el valor deseado.
La manipulación de los controladores temperatura y presión de la columna A no se incluyen en
la secuencia de ajuste de los lazos de control, ya que estos pasos desestabilizan los controladores de
nivel lo que puede traer consecuencias negativas en los equipos por falta de ujo, por ejemplo, en la
bomba de reujo de la columna A. Se planteó trabajar con el perl de temperaturas de la columna
A desplazado a mayores temperaturas el cual se restablece una vez restaurado el ujo proveniente
del área de reacción. Teniendo en cuenta cada uno de los pasos descritos anteriormente se dene
como secuencia de ajuste las siguientes acciones: cambio a modo manual de los controladores de
ujo de vapor de alimentación a la columna A (FIC-103: 50% de apertura), presión del D-101 (PIC-
101: 50% de apertura) y producto de tope de la columna B (FIC-107: 43% de apertura).
60
50
40
30
20
10
0
100
90
80
70
60
50
Tiempo (h)
0 2 4 6 8 10 12 14 16 18 20 22
Figura 13. Respuesta del nivel de tanque de reujo de la columna A con
5500 kg/h de ujo de agua de enfriamiento del condensador
Es importante mencionar que el modelo desarrollado en esta investigación se basa en el sis-
tema automático de control actual de la planta de polietileno lineal, razón por la cual se recomienda
evaluar una propuesta de sistema de control como alternativa a la actual que pueda compensar más
ecientemente las perturbaciones en las variables, así como también un estudio anexo que involu-
cre el área de reacción y el procedimiento de parada de reacción.
79
Ingeniería 29 (2): 59-80, julio-diciembre, 2019. ISSN: 2215-2652. San José, Costa Rica. DOI 10.15517/ri.v29i2.32611
4. CONCLUSIONES
La columna A es la más afectada al realizar la interrupción del ujo de C
2
H
4
y C
4
H
8
, pues
se desestabiliza el control del nivel del rehervidor y condensador. El perl de temperatura
aumenta 14 ºC aproximadamente en las primeras etapas de la torre.
En la columna B, solo se ve afectado el nivel del tanque de reujo, mientras que la
columna C logra estabilizarse por sola. Se propusieron diferentes acciones para estabilizar
los controladores del área de recuperación de solvente: tanto el controlador de ujo de vapor
de alimentación a la columna A como el de presión del D-101 se cambian a modo manual con
apertura de la válvula en 50%.
Para estabilizar los controladores de nivel se debe cambiar a modo manual el controlador
del producto de tope de la columna B con 43% de apertura de válvula. Ajustar el perl de
temperaturas de la columna A, manipulando el ujo de vapor del rehervidor, desestabiliza
los lazos de nivel, estos pueden ajustarse nuevamente manipulando el controlador el ujo de
agua de enfriamiento del condensador, pero es un proceso lento que representaría problemas
operacionales en la columna A.
REFERENCIAS
Abdullah, Z., Ahmad, Z. y Aziz, N. (Agosto 2009). MIMO Neural Network Model for Pilot Plant Distillation
Column. En 10th International Symposium on Process Systems Engineering, Salvador de Bahía, Brasil.
Al Ghafri, S. Z., Maitland, G. C. y Trusler, J. M. (2014). Experimental and modeling study of the phase
behavior of synthetic crude oil + CO
2
. Fluid Phase Equilibria, 365, 20-40.
Batres, R. (2013). Simulation-based Planning of Shutdown Operations. Procedia Computer Science, 22,
1294-1302.
Carlson, E. (1996). Don’t gamble with physical properties for simulations. Chemical Enginering Progress,
20, 36-46.
Fabro, J., Arruda, L. y Neves, F. (2005). Startup of a distillation column using intelligent control techniques.
Computers & Chemical Engineering, 30, 309-230.
Gagieva, S. C., Kurmaev, D. A., Tuskaev, V. A., Zubkevich, S. V., Borissova, A. O., Fedyanin, I. V. et al.
(2015). Preparation of linear low-density polyethylene from ethylene by tandem catalysis of titanium(IV)
and nickel(II) non-metallocene catalysts. Polyhedron, 98, 131-136.
Hosgor, E., Kucuk, T., Oksal, I. N. y Kaymak, D. B. (2014). Design and control of distillation processes for
methanol–chloroform separation. Computers & Chemical Engineering, 67, 166-177.
Kusolsongtawee, T. y Bumroongsri, P. (2017). Optimización of energy consumption in gas-phase polyme-
rization process for linear low density polyethylene production. Energy Procedia, 138, 772-777.
Leitao, P., Colombo W. y Karnouskos S. (2015). Industrial automation based on cyber-physical systems
technologies: Prototype implementations and challenges. Computers in Industry, 35 (8), 11-25.
RAMOS, GUTIÉRREZ Y MARRUFO: Análisis del control dinámico de un tren de destilación...
80
Luyben, W. (2012). Use of dynamic simulation for reactor safety analysis. Computers & Chemical Engi-
neering, 40, 97-109.
Luyben, W. (2013). Distillation design and control using AspenTM Simulation. Nueva York: John Wiley &
Sons, Inc.
Malpass, D. (2010). Introduction to Industrial Polyethylene: Properties, Catalysts, and Processes. Nueva
York: John Wiley & Sons, Inc.
Niculescu, A., Constantin, T., Ana, G. y Draghia, M. (2017). Dynamic simulation of a multicomponent dis-
tillation column for D-T separation. Fusion Engineering and Design, 124, 752-756.
Sabetzadeh, M., Bagheri, R. y Masoomi, M. (2016). Effect of nanoclay on the properties of low density
polyethylene/linear low density polyethylene/thermoplastic starch blend lms. Carbohydrate Polymers,
141, 75-81.
Taqvi, S. A., Tufa, L. D. y Muhadizir, S. (2016). Optimization and Dynamics of Distillation Column Using
Aspen Plus®. Procedia Engineering, 148, 978-984.
Whiteley, K. (2012). Ullmann’s Encyclopedia of Industrial Chemistry. Nueva York: John Wiley & Sons, Inc.
Xu, S. y Bao, J. (2017). Distributed control of plantwide chemical processes. Journal of Process Control,
19 (10), 1671-1687.
Ziegler, J. G. y Nichols, N. B. (1993). Optimum settings for automatic controllers. Journal of Dynamic Sys-
tems, Measurement, and Control, 115, 220-222.