Modeling Eco-Industrial Networks—A Representative Literature Review and Design Requirements

The development of industrial clusters has, until recently, been guided solely by economic considerations. Recent advances in measuring and assessing the environmental and social consequences of industrial production have revealed the impact of human activity on the environment and society. Eco-Industrial Networks (EINs) are being viewed as a possible solution for reducing the environmental consequences of industrial production. In this paper, we review existing approaches used for analyzing and designing EINs. Based on salient works, design requirements for models to facilitate the organization of EINs are presented and discussed.


Introduction
Increased attention is being paid to ecosystems and the consumption of resources. As a result, governments and commercial organizations are supporting shifts towards sustainable development. There are many examples of this: companies are being held liable for their environmental impact in many regions all over the world (Abdullah et al. 2019;Roberts 2004) and some countries, for instance China, have gone a step further by promoting a comprehensive legal strategy called circular economy (Yuan et al. 2006). When proposing the concept of industrial ecology, Frosch and Gallopoulos (1989) suggested a holistic approach for companies to efficiently achieve improvements in all three dimensions of sustainable development, economy, environment, and society, suggesting that industrial systems should operate like natural eco-systems (Frosch 1994;Allenby 1992;Jelinski et al. 1992). The waste and byproducts of one company could be the inputs of another (Abdullah et al. 2019). This can be done by setting up Eco-Industrial Parks (EIPs) or Networks (EINs), which involve the cooperation of companies and communities sharing and using their resources and byproducts, synergistically reducing waste. Economies of scale can help to achieve economic improvements (Tudor et al. 2007). The case of Kalundborg in Denmark is one of many promising examples (Ehrenfeld 1997;Bain et al. 2010). However, Drexhage and Murphy (2012) and Gu et al. (2013) claim that this concept exists mainly in theory and that successful cases are not set up from scratch, however have developed in response to fortunate circumstances. Although there is governmental and private interest for designing and operating industrial parks or networks, the literature does not reveal much in the way of mathematical and computational modeling for the improvement and design of corporative networks. The development of such models can help to support the process of turning the key concepts of the idea of industrial ecology by providing optimal decisions, assessing patterns, and investigating key factors in the design and operation of EINs (Gu et al. 2013). This paper aims to map the mathematical approaches used for EIN design and indicate how they may be used for designs in the future.

Literature Review and Modeling Requirements
Eco-Industrial Networks and related concepts have been identified by many different authors over the last two decades (Gu et al. 2013;Veiga and Magrini 2009;Lowe 2001;Schlarb 2001;Chertow 2000;Côté and Cohen-Rosenthal 1998;Côté and Hall 1995). The principal related definitions are: Eco-industrial park (EIP): "An eco-industrial park or estate is a community of manufacturing and service businesses located together on a common property" (Lowe 2001). Eco-industrial Networks: EIPs without geographic proximity (Ehrenfeld and Chertow 2002). The terms "industrial symbiosis networks" (Domenech and Davies 2011), and "zero waste networks" (Curran and Williams 2012) are interchangeable. The participants of a geographically spread virtual EIP are EINs. Roberts defined EINs as "networks of EIPs at national or global levels" (Roberts 2004).
An appraisal of the literature on the application of modeling methods and optimization approaches used in the field of Industrial Ecology as related to EINs has been conducted. The search was conducted with Google Scholar primarily, using terms such as: "EIN", eco industrial parks", "circular economy" with "mathematical model", "optimization" etc. From the resultant list, salient articles were identified based on the type of solution methodology used, and these have been reported in this paper. These approaches are clustered (Table 14.1) based on the method used and classified according the discussion in the previous paragraphs. Some proposed models may be part of more than one cluster. The cluster bi-level fuzzy optimization requires, for example, a fuzzy optimization, and mixed-integer linear programming (MILP) or mixed-integer non-linear programming (MINLP), which are clusters themselves. The clusters represent main practices of approaching modeling of industrial ecology. Walther et al. (2008), and Miettinen (1999) develop a negotiation algorithm for the coordination of material flow in recycling networks. The idea of industrial ecology has not been mentioned in these publications; however, based on mathematical models and an interactive negotiation algorithm, new symbiosis can be created. In  (2002), Duchin (1992) Material flow analysis Bringezu et al. (1997), Bringezu and Moriguchi (2002), Bringezu and Moriguchi (2018), Lee et al. (2006) Mixed-integer linear programming Boix et al. (2015), Gonela and Zhang (2014), Karlsson and Wolf (2007) (2017) Agent-based modeling Romero and Ruiz (2014), Romero and Ruiz (2013) order to solve the optimization model, Lagrange relaxation is applied. While the original objective function contains economic measurements, the Lagrange relaxation allows a variation of the recycling rate and thus accounts for environmental issues. Penalty functions are a common multi-objective optimization method (Miettinen 1999). Pishvaee et al. (2009) provide a meta-investigation of modeling approaches for reverse and integrated networks considering uncertainty. While the subject matter of this investigation is a sustainable reverse logistics network, environmental targets or the idea of industrial ecology have not been applied. Different methodologies have been developed in the literature in order to solve multiobjective optimization models (Walther et al. 2008). A successful application of the NIMBUS (non-differentiable interactive multi-objective bundle-based optimization system) method to optimization of eco-industrial parks has been proposed by Gu et. al. (2013); they apply the whole process of an interactive multi-objective optimization to both eco-industrial park design and optimization; this composition of computational and mathematical modeling considers multiple waste product flows. Since this tool is web-based, it is considered to have a high usability. Erol and Thöming (2005) consider multiple stakeholders with an interactive negotiation framework, neglecting social performance; uncertain behavior is also not considered. This methodology supports the improvement and design of eco-industrial parks. Li et al. (2009) consider chemical processes in general for industrial ecology; they apply TOPSIS (technique for order preference by similarity to ideal solution) and solve with an NSGA-II (nondominated sorting generic algorithm). Erol and Thöming (2005) combine the simultaneous analysis of environmental impact sensitivity (SAEIS) with multi-objective optimization performed by MINLP; they model the trade-off between economy and environment under consideration of life cycle analysis (LCA) factors. Li et al. (2009) provide the basis for this approach by illustrating the application of LCA to process optimization. The interactive surrogate worth trade-off method (ISWT) has only been applied to power plants, not to EIPs. Taskhiri et al. (2011) suggest a model to achieve a compromise among the potentially conflicting fuzzy goals of the various EIP stakeholders. Unlike the following approach, this mathematical optimization model does not consider a hierarchical structure. Karlsson and Wolf (2007) consider the hierarchy of decision-making in an eco-industrial park using a bi-level fuzzy optimization; they take into account the participating plants by means of an individual fuzzy cost goal while the upper level and overall goal of an EIP authority is the minimization of resource consumption and generation of waste. Their model thus includes environmental and economic targets but does not consider social issues.
Using a fuzzy function, lower and upper boundaries are included and provide a range in which alternative economic outcomes are acceptable for participants. The bi-level especially considers the hierarchy of stakeholders. It applies the Stackelberg Game to mathematical optimization. The basic idea of applying game theory approaches has also been investigated by Chew et al. (2009). This multi-objective bi-level optimization has also been applied to other problems such as transport planning and management problems in the past, and Kastner et al. (2015), and Yin (2002) use a nonlinear solver to find an optimal solution of an example case. The application of the max-min-concept seeks to maximize the satisfaction of the least satisfied company. Huo and Chai (2008) set up a simulation to understand evolution of industrial ecology patterns and provide new implications for design, improvement, and prediction of structural evolutions. They investigate patterns and apply evolutionary principles as well as nonlinear partial differential equations with boundary conditions and thus computationally implement interacting organisms. Evolutionary algorithms solve many nonlinear programs, and are can be used to solve multi-objective optimization problems (Zitzler and Thiele 1999). However, other than Qu et al. (2014), most of the nonlinear programs have an underlying mathematical model to be solved. Ruiz (2013, 2014) propose the application of agent-based modeling to the optimization and design of eco-industrial parks. Single companies are implemented as agents with an individual behavior and an individual economic and ecological goal (Bichraoui et al. 2013). Among other notable approaches, LCA is product based and not based on company level approach for handling eco-industrial models. However, Tong et al. (2013) applies the LCA to a system for water reuse in an industrial park. In order to determine the correct partners for increasing competitive advantage, many mathematical programming models, such as linear programming (Pan 1998), stochastic integer programming (Afshari et al. 2016), and multi-objective programming (Huang et al. 2010) have been proposed.

Analysis and Requirements
Based on the literature review, eight requirements for modeling EINs can be identified. Models should be able to include (i) Economic objectives; (ii) Environmental objectives; (iii) Social objectives; (iv) Multiple flows; (v) Multiple stakeholders; (vi) Negotiation of alternatives; (vii) Uncertainty; (viii) Ability to generate optimal solutions; and (ix) Usability. The relative importance of the requirements is dependent on the particular situation, but a comprehensive model for an EIN should address all these requirements.
The development potential of the different approaches is shown in Table 14.2 by four categories: A, B, C and D. Category A: for analysis usefulness; B: for improvements; C: for extending capabilities; and D: for design of Eco-Industrial Networks. From the table, if the goal of the model being developed relates to metrics associated with multiple flows, Input-Output analysis, Material Flow and Lagrangian Penalty based methods can be used for the purpose. Similarly, based on the development potential, Input-Output analysis, Material Flow can be used for analyzing exchanges in Eco-Industrial Networks. Likewise, Mixed Integer Programming, Agent based approaches can be used for designing new EINs. Multicriteria methods, Fuzzy approaches and Bilevel programming have been used for improving existing networks, and, Evolutionary and System Dynamics methods can be used to extend/augment existing EINs either by extending the network size or by including additional materials.

Summary
Some special patterns have been discovered in this review of mathematical approaches for modeling and designing EINs. There is no current method for achieving an optimized decision for creating new eco-industrial parks and networks. The only approach providing such an idea has been proposed by Ruiz (2013, 2014). While nearly every modeling approach that has been applied by a publication in the field of industrial ecology directly considers economic and ecological performance indicators in the objective functions, social performance has not been addressed. It is noteworthy that most of the publications only consider a single flow of material in a network; useful by-products are rarely considered explicitly by any model.
Based on the review, several design requirements for modeling EINs are determined. The domain requirements and the modeling techniques best suited for each domain are mapped. Future work on the use of these requirements for designing EINs for developing economies is ongoing.  Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license and indicate if changes were made. The images or other third party material in this chapter are included in the chapter's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.