Optimization of wind farm portfolios for minimizing overall power fluctuations at selective frequencies – a case study of the Faroe Islands
 ^{1}Faculty of Science and Technology, University of the Faroe Islands, Tórshavn, Faroe Islands
 ^{2}Department of Wind Energy, Technical University of Denmark, Roskilde, Denmark
 ^{1}Faculty of Science and Technology, University of the Faroe Islands, Tórshavn, Faroe Islands
 ^{2}Department of Wind Energy, Technical University of Denmark, Roskilde, Denmark
Correspondence: Turið Poulsen (turidp@setur.fo)
Hide author detailsCorrespondence: Turið Poulsen (turidp@setur.fo)
Hourly modeled wind turbine power output time series – modeled based on outputs from the mesoscale numerical weather prediction system Weather Research and Forecasting Model (WRF) – are used to examine the spatial smoothing of various wind farm portfolios located on a complex isolated island group with a surface area of 1400 km^{2}. Power spectral densities (PSDs), hourly stepchange functions, and duration curves are generated, and the 5th and 95th percentiles and the standard deviations of the hourly stepchange functions are calculated. The spatial smoothing is identified from smaller highfrequency PSD amplitudes, lower hourly fluctuations, and more flat duration curves per installed wind power capacity, compared with single wind turbine outputs. A discussion on the limitation of the spatial smoothing for the region is included, where a smoothing effect is observed for periods of up to 1–2 d, although it is most evident at higher frequencies. By maximizing the smoothing effect, optimal wind farm portfolios are presented with the intention of minimizing overall wind power fluctuations. The focus is mainly on the smoothing effect on the 1–3 h timescale, during which the coherency between wind farm power outputs is expected to be dependent on how the regional weather travels between local sites, thereby making optimizations of wind farm portfolios relevant – in contrast to a focus on either lower or higher frequencies on the scale of days or minutes, respectively, during which wind farm power output time series are expected to be either close to fully coherent due to the same weather conditions covering a small region or not coherent as the turbulences in separate wind farm locations are expected to be uncorrelated. Results show that an optimization of the wind farm capacities at 14 predefined wind farm site locations has a minimal improvement on the hourly fluctuations compared with a portfolio with equally weighted wind farm capacities. However, choosing optimized combinations of individual wind farm site locations decreases the 1–3 h fluctuations considerably. For example, selecting a portfolio with four wind farms (out of the 14 predefined wind farm site locations) results in 15 % lower 5th and 95th percentiles of the hourly stepchange function when choosing optimal wind farm combinations compared with choosing the worst wind farm combinations. For an optimized wind farm portfolio of seven wind farms, this number is 13 %. Optimized wind farm portfolios consist of distant wind farms, while the worst portfolios consist of clustered wind farms.
Most nations strive to make their electricity generation less dependent on fossil fuels. Wind power is an attractive solution, being environmentally friendly, mature, and affordable. However, an increase in the wind power contribution to the power grid can cause challenges as wind varies across various timescales. Measures are required to mitigate natural wind power fluctuations in order to balance and keep the grid stabilized. One possible measure is spatial distribution of wind farm siting, see, e.g., Beyer et al. (1990) or Katzenstein et al. (2010). Isolated grids face special challenges, as they usually span across small areas. Therefore, the available wind resource is partly or fully covarying, which poses a serious challenge for maintaining the energy balance of the grid.
The wind power fluctuations in a confined isolated region are examined in this study, in which the Faroe Islands are used as a case study. The Faroe Islands are a small isolated mountainous island group (1400 km^{2}) in the north east of the Atlantic Ocean, ∼ 300 km away from any mainland; see Fig. 1. The topography is complex, and the climate is windy. Battery systems have been installed in the Faroe Islands and can handle wind variability up to timescales of several minutes. Manual setpoint control of hydropower and fossil power plants is presently used to compensate for wind variability on timescales from several minutes to several days. In the future, the contribution from wind power will increase rapidly, and the present manualcontrol regime will be put under pressure. There is therefore considerable local interest in limiting the inherent wind variability on hourly timescales in the Faroe Islands.
The aim of the study is to analyze and optimize the smoothing effect caused by the nature of dispersed wind farm sites in the region. Wind farm portfolios are optimized with the objective of minimizing the 2–3 h power output fluctuations.
There are numerous studies in the literature on analyzing the smoothing effect of spatially distributed wind farm sites on various time and spatial scales, applying different techniques, e.g., correlation analysis (Beyer et al., 1993; Katzenstein et al., 2010; Pearre and Swan, 2018; Giebel, 2000), stepchange analysis (Beyer et al., 1990, 1993; Katzenstein et al., 2010; Pearre and Swan, 2018; Giebel, 2000), visualization of duration curves (Giebel, 2000; Barasa and Aganda, 2016; Katzenstein et al., 2010), and analysis in the frequency domain (Beyer et al., 1990, 1993; Nanahara et al., 2004a, b; Katzenstein et al., 2010; Poulsen and Beyer, 2020). One advantage of frequency analysis is the ability to characterize the extent of the smoothing effect for many frequencies. The reduction in the total variance of aggregated wind farm data can be small even though a significant smoothing effect occurs at higher frequencies, as observed by Beyer et al. (1990). In contrast, time domain analysis – such as stepchange analysis and duration curves – may be more intuitive to interpret. Other examples of studies analyzing the smoothing of spatially distributed wind farms in the time domain include Palutikof et al. (1990), Wan et al. (2003), and Giebel (2001).
The study by Frank et al. (2021) is an example of a recent study examining the smoothing effect of the natural variability in both wind and solar power. They analyzed daily data from 11 European countries. Although their focus was both on wind and solar power, the smoothing effect of spatially distributed wind power can be seen in their results, based on the low correlation coefficients and low joint probabilities of the minimum and maximum wind power between countries. In addition, they optimize the installation ratio between solar and wind power in order to reduce extreme fluctuations. Even though their results cannot be directly compared with the results in this study, their study demonstrates the potential of balancing the natural fluctuations from renewable energy resources by considering the configuration of the resources.
Reichenberg et al. (2014) and Cassola et al. (2008) are examples of studies focusing on minimizing fluctuations in the total wind power production through wind farm optimization. Reichenberg et al. (2014) developed a method for optimizing the locations of wind farms by minimizing the coefficient of variation of the total wind power output time series. Their method is based on sequential optimization of site localization, and it is applied in the Nordic countries and Germany. Cassola et al. (2008) proposed a procedure for optimally distributing the relative size of considered wind farms either by minimizing the variation in the total wind power output time series or by minimizing the ratio between the variation and the total wind energy production. Their method was tested with regards to wind power data at 10 locations in Corsica (France), a complex island in the Mediterranean Basin stretching 175 km in the latitude direction and 80 km in the longitude direction with a mountain chain crossing the islands in the north–south direction and peaks higher than 2000 m. However, the results from these studies cannot be directly compared to the results obtained in this study, as the methodologies are different as well as the geography.
The results in the present study are based on hourly modeled wind turbine power output data, which are analyzed by generating power spectral densities (PSDs), hourly stepchange functions, and duration curves for individual and lumped power output time series. Moreover, the 5th and 95th percentiles and the standard deviations of the hourly stepchange functions are extracted for quantitatively comparing the results. The PSDs are further used to optimize wind farm portfolios by identifying the wind farm combinations with the lowest 2–3 h wind power fluctuations.
The methodology applied in the presented wind farm portfolio optimization study is unconventional – using spectral results as the objective function for minimizing the overall wind power fluctuations, a powerful methodology in the sense that the optimization focuses solely on the preferred frequencies. No such optimization study has ever been conducted for the Faroe Islands region.
The datasets used in this study are introduced in Sect. 2, the methods applied are described in Sect. 3, the results are presented in Sect. 4, and a summary and conclusions are given in Sect. 5.
1.1 A note on ignoring the wind farm smoothing effect – paper assumption
In this study, turbine outputs are analyzed, thus neglecting the potential smoothing effects within wind farms. Various studies investigate the coherence function between time series on a wind farm scale (Schlez and Infield, 1998; ViguerasRodríguez et al., 2012; Vincent et al., 2013); the coherence function can be used as a measure of the smoothing effect between time series as a function of frequency. If the coherence is high, the smoothing effect is small in the considered frequency range. If the coherence is small, a smoothing effect occurs.
The distance between the two farthest turbines in the largest wind farm in the Faroe Islands at the time of the preparation of this paper – with thirteen 0.9 MW turbines – is about 670 m. For a wind speed of 10 m s^{−1} (approximate average values as observed for the Faroe Islands region by Poulsen et al., 2021), a distance of 670 m, and a frequency of (2 h)^{−1}, the analytical coherence model presented by ViguerasRodríguez et al. (2012) yields squared coherence values of 0.81 and 0.92 in the lateral and longitudinal directions, respectively, and the wind farm smoothing effect is therefore expected to be small for frequencies of up to (2 h)^{−1}. For larger wind farms, a wind farm smoothing effect is expected. Although this is outside the scope of this study, a more detailed study could apply methods such as those presented by Nørgaard and Holttinen (2004) or Sørensen et al. (2008) when scaling up turbine outputs to large wind farms.
Hourly Weather Research and Forecasting Model (WRF)generated wind speed data are the primary dataset used throughout this study. Time series are extracted at locations with favorable wind farm conditions. In addition, available empirical data are used for comparison.
2.1 WRFgenerated wind speed data
WRFgenerated data were simulated by Kjeller Vindteknikk and made available to us by SEV (the local power company of the Faroe Islands). The model setup is documented in Haslerud (2019). Multiple output parameters were generated, including the wind speed data used in this study, validated by Poulsen et al. (2021). WRF version v3.8.1 was used for the simulations, with 51 vertical levels (8 in the lower 200 m). Three nested domains were used, with horizontal resolutions of 500×500, 1500×1500, and 7500×7500 m, for the innermost, middle, and outer domains, respectively. The innermost domain covered the entire Faroe Islands. Data were stored every hour during the period from 1 July 2016 to 30 June 2018. The Thompson scheme, the Mellor–Yamada–Janjić scheme, and the NOAH scheme were applied for microphysics, boundary layer mixing, and the surface, respectively. Every 3 h, ERA5 reanalysis data with a horizontal resolution of approximately 0.25 ^{∘} were used as boundary condition (available from the European Centre for MediumRange Weather Forecasts, ECMWF).
2.2 Empirical data
During the period from 1 July 2016 to 30 June 2018, hourly power production data from three wind farms operating on the Faroe Islands are available. The hub heights of all turbines (except one smaller) are 45 m a.g.l. In addition, 10 min averaged wind speed measurements from an additional site at a height of 52.8 m a.g.l. during the period from 21 July 2016 to 30 June 2018 are available. The locations of the four sites are pinpointed in Fig. 2.
2.3 Favorable wind farm site locations
Magnussen (2017) presented a map with current and potential wind farm site areas on the Faroe Islands; see Fig. 2. The areas were selected based on modeled wind resource, distance from the road, distance from the highvoltage grid, terrain slope, populated areas, and presence of water and rivers. These considerations were weighted and merged. Based on additional considerations, as expected turbulent areas, some locations were rejected, and a selection of favorable wind farm sites was composed.
Data at these favorable locations are used throughout the study.
Modeled wind turbine power output data at a height of 45 m a.g.l. from the period of July 2016 to June 2018 are examined. The data processing of the time series is explained in the following subsection. The applied analytical methods are defined in the second and third subsections.
3.1 Data processing
WRFgenerated wind speed time series for selected sites are vertically interpolated to a height of 45 m a.g.l. This height is chosen because all operating wind turbines in the Faroe Islands at the time of the preparation of this paper had a hub height of 45 m a.g.l. The interpolated time series are subsequently corrected for the aliasing effect, as elaborated on in the Appendix.
The wind speed time series are modeled to power output time series using the power curve of an Enercon E44 wind turbine with a storm control function (Enercon, 2012) and a rated power of 0.9 MW. This turbine model is chosen because most of the currently operating wind turbines in the Faroe Islands are of the type Enercon E44.
To be able to compare results, all power output time series are normalized by their rated wind power capacity. When time series are aggregated, the lumped time series is computed prior to normalization.
3.2 Approaches for the characterization of the power output time series
In this study, three analyzing techniques are applied to the time series – spectral analysis, stepchange analysis, and the generation of duration curves – each of which is elaborated on in the following subsections. The first two characterize the time series fluctuations, while the latter characterizes the distribution of the power generation.
3.2.1 Spectral analysis
PSDs are generated for normalized power output time series using discrete Fourier transform. Only the amplitudes corresponding to the positive frequencies of the discrete Fourier transform are extracted. Thus, the integral of the raw spectra with respect to frequency yields about half of the variance of the time series. Due to the inherent uncertainty of raw spectral estimates, smoothed PSDs are generated by dividing the time series into chunks, calculating the spectral results for each chunk, and then averaging over all chunks. For each chunk, the average value is subtracted, and a hamming window is applied. A 50 % overlap between chunks is used.
The length of the chunks is a compromise between the accuracy of the PSD estimates (smaller chunks, i.e., more chunks) and the frequency resolution and the lowest resolvable frequency (longer chunks). In this study, a length of 256 data points was chosen (10 d and 16 h), giving 135 overlapping chunks for the 2yearlong hourly time series. The PSD estimates will therefore be generated for frequencies between (256 h)^{−1} (thus, including PSD estimates for the 3–4 d period of the timescale of migratory lowpressure systems at mid and high latitudes) and the Nyquist frequency of (2 h)^{−1} with a resolution of (256 h)^{−1}.
3.2.2 Stepchange analysis
The stepchange function of the wind power time series is calculated as given in Eq. (1).
where P_{t} is the wind power production at time t, P_{t+1} is the wind power production of the consecutive time step, and ΔP_{t} is the corresponding change in the wind power production between these two time steps. The distribution of the wind power stepchange function is presented with its probability density function (pdf). In addition, the 5th and 95th percentiles and the standard deviation of the stepchange function are identified. Although the hourly stepchange function is not normally distributed, the standard deviation is a measure for characterizing the hourly fluctuations of the time series. The higher the standard deviation, the more frequently large hourly fluctuations occur, and vice versa. The 5th and 95th percentiles are measures characterizing the more extreme hourly fluctuations, whereby the hourly fluctuations are beyond either of these values 10 % of the time.
3.2.3 Extraction of power duration curves
The power output time series data are sorted according to descending magnitude to extract the fraction of exceedance of the power levels, presenting the power duration curve.
3.3 Optimization method
By optimizing wind farm portfolios with the aim of minimizing the fluctuations of the total time series, the most stable wind farm configurations are derived. For this purpose, modeled wind turbine power output time series are used, positioned at the favorable wind farm site locations shown in Fig. 2, and modeled as described in Sect. 3.1. The overall power output time series of the wind farm portfolio, P_{lumped}, is generated by aggregating the turbine power output time series of the combination of wind farms in the given portfolio, as demonstrated in Eq. (2).
where P_{i}(t) is the modeled wind turbine power output time series at location i, and n is the number of considered wind farm site locations in the given portfolio. P_{inst,i} is the relative wind farm capacity at location i with respect to the total capacity of the wind farm portfolio. The relative capacity is used in order to normalize the aggregated time series, so it represents wind power output data per installed wind power capacity, as mentioned in Sect. 3.1. Thus, by definition, the sum of P_{inst,i} must equal one:
where P_{inst,i} is always greater than or equal to zero:
As PSD estimates represent time series fluctuations with respect to frequency, the objective function in the optimization is set to minimize the integral of the PSD estimates of the aggregated power output time series from the wind farm portfolio (Eq. 2) for frequencies between (2 h)^{−1} and (3 h)^{−1}:
with the constraints given in Eq. (3) and (4). Thus, the optimization minimizes the 2–3 h fluctuations of the wind farm portfolio power output time series. These are the highest resolvable spectral frequencies when working with hourly data. Lower frequencies could also be considered but are not chosen because larger periods, on the scale of days, are expected to be close to fully correlated across the Faroe Islands, as the same weather systems are expected to cover the small island region (1400 km^{2}). This is discussed further in Sect. 4.3 and 4.5. The optimized parameters are set to be one of the following:
 1
the normalized wind power capacities, P_{inst,i}, where the locations of the wind farms in the portfolio, i, are fixed;
 2
the locations of the wind farms, i, where the total number of wind farms in the portfolio, n, is fixed, and the normalized wind power capacities are equal at all site locations ($\mathrm{1}/n$).
Optimized portfolios using the first mentioned parameters are presented in Sect. 4.2, while optimized portfolios using the latter mentioned parameters are presented in Sect. 4.4.
All results represent normalized power output time series with respect to the total rated power. This enables comparisons between results.
First, statistical characteristics of modeled turbine power output time series from WRFgenerated wind speeds are compared with those of measured data. In the subsequent subsection, the modeled turbine power output time series are used to optimize wind farm capacities at selective wind farm site locations, with the objective of minimizing the highfrequency fluctuations of the aggregated time series. Limitations of the spatial wind farm smoothing effect in the small island system are discussed in the third subsection. An optimization of the geographical distribution of wind farms is presented in the fourth subsection. The last subsection includes a discussion on the potential optimization effects for minimizing the overall power output time series for a variety of frequency ranges.
4.1 Comparing statistical characteristics of WRF modeled and measured dataset
To examine the validity of the modeled turbine power output data based on WRFgenerated wind speeds, these time series are compared to actual measurements. Three wind farms were operating on the Faroe Islands during the period of July 2016 to June 2018, the period of the WRFgenerated data. One additional measured time series is gained by converting meteorological data into turbine output data.
PSD, pdf of hourly stepchange functions, and duration curves for sitespecific and lumped wind farm power output data are calculated and depicted in Fig. 3; the lumped time series weights the installed capacity at the four site locations equally. The sitespecific PSD and pdf display similar behaviors as observed in the measurements. However, the PSD amplitudes at the highest frequencies are somewhat overestimated at sites B and C. This is also clear from the 5th and 95th percentiles and the standard deviations of the hourly stepchange functions presented in Table 1. The duration curves have evident sitespecific deviations.
4.2 Optimization of wind farm capacities
By optimizing the wind power capacity of scattered wind farms with the aim of minimizing the fluctuations of the total time series, the most stable wind farm configuration is derived. For this purpose, modeled wind turbine output data from predefined locations are used – the favorable wind farm locations in the region as selected by Magnussen (2017); see Sect. 2.3. The objective of the optimization is to minimize the 2–3 h fluctuations; see Sect. 3.3. The optimization algorithm is conducted for the following three cases:

Case I – all 14 favorable wind farm site locations mapped in Fig. 2 are considered;

Case II – the two southernmost wind farm site locations are excluded;

Case III – upper and lower boundaries for the installed wind power capacity are given for each of the 14 wind farms.^{1}
The optimized wind farm configurations for the three cases are presented in Fig. 4. The distributions of the optimized wind farm capacities are logical, with generally more capacities at more remote sites and fewer capacities for closely clustered wind farm sites. However, the effect of the optimization is limited compared with if the wind power capacities were to be equally distributed over all considered sites. The reduction in the 5th and 95th percentiles and standard deviations of the hourly stepchange functions is 2 % or less for the three optimized cases. It may be speculated that the observed limited hourly smoothing effect is because the improvement that the optimization adds is small compared with the smoothing that already occurred by considering the equally weighted spatial distribution of the wind farms.
To investigate this further, small wind farm portfolios in the region are examined by looking at combinations of three wind farms. For the combination of the southernmost site and the two clustered northernmost sites, optimized wind farm capacity distribution places around 27 % of the total power capacity of the wind farm portfolio at each of the northernmost sites and almost 47 % at the southernmost site. The distribution is logical, with lower capacities at both of the clustered wind farm sites and higher capacity at the distant site. If the two clustered sites were fully correlated on a 2–3 h scale while the distant site was uncorrelated, the capacity distribution would be expected to be 25 %, 25 %, and 50 %, respectively, which is close to what is observed for this setup. For the optimized portfolio with this setup, the 5th and 95th percentiles and the standard deviation of the hourly power output stepchange function are 4 %–5 % lower compared with if the wind power capacity were to be equally distributed over the three sites. Comparing the optimized wind farm portfolio to a portfolio in which 40 % of the total wind power capacity is placed at each of the two northernmost wind farms and 20 % is placed at the southernmost wind farm – a rather illogical distribution in terms of spatial smoothing of the power capacity – the optimized wind farm portfolio contains 12 %–13 % lower 5th and 95th percentiles and standard deviation. This indicates that the power capacity distribution of the wind farms matters.
Going back to Case I (all 14 favorable wind farm sites in the region), although the optimized portfolio has only limited effect on the smoothing effect compared with the portfolio in which the power capacities are equally distributed, comparing the overall power outputs of the optimized distribution of Case I to the power output of the equally distributed capacity of Case III with the constraint of the given boundary conditions, the 5th and 95th percentiles and the standard deviation of the hourly stepchange function are 4 % lower. Likewise, by comparing the power output of the optimized distribution of Case I to a setup where 11 % of the total wind power capacity is placed at the four closest clustered sites, as well as at the two northernmost sites, and 4.5 % at the rest of the sites (with the intent to distribute the capacity of the wind farms poorly), the 5th and 95th percentiles and the standard deviation of the hourly stepchange function are 6 % and 8 % lower, respectively. Again, this indicates that the wind power capacity distribution matters.
It is expected that the results are different for other case studies. The magnitude of the smoothing due to the spatial distribution of the wind farm capacities depends on the magnitude of the coherences and the details of the shape of the PSD curves. Other locations with, e.g., simpler terrain are expected to have larger coherencies. In contrast, larger regions have the advantage of being able to have scattered wind farms to a larger extent.
4.2.1 Comparison of optimized portfolios according to different conditions
PSD, pdf of the hourly stepchange functions, and the duration curves of the optimized lumped time series are displayed in Fig. 5. For reference, results from single turbine outputs are superimposed. The 5th and the 95th percentiles and the standard deviations of the stepchange functions are given in Table 2.
No clear distinction can be observed in the pdf of the stepchange functions and the duration curves of the three cases. However, a smoothing effect is observed compared with single turbine outputs with lower hourly stepchanges and less frequent intervals with zero and rated power production.
The spectra for all cases are equivalent for low frequencies and similar to those of the single turbine outputs, while the energy content at higher frequencies differs. A smoothing effect for the combined time series (Cases I–III) compared with the single turbine outputs is most evident at the highest frequencies but noticeable for periods of up to 1–2 d.
Out of the three cases, Case II (the exclusion of the two southernmost site locations) has the heaviest spectral tail, followed by Case III (bound by current and committed wind farm power capacities as distributed in the study by Tróndheim et al., 2021) and finally Case I (all site locations). The 5th and 95th percentiles and the standard deviations of the hourly stepchange functions of the time series show the same trend, with the highest values for Case II, followed by Case III and then Case I. Excluding the two southernmost site locations increases the 5th and 95th percentiles and the standard deviation of the hourly stepchange function by 7 %, while setting boundaries to the lower and upper wind farm capacities increases the percentiles and standard deviation by 3 %.
4.2.2 Sensitivity analysis on the PSD of power outputs with the inclusion of remote sites
As observed in Fig. 5 and Table 2, excluding the two time series of the southernmost island increases the optimized lumped time series fluctuation for the highest frequency range considerably. In other words, the inclusion of the two southernmost site locations smoothes the total time series.
To test if the smoothing effect originates from the inclusion of two additional sites in the given region or whether it is because the two sites are further away, PSDs for various optimized site combinations are generated and displayed in Fig. 6:

including all 14 sites (Case I);

excluding the two southernmost sites (Case II);

all possible combinations of 12 sites with the constraint that 2 out of the 12 sites are those located on the southernmost islands.
Case II has higher PSD amplitudes in the highest frequency range compared with other combinations of the 12 site locations that include the 2 southernmost site locations. This indicates that the smoothing effect for the highest frequencies is more pronounced for remote sites.
4.3 Limitation of spatially distributed wind farms
To test the limitations of the smoothing effect from spatially distributed wind farms in the Faroe Islands, PSDs are generated for multiple combinations of up to 14 turbine power output time series – the 14 favorable site locations marked in Fig. 2. Each computation considers an equal distribution of the rated power over all considered sites. Results are displayed in Fig. 7. As the number of lumped time series, n, increases, the PSD amplitudes at the highest frequencies decrease. The smoothing effect is evident for periods of up to 1–2 d but more evident for the higher frequencies. The smoothing effect observed when adding one additional wind farm to the lumped time series becomes less pronounced when n is large. These characteristics are similar to results observed by, e.g., Katzenstein et al. (2010) and Beyer et al. (1993), who analyzed the variability of interconnected wind power time series spatially dispersed in the area of Texas and northwest Germany, respectively.
The PSDs of the total power output when combining n wind farms differ for the various combinations of the individual wind farms. This will be analyzed further in the following subsection.
4.4 Optimization of wind farm positioning
As observed in Fig. 7, the same number of lumped wind farm time series can display different PSDs, depending on which individual time series are aggregated.
By extracting wind farm combinations with the smallest integrals of the PSD with respect to frequencies between (3 h)^{−1} and (2 h)^{−1}, the wind farm combinations with the lowest fluctuations are obtained. Figures 8 and 9 give examples of optimized wind farm combinations for n=4 and n=7, respectively. For reference, the worst wind farm combinations are also given – those with the highest integrals. It is observed that the wind farms appearing in the optimized combinations, i.e., with the smallest integrals, are scattered, while the opposite is observed for the worst combinations.
The 5th and 95th percentiles and the standard deviation of the hourly stepchange function of the lumped time series for each example are also given in Fig. 8 and 9. Both the percentiles and the standard deviations of the optimized combinations are considerably lower compared with the worst combinations, around 15 % and 18 %, respectively, for n=4 and around 13 % and 15 %, respectively, for n=7. In addition, the percentiles and standard deviation of the hourly stepchange function of the lumped time series for all combinations of n<14 wind farms are given in Table 3 together with the corresponding percentiles and standard deviation of the lumped time series with the smallest integral.
It can be concluded that the combinations of individual site locations when building n wind farms in the Faroe Islands have a considerable impact on the hourly fluctuations of the total power output time series – information that ought to be of interest to operators. Wind farm portfolios with distant sites are preferred in order to have the most stable lumped wind power time series.
4.5 Coherency between pairs of wind farm power output datasets at various frequency ranges
The focus of this wind farm optimization study has been to minimize the 2–3 h fluctuations of the overall power output time series. As argued in the introduction and observed in the results, optimized wind farm portfolios yield lower fluctuations at these frequencies. For longer periods, longer than a few days, the results in Sect. 4.3 established that there is a limited smoothing effect for any wind farm combinations in the region. The focus of this subsection is on the smoothing effect for frequencies in between these two.
High wind speeds are associated with high standard deviations of the time series, thus high spectral amplitudes, as the integral of the spectral amplitudes over frequency equals about half of the variance of the time series. As the spectral amplitudes are higher for the lower frequencies in the considered frequency range, it can be expected that the PSDs with the smallest spectral amplitudes at the lower frequencies are correlated with time series with smaller standard deviations and smaller wind speeds, thus smaller capacity factors. This is not desirable, as minimizing the fluctuations of the overall wind power output time series should not be at the cost of lower power production. Therefore, instead of looking at the PSD, this subsection will focus on the coherence function between pairs of wind turbine power output time series.
There is a connection between the coherency between time series and the smoothing effect of their combined time series. The higher the coherence, the lower the smoothing. For a coherence value of one, there is no smoothing of the combined time series, while a coherence of zero means that the time series are uncorrelated, and the fluctuations of the combined time series will be lower with respect to the produced power.
Figure 10 displays the squared coherence functions of the turbine power output time series for all possible site pairs out of the 14 favorable wind farm locations marked in Fig. 2 with respect to intersite distances. At the lower frequencies, the coherences are high, consistent with the fact that the same weather systems travel across the entire region. As the frequencies increase, the coherence values decrease, which is why in Fig. 7 in Sect. 4.3 it is observed that the smoothing effect from combining wind farms is more evident for higher frequencies. Another general characteristic that can be observed in Fig. 10 is that the closer the intersite distances are, the higher their coherence values are, which is consistent with the fact that the wind flow has had less time to change when traveling between two closer wind farm sites, which is why in Sect. 4.4 the optimized wind farm portfolio is a combination of distant wind farms.
To examine the potential smoothing effect on fluctuations at selective frequencies when combining wind farm power output time series, the average squared coherences between the selected frequency ranges are extracted from Fig. 10 and displayed in Fig. 11 with respect to the intersite distances between site pairs. Reasonably, the same general characteristics can be observed: higher coherence values for closer wind farms, higher coherence values at lower frequencies, and decreasing coherence values with increasing frequencies. The average squared coherence values for $(\mathrm{3}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}<f<(\mathrm{2}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}$ decrease rapidly with intersite distance and are already below 0.02 for all pairs of time series with intersite distances larger than 10 km. This means that in order to minimize the 2–3 h fluctuations of aggregated wind farm power output time series, all wind farms in the portfolio should be placed at least 10 km apart from each other. In order to do so in a small region and at the same time include several wind farms in the portfolio, most of the spatial area in the region has to be utilized. Although comparably less, the average squared coherence for $(\mathrm{9}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}<f<(\mathrm{6}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}$ also decreases rapidly with intersite distance, reaching a value of around 0.05 for an intersite distance of 30 km. This means that there are still possibilities for a maximum smoothing effect on the 6–9 h fluctuations for optimized wind farm portfolios consisting of a few wind farms, keeping in mind that the largest intersite distances are around 90 km.
The average squared coherence values for $(\mathrm{24}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}<f<(\mathrm{16}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}$ and $(\mathrm{72}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}<f<(\mathrm{48}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}$ decrease comparably more slowly with respect to intersite distance, and the values are scattered. The lowest average squared coherence values for $(\mathrm{24}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}<f<(\mathrm{16}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}$ and $(\mathrm{72}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}<f<(\mathrm{48}\phantom{\rule{0.125em}{0ex}}\mathrm{h}{)}^{\mathrm{1}}$ are 0.16 and 0.36, respectively. Although a smoothing effect at these frequencies should be observed for combinations of two time series with low coherence values, combining more than a couple of wind farms into a wind farm portfolio is assumed to result in a limited smoothing effect at these frequency ranges in the given region. Moreover, optimized wind farm portfolios consisting of more than a couple of wind farms are speculated to be determined by the standard deviation of the power output time serious at the expense of the total power production rather than by a combination of wind farms with low coherence values. However, the trend of a decrease in the average squared coherence values with intersite distance is noticeable. By extrapolating, it is expected that optimized wind farm portfolios consisting of multiple wind farms could have a high daily smoothing effect across larger regions without a reduction in the total power production.
The identification of the best distribution of wind farm production within a region has mainly been done in terms of identifying the energy yield. However, the present study discusses the use of knowledge on the underlying spatial–temporal characteristics of the governing wind field to assist in the search on system configurations for assuring a reduced variability in the power generation. The natural smoothing effect from spatially distributed wind farm sites in a small island system with a complex terrain is investigated in this study, using the Faroe Islands as a case study. For this, hourly modeled wind turbine power output data 45 m a.g.l. are analyzed. The modeled data are converted from WRFgenerated wind speeds at 14 favorable wind farm site locations as selected by Magnussen (2017). All results are presented per installed wind power capacity, i.e., normalized with respect to the rated power.
PSD, hourly stepchange functions, and duration curves are generated, and the 5th and 95th percentiles and the standard deviations of the hourly stepchange functions are extracted. As expected based on the literature, the smoothing from lumped power output time series is evident, with smaller highfrequency PSD amplitudes, lower hourly fluctuations, and fewer periods with zero and rated power production compared with single wind turbine outputs.
By optimally distributing wind farms, spatial smoothing can be maximized, resulting in less regulation needed from the operator. The focus of this study is on the smoothing effect on the 1–3 h scale. With the composition of the local power grid, it is considered that limiting the inherent wind variability on the hourly timescale is of local interest, especially as installed wind power is expected to increase rapidly in the future. In addition, it is argued that the hourly coherency between wind farm power outputs in a small region is expected to be dependent on how the regional weather travels between local sites, while, e.g., pairs of wind farm power output data are expected to be uncorrelated at the scale of seconds or minutes and close to fully coherent at the scale of days due to the same weather regime being present across the country. This study presents optimized wind farm portfolios for the Faroe Islands, where the objective function is set to be the integral of the PSD for frequencies between (3 h)^{−1} and (2 h)^{−1}, thus minimizing 2–3 h fluctuations.
Wind farm capacities at 14 predefined wind farm site locations are optimized. The results show that the wind farm capacities at remote sites should be higher and that the wind farm capacities at clustered sites should be lower. However, the optimization has only a small influence on the hourly fluctuations compared with if the wind farm capacities at the predefined wind farm site locations were equally distributed over all wind farms. The decrease in the 5th and 95th percentiles and standard deviation of the hourly stepchange function is about 2 %.
The optimization algorithm is performed for two additional cases: Case II – excluding two distant site locations, located ≥25 km apart from the rest of the sites; and Case III – setting upper and lower boundaries for the installed wind power capacity at each site. Also here, the optimization has only a small influence on the hourly fluctuations compared with if the wind farm capacities at the predefined wind farm site locations were equally distributed over all wind farms. However, comparing Case II to the first case, the 5th and 95th percentiles and standard deviation of the hourly stepchange function are increased by 7 %. The increase is found to be more pronounced when the two distant sites are excluded compared with if any other two site locations were to be excluded instead, demonstrating the importance of the smoothing effect from distant sites.
As the Faroe Islands consist of a limited spatial area surrounded by ocean far from any other land, the possibilities for the spatial distribution of wind farms are limited. To examine the achievable smoothing effect, PSDs are generated for up to 14 aggregated wind turbine time series. Results show less smoothing as more wind farms are integrated. The smoothing effect is most evident at the highest frequencies but noticeable for periods of up to 1–2 d. It is also seen that the highfrequency fluctuations are highly dependent on which of the individual site locations are considered. Optimized wind farm combinations, in order to minimize the 2–3 h power output fluctuations, are wind farms that are located distant from each other, while the worst wind farm combinations consist of clustered wind farms. When building wind farms in 4 out of the 14 favorable wind farm areas, optimized wind farm positioning decreases the 5th and 95th percentiles and the standard deviation of the hourly stepchange function by 15 % and 18 %, respectively, compared with the worst wind farm combinations.
By examining the coherency of the power output data between pairs of wind farms, it is concluded that a general characteristic in order to reach a maximum smoothing effect on the hourly scale is to place wind farms ≥10 km apart from each other. Although the focus of this study is primarily on the hourly scale, the coherency between pairs of wind farms focusing on a variety of lowerfrequency ranges is also presented and used to discuss the potential optimal smoothing effect for the respective frequencies. As hypothesized, coherency on the scale of days is high for all pairs of wind farms, implying that it is inefficient to conduct the optimization algorithm with a focus on minimizing dailyscale fluctuations. However, a decreasing trend in the daily coherency is observed with respect to distance between site pairs, indicating that an optimization algorithm minimizing dailyscale fluctuations could be applied for larger regions.
The importance of choosing the best wind farm sites is emphasized in this study in order to naturally balance wind power fluctuations. This feature should be of interest to an operator, as smoother time series result in a lower operating effort for the power grid.
This study exclusively considers wind power data 45 m a.g.l., being the hub height of all operating turbines in the Faroe Islands at the time of the preparation of this paper. However, current and future planned wind farms in the region consist of taller wind farms. In fact, six turbines with hub heights of 92 m a.g.l. have recently been installed and started producing electricity for the power grid. The presented work cannot be generalized without further calculations to wind power production at this or other heights, as wind is a function of height, and so are the time series fluctuations and their spatial coherencies. A suggestion for a future study is to include diverse wind turbine models in the analysis in order to investigate if and how these affect the results.
The WRFgenerated wind speed dataset used in this study displays artificial highfrequency fluctuations from the aliasing effect (Poulsen et al., 2021). Thus, before usage, these time series must be corrected: for each WRFgenerated wind speed time series, the power spectral density is generated and corrected for the aliasing effect using the method given in Kirchner (2005), with the parameters fc=1 and $f\mathit{\_}\mathrm{limit}=\frac{\mathrm{1}}{\mathrm{1}\mathrm{day}}$; see Fig. A1. By reversing the corrected spectral calculations and preserving the phase of the discrete Fourier transform of the original time series, a new dealiased time series is generated. Finally, the average wind speed value of the original time series is added to the dealiased time series. Figure A2 displays an example of a scatterplot between the WRFgenerated wind speed time series before and after dealiasing of the time series.
The results presented in this article are generated using the software MATLAB. The mfiles can be found here: https://doi.org/10.5281/zenodo.7317457 (Poulsen, 2022). In addition, for dealiasing of the WRFgenerated wind speed data, the mfile presented by Kirchner (2005, https://doi.org/10.1103/PhysRevE.71.066110) was used.
The data used in this paper are not publicly accessible but provided to the researchers at the University of the Faroe Islands for research purposes only by the local power company in the Faroe Islands, SEV. This includes both the WRFsimulated data and the measurements. The WRF data were simulated by Kjeller Vindteknikk with SEV as their client.
TP generated the results and wrote the initial manuscript. HGB supervised. All authors provided important input through discussion, feedback, and revision of the manuscript.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The authors are grateful to the local power company SEV for giving access to measurements and WRFgenerated data. We thank Kjeller Vindteknikk for their fast response to our questions regarding the WRFgenerated data. Turið Poulsen's PhD project was funded through the University of the Faroe Islands by research funding from Equinor.
This research has been supported in part by research funding from Equinor for Turið Poulsen through the University of the Faroe Islands.
This paper was edited by Mingming Zhang and reviewed by two anonymous referees.
Barasa, M. and Aganda, A.: Wind power variability of selected sites in Kenya and the impact to system operating reserve, Renew. Energ., 85, 464–471, https://doi.org/10.1016/j.renene.2015.05.042, 2016. a
Beyer, H. G., Luther, J., and SteinbergerWillms, R.: Fluctuations in the combined power output from geographically distributed grid coupled wind energy conversion systems – an analysis in the frequency domain, Wind Engineering, 14, 179–192, 1990. a, b, c, d
Beyer, H. G., Luther, J., and SteinbergerWillms, R.: Power fluctuations in spatially dispersed wind turbine systems, Sol. Energy, 50, 297–305, https://doi.org/10.1016/0038092X(93)90025J, 1993. a, b, c, d
Cassola, F., Burlando, M., Antonelli, M., and Ratto, C. F.: Optimization of the Regional Spatial Distribution of Wind Power Plants to Minimize the Variability of Wind Energy Input into Power Supply Systems, J. Appl. Meteor. Climatol., 47, 3099–3116, https://doi.org/10.1175/2008JAMC1886.1, 2008. a, b
Enercon: Enercon Product Overview, ENERCON GmbH, Aurich, Germany, April 2012, https://pdf.directindustry.com/pdf/enercon/enerconproductoverview/20877243513.html (last access: 16 October 2022), 2012. a
Frank, C., Fiedler, S., and Crewell, S.: Balancing potential of natural variability and extremes in photovoltaic and wind energy production for European countries, Renew. Energ., 163, 674–684, https://doi.org/10.1016/j.renene.2020.07.103, 2021. a
Giebel, G.: Equalizing effects of the wind energy production in Northern Europe determined from reanalysis data, Denmark, Forskningscenter Risoe, RisoeR, No. 1182(EN), https://backend.orbit.dtu.dk/ws/portalfiles/portal/7729010/ris_r_1182.pdf (last access: 16 October 2022), 2000. a, b, c
Giebel, G.: On the benefits of distributed generation of wind energy in Europe, PHD thesis from the Carl von Ossietzky Universität Oldenburg, Fortschr.Ber., VDI Reihe 6, Nr. 444, Düsseldorf, VDI Verlag 2001, ISBN 3183444062, 2001. a
Haslerud, A. S.: Faroe Islands – Detailed wind maps, availability limited to client, Kjeller Vindteknikk, Lillestrøm, Norway, Tech. Rep. KVT/ASH/2019/R032, 2019. a
Katzenstein, W., Fertig, E., and Apt, J.: The variability of interconnected wind plants, Energ. Policy, 38, 4400–4410, https://doi.org/10.1016/j.enpol.2010.03.069, 2010. a, b, c, d, e, f
Kirchner, J. W.: Aliasing in $\mathrm{1}/{f}^{\mathit{\alpha}}$ noise spectra: Origins, consequences, and remedies, Phys. Rev. E, 71, 066110, https://doi.org/10.1103/PhysRevE.71.066110, 2005. a, b, c
Magnussen, J. P.: Vindmyllustaðseting – val av økjum til vindmyllulundir í Føroyum & dømi um staðsetingar, Orkudeildin á Umhvørvisstovuni, Argir, Faroe Islands, Tech. Rep., 2017. a, b, c, d, e, f
Nanahara, T., Asari, M., Maejima, T., Sato, T., Yamaguchi, K., and Shibata, M.: Smoothing effects of distributed wind turbines, Part 2, Coherence among power output of distant wind turbines, Wind Energy, 7, 75–85, https://doi.org/10.1002/we.108, 2004a. a
Nanahara, T., Asari, M., Sato, T., Yamaguchi, K., Shibata, M., and Maejima, T.: Smoothing effects of distributed wind turbines, Part 1, Coherence and smoothing effects at a wind farm, Wind Energy, 7, 61–74, https://doi.org/10.1002/we.109, 2004b. a
Nørgaard, P. and Holttinen, H.: A MultiTurbine Power Curve Approach, proceedings of Nordic Wind Power Conference, Gothenburg, 1–2 March 2004, 1–5, http://lib.tkk.fi/Diss/2004/isbn9513864278/article4.pdf (last access: 16 October 2022), 2004. a
Palutikof, J., Cook, H., and Davies, T.: Effects of geographical dispersion on wind turbine performance in England: A simulation, Atmospheric Environment, Part A, General Topics, 24, 213–227, https://doi.org/10.1016/09601686(90)90458Y, 1990. a
Pearre, N. S. and Swan, L. G.: Spatial and geographic heterogeneity of wind turbine farms for temporally decoupled power output, Energy, 145, 417–429, https://doi.org/10.1016/j.energy.2018.01.019, 2018. a, b
Poulsen, T.: Code used to generate optimized wind farm portfolios for minimizing overall power fluctuations, Zenodo [code], https://doi.org/10.5281/zenodo.7317457, 2022. a
Poulsen, T. and Beyer, H. G.: Spectral characteristics and spatial smoothing of wind power – a case study of the Faroe Islands, Meteorol. Z., 29, 427–438, https://doi.org/10.1127/metz/2020/1041, 2020. a
Poulsen, T., Niclasen, B., Giebel, G., and Beyer, H. G.: Validation of WRF generated wind field in complex terrain, Meteorol. Z., 30, 413–428, https://doi.org/10.1127/metz/2021/1068, 2021. a, b, c
Reichenberg, L., Johnsson, F., and Odenberger, M.: Dampening variations in wind power generation – the effect of optimizing geographic location of generating sites, Wind Energy, 17, 1631–1643, https://doi.org/10.1002/we.1657, 2014. a, b
Schlez, W. and Infield, D.: Horizontal, Two Point Coherence for Separations Greater Than the Measurement Height, Bound.Lay. Meteorol., 87, 459–480, https://doi.org/10.1023/A:1000997610233, 1998. a
Sørensen, P., Cutululis, N. A., ViguerasRodríguez, A., Madsen, H., Pinson, P., Jensen, L. E., Hjerrild, J., and Donovan, M.: Modelling of power fluctuations from large offshore wind farms, Wind Energy, 11, 29–43, https://doi.org/10.1002/we.246, 2008. a
Tróndheim, H. M., Niclasen, B. A., Nielsen, T., Faria da Silva, F., and Bak, C. L.: 100 % Sustainable Electricity in the Faroe Islands: Expansion Planning through Economic Optimisation, IEEE Open Access Journal of Power and Energy, 8, 23–34, https://doi.org/10.1109/OAJPE.2021.3051917, 2021. a, b, c
ViguerasRodríguez, A., Sørensen, P., Viedma, A., Donovan, M., and Gómez Lázaro, E.: Spectral coherence model for power fluctuations in a wind farm, J. Wind Eng. Ind. Aerod., 102, 14–21, https://doi.org/10.1016/j.jweia.2011.12.006, 2012. a, b
Vincent, C., Larsén, X., Larsen, S., and Sørensen, P.: Crossspectra over the sea from observations and mesoscale modelling, Bound.Lay. Meteorol., 146, 297–318, https://doi.org/10.1007/s1054601297541, 2013. a
Wan, Y., Milligan, M., and Parsons, B.: Output Power Correlation Between Adjacent Wind Power Plants, J. Sol. Energ. Engin., 125, 551–555, https://doi.org/10.1115/1.1626127, 2003. a
Upper bounds are defined based on the sitespecific areas as marked in Magnussen (2017). For every 30 000 m^{2} (3⋅44 m$\cdot \mathrm{5}\cdot \mathrm{44}$ m ∼ 30 000 m^{2}), 0.9 MW is added to the possible upper boundary. The rotor diameter of an Enercon E44 turbine is 44 m, and it currently is the most frequently used turbine model on the Faroe Islands. However, none of the optimized values exceed the upper boundaries, making the upper boundaries redundant. The lower boundaries are set to the capacity of the current and committed wind farms on the Faroe Islands as distributed in the study by Tróndheim et al. (2021), with a total of 60.63 MW. Values for both the upper and lower boundaries are normalized by a total power capacity of 168 MW, being the proposed wind power capacity in 2030 by Tróndheim et al. (2021).