Layout Optimization Process to Minimize the Cost of Energy of an Offshore Floating Hybrid Wind–Wave Farm

Offshore floating hybrid wind and wave energy is a young technology yet to be scaled up. A way to reduce the total costs of the energy production process in order to ensure competitiveness in the sustainable energy market is to maximize the farm’s efficiency. To do so, an energy generation and costs calculation model was developed with the objective of minimizing the technology’s Levelized Cost of Energy (LCOE) of the P80 hybrid wind-wave concept, designed by the company Floating Power Plant A/S. A Particle Swarm Optimization (PSO) algorithm was then implemented on top of other technical and decision-making processes, taking as decision variables the layout, the offshore substation position, and the export cable choice. The process was applied off the west coast of Ireland in a site of interest for the company, and after a quantitative and qualitative optimization process, a minimized LCOE was obtained. It was then found that lower costs of ~73% can be reached in the short-term, and the room for improvement in the structure’s design and materials was highlighted, with an LCOE reduction potential of up to 32%. The model serves usefully as a preliminary analysis. However, the uncertainty estimate of 11% indicates that further site-specific studies and measurements are essential.


Introduction
One of the main conditioning issues of the planet's climate since the industrial revolution is global warming, a consequence of greenhouse gases in the atmosphere. More than 97% of the scientific community agree on the fact that human emissions are the main reason of this phenomenon [1]. Consequently, since the late 1990s, several efforts have been launched at the international level. Among other targets, middle-term national green new deals point at a transition to a full sustainable energy generation system as it is one of the main sources of these greenhouse gas emissions, representing around 25% of them according to the Intergovernmental Panel on Climate Change (IPCC) [2]. equations involved in the optimization process, focus of this study, and the constraints associated to it are presented. Last, the used optimization algorithm is concisely explained.
As shown in Figure 1, the tool is composed of four main modules: the wind and waves power generation models, the costs calculation model and finally, the overall optimization algorithm that modifies the decision variables, thus improving the output of the model in an iterative way. The waves and wind power production modules aim at simulating the physical and technical phenomena that result in a certain Annual Energy Production (AEP), including the interaction between platforms (wind wakes and waves shadows). This is the source of most of the nonlinearities of the model. Their detailed description can be found in [17] together with an outline of FPP's P80 technology. The tool presented in the diagram of Figure 1 consists of a lifetime energy generation and costs calculation of an offshore hybrid farm project. The farm is composed of several P80 platforms that are connected to an offshore floating electrical substation through inter-array cables; from the offshore substation, an export cable transmits the electricity produced in the farm to an inland electrical substation which is connected to the grid; furthermore, a harbor is used to carry out the marine operations in the farm.

The Costs Model
The costs model is a series of economical calculations, mostly linear, that use these AEP results as inputs to find the costs of the project in relation to the energy production by means of the LCOE indicator. The costs model includes the costs due to the capital expenses (CAPEX); the costs derived from the operational expenditures (OPEX); and those associated to the energy losses due to the generators' downtime periods and transmission of the electricity. The inputs to the costs model are the gross AEP from the wind and the waves after the wakes and shadow losses, the values related to the site characterization (bathymetry and coordinates, and the inland substation, harbor, and shore positions), and more than 140 costs related figures. The output is the LCOE, introduced below.

Levelized Cost of Energy (LCOE)
The LCOE, expressed in A C/MWh, indicates the average cost of the energy produced throughout the lifetime of a project, and therefore the minimum average price at which the energy must be sold to make the project feasible [16]. Unlike other simpler indicators, the LCOE takes into account the time value of money, as it weights the annual expenses TC t in each year t by means of a discount rate r, which is determined by the investors according to the inflation and the estimated risk of the project.
It is arguably the most useful indicator for seeking financing as it is ideal to compare the costs of different technologies [18]. The LCOE is calculated in the model as follows, where q t represents the yearly energy production in period t (obtained in the wind and waves models), and LT the lifetime of the project. The annual expenses included in TC t are presented below. As stated above, in this paper the LCOE values are expressed as dimensionless indexed numbers where the relative change is well represented.

Capital Expenditures (CAPEX)
The capital expenditures of a project involve all the investment, administrative, acquisition and installation related costs before and along its lifetime. They are calculated as follows, The categories considered in the model are design and project management expenses (C PM ), wind turbine acquisition cost (C WTG ), P80 costs (C P80 ), mooring costs (C Moor ), costs associated to the electrical infrastructure (C EI ), installation costs (C Inst ), and decommissioning costs (C Decom ); all of them expressed in A C.
• The "design and project management" costs are applied in year zero and they are divided into concession (C concession ) and design (C design ) costs: C concession = C con., cap · Cap Tot (4) • The acquisition cost of the wind turbines, one of the main cost components of a wind farm project, is calculated as a linear function of the wind turbine capacity: where is the total capacity of the farm. C WT cap. [A C/MW] is the considered cost per capacity.
• The acquisition costs of the P80 platform include the costs associated to the platform's structure and the waves energy converter (absorbers and power take-off system). They are calculated as follows, where C plat [A C/P80] is the cost of a P80 platform. n plat is the number of platforms.
are costs associated to condition monitoring and the SCADA systems.
• The mooring costs (C Moor , in A C) are divided into the anchor acquisition costs, which include the six anchor kits needed in each P80 unit, and the mooring line costs, which are considered to increase linearly with depth: depth under platform i.
• The investment costs associated to the electrical infrastructure are composed of those associated to the inter-array cables (C MV ), the export high voltage cable (C HV ), the offshore substation (C sub. ), and some additional costs associated to other devices (C el. other ): The cables costs are modeled as a linear function of their length, the additional costs depend on the number of platforms, and the substation's costs are calculated as follows [19]: Tot (10) • The installation costs are those expenses associated to both the human workforce and the technology used in the transportation, assembly, connection, and testing of all the infrastructure involved in the construction of the farm up to its commissioning: where the presented elements describe the various installation processes involved at the beginning of the project. They correspond to the following.
is the cost of transportation of the components to the assembly site.
assembly costs.
costs of installation of the electrical infrastructure.
costs of on-site platforms installation.
commissioning or connection costs.
Most of these costs, not detailed here for the sake of brevity, but specified in [17], are fixed or of linear nature depending on variables such as the number of elements (e.g., platforms), length of cable used, cable laying speed, cruise speed of cable lying, and so on, multiplied by the unitary cost of those items.
• Finally, the decommissioning costs are calculated. They are applied in the last year of the project's lifetime. The costs due to the removal of all the installed elements may be expressed as representing the various decommissioning categories, namely, costs related to Once again, most of these costs are of linear nature.

Operational Expenditures (OPEX)
The Operational Expenditures include the costs associated to all types of maintenance and any other operation along the lifetime of a project (e.g.,: inspections, checks, testing, cleaning, painting, replacement of components, etc.), being often divided into planned and unplanned or unscheduled O&M [20]. These maintenance operations have an associated energy generation downtime, as they often involve a failure or require shutting down the generators. The generation downtime is thus modeled together with the O&M costs as part of the OPEX.
As done when modeling the CAPEX, and more relevantly in this case, due to its distribution along the years, the annual costs associated to the OPEX are calculated individually every year, so they can be discounted in the LCOE calculation as explained above. Thus, for each year, i, in the project's lifetime, the maintenance operations (O&M i , in A C) and their associated downtime (DT i , in hours) are calculated as follows.
• Planned O&M and Downtime: there are different types of planned O&M operations: yearly, intermediate, and major operations; each of which has an associated frequency, cost and downtime period. Depending on the planned operations in each year i, the associated costs and downtime periods are expressed as follows, where, • Unplanned O&M and Downtime. In the case of the operations due to unexpected or sudden breakdowns, the costs and their associated downtimes are implemented by sorting the failures into five different severity levels: remote reset, minor fault/ inspection, medium fault (module replacement), large fault (offshore repair), and major fault (harbor repair). Each one of these categories presents a certain occurrence probability that is translated in the model as a number of incidents per year. The yearly unplanned O&M costs and downtime are thus implemented in the model as shown below,

Optimization
Given the wide range of variables involved in the model, the nature of many of the parameters (some of which are discrete and cannot be easily described by explicit functions and/or involve specific processes for their calculation), the spatial constraints, and the lack of convexity (the presence of several minima), it is not feasible to carry out a classical analytical or gradient-based optimization [11,12]. For this reason, it is necessary to use meta-heuristic methods instead, for example evolutionary algorithms. They are calculation methods that take metaphorical inspiration on certain processes that occur in nature [10]. They work by improving the obtained results in an iterative way, modifying the decision variables in each iteration according to the algorithm's working principle or logic.
The Genetic Algorithm (GA) and the Particle Swarm Optimization (PSO) are among the most used methods. The former, extensively used in the wind energy field [11,13,14], emulates the principles of evolution theory by combining and mutating offspring solutions that give birth to new ones. With regard to the PSO, it is based on bird flocking, fish schooling, or swarming theory, also having strong ties with the GA [15]. Additionally, in the optimization field, and taking into consideration the description of the problem, both single and multi-objective algorithms may be used. As an example of a multi-objective algorithm, the Agent Swarm Optimization (ASO), originated in the water supply field, shall be mentioned. ASO integrates behaviors of the two described algorithms, among others, preventing it from terminating in local minima, thus increasing the efficiency and, as a result, the total time devoted to solve the problem [21,22].
In this work, as a first approximation, it was used the classical PSO in a single objective optimization setting. The algorithm was set up to obtain the minimum LCOE in any site with its specific inputs. The coding was carried out in MATLAB with inputs provided by Excel data sheets following three objectives: first, minimizing as much as possible the many uncertainties involved in this complex calculation, naturally still assuming some limitations; second, reducing its running time, especially for the sake of the optimization process; and, third, developing it in a user-friendly manner.

Particle Swarm Optimization
The PSO is primarily based on the movement of the different individuals (in this case known as particles, each of which represents a set of decision variables) with respect to the rest of the members (population) of the swarm [15]. The new variables used in each iteration shift from their predecessors according to three different components [23]: the current (time k) position of a particle i, namely, x i (k); the best position ever achieved by the particle, x p best i ; and the global best position, x n best , achieved by any of the particles in the swarm (or neighborhood). The current position of the particle is updated through addition of its velocity, calculated as In this expression, ω is the inertia parameter, c 1 the cognitive/own-best position/self-adjustment weight, and c 2 the neighborhood/social/group adjustment weight, the factors that influence the particle's motion; r 1 and r 2 are random numbers within [0, 1] to represent nature's inexact condition.
Four main parties are involved in the way the optimization process is built: the objective function (minimization of the LCOE, already described in Section 2.1.1), the decision variables (subject of change to improve the result), the decision variables' associated constraints, and the penalties (increase of the objective value due to the infringement of the constraints).

Decision Variables
Among the wide range of parameters that could be subject of change to minimize the LCOE, in the last version of the optimization algorithm eight parameters, divided in three groups, were chosen as decision variables, as described below.
• First, the decision variables associated to the positions of the platforms: arguably, the most relevant decision variables of the model. With a set number of platforms in the farm, different ways of arranging the turbines were considered. At first, the platforms were allowed to be located independently within the limited site's area. However, rather messy and unorganized arrangements were obtained, and this option was soon discarded due to consequent negative visual impact and the operational difficulties that this may involve (e.g., definition of boat routes and safe access areas). Moreover, this arrangement method involved having 2 decision variables per platform (x and y coordinates), which increased the computational time and decreased the efficiency of the algorithm significantly. To solve these issues, the platforms were finally arranged in rows and it was resolved that a total of five decision variables were necessary to define completely the positions of all the platforms: these are the number of rows (d. var 1 ); the separation between rows (d. var 2 ); the separation between platforms in each row (d. var 3 ); the row offset (d. var 4 ), i.e., the relative positional shift along the rows' direction of consecutive rows; and the orientation angle of the rows (d. var 5 ). As it can be seen in Figure 2, the platforms are located around the site's center coordinates, so the rows are distributed equally to each side. These input parameters or decision variables are then converted into the platforms coordinates inside the function: first, the number of platforms and rows is read, and therefore the rows are created with their specified spacing (d. var 2 and d. var 3 ); then, the uneven extra platforms are successively added from the last row towards the first one; finally, the orientation angle is applied by multiplying each platform's coordinates by a a two-dimensional rotation matrix (R θ ). • The second group of decision variables (d. var 6 and d. var 7 ) are the coordinates of the offshore substation, specified as x and y distances with respect to the center of the site O. • Finally, the last decision variable (d. var 8 ) is the farm's export cable choice among the available ones, which has an impact in the costs and transmission losses as the chosen cable has its own specifications.
A number of other parameters that a priori seemed subject of becoming decision variables were eventually disregarded as such due to various reasons. These parameters are the following.

•
Harbor and inland substation choice among the available ones in the area: even though the LCOE output is definitely influenced by this decision, there are many external factors that condition it as this choice highly depends on the state legislation and planning; it also depends on particular characteristics of both the harbour and substation: their size, availability, grid strength, and so on. Therefore, it would be unrealistic to model these factors so it was considered that it should be the user who makes this decision, after a site-specific study.

•
The number of platforms: similarly, this factor depends mainly on the region and country where the farm would be installed. Factors like the strength of the grid or the need of power in the surrounding area will determine the total capacity of the farm, and thus it is something that must be chosen manually after studying these particularities.

•
WTG size: One of the initial decision variables considered. However, as the P80 is currently designed to host a 8 MW WTG, it was finally decided to develop the model with it. Another type of WTG would imply different input costs, a different hub height and therefore different wind speed inputs, and probably also a different scale of the platforms and WECs.

•
Lifetime: This parameter depends entirely on the design of the structures involved in the project, not on the project planning itself. As the P80 design is given for granted in the model with its consequent costs and characteristics, optimizing the lifetime is out of the scope of the model. In fact, the model would tend to increase the lifetime indefinitely in order to share the unrealistically constant investment costs with a higher amount of energy production along the years, thus reducing the LCOE as the lifetime increases. It was also considered to model an aging rate, but this was regarded as being too uncertain due to the lack of information about the real aging behavior of the P80.

Constraints and Penalties
All the decision variables have constraints that are directly or indirectly associated to them due to the fact that they often cannot adopt any size, depending on their nature; even third party variables, which depend on the values of these decision variables, may also have their own constraints. Below, the constraints involved in the model are listed: • Minimum separation between platforms: they cannot be placed closer than a certain distance (because of the turbulence at the wake of the turbines, the mooring lines length, etc.).

•
Site delimitation: the available space cannot be infinite as the model considers a single wind and waves climate for the whole site. The site area is defined with a radius of 10 km around the center coordinates.
These two first constraints affect directly the coordinates of the platforms, and indirectly the decision variables that determine their arrangement (d. var 1 to d. var 5 ).

•
Depth range: which is determined by the depth boundaries of the mooring system (from 40 to 200 m, for the case study presented later). The depth range may result in restricted areas where the seafloor is either too shallow or too deep.

•
Minimum separation between platforms and substation: similarly to the first constraint, due to security and operational reasons.

•
Rated power of the export cable: so the algorithm chooses an export cable that can bear the farm's maximum output power.
On top of these constraints, the eight decision variables specified have associated boundary parameters that define the maximum and minimum values they can adopt. Technically, this would not be necessary as long as the constraints are fulfilled. However, setting these values is convenient so as to minimize the computing time (e.g., setting the orientation angle from −180 • to +180 • ), so that the optimizer does not try values out of this range.
As the PSO is a single-objective algorithm, all the presented constraints have associated penalty values that prevent the algorithm from choosing a non-feasible result as an optimal solution. These penalties are modeled as a function of the degree of unfulfillment of the constraint, i.e., the further the variable is from the constraint threshold, the higher the objective function (the LCOE) will become, thus the algorithm will be encouraged to progress towards feasible solutions.

Performance Improvement
Several parameters can be modified in order to optimize the performance of the algorithm. This performance is basically determined by two aspects: the computing speed of the algorithm (how long it takes to find an optimized minimum), and the ability to avoid the choice of a local minimum as the final solution. It has been tried to deal with these issues as explained below: Once the algorithm has found what looks like the surrounding area of a minimum, it often starts to modify the values of the particles (decision variables) in a very small scale. Technically, this still improves the output result, but it actually does so slowly and irrelevantly, as the small changes may improve the output value in very small values. Apart from not really making a difference and preventing the algorithm from finishing the optimization quickly, the resolution of the results at that stage is too precise given the model's uncertainties and error tolerance. To prevent this from happening several actions were tried: First, all the decision variables (not only the integers) were rounded to a certain decimal base that was regarded as a reasonable resolution according to its nature. For example, the orientation angles were rounded to integers, having a change resolution of 1 • , or the separation between platforms were rounded to the nearest hundred meters. However, it was found that this solution resulted in a termination of the optimization too prematurely, as the algorithm would often see no change in the objective function when modifying the values of the inputs, sometimes preventing it from finding and jumping to a new better local minimum.
Second, the rounding was also applied to the LCOE output value, which was found to provide better results. This solution satisfactorily stops the algorithm from trying to find irrelevant better solutions for too long, however still with a risk of ending the optimization too soon in a non-global minimum. Something similar is achieved when modifying the option TolFun of the algorithm, which determines the tolerance in the change of the objective function for the algorithm to consider it as an improvement.
The third modification that can be made is adjusting the values of the previously described inertia weight (InertiaRange option in MATLAB), which can also be dynamic (i.e., it can change along the optimization loops), the self-adjustment weight (SelfAdjustmentWeight), and the social adjustment weight (SocialAdjustmentWeight). Optimum values for them and for the amount of particles in each iteration (SwarmSize) were found in literature [24] given the number of input variables, although this also very much depends on the structure of the function being optimized. Three sets of optimum parameters are given as in Table 1. Being the problem dimensions of 10, the case found in literature closest to the eight decision variables of this optimization case. After implanting these changes, an improvement was observed in the results but with an increase in the optimization time, as presented in Section 3.2.2.
The last modification made to improve the performance of the algorithm was running successive optimizations, defining each time the optimum result found in the previous optimization process as one of the initial particles/decision variables (InitialSwarmMatrix) of the next loop. A slight improvement in the final value was found as a result. However, again the risk of limiting the search area of the algorithm to a local minimum exists. The methodology followed to mitigate this type of error is detailed in the next section.

Case Study and Results
The model was applied to minimize the LCOE of a farm in Sligo bay, at the northwest of the Republic of Ireland (Figure 3), an area that presents appropriate conditions for the P80 concept. Below, the application is introduced and subsequently the details regarding the optimization process carried out are presented.

General Description
The Republic of Ireland is one of the countries in Europe with the highest share of wind power production, entailing a 29% of the country's electricity demand by 2018. Furthermore, the State plans to increase the share of renewable generation from 30% to 70% by adding 12 GW of clean energy capacity by 2030, 3.5 GW of which are reserved for offshore sources [25,26]. In this sense, the new Renewable Electricity Support Scheme (RESS) includes Floating Feed-in Premiums, a relevant feature for the P80 [27]. Focusing in the Sligo area, even though it is not one of the regions with the highest electricity demand in the country, the phasing out process of nonrenewable energy generation, the ongoing North-West project (that is supposed to strengthen the grid in the area [28]), and the mentioned optimal wind and waves conditions, make the bay an appropriate area for a P80 farm.
The region presents a number of substations that may serve as a connection of the farm to the electrical grid. Most of them have a rated voltage of 110 kV, except for the 220 kV substation placed in Srananagh, 10 km away from Sligo town. This was the substation chosen for the project as it allows the connection of a big generation as the one modeled in this work. The area also presents two harbors, from which the Sligo harbor was chosen due to its bigger size that makes it more suitable for the installation, O&M, and decommissioning. Regarding the environmental impact, the only protected areas in the bay are the estuaries of the rivers; as for the issues involving the main stakeholders in the area, neither big cities nor major shipping routes can be found, so the visual and operational impacts of the farm are minor issues in this case.
After analyzing the area's grid strength using information from the substation's characteristics, an optimum farm size of 25 platforms was found [17]. A total of 16 possible sites, separated from each other by an average distance of 20 km, was considered throughout the bay, as it can be assumed that the climate is constant in each one of these sites. Six of them (1, 2, 3, 14, 15, and 16) were discarded prior to the optimization process due to factors like the depth constraints, and the proximity to shore (see Figure 4).
In terms of the energy yield calculation, the 8MW wind turbine developed in EU FP7's project, LEANWIND [29], with a 110 m hub height, was chosen. Accordingly, the wind data at the given height was obtained from the database of the Global Wind Atlas 3.0 [30], and that of the waves was obtained from the Danish Hydraulics Institute (DHI) [31]. Non-surprisingly, southwestern wind and waves resources were found to prevail, presenting an increase in energetic resource with distance to shore, as a trade-off with the expected growth of both the capital and operational expenses. Therefore, the LCOE minimization problem first requires the selection of a site for the farm, and then an optimization of the layout (arrangement of the platforms).

Processes Involved in the Optimization Model
In this section, the values chosen for the problem's variables are justified. Afterwards, the steps used along the optimization process are presented.

Definition of Constraints and Parameters
• Two main factors influence the minimum separation between the platforms: the turbulence at the wake of the wind turbines with the subsequent undesired loads in the neighboring turbines, and the horizontal length of the mooring lines which, given their dynamic nature, cannot intersect each other. Taking into account these considerations, the constraint is modeled as follows.
• It was decided to use a minimum separation between the platforms and the offshore substation of 1000 m. Given that the mooring system of the substation is considered to be much more simple and static than that of the P80s, this separation is more than enough. Furthermore, it is assumed that no further interaction between substation and generators exists in terms of wind blocking and waves shadow.

•
Sixteen types of HV export cables were modeled in this case analysis, each of them with a nominal current (I r ), section (Sec), nominal voltage (V r ), rated power (S r ), and metric cost. All these cables are subject to be chosen as long as they fulfill the maximum power output constraint.
The optimizer thus intends to find a balance between the total cost and transmission losses.

•
The site delimitation was set to a circular area with a radius of 10 km around the site's center coordinates, due to the fact that the wind and waves resources cannot be considered constant at higher distances. Furthermore, more extensive areas are not necessary for big farms like this one. • Finally, the discount rate value that has been considered according to the type of technology and which has been used throughout the optimization process in this case analysis, is that of 9% [33], although this is one of the parameters with the highest level of uncertainty.

Optimization Process Methodology
The optimization process to choose the best site and determine the optimized layout was carried out in three consecutive steps, presented with detail below.

•
Step 1: finding the best sites The first part of the process had the goal of discarding the least feasible sites considered up to this point. The model was run for every site with the same identical decision variables and layout (shown in Figure 5), composed of 25 platforms. This way, the long optimization process could be avoided in many of these sites. The site's mean depth and the AEP wind , AEP waves , and LCOE values obtained in the different sites are shown in Table 2. The LCOE values are presented as dimensionless numbers relative to the most optimal result obtained throughout the optimization process, which was given a reference value of 100. As it can be deducted from the results, Sites 4, 7, 8, and 11 (in green) were found to provide the best figures due to a better balance between the energy production and the costs associated to the distance from shore. It is easy to realize that these four sites are the ones placed closer to shore as the operational, installation, and electric costs are lower, and in these sites, the seafloor is more shallow (although this last factor does not have such a big impact in the LCOE). On the other hand, even though both the wind and waves energy production in the rest of the sites is generally higher (with the exception of Site 4), this fact is not enough to balance their longer distance to shore, which increases their costs significantly. This is specially observed in the sites marked in red (Site 6, 10, and 13), the ones placed the furthest from shore. Given the significant difference in the LCOE value, only the sites displayed in green in Table 2 were chosen for the second step of the process.
As discussed more in depth further in the paper, a conclusion that may be extracted from these preliminary results is that in this area the climate's power growth with distance is not steep enough for the technology's costs increase. Also, note the big difference obtained between the wind and waves energy generation. In every site, the latter only represents 10-15% of that of the wind.

•
Step 2: optimization and choice of the best site The second step of the process consisted in carrying out the first optimization in the selected best sites. From the results obtained in the previous step, it seems like Site 7 is the optimum as its LCOE value is approximately 2 points lower than that of the second best site (Site 4). However, given that the energy production of this second one is higher (better wind resource), the optimization of the layout could have a more relevant effect in the minimization of its LCOE.
A new site dismissal was made. Given the nature of the model, which takes into account the linear distances (between site, harbor, onshore substation and shore connection point), a modification should be made in Site 11. As it is located behind a cape, the length of the export cable and the distances of the O&M and installation operations are actually higher than the ones considered in the first approach as their path would need to stretch around the prominent land. These factors, together with the fact that its LCOE is one of the highest in the selected group (and it cannot be minimized too much since the site's wind and waves power is not the best), were considered enough to also discard Site 11 for this step.
After a few prior optimization attempts, it was noticed that the integer decision variables of the optimization problem presented some issues, as the algorithm would get trapped in a localized minimum. Two integer decision variables are part of the optimization: the number of rows and the choice of the HV export cable. Regarding the latter, after running the model several times, it was observed that given the orders of magnitude of the distances to shore and the power generation, one of the cables was providing the best results in every case: the 220 kV, 800 mm 2 , and 775 A cable. Therefore, this decision variable was eliminated and this specific cable was set as a fixed parameter in every case.
As for the number of rows, the solution was not so simple. Five different optimizations for each number of rows were carried out in one of the three sites: Site 4. This way, similarly as it was done when discarding sites in Step 1, the best numbers of rows were found and the optimization process with other less optimal numbers of rows in the rest of the sites was avoided. An example of these types of successive optimizations can be seen in Figure 6, where the graph on the left depicts the improvement in the LCOE in each optimization process until it stabilizes, and the graph on the right shows the layout that provides the newest value of LCOE. Figure 6. Example of simulation for eight rows (at the left, LCOE obtained in each iteration; at the right, the optimal layout found by the algorithm). Therefore, given the consistent and better results obtained in Site 4 for the numbers of rows equal to 5, 6, 7, and 8 (although this last option provided quite unorganized layouts), the same simulations were then run for the rest of the sites. The average results of the five optimizations made in each case are shown in Table 3: Clearly, it is Site number 7 the one that presents the best results with 5 and 7 rows. In all the cases, this first optimization process reduced the LCOE by about 2 points from the values obtained in the previous step. As can be deducted from the table, even though the energy resource is better in Site 4, the optimization of the layout did not mean a significant improvement in the LCOE value with respect to the other two sites, and Site 7 still resulted to be the one with the lowest energy cost. Once the best site has been chosen, the final layout is obtained in the last step.

•
Step 3: final optimization In the last step of the process, the optimizations for Site 7 in the cases of 5 and 7 rows were carried out. In each case, this was done in two ways to make sure that all the results span was covered. Each one of these types of optimization was run 10 times.
In the first type of optimization, the default PSO parameters of inertia, self, and social adjustment weights were used (dynamic ω = 0.1 − 1.1, c 1 = 1.49 and c 2 = 1.49). Furthermore, the five results obtained in the previous step were given to the algorithm as an initial solution (InitialSwarmMatrix). This could have two effects: either it makes the optimization start from a good solution (if it is considered to be close to the global optimal), and therefore it allows it to further improve it, or it prevents the algorithm from finding other existent minimums that may be better.
Due to this second possibility, an alternative way to carry out the optimization was attempted: one of the optimal sets of PSO parameters proposed in literature (shown in Table 1) were used instead: S = 63, ω = 0.6571, c 1 = 1.6319 and c 2 = 0.6239. In this case, no initial sets of decision variables were given to the algorithm, leaving this option to be randomly decided by the PSO as it has been done in Step 2.
After a long process that involved running the optimization ten times for each one of the four cases (each optimization taking an average of 30 to 60 min in the case of the ones with the default variables and an hour to an hour and a half for those using the recommended parameters), the five best values obtained in each case are the following.
It is quite clear how the optimizations run with the default parameters that make use of the initial swarm from Step 2 provided results that slightly improved those from the previous step. However, when they are compared to the optimizations that make use of the recommended parameters, it is evident how the former are stuck in a localized minimum that may be improved.

Results
In the view of the results obtained and presented in the previous section, several final layouts were considered candidates to be chosen as the most convenient. Even though reducing the value of the LCOE of the project is the ultimate goal, other considerations were also taken into account.
Among the layouts that provide the best LCOE output, two clear groups of very similar arrangements were found. Within each of these two groups, the different solutions differ in small variations in the orientation angle or separation between platforms which result in the very small difference of LCOE observed in Table 4. These two main groups of arrangements are depicted in Figures 7 and 8.
After considering the pros and cons of both layouts, it was decided that it is the arrangement depicted in Figure 7, with five rows, the one that would be more convenient in a real case. On the one hand, the small difference of LCOE between both cases (0.07 points) is, most likely, below the error threshold of the model, and therefore it is not regarded as a determining factor. On the other hand, it was considered that the second arrangement is quite unorganized, and therefore more messy for the visual impact. It also makes it more difficult to establish routes for vessels operating in the farm.   The input variables of the selected layout are shown in Table 5. In a first intuitive outlook of the the optimized input variables, it is easy to notice that these are quite logical. First of all, for every layout solution, the position of the offshore substation is always located in the fourth quadrant of the site (at the southeast), as it is oriented towards the direction where the inland substation is located, minimizing the export cable length and cost, as seen in Figure 9. The distance from the farm to the onshore substation (~80 km) and the harbor (in the order of 60 km) can be observed, being that to the closest land~15 km. Again, note that the land closest to the farm is quite uninhabited, only finding very small villages in the surroundings, like Glencolmcille, with only a few hundreds of inhabitants. Figure 9. Location of the farm in the area. In red: the subsea stretch of the export cable. In orange: its inland stretch.
A more close-up look to the site, depicted in Figure 10, reveals more information about the farm. Given the influence of the inter-array cables cost, the offshore substation is quite centered in the site, as the algorithm finds a balance between the costs and transmission losses of these two types of cables (the influence of the different depths in the mooring cost should not be forgotten either, although its cost contribution is quite low). Second, and with a higher influence in the LCOE output, it makes sense that the orientation of the rows and the position of the WTGs tend to minimize the wind wakes when the main directions of the wind are taken into consideration. For this reason, the rows directions have resulted to be perpendicular to the main wind direction observed in the wind rose. Given the presented arrangement of the wind farm, the energy values obtained in the model are those shown in Table 6. It is quite evident how the algorithm tends to minimize the wind wakes losses given its low values, one of the advantages of offshore areas with low spatial restrictions. Furthermore, as mentioned above, the waves shadow range is much smaller than that of the wind wakes and this is shown in the output values. In this case, the result is an energy loss due to the wave's shadows ten times smaller in proportional terms than that caused by the WTG's wakes. The rest of the parameters depicted are not as dependent on the layout, but rather on the site's location and the technology.

Conclusions
The main objective of the study was to develop a tool integrated by various (technical, mathematical, and decision) processes that allows to minimize the average lifetime costs with respect to the energy generation of an offshore farm composed of the Floating Power Plant's hybrid wind and wave concept. The tool consists of waves and wind energy generation and costs models, which was then used to carry out a layout optimization process in order to minimize the LCOE of the hybrid wind and waves P80 concept.
The focus of this paper is the undertaking of optimization process that carry out the cost minimization by modifying certain decision variables. Given the problem's non-convex nature and the presence of discrete parameters, a metaheuristic algorithm was needed. For this preliminary study, the Particle Swarm Optimization (PSO), based on animal swarm behaviors, was chosen. The definition of the decision variables was one of the main determinants of the optimization process, as they depend on the physics and logic after which the model is built, and the nature and number of variables also condition the optimization process' success. The decision variables that were finally implemented are the platforms' positions (arranged in rows for visual and operational considerations), substation position, and the export cable choice. Furthermore, constraints depending on the variables nature, and their correspondent penalties were implemented. The algorithm is able to find optimized solutions but a few issues were found regarding integer input variables, which were addressed in the specific case study.
The case study was carried out off the west coast of Ireland, an area with a convenient wind and waves climate. A farm size of 25 platforms and a grid connection through an onshore 220 kV substation were chosen by means of a study of the area's grid strength. After a qualitative and quantitative process of elimination that involved three phases, Site 7 was chosen, with an arrangement of five rows oriented from northwest to southeast. It was found that given the amount of sites studied and possible distributions, it was convenient to follow a plan to guide the tool along the optimization process. In this regard, it was concluded that there is an interesting path that can be followed to automatize this type of optimization process as it was found that the PSO has some limitations in terms of time consumption and finding consistent solutions. This issue may be solved in the future with other types of meta-heuristic algorithms, most likely multi-objective ones such as the ASO.
After initial sensitivity and uncertainties analyses, the results were put in context, finding that the factor with the highest influence in the final cost, specially considering its high degree of uncertainty, is the discount rate, with a variation effect on the LCOE of approximately 14%, for a 2% change of the parameter (from its base 9% value to a 7%). Furthermore, situations with certain modifications were studied. A plausible short-term scenario with certain parameter variations was calculated, obtaining a reduction of 74%, an LCOE range significantly close to that of bottom-fixed wind farms currently being auctioned (60 to 75 A C/MWh). It is considered that the technology presents good potential as there is room for costs reduction, specially in the structure's material (first cost component of floating concepts).  Acknowledgments: Thanks to Sarah Thomas, Pilar Heras, Anders Køhler, Morten Kramer, and Chris McConville for offering help and advice and providing the information needed from FPP. I would also like to thank Niels Mortensen, Daniel Hermosilla, Andreas Bechmann, and Nikolay Dimitrov, from DTU Wind Energy, for offering their knowledge and support. Special gratitude to Alfredo Peña, who granted me access to his work and let me use it in the tool, Neil Davis, who provided crucial wind data for the model, and Hugo de Sevin and Gyde Ohlsen, for allowing me to have a basis from where to build this work and take inspiration.

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

Abbreviations
The following abbreviations are used in this manuscript.