Multistage Bounded Evolutionary Algorithm to Optimize the Design of Sustainable Photovoltaic (PV) Pumping Irrigation Systems with Storage

: Small o ﬀ -grid photovoltaic (PV) pumping irrigation systems with storage tanks are an environmentally friendly, cost e ﬀ ective and e ﬃ cient way of taking advantage of solar energy to irrigate crops, and they are increasingly being used today. However, ﬁnding the optimal design of this type of system is cumbersome since there are many possible designs. In this work, a new heuristic method based on the hybrid approach, which uses search space reduction, is developed and adapted to the optimal design for this type of PV irrigation system. At di ﬀ erent stages, the proposed approach iteratively combines a bounding strategy based on the application of engineering rules with the aim of reducing the search space with a genetic algorithm to ﬁnd the optimal design within this search space. The proposed methodology was applied to minimize the cost of a benchmark case study consisting of a real farm placed in the province of Almer í a (Spain). The proposed methodology was able to provide a faster and an accurate convergence due to the reduction of the search space. This fact led to a reduced total cost of the system. This study proved that the most sensitive variables were the number of modules and the type of pump, whereas the diameter of the pipe and volume of the storage tank remained more stable.


Introduction
The use of renewable energy, such as solar energy, is a suitable measure for coping with the scarcity of fossil fuels and their as sociated problems, such as the emission of Greenhouse Gases (GHG) and their growing impact on global warming and climate change. The use of off-grid photovoltaic pumping irrigation systems has increased recently not only because of their environmental benefits but also because of their economic profitability because the cost of conventional electrification in isolated rural areas is rising, whereas the price of PV cells has dropped in recent decades.
The possibility of storing both water and energy in an elevated reservoir or in the soil is another valuable feature of photovoltaic irrigation systems, although in some cases, battery banks are used to store the energy [1]. In addition, in many crops, water consumption is strongly dependent on radiation [2,3], thus PV energy production and energy consumption can be synchronized so the design of these systems can be simpler and more economically and environmentally efficient.
Different typologies of PV irrigation systems exist, namely those that pump water to an elevated tank or reservoir and direct pumping irrigation systems that transport water directly to the irrigation pipe distribution system. The latter type is less common and more complex to operate, primarily because energy production and consumption need to be matched perfectly [4][5][6] as the capacity of the soil profile to store water is limited, especially in sand dominant or coarse soils. In addition, the flowrate of these systems is usually constant because it depends on the discharge and number of the emitters of the irrigation units, so adopting strategies to adapt the flowrate to the available energy is much more complicated. Different authors have comprehensively analyzed this kind of system in order to propose operation strategies able to match energy supply and demand and optimize their performance [7][8][9][10].
The first type of PV irrigation system mentioned is the most common because the bias between the energy supply and demand can be balanced by storing water in an elevated storage tank or reservoir [11]. Despite their advantages, these systems also feature some drawbacks, such as the reservoir investment costs and higher evaporation losses from the reservoir [12].
The sizing of these systems must consider the different parts that comprise the system itself, which makes finding an optimal design a difficult task because there are many options to choose from and all the components of the system are interrelated.
Many authors have developed simulation models of the PV irrigation systems and have tried to find an optimal design by performing a sensitivity analysis. The said procedure is simple, but the number of iterations should be very high in order to find a suitable solution. For example, Glasnovic and Margeta [4] developed a simulation model in order to find the optimal solution by performing a sensitivity analysis. Yahyaoui et al. [13] developed an algorithm for sizing PV irrigation systems for tomato crops. Campana et al. [14] followed a similar approach. Bouzidi [15] also developed a simulation model of the system and used two optimization criteria: The concept of loss of power supply probability (LPSP) for the reliability as sessment and the life cycle cost (LCC) for the economic evaluation. Maleki and Askarzadeh [16] and Maleki and Pourfayaz [17] also used the LPSP concept for different kinds of hybrid PV systems. Hamidat and Benyoucef [11] proposed a systematic procedure for sizing PV irrigation systems with storage.
Due to the complexity of the problem, the use of powerful and robust optimization methods is needed for determining an optimal solution, both from an economic and environmental perspective. More specifically, the use of heuristic methods is the most suitable approach for this problem as heuristic approaches can handle complex global optimization problems. Heuristic methods can be defined as a special type of optimization method based on simple rules, which on many occasions mimic some biological or physical processes, capable of finding a "satisfactory" solution. Albeit this solution is not guaranteed to be optimal [18], it is close enough to the optimum from a practical standpoint.
Several heuristic optimization methods have been applied to optimize the design of different types of PV powered systems. Among them, Maleki et al. [19] applied different PSOs (particle swarm optimization algorithms) to the optimum design of hybrid (PV/wind/Batt) systems. Maleki and Pourfayaz [17] also applied three classic heuristic algorithms and four new metaheuristics to the optimal design of a PV/WT/battery hybrid. Sanchez et al. [20] evaluated the performance of the Particle Swarm Optimizer (PSO) method for the optimal design of a residential hybrid wind-photovoltaic-fuel power plant in Mexico.
Regarding the specific case of stand-alone PV pumping systems with water storage, Bakelli and Kaabeche [21] proposed different metaheuristics, such as ant lion optimizer (ALO), bat algorithm (BA), cuckoo search algorithm (CSA) and gray wolf optimization algorithm (GWO). They concluded that the GWO algorithm outperformed the other algorithms. Olcan [22] developed a multi-objective model that considered both the reliability and cost as objectives. On the other hand, EL-Shimy [23] preferred to apply an efficient non-conventional direct method.
Genetic or evolutionary algorithms are perhaps the most well known and widely used heuristic methods and they have also been applied to the PV pumping design problem. For example, Dufo-López and Bernal-Agustín [24] and Carroquino et al. [25] used a Genetic Algorithm (GA) to design a hybrid PV-diesel system. Ma et al. [26] also used a GA to optimize the design of a solar photovoltaic system with pumped storage. Mérida et al. [27] proposed a multiobjective approach and developed a model called MOPISS that used a Non-dominated Sorting Genetic Algorithm (NSGA-II) for the optimal sizing of a PV irrigation system. Despite their valuable advantages, the accuracy of heuristic methods and their convergence speed decrease as the size of the optimization problem increases. This inefficiency occurs because general purpose heuristic algorithms must explore a vast search space, and thereby devote a great effort to evaluating unfeasible solutions. To overcome this drawback, Reca et al. [28] proposed a new hybrid evolutionary algorithm that was applied to the optimal design of a pipe network. This approach reduced the search space thus increasing both the convergence speed and accuracy of the solution.
In this study, a new heuristic method based on the hybrid approach using search space reduction is developed and adapted to the optimal design of stand-alone PV irrigation pumping systems with a storage tank. The proposed approach is based on iteratively reducing the search space by applying engineering-based rules to the problem with the aim of reducing the number of possible scenarios to analyze. The proposed methodology has been applied to minimize the cost of a benchmark case study. The proposed case study consisted in finding the optimal design of a stand-alone PV irrigation system for an olive orchard farm located in Tabernas (Almería).

Overall Description of the Optimal Design Approach
The system is composed of a photovoltaic array, a frequency converter, a variable speed pump and an elevated tank or reservoir that stores water with enough elevation to distribute it to the irrigation system. A scheme of the PV pumping system analyzed in this work is illustrated in

117
The decision variables that must be optimized are the following:    The decision variables that must be optimized are the following:

1.
Diameter of the delivery pipe 2.
Capacity of the elevated tank 4.
Number of PV modules.
The proposed optimal design method is composed of a simulation model of the PV irrigation pumping system and an optimization module based on a heuristic approach. These models are described hereafter.

PV Subsystem
The PV subsystem model estimated the net photovoltaic power provided by the PV generator (P PV ) depending on the calculated incoming irradiance. An isotropic distribution of the diffuse radiation [29] was supposed in order to calculate the temporal distribution of the irradiance, I(t). The model by Collares-Pereira and Rabl [30] was applied to disaggregate direct and diffuse irradiance from the global daily radiation in MJ m −2 . A detailed description of the calculation procedure is further described in previous works [7,31].

Pumping Subsystem
The pumping system can be described by both the performance curves and the system curve. In this work, the Head-Discharge (H-Q) and Power-Discharge (P-Q) performance curves of the pump for its nominal rotational speed are given as the input data of the problem. These curves can be fitted to second degree polynomial functions and are shown in the following set of equations.
where: H is the pumping head (m), Q is the pump discharge (m 3 s −1 ) and P is the pumping power.
A to C and E to G are regression coefficients. The performance of the pump depends on not only the type of pump but also on its working speed. The performance of a variable speed pump can be modelled by applying affinity laws that are shown in the following set of equations: where: N 1 and N 2 are two specific pump speeds (in rpm) 1 and 2, f 1 and f 2 are the electric frequency, Q 1 and Q 2 are the pump discharges, H 1 and H 2 are the pumping heads and P 1 and P 2 are the pumping power for these two specific pump speeds. λ is the pump speed ratio. The system curve (H) is given in Equation (3): where: ∆z is the static head, i.e., elevation change between the water level in the suction of the pump and the water level in the elevated tank (m), h f is the friction head, i.e., the total head losses in the delivery pipe (m). Pressure losses were calculated by applying the Darcy-Weisbach equation [32].
where: f p is the friction factor of the pipe, which in turn depends on the relative roughness of the pipe and the Reynolds number, L is the length of the pipe, D is the diameter and Q is the total discharge. The net power transferred to the water by the pump is equal to: where: γ is the specific weight of the water (N m −3 ).
The total PV power required to drive the pump is estimated by applying Equation (6): where: η FC is the efficiency of the converter, η AM is the efficiency of the asynchronous motor and η P is the pump efficiency.

Water Balance
The irrigation water demands were calculated using the methodology proposed by FAO [33]. The crop evapotranspiration (ET c ) was calculated applying Equation (7): where: ET o is the reference evapotranspiration and K c is the crop coefficient and K r is a reduction coefficient for sparse crops with limited canopy cover. In this study, ET o was computed applying the Penman-Monteith formula. The net irrigation requirements, N n , (mm) can be estimated as the difference between ET c and the effective precipitation [34], meaning rainfall that is effectively used by the crop (P e ).
Finally, the gross irrigation water requirements (N g , in mm) can be calculated as a function of the N n by applying an irrigation efficiency (η a ) in order to take into account the water losses in the system (See Equation (8)).
A model to calculate irrigation water requirements previously developed by the authors and applied in previous works was used in this research.
The water volume (V Ri , in m 3 ) needed to irrigate a farm with an area S (ha) during a specific period i can be calculated as a function of the gross irrigation requirements for that period (N gi ) as shown in Equation (9): The water volume pumped to the reservoir (V Pi , in m 3 ) for a specific period i can be calculated by multiplying the pumping flowrate (Q i , in m 3 /s) by the pumping time (T i , in s). The volume of water stored in the tank can be modelled by applying a water balance model: where: V i and V i-1 are the volume of water stored in the reservoir in periods i and i-1, respectively, ID i is the irrigation deficit in period i. The volume stored in the reservoir must be lower than the maximum capacity of the tank.

Statement of the Optimization Problem
A multistage bounded genetic algorithm module was applied to optimize the design of the PV pumping system. The GA algorithm developed in this work is based on the GENOME (University of Almería, Almería, Spain) software model developed by Reca and Martínez [35], which was later converted in B-GENOME [28], an improved bounded version of GENOME that was applied to optimize the cost of a pipe network. This code was modified and adapted with the aim of optimizing the above-mentioned PV pumping design problem (BPV-GENOME). The latter model was developed using the VBA programming language and embedded in an EXCEL © (Microsoft Corporation, Redmond, Washington, EE.UU) workbook. B-GENOME used an efficient integer-coding scheme. Each individual (possible solution) is described by a chromosome consisting of a vector of n discrete variables (genes). Genes were codified by an integer value ranging from one to n i (maximum number of possible values as signed to the variable i, i.e., alleles). This coding scheme was also adopted in BPV-GENOME, in which case the chromosome is composed of four genes, one per decision variable. The first gene codifies the diameter of the pipe, the second gene codifies the type of pump, the third one the capacity of the tank and the last one the number of PV modules. The number of alleles per gene depends on the input data of the problem. The number of alleles of the first and second gene depends on the total number of pipe diameters in the input pipe database and the number of pumps in the pump database, respectively. The third and fourth variables have a few alleles that results from a predefined minimum and maximum value of each variable (volume of the tank and number of PV modules, respectively) and a variable step. The bounding strategy used in this work is described hereafter.

Multistage Bounding Strategy
A multistage bounding strategy is proposed in this study with the aim of reducing the search space and improving the efficiency of the GA algorithm. The initial set-up of the limits of the variables is done by applying engineering-based rules. In this initial stage, the aim is to apply relatively wide bounds, so that the resulting reduced search space includes not only the problem solution but also its surroundings, which may provide useful information. Individuals that are clearly beyond the limits and drastically violate the problem's restrictions are simply discarded to reduce the complexity of the problem. The steps in this initial bounding are described as follows: (a) Bounding the diameter series: Maximum and minimum reasonable pipe diameters are estimated depending on the flow range required by the irrigation system. The first step is to calculate the fictitious continuous flow (q), which is defined as the flowrate required to supply the maximum crop water requirement if the pumping system operates continuously during the whole day. The maximum flow is obtained multiplying the fictitious continuous flow by a degree of freedom (DF), which indicates the number of times that the system would be able to supply the required irrigation water in a day. The higher the selected value of the degree of freedom, the wider the bound for the maximum flow (see Equation (11)).
where: q (m 3 /s) is the continuous flow, DF is the degree of freedom (dimensionless) and Q M (m 3 /s) is the maximum pumped flow.
To calculate the diameter bounds, velocity restrictions are imposed on the flow in the pipes. Excessively high flow velocities in the pipe may produce pipe erosion, loud noise, undesired effects of hydraulic transients and excessive head losses. On the contrary, very low velocities may produce sedimentation and oversizing of the system. Maximum (D M ) and minimum (D m ) diameter sizes can be calculated by applying these velocity restrictions in the pipe. D M is limited to the diameter able to supply the maximum flow with the minimum admissible flow velocity in the pipe. D m is defined as the size of the pipe that would supply the continuous flow with the maximum admissible velocity in the pipe (see Equation (12)). Diameters in the pipe database out of this range are neglected.

238
(rotating speed). Pumps above this limit would produce operating points far beyond the upper limit.

239
The lower envelope requires that the pump curve falls above the lower operating point operating at   Extreme maximum and minimum system curves are defined for the maximum pumped flow (Q M ) and continuous flow (q), respectively. These curves were computed by applying Equations (3) and (4). Two extreme operating points, (Q M , H M ) and (q, H m ), are thus defined. These points are used to limit the range of admissible pumps. Admissible pumps are those whose performance curves lie within two envelopes. The upper envelope requires that the performance H-Q curve of the pump falls below the maximum operating point with the pump operating at the minimum frequency (rotating speed). Pumps above this limit would produce operating points far beyond the upper limit. The lower envelope requires that the pump curve falls above the lower operating point operating at its maximum speed. Pumps with H-Q curves below this limit are clearly unable to supply the flow and energy demands.
(c) Bounding the capacity of the storage tank: Maximum and minimum capacity levels are defined. The minimum capacity of the tank was set to the volume of water needed to irrigate for one day. This volume makes it possible to have a minimum water storage in case the PV pumping system fails. The maximum capacity was set to the volume equal to 30 days of water consumption. Higher volumes are considered excessive. The volume step is equal to one day of storage so that 30 possible storage volumes (multiple of the minimum volume to store one day of operation) are considered in this initial bounding.
(d) Bounding the PV power: The maximum and minimum number of PV modules are calculated depending on the power required to drive the most and least powerful pumps in the selected series of admissible pumps. Required PV power for both cases can be calculated by applying the appropriate inverter and motor efficiencies (Equation (6)). If the difference between the maximum and minimum number of modules is lower than 50, the step equals one module. The maximum number of steps was set to 50 so if the difference is higher, the step can be higher than one module.
After this initial bounding is performed, the GA algorithm is applied recursively (in this work ten runs of the GA were performed). The solutions of this first GA phase provide the range of solutions that serve to refine the search space and redefine the problem limits in the subsequent stage. Since the solution space is now narrower, the number of PV modules and capacity of the tank steps are now smaller. After a second round of GA analyses, a new more precise solution with higher resolution can be found. In this work, two stages were performed, although more iterative stages can be also conducted.

Genetic Algorithm
Genetic algorithms (GA) are heuristic methods that mimics the evolutionary mechanisms of natural selection and genetics [18]. Genetic algorithms are based on the effective optimization model that has naturally evolved for dealing with large, highly complex systems.
In the GA formulation, the initial population is obtained randomly. Using a random generator, integer values are generated, ranging from 1 to n j , where n j is the number of alleles of each gene. Then, population evolves by undergoing an iterative reproductive cycle in subsequent generations. This evolutive cycle is composed of three basic operators: Selection, crossover and mutation. There are many possible methods to perform these basic operations: The selection operator used in this work is called steady-state-delete-worst plan. According to this plan, a new individual is inserted in the population if its fitness exceeds that of the least fit member of the parent population. That least fitted individual is consequently removed and replaced by the offspring. The crossover process involves the exchange of information between two parent chromosomes to produce two new offspring chromosomes that inherit their characteristics. This process is controlled by the input parameter p cross , which defines the probability of crossing two chromosomes. In this work, a uniform crossover procedure was applied, meaning that the parents' chromosomes exchange information gene to gene. A gene crossing rate (r cross ) parameter is defined with the aim of defining the probability of exchanging a specific gene. Finally, a mutation probability (P mut ) is defined to specify the mutation rate. The mutation operator swaps the allele of the selected gene with a randomly generated value.
A population size of 50 individuals and 50 generations was selected in this work because it provided a good trade-off between convergence speed and accuracy. The remaining GA parameters and options used in this work are summarized in Table 1 and were recommended in other works [28,35]. Table 1. Input parameters for the genetic algorithm.

Objective Function
For every generation, the fitness of every individual is evaluated. This fitness is defined as the sum of the cost of the system (cost of the delivery pipe, cost of the pump, cost of the tank and cost of the PV plant) plus a penalty function to consider the supply deficit. A supply deficit occurs when one specific induvial (system design) is unable to supply the target water demands for irrigation (see Equation (13)).
where: FF = fitness function, C D = cost of the pipe, C P = cost of the pump, C S = cost of the storage tank, C PV = cost of the PV plant, K P = penalty multiplier and ID = cumulative water deficit in the system (sum of deficits for every study period). The penalty multiplier (K P ) may affect the convergence and accuracy of the GA algorithm, so the value of this parameter should be adjusted. Some researchers recommend different constraint-handling techniques, such as the use of variable values or self-adaptive penalty functions [36]. However, in this work, a high value (10 9 €/m) was selected for this parameter to discard infeasible solutions that imply an insufficient irrigation water supply.

Structure of the BPV-GENOME Model
The flowchart of the BPV-GENOME (Bounded Photovoltaic irrigation-GENOME) model is depicted in Figure 3.
where: FF = fitness function, CD = cost of the pipe, CP = cost of the pump, CS = cost of the storage tank, 294 CPV = cost of the PV plant, KP = penalty multiplier and ID = cumulative water deficit in the system 295 (sum of deficits for every study period).

296
The penalty multiplier (KP) may affect the convergence and accuracy of the GA algorithm, so the 297 value of this parameter should be adjusted. Some researchers recommend different constraint-

298
handling techniques, such as the use of variable values or self-adaptive penalty functions [36].

299
However, in this work, a high value (10 9 €/m) was selected for this parameter to discard infeasible 300 solutions that imply an insufficient irrigation water supply.

306
The model BPV-GENOME is composed of several modules: An input/output module, a GA  The model BPV-GENOME is composed of several modules: An input/output module, a GA module, a PV pumping simulation model and a bounding module. The model was programmed in the EXCEL © spreadsheet platform using the VBA programming language. The input module reads the general input data of the problem and the specific databases, i.e., the pipe database and the pump database. These databases are stored in EXCEL © sheets. These databases contain the basic characteristics of the pipe and pump series and their as sociated costs. The output module writes the final solution found by the model (optimal combination of pipe diameter, pump unit, storage capacity and PV peak power) and the best fitness function value for every generation in a results spreadsheet.

Case Study
A specific case study was used in order to test the methodology. This case study consisted in finding the optimal design of a PV irrigation system for an olive farm in Tabernas, in the province of Almería (Southeast of Spain). A description of the area and the characteristics of the farm can be found in Zavala et al. [37].
The location of the PV plant is 37 • 06' north latitude and 02 • 18' west longitude and its elevation is 435 meters above sea level. A fixed PV frame was considered with a tilt angle of 20 • , the azimuth is 0 • and the albedo is 0.2. Conventional crystalline cells were used in this study. The peak power per module was 275 W and the cost of the PV plant was supposedly 0.5 €/PW.
The calculated irrigation water requirements are given in Table 2 for the whole farm. The water is pumped from a well to an elevated tank with a static head of 38 m. The total length of the pipe from the well to the tank is 420 m.
The pipe series is composed of 14 PVC pipes with different diameter sizes taken from a major manufacturer. The working pressure of the selected pipes was 10 at m. The relative roughness of the pipe was 0.0015 mm. The pipe database is supplied as a supplementary material.
The pump database is composed of 229 different pump units. This database includes submerged pumps of 4 and 6 inches. These pumps differ in the impeller type and the number of impellers per pump unit. This database was taken from an international pump manufacturer and it is also supplied as a supplementary material.
The overall efficiency of the system, comprising the pump, inverter and PV array efficiencies was considered equal to 0.5. In this research, a constant pump efficiency was assumed.
For the storage tank, a reservoir constructed of excavated natural earth and lined with a plastic sheet was considered. The cost function of this type of reservoir was obtained by fitting unitary cost data provided by an experienced consulting company (see Figure 4).

Bounding Results
After applying the proposed bounding algorithm in this first stage to the optimal design of the case study, the following ranges of variables are obtained. The selected pipe series was reduced to a total of 11 admissible diameter pipes (out of 14 total diameters), ranging from a minimal inner diameter of 57 mm (PVC 63) to a maximum diameter of 290.8 mm (PVC 315). The series of admissible pumps was reduced to 161 pump units (out of 229). Regarding the capacity of the storage tank, the volume was discretized in 30 steps of 400 m 3 , which represents the volume required to supply the irrigation amount of one day of operation during the peak period. The capacity ranged from a minimum of 400 m 3 to a maximum of 12,000 m 3 (volume needed for 30 days of operation). The number of PV modules was limited to 50 cases, ranging from 0 to 1050 PV modules with a step of 21 modules.
The chromosome of every individual is then composed of an array of four genes (four integer numbers). The first gene was codified by an integer number ranging from 1 to 11 (selected pipe diameter series); the second one from 1 to 161 (selected pump series); the third one ranged from 1 to 30 (volume steps) and the last one from 1 to 50 (number of PV arrays steps). The limited search space had a total size of 2,656,500 combinations.

Fitness Function
The GA analysis performed in this first stage yielded the results shown in Figure 5 and Table 3. Figure 5 shows the evolution of the fitness function for ten iterative simulations. Convergence was quickly reached and the optimal fitness function value (optimal cost) ranged from a minim value of 34,142.2 € to a maximum of 41,731.7 €. The average minimum cost in the ten runs was 36,726.3 € with a coefficient of variation of 5.2%.

First Stage Results
The diameter of the pipe that appeared most frequently was 101.6 (8 out of 10 times). The type of pump varied among the iterations. Pumps 123 and 125 were chosen three times and other ones were selected once or twice at most, meaning the choice of pump is a sensitive variable. The number of PV modules ranged from 42 (11.55 kWp) to 105 (28.87 kWp). The number of modules that were most repeated (seven times out of ten) was 63 (17.32 kWp). The storage capacity was always 800 m 3 , which is close to the minimum storage capacity.

Bounding Results
A second stage was performed to improve the accuracy and resolution of the solution. The series of pipe diameters and pumps were not modified in this stage. Regarding the storage volume, the variable was discretized in 50 possible values, starting from 400 and continuing to 1400 m 3 with a step of 20 m 3 . The center of this interval coincides with the solution provided in the first stage (800 m 3 ). Applying a similar procedure, the number of modules varied from 40 to 140 with a step of two modules.

Fitness Function
The GA analysis performed in this stage yielded the results shown in Figure 6 and Table 4. Figure 6 shows the evolution of the fitness function for the ten iterative simulations. Convergence was reached fast after the 50 generations. The optimal fitness function value (optimal cost) ranged from a minim value of 30,994.5 € (almost 10% lower than in the first stage) to a maximum of 32,830.5 €. The average minimum cost in the ten runs was 31,773.8 € (an improvement of almost 14% of the first stage) with a coefficient of variation of only 2%. Considerable cost reductions were achieved in performing this second stage.

Second Stage Results
As expected, the optimal diameter of the pipe was 101.6 mm, precisely as in the first stage. The type of pump varied among iterations. Pump 123 was chosen three times; 167 and 125, twice; and pumps 166, 124 and 126, once. Thus, the choice of pump is a more sensitive variable. The number of PV modules ranges from 42 to 52. The number of modules that were selected more times (four out of ten) was 42 (11.55 kWp). The storage capacity ranged from 420 m 3 to 460 m 3 , although the most common capacity was 420 (6 out of ten).

404
Due to the stability and sufficiently high resolution of the solutions found in the second stage, a 405 third stage was considered unnecessary. The optimal design of the system was that which provided 406 the minimum cost in the second stage (see Table 5).

407
The proposed methodology was able to provide a much faster and accurate convergence to the 408 optimum because of the considerable reduction of the search space achieved.
409 Figure 6. Evolution of the fitness function during the first stage for ten iterative runs.

Final Optimal Design
Due to the stability and sufficiently high resolution of the solutions found in the second stage, a third stage was considered unnecessary. The optimal design of the system was that which provided the minimum cost in the second stage (see Table 5).
The proposed methodology was able to provide a much faster and accurate convergence to the optimum because of the considerable reduction of the search space achieved.
It is relevant to note that in the minimum cost solution found, the storage volume tends to be close to the minimum because, as the storage capacity increases, the raise in the cost of the tank exceeds the reduction in the pumping and PV costs. However, this low storage capacity causes much of the PV energy to be wasted, thereby yielding a lower Energy Use Efficiency (EUE). In the following section, a more in-depth analysis of the resulting EUE is conducted and appropriate measures for improving this energy efficiency are proposed.

Energy Use Efficiency (EUE)
With the aim of evaluating the performance of the system from an environmental standpoint, an efficiency index was calculated. The energy use efficiency (EUE) index is defined as the fraction of the maximum energy potentially generated by the PV system that is effectively used (see Equation (14)). The effectively used energy is computed as the potential energy gained by the water pumped to the tank.
Higher values of EUE implies a more efficient use of the solar energy. Energy losses include the energy potentially produced but not instantaneously consumed, energy losses in the motor, inverter and pump and in the delivery pipe.
where: E P is the energy effectively used to pump water to the tank (kWh), E PV is the total energy produced by the PV system (kWh), calculated by the model as described in Section 2.2.1. E P can be calculated by applying Equation (15).
where: V P is the total volume of water pumped (m 3 ) and H is the pumping head (m). The volume of water pumped for the optimal solution was 39,484.9 m 3 . The increase in the water level between the well and the tank was 38 m. The increasing of the potential energy of the water accumulated in the tank (net energy received by the water) was 4088.660063 kWh. The potential energy produced by the PV plant as a function of the irradiance was 21,657.46 kWh. The resulting EUE for the minimum cost design was thus 18.88%. These results are similar to those found by Zavala et al. [37] for a direct pumping system without storage.
As discussed before, EUE is low, which implies that much of the potential PV energy is wasted because the tank is often full due to low storage capacity. The EUE could be increased if the storage capacity were increased. To illustrate this scenario, we simulated the system performance for increased values of storage capacity. For each storage value, we determined the minimum required peak power of the PV plant and the EUE was then calculated. Figure 7 depicts these relationships. It can be seen that as the storage volume increases from 420 m 3 (economically optimal design), the peak power of the PV plant decreases and the EUE rises. To illustrate this effect, for a storage capacity of 10,000 m 3 , the peak power of the plant reduces almost half its value for 420 m 3 and EUE increases up to 40%.

448
According to these results, minimizing the system cost should not be the only objective to 449 optimize the performance of the system, increasing the EUE should be considered as well. A According to these results, minimizing the system cost should not be the only objective to optimize the performance of the system, increasing the EUE should be considered as well. A multiobjective optimization approach with two objectives, i.e., minimizing cost and maximizing EUE, is strongly recommended for future research in order to obtain balanced designs with the optimal trade-off between efficiency and cost.

Conclusions
In this study, a new heuristic method based on a hybrid approach and a search space reduction has been developed and adapted to the optimal design of stand-alone PV irrigation pumping to a storage tank. The proposed approach is based on iteratively bounding the search space by applying engineering-based rules to the problem with the aim of reducing the number of possible scenarios to analyze and applying a genetic algorithm to the reduced search space. This methodology has been applied to find the optimal design of a benchmark PV irrigation system with storage for a real farm located in Tabernas in the province of Almería (Spain).
A considerable reduction of the search space was accomplished with the proposed methodology, which also provided a much faster and accurate convergence.
The total cost of the PV irrigation system was considerably reduced. This study proved that the most sensitive variables were the number of modules and the type of pump, whereas the diameter of the pipe and volume of the storage tank remained more stable.
It has been proven that minimizing the system cost may result in low energy use efficiency (EUE). To overcome this situation, a multiobjective optimization approach considering two objectives, i.e., minimizing cost and maximizing EUE, is strongly recommended for future research in order to obtain balanced designs with the optimal trade-off between efficiency and cost.