Quasi-Z-Source Inverter-Based Photovoltaic Power System Modeling for Grid Stability Studies

: Quasi-Z-source inverters (qZSIs) are becoming a powerful power conversion technology in photovoltaic (PV) power systems because they allow energy power conversion in a single stage operation. However, they can cause system resonances and reduce system damping, which may lead to instabilities. These stability problems are well known in grid-connected voltage source converter systems but not in quasi-Z-source inverter (qZSI)-based PV power systems. This paper contributes with Matlab/Simulink and PSCAD/EMTDC models of qZSI-based PV power systems to analyze transient interactions and stability problems. These models consider all power circuits and control blocks of qZSI-based PV power systems and can be used in sensitivity studies on the inﬂuence of system parameters on stability. PV power system stability is assessed from the proposed models. The causes of instabilities are analyzed from numerical simulations and possible solutions are proposed.


Introduction
Penetration of wind and solar PV energies in grids has been rapidly growing in the last few years [1][2][3][4].PV power systems have different configurations achieved in a single or two stages [5,6], with the two-stage DC/DC and DC/AC voltage source inverter (VSI) topology being the most common configuration.Nevertheless, single stage topologies with a centralized inverter are gaining ground in a wide range of applications because they are simpler and cheaper than two stage topologies [6][7][8][9][10].Several impedance-source-networkbased power converters with different converter topologies that overcome the drawbacks of traditional converter topologies are reported in research works.A comprehensive review of them, together with their main topologies, modeling, control, and modulation techniques from an engineering application point of view are presented in [8,9].A complete analytical comparison of impedance-source-network-based power converters related to their passive component count and semiconductor stress, and a comparison between two and multilevel buck-boost inverters based on the already existing impedance-source networks are provided in [10].Recent solutions based on these networks have been extended to various fields of application such as Z-source inverters and quasi-Z-source inverters (qZSIs) for grid integration into low-voltage energy sources [10].In particular, introduction of qZSIs in PV power systems improves their efficiency and reliability because they can make both buck or boost DC voltage and energy power conversion in a single-stage operation [11][12][13][14][15][16][17].The traditional qZSI based on an impedance and semiconductor network coupling the inverter to the dc source is studied in [11][12][13][14].Four new qZSIs with several improvements and no disadvantages compared to ZSIs are presented and theoretically studied in [15].In addition, [16,17] examine a new qZSI topology, known as three-level neutral-point-clamped

•
Simulation tools, fully developed for implementation, based on Matlab/Simulink and PSCAD/EMTDC small-signal state-space averaged models of qZSI-based PV power systems which include all their components and controls.

•
An exhaustive study about the influence of qZSI-based PV power system parameters on stability.

•
Solutions to improve qZSI-based PV power system stability.
These contributions can be a useful tool and a guide on qZSI-based PV power system stability works for designers, engineers and researchers.

QZSI-Based PV Power Systems
The circuit and block diagram of qZSI-based PV power systems is illustrated in Figure 1.It consists of a DC bus where an N p × N s (number of PV panels in parallel × number of PV panels in series) PV array is connected with a shunt capacitor C p and a DC cable, characterized with a resistance R c , to the qZSI, which boosts the output voltage v dc at the VSI terminals.The maximum power point (MPP) tracking (MPPT) control is used to obtain the maximum power from the PV panel.The PV and current control loops ensure power flow from the PV panel to the AC grid.Moreover, the duty cycle control of the qZSI allows the DC-link peak voltage v dc,p to be fixed.
The Matlab/Simulink and PSCAD models of the qZSI-based PV power system are presented in the following Sections to assess system stability.

Matlab/Simulink Model of the qZSI-Based PV Power System
The Matlab/Simulink small-signal state-space averaged model of the qZSI-based PV power system in Figure 2 is derived.It is used to study the dynamic behavior of the PV power system in order to assess PV system stability and analyzing possible solutions to mitigate instability problems.The dynamic response of the PV power system is obtained by simulations, and the state-space matrix A and eigenvalues λi = σi ± j•ωi of the state-space modes of the PV power system can also be obtained from the Matlab/Simulink small-signal state-space averaged model of the qZSI-based PV power system.Moreover, the participation factors of each state-space variable of each mode can be calculated as pki = ϕki•ψik, where ϕki and ψik are the elements of the right and left eigenvectors of the state-space matrix A, respectively.This Section is divided into Subsections introducing the Matlab/Simulink model of the qZSI-based PV power system subsystems in Figure 2. In order to obtain these subsystem models, the small-signal state-space averaged equations of the circuits in Figure 1 are derived according to the state-space equation form where A, B, C, and D are the state-space, input, output, and feedthrough matrices, x is the state vector, u is the input vector, and y is the output vector.

PV Installation Modeling
The small-signal state-space model of the PV installation, formed by the PV panel with the MPPT control, the shunt capacitor Cp and the resistance Rc of the cable connecting the PV panel and the qZSI (see Figure 1), is derived.
The equation that relates the output current ipv and the terminal voltage vpv in the PV panel is expressed as [2], where IpvL and Ipv0 are the photovoltaic and the saturation currents of the PV cell, Rs is the equivalent series resistance of the PV cell, n is the diode quality factor and VT = k•T/q is the thermal constant of the PV cell, with k being the Boltzman constant (1.38 × 10 −23 J/K), q the Coulomb constant (1.602 × 10 −19 C), and T the PV cell temperature in Kelvin.The photovoltaic current IpvL is proportional to the irradiance level G and is usually referred to as the rated irradiation (i.e., G = 1 Sun = 1000 W/m 2 ).The photovoltaic and saturation currents, IpvL and Ipv0, also depend on the temperature T, which is 25 °C.

Matlab/Simulink Model of the qZSI-Based PV Power System
The Matlab/Simulink small-signal state-space averaged model of the qZSI-based PV power system in Figure 2 is derived.It is used to study the dynamic behavior of the PV power system in order to assess PV system stability and analyzing possible solutions to mitigate instability problems.The dynamic response of the PV power system is obtained by simulations, and the state-space matrix A and eigenvalues λ i = σ i ± j•ω i of the statespace modes of the PV power system can also be obtained from the Matlab/Simulink small-signal state-space averaged model of the qZSI-based PV power system.Moreover, the participation factors of each state-space variable of each mode can be calculated as p ki = φ ki •ψ ik , where φ ki and ψ ik are the elements of the right and left eigenvectors of the state-space matrix A, respectively.This Section is divided into Subsections introducing the Matlab/Simulink model of the qZSI-based PV power system subsystems in Figure 2. In order to obtain these subsystem models, the small-signal state-space averaged equations of the circuits in Figure 1 are derived according to the state-space equation form where A, B, C, and D are the state-space, input, output, and feedthrough matrices, x is the state vector, u is the input vector, and y is the output vector.

PV Installation Modeling
The small-signal state-space model of the PV installation, formed by the PV panel with the MPPT control, the shunt capacitor C p and the resistance R c of the cable connecting the PV panel and the qZSI (see Figure 1), is derived.
The equation that relates the output current i pv and the terminal voltage v pv in the PV panel is expressed as [2], where I pvL and I pv0 are the photovoltaic and the saturation currents of the PV cell, R s is the equivalent series resistance of the PV cell, n is the diode quality factor and V T = k•T/q is the thermal constant of the PV cell, with k being the Boltzman constant (1.38 × 10 −23 J/K), q the Coulomb constant (1.602 × 10 −19 C), and T the PV cell temperature in Kelvin.The photovoltaic current I pvL is proportional to the irradiance level G and is usually referred to as the rated irradiation (i.e., G = 1 Sun = 1000 W/m 2 ).The photovoltaic and saturation currents, I pvL and I pv0 , also depend on the temperature T, which is 25   Figure 3a shows the ipv-vpv characteristics (2) of a PV panel based on 60 W PV Solarex MSX60 modules [2].The influence of Np and Ns on the MPP is also illustrated in Figure 3a for three PV array configurations.
The small-signal model of the PV panel is obtained from the linearization of ipv-vpv characteristics (2) around the MPP (Vpv = Vmpp and Ipv = Impp) [5], where the incremental symbol ∆ identifies the small-signal variables and  Figure 3a shows the i pv -v pv characteristics (2) of a PV panel based on 60 W PV Solarex MSX60 modules [2].The influence of N p and N s on the MPP is also illustrated in Figure 3a for three PV array configurations.
The small-signal model of the PV panel is obtained from the linearization of i pv -v pv characteristics (2) around the MPP (V pv = V mpp and I pv = I mpp ) [5], where the incremental symbol ∆ identifies the small-signal variables and Energies 2021,  The small-signal circuit of the PV panel (3), shunt capacitor, and cable is shown in Figure 3b and the small-signal state-space model becomes The MPPT control generates the PV panel reference voltage by imposing the maximum value of the instantaneous power from the PV panel, ppv = ipv•vpv.Thus, the MPP can be mathematically characterized as [5] MPP MPP 0 0 .
According to (6), the sum of the actual conductance and the incremental conductance at the MPP (Vmpp, Impp) is equal to zero.The sum Σmpp is related to the PV voltage around the MPP by means of the linear I-V characteristics of the PV panel (3) as The small-signal circuit of the PV panel (3), shunt capacitor, and cable is shown in Figure 3b and the small-signal state-space model becomes The MPPT control generates the PV panel reference voltage by imposing the maximum value of the instantaneous power from the PV panel, p pv = i pv •v pv .Thus, the MPP can be mathematically characterized as [5] dp pv dv pv MPP = According to (6), the sum of the actual conductance and the incremental conductance at the MPP (V mpp , I mpp ) is equal to zero.The sum Σ mpp is related to the PV voltage around the MPP by means of the linear I-V characteristics of the PV panel (3) as where R mpp = V mpp/ I mpp is the incremental resistance of the linear i pv -v pv characteristics of the PV panel at the MPP (3).Considering the previous characterization of the MPP, the PV panel reference voltage is obtained from a PI control as the small-signal relation where k m p and k m i are the PI control proportional and integral gains of the MPPT control and k m = −2/(V mpp •R mpp ) [5].
According to ( 8), the small-signal state-space model of the MPPT control is represented by Figure 3c,d shows the influence of N p , N s and G on the MPP of the PV panel (i.e., on I mpp , V mpp and R mpp ).The MPP corresponding to the PV array configuration N p = 55 and N s = 42 and the operating point T 1 = 25 • C with G = 0.5 Sun and 0.75 Sun labeled with the blue and red dots, respectively.The current I mpp is affected by N p and by G, and the increase of both variables leads to the increase of I mpp .Variations of I mpp can have a large impact on the performance of the qZSI-PV-based power system, even leading the system to instability in some conditions (see Section 5).The voltage V mpp is mainly affected by N s , and the increase of this variable results in the increase of V mpp .The number of PV panels in series N s is usually fixed according to the rated voltage of the PV power system, and low variations of V mpp due to other causes (e.g., the irradiance level G) are controlled by means of the duty cycle control of the qZSI.The incremental resistance R mpp = V mpp /I mpp varies in accordance with the variations of V mpp and I mpp .

PV Panel Voltage Control and Grid-Connected VSI Modeling
The small-signal state-space model of the PV panel voltage control, VSI current control, and grid-connected VSI is presented in Figure 1.The grid d-frame reference current i dr is generated by the PI-based control of the PV panel voltage.The grid q-frame reference current i qr is set to zero considering that the VSI operates at unity power factor [6,8].The grid dq-frame reference voltage v dqr is provided by the PI-based control of the VSI.The VSI PLL is not considered in the study.
The grid d-reference current generated by the PV panel voltage control is expressed as the small-signal relation Energies 2021, 14, 508 7 of 16 where k pv p and k pv i are the PI control proportional and integral gains of the PV panel voltage control.According to (10), the small-signal state-space model of the PV panel voltage control is expressed as The output d-reference voltage generated by the PI compensator of the dq-frame VSC current control is expressed as the small-signal relation where k cc p and k cc i are the proportional and integral gains of the PI compensator.According to (12), the small-signal state-space model of the PI compensator is represented by The small-signal model of the grid-connected VSI is obtained from the d-current dynamics and the grid d-reference voltage v dr generated by the dq-frame VSI current control as where L f is the inductance of the converter filter.According to (14), the small-signal state-space model of the grid-connected VSI is represented by (15)

qZSI Modeling
This Subsection presents the small-signal state-space averaged model of the qZSI in Figure 1.The qZSI works in two different operational states during one switching cycle T: The shoot-through and the non-shoot-through states (see Figure 4) [11].The first interval is T 0 , where the inverter bridge is equivalent to short circuit.The second interval is T 1 , where the inverter is equivalent to a current source, which represents VSI consumption.The shoot-through interval is defined by the duty cycle d = T 0 /T and the non-shoot-through interval by T 1 /T = 1 − d.The qZSI control of the duty cycle is also plotted in Figure 1.This control allows the adjustment of the dc-link voltage v dc [11].

Power Circuit Model
The equations during each switching interval of the qZSI can be deduced from the equivalent circuits of the qZSI in Figure 4 [11].• Shoot-through state: The matrices of the state-space model in this interval are ( ) ( ) • Non-shoot-through state: The matrices of the state-space model in this interval are ( ) ( ) The small-signal state-space averaged model of the qZSI is finally expressed from ( 16) and ( 17) as [11]

•
Shoot-through state: The matrices of the state-space model in this interval are • Non-shoot-through state: The matrices of the state-space model in this interval are The small-signal state-space averaged model of the qZSI is finally expressed from ( 16) and ( 17) as [11] d dt where D is the steady state duty cycle and , and V C2 being the steady state values of the state variables.The qZSI input output current i dc is obtained from the DC/AC VSI power balance, which, assuming a lossless VSI, can be written as Energies 2021, 14, 508 9 of 16 The power balance ( 19) is expressed as the following small-signal relation: where m d0 = V d/ V dc is the steady-state operation point of the modulation function, G dc = 1/R dc = −P/V 2 dc , P is the active power flowing through the VSI from the DC to the AC side (see Figure 1), V d is the steady = state converter voltage, and V dc and I dc are the steady-state qZSI output voltage and current, respectively.The Matlab/Simulink model of this small-signal relation of the power balance is presented in Figure 2b.

Duty Cycle Control Model
The duty cycle control of the qSZI-based PV power system in Figure 1 is based on the adjustment of the dc-link peak voltage v dc,p through a PI-based control of the dc-link voltage reference V * dc, p together with a P-based control for the inductor-L 2 current.Finally, in order to smooth the duty cycle delivered to the VSI, a low-pass-filter (LPF) tuned to f c = 25 Hz is used.
The control laws of the dc-link voltage and inductor-L 2 current loops are expressed as the small-signal relation where k dc p and k dc i are the proportional and integral gains of the dc-link voltage PI-based control, and k L p is the proportional gain of the inductor-L 2 current P-based control.The small-signal averaged expression of the dc-link peak voltage is derived as [10] The small-signal state-space model of the duty cycle d r is derived from ( 21) and (22) as The small-signal state-space model of the LPF is represented by The Matlab/Simulink model of the duty cycle control is presented in Figure 2c.

PSCAD/EMTDC Model
PSCAD implementation of the qZSI-based PV power system is illustrated in Figure 5 according to the power and control systems of the qZSI-based PV power system in Figure 1 and the qZSI-based PV power system modeling in Section 3. The PSCAD model, composed of the power modules of the PV panel, qZSI, VSI, and AC grid, is shown in Figure 5a.Furthermore, Figure 5b presents the three control loops included in the qZSI-based PV power system, i.e., the PV voltage (left-top), current (right-top), and duty cycle (bottom) control loops.
according to the power and control systems of the qZSI-based PV power system in Figure 1 and the qZSI-based PV power system modeling in Section 3. The PSCAD model, composed of the power modules of the PV panel, qZSI, VSI, and AC grid, is shown in Figure 5a.Furthermore, Figure 5b

Applications
Simulations with the Matlab/Simulink and PSCAD/EMTDC models of the qZSI-based PV power system in Figure 1 and Table 1 data are shown in Figure 6a.The PV panel supplies a 138 kW 1 kV DC PV power system connected to a strong 400 V AC grid.Note that small-scale (Ppv < 10 kW) PV power systems usually work with DC voltages below 1 kV [7,13,16], but large-scale PV power systems (Ppv > 10 kW) tend to work with DC voltages in the range of 1 kV to 1.5 kV to reduce the initial PV array cost and total cabling required [7,11,19,25]

Applications
Simulations with the Matlab/Simulink and PSCAD/EMTDC models of the qZSIbased PV power system in Figure 1 and Table 1 data are shown in Figure 6a.The PV panel supplies a 138 kW 1 kV DC PV power system connected to a strong 400 V AC grid.Note that small-scale (P pv < 10 kW) PV power systems usually work with DC voltages below 1 kV [7,13,16], but large-scale PV power systems (P pv > 10 kW) tend to work with DC voltages in the range of 1 kV to 1.5 kV to reduce the initial PV array cost and total cabling required [7,11,19,25].In the application, the PV panel has 55 arrays connected in parallel, each of which consists of 42 PV Solarex MSX60 modules connected in series.Key data of these modules at 1 Sun are V oc = 21 V, I sc = 3.74 A, V max = 17.1 V, I max = 3.5 A, P max = 59.9 W [2].An irradiance level of 0.5 Sun and MPP voltage and current equal to V mpp = 16.73V and I mpp = 1.77A, which corresponds to a PV panel voltage, current and active power V pv = 16.73 × 42 = 702.9V, I pvs = 1.77 × 55 = 97.35A and P pv = 702.9× 97.35 = 68.2kW, respectively, are considered in the studied cases.It must also be noted that traditional IGBTs at the 10 kHz switching frequency f sw are set in the application, just as in [11,15,20].It is well known that the use of SiC MOSFETs or Si MOSFETs would allow the switching frequency to be increased up to the 25-65 kHz range to reduce passive component size and improve system efficiency [17,21,25].According to this, it would be interesting to further develop studies about the influence of switching frequency on dynamic behavior of qZSI-based PV power systems.
The Matlab/Simulink small-signal model of the qZSI-based PV power system model depends on the steady-state stable operation point values (i.e., V mpp = 702 V, I mpp = 94.2A, I L1 = I L2 = 94.1 A, V c1 = 760 V, V c2 = 63 V and D = 0.06), as obtained from PSCAD/EMTDC time domain simulations.

Circuit and Control
Parameters Data 0667 Ω, 10 mF qZSI source network  The dynamic response of the dc-link voltage vdc and the grid d-frame current id is studied when the PV panel current Ipvs is increased by 2 A and subsequently decreased by 1 A at 0.5 s and 1.5 s, respectively.The fair accuracy between the qZSI-based PV power system models is worth noting by comparing the dynamic response of the Matlab/Simulink and PSCAD/EMTDC simulations.It can also be noted that the PV power system is stable at all operating points.As an example, the PV power system eigenvalues of the operating point before 0.5 s, which is called reference operating point (ROP), are shown in Figure 6b.It is observed that the PV power system has no right hand plane (RHP) eigenvalues at the ROP.The same is true for the other two operating points.The dynamic response of the dc-link voltage v dc and the grid d-frame current i d is studied when the PV panel current I pvs is increased by 2 A and subsequently decreased by 1 A at 0.5 s and 1.5 s, respectively.The fair accuracy between the qZSI-based PV power system models is worth noting by comparing the dynamic response of the Matlab/Simulink and PSCAD/EMTDC simulations.It can also be noted that the PV power system is stable at all operating points.As an example, the PV power system eigenvalues of the operating point before 0.5 s, which is called reference operating point (ROP), are shown in Figure 6b.It is observed that the PV power system has no right hand plane (RHP) eigenvalues at the ROP.The same is true for the other two operating points.
The sensitivity of the circuit and control parameters on the qZSI-based PV power system eigenvalues and stability is numerically analyzed from the Matlab/Simulink smallsignal state-space averaged model (see results in Table 2).It can be observed that several parameters affect to some extent (i.e., low, medium, or high sensitivity) the locus eigenvalues but not all of them impact on system stability.In order to analyze this issue, Figure 7 shows the locus of the ROP eigenvalues with low, medium, and high sensitivity (labeled in blue and red as in Figure 6b) with respect to the circuit and control qZSI-based PV power system parameters.These parameters are varied according to 10% increments above and below the ROP circuit and control parameter values.The arrows indicate the growing value of the circuit and control parameters.The circuit and control parameters with no impact on the eigenvalues are omitted in Figure 7.It is observed that only the red eigenvalue is related to system instability because it could be shifted to the RHP by circuit and control parameter variations.Thus, only the circuit and control parameters that affect this eigenvalue can impact to some extent (low, medium, or high sensitivity) on system stability (i.e., C p , R c , L 1 , L 2 , C 1 , C 2 , I pvs and V * dc, p ).The sensitivity level depends on how the red eigenvalue approaches close to the RHP.It was numerically verified that only the external circuit parameter I pvs = 2•I mpp in Figure 5a can cause system instability.This is so because it is related to the active power P flowing through the VSI (large I pvs values lead to high P values) and this power affects system stability.The small-signal DC input current and voltage of the VSI, i.e., the qZSI output current and voltage i dc and v dc , are related by means of the virtual conductance G dc = -P/V 2 dc (20), which has a negative value for inverter operation of the VSI (i.e., for P > 0).In this situation, the VSI has a non-passive behavior at the DC side, which reduces the damping of qZSI-based PV power systems at resonances, and this can lead to system instability.This external parameter I pvs increases when the irradiance level G and the number of PV panels in parallel N p increase (see Figure 3c).Therefore, these factors increase the active power P flowing through the VSI and enlarge the value of G dc , compromising qZSI-based PV power system stability.It can be observed from Figure 7 that low values of L 2 , C 2 , k pv p or high values of L 1 , C 1 , R c , C p can move the system away from the RHP and help make the system more stable.It must also be noted that any factor that reduces the value of G dc in (20), such as a large steady state qZSI output voltage V dc , may also improve qZSI-based PV power system stability.For this reason, the increase of V * dc, p moves the system away from the RHP in Figure 7.Comparison of the above conclusions with the qZSI [9] and qZSI-based PV power system [22], ref. [23] stability studies in the literature shows that:

•
The influence of the inductors L 1 and L 2 and the capacitors C 1 and C 2 on qZSI dynamic behavior [9] can be extended to qZSI-based PV power system dynamic behavior: Inductors and capacitors have a high and low impact on the dynamic response, respectively.However, unlike [9], the present study differentiates the value of both inductors and concludes that the increase of L 1 may improve qZSI-based PV power system stability.

•
The detailed analysis of the present study shows the impact of the qZSI-based PV power system operating point on instabilities, in particular high and low values of the active power P flowing through the VSI and the steady-state qZSI output voltage V dc , respectively.The latter is directly related to the negative impact of high values of the duty cycle control in [22,23].

•
The impact of the AC network on qZSI-based PV power systems [22,23] is not studied in this paper because a strong AC grid is considered in the application.However, the proposed Matlab/Simulink and PSCAD/EMTDC small-signal state-space averaged models can be used to analyze this concern.
Energies 2021, 14, x FOR PEER REVIEW 14 of 17 both inductors and concludes that the increase of L1 may improve qZSI-based PV power system stability.

•
The detailed analysis of the present study shows the impact of the qZSI-based PV power system operating point on instabilities, in particular high and low values of the active power P flowing through the VSI and the steady-state qZSI output voltage Vdc, respectively.The latter is directly related to the negative impact of high values of the duty cycle control in [22,23].

•
The impact of the AC network on qZSI-based PV power systems [22,23] is not studied in this paper because a strong AC grid is considered in the application.However, the proposed Matlab/Simulink and PSCAD/EMTDC small-signal state-space averaged models can be used to analyze this concern.= +20%), and the increase of C1 from 3 mF to 3.9 mF (ΔC1 = +30%) restore system stability after the qZSI-based PV power system becomes unstable due to the step-up of the irradiance level.Although not shown for the sake of space, the influence of the other parameters (i.e., Np, C2, pv p k , L1, Rc and Cp) on stability was also validated.PSCAD/EMTDC simulations are shown in Figure 8 to verify the above conclusions.The dynamic response of the currents i L1 , i dc , and i d to G, L 2 , V * dc, p and C 1 variations is plotted.It can be noted that the increase of G from 0.5 Sun to 0.75 Sun causes system instability.On the other hand, the decrease of L 2 from 0.3 mH to 0.225 mH (∆L 2 = −25 %), the increase of V * dc, p from 800 V to 960 V (∆V * dc, p = +20%), and the increase of C 1 from 3 mF to 3.9 mF (∆C 1 = +30%) restore system stability after the qZSI-based PV power system becomes unstable due to the step-up of the irradiance level.Although not shown for the sake of space, the influence of the other parameters (i.e., N p , C 2 , k pv p , L 1 , R c and C p ) on stability was also validated.

Conclusions
This first part of the paper contributes to qZSI-based PV power system stability studies with two simulation tools based on a Matlab/Simulink and a PSCAD model of qZSI-based PV power systems.Both models consider all the main components and controls.They make it possible to analyze system stability based on system eigenvalues and dynamic numerical simulations, respectively.The second part of the paper (Section 5) contributes with the stability study performed from the above models and the discussion about the impact of parameter variations on system stability.Direction on the qZSI-based PV power system stability is provided to designers, engineers, and researchers.The following recommendations and solutions are proposed:

•
The operating point is the main factor related to system stability.In particular, o The increase of active power flowing through the VSI, which directly depends on the irradiance level and the number of PV panels in parallel, can lead to system stability loss.o The dc-link peak voltage also affects system stability.According to this, the increase of the dc-link peak voltage reference may be used to enhance stability because it mitigates the impact of active power flowing through the VSI.

•
The qZSI component design is important for upgrading PV power system stability.In particular, low values of L2, C2 and high values of L1, C1 can help to make the system more stable.

•
Other PV power system parameters such as the PV panel shunt capacitor Cp and the DC cable Rc could help to improve stability if they take high values.

Conclusions
This first part of the paper contributes to qZSI-based PV power system stability studies with two simulation tools based on a Matlab/Simulink and a PSCAD model of qZSI-based PV power systems.Both models consider all the main components and controls.They make it possible to analyze system stability based on system eigenvalues and dynamic numerical simulations, respectively.The second part of the paper (Section 5) contributes with the stability study performed from the above models and the discussion about the impact of parameter variations on system stability.Direction on the qZSI-based PV power system stability is provided to designers, engineers, and researchers.The following recommendations and solutions are proposed:

•
The operating point is the main factor related to system stability.In particular, The increase of active power flowing through the VSI, which directly depends on the irradiance level and the number of PV panels in parallel, can lead to system stability loss.The dc-link peak voltage also affects system stability.According to this, the increase of the dc-link peak voltage reference may be used to enhance stability because it mitigates the impact of active power flowing through the VSI.

•
The qZSI component design is important for upgrading PV power system stability.In particular, low values of L 2 , C 2 and high values of L 1 , C 1 can help to make the system more stable.

•
Other PV power system parameters such as the PV panel shunt capacitor C p and the DC cable R c could help to improve stability if they take high values.
• PV installation, qZSI, PV control, and D control parameters have little influence on system stability.

Figure 3 .
Figure 3. Study of the PV panel (T1 = 25 °C): (a) ipv-vpv characteristics; (b) small-signal equivalent circuit; (c) influence of Np and G on MPP; (d) influence of Ns and G on MPP.

Figure 3 .
Figure 3. Study of the PV panel (T 1 = 25 • C): (a) i pv -v pv characteristics; (b) small-signal equivalent circuit; (c) influence of N p and G on MPP; (d) influence of N s and G on MPP.

Figure 4 .
Figure 4. Equivalent circuit of qZSI with continuous current input: (a) Shoot-through state; (b) non-shoot-through state.

Figure 4 .
Figure 4. Equivalent circuit of qZSI with continuous current input: (a) Shoot-through state; (b) non-shoot-through state.
. In the application, the PV panel has 55 arrays connected in parallel, each of which consists of 42 PV Solarex MSX60 modules connected in series.Key data of these modules at 1 Sun are Voc = 21 V, Isc = 3.74 A, Vmax = 17.1 V, Imax = 3.5 A, Pmax = 59.9 W [2].An irradiance level of 0.5 Sun and MPP voltage and current equal to Vmpp = 16.73V and Impp = 1.77A, which corresponds to a PV panel voltage, current and active power Vpv = 16.73 × 42 = 702.9V, Ipvs = 1.77 × 55 = 97.35A and Ppv = 702.9× 97.35 = 68.2kW, respectively, are considered in the studied cases.It must also be noted that traditional

Figure 6 .
Figure 6.Study of the qZSI-based PV power system with the Matlab/Simulink and PSCAD/EMTDC models: (a) Time domain simulations.(b) Eigenvalues.

Figure 6 .
Figure 6.Study of the qZSI-based PV power system with the Matlab/Simulink and PSCAD/EMTDC models: (a) Time domain simulations.(b) Eigenvalues.

Figure 7 .
Figure 7. Root locus of the qZSI-based PV power system eigenvalues.PSCAD/EMTDC simulations are shown in Figure 8 to verify the above conclusions.The dynamic response of the currents iL1, idc, and id to G, L2, * , dc p V and C1 variations is plotted.It can be noted that the increase of G from 0.5 Sun to 0.75 Sun causes system instability.On the other hand, the decrease of L2 from 0.3 mH to 0.225 mH (ΔL2 = −25 %), the increase of * , dc p V from 800 V to 960 V (Δ * , dc p V = +20%), and the increase of C1 from 3 mF

Figure 7 .
Figure 7. Root locus of the qZSI-based PV power system eigenvalues.

Table 2 .
Sensitivity analysis of the qZSI-based PV power system eigenvalues and stability to circuit and control parameters.