The optimization of spare-part logistics and maintenance scheduling constitutes a fundamental challenge in the management of floating offshore wind (FOW) farms. This research introduces a model designed to explore the complexities of this industrial problem, providing a robust framework for decision-making. Integral to the model’s conception is the acknowledgment of the holistic nature of spare part distribution, which is influenced by fluctuating economic conditions, such as price indices and industry demand. These factors, coupled with the critical need for timely part delivery to prevent turbine downtime, add layers of complexity to the logistics process.
2.1. Problem Statement
The model implemented in this research aims to optimize spare-part logistics and maintenance decisions within a network of geographically distributed offshore wind turbines.
The logistics of spare parts are intrinsically complex due to several factors, prominently including the stochastic nature of part failures and the subsequent need for replacements. As depicted in
Figure 1, the model developed under this research centralizes its analysis on a maintenance center/warehouse (fed by different geographically distributed part suppliers) and the respective connections to multiple offshore turbines. These connections are not merely physical but also temporal, with routes from parts suppliers to maintenance centers and turbines being quantified in hours.
In
Figure 2, two significant dependencies under the proposed framework are highlighted: the unpredictable nature of part supply, which is contingent on factors like price indices and industry demand, and maintenance operations, which are greatly influenced by weather conditions. These dependencies underscore the complexity of the model, which must account for these uncertainties through its optimization equations and the respective variables assumed.
As stated, maintenance operations within FOW farms are heavily contingent upon variable and often adverse weather conditions, which can significantly disrupt planned activities. The proposed model incorporates a dynamic approach to maintenance scheduling that is sensitive to these environmental uncertainties. By integrating the influence of meteorological data, the model is capable of adapting maintenance strategies to improve the timing and execution of operations.
The architecture of the proposed model is modular, encompassing independent yet interrelated components that handle various operational aspects of the maintenance process. These modules, responsible for predicting failures, predicting electricity prices, and forecasting metocean conditions, collectively inform a layered decision-making framework. This framework, schematically presented in
Figure 3, distinguishes between tactical short-term actions, which respond to immediate operational needs, and strategic long-term planning, which shapes the broader scope of maintenance strategies.
The hypotheses of the model are grounded in the challenges faced by floating offshore wind farms. These turbines operate under harsh environmental conditions, and the supply chain for spare parts stretches across vast distances, complicating the logistics and elevating the risk of extended downtimes:
Temporal flexibility in maintenance scheduling: The model assumes that maintenance activities can be scheduled flexibly over short-term and long-term periods. The use of nested GAs allows the interaction between these two temporal dimensions. By analyzing both short-term and long-term spectra and considering their mutual influences, the model hypothesizes that decisions made in one timeframe have significant implications for the other.
Predictive maintenance capability: The model operates under the hypothesis that maintenance needs can be anticipated through the analysis of operational data. This is incorporated through component condition thresholds intervals between inspections as variables, seeking to preempt maintenance issues before they lead to component failure.
Impact of Maintenance on Energy Production and Revenue: The model hypothesizes that maintenance operations have a direct impact on the energy production efficiency of FOW farms; furthermore, it takes into consideration the electricity price forecast to propose the execution of maintenance procedures during periods of lower electricity market prices.
Operational and Environmental Constraints: The inclusion of variables such as weather conditions and logistic considerations (e.g., availability of vessels for tow-to-port operations and the existence of spare parts at the maintenance center) implies that the implemented modeling of the operational efficiency of FOW farms is heavily influenced by environmental conditions and logistical factors.
As any model that simulates a certain reality, the proposed model has some limitations that need to be acknowledged:
Data Dependence: The model relies heavily on the quality and availability of input data, such as accurate failure detection, electricity price forecasts, and metocean forecasts. If these data are incomplete, inaccurate, or outdated, the model’s outputs will be compromised.
Environmental and Regulatory Changes: Environmental policies and regulations impacting O&M activities in FOW farms are subject to change, which may not be fully accounted for in a model that is based on current standards and practices.
Cost Model Simplifications: The model simplifies the complex cost structures associated with maintenance activities of OSW farms, which could potentially overlook factors such as varying labor costs, penalties for downtime, and the financial impact of unexpected delays.
Interaction Effects: The model considers the interdependencies between short-term and long-term maintenance decisions. However, capturing the full spectrum of interaction effects accurately is challenging and may not reflect all the complexities of real-world operations.
Complexity in Real-World Application: While the model may perform well under simulated or controlled conditions, translating its recommendations into practical, real-world FOW farm management involves logistical complexities and may encounter unforeseen implementation challenges.
Acknowledging its limitations, the model is intentionally designed with modularity to ensure flexibility and ease of updates. To put the developed model up to test, the present research will use as case study the Windfloat Atlantic floating OSW farm—as partially revealed in
Figure 1. This farm, implemented off the north coast of Portugal back in 2019, was one of the first grid-connected farms of its kind. It is a pre-commercial 3-turbine farm—from the manufacturer Vestas, with a total output of 25 MW—and is expected to operate for the next 25 years. Its pioneering status as a pre-commercial venture offers valuable insights into the unique challenges and opportunities of floating OSW farms, making it an exemplary testbed for the model’s strategic maintenance optimization features.
2.2. Decision-Making Support Module
A nested GA framework is developed under this research to address the complex and interdependent decision-making processes in OSW farm maintenance. The traditional GA-based approach is expanded into a nested structure to capture the interplay between short-term operational decisions and long-term strategic planning. The developed model is built within the Python environment, leveraging the DEAP (Distributed Evolutionary Algorithms in Python) library [
17]. The DEAP library provides a versatile toolbox for evolutionary computation, enabling the implementation of complex genetic operations such as mutation, crossover, and gene correction within a constraint optimization context.
For the outer layer, known as the “Captain GA”, the ‘eaMuPlusLambda’ evolutionary strategy is employed, which uses a combination of the best individuals (µ) and the offspring (λ) to create a new generation, thereby maintaining elite solutions while encouraging diverse exploratory steps through crossover and mutation [
17]. This approach promotes robust search capability within the solution space and is particularly effective when contending with the stochastic nature of the failure models in OSW farms.
Within the proposed framework, schematically shown in
Figure 4, the inner GA, also referred to as “Stewart GA”, is dedicated to optimizing immediate maintenance activities that are sensitive to operational exigencies and rapidly changing environmental and meteorological conditions. The resulting decisions are then integrated as parameters into the outer GA, referred to as the “Captain GA”, which is tasked with improving strategic maintenance parameters that secure long-term operational efficiency and farm viability.
For both the “Captain GA” (outer layer) and the “Stewart GA” (inner layer), the ‘eaMuPlusLambda’ algorithm was utilized. The GAs implemented within this framework are also characterized by the additional features:
Mutation: This operation integrates randomness into the gene pool, creating opportunities for novel solutions to arise.
Crossover: By combining traits from paired individuals, this process encourages the propagation of advantageous genetic information.
Gene correction: This mechanism ensures that solutions adhere to the specific operational constraints of OSW maintenance, thus maintaining the relevance and applicability of the developed model.
Elitism: Within each model run, an elitism strategy is deployed to preserve top solutions, ensuring the retention of high-quality genetic material for subsequent generations.
Hall of Fame: To accommodate the stochastic nature of the failure module (presented next), a hall of fame across multiple GA runs is maintained. This repository of best solutions stores the best solutions for varying parameter sets, providing a robust basis for the model’s conclusions.
Average of results: By averaging outcomes from multiple simulations, the randomness of failure occurrences is mitigated, leading to a more consistent and reliable assessment of short- and long-term maintenance strategies.
The primary objective of the “Captain GA” is to enhance the long-term profitability of the FOW farm by optimizing maintenance parameters. It achieves this by adjusting the interval between inspections and setting precise repair thresholds for each subsystem component. Its parameters include:
Interval between Inspections: This variable, ranging from 1 to 12 months, determines the frequency of maintenance checks, which is crucial for early fault detection and reduction of unscheduled downtimes.
Inspection Reliability Thresholds: These thresholds, ranging from 0 to 100%, represent the condition at which a certain turbine part is deemed to require repair or replacement based on visual inspections. These thresholds are independently considered for each part analyzed and are detailed below.
The fitness function for “Captain GA” aims to maximize the revenue over a long-term horizon. It is formulated to account for operational earnings while deducting the costs associated with maintenance activities. Mathematically, it can be represented by a function that calculates the gross profit throughout the entire period of analysis, considering the income from energy production minus the costs from inspections, repairs, and associated logistics throughout that long-term period. The function can be expressed as follows:
The “Captain GA” operates under the assumption that shorter intervals between inspections (IbIs) increase the likelihood of fault detection before total loss of reliability of a certain mechanical part of the turbine, albeit at a higher operational cost. The repair thresholds are set based on the damage estimated by technicians from in situ inspections, a proxy for each component’s health. Thus, finding a balance between the frequency of inspections and the repair thresholds is critical, as overly frequent inspections can lead to unnecessary expenditure, while infrequent inspections might result in missed detections and subsequent catastrophic failures
Conversely, the “Stewart GA” focuses on optimizing maintenance procedures within a short-term window, primarily the upcoming thirty days from the date at which the failure is detected. Its variables are:
Best Hourly Slot for Maintenance: This variable selects the optimal hour within a 30-day span for conducting maintenance, balancing operational availability with the least disruption to energy production, and, particularly, total revenue during that period.
Component Availability: A binary variable indicating the availability of the necessary spare part in the warehouse, directly influencing the decision to undertake immediate maintenance action or to wait for the part shipment and arrival.
The fitness function for the “Stewart GA” is designed to maximize the profit from each maintenance procedure conducted within the short-term period. This function considers the immediate operational revenue, adjusted for the potential costs of maintenance and the impact of component availability, as presented next:
where E
prod is the hourly electricity production vector for the 30-day timespan analyzed, E
price is the hourly predicted electricity price for this 30-day timespan, C
vessel is the renting cost of the vessel required for each maintenance procedure, C
HR is the cost of technical workers required for each maintenance procedure, C
maint is the cost of executing each maintenance procedure, and C
ware is the cost associated with storing one component to then be used as replacement—in this respect, if a spare part is not available in the warehouse, a waiting time for ordering and shipment of that specific part is considered, which reduces the flexibility of the farm owner to decide which hourly-slot is the best for repair. C
inspection is the cost of inspecting the turbine, according to the IbI established on the “Captain GA”.
These two optimization equations at the core of the proposed model are critical, as they weigh the costs of maintenance and inventory against the penalties of asset downtime. By capturing these factors, the model searches for a strategy that optimizes financial and operational efficiency.
In practice, the model operates by evaluating the operational condition of each turbine over time, factoring in the interval between inspections (IbI) and the set inspection reliability thresholds from the “Captain GA”, alongside the failure predictions from the failure module. It systematically proposes and analyzes maintenance actions in response to each identified failure event. Turbines are presumed to be equipped with Structural Health Monitoring (SHM) systems that detect actual failures as they occur, an assumption that aligns with studies indicating SHM systems can allow for less frequent inspections without sacrificing revenue over the project’s lifespan [
18]. However, SHM’s detection capabilities are considered, under the present research, to only detect already existing failures. In contrast, regular inspections aim to preempt potential failures by measuring a component’s reliability against established thresholds. If reliability is found to be below these thresholds during an inspection, a maintenance action is scheduled, which is aimed at fully restoring the component’s reliability. A depiction of an example for the timeline of a turbine under the developed model is proposed in
Figure 5.
2.3. Failure Detection Module
The developed model integrates within a module designed to forecast component failures. Given the absence of real-world failure data for offshore wind farms, the module relies on failure rates and operational cost data derived from the research by James Carroll et al. [
19], distributing failures stochastically over the entire duration of the analysis period considered for model simulation. Carroll et al. [
19] attribute the severity of the failures of each subsystem based on the cost range of correcting a certain failure. The failure distribution output of the farm tested under the present research is obtained using a Monte Carlo simulation—a computational algorithm that relies on repeated random sampling to obtain numerical results, which in this case, predicts the timing of failures. The module operates under the assumption that failures follow an exponential distribution, which is a common hypothesis for failure analysis in the absence of more detailed data. The failure rates per component for minor repairs, major repairs, and replacements are presented below in
Table 1.
The module uses these rates to generate a structured dataset that describes the failure events’ timings. A partial example of a list of simulated failures obtained from the proposed module is provided below in
Table 2.
Table 2 presents an hourly-based schedule of the main farm events, including type and severity of failures per turbine, as well as the inspection dates.
Given the stochastic nature of the exponential distribution, the reliability of the predictions from a single run of the model would be limited due to the inherent randomness of the process. To address this variability and enhance the robustness of the model’s forecasts, multiple instances of the model are run. Each iteration of the model uses a distinct sequence of failures based on the same underlying probability distributions, but with different random seeds, leading to varied outcomes. By aggregating the results from numerous simulations, the model averages out the randomness associated with individual runs. Consequently, the multiple runs collectively form a more comprehensive and accurate picture of the expected maintenance requirements over the analysis period, enabling more informed and strategic decision-making for maintenance scheduling.
2.4. Electricity Prices and Metocean Forecast Modules
Within the scope of this research, two distinct but interconnected modules have been integrated to process critical environmental data and energy market dynamics. The model assimilates historical electricity price data from 2017 to 2022 obtained from the OMIE database, which is the entity in charge of the Iberian energy market (MIBEL) [
20]. These historical data are treated within the model as a predictive estimate, forecasting prices 30 days into the future to simulate a realistic decision-making environment. Such an approach ensures that maintenance scheduling will align with periods of maximum energy price to optimize revenue.
For metocean conditions, the model integrates wind and wave data acquired from the ERA5 reanalysis dataset, corresponding to the geographic coordinates nearest to the offshore wind farm. ERA5 provides access to a dataset of weather station data collected since 1950 to deliver hourly estimates of atmospheric, oceanic, and land climate variables over a global grid, with each cell spanning across 30 × 30 km [
21]. For the purpose of this study, the following key parameters are used:
Significant wave height, providing an average measure of the highest one-third of waves.
The u component (eastward) and v component (northward) of the 10 m wind. These values have been adjusted to turbine height by extrapolating the wind speed from the standard 10 m height to 119 m, using the logarithmic wind profile equation, to reflect the actual hub height of the wind turbine considered.
The absence of public data on the Vestas 8.3 MW turbine was overcome using the DTU 10 MW reference turbine [
22,
23]. This reference turbine serves as a scientific standard, allowing for the comparison of research outcomes between different research groups. By using the same turbine model, results can be consistently compared across different research. Employing the reference turbine’s power curve (see
Figure 6) in conjunction with real-time electricity pricing and metocean data enables an accurate hourly forecast of energy production and associated revenue. This becomes a feedback mechanism for the model’s decision-making, informing the selection process of optimal maintenance timing to enhance the operational profitability of each turbine.
The model also accounts for production halts during major repairs and anticipates downtime when spare parts are in transit. This allows for the preemptive shutdown of turbines to prevent compounding damage, a measure that safeguards other components and overall turbine integrity.
2.5. Cost and Logistic Assumptions
This section discusses and presents the cost and logistical assumptions assumed in this research for maintenance operations that require towing turbines to port (tow-to-port). The data draw on insights from Carroll et al. [
19] work and on the internal databases of WavEC, reflecting the direct and indirect expenses associated with different maintenance activities—see
Table 3.
The maintenance costs for the key subsystems requiring tow-to-port procedures are classified into minor repairs, major repairs, and major replacements. Each category assumes distinct costs, reflecting the severity and complexity of the maintenance required. The expenditure associated with minor repairs is comparatively low, while major repairs and replacements incur significantly higher costs due to their comprehensive nature and the need for specialized equipment and labor. For the present research, it is assumed that TTP replacement events can be produced in just 72 h. Nonetheless, this value can be highly optimistic, particularly for newer farms or farms that do not have support ports in the vicinities. For example, recently, at the Kincardine OSW farm, which is floating and also counts on the Windfloat semi-submersible platform, total time for one TTP operation was of an astounding 94 days [
25]. This operational experience highlights the need for optimizing O&M actions for floating OSW farms, as well as the need for adequate support infrastructure to accommodate the exigent needs of newer O&M strategies.
The logistic component of the model considers the type of vessel required for each maintenance activity—CTV (Crew Transfer Vessel), SOV (Service Operation Vessel), and DOCK (Docking Services). The choice of vessel is dictated by the specific maintenance operation, with DOCK being essential for major replacements (assuming the use of an SOV to execute the towing process). The associated costs are influenced by factors such as rental rates, human resources, and the duration of the maintenance procedure. The costs of renting vessels and the potential downtime must be factored into the overall budgeting and scheduling of maintenance operations. The research assumes at this stage that vessels are always available, although this may not reflect the real-world scenario and could introduce a level of uncertainty in the logistics planning.
Table 4 presents the vessels required for each failure type and for each turbine part considered under this research.
The model accounts for hourly rates and daily costs for each vessel type. These rates include not only the base rental cost but also fuel charges and other operational expenses. For instance, the CTV has a lower hourly rate compared to the SOV and DOCK, but its operational costs per day might be less favorable when considering longer maintenance windows and conjoined maintenance operations. These values are presented in
Table 5.
The model also incorporates maximum operational wind speeds and wave heights for each vessel type, which are essential for ensuring safety and feasibility of maintenance operations. These thresholds are aligned with the industry practices and were retrieved from WavEC’s database. If metocean conditions exceed these thresholds, maintenance activities will be delayed, which will consequently impact the overall scheduling and costs—see
Table 6.
Lastly, the availability of spare parts in the warehouse and the lead time for ordering components from manufacturers are crucial parameters for the model (in fact, the availability of each part in the warehouse is one of the genes assumed for the “Stewart GA”). These factors influence maintenance scheduling, as the lack of immediate part availability could lead to extended downtimes, especially in the event of unscheduled maintenance requirements. The integration of these cost and logistic assumptions provides realistic estimates of the direct costs associated with specific maintenance activities.
2.6. Genetic Algorithm Parameters/Simulation Parameters
In this research, two distinct sets of parameters, which frame the evolutionary search process within each algorithm’s loop, were defined for the “Captain GA” and the “Stewart GA”. These are summarized below in
Table 7.
The parameter “Number of Averages” is unique to the “Captain GA”, reflecting the emphasis of the developed model on long-term strategic decision-making by averaging results over multiple simulations to address the stochastic nature of the failure generation model.
The chromosomes of both genetic loops present the blueprints of the algorithms, being formed by the decision variables essential to maintenance scheduling and decision-making. The “Captain GA” chromosome is constituted by parameters that influence long-term strategic decisions, including the interval between inspections (IbI) and repair thresholds for various components, where each gene corresponds to a vital decision variable for long-term planning. Specifically, this chromosome can be represented as:
where the IbI ranges from 1 to 12 months, and the repair thresholds for each of the
= 5 turbine parts considered may vary from 10 to 90%. The “Stewart GA” chromosome encodes operational decisions such as the optimal maintenance hour within a short-term horizon and the availability of spare parts in the warehouse. It can be represented as:
where the optimal maintenance hourly slot can take a value between zero and 720 (7 days), and the spare part availability is a binary variable, which may take a value of zero (no spare part available) or one. These structures allow the algorithms to efficiently seek optimal maintenance parameters that maximize the short-term and long-term profitability of the farm.