Next Article in Journal
Analysis of a Costly Fiberglass-Polyester Air Filter Fire
Previous Article in Journal
A Dual-Harmonic Pole-Changing Motor with Split Permanent Magnet Pole
Previous Article in Special Issue
A Novel Interval Programming Method and Its Application in Power System Optimization Considering Uncertainties in Load Demands and Renewable Power Generation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Convex-Optimization-Based Power-Flow Calculation Method for Offshore Wind Systems

1
PowerChina Huadong Engineering Corporation Limited, Hangzhou 311122, China
2
College of Electrical Engineering, Zhejiang University, Hangzhou 310027, China
*
Author to whom correspondence should be addressed.
Energies 2022, 15(20), 7717; https://doi.org/10.3390/en15207717
Submission received: 23 August 2022 / Revised: 22 September 2022 / Accepted: 28 September 2022 / Published: 19 October 2022

Abstract

:
Offshore wind farms have boomed worldwide due to the sustainability of wind power and ocean resources. Power grid companies should consider the wind power consumption problem with more power generated. Power-flow calculation is the most fundamental tool in energy management. This paper proposes the convex-relaxation-based method for offshore wind farms’ power flow. In this method, the traditional equations’ problem solving is transferred into standard convex optimization, which can be solved efficiently with unique optimum. Second-order cone relaxations are imposed to describe the quadratic relationship. The exactness of the relaxation is guaranteed with the special definition of the objective function.The superiority of the proposed method is tested on the case study, for which a computational efficiency improvement is shown. Moreover, the reliability of the power-flow results is verified within the precision tolerance.

1. Introduction

Power-flow analysis is to obtain the voltage and power information when the system’s topology is fixed. With this information, the power system’s planning, operation, and system-stability management states can be guided and decided [1]. Offshore wind energy resources are abundant and stable, and the global wind power elaboration shows a trend of developing from onshore to offshore [2].
Calculating the power flow of offshore wind systems helps to regulate and dispatch the energy. Since offshore wind systems are mostly in radial topology or weak-meshed topology, it is necessary to evaluate the methods of solving power flow. The main difficulty of the power-flow problem is solving a series of non-linear equations that describe the physical connection of the offshore wind system. The history of power flow can be traced back to the mid-20th century. In the early 1950s, power-flow calculation was carried out through the physical simulation of the DC computing table and the AC computing table. With the development of computer science, power-flow calculation transferred from analog technology to digital technology [3]. From then on, power-flow analysis developed with memory requirements, calculation speed and precision, and convergence. Thus, various kinds of power-flow calculation algorithms are mushrooming around the world [4].
Initially, the matrix-based Gauss–Seidel method [5] was used in power-flow calculation. This method only requires limited storage but shows a slow convergence rate. Later, the impedance matrix-based iterative method has been developed [6,7]. In comparison, the computational memory’s consumption of this method limits the system’s scale. As the power system becomes more and more complicated, this method cannot be applied anymore. At the same time, the Newton–Raphson method [8] has been created for power-flow analysis. Moreover, the PQ decomposition method [9], the artificial intelligence-based method [10], the backward/forward sweep method [11], and the approximated linear method [12] have also been applied in power-flow calculation. The stability of considered power flow has been proposed in [13]. The power-flow problem considering renewable energy has been proposed in [14]. However, the heuristic methods require much storage and computational time in the calculation [15]. The approximated linear equations will result in non-negligible errors [16], which is regarded as DC power flow. The Newton-series method is most sensitive to the initial value [17] since the essence of power-flow calculation is to solve a series of non-linear equations [18]. Sometimes, the solution of a specific power flow is not unique. Therefore, secondary processing is also required [19].
Another way to analyze power flow is to solve this through an optimization process since optimal power flow is the power flow solution with specific objectives. Among the optimization formats, convex optimization has strong points considering the character of the global optimum [20]. Therefore, the combination of convex optimization and power flow can guarantee a unique solution. Moreover, the convex optimization shows high efficiency in solving [21]. While the power flow constraints are non-linear, many relaxation methods should be implemented to relax the non-convex equations into convex ones, such as convex hull relaxations [22], convex-concave procedure technique [23], and semi-definite relaxation [24], second- order-cone relaxation [25] and so on. With the convex relaxation implemented, it is meaningful to check the relaxation gap. The large gaps will result in an infeasible optimal solution of the optimization before relaxations [26]. Therefore, to take advantage of convex optimization, exact relaxations should be guaranteed.
In this paper, the power-flow analysis of offshore wind farms is implemented. Convex optimization is combined with power-flow analysis, which helps to enhance the solving efficiency and guarantee a unique solution. The exactness of the relaxation is proved with the specific configuration of the objective function. The real wind farm system is tested in the case study to evaluate the performance of the proposed method with other methods, in which the superiority of the proposed method is demonstrated.
The main technical contributions made in this work can be summarized as follows:
(1)
A convex-relaxation-based method for power-flow calculation of offshore wind farms has been proposed.
(2)
The exactness of the relaxation can be guaranteed due to the design of the optimization objectives.
(3)
The unique solution of the power flow can be guaranteed with the strong points of the convex optimization.
This paper is organized as follows: In Section 2, the formulation of the power-flow problem is described. In Section 3, convex relaxations are imposed to convexify the power-flow equations, the exactness of which is proved mathematically. In Section 4, the proposed method is applied in a real case study. In Section 5, the main findings of the paper are summarized.

2. Power-Flow Calculation Problem

In this section, the original power flow will be demonstrated from the aspects of the notations, formulations, and descriptions.

2.1. Notations

The offshore wind system is a typical power system. In this paper, we make the following definitions.
  • Variables
    Define the voltage of node i with the complex variable V i , where V i = U i θ i . U i and θ i are the real variables, which denote the amplitude and the phase angle of node i, respectively. The active power and reactive power of node i are defined as the real variables P i and Q i , respectively. The apparent power of node i is defined as the complex variable S i . Similarly, the corresponding voltage and power variables are defined for node j as V j = U j θ j , S j = P j + i Q j . The injected current of node i is set as I i , where I i is the complex variable. For the transmission line ij, the transmitted active power and reactive power are shown with the real variables P i j and Q i j . For the apparent power of ij, the complex variables are S i j , where S i j = P i j + i Q i j .
  • Constants
    For the transmission line ij, the conductance and susceptance of the transmission line are defined as the real constants G i j and B i j , where the complex constant admittance is Y i j = G i j + i B i j . It is important to mention that the transmission line in this paper is set as the in-line equivalent circuit, where the susceptance to the ground is neglected.
  • Operator
    * is the conjugate operator of complex variables and constants.
The above notations are summarized in Figure 1.
A summary table of the notations has been added in the end of the paper.

2.2. Power Flow Formulation

For each node i in the power system, the relationship between the injected current and node’s voltage can be represented as
I i = j = 1 n Y i j V j ( i , j = 1 , 2 , . . . , n ) ,
where n is the total number of nodes in the system. This is the basic principle in the power system. The injected current of each node equals the total current of each connected node.
To solve the power-flow problem, the apparent power of the node and corresponding current should be used. Due to the physical principle, the relationship is shown as
I i = P i i Q i V i * ( j = 1 , 2 , , n ) ,
in mathematics. By substituting (2) into (1), the basic power flow formulation is obtained. It will be
P i i Q i V i * = j = 1 n Y i j V j ( j = 1 , 2 , , n ) .
In (3), there are n complex non-linear equations. In the power system, by constructing the equation with each node’s current, we can connect the branch’s information with the nodes’. The different ways to deal with this equation constructs various methods to solve power flow equations.

2.3. Description of Power Flow Analysis

In the power system, the operation state of each node i can be represented with a set of four variables: U i , θ i , P i , Q i . Therefore, for the power system with n nodes, there are 4 n variables in total. While not all of them are unknown, some of them is known before calculation.
To simplify the power-flow equations, the real parts and imaginary parts are decoupled individually. In the power-flow calculation (3), the n complex equations can be rewritten as 2 n real equations by decoupling the active power and reactive power, respectively. So, only 2 n variables can be obtained by solving power-flow equations with the other 2 n variables fixed. The nodes in the power system can be classified into three types: PQ node, PV node and slack node, according to the operation state. Whether the nodes’ information is known or unknown depends on the type of the nodes. More specifically, three types of nodes are defined as follows:
  • PQ node
    PQ nodes that characterize the known active power and the reactive power’s value. Additionally, the power-flow calculation aims to solve the node’s voltage amplitude and phase angle. The nodes connected with generators are always regarded as the PQ node in the power system. In offshore wind systems, the node of wind turbines is PQ nodes.
  • PV node
    The PV node characterizes the known active power and the phase angle’s value. Additionally, the power-flow calculation aims to solve the node’s voltage amplitude and reactive power. This node is embedded with reactive power sources to set the specific voltage amplitude in operation. The power plant buses with certain reactive power reserves are the PV nodes. In the offshore wind system, there are few nodes connected with loads. Sometimes, the storage units are bidden by the wind system. Thus, these nodes can be regarded as PV nodes. To guarantee the electricity’s quality of wind farms, there are regulations for the voltage limit of each node.
  • Slack node
    Only one node can be set as the slack node in power-flow calculation. This node regulates the reference phase angle of the voltage in the system. The node’s active and reactive power should be solved for the slack node. The slack node always poises the power balance of the whole system. Mostly, the frequency adjustment bus of the power plant is set as the slack node. In an offshore wind system, the grid node that transfers the energy between the power system and the offshore wind system can be regarded as the slack node.
To sum up, the power-flow analysis is used to solve Equation (3). In this equation, considering the features of the nodes PQ node, PV node, and slack node, there are 2 n known variables and 2 n unknown variables. The target of the power-flow analysis for the offshore wind system is to solve the unknown variables efficiently and accurately.

3. Convex Optimization Based Method

Since Equation (3) is non-convex, if the convex relaxation is implemented, the relaxations should be taken as precautions. In this section, the variables’ relaxation is defined firstly. Thereafter, the convex relaxation is imposed for the power-flow equations, which are used to convert the power-flow problem into a convex optimization.

3.1. Variables Relaxations

The first step to relaxing the non-linear equations is by variables’ relaxations. In the relationship between the injected power of each node and the voltage, there will be square variables of voltages. In the variables’ relaxation, we will focus on them. For two arbitrary connecting nodes ij in the offshore wind system, the variables x i , x j , x i j 1 , and x i j 2 are defined as
x i = U i 2 ,
x j = U j 2 ,
x i j 1 = U i U j sin ( θ i θ j ) ,
x i j 2 = U i U j cos ( θ i θ j ) .
The defined variables are the quadratic format of voltage variables, which describes the squared voltages of nodes, the real and imaginary of V i V j .
According to the definition above, we reconstruct some indirect variables. We define three variables x i j , y i j , and z i j , referring to the value of x i , x j , x i j 1 , and x i j 2 , that is
x i j = x i j 1 ,
y i j = x i x j ,
z i j = x i + x j 2 x i j 2 .
From the above equations, it can be known that the variables x i , x j , x i j 1 , and x i j 2 can be replaced by the linear representations of variables x i j , y i j , and z i j .
When comparing (4)–(7) with (8)–(10), it can be known that with the value of x i j , y i j , and z i j we can have endless sets of value x i , x j , x i j 1 , and x i j 2 as
1 1 0 1 1 2 x i x j x i j 2 = y i j z i j
When solving (11), we can see that there is an infinite set of solutions. By inducing an arbitrary variable ξ , we can obtain a set of solutions, for example,
x i = z i j 6 + y i j 2 + ξ i j x j = z i j 6 y i j 2 + ξ i j x i j 1 = x i j x i j 2 = z i j 3 + ξ i j
With different ξ , there will be various relationship representations between the variables x i , x j , x i j 1 , x i j 2 and x i j , y i j , and z i j .

3.2. Convex Relaxation

The power-flow equations can be rewritten with the relaxed variables x i j , y i j and z i j . In the power system, considering the relationship of the transmission line’s active power and reactive power with x i j , y i j , and z i j , we have
P i j P j i Q i j Q j i = B i j G i j 2 G i j 2 B i j G i j 2 G i j 2 G i j B i j 2 B i j 2 G i j B i j 2 B i j 2 x i j y j i z i j .
From Equation (13), the non-linear power-flow equations are converted into linear ones with the relaxed variables. This equation explains the power-flow equations from the aspects of the active power and the reactive power. The decoupling of this step converts the complex variables into real variables.
For the voltage constraints of the real operation state, we have the corresponding constraints for x i j , y i j , and z i j . In the offshore wind systems, we have some limitations for each node’s voltages. That is, to ensure the effectiveness of the solution, both the voltage’s amplitude and phase angle are limited. That is,
U i min U i U i max ,
U j min U j U j max ,
θ i j min θ i θ j θ i j max .
The superscript min and max are the lower and upper bounds of the variables, which are fixed as the constants due to the stable operation of the system. With the ranges in (14)–(16), the value of variables x i j , y i j , and z i j can be limited as well.
min ( U i U j sin ( θ i θ j ) ) x i j max ( U i U j sin ( θ i θ j ) ) ,
min ( U i 2 U j 2 ) y i j max ( U i 2 U j 2 ) ,
min ( U i 2 + U j 2 2 U i U j cos ( θ i θ j ) ) z i j max ( U i 2 + U j 2 2 U i U j cos ( θ i θ j ) ) .
The min and max symbol mentioned above are the minimize and maximize function. With the fixed value ranges in (14)–(16), the minimization and maximization will be fixed with constants accordingly.
Moreover, due to the power balance law, the active power and reactive power of the node should be equal to the sum of the total power from connected branches.
P i = i j E P i j ,
Q i = i j E Q i j .
where E denotes the set of all the branches in the power system.
From (4)–(7), we can know that there is the following relationship among variables x i , x j , x i j 1 , x i j 2 .
x i x j = x i j 1 2 + x i j 2 2
This is a quadratic constraint with squared terms and the product term. This equation constructs the non-linearity and non-convexity of the power-flow problem. Thus, convex relaxations are imposed as follows.
To make this optimization problem a convex optimization problem, we relax the equality constraints into inequality constraints as a second-order cone format. By substituting (12) into convex-relaxed (22), (22) can be represented as
x i j 2 + 1 4 y i j 2 + 1 12 z i j 2 z i j ξ i j
A schematic diagram of the second-order-cone relaxations is shown in Figure 2. The original constraints in (22) can be described as a surface of cones, as shown in the left picture in Figure 2. With relaxations, the constraints’ region is enlarged, as shown in the whole space contained in the cone in the right picture in Figure 2. Obviously, when the optimal solution is obtained on the surface of the cone, we can consider that the relaxation is exact. That is to say, although the solution is obtained by relaxation, the optimal solution is the same as the one obtained before relaxation.
Thus, the power-flow problem can be described as
min i j E ξ i j   subject to ( 13 ) , ( 17 ) ( 21 ) , ( 23 )
with convex relaxations.
This optimization problem becomes a problem relating to the variables x i j , y i j , z i j and ξ i j . Since ξ i j is a variable with a random possible value, this relaxation always tends to be exact. The exactness of the relaxation will be discussed in the next section.
In Equation (24), both the equality and inequality constraints are in a convex format. The objective function is a linear representation of the variable ξ i j , which is also a convex function. Thus, optimization (24) is a standard convex optimization problem. Due to the strong points of convex optimization, (24) can be solved efficiently with one unique optimum. That is, only the power flow solution, which meets the power system’s physical requirements, can be obtained with (24). Compared with the method in [23], the (24) shows the method for power-flow calculation, in which the second-order cone relaxation is concluded, while in [23], the method for power flow optimization is proposed, whose aim is to optimize the power loss. To apply the method in [23] for power-flow calculation, the objective function should be changed as a constant. Moreover, in [23], the convex relaxation is based on the convex-concave procedure technique.

3.3. Exactness of Relaxations

In Section 3.2, it can be shown that the power-flow problem can be transferred into a convex optimization with relaxations. Since the convex relaxations have been implemented, it is essential to discuss the exactness. The exactness refers to the fact that the optimal solution is solved when (23) has an equality operator. That is the same as the one before convex relaxations. If the relaxation is not exact, the optimal value of the objective function could not be obtained.
So in this section, we try to prove the exactness of the relaxation mathematically. The following theorem is given.
Theorem 1.
For ξ i j with any value, the second-order-cone relaxation (23) will always be exact.
Proof of Theorem 1. 
The proof is carried out by the contradiction.
Assume after convex optimization (24), a ξ i j * is obtained on the transmission line ij when
x i j 2 + 1 4 y i j 2 + 1 12 z i j 2 < z i j ξ i j * .
Therefore, the gap of (25) can be recorded as ϵ i j * , where
ϵ i j * = z i j ξ i j * x i j 2 1 4 y i j 2 1 12 z i j 2 .
Due to (25), it can be known that
ϵ i j * > 0 .
So there must be an arbitrary constant ϵ i j , where
ϵ i j * > ϵ i j 0 .
In the power-flow calculation, ϵ i j is a defined variables with no real physical meaning. As shown in the optimization, there are no specific limitations for this. Therefore, both of ϵ i j * and ϵ i j can be the optimized solution. Since in convex optimization (24), the local optimum will be the global one. Thus, there will be only one unique solution, and the optimum of ϵ i j * is refuted.
That implies that, unless the exactness of the relaxation is obtained, the optimal solution of convex optimization will not be obtained. This finishes the proof. □
Since the exactness of the relaxations can be ensured, with the value of ξ , x i j , y i j , and z i j from the optimal solution in (24), the value of voltage variables x i , x j , x i j 1 , and x i j 2 can be recovered with (12).

3.4. Solving Process

The whole solving process of the proposed convex-relaxation-based method for power flow can be summarized in the following Algorithm 1.
Algorithm 1: Solve power-flow problem by the proposed method
Input: ( P i , Q i ) of the PQ nodes; ( U i , P i ) of the PV nodes
Output: ( U i , θ i ) of the PQ nodes; ( Q i , θ i ) of the PV nodes
1. Define variables x i , x j , x i j 1 , x i j 2 ; define variables x i j , y i j , z i j ;
2. Introduce the variable ξ i j ;
 represent x i , x j , x i j 1 , x i j 2 with x i j , y i j , z i j and ξ i j ;
3. Construct the optimization problem with limited ranges and set the objective function;
4. Relax the non-convex equality constraints with second-order-cone relaxation;
5. Solve the optimization problem and recover the electrical variables from x i , x j , x i j 1 , x i j 2 ,
x i j , y i j , z i j and ξ i j .
From the above, the whole process can be divided into five steps.
  • Step 1: The first step focuses on the variables’ relaxations. As stated in (4)–(7), the variables x i , x j , x i j 1 , and x i j 2 are defined. Afterwards, the variables x i j , y i j , and z i j are defined according to x i , x j , x i j 1 , and x i j 2 in (8)–(10).
  • Step 2: The slack variable ξ i j is introduced in (12). Thus, the variables x i , x j , x i j 1 , and x i j 2 can be represented with x i j , y i j , z i j , and ξ i j .
  • Step 3: With the newly-defined variables, the power-flow problem can be represented in an optimization format, where the defined variables should be limited according to the ranges of voltage. Moreover, the power flow equation should be satisfied as an linear constraint. The objective function is set with the slack variable ξ i j .
  • Step 4: Relax the non-convex constraint (22) with the second-order-cone relaxation as (23).
  • Step 5: Solve the convex optimization problem (24) and obtain the value of x i j , y i j , z i j and ξ i j . Then recover x i , x j , x i j 1 , and x i j 2 with (12). It needs to be mentioned that the active power and reactive power can be obtained from the optimization (24). Additionally, there is no need to recover them.

4. Case Study

In this section, a real offshore wind farm will be implemented as the case study. In this wind farm, there are 103 wind turbines in total. The topology of the wind farm is shown in Figure 3, which is radial.
The Arabic numerals shown in Figure 3 denote the number of wind turbines. The differences between cables come from the distance to the offshore boosters. The nearer to the booster, the more over-current should be considered. The information on the connected wind turbines and their offshore submarine cables is shown in Table 1.
To test the effect of the proposed method, we evaluate the results from these aspects: 1. the exactness of the relaxations; 2. a comparison of the solving efficiency and solving accuracy with other methods.
The simulation is implemented by the Matlab program in the computer with an Intel 7th gen CPU. By constructing the optimization problem (24) in Yalmip [27], the power-flow problem can be solved by the optimizer toolbox Cplex [28].
Firstly, to test the exactness of the relaxation, for each offshore cable, we define the relaxation gap as
ϵ i j = | x i j 2 + 1 4 y i j 2 + 1 12 z i j 2 z i j ξ i j | .
By performing sampling calculation, we randomly select ten branches each time and calculate the average gap value of their gaps by
ϵ i j ¯ = | x i j 2 + 1 4 y i j 2 + 1 12 z i j 2 z i j ξ i j | 10 .
The values of ϵ i j ¯ changing in each sample are shown in the Figure 4.
From this figure, we can know that by setting the calculation accuracy to 10 6 , each sample of the test meets the requirement. Although the gaps vary in the test, all of the gaps are very small. By carefully checking the gap of each branch, it can be found that all of them are less than 10 6 . The accuracy standard is set as 10 6 in the solver considering the computational efficiency.
Secondly, to test the solving efficiency, we compare the proposed method with the Gauss–Seidel Method. The test results are shown in Table 2. The iteration times of the algorithm and the computational time are compared.
From the comparison, it can be found that the proposed method requires fewer iterations and thus has higher computational efficiency. The solution results of each node and branch’s power flow have been compared as well, which demonstrates the same solution.
Moreover, we extend the case study into a bigger wind farm, which is a combination of three proposed wind farms whose topology is shown in Figure 5.
In the case study, we calculate the power flow of this big system with +the Gauss–Seidel Method and the proposed method. To quantify the solution differences of these two methods, we define the variable τ as
τ = i = 1 , 2 , 3 . . . k | P i proposed method P i Gauss Seidel method | + | Q i proposed method Q i Gauss Seidel method | .
In this equation, k indicates the total nodes of the whole system. The result of the case study is shown in the Table 3.
From the Table 3, it can be found that when the system becomes bigger, the solving efficiency is improved. While the proposed method is still faster than the Gauss–Seidel method, the value of τ equals zero under the tolerance of 10 6 , which demonstrates the same electrical solution of the two methods. In conclusion, the proposed convex relaxation shows exact relaxations in the power-flow calculation, which guarantees the correctness of the solution. Moreover, the solving efficiency is relatively high compared with the existing method.

5. Conclusions

In this paper, a convex-optimization-based method for power-flow calculation of offshore wind systems has been proposed. The convex relaxation is implemented with an arbitrary defined variable, of which the exactness is proved. The method is applied in a real offshore wind farm. In the case study, both the solving efficiency and the correctness of the solution are verified. Additionally the superiority of the method is performed. The main findings of this manuscript can be summarized as
  • The computational time of the proposed method is significantly improved compared with the Gauss–Seidel method.
  • The average gap of the proposed convex relaxation is within the set accuracy standard.
In our future work, the proposed method will be compared with more existing methods, such as the fast decoupled load-flow method, the direct-current load-flow method, and the deep-learning-based artificial intelligence method. Moreover, technical improvements should be implemented to boost the computational efficiency further.

Author Contributions

Conceptualization, Y.C. and Q.Y.; methodology, Y.C.; software, Y.C. and H.L.; validation, Y.C., H.X., H.Q. and H.L.; formal analysis, Y.C. and Q.C.; writing—original draft preparation, Y.C.; writing—review and editing, H.L., H.Q. and H.X.; and supervision, Q.Y. and Q.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded National Key R&D program of China (2021YFB2400600) and Zhejiang Postdoctoral Preferential Funding Project (202211).

Conflicts of Interest

Yuwei Chen, Haifeng Qi, Hongke Li, Han Xu, and Qing Chen are employed by the PowerChina Huadong Engineering Corporation Limited. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Notations

The following notations are used in this manuscript:
V i = U i θ i Voltage of node i
U i Amplitude of V i
θ i Phase angle of the V i
P i Active power of node i
Q i Reactive power of node i
S i = P i + i Q i Apparent power of node i
I i injected current of node i
P i j Transmitted active power of the transmission line ij
Q i j Transmitted reactive power of the transmission line ij
S i j = P i j + i Q i j Transmitted apparent power of the transmission line ij
G i j Conductance of the transmission line ij
B i j Susceptance of the transmission line ij
Y i j = G i j + i B i j Admittance of the transmission line ij
*The conjugate operator

References

  1. Athay, T. An Overview of power-flow analysis. In Proceedings of the 1983 American Control Conference, San Francisco, CA, USA, 22–24 June 1983; pp. 404–410. [Google Scholar] [CrossRef]
  2. Esteban, M.D.; Diez, J.J.; López, J.S.; Negro, V. Why offshore wind energy? Renew. Energy 2011, 36, 444–450. [Google Scholar] [CrossRef] [Green Version]
  3. Ward, J.B. Equivalent circuits for power-flow studies. Electr. Eng. 1949, 68, 794. [Google Scholar] [CrossRef]
  4. Zhang, W.; Tolbert, L. Survey of reactive power planning methods. In Proceedings of the IEEE Power Engineering Society General Meeting, San Francisco, CA, USA, 12–16 June 2005; Volume 2, pp. 1430–1440. [Google Scholar] [CrossRef] [Green Version]
  5. Ahmadi, A.; Manganiello, F.; Khademi, A.; Smith, M.C. A Parallel Jacobi-Embedded Gauss–Seidel Method. IEEE Trans. Parallel Distrib. Syst. 2021, 32, 1452–1464. [Google Scholar] [CrossRef]
  6. Brown, H.; Carter, G.; Happ, H.; Person, C. Power flow solution by impedance matrix iterative method. IEEE Trans. Power Appar. Syst. 1963, 82, 1–10. [Google Scholar] [CrossRef]
  7. Cheng, C.S.; Shirmohammadi, D. A three-phase power flow method for real-time distribution system analysis. IEEE Trans. Power Syst. 1995, 10, 671–679. [Google Scholar] [CrossRef]
  8. Sereeter, B.; Vuik, C.; Witteveen, C. On a comparison of Newton–Raphson solvers for power-flow problems. J. Comput. Appl. Math. 2019, 360, 157–169. [Google Scholar] [CrossRef]
  9. Ao, L.; Cheng, B.; Li, F. Research of power flow parallel computing based on MPI and PQ decomposition method. In Proceedings of the 2010 International Conference on Electrical and Control Engineering, Wuhan, China, 25–27 June 2010; IEEE: Piscataway, NJ, USA, 2010; pp. 2925–2928. [Google Scholar] [CrossRef]
  10. Yan, Z.; Xu, Y. Real-time optimal power flow: A lagrangian based deep reinforcement learning approach. IEEE Trans. Power Syst. 2020, 35, 3270–3273. [Google Scholar] [CrossRef]
  11. Shen, T.; Li, Y.; Xiang, J. A graph-based power flow method for balanced distribution systems. Energies 2018, 11, 511. [Google Scholar] [CrossRef] [Green Version]
  12. Tang, K.; Dong, S.; Song, Y. Successive-intersection-approximation-based power flow method for integrated transmission and distribution networks. IEEE Trans. Power Syst. 2020, 35, 4836–4846. [Google Scholar] [CrossRef]
  13. Veerasamy, V.; Wahab, N.I.A.; Ramachandran, R.; Othman, M.L.; Hizam, H.; Devendran, V.S.; Irudayaraj, A.X.R.; Vinayagam, A. Recurrent network based power flow solution for voltage stability assessment and improvement with distributed energy sources. Appl. Energy 2021, 302, 117524. [Google Scholar] [CrossRef]
  14. Khan, I.U.; Javaid, N.; Gamage, K.A.; Taylor, C.J.; Baig, S.; Ma, X. Heuristic algorithm based optimal power flow model incorporating stochastic renewable energy sources. IEEE Access 2020, 8, 148622–148643. [Google Scholar] [CrossRef]
  15. Schloss, P.D.; Westcott, S.L. Assessing and improving methods used in operational taxonomic unit-based approaches for 16S rRNA gene sequence analysis. Appl. Environ. Microbiol. 2011, 77, 3219–3226. [Google Scholar] [CrossRef] [Green Version]
  16. Wood, S.N. Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J. R. Stat. Soc. Ser. B Stat. Methodol. 2011, 73, 3–36. [Google Scholar] [CrossRef] [Green Version]
  17. Sadjina, S.; Kyllingstad, L.T.; Skjong, S.; Pedersen, E. Energy conservation and power bonds in co-simulations: Non-iterative adaptive step size control and error estimation. Eng. Comput. 2017, 33, 607–620. [Google Scholar] [CrossRef] [Green Version]
  18. Sur, U.; Sarkar, G. Existence of Explicit and Unique Necessary Conditions for Power Flow Insolvability in Power Distribution Systems. IEEE Syst. J. 2019, 13, 702–709. [Google Scholar] [CrossRef]
  19. Park, S.; Zhang, R.Y.; Lavaei, J.; Baldick, R. Uniqueness of Power Flow Solutions Using Monotonicity and Network Topology. IEEE Trans. Control Netw. Syst. 2021, 8, 319–330. [Google Scholar] [CrossRef]
  20. Boyd, S.; Vandenberghe, L.; Faybusovich, L. Convex Optimization. IEEE Trans. Autom. Control 2006, 51, 1859. [Google Scholar] [CrossRef]
  21. Bertsekas, D. Convex Optimization Theory; Athena Scientific: Nashua, NH, USA, 2009; Volume 1. [Google Scholar]
  22. Goodarzi, S.; Gitizadeh, M.; Reza Abbasi, A. Efficient linear network model for TEP based on piecewise McCormick relaxation. IET Gener. Transm. Distrib. 2019, 13, 5404–5412. [Google Scholar] [CrossRef]
  23. Chen, Y.; Li, H.; Chen, Q.; Xie, R.; Wang, X. Convex-Concave Procedure Based Method for Optimal Power Flow of Offshore Wind Farms. Front. Energy Res. 2022, 10, 1056. [Google Scholar] [CrossRef]
  24. Bai, X.; Wei, H.; Fujisawa, K.; Wang, Y. Semidefinite programming for optimal power-flow problems. Int. J. Electr. Power Energy Syst. 2008, 30, 383–392. [Google Scholar] [CrossRef]
  25. Chen, Y.; Xiang, J.; Li, Y. Socp relaxations of optimal power-flow problem considering current margins in radial networks. Energies 2018, 11, 3164. [Google Scholar] [CrossRef] [Green Version]
  26. Chen, Y.; Xiang, J.; Li, Y. A quadratic voltage model for optimal power flow of a class of meshed networks. Int. J. Electr. Power Energy Syst. 2021, 131, 107047. [Google Scholar] [CrossRef]
  27. Johan, L. YALMIP: A toolbox for modeling and optimization in MATLAB. Skelet. Radiol. 2011, 41, 287–292. [Google Scholar] [CrossRef]
  28. Zuo, X.; Zhu, C.; Huang, C.; Xiao, Y. Using AMPL/CPLEX to model and solve the electric vehicle routing problem (EVRP) with heterogeneous mixed fleet. In Proceedings of the 2017 29th Chinese Control And Decision Conference (CCDC), Chongqing, China, 28–30 May 2017; pp. 4666–4670. [Google Scholar] [CrossRef]
Figure 1. The notations of power system.
Figure 1. The notations of power system.
Energies 15 07717 g001
Figure 2. The second-order-cone relaxation.
Figure 2. The second-order-cone relaxation.
Energies 15 07717 g002
Figure 3. The topology of the offshore wind farm.
Figure 3. The topology of the offshore wind farm.
Energies 15 07717 g003
Figure 4. The value of ϵ i j ¯ changing in each sample.
Figure 4. The value of ϵ i j ¯ changing in each sample.
Energies 15 07717 g004
Figure 5. The Topology of the offshore wind farm with a larger number of bars.
Figure 5. The Topology of the offshore wind farm with a larger number of bars.
Energies 15 07717 g005
Table 1. The information of the connected wind turbines and their offshore submarine cables.
Table 1. The information of the connected wind turbines and their offshore submarine cables.
Loop No.Wind Turbine No.Length of Cables (km)Cross-Section of Cables (mm)
11812,7843 × 300
21216,9823 × 300
3617,8963 × 300
43140103 × 500
54656253 × 500
67216,9873 × 300
77720573 × 500
82410,2473 × 300
93819873 × 500
103923673 × 500
115310,7823 × 500
126616,4423 × 300
138412,7833 × 500
146010,2573 × 300
159713,5783 × 300
1610311,2643 × 300
Table 2. Power flow results’ comparison.
Table 2. Power flow results’ comparison.
MethodComputational TimeIteration Times
Gauss–Seidel1.7 s486
Proposed Method0.68 s4
Table 3. Power flow results’ comparison.
Table 3. Power flow results’ comparison.
MethodComputational TimeIteration Times τ
Gauss–Seidel1.27 s384
Proposed Method0.82 s21
0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Y.; Qi, H.; Li, H.; Xu, H.; Yang, Q.; Chen, Q. Convex-Optimization-Based Power-Flow Calculation Method for Offshore Wind Systems. Energies 2022, 15, 7717. https://doi.org/10.3390/en15207717

AMA Style

Chen Y, Qi H, Li H, Xu H, Yang Q, Chen Q. Convex-Optimization-Based Power-Flow Calculation Method for Offshore Wind Systems. Energies. 2022; 15(20):7717. https://doi.org/10.3390/en15207717

Chicago/Turabian Style

Chen, Yuwei, Haifeng Qi, Hongke Li, Han Xu, Qiang Yang, and Qing Chen. 2022. "Convex-Optimization-Based Power-Flow Calculation Method for Offshore Wind Systems" Energies 15, no. 20: 7717. https://doi.org/10.3390/en15207717

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