Next Article in Journal
Artificial Intelligence-Assisted RFID Tag-Integrated Multi-Sensor for Quality Assessment and Sensing
Previous Article in Journal
Relative Motion Estimation Algorithm for Noncooperative Targets Considering Multiple Solutions of Rotational Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multi-Level Operation Method for Improving the Resilience of Power Systems under Extreme Weather through Preventive Control and a Virtual Oscillator

1
Electric Power Research Institute of State Grid Henan Electric Power Company, Zhengzhou 450000, China
2
College of New Energy, Harbin Institute of Technology at Weihai, Weihai 264209, China
3
College of Lilac, Harbin Institute of Technology at Weihai, Weihai 264209, China
4
School of Electrical Engineering and Automation, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(6), 1812; https://doi.org/10.3390/s24061812
Submission received: 23 January 2024 / Revised: 24 February 2024 / Accepted: 6 March 2024 / Published: 12 March 2024
(This article belongs to the Section Electronic Sensors)

Abstract

:
This paper proposes a multi-level operation method designed to enhance the resilience of power systems under extreme weather conditions by utilizing preventive control and virtual oscillator (VO) technology. Firstly, a novel model for predicting time intervals between successive failures of the power system during extreme weather is introduced. Based on this, this paper proposes a preventive control method considering the system ramping and transmission constraints prior to failures so as to ensure the normal electricity demand within the system. Further, a VO-based adaptive frequency control strategy is designed to accelerate the regulation speed and eliminate the frequency deviation. Finally, the control performance is comprehensively compared under different experimental conditions. The results verify that the method accurately predicted the time of the line fault occurrence, with a maximum error not exceeding 3 min compared to the actual occurrence; also, the virtual oscillator control (VOC) strategy outperformed traditional droop control in frequency stabilization, achieving stability within 2 s compared to the droop control’s continued fluctuations beyond 20 s. These results highlight VOC’s superior effectiveness in frequency stability and control in power systems.

1. Introduction

In recent years, the escalating severity of climate change has led to an increased frequency of extreme weather events, critically impacting the operational stability of power systems [1]. For example, Texas experienced cold weather conditions in 2021, leading to the most significant load shedding in the state’s history. Forty states in the US have reported longer power outages due to destructive storms, and temperatures reaching up to 40 degrees in some southern provinces of China have increased the risk to power systems due to increased electricity demand [2]. In light of the growing challenges posed by extreme weather, the study of power system resilience has become increasingly critical. Resilience refers to the system’s ability to maintain operational functionality and rapidly restore system performance when facing low-probability extreme events [3]. Furthermore, confronted with challenges to power systems’ resilience, it is important to develop effective methods for assessing the risk of failures and implement measures to enhance grid stability and resilience.
Numerous studies have been conducted to enhance the resilience of power systems under extreme weather conditions. Some studies have focused on evaluations of system resilience. These evaluations employ various quantification techniques to identify the strengths and weaknesses of power systems and propose resilience enhancement strategies. For instance, one study introduced a scoring matrix to evaluate system functions from different perspectives [4], while another proposed the resilience “trapezoid” to capture critical system degradation and recovery features across different event phases [5]. Additionally, resilience has been quantified as the difference between the fully functioning system’s capacity and its post-event capacity, offering a more nuanced understanding of resilience impacts [6].
In the realm of operational and planning strategies to enhance resilience, a lot of approaches have been explored. Studies have proposed probabilistic proactive generation redispatch strategies [7] and optimized recovery strategies considering distribution system reconfiguration [8]. A mixed-integer linear programming (MILP)-based method for generation redispatch during ice storms aims to strengthen resilience [9], while a three-stage decision-making methodology evaluates weighted pre-event and post-event performance losses, addressing the spatial uncertainty in fault probabilities [10]. This is complemented by innovative applications of battery energy storage systems (BESSs) for earthquake resilience [11] and the use of machine learning for fault prediction and resilience enhancement, like a distributed offline–online architecture [12] and a multiagent-based hybrid soft actor–critic (HSAC) algorithm [13]. Preventive-action-based strategies, such as multi-sensor prediction-based wide-area monitoring and control [14] and Monte Carlo simulation-based proactive unit commitment frameworks [15], are also gaining prominence in the resilience enhancement discourse.
Furthermore, the integration of renewable energy sources (RESs) has been identified as a crucial factor in augmenting power system resilience. Studies have demonstrated the efficacy of microgrid islanding in creating more resilient distribution systems [16]. A multistage restoration approach that utilizes heuristic optimization algorithms for system reconfiguration and DG islanding has been proposed [17], emphasizing the potential of RESs in emergency response scenarios. This includes approaches like transmission line switching [18], load redispatching [19], and resilience-based unit commitment [20], which collectively underscore the growing importance of RESs in enhancing the resilience of power systems against extreme weather impacts.
The aforementioned methods are crucial for enhancing the resilience of a power system, yet they have limitations. On the one hand, under extreme weather conditions, the power system may experience consecutive failures, potentially leading to insufficient local supply capabilities or even the emergence of localized power supply islands [21,22,23]. In such scenarios, the combined ramping capability of the generators within the local supply area and the power transmission capacity from the external grid may not be sufficient to meet the normal electricity demand within the area. However, current research primarily focuses on emergency control after power system failures, with little consideration given to the impact of potential consecutive failures before they occur. This oversight is particularly evident in the formulation of emergency dispatch strategies for power systems, where there is a lack of comprehensive integration of system ramping and transmission constraints. The system ramping constraint is critical, as it delineates the rate at which a power generator can regulate its output to match fluctuations in electricity demand, given the physical limitation on the rapidity of power output adjustment. Similarly, the transmission constraint encapsulates the capacity limitation of the power transmission network, encompassing a high-voltage transmission line and other grid components, which are pivotal in ensuring the safe and efficient transportation of electrical power without causing overheating or unacceptable voltage drops. It is important that an emergency dispatch strategy accounts for these two constraints to ensure the operational feasibility and resilience of a power system against potential failures.
On the other hand, extreme weather can lead to fragile, isolated local power supply areas with insufficient resources, poor inertia, and low disturbance resistance, threatening frequency stability. Advanced frequency regulation methods are required for their stability. The increasing presence of power electronic interfaces in power grids has enhanced regulation flexibility. However, traditional controls like droop and virtual synchronous machine controls in power electronics, such as energy storage units, are slow to respond to disturbances, particularly in isolated power areas. Recently, virtual oscillator control (VOC) has become important for frequency regulation in power systems rich in power electronics [24,25]. VOC, using the digital control of converters, simulates LC second-order oscillation dynamics. Research on VOC aims to improve inverter performance and practical VOC applications, such as power scheduling in grid-connected modes, islanding support, and enhanced disturbance management. A challenge is balancing VOC’s preference for larger inertia for angular frequency stability against the need for smaller inertia for a better active power response. Thus, a flexible, adjustable VOC beneficial for the frequency response is preferred. Additionally, developing a VOC-based controller that minimizes frequency errors is necessary.
To address the aforementioned issues, this paper proposes a multi-level operation method for improving the resilience of power systems under extreme weather through preventive control and a VO. The new contributions can be summarized as follows:
  • This paper introduces a novel method of preventive control for power systems, particularly focusing on the power constraints between successive failures. This method is significant, as it considers both the combined ramping capability of generators within the local supply area and the power transmission capacity from the external grid before the occurrence of successive failures. This approach is designed to ensure that the normal electricity demand within the local supply area is met, even in the face of potential disruptions. Such a preventive control mechanism is innovative in its comprehensive approach to power system resilience, especially in scenarios of extreme weather and high-risk conditions for cascading failures.
  • The paper proposes a VOC-based adaptive frequency control strategy for power systems. A key element of innovation here is the dynamic adjustment of the VOC loop coefficients according to the system’s operating conditions during transient states. This adaptive strategy marks a significant advancement over traditional VOC implementations, where coefficients are usually static. By dynamically adjusting the VOC loop coefficients, this method enhances the frequency recovery rate, offering a more responsive and efficient approach to frequency regulation under varying system conditions.
The rest of this paper is organized as follows. In Section 2, an upper-level preventive control method that considers successive failures is presented. Section 3 presents a lower-level control method based on a VO to accelerate the frequency regulation speed. A case study is provided in Section 4. Section 5 concludes the paper.

2. Upper-Level Preventive Control Considering Successive Failures

Extreme weather can potentially lead to the successive tripping of multiple power lines, and the whole power system might be partitioned into relatively independent regions and even several isolated regions. In this circumstance, some regions may not have an abundant power supply due to the ramping limits of generators and the limited power transmission capacity of lines. This upper-level control framework aims to ensure that the combined capability of the ramping up of generators within the local supply region and the power transmission capacity from the external grid are greater than the load within the local supply region.
During the fault process, the transmission system needs to redispatch generators after the occurrence of faults in order to meet the security conditions of the system under successive failures. The objective function is firstly established according to the dynamic situation of the power system, as shown in (1). The objective function aims to minimize the cost of electricity generation as much as possible while also ensuring minimal variation in the system’s states (voltage and phase angle) before and after any faults.
min C = m = 1 N G c G m P G m , t + m = 1 N E c E m P E m , t + z t h ( x t )
where P G m , t is the active power of the m-th generator at time t. c G m is the cost coefficient of the m-th generator. N G is the number of generators. P E m , t is the active power of the m-th energy storage unit at time t. c E m is the cost coefficient of the m-th energy storage unit. N E is the number of energy storage units. z t represents the measurement vector including the bus voltage magnitude V i , t ( i Ω B , and Ω B is the set of system buses) and phase angle deviation θ i j ( i Ω B , j Ω B ) before the occurrence of the fault. h ( x t ) is the fault state vector including the bus voltage magnitude V i , t ( i Ω B , and Ω B is the set of system buses) and phase angle deviation θ i j ( i Ω B , j Ω B ) after the occurrence of the fault.
Furthermore, in the method of preventive control, it is necessary to consider the following constraints:
(1)
Power flow constraint: The power system must satisfy the power flow equations, and the line flow should be within the prescribed limits [26].
P G i , t all P D i , t all = j Ω B V i , t V j , t ( G i j cos θ i j , t + B i j sin θ i j , t ) Q G i , t all Q D i , t all = j Ω B V i , t V j , t ( G i j sin θ i j , t B i j cos θ i j , t )
P i j , t = ( V i , t ) 2 G i j + V i , t V j , t ( G i j cos θ i j , t + B i j sin θ i j , t ) Q i j , t = ( V i , t ) 2 B i j + V i , t V j , t ( G i j sin θ i j , t B i j cos θ i j , t ) P i j , min P i j , t P i j , max Q i j , min Q i j , t Q i j , max
where P G i , t all and Q G i , t all are the active and reactive power injected into bus i of all kinds of generation systems, including generators, energy storage units, wind power systems, etc. P D i , t all and Q D i , t all are the active and reactive power demands at bus i of all kinds of loads. P i j , t and Q i j , t are the active and reactive power from the i-th bus to the j-th bus. G i j and B i j represent the admittance parameter between transmission line ij. P i j , max and P i j , min represent the active power limits of transmission line ij. Q i j , max and Q i j , min denote the reactive power limits of transmission line ij.
(2)
System reserve constraint: The system needs positive and negative spinning reserve capacity to mitigate the impact of uncertainties in wind power and load forecasting, thus ensuring the security and stability of the system operation.
i = 1 N G R G i , t up + i = 1 N E R E i , t up ε L P Load , t + ε W P fW , t i = 1 N G R G i , t down + i = 1 N E R E i , t down ε L P Load , t + ε W P fW , t
where ε L and ε W represent the proportional forecasting errors of the system load and wind power, respectively. R G i , t up and R G i , t down denote the positive/negative spinning reserve capacity of generator i at time t. R E i , t up and R E i , t down denote the positive/negative spinning reserve capacity of energy storage unit i at time t. P Load , t denotes the total load within the system, while P fW , t denotes the forecasted wind power output.
(3)
Constraints on the operation of the generator: These encompass the constraints on the output, the ramping constraints, the startup and shutdown times, and the reserve constraints of generators.
P G i , min P G i , t P G i , max
T G i , t on T G i , min on T G i , t off T G i , min off P G i , t + 1 P G i , t v Ramp , G i up Δ t P G i , t P G i , t + 1 v Ramp , G i down Δ t
P G i , t + R G i , t up P G i , max 0 R G i , t up v Ramp , G i up Δ t P G i , t R G i , t down P G i , min 0 R G i , t down v Ramp , G i down Δ t
where P G i , max and P G i , min , respectively, denote the upper and lower limits of the active power output of generator i. T G i , t on represents the startup time of generator i during time period t, while T G i , t off represents the shutdown time of generator i during time period t. T G i , min on denotes the minimum allowable startup time for generator i, and T G i , t off represents the minimum allowable shutdown time for generator i. v Ramp , G i up and v Ramp , G i down represent the upward and downward ramp rates of generator i, respectively. Δ t denotes the minimum scheduling period.
(4)
Constraints on the operation of the energy storage unit: These encompass constraints on the reserve and capacity of energy storage units.
0 R E i , t up P N , E i P E i , t 0 R E i , t down P E i , t + P N , E i
E min , i E t , i ( P E i , t + R E i , t up ) / η i E max , i E min , i E t , i ( P E i , t R E i , t down ) / η i E max , i
where P N , E i represents the rated power of energy storage unit i, while E min , i and E max , i represent the minimum and maximum energy capacities of storage unit i. η i represents the charge or discharge efficiency of storage unit i.
(5)
Power constraint between successive failures: Following a power system failure, isolated operation in certain regions or limited interconnection with the rest of the large power grid may ensue. In such scenarios, the affected area may experience insufficient generation capacity. In particular, in the intervals between consecutive failures, without a well-orchestrated generation plan, generators within the region may struggle to swiftly adjust their output to anticipate the next contingency. Additionally, due to transmission constraints, the remainder of the large power grid may be unable to supply adequate power to this area, thereby increasing the risk of widespread power outages in the power system.
We utilize Figure 1 for a more intuitive explanation. Assuming that Node 9 bears a substantial load at time t 0 while Node 3 exhibits a limited generation level at the same time, at this juncture, the load at Node 9 will be sustained by Lines 3–9, 6–9, and 8–9. Furthermore, employing the approach outlined in Section 2, if it is possible to forecast the occurrence of faults on Lines 6–9 and 8–9, and the time of occurrence for each fault is denoted by t 1 and t 1 + K Δ t , respectively, it becomes evident that, following the fault on Line 6–9, the load at Node 9 will be supported by Lines 3–9 and 8–9. Subsequently, after the fault on Line 8–9, the load at Node 9 can only be supplied by Line 3–9. If the ramping capability of the generator at Node 3 is insufficient to meet the load requirement at Node 9 at time t 1 + K Δ t , it becomes necessary for Node 9 to curtail a portion of its load to maintain supply–demand equilibrium. One potential resolution is to preemptively increase the output of the generator at Node 3 upon forecasting the fault occurrences on Lines 6–9 and 8–9 at times t 1 and t 1 + K Δ t , thereby mitigating Node 9’s reliance on Lines 6–9 and 8–9 for supply and ensuring that the generator’s output at Node 3 is able to ramp up to meet the load requirements at times t 1 and t 1 + K Δ t .
Based on the aforementioned analysis, the power constraint between successive failures is set as follows:
i Ω ISO ( P G i , t + R G i , t up ) + i j Ω CON P i j , max i Ω ISO P D i , t all i Ω ISO ( P G i , t + 1 + R G i , t + 1 up ) + i j Ω CON P i j , max i Ω ISO P D i , t + 1 all i Ω ISO ( P G i , t + K + R G i , t + K up ) + i j Ω CON P i j , max i Ω ISO P D i , t + K all
where K represents the time interval between adjacent failures. Ω ISO denotes the set of nodes in the local power supply area, while Ω CON signifies the set of lines connected to the local power supply area.
The equation illustrates that, following a failure, the power system may be partitioned into relatively independent regions. It is imperative to ensure that, during the occurrence of a failure, the load within the local power supply area remains unaffected. Specifically, the combined ramping capability of the generators within the local supply area and the power transmission capacity from the external grid must be sufficient to meet the normal electricity demand within the area.
(6)
Constraint on renewable energy output: The actual dispatch output of wind power must not exceed the forecasted wind power output.
0 P sW , t P fW , t
where P sW , t represents the actual output of wind power.
To sum up, the power flow and generator/energy storage system operation constraints are critical for ensuring that the physical operation of the power system adheres to technical specifications and limits. The constraints of the system reserve and power stability between failures are vital for preparing the system against uncertainties (like wind power and load forecasting errors) and for maintaining stability and continuous operation in the aftermath of system failures, ensuring that there is sufficient generation capacity and strategic operation planning to avoid widespread outages. The constraint on renewable energy output focuses on aligning the actual renewable power dispatch with forecasted outputs. By combining the above six constraints and the objective function, the upper-level preventive control scheme considering successive failures can be obtained.
Notably, in the preventive control optimization process, we utilize a differential evolution–particle swarm optimization (DE-PSO) algorithm [27]. The DE-PSO algorithm, merging differential evolution (DE) and particle swarm optimization (PSO), excels in optimizing complex problems without easily getting trapped in local optima. Its hybrid nature combines DE’s global search with PSO’s swift convergence, enabling the effective exploration and exploitation of the search space. This makes DE-PSO particularly adept at handling problems with numerous constraints, as it can adaptively balance between exploring new areas and refining existing solutions. Thus, DE-PSO stands out for its ability to navigate challenging optimization landscapes and find high-quality solutions amidst a multitude of constraints. Considering the dynamic and often rapidly changing conditions of power systems during extreme weather events, the optimization is set to recur at intervals of 15 min. This interval is optimal for updating the system’s preventive strategies in response to evolving environmental and operational conditions, thereby enhancing the system’s adaptability and resilience against cascading failures. Testing indicates that the PSO algorithm can complete an optimization cycle in approximately 10 s, and this duration is acceptable and ensures a balance between computational expediency and precision.

3. Lower-Level Control Based on Virtual Oscillator

On short time scales, the stochastic nature of the load and renewable energy sources may influence the stability of the system frequency, especially when the power system is divided into relatively independent regions due to faults. This section (Section 3) explores the rapid adjustment characteristics of power electronic devices of energy storage devices and proposes a frequency control method for power systems based on a VO.
Firstly, the integration of VOC into a multi-level frequency control architecture marks a significant advancement for power systems, particularly those rich in power electronics. This approach represents a pioneering step in harmonizing VOC strategies with contemporary power system control requirements. Secondly, we developed an innovative frequency regulation method using the VOC strategy. This method utilizes a rapid time-domain model in conjunction with coupling droop characteristics, enabling VOC-controlled inverters to measure frequency without the need for a Phase-Locked Loop (PLL). This advancement significantly reduces the computational load, enhancing the efficiency of the system. Thirdly, this manuscript details a novel methodology for optimizing VOC parameters, including virtual capacitance and inductance. This optimization is crucial for dynamically adjusting system parameters in response to real-time operational conditions, thereby improving the frequency regulation performance of power systems [28].

3.1. Characteristic Analysis of VOC Inverter

Figure 2 depicts a typical VOC inverter of the Van der Pol (VdP) type, comprising a VdP type of VO, DG, and an inverter. The circuit model of the VdP type of VO consists of three parallel parts: (1) an LC resonance loop that can generate a resonant frequency and whose oscillation characteristics are reflected in the bidirectional energy flow between a virtual capacitor and a virtual inductor; (2) a negative impedance, which is essential to the system’s transient process and closely related to the inertial frequency regulation scale; (3) a voltage-controlled current source.
The dynamics of VOC are given by [29]
d i L d t = v c L d v c d t = α C v c 3 + σ C v c 1 C i L κ i C i
where i L is the inductor current. v c is the capacitor voltage. C , L , and σ are the virtual capacitance, inductance, and conductance, and the nonlinear voltage-dependent current source is parameterized by α . The controller is interfaced via voltage and current scaling values, κ v and κ i .
From Figure 2, the output voltage command is
v = κ v v C cos Φ ε i L sin Φ
where ε = L / C , and Φ is a fixed user-defined angle, typically set to either 0 or π / 2 , which specifies the type of droop-like characteristic that the VO-controlled inverter exhibits in the steady state.
Further, assuming Φ = π / 2 , an averaged model is utilized by averaging the trajectory of the oscillator over one ac cycle, and then we can obtain the following amplitude and frequency dynamics:
d θ ¯ *   d t = κ v κ i 6 C V ¯ 2 P
where P is the average three-phase real power over one cycle ω * / 2 π .
Given (14) and the fact that θ t + ω t = θ * t + ω * t and d θ / d t = 0 in the steady state, we have
ω eq = ω * κ v κ i 2 C V ¯ oc 2 P eq
This VOC inverter uses symmetrical pulse-width modulation (SPWM) in the control. A reference signal (usually a sine-wave signal) is compared with a high-frequency triangular-wave signal, and the on-time of the switching device is controlled according to the comparison results.
Assuming that the reference signal is V ref and the high-frequency triangular-wave signal is V tri , their comparison yields a result. If V ref > V tri , the switching device conducts; if V ref < V tri , the switching device turns off.
By controlling the conduction time of the switching devices, the output voltage can be regulated. The longer the conduction time, the higher the output voltage; the shorter the conduction time, the lower the output voltage. Specifically, an SPWM signal can be derived from the comparison of V ref and V tri to control the conduction time of the switching devices. This modulation signal can be described by the following formula:
D = V ref V tri × 100 %
where D is the duty cycle representing the percentage of conduction time in the total period, V ref represents the reference signal, and V tri represents the triangular-wave signal. By controlling this duty cycle, the precise adjustment of the output voltage can be achieved, thus realizing SPWM control.
The converters based on VOC exhibit deviations from their nominal responses, which vary depending on the active power. These deviations are characterized by the slope of the droop curves, indicating that VOC droop characteristics are embedded in its nonlinear dynamics model. As a result, VOC has the primary frequency regulation characteristics of droop control. The VOC time-domain model is rapid, and when combined with its coupling droop characteristics, it enables the measurement of frequency without requiring the use of a Phase-Locked Loop (PLL). This feature significantly reduces the calculation burden of the controller. Additionally, the primary frequency regulation scale of VOC is shorter than that of droop control, allowing for faster attainment of the steady state.

3.2. Control Method

3.2.1. Overall Framework

This paper proposes a multi-level frequency control architecture suitable for a VOC strategy in power systems, as depicted in Figure 3. The proposed architecture consists of two levels: lower-level primary frequency control and upper-level secondary frequency control. At the lower level, a multistage frequency regulation architecture is used with VOC to achieve frequency regulation. Additionally, a virtual link is added to compensate for lower inertia power, and the dynamic adjustment of the virtual capacitance is performed to optimize the inertia of the frequency regulation dimension. To facilitate information interaction between distributed generation units, network layer communication is utilized.

3.2.2. Determination of Virtual Capacitance

Figure 4 shows the transient frequency characteristics of the power system before and after optimal inertia frequency regulation optimization. The red curve is the frequency characteristics before optimization, and the blue curve is the frequency characteristics after optimization. Time scales t 1 and t 3 are the frequency dip and frequency recovery scales before optimization, while time scales t 2 and t 4 are the frequency dip and frequency recovery scales after optimization, respectively. At the scale of t 1 or t 2 , the frequency change rate rises sharply after the disturbance appears and gradually decreases within this time scale. However, at the t 3 or t 4 scale, the frequency change rate gradually increases. This paper expects a longer frequency drop process with a smaller frequency deviation, that is, t 2 > t 1 and Δ f 2 < Δ f 1 . Also, the frequency recovery process is expected to be faster, that is, t 4 < t 3 .
Based on the analysis of the influence of VOC parameters on the control transient characteristics in Appendix B, it was found that virtual capacitance plays a crucial role and has a proportional relationship with the transient response time. In order to achieve the ideal optimization effect, the value of virtual capacitance should be larger during the frequency drop stage but smaller during the frequency recovery stage. To meet the requirements for inertia frequency regulation optimization, an adaptive virtual inertia control strategy based on the dynamic change in VOC parameters is proposed in this article, which is
C i n e r t i a = k J ω ω 0 ω ω 0 d Δ ω d t C c o n = C 0 + C i n e r t i a L = 1 C ω 0 2
where C i n e r t i a is the value of the virtual compensation capacitance. C 0 is the initial value of the virtual capacitor. C c o n is the final virtual capacitance value. Δ ω is the frequency variation. ω 0 is the rated angular frequency. k J is the inertia compensation coefficient.
By substituting (16) into Equation (15), the following equations can be obtained:
2 C 0 + C i n e r t i a V ¯ oc 2 k v k i ω ω * = P out
2 C i n e r t i a V ¯ oc 2 k v k i ω ω * = P out 2 C 0 V ¯ oc 2 k v k i ω ω *
Equation (17) represents the active power frequency sag equation of VOC, with the right side showing the coupling effect between active power and frequency, while the left side represents the superimposed virtual inertia power. The adaptive control method of virtual inertia allows for the automatic adjustment of the virtual capacitance based on changes in frequency variation. When the frequency change rate is high, a higher value of virtual capacitance can be set to increase the transient response time and reduce the frequency change rate simultaneously. As the frequency approaches stability, the frequency change rate gradually decreases, and the value of the virtual capacitance decreases accordingly, thus accelerating the transient response speed of the system.

3.2.3. Determination of Virtual Inductor

In Appendix B, we analyze the coupling droop characteristics of the VOC inverter, and here, we assume that the coupling droop relation of active power and frequency is
ω i = ω * m P ω P ¯ i m P ω = k v k i 2 C c o n V ¯ oc 2
where m P ω is the droop coefficient. ω i is the angular frequency for node i. P ¯ i is the active power of node i. ω * is the rated angular frequency.
As shown in (19), the droop coefficient of VOC is determined by the VOC parameters. Therefore, dynamically adjusting the VOC parameters can be used to adjust the frequency characteristics of a power system. Based on sparse communication networks, the upper-level control should be incorporated into dynamically adjustable VOC parameters. Assuming that the power-generating unit acts as a node, the adjacency matrix A forms a communication network connection between nodes i and j, with its elements a i j = a j i > 0 describing the network’s communication status. The upper-level frequency control based on VOC and consensus theory is as follows:
ω i = ω * + k v k i 2 C c o n V ¯ oc 2 P ¯ i Ω i
d Ω i d t = g f ω i ω * j = 1 n a i j Ω i Ω j
where Ω i is the control variable of the frequency controller. g f is the positive gain that affects the frequency adjustment speed of node i. a i j is the positive gain that affects the power-sharing rate between nodes i and j.
According to consensus theory, the state variables of both Ω i and Ω j will eventually converge to the same value. Specifically, each node of the droop characteristic curve will move by the same amount, ω i ω 0 , as depicted in Figure 5. Therefore, there is an agreement in the change in droop characteristics across all nodes.
When the system reaches a steady state, the derivative on the left-hand side of (21) becomes zero. This is also the point where the frequency deviation in (21) equals zero, which means that the system frequency and the rated angular frequency are equal to each other. Therefore, the system regulation achieves no difference, i.e., ω i = ω * . To control the upper-level controller, the VOC parameter is introduced, and the virtual capacitor C remains unchanged. On the other hand, the virtual inductor L can be adjusted, which results in the system frequency being determined by the adjustable virtual inductor L a d j . Based on the adjustable virtual inductor L a d j , the system frequency is
ω i = 1 L a d j C c o n
Assuming that the parameters of the LC resonance loop are L c o n and C c o n when the system frequency is the rated value, by substituting (22) into (20), one can obtain
L a d j = 1 / C c o n 1 L c o n C c o n + k v k i 2 C V ¯ i 2 P ¯ i Ω i 2
As can be seen from (23), the control variable Ω i of the controller is absorbed into the coupling droop characteristics of VOC through the adjustable inductance of the VO controller. When the system reaches a steady state, the adjustable inductance L a d j value will tend to be L c o n , and the frequency ω i will also tend to be the rated frequency ω * .

4. Case Study

4.1. System Description

This paper utilized the IEEE-39 power system in the simulation experiments. This system comprises 38 transmission lines, 10 generators, 3 wind farms, and 6 energy storage units. A detailed visualization of the power system is presented in Figure 6. The parameters of the generators are listed in Table 1. It is assumed that the wind speed in the wind farm follows a Weibull distribution, and the scale factor, shape factor, cut-in wind speed, rated wind speed, and cut-out wind speed are listed in Table 2. The parameters of the energy storage units are listed in Table 3. VOC is employed in the inverters of the energy storage units, and they are coordinated through the cyber network layer shown in Figure 7. The parameter settings of the VO controllers for the energy storage units are listed in Table 4.

4.2. Verification of Upper-Level Preventive Control

Suppose that high wind occurs at 9:00 throughout the power system, causing Line 14–15, Line 9–39, and Line 3–4 trips in succession. This results in the isolated operation of the local power grid, as shown in the following Figure 8. The occurrence process of successive failures that is measured in the actual power system is as follows:
(1)
At 9:28, Line 14–15 is broken and fails to reclose;
(2)
At 10:02, Line 9–39 is broken and fails to reclose;
(3)
At 10:42, Line 3–4 is broken and fails to reclose, and the whole power system is separated into two regions, as shown in Figure 8.
Figure 8. Separated power system after successive failures.
Figure 8. Separated power system after successive failures.
Sensors 24 01812 g008
To verify the effectiveness of the model of successive failures in the power system, this whole successive failure process is modeled by the proposed method according to Equations (A2) and (A3).
In Equations (A2) and (A3), the expectation value of the Poisson distribution λ is a key parameter. It represents the average rate of occurrences (such as line outages) per unit time. To determine λ precisely, historical data on the frequency of line outages, particularly under specific weather conditions, are required. These data form the basis for a reliable estimate of λ . The calculation of λ involves the following two steps:
(1)
Data Collection: Gather historical data on line outages. These data should include the number of outages that occurred over a significant period, along with the duration of each observation period. It is essential to consider the specific weather conditions during which these outages occurred, as λ may vary with different environmental factors.
(2)
Calculating the Average Rate: Divide the total number of observed outages by the total observation time to calculate the average rate of outages per unit time. The formula for this calculation is as follows:
λ = T o t a l   n u m b e r   o f   o u t a g e s T o t a l   o b s e r v a t i o n   t i m e
After the determination of λ , the probability distribution and probability density function of the lines can be calculated. As a result, Figure 9 and Figure 10 show the probability distribution and probability density function of Line 3–4, Line 9–39, and Line 14–15. Furthermore, the time interval between successive failures is obtained, as shown in Figure 11. Given a confidence level of 1 α = 0.95 , the failure time of Line 14–15 is calculated as 9:30, the time interval t L , 1 α between the successive failures of Line 14–15 and Line 9–39 is calculated as 30 min according to Equation (A7), and the time interval t L , 1 α between the successive failures of Line 9–39 and Line 3–4 is calculated as 45 min according to Equation (A7). After the successive failures of Line 14–15, Line 9–39, and Line 3–4, the whole power system is separated into two regions. It is evident that the proposed method accurately predicts the time of the line fault occurrence, with a maximum error not exceeding 3 min in comparison to the actual occurrence. Therefore, the proposed method is capable of providing a relatively precise estimation of the occurrence time of transmission line faults in the power system, thus furnishing essential prospective fault information for subsequent preventive control measures.
Furthermore, based on the predicted fault occurrence time of the line, the method described in Section 2 is applied for the preventive control of the power system, with a time interval of 15 min between adjacent schedules. Figure 12 depicts the output of generators and energy storage units from 9:00 to 10:45, while Figure 13 illustrates the active power flows in the transmission lines connecting Region A and Region B. It is observed that at 9:00, the flow through Line 14–15 remained relatively high, reaching 32.67 MW, whereas at the scheduling time of 9:15, prior to the fault occurrence on Line 14–15, the flow through Line 14–15 had already decreased to 13.46 MW. Consequently, when the fault occurred on Line 14–15 at 9:28, the power system was not significantly affected by the sudden outage of Line 14–15. Furthermore, at various times from 9:00 to 10:30 before the fault occurred on Line 3–4, the proposed method gradually reduced the flow through this line by appropriately scheduling the generators and energy storage, thereby incrementally increasing the output of internal generation resources in Region B. This ensured that, in the event of a sudden fault on Line 3–4, the internal load in Region B could still receive continuous power supply, thereby averting curtailment resulting from the insufficient ramping and adjustment capability of the internal generation resources, thereby avoiding the phenomenon of power limitation within the region. Consequently, the proposed method adequately considers the impact of line faults on the reduced connectivity of the entire power system, proactively adjusting the output of the generators and energy storage units to reduce the power of the transmission lines connecting the two regions, thereby preventing severe oscillations and widespread transfer of power when a line fault occurs, thereby ensuring the continuity and sustainability of the power supply and enhancing the stability and security of power system operations.

4.3. Verification of Virtual-Oscillator-Based Lower-Level Control

In this section, we evaluate the efficacy of the subordinate control approach. In order to align with the aforementioned upper-level preventive control outcomes, we considered the dispatch results of the power system at 9:45 as the initial state for subordinate control and conducted the validation of the proposed voltage- and frequency-based lower-level control method using Region B as the test system. Initially, two frequency events are defined, as described below.
Scenario I: The initial active power of the three energy storage units at Nodes 4, 8, and 14 of Region B are 216.03 MW, 226.24 MW, and 310.31 MW, respectively, which are the dispatch results of the power system at 9:45. The output of the wind farm at Node 9 changes from 400 MW to 300 MW suddenly at t = 100 s to emulate the frequency event.
Scenario II: The initial active power of the three energy storage units at Nodes 4, 8, and 14 of Region B are 216.03 MW, 226.24 MW, and 310.31 MW, respectively. The output of the wind farm at Node 9 changes from 300 MW to 335 MW suddenly at t = 100 s to emulate the frequency event.
Figure 14 shows the frequency and active power curves of the energy storage units under Scenario I using the proposed lower-level control method. The initial frequency of the microgrid system is 50.00 Hz. Then, the frequency drops dramatically at t = 100 s, and the frequency nadirs for the three energy storage units are about 49.885 Hz, 49.865 Hz, and 49.88 Hz, respectively. The frequency starts to recover after the frequency nadir points, and the frequencies of the three energy storage units are basically the same after the transient dynamic processes. The system frequency recovers to 49.92 Hz, 49.96 Hz, and 50.00 Hz at about t = 101.8 s, t = 104.5 s, and t = 120 s, respectively. As for the active power of the energy storage units, the initial values are 216.03 MW, 226.24 MW, and 310.31 MW, respectively. At t = 100 s, the active power of the energy storage units increases dramatically, and the stable values of the active power for the three energy storage units are about 236 MW, 248 MW, and 340 MW, respectively. Similarly, Figure 15 shows the frequency and active power curves of the energy storage units under Scenario II using the proposed lower-level control method. In this scenario, the initial frequency of the power system is set at 50 Hz. At t = 100 s, the frequency of the three energy storage units suddenly rises to approximately 50.05 Hz. This frequency level, however, is not exceeded, and the frequency then slowly decreases until it returns to the original frequency at around t = 120 s. The active power values of the energy storage units decrease significantly at t = 100 s, and it takes approximately 1 s to stabilize.
Furthermore, this study examines the performance of energy storage in other initial states of lower-level control. Specifically, one frequency event is defined, which is described below.
Scenario III: The initial active power of the three energy storage units at Nodes 4, 8, and 14 of Region B is 246 MW, 200 MW, and 273 MW, respectively. The output of the wind farm at Node 9 changes from 400 MW to 300 MW suddenly at t = 100 s to emulate the frequency event.
Figure 16 shows the frequency and active power curves of the energy storage units under Scenario III using the proposed lower-level control method. At t = 100 s, the frequency drops sharply and reaches its lowest point. The frequency begins to recover after reaching the lowest frequency, and the three energy storage units basically return to the original frequency at 120 s. The active power of the energy storage units increases significantly and begins to stabilize after increasing to the maximum value. Scenario III exhibits dynamic processes that are strikingly similar to those in Scenario I and Scenario II, indicating that the proposed approach can effectively regulate energy storage in different initial states. Furthermore, the proposed control method yields satisfactory and stable results in variable control experiments.

4.4. Comparison with Droop Control in Lower-Level Control

The output characteristics of the traditional converters of the energy storage units generally simulate the characteristics of a synchronous generator. Depending on the degree of accuracy of the simulation of the synchronous generator, the most representative ones are droop control and VSG control. In order to analyze the difference between the VOC strategy and a traditional control strategy, the droop control strategy has been selected as the compared method in this paper. Droop control is designed to simulate the steady-state droop characteristics of a synchronizer and is a type of proportional control in the phasor domain. Droop control serves as a pertinent comparison to VOC due to its contrasting approach. Characterized by decentralized control, simplicity, and proven performance under varied conditions, droop control serves as a robust benchmark to evaluate the advanced capabilities of VOC. This comparison is instrumental in highlighting the advantages of VOC, such as faster response, enhanced stability, and adaptability to modern grid demands.
In the simulation example, the VOC strategy of the energy storage units in Figure 6 is replaced with the droop control strategy. Meanwhile, considering the differences between the VOC strategy and the droop control strategy, the role of the sparse communication network in the coordination between the VOC strategy and the droop control strategy is temporarily ignored in this paper. That is, when the communication network is removed, all of the energy storage units adopt local control. To compare the VOC strategy and droop control strategy, Scenario IV and Scenario V, described below, were selected to test the two controllers.
Figure 17 and Figure 18 show the active power curves of the energy storage units under Scenario I using both methods. Figure 19 and Figure 20 show the frequency and active power curves of the energy storage units under Scenario II using both methods. By comparing the simulation results, it can be seen that the change process of frequency and active power in the VOC strategy are much faster than those in the droop control strategy. Additionally, the droop control strategy takes a longer time to return to the new stable value, and the maximum overshoot and adjustment time are significantly greater than those of the VOC strategy. Taking Figure 17 and Figure 18 as an example, the VOC strategy is able to stabilize the frequency within 2 s, in contrast to the droop control strategy, where the frequency remains fluctuating even after 20 s of adjustment. Moreover, the VOC strategy effectively avoids overshooting during the adjustment process. As illustrated in Figure 17, under the VOC strategy, the steady-state frequency is approximately 49.85 Hz, with the lowest frequency point being slightly below 49.85 Hz. In contrast, as shown in Figure 18, with the droop control strategy, the system’s frequency fluctuates within a range of about 49.5 Hz to 50.15 Hz, indicating a lower stability in frequency control. This comparison highlights the superior performance of the VOC strategy in maintaining frequency stability and preventing excessive fluctuations, making it a more efficient and reliable control method for power systems.

4.5. Influence of the Adaptive Control Parameters in Lower-Level Control

In this paper, an adaptive frequency regulation method with variable parameters is proposed. To clearly capture the indicators measuring the inertial frequency regulation response with different selections of parameters, the frequencies under Scenario IV and Scenario V, described below, using fixed control parameters and adaptive control parameters are compared.
Scenario IV: The initial active power of the three energy storage units at Nodes 4, 8, and 14 of Region B is 216.03 MW, 226.24 MW, and 310.31 MW, respectively. The output of the wind farm at Node 9 changes from 400 MW to 332 MW suddenly at t = 100 s to emulate the frequency event.
Scenario V: The initial active power of the three energy storage units at Nodes 4, 8, and 14 of Region B is 216.03 MW, 226.24 MW, and 310.31 MW, respectively. The output of the wind farm at Node 9 changes from 300 MW to 345 MW suddenly at t = 100 s to emulate the frequency event.
Figure 21 shows the frequency curves under Scenario IV using the fixed control parameters and the adaptive control parameters. Figure 22 shows the frequency curves under Scenario V using the fixed control parameters and the adaptive control parameters.
It can be observed that both the VOC controllers with fixed and adaptive control parameters enter the transient process at the same time. However, the performance of these two control strategies differs significantly. VOC with adaptive control parameters is able to restore the rated frequency in a shorter period compared to the VOC controller with fixed control parameters. Specifically, the VOC controller with adaptive control parameters takes about 15~20 s to restore the rated frequency, while the VOC controller with fixed control parameters takes much longer. This suggests that the adaptive virtual inertia control strategy is more effective than the traditional VOC strategy in regulating the frequency of the power system. Through the incorporation of virtual inertia elements into the control strategy, VOC controllers with adaptive control parameters are able to respond faster to changes in power fluctuations, resulting in a shorter adjustment time in the frequency response. This advantage is particularly pronounced under transient conditions when the system experiences sudden changes in power demand or supply.

5. Conclusions

This paper develops a multi-level operation method for improving the resilience of power systems under extreme weather conditions by incorporating preventive control and virtual oscillator (VO) technology. The method’s ability to precisely estimate the occurrence time of transmission line faults is a significant advancement, enabling the implementation of effective preventive controls. By considering the impact of line faults on system connectivity, the method proactively adjusts the output of generators and energy storage units, effectively reducing the power in transmission lines and preventing severe oscillations during faults. This ensures the continuity and sustainability of the power supply, thereby enhancing the overall stability and security of power system operations. In addition, the method demonstrates robust performance in regulating energy storage under various conditions and shows superior stability and speed in control experiments, particularly highlighting the effectiveness of VOC over traditional droop control. Specifically, compared with droop control, VOC demonstrates the ability to stabilize frequency more quickly, effectively avoids overshooting during the adjustment process, and achieves greater stability in frequency control. These results emphasize the method’s potential as a comprehensive solution for enhancing power system resilience in the face of increasing extreme weather events.
Although the proposed multi-level operation method enhances the resilience of power systems, we must admit that there are challenges in ensuring the stability of adaptive VOC under all possible circumstances. Specifically, our approach may struggle in situations like the tripping of a major generating station. To address this issue, future work could explore new control strategies to improve the system’s adaptability and resilience to sudden events. This could include the development of more advanced predictive models capable of the real-time monitoring and prediction of the behavior of power systems and the automatic adjustment of the VOC strategy upon detecting potential instability signals. For example, by integrating machine learning, it would be possible to design control systems that learn from past failure events and quickly adapt to new scenarios.

Author Contributions

Conceptualization, J.H.; data curation, L.X.; formal analysis, C.L., C.T. and G.Z.; investigation, Z.F.; methodology, C.L., D.Z., J.H., L.L. and G.Z.; software, J.H. and C.W.; validation, L.X.; writing—original draft, J.H. and L.X.; writing—review and editing, J.H. and C.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Modeling of Successive Failures in Power Systems

Under extreme weather conditions, the failure probability of a transmission line can be expressed using a Poisson model. This is because the Poisson model is a tool for characterizing stochastic and rare occurrences of events. In extreme weather conditions, the occurrence of transmission line faults may exhibit stochastic and rare characteristics. The Poisson model is well suited to capturing this stochastic and rare behavior with a certain degree of accuracy [30,31] and is therefore employed to describe the probability of transmission line faults under extreme weather conditions. Given a specific extreme weather event, the probability of a single-component failure, taking a line outage as an example, is described by a Poisson model as follows [32]:
P X = k = λ k k ! e λ E X = λ
where k represents the number of line failures; λ is the expectation value of the Poisson distribution.
Subsequently, we can formulate the Poisson process using the above as follows:
P X t = k = ( λ t ) k k ! e λ t
where P X t = k is the probability of k line outages in a time period 0 , A .
Let Y be the time interval between the k th and k + 1 th trips. The probability distribution and probability density function of Y are
P Y t = 1 P ( Y > t ) = 1 P X t = 0 = 1 e λ t
p t = λ e λ t t 0 0 t < 0
Let t i and t j be the tripping time of two lines; then, their probability density functions are
f 1 t = λ 1 e λ 1 t t 0 0 t < 0
f 2 t = λ 2 e λ 2 t t 0 0 t < 0
Then, the probability density function of the time interval between t i and t j , which is t = t i t j , can be written as
f t = + f x , x t d x = + f 1 x f 2 x t d x = λ 1 λ 2 λ 1 + λ 2 e λ 2 t t < 0 λ 1 λ 2 λ 1 + λ 2 e λ 1 t t 0
Further, the probability density function and probability distribution of t are
f t = λ 1 λ 2 λ 1 + λ 2 e λ 1 t + e λ 2 t
F t = 1 λ 2 λ 1 + λ 2 e λ 1 t λ 1 λ 1 + λ 2 e λ 2 t
Utilizing the interval estimation technique, one can determine the prediction interval for the duration between two line failures t L , 1 α , t U , 1 α with a confidence level set at 1 α . Herein, t L , 1 α serves as the gap between failures.

Appendix B. The Influence of VOC Parameters on Transient Characteristics

This part focuses on how the VOC parameters affect the inverter output transient process. The transient process refers to the flow of energy between the energy storage element and resistance. As such, it is predominantly influenced by the circuit’s initial energy, structure, and the parameters of the circuit elements, rather than external signals like the power supply. This means that the circuit structure can be simplified and analyzed with greater accuracy. To achieve this simplification, the circuit model of VOC is represented as a second-order circuit in parallel with RLC without power excitation, as shown in Figure A1.
Figure A1. Simplified RLC parallel circuit.
Figure A1. Simplified RLC parallel circuit.
Sensors 24 01812 g0a1
When i L ( 0 ) = I L 0 > 0 , the switch is set to “on” at t = 0 s. According to circuit theory, one can obtain
d 2 i L d t 2 1 R C d i L d t + 1 L C i L = 0
with the initial conditions
i L ( 0 ) = i L ( 0 + ) = I L 0 d i L d t t = 0 + = 1 L u ( 0 + ) = 1 L u ( 0 ) = 0
Suppose that the circuit has a zero-input response; that is, the forced component of u c is zero, and the characteristic equation and characteristic root of the differential equation are
p 2 1 R C p + 1 L C = 0
p = 1 2 R C ± 1 2 R C 2 1 L C
Let γ = 1 / ( 2 R C ) and ω * = 1 / L C ; the two root characteristics can be expressed as
p 1 = γ + γ 2 ( ω * ) 2 p 2 = γ γ 2 ( ω * ) 2
The above formula shows that the characteristic roots p 1 and p 2 are determined by the component parameters in the circuit ( R < 0 , L > 0 , C > 0 ) . p 1 and p 2 can be double roots, distinct real roots, or conjugate complex roots. In order to explore the influence of circuit parameters on the transient process, it is necessary to derive the time constant of the transient component.
Condition 1: Let γ > ω * ; the circuit structure parameters would satisfy R > 0.5 λ . In this circumstance, p 1 and p 2 are two different negative real roots, and the general solution of the equation is i L = A 1 e p 1 t + A 2 e p 2 t . With the initial conditions, one can obtain
A 1 = p 2 p 2 p 1 I L 0 A 2 = p 1 p 2 p 1 I L 0
It is obvious that A 1 > 0 , A 2 < 0 , A 1 > A 2 , and we can know from (A13) that p 1 < 0 , p 2 < 0 , p 1 < p 2 .
Condition 2: Let γ = ω * ; the circuit structure parameters would satisfy R = 0.5 λ . In this circumstance, p 1 and p 2 are equal and are the negative real roots, and the circuit is in a state between that with oscillations and that without oscillations.
Condition 3: Let γ < ω * ; the circuit structure parameters would satisfy R < 0.5 λ . In this circumstance, p 1 and p 2 are a pair of conjugate complex roots, and the general solution of the equation is i L = A e γ t sin ω d t + θ . This circuit is a process of vibration attenuation, and the attenuation rate is determined by the attenuation coefficient γ ; also, the time constant is τ = 1 / γ .
In view of the above three conditions, the transient process of Condition 3 is the wanted oscillating structure, and the time constant is τ = 1 / γ . Because the rated angular frequency ω * is constant, the time constant τ is inversely proportional to γ . In other words, τ is proportional to C.
Through considering Condition 3, we know that a small λ 0 results in a nearly sinusoidal output voltage in the steady state with a lower frequency deviation from the nominal value ω * and reduced harmonic content in the output terminal voltage. However, a small λ 0 slows down the dynamic response of the VOC inverter. Therefore, in the design of VOC parameters, the control parameters L and C should be reasonably designed according to the control requirements. The design principles of L and C are
max C | Δ ω | max min , C δ 3 : 1 min C C t rise max L = 1 C ω * 2
C | Δ ω | max min = 1 2 | Δ ω | max V ¯ oc   V ¯ min Q ¯ rated   P ¯ rated   C t rise max = t rise max 6 V ¯ oc V ¯ min   V ¯ oc 2 V ¯ oc 2 V ¯ min 2 C δ 3 : 1 min = 1 8 ω * δ 3 : 1 max V ¯ oc V ¯ min V ¯ oc 2 V ¯ oc 2 V ¯ min 2
where t rise   max is the maximum permissible rise time. δ 3 : 1 max is the maximum permissible ratio. V ¯ oc is the open-circuit voltage. P ¯ rated   is the rated real power. V ¯ min   is the voltage at the rated power. Q ¯ rated   is the rated reactive power. | Δ ω | max is the maximum frequency offset. t rise max is the rise time. δ 3 : 1 max is the ratio of the third-to-first harmonic.
It can be seen that the optimal parameters of VOC are affected by the operating state of the system, and it is necessary to adjust the VOC parameters in a timely manner according to the changes in the operating state of the system to achieve the optimal frequency regulation performance. It should be noted that Formula (A14) only provides the adjustment ranges of L and C, without determining the optimal control parameters. Therefore, an adaptive multi-level frequency regulation method and VOC parameter design method are proposed.

References

  1. Han, J.; Lyu, W.; Song, H.; Qu, Y.; Wu, Z.; Zhang, X.; Li, Y.; Chen, Z. Optimization of Communication Network for Distributed Control of Wind Farm Equipped with Energy Storages. IEEE Trans. Sustain. Energy 2023, 14, 1933–1949. [Google Scholar] [CrossRef]
  2. Pan, W.; Li, Y. Improving Power Grid Resilience Under Extreme Weather Conditions With Proper Regulation and Management of DERs—Experiences Learned From the 2021 Texas Power Crisis. Front. Energy Res. 2022, 10, 921335. [Google Scholar] [CrossRef]
  3. Bie, Z.; Lin, Y.; Qiu, A. Concept and Research Prospects of Power System Resilience. Autom. Electr. Power Syst. 2015, 39, 1–9. [Google Scholar]
  4. Roege, P.E.; Collier, Z.A.; Mancillas, J.; McDonagh, J.A.; Linkov, I. Metrics for energy resilience. Energy Policy 2014, 72, 249–256. [Google Scholar] [CrossRef]
  5. Panteli, M.; Mancarella, P.; Trakas, D.N.; Kyriakides, E.; Hatziargyriou, N.D. Metrics and Quantification of Operational and Infrastructure Resilience in Power Systems. IEEE Trans. Power Syst. 2017, 32, 4732–4742. [Google Scholar] [CrossRef]
  6. Reed, D.A.; Kapur, K.C.; Christie, R.D. Methodology for Assessing the Resilience of Networked Infrastructure. IEEE Syst. J. 2009, 3, 174–180. [Google Scholar] [CrossRef]
  7. Mahzarnia, M.; Moghaddam, M.P.; Baboli, P.T.; Siano, P. A Review of the Measures to Enhance Power Systems Resilience. IEEE Syst. J. 2020, 14, 4059–4070. [Google Scholar] [CrossRef]
  8. Arjomandi-Nezhad, A.; Fotuhi-Firuzabad, M.; Moeini-Aghtaie, M.; Safdarian, A.; Dehghanian, P.; Wang, F. Modeling and Optimizing Recovery Strategies for Power Distribution System Resilience. IEEE Syst. J. 2021, 15, 4725–4734. [Google Scholar] [CrossRef]
  9. Yan, M.; Ai, X.; Shahidehpour, M.; Li, Z.; Wen, J.; Bahramira, S.; Paaso, A. Enhancing the Transmission Grid Resilience in Ice Storms by Optimal Coordination of Power System Schedule with Pre-Positioning and Routing of Mobile DC DE-Icing Devices. IEEE Trans. Power Syst. 2019, 34, 2663–2674. [Google Scholar] [CrossRef]
  10. Noebels, M.; Quirós-Tortós, J.; Panteli, M. Decision-Making Under Uncertainty on Preventive Actions Boosting Power Grid Resilience. IEEE Syst. J. 2022, 16, 2614–2625. [Google Scholar] [CrossRef]
  11. Nazemi, M.; Moeini-Aghtaie, M.; Fotuhi-Firuzabad, M.; Dehghanian, P. Energy Storage Planning for Enhanced Resilience of Power Distribution Networks Against Earthquakes. IEEE Trans. Sustain. Energy 2020, 11, 795–806. [Google Scholar] [CrossRef]
  12. Karim, M.A.; Currie, J.; Lie, T.T. Distributed Machine Learning on Dynamic Power System Data Features to Improve Resiliency for the Purpose of Self-Healing. Energies 2020, 13, 3494. [Google Scholar] [CrossRef]
  13. Kamruzzaman, M.; Duan, J.; Shi, D.; Benidris, M. A Deep Reinforcement Learning-based Multi-Agent Framework to Enhance Power System Resilience using Shunt Resources. IEEE Trans. Power Syst. 2021, 36, 5525–5536. [Google Scholar] [CrossRef]
  14. Musleh, A.S.; Khalid, H.M.; Muyeen, S.M.; Al-Durra, A. A Prediction Algorithm to Enhance Grid Resilience Toward Cyber Attacks in WAMCS Applications. IEEE Syst. J. 2019, 13, 710–719. [Google Scholar] [CrossRef]
  15. Wang, Y.; Huang, L.; Shahidehpour, M.; Lai, L.L.; Yuan, H.; Xu, F.Y. Resilience-Constrained Hourly Unit Commitment in Electricity Grids. IEEE Trans. Power Syst. 2018, 33, 5604–5614. [Google Scholar] [CrossRef]
  16. Liu, X.; Shahidehpour, M.; Li, Z.; Liu, X.; Cao, Y.; Bie, Z. Microgrids for Enhancing the Power Grid Resilience in Extreme Conditions. IEEE Trans. Smart Grid 2017, 8, 589–597. [Google Scholar]
  17. Wang, F.; Chen, C.; Li, C.; Cao, Y.; Li, Y.; Zhou, B.; Dong, X. A Multi-Stage Restoration Method for Medium-Voltage Distribution System With DGs. IEEE Trans. Smart Grid 2017, 8, 2627–2636. [Google Scholar] [CrossRef]
  18. Dehghanian, P.; Aslan, S.; Dehghanian, P. Quantifying Power System Resilience Improvement Using Network Reconfiguration. In Proceedings of the IEEE 60th International Midwest Symposium on Circuits and Systems (MWSCAS), Boston, MA, USA, 6–9 August 2017; pp. 1558–3899. [Google Scholar]
  19. Huang, G.; Wang, J.; Chen, C.; Qi, J.; Guo, C. Integration of Preventive and Emergency Responses for Power Grid Resilience Enhancement. IEEE Trans. Power Syst. 2017, 32, 4451–4463. [Google Scholar] [CrossRef]
  20. Li, Q.; Zhang, X.; Guo, J.; Shan, X.; Wang, Z.; Li, Z.; Tse, C.K. Integrating Reinforcement Learning and Optimal Power Dispatch to Enhance Power Grid Resilience. In IEEE Transactions on Circuits and Systems II: Express Briefs; IEEE: Piscataway, NJ, USA, 2021; Volume 69, pp. 1402–1406. [Google Scholar]
  21. Wei, X.; Liang, Y.; Guo, L.; Zhang, D.; Jia, Y. Research on fault location method of distribution network based on Lambda algorithm. Electr. Meas. Instrum. 2022, 59, 71–76. [Google Scholar]
  22. Duanmu, J.C.; Yuan, Y. Calculation method of outage probability of distribution network based on real-time failure rate of equipment. Electr. Meas. Instrum. 2023, 60, 26–32. [Google Scholar]
  23. Zhang, L.; Dou, X.; Zhao, W.; Hu, M. Islanding detection for active distribution network considering demand response. Electr. Meas. Instrum. 2019, 56, 63–68; 75. [Google Scholar]
  24. Lu, M.; Dhople, S.; Johnson, B. Benchmarking Nonlinear Oscillators for Grid-Forming Inverter Control. IEEE Trans. Power Electron. 2022, 37, 10250–10266. [Google Scholar] [CrossRef]
  25. Lu, M. Virtual Oscillator Grid-Forming Inverters: State of the Art, Modeling, and Stability. IEEE Trans. Power Electron. 2022, 37, 11579–11591. [Google Scholar] [CrossRef]
  26. Nguyen, N.; Almasabi, S.; Bera, A.; Mitra, J. Optimal power flow incorporating frequency security constraint. IEEE Trans. Ind. Appl. 2019, 55, 6508–6516. [Google Scholar] [CrossRef]
  27. Hao, Z.; Guo, G.; Huang, H. A Particle Swarm Optimization Algorithm with Differential Evolution. In Proceedings of the 2007 International Conference on Machine Learning and Cybernetics, Hong Kong, China, 19–22 August 2007; pp. 1031–1035. [Google Scholar] [CrossRef]
  28. Mazumder, S.K.; Greidanus, M.D.R.; Liu, J.; Mantooth, H.A. Vulnerability of a VOC-Based Inverter Due to Noise Injection and Its Mitigation. IEEE Trans. Power Electron. 2022, 38, 1445–1450. [Google Scholar] [CrossRef]
  29. Johnson, B.B.; Sinha, M.; Ainsworth, N.G.; Dorfler, F.; Dhople, S.V. Synthesizing virtual oscillators to control islanded inverters. IEEE Trans. Power Electron. 2015, 31, 6002–6015. [Google Scholar] [CrossRef]
  30. Wang, J.; Heydt, G.T. Application of Poisson models to electric power system reliability evaluation. IEEE Trans. Power Syst. 2006, 21, 1493–1500. [Google Scholar]
  31. Hu, Y.; Wang, W.; Liu, Y. Short-term load forecasting using a hybrid model based on Poisson regression and artificial neural network. Appl. Energy 2017, 185, 1491–1500. [Google Scholar]
  32. Zhang, Z.; Cui, W.; Niu, S.; Wang, C.; Xue, C.; Ma, X.; Liu, S.; Bie, Z. Preventive control of successive failures in extreme weather for power system resilience enhancement. IET Gener. Transm. Distrib. 2022, 16, 3245–3255. [Google Scholar] [CrossRef]
Figure 1. The pre-control process of power flow in a transmission section.
Figure 1. The pre-control process of power flow in a transmission section.
Sensors 24 01812 g001
Figure 2. The VdP type of VO.
Figure 2. The VdP type of VO.
Sensors 24 01812 g002
Figure 3. Multi-level control architecture of power system.
Figure 3. Multi-level control architecture of power system.
Sensors 24 01812 g003
Figure 4. Comparison of transient frequency characteristics of power system before and after ideal inertia optimization.
Figure 4. Comparison of transient frequency characteristics of power system before and after ideal inertia optimization.
Sensors 24 01812 g004
Figure 5. The secondary frequency regulation process based on the VOC droop curve.
Figure 5. The secondary frequency regulation process based on the VOC droop curve.
Sensors 24 01812 g005
Figure 6. Topology diagram of power system.
Figure 6. Topology diagram of power system.
Sensors 24 01812 g006
Figure 7. Cyber connection of the energy storage units.
Figure 7. Cyber connection of the energy storage units.
Sensors 24 01812 g007
Figure 9. Probability distribution of Line 3–4, Line 9–39, and Line 14–15.
Figure 9. Probability distribution of Line 3–4, Line 9–39, and Line 14–15.
Sensors 24 01812 g009
Figure 10. Probability density function of Line 3–4, Line 9–39, and Line 14–15.
Figure 10. Probability density function of Line 3–4, Line 9–39, and Line 14–15.
Sensors 24 01812 g010
Figure 11. Time interval between successive failures of Lines 3–4 and 9–39 and Lines 9–39 and 14–15.
Figure 11. Time interval between successive failures of Lines 3–4 and 9–39 and Lines 9–39 and 14–15.
Sensors 24 01812 g011
Figure 12. Output of generators and energy storage units in different time periods.
Figure 12. Output of generators and energy storage units in different time periods.
Sensors 24 01812 g012
Figure 13. Active power flows in the transmission lines connecting Region A and Region B.
Figure 13. Active power flows in the transmission lines connecting Region A and Region B.
Sensors 24 01812 g013
Figure 14. Frequency and active power curves of the energy storage units under Scenario I using the proposed lower-level control method.
Figure 14. Frequency and active power curves of the energy storage units under Scenario I using the proposed lower-level control method.
Sensors 24 01812 g014
Figure 15. Frequency and active power curves of the energy storage units under Scenario II using the proposed lower-level control method.
Figure 15. Frequency and active power curves of the energy storage units under Scenario II using the proposed lower-level control method.
Sensors 24 01812 g015
Figure 16. Frequency and active power curves of the energy storage units under Scenario III using the proposed lower-level control method.
Figure 16. Frequency and active power curves of the energy storage units under Scenario III using the proposed lower-level control method.
Sensors 24 01812 g016
Figure 17. Frequency and active power curves of the energy storage units under Scenario I using the proposed lower-level control method without the communication network.
Figure 17. Frequency and active power curves of the energy storage units under Scenario I using the proposed lower-level control method without the communication network.
Sensors 24 01812 g017
Figure 18. Frequency and active power curves of the energy storage units under Scenario I using the droop-based primary frequency control method.
Figure 18. Frequency and active power curves of the energy storage units under Scenario I using the droop-based primary frequency control method.
Sensors 24 01812 g018
Figure 19. Frequency and active power curves of the energy storage units under Scenario II using the proposed lower-level control method without the communication network.
Figure 19. Frequency and active power curves of the energy storage units under Scenario II using the proposed lower-level control method without the communication network.
Sensors 24 01812 g019
Figure 20. Frequency and active power curves of the energy storage units under Scenario II using the droop-based primary frequency control method.
Figure 20. Frequency and active power curves of the energy storage units under Scenario II using the droop-based primary frequency control method.
Sensors 24 01812 g020
Figure 21. Frequency curves under Scenario IV using the fixed control parameters and the adaptive control parameters.
Figure 21. Frequency curves under Scenario IV using the fixed control parameters and the adaptive control parameters.
Sensors 24 01812 g021
Figure 22. Frequency curves under Scenario V using the fixed control parameters and the adaptive control parameters.
Figure 22. Frequency curves under Scenario V using the fixed control parameters and the adaptive control parameters.
Sensors 24 01812 g022
Table 1. Parameters of the generators.
Table 1. Parameters of the generators.
Parameter SymbolParameter Name and UnitValue
Node 30Node 31Node 32Node 33Node 34Node 35Node 36Node 37Node 38Node 39
P Upper limit of active power (MW)10406467256525086875805648651100
P Lower limit of active power (MW)0000000000
Q Upper limit of reactive power (Mvar)400300300250167300240250300300
Q Lower limit of reactive power (Mvar)140−10015000−10000−150−100
Table 2. Parameters of the wind farms.
Table 2. Parameters of the wind farms.
Parameter SymbolParameter Name and UnitValue
Node 9Node 18Node 23
P Capacity (MW)500400450
Q Scale factor10.78.710
LfShape factor3.974.74.2
LgCut-in wind speed (m/s)3.53.53.5
CfRated wind speed (m/s)151515
| Δ ω | max Cut-out wind speed (m/s)252525
Table 3. Parameters of the energy storage units.
Table 3. Parameters of the energy storage units.
Parameter SymbolParameter Name and UnitValue
Node 3Node 4Node 8Node 14Node 22Node 25
P Rated active power (MW)200300250350200100
Q Rated reactive power (Mvar)203025352010
Lf(μH)300240400300240400
Lg(μH)300240400300240400
Cf(μF)202020202020
Rd(Ω)1.51.51.51.51.51.5
ω 0 Rated angular velocity (rad/s) 100 π 100 π 100 π 100 π 100 π 100 π
| Δ ω | max Maximum angular velocity offset (rad/s) π π π π π π
Table 4. Parameter settings of the VO controllers for the energy storage units.
Table 4. Parameter settings of the VO controllers for the energy storage units.
Parameter SymbolParameter NameValueUnit
k v Voltage proportionality factor110-
σ Negative impedance5.0868Ω−1
α Cubic current source coefficient3.414A/V3
L0 58.9μH
C0 119.3mF
k i 1 The first current scaling factor0.152-
k i 2 The second current scaling factor0.304-
k i 3 The third current scaling coefficient0.456-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, C.; Zhang, D.; Han, J.; Tian, C.; Xie, L.; Wang, C.; Fang, Z.; Li, L.; Zhang, G. A Multi-Level Operation Method for Improving the Resilience of Power Systems under Extreme Weather through Preventive Control and a Virtual Oscillator. Sensors 2024, 24, 1812. https://doi.org/10.3390/s24061812

AMA Style

Li C, Zhang D, Han J, Tian C, Xie L, Wang C, Fang Z, Li L, Zhang G. A Multi-Level Operation Method for Improving the Resilience of Power Systems under Extreme Weather through Preventive Control and a Virtual Oscillator. Sensors. 2024; 24(6):1812. https://doi.org/10.3390/s24061812

Chicago/Turabian Style

Li, Chenghao, Di Zhang, Ji Han, Chunsun Tian, Longjie Xie, Chenxia Wang, Zhou Fang, Li Li, and Guanyu Zhang. 2024. "A Multi-Level Operation Method for Improving the Resilience of Power Systems under Extreme Weather through Preventive Control and a Virtual Oscillator" Sensors 24, no. 6: 1812. https://doi.org/10.3390/s24061812

APA Style

Li, C., Zhang, D., Han, J., Tian, C., Xie, L., Wang, C., Fang, Z., Li, L., & Zhang, G. (2024). A Multi-Level Operation Method for Improving the Resilience of Power Systems under Extreme Weather through Preventive Control and a Virtual Oscillator. Sensors, 24(6), 1812. https://doi.org/10.3390/s24061812

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