Grid Interconnection Modeling of Inverter Based Resources (IBR) Plant for Transient Analysis

: The increase in penetration levels of inverter-based resources (IBRs) is changing the dynamic performance of power grids of different parts of the world. IBRs are now being more and more integrated into the grid at a single connection point as an IBR plant. Due to the complex nature and dynamicity of each inverter model, it is not realistic to build and analyze full complex models of each inverter in the IBR plant. Moreover, simulating a large plant including detailed models of all the IBRs would require high computing resources as well as a long simulation time. This has been the main issue addressed in the new IEEE Std 2800-2022. This paper proposes a novel approach to model an IBR plant, which can capture the transient nature at the plant level, detailed IBR control at the inverter level, interactions of multiple IBR groups in a plant structure, and a collector system connecting the IBRs to the grid. The IBRs in the plant use a voltage source inverter topology combined with a grid-connected ﬁlter. The control structure of the IBR includes a cascaded loop control where an inner current control and outer power control are designed in the dq-reference frame, and a closed-loop phase-locked loop is used for the grid synchronization. The mathematical study is conducted ﬁrst to develop aggregated plant models considering different operating scenarios of active IBRs in an IBR plant. Then, an electromagnetic transient simulation (EMT) model of the plant is developed to investigate the plant’s dynamic performance under different operating scenarios. The performance of the aggregated plant model is compared with that of a detailed plant model to prove the effectiveness of the proposed strategy. The results show that the aggregated EMT simulation model provides almost the same result as the detailed model from the plant perspective while the running time/computation burden is much lower.


Introduction
Renewable energy resources (RESs), such as wind turbines and solar photovoltaic (PV) modules, are typically connected to the grid using inverters, and are therefore named inverter-based resources (IBRs) [1].With the fast growth of renewable energy, solar and wind in plant structures are emerging rapidly worldwide, which presents an urgent need for an accurate IBR plant model that can be used to investigate the integration of IBR plants with the grid and which drives a new IEEE Std 2800-2022 published on 22 April 2022 [2].An IBR plant consists of a huge number of IBR units connected together via a collector system at a single location [3], and, after that, the aggregated power is transported through a transmission line and sent to the grid.An IBR operation is different and complicated compared to that of a synchronous generator as it has more complicated dynamics.When hundreds of IBRs are tied together and operate as a plant, the transient complexity increases enormously, and its impact on the grid can vary based on different IBR configurations in the plant [4].
In order to analyze the impact of an IBR plant on the bulk power system (BPS), it is important to develop a transient simulation model that is computationally affordable and can accurately portray the control of the IBRs at the inverter level, as well as the interconnection dynamics of the plant with the grid.Most developments of the IBR plant model available in the literature are related to wind power plants (WPP) [5][6][7].In [8], an equivalent impedance model of a WPP collector system is proposed for power system planning studies, where the wind turbines are set up in series and parallel arrangements.However, detailed transient modeling of IBRs within the plant is absent in [8], and no EMT simulation is performed.Another WPP equivalent model is proposed in [9], which reduces the plant model containing a fixed number of 168 wind turbines into an equivalent model with one or two machines connected to the grid via one or two lines representing the equivalent collector system.In [10], the impact of transient output characteristics of a WPP equivalent model is analyzed and compared with an actual WPP, in which it is found that the model does not provide an accurate response when the output power is larger than the grid short circuit capacity.A PV plant equivalent model is proposed in [11], which uses a model with an aggregated PV array, dc/dc converter, dc/ac inverter, gridconnected filter, and converter controllers.A closed-loop parameter identification approach is proposed in [11] to tune the parameters of the aggregated PV plant and controller models to minimize the error between the measured and model output powers of the plant.Another equivalent model of a PV power plant is presented in [12], which represents 100 PV inverters.The focus of the model is only to mimic the dynamics of the active and reactive power of the actual plant at the POI (point of interconnection).In [13,14], model reduction approaches of multiple IBRs are proposed, with the primary focus on developing an aggregate droop control equivalence for a group of droop-controlled inverters.In [15], a coherency equivalence method is proposed to develop an aggregate model for a group of modular multilevel converters (MCCs) used in large-scale renewable energy systems, with focus mainly on developing an aggregate equivalent model for the power circuit portions of the MMCs.Purba et.al., proposed a reduced-order IBR plant model in [16,17], in which the collector system impact is not considered and IBRs are assumed to be connected in parallel and then tied to the grid.A scaling factor is used in [16] to determine the aggregated plant model parameters, including aggregated grid-connected filter and controller parameters of the plant.In all aggregate plant models developed above, the characteristics at the individual IBR level, including its controller, are unable to simulate and evaluate and the EMT impacts for grid integration of the IBR plant.They are unable to be correctly captured because the differences and variations of collector systems and transmission lines are not properly and fully addressed.A comparison of properties of the existing models and the proposed model is provided in Table 1.This article develops an aggregated IBR plant model that can simulate the dynamic performance at the IBR level as well as the grid interconnection level of the plant, and is fundamentally different from the aggregated plant models developed previously by others.The significant and original contributions of this article comprise: (1) a detailed mathematical proving process to prove the theoretical foundation of the proposed aggregated IBR plant model development under different IBR plant configurations, (2) a novel aggregated IBR plant model that can simulate transient performance at the IBR level, evaluate the detailed and complete design of the IBR controller and grid-connected filter, and investigate the interactions of IBR groups in a plant structure, (3) a strategy that carefully addresses the collector system impact in diverse parallel and series connection scenarios, (4) an aggregate model that can be used to evaluate the transient performance for the grid interconnection of the plant based on the designed IBR system and controller parameters of individual IBRs, instead of scaled and modified parameters used in the aggregate model, and (5) a process to evaluate the effectiveness of the proposed aggregate IBR plant model through a comparison study with the detailed IBR plant model in a recursive way.Details of the article are presented in the following sections.

IBR Plant
The IBR family typically includes solar PV generators, wind turbine (WT) generators such as doubly fed induction generators (Type-3), permanent magnet synchronous generators (Type-4), and energy storage systems [26].Except for the Type-3 WT (not the focus of this paper), IBRs are connected to the BPS through a full-size dc/ac converter system (Figure 1a), which includes four main components: (i) a converter at the source-side, (ii) a capacitor at the dc-link, (iii) a converter at the grid-side, and (iv) a grid connected filter [27][28][29].The converter at the source-side controls the energy supply of solar or wind, or performs energy management of the energy storage system, and is an ac/dc type for a Type-4 WT and a dc/dc type for a solar PV and energy storage system.The grid-side converter regulates reactive power flow at the point of common coupling (PCC) and maintains a fixed dc-link voltage and is always a dc/ac type.An IBR plant usually consists of a collection of closely located IBRs that are electrically linked together to generate electricity.With the declining cost of solar panels and wind turbines, the number of IBR plants are increasing quickly.Utility-scale wind power plants (WPPs) and solar power plants (SPPs) have been built with hundred or thousand MW scale worldwide.Figure 1b shows a generalized connection diagram of a grid connected onshore WPP, where k represents the parallel feeder index, and i represents the index of a series path in each parallel feeder.Typically, the output end voltage of each WT generator is 690 V, and all the WTs within a WPP are connected to a medium-voltage (MV) substation (e.g., 34.5k V/120 kV) via a collector system.The connection point is popularly known as the POM (point of measurement).The POM is then connected to the bulk power system (BPS) via a transmission line, which can be up to several hundred miles long.This connection point is called the POI (point of interconnection).Similarly, a utility-scale SPP is built with hundreds to thousands of PV panel strings.A string-inverter connects each string with the collector system.The collector system is connected to the POM via a step-up transformer.Finally, a transmission line connects the SPP with the POI via an interconnection transformer.

Control of an IBR Plant
An IBR operating in an IBR plant can be assumed to have a constant dc-link voltage with a controller for the grid-side inverter, which consists of a back-to-back inner current control loop and an outer power (active/reactive) or bus voltage control loop.The control loops are usually designed in the d-q reference frame.The power control reference signals (P * PCC and Q * PCC) are generated in accordance with the grid requirements (either at the POM or POI) and sent to the outer control loop.The outer control loop produces d-and qaxis current reference signals and sends them to the inner control loop.The inner control loop produces d-q controlled output voltages, denoted as v * d_inv and v * q_inv.The controlled voltages are injected to the power grid.The relation between the inverter end voltage (vdq_inv) and the controller output voltage (v * dq_inv) can be written as follows [30]: An IBR plant usually consists of a collection of closely located IBRs that are electrically linked together to generate electricity.With the declining cost of solar panels and wind turbines, the number of IBR plants are increasing quickly.Utility-scale wind power plants (WPPs) and solar power plants (SPPs) have been built with hundred or thousand MW scale worldwide.Figure 1b shows a generalized connection diagram of a grid connected onshore WPP, where k represents the parallel feeder index, and i represents the index of a series path in each parallel feeder.Typically, the output end voltage of each WT generator is 690 V, and all the WTs within a WPP are connected to a medium-voltage (MV) substation (e.g., 34.5k V/120 kV) via a collector system.The connection point is popularly known as the POM (point of measurement).The POM is then connected to the bulk power system (BPS) via a transmission line, which can be up to several hundred miles long.This connection point is called the POI (point of interconnection).
Similarly, a utility-scale SPP is built with hundreds to thousands of PV panel strings.A string-inverter connects each string with the collector system.The collector system is connected to the POM via a step-up transformer.Finally, a transmission line connects the SPP with the POI via an interconnection transformer.

Control of an IBR Plant
An IBR operating in an IBR plant can be assumed to have a constant dc-link voltage with a controller for the grid-side inverter, which consists of a back-to-back inner current control loop and an outer power (active/reactive) or bus voltage control loop.The control loops are usually designed in the d-q reference frame.The power control reference signals (P * PCC and Q * PCC ) are generated in accordance with the grid requirements (either at the POM or POI) and sent to the outer control loop.The outer control loop produces dand q-axis current reference signals and sends them to the inner control loop.The inner control loop produces d-q controlled output voltages, denoted as v * d_inv and v * q_inv .The controlled voltages are injected to the power grid.The relation between the inverter end voltage (v dq_inv ) and the controller output voltage (v * dq_inv ) can be written as follows [30]: where k PWM is the multiplication constant between inverter end voltage and the controller output voltage caused by the PWM (pulse width modulation) process [30].Figure 2 illustrates the overall plant control structure, as well as its interaction with individual IBR level controller.
Energies 2022, 15, x FOR PEER REVIEW 5 of 20 where kPWM is the multiplication constant between inverter end voltage and the controller output voltage caused by the PWM (pulse width modulation) process [30].Figure 2 illustrates the overall plant control structure, as well as its interaction with individual IBR level controller.

Approximate EMT Circuit Models of IBR Plant Connected to the Grid
An IBR is usually interfaced to the collector system through a passive filter.The filter could be L, LC, or LCL type.Figure 3 shows the circuit diagram of an IBR connected to the grid via an LCL-filter.In the figure, RfI and LfI represent the inverter-side filter resistance and inductance, RfG and LfG denote the grid-side filter resistance and inductance, C stands for the filter capacitance, and Rg and Lg represent the equivalent grid resistance and inductance from the PCC.In Figure 3, vα,β,γ_inv, vα,β,γ_PCC, iα,β,γ_inv, iα,β,γ, vc_α,β,γ, and vα,β,γ  The equation of voltage balance across the inductor of the inverter-side can be written as, The current equation of the LCL capacitor is,

Approximate EMT Circuit Models of IBR Plant Connected to the Grid
An IBR is usually interfaced to the collector system through a passive filter.The filter could be L, LC, or LCL type.Figure 3 shows the circuit diagram of an IBR connected to the grid via an LCL-filter.In the figure, R fI and L fI represent the inverter-side filter resistance and inductance, R fG and L fG denote the grid-side filter resistance and inductance, C stands for the filter capacitance, and R g and L g represent the equivalent grid resistance and inductance from the PCC.In Figure 3, v α,β,γ_inv , v α,β,γ_PCC , i α,β,γ_inv , i α,β,γ , v c_α,β,γ , and  v α,β,γ represent the inverter terminal voltage, PCC voltage, current through the R fI and L fI , current injected to the grid, LCL filter capacitor voltage, and grid voltage, respectively.where kPWM is the multiplication constant between inverter end voltage and the controller output voltage caused by the PWM (pulse width modulation) process [30].Figure 2 illustrates the overall plant control structure, as well as its interaction with individual IBR level controller.

Approximate EMT Circuit Models of IBR Plant Connected to the Grid
An IBR is usually interfaced to the collector system through a passive filter.The filter could be L, LC, or LCL type.Figure 3 shows the circuit diagram of an IBR connected to the grid via an LCL-filter.In the figure, RfI and LfI represent the inverter-side filter resistance and inductance, RfG and LfG denote the grid-side filter resistance and inductance, C stands for the filter capacitance, and Rg and Lg represent the equivalent grid resistance and inductance from the PCC.In Figure 3, vα,β,γ_inv, vα,β,γ_PCC, iα,β,γ_inv, iα,β,γ, vc_α,β,γ, and vα,β,γ  The equation of voltage balance across the inductor of the inverter-side can be written as, The current equation of the LCL capacitor is, The equation of voltage balance across the inductor of the inverter-side can be written as, The current equation of the LCL capacitor is, The equation of voltage balance across the inductor of the grid-side can be written as, The equation of voltage balance across the impedance of the grid can be written as Figure 4 shows the equivalent circuit of an IBR plant.In the figure, IBRs are coupled to the power grid through a collector system.R CSi , and L CSi correspond to the collector system resistance and inductance.To analyze the plant structure and keep the collector system into consideration, four cases are developed in this paper, which are described as follows: 1.
The collector system impedance is neglected and all IBRs in the IBR plant have the same filter parameter. a.
All IBRs in the IBR plant have the same output terminal voltage.b.
All IBRs in the IBR plant have different output terminal voltages. 2.
The collector system impedance is considered, and all IBRs in the IBR plant have distinct filter characteristics. a.
Each IBR in the IBR plant is coupled to the POM in parallel through different collector impedances; therefore, they output different terminal voltage.b.
Each IBR in the IBR plant is connected in a series and parallel combination to the POM through a collector system (depicted in Figure 1b), and they output different terminal voltage.
Energies 2022, 15, x FOR PEER REVIEW 6 of 20 Figure 4 shows the equivalent circuit of an IBR plant.In the figure, IBRs are coupled to the power grid through a collector system.RCSi, and LCSi correspond to the collector system resistance and inductance.To analyze the plant structure and keep the collector system into consideration, four cases are developed in this paper, which are described as follows: 1.The collector system impedance is neglected and all IBRs in the IBR plant have the same filter parameter.
a.All IBRs in the IBR plant have the same output terminal voltage.b.All IBRs in the IBR plant have different output terminal voltages.
2. The collector system impedance is considered, and all IBRs in the IBR plant have distinct filter characteristics.
a.Each IBR in the IBR plant is coupled to the POM in parallel through different collector impedances; therefore, they output different terminal voltage.b.Each IBR in the IBR plant is connected in a series and parallel combination to the POM through a collector system (depicted in Figure 1b), and they output different terminal voltage.
A grid-connected IBR plant including IBRs, IBR filters, collector impedances, and grid impedance.

A Grid-Connected IBR Plant including IBRs, IBR Filters, Collector Impedances, and Grid Impedance
To identify IBRs in the IBR plant, superscript k is used, where k = 1, 2,..., n. n defines the total number of online IBRs in the IBR plant.

Theorem 1. If all IBRs online in the IBR plant have identical filter characteristics and terminal voltage, then each IBR is assumed to be linked to the BPS via a changeable grid impedance, which is equal to the grid impedance (when one IBR is connected) multiplied by the total number of online IBRs (n) in the IBR plant.
Proof.In this case, the collector system impedance is ignored, therefore, the PCC voltage of each IBR would be the POM voltage (according to Figure 4).As each IBR outputs the same voltage and has identical filter parameters, each IBR will have a similar voltage and current, as shown in Equations ( 2)-( 4), and will deliver the same current to the grid.Then,  Theorem 1.If all IBRs online in the IBR plant have identical filter characteristics and terminal voltage, then each IBR is assumed to be linked to the BPS via a changeable grid impedance, which is equal to the grid impedance (when one IBR is connected) multiplied by the total number of online IBRs (n) in the IBR plant.
Proof.In this case, the collector system impedance is ignored, therefore, the PCC voltage of each IBR would be the POM voltage (according to Figure 4).As each IBR outputs the same voltage and has identical filter parameters, each IBR will have a similar voltage and current, as shown in Equations ( 2)-( 4), and will deliver the same current to the grid.Then, under this case, we have As a result, the equation of voltage balance across the impedance of the grid is The equivalent circuit of an IBR plant, looking from an IBR point of view, is developed from Equation ( 7) and shown in Figure 5.In the figure, the IBR is coupled to the BPS through a virtually changeable transmission line impedance.The impedance of the transmission line can be calculated by multiplying the grid impedance by the number of online IBRs in the IBR plant.
Energies 2022, 15, x FOR PEER REVIEW 7 of 20 The equivalent circuit of an IBR plant, looking from an IBR point of view, is developed from Equation ( 7 Lemma 1.If all IBRs in an IBR plant have identical filter characteristics and terminal voltage, then, from the grid point of view, the IBR plant will be assumed to be interfaced to the BPS through a virtually changeable filter impedance, which can be calculated by dividing the individual IBR filter impedance by the total number of online IBRs in the IBR plant. Proof.Define the current flowing through the transmission line of the grid as iα,β,γ_g.Then, (7) can be rewritten as, Similarly ( 6), ( 2), (3), and (4) can be rewritten by using iα,β,γ_g as An IBR plant (consists of n IBRs) equivalent circuit, looking from the grid point of view, is developed from 8 to 11 and shown in Figure 6.The figure shows that the IBR plant is tied to the BPS through a virtually changeable filter impedance.The filter impedance is determined by the parallel combination of filter impedances of all IBRs online in the IBR plant.The equivalent circuit outputs the same voltage as that of individual IBR and injects the sum of currents of IBRs online in the plant.□ Lemma 1.If all IBRs in an IBR plant have identical filter characteristics and terminal voltage, then, from the grid point of view, the IBR plant will be assumed to be interfaced to the BPS through a virtually changeable filter impedance, which can be calculated by dividing the individual IBR filter impedance by the total number of online IBRs in the IBR plant.
Proof .Define the current flowing through the transmission line of the grid as i α,β,γ_g .Then, (7) can be rewritten as, Similarly ( 2)-( 4), ( 6) can be rewritten by using i α,β,γ_g as An IBR plant (consists of n IBRs) equivalent circuit, looking from the grid point of view, is developed from 8 to 11 and shown in Figure 6.The figure shows that the IBR plant is tied to the BPS through a virtually changeable filter impedance.The filter impedance is determined by the parallel combination of filter impedances of all IBRs online in the IBR plant.The equivalent circuit outputs the same voltage as that of individual IBR and injects the sum of currents of IBRs online in the plant.

Case 2: IBRs with Identical Filter Parameters but Different Terminal Voltages
Assume all IBRs have different output terminal voltages.As the collector system impedance is neglected, PCC voltages of all IBRs are the same and equivalent to the POM voltage.

Theorem 2. If all IBRs online in an IBR plant have identical filter parameters but output distinct terminal voltages, then, looking from the grid side, the IBR plant is assumed to be coupled to the BPS by a virtual changeable filter impedance, which can be calculated by dividing the individual IBR filter impedance by the total number of online IBRs in the IBR plant. The IBR plant outputs the average terminal voltage of all IBRs online in that plant.
Proof.As all IBRs have identical filter parameters but output distinct terminal voltages, Equations ( 2)-( 5) can be rewritten as follows: , , , , , , , , , , , , , , For all IBRs, we obtain, Rewrite ( 16), ( 17), (18), and (15) as  Theorem 2 .If all IBRs online in an IBR plant have identical filter parameters but output distinct terminal voltages, then, looking from the grid side, the IBR plant is assumed to be coupled to the BPS by a virtual changeable filter impedance, which can be calculated by dividing the individual IBR filter impedance by the total number of online IBRs in the IBR plant.The IBR plant outputs the average terminal voltage of all IBRs online in that plant.
Proof.As all IBRs have identical filter parameters but output distinct terminal voltages, Equations ( 2)-( 5) can be rewritten as follows: For all IBRs, we obtain, Rewrite ( 15)- (18) as Therefore, an IBR plant equivalent circuit (consist of n IBRs) can be drawn by Equations ( 19)-( 22) and shown in Figure 6.In the figure, an IBR plant is coupled to the grid via a virtually changeable filter impedance, which is equivalent to the parallel combination of filter impedances of all IBRs online in the IBR plant.The equivalent circuit outputs the average terminal voltages and injects the sum of currents of all IBRs online in the IBR plant.
Lemma 2. If all IBRs online in an IBR plant have identical filter characteristics but distinct terminal voltages, then, from the individual IBR point of view, the IBR plant can be assumed to be coupled to the BPS through a virtually changeable grid impedance, which is equivalent to the product of grid impedance and the number of IBRs online in the IBR plant.The IBR plant provides the average terminal voltage of all IBRs in that plant.
Proof.By rewriting ( 15)-( 18) as: An IBR plant equivalent circuit can be drawn from Equations ( 23)-( 26) and looks the same as in Figure 5.In the figure, an IBR plant can be assumed to be coupled to the BPS through a changeable transmission line impedance, which is equivalent to the product of grid impedance and the number of online IBRs in the IBR plant.The equivalent circuit outputs the average terminal voltages and injects the sum of currents of all IBRs online in the IBR plant.To develop an EMT model for an IBR plant, either the time domain (used in Cases 1 and 2) or the Laplace transform technique could be considered.Both methods provide very complicative state-space system equations and IBR plant models, which are not easy to simulate or evaluate.Instead of developing a detailed EMT model, an approximate EMT model of an IBR plant is proposed in this section based on Lemma 2 (developed for Case 2) as follows: (1) An average IBR represents an IBR plant consist of n IBRs, which outputs the average terminal voltage of all IBRs in the IBR plant as follows: (2) The filter characteristics of the average IBR are given by averaging the values of filter parameters of all the IBRs in the IBR plant as follows: (3) The collector system parameters(R CS_avg and L CS_avg ) of the average IBR are given by averaging the collector system parameters (R k CS and L k CS ) connecting all the IBRs to the POM as follows The approximate equivalent EMT circuits of an IBR plant are obtained from Equations ( 27)-( 29) and Figures 6 and 7, and they are shown in Figure 8a,b.

2, 15, x FOR PEER REVIEW
10 of 20 (3) The collector system parameters ( ) connecting all the IBRs to the POM as follows The approximate equivalent EMT circuits of an IBR plant are obtained from Equations ( 27)-( 29) and Figures 6 and 7   Energies 2022, 15, x FOR PEER REVIEW 10 ) connecting all the IBRs to POM as follows The approximate equivalent EMT circuits of an IBR plant are obtained from Eq tions ( 27)-( 29) and Figures 6 and 7    An approximate EMT IBR plant model is proposed in this section based on theorems and lemmas developed for Cases 1 to 3 as follows: (1) An average IBR represents all IBRs in an IBR plant.The average IBR outputs the average terminal voltages of all IBRs in the IBR plant as shown in ( 27). ( 2) The filter characteristics of the average IBR are given by averaging the values of filter parameters of all the IBRs in the IBR plant as shown in (28).(3) Similar to Figure 1b, assume an IBR plant has K parallel feeders, and each parallel feeder has N k number of series paths.Equation (30a) provides the approximate and aggregated average impedance of collector system of each feeder.Therefore, the overall collector system (with K parallel feeders) parameters of an IBR plant are obtained from Equation (30b).
where, R k i and L k i denote the resistance and inductance of the ith series path within the kth parallel feeder, and R k_avg and L k_avg stand for the approximate and aggregated average resistance and inductance of the kth parallel feeder.(4) By replacing the parameters developed in ( 27), (28), and (30) into Figure 7b, the equivalent and approximate EMT circuit model is obtained.

EMT Simulation Model of IBR Plant Grid Interconnection
According to the theoretical study presented in Section 3, the aggregated EMT simulation models of an IBR plant are developed in this section.The development of the aggregated EMT simulation models is presented as correspondence to Cases 1 to 4, from simple to complex/complete as follows.

Approximated EMT Model of Case 1
The aggregated IBR plant model consists of a power circuit model and a controller model.The power circuit model is obtained from Figure 4, in which the POM voltage is applied directly to each individual IBR, while the current flowing into the grid at the POM is the individual IBR current multiplied by n, as shown in Figure 7.For this equivalent plant model, the original controller model of individual IBR can be applied directly without the need for an equivalent plant controller model as developed previously and reported in the literature [27, 30,31].The IBR controller includes an inner current loop controller and an outer power loop controller (Figure 2).For the conditions of Case 1, each IBR outputs the same voltage to the grid at the inverter terminal, which can also be translated to the power references to the controller, which are the same for all the IBRs.Thus, a complete aggregated IBR plant model, including the power circuit and controller, is obtained in Figure 8, where we have the voltage and current relations at the PCC of the IBR and the POM of the grid as follows

Approximated EMT Model of Case 2
The aggregated IBR plant model in this case also consists of a power circuit model and a controller model.In this case, each IBR within an IBR plant has identical filter characteristics.Thus, the power circuit model should be the same as that shown in Figure 8, except that the IBR is represented by an average IBR.According to Lemma 2, the average IBR in Case 2 is considered as being connected to the grid and as having the average of the output voltages and currents of all the IBRs connected online.This means that the power references of IBR controllers are different for all the IBRs.Therefore, the reference power to the average IBR at any time instant is the average of the power references of all the IBRs.In this way, a complete aggregated IBR plant can also be modeled using Figure 7, in which, the power reference of the average IBR at any time instant is the average of the power references of all the IBRs, i.e., where, P * PCC_i (t) and Q * PCC_i (t) are the active and reactive power references of each individual IBR (i = 1, . . ., n).

Approximated EMT Model of Case 3
Case 3 is similar to Case 2. However, there are two major differences between the two cases.First, in Case 3, the parameters of the IBR grid-connected filters are different.Therefore, average values of all the grid-connected filters are used in the average IBR model.In addition, in Case 3, each IBR is connected to the POM through a different collector system impedance.In the aggregated model, this is represented as an average collector resistance R CS_avg and an average collector inductance L CS_avg that are calculated via (29).By considering these two factors, a complete aggregated IBR plant model is developed, as shown in Figure 9, in which the POM' represents the POM referred to the average IBR model side, and the current and voltage relation between the POM and POM' is illustrated by

Approximated EMT Model of Case 3
Case 3 is similar to Case 2. However, there are two major differences between the two cases.First, in Case 3, the parameters of the IBR grid-connected filters are different.Therefore, average values of all the grid-connected filters are used in the average IBR model.In addition, in Case 3, each IBR is connected to the POM through a different collector system impedance.In the aggregated model, this is represented as an average collector resistance _ CS avg R and an average collector inductance L that are calculated via (29).By considering these two factors, a complete aggregated IBR plant model is developed, as shown in Figure 9, in which the POM` represents the POM referred to the average IBR model side, and the current and voltage relation between the POM and POM' is illustrated by In the aggregated IBR plant model, the power reference of the average IBR at any time instant is the average of the power references of all the IBRs, as shown in (32).

Approximated EMT Model of Case 4
The aggregated IBR plant model of Case 4 is the same as the aggregated IBR plant model of Case 3, except that the average collector system resistance and inductance are calculated by using (30) instead of (29).Thus, Figure 9 can be used to represent the aggregated IBR plant model of Case 4, in which the average collector system resistance and inductance are calculated by using (30).The aggregated model of Case 4 can also be applied to Cases 1 to 3 by considering each as a special case of Case 4.

Detailed IBR Plant Model in a Recursive Way
A detailed IBR plant EMT model is important for the validation of the aggregated In the aggregated IBR plant model, the power reference of the average IBR at any time instant is the average of the power references of all the IBRs, as shown in (32).

Approximated EMT Model of Case 4
The aggregated IBR plant model of Case 4 is the same as the aggregated IBR plant model of Case 3, except that the average collector system resistance and inductance are calculated by using (30) instead of (29).Thus, Figure 9 can be used to represent the aggregated IBR plant model of Case 4, in which the average collector system resistance and inductance are calculated by using (30).The aggregated model of Case 4 can also be applied to Cases 1 to 3 by considering each as a special case of Case 4.

Detailed IBR Plant Model in a Recursive Way
A detailed IBR plant EMT model is important for the validation of the aggregated IBR plant model.However, a detailed model with hundreds of IBR units within an IBR plant is almost impossible to build and simulate.To overcome the challenge, a recursive way to build the detailed IBR plant model is developed.The recursive way first started with a detailed IBR plant model containing a certain number of IBR units (for example 10 units), that was then compared with an aggregated IBR plant model with the same filter and collector system configuration of the 10 IBR units.After the validation, the number of IBRs in the IBR plant was increased to 100 units that were divided into 10 groups, including a collector system connecting the 10 groups, with each group modeled by an aggregated IBR model containing 10 IBR units.In this way, a larger IBR plant system with thousands of IBRs can be achieved if needed.This proposed strategy can also be used to evaluate different transient interactions of IBR groups in a plant structure.

Approximated IBR Plant Grid Interconnection Model Simulation and Its Validation with Detailed Recursive Model
To evaluate the equivalent and aggregated IBR plant models developed in Section 4, an EMT simulation model of an IBR plant is built in MATLAB (Figure 10a).The aggregated IBR plant model consists of an average IBR model.The average IBR is coupled to a step-up transformer (690 V/25 kV) through an LCL filter.The LCL filter can be conveniently changed to either an L or LC filter.After the step-up transformer, the average IBR model is coupled to the POM through an equivalent aggregated collector system (Section 4), a 25 kV/120 kV step-up transformer, and then to the POI via a transmission line.At the POI, the grid voltage is 120 kV.The IBR unit is then connected to a renewable energy source (RES), such as a photovoltaic array or a wind turbine, on the left side, and a grid-connected inverter on the right side.According to an active power reference (P*), the RES controller controls the active power extracted from the RES and delivers it to the BPS.The grid-connected inverter controller maintains the IBR interconnection to the grid.The configuration parameter of the aggregated IBR plant model is given in Table 2.

EER REVIEW
13 of 20

Approximated IBR Plant Grid Interconnection Model Simulation and Its Validation with Detailed Recursive Model
To evaluate the equivalent and aggregated IBR plant models developed in Section 4, an EMT simulation model of an IBR plant is built in MATLAB (Figure 10a).The aggregated IBR plant model consists of an average IBR model.The average IBR is coupled to a step-up transformer (690 V/25 kV) through an LCL filter.The LCL filter can be conveniently changed to either an L or LC filter.After the step-up transformer, the average IBR model is coupled to the POM through an equivalent aggregated collector system (Section 4), a 25 kV/120 kV step-up transformer, and then to the POI via a transmission line.At the POI, the grid voltage is 120 kV.The IBR unit is then connected to a renewable energy source (RES), such as a photovoltaic array or a wind turbine, on the left side, and a gridconnected inverter on the right side.According to an active power reference (P*), the RES controller controls the active power extracted from the RES and delivers it to the BPS.The grid-connected inverter controller maintains the IBR interconnection to the grid.The configuration parameter of the aggregated IBR plant model is given in Table 2. Figure 10b shows the inverter controller configuration of the average IBR model, which has a current controller in the inner loop and a dc-link voltage controller and reactive power or PCC voltage controller in the outer loop.The controllers are designed according to the well-known standard vector control strategy [27,28].The PI controller gains  Figure 10b shows the inverter controller configuration of the average IBR model, which has a current controller in the inner loop and a dc-link voltage controller and reactive power or PCC voltage controller in the outer loop.The controllers are designed according to the well-known standard vector control strategy [27,28].The PI controller gains are determined based on the transfer functions and tuned by following the standard bode plot method [32].In addition, a limiting block is added in the inner current control loop to prevent overcurrent [29,33].Similarly, a voltage limiting block is connected to limit the voltage produced by the inner-loop current controller according to [34].This will ensure the IBR operation within the linear modulation boundary.A saturation block is added to the output of all PI controllers to limit the integral control part without crossing the maximum IBR rated current, IBR output voltage, and reactive power limits.Thus, Figure 10 shows the state-of-the-art IBR control technology used in the industry.The aggregated model is used to simulate the IBR plant transient performance for Cases 1 to 4, and it is validated via the comparison with the detailed recursive plant model as follows:

Approximated IBR Plant Model Evaluation of Case 1
For Case 1, the impedance of the collector system is not considered and the parameters for all the IBRs in the plant are the same.The total number of IBRs is 600.Regarding the LCL filter parameters of an IBR, the capacitance is 25 µF while the inductance is 0.2 mH, and the resistance of the filter inductor is 0.0015 Ω for both the inverter-and grid-side inductors.Moreover, the IBRs are synchronized and have same PCC voltage.Thus, the aggregated and detailed IBR plant models should generate the same results, according to the studies presented in Sections 3 and 4. Figure 11 shows a comparison of the EMT simulation using the two models, in which the active power reference changes from 0 kW to 500 kW at 0.4 s, and the reactive power reference changes from 20 kVAR to 40 kVAR at 0.5 s.Table 3 gives a comparison of the current and voltage total harmonic distortions (THDs) by using the aggregated and detailed IBR plant models.Both show that the aggregated and detailed models generate the same results, validating the theoretical study given in Sections 3 and 4.

Approximated IBR Plant Model Evaluation of Case 2
All IBRs within the IBR plant have the same LCL filter parameters in this case, but they are not synchronized and deliver different output voltages.Hence, the power reference of each IBR is different in the detailed model.The power reference of the average IBR in the aggregated model is the average of the power references of all the IBRs, according to (32). Figure 12 shows the active and reactive power at the POM of the IBR plant.The THD comparison by using the aggregated and detailed IBR plant models for Case 2 is given in Table 3. Unlike Case 1, the results obtained by using the aggregated and detailed models are not the same, which is consistent with the theoretical study shown Section 3.However, the difference between the aggregated and detailed models is very small, demonstrating the effectiveness of the proposed aggregated model.

Approximated IBR Plant Model Evaluation of Case 3
In this case, all the IBRs within the IBR plant have different LCL filter parameters, deliver different IBR output voltages, and are connected in parallel with different collector impedances to the POM.For the detail model, a perturbation of 20% on filter inductance and resistance is applied to obtain a random value of the filter parameters.A similar process is applied to the collector impedance.This represents the parallel connected IBRs having different collector lengths and having different filters.For the aggregated model, according to Section 4.3, the filter parameters of the average IBR are the average of all IBR filter parameters calculated using (28), and the collector resistance and inductance are the average of the collector resistance and inductance of all the IBRs calculated using (28).All IBRs within the IBR plant have the same LCL filter parameters in this case, but they are not synchronized and deliver different output voltages.Hence, the power reference of each IBR is different in the detailed model.The power reference of the average IBR in the aggregated model is the average of the power references of all the IBRs, according to (32). Figure 12 shows the active and reactive power at the POM of the IBR plant.The THD comparison by using the aggregated and detailed IBR plant models for Case 2 is given in Table 3. Unlike Case 1, the results obtained by using the aggregated and detailed models are not the same, which is consistent with the theoretical study shown Section 3.However, the difference between the aggregated and detailed models is very small, demonstrating the effectiveness of the proposed aggregated model.

Approximated IBR Plant Model Evaluation of Case 2
All IBRs within the IBR plant have the same LCL filter parameters in this case, but they are not synchronized and deliver different output voltages.Hence, the power reference of each IBR is different in the detailed model.The power reference of the average IBR in the aggregated model is the average of the power references of all the IBRs, according to (32). Figure 12 shows the active and reactive power at the POM of the IBR plant.The THD comparison by using the aggregated and detailed IBR plant models for Case 2 is given in Table 3. Unlike Case 1, the results obtained by using the aggregated and detailed models are not the same, which is consistent with the theoretical study shown Section 3.However, the difference between the aggregated and detailed models is very small, demonstrating the effectiveness of the proposed aggregated model.

Approximated IBR Plant Model Evaluation of Case 3
In this case, all the IBRs within the IBR plant have different LCL filter parameters, deliver different IBR output voltages, and are connected in parallel with different collector impedances to the POM.For the detail model, a perturbation of 20% on filter inductance and resistance is applied to obtain a random value of the filter parameters.A similar process is applied to the collector impedance.This represents the parallel connected IBRs having different collector lengths and having different filters.For the aggregated model, according to Section 4.3, the filter parameters of the average IBR are the average of all IBR filter parameters calculated using (28), and the collector resistance and inductance are the average of the collector resistance and inductance of all the IBRs calculated using (28).Again, the power reference of the average IBR in the aggregated model is the average of the power references of all the IBRs, according to (32). Figure 13 shows the EMT simulation results of the active and reactive powers at the POM obtained by using the detailed and aggregated models, and Table 3 gives the THD comparison.As expected, the difference between results obtained by using the aggregated and detailed models is higher for Case 3 than that for Case 2. However, the difference is still very small, demonstrating the effectiveness of the proposed aggregated model for Case 3.
Energies 2022, 15, x FOR PEER REVIEW 16 of 20 Again, the power reference of the average IBR in the aggregated model is the average of the power references of all the IBRs, according to (32). Figure 13 shows the EMT simulation results of the active and reactive powers at the POM obtained by using the detailed and aggregated models, and Table 3 gives the THD comparison.As expected, the difference between results obtained by using the aggregated and detailed models is higher for Case 3 than that for Case 2. However, the difference is still very small, demonstrating the effectiveness of the proposed aggregated model for Case 3.

Approximated IBR Plant Model Evaluation of Case 4
This case is like Case 3, except that the IBRs are connected to the POM in both parallel and series ways.Therefore, the approximate collector resistance and inductance of the aggregated model are obtained using (30) instead of (29). Figure 14 presents the EMT simulation results of the active and reactive powers at the POM obtained by using the detailed and aggregated models, and Table 3 gives the THD comparison.For this case, a comparison of the voltage and current waveforms using the detailed and aggregated models are also shown in Figure 14.As it can be seen, the difference between results obtained by using the aggregated and detailed models is still small, illustrating the effectiveness of the proposed aggregated model development for Case 4.

Approximated IBR Plant Model Evaluation of Case 4
This case is like Case 3, except that the IBRs are connected to the POM in both parallel and series ways.Therefore, the approximate collector resistance and inductance of the aggregated model are obtained using (30) instead of (29). Figure 14 presents the EMT simulation results of the active and reactive powers at the POM obtained by using the detailed and aggregated models, and Table 3 gives the THD comparison.For this case, a comparison of the voltage and current waveforms using the detailed and aggregated models are also shown in Figure 14.As it can be seen, the difference between results obtained by using the aggregated and detailed models is still small, illustrating the effectiveness of the proposed aggregated model development for Case 4.

Performance Evaluation of the Approximate IBR Plant Model during Abnormal Operation Conditions
In this subsection, the performance of the proposed IBR plant during a fault event in the grid is analyzed.A real world event reported by Electric Reliability Council of Texas (ERCOT) in [35] is taken into account, where a type-4 WPP is connected to the grid, and a three phase to ground fault occurs, resulting abnormal operation in the WPP.To perform the analysis, the most comprehensive model, with the IBRs arranged in series and parallel combinations (explained in case 4), is considered.The IBR plant consists of 600 online IBRs connected via a 25k V collector system, and then interfaced to the grid via two 120 kV transmission lines.While operating normally, the short circuit ratio at the point of measurement (POM) is 3.5, which means the grid strength is in the medium range.A three phase to ground fault is introduced to one of the transmission lines, which means one of the transmission lines is disconnected from the grid.Due to the event, the impedance of the transmission line increases, and the short circuit ratio at POM decreases to 1.75 (weak grid).In this situation, the response of the IBR plant depends on the power supply scenario of the IBR plant.In Figure 15a, the response of the IBR plant at the POM is shown when the IBR plant's real and reactive power supply to the grid is low.The result is also compared with the detailed model to demonstrate the proposed model performance.In Figure 15b, the response with the higher real and reactive power supply is provided.In this case, an undamped oscillation can be noticed after the fault is introduced.This oscillation is denoted as sub-synchronous oscillation, and any further change in the plant, as well as in the grid, can lead to IBRs being tripped, resulting power loss of the IBR plant.A similar result, found in the detailed model simulation, is also shown in Figure 15b.These results are consistent with the results presented in the real case results of [35].

Performance Evaluation of the Approximate IBR Plant Model during Abnormal Operation Conditions
In this subsection, the performance of the proposed IBR plant during a fault event in the grid is analyzed.A real world event reported by Electric Reliability Council of Texas (ERCOT) in [35] is taken into account, where a type-4 WPP is connected to the grid, and a three phase to ground fault occurs, resulting abnormal operation in the WPP.To perform the analysis, the most comprehensive model, with the IBRs arranged in series and parallel combinations (explained in case 4), is considered.The IBR plant consists of 600 online IBRs connected via a 25k V collector system, and then interfaced to the grid via two 120 kV transmission lines.While operating normally, the short circuit ratio at the point of measurement (POM) is 3.5, which means the grid strength is in the medium range.A three phase to ground fault is introduced to one of the transmission lines, which means one of the transmission lines is disconnected from the grid.Due to the event, the impedance of the transmission line increases, and the short circuit ratio at POM decreases to 1.75 (weak grid).In this situation, the response of the IBR plant depends on the power supply scenario of the IBR plant.In Figure 15a, the response of the IBR plant at the POM is shown when the IBR plant's real and reactive power supply to the grid is low.The result is also compared with the detailed model to demonstrate the proposed model performance.In Figure 15b, the response with the higher real and reactive power supply is provided.In this case, an undamped oscillation can be noticed after the fault is introduced.This oscillation is denoted as sub-synchronous oscillation, and any further change in the plant, as well as in the grid, can lead to IBRs being tripped, resulting power loss of the IBR plant.A similar result, found in the detailed model simulation, is also shown in Figure 15b.These results are consistent with the results presented in the real case results of [35].

Conclusions
This paper presents an aggregated dynamic transient simulation model of an IBR plant.The aggregated model has incorporated different aspects of an IBR plant, as well as IBR configurations within the plant, such as collector system impacts, grid-connected filters of the IBR inverter, detailed control configurations at the inverter level, series/parallel orientation of IBRs within the plant, etc.The model development is supported by a detailed mathematical analysis and proving process.The developed IBR plant model can be used to simulate transient performance at the individual IBR level, evaluate the detailed and complete design of the IBR controller and grid-connected filter, and evaluate the transient interactions of multiple IBR groups in a plant structure.At the same time, the model can be used to evaluate the transient performance for the grid interconnection of the plant based on the IBR system and controller parameters of individual IBRs, instead of scaled and modified parameters that were normally needed to fit an aggregated plant model.The proposed model development has carefully addressed the collector system impact in diverse parallel and series connection scenarios.The model is verified through a comparison study with a detailed IBR plant model in a recursive way.The study and evaluation results show that the proposed IBR plant model is effective and accurate in investigating the interconnection of an IBR plant, with the grid as well as IBR controllers, at the inverter level in a plant configuration scenario.

Conclusions
This paper presents an aggregated dynamic transient simulation model of an IBR plant.The aggregated model has incorporated different aspects of an IBR plant, as well as IBR configurations within the plant, such as collector system impacts, grid-connected filters of the IBR inverter, detailed control configurations at the inverter level, series/parallel orientation of IBRs within the plant, etc.The model development is supported by a detailed mathematical analysis and proving process.The developed IBR plant model can be used to simulate transient performance at the individual IBR level, evaluate the detailed and complete design of the IBR controller and grid-connected filter, and evaluate the transient interactions of multiple IBR groups in a plant structure.At the same time, the model can be used to evaluate the transient performance for the grid interconnection of the plant based on the IBR system and controller parameters of individual IBRs, instead of scaled and modified parameters that were normally needed to fit an aggregated plant model.The proposed model development has carefully addressed the collector system impact in diverse parallel and series connection scenarios.The model is verified through a comparison study with a detailed IBR plant model in a recursive way.The study and evaluation results show that the proposed IBR plant model is effective and accurate in investigating the interconnection of an IBR plant, with the grid as well as IBR controllers, at the inverter level in a plant configuration scenario.

Figure 1 .
Figure 1.A generalized diagram of (a) a full-scale IBR system and (b) a grid connected WPP.

Figure 1 .
Figure 1.A generalized diagram of (a) a full-scale IBR system and (b) a grid connected WPP.

Figure 2 .
Figure 2. Plant level and IBR level controller structure of an inverter interfaced IBR plant.

Figure 2 .
Figure 2. Plant level and IBR level controller structure of an inverter interfaced IBR plant.

Figure 2 .
Figure 2. Plant level and IBR level controller structure of an inverter interfaced IBR plant.

3. 1 .
A Grid-Connected IBR Plant Including IBRs, IBR Filters, Collector Impedances, and Grid Impedance To identify IBRs in the IBR plant, superscript k is used, where k = 1, 2,..., n. n defines the total number of online IBRs in the IBR plant.

Figure 5 .
Figure 5. Equivalent circuit diagram of an IBR plant looking from the IBR side.

Figure 5 .
Figure 5. Equivalent circuit diagram of an IBR plant looking from the IBR side.

Figure 6 .
Figure 6.Equivalent circuit diagram of an IBR plant looking from the grid side.

Figure 6 .
Figure 6.Equivalent circuit diagram of an IBR plant looking from the grid side.

3. 2 .
Case 2: IBRs with Identical Filter Parameters but Different Terminal Voltages Assume all IBRs have different output terminal voltages.As the collector system impedance is neglected, PCC voltages of all IBRs are the same and equivalent to the POM voltage.

3. 3 .
Case 3: IBRs with Different Filter Parameters and Different Terminal Voltages-Each IBR Is Coupled in Parallel to the POM through a Different Collector Impedance

Figure 7 .
Figure 7. Equivalent circuit diagram of an IBR plant: (a) looking from the IBR side, (b) looking from the average IBR viewpoint.

3. 4 .Figure 7 .
Figure 7. Equivalent circuit diagram of an IBR plant: (a) looking from the IBR side, (b) looking from the average IBR viewpoint.

3 )
The collector system parameters (

Figure 7 .
Figure 7. Equivalent circuit diagram of an IBR plant: (a) looking from the IBR side, (b) looking f the average IBR viewpoint.

3. 4 .
Case 4: IBRs with Different Filter Parameters and Different Terminal Voltages-Each IBR Is Coupled in a Combination of Series and Parallel to the POM through a Collector System (As Depicted in Figure 1b)

)t
are the active and reactive power references of each individual IBR (i = 1, …, n).

Figure 10 .
Figure 10.Configuration of the aggregated IBR plant model: (a) IBR plant structure, (b) IBR control structure.

Figure 10 .
Figure 10.Configuration of the aggregated IBR plant model: (a) IBR plant structure, (b) IBR control structure.

Figure 11 .
Figure 11.Active and reactive power at the POM for Case 1 using aggregated and detailed models: (a) Active power, (b) Reactive Power.

Figure 12 .
Figure 12.Active and reactive power at the POM for Case 2 using aggregated and detailed models: (a) Active power, (b) Reactive Power.

Figure 11 .
Figure 11.Active and reactive power at the POM for Case 1 using aggregated and detailed models: (a) Active power, (b) Reactive Power.

Figure 11 .
Figure 11.Active and reactive power at the POM for Case 1 using aggregated and detailed models: (a) Active power, (b) Reactive Power.

Figure 12 .
Figure 12.Active and reactive power at the POM for Case 2 using aggregated and detailed models: (a) Active power, (b) Reactive Power.5.3.Approximated IBR Plant Model Evaluation of Case 3 In this case, all the IBRs within the IBR plant have different LCL filter parameters, deliver different IBR output voltages, and are connected in parallel with different collector impedances to the POM.For the detail model, a perturbation of 20% on filter inductance

Figure 12 .
Figure 12.Active and reactive power at the POM for Case 2 using aggregated and detailed models: (a) Active power, (b) Reactive Power.

Figure 13 .
Figure 13.Active and reactive power at the POM for Case 3 using aggregated and detailed models: (a) Active power, (b) Reactive Power.

Figure 13 .
Figure 13.Active and reactive power at the POM for Case 3 using aggregated and detailed models: (a) Active power, (b) Reactive Power.

Figure 14 .
Figure 14.Active and reactive power and current at the POM for Case 4 using aggregated and detailed models: (a) Active power, (b) Reactive Power, (c) Three-phase current plus a comparison of Phase a current.

Figure 14 .
Figure 14.Active and reactive power and current at the POM for Case 4 using aggregated and detailed models: (a) Active power, (b) Reactive Power, (c) Three-phase current plus a comparison of Phase a current.

Figure 15 .
Figure 15.IBR plant response in case of three phase to ground fault in transmission line (a) strong grid scenario, (b) weak grid scenario.

Figure 15 .
Figure 15.IBR plant response in case of three phase to ground fault in transmission line (a) strong grid scenario, (b) weak grid scenario.

Table 1 .
Summary of IBR models presented in literature.

Table 2 .
System parameters of the IBR plant model.

Table 3 .
THD comparison between aggregate and detail model at POM.