3.1. System Description and Mathematical Model of the Single VSC
In this paper, the transient stability of the multi-VSC parallel system is targeted. To coherently derive a reduced order model for the multi-VSC parallel system, taking a single VSC grid-connected system as an example, the topology and control strategy of the adopted grid-following VSC are introduced, and the reduced-order model of the single VSC is obtained.
Figure 3 represents the simplified diagram and the adopted control method for each VSC. For VSC in grid-following mode, there is usually a large capacitor or a control part adopted to stabilize the DC voltage, so the DC side can be considered an ideal voltage source, the voltage of which can be represented as
On the AC side, the VSC is filtered by the inductor (
). After that, the VSC is connected to the weak grid by line impedance (
). The weak grid consists of the ideal voltage source (
) and grid impedance (
). The amplitude and phase of the grid voltage
can be represented as
and
, respectively. In terms of the control strategy, the VSC measures the three-phase current (
) and the point of common coupling (PCC) voltage (
) to the PLL to obtain a phase
. Through Park transformation,
can be transferred to the active current (
) and the reactive current (
). The active current reference is
, and the reactive current reference is
. The current is regulated by the PI controller to achieve zero-steady-state error.
Figure 4 shows the control principle of the PLL, where
and
represent the proportional and integral parameters of the PI controller. Due to the presence of PI control and frequency integration, the disturbance of
will be transmitted to
through two integrators. Therefore, the PLL cannot achieve instant phase locking. In addition, due to the presence of line impedance and weak grid impedance, there is complex coupling between different parallel VSCs as well as between the VSC and the grid, making the transient stability analysis of multi-VSC systems more complex.
To facilitate modeling, the amplitude and phase representation is used for impedance, voltage, and current, which can be represented as
The utility grid composed of an ideal voltage source and grid impedance in series can be converted into the current source form, where
The bus voltage can be represented as
Thus, the voltage of the PCC can be represented as
For simplification, define
Combining (18)–(24), the projection of PCC voltage on the phase-locked angle can be derived as
where
is the angular frequency and
is the rated angular frequency value.
indicates the deviation between the system frequency and the rated frequency. Power angle (
) can be defined as
The control principle of PLL is to achieve phase lock using a PI controller to make the projection value of the q-axis PCC voltage zero. The q-axis voltage can be derived as follows:
where
represents taking the imaginary part. Generally, the system frequency is relatively close to the rated frequency, thus
When displays the inductive characteristic, and . When displays the capacitive characteristic, and
To simplify the expression and highlight the components of the q-axis voltage projection value, define the constant component as
, thus the
q-axis PCC voltage can be represented as
According to the PLL control principle shown in
Figure 4, the state-space model can be represented as
where
is the space variable introduced by the PLL integrator controller.
According to (26), considering the PLL control strategy, expressions about power angle can be derived as
Combining (29) and (31), this can be derived:
Combining (30)–(32), the reduced-order big-signal model for the single VSC system can be developed as
The big-signal model can be represented in
Figure 5.
Trigonometric functions in large-signal models pose challenges for transient stability analysis. Consequently, Taylor expansion is applied to derive the second-order and cubic-order truncated models based on the reduced-order large-signal model.
The state variable can be expressed as a combination of the steady-state quantity and the incremental value, as
, where
and
represent the steady-state and the incremental value, respectively. Taylor’s expression about the trigonometric function
is as follows:
where coefficients are
.
Combining (33) and (34) by truncating the system model to the quadratic term from the steady-state operating point, the reduced-order quadratic truncation model for a single VSC grid-connected system can be derived as follows:
Similarly, by truncating the system model to the cubic term, the reduced-order cubic truncation model for a single VSC grid-connected system can be derived as follows:
Thus, the reduced-order large-signal model, along with the corresponding quadratic and cubic truncation models, for the single VSC grid-connected system is derived.
In the case of single VSC grid-connected systems, the reduced-order large-signal model is second-order, enabling direct application of mature phase diagram or inverse trajectory methods for transient stability analysis. However, as the number of VSC devices connected to the system increases, the number of system nodes grows, and the system order increases significantly. Consequently, the phase diagram method and inverse trajectory method used for analyzing low-order systems cannot be directly applied to assess the transient stability of the system.
3.3. Reduced-Order Derivation for the Multi-VSC Parallel System
Currently, most renewable energy sources are connected to the power grid through grid-following VSCs. Thus, this paper primarily focuses on modeling and transient stability analysis for scenarios involving multiple grid-following VSCs connected to the grid. Since the VSCs employ a similar control strategy, the circuit model and control strategy of multi-VSCs connected to the grid will not be extensively discussed further. Each grid-connected VSC will be labeled sequentially as
,
. The symbol representation of circuit parameters and control parameters in a multi-VSC grid-connected system follows the same as those defined in a single VSC grid-connected system.
and
respectively denote the filtering inductance and line impedance associated with the
k-th VSC circuit. Other parameters follow a similar pattern and will not be reiterated. To simplify the circuit topology of a multi-VSC grid-connected system, a similar circuit transformation is employed to convert the series combination of the ideal voltage source and grid impedance into a current source form.
Figure 8 presents the corresponding schematic for the simplified circuit topology of the multi-VSC grid-connected system.
Bus voltage in the multi-VSC grid-connected system can be derived as
The measured voltage at the PCC of each VSC, which is used for the PLL of each VSC itself, can be expressed as
Each VSC’s PLL is regulated by locally measuring voltage information. As a result of different line impedances, the active and reactive current reference values for VSCs and angles detected by VSCs are different.
Figure 9 displays the phase angle diagram for the grid-connected system with multiple VSCs.
For simplification, define
Thus, the voltage of PCC can be derived as
The control principle of the PLL in each VSC involves using a PI controller to make the
q-axis voltage zero, ultimately determining the phase at PCC. Therefore, the q-axis voltage at each PCC can be expressed as
Thus, the
q-axis voltage can be derived as
where
represents the power angle of the
k-th VSC. According to (42), the PCC voltage of the
k-th VSC is influenced not only by its own loop’s current reference value as well as the grid current and grid impedance but also by the current and impedance on other VSC lines. It is illustrated that coupling of the frequency and phase between multiple VSCs occurs and that the complexity of the whole system is dramatically improved.
In the multi-VSC grid-connected system, considering the control principle of PLL, similar state-space model can be developed for the
k-th VSC:
where
represents the state variable introduced by the integrator controller in the k-th VSC.
Taking the example of a grid-connected system comprising two VSCs, the corresponding large-signal block diagram is depicted in
Figure 10.
The multi-VSC grid-connected system, as depicted in
Figure 10, exhibits high-order nonlinear characteristics and involves complex phase angle nonlinear coupling in its large signal model. Their inherent traits pose challenges in directly analyzing the transient stability of the multi-VSC parallel system.
3.4. System Transient Stability Analysis Based on CFND
This paper illustrates a three-VSC parallel system to analyze the transient stability. All state variables can be expressed as steady-state quantities (denoted by the symbol
) superimposed on the variable quantities (denoted by
), such as
. Taylor expansion can be utilized for a model approximation when dealing with nonlinear functions, such as trigonometric functions, in large-signal models. Due to the comparison of the ROAs obtained from the single VSC reduced-order model and the quadratic and cubic truncation models in the previous section, it is shown that truncating to the quadratic term is sufficient to reflect the characteristics of the original system. Therefore, in the analysis of the multi-VSC grid-connected system, the truncated to quadratic Taylor expansion is used for model approximation. The relevant expansion equations are
In general, the frequency difference of VSCs in the system is slight. For simplification, if
, it can be derived that
. When
is inductive,
, where
is the equivalent inductance. When
is capacitive,
, where
is the equivalent capacitance [
22]. For the sake of simplification, this paper mainly addresses the transient stability problem with slow outer loop dynamic characteristics, thus ignoring the influence of the fast current inner loop and the line dynamics [
17]. Combining (42), (44) and (45), the expression of the PCC point voltage truncated to the quadratic term after Taylor transformation is as follows:
Combining (43) and (46), the nonlinear model truncated to the quadratic term of the three-VSC parallel system is developed. The form of the model is (2), and the proposed nonlinear decoupling method is utilized for stability analysis. The parameters adopted by the system are shown in
Table 1. The state variables of the system can be expressed as
To verify the proposed method’s effectiveness, transient stability assessments are conducted for three cases with different operation points.
CASE I: For the case when the grid voltage amplitude suddenly drops from 311 V (1.0 pu) to 31.1 V (0.1 pu), the state equations after decoupling can be obtained, and the initial operation point (
) can be calculated, as shown in the
Table 2.
Since most of the coefficients in decoupled equations are complex, separating the real and imaginary parts for stability analysis in the real number domain is helpful. For quadratic differential equations, the mature trajectory reversing method and the phase diagram can be used to obtain the ROA.
The ROA of the coupling pair
is shown inside the blue curve in
Figure 11a, and the phase diagram is shown in
Figure 11b. The system gradually stabilizes when the initial point lies within the blue curve range. Conversely, when the initial point is outside the ROA, the system diverges and cannot be stabilized. As
Figure 11 shows, the two initial points are marked by red triangles outside the stability region, so the analysis illustrates that the system faces large-signal instability when the voltage drops from 311 V to 31.1 V.
Furthermore, the phase diagrams of the coupling pairs
and
,
and
, as well as the corresponding initial points, are shown in
Figure 12. The phase diagrams illustrate that these two pairs of state variables are globally stable.
CASE II: For the situation when the grid voltage amplitude drops from 40 V to 31.1 V, the system parameters and steady-state operation point are the same as in case I. Therefore, the ROAs corresponding to the two cases are consistent.
As shown in
Figure 11a,b, the green triangles indicate the initial points when the grid voltage amplitude drops from 40 V to 31.1 V. It can be concluded that the initial points are inside the ROA. Thus, there will be no instability when the grid amplitude drops from 40 V.
CASE III: Compared with the transient process of sudden voltage drop from 311 V to 31.1 V, the model parameters (voltage amplitude), the steady-state operation point, and the initial operation point significantly change when the grid voltage rises abruptly from 31.1 V to 311 V. The eigenvalues corresponding to the model are all real roots. The nonlinear decoupling method combines the two pairs of state variables with large coupling factors into coupled pairs. The remaining variables with limited coupling factors are kept as isolated variables. The phase diagrams of the coupling pairs
and
and
are shown in
Figure 13. The phase diagrams of these two coupling pairs are stable over a wide range, which prevents instability. The red triangles are still used to mark the initial positions of the corresponding state variables.
For isolated variables
and
, the diagram of
can be calculated to view the ROA. As shown in
Figure 14a, the curve intersects the axis at points
a and
b, and the red triangle marks the initial point. It can be easily proved that the ROA of
is
, which means that when the initial point is in the interval, the state variable returns to the stable operating point. According to
Figure 14b, the ROA of the state variable
is
, and the initial point lies within the region. Hence, for a sudden voltage rise from 31.1 V to 311 V, all six variables are in stable regions, so the system is robust against large-signal instability.
A comparison of the above three cases indicates that large-signal instability occurs in the system when the grid voltage amplitude drops from 311 V to 31.1 V, but the system can maintain stability in the process of recovery from 31.1 V to 311 V. The cases also confirm that the large-signal stability of the nonlinear system needs to be analyzed by combining the initial and operation points.