Power-Flow Simulations for Integrating Renewable Distributed Generation from Biogas, Photovoltaic, and Small Wind Sources on an Underground Distribution Feeder

: The rapid expansion of distributed generation leads to the integration of an increasing number of energy generation sources. However, integrating these sources into electrical distribution networks presents speciﬁc challenges to ensure that the distribution networks can effectively accommodate the associated distributed energy and power. Thus, it is crucial to evaluate the electrical effects of power along the conductors, components, and loads. Power-ﬂow analysis is a well-established numerical methodology for assessing parameters and quantities within power systems during steady-state operation. The University of S ã o Paulo’s Cidade Universit á ria “Armando de Salles Oliveira” (CUASO) campus in S ã o Paulo, Brazil, features an underground power distribution system. The Institute of Energy and Environment (IEE) leads the integration of several distributed generation (DG) sources, including a biogas plant, photovoltaic installations, and a small wind turbine, into one of the CUASO’s feeders, referred to as “USP-105”. Load-ﬂow simulations were conducted using the PowerWorld TM Simulator v.23, considering the interconnection of these sources. This dataset provides comprehensive information and computational ﬁles utilized in the simulations. It serves as a valuable resource for reanalysis, didactic purposes, and the dissemination of technical insights related to DG implementation.


Summary
The rapid expansion of distributed generation (DG) throughout the 21st century is one of the many outcomes of the profound transformations that the organization of the electrical industry has been undergoing in nearly all countries.
Following a century dominated by vertically integrated generation, transmission, and distribution systems constituting regulated natural monopolies, structural and regulatory changes were initiated in the last few decades.These changes were driven by advances in information technology, online data management, and the automation and control of the constituent units comprising the electrical systems.Political decisions adopted by various countries or regions aimed at market liberation and the introduction of competitive mechanisms for the generation and commercialization of electric energy.
A new cycle of structural changes has emerged in the current century, propelled by a global movement to respond effectively to a set of concerns and opportunities, including environmental protection on a global scale related to climate change, the reduction in greenhouse gas emissions, and decreased dependency on fossil fuels, coupled with the Data 2023, 8, 150 2 of 12 expansion of energy production and consumption in a manner that is attentive to environmental issues.In addition, the potential for maintaining or even improving access to energy services such as lighting, thermal comfort, motor force, and mobility, among others, can be achieved through the more efficient consumption of energy as a consequence of the introduction of technologies, processes, and demand management systems.Demand Side Management (DSM), for instance, plays an important role, allowing for a reduction in the consumption intensity of primary resources in generation.Technological advances, particularly wind and photovoltaics, have also played a key role when substituting expensive fossil-derived energy.Furthermore, investments in electricity transmission and distribution systems could be postponed or even avoided.
Smaller distributed wind or rooftop photovoltaic installations, combined or not with storage systems, are characterized as Distributed Energy Resources (DER).New challenges arise for distribution system planners as more DER are interconnected to the distribution network, which must handle variable power and operate differently either to supply power to normal consumers or to host excess power injected by DER.Despite distributed sources possibly improving end-user power quality and reliability [1], the concepts of reliability and the learned operation methodology must be fulfilled while respecting the technical requirements and limits, otherwise the photovoltaic DG power fed into the network may lead to technical issues such as overvoltage [2].
When considering electrical energy generation, transmission, and distribution, it is imperative to evaluate the electrical effects of the power along conductors, components, and loads.Load-flow (or power-flow) is a well-established numerical methodology for analyzing parameters and quantities within power systems during steady-state operation.Primarily, power-flow determines voltages, currents, angles, and active and reactive powers.These calculations can be executed using various tools [3], and in [4] a benchmark comparing performance between various computational tools [5][6][7][8][9] for analyzing power grids is presented.
A power-flow study usually employs moderately simplified modelling of the system s components derived from a single-line diagram and uses the per-unit (p.u.) system to determine the operating point.This approach considers the balance of load and generation, mainly when dealing with multiple generators, as in DG conditions.
The University of São Paulo's campus, Cidade Universitária "Armando de Salles Oliveira" (CUASO), located in the city of São Paulo, Brazil, features an underground power distribution system.The Institute of Energy and Environment (IEE) is at the forefront of a campus project aimed at integrating multiple DG sources on the campus, including a biogas plant, photovoltaic (both rooftop and ground) installations, and a small wind turbine.These DG installations are predominantly concentrated on a particular feeder referred to as "USP-105", one of the five circuits constituting the CUASO distribution network.
A hosting capacity study, which outlines a methodology to estimate the potential accommodation of photovoltaic installations within the USP-105 circuit, is detailed in [10].Figure 1 illustrates a map of CUASO, depicting the main feeder medium voltage path of USP-105, its associated laterals, and the position of twelve load centers comprising educational and research facilities served by the USP-105 circuit.These load centers include the Institute of Energy and Environment (IEE, São Paulo, Brazil), PUSP-C (Prefecture of the USP Campus, São Paulo, Brazil), NUCEL (Cellular and Molecular Therapy Center, São Paulo, Brazil), HU (University Hospital, São Paulo, Brazil), FOFITO (Department of Physiotherapy, Speech Therapy and Occupational Therapy, São Paulo, Brazil), ICB "I", "II", "III" and "IV" (Institute of Biomedical Sciences, São Paulo, Brazil), FO (School of Dentistry, São Paulo, Brazil), FMVZ (School of Veterinary Medicine and Animal Science, São Paulo, Brazil), and IB (Institute of Biosciences, São Paulo, Brazil).This dataset comprises essential information, intricate details, and source files used to conduct power-flow studies of the USP-105 underground feeder within the distribution medium voltage network at the University of São Paulo campus, considering the integration of existing and foreseen DG sources.The power-flow simulations were executed using the PowerWorld TM Simulator, v.23 [11][12][13][14][15].This dataset provides valuable insights for reanalysis, didactic applications, and the dissemination of technical knowledge regarding the practical implementation of DG.

Data Description
The distribution system´s topology, including cables, transformers, loads, and generators, was configured within the PowerWorld TM simulator for processing the powerflow calculations.The computational files associated with the USP-105 circuit were then included in the dataset, titled as follows:  This dataset comprises essential information, intricate details, and source files used to conduct power-flow studies of the USP-105 underground feeder within the distribution medium voltage network at the University of São Paulo campus, considering the integration of existing and foreseen DG sources.The power-flow simulations were executed using the PowerWorld TM Simulator, v.23 [11][12][13][14][15].This dataset provides valuable insights for reanalysis, didactic applications, and the dissemination of technical knowledge regarding the practical implementation of DG.

Data Description
The distribution system s topology, including cables, transformers, loads, and generators, was configured within the PowerWorld TM simulator for processing the power-flow calculations.The computational files associated with the USP-105 circuit were then included in the dataset, titled as follows:

Modelling
The current section presents information regarding the modelling process, including the medium-voltage circuit, loads, DG sources, and the generators.

Modelling
The current section presents information regarding the modelling process, including the medium-voltage circuit, loads, DG sources, and the generators.The USP-105 is a radial circuit with a rated voltage of 13.8 kV and a length of 3.5 km.The cross-sectional area of the main feeder cable is 3 × 240 mm 2 , and its laterals use either 3 × 70 mm 2 or 3 × 35 mm 2 cables.The main insulation is composed of a polyethylene compound (XLPE) capable of operating at a temperature of 90 °C.
The medium-voltage cable impedances are provided in Table 1, and the ampacity values were reevaluated according to the installation criteria outlined in the IEC 60502-2 Standard (for tripolar cables in buried earthenware ducts).The USP-105 is a radial circuit with a rated voltage of 13.8 kV and a length of 3.5 km.The cross-sectional area of the main feeder cable is 3 × 240 mm 2 , and its laterals use either 3 × 70 mm 2 or 3 × 35 mm 2 cables.The main insulation is composed of a polyethylene compound (XLPE) capable of operating at a temperature of 90 • C.
The medium-voltage cable impedances are provided in Table 1, and the ampacity values were reevaluated according to the installation criteria outlined in the IEC 60502-2 Standard [16] (for tripolar cables in buried earthenware ducts).Based on the distances of the sections depicted in the single-line diagram (Figure 2) and the impedance values (Table 1), the impedances of each section were calculated.For p.u. calculations, a base power of 1 MVA and a base voltage of 13.8 kV were adopted, leading to the base impedance of 190.44 Ω.Consequently, the calculated modelled impedances (both resistance and reactance) are presented in Table 2.The distribution transformers were modelled using their nominal ratios (13.8 kV/220 V) and the series resistances and reactances according to the rated power, as presented in Table 3.The twelve loads positioned along the lines were modelled using their active and reactive profile curves, sampled every 15 min across a 24 h day.These profiles were obtained for high load and light load conditions.The high load condition refers to weekdays in the summer when the campus experiences peak energy and power consumption.Conversely, the light load condition refers to weekend days in the winter when energy demand is lower.
The electric distribution system operation was performed by the University of São Paulo at the CUASO.Load measurement data were acquired through its supervisory system from March 2020 to August 2021.The heavy and light load months were identified as March (summer, during classes) and July (winter, no-classes period), respectively.The load profiles employed in the power-flow simulations were derived from linearizing the averaged curves obtained from each set of active and reactive power measurements, considering heavy or light load conditions.Figure 3 illustrates the modelling process outcomes for two distinct sets of load profiles derived from the acquired curves: the high load active power of the HU consumer and the light load reactive power of the IB consumer.The electric distribution system operation was performed by the University of São Paulo at the CUASO.Load measurement data were acquired through its supervisory system from March 2020 to August 2021.The heavy and light load months were identified as March (summer, during classes) and July (winter, no-classes period), respectively.The load profiles employed in the power-flow simulations were derived from linearizing the averaged curves obtained from each set of active and reactive power measurements, considering heavy or light load conditions.Figure 3 illustrates the modelling process outcomes for two distinct sets of load profiles derived from the acquired curves: the high load active power of the HU consumer and the light load reactive power of the IB consumer.The 15 min hourly modelled curves for all the loads are presented in the dataset.These curves are used in the power-flow simulations and can be identified by the following file names: During the PowerWorld TM simulations, the load profiles were replicated using the time step simulation tool, specifying the active and reactive powers at 15 min intervals throughout the 24 h of the day.The loads' variations were incorporated in the time step simulation files named in the dataset, as follows: Figure 4 shows examples of the modelled active and reactive power profiles for the HU consumer under both heavy and light load conditions.The 15 min hourly modelled curves for all the loads are presented in the dataset.These curves are used in the power-flow simulations and can be identified by the following file names: During the PowerWorld TM simulations, the load profiles were replicated using the time step simulation tool, specifying the active and reactive powers at 15 min intervals throughout the 24 h of the day.The loads' variations were incorporated in the time step simulation files named in the dataset, as follows: Figure 4 shows examples of the modelled active and reactive power profiles for the HU consumer under both heavy and light load conditions.

DG Sources and Generator Modelling
The distribution system was supplied by the public power network at the connection substation and named ETD-USP at the 88 kV level, which was modelled as a slack bus/generator used as a reference bus 13.8 kV/1 p.u.
Fifteen-minute hourly power profiles were modelled for ten existing and planned DG sources (G1 to G10 in Figure 2).These profiles, illustrated in Figure 5, consider three distinct types of generation sources:
The 15 min hourly modelled DG profiles for all the sources are in the dataset placed within the file named "DG Table 24 h.csv".
Similar to the load modelling, the DG generation curves were simulated using the time step simulation tool in PowerWorld TM , in which the active powers were specified at 15 min intervals throughout the 24 h of the day.The generation profiles were incorporated in the same time step simulation files described in Section 2.1.2,whose file names in the dataset are as follows:

DG Sources and Generator Modelling
The distribution system was supplied by the public power network at the connection substation and named ETD-USP at the 88 kV level, which was modelled as a slack bus/generator used as a reference bus 13.8 kV/1 p.u.
Fifteen-minute hourly power profiles were modelled for ten existing and planned DG sources (G1 to G10 in Figure 2).These profiles, illustrated in Figure 5, consider three distinct types of generation sources:
The 15 min hourly modelled DG profiles for all the sources are in the dataset placed within the file named "DG Table 24 h.csv".
Similar to the load modelling, the DG generation curves were simulated using the time step simulation tool in PowerWorld TM , in which the active powers were specified at 15 min intervals throughout the 24 h of the day.The generation profiles were incorporated in the same time step simulation files described in Section 2.1.2,whose file names in the dataset are as follows:

Results
Figure 6 provides, for illustrative purposes, screen captures from PowerWorld TM demonstrating the power-flow balance under two distinct conditions.

Results
Figure 6 provides, for illustrative purposes, screen captures from PowerWorld TM demonstrating the power-flow balance under two distinct conditions.
Since one can use power-flow simulations for countless purposes, the current dataset offers a collection of organized results, which can be found in the following dataset files:
Figure 7 illustrates some of the analyses derived from the simulation results stored in the aforementioned files, including the voltage profile along the feeder, the calculation of the power losses, and the determination of the maximum loading of the feeder.Since one can use power-flow simulations for countless purposes, the current dataset offers a collection of organized results, which can be found in the following dataset files: • Results Light Load.csv; • Results Heavy Load.csv.
Figure 7 illustrates some of the analyses derived from the simulation results stored in the aforementioned files, including the voltage profile along the feeder, the calculation of the power losses, and the determination of the maximum loading of the feeder.

Materials and Methods
The dataset contains proprietary files from the power system simulation software PowerWorld TM , version 23.The extension of these files are: .pwb,.pwdand .tsb.
The modelling process for representing the electrical loads along the medium-voltage distribution system relied on measurements of active and reactive powers obtained from a SCADA system implemented by the University of São Paulo at CUASO.
The supervisory system was built upon the Elipse TM SCADA platform [18], a widely used process management tool that allows the integration of energy meters or other measurement equipment that uses standardized communication protocols.
The preferred means of communication used for the traffic of information in CUASO are the Internet, through the MODBUS TCP/IP communication protocol, and a private fiber-optical network using the MODBUS RTU serial protocol.
Figure 8 showcases a few screens of the supervisory system in operation.

Materials and Methods
The dataset contains proprietary files from the power system simulation software PowerWorld TM , version 23.The extension of these files are: .pwb,.pwdand .tsb.
The modelling process for representing the electrical loads along the medium-voltage distribution system relied on measurements of active and reactive powers obtained from a SCADA system implemented by the University of São Paulo at CUASO.
The supervisory system was built upon the Elipse TM SCADA platform [19], a widely used process management tool that allows the integration of energy meters or other measurement equipment that uses standardized communication protocols.
The preferred means of communication used for the traffic of information in CUASO are the Internet, through the MODBUS TCP/IP communication protocol, and a private fiber-optical network using the MODBUS RTU serial protocol.
Figure 8 showcases a few screens of the supervisory system in operation.
The information was acquired using energy meters and/or protection relays installed within the power supply cabinets and was transmitted to the system's server.Physical measurements could be recorded in the supervisory system at a maximum rate of 1 register per second.However, for this study, the measurements were evaluated and modelled on a 15 min interval basis.
The CSV files provided in the dataset employ a semicolon symbol (;) as the column separator, while a dot symbol (.) is used as the decimal separator.The first row of each CSV file corresponds to the header row, helping in the identification of the data.The information was acquired using energy meters and/or protection relays installed within the power supply cabinets and was transmitted to the system's server.Physical measurements could be recorded in the supervisory system at a maximum rate of 1 register per second.However, for this study, the measurements were evaluated and modelled on a 15 min interval basis.
The CSV files provided in the dataset employ a semicolon symbol (;) as the column separator, while a dot symbol (.) is used as the decimal separator.The first row of each CSV file corresponds to the header row, helping in the identification of the data.

2. 1 . 1 .
Medium-Voltage Modelling The USP-105 single-line diagram is depicted in Figure 2. The simulation model includes the medium-voltage buses and cables and the low-voltage buses and their corresponding transformers once the loads and DG installations have become interconnected within the low-voltage distribution network.Data 2023, 8, x FOR PEER REVIEW 4 of 12

2. 1 . 1 .
Medium-Voltage ModellingThe USP-105 single-line diagram is depicted in Figure2.The simulation model includes the medium-voltage buses and cables and the low-voltage buses and their corresponding transformers once the loads and DG installations have become interconnected within the low-voltage distribution network.

Figure 3 .
Figure 3. Examples of modelled load curves from measured data.

Figure 3 .
Figure 3. Examples of modelled load curves from measured data.

Figure 6 .
Figure 6.Screen captures of the power-flow calculation software PowerWorld TM .(a) Light load condition, no DG sources.(b) Heavy load condition, all DG sources interconnected.Figure 6. Screen captures of the power-flow calculation software PowerWorld TM .(a) Light load condition, no DG sources.(b) Heavy load condition, all DG sources interconnected.

Figure 6 .
Figure 6.Screen captures of the power-flow calculation software PowerWorld TM .(a) Light load condition, no DG sources.(b) Heavy load condition, all DG sources interconnected.Figure 6. Screen captures of the power-flow calculation software PowerWorld TM .(a) Light load condition, no DG sources.(b) Heavy load condition, all DG sources interconnected.

Figure 8 .
Figure 8. Screen captures of the SCADA implemented for CUASO distribution system operation and supervision.

Figure 8 .
Figure 8. Screen captures of the SCADA implemented for CUASO distribution system operation and supervision.

Table 2 .
Modelled impedances of the medium-voltage sections.