Realistic Optimization of Parallelogram-Shaped Offshore Wind Farms Considering Continuously Distributed Wind Resources

: Offshore wind power plants are becoming a realistic option for the renewable production of electricity. As an improvement tool to the proﬁtability of OWFs, this work presents the ﬁrst complete non-genetic (and non-binary) evolutionary algorithm to optimize the location, size and layout of a parallelogram-shaped offshore wind farm, as the arrangement that is becoming an standard for offshore wind farms. It has been tested in the HR-I site. Most relevant economic data inﬂuencing the investment proﬁtability have been taken into account. In addition, the paper introduces a new approach to offshore wind farm optimization based on a continuous behaviour of varying wind conditions, which allows a more realistic estimation of the energy produced. The proposed optimization approach has been tested based on the available information from HR-I. Obtained solutions present similar values to the actual offshore wind farm in terms of investment and annual energy produced, but differs with respect to the optimal orientation and proﬁtability. The contributions of this paper are: it details the ﬁrst method to interpolate a continuous distribution of wind rose and Weibull parameters; it presents the ﬁrst algorithm to obtain a realistic optimal solution to the location+sizing+micro-siting problem for regular arrangements; it is prepared to work with the most complete set of economic, bathymetric, and wind data.


Introduction
Offshore wind technology is becoming established as a realistic option for the production of electricity from renewable sources. The accumulated power capacity at the end of 2019 was 28.3 GW, of which 22 GW corresponds to Europe [1]. This represents only 0.3% of global energy production, albeit swith a growing role in the world electricity market.
The knowledge acquired in the development of previous power plants, and the determination of the sources and severity of the contingencies that delay the development of the project, make offshore wind farm (OWF) projects an investment that is less and less risky and therefore more attractive to potential developers.
However, it is yet one of the most expensive technologies in terms of cost-per-MW installed, and a thorough profitability analysis is necessary before considering the realization of an offshore wind farm. In this analysis, the site conditions must be taken into account in order to refine the cost of the investment and deduce the best arrangement of the turbines to maximize the production of electrical energy and minimize the investment.

Genetic and Non-Genetic Evolutive Algorithms
There are numerous papers that propose an algorithm to solve the micro-siting problem and obtain the optimal arrangement of the turbines in order to maximize the annual energy production (AEP). Serrano in [2] and Herbert-Acero in [3], among others, conduct a deep review of the literature about the wind farm optimization problem (WFOP). It should be noted that the vast majority of previous works that have tackled the problem of optimal micro-positioning of wind turbines have resorted to meta-heuristic optimization techniques to solve the problem. This is because it is a problem with multiple optimal solutions (convexity), and it is not possible to describe its objective function in a completely analytical manner, which makes the problem non-derivable, thus preventing the use of classical gradient-based optimization techniques.
Kirchner-Bossi [4] uses two GA parametrizations to determine the position of the turbines that optimize the energy produced, assuming that the turbines are restricted to the area occupied by Horns Rev I (HR-I) and for Princess Amalia. Gao et al. [5] presented a multi-population genetic algorithm for a case study optimization applied to an offshore wind farm in Hong Kong. A similar study was also performed by Pookpunt and Ongsakul [6], in this case for an offshore wind farm in Thailand, by means of a binary particle swarm algorithm. Chowdhury et al. [7] introduced a particle swarm algorithm to optimize the location of wind turbines and considered different wind conditions over the study area. Furthermore, a particle swarm optimization algorithm was proposed by Hou et al. [8] to address the design of large offshore wind farms. Salcedo et al. [9] relied on a novel coral reef algorithm with the objective of maximizing the output of offshore wind projects. Feng in [10] uses a random search algorithm to maximize the power production in HR1. A hybrid method based on a genetic algorithm (GA) combined with local search was proposed by Mittal et al. [11]. A newly reported mutation operator was developed by Mayo with the objective of maximizing energy yield [12]. In a similar way, but referring to large wave farm consisting of fully submerged wave energy converters, Neshat et al. [13] use different approaches to evolutionary algorithms to maximize the harnessed power. Five of the approaches fall into the group of continuous methods, like the one used in the present work, and four of them in that of the discrete methods, of which GAs are its greatest exponent.
After reviewing these works and other similar ones, we found that they are mostly based on two approaches to model the WFOP: • Defining a discrete computational domain in the form of a grid in which the turbines can be located in the centre of each cell. This is the most popular method. For example, Kunakote et al. in [14] conduct a comparative performance of twelve metaheuristics to optimize the layout of a wind farm after discretizing the studied area into 10 × 10 grids. Binary reasoning is applied to discriminate whether each cell includes a turbine or not. • Assuming a predefined number of wind turbines in the wind farm and optimizing the value of the two geographic coordinates of each wind turbine.
These approaches generally result in power plants with very irregular layouts. This characteristic is not desirable in real OWFs due to the problems that it originates in navigation and the access to turbines for maintenance and repair.
Despite this, little work has been done to limit solutions to OWFs with evenly distributed turbines [15][16][17][18], which is a solution that is increasingly being adopted by designers of real OWFs. Among them, parallelogram-shaped wind farms are a particular and frequent case. In cases with a very limited concession area where the optimal solution will hardly have a regular shape, Serrano et al. in [15] uses a search algorithm for a parallelogram-shape solution to determine the distance between turbines, orientation and angle of the parallelogram that maximize energy production. Neubert in [16] uses a deterministic algorithm to investigate a discrete number of cases with the motivation of achieving an appealing visual impact and presumably reduce the infrastructure cost for cables. Feng in [17] uses a parallelogram shape solution in a first stage of the WFOP to deduce the best orientation and shape of the OWF, before performing the complete layout optimization. Finally, Stanley in [18] uses five parameters to define the solutions bounded in a concession area with the benefit of reducing the computational expense. For other renewable sources, a continuous exploration of the position and orientation of a solar array is presented in [19].
1.1.1. Contribution 1. Non-Genetic Evolutionary Algorithm Searching in a Vast Concession Area By exploiting the complexity reduction derived from considering parallelogram-shape layout solutions (with evenly distributed turbines), this work introduces a complete WFOP that searches the optimal values for position, size and layout of the OWF in order to maximize the project profitability. Contrary to existing works based on a GA, the proposed method can search within a vast concession area without compromising the computation time. Such a complete search spectrum is only possible if the idea of using GAs to independently position each turbine is abandoned. Instead, this work proposes a non-genetic evolutionary algorithm (NGEA) that works with eight continuous variables regarding the centre coordinates, number of rows, number of turbines per row, distances between rows and between turbines in a row, row orientation and parallelogram angle.
As a relative drawback, the turbine locations cannot freely be located in the OWF area and are forced to be uniformly distributed, thus reducing the spectrum of possible candidates to the optimum. However, optimizations made with and without the restriction of uniformly distributed turbines show a very small difference in the resulting net energy yield, less than 1% [4,16,17].

Contribution 2. Analysis of the Optimal Orientation for Rows and Columns
The proposed method has been tested in the popular HR-I site (see Figure 1). This OWF has been used in many research works as a testing bench to prove the validity of a proposed optimization algorithm, or as a method to evaluate the wake losses. For the first objective, the output of most of the tested algorithms was a layout with turbines irregularly distributed along the concession area, and hence no useful information resulted about the preferential directions of the turbine alignments. However, few research papers have been found discussing whether the choice of orientation and inclination of the concession area was appropriate. Gerdes in [20] stated that the turbine spacing was chosen as 560 m both between rows and along a row (equivalent to seven rotor diameter of the upwind turbine (D)), based on a wide wind direction distribution that revealed no clear preferable orientation and demanded an equal-distance spacing for both wind farm axes. The report [21] exposed that the modeled wake losses varied by only 0.11% over the full range of grid rotation, which was relatively unimportant in determining the total plant losses.
The present work analyzes the optimal values for orientation and inclination for a parallelogram-shaped OWF and compares the result to the actual layout. An additional advantage of using uniformly distributed turbines is that two of the variables defining the optimum, orientation, and inclination provides the designer with information about the preferred directions to align the turbines with the aim of being used in future for the selection of the concession area that is most attractive to developers. If the concession area is large enough, albeit with limitations in the number of turbines or MW to be installed, the algorithm also finds the optimal solution by completely defining its location, distance, orientation, inclination and layout (distribution in rows and columns).

Wake Loss Model
Regarding the models used for the evaluation of the energy produced, Archer et al. in [22] reviews and compares six popular analytical wake loss models. Among them, the vast majority of published works use the Jensen model [23] (also named as the PARK model) to calculate the wake effect at distances in excess of three turbine [24] diameters. This is the analytical model used in commercial software packages such as WAsP, WindPro and Meteodyn to evaluate the production of wind farms. In addition to the analytical models, there are others based on computational fluid dynamics (CFD) that allow a more accurate assessment of the effect of wakes. Walker et al. in [25] benchmarks the performance of four commonly used wake models against production data from five offshore wind farms. Nevertheless, the very high computational cost required makes their applicability to the optimal wind turbine micro-siting problem unfeasible in practical terms.
On the other hand, due to unavailability of the complete wind series or difficulty in their treatment, wind roses or Weibull parameters distributed in few sectors are generally used. The consideration of a discrete wind rose leads to the fact that the wake effect for incoming wind directions close to the transition zone between two sectors is not properly evaluated, which in turn causes discontinuities and unrealistic estimates of the produced energy. This aspect is even more important in the case of large offshore wind farms based on parallelogram-shaped layouts, since the errors made by concentrating all the probability on the bisector of a wind rose sector are amplified due to the high degree of symmetry, being successively reproduced for all the rows or columns.

Contribution 3. Continuous Interpolation of Sectorized Wind Data
For this reason, it is necessary to transform a sectorized wind rose into a continuous probability distribution, and analogously for the Weibull parameters. In the reviewed examples [26], the estimated distribution matches the available wind rose fairly well. However, instead of looking for a high match, a more realistic distribution will be estimated that yields the same mean value or the same cumulative probability at each sector.
To sum up, the literature review shows the existence of two gaps in the current body of knowledge on the problem of optimal micro-siting of OWFs. On the one hand, there are no previous works published that deal efficiently with the problem of optimal micro-siting of offshore wind farms based on parallelogram-shape layout, which is also one of the most common types of configurations. On the other hand, all previous works on optimal wind farm design have considered a discrete wind rose (usually composed of 8, 12 or 16 sectors), which can lead to errors in the evaluation of the wake effect with respect to a realistic situation in which the wind direction has a continuous evolution. Therefore, the present work introduces two significant improvements by developing an efficient optimization approach for parallelogram-shaped offshore wind farm designs and by introducing a new realistic model for evaluating the energy produced based on the continuous evolution of the wind resource.
The remaining of the paper is structured as follows: • Section 2 lists the expressions for two of the most popular economic indicators used to assess the profitability of a investment, presents a brief but complete description of the Park model to evaluate the wind speed deficit, introduces a new method to interpolate a continuous distribution of the wind rose and the Weibull parameters when data for a reduced number of sectors are given, and presents the basis of the NGEA used to solve the WFOP. • Section 3 presents and discusses the results obtained after testing it in the HR-I site; • Finally last Section exposes the conclusions extracted from the work.

Economic Indicators
Rodrigues in [27] explains different economic functions to evaluate the profitability of a prospective investment over time. Among them, the levelized cost of energy (LCOE), the internal rate of return (IRR), the discounted payback time, the return on investment, or the benefit to cost ratio are used as measurements of the operation risk and the investment quality. On the other hand, the net present value (NPV) is suitable to choose between mutually exclusive options for OWF layouts and hence should be used as objective function in optimization algorithms, although limiting the maximum investment to a certain value.
Since NPV is the difference between the present value of a company's cash inflows and the present value of cash outflows over a given time period, it can be expressed as: where T is the utility life-time, AEP is the annual net energy production (after discounting the total losses due to wake interferences and in the electrical infrastructure), p kwh and ∆p kWh are the energy price and its annual increment, OPEX and ∆OPEX are the operation and maintenance cost and its increment, CAPEX is the total investment cost and D ec is the on-time decommissioning cost after deducing the residual price of the plant. Finally, r is the discount rate, obtained from the nominal interest rate r i and the inflation in f as The annuity factor (a) can be defined as and LCOE can be expressed as In some research works, the discount rate is assumed to be similar to the weighted average cost of capital (WACC) and slightly higher to adjust for risk. WACC is more difficult to obtain [28], although reference values are given in [29] for Denmark, oscillating between 7.4% in 2001 and 8.9% in 2013. In this work, r i = 9.4% and in f = 1.5% have been used, yielding r = 7.78%.
The main drawbacks of NPV are: • A positive value is indicative of a profitable investment but does not give an idea of its quality since, in general, NPV will increase with the OWF size. It cannot be used to compare two OWFs with different sizes, locations or bid prices. • The value of LCOE is strongly influenced by the discount rate, while the value of IRR is influenced by the sale price of the energy produced. The calculation of NPV is determined by both data, and the variation in one of them will not only produce an erroneous value in the objective function but will also distort the selection between individuals. For example, if a too low energy price or a too high discount rate have been taken, then the best solution found will tend to reduce the distances between turbines more than what is ideal.
A value of 58 €/MWh can be assigned to the energy price, as the sum of the feed-in tariff and the market price (at least during 10 years) [20,30]. Nevertheless, with this energy price, profitability is doubtful. As a consequence, LCOE has been chosen instead of NPV or IRR as a objective function, since this indicator does not depend on the energy price.

Annual Energy Production (AEP)
The annual energy produced by a wind power plant can be obtained by adding together the energy, E k , produced by each of the N t turbines, for every wind direction ϕ and for every wind speed u. Throughout a year (T y = 8766h), AEP can be expressed as: where u cut−in is the cut-in speed, u cut−out is the cut-out speed, P c is the electrical energy at a certain wind speed given by the turbine power curve and f r(ϕ, u) is the density of probability at speed ν and wind direction ϕ. In this expression, the power produced by a turbine k of the array depends of the effective speed of the air flow that reaches the k-th turbine, ν k . Due to the wake effect, this effective speed is smaller than the free-flow air speed u. There are different methods to estimate the wind speed deficit caused by the wakes. The simplest and most popular wake model used by researchers in WFOP is the PARK model [31]. It was first proposed by Jensen and Katic [23], and starting from a diameter of the wake (D w ) equal to the rotor diameter of the upwind turbine (D), it assumed a expansion of the wake behind the upwind wind turbine linear with the projection of the distance between turbines onto the wind direction (s) (see Figure 2) The proportionality factor is the wake decay constant (k w ), which can be expressed as a function of the tower hub height (h) and the surface roughness length (z 0 ) For the first evaluations, a value of 0.0525 was used, obtained from z 0 = 0.005 (as usual for offshore sites) and h = 70. However, the resulting park efficiency for the HR-I layout yielded 91.5%, higher than the observed one (89% in [32]). In fact, the value obtained from (7) should be adjusted to take into account the turbulence intensity TI, considering that the lower the stability, the higher the k w [33]. When TI is not available, authors fit k w so that the results obtained from the wake loss model match the observed values, or more frequently, they use k w = 0.04 for offshore sites [34]. This is the value used in this work.
According to the momentum theory, the wake expansion produces a wind speed deficit δ that depends on the upwind thrust coefficient where v is the wake velocity at the downstream position and v ∞ is the undisturbed wind speed.
In the case of turbines that are perturbed by more than one upwind turbines, the effect of the wakes must be added to obtain the effective wind speed deficit. Different methods of aggregation have been suggested [31], and among them, the preferred one is the root sum of squares which obtains the total wind speed deficit at any turbine as the root mean square of the deficits at downstream place due to the single wake from the N upstream turbines. In the work of Tian et al. [35], this combination method is the one that exhibits the best match with the observed power data. The previous expression can be easily modified to account for the situation in which the downwind turbine is only partially affected by the wake. Figure 2 illustrates this situation that allows the calculation of the effective surface of the downwind rotor affected by the upstream weak (A e f f ) once the projection s and the transversal separation between turbines (y) are known.
w +y 2 2 y R t (11) and expression (9) for the wind speed deficit turns into where A e f f i is the surface of the partial wake of the upwind turbine i affecting the downwind turbine under study.

Continuous Distribution of Magnitudes Given by the Wind Rose and Averaged Values
Inspecting (5) makes it clear that accurate data for the density of probability f r(u, ϕ) are mandatory for an adequate AEP evaluation. Complete wind data sets are frequently unavailable to researchers or for a first site assessment. It is more common to have such data as values collected and summarized in a wind rose or as a set of mean values for each of the ns sectors under study.
In the case of HR-I, several sources can be found with different values, although they share the existence of a prevailing wind coming from NW. This work will use the ones obtained from [36], represented in Table 1, resulting in the wind rose of Figure 3. It is necessary to adjust the scale parameter A from the measurement height (z re f = 62 m) to the hub height (z hub = 70 m) due to the wind shear effect. The relationship between scale factors, and in general between wind speeds, at different heights is given by As previously mentioned, the usual value of z 0 for offshore sites is 0.005, which is also consistent with the wind profiles presented in [36].  Figure 3. Wind rose obtained from [36]. Only values between cut-in speed and cut-out speed are represented.

3%
Logically, these data do not present a continuous evolution, and in general there appear abrupt changes in the edge from one sector to another. This unrealistic distribution prevents a proper evaluation of produced energy, leading to inconsistent situations and solutions that are far from the real one. It is therefore necessary to obtain a distribution that, although estimated, is more realistic and consistent with the continuous distributions of parameters and probabilities that are assumed for the site. Feng [26] proposed a way to obtain a continuous distribution starting from cubic spline interpolations. The result effectively describes an evolution that passes through the bisectors of the ns sectors, but in general it will give rise to distributions that do not have the same mean value (in the case of the Weibull parameters), or the same cumulative probability (in the case of the wind rose) as the starting data.
Instead, this article proposes a way to obtain a continuous distribution from the averaged data for the different sectors, while respecting the aforementioned concerns. In the following, we will assume that the average value for a certain magnitude x is available for the ns sectors. At the end of this section, we will recover the Weibull parameters as this magnitude. For simplicity, we define a sector s i as although the explanation can be easily extended to s i = 2π ns (i − 1.5), 2π ns (i − 0.5) in order to centre the first sector at the origin (0 deg). The average value for a sector i can be expressed as: where x can be considered as unavailable and hence unknown.
Let us define the uniform knot vector, and a vector of ns + 1 ordinates where for convenience, we can consider or defineȳ 0 =x 0 := 0. A cubic spline S(ϕ) can be obtained as a piecewise function defined from a series of polynomials passing through the ns + 1 points (t i , y i ), i.e., satisfying On the other hand, since S(t) is a cubic (i.e., smooth) spline, its derivative exits, is finite and continuous ∀ϕ ∈ [0, 2π]. By calling it yields that, along a sector i, the mean value forx is with the result that thatx is a coherent estimation of the unknown x.
Consequently, for a certain magnitude x, by finding the derivative of the cubic spline of the sum of its average values, it is possible to obtain a continuous distribution of this magnitude while keeping its mean along a sector. The procedure is: • Calculate the accumulated sum of the average values from (17). • Obtain the cubic spline S satisfying (19). • An estimation of x can be found as the derivative of S.
This procedure can be applied to Weibull parameters, starting from their mean values for each of the sectors. Figure 4 shows the original discrete values for the Weibull parameters, together with their obtained continuous distributions.  The case of the wind rose is conceptually different, because the available data are the probability PS i that the wind comes from sector i, and the goal is to obtain a continuous distribution with the same accumulated probability at the end of each sector. However, the procedure is analogous, with the only difference of defining y i as with where p(ϕ) is the probability density that the wind blows in the direction ϕ and PS i is the probability that the wind comes from sector i, defined by the wind rose. The estimated probability densityp is the derivative of the cubic spline, and satisfies that its aggregated value along sector i is Therefore,p is a coherent estimation for the probability density.
The procedure in the case of calculating a continuous probability distribution of the wind speed is: • Obtain PS i from the wind rose. • Calculate the accumulated sum of the probabilities along the sectors from (22). • Obtain the cubic spline S. • An estimation of p can be found as the derivative of S.
With this procedure, a continuous distribution is obtained in the interval (0, 2π) as shown in the central plot of Figure 5. However, a traditional spline interpolation does not take into account that p(0) = p(2π). To fix this, the procedure has been repeated ns = 12 times, exploiting the periodicity and starting the interpolation at different angles (0, 30 deg, 60 deg . . . , 330 deg). Then, the results have been averaged to give rise to the final distribution plotted at the right of Figure 5.

Description of the Non-Genetic Evolutionary Algorithm
The aim of the WFOP is to obtain the optimal position and layout for an OWF in the concession area, such that it optimizes a certain objective function, which in this work is the LCOE. Unlike almost all of the articles reviewed, it does not seek to position each of the turbines independently [5][6][7][9][10][11][12]14,26,27], which would lead to irregular layouts. Instead, in this work, the individuals of each generation are parallelogram-shaped OWFs in which the turbines are evenly distributed. It is defined by the following 8-tuple: coordinates X and Y of the centre; number of rows (n r ); number of turbines per row (n tr ); distance between rows, in diameters (d r ); distance between turbines in a row, in diameters (d tr ); orientation of an array with respect to North, in deg, sense CW (θ); parallelogram angle; equal to 90 deg for rectangle (φ). Figure 6 represents these decision variables.
(X c , Y c , n r , n tr , d r , d tr , θ, φ). Therefore, the proposed non-genetic evolutionary algorithm (NGEA) must obtain these eight decision variables defining the OWF that minimizes the LCOE. Among them, six are real numbers and two (number of rows and columns) are discrete values. These two variables do not add additional complication because they can easily be converted to continuous values. Thus, the algorithm can be considered to work with continuous values.
In general, the following limitations apply to the presented algorithm: • All of the turbines are of the same type and the rotors are at the same height. Therefore, it is only valid for offshore sites. • Turbines are uniformly distributed forming a parallelogram-shaped OWF. All of the rows have the same number of turbines. • The cable size and the substation position are optimized, but not the connection between turbines; i.e, all of the turbines of a row must be connected from the last one to the first one following a radial cable arrangement as in HR-I [37]. • The HV cable trajectory from the offshore substation to the coast transition is not optimized.
For this particular case of study, the following constraints apply: • A maximum of 80 turbines are allowed. • The turbine positions are limited to the shallow waters plotted in clear blue in Figure 1 (Note: the shallow waters at the right of the dashed line is considered as forbidden zone and hence disregarded). The search area in this study is limited to a reduced space, but in general the search area can be enlarged to thousands of km 2 , without compromising the computation time. A GA like the one used by most authors working on WFOP is not appropriate for handling continuous variables. In this sense, an NGEA operating with real values has been programmed to determine the eight decision variables that define the optimal parallelogramshaped OWF. It has basically the same structure as a typical GA, although adapting its crossover and mutation operators to work with continuous variables. The algorithm has been programmed following the operation sequence of Figure 7, where the term crossover is substituted by recombination, as usually preferred in NGEA to design the operator that combines the parents' genotype to generate new offspring. Its operative specifications are: size of the population, 150; maximum number of generations, 250; repetitions to finish, 40; individuals in the initial population, 100; crossing probability, 80%; mutation probability, 12%.
• Turbine cost and characteristics.
• Foundations. Depth charts and type of soil

Anticipation of the Optimal Orientation and Inclination Angle
Before launching the algorithm that performs a complete search for the optimal solution in terms of location, dimensions and layout, an exhaustive scan has been carried out to determine the optimal orientation and inclination for a given number of turbines and separation between them. This task is not mandatory, but it will provide the designer with certain hints about the validity of the complete algorithm solutions.
The objective is to obtain the energy lost due to the wake effect, assuming the following conditions:

•
The wind rose and the Weibull parameters are those specified for the HR-I site during the period 1999-2002, obtained from [36] and reproduced in Table 1. Data have been interpolated according to the method explained in Section 2.3 to estimate a continuous distribution of the density of probability and Weibull parameters. • The number of rows and turbines per row are n r = 10 and n tr = 8, respectively. • The separation between turbines of a row as well as between rows is d r = d tr = 7 D.

•
The turbines are Vestas V80, 70 m high tower. Power and thrust data were obtained from [38,39] and are shown in Table 2. Under these conditions, the energy deficit caused by the wake effect has been calculated for each of the possible orientations of the arrays between 0 and 180 deg, as well as for different angles of the parallelogram, between 45 and 135 deg, both with a discretization of 1 deg. The resulting AEP is shown in Figure 8. Maximum values of the produced energy are obtained for an orientation of θ = 65 deg, and an inclination of φ = 91 deg. The actual layout orientation is considerable different, θ = 173 deg, φ = 83 deg. Layouts producing 99.9% of the maximum energy yield are contoured in blue, highlighting that there is a large region of values for which the energy yield remains almost in its highest value. The right-hand plot shows the AEP dependence on θ for different values of φ. In fact, the maximum variation in AEP is 4 GWh, that represents 0.55% of the maximum yield, a reduction of 0.30 €/MWh in the LCOE, and an increase of 1.65 M€ in the NPV. This visualizes the importance of an adequate selection of these parameters. As deduced from the wind rose of Figure 5, the most frequent wind directions are centered at the WWN direction, with very frequent winds coming from West, for which the actual OWF layout is very sensitive in terms of wind speed deficit. Thus, higher OWF efficiencies are reasonably expected for a layout arranged to alleviate the wake effect for these wind directions (see central layout of Figure 9). The same evaluation was performed in other three situations, whose results are presented in Figure 10. The left-hand plot represents the results for a similar layout with half of the arrays and half of the columns (5 × 4). High similarity with the left-hand plot is observed, with a slight difference in the maximum location. This deviation is negligible because the zone contoured in blue (indicating 99.9% of the maximum energy) is very similar to that of the original 10 × 8 layout, and any solution inside this contour will give rise to the highest (or almost) values for AEP.

Continuous. Original 10x8
The central plot chart shows an analogous representation for a quincunx layout (see the right-hand layout of Figure 9), in which the even rows are displaced by d q tr /2, where d q tr is the distance between turbines (in a row) in this quincunx layout. In order to occupy the same concession area, d q tr /2 is obtained from the distance between turbines in a row, in diameters (d tr ), as d q tr = d tr * n tr − 1 n tr − 0.5 = 7 · 7/7.5.
Results are not better than the original parallelogram-shaped OWF, and it is not expected that using a quincunx layout will increase the energy yield if the distance between turbines is reduced according to (26).
Finally, the right-hand plot of Figure 10 corresponds to evaluations using the original sectorized data from sectors, without performing any interpolation. Although a logical similarity exists between this representation and that from Figure 8, there are significant differences in the maximum magnitude and location. Even when the actual values of the wind resources are unknown, the continuous distribution is more likely to behave as the real one, and the sectorized solution will give rise to erroneous results. Figure 10. Comparison of the energy produced in different scenarios. Contour lines for constant levels of energy deficit for a simplified 5 × 4 layout (left-hand plot), using the original discretized data (central plot) and for a quincunx layout (right-hand plot). The optimal values are pointed out with red asterisks, and zones with 99.9% of the maximum yield are contoured in blue.

Possibility to Limit the Speed Range for the Deficit Calculation
During the evaluation of the energy deficit of the possible layout, the complete set of speeds is evaluated. Figure 11 plots the contribution of every speed and wind direction to the final produced energy, evidencing that only wind speeds between 5 m/s and 15 m/s are of consideration for the deficit evaluation. At lower wind speeds, the extracted power is reduced. At the other end, at speeds higher than 15 m/s, C t is lower than 0.25 (see Table 2) and, from expression (8), the impact of wakes is reduced. In addition, the slope of the curve v-P is low at these high speeds, and moderate speed reductions will not provoke significant power deficits.  Figure 11. Contribution of different wind speeds and directions to the produced power and energy. Electric power produced by the OWF for every wind speed and wind direction (left-hand plot) and energy taking into account the density of probability (right-hand plot). Contour lines are added for a value = 1% of the maximum one.

Results of the Layout Optimization Search
The algorithm shown in Section 2.4 was performed ten times to determine the best layout for a wind farm in the site where HR-I is.
The first five optimization searches were made with the constraint of limiting the turbines to the clear area corresponding to shallow waters (<10 m) of Figure 1, excluding the extraction areas and sensitive zones close to the coast. Results are presented in Figure 12. As can be seen, there is no dispersion in the location of the OWF centre (left-hand plot). The orientation and inclination are also fairly constant (central plot). The right-hand plot reflects that solutions with very long distances are preferred, since they reduce the wake effect (the maximum distance between turbines was limited to 20 D). However, concession areas are usually reduced to smaller zones, and developers are forced to reduce the separation between turbines.
According to this concern, distances were fixed to 7D, a new battery of five searches was performed, and the results are presented in Figure 13. Again, the center is well located in the middle of the concession area. Orientation values vary between 55 deg and 95 deg, while the parallelogram angles are between 85 and 95 deg. The right-hand plot shows the most frequent layouts and the mean of LCOE for every layout. The layout of 16 × 5 is the most frequent one, with the lower (i.e., best) values for the objective function LCOE. Some solutions exist with 10 × 8, 11 × 7 and 13 × 6, these two last ones including less turbines than allowed by the concessionary. These solutions hardly appear at the end of the search process if NPV is chosen as the objective function. The evolution of the fitness at each generation is displayed in Figure 14, showing that the optimum is usually reached in less than 100 generations. The computation time is also presented, with an average of 2100 s per optimization search for long distances between turbines and 2700 s for normal separations.  Figure 15 represents the evolutions of the most frequent orientations for the rows and columns. The column orientation is obtained by subtracting the orientation and the inclination (see Figure 9). The preferred alignment directions are around 60 deg and 150 deg, trying to avoid the wind directions with more energy deficit (90 deg and 170 deg from Figure 11). Results for the five optimization searches are represented with different colors. The plot at the left represents the values when no separation restriction is imposed to the separation between rows and columns. In the right-hand plot, a separation of 7D is imposed.
Finally, Table 3 lists the decision variables of the final optimum for both types of searches and compares them with the one of the actual HR-I OWF, together with the breakdown of estimated costs.  [40]. Starting from these, the estimated cost for the actual layout is EUR 293.5 M, which means 1.83 M EUR/MW. The real project cost, obtained from [20], was EUR 278 M. The third column corresponds to the optimization result when distances between turbines are fixed and equal to 7D between rows, and between turbines in a row. Optimum orientations for rows and columns are different from the actual ones and similar to the central plot of Figure 9. However, the main difference lies in the number of rows, which in this case increases to 16 ). With regard to the economic indicators, the investment cost hardly varies from one solution to another. An exception is the significant reduction in the collection MV electrical infrastructure for the Optim 7D, due to grouping of sixteen turbines in the same row, for which the cable section can be optimized. For any of the solutions, the values for LCOE indicate that the price of the energy should be sufficiently higher than 60 EUR/MWh to make this project attractive to developers.

Conclusions and Discussion
There are numerous research works carried out on the energy losses suffered by the wake effect, especially relative to the HR-I OWF. Genetic algorithms also abound that obtain the layout that maximizes the return on investment, although none takes into account all of the factors that affect the profitability of the project. In addition, they give rise to configurations that are not desirable due to their lack of uniformity in the location of the turbines.
In this work the result of a complete optimization algorithm is presented: • that takes into account all of the factors (or at least more than any other paper found) required for the economic evaluation; • that fully defines the position, size and configuration of the optimal OWF with regularly distributed turbines and shaped like a parallelogram; • that, alternatively, indicates what should be the orientation and inclination of the concession area to make it more attractive to potential developers of wind farms.
With the fixed electricity production feed-in tariff provided for grid connections at that moment equal to 0.058 EUR/kWh [20,30], the project is not profitable except for very low interest rates. In the case of uncertain values for the energy price, the most suitable indicator to evaluate the project profitability is the LCOE.
From the exploration of the results obtained with a wide concession area, a tendency is observed for the algorithm to search for very dispersed OWFs, with little exploitation of the available surface. If, to reduce the occupied area, the distance between turbines is limited to, for example, seven diameters (as is the case with HR-I), a solution search algorithm with regular patterns offers the optimal orientation, inclination and layout. In this case, and assuming the rows and columns configurations of (4 × 20, 5 × 16, 8 × 10, 10 × 8, 16 × 5, 20 × 4), an exhaustive search only needs 180 × 90 × 6 90,000 evaluations to calculate the optimal configuration, so that the use of a random search algorithm is meaningless. With regard to the orientation and inclination, results show that optimal solutions tend to avoid the more frequent wind direction given by the wind rose.
However, if the algorithm also has to find the exact location, orientation, inclination and optimal separation distances between turbines (respecting a more or less wide concession area), the number of possibilities increases to the extent of requiring an efficient algorithm that can handle these variables. Under these requirements, a traditional GA is not appropriate. Instead, an NGEA that operates with real values has been used to determine the decision variables that define a parallelogram-shaped OWF. The inspection of the results shows that other regular layouts were possible with a better site exploitation.
In order to have a confident estimation of wind resources, a novel method of interpolating the data given for a reduced number of sectors is presented. Consequently, the discontinuous values describing the wind direction probability and Weibull parameters are interpolated along 360 deg using cubic splines, with the concern of maintaining the accumulated probability (in the case of the wind rose) or the average value (in the case of the Weibull parameters) at each sector. Obtained results show significant differences between the optimal orientation values when the algorithms uses continuous distributions for the wind resources and when it uses sectorized data.
Additional contributions from this project are: • discussion on the optimal orientation of the concession area for HR-I, • recommendation to reduce the range of speeds to be studied when evaluating the orientation and inclination of the OWF during the optimization algorithm execution, • possibility to deduce the optimal orientation and inclination replacing the original layout into a scaled version of the OWF, with half of the rows and half of the turbines per row, • the evaluation of a quincunx configuration, which does not entail an increase in the energy produced, mainly if the position of the turbines must be limited to a certain concession area.