^{1}

^{1}

^{1}

In recent years, power generation using renewable energy sources has been increasing as a solution to the global warning problem. Wind power generation can generate electricity day and night, and it is relatively more efficient among the renewable energy sources. The penetration level of variable-speed wind turbines continues to increase. The interconnected wind turbines, however, have no inertia and no synchronous power. Such devices can have a serious impact on the transient stability of the power grid system. One solution to stabilize such grid with renewable energy sources is to provide emulated inertia and synchronizing power. We have proposed an optimal design method of current control for virtual synchronous generators. This paper proposes an optimal control method that can follow the virtual generator model under constrains. As a result, it is shown that the proposed system can suppress the peak of the output of semiconductor device under instantaneous output voltage drop.

At present day, the introduction of renewable energy that can suppress global warming and gas emissions that causes global warming is progressing. Renewable energy includes “solar power”, “wind power”, “ocean energy”, “hydropower”, “geothermal”, and “biomass”, etc. In particular, wind power can generate electricity day and night and it is a relatively efficient power generation method among the renewable sources. According to the Global Wind Report 2018, the installed amount of wind power in the world was 540 GW in 2017 and 591 GW in 2018, and the installed amount is increasing [

In [

In [

There are many studies applying such advanced control methods to efficiently obtain wind power.

Generally, synchronous generators are used in power systems to supply electrical energy. The inertia and synchronization power of the synchronous generator contribute to the stability of the commercial power system. On the other hand, it is pointed out that grid-connected inverters usually have no inertia and no synchronizing power, and thus, if the penetration level of renewable energy increases, the ability to stabilize the power system may be insufficient.

As a solution to this problem, a method has been proposed in which inertia and synchronous power are virtually provided to the grid-connected inverter [

Application of the virtual synchronous generator control to a wind diesel system with energy storage system (ESS) has been proposed in [

In [

As a result, the interconnected inverter of each wind power generator is expected to have the inertia and the synchronous power virtually and supplement the lack of inertia and the synchronizing power of the system. A converter with virtual inertia control uses current control (voltage control) to realize the virtual inertia. In the previous researches, PI control is often used after converting to a rotating coordinate system. An optimal current control method was proposed in [

Such equipment typically requires constraint to prevent excess output over the rated capacity in the event of a grid system failure. If such a system demands strict control performance without considering the saturation of the control input, the control performance may be degraded or unstable. To deal with such problems, a model-recovery-anti-windup compensator scheme of the electro hydro servo-system is proposed to suppress the control saturation caused by external load disturbance and modeling uncertainly [

In previous research, we have proposed discrete-time model following control for optimal current control that changes the state feedback gain for each sample time for a virtual synchronous generator model with a inertia and synchronous power [

In this paper, we propose an optimal current control method with suppressed inverter side filter current, the inverter side filter voltage, and the inverter side filter power, and consider the optimized weight of the evaluation function.

This paper describes the proposed control system for the grid-connected wind power generator system in Chapter 2. Chapter 3 describes the simulation result. Finally, we summarize the results in Chapter 4.

We consider a large-scale wind power generator system with a grid-connected inverter as shown in

A wind generator is connected to the grid via AC/DC converter, grid-connected inverter, LCL filter, and grid-connected transformer.

In this paper, we focus on the interconnected inverter and the system, assuming that the wind turbine of the wind generator is rotating under the rated constant wind speed.

An equivalent circuit of the interconnected inverter with an LCL filter connected to the system is shown (see

Since this inverter is connected to a balanced three-phase power system,

The grid system model is composed of equivalent system impedance and infinite bus. The DC link capacitor and the full-bridge inverter are simulated with an ideal voltage source and are connected to the system via an LCL filter. The following differential equation is defined for the system model to which the interconnection inverter is connected, as the control plant model.

where, the filter output current vector i_{po} is composed of zero-phase, negative-phase, and positive-phase components of the three-phase current, and each component is an instantaneous current normalized by the rated current. Similarly, i_{i} is the input current vector of the LCL filter. v_{c} standardizes the instantaneous value vector of the star-connected capacitance voltage normalized by the rated voltage. Similarly, v_{i} is the filter input voltage vector, v_{o} is the output voltage vector, and v_{g} is the infinite bus voltage vector. C_{Y} is the matrix of star-connected capacitances multiplied by the rated impedance, and l_{1} and l_{2} are the input and output inductance matrices of the filter normalized by the rated impedance. T_{1} and T_{2} are the input and output time constant matrices of the filter. The above matrix is defined by the following equation.

where, the rated capacity and rated voltage is the self-capacity of the virtual generator, the rated capacity is S_{base} [VA] and the rated phase voltage is V_{baseY} [V]. The three-phase to two-phase conversion matrix is_{LCL}_{1} [Ω] and R_{LCL}_{2} [Ω], Inverter-side inductance of LCL filter in scalaris L_{LCL}_{1} [H], Grid-side LCL filter inductance in scalaris L_{LCL}_{2} [H], Grid system inductance component in scalaris L_{grid} [H], and the star-connected capacitance in scalar is C_{LCL} [F]. A balanced three-phase circuit is assumed, so the parameters of the three-phase circuit are the same for each phase.

The sum of Equations (1) and (2) gives the following equation:

The voltage vector of the infinite bus is given by the following equation.

where, v is a unit vector of the system voltage. The magnitude of the system voltage |v_{g}|(= C_{pvg}) is a constant scalar quantity. ω_{base} [rad/sec] is a scalar quantity and the rated system frequency, and S_{αβ} is the following oblique matrix.

Assuming that the input voltage vector of filter v_{i} and the grid system voltage vector v_{g} are constant at the time interval T_{s}, the values at the discrete time k are defined as v_{ik} and v_{gk}. Then, Equations (3) to (8) can be represented by the following discrete-time equations.

where i_{pok}, v_{ck}, i_{ik}, and v_{k} are discrete-time vectors at k-th time of the continuous time vector i_{po}, v_{c}, i_{i}, and v at the time interval T_{s}, A_{pio}, A_{pii}, and A_{pvg} are the transition matrices obtained from the differential Equations (6), (3), and (4), respectively, and are given as follows.

By combining Equations (10) to (14), the equation of the controlled model can be expressed as follows.

where, the state vector is composed of _{k}, and the input vector is the filter input voltage vector v_{ik}.

The output voltage v_{o} is given by the following equation from Equations (1) and (2).

Therefore, the output voltage v_{ok} at the time k is obtained by the following equation.

Therefore, from Equations (19) and (14), the output voltage is obtained from the state vectors x_{p}_{1k} and v_{k}.

The swing equation of the virtual generator is given as follows.

Here, M_{g} [s] is a scalar and is the inertia constant of the virtual generator (twice the stored energy constant). T_{g} [s] is a scalar quantity and the time constant of the mechanical system, ω_{g} [rad/sec] is a scalar quantity and the rotation speed of the virtual synchronous generator, τ_{i} [pu] is the input torque as a scalar quantity, τ_{o} [pu] is the output torque as a scalar quantity. These to rquesare, when P_{i} [pu] is input as a scalar quantity and P_{o} [pu] is output as a scalar quantity, _{base} [rad/sec] is equal to the grid rated frequency because the virtual synchronous generator is a 2-pole virtual synchronous generator.

The equivalent circuit equation of the virtual generator is as follows.

where, T_{vo} [s] is the time constant matrix of the electric system of the virtual generator, l_{g} is the inductance matrix, i_{vo} is the output current vector of the virtual generator, and e_{g} is the no-load voltage vector. These vectors are composed of zero-phase, positive-phase, and negative-phase components in the grid-connected inverter system. e_{g} is given by the following equation.

where, e is the unit vector of the no-load voltage vector, and |e_{g}| [pu] is the magnitude of the no-load voltage. |e_{g}| is given by the following equation.

_{avg}, is the output voltage regulator gain, _{o}| [pu] is the magnitude of the output voltage. _{o}| and keeps the output voltage at the target value.

where, K_{gov} is the speed adjustment gain.

The input active power P_{i} of the virtual generator is composed of the output power reference

Assuming that ω_{g}, v_{o}, P_{o}, and P_{i} are constant as ω_{gk}, v_{ok}, P_{ok}, and P_{ik} during the sampling time T_{s}, the discrete-time system of Equation (20) to (25) can be expressed by the following equation.

where τ_{ik} and τ_{ok} are continuous-time scalar variables of τ_{i} and τ_{o} at discrete time at k, and v_{ok}, e_{gk}, i_{vok}, and e_{k} are continuous-time vectors of v_{o}, e_{g}, i_{vo}, and e at time k. A_{vw}, A_{vi}, and A_{ve} are the transition matrices obtained by differential Equations (20), (21), and (22), respectively, and are defined by the following equations.

According to the differential equation solution formula, B_{vw} and B_{vi} are defined as follows.

The control system is a hybrid system composed of three elements. This system is composed of three blocks: a continuous-time full-bridge inverter circuit connected to the grid system via an LCL filter, a discrete-time virtual generator model, and a state feedback system.

The virtual generator model is a time-varying vibration system. The model is not reachable. Therefore, the control system is derived by dividing the state space of the control system into reachable space and the non-reachable space.

The overall system equation is given by

Equation (32), the state equation, is divided into reachable space and others as follows.

Here, the state vector of the reachable space is

where,

(0 in the above equation is a zero matrix, 1 is a unit matrix)

We have the quadratic evaluation function as follows

where, state weight Q, state operation weight S and operation weight R are given by positive definite matrix.

In this case, the optimal input is:

where P is the solution of the following discrete-time Riccati equation:

where, the matrix P is divided into small matrices.

where

Since the lower submatrix of the B matrix in Equation (32) is a zero matrix, Equation (38) becomes as follows.

Therefore, if the solutions of P_{11} and P_{12} satisfying the following equations are obtained, it is enough to obtain Equation (43).

Since A_{11} and B_{1} are reachable, there is P_{11} that satisfies Equation (42).

P_{12} can also be obtained from the following equation.

Since the virtual generator model is a time-varying system, the optimal state feedback gain is derived for each sampling time.

The rated output of the virtual generator model is 1.2 MVA, the rated voltage is 400 V, and the rated frequency of the system is 50 Hz. The synchronous impedance is 2.0 pu, the synchronous impedance resistance component is 7.5 × 10^{−2} pu (0.01Ω), the inertia constant is 3.5 sec, and the mechanical time constant is 100 sec because the friction loss is assumed to be small. The initial internal phase difference angle cosδ = 0.83 and the reference back electromotive force_{gov}: 1/0.05 ω_{base}, voltage adjustment gain K_{avg}: 0.5.

The weight of the evaluation function of the optimal control system is optimized by a multi-objective genetic algorithm, where, the operation weight R = r1, the state weight is_{ii}, and the error weight q_{err} (“1” is unit matrix, “0” is zero matrix).

The weight s_{q} is for the reactive power of the input and the weight for the active power is s_{p}, and the lower submatrix S_{2} is defined as a zero matrix.

(1 in the above equation is the unit matrix)

We define the following objective function to be minimized for searching for Pareto optimal parameters using a multi-objective genetic algorithm.

< >

1) Average root-mean-square of the error current

2) Maximum peak of the filter input voltage

3) Maximum peak of the filter input current

4) Inductance of the inverter side L_{LCL}_{1} [H]

5) Inductance of the grid side L_{LCL}_{2} [H]

6) Capacitance of the filter C_{LCL} [F]

where, the output current error is given by the following equation.

The primary purpose of the control is to minimize the error between the output current and the model output current. The purpose of this evaluation is (1).

On the other hand, a trade-off relationship is expected between the error current and the filter input voltage. Therefore, at the same time, (2) is added as the purpose for minimizing the input voltage.

The purpose (3) is to minimize the output of the inverter because the capacity of the inverter is limited.

In the proposed inverter, a filter is placed between the grid system and the inverter optimal controlled for an arbitrary system. Cost reduction can be expected if the inductance and capacitance of the designed filter are small. For this purpose, objectives (4) to (6) are added.

In consideration of the above, we solve the optimization problem with 8 objective functions.

The following eight design variables are applied to multi-objective genetic optimization as genes.

< >

1) Inductance of the inverter side L_{LCL}_{1} [H]

2) Inductance of the grid side L_{LCL}_{2} [H]

3) Cut-off frequency ω_{cutoff} [rad/sec]

4) Error weight q_{err}

5) Weight of the filter input current q_{ii}

6) Weight of the filter input voltage r

7) Weight of the filter input reactive power s_{q}

8) Weight of the filter input active power s_{p}

Individual with large errors are not design solutions that can be selected, so individual with errors exceeding 10% of the rating are rejected. Individual with a maximum input current exceeding 60 [pu] are discarded. Since the cut-off frequency is based on the rated frequency and the switching frequency is 7.4 kHz, solutions of 74th order or higher are discarded.

Based on [

Step 1: Initialize

Give the individual i_{i} (60 individuals) as element of the external set

Step 2: Fitness assignment

The fitness F(i_{i}) of individual i_{i} is following equation:

where the strength S(i_{i}) of individual I_{i} is the number of solutions it dominates:

where the number of elements of the set A is expressed as num (A) and the equation

Step 3: Selection

Delete one of the closest individuals in the objective function space. Good individuals are selected for the next generation of population set using binary tournament selection.

Step 4: Recombination and mutation operators

Step 5: Algorithm termination judgment

Return to step 2, if not the last generation (60 generations).

O_{j} (I_{i}) is the objective function with objective number j, the ith individual with a set of genes {1..8} is i_{i}.

Here, objective functions (1) to (3) were evaluated by Runge-Kutta method for the initial time-varying of each individual (design solution) in the ideal hybrid system of continuous time system (solution time step: 0.01Ts) and discrete time system with sample time T_{s}.

The markers in

Referring to

The filter input current tends to increase as the error weight increases and the output current error decreases. It seems that individuals are gathering to one place where the input current increases by increasing the error weight and the effect of decreasing the input current by increasing the input current weight is cancelled each other.

On the other hand, the active power weight and the reactive power weight are concentrated in one place, but the distribution is wide (_{q} snd s_{p} is insufficient.

From the figure, it can be seen that the higher the cut-off frequency, the lower the peak current becomes in the Pareto curve. The absence of a design solution of order 11 or higher is considered to be the result of inferior solution discard of large error current.

From the Pareto optimal solution, the filter and load are ω_{base} L_{LCL}_{1}: 2.36 × 10^{−2} pu, ω_{base} L_{LCL}_{2}: 4.41 × 10^{−2} pu, cut-off frequency 10.2 × ω_{base} [rad/sec], error weight q_{err}: 1.44× 10^{13}, input current weight q_{ii}: 1.86 × 10^{9}, input voltage weight r: 14.0, active power weight s_{p}: 5.60 × 10^{8}, reactive power weight s_{q}: 1.01 × 10^{9}.

Since the input of the model changes, the output of the inverter also fluctuates step by step (

the waveform follows the target waveform with an error of 0.5% or less (

The results of the instantaneous voltage drop using the Pareto-optimized load described in the previous section are shown.

In the sequence of the instantaneous voltage drop, it is assumed that the voltage drops by 10% at 4 seconds and returns after 0.1 seconds. The system

parameters are the same as in the previous section.

As seen from

The maximum peak was 3.5pu without weight (_{ii}, s_{p}, and s_{q}, it was affected by the weight, and the maximum was 2.50 pu, but was reduced to 1.85 pu (

It is shown that the inverter rating can be reduced by optimizing the load of the evaluation function while achieving the same performance.

In this paper, as a current control method for a virtual generator model with synchronization power, we proposed the application of discrete-time model following control considering output suppression by optimizing the gain of a multivariable evaluation function.

Since the proposed control method uses instantaneous value control, it can follow the virtual generator model well.

Since a multivariable evaluation function was used, not only the tracking error but also the suppression of the input voltage and input current from the inverter device to the filter could be considered in the control system.

In this control method, good tracking performance and output suppression can be expected by deriving the feedback gain for each sample time even when the frequency fluctuation is large.

Therefore, it is concluded that the proposed method is useful for realizing virtual generator models that contribute to enhancing the grid system stability.

Grid stabilization depends on the virtual synchronous generator model performance. As a future study, we are planning to propose new design method for the condition of output suppression in the virtual generator model supporting and stabilizing the grid system.

Supported by JSPS KAKENHI Grant Number JP17K06289.

The authors declare no conflicts of interest regarding the publication of this paper.

Umemura, A., Takahashi, R. and Tamura, J. (2020) A Study of Discrete-Time Optimum Current Controller for the Virtual Synchronous Generator under Constraints. Smart Grid and Renewable Energy, 11, 1-20. https://doi.org/10.4236/sgre.2020.111001

[

[

[

[

[

[

[

[

[

[

[

V_{inva}, V_{invb}, V_{invc}: the filter input voltage (a-b-c phase) [V]

L_{LCL}_{1}: the inverter-side inductance of the LCL filter [H]

R_{LCL}_{1}: the inverter-side inductance loss of the LCL filter [Ω]

C_{LCL}: the star connected capacitance of the LCL filter [F]

L_{LCL}_{2}: the grid-side inductance of the LCL filter [H]

R_{LCL}_{2}: the grid-side inductance loss of the LCL filter [Ω]

L_{grid}: the grid system inductance component [H]

R_{grid}: the grid system loss [Ω]

V_{grida}, V_{gridb}, V_{gridc}: the grid system voltage (a-b-c phase) [V]

V_{cYa}, V_{cYb}, V_{cYc}: the capacitance voltage (a-b-c phase) [V]

V_{outa}, V_{outb}, V_{outc}: the output voltage of the LCL filter (a-b-c phase) [V]

I_{inva}, I_{invb}, I_{invc}: the input current of the LCL filter (a-b-c phase) [A]

I_{outa}, I_{outb}, I_{outc}: the output current of the LCL filter (a-b-c phase) [A]

S_{base}: the rated capacity of the virtual generator [VA]

V_{baseY}: the rated voltage of the virtual generator [V]

ω_{base}: the rated system frequency [rad/sec]

ω_{cutoff}: the cutoff frequency of the LCL filter [rad/sec]

T_{s}: sampling time [sec]

k: discrete-time

i_{po}: the normalized LCL filter output current vector [pu]

i_{i}: the normalized LCL filter input current vector [pu]

v_{c}: the normalized LCL filter capacitance voltage vector [pu]

v_{o} : the normalized LCL filter output voltage vector [pu]

v_{g} : the normalized grid system voltage vector [pu]

v: the grid system voltage unit vector

v_{i}: the normalized LCL filter input voltage vector [pu]

T_{gd}: the time constant matrix of the grid system [sec]

l_{gd}: the inductance matrix of the grid system

T_{1}: the time constant matrix of the inverter-side LCL filter [sec]

l_{1}: the inductance matrix of the inverter-side LCL filter

T_{2}: the time constant matrix of the grid-side LCL filter [sec]

l_{2}: the inductance matrix of the grid-side LCL filter

c_{Y}: the capacitance matrix of the LCL filter

ω_{g}: the rotation speed of the virtual generator [rad/sec]

P_{i}, P_{o}: the input power and the output power of the virtual generator [pu]

τ_{i}, τ_{o}: the input torque and the output torque of the virtual generator [pu]

T_{g}: the time constant of the virtual generator [sec]

M_{g}: the inertia constant of the virtual generator [sec]

i_{vo}: the output current vector of the virtual generator [pu]

e_{g}: the normalized no-load voltage vector of the virtual generator [pu]

e: the no-load voltage unit vector of the virtual generator [pu]

T_{vo}: the time constant matrix of the virtual generator [sec]

l_{g}: the inductance matrix of the virtual generator

K_{avg}: the voltage regulator gain of the virtual generator

K_{gov}: the speed adjustment gain of the virtual generator

ε_{kiok}: the output current error discrete-time vector

z_{k}: the output current error evaluation state discrete-time vector

q_{ii}: the filter input current weight

q_{err}: the error weight

r: the filter input voltage weight

s_{p}: the filter input active power weight

s_{q}: the filter input reactive power weight

I: the population set

i_{i}: the ith individual (design solution)

F(i_{i}): the fitness function of the ith individual

S(i_{i}): the strength function of the ith individual

O_{j}(i_{i}): the jth objective function of the ith individual

1) Average root-mean-square of the error current

2) Maximum peak of the filter input voltage

3) Maximum peak of the filter input current

4) Inductance of the inverter side L_{LCL}_{1} [H]

5) Inductance of the grid side L_{LCL}_{2} [H]

6) Capacitance of the filter C_{LCL} [F]

where, the output current error is given by the following equation.

The primary purpose of the control is to minimize the error between the output current and the model output current. The purpose of this evaluation is (1).

On the other hand, a trade-off relationship is expected between the error current and the filter input voltage. Therefore, at the same time, (2) is added as the purpose for minimizing the input voltage.

The purpose (3) is to minimize the output of the inverter because the capacity of the inverter is limited.

In the proposed inverter, a filter is placed between the grid system and the inverter optimal controlled for an arbitrary system. Cost reduction can be expected if the inductance and capacitance of the designed filter are small. For this purpose, objectives (4) to (6) are added.

In consideration of the above, we solve the optimization problem with 8 objective functions.

The following eight design variables are applied to multi-objective genetic optimization as genes.

< >

1) Inductance of the inverter side L_{LCL}_{1} [H]

2) Inductance of the grid side L_{LCL}_{2} [H]

3) Cut-off frequency ω_{cutoff} [rad/sec]

4) Error weight q_{err}

5) Weight of the filter input current q_{ii}

6) Weight of the filter input voltage r

7) Weight of the filter input reactive power s_{q}

8) Weight of the filter input active power s_{p}

Individual with large errors are not design solutions that can be selected, so individual with errors exceeding 10% of the rating are rejected. Individual with a maximum input current exceeding 60 [pu] are discarded. Since the cut-off frequency is based on the rated frequency and the switching frequency is 7.4 kHz, solutions of 74th order or higher are discarded.

Based on [

Step 1: Initialize

Give the individual i_{i} (60 individuals) as element of the external set

Step 2: Fitness assignment

The fitness F(i_{i}) of individual i_{i} is following equation:

where the strength S(i_{i}) of individual I_{i} is the number of solutions it dominates:

where the number of elements of the set A is expressed as num (A) and the equation

Step 3: Selection

Delete one of the closest individuals in the objective function space. Good individuals are selected for the next generation of population set using binary tournament selection.

Step 4: Recombination and mutation operators

Step 5: Algorithm termination judgment

Return to step 2, if not the last generation (60 generations).

O_{j} (I_{i}) is the objective function with objective number j, the ith individual with a set of genes {1..8} is i_{i}.

Here, objective functions (1) to (3) were evaluated by Runge-Kutta method for the initial time-varying of each individual (design solution) in the ideal hybrid system of continuous time system (solution time step: 0.01Ts) and discrete time system with sample time T_{s}.

The markers in

Referring to

The filter input current tends to increase as the error weight increases and the output current error decreases. It seems that individuals are gathering to one place where the input current increases by increasing the error weight and the effect of decreasing the input current by increasing the input current weight is cancelled each other.

On the other hand, the active power weight and the reactive power weight are concentrated in one place, but the distribution is wide (_{q} snd s_{p} is insufficient.

From the figure, it can be seen that the higher the cut-off frequency, the lower the peak current becomes in the Pareto curve. The absence of a design solution of order 11 or higher is considered to be the result of inferior solution discard of large error current.

From the Pareto optimal solution, the filter and load are ω_{base} L_{LCL}_{1}: 2.36 × 10^{−2} pu, ω_{base} L_{LCL}_{2}: 4.41 × 10^{−2} pu, cut-off frequency 10.2 × ω_{base} [rad/sec], error weight q_{err}: 1.44× 10^{13}, input current weight q_{ii}: 1.86 × 10^{9}, input voltage weight r: 14.0, active power weight s_{p}: 5.60 × 10^{8}, reactive power weight s_{q}: 1.01 × 10^{9}.

Since the input of the model changes, the output of the inverter also fluctuates step by step (

the waveform follows the target waveform with an error of 0.5% or less (

The results of the instantaneous voltage drop using the Pareto-optimized load described in the previous section are shown.

In the sequence of the instantaneous voltage drop, it is assumed that the voltage drops by 10% at 4 seconds and returns after 0.1 seconds. The system

parameters are the same as in the previous section.

As seen from

The maximum peak was 3.5pu without weight (_{ii}, s_{p}, and s_{q}, it was affected by the weight, and the maximum was 2.50 pu, but was reduced to 1.85 pu (

It is shown that the inverter rating can be reduced by optimizing the load of the evaluation function while achieving the same performance.

In this paper, as a current control method for a virtual generator model with synchronization power, we proposed the application of discrete-time model following control considering output suppression by optimizing the gain of a multivariable evaluation function.

Since the proposed control method uses instantaneous value control, it can follow the virtual generator model well.

Since a multivariable evaluation function was used, not only the tracking error but also the suppression of the input voltage and input current from the inverter device to the filter could be considered in the control system.

In this control method, good tracking performance and output suppression can be expected by deriving the feedback gain for each sample time even when the frequency fluctuation is large.

Therefore, it is concluded that the proposed method is useful for realizing virtual generator models that contribute to enhancing the grid system stability.

Grid stabilization depends on the virtual synchronous generator model performance. As a future study, we are planning to propose new design method for the condition of output suppression in the virtual generator model supporting and stabilizing the grid system.

Supported by JSPS KAKENHI Grant Number JP17K06289.

The authors declare no conflicts of interest regarding the publication of this paper.

Umemura, A., Takahashi, R. and Tamura, J. (2020) A Study of Discrete-Time Optimum Current Controller for the Virtual Synchronous Generator under Constraints. Smart Grid and Renewable Energy, 11, 1-20. https://doi.org/10.4236/sgre.2020.111001

V_{inva}, V_{invb}, V_{invc}: the filter input voltage (a-b-c phase) [V]

L_{LCL}_{1}: the inverter-side inductance of the LCL filter [H]

R_{LCL}_{1}: the inverter-side inductance loss of the LCL filter [Ω]

C_{LCL}: the star connected capacitance of the LCL filter [F]

L_{LCL}_{2}: the grid-side inductance of the LCL filter [H]

R_{LCL}_{2}: the grid-side inductance loss of the LCL filter [Ω]

L_{grid}: the grid system inductance component [H]

R_{grid}: the grid system loss [Ω]

V_{grida}, V_{gridb}, V_{gridc}: the grid system voltage (a-b-c phase) [V]

V_{cYa}, V_{cYb}, V_{cYc}: the capacitance voltage (a-b-c phase) [V]

V_{outa}, V_{outb}, V_{outc}: the output voltage of the LCL filter (a-b-c phase) [V]

I_{inva}, I_{invb}, I_{invc}: the input current of the LCL filter (a-b-c phase) [A]

I_{outa}, I_{outb}, I_{outc}: the output current of the LCL filter (a-b-c phase) [A]

S_{base}: the rated capacity of the virtual generator [VA]

V_{baseY}: the rated voltage of the virtual generator [V]

ω_{base}: the rated system frequency [rad/sec]

ω_{cutoff}: the cutoff frequency of the LCL filter [rad/sec]

T_{s}: sampling time [sec]

k: discrete-time

i_{po}: the normalized LCL filter output current vector [pu]

i_{i}: the normalized LCL filter input current vector [pu]

v_{c}: the normalized LCL filter capacitance voltage vector [pu]

v_{o} : the normalized LCL filter output voltage vector [pu]

v_{g} : the normalized grid system voltage vector [pu]

v: the grid system voltage unit vector

v_{i}: the normalized LCL filter input voltage vector [pu]

T_{gd}: the time constant matrix of the grid system [sec]

l_{gd}: the inductance matrix of the grid system

T_{1}: the time constant matrix of the inverter-side LCL filter [sec]

l_{1}: the inductance matrix of the inverter-side LCL filter

T_{2}: the time constant matrix of the grid-side LCL filter [sec]

l_{2}: the inductance matrix of the grid-side LCL filter

c_{Y}: the capacitance matrix of the LCL filter

ω_{g}: the rotation speed of the virtual generator [rad/sec]

P_{i}, P_{o}: the input power and the output power of the virtual generator [pu]

τ_{i}, τ_{o}: the input torque and the output torque of the virtual generator [pu]

T_{g}: the time constant of the virtual generator [sec]

M_{g}: the inertia constant of the virtual generator [sec]

i_{vo}: the output current vector of the virtual generator [pu]

e_{g}: the normalized no-load voltage vector of the virtual generator [pu]

e: the no-load voltage unit vector of the virtual generator [pu]

T_{vo}: the time constant matrix of the virtual generator [sec]

l_{g}: the inductance matrix of the virtual generator

K_{avg}: the voltage regulator gain of the virtual generator

K_{gov}: the speed adjustment gain of the virtual generator

ε_{kiok}: the output current error discrete-time vector

z_{k}: the output current error evaluation state discrete-time vector

q_{ii}: the filter input current weight

q_{err}: the error weight

r: the filter input voltage weight

s_{p}: the filter input active power weight

s_{q}: the filter input reactive power weight

I: the population set

i_{i}: the ith individual (design solution)

F(i_{i}): the fitness function of the ith individual

S(i_{i}): the strength function of the ith individual

O_{j}(i_{i}): the jth objective function of the ith individual