A Hybrid Multi-Objective Optimization Framework for Preliminary Process Design Based on Health, Safety and Environmental Impact

: Due to increasingly stringent legal requirements and escalating environmental control costs, chemical industries have paid close attention to sustainable development without compromising their economic performance. Thus, chemical industries are in need of systematic tools to conduct sustainability assessments of their process / plant design. In order to avoid making costly retroﬁts at later stages, assessments during the preliminary design stage should be performed. In this paper, a systematic framework is presented for chemical processes at the preliminary design stage. Gross proﬁt, Health Quotient Index (HQI), Inherent Safety Index (ISI) and the Waste Reduction (WAR) algorithm are used to assess the economic performance, health, safety and environmental impact of the process, respectively. The fuzzy optimization approach is used to analyse the trade-o ﬀ among the four aspects simultaneously, as they often conﬂict with each other. Deviation between the solution obtained from mathematical optimization model and process simulator is determined to ensure the validity of the model. To demonstrate the proposed framework, a case study on 1, 4-butanediol production is presented.


Introduction
Process design is a core element in the field of chemical engineering.It can be considered a centre point, bringing together all chemical engineering components as a whole.Process design is associated with creating processes or improving existing processes.An integral part of process design is process synthesis.Process synthesis is defined as "the discrete decision-making activities of conjecturing (1) which of the many available component parts one should use, and (2) how they should be interconnected to structure the optimal solution to a given design problem."[1].The field of process synthesis has seen significant developments since its inception in the 1960s led by the late Roger W. H. Sargent [2].Most notably, key pioneering contributions have been an integral part of establishing what process synthesis is and what it entails.For instance, Nishida et al. [3] provided an important overview of developments within the boundary of process synthesis and defined process synthesis as "an act of determining the optimal interconnection of processing units as well as the optimal type and design of the units within a process system".Meanwhile, the Onion model [4] was proposed as a systematic overview and guide to process synthesis thinking.The Onion model emphasized that a reactor is designed first followed by separation and recycle streams, heat recovery systems and utility systems [4].Foo and Ng [5] then extended the Onion model further by incorporating material recovery and treatment systems (see Figure 1).On the other hand, Douglas' model is another well-accepted decision hierarchy approach that was proposed for process synthesis in the late 80s [6].
Processes 2019, 7, x FOR PEER REVIEW 2 of 19 proposed as a systematic overview and guide to process synthesis thinking.The Onion model emphasized that a reactor is designed first followed by separation and recycle streams, heat recovery systems and utility systems [4].Foo and Ng [5] then extended the Onion model further by incorporating material recovery and treatment systems (see Figure 1).On the other hand, Douglas' model is another well-accepted decision hierarchy approach that was proposed for process synthesis in the late 80s [6].To date, there has been a vast number of developed process synthesis methodologies.Often, these methodologies pay close attention to the technical and economic performance of a process design.Meanwhile, other aspects such as safety, occupational health and environmental impacts are typically considered at the later/final stages of design [7,8].However, several papers argue that these aspects should be considered at the preliminary stage of process synthesis as the cost of process improvement and operational risks can be significantly lowered compared to at the later stages (detailed design) [9,10].In this respect, several alternative methods for evaluating occupational health, safety and environmental impacts have been presented.
In the area of safety evaluations, Edwards and Lawrence [11] developed the earliest method for assessing inherent safety in a given process.This method is known as Prototype Index of Inherent Safety (PIIS).PIIS ranks the inherent safety level of alternative chemical process routes based on main reactions and parameters such as pressure, temperature, yield, heat of reaction, inventory, flammability, toxicity and explosiveness.Despite considering the mentioned parameters, PIIS focuses solely on the main reaction of a process and not the other parts of the process.In view of this, Heikkila [12] proposed an alternative approach called the Inherent Safety Index (ISI).ISI considers the same factors as in the PIIS along with additional ones such as corrosions, side reactions, and To date, there has been a vast number of developed process synthesis methodologies.Often, these methodologies pay close attention to the technical and economic performance of a process design.Meanwhile, other aspects such as safety, occupational health and environmental impacts are typically considered at the later/final stages of design [7,8].However, several papers argue that these aspects should be considered at the preliminary stage of process synthesis as the cost of process improvement and operational risks can be significantly lowered compared to at the later stages (detailed design) [9,10].In this respect, several alternative methods for evaluating occupational health, safety and environmental impacts have been presented.
In the area of safety evaluations, Edwards and Lawrence [11] developed the earliest method for assessing inherent safety in a given process.This method is known as Prototype Index of Inherent Safety (PIIS).PIIS ranks the inherent safety level of alternative chemical process routes based on main reactions and parameters such as pressure, temperature, yield, heat of reaction, inventory, flammability, toxicity and explosiveness.Despite considering the mentioned parameters, PIIS focuses solely on the main reaction of a process and not the other parts of the process.In view of this, Heikkila [12] proposed an alternative approach called the Inherent Safety Index (ISI).ISI considers the same factors as in the PIIS along with additional ones such as corrosions, side reactions, and inventory for both inside and outside the battery limits, types of equipment and process structure.Later, Palaniappan et al. [13] added a few supplementary indexes such as worst reaction index, overall chemical index and total chemical index to the previously proposed ISI to form i-Safe.
Aside from safety, various methods have been developed to evaluate occupational health.For example, Hassim and Edwards [14] developed Process Route Healthiness Index (PRHI) to rank the process alternatives based on the potential of working activities and process conditions that may harm workers.However, PRHI is not suitable for assessment in the preliminary design stage, as it requires complete process information (e.g., points for manual handling etc.).In this respect, Hassim et al. developed the Inherent Occupational Health Index (IOHI) [10], which is suitable for application at the initial stage of the process research and development (R&D).This work was then extended to a more detailed assessment called the Health Quotient Index (HQI), which uses detailed process information available in the preliminary design stage [15].HQI is able to rank the process alternatives based on health risk values from the fugitive emissions and to calculate the risk of a given process.
Apart from occupational health, there are many well-established approaches reported in the literature on environmental impact assessments.The Environmental Hazard Index (EHI) [16] and Atmospheric Hazard Index (AHI) [17] are among the early methods that assess inherent environmental performance of chemical process routes.Later, Gunasekera and Edwards [18] integrated EHI with AHI into a new method which is known as Inherent Environmental Toxicity Hazard (IETH).IETH can estimate the inherent environmental friendliness of a chemical plant in all media including air, soil and aquatic due to total loss of contaminants.Cabezas et al. [19] then introduced the generalized Waste Reduction (WAR) algorithm.This algorithm determines the total environmental impacts based on the Potential Environmental Index (PEI) balance, which assigned environmental impact values to different pollutants.Young and Cabezas [20] extended PEI balance to integrate the energy consumption of the chemical process into the environmental evaluation.Similarly, Andiappan et al. [21] presented the incremental environmental burden assessment (IEBA), which used the concept of economic potential assessment method developed by Seider et al. [22] to determine the new environmental burden for each process.
It is essential to note that the abovementioned evaluation methods focus only on assessing the performance of a process design based on just one aspect (e.g., safety or occupational health or environmental impact) and not addressing them simultaneously.In this respect, there have been attempts to develop process synthesis approaches that address safety, occupational health and environmental impacts simultaneously.For instance, Azapagic et al. [8] and Othman et al. [23] presented sustainable assessment and design selection approaches that consider economic, environmental, and social aspects.Besides, Al-Sharrah et al. [24] presented a multi-objective optimization model which considers environmental impact, economic performance and operational risk simultaneously for the petrochemical industry.Liew et al. [25] developed a systematic approach to the screening of sustainable chemical reaction pathways at the research and development (R&D) stage.Based on the proposed approach, fuzzy optimization is used to trade-off economic performance, health, safety and environmental impact.Similarly, Ng et al. [26] presented a multi-objective process synthesis approach which considered economic performance, health, safety and environmental impact for biorefineries.Following this, a visualization tool called the Piper diagram was also proposed for considering economic, safety and environmental aspects simultaneously [27].However, due to its graphical nature, the Piper diagram is unable to consider process optimization.
It is clear that the abovementioned work presents approaches to evaluate safety, occupational health and environmental impacts at the preliminary stage of design.However, these approaches focus on preliminary screening without considering operating conditions in a given process.The choice of process operating conditions has a process direct impact on screening decisions and this should be given closer attention.To address this, process simulators can be used to analyse operating conditions for a given process.Process simulators can be used to represent chemical processes in terms of mathematical models and solve them to attain insights on their performance [28].Process simulation is an essential and complementary task of process synthesis, as it predicts how a process design would behave under defined operating conditions [29,30].Despite the benefits of process simulators, they suffer from various limitations.For instance, simulators are limited to considering a single objective (e.g., product yield or production rate) at a given run.This limitation disables simulators from finding a trade-off between multiple aspects such as process economics, occupational health, safety and environmental impacts for a given process design.In order to consider multiple aspects simultaneously, multi-objective optimization is required.
As such, this paper presents a systematic framework which combines the benefits from both process simulation and multi-objective optimization to address economic, environmental, health and safety aspects simultaneously at the preliminary design stage.This paper is structured as follows: First, a formal problem statement is given in Section 2. A systematic framework for preliminary process design is presented in Section 3; To illustrate the proposed framework, a case study on the screening of 1,4-butanediol production processes is presented in Section 4; and the work is finally concluded in Section 5.

Problem Statement
The problem addressed in this work is stated as follows; given a set of alternative processes k to produce chemical product (in a given output stream) p.Each alternative process k has a set of unit operation j with operating capacities, x j,k and a set of (output) stream p to or from unit operation j represented by matrix a p,j,k .Alongside this, each alternative process k differs in economic performance, health, safety and environment impacts.In this work, four objectives are considered for performing a preliminary evaluation and screening of alternative process k.The economic performance is determined based on the gross profit of alternative process k (GP).The health impacts of process alternative k are evaluated using Health Quotient Index (HQI).Meanwhile, the safety and environmental impact of each alternative process k are determined via Inherent Safety Index (ISI) and Potential Environmental Index (WAR) respectively.
This work combines the advantages of process simulation, input-output modelling (IOM) and fuzzy optimization to determine the trade-off between the GP, HQI, ISI and WAR objectives.In fuzzy optimization, the degree of satisfaction (λ k ) is introduced to quantify the degree of satisfaction for four objectives in each alternative process k.Following this, the alternative process k with the highest λ k will be selected.A detailed description of the proposed framework to address the stated problem is presented in the following section.

Framework for Preliminary Process Design
Figure 2 presents the proposed framework for preliminary process design, which considers economic performance, health, safety and environment impacts simultaneously.As shown, the proposed framework begins by identifying the product i that is expected to be produced.Next, alternative process k which produce product i is determined.In order to analyse the performance of each alternative process k, process simulation tools (e.g., Aspen HYSYS, SuperPro Designer, PRO/II, etc.) can be used to simulate the process.Based on the simulation results, process data such as mass flow rates and energy requirements for each process unit j can be determined.These data are then used to develop an input-output model for each process.Input-output model was presented by Leontief [31,32] to analyse the relationship among the raw materials requirement, goods production and the exchange of materials within different economic sectors.Tan et al. [33] extended the usage of input-output model in life cycle assessment and ecological footprint analysis.Aviso et al. [34] integrated the approach with optimization framework to analyse the eco-industrial supply chain under a water footprint constraint.Recently, the input-output model has been used to describe the material and energy balance of processes in a system [35].Based on the approach [35], the optimal operational adjustment in multi-functional energy systems in response to process inoperability can be determined.Kasivisvanathan et al. [36] also used the input-output model with robust optimization for process synthesis and design of multi-functional energy systems with uncertainties.Most recently, Foong et al. [37] extended the use of the input-output model for sustainable oil palm plantation development.A detailed description on the procedure to develop an input-output model is presented with an example in Andiappan et al. [38].The input-output model for each alternative process k, given as in Equation ( 1) Processes 2019, 7, x FOR PEER REVIEW 5 of 19 plantation development.A detailed description on the procedure to develop an input-output model is presented with an example in Andiappan et al. [38].The input-output model for each alternative process k, given as in Equation ( 1) In Equation ( 1), ap,j,k represents the matrix of input or output mass flow rate of stream p to/from unit operation j in process k. xj,k represents the capacity or size of the unit operation j in process k and yp,k is the net output of stream p in process k.A positive value of yp,k indicates that it is an output stream which is either product or effluent stream whereas a negative value indicates that it is an input stream and a value of zero denotes that it is an intermediate stream.In Equation (1), a p,j,k represents the matrix of input or output mass flow rate of stream p to/from unit operation j in process k. x j,k represents the capacity or size of the unit operation j in process k and y p,k is the net output of stream p in process k.A positive value of y p,k indicates that it is an output stream which is either product or effluent stream whereas a negative value indicates that it is an input stream and a value of zero denotes that it is an intermediate stream.
Meanwhile, the lower and upper limits of the capacity of unit operation j in process k are presented in Equation ( 2 In addition to Equations ( 1) and ( 2), the economic, environmental, health and safety impact for each process k is determined.Firstly, in this work, economic performance for each process k is measured via gross profit as shown in Equation (3): where y Prod p,k and y Raw p,k are the mass flow rate of the products and raw materials in process k respectively; C Prod p,k and C Raw p,k are the product price and cost of raw materials in process k respectively; AOT k is the annual operating time for process k.Note that y Raw p,k is obtained from the input-output model as a negative value.This is because it is an input stream into the system, as stated previously in Equation (1).
On the other hand, the environmental impact of each process k is also evaluated.For this work, the WAR algorithm is adapted to evaluate the environmental impact because it can determine the average possible impact of a chemical process on the environment, based on environmental impact values of different pollutants and their respective mass flows.According to the WAR algorithm, the total environmental impact generated by process k, WAR k can be expressed as in Equation (4) as follows.
where i i where y Elec i,k is the power consumption of process k; PEI a and PEI Elec are the score of the potential environmental impact of chemical component a and electricity respectively; α l is the weighting factor of the impact at category l; PEI a,l and PEI Elec l are the potential environmental impact score of chemical component a and electricity at each of the category l respectively.Note that a total of eight categories are considered for PEI, that is, human toxicity potential by exposure, both dermal and inhalation (HTPE), human toxicity potential by ingestion (HTPI), aquatic toxicity potential (ATP), terrestrial toxicity potential (TTP), ozone depletion potential (ODP), global warming potential (GWP), acidification potential (AP) and photochemical oxidation potential (PCOP).
In order to evaluate the health impact of process k, the Hazard Quotient Index (HQI) is selected for this work.This is because HQI can be used for simple process flow diagrams (PFDs) containing limited information, namely process drawings and process descriptions.Such method is suitable for the case of preliminary process design optimization as it allows for the comparison of alternative processes by ranking them based on the risk value based on minimal available information.The calculation of HQI consists of four parts, i.e., estimation of fugitive emissions, air volumetric flow rate, airborne chemical concentration and the health risk (HQI) value.Fugitive emission of each of the chemical components a in process k, m a,k can be calculated via Equation (10): where x j,k , w a,j,k and FE p,j,k are the capacity of the unit operations j, weight composition of chemical component a and the estimated fugitive emission rates in stream p to or from unit operation j respectively.The pre-calculated fugitive emission rates database for the unit operation stream can be obtained from Hassim et al. [15].The air volumetric flow rate is determined by the following Equations ( 11)-( 14).
where A T k and A j,k are the total process floor area and the floor area of each unit operations j in process k respectively; s k is the side length of the process k; A c k is the cross-sectional area of the process; h k is the average height of the main unit operations in process k; Q k is the air volumetric flow rate; v is the wind speed.The average concentration of the chemical components in the air at the downwind edge of the plot area, C a,k can be determined using Equation (15): HQI of component a in process k can be calculated using Equations ( 16) and ( 17): ) where HQI a,k is the HQI of each chemical components a; C EL a is a constant which represents the threshold limit of the chemical component a; HQI k is the total HQI of process k.
The safety impact of a process is assessed via Inherent Safety Index (ISI).ISI is selected because it is a simple scoring method that can be incorporated into an optimization framework.Moreover, ISI considers all parts of the process and equipment, unlike the PIIS method.The total ISI score for The sub-index for Chemical ISI, I k CI is expressed as below: where I RM, max are the sub-index for the factor of the heat of main reaction, heat of side reaction, chemical interaction, flammability, explosiveness, toxic exposure and corrosiveness of the chemical components present in the process k respectively.The sub-index for process ISI, I k PI is given as: where I I k , I T, max k , I P, max k , I EQ, max k and I ST, max k are the sub-index for the factor of inventory, process temperature, process pressure, equipment safety and safe process structure, respectively.The calculations for both I k CI and I k PI are performed on the basis of the worst-case scenario.Among all chemical components present in the process, the scores of the chemical with the most severe hazard in terms of flammability, explosiveness and toxicity are used in the I k CI calculation.Besides, the highest temperatures of the main and side reaction and the worst possible chemical substance interaction in the process are considered.Meanwhile, the maximum expected values for inventory, process temperature, and pressure and the worst process structure are taken into account for the calculation of I k PI .The sub-index for inventory is influenced by the total output of the process, y Out k which is calculated using Equation ( 21): where y Out p,k is the mass flow rate of the output stream for product and effluent streams.The score of the sub-index for inventory is assigned using Equations ( 22) and ( 23) [39].In Equation (22), L k,r and U k,r are the lower and upper bounds of a given criteria r.The criteria r can often be in the form of a range.For instance, r = 1 may refer to the criteria where mass flowrate is between the range of 300 -500 kg/hr.In this context, 300 kg/hr would be L k,r and 300 kg/hr would be U k,r .If y Out p,k falls between this range r = 1 or meets this criteria, Equation (22) will assign a specific score S k,r , which is activated by binary variable b k,r . ( Although several objectives are considered simultaneously in this work, it is important to note that these objectives conflict with each other.To address such conflict, this work employs a multi-objective approach known as fuzzy optimization.Fuzzy optimization is a simple multi-objective optimization approach that was founded upon the fuzzy decision-making approach introduced by Bellman and Zadeh [40].Zimmermann [41] then extended the fuzzy approach to deal with linear and non-linear programming problems that contain multiple objectives.In this approach, a continuous interdependence variable, λ, which is also known as the degree of satisfaction, is introduced.Every fuzzy constraint will be satisfied partially at least to λ.Thus, the multi-objective functions in the optimization can be integrated into a single objective function within the optimization framework.
Fuzzy optimization has been widely used to determine the optimum process alternative that considers all aspects simultaneously based on the pre-defined limits.The fuzzy optimization model is expressed as: where integer variable, b k is used to indicate the existence (or absence) of λ k of process k.Note that the formulation in Equation ( 24) results in the model being non-linear.The optimization objective of λ is maximized, subject to the predefined upper and lower bounds.All flexible targets (GP, HQI, ISI and WAR) are predefined as fuzzy goals which are given by a linear membership function bounded by the upper (GP U , HQI U , ISI U , WAR U ) and lower limits (GP L , HQI L , ISI L , WAR L ). λ of each processes k is determined using Equations ( 26)-( 29): Following this, the model is solved via fuzzy optimization to select a more sustainable process configuration based on the aforementioned four objectives.Based on the optimized solution, the operating capacities, x j,k of the selected process configuration are then examined.If the operating capacities do not differ from the values used in the preliminary simulation, then the selected process k can be recommended for the next process design phase.However, if the operating capacities do differ from the preliminary simulation, the selected process k would be re-simulated in accordance to the optimized operating capacities.Once the selected process k is re-simulated, the operating capacities are checked.In the event where the operating capacities in the simulation exceed an allowable error margin (i.e., ≥10%), the preliminary simulation must be revisited for troubleshooting.For cases where the operating capacities are within the allowed error margin, the re-simulated conditions for process k can be recommended for the next process design phase.Note that the error margin can be changed based on several factors such as type of industry, experience and knowledge of the decision-maker, and type of equipment used.

Case Study
To illustrate the proposed framework, a case study on 1,4-butanediol (C 4 H 10 O 2 ) production process is presented.1,4-butanediol is a colorless and non-corrosive solution which is widely used as a solvent in the industry to manufacture elastic fibers, technical plastics and polyurethanes [42].Based on literature review, Reppe [42] and Davy [43,44] processes are the two common and well-established alternative processes available to produce 1,4-butanediol.Table 1 shows the reaction pathway for both processes.As shown in Table 1, the raw materials needed for Reppe process are formaldehyde (C 2 H 2 ), acetylene (CH 2 O) and hydrogen (H 2 ).Formaldehyde and acetylene are fed into two batch reactors which operate at 5 bar, 80 • C and are arranged in parallel.Formaldehyde is reacted with acetylene in the presence of copper, bismuth and silicon dioxide supporting catalyst to produce 1,4-butynediol as an intermediate product with a conversion rate of 60%.Aqueous sodium hydroxide and hydrochloric acid are used to maintain the pH of the medium in the reactor within the range of pH 6 to 8. Based on the reaction, process synthesis is started by establishing the remaining unit operations in the process.For example, the product from the batch reactors is sent to a buffer tank and is stored for approximately 1-30 h before feeding to a distillation column for further separation.After separation, the bottom product stream of the distillation column is rich in 1,4-butynediol and is fed to a trickling and hydrogenation reactor for further reaction to produce 1,4-butanediol as the main product with a conversion rate of 90%.The operating condition of this reactor is 160 • C and 150 bar. Figure 3 shows the process flow diagram for the Reppe process.
For the Davy process, the raw materials are dimethyl maleate (DMM), C 6 H 8 O 4 and hydrogen, H 2 .DMM is completely vaporised with H 2 before feeding into the first fixed bed reactor.The fixed bed reactor operates at 15 bar, 150 • C and consists of two bed of catalysts which are palladium on alumina and copper-zinc oxide.In the reactor, DMM is fully converted into dimethyl succinate (DMS) through hydrogenation process in the presence of palladium on alumina as catalyst in the first reaction stage.In the second reaction stage, DMS is mainly converted into gammabutyrolactone (GBL) with trace amounts of tetrahydrofuran (THF) and 1,4-butanediol in the presence of copper-zinc oxide as catalyst.The conversion rate of the second reaction is 97%.Based on these two reactions, process synthesis is conducted to list all the other required unit operations.For instance, the vapour-phase product from the reactor is condensed at 30 • C and is fed into a gas-liquid separator to remove the unreacted H 2 .The unreacted H 2 is then being recycled to the fixed bed reactor for further reaction.Then, the GBL rich stream from the first reactor is fed into the second hydrogenation reactor for further reaction to produce 1,4-butanediol with trace amounts of methanol and THF as by-products.The operating condition of the second reactor is 210 • C and 75 bar.The catalyst used in the second hydrogenation reactor is copper-zinc oxide and the conversion of GBL to 1,4-butanediol is around 95%.The PFD for the Davy process is shown in Figure 4.The annual production rates for both Reppe and Davy process are assumed to be 60,000 tonnes of 1,4-butanediol, operating at 8000 h/y.The acceptable range of fluctuation on the amount of the product output given by the decision maker is within ±30%.The average height of the main unit operations in the process is assumed to be below 7 m [45].Wind speed is assumed to be 4 m/s, which is a typical value for outdoor facilities, since local average wind speed is not available [45].
Based on the information mentioned above, process simulation can be performed to analyse the performance for both Reppe and Davy processes.In this work, both processes are simulated with the aid of commercial process simulation software, Aspen HYSYS version 8.8 (Aspentech, Bedford, MA, USA, 2014) [46] (see simulation flowsheets in Figures S1,S2 and simulation files in Supplementary Materials respectively).Settings and parameters used for both simulations are summarized in Tables S1 and S2 (Supplementary Materials) respectively.
Based on the simulation result, data such as input and output mass flow rates and energy consumption are then extracted.The extracted data was then used to develop the input-output models for both processes (see Tables S3 and S4 in Supplementary Materials).
In the following step, the input-output models are evaluated based on economic, environmental, health and safety aspects.For the environmental aspects, the weighting factor of the impact at category l,α l is assumed as one.This indicates that the impact for each category is equally important.Meanwhile, Table 2 shows the information obtained for both Reppe and Davy processes prior to the optimisation step.As shown in Table 2, the Davy process has higher initial GP and WAR; however, with a lower score of HQI and ISI, compared with the Reppe process.The latter has a lower initial GP mainly due to its higher cost of raw materials.Its higher HQI score is mainly due to the presence of harmful chemical components in the process, such as formaldehyde and 1,4-butynediol.Note that the threshold limit (TLV) of formaldehyde and 1,4-butanediol are very low, that is, 1.228 mg/m 3 and 0.5 mg/m 3 respectively.Besides, the higher ISI score for the Reppe process is mainly due to its highly exothermic side reaction involving formaldehyde, acetylene and hydrogen to form 2-propanol.Note that the Reppe process has a lower score of WAR (i.e., it is safer) because it converts the hazardous raw material (formaldehyde) to a product (1,4-butanediol) with a lower environmental impact.As such, the selection between these two processes is complex, particularly when four different aspects are considered simultaneously.In order to determine the optimum process with consideration of all aspects simultaneously, fuzzy optimisation is used.Based on Equations ( 26)- (29), it is noted that lower and upper limits for GP, HQI, ISI and WAR are required.To obtain these values, GP k is first maximized to determine the upper limit GP U .The corresponding values for HQI k , ISI k and WAR k are used as upper limits in HQI U , ISI U and WAR U respectively.Following this, WAR k is minimized to determine WAR L .The corresponding GP k , ISI k and WAR k values were taken as lower limits GP L , HQI L and ISI L respectively.The obtained values for GP U , GP L , HQI U , HQI L , ISI U , ISI L , WAR U and WAR L are 185.112× 10 6 USD/y, 22.817 × 10 6 USD/y, 0.593, 0.017, 36, 27, 50.519 and 1.045 respectively.The model is then solved by maximizing λ in Equation (24) with constraints in Equations ( 1)-( 23) and ( 26)- (29).The model was solved using LINGO version 14.0 (Lindo Systems, Chicago, IL, USA, 2015) on a Lenovo P700 with 8 GB RAM and Intel®Core™ i7, 2.60 GHz Processor.The mixed integer non-linear programming (MINLP) model consists of 187 variables, 6 integers and 341 constraints.
Based on the results obtained, the Davy process was selected as it has a higher value of λ, that is, 0.467, as compared to that of the Reppe process (0.111).The production of 1,4-butanediol after optimisation is scaled down from 7.9 ton/h to 7.64 ton/h, which is within the acceptable range for the amount of product output.Table 3 shows the comparison of the four aspects before and after optimisation for the Davy process.After optimisation, it is observed that the operating capacity of the process has been scaled down by 3.3% from the target production to obtain a trade-off among economic, health, safety and environmental aspects.The reason for this is that, in order to achieve the trade-off solution shown in Table 3, the operating capacity for the Davy process has to be scaled down by 3.3%.From here, the new scaled down operating capacity of the Davy process is re-simulated to study the deviations.The net output of the process streams obtained from the developed model has a slight difference to the solution generated from process simulator with a deviation of less than 10%.This is mainly due to the different models used in generating the solution.The solution obtained by the mathematical model is linear, while that used by the process simulator is in rigorous mode.Since the deviation obtained is small, the selected Davy process with its optimized operating capacity can be recommended for the next stage of design.

Conclusions
In conclusion, this work presented a systematic framework to screen and select a sustainable chemical process at the preliminary design stage.The presented framework consists of three main tasks.First, process simulation was carried out to analyse the performance of each alternative process.Based on the simulation, process data such as mass flow rates and energy requirements for each process unit were determined.These data were then used in the next task, which is to develop an input-output model for each process.Each input-output model was formulated to address four conflicting objectives.These four objectives, namely economic performance, health, safety and environment aspects were assessed using gross profit (GP), the Health Quotient Index (HQI), the Inherent Safety Index (ISI) and the WAR algorithm respectively.Following this, fuzzy optimization was used in the third task to optimize the abovementioned aspects simultaneously and to select the most sustainable process.A sustainable process in this context refers to the process with the highest degree of satisfaction (λ) among all four aspects considered.The proposed framework was illustrated with a case study on 1,4-butanediol production, where the proposed framework was used to determine the most sustainable process to produce 1,4-butanediol.Two processes, namely the Reppe and Davy processes were considered.Each process was simulated using Aspen HYSYS Version 8.8 and subsequently underwent input-output modelling.The developed input-output models were then optimized using fuzzy optimization.Results from the case study indicate that the Davy process is more sustainable compared to the Reppe process.In particular, the Davy process had the highest λ value.Furthermore, the results suggest that the operating capacity of the Davy process should be scaled down by 3.3% to meet the trade-off scores determined by the λ value obtained.Based on this recommendation, the Davy process was re-simulated.It was found that there was a small percentage of deviations between the mathematical and process simulation models.This proved that the developed model is feasible for determining a more sustainable process configuration.For future work, a selection of different unit operations may be incorporated into the model to determine their impact on flowsheet selection.

EP L
Lower limit for gross profit EP U  Upper limit for gross profit WAR L Lower limit for environmental impact WAR U Upper limit for environmental impact ISI L  Lower limit for inherent safety index ISI U  Upper limit for inherent safety index HQI L  Lower limit for health quotient index HQI U  Upper limit for health quotient index PEI a Score of potential environmental impact for chemical component a PEI a,l Score of potential environmental impact for chemical component a at each of the category l PEI Elec Score of potential environmental impact for electricity α l Weighting factor of the impact at category l FE i,j,k Estimated fugitive emission rates in stream i to or from unit operation j in process k Existence of alternative process k

Figure 1 .
Figure 1.The extended Onion Model [5].Reproduced with permission from Foo, D.C.Y.; Ng, D.K.S., Process Integration for Cleaner Process Design; In Handbook of Process Integration: Minimisation of Energy and Water Use, Waste and Emissions; published by Elsevier Science, 2013.

Figure 1 .
Figure 1.The extended Onion Model [5].Reproduced with permission from Foo, D.C.Y.; Ng, D.K.S., Process Integration for Cleaner Process Design; In Handbook of Process Integration: Minimisation of Energy and Water Use, Waste and Emissions; published by Elsevier Science, 2013.

Figure 2 .
Figure 2. Systematic design framework for preliminary process design.Figure 2. Systematic design framework for preliminary process design.

Figure 2 .
Figure 2. Systematic design framework for preliminary process design.Figure 2. Systematic design framework for preliminary process design.
k are the input and output rates of impact in the process k, i (ep),out k is the power consumption for process k.The value of i j,k PEI a ∀k k is represented by ISI k .ISI k is calculated by summing the Chemical ISI, I k CI and Process ISI, I k PI as shown in Equation (18):

Figure 3 .
Figure 3. Process flow diagram of the Reppe Process.

Figure 3 .
Figure 3. Process flow diagram of the Reppe Process.

Figure 4 .
Figure 4. Process flow diagram of the Davy Process.Figure 4. Process flow diagram of the Davy Process.

Figure 4 .
Figure 4. Process flow diagram of the Davy Process.Figure 4. Process flow diagram of the Davy Process.

k
Total process floor area in process k A j,kThe floor area of each unit operations j in process k s kSide length of the process k A c kCross-sectional area of process k h kAverage height of the main unit operations in process kQ k Air volumetric flow rate in process k V Wind speed C a,kAverage concentration of chemical components a in the air at downwind edge of the plot areaC EL a Threshold limit of the chemical component a L k,rLower bound for a given criteria r in process k U k,rUpper bound for a given criteria r in process k S k,r Score of sub-index for inventory in process k based on a given criteria r Variables x j,k Capacity or size of the unit operation j in process k y p,k Net input or output of stream p in process k y Out p,k Net output of stream p in process k y In p,k Net input of stream p in process k y Elec p,k Power consumption of process k w a,j,k Weight composition of chemical component a from unit j in process k m a,k Fugitive emission of each of the chemical components a in process k GP k Gross profit for process k WAR k Total environmental impact generated by process k i (cp),out k Output rates of impact in the process k i (cp),in k Input rates of impact in the process k i (ep),out k Impact of power consumption for process k I k CI Chemical ISI score of process k I k PI Process ISI score of process k I RM, max k Sub-index for the factor of the heat of main reaction in process k I RS, max k Sub-index for the factor of the heat of side reaction in process k I INT, max k Sub-index for the factor of chemical interaction in process k I FL, max k Sub-index for the factor of flammability in process k I EX, max k Sub-index for the factor of explosiveness in process k I TOX, max k Sub-index for the factor of toxic exposure in process k I COR, max k Sub-index for the factor of corrosiveness of the chemical components in process k I I k Sub-index for the factor of inventory in process k I T, max k Sub-index for the factor of process temperature in process k I P, max k Sub-index for the factor of process pressure in process k I EQ, max k Sub-index for the factor of equipment safety in process k I ST, max k Sub-index for the factor of safe process structure in process k ISI k Total ISI score of process k HQI a,k HQI of each chemical components a in process k HQIk Total HQI score of process k λ k Degree of satisfaction for multiple objectives in each alternative process k b k

Table 2 .
Data obtained for Reppe and Davy processes before optimisation.

Table 3 .
Comparison of the four aspects for the Davy process before and after optimisation step.