Next Article in Journal
Economic Analysis of Special Rate for Renewable Energy Based on the Design of an Optimized Model for Distributed Energy Resource Capacities in Buildings
Previous Article in Journal
Interval Load Flow for Uncertainty Consideration in Power Systems Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cable Connection Optimization for Heterogeneous Offshore Wind Farms via a Voronoi Diagram Based Adaptive Particle Swarm Optimization with Local Search

1
School of Computer Science, University of Electronic Science and Technology of China, Zhongshan Institute, Zhongshan 528402, China
2
School of Computer Science and Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
3
SEWPG European Innovation Center, 8000 Aarhus, Denmark
4
Department of Operation, University of Groningen, 9747 Groningen, The Netherlands
5
Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
6
Center of Electric Power and Energy, Department of Electrical Engineering, Technical University of Denmark, 2800 Lyngby, Denmark
7
School of Electrical Engineering and Telecommunications, University of New South Wales, Sydney, NSW 2052, Australia
*
Author to whom correspondence should be addressed.
Energies 2021, 14(3), 644; https://doi.org/10.3390/en14030644
Submission received: 16 December 2020 / Revised: 10 January 2021 / Accepted: 22 January 2021 / Published: 27 January 2021
(This article belongs to the Section A3: Wind, Wave and Tidal Energy)

Abstract

:
Offshore wind energy, as one of the featured rich renewable energy sources, is getting more and more attention. The cable connection layout has a significant impact on the economic performance of offshore wind farms. To make better use of the wind resources of a given sea area, a new method for optimal construction of offshore wind farms with different types of wind turbines has emerged in recent years. In such a wind farm, the capacities of wind turbines are not identical which brings new challenges for the cable connection layout optimization. In this work, an optimization model named CCLOP is proposed for such wind farms. The model incorporates both the cable capital cost and the cost of power losses associated with the cables in its objective function. To get an optimized result, a Voronoi diagram based adaptive particle swarm optimization with local search is proposed and applied. The simulation results show that the proposed method can help find a solution that is 12.74% outperformed than a benchmark.

Graphical Abstract

1. Introduction

With the increasing demand for green energy, the use of rich renewable energy sources such as wind energy has drawn a wide attention [1]. Compared with onshore wind farms, offshore wind farms (OWFs) possess superiorities involving higher wind speed, rich wind energy resources, less turbulence intensity, good visualization, noise-free to residents, etc. [2]. However, the immense investment, including capital expenditure (CAPEX) and the operational cost of OWFs, is much higher than that of onshore wind farms, which seriously hinders their commission and development [3]. Therefore, in order to obtain an OWF with high investment benefits, researchers are devoted to studying the cable connection layout optimization problem (CCLOP) of OWFs [4]. Based on the concept of the minimum spanning tree (MST) in graph theory, the greedy algorithm has been used to solve the CCLOP [5]. In [5], the whole wind farm was divided into several regions using the fuzzy c-means algorithm, then the Prim algorithm is used to solve the CCLOP for each area. Due to the complexity of the CCLOP, metaheuristic algorithms were also adopted to solve the CCLOP [6,7]. In [6], the optimization model of the CCLOP is established considering the voltage level, quantity, and location of the offshore substation (OS), then the model was solved by the genetic algorithm. In [7], a new model based on the concept of MST was proposed and solved by the adaptive particle swarm optimization (APSO) algorithm. The simulation result showed that the APSO outperformed the result obtained by the prime algorithm and the proposed dynamic MST algorithm. However, the influence of the cable energy losses cost (CELC) [8] on the cable connection layout is not considered in the above work.
With the increasing capacity and distance to shore of offshore wind farm, the CELC becomes increasingly significant for CCLOP [8]. To consider this factor, some studies have been conducted. In [9], a hybrid immune genetic algorithm was proposed to solve the CCLOP. The cost of more components such as transformers and substations was considered simultaneously with the cables’ costs and solved by an improved genetic algorithm in [10]. In addition to using heuristic and metaheuristic algorithms, mathematical programming was also applied in solving CCLOP considering CELC’s impact [11,12]. In [11], the ordinal optimization and mixed-integer programming method were combined to solve the CCLOP and the results demonstrated that combined method obtains a better solution than mixed-integer programming. The CCLOP was formulated as mixed-integer linear programming problem in [12]. It was designed to ensure an uncrossed cable connection layout so as to eliminate the extra cots on installation as well as operations and maintenance.
On the other hand, existing OWFs are always composed of the same type of wind turbine (WT) whereas several recent studies have shown the benefits of constructing an OWF with diverse types of WTs [13,14]. Therefore, with the willingness of maximizing the use of the resources of the sea area, reducing the levelized cost of energy of the OWF, a new type of OWF, namely, an offshore wind farm with different wind turbine types (OWFD) came into being and attracted the worldwide attention. Reference [15] is one of the earliest studies on solving the wind farm layout optimization problem for an OWFD where the WTs positions and the type on each position were optimized simultaneously to maximize the total power generation. Similarly, the WTs’ positions and the type based on the Jensen wake model to minimize the levelized cost of energy considering the cost of the foundation system, offshore electrical system, offshore transportation and installation in [16]. A 3D wake model for complex terrain was established in [17] while the WT micro-siting problem was solved by the greedy algorithm considering installing WTs with different hub heights.
However, the above literature only concerns the micro-siting optimization of the OWFD. To our knowledge, our work is the first to put forward a method for solving the CCLOP of the OWFD. The CCLOP of OWFD aims to optimize the cable connection while meets all the technical and engineering constraints, which is an intractable combinatorial optimization problem. The exact algorithms face the challenges of ensuring a feasible solution to such a problem [18]. In comparison, meta-heuristic algorithms find suboptimal solutions in an acceptable time, and have been adopted in many works [19,20]. Ref. [19] proposed a technique of coding the topology to a binary string and apply genetic algorithm to solve the CCLOP considering both middle and high-voltage alternating-current grids. Considering the wake effect and accurate power losses calculation method, ref. [20] proposed solved by an APSO to solve the CCLOP.
According to our previous work, the proximity relationship between WTs is important in obtaining an optimal cable connection layout. To give gain a better judgement of the the proximity relationship between the vertices (i.e., OS or WTs) and thus find a better solution, the Voronoi distance [21,22] is introduced to define the proximity relationship between vertices, and the two adjacent vertices are connected preferentially to improve the searching efficiency of the algorithm in this paper. Voronoi distance has achieved good results in vehicle routing problem optimization [23], and our work is the first to incorporate the Voronoi distance into the CCLOP model.
The main contribution of this work is summarized in the following:
(1)
Considering the influence of CELC on the economic performance, a CCLOP model of OWFD is established. To the best of our knowledge, this is the first work that proposes an optimization model for solving the CCLOP of OWFD.
(2)
To resolve the CCLOP, a novel V-APSO-LS algorithm is proposed, which adopts a K-ring shaped Voronoi neighbor (KSVN) and a local search strategy (LSS) to enhance the optimization ability of APSO.
(3)
The proposed algorithms are well verified via two case studies.
The remaining of this paper is organized as follows: Section 2 describes the CCLOP of OWFD. Section 3 provides a V-APSO-LS algorithm to solve the model. A wind farm with 72 WTs is chosen as the study case to demonstrate the proposed method and its strategies in Section 4. Some conclusions and future work on OWFD are given in Section 5.

2. Mathematical Models

In this section, the assumptions associated with the model are introduced at first. A new CCLOP model is specified which is established according to the characteristics of OWFD.

2.1. Assumptions

To make the paper more readable, some corrections have been made for the assumption section in the following.
The type of WT in each position within the wind farm is given which is also the truth that the CCLOP was solved after the micro-siting of wind turbines were decided.
To simply the cable length calculation, the distance between the OS and the WTs is a two-dimensional Euclidean distance. In reality, this length is always modified by a redundancy factor.
All WTs are assumed to be operated at 1 p. u. voltage and the power factor is assumed to be 0.95.
Based on the above assumptions, this article will establish an optimization model of CCLOP based on OWFD.

2.2. Optimization Model

2.2.1. Objective Function

In graph theory, a spanning tree is a subgraph of an undirected graph G [7]. The subgraph contains all vertices and some branches in G. All vertices have paths to communicate with each other, but there are no loops. The spanning tree of a weighted G can be expressed as:
G T = ( V ,   B T ,   W T ) ,   G T G ,   B T B ,   W T W
where V represents the set of all vertices, B is the set of all branches connecting the vertices in V. W is the set of weights corresponding to all branches in B. G T is a subgraph of G, representing a spanning tree of G. Similarly, B T is the set of the branches in B while W T is the weight corresponding to all branches in B T .
According to the definition of the spanning tree shown above, we know that if the OS (represented by V O S ) and WTs (represented by V W T ) are regarded as the vertices of G, the cost of the cable connecting the two vertices (between OS and WT or WT and WT) is the weight of the branch. Then the CCLOP problem can be transformed into a mathematical problem of finding a G T with minimum total weight in a weighted G. The objective function of the CCLOP problem includes two parts: CAPEX of cables and CELC, of which the detailed mathematical expressions are discussed in the following sections.
The cost of CAPEX of branch m with cable type p is calculated as:
Z CAPEX , G T , m , p = C p L G T , m φ G T , m , p
Then, the total cost of CAPEX of all cables can be formulated as:
Z CAPEX , G T = m = 1 M p = 1 P Z CAPEX , G T , m , p
In addition, the power losses along the cables reduce the profits of the wind farm owner and appear all over the lifetime of the offshore wind farm. In order to account for it, the net present value of the CELC during the wind farm’s lifetime is incorporated into the objective function as follows:
Z CELC , G T , m , p = ( I G T , m , p ) 2 R p L G T , m φ G T , m , p
where I G T , m , p is the sum of currents flowing through the branch m. In a homogeneous wind farm, the power through the cable is usually calculated according to the total number of WTs that the cable undertaken which has been adopted in [7,12,20,24]. In an OWFD, the rated powers of WTs are not identical. In order to calculate the power flow in each cable when all the WTs are in full load condition, the WTs’ types associated with the cable should be identified. A simple example is shown to explain how the I G T , m , p is calculated in such a case.
Supposing that there are two offshore wind farms with the same layout (one OS and four WTs are with the same coordinates). One is the homogeneous wind farm which composed of the same type of WTs while the other one is the OWFD. All the WTs in each wind farm are operated in full load condition. The rated power for type one and two WT is 2 MW and 3 MW, respectively. Neglecting the power losses, the actual power along each branch are compared in Figure 1.
In Figure 1a, there are four WTs transmitting their energies to cable (1,2), which has the same type. Therefore, the actual carrying capacity of cable (1,2) in Figure 1a is 4 × 2 = 8 MW. In Figure 1b, the WTs transmitting energies to cables (1,2) include two type 1 WTs and two type 2 WTs. Therefore, the actual power reached at the cable (1,2) in Figure 1b is 2 + 3 + 3 + 2 = 10 MW. It can be seen that when there is only one type of WTs in the wind farm, it can be obtained by the product of the number of WTs transmitting energies to the cable and the rated WT’s power. However, if the WTs in the wind farm are not identical, the corresponding WT types and number of each type of WTs related to the cable should be identified so that the actual power of each cable can be calculated. As a result, Equation (5) is modified to resolve the above issue:
I G T , m , p = h = 1 H G T , m I rate , q δ G T , q , h
where H G T , m is the number of the WTs that transmit their power through the branch m of G T while I rate , q is decided by the WT type q, the mathematical expression is written as:
I rate , q = P rate , q 3 c o s ψ U rate , q
Then the CELC of all cables can be written as:
Z CELC , G T = 3 × ξ = 1 T ( 1 + η ) ξ ψ Δ C loss ( m = 1 M p = 1 P Z CELC , G T , m , p )
Finally, the overall objective function is defined as:
Z min , total = min ( Z CAPEX , G T + Z CELC , G T )

2.2.2. Constraints

The main constraints involve in this work are summarized in the following:
G T G
I G T , m , p I max , p ,     m = 1 , 2 , , M   ,     p = 1 , 2 , , P
1 F G T , V OS , in 10 ,   V OS V
1 F G T , V WT , m , in 2   ,   V WT ,   m V WT   ,     m = 1 , 2 , , M
F G T , V WT , m , out = 1 ,   V WT ,   m V WT   ,     m = 1 , 2 , , M
  B m 1 B m 2 = ϕ ,     B m 1 B T ,   B m 2 B T ,   m 1 m 2 [ 1 , M ]
Equation (9) is to ensure that G T is a spanning tree of G. Based on the principle of equipment protection, the actual current flowing through each cable cannot exceed its maximum current carrying capacity which is constrained in Equation (10). From the engineering point of view, some constraints are made from Equations (11)–(13) which represents the maximum number of cables that an offshore substation can hold, the number of incoming cables for each WT, and the number of outgoing cables for one WT respectively. The crossed cables would not only endanger the operation of OWF but also increase the installation cost. To ensure an uncrossed cable connection layout, Equation (14) is added to this work.

3. Proposed Optimization Method

In this section, a novel V-APSO-LS algorithm is proposed for solving the above intractable problem. The canonical APSO algorithm structure is specified at first. Secondly, a simple example is elaborated to show how the k-ring shaped Voronoi neighbor-based model is created. Then, the local search enhanced APSO is shown. The overall framework of the V-APSO-LS is specified in the end.

3.1. APSO Algorithm

As one of the population-based methods, particle swarm optimization was firstly proposed by Kennedy and Eberhart for solving a continuous nonlinear problem [25,26]. It is a random search algorithm based on group cooperation, which is developed by simulating the foraging behavior of birds. In this algorithm, particles have two properties: speed and position. The speed of the particles represents the speed of movement, and the position of the particles represents the current solution. The algorithm finds the optimal solution by updating the speed and position of the particles, which can be expressed as:
v i , t + 1 = ω · v i , t + r 1 · rand ( ) · ( LB i , t x i , t ) + r 2 · rand ( ) · ( GB i , t x i , t )
x i , t + 1 = x i , t + v i , t + 1
where the first term in Equation (15) is a weighting factor, the second term is a cognitive learning part, and the third term is a social learning part. The later research [27] found that the larger the w is, the stronger the global optimization searching ability that the algorithm has, and thus shows a weaker the local optimization searching ability; if it gets smaller, the algorithm’s local optimization ability is stronger, and the global optimization ability is weaker. It can be seen that the global and local optimization searching performance of the algorithm can be adjusted by tuning the parameters. In the light of this, an APSO with a non-linear reduction of inertia weight was proposed in [27]. The weight was changed according to the following rule:
ω = 1 1 + 1 . 5 e 2 . 6 μ [ 0 . 4 , 0 . 9 ] ,     μ [ 0 , 1 ]
μ = d avg d min d max d min
According to our previous work [7,20,28], APSO with initial values of r 1   =   2 , r 2   =   2 and ω   =   0 . 9 has shown good performance in solving CCLOP compare with MST. Furthermore, as indicated in [29], particle swarm optimization has higher computation efficiency in solving non-linear problem compared with other meta-heuristic algorithms. Therefore, the APSO is selected as the core algorithm to optimize CCLOP in this work.

3.2. V-APSO-LS Algorithm

In this section, a new coding and decoding strategy based on k-ring shaped Voronoi neighbor and APSO, named V-APSO-LS, is proposed. In order to further enhance the local search ability of the algorithm, a random local search strategy is incorporated and specified at last.

3.2.1. K-Ring Shaped Voronoi Neighbor

A Voronoi diagram is a typical metric space decomposition method [22,30,31,32]. According to a nearest-neighbor rule, each discrete point in the Voronoi diagram is associated with the region of the plane close to it. Each discrete point corresponds to only one region. Due to the bisected features of the Voronoi diagram in spatial decomposition, it is often applied to solve the problems such as nearest points, adjacencies, and proximity. In 1991 Aurenhammer et al. proposed the definitions of the Voronoi diagram, Voronoi region, and Voronoi branches [30,31,32]. In 2017, Fang proposed the definition of k-ring shaped Voronoi neighbor [22]. In the following, the detailed definitions related to the Voronoi diagram are specified.
Definition 1. 
The discrete points distributed on the Voronoi diagram are defined as the Voronoi site. Hence, the WTs and OS are Voronoi sites in this work.
Definition 2. 
The corresponding Voronoi regions of the site E 1 and E 2 are VR ( E 1 ) and VR ( E 2 ) respectively. If there is a Voronoi edge between VR ( E 1 ) and VR ( E 2 ) , then E 1 and E 2 are Voronoi neighbors. For instance, there are totally 29 Voronoi sites inFigure 2 while E 1 and E 2 are one pair of the Voronoi neighbors.
Definition 3. 
The Voronoi distance between E 1 and E 2 is defined as the number of Voronoi branch intersected by the line connecting E 1 and E 2 . For example, the Voronoi distance between E 1 and E 2 in Figure 2 is one.
Definition 4. 
The k-ring shaped Voronoi neighbors set of E 1 is referring to the set of k points with Voronoi distance from E 1 , denoted as KN ( E 1 ,   k ) . For example, KN ( E 1 ,   1 ) , KSVN ( E 1 ,   2 ) and KN ( E 1 ,   3 ) is the 1-ring, 2-ring, 3-ring Voronoi neighbors sets of E 1 which are shown in Figure 2, respectively.
As can be seen from Figure 2, the metric space is divided by the Voronoi diagram. Instead of judging two points’ proximity relationship by the Euclidean distance, the Voronoi distance is obtained based on redefining the proximity relationship between the sites, which is a more scientific method of judging the proximity relationship [23]. The reason is that a shorter Euclidean distance two sites does not also mean that their proximity relationship is closer. A simple example is given in this work to illustrate the abovementioned situation, as shown in Figure 3.
In Figure 3, there is a wind farm with 1 OS and 5 WTs (totally 2 types of WTs). It can be seen that the Euclidean distance between the 6th node to the OS is closer than the distance between the 3rd node to the OS. However, the 3rd node is the 1st order neighbor of OS while the 6th node is the 2nd order neighbor of the 1st node. Hence, the 3rd node is a closer neighbor of OS compared with the 6th node. Similarly, the 4th node and 5th node are also the closer neighbors of OS compared with the 6th node. It can be seen that in generating the cable connection layout, when the WT connected to the OS needs to be selected, the priority of the 3rd node should be higher than the 6th node, and the 6th node should be connected to the 4th site or 5th node.
It can be seen that the Voronoi diagram judge the adjacent relationship between OS and WT, and between WT and WT in a wind farm more scientifically, which is conducive to selecting a better branch for cable connection in the layout process. The next section will introduce the KSVN based codec scheme, and elaborate on how to apply KSVN to cable connection layout optimization in this article.

3.2.2. Encoding and Decoding of Particle

Encoding of Particle
According to the model shown in Section 2, the algorithm in this paper does not only need to get the basic cable connection information, but also needs to consider the further influence of energy loss on the selection of cable type on the layout. Therefore, the particle code must contain two parts of information, one is the cable connection layout information, the other is the information of each cable type. Based on the above reasons, assuming the number of WTs is m, the particle i coding in this paper is shown in Figure 4.
In Figure 4, S 1 S M represent the connection index which belongs to [1,m × (Mm + 1)] and m = 1 , 2 , , M while S M + 1 S 2 M are the Cable Type Incremental Number (CTIN) index which are in the range of [1, P].
Decoding of Particle
Particle decoding includes two stages, the first stage is to get the cable connection information, and the second stage is to select the cable model according to the current cable connection information. In the whole decoding process, V-APSO-LS needs to use 6 basic sets and 3 matrices as follows:
Set A:
Including the vertices which have already been connected in MST.
Set B:
Including the vertices which have not yet been connected in MST.
Set C:
Including the length of each branch connecting vertices in set A.
Set D:
Including actual current of each branch in MST.
Set E:
Including the minimal cable type of each branch in MST.
Set F:
Including the final cable type of each branch in MST.
Adjacency matrix AM: Including all the Voronoi distance and the geometrical distance between each pair of adjacent stochastic vertices which are from set A and set B. Where set A is the index of the column and set B is the index of the row.
Sorted-adjacency matrix SAM: All the branches in the adjacency matrix AM are sorted from small to large according to the Voronoi distance. If the distance of two sides is equal, the branches are sorted from small to large according to the length of the branch. Finally, the sorted-adjacency matrix SAM is obtained.
Cable matrix: Including the current carrying capacity, resistance and price per unit length for each type of cable, the cable data is arranged according to the ampacity of the cable from small to large.
In the first stage, set A only contains OS, set B contains all WTs, while set C, D, E and F are empty. First of all, according to the vertices of set A and B, the adjacency matrix is generated, and the branches in the adjacency matrix are sorted to get the sorted adjacency matrix. Then, V-APSO-LS selects a new branch according to the connection index of particle i and sorted adjacency matrix, and puts the branch related information into set C. At the same time, the vertices corresponding to the branch are deleted from set B and then put into set A. This phase will not end until set B is empty. It is worth noting that in the selection process of vertex and its corresponding branch, if the cable connection layout after the access of the branch does not meet the constraint from Equation (10) to Equation (14), the next branch in the sorted adjacency matrix is selected in order until the constraint is met. It should be noted that the new branch cannot cross the other branches in the spanning tree in each step which is just by the method specified in [2].
In the second stage, the actual power through each branch (i.e., cable) can be obtained according to the information in set C. The current passing through each branch (i.e., cable) will be calculated by Equations (5) and (6), and the current value is put into set D accordingly. Further, according to the current value in set D, the cable with the minimum current carrying capacity while meeting the constraint of Equation (10) will be selected [7], and put into set E. Finally, the selected cable type of each branch is decided according to the sum of the corresponding index in set E and the CTIN index of particle i which is the information that put into the set F. It is worth noting that if the added value exceeds the maximum index of the optional cable type, then the cable type will be set to the maximum cable type index.
To sum up, after two stages of decoding, the obtained information in set C and F can constitute a cable connection layout. Combined with set D, the corresponding cost of the cable connection layout can be calculated. In order to better explain the two stages of particle coding, an example is given in the following. It is assumed that there are totally six WTs with two types in the reference wind farm, the position of particle i is assumed to be [2,4,1,1,2,0,1,2,1,0]. There are three types of optional cable models (the larger the cable index, the greater the carrying capacity of the cable). Among them, the first five positions of the particle represent the connection index, and the last five positions represent the CTIN index. The detailed decoding process of particle i is shown in Figure 5 and Figure 6 which represent the first and the second stages, respectively.
Assuming the position of particle i has been decoded the 1st dimension, and the current set A is [1,2], set B is [3,4,5,6], and set C is {{(1, 2), 1.00}}. Then, decoding the 2nd dimension of the particle i is shown in Figure 5. First, according to set A and set B, we get the adjacency matrix AM. Second, sorting all the branches in the AM from small to large according to the Voronoi distance, and then we get the sorted-adjacency matrix SAM. Third, because the value of the second dimension of particle i is four, the 4th branch (1,3) of the SAM is pointed. Lastly, assuming that the branch (1,3) satisfies the constraints from Equation (10) to Equation (14), then the relevant information of this branch will be put into set C. So the information in set C will be updated to {{(1,2),1.00}, {(1,3),2.00}}. Meanwhile, set A and set B are updated to [1,2,3] and [4,5,6] respectively. Similarly, the 3rd to 5th dimensions of the particle i can also be decoded according to the above steps and then set A will be created as [1,2,3,4,6,5] and set C as {{(1,2),1.00}, {(1,3),2.00}, {(1,4),0.50}, {(4,6),0.71}, {(4,5),1.12}}.
Next, the second stage of decoding is executed, as shown in Figure 6. First, according to the cable connection information of set C, set D that includes the current passing through each branch is calculated by Equations (5) and (6). Second, assuming set D is obtained as [1.0,1.0,3.2,1.2,1.0], under the rule of selecting the minimum current-carrying cable type, set E [1,1,2,1,1] is calculated. Third, the last five dimensions of the position of particle i are [0,1,2,1,0], as introduced above it should be added with the information in set E, so that the set F is initialized as [1,2,4,2,1]. However, since there is no cable with cable type 4, cable type 3 will be picked up and then set F becomes [1,2,3,2,1]. Finally, set C, set F and set D can get the cable connection layout corresponding to the particle i’s position and its related information.

3.2.3. Local Search Strategy

LSS is often used to improve the local optimization ability of the intelligent algorithm [33,34]. LSS directly optimizes the global optimal particles, and the updated global optimal particles can enhance the social learning ability of particles, as shown in Equation (19):
v i , t + 1 = ω · v i , t + r 1 · rand ( ) · ( LB i , t x i , t ) + r 2 · rand ( ) · ( GB local , i , t x i , t )
where GB local , i , t refers to the best position of all particles after using LSS at iteration t. It is not hard to see that Equation (19) is the modification of Equation (15). According to the characteristics of OWFD and APSO, a new LSS is proposed in this paper. As far as we know, previous scholars mainly applied LSS for OWF micro-siting [35,36,37], then our work is the first to present a LSS to solve the CCLOP.
The mechanism of the proposed LSS can be explained as: randomly change the value of one dimension of the global optimal particle position. If the cable connection layout corresponding to the new particle benefits the objective function, it will replace the original global optimal particle, otherwise, it will be re-executed until the number of local search iteration reaches the predefined limit. Since the particle position includes the connection indexes and the CTIN indexes, the proposed LSS can help optimize the cable connection and cable selection at the same time. Continue with the same example as shown in Figure 6, two cases of particle updating using LSS are shown in Figure 7.
As can be seen in Figure 7, there are in total two cases. In the first case, it is assumed that the LSS randomly selects the 5th connection index and changes its value to 1, and then the particle 2 positon [2,4,1,1,1,0,1,2,1,0] is obtained. Further, decoding the particle 2 positon, the cable connection layout of particle 2 is shown in the bottom left of Figure 7. It can be seen from the figure that the cable connection layouts of particle 1 and particle 2, after triggering the LSS, the branch (4,5) becomes branch (1,5). In the second case, supposing that the LSS randomly selects the 1st CTIN index, and changes the value on this position to 1 to get the new particle 3 position [2,4,1,1,2,1,1,2,1,0]. Similarly, decoding the particle 3 positon, the cable connection layout of particle 3 is shown in the bottom right of Figure 7. It can be seen from the cable connection layout after two particle codec that after LLS, the cable type of branch (1,2) changes from type 1 to type 2.
It is worth noting that the above two cases only generate candidate particles. The information of the global optimal particles will never update except that the LSS helps generate a better solution. Therefore, this strategy can not only retain the original global search ability of APSO, but also improve the local searching ability of the algorithm. The performance of LSS in finding a better solution will be further demonstrated in the simulation.

3.3. Optimization Framework

In this paper, we use the V-APSO-LS algorithm to solve the CCLOP for an OWFD. The optimization framework is shown in Figure 8. As can be seen in Figure 8, the algorithm starts with initializing Voronoi distance matrix which is generated as described in Section 3.2.1. According to the locations of the WTs and OS, the particle population and its related parameters are initialized in the next step. After importing such information into fitness function, the fitness value of each particle is obtained followed by the procedure in Section 3.2.2, which is the basis of the later fitness comparison. Further, the position and velocity of particles and their corresponding fitness values are updated, and the current global optimal solution will be updated via the LSS in Section 3.2.3. The above steps will not be terminated until the stop criteria are met. Finally, the optimal cable layout and the corresponding costs will be obtained.

4. Case Study

In this section, a reference OWFD is introduced at first. Then, two scenarios are presented with detailed results comparison and discussion at the end. In addition, the simulation software is MatLab, and the experimental simulation environment is Windows 10 OS, Intel Corei7-6800K [email protected] GHz, 16 GB RAM.

4.1. Reference Wind Farm

In this work, an offshore wind farm with one OS and 72 WTs was selected. There are three types of WTs in this wind farm, with the rated power of 3 MW, 3.3 MW and 6.4 MW respectively. The positions of the OS and WTs are already known. The technical and economic information related to the wind farm is shown in Table A1 while the cable specifications are shown in Table A2. In addition, all cables are 3-core cross-linked polyethylene AC cables.

4.2. Simulation Results and Discussion

In this section, two scenarios are presented to validate the proposed method. According to our many trials, when the number of iterations is set to 150 and the population size is set to 500, the algorithm can solve CCLOP of OWFD better. Then, to ensure the comparison fairness, the number of iterations of each algorithm is set to 150 and the population size of each algorithm is set to 500, while each algorithm is run 10 times independently. The best solution among the 10 trials is selected as the optimal cable connection scheme, the result of which will be discussed and analyzed at last.

4.2.1. Scenario I

In this scenario, two optimization models for solving CCLOP are adopted and solved by the proposed V-APSO-LS. The first model (Model 1) is a cable connection layout scheme constructed on the minimum cross-sectional area cable selection rule, that is, cables with the small cross-sectional areas but satisfying Equation (11) are selected as far as possible for layout which is the model proposed in [7]. The second model (Model 2) further considers the impact of the CELC on the cable type selection, which is the model proposed in Section 2. The simulation results are shown in Figure 9.
It can be seen in Figure 9, no crossed cables appear in either the layout obtained by Model 1 or Model 2. Six kinds of cables are used for Model 1 case while in Model 2 seven kinds of cables are used. Compared with Figure 9a, there are fewer blue cables in Figure 9b which shows that Model 2 chooses more cables with the larger cross-sectional area. The reason is that Model 2 considering the influence of CELC. To highlight the difference, the branches with the larger cross-sectional area are shown in Figure 10.
The statistic results for running each model using V-APSO-LS are shown in Figure 11. As can be seen from Figure 11, the optimal solutions of Model 1 and Model 2 appear in the 8th and 2nd experiments respectively. Compared to Model 1, Model 2 gives a better cable connection layout. In order to better understand the advantages and disadvantages of the two models, the detailed cost information of the optimal cable connection layout obtained by each model is shown in Table 1.
It can be seen from Table 1 that the cable procurement cost of Model 2 is 3030.80 kEUR more expensive than that of Model 1. However, the energy consumption cost of Model 2 is only 10,272.63 kEUR, which is 38.36% less than that of Model 1. As a result, the total cost of Model 2 is 9.43% lower than that of Model 1. Although more expensive cables are purchased in Model 2, the savings in CELC lead to a lower overall cost compared with the situation of using Model 1.

4.2.2. Scenario II

In order to verify the proposed V-APSO-LS, and the effect of using KSVN and LSS, the other three algorithms are used to solve Model 2 and compared with V-APSO-LS. The first algorithm, Voronoi diagram-based adaptive particle swarm optimization (V-APSO) is constructed by V-APSO-LS without using LSS. The second adaptive particle swarm optimization minimum spanning tree (APSO-MST) algorithm is from [7]. The third algorithm, Voronoi diagram-based wild goats algorithm with local search (V-WGA-LS), is to replace the APSO in V-APSO-LS with wild goats algorithm proposed in [38], but still, retain KSVN and LSS. The 10 trials results of the four algorithms are shown in Figure 12, and the average computational times of V-APSO-LS, V-APSO, APSO-MST and V-WGA-LS are 1284 s, 1176 s, 1522 s, 2452 s, respectively.
As can be seen from Figure 12, the optimal solutions of V-APSO-LS, V-APSO, APSO-MST and V-WGA-LS appeared at the 2nd, 3rd, 6th and 5th times respectively. The results obtained by the APSO-MST are the worst. In contrast, the V-APSO finds the best solution. The optimized layout obtained by V-APSO, APSO-MST and V-WGA-LS are shown in Figure 13. Then, the fitness value corresponding to each iteration is shown in Figure 14.
As can be seen from Figure 14, the result obtained by the V-APSO with KSVN is better than the result obtained by the APSO-MST which shows the rationality of using the Voronoi distance to judge the proximity relations between points in the layout. It is worthy of note that after the 10th iteration, the convergence rate of V-APSO becomes slower. Furthermore, V-APSO-LS with LSS is much better than V-APSO after 10 iterations. It is verified that the LSS proposed in this paper can greatly enhance the searching ability of the APSO.
In addition, after the 10th iteration, the trace of the convergence of V-WGA-LS is similar to that of V-APSO-LS. However, the convergence speed of V-APSO-LS is faster while retaining a better solution which shows the superiority of the APSO than the wild goats algorithm. The detailed cost comparison of the optimal cable connection layout obtained by each algorithm is shown in Table 2.
As can be seen in Table 2, compared with the best, mean and worst solutions of the algorithms respectively, the CAPEXes obtained by the APSO-MST are always the lowest, whereas due to the higher CELCs, the total costs are the largest. Therefore, the accuracy and stability of the V-APSO-LS are better than the other three algorithms. Furthermore, the V-APSO-LS achieves the best solution which is 8.04%, 12.74% and 5.35% lower than the solutions found by the V-APSO, APSO-MST, and V-WGA-LS respectively.

5. Conclusions

In this paper, the influence of energy loss on the overall cable connection structure economy is fully considered. Aiming at minimizing the total cost of OWFD, a CCLOP model based on OWFD is established. In order to get an optimized result, this paper proposes a V-APSO-LS which takes APSO as the core and adopts the KSVN and LLS to enhance the optimized performance. In the simulation experiment of Scenario I, this paper proves that fully considering the impact of CELC on the overall cable connection layout is conducive to reducing the total cost of OWFD during the operation period. Furthermore, in scenario II, the layout scheme of V-APSO-LS is superior to V-APSO, APSO-MST and V-WGA-LS, which proves the effectiveness of the algorithm and the optimization strategy proposed in this paper. Therefore, this paper provides effective theoretical guidance and technical support for the cable connection layout optimization of OWFD.
In the future, CCLOP of large-scale OWFD can be established and optimized considering the number and location of offshore substations, and then compared with the used industrial layout. At the same time, more optimization strategies is promising to be adopted for improve the ability of the intelligent algorithm in solving the complex combinational optimization tasks in terms of offshore wind energy.

Author Contributions

Writing—original draft preparation, Y.Q.; writing—review and editing, P.H.; review and editing, G.L., Z.Y., G.Y. and Z.D.; writing—original draft preparation, R.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Key Project in Higher Education of Guangdong Province, China under grant No. 2020ZDZX3030, Young Innovation Talents Project in Higher Education of Guangdong Province, China under under No. 2018KQNCX333, National Science Foundation of China under grants 52077213 and 62003332, Natural Science Foundation of Guangdong (No. 2018A030313755).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

This research was funded by Key Project in Higher Education of Guangdong Province, China under grant No. 2020ZDZX3030, Young Innovation Talents Project in Higher Education of Guangdong Province, China under under No. 2018KQNCX333, National Science Foundation of China under grants 52077213 and 62003332, Natural Science Foundation of Guangdong (No. 2018A030313755).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

OWFOffshore Wind Farm H G T , m Number of the WTs that transmit their power through the branch m of G T
OWFDOffshore Wind Farm with Diverse Wind Turbine TypeshIndex of the WTs that transmit their power through one branch
OSOffshore Substation φ G T , m , p If branch m of G T uses type p cable, φ G T , m , p = 1 , otherwise φ G T , m , p = 0 .
WTWind Turbine δ G T , q , h If the hth WT transmits its power through the branch m of G T is WT type q, δ G T , q , h = 1 , otherwise δ G T , q , h = 0 .
CAPEXCapital ExpenditureTLifetime of the wind farm [Year]
CELCCable Energy Losses Cost C loss Unit cost of energy loss of cable [EUR/MWh]
CCLOPCable Connection Layout Optimization Problem ψ Δ Duration time of peak energy loss [h]
MSTMinimum Spanning Tree ψ Power factor
APSOAdaptive Particle Swarm Optimization η Interest rate
APSO-MSTAdaptive Particle Swarm Optimization Minimum Spanning Tree Z CAPEX , G T , m , p Cost of CAPEX of branch m with cable type p in G T [EUR]
KSVNK-ring Shaped Voronoi Neighbor Z CAPEX , G T Total cost of CAPEX in G T [EUR]
LSSLocal Search Strategy Z CELC , G T , m , p Cable energy losses of branch m with cable type p in G T [MWh]
CTINCable Type Incremental Number Z CELC , G T Cost of CELC in G T [EUR]
V-APSO-LSVoronoi Diagram based Adaptive Particle Swarm Optimization with Local Search Z min , total Total cost for the multiple wind turbine types offshore wind farm [EUR]
V-APSOVoronoi Diagram based Adaptive Particle Swarm Optimization V WT ,   m The vertex m which notes WT m in V WT
V-WGA-LSVoronoi Diagram based Wild Goats Algorithm with Local Search F G T , V OS , in Number of Feeders in G T
Gan undirected graph F G T , V WT , m , in Number of input cables of WT m in G T
G T Sub-graph belongs to G F G T , V WT , m , out Number of output cables of WT m in G T
VVertex set in G B m 1 , B m 2 Branch m1 and branch m2 in B T
BAll branches that connect V in GiIndex of particles
B T All branches that connect V in G T T ite Number of iterations
WWeight of each branch in GtIndex of iterations
W T Weight of each branch in G T x i , t , x i , t + 1 Position of particle i at iteration t and t + 1
V O S The vertex which notes the OS in V v i , t , v i , t + 1 Velocity of particle i at iteration t and t + 1,
V W T The vertex set which notes WTs in V L B i , t Best position of particle i at iteration t
MNumber of wind turbines G B i , t Best position of all particles at iteration t
PNumber of cable types ω Inertia weight
QNumber of WT types r 1 , r 2 Learning factors
MIndex of a branch in G T ; Index of WT in V W T ; Index of the dimension of the particle position μ Evolutionary factor
RIndex of wind farm operation period d a v g Mean distance from the global best particle to the others
PIndex of cable type d m i n , d m a x Minimal and maximum mean distances from one particle to the others
QIndex of WT type E 1 , E 2 Site E 1 and site E 2 in the Voronoi diagram
C p Unit cost of cable type p (EUR/m) VR ( E 1 ) ,   VR ( E 2 ) Voronoi regions of site E 1 and site E 2
R p Unit resistance of cable type p (𝛀/km) KN ( E 1 ,   k ) The k-ring shaped Voronoi neighbors set of site E 1
I max , p Current-carrying capacity of cable type p (A) GB local , i , t Best position of all particles after using LSS at iteration t
P rate , q Rated active power of WT type q (MW) T LS Maximum attempts of local search
U rate , q Rated voltage of WT type q (kV)tlsThe tlsth attempt of T LS
I rate , q Rated current of WT type q (A)AMAdjacency matrix
L G T , m Length of the branch m in G T (m)SAMSorted-adjacency matrix
I G T , m , p Actual current in branch m with cable type p in G T (A)

Appendix A

Table A1. Other basic technical and economic data of the wind farm.
Table A1. Other basic technical and economic data of the wind farm.
ItemValueItemValue
U rate , q 66.0 kV η 0.049
T25 Year ψ Δ 2608 h
cos   ψ 0.95 C loss 109.055 EUR/MWh
Table A2. Basic technical and economic data for 66 KV, XRUHAKXS cable.
Table A2. Basic technical and economic data for 66 KV, XRUHAKXS cable.
Type No.Sectional AreaPrice (EUR/m) Resistance   ( Ω / k m ) Ampacity (A)
T195233.6340.2500260
T2120251.3400.1458315
T3150271.3550.1250350
T4185301.6330.0970428
T5240343.8440.0729440
T6300382.4620.0600520
T7400443.9180.0437600

References

  1. Chen, Y.; Dong, Z.; Meng, K.; Luo, F.; Yao, W.; Qiu, J. A novel technique for the optimal design of offshore wind farm electrical layout. J. Mod. Power Syst. Clean Energy 2013, 1, 258–263. [Google Scholar] [CrossRef] [Green Version]
  2. 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]
  3. Hou, P.; Hu, W.; Soltani, M.; Chen, C.; Chen, Z. Combined optimization for offshore wind turbine micro siting. Appl. Energy 2017, 189, 271–282. [Google Scholar] [CrossRef]
  4. Hou, P.; Hu, W.; Soltani, M.; Chen, Z. Optimized placement of wind turbines in large-scale offshore wind farm using particle swarm optimization algorithm. IEEE Trans. Sustain. Energy 2015, 6, 1272–1282. [Google Scholar] [CrossRef] [Green Version]
  5. Ling-ling, H.; Ning, C.; Hongyue, Z.; Yang, F. Optimization of large-scale offshore wind farm electrical collection systems based on improved FCM. In Proceedings of the International Conference on Sustainable Power Generation and Supply (SUPERGEN 2012), Hangzhou, China, 8–9 September 2012; pp. 1–6. [Google Scholar]
  6. Huang, L.; Yang, F.; Guo, X. 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, Supergen, Nanjing, China, 6–7 April 2009; pp. 1–4. [Google Scholar]
  7. Hu, W.; Chen, Z.; Hou, P. Optimisation for offshore wind farm cable connection layout using adaptive particle swarm optimisation minimum spanning tree method. IET Renew. Power Gener. 2016, 10, 694–702. [Google Scholar] [CrossRef]
  8. Chen, Y.; Dong, Z.Y.; Meng, K.; Luo, F.; Xu, Z.; Wong, K.P. Collector system layout optimization framework for large-scale offshore wind farms. IEEE Trans. Sustain. Energy 2016, 7, 1398–1407. [Google Scholar] [CrossRef]
  9. Li, D.D.; He, C.; Fu, Y. Optimization of internal electric connection system of large offshore wind farm with hybrid genetic and immune algorithm. In Proceedings of the International Conference on Electric Utility Deregulation and Restructuring and Power Technologies, Nanjing, China, 6–9 April 2008; pp. 2476–2481. [Google Scholar]
  10. Gonzalez-Longatt, F.M.; Wall, P.; Regulski, P.; Terzija, V. Optimal electric network design for a large offshore wind farm based on a modified genetic algorithm approach. IEEE Syst. J. 2012, 6, 164–172. [Google Scholar] [CrossRef]
  11. Lumbreras, S.; Ramos, A.; Sanchez-Martin, P. Offshore wind farm electrical design using a hybrid of ordinal optimization and mixed-integer programming. Wind Energy 2015, 18, 2241–2258. [Google Scholar] [CrossRef]
  12. 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]
  13. Chowdhury, S.; Zhang, J.; Messac, A.; Castillo, L. Optimizing the arrangement and the selection of turbines for wind farms subject to varying wind conditions. Renew. Energy 2013, 52, 273–282. [Google Scholar] [CrossRef]
  14. Abdulrahman, M.; Wood, D. Investigating the Power-COE trade-off for wind farm layout optimization considering commercial turbine selection and hub height variation. Renew. Energy 2017, 102, 267–278. [Google Scholar] [CrossRef]
  15. Chowdhury, S.; Zhang, J.; Messac, A.; Castillo, L. Unrestricted wind farm layout optimization (UWFLO): Investigating key factors influencing the maximum power generation. Renew. Energy 2012, 38, 16–30. [Google Scholar] [CrossRef]
  16. Feng, J.; Shen, W.Z. Design optimization of offshore wind farms with multiple types of wind turbines. Appl. Energy 2017, 205, 1283–1297. [Google Scholar] [CrossRef]
  17. Chen, K.; Song, M.X.; Zhang, X.; Wang, S.F. Wind turbine layout optimization with multiple hub height wind turbines using greedy algorithm. Renew. Energy 2016, 96, 676–686. [Google Scholar] [CrossRef]
  18. Gong, X.; Kuenzel, S.; Pal, B.C. Optimal wind farm cabling. IEEE Trans. Sustain. Energy 2018, 9, 1126–1136. [Google Scholar] [CrossRef]
  19. Dahmani, O.; Bourguet, S.; Machmoum, M.; Guerin, P.; Rhein, P.; Josse, L. Optimization of the connection topology of an offshore wind farm network. IEEE Syst. J. 2015, 9, 1519–1528. [Google Scholar] [CrossRef]
  20. Jin, R.; Hou, P.; Yang, G.; Qi, Y.; Chen, C.; Chen, Z. Cable routing optimization for offshore wind power plants via wind scenarios considering power loss cost model. Appl. Energy 2019, 254, 113719. [Google Scholar] [CrossRef]
  21. Duczmal, L.H.; Moreira, G.J.P.; Burgarelli, D.; Takahashi, R.H.C.; Magalhaes, F.C.O.; Bodevan, E.C. Voronoi distance based prospective space-time scans for point data sets: A dengue fever cluster analysis in a southeast Brazilian town. Int. J. Health Geogr. 2011, 10, 29. [Google Scholar] [CrossRef] [Green Version]
  22. Fang, Z.; Tu, W.; Li, Q.; Shaw, S.-L.; Chen, S.; Chen, B.Y. A Voronoi neighborhood-based search heuristic for distance/capacity constrained very large vehicle routing problems. Int. J. Geogr. Inf. Sci. 2013, 27, 741–764. [Google Scholar] [CrossRef]
  23. Fang, Z.; Li, Q.; Shaw, S.-L.; Chen, B.; Chen, B. A bi-level Voronoi diagram-based metaheuristic for a large-scale multi-depot vehicle routing problem. Transp. Res. Part Logist. Transp. Rev. 2014, 61, 84–97. [Google Scholar] [CrossRef]
  24. Fischetti, M.; Pisinger, D. Optimizing wind farm cable routing considering power losses. Eur. J. Oper. Res. 2018, 270, 917–930. [Google Scholar] [CrossRef] [Green Version]
  25. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the ICNN’95—International Conference on Neural Networks, Perth, WA, Australia, 27 November–1 December 1995; Volume 4, pp. 1942–1948. [Google Scholar]
  26. Kennedy, J. The particle swarm: Social adaptation of knowledge. In Proceedings of the 1997 IEEE International Conference on Evolutionary Computation (ICEC ’97), Indianapolis, IN, USA, 13–16 April 1997; pp. 303–308. [Google Scholar]
  27. Zhan, Z.-H.; Zhang, J.; Li, Y.; Chung, H.S.-H. Adaptive particle swarm optimization. IEEE Trans. Syst. Man Cybern. Part (Cybern.) 2009, 39, 1362–1381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Hou, P.; Hu, W.; Zhang, B.; Soltani, M.; Chen, C.; Chen, Z. Optimised power dispatch strategy for offshore wind farms. IET Renew. Power Gener. 2016, 10, 399–409. [Google Scholar] [CrossRef]
  29. Rania, H.; Babak, C.; Olivier, W. A comparison of particle swarm optimization and the genetic algorithm. In Proceedings of the 46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, Austin, TX, USA, 18–21 April 2005; pp. 1–13. [Google Scholar]
  30. Aurenhammer, F. Voronoi diagrams—A survey of a fundamental geometric data structure. ACM Comput. Surv. 1991, 23, 345–405. [Google Scholar] [CrossRef]
  31. Stoyan, D. Spatial tessellations: Concepts and applications of voronoi diagrams. Biom. J. 2010, 36, 146. [Google Scholar] [CrossRef]
  32. Chen, J.; Li, C.; Li, Z.; Gold, C. A Voronoi-based 9-intersection model for spatial relations. Int. J. Geogr. Inf. Sci. 2001, 15, 201–220. [Google Scholar] [CrossRef]
  33. Fischetti, M.; Monaci, M. Proximity search heuristics for wind farm optimal layout. J. Heuristics 2016, 22, 459–474. [Google Scholar] [CrossRef] [Green Version]
  34. Wagner, M.; Day, J.; Neumann, F. A fast and effective local search algorithm for optimizing the placement of wind turbines. Renew. Energy 2013, 51, 64–70. [Google Scholar] [CrossRef] [Green Version]
  35. Wan, C.; Wang, J.; Yang, G.; Gu, H.; Zhang, X. Wind farm micro-siting by Gaussian particle swarm optimization with local search strategy. Renew. Energy 2012, 48, 276–286. [Google Scholar] [CrossRef]
  36. Abdelsalam, A.M.; El-Shorbagy, M.A. Optimization of wind turbines siting in a wind farm using genetic algorithm based local search. Renew. Energy 2018, 123, 748–755. [Google Scholar] [CrossRef]
  37. Carrabs, F.; Cerrone, C.; Pentangelo, R. A multiethnic genetic approach for the minimum conflict weighted spanning tree problem. Networks 2019, 74, 134–147. [Google Scholar] [CrossRef]
  38. Shefaei, A.; Mohammadi-Ivatloo, B. Wild goats algorithm: An evolutionary algorithm to solve the real-world optimization problems. IEEE Trans. Ind. Inform. 2018, 14, 2951–2961. [Google Scholar] [CrossRef]
Figure 1. A simple example of Energy calculation when all WTs are operated in full load condition. (a) Energy distribution for the homogeneous wind farm; (b) Energy distribution for the OWFD.
Figure 1. A simple example of Energy calculation when all WTs are operated in full load condition. (a) Energy distribution for the homogeneous wind farm; (b) Energy distribution for the OWFD.
Energies 14 00644 g001
Figure 2. Voronoi diagram.
Figure 2. Voronoi diagram.
Energies 14 00644 g002
Figure 3. Voronoi diagram of a multiple wind turbine types offshore wind farms.
Figure 3. Voronoi diagram of a multiple wind turbine types offshore wind farms.
Energies 14 00644 g003
Figure 4. Encoding of particle i.
Figure 4. Encoding of particle i.
Energies 14 00644 g004
Figure 5. Decoding of particle (stage 1).
Figure 5. Decoding of particle (stage 1).
Energies 14 00644 g005
Figure 6. Decoding of particle (stage 2).
Figure 6. Decoding of particle (stage 2).
Energies 14 00644 g006
Figure 7. LSS for OWFD.
Figure 7. LSS for OWFD.
Energies 14 00644 g007
Figure 8. Optimization framework of V-APSO-LS.
Figure 8. Optimization framework of V-APSO-LS.
Energies 14 00644 g008
Figure 9. Optimized cable connection scheme. (a) Optimized layout of Model 1 obtained by V-APSO-LS; (b) Optimized layout of Model 2 obtained by V-APSO-LS.
Figure 9. Optimized cable connection scheme. (a) Optimized layout of Model 1 obtained by V-APSO-LS; (b) Optimized layout of Model 2 obtained by V-APSO-LS.
Energies 14 00644 g009
Figure 10. Optimized layout of Model 2 and highlighted the cables which changed the cross-section.
Figure 10. Optimized layout of Model 2 and highlighted the cables which changed the cross-section.
Energies 14 00644 g010
Figure 11. Ten trails for layout optimization of Model 1 and Model 2 obtained by V-APSO-LS.
Figure 11. Ten trails for layout optimization of Model 1 and Model 2 obtained by V-APSO-LS.
Energies 14 00644 g011
Figure 12. Ten trails for layout optimization of Model 2 obtained by V-APSO-LS, V-APSO, APSO-MST and V-WGA-LS.
Figure 12. Ten trails for layout optimization of Model 2 obtained by V-APSO-LS, V-APSO, APSO-MST and V-WGA-LS.
Energies 14 00644 g012
Figure 13. Optimized cable connection scheme. (a) Optimized layout obtained by V-APSO; (b) Optimized layout obtained by APSO-MST; (c) Optimized layout obtained by V-WGA-LS.
Figure 13. Optimized cable connection scheme. (a) Optimized layout obtained by V-APSO; (b) Optimized layout obtained by APSO-MST; (c) Optimized layout obtained by V-WGA-LS.
Energies 14 00644 g013
Figure 14. Total cost corresponding to each generation.
Figure 14. Total cost corresponding to each generation.
Energies 14 00644 g014
Table 1. Simulation results of Model 1 and Model 2.
Table 1. Simulation results of Model 1 and Model 2.
ItemModel 1Model 2
CAPEX (kEUR)18,983.4922,014.29
CELC (kEUR)16,666.4710,272.63
Total cost (kEUR)35,649.9632,286.92
Table 2. The best, mean and worst solutions obtained by different algorithms in 10 trials.
Table 2. The best, mean and worst solutions obtained by different algorithms in 10 trials.
ItemCAPEX (kEUR)CELC (kEUR)Total Cost (kEUR)
Best SolutionV-APSO-LS22,014.2910,272.6332,286.92
V-APSO21,081.6314,026.2335,107.86
APSO-MST19,791.8017,208.9737,000.77
V-WGA-LS21,869.7012,241.7734,111.47
Mean SolutionV-APSO-LS22,970.3410,643.3033,613.64
V-APSO21,015.1915,862.9936,878.18
APSO-MST19,622.7918,700.8038,323.59
V-WGA-LS22,065.2012,412.2934,477.49
Worst SolutionV-APSO-LS23,923.7510,514.2034,437.95
V-APSO22,823.7516,742.7439,566.49
APSO-MST19,360.4820,470.1039,830.58
V-WGA-LS22,324.5312,855.7235,180.25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qi, Y.; Hou, P.; Liu, G.; Jin, R.; Yang, Z.; Yang, G.; Dong, Z. Cable Connection Optimization for Heterogeneous Offshore Wind Farms via a Voronoi Diagram Based Adaptive Particle Swarm Optimization with Local Search. Energies 2021, 14, 644. https://doi.org/10.3390/en14030644

AMA Style

Qi Y, Hou P, Liu G, Jin R, Yang Z, Yang G, Dong Z. Cable Connection Optimization for Heterogeneous Offshore Wind Farms via a Voronoi Diagram Based Adaptive Particle Swarm Optimization with Local Search. Energies. 2021; 14(3):644. https://doi.org/10.3390/en14030644

Chicago/Turabian Style

Qi, Yuanhang, Peng Hou, Guisong Liu, Rongsen Jin, Zhile Yang, Guangya Yang, and Zhaoyang Dong. 2021. "Cable Connection Optimization for Heterogeneous Offshore Wind Farms via a Voronoi Diagram Based Adaptive Particle Swarm Optimization with Local Search" Energies 14, no. 3: 644. https://doi.org/10.3390/en14030644

APA Style

Qi, Y., Hou, P., Liu, G., Jin, R., Yang, Z., Yang, G., & Dong, Z. (2021). Cable Connection Optimization for Heterogeneous Offshore Wind Farms via a Voronoi Diagram Based Adaptive Particle Swarm Optimization with Local Search. Energies, 14(3), 644. https://doi.org/10.3390/en14030644

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