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 6–7 (; ).
, 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 1–3.
FIGURE 2
TABLE 1
| Bus | (kV) | Type | |||
|---|---|---|---|---|---|
| 0 | 15.000 | AC | |||
| 1 | 0.400 | AC | |||
| 2 | 0.400 | AC | |||
| 3 | 0.630 | DC | |||
| 4 | 0.800 | DC | |||
| 5 | 0.860 | DC | |||
| 6 | 0.800 | DC | |||
| 7 | 0.230 | AC | |||
| 8 | 0.630 | AC |
| Line | Length (km) | (/km) | (/km) | ||
|---|---|---|---|---|---|
| 0 | (2,1) | 0.15 | 0.5 | 0.35 | |
| 1 | (4,6) | 0.20 | 0.5 | – |
| Line | (nF/km) | (kA) | Type | ||
|---|---|---|---|---|---|
| 0 | (2,1) | 100 | 100 | AC | |
| 1 | (4,6) | – | 100 | DC |
| Transformer | (kVA) | (%) | (%) | ||
|---|---|---|---|---|---|
| 0 | (1,0) | 250 | 1.5 | 1.0 |
| Converter | (kVA) | Control | |||
|---|---|---|---|---|---|
| 0 | (3,2) | 20 | 0.99 | Grid following | |
| 1 | (3,4) | 20 | 0.99 | Grid forming | |
| 2 | (4,5) | 30 | 0.86 | Grid forming | |
| 3 | (4,7) | 12 | 0.99 | Grid forming | |
| 4 | (8,4) | 5 | 0.99 | Grid following |
| Generator | Bus | (kW) | (kW) | (kW) | (kW) |
|---|---|---|---|---|---|
| 0 | 3 | 22.14 | 0.5 | 0 | 0 |
| 1 | 8 | 4.20 | 0.0 | 10 | −10 |
| 2 | 2 | 0.00 | 0.0 | 10 | −10 |
| 3 | 7 | 0.00 | 0.0 | 10 | −10 |
| Load | Bus | (kW) | (kW) | ||
|---|---|---|---|---|---|
| 0 | 6 | 5 | 0 | ||
| 1 | 7 | 4 | 0 |
| Storage | Bus | (kW) | |||
|---|---|---|---|---|---|
| 0 | 5 | 25 |
TIGON CEDER electrical data.
TABLE 2
| Generator | Bus | (€) | (€) | (€/year) |
|---|---|---|---|---|
| 0 | 3 | 23,000 | 100 | 200 |
| 1 | 8 | 10,000 | 150 | 600 |
| Generator | Bus | (€/kWh) | (%) | (kgCO2/kWh) |
|---|---|---|---|---|
| 0 | 3 | 0.003 | 33.5 | 0.03500 |
| 1 | 8 | 0.008 | 42.0 | 0.00464 |
Techno-economic information of TIGON CEDER generators.
TABLE 3
| (€) | (€) | (€/year) | (%) | (years) | (€/kWh) |
|---|---|---|---|---|---|
| 195,500 | 28,900 | 1,400 | 1 | 25 | 0.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
| Equipment | Investment cost (€) | Residual value (€) |
|---|---|---|
| DC PV converter | 23,000 | 2,500 |
| AC PV converter | 22,000 | 2,000 |
| Wind turbine converter | 10,000 | 500 |
| Battery | 34,000 | 400 |
| Battery converter | 25,000 | 3,000 |
| AC load converter | 31,000 | 2,500 |
| Wiring | 15,500 | 4,000 |
| Cabins | 18,000 | 18,000 |
| Labor | 17,000 | 0 |
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) |
|---|---|---|---|---|---|
| 0 | H1 | 0.00000 | 0.00000 | 14.82807 | 0.06733 |
| 1 | H1 | 0.00000 | 0.00000 | 0.39542 | 0.06733 |
| 2 | H1 | 0.00000 | −0.00074 | 0.39542 | 0.06733 |
| 3 | H1 | 19.75450 | 0.00000 | 0.62046 | 0.00000 |
| 4 | H1 | 0.00000 | 0.00000 | 0.80000 | 0.00000 |
| 5 | H1 | −12.50000 | 0.00000 | 0.86000 | 0.00000 |
| 6 | H1 | −5.00000 | 0.00000 | 0.79937 | 0.00000 |
| 7 | H1 | −4.00000 | 0.00000 | 0.23000 | 0.00000 |
| 8 | H1 | 4.06289 | 0.00000 | 0.62046 | 0.00000 |
| 0 | H2 | −0.05089 | 0.00000 | 14.99987 | 0.59858 |
| 1 | H2 | 0.00000 | 0.00000 | 0.40000 | 0.59858 |
| 2 | H2 | 0.00000 | −0.00075 | 0.40001 | 0.59860 |
| 3 | H2 | 19.80580 | 0.00000 | 0.63000 | 0.00000 |
| 4 | H2 | 0.00000 | 0.00000 | 0.80000 | 0.00000 |
| 5 | H2 | −12.50000 | 0.00000 | 0.86000 | 0.00000 |
| 6 | H2 | −5.00000 | 0.00000 | 0.79937 | 0.00000 |
| 7 | H2 | −4.00000 | 0.00000 | 0.23000 | 0.00000 |
| 8 | H2 | 4.06300 | 0.00000 | 0.63000 | 0.00000 |
| 0 | H3 | 0.00000 | 0.00000 | 14.95132 | 0.06573 |
| 1 | H3 | 0.00000 | 0.00000 | 0.39870 | 0.06573 |
| 2 | H3 | 0.00000 | −0.00075 | 0.39870 | 0.06573 |
| 3 | H3 | 20.00000 | 0.00000 | 0.62046 | 0.00000 |
| 4 | H3 | 0.00000 | 0.00000 | 0.80000 | 0.00000 |
| 5 | H3 | −12.50000 | 0.00000 | 0.86000 | 0.00000 |
| 6 | H3 | −5.00000 | 0.00000 | 0.79937 | 0.00000 |
| 7 | H3 | −4.00000 | 0.00000 | 0.23000 | 0.00000 |
| 8 | H3 | 3.81736 | 0.00000 | 0.62046 | 0.00000 |
| 0 | H4 | −2.49415 | 0.00000 | 14.75867 | 0.01449 |
| 1 | H4 | 0.00000 | 0.00000 | 0.39361 | 0.01460 |
| 2 | H4 | 0.00000 | 0.00166 | 0.39408 | 0.01545 |
| 3 | H4 | 22.14000 | 0.00000 | 0.62046 | 0.00000 |
| 4 | H4 | 0.00000 | 0.00000 | 0.80000 | 0.00000 |
| 5 | H4 | −12.50000 | 0.00000 | 0.86000 | 0.00000 |
| 6 | H4 | −5.00000 | 0.00000 | 0.79937 | 0.00000 |
| 7 | H4 | −4.00000 | 0.00000 | 0.23000 | 0.00000 |
| 8 | H4 | 4.20000 | 0.00000 | 0.62046 | 0.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) |
|---|---|---|---|---|---|---|
| 0 | H1 | 0.000 | 14.790 | H2 | −0.050 | 14.990 |
| 1 | H1 | 0.000 | 0.399 | H2 | 0.000 | 0.399 |
| 2 | H1 | 0.000 | 0.399 | H2 | 0.000 | 0.399 |
| 3 | H1 | 19.750 | 0.612 | H2 | 19.800 | 0.611 |
| 4 | H1 | 0.000 | 0.799 | H2 | 0.000 | 0.808 |
| 5 | H1 | −12.500 | 0.859 | H2 | −12.500 | 0.862 |
| 6 | H1 | −5.000 | 0.799 | H2 | −5.000 | 0.799 |
| 7 | H1 | −4.000 | 0.230 | H2 | −4.000 | 0.230 |
| 8 | H1 | 4.100 | 0.629 | H2 | 4.102 | 0.630 |
| 0 | H3 | 0.000 | 14.920 | H4 | −2.490 | 14.880 |
| 1 | H3 | 0.000 | 0.398 | H4 | 0.000 | 0.399 |
| 2 | H3 | 0.000 | 0.398 | H4 | 0.000 | 0.399 |
| 3 | H3 | 20.000 | 0.613 | H4 | 22.140 | 0.632 |
| 4 | H3 | 0.000 | 0.806 | H4 | 0.000 | 0.805 |
| 5 | H3 | −12.500 | 0.860 | H4 | −12.500 | 0.865 |
| 6 | H3 | −5.000 | 0.801 | H4 | −5.000 | 0.799 |
| 7 | H3 | −4.000 | 0.230 | H4 | −4.000 | 0.230 |
| 8 | H3 | 3.800 | 0.628 | H4 | 4.200 | 0.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 (%) | H1 | H2 | H3 | H4 |
|---|---|---|---|---|
| Active power | 0.10 | 0.31 | 0.05 | 0.02 |
| Voltage | 0.56 | 0.55 | 0.44 | 0.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
| Scenario | KPI1 (kWh) | KPI2 (kgCO2) | KPI3 (%) | KPI4 (kW) |
|---|---|---|---|---|
| Baseline | 80,425 | 2,346 | 102.01 | 25 |
| No battery | 80,425 | 2,346 | 102.01 | 0 |
| Battery flexibility | 80,425 | 2,346 | 102.01 | 25 |
| VB flexibility | 80,425 | 2,346 | 102.01 | 25 |
| Scenario | KPI5 (€) | KPI6 (€) | KPI7 (years) | KPI8 (€/kWh) |
|---|---|---|---|---|
| Baseline | 256,824.77 | 260,635.35 | Never | 145.73 |
| No Battery | 256,824.77 | 204,286.56 | 21 | 113.92 |
| Battery | 301,874.24 | 260,635.35 | 24 | 145.73 |
| VB flexibility | 301,874.24 | 204,286.56 | 17 | 113.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 |
Footnotes
2.^https://www.ciemat.es/portal.do?IDM=6&NM=1
3.^https://www.sistemaelectrico-ree.es/informe-del-sistema-electrico/mercados/servicios-ajuste/energias-precios-balance
References
1
AbadieL.ChamorroJ. (2019). Levelized cost of electricity: key drivers and valuation methods. Dyna94, 656–661. 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, 840–848. 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, 28–36. 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, 1–36. 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, 1243–1257. 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), 313–318.
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, 78–94. 10.1109/mpae.2007.376583
12
KatiraeiF.IravaniR.HatziargyriouN.DimeasA. (2008). Microgrids management. IEEE power energy Mag.6, 54–65. 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, 305–308. 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), 139–144.
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, 6569–6582. 10.1109/tii.2021.3058218
19
MontesA.CastroJ. (1995). Solving the load flow problem using gröbner basis. ACM SIGSAM Bull.29, 1–13. 10.1145/216685.216686
20
PapapetrouM.CipollinaA.La CommareU.MicaleG.ZaragozaG.KosmadakisG. (2017). Assessment of methodologies and data used to calculate desalination costs. Desalination419, 8–19. 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, 533–540. 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, 902–906. 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, 166502–166517. 10.1109/access.2021.3135083
26
SaeedifardM.GraovacM.DiasR.IravaniR. (2010). “DC power systems: challenges and opportunities,” in IEEE PES general meeting (IEEE), 1–7.
27
SalamA.MohamedA.HannanM. (2008). Technical challenges on microgrids. ARPN J. Eng. Appl. Sci.3, 64–69.
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, 1303–1319. 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, 6940–6949. 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, 25–57. 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, 76–83. 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, 37–44. 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
© 2024 Martín-Crespo, Hernández-Serrano, Izquierdo-Monge, Peña-Carro, Hernández-Jiménez, Frechoso-Escudero and Baeyens.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alejandro Martín-Crespo, alemar@cartif.es
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.