Small-Signal Modeling and Parameter Optimization Design for Photovoltaic Virtual Synchronous Generator

: With the continuous proliferation of renewable energy generation, distributed photovoltaic inverters operating at a maximum power point reduce the inertia of power systems, degrading system frequency stability and potentially causing severe oscillations in systems after being disturbed. The virtual synchronous generator (VSG) control method, which causes photovoltaic inverters to possess inertia and damping, now plays an important role in the ﬁeld of distributed generation. However, while introducing the advantages of synchronous machines, problems with oscillations are also introduced and the stochastic ﬂuctuation characteristic of photovoltaics results in the stochastic drifting of the operating point. This paper presents an adaptive controller parameter design method for a photovoltaic-VSG (PV-VSG) integrated power system. Firstly, a small-signal model of the PV-VSG is built and a state space model is deduced. Then, the small-signal stability and low frequency oscillation characteristics of the photovoltaic power generation system are analyzed. Finally, considering the limitations of system oscillations and the stochastic drifting of the operating point, a global optimization design method for controller parameters used to improve system stability is proposed. The time domain simulation shows that an optimized PV-VSG could provide su ﬃ cient damping in the case of photovoltaic power output changes across a wider range.


Introduction
Resource shortages around the world are becoming more and more serious, and the development of renewable power generation technologies is extremely urgent. As a clean and renewable energy source, solar power plays a more and more important role in China's energy strategy and its installed capacity is increasing day by day [1,2]. However, since photovoltaic inverters based on power electronic inverter interfaces have almost no rotational kinetic energy, the inertia of these systems will also decrease when they are connected to the power grid on a large scale. When the system power fluctuates, the frequency may oscillate or even cause the frequency stability to collapse, which seriously threatens the stable operation of the system. The impact of large-scale solar power access on grid stability will be further exposed and deteriorated, and system stability faces unprecedented challenges [3][4][5].
As an important factor affecting the dynamic characteristics of the frequency, system inertia has become a key point of research regarding the stable operation of a high-permeability photovoltaic microgrid. The low inertia caused by photovoltaic inverters is gaining more and more attention [6][7][8].
In order to improve the frequency dynamic stability of photovoltaic systems, it is particularly important to improve the control performance of the inverter. When a power system is in the electromechanical transient transition process, the mechanical inertia of the synchronous generator rotor provides transient energy, suppresses the change in the grid frequency, and strives for a response time for regulation, which plays a vital role in grid stability. Following the principle of its function, the idea of virtual synchronous generator control for renewable energy has been proposed by scholars and has received extensive attention in academic circles and industry in recent years [9][10][11][12]. The virtual synchronous generator (VSG) introduces a rotor equation to make the inverter similar to a traditional synchronous generator in both the operating mechanism and the external characteristics, which can provide significant inertial support for the system [13,14].
The characteristics of the VSG meet the needs of future grids for distributed photovoltaic (PV) power access. Thus, virtual synchronous machine technology provides a new idea and technical route for solving many problems brought about by large-scale photovoltaic power supply access to power systems. In the existing literature, virtual synchronous machine topology and control structures are mainly used for cases in which the DC side is an energy storage battery, that is, the virtual prime mover simulated on the DC side is a kind of "power source with constant voltage and unlimited capacity". The DC source of a traditional VSG is mainly used for energy storage batteries, which ignores the influence of DC side dynamic characteristics to some extent [9,12]. In general, distributed photovoltaic power output is susceptible to external environmental factors and has strong randomness, volatility, and finiteness [15][16][17]. Hence, the traditional VSG has a limited application range and is difficult to apply to distributed PV access.
For this reason, the PV-VSG has been proposed by scholars as a way to apply virtual synchronous machine control to photovoltaic inverters to cause them to have better inertia and damping [18][19][20][21]. However, most of existing studies mainly simulate different PV-VSG control strategies and do not analyze the PV-VSG grid-connected system from the perspective of system stability. The parameter design is mainly designed for the operating point at a single output level. In [22], which combs VSG technology and considers the dynamic characteristics of a primary source, a PV-VSG and its control strategies are proposed. However, this work only simulates the effect of a PV-VSG and does not analyze the stability of the system and the selection of controller parameters. Virtual synchronous machine control allows the photovoltaic power station to participate in power grid frequency modulation and improve the system frequency stability. At the same time, it opens the channel that interacts with the active power of the system under disturbance, and has a non-negligible influence on the oscillation of the power system [23,24].
In this paper, an accurate small-signal model of a PV-VSG is established and the influence of controller parameters on the small-signal stability of the system is analyzed, and based on this, a globally optimized parameter design method is proposed. Firstly, with full consideration of a detailed model of each important link in the photovoltaic power generation system, a small-signal model of a PV-VSG that is suitable for transient stability analysis is established. Then, the small-signal stability and low frequency oscillation characteristics of the photovoltaic power generation system are analyzed. Finally, considering the limitations of system oscillations and the stochastic drifting of the operating point, a global optimization design method for controller parameters used to improve system stability is proposed. A time domain simulation shows that the optimized PV-VSG could provide sufficient damping in the case of photovoltaic power output changes across a wider range.

Small-Signal Model of a PV-VSG Integrated Power System
The photovoltaic power generation system consists of photovoltaic cell components, energy storage capacitors, inverters, and controllers. The photovoltaic power generation system structure is shown in Figure 1.

PV Cell Model
The exact model of the photovoltaic cell module can be established according to the theory of electronics, but the exact model is very complicated, and some parameters are difficult to measure directly, which is not convenient for research and applications. The traditional practical engineering model is utilized as the PV cell model. The model uses four factory electrical parameters of the photovoltaic cell, namely, short-circuit current where e is the natural logarithmic base, approximately 2.71828. The compensation coefficients  ,  , and  are 0.0025 °C, 0.5 2 m W , and 0.00288 °C, respectively.

DC Link Model
The mathematical model of the intermediate capacitor is as shown in Equation (10),

PV Cell Model
The exact model of the photovoltaic cell module can be established according to the theory of electronics, but the exact model is very complicated, and some parameters are difficult to measure directly, which is not convenient for research and applications. The traditional practical engineering model is utilized as the PV cell model. The model uses four factory electrical parameters of the photovoltaic cell, namely, short-circuit current I sc0 , open circuit voltage U oc0 , and current I m0 and voltage U m0 , when the photovoltaic cell obtains maximum power under standard conditions and test parameters (S 0 = 1000W/m 2 , T 0 = 25°C). The U-I equation, Equation (1), under non-standard conditions is The short-circuit current I sc , open-circuit voltage U oc , and current at maximum power I m and voltage U m under non-standard conditions can be obtained by Equations (4)-(9), i.e., where e is the natural logarithmic base, approximately 2.71828. The compensation coefficients α, β, and γ are 0.0025 • C, 0.5 m 2 /W, and 0.00288 • C, respectively.

DC Link Model
The mathematical model of the intermediate capacitor is as shown in Equation (10), where C dc is the capacitance value of the intermediate capacitor, V dc is the DC voltage of the capacitor, I pv is the current from the PV cells to the intermediate capacitor, and P g is the active power from the intermediate capacitor to the grid-side converter.

Inverter and Controller Model
In the photovoltaic power generation system, the topology of the voltage-type three-phase PWM (Pulse Width Modulation) inverter controlled by sinusoidal pulse width modulation (SPWM) is shown in Figure 2. Considering that the research object of this paper is mainly the electromechanical frequency band, the high frequency dynamic process in the system was ignored to simplify the research.

Inverter and Controller Model
In the photovoltaic power generation system, the topology of the voltage-type three-phase PWM (Pulse Width Modulation) inverter controlled by sinusoidal pulse width modulation (SPWM) is shown in Figure 2. Considering that the research object of this paper is mainly the electromechanical frequency band, the high frequency dynamic process in the system was ignored to simplify the research.
-Reset  . θ g , the frequency of the equivalent internal potential of VSG; ω n , the base frequency of the grid; P re f and Q re f , the set values of active and reactive which are used to control the output power of the inverter; J, the inertia of the virtual synchronous machine; D, the damping coefficient; K, the gain factor of the reactive controller.

The dynamic equation of the VSG control method is
where P e and Q e are the output active and reactive powers of the inverter, and the calculations are In a traditional virtual synchronous generator control, the virtual synchronous machine topology and the control structure are mainly used for the cases in which the DC side is an energy storage battery, that is, the virtual prime mover simulated on the DC side is a kind of "power source with constant voltage and unlimited capacity". However, the dynamic characteristics of photovoltaic power supply are more complicated. The available output range of photovoltaic power supply is [0, P max ], and there are stable operation areas and unstable operation areas. The stable operation area of the photovoltaic power supply can also be determined by the DC voltage, that is, the DC voltage (V dc−mpp ) corresponding to the MPP (Maximum Power Point) point is the minimum DC voltage that can be achieved by the photovoltaic power supply. If the DC voltage continues to decrease, it will travel to the unstable region. Thus, in order to fully consider the dynamic characteristics of the photovoltaic power supply, and regardless of the external power demand, it should be guaranteed to operate in a stable area.
As is shown in Figure 3, V dc−mpp is the maximum power point and P max is the corresponding output power which could be obtained by the maximum power tracing algorithm. Considering the stable area of the PV-VSG, the upper limit of ∆X is set as 0, and thus the P re f will lower than P max , and the PV-VSG operate within a stable area. The control equation of the P re f and V re f can be obtained as The dynamic equation is where K iv is the integral coefficient of the PI controller and V dc−mpp is the maximum power point of the photovoltaic component.

The Linearization Model
The PV-VSG dynamic model as given in Equations (1)- (14) is linearized to obtain a small-signal model of a photovoltaic battery with VSG control.

The Linearization Model
The PV-VSG dynamic model as given in Equations (1)- (14) is linearized to obtain a small-signal model of a photovoltaic battery with VSG control.
In Equation (15), [∆U gx , ∆U gy ] T is the PV-VSG bus voltage, and ∆I P = [∆I x , ∆I y ] T is the PV-VSG bus injection current. The matrices in the formula are where the elements in the matrix are By combining Equations (15)-(23), a small-signal state space model of the system shown in Figure 1 can be obtained, i.e., where

is the input vector; and
are the state matrix, input matrix, output matrix, and direct transfer matrix, respectively. Equation (24) is a linearized model of a PV-VSG grid-connected system.

Analysis of Eigenvalues of PV-VSG Systems
The system shown in Table 1 is taken as an example for small-signal analysis of the PV-VSG. The initial conditions of the PV battery are S = 1000, W/m 2 , and T = 25 • C.
According to the linearization model of the PV-VSG, small-signal stability analysis was performed on the test system and the eigenvalues of the state matrix are shown in the Table 2. As can be seen from the table, the damping of modes 2-3 is small and subsynchronous oscillation is prone to occur. Modes 4-5 are low-frequency oscillation modes. Table 2. Eigenvalues of the studied PV system.

Mode
Eigenvalue Frequency/Hz Damping Ratio In order to analyze the stability range of the system and the variation of the eigenvalues, the eigenvalue trajectories of the PV-VSG generation system have been calculated, as shown in Figure 4.
(1) Active power loop parameters. The active power loop parameters include D and J, which simulate the characteristics of the generator and provide inertia and damping for the system. It can be seen from Figure 4a that the change in D mainly affects λ 2,3 and has little effect on the other modes. As D increases, λ 2,3 moves toward the stable area and the imaginary part decreases significantly. The change in modes λ 1 and λ 4,5 is not obvious, so the increase in D is beneficial to the attenuation of the mode. It can be seen from Figure 4b that the change in J mainly affects λ 2,3 ,λ 4,5 , and has little effect on the other modes. As J increases, λ 2,3 moves toward the unstable area and the imaginary part increases significantly. The change in mode λ 1 is not obvious. However, the increase also leads to the movement of mode λ 4,5 towards the stable area, but this is less obvious. Thus, the increase in J should be limited to a certain range.
(2) Voltage regulation parameters. The voltage regulation parameters include D q and K, which simulate the characteristics of the excitation regulation and adjust the voltage at the point of interconnection. It can be seen from Figure 4c that the change in D q mainly affects λ 1 . As D increases, λ 1 moves toward the stable area. The imaginary part of mode λ 2,3 decreases slightly and the change in λ 4,5 is not obvious. It can be seen from Figure 4d that the change in K mainly affects λ 1 and λ 2,3 . As K increases, mode λ 1 moves toward the unstable area and mode λ 2,3 moves towards the stable area. The change in mode λ 4,5 is not obvious.
(3) PV power regulation parameter. The voltage regulation parameters include K iv and K pv , which generate the P re f based on the output characteristics of the photovoltaic battery and keep the photovoltaic generation system working around the maximum power point. It can be seen from Figure 4e that the change in P iv mainly affects λ 4,5 . As P iv increases, mode λ 4,5 gradually moves from the unstable region toward the stable region, and then the imaginary part begins to increase. It can be seen from Figure 4f that the change in P kv mainly affects λ 2,3 and λ 4,5 . As P kv increases, the imaginary part of the mode λ 4,5 gradually decreases and then begins to move from the stable region toward the unstable region. Mode λ 2,3 gradually moves from the unstable area toward the stable area.

Optimization Method for Controller Parameters of PV-VSG
In order to improve the stability and stability margin of the photovoltaic power generation system, controller parameters can be designed. According to the influence law of the characteristic value on the transient characteristics of the system, the parameter selection of the controller should make the absolute value of the real part of the eigenvalue as large as possible to shorten the transient process and return it to the steady state. Considering that the frequency of the power system is 50 Hz, it is desirable that the oscillation frequency of the oscillation mode be much smaller than 50 Hz. In order to reduce the number of oscillation cycles, it is desirable that the damping of each mode be relatively large.
The traditional controller parameter design method is to determine each controller parameter according to the analysis of the eigenvalue trajectory of the system. However, the traditional design method lacks scientific and reasonable quantitative design standards and cannot provide an optimal parameter combination in the global scope for the characteristics of system parameter coupling. Based on the general rules of parameter design, this paper utilizes a quantitative model to reflect the comprehensive benefits of different parameter combinations. In summary, the control parameters should be chosen such that the real part of the eigenvalue is negative and its absolute value is large; the absolute value of the eigenvalue's imaginary part is small. Therefore, the objective function of control parameter optimization is where i  is the i th eigenvalue of the system; Re( i  ) and Im( i  ) are the real and imaginary parts of i  , respectively; and i w and i p are the weights of the real part of i  , respectively.
The initial conditions of the photovoltaic power generation system under different lighting conditions and temperature are different, which affects the accuracy of the small-signal model. Thus, the linearization is performed under different initial conditions. The objective function of the control parameter optimization is

Optimization Method for Controller Parameters of PV-VSG
In order to improve the stability and stability margin of the photovoltaic power generation system, controller parameters can be designed. According to the influence law of the characteristic value on the transient characteristics of the system, the parameter selection of the controller should make the absolute value of the real part of the eigenvalue as large as possible to shorten the transient process and return it to the steady state. Considering that the frequency of the power system is 50 Hz, it is desirable that the oscillation frequency of the oscillation mode be much smaller than 50 Hz. In order to reduce the number of oscillation cycles, it is desirable that the damping of each mode be relatively large.
The traditional controller parameter design method is to determine each controller parameter according to the analysis of the eigenvalue trajectory of the system. However, the traditional design method lacks scientific and reasonable quantitative design standards and cannot provide an optimal parameter combination in the global scope for the characteristics of system parameter coupling. Based on the general rules of parameter design, this paper utilizes a quantitative model to reflect the comprehensive benefits of different parameter combinations. In summary, the control parameters should be chosen such that the real part of the eigenvalue is negative and its absolute value is large; the absolute value of the eigenvalue's imaginary part is small. Therefore, the objective function of control parameter optimization is where λ i is the i th eigenvalue of the system; Re(λ i ) and Im(λ i ) are the real and imaginary parts of λ i , respectively; and w i and p i are the weights of the real part of λ i , respectively. The initial conditions of the photovoltaic power generation system under different lighting conditions and temperature are different, which affects the accuracy of the small-signal model. Thus, the linearization is performed under different initial conditions. The objective function of the control parameter optimization is where U j is the objective function of the control parameter optimization in the j th initial condition. In order to ensure the stability of the system, the value of the controller parameters and the real part of the eigenvalue must be limited, i.e., where J max , K max , D qmax , K ivmax , and K pvmax are the upper limit stability thresholds which can be obtained according to the eigenvalue trajectories of the PV-VSG generation system. Considering system dynamic performance, D max can be obtained by maintaining the damping of the open loop transfer function of the VSG. C th is the stable threshold of the real part of the eigenvalue, which is a negative number, and the larger the absolute value, the more stable the stability of the system. The controller parameter optimization design goal is transformed into searching for the optimal combination of all parameters satisfying the constraint within the entire parameter range, so that the objective function U is maximized.

Simulation and Results
In the test system, D max , J max , K max , D qmax , K ivmax , K pvmax , and C th are taken to be 30, 10, 20, 10, 20, 20, and −2, respectively. The optimal combination of parameters according to the global optimal design method is shown in Table 3, and the corresponding system eigenvalues are shown in Table 4.  Table 4. Eigenvalues of PV-VSG using initial parameters and optimized parameters. To further demonstrate the effectiveness of the proposed parameter design method, a time domain simulation has been performed. A contrast method in [25] is utilized to highlight the effectiveness of the proposed method. The comparison method is designed based on the small-signal model when the light intensity is 750.

Frequency/Hz Damping Ratio
Case A: First let the system work at the operating point where the light intensity is 1000; at this time the PV-VSG output active power is 1 p.u. At 1 s, the light intensity of the PV-VSG is increased by 5% and reverted back to the initial values 20 ms later. Figure 5 shows the active power of the PV-VSG using the initial parameters, the proposed method, and the contrast method. Because the robustness has been considered in the parameter design, it is not surprising that the damping ratios of the proposed method are considerably better than those of the other method.
Energies 2020, 13, 398 11 of 15 VSG using the initial parameters, the proposed method, and the contrast method. Because the robustness has been considered in the parameter design, it is not surprising that the damping ratios of the proposed method are considerably better than those of the other method.  Case B: First let the system work at the operating point where the light intensity is 1000; at this time the PV-VSG output active power is 1 . . At 1 s, the light intensity of the PV-VSG is increased by 5% and reverted back to the initial values at 5 s. Figure 6 shows the active power of the PV-VSG using the initial parameters, the proposed method, and the contrast method. It can be seen that after the application of the optimal parameters, the active power is quickly attenuated and the PV-VSG switches to the new operating point with a faster speed and less oscillation. Case C: The fluctuation of solar power output is used to test the performance of the proposed method. At 1 s, the light intensity of PV-VSG starts to change randomly by 5%. In Figure 7, the blue, red, and black curves represent the dynamic active powers and DC voltage of the PV-VSG with the initial parameters, with the proposed method, and with the contrast method, respectively. From the figure, it can be observed that the red curves have smaller oscillation magnitudes than the blue and black curves, proving the increased robustness of the proposed method. Case B: First let the system work at the operating point where the light intensity is 1000; at this time the PV-VSG output active power is 1 p.u. At 1 s, the light intensity of the PV-VSG is increased by 5% and reverted back to the initial values at 5 s. Figure 6 shows the active power of the PV-VSG using the initial parameters, the proposed method, and the contrast method. It can be seen that after the application of the optimal parameters, the active power is quickly attenuated and the PV-VSG switches to the new operating point with a faster speed and less oscillation.
Energies 2020, 13, 398 11 of 15 VSG using the initial parameters, the proposed method, and the contrast method. Because the robustness has been considered in the parameter design, it is not surprising that the damping ratios of the proposed method are considerably better than those of the other method. Case B: First let the system work at the operating point where the light intensity is 1000; at this time the PV-VSG output active power is 1 . . At 1 s, the light intensity of the PV-VSG is increased by 5% and reverted back to the initial values at 5 s. Figure 6 shows the active power of the PV-VSG using the initial parameters, the proposed method, and the contrast method. It can be seen that after the application of the optimal parameters, the active power is quickly attenuated and the PV-VSG switches to the new operating point with a faster speed and less oscillation. Case C: The fluctuation of solar power output is used to test the performance of the proposed method. At 1 s, the light intensity of PV-VSG starts to change randomly by 5%. In Figure 7, the blue, red, and black curves represent the dynamic active powers and DC voltage of the PV-VSG with the initial parameters, with the proposed method, and with the contrast method, respectively. From the figure, it can be observed that the red curves have smaller oscillation magnitudes than the blue and black curves, proving the increased robustness of the proposed method. Case C: The fluctuation of solar power output is used to test the performance of the proposed method. At 1 s, the light intensity of PV-VSG starts to change randomly by 5%. In Figure 7, the blue, red, and black curves represent the dynamic active powers and DC voltage of the PV-VSG with the initial parameters, with the proposed method, and with the contrast method, respectively. From the figure, it can be observed that the red curves have smaller oscillation magnitudes than the blue and black curves, proving the increased robustness of the proposed method.  Case D: At 1 s, the light intensity of PV-VSG decreases from 1000 to 300. Then, the light intensity increases from 300 to 750 at 4 s. The parameter design method proposed in this paper takes into account the difference in operating points caused by light intensity and temperature. Great damping is achieved for different output levels and is not designed for a single operating point. The time domain simulation results using the initial parameters, the proposed method, and the contrast method are shown in Figure 8. When the output changes and the system operating point is far away from the initial operating point of the initial parameter controller design, the proposed method has a better suppression effect on oscillation than the other method. Case E: A three-phase temporary fault is applied at the bus of PV-VSG at t = 1 s and is cleared 20 ms later. The dynamic responses of the PV-VSG are shown in Figure 9. The four blue, red, and black curves represent the dynamic active powers and DC voltage of the PV-VSG with the initial parameters, with the proposed method, and with the contrast method, respectively. It can be seen from the figure that the proposed method gives a better damping performance under fault conditions. Case D: At 1 s, the light intensity of PV-VSG decreases from 1000 to 300. Then, the light intensity increases from 300 to 750 at 4 s. The parameter design method proposed in this paper takes into account the difference in operating points caused by light intensity and temperature. Great damping is achieved for different output levels and is not designed for a single operating point. The time domain simulation results using the initial parameters, the proposed method, and the contrast method are shown in Figure 8. When the output changes and the system operating point is far away from the initial operating point of the initial parameter controller design, the proposed method has a better suppression effect on oscillation than the other method.
Energies 2020, 13  Case D: At 1 s, the light intensity of PV-VSG decreases from 1000 to 300. Then, the light intensity increases from 300 to 750 at 4 s. The parameter design method proposed in this paper takes into account the difference in operating points caused by light intensity and temperature. Great damping is achieved for different output levels and is not designed for a single operating point. The time domain simulation results using the initial parameters, the proposed method, and the contrast method are shown in Figure 8. When the output changes and the system operating point is far away from the initial operating point of the initial parameter controller design, the proposed method has a better suppression effect on oscillation than the other method. Case E: A three-phase temporary fault is applied at the bus of PV-VSG at t = 1 s and is cleared 20 ms later. The dynamic responses of the PV-VSG are shown in Figure 9. The four blue, red, and black curves represent the dynamic active powers and DC voltage of the PV-VSG with the initial parameters, with the proposed method, and with the contrast method, respectively. It can be seen Case E: A three-phase temporary fault is applied at the bus of PV-VSG at t = 1 s and is cleared 20 ms later. The dynamic responses of the PV-VSG are shown in Figure 9. The four blue, red, and black curves represent the dynamic active powers and DC voltage of the PV-VSG with the initial parameters, with the proposed method, and with the contrast method, respectively. It can be seen from the figure that the proposed method gives a better damping performance under fault conditions. It is clear from Figure 9 that the proposed method damps oscillations within 200 ms, thereby confirming the robustness of the proposed methodology under fault conditions. Energies 2020, 13,398 13 of 15 It is clear from Figure 9 that the proposed method damps oscillations within 200 ms, thereby confirming the robustness of the proposed methodology under fault conditions.

Conclusions
For the PV-VSG integrated power system, this paper presents an adaptive controller parameter design method. A small-signal model of a PV-VSG is built and the small-signal stability and oscillation characteristics of the photovoltaic power generation system are analyzed. Considering the limitations of system oscillations and the stochastic drifting of the operating point, controller parameters are designed using a global optimization method to improve system stability. The optimized parameters are used to maintain a good damping performance when the output of the PV cells randomly drifts at the operating point. The time domain simulation results of the PV-VSG integrated power system prove the effectiveness and adaptability of the proposed parameter design method.
In this paper, the power system into which PV-VSG is incorporated is an infinite system. In an actual system, the interaction between PV-VSG and other dynamic components in the system such as generators and HVDC (High Voltage Direct Current), etc., may cause oscillation. Studying the dynamic interaction between PV-VSG and complex systems and the suppression of oscillations will be the focus of future research.

Conflicts of Interest:
No conflict of interest exists in the submission of this manuscript and the manuscript has been approved by all authors for publication. We would like to declare on behalf of my co-authors that the work described was original research that has not been published previously and is not under consideration for publication elsewhere, in whole or in part.

Conclusions
For the PV-VSG integrated power system, this paper presents an adaptive controller parameter design method. A small-signal model of a PV-VSG is built and the small-signal stability and oscillation characteristics of the photovoltaic power generation system are analyzed. Considering the limitations of system oscillations and the stochastic drifting of the operating point, controller parameters are designed using a global optimization method to improve system stability. The optimized parameters are used to maintain a good damping performance when the output of the PV cells randomly drifts at the operating point. The time domain simulation results of the PV-VSG integrated power system prove the effectiveness and adaptability of the proposed parameter design method.
In this paper, the power system into which PV-VSG is incorporated is an infinite system. In an actual system, the interaction between PV-VSG and other dynamic components in the system such as generators and HVDC (High Voltage Direct Current), etc., may cause oscillation. Studying the dynamic interaction between PV-VSG and complex systems and the suppression of oscillations will be the focus of future research.