1. Introduction
Grid-following (GFL) inverters are the fundamental elements of inverter-based distributed energy resources in the future, modern power systems. Currently, GFL inverters are increasingly adopted as the power conversion interface for grid integration of renewable energy resources [
1]. In this regard, different studies are being conducted on current, power, and the DC-link voltage control of GFL inverters, with the focus on grid synchronization, especially in weak grid conditions with/without phase locked loop (PLL) [
2,
3,
4,
5,
6,
7], the stability of the phase locked loop (PLL), faults right through the capability of th eGFL inverters and their operation under faulty and unbalanced grid conditions [
8,
9,
10,
11], reliability enhancement and cost minimization with a reduction in the number of the required sensors in the control loops [
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22].
The grid integration of a GFL inverter is mainly done through L- or LCL-type filters to limit the harmonic current injection. An L-type filter has a weak harmonic attenuation characteristic, compared to the LCL filters. However, control with an L-type filter is simple, which makes it the preferred option for the grid integration of the GFL inverters. Typically, in the control of the GFL inverters with an L-type filter, at least two AC voltages, two AC currents and one DC-link voltage sensors are required. On the other hand, an LCL filter has superior attenuation characteristics compared to L-type filters, but at the cost of a complex control structure to damp its resonance characteristics. To add more damping to the system, further sensors are required in the control structure.
In the case of sensor failures, the inverter usually loses its stability. Thus, proposing a control structure with a reduced number of the sensors is a key solution in improving the reliability of GFL inverters. In addition, the need to use a reduced number of sensors considerably decreases the cost and the complexity of the system. In this regard, different studies in the literature have tried to propose a control structure with a minimum number of the required sensors [
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22]. Specifically, active damping and control of the GFL inverters with an LCL filter requires the use of multiple signals, which necessitates the use of a reduced number of sensors for the sake of cost and complexity reduction. In this respect, using state observer-based control schemes are presented in the existing studies in the literature [
12,
13,
14,
15]. While the studies in [
12,
13] use only grid-side current sensors and grid voltage, the presented methods in [
14,
15] also eliminated the necessity for a grid voltage sensor. The conventional synchronization schemes of GFLs mostly rely on the use of measured grid-voltages. To remove the grid-voltage measurement from the control structure, virtual-flux observer-based [
19] and grid voltage observer-based [
16,
17,
18] methods are presented. However, current sensors are usually more expensive and their accuracy has more impact on the system-control performance. Therefore, the elimination of the grid current sensors has considerable effect on the cost reduction, reliability enhancement and fault tolerant capability under sensor failure. An observer-based robust current sensor-less control scheme is presented in [
20] for the inverter-based multi-terminal direct current transmission system, which eliminates the current control loops and current measurements. In addition, reference [
21] proposed a new control scheme for a double-fed induction generator, which substitutes the stator-current regulators for the rotor-current regulators, and eliminates the necessity for rotor-current measurement. A current sensor-less control scheme is presented in [
22], for the control of a permanent magnet synchronous motor. Despite the different available studies, the presented methods are not applicable for GFL inverters.
To the best of the authors’ knowledge, measuring the output current of the inverter is required in all of the presented control schemes in the literature, and there is no grid-side current sensor-less control scheme for GFL inverters. To cover this gap, this paper proposes a new scheme that controls the GFL inverter without requiring any grid-side current measurement. In the proposed control scheme, only the grid-side and DC-link voltages are used in the control system. Considering the obtained dynamic interaction between the grid-side current and the DC-link voltage, the grid-side currents are estimated via an observer-based controller. Accordingly, the measured currents are used in the control loop of the system. The novelties and contributions of the paper are summarized as follows.
Proposing an observer for the estimation of the grid-side currents of a GFL-inverter considering the nonlinear dynamics of DCLV;
Including the nonlinear model of DC-link voltage dynamic and the filter inductance instantaneous power in the estimation process;
Providing the well-proved eigenvalue analysis to verify the overall stability of the GFL inverter control system with the proposed observer.
Using the proposed scheme reduces the overall cost of manufacturing the GFL inverter and enhances its reliability and fault-tolerant capability.
The different parts of the paper are organized in seven complementary sections. In
Section 2, the modeling of a GFL inverter with emphasize on the DCLV dynamic is presented. Conventional control of GFL inverter is explained in
Section 3. III. The proposed scheme for the estimation of the grid-side currents is described in
Section 4. The overall state-space model of the GFL inverter, with both the conventional and the proposed control schemes, is obtained in
Section 5, and the stability of the overall system is studied in possible operating points. Various case studies and time-domain simulations are presented in
Section 6, to show the effectiveness of the proposed scheme under different scenarios and dynamic conditions. Finally, the paper is concluded in
Section 7.
2. Modeling of Grid-Following Inverter with L-Type Filter
Figure 1 shows a three-phase diagram of a GFL inverter with an L-type filter which delivers the power from the primary source to the grid. The state-space dynamic equation of the system consists of the dynamics of the grid-side currents and DC-link voltage in the “
abc” frame as follows:
In these equations and the associated
Figure 1,
is the output filter,
is the DC-link capacitor,
is the DC-link voltage,
is the power generated by the primary power source,
is the VSC terminal voltage,
is the grid (PCC) voltage, and
is the grid-side three phase currents. To remove the nonlinearity of (2), it is a common practice to use the DC-link energy
instead of the DC-link voltage
as a state variable. Using a phase-locked loop (PLL) based on
and Park transformation, (1) and (2) can be shown by (3) in the “dq” frame representation form:
where
,
and
are state variables;
−
, and
−
are the inverter and PCC voltages in the
d-
q frame, respectively. It can be seen that the DC-link voltage dynamic is still nonlinear, and (4) shows the dynamic equation, by substituting the grid-voltages instead of the terminal voltages:
where
is the instantaneous power of output filter, which is further expressed in (5):
The first part models the change in the energy stored in the inductance of the filter and the second term models the conduction and switching losses in the inverter and output filter. To consider the worse possible condition from a stability perspective, the loss parts are ignored. Substituting (5) into (4) and with some simplification, the dynamic equation is as given in (6):
In addition, considering the zero value of
in steady state due to the operation of PLL, the output active and reactive power of the inverter are as given in (7):
It should be pointed out that the nonlinear terms in the right-hand side of (6) complicate the control design for this system. In order to alleviate the complexity of the nonlinear DC-link voltage dynamic, either a simplified model without non-linear parts or a linearized model with non-linear parts can be used.
The simplified model is as shown in (8):
Accordingly, the closed form of the state equations are as follows:
In the linearized model, the nonlinear terms are linearized at the operating point (
) as follows:
According to (7) the operating point values of
and
are:
Substituting (3) and (12) into (11) leads to the following linearized model of DC-link voltage dynamic:
By defining
,
and simplifying (13), the linearized model of DC-link voltage is represented by (13):
Despite the higher accuracy, the system model of (14) is dependent on the operating point, and is more complicated compared to (8). Based on (14), the closed form of the state-space equations is as follows:
where,
It is worth mentioning that the model of (9) is the same as the model of (15), but in a special operating point, i.e., ().
3. Conventional Control of Grid-Following Inverters
The two main control purposes of a GFL-inverter are (i) delivery of the produced active power by the primary source to the grid and the exchange of the desired reactive power to the grid, according to the grid-codes, (ii) control of the DC-link voltage in the desired reference value. Both of the requirements should be met with a proper response time and stability margins.
According to (7), it is general practice to control the and with and , respectively. In addition, according to (6), control of DC-link voltage at the predetermined value drives to ; and, accordingly, the desired is delivered to the grid. Therefore, the control goals would be control of the reactive power and DC-link voltage, which are done via control of and . In this cascade control structure, the current control loops of and are inner loops, and the reactive power and DC-link voltage loops are outer control loops.
3.1. Current Control Loop
As shown by (3), there is a tight coupling between the
and
dynamics, and the grid voltages also affect the system dynamics. Hence, the proper decoupling actions and the feed-forward-based disturbance rejection approaches are required. Equation (15) provides the necessary control effort to cover both demands.
where,
and
are intermediate control variables, which are introduced later on in (19). Substituting (17) into (3) yields to the resulted current dynamics of (16):
By defining the intermediate control variables of (17), the current control loop transfer functions are as given by (18), in which
is the current control loop bandwidth (BW):
usually is chosen at least 10 times smaller than that of the switching frequency of VSC, to preserve the proper system stability in the presence of the time delay in the current control loop, imposed by computations and PWM switching in the control loop.
The current references of and are obtained from the reactive power and DC-link control loops, respectively.
3.2. Reactive Power Control Loop
The reactive power control loop determines
using the reactive current setpoint of
based on (21):
According to (7), using (21) and even without
,
is driven to
. However, uncertainty in the system, the delays in the current control loop, and the simplification in the system model may lead to a steady state error in controlling
. To guarantee a zero steady-state error, a PI controller is also used, which provides the
by (22):
The bandwidth of the reactive power control loop is chosen at least four times lower than that of the control loop.
3.3. DC-Link Control Loop
Considering the three dynamic models for the DC-Link presented in Equations (3), (8) and (14), three schemes can be used for DC-link control [
23], which are:
Using the simplified DC-Link model of (8) and designing a PI controller. This method is the prevalent control design scheme in the literature, known as “conventional DC-link control scheme”;
Using the operating point-dependent model of (14) and designing an adaptive DC-link controller;
Using the nonlinear dynamic model of (3) and designing a nonlinear DC-link controller for a complicated system.
Among the above options, using the simplified DC-link model of (8) results in a simplified control design. However, it imposes restrictions on the bandwidth of the DC-link control loop, especially in rectifying mode. The deficiency of the restricted bandwidth, and accordingly the poor disturbance rejection capability, can be solved by using a larger DC-link capacitor. As a result, this method is a commonly used approach in the literature for finding the “d” axis current reference. According to (8), the DC-link voltage dynamic depends on both
and
. To cancel the impact of
, it is considered as a measurable disturbance and compensated for by feedforward control. Finally,
is determined as given by (21):
It is worth mentioning that the bandwidth of the DC-link control loop is selected at least four times lower that of control loop. In addition, reducing the bandwidth improves the system stability, considering the uncertainty in the model of (8).
4. Proposed Observer and Control Structure
According to the control structure of the GFL-inverter, the grid-side currents, the grid voltages and the DC-link voltage are required for control of the inverter. By analyzing the system dynamic equation of (3), it can be seen that there is a coupling between , , and . Then, theoretically, it is possible to estimate the grid currents (, and ) by measuring the DC-link voltage. This is the key idea of this paper in developing an observer for the grid-side current estimation. In this respect, a well-adapted Luenberger observer is used in this paper, as the system is observable. In the following subsection, the fundamentals of the Luenberger observer are presented.
4.1. Fundamentals of the Luenberger Observer
Generally, the Luenberger observer is designed for linear systems with the given block diagram of
Figure 2. This figure consists of three main blocks of (i) the real system, (ii) a dynamic model of the system, and (iii) Matrix L. The goal of the observer is the accurate estimation of the state variables of the system. The best dynamic response and accurate estimation is achieved if the dynamic model of the system accurately models the real system. In other words, it is proved that by providing the accurate state-space dynamic equation of the system, the correct initial conditions, and the inputs of the real system into the simulated state-space equation, the response of the simulated system, including the output and state variables, are equal to those of the real system.
However, the model of the system may have some uncertainties. In addition, the initial values of the state variables in the real system are not available, and using them in the simulated system is not feasible. To compensate for this deficiency, the error between the measurable output of the real system and that of the simulated system is calculated and multiplied with a matrix, i.e., observer matrix L. The results are fed back to the simulated system, which finally reduces the error values. In an observable system, it is proved that all of the errors between the state variables of the real system and the simulated system are driven to zero by choosing a proper value for matrix L. Accordingly, the time-domain waveforms of the state variables of the simulated system become equal to those of the real system. By this procedure, the state variables of the real systems are obtained without using a direct measurement tool. The details of the Luenberger observer design procedure can be found in [
24], which is not further elaborated for the sake of brevity. The application of the Luenberger observer for the control of the GFL-inverter is provided in the following subsection.
4.2. Estimation of the Grid-Side Currents Using Luenberger Observer
In this subsection, the presented Luneburg observer is applied on a GFL inverter to estimate the grid-side currents. As mentioned before, this observer is appliable on linear system, and the linearized equations of the system are required for the design procedure. In this regard, considering the nonlinear dynamics of the DC-link control, the models of either (8) or (13) can be used in the observer-design procedure. In the grid current estimation, there are two alternatives of:
Using the system model in (9) for the observer design, which includes the simplified linear DC-link model of (8);
Using the system model in (15) for the observer design, which includes the linearized operating point-dependent DC-link model of (13).
As a rule of thumb, it is clear that using (9) leads to a simple design process; however, the approximation of (9) may lead to an undesirable response. On other side, (15) is an operating point-dependent model, and consequently, an adaptive observer, such as a gain-scheduling observer, should be designed. To make a compromise, the observer is designed for the middle operating point of
in this paper, which still can be used for all of the operating points. Therefore, it imposes less approximation compared to (9) on average. Based on the Luenberger observer design scheme, the state-space equation of the observer consists of (1) state-space equation of the main system and (2) a compensating term. Considering the GFL-inverter linearized model of (15) in the middle operating point
, the observer state-space equations can be written as (24):
In this equation, the last term is the compensating term; other terms are a model of the system at , the matrix is the observer matrix with three elements, is the estate vector, and is the estimated state vector. In , and , are, respectively, the estimated values of the grid current in d and q axes, and is the estimated value of the DC-link energy. As the DC-link energy is directly measured, its estimation is not used in the control structure. Inputs of the observer are , , , (if available), and the control variables of and . Based on the Luenberger observer design scheme, the observer matrix must be determined in such a way that all of the eigenvalues of the matrix become stable.
It is worth noting that the first proposed observer in (24) is a third-order system. Considering the order of the system, which is three, and the measured value of
, it is possible to estimate the
and
by a second-order observer (reduced order observer). This may lead to a lower computational burden in the control loop. However, to include the nonlinear DC-link model in the observer implementation, it is necessary to use the third-order observer, which considerably increases the performance. Moreover, although matrix
and the compensating term of the observer are designed based on a linearized model of (24), the exact nonlinear model of the system can be used in the final form of the observer. Accordingly, the observer of (24) is replaced by (25), which includes the nonlinear DC-link dynamics of the system:
In summary, first, a constant observer matrix is designed, based on the system model of (15) in middle operating point , and then it is used in the compensating term beside the exact nonlinear model of the system to complete the observer dynamic, as in (25).
Figure 3 shows the block diagram of the overall system, including the controller, observer and the GFL-inverter power circuit. It can be seen that only the grid voltages and the DC-link voltage are measured, and using the proposed scheme, the estimated values of the grid currents (
,
) are used as feedback signals in the control loop, instead of measuring (
,
). In this structure, the three current sensors of the GFL-inverter are eliminated, and the terminal voltages description of the inverter will be as shown by (26):
5. Stability Analysis of the GF-VSC with the Proposed Observer
The stable operation with a high-enough stability margin is the preliminary requirement of a given control system. In this section, the stability of the GFL-inverter with the proposed observer and control structure is analyzed, using the eigenvalue theory. To cover all of the dynamics of the system, the dynamic model of the grid-side current and the DC-link voltage of (14), the observer dynamic of (25), and the controllers of (21) and (23) are integrated and used to obtain the overall system model. It should be mentioned that, for stability analysis by eigen values, a linearized model of the overall system is needed. A linearized state-space model of the primary system was presented in (15), also, the linearized model of the observer in (25) is given by (27):
Therefore, it is only needed to extract the state-space model of controllers, to complete the state-space model of the overall system shown in
Figure 3. According to the block diagram in
Figure 3, the two state variables of
,
are defined to represent the controller in the state-space equations as follows:
Furthermore, based on the block diagram of the system in
Figure 3, the control variables of
and
are as given in (29):
In which, the current references of
and
are given in (30):
Substituting (30) into (29), and with some simplification, the result is as given in (31):
where,
It is worth noting that the inputs of the system
and
do not affect the stability of the system, and they are eliminated in (31) for the stability analysis. Applying the control variables of (31) into the linearized model of the system and observer in (15) and (27), then by merging with the controller-state equations of (28), the resultant linearized model of the overall system is given in (33):
The closed form of (31) is written in (34):
The stability of the overall control system is dependent on the eigenvalues of the state matrix of (34). Given that and are dependent on the operating points of GFL inverter (), the eigenvalues of the state matrix of (34) should be analyzed in a full range of operating points. In addition, the eigenvalues of the conventional control system without an observer are presented for comparison purposes.
In the case study and numerical analysis, a sample GFL-inverter, with the characteristic parameters given in
Table 1, is used in the following.
Figure 4 shows the poles of the GFL inverter with a conventional control system without using an observer in different operating points. From the figure, it can be seen that the system is a fifth-order system, and increasing the
would shift the poles to the more negative values. In addition,
does not affect the poles, and hence, the poles are shown only for
. Furthermore, all of the poles have a damping coefficient of
regardless of the operating points.
In the design of matrix
, it is possible to place the observer poles of
at any of the desired values. The negative poles with higher values provide a higher bandwidth and increased observer speed, which consequently decreases the estimation error in the presence of uncertainty in the model of the system. However, it may deteriorate the stability of the closed loop system in the presence of uncertainty. To further evaluate the effect of the observer poles on the stability of the overall system in various operating points,
are illustrated in
Figure 5. Four cases for the poles of the observer are examined. Suppose that the four-observer matrix is designed in such a way that the poles of the observer are as follows.
In this study, the position of the poles ([
,
,
]) are chosen depending on the current control-loop bandwidth (
), and
, and
are considered similar to each other.
Figure 5 shows the eigenvalues of the closed loop system for the four observers of Observer#1-to- Observer#4 in the six operating points under the full range of active power and reactive power production, as given in (36):
From
Figure 5, it can be seen that the overall system is stable for all four observers in all six operating points. In addition, increasing the observer bandwidth moves the poles to the far away points in the s-plane, but as expected, decreases the poles damping coefficient
. This is due to the difference between the model of (24), used for the design of matrix
L, and the exact model of system in (3). Moreover, it decreases the delay margin of the system, which, considering the time delay imposed in the control loop by VSI, may destabilize the system. On the other hand, reducing the observer bandwidth may increase the estimation error in the presence of the uncertainty. Accordingly, a compromise is needed in the selection of the observer bandwidth and the system stability. In this regard, the third observer, whose corresponding matrix
L is presented in
Table 1, is chosen as the final observer.
6. Case Studies, Numerical Results, and Discussions
The performance of the proposed observer and controller is evaluated through the various case studies on the sample test system in this section. The complete characteristics of the system are presented in
Table 1.
To evaluate the performance of the system under different operating conditions, different values and changes are applied to , , and Both the actual and estimated grid-side currents, the DC-link voltage, and the active and reactive powers are presented. In addition, the system performance in the presence of uncertainty in the output filter inductance is examined. The detail of the scenarios is as follows.
Initially, the values are set at
10 kW,
, and
. A step-up and a step-down command are given to
, respectively. Then, a step-up and a step-down command are applied to
, respectively. Furthermore,
is gradually changed from 10 kW to 0 kW. While
is 0 kW, another step-up and a step-down command is given to
and
. The changes in the reference values (
,
, and
) are shown in
Figure 6. To make the simulation tests close to the practical experimental tests, zero order hold sampling, PWM and a computation delay of
are also intended in the simulations. The simulation results are presented in
Figure 6 and
Figure 7.
Figure 6 shows that using the proposed control scheme, including the proposed observer, DC-link voltage and
of the GFL inverter, properly tracks their reference values with a fast dynamic response and a high stability margin. In other words,
Figure 7 shows that using the observer in the control loop of the GFL inverter instead of the grid-side current measurement, does not affect the dynamic performance of the overall GFL inverter. Also, as shown in
Figure 7, the estimated currents
and
are driven to their actual values very fast in startup and under the reference values variation. In addition, for better vision, the zoomed view of deriving
and
to their actual values in the startup time are shown in
Figure 6c,d, which demonstrates the proper performance of the proposed observer.
Moreover, to show the robustness under the presence of uncertainties, a 20% uncertainty in the output filter is considered. It is supposed that the actual
is 20% larger than the value used in the design of the observer and the control loop. In addition, as another scenario, it is supposed that the actual
is 20% smaller than the value used in the design of the observer and the control loop. The simulation results in
Figure 8 show that such large uncertainties in the
value partially affect the system performance. It negligibly decreases the speed of the system in tracking the reference values of reactive power and DC-link voltage. Such minor effects can be ignored, since the important features of the system, including the tracking capability, dynamic performance, high functionality and stability margin of the system, are preserved even under the presence of the uncertainty in the model of the system.
The different case studies and scenarios and the obtained results demonstrate that the proposed control scheme effectively provides proper performance even with a reduced number of measurement signals. This provides the manufacturers of the GFL inverter with an opportunity to reduce the cost and increase the reliability of the inverters by reducing the number of required sensors.
7. Conclusions
In this paper, the control of a GF-inverter with an L-type filter is studied and the idea of estimating the output filter current was proposed. In the proposed idea, the output filter currents were estimated with an observer-based controller, which eliminates the need to use the grid-side current sensors. Using a lower number of sensors, the total cost of the system is reduced, and the reliability and fault-tolerant capability of the GF inverter are increased. To develop the analytical formulation of the proposed idea, modeling of the GF inverter with a focus on the DC-link voltage dynamic model is presented. It is shown that the output filter currents are observable and can be estimated by measuring the DC-link voltage and employing a linear observer. In the design of the observer matrix, the linearized model of the non-linear DC-link voltage dynamic is used, which provides a simple constant observer with proper accuracy. To increase the performance in estimating the real system dynamics, the nonlinear model of DC-link voltage is employed in the observer implementation. Moreover, the stability of overall system is evaluated by using the system state-space equations in different operating points through eigenvalue analysis. Various time domain simulation in MATLAB/SIMULINK and case studies verify the proper suitable tracking, estimation functionality, dynamic performance, and stability margin of the proposed control system in all operating pints even under presence of uncertainty in output filter inductance.