Next Article in Journal
Time-Dependent Seismic Hazard Analysis for Induced Seismicity: The Case of St Gallen (Switzerland), Geothermal Field
Next Article in Special Issue
Multi-Step Energy Demand and Generation Forecasting with Confidence Used for Specification-Free Aggregate Demand Optimization
Previous Article in Journal
Increased Internal Combustion Engine Efficiency with Optimized Valve Timings in Extended Stroke Operation
Previous Article in Special Issue
Price Forecasting for the Balancing Energy Market Using Machine-Learning Regression
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

OptiMEMS: An Adaptive Lightweight Optimal Microgrid Energy Management System Based on the Novel Virtual Distributed Energy Resources in Real-Life Demonstration

by
Angelina D. Bintoudi
*,
Lampros Zyglakis
,
Apostolos C. Tsolakis
,
Paschalis A. Gkaidatzis
,
Athanasios Tryferidis
,
Dimosthenis Ioannidis
* and
Dimitrios Tzovaras
Centre for Research and Technology—Hellas, Information Technologies Institute, P.O. Box 60361, 5700 Thessaloniki, Greece
*
Authors to whom correspondence should be addressed.
Energies 2021, 14(10), 2752; https://doi.org/10.3390/en14102752
Submission received: 2 April 2021 / Revised: 20 April 2021 / Accepted: 8 May 2021 / Published: 11 May 2021
(This article belongs to the Special Issue Control and Optimization of Renewable Energy Systems)

Abstract

:
As microgrids have gained increasing attention over the last decade, more and more applications have emerged, ranging from islanded remote infrastructures to active building blocks of smart grids. To optimally manage the various microgrid assets towards maximum profit, while taking into account reliability and stability, it is essential to properly schedule the overall operation. To that end, this paper presents an optimal scheduling framework for microgrids both for day-ahead and real-time operation. In terms of real-time, this framework evaluates the real-time operation and, based on deviations, it re-optimises the schedule dynamically in order to continuously provide the best possible solution in terms of economic benefit and energy management. To assess the solution, the designed framework has been deployed to a real-life microgrid establishment consisting of residential loads, a PV array and a storage unit. Results demonstrate not only the benefits of the day-ahead optimal scheduling, but also the importance of dynamic re-optimisation when deviations occur between forecasted and real-time values. Given the intermittency of PV generation as well as the stochastic nature of consumption, real-time adaptation leads to significantly improved results.

1. Introduction

Microgrids (MG), foreseen as the avant-garde building blocks of the future Distribution Network (DN), have attracted much attention from various research fields, ranging from power electronics [1] and control systems [2] to machine learning [3]. As such, it is essential to move past simulations [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23] and hardware-in-the-loop (HIL) configurations [24,25,26,27,28,29,30] towards the evaluation of technologies in real-life applications [31,32,33,34]. In general, MG operation must be technically reliable and economically feasible, regardless of the stochastic nature of Renewable Energy Sources (RES), or the connection status with the DN [16,35]. The MG day-ahead scheduling problem essentially consists of the implementation of an optimised energy management strategy. Optimisation techniques lead to improved MG operation compared to rule-based approaches [36], which, though simple and direct, limit the MG operation [24]. The optimisation objectives may vary: cost [37] or profit-efficient MG operation [38], via stochastic programming [39,40], or a combination of integrated Two-Stage Stochastic Programming (ITSP) and Worst-Case Value at Risk (VCVaR) [18]; social benefit maximisation via Mixed Integer Point Solver (MIPS) [14]; extension of equipment’s lifetime [8,41]; load feeding reliability, utilisation of resources, e.g., RES production [42]. They are solved using common optimisation techniques, e.g., Dynamic, Linear/Mixed-Integer Linear Programming (LP/MILP) [43]; Nonlinear (NLP)/Quadratic Programming [15,34]; or heuristics ones, such as Genetic Algorithm [44], Isolation Niche Immune Genetic Algorithm [45], Non-Dominated Sorting Genetic Algorithm type II [46], Niching Evolutionary Algorithm [47], Particle Swarm Optimisation [37,48], Ant Colony Optimisation [49], or even Binary Backtracking Search Algorithm [20] and Glowworm Swarm Optimisation [12]. Heuristic algorithms are rarely applied beyond simulation. Due to their structure, they are susceptible to local minima entrapment [50]. Moreover, they are quite complex to implement, leading to long computation time. Ergo, they are not quite deemed suitable for real-time operation but rather for cases of long-term planning and assessment. Simpler, deterministic approaches, such as LP/MILP have been identified as more suited for real-site demonstration [9,13,26,51,52].
As presented in Table 1, there is a significant gap in the literature regarding actual deployment of the introduced solutions, either in terms of Day-Ahead (DA) and Real-Time (RT) optimisation schemes, which could very well lead to completely different results than the ones extracted from simulations, or even HIL setups. In general, simulations and HIL setups are good at understanding system margins, making design trade-offs and verification. However, each simulation software makes a series of assumptions in regards to the equipment modelling and operating conditions. Similarly, small-scale experimental test-beds, though closer to reality, still constitute a controlled environment, which, especially in the presence of purely stochastic phenomena (such as loading and weather conditions), cannot describe reality accurately [53]. For these reasons, this work aims towards presenting the real-life application of a microgrid optimal scheduler, where real weather conditions and real home users are involved for a period of 9 days.
Aiming to create a real-world applied optimal scheduler, in this paper, an Optimal Energy Management System (OptiMEMS) has been developed, deployed and assessed upon a residential MG. OptiMEMS framework includes four components: an MILP-based optimisation engine that produces the optimal day- or steps-ahead MG schedules; a real-time supervising tool that monitors and validates the optimal schedule application; a machine-learning-based load consumption forecasting tool and a hybrid deterministic/stochastic PV generation predicting tool (Figure 1). The optimised scheduler solves an adjusted Unit Commitment (UC) problem, tailored for MGs operating in grid-connected mode, with the objective being the minimisation of the operation cost. Finally, a customised module is responsible for applying the optimised schedule whenever it is allowed and, for triggering its recalculation in real-time, in case of significant divergence between real load consumption, RES generation and the forecasted ones. The application of the OptiMEMS is realised on a fully RES-based, grid-connected AC microgrid [54].
This paper is organised as follows. In Section 2.1, the problem formulation is analysed. In Section 2.2, the System Architecture is presented. In Section 3.1, the experimental setup is described, and in Section 3.2, the results are presented. Finally, Section 4 is devoted to conclusions and future work.
Table 1. Day-ahead MG optimisation approaches.
Table 1. Day-ahead MG optimisation approaches.
Ref No.ObjectiveTechniqueMG ModeDERsImplementationSchedule
[4]EM - min costMILPIslanded Grid ConnectedPV, Wind Turbines, Fuel Cells, Micro Turbines, Diesel Generators, ESSSimulationsDA RT
[5]EM - min lossesDynamic ProgrammingGrid ConnectedWind Turbines, ESS, PV, ESSSimulationDA
[6]EM - min costMILPGrid ConnectedPV, Wind Turbines, ESSSimulationsDA
[7]EMMILP/NLPIsolatedPV, Wind Turbines, Fuel Cells, Micro Turbines, Diesel Generators, ESSSimulationsDA
[8]multi-objectiveMILPIsolatedPV, Diesel Generators, ESSSimulationsDA
[10]min cost (operation)MILPGrid ConnectedPV, Wind Turbines, ESSSimulationsDA
[11]min cost (operation)MILPGrid ConnectedPV, Wind Turbines, ESSSimulationsDA
[12]min losses + min fuel consumptionHeuristicGrid ConnectedPV, Micro TurbinesSimulationsDA
[14]Social benefit maxQuadratic ProgrammingGrid ConnectedPV, Wind Turbines, Fuel Cells, Micro Turbines, Pico Hydel ESSSimulationsDA RT
[15]min cost (operation)Quadratic ProgrammingGrid ConnectedPV, Wind Turbines, ESSSimulationsDA RT
[16]min cost (operation)Stochastic/Robust optimisationGrid ConnectedPV, Wind Turbines, Diesel Generators, ESSSimulationsDA
[17]min cost (operation)Quadratic ProgrammingGrid ConnectedPV, ESSSimulationsDA
[18]min costITSP and WCVaRGrid ConnectedPV, Wind Turbines, Micro Turbines, CHP, ESSSimulationsDA
[20]EMHeuristicGrid ConnectedPV, Wind Turbines, Fuel Cells, Diesel Generators, ESSSimulationsDA
[21]min cost (operation)MILPIsolatedHeavy/Light fuel units, PV, Wind Turbines, ESSSimulationsDA
[22]min cost (operation)MILPIslandedPV, Wind Turbines, Diesel Generators, ESSSimulationsDA Intra-DA
[23]min cost—max ESS utilisation efficiencyMILPGird ConnectedPV, Wind Turbines, Micro Turbines, Diesel Generators, ESSSimulationsHA MA
[24]min costDynamic ProgrammingGrid ConnectedPV, ESSHILDA
[25]min cost (operation)Rule-basedIsolatedPV, Fuel CellsHILDA
[26,27]min cost (operation)LP/MILPIslanded Grid ConnectedPV, Wind Turbines, ESSHILDA
[30]min cost (operation)MILPIslanded Grid ConnectedPV, Wind Turbines, Fuel Cells, ESSHILDA
[32]min costMILPIsolatedPV, Wind Turbines, microCHP, ESSReal-life PrototypeDA
[33]min cost (operation + emissions)Rule-basedIsolatedPV, Diesel GeneratorsSimulationsYearly
[34]min cost (operation + emissions)NLP/Sequential Quadratic ProgrammingIslanded Grid ConnectedPV, Wind Turbines, Micro Turbines, CHP, ESSHILDA
[36]EM - min costRule-basedGrid ConnectedPV, Wind Turbines, Micro Turbines, ESSSimulations + HILDA RT
[39]min cost (operation)Stochastic ProgrammingGrid ConnectedPV, Wind Turbines, Micro Turbines, ESSSimulationsDA
[37,48]min cost (operation)HeuristicGrid ConnectedPV, Wind Turbines, CCHP, ESSSimulationsDA RT
[38]max profittwo-stage Stochastic ProgrammingGrid ConnectedPV, Wind Turbines, ESSSimulationsDA
[40]max power productionMINLPGrid ConnectedPV, ESSSimulations + HILDA RT
[41,51]min cost (operation + extension of hybrid ESS lifetime)MILPIsolatedPV, Diesel Generators, ESSSimulations + HILDA
[43]min fuel consumption + min cost (operation)LP/MILPIsolatedPV, Diesel Generators, ESSSimulationsDA
[44]min cost (operation)HeuristicIsolatedPV, Wind Turbines, Diesel Generators, ESSSimulationsDA
[45]min cost (operation)HeuristicIslanded Grid ConnectedPV, Wind Turbines, HT, Fuel Cells, ESS, GTSimulationDA
[46]min embodied energy + min LPSPNLP/Sequential Quadratic ProgrammingGrid ConnectedPV, Wind Turbines, ESSSimulationshours-ahead
[47]min cost (operation + emission)HeuristicGrid ConnectedESSSimulationshours-ahead
[49]min costHeuristicGrid ConnectedPV, Wind Turbines, Micro Turbines, ESSSimulations + HILDA RT
[52]max profitMILPGrid ConnectedPV, ESSSimulations + HILDA RT

2. Materials and Methods

2.1. Optimisation Problem Formulation

2.1.1. Virtual Distributed Energy Resources

In the context of OptiMEMS, the optimisation module is responsible for solving an adjusted UC problem tailored to MG requirements. Its basis is a novel concept called: virtual Distributed Energy Resource (vDER). Using the Passive Sign Convention, a vDER can be a power source (provider: negative power) or a power load (consumer: positive power). Therefore, in the present optimisation problem, the vDERs “ESS charge” and “Grid Export” are positive, while “PV”, “ESS Discharge” and “Grid Import” are negative; this information is given on the superscript symbols. The accumulated output of all vDERs must satisfy the load demand in each time slot.
Regarding the formulation of vDERs, unidirectional DERs, such as PV modules, remain as they were, whereas DERs with a bidirectional energy flow, such as ESSs, are split into two complementary virtual ones. In this study, the formulation of the vDERs concerns a grid-connected MG, consisting of PV modules, ESS modules (assigned index S) and the Point of Common Coupling (PCC), with the latter being treated as a bidirectional DER (assigned index G for Grid). For example, ESS Discharge and Grid Import are modelled as vDERs that provide power, and ESS Charge and Grid Export as vDERs that consume power. During grid-connected operation, aiming towards PV production maximisation, PV modules can be modelled as negative loads: forecasted production is subtracted from the forecasted consumption, leading to a new artificial curve of “Net Demand”. Generally, if PV curtailment is possible, PV modules are modelled as an additional vDER (negative power). In the following paragraphs, the latter, more general approach is employed.

2.1.2. Adjusted Unit Commitment Problem Formulation

Since the distances among the physical assets of the considered MG are limited, the physical aspects of the MG circuitry can be omitted and thus, the need for any kind of optimal power flow related constraints. Therefore, the single bus model approach of the UC can be adopted, without loss of reality [49,55]. It is noted that all the variables and symbols used in the following paragraphs are given in the Abbreviations section at the end of the document. The problem contains two sets of optimisation variables, shown below: the continuous vDER energies as a function of time and a set of auxiliary integer (binary) optimisation variables modelling the mutually exclusive operation of the complementary vDERs related to PCC and batteries. In the general case of an MG consisting of N batteries, M PV modules and a PCC, the optimisation variables for a time horizon T are:
E S ( t ) = E S 1 + ( t ) E S N + ( t ) E S 1 ( t ) E S N ( t ) = E S + ( t ) E S ( t )
E G ( t ) = E G ( t ) E G + ( t )
E P V ( t ) = E P V 1 ( t ) E P V M ( t )
β S ( t ) = β S 1 + ( t ) β S N + ( t ) β S 1 ( t ) β S N ( t ) = β S + ( t ) β S ( t )
β G ( t ) = β G ( t ) β G + ( t )
The optimisation objective, expressed by Equation (3), is to minimise the operation cost of a grid-connected MG. The main beneficiary can either be the MG owner or an aggregator depending on the proper energy price: a Real-Time Pricing in the case of the MG owner and the System Marginal Price for the aggregator. In this work, the MG owner is chosen as the end-user.
O F = min t = 1 T j G · ( E S ( t ) · L C O E S ) E G ( t ) · R T P ( t ) + E P V ( t ) · L C O E P V
The costs that have been taken into account are: a) the energy importing cost, b) the energy exporting price and c) the Levelised Cost of Energy (LCOE) for the PV and the ESS modules. Regarding costs (a) and (b), an RTP scheme has been used in order to generalise the problem formulation, which is defined as R T P ( t ) = [ R T P ( t ) R T P ( t ) + ] , t [ 1 , T ] . The LCOE calculations are done as in [56], since the usage of more elaborate schemes is considered to be out of the scope of the current study. All constraints are time-dependent since they should be applied in every time instant during the whole optimisation time horizon. Hence, for all the following equations: t [ 1 , T ] .
The first set of constraints in (4) expresses the equilibrium between energy supply and demand:
j G · ( E S ( t ) · j S ) + E G ( t ) · j G + E P V ( t ) · j P V = E L ( t )
The second constraint set (5) describes the vDERs minimum/maximum allowed limits, multiplied (using the ⊙ Hadamard, or the element-wise product) with the respective operation variables:
E S m i n β S + ( t ) E S + ( t ) E S m a x β S + ( t ) E S m i n β S ( t ) E S ( t ) E S m a x β S ( t ) E G m i n β G ( t ) E G ( t ) E G m a x β G ( t ) E P V m a x ( t ) E P V ( t ) 0
The set of constraints in (6) and (7) ensure that at all times, the energy stored in each battery follows the natural continuity and also, that the energy during both charge and discharge operation, is defined by preset limits, according to charge and discharge C-rates:
E S m i n τ = 0 t 1 E S + ( τ ) + E S ( τ ) + E S + ( t ) + E S ( t ) E S m a x
where, for τ = 0 , the sum E S + ( τ ) + E S ( τ ) corresponds to the initially stored energy ( E S ( 0 ) ) for the N microgrid batteries right before the beginning of the optimisation horizon start, and these values are given as input to the algorithm.
0 E S + ( t ) R C + C a p D o D · r t 0 E S ( t ) R C C a p D o D · ( r t )
The final constraint (8) dictates the mutually exclusive operation of the PCC and battery-related vDERs as follows:
β S + ( t ) + β S ( t ) = β G ( t ) + β G + ( t ) = j T
Since all optimisation expressions are linear and the optimisation variables are either continuous or binary, the optimisation problem is formed as an MILP problem. This is considered an asset, because LP/MILP problems are solved much faster than their non-linear equivalents. Moreover, re-running the optimisation algorithm can be achieved quickly, thus allowing for real-time applications. Since the OptiMEMS framework is accompanied by two auxiliary machine-learning based modules responsible for the forecast of the upcoming consumption and generation (presented briefly in the following subsection), there is no need to employ the Stochastic Optimisation approach since stochasticity is tackled in these two separate software modules.

2.2. System Architecture

2.2.1. Optimal Scheduling Engine for Microgrids

The OptiMEMS system architecture is depicted in Figure 1. It consists of software modules interconnected in the OptiMEMS platform, which interacts with web-based services in the cloud to fetch necessary input time series and also, with the MG assets for the setpoint signal dispatch and monitoring. The OptiMEMS platform is comprised of the following four sub-components: (a) one core module called Optimal Scheduling Engine for Microgrids (OSEM), which implements the Adjusted Unit Commitment Problem described in the previous section, (b) two forecasting engines providing the necessary production and generation time series forecasts to OSEM and finally, (c) the Real-Time Validator/Applicator of Optimal Schedule (RT-VAOS) module, which is responsible for the MG monitoring, dispatch of the produced setpoints and, whenever necessary, the trigger of recalculation of the optimal schedule. Automation-wise, all OptiMEMS components operate on a 24 h basis automatically, with 15-min granularity without human interference. In the following paragraphs, each module is described in terms of functional details, implementation and data flows.
The main component of OptiMEMS is the Optimal Scheduling Engine for MGs (OSEM) module, which formulates the optimisation problem described in the previous paragraph. OSEM is designed to produce an optimal schedule for either the whole day (24 h, 15 min of resolution), or upon request, from the next quarter of an hour till the end of the day. The OSEM engine is purely developed in Python 3.7. The mathematical formulation presented in Section 2.1 is realised with the help of the open-source library “PuLP” [57], while the MILP problem solver is the C++ open-source COIN-OR Branch and Cut (CBC) solver [58], characterised by high solving accuracy and computational speed. The choices of PuLP and CBC systems do not affect the functionalities of the OSEM module. Other optimisation frameworks, e.g., the usage of GUROBI, GAMS or MATLAB have also been tried; however, the authors chose the creation of an open-source optimisation engine.
The OSEM module needs three input variables: (a) the day-ahead RTP, (b) the load consumption and (c) PV production forecast time series, both with a predefined time granularity (e.g., quarterly, hourly etc.). When OptiMEMS runs at midnight, it receives the day-ahead forecasts, while, in case a recalculation of the optimised schedule is requested by RT-VAOS due to deviations, it receives the short-term forecasts, which correspond to a time series starting from the upcoming time slot up until the end of the day, and they are both updated periodically every hour so that even the latest deviations from normality are taken into consideration. The output of OSEM is a set of time series containing setpoints in terms of energy to be provided or consumed by the MG vDERs at each time slot of the optimisation horizon with the same granularity (e.g., quarterly), i.e., either day-ahead or for the remainder of the day in case recalculation is triggered.

2.2.2. Real-Time Validation and Application of Optimal Schedule

The module responsible for the application of day-ahead scheduling and coordination of its potential calibration is called “Real-Time Validation and Application of Optimal Schedule” (RT-VAOS). The term “real-time” is used in this work in order to describe the fact that at the moment of measurements logging (monitoring), control actions (either optimal schedule application or schedule bypass) are applied upon the microgrid inverters. RT-VAOS has a one-minute resolution, meaning that each minute measurements are monitored and control actions (following either optimal schedule or safety limits) are applied upon the microgrid inverters. The key feature of RT-VAOS is the use of real-time measurements to make the MG less vulnerable to possible errors deriving from forecasting tools. This component validates, in real-time, whether the current MG production and consumption are aligned or not (binary logic) with the pre-issued forecasted time series, within some tolerance limits. Their calculation has been performed via trial-and-error for various loading cases of the examined MG, leading to the creation of a customised “Tolerance Limits Look-up Table” (TL-LuT). In case of absence of significant deviations and if the MG DERs are capable of meeting the required energy demand, then appropriate setpoints are forwarded to the DER inverters. More specifically, RT-VAOS evaluates a series of predefined different scenarios in order to dispatch the proper signals. These scenarios are formed by the amount and type of vDERs defining the MG. RT-VAOS assigns each segment of the optimisation schedule in a specific case and calculates the forwarded setpoints to MG assets, which in all cases respect the battery’s safety limits.
In case of significant deviation of actual real-time overall load and/or PV measurements from the forecasted values, RT-VAOS enforces the “MG-AUTO” mode and consequently, demands recalculation of the optimised vDER scheduling for the remainder of the day, based on short-term forecasts. “MG-AUTO” operates the MG PV modules in their MPP, satisfies all connected loads and, whenever needed, charges the MG ESS, keeping them constantly fully charged. Any excess of energy is injected uncontrollably back to the grid. Figure 2 depicts the RT-VAOS algorithm flow chart.
Finally, RT-VAOS calculates the actual operating cost derived by the application of the optimal schedule, which, at the end of the day, is compared to the expected cost. This comparison is used for readjusting TL-LuT if necessary.

2.2.3. Load and PV Forecasting Engines

The MG load characteristics are highly related to seasonality, user behaviour and weather conditions. These parameters make load forecasting a challenging task, and extensive research has been applied to this area, a comprehensive review of which can be found in [3]. There is a rivalry between deterministic and machine-learning based approaches, with the latter prevailing in cases of limited knowledge, regarding each load, their types, the precise scheduling of MG users, etc. Thus, and due to the proposed OptiMEMS requirements, a machine-learning load-forecasting engine has been designed that exploits one-year-long real-time energy consumption measurements from an extensive smart-metering system of prime importance for a higher accuracy of the prediction. The engine produces a 24 h ahead load forecast with a 15-min resolution (long-term) once per day, at midnight. Day-ahead load forecast is considered a multi step time series forecasting problem; thus, multiple models for each time horizon are developed (Direct Multi-Step Forecast). Except for the load consumption, time features (hour, day, month, year) and temperature are also used as additional features in order to capture the seasonality and improve the efficiency of the module. Furthermore, this module has the functionality to run upon request in shorter time intervals (short-term, every quarter of an hour), thus reducing the mismatches between forecasted and real-time measurements. The load forecasting models are produced by using gradient boosting trees via the extreme gradient boosting (XGBoost) system [59]. For more information on the XGBoost model design, readers are encouraged to refer to the work of Xue et al. [60]. The design, implementation and testing of the load forecasting engine produced for the needs of OptiMEMS have been submitted by the authors for publication, and by the time this paper is written, it will be under review.
Regarding the PV forecast, a hybrid deterministic/machine-learning-based approach is adopted, as presented by the authors in previous works [61,62]. Briefly, the forecast engine is based on the open-source python library PVlib [63], adjustable to every PV installation, while two data-driven machine-learning models correct the forecasts PV production according to the expected weather forecast. More precisely, the first step is to feed a numerical weather prediction (NWP) to the PVlib engine in order to calculate the long-term (day-ahead) or short-term (few hours ahead) PV generation time series, with a 15-min time interval. The production curve is further improved by being fitted into two models that were trained utilising past data of actual and forecasted power generation. The first model takes into consideration the PV forecast errors with respect to the solar angles assuming ideal weather conditions, i.e., complete sunshine. The latter incorporates weather information concerning clouds, analysing past errors on different levels of cloudiness. Hence, the weather stochasticity and thus the PV volatile generation are taken into account within this sub-module. In terms of training the PV forecast models, a set of year-long data is required, which must contain an energy production time series from the PV modules of the selected MG. A detailed analysis regarding the hybrid deterministic/machine-learning-based PV forecasting module can be found in the recently published works of the authors of [61,62].

3. Experimental Validation

3.1. Experimental Setup

The application of the proposed OptiMEMS was realised on the CERTH/ITI SmartHome [54], a small-scale MG with AC-coupled DERs. This pilot is a real-life residential establishment with commercial (e.g., inverters) and custom equipment (e.g., smart meters, etc.) where various smart grid related technologies are evaluated for R&D purposes.

3.1.1. Microgrid Infrastructure

The MG (see Figure 3) consists of 9.57 kWp PV modules on a 10 kW inverter, a 5 kWh Li-ion battery with 3 single-phase 3 kVA inverters/chargers configured as Master-Slave-Slave in order to form a three-phase system, a variety of over 50 kWp controllable and non-controllable loads and a smart LV bus (fully controllable PCC). The system was configured to allow real-time monitoring and control on multiple levels using various wireless and wired communication protocols (e.g., MODBUS, TCP/IP). A detailed description of the CERTH/ITI SmartHome Microgrid can be found in [64]. This MG configuration allows control over the DER inverters as dictated by the proposed OptiMEMS under real conditions. Considering that the examined MG has AC-coupled DERs, and in order to highlight the performance of the algorithm without any load or PV curtailment, the mainly controlled assets are the battery inverters. A Raspberry Pi 3b running RT-VAOS was installed on-site in order to monitor and control the battery/PV inverters, along with the PCC.

3.1.2. Local Customisation—Web Services

As shown in Figure 1, both weather and RTP time series were fetched through online web-based Java APIs using RESTful services. The weather API collects current measurements every 15 min and the day-ahead predicted conditions every hour in order to be used as input at the two forecasting engines. The necessary variables are local temperature, irradiance (Diffuse Horizontal Irradiation—DHI, Direct Normal Irradiation—DNI, Global Horizontal Irradiation—GHI), relative cloud coverage in various atmospheric altitudes, relative humidity and wind speed. To improve prediction, two different weather forecasts were used (weatherbit.io and darksky.net), with empirically selected levels of trust, at 30% and 70%, respectively. The daily RTP scheme fed to the OSEM module was fetched by a German Utility Provider, namely Stadtwerk Haßfurt (eex.stwhas.de).

3.1.3. Experimental Scenarios

Towards providing not only proof-of-concept but also real-life experimental results, OptiMEMS was deployed on the aforementioned infrastructure over a period of nine days, during which the following scenarios were tested:
  • Baseline Operation: A fully automatic state, without any optimisation scheme applied, aiming at maximising PV generation and keeping batteries fully charged.
  • Scenario A—Optimal Day-Ahead Scheduling (“Opt Mode”): Right before midnight, OSEM creates the day-ahead optimal schedule as described in Section 2.2.1. Within the day and if the scheduling is applicable, the RT-VAOS follows that schedule and sends commands every minute to each asset with the appropriate setpoints.
  • Scenario B—Adaptive Optimal Day-Ahead Scheduling (“reOpt Mode”): As in Scenario A, OSEM creates the day-ahead schedule for the examined day. RT-VAOS monitors the status of MG DERs. When an out-of-limits deviation occurs between actual and forecasted load, or PV generation, RT-VAOS triggers a recalculation of the schedule, for the remainder of the day (see Section 2.2.2).
Towards presenting financial metrics to better assess the performance of the proposed framework, Scenarios A and B are compared to each other, using the Baseline Operation as reference. Since OptiMEMS is deployed in a real-life establishment, it should be noted that it is not possible to apply Opt and reOpt mode on the same day. Thus, two different sets of days have been chosen; in particular, Scenario A (Opt) was applied during 15–17th of July, and Scenario B (reOpt), on 19, 22–23, 25–26 and 29th. In order to be able to extract comparable results in regards to the performance of OptiMEMS, scenarios A and B were deployed in days with similar weather and operational conditions. It should be noted that during all experiments, CERTH/ITI SmartHome users were unaware of the particular days that the experiment was taking place; thus, the load profile remained intact. Therefore, it can be considered that the OptiMEMS prototype has been successfully demonstrated in an operational environment (Technology Readiness Level—TRL: 7).

3.2. Experimental Results

To comprehend how OptiMEMS works, the MG scheduling and actual operation are demonstrated for two sample summer days: Figure 4 and Figure 5 depict the MG operation under Experimental Scenario A on the 15th of July, while Figure 6 and Figure 7 show the MG operation under Scenario B on 22 July. To highlight the importance of MG operation monitoring and the adaptive nature of the optimal MG scheduling provided in Scenario B, the following interesting remark is made: the drop on the PV production curve shown with the black line in Figure 4b corresponds to the actual production, which was not accurately predicted by the PV forecasting engine in the day ahead horizon. Under Scenario A operation, the deviations between the PV forecast and actual production force the MG to bypass the originally optimal schedule calculated in the beginning of the day and turn to MG-AUTO during the slots of the deviations. This fact can be observed in Figure 5a, which shows the battery operation of that day: around the times 16:00–17:30, the battery was originally scheduled to charge (shown by the dashed black lines), however, due to the PV output drop, the OptiMEMS engine does not allow to proceed to charging so that unexpected expenses from excessive energy imports are avoided. At first glance, it is evident that in reOpt mode, the MG vDERs follow the dictated optimised schedule closer. In addition, in Figure 7a,b, the differentiation between Opt and reOpt modes is illustrated: until the first triggering at 09:30 a.m., i.e., request for a new optimised schedule (orange lines), the optimal schedule calculated at the beginning of the day (black dashed line) is identical to that actually applied (continuous grey line); however, after re-optimisation is triggered, the dispatched schedules diverge significantly. Overall, for these days, during Scenario A, the optimal scheduling was ignored by the RT-VAOS module and the MG was set to MG-AUTO mode in a total of 14–15-min slots. On the other hand, in Scenario B, the optimal schedule was not followed in only 11–15-min slots. It is interesting that the merits of reOpt mode are not entirely based upon “missing” fewer time slots but on the calculation of more accurate scheduling. This improved performance lies predominantly on the accuracy of the two short-term forecasts and is quantified using the Mean Absolute Error (MAE) of each forecasting engine (long-term and short-term forecasts for Scenario A and B, respectively) in Table 2. Indeed, the average MAE for long-term load and PV forecast are 0.268 and 0.293, respectively, while in the case of short-term forecasts, they become 0.124 and 0.107, i.e., relative reductions of MAE equal to ( 0.268 0.124 ) / 0.268 × 100 % = 53.73 % and ( 0.293 0.107 ) / 0.293 × 100 % = 63.48 % .
In order to assess the overall performance of Opt and reOpt Modes throughout the experiments, in Figure 8a, a key performance indicator (KPI) named “MGA” (%) is defined as the relative error of actually achieved to expected daily profits, with respect to baseline daily costs (9). As observed, when the MG operated in Opt mode, this KPI is always negative, meaning that it is impossible to reach the initially set goal since time slots are being “lost” when the schedule is disregarded (MG-AUTO). On the contrary, on certain days, reOpt mode may lead to positive values of the MGA KPI, meaning that the actually achieved daily profit was eventually greater than what originally predicted, a fact justified by the improved accuracy of both short-term forecasts as this derives from their re-calibration upon real-time measurements. This observation highlights the superiority of reOpt mode over the Opt mode. This statement is further strengthened by the analysis of daily cost reductions relative to the baseline daily costs (Figure 8b). Indeed, in both Opt and reOpt modes, the overall MG cost is less than what it would be without any optimisation. Specifically, Opt mode leads to approximately 5.5% daily cost reduction, whereas by employing reOpt mode, cost reduction fluctuates between 6% and 13%, with an average value of 9%. Additionally, for all Scenario B days, the daily cost reduction that would be reached if Opt mode was applied is calculated for comparison reasons. For example, during the demonstrated day with Scenario B, the daily profit was reached €0.198 (7.11% relative profit), while if Scenario A was applied on the same day, the daily profit would be equal to €0.09 (3.25% relative profit). Generally, as seen in Figure 8b, reOpt mode leads to 1.5–2 times greater relative cost reductions compared to Opt mode.
M G A = Π a Π e Π a · 100

4. Conclusions

This work presented the design and application of an optimal scheduling scheme, tailored to grid-connected, fully RES-based MG, with AC-coupled DERs. The core OSEM module is accompanied by load and PV forecasting engines, cloud services APIs and the RT-VAOS module controls of commercially available inverters. The OptiMEMS framework has been applied successfully and continuously on a real-life MG over nine summer days, and the results of its application are presented. Initial findings are very encouraging with respect to minimising daily cost during grid-connected mode and, at the same time, thanks to the developed real-time monitoring system, ensuring proper MG operation even in the case of deviations between forecasted and real-time measurements. A comparative analysis of the two possible operation modes, day-ahead-Opt mode/intra-day-reOpt mode clearly indicate the predominance of the adaptive scheduling logic. Another significant merit of the proposed methodology is the applicability of the overall framework on real world MGs due to (a) OptiMEMS coupling with commercial inverters, (b) limited infrastructure prerequisites (internet access, one computer and a Raspberry Pi) and minimal computational requirements for all developed software.
Furthermore, OptiMEMS is easily expandable, not only in terms of numbers of controllable MG DERs but also through additional functionalities, such as the inclusion of optimal vehicle-to-grid operation, a feature soon to be available by the authors as an expansion pack. Future work involves optimal load shedding capabilities for islanded operation and testing more optimisation techniques (e.g., metaheuristic techniques, or robust optimisation) for further improving the final outcome. By integrating these extra features, OptiMEMS efficiency is expected to increase, providing an even more effective tool for MG management. Finally, it is noteworthy that throughout the testing phase, it became evident that the accuracy, and thus, the overall framework effectiveness heavily depends on accurate forecasting engines. Thus, relevant research endeavours should intensify.

Author Contributions

Conceptualisation, A.D.B.; methodology, A.D.B., L.Z. and A.C.T.; formal analysis, A.D.B.; investigation, A.D.B., L.Z. and A.C.T.; experimental setup, L.Z., A.D.B. and A.C.T.; data curation, L.Z.; software, A.D.B. and L.Z.; validation, A.D.B., L.Z., A.C.T. and P.A.G.; supervision, A.C.T. and D.I.; funding acquisition, A.T., D.I. and D.T.; project administration, A.C.T. and A.T.; resources, D.T.; writing—original draft, A.D.B., L.Z., A.C.T. and P.A.G.; writing—review and editing, A.D.B., L.Z., A.C.T., P.A.G. and A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work has partially received funding from the EU H2020 R&I programme through the inteGRIDy project (grant agreement: No 731268), and the General Secretariat for Research and Innovation of Greece through the 3DMicroGrid project via the ERANETMED initiative (identification: ERANETMED-energy-11-286).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following nomenclature abbreviations are used in this manuscript:
TTotal duration of the optimisation horizon as number of t steps.
r t Time resolution considered, given the time slot duration in minutes divided by 60.
j T Auxiliary unit vector T-sized, respectively.
E S ( t ) Energy setpoints for N Energy Storage Systems (ESSs) sized 2 × N (in Wh).
E S + ( t ) , E S ( t ) Energy setpoint N-sized vectors for the charging/discharging vDERs of N ESSs (in Wh).
E S ( 0 ) The initially stored energy of the batteries (N-sized vector in Wh).
E S m i n , E S m a x ESS minimum/maximum stored energy (N-sized vectors in Wh).
E G ( t ) Energy setpoint ( [ E G ( t ) E G + ( t ) ] ( t ) ) for the import/export at MG PCC (2-sized vector in Wh).
E G m i n , E G m a x Minimum/maximum allowed energy exchange with the distribution network (T-sized vectors in W).
E P V ( t ) Energy setpoints ( E P V i ( t ) , i [ 1 , M ] ) for the M PV modules of the MG (M-sized vectors in Wh).
E P V m a x ( t ) Maximum Power Point of M PV modules (M-sized vector in Wh).
E L ( t ) MG load demand (in Wh).
β S ( t ) Index of the N ESS operation ( 2 × N ).
β S + ( t ) , β S ( t ) Operation Index (0: not operating, 1:operating) for the N charging/discharging vDERs (N-sized vectors).
β G ( t ) Binary 2-sized vector for the operation of importing/exporting ( [ β G ( t ) β G + ( t ) ] ( t ) ) vDERs.
j S , j G , j P V Auxiliary unit vectors, N/2/M-sized, respectively.
L C O E S , L C O E P V LCOE values for the N ESSs and M PV modules (in €/Wh) N/M-sized, respectively.
R C + , R C Charge/discharge C-Rates N-sized vectors, for N ESSs.
C a p Nominal N-sized vector capacity for N MG ESSs (in Wh).
D o D N-sized vector for the Depth-of-Discharge for N MG ESSs (in %).
R T P ( t ) Import ( R T P ( t ) )/export ( R T P ( t ) + ) real-time pricing vector ( T × 2 -sized vector in €/Wh).
Π e Expected profit by the microgrid throughout T optimisation horizon at the beginning of the scheduling (in €).
Π a Achieved profit by the microgrid throughout T optimisation horizon after the scheduling horizon has passed (in €).
M G A KPI “Microgrid Achievement”.

References

  1. Pulcherio, M.C.; Renjit, A.A.; Illindala, M.S.; Khalsa, A.S.; Eto, J.H.; Klapp, D.A.; Lasseter, R.H. Evaluation of Control Methods to Prevent Collapse of a Mixed-Source Microgrid. IEEE Trans. Ind. Appl. 2016, 52, 4566–4576. [Google Scholar] [CrossRef]
  2. Meng, L.; Sanseverino, E.R.; Luna, A.; Dragicevic, T.; Vasquez, J.C.; Guerrero, J.M. Microgrid supervisory controllers and energy management systems: A literature review. Renew. Sustain. Energy Rev. 2016, 60, 1263–1273. [Google Scholar] [CrossRef]
  3. Amasyali, K.; El-Gohary, N.M. A review of data-driven building energy consumption prediction studies. Renew. Sustain. Energy Rev. 2018, 81, 1192–1205. [Google Scholar] [CrossRef]
  4. Jiang, Q.; Xue, M.; Geng, G. Energy Management of Microgrid in Grid-Connected and Stand-Alone Modes. IEEE Trans. Power Syst. 2013, 28, 3380–3389. [Google Scholar] [CrossRef]
  5. Levron, Y.; Guerrero, J.M.; Beck, Y. Optimal Power Flow in Microgrids With Energy Storage. IEEE Trans. Power Syst. 2013, 28, 3226–3234. [Google Scholar] [CrossRef] [Green Version]
  6. Malysz, P.; Sirouspour, S.; Emadi, A. An Optimal Energy Storage Control Strategy for Grid-connected Microgrids. IEEE Trans. Smart Grid 2014, 5, 1785–1796. [Google Scholar] [CrossRef]
  7. Olivares, D.E.; Cañizares, C.A.; Kazerani, M. A Centralized Energy Management System for Isolated Microgrids. IEEE Trans. Smart Grid 2014, 5, 1864–1875. [Google Scholar] [CrossRef]
  8. Chalise, S.; Tonkoski, R. Day ahead schedule of remote microgrids with renewable energy sources considering battery lifetime. In Proceedings of the 11th IEEE/IAS International Conference on Industry Applications, Juiz de Fora, Brazil, 7–10 December 2014; pp. 1–5. [Google Scholar] [CrossRef]
  9. Klanšek, U. A comparison between MILP and MINLP approaches to optimal solution of Nonlinear Discrete Transportation Problem. Transport 2015, 30, 135–144. [Google Scholar] [CrossRef] [Green Version]
  10. Pegueroles-Queralt, J.; Igualada-Gonzalez, L.; Corchero-Garcia, C.; Cruz-Zambrano, M.; del Rosario-Calaf, G. Coordination of control and energy management methods for microgrid systems. In Proceedings of the IEEE PES Innovative Smart Grid Technologies Conference Europe, Instabul Turkey, 12–15 October 2014; pp. 1–6. [Google Scholar]
  11. Kanwar, A.; Rodríguez, D.I.H.; von Appen, J.; Braun, M. A Comparative Study of Optimization-and Rule-Based Control for Microgrid Operation. In Proceedings of the 2015 Power and Energy Student Summit (PESS), Dortmund, Germany, 13–14 January 2015; pp. 1–6. [Google Scholar]
  12. Graditi, G.; Silvestre, M.L.D.; Gallea, R.; Sanseverino, E.R. Heuristic-Based Shiftable Loads Optimal Management in Smart Micro-Grids. IEEE Trans. Ind. Inform. 2015, 11, 271–280. [Google Scholar] [CrossRef]
  13. Wierzbowski, M.; Lyzwa, W.; Musial, I. MILP model for long-term energy mix planning with consideration of power system reserves. Appl. Energy 2016, 169, 93–111. [Google Scholar] [CrossRef] [Green Version]
  14. Kaur, J.; Sood, Y.R.; Shrivastava, R. A two-layer optimization approach for renewable energy management of green microgrid in deregulated power sector. J. Renew. Sustain. Energy 2017, 9, 065905. [Google Scholar] [CrossRef]
  15. Choi, S.; Min, S. Optimal Scheduling and Operation of the ESS for Prosumer Market Environment in Grid-Connected Industrial Complex. IEEE Trans. Ind. Appl. 2018, 54, 1949–1957. [Google Scholar] [CrossRef]
  16. Li, Y.; Zhao, T.; Wang, P.; Gooi, H.B.; Ding, Z.; Li, K.; Yan, W. Flexible scheduling of microgrid with uncertainties considering expectation and robustness. IEEE Trans. Ind. Appl. 2017, 54, 3009–3018. [Google Scholar] [CrossRef]
  17. Paul, T.G.; Hossain, S.J.; Ghosh, S.; Mandal, P.; Kamalasadan, S. A Quadratic Programming Based Optimal Power and Battery Dispatch for Grid-Connected Microgrid. IEEE Trans. Ind. Appl. 2018, 54, 1793–1805. [Google Scholar] [CrossRef]
  18. Ji, L.; Huang, G.; Xie, Y.; Zhou, Y.; Zhou, J. Robust cost-risk tradeoff for day-ahead schedule optimization in residential microgrid system under worst-case conditional value-at-risk consideration. Energy 2018, 153, 324–337. [Google Scholar] [CrossRef]
  19. Sharma, R.; Mudaliyar, S.; Mishra, S. Power management and economic load dispatch based control of hybrid PV-battery-diesel standalone AC system. In Proceedings of the 2018 IEEMA Engineer Infinite Conference (eTechNxT), New Delhi, India, 13–14 March 2018; pp. 1–6. [Google Scholar] [CrossRef]
  20. Abdolrasol, M.G.M.; Hannan, M.A.; Mohamed, A.; Amiruldin, U.A.U.; Abidin, I.B.Z.; Uddin, M.N. An Optimal Scheduling Controller for Virtual Power Plant and Microgrid Integration Using the Binary Backtracking Search Algorithm. IEEE Trans. Ind. Appl. 2018, 54, 2834–2844. [Google Scholar] [CrossRef]
  21. Psarros, G.N.; Karamanou, E.G.; Papathanassiou, S.A. Feasibility Analysis of Centralized Storage Facilities in Isolated Grids. IEEE Trans. Sustain. Energy 2018, 9, 1822–1832. [Google Scholar] [CrossRef]
  22. Qiu, H.; Gu, W.; Xu, Y.; Zhao, B. Multi-Time-Scale Rolling Optimal Dispatch for AC/DC Hybrid Microgrids With Day-Ahead Distributionally Robust Scheduling. IEEE Trans. Softw. Eng. 2019, 10, 1653–1663. [Google Scholar] [CrossRef]
  23. Yang, F.; Feng, X.; Li, Z. Advanced Microgrid Energy Management System for Future Sustainable and Resilient Power Grid. IEEE Trans. Ind. Appl. 2019, 55, 7251–7260. [Google Scholar] [CrossRef]
  24. Riffonneau, Y.; Bacha, S.; Barruel, F.; Ploix, S. Optimal Power Flow Management for Grid Connected PV Systems With Batteries. IEEE Trans. Sustain. Energy 2011, 2, 309–320. [Google Scholar] [CrossRef]
  25. Cingoz, F.; Elrayyah, A.; Sozer, Y. Optimized Resource Management for PV-Fuel-Cell-Based Microgrids Using Load Characterizations. IEEE Trans. Ind. Appl. 2016, 52, 1723–1735. [Google Scholar] [CrossRef]
  26. Luna, A.C.; Diaz, N.L.; Graells, M.; Vasquez, J.C.; Guerrero, J.M. Mixed-Integer-Linear-Programming-Based Energy Management System for Hybrid PV-Wind-Battery Microgrids: Modeling, Design, and Experimental Verification. IEEE Trans. Power Electron. 2017, 32, 2769–2783. [Google Scholar] [CrossRef] [Green Version]
  27. Luna, A.C.; Meng, L.; Diaz, N.L.; Graells, M.; Vasquez, J.C.; Guerrero, J.M. Online Energy Management Systems for Microgrids: Experimental Validation and Assessment Framework. IEEE Trans. Power Electron. 2018, 33, 2201–2215. [Google Scholar] [CrossRef]
  28. Yi, Z.; Dong, W.; Etemadi, A.H. A Unified Control and Power Management Scheme for PV-Battery-Based Hybrid Microgrids for Both Grid-Connected and Islanded Modes. IEEE Trans. Smart Grid 2018, 9, 5975–5985. [Google Scholar] [CrossRef]
  29. Kotsampopoulos, P.; Lagos, D.; Hatziargyriou, N.; Faruque, M.O.; Lauss, G.; Nzimako, O.; Forsyth, P.; Steurer, M.; Ponci, F.; Monti, A.; et al. A Benchmark System for Hardware-in-the-Loop Testing of Distributed Energy Resources. IEEE Power Energy Technol. Syst. J. 2018, 5, 94–103. [Google Scholar] [CrossRef]
  30. Morais, H.; Kádár, P.; Faria, P.; Vale, Z.A.; Khodr, H. Optimal scheduling of a renewable micro-grid in an isolated load area using mixed-integer linear programming. Renew. Energy 2010, 35, 151–156. [Google Scholar] [CrossRef] [Green Version]
  31. Barnes, M.; Kondoh, J.; Asano, H.; Oyarzabal, J.; Ventakaramanan, G.; Lasseter, R.; Hatziargyriou, N.; Green, T. Real-world microgrids-an overview. In Proceedings of the 2007 IEEE International Conference on System of Systems Engineering, San Antonio, TX, USA, 16–18 April 2007; pp. 1–8. [Google Scholar]
  32. Silvente, J.; Graells, M.; Espuña, A.; Salas, P. An optimization model for the management of energy supply and demand in smart grids. In Proceedings of the 2012 IEEE International Energy Conference and Exhibition (ENERGYCON), Florence, Italy, 9–12 September 2012; pp. 368–373. [Google Scholar] [CrossRef]
  33. Pell, S.; Turcotte, D.; Colgate, G.; Swingler, A. Nemiah Valley Photovoltaic-Diesel Mini-Grid: System Performance and Fuel Saving Based on One Year of Monitored Data. IEEE Trans. Sustain. Energy 2012, 3, 167–175. [Google Scholar] [CrossRef]
  34. Aluisio, B.; Dicorato, M.; Forte, G.; Trovato, M. An optimization procedure for MG day-ahead operation in the presence of CHP facilities. Sustain. Energy Grids Netw. 2017, 11, 34–45. [Google Scholar] [CrossRef]
  35. Lai, K.; Illindala, M.S.; Haj-ahmed, M.A. Comprehensive protection strategy for an islanded microgrid using intelligent relays. In Proceedings of the 2015 IEEE Industry Applications Society Annual Meeting, Addison, TX, USA, 18–22 October 2015; pp. 1–11. [Google Scholar] [CrossRef]
  36. Marzband, M.; Sumper, A.; Ruiz-Álvarez, A.; Domínguez-García, J.L.; Tomoiagă, B. Experimental evaluation of a real time energy management system for stand-alone microgrids in day-ahead markets. Appl. Energy 2013, 106, 365–376. [Google Scholar] [CrossRef]
  37. Bao, Z.; Zhou, Q.; Yang, Z.; Yang, Q.; Xu, L.; Wu, T. A Multi Time-Scale & Multi Energy-Type Coordinated Microgrid Scheduling Solution—Part I: Model & Methodology. IEEE Trans. Power Syst. 2015, 30, 2257–2266. [Google Scholar] [CrossRef]
  38. Nguyen, D.T.; Le, L.B. Optimal Bidding Strategy for Microgrids Considering Renewable Energy and Building Thermal Dynamics. IEEE Trans. Smart Grid 2014, 5, 1608–1620. [Google Scholar] [CrossRef]
  39. Ebrahimi, M.R.; Amjady, N. Adaptive robust optimization framework for day-ahead microgrid scheduling. Int. J. Electr. Power Energy Syst. 2019, 107, 213–223. [Google Scholar] [CrossRef]
  40. Conte, F.; Massucco, S.; Saviozzi, M.; Silvestro, F. A Stochastic Optimization Method for Planning & Real-Time Control of Integrated PV-Storage Systems: Design & Experimental Validation. IEEE Trans. Softw. Eng. 2018, 9, 1188–1197. [Google Scholar] [CrossRef]
  41. Oriti, G.; Julian, A.L.; Anglani, N.; Hernandez, G.D. Novel Economic Analysis to Design the Energy Storage Control System of a Remote Islanded Microgrid. IEEE Trans. Ind. Appl. 2018, 54, 6332–6342. [Google Scholar] [CrossRef]
  42. Fathima, A.H.; Palanisamy, K. Optimization in microgrids with hybrid energy systems—A review. Renew. Sustain. Energy Rev. 2015, 45, 431–446. [Google Scholar] [CrossRef]
  43. Dolara, A.; Grimaccia, F.; Magistrati, G.; Marchegiani, G. Optimization Models for Islanded Micro-Grids: A Comparative Analysis between Linear Programming and Mixed Integer Programming. Energies 2017, 10, 241. [Google Scholar] [CrossRef]
  44. Vrettos, E.I.; Papathanassiou, S.A. Operating Policy and Optimal Sizing of a High Penetration RES-BESS System for Small Isolated Grids. IEEE Trans. Energy Convers. 2011, 26, 744–756. [Google Scholar] [CrossRef]
  45. Liao, G. The optimal economic dispatch of smart Microgrid including Distributed Generation. In Proceedings of the 2013 International Symposium on Next-Generation Electronics, Kaohsiung, Taiwan, 25–26 February 2013; pp. 473–477. [Google Scholar]
  46. Abbes, D.; Martinez, A.; Champenois, G. Eco-design optimisation of an autonomous hybrid wind-photovoltaic system with battery storage. IET Renew. Power Gener. 2012, 6, 358–371. [Google Scholar] [CrossRef]
  47. Conti, S.; Nicolosi, R.; Rizzo, S.A.; Zeineldin, H.H. Optimal Dispatching of Distributed Generators and Storage Systems for MV Islanded Microgrids. IEEE Trans. Power Deliv. 2012, 27, 1243–1251. [Google Scholar] [CrossRef]
  48. Bao, Z.; Zhou, Q.; Yang, Z.; Yang, Q.; Xu, L.; Wu, T. A Multi Time-Scale & Multi Energy-Type Coordinated Microgrid Scheduling Solution—Part II: Optimization Algorithm & Case Studies. IEEE Trans. Power Syst. 2015, 30, 2267–2277. [Google Scholar] [CrossRef]
  49. Marzband, M.; Yousefnejad, E.; Sumper, A.; Domínguez-García, J.L. Real time experimental implementation of optimum energy management system in standalone Microgrid by using multi-layer ant colony optimization. Int. J. Electr. Power Energy Syst. 2016, 75, 265–274. [Google Scholar] [CrossRef] [Green Version]
  50. Mirjalili, S. Evolutionary Algorithms and Neural Networks: Theory and Applications; Springer: Cham, Switzerland, 2019; pp. 87–104. [Google Scholar] [CrossRef]
  51. Anglani, N.; Oriti, G.; Colombini, M. Optimized Energy Management System to Reduce Fuel Consumption in Remote Military Microgrids. IEEE Trans. Ind. Appl. 2017, 53, 5777–5785. [Google Scholar] [CrossRef]
  52. Conte, F.; D’Agostino, F.; Pongiglione, P.; Saviozzi, M.; Silvestro, F. Mixed-Integer Algorithm for Optimal Dispatch of Integrated PV-Storage Systems. IEEE Trans. Ind. Appl. 2019, 55, 238–247. [Google Scholar] [CrossRef]
  53. AhmadiAhangar, R.; Rosin, A.; Niaki, A.N.; Palu, I.; Korõtko, T. A review on real-time simulation and analysis methods of microgrids. Int. Trans. Electr. Energy Syst. 2019, 29, e12106. [Google Scholar] [CrossRef] [Green Version]
  54. Digital Innovation Hubs. Smart Specialisation Platform. Available online: https://s3platform.jrc.ec.europa.eu/digital-innovation-hubs-tool (accessed on 7 September 2020).
  55. Van Ackooij, W.; Danti Lopez, I.; Frangioni, A.; Lacalandra, F.; Tahanan, M. Large-scale unit commitment under uncertainty: An updated literature survey. Ann. Oper. Res. 2018, 271, 11–85. [Google Scholar] [CrossRef] [Green Version]
  56. Lotfi, H.; Khodaei, A. Levelized cost of energy calculations for microgrids. In Proceedings of the 2016 IEEE Power and Energy Society General Meeting (PESGM), Boston, MA, USA, 17–21 July 2016; pp. 1–5. [Google Scholar] [CrossRef]
  57. Mitchell, S.; Consulting, S.M.; Dunning, I. PuLP: A Linear Programming Toolkit for Python. Available online: http://www.optimization-online.org/DB_FILE/2011/09/3178.pdf (accessed on 1 September 2020).
  58. Forrest, J.; Lougee-Heimer, R. COIN-OR Branch and Cut Solver. Available online: https://projects.coin-or.org/Cbc (accessed on 1 September 2020).
  59. Chen, T.; Guestrin, C. XGBoost: A Scalable Tree Boosting System. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016; pp. 785–794. [Google Scholar] [CrossRef] [Green Version]
  60. Xue, P.; Jiang, Y.; Zhou, Z.; Chen, X.; Fang, X.; Liu, J. Multi-step ahead forecasting of heat load in district heating systems using machine learning algorithms. Energy 2019, 188, 116085. [Google Scholar] [CrossRef]
  61. Timplalexis, C.; Bezas, N.; Bintoudi, D.A.; Zyglakis, L.; Pavlopoulos, V.; Tsolakis, C.A.; Krinidis, S.; Tzovaras, D. A hybrid physical/statistical day-ahead direct PV forecasting engine. In Proceedings of the Mediterranean Conference on Power Generation, Transmission, Distribution and Energy Conversion (MEDPOWER 2020), Pafos, Cyprus, 9–11 November 2020. [Google Scholar]
  62. Salamanis, A.I.; Xanthopoulou, G.; Bezas, N.; Timplalexis, C.; Bintoudi, A.D.; Zyglakis, L.; Tsolakis, A.C.; Ioannidis, D.; Kehagias, D.; Tzovaras, D. Benchmark Comparison of Analytical, Data-Based and Hybrid Models for Multi-Step Short-Term Photovoltaic Power Generation Forecasting. Energies 2020, 13, 5978. [Google Scholar] [CrossRef]
  63. Holmgren, W.F.; Hansen, C.W.; Mikofski, M.A. pvlib python: A python package for modeling solar energy systems. J. Open Source Softw. 2018, 3, 884. [Google Scholar] [CrossRef] [Green Version]
  64. Tsolakis, A.C.; Bintoudi, A.D.; Zyglakis, L.; Zikos, S.; Timplalexis, C.; Bezas, N.; Kitsikoudis, K.; Ioannidis, D.; Tzovaras, D. Design and Real-life Deployment of a Smart Nanogrid: A Greek Case Study. In Proceedings of the 2020 IEEE International Conference on Power and Energy (PECon), Penang, Malaysia, 7–8 December 2020; pp. 321–326. [Google Scholar]
Figure 1. System Architecture of the OptiMEMS Framework.
Figure 1. System Architecture of the OptiMEMS Framework.
Energies 14 02752 g001
Figure 2. Logical flow chart of the RT-VAOS module.
Figure 2. Logical flow chart of the RT-VAOS module.
Energies 14 02752 g002
Figure 3. SmartHome microgrid infrastructure.
Figure 3. SmartHome microgrid infrastructure.
Energies 14 02752 g003
Figure 4. Scenario A—Optimisation mode on 15/7: consumption/generation.
Figure 4. Scenario A—Optimisation mode on 15/7: consumption/generation.
Energies 14 02752 g004
Figure 5. Scenario A—Optimisation mode on 15/7: BESS/grid operation.
Figure 5. Scenario A—Optimisation mode on 15/7: BESS/grid operation.
Energies 14 02752 g005
Figure 6. Scenario B—Re-optimisation mode on 22/7: consumption/generation.
Figure 6. Scenario B—Re-optimisation mode on 22/7: consumption/generation.
Energies 14 02752 g006
Figure 7. Scenario B—Re-Optimisation mode on 22/7: BESS/grid operation.
Figure 7. Scenario B—Re-Optimisation mode on 22/7: BESS/grid operation.
Energies 14 02752 g007
Figure 8. Comparative performance graphs during scenarios A (grey) and B (black).
Figure 8. Comparative performance graphs during scenarios A (grey) and B (black).
Energies 14 02752 g008
Table 2. MAE of load and PV forecasting tools (kW per hour quarter “kWh-q”).
Table 2. MAE of load and PV forecasting tools (kW per hour quarter “kWh-q”).
Opt Dates15/716/717/7 Mean Values
MAE Load0.1920.2320.379 0.268
MAE PV0.0920.2650.523 0.293
reOpt Dates19/722/723/725/726/729/7Mean Values
MAE Load0.1530.1260.1240.1140.1310.0970.124
MAE PV0.1110.0490.1390.1220.1250.1010.107
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bintoudi, A.D.; Zyglakis, L.; Tsolakis, A.C.; Gkaidatzis, P.A.; Tryferidis, A.; Ioannidis, D.; Tzovaras, D. OptiMEMS: An Adaptive Lightweight Optimal Microgrid Energy Management System Based on the Novel Virtual Distributed Energy Resources in Real-Life Demonstration. Energies 2021, 14, 2752. https://doi.org/10.3390/en14102752

AMA Style

Bintoudi AD, Zyglakis L, Tsolakis AC, Gkaidatzis PA, Tryferidis A, Ioannidis D, Tzovaras D. OptiMEMS: An Adaptive Lightweight Optimal Microgrid Energy Management System Based on the Novel Virtual Distributed Energy Resources in Real-Life Demonstration. Energies. 2021; 14(10):2752. https://doi.org/10.3390/en14102752

Chicago/Turabian Style

Bintoudi, Angelina D., Lampros Zyglakis, Apostolos C. Tsolakis, Paschalis A. Gkaidatzis, Athanasios Tryferidis, Dimosthenis Ioannidis, and Dimitrios Tzovaras. 2021. "OptiMEMS: An Adaptive Lightweight Optimal Microgrid Energy Management System Based on the Novel Virtual Distributed Energy Resources in Real-Life Demonstration" Energies 14, no. 10: 2752. https://doi.org/10.3390/en14102752

APA Style

Bintoudi, A. D., Zyglakis, L., Tsolakis, A. C., Gkaidatzis, P. A., Tryferidis, A., Ioannidis, D., & Tzovaras, D. (2021). OptiMEMS: An Adaptive Lightweight Optimal Microgrid Energy Management System Based on the Novel Virtual Distributed Energy Resources in Real-Life Demonstration. Energies, 14(10), 2752. https://doi.org/10.3390/en14102752

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop