Evaluation of the Specific Energy Consumption of Sea Water Reverse Osmosis Integrated with Membrane Distillation and Pressure–Retarded Osmosis Processes with Theoretical Models

In this study, theoretical models for specific energy consumption (SEC) were established for water recovery in different integrated processes, such as RO-PRO, RO-MD and RO-MD-PRO. Our models can evaluate SEC under different water recovery conditions and for various proportions of supplied waste heat. Simulation results showed that SEC in RO increases with the water recovery rate when the rate is greater than 30%. For the RO-PRO process, the SEC also increases with the water recovery rate when the rate is higher than 38%, but an opposite trend can be observed at lower water recovery rates. If sufficient waste heat is available as the heat source for MD, the integration of MD with the RO or RO-PRO process can significantly reduce SEC. If the total water recovery rate is 50% and MD accounts for 10% of the recovery when sufficient waste heat is available, the SEC values of RO, RO-PRO, RO-MD and RO-MD-PRO are found to be 2.28, 1.47, 1.75 and 0.67 kWh/m3, respectively. These critical analyses provide a road map for the future development of process integration for desalination.


Introduction
Climate change and water scarcity are the two pervasive problems posing serious threats to people around the world [1]. To solve these problems, seawater desalination technology for potable water has been rapidly developed. Reverse osmosis (RO) is currently the most common membrane process for desalination since high-performance RO membranes and modules have been well established. However, this process requires high electric energy consumption and the discharge of concentrated brine solution, which have confined the application of RO systems for years [2][3][4][5][6][7]. Thus, developing methods to reduce RO energy consumption and brine disposal is crucial for process efficiency and sustainability. One of the techniques that has received extensive attention is process integration, or the combination of RO with pressure-retarded osmosis (PRO) and/or membrane distillation (MD). This integrated approach can reduce the energy requirement for water recovery and enhance the water recovery rate [8].
PRO is an energy-harvesting process that converts the osmotic pressure of a saline solution to hydraulic pressure. In a PRO system, water from a low-salinity solution

Materials and Methods
In this study, Equation (1) is used to indicate the energy consumed per unit of water production [2,3]. SEC = power consumption (kW) water production rate ( m 3 h ) This study integrated RO with other unit operations, such as ERD, MD, and PRO, with the aim of analyzing SEC variations in different integrated processes.

SEC of the RO Process
A schematic diagram of the stand-alone RO process is shown in Figure 1. Since the suspended particles, microbes, and organic and inorganic matter in seawater can block the RO membrane and lead to fouling, the raw seawater feed solution is often pretreated with UF/MF prior to RO. The SEC for pretreatment per unit of RO permeate (SEC Pre ) is described in Equation (2). SEC Pre = SEC pre,sw Q SW Q P,RO = SEC pre,sw Y r (2) where SEC Pre,sw is the pretreatment energy consumption per unit seawater feed solution, Q SW is the volumetric flow rate of the seawater feed solution (m 3 /h), Q P,RO is the RO permeate flow rate (m 3 /h), and Y r = Q P,RO Q SW is the RO water recovery rate.

Materials and Methods
In this study, Equation (1) is used to indicate the energy consumed per unit of water production [2,3].

SEC = power consumption(kW)
water production rate ( m 3 h ) (1) This study integrated RO with other unit operations, such as ERD, MD, and PRO, with the aim of analyzing SEC variations in different integrated processes.

SEC of the RO Process
A schematic diagram of the stand-alone RO process is shown in Figure 1. Since the suspended particles, microbes, and organic and inorganic matter in seawater can block the RO membrane and lead to fouling, the raw seawater feed solution is often pretreated with UF/MF prior to RO. The SEC for pretreatment per unit of RO permeate (SECPre) is described in Equation (2). SEC Pre = SEC pre,sw Q SW Q P,RO = SEC pre,sw Y r (2) where SECPre,sw is the pretreatment energy consumption per unit seawater feed solution, QSW is the volumetric flow rate of the seawater feed solution (m 3   Glueckstern and Priel suggested that when ultrafiltration is used prior to RO, SECPre,sw is 0.095 kWh/m 3 [29]. Pretreated seawater is pressurized via a high-pressure pump (HP) to the RO operating pressure (PRO) and then enters the RO unit. The minimum operating pressure required for RO is the osmotic pressure of the brine. In this study, the RO operating condition is set equal to the minimum operating pressure at each water recovery rate (i.e., osmotic pressure of brine in RO) for determining the SEC. The SEC of the stand-alone RO process (SECRO) was evaluated by Equation (3) [2,3].
where πsw is the osmotic pressure of the feed solution (bar), Rt is the salt rejection rate of the membrane, and ηp is the HP efficiency. The osmotic pressure can be calculated by the van 't Hoff equation as below [30]: where π is the osmotic pressure (bar), i is the van 't Hoff factor (i = 1.9 when the solute is NaCl [30]), C is the salt concentration (M), R is the ideal gas constant, and T is the temperature (K). In this study, the temperature is set to 298 K, unless a different value is specified. The RO-ERD process is presented in Figure 2. ERD is utilized to recover the remaining mechanical energy in the brine solution; herein, a PX was selected as the ERD. After seawater pretreatment, part of the feed solution enters the RO unit via HP, and the other Glueckstern and Priel suggested that when ultrafiltration is used prior to RO, SEC Pre,sw is 0.095 kWh/m 3 [29]. Pretreated seawater is pressurized via a high-pressure pump (HP) to the RO operating pressure (P RO ) and then enters the RO unit. The minimum operating pressure required for RO is the osmotic pressure of the brine. In this study, the RO operating condition is set equal to the minimum operating pressure at each water recovery rate (i.e., osmotic pressure of brine in RO) for determining the SEC. The SEC of the stand-alone RO process (SEC RO ) was evaluated by Equation (3) [2,3].
where π sw is the osmotic pressure of the feed solution (bar), R t is the salt rejection rate of the membrane, and η p is the HP efficiency. The osmotic pressure can be calculated by the van 't Hoff equation as below [30]: where π is the osmotic pressure (bar), i is the van 't Hoff factor (i = 1.9 when the solute is NaCl [30]), C is the salt concentration (M), R is the ideal gas constant, and T is the temperature (K). In this study, the temperature is set to 298 K, unless a different value is specified. The RO-ERD process is presented in Figure 2. ERD is utilized to recover the remaining mechanical energy in the brine solution; herein, a PX was selected as the ERD. After seawater pretreatment, part of the feed solution enters the RO unit via HP, and the other part enters the RO unit via the PX and booster pump (BP). The brine solution is then discharged after the available brine energy is transferred via the PX at atmospheric pressure. The SEC of the RO-ERD process (SEC RO-ERD ) can be calculated as in Equation (5) [2,3]: where η PX is the PX efficiency.
where ηPX is the PX efficiency.  Figure 3 shows the flow chart of the RO-PRO process. The brine solution, which served as the DS for the PRO unit, is discharged from the RO system and releases pressure via the PX1 to PRO operating pressure (PPRO). Previous studies often used low-saline water as the FS for PRO; however, in practice, urban wastewater should be considered as the PRO FS [27]. In this study, a 0.01-M NaCl solution is used to simulate the FS of the PRO unit. Considering the environmental impact after a large amount of brine discharge, the RO brine is assumed to be diluted to the original seawater concentration after passing through the PRO unit. The SEC of RO-PRO (SECRO-PRO) can be calculated as in Equation (6

SEC of the RO-PRO Process
where W pump,RO-PRO and W pump,RO are the energy consumption levels of the HP in RO-PRO and stand-alone RO processes, respectively. W pump,RO can be calculated in Equation (7).   Figure 3 shows the flow chart of the RO-PRO process. The brine solution, which served as the DS for the PRO unit, is discharged from the RO system and releases pressure via the PX1 to PRO operating pressure (P PRO ). Previous studies often used low-saline water as the FS for PRO; however, in practice, urban wastewater should be considered as the PRO FS [27]. In this study, a 0.01-M NaCl solution is used to simulate the FS of the PRO unit. Considering the environmental impact after a large amount of brine discharge, the RO brine is assumed to be diluted to the original seawater concentration after passing through the PRO unit. The SEC of RO-PRO (SEC RO-PRO ) can be calculated as in Equation (6)

SEC of the RO-PRO Process
where ηPX is the PX efficiency.  Figure 3 shows the flow chart of the RO-PRO process. The brine solution, which served as the DS for the PRO unit, is discharged from the RO system and releases pressure via the PX1 to PRO operating pressure (PPRO). Previous studies often used low-saline water as the FS for PRO; however, in practice, urban wastewater should be considered as the PRO FS [27]. In this study, a 0.01-M NaCl solution is used to simulate the FS of the PRO unit. Considering the environmental impact after a large amount of brine discharge, the RO brine is assumed to be diluted to the original seawater concentration after passing through the PRO unit. The SEC of RO-PRO (SECRO-PRO) can be calculated as in Equation (6).

SEC of the RO-PRO Process
where W pump,RO-PRO and W pump,RO are the energy consumption levels of the HP in RO-PRO and stand-alone RO processes, respectively. W pump,RO can be calculated in Equation (7).  Previous studies indicated that the maximum power density can be obtained when the operating pressure equals half of the transmembrane osmotic pressure [31]. Thus, P PRO is set accordingly in Equation (8).
where π RO−PRO,DS and π PRO,FS are the osmotic pressures of the DS and FS, respectively. The recovery of the mechanical energy of the RO brine after depressurization from P RO to P PRO via PX1 ( . W pump,PX1 ) can be evaluated by Equation (9).
where Q B,RO is the RO brine flow rate. The recovery of mechanical energy at the PRO DS outlet after depressurization from P PRO to atmospheric pressure via PX2 ( . W pump,PX2 ) can be evaluated by Equation (10).
where Q PRO,O and Q P,PRO are the flow rate at the PRO DS outlet and permeation rate through the PRO membrane, respectively. After combining Equation (6) and (10), the SEC of the RO-PRO process can be presented as shown in Equation (11).

SEC of the RO-MD Process
.
where Q P,MD is the MD permeate flow rate and Q B,MD is the brine flow rate discharged from MD. Assuming there is sufficient waste heat for the MD unit, SEC RO-MD described in Equation (12) can be converted to the form shown in Equation (14).
W pump,PX1 is the mechanical energy recovered by PX1, W pump,PX2 is the energy recovered by PX2, PRO is the RO operating pressure, and Qsw is the volumetric flow rate at the seawater inlet. Previous studies indicated that the maximum power density can be obtained when the operating pressure equals half of the transmembrane osmotic pressure [31]. Thus, PPRO is set accordingly in Equation (8).
where π RO-PRO,DS and π PRO,FS are the osmotic pressures of the DS and FS, respectively. The recovery of the mechanical energy of the RO brine after depressurization from PRO to PPRO via PX1 (W pump,PX1 ) can be evaluated by Equation (9).
where QB,RO is the RO brine flow rate. The recovery of mechanical energy at the PRO DS outlet after depressurization from PPRO to atmospheric pressure via PX2 (W pump,PX2 ) can be evaluated by Equation (10).
where QPRO,O and QP,PRO are the flow rate at the PRO DS outlet and permeation rate through the PRO membrane, respectively. After combining Equation (6) and Equation (10), the SEC of the RO-PRO process can be presented as shown in Equation (11). Figure 4 illustrates the RO-MD process. RO brine is first depressurized to PMD via PX3 and then transported to the MD unit. The SEC of RO-MD (SEC RO-MD ) and the mechanical energy recovered via PX3 (W pump,PX3 ) can be presented by Equations (12) and (13), respectively.

SEC of the RO-MD process
W pump,PX3 =η PX (P RO -P MD )Q B,RO (13) where QP,MD is the MD permeate flow rate and QB,MD is the brine flow rate discharged from MD. Assuming there is sufficient waste heat for the MD unit, SEC RO-MD described in Equation (12) can be converted to the form shown in Equation (14).
Membranes 2022, 12, 432 is the MD water recovery rate. If the waste heat supply for MD is not sufficient and a supplementary heat source is then required, the MD thermal energy consumption (SEC Thermal-MD ) and total SEC of RO-MD (SEC Thermal RO-MD ) can be expressed as shown in Equations (15) and (16).
SEC Thermal RO-MD = SEC RO-MD +SEC Thermal-MD (16) where X is the proportion of waste heat supplied to the total thermal energy consumption of MD and EE is the energy efficiency of MD, which is defined as the percentage of the thermal energy associated with liquid evaporation [13]. ρ P,MD is the density of water, and ∆H vap is the enthalpy of vaporization.

SEC of the RO-MD-PRO Process
The power density of PRO can be enhanced by increasing the DS concentration; therefore, utilizing the high-concentration RO-MD brine discharge as the PRO DS can improve power production. Figure 5 is the schematic diagram of the integrated RO-MD-PRO process. MD brine is pressurized in the PRO unit as the DS. It is assumed that PRO DS will be diluted to the initial seawater concentration and that the mechanical energy will be transferred to the RO feed solution via PX5. The SEC of RO-MD-PRO (SEC RO-MD-PRO ) and the recovery mechanical energy from PX4 ( . W pump,PX4−recover ) can be expressed as shown in Equations (17) and (18).
where . W pump,PX4-reuse is the mechanical energy transferred from the RO brine to PRO DS via PX4, and . W pump,PX5 is the mechanical energy recovered from the PRO DS outlet as the DS depressurized from PPRO to atmospheric pressure via PX5. By combining Equations (17)- (20) and assuming that there is sufficient waste heat to meet the thermal energy requirement of MD (X = 1), the SEC of the RO-MD-PRO process can be obtained.  By combining Equations (17)- (20) and assuming that there is sufficient waste heat to meet the thermal energy requirement of MD (X = 1), the SEC of the RO-MD-PRO process can be obtained.

The Fractional Energy Savings
Based on the RO-ERD process, the fractional energy savings (FES) in water production based on hybrid processes can be defined as shown in Equation (22).
When SEC integrated process is zero, the FES equals 100%. Under these circumstances, energy consumption and energy generation are equal, i.e., no further energy supply is required for this process. However, when the FES is greater than 100%, the energy generation is larger than the energy consumption, and SEC integrated process for this process is negative, and vice versa.

Results and Discussion
Herein, the SEC of RO, RO-PRO, RO-MD, and RO-MD-PRO were analyzed under different operating conditions. These calculated SECs were further compared with data from previous research articles. Table 1 lists the parameters used in stand-alone RO and/or RO-ERD process simulations. Figure 6 shows the variations in SEC RO and SEC RO-ERD at different water recovery rates. The SEC of the pretreatment unit, SEC pre , is also shown in Figure 6. As the water recovery rate decreases, the amount of seawater to be pretreated increases, which leads to an increase in SEC pre . SEC RO is much higher than SEC pre ; therefore, the influence of SEC pre on SEC is small in the stand-alone RO process. The SEC RO reaches a minimum value of 4.0 kWh/m 3 when the water recovery rate is 51%. Additionally, SEC RO-ERD is lower than SEC RO ; hence, the addition of the pretreatment unit has a pronounced influence on the total SEC. In the RO-ERD process, SEC is maintained at approximately 2.0 kWh/m 3 when the water recovery rate ranges from 27-36%.

RO-PRO Process
The PRO unit provides both power generation and a reduction in the brine discharge concentration. Figure 7 shows the energy consumption of RO and energy generation of PRO in the RO-PRO process. When the water recovery rate is higher than 29%, both power generation and energy consumption increase as the recovery rate increases; however, the rate of increase of the latter is significantly larger. The results show that RO energy consumption is more dominant than PRO energy generation in the RO-PRO process.

RO-PRO Process
The PRO unit provides both power generation and a reduction in the brine discharge concentration. Figure 7 shows the energy consumption of RO and energy generation of PRO in the RO-PRO process. When the water recovery rate is higher than 29%, both power generation and energy consumption increase as the recovery rate increases; however, the rate of increase of the latter is significantly larger. The results show that RO energy consumption is more dominant than PRO energy generation in the RO-PRO process.  Figure 8 shows a comparison of SEC between the RO-ERD and RO-PRO proce When the water recovery rate is low, the SEC of these two processes decreases as the w recovery rate increases. Notably, the increased water recovery rate will reduce the amo of seawater required, thereby decreasing SECpre. However, when the water recovery is high, SEC increases as the water recovery rate increases. Since the RO operating p sure increases as the water recovery rate increases, HP energy consumption results i increase in the total SEC. When the RO water recovery rate is 30%, the RO-ERD pro yields a minimum SEC of 1.97 kWh/m 3 , and the RO-PRO process has a minimum SE Comparison of RO energy consumption and PRO energy production in the RO-PRO process. Figure 8 shows a comparison of SEC between the RO-ERD and RO-PRO processes. When the water recovery rate is low, the SEC of these two processes decreases as the water recovery rate increases. Notably, the increased water recovery rate will reduce the amount of seawater required, thereby decreasing SEC pre . However, when the water recovery rate is high, SEC increases as the water recovery rate increases. Since the RO operating pressure increases as the water recovery rate increases, HP energy consumption results in an increase in the total SEC. When the RO water recovery rate is 30%, the RO-ERD process yields a minimum SEC of 1.97 kWh/m 3 , and the RO-PRO process has a minimum SEC of 1.33 kWh/m 3 when the RO water recovery rate is 38% [32]. Comparing the two processes, RO-PRO has a lower SEC due to the energy generation by PRO and obtains a higher water recovery rate at the minimum SEC.
cess. Figure 8 shows a comparison of SEC between the RO-ERD and RO-PRO processes. When the water recovery rate is low, the SEC of these two processes decreases as the water recovery rate increases. Notably, the increased water recovery rate will reduce the amount of seawater required, thereby decreasing SECpre. However, when the water recovery rate is high, SEC increases as the water recovery rate increases. Since the RO operating pressure increases as the water recovery rate increases, HP energy consumption results in an increase in the total SEC. When the RO water recovery rate is 30%, the RO-ERD process yields a minimum SEC of 1.97 kWh/m 3 , and the RO-PRO process has a minimum SEC of 1.33 kWh/m 3 when the RO water recovery rate is 38% [32]. Comparing the two processes, RO-PRO has a lower SEC due to the energy generation by PRO and obtains a higher water recovery rate at the minimum SEC.

RO-MD Process
In the RO-MD process, MD can enhance the overall water recovery rate. This section will discuss the operating parameters associated with the minimum SEC and the influence of the heat source for the MD unit on SEC. Table 2 lists the parameters of the simulation. Since the RO water recovery rate (Y r ) is lower than 50%, Y r is set to 40% or 50% in the simulation. Considering the saturated concentration of NaCl, the overall water recovery rate (Y r + Y m ) of seawater in the integrated process should be less than 85%. Moreover, the liquid entry pressure (LEP) of the MD membrane should also be taken into consideration; thus, the MD operating pressure is set to 3 bar [33]. For the estimation of the MD heat requirement, EE is set to 60% [27].  Figure 9 displays the SEC RO-MD values at different water recovery rates. When the RO water recovery rate remains constant and there is sufficient waste heat as the MD heat source (Equation (15), X = 1), SEC RO-MD can be reduced as the MD water recovery rate increases. When the MD water recovery rate increases from 1% to 45%, SEC RO-MD decreases from 2.14 to 1.03 kWh/m 3 . If no waste heat supplies the MD process, SEC will largely increase as the water recovery rate increases due to the high energy consumption of the enthalpy of vaporization from the MD process. When the RO recovery rate is maintained at 40% and the MD water recovery rate is 45%, SEC RO-MD can be as high as 584 kWh/m 3 . Because the amount of waste heat supplied by industrial sites is confined, the water recovered by the MD unit in large-scale seawater operations is less than 10% of all water recovered. Thus, both the gray and blue rectangular areas in Figure 9 reflect practical conditions. Figure 9 displays the SECRO-MD values at different water recovery rates. When the RO water recovery rate remains constant and there is sufficient waste heat as the MD hea source (Equation (15), X = 1), SECRO-MD can be reduced as the MD water recovery rate in creases. When the MD water recovery rate increases from 1% to 45%, SECRO-MD decrease from 2.14 to 1.03 kWh/m 3 . If no waste heat supplies the MD process, SEC will largely in crease as the water recovery rate increases due to the high energy consumption of the enthalpy of vaporization from the MD process. When the RO recovery rate is maintained at 40% and the MD water recovery rate is 45%, SECRO-MD can be as high as 584 kWh/m 3 Because the amount of waste heat supplied by industrial sites is confined, the water re covered by the MD unit in large-scale seawater operations is less than 10% of all wate recovered. Thus, both the gray and blue rectangular areas in Figure 9 reflect practical con ditions. Figure 9. SEC of the RO-MD process with or without waste heat contribution. Figure 9 indicates that waste heat can largely influence SECRO-MD. If there is not suffi cient waste heat supply, electricity, steam or other heat sources should be provided fo Figure 9. SEC of the RO-MD process with or without waste heat contribution. Figure 9 indicates that waste heat can largely influence SEC RO-MD . If there is not sufficient waste heat supply, electricity, steam or other heat sources should be provided for the MD process. Figure 10 shows the relationship between the total water recovery rate and SEC RO-MD at different percentages of the waste heat supply when the RO water recovery rate is 40%. If the waste heat could provide 90% or 95% of the required heating energy (X = 0.9 or 0.95), SEC RO-MD could be reduced to 25.5 or 13.61 kWh/m 3 , respectively, when the total water recovery rate is 51%. the MD process. Figure 10 shows the relationship between the total water recovery rate and SECRO-MD at different percentages of the waste heat supply when the RO water recovery rate is 40%. If the waste heat could provide 90% or 95% of the required heating energy (X = 0.9 or 0.95), SECRO-MD could be reduced to 25.5 or 13.61 kWh/m 3 , respectively, when the total water recovery rate is 51%. Figure 10. SEC of the RO-MD process with or without a sufficient waste heat source for MD. Figure 11 displays an SEC comparison between the RO-MD-PRO and RO-PRO processes. It is assumed that there is sufficient waste heat supply for MD operations. The straight line indicates SECRO-MD-PRO at different water recovery rates. This figure also dis- Figure 10. SEC of the RO-MD process with or without a sufficient waste heat source for MD. Figure 11 displays an SEC comparison between the RO-MD-PRO and RO-PRO processes. It is assumed that there is sufficient waste heat supply for MD operations. The straight line indicates SEC RO-MD-PRO at different water recovery rates. This figure also displays SEC RO-ERD and SEC RO-PRO when the RO water recovery rate is 40%. When the RO water recovery rate is held constant, SEC RO-MD-PRO decreases as the MD water recovery rate increases. Moreover, when the total water recovery rate is 64%, SEC RO-MD-PRO changes from positive to negative, which suggests that energy generation by PRO is larger than energy consumption by pumping. Figure 10. SEC of the RO-MD process with or without a sufficient waste heat source for MD. Figure 11 displays an SEC comparison between the RO-MD-PRO and RO-PRO processes. It is assumed that there is sufficient waste heat supply for MD operations. The straight line indicates SECRO-MD-PRO at different water recovery rates. This figure also displays SECRO-ERD and SECRO-PRO when the RO water recovery rate is 40%. When the RO water recovery rate is held constant, SECRO-MD-PRO decreases as the MD water recovery rate increases. Moreover, when the total water recovery rate is 64%, SECRO-MD-PRO changes from positive to negative, which suggests that energy generation by PRO is larger than energy consumption by pumping. Figure 11. SEC of the RO-MD-PRO process at various MD water recovery rates and total water recovery rates. Figure 11. SEC of the RO-MD-PRO process at various MD water recovery rates and total water recovery rates. Figure 12 illustrates the FES for different water recovery rates. When the MD heat requirement is totally met by waste heat, the FES of RO-MD, RO-PRO, and RO-MD-PRO increases as the total water recovery rate increases since MD can enhance water recovery and PRO can extract osmotic pressure energy from brine discharge. At the same water production rate, when the water recovery rate increases, the amount of brine will decrease, but the brine concentration will increase. According to Equations (7) and (9), these two phenomena have opposite effects on PRO energy generation. Therefore, the increase in the rate of FES for RO-PRO is lower than that for RO-MD and RO-MD-PRO.  Figure 12 illustrates the FES for different water recovery rates. When the MD heat requirement is totally met by waste heat, the FES of RO-MD, RO-PRO, and RO-MD-PRO increases as the total water recovery rate increases since MD can enhance water recovery and PRO can extract osmotic pressure energy from brine discharge. At the same water production rate, when the water recovery rate increases, the amount of brine will decrease, but the brine concentration will increase. According to Equations (7) and (9), these two phenomena have opposite effects on PRO energy generation. Therefore, the increase in the rate of FES for RO-PRO is lower than that for RO-MD and RO-MD-PRO.  Table 3 lists the SEC values of RO-integrated processes reported in previous studies and in this study. The results indicate that previous studies mainly focused on SEC at specific water recovery rates, and few have evaluated the impact of the MD water recov-  Table 3 lists the SEC values of RO-integrated processes reported in previous studies and in this study. The results indicate that previous studies mainly focused on SEC at specific water recovery rates, and few have evaluated the impact of the MD water recovery rate on SEC. In this study, SEC at a total water recovery rate ranging from 10% to 85% was studied; moreover, the conditions associated with different MD water recovery rates were discussed. Thus, this work can provide a broader range of evaluations for process integration. Compared with Wan and Chung's results [27], SEC RO-PRO in this work is slightly higher than theirs (0.28 kWh/m 3 ) when the RO water recovery rate is 50%. This difference is mainly attributed to the different efficiency settings of the pump and PX. Compared with the work of Kim et al. [25], SEC RO-MD-PRO in this study is lower. The difference is that this study assumes that the MD heat source is completely provided by waste heat. When the ratio of MD water recovery to total water recovery increases, the SEC of the RO-MD-PRO process can be reduced.

Conclusions
Based on the principles of energy consumption during fluid transport, fluid mechanical energy recovery, MD thermal energy demand and PRO energy generation, this study established theoretical SEC models for water recovery by integrating RO with other unit operations, such as MD and PRO. The purpose of this study was to analyze the effect of the total water recovery rate on SEC in different integrated processes, such as RO-PRO, RO-MD and RO-MD-PRO. From this study, the following conclusions can be drawn.
(1) The brine flow rate decreases as RO water recovery increases. When the water recovery rate is greater than 30%, recyclable mechanical energy in the PX begins decreasing in availability, and SEC starts increasing. The minimum SEC is 1.97 kWh/m 3 at a recovery rate of 30% for the RO/ERD process. (2) For the RO-PRO process, SEC reaches a minimum value of 1.33 kWh/m 3 at a recovery rate of 38%. RO-PRO can give a lower SEC than RO due to the energy generation by PRO, and a higher water recovery rate is obtained at the minimum SEC. (3) For the RO-MD process, when the RO water recovery rate remains constant and there is sufficient waste heat as the MD heat source, SEC can be reduced as the MD water recovery rate increases. If the water recovery rate of RO is fixed at 40% and the total water recovery rate is 85%, SEC is 1.03 kWh/m 3 . (4) For the RO-MD-PRO process and an RO water recovery rate that is constant at 40%, the energy consumption due to pumping and energy generation by PRO reaches a balance at a total water recovery rate of 64%, which means that SEC is zero under these conditions. When the total recovery rate exceeds 64%, the FES of the integrated process is greater than 100%.
(5) The limit of the water recovery rate for SWRO is generally 50%. Thus, if we assume that the RO water recovery rate is 50%, the recovery rate of water for MD is 10%, and sufficient waste heat is available as a heat source for the MD unit; the SECs of the RO, RO-PRO, RO-MD and RO-MD-PRO processes are found to be 2.28, 1.47, 1.75, and 0.67 kWh/m 3 , respectively. The corresponding FES values of the integrated processes are 37%, 23% and 70% when compared with the baseline RO process.
It is noted that this study only evaluated the energy consumption for water production from integrated RO processes. When evaluating the water production cost in future engineering applications, further consideration of capital costs, such as land, construction, equipment, module, maintenance and operation costs, is needed.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: Nomenclature