Optimal heterogeneity in a simplified highly renewable European electricity system
The resource quality and the temporal generation pattern of variable renewable energy sources vary significantly across Europe. In this paper spatial distributions of renewable assets are explored which exploit this heterogeneity to lower the total system costs for a high level of renewable electricity in Europe. Several intuitive heuristic algorithms, optimal portfolio theory and a local search algorithm are used to find optimal distributions of renewable generation capacities that minimise the total costs of backup, transmission and renewable capacity simultaneously. Using current cost projections, an optimal heterogeneous distribution favours onshore wind, particularly in countries bordering the North Sea, which results in average electricity costs that are up to 11% lower than for a homogeneous reference distribution of renewables proportional to each country’s mean load. The reduction becomes even larger, namely 18%, once the transmission capacities are put to zero in the homogeneous reference distribution. Heuristic algorithms to distribute renewable capacity based on each country’s wind and solar capacity factors are shown to provide a satisfactory approximation to fully optimised renewable distributions, while maintaining the benefits of transparency and comprehensibility. The sensitivities of the results to changing costs of solar generation and gas supply as well as to the possible cross-sectoral usage of unavoidable curtailment energy are also examined.
keywords:large-scale integration of renewables, system design, renewable energy networks, wind power generation, solar power generation, levelised system cost of electricity, Europe
vresVRESVariable Renewable Energy Sources \newacronymreaREARenewable Energy Atlas \newacronymncepNCEPAmerican National Centers for Environmental Prediction \newacronymcfsrCFSRClimate Forecast System Reanalysis \newacronymacACAlternating Current \newacronymdcDCDirect Current \newacronymlcoeLCOELevelised Cost of Electricity \newacronymhvdcHVDCHigh Voltage Direct Current \newacronymisetIWESFraunhofer-Institut für Windenergie und Energiesystemtechnik \newacronymgaGAGenetic Algorithm \newacronymgasGASGreedy Axial Search \newacronymcsCSCuckoo Search \newacronymdeDEDifferential Evolution \newacronymCapExCapExcapital expenditures \newacronymOpExOpExoperational expenditures
The ambitious renewable energy targets set by European governments eu2050 imply that the share of renewables in electricity generation will increase significantly in the years to come. At present, the leading renewable technologies are wind, solar photovoltaics (PV) and hydroelectricity, of which only wind and solar PV have the potential for large scale expansion. The uneven distribution of wind and solar resources across the continent raises the question of how best to exploit these heterogeneous resources. If wind and solar generation capacities are concentrated in those countries with the best resources, this may increase demand for transmission and increase energy imbalances between countries; if wind and solar generation are distributed homogeneously, then the best renewable resources will not be fully used and total system costs may be higher than the heterogeneous optimum. In this paper, the consequences of heterogeneity for the whole electricity system, including backup generation and transmission, will be quantified.
Since wind and solar PV are both \glsvres, backup generation is needed if the electrical demand is to be met at all times. Backup generation introduces additional system costs, which depend on the mismatch between \glsvres generation and load. Using the degrees of freedom associated with the choice of the capacity distributions of \glsvres for each country, it is possible to smooth out the aggregated temporal generation pattern or even shape it towards the load pattern. As a result, the mismatch and thus the backup requirements is lowered. To decrease the dimensionality of the problem, renewable assets can be assigned homogeneously, proportional to the mean load of each country, with a uniform wind-to-solar mixing factor. This approach is demonstrated in Heide2010 ; Heide2011 , where optimal wind-to-solar mixes for Europe are found that minimise balancing and storage costs. Further reductions in backup requirements are possible by extending the transmission network to enable more energy exchange between the countries Rodriguez2014 ; sarah . The implications for total system costs of different homogeneous renewable penentrations, wind-solar mixes and transmission levels were considered in Sensitivity , where the cost-optimal design was found to consist of a renewable energy penetration of 50% and a wind fraction of 94%. Other relevant research on the advantages of grid extensions for the integration of renewables, including reduced variability and smaller forecast errors, can be found in Tradewind ; WE:WE410 ; Papae ; Schaber ; Schaber2 ; Egerer ; Brown .
In this paper the consequences of moving from a homogeneous spatial distribution of \glsvres and a uniform wind-to-solar mixing factor to a cost-optimal placement of \glsvres capacities around Europe are explored. The distribution of \glsvres plants is determined by at least two considerations. The first consideration is the geographical variation of the \glsvres quality. The resource quality is quantified through the capacity factor (CF) defined as
The capacity factor is a number between 0 and 1, where 0 means no generation and 1 means maximum generation at all times. Capacity factors for the European countries for onshore wind and solar PV are calculated using (1) and listed in Table 2. The second consideration is the geographical variation of the temporal generation pattern for a given \glsvres type. This effect is particularly important for wind since Europe is large compared to the correlation length of wind of 600 km npg-15-803-2008 ; Widen2011 ; 1748-9326-10-4-044004 , and wind therefore benefits from smoothing effects across the continent.
With these points in mind, the optimal heterogeneous spatial layouts of wind and solar PV across Europe is investigated and compared to the homogeneous layouts. The main point of comparison is the average cost of electricity, which is composed of the \glsvres, backup and transmission costs. Different approaches to cope with the resulting large number of degrees of freedom are considered. In the literature a common approach for heterogeneous systems is to use linear programming to optimise generation and transmission capacities simultaneously Czisch ; Scholz ; Hagspiel ; Schlachtberger2016 , but this has the drawback that only a selection of representative weather conditions can be considered before computation times become infeasible. This makes the results susceptible to over-tuning to the weather selection. Other groups have used genetic algorithms to optimise generation, storage and transmission over a full year in Australia Elliston and over three years in Europe, the Middle East and North Africa Bussar201440 . In this paper a novel local search algorithm was found to be most effective given the size and non-linear formulation of the optimisation problem, allowing 8 years of hourly weather to be considered.
A downside of pure optimisation approaches is that one loses an understanding of why particular solutions are optimal. This makes it hard to justify investment strategies to policy makers and to the public. To counter this downside, more intuitive heuristic methods are developed here to construct layouts based on knowledge of resource quality, which are then compared to layouts obtained through optimisation. Distributions proportional to capacity factors (similar to the approach in Schaber2 ) and distributions based on optimal portfolio theory that reduce risk, or standard deviation, of the in-feed (similar to approaches in HOL08 ; ROQ10 ; ROM11 ; THO16 ) are considered and compared.
This paper is organised as follows: Section 2 discusses the general modelling of the simplified European electricity system and the key infrastructure measures. Section 3 describes the construction of heterogeneous layouts. In Section 4 the performance of the different layouts and the resulting renewable penetrations for individual European countries are discussed. Section 5 contains an analysis of the sensitivity of the results to variations in component costs. We conclude the paper with a discussion on the results and an outlook on future research.
|Set of nodes|
|Mismatch (VRES generation minus load)|
|Generation of wind, solar or backup|
|Total renewable generation|
|Net power balance|
|Wind, solar or backup capacity|
|Transmission capacity for link|
|Power flow on link|
|CF||Wind/solar capacity factor|
|Average value of|
2 Methods I: general modeling
2.1 Renewable resource assessment
Realistic time series describing the country-specific wind and solar PV power generation and the load are the starting point of the advocated weather-driven modelling of a simplified networked European electricity system. The utilized data set has been released from the Fraunhofer Institute for Wind Energy and Energy System Technology (formerly ISET, now IWES) bofinger . This data set covers the eight-year period from January 2000 to December 2007, has a temporal resolution of one hour and a spatial resolution of km over all of Europe. Fixed country-specific capacity layouts have been used to first convert the weather data into onshore wind and solar power generation, and then to aggregate the latter over each of the 30 European countries; off-shore wind power generation is not considered. The country-specific load time series have been obtained from publicly available sources, extrapolated to cover missing data, and detrended from an annual growth of around 2% to their year 2007 values. For more specific details see bofinger ; Heide2010 . A good alternative description of the conversion modelling is given in AND15 ; STA16 ; PFE16 .
The obtained wind and solar PV power generation time series have been rescaled to the capacity factors (CFs) from 2014. The latter have been determined in accordance with equation 1 from the EuroStat data for the installed capacities and the total generation for the year 2014 eurostat1 ; eurostat2 ; barometer . The resulting CFs for each country and each technology are listed in Table 2. For some of the countries (particularly smaller countries) no data was available or the calculated result was too uncertain because of too little or no installed capacity. For these countries the CFs are calculated as an average value from surrounding countries. These cases are marked by a star. Some countries with an already high installed capacity have a relatively low capacity factor compared to validated results from AND15 . For them the CFs have been raised by a small factor: 8% in Germany, 4% for wind in Spain, 4% for solar in Italy and 2% for wind in Great Britain. The final capacity factors presented in Table 2 are in accordance with STA16 ; PFE16 , which presents a critical assessment of current and future national capacity factors. Capacity factors for wind are likely to rise further in the future because of re-powering of wind turbines with more efficient, modern turbines at higher hub heights irena2 .
2.2 The electricity network
The European electricity network is modelled as a simplified 30-node model, where each node represents a country. For each node the generation from \glsvres (see Table 1 for a summary of nomenclature),
can be expressed through two parameters. The penetration determines the amount of renewable energy generated relative to the mean load of the node,
while the mixing parameter fixes the wind-to-solar ratio,
Other forms of renewable power generation are neglected in this simplistic modelling approach.
The nodal difference between \glsvres generation and load
is called the mismatch. To avoid power outages, the demand must be met at all times. Since storage is not considered, any power deficits must be covered by backup generation. Dispatchable resources are not modelled explicitly, but are considered as part of the backup generation. If , excess energy must be curtailed, while if backup generation is needed. Together the two terms form the nodal balancing . It is possible to lower the balancing needs with transmission. Nodes with excess generation export energy , allowing nodes with an energy deficit to import energy to (partly) cover their energy deficit. The nodal injection, , is denoted . This leads to the nodal balancing equation,
The vector of nodal injections is called the injection pattern, and fullfills . The actual imports and exports, and thus the injection pattern, depend on the dispatch of the nodal balancing. The synchronised balancing scheme,
where all nodes are curtailing/generating backup synchronously (relative to ), fulfills two top priorities: it minimises the total backup generation for each time step and it minimises the overall backup capacity Rolando2015 . This stylised synchronised balancing scheme has also been chosen in view of the layout optimisation, since the computational time for an update step is much smaller than for other dispatch schemes, like for example the localised flow scheme used in two previous publications Rodriguez2014 ; sarah .
The linear relationship follows from the DC approximation, which is known to be a good approximation for high-voltage flows. For the Power Transfer Distribution Factors we have assumed unit susceptances Rolando2015 , allowing its construction from the Moore-Penrose pseudo inverse of the underlying network Laplacian.
2.3 Infrastructure measures
Following Sensitivity , the energy system cost is calculated based on a few key measures. Besides the cost of the \glsvres capacities, and , costs for the backup system and the transmission network are included. The backup system cost is split into two components, the cost of backup capacity and the cost of backup energy . The backup capacity cost covers expenses related to construction and to keeping the power plants online while the backup energy cost accounts for actual fuel costs. Expressed in units of the average annual load, the backup energy is given by
In principle, the backup capacity is fixed by a single extreme event. However with this definition, the results will be highly coupled to the particular data set used. To decrease the coupling, the 99% quantile is used rather than the maximum value,
where is the time sampled distribution of backup generation and . With this choice, the backup system will be able to fully cover the demand 99% of the time. The remaining 1% is assumed to be covered by unmodelled balancing initiatives, e.g. demand side management. Given the nodal values , the overall backup capacity
is calculated by summation.
In analogy, the transmission capacity is defined so that the flow is met 99% of the time. Transmission can be positive and negative, but since links are assumed bidirectional, only the magnitude (not the sign) of the flow is to be considered. Hence
where is the time sampled distribution of absolute flows and . Since the link length varies, is not calculated directly by summation, but instead as a weighted sum,
where denotes the length of link . Link lengths are estimated as the distance between the country capitals.
In this paper will be expressed in units of average annual load, in units of average hourly load and in units of average hourly load megametre.
2.4 Cost modelling
Cost assumptions for the elements of an electricity system vary greatly across the literature. In this study, the cost assumptions published by Sensitivity have been adapted with a single modification. The cost of solar has been reduced by 50% in accordance with near future solar PV panel price projections irena . The resulting estimates are listed in Table 3. In general, the cost assumptions are in the low end for \glsvres which reflects the expectation that the cost of \glsvres will go down in the future as the penetration increases. Backup generation is priced based on the cost of Combined Cycle Gas Turbines (CCGTs).
From the \glsvres penetration, the mixing factor and the mean load, the mean generation of each node can be calculated. Dividing by the associated capacity factor, the capacity is obtained. Except for transmission capacity, the present value of each element can be calculated directly as
where is the rate of return assumed to be 4% per year. The transmission capacity cannot be translated directly into cost as the cost depends on the length and the type of the link. Link costs are assumed to be 400€ per km per MW for AC links and 1,500€ per km per MW for HVDC links. For HVDC links, an additional cost of 150,000€ per MW per converter station pair (one at each end) is added McKinsey ; Schaber ; Schaber2 . The layout of AC and HVDC lines has been constructed by Rodriguez2014 according to the existing European network reported by ENTSO-E for the year 2011 ENTSOEcaps and new predicted lines until 2014 NorNed ; BritNed . It is shown in Figure 10.
To allow for comparison of different system layouts, the \glslcoe is a convenient measure Sensitivity ; UEC13 ; HIR15 . The \glslcoe is the cost that every generated unit of energy consumed during the lifetime of the project has to match the present value of investment Short1995 ,
Since the life time of the system elements differs, the \glslcoe is evaluated separately for each system element from each respective present value. The \glslcoe for the complete system is calculated by summation. Life times of 25 years for solar PV and onshore wind, 30 years for CCGT plants and 40 years for transmission infrastructure were assumed.
3 Methods II: heterogeneous layouts
The simplest way to distribute the renewable resources is to assign them homogeneously (relative to the mean load of the node) so that and . This homogenous layout is denoted as HOM. However this assignment might not be ideal since the capacity factors vary significantly between the nodes. Three heuristic schemes and a straightforward optimisation for the construction of heterogeneous layouts will be presented in the following four subsections. The naming of the distribution algorithms is summarised in Table 4.
3.1 Heuristic layout I: CF proportional (CFprop)
An intuitive first approach, called CFprop, is to assign resources proportional to the CF, or more general to the CF raised to an exponent . For a wind-only layout, the nodal renewable penetrations are given by
where is the overall penetration assumed to be 1. An equivalent expression for the solar-only layout is obtained by the substitution . Examples for = 1 are shown in Figure 1a for the wind- and solar-only layouts. In the layout illustrations, each bar represents a country.
CFprop layouts for any value of can be constructed as a linear combination of the wind and solar only layouts with
For practical reasons, it is not possible to realise extremely heterogeneous layouts. On the one hand the geographical potentials for \glsvres installations in countries with good renewable resources may be a limiting factor. On the other hand countries with poor renewable resources may not want to become too dependent on imports. To constrain heterogeneity, the heterogeneity parameter K is introduced by requiring
With this definition, K = 1 corresponds to a homogeneous layout while K = represents unconstrained heterogeneity. For the CFprop layouts, each value of K translates into an -dependent value of . For a given value of , the corresponding value is found by increasing until the first country violates equation 20. At the mix the values correspond to , respectively.
3.2 Heuristic layout II: extreme K-constrained (CFmax)
Although the overall capacity factor of a CFprop layout for 0 is higher than the capacity factor of the homogeneous layout, it is possible to achieve an even higher capacity factor without violating the constraints in equation 20. In the wind- and solar-only cases, the capacity factor is maximised by assigning = K to the countries with the highest capacity factor and to the remaining countries, except for a single in-between country which is fixed by the constraint
The wind- and solar-only cases of the CFmax layout constrained by are shown in Figure 1b. Similar to the CFprop layouts, the CFmax layouts for arbitrary values can be constructed as linear combinations (18) of the wind- and solar-only layouts.
3.3 Heuristic layout III: OPT
The optimal portfolio theory (OPT) is well known in mathematical finance MAR52 . It discusses different assets obtained from the tradeoff between maximizing their return and minimizing their risk. This concept has also been applied to find optimal deployment of wind and solar energy resources in large-scale energy systems ROQ10 ; ROM11 ; THO16 , where the overall capacity factor has been treated as the return and the variance of the renewable power generation as the risk. In modified form, we will use OPT to further explore \glsvres capacity layouts over Europe with low system cost of electricity.
The overall capacity factor of a wind-only () or solar-only () layout is defined as
represents the overall installed capacity. The overall capacity factor is a useful measure of return as this is of high importance for investors of renewable generation capacity. Investors seek to minimise the overall capacity investment, which corresponds to maximising the overall capacity factor.
OPT’s second measure is risk, for which we select the relative standard deviation of the overall mismatch
based on the country-specific mismatches (6). The smaller the risky standard deviation becomes the more likely is the reduced need for a backup infrastructure, which an investor tries to minimise HOL08 .
A possible heterogeneous wind- or solar-only capacity layout is sampled from a Monte Carlo procedure. The country-specific renewable penetrations are randomly and independently drawn from a Beta distribution
defined on the compact support (20). is the Gamma function. The two shape parameters and are determined by requiring and by envoking the maximum entropy principle HON98 to maximal smear out the Beta distribution over the interval (20). For the two parameters result in , , and for they are , . A capacity layout sampled with this procedure does not necessarily meet the requirement (21). For such cases, all are uniformly rescaled upwards or downwards until the requirement is fulfilled. During the rescaling some of the penetration parameters hit the constraint (20), and are then frozen for the remainder of the rescaling procedure.
The wind-only and solar-only portfolios for are shown in Figure 2 in blue and yellow, respectively, with the overall mismatch measure on the first axis and the overall capacity factor measure on the second axis. Each of the portfolios consists of 100000 layouts. Due to the elongated shape of the portfolios there is no clear extended Pareto front in the upper left corners. The Pareto front defines a line, for which at the same time the standard deviation of the overall mismatch (risk) can not be reduced further for a fixed overall capacity factor and the overall capacity factor (return) can not be increased further for a fixed standard deviation of the overall mismatch. For both portfolios we identify a single point to characterise minimum risk and maximum return. This is done by extracting a subset of the points which are simultaneously a part of the top 200 capacity factors and bottom 200 standard deviations. For this leaves a sample of 28 layouts for wind and 67 layouts for solar to average and to calculate the respective new overall capacity factor and new overall standard deviation. The resulting points are plotted in white on top of the portfolios. The layouts of these two Pareto optimal points are shown in Figure 1c.
In order to find an optimal combined layout, we interpolate between the Pareto-optimal wind-only and solar-only layouts according to (18) and (19). This interpolation conserves the constraint (21) and results in the line shown in Figure 2. Apparently some of the interpolated layouts are able to reduce the standard deviation of the global mismatch further. The interpolated layout marked with a black dot comes with the mixing parameter .
3.4 Optimised layouts
The full optimisation of the layouts is considered, with the objective to minimise the \glslcoe with respect to the 60 variables for the countries. Given the high dimensionality of the search space, a number of optimisation algorithms were tested including the Nelder-Mead method nelder , simulated annealing sa , genetic algorithms ga and cuckoo search cs . It was found that the continuous enforcement of the normalisation criterion (21) generally decreased the performance of the tested algorithms, and for that reason a new hybrid algorithm was developed to address this problem. While being a classical greedy algorithm in the sense that the locally optimal choice is always taken, the renormalisation problem was circumvented by moving only along the axial directions. The algorithm has been denoted Greedy Axial Search (GAS).
When a solution is renormalised, all values are scaled either up or down. Therefore, it is possible that some values end up outside the boundary (20). The values are fixed at the boundary and the rescaling is only applied to the remaining free values. In general this approach is problematic since it can change the direction of the search. This is circumvented by holding the specific value constant that is considered during the step up/down procedure along a given axis. In this way only some values are scaled down/up and the feasibility of moving up/down along the considered axis can be determined. This is the underlying principle of \glsgas.
As any greedy algorithm, the \glsgas algorithm works by taking the locally optimal choice. Hence the feasibility for each direction is evaluated, but only the best choice is accepted. This process is repeated until a convergence criterion is fulfilled. At this point the step size is reduced and the iterative optimisation procedure repeated until the step size drops below some tolerance. The algorithm structure is sketched in Algorithm 1. The StepUp and StepDown subroutines generate new solutions by stepping a solution (first argument) up/down along axis (second argument) with some step size (third argument) after which the solution is renormalised as described above. Values of , and were found to be appropriate.
All optimised layouts have been obtained using the GAS routine. These layouts will be denoted GAS layouts. Constraining the transmission and thereby reducing the transmission capacity can lead to an overall lower LCOE. This is discussed in Section 4.4. The layouts resulting from this additional optimisation will be denoted GAS* layouts.
|HOM||Homogeneous distribution proportional to the mean of each country’s load|
|CFprop||Distribution proportional to a power of the capacity factor|
|CFmax||Assignment to each country extremised within depending on|
|OPT||Distribution using Optimal Portfolio Theory|
|GAS||Distribution optimised using Greedy Axial Search algorithm|
|GAS*||As GAS, but with optimally constrained transmission|
|GASnoT||As GAS but with no transmission between countries, so that each country is self-sufficient at all times|
The optimal heuristic layouts CFprop, CFmax, OPT as well as the optimized layouts GAS will be discussed in the next three subsections, first for , then for , and finally for . The fourth subsection focuses on the transmission capacities.
By construction, the layouts CFprop, CFmax and OPT become identical and homogeneous for . Due to Eq. (20), their respective renewable penetrations are . Moreover, according to Eq. (19) their renewable mixes also turn out to be independent of the country index. For these strictly homogeneous layouts Figure 3 shows the dependence of the key infrastructure measures on as the blue curves. For the backup energy and backup capacity, the optimal mixing parameters are located around , which is slightly larger than the values found by Heide2010 ; Heide2011 . For the transmission capacity, the minimum occurs around . The main measure of interest, the \glslcoe, has a minimum at . The high cost at is caused by a combination of high backup energy/capacity costs and the fact that the CF of solar is generally lower than for onshore wind. The cost of producing one unit of energy is thus higher for solar than for onshore wind even though the specific \glsCapEx is lower for solar.
The homogeneous layout producing the minimum LCOE at is denoted as the ’HOM’ layout. It is illustrated in Figure 4a. Its total LCOE amounts to 59.7 €/MWh. The componentwise LCOE corresponding to the wind, solar, backup and transmission parts are listed in the third column of Table 5 and graphed as the second bar in Figure 5. Wind power dominates the overall LCOE. Its contribution amounts to 61%, and is followed by 21% from backup, 10% from solar and 8% from transmission.
Contrary to the HOM layout, the GAS layout is no longer strictly homogeneous. Of course, all renewable penetrations are still equal to , but as a result of the optimisation the wind-solar mixing parameters become heterogeneous. This is illustrated in Figure 4b. Two-thirds of the countries are wind-only with . The remaining countries have a significant share of solar. For some of those this was to be expected. Spain, Greece, Italy, Romania and Serbia have very large solar capacity factors. See again Table 2. However, other solar-rich countries, like Portugal, Bulgaria, Bosnia and Croatia, are not amongst them. Instead, Germany is also assigned a significant share of solar, although its solar capacity factor is only average. By taking a closer inspection of Table 2 we discover the following empirical finding for the GAS layout: all countries with come with a ratio between their solar and wind capacity factor which is smaller than . The countries with have a larger ratio , except for the three smallest countries Estonia, Latvia and Luxembourg.
Compared to the HOM layout, the -heterogeneity of the GAS layout is able to reduce the total LCOE by 3%. This is mostly a consequence of the reduced combined component costs for wind and solar power. Note, that the overall mixing parameter has also slightly reduced from 0.90 (HOM) to 0.84 (GAS). See the fourth column of Table 5 and the third bar of Figure 5. The costs for backup and transmission have not changed much; which is also apparent from the rightmost panel of Figure 3.
All layouts discussed so far include the transmission infrastructure. It is also interesting to compare them to an optimised layout without transmission. No exports and imports would then be possible and the injection pattern would always be zero. No transmission investment would be needed and the respective componentwise LCOE would be zero. However, the countries then have to balance their mismatches all by themselves, and this in turn requires more backup infrastructure with higher respective componentwise LCOE. For the GAS layout without the transmission infrastructure, which for clarity we denote as GASnoT, the total LCOE turns out to be 64.5 €/MWh. Compared to the HOM layout, the combined LCOE components for wind and solar power generation are almost the same, but the increase of the LCOE components for the backup power generation and capacity is significantly larger than the disappearance of the transmission component. See again Figure 3, Table 5 and Figure 5. The total LCOE of the GASnoT layout is 8% and 11.5% larger than for the HOM and GAS layout respectively. This clearly demonstrates the benefit of transmission Rodriguez2014 ; Sensitivity .
The GAS and GASnoT layouts are obtained from two independent optimisation efforts. This explains why the two layouts are actually quite different in the distribution of the wind and solar resources. Figure 4c illustrates the resulting wind-solar mixing parameters for the GASnoT layout. Contrary to the more extreme GAS layout, the majority of the countries comes with a mix below and well above . Only the most northern countries turn out to be wind-only. However, on average the mixing parameter for the GASnoT layout is again close to for the GAS layout.
More heterogeneity is introduced once is chosen to be larger than one. Figures 6a-c illustrate the optimal heuristic CFprop, CFmax and OPT layouts for . Their respective values are 0.86 – 0.87 (see Table 5), and have been fixed by minimising the LCOE (see Figure 3d). The general -dependence of the other infrastructure measures are illustrated in Figure 3a-c. The backup energies required for the three layouts are quasi identical, and no difference is seen to the HOM layout. Also the backup capacities are almost identical for the three layouts, and are slightly less than for the HOM layout. Differences are observed for the transmission capacities. The CFprop layout comes with the smallest transmission capacities, followed by the OPT layout. The CFmax layout has the largest transmission capacities because its heterogeneity is the largest. All layouts are found to have larger transmission capacities than the respective layouts.
The total LCOE of the three heuristic layouts are inbetween 57.2 – 57.9 €/MWh. See columns 6-8 in Table 5 and bars 5-7 in Figure 5. This is very close to the value 57.8 €/MWh found for the GAS layout. In this respect, the larger heterogeneity of the layouts do not represent a clear cost advantage when compared to the GAS layout, which is homogeneous in the renewable penetration parameters . The situation changes once the optimised GAS layout is considered, which is exemplified in Figure 6d. It exploits the wind resources over Europe in a more efficient way and reduces the wind component in the LCOE; consult column 9 of Table 5 and bar 8 in Figure 5. This reduces the total LCOE to 55.3 €/MWh.
The overall renewable penetration of the GAS layout is ; consult again Equation (21). However, the individual renewable penetration parameters now scatter within . As can be seen in Figure 6d, their distribution is extremely heterogenous. For half of the countries they are either or , and for the other countries just somewhere in-between. A more carefull inspection reveals an approximate heuristic law, which expresses the renewable penetration parameters
as a continuous and piece-wise linear function of an effective capacity factor
A least-square fit is shown in Figure 7.
The overall mixing parameter of the GAS layout is almost the same as for the GAS layout. Both layouts also have in common that 20 out of the 30 countries come with . The five largest of the countries with a non-zero solar component are also identical.
It is worth to take again a quick look at Figure 2. It shows that for the optimal CFprop, CFmax and OPT layouts have more or less the same close-to-minimum standard deviation of the overall mismatch (24) as the optimised GAS layout. This indicates that a minimised mismatch standard deviation serves as a good measure to determine an optimal infrastructure HOL08 . However, it is still a rough measure, since it does not allow to finetune the minimum-cost infrastructure.
For the GAS algorithm has more freedom to optimise the heterogeneous layout and to reduce the overall LCOE, see (20). The resulting layout is depicted in Figure 8. It has some similarity to the GAS layout, but of course the GAS layout is even more extreme. Its overall wind-solar mixing parameter is almost the same as for the counterpart. The overall cost reduction turns out to be small. As can be seen in Table 5, the total LCOE for the and GAS layouts are 55.3 and 54.5€/MWh, respectively. This small cost reduction is mainly caused by the opportunity to allocate more wind resources to the sites with a very high capacity factor, and it is weakened to some extend by slightly increased costs for the transmission component; compare column 14 with column 9 in Table 5.
Bulk results for the optimal heuristic layouts CFprop, CFmax and OPT are also listed in Table 5 and Figure 5. Their layouts are also found to be wind-dominated, with nearly the same values as for the respective GAS layout. The LCOE for these three heuristic layouts are larger than for the GAS layout. This of course was to be expected. However, their LCOE also turn out to be slightly larger than for the less heterogeneous GAS layout.
Another reason that the GAS optimisation might have been better than the heuristic layouts is that the GAS algorithm sees not just the capacity factors at each site, like the heuristic layouts, but also the geographical variation of the temporal generation pattern, which the GAS algorithm can exploit to shape the \glsvres generation pattern towards the load. However if this was the reason, the backup generation costs would have decreased from the heuristic to the GAS layout, which they do not. This suggests that the GAS optimisation’s success really lies with the free exploitation of capacity factors.
4.4 Transmission capacities
So far, only the total contribution of the transmission capacities to the overall LCOE have been discussed for various system layouts in Table 5 and Figure 5. Its geographic distribution has not yet been specified. This will be done in this Subsection, but not right away. At first we will investigate a procedure which further reduces the overall LCOE by reducing the transmission capacities to some extend.
The transmission capacities defined in Equation (13) have been derived from unconstrained power flows. They are determined by the most extreme flow events, which typically occur between countries with a large energy deficit and others with a large excess. These events are not expected to overlap with other extreme events when all countries face a large energy deficit. The latter determine the required backup capacities. Consequently, it can be expected that a modest reduction of the total transmission capacities will not, or at least not much, affect the total backup capacities and the total backup energy, and will lower the overall LCOE.
The objective is to minimise the expression in the first line, taking into account the two constraints of the second and third line. denotes the constrained transmission capacity of line . In the limit of unconstrained flows, where the second constraint can be discarded, the objective (28) can be rewritten as with the method of Lagrange multipliers and leads to the solution (8). For the following we will downscale the unconstrained transmission capacities from (13) by a uniform scaling parameter to obtain the constrained transmission capacities
Figure 9 illustrates the dependence of the LCOE on the transmission constraints by taking the unconstrained transmission capacities of the GAS layout and scaling them down by the uniform factor . At first, as decreases, the LCOE also decreases. A minimum is found at . For the and GAS layouts the minimum is found at the optimal values and , respectively. If the transmission capacities are downscaled further the LCOE starts to increase again due to increasing requirements for backup energy and backup capacity.
Table 5 lists also the modified GAS layouts resulting from the optimal scaling parameters. For clarity, we denote them as GAS layouts. Compared to the GAS layouts, the transmission contribution to the total LCOE is reduced and the backup contributions are slightly increased. The wind and solar components of the GAS and GAS layouts are of course identical. Compared to the GAS layout, the total LCOE of the GAS layout is reduced by 1.2€/MWh in absolute units and by 2.1% in relative units. For and the reductions are 2.7% and 2.8%, respectively. The reductions are also illustrated in Figure 5.
The geographic distribution of the transmission capacities for the GAS layout is shown in Figure 10. The transmission capacities are not homogeneously distributed across the network. By far the strongest links are attached to Spain and Great Britain, which are the two largest countries with severe renewable excess generation. Links to their second neighbours with big deficits in renewable power generation, in particular Germany and Italy, also turn out to be quite strong. The more expensive HVDC transmission lines are utilised less extensively.
5 Sensitivity analysis
5.1 Reduced solar cost
For the optimised GAS layouts as well as for the heuristic CFprop, CFmax and OPT layouts the optimal mixing parameter minimising the overall costs is located in the wind dominated region. This is a consequence of the substantially higher costs of solar generation compared to wind. The future price development of solar photovoltaic systems is rather uncertain. To analyse the sensitivity to future price drops in solar cost, we calculate optimised layouts for solar cost reductions of 25%, 50% and 75%. Cost reductions could come from improved production processes, or alternatively from increasing capacity factors. Based on data from AND15 , the capacity factor can be increased by up to 40% by applying dual axis tracking compared to the fixed position installation assumed in Table 2, which may offset the higher capital costs of such systems. In addition, studies on increasing the energy conversion efficiency are still being conducted. A recent study suggests a huge decrease in the total system cost of PVs in a far future system irena2 .
The resulting GAS portfolios are visualized in Figure 11. Not surprisingly we find that a decrease in solar cost leads to a continuously increase in totally installed solar capacity. This increase is not found to be equal at all nodes. The main solar electricity supplier, Spain, initially increases its solar capacity, but for the more extreme price reductions decreases it again. It seems more efficient to shift the production to other sites. Spain is the clear leader in terms of solar generation for large solar costs. However, in the case of 50% solar cost reductions Germany almost produces equal amounts as Spain. This might not appear to be intuitive from the figure as the renewable penetration of Germany is always smaller than for Spain, but the mean load of Germany is more than twice as large as the one of Spain. In the 75% scenario Germany passes Spain and becomes the main producer of solar power. In this most extreme scenario almost all countries deploy solar resources.
We illustrate the change in the associated LCOE due to the cost reductions in the cases of K = 1, 2 and 3 in Figure 12. For all cases of heterogeneity the associated total European LCOE drops steadily for decreasing solar costs. For a reduction of the solar cost by 25% the optimal mix is shifted from above to below this point, and the LCOE values drop by around 2€. As the solar cost is reduced by 50% the optimal mix drops further and lies between 0.6 and 0.7. For the LCOE is reduced by almost 5€ compared to the reference scenario. When reducing the cost of solar by 75%, solar becomes much cheaper than wind, and the optimal mix is shifted below , indicating a dominant share of solar. Compared to the reference scenario, the LCOE dropped by around 9€ for the case . We have to be aware that such large cost reductions for solar photovoltaic systems might not be plausible. A cost reduction is mostly to be expected from material and production costs but not from installation costs.
5.2 Increased backup cost
The future price developments of fossil fuels, which are likely to increase, will affect the cost of electricity. An increase in the cost of gas used by the CCGT generators leads to an increase in the variable operational expenses associated with backup generation. In principle this will also affect the structure of the optimised layouts, but we expect the structural change to be very small. As Figure 3a reveals, the mixing parameters -0.84 of the optimised GAS layouts also produce the minimum of the backup energy. Consequently, the structure of the layouts will more or less not change, but of course their LCOE will increase as the gas price increases. This increase is linear. For the GAS layout an increase in backup fuel price to 150% leads to a LCOE of 59.0€/MWh, which is an increase of 6.8%. An increase to 200% of the gas price results in a LCOE of 62.8€/MWh, which equals an increase of 13.6%.
The increased backup costs can to some degree be counterbalanced by the sale of curtailment energy. So far we have assumed that curtailed electricity is wasted renewable production. Selling the curtailment energy to other energy sectors like the heating and transportation sector is a promising possibility. The resulting decrease in LCOE depends on the selling price and the amount of electricity sold. Since we are discussing an all-European renewable penetration of throughout this paper, the total amount of curtailment energy is identical to the backup energy. Assuming to sell of it at a price of 80€/MWh, the LCOE of the GAS layout is reduced to 50.2€/MWh, which is a decrease of 9.2%. Note however, that the sale of curtailment energy might have a slightly bigger impact on the structural change of the optimised GAS layouts than increased backup costs. Since, again, the amount of curtailment energy is equal to the backup energy, Figure 3a also illustrates the dependence of the curtailment energy on the mixing parameter . For parameter values below -0.84 the curtailment energy increases strongly. Consequently, when taking the sale of curtailment energy into account, a proper layout optimisation will shift to some degree towards smaller mixing parameters.
5.3 Interpolations towards more and less heterogeneity
As the heterogeneity parameter changes from to 2 and 3, the LCOE of the optimised GAS layouts has decreased further; consult again Table 5 and Figure 5. It is quite natural to ask how much further the LCOE might decrease as gets even larger. The answer is shown in Figure 13. The LCOE decreases continuously with increasing heterogeneity. However, the benefit of increased heterogeneity becomes smaller and smaller. The increasing cost of transmission leads to a point where it is almost no longer economic beneficial to increase the heterogeneity. The LCOE of 54.5€/MWh for the GAS layout is already very close to the asymptotic value of 54€/MWh for very large .
On the contrary, it might be more politically correct to reduce the heterogeneity. If the optimised GAS layouts were to represent the minimum of a rather shallow cost landscape, then other, more homogeneous layouts could be found in their vicinity without increasing the LCOE too much. Unfortunately, the search space for the exploration is high-dimensional, 60-dimensional to be more precise, as each of the 30 countries comes with its two variables and . If for each variable we were to test two smaller and two bigger values around its GAS value, we would end up in testing layout explorations. This is infeasible. Instead, we explore simple one-parameter interpolations between the heterogenous GAS layouts and the homogeneous HOM layout:
The interpolation parameter is confined to . A value of represents the GAS layout while reproduces the homogeneous layout. Figure 14 illustrates the LCOE of the interpolated layouts. The dependence on turns out to be almost linear. It is only weakly convex. This might indicate that the cost landscape around the GAS minimum is not flat, and that it might not be possible to find more homogeneous layouts without increasing the LCOE too much.
6 Discussion and outlook
In this paper the heterogeneity of renewable resources in different countries has been explored, but the distribution of wind and solar capacities within each country was fixed. Further heterogeneity of renewables, particularly wind, could be exploited by fine-tuning the distribution of renewables within each country, or by using a finer-scale model of Europe that exposes the locations with high capacity factors. In a recent paper Hoersch2016 it was shown that the \glsvres costs in a heterogenous optimisation are up to 10% lower when using a 362 node model of Europe compared to a one-node-per-country model with 37 nodes, because the better exploitation of good sites offsets the increased exposure of grid bottlenecks within each country.
Only three generation technologies were considered here: solar PV, onshore wind and natural gas. The inclusion of offshore wind might not improve system costs, given its high LCOE, but the LCOE may be offset by the system benefit of its steadier feed-in profile. In addition, offshore offers other benefits compared to onshore wind which are not accounted for by the cost optimisation, such as higher rates of public acceptance. Given that offshore wind is geographically concentrated along the coastlines of countries, a finer-resolution grid model would be advisable to fully assess the integration of offshore wind.
Modelling hydroelectricity, which already supplies 17% of Europe’s electricity, would reduce the costs of backup energy and provide extra flexibility to integrate the \glsvres. Similarly, the incorporation of storage or the use of flexibility from the electrification of transport and heating may alow \glsvres to be balanced more locally, favouring homogeneous solutions.
Finally, while the cost reduction is a strong argument for a heterogeneous \glsvres layout, the realisation might be a political challenge. Since the optimal placing of resources was derived from a system perspective, a realisation would require full collaboration from all countries. Countries with low capacity factors would no longer be self sufficient, while countries with high shares of renewables, such as the countries bordering the North Sea with good wind sites, may encounter problems finding enough sites or with public acceptance.
An unequal distribution of renewable energy generation also raises the question of who should pay for the generation and transmission assets. Current market conditions do not allow renewable generators to recover their capital costs from the energy-only market, forcing countries to subsidise the expansion of renewables. A highly heterogeneous system would therefore require a system for countries to compensate each other for their renewable imbalances. Recent work on the allocation of network flows to users in highly renewable networks Brown2014 ; Tranberg may provide the basis for an equitable distribution of such costs in a highly heterogeneous system.
In this paper the cost-optimal spatial distribution of \glsvres in a simplified European electricity system has been investigated for the case where the mean \glsvres generation equals the mean load (). A heterogenous distribution of wind and solar capacities has been shown to result in an average electricity cost that is up to 11% lower than a homogeneous distribution of renewables proportional to each country’s mean load. This is because the capital costs of wind and solar dominate the total system costs, and allowing the system to build more \glsvres in countries with better capacity factors means that fewer wind turbines and solar panels need to be built in order to produce the same amount of energy.
If the heterogeneity parameter , which controls the maximum and minimum levels of renewables generation in each country compared to its mean load, is gradually relaxed from (homogeneous) to larger values (heterogenous) then there is a clear trend of cost reduction, which is steepest for smaller values of and flattens out above . This has the important policy consequence that Europe can profit from the benefits of heterogeneity without allowing renewable imbalances between countries to become excessive.
The optimal mixing parameter between wind and solar is remarkably robust as the heterogeneity is increased, favouring a high proportion of wind of between 80% and 90% in the \glsvres mix. The mixing parameter is, however, sensitive to the relative capital costs of wind and solar, dropping to between 60% and 70% as solar capital costs are decreased by 50% compared to the default cost assumption.
While the best results in terms of low total system costs have been obtained here by explicit optimisation, heuristic methods for heterogeneously distributing wind and solar capacities, based for example on capacity factors, produce results that have costs only a few percent higher than the optimal systems. Given the increased comprehensibility and transparency that heuristic methods provide, this may be a price worth paying for policy makers.
Tom Brown is funded by the CoNDyNet project, which is supported by the German Federal Ministry of Education and Research under grant no. 03SF0472C. The responsibility for the contents lies solely with the authors.
-  European Commission. A roadmap for moving to a competitive low carbon economy in 2050. Technical report, EC, March 2011.
-  Heide, D., Von Bremen, L., Greiner, M., Hoffmann, C., Speckmann, M., and Bofinger, S. Seasonal optimal mix of wind and solar power in a future, highly renewable Europe. Renewable Energy, 35:2483–2489, 2010.
-  Heide, D., Greiner, M., Von Bremen, L., and Hoffmann, C. Reduced storage and balancing needs in a fully renewable European power system with excess wind and solar power generation. Renewable Energy, 36:2515–2523, 2011.
-  Rodriguez, R.A., Becker, S., Andresen, G., Heide, D., and Greiner, M. Transmission needs across a fully renewable European power system. Renewable Energy, 63:467–476, 2014.
-  Becker, S., Rodríguez, R.A., Andresen, G.B., Schramm, S., and Greiner, M. Transmission grid extensions during the build-up of a fully renewable pan-European electricity supply. Energy, 64:404–418, 2014.
-  Rodríguez, R.A., Becker, S., and Greiner, M. Cost-optimal design of a simplified, highly renewable pan-European electricity system. Energy, 83:658 – 668, 2015.
-  Van Hulle, F., Tande, J.O., Uhlen, K., Korpås, M., Meibom, P., and …Völker, J. Integrating wind: Developing Europe’s power market for the large-scale integration of wind power. Technical report, European Wind Energy Association (EWEA), 2009.
-  Holttinen, H., Meibom, P., Orths, A., Lange, B., O’Malley, M., Tande, J.O., Estanqueiro, A., Gomez, E., Söder, L., Strbac, G., Smith, J.C., and Van Hulle, F. Impacts of large amounts of wind power on design and operation of power systems, results of IEA collaboration. Wind Energy, 14:179–192, 2011.
-  Papaemmanouil, A., Tuan, L.A., Andersson, G., Bertling, L., and Johnsson, F. A cost-benefit analysis of transmission network reinforcement driven by generation capacity expansion. In Power and Energy Society General Meeting. IEEE, 2010.
-  Schaber, K., Steinke, F., and Hamacher, T. Transmission grid extensions for the integration of variable renewable energies in Europe: Who benefits where? Energy Policy, 43:123 – 135, 2012.
-  Schaber, K., Steinke, F., Mühlich, P., and Hamacher, T. Parametric study of variable renewable energy integration in Europe: Advantages and costs of transmission grid extensions. Energy Policy, 42:498–508, 2012.
-  Egerer, J., Lorenz, C., and Gerbaulet, C. European electricity grid infrastructure expansion in a 2050 context. In 10th International Conference on the European Energy Market. IEEE, 2013.
-  Brown, T., Schierhorn, P., Tröster, E., and Ackermann, T. Optimising the European transmission system for 77% renewable electricity by 2030. IET Renewable Power Generation, 10:3–9, 2016.
-  Kiss, P. and Jánosi, I.M. Limitations of wind power availability over Europe: a conceptual study. Nonlinear Processes in Geophysics, 15:803–813, 2008.
-  Widen, J. Correlations between large-scale solar and wind power in a future scenario for Sweden. IEEE Transactions on Sustainable Energy, 2:177–184, 2011.
-  St. Martin, C.M., Lundquist, J.K., and Handschy, M.A. Variability of interconnected wind plants: correlation length and its dependence on variability time scale. Environmental Research Letters, 10:044004, 2015.
-  Czisch, G. Szenarien zur zukünftigen Stromversorgung. PhD thesis, Universität Kassel, 2005.
-  Scholz, Y. Renewable energy based electricity supply at low costs - Development of the REMix model and application for Europe. PhD thesis, Universität Stuttgart, 2012.
-  Hagspiel, S., Jägemann, C., Lindenburger, D., Brown, T., Cherevatskiy, S., and Tröster, E. Cost-optimal power system extension under flow-based market coupling. Energy, 66:654–666, 2014.
-  Schlachtberger, D., Brown, T., Schramm, S., and Greiner, M. The price of non-cooperation in a highly renewable European electricity system. In Proceedings of SDEWES 2016, 2016.
-  Elliston, B., MacGill, I., and Diesendorf, M. Least cost 100% renewable electricity scenarios in the Australian National Electricity Market. Energy Policy, 59:270 – 282, 2013.
-  Bussar, C., Moos, M., Alvarez, R., Wolf, P., Thien, T., Chen, H., Cai, Z., Leuthold, M., Sauer, D.U., and Moser, A. Optimal allocation and capacity of energy storage systems in a future European power system with 100% renewable energy generation. Energy Procedia, 46:40 – 47, 2014.
-  Holttinen, H., Milligan, M., Kirby, B., Acker, T., Neimane, V., and Molinski, T. Using standard deviation as a measure of increased operational reserve requirement for wind power. Wind Engineering, 32:355–378, 2008.
-  Roques, F., Hiroux, C., and Saguan, M. Optimal wind power deployment in Europe - a portfolio approach. Energy Policy, 38:3245–3256, 2010.
-  Rombauts, Y., Delarue, E., and D’haeseleer, W. Optimal-portfolio-theory-based allocation of wind power - taking into account cross-border transmission-capacity constraints. Renewable Energy, 36:2374–2387, 2011.
-  Thomaidas, N.S., Santos-Alamillos, F.J., Pozo-Vazquez, D., and Usaola-Garcia, J. Optimal management of wind and solar energy resources. Computers & Operations Research, 66:284–291, 2016.
-  Bofinger, S., Von Bremen, L., Knorr, K., Lesch, K., Rohrig, K., Saint-Drenan, Y-M., and Speckmann, M. Raum-zeitliche Erzeugungsmuster von Wind- und Solarenergie in der UCTE-Region und deren Einfluss auf elektrische Transportnetze. Abschlussbericht, Institut für Solare Energieversorgungstechnik, ISET e.V., 2008.
-  Andresen, G.B., Sœndergaard, A.A., and Greiner, M. Validation of Danish wind time series from a new global renewable energy atlas for energy system analysis. Energy, 93:1074–1088, 2015.
-  Staffell, I. and Pfenninger, S. Using bias-corrected reanalysis to simulate current and future wind power output. Energy, 114:1224–1239, 2016.
-  Pfenninger, S. and Staffell, I. Long-term patterns of European PV output using 30 years of validated hourly reanalysis and satellite data. Energy, 114:1251–1265, 2016.
-  Eurostat, The Statistical Office of the European Union. Energy statistics: supply, transformation and consumption, Table nrg_105a, 2016.
-  Eurostat, The Statistical Office of the European Union. Infrastructure: electricity: annual data, Table nrg_113a, 2016.
-  EurObserv’ER. Photovoltaic Barometer, 2016.
-  International Renewable Energy Agency (IRENA). The Power to Change: Solar and Wind Cost Reduction Potential to 2025. Technical report, 2016.
-  Rodriguez, R.A., Dahl, M., Becker, S., and Greiner, M. Localized vs. synchronized exports across a highly renewable pan-European transmission network. Energy, Sustainability and Society, 5:21, 2015.
-  International Renewable Energy Agency (IRENA). Renewable energy technologies: Cost analysis series, Volume 1: Power Sector, Issue 4/5, Solar Photovoltaics. Technical report, 2012.
-  McKinsey. RoadMap 2050: A Practical Guide to a Prosperous, Low-Carbon Europe. Technical report, European Climate Foundation, 2010.
-  European Transmission System Operators. Indicative values for Net Transfer Capacities (NTC) in Continental Europe. Technical report, ENTSO-E, 2011.
-  Skog, J.E., Koreman, K., Pääjärvi, B., Worzyk, T., and Andersröd, T. The NorNed HVDC cable link: a power transmission highway between Norway and the Netherlands. Technical report, ABB, 2007.
-  BritNed Development Limited. BritNed Construction. http://www.britned.com/.
-  Ueckerdt, F., Hirth, L., Luderer, G., and Edenhofer, O. System LCOE: What are the costs of variable renewables? Energy, 63:61–75, 2013.
-  Hirth, L., Ueckerdt, F., and Edenhofer, O. Integration costs revisited – An economic framework for wind and solar variability. Renewable Energy, 74:925–939, 2015.
-  Short, W., Packey, D., and Holt, T. A Manual for the Economic Evaluation of Energy Efficiency and Renewable Energy Technologies. Technical report, National Renewable Energy Laboratory, 1995.
-  Markowitz, L. Portfolio selection. Journal of Finance, 7:77–91, 1958.
-  Honerkamp, J. Statistical Physics – An Advanced Approach with Applications. Springer Verlag Berlin Heidelberg, 1998.
-  Nelder, J. A. and Mead, R. A simplex method for function minimization. The Computer Journal, 7:308–313, 1965.
-  Kirkpatrick, S., Gelatt, C.D., and Vecchi, M.P. Optimization by simulated annealing. Science, 220:671–680, 1983.
-  Goldberg, D.E. Genetic Algorithms in Search, Optimization and Machine Learning. Addison-Wesley Longman Publishing Co., Inc., Boston, MA, USA, 1st edition, 1989.
-  Yang, X-S. and Deb, S. Cuckoo search via Lévy flights. In 2009 World Congress on Nature & Biologically Inspired Computing (NaBIC), pages 210–214, 2009.
-  Hörsch, J. and Brown, T. The role of spatial scale in joint optimisations of generation and transmission for european highly renewable scenarios. In Proceedings of 14th International Conference on the European Energy Market (EEM 2017), 2017.
-  Brown, T. Transmission network loading in Europe with high shares of renewables. IET Renewable Power Generation, 9:57–65, 2016.
-  Tranberg, B., Thomsen, A.B., Rodriguez, R.A., Andresen, G.B., Schäfer, M., and Greiner, M. Power flow tracing in a simplified highly renewable European electricity network. New Journal of Physics, 17:105002, 2016.