A multi-objective optimisation approach applied to offshore wind farm location selection

This paper compares the three state-of-the-art algorithms when applied to a real-world case of the wind energy sector. Optimum locations are suggested for a wind farm by considering only Round 3 zones around the UK. The problem comprises of some of the most important techno-economic life cycle cost-related factors, which are modelled using the physical aspects of each wind farm location (i.e., the wind speed, distance from the ports, and water depth), the wind turbine size, and the number of turbines. The model is linked to NSGA II, NSGA III, and SPEA 2 algorithms, to conduct an optimisation search. The performance of these three algorithms is demonstrated and analysed, so as to assess their effectiveness in the investment decision-making process in the wind sector, more importantly, for Round 3 zones. The results are subject to the specifics of the underlying life cycle cost model.


Introduction
In the last few decades, a necessity to reduce carbon emissions has been raised after concerns of the global warming effect that causes rapid changes in the environment. In fact, electricity production was found to be responsible for at least 24% of the total greenhouse emitted gases, in 2013 (Lin and Chen 2013 led to further alternative ways to generate energy for more sustainable solutions. According to the 20-20-20 target on reducing carbon emissions and the new Climate Conference in Paris (COP 21) on keeping the global warming temperature below 2 • C, it is important to contribute to the renewable energy (RE) investment growth in the UK by making the investments more attractive, information-rich and less risky (BEC CREW 2015).
Wind energy is one of the fastest growing forms of RE in the UK; however, since structural material prices have significantly increased over the last years, it has a direct impact on larger scale wind projects, the overall cost of turbines, and their operational and maintenance costs (European Observation 2011; Lin and Chen 2013; Mytilinou et al. 2015). The UK technology roadmap highlights that the offshore wind costs need to be reduced to £100 per MWh by 2020 and a greater confidence over financial motivations is required (Department of Energy and Climate Change 2011). The location of a wind farm and the type of support structure have great impacts on the installation costs. The most important costs in an offshore wind farm can be found in HM Government (2013).
The Crown Estate released Round 3 leases and provided 9 new considerably larger zones than Rounds 1 and 2; offshore wind farm zones will include up to 32 GW of power capacity. The new leases encourage larger scale investment plans and bigger wind turbines. The new zones include locations further away from the shore and deeper waters which could be more challenging (Department of Energy and Climate Change 2011; Renewables First 2017; The Crown Estate 2010a, b, 2013).
Decision making for offshore wind energy investment is governed by a variety of criteria as can be found in (Kolios et al. 2010(Kolios et al. , 2014(Kolios et al. , 2016Lozano-Minguez et al. 2011;Mar- One of the most important decisions arising when starting a new investment is the selection of a suitable offshore location and always requires extended effort. A methodology is proposed to help the decision-making process at these first stages of a wind farm investment considering the Round 3 zoned in the UK. Through this methodology, the location based on different physical aspects is selected. The aim of this paper is to discover the optimum offshore Round 3 location based on financial costs and demonstrate the effectiveness of the underlying methods. The wind energy project costs associated with the design specifications of the site in conjunction with the turbine type and number of turbines should be numerically captured. This considers physical aspects unique to each offshore location, such as wind speed, water depth, and distance from appropriate construction ports, and will be modelled using life cycle cost (LCC) analysis. The best location should be discovered by considering the conflicting nature of the cost elements, so as to reduce overall cost at the early stages of a wind energy investment. Multi-objective optimisation algorithms will be coupled to the aforementioned model and used to reveal the interplay among the cost elements. To increase the effectiveness of the optimisation process and the diversity of the results, three different algorithms will be employed.
The contribution to knowledge is as follows. First, the combination of a newly developed prototype framework that includes economic modelling and optimisation process is assessed, so as to select the optimum offshore Round 3 location of a wind farm in the UK. Second, a set of non-dominated optimal solutions is suggested using the prototype framework. Both are expected to assist project developers and researchers at the first stages of the design of a wind farm.
The suggested methodology is depicted in Fig. 1. It comprises of individual modules, coupled with generic interfaces, so as to enable incremental development. Here, the physical aspects of each wind farm location (i.e., the wind speed, distance from the ports and water depth), the wind turbine size, and the number of turbines are considered, to manually select the optimum, economically efficient, and viable option (Shafiee et al. 2015).
The structure of this paper is as follows. A literature review in LCC analysis and the main phases will be explained. A description of the optimisation process and categories, and their benefits along with the most commonly used optimisation algorithms will be presented. Next, the methodology will follow. Finally, the results will be analysed and followed by a discussion. Future avenues will be drawn in the conclusions.

Multi-objective optimisation
In real-world cases, multiple and conflicting objectives have to be improved simultaneously and multi-objective optimisation techniques have to be used, for example, minimising time versus energy efficiency (Branke et al. 2008).
During MOO, the decision space, the hyperplane that combines all the decision variables, is searched by evaluating the constraints and objectives. At the end of the process, a set of solutions is obtained and they are at least as many as the considered objectives. This is frequently called the Pareto Optimal Set or Pareto Front, where it is not possible to improve any objective without compromising in any of the others (Adinolfi et al. 2015).
MOO assists decision makers in appreciating the trade-off among conflicting objectives, before selecting the optimum solution for implementation, while understanding the interplay among the considered objectives. Some MOO studies can be found in Table 1.

SPEA
Strength Pareto Evolutionary Algorithm (SPEA) is an evolutionary based algorithm and it is also an MOO algorithm. In other words, SPEA is an Evolutionary Multiple Objective algorithm or EMO. SPEA is closely related to other evolutionary algorithms such as the NSGA, Vector-Evaluated Genetic Algorithm (VEGA), and Pareto Archived Evolution Strategy (PAES). SPEA has two versions, i.e., SPEA and SPEA 2, an extension of the former. More extensions can be found under the name SPEA+ and iSPEA. The aim of this algorithm is to locate and maintain a collection of non-dominated solutions (Pareto front) by examining thoroughly the search area by following an evolutionary procedure.
SPEA discovers and maintains a set of Pareto optimal solutions. An evolutionary process is used to investigate through the search space. During the selection process, a utility function is used, where an assessment method for dominance is combined with a density estimator. An archive of the Pareto Front is kept separately from the population of potential solutions used in the evolutionary process, thus demonstrating a form of elitism (Brownlee 2011).
SPEA 2 was selected for this work due to its suitability for MOO problems (Brownlee 2011;Nalianda 2012). SPEA is appropriate for combinatorial and continuous function MOO problems. A binary representation can be used for continuous function optimisation problems along with classical genetic operators such as one-point crossover and point mutation. In SPEA 2, the size of the archive is commonly smaller than the size of the population. It is possible to implement optimisation in calculations of density (of the revealed solutions) and Pareto dominance. More can be found in Zitzler et al. (2001).

NSGA
NSGA stands for non-dominated sorting genetic algorithm and it is also an MOO algorithm and an EMO. Currently, there are three versions of the code: NSGA, NSGAII, and NSGAIII, among which the last two will be considered in this study. NSGAII has been employed on a number of optimisation problems, mostly with two objectives, whereas NSGAIII is expected to be more appropriate when the number of objectives increases (Yuan et al. 2014), as in this research.
For comparison purposes, the performance of all three selected optimisers is depicted in Fig. 2, where they were applied to a benchmark test function with three objectives. Ideally, the trade-off should be as much dense and wide as possible, so as to uniformly cover the performance of the objective space. In Fig. 2, the shape of a reference problem (which was not as complex as the developed LCC) is illustrated. The application of these algorithms in the optimum selection of the wind farm location, as a class of problems, has never been attempted before. As part of the contribution to knowledge and to demonstrate by evidence the effectiveness of the optimisers to deliver non-dominated solutions in real-world modelling, all three optimisers were trialled. The optimisers were considered to compare the quality of the solutions, irrespectively of their internal functions. All algorithms can yield satisfactory results and are tested in different sectors according to the literature. Frequently, using one of them is acceptable and has delivered satisfactory results. As expected, there are advantages and disadvantages for each method, which will be investigated further below. Hence, employing all three methods will highlight the differences among them, as certain algorithms behave better in certain problems (Wolpert and Macready 1997).
The class of NSGA algorithms was selected in this study, because it is suitable for MOO problems. NSGA is appropriate for continuous function MOO problems. A binary representation of a decision variable can be used along with classical genetic operators such as one-point crossover and point mutation. A real-valued representation is recommended for continuous function optimisation problems, which consequently requires specific genetic operators such as Simulated Binary Crossover (SBX) and polynomial mutation (Deb and Agrawal 1994;Jain and Deb 2014;Nalianda 2012). The number of divisions needs to be set in the NSGAIII algorithm. The divisions are a mechanism NSGAIII uses to control the spacing of a reference point as it progresses through the loops. It is suggested by Chiang (2014), Deb and Jain (2014) and Jain and Deb (2014) to use 12 divisions, where the domain of each of the many objectives is separated into regions using reference points to define normalised hyperplanes, for greater effectiveness in the discovery process. Both optimisers are considered as the state-of-the-art multi-objective algorithms, and they are selected for this study. In general, the optimisation problem can be solved through genetic algorithms, Tabu Search, Simulated Annealing, and other Heuristic methods (Momoh and Reddy 2014). A few methods that can also be applied to this study are the following: multi-objective evolutionary algorithm based on decomposition (MOEA/D), covariance matrix adaptation evolution strategy (CMAES), the third evolution step of generalized differential evolution (GDE3), indicator-based evolutionary algorithm (IBEA), optimized multi-objective particle swarm optimiser (OMOPSO), speed-constrained multi-objective PSO (SMPSO), etc. (Hadka 2015). However, the SPEA2 algorithm has not been considered yet for that type of applications.
Evolution and genetic algorithms are used throughout the literature in many energy-related sectors. Karimi presented an approach that links a multi-objective genetic algorithm to the design of a floating wind turbine. By varying nine design variables related to the structural characteristics of the support structure, multiple concepts of support structures were modelled and linked to the optimiser. The aim is to minimise the economic costs (as a combination of the mooring system, the anchor, and the offshore floating platform) and maximise the turbine performance (standard deviation of nacelle acceleration. The Pareto Front contained a wide number of results that reflect solutions of either the single-body platforms or tension-leg platforms or multi-body platforms. In the future steps, it was suggested to optimise the levelised cost of energy and to consider a different parameterisation scheme both of which could extend to the present research (Karimi et al. 2017).
Yan suggests a methodology that includes genetic algorithms and the analytical hierarchy process decision-making method in a study related to green suppliers (Yan 2009). The design of a new optimisation algorithm is proposed in Saavedra-Moreno et al. (2011), to optimise the layout of turbines in a wind farm. New aspects were considered for the objectives, such as the shape of the wind farm, a range of costs (expressed in a benefit for investment), and orography. The optimisation algorithm is based on evolutionary algorithms and is seeded by a greedy approach, where experimental comparisons have been demonstrated.
Five different types of optimisation algorithms were used in Elkinton et al. (2008) to optimise the layout of a wind farm. As it was also mentioned in Cagan et al. (2002), solutions of higher quality can be discovered using algorithms with stochastic elements. However, these operate at a much slower speed, than the deterministic counterparts. In fact, in Elkinton et al. (2008), combining genetic algorithms with heuristics was more effective and faster than using only one of them. In this particular instance, it was suggested to use layout optimisation in small areas, or, at least, focus in good areas before launching a bigger optimisation case.
Another approach, presented in Wan et al. (2012), was used to optimise the layout of a wind farm (micro-siting optimisation: choosing the type and location of wind turbines) by considering continuous space and using particle swarm optimisation techniques. A special local search scheme was also introduced in the optimisation algorithm, to successfully speed up the process, where realistic solutions were discovered that delivered more electricity. In Papatheou et al. (2015), an evolutionary optimisation algorithm is used in the area of supervisory control and data acquisition or SCADA in a Swedish offshore wind location called Lillgrund for monitoring purposes to optimise and predict the energy production from every wind turbine in a wind farm considering related data from other wind turbines. Finally, in González et al. (2010), evolutionary algorithms are applied in a wind farm-related optimisation problem. The configuration of the layout of the turbines is optimised based on a cost model. The problem is very complex, since the optimum layout has been troubling many specialists for years. The suitability of the suggested evolutionary techniques is proven in the study.

Life cycle cost modelling
Life cycle cost was established by the Federal Energy Management Program to evaluate economically energy, water conservation, and RE projects for federal facilities. A guide was created in a handbook in Fuller and Petersen (1996), where the methodology and criteria were evaluated and presented. LCC analysis can evaluate and suggest cost reductions throughout a project's life. The outcome of the analysis can provide useful information in an investment and can direct decision making from the initial stages of a new project.
Calculating the LCC of an offshore wind project involves five stages from the predevelopment to the decommissioning phase, and there is not a common universal reference point for wind projects. In Shafiee et al. (2015), a parametric whole life cost framework for an offshore wind farm and a cost breakdown structure is presented and analysed. LCC analysis is essential for the insurers, wind farm operators, and investors to ensure a cost-efficient long and profitable investment plan to produce power.
LCC analysis gains more ground over the years because of the larger scale in wind projects. For example, Nordahl (2011) studied the advantages and disadvantages of the transition to offshore wind and proposed an LCC model of an offshore wind development. However, the study mainly centres in a simplified model and especially the operation and maintenance stage of the LCC analysis and it is suggested that there can be a further full-scale LCC framework in the future.
A study that states legal, financial, etc. related problems of wind turbines that are positioned both onshore and offshore is described in Angelakoglou et al. (2014) and involves an assessment using the LCC analysis methodology to select the best option, environmentally, energywise, and economically. Martínez et al. (2009) offer further insight in LCC analysis for a wind turbine throughout its whole life beginning from the manufacturing and installation to the decommissioning. This aims to quantify the impact of each stage along with important aspects such as manufacturing, transportation to site, and material waste.
The importance of the LCC analysis and an economic life cost-related model in three different offshore floating wind devices was presented in Laura and Vicente (2014). This study aims to develop a framework and minimise some of the most important costs in a floating type of turbine making the floating devices a more attractive investment. In the same direction, Myhr et al. (2014) also studies five different types of floating wind turbine concepts in one offshore location and compares them using LCC-related features. The water depth was found to have the most significant impact on the total cost.
According to Shafiee et al. (2015), LCC analysis includes costs calculated from several stages of the wind project such as the predevelopment and consenting (C P&C ), production and acquisition (C P&A ), installation and commissioning (C I&C ), operation and maintenance (C O&M ), and decommissioning and disposal (C D&D ) stage. Since foundations and support structures moved further towards deeper waters, Round 3 locations moved further away from the shore, and larger scale wind turbines are now becoming more common along with transportation and delays especially through vessels, the cost has been increased considerably. This study follows the LCC framework that was developed in Shafiee et al. (2015) and investigates different Round 3 offshore locations using it as common ground to select the most costefficient one. Through the study, the main aspects that impact the final costs are discussed.

Selecting an offshore location
To the author's knowledge, there are no studies that combine the concept of LCC with optimisation techniques with a focus on the individual LCC costs, to find the optimum offshore location for wind farm projects. Moreover, there are no studies that consider objectives based on economic figures and select the optimum Round 3 offshore location in the UK. In fact, for the selection of the location, there is very limited work accessible and with a small amount of focused and related criteria on this topic.
In the literature, many RE and location selection studies can be found, but the findings and the formulation of the problems provided follow a different direction. For example, using goal programming, the offshore location for a wind farm installation was selected in Jones and Wall (2015). The study involves Round 3 locations in the UK, while shows its flexibility to combine decision-making methods. This work shows the energy production, costs, and multi-criteria nature of the problem also considering some important factors related to environmental, social, technical, and economic aspects. The LCC analysis and their formulation into an MOO problem were not employed in the application. According to Jones and Wall (2015), multi-objective modelling techniques for both onshore and offshore wind farm are quite underdeveloped. The present study focuses on a methodology to fill this gap by linking MOO with LCC as objective functions and compares optimisation algorithms to select the optimum solution.
Another study on offshore locations is provided in Cradden et al. (2016) for an RE platform using multiple criteria and geographical information systems (GIS). A range of problems that exist around offshore RE platforms have been reviewed and a combination of criteria has been selected for the Atlantic facing shores in Europe. The potential risks were studied and it was found that the extreme weather conditions show the necessity of a compromise between the designing costs and the extra energy production. Very important factors were also the lack of ports with suitable available construction infrastructure that results in under-exploited sites, access problems, and weather window conditions, even during the summer months (Cradden et al. 2016). Although the study is very thorough, it is mostly focused on environmental, geographical, and weather issues that are out of scope because of the economic nature of the objectives and related criteria.

A framework for the optimisation of deployment sites for Round wind farms in the UK
In Espinosa (2014) and Shafiee et al. (2015), a whole LCC formulation is provided and this study follows the steps and phases of the analysis for the optimisation problem. Assump-tions and useful data in the modelling of the problem can be found in references 4COffshore (2017c), Dicorato et al. (2011), Espinosa (2014, Laura and Vicente (2014), Shafiee et al. (2015), The Crown Estate (2017) and Wind Energy The Facts (2017). Based on the previous work, a new model was developed so as to be coupled with the optimisation algorithms and drive the optimisation search. The framework described in this section assesses the effectiveness of a suggested methodology to discover the optimum Round 3 offshore locations in the UK and improve the decision-making processes. Conceptually, the framework comprises of a model and an optimisation algorithm. The flow chart in Fig. 1 shows the optimisation model that includes seven objectives, four LCC-related objectives, described in Shafiee et al. (2015), and three additional objectives. Optimising seven objective functions at the same time, which are conflicting (from the mathematical formulation below), classifies the problem as multi-objective and it is considered rather complicated because of the interplay of the objectives and the nature of the variables.
The LCC model described in Shafiee et al. (2015) is used in this study and provided below in detail.
The LCC is calculated as follows: Life cycle cost C P&C : Predevelopment and consenting cost C P&A : Production and acquisition cost C I&C : Installation and commissioning C O&M : Operation and maintenance cost C D&D : Decommissioning and disposal C P&C = C projM + C legal + C surveys + C eng + C contingency C projM : Project management cost C legal : Legal cost C surveys : Survey cost C eng : Engineering activities cost C contingency : Contingency cost C P&A = C WT + C SS + C PTS + C monitoring C WT : Wind turbine procurement cost C SS : Support structure/foundation cost C PTS : Electricity transmission cost C monitoring : Monitoring cost C I&C = C I&C-port + C I&C-comp + C comm + C I&C-ins C I&C-port : Port-related cost More can be found in Shafiee et al. (2015). The present study only considers the first four life cycle costs (depicted in Fig. 3) in the proposed methodology, because the aim is to understand the interplay between the CAPEX and OPEX costs to improve the decision-making process. Both these costs are considered to drive investment decisions, so as to prevent potential risks and issues when beginning a wind project by the developers.
As it is depicted in Fig. 1, the first four objectives are the costs from the LCC analysis. More specifically, the present model includes the predevelopment and consenting, production and acquisition, installation and commissioning, and finally operation and maintenance costs. The decommissioning and disposal cost is not considered at this stage. All of the cost-related objectives are minimised, as shown in Fig. 1. The mapping between the variables and the objectives that are estimated using the LCC analysis is depicted in Fig. 4. In this representation, the number of turbines, the distance from the ports, the water depth, and the power rate, which is determined by the wind turbine size shown in the specifications in Table 3, are the decision variables, as shown in Fig. 1.
The last three objectives are the number of turbines (NWT), the power that is extracted (P) from each offshore site and the total installed capacity (TIC), which are minimised, maximised, and maximised, respectively. The power extracted is calculated by the specific mean annual wind speed of each location along with the characteristics of each wind turbine both of which are considered inputs (listed in Table 2). The TIC is calculated by the number of turbines and the rated power of each of them.
The power extracted in this optimisation model is maximised and it is calculated for each site and each wind turbine, respectively, from P = 1 2 AC p ρu 3 A: Area of the wind turbine ρ: Air density C p : Power coefficient  (2017). The list contains designated, appropriate, and sufficient construction ports that are suitable for the installation, manufacturing, and maintenance works for wind farms. New ports are agreed to be built for the conveniences of new wind farms. This study will consider only the parts and sites depicted in Fig. 5 For the distances from the ports calculation, QGIS was used. QGIS is an Open Source licensed Geographic Information System (GIS), which is a part of the Open Source Geospatial Foundation (OSGeo) (QGIS). The distances were calculated by the assumption of the nearest port to the individual wind farm, in a straight line. The specifications in Table 2 were acquired from 4COffshore (2017c); for each offshore location, a special profile was created including the coordinated, the distances to the shore and port, annual wind speed, and average site water depth. Table 2 shows among various data, the locations that each of these zones contains. Each location correlates with their specific data used in this problem. Table 3 lists the specifications of the turbines that are considered in this study. These were extracted from reference turbines in 4COffshore (2017d).
The optimisation problem formulates as follows: where the objectives are described above. The site index and turbine-type index are specified in Tables 2 and 3, respectively. In the formulation, minimising the negative TIC and P is equivalent to their maximisation. The number of turbines was deliberately selected both as variable and as objective so as to minimise the CAPEX.
The optimisation modelling has been completed using the library platypus in python (Hadka 2015) and the selected optimisation algorithms are NSGA II, NSGA III with 12 divisions, and SPEA 2. For all the algorithms, the default implementations were employed and the stopping criteria were set to 10,000 iterations.  . 6 Scatter plot matrix among the seven objectives, i.e., C A&P , C P&C , C I&C , C O&M , NWT, P, and TIC using NSGAII In this study, a comparison among the pareto front solutions from each optimiser is conducted, where each optimiser delivered a non-dominated trade-off with the best possible solutions for offshore wind farm locations.

Results
The pareto front solutions from every algorithm and their trade-offs are presented. Although the variables are discrete in nature, the revealed trade-offs appear to be continuous. Each pareto front solution consists of the wind farm location, number of turbines (NWT), type of turbine, power extracted (P) for the specific site, total installed capacity (TIC) (i.e., the total capacity of the wind farm), and the life cycle costs (i.e., C P&A , C P&C , C I&C , and C O&M ).
A scatter plot matrix was considered as the best option to visualise a seven-dimensional problem and it is depicted from Figs. 6,7,8,9,10,and 11. In Figs. 6,8, and 10, seven objectives, i.e., C P&A , C P&C , C I&C , C O&M , NWT,P,and TIC,are illustrated,and in Figs. 7,9,and 11, the CAPEX, OPEX, NWT, P, and TIC are also provided. This is part of the Pareto Front, as it was discovered by the optimisers. In the main diagonal of the figure, the histograms represent the concentration of points in ten buckets of equal size. All the trade-offs are continuous, which means that there is not any (discontinuous) gap in performance.

NSGA II
In Figs. 6 and 7, the revealed trade-offs demonstrate that non-dominated solutions were discovered. However, certain areas were not explored thoroughly. For instance, the plot of 'P versus NWT' illustrates that there is a strong upper limit, as expected from the nature of the model. In addition, in the same plot, it is obvious that the majority of solutions are Fig. 7 Scatter plot matrix among CAPEX, OPEX, NWT, P, and TIC using NSGAII focused on the upper bound, whereas below that, the number of points explored is scarce. This practically means that the optimiser quickly identified chromosomes with good performance and used the computational budget towards exploring along this direction. This pattern applies to most of the plots. Plots whose performance has a pattern similar to 'C O&M versus C I&C ' also demonstrate that certain performance routes were discovered and followed throughout the optimisation search. This explains the solid legs within the swarm of solutions. The bottom area of 'P versus C P&A ' illustrates that a few areas of dominated performance (sporadic solutions) were discovered, but were not further investigated. Finally, 'C P&C versus C P&A ', 'TIC versus C I&C ', and 'C O&M versus NWT' vary in harmony.
By interactively investigating the results, the following behaviour was identified. Each of the legs of the plot 'TIC versus C O&M ' is linked to a particular wind turbine. Then, by varying the combination of the number of turbines and the site, a wide range of values of the cost element can be obtained. This is more obvious in the PFs discovered by NSGAIII because of the continuous points discovered for the same type of plot, whereas this would be very hard to spot in the results from NSGAII. This feature could be integrated into a process to detect any hidden performance relationships.

NSGA III
Compared to NSGAII, NSGAIII was also found to behave better when the number of objectives increases (relative to the other two optimisers), as expected by definition. Hence, the final results appear more uniform, as shown in Figs. 8 and 9. Consequently, the trade-offs for certain pairs of objectives are much more complete (in terms of the distance between any two points) and richer (in terms of a number of points within a relatively narrow area). It seems that the chromosomes covered a wider spectrum of solutions. In NSGA III, the 'NWT versus C O&M ' present the most straightforward relationship as they are approaching an exponential trend. As before, the same three combinations clearly vary in harmony. In general, this demonstrates the suitability of NSGAIII in such problems (and will be considered as the main optimiser in the next stages of the research).

SPEA2
The performance of SPEA2 is depicted in Figs. 10 and 11. Fundamentally, the same characteristics can be observed. However, SPEA2 stands between NSGAII and NSGAIII, in terms of trade-off findings. It discovered more diverse solutions than NSGAII but less diverse than NSGAIII.

Discussion
Optimum locations for a wind farm have been discovered using three different MOO algorithms, NSGA II, NSGA III, and SPEA 2 and using the LCC analysis, to achieve costefficient solutions. The results follow a consistent trend and they seem to be in relative agreement. TIC, C P&A , C I&C , and C P&A vary in harmony, as shown in Figs. 6,8,and 10,which Fig. 8 Scatter plot matrix among the seven objectives, i.e., C A&P , C P&C , C I&C , C O&M , NWT, P, and TIC using NSGAIII Fig. 9 Scatter plot matrix among CAPEX, OPEX, NWT, P, and TIC using NSGAIII suggests to be compound into a single objective in the future, for simplicity. In addition, C O&M and NWT also vary in harmony by following a parabolic trend. In the same figures, most of the discovered solutions live at the lower end of the range of C P&A and C P&C . This points out future research directions, so as to identify areas, where the performance remains constant.
In the histogram of C O&M , many points of the optimal revealed behaviour reside at the lower bound of the range of the objective. Consequently, this is also noted in the concentration of points in the CAPEX and OPEX. This was the easiest to discover relative to the other objectives. Because the shape of the trade-offs among C P&A , C P&C , and C I&C is similar and the same holds of the trends of the histograms,  it is sensible to consider the CAPEX objective in the future steps of the research (instead of all these three individual objectives).
In SPEA2 and NSGAII results, the sparse dots at the end of each graph represent that the method (i.e., the search patterns employed by those algorithms) has not discovered a wide range of solutions in those areas.
Several results show a relative harmony between each other and some others show the conflicting performance between the objectives. Six offshore locations scored higher than 10% in the frequency graph, i.e., the Seagreen Alpha and Bravo (in Firth of Forth), Teesside C and D (in Dogger Bank), Rampion (Hastings), and the Celtic Array South West Potential development Area, which represent the Irish Sea (Celtic Array). A frequency graph has been created in Fig.  12 to summarise and depict the % frequency of each offshore location that appears in the pareto front solutions. In general, the Firth of Forth, Dogger Bank, and Irish Sea (Celtic Array) zones reached the highest score for all optimisers.
The Teesside C scored first and it is closely followed by Seagreen Bravo when using NSGAII. The Seagreen Alpha in Firth of Forth seems to have the highest scores of all when using the SPEA2 algorithm, and thus, it is considered one of the best options for a wind farm project, as it was identified by all three optimisers. Finally, Seagreen Alpha was scored first for the NSGAIII. In general, the Hornsea Projects family has been found to be suitable in only a small amount of appearances in the pareto front.
The performance of the algorithms has shown that NSGAIII demonstrated its suitability in multi-objective problems as its results appear to be more uniform and clear because of its main design, compared to the other optimisers. Therefore, the trade-offs for certain pairs of objectives are more complete, wider, and richer, in terms of a number of points. In general, the patterns of the revealed trade-off are very clear and distinct for the results from NSGAIII with 12 divisions.
The non-dominated results demonstrate that employing MOO algorithms was a sensible choice, so as to complement the process of wind farm location selection. For example, both Moray Firth Western Development Area and Seagreen Alpha were found to be in the trade-off at least once by the optimisers. For the developer, this means that it is equally cost efficient to choose either location. However, the latter has appeared significantly more frequent than the former. The developer could accordingly allocate the development  budget in different development phases, as required. More specifically, the Seagreen Alpha provides many more optimum options that are equally efficient (involving the number and type of turbines in each solution) than the Moray Firth Western Development Area, which gives the flexibility to invest more money in the installation or the maintenance stage of the project. An application has been submitted and consented for the Firth of Forth Seagreen Alpha and Bravo, which are at the first phase of the Round 3 Firth of Forth development according to 4COffshore (2017f), RYA (2017). Both Seagreen Alpha and Bravo faced some engineering and environmental problems and developed the project accordingly to accommodate such constraints that appeared in their offshore locations (4COffshore 2017f). Teesside C and D in Dogger Bank applications were submitted together by Forewind. Unfortunately, the projects have been cancelled according to 4COffshore (2017a4COffshore ( , 2017b. Rampion (Hastings) is currently under construction according to 4COffshore (2017e). Finally, Celtic Array South West Potential development Area was also cancelled (4COffshore 2014).
Probably, the six primary selected offshore locations have been either cancelled, consented or just submitted and only one of them is under construction, because other factors are also involved in the selection process. For example, the cancellation in all cases happened, because the Crown Estate asked the developers to revise all the terms of their agreement, while the project transitioned toward the development phase (4COffshore 2017e). Other reasons could also result in cancellation such as environmental problems, legal and procurement restrictions, etc.
By relating some of the most important techno-economic LCC factors to the physical aspects of each wind location (i.e., the wind speed, distance from the ports, and water depth), the wind turbine siz,e and the number of turbines, it is possible to discover a more cost-efficient solution. By comparing the outcomes of the three algorithms and suggesting the most suitable locations, useful insights are provided for both industrial and educational purposes in the wind sector for future investments.
In Fig. 13, the performance of each optimisation algorithm in CAPEX versus OPEX is shown. As expected, most of the solutions gather at the lower left end. NSGAII and SPEA2 discovered and followed certain trends for a particular wind turbine, as shown by the various legs in the figure. Because of the design of NSGAIII, the discovered solutions are more spread throughout the objective space. Table 4 shows the top three locations discovered by each optimiser in terms of frequency in the Pareto Front from Fig. 12. The listed solutions were selected manually so as to demonstrate the conflicting nature of the objectives. The table numerically demonstrates the conflict among the objectives in the trade-off, for instance, comparing Seagreen Alpha and Seagreen Bravo using SPEA 2 algorithm, when the NWT increases then P decreases. In the results of the same algorithm, including Rampion (Hastings), when the NWT increases, TIC reduces (whereas previously, it increased). This demonstrates the conflicting nature of the results and non-linear relationships, which will be further investigated in the future.
From the obtained results, the average savings were calculated to assess the performance of the approach and the framework. An average savings formula was used (i.e., Average savings = Maximum cost−Average cost Maximum cost ) to calculate percentages for CAPEX and OPEX using NSGA II, NSGA III, and SPEA 2. NSGA II reached 77 and 66.7% savings for CAPEX and OPEX, respectively. NSGA III found 69.1 and 59.2% for CAPEX and OPEX, respectively. Finally, SPEA 2 shows 70.1 and 55.8% savings for CAPEX and OPEX, respectively. It appears that the optimisation approach can discover solutions with lower CAPEX deviation than OPEX. The high percentages were expected because of the large number of conflicting objectives which forces the optimisers to discover a great number of solutions in the Pareto Front.

Conclusions
This study successfully demonstrated by example the effectiveness of the newly developed optimisation process and delivered satisfactory outcomes for the most suitable and cost-efficient offshore wind farm Round 3 locations. A comparison has been presented among the three state-of-the-art algorithms (i.e., NSGA II, NSGA III, and SPEA 2), which were applied to a wind energy real-world case. The comparison and the useful outcomes on their performance have been illustrated and discussed. The optimum locations for a wind farm have been suggested by considering the significant input of the LCC analysis. Six sites were suggested (with the frequency of appearance higher than 10% in the parent front).
The results follow a relatively similar and consistent trend. The performance of the algorithms has shown that the NSGAIII demonstrated its suitability in multi-objective problems as its results appear to be more uniform and clear because of its design compared to the other optimisers.
The limitations of this work are related to the LCC calculations and the associated assumptions that had to be made in the development of the model. To get more accurate results, more precise data are required to validate and calibrate the LCC, which could refine the results of this work. Many data and assumptions taken for this study have been obtained and chosen from wind-related databases and crown estate reports, to reach the real-world values. However, research and surveys are important for every individual site to have more accurate inputs. Data acquisition is the hardest part as it is impossible for a developer to proceed without a project plan. Here, only one type of foundation was considered, the jacket structure. The impact of Net Present Value on the economic objectives (i.e., CAPEX and OPEX) has not been considered and will be further investigated in the future.
The revealed outcomes will have an important impact on a possible extension of the Round 3 zones in the future of the UK and will help decision makers for their next cost efficient investment move. The proposed framework could also be applied to other sectors to increase investment confidence and provide optimum solutions. For example, the installation of floating offshore wind and wave devices could be benefited by the framework, where the optimum locations can be suggested according to cost and operational aspects for each technological need.