1. Introduction
Renewable energies have become a replacement for conventional means of electrical power generation [
1,
2,
3]. Specifically, photovoltaic systems are of particular interest and therefore have been widely studied [
1]. The incorporation of renewable energy sources (RES) has helped to distribute clean energy and increase the proliferation of microgrids with everlasting power when compared with traditional fuel ones [
4].
Due to the low inertia of microgrids, they are greatly affected by power unbalancing between generation and consumption, leading to amplitude and frequency fluctuations [
5,
6]. Frequency variations—up to 100%—have been reported in applications such as aircraft, electrical ships, ferries, vessels, and seaports, among others [
2,
4]. This is the reason power converters need to employ algorithms capable of addressing these conditions and, subsequently, be able to integrate with systems such as those aforementioned.
To inject power energy into the AC grid, two approaches can be employed: (i) include a transition DC–DC converter and then a DC–AC inverter or (ii) only include a DC–AC stage, which reduces losses since one stage is removed [
7]. In other words, maximum power extraction and injection can be carried out with a direct DC/AC converter topology [
8,
9]. Therefore, the DC-link voltage is managed directly by the maximum power point tracking (MPPT) algorithm, e.g., a P&O algorithm [
8,
10,
11,
12]. This means a simplification in the topology used for power conversion since fewer elements and less physical space is used [
13].
In addition to the maximum power point tracking, power injection into the grid also requires synchronization algorithms, such as Phase-Locked Loop (PLL), to accurately track the grid frequency and phase shifting, as shown in references [
14,
15,
16,
17]. In this paper, the proposed control strategy employs a PLL algorithm, similar to that presented in [
14], in order to allow the converter to operate effectively under a wide range of frequencies, such as the conditions of microgrids herein employed.
Frequency and amplitude variations are part of the operation of microgrids and weak-grid systems [
14]; therefore, power converters connected to these kinds of networks require controllers capable of avoiding a fall into overmodulation or even instability [
4]. Research conducted in [
18] has demonstrated that variations in microgrid voltage can result in variations in the operating region. Specifically, when the frequency increases, the operating region is reduced due to an increase in passive filter impedance. On the other hand, under voltage swell conditions, the operating region is also decreased due to a reduction in the voltage difference between the injected voltage (
voabc, due to the DC voltage) and the grid voltage
vsabc, as described in [
19,
20]. It is crucial to note that power converters connected to microgrids may fall out of the operating region, and as such, power converter control must be designed to work under these conditions.
Due to the converter operating region studied, some of the latest control strategies propose a limit control on the Active/Reactive Power (PQ) plane operating region under unfavorable conditions [
20,
21]. These papers propose a way to overcome this issue by modifying the operating region in order to avoid overmodulation, which is the result of working out of the operating region.
As previously mentioned, the operating region can be reduced due to the operating conditions of a microgrid; for example, an increase in voltage frequency and/or swell amplitude variation. Therefore, in this work, a control method to overcome the operating region issue is presented; specifically, overmodulation is overcome by including the modification of the power factor and increasing the DC-link voltage whilst trying to maintain solar energy injection as high as possible. The proposed algorithm is designed for a single-stage DC/AC converter, and thus, increasing the DC voltage requires stepping away from the maximum power point (MPP) and should therefore be avoided as much as possible and only employed when necessary to overcome overmodulation, which may lead to instability. As a result, an overmodulation protection algorithm is proposed that allows the converter to continue operating over a wide frequency range and sag/swell conditions, at the cost of briefly losing the MPP of the photovoltaic system and the unit power factor.
The proposed control strategy is simulated and applied to a photovoltaic system using PSim© software. Finally, several experimental tests are carried out to validate the model by means of a laboratory prototype setup and the use of a TMS320F28335 digital signal controller (DSC) board from Texas Instruments to visualize some of the internal control variables through the graphics tool of the Code Composer Studio software. The most relevant conclusions of the work are also included.
2. Power Converter Model
A PV system connected to the point of common coupling (PCC) using a Voltage Source Inverter (VSI) configuration is shown in
Figure 1. It can be mathematically modeled through the Kirchhoff laws.
Employing the SPWM modulation technique, with injected voltage
voabc and injected current
isdc, can be defined as:
where
Gac represents the modulation technique gain and
mabc, the modulation function for each phase.
Analyzing the circuit in
Figure 1, the voltage equation on the AC side results in:
On the other hand, on the DC side, it can be found that
idc =
ic +
isdc, and therefore the DC-link capacitor model including (1) results in:
Thus, Equations (2) and (3) model the power system dynamically. However, to work with continuous variables, the
dq Park Transformation can be employed. Then, over the continuous variables system, it is possible to apply classical control theory to manage the currents and voltages. Consequently, by applying the Park transform to Equations (2) and (3), it is obtained:
where
W =
, ω = 2π
f, and
f is the grid frequency.
Therefore, the nonlinear dynamic model of the system in
Figure 1 is depicted in Equations (4) and (5).
3. Operation Region
In most grid-tied applications, the voltage maintains constant parameters (frequency and amplitude) and, consequently, the control is designed under these circumstances. However, in microgrids, the voltage frequency and amplitude may vary considerably due to their small electric system inertia [
22]. These variations affect the converters and their respective operating regions [
23].
A study of the operating region allows us to obtain the set of operation points where no overmodulation is reached. In addition, it helps to mark off the zone where the power converter engages overmodulation and, therefore, the model in (4) and (5) is not valid. Such a region is computed when the converter is in steady state, i.e., the derivations in (4) and (5) are equal to zero because the model is in
dq reference frame. Thus, the operating region is found from the following equations:
In addition, the equation that relates the power factor (
pf) and the currents in the synchronous axis can be stated as:
Subsequently, it is possible determine a relationship for
mdq, as a function of the state variables and system parameters, as:
where the
q-axis current is:
Then, to obtain the operating region, a sweep for values of
vdc and
idq is performed. These values are entered into a Matlab
® script and all possible values of modulators are obtained, delivering the overmodulation information. One way to see the region graphically is by considering a valid model (VM) and a not valid model (NVM), which represent no overmodulation and with overmodulation, respectively. Therefore, to avoid this issue, and considering a carrier from −0.5 to 0.5 p.u., the modulator amplitude condition that must be met is the following:
To graph the operating region considering the restriction in (12), the parameters of the system must be defined, which are listed in
Table 1.
The simulation is performed in Matlab
® to obtain the 3D plot of the operating region, presented in
Figure 2. The plot shows the overmodulation zone in red, and in blue when Equation (12) is satisfied, i.e., where the model is valid.
Other interesting cases to analyze involve how the operating region behaves under different power factors. The main studied cases are: (i) unity, (ii) inductive
pf of 0.93, and (iii) capacitive
pf of 0.93. The simulation results of the three cases are presented in
Figure 3, operating with a fixed frequency of 50 Hz. It can be observed that the converter has the greatest valid model (VM) region when the
pf is inductive, a smaller region for unity
pf, and the smallest VM region when capacitive
pf is imposed.
The effect of changing the
pf on the power converter operating region is summarized in
Table 2, where the unity
pf is taken as the reference. It is easy to note that the operating region is increased by around 12% for an inductive
pf and reduced by around 12% for a capacitive
pf. This result presents an opportunity to increase the operating region.
It is important to highlight that frequency affects the operating region because of the natural frequency response of inductances and capacitances.
Figure 4 shows the results considering a unitary
pf, tested for different frequencies: 30 Hz, 50 Hz, and 100 Hz. It can be observed that the greater the frequency value, the smaller the valid model operating region obtained. Therefore, with a lower frequency and more inductive
pf, a greater linear operation region can be obtained for the inverter. However, in this case, frequency is a disturbance and cannot be employed to enlarge the operating region, but the control can be aware that increasing frequency will decrease the operating region.
4. Phase-Locked-Loop
To ensure proper current control, the power converter needs to be synchronized with the grid mains, and the Phase Lock Loops (PLL) algorithm is one of the solutions to achieve phase tracking [
14].
The synchronizer used in this work seeks to maintain N samples per period independent of the grid frequency (external signal), i.e., the algorithm constantly changes the sampling time to track the grid frequency.
This PLL operates based on two variables: an external variable,
ug, and an internal variable. The external variable is sensed from the three-phase voltage and can be mathematically represented as:
Similarly, the internal variable has a behavior that can be described as:
where cosine(∙) represents the
N samples discretization of the cosine function, such that:
Each sampling time updates the values of
,
, and
according to:
with
,
, and
. Therefore, for the divisions to be integer numbers, this synchronizer requires a quantity of samples for periods in multiples of three for its correct operation.
The dot product between the external and internal signals gives the result in (19) as:
Moreover, to achieve synchronism, (19) requires that
2πn0(k)/N = θg(t) − 2mπ, or in other words,
u(
k) = 0, with
:
To track the grid frequency, a discrete PI controller is employed [
24]:
The algorithm, shown in
Figure 5, is tested by means of computer simulation using PSim
®, as shown in
Figure 6 and
Figure 7, such that there are
N = 204 samples per period. The PI controller gains
kp = 2.494 × 10
−4 and
a = 0.98.
In
Figure 6, it can be seen how the PLL correctly follows the frequency of the grid, increasing the value of N until it achieves 204. This is how the sinusoidal signal is constructed, and at 0.2 s and 0.35 s, it can be observed how the variation in frequency affects the signal, but the PLL instantly adjusts the new value of
Ts. In the same way, the synchronization of each phase of the grid can be observed in
Figure 7. It is important to note that the amplitude of the grid is divided by 311 for better comparison with the modulating reference signal, which uses the implemented algorithm as a speed reference. Throughout the simulation, it can be seen that the reference is in synchronization with the grid, including during the perturbation time.
5. MPPT-Based Perturbation and Observation
Electric power generation systems based on PV modules require DC/AC power converters to interface between the DC photovoltaic array and the AC electrical grid to inject the conditioned power. In some applications, this conversion is completed in two stages: first, a DC/DC conversion to extract the maximum power from the PV array and to increase the DC voltage on the DC-link capacitor, then a DC/AC inverter. However, in many photovoltaic generation systems, the implementation of the first DC/DC stage is not desired because it requires additional control and increases losses that can be avoided [
25,
26,
27]. Therefore, with a suitable design of the PV array, it is possible to reach a voltage level sufficient to feed the DC link. As result, the overall system increases in efficiency and only requires regulating the DC voltage based on environmental conditions. Therefore, in this research, a single-stage DC/AC converter is selected, which also tracks the MPP to maximize the power injected into the grid.
There are different MPPT algorithms, but the most commonly used method is P&O due to the simplicity of its implementation [
28]. Thus, the behavior of the P&O with variable step and its behavior when implemented in a power converter under distorted grid frequency is studied.
The principle of operation of the P&O algorithm is based on manipulating the DC-link voltage by increasing or decreasing DC-link voltage (
vdc) as a function of power variation to track the MPP. If the path specifies an increase in voltage, because increasing the voltage increases the power, the P&O algorithm will continue to do so until the power starts to decrease, at which point the algorithm takes the opposite action by reducing the voltage at the solar array. Additionally, a minimum and maximum voltage variation (
δvmin and
δvmax) are defined based on the point of operation on the P-V curve of the array. In this algorithm, it is necessary to calculate the power on each iteration and analyze it with respect to the previous one if the actual power,
p(
k), is greater or lesser than the previous
p(k − 1
). If the difference between
p(k) and
p(k − 1
) is greater than a tuning parameter
M (where
M depends on the array power and desired separation), the maximum voltage variation
δvmax is used. If the difference is less than M, it means the MPP has been reached and the minimum voltage variation
δvmin is used, as shown in
Figure 8.
As the MPPT algorithm runs in the same script as the converter’s main controller, both share the sampling time,
Ts, which, as previously mentioned, varies depending on the grid frequency. In order to minimize the voltage step variation per time caused by the variable sampling frequency, this paper proposes normalizing the voltage step variation as a function of the sampling time, as follows:
If this normalization is neglected, the MPPT algorithm exhibits oscillations, as shown in
Figure 9. However, the proposed normalization in (22) mitigates the effect of the variable sampling frequency, as also shown in
Figure 9.
6. Power Converter Control
After defining the MPPT, the VSI controller is studied in this section. The injection strategy is separated into three parts, as outlined below.
6.1. Non-Linear Feedback Control
From (4), the
dq-axis current dynamics can be expressed mathematically as:
It can be observed from Equations (23) and (24) that the system is nonlinear and coupled. Therefore, a nonlinear control strategy that employs state feedback is chosen. The exact linearization method is used by introducing two new variables,
w1 and
w2, as inputs to the power converter and by imposing that:
Thus, the Transfer Functions that relate
isd,
isq,
w1, and
w2, respectively, are integrators. By substituting
w1 and
w2 into Equations (23) and (24), the modulators are controlled to ensure a linear response within the boundaries of the operating region. Solving Equations (23) and (24) defines the modulators as functions of state variables, with new inputs
w1 and
w2 and the disturbances as follows:
It is important to note that these modulators enable the current response to be linear and decoupled within the operating region, thereby allowing for the easy definition of single-input single-output controllers.
6.2. PI Control of Currents
The control of the
isdq current is based on a PI controller that approximates the plant as a pure integrator, as shown in Equation (25). The PI controller has the following ideal transfer function:
Thus, the closed-loop transfer function results:
The dumping factor, ξ, is set to 0.707 through analytical tuning to achieve an overshoot of less than 5%. The natural frequency, ω
n, is set to 2π/
Tp, where
Tp =
k(L/R) and
k is a factor that represents the response of the closed-loop system. For
k = 0.5, the closed loop is faster, and the controller is designed with gains
KcPIi = 740 and
TiPIi = 2.7 × 10
3. However, the zero, caused by (
Tis + 1), can increase the overshoot, but this has been resolved by incorporating a filter reference (
Figure 10a).
6.3. Active Power PI Control Design
Additionally, to control the currents, the converter requires an outer power control loop, which directly impacts the ability to regulate the energy injected into the grid, the imposed power factor, and the required reactive power reference.
From a power balance of the system, the active power can be defined as ps = ppv − pRL − pCdc, where pRL is the power losses in the RL filter, pCdc is the power in the DC-link capacitor, ppv is the power supplied by the PV array, and ps is the active power injected into the grid.
For the active power control design, the plant transfer function in the Laplace domain is used:
HCdc(
s) = 2/(
Cdcs), which can be obtained from the expression of the energy in a capacitor [
29]. Given that the plant has only one pole, it is possible to use a PI control to achieve zero steady-state error. Therefore, the closed-loop transfer function for the power PI controller results in (30). Applying analytical tuning with ξ = 0.707,
Cdc = 2.35 mF, and
ωnPIi is at least 12 times faster than the current
ωnPIp. The gains for the power PI controller can be obtained as
KcPip =
Cdcξω
n and
TiPIp =
2 KcPip/(
Cdcξω
n2).
6.4. Reactive Power Control
Even though the main characteristic of a PV system is active power injection, related to the
d-axis current, it is possible to provide the power control system with reactive power regulation through the
q-axis current, as shown in
Figure 10, [
30]. The reactive power reference can be obtained as a function of the active power reference and the desired power factor, as:
From (31),
can be defined as a constant value,
Kpf; thus, (31) results in:
The block diagram presented in
Figure 10 shows the voltage/power control loop where the active and reactive power references are obtained. These are later introduced in the internal control loop as
dq-axis current references.
6.5. Overmodulation Control Based on Operation Region
One important drawback of a power converter connected to a microgrid is that the operating region (amplitude and frequency) may vary. In this paper, an overmodulation control is proposed to mitigate the effects of this reduction in the operating region. The results, shown in
Figure 3 and
Figure 4, conclude that as the frequency increases and the voltage amplitude increases, the operating region is reduced. Therefore, the proposed overmodulation control aims to counteract these effects. The proposed control diagram is simplified in
Figure 11 by introducing the main steps the algorithm followed to avoid overmodulation.
As presented in
Section 3, the valid area of the operating region is always greater when increasing the DC-link voltage (
vdc). This is independent of variations in the grid frequency, amplitude, or a non-unity power factor. As can be seen in
Figure 3 and
Figure 4, an increment of
vdc allows for a larger valid operating region. As the MPPT is based on varying the voltage in
vdc, when overmodulation occurs, the MPPT will be momentarily disconnected.
On the other hand, the effects of the power factor in the operating region, presented in
Figure 3, show that to relieve the converter in an overmodulation situation, changing from an operating state with a unity power factor to an inductive power factor allows for obtaining a greater region for the converter. Therefore, another part of the proposal to counteract overmodulation is to control the inverter between unity and inductive power factor to enlarge the VM region.
From the previous analysis, to enlarge the operating region in the proposed algorithm, both effects will be used, i.e., control over the DC-link voltage vdc and power factor pf.
The proposed algorithm, shown in
Figure 12 and
Figure 13, is separated into two main steps, explained in detail later. To avoid entering the overmodulation control loop in transient time, where the modulator may be out of range momentarily, the following sentence is first encountered:
where ε = [ε
1 ε
2]
T and ε
1, ε
2 are parameters to be syntonized, and for the experimental setup, a value of 5 × 10
−3 was used. If the condition in (33) is met, it can be said that the system is in a steady-state condition.
In order to avoid changing from step to step at all time, the algorithm will pass to the next step only if the “count” reaches “Time_wait”, raising the Action Time Control (ATC) flag “Flag_ATC” (see
Figure 12) and giving enough time to reach the steady-state condition, avoiding oscillations between the steps and cases (presented later), where “Time_wait” is a tuning parameter.
Given that the proposed control strategy is based on the exact value of the modulator
Mdq, a low-pass filter must be implemented to attenuate undesirable noise for this variable:
Using
N = 204 from PLL algorithm, it is obtained:
Overmodulation Control Step by Step
First Step: The first step looks for overmodulation, (12), while the power converter is in normal operation. If the filtered signal Mf is greater than Mmax, the first action is to enlarge the operating region by changing the power factor to 0.96 inductive, which gives room for the second step.
Second Step: The second step is activated only when the First Step has already been activated (in
Figure 13 related to “Step == 1”), and it is divided into three cases:
- ○
Case 1: If the filtered
Mf is greater than the maximum value
Mmax, then the worst state is happening. MPPT is deactivated to increase the voltage and therefore increase the operating region, such that:
In addition, the power factor is kept at 0.96 inductive.
- ○
Case 2: When the filtered Mf is within the range [Msub-max∙Mmax, Mmax], where Msub-max is an adjustable parameter of sub maximum. For this work, a value of Msub-max = 0.9 is used. Then, the system is kept in the actual position since it is possible that the problem (over voltage or over frequency) is still ongoing.
- ○
Case 3: If the filtered Mf is below Msub-max∙Mmax, then it is said that the system is out of trouble and can return to normal operation. Therefore, the MPPT is reactivated, the power factor is set to unity, and the “look up Step 1” is activated.
Thus, the aforementioned algorithm is able to deal with the overmodulation caused by grid disturbances.
6.6. Overmodulation Control Based on Operation Region
The primary advantages of the algorithm lie in its ease of implementation and the ability to be integrated with other injection-control techniques, with the goal of delivering greater stability to the grid. Conversely, a major disadvantage in its implementation is the high total computation time required for all associated strategies in systems characterized by imbalanced amplitude and a wide frequency range. However, the proposed strategy is characterized by low cost, which served as a focal point for improvement during the implementation process.
Table 3 illustrates the advantages and disadvantages of comparing the proposed algorithm with others found in the literature.
In [
31], The paper discusses different overmodulation strategies used in modern space vector PWM converters, analyzing their voltage transfer ratio (VTR) and low-order harmonics. The complexity of the modulation algorithm is elaborated, and the accuracy and linearization of the VTR are obtained through Fourier transform analysis. The current paper is appended to the table with the primary objective of categorizing all the strategies in comparison with the proposed approach presented in this study.
In [
32], an improved overmodulation strategy is proposed for a 3L-NPC-VSI in a maglev traction system. The strategy, based on the minimum amplitude error method, increases DC-bus voltage utilization and adopts a virtual space vector modulation strategy for NP voltage balance, CMV suppression, and leakage current suppression.
In [
33], a simple formulation of the virtual-vector-based PWM strategy for multilevel three-phase neutral-point-clamped (NPC) DC-AC converters is presented. The algorithm maintains capacitor voltage balance in every switching cycle and is applicable to any number of levels and operating conditions, including overmodulation. Simulation results confirm the effectiveness of the proposed formulation. This paper provides a comprehensive, effective, and simple PWM formulation for multilevel, three-phase NPC DC–AC converters.
The aforementioned studies are summarized in
Table 3.
8. Conclusions
Microgrids have been shown to exhibit variations in frequency and amplitude due to their low inertia. One of the main problems related to voltage variation is the change in the operating region, which means that the point of operation may fall out of the valid region. Therefore, in this work, an algorithm is proposed to deal with these changes by employing the characteristics of the power converter such as its operating region variation as a function of the power factor and the DC voltage. The control algorithm first includes reactive power to avoid overmodulation while still trying to keep the MPPT working to improve overall system efficiency. However, if this is not sufficient, the next step is to deactivate the MPPT algorithm to increase the DC voltage and therefore increase the operating region, trying to avoid overmodulation. Thus, maximum power extraction is sacrificed for a short period of time.
The algorithm, based on the aforementioned criteria, in conjunction with a nonlinear-based control for the dq-axis currents and PI strategy for the DC voltage, can work in an extended operating region based on the model certainties. The entire system response has satisfactory results in simulations and experiments, extending the operating region by employing the proposed algorithm by around 12% using the first step of reactive power control. Therefore, power converters can be employed in microgrids where amplitude and frequency can vary notably. For future work, it is of interest to implement the proposal in microgrids with an unbalanced system in voltage amplitude. In this way, this proposal contributes to keeping microgrid power generation systems in operation under the most unfavorable conditions, without the need to shut down power due to disturbances.