Next Article in Journal
Prospective of Upfront Nitrogen (N2) Removal in LNG Plants: Technical Communication
Next Article in Special Issue
Evaluation of Flicker of Light Generated by Arc Furnaces
Previous Article in Journal
Model-Based Contract Design for Low Energy Waste Heat Contracts: The Route to Pricing
Previous Article in Special Issue
Modeling Electricity Price and Quantity Uncertainty: An Application for Hedging with Forward Contracts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Wind Farm Cable Connection Layout Optimization with Several Substations

by
Adelaide Cerveira
1,*,
Eduardo J. Solteiro Pires
2 and
José Baptista
2
1
Department of Mathematics, School of Science and Technology of University of Trás-os-Montes and Alto Douro and INESC-TEC, UTAD’s Pole, 5000-801 Vila Real, Portugal
2
Department of Engineering, School of Science and Technology of University of Trás-os-Montes and Alto Douro and INESC-TEC, UTAD’s Pole, 5000-801 Vila Real, Portugal
*
Author to whom correspondence should be addressed.
Energies 2021, 14(12), 3615; https://doi.org/10.3390/en14123615
Submission received: 13 May 2021 / Revised: 10 June 2021 / Accepted: 12 June 2021 / Published: 17 June 2021
(This article belongs to the Collection Featured Papers in Electrical Power and Energy System)

Abstract

:
Green energy has become a media issue due to climate changes, and consequently, the population has become more aware of pollution. Wind farms are an essential energy production alternative to fossil energy. The incentive to produce wind energy was a government policy some decades ago to decrease carbon emissions. In recent decades, wind farms were formed by a substation and a couple of turbines. Nowadays, wind farms are designed with hundreds of turbines requiring more than one substation. This paper formulates an integer linear programming model to design wind farms’ cable layout with several turbines. The proposed model obtains the optimal solution considering different cable types, infrastructure costs, and energy losses. An additional constraint was considered to limit the number of cables that cross a walkway, i.e., the number of connections between a set of wind turbines and the remaining wind farm. Furthermore, considering a discrete set of possible turbine locations, the model allows identifying those that should be present in the optimal solution, thereby addressing the optimal location of the substation(s) in the wind farm. The paper illustrates solutions and the associated costs of two wind farms, with up to 102 turbines and three substations in the optimal solution, selected among sixteen possible places. The optimal solutions are obtained in a short time.

1. Introduction

1.1. General Context and Motivation

Climatic changes and the population’s awareness of pollution have contributed to policymakers changing their policies on energy production. On the other hand, we have seen an exponential increase in the world population resulting in high electricity consumption through the industrial and domestic load. Governments have cherished clean energies to the detriment of fossil energies that contribute to the planet’s pollution. Wind energy has benefited from these policies. Wind energy has been a considerable investment for almost all developed countries in the last two decades, with a strong bet on onshore wind farms. This commitment to strengthening renewable energies with a particular focus on wind power will continue in the coming decades. An example of this strategy is the Renewable Energy Directive, Directive (EU) 2018/2001, (RED II), which established a common framework for the promotion of energy from renewable sources in the EU and set a binding target of 32% for the overall share of energy from renewable sources in the EU’s gross final consumption of energy in 2030 [1]. In this context, the design of large wind farms using several substations is becoming increasingly common, bringing more efficient and more profitable exploration models to the new wind farm concessions.

1.2. Related Work

The design of wind turbines has been approached by many researchers using classical methods [2,3,4,5,6,7,8] and meta-heuristics [9,10,11,12,13]. However, works considering more than one station are scarce [14,15,16,17,18,19,20,21,22]. Moreover, they use meta-heuristics or optimization in two phases, which do not guarantee obtaining optimal solutions.
Chen et al. [14] use a fuzzy clustering algorithm to assign the turbines to a substation. The optimization of the cable layout is carried on using the minimum spanning tree algorithm. Considering meta-heuristics approaches and the optimization in two independent phases could lead to suboptimal solutions. In their work, a wind farm with three substations and 50 turbines was presented.
Wang et al. [15,23] use an integrated design method to minimize the total cost considering the substation location, connection topology, and cable cross-sections. The method uses an evolutionary algorithm to determine the substations coordinates, the substation associated with each turbine, and the cable layout. They applied the algorithm in a wind farm with two substations and up to 66 turbines.
Pillai et al. [16] use an approach for a wind farm’s cable design layout. To solve the problem, they divide it into two phases. The first phase places the substation using the capacitated clustering algorithm. In the second phase, for each subproblem, they use a mixed integer linear program (MILP) to find out the cable connection layout and cable type to be installed. The MILP takes into account initial solutions provided by some heuristics.
Huang et al. [17] propose a genetic algorithms to optimize the Horn Rev wind farm. They consider the transformers connected to wind turbines, substations costs, and cable costs.
Zuo et al. [18] analyzed the economic and reliability benefits of the collector system topology for offshore wind farms’ repowering and expansion. The method was based on cross-substation incorporation and radial topology to link repowered old wind farms and new ones. The multilayer optimization consists of an offshore substation (OS) refinement layer, an offshore wind farm partition layer, and an intra-zone cable connection layer. The collector system connection design was based on the predetermined WT selections and layouts, where the WTs were uniformly distributed by the substations. Moreover, the authors considered the repowered farms as one wind turbine. In each zone, the cable connection was based on the minimum spanning tree. A fuzzy clustering technique was used to determine the substation locations.
Dutta and Overbye [19] proposed a clustering-based algorithm for the cable layout of a large-scale wind power plant. Comparison of the proposed method with the radial feeder cable configuration shows that real power losses in the collector system are lowered and greater reliability is achieved with the proposed design. In this work, three different cable layout configurations have been compared.
Wu and Wang [21] use an ant colony algorithm to reduce wind farms’ construction costs and to improve the reliability of the collector system. They use the K-means clustering algorithm to partition the wind turbines into several groups to find the shorter collection lines in a radial configuration for the wind farm. They presented a problem with four substations and regions.
Pérez-Rúa et al. [22] propose a mixed-integer linear programming to optimize the cable system of offshore wind farms considering economic costs and losses power. The cable layout encompasses the interconnection between wind turbines (WTs) and transmission systems to couple Offshore Substations (OSSs) to the Onshore Connection Point. The work considers three case studies with up to 175 turbines and eight possible locations for substations, of which two are intended to be chosen. The computational time ranges from 1.48 h to 23.29 h.

1.3. Contribution and Document Structure

This work proposes an exact methodology to solve the optimization wind farm layout problem with several substations. Additionally, restrictions are used to limit connections on walkways. This paper’s main contribution is to efficiently solve the wind farm layout optimization, considering the topology and cable connection simultaneously, with several substations and many wind turbines. The proposed integer linear programming model determines the topology, selecting the connections and the substations that should be in the optimal solution and optimal cable to be used, minimizing energy losses and cable installation costs, all at the same time. The model also determines the optimal set of substations to be in the optimal solution.
The paper is divided into five sections. Section 2 describes the model of the electrical power grid underlying the wind farms. In Section 3, an integer linear programming model to optimize the wind farm layout considering several substations is presented. Section 4 presents and discusses the obtained results. Finally, Section 5 draws the main conclusions.

2. Electrical Grid Power Flow

The steady-state analysis of the power flow in the wind farm distribution network, where each wind turbine represents a node, assumes a fundamental role in this kind of research. It is essential to know all the network parameters to develop an equivalent model for all the constituent elements. Usually, this collector system is established with voltage levels between 20 kV and 30 kV. Since these networks have a radial structure, it is crucial to calculate the power flow circulating in each branch. Some research papers address exactly the problem of power flow in radial networks [5,24,25,26]. In addition, research such as [27], dedicated to the problem of reactive energy optimization in radial networks, present contributions to the calculation of the load flow also adapted to the radial feature of these networks.
As referred to in Reference [5], in the study of wind farms collecting system, the short line model must be used. This is due to the fact that the networks that connect the various elements of a wind farm are not longer than hundreds of meters or a few kilometers, where the R/X ratio is high. Considering the short line model, several simplifications can be made in the model, and the cables’ shunt admittance can be neglected. Therefore, the network branches can be represented by the model of Figure 1, where active ( P n ) and reactive ( Q n ) power flows between the buses are calculated by using Equations (1) and (2), respectively. The bus voltage is calculated using Equation (3) [28]. In this model, the turbines are represented as injector points of active and reactive power in the collecting system.
P n = P n + 1 + P wt n R n , n + 1 · P n + 1 2 + Q n + 1 2 | V n + 1 | 2 ,
Q n = Q n + 1 + Q wt n X n , n + 1 · P n + 1 2 + Q n + 1 2 | V n + 1 | 2 ,
| V n | 2 = | V n + 1 | 2 2 ( R n , n + 1 P n + 1 + X n , n + 1 Q n + 1 ) + ( R n , n + 1 2 + X n , n + 1 2 ) . P n + 1 2 + Q n + 1 2 | V n + 1 | 2 .
where R n , n + 1 and X n , n + 1 represents the branch resistance and reactance between buses n and n + 1 , respectively. Considering
I n + 1 2 = P n + 1 2 + Q n + 1 2 | V n + 1 | 2 ,
the active and reactive branch losses are given by Equations (5) and (6), respectively.
P loss ( n , n + 1 ) = R n , n + 1 · I n + 1 2 ,
Q loss ( n , n + 1 ) = X n , n + 1 · I n + 1 2 .
Therefore, the total power losses, P loss total , can then be obtained by adding all losses from all line sections as given by
P loss total = k B P loss ( k , k + 1 ) + Q loss ( k , k + 1 )
where B is the set of the buses in the network.
The resistance R n , n + 1 and reactance X n , n + 1 values between buses depend on the cable characteristics, namely: resistance R, inductance L, and the maximum current I z that can support. This latter value limits the number of turbines connected through a cable. In fact, the rated current drawn by each turbine, with a rated power P r and a U interconnection grid voltage, is defined by:
I r = P r 3 · U · cos φ
where the value of cos φ is the turbines’ power factor. For example, assuming that P r = 2 MW and U = 20 kV, the rated current drawn by each turbine is I r = 57.735 A. The I r value restricts the cable type used in a connection depending on the number of downstream wind turbines.
Consider the wind farm example shown in Figure 2, with two substations, nodes 19 and 20, and eighteen wind turbines, nodes 1 to 18. This layout has four branch lines connected to each substation. The total current reaching a substation from a branch line is the sum of the currents drawn by all turbines connected through this branch. For example, the branch starting in connection ( 19 , 3 ) supports three wind turbines (including turbine 3), and so, the current crossing this cable is 3 × I r = 173.22 A. The branch starting in connection ( 20 , 12 ) supports five wind turbines (including turbine 12), and so, the current crossing this cable is 5 × I r = 288.68 A.
In this work, the available cable types are presented in Table 1. So, in the connection ( 19 , 3 ) of Figure 2 a cable type less than k = 2 cannot be used, which has I z = 207 A, and in connection ( 20 , 12 ) a cable type less than k = 5 cannot be used, which has I z = 313 A.
Among all the available cable types, the one with the highest I z value is k = 10 with I z 10 = 585 A. So, in the case that I r = 57.735 A, any branch line cannot have more than I z 10 I r = 585 57.7 = 10 wind turbines, where a denotes the maximum integer not greater than a.

3. Mathematical Formulation

This section presents an Integer Linear Programming (ILP) model to obtain the optimal wind farm layout considering the topology and cables connection simultaneously for wind farms with several substations.

3.1. Data Sets and Parameters

Consider the node-sets N = { 1 , , n } , corresponding to wind turbines location and S = { n + 1 , , n + s } corresponding to substation locations. The goal is to obtain the wind-farm connection layout, i.e., a spanning forest of a graph, G = ( V , A ) with V = N S and A = { ( i , j ) : i S , j N } { ( i , j ) : i N , j N , i j } where the root of each spanning trees is a node belonging to set S.
In the layout problem is considered the set of available cable types K = { 1 , , k } presented in Table 1. The maximum current intensity I z k that each cable type k K can support bounds the number of downstream wind turbines for this cable by I z k I r . So, the maximum number of wind turbines in any branch line is:
Q = max k K I z k I r .
Therefore, the number of downstream wind turbines that each connection ( i , j ) can support, being node i in the substation side, is:
Q ( i ) = Q , i S Q 1 , i N .
Another parameter to consider is the F parameter that represents the maximum number of connections that can be linked to the substation, called feeders.

3.2. Costs

The total cost is given by the sum of the costs of active losses, c p , and reactive losses, c q , during the expected wind farm lifetime, and the infrastructure cost, c I , which includes the cable costs and digging cost.
For all pairs of nodes ( i , j ) A it is known the distance i j , in meters, between them. By using Equations (5) and (6), the cost of the active and reactive losses, during the wind farm lifetime, in a connection ( i , j ) supporting t downstream turbines, using a three-phase cable of type k, are, respectively,
c p i j k t = 3 · h · c e p j i · R k 10 3 · t 2 · l f 2 · I r 2 ,
and
c q i j k t = 3 · h · c e q j i · ω · L k 10 6 · t 2 · l f 2 · I r 2
where h is the number of hours during the expected wind farm lifetime, c e p is the active energy cost, c e q is the reactive energy cost, and l f is the load factor which reflects the real operating conditions during the wind farm lifetime; ω is the angular frequency.
The infrastructure cost to make a cable connection ( i , j ) using a cable of a type k is
c I i j k = ( D + 3 · C k ) · i j ,
where D is the digging cost (EUR by meter) and 3 · C k is the three-phase cable cost (EUR by meter).
Following [3], a preprocessing calculus is performed to determine the optimal cable type, k i j t , for a connection ( i , j ) supporting the current of t downstream wind turbines,
k i j t = arg min k K : t · I r I z k c I i j k + c p i j k t + c q i j k t ,
and the correspondent cost,
T i j t = c I i j k i j t + c p i j k i j t t + c q i j k i j t t
is the minimum cost for the connection ( i , j ) with t downstream turbines.

3.3. Decision Variables

The decision variables are:
  • For all ( i , j ) A , binary variables x i j t taking value 1 if the nodes i and j are connected (being node i on the substation side) and supports the current of t downstream wind turbines (including the one located in j); otherwise, it takes value zero.
  • For all i S , binary variables w i taking value 1 if substation i is in the solution; otherwise, it takes value zero.

3.4. Cable Connection Layout Model

The ILP model to optimize the wind farm layout considering multiple substations, WFLMS, is given by:
min ( i , j ) A t = 1 Q ( i ) T i j t · x i j t
s u b j e c t   t o
i S j N t = 1 Q t · x i j t = n
j N t = 1 Q x i j t F · w i , i S
i V t = 1 Q ( i ) x i j t = 1 , j N
i V t = 1 Q ( i ) t · x i j t 1 = i N t = 1 Q ( i ) t · x j i t , j N
x i j t { 0 , 1 } , ( i , j ) A , t = 1 , , Q ( i )
w i { 0 , 1 } , i S
The objective function (15) minimizes the total cost layout. Equation (16) guarantees that the network connects n wind turbines. Constraints (17) impose the maximum number of feeders, F, that can enter each substation. Constraints (18) guarantee that each wind turbine j N has one incoming connection. Constraints (19) are the flow conservation constraints and guarantee that for each wind turbine j N , if there exists an incoming connection supporting t downstream wind turbines, then the outgoing connections from turbine j must support t 1 downstream wind turbines. Finally, constraints (20) and (21) are the variable domain constraints.

3.5. Bounding Connections

If it is desired to bound the number of connections between a set of turbines, B, to the remaining wind turbines and substations, the following constraint must be included,
j B i V , i B t = 1 Q ( i ) x i j t + j B i N , i B t = 1 Q ( j ) x j i t L
where L is the maximum number of links between set turbines, B, to the remaining wind turbines and substations.

3.6. Substations Selection

The model WFLMS enables one to determine also the set of substations to include in the solution. So, it can be used to determine the best substation locations as long as a discrete set of substation locations is available. Furthermore, if a maximum number, M, of substations to be present in the solution is predefined, only Constraint (23) needs to be included.
i S w i M .

3.7. General Comments

It is also possible to incorporate in model WFLMS the topology of the wind turbines farm, choosing a given number of turbines from a set of possible locations of wind turbines. This is done by enlarging the turbine’s set, N, to the set of possible locations of the turbines from which it is intended to choose n turbines. In terms of the optimization model, it is not necessary to change the model, it is only necessary to extend the data files.
To take into account the real landscape and forbidden zones, it is necessary to compute the distances between the nodes in order to contemplate the ground situation.

4. Results and Discussion of the Case Studies

This section presents and discusses the obtained results using the proposed model to optimize the layout of two wind farms, WF-102-S2 and WF-74-S3, generated according to examples found in [29,30]. The first one has two substations and 102 wind turbines, while the second one has 74 wind turbines and sixteen possible sites for a substation, of which a maximum of three are to be chosen.
In all cases, the ten cable types presented on Table 1 were considered. Additionally, it was considered: h = 24 × 365 × 20 as the number of hours during the expected wind farm lifetime, assuming that it is 20 years; c e p = 102.52 × 10 6 EUR/Wh as the cost of active energy; c e q = 51.26 × 10 6 EUR/Wh as the cost of reactive energy; l f = 0.35 as the load factor, which reflects the real operating conditions during the wind farm lifetime, and is the ratio between the generated current and the maximum current that can be generated; ω = 100 π rad/s as the angular frequency.
The optimization models were constructed using FICO Xpress Mosel (Xpress Mosel Version 4.8.0), and then they were solved with FICO Xpress Optimizer. Computations were performed on a computer Intel(R) Core(TM) i7-8550U CPU @ 1.80 GHz 1.99 GHz with 8GB RAM and 64 bits.
The results of each wind farm are analyzed in the following sections.

4.1. WF-102-S2 Wind Farm

The first case study is the WF-102-S2 wind farm with two substations and 102 wind turbines with P r = 2 MW of rated power, interconnected by a U = 20 kV grid. With these parameters, the rated current drawn by each turbine is I r = 57.735 A, and the maximum number of wind turbines per branch line is Q = 10 . This value is due to the fact that the maximum current intensity that the available cables can support is I z = 585 , Table 1, and so Q = 585 57.7 = 10 , Equation (9).
The coordinates of the wind turbines and substations are in Table A1.
Two scenarios are considered for this wind farm: the wind farm original, WF-102-S2, and the wind farm WF-102-S2W, which includes a limit in the sidewalk connection.
For the first scenario, WF-102-S2, the obtained cable connection layout is presented in Figure 3, the costs in the optimal solution are presented in Table 2, and information about optimal cable connections is presented in Table 3.
The correspondent model WFLMS for this scenario has 207 constraints and 94,758 variables, and the processing time to obtain the optimal solution was 79 s.
The total cost obtained is EUR 10,475,772.3, where 60.5% is the infrastructure cost, corresponding to EUR 6,332,484.5, 23.0% is the active losses cost, corresponding to EUR 2,410,659.4, and 16.5% is the reactive losses cost, corresponding to EUR 1,732,628.4. The highest amount corresponds to the infrastructure cost, and the smallest part is the reactive losses cost during wind farm lifetime. Table 3 shows information about the cable types being used in the optimal solution. In the optimal solution, only five different types of cables are used: type 3 in 41 connections, type 4 in 23 connections, type 7 in 18 connections, type 8 in 6 connections, and type 10 in 14 connections.
The second scenario, WF-102-S2W, is obtained by limiting to two the connections passing on the walkway, adding constraint (22) with L = 2 in the model WFLMS.
The optimal connection layout for WF-102-S2W is presented in Figure 4 and the costs in the optimal solution are presented in Table 2.
The correspondent model has 207 constraints and 94,758 variables, and the processing time to obtain the optimal solution was 85.3 s.
In the optimal solution, only five different types of cables are used: type 3 in 39 connections, type 4 in 22 connections, type 7 in 17 connections, type 8 in 10 connections, and type 10 in 14 connections, Table 3.
The total cost is EUR 10,667,045.8, where: 57.3% is the infrastructure cost, corresponding to EUR 6,112,964.7; 23.7% is the active losses cost, corresponding to EUR 2,527,016.9; and 19.0% is the reactive losses cost, corresponding to EUR 2,027,064.2.
There are two sectors in both optimal wind farm solutions, WF-102-S2 and WF-102-S2-w: one has an installed capacity of 100 MW, with 50 wind turbines connected to substation 0 1 and the other one has 104 MW of installed capacity, with 52 turbines connected to substation 0 2 . There are ten branch lines linked to each substation, and the number of wind turbines in each branch line ranges between four and eight turbines. In the optimal solutions, only five different types of cables are used 3, 4, 7, 8, and 10, as shown in Table 3.
Comparing both scenarios’ optimal solutions, in Table 2 and Table 3, as expected, the total cost corresponding to the optimal solution when limiting the number of connections to a subset of turbines, study case WF-102-S2-W, is higher. This phenomenon is due to the increasing of active and reactive loss costs that are not compensated by the decreasing observed in the infrastructure cost. Note that these changes in costs result from the use of lower types of cables, which are cheaper but have higher energy losses.

4.2. WF-74-S3 Wind Farm

The second wind farm, WF-74-S3, is formed by 74 turbines with P r = 2 MW of rated power, interconnected by a U = 20 kV grid. With these parameters, the rated current drawn by each turbine is I r = 57.735 A, and the maximum number of wind turbines per branch line is Q = 10 . Sixteen possible positions for substations, distributed in a grid of 4 × 4-type points over the wind farm, are considered: O 1 , O 2 , ..., O 16 . The coordinates of the wind turbines and substations are presented in Table A2.
The goal is to optimize the cable layout, choosing at most three of the available substations. To solve this problem, the model WFLMS including constraint (23) with M = 3 is considered.
The optimal connection layout is presented in Figure 5, the costs in the optimal solution are shown in Table 4, and information about optimal cable connections is presented in Table 5.
The model has 166 constraints and 60,474 variables, and the processing time to obtain the optimal solution was 158.4 s.
The optimal wind farm has three sectors: one has an installed capacity of 26 MW, with 18 wind turbines connected to the substation O 1 ; the other one has 50 MW of installed capacity, with 25 turbines connected to the substation O 8 ; and the last one has 31 wind turbines connected to the substation O 14 , having 62 MW.
There are seven branches connected to substation O 1 , in which the number of turbines varies from one to five. Substation O 8 is linked to eight branch lines; it also has a number of turbines ranging between one and five. There are nine branch lines linked to the substation O 14 , in which the number of turbines varies between one and six. In the optimal solutions, only five different types of cables are used: type 3 in 30 connections, type 4 in 21 connections, type 7 in 13 connections, type 8 in 6 connections, and type 10 in 4 connections, as shown in Table 5.
The total cost is EUR 2,908,787.57, where 64.9% is the infrastructure cost, corresponding to EUR 1,887,148.4, 22.7% is the active losses cost, corresponding to EUR 660,641.2, and 12.4% is the reactive losses cost, corresponding to EUR 360,997.9, Table 4.
Once again, the highest amount corresponds to the infrastructure cost and the lowest amount is the reactive losses cost during wind farm lifetime.
Comparing this approach with approaches in the literature, Reference [22], the presented model is much faster at finding the optimal solution.

5. Conclusions

This work presents an ILP model to solve a cable connection layout considering wind farms with several substations. The model was applied to two wind farms with up to 102 turbines. In the wind farm with 102 turbines, there are two substations, and two variants of this case study are considered: one in which the original model WFLMS is considered and another in which a new model is considered where, to meet the solution installed on the ground, the number of connections between a set of turbines and the rest of the wind farm is limited. An additional constraint was considered in the latter case study. The proposed model was able to determine the optimal solution in both variants in a very short time, around 80 s.
In a wind farm with 74 turbines, sixteen possible locations for substations are considered, and the goal is to optimize the layout of the wind farm with at most three substations. To contemplate this new issue, it is only necessary to include a new constraint (23), to the initially proposed model, WFLMS. The optimal solution for this case study was obtained in only 158.4 s.
To contemplate this new issue, a new constraint (23), to the initially proposed model, WFLMS, was included. The optimal solution for this case study was obtained in only 158.4 s.
The results are promising. The optimal solutions are obtained in a short time and, if desired, the number of possible locations for the substations can be increased, bringing the discretized problem closer to the problem of choosing the best locations globally across the park.
The proposed model can also consider forbidden zones and can be adapted to optimize also the topology of the wind turbines farm, choosing a given number of turbines from a set of possible locations of wind turbines. Furthermore, compared with the literature, the presented model is fast at finding the optimal solution.

Author Contributions

Conceptualization, A.C., E.J.S.P. and J.B.; Formal analysis, E.J.S.P.; Investigation, A.C., E.J.S.P. and J.B.; Methodology, A.C. and E.J.S.P.; Resources, J.B.; Validation, J.B.; Writing—review & editing, A.C., E.J.S.P. and J.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Wind Farm Coordinates

This section presents the wind farms’ coordinates WF-102-S2 and WF-74-S3 in Table A1 and Table A2. Coordinates for wind farms WF-102-S2 and WF-74-S3 are expressed in WGS84 and Cartesian Coordinates, respectively. Substations and wind turbines are labeled on column “No”. Columns “Latitude” and “Longitude” show the correspondent coordinates.
Table A1. WF-102-S2 wind farm coordinates (WGS84).
Table A1. WF-102-S2 wind farm coordinates (WGS84).
No.LatitudeLongitudeNo.LatitudeLongitudeNo.LatitudeLongitude
0 1 54.04470000−3.502550003454.08490000−3.574716676954.09626667−3.63493333
0 2 54.07951667−3.575116673554.08981667−3.585100007054.10030000−3.64470000
154.02993333−3.447316673654.10518333−3.617616677154.10380000−3.65460000
254.03440000−3.455883333754.10988333−3.627533337254.01888333−3.50275000
354.03886667−3.464466673854.11658333−3.641716677354.02335000−3.51133333
454.04331667−3.473050003954.02441667−3.475033337454.02781667−3.51991667
554.04778333−3.481633334054.02886667−3.483616677554.03226667−3.52850000
654.05225000−3.490216674154.03333333−3.492200007654.03673333−3.53708333
754.05671667−3.498800004254.03780000−3.500783337754.04118333−3.54568333
854.06116667−3.507400004354.04226667−3.509366677854.04565000−3.55426667
954.06563333−3.515983334454.04671667−3.517950007954.05010000−3.56286667
1054.07010000−3.524566674554.05118333−3.526550008054.05475000−3.57295000
1154.07455000−3.533166674654.05563333−3.535133338154.05955000−3.58366667
1254.07901667−3.541766674754.06010000−3.543733338254.06296667−3.59505000
1354.08361667−3.550650004854.06836667−3.562233338354.06450000−3.60870000
1454.08823333−3.559533334954.07218333−3.572150008454.06603333−3.62235000
1554.09283333−3.568433335054.07598333−3.582066678554.06756667−3.63601667
1654.09743333−3.577316675154.07980000−3.591966678654.09203333−3.64363333
1754.11125000−3.604000005254.08360000−3.601900008754.09580000−3.65346667
1854.11585000−3.612900005354.10071667−3.626233338854.01613333−3.51660000
1954.12045000−3.621800005454.11020000−3.648166678954.02058333−3.52520000
2054.12536667−3.631283335554.02165000−3.488883339054.02505000−3.53376667
2154.02716667−3.461166675654.02611667−3.497466679154.02950000−3.54235000
2254.03163333−3.469750005754.03058333−3.506066679254.03396667−3.55095000
2354.03610000−3.478333335854.03503333−3.514650009354.03841667−3.55953333
2454.04056667−3.486916675954.03950000−3.523233339454.04180000−3.56891667
2554.04503333−3.495500006054.04395000−3.531816679554.04603333−3.57950000
2654.04948333−3.504083336154.04841667−3.540400009654.04980000−3.59058333
2754.05395000−3.512666676254.05288333−3.549000009754.05308333−3.60208333
2854.05840000−3.521266676354.05733333−3.557583339854.05588333−3.61396667
2954.06286667−3.529866676454.06156667−3.567583339954.05816667−3.62616667
3054.06733333−3.538450006554.06580000−3.5775833310054.05991667−3.63860000
3154.07178333−3.547050006654.06933333−3.5883666710154.06475000−3.64645000
3254.07616667−3.556266676754.07296667−3.6009500010254.08780000−3.65233333
3354.08053333−3.565483336854.07660000−3.61353333
Table A2. WF-74-S3 wind farm coordinates (Cartesian).
Table A2. WF-74-S3 wind farm coordinates (Cartesian).
No.LatitudeLongitudeNo.LatitudeLongitudeNo.LatitudeLongitude
O 1 0639.10640.6151800.01410.0453000.03000.0
O 2 0639.11546.7161800.01812.0463000.03410.0
O 3 0639.12452.8171375.01812.0473000.03812.0
O 4 0639.13358.9181000.02225.0483800.03410.0
O 5 1542.20640.6190600.01410.0493800.03812.0
O 6 1542.21546.7200187.51000.0503800.01000.0
O 7 1542.22452.8210600.01812.0513000.01812.0
O 8 1542.23358.9220187.51410.0523375.01812.0
O 9 2445.30640.6230187.51812.0533375.02225.0
O 10 2445.31546.7240187.52225.0543375.02600.0
O 11 2445.32452.8251375.02225.0553800.02600.0
O 12 2445.33358.9261800.02225.0563800.00187.5
O 13 3348.40640.6271800.02600.0573800.01410.0
O 14 3348.41546.7281375.03000.0583800.03000.0
O 15 3348.42452.8291000.03000.0593375.03000.0
O 16 3348.43358.9300600.03000.0603000.02225.0
010187.50187.5310187.53000.0613000.02600.0
020187.50600.0320187.53410.0623000.01410.0
030600.01000.0331375.03410.0633000.01000.0
040600.00187.5341000.03410.0643000.00600.0
051000.00187.5351375.03812.0653800.00600.0
061000.00600.0361000.03812.0663375.00187.5
071375.01000.0370600.03812.0672600.02225.0
081375.00600.0381800.03410.0682187.52225.0
091800.00600.0391800.03812.0692187.51812.0
101800.00187.5402187.53812.0702600.01410.0
112600.00187.5412187.53410.0712187.51410.0
122187.50600.0422187.53000.0722600.01000.0
132187.51000.0432187.52600.0733375.01000.0
141800.01000.0442600.03000.0742600.01812.0

References

  1. Directive (EU) 2018/2001 on the promotion of the use of energy from renewable sources. Off. J. Eur. Union 2018, L 328, 209.
  2. Cerveira, A.; Baptista, J.; Solteiro Pires, E.J. Optimization Design in Wind Farm Distribution Network. In International Joint Conference SOCO’13-CISIS’13-ICEUTE’13; Herrero, Á., Baruque, B., Klett, F., Abraham, A., Snášel, V., de Carvalho, A.C., Bringas, P.G., Zelinka, I., Quintián, H., Corchado, E., Eds.; Springer International Publishing: Cham, Switzerland, 2014; pp. 109–119. [Google Scholar]
  3. Cerveira, A.; de Sousa, A.; Solteiro Pires, E.J.; Baptista, J. Optimal Cable Design of Wind Farms: The Infrastructure and Losses Cost Minimization Case. IEEE Trans. Power Syst. 2016, 31, 4319–4329. [Google Scholar] [CrossRef]
  4. Wędzik, A.; Siewierski, T.; Szypowski, M. A new method for simultaneous optimizing of wind farm’s network layout and cable cross-sections by MILP optimization. Appl. Energy 2016, 182, 525–538. [Google Scholar] [CrossRef]
  5. Cerveira, A.; Baptista, J.; Solteiro Pires, E.J. Wind farm distribution network optimization. Integr. Comput. Aided Eng. 2016, 23, 69–79. [Google Scholar] [CrossRef]
  6. Fischetti, M.; Pisinger, D. Mixed Integer Linear Programming for new trends in wind farm cable routing. Electron. Notes Discret. Math. 2018, 64, 115–124. [Google Scholar] [CrossRef] [Green Version]
  7. Bentz, C.; Costa, M.C.; Hertz, A.; Poirion, P.L. Cabling optimization of a windfarm and capacitated K-Steiner tree. In Proceedings of the PGMO-COPI’14 Gaspard Monge Program for Optimization—Conference on Optimization Practices in Industry, Palaiseau, France, 28–30 October 2014; p. 4. [Google Scholar]
  8. Hertz, A.; Marcotte, O.; Mdimagh, A.; Carreau, M.; Welt, F. Optimizing the Design of a Wind Farm Collection Network. INFOR Inf. Syst. Oper. Res. 2012, 50, 95–104. [Google Scholar] [CrossRef]
  9. Baptista, J.; Solteiro Pires, E.J.; Cerveira, A. Wind Farm Distribution Network Optimization Based in a Hierarchical GA. In Proceedings of the ISAP2013—Seventeenth International Conference on Intelligent System Applications to Power Systems, Tokyo, Japan, 1–4 July 2013; p. 6. [Google Scholar]
  10. Veeramachaneni, K.; Wagner, M.; O’Reilly, U.M.; Neumann, F. Optimizing energy output and layout costs for large wind farms using particle swarm optimization. In Proceedings of the 2012 IEEE Congress on Evolutionary Computation, Brisbane, Australia, 10–15 June 2012; pp. 1–7. [Google Scholar] [CrossRef]
  11. Qi, Y.; Hou, P.; Yang, L.; Yang, G. Simultaneous Optimisation of Cable Connection Schemes and Capacity for Offshore Wind Farms via a Modified Bat Algorithm. Appl. Sci. 2019, 9, 265. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, L.; Wu, J.; Han, R.; Wang, T. Minimizing Energy Loss by Coupling Optimization of Connection Topology and Cable Cross-Section in Offshore Wind Farm. Appl. Sci. 2019, 9, 3722. [Google Scholar] [CrossRef] [Green Version]
  13. Khanali, M.; Ahmadzadegan, S.; Omid, M.; Nasab, F. Optimizing layout of wind farm turbines using genetic algorithms in Tehran province, Iran. Int. J. Energy Environ. Eng. 2018, in press. [Google Scholar] [CrossRef] [Green Version]
  14. Srikakulapu, R.; U, V. Optimized design of collector topology for offshore wind farm based on ant colony optimization with multiple travelling salesman problem. J. Mod. Power Syst. Clean Energy 2018, 6, 1181–1192. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, L.; Wu, J.; Wang, T.; Han, R. An optimization method based on random fork tree coding for the electrical networks of offshore wind farms. Renew. Energy 2020, 147, 1340–1351. [Google Scholar] [CrossRef]
  16. Pillai, A.; Chick, J.; Johanning, L.; Khorasanchi, M.; de Laleu, V. Offshore wind farm electrical cable layout optimization. Eng. Optim. 2015, 47, 1689–1708. [Google Scholar] [CrossRef] [Green Version]
  17. Lingling, H.; Yang, F.; Xiaoming, G. Optimization of electrical connection scheme for large offshore wind farm with genetic algorithm. In Proceedings of the 2009 International Conference on Sustainable Power Generation and Supply, Nanjing, China, 6–7 April 2009; pp. 1–4. [Google Scholar] [CrossRef]
  18. Zuo, T.; Zhang, Y.; Meng, K.; Tong, Z.; Dong, Z.Y.; Fu, Y. Collector System Topology Design for Offshore Wind Farm’s Repowering and Expansion. IEEE Trans. Sustain. Energy 2021, 12, 847–859. [Google Scholar] [CrossRef]
  19. Dutta, S.; Overbye, T.J. A clustering based wind farm collector system cable layout design. In Proceedings of the 2011 IEEE Power and Energy Conference at Illinois, Urbana, IL, USA, 25–26 February 2011; pp. 1–6. [Google Scholar] [CrossRef]
  20. Dutta, S.; Overbye, T. Optimal Wind Farm Collector System Topology Design Considering Total Trenching Length. IEEE Trans. Sustain. Energy 2012, 3, 339–348. [Google Scholar] [CrossRef]
  21. Wu, Y.; Wang, Y. Collection line optimisation in wind farms using improved ant colony optimisation. Wind Eng. 2021, 45, 589–600. [Google Scholar] [CrossRef]
  22. Pérez-Rúa, J.A.; Stolpe, M.; Cutululis, N.A. Integrated Global Optimization Model for Electrical Cables in Offshore Wind Farms. IEEE Trans. Sustain. Energy 2020, 11, 1965–1974. [Google Scholar] [CrossRef]
  23. Wang, L.; Wu, J.; Tang, Z.; Wang, T. An Integration Optimization Method for Power Collection Systems of Offshore Wind Farms. Energies 2019, 12, 3965. [Google Scholar] [CrossRef] [Green Version]
  24. Chen, T.; Chen, M.; Hwang, K.; Kotas, P.; Chebli, E.A. Distribution system power flow analysis—A rigid approach. IEEE Trans. Power Deliv. 1991, 6, 1146–1152. [Google Scholar] [CrossRef]
  25. Moon, Y.H.; Choi, B.K.; Cho, B.H.; Kim, S.H.; Ha, B.N.; Lee, J.H. Fast and reliable distribution system load flow algorithm based on the Y/sub BUS/ formulation. In Proceedings of the 1999 IEEE Power Engineering Society Summer Meeting, Conference Proceedings (Cat. No.99CH36364), Edmonton, AB, Canada, 18–22 July 1999; Volume 1, pp. 238–242. [Google Scholar] [CrossRef]
  26. Gómez Expósito, A.; Romero Ramos, E. Reliable load flow technique for radial distribution networks. IEEE Trans. Power Syst. 1999, 14, 1063–1069. [Google Scholar] [CrossRef]
  27. Baran, M.E.; Wu, F.F. Optimal capacitor placement on radial distribution systems. IEEE Trans. Power Deliv. 1989, 4, 725–734. [Google Scholar] [CrossRef]
  28. Haque, M. Efficient load flow method for distribution systems with radial or mesh configuration. IEE Proc. Gener. Transm. Distrib. 1996, 143, 33–38. [Google Scholar] [CrossRef]
  29. Seafish. Kingfisher Wind Farms Chart: Walney 1–3 Offshore Wind Farm. 2021. Available online: https://kis-orca.org/downloads/ (accessed on March 2021).
  30. Pricope, R.; Neagu, B. Optimal Reconfiguration of a Wind Farm Power Distribution Network. Bul. Inst. Politeh. Din Iasi 2015, LXI, 91–101. [Google Scholar]
Figure 1. Line diagram of a radial distribution system.
Figure 1. Line diagram of a radial distribution system.
Energies 14 03615 g001
Figure 2. Wind farm layout example.
Figure 2. Wind farm layout example.
Energies 14 03615 g002
Figure 3. Optimal cable connection layout for WF-102-S2 wind farm.
Figure 3. Optimal cable connection layout for WF-102-S2 wind farm.
Energies 14 03615 g003
Figure 4. Optimal cable connection layout for WF-102-S2-W wind farm.
Figure 4. Optimal cable connection layout for WF-102-S2-W wind farm.
Energies 14 03615 g004
Figure 5. Optimal cable connection layout for the WF-74-S3 wind farm.
Figure 5. Optimal cable connection layout for the WF-74-S3 wind farm.
Energies 14 03615 g005
Table 1. Unipolar cable characteristics (LXHIOV) 18/30 kV.
Table 1. Unipolar cable characteristics (LXHIOV) 18/30 kV.
TypeSectionInductanceElectrical ResistanceMax. CurrentPrice
k(mm 2 )L (mH/km)R ( Ω / km) I z (A)C (EUR/m)
10500.620.641016906.80
20700.590.443020707.12
30950.570.320024707.98
41200.550.253028108.70
51500.540.206031312.77
61850.530.164035413.23
72400.500.125040814.89
83000.490.100045817.50
94000.470.077851921.09
105000.460.060558523.77
Table 2. Optimal costs: WF-102-S2 and WF-102-S2-W wind farm.
Table 2. Optimal costs: WF-102-S2 and WF-102-S2-W wind farm.
Case StudyCosts
C I C p C q Total
WF-102-S26,332,484.52,410,659.41,732,628.410,475,772.3
WF-102-S2-W6,112,964.72,527,016.92,027,064.210,667,045.8
Table 3. Optimal solution cable type: WF-102-S2 and WF-102-S2-W wind farm.
Table 3. Optimal solution cable type: WF-102-S2 and WF-102-S2-W wind farm.
Case Study# Feeders/SubstationConnections/Cable Type
O 1 O 2 347810
WF-102-S210103922171014
WF-102-S2-W1010412318614
Table 4. Optimal costs: WF-74-S3.
Table 4. Optimal costs: WF-74-S3.
Case StudyCosts
C I C p C q Total
WF-74-S31,887,148.4660,641.2360,997.92,908,787.5
Table 5. Optimal solution cable type: WF-74-S3 wind farm.
Table 5. Optimal solution cable type: WF-74-S3 wind farm.
Case Study# Feeders/Substation# Connections/Cable Type (k)
O 1 O 8 O 14 347810
WF-74-S378930211364
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cerveira, A.; Pires, E.J.S.; Baptista, J. Wind Farm Cable Connection Layout Optimization with Several Substations. Energies 2021, 14, 3615. https://doi.org/10.3390/en14123615

AMA Style

Cerveira A, Pires EJS, Baptista J. Wind Farm Cable Connection Layout Optimization with Several Substations. Energies. 2021; 14(12):3615. https://doi.org/10.3390/en14123615

Chicago/Turabian Style

Cerveira, Adelaide, Eduardo J. Solteiro Pires, and José Baptista. 2021. "Wind Farm Cable Connection Layout Optimization with Several Substations" Energies 14, no. 12: 3615. https://doi.org/10.3390/en14123615

APA Style

Cerveira, A., Pires, E. J. S., & Baptista, J. (2021). Wind Farm Cable Connection Layout Optimization with Several Substations. Energies, 14(12), 3615. https://doi.org/10.3390/en14123615

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop