1. Introduction
Climate change has a great impact on the water supply sector of many regions worldwide, such as in southern Europe, where the hydrological stress is expected to shortly increase due to this phenomenon and some areas are already facing serious water problems indeed [
1]. One of these is the province of Almería (southeast of Spain), which is identified as one of the driest regions in the continent, but has one of the largest agricultural production systems. What is more, the main driving force of the economy in this province is agriculture, with around 30,000 ha of greenhouse crop production [
2]. The development of this industry has been tied, for a long time, to the depletion of freshwater reservoirs (despite the efficient management of this resource that has been performed in this sector [
3]). This fact made the case for the inclusion of alternative water sources such as desalination plants, thus enhancing the availability of fresh water for the sustainability of the Almerian agricultural system [
4]. As a consequence, the current panorama in Almería can be visualized as a distributed water network composed of (i) producers, based on traditional (wells and water public utility network) and non-traditional sources (desalination facilities and other non-conventional systems), and (ii) consumers, such as industries related to agriculture and greenhouses.
As in any other kind of distribution network or multi-agent system, this requires integral and optimal management [
5,
6,
7] and, until now, different approaches dealing with this issue have been introduced and formulated in the literature. For example, in the studies of Ocampo Martinez et al. [
8] and Pascual et al. [
9], Model Predictive Control (MPC) formulations were proposed for the efficient management of the urban water cycle of Barcelona (Spain) in terms of operating costs. A similar approach was employed by Lopez Farias et al. [
10] but, in this case, to improve the forecasting of the control method. In the work presented in Ref. [
11], a distributed MPC approach, which was aimed at decreasing the required resolution time and computational cost, was put forward for the same problem. Another interesting strategy was presented in Ref. [
12], where an MPC controller was in charge of optimizing the energy–water nexus in urban water networks. Furthermore, a scheduling method was formulated by Zhang et al. [
13], in which the distributed water network was modeled and divided into three levels: water supply source, water station, and water user. Then, the developed model was used to implement effective scheduling strategies for the guidance of regional water distribution systems.
Although all the aforementioned studies presented satisfactory results, the proposed techniques are focused on optimizing the transport water network by minimizing the operational costs. However, in most cases, not only the water needs to be optimally managed but also energy trading. The presence of desalination plants and other non conventional water sources introduces new agents of intermittent nature in the problem, as they are normally powered by renewable energy to improve their efficiency [
14,
15] and to reduce water costs [
16]. In addition, adequate storage systems [
17] and the use of the public utility network as a backup are required for their continuous operation. These elements, together with the related industries and greenhouses’ necessities, make the whole system to constitute an agro-industrial district [
18] in which the generation, storage, and distribution of heterogeneous resources must be optimally performed for its efficient proper exploitation.
Over the last few years, some studies dealing with the optimal management of multi-agents systems that include non-conventional water sources, such as desalination plants, have been published. Most of these analyses are based on the typical approach for an energy hub (EH), which relies on a simplified modeling methodology that represents the interactions given inside manifold systems, attending to its input–output configuration [
19]. Gharffarpur et al. [
20] proposed a scheduling method to manage an isolated district, including a Reverse Osmosis (RO) plant, renewable energy generation, and electricity, heat, and water demands. However, their formulation did not consider the connection to the public utility network and the objective of the scheduling method was not aimed at improving operating costs, but it was devoted to maintaining some levels of resilience. Authors in Ref. [
21] introduced an alternative modeling methodology to integrate energy and water systems that include multiple energy sources. This was then used to pose a Mixed Integer Nonlinear optimization Problem (MINLP) and to perform the optimal management of the different systems comprising a shale-gas production plant, but no connection to the public network was taken into account either. In the work presented in Ref. [
22], a similar procedure was also proposed to manage an isolated micro-grid located in an island, considering the optimal dispatch of water and energy. In this case, the planning method was tasked with minimizing the environmental pollutants as well as the operation and investment costs. Once more, the approach is only valid for off-grid regions (without public utility network connections) such as islands. A more complete approach was recently published by the authors in Ref. [
23], in which instance the multi-agent system also included issues related to energy and water, but both conventional (i.e., a well) and non-conventional water sources (i.e., a desalination plant) were considered. The management method was also based on other paradigms of energy hubs’ and the planning problem included a multi-objective optimization algorithm that minimized, at the same time, the energy and the fresh water extracted from underground reservoirs. Although in this case a connection to the public utility network was included, the sale of surplus resources was not considered in the problem, which might prove to be a determining factor for the correct economical exploitation of agro-industrial districts.
Accordingly, based on the above review, the main gaps identified in the literature that are the basis of this manuscript are listed as follows:
At present, the management strategies presented for distributed water networks are mainly focused on the optimal performance of the transportation system, without contemplating the sources of water, and without considering any other resource apart from water in the problem.
Regarding the strategies based on multi-agent systems, they are mainly focused on isolated plants. In this way, trading with the public utility network is not usually considered, neither are the sale of surplus resources nor backup systems to carry out continuous operations, which can be fundamental in agro-industrial districts to maintain the desired level of quality and productivity.
To cover the above issues, this paper extends the case study addressed in our previous work Ref. [
18] by including all the plants related to water issues (as well as the exchange of resources among them) contemplated in the CHROMAE research project (“Control and Optimal Management of Heterogeneous Resources in Agroindustrial Production Districts Integrating Renewable Energies”,
www2.ual.es/chromae (accessed on 25 January 2021)). The resulting system is in fact an agro-industrial district in which several water sources and resources of different nature must be managed. In addition, the connection to the water public utility network and the sale of surpluses is considered. Unlike the approach presented in Ref. [
18], which was focused on the development of an enabling platform for the management of the distributed facilities, and presented a basic MPC strategy for solving the problem, in this work we use the approach presented in Ref. [
24] but adding new capabilities to the energy hub model. For example, variable dependent loads are considered by introducing a term in the model that reflects the variable production of the facilities. This is especially important to the operation of desalination facilities in these kind of environments as it allows the management strategy to adapt their production to the operating conditions at each sampling time (water needs, availability of resources, etc.), and therefore to perform an optimal dispatch. By using real historical meteorological data and a tool developed for such problems (ODEHubs), a simulation is presented to exemplify this fact and then compared to a manual operation, evidencing the benefits achieved by the proposed strategy.
The rest of the paper is arranged as follows.
Section 2 presents the case study, the energy hub modeling approach, and the particularization of this model for the analyzed system.
Section 3 presents the simulation results and the comparison with the manual strategy. Finally,
Section 4 discusses about the main outcomes of the work.
3. Results and Discussion
The analysis presented in the following subsections consists of a comparison between the operation of the EH, over a week in March 2014, under two different conditions placed on the desalination plant: Case 1 considers the operation point, in terms of the desalinized flow of water (
), as a decision variable; whereas in Case 2 that variable is constrained to force that the plant work at its maximum capacity whenever it is activated (see
Table 3), which is representative of what a human operator would do. ODEHubs was used to solve the optimization problem presented in Ref. [
30] and above adapted for the district described in
Section 2.1. The toolbox was configured to use a scheduling mode strategy (fixed horizon, with deterministic and measurable disturbances) and the solver
intlinprog [
32] and it was executed on an Intel
® Core
TM i7-6700K 4GHz CPU, taking around 96 s to simulate each of the two cases considered.
Figure 7 and
Figure 8 correspond to Case 1 and
Figure 9 and
Figure 10 correspond to Case 2, in which the results are arranged in the same way, that is, using a think solid line to represent each of the hourly demand profiles (outputs of the energy hub that form
) together with a colored stacked bar graph showing the sources (elements of
) that meet those demands. The solid thin lines depicting the evolution of the sold resources represent in fact the accumulated values of
and
, respectively, hence the actual values of
and
can be obtained by subtraction. All these elements are expressed in terms of power or flow according to the scale of the left vertical axis, whereas the right one is only employed for the dashed lines that symbolize the state of charge of each of the storage systems (elements of
), in terms of energy, volume, or mass. Note that both the charge (
) and discharge flows (
) are deliberately missing in the said figures, since they can be deducted from either the slope of the dashed lines or the difference between the colored stacked bar graph and the stacked solid lines. The results have been split into non-dependent outputs (
Figure 7 and
Figure 9), i.e., those ones related to loads or demands that cannot be controlled and would need to be met ad hoc; and dependent outputs (
Figure 8 and
Figure 10), which hinge on the on/off state of some of the devices.
Regarding the similarities in both cases, the broad strategy to schedule the dispatch of resources over the week consists in making use of the solar energy to yield either energy or water from dawn to dusk at a lower cost. This is close to zero, since radiation was considered to be freely available (
), but not exactly zero because sometimes using electricity or biomass compensates the cost of acquiring water from the public network, which is the case for the nanofiltration plant on 16–20 March (
Figure 8 and
Figure 10). In Case 1 (but not in Case 2), this also happens to the desalination plant on 18 March night (
Figure 9) and on 17’s dawn and dusk, when the modules were allowed to operate at partial load. As a result that the electricity price tends to be higher around midday, except for 20 March (see
Figure 5), and the electricity demanded by the facilities (CIESOL and the greenhouse) remains quite stable at about 6 kW, most of the electricity produced by the photovoltaic field is either directly sold or stored to be sold straightaway after noon, when the price is still profitable (see
Figure 7 and
Figure 9). That is also the reason why no electricity is kept stored at night, given the abundance of solar radiation the following days and the lower price at that time.
and
are managed analogously, but they show more variability in the demand. Water (
) is required especially during the working time and when the evapotranspiration rises (as outdoor temperature and irradiance do), which happens to be higher on 17–20 March (hence, storage is required during the previous night); whereas the thermal needs (
), which are nil on 16 March, exhibit a quite irregular pattern.
On the contrary, the differences found between Case 1 and Case 2 are mainly justified from
Figure 8 and
Figure 10. Since the desalination plant is no longer able to work at a partial load in Case 2, the amount of excess water that would be produced in comparison to Case 1 does not compensate the additional energy required for it. In other words, as the levels of irradiance and stored energy close to the dawn and the dusk are not enough to cover the needs of the plant operating at its maximum capacity, owing to physical constraints, it is preferable to use water from either the public network or the nanofiltration plant. This fact is noticeable in the increase of water coming from these sources in
Figure 8 and in the amount of electric energy supplied to
on 20–22 March (
Figure 10). For the same reason,
Figure 10 shows a minor usage of
, in comparison to the same variable in
Figure 9. Finally, why a first glance might look like a baffling performance on 18 March, when
is no longer used to store heat, is explained because in Case 1 the biomass boiler is used to feed the desalination plant (see
Figure 9), whereas in Case 2 it is used to feed the facilities, but the amount of consumed biomass is altogether the same, as summarized in
Table 6.
In addition to the above discussion and figures, the accumulated amounts of resources involved in each case are summarized in
Table 6 for vectors
,
, and
, as well as the total cost of purchasing or selling resources according to vectors
and
. Note that Case 1 and Case 2 differ in
,
,
, and
, which depend on the decision variables that minimize the operation cost of the plant. As shown in
Figure 7 and
Figure 8, the constraint placed on the desalination plant in Case 2 produces a shift of supply sources and more water is acquired via nanofiltration or the public network. Thus, the flexible operation of Case 1 results in a higher profit of 4.68%, a decrease in the water consumption from the public grid of 1.8 m
3 (58.1%), and an increase in the desalination plant’s production of 20.3%, in comparison to Case 2. Note that, in addition to the economic benefits, if this is extrapolated to the annual operation of the entire Almerian region, it would have a significant impact on the environment owing to the amount of water that would not be extracted from the sweetwater reservoirs.
4. Conclusions
This paper compares two cases of water management where the operation of a desalination plant is constrained to work similarly to a manual operation mode, in contrast to its flexible use, in which the amount of distilled flow is adapted to the consumption needs. The simulations performed on a realistic test-bed plant, which can be defined as an energy hub, show that those constraints actually make the system to economically underperform. It also proves ODEHubs to be suitable tool for resource scheduling problems, since the results for both cases are coherent with respect to the intuitively expected behavior of the controller regarding the economic objective and they constitute a paradigm of device-dependent variable loads, whose theoretical framework was presented in Ref. [
24] but not exemplified by a real-world study as the one above-presented. Although the scheduling strategy was used in this paper to demonstrate the differences in the flexible operation—deterministic scenario where all the disturbances are beforehand known, without a feedback loop—ODEHubs allows the users to simulate MPC-based strategies with receding horizon, which would take into account the presence of uncertainty in more realistic scenarios.
Some related studies have shown a similar performance when operating distillation facilities in terms of economic costs, as in Ref. [
33], in which the cost per unit of demanded water was 0.44 EUR/m
3. If compared with the results obtained, considering the same electricity and water prices as in Ref. [
33] together with the material and energetic needs in
Table 1, cost results in 0.12 EUR/m
3 for Case 1 and 0.14 EUR/m
3 for Case 2. This difference is due to the fact that in this study the distillation unit is enforced to yield water at the optimal operation point. Note that, in both studies, the thermal cost is neglected because it is assumed to be met by the field of solar collectors.
On the other hand, these analyses and methods are applicable to the wider Almerian ecosystem, in order to manage sets of greenhouse and water producers distributed over certain area. The modifications to be performed on the ODEHubs’s model and on the code itself have been already posed and are under implementation within a framework based on game theory in which each of the prosumers constitute different players with, at times, opposing interests.