Development and validation of a real-time testbed for renewable energy integration studies: South African grid code case study

Gideon Daniel Joubert (Department of Electrical, Electronic and Computer Engineering, Cape Peninsula University of Technology – Bellville Campus, Cape Town, South Africa)
Atanda Kamoru Raji (Department of Electrical, Electronic and Computer Engineering, Cape Peninsula University of Technology – Bellville Campus, Cape Town, South Africa)

Journal of Engineering, Design and Technology

ISSN: 1726-0531

Article publication date: 14 April 2023

926

Abstract

Purpose

Despite South Africa’s ailing electrical grid, substantial renewable energy (RE) integration is planned for the country. As grid-integrated RE affects all grids differently, this study aims to develop an adaptable grid code-guided renewable power plant (RPP) control real-time simulation testbed, tailored to South African grid code requirements to study grid-integrated RE’s behaviour concerning South Africa’s unique conditions.

Design/methodology/approach

The testbed is designed using MATLAB’s Simulink and live script environments, to create an adaptable model where grid, RPP and RPP guiding grid codes are tailorable. This model is integrated with OPAL-RT’s RT-LAB and brought to real-time simulation using OPAL-RT’s OP4510 simulator. Voltage, frequency and short-circuit event case studies are performed through which the testbed’s abilities and performance are assessed.

Findings

Case study results show the following. The testbed accurately represents grid code voltage and frequency requirements. RPP point of connection (POC) conditions are consistently recognized and tracked, according to which the testbed then operates simulated RPPs, validating its design. Short-circuit event simulations show the simulated wind farm supports POC conditions relative to short-circuit intensity by curtailing active power in favour of reactive power, in line with local grid code requirements.

Originality/value

To the best of the authors’ knowledge, this is the first design of an adaptable grid code-guided RPP control testbed, tailored to South African grid code requirements in line with which RPP behavioural and grid integration studies can be performed.

Keywords

Citation

Joubert, G.D. and Raji, A.K. (2023), "Development and validation of a real-time testbed for renewable energy integration studies: South African grid code case study", Journal of Engineering, Design and Technology, Vol. ahead-of-print No. ahead-of-print. https://doi.org/10.1108/JEDT-07-2022-0380

Publisher

:

Emerald Publishing Limited

Copyright © 2023, Gideon Daniel Joubert and Atanda Kamoru Raji.

License

Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial & non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence may be seen at http://creativecommons.org/licences/by/4.0/legalcode


Nomenclature

F

= frequency, Hz;

HFRT

= high-frequency ride-through;

HVRT

= high-voltage ride-through;

I

= current, A;

LFRT

= low-frequency ride-through;

LVRT

= low-voltage ride-through;

MATLAB

= programming platform used for designing and analyzing systems;

OP4510

= OPAL-RT real-time simulator;

OPAL-RT

= a leading developer of real-time simulation software and simulators;

P

= power, W;

pu

= per unit system, quantities are expressed as a ratio relative to a common base;

POC

= Point of connection;

Q

= reactive power, var;

RE

= renewable energy;

rad

= radians, a unit of angle;

RPP

= renewable power plant;

S

= apparent power, VA;

SC_Target

= model console subsystem;

TCP/IP

= connection between computer and simulator;

V

= voltage, V; and

Z

= impedance, ohm.

Symbols

|u|

= absolute value;

θ

= angle, degrees; and

π

= pi, 3.14159.

Subscripts

abc

= three-phase value;

active

= active component of a value;

B3

= Bus 3/POC;

b

= base value;

in

= input;

new

= modified value;

nu

= upper grid code boundary value of an RPPs continuous operating range;

nl

= lower grid code boundary value of an RPP’s continuous operating range;

out

= output;

old

= original value;

peak

= peak value;

ph

= phase measurement;

reactive

= reactive component of a value;

spreadsheet

= imported value;

wt

= wind turbine.

1. Introduction

Globally, an energy shift is underway where renewable energy (RE) alternatives are increasingly being favoured over their fossil fuel-driven counterparts, in a bid to reduce both harmful emissions, and reliance on exhaustible energy resources. Leading the way for the electrical generation industry is hydropower, wind, biomass and solar technologies, owing to their versatility and widespread application. These applications include small-scale implementation of technologies such as the vertically orientated Savonius-style wind turbine optimized by Roy et al. (2018) for agricultural, household and rural applications. The study of Goswami and Das (2020) investigates another such implementation of a waste heat-driven biomass engine, capable of charging a 12 V battery, making it a valuable and implementable design for various small-scale applications.

RE generation’s more significant contribution to the displacement of its fossil fuel-driven counterparts however comes in the form of large-scale grid integration, bringing with it added challenges to existing electrical grids. Li et al. (2017) investigate this, noting the most significant challenges are increased system voltage levels, increased short-circuit currents, deterioration in power supply reliability and power quality, highlighting the importance of a strong stable grid for successful large-scale RE integration. A major contributing factor is noted to be RE’s lacking contribution to system inertia, which refers to the loss of conventional generation’s spinning reserves. Saha et al. (2023) investigate the impact of this, highlighting the effect loss of system inertia has on grid frequency stability, and the added challenges this poses to system operators in terms of maintaining system frequency within prescribed limits. Concerns surrounding the loss of system inertia are supported by Mararakanye and Bekker (2019), also noting that loss of inertia implications greatly vary between power grids, making it essential to consider the circumstances of a specific region.

This paper focuses on the grid integration of RE in South Africa, which has recently experienced several challenges surrounding supply reliability of the country’s energy needs, leading to a decline in the average Energy Availability Factor (EAF) down to only 65% in 2020 (Hanto et al., 2022). This follows years of infrastructure neglect and lack of maintenance leading Eskom, the country’s primary energy provider, to battle unplanned breakdowns of its old and poorly maintained power plants to remain below 9,500 MW, at which point planned load shedding is initiated to maintain stability of the country’s strained electrical grid (Statistics South Africa, 2018, 2019, 2020; Eskom, 2019b, 2019a; Eskom generation division, 2019; Maune, 2019; Creamer, 2021; Ayamolowo et al., 2022). As part of the solution to the ongoing energy crises South Africa is planning to significantly ramp up its use of RE for electricity production, and therefore developed RE specific grid codes to potentially limit added instabilities brought about by its integration. RE integration plans are reflected by initiatives such as the 2030 Integrated Resource Plan (IRP), setting a renewable generation addition target of 14,725 MW by 2030, 6,422 MW of which are already operational and around 4,100 MW of which are already connected and grid integrated (Department of Energy, 2010, 2011; National Energy Regulator of South Africa (NERSA), 2012; Department of Energy and Department of National Treasury, 2018; RSA Department of Energy, 2018; Department of Mineral Resources and Energy, 2019; Statistics South Africa, 2020; Ayamolowo et al., 2022; Hanto et al., 2022).

Considering South Africa’s unique implementation of RE as aid for their already struggling electrical grid, and the foreseeable challenges associated with large-scale RE integration, South Africa is set to greatly benefit from tailored RE integration studies. Given the obligation of grid-connected RE to adhere to renewable power plant (RPP)-specific grid code requirements in South Africa, it would be advantageous to the accuracy of RE integration studies to include grid code requirements as part of such simulations. Paquin et al. (2007) developed a MATLAB-integrated OPAL-RT real-time simulation model to perform detailed studies of doubly fed induction generator wind farms connected to a 24-bus electrical network, noting the ability of such a platform to study various aspects of RE’s performance. Merabet et al. (2014) developed a wind turbine-generator emulating system, implementing a MATLAB-integrated OPAL-RT real-time model to simplify control scheme development of both generator and grid-side converters. Li et al. (2016) developed an analogue digital hybrid MATLAB integrated OPAL-RT real-time simulation platform for single- or multi-generation unit grid connected behavioural simulations, noting the feasibility of such platforms for this purpose. Kemal et al. (2017) developed a model using MATLAB/Simulink and OPAL-RT to demonstrate the abilities of these tools to perform hardware-in-the-loop simulations when studying various aspects of the RE integration process.

Although several real-time simulation platforms to study various aspects of the RE integration process have been developed, none presents a tailorable model incorporating grid code requirements against which simulated RPPs can be operated when performing RE grid integration studies. This study, therefore, develops a versatile grid code-guided RPP control testbed tailored to South African grid code requirements for performing RE grid integration studies.

It furthermore became apparent that MATLAB-integrated OPAL-RT models are favoured for their power generation real-time simulation abilities. Noureen et al. (2018) investigate and summarize the strengths and abilities of OPAL-RT’s renowned simulators, and the various advantages of implementing MATLAB as front-end software. The validity of this combination for power generation modelling is corroborated by Singh et al. (2015), modelling the New England 39 bus, and WSCC 9 bus systems to assess the performance and limitations of both model and simulator under various test cases. Given this, the testbed of this study will implement MATLAB Simulink for designing the model, which will then be brought to real-time simulation using OPAL-RT’s RT-LAB, and OP4510 real-time simulator. The design incorporates grid code requirements against which a simulated RPP connected to a grid representing network is operated, concerning point of connection (POC) conditions. This method of operation then allows the strain placed on RPPs by grid code requirements to be studied, in addition to the simulation abilities offered by other RE grid integration platforms reviewed.

Though tailored to South African specifications for this study, the testbed’s design will also allow the user to add MATLAB live script contained grid code requirements, replay any grid representing spreadsheet contained data or insert any type and size of RPP, enabling the testbed’s applications to reach well beyond that of South Africa. The option of integrating additional grid code requirements will then allow the testbed to be implemented for studies like that of Al Riyami et al. (2016), simulating a 50 MW Dhofar network integrated wind farm to evaluate new RPP grid code requirements for the region. Studies such as that of Abdalla (2016, 2018), in which technical requirements are reviewed for PV systems connected to Egypt’s electrical network, can then also benefit from adopting such a testbed to perform grid-connected RPP behavioural studies, implementing the technical specifications it reviewed.

The paper is structured as follows. Section 2 summarizes the relevant South African grid code requirements integrated as part of the testbed’s design. Section 3 discusses the model development process, including the circuit sections and their respective inputs and outputs. Case study results and the discussion thereof are covered in Section 4 before Section 5 concludes and offers suggestions for future work. Based on the outcome of this study, its main contribution is an adaptable real-time grid code-guided RPP control testbed, to be used for RE grid integration studies. The testbed may then be tailored and implemented by students and researchers alike to simulate the behaviour of grid connected RPPs at a specific POC, under various replayed or recreated conditions.

2. A review of testbed integrated South African renewable power plant grid codes

Grid codes are important to the testbed’s development as their requirements govern conditions during which RPPs should remain connected, allowing the demands placed on RPPs by these requirements during disturbances to be studied. As South Africa has already developed RPP-specific grid codes with future stability and reliability in mind, the testbed’s development proceeded to incorporate these as what would be required of RPPs to ensure future stability and reliability leading up to the integration of 14,725 MW renewables by 2030.

South Africa’s RPP grid code requirements however accommodate different types and sizes of RPPs, an ability that had to be incorporated into the testbed. Determining the requirements of a specific RPP is its category based on size, classified as follows (Yan et al., 2015; National Energy Regulator of South Africa, 2016; Etxegarai et al., 2017; Department of Energy and Department of National Treasury, 2018; RSA Department of Energy, 2018; Mararakanye and Bekker, 2019; Mchunu, NERSA and Eskom Transmission Division, 2019):

  • Category A1: 0 to 13.8 kVA;

  • Category A2: > 13.8 kVA to < 100 kVA;

  • Category A3: 100 kVA to < 1 MVA;

  • Category B: 1 MVA to < 20 MVA; and

  • Category C: 20 MVA or greater.

With an RPP’s details known, its grid-connected requirements are identified from the South African RPP grid code document. Although grid codes include a wide range of support and ride-through specifications, the testbed’s focus will fall on incorporating grid code voltage and frequency specifications, being fundamental requirements. These are given by grid codes in the form of graphs, such as that in Figure 1, illustrating requirements for non-synchronous Category C RPPs (Mchunu, NERSA and Eskom Transmission Division, 2019).

Considering Figure 1, the voltage and frequency requirement graphs can essentially be divided into three areas of RPP operation. The first of these is the “Continuous Operating Range” respectively indicated on both graphs. When in this region of operation, it is expected that the RPP remains connected indefinitely and irrespective of the number of fluctuations within this band. The second area falls outside the continuous operating band, but within the outer boundaries of each graph. This area is known as the fault ride-through region of operation, where it is expected that the RPP remains connected, riding through the fault in support of the grid, thereby allowing time for a fault to potentially clear. The final area of operation is the trip region outside the respective graph boundaries, at which point it is expected of an RPP to disconnect from the grid to prevent or limit any potential damage (Mchunu, NERSA and Eskom Transmission Division, 2019).

3. Model development

The section focuses on the model developed as a real-time testbed, elaborating on the inputs, sections of the main circuit, subcircuits responsible for validating South African RPP grid codes in terms of voltage and frequency and the outputs generated to allow assessment of the studied RPP’s performance. As this study develops a real-time testbed using a MATLAB integrated OPAL-RT model, inputs and outputs will need to be split from the main circuit, allowing the model to be excecated in real-time using RT-LAB and the OPAL-RT OP4510 real-time simulator. A TCP/IP cable is used to establish a physical wired connection using the local network between the equipment, while communication is set up and configured using the host computer and RT-LAB software (Buraimoh and Davidson, 2022). The split layout of the model design is given by the diagram in Figure 2, which also includes the respective host computer and OP4510 real-time simulator used in the lab.

Although designed with versatility and adaptability in mind, this study’s simulated RPP will be limited to that of a wind farm operated at optimal wind speed for testbed testing purposes. Furthermore, model-integrated grid codes only include voltage and frequency requirements set for RPPs in South Africa, while fabricated data will be replayed during test simulations to achieve desired POC conditions for testing. Results and conclusions drawn are therefore oriented towards testing of the testbed and are less focussed on conducting an RE grid integration study as part of this paper. Figure 3 follows, illustrating the testbed’s main circuit and flow diagram.

3.1 Data imports

Part of developing the circuit was to make it adaptable, allowing any type and size of RPP to be inserted and studied. As different category RPPs need to adhere to different requirements, variable inputs had to be obtained based on the simulated RPP, which include the normal upper (Vnu) and lower (Vnl) operating voltage boundaries. For increased versatility, the circuit uses a single source with inputs “Vpu spreadsheet” and “f spreadsheet” to represent the larger grid network, allowing the testbed to replay previously recorded grid conditions accurately and easily as compared to an elaborate static grid representing circuit. All inputs additionally pass through an OpComm block indicated in Figure 3, allowing the input data housed on the host computer to be imported and used on the physically separated real-time simulator. The subsections, inputs and outputs can then further be defined below.

3.1.1 Vnu and Vnl.

The “normal upper” (Vnu) and “normal lower” (Vnl) voltage boundary inputs represent the continuous operating voltage limits and are specified according to an RPP’s class and POC voltage. For the Eskom Sere representing wind farm modelled in the study, these are respectively represented by the upper and lower boundaries of the “Continuous operating range” indicated as Area A in Figure 1, discussed in Section 2. A MATLAB live script is also implemented, requiring the user to select the simulated RPP’s category and POC voltage, after which “if” functions select and output the appropriate “Vnu” and “Vnl” values.

3.1.2 Vpu spreadsheet and f spreadsheet.

Vpu spreadsheet” and “f spreadsheet”, respectively, represent the pu voltage and frequency in Hz inputted to the grid representing source as a value concerning time. This data is housed and imported from an Excel document, making it easy to replay past measured grid data and events, or input custom data for testing purposes. As this study aims to demonstrate the abilities and accuracy of the developed testbed, custom data are used to recreate the desired scenarios of the respective case studies.

3.2 Model components

3.2.1 Grid representing source.

The “grid representing source” is seen as the first significant component forming part of the grid section labelled in Figure 3 and uses the imported “Vpu spreadsheet” and “f spreadsheet” data to generate a usable three-phase output to the circuit. To achieve this, a base voltage and nominal frequency of 25 kV at 50 Hz is selected, from which this source block then generates the three-phase output according to equation (1), given by Dixon (2019):

(1) Vpu=VoutVbase1

Consequently, when applied to the grid representing source, Vpu = Vpu spreadsheet, Vbase1 =the transformer’s primary side voltage, and Vout = the grid representing source output voltage.

Therefore,

(2) Vout=Vpu spreadsheed×Vbase1

However, as in the end, the per unit voltage is of interest at the POC, being the primary input to the grid code validation subsystems, it can be stated that the grid representing source output is simply given by equation (3):

(3) Vout(pu)=Vpu spreadsheed

3.2.2 Transformer.

Connected to the source through Bus 1 is an ideal 25 kV/690 V step-down wye-delta transformer, of which the secondary voltage was chosen to match that of the simulated RPP. In terms of per unit values, the transformer’s circuit impedance (ZTnew) can be represented by equation (4), as given by Müller et al. (2009):

(4) ZTnew=Zold(SB_newSB_old)(VB_oldVB_new)2

where SB_new and VB_new represent chosen base circuit values, while SB_old, VB_old and Zold, respectively, represent the transformer voltage, power and impedance rating.

3.2.3 Load.

Connected to the transformer through Bus 2 is a 100 MW grounded static local load. This can be modelled concerning per unit values using equation (5) as in Müller et al. (2009):

(5) Zloadpu= LoadZB2

where ZB2 is given by equation (6):

(6) ZB2=VB2SB

And, VB2 = the secondary transformer voltage chosen as base value, while SB = the base power rating chosen for the circuit.

3.2.4 Renewable power plant.

Given that the RPP modelled in Figure 3 is based on Eskom’s Sere wind farm, it consists of 46 × 2.3 MW asynchronous wind turbines with an output voltage of 690 V. To maintain consistency in validating the respective model and validation subsystems, wind speed is assumed constant at 8 m/s for this study. Like the grid representing source, the RPP’s voltage (Vpu(RPP)) can be represented in terms of the per unit system by equation (7), as given in Dixon (2019):

(7) Vpu(RPP)=Vout(RPP)Vbase2

where Vout(RPP) = generated RPP terminal voltage, and Vbase2 = the secondary voltage of the transformer chosen as base value.

3.2.5 Point of connection.

The POC labelled in Figure 3 represents the connection point of the studied RPP to the rest of the network and is of particular significance, being that all measurements are obtained from this bus unless specified otherwise. This includes the complex per unit voltage and current measurements used by the voltage and frequency validation subsystems. Consequently, considering Figure 3, power flows through the POC feeding the load and/or grid representing source depending on simulated conditions. POC per unit current (Ipu(POC)) and voltage (Vpu(POC)) can therefore be given by equations (8) and (9):

(8) Ipu(POC)=Ipu(Load)Ipu(Grid)
(9) Vpu(POC)=Ipu(POC)×Zcircuit

where Ipu(Load) is the current flowing to the load, Ipu(Grid) is the current flowing from or to the grid representing source and Zcircuit is the impedance of the circuit concerning the POC.

3.2.6 Voltage validation subsystem.

The voltage validation subsystem labelled in Figure 3 consists of three, respective, parts, namely, “Detect Event”, “Measure Event Time” and “Voltage Violation”. The first of these, “Detect Event”, use “Vnu”, “Vnl” and the measured POC bus voltage “V_POC” as inputs to detect when POC conditions deviate outside normal Vnu/Vnl operating boundaries. When a voltage event is present, a “1” is generated at the “LV/HV Event” output and passed on to the adjacent “Measure Event Time” section. Using the “LV/HV Event” output as input to the “Measure Event Time” section, the duration of the event is measured and outputted in seconds at this section’s “Event Time” output. Using this as input, the “Voltage Validation” section continually compares the event duration to the grid code voltage requirement graph boundaries of Figure 1 using lookup tables, allowing the maximum (Vu out) and minimum (Vl out) instantaneous POC voltage boundaries to be generated as outputs. Comparing “Vu out” and “Vl out” to V_POC then allows grid code voltage violations to be identified, at which time a “1” is generated at the section’s “V Violation” output.

3.2.7 Frequency validation subsystem.

The purpose of the frequency validation subsystem labelled in Figure 3 is to validate compliance set for RPPs by South African grid codes concerning frequency. As frequency requirements remain constant irrespective of RPP class and POC voltage, the frequency validation subsystem does not require external parameters to be imported into the model. Frequency can however not directly be measured at the POC given the type of simulation, and therefore had to be obtained using the complex POC voltage (V_POC), also used as input to the voltage validation subsystem. Considering the tools available in the MATLAB library, POC frequency could be determined using a “frequency (Phasor)” block, if the phase angle in radians “phase(rad)” is known. Using V_POC, this was then determined. First, the complex POC voltage (V_POC) measured at the POC in the format given by equation (10) had to be converted to a more appropriate magnitude (|u|) angle (θ) format, respectively, using equations (11) and (12), as given by the Centre for Innovation in Mathematics Teaching (2009):

(10) Vpoc=a+bi
(11) |u|=|a2+b2|
(12) θ=arctanba

The angle (θ) in degrees could then be converted to radians using equation (13), and used as input to the “frequency (Phasor)” MATLAB library block with output F_POC representing the POC frequency (Centre for Innovation in Mathematics Teaching, 2009):

(13) rad=π180

Next, the grid code frequency requirements of Figure 1 were considered, and a circuit was created implementing comparators to compare F_POC to the constant frequency boundaries of Figure 1, also applying the appropriate time delays where necessary. The comparators will consequently output a “1” when their respective frequency boundary is violated, which after having been delayed according to the graph of Figure 1 is passed on using the “F Violation” output of the frequency validation subsystem.

3.2.8 Renewable power plant data acquisition subsystem.

The RPP data acquisition subsystem was included to obtain the data needed for RPP short-circuit testing specified by South African grid codes. Active (I_active) and reactive (I_reactive) current consequently had to be obtained using the available complex three-phase per unit voltage “V_POC” and current “I_POC” available. To achieve this active (P) and reactive (Q) power were first determined using a MATLAB library circuit, and the known relation given by equations (14) and (15), as in Müller et al. (2009):

(14) P=3×Vph×Iph×cosϕ

And

(15) Q=3×Vph×Iph×sinϕ

Using the active (P) and reactive (Q) power outputs, the active (I_active) and reactive (I_reactive) current components were determined by adapting equation (16) to find equations (17) and (18) (All about circuits, 2020):

(16) P=I×V
(17) Iactive=PV_POC
(18) Ireactive=QV_POC

3.3 Circuit outputs

Voltage simulation results and measurements are obtained and logged using the MATLAB scopes displayed in Figure 4. As these measurements are obtained from the SC_Target subsystem, the data of all scopes passes through an OpComm block before entering the respective scope. The scopes and their respective inputs discussed throughout can then further be defined as follows.

3.3.1 Voltage measurements scope.

The Voltage measurements scope of Figure 4 is used to output three graphs, the first displaying the imported per unit spreadsheet voltage “Vpu spreadsheet”, the second the voltage validation subsystem output “V Violation”, while the third gives the per unit POC voltage “V_POC” alongside the respective voltage boundary limits, “Vnu”, “Vnl”, “Vu out” and “Vl out” given by Figure 1.

3.3.2 Frequency measurements scope.

The frequency measurements scope of Figure 4 again outputs three graphs, the first displaying the imported spreadsheet frequency “f spreadsheet” in Hz, the second the frequency validation trip signal “F Violation”, while the third gives the POC frequency “F_POC” alongside the constant normal, and frequency ride-through boundary limits, “Fnu”, “Fnl”, “Lower f limits” and “Upper f limits” given by Figure 1.

3.3.3 Console subsystem (SC) test scope.

The SC Test scope of Figure 4 is used to output four graphs, namely, the active (I_active), reactive (I_reactive) and peak (I_peak) currents, determined by the RPP data acquisition subsystem, in addition to the final graph used for the POC voltage (V_POC), and RPP trip input (Trip in).

3.3.4 Wind turbine measurements scope.

The final wind turbine measurements scope of Figure 4 is used to display four graphs, indicating the per unit wind turbine voltage (Vabc_wt), current (Iabc_wt), active power generated in MW (P_out), reactive power generated in MVAR (Q_out) and the trip input signal (Trip in). These signals are all obtained or calculated using direct wind turbine measurement outputs, included to illustrate the load put on the simulated RPP under recreated conditions.

4. Case study results and discussion

Case studies were selected aimed at testing the ability of the model to appropriately respond to voltage and frequency events when feeding test data with a known outcome. This way, the effectiveness and accuracy of the respective circuit sections could be assessed while illustrating their ability to respond to replayed test conditions as expected.

4.1 Voltage validation subsystem testing

To test the ability of the voltage validation subsystem to assess and respond in line with RPP grid code voltage requirements, three case studies were carried out mimicking events that would cause POC conditions to operate within the three regions of RPP operation. These include the continuous operating region indicated by the blue-dotted boundaries of each graph, the fault ride-through region given by the green-dotted boundaries and the trip region outside the green-dotted ride-through boundaries. Looking at Figures 5 to 7, it should be noted that the green-dotted ride-through boundaries changes with time, as it adjusts from the initiation of a fault concerning the voltage ride-through requirement graph of Figure 1. From the results, it is expected to observe that the RPP continues operation while operating within the respective continuous and ride-through regions, resetting for faults with quick recovery times and tripping once ride-through boundaries are violated. Figures 5 to 7 illustrate the results obtained from the respective voltage case study simulations.

In Case study 1 to 3, Figures 5 to 7, the various degrees of voltage spikes and dips introduced to assess voltage validation subsystem response is observed. From the results, it is noted that some distortion in V_POC occurs at the initiation and termination of events, as the RPP compensates and adjusts to the abrupt V_POC changes introduced. In terms of the voltage validation subsystem’s response, it is observed that while V_POC remain within the blue-dotted continuous operating region, no trip is generated and the green-dotted fault ride-through boundaries does not respond. A fault ride-through boundary response indicating abnormal V_POC conditions is observed in Case studies 2 and 3, as V_POC exits the blue-dotted continuous operating region. The voltage validation subsystem responds by adjusting green HVRT/LVRT boundaries, in line with Figure 1’s grid code voltage requirements it represents. For Case study 2 the fault is short lived, and V_POC recovers in time to reset the subsystem and avoid a trip. For Case study 3 the fault persists, violating the green-dotted fault ride-through boundaries, at which point a trip is generated. Based on the voltage validation subsystem’s response to Case study 1 to 3’s events, it is concluded that POC conditions were effectively assessed, that the subsystem adapted HVRT/LVRT boundaries in line with Figure 1’s grid code voltage requirements, and that the subsystem was able to operate the simulated RPP accordingly.

4.2 Frequency validation subsystem testing

Like that of voltage, the frequency validation subsystem’s ability to respond to the three regions of operation will be assessed through three case studies mimicking events causing conditions to operate the RPP within these regions. These include the continuous operating region indicated by the blue-dotted boundaries of each graph, the fault ride-through region given by the adjusting green dotted boundaries, and the trip region outside the green-dotted ride-through boundaries. Similarly, it is expected of the RPP to remain connected indefinitely whilst within the continuous operating region, ride-through the disturbance while within the green-dotted ride-through boundaries, and trip once these are violated. Figures 8 to 10 illustrate the results obtained from the respective frequency case study simulations.

Like voltage, Case studies 4 to 6, Figures 8 to 10 introduced various degrees of frequency spikes and dips to assess frequency validation subsystem response to F_POC changes. From the results, the frequency validation subsystem is seen not to respond to F_POC variations within the blue-dotted continuous operating region. A fault ride-through response indicating abnormal F_POC conditions is observed for Case studies 5 and 6, as F_POC exits the blue-dotted continuous operating region. This is in the form of adjusting green HFRT/LFRT boundaries, in line with Figure 1’s grid code frequency requirements it represents. For Case study 5, although coming close, F_POC does not violate the green fault ride-through boundaries, recovering in time to reset the subsystem and avoid a trip. In Case study 6, the fault was severe enough to violate green-dotted ride-through boundaries, at which point a trip was generated. Based on the frequency validation subsystem’s response to Case study 4 to 6’s events, it is concluded that POC conditions were effectively assessed, that the subsystem adapted HFRT/LFRT boundaries in line with Figure 1’s grid code frequency requirements, and that the subsystem was able to operate the simulated RPP accordingly.

4.3 Short-circuit testing

Case study 7 was performed in response to mandatory short-circuit testing requirements of grid connected RPPs in South Africa. Specifications mandate the introduction of short circuits leading up to a 500 millisecond 20%, 30%, 50%, and 80% voltage drop at the POC, monitoring the RPP's current flow. During such a simulation, it is expected to observe a decrease in the active current flow in favour of an increased reactive current, thereby supporting recovery of the POC voltage. Furthermore, it is also expected that this response will become more prominent as short-circuit severity increase from 20% up to 80% voltage drop. Figure 11 follows, respectively showing the RPP's active and reactive current response for the conditions described.

Case study 7, Figure 11 confirms that once the respective short circuits are initiated at t = 0s, the RPP responds by significantly reducing the active current generated, in favour of reactive current, becoming more prominent with short-circuit intensity. For the yellow 80% short-circuit simulation, a trip is also generated as frequency ride-through requirements are violated. Case study 7 then shows the testbed to effectively respond during short circuits in line with what was expected, giving insight into the support demands placed on RPPs during such conditions.

5. Conclusion

In this paper, a versatile grid code-guided RPP control testbed is developed, tailored to South African grid code requirements for performing RE grid integration studies. The development of such a tool stems from the fact that RE integration affects all networks differently, making tailored simulations essential to the accuracy and quality of results produced. For this reason, the testbed was developed so that grid side, RPP and grid codes guiding RPP operation are tailorable, which in the case of this study, was tailored to South African conditions. Voltage, frequency and short-circuit event case studies were performed, evaluating the testbed’s design and abilities. Voltage and frequency case studies aimed at assessing testbed design validated its ability to accurately track POC conditions, allowing RPP operation in the continuous, fault ride-through and trip region to be identified, in line with which the testbed then operated the simulated RPP. Short-circuit case studies aimed at illustrating the testbed’s abilities as RE behavioural studies tool showed the testbed to provide valuable insight into grid code imposed demands, and grid integrated RPP behaviour. This includes that the simulated wind farm will adapt its output to provide reactive current support at the POC relative to short-circuit intensity. During the 80% short-circuit simulation, however, a grid code frequency violation causing the testbed to disconnect the RPP was observed, as the simulated wind farm was unable to provide sufficient POC support within prescribed grid code limits at this time. Studying an RPP’s support abilities, grid connected behaviour, fault response and the strain placed on RPPs by grid code requirements, is among the intended uses of the developed testbed ultimately forming part of RE grid integration studies to ensure the success of RE integration goals. Future work may then involve tailoring the testbed to the specifications of another country or region, as well as expanding the testbed’s grid integrated parameters, thereby allowing it to be implemented for broader RE integration study purposes.

Figures

Non-synchronous Category C RPP voltage and frequency ride-through requirements

Figure 1.

Non-synchronous Category C RPP voltage and frequency ride-through requirements

Testbed design layout

Figure 2.

Testbed design layout

MATLAB SM target sub-system flow diagram and circuit

Figure 3.

MATLAB SM target sub-system flow diagram and circuit

Testbed measurement scopes

Figure 4.

Testbed measurement scopes

Case study 1, normal voltage fluctuations real-time results

Figure 5.

Case study 1, normal voltage fluctuations real-time results

Case study 2, fault ride-through real-time results

Figure 6.

Case study 2, fault ride-through real-time results

Case study 3, trip region real-time results

Figure 7.

Case study 3, trip region real-time results

Case study 4, normal frequency fluctuations real-time results

Figure 8.

Case study 4, normal frequency fluctuations real-time results

Case study 5, fault ride-through real-time results

Figure 9.

Case study 5, fault ride-through real-time results

Case study 6, trip region real-time results

Figure 10.

Case study 6, trip region real-time results

Case study 7, RPP active and reactive current response

Figure 11.

Case study 7, RPP active and reactive current response

References

Abdalla, O.H. (2016), “Technical design specifications and criteria for integrating PV systems into distribution networks in Egypt”, The 1st FUE International Conference on New Energy and Environmental Engineering, pp. 11-114, available at: www.academia.edu/download/57962711/Solar_ssPV_Paper.pdf

Abdalla, O.H. (2018), “Technical requirements for connecting medium and large solar power plants to electricity networks in Egypt”, Journal of the Egyptian Society of Engineers, Vol. 57 No. 1, pp. 25-36, available at: www.works.bepress.com/omar/103/

Al Riyami, H.A.S., Al Busaidi, A.G., Al Nadabi, A., Abdalla, O.H., Al Manthari, K., Hagenkort, B. and Fahmi, R. (2016), “Grid code compliance for integrating 50 MW wind farm into Dhofar power grid”, The Selected Works of Omar H. Abdalla, pp. 152-161, available at: www.works.bepress.com/omar/25/

All about circuits (2020), “True, reactive, and apparent power, alternating current (AC)”, available at: www.allaboutcircuits.com/textbook/alternating-current/chpt-11/true-reactive-and-apparent-power/ (accessed 20 February 2023).

Ayamolowo, O.J., Manditereza, P.T. and Kusakana, K. (2022), “South Africa power reforms: the path to a dominant renewable energy-sourced grid’, in Elsevier ltd”, Energy Reports, Vol. 8, pp. 1208-1215, doi: 10.1016/j.egyr.2021.11.100.

Buraimoh, E. and Davidson, I.E. (2022), “Laboratory procedure for real-time simulation experiment of renewable energy systems on OPAL-RT digital simulator”, 10th IEEE International Conference on Smart Grid (icSmartGrid), pp. 221-226, doi: 10.1109/icSmartGrid55722.2022.9848570.

Centre for Innovation in Mathematics Teaching (2009), “Complex numbers”, available at: www.cimt.org.uk/projects/mepres/alevel/fpure_ch3.pdf (accessed 20 February 2023).

Creamer, T. (2021), “With South Africa set to breach 1000 hours of load-shedding in 2021, concerns grow about finance for maintenance”, available at: www.engineeringnews.co.za/article/with-south-africa-set-to-breach-1-000-hours-of-load-shedding-in-2021-concerns-grow-about-maintenance-finance-and-skills-2021-11-08 (accessed 1 June 2022).

Department of Energy (2010), “Executive summary of the draft integrated electricity resource plan for South Africa – 2010 to 2030”, available at: www.doe-irp.co.za/content/Executive_Summary_Draft_IRP2010_22Oct2010.pdf (accessed 15 May 2022).

Department of Energy (2011), “Integrated resource plan for electricity 2010-2030”, available at: www.energy.gov.za/IRP/irpfiles/IRP2010_2030_Final_Report_20110325.pdf (accessed 15 May 2022).

Department of Energy and Department of National Treasury (2018), “Independent power producers procurement programme (IPPPP), an overview”, available at: www.ipp-projects.co.za/Publications/GetPublicationFile?fileid=71515393-218c-e811-9489-2c59e59ac9cd&fileName=20 (accessed 22 May 2022).

Department of Mineral Resources and Energy (2019), “Integrated resource plan 2019”, Department of Energy Republic of South Africa, available at: www.energy.gov.za/IRP/2019/IRP-2019.pdf (accessed 1 June 2022).

Dixon, A. (2019), “Per unit systems for frequency analysis”, Modern Aspects of Power System Frequency Stability and Control, pp. 39-45, doi: 10.1016/b978-0-12-816139-5.00003-5.

Eskom (2019a), “Camden power station”, Eskom heritage, available at: www.eskom.co.za/sites/heritage/Pages/Camden.aspx (accessed 3 June 2022).

Eskom (2019b), “Eskom transmission development plan 2019 to 2028 (public version)”, Transmission Development Plan, available at: www.eskom.co.za/Whatweredoing/TransmissionDevelopmentPlan/Documents/2019-2028PublicTDPreport1.pdf (accessed 3 June 2022).

Eskom generation division (2019), “Eskom power stations”, doi: 10.1016/s1359-6128(05)00671-3.

Etxegarai, A., Eguia, P., Torres, E., Buigues, G. and Iturregi, A. (2017), “Current procedures and practices on grid code compliance verification of renewable power generation”, Renewable and Sustainable Energy Reviews, Vol. 71, pp. 191-202, doi: 10.1016/j.rser.2016.12.051.

Goswami, R. and Das, R. (2020), “Waste heat recovery from a biomass heat engine for thermoelectric power generation using two-phase thermosyphons”, Renewable Energy, Elsevier Ltd, Vol. 148, pp. 1280-1291, doi: 10.1016/j.renene.2019.10.067.

Hanto, J., Schroth, A., Krawielicki, L., Oei, P.Y. and Burton, J. (2022), “South Africa’s energy transition – unraveling its political economy”, Energy for Sustainable Development, Elsevier Inc, Vol. 69, pp. 164-178, doi: 10.1016/j.esd.2022.06.006.

Kemal, M.S., Pedersen, R., Iov, F. and Olsen, R. (2017), “DiSC-OPAL: a simulation framework for real-time assessment of distribution grids”, 2017 Workshop on Modeling and Simulation of Cyber-Physical Energy Systems (MSCPES), pp. 1-5. doi: 10.1109/mscpes.2017.8064541.

Li, Y., Shi, W., Zhang, X. and Tang, W. (2016), “Renewable energy generation power in the loop real time simulation platform research”, 2016 19th International Conference on Electrical Machines and Systems (ICEMS), The Institute of Electrical Engineers of Japan, available at: www.ieeexplore-ieee-org.libproxy.cput.ac.za/document/7837493

Li, Y., Tian, X., Liu, C., Su, Y., Li, L., Zhang, L., Sun, Y. and Li, J. (2017), “Study on voltage control in distribution network with renewable energy integration”, 2017 IEEE Conference on Energy Internet and Energy System Integration (EI2), doi: 10.1109/EI2.2017.8245755.

Mararakanye, N. and Bekker, B. (2019), “Renewable energy integration impacts within the context of generator type, penetration level and grid characteristics”, Renewable and Sustainable Energy Reviews, Elsevier Ltd., Vol. 108, pp. 441-451, doi: 10.1016/j.rser.2019.03.045.

Maune, B. (2019), “Load shedding: timeline of Eskom’s battle to keep the lights on”, The South African, 10 December 2019, available at: www.thesouthafrican.com/news/eskom-load-shedding-timeline-since-2007/ (accessed 5 May 2022).

Mchunu, T, NERSA and Eskom Transmission Division (2019), “Grid connection code for renewable power plants connected to the electricity transmission system or the distribution system in South Africa”, Eskom Transmission Division, available at: www.nersa.org.za/Admin/Document/Editor/file/Electricity/TechnicalStandards/NewableEnergy/GridConnectionCodeforRenewablePowerPlantsRPPs.pdf%0A

Merabet, A., Tawfique, K.A., Islam, M.A., Enebeli, S. and Beguenane, R. (2014), “Wind turbine emulator using OPAL-RT real-time HIL/RCP laboratory”, Proceedings of the International Conference on Microelectronics, ICM, IEEE, March 2015, pp. 192-195, doi: 10.1109/ICM.2014.7071839.

Müller, I., Blauth, Y.B., Pereira, C.E., Gabiatti, G. and Bonan, G. (2009), “Simplifying the design of parallel-connected UPS inverters by means of the per unit system”, 2009 Brazilian Power Electronics Conference, COBEP2009, pp. 544-550, doi: 10.1109/COBEP.2009.5347691.

National Energy Regulator of South Africa (2016), “Grid connection code for renewable power plants (RPPs) connected to the electrical transmission system or the distribution system in South Africa”, doi: 10.1109/powerafrica.2017.7991291.

National Energy Regulator of South Africa (NERSA) (2012), “Electricity supply statistics for South Africa 2012”, doi: 10.1787/data-00462-en.

Noureen, S.S., Roy, V. and Bayne, S.B. (2018), “An overall study of a real-time simulator and application of RT-LAB using MATLAB SimPowerSystems”, 2017 IEEE Green Energy and Smart Systems Conference, IGESSC 2017, pp. 1-5, doi: 10.1109/IGESC.2017.8283453.

Paquin, J.N., Moyen, J., Dumur, G. and Lapointe, V. (2007), “Real-time and off-line simulation of a detailed wind farm model connected to a multi-bus network”, 2007 IEEE Canada Electrical Power Conference, EPC 2007, pp. 145-152, doi: 10.1109/EPC.2007.4520321.

Roy, S., Das, R. and Saha, U.K. (2018), “An inverse method for optimization of geometric parameters of a Savonius-style wind turbine”, Energy Conversion and Management, Elsevier, Vol. 155, pp. 116-127, doi: 10.1016/j.enconman.2017.10.088.

RSA Department of Energy (2018), ‘2018 South African Energy Sector Report’, in Ratshomo, K. (Ed.), Department of Energy, doi: 10.2172/1430268 (accessed 20 May 2022).

Saha, S., Saleem, M. and Roy, T. (2023), “Impact of high penetration of renewable energy sources on grid frequency behaviour”, International Journal of Electrical Power and Energy Systems, Vol. 145, doi: 10.1016/j.ijepes.2022.108701.

Singh, S.K., Padhy, B.P., Chakrabarti, S., Singh, S.N., Kolwalkar, A. and Kelapure, S.M. (2015), “Development of dynamic test cases in OPAL-RT real-time power system simulator”, 2014 18th National Power Systems Conference, NPSC 2014, IEEE, doi: 10.1109/NPSC.2014.7103856.

Statistics South Africa (2018), “Electricity generated and available for distribution (preliminary)”, available at: www.statssa.gov.za/publications/P4141/P4141October2018.pdf (accessed 5 June 2022).

Statistics South Africa (2019), “Electricity generated and available for distribution”, available at: www.statssa.gov.za/publications/P4141/P4141November2019.pdf (accessed 5 June 2022).

Statistics South Africa (2020), “Key findings: P4141 – Electricity generated and available for distribution, November 2019”, available at: www.statssa.gov.za/?page_id=1856&PPN=P4141&SCH=7534 (accessed 5 June 2022).

Yan, R., Saha, T.K., Modi, N., Masood, N.A. and Mosadeghy, M. (2015), “The combined effects of high penetration of wind and PV on power system frequency response”, Applied Energy, Elsevier, Vol. 145, pp. 320-330, doi: 10.1016/j.apenergy.2015.02.044.

Acknowledgements

Funding: Article Processing Charges for the publishing of this article was covered by the Cape Peninsula University of Technology, CPUT - SANLiC agreement with Emerald.

Corresponding author

Gideon Daniel Joubert can be contacted at: gd.joubert2@gmail.com

Related articles