Source: http://docplayer.es/8274029-Uso-de-algoritmos-geneticos-para-resolver-el-modelo-determinista-y-estocastico-para-el-diseno-de-una-red-de-recogida-de-residuos.html
Timestamp: 2018-03-19 05:52:41+00:00

Document:
Uso de algoritmos genéticos para resolver el modelo determinista y estocástico para el diseño de una red de recogida de residuos - PDF
Download "Uso de algoritmos genéticos para resolver el modelo determinista y estocástico para el diseño de una red de recogida de residuos"
Joaquín Salas Contreras
1 International Conference on Industrial Engineering & Industrial Management - CIO Uso de algoritmos genéticos para resolver el modelo determinista y estocástico para el diseño de una red de recogida de residuos Miguel Ortega-Mier, Joaquín Delgado Hipólito, Álvaro García-Sánchez Dpto. de Ingeniería de Organización, Administración de Empresas y Estadística. Escuela Técnica Superior de Ingenieros Industriales. Universidad Politécnica de Madrid. Calle José Gutiérrez Abascal, Madrid. Resumen Unos de los problemas que aparecen en la gestión de las cadenas de suministro inversas es el de las decisiones relativas a la localización de plantas de tratamiento, centros de recogida y de transferencia. En esta comunicación se presentan un modelo determinista MIP y el modelo estocástico derivado para ayudar a la decisión relativa a la localización de una planta de tratamiento y los centros de transferencias necesarios para reducir los costes totales de funcionamiento de un sistema de recogida de residuos. Estos modelos se intentan resolver de forma exacta con CPLEX 10.0 y dada la difi cultad se plantea un enfoque distinto utilizando algoritmos genéticos y programación lineal continua. Palabras clave: recogida residuos, logística inversa, algoritmos genéticos. 1. Introducción Uno de los problemas que aparecen en la gestión de las cadenas de suministro es la localización de las instalaciones clave de dicha cadena (plantas de producción, almacenes, etc.) En el área de la logística inversa este problema de diseño se traduce en las decisiones acerca de la localización de los diferentes centros de recogida, centros de transferencia o plantas de tratamiento. Existen muchos modelos matemáticos que intentan facilitar la toma de estas decisiones a ese respecto. Cuando las implicaciones de estas decisiones son a largo plazo es muy difícil de estimar, con precisión, la evolución de los valores de todos los parámetros y las variables de los modelos correspondientes. Esta situación es muy frecuente cuando se construyen modelos de logística inversa. Este es el caso por ejemplo de la evolución de los costes de los combustibles o la cantidad de residuo generado a lo largo de los años. En el siguiente apartado se presenta un problema de localización de instalaciones de recogida y tratamiento de residuos. En los siguientes apartados se ofrece un modelo MIP determinista que resuelve el problema previamente anunciado y el modelo estocástico asociado. Tras comentar la dificultad asociada a la resolución de dicho modelo para casos de gran tamaño, en el último apartado se presenta una metodología distinta basada en los algoritmos genéticos. 2. Presentación del problema El problema que se presenta, muy habitual cuando se quiere diseñar una red de recogida de residuos (plásticos, RSU, etc.) es el siguiente. Dados una serie de municipios en los que se genera residuo, se trata de determinar en qué municipio colocar la planta de tratamiento a la cual enviar dichos residuos y los posibles centros de transferencias donde llevar el residuo desde cada municipio que, una vez compactado, se transportar hasta la planta de tratamiento.
2 1444 Quantitative Methods El objetivo es incurrir en el menor coste posible durante un horizonte contemplado de varios periodos, pudiendo colocar la planta de tratamientos y los centros de transferencia en cualquier municipio. Se admite que en ese horizonte no tiene sentido cambiar la localización de la planta pero quizás sí el número o la localización de los centros de transferencia. Este problema surge principalmente cuando la red de recogida no existe y, por lo tanto, se tiene que diseñar el sistema de recogida partiendo de cero. La recogida del residuo se realiza con unos camiones recolectores que están preparados para recoger el residuo desde los contenedores de la calle y que transportan el residuo sin compactar al centro de transferencia correspondiente. Una vez en el centro de transferencias, el residuo se compacta y, cuando un remolque está lleno, es transferido en trailer a la planta de tratamiento. En este problema, cobra mucha importancia la relación entre los costes de recogida del camión recolector y los costes de transferencia. La relación entre estos costes de transporte y recogida, junto con los costes fijos de los centros de transferencia, llevarán a soluciones en las que se utilizan más o menos centros de transferencia. En la figura 1 (izda.) se muestra el conjunto de municipios en los que se genera residuo. Una posible solución del problema se ofrece en la misma figura (dcha); en este caso ya se ha localizado la planta de tratamiento (cruz) y se ha utilizan tres centros de transferencia (triángulos). Figura 1. Situación inicial del problema. La planta de tratamiento y los centros están por elegir (izda.) Solución del problema, la planta de tratamiento y los centros de transferencia están elegidos (dcha.) 3. El modelo de localización de una planta de tratamiento y los centros de transferencia asociados: STPNTCLP (Single Treatment Plant and Necessary Transfer Centers Location Problem) Se ha desarrollado un modelo matemático de localización dinámica asociado al problema presentado en el apartado anterior, que se ha llamado Single Treatment Plant and Necessary Transfer Centers Location Problem (STPNTCLP). Por un lado, en la literatura existen numerosos artículos que tratan de localización dinámica, (Saldanha, 1998; Antunes, 2000; Hinojosa, 2000; Bose, 2003) y; por otro, sobre modelos de logística inversa (Fleischmann, 2001). En el modelo STPNTCLP se admiten las siguientes hipótesis: El modelo es dinámico. Se contempla un horizonte dividido en varios periodos (para este problema, años). La localización de la planta de tratamiento es fija para todo el horizonte contemplado. El
3 International Conference on Industrial Engineering & Industrial Management - CIO coste fijo de la planta de tratamiento se supone independiente de dónde se localiza. Los centros de transferencia se pueden cambiar cada periodo si eso permite reducir los costes totales de la gestión. Tienen un coste fijo por periodo asociado por su uso. Los costes de apertura y cierre no son relevantes.. No se tienen en cuenta costes de recogida del residuo dentro de cada municipio, se admite que no tienen relevancia ya que pueden ser los mismos independiente de la decisión que se tome. Esta hipótesis es aceptable cuando los municipios son urbanos, y no tanto cuando un municipio esté compuesto por poblaciones o barrios dispersos. Los costes de transporte del residuo y los de transferencia son proporcionales a las distancias existentes entre el lugar de origen y destino de los mismos. Esto supone que no se tienen en cuenta rutas, lo cual se puede admitir, dado que este problema se considera estratégico. Una vez decididas las localizaciones sí se pueden calcular las rutas correspondientes.en las siguientes tablas, se puede observar cuáles son los índices, parámetros y variables utilizados en el modelo. Tabla 1. Índices del problema STPNTCLP Índices Rango i Municipios productores de residuo 1...M j Posibles centros de transferencia 1...M k Municipios candidatos para la planta de tratamiento 1...M t Periodos 1...T Tabla 2. Parámetros del problema STPNTCLP A it Producción de residuo de cada municipio en el año t (en Tm) Capacidad anual de un centro de transferencia j (en CAP j Tm) Coste medio de transferencia del residuo compactado entre un centro de ctr transferencia y la planta de tratamiento (en /km Tm) Coste medio de transporte del residuo recogido en cada municipio y llevado al cre centro de transferencia (en /km Tm) F j Costes fijos del centro de transferencia j por periodo(en ) Tabla 3. Definición de variables del problema STPNTCLP Variable z k y jt x ijt h jkt Valor 0 (1) si la planta de tratamiento está (no está) en el municipio k Tipo Binaria Valor 0 (1) si la instalación j está cerrada (abierta) en el periodo t Binaria Cantidad de residuo llevado del municipio i al centro de transferencia j el año t Cantidad de residuo compactado llevado del centro de transferencia j a la planta de tratamiento situada en el municipio k Continua Continua En la figura 2 se refleja claramente el significado de cada una de las variables utilizadas en el modelo.
4 1446 Quantitative Methods Figura 2. Significado de las variables del modelo El modelo matemático correspondiente al STPNTCLP se expresa de la siguiente forma: (1) sujeto a: (2) (3) (4) (5) (6) (7) 3.1. Solución exacta Este problema es NP-HARD (Garey, 1979) y resulta difícil resolverlo de forma exacta. Se han intentado resolver varios casos de complejidad creciente utilizando CPLEX En la tabla 4 se puede observar que, si bien para pequeños problemas se llega pronto a la solución óptima, cuando los problemas tienen mayor dimensión crece el tiempo de resolución y llega un momento en el que no se pueden resolver. Tabla 4. Tiempo de resolución y cercanía al óptimo de varios ejemplos del problema STPNTCLP 4. Modelo STPNTCLP estocástico El modelo de localización dinámica, planteado primero de forma determinista, se presenta ahora en su variante estocástica para distintos escenarios posibles.
5 International Conference on Industrial Engineering & Industrial Management - CIO En cualquier modelo estocástico existen parámetros con incertidumbre, muy usuales cuando existen decisiones de diseño cuyas implicaciones son de varios años. En este problema los parámetros con incertidumbre son: la cantidad de residuos generados en el futuro en cada municipio y el coste del combustible, que influyen en los costes de recogida y transferencia. El problema estocástico se resuelve con un enfoque bi-etápico. En la primera fase, (con incertidumbre) se obtienen los valores de las variables de decisión de diseño ( dónde y cuándo, dónde colocar la planta de tratamiento, dónde colocar los centros de transferencia y cuándo utilizarlos). El resto de variables (operativas) dependen para cada escenario (es decir de la segunda etapa). En las tres tablas siguientes se presentan los índices, parámetros y variables del modelo. Tabla 5. Índices del problema STPNTCLP estocástico Índices Rango i Municipios productores de residuo 1...M j Posibles centros de transferencia 1...M k Municipios candidatos para la planta de tratamiento 1...M t Periodos 1...T s Escenarios 1...S Existen parámetros dependientes del escenario contemplado (producción de residuo en cada municipio y los costes de recogida y transferencia). El resto se suponen independientes del escenario. Además hay que añadir un parámetro nuevo: la probabilidad asociada a cada escenario. Parámetros A it s CAP j Tabla 6. Parámetros del problema STPNTCLP estocástico Producción de residuo de cada municipio en el año t (en Tm) Capacidad anual de un centro de transferencia j (entm) ctr s Coste medio de transferencia del residuo compactado entre un centro de transferencia y la planta de tratamiento (en /km Tm) cre s Coste medio de transporte del residuo recogido en cada municipio y llevado al centro de transferencia (en /km Tm) F j Costes fijos del centro de transferencia j (en ) Probabilidad del escenario s p s Tabla 7. Definición de variables del problema STPNTCLP estocástico Variable Tipo Etapa z k Valor 0 (1) si la planta de tratamiento está (no está) Binaria en el municipio k y Primera jt Valor 0 (1) si la instalación j está cerrada (abierta) Binaria en el periodo t x ijt Cantidad de residuo llevado del municipio i al centro Continua de transferencia j en el periodo t h jkt Cantidad de residuo compactado llevado del centro Segunda de transferencia j a la planta de tratamiento situada Continua en el municipio k en el periodo t
6 1448 Quantitative Methods El modelo matemático se expresa de la siguiente forma: (8) sujeto a: (9) (10) (11) (12) (13) (14) Si el problema determinista es de difícil solución, el estocástico asociado es todavía más difícil ya que su complejidad es aún mayor. 5. Enfoque algoritmos genéticos Debido a la dificultad de resolución de las variantes deteriminista y estocástica del problema STPNTCLP (como se ha comentado en los apartados anteriores), y también a que, generalmente, para problemas reales con incertidumbre hablar de óptimo es un poco arriesgado y que es aceptable llegar a soluciones buenas, se ha planteado la resolución del problema utilizando algoritmos genéticos. Este enfoque mediante algoritmos genéticos, que se puede utilizar tanto para el problema determinista como para el estocástico del STPNTCLP, se presenta en esta comunicación aplicado al determinista (más sencillo). A su vez este mismo enfoque deriva en dos variantes. Ambas se presentan a continuación Variante 1 del enfoque mediante AG. Función de supervivencia calculada con un modelo LP En este enfoque basado en los algoritmos genéticos se van a presentar sólo aquellos aspectos específicos y más relevantes para este problema concreto (individuos, cálculo de la función de supervivencia, cruces, mutaciones, soluciones no factibles). El resto de aspectos son análogos a los que aparecen en cualquier enfoque basado en algoritmos genéticos Individuo Un individuo está compuesto de dos cromosomas (ver figura 3). El primer cromosoma (la matriz binaria) expresa si un centro de transferencia que esté en el municipio i ( fila i) se abre (1) o no (0) en el periodo t (columna t). El segundo cromosoma es un número entre (1 y N)
7 International Conference on Industrial Engineering & Industrial Management - CIO que indica en qué municipio está la planta de tratamiento (en la figura el número 23). Como se puede observar, el individuo contiene toda la información relativa a las variables de decisión de diseño (dónde ubicar las instalaciones y cuándo abrirlas). Figura 3. Individuo (izda) y método de cálculo de la función de supervivencia Cálculo de la función de supervivencia Con la información de un individuo es sencillo calcular el coste asociado a dicha solución. Habría que resolver el modelo STPNTCLP representado en las ecuaciones (1-7), pero teniendo en cuenta que las variables y y z ahora son parámetros de un modelo que ya no tiene variables jt k binarias (LP). Este modelo es de fácil y rápida resolución. Una vez obtenido el coste de la solución es sencillo calcular la función de supervivencia haciendo el inverso del coste Cruces Se han definido tres tipos de cruces, en función de la información que intercambian las soluciones: 1. Intercambio de la información relativa a los periodos (columnas). Aleatoriamente se hace un corte en las matrices de ceros y unos (primer cromosoma) y se construyen los individuos hijos a partir del segundo cromosoma y la primera parte del primer cromosoma de un padre junto con la información correspondiente a la otra mitad del primer cromosoma del otro padre. De forma análoga se construye el otro hijo. Figura 4. Cruce de dos individuos compartiendo información relativa a los periodos 2. Intercambio de la información relativa a los centros de transferencia. Ahora, análogamente al cruce anterior, el corte de las matrices se hace horizontal provocando un intercambio de
8 1450 Quantitative Methods la información relativa al uso de los centros de transferencia (filas). Figura 5. Cruce de dos individuos compartiendo información relativa a los centros de transferencia 3. Intercambio de la información relativa a la localización de la planta de tratamiento. En este caso, los padres cruzan el segundo cromosoma, sin tocar el primero. Figura 6. Cruce de dos individuos compartiendo información relativa a la localización de la planta Estos tres cruces son independientes, y sólo se puede realizar un cruce en el mismo momento. Aleatoriamente en cada instante se decide qué tipo de cruce realizar Mutaciones Las mutaciones se producen cambiando un cero por un uno (o viceversa) en el primer cromosoma; o aleatoriamente cambiando el valor del segundo cromosoma Soluciones no factibles Al realizar los cruces entre individuos y las mutaciones, que se realizan de forma aleatoria, se pueden generar soluciones no factibles. Para convertir estos nuevos individuos en factibles se tiene que implementar un procedimiento que arregle la infactibilidad abriendo centros de transferencia (cambiando ceros por unos en el primer cromosoma) Variante 2 del enfoque mediante AG. Función de supervivencia calculada con un modelo MIP Se plantea una variante al enfoque presentado en el apartado anterior. Las diferencias son las siguientes: Los individuos están formados sólo por un único cromosoma, el relativo a los centros de transferencia (matriz binaria). La función de supervivencia se calcula resolviendo un modelo asociado que en este caso es un modelo MIP en el que existen variables K variables binarias (z ) k En esta variante se busca aprovechar las ventajas relativas a la existencia de un único cromosoma ) k
9 International Conference on Industrial Engineering & Industrial Management - CIO y la simplificación del algoritmo, a consta de las desventajas relativas a la complicación del cálculo de la función de supervivencia (ahora un modelo MIP, aunque de pocas variables binarias). Figura 7. Variante 2 del enfoque AG (individuo y procedimiento de cálculo de la función de supervivencia) 6. Conclusiones En este artículo se ha presentado un problema de localización de las instalaciones de transferencia y tratamiento en un sistema de recogida de residuos. Junto al problema también se han presentado los modelo MIP deterministas y estocásticos correspondientes (STPNTCLP) que ayuda a la toma de decisiones en el problema anterior. Además se ha planteado un enfoque basado en algoritmos genéticos que intenta resolver los modelos anteriores dada la dificultad de la resolución de forma exacta. Referencias Antunes, A.; Peeters, D. (2000). A dynamic optimization model for school network planning. Socio-Economic Planning Sciences, vol. 34, pp Booker, L.B.; Goldberg, D.E.; Holland, J.H. (1989). Classifier systems and genetic algorithms. Artificial Intelligence, vol 40, Issues 1-3, pp Bose, I.; Eryarsoy, E.; He, L. (2003). Multi-period design of survivable wireless access networks under capacity constraints. Decision Support Systems. Hinojosa, Y.; Puerto, J.; Fernández, F.R. (2000). A multiperiod two-echelon multicommodity capacitated plant location problem. European Journal of Operational Research, vol. 123, pp Fleischmann, M. et al., (2001). The impact of product recovery on logistics network design. Production and Operations Management, vol 10, pp Saldanha, F.; Captivo, M.E. (1998). A heuristic approach for the discrete dynamic location problem. Location Science, vol. 6, pp Garey, M.R.; Johnson, D.S. (1979). Computers and Intractability: A Guide to the Theory of NP- Completeness. W.H. Freeman and Company.
POLITÉCNICA POLITÉCNICA
POLITÉCNICA NÚMERO ESPECIAL XI Congreso de Ingeniería de Organización International Conference on Industrial Engineering and Industrial Management D IRECCION Y O RGANIZACION POLITÉCNICA UNIVERSIDAD POLITÉCNICA
La Técnica Enumerativa de Programación Dinámica en los Problemas de Producción e Inventario y los Algoritmos Genéticos
La Técnica Enumerativa de Programación Dinámica en los Problemas de Producción e Inventario y los Algoritmos Genéticos AUTOR ING. MIGUEL JIMÉNEZ CARRIÓN M.Sc. Profesor Principal a DE, adscrito al Departamento
Algoritmo evolutivo para el problema de la ubicación de productos en las estanterías de una gran superficie
X Congreso de Ingeniería de Organización Valencia, 7 y 8 de septiembre de 2006 Algoritmo evolutivo para el problema de la ubicación de productos en las estanterías de una gran superficie José M. Albarracín
El Stroke y la Matriz de Operaciones y Materiales, nuevo enfoque para resolver el problema GMOP
5 th International Conference on Industrial Engineering and Industrial Management XV Congreso de Ingeniería de Organización Cartagena, 7 a 9 de eptiembre de 2011 El troke y la Matriz de Operaciones y Materiales,
Sujeta a: 1.1.1 Modelo Propuesto. Min (f 1, f 2 ) f 1 =
Un Modelo bi-objetivo para la Cadena de Suministro de un Sistema de Manufactura. López A. 1 Resumen. Este artículo presenta el modelo matemático de un problema de diseño de cadena de suministro directa
Escuela Técnica Superior de Ingenieros Industriales. Máster Universitario en Ingeniería Avanzada de Producción, Logística y Cadena de Suministro
Escuela Técnica Superior de Ingenieros Industriales Máster Universitario en Ingeniería Avanzada de Producción, Logística y Cadena de Suministro Tesina de Máster Resolución de modelos de planificación de
ANÁLISIS ENTRE DISTINTAS ALTERNATIVAS DE RECUPERACIÓN DE VALOR A TRAVÉS DE LA SIMULACIÓN
X Congreso de Ingeniería de Organización Valencia, 7 y 8 de septiembre de 2006 ANÁLISIS ENTRE DISTINTAS ALTERNATIVAS DE RECUPERACIÓN DE VALOR A TRAVÉS DE LA SIMULACIÓN Isabel Fernández 1, Paolo Priore
Planificación Agregada Integrada de la Empresa: una propuesta para la clasificación de problemas *
4 th International Conference on Industrial Engineering and Industrial Management XIV Congreso de Ingeniería de Organización Donostia- San Sebastián, September 8 th -10 th 2010 Planificación Agregada Integrada
Fundamentos de Matemática Aplicada. (Prácticas)
Fundamentos de Matemática Aplicada (Prácticas) Damián Ginestar Peiró UNIVERSIDAD POLITÉCNICA DE VALENCIA 1 Índice general 1. Matrices dispersas 3 1.0.1. Esquemas de almacenamiento.............. 3 1.0.2.

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 Resolución