
Referring tothe study of epidemic mathematical models, this manuscript presents a noveldiscrete-time COVID-19 model that includes the number of vaccinated individuals as an additional state variable in the system equations. The paper shows that the proposed compartment model, described by difference equations, has two fixed points, i.e., a disease-free fixed point and an epidemic fixed point. By considering both the forward difference system and the backward difference system, some stability analyses of the disease-free fixed point are carried out.In particular, for the backward difference system a novel theorem is proved, which gives a condition for the disappearance of the pandemic when an inequality involving some epidemic parameters is satisfied. Finally, simulation results of the conceived discrete model are carried out, along with comparisons regarding the performances of both the forward difference system and the backward difference system.
Citation: A Othman Almatroud, Noureddine Djenina, Adel Ouannas, Giuseppe Grassi, M Mossa Al-sawalha. A novel discrete-time COVID-19 epidemic model including the compartment of vaccinated individuals[J]. Mathematical Biosciences and Engineering, 2022, 19(12): 12387-12404. doi: 10.3934/mbe.2022578
[1] | Sue Wang, Chaohong Zhou, Saleem Riaz, Xuanchen Guo, Haider Zaman, Alsharef Mohammad, Ahmad Aziz Al-Ahmadi, Yasser Mohammed Alharbi, NasimUllah . Adaptive fuzzy-based stability control and series impedance correction for the grid-tied inverter. Mathematical Biosciences and Engineering, 2023, 20(2): 1599-1616. doi: 10.3934/mbe.2023073 |
[2] | Pan Wang, Xuechen Li, Qianqian Zheng . Synchronization of inertial complex-valued memristor-based neural networks with time-varying delays. Mathematical Biosciences and Engineering, 2024, 21(2): 3319-3334. doi: 10.3934/mbe.2024147 |
[3] | Qi Wang, Yufeng Guo, Dongrui Zhang, Yingwei Wang, Ying Xu, Jilai Yu . Research on wind farm participating in AGC based on wind power variogram characteristics. Mathematical Biosciences and Engineering, 2022, 19(8): 8288-8303. doi: 10.3934/mbe.2022386 |
[4] | Yuhan Su, Shaoping Shen . Adaptive predefined-time prescribed performance control for spacecraft systems. Mathematical Biosciences and Engineering, 2023, 20(3): 5921-5948. doi: 10.3934/mbe.2023256 |
[5] | P. E. Greenwood, L. M. Ward . Rapidly forming, slowly evolving, spatial patterns from quasi-cycle Mexican Hat coupling. Mathematical Biosciences and Engineering, 2019, 16(6): 6769-6793. doi: 10.3934/mbe.2019338 |
[6] | Wenjing Wang, Jingjing Dong, Dong Xu, Zhilian Yan, Jianping Zhou . Synchronization control of time-delay neural networks via event-triggered non-fragile cost-guaranteed control. Mathematical Biosciences and Engineering, 2023, 20(1): 52-75. doi: 10.3934/mbe.2023004 |
[7] | Zhanying Tong, Yingying Zhou, Ke Xu . An intelligent scheduling control method for smart grid based on deep learning. Mathematical Biosciences and Engineering, 2023, 20(5): 7679-7695. doi: 10.3934/mbe.2023331 |
[8] | Yu Lei, Zhi Su, Xiaotong He, Chao Cheng . Immersive virtual reality application for intelligent manufacturing: Applications and art design. Mathematical Biosciences and Engineering, 2023, 20(3): 4353-4387. doi: 10.3934/mbe.2023202 |
[9] | Yixin Zhuo, Ling Li, Jian Tang, Wenchuan Meng, Zhanhong Huang, Kui Huang, Jiaqiu Hu, Yiming Qin, Houjian Zhan, Zhencheng Liang . Optimal real-time power dispatch of power grid with wind energy forecasting under extreme weather. Mathematical Biosciences and Engineering, 2023, 20(8): 14353-14376. doi: 10.3934/mbe.2023642 |
[10] | Zhen Yang, Zhengqiu Zhang, Xiaoli Wang . New finite-time synchronization conditions of delayed multinonidentical coupled complex dynamical networks. Mathematical Biosciences and Engineering, 2023, 20(2): 3047-3069. doi: 10.3934/mbe.2023144 |
Referring tothe study of epidemic mathematical models, this manuscript presents a noveldiscrete-time COVID-19 model that includes the number of vaccinated individuals as an additional state variable in the system equations. The paper shows that the proposed compartment model, described by difference equations, has two fixed points, i.e., a disease-free fixed point and an epidemic fixed point. By considering both the forward difference system and the backward difference system, some stability analyses of the disease-free fixed point are carried out.In particular, for the backward difference system a novel theorem is proved, which gives a condition for the disappearance of the pandemic when an inequality involving some epidemic parameters is satisfied. Finally, simulation results of the conceived discrete model are carried out, along with comparisons regarding the performances of both the forward difference system and the backward difference system.
To achieve the "double carbon goal" and develop green and low-carbon energy, distributed power generation technology has become an essential component of the new energy industry [1,2]. However, grid-tied inverters realize the integration of distributed power generation with the grid, exhibiting low inertia and weak damping, which is insufficient to provide the grid system's frequency and voltage support capability. The system power and frequency may lose self-regulation capability under the condition of disturbance, which seriously threatens the operation characteristics of the new power system [3,4]. Therefore, it is crucial to control the grid-tied inverter effectively. Several strategies were proposed to control VSG technology [5,6], focusing on frequency and voltage regulation features for the synchronous generator (SG). With extra-operational characteristics, the grid-tied inverter can maintain the grid power balance and frequency stability, improving the stability and safety of grid-tied system operation [7,8].
From the parameter optimization study for the GTVSG in [9,10], the rotational inertia (represented by J) and the damping factor (represented by D) determine the output response characteristics of the active power. The larger the rotational inertia J, the longer the regulation time, raising the active power overshoot and deteriorating system stability. On the other hand, when the damping factor D is larger, the active power overshoot decreases, the damping effect becomes better, and the dynamic response slows down. However, the independently increasing D seriously affects the frequency modulation characteristics of VSG active output power, leading to poor steady-state performance [11,12]. Therefore, the unreasonable selection of parameters makes the VSG prone to active power oscillations and poor active power transient and steady-state characteristics [13].
In [14] authors investigate the transient response features of GTVSG systems and suppress the KVA output power oscillations by adjusting the inertial damping. Similarly, [15] proposes an adaptive virtual inertia control approach for SG, taking advantage of the association between frequency fluctuation and virtual inertia and adjusting the inertia to enhance the microgrid frequency's dynamic characteristics. A virtual inertia adaptive control strategy for microgrids based on an improved bang-bang control strategy is proposed in [16], which uses the tunability of virtual inertia to reduce the VSG dynamic frequency deviation and improve frequency stability. To solve the problem of large fluctuations in frequency and power caused by excessive load fluctuations, reference [17] proposes a control strategy with adaptive adjustment of rotational inertia; the control strategy improves the system damping ratio and effectively suppresses the power overshoot and oscillation during transients.
The related literature work aimed to suppress the power fluctuation and frequency oscillation and improve the power output characteristics by adaptively changing the inertia and damping factor. However, in GTVSG, there is a contradiction between the active output power's transient and steady-state performance. In order to cope with these issues, according to the relationship between the system frequency change rate and inertia damping, authors in [10,18] proposed an optimal method for designing virtual inertia and damping factor for reduction of the output power and frequency fluctuations, improves dynamic power characteristics while minimizing active power steady-state errors. In [13,19], a differential compensation link is added to the forward channel of the VSG control structure. The pole-zero position of the active closed-loop transfer function is changed by adjusting the compensation coefficient, thus, improving the dynamic performance of the active power and ensuring the output power steady-state control accuracy. For the problem of mutual influence between the primary frequency modulation characteristics and damping characteristics in VSG control, the damping feedback link is introduced [20,21]. It achieves the decoupling of frequency regulation and damping in the steady-state case and solves the problem of the trade-off between the transient and steady-state characteristics of the VSG. In [22], a virtual damped feedforward VSG control mechanism is proposed that decouples the VSG frequency modulation characteristics from the damped power connection and enhances the steady-state features of the active output power when the grid frequency deviates. By adding a correction differential link to the VSG inertia damping characteristics at various stages, enhanced inertia damping characteristics are achieved [23].
The work presented in the literature has improved the steady-state performance of active output power, solving the contradiction between the transient and the steady-state performance of the active output power. However, when designing the parameters, there is still the problem that the damping factor affects the steady-state power accuracy and is too complicated to analyze the VSG active power output performances.
This paper proposes an active power control decoupling strategy for GTVSG based on DPDIDC. By connecting a PD control link in series with the VSG inertia forward channel and a PD control link in series with the VSG damping feedback channel, the VSG system structure is changed to optimize the performance of active output power. First of all, the power-frequency response of the GTVSG system is evaluated, and pitfalls in the transient and steady-state performance of the active output power are discovered using the GTVSG control concept. Secondly, pointing out the contradictory problems between the stability and transient features of the active output power. Then, a GTVSG, active decoupling control scheme, based on DPDIDC is proposed. The closed-loop transient and steady-state characteristics of the improved GTVSG system are analyzed. The parameter constraints of DPDIDC are given and finally, the correctness and superiority of the control scheme are verified by simulation.
The GTVSG central circuit topology is shown in Figure 1. eabc, uabc and ioabc are the VSG three-phase circuit potential, VSG output terminal voltage, and VSG output current; iabc represents the three-phase AC current flowing through inductor; ic-abc denotes capacitance current. The LCL filter is chosen to filter the inverter output harmonics, L1, 2 and R1, 2 are the filter inductor and parasitic resistance, C is the filter capacitor; Lg and Rg are the line inductor and resistance, Pref, Qref, ωo and UN are the reference values of the active output power, reactive power, the rated values of the angular frequency and voltage; Pe and Qe are the effective values of the VSG output active power and reactive power; θ and Em are the phase and amplitude of the active and reactive loop output reference voltages from the VSG control algorithm.
According to Figure 1, the VSG system control structure mainly includes four parts: voltage and current acquisition and processing power calculation, VSG control algorithm module (active-frequency control and reactive-voltage control), voltage and current double closed-loop control and modulation strategy.
The main control principle: firstly, the power is calculated by collecting the capacitor voltage and VSG output current, and obtain the output active power, reactive power, angular frequency and voltage. Then combined with the reference value of output active power and reactive power, the rated value of angular frequency and voltage, the phase and amplitude of the output reference voltage are calculated through the active power control and reactive power control. The output voltage is controlled by voltage and current double loops, the AC voltage signal is tracked, and finally the inverter is controlled by modulating the drive signal.
Drawing on the mathematical model of synchronous generators, while simulating the inertia and damping characteristics of synchronous generators, the active power controller of VSG consists of two parts: the rotor equation of motion and the active-frequency droop control. The VSG rotor equation of motion which is obtained as follows [24]:
{JωodΔωdt=Pm−Pe−D(ω−ωo)dθdt=ω | (1) |
where J is the rotational inertia; D is the damping factor; Pm is the VSG virtual mechanical power; Pe is the VSG output active power; ω and ωo are the effective and rated values of the VSG angular frequency; and θ is the phase angle of the VSG. VSG uses P-ω droop control to simulate the governor of the synchronous generator, to achieve the frequency modulation characteristic of the GTVSG system [9], the expression is:
Pm=Pref−Kf(ω−ωo) | (2) |
where Kf is the droop factor of active power; Pref is the reference value of VSG active power output.
The active-frequency control expression is obtained as:
{JωodΔωdt=Pref−Pe−(D+Kf)(ω−ωo)dθdt=ω | (3) |
Therefore, the VSG active-frequency control block diagram can be obtained as shown in Figure 2.
The droop part of the controller implements the primary frequency regulation function, and the rotor motion characteristics realize the inertia damping support. The controller outputs the frequency of the reference voltage signal, which is integrated to form the phase angle.
VSG achieves output reactive power control by regulating the interrelationship between reactive power and voltage through the reactive power droop factor. So the VSG uses Q-U droop control to simulate the regulating characteristics of the synchronous generator [9], then the voltage-reactive power droop relationship equation is:
Qref−Qe=Kq(UN−U) | (4) |
where Em is the VSG output voltage; Kq denotes droop factor of KVAR; Qe and Qref represent the actual and reference values for the VSG output; UN is the VSG rated voltage value. The dynamic response of reactive-voltage regulation can be described using first-order inertia links.
The introduction of the 1/Kis link makes the dynamic response of the excitation regulation inertial in nature and does not undergo step changes. This feature is closer to the synchronous motor in transient characteristics to prevent the power electronics from being too sensitive to reactive load changes in reactive load, causing system instability. Its reactive-voltage control expression is:
Em=1KiS[Qref−Qe+Kq(UN−U)] | (5) |
where the Ki is the factor of relevant power integration. Therefore, the VSG reactive-voltage control block diagram can be obtained as shown in Figure 3.
VSG reactive-voltage control enables the inverter to have the function of voltage regulation by adjusting the voltage amplitude. After reactive-voltage control to obtain the amplitude of the reference voltage, combined with the phase obtained by active-frequency control can synthesize the reference voltage eabc=[emaembemc]T, the expression is:
{ema=Emsinθema=Emsin(θ−2/23π3π)ema=Emsin(θ+2/23π3π) | (6) |
To maintain the output voltage stability of the virtual synchronous generator, tracking control of the voltage is required. The PI controller is used here to achieve tracking control of the voltage signal. In VSG voltage-current dual closed-loop control, the outer-loop voltage feedback is generally selected as the capacitor voltage and the inner-loop feedback is selected as the capacitor current. Figure 4 shows the VSG voltage and current double closed loop control block diagram.
where Ucref (s) is capacitance voltage command value; GuPI (s) is voltage outer loop regulator; GiPI (s) is current internal loop regulator; U(s) is the VSG output capacitor voltage; i(s) is the inductor current; ic(s) is the capacitor current; 1/sL1, 1/sL2 is the filter inductor equivalent impedance; 1/sC is the filter capacitor equivalent impedance; io(s)is the inverter output current; KPWM is the inverter equivalent gain, using PWM modulation, KPWM = l.
As the LCL filter generates resonance problem, it can add damping to the system, so in the VSG voltage and current dual-loop control, the output capacitor voltage is selected as the voltage outer-loop feedback quantity, the capacitor current is selected as the current inner-loop feedback quantity, which not only inhibits resonance, but also solves the output coupling problem. At the same time, the capacitive current inner-loop control as the feedback has more negligible output impedance, insensitive to the disturbance changes of load, and has a stronger carrying capacity [25]. Its voltage and current dual-loop both use the PI controllers to achieve the purpose of tracking and control the voltage signal.
According to the VSG active-frequency control expression, the power-frequency characteristic of its system is mainly reflected in the P-ω characteristic. Then the expression of the P-ω characteristic is:
GP−ω(s)=ω−ωoPref−Pe=ΔωΔPe=1Jωos+D+Kf | (7) |
According to Eq (7), the parameters affecting the system performance are J, D, and Kf, the power-frequency transient characteristics of the system are related to J, D, and Kf. In contrast, the power-frequency steady-state characteristics are related to D and Kf. The system can control Δω by adjusting J and D to ensure the stable output of the system power-frequency output response. Since D and Kf have the same effect on the system characteristics, generally, the fixed value of Kf is 100% of the output power variation, and the angular frequency variation is Δω = 1%ωo [26,27]. Therefore, its effect on the VSG active loop power-frequency response is analyzed by varying the magnitude of D.
As we can clearly see the Figure 5 shows the VSG active-power control structure and terminologies of the GTVSG closed-loop control block diagram.
KP is the gain factor of the relationship between the active output power Pe and the phase angle θ in Figure 5. According to the VSG line power transmission principle, the GTVSG system output power is:
{Pe=3EmUgXθQe=(3Em−Ug)X⋅Ug | (8) |
where Ug is the grid voltage amplitude, and X is the GTVSG output equivalent impedance, according to Eq (8): KP = 3EmUg /X [24]. The terminologies used in the above Figure 5 for the GTVSG active closed-loop control block diagram, the following expression can be drawn for a grid-connected inverter output active power Pe such that:
Pe=KPJωos2+(D+Kf)s+KPPref+KP(Jωos+D+Kf)Jωos2+(D+Kf)s+KP(ωo−ω)=GPe−PrefPref+GPe−ω(ωo−ω) | (9) |
From Eq (9), the active output power Pe is affected by the VSG active power reference value and the fluctuation of the grid frequency. Thus, to recapitulate the final value theorem, the stable active power Peo for a grid-connected operation can be described as follows:
Peo=lims→0[GPe−PrefPref+GPe−ω(ωo−ω)]=Pref+(D+Kf)(ωo−ω) | (10) |
According to Eq (10), the steady-state value of active output power Peo is affected by two conditions: most important one is the change of active power reference value Pref, while the other one is the stable power deviance (D + Kf) (ωo - ω) due to frequency disturbance.
When the frequency fluctuations are higher, it is noticeable that the relevant steady-state error is larger. It is noticeable that parameters D and Kf mainly affect the steady-state enactment of active power Peo. While keeping Kf is fixed, the greater the damping factor D is, the higher frequency disturbance capability is. On the other hand the weaker the frequency modulation characteristics of the VSG active output power means the larger the error. Therefore, the VSG damping characteristics affect the frequency modulation capability of the VSG system, which affects the steady-state performance of the active output power. The step response of different D on the anti-frequency disturbance capability of VSG steady-state active output power is shown in Figure 6.
In Figure 6, when the active power reference value is 0, the frequency fluctuations (frequency deviation of 0.01 Hz), the damping power D seriously affects the frequency modulation accuracy of the VSG output steady-state active power. While the primary frequency modulation capability of the VSG system is very slow, and the stable output active power varies significantly with the increase of D. Hence, the steady-state performance of the active output power goes worse according.
The overall closed-loop small-signal transfer function may be used to study active-loop power-frequency dynamic response features of GTVSG systems [19]. In Figure 5, by expressing the relationship between the active output power and its reference value, the active output power of closed-loop small-signal transfer function is obtained as shown in Eq (11). By relating the angular frequency and the active power reference value, we can easily get the output angular frequency closed-loop small-signal transfer function as shown in Eq (12).
GPe=ΔPeΔPref=KPJωos2+(D+Kf)s+KP | (11) |
Gω=ΔωΔPref=sJωos2+(D+Kf)s+KP | (12) |
According to Eqs (11) and (12), the small signal transfer function for closed-loop active power is a typical second-order system. The VSG KVA power can be deduced with the help of pole placement such as:
s1,2=−(D+Kf)±√(D+Kf)2−4JωoKP2Jωo | (13) |
From Eq (13), the J and D affect the poles s1 and s2 of the system active closed loop. The unreasonable selection of parameters affects the distribution of zero poles, which leads to the degradation of the transient performance of the GTVSG system in the dynamic adjustment process of output active power and frequency.
According to the expression of the transfer function in Eq (11), the step response can be obtained by taking the inverse Laplace transform:
ΔPe(t)=(1−1√1−(D+Kf)24JωoKPe−D+Kf2Jωotgsin(√KPJωo√1−(D+Kf)24JωoKPt+arccosD+Kf2√JωoKP))ΔPref(t) | (14) |
According to Eq (14), the peak time tpp, regulation time tsp, and the power overshoot σp% of the active step response can be expressed as:
{tpp=2πJωo√8JωoKP−(D+Kf)2tspe=8JωoD+Kfσp%=e−πD+Kf√4JωoKP−(D+Kf)2 | (15) |
where above expressions demonstrate that the tpp, tsp, and σp% in Eq (15), a graph of the relationship between the transient performance for an active power and the parameters J, D are drawn according to simulation. Furthermore we have analyzed the effects of J and D on the transient performance index. The variation trend of active output power transient performance index with different J and D can be seen in the following two figures (see Figures 7 and 8(a)–(c)).
Figures 7 and show that when J is larger, the damping ratio of the system is smaller, the peak time tpp, regulation time tsp and overshoot σp% increase, the power response slows down, and the dynamic performance and stability become worse; When D is larger, the damping ratio of the system gradually increases, the peak time tpp of active power increases, the regulation time tsp and overshoot σp% decrease, and the damping effect becomes better, the power oscillation slows down, the dynamic response slows down, and the stability becomes better.
Next, the transient performance of the system affected by J and D is simulated to verify the correctness of the above theoretical analysis. Firstly, when the active power reference value is 10kW, observe the effect of J and D on the active output power step response. The step response of active power with different J and D is shown in Figure 9.
Figure 9 shows that when J increases, the peak time observed from step response grows from 1.5 to 1.8 s, the regulation time increases from 4.8 to 5.7 s, and the power overshoot increases from 44 to 53%. When D increases, the peak time of active step response increases from 1.6 to 1.8 s, the regulation time decreases from 6.3 to 4.3 s, and the power overshoot decreases from 54 to 25%. The simulation results produced are in line with the theoretical results, verifying the variation trend of active output power transient performance index with different J and D described in this paper. The comparison of the system power transient performance indicators is shown in Table 1.
Transient performance indicators | tpp (s) | tsp (s) | σp% | |
D constant | J = 10 | 1.5 | 4.8 | 44 |
J = 15 | 1.8 | 5.7 | 53 | |
J constant | D = 10 | 1.6 | 6.3 | 54 |
D = 20 | 1.8 | 4.3 | 25 |
Similarly, take the inverse Laplace transform of Eq (12), and obtain the unit step response of frequency as:
Δω(t)=((√KPJωo/√1−(D+Kf)24JωoKPe−D+Kf2Jωot)gsin(√KPJωo√1−(D+Kf)24JωoKPt))ΔPref(t) | (16) |
According to Eq (16), the peak time tpω, regulation time tsω and power overshoot σω% of the system frequency step response can be expressed as:
{tpω=2JωoKParctan√4JωoKP−(D+Kf)2D+Kf√4JωoKP−(D+Kf)2tsω=8JωoD+Kfσω%=√4KPJωoe−πD+Kf√4JωoKP−(D+Kf)2arctan√4JωoKP−(D+Kf)2D+Kf | (17) |
Similarly, according to the expressions of tpω, tsω and σω% in Eq (17), a graph of the relationship between the transient performance index of the system angular frequency and J, D is drawn by simulation. Also we have analyzed the effects of J and D on the transient performance index of the angular frequency. The variation trend of the frequency transient performance index with different J and D is shown in Figures 10 and 11(a)–(c).
Figures 10 and 11 show that J can suppress the frequency variation and determines the oscillation frequency in the active output power response. When J is larger, the peak time tpω and regulation time tsω lengthen, the overshoot σp% decreases, and the system frequency oscillation time is too long, which is not conducive to system stability; D can reduce the amount of frequency deviation and determine the active power response oscillation decay rate. When D is larger, the peak time tpω, regulation time tsω and overshoot σω% decrease, and the dynamic performance becomes better.
When the system frequency rating is 50 Hz, observe the effect of J and D on the output frequency step response. The step response of frequency with different J and D is shown in Figure 12.
Figure 12 shows that when J increases, the peak time of the frequency step response increases from 0.7 to 0.9 s, the regulation time increases from 5.2 to 6.3 s, and the frequency overshoot decreases from 0.30 to 0.25%. When D increases, the peak time of the frequency step response decreases from 0.9 to 0.7 s, the regulation time decreases from 6.3 to 4.4 s, and the frequency overshoot decreases from 0.3 to 0.23%. The simulation results are consistent with the theoretical analysis results, which verifies the variation trend of the frequency transient performance index with different J and D described in this paper. The comparison of the system frequency transient performance indicators is shown in Table 2.
Transient performance indicators | tpp (s) | tsp (s) | σp% | |
D constant | J = 10 | 0.7 | 5.2 | 0.30 |
J = 15 | 0.9 | 6.3 | 0.25 | |
J constant | D = 10 | 0.9 | 6.3 | 0.30 |
D = 20 | 0.7 | 4.4 | 0.23 |
According to the analysis of the system active loop power-frequency response in Sections 3.1 and 3.2. It is obvious that the choice of J and D for the output active power has a great influence on the characteristics response, when VSG system is connected to grid. In GTVSG control strategy has a contradiction in the optimal choice of parameters for the regulation of the transient and steady-state features of the power (KVA). While the transient response of active power and frequency is associated with both J and D. We can easily chose smaller J and larger D which is advantageous to the transient characteristics of the system. However, the steady-state response of active power is only linked to D. When D increases, it means the rise in frequency deviation, the GTVSG steady-state power error increases, which is not conducive to frequency steady-state regulation, and the active output power steady-state performance deteriorates. In summary, the selection of J and D cannot satisfy the requirements of the transient as well as steady-state characteristics of the active output power, so the structure of the VSG needs to be reasonably designed to improve and satisfy the performance requirements of the system. Therefore, to address these issues, the paper mainly introduces PD control to the inertia damping channel, which varies the GTVSG closed-loop construction. That is decoupling control on the transient and stable performance of GTVSG system.
According to the classical control theory, PD control and feedback control are commonly used to improve the performance of second-order systems. PD control produces a correction effect in advance of the position error of the output signal, so as to achieve the purpose of improving the system performance. It is equivalent to increasing the system damping, adding a closed-loop zero to the system, so that the system step response overshoot is reduced and the setting time is shortened, and does not affect the constant steady-state error of the system. Due to the adoption of differential terms, larger open-loop gain is allowed to be selected, so the steady-state error can be reduced under ensure dynamic performance conditions [28]. Therefore, PD control is introduced in the VSG inertia forward channel and the damping feedback channel respectively. The VSG control block diagram based on dual PD control inertia damping is shown in Figure 13.
According to Figure 13, the expression of GTVSG active output power Pemod based on duplex PD control inertia damping is:
Pemod=KP+KKPs(Jωo+D)s2+(Kf+KKP)s+KPPref+KP(Jωo+D)s+KfKP(Jωo+D)s2+(Kf+KKP)s+KP(ωo−ω) | (18) |
In a view of Eq (18), the GTVSG system with dual PD control enhances a closed-loop response. The zero pole placement could be adjusted and controlled by the differential gain K, thus satisfying the requirements of the system transient characteristics index.
According to the final value theorem, the active output power Peomod expression of the GTVSG system in the steady-state case is:
Peomod=Pref+Kf(ωo−ω) | (19) |
According to Eq (19), the active output power steady-state value Peomod consists of active power reference value Pref and primary frequency modulation link Kf (ωo-ω). The VSG damping characteristics do not affect the primary frequency modulation characteristics of the active power control. There is no steady-state error caused by D during GTVSG steady-state operation, and improves the stable performance of the KVA power. Without introducing other parameters, PD control is added to the feedback channel to improve the steady-state performance of the active power. Based on the inertia damping control structure with dual PD control, the step response of D on the improved VSG steady-state output active power against frequency disturbance is obtained as shown in Figure 14.
Figure 14 shows that, after adding the PD control, observing the influence of D on the improved GTVSG steady-state output active power. Similarly, when the active power reference value is set to 0, after the VSG inertia damping control structure changes, when the frequency fluctuates, the output of the steady output power changes little. The D has no effect of a frequency regulation for the VSG output steady-state output power. We know that the output power is dynamically adjusted with the help of frequency regulation characteristics. Also consecutively reaching to a fixed value and then steady-state output, which can realize the steady-state operation of the active power. Simultaneously, with the increase of D, the active power oscillation slows down and improves the system transient features.
We can easily adjust the pole placement with the choice of K value of grid-connected active power that can be changed to improve the transient characteristics of grid-connected active power. Referring to the system simulation parameters in Table 3, when J = 2, D = 15 and K = 0–0.2, the GTVSG active closed-loop zero-pole placement is obtained which is illustrated in Figure 15.
Parameters | Value | Parameters | Value |
DC voltage Udc | 800 V | Line inductors Lg | 0.2 mH |
Rated voltage UN | 311 V | Line resistance Rg | 0.1 Ω |
Grid voltage Ug | 311 V | Rotational inertia J | 6 kg·m2 |
Rated frequency fo | 50 Hz | Damping factor D | 20 |
Filter inductors L1 | 5 mH | Active droop Kf | 20 |
Filter inductors L2 | 5 mH | Reactive droop Kq | 30 |
Filter capacitor C | 20 µF | Reactive power integration Ki | 50 |
Filter resistors R1 | 0.01 Ω | Active power reference Pref | 10 kW |
Filter resistors R2 | 0.01 Ω | Reactive power reference Qref | 0 Var |
From Figure 15, as the K increases further, the zero z of the closed-loop system transfer function rapidly approaches the imaginary axis. The two poles s1 and s2 approach the imaginary axis, gradually approach the real axis and meet at the real axis. The system enters the over-damping state, the damping ratio increases, the power oscillation amplitude decreases, and the overshoot decreases. Therefore, the added differential term has the same effect as D. The transient response of the system can be improved by increasing the K without affecting the steady-state evaluation of the system.
By improving the GTVSG active link loop, it is found that the selection of J and D can take into account the requirements of both the transient and steady-state characteristics of the active output power. Simultaneously, dual PD control increases the selection range of parameters of the inertial damping. The choice of the K has an essential impact on the active power performance. By adjusting the zero-pole position of the grid-connected system, the overall characteristics of the system are improved, so as to realize the decoupling control of the steady-state and transient characteristics of the active power.
According to Eq (18), based on the dual PD control inertia damping, the natural oscillation frequency ωnmod and damping ratio ζmod of the GTVSG transfer function can be obtained:
{ωnmod=√KPJωo+Dζmod=KKP+Kf2√KP(Jωo+D) | (20) |
In Eq (20), dynamic adjustment of K and J can effectively realize the damping regulation of VSG active closed-loop system. Increasing the K can raise the damping ratio of the system and satisfy the constraints of stability margin, overshoot and damping ratio.
1) Stability margin
The open-loop transfer function Gmod(s) of the active power control loop based on the dual PD control is:
Gmod(s)=KKPs+KP(Jωo+D)s2+Kfs | (21) |
Convert Eq (21) to a standard Type Ⅰ system:
Gmod(s)=KP/Kf(Ks+1)s(Jωo+DKfs+1) | (22) |
when the open-loop gain of the system is H = KP / Kf and the time constant is T = (Jωo + D) / Kf, then Eq (22) can be transformed into:
Gmod(s)=H(Ks+1)s(Ts+1) | (23) |
According to Eq (23), the open-loop frequency characteristics of the system can be obtained:
{Gmod(jω)=H(jωK+1)jω(jωT+1)A(ω)=H√(ωK)2+1ω√(ωT)2+1φ(ω)=−90∘+arctgωK−arctgωT | (24) |
where A(ω) is the amplitude and ϕ(ω) is the phase angle. According to classical control theory, the stability of a system in the frequency domain can be measured in terms of stability margins (phase angle margin γ and amplitude margin h) [28]. Thus the phase angle margin γ and amplitude margin h can be expressed as:
{h=1A(ωc)=+∞γ=180∘+φ(ωc)=90∘+arctgωcK−arctgωcT | (25) |
where ωc is the cut-off frequency, the phase angle margin γ becomes larger with the increase of the K, which is beneficial to the system's stability. The cut-off frequency ωc of the system:
A(ωc)=H√(ωcK)2+1ωc√(ωcT)2+1=1 | (26) |
Equation (26) can be solved as follows:
ωc=√(HK)2+1+√(2HT)2+(1−(HK)2)2√2T | (27) |
In Eq (27), dωc / dK < 0, so the system cut-off frequency decreases with the increase of the K. The reduction of the system cut-off frequency improves the system's ability to suppress frequency oscillations and improves the stable response of the KVA power. To obtain satisfactory performance, the phase angle margin γ and amplitude margin h should satisfy the following conditions:
{30∘⩽γ⩽60∘h⩾10dB | (28) |
The value range of K can be obtained as:
{h=1A(ωc)=+∞⩾10dBJωo+D−√3Kfωc(√3(Jωo+D)+Kf)⩽K⩽3(Jωo+D)−√3Kfωc(√3(Jωo+D)+3Kf) | (29) |
2) Damping ratio
When designing the PD control system, appropriately select a higher open-loop gain to not only improve the system transient characteristics, but also compensate for the loss of steady-state error. To satisfy the transient performance index with moderate overshoot and short setting time, the damping ratio ζmod should fulfill the following conditions:
0.6⩽ζmod=KKP+Kf2√KP(Jωo+D)⩽0.8 | (30) |
Therefore, according to the system constraints, the value range of K and the corresponding parameters, the range of K and J are obtained as:
{0.01⩽K⩽0.20.2⩽J⩽15 | (31) |
We can testify the rationality and superiority of the GTVSG system with dual PD control inertia damping proposed in the paper by using Matlab/Simulink. We have to build the GTVSG central circuit topology and control simulation model, where the system central circuit and control parameters are shown in Table 3.
To verify the correctness of the control scheme, the effects of D and K on the transient and steady-state performance of the GTVSG active power based on dual PD control are analyzed by simulation. In analyzing the active output power response of the GTVSG system under disturbance, this paper follows two cases: (ⅰ) the step response change of VSG output power and output frequency when the active power reference value changes abruptly; (ⅱ) the step change of VSG active output power when the grid frequency changes abruptly.
Under the VSG grid-connection, the active power reference value Pref is 10 kW. While the reference value Pref rises from 10 kW to 20 kW at 0.5 s. It can easily observed that the active output power response of the GTVSG system in four cases (J = 5, D = 0; J = 5, D = 20; J = 5, D = 0, K = 0.05; J = 5, D = 20, K = 0.05) is different. The simulation waveform diagram in Figure 16, the corresponding simulation waveform diagram of output frequency in Figure 17.
In Figures 16 and 17, when the active power reference value increases, the VSG active power and frequency will step up. After dynamic adjustment, the active power output reaches the reference value and the frequency returns to the rated operating state, which is aligned with the grid frequency. Through the comparison of the four cases, it is found that with the increase of D, the transient performance of the system output power and frequency is better, the peak of active power fluctuation becomes smaller, the time to reach the steady-state is shorter, and the power oscillation is reduced, which helps to improve the transient performance of the system output response. With the dual PD control strategy and the introduction of K, the active power changes smoothly, which vastly reduces the overshoot and regulation time of the system, accelerates the dynamic response, and enters the stable state quickly. At the same time, K helps to improve the frequency modulation characteristics of the system, and the transient performance of the output power and frequency also is greatly improved.
In Figures 18 and 19, under the VSG grid connection, the simulation waveform of VSG active output power when D and K are adjusted for a 0.05 Hz increase in grid frequency at 0.4 s.
In Figure 18, when D increases, reducing the power oscillation, because it affects the frequency modulation characteristics of the power, the active output power is lower than the regulated steady-state value, which generates a steady-state error and affects the steady-state performance of the VSG active power. After adopting the inertia damping control strategy with dual PD control, the output power is consistent with the steady-state power when reaching the steady -state, and the active power steady-state error is 0. Meanwhile, the power has no oscillation, the regulation time and overshoot are reduced, and the system soon enters the steady-state. Therefore, the inertia damping control strategy based on dual PD control not only eliminates the steady-state error, but also improves the steady-state performance of VSG active power as well as enhances the transient performance of GTVSG active power.
The simulation waveform of VSG active output power under grid frequency mutation with different K in Figure 19. By comparison, it is found that when the K is increased in a particular range, it not only does not disturb the stable response of GTVSG active power, but also inhibits the power oscillation. The active power overshoot and regulation time are reduced, it quickly enters stability, speeding up the response, and the transitory performance of VSG active power is greatly improved.
The power response performance of the system is compared between the differential compensation virtual inertia control strategy and the dual PD inertia damping control strategy to verify the superiority of the proposed control method. The comparison charts of the simulation results are shown in Figures 20–22.
1) Comparison of system transient performance
The comparison shows that after adding D, the transient performance of the system output power and frequency is better, the peak of active power fluctuation becomes smaller, the time to reach the steady state is shorter, and the power oscillation is reduced, which helps to improve the transient performance of the system output response. Therefore, the dual PD control strategy of adding D and K is more helpful to improve the frequency modulation characteristics of the system, and the transient performance of the system output power and frequency is greatly improved.
2) Comparison of system steady state performance
The comparison shows that the system output active power of differential compensation virtual inertia control strategy still has steady-state error, which affects the steady-state performance of VSG active power. After adopting the inertia damping control strategy with dual PD control, the output power is consistent with the steady-state power when reaching the steady-state, and the active power steady-state error is 0. The steady-state performance of VSG active power is greatly improved.
This paper presented a decoupling control for GTVSG founded on DPDIDC, which solves stability and transient response control problems. The theoretical analysis and simulation results show that: 1) The smaller J and the larger D are beneficial to the transient characteristics of the system. Still, too large D will increase the frequency deviation, GTVSG steady-state power error increases, and the active output power steady-state performance deteriorates. The selection of J and D is contradictory to the optimal control of active power transient and steady-state characteristics. Thus it cannot satisfy the requirements of active output power transient characteristics and steady-state characteristics simultaneously. 2) Using the proposed control strategy, the VSG inertia forward channel and the damping feedback channel are connected with PD control. With fewer parameter variables introduced, changing the structure of the GTVSG system eliminates the steady-state error in power and improves the transient characteristics of the system in both frequency and power. 3) The parameters are optimized according to the performance index requirements of the dual PD control active output power, obtaining the constraints on the parameters and the range of values.
This research work was supported by Taif University Researchers Supporting Project number (TURSP-2020/144), Taif University, Taif, Saudi Arabia.
The authors declare there is no conflict of interest.
[1] |
H. W. Hethcote, The mathematics of infectious diseases, SIAM Rev., 42 (2000), 599–653. https://doi.org/10.1137/S0036144500371907 doi: 10.1137/S0036144500371907
![]() |
[2] |
T. T. Marinov, R. S. Marinova, Dynamics of COVID-19 using inverse problem for coefficient identification in SIR epidemic models, Chaos Solit. Fract., 5 (2020). https://doi.org/10.1016/j.csfx.2020.100041 doi: 10.1016/j.csfx.2020.100041
![]() |
[3] |
J. T. Wu, K. Leung, M. Bushman, N. Kishore, R. Niehus, P. M. de Salazar, et al., Estimating clinical severity of COVID-19 from the transmission dynamics in Wuhan, China, Nat. Med., 26 (2020), 506–510. https://doi.org/10.1038/s41591-020-0822-7 doi: 10.1038/s41591-020-0822-7
![]() |
[4] |
S. Mangiarotti, M. Peyre, Y. Zhang, M. Huc, F. Roger, Y. Kerr, Chaos theory applied to the outbreak of COVID-19: An ancillary approach to decision making in pandemic context, Epidem. Infect, 148 (2020), 1–29. https://doi.org/10.1017/S0950268820000990 doi: 10.1017/S0950268820000990
![]() |
[5] |
S. Gounane, Y. Barkouch, A. Atlas, M. Bendahmane, F. Karami, D. Meskine, An adaptive social distancing SIR model for COVID-19 disease spreading and forecasting, Epidem. Methods, 10 (2021), 20200044. https://doi.org/10.1515/em-2020-0044 doi: 10.1515/em-2020-0044
![]() |
[6] |
G. Giordano, F. Blanchini, R. Bruno, P. Colaneri, A. Di Filippo, A. Di Matteo, et al., Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy, Nat. Med., 26 (2020), 855–860. https://doi.org/10.1038/s41591-020-0883-7 doi: 10.1038/s41591-020-0883-7
![]() |
[7] |
A. Ajbar, R. T. Alqahtani, M. Boumaza1, Dynamics of an SIR-Based COVID-19 Model With Linear Incidence Rate, Nonlinear Removal Rate, and Public Awareness, Front. Phys., (2021). https://doi.org/10.3389/fphy.2021.634251 doi: 10.3389/fphy.2021.634251
![]() |
[8] |
P. Kumara, V. S. Erturk, M. Murillo-Arcila, A new fractional mathematical modelling of COVID-19 with the availability of vaccine, Results Phys., 24 (2021), 2211–3797. https://doi.org/10.1016/j.rinp.2021.104213 doi: 10.1016/j.rinp.2021.104213
![]() |
[9] |
N. Gozalpour, E. Badfar, A. Nikoofard, Transmission dynamics of novel coronavirus SARS-CoV-2 among healthcare workers, a case study in Iran, Nonlinear Dynam., 105 (2021), 3749–-3761. https://doi.org/10.1007/s11071-021-06778-5 doi: 10.1007/s11071-021-06778-5
![]() |
[10] |
E. Badfar, E. J. Zaferani, A. Nikoofard, Design a robust sliding mode controller based on the state and parameter estimation for the nonlinear epidemiological model of Covid-19, Nonlinear Dynam., (2021), 5–-18. https://doi.org/10.1007/s11071-021-07036-4 doi: 10.1007/s11071-021-07036-4
![]() |
[11] |
A. Rajaei, M. Raeiszadeh, V. Azimi, M. Sharifi, State estimation-based control of COVID-19 epidemic before and after vaccine development, J. Pro. Control, 102 (2021), 1–14. https://doi.org/10.1016/j.jprocont.2021.03.008 doi: 10.1016/j.jprocont.2021.03.008
![]() |
[12] |
M. De la Sen, A. Ibeas, R. Nistal, About partial reachability issues in an SEIR epidemic model and related infectious disease tracking in finite time under vaccination and treatment controls, Discrete Dynam. Nat. Soc., (2021). https://doi.org/10.1155/2021/5556897 doi: 10.1155/2021/5556897
![]() |
[13] |
M. De la Sen, A. Ibeas, On an SE(Is)(Ih)AR epidemic model with combined vaccination and antiviral controls for COVID-19 pandemic, Adv. Difference Equ., 92 (2021). https://doi.org/10.1186/s13662-021-03248-5 doi: 10.1186/s13662-021-03248-5
![]() |
[14] |
S. Zhai, G. Luo, T. Huang, X. Wang, J. Tao, P. Zhou, Vaccination control of an epidemic model with time delay and its application to COVID-19, Nonlinear Dynam., 106 (2021), 1279–1292. https://doi.org/10.1007/s11071-021-06533-w doi: 10.1007/s11071-021-06533-w
![]() |
[15] |
E. Hwang, Prediction intervals of the COVID-19 cases by HAR models with growth rates and vaccination rates in top eight affected countries: Bootstrap improvement, Chaos Solit. Fract., 155 (2022). https://doi.org/10.1016/j.chaos.2021.111789 doi: 10.1016/j.chaos.2021.111789
![]() |
[16] |
P. Mahmood, M. Saeed, Stability of the equilibria in a discrete-time sivs epidemic model with standard incidence, Filomat, 33 (2019), 2393–2408. https://doi.org/10.1016/j.chaos.2021.111789 doi: 10.1016/j.chaos.2021.111789
![]() |
[17] |
M. De la Sen, S. Alonso-Quesada, A. Ibeas, On a Discrete SEIR Epidemic Model with Exposed Infectivity, Feedback Vaccination and Partial Delayed Re-Susceptibility, Mathematics, 9 (2021), 5–9. https://doi.org/10.3390/math9050520 doi: 10.3390/math9050520
![]() |
[18] |
M. De la Sen, S. Alonso-Quesada, A. Ibeas, R. Nistal, On a Discrete SEIR Epidemic Model with Two-Doses Delayed Feedback Vaccination Control on the Susceptible, Vaccines, 9 (2021). https://doi.org/10.3390/vaccines9040398 doi: 10.3390/vaccines9040398
![]() |
[19] |
Y. Omae, Y. Kakimoto, M. Sasaki, J. Toyotani, K. Hara, Y. Gon, et al., SIRVVD model-based verification of the effect of first and second doses of COVID-19/SARS-CoV-2 vaccination in Japan, Math. Biosci. Eng., 19 (2021), 1026–1040. https://doi.org/10.3934/mbe.2022047 doi: 10.3934/mbe.2022047
![]() |
[20] | N. Djenina, I. Rezzoug, A. Ouannas, T-E. Oussaeif, Giuseppe Grassi, A new COVID-19 pandemic model including the compartment of vaccinated individuals: Global stability of the disease-free fixed point, Submitted to CMMM, 2022 (2022). |
[21] |
P. van den Driesschea, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002). https://doi.org/10.1016/S0025-5564(02)00108-6 doi: 10.1016/S0025-5564(02)00108-6
![]() |
[22] |
S. Elaydi, An introduction to difference equations, Springer SBM, 3 (2005). https://doi.org/10.1007/0-387-27602-5 doi: 10.1007/0-387-27602-5
![]() |
1. | You-Kun Tai, Kuo-Ing Hwu, A Novel ZVS/ZCS Push-Pull LC Resonant DC-DC Converter for Energy Sources, 2023, 16, 1996-1073, 2892, 10.3390/en16062892 | |
2. | Mohammed Amine Hoummadi, Manale Bouderbala, Hala Alami Aroussi, Badre Bossoufi, Najib El Ouanjli, Mohammed Karim, Survey of Sustainable Energy Sources for Microgrid Energy Management: A Review, 2023, 16, 1996-1073, 3077, 10.3390/en16073077 | |
3. | Zaher Lamaouche, Abdelhamid Hamadi, Auguste Ndtoungou, Kettly Gustave, Kamal Al-Haddad, 2023, Cascade Droop-Virtual Synchronous Generator Control Approach to Improve Frequency Deviation, 979-8-3503-3182-0, 1, 10.1109/IECON51785.2023.10312210 |
Transient performance indicators | tpp (s) | tsp (s) | σp% | |
D constant | J = 10 | 1.5 | 4.8 | 44 |
J = 15 | 1.8 | 5.7 | 53 | |
J constant | D = 10 | 1.6 | 6.3 | 54 |
D = 20 | 1.8 | 4.3 | 25 |
Transient performance indicators | tpp (s) | tsp (s) | σp% | |
D constant | J = 10 | 0.7 | 5.2 | 0.30 |
J = 15 | 0.9 | 6.3 | 0.25 | |
J constant | D = 10 | 0.9 | 6.3 | 0.30 |
D = 20 | 0.7 | 4.4 | 0.23 |
Parameters | Value | Parameters | Value |
DC voltage Udc | 800 V | Line inductors Lg | 0.2 mH |
Rated voltage UN | 311 V | Line resistance Rg | 0.1 Ω |
Grid voltage Ug | 311 V | Rotational inertia J | 6 kg·m2 |
Rated frequency fo | 50 Hz | Damping factor D | 20 |
Filter inductors L1 | 5 mH | Active droop Kf | 20 |
Filter inductors L2 | 5 mH | Reactive droop Kq | 30 |
Filter capacitor C | 20 µF | Reactive power integration Ki | 50 |
Filter resistors R1 | 0.01 Ω | Active power reference Pref | 10 kW |
Filter resistors R2 | 0.01 Ω | Reactive power reference Qref | 0 Var |
Transient performance indicators | tpp (s) | tsp (s) | σp% | |
D constant | J = 10 | 1.5 | 4.8 | 44 |
J = 15 | 1.8 | 5.7 | 53 | |
J constant | D = 10 | 1.6 | 6.3 | 54 |
D = 20 | 1.8 | 4.3 | 25 |
Transient performance indicators | tpp (s) | tsp (s) | σp% | |
D constant | J = 10 | 0.7 | 5.2 | 0.30 |
J = 15 | 0.9 | 6.3 | 0.25 | |
J constant | D = 10 | 0.9 | 6.3 | 0.30 |
D = 20 | 0.7 | 4.4 | 0.23 |
Parameters | Value | Parameters | Value |
DC voltage Udc | 800 V | Line inductors Lg | 0.2 mH |
Rated voltage UN | 311 V | Line resistance Rg | 0.1 Ω |
Grid voltage Ug | 311 V | Rotational inertia J | 6 kg·m2 |
Rated frequency fo | 50 Hz | Damping factor D | 20 |
Filter inductors L1 | 5 mH | Active droop Kf | 20 |
Filter inductors L2 | 5 mH | Reactive droop Kq | 30 |
Filter capacitor C | 20 µF | Reactive power integration Ki | 50 |
Filter resistors R1 | 0.01 Ω | Active power reference Pref | 10 kW |
Filter resistors R2 | 0.01 Ω | Reactive power reference Qref | 0 Var |