1. Introduction
Inverter-interfaced distributed energy resources (DERs) are increasingly common in three-phase distribution networks, forming microgrids at medium or low voltage. The microgrid concept offers numerous benefits, including the capacity to integrate renewable energy sources near the loads, the integration of energy storage devices, and the possibility of islanded operation [
1]. The DER units in a microgrid also can provide support for a host of power quality issues. For example, when the microgrid operates in a grid-connected mode, the sharing of reactive power by the DERs contributes to the local voltage regulation [
2]. Among the challenges in power quality management, the compensation of imbalance is a significant one where DERs can make a substantial contribution. Unbalanced loads are characteristic of low-voltage networks and demand asymmetric currents that, in turn, can considerably affect the symmetry of the microgrid’s voltage. This may be exacerbated during asymmetric faults. In such cases, the DER units are required to remain in operation, at least temporarily, compensating some of the imbalance [
3]. In [
4], a current reference strategy including an explicit negative sequence component is suggested to contribute to voltage equalization during asymmetric fault ride-through. The compensation of imbalance is particularly important in the islanded mode of operation [
5,
6,
7,
8] where the DERs must supply the whole of the imbalance currents.
The implementation of appropriate control algorithms for unbalanced voltage support in unbalanced grids requires the decoupling of the so-called positive and negative sequences of asymmetrical three-phase voltages and currents. Typically, the voltage sequence components are used for current reference generation, and the current sequence components are used directly for feedback in the inner current control loop; see, e.g., [
8]. Thus, the dynamic response of the current sequence separation algorithm influences the system’s response and should be considered in the control design for optimal performance. However, this aspect is, to the best of our knowledge, overlooked in the literature.
The concept of positive and negative sequence components in asymmetric three-phase systems was introduced by Fortescue [
9] more than one hundred years ago, and stems from the fact that real sinusoidal signals can be expressed in terms of two complex exponentials of positive and negative frequencies, i.e.,
and
via Fourier series. Hence, one of the main characteristics of the Fortescue transformation is that it only applies to the steady-state regime. Afterwards, Clarke [
10] introduced the idea of instantaneous complex vector transformation applicable to time-domain signals in transient regimes. Although this transformation does not include an explicit sequence separation, both Fortescue and Clarke transformations are strongly interconnected [
11,
12]; the former is applicable to signals in phasor representation, whereas the latter in the time domain. As a matter of fact, the Clarke transformation [
10] yields a representation of the three three-phase signals in a coordinate system aligned with the direction of its average (or common mode) and a plane orthogonal to this common mode, the so-called
plane [
13]. The electrical quantities of interest can then be represented as a vector,
, or in complex representation
. In this plane, the pure positive and negative sequences of a stationary three-phase system would be fixed magnitude complex signals that rotate in opposite directions at a fixed (synchronous) frequency. However, in this representation, both positive and negative sequences are condensed in a single quantity. The Clarke transformation leads to a more general representation of arbitrary signals that may include transients, harmonics, and noise. Unfortunately, there is no bijective transformation that separates the positive and negative sequence components from the resulting complex signal. Therefore, the positive and negative sequences can only be estimated by a dynamic algorithm, such as filters or observers.
Possibly the best-known algorithm for sequence separation is the delayed signal cancellation method [
14]. This method uses a linear combination of the complex representation of the electrical quantities of interest
with a version of itself delayed by a quarter of its fundamental period (
). The only dynamic element on this filter is delay, and the rest are algebraic calculations. Therefore, this algorithm does not perform well with distorted signals or during transients. A family of methods that do incorporate filtering in the classical sense stems from the synchronous reference frame phase-locked loop (SRF-PLL) of [
15]. To better handle unbalance, the use of two synchronous rotating frames (in opposite directions) and a decoupling network has been proposed in the double synchronous reference frame PLL (DSRF-PLL) [
16]. The decoupled double synchronous reference frame (DDSRF) algorithm can be used for sequence separation [
4]. Another family of synchronization methods is based on the second-order generalized integrator (SOGI-PLL) [
17], originally devised to make use of the advantages of complex signal processing on single-phase (real) signals. The second-order generalized integrator is an adaptive band-pass filter centered at the fundamental frequency, providing additional filtering with respect to SRF-based methods and improving its handling of distorted signals. By applying the SOGI method to both components
in an unbalanced system and using a decoupling network, the double second-order generalized integrator PLL (DSOGI-PLL) [
18] and its frequency-locked loop variant DSOGI-FLL [
19] can separate the positive and negative sequences effectively. For these reasons, the DSOGI has become a well-established method for sequence separation and synchronization in distorted and weak grids [
20] and can be used as a benchmark for other methods. However, all the DDSRF- and DSOGI-based methods rely on coordinate rotations, filters, and linear combinations of the signals, resulting in deterministic implementations. The random nature of the sequence separation problem has been acknowledged in [
21], which has led to the utilization of stochastic estimation techniques. In [
21], a complex-valued Kalman filter is used to estimate the positive sequence and negative sequences. Moreover, the order of the model is expanded to estimate harmonic components. In [
22,
23], the problem of sequence separation is tackled together with the frequency estimation problem by implementing an extended Kalman filter (EKF). The idea of utilizing a complex-valued model was further explored in [
24], where an extended model was developed and the estimation was carried out using a complex-valued EKF. However, the EKF approach in real-time closed-loop applications presents a high computational load and results in an uncertain dynamic response.
In many of the techniques described above, the real and imaginary components of the complex variables that arise from the transformations are treated separately, favoring vector
(real) representation. This approach disregards the compactness and analytical advantages of dealing directly with complex-valued variables. The processing of complex signals (
) has many advantages over the processing of real signals and this fact is used extensively in many sub-fields of electrical and electronic engineering, for example, in signal processing and communications (see, e.g., [
25] and the references therein) and in drive control [
26]. On the other hand, the advantage of using complex signals in sequence separation has been recognized in [
21], where the real and imaginary parts of the process noise and measurement noise are assumed mutually independent. In this paper, the solution of the sequence estimation by a complex time-varying Kalman filter (KF) proposed in [
21] is reviewed in detail in terms of the statistical properties of the process noise and measurement noise allowing for a straightforward implementation of the complex-valued Kalman filter, in which the transpose operator is simply replaced by the hermitian (conjugated transpose) operator. Moreover, in this paper, a linear transformation of the state variable is proposed, resulting in a stationary solution of the Kalman filter [
27]. This complex-valued stationary KF for sequence separation is the main contribution of the article and leads to a sequence estimation algorithm with similar performance to the conventional (real valued) time-varying KF, but at a fraction of the computational cost. Additionally, this stationary filter results in a constant and predictable dynamic performance, facilitating its incorporation into the control design. The benefits of this property for modern control design are illustrated in a numerical example.
3. Modeling and Estimation of the Complex Positive and Negative Sequence Components
The asymmetrical three-phase signals have been expressed as a complex variable in a
reference frame in (
6). Here,
is the observed or output variable, and
,
are its positive and negative sequence components to be estimated, i.e., the state of the system to be estimated. In other words, (
6) can be thought of as the output equation of a dynamic system. Unfortunately, the term
provides a time-varying characteristic to the dynamic model. Another difficulty to overcome is that the sequence components, i.e., the state of the system, depend on unknown conditions, like load imbalance, for example. Therefore, it may be concluded that an accurate dynamic model is not readily available for the system of interest. However, it can be safely assumed that
is a constant variable at steady state. Hence, its dynamic behavior may be modeled as an uncorrelated Gaussian random walk:
with
, and
a zero mean circularly symmetric (i.e., proper) complex Gaussian random variable with variance
,
.
In the following, estimation algorithms that consider both the simplified dynamic model (
9) and its time-varying output Equation (
6) are proposed. Since the original problem is formulated in the complex plane, the use of the complex-valued Kalman filter (CKF) [
21,
29] is also considered in this section as a more compact alternative. Finally, based on the representation of the system in complex form, a linear transformation is used, assigning the time-varying term to the actual signal to be estimated, resulting in a time-invariant model. By applying this transformation, the stationary complex-valued Kalman filter (SCKF) is deduced. This implementation greatly reduces the computational load of the sequence separation without loss of accuracy. These two algorithms are detailed next.
3.1. Complex-Valued Kalman Filter (CKF)
The dynamic model (
9) and the output Equation (
6) correspond to a complex-valued state-space formulation. This representation can be used, for instance, in the formulation of a time-varying linear quadratic estimator problem, with the well-known time-varying KF as the optimal solution [
30]. For complex linear systems, the CKF can be formulated using the well-known structure of the real-valued KF [
31] when the initial condition. The process noise and the measurement noise are all proper complex random variables. In this case, the transpose operator from the real-valued KF is simply replaced by the hermitian (conjugate transpose) operator; see e.g., [
29]. More details about the derivation of the KF and the complex-valued KF are provided in
Appendix A and
Appendix B, respectively.
For implementation, the output Equation (
6) can be written as
With the state and output equations in (
9) and (
10), respectively, the corresponding CKF is obtained:
In the algorithm, the time-variant part of , can be assumed to be known as the grid frequency is known by the local controllers, and the phase of the frame is arbitrary.
3.2. Stationary Complex Kalman Filter (SCKF)
In this section, the main results of the article are presented. Here, a change in variables to the positive and negative sequence in state vector
is introduced, given by (
9) and (
10), that yields a time-invariant model. For estimation purposes, this model simplifies the implementation and reduces the computational load of the associated Kalman filter. To that end, a new complex-valued random variable
is estimated instead of
. This new variable is defined as
By inserting (
12) to output Equation (
10), the new output Equation (
13) is obtained:
Similarly, transformation (
12) can be applied to the time-evolution of the variables to be estimated in (
9), yielding
As mentioned before, variable
is known for each time instant
k and its time domain model is given by
. Therefore, the time evolution of transformation
can be written as
and the time evolution of
is given by
where
. Again, since
is circularly symmetric and transformation
is a rotation in the complex plane, the vector of uncertainties preserves the original properties
.
Finally, the CKF algorithm can be evaluated with state and output Equations (
16) and (
13), respectively, resulting in
The matrices
and
are time-invariant and the pair
is detectable; therefore, the convergence of the complex Riccati equation associated to the estimation problem is ensured. Consequently, one can implement the stationary solution of the CKF (SCKF) without loss of estimation accuracy and at low computational complexity. Here, the time index “
k” is replaced by the symbol “∞” to denote the quantities that converge (i.e., the filtered covariance and the Kalman gain) for the stationary filter. In order to obtain the corresponding SCKF algorithm, the following discrete-time complex algebraic Riccati equation should be solved:
To solve the algebraic Riccati equation, numerical packages are widely available. One alternative that can handle complex entries is MATLAB R2023b (The MathWorks, Inc., Natick, MA, USA), using command
idare(·). Then, given the numerical value of
, the stationary complex-valued Kalman gain can be found evaluating the following:
With the stationary complex-valued Kalman gain obtained, the implementation of the SCKF for
is straightforward, and can be written as
This is a recursive time-invariant algorithm with constant matrices
. Finally, solving (
20) for each time instant
k, the complex vector
is found and the estimated positive and negative sequences are obtained by evaluating (
12) for the estimated values. That is,
For the sequence estimation problem, transformation (
12) is a simple and effective way to express the time-varying nature of the problem. As a consequence, the estimation problem results to be time-invariant in the new state space. Hence, the stationary solution of the complex-valued Kalman filter can be obtained, simplifying the implementation and significantly reducing computational cost, without loss of accuracy. The straightforward SCKF for the sequence separation problem is expressed in Algorithm 1.
Algorithm 1: Proposed SCKF Algorithm. |
|
4. Simulation Results
To evaluate the performance of the proposed SCKF and to compare it to the performance obtained by the well-established DSOGI method [
18,
19,
20], an open-loop sequence separation problem is simulated. A brief explanation of the DSOGI algorithm is included in
Appendix C, and the interested reader can find further details in [
20].
The simulated system for this test is depicted in
Figure 1. It consists of a symmetrical three-phase voltage source of 111 V phase—neutral, an islanding transformer connected in Delta–Delta, with turn ratio
, and a three-phase resistive load with
per phase. This load is intentionally unbalanced to produce asymmetrical three-phase current waveforms by disconnecting load phase
b.
For the implementation of the SCKF estimator, the values of the state and measurement noise covariance matrices,
and
are required. These serve as tuning parameters for the simulation. In the interest of fairness, the same values used for the experimental testing discussed in the next section are also used in the simulation. These values can be found in
Table 1. Noise covariance matrices were selected by trial and error. On the other hand, the DSOGI algorithm only has one tuning parameter, and value gain
used in [
18,
19,
20] is considered.
A thousand Monte Carlo simulations were carried out, which is equivalent to measuring the same event a thousand times.
Figure 2 shows the mean of the Monte Carlo simulation results. To examine and compare the estimation performance of each method, the separated sequence currents are presented in positive and negative rotating
frames, respectively. These results show that, although both algorithms are efective in achieving sequence separation, the SCKF has a faster response than DSOGI.
5. Experimental Validation
To confirm the aforementioned results, experimental tests are run in real time to evaluate the performance of the SCKF proposed in
Section 3.2. Two different case studies are considered to validate the effectiveness of our proposal: an open-loop sequence separation problem and a closed-loop current control application. In the first case, a three-phase resistive load connected to the grid is severely unbalanced, resulting in a fast transition from symmetric to asymmetric currents replicating experimentally the simulation of the previous section using a balanced star-connected resistive load of
, 1 kW per phase and a circuit breaker in phase
b to produce the imbalance. A resistive load is chosen because it demands the fastest change in load current, making it the most stringent test for assessing the performance of sequence separation algorithms. In this case, the grid is physically emulated by a three-phase
Chroma 61704 AC power supply (Chroma ATE Inc., Taoyuan, Taiwan) that produces very clean sinusoidal signals and an isolation transformer connected in Delta–Delta with turn ratio
. The setup is instrumented with LEM current transducers, and the current measurement signals are sampled at
s by a dSPACE 1103 (dSPACE GmbH, Paderborn, Germanny) that performs all the calculations in real time. A general-purpose I/O signal of the dSPACE is used to control the circuit breaker via a chain of relays. For this experimental case study, the two discussed complex versions of the recursive algorithm (CKF and SCKF) are used to separate the positive and negative sequences of the currents to contrast their estimation performance and also their computational burden. In addition, the equivalent real-valued KF is used for comparison purposes. The description of the real-valued KF can be found in
Appendix A. In the second case, the estimation of the SCKF-based algorithm is used as feedback in a positive and negative sequence control scheme of an
-filtered DER unit, and the control performance is compared with the direct feedback of the asymmetrical variables. The DER consists of a 6 kW Semikron laboratory 2-level voltage source inverter (2L-VSC) (Semikron GmbH, Nuremberg, Germany) properly interfaced to be modulated by external signals. The primary energy resource (ER) is emulated using a programmable DC power supply (model N8762A of Agilent Technologies, Inc., Santa Clara, CA, USA), set at
. Finally, the same isolating transformer and three-phase load as in the previous case study are used. For this experimental case study, the estimation, the control and the space-vector PWM modulation are implemented in a dSPACE 1103 with a sampling time of
s, which is also the period of the PWM carrier for close loop control.
Figure 3 shows a picture of this experimental setup.
5.1. Case Study 1: Positive and Negative Sequence Estimation of an Asymmetrical Signal
This case study evaluates the estimation performance of the estimation algorithms described in
Section 3. In particular, it is of interest to demonstrate that the proposed time-varying CKF exhibits performance similar to that of the conventional Kalman filter when the latter is implemented with an extended state vector of real variables (see
Appendix A for details). These results help to validate the assumption of circular symmetric noises underlying the CKF formulation. Furthermore, it is of interest to show that the proposed SCKF exhibits almost identical estimates to the time-varying CKF after the initial transient of the Kalman gain matrix passes. As a result, the SCKF enables comparable performance with a significantly reduced computational load of the estimator. The experimental setup corresponds to the circuit depicted in
Figure 1 and its parameters are the same as in the simulations of the previous section. In this experiment, the load is intentionally unbalanced to produce asymmetrical three-phase current waveforms by opening a single-phase breaker connected in phase
b.
The implementation of CKF and SCKF estimators requires the knowledge of state and measurement noise covariance matrices,
and
. Additionally, these algorithms require to be initialized with the expected value of the initial state and covariance matrices at time instant
. The values used in each algorithm are given in
Table 1. Notice that for the stationary filter, i.e., SCKF, the stationary Kalman gain
can be directly computed from matrices
,
and the matrices that define the dynamic system (
and
, see (
18) and (
19)). Therefore,
does not depend on
. To examine and compare the estimation performance of each method, the separated sequence currents are presented in positive and negative rotating
frames, oriented with the supply voltage as
, respectively.
The experimental results are shown in
Figure 4. The three-phase load is balanced, with Br
1 closed, hence demanding symmetric currents. The estimation is carried out using the three algorithms at
s with the same initial value (either in real or complex form). For KF and CKF algorithms, the estimated states are almost identical, validating the assumption of circular symmetric complex noises in the CKF formulation. On the other hand, the dynamic response of the proposed SCKF differs slightly for the first few samples due to the difference in the initial Kalman gain of the time-varying algorithms (KF and CKF) with respect to the steady state
used in the SCKF. The three algorithms converge to the same Kalman gains after
ms, making their time responses almost indistinguishable. At
s, the breaker Br
1 is opened and the load in phase-
b is disconnected, resulting in a considerable imbalance. As shown in
Figure 4, the three implementations of the KF yield almost identical state estimates during this large transient, converging to the same steady-state values within (approximately) a third of the fundamental period. Therefore, the use of the proposed SCKF for sequence estimation can effectively replace the time-varying algorithms, KF and CKF. Consequently, positive and negative sequence estimation can be easily and efficiently implemented without loss of performance. In terms of computational complexity, first, the structure of the problem of interest is considered. That is, it is taking into account that some elements are equal to 1 when involved in multiplications, that some elements are real when involved in additions and multiplications, and that some matrices are symmetric or Hermitian when computed. In addition, it is considered that a real sum and a real multiplication are both equivalent to one floating point operation (FLOP), a complex sum equals 2 FLOPs and a complex multiplication equals 6 FLOPs. Then, the real-valued KF (see [
21]) results in 213 FLOPs per iteration, the complex-valued KF in 74 FLOPs, whilst the proposed stationary CKF in only 50 FLOPs, which corresponds to a reduction of
of the complexity with respect to the real-valued KF and a reduction of
of the complexity with respect to the complex-valued KF.
5.2. Case Study 2: The Use of the Estimated Sequences on a Closed-Loop Strategy
This case study explores the application of the proposed SCKF in a closed-loop control strategy that uses the estimated positive and negative sequences as feedback variables. The system under study is shown in
Figure 5, where the
filtered stand-alone DER is subjected to an unbalanced load disturbance. The main goal of this case study is to confirm that the utilization of the state estimations (i.e., positive and negative sequences) that are obtained from the proposed SCKF in a real closed-loop application can yield good control performance in the presence of unbalanced loads. For comparison purposes, a closed-loop strategy based on the direct feedback of the asymmetrical measured variables, and otherwise similar control strategy, is used. The experimental set-up in
Figure 5 has the following parameters:
,
and
. A second output filter inductance is provided by the leakage inductance of the transformer, approximately
, but for control purposes, it can be considered part of the load of the system.
The two control strategies that are compared in this case study are based on state feedback. These are briefly described in the following paragraph.
The first control strategy uses the estimated positive and negative sequences from the proposed SCKF as feedback, and it can be summarized by the control diagram in
Figure 6. In this scheme, the proposed SCKF algorithm follows the design presented in
Table 1. In
Figure 6,
and
are the positive and negative sequence state variables,
and
are the positive and negative sequences of voltage
modulated by the VSC. Currents
and
are demanded by the load and are considered as disturbances for the system. Finally,
,
are the voltages at the shunt capacitor which are the output variables to be controlled. On the other hand, a direct or naive control strategy, without the sequence separation, is illustrated in
Figure 7. Here, the state is given by
, the control action is
, the disturbance is given by load current
and the output to be controlled is
.
Both control strategies follow a scheme based on state feedback that includes the set point of the variable to be controlled (
) and the compensation of the measured disturbance (
). It is worth noting that in this case study, where the capacitor voltage is the controlled variable, orientation angle
is directly obtained using the integration of frequency, and its phase is automatically adjusted by controlling the imaginary component
to a zero value. Feedback gains
are designed as a linear quadratic regulator, and the matrices that map the set point (
) and disturbance (
) to the input variable are obtained following the analytical results in [
27,
32]. For the separated sequence control, the gain matrices differ only in the signs of some of its entries as summarized in
Table 2. For direct control, the gain matrices are given by
,
and
. Note that this controller, shown in the scheme of
Figure 7, is similar to the control of the positive sequence in the scheme of
Figure 6, but the asymmetrical signal is used as feedback. Consequently, it is expected that the main difference between the closed-loop performances of both control methods would be manifested in the negative sequence of the variable to be controlled,
.
The experimental response of the closed-loop current control in response to a step change from balanced to unbalanced load is shown in
Figure 8, where the magnitude of the current and voltage vectors,
and
, and the actual currents and voltages in
coordinates are shown for both methods. The set-ups of
Figure 6 and
Figure 7 start with the balanced three-phase load connected, and the voltages references are
and
. At
s, the phase
b load is disconnected, resulting in an asymmetrical output current,
. Consequently, the DER unit’s output filter is subject to asymmetrical disturbance
that should be properly compensated at input
. As it can be seen in
Figure 8, the effect of load imbalance appears as an oscillation at the double of the fundamental frequency in the
frame for both methods, but the benefit of actively regulating both sequences is evident in the results shown in
Figure 8f. The impact of the imbalance on the controlled variable,
, is significantly reduced in this case. This confirms that the proposed SCKF can be used for sequence separation in closed-loop control, reducing the impact of the load imbalance on control performance.
The second closed-loop experimental result shows that by using the feedback of both estimated sequences, as shown in
Figure 6, the negative sequence of the voltage can be regulated with independent reference
. This is of particular importance in applications where a DER unit is required to share the imbalance of the loads. The result is shown in
Figure 9, where the initial condition corresponds to the same steady state as that at the end of the previous result, i.e., with
and
and unbalanced load. At
s, the reference of the negative voltage sequence is changed to
and an intentional imbalance at the voltage reference is introduced. As it can be seen in the measurement of
in
Figure 9d, both voltage sequences follow the change on demand with good dynamic response and high level of decoupling, showing that the proposed SCKF is useful for separate sequence control. It is important to note that the response exhibits a steady-state tracking error, which is an inherent characteristic of the proportional control utilized for testing and not a consequence of the proposed sequence separation algorithm.
In the above experimental results, the output filter parameters are assumed to be accurately known and the nonlinearities of the inverter are neglected. Those assumptions may be questionable for applications requiring high-performance closed-loop control. However, one significant advantage of the proposed SCKF is its predictable frequency response, which can be considered in the control design process. This allows for a systematic way of dealing with the above-mentioned modeling errors, for example, limiting the closed-loop bandwidth, performing sensitivity analysis or including integral action in the controller. Nevertheless, such control-tuning procedures are out of the scope of this paper.