ORIGINAL RESEARCH article

Front. Energy Res., 30 August 2024

Sec. Smart Grids

Volume 12 - 2024 | https://doi.org/10.3389/fenrg.2024.1399114

AC/DC optimal power flow and techno-economic assessment for hybrid microgrids: TIGON CEDER demonstrator

  • 1. Centro Tecnológico CARTIF, Boecillo, Spain

  • 2. Centro de Desarrollo de Energías Renovables (CEDER)—Centro de Investigaciones Energéticas Medioambientales y Técnológicas (CIEMAT), Soria, Spain

  • 3. Departamento de Ingeniería Eléctrica, Escuela de Ingenierías Industriales, Universidad de Valladolid, Valladolid, Spain

  • 4. Instituto de las Tecnologías Avanzadas de la Producción, Universidad de Valladolid, Valladolid, Spain

Abstract

In recent years, the interest in electric direct current (DC) technologies (such as converters, batteries, and electric vehicles) has increased due to their potential in energy efficiency and sustainability. However, the vast majority of electric systems and networks are based on alternating current (AC) as they also have certain advantages regarding cost-effective transport and robustness. In this paper, an AC/DC optimal power flow method for hybrid microgrids and several key performance indicators (KPIs) for its techno-economic assessment are presented. The combination of both calculations allows users to determine the viability of their hybrid microgrids. AC/DC networks have been modeled considering their most common elements. For the power flow method, polynomial optimization is formulated considering four different objective functions: the minimization of energy losses, voltage deviation, and operational costs and the maximization of the microgrid generation. The power flow method and the techno–economic analysis are implemented in Python and validated in the Centro de Desarrollo de Energías Renovables (CEDER) demonstrator for TIGON. The results show that the calculated power flow variables and those measured at CEDER are practically the same. In addition, the KPIs are obtained and compared for four operating scenarios: baseline, no battery, battery flexibility, and virtual battery (VB) flexibility. The last scenario results in the most profitable option.

1 Introduction

The global shift toward decarbonization has propelled significant transformations in the design, operation, and management of electric grids. The urgent need to mitigate climate change has led to the adoption of renewable energy sources and the phasing out of fossil fuel-based power generation, resulting in a paradigm shift in the electricity sector. The integration of variable renewable sources, such as solar and wind, presents unique challenges due to their intermittent nature and geographical distribution. As a result, electric grids have witnessed a remarkable transition toward more dynamic, flexible, and intelligent systems.

Microgrids (; ; ; ; ) are localized and self-contained electricity distribution systems. They have gained prominence due to their ability to effectively integrate distributed energy resources (DERs). DERs include small-scale renewable energy installations, energy storage systems, and demand response capabilities. Microgrids provide an innovative solution to enhance the resilience, reliability, and sustainability of the electric grid at a smaller scale while offering opportunities for local energy generation, utilization, and management. In general, microgrids can work as an island or connected to the main power network, which acts as an external grid ().

Furthermore, the choice between alternating current (AC) and direct current (DC) elements within microgrids has become a subject of considerable interest (). While AC has historically been the dominant standard for power transmission and distribution, recent advancements in DC technologies, such as state-of-the-art batteries and electric vehicles (EVs), have brought attention to these systems (; ; ).

In microgrids, where local generation and consumption are tightly integrated, DC elements offer several benefits (; ; ; ). First, DC distribution systems enable higher efficiency in the utilization of renewable energy sources. Most renewable energy technologies, such as solar panels and batteries, inherently generate and store DC power. By directly integrating these DC sources into the microgrid without the need for AC/DC conversions, energy losses associated with multiple conversions can be minimized, resulting in improved overall system efficiency. Moreover, DC systems offer increased flexibility for the integration of emerging technologies. As the demand for EVs grows, the DC charging infrastructure becomes crucial. DC microgrids can seamlessly accommodate EV charging stations without the need for additional power conversion equipment, reducing infrastructure costs and improving charging efficiency ().

Despite that, AC technologies still offer certain advantages that make them relevant and preferred in specific aspects of microgrid design and operation. For instance, high-voltage AC transmission systems inherently offer better voltage control than DC through reactive power, they are easier to isolate and interrupt in the case of faults, and the existing infrastructure is more abundant. In the end, both AC and DC technologies have their unique advantages and limitations, and their selection should be based on the careful evaluation and analysis of the specific circumstances. In this context, hybrid AC/DC microgrids emerge as a suitable solution for the transition to an electricity system with reduced or zero greenhouse gas emissions, taking advantage of the benefits of both forms of electricity current. One example of a project that seeks to maximize the benefits of these networks is the Horizon 2020 European project called Towards Intelligent DC-based hybrid Grids Optimizing the Network performance (TIGON). It also aims to improve the reliability, resilience, performance, and cost efficiency of hybrid AC/DC grids.

The main contribution of this article is the development of a procedure to study and evaluate the correct operation and the technical and economic feasibility of hybrid microgrid installations. The developed procedure consists of two components. The first component is a power flow calculation method for hybrid AC/DC microgrids based on optimization. The power flow can be performed by choosing among four different cost functions, depending on the objective to be achieved. The method is based on the recent literature on both AC and DC load flows (, ; ; ) but differs from them by including a larger number of elements and a different formulation. The second component is a techno-economic evaluation based on key performance indicators (KPIs). In this case, several references were used (; ; ; ) although the methodology of this article differs as it is specifically adapted to hybrid AC/DC microgrids. Another contribution of the article is the validation of the developed procedure in a real hybrid microgrid located in the facilities of the Centro de Desarrollo de Energías Renovables (CEDER), which is part of the EU-funded TIGON project dedicated to the demonstration of innovations in hybrid microgrids for greener, more resilient, and safer power grids. The measurable variables in the CEDER hybrid AC/DC microgrid have been compared with the values obtained in the power flow simulation for validation, demonstrating the accuracy and validity of the developed procedure.

The remainder of this article is organized as follows: Section 2 presents the optimal power flow formulation and the possible objective functions that can be used for optimization. Section 3 presents the KPIs used for the techno-economic assessment. Section 4 details the characteristics of the TIGON CEDER demonstrator. Section 5 presents the experimental and simulation results for different operating scenarios. The conclusion is given in Section 6. Finally, the microgrid model used for the AC/DC optimal power flow and the techno-economic analysis is detailed in Appendix A.

2 AC/DC optimal power flow

The AC/DC optimal power flow allows us to study the feasibility of the microgrid operation, self-consumption capability, load supply, and power losses.

2.1 Formulation

The nomenclature used in the formulation of the problem is detailed at the end of the paper. All electrical variables are represented in phasor form. Consider an electrical network whose topology is represented by a graph , where is the set of buses (vertices) and is the set of lines (edges). The lines are unordered pairs of buses , where and are the pair of buses connected by the line.

A bus is adjacent to another bus if there is a line connecting them, i.e., if . The set of buses adjacent to bus is denoted by and is defined as follows:

The state of the network variables is physically related (). First, the apparent power at each bus is expressed using Equation 2, where is the complex conjugate of .

is an aggregation of currents, as stated in Equation 3.

Both and can be calculated using Ohm’s law, as shown in Equations 4, 5, respectively.

Power flow equations are obtained combining all the previous expressions. In the case of AC buses, the resulting expressions are Equations 67 (; ).

, which is the power injected or transferred from converters, is added to the active power load flow equation. In the case of DC buses, the expressions are Equations 8, 9.

In any case, and are the sum of generation and demand (Equations 10, 11).

At all times, and are limited to the nominal power and the minimum power of the generators (Equations 12, 13, respectively).

If two buses are connected through a converter, the expression that describes the power exchange between them is given by Equation 14.

The maximum and minimum voltage limits of buses are maintained, as shown in Equation 15.

When grid-forming mode is activated in a converter, the voltage at the output bus of the converter is set to its nominal value, as described in Equation 16.

When the converter is in grid-following mode, this restriction is not considered.

The total current that lines can transport is limited by Equation 17.

Transformers cannot exceed their nominal power when operating, as shown in Equations 18, 19.

In the case of converters, this limitation is expressed using Equation 20.

2.2 Optimization

Four different objective functions

can be chosen for minimization in the AC/DC optimal power flow, depending on the objective to be achieved.

  • 1. H1: Total active power generated. This function focuses on reducing energy losses (Equation 21).

  • 2. H2: Bus voltage deviation from their nominal value. This function focuses on achieving grid stability (Equation 22).

  • 3. H3: Total amount of operational costs associated with each generator. This function focuses on achieving economic savings (Equation 23).

  • 4. H4: Active power microgrid generation. This function focuses on making the highest possible use of the available generation resources in the microgrid (Equation 24). The objective function is set to be negative in order to calculate maximization.

The variables that are optimized when performing the power flow are , , , , and .

3 Techno-economic assessment

The techno-economic assessment presented in this paper consists of the calculation of eight KPIs. They allow the evaluation of a microgrid in terms of costs, energy generation, storage capabilities, and financial feasibility. Moreover, it provides critical insights for making informed decisions and maximizing the overall performance of the microgrid. The KPIs, presented below, are divided into two categories: technical and economic.

3.1 Technical KPIs

  • KPI1: Electrical energy generated. This indicator represents the amount of electrical energy generated per year.

  • KPI2: CO2emissions. This indicator represents the total CO2 emitted by all energy carriers associated with the primary energy use in the microgrid per year.

  • KPI3: Self-consumption percentage. This indicator represents the amount of energy obtained from the generators of the microgrid in relation to the energy used by the loads.

  • KPI4: Storage flexibility. This indicator represents the total flexible power available in the microgrid due to the storage systems.

3.2 Economic KPIs

  • KPI5: Total life cycle income. This indicator represents the total income earned by the microgrid over its useful life. This index includes income from energy production and flexibility management.

  • KPI6: Total life cycle cost. This indicator represents the total cost incurred by the microgrid over its useful life. The index includes the investment, operation, and management costs and residual values of the microgrid and generators.

  • KPI7: Payback. This indicator represents the period of time required to recover the capital investment of the microgrid ().

  • KPI8: Levelized cost of energy (LCOE). This indicator represents the price at which the generated electricity should be sold to break even by the end of the useful life of the microgrid (; ).

4 TIGON CEDER demonstrator

CEDER is the acronym for Center for the Development of Renewable Energies1. It is located in Lubia (Soria, Spain) and serves as a national center for energy research that belongs to the Center for Energy, Environmental, and Technological Research (CIEMAT)2, a public research organization under the Ministry of Science and Innovation. This facility covers an area of 640 ha (13,000 m2 built) and features a smart microgrid (electrical and thermal) operated and managed in real time (see Figure 1, left).

FIGURE 1

The Spanish demonstrator of TIGON is installed at CEDER (see

Figure 1

, right) and consists of the following elements:

  • 1. Transformer station: 15 kVAC–400 VAC.

  • 2. Small wind turbine: A three-bladed, horizontal axis wind turbine with a nominal power of 4.2 kW (Ryse E-5).

  • 3. Photovoltaic (PV) system: Comprising 3 strings with 18 modules (URECO) of 410 W each, amounting to a total capacity of 22.14 kW.

  • 4. NMC batteries: 3 modules, each with 80 cells (50 Ah and 3.6 V per cell).

  • 5. Programmable AC loads: Three programmable AC loads of 2.9 kW each.

  • 6. DC loads: Three adjustable DC loads of 4 kW each.

A schematic diagram of the CEDER demonstrator is shown in Figure 2. It consists of AC loads, wind turbines, transformers, and DC sections in the network. The characteristics of the elements used to test the microgrid are given in Tables 13.

FIGURE 2

TABLE 1

Bus (kV)Type
015.000AC
10.400AC
20.400AC
30.630DC
40.800DC
50.860DC
60.800DC
70.230AC
80.630AC
LineLength (km) (/km) (/km)
0(2,1)0.150.50.35
1(4,6)0.200.5
Line (nF/km) (kA)Type
0(2,1)100100AC
1(4,6)100DC
Transformer (kVA) (%) (%)
0(1,0)2501.51.0
Converter (kVA)Control
0(3,2)200.99Grid following
1(3,4)200.99Grid forming
2(4,5)300.86Grid forming
3(4,7)120.99Grid forming
4(8,4)50.99Grid following
GeneratorBus (kW) (kW) (kW) (kW)
0322.140.500
184.200.010−10
220.000.010−10
370.000.010−10
LoadBus (kW) (kW)
0650
1740
StorageBus (kW)
0525

TIGON CEDER electrical data.

TABLE 2

GeneratorBus (€) (€) (€/year)
0323,000100200
1810,000150600
GeneratorBus (€/kWh) (%) (kgCO2/kWh)
030.00333.50.03500
180.00842.00.00464

Techno-economic information of TIGON CEDER generators.

TABLE 3

(€) (€) (€/year) (%) (years) (€/kWh)
195,50028,9001,4001250.145

Techno–economic information of the TIGON CEDER microgrid.

and have been fixed at 5% above and below the nominal value, respectively. of generators 0 and 1 have been obtained from the study by ). Generators 2 and 3 model the reactive power management of converters 0 and 3, respectively. The CEDER microgrid has one connection to an external grid located on bus 0. The frequency is 50 Hz ( rad/s).

and have been calculated as the sum of the investment cost and the residual value of the microgrid equipment and labor, as given in Table 4.

TABLE 4

EquipmentInvestment cost (€)Residual value (€)
DC PV converter23,0002,500
AC PV converter22,0002,000
Wind turbine converter10,000500
Battery34,000400
Battery converter25,0003,000
AC load converter31,0002,500
Wiring15,5004,000
Cabins18,00018,000
Labor17,0000

Investment costs and residual values of TIGON CEDER equipment and labor.

5 Results and discussion

5.1 AC/DC optimal power flow

The AC/DC optimal power flow simulation has been implemented in Python due to its wide range of available open-source optimization libraries. In this paper, the optimization has been performed using CasADi (), which is a software library equipped with specific tools focused on the modeling, optimization, and control of nonlinear dynamic systems. CasADi is widely used to define both mathematical models and constraints involved and allows us to utilize different solvers in order to optimize the problem. In this study, the solver IPOPT () has been used. IPOPT applies sequential quadratic programming (SQP) to solve constrained nonlinear optimization problems, which is the case in the AC/DC optimal power flow. The calculation time is not significant: it is only a few seconds.

Four scenarios have been tested, one per objective function. In all scenarios, TIGON CEDER storage has been considered to act as a load that consumes electricity at half its nominal power (/2), and the external grid has been limited to power consumption (if required). The results obtained in each scenario are given in Table 5.

TABLE 5

Bus Scenario (kW) (kVA) (kV) (rad)
0H10.000000.0000014.828070.06733
1H10.000000.000000.395420.06733
2H10.00000−0.000740.395420.06733
3H119.754500.000000.620460.00000
4H10.000000.000000.800000.00000
5H1−12.500000.000000.860000.00000
6H1−5.000000.000000.799370.00000
7H1−4.000000.000000.230000.00000
8H14.062890.000000.620460.00000
0H2−0.050890.0000014.999870.59858
1H20.000000.000000.400000.59858
2H20.00000−0.000750.400010.59860
3H219.805800.000000.630000.00000
4H20.000000.000000.800000.00000
5H2−12.500000.000000.860000.00000
6H2−5.000000.000000.799370.00000
7H2−4.000000.000000.230000.00000
8H24.063000.000000.630000.00000
0H30.000000.0000014.951320.06573
1H30.000000.000000.398700.06573
2H30.00000−0.000750.398700.06573
3H320.000000.000000.620460.00000
4H30.000000.000000.800000.00000
5H3−12.500000.000000.860000.00000
6H3−5.000000.000000.799370.00000
7H3−4.000000.000000.230000.00000
8H33.817360.000000.620460.00000
0H4−2.494150.0000014.758670.01449
1H40.000000.000000.393610.01460
2H40.000000.001660.394080.01545
3H422.140000.000000.620460.00000
4H40.000000.000000.800000.00000
5H4−12.500000.000000.860000.00000
6H4−5.000000.000000.799370.00000
7H4−4.000000.000000.230000.00000
8H44.200000.000000.620460.00000

Simulation results for different scenarios.

In all scenarios, all the power demanded by loads and storage is delivered. In none of them exists the demand for reactive power, so all the reactive power in the microgrid is generated in the AC section at bus 2 because of the transformer and AC line reactance.

In scenarios H1 and H3, the exact active power is generated to supply the demand and compensate for losses. The difference between them is that the wind turbine reduces its active power generation in scenario H3 because this technology has higher operational costs than the PV, whereas more power is saved in scenario H1.

In scenario H2, at all buses is the nominal value or very close to it. To achieve it, the generation of active power and the consumption of the external grid are precisely adjusted.

In scenario H4, the PV and wind turbine generate the maximum active power, which is then supplied to the external grid.

All the scenarios have been recreated in the real environment of the TIGON CEDER microgrid. The two quantities accessed for measurement are active power and voltage. The values obtained in the measurements taken for each scenario are shown in Table 6.

TABLE 6

Bus Scenario (kW) (kV)Scenario (kW) (kV)
0H10.00014.790H2−0.05014.990
1H10.0000.399H20.0000.399
2H10.0000.399H20.0000.399
3H119.7500.612H219.8000.611
4H10.0000.799H20.0000.808
5H1−12.5000.859H2−12.5000.862
6H1−5.0000.799H2−5.0000.799
7H1−4.0000.230H2−4.0000.230
8H14.1000.629H24.1020.630
0H30.00014.920H4−2.49014.880
1H30.0000.398H40.0000.399
2H30.0000.398H40.0000.399
3H320.0000.613H422.1400.632
4H30.0000.806H40.0000.805
5H3−12.5000.860H4−12.5000.865
6H3−5.0000.801H4−5.0000.799
7H3−4.0000.230H4−4.0000.230
8H33.8000.628H44.2000.625

CEDER microgrid measurements for different scenarios.

The measures obtained are very similar to the calculated values in the simulations, as shown in Figures 3, 4. In Table 7, the mean relative error (MRE) has been calculated for both active power and voltage using Equation Equation 33.

FIGURE 3

FIGURE 4

TABLE 7

MRE (%)H1H2H3H4
Active power0.100.310.050.02
Voltage0.560.550.440.80

Mean relative error of measurements.

The largest differences are observed between the expected and actually measured PV voltages, with a maximum relative error of 3.11%. This occurs because the PV voltage (bus 3) varies considerably as the delivered active power changes. Small differences in accuracy arise from the monitoring devices and the difficulty of precisely obtaining at the same time the proposed values of active power generation from both the wind turbine and PV due to their inherent variability depending on the weather conditions and technical restrictions.

The results confirm that the AC/DC optimal power flow works correctly and accurately enough to assess the operation of the microgrid.

5.2 Techno-economic assessment

The techno-economic assessment calculations have been implemented in Python. Again, the calculation time is very short, only a few seconds. Four scenarios in which the TIGON CEDER microgrid could operate have been studied: baseline, no battery, battery flexibility, and virtual battery (VB) flexibility.

In the baseline scenario, the microgrid elements have been considered with the same characteristics presented in Section 4. Nevertheless, in this situation, which is real, the battery is being used only for performing tests and research experiments. For this reason, we have considered three more scenarios in which the microgrid could be more profitable. In the no-battery scenario, the battery and its converter have been removed from the microgrid, along with their investment costs and residual values. In the battery flexibility scenario, the battery has not been eliminated, but it has been deemed to be used in the Spanish upward tertiary regulation market. Batteries are loads with inherent electrical flexibility as they can be charged and discharged at the most convenient time, keeping the daily balance of the generated and consumed electricity unchanged. The upward tertiary regulation market has been chosen because it is the balance market with the highest average price in Spain in 2022: 224.17 €/MWh3. The power of the battery is very small to participate in the Spanish tertiary regulation market on its own as it is necessary to make bids of at least 1 MW (). Therefore, it needs to be part of an aggregation that participates in the market as a unitary market agent, following the methodology explained by . Lastly, in the VB flexibility scenario, the battery has been replaced by a VB consisting of an aggregation of thermostatically controlled loads (TCLs) already in place, with the same nominal power . These TCLs have the capacity to store thermal energy equivalent to a given electricity consumption. They could be heat pumps or water heaters, for instance. The aggregation of TCLs in the VB has been achieved using the method proposed by . This method is a control algorithm executed in real time and consists of three steps: Check of TCLs, Aggregation, and Priority Control. As a result of these stages, the controller decides which TCLs should be switched on or off to meet the stated power requirement. For both flexibility scenarios, storage participation in the market for 1 h per day is assumed to operate at its nominal power and the aggregation shares its benefits proportionally between its loads. This results in an additional flexibility income of 2,045.55 € per year.

The KPI values obtained for each scenario are given in Table 8.

TABLE 8

ScenarioKPI1 (kWh)KPI2 (kgCO2)KPI3 (%)KPI4 (kW)
Baseline80,4252,346102.0125
No battery80,4252,346102.010
Battery flexibility80,4252,346102.0125
VB flexibility80,4252,346102.0125
ScenarioKPI5 (€)KPI6 (€)KPI7 (years)KPI8 (€/kWh)
Baseline256,824.77260,635.35Never145.73
No Battery256,824.77204,286.5621113.92
Battery301,874.24260,635.3524145.73
VB flexibility301,874.24204,286.5617113.92

KPI results.

KPI2 (CO2) proves that carbon dioxide emissions are very low. This is caused by the usage of renewable generation technologies, specifically PV and wind turbines. KPI3 (self-consumption percentage) shows that the microgrid produces more energy than required by the loads per year, thus compensating for energy losses.

KPI7 (payback) shows that in the baseline scenario, the microgrid is not profitable. It makes sense as the TIGON CEDER microgrid is a small-scale network where the elements are high-cost prototypes. The battery, which is the storage element of the grid, is key for increasing or decreasing costs and incomes. The battery flexibility scenario demonstrates that using the battery for making bids in the Spanish electricity balance markets enhances the incomes and makes the microgrid profitable. Nevertheless, the flexibility remuneration is not high enough to make the microgrid payback lower than that when the battery is removed, i.e., the no-battery scenario. In the no-battery scenario, costs are highly reduced, as well as in KPI8 (LCOE). This is at the expense of reducing KPI4 (storage flexibility) to zero. The most profitable scenario is VB flexibility as it combines the additional income obtained by participation in the tertiary regulation market with the cost savings caused by the disappearance of the physical battery. KPI7 (payback) is reduced to 17 years and KPI8 (LCOE) to 113.92 €/kWh, whereas KPI4 (storage flexibility) remains at 25 kW.

6 Conclusion

In this paper, an AC/DC optimal power flow and a techno-economic assessment have been presented with the aim of helping users evaluate the operation and viability of their hybrid AC/DC microgrids. In addition, both calculations can be used for proposing improvements and new investments. The AC/DC optimal power flow is a useful technique for checking the correct operation of hybrid microgrids under specific instantaneous conditions, while the techno-economic assessment allows us to verify their performance in the long term.

Both methodologies have been tested on the TIGON CEDER microgrid, which has been described in this paper. The AC/DC optimal power flow has been validated using real measurements and considering different objective function scenarios. Among the four techno-economic scenarios, the VB flexibility scenario is the most profitable and offers the fastest return on investment. Based on the results of flexibility scenarios, we encourage the stakeholders in the Spanish electricity system to increase the remuneration of flexibility to increase the penetration of renewable energy sources and advance the transition to a more efficient AC/DC hybrid power system with low greenhouse gas emissions.

In future work, we will apply the AC/DC optimal power flow and the techno-economic assessment to other microgrids, especially those primarily aimed at supplying electricity to residential consumers and industries. We also hope that the two developed techniques and the TIGON CEDER microgrid example can be used by other companies and research institutions in their developments.

Statements

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

AM-C: conceptualization, data curation, formal analysis, investigation, methodology, software, validation, visualization, writing–original draft, and writing–review and editing. AH-S: conceptualization, formal analysis, investigation, methodology, project administration, software, and writing–review and editing. ÓI-M: conceptualization, data curation, formal analysis, investigation, project administration, resources, validation, visualization, and writing–review and editing. PP-C: investigation, resources, and writing–review and editing. ÁH-J: investigation, resources, and writing–review and editing. FF-E: conceptualization, methodology, supervision, and writing–review and editing. EB: conceptualization, methodology, supervision, and writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research received funding from the European Union’s Horizon 2020 TIGON project under grant agreement no. 957769.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Nomenclature

Bus voltage phase (rad)
Converter performance
Simulated value
Set of buses
Microgrid topology network graph
Set of lines
Microgrid angular frequency (rad/s)
Measured value
Line susceptance to ground
Line capacitance (F)
Line conductance
Generator capacity factor (%)
Real part of bus voltage (V)
Electricity price (currency/kWh)
Imaginary part of bus voltage (V)
Flexibility income (currency)
Generator CO2 equivalent (kgCO2/kWh)
Subindices denoting buses
Bus injected current (A)
Current flowing to the ground between i and k (A)
Current flowing from i to k or vice versa (A)
Line current limit (A)
Maximum current flowing from i to k or vice versa (A)
Microgrid investment cost (currency)
Generator investment cost (currency)
Imaginary unit
Generator maintenance cost (currency/year)
Number of measurements
Subindex denoting years
Generator operational cost (currency/kWh)
Microgrid operation and maintenance cost (currency/year)
Bus total active power (W)
Converter injected or transferred active power (W)
Generator active power (W)
Generator minimum active power (W)
Generator nominal active power (W)
Load active power (W)
Storage nominal power (W)
Bus total reactive power (VAr)
Generator reactive power (W)
Generator minimum active power (VAr)
Generator nominal active power (VAr)
Load reactive power (W)
Discount rate (%)
Line resistance
Microgrid residual value (currency)
Generator residual value (currency)
Bus apparent power (VA)
Base power (VA)
Line susceptance
Transformer or converter nominal power (VA)
Microgrid useful life (years)
Bus voltage (V)
Base voltage (V)
Transformer zero-sequence relative short-circuit voltage percentage (%)
Transformer nominal voltage (V)
Bus maximum voltage (V)
Bus minimum voltage (V)
Bus nominal voltage (V)
Transformer-resistive part of the zero-sequence relative to short-circuit voltage percentage (%)
Line reactance
Line admittance
Base impedance
Line impedance

References

  • 1

    AbadieL.ChamorroJ. (2019). Levelized cost of electricity: key drivers and valuation methods. Dyna94, 656661. 10.6036/9223

  • 2

    Agundis-TinajeroG.AldanaN. L. D.LunaA. C.Segundo-RamírezJ.Visairo-CruzN.GuerreroJ. M.et al (2018). Extended-optimal-power-flow-based hierarchical control for islanded ac microgrids. IEEE Trans. Power Electron.34, 840848. 10.1109/tpel.2018.2813980

  • 3

    Agundis-TinajeroG.Segundo-RamírezJ.Visairo-CruzN.SavaghebiM.GuerreroJ. M.BarocioE. (2019). Power flow modeling of islanded ac microgrids with hierarchical control. Int. J. Electr. Power and Energy Syst.105, 2836. 10.1016/j.ijepes.2018.08.002

  • 4

    AlexanderC. K.SadikuM. N. (2013). Fundamentos de circuitos elétricos. McGraw-Hill Interamericana de España S.L.

  • 5

    AnderssonJ. A. E.GillisJ.HornG.RawlingsJ. B.DiehlM. (2019). CasADi – a software framework for nonlinear optimization and optimal control. Math. Program. Comput.11, 136. 10.1007/s12532-018-0139-41007/s12532-018-0139-4

  • 6

    AshiqueR. H.SalamZ.AzizM. J. B. A.BhattiA. R. (2017). Integrated photovoltaic-grid dc fast charging system for electric vehicle: a review of the architecture and control. Renew. Sustain. Energy Rev.69, 12431257. 10.1016/j.rser.2016.11.245

  • 7

    BaldwinS. (2006). Carbon footprint of electricity generation. Parliam. Office Sci. Technol.268.

  • 8

    ChopraS.VanaprasadG. M.TinajeroG. D. A.BazmohammadiN.VasquezJ. C.GuerreroJ. M. (2022). Power-flow-based energy management of hierarchically controlled islanded ac microgrids. Int. J. Electr. Power and Energy Syst.141, 108140. 10.1016/j.ijepes.2022.108140

  • 9

    CuiY. (2017). “Static parameter detection and identification of power grid branch based on PMU measurement,” in 2016 4th international conference on machinery, materials and information technology applications (Atlantis Press), 313318.

  • 10

    FotopoulouM.RakopoulosD.TrigkasD.StergiopoulosF.BlanasO.VoutetakisS. (2021). State of the art of low and medium voltage direct current (DC) microgrids. Energies14, 5595. 10.3390/en14185595

  • 11

    HatziargyriouN.AsanoH.IravaniR.MarnayC. (2007). Microgrids. IEEE power energy Mag.5, 7894. 10.1109/mpae.2007.376583

  • 12

    KatiraeiF.IravaniR.HatziargyriouN.DimeasA. (2008). Microgrids management. IEEE power energy Mag.6, 5465. 10.1109/mpe.2008.918702

  • 13

    KiranD. (2022). Principles of economics and management for manufacturing engineering. Butterworth-Heinemann.

  • 14

    LasseterR. H. (2002). “Microgrids,” in 2002 IEEE power engineering society winter meeting. Conference proceedings (Cat. No. 02CH37309) (IEEE), 1, 305308. 10.1109/pesw.2002.985003

  • 15

    MarnayC.ChatzivasileiadisS.AbbeyC.IravaniR.JoosG.LombardiP.et al (2015). “Microgrid evolution roadmap,” in 2015 international symposium on smart electric distribution systems and technologies (EDST) (IEEE), 139144.

  • 16

    Martín-CrespoA.BaeyensE.Saludes-RodilS.Frechoso-EscuderoF. (2023). Aggregated demand flexibility prediction of residential thermostatically controlled loads and participation in electricity balance markets. arXiv preprint arXiv:2308.04853.

  • 17

    Martín-CrespoA.Saludes-RodilS.BaeyensE. (2021). Flexibility management with virtual batteries of thermostatically controlled loads: real-time control system and potential in Spain. Energies14, 1711. 10.3390/en14061711

  • 18

    MohammedS. A. Q.JungJ.-W. (2021). A state-of-the-art review on soft-switching techniques for DC–DC, DC–AC, AC–DC, and AC–AC power converters. IEEE Trans. Industrial Inf.17, 65696582. 10.1109/tii.2021.3058218

  • 19

    MontesA.CastroJ. (1995). Solving the load flow problem using gröbner basis. ACM SIGSAM Bull.29, 113. 10.1145/216685.216686

  • 20

    PapapetrouM.CipollinaA.La CommareU.MicaleG.ZaragozaG.KosmadakisG. (2017). Assessment of methodologies and data used to calculate desalination costs. Desalination419, 819. 10.1016/j.desal.2017.05.038

  • 21

    PiresV. F.PiresA.CordeiroA. (2023). DC microgrids: benefits, architectures, perspectives and challenges. Energies16, 1217. 10.3390/en16031217

  • 22

    QuezadaV. M.AbbadJ. R.RomanT. G. S. (2006). Assessment of energy distribution losses for increasing penetration of distributed generation. IEEE Trans. power Syst.21, 533540. 10.1109/tpwrs.2006.873115

  • 23

    RaufS.WahabA.RizwanM.RasoolS.KhanN. (2016). Application of dc-grid for efficient use of solar PV system in smart grid. Procedia Comput. Sci.83, 902906. 10.1016/j.procs.2016.04.182

  • 24

    Red Eléctrica de España (2021). Guía descriptiva. Ser Proveed. Serv. balance. Versión 03. (in Spanish).

  • 25

    SaeedM. H.FangzongW.KalwarB. A.IqbalS. (2021). A review on microgrids’ challenges and perspectives. IEEE Access9, 166502166517. 10.1109/access.2021.3135083

  • 26

    SaeedifardM.GraovacM.DiasR.IravaniR. (2010). “DC power systems: challenges and opportunities,” in IEEE PES general meeting (IEEE), 17.

  • 27

    SalamA.MohamedA.HannanM. (2008). Technical challenges on microgrids. ARPN J. Eng. Appl. Sci.3, 6469.

  • 28

    SamperioÁ. (2023). Sparse recovery of an electrical network. arXiv preprint arXiv:2304.06676.

  • 29

    SartoriD.CatalanoG.GencoM.PancottiC.SirtoriE.VignettiS.et al (2014). Guide to cost-benefit analysis of investment projects. Economic appraisal tool for cohesion policy 2014-2020, European Commission. Directorate-General for Regional and Urban policy. Tech. Rep.

  • 30

    ShaoS. J.AgelidisV. G. (2010). Review of DC system technologies for large scale integration of wind energy systems with electricity grids. Energies3, 13031319. 10.3390/en3061303

  • 31

    ShiR.SemsarS.LehnP. W. (2017). Constant current fast charging of electric vehicles via a DC grid using a dual-inverter drive. IEEE Trans. Industrial Electron.64, 69406949. 10.1109/tie.2017.2686362

  • 32

    TinajeroG. D. A.NasirM.VasquezJ. C.GuerreroJ. M. (2021). Comprehensive power flow modelling of hierarchically controlled AC/DC hybrid islanded microgrids. Int. J. Electr. Power and Energy Syst.127, 106629. 10.1016/j.ijepes.2020.106629

  • 33

    WächterA.BieglerL. T. (2006). On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming. Math. Program.106, 2557. 10.1007/s10107-004-0559-y

  • 34

    WangP.GoelL.LiuX.ChooF. H. (2013). Harmonizing AC and DC: a hybrid AC/DC future grid solution. IEEE Power Energy Mag.11, 7683. 10.1109/mpe.2013.2245587

  • 35

    ZubietaL. E. (2016). Are microgrids the future of energy? DC microgrids from concept to demonstration to deployment. IEEE Electrification Mag.4, 3744. 10.1109/mele.2016.2544238

Appendix A

Microgrid Model

The model representation of grid lines (Figure A1 (; )) is used for performing the AC/DC optimal power flow. In the case of DC lines, and are not considered, as the reactive part of impedances does not affect electric charges when they move always in the same direction over time.

FIGURE A1

is the inverse of the line impedance , which is expressed by Equation A1.

Both and are directly proportional to the line length, and, in the case of , also the frequency. In the case of DC lines, is not considered because it is a reactive parameter. For AC cases, is obtained using Equation A2.

is directly proportional to the line length. Each line has a current limit .

Voltage at each bus can be expressed by its real and imaginary part, as in Equation A3

The bus voltage phase can be calculated by Equation A4.

Each bus of the grid can have a different reference voltage. In order to ease load flow calculations in grids where voltage and current transformations occur through converters or transformers, a base power and a base voltage are set. Base impedance is calculated as in Equation A5.

All magnitudes in the grid model are divided by their base values, so they can be expressed per unit (p.u.). Maximum and minimum voltages at each bus are specified ( and , respectively).

Transformers in the grid are modelled in the same way as lines, but without taking into account because it is considered to have a minimal impact on the total transformer impedance. The magnitude and the resistive part of their are calculated with Equations A6, A7.

Besides transformers, converters are the other grid elements that changes voltage levels between buses. There exists AC/DC, DC/AC, AC/AC and DC/DC converters (). Their characteristics needed to perform load flow calculations are nominal power , performance and control mode (grid-following or grid-forming). The reactive power at the output of DC/AC and AC/AC converters is included in the grid model of this paper as a reactive power generator in order to recreate the reactive power management operation of these converters.

All generators existing in the microgrid are characterized by their nominal active power , minimum active power , nominal reactive power and minimum reactive power . Regarding loads, they are defined by two parameters: demanded active power and demanded reactive power . Supplying and is considered mandatory. In the case of storages, they can act as generators of loads in the power flow. Their defining parameter is .

If present, external grids are considered to be capable of providing and consuming active power and reactive power with no limits, but not at the same time.

For the techno-economic assessment, additional information about generators is considered. This includes their investment cost , residual value , maintenance cost , operational cost , capacity factor and CO2 equivalent . is the percentage of power produced by the generator with respect to the power that could have been produced if the generator had operated at maximum rate ().

Regarding the rest of the microgrid, investment cost , residual value , operation and maintenance cost and useful life are used. In addition, the price at which the microgrid sells electricity ,and the discount rate are included in the model. Flexibility income is considered in the flexibility scenario of the techno-economic analysis.

Summary

Keywords

AC/DC optimal power flow, hybrid microgrids, key performance indicators, techno-economic assessment, polynomial optimization, Python

Citation

Martín-Crespo A, Hernández-Serrano A, Izquierdo-Monge Ó, Peña-Carro P, Hernández-Jiménez Á, Frechoso-Escudero F and Baeyens E (2024) AC/DC optimal power flow and techno-economic assessment for hybrid microgrids: TIGON CEDER demonstrator. Front. Energy Res. 12:1399114. doi: 10.3389/fenrg.2024.1399114

Received

11 March 2024

Accepted

31 July 2024

Published

30 August 2024

Volume

12 - 2024

Edited by

George Tsekouras, University of West Attica, Greece

Reviewed by

Pavlos Georgilakis, National Technical University of Athens, Greece

Roberto Faranda, Polytechnic University of Milan, Italy

Updates

Copyright

*Correspondence: Alejandro Martín-Crespo,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics