تعداد نشریات | 43 |

تعداد شمارهها | 1,612 |

تعداد مقالات | 13,135 |

تعداد مشاهده مقاله | 28,991,189 |

تعداد دریافت فایل اصل مقاله | 11,583,773 |

## Thermodynamic simulation and economic modeling and optimization of a multi generation system partially fed with synthetic gas from gasification plant | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Gas Processing Journal | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

مقاله 5، دوره 5، شماره 2، اسفند 2017، صفحه 49-68 اصل مقاله (1.25 M)
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

نوع مقاله: Research Article | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

شناسه دیجیتال (DOI): 10.22108/gpj.2018.110943.1025 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

نویسندگان | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Mohammad Rahimi^{*} ؛ Mohammad Hosein Hamedi؛ Majid Amidpour
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

^{}Faculty of Mechanical Engineering-Energy Division, K.N. Toosi University of Technology | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

چکیده | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

This paper presents thermodynamic simulation, economic modeling and annual profit optimization of a multi generation system which produces both power and fresh water. The fuel of the combined system is natural gas plus synthesis gas which is produced in biomass gasification reactor. In order to evaluate thermodynamic performance of the biomass gasification reactor, visual simulation software was developed in C# programming language. The multi generation system is analyzed both with inlet air cooling and without inlet air cooling. The final results show that the total cost of produced power is 0.0286 $/kWh and total cost of produced water is 0.7408286 $/m^{3}. Also the total annual profit which comes from selling power and water to the market is 35.103 M$ and the CHP efficiency is 67.08. Optimization of the configuration is carried out once the simulation phase is finished. The optimization results in 10.5% increase in total annual profit and 6.6% increase in CHP efficiency. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

کلیدواژهها | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Synthetic gas؛ Desalination system؛ Power and water cost؛ Net annual profit؛ Gasification؛ Genetic Algorithm | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

اصل مقاله | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Multi generation thermal systems have drawn great attention nowadays. Multi generation means the combined production of heat, power, water, cooling, liquid fuel, etc., for consumption within a site. Heat can have several uses. For example, it can be used as the motive steam for thermal desalination systems (MSF or MED) or as the source of heat for absorption chiller. Whilst power and heat are provided by the turbine and the exhaust gases refrigeration could be obtained in two different ways, either by using an absorption system in combination with low grade heat or by using an electrically driven compression system. The use of one or another will depend on the process heat/power ratio needs and specific site characteristics. Thermal desalination is among the most useful applications of multi generation. Many researchers have studied thermal desalination from thermodynamic and economic points of view. Sayyaddi et al. (Sayyaadi&Ghorbani, 2018) introduced a systematic approach for the design of Stirling-desalination system which was found to be a reliable option for the small-scale power-water production. The proposed system could deliver 2.58 kW of the electric power as well as 23.3 m Shih (Shih, 2005) evaluated the technologies of thermal desalination using low-grade heat present in a sulfuric acid plant, Mabrouk (Mabrouk, 2013) explored a techno-economic comparison between long tube (LT) and cross tube (CT) bundles of MSF evaporator for a unit production of equal and greater than 20 MIGD, Fiorini and Sciubba (Fiorini&Sciubba, 2005) adapted a modular simulation code, CAMEL™, developed by the University of Roma1, to include the capability to perform a thermo economic analysis of a MSF desalination plant (in addition to the thermodynamic and exergetic analyses) and Nafey et al. (Nafey, Fath, &Mabrouk, 2006) did a number of comparisons for Multi Effect Evaporation (MEE) and hybrid Multi Effect Evaporation-Multi Stage Flash (MEE-MSF) systems using the exergy and thermo economic analysis. Recently, some studies have been performed to produce fresh water with the use of renewable energy resources. Mentis et al. (Mentis et al., 2016) developed a tool for designing and optimally sizing desalination and renewable energy units. Ghaffour et al. (Ghaffour et al., 2014) worked on developing new desalination processes, adsorption desalination (AD) and membrane distillation (MD), which can be driven by waste heat, geothermal or solar energy. They constructed a demonstration solar-powered AD facility. A life cycle assessment showed that its specific energy consumption is less than 1.5 kWh per cubic meter of desalinated water, which is far less than the energy consumption of conventional desalination methods. On the other hand, in the field of desalination system’s parameter and configuration optimization, some research activities have been done. Kwon et al. (Kwon, Won, & Kim, 2016) developed a superstructure model using Mixed Integer Linear Programming to determine the optimal configuration of a renewable-based power supply.
Diverse economic factors such as transmission and reclamation costs were considered to ensure minimal cost while satisfying electricity demands, Ansari et al. (Ansari, Sayyaadi, &Amidpour, 2011) considered a typical 1000 MW Pressurized Water Reactor (PWR) nuclear power plant coupled to a multi effect distillation desalination system with a thermo-vapor compressor (MED–TVC) for optimization. Shakib et al. (Shakib, Amidpour, &Aghanajafi, 2012) did an optimization study for a combined system of gas turbine, Heat Recovery Steam Generator and MED desalination unit in view of three approaches, Kamali et al. (Kamali, Abbassi, SadoughVanini, &SaffarAvval, 2008; Kamali&Mohebinia, 2008), did a parametric optimization analysis of a multiple effect desalination system with thermal vapor compression (MED-TVC) process to increase gain output ratio (GOR), Ameri et al. (Ameri, Mohammadi, Hosseini, &Seifi, 2009) studied the effects of different design parameters such as number of evaporation effects, inlet steam pressure, temperature difference of the effects, etc. on MED system specifications, Mehrpooya et al. (Mehrpooya, Ghorbani, Jafari, Aghbashlo, &Pouriman, 2018) investigated a novel hybrid model based on neural network. The proposed model was a combination of Group Method of Data Handling type neural networks and Genetic Algorithm. The Genetic algorithm was used to optimize the correlation parameters to improve the accuracy of model, Agashichev and El-Nashar (Agashichev& El-Nashar, 2005) developed a system of models for the techno-economic evaluation of a triple hybrid, reverse osmosis (RO), multistage flush (MSF) and power generation process. Researchers, nowadays, have been doing extensive studies about the systems that can produce more than two or even three forms of useful products. In fact, the trend is also accelerating toward using renewable forms of input energy. Malik et al. (Malik, Dincer, & Rosen, 2015) developed a renewable energy-based multi-generation system and studied it both energetically and exergetically. They employed Two renewable sources of energy, biomass and geothermal, to deliver five useful outputs. They found that the energy and exergy efficiency of the entire system is 56.5% and 20.3% respectively. Shariatiniasar et al. (ShariatiNiasar et al., 2017) studied a cogeneration system of four useful outputs including power, heating, cooling and liquid fuels with the use of gasification of coal. Ghorbani et al. evaluated an integrated system for co-production of LNG and NGL, based on MFC and absorption refrigeration systems.They found the highest and the lowest exergy destruction parts of the system and concluded that the forth compressor has the highest exergy destruction cost. Salehi et al. did an optimization on an integrated heat and power system which were part of a distillation column sequence. The results showed that a large amount of power can be produced between the columns due to having high flow rate flows between the columns. Huang et al. (Huang et al., 2013) did a simulation and techno-economic analysis of small scale biomass trigeneration system. The study investigates the impact of different biomass feedstock on the performance of trigeneration plant. The results specified the maximum efficiencies and the best breakeven electricity selling prices. This paper presents a new evaluation system for the comparison of different configurations of combined cycles of power, cooling and desalinated water production. A complete thermodynamic and economic modeling and optimization procedure is used for this purpose.
The multi generation system is composed of a gas turbine which is a Siemens V94.2 (with nominal power output of 148.8 Mw in 15 °C inlet air temperature) and the desalination system which is chosen to be MSF (Multi Stage Flash) system. The high pressure steam produced in HRSG is used for power production in the back pressure steam turbine, but the low pressure steam, which is saturated, is mixed with the low pressure steam exiting the back pressure steam turbine to feed the brine heater of MSF system. The schematic of the system and main inputs are shown in figure 1 and table 1 respectively. Simulation is done both with inlet air cooling and without inlet air cooling.
In this section the thermodynamic modeling of the above mentioned configuration is outlined. The first stage in every modeling process is to identify the inputs and the outputs of the problem.Generally speaking, the inputs fall into two categories:system parameters(which are considered fixed values during modeling and simulation)and decision variables(which
are changeable within a certain limit). The outputs are actually the results obtained from inputs, using the fundamental laws (in this paper, the first and the second laws of thermodynamics). System parameters as the major assumptions that are considered for modeling of the system are mentioned in table 2. For example, in this work, a Siemens V94.2 gas turbine was chosen and the ambient design temperature decided to be 35 C Other inputs necessary for carrying out the modeling of the configuration were presented in the previous sections, among them the pressure and temperature of High Pressure (HP) and Low Pressure (LP) steam of HRSG and the approach and pinch temperature are of great importance. Dependent variables are those variables that will be obtained after running the simulation code. They are actually the outputs of the system and their values are dependent on the values of parameters and decision variables of the system. Major dependent variables are as follows: 1- Power produced by the gas turbine 2- Power produced by the steam turbine 3- The net power output of the cycle 4- The mass flow rate of the fuel 5- The volume flow rate of desalinated water 6-Total capital investment of power production 7- Total capital investment of water production 8- The net annual profit of the plant
In order to evaluate thermodynamic performance of the biomass gasification reactor, visual simulation software was developed in C# programming language. The software is composed of different tabs which have specific purposes and perform required calculations. Two modes of simulations are available for both fixed and fluidized bed reactors; both of them solve a set of non-linear equations to get the desired outputs. In the first mode, the reaction temperature is given and the program will calculate the required amount of air to satisfy energy balance in the reactor. In the second mode, the amount of air injected to the gasifier is given and the program uses an iterative approach to calculate final reaction temperature. In this mode, giving an initial reaction temperature is necessary for starting the calculations. The software uses extended Newton-Raphson approach to solve the set non-linear equations and obtain the producer gas composition. First, the biomass inlet characteristics (atomic composition, LHV, moisture, biomass flow, etc.), environmental conditions (temperature and pressure), gasifier parameters (gasifier temperature or mole of injected air) and other parameters (gasifier operating pressure, heat loss) are introduced into the program. In second stage, the product gas composition is calculated by stoichiometric equilibrium model. Flowchart of the calculation approach in the first and second mode of computation can be observed in figure 2 and 3. Procedure of formulation for the second mode of simulation (air injected to the gasifier given) with the use of stoichiometric method is as follows, it should be noted that formulation of the first mode of simulation (reaction temperature given) is quite similar with marginal changes: First off, starting from the mass fractions of carbon, hydrogen, oxygen, nitrogen and sulfur (CHONS) in the biomass and the relative mass of the moisture, the substitution fuel and the molar water content can be evaluated. In the second stage, the composition of the producer gas is estimated, using the initial value of reaction temperature and calculation of the equilibrium constants. Then the reaction temperature corresponding to the actual producer gas composition is calculated, equating the enthalpy of the entering biomass and moisture and the enthalpy of the producer gas. Using the calculated reaction temperature, the input of the next composition calculation is formed, iterating the process until chemical and thermodynamic equilibrium have been reached.
Once the final producer gas composition and its corresponding reaction temperature are obtained, other outputs of the gasification process including the heating value of the producer gas, Cold Gas Efficiency, etc. can be derived. Estimating the composition of producer gas is based on chemical equilibrium between different species, neglecting tar content in the producer gas. The reaction, in its general form, can be written as (Melgar, Pérez, Laget, &Horillo, 2007): The variable x corresponds to the molar quantity of air used during the gasifying process and is one of the inputs of the simulation. The value of m, p, q, and r can be calculated from weight percent of Hydrogen, Carbon, Oxygen, Nitrogen, Sulfur and their molecular weights. Also, from molecular weight of biomass and water and the relative moisture of biomass, the value of ω can be calculated. Writing the atomic balance for C, H, O, N and S, respectively and assuming that no oxygen will be present in the producer gas, following six equations are derived (Melgar et al., 2007).
In order to solve the above system of equations, to find out eight unknown variables, two more equations are needed. The first one is reduction of hydrogen to methane in reduction zone. The second one is known as the water gas shift reaction, which is the equilibrium between CO and H2 in the presence of water.
The corresponding equilibrium constants of the above mentioned equations can be obtained from either the molar composition of syngas or Gibbs free energy. If the second is substituted form of equilibrium constant equation (Gibbs free energy) into the first one (molar composition), the complementary equations will be found (Melgar et al., 2007):
In which Gibbs free energy can be calculated from (Melgar et al., 2007):
Thermodynamic properties are extracted from NIST-JANAF thermochemical tables (Chase, 1998). Two of the mentioned system of equations has non-linear structure, so an extended Newton- Raphson scheme is employed in order to solve the system of equations. Once the above mentioned system of equations is solved, the syngas composition will be determined at initial reaction temperature. Knowing the syngas composition, corresponding reaction temperature can be computed, which is, in turn, the initial reaction temperature of next iteration. The reaction temperature of next iteration can be estimated using the first law of thermodynamic according to the following equations (Melgar et al., 2007):
There are two approaches for Gas Turbine modeling. The first one is using classical thermodynamics laws, namely the first and second laws of thermodynamics and utilizing the concepts of isentropic efficiency, etc.(Bejan A, 1996). In this way, some assumptions are needed. For the sake of simplicity, the fuel is considered Methane, physical properties of all streams are calculated in mean inlet and outlet temperature and the air and combustion products are treated as ideal gases. Referring to figure 2, the temperature of air leaving the compressor in ideal condition would be (Bejan A, 1996):
In which, K is the ratio of the specific heat capacity at constant pressure to the specific heat capacity at constant volume. Specific heat capacity of every component of air (or any other ideal gas mixture) can be calculated according to the following equation in which a, b, c and d are constant coefficients for each component (Bejan A, 1996).
Actual air temperature leaving the compressor can be calculated utilizing compressor isentropic efficiency (Sonntag, Borgnakke, Van Wylen, & Van Wyk, 1998):
In a similar way, the ideal and actual temperature of combustion products leaving the turbine can be computed (Sonntag et al., 1998):
Knowing the composition of ambient air (mentioned in table 5) and for complete combustion of methane, the chemical equation takes the form (Bejan A, 1996):
In which is fuel to air mole fraction. By balancing the two sides of equation, the mole fraction of the products components will be obtained (Bejan A, 1996).
So the molar analysis of combustion products is fixed once the fuel to air ratio ( ) has been determined. The fuel to air ratio ( ) can be obtained from an energy rate balance around combustion chamber (Bejan A, 1996):
Since , the fuel and air mass flow rates are related by:
This approach, for modeling a certain type of gas turbine, like V94.2, will lead to small deviation from the actual performance of the gas turbine, since the manufacturer has probably used special assumptions for modeling and constructing the gas turbine. The second approach which is more accurate than the first is to use the manufacturer’s graphs and tables and use regression if required to get our desired outputs. After finding the temperature, flow rate and composition of the flue gas leaving the gas turbine, the next step is finding the enthalpy and entropy of the water streams. But first it is needed to do the pressure analysis of the water side of the cycle.
According to the assumptions for the values of pressure drop in various portions of HRSG and the known pressures of the cycle, which were mentioned before, the pressure of every single stream can be found (Sonntag et al., 1998).
It should be noted that for finding steam properties, XSteam code, which is available at (Holmgren, 2006), has been used.
Having completed the pressures analysis and before finding the mass flow rates, the enthalpies and entropies of the all streams were found. For some streams, it is completely straightforward, since the pressure and temperature of that stream is known. Using XSteam code, it is easy to find the enthalpy and entropy of all streams.
For some other streams, for example the exit of pumps or steam turbine it is necessary to use the isentropic efficiency formula. For pump1 it is as follows (Sonntag et al., 1998):
Since h The isentropic efficiency formula for steam turbine is (Sonntag et al., 1998):
And by the use of the mentioned formula, the enthalpy, entropy, temperature and wetness of the leaving streams of steam turbines can be found. Generally, to find out the steam flow rates, a proper control volume around one component or two components, depending on the configuration of the cycle should be considered. Next, by employing the laws of conservation of mass and energy and simultaneous solving of the equations, the desired mass flow rate will be obtained. To find HP steam mass flow rate, it is needed to consider a control volume encompassing both super heater and evaporator1. Writing the conservation of energy for this control volume will lead to:
T For finding LP mass flow rate, a control volume around evaporator1 would be adequate.
In order to find the power output (for steam turbine) or power input (for pumps), the steam (or water) flow rates and enthalpies upstream and downstream of the component as well as the mechanical efficiencies are needed. The former was calculated in the previous sections. The latter is considered a fixed parameter and was mentioned in table1 for each part. With reference to figure1, employing two equations as follows will lead to the power output of the steam turbine and power input of water pump:
For evaluation of thermal performance, a mathematical model is developed by applying mass and energy conservation laws to the flashing stages and condenser (Hisham T. El-Dessouky, 2002). The final objective is to obtain the total produced desalinated water per day. For this purpose the mass flow rate and the temperature and pressure of the motive steam is needed. These values were found in the previous section. In this study, brine circulation MSF process has been chosen. The following assumptions are considered in this regard: Distillate product is salt free, Specific heat at constant pressure, Cp, for all liquid streams, brine, distillate and seawater is constant and equal to 4.18 kJ/kg C, Sub cooling of condensate or superheating of heating steam has negligible effect on the system energy balance and The heat losses to the surroundings are negligible because the flashing stages and the brine heater are usually well insulated. Schematic of the brine circulation MSF process is shown in Figure 3 below. In addition to the above mentioned assumptions, some key parameters of the MSF system should be known which were indicated in table 2. The procedure of modeling is as follows:The overall material balance equation of the system can be arranged to obtain the expression for the total feed flow rate in terms of the distillate flow rate which is given by equation 41 (Hisham T. El-Dessouky, 2002):
Where M is the mass flow rate and the subscript b, d and f defines the brine, distillate, and feed and X is the salt concentration. This equation assumes that the distillate is salt free. The temperature distribution in the MSF system is defined in terms of four temperatures; these are the temperature of the steam, T A linear profile for the temperature is assumed in the stages and the condensers, the temperature drop per stage, ΔT, is ΔT = (T By performing an energy balance on stage i, Assuming the temperature difference, T This gives the general relation for the seawater temperature in the rejection section T Using the conservation of energy within each stage and some mathematical work, the total distillate flow rate is obtained by summing the values of Di for all stages (Hisham T. El-Dessouky, 2002). (D
In which M
Where C The Gain Output Ratio (GOR) is a measure of water produced relative to steam consumed. Specifically, Gain Output Ratio (GOR) is defined as: (kilograms desalinated water produced) / (kilograms steam condensed)
The main objective of economic modeling is to calculate the cost of produced power and desalinated water as well as the total annual profit of the plant. Cost of produced power and water is a function of the net power and water output, capital investment of the plant, fixed and variable operating and maintenance cost, fuel price, etc. since each system is different from other systems from several aspects, its cost of power and water would be different. Below, the process of calculating the two is outlined.
The cost of produced power is primarily composed of four parts. All of which would be in $/kWh.
Capital investment of a thermal system is composed of several parts, including purchase equipment costs of specialized equipments (gas turbine, steam turbine, HRSG, condenser, gasifier, etc.), purchase equipment costs of other equipments (pumps, tanks, cooling tower, etc.), civil and structural works, mechanical works (equipment erection and piping), electrical and wiring works and plant startup and engineering services. Purchase equipment costs of the plant components are listed in table 3. Cost data are often presented as cost versus capacity charts, or expressed as a power law of capacity(Smith, 2005).
Where C Such data can be brought up-to-date and put on a common basis using cost indexes. Commonly used indices are Marshall and Swift, published in Chemical Engineering magazine
In which CAP is the total money invested for power generation components in million dollars, CRF is the Capital Recovery Factor and is obtained by the following formula (Smith, 2005):
“ ” and “ ” are respectively the number of years that the plant is in operation and interest rate. P is the net power output of the plant in Mw and is the availability of the power generation plan which is defined as the ratio of the total days in year which the plant is in operation and produce power to the total days of a year. Refer to table 7 for the detail results.
Fixed operating and maintenance cost of different power generation plants usually ranges from 10 to 40 $/KW-year. For the configuration presented in this paper, this cost is as follows ("Thermo flow (GTpro module)"):
Since all the costs should be congruent ($/kWh) the following formula is used:
This cost is usually expressed in $/kWh. For a combined cycle power plant, similar to the herein scenario, a reasonable value is 0.002 $/kWh according to the information of the plants currently in operation in Iran.
Considering the price of the fuel (methane) to be 0.003 $/MJ LHV and using the following formula, the power cost, resulted from fuel price, can be calculated (Smith, 2005):
HR is the Heat Rate in MJ/kWh and is calculated according to the following formula (Sonntag et al., 1998):
is the fuel mass flow rate in Kg/s and LHV is the fuel Low Heating Value in KJ/Kg and P is the net power output in MW. Thus the total cost of produced power in $/kWh is the sum of the four parts mentioned above.
The cost of produced water is composed of three parts. All parts should be in $/m
Capital investment of a desalination system, like power generation system, is composed of several parts, the most important of which are brine heater, flashing stages, transferring pumps and related piping and civil and structural works. Knowing the total capital investment of the desalination system in M$, the following formula is used to find the corresponding cost term:
The fixed operating cost of a desalination system, according to the technical reports, is 0.1159 $/m The total cost of produced water is the sum of the three parts mentioned above.
Having finished the process of thermodynamic and economic modeling, the final results can be presented. In the table 3, the major data for economic calculations are mentioned. In table 4, the detailed results of thermodynamic and economic modeling are shown. As it can be seen from table 4, the total cost of produced power is 0.0286 $/kWh and total cost of produced water is 0.7408286 $/m In table 5, a precise comparison between the outputs of the combined system, with and without inlet air cooling is shown. CHP efficiency, Total annual profit and the total cost of produced power and water are among the comparison parameters. Table 5 shows that the configuration with inlet air cooling, in which the low pressure steam generated in HRSG, plus the low pressure steam leaving the back pressure steam turbine is used to feed the brine heater, has the higher total annual profit. This configuration is optimized in the next section. The details of the optimization process are given in the next section.
In order to achieve the optimal value of decision variables, an optimization algorithm should be employed. Although gradient descent methods are the most elegant and precise numerical methods to solve optimization problems, however, they have the possibility of being trapped at local optimum points depending on the initial guess of solution. Stochastic optimization methods such as genetic algorithm (GA) and Particle Swarm Optimization (PSO) seem to be promising alternatives for optimization problems similar to this paper’s configuration. In general, they are robust search and optimization techniques, able to cope with ill-defined problem domain such as multimodality, discontinuity and time-variance (Shakib et al., 2012). GA is a population based optimization technique that searches the best solution of a given problem based on the concepts of natural selection, genetics and evolution (Holland, 1992). PSO is a heuristic population based optimization algorithm simulating the movement and flocking of birds (Modares & Naghibi Sistani, 2011). In this work, genetic algorithm has been chosen as the optimization method with an economic objective function, namely the total annual profit.
The first stage in an optimization problem is to fully specify system parameters (which are considered fixed values during optimization), decision variables (which are changeable within a certain limit) and dependent variables (which are actually the outputs of the problem). During the optimization process, objective function, that is the most important dependent variable, should be maximized (or minimized, depending on its nature) with changing the decision variables within their limits. These limits are dependent on the physical, mechanical and thermodynamic constraints. In the optimization phase, the best values of these variables for satisfying the objective function will be chosen by the optimizer. In this paper, in order to optimize the configuration, the following variables are chosen as decision variables: 1- The pressure of High Pressure (HP) steam of HRSG 2- The temperature of High Pressure (HP) steam of HRSG 3- The pressure of Low Pressure (LP) steam of HRSG 4- Pinch temperature difference of the first evaporator 5- Approach temperature of the first evaporator 6- Pinch temperature difference of the second evaporator 7- Approach temperature of the second evaporator In table 6, the range of change of decision variables is shown. The major system parameters and dependent variables of this configuration were mentioned in previous sections.
Once the thermodynamic and economic model of the combined system is built, the corresponding code (in MATLAB programming language) is created, the decision variables are decided and the physical constraints are exerted, the optimization process starts using the genetic algorithm toolbox of MATLAB. The optimum condition, in which the annual net profit is maximized, is achieved after several iterations. Table 8 presents the values of decision variables as well as dependent variables in the base case and optimum case. The table shows that the optimization, results in 10.5% increase in total annual profit and 6.6% increase in CHP efficiency.
Economic parameters (for example fuel price, power sale price, water sale price, utilization years, etc.) have a great impact on the final results, including the net annual profit. In the modeling phase as well as the optimization phase, specific values for economic parameters are chosen, for example, the fuel price considered to be 3 $/GJ (LHV), Interest rate to be 15% and so on (refer to table 6). Furthermore, thermodynamic variables (for example pressure of HP and LP steam) have the similar effects on final results and are of great importance. In this section, the impact of important economic and thermodynamic parameters on the cost of produced power and the net annual profit of the optimum case is evaluated. In the figure 6 and figure 7, the cost of produced power and the net annual profit as a function of fuel price and utilization years are shown respectively. The range of change of fuel price is from 1 $/GJ (LHV) to 6 $/GJ (LHV). As it can be seen, for fuel price of 1 $/GJ (LHV) and utilization years of 10, 15 and 20, the cost of produced power are 1.85, 1.71 and 1.66 cent/kWh respectively. It is predictable that by increasing the years of utilization of the combined system, the cost of produced power will decrease. It is also observable that the annual profits, for fuel price of 1 $/GJ (LHV) and utilization years of 10, 15 and 20, are 50.55, 54.25 and 55.72 M$ respectively. These figures also show that, for utilization years of n = 15, increasing the fuel price to 6 $/GJ (LHV) will increase the cost of produced power to 3.94 Cent/kWh and decrease the net annual profit to 23.26 M$. Figure 8 shows the annual profit as a function of PSP (power sale price) and WSP (water sale price). Apparently, increasing the two will increase annual profit with a trend shown in the figure. As it is shown in the figure 8, for power sale price of 3.5 cent/kWh, and water sale price of 0.5, 0.8, 1.11 and 1.5 $m
Next figure, figure 9, shows the annual profit as a function of utilization years and interest rate. It is observed that increasing the interest rate will highly decrease the total annual profit but increasing the utilization years, especially from 20 to 30, will have a marginal increase effect on annual profit. In fact, for 10 years of utilization, increasing the interest rate from 10% to 25% will decrease the annual profit from 42.92 M$ to 27.51M$, but increasing the years of utilization from 20 to 30 will increase the annual profit from 48.87 M$ to 50.36 M$ (for fixed interest rate of 10%). In figure 10, the annual profit as a function of pressure and temperature of HP steam is shown. As it can be seen, increasing the steam pressure initially increases the annual profit, but further increase in steam pressure will no longer increase the annual profit, but instead, marginally decrease it. In fact, figure 10 shows that increasing the high pressure steam pressure form 50 bar to 75 bar (for high pressure steam temperature of 450 ˚C) will increase the annual profit for 135000 $, but increasing this pressure from 75 bar to 100 (again for high pressure steam temperature of 450 ˚C) will conversely decrease the annual profit for 4000$. On the other hand, increasing the high pressure steam from 450 ˚C to 510 ˚C will increase the annual profit on almost a regular basis. Figure 11 demonstrates the effect of LP steam pressure on the annual profit for three HP steam temperatures. It is observed that as a result of increasing LP steam pressure, the annual profit will decrease for all three HP steam temperatures. It is observed in figure 11 that the rate of decrease is quite the same for all three HP steam temperatures. It is observable that by increasing the LP steam pressure from 2 bar to 4 bar (for high pressure steam temperature of 450 ˚C), there would be a decrease of 2,274,000 $ in annual profit which a considerable value. On the other hand, for fixed LP pressure of 2 bar, increasing the HP steam temperature from 450 ˚C to 510 ˚C will only increase the annual profit for 275,500 $.
In this paper, a novel thermal system for combined production of power and desalinated water was modeled and analyzed from both thermodynamic and economic points of view. The impact of inlet air cooling on thermodynamic and economic performance of the configuration was also investigated. Optimization of the configuration was carried out next. The most important outputs of the modeling were the net power output, total desalinated water produced, CHP efficiency, total power cost, total water cost and finally net annual profit. The results of modeling showed that the configuration had a relatively high total annual profit and the CHP efficiency. In the next section, the configuration was optimized through genetic algorithm method. Total annual profit of the combined system was chosen as objective function. The optimization process resulted in 10.5% increase in total annual profit and 6.6% increase in CHP efficiency. Evaluation of the impact of important economic and thermodynamic parameters on objective function was done in the last section. It showed the effects of fuel price, power sale price, water sale price, utilization years, interest rate and HP and LP conditions on net annual profit and cost of produced power.
^{2}) adjacent fins of exchanger
^{2})
_{ }equipment cost with capacity Q_{B }(Base capacity)
_{ }Specific heat capacity at constant pressure (KJ/Kg)
Gibbs free energy (KJ/Kmol)
Molar enthalpy (KJ/Kmol) enthalpy of formation (KJ/Kmol)
Mass flow rate (Kg/s)
^{th} component of producer gas
Time rate of heat (KJ/s)
Time rate of work (KJ/s)
Subscript
Greek letters
latent heat
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

مراجع | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Agashichev, S. P., & El-Nashar, A. M. (2005). Systemic approach for techno-economic evaluation of triple hybrid (RO, MSF and power generation) scheme including accounting of CO2 emission. Energy, 30(8), 1283-1303. doi:http://dx.doi.org/10.1016/j.energy.2004.02.002
Al-Hengari, S., El-Bousiffi, M., & El-Mudir, W. (2005). Performance analysis of a MSF desalination unit. Desalination, 182(1–3), 73-85. doi:http://dx.doi.org/10.1016/j.desal.2005.03.010
Alasfour, F. N., Darwish, M. A., & Bin Amer, A. O. (2005). Thermal analysis of ME—TVC+MEE desalination systems. Desalination, 174(1), 39-61. doi:http://dx.doi.org/10.1016/j.desal.2004.08.039
Alhazmy, M. M. (2014). Economic and thermal feasibility of multi stage flash desalination plant with brine–feed mixing and cooling. Energy, 76(0), 1029-1035. doi:http://dx.doi.org/10.1016/j.energy.2014.09.022
Ameri, M., Mohammadi, S. S., Hosseini, M., & Seifi, M. (2009). Effect of design parameters on multi-effect desalinationsystem specifications. Desalination, 245(1–3), 266-283. doi:http://dx.doi.org/10.1016/j.desal.2008.07.012
Ansari, K., Sayyaadi, H., & Amidpour, M. (2011). A comprehensive approach in optimization of a dual nuclear power and desalination system. Desalination, 269(1–3), 25-34. doi:http://dx.doi.org/10.1016/j.desal.2010.12.035
Bejan A, T. G., Moran M. (1996). Thermal Design and Optimization: Wiley.
Bridgwater, A. V., Toft, A. J., & Brammer, J. G. (2002). A techno-economic comparison of power production by biomass fast pyrolysis with gasification and combustion. Renewable and Sustainable Energy Reviews, 6(3), 181-246. doi:https://doi.org/10.1016/S1364-0321(01)00010-7
Carapellucci, R., & Giordano, L. (2013). A comparison between exergetic and economic criteria for optimizing the heat recovery steam generators of gas-steam power plants. Energy, 58, 458-472.
Chase, M. (1998). NIST—JANAF Thermochemical Tables (Journal of Physical and Chemical Reference Data Monograph No. 9). American Institute of Physics.
El-Sayed, Y. M. (2013). The thermoeconomics of energy conversions: Elsevier.
Fiorini, P., & Sciubba, E. (2005). Thermoeconomic analysis of a MSF desalination plant. Desalination, 182(1–3), 39-51. doi:http://dx.doi.org/10.1016/j.desal.2005.03.008
Ghaffour, N., Lattemann, S., Missimer, T., Ng, K. C., Sinha, S., & Amy, G. (2014). Renewable energy-driven innovative energy-efficient desalination technologies. Applied Energy, 136, 1155-1165. doi:http://dx.doi.org/10.1016/j.apenergy.2014.03.033
Ghorbani, B., Mehrpooya, M., & Sadeghzadeh, M. (2018). Developing a tri-generation system of power, heating, and freshwater (for an industrial town) by using solar flat plate collectors, multi-stage desalination unit, and Kalina power generation cycle. Energy Conversion and Management, 165, 113-126.
Hisham T. El-Dessouky, H. M. E. (2002). Fundamentals of Salt Water Desalination: ELSEVIER.
Holland, J. H. (1992). Adaptation in natural and artificial systems: MIT Press.
Holmgren, M. (2006). XSTEAM FOR MATLAB. www.x-eng.com.
Huang, Y., Wang, Y. D., Rezvani, S., McIlveen-Wright, D. R., Anderson, M., Mondol, J., . . . Hewitt, N. J. (2013). A techno-economic assessment of biomass fuelled trigeneration system integrated with organic Rankine cycle. Applied Thermal Engineering, 53(2), 325-331. doi:http://dx.doi.org/10.1016/j.applthermaleng.2012.03.041
Kafi, F., Renaudin, V., Alonso, D., & Hornut, J. M. (2004). New MED plate desalination process: thermal performances. Desalination, 166(0), 53-62. doi:http://dx.doi.org/10.1016/j.desal.2004.06.058
Kahraman, N., & Cengel, Y. A. (2005). Exergy analysis of a MSF distillation plant. Energy Conversion and Management, 46(15–16), 2625-2636. doi:http://dx.doi.org/10.1016/j.enconman.2004.11.009
Kamali, R. K., Abbassi, A., Sadough Vanini, S. A., & Saffar Avval, M. (2008). Thermodynamic design and parametric study of MED-TVC. Desalination, 222(1–3), 596-604. doi:http://dx.doi.org/10.1016/j.desal.2007.01.120
Kamali, R. K., & Mohebinia, S. (2008). Experience of design and optimization of multi-effects desalination systems in Iran. Desalination, 222(1–3), 639-645. doi:http://dx.doi.org/10.1016/j.desal.2007.01.182
Kwon, S., Won, W., & Kim, J. (2016). A superstructure model of an isolated power supply system using renewable energy: Development and application to Jeju Island, Korea. Renewable Energy, 97, 177-188. doi:http://dx.doi.org/10.1016/j.renene.2016.05.074
Mabrouk, A.-N. A. (2013). Techno-economic analysis of tube bundle orientation for high capacity brine recycle MSF desalination plants. Desalination, 320(0), 24-32. doi:http://dx.doi.org/10.1016/j.desal.2013.04.014
Malik, M., Dincer, I., & Rosen, M. A. (2015). Development and analysis of a new renewable energy-based multi-generation system. Energy, 79, 90-99.
Mehrpooya, M., Ghorbani, B., Jafari, B., Aghbashlo, M., & Pouriman, M. (2018). Modeling of a single cell micro proton exchange membrane fuel cell by a new hybrid neural network method. Thermal Science and Engineering Progress, 7, 8-19.
Melgar, A., Pérez, J. F., Laget, H., & Horillo, A. (2007). Thermochemical equilibrium modelling of a gasifying process. Energy Conversion and Management, 48(1), 59-67. doi:http://dx.doi.org/10.1016/j.enconman.2006.05.004
Mentis, D., Karalis, G., Zervos, A., Howells, M., Taliotis, C., Bazilian, M., & Rogner, H. (2016). Desalination using renewable energy sources on the arid islands of South Aegean Sea. Energy, 94, 262-272. doi:http://dx.doi.org/10.1016/j.energy.2015.11.003
Modares, H., & NaghibiSistani, M. B. (2011). Solving nonlinear optimal control problems using a hybrid IPSO SQP algorithm. Engineering Applications of Artificial Intelligence, 24, 476–484.
Nafey, A. S., Fath, H. E. S., & Mabrouk, A. A. (2006). Thermo-economic investigation of multi effect evaporation (MEE) and hybrid multi effect evaporation—multi stage flash (MEE-MSF) systems. Desalination, 201(1–3), 241-254. doi:http://dx.doi.org/10.1016/j.desal.2005.09.044
Nisan, S., & Dardour, S. (2007). Economic evaluation of nuclear desalination systems. Desalination, 205(1–3), 231-242. doi:http://dx.doi.org/10.1016/j.desal.2006.05.014
Salimi, M., & Amidpour, M. (2017). Investigating the integration of desalination units into cogeneration systems utilizing R-curve tool. Desalination, 419, 49-59.
Sayyaadi, H., & Ghorbani, G. (2018). Conceptual design and optimization of a small-scale dual power-desalination system based on the Stirling prime-mover. Applied Energy, 223, 457-471.
Sayyaadi, H., & Mehrabipour, R. (2012). Efficiency enhancement of a gas turbine cycle using an optimized tubular recuperative heat exchanger. Energy, 38(1), 362-375.
Shakib, S. E., Amidpour, M., & Aghanajafi, C. (2012). Simulation and optimization of multi effect desalination coupled to a gas turbine plant with HRSG consideration. Desalination, 285(0), 366-376. doi:http://dx.doi.org/10.1016/j.desal.2011.10.028
Shariati Niasar, M., Amidpour, M., Ghorbani, B., Rahimi, M.-J., Mehrpooya, M., & Hamedi, M.-H. (2017). Superstructure of Cogeneration of Power, Heating, Cooling and Liquid Fuels Using Gasification of Feedstock with Primary Material of Coal for Employing in LNG Process. Gas Processing, 5(1), 1-23. doi:10.22108/gpj.2017.100034.1000
Shih, H. (2005). Evaluating the technologies of thermal desalination using low-grade heat. Desalination, 182(1–3), 461-469. doi:http://dx.doi.org/10.1016/j.desal.2005.02.038
Smith, R. (2005). Chemical process: design and integration: John Wiley & Sons.
Sonntag, R. E., Borgnakke, C., Van Wylen, G. J., & Van Wyk, S. (1998). Fundamentals of thermodynamics (Vol. 6): Wiley New York.
Thermoflow (GTpro module). Thermoflow Inc.
Yassin, L., Lettieri, P., Simons, S. J. R., & Germanà, A. (2009). Techno-economic performance of energy-from-waste fluidized bed combustion and gasification processes in the UK context. Chemical Engineering Journal, 146(3), 315-327. doi:https://doi.org/10.1016/j.cej.2008.06.014 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

آمار تعداد مشاهده مقاله: 958 تعداد دریافت فایل اصل مقاله: 334 |