Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Combined Effect of Ferrous Ion and Biochar on Cadmium and Arsenic Accumulation in Rice
Next Article in Special Issue
Multi-Objective Optimal Reactive Power Planning under Load Demand and Wind Power Generation Uncertainties Using ε-Constraint Method
Previous Article in Journal
Finite Element Method and Cut Bar Method-Based Comparison Under 150°, 175° and 310 °C for an Aluminium Bar
Previous Article in Special Issue
Field Optimization and Electrostatic Stress Reduction of Proposed Conductor Scheme for Pliable Gas-Insulated Transmission Lines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved Mixed AC/DC Power Flow Algorithm in Hybrid AC/DC Grids with MT-HVDC Systems

by
Fazel Mohammadi
1,*,
Gholam-Abbas Nazri
2 and
Mehrdad Saif
1
1
Electrical and Computer Engineering (ECE) Department, University of Windsor, Windsor, ON N9B 1K3, Canada
2
Electrical and Computer Engineering, Wayne State University, Detroit, MI 48202, USA
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(1), 297; https://doi.org/10.3390/app10010297
Submission received: 30 November 2019 / Revised: 18 December 2019 / Accepted: 19 December 2019 / Published: 31 December 2019
(This article belongs to the Special Issue Integration of High Voltage AC/DC Grids into Modern Power Systems)

Abstract

:
One of the major challenges on large-scale Multi-Terminal High Voltage Direct Current (MT-HVDC) systems is the steady-state interaction of the hybrid AC/DC grids to achieve an accurate Power Flow (PF) solution. In PF control of MT-HVDC systems, different operational constraints, such as the voltage range, voltage operating region, Total Transfer Capability (TTC), transmission reliability margin, converter station power rating, etc. should be considered. Moreover, due to the nonlinear behavior of MT-HVDC systems, any changes (contingencies and/or faults) in the operating conditions lead to a significant change in the stability margin of the entire or several areas of the hybrid AC/DC grids. As a result, the system should continue operating within the acceptable limits and deliver power to the non-faulted sections. In order to analyze the steady-state interaction of the large-scale MT-HVDC systems, an improved mixed AC/DC PF algorithm for hybrid AC/DC grids with MT-HVDC systems considering the operational constraints is developed in this paper. To demonstrate the performance of the mixed AC/DC PF algorithm, a five-bus AC grid with a three-bus MT-HVDC system and the modified IEEE 39-bus test system with two four-bus MT-HVDC systems (in two different areas) are simulated in MATLAB software and different cases are investigated. The obtained results show the accuracy, robustness, and effectiveness of the improved mixed AC/DC PF algorithm for operation and planning studies of the hybrid A/DC grids.

1. Introduction

Due to the recent developments in the power electronics technology, Voltage-Sourced Converter (VSC)-High Voltage Direct Current (HVDC) systems have solved the problem of bidirectional Power Flow (PF) in hybrid AC/DC grids [1,2,3]. MT-HVDC systems are capable of controlling the active and reactive power, independently. One of the important considerations to control the HVDC systems is that the V D C -control and P -control VSC-HVDC stations should be capable of operating in inverter and rectifier modes, respectively [4,5,6]. The main purpose of applying different control strategies in MT-HVDC systems is to achieve a precise and secure control mode for MT-HVDC systems without violating the operational constraints. For stable operation and active and reactive PF, MT-HVDC systems need to maintain the DC voltage and frequency within the operating limits [1,2,3,4].

1.1. AC/DC Power Flow for Hybrid AC/DC Grids with MT-HVDC Systems

There have been some relevant surveys about power system operation considering PF problem solutions [7,8]. In traditional AC systems, the PF can be controlled through three hierarchical levels which are, HLI (generation level), HLII (generation and transmission levels), and HLIII (generation, transmission, and distribution levels). Considering hierarchical levels and the market analysis, the Independent Electricity System Operator (IESO) can only control the generated power on the power system, and based on that, control the PF to the next hierarchical level(s). The integration of MT-HVDC systems to the existing AC grids leads to increasing the region of controllability of the hybrid AC/DC grids. This fact is due to the capability of controlling active and reactive PF by each converter station in MT-HVDC systems [1,4,9]. Therefore, MT-HVDC systems can change the PF patterns, and from the IESO and Transmission System Operator (TSO) perspectives, these changes in the pattern may cause significant issues in hybrid AC/DC grids. Moreover, due to the temporary or permanent outages of the components in hybrid AC/DC grids, the dynamic behavior of nonlinear components can change the PF and enhance the risk of instability in the entire hybrid AC/DC grids [9,10,11].

1.2. Solving the Power Flow Problem for Hybrid AC/DC Grids with MT-HVDC Systems

The unified method and sequential method are the two well-known methods to solve the PF problem for hybrid AC/DC grids with MT-HVDC systems.

1.2.1. The Unified Method

The unified method solves the PF problem for the entire hybrid AC/DC grids using a modified Jacobian matrix [12,13,14], where all the AC and DC variables, such as the impact of DC links in the Jacobian matrix, can be calculated in each iteration process. Many techniques are proposed and developed to improve the efficiency of the unified method, but the main drawback of those methods is neglecting the impact of droop parameters’ settings on the AC/DC PF [12,13,14]. Another issue of the unified method to solve the PF problem is that it needs an alternation of the extension of an existing AC PF.

1.2.2. The Sequential Method

The sequential method solves the AC/DC PF equations sequentially, one after another, in each iteration [15,16]. The main advantage of the sequential method is to make the solution easy to combine the DC PF to the AC PF solution, and it can be implemented easily when the extension of an existing AC PF is needed. In [17], a numerical method based on the Newton-Raphson algorithm to calculate the converter station losses is proposed. A detailed steady-state model of the converter station to solve the AC/DC PF problem sequentially based on the Newton-Raphson technique considering converter station losses and reference power is proposed in [18,19]. An algorithm with per-unit conversion and changing the bus numbers to simplify solving the PF problem when multiple DC lines and converter station outages occur is developed in [18]. The main drawback of this method is neglecting the AC grid connection in the problem. In [20], a detailed model of the converter station with AC/DC PF equations, including converter station losses, for solving the PF sequentially is proposed. In [21], a method to solve the PF problem through the Gauss-Seidel method is developed.

1.3. AC/DC Power Flow for MT-HVDC Systems Considering Droop Parameters

Controlling the droop parameters in MT-HVDC systems has a significant impact on the PF of AC and DC grids after an outage [4]. In [16], the concept of distributed DC voltage control with the droop parameters on the PF problem in MT-HVDC systems is considered. In [22], a methodology to determine the mean voltage instead of a single slack converter station in MT-HVDC systems to solve the DC PF is proposed. By interconnecting MT-HVDC systems to the large-scale AC grids, solving the AC/DC PF goes through a complicated process, and a combined solution of AC/DC PF considering all the system’s variables and constraints is required.
To address all the above-mentioned challenges, a mixed AC/DC PF algorithm for hybrid AC/DC grids with MT-HVDC systems is proposed in this paper. The proposed strategy is a fast and accurate method, which is capable of optimizing the AC/DC PF calculations. Except for the high accuracy and optimized performance, considering all operational constraints and control objectives of the integration of MT-HVDC systems into the large-scale AC grids is the other contribution of this paper. The calculated results by the mixed AC/DC PF problem can be used for the planning, scheduling, state estimation, small-signal stability analyses. The mixed AC/DC PF algorithm is applied to a five-bus AC grid with a three-bus MT-HVDC system and the modified IEEE 39-bus test system with two four-bus MT-HVDC systems (in two different areas) which are simulated in MATLAB software. To check the performance of the mixed AC/DC PF algorithm, different cases are considered. The obtained results show the accuracy, robustness, and effectiveness of the mixed AC/DC PF algorithm.

2. Principles of Power Flow in Power Systems

2.1. AC Grid Power Flow

The main objective of the AC PF is to determine the magnitude and angle at each bus of the AC grids. In order to analyze the AC PF in power systems, the following assumptions are considered.
  • The AC transmission networks have fast dynamics compared to the other components. In this regard, AC transmission networks can be represented by algebraic equations.
  • Each transmission line and transformer is modeled by an equivalent π model.
  • The power in AC grids is balanced.
  • The positive sequence parameters on a per-phase basis are assumed.
Accordingly, the network equations can be written as follows:
[ I ¯ 1 I ¯ 2 I ¯ i I ¯ n A C ] = [ Y ¯ 11 Y ¯ 21 Y ¯ i 1 Y ¯ n A C 1 Y ¯ 12 Y ¯ 22 Y ¯ i 2 Y ¯ n A C 2 Y ¯ 1 i Y ¯ 2 i Y ¯ i i Y ¯ n A C i Y ¯ 1 n A C Y ¯ 2 n A C Y ¯ i n A C Y ¯ n A C n A C ] [ V ¯ 1 V ¯ 2 V ¯ i V ¯ n A C ]
where I ¯ i and V ¯ i   = V i e j θ i are the injected current and voltage at the i th node, Y ¯ i i is the self-admittance at the i th node, Y ¯ i j is the mutual admittance between nodes i and j , and n A C represents the total number of buses in the AC grid.
The voltage of each node can be used for solving the equation of [ I ¯ ] [ Y b u s ¯ ] [ V ¯ ] = 0 . The injected current at the i th node is related to the injected power and bus voltage at that node and it can be calculated as follows:
I ¯ i = P i j Q i V ¯ i
Considering the operational constraints in power systems, the PF problem becomes nonlinear, and the nodes in AC grids can be classified into four types:
  • Slack bus where the voltage magnitude ( V i ) and angle ( θ i ) are determined.
  • PV bus where the active power ( P i ) injected to the grid and the voltage magnitude ( V i ) are known.
  • PQ bus where the active power and reactive power ( Q i ) injected into the grid are known.
  • MT-HVDC Point of Common Coupling (PCC) bus where different constraints based on the control mode of each converter station are applied.
From Equation (1), the injected current at the i th node (without considering the PCC buses of MT-HVDC systems) can be written as follows:
I ¯ i = j = 1 n A C Y i j ¯ V j ¯
Therefore, the active and reactive power injected to the grid can be derived as follows:
P i = j = 1 n A C V i V j [ G i j cos ( θ i θ j ) + B i j sin ( θ i θ j ) ]
Q i = j = 1 n A C V i V j [ G i j sin ( θ i θ j ) B i j cos ( θ i θ j ) ]
where Y ¯ i j = G i j + j B i j .
It should be noted that a generator can be either a slack bus or a PV bus, and a load can be a PQ bus with known active and reactive power.

2.2. MT-HVDC Systems Power Flow

The main aim of MT-HVDC systems PF is to determine the DC voltage magnitude at each converter station and the PF within the DC grid. Considering n converter stations in MT-HVDC systems, the following items should be considered to derive the PF equations for the i th converter station of MT-HVDC systems.
  • Interface of the converter station with the AC grid
  • AC side of the converter station
  • Interface of the converter station AC and DC sides
  • DC side of MT-HVDC systems
  • Control modes of the converter station

2.2.1. Interface of the Converter Station with the AC Grid

The PF equations of at the PCC bus of the i th converter station interfacing with the AC grid is as follows:
P g i = j = 1 n A C V g i V j [ G i j cos ( θ g i θ j ) + B i j sin ( θ g i θ j ) ]
Q g i = j = 1 n A C V g i V j [ G i j sin ( θ g i θ j ) B i j cos ( θ g i θ j ) ]
where V ¯ g i = V g i e j θ g i .

2.2.2. AC Side of the Converter Station

The PF equations of at the PCC bus of the i th converter station for the asymmetric bipolar MT-HVDC systems is as follows:
P g i P g i p P g i n = 0
Q g i Q g i p Q g i n = 0
where P g i p , P g i n , Q g i p , and Q g i n are the active and reactive power of the positive and negative poles of the i th converter station, respectively.
Therefore, the apparent power at the PCC bus of the i th converter station by the positive and negative poles can be written as follows:
S g i p = V ¯ g i I ¯ g i p * = P g i p + j Q g i p
S g i n = V ¯ g i I ¯ g i n * = P g i n + j Q g i n
Applying KCL at the PCC bus of the i th converter station leads to deriving the following equation.
I ¯ g i = I ¯ g i p + I ¯ g i p
Using KVL between the PCC bus and the converter terminal of the i th converter station leads to deriving the following equations.
I ¯ g i p = V ¯ t i p V ¯ g i Z i p = ( V ¯ t i p V ¯ g i ) ( G i p + j B i p )
I ¯ g i n = V ¯ t i n V ¯ g i Z i n = ( V ¯ t i n V ¯ g i ) ( G i n + j B i n )
where Z i p = R i p + j X i p and Z i n = R i n + j X i n . Also, V ¯ t i p = V t i p e j θ t i p and V ¯ t i n = V t i n e j θ t i n are the voltage of the positive and negative pole converter terminal, respectively.
Substituting Equations (13) and (14) into Equations (10) and (11), and also separating the active and reactive power, the following equations are obtained.
P g i p = V g i [ V g i G i p + V t i p { G i p cos ( θ g i θ t i p ) + B i p sin ( θ g i θ t i p ) } ]
Q g i p = V g i [ V g i B i p + V t i p { G i p sin ( θ g i θ t i p ) B i p cos ( θ g i θ t i p ) } ]
P g i n = V g i [ V g i G i n + V t i n { G i n cos ( θ g i θ t i n ) + B i n sin ( θ g i θ t i n ) } ]
Q g i n = V g i [ V g i B i n + V t i n { G i n sin ( θ g i θ t i n ) B i n cos ( θ g i θ t i n ) } ]

2.2.3. Interface of the Converter Station AC and DC Sides

Regardless of the converter station switching losses, the active power balance between converter station AC and DC sides is as follows:
P t i p = V D C i p I D C i p
P t i n = V D C i n I D C i n
where V D C i p , I D C i p , V D C i n , and I D C i n are the DC bus voltage and current of the positive and negative pole of the i th converter station.
Also, the active power at the AC-side terminal of the positive and negative pole of the converter stations can be written as follows:
P t i p = V t i p [ V t i p G i p V g i { G i p cos ( θ t i p θ g i ) + B i p sin ( θ t i p θ g i ) } ]
P t i n = V t i n [ V t i n G i n V g i { G i n cos ( θ t i n θ g i ) + B i n sin ( θ t i n θ g i ) } ]

2.2.4. DC Side of MT-HVDC Systems

When the current flowing through the DC-bus capacitor is zero, I D C i p = I L i p and I D C i n = I L i n , where I L i p and I L i n are the current of the positive and negative DC link connected to the i th h converter station. The injected current can be determined by the conductance matrix of the DC grid. The DC current injection at the positive and negative poles of the DC bus is as follows:
[ I D C 1 p I D C i p I D C N p I D C 1 n I D C i n I D C N n ] = [ G p p 11 G p p i 1 G p p N 1 G n p 11 G n p i 1 G n p N 1 G p p 1 i G p p i i G p p N i G n p 1 i G n p i i G n p N i G p p 1 N G p p i N G p p N N G n p 1 N G n p i N G n p N N G p n 11 G p n i 1 G p n N 1 G n n 11 G n n i 1 G n n N 1 G p n 1 i G p n i i G p n N i G n n 1 i G n n i i G n n N i G p n 1 N G p n i N G p n N N G n p 1 N G n p i N G n p N N ] [ V D C 1 p V D C i p V D C N p V D C 1 n V D C i n V D C N n ]
Hence, the injected current at the i th bus in the DC side is as follows:
I D C i p = j = 1 N ( G p p j i V D C j p + G p n j i V D C j n )
I D C i n = j = 1 N ( G n p j i V D C j p + G n n j i V D C j n )
where N represents the total number of buses in the DC grid. In addition, G p p j i , G n n j i , G p n j i , and G n p j i are the conductance of the positive and negative DC link between nodes i and j , respectively.

2.2.5. Control Modes of the Converter Station

There are different modes of operation for each converter station in MT-HVDC systems [4], and based on them, the steady-state equations in the form of equality constraints can be obtained. Also, there are some inequality constraints due to the limits, which are imposed by the converter station voltage and current ratings, as follows:
V D C i min V D C i V D C i max
V i min V i V i max
I g i I g i max
I D C i I D C i max

3. Mixed AC/DC Power Flow Algorithm

The mixed AC/DC PF algorithm, which is an improved sequential AC/DC PF algorithm [19], can be used to obtain the initial operating points to analyze the dynamics of the hybrid AC/DC grids by solving the AC and DC PF sequentially and keeping both the converter station power and voltage at each node constant.
In order to implement the mixed AC/DC PF algorithm, the per-unit conversion should be performed for the entire system. Each converter station is connected to both AC and DC grids. The AC side of the converter station is modeled by a voltage source connected to the AC bus through a phase reactor ( Z C = R C + j ω L C ), a capacitor ( Z F = j ω C ), and a transformer ( Z T R = R T R + j ω L T R ). Also, the DC side of the converter station is connected to the DC grid. The power losses at the converter station are considered as a quadratic function of the converter station Root Mean Square (RMS) AC current as follows [23]:
P L o s s = a + b I + c I 2
where a , b , and c are the loss coefficients.
It should be noted that,
P g i + P D C i + P L o s s i = 0
Assume an MT-HVDC system with n converter stations and n D C DC lines. For AC/DC PF analysis, at least one converter station should be capable of controlling V D C (DC-slack bus) in the entire DC grids, and the rest of the converter stations control the active power. The detailed explanations of the mixed AC/DC PF are given as follows:
  • Start by an initial guess of the active power injected to the AC grid by the DC-slack converter station at the l th iteration ( P g n s ( l ) ), where n s represents the index of the DC-slack bus.
  • Transform all converter stations which are connected to the j th bus to PV or PQ buses based on their control modes (including droop-based control strategies) and solve the AC PF. At the l th iteration, the active power injected by all the non-slack converter stations is constant, while P g n s ( l ) changes.
  • Calculate the converter station losses using Equation (30) considering the active power injected by the i th converter to the AC grid ( P g i ( l ) ) and the active power injected by the i th converter station to the DC grid ( P D C i ( l ) ). In this step, the AC/DC connections and the converter stations’ limits should be considered.
  • Solve the DC PF for the DC grid using the Newton-Raphson method. In this step, the DC-slack bus regulates the DC voltage ( V D C n s ( 0 ) ) initially, and the DC buses determine the active power injection ( P D C i ( l ) ). Therefore, the DC voltage at each bus and also P D C n s ( l ) are calculated.
  • Compute DC-slack and droop buses iteration ( k ). As a new value of P g n s ( l + 1 ) is calculated, considering the converter station losses,
    • Initialization: P g n s ( k = 0 ) = P g n s ( l ) .
    • Solve the branch j , n s g , n s considering V j n s , θ j n s , Q j n s , and P g n s ( k ) using the Newton-Raphson method.
    • Obtain the new value of P g n s ( k + 1 ) with P D C n s and P L o s s n s using Equations (30) and (31).
    • If | P g n s ( k + 1 ) P g n s ( k ) | < ε , stop the calculations. Otherwise, k = k + 1 and return to step II. The output is P j n s ( k + 1 ) .
  • Check the convergence criterion. If | P g n s ( k + 1 ) P g n s ( k ) | < ε , stop the calculations. Otherwise, k = k + 1 and return to step 1.
It should be noted that all linear and nonlinear variables are considered in the mixed AC/DC PF algorithm. Based on the topology of the hybrid AC/DC grids, the rating and length of each transmission line are determined. Also, the droop parameters and reference voltage have a direct impact on the reference power.

4. Results and Discussions

For validation and to demonstrate the performance of the mixed AC/DC PF algorithm, a five-bus AC grid with a three-bus MT-HVDC system is simulated in MATLAB software, as shown in Figure 1. The test system is composed of both AC and DC grids. The data of the system is provided in Table 1, Table 2, Table 3, Table 4 and Table 5.
The parameters of the converter stations for the PF calculations are as follows:
  • Converter Station #1: P Q control mode, P g = 60   MW and Q g = 40   MVAR
  • Converter Station #2: V D C V control mode, V D C = 1   p . u . and V = 1   p . u .
  • Converter Station #3: P Q control mode, P g = 35   MW and Q g = 5   MVAR
It should be noted that S B = 100   MVA , V D C B = 345   kV , and V B = 345   kV are the base values for the hybrid AC/DC grids per-unit system. Also, the links between the two nodes on the DC grids are bipolar.
The simulations are accomplished using a laptop with the Intel Core i7-8550U processor at 1.80 GHz clock speed and 12-GB of RAM.

4.1. Case 1: AC Power Flow without MT-HVDC Systems

In the first case, the AC PF calculations are performed without considering the DC system. The algorithm is converged in 0.44 s and in three iterations. Table A1 and Table A2 show the results of AC PF without DC grids.

4.2. Case 2: AC/DC Power Flow Considering MT-Systems with Constant Active Power and DC Voltage

In the second case, the algorithm is applied to solve the AC/DC PF problem for the studied system considering that converter stations #1 and #3 are operated in constant P -mode and converter station #2 is operated in constant V D C -mode. The algorithm is converged in 0.82 s and in three iterations. Table A3, Table A4, Table A5, Table A6 and Table A7 show the obtained results of the mixed AC/DC PF calculations in Case 2.
Compared to Case 1, in this case, both the active and reactive PF on each line have decreased and accordingly, the total power losses are decreased. It is observed that the generator connected to bus 2 injects its maximum active power. In addition, compared to the previous case, the total injected reactive power by the two generators is decreased.

4.3. Case 3: AC/DC Power Flow Considering Converter Station Outage

In the third case, the impact of the converter station outage on the AC/DC PF results is analyzed. Table A8, Table A9, Table A10 and Table A11 demonstrate the results of PF calculations both in all AC and DC buses and in branches, in case of the outage of converter station #1. Based on the obtained results, the converter station #3 is capable of operating in constant P -mode. It should be noted that the algorithm is converged in 0.55 s and in three iterations.
In this case, the generator connected to bus 2 injects 40 MW active power to the grid (maximum active power capability) and the total injected reactive power this generator is approximately tripled. The outage of converter station #1 is led to an increase in the active and reactive PF of the AC lines. Therefore, the total power losses on the AC grids are increased. Due to the fact that converter station #1 is not connected to the grid, the summation of the PF between the DC lines (1-2 and 1-3) is zero and consequently, the overall DC power losses are decreased.

4.4. Case 4: AC/DC Power Flow Considering Droop Control Strategy and Converter Outage

As stated in Section 1.3, the droop control strategy is an efficient way of controlling MT-HVDC systems that can improve the PF of the AC and DC grids after an outage. In the fourth case, it is assumed that all converter stations in the studied system are equipped with the droop controllers. Table A12 illustrates the detailed information of the droop settings for each converter station. Table A13, Table A14, Table A15 and Table A16 show the results of PF calculations both in all AC and DC buses and in branches in case of the outage of converter station #1 and considering droop control settings provided in Table A12. It should be noticed that the algorithm is converged in 0.45 s and in three iterations.
Compared to Case 3, in this case, the contribution of the generator connected to bus 2 to power generation is decreased. As the outage of converter station #1 is considered, compared to Case 3, the PF on the DC lines is decreased. It should be noted that as the performance of each converter station is based on the pre-determined droop parameter, the PF on the DC lines should be in a way that the total power losses on the DC lines become zero.

4.5. Case 5: AC/DC Power Flow Considering Changes of Droop Parameters and Converter Outage

In the fifth case, to demonstrate the impact of changes in the droop parameters of each converter station, the same test as Case 4 is evaluated by changing the droop parameters as shown in Table A17. Table A18, Table A19, Table A20 and Table A21 depict the results of PF calculations both in all AC and DC buses and in branches in case of the outage of converter station #1 and considering the changes in droop parameters according to Table A17. It is worth mentioning that the algorithm is converged in 0.43 s and in three iterations.
Changing the droop parameters is led to an increase in the total active power generated by the generator connected to bus 2. As a consequence, there are some slight changes of the PF on the AC lines. However, compared to the previous case, those changes are negligible. In addition, changing the droop parameters is causes an increase in the DC PF on the DC lines (compared to Case 4), but as the converter station #1 is not connected to the grid, the summation of the PF between the DC lines is equal to zero.

4.6. Case 6: AC/DC Power Flow Considering Converter Station Limits and Converter Outage

In the sixth case, the impact of the converter station limits on the AC/DC PF solution is investigated. When an active power set-point of a converter station equipped with the P controller is outside of the P-Q capability chart (P-Q capability chart shows the possible operation points.), the active power order should be reduced to comply with the predefined limit. Similarly, when a reactive power set-point of a converter station equipped with the Q controller is outside of the P-Q capability chart, the reactive power order should be reduced to comply with the predefined limit, subject to not reaching the active power limit. In addition, when a converter station equipped with the V D C -droop controller reaches its limit, the converter station should be set to a constant P injection equal to the maximum active power limit of the converter station. In addition, when a converter station equipped with the V -droop controller reaches its limit, the converter station should be set to a constant Q injection based on the predefined limit. For both converter stations equipped with the Q controller and V D C -droop controller, the priority is given to active power over reactive power, when enforcing the limits. It should be noted that however, all the DC-slack buses are disregarded from the analysis, they are rechecked at the end of AC/DC PF calculations.
According to the above explanations, it is assumed that the converter station #1 reaches the reactive power limit and the converter station control is changed from constant V -mode to constant Q -mode. Table A22, Table A23, Table A24 and Table A25 illustrate the results of PF calculations in all both AC and DC buses and branches after enforcing the converter station’s current and voltage limits. Meanwhile, the algorithm is converged in 0.60 s and in three iterations. After three iterations, the reactive power set-point of the converter station #1 is changed from −40 MVAR to −36.25 MVAR.
To demonstrate the impact of changing the set-point of the converter station, the active power set-point of the converter station #1 is set to −130 MW so that simultaneously both active and reactive power violate their limits. Table A26, Table A27, Table A28 and Table A29 show the results of PF calculations both in all AC and DC buses and in branches after changing the active power set-point of the converter station #1 to −130 MW. It should be noted that the algorithm is converged in 0.71 s and in three iterations.
After three iterations, the active and reactive power set-points of the converter station #1 are changed from −130 MW to −121.21 MW, and −40 MVAR to −8.96 MVAR, respectively. In this case, the total active power generated by the generation units is increased (tangible changes are related to the one connected to bus 2). The overall power losses on the AC lines are decreased. As the converter station #1 is connected to the grid and its set-points are reached their maximums limits, the active power losses (correspond to Z I 2 ) are increased. Therefore, the PF on the DC lines is increased.

4.7. Case 7: AC/DC Power Flow for the Large-Scale Hybrid AC/DC Grids with MT-HVDC Systems

In order to validate and show the performance of the mixed AC/DC PF algorithm for the large-scale hybrid AC/DC grids with the integration of MT-HVDC system, the modified IEEE 39-bus test system with two four-bus MT-HVDC systems (with different colors) is simulated in MATLAB software, as shown in Figure 2. The data of the MT-HVDC systems are provided in Table A30, Table A31 and Table A32.
As it can be observed from Figure 2, the IEEE 39-bus test system is divided into three areas, in which area 1 comprises of three coherent generators (G8, G9, and G10), area 2 comprises of three coherent generators (G1, G2, and G3), and area 3 comprises of four coherent generators (G4, G5, G6, and G7). The main corridors to interconnect those three areas are given as follows:
  • Area 1 to Area 2: Lines from bus 1 to bus 39, and from bus 3 to bus 4
  • Area 1 to Area 3: Lines from bus 3 to bus 18, and from bus 27 to bus 17
  • Area 2 to Area 3: Line from bus 14 to bus 15
The parameters of the converter stations for the PF calculations are as follows:
  • Converter Station #1: P Q control mode, P g = 60   MW and Q g = 20   MVAR
  • Converter Station #2: V D C V control mode, V D C = 1   p . u . and V = 1   p . u .
  • Converter Station #3: P Q control mode, P g = 40   MW and Q g = 40   MVAR
  • Converter Station #4: P Q control mode, P g = 40   MW and Q g = 40   MVAR
  • Converter Station #5: V D C V control mode, V D C = 1   p . u . and V = 1   p . u .
  • Converter Station #6: P Q control mode, P g = 40   MW and Q g = 30   MVAR
  • Converter Station #7: P Q control mode, P g = 40   MW and Q g = 20   MVAR
  • Converter Station #8: P Q control mode, P g = 40   MW and Q g = 40   MVAR
From the operation and planning perspectives, the outage of the mentioned lines can cause entirely disconnection of the two areas from each other and accordingly, prevention of power exchange between two areas, operating in islanded mode, and instability of the AC grids. Therefore, there is a need for strengthening the power transmission lines among the mentioned areas. To do so, two four-bus MT-HVDC systems are considered to interconnect the mentioned areas together. Areas 1 and 3 are interconnected with four converter stations (CS1, CS2, CS3, and CS4) and areas 2 and 3 are interconnected with four converter stations (CS5, CS6, CS7, and CS8). Hence, in case of AC lines outage, the power can be transferred via DC links. Table A33 and Table A34 show the results of AC PF on the IEEE 39-bus test system without DC grids (Case 7-1).
Table A35, Table A36, Table A37 and Table A38 show the obtained results of the mixed AC/DC PF calculations when MT-HVDC systems are connected to the IEEE 39-bus test system (Case 7-2). In a general view, due to the fact that by interconnecting the MT-HVDC systems to the grid, more power is required, the total generated active and reactive power by the generators should increase and decrease, respectively. There is no connection between the MT-HVDC systems and any of the generators. Therefore, only the voltage magnitude and angles are changed. This should be considered that each converter has its own power losses. Therefore, compared to the AC PF, the total power losses on the AC lines in the mixed AC/DC PF is increased by 0.11 MW. There are some slight changes in the total active power generated by the generator #31. However, the generated reactive power by all the generators is changed. The changes of the PF on the AC lines 1–39, 3–4, 3–18, 4–14, 14–15, 15–16, and 17–27 are significant since they are directly connected to the MT-HVDC systems. To check the changes of the PF on the other lines, the connected lines to the mentioned lines are checked and it is observed that there are some slight changes of the PF on the lines 4-5 and 13-14. Based on the obtained results, the slack converters are capable of transferring power on the DC lines based on the total demand. It is also noticed that compared to the other DC lines, the changes of the PF on the lines 5–6, 5–7, and 5–8 are considerable and the changes of the PF on the DC lines 6–7 and 6–8 are almost zero.
To minimize the power losses, in the next case, it is considered that the converters are equipped with the droop controllers. Table A39 illustrates the detailed information of the droop settings for each converter station. Table A40, Table A41, Table A42 and Table A43 illustrate the results of mixed AC/DC PF considering droop control settings provided in Table A39 (Case 7-3). Analyzing the obtained results shows that the changes of the PF on the AC lines 1–39, 2–3, and 17–27 are almost negligible but the total power losses on the AC lines are decreased. The DC lines 1–4 and 2–3 have the most power losses. In addition, no power is transferred on lines 6–7 and 6–8. Therefore, they can be considered as reserve lines for planning purposes. It is worth mentioning that compared to Case 7-2, the total power losses on the DC lines are decreased.
To analyze the impact of MT–HVDC systems on the PF analysis, AC lines 3–4, 3–18, 4–14, and 14–15 as the main interconnected corridors between each of area are disconnected. Table A44 and Table A45 show the results of AC PF without DC grids in the case of disconnecting the AC lines 3–4, 3–18, 4–14, and 14–15 (Case 7-4). The disconnection of the mentioned lines is caused that some of the lines reach their maximum transfer power capabilities. The changes of the PF on the AC lines 4–5, 5–6, 6–11, and 10–11 are drastically changed. However, the changes of the PF on the AC lines 10–13 and 13–14 are deceased. The rest of the lines have either no or very slight changes on the PF.
To improve the reliability of the hybrid AC/DC grids, the MT-HVDC systems are connected to compensate for the disconnections of the AC lines 3–4, 3–18, 4–14, and 14–15. Table A46, Table A47, Table A48 and Table A49 show the obtained results of the mixed AC/DC PF calculations when MT-HVDC systems are connected to the IEEE 39-bus test system considering the disconnections of the AC lines 3–4, 3–18, 4–14, and 14–15 (Case 7-5). Based on the obtained results, the significant changes on the PF are related to the AC lines 1–2, 1–39, 10–13, 13–14, and 17–27 (increase) and 2–25, 5–6, 6–11, 10–11, 15–16, 17–18, and 26–27 (decrease). The total power losses on the AC lines are closed to the previous case. In addition, DC lines 5–6, 5–7, and 5–8 are the DC lines with high power losses.
In order to minimize the power losses, it is assumed that the converter stations are equipped with the droop controllers as Table A39. Table A50, Table A51, Table A52 and Table A53 show the obtained results of the mixed AC/DC PF calculations when MT-HVDC systems are connected to the IEEE 39-bus test system considering the droop parameters and disconnections of the AC lines 3–4, 3–18, 4–14, and 14–15 (Case 7-6). The obtained results show that compared to the previous case, only the PF of the AC lines 4–5, 5–6, 15–16, 17–18, and 26–27 are increased and the total power losses both on the AC and DC lines are decreased.

5. Conclusions

In this paper, a mixed AC/DC Power Flow (PF) algorithm for the steady-state interaction of the large-scale MT-HVDC systems is investigated. This algorithm is an improved sequential AC/DC PF algorithm, which uses the Newton-Raphson method to solve the DC PF problem. Different operational constraints and control strategies along with contingency analysis in the hybrid AC/DC grids are considered in this study. Fast convergence and high accuracy are the main advantages of the mixed AC/DC PF algorithm. In addition, it is a powerful tool for sensitivity analysis and congestion management in power systems. Various cases are studied in this paper to evaluate the performance of the mixed AC/DC PF algorithm. The obtained results demonstrate the robustness and effectiveness of the mixed AC/DC PF algorithm for power system operation and planning studies.

Author Contributions

F.M. was responsible for methodology, collecting resources, data analysis, writing—original draft preparation, and writing—review and editing. G.-A.N. and M.S. were responsible for the supervision, and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

M. Saif acknowledges the financial support of Natural Sciences and Engineering Research Council (NSERC) of Canada, as well as National Science Foundation of China under grant number 61873144.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Note: P and Q are in M W and M V A R , respectively.
Table A1. Results of AC PF calculations in Case 1. AC bus data.
Table A1. Results of AC PF calculations in Case 1. AC bus data.
BusVoltageGenerationLoad
| V |   ( p . u . ) θ   ( ° ) P   Q P   Q
11.0600.000131.1290.82
21.000−2.06140.00−61.5920.0010.00
30.987−4.63745.0015.00
40.984−4.95740.005.00
50.972−5.76560.0010.00
Total:171.1229.23165.0040.00
Table A2. Results of AC PF calculations in Case 1. AC branch data.
Table A2. Results of AC PF calculations in Case 1. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
11289.3374.00−86.85−72.912.4867.46
21341.7916.82−40.27−17.511.5184.55
32324.47−2.52−24.11−0.350.3601.08
42427.71−1.72−27.25−0.830.4611.38
52554.665.56−53.44−4.831.2153.65
63419.392.86−19.35−4.690.0400.12
7456.600.52−6.56−5.170.0430.13
Total:6.12318.37
Table A3. Results of the mixed AC/DC PF calculations in Case 2. AC bus data.
Table A3. Results of the mixed AC/DC PF calculations in Case 2. AC bus data.
BusVoltageGenerationLoad
| V |   ( p . u . ) θ   ( ° ) P   Q P   Q
11.0600.000133.6484.32
21.000−2.38340.00−32.8420.0010.00
31.000−3.89545.0015.00
40.996−4.26240.005.00
50.991−4.14960.0010.00
Total:173.6451.48165.0040.00
Table A4. Results of the mixed AC/DC PF calculations in Case 2. AC branch data.
Table A4. Results of the mixed AC/DC PF calculations in Case 2. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
11298.3871.37−95.66−69.592.7178.15
21335.2612.96−34.20−15.081.0623.19
32313.25−6.22−13.142.570.1160.35
42417.08−5.18−16.891.740.1810.54
52525.33−1.85−25.07−0.350.2570.77
63423.094.64−23.04−6.470.0570.17
745−0.07−0.270.07−4.650.0040.01
Total:4.39413.18
Table A5. Results of the mixed AC/DC PF calculations in Case 2. DC bus data.
Table A5. Results of the mixed AC/DC PF calculations in Case 2. DC bus data.
Bus DCBus ACVoltage Magnitude (p.u.)Active Power (MW)
121.008−58.627
231.00021.901
350.99836.186
Table A6. Results of the mixed AC/DC PF calculations in Case 2. Converter station data.
Table A6. Results of the mixed AC/DC PF calculations in Case 2. Converter station data.
Bus DCBus InjectionConverter VoltageTotal Loss
P Q | V |   ( p . u . ) θ   ( ° ) P  
1−60.00−40.000.890−13.0171.37
220.767.141.007−0.6551.14
335.005.000.9951.4421.19
Total:3.70
Bus DCConverter PowerFilterTransformer LossReactor LossConverter Loss
P Q Q P Q P Q P
1−59.22−32.63−8.120.085.830.019.661.29
220.76−0.65−9.020.010.540.000.701.14
335.02−0.37−8.830.021.430.002.031.17
Total:0.117.800.0112.393.60
Bus DCGrid PowerTrans. Filter PowerFilterConv. Filter PowerConv. Power
P Q P Q Q Q P Q
1−60.00−40.00−59.92−34.17−8.12−42.29−59.92−32.63
220.767.1420.767.68−9.02−1.3520.76−0.65
335.005.0035.026.43−8.83−2.4035.02−0.37
Table A7. Results of the mixed AC/DC PF calculations in Case 2. DC branch data.
Table A7. Results of the mixed AC/DC PF calculations in Case 2. DC branch data.
Branch From Bus To Bus Power Losses
From To P   P   P  
11230.66−30.420.24
2238.52−8.500.02
31327.96−27.680.28
Total:0.54
Table A8. Results of the mixed AC/DC PF calculations in Case 3. AC bus data.
Table A8. Results of the mixed AC/DC PF calculations in Case 3. AC bus data.
BusGenerationLoad
P   Q P   Q
1133.9384.93
240.00−90.4820.0010.00
345.0015.00
440.005.00
560.0010.00
Total:173.93−5.55165.0040.00
Table A9. Results of the mixed AC/DC PF calculations in Case 3. AC branch data.
Table A9. Results of the mixed AC/DC PF calculations in Case 3. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
11284.9575.29−82.56−74.502.3867.16
21348.989.64−47.17−9.501.8195.46
32334.58−12.23−33.8010.570.7802.34
42434.13−9.79−33.408.010.7352.21
52533.85−3.96−33.392.380.4611.38
634−1.6913.761.71−15.690.0220.07
745−8.312.688.39−7.380.0770.23
Total:6.28318.85
Table A10. Results of the mixed AC/DC PF calculations in Case 3. DC bus data.
Table A10. Results of the mixed AC/DC PF calculations in Case 3. DC bus data.
Bus DCBus InjectionTotal Loss
P Q P  
10.000.000.00
2−37.6529.841.22
335.005.001.19
Total:2.41
Table A11. Results of the mixed AC/DC PF calculations in Case 3. DC branch data.
Table A11. Results of the mixed AC/DC PF calculations in Case 3. DC branch data.
Branch From Bus To Bus Power Losses
From To P   P   P  
112−10.6710.700.03
22325.73−25.550.17
31310.67−10.630.04
Total:0.24
Table A12. Droop control settings for each converter station.
Table A12. Droop control settings for each converter station.
Converter StationDroop Parameter P D C   ( M W ) V D C   ( p . u . )
10.005−58.62741.0079
20.00721.90131.0000
30.00536.18560.9778
* P D C and V D C are the reference power and reference DC voltage, respectively.
Table A13. Results of the mixed AC/DC PF calculations in Case 4. AC bus data.
Table A13. Results of the mixed AC/DC PF calculations in Case 4. AC bus data.
BusGenerationLoad
P   Q P   Q
1133.3984.72
240.00−81.0520.0010.00
345.0015.00
440.005.00
560.0010.00
Total:173.390.67165.0040.00
Table A14. Results of the mixed AC/DC PF calculations in Case 4. AC branch data.
Table A14. Results of the mixed AC/DC PF calculations in Case 4. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
11290.3473.70−87.83−72.542.5107.53
21343.0511.02−41.59−11.961.4564.37
32325.47−9.78−25.057.060.4261.28
42428.39−7.76−27.885.290.5031.51
52553.97−0.97−52.811.521.1653.50
63418.0710.21−18.03−12.070.0450.14
7455.911.77−5.87−6.520.0430.13
Total:6.14818.46
Table A15. Results of the mixed AC/DC PF calculations in Case 4. DC bus data.
Table A15. Results of the mixed AC/DC PF calculations in Case 4. DC bus data.
Bus DCBus InjectionTotal Loss
P Q P  
10.000.000.00
2−3.5720.311.13
31.335.001.11
Total:2.24
Table A16. Results of the mixed AC/DC PF calculations in Case 4. DC branch data.
Table A16. Results of the mixed AC/DC PF calculations in Case 4. DC branch data.
Branch From Bus To Bus Power Losses
From To P   P   P  
112−0.720.720.00
2231.72−1.720.00
3130.72−0.720.00
Total:0.00
Table A17. Changes in the droop parameters of each converter station.
Table A17. Changes in the droop parameters of each converter station.
Converter StationDroop Parameter
10.0010
20.0014
30.0010
Table A18. Results of the mixed AC/DC PF calculations in Case 5. AC bus data.
Table A18. Results of the mixed AC/DC PF calculations in Case 5. AC bus data.
BusGenerationLoad
P   Q P   Q
1133.3884.72
240.00−81.2220.0010.00
345.0015.00
440.005.00
560.0010.00
Total:173.383.50165.0040.00
Table A19. Results of the mixed AC/DC PF calculations in Case 5. AC branch data.
Table A19. Results of the mixed AC/DC PF calculations in Case 5. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
11290.2473.73−87.74−72.582.5087.52
21343.1411.00−41.68−11.921.4614.38
32325.62−9.82−25.197.110.4301.29
42428.48−7.80−27.975.340.5071.52
52553.64−1.02−52.491.541.1513.45
63417.7510.26−17.71−12.120.0440.13
7455.681.78−5.64−6.540.0410.12
Total:6.14218.41
Table A20. Results of the mixed AC/DC PF calculations in Case 5. DC bus data.
Table A20. Results of the mixed AC/DC PF calculations in Case 5. DC bus data.
Bus DCBus InjectionTotal Loss
P Q P  
10.000.000.00
2−4.1120.451.13
31.875.001.11
Total:2.24
Table A21. Results of the mixed AC/DC PF calculations in Case 5. DC branch data.
Table A21. Results of the mixed AC/DC PF calculations in Case 5. DC branch data.
Branch From Bus To Bus Power Losses
From To P   P   P  
112−0.880.880.00
2232.11−2.110.00
3130.88−0.880.00
Total:0.00
Table A22. Results of the mixed AC/DC PF calculations in Case 6-1. AC bus data.
Table A22. Results of the mixed AC/DC PF calculations in Case 6-1. AC bus data.
BusGenerationLoad
P   Q P   Q
1133.3884.72
240.00−81.2220.0010.00
345.0015.00
440.005.00
560.0010.00
Total:173.383.50165.0040.00
Table A23. Results of the mixed AC/DC PF calculations in Case 6-1. AC branch data.
Table A23. Results of the mixed AC/DC PF calculations in Case 6-1. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
11290.2473.73−87.74−72.582.5087.52
21343.1411.00−41.68−11.921.4614.38
32325.62−9.82−25.197.110.4301.29
42428.48−7.80−27.975.340.5071.52
52553.64−1.02−52.491.541.1513.45
63417.7510.26−17.71−12.120.0440.13
7455.681.78−5.64−6.540.0410.12
Total:6.14218.41
Table A24. Results of the mixed AC/DC PF calculations in Case 6-1. DC bus data.
Table A24. Results of the mixed AC/DC PF calculations in Case 6-1. DC bus data.
Bus DCBus InjectionTotal Loss
P Q P  
10.000.000.00
2−4.1120.451.13
31.875.001.11
Total:2.24
Table A25. Results of the mixed AC/DC PF calculations in Case 6-1. DC branch data.
Table A25. Results of the mixed AC/DC PF calculations in Case 6-1. DC branch data.
Branch From Bus To Bus Power Losses
From To P   P   P  
112−0.880.880.00
2232.11−2.110.00
3130.88−0.880.00
Total:0.00
Table A26. Results of the mixed AC/DC PF calculations in Case 6-2. AC bus data.
Table A26. Results of the mixed AC/DC PF calculations in Case 6-2. AC bus data.
BusGenerationLoad
P   Q P   Q
1135.6383.69
240.00−61.3820.0010.00
345.0015.00
440.005.00
560.0010.00
Total:175.6322.31165.0040.00
Table A27. Results of the mixed AC/DC PF calculations in Case 6-2. AC branch data.
Table A27. Results of the mixed AC/DC PF calculations in Case 6-2. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
112113.5967.07−110.41−63.903.1809.54
21322.0416.62−21.42−20.080.6151.84
323−7.870.697.91−4.570.0420.12
4240.180.16−0.18−4.130.0030.01
52516.890.63−16.78−3.250.1160.35
63448.34−3.45−48.102.160.2340.70
7458.28−3.02−8.22−1.750.0560.17
Total:4.24612.73
Table A28. Results of the mixed AC/DC PF calculations in Case 6-2. DC bus data.
Table A28. Results of the mixed AC/DC PF calculations in Case 6-2. DC bus data.
Bus DCBus InjectionTotal Loss
P Q P  
1−121.218.961.69
279.82−13.101.38
335.005.001.19
Total:4.26
Table A29. Results of the mixed AC/DC PF calculations in Case 6-2. DC branch data.
Table A29. Results of the mixed AC/DC PF calculations in Case 6-2. DC branch data.
Branch From Bus To Bus Power Losses
From To P   P   P  
11273.60−72.241.36
223−8.978.990.02
31345.92−45.180.74
Total:2.12
Table A30. Parameters of the converter stations in per-unit.
Table A30. Parameters of the converter stations in per-unit.
Converter Station S N   ( M V A ) R T R X T R B F R T R X T R
11000.00150.11210.08870.00010.1643
21000.00150.11210.08870.00010.1643
31000.00150.11210.08870.00010.1643
41000.00150.11210.08870.00010.1643
51000.00150.11210.08870.00010.1643
61000.00150.11210.08870.00010.1643
71000.00150.11210.08870.00010.1643
Table A31. Power losses coefficients of the converter stations.
Table A31. Power losses coefficients of the converter stations.
Converter Station a b c r e c c i n v
11.1030.8872.8854.371
21.1030.8872.8854.371
31.1030.8872.8854.371
41.1030.8872.8854.371
51.1030.8872.8854.371
61.1030.8872.8854.371
71.1030.8872.8854.371
Table A32. Parameters of the DC lines.
Table A32. Parameters of the DC lines.
FromTo R D C   ( p . u . ) V D C   ( k V )   P D C   ( M W )
120.0520345100
140.0520345100
230.0520345100
340.0520345100
560.0730345100
570.0730345100
580.0730345100
670.0730345100
680.0730345100
Table A33. Results of AC PF calculations in Case 7-1. AC bus data.
Table A33. Results of AC PF calculations in Case 7-1. AC bus data.
BusVoltageGenerationLoad
| V |   ( p . u . ) θ   ( ° ) P   Q P   Q
11.039−13.53797.6044.20
21.048−9.785
31.031−12.276322.002.40
41.004−12.627500.00184.00
51.006−11.192
61.008−10.408
70.998−12.756233.8084.00
80.998−13.336522.00176.60
91.038−14.1786.50−66.60
101.018−8.171
111.013−8.937
121.001−8.9998.5388.00
131.015−8.930
141.012−10.715
151.016−11.345320.00153.00
161.033−10.033329.0032.30
171.034−11.116
181.032−11.986158.0030.00
191.050−5.410
200.991−6.821680.00103.00
211.032−7.629274.00115.00
221.050−3.183
231.045−3.381247.5084.60
241.038−9.914308.60−92.20
251.058−8.369224.0047.20
261.053−9.439139.0017.00
271.038−11.362281.0075.50
281.050−5.928206.0027.60
291.050−3.170283.5026.90
301.050−7.370250.00161.76
310.9820.000677.87221.579.204.60
320.984−0.188650.00206.96
330.997−0.196632.00108.29
341.012−1.631508.00166.69
351.0491.777650.00210.66
361.0644.468560.00100.16
371.028−1.583540.00−1.37
381.0273.893830.0021.73
391.030−14.5351000.078.471104.00250.00
Total:6297.871274.926254.231387.10
Table A34. Results of AC PF calculations in Case 7-1. AC branch data.
Table A34. Results of AC PF calculations in Case 7-1. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
112−173.7−40.31174.68−24.360.97811.48
213976.10−3.89−76.03−74.750.0661.65
323319.9188.59−318.58−100.881.33515.51
4225−244.5982.97248.93−93.844.3375.33
5230−250.00−147.20250.00161.760.00014.56
63437.34113.06−37.13−132.590.2083.40
7318−40.76−14.5940.78−7.940.0170.21
845−197.45−4.09197.76−4.520.3094.95
9414−265.42−47.32265.9942.480.5719.22
1056−536.94−43.11537.5146.160.5737.45
1158339.1847.64−338.24−49.390.93313.07
1267453.8281.55−452.56−73.591.26119.33
13611−322.65−38.85323.3833.140.7248.48
14631−668.67−88.85668.67216.970.000128.12
1578218.76−10.41−218.564.840.1922.21
168934.81−132.06−34.4897.720.3245.11
1793927.98−31.12−27.97−96.780.0180.44
181011327.9073.37−327.46−76.180.4384.71
191013322.1037.49−321.69−40.650.4074.38
201032−650.00−110.87650.00206.960.00096.1
211211−4.06−42.254.0943.040.0290.79
221213−4.47−45.754.5146.680.0340.93
231314317.18−6.03−316.30−1.800.8799.87
24141550.31−40.68−50.263.660.0530.64
251516−269.74−156.66270.56147.330.8258.61
261617224.02−42.54−223.6832.500.3384.29
271619−451.30−54.20454.3858.753.07837.52
281621−329.6014.44330.42−27.740.82113.86
291624−42.68−97.3342.7190.630.0300.59
301718199.0411.05−198.78−22.060.2613.06
31172724.64−43.56−24.629.230.0160.21
321920174.73−9.170−174.5113.480.2184.30
331933−629.11−49.58632.00108.292.89458.71
342034−505.49−116.48508.00166.692.51150.21
352122−604.42−87.26607.21108.152.78348.7
36222342.79041.88−42.77−61.750.0250.40
372235−650.00−150.04650.00210.660.00060.63
382324353.84−0.50−351.311.572.52940.24
392336−558.57−22.35560.00100.161.43077.82
40252665.41−18.81−65.29−39.040.1261.27
412537−538.3465.45540.00−1.371.65764.08
422627257.3068.21−256.38−84.730.9209.66
432628−140.82−21.21141.61−56.360.7888.69
442629−190.19−24.96192.10−67.791.91420.98
452829−347.6128.76349.16−39.441.55616.78
462938−824.7780.33830.021.735.234102.06
Total:43.6401000.61
Table A35. Results of Mixed AC/DC PF calculations in Case 7-2. AC bus data.
Table A35. Results of Mixed AC/DC PF calculations in Case 7-2. AC bus data.
BusVoltageGenerationLoad
| V |   ( p . u . ) θ   ( ° ) P   Q P   Q
11.047−12.43997.6044.20
21.051−8.959
31.034−11.548322.002.400
40.998−12.879500.00184.00
51.001−11.443
61.004−10.675
70.994−12.969233.8084.00
80.994−13.517522.00176.60
91.037−13.8046.50−66.60
101.012−8.658
111.009−9.390
121.000−10.0758.53088.00
131.008−9.464
141.000−11.243
151.016−10.879320.00153.00
161.035−9.343329.0032.30
171.040−10.213
181.038−11.041158.0030.00
191.051−4.727
200.991−6.135680.00103.00
211.034−6.945274.00115.00
221.051−2.507
231.046−2.705247.5084.60
241.040−9.223308.60−92.20
251.060−7.519224.0047.20
261.056−8.561139.0017.00
271.043−10.465281.0075.50
281.052−5.060206.0027.60
291.051−2.306283.5026.90
301.050−6.550250.00147.24
310.9820.000691.88240.399.204.600
320.984−0.633650.00232.33
330.9970.488632.00103.03
341.012−0.946508.00164.25
351.0492.449650.00204.75
361.0645.139560.0096.81
371.028−0.744540.00−10.62
381.0264.752830.0014.76
391.030−13.7881000.054.671104.00250.00
Total:6311.881247.616254.231387.10
Table A36. Results of Mixed AC/DC PF calculations in Case 7-2. AC branch data.
Table A36. Results of Mixed AC/DC PF calculations in Case 7-2. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
112−161.82−30.43162.66−36.590.8389.85
2139104.2226.23−104.08−103.60.1403.51
323333.4580.70−332.03−92.241.41416.43
4225−246.1189.31250.55−100.124.4395.45
5230−250.00−133.42250.00147.240.00013.82
634123.06160.38−122.51−174.320.5448.92
7318−73.03−30.5473.088.300.0590.71
845−196.54−21.20196.8612.780.3125.00
9414−220.94−8.49221.341.0300.3926.33
1056−521.95−55.55522.5058.320.5497.14
1158325.0942.76−324.23−45.370.86312.09
1267439.6879.60−438.49−72.561.19418.31
13611−279.50−37.36280.0529.730.5496.43
14631−682.68−100.56682.68235.790.000135.24
1578204.69−11.44−204.525.680.1701.95
16896.75−136.92−6.420102.830.3265.15
17939−0.08−36.230.08−91.730.0070.19
181011309.1859.26−308.79−62.530.3894.18
191013340.8274.67−340.34−76.980.4775.13
201032−650.00−133.93650.00232.330.00098.4
211211−28.72−31.9828.7532.800.0300.81
221213−25.68−30.3925.7131.080.0260.70
231314314.6345.89−313.73−53.120.90410.14
241415−35.63−89.1635.7553.330.1131.36
251516−315.75−176.33316.86169.981.11411.63
261617177.78−80.31−177.5468.940.2423.07
271619−451.33−47.15454.3951.363.05937.28
281621−329.6219.33330.44−32.750.82013.84
291624−42.68−94.1642.7187.390.0280.55
301718191.326.86−191.08−18.300.2382.79
31172726.22−35.79−26.201.0900.0120.16
321920174.72−6.89−174.511.180.2184.29
331933−629.11−44.47632.00103.032.88658.55
342034−505.50−114.18508.00164.252.50350.07
352122−604.44−82.25607.21102.852.77048.47
36222342.7941.60−42.76−61.500.0250.39
372235−650.00−144.44650.00204.750.00060.31
382324353.84−3.95−351.314.812.52240.13
392336−558.57−19.16560.0096.811.42877.66
40252663.79−21.80−63.68−36.420.1181.19
412537−538.3474.72540.00−10.621.65864.10
422627255.6859.54−254.8−76.590.8889.32
432628−140.79−18.17141.58−59.770.7908.70
442629−190.22−21.95192.14−71.241.91521.00
452829−347.5832.17349.13−42.931.55516.77
462938−824.7787.27830.0014.765.232102.02
Total:43.7561009.53
Table A37. Results of the mixed AC/DC PF calculations in Case 7-2. DC bus data.
Table A37. Results of the mixed AC/DC PF calculations in Case 7-2. DC bus data.
Bus DCBus InjectionTotal Loss
P Q P  
1−40.00−20.001.24
2−45.8725.631.25
340.0040.001.24
440.0040.001.24
5−128.03−141.242.44
640.0030.001.22
740.0020.001.20
840.0040.001.24
Total:11.07
Table A38. Results of the mixed AC/DC PF calculations in Case 7-2. DC branch data.
Table A38. Results of the mixed AC/DC PF calculations in Case 7-2. DC branch data.
Branch From Bus To Bus Power Losses
From To P   P   P  
112−2.192.190.00
21440.95−40.510.44
32342.43−41.960.47
4340.72−0.720.00
55641.86−41.220.64
65741.85−41.210.64
75841.87−41.230.64
867−0.010.010.00
9680.01−0.010.00
Total:2.83
Table A39. Droop control settings for each converter station.
Table A39. Droop control settings for each converter station.
Converter StationDroop Parameter P D C   ( M W ) V D C   ( p . u . )
10.0017−50.00001.0079
20.0017−35.00001.0000
30.001760.00001.0000
40.001760.00001.0000
50.0017−35.00001.0000
60.001740.00001.0000
70.001740.00001.0000
80.001740.00001.0000
* P D C and V D C are the reference power and reference DC voltage, respectively.
Table A40. Results of Mixed AC/DC PF calculations in Case 7-3. AC bus data.
Table A40. Results of Mixed AC/DC PF calculations in Case 7-3. AC bus data.
BusVoltageGenerationLoad
| V |   ( p . u . ) θ   ( ° ) P   Q P   Q
11.047−12.50497.6044.20
21.051−9.194
31.034−11.802322.002.40
40.997−12.917500.00184.00
51.001−11.426
61.004−10.647
70.994−12.949233.8084.00
80.994−13.501522.00176.60
91.037−13.8236.50−66.60
101.012−8.558
111.009−9.310
121.000−9.9268.53088.00
131.008−9.343
141.000−11.08
151.016−11.186320.00153.00
161.034−9.727329.0032.30
171.04−10.668
181.037−11.483158.0030.00
191.051−5.110
200.991−6.519680.00103.00
211.034−7.328274.00115.00
221.051−2.889
231.046−3.087247.5084.60
241.040−9.607308.60−92.20
251.060−7.783224.0047.20
261.055−8.920139.0017.00
271.042−10.869281.0075.50
281.052−5.418206.0027.60
291.051−2.664283.5026.90
301.050−6.785250.00146.82
310.9820.000689.99240.449.2004.60
320.984−0.532650.00232.48
330.9970.105632.00103.62
341.012−1.329508.00164.52
351.0492.067650.00205.41
361.0644.758560.0097.19
371.028−1.007540.00−9.75
381.0274.395830.0015.39
391.030−13.8301000.0054.041104.00250.00
Total:6309.991250.166254.231387.10
Table A41. Results of Mixed AC/DC PF calculations in Case 7-3. AC branch data.
Table A41. Results of Mixed AC/DC PF calculations in Case 7-3. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
112−154.00−31.28154.76−36.70.7598.91
2139102.5227.08−102.38−104.510.1383.46
323335.8681.69−334.43−92.981.43516.67
4225−240.6288.03244.87−99.074.2535.23
5230−250.00−133.01250.00146.820.00013.8
634104.84161.43−104.34−176.120.4988.17
7318−46.55−30.8546.578.230.0260.32
845−203.98−20.74204.3212.720.3365.38
9414−247.94−7.14248.431.330.4947.97
1056−529.57−55.39530.1358.380.5657.35
1158325.2542.68−324.39−45.260.86412.1
1267441.2479.61−440.04−72.451.20318.44
13611−290.58−36.72291.1729.610.5936.94
14631−680.79−101.27680.79235.840.000134.58
1578206.24−11.55−206.075.830.1721.98
16898.45−137.17−8.12103.120.3285.18
179391.62−36.52−1.62−91.440.0070.18
181011317.5759.28−317.16−62.330.4094.4
191013332.4374.78−331.97−77.320.4554.9
201032−650.00−134.06650.00232.480.00098.41
211211−25.96−31.9725.9932.720.0270.75
221213−24.57−30.4224.5931.090.0250.67
231314307.3846.23−306.52−53.900.8649.7
2414152.66−91.24−2.5655.210.0961.16
251516−300.93−178.21301.97171.11.04010.86
261617192.66−79.73−192.3968.820.2773.52
271619−451.33−47.94454.3952.193.06137.3
281621−329.6218.78330.44−32.190.82013.84
291624−42.68−94.5242.7187.760.0280.56
301718188.286.71−188.04−18.230.2312.7
31172720.60−35.53−20.590.800.0090.12
321920174.72−7.14−174.5011.440.2184.29
331933−629.11−45.05632.00103.622.88758.57
342034−505.50−114.44508.00164.522.50450.08
352122−604.44−82.81607.21103.442.77148.5
36222342.7941.63−42.76−61.530.0250.39
372235−650.00−145.07650.00205.410.00060.34
382324353.84−3.56−351.314.442.52340.14
392336−558.57−19.51560.0097.191.42877.67
40252669.47−21.98−69.33−35.990.1391.41
412537−538.3473.85540.00−9.751.65864.1
422627261.3459.66−260.41−76.300.9269.72
432628−140.79−18.45141.58−59.460.7898.70
442629−190.22−22.23192.13−70.921.91521.0
452829−347.5831.86349.13−42.611.55516.77
462938−824.7786.63830.0015.395.232102.03
Total:43.5831009.26
Table A42. Results of the mixed AC/DC PF calculations in Case 7-3. DC bus data.
Table A42. Results of the mixed AC/DC PF calculations in Case 7-3. DC bus data.
Bus DCBus InjectionTotal Loss
P Q P  
1−56.26−20.001.30
2−42.0025.611.23
345.8640.001.26
446.1240.001.26
5−55.43−143.802.04
616.5130.001.16
716.5320.001.14
816.4840.001.19
Total:10.58
Table A43. Results of the mixed AC/DC PF calculations in Case 7-3. DC branch data.
Table A43. Results of the mixed AC/DC PF calculations in Case 7-3. DC branch data.
Branch From Bus To Bus Power Losses
From To P   P   P  
1125.27−5.260.01
21449.69−49.030.65
32346.03−45.470.56
434−1.651.650.00
55617.79−17.670.12
65717.79−17.670.12
75817.79−17.670.12
8670.000.000.00
9680.000.000.00
Total:1.58
Table A44. Results of Mixed AC/DC PF calculations in Case 7-4. AC bus data.
Table A44. Results of Mixed AC/DC PF calculations in Case 7-4. AC bus data.
BusVoltageGenerationLoad
| V |   ( p . u . ) θ   ( ° ) P   Q P   Q
11.043−15.38397.6044.20
21.058−11.984
31.055−14.490322.002.40
40.932−16.184500.00184.00
50.963−12.185
60.971−10.906
70.961−13.647233.8084.00
80.960−14.378522.00176.60
91.023−15.6246.50−66.6
100.995−6.307
110.985−7.818
120.975−7.1658.5388.00
130.993−6.380
140.994−6.385
151.016−14.260320.00153.00
161.033−12.688329.0032.30
171.035−13.530
181.032−14.217158.0030.00
191.050−8.066
200.991−9.477680.00103.00
211.032−10.284274.00115.00
221.050−5.839
231.045−6.037247.5084.60
241.038−12.569308.60−92.20
251.064−10.503224.0047.20
261.055−11.705139.0017.00
271.040−13.691281.0075.50
281.052−8.202206.0027.60
291.051−5.447283.5026.90
301.050−9.591250.00108.31
310.9820.000683.41362.789.204.60
320.9841.860650.00313.16
330.997−2.849632.00107.87
341.012−4.287508.00166.49
351.049−0.880650.00210.19
361.0641.812560.0099.90
371.028−3.746540.00−26.63
381.0271.612830.0016.42
391.030−16.1971000.00126.311104.00250.00
Total:6303.411484.806254.231387.10
Table A45. Results of Mixed AC/DC PF calculations in Case 7-4. AC branch data.
Table A45. Results of Mixed AC/DC PF calculations in Case 7-4. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
112−160.79−56.77161.63−10.460.8439.90
213963.1912.57−63.13−91.600.0631.57
323323.21−12.20−322.00−2.401.21314.09
4225−234.85118.79239.31−129.744.4585.48
5230−250.00−96.12250.00108.310.00012.19
6340.000.000.000.000.0000.00
73180.000.000.000.000.0000.00
845−500.00−184.00502.60213.482.59541.53
94140.000.000.000.000.0000.00
1056−819.04−214.86820.59230.881.54420.08
1158316.451.38−315.58−2.930.86412.10
1267490.2481.71−488.66−68.051.57924.20
13611−636.61−105.52639.69128.363.08436.12
14631−674.21−207.08674.21358.180.000151.11
1578254.86−15.95−254.5812.010.2823.25
168948.16−185.68−47.40160.290.76312.05
1793940.90−93.69−40.87−32.100.0250.63
181011617.69175.82−616.01−165.001.67217.97
19101332.3129.83−32.31−36.940.0090.09
201032−650.00−205.65650.00313.160.000107.51
21121123.71−35.79−23.6836.640.0310.85
221213−32.24−52.2132.3153.950.0641.74
2313140.00−17.010.000.000.0010.01
2414150.000.000.000.000.0000.00
251516−320.00−153.00321.07146.271.07411.21
261617173.51−42.68−173.3130.960.2062.62
271619−451.30−53.64454.3858.173.07737.50
281621−329.6014.83330.42−28.140.82113.85
291624−42.68−97.0842.7190.370.0300.58
301718158.1717.88−158.00−30.000.1681.96
31172715.14−48.84−15.1214.440.0150.20
321920174.73−8.99−174.5113.300.2184.30
331933−629.11−49.17632.00107.872.89458.70
342034−505.49−116.30508.00166.492.51050.20
352122−604.42−86.86607.21107.732.78248.69
36222342.7941.86−42.77−61.730.0250.40
372235−650.00−149.59650.00210.190.00060.60
382324353.84−0.77−351.311.832.52940.23
392336−558.57−22.10560.0099.901.43077.80
40252675.03−8.33−74.86−49.500.1731.74
412537−538.3490.86540.00−26.631.66164.23
422627266.8774.07−265.88−89.940.99210.42
432628−140.79−18.90141.58−58.950.7898.70
442629−190.21−22.67192.13−70.411.91520.99
452829−347.5831.35349.14−42.101.55516.78
462938−824.7785.61830.0016.425.232102.03
Total:49.1861106.20
Table A46. Results of Mixed AC/DC PF calculations in Case 7-5. AC bus data.
Table A46. Results of Mixed AC/DC PF calculations in Case 7-5. AC bus data.
BusVoltageGenerationLoad
| V |   ( p . u . ) θ   ( ° ) P   Q P   Q
11.046−5.46197.6044.20
21.0601.265
31.063−0.937322.002.40
40.934−16.570500.00184.00
50.969−12.289
60.977−11.124
70.967−13.306233.8084.00
80.966−13.764522.00176.60
91.025−10.9116.50−66.60
101.003−7.787
110.994−8.935
121.000−9.2468.5388.00
131.002−8.201
141.000−8.935
151.0324.023320.00153.00
161.0455.368329.0032.30
171.0514.342
181.0503.841158.0030.00
191.0559.953
200.9938.557680.00103.00
211.0417.738274.00115.00
221.05512.140
231.05011.944247.5084.60
241.0495.488308.60−92.20
251.0703.371224.0047.20
261.0644.157139.0017.00
271.0533.112281.0075.50
281.0577.633206.0027.60
291.05410.373283.5026.90
301.0503.655250.0099.19
310.9820.000701.04343.289.204.60
320.9840.317650.00277.42
330.99715.160632.0078.96
341.01213.742508.00153.10
351.04917.078650.00177.74
361.06419.762560.0081.49
371.02810.095540.00−55.46
381.02617.419830.00−5.42
391.030−8.8491000.00115.541104.00250.00
Total:6321.041265.846254.231387.10
Table A47. Results of Mixed AC/DC PF calculations in Case 7-5. AC branch data.
Table A47. Results of Mixed AC/DC PF calculations in Case 7-5. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
112−314.89−27.33318.06−12.823.17637.29
2139257.2923.13−256.65−87.880.64316.07
323282.95−55.56−282.0037.600.94711.00
4225−351.01155.69360.37−160.759.35711.50
5230−250.00−87.31250.0099.190.00011.88
6340.000.000.000.000.0000.00
73180.000.000.000.000.0000.00
845−540.00−204.00543.03240.393.03448.54
94140.000.000.000.000.0000.000
1056−758.72−243.95760.07257.421.35217.58
1158215.683.56−215.29−11.800.3985.57
1267396.7885.56−395.74−80.251.04215.98
13611−465.01−158.13466.76165.201.75520.55
14631−691.84−184.85691.84338.680.000153.82
1578161.94−3.75−161.82−2.250.1121.29
1689−144.89−162.55145.92141.151.03416.33
17939−152.42−74.55152.65−46.590.2225.56
181011480.08167.66−479.05−163.811.03411.11
191013169.926.61−169.8−12.690.1151.24
201032−650.00−174.27650.00277.420.000103.15
211211−12.281.4612.29−1.390.0020.07
221213−42.29−15.3142.3216.200.0330.89
231314127.48−3.51−127.33−12.110.1461.64
2414150.000.000.000.000.0000.00
251516−280.00−123.00280.77112.640.7738.07
261617214.10−92.92−213.7682.520.3414.34
271619−451.47−14.50454.4617.422.98836.42
281621−329.7141.95330.53−55.830.82013.83
291624−42.70−79.4842.7272.440.0210.41
301718118.09−3.51−118.00−10.000.0891.04
311727135.67−39.01−135.456.390.2222.96
321920174.693.61−174.470.650.2164.26
331933−629.14−21.03632.0078.962.85657.93
342034−505.53−103.65508.00153.102.47249.45
352122−604.53−59.17607.2478.502.71447.49
36222342.7640.27−42.73−60.330.0240.38
372235−650.00−118.77650.00177.740.00058.97
382324353.82−19.77−351.3219.762.49939.75
392336−558.58−4.49560.0081.491.41577.00
402526−46.04−6.6646.12−53.070.0750.76
412537−538.33120.21540.00−55.461.67564.75
422627145.8858.48−145.55−81.890.3273.43
432628−140.68−9.26141.48−69.690.7978.79
442629−190.31−13.14192.24−81.251.92521.10
452829−347.4842.09349.03−53.071.55416.76
462938−824.77107.42830.00−5.425.231102.00
Total:53.4361110.95
Table A48. Results of the mixed AC/DC PF calculations in Case 7-5. DC bus data.
Table A48. Results of the mixed AC/DC PF calculations in Case 7-5. DC bus data.
Bus DCBus InjectionTotal Loss
P Q P  
1−40.00−20.001.25
2−46.0474.151.42
340.0040.001.23
440.0040.001.24
5−127.33−12.111.75
640.0030.001.22
740.0020.001.2
840.0040.001.24
Total:10.55
Table A49. Results of the mixed AC/DC PF calculations in Case 7-5. DC branch data.
Table A49. Results of the mixed AC/DC PF calculations in Case 7-5. DC branch data.
Branch From Bus To Bus Power Losses
From To P   P   P  
112−2.202.200.00
21440.95−40.510.44
32342.43−41.960.47
4340.73−0.730.00
55641.86−41.220.64
65741.85−41.210.64
75841.87−41.230.64
867−0.010.010.00
9680.01−0.010.00
Total:2.83
Table A50. Results of Mixed AC/DC PF calculations in Case 7-6. AC bus data.
Table A50. Results of Mixed AC/DC PF calculations in Case 7-6. AC bus data.
BusVoltageGenerationLoad
| V |   ( p . u . ) θ   ( ° ) P   Q P   Q
11.049−8.26797.6044.20
21.062−2.889
31.066−5.034322.002.40
40.932−16.723500.00184.00
50.967−12.294
60.975−11.059
70.965−13.409233.8084.00
80.965−13.950522.00176.60
91.025−12.2406.50−66.60
101.002−7.176
110.993−8.494
121.000−8.5898.5388.00
131.001−7.424
141.000−7.736
151.03−3.037320.00153.00
161.044−1.572329.0032.30
171.050−2.488
181.048−3.091158.0030.00
191.0543.016
200.9931.618680.00103.00
211.0400.800274.00115.00
221.0545.205
231.055.008247.5084.60
241.048−1.452308.60−92.20
251.069−1.128224.0047.20
261.063−1.501139.0017.00
271.052−3.084281.0075.50
281.0561.978206.0027.60
291.0544.719283.5026.90
301.050−0.506250.0082.91
310.9820.000696.02347.709.204.60
320.9840.934650.00280.53
330.9978.223632.0080.97
341.0126.803508.00154.03
351.04910.145650.00180.00
361.06412.829560.0082.78
371.0285.600540.00−52.18
381.02711.767830.00−3.06
391.030−10.9161000.0099.461104.00250.00
Total:6316.021253.146254.231387.10
Table A51. Results of Mixed AC/DC PF calculations in Case 7-6. AC branch data.
Table A51. Results of Mixed AC/DC PF calculations in Case 7-6. AC branch data.
Branch From Bus InjectionTo Bus InjectionPower Losses
From To P   Q P   Q P   Q
112−254.22−39.12256.28−14.612.05624.14
2139202.7534.92−202.32−105.310.42610.66
323277.03−56.22−276.1337.600.90410.50
4225−283.31142.35289.70−151.106.3847.84
5230−250.00−71.52250.0082.910.00011.39
6340.000.000.000.000.0000.00
73180.000.000.000.000.0000.00
845−556.28−204.00559.49243.273.21151.37
94140.000.000.000.000.0000.00
1056−800.66−243.89802.15259.241.49619.45
1158241.170.62−240.67−7.410.4986.97
1267425.5484.41−424.34−76.761.19318.29
13611−540.87−153.36543.18166.982.31127.07
14631−686.82−190.29686.82343.100.000152.81
1578190.54−7.24−190.391.780.1561.79
1689−90.94−170.9691.73145.680.78512.39
17939−98.23−79.0898.32−45.230.0942.36
181011548.16170.50−546.84−163.591.31814.17
191013101.846.52−101.80−13.390.0420.45
201032−650.00−177.03650.00280.530.000103.50
211211−3.663.413.66−3.400.0000.01
221213−47.06−14.1047.1015.160.0391.06
23131454.70−1.77−54.67−15.170.0270.31
2414150.000.000.000.000.0000.00
251516−303.49−123.00304.38113.910.8919.30
261617190.48−88.27−190.2077.070.2753.50
271619−451.46−17.25454.4520.272.99336.48
281621−329.7040.04330.52−53.880.81913.83
291624−42.70−80.7242.7273.700.0210.42
301718141.60−3.02−141.47−10.000.1281.49
31172765.09−34.05−65.04−0.740.0530.71
321920174.692.73−174.471.530.2164.26
331933−629.14−23.00632.0080.972.85857.97
342034−505.53−104.53508.00154.032.47549.50
352122−604.52−61.12607.2480.552.71847.57
36222342.7640.38−42.74−60.430.0240.38
372235−650.00−120.93650.00180.000.00059.07
382324353.82−18.44−351.3218.502.50039.78
392336−558.58−5.73560.0082.781.41677.05
40252624.63−12.96−24.61−47.170.0250.26
412537−538.33116.86540.00−52.181.67364.68
422627216.6054.67−215.96−74.760.6386.70
432628−140.69−10.32141.49−68.530.7968.78
442629−190.30−14.18192.23−80.071.92321.09
452829−347.4940.93349.04−51.881.55416.76
462938−824.77105.05830.00−3.065.230101.99
Total:50.1661098.10
Table A52. Results of the mixed AC/DC PF calculations in Case 7-6. DC bus data.
Table A52. Results of the mixed AC/DC PF calculations in Case 7-6. DC bus data.
Bus DCBus InjectionTotal Loss
P Q P  
1−56.28−20.001.32
2−42.1977.311.42
345.8740.001.25
446.1240.001.26
5−54.67−15.171.29
616.5130.001.16
716.5320.001.14
816.4940.001.19
Total:10.03
Table A53. Results of the mixed AC/DC PF calculations in Case 7-6. DC branch data.
Table A53. Results of the mixed AC/DC PF calculations in Case 7-6. DC branch data.
Branch From Bus To Bus Power Losses
From To P   P   P  
1125.27−5.260.01
21449.69−49.030.65
32346.03−45.470.56
434−1.651.650.00
55617.79−17.670.12
65717.79−17.670.12
75817.79−17.670.12
8670.000.000.00
9680.000.000.00
Total:1.58

References

  1. Mohammadi, F.; Nazri, G.-A.; Saif, M. A Bidirectional Power Charging Control Strategy for Plug-in Hybrid Electric Vehicles. Sustainability 2019, 11, 4317. [Google Scholar] [CrossRef] [Green Version]
  2. Bahrman, M.P.; Johnson, B.K. The ABCs of HVDC Transmission Technologies. IEEE Power Energy Mag. 2007, 5, 32–44. [Google Scholar] [CrossRef]
  3. Flourentzou, N.; Agelidis, V.G.; Demetriades, G.D. VSC-Based HVDC Power Transmission Systems: An Overview. IEEE Trans. Power Electron. 2009, 24, 592–602. [Google Scholar] [CrossRef]
  4. Mohammadi, F.; Nazri, G.-A.; Saif, M. An Improved Droop-Based Control Strategy for MT-HVDC Systems. Electronics 2020, 9, 87. [Google Scholar]
  5. Xu, L.; Fan, L. Impedance-Based Resonance Analysis in a VSC-HVDC System. IEEE Trans. Power Deliv. 2013, 28, 2209–2216. [Google Scholar]
  6. Xu, L.; Fan, L.; Miao, Z. DC Impedance-Model-Based Resonance Analysis of a VSC-HVDC System. IEEE Trans. Power Deliv. 2015, 30, 1221–1230. [Google Scholar] [CrossRef]
  7. Pandya, K.S.; Joshi, S.K. A Survey of Optimal Power Flow Methods. J. Theor. Appl. Inf. Technol. 2008, 4, 450–458. [Google Scholar]
  8. Frank, S.; Steponavice, I.; Rebennack, S. Optimal Power Flow: A Bibliographic Survey II: Non-Deterministic and Hybrid Methods. Energy Syst. 2012, 3, 259–289. [Google Scholar] [CrossRef]
  9. Mohammadi, F.; Nazri, G.-A.; Saif, M. A New Topology of a Fast Proactive Hybrid DC Circuit Breaker for MT-HVDC Grids. Sustainability 2019, 11, 4493. [Google Scholar] [CrossRef] [Green Version]
  10. Mohammadi, F.; Zheng, C. Stability Analysis of Electric Power System. In Proceedings of the 4th National Conference on Technology in Electrical and Computer Engineering, Tehran, Iran, 27 December 2018. [Google Scholar]
  11. Mohammadi, F.; Nazri, G.-A.; Saif, M. A Fast Fault Detection and Identification Approach in Power Distribution Systems. In Proceedings of the 5th International Conference on Power Generation Systems and Renewable Energy Technologies (PGSRET), Istanbul, Turkey, 26–27 August 2019. [Google Scholar]
  12. Arrillaga, J.; Bodger, P. Integration of HVDC Links with Fast-Decoupled Load-Flow Solutions. Proc. Inst. Electr. Eng. 1977, 124, 463–468. [Google Scholar] [CrossRef]
  13. Baradar, M.; Ghandhari, M.; van Hertem, D. The Modeling Multi-Terminal VSC-HVDC in Power Flow Calculation using Unified Methodology. In Proceedings of the 2nd IEEE PES International Conference and Exhibition on Innovative Smart Grid Technologies, Manchester, UK, 5–7 December 2011. [Google Scholar]
  14. Baradar, M.; Ghandhari, M. A Multi-Option Unified Power-Flow Approach for Hybrid AC/DC Grids Incorporating Multi-terminal VSC-HVDC. IEEE Trans. Power Syst. 2013, 28, 2376–2383. [Google Scholar] [CrossRef]
  15. El-Hawary, M.E.; Ibrahim, S.T. A New Approach to AC-DC Load Flow Analysis. Electr. Power Syst. Res. 1995, 33, 193–200. [Google Scholar] [CrossRef]
  16. Beerten, J.; van Hertem, D.; Belmans, R. VSC-MTDC Systems with a Distributed DC Voltage Control-A Power Flow Approach. In Proceedings of the 2011 IEEE Trondheim PowerTech, Trondheim, Norway, 19–23 June 2011. [Google Scholar]
  17. Haileselassie, T.M.; Uhlen, K. Power Flow Analysis of Multi-Terminal HVDC Networks. In Proceedings of the 2011 IEEE Trondheim PowerTech, Trondheim, Norway, 19–23 June 2011. [Google Scholar]
  18. Beerten, J.; Cole, S.; Belmans, R. Implementation Aspects of a Sequential AC/DC Power Flow Computation Algorithm for Multi-Terminal VSC-HVDC Systems. In Proceedings of the 9th IET International Conference on AC and DC Power Transmission, London, UK, 19–21 October 2010. [Google Scholar]
  19. Beerten, J.; Cole, S.; Belmans, R. A Sequential AC/DC Power Flow Algorithm for Networks Containing Multi-Terminal VSC-HVDC Systems. In Proceedings of the IEEE PES General Meeting, Providence, RI, USA, 25–29 July 2010. [Google Scholar]
  20. Beerten, J.; Cole, S.; Belmans, R. Generalized Steady-State VSC MTDC Model for Sequential AC/DC Power Flow Algorithms. IEEE Trans. Power Syst. 2012, 27, 821–829. [Google Scholar] [CrossRef] [Green Version]
  21. Gonzalez-Longatt, F.; Roldan, J.; Charalambous, C.A. Power Flow Solution on Multi-Terminal HVDC Systems: Supergrid Case. In Proceedings of the International Conference on Renewable Energies and Power Quality, Santiago de Compostela, Galicia, Spain, 28–30 March 2012. [Google Scholar]
  22. Hendriks, R.L.; Paap, G.C.; Kling, W.L. Control of a Multi-Terminal VSC Transmission Scheme for Interconnecting Offshore Wind Farms. In Proceedings of the European Wind Energy Conference and Exhibition, Milan, Italy, 7–10 May 2007. [Google Scholar]
  23. Yang, Z.; Zhong, H.; Bose, A.; Xia, Q.; Kang, C. Optimal Power Flow in AC-DC Grids with Discrete Control Devices. IEEE Trans. Power Syst. 2018, 33, 1461–1472. [Google Scholar]
Figure 1. Single-line diagram of a three-bus MT-HVDC system.
Figure 1. Single-line diagram of a three-bus MT-HVDC system.
Applsci 10 00297 g001
Figure 2. Single-line diagram of the modified IEEE 39-bus test system with two four-bus MT-HVDC systems.
Figure 2. Single-line diagram of the modified IEEE 39-bus test system with two four-bus MT-HVDC systems.
Applsci 10 00297 g002
Table 1. Parameters of the AC system.
Table 1. Parameters of the AC system.
BusType V   ( p . u . ) θ   ( ° )   P G   ( M W ) Q G   ( M V A R ) P D   ( M W ) Q D   ( M V A R )
1Slack1.0600.000.000.00
2PV1.00040.0020.0010.00
3PQ45.0015.00
4PQ40.005.00
5PQ60.0010.00
Table 2. Parameters of the AC lines.
Table 2. Parameters of the AC lines.
FromTo R   ( p . u . ) X   ( p . u . )   B   ( p . u . ) S N   ( M V A )
120.020.060.06100
130.080.240.05100
230.060.180.04100
240.060.180.04100
250.040.120.03100
340.010.030.02100
450.080.240.05100
Table 3. Parameters of the converter stations in per-unit.
Table 3. Parameters of the converter stations in per-unit.
Converter Station S N   ( M V A ) R T R X T R B F R T R X T R
11000.00150.11210.08870.00010.1643
21000.00150.11210.08870.00010.1643
31000.00150.11210.08870.00010.1643
Table 4. Power losses coefficients of the converter stations.
Table 4. Power losses coefficients of the converter stations.
Converter Station a b c r e c c i n v
11.1030.8872.8854.371
21.1030.8872.8854.371
31.1030.8872.8854.371
Table 5. Parameters of the DC lines.
Table 5. Parameters of the DC lines.
FromTo R D C   ( p . u . ) V D C   ( k V )   P D C   ( M W )
120.0260345100
130.0365345100
230.0260345100

Share and Cite

MDPI and ACS Style

Mohammadi, F.; Nazri, G.-A.; Saif, M. An Improved Mixed AC/DC Power Flow Algorithm in Hybrid AC/DC Grids with MT-HVDC Systems. Appl. Sci. 2020, 10, 297. https://doi.org/10.3390/app10010297

AMA Style

Mohammadi F, Nazri G-A, Saif M. An Improved Mixed AC/DC Power Flow Algorithm in Hybrid AC/DC Grids with MT-HVDC Systems. Applied Sciences. 2020; 10(1):297. https://doi.org/10.3390/app10010297

Chicago/Turabian Style

Mohammadi, Fazel, Gholam-Abbas Nazri, and Mehrdad Saif. 2020. "An Improved Mixed AC/DC Power Flow Algorithm in Hybrid AC/DC Grids with MT-HVDC Systems" Applied Sciences 10, no. 1: 297. https://doi.org/10.3390/app10010297

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