P. Palensky
Please Note
51 records found
1
Cyber Security of Digital Substations
Machine Learning-based Deep Packet Inspection for Detecting Cyber Attacks on IEC 61850 GOOSE
This work proposes a Hardware-in-the-Loop (HiL) test setup to execute various GOOSE cyber attacks and thereby simulate a hacker's actions. This setup consists of a simple power system simulated on a Real-Time Digital Simulator (RTDS), a physical Intelligent Electronic Device (IED), and a communication network connecting all. The simulated power system communicates node voltages and breaker currents to the IED via IEC 61850 Sample Values (SV), and the IED responds by sending GOOSE traffic. An additional workstation is connected to the communication network to launch cyber attacks that cause physical impact on the simulated power system.
Secondly, the HiL setup is used to evaluate which alterations to the GOOSE packet will result in a physical impact on the simulated power system. Several attributes in the GOOSE PDU are modified, and together with changes in AllData (for circuit breaker tripping), the circuit breaker in the simulated power system should be tripped. An attempt is also made to block legitimate traffic during a fault, with an attack. Based on these findings, a cyber-physical dataset was constructed containing GOOSE communication network traffic recorded during normal operation, faults, and the examined cyber attacks that yielded physical impact.
Furthermore, an anomaly-based deep packet inspection (DPI) intrusion detection system (IDS) is proposed for the mitigation of cyber attacks. This DPI-IDS uses features from the GOOSE PDU attributes and a long short-term memory (LSTM) model to distinguish GOOSE packets from normal operation, faults, and cyber attacks. The LSTM's hyperparameters were optimized, and the complete DPI model was trained on primarily GOOSE traffic from normal operation and fault conditions. The performance of the DPI-IDS on the collected dataset was evaluated using several metrics. For all attacks in the dataset, the performance is evaluated separately to identify which attacks the DPI-IDS model performs best for.
The goal of the DPI-IDS is to classify legitimate traffic from malicious traffic. Normal operation traffic and traffic during faults should be classified correctly as legitimate traffic. Correct classification of malicious traffic would cause the traffic to be flagged, indicating to an operator to take action. Overall, the results of legitimate traffic identification (normal operation and faults) show that the DPI-IDS performs well on separating these two classes. However, the classification of malicious traffic is more difficult, due to the limited availability of malicious traffic in the training data. This underlines the importance of developing an effective mitigation strategy for cyber attacks on GOOSE communication traffic. ...
This work proposes a Hardware-in-the-Loop (HiL) test setup to execute various GOOSE cyber attacks and thereby simulate a hacker's actions. This setup consists of a simple power system simulated on a Real-Time Digital Simulator (RTDS), a physical Intelligent Electronic Device (IED), and a communication network connecting all. The simulated power system communicates node voltages and breaker currents to the IED via IEC 61850 Sample Values (SV), and the IED responds by sending GOOSE traffic. An additional workstation is connected to the communication network to launch cyber attacks that cause physical impact on the simulated power system.
Secondly, the HiL setup is used to evaluate which alterations to the GOOSE packet will result in a physical impact on the simulated power system. Several attributes in the GOOSE PDU are modified, and together with changes in AllData (for circuit breaker tripping), the circuit breaker in the simulated power system should be tripped. An attempt is also made to block legitimate traffic during a fault, with an attack. Based on these findings, a cyber-physical dataset was constructed containing GOOSE communication network traffic recorded during normal operation, faults, and the examined cyber attacks that yielded physical impact.
Furthermore, an anomaly-based deep packet inspection (DPI) intrusion detection system (IDS) is proposed for the mitigation of cyber attacks. This DPI-IDS uses features from the GOOSE PDU attributes and a long short-term memory (LSTM) model to distinguish GOOSE packets from normal operation, faults, and cyber attacks. The LSTM's hyperparameters were optimized, and the complete DPI model was trained on primarily GOOSE traffic from normal operation and fault conditions. The performance of the DPI-IDS on the collected dataset was evaluated using several metrics. For all attacks in the dataset, the performance is evaluated separately to identify which attacks the DPI-IDS model performs best for.
The goal of the DPI-IDS is to classify legitimate traffic from malicious traffic. Normal operation traffic and traffic during faults should be classified correctly as legitimate traffic. Correct classification of malicious traffic would cause the traffic to be flagged, indicating to an operator to take action. Overall, the results of legitimate traffic identification (normal operation and faults) show that the DPI-IDS performs well on separating these two classes. However, the classification of malicious traffic is more difficult, due to the limited availability of malicious traffic in the training data. This underlines the importance of developing an effective mitigation strategy for cyber attacks on GOOSE communication traffic.
MMC Control in HVDC System connected to Offshore wind farms
Control Strategy for asymmetric offshore AC faults
Current TSA methods mainly rely on passive clustering, focusing on the statistical proximity between the representative periods and the full-space time-series. However, this approach does not guarantee a satisfactory solution for the final expansion planning problem in general cases, even with predefined extreme periods implemented such as days with maximum load and minimum available RES. The operation of the power system and extreme conditions are highly sensitive to the specific power system configurations, making the standalone TSA method unreliable in practical applications.
To improve the time-series aggregation in terms of the power system operation, firstly, the methods of assessing the performance of the estimated investment decision are explored directly in terms of the objective function. By introducing the full-space operational cost model for the power system with the estimated investment decision, the operational cost error made by representative periods can be obtained, referred to as the operational estimation error. The actual objective difference between the estimated investment decision and the optimal investment decision found by the full-space expansion planning model can also be evaluated, denoted as the optimality gap. It is found that the optimality gap is bounded by the difference in the operational estimation error of representative periods for the power system with the two investment decisions. As the operations of the power system with the estimated and optimal investment decisions are better estimated, the simplified model can provide a closer investment decision.
Subsequently, looking into the operational estimation error, the performance of representative periods in estimating the full-space operational cost is highly unevenly distributed among the full-space time-series. Representative periods fail to accurately estimate a minor portion of the full-space time-series, causing extremely high operational estimation errors that contribute to the majority of the total operational estimation error. This uneven distribution remains as the number of representative periods increases, indicating that standalone TSA methods are not capable of capturing the extreme conditions of the specific power system operation.
Therefore, to improve the time-series aggregation in terms of its ability to better estimate the operation of the power system, bad-performing representative periods and original periods with high operational estimation error can be identified and prioritized, forming a feedback enhancement loop. Case studies show that the feedback enhancement with re-clustering on the bad-performing representative periods improves the optimality gap by more than 50% compared with the standard mean-based clustering method. ...
Current TSA methods mainly rely on passive clustering, focusing on the statistical proximity between the representative periods and the full-space time-series. However, this approach does not guarantee a satisfactory solution for the final expansion planning problem in general cases, even with predefined extreme periods implemented such as days with maximum load and minimum available RES. The operation of the power system and extreme conditions are highly sensitive to the specific power system configurations, making the standalone TSA method unreliable in practical applications.
To improve the time-series aggregation in terms of the power system operation, firstly, the methods of assessing the performance of the estimated investment decision are explored directly in terms of the objective function. By introducing the full-space operational cost model for the power system with the estimated investment decision, the operational cost error made by representative periods can be obtained, referred to as the operational estimation error. The actual objective difference between the estimated investment decision and the optimal investment decision found by the full-space expansion planning model can also be evaluated, denoted as the optimality gap. It is found that the optimality gap is bounded by the difference in the operational estimation error of representative periods for the power system with the two investment decisions. As the operations of the power system with the estimated and optimal investment decisions are better estimated, the simplified model can provide a closer investment decision.
Subsequently, looking into the operational estimation error, the performance of representative periods in estimating the full-space operational cost is highly unevenly distributed among the full-space time-series. Representative periods fail to accurately estimate a minor portion of the full-space time-series, causing extremely high operational estimation errors that contribute to the majority of the total operational estimation error. This uneven distribution remains as the number of representative periods increases, indicating that standalone TSA methods are not capable of capturing the extreme conditions of the specific power system operation.
Therefore, to improve the time-series aggregation in terms of its ability to better estimate the operation of the power system, bad-performing representative periods and original periods with high operational estimation error can be identified and prioritized, forming a feedback enhancement loop. Case studies show that the feedback enhancement with re-clustering on the bad-performing representative periods improves the optimality gap by more than 50% compared with the standard mean-based clustering method.
This thesis proposes using a black box approach and machine learning models to predict the hysteresis factor for SOC correction more accurately. The research aims to enhance SOC accuracy in real-world EV environments, beyond laboratory settings. Improved SOC estimation through machine learning could significantly advance battery management systems (BMS), yielding both economic and environmental benefits by increasing EV efficiency. The study addresses the existing gaps in SOC estimation for voltage hysteresis, considering the vehicle as a system by incorporating driving cycles, BMS technical limitations, and cell chemistry. By improving EV reliability and efficiency, this research aims to promote broader acceptance and contribute to a sustainable future in personal transportation.
...
This thesis proposes using a black box approach and machine learning models to predict the hysteresis factor for SOC correction more accurately. The research aims to enhance SOC accuracy in real-world EV environments, beyond laboratory settings. Improved SOC estimation through machine learning could significantly advance battery management systems (BMS), yielding both economic and environmental benefits by increasing EV efficiency. The study addresses the existing gaps in SOC estimation for voltage hysteresis, considering the vehicle as a system by incorporating driving cycles, BMS technical limitations, and cell chemistry. By improving EV reliability and efficiency, this research aims to promote broader acceptance and contribute to a sustainable future in personal transportation.
Low cost power quality measuring unit for household usage and small to enterprise scale installations
Designing a Low Cost Power Quality Analyzer
Demand response in a container terminal
A stochastic optimization of the operational planning considering energy consumption
Facilitating flexibility trading between asset owners and system operators
Creating a protocol for flexibility exchange between the grid operator and flexible assets
Defense Against Malicious Parameter Identification
System Immersion Coding and Hybrid Multiplicative Watermarking
The system immersion coding method, derived from the privacy solution in federated learning, is adapted to enhance its capability to detect malicious attacks by merging the input information and defend malicious parameter identification by leveraging its privacy-preserving properties. This method involves mapping the plant output into a higher-dimensional space and introducing carefully defined noise, which can create arbitrarily large disturbances without compromising performance. The introduced disturbance disrupts the attacker's parameter estimation. Theoretical conditions are provided to discuss the detection performance of replay attacks, control-signal-injection zero-dynamics attacks, and sensor-signal-injection zero-dynamics attacks. However, we also identify that the system immersion coding method is vulnerable to known-plaintext attacks.
Watermarking is a promising active diagnosis technique for the detection of highly sophisticated attacks. Motivated by the computational hardness problems of cryptography analysis, we propose a hybrid multiplicative watermarking scheme as an active diagnosis technique. In this scheme, watermarking parameters are periodically updated based on the dynamics of unobservable states in specifically designed piecewise affine (PWA) hybrid systems. We conduct a theoretical analysis to assess the impact of this scheme on closed-loop performance, demonstrating its stability preservation. We also provide conditions to detect replay attacks and control-signal-injection zero-dynamics attacks. Furthermore, we demonstrate that the proposed approach makes it challenging for an eavesdropper to reconstruct watermarking parameters, considering both computational complexity and systems theoretic perspectives. ...
The system immersion coding method, derived from the privacy solution in federated learning, is adapted to enhance its capability to detect malicious attacks by merging the input information and defend malicious parameter identification by leveraging its privacy-preserving properties. This method involves mapping the plant output into a higher-dimensional space and introducing carefully defined noise, which can create arbitrarily large disturbances without compromising performance. The introduced disturbance disrupts the attacker's parameter estimation. Theoretical conditions are provided to discuss the detection performance of replay attacks, control-signal-injection zero-dynamics attacks, and sensor-signal-injection zero-dynamics attacks. However, we also identify that the system immersion coding method is vulnerable to known-plaintext attacks.
Watermarking is a promising active diagnosis technique for the detection of highly sophisticated attacks. Motivated by the computational hardness problems of cryptography analysis, we propose a hybrid multiplicative watermarking scheme as an active diagnosis technique. In this scheme, watermarking parameters are periodically updated based on the dynamics of unobservable states in specifically designed piecewise affine (PWA) hybrid systems. We conduct a theoretical analysis to assess the impact of this scheme on closed-loop performance, demonstrating its stability preservation. We also provide conditions to detect replay attacks and control-signal-injection zero-dynamics attacks. Furthermore, we demonstrate that the proposed approach makes it challenging for an eavesdropper to reconstruct watermarking parameters, considering both computational complexity and systems theoretic perspectives.
In this work, a blockchain-based application for the monitoring and control of a feeder in the Low-Voltage (LV) distribution grid is designed and tested. A smart contract is created and deployed in a private Ethereum blockchain utilizing the Proof of Authority (PoA) consensus mechanism. The blockchain application enhances the cyber security of the LV distribution system in three ways. First, it detects cyber attacks targeting DERs by comparing the setpoints received by prosumers with smart meter measurements. Second, it prevents cyber attacks by enabling the exchange of measurements and setpoints on-chain and by preventing unreliable prosumers from participating in the voltage regulation market. Third, it mitigates the effects of cyber attacks on the steady-state voltage magnitudes by enforcing a novel voltage regulation mechanism, in which a new metric is proposed to quantify the power-to-voltage relationship while considering the location of the power exchange.
The efficacy of the blockchain application is tested in a co-simulation environment together with a modeled LV distribution network, simulated in DigSILENT PowerFactory. The distribution network model is first used to assess the impact of cyber attacks manipulating the setpoints of Battery Energy Storage Systems (BESSs), which have been identified as the most critical DERs. The simulation results demonstrate that the considered cyber attacks can force the disconnection of inverters by causing violations of the acceptable steady-state voltage magnitudes. One of the scenarios demonstrates that a cyber attack targeting half of the BESSs in a feeder can lead to the collapse of the voltage, causing a local outage. Finally, the results of the co-simulation of the blockchain-based monitoring and control system, achieved by the Open Platform Communications Unified Architecture (OPC UA) communication protocol and by a series of clients managing the data streams, demonstrate its efficacy in detecting cyber attacks and mitigating their impact on the voltage magnitude across the feeder, thus reducing the number of disconnected DERs.
...
In this work, a blockchain-based application for the monitoring and control of a feeder in the Low-Voltage (LV) distribution grid is designed and tested. A smart contract is created and deployed in a private Ethereum blockchain utilizing the Proof of Authority (PoA) consensus mechanism. The blockchain application enhances the cyber security of the LV distribution system in three ways. First, it detects cyber attacks targeting DERs by comparing the setpoints received by prosumers with smart meter measurements. Second, it prevents cyber attacks by enabling the exchange of measurements and setpoints on-chain and by preventing unreliable prosumers from participating in the voltage regulation market. Third, it mitigates the effects of cyber attacks on the steady-state voltage magnitudes by enforcing a novel voltage regulation mechanism, in which a new metric is proposed to quantify the power-to-voltage relationship while considering the location of the power exchange.
The efficacy of the blockchain application is tested in a co-simulation environment together with a modeled LV distribution network, simulated in DigSILENT PowerFactory. The distribution network model is first used to assess the impact of cyber attacks manipulating the setpoints of Battery Energy Storage Systems (BESSs), which have been identified as the most critical DERs. The simulation results demonstrate that the considered cyber attacks can force the disconnection of inverters by causing violations of the acceptable steady-state voltage magnitudes. One of the scenarios demonstrates that a cyber attack targeting half of the BESSs in a feeder can lead to the collapse of the voltage, causing a local outage. Finally, the results of the co-simulation of the blockchain-based monitoring and control system, achieved by the Open Platform Communications Unified Architecture (OPC UA) communication protocol and by a series of clients managing the data streams, demonstrate its efficacy in detecting cyber attacks and mitigating their impact on the voltage magnitude across the feeder, thus reducing the number of disconnected DERs.
Security of supply during the energy transition
Evaluating the adequacy of the Dutch energy system during the transition to a zero-carbon energy system in a realistic scenario
The way the current Dutch electricity market is designed will lead to a drop in electricity prices on the wholesale market, when the market is penetrated with a high concentration of renewable energy sources. So the electricity market mechanism needs to be reformed to have a price determination system better suited for energy sources with low marginal costs. Furthermore, a capacity remuneration system needs to be implemented to create enough investment incentive in controllable energy capacity to maintain a high level of security of supply in the Dutch energy system during the energy transition.
There exist many different forms of capacity remuneration mechanisms, but the strategic reserve and capacity markets are examined to check which one would work best in the Dutch energy system. In a strategic reserve a central authority sets the capacity volume, which will be contracted by an operator. This operator will dispatch the contracted capacity when needed. In a capacity market a central authority determines the amount of capacity each consumer should acquire through buying capacity credits on the market.
To determine which capacity remuneration mechanism would deliver the highest level of security of supply, simulations of the Dutch energy system are run. The simulations were done using the models AMIRIS and EMLab, which are coupled to give a more accurate representation of reality. The coupling takes place in the Spine Toolbox. The capacity remuneration mechanisms were created in Python and can be activated as modules in the Spine Toolbox.
The vertices are sketched of the possible future energy scenario of the Netherlands in 2050. These vertices are the so called scenarios, regional governance, national governance, European CO2-governance and international governance. The simulations were run with a scenario in which the energy demand stays equal and with a scenario where the energy demand descends similar to the international governance scenario. The international governance scenario was chosen to base the simulations on, because the simulations make profit driven decisions, similar to the international governance scenario.
The results of the simulations show that, firstly, without any governmental interference, a (near) zero-carbon energy system will not be achieved in 2050. Secondly, when running the energy-only market in a realistic scenario, no shortage hours will occur and the system had sufficient generation capacity for supplying the demand, even in a "dunkelflaute" scenario. Thirdly, the basic version of the strategic reserve offered the least costly solution for society for providing sufficient security of supply, so no significant shortage periods will occur in a realistic scenario. Finally, the best option for maintaining security of supply in the Netherlands during the transition to a zero-carbon energy system is the implementation of the yearly capacity market. ...
The way the current Dutch electricity market is designed will lead to a drop in electricity prices on the wholesale market, when the market is penetrated with a high concentration of renewable energy sources. So the electricity market mechanism needs to be reformed to have a price determination system better suited for energy sources with low marginal costs. Furthermore, a capacity remuneration system needs to be implemented to create enough investment incentive in controllable energy capacity to maintain a high level of security of supply in the Dutch energy system during the energy transition.
There exist many different forms of capacity remuneration mechanisms, but the strategic reserve and capacity markets are examined to check which one would work best in the Dutch energy system. In a strategic reserve a central authority sets the capacity volume, which will be contracted by an operator. This operator will dispatch the contracted capacity when needed. In a capacity market a central authority determines the amount of capacity each consumer should acquire through buying capacity credits on the market.
To determine which capacity remuneration mechanism would deliver the highest level of security of supply, simulations of the Dutch energy system are run. The simulations were done using the models AMIRIS and EMLab, which are coupled to give a more accurate representation of reality. The coupling takes place in the Spine Toolbox. The capacity remuneration mechanisms were created in Python and can be activated as modules in the Spine Toolbox.
The vertices are sketched of the possible future energy scenario of the Netherlands in 2050. These vertices are the so called scenarios, regional governance, national governance, European CO2-governance and international governance. The simulations were run with a scenario in which the energy demand stays equal and with a scenario where the energy demand descends similar to the international governance scenario. The international governance scenario was chosen to base the simulations on, because the simulations make profit driven decisions, similar to the international governance scenario.
The results of the simulations show that, firstly, without any governmental interference, a (near) zero-carbon energy system will not be achieved in 2050. Secondly, when running the energy-only market in a realistic scenario, no shortage hours will occur and the system had sufficient generation capacity for supplying the demand, even in a "dunkelflaute" scenario. Thirdly, the basic version of the strategic reserve offered the least costly solution for society for providing sufficient security of supply, so no significant shortage periods will occur in a realistic scenario. Finally, the best option for maintaining security of supply in the Netherlands during the transition to a zero-carbon energy system is the implementation of the yearly capacity market.
The objective of this thesis is to explore the potential of modern signal processing and machine learning (ML) techniques for the development of an explainable classification scheme of faults in the distribution network. ML explainability is a huge part of the research held for this thesis, as it is very important for the fault analysis department to understand the model and for the grid operators to trust the results of the model. The problem is divided into two parts. The first part concerns the construction of an explainable ML-based classification model that can accurately differentiate between types of stable faults. The second part concerns the identification of suitable criteria for single-phase and multi-phase fault stability. Combined with a fault classifier based on the first part research, this will lead to a classification and location scheme that can incorporate all fault types.
First, a literature study is done. The pre-modeling explainability stage is the most suitable one and based on that, the Short-time Fourier Transform (STFT) is chosen as the signal processing technique since it can lead to more explainable features than other techniques, e.g. the Wavelet Transform (WT). Also, the Support Vector Machines (SVMs) are selected as the supervised learning technique, as it performs well in classification problems and it is not a method as complicated as for example the Artificial Neural Network (ANN). Next, three different feature sets are shaped appropriately from the three-phase fault current and voltage signals, which are the inputs for the SVM classifier model, so then this model can be tuned and evaluated. The model with the symmetrical feature set is selected serving suitably both the performance and explainability requirements. Following that, a set of stability rules are developed to better identify the stable faults that previously were classified as unstable. These rules are tuned based on a set of initially labeled unstable faults. In the end, the ML classification and the stability model combined with section selection rules constitute the final classification scheme. This is tested on hidden data and the results indicate good performance in classifying faults, locating faults and identifying stable faults that initially would be classified as unstable. In particular, the performance of the proposed scheme on that test set has f1-score of 95.1% exceeding 93.3% of the current algorithm tested on the same data, errors regarding the fault loop-impedance/reactance are located within the acceptable limits and 4 stable faults correctly identified that previously were seen as unstable. ...
The objective of this thesis is to explore the potential of modern signal processing and machine learning (ML) techniques for the development of an explainable classification scheme of faults in the distribution network. ML explainability is a huge part of the research held for this thesis, as it is very important for the fault analysis department to understand the model and for the grid operators to trust the results of the model. The problem is divided into two parts. The first part concerns the construction of an explainable ML-based classification model that can accurately differentiate between types of stable faults. The second part concerns the identification of suitable criteria for single-phase and multi-phase fault stability. Combined with a fault classifier based on the first part research, this will lead to a classification and location scheme that can incorporate all fault types.
First, a literature study is done. The pre-modeling explainability stage is the most suitable one and based on that, the Short-time Fourier Transform (STFT) is chosen as the signal processing technique since it can lead to more explainable features than other techniques, e.g. the Wavelet Transform (WT). Also, the Support Vector Machines (SVMs) are selected as the supervised learning technique, as it performs well in classification problems and it is not a method as complicated as for example the Artificial Neural Network (ANN). Next, three different feature sets are shaped appropriately from the three-phase fault current and voltage signals, which are the inputs for the SVM classifier model, so then this model can be tuned and evaluated. The model with the symmetrical feature set is selected serving suitably both the performance and explainability requirements. Following that, a set of stability rules are developed to better identify the stable faults that previously were classified as unstable. These rules are tuned based on a set of initially labeled unstable faults. In the end, the ML classification and the stability model combined with section selection rules constitute the final classification scheme. This is tested on hidden data and the results indicate good performance in classifying faults, locating faults and identifying stable faults that initially would be classified as unstable. In particular, the performance of the proposed scheme on that test set has f1-score of 95.1% exceeding 93.3% of the current algorithm tested on the same data, errors regarding the fault loop-impedance/reactance are located within the acceptable limits and 4 stable faults correctly identified that previously were seen as unstable.
Modelling Hydrogen in Power Systems
Optimisation for Investment and Operational Models
The hydrogen pathway as aforementioned is divided over three technologies: hydrogen generation with means of water electrolysis, also known as 'green hydrogen', storage in compression vessels and reconversion of hydrogen into electricity in the form of a fuel cell technology (also known as Power-to-Gas).
The research focuses on identifying technical parameters and operational policies of the water electrolysis systems that can be translated into optimisation constraints, assessing the level of detail required to create an accurate optimisation model. A generic model is developed that can be scaled for further research, making different case studies and sizing possible. The research compares the performance of the models in terms of accuracy to the computational burden. The comparison is done for the level of detail and complexity added to the model.
After a literature review of technical parameters and operational policies regarding the technologies, two models were created in a mathematical framework. The two models proposed were Linear Programming (LP) and a Mixed-Integer Programming (MIP) Model. On the LP model 6 different sensitivity analysis has been performed, to be precise on Capital Expenditures (CAPEX), efficiency, lifetime, ramping rates, interest rates and finally different time horizons. The outcome of these analyses is that the technology mix can best be used in a combined manner, whereby each component of the mix contributes towards minimising the objective value: the Total Annualised Cost.
Lastly the two models are compared with different types of configurations, each with a different set of constraints. The constraints to be modelled were: minimum uptime and downtime, start-up costs, degradation due to cycling and finally the part-load operation. ...
The hydrogen pathway as aforementioned is divided over three technologies: hydrogen generation with means of water electrolysis, also known as 'green hydrogen', storage in compression vessels and reconversion of hydrogen into electricity in the form of a fuel cell technology (also known as Power-to-Gas).
The research focuses on identifying technical parameters and operational policies of the water electrolysis systems that can be translated into optimisation constraints, assessing the level of detail required to create an accurate optimisation model. A generic model is developed that can be scaled for further research, making different case studies and sizing possible. The research compares the performance of the models in terms of accuracy to the computational burden. The comparison is done for the level of detail and complexity added to the model.
After a literature review of technical parameters and operational policies regarding the technologies, two models were created in a mathematical framework. The two models proposed were Linear Programming (LP) and a Mixed-Integer Programming (MIP) Model. On the LP model 6 different sensitivity analysis has been performed, to be precise on Capital Expenditures (CAPEX), efficiency, lifetime, ramping rates, interest rates and finally different time horizons. The outcome of these analyses is that the technology mix can best be used in a combined manner, whereby each component of the mix contributes towards minimising the objective value: the Total Annualised Cost.
Lastly the two models are compared with different types of configurations, each with a different set of constraints. The constraints to be modelled were: minimum uptime and downtime, start-up costs, degradation due to cycling and finally the part-load operation.
the longer-term energy storage systems are not fully represented since, for existing storage systems, the self-serving nature of these leads to participation in multiple energy markets. This is because participating in other markets, like the balancing markets, can lead to higher overall profits than a storage system only participating in the wholesale market.
This thesis investigates different energy storage technologies and multiple prominent storage applications for grids. Furthermore, an overview of the European energy markets will be examined, and different design options will be discussed. These markets include frequency containment reserve (FCR), frequency regulation reserves (aFRR/mFRR) and the wholesale markets. The review of storage technologies, applications, and available markets has led to the development and simulation of single-purpose energy storage models fulfilling grid applications.
By combining the specific purpose models, a complete energy market and energy storage model representation could be created. The model created is unique since the complete energy system model allows energy storage systems to optimally dispatch over multiple markets while at the same time also influencing these markets. Multiple cases were investigated using this model, such as the influence of increasing storage capacity on the wholesale and balancing market and the influence of storage systems just performing one service, so only regulation, arbitrage or peak-shaving. Based on the model results, recommendations are made on improving the current energy market designs and how to better represent storage systems in existing energy system models. ...
the longer-term energy storage systems are not fully represented since, for existing storage systems, the self-serving nature of these leads to participation in multiple energy markets. This is because participating in other markets, like the balancing markets, can lead to higher overall profits than a storage system only participating in the wholesale market.
This thesis investigates different energy storage technologies and multiple prominent storage applications for grids. Furthermore, an overview of the European energy markets will be examined, and different design options will be discussed. These markets include frequency containment reserve (FCR), frequency regulation reserves (aFRR/mFRR) and the wholesale markets. The review of storage technologies, applications, and available markets has led to the development and simulation of single-purpose energy storage models fulfilling grid applications.
By combining the specific purpose models, a complete energy market and energy storage model representation could be created. The model created is unique since the complete energy system model allows energy storage systems to optimally dispatch over multiple markets while at the same time also influencing these markets. Multiple cases were investigated using this model, such as the influence of increasing storage capacity on the wholesale and balancing market and the influence of storage systems just performing one service, so only regulation, arbitrage or peak-shaving. Based on the model results, recommendations are made on improving the current energy market designs and how to better represent storage systems in existing energy system models.
The LEPP to 2050 is made possible by normalising the selected markets. As a result, the conditions unprecedented in Dutch data are covered in the normalised and combined data set. Generally, under scenarios with high proportions of hydrogen power in the energy portfolio, the cWGAN-GP model projects that DAM price series would have slightly lower mean and daily standard deviation than the 2019 level. Whereas much lower mean and daily standard deviation are projected when natural gas is still the fuel of the most frequent final generating technology. To explore the possible application of the projector model, the German DAM prices series in 2019 have been projected and evaluated, and the projections under Dutch 2050 energy scenarios have been used in calculating the generic profit potential of energy storage.
Five findings can be summarised from the main results. Firstly, from a literature survey and importance analyses, seven features are shown relevant to the DAM price in the combined data set, namely month of the year, day of the week, total hourly load forecast, national daily mean temperature, fuel cost of the most frequent final generating technology, hourly renewable power generation forecast and total installed renewable power capacity. Secondly, it has been found that two of the four proposed market state normalisation solutions, the Renewable Scarcity Factor (RSF) and the Renewable-Load Ratio (RLR) help the cWGAN-GP model strike a balance between price value distribution and hourly inter-dependencies. Thirdly, in this LEPP study, the cWGAN-GP model performs better than the Conditional Variational Auto-Encoder (CVAE) and multivariate Gaussian distribution (mGaus) models. Compared with the two alternatives, the cWGAN-GP model produces samples in better quality while remaining sensitive to temporal conditions. Fourthly, projections by the cWGAN-GP model are more realistic than those made by the Energy Transition Model (ETM), with price values varying continuously in smooth boundaries. Finally, the fuel cost of the most frequent final generating technology is found critical to LEPP. The annual mean and daily standard deviation of the DAM price series are expected to rise significantly when natural gas is mostly replaced by hydrogen power in the national energy portfolio. ...
The LEPP to 2050 is made possible by normalising the selected markets. As a result, the conditions unprecedented in Dutch data are covered in the normalised and combined data set. Generally, under scenarios with high proportions of hydrogen power in the energy portfolio, the cWGAN-GP model projects that DAM price series would have slightly lower mean and daily standard deviation than the 2019 level. Whereas much lower mean and daily standard deviation are projected when natural gas is still the fuel of the most frequent final generating technology. To explore the possible application of the projector model, the German DAM prices series in 2019 have been projected and evaluated, and the projections under Dutch 2050 energy scenarios have been used in calculating the generic profit potential of energy storage.
Five findings can be summarised from the main results. Firstly, from a literature survey and importance analyses, seven features are shown relevant to the DAM price in the combined data set, namely month of the year, day of the week, total hourly load forecast, national daily mean temperature, fuel cost of the most frequent final generating technology, hourly renewable power generation forecast and total installed renewable power capacity. Secondly, it has been found that two of the four proposed market state normalisation solutions, the Renewable Scarcity Factor (RSF) and the Renewable-Load Ratio (RLR) help the cWGAN-GP model strike a balance between price value distribution and hourly inter-dependencies. Thirdly, in this LEPP study, the cWGAN-GP model performs better than the Conditional Variational Auto-Encoder (CVAE) and multivariate Gaussian distribution (mGaus) models. Compared with the two alternatives, the cWGAN-GP model produces samples in better quality while remaining sensitive to temporal conditions. Fourthly, projections by the cWGAN-GP model are more realistic than those made by the Energy Transition Model (ETM), with price values varying continuously in smooth boundaries. Finally, the fuel cost of the most frequent final generating technology is found critical to LEPP. The annual mean and daily standard deviation of the DAM price series are expected to rise significantly when natural gas is mostly replaced by hydrogen power in the national energy portfolio.
Congestion forecasting using a custom loss function
Application to congestion mitigation on substation Middelharnis
In this thesis, the use of a a custom loss function is explored, which incorporates the different costs associated with congestion mitigation. This function aims to improve upon existing congestion forecasts by having a loss function in line with the congestion forecasts goal. For this research, the case of substation Middelharnis is used. This substation encounters congestion due to connected wind parks and distributed photovoltaics. The congestion encountered will be solved by grid expansion in 2024, but until then will have to be mitigated using the redispatching market GOPACS.
A custom function was generated using two cost components: a congestion fee, which needs to be paid if wind parks are disconnected and GOPACS costs, which are the price of buying flexible power on the redispatch market. The function takes into account two-time horizons: costs associated with buying GOPACS power day-ahead acting on a prediction of the load, and the resulting cost on the day itself if any remaining congestion was resolved via the congestion fee. The contribution of this function is that it depends on the difference between the costs of the prediction and the cost of the realization, instead of the cost of the difference between prediction and realization.
Four models were trained, each with a different loss function: MSE, Pinball, and Cost. The cost model is trained twice and for each model, a different value for GOPACS cost is used. The results show that the cost models outperform the Pinball and MSE-trained models, when the performance is evaluated on the final cost metric. However, if the GOPACS price is high, the cost model becomes conservative in predicting congestion peaks. It is not expected that this will be an issue in the future, as with more renewable energy, the energy price and subsequently the GOPACS price will be lower. ...
In this thesis, the use of a a custom loss function is explored, which incorporates the different costs associated with congestion mitigation. This function aims to improve upon existing congestion forecasts by having a loss function in line with the congestion forecasts goal. For this research, the case of substation Middelharnis is used. This substation encounters congestion due to connected wind parks and distributed photovoltaics. The congestion encountered will be solved by grid expansion in 2024, but until then will have to be mitigated using the redispatching market GOPACS.
A custom function was generated using two cost components: a congestion fee, which needs to be paid if wind parks are disconnected and GOPACS costs, which are the price of buying flexible power on the redispatch market. The function takes into account two-time horizons: costs associated with buying GOPACS power day-ahead acting on a prediction of the load, and the resulting cost on the day itself if any remaining congestion was resolved via the congestion fee. The contribution of this function is that it depends on the difference between the costs of the prediction and the cost of the realization, instead of the cost of the difference between prediction and realization.
Four models were trained, each with a different loss function: MSE, Pinball, and Cost. The cost model is trained twice and for each model, a different value for GOPACS cost is used. The results show that the cost models outperform the Pinball and MSE-trained models, when the performance is evaluated on the final cost metric. However, if the GOPACS price is high, the cost model becomes conservative in predicting congestion peaks. It is not expected that this will be an issue in the future, as with more renewable energy, the energy price and subsequently the GOPACS price will be lower.
The integration between these three energy sectors is assumed to be chiefly driven by the operation of electrolysers and fuel cells. Two other main assumptions set the base for this work: first, the gas sector is assumed to be entirely repurposed to operate with hydrogen and, secondly, the waste heat coming from fuel cells and electrolysers is assumed to be the main thermal energy input of a district heating network. An electricity and a hydrogen market have also been modelled to simulate the interaction of this region with the external grid.
The analysis is carried out by means of a linear optimization algorithm coded using oemof, an open source python package for multi-energy system modelling and optimization. Most of the input data (i.e. energy demand and generation profiles) comes from the Integrale Ifrastructuurverkenning 2030-2050 study from TenneT, Gasunie and the Dutch DSOs, which has been regionalized for the Drechsteden region is order to optimize the investments on the main energy assets (transformers, hydrogen substations, electrolysers, fuel cells, batteries) needed to run the future energy system.
An optimal system configuration is also calculated for a "reinforcement" scenario, in which energy sectors remain independent, and a comparison with the system integration scenario is presented. The overall system costs appear to be only about 5% lower for the system integration scenario. The costs imputable to assets is higher with system integration due to the addition of expensive hydrogen conversion assets, but, in this scenario, the operation on the energy market (driven chiefly by hydrogen exports) is more advantageous.
Finally, a sensitivity analysis on the share of heat demand to be satisfied by district heating and on the price of batteries is carried out, along with an investigation on the effect of adding thermal storage to the system. ...
The integration between these three energy sectors is assumed to be chiefly driven by the operation of electrolysers and fuel cells. Two other main assumptions set the base for this work: first, the gas sector is assumed to be entirely repurposed to operate with hydrogen and, secondly, the waste heat coming from fuel cells and electrolysers is assumed to be the main thermal energy input of a district heating network. An electricity and a hydrogen market have also been modelled to simulate the interaction of this region with the external grid.
The analysis is carried out by means of a linear optimization algorithm coded using oemof, an open source python package for multi-energy system modelling and optimization. Most of the input data (i.e. energy demand and generation profiles) comes from the Integrale Ifrastructuurverkenning 2030-2050 study from TenneT, Gasunie and the Dutch DSOs, which has been regionalized for the Drechsteden region is order to optimize the investments on the main energy assets (transformers, hydrogen substations, electrolysers, fuel cells, batteries) needed to run the future energy system.
An optimal system configuration is also calculated for a "reinforcement" scenario, in which energy sectors remain independent, and a comparison with the system integration scenario is presented. The overall system costs appear to be only about 5% lower for the system integration scenario. The costs imputable to assets is higher with system integration due to the addition of expensive hydrogen conversion assets, but, in this scenario, the operation on the energy market (driven chiefly by hydrogen exports) is more advantageous.
Finally, a sensitivity analysis on the share of heat demand to be satisfied by district heating and on the price of batteries is carried out, along with an investigation on the effect of adding thermal storage to the system.
Due to dependence on threshold values that vary across the literature and the limitations associated with calculating thresholds as a percentage of installed capacity, it was demonstrated that binary ramp definitions are not ideal and result in under-reporting. On the other hand, the wavelet approach extracts ramp events from the generation using statistically determined threshold values. As a result, the problem of under-detection of ramp events is mitigated. The proposed approach of "significant ramps" allows the evaluation of which ramp events are important and which are far less disruptive and may be ignored.
It was demonstrated that anti-correlation between wind and solar resources alone is not adequate to promise a smoother output as it does not provide sufficient information about ramp events. Anti-correlations at shorter time resolutions, such as 15 minutes or an hour, could be preferable. While seasonal anti-correlation may benefit national system adequacy, it does not benefit daily ramping events.
The optimal wind-PV capacity size for decreasing the total number of ramps was such that wind turbines filled the grid capacity, as solar power would result in extra ramps. It was observed that solar over-planting leads to a significantly increased number of ramp events, whereas wind over-planting results in a minimal change in ramp events. A penalty price was proposed to internalize the severity of ramp events, which could influence the choice between wind and solar over-planting. A solution was presented to mitigate ramp incidents in a hybrid power plant using a battery which was found to be more effective and/or more economical in minimizing ramps compared to over-planting.
The proposed "significant wavelet ramp approach" is shown to be a useful metric for characterizing wind-based hybrid power plant ramp occurrences. For a future in which variable renewable energy sources account for a substantial portion of the energy mix, it is proposed that demand information be considered when defining ramp events. More attention must be paid to power ramp occurrences, either by penalizing ramps or enforcing tougher grid codes. The ramp events must be considered at the sizing and development stage, with the possibility of including a ramp-mitigating battery strategy. A thorough examination of ramp events in hybrid power plants demonstrates the importance of minimizing and managing ramp events for both the system operator and the producer.
...
Due to dependence on threshold values that vary across the literature and the limitations associated with calculating thresholds as a percentage of installed capacity, it was demonstrated that binary ramp definitions are not ideal and result in under-reporting. On the other hand, the wavelet approach extracts ramp events from the generation using statistically determined threshold values. As a result, the problem of under-detection of ramp events is mitigated. The proposed approach of "significant ramps" allows the evaluation of which ramp events are important and which are far less disruptive and may be ignored.
It was demonstrated that anti-correlation between wind and solar resources alone is not adequate to promise a smoother output as it does not provide sufficient information about ramp events. Anti-correlations at shorter time resolutions, such as 15 minutes or an hour, could be preferable. While seasonal anti-correlation may benefit national system adequacy, it does not benefit daily ramping events.
The optimal wind-PV capacity size for decreasing the total number of ramps was such that wind turbines filled the grid capacity, as solar power would result in extra ramps. It was observed that solar over-planting leads to a significantly increased number of ramp events, whereas wind over-planting results in a minimal change in ramp events. A penalty price was proposed to internalize the severity of ramp events, which could influence the choice between wind and solar over-planting. A solution was presented to mitigate ramp incidents in a hybrid power plant using a battery which was found to be more effective and/or more economical in minimizing ramps compared to over-planting.
The proposed "significant wavelet ramp approach" is shown to be a useful metric for characterizing wind-based hybrid power plant ramp occurrences. For a future in which variable renewable energy sources account for a substantial portion of the energy mix, it is proposed that demand information be considered when defining ramp events. More attention must be paid to power ramp occurrences, either by penalizing ramps or enforcing tougher grid codes. The ramp events must be considered at the sizing and development stage, with the possibility of including a ramp-mitigating battery strategy. A thorough examination of ramp events in hybrid power plants demonstrates the importance of minimizing and managing ramp events for both the system operator and the producer.