Operating Energy Needed for Desalination Systems in Cogeneration Plants

Abstract

This study investigated the energy requirement for running desalination units coupled to cogeneration plants. Various cogeneration systems were explored using power- and heat-allocated approaches. The specific work and heat necessary for operating different desalination systems were determined. The investigation revealed that the specific work and heat remain consistent regardless of the desalination daily capacity. It was observed that the energy demand for operating a desalination system mainly relies on power plant efficiency. The investigation revealed that the energy demand for a plain multi-effect desalination system was lower than that for multi-effect desalination with thermal vapor compression. Additionally, the energy requirement for a multi-effect desalination system with preheaters was lower than that for plain multi-effect desalination. Comparisons also indicated that the energy demand of multi-stage flash exceeds that of different multi-effect desalination systems. Based on the primary thermal energy input, a universal performance ratio was used to evaluate the desalination unit performance. Furthermore, a new correlation was proposed to predict the universal performance ratio.

1. Introduction

Freshwater demand is a critical concern in Saudi Arabia and worldwide. Desalination is a reliable solution to tackle water shortage problems. Based on the Shared Socioeconomic Pathway (SSP2) study, Gao et al. [1] concluded that, by 2050, desalination will become increasingly feasible for developing nations. Desalination is an energy-consuming process. Consequently, the adaption of cogeneration, which combines power generation and desalination, is used to maximize the system’s energy utilization efficiency. The US Department of Energy recommended using renewable energy for running cogeneration plants to minimize the operational expenses of desalination water production [2].

The primary factor in selecting the appropriate cogeneration systems lies in the economics of water production. Evaluating the performance of cogeneration plants, where electric and thermal energy are involved, is a complex task and is not easy to perform. Detailed reviews of thermo-economic models assessing the cost of energy required for desalination processes are presented in [3,4]. Thermo-economic analyses of desalination systems were carried out by [5,6,7,8] based on the economic model established by Kavvadias and Khamis [9]. The study of Khan et al. [5] indicated differing cost trends for the examined desalination technologies. Meanwhile, Haya et al. [6] suggested that using nuclear energy in cogeneration was a promising technique to reduce water production expenses. Moreover, Sadeghi et al. [7] proposed a hybrid reverse osmosis multi-effect desalination (RO-MED) system with an optimal ratio of 0.7 to reduce water production costs.

Major desalination processes include thermal and electric-driven systems. Comparisons of the energy required to operate these systems are not a straightforward task [10,11,12,13]. According to Altman et al. [10], comparing the thermal energy requirements of multi-effect desalination (MED) and multi-stage flash (MSF) systems with the electric energy needs of reverse osmosis (RO) systems could not lead to definitive conclusions. As reported by [10,14], there is no standard technique to estimate the cost of freshwater production in cogeneration plants. Consequently, distributing the energy cost between power and water production remains a subject of debate. ElNashar [12] pointed out that existing methods used for cost allocation in cogeneration plants presented significant challenges.

The feasibility of the cogeneration systems can be judged based on the produced freshwater levelized cost. The type and the quantity of energy required to operate desalination systems represent important factors. Various energy demand approaches were explored, including exergy and energy methods [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25]. An exergy approach was presented by [10,16,21,22,26]. Altmann et al. [10] developed a theoretical exergy-based approach to determine the energy consumption of 48 different configurations of cogeneration plants. In this investigation, the heat and power requirements of these systems were traced back to their primary form. Shahzad et al. [21,22] and Ng et al. [26] introduced an exergy-based analysis to determine the quantity of fuel needed by the desalination units in a cogeneration plant. Shahzad et al. [21,22] and Ng et al. [26] assessed different desalination methods using a universal performance ratio (UPR) based on the type of primary energy input.

Two main approaches can be employed to determine the energy requirements of desalination systems: the power-allocated method (PAM) [11,15,16,17,18,23,24,25] and the heat-allocated method (HAM) [19,20]. In the PAM approach, steam extraction from the turbine to heat the thermal desalination system decreases the cogeneration power plant’s capacity to generate electricity. The reduction in the cogeneration electric power output is considered as the energy demand to run the desalination system. The HAM method considers the additional heat required to substitute the steam thermal energy extracted to heat the desalination system as the energy requirement of the desalination system. This additional heat is needed to maintain the stand-alone power plant’s full load output.

Zeitoun et al. [27] developed thermo-economic models for cogeneration power and desalination systems. In these models, the levelized cost of water was estimated for nine configurations where three types of desalination systems were connected with three types of power plants. In this study, the energy consumption of the desalination systems was estimated using the PAM and HAM procedures. The main key findings of this investigation were as follows:

• Applying the HAM method for estimating the levelized cost of water narrowed the gap between the cost of water produced by the MED and the seawater reverse osmosis (SWRO) systems.
• The levelized cost of water for the simple MED was lower than that for multi-effect distillation with thermal vapor compression (MED-TVC).
• The profit margin for cogeneration combined power plants powered by natural gas was higher than that for nuclear power plants.

The previous investigation of Zeitoun et al. [27] raised many important questions, including the impact of using highly efficient power plants, e.g., supercritical steam power plants, on the water levelized cost, if the simple MED competes with MED-TVC, and the impact of using preheaters with the simple MED.

The current study focused on the energy consumption of thermal desalination systems integrated with steam, nuclear, combined, and supercritical pressure power plants. The desalination systems investigated included plain MED, MED with preheaters, MED-TVC, and once-through MSF systems. The energy consumption of these desalination systems was compared to the SWRO energy requirement to explore the reason for RO spread adoption. According to published statistics, most of the desalination plants newly installed worldwide are RO-based [8].

The procedure of the current study involved the following steps:

-

Analyze the performance of a stand-alone power plant to determine its full power output and overall efficiency for a given amount of heat added;

-

Analyze the performance of a power plant connected to a desalination system to determine the loss of its power generation for the same amount of heat added to determine the loss in power.

-

Analyze the performance of a power plant connected to a desalination system to determine the additional heat required to maintain the full power output.

-

Use the above results to determine the energy requirement to run desalination systems based on the PAM and HAM methods.

2. Stand-Alone Power Plant Simulation

In the current investigation, supercritical power plants (ScPPs) connected to plain MED, multi-effect with preheaters (MED-PH), MED-TVC and once-through MSF desalination systems were examined. The configurations of these cogeneration plants are shown in Figure 1, Figure 2, Figure 3 and Figure 4. The configurations presented by Zeitoun et al. [27] were revisited to determine the energy requirements for operating MED, MED-TVC, and MSF systems connected to a steam power plant (SPP), nuclear power plant (NPP), and combined power plant (CPP). This study extended the cases analyzed by [27] by including the MED-PH configuration. The figures illustrating these configurations are shown in Appendix A. The complete list of cogeneration configurations investigated is listed in

Table 1. Investigated configurations of cogeneration power plants/desalination systems.

Power Plant	Desalination Systems
	Plain MED	MED with Preheaters	MED-TVC	OT MSF
Steam power plant	SPP-MED	SPP-MED-PH	SPP-MED-TVC	SPP-MSF
Nuclear power plant	NPP-MED	NPP-MED-PH	NPP-MED-TVC	NPP-MSF
Combined power plant	CPP-MED	CPP-MED-PH	CPP-MED-TVC	CPP-MSF
Supercritical power plant	ScPP-MED	ScPP-MED-PH	ScPP-MED-TVC	ScPP-MSF

.

A plain MED unit connected to a ScPP is depicted in Figure 1. The supercritical pressure boiler supplies steam to a turbine of three stages. The steam undergoes reheating after expansion in the high-pressure turbine, as shown in Figure 1. The power plant is equipped with nine feed-water heaters; one of the open type and eight of the closed type, as shown in Figure 1. The pressure distribution along the feed-water heaters was optimized as reported in [28]. The parameters used to simulate the ScPP are detailed in

Table 2. Stand-alone ScPP simulation data.

ScPP Input Data			Components’ Efficiencies
	Pres. (bar)	Temp. (°C)	ηt	90%
Boiler outlet	330	610	ηb	95%
Condenser	0.086	43	ηp	90%
Seawater		33	ηm	95%
Reheater outlet	45	630
CFWH1	68.38	284.3	ScPP Output Data
CFWH2	45	257.5	ηov	44.19%
CFWH3	28.29	230.7	HRov	8146	kJ/kWh
CFWH4	16.83	203.9	Wfull	1325.78	MW
OFWH	9.356	177
FWH5	4.787	150.2
CFWH6	2.21	123.4
CFWH7	0.8969	96.6
CFWH8	0.3092	69.8

. Thermal simulation of the stand-alone ScPP was performed using EES software V11.725 [29]. The stand-alone simulation was conducted to obtain the overall characteristics of the power plant under full load conditions.

For the investigated power plants, the following input data were maintained:

• The heat added to operate the different examined power plants was maintained at Q_H = 3000 MW.
• The rise in the temperature of the condenser cooling water was assumed to be 5 °C, and the terminal temperature difference at the condenser exit was assumed to be 5 °C.
• The efficiencies of the components of ScPP are listed in Table 2.

A thermodynamic model using the mass and energy conservation equations of the ScPP components was developed. The overall balances of mass and energy for the power plant were checked for the solved conditions. At the exit of the boiler of the ScPP, the pressure and temperature were maintained at 330 bar and 610 °C, respectively, see Table 2. The pressures along the components of the ScPP are listed in Table 2. The steam was reheated to 630 °C after the expansion in the high-pressure turbine. The simulation results of the stand-alone ScPP plant are also listed in Table 2, where the overall efficiency is 44.19%, and full output is 1325.78 MW.

3. Cogeneration Plant Simulation

In the cogeneration plant simulation, thermodynamic models of plain MED, MED-PH, MED-TVC, and MSF systems [30,31] were integrated into an ScPP model for the configurations shown in Figure 1, Figure 2, Figure 3 and Figure 4. The validation of the desalination models was discussed in [27]. In addition, a model of MED-PH was integrated into the models of SPP, CPP, and NPP plants presented in [27], and their configurations are shown in Appendix A. The main specifications and parameters of the investigated desalination systems are listed in

Table 3. Desalination system characteristics.

Cogeneration Plant	N	Heating Steam		Primary Steam for TVC		Simulation Obtained Data
		P, Bar	Tsat, °C	P, Bar	Tsat, °C	TBT, °C	GORth	QHAM kWh/m3	WPAM kWh/m3	UPR	ηD	ηov	UPRPr	$$\begin{gathered} \frac{∆\mathit{U}\mathit{P}\mathit{R}}{\mathit{U}\mathit{P}\mathit{R}} \\ \% \end{gathered}$$
SPP-MED	8	0.41	76.5			70.3	5.98	24.49	9.64	27.21	0.085	0.394	27.69	1.8
NPP-MED		0.416	76.9			70.6	5.97	27.83	9.90	23.95	0.087	0.356	24.32	1.5
CPP-MED		0.236	63.7			59.1	6.27	13.69	7.04	48.66	0.064	0.514	50.69	4.2
ScPP-MED		0.31	69.8			66.5	6.19	17.25	7.62	38.63	0.068	0.442	39.99	3.5
SPP-MED-PH	8	0.41	76.5			70.3	7.55	19.61	7.75	33.97	0.085	0.394	34.97	2.9
NPP-MED-PH		0.416	76.9			70.6	7.55	22.15	7.88	30.08	0.087	0.356	30.77	2.3
CPP-MED-PH		0.236	63.7			59.1	7.58	11.49	5.91	57.99	0.064	0.514	61.27	5.7
ScPP-MED-PH		0.31	69.8			66.5	7.57	14.2	6.28	46.92	0.068	0.442	48.95	4.3
SPP-MED-TVC	10 NTVC = 7	0.312	70	3.99	143.5	67.3	12.7	28.49	11.25	23.39	0.237	0.394	23.93	2.3
NPP-MED-TVC		0.312	70	4.1	144.6	67.3	13.13	31.49	11.21	21.16	0.241	0.356	22.41	5.9
CPP-MED-TVC		0.312	70	6.8	164	67.3	12.9	30.38	15.62	21.93	0.301	0.514	25.62	16.8
ScPP-MED-TVC		0.312	70	4.79	150.2	67.3	12.68	27.12	11.98	24.57	0.253	0.442	25.17	2.4
SPP-MSF	40	3.99	143.5			120	10.52	34.18	13.46	19.55	0.237	0.394	17.51	−10.4
NPP-MSF		4.1	144.6			120	10.52	38.13	13.57	17.53	0.241	0.356	15.56	−11.2
CPP-MSF		6.8	164			120	10.52	38.6	19.85	17.31	0.301	0.514	17.97	3.8
ScPP-MSF		4.79	150.2			120	10.52	32.62	14.41	20.49	0.253	0.442	18.34	−10.5

.

The investigated plain MED and MED-PH systems consisted of eight effects, as shown in Figure 1 and Figure 2. The steam was withdrawn at the ninth feed water heater to heat the MED and MED-PH systems, as shown in Figure 1 and Figure 2. The thermal simulation of the ScPP coupled to a MED or a MED-PH was conducted for different desalination daily capacities using EES software [29]. The energy requirement to run the MED and MED-PH was determined according to the PAM and HAM approaches discussed before. Figure 3 shows the investigated MED-TVC system. The MED-TVC consisted of 10 effects where the TVC extracted vapor from the seventh effect. The heating vapor was fed to the first effect at 70 °C. The examined once-through MSF, shown in Figure 4, consisted of 40 stages where the top brine temperature was maintained at 120 °C and the TTD was fixed at 7 °C at the MSF exit [31]. The thermal gain output ratio of the simulated systems is included in Table 3 The thermal gain output ratio GOR_th, based on the energy ratio, was estimated as follows:

$${GOR}_{th}=\frac{m_d h_{fgN}}{m_{hs}{∆h}_{hs}}$$

(1)

where m_d is the desalination capacity in kg/s, m_hs is the steam withdrawn to heat the thermal desalination system, h_fgN is the latent heat at the last effect or the last stage, and ∆h_hs is the drop in the enthalpy of the extracted steam. The data obtained for the cases investigated by Zeitoun et al. [27] are also included in Table 3.

3.1. Desalination Energy Requirement According to PAM

In the PAM approach, the decrease in the power plant net output was considered as the energy necessary to run the desalination system. The fuel energy (Q_H) supplied to the boiler was kept at a rate of 3000 MW. The power requirement to run the desalination system W_D included the power lost due to the steam extraction and the power required to run the pumps of the desalination system [27]:

$$W_D=∆W+W_{pump}$$

(2)

where ΔW is the drop in the cogeneration plant output power [27];

$$W_D=W_{full}-W_p+W_{pump}$$

(3)

W_full represents the full power of the stand-alone plant, and W_p represents the power developed by the cogeneration plant. W_pump is the power required to circulate the seawater and to extract the produced fresh water and concentrate brine from the desalination unit; see Appendix B for details.

The specific work W_PAM (kWh/m³), which is the work needed to yield one cubic meter of fresh water, can be estimated by dividing the daily energy consumption in kWh by the daily desalination capacity:

$$W_{PAM}=24W_D/m_D$$

(4)

The equivalent specific heat Q_PAM (kWh/m³), which is the heat needed to produce one cubic meter of fresh water, was determined from the following:

$$Q_{PAM}=W_{PAM}/{\mathit{\eta }}_{ov}$$

(5)

where η_ov is the stand-alone power plant overall efficiency at the full load. The specific work and heat required to run the different investigated desalination systems, using the PAM method, are shown in Figure 5 and Figure 6, respectively. W_PAM data for the investigated configuration are also listed in Table 3.

3.2. Desalination Energy Requirement According to HAM

In the HAM method, the ScPP output connected to a desalination unit was maintained at the designed full load W_full, (Table 2). Consequently, the heat input to the cogeneration plant was increased by ∆Q to maintain the full load. This heat increase is needed to compensate for the thermal energy of the steam pulled out to heat the desalination unit. It was assumed that the boiler of the ScPP could produce additional steam to cover the heating steam needs of the connected desalination system. The heat needed to run the desalination system includes the heat increase and the heat corresponding to the pumping work of the desalination system:

$${\mathrm{Q}}_D=∆Q+W_{pump}/{\mathit{\eta }}_{ov}$$

(6)

The specific heat Q_HAM (kWh/m³), which is the heat needed to produce one cubic meter of fresh water, was estimated by dividing the daily heat consumption in kWh by the daily desalination capacity:

$$Q_{HAM}=24 Q_D/m_D$$

(7)

The specific equivalent work W_HAM (kW h/m³), which is the work needed to produce one cubic meter of fresh water, can be estimated as follows:

$$W_{HAM}=Q_{HAM} {\mathit{\eta }}_{ov}$$

(8)

The specific heat and work required to run the different investigated desalination systems using the HAM approach are also presented in Figure 5 and Figure 6, respectively. Q_HAM data for the investigated cases are also listed in Table 3.

3.3. Results for Energy Requirements of Desalination Systems

The data shown in Figure 5 and Figure 6 and listed in Table 3 indicate that the PAM and HAM approaches lead to the same results, and the specific heat and work needed to operate the examined desalination systems are independent of the systems’ daily desalination capacities. The obtained data revealed that the energy needed to run the MED-PH system is the lowest and the energy needed to operate the MSF system is the highest. Extracting steam with high pressure, as in MSF and MED-TVC, reduces the power output of the low-pressure steam turbines, leading to a significant increase in the energy needed to operate the MSF and MED-TVC systems. The energy needed to run the desalination system does not mainly depend on the GOR_th of the desalination system. The MED-TVC system, of high GOR_th, needs more energy than the simple MED system of low GOR_th. The GOR_th of the desalination system does not separately indicate the efficiency of the desalination unit. The overall efficiency of the power plant represents a main parameter in determining the efficiency of the desalination plant.

Figure 7 and Figure 8 show the specific work and heat needed to operate different desalination systems. The data for the needed work of Figure 7 falls near the range reported by Al-Karaghouli and Kazmerski [32]. They reported 14.45–21.35 kWh/m³ for MED-TVC and 19.6–27.25 kWh/m³ for MSF. These values were estimated based on 30% overall power plant efficiency. Figure 8 shows that the heat needed to run the desalination system decreases as the overall efficiency of the plant increases. However, the trend is different for MED-TVC and MSF connected to a CPP. There are two reasons for this problem: the high pressure of the extracted steam and the limited amount of steam in the CPP. These two factors dramatically affect the performance of the CPP connected to MED-TVC or MSF.

4. Universal Performance Ratio

The universal performance ratio (UPR) presented in [21,22] can be used to indicate the performance of a thermal desalination unit connected to a power plant. The UPR represents the ratio between the heat required to evaporate the desalinated water production rate and the thermal primary form of energy required to run the desalination system. UPR was reduced to the following form in the current investigation:

$$UPR=1000 h_{fgN}/3600 Q_{HAM}$$

(9)

Q_HAM can be replaced by Q_PAM. The obtained UPR data are included in Table 3. Figure 9 shows the UPR for different desalination systems. Figure 9 indicates that the simple MED without or with preheaters, MED or MED-PH, performs better than MED-TVC, which was also concluded by [10]. Removing the TVC will also reduce the capital investment in the desalination system. The UPR of MSF is the poorest compared to simple MED, MED-PH, and MED-TVC. The reason for the low UPR of MSF and MED-TVC is the high reduction in the power plant output due to the high pressure of the extracted steam. Figure 9 indicates that the UPR mainly depends on the power plant’s overall efficiency. For simple MED and MED-PH, the UPR significantly increases from NPP to CPP, i.e., from a low-efficiency to high-efficiency plant. The UPR of the desalination systems connected to the CPP is the best due to the CPP’s high efficiency. However, when a CPP is connected to MED-TVC or MSF, the UPR decreases. As discussed earlier, there are two reasons for this problem: the high pressure of the extracted steam required to heat the desalination system and the limited quantity of steam in the CPP.

As reported by [21,22], the minimum work for salt separation is W_min = 0.78 kWh/m³. The theoretical limit of a desalination system connected to the CPP can be estimated as follows:

$${UPR}_{TL}=1000 h_{fgN} {\mathit{\eta }}_{ov}/3600 W_{min}$$

(10)

where η_ov is the overall efficiency of the CPP at full load. As mentioned before, the energy required to run the desalination system was estimated by [10,21,22] using the exergy approaches. Figure 10 illustrates a comparison between the theoretical limit, the current data, and the data of [10,21,22]. The energy required to run SWRO, in the current investigation, is based on a specific work consumption of 3.5 kWh/m³, as obtained from [33]. The data in Figure 10 indicate the large gap between the theoretical limit and the actual cases.

Investigating the UPR data in Figure 9 or Table 3 indicates that the UPR depends on the GOR_th of the desalination system and the overall efficiency of the stand-alone power plant. It is known that the performance of the cogeneration system deteriorates as the pressure of the extracted steam increases. Consequently, a new parameter was introduced to design a correlation for UPR prediction. This parameter represents the efficiency of the Rankine cycle related to the desalination process. This efficiency η_D represents the ratio between the lost expansion turbine work and the thermal energy extracted from the heating steam,

$${\mathit{\eta }}_D=\frac{h_a-h_b}{h_a-h_{fr}}$$

(11)

where h_a is the enthalpy of the extracted steam to heat the desalination system, h_b is the enthalpy of the low-pressure turbine exit, and h_fr is the enthalpy of the saturated liquid returned to the power plant cycle. The following correlation was introduced to predict the UPR:

$$UPR={GOR}_{th}{\mathit{\eta }}_{ov}/{\mathit{\eta }}_D$$

(12)

The predicted UPRs of this correlation are listed in Table 3. The differences between predicted UPR_pr and system UPR are also listed in Table 3. The comparison between the predicted UPR_pr and UPR data is shown in Figure 11. For MED and MED-PH systems, the predictions of the new correlation are reasonably good and fall within +1.5 and 5.7% of the system UPR. For MED-TVC and MSF, the predictions of the correlation are reasonable and fall within −10% and +16.8% of system UPR. The differences between predicted and actual UPR can be attributed to the extracted steam conditions, which are not completely dry saturated. The extracted steam was either superheated or wet vapor.

5. Conclusions

The energy required to operate desalination systems coupled to different power plants was investigated using the power- and heat-allocated methods. Both procedures led to the same results. The obtained data revealed the following:

• The specific heat and work needed to operate the examined desalination systems are independent of daily system desalination capacity.
• The MED-PH system needs the lowest amount of energy to operate.
• The MSF system needs the highest amount of energy to operate.
• The MED-TVC system needs more energy than the simple MED system.
• Extracting steam with high pressure, as in MSF and MED-TVC, reduces the output power of the low-pressure turbines, leading to a significant increase in the energy required to run the desalination systems.
• The GOR_th of the desalination system, in cogeneration plants, does not independently indicate the effectiveness of the desalination plant. The overall efficiency of the power plant is a main parameter in determining the effectiveness of the desalination plant.

The universal performance ratio, based on the thermal primary form of energy input, was used to evaluate the systems’ performance. A new correlation was introduced to predict the universal performance ratio. The predicted value of the new correlation was good for the plain MED and MED-PH systems and reasonable for the MED-TVC and MSF systems. For the MED and MED-PH systems, the predictions of the UPR correlation fell within +1.5 and 5.7% of the actual UPR. For the MED-TVC and MSF, the predictions of the UPR correlation fell within −10% and +16.8% of system UPR.