Modeling simple trigeneration systems for the ... - Semantic Scholar

3 downloads 62113 Views 588KB Size Report
Dec 9, 2011 - Umberto software is specifically designed to analyze the distribution of material and energy ... bookkeeping and cost accounting, allowing the setup of complex .... to evaluate the environmental loads associated with a product,.
Environmental Modelling & Software 30 (2012) 71e80

Contents lists available at SciVerse ScienceDirect

Environmental Modelling & Software journal homepage: www.elsevier.com/locate/envsoft

Modeling simple trigeneration systems for the distribution of environmental loads Monica Carvalho a,1, *, Miguel A. Lozano a, Luis M. Serra a, Volker Wohlgemuth b a

Group of Thermal Engineering and Energy Systems (GITSE), Aragon Institute of Engineering Research (I3A), Department of Mechanical Engineering, University of Zaragoza, CPS de Ingenieros (Edf. Agustin de Betancourt), C/Maris de Luna, s/n, 50018, Zaragoza, Spain b HTW Berlin, University of Applied Sciences, Industrial Environmental Informatics Unit, Wilhelminenhofstr. 75A, 12459 Berlin, Germany

a r t i c l e i n f o

a b s t r a c t

Article history: Received 4 May 2011 Received in revised form 1 November 2011 Accepted 12 November 2011 Available online 9 December 2011

Integration of thermoeconomics and Life Cycle Analysis was carried out within the framework of an Environmental Management Information System. This combined approach identified where environmental loads were generated and tracked environmental loads throughout the system, allowing for a more precise understanding of operational activities. A trigeneration system was modeled, providing electricity, heat, and cooling to a building. The trigeneration system consists of a cogeneration module, auxiliary boiler, absorption chiller and electrical chiller. The trigeneration system model is flexible, as it allows electricity from/to the electric grid to be purchased/sold, and part of the cogenerated heat to be wasted. Umberto software is specifically designed to analyze the distribution of material and energy resources throughout a productive system. The software is based on Petri networks, double-entry bookkeeping and cost accounting, allowing the setup of complex systems and also a combined material, energy and inventory calculation. An assistant was built to include the tracking of emissions through the application of algebra and rules similar to those used in thermoeconomic analysis. It is possible to evaluate the environmental impact in terms of the consumption of natural resources and generation of emissions in the system, from the input of natural resources to the output of the final products. Network parameters were used to calculate the emissions associated with the operation of the system. The issue of allocating environmental loads was introduced and two scenarios for each operational mode were compared: the trigeneration system vs. a conventional energy supply system in which electricity was produced in a representative coal power plant. In this case the trigeneration system operated with significant reduction of the CO2 emitted into the atmosphere. Ó 2011 Elsevier Ltd. All rights reserved.

Keywords: Trigeneration CO2 Material flow networks Environmental costs Umberto software Environmental Management Information System

Software availability Umberto software: Developed by ifu Hamburg GmbH. Contact address: ifu Hamburg GmbH, Max-Brauer-Allee 50, D-22765 Hamburg, Germany. Tel.: þ49 404800090; fax: þ49 404800922. [email protected], http://www.umberto.de/en/ LINGO software: Developed by LINDO systems Inc. Contact address: 1415 North Dayton St, Chicago, IL, 60622, United States. Tel.: þ1 312 988 7422. [email protected], www.lindo.com

* Corresponding author. Present address: MIRARCO, Laurentian University, 935 Ramsey Lake Road, Sudbury, ON P3E 2C6, Canada. Tel.: þ1 705 675 1151x6613; fax: þ1 705 675 4838. E-mail addresses: [email protected] (M. Carvalho), [email protected] (M.A. Lozano), [email protected] (L.M. Serra), [email protected] (V. Wohlgemuth). 1 Tel.: þ34 976 76 1000x5258; fax: þ34 976 76 2616. 1364-8152/$ e see front matter Ó 2011 Elsevier Ltd. All rights reserved. doi:10.1016/j.envsoft.2011.11.005

SimaPro software: Developed by PRé Consultants. Contact address: Printerweg 18, 3821 AD Amersfoort, The Netherlands. Tel.: þ31 33 450 4010. [email protected], http:// www.pre-sustainability.com/

1. Introduction Residential and tertiary sectors are responsible for more than 40% of final energy consumption in the European Community (Directive COM 2002/91/EC, 2002). The tertiary sector includes different types of buildings (hospitals, schools, hotels, etc.) with a great variety of uses and energy services (heating, cooling, and electricity). European research projects, such as CHOSE (2001), TRIGEMED (2003), and Summerheat (2008), concluded that there is a significant technical potential for the implementation of trigeneration in the residential and tertiary sector of countries in the Mediterranean area. In these

72

M. Carvalho et al. / Environmental Modelling & Software 30 (2012) 71e80

Nomenclature AB AC C CM COP Ed Ep Er Es EM ExC EC Fc Fa L MFN P pep

auxiliary boiler absorption chiller operation mode (followed by subscript) cogeneration module coefficient of performance electricity demand electricity purchased from the grid electricity input to electrical chiller electricity sold to the grid (cogenerated) environmental loads (followed by subscript) example of operation mode (followed by subscript) electrical chiller fuel for the cogeneration module (natural gas) fuel for the auxiliary boiler (fuel oil) loss node Material flow networks Purchase node price of purchased electricity

countries the need for heating is restricted to few winter months, limiting the application of cogeneration systems. There is, however, a significant need for cooling during the summer period. One solution is the use of absorption chillers for cooling. By combining cogeneration with absorption chillers (trigeneration), the energy demand can be extended into the summer months to match cooling loads. Advantages of trigeneration systems in buildings have been demonstrated in literature, as the improved use of fuel is associated with economic savings and sparing of the environment, as less fuel is consumed and consequently less pollution is generated (Maglorie et al., 2002; Chicco and Mancarella, 2008). As sustainability-related issues such as energy consumption and environmental impact become a more integrated part of operational and long-term planning decisions in energy systems, simulation modeling and analysis tools are needed to aid in the decision making process. Thermoeconomics is a potent tool for energy analysis and has been used to support the design, synthesis and operation of energy systems by providing crucial information not available through conventional analyses. Thermoeconomics combines economic and thermodynamic analysis with the purpose of revealing opportunities of energy and cost savings when designing and operating energy conversion systems (El-Sayed and Evans, 1970; Gaggioli, 1983, ElSayed and Gaggioli, 1989; Lozano and Valero, 1993; Serra et al., 2009). The basic concept of thermoeconomic analysis is the energy cost, understood as the amount of energy resources consumed for obtaining a piece of equipment, a flow or a commodity. Hence, the cost of a flow in a plant represents the amount of resources that have to be supplied to the overall system to produce this flow. Unit costs are used by thermoeconomic cost accounting theories for rational price assessment, and allow us to follow the cost formation process throughout the system, from energy resources to final products. Thermoeconomic costs can be expressed in monetary or energy units. The necessity of considering the environment as an additional design factor is an increasing demand due to the uprise of environmental conscience and the requirements to reduce the environmental impact of modern society. Life Cycle Analysis (LCA) provides a more global perspective of environmental loads and has the potential to fulfill the need for an adequate design tool for energy supply systems (Guinée, 2002). LCA is an objective process to evaluate the environmental loads associated with a product, process, or activity (Awasthi and Chauhan, 2011; Carvalho et al.,

pes pfa pfc Q Qa Qc Qcc Ql Qr Qd R Rd Re Rq rql S Wc Wcc

a h

price of electricity sold to the grid price of fuel oil (for auxiliary boiler) price of natural gas (for cogeneration module) heat distribution node heat produced by auxiliary boiler heat produced by cogeneration module heat internally consumed wasted heat heat input to absorption chiller heating demand cooling node cooling demand cooling produced by electrical chiller cooling produced by absorption chiller cost of waste heat sale node electricity produced by cogeneration module electricity internally consumed efficiency (cogeneration module) efficiency (boiler)

2011a,b; Turner et al., 2011). LCA also identifies and quantifies the use of mass and energy as well as the emissions to the environment, determining the impact of the use of resources and emissions, allowing for evaluation and implementation of strategies of environmental improvement. The life cycle or cradle-to-grave impacts include those resulting from extraction of raw materials, fabrication of the product, transportation or distribution of the product to the consumer, use of the product by the consumer, and disposal or recovery of the product after its useful life. Herein, the environmental cost of a flow represents the amount of environmental loads that have been generated in the overall system to produce this flow. Thermoeconomic and LCA techniques are both based on the premise that all of the resources required for producing a good or service need to be accounted for. Thermoeconomics is usually applied to energy conversion systems and the limits of the system are those of the associated plant. However, there is no constraint that impedes widening the limits of analysis to include the well or the mine from where the natural resources were extracted. Thus, merging thermoeconomics and LCA methodologies provides a global perspective of a complex system via an integrated analysis of energy, economics and environment. Generally, the analyzed system in LCA is treated as a black box from which only its inputs and outputs are measurable, without further knowledge of the inner structure. Applying the philosophy of thermoeconomics to energy systems opens this black box and unravels the process of environmental burden formation, which is where the importance of combining thermoeconomics with LCA lies (Carvalho, 2011). Similar to the cost formation process in thermoeconomics, it will be possible to evaluate the process of formation of the environmental impact linked with consumption of natural resources and distribution of environmental loads throughout the system e i.e., from the input of natural resources to the output of final products. Integration of thermoeconomics and LCA was carried out within the framework of an Environmental Management Information System (EMIS). EMIS are designed to detect, evaluate and prevent a wide range of environmental dangers and stresses. In more concrete terms, EMIS consist of computer programs that support management by collecting, documenting and evaluating all relevant data about an enterprise’s interaction with its environment and plan, initiate and control all activities related to environmental

M. Carvalho et al. / Environmental Modelling & Software 30 (2012) 71e80

protection (Pokorný, 2006; Wohlgemuth et al., 2006; Díez and McIntosh, 2009.). Combination of different techniques or analyses is not uncommon, addressing the shortcomings of individual analyses and allowing for a more complete result, such as the combination of Geographical Information System analysis, System Dynamics model and 3D visualization for sustainability assessments (Xu and Coors, 2011) or the combination of PrincipaleAgent models with Linear Programming in order to deal with the design of environmental regulation (Viaggi et al., 2009). Umberto software (2006) is an EMIS and has been specifically designed for analyzing the distribution of material and energy resources throughout a productive system. Umberto can be used as an instrument for analyzing complex networks of materials and energy flows, stocks, as well as material and energy transformations. This study modeled a simple trigeneration system with the Umberto software, incorporating environmental information on the consumption of resources. Even if a network can provide data on the level of material and energy flows, it is necessary to extend the boundaries because life cycle assessments require a representation of the entire life cycle of products and services, including raw material extraction, distribution, use phase, and waste disposal (Möller, 2010). Such an extension is interesting for companies, as the concept of resource productivity is a new administrative approach to deal with sustainability challenges (Porter and van der Linde, 1995). The combined methodology allowed consumers of electricity, heat, and cooling produced by trigeneration systems to know the unit environmental loads (equivalent to the thermoeconomic unit costs) associated with the consumption of each energy service. Incorporating economic and environmental criteria into the design and planning process should strive toward the following: (1) increased efficiency in the use of energy and materials; (2) reduction of energy cost of final products; and (3) reduction of environmental burden. The allocation of environmental loads in multiproduct processes is a much debated and critical issue. This is the case for trigeneration plants, and in general, for any energy plant producing more than one useful energy flow. Through a detailed examination of the operation modes of a trigeneration system, an allocation proposal is presented, providing better insight on the characteristics of a trigeneration system (exposing the distribution of environmental loads throughout the trigeneration system). The environmental allocation proposal provided energy services with fewer environmental loads than those associated with separate production,2 which is very important on a consumer-oriented basis. If the consumers of the trigeneration system assess that allocation was fair, their buy-in is more likely to occur. A fair apportionment will contribute to the acceptance of the more complex but more efficient trigeneration systems by consumers, which is essential for the success of such systems when they are oriented to multiple users.

The technology behind trigeneration is fundamentally based on the coupling of a cogeneration module with an absorption chiller. The simple trigeneration system analyzed herein is complemented by the usual equipment present in a conventional plant: a hot water boiler and an electrical chiller (Lozano et al., 2009a). Both technologies are used to guarantee supply and also to avoid oversizing the cogeneration module and the associated absorption chiller. The idea is that the cogeneration module, jointly with the absorption chiller, satisfies the average thermal demand for heat and cooling, while the conventional units (boiler and electrical chiller) are utilized in an auxiliary way to make up for the demand peaks. Therefore, supply is guaranteed and the installation is reliable, since the existence of conventional equipment assures the satisfaction of the thermal demand (Lozano et al., 2009b). Trigeneration systems have wide range of applications: single residential applications (Wang et al., 2008), industry (Colonna and Gabrielli, 2003), or buildings (Marimón et al., 2011). 2.1. Structure and operation The purpose of the trigeneration system (Fig. 1) was to meet the demand of different energy services (electricity, Ed; heating, Qd; and cooling, Rd) of a consumer center. Fig. 1 shows a diagram of the analyzed trigeneration system, with internal and product flows. The simple trigeneration system consists of the following productive units: a cogeneration module CM (providing heat, Qc, and electricity, Wc), an auxiliary boiler AB (providing heat, Qa), an absorption chiller AC (providing cooling, Rq, and driven by heat, Qr) and an electrical chiller EC (providing cooling, Re, and driven by electricity, Er). Table 1 presents the technical parameters of the simple trigeneration system. Efficiency coefficients were obtained from equipment catalogs and consultations with manufactures. For each equipment, Pnom is the nominal capacity (power of its main product). Taking the cogeneration module as an example, electricity is the main product and Pnom ¼ Wc nom ¼ 350 kW. To produce Wc kW of electricity, Fc ¼ (1/aw)Wc kW of natural gas will be consumed, recovering Qc ¼ (aq/aw)Wc kW of heat. It was considered that the efficiency coefficients were constant and independent from the production P  Pnom of the equipment. Demands were always met either by the productive units of the trigeneration system or with the help of purchased electricity from the electric grid (Ep). The possibilities also existed that a fraction (Ql > 0) of the cogenerated heat could be wasted, and the electricity

2. Simple trigeneration system Trigeneration is the combined production of heat, cooling and electricity from the same source of energy. The benefits of trigeneration arise from the comprehensive integration of the processes and technologies used to supply energy in all necessary forms to meet defined energy demand profiles. Trigeneration technologies have socio-economic and environmental benefits that relate to their efficient use of energy resources, typically resulting in the reduction of operation costs and in environmental benefits (reduced carbon emissions).

2 Electricity purchased from the grid, heat produced by a boiler, and cooling produced by an electric chiller driven by purchased electricity.

73

Fig. 1. Simple trigeneration system with flows.

74

M. Carvalho et al. / Environmental Modelling & Software 30 (2012) 71e80

Table 1 Technical parameters of the simple trigeneration system. Transition/equipment

Efficiency coefficient

Nominal capacity (kW)

Cogeneration module

aw h Wc/Fc ¼ 0.35 aq h Qc/Fc ¼ 0.40 hq h Qa/Fa ¼ 0.80

Wc

Auxiliary boiler Absorption chiller Electrical chiller

COPq h Rq/Qr ¼ 0.625 COPe h Re/Er ¼ 5.0

nom ¼ 350

Qa nom ¼ 400 Rq nom ¼ 250 Re nom ¼ 250

could be sold to the market (Es). Fc and Fa refer to the fuel utilized by the cogeneration module and the auxiliary boiler, respectively. Environmental loads were generated in the cogeneration module, which operates on natural gas and in the auxiliary boiler, which operate on fuel oil. When electricity was purchased or sold from/to the electric grid, environmental loads were also considered. The minimal operation CO2 emissions of the system were evaluated by a linear programming model, as explained in section 2.3. The model was applied to very different demands of energy services and the results show how different the operation modes for a simple trigeneration system can be. A greater sophistication of the optimization model, using non linear production restrictions taking into account how ambient conditions affect capacity and efficiency coefficients and binary variables limiting both the minimum load of the equipments and the on/off status, would provide more precise results but, generally, the achieved conclusions would still prevail. The approach presented in this paper can be used for whole system design by considering, in the optimization model and analysis methodology, the life cycle environmental loads of the system: a) the environmental loads embedded in the selected equipment for the plant and b) the environmental loads associated with the operation of the plant. In this paper, only the second part (b) was studied: the minimization of the environmental loads associated with the operation of a previously dimensioned trigeneration system. Minimization of environmental loads associated with the whole system should also be considered as an obligatory prior step of project and design methodologies for new trigeneration systems (Carvalho et al., 2011a,b). 2.2. Environmental loads of fuels and electricity SimaPro (2008) is a specialized LCA tool, one of the most widely used LCA softwares, used by major industries and consultants, through to research institutes and universities, and was utilized to calculate the impact associated with the operation of the system (consumption of resources). This was possible because SimaPro includes several inventory databases with thousands of processes and several well-known impact assessment methods. The system interacted with the economic environment (market) through the purchase of natural gas, fuel oil, and electricity from the grid, as well as through the sale of cogenerated electricity to the grid. LCA analyzes the environmental impacts associated with a process or product from ‘the cradle to the grave’, which begins with the gathering of raw materials from the earth to create the product/service and ends at the point when all materials are returned to the earth (SAIC, 2006). 2.2.1. Natural gas Natural gas was characterized by utilizing the related emissions of combustion of natural gas, from the IDEMAT database (IDEMAT, 2001), and the total aggregated system inventory for a natural gas consumer in Spain, from the Ecoinvent database (Ecoinvent, 2007). The CO2 emissions associated with the consumption of natural gas in Spain were obtained by utilizing SimaPro, calculated as EMfc ¼ 0.272 kg CO2 per kWh of consumed natural gas.

2.2.2. Fuel oil Fuel oil was characterized by an inventory module (extraction, production at refinery and transportation from refinery to an average European end user) and related emissions of controlled burning, from the Ecoinvent database. The CO2 emissions associated with the consumption of fuel oil were obtained by utilizing SimaPro, calculated as EMfa ¼ 0.305 kg CO2 per kWh of consumed fuel oil. 2.2.3. Electricity The CO2 emissions associated with the local electricity were also calculated by SimaPro, utilizing the Ecoinvent database and considering that all electricity originated from a single-fuel representative coal power plant (EMep ¼ 1.020 kg CO2/kWh). The evaluation of CO2 emissions in a representative coal power plant was utilized as an example here, in order to illustrate the application of the proposed methodology. Carvalho (2011) analyzed other environmental indicators as well as other types of power plants, and different national electricity mixes. Using coal power plants, different operation modes of the system appear when varying the energy services demands. This allowed for a clearer and complete validation of the proposed methodology. 2.3. Optimal operation modes A linear programming model was solved in order to obtain the optimal operation mode from an environmental viewpoint. The environmental analysis considered that the only significant variable environmental loads were electricity, natural gas and fuel oil. No environmental burden was associated with the waste of cogenerated heat, i.e., EMql ¼ 0. The objective function to be minimized was the operation variable emissions (in kg CO2/h):

Operation variable emissions ¼ EMfc Fc þ EMfa Fa þ EMep Ep  EMes Es þ EMql Ql

(1)

Cogenerated electricity sold to the grid was considered to have the same environmental load as that of electricity purchased from the grid (EMes ¼ EMep). The concept of avoided emissions is presented as the emissions avoided elsewhere with the production of electricity by the cogeneration module, consequently avoiding the purchase of electricity from the grid. Equation (1) was subject to restrictions of capacity limit and equipment efficiency as well as balance equations, previously presented in Lozano et al. (2009a). Results were obtained by utilizing the modeling language and optimizer Lingo (2008). Lingo is a commercial software package for solving optimization problems that uses the branch and bound solver to enforce any integer restrictions contained in a model. The advanced capabilities of Lingo such as cut generation, tree reordering, advanced heuristic and presolve strategies were used as needed. Given the energy demands to be satisfied, according to the different operation modes, Lingo solved the previous model and determined the feasible operation mode with the minimum operation variable emissions. The resulting feasible operation states could be classified into nine different operation modes, based on the values of purchased electricity (Ep), sold electricity (Es), auxiliary heat (Qa) and waste heat (Ql). These operation modes corresponded to different demand of the energy services of the consumer center and are shown in Table 2. A summary of results (demand, flows, and hourly cost) obtained with Lingo for four examples ExC1, ExC3, ExC7 and ExC9 that corresponded to different operation modes (C1, C3, C7 and C9) is presented in Table 3. For each different example, the minimum

M. Carvalho et al. / Environmental Modelling & Software 30 (2012) 71e80

emissions of satisfying the energy service demand of the consumer center was reached in a different operation mode, which exchanged energy flows with the environment and utilized the productive capacity of the installed equipment. Table 3 shows the energy flows obtained and the operation variable CO2 emissions. 3. Environmental management information system 3.1. Umberto software for material and energy flow analysis Software support for Material Flow Networks (MFN) must meet a number of requirements. For instance, the software must be capable of modeling complex production processes from different fields of application, such as chemistry or engineering, and must offer flexible data management for updating and extending the model. Different ways of data interpretation and display must also be supported. The Umberto software, which entered the market in the mid-1990s, was the first suite of programs that tried to meet these goals. Umberto software is a powerful and user-friendly material flow analysis tool. It runs under the Microsoft Windows family of operating systems and meets software ergonomic standards (Wohlgemuth et al., 2006). The inclusion of environmental information on the usage and consumption of resources into Umberto allowed for the evaluation of the distribution of the environmental loads associated with each flow of the system. Additionally, integrated models of energy flows facilitated a better understanding of the assignment of environmental and economic costs to the internal and final products of the trigeneration system. The modeling of material flows in multi-stage production systems initially focused on absolute flows of companies and supply chains. The objective is to trace the material flows within a company or between different companies within a value chain. In MFN, the term material refers to substances and energy, meaning there is virtually no distinction between substances and energy. Based on the concept of material flow networks, the powerful calculation algorithm of Umberto allows for the determination of all material and energy flows in the system under study. According to Wohlgemuth et al. (2006), the most attractive feature of MFN is the possibility to combine the compilation of ecobalances for an industrial plant with an analysis of material flows associated with given products or services. An eco-balance refers to the consumption of energy and resources and the pollution caused by the production cycle of a product. An advantage of the MFN approach resides in its gradual modeling approach, starting from a very basic model of few processes with simple specifications, the model can be extended step by step to include further processes, sites, more complex specifications, costs, etc. (Viere et al., 2010). Umberto software allows for the visualization of processes, units and flows, carrying out mass and energy balances and analyzing from an environmental point of view the emissions generated. Petri Nets, double-entry bookkeeping, and cost accounting are the basis of Umberto software, allowing the setup of complex systems and also

75

Table 3 Energy flows and variable CO2 emissions. ExC1

ExC3

ExC7

ExC9

Ed Qd Rd

kW kW kW

400 400 400

400 100 100

200 600 100

200 100 100

Ep Es Fc Fa Wc Qc Wcc Er Ql Qcc Qa Qr Rq Re

kW kW kW kW kW kW kW kW kW kW kW kW kW kW

100 0 1000 300 350 400 350 50 0 400 240 240 150 250

50 0 1000 0 350 400 350 0 140 260 0 160 100 0

0 130 1000 250 350 400 220 20 0 400 200 0 0 100

0 150 1000 0 350 400 200 0 140 260 0 160 100 0

CO2 emissions Operation mode

kg CO2/h

465.50 C1

323.00 C3

215.65 C7

119.00 C9

a combined material, energy and inventory calculation. Umberto models consist of places, transitions, and arrows (directed graphs). 3.1.1. Places An important function of places is that they delimit the system from its environment; they are points of contact with the world. The inputs of the simple trigeneration system (Fig. 4, light gray circles) were the consumption of fuel by the cogeneration module (Fc) and auxiliary boiler (Fa), and the electricity purchased from the grid (Ep). The outputs of the system (dark gray circles in Fig. 4) were the demands of electricity (Ed), heat (Qd), and cooling (Rd). Freedom was available to the consumer to decide how the system operated: wasting heat permitted the operation of the cogeneration module to match the demand of the consumer center and the sale of surplus cogenerated electricity permitted to realize profit. Therefore two more outputs of the system were waste heat (Ql) and the cogenerated electricity sold to the grid (Es). The place “Environmental loads” accounted for the environmental loads originating from the consumption of natural gas and fuel oil, and from the purchase or sale of electricity. The two “Environmental loads” outputs seen in Fig. 2 are duplicate places. If an arrow leads to a place far away, the graphical display might become incomprehensible. Therefore the “Environmental loads” place was duplicated and the copy was positioned in the vicinity of transition P. All emissions go into the atmosphere, but Umberto software tracks the contribution of each transition to account for its share of emissions. 3.1.2. Transitions Each piece of equipment was modeled as a transition. A slightly more complex but more flexible method to specify transitions was

Table 2 Operation modes.

Yes auxiliary boiler No waste heat No auxiliary boiler No waste heat No auxiliary boiler Yes waste heat

Yes purchased electricity

No purchased electricity

No sold electricity

No sold electricity

No purchased electricity Yes sold electricity

C1

C4

C7

C2

C5

C8

C3

C6

C9

76

M. Carvalho et al. / Environmental Modelling & Software 30 (2012) 71e80

Fig. 2. Umberto model of simple trigeneration system.

applied, utilizing expressions to describe the relationships between input and output flows of a transition, using the technical parameters shown in Table 1. To guarantee that the network could be calculated in both directions, the user-defined functions also considered the inverse form. Environmental load values were incorporated into Umberto as net parameters, which are parameters valid not only locally for a single transition but for an entire network or for subnets. Net parameters can be used in the functions of the transition specifications in the same way as transition parameters. Branching and merging points S (Sale node), L (Waste heat node), P (Purchase node), Q (Heat distribution node), and R (Cooling node) were also modeled as transitions. Branching and merging points can be interpreted as decision points, in which possibilities are reflected. Point S refers to the possibility of selling cogenerated electricity to the grid; point L refers to the possibility of wasting part of the cogenerated heat; point P refers to the possibility of purchasing electricity from the grid; point Q refers to the possibility of operating the auxiliary boiler, and point R adds the contributions of the chillers to satisfy the cooling demand.

Umberto offers several interfaces to other programs, one of which allows the specification of transitions. In combination with complex algorithms it can embed sub-models into MFN that have been developed outside of Umberto. This can be done through a script written in any of the languages supporting Microsoft’s Active Scripting Architecture. The new functions and extensions were implemented within the menu structure, utilizing structural language XML with code/logic J#. XML was chosen because it is human-legible and reasonably clear, easy to create, and the user has the advantage of acting independently of software (data can be moved through software upgrades and even to different software products). J# was used basically for its text highlighting abilities and because it uses Microsoft’s Common Language Runtime. The assistant was an application that collected data of the calculated flows to carry out cost accounting. The assistant was necessary because Umberto calculates flows and costs simultaneously, and the implementation of thermoeconomic equations required the flows to be previously calculated. The assistant was validated with economic costs (correctly reproducing thermoeconomic cost results published in Lozano et al. (2009a,c). By changing ‘market prices’ to ‘environmental loads’, the assistant turned to an environmental perspective, giving the assistant flexibility to support calculations regarding environmental loads or economic costs. Balances were formulated in each transition and external resources used in the production process were valued by the environmental burden caused. Applying the condition of cost conservation to each transition yields: Equipment

CM :

EMfc Fc ¼ EMwc Wc þ EMqc Qc

(2)

AB :

EMfa Fa ¼ EMqa Qa

(3)

AC :

EMqr Qr ¼ EMrq Rq

(4)

EC :

EMer Er ¼ EMre Re

(5)

Branching and merging points

3.2. Environmental cost accounting

S:

EMwc Wc ¼ EMwcc Wcc þ EMes Es

(6)

As previously mentioned, thermoeconomic analysis (cost accounting) can be combined with LCA, and both techniques can be integrated into an EMIS. In the particular case of Umberto, it was necessary to build an assistant to integrate the philosophy/ methodology utilized in energy cost analysis (thermoeconomics) with the evaluation of environmental loads. The concept of cost can involve different magnitudes, as for example, environmental loads. Environmental costs can be understood as a category of cost (according to the consumption of natural resources and generation of environmental loads in order to obtain a flow). For the implementation of the environmental allocation method based on thermoeconomics, an assistant was created in Umberto software. The assistant performed calculations of environmental loads of internal flows and products after network calculation. The input monitor of Umberto software was used to model the system with greater flexibility, allowing modifications in operation modes to be carried out easily. The flows that defined each operation mode were established in the input model, and then Umberto software calculated the remaining flows.

P:

EMwcc Wcc þ EMep Ep ¼ EMer Er þ EMed Ed

(7)

L:

EMqc Qc þ EMql Ql ¼ EMqcc Qcc

(8)

R:

EMrq Rq þ EMre Re ¼ EMrd Rd

(9)

Q:

EMqcc Qcc þ EMqa Qa ¼ EMqr Qr þ EMqd Qd

(10)

Considering that the operation state of the plant was known, then all energy flows, environmental loads for fuel and electricity and the environmental load entailing waste heat were also known. Here it was considered that EMql ¼ 0 because the objective was to assess all environmental loads to useful final products. EMql ¼ 0 does not mean an absence of environmental impact, but that there was no environmental burden associated with the act of wasting heat. Consequently, there were 12 unit environmental loads of internal flows and final products to be calculated: EMwc, EMwcc, EMer, EMed, EMqc, EMqcc, EMqa, EMqr, EMqd, EMrq, EMre, and EMrd. As

M. Carvalho et al. / Environmental Modelling & Software 30 (2012) 71e80

Fig. 3. Control volume accounting for the interaction of the cogeneration module with the environment.

the system was described using 9 equations with 12 unknowns, 3 auxiliary equations were needed. It was considered that the unit environmental load of several flows obtained from a homogeneous flow was the same. Applying this rule to branching points P and Q, two auxiliary equations were obtained:

P:

EMer ¼ EMed

(11)

Q:

EMqr ¼ EMqd

(12)

This means that the environmental loads assessed to the electricity internally consumed in the electrical chiller (Er) and the electricity provided to the consumer center (Ed) are the same. The environmental loads assessed to the heat internally consumed in the absorption chiller (Qr) and the heat provided to the consumer center (Qd) follow the same rule. The third auxiliary equation must define how the environmental loads generated in the cogeneration module are attributed to its products: heat and electricity. Allocation is a very important issue when apportioning environmental loads to multiproduct systems, to ensure each party is credited with their appropriate share. Research on allocation of emissions and environmental burdens will allow the environmental benefits of co- and tri-generation technologies (adequately designed and operated) to be better understood and exploited (Rosen, 2008; Abusoglu and Kanoglu, 2009; Carvalho, 2011). Umberto software supports the consideration of different approaches to the allocation issue, and even supports the use of scripts to attach complex rules or models to a transition. However, the issue of allocation was not addressed in depth in this paper. A deep discussion on the allocation of environmental loads is presented in Carvalho (2011). Different allocation methods of environmental loads to electricity and heat products (third auxiliary equation for the analyzed system) were found in literature (Phylipsen et al., 1998).

77

However, the main issue found during the utilization of such simple methods focused on the immediate products of the cogeneration module, Qc and Wc, not accounting for possible different destinations or uses of Qc and Wc nor for the interaction of the cogeneration module with its productive environment. Since traditional solutions to the problem were unsatisfactory, the authors proposed an innovative allocation method in Carvalho et al. (2010), which was implemented in Umberto software through the aforementioned assistant, allocating emissions to the consumed products of the cogeneration module (Wcc and Qcc), in proportion to the emissions generated by their alternative production:

EMqcc EMqa ¼ EMwcc EMep

(13)

EMep being the environmental loads corresponding to electricity purchased from the electric grid, and EMqa the environmental loads associated with the heat produced in the auxiliary boiler (EMqa ¼ EMfa/hq). Equation (13) could be applied directly to all operation modes, as it was previously established that EMep ¼ EMes. The choice of such control volume allows the benefits of selling electricity and the inefficiency of wasting heat to be both distributed between heat and electricity internally consumed. When considering different equipment, activities, and options included in the trigeneration system, the assignment of unit costs should rather consider the products of the cogeneration module that were consumed (Wcc and Qcc). In this way, as shown in Fig. 3, the emissions associated with the operation of the cogeneration module, the sale of electricity to the grid, and waste heat (EMfcFc  EMesEs þ EMqlQl) were distributed between heat and electricity internally consumed (EMwccWcc and EMqccQcc). 4. Results and discussion Based on the flow quantities entered in the input monitor (energy demands and environmental loads associated with operation), the system calculated the flows of the entire network using the transition specifications. For comparison, a reference system (Fig. 4) was created, in which all energy demands were satisfied in a conventional way (electricity was purchased from the grid for Ed and Rd, and Qd was satisfied by an auxiliary boiler). Table 4 shows the emissions for the trigeneration and conventional systems. Umberto software provided the overall amount of emissions associated with the operation of the trigeneration system. It was verified that trigeneration technology reduced significantly the kg of CO2 emitted into the atmosphere. Table 5 shows the unit environmental loads associated with the operation of the simple trigeneration system for all operation examples.

Fig. 4. Conventional system.

78

M. Carvalho et al. / Environmental Modelling & Software 30 (2012) 71e80

Table 4 CO2 emissions of trigeneration system vs. conventional system.

Conventional Trigeneration

kg CO2/h kg CO2/h

Table 6 Energy resources and CO2 emissions assessed to final products.

ExC1

ExC3

ExC7

ExC9

642.12 465.50

466.53 323.00

453.18 215.65

262.53 119.00

Comparison of emissions between examples ExC1 and ExC3 gave indications on what occurred when some heat was wasted. EMwc and EMwcc increased their value in ExC3, reflecting the inefficiency of wasting heat, and consequently, EMed had a higher value. EMqd presented a lower value in ExC3, promoting the consumption of waste heat. The comparison between examples ExC1 and ExC7 gave indications on the behavior of the system when electricity was sold to the grid. In ExC7, EMwc presented a higher value and EMwcc presented a lower value than in ExC1. The sale of electricity with lower emissions (but evaluated as having higher emissions) consequently lowered the cost of EMwcc. The benefits of the sale of electricity were positively reflected on the values of EMqc and EMqcc, and ultimately on the final emissions of Ed, Qd and Rd, which were lower. In ExC9, EMed increased reflecting the waste of heat, but with sale of electricity, EMed was still environmentally sounder than EMep. The sale of electricity benefited all final energy services, resulting in lower emission values when comparing ExC3 and ExC9. Internal flows too, were lower in ExC9. When comparing ExC7 and ExC9, it could be seen that EMqc had a lower value, resulting in a lower value for EMqd which should promote consumption of otherwise wasted heat. EMwc and EMwcc had increased values which were translated into higher emissions for EMed and EMrd. The benefits as well as the penalties of the system were reflected in all energy services produced in the cogeneration module. Table 6 shows the breakdown of energy resources and CO2 emissions assessed to final products (electricity, heat and cooling demands). Umberto provided the emissions associated with the consumption of each energy service, and the origin of the emissions, tracking the contributions of electricity purchased from the grid, and fuels of the cogeneration module and auxiliary boiler. Umberto software uses diagrams to display materials, energy and cost flows. Figs. 5e7 show the breakdown of the emissions associated with the consumption of energy services in case ExC1. A scaling variant (Maximum flow) was utilized for visualization purposes, so the maximum width in the network was limited to each energy service. All other widths were calculated proportionally.

Ep (kWh)

Fc (kWh)

Fa (kWh)

CO2 (kg CO2)

ExC1

Ed Qd Rd Total

88.89 e 11.11 100.00

622.72 187.15 190.13 1000.00

e 187.50 112.50 300.00

260.05 108.09 97.36 465.50

ExC3

Ed Qd Rd Total

50.00 e e 50.00

700.56 115.17 184.27 1000.00

e e e e

241.55 31.33 50.12 323.00

ExC7

Ed Qd Rd Es Total

e e e e e

400.32 299.44 40.03 260.21 1000.00

e 250.00 e e 250.00

108.89 157.70 10.89 61.83 215.65

ExC9

Ed Qd Rd Es Total

e e e e e

400.32 115.17 184.27 300.24 1000.00

e e e e e

108.89 31.32 50.12 71.33 119.00

Fig. 5 shows the contributors to the environmental loads of the electricity demand (Ed), which is satisfied by the generation of electricity in the cogeneration module (operating with fuel Fc), and is complemented by the purchase of electricity from the grid (Ep). Fig. 6 shows the contributors to the environmental loads of the heat demand (Qd). The cogeneration module generates heat (operating with fuel Fc) and is complemented by heat produced by the auxiliary boiler Qa (which operates with fuel Fa). Fig. 7 shows the contributors to the environmental loads of the cooling demand Rd. Electricity purchased Ep from the grid and electricity generated Wc by the cogeneration module contribute to operate the electrical chiller. The auxiliary boiler operates on heat produced by the cogeneration module (Qc) plus heat produced by the auxiliary boiler (Qa). It is important to highlight the fact that the integration proposed in this paper (LCA and thermoeconomics into an EMIS) supports consideration of any environmental indicator, i.e., Eco-indicator 99 points, ozone layer depletion values (kgeq CFC-11 into air),

Table 5 CO2 emissions (kg CO2/kWh) associated with the operation of the trigeneration system. ExC1

ExC3

ExC7

ExC9

EMed EMqd EMrd

0.6502 0.2702 0.2434

0.6597 0.2273 0.3638

0.3773 0.2211 0.0755

0.4004 0.1497 0.2395

EMes EMwc EMqc EMwcc EMer EMqcc EMqa EMqr EMrq EMre

e 0.5445 0.2035 0.5445 0.6502 0.2035 0.3813 0.2702 0.4323 0.1300

e 0.6083 0.1478 0.6083 e 0.2273 e 0.2273 0.3638 e

1.0200 0.6160 0.1410 0.3773 0.3773 0.1410 0.3813 e e 0.0755

1.0200 0.6660 0.0973 0.4004 e 0.1497 e 0.1497 0.2395 e

Fig. 5. Emissions associated with the consumption of Ed in example ExC1.

M. Carvalho et al. / Environmental Modelling & Software 30 (2012) 71e80

79

term ecologically sustainable production, particularly in changing climatic conditions. From an environmental viewpoint, the purpose of installing a trigeneration system was to provide a reduction in the amount of environmental loads produced to attend a consumer center, which was achieved here. 5. Summary and conclusions

Fig. 6. Emissions associated with the consumption of Qd in example ExC1.

ecotoxicity values (kgeq triethylene glycol into water and soil), etc. Information shown in Tables 5 and 6 and Figs. 5e7 allows for the study of the distribution of consumed resources, of the cost formation process, and the generation of environmental loads throughout the productive system. This study emphasized the importance of integrating environmental information (obtained with the application of LCA techniques) and the philosophy of thermoeconomic analysis in Environmental Information Management Systems. Such integration was easily accomplished through an assistant implemented in Umberto. The assistant facilitated the registration and tracking of environmental impacts generated in each piece of equipment as well as the assessment of environmental loads to each output and internal flow, with the possibility of considering different approaches to the allocation issue. This combined approach allowed for a more precise understanding of operational activities. International markets are increasingly seeking an indication of environmental performance and the carbon emissions associated with products and services. The responsibility exists to ensure long

A trigeneration system was modeled in Umberto software, with the flexibility to purchase/sell electricity from/to the electric grid. The possibilities of wasting part of the cogenerated heat and operating an auxiliary boiler also existed. The network was calculated for specific energy service demands, identifying where emissions were generated, and tracking the emissions throughout the system. This study case concentrated on the issue of climate change and therefore considered the emissions of CO2. Flow analysis of individual production steps specific to operation made it possible to study the operational activities more precisely. This paper demonstrates that LCA can and has been formally combined with thermoeconomic analysis. The result is an ability to take thermoeconomics and LCA e and their tradeoff relationships e into account in product/process design decision making. This combined approach was integrated into an EMIS, identifying where environmental loads were generated and tracked their distribution to the final products of trigeneration systems. In an attempt to address the ongoing debate, an innovative environmental allocation method was proposed. The allocation proposal for trigeneration systems considered that environmental loads of the cogeneration module were distributed among the consumers of the final products, resulting in overall reduced emissions derived from the combined production. Such reductions were evaluated in proportion to the emissions associated with obtaining heat and electricity separately via conventional systems. The assistant developed also provided Umberto with the capability of calculating economic costs. A potent tool was obtained for energy, economic, and environmental analysis, to study the distribution of resources throughout productive systems as well as the formation process of costs and environmental loads. Effective environmental related strategies connect the reduction of emissions with a system’s operational strategy (consumption of resources). Therefore the usage of EMIS with double input from thermoeconomics and LCA tools could be promoted to (1) analyze the distribution of material and energy resources throughout a productive system; (2) serve the numerical registration and interpretation of environmental effects; (3) calculate unit economic costs and study the cost formation process; (4) identify the most environmentally beneficial among competing technologies; and (5) allow an emission-efficient economy to develop. By incorporating environmental information on the usage and consumption of resources into Umberto software, the approach of material flow networks gave insight on the environmental loads associated with each flow of the system. Thus, the consumers of a productive system can be made aware of the environmental loads associated with the consumption of each product (either internal or final). This information can be very useful for the introduction of strategies oriented to changes and improvements in the design and operation of productive systems as well as in consumption patterns and resource conservation, contributing to the development of a more sustainable economy. Acknowledgments

Fig. 7. Emissions associated with the consumption of Rd in example ExC1.

This work was developed within the framework of research project ENE2007-67122, funded in part by the Spanish Government (Energy program) and the European Union (FEDER program).

80

M. Carvalho et al. / Environmental Modelling & Software 30 (2012) 71e80

Special thanks are extended to the Umberto Competence Centre at HTW Berlin for the support given to Monica Carvalho during her research stay. References Abusoglu, A., Kanoglu, M., 2009. Allocation of emissions for power and steam production based on energy and exergy in diesel engine powered cogeneration. Energy & Fuels 23, 1526e1533. Awasthi, A., Chauhan, S.S., 2011. Using AHP and Dempster-Shafer theory for evaluating sustainable transport solutions. Environmental Modelling & Software 26 (6), 787e796. Carvalho, M., et al., 2010. Distribution of resources and allocation of environmental loads applied to trigeneration systems. In: Proceedings of the 23nd International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact of Energy Systems, Lausanne, Switzerland. Carvalho, M., 2011. Thermoeconomic and environmental analyses for the synthesis of polygeneration systems in the residential-commercial sectorÓ. Ph.D. thesis, University of Zaragoza, Spain. (accessed 01.11.11). Carvalho, M., Serra, L.M., Lozano, M.A., 2011a. Geographic evaluation of trigeneration systems in the tertiary sector. Effect of climatic and electricity supply conditions. Energy 36 (4), 1931e1939. Carvalho, M., Serra, L.M., Lozano, M.A., 2011b. Optimal synthesis of trigeneration systems subject to environmental constraints. Energy 36 (6), 3779e3790. Chicco, G., Mancarella, P., 2008. Assessment of the greenhouse gas emissions from cogeneration and trigeneration systems. Part I: Models and indicators. Energy 33 (3), 410e417. CHOSE, 2001. Energy savings by combined heat cooling and power plants (CHCP) in the Hotel Sector. Final Report, SAVE Project. (accessed 01.11.11). Colonna, P., Gabrielli, S., 2003. Industrial trigeneration using ammoniaewater absorption refrigeration systems (AAR). Applied Thermal Engineering 23 (4), 381e396. COM 2002/91/EC, Directive on the energy performance of buildings. Official Diary of the EC, 4.1.2003, L1, 65e71. Díez, E., McIntosh, B.S., 2009. A review of the factors which influence the use and usefulness of information systems. Environmental Modelling & Software 24 (5), 588e602. Ecoinvent v2.0 Database, 2007. Swiss Centre for Life Cycle Inventories, Dübendorf, Switzerland. El-Sayed, Y.M., Evans, R.B., 1970. Thermoeconomics and the design of heat systems. Transactions of ASME Journal of Engineering for Power 92, 27e34. El-Sayed, Y.M., Gaggioli, R.A., 1989. A critical review of second law costing methods: Parts I and II. ASME Journal of Energy Resource Technology 111, 1e15. Gaggioli, R., 1983. Second law analysis for process and energy engineering. In: Efficiency and Costing, ACS Symposium Series, vol. 235, pp. 3e50. Guinée, J.B. (Ed.), 2002. Handbook on Life Cycle Assessment: Operational Guide to the ISO Standards. Kluwer Academic Publishers, Dordrecht; Boston. IDEMAT Database, Section for Environmental Product Development, Faculty of Industrial Design, 2001. Delft University of Technology, Delft, The Netherlands. (accessed 01.11.11). Lingo, 2008. The Modeling Language and Optimizer. Chicago, USA (accessed 01.11.11). Lozano, M.A., Valero, A., 1993. Theory of the exergetic cost. Energy 18 (9), 939e960. Lozano, M., Carvalho, M., Ramos, J.C., Serra, L.M., 2009a. Thermoeconomic analysis of simple trigeneration systems. International Journal of Thermodynamics 12 (3), 147e153.

Lozano, M.A., Ramos, J.C., Carvalho, M., Serra, L.M., 2009b. Structure optimization of energy supply systems in tertiary sector. Energy and Buildings 41 (10), 1063e1075. Lozano, M.A., et al., 2009c. Energy cost analysis of simple trigeneration systems under variable operation conditions. In: Proceedings of the 22nd International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact of Energy Systems Foz do Iguaçu, Brazil. Maglorie, P., et al., 2002. Economie d’énergie en trigénération. International Journal of Thermal Sciences 41, 1151e1159. Marimón, M.A., Arias, J., Lundqvist, P., Bruno, J.C., Coronas, A., 2011. Integration of trigeneration in an indirect cascade refrigeration system in supermarkets. Energy and Buildings 43 (6), 1427e1434. Möller, A., 2010. Material and energy flow-based cost accounting. Chemical Engineering & Technology 33 (4), 567e572. Phylipsen, G., et al., 1998. Handbook on International Comparisons of Energy Efficiency in the Manufacturing Industry. Utrecht University Press, the Netherlands. Pokorný, J., 2006. Database architectures: current trends and their relationships to environmental data management. Environmental Modelling & Software 21 (11), 1579e1586. Porter, M., van der Linde, C., 1995. Green and competitive: ending the stalemate. Harvard Business Review 73, 120e134. Rosen, M., 2008. Allocating carbon dioxide emissions from cogeneration systems: descriptions of selected output-based methods. Journal of Cleaner Production 16 (2), 171e177. SAIC e Scientific Applications International Corporation, 2006. Life Cycle Assessment: Principles and Practice. National Risk Management Research Laboratory, U.S. Environmental Protection Agency, Cincinnati, Ohio, US. Serra, L.M., et al., 2009. Polygeneration and efficient use of natural resources. Energy 34 (5), 575e586. Simapro, 2008. Life Cycle Assessment Software. PRé Consultants, Amersfoort, The Netherlands. (accessed 01.11.11). Summerheat Project, 2008. Cooling in Summer with Heat from CHP Technologies. Euro Heat & Power English Edition, vol. 5. . http://www.eu-summerheat.net/, III/2008, pp. 22e25. TRIGEMED, 2003. Promotion of Tri-Generation Technologies in the Tertiary Sector in Mediterranean Countries. Survey Report, SAVE Project. . http://www.managenergy.net/resources/368. Turner, G.M., Hoffman, R., McInnis, B.C., Poldy, F., Foran, B., 2011. A tool for strategic biophysical assessment of a national economy e the Australian stocks and flows framework. Environmental Modelling & Software 26 (9), 1134e1149. Umberto, version 5.0, ifu Hamburg GmbH and ifeu GmbH, Hamburg. , 2006 Germany (accessed 01.11.11). Viaggi, D., Bartolini, F., Raggi, M., 2009. Combining linear programming and principaleagent models: an example from environmental regulation in agriculture. Environmental Modelling & Software 24 (6), 703e710. Viere, T., et al., 2010. J. Verbund-simulation e strategic planning and optimization of integrated production networks. Chemical Engineering & Technology 33 (4), 582e588. Wang, J.-J., Jing, Y.-Y., Zhang, C.-F., Shi, G.-H., Zhang, X.-T., 2008. A fuzzy multi-criteria decision-making model for trigeneration system. Energy Policy 36 (10), 3823e3832. Wohlgemuth, V., et al., 2006. Combining discrete event simulation and material flow analysis in a component-based approach to industrial environmental protection. Environmental Modelling & Software 21 (11), 1607e1617. Xu, Z., Coors, V., 2011. Combining system dynamics model, GIS and 3D visualization in sustainability assessment of urban residential development. Building and Environment 47, 272e287.