Abstract
For the planning, operation and control of multi-terminal voltage source converter (VSC) based high-voltage direct current (HVDC) (VSC-MTDC) systems, an accurate power flow formulation is a key starting point. Conventional power flow formulations assume the constant frequencies for all asynchronous AC systems. Therefore, a new feature about the complex coupling relations between AC frequencies, DC voltages and the exchanged power via VSC stations cannot be characterized if VSC-MTDC systems are required to provide cross-regional frequency responses. To address this issue, this paper proposes a comprehensive frequency-dependent power flow formulation. The proposed approach takes the frequencies of asynchronous AC systems as explicit variables, and investigates the novel bus models of the interlinking buses of VSC stations. The proposed approach accommodates different operation modes and frequency droop strategies of VSC stations, and considers the power losses of VSC stations. The effectiveness and generality of the developed approach are validated by a 6-terminal VSC-HVDC test system. The test system presents the characteristics of the coexistence of numerous VSC operation modes, the absence of slack buses in both AC and DC subsystems, and diversified grid configurations such as point-to-point integration of renewable energy sources and one AC system integrated with multiple VSC stations.
Keywords
Voltage source converter based high-voltage direct current (VSC-HVDC) ; multi-terminal VSC based HVDC (VSC-MTDC) ; cross-regional frequency response ; power flow formulation ; VSC operation mode
OVER the years, the increasing electricity demand and the renewable energy integration have greatly promoted the rapid development of high-voltage direct current (HVDC) technologies [
As the local synchronous generators (SGs) are replaced by feed-in DC transmissions, a larger rate of change of frequency (RoCoF), a worse frequency nadir and longer frequency recovery time may occur in one single synchronous system with the decreased rotating inertia and frequency regulation capacity [4]-[6]. Thus, a cross-regional frequency response is imperative for VSC-MTDC systems.
To achieve the cross-regional frequency responses and unbalanced power allocation among different asynchronous AC grids, droop-controlled strategies are usually adopted by VSC stations [7]-[10]. In such a system, the frequency of each asynchronous AC system varies under different operation conditions, which results in that the units or VSC stations may change the output power or exchanged power according to their frequency regulation strategies. As a result, a salient feature is introduced to the VSC-MTDC systems that AC frequencies, DC voltages, the power injected into the point of common connection (PCC) of VSC stations and the output power of units are coupled together. Therefore, an investigation into this prominent coupling feature is necessary to accurately analyze the steady state of the entire VSC-MTDC systems.
The conventional power flow approaches for AC grids [
To investigate a power flow formulation incorporating the above-mentioned frequency coupling feature, [19]-[22] indirectly couple the AC frequency and DC bus voltages by controlling the exchanged power via converters. However, the control strategy used in these studies applies to a certain VSC operation mode, i.e., active power is controlled in the active current channel of VSC controllers, which cannot be generalized to other possible VSC operation modes in meshed MTDC systems. Also, simple network configurations are designed in these studies, e.g., only one AC subsystem integrates with one DC subsystem [19]-[21], or two AC (or DC) subsystems integrate with one DC (or AC) subsystem via one converter separately [
Therefore, a novel comprehensive approach is developed in this paper to formulate the power flow of VSC-MTDC systems considering the cross-regional frequency regulations. The contributions of this paper are as follows: ① different operation modes of VSC stations and versatile frequency droop strategies are accommodated by the proposed approach; ② new bus types are defined to deal with the difficulty that almost all variables related to the PCC and DC interlinking buses of VSC stations are unknown before the PFC; ③ additional power flow variables, i.e., AC system frequencies and the losses of VSC stations are introduced, and the corresponding extra mismatch equations are formulated to guarantee that the number of unknown variables is equal to the number of equations; ④ the challenges are overcame by the proposed power flow approach such as the absence of slack buses in both AC and DC subsystems, and different network configurations, e.g., one AC subsystem integrated with multiple VSC stations.
The remainder of this paper is organized as follows. Section II introduces the model of VSC stations and defines novel bus types considering converter losses and frequency droop strategies. Section III investigates the power flow formulation of VSC-MTDC systems. Section IV presents the case study to demonstrate the effectiveness of the proposed approach. Finally, the conclusions are drawn in Section V.
In the steady-state analysis, a VSC station can be represented by a controllable AC voltage source and a controllable DC current source, as shown in

Fig. 1 Steady-state model of a VSC station.
The PCC and DC interlinking buses are directly involved in the PFC, while the internal buses F and B 1 account for the calculation of VSC station losses. The current injected into the PCC is calculated by:
(1) |
where Ps and Qs are the active power and reactive power injected into the PCC bus, respectively; and Vs and δs are the voltage magnitude and phase angle of the PCC bus, respectively.
The bus voltage of an AC filter Vf , the passing current of a phase reactor Ic , and the power injected into the DC interlinking bus Pc , DC are obtained by (2)-(4), respectively.
(2) |
(3) |
(4) |
where PcL is the total loss of the VSC station, which consists of two parts. The first part is the transformer and phase reactor losses PcL 1, and the second part is the converter loss PcL 2. Equations (6) and (7) indicate that the converter losses are related to the magnitude of the passing currents.
(5) |
(6) |
(7) |
where acL , bcL and ccL are the no-load converter losses, linear and quadratic coefficients of converter losses, respectively.
For the modular multilevel converter (MMC), there is no need to use an AC voltage filter, so the corresponding steady-state model can be simplified as (1), (4)-(7) with .

Fig. 2 Operation modes of a VSC station.
Combining different references from active and reactive current channels, different modes of VSC stations are determined such as P-VAC
, P-Q, VDC
-Q, and VDC
-VAC
. Specifically, VSC stations can also operate in a VAC
-F mode if they are connected to an islanded AC system or an RES power plant, as described in

Fig. 3 VAC -F operation mode of a VSC station.
Based on the different operation modes of VSC stations, versatile frequency regulation strategies can be designed. In this paper, different frequency-related droop strategies are adopted to enable the units in one subsystem providing cross-regional responses to the frequency events in other asynchronous AC subsystems. A VDC -ωs or Ps -ωs droop strategy allows VSC stations to change the DC bus voltage or the power injected into the PCC by perceiving the frequency deviations of the connected AC subsystem with frequency events, and a ωs -VDC or ωs -Ps droop strategy allows VSC stations to change the frequencies of AC subsystems, in which the units can provide frequency responses based on only local measurements, by perceiving the deviations of DC bus voltage or the power injected into the PCC.
Since the converter losses and droop coupling variables of VSC stations are unknown before PFC, four variables of the PCC bus (active power Ps , reactive power Qs , voltage magnitude |Vs | and phase angle δs ), two variables of the DC interlinking bus (active power PDC and bus voltage VDC ), and AC frequency ωs cannot be pre-specified, which contradict the assumption that half of the variables of traditional bus types are known before the PFC such as P-Q, P-V and slack buses in AC systems, or constant V and P buses in DC systems. Thus, new bus types need to be defined.
Taking the converter losses, different VSC operation modes and control strategies into consideration, Table I summarizes the bus type classification of the PCC and DC interlinking buses. In Cases 1-5 and 6-12, VSC stations adopt constant and droop-controlled strategies, respectively. As mentioned before, VSC stations typically operate in a VAC -F mode when they are connected to an islanded AC system or an RES station, so the PCC acts as a V-δ bus in Cases 1 and 6. The bus type classification of PCC and DC interlinking buses is as follows.
In Cases 1-3 and 6-12, since the power or/and voltage of the PCC bus is/are unknown, the converter losses cannot be obtained by (1)-(7), and the power injected into the DC interlinking bus is unknown before the PFC. Therefore, if the VSC stations do not adopt a constant VDC control strategy, the corresponding DC bus is neither a constant P bus nor a VDC bus. However, DC power mismatch equations are necessary to solve the DC power flow variables. Hence, a novel bus type has to be defined, i.e., DC tie bus, to classify the DC interlinking buses in Cases 1-3 and 6-12. The DC tie bus shares the similar power mismatch equation as the constant P bus, and the power injected into the DC interlinking bus needs to be updated during the PFC iterations.
If a VSC station operates in a constant VDC mode or a droop-controlled mode (Cases 4-5 and 7-12), the corresponding PCC bus is neither a P-Q bus nor a P-VAC bus because the converter losses and the active power injected into the PCC are unknown before the PFC. However, AC power mismatch equations are necessary to solve the AC power flow variables. Hence, two novel bus types, i.e., P-V tie bus and P-Q tie bus, are defined to classify the PCC bus if the VSC station adopts VAC and QAC control strategies in the reactive current channel, respectively. The P-V tie bus and P-Q tie bus share the similar power mismatch equations as the P-V bus and P-Q bus, respectively, and the active power injected into the PCC bus needs to be iteratively updated during the PFC iterations.
As described in Section II, since the output power of droop-controlled units and the power injected into the PCC or the DC voltage of droop-controlled VSC stations are coupled with AC frequencies, the frequencies of different asynchronous systems inevitably become variables. In addition, the converter losses cannot be pre-specified before PFC. To update the power injected into the PCC and the DC interlinking buses, additional state variables are introduced to represent the converter losses.
To solve these additional state variables, extra mismatch equations need to be formulated. The mismatch equations of the entire VSC-MTDC system are extended to four groups: AC power mismatch equations, DC power mismatch equations, active power-balance mismatch equations and droop-controlled mismatch equations.
A unified power flow formulation is described mathematically by a set of linear and nonlinear equations f ( x ) in a polar form, as shown in (8).
(8) |
(9) |
where f is the vector containing mismatch equations including AC mismatch equation f AC , DC mismatch equation f DC and VSC station mismatch equation f ACDC ; x is the vector containing the iterative power flow variables, which are divided into 3 parts, i.e., AC iterative variable x AC , DC iterative variable x DC and converter iterative variable x cL ; x AC contains AC bus voltage phase angle δ , voltage magnitude | V | and system frequency ω s ; x DC contains DC bus voltage V DC ; and x cL contains converter loss P cL .
Collecting all the AC active power mismatch equation Δ P AC , reactive power mismatch equation Δ Q AC , DC power mismatch equation Δ P DC , active power-balance mismatch equation Δ P ACDC and droop-controlled mismatch equation Δ D ACDC , yields the detailed power flow equation in (10) expressed by a modified Jacobian matrix J . The detailed expressions of these mismatch equations are given in the following subsections.
(10) |
Equation (10) is capable of analyzing VSC-MTDC systems with arbitrary network topologies and versatile control strategies of VSC stations. If multiple VSC stations are integrated into one AC system, all except one of the corresponding connected DC bus voltages need to be taken out of the matrix in (10) to maintain an equal number of variables and equations.
For a general AC bus, the power mismatch equations are:
(11) |
where PGi and QGi are the generated active power and reactive power, respectively; Psi and Qsi are the active power and reactive power injected into the PCC bus i from the connected VSC station, respectively, which are equal to zero if the bus i is not connected to a VSC station; PLi and QLi are the load active power and reactive power, respectively; and Pcalc,i and Qcalc,i are the calculated active power and reactive power at bus i, respectively, which can be obtained by:
(12) |
where |Vi | and δi are the voltage magnitude and phase angle of bus i, respectively; θij and Yij are the phase of line impedance and admittance between bus i and bus j, respectively; and NAC is the number of AC buses.
For VSC-MTDC systems, AC buses
B
AC
are classified into 6 types: P-V buses
B
PV
, P-Q buses
B
PQ
, V-δ buses
B
Vδ
, AC generator droop buses
B
ACd
, P-V tie buses
B
PVt
, and P-Q tie buses
B
PQt
. The denotation
B
represents the set of buses. The first two types of buses are the same with the conventional AC power flow approach [
(13) |
where |V 0| and ωs 0 are the nominal voltage magnitude and frequency, respectively; PGi, 0 and QGi, 0 are the nominal active power and reactive power, respectively; and mpi and mqi are the frequency and voltage droop coefficients, respectively.
Therefore, the mismatch equations of AC subsystems can be summarized by (14), where the detailed expression of f AC is shown in (15).
(14) |
(15) |
2) Mismatch Equations of DC Subsystem
For a general DC bus, the power mismatch equation is:
(16) |
where PG,DCi is the generated active power; Pc,DCi is the active power injected into the DC interlinking bus i from the connected VSC station, which equals to 0 if the bus i is not connected with a VSC station; PL,DCi is the active power demand; and Pcalc,DCi is the calculated active power at the bus i that can be obtained by:
(17) |
where VDCi , GDC,ij and NDC are the DC voltage of the DC bus i, the conductance between the DC buses i and j, and the number of DC buses, respectively.
DC buses
B
DC
are classified into 4 types: VDC
buses
B
V
, P buses
B
P
, DC tie buses
B
Pt
, and DC generator droop buses
B
DCd
. VDC
buses and P buses are the same with the conventional DC power flow approach [
(18) |
where VDCi, 0 and PG,DCi, 0 are the nominal DC voltage and active power, respectively; and mdcpi is the voltage droop coefficient.
Therefore, the mismatch equations of DC subsystems can be summarized by (19), where the detailed expression of f DC is described by (20).
(19) |
(20) |
3) Mismatch Equations of VSC Stations
The mismatch equations of VSC stations consist of two parts: active power-balance mismatch equations and droop-controlled mismatch equations.
The active power-balance mismatch equation is shown in (21), and the corresponding detailed expression is described in (22).
(21) |
(22) |
where B VSC is the PCC or DC interlinking buses that are connected to VSC stations; and Psi , PcLi , Pc , DCi are the power injected into the PCC calculated by (12), the converter losses calculated by (5), and the power injected into the DC interlinking bus calculated by (17), respectively.
The droop-controlled mismatch equations of the VSC stations adopting droop strategies are presented in (23). Different types of droop relations such as VDC -ωs , Ps -ωs , ωs -VDC , ωs -Ps , VDC -PDC and VDC -IDC , are expressed by (24)-(29), respectively.
(23) |
(24) |
(25) |
(26) |
(27) |
(28) |
(29) |
where
B
Pt
is the PCC buses or the DC interlinking buses that are connected to the VSC stations adopting droop strategies,
B
Pt
⊆
B
VSC
; VDCi
, Psi
, ωsi
are the DC bus voltage, the power injected into the PCC and AC system frequency at the interlinking bus i, respectively; VDCi
,0, Psi
,0, ωsi
,0 are the nominal DC bus voltage, the power injected into the PCC and the AC system frequency at the interlinking bus i, respectively; and kVdcωi
, kPωi
, kωVdci
, kωPi
, kVdcPdci
, kVdcIdci
are the corresponding coefficients of VDC
-ωs
, Ps
-ωs
, ωs
-VDC
, ωs
-Ps
, VDC
-PDC
and VDC
-IDC
droop relations, respectively. In addition, VSC stations can also adopt nonlinear droop relations between VDC
or Ps
and ωs
, PDC
or IDC
[
Therefore, the mismatch equations of VSC stations are summarized by (30), where the detailed expression of f ACDC is described by (31).
(30) |
(31) |
The flow chart of the proposed power flow formulation in VSC-MTDC systems is depicted in

Fig. 4 Flow chart of the proposed power flow formulation of VSC-MTDC systems with cross-regional frequency responses.
The proposed power flow formulation is solved by the numerical methods such as a globally convergent trust region method [
The variables and mismatch equations for solving the power flow of VSC-MTDC systems are summarized in Table II, where NP , NPt and NDCd are the numbers of P, DC tie and DC droop buses, respectively; NPV , NPQ , NPVt , NPQt , NACd are the numbers of P-V, P-Q, P-V tie, P-Q tie and AC droop buses, respectively; and NVSC and NINT are the numbers of VSC stations and asynchronous AC systems, respectively.
A 6-terminal VSC-HVDC test system is applied in this paper, as described in

Fig. 5 6-terminal VSC-HVDC test system.
The parameters of the test system can be found in Table III and Table IV. The resistance and reactance of AC cables are RAC = 0.01 Ω/km and XAC = 0.04 Ω/km (nominal frequency is 50 Hz), respectively. The resistance of DC cables RDC is equal to 0.02 Ω/km. For VSC stations, Zc = 0.0001+j0.15, ZT = 0.001+j0.1. VSC stations are represented by an average model, and their AC filters are ignored. The base capacity of the entire system Sbase is equal to 300 MVA. The base AC and DC voltages are 220 kV and 400 kV, respectively. For each droop generator bus connected to SGs and WTs, mp = 20, mq = 0.
The seven cases listed in Table V are studied. Case 1 provides a benchmark that the units and VSC stations do not provide any frequency responses, and AC frequencies are set to the nominal values. In Cases 2.1, 2.2, 3.1, 3.2, 4.1 and 4.2, different droop strategies are adopted. In AC1, VSC station 1 adopts the VDC -ωs (Cases 2.1 and 2.2), Ps -ωs (Cases 3.1 and 3.2) or VDC -PDC (Cases 4.1 and 4.2) droop strategy, and SG1 adopts the Pm -ωs droop strategy, where Pm is the input mechanical power of SGs. In AC2-AC4, VSC stations adopt the ωs -VDC droop strategy and WFs adopt the Pe -ωs droop strategy, where Pe is the output electromagnetic power. In AC5, VSC stations adopt the VDC -ωs droop strategy, and SG5 adopts the Pm -ωs droop strategy. Different locations of load changes are also considered. In Cases 2.1, 3.1 and 4.1, the load power is increased by 0.2 p.u. at AC5 bus 3 (ΔPL 51 =0.2 p.u.). In Cases 2.2, 3.2 and 4.2, the load power is increased by 0.2 p.u. at AC1 bus 2 (ΔPL 1 =0.2 p.u.). In addition, Case 1 provides the initial references to all droop-controlled units and VSC stations in Cases 2.1-4.2.
According to the bus type classification introduced in Section II, the types of the PCC buses that are connected to the VSC stations 1-6 are P-V tie, V-δ, V-δ, P-V tie, P-Q tie and V-δ, respectively. Taking Case 2.1 as an example, Table VI compares the results of the proposed power flow approach and the steady-state results obtained by the time-domain simulation. In Table VI, the maximum errors of voltage magnitude, phase angle and system frequency are less than 0.0017, 0.0007 and 0.00006, respectively. The closely-matched results validate the effectiveness of the proposed approach.
C. Impact of Droop Coefficients and Relations in VSC Stations on Cross-Regional Frequency Regulation
To investigate the impact of droop-controlled VSC stations on the power flow and load allocation in VSC-MTDC systems, different droop coefficients of ωs
-VDC
(VSC stations 2, 3 and 6) and VDC
-ωs
(VSC stations 1, 4 and 5) are compared in Figs. 6-9. With the increase of kωVdc
and kVdcω
, the system frequency deviations of AC1 and AC5 are reduced, and those of AC2-AC4 are increased as shown in

Fig. 6 Frequency of AC subsystems with different kωVdc in Case 2.1 (ΔPL 51 = 0.2 p.u.).

Fig. 8 Frequency of AC subsystems with different kVdcω in Case 2.1 (ΔPL 51 = 0.2 p.u.).

Fig. 7 Generated power of SGs and WFs with different kωVdc in Case 2.1 (ΔPL 51 = 0.2 p.u.).

Fig. 9 Generated power of SGs and WFs with different kVdcω in Case 2.1 (ΔPL 51 = 0.2 p.u.).
In
VSC stations adopting different droop strategies may present the distinctive cross-regional frequency regulation capabilities in VSC-MTDC systems.
The VSC station adopting a VDC -ωs droop strategy takes both AC and DC variables into account. The units in the corresponding connected AC subsystem can provide responses to the frequency events in other asynchronous AC subsystems, and the units in other AC subsystems can also participate in the frequency regulations of this AC subsystem. In Case 2.2, with the increase of kVdcω of VSC station 1, the units in other subsystems increase their output power and provide the responses to the frequency event in AC1, as shown in Figs. 10 and 11.
The VSC station adopting a Ps
-ωs
droop strategy takes only AC variables into account. The units in the corresponding connected AC subsystem cannot respond to the frequency events in other AC subsystems, but the units in other AC subsystems are capable of providing frequency responses to this subsystem. In Cases 3.1 and 3.2, the VSC station 1 adopts a Ps
-ωs
droop strategy. In Case 3.1, SG1 does not provide any response to the frequency event in AC5 no matter what value kPω
of the VSC station 1 is taken, as shown in

Fig. 11 Generated power of SGs and WFs with different kVdcω in Case 2.2 (ΔPL 1 = 0.2 p.u.).

Fig. 12 Frequency of AC subsystems with different kPω in Case 3.1 (ΔPL 51 = 0.2 p.u.).

Fig. 13 Frequency of AC subsystems with different kPω in Case 3.2 (ΔPL 1 = 0.2 p.u.).

Fig. 14 Generated power of SGs and WFs with different kPω in Case 3.2 (ΔPL 1 = 0.2 p.u.).
The VSC station adopting a VDC
-PDC
droop strategy takes only DC variables into account. In Cases 4.1 and 4.2, the VSC station 1 adopts a VDC
-PDC
droop strategy. In Case 4.1, a frequency event occurs in AC5. With the increase of |kVdcPdc
| of VSC station 1, only the frequency of the AC1 is restored, but the frequency of the AC5 is deteriorated even though the WFs in the AC2-AC4 increase their output power, as shown in

Fig. 15 Frequency of AC subsystems with different kVdcPdc in Case 4.1 (ΔPL 51 = 0.2 p.u.).

Fig. 16 Generated power of SGs and WFs with different kVdcPdc in Case 4.1 (ΔPL 51 = 0.2 p.u.).

Fig. 17 Frequency of AC subsystems with different kVdcPdc in Case 4.2 (ΔPL 1 = 0.2 p.u.).
Compared to the Ps
-ωs
droop strategy in

Fig. 10 Frequency of AC subsystems with different kVdcω in Case 2.2 (ΔPL 1 = 0.2 p.u.).
A central challenge of analyzing VSC-MTDC systems lies in the formulation and computation of power flow with the coupling relations between AC frequencies, DC voltages, the power injected into the PCC of VSC stations and the output power of the units. This paper revisits the conventional assumptions associated with power flow problems, and introduces a suitable approach of frequency-dependent power flow formulations. The proposed approach can be generalized to analyze the steady states of VSC-MTDC systems with arbitrary system topologies, operation modes and the frequency droop strategies adopted by VSC stations.
The time-domain simulation results obtained by MATLAB/Simulink validate that the proposed approach can precisely analyze the power flow under different operation conditions. Simulation results reveal that the VDC -ωs droop strategy adopted by VSC stations presents the best cross-regional frequency regulation performance compared with the Ps -ωs and VDC -PDC droop strategies. The results also illustrate that larger droop coefficients of VDC -ωs and Ps -ωs may enhance the cross-regional frequency regulation performance. In contrast, the VDC -PDC droop strategy may deteriorate the cross-regional frequency regulation in VSC-MTDC systems because it focuses only on the power allocation of DC systems, regardless of the frequency performances of the AC systems.
The proposed model can further aid the design, operation and planning of VSC-MTDC and hybrid LCC-VSC HVDC systems. The operation modes and droop coefficients of VSC stations can be optimized to trade off the economic dispatch and operation stability of the entire system. In addition, improved trust-region and Newton-Raphson methods can be investigated to solve the power flow problem of larger systems with complicated topologies.
REFERENCES
E. Pierri, O. Binder, N. G. Hemdan et al., “Challenges and opportunities for a European HVDC grid,” Renewable and Sustainable Energy Reviews, vol. 70, pp. 427-456, Apr. 2017. [百度学术]
N. Flourentzou, V. G. Agelidis, and G. D. Demetriades, “VSC-based HVDC power transmission systems: an overview,” IEEE Transactions on Power Electronics, vol. 24, no. 3, pp. 592-602, Mar. 2009. [百度学术]
T. An, G. Tang, and W. Wang, “Research and application on multi-terminal and DC grids based on VSC-HVDC technology in China,” High Voltage, vol. 2, no. 1, pp. 1-10, Jun. 2017. [百度学术]
E. Ørum, M. Kuivaniemi, M. Laasonen et al., “Future system inertia,” [Online]. Available: https://docstore.entsoe.eu/Documents/Publications/SOC/Nordic/Nordic_report_Future_System_Inertia.pdf [百度学术]
L. Che, X. Shen, and M. Shahidehpour, “Primary frequency response based rescheduling for enhancing microgrid resilience,” Journal of Modern Power Systems and Clean Energy, vol. 7, no. 4, pp. 696-704, Jul. 2019. [百度学术]
Z. Wu, W. Gao, T. Gao et al., “State-of-the-art review on frequency response of wind power plants in power systems,” Journal of Modern Power Systems and Clean Energy, vol. 6, no. 1, pp. 1-16, Jan. 2018. [百度学术]
K. Rouzbehi, A. Miranian, J. I. Candela et al., “A generalized voltage droop strategy for control of multiterminal DC grids,” IEEE Transactions on Industry Applications, vol. 51, no. 1, pp. 607-618, Jan. 2015. [百度学术]
M. A. Abdelwahed and E. F. El-Saadany, “Power sharing control strategy of multiterminal VSC-HVDC transmission systems utilizing adaptive voltage droop,” IEEE Transactions on Sustainable Energy, vol. 8, no. 2, pp. 605-615, Apr. 2017. [百度学术]
W. Wang, Y. Li, Y. Cao et al., “Adaptive droop control of VSC-MTDC Ssystem for frequency support and power sharing,” IEEE Transactions on Power Systems, vol. 33, no. 2, pp. 1264-1274, Mar. 2018. [百度学术]
Y. Li, Z. Xu, J. Østergaard et al., “Coordinated control strategies for offshore wind farm integration via VSC-HVDC for system frequency support,” IEEE Transactions on Energy Conversion, vol. 32, no. 3, pp. 843-856, Sept. 2017. [百度学术]
J. J. Grainger and W. D. Stevenson, Power System Analysis, New York, USA: McGraw-Hill, 1994. [百度学术]
H. Fudeh and C. M. Ong, “A simple and efficient AC-DC load-flow method for multiterminal DC systems,” IEEE Transactions on Power Apparatus and Systems, vol. PAS-100, no. 11, pp. 4389-4396, Nov. 1981. [百度学术]
T. Smed, G. Andersson, G. B. Sheble et al., “A new approach to AC/DC power flow,” IEEE Transactions on Power Systems, vol. 6, no. 3, pp. 1238-1244, Aug. 1991. [百度学术]
C. Liu, B. Zhang, Y. Hou et al., “An improved approach for AC-DC power flow calculation with multi-infeed DC systems,” IEEE Transactions on Power Systems, vol. 26, no. 2, pp. 862-869, May 2011. [百度学术]
J. Beerten, S. Cole, and R. Belmans, “Generalized steady-state VSC MTDC model for sequential AC/DC power flow algorithms,” IEEE Transactions on Power Systems, vol. 27, no. 2, pp. 821-829, May 2012. [百度学术]
M. Baradar and M. Ghandhari, “A multi-option unified power flow approach for hybrid AC/DC grids incorporating multi-terminal VSC-HVDC,” IEEE Transactions on Power Systems, vol. 28, no. 3, pp. 2376-2383, Aug. 2013. [百度学术]
W. Wang and M. Barnes, “Power flow algorithms for multi-terminal VSC-HVDC with droop control,” IEEE Transactions on Power Systems, vol. 29, no. 4, pp. 1721-1730, Jul. 2014. [百度学术]
S. Khan and S. Bhowmick, “A comprehensive power-flow model of multi-terminal PWM based VSC-HVDC systems with DC voltage droop control,” International Journal of Electrical Power and Energy Systems, vol. 102, pp. 71-83, Nov. 2018. [百度学术]
A. A. Hamad, M. A. Azzouz, and E. F. El-Saadany, “A sequential power flow algorithm for islanded hybrid AC/DC microgrids,” IEEE Transactions on Power Systems, vol. 31, no. 5, pp. 3961-3970, Sept. 2016. [百度学术]
A. Eajal, M. A. Abdelwahed, E. El-Saadany et al., “A unified approach to the power flow analysis of AC/DC hybrid microgrids,”" IEEE Transactions on Sustainable Energy, vol. 7, no. 3, pp. 1145-1158, Jul. 2016. [百度学术]
M. E. Nassar, A. A. Hamad, M. Salama et al., “A novel load flow algorithm for islanded AC/DC hybrid microgrids,” IEEE Transactions on Smart Grid, vol. 10, no. 2, pp. 1553-1566, Mar. 2019. [百度学术]
E. Aprilia, K. Meng, M. Al Hosani et al., “Unified power flow algorithm for standalone AC/DC hybrid microgrids,” IEEE Transactions on Smart Grid, vol. 10, no. 1, pp. 639-649, Jan. 2019. [百度学术]
W. Zhang, A. M. Cantarellas, J. Rocabert et al., “Synchronous power controller with flexible droop characteristics for renewable power generation systems,” IEEE Transactions on Sustainable Energy, vol. 7, no. 4, pp. 1572-1582, Oct. 2016. [百度学术]
F. Mumtaz, M. Syed, M. Al Hosani et al., “A novel approach to solve power flow for islanded microgrids using modified newton raphson with droop control of DG,” IEEE Transactions on Sustainable Energy, vol. 7, no. 2, pp. 493-503, Apr. 2016. [百度学术]
M. M. A. Abdelaziz, H. E. Farag, E. F. El-Saadany et al., “A novel and generalized three-phase power flow algorithm for islanded microgrids using a newton trust region method,” IEEE Transactions on Power Systems, vol. 28, no. 1, pp. 190-201, Feb. 2013. [百度学术]
A. Elrayyah, Y. Sozer, and M. E. Elbuluk, “A novel load-flow analysis for stable and optimized microgrid operation,” IEEE Transactions on Power Delivery, vol. 29, no. 4, pp. 1709-1717, Aug. 2014. [百度学术]
R. T. Pinto, P. Bauer, S. F. Rodrigues et al., “A novel distributed direct-voltage control strategy for grid integration of offshore wind energy systems through MTDC network,” IEEE Transactions on Industrial Electronics, vol. 60, no. 6, pp. 2429-2441, Jun. 2013. [百度学术]
A. R. Conn, N. I. Gould, and P. L. Toint, “Trust region methods,” MOS-SIAM Series on Optimization. [Online]. Available: https://doi.org/10.1137/1.9780898719857 [百度学术]