Next Article in Journal
Towards a Digital Product Passport Fit for Contributing to a Circular Economy
Next Article in Special Issue
Enhanced Stability Criteria of Network-Based Load Frequency Control of Power Systems with Time-Varying Delays
Previous Article in Journal
Research on Zero-Voltage Ride Through Control Strategy of Doubly Fed Wind Turbine
Previous Article in Special Issue
A Novel Structure of a Power System Stabilizer for Microgrids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Load Frequency Control of Multi-Region Interconnected Power Systems with Wind Power and Electric Vehicles Based on Sliding Mode Control

1
College of Electrical Engineering, Shanghai Dianji University, Shanghai 201306, China
2
College of Electrical and Electronic Engineering, North China Electric Power University, Beijing 102206, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(8), 2288; https://doi.org/10.3390/en14082288
Submission received: 6 March 2021 / Revised: 5 April 2021 / Accepted: 12 April 2021 / Published: 19 April 2021
(This article belongs to the Special Issue Advances in Power System Stability and Control)

Abstract

:
In recent years, wind power systems have been used extensively, which not only improve the efficiency of current conventional power generation systems, but also can save traditional fossil fuel resources. However, considering the instability of wind power, after being grid connected, it can easily cause an impact on the stability of the grid operation. Considering the above problems, this paper considers to make full use of the energy storage part of electric vehicles (EVs) to increase the stability of grid operation. Based on the mathematical model, this paper studies the load frequency control (LFC) problem of a multi-region interconnected power system with wind power and EVs. First, since the system states are difficult to be monitored, a state observer is designed to estimate the state. Based on this, the integral sliding mode controller (SMC) is designed to realize the LFC of the interconnected power system. Meanwhile, to obtain better control performance, this paper further analyzes and optimizes the controller parameters based on Lyapunov stability theory. At last, simulations are carried out for the power systems with two regions in Simulink. The results show that the designed controllers are effective to compensate the load demand disturbances. In addition, it is demonstrated that the battery storage of EVs can play the role of peak-shaving and valley-filling in LFC.

1. Introduction

With the development of grid technology, communication technology is gradually integrated into the grid. As the connection between different regions is strengthened, a multi-regional interconnected power system is gradually formed in [1,2]. Multi-region interconnection can realize the economic operation of power system by coordinating different generation costs among different regions. Meanwhile, the interconnected areas can support each other with the help of tie-line, which improves the reliability of power supply. However, the multi-regional interconnection of power systems also has some unavoidable problems. Failures in the operation of power generation equipment and instability of renewable energy sources tend to cause fluctuations in the generation process. The fluctuations may affect the stability of the entire power system frequency through the tie-line [3,4]. Therefore, under the complex environment, more stringent standards are proposed for the stable operation of modern power systems. For power generation, integration of various energy sources is to be accomplished. Meanwhile, it is essential to ensure the safety and the efficiency of the transmission in electricity. More importantly, to ensure the control performance, it demands higher requirements for LFC in interconnected power systems.
With the increasing maturity of the wind power technology, wind power has been widely used in many forms of renewable energy power generation. However, the randomness and uncertainty of wind energy tends to cause wind abandonment and load shedding. This phenomenon causes the imbalance of supply and demand of the system power, which may bring in the deviation of the power system frequency. With the large-scale wind power integration into the grid, it will affect the stability of the grid frequency from a certain extent [5,6,7]. Therefore, some scholars have researched the LFC of the power system after the grid connection of wind power in recent years. In order to track the stochastic fluctuations of wind power output, a model for decentralized load frequency prediction control of interconnected grids containing the wind farms was proposed in [8]. In [9], a decentralized SMC is designed to suppress the fluctuation of load frequency based on the construction of a power system model containing wind power. A robust LFC strategy based on event-driven communication is proposed in [10]. While reducing the transmission volume of network communication, the frequency stability of the power system is ensured. In [11], an online reinforcement learning method based on an adaptive dynamic model is proposed. It is shown that the method can suppresses the uncertainty caused by the large-scale access of wind power effectively.
Considering the uncertainty of wind power forecasting, a distributed demand-side management approach for smart grids is proposed in [12]. Through game-theoretic approach, the cost of consumer participation in demand-side management programs is minimized. According to the research mentioned in the above literature, it is mainly to adopt appropriate control strategies to solve the effect of wind power uncertainty on the system load frequency. However, considering the limitations of only relying on the controller to ensure system stability, we further consider supplementing the energy storage component to improve the robustness of the control system.
In addition, with the popularization of the concept of environmental protection, the market share of the electric vehicles (EVs) is rising year by year. According to statistics, except for the demand of daily traffic, nearly 90% of the EVs are not used. Therefore, LFC of power system with EVs has become a key research topic for researchers in recent years. The scale effect of EVs concentration is analyzed in [13]. It is verified that EVs energy storage can provide backup energy storage for the systems. In [14], a LFC model for the implementation of the controlled energy dynamic changes in EVs is constructed. The simulation shows that the EVs are able to switch between power and controllable load. Although the literature [13,14] has investigated the mode of electric vehicles, the model of wind power and electric vehicles both connected to the grid is still open. To solve the problems of inaccurate mode of EVs entry and high anti-interference requirements, a linear self-interference control method is proposed in [15]. Considering the influence of the droop control characteristics of EVs, a combined optimization method of EVs-assisted frequency control and secondary frequency control of conventional units is proposed in [16]. The literatures [15,16] focus on control performance of electric vehicles after grid integration. However, in the system that contains wind power and EVs, more effective control strategies need to be considered to improve the system performance.
Through the above analysis, it is found that when considering the system containing both wind power and electric vehicle, the problems of system model building and control strategy selection still need to be paid more attention. Meanwhile, the point should also be noted that, to reduce fluctuations in the frequency of the load, the traditional energy storage methods such as batteries are used in general. It not only increases the economic cost of grid construction, but also fails to make full use of the energy storage capacity of idle EVs. In addition, when the system suffers from load demand disturbances, the load frequency is also affected and fluctuates. For complex systems with perturbations, sliding mode control can effectively overcome effects of perturbations.
Therefore, to reduce the deviation of load frequency after area interconnection, it is necessary to study the LFC of the complex system containing wind power and EVs using the sliding mode control.
Motivated by the above discussion, the contribution of this work is specified as follows:
i.
After the wind power is connected to the grid, to solve the impact of its uncertainty on the system frequency stability, the energy storage of idle EVs is considered as the buffer in this paper. To study the LFC of grid, when the system suffers from load demand perturbation, the model of the system containing wind power and EVs is developed.
ii.
In practice, the state is difficult to be grasped due to the limitation of equipment. Based on the model built in (I), the state observer is designed in this paper.
iii.
SMC is an effective control method to deal with system disturbance. Therefore, in this paper, to achieve the goal of load frequency control, the integral sliding mode controller is designed.
iv.
By combining the Lyapunov stability theory, the parameters of the controller are further optimized. In this way, better control performance of both the error system and the observer is ensured in the paper.
In this paper, the structure of the paper can be organized into five sections. First, we describe the background and motivation of this paper in Introduction of the Section 1. Then, the Section 2 presents description of the model building. In the Section 3, we focus on the theoretical analysis. The Section 4 includes the simulation results and the analysis. Finally, through simulation verification, we summarize the conclusions in Section 5.

2. System Model and Description

2.1. Wind Power and EVs Model

For wind power systems, the main issue is the active power output. Therefore, the wind power dynamics model is established as follows [17]:
Δ P ˙ WTG i = 1 T WTG i Δ P w i 1 T WTG i Δ P WTG i
where Δ P WTG i is the wind turbine output power deviation, T WTG i is the wind turbine time constant, and Δ P w - i is the power of wind.
Normally, EVs tend to play the role of getting electricity from the grid. Benefit from the development of V2G technology, EVs energy storage enables participation in LFC. The main idea is to make full use of the battery energy when the EVs are non-working; to use it as a buffer for the fluctuations caused when renewable energy is integrated into the grid. But EVs cannot be integrated into the grid randomly. It not only needs to be managed centrally for EVs energy storage, but also combined with intelligent control between the grid. The control input of the EV part of this paper is expressed as:
u e i ( t ) = α e i u ( t )
where α e i , u ( t ) denotes EVs scale factor, control input, respectively.
To model EVs, the literature [18] focuses on the analysis of when EVs are connected to the grid. Meanwhile, to prevent the EVs charging and discharging simultaneously, the paper supplements the constraint. Instead, by establishing mathematical models, ref. [19] adequately equates the output power of EVs as energy storage part. To realize the control objective, considering the partial output of EVs, the EVs dynamics model is written in the following form [19]:
Δ P ˙ e i = 1 T e i Δ P e i + α e i k e i T e i u ( t )
where Δ P e i is the deviation of EVs output power, k e i is EVs gain factor, and T e i is EVs time instants.

2.2. Load Frequency Control Model with Wind Power and EVs

Interconnected power systems usually consist of multiple control areas. Each area is connected to another area by tie lines to overcome the mismatch between generation and demand. In traditional LFC systems, automatic generation control acts as the secondary controller to control the grid frequency. To solve the instability caused by the wind power in the grid in recent years, some new topology structure of the interconnected power system is established after the introduction of electric vehicles, which is shown in Figure 1.
Under the stable operation of the power system, the system is linearized and simplified. In this model, a generator is used to represent the overall performance of the generator set. In addition, the time delay phenomenon of the control signal is ignored. Therefore, mathematical models containing the wind power and EVs power systems are shown in Figure 2.
According to the system model established in Figure 2, the mathematical model of the system can be expressed as:
Δ f ˙ i = D i M i Δ f i + 1 M i Δ P m i 1 M i Δ P tie i + 1 M i Δ P WTG i + 1 M i Δ P e i 1 M i Δ P d i Δ P ˙ m i = 1 T ch i Δ P m i + 1 T ch i Δ P v i Δ P ˙ v i = 1 R i T g i Δ f i 1 T g i Δ P v i + α g i T g i u ( t ) Δ P ˙ tie i = 2 π j = 1 , i 1 N T i j Δ f i δ i = β i Δ f i + Δ P tie i Δ P ˙ WTG i = 1 T WTG i Δ P w i 1 T WTG i Δ P WTG i Δ P ˙ e i = 1 T e i Δ P e i + α e i k e i T e i u ( t )
where Δ f i is the frequency deviation, Δ P m i is the generator power deviation, Δ P v i is the control valve position deviation, Δ P tie i is the tie-line power exchange, δ i is the area control deviation, M i is the rotational inertia of the generator set, D i is the load damping coefficient, T g i is the governor time constant, T ch i is the turbine time constant, R i is the governor droop characteristic, β i is the frequency deviation factor, Δ P d i is the load demand disturbance, and α g i is the turbine proportionality factor. The summary table of nomenclature is shown in Appendix A.
By simplifying the above equation of state, the system can be further expressed as:
x ˙ ( t ) = A x ( t ) + B u ( t ) + D f ( t ) , y ( t ) = C x ( t ) ,
with
x ( t ) = x 1 ( t )   x 2 ( t )   x 3 ( t )   x n ( t ) T , y ( t ) = y 1 ( t )   y 2 ( t )   y 3 ( t ) y n ( t ) T ,
x i ( t ) = Δ f i   Δ P m i   Δ P v i   Δ P tie i   δ i   Δ P WTG i   Δ P e i T , y i ( t ) = Δ f i   δ i T , f ( t ) = Δ P d i   Δ P w i T ,
A i i = D i M i 1 M i 0 1 M i 0 1 M i 1 M i 0 1 T ch i 1 T ch i 0 0 0 0 1 R i T g i 0 1 T g i 0 0 0 0 2 π j = 1 , i 1 N T i j 0 0 0 0 0 0 β i 0 0 0 0 0 0 0 0 0 0 0 1 T WTG i 0 0 0 0 0 0 0 1 T e i , A i j = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 π j = 1 , i 1 N T i j 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ,
B = 0 0 α g i T g i 0 0 0 α e i k e i T e i , C = 1 0 0 0 0 0 0 0 0 0 0 1 0 0 , D = 1 M i 0 0 0 0 0 0 0 0 0 0 0 1 T WTG i 0 ,
where x ( t ) is system state, A is state matrix, B is system input matrix, C is output matrix, y ( t ) is system output matrix, u ( t ) is control signal, and D is coefficient matrix.
In this paper, we consider the effect of the disturbance component on the system. The disturbance part mainly has two parts: load demand disturbance and wind power magnitude. To help the derivation of the subsequent theory, assume that the perturbation considered in this paper is bounded, and it can be written as:
f ( t ) L i ,
where L i is a known constant.

3. Model Characterization

In this section, based on the model, we focus on the theoretical analysis of the state observer, the SMC, and the stability of the system. The flow chart is shown in Figure 3.

3.1. Observer Design

Since the integration of wind power and EVs, the complexity of the system model becomes complex. Meanwhile, due to the limitation of monitoring equipment and cost in engineering, the system state variable may not be accurately measured. To solve this problem, the state observer is designed as follows:
x ^ ˙ ( t ) = A x ( t ) + B u ( t ) + O ( y ( t ) y ^ ( t ) ) , y ^ ( t ) = C x ^ ( t ) ,
where O is the observer gain matrix.
Considering state observer (6), the following observer-based integral sliding surface is designed:
s ( t ) = S x ^ ( t ) S 0 t ( A B K ) x ( τ ) d τ ,
where S and K are matrices of appropriate dimensions to be designed.
Define the error variable e ( t ) = x ( t ) x ^ ( t ) . When the system state reaches the sliding surface, it can be obtained that:
s ˙ ( t ) = 0 ,
s ˙ ( t ) = S B u ( t ) + S O C e ( t ) + S B K x ^ ( t ) ,
Thus, the equivalent control can be written as:
u eq ( t ) = K x ( t ) ( S B ) 1 S O C e ( t ) .
Substituting (10) into (6), the dynamic equation of the system is obtained as follows:
x ^ ˙ ( t ) = ( A B K ) x ^ ( t ) + S ^ O C e ( t ) ,
where S ^ = I B ( S B ) 1 S .
According to (5), (6), and (10), the error system is obtained as follows:
e ˙ ( t ) = ( Α O C ) e ( t ) + F f ( t ) + B ( u ( t ) u eq ( t ) ) .
Based on the above analysis, the controller u ( t ) is designed such that the LFC of the system (5) can be achieved.

3.2. Reachability Analysis

To achieve the reachability of the sliding mode surface, a SMC is designed as follows:
u ( t ) = K x ^ ( t ) ( S B ) 1 S O C e ( t ) μ ( S B ) 1 s ( t ) ν ( S B ) 1 sgn ( s ( t ) ) ,
where μ > 0 and ν > 0 are positive control parameters to be designed.
The Lyapunov function is chosen as:
V ( t ) = 1 2 s T ( t ) s ( t ) .
The derivative of (14) can be written as:
V ˙ ( t ) = s T ( t ) s ˙ ( t ) s ( t ) ,
Noting that s ( t ) sgn ( s ( t ) ) = s ( t ) 1 s ( t ) and s T ( t ) s ( t ) s ( t ) 2 . It can be obtained from (15) that
V ˙ ( t ) = s T ( t ) s ( t ) ( S x ^ ˙ ( t ) S ( A B K ) x ^ ( t ) ) μ V ( t ) ν .
Considering (16), one further has
V ( t ) ( V ( 0 ) + ν μ ) e μ t ν μ .
Therefore, under controller (13), when t t * = μ 1 ln μ s ( 0 ) + ν ν is satisfied, the observer system trajectory can reach the sliding surface in finite time and remain there.

3.3. Stability Analysis

The sliding mode dynamics of the systems (5) can be described as an augmented system as follows:
e ˙ ( t ) x ^ ˙ ( t ) = A O C 0 S ^ O C A B K e ( t ) x ^ ( t ) + D 0 f ( t ) + 0 μ ν ( S B ) 1 s ( t ) + 0 ν B ( S B ) 1 sgn ( s ( t ) ) .
Ignoring the variable structure part, system (18) can be expressed as:
e ˙ ( t ) x ^ ˙ ( t ) = A O C 0 S ^ O C A B K e ( t ) x ^ ( t ) + D 0 f ( t ) .
According to the Lyapunov theory, if there exist positive definite matrices Q 1 , Q 2 , P 1 , and P 2 that satisfy the following LMIs:
s y m ( P 1 ( A O C ) ) + Q 1 < 0 , s y m ( ( A B K ) P 2 ) + Q 2 < 0 ,
then the stability of the system (19) can be guaranteed.
Considering the system (18), one has
e ( t ) = e ( Α O C ) t e ( 0 ) + 0 t e ( A O C ) ( t τ ) D f ( τ ) d τ , x ^ ( t ) = e ( A B K ) t x ^ ( 0 ) + 0 t e ( A B K ) ( t τ ) ( - μ B ( S B ) 1 s ( t ) ν B ( S B ) 1 sgn ( s ( t ) ) + S ^ O C e ( τ ) ) d τ .
Since A O C and A B K are Hurwitz matrix, there exist positive values λ e > 0 , λ o > 0 , M 1 and M 2 such that (22) holds. Moreover, the equation (21) can be rewritten as:
e ( t ) M 1 e λ e t + M 1 0 t e λ e ( t τ ) D f ( τ ) d τ , x ^ ( t ) M 2 e λ o t + M 2 0 t e λ o ( t τ ) ( μ B ( S B ) 1 ( s ( 0 ) + ν μ e μ τ ν μ ) + ν B ( S B ) 1 + S ^ O C e ( τ ) ) d τ .
Considering f ( t ) L i . By combing (17) with (22), one further has
e ( t ) M 1 e λ e t + M 1 0 t e λ e ( t τ ) D f ( τ ) d τ , x ^ ( t ) M 2 e λ o t + M 2 0 t e λ o ( t τ ) ( μ B ( S B ) 1 ( s ( 0 ) + ν μ e μ τ ν μ ) + ν B ( S B ) 1 + S ^ O C e ( τ ) ) d τ .
Taking the limit of (23), one has
lim t e ( t ) M 1 D L λ e , lim t x ^ ( t ) = M 1 M 2 S ^ O C D L λ o λ e .
From the definition e ( t ) = x ( t ) x ^ ( t ) , (24) can be expressed as:
lim t x ( t ) = M 1 D L λ e ( 1 + M 2 S ^ O C λ o ) .
Therefore, when LMI (20) is satisfied, the SMC (13) can guarantee the asymptotic stability of the closed-loop system (18). Moreover, the system state is within the region of the equilibrium with the radius M 1 D L λ e ( 1 + M 2 S ^ O C λ o ) .

3.4. Optimize the Control Parameters

Considering the SMC (13), when selecting larger values of μ and β , the system state can reach the sliding mode surface faster. At the same time, it is able to reduce the time of load frequency oscillation. However, when the system state reaches the sliding surface with a high speed, the deviation of the load frequency will become larger simultaneously.
In addition, the Lyapunov function of the system (19) is designed as:
V ( t ) = e T ( t ) P 1 e ( t ) ,
Taking the derivative of (26), one has
V ˙ ( t ) = e T ( t ) ( s y m ( P 1 ( A O C ) ) + Q 1 ) e ( t ) .
According to the Lyapunov stability theory, when the eigenvalue of Q 1 take larger values, it will enable the error system to converge quickly and improve the dynamic performance of the error system. However, a larger value of the eigenvalue λ e of Q 1 will contribute to larger parameter for the observer gain O . It is seen from the observer
x ^ ˙ ( t ) = ( A B K ) x ^ ( t ) + S ^ O C e ( t )
that the dynamic performance of the observer will be damaged. To improve the dynamic performance of the closed-loop system, a compromise solution is proposed.
Combining (17) and (26), it can be obtained that
V ˙ ( t ) e T ( t ) Q 1 e ( t ) λ m i n ( Q 1 ) e ( t ) 2 .
Considering the LMIs (20), with the increasing eigenvalues of Q 1 , λ e also increases. The parameters Q 2 and λ o have the same property. Therefore, by optimizing the trace of Q 1 and Q 2 , the larger λ e and λ o can be obtained. At the same time, the smaller value of S ^ O C can be guaranteed by optimizing the trace of the matrix O C .
From the discussion above, the matrix O , K and P 1 , P 2 , Q 1 , Q 2 can be solved by the following optimization problem:
max Q 1 ,   Q 2 t r ( Q 1 ) , t r ( Q 2 ) ; min O   t r ( O C ) s . t . s y m ( P 1 ( A ( k ) O C ) ) + Q 1 < 0 , s y m ( ( A ( k ) B ( k ) K ) P 2 ) + Q 2 < 0 ,
where k = 1 ,   2   ,     , M .
Through the above analysis, in the process of optimization of the controller parameters, we have supplemented the maximum values of the trace of Q 1 Q 2 and the maximum value of the trace of O C , denoted by max t r ( Q 1 ) ,   t r ( Q 2 ) and min t r ( O C ) , respectively. In this case, we can use the minimum value of t r ( Q 1 ) ,   t r ( Q 2 ) equivalent to the maximum value of t r ( Q 1 ) ,   t r ( Q 2 ) . In this way, after combining min t r ( O C ) , we use the “mincx” function in the LMI toolbox to solve the controller parameters.

4. Simulation Analysis

In this section, we consider the interconnected power systems with two regions containing wind power and EVs. The designed SMC is applied to the system model for LFC. For the conventional PID (Proportion Integral Differential) controller, by comparing the systems contain EVs or not, the simulations are given.
Assume the system subjected to perturbations from changes in load demand at t 0 . The magnitudes are Δ P d 1 = 0.01 pu and Δ P d 2 = 0.02 pu , respectively. In the two regions, the turbine scale factor is α g 1 = 0.9 , α g 2 = 0.6 , EVs scale factor is α e 1 = 0.1 , α e 2 = 0.4 , and wind power is Δ P w 1 = 0.04 pu , Δ P w 2 = 0.05 pu . According to the topology and the analysis in the mathematical model, the generator units are equivalent in each region in the simulation process. After idealized values for the model parameters, the system parameters of the two areas are shown in Table 1.

4.1. PID Control

The PID controller parameters are shown in Table 2. The simulation results of the two regional power systems are shown in Figure 3.
As shown in Figure 4, the PID controller is able to stabilize the load frequency deviation at t = 50   s . It is shown from the simulation that the system containing EVs can take less time to achieve stability. In addition, the system containing EVs also has less fluctuations. The effectiveness of utilizing energy storage of EVs to participate in frequency control is demonstrated.

4.2. Integral Sliding Mode Control

Next, the SMC will be used for simulation. The matrices in the two area controllers are as follows:
S 1 = 0 0 1 0 0 0 17 S 2 = 0 0 1 0 0 0 18 , μ 1 = 5 , μ 2 = 9 , ν 1 = 0 . 1 , ν 2 = 0 . 1 .
By solving (20), the sliding surface for region 1 and region 2 are obtained as:
s 1 ( t ) = 53.6373 Δ f ^ 1 85.3485 Δ P ^ m 1 37.1822 Δ P ^ v 1 + Δ P ^ ˙ v 1 + 23.6768 Δ P ^ tie 1 40.4814 δ ^ 1 + 99.7098 Δ P ^ WTG 1 114.8640 Δ P ^ e 1 + 17 Δ P ^ ˙ e 1 ,
s 2 ( t ) = 56.7137 Δ f ^ 1 89.5246 Δ P ^ m 1 39.0279 Δ P ^ v 1 + Δ P ^ ˙ v 1 + 24.2960 Δ P ^ tie 1 42.6785 δ ^ 1 + 104.6242 Δ P ^ WTG 1 120.5132 Δ P ^ e 1 + 18 Δ P ^ ˙ e 1 ,
The sliding controller for region 1 and region 2 are obtained as:
u 1 ( t ) = 3.9488 Δ f ^ 1 6.5906 Δ P ^ m 1 1.9060 Δ P ^ v 1 + 1.8283 Δ P ^ tie 1 3.1260 δ ^ 1 7.6996 Δ P ^ WTG 1 7.5571 Δ P ^ e 1 1.2406 ( Δ f 1 Δ f ^ 1 ) + 0.0186 ( δ 1 δ ^ 1 ) 0.2066 s ( t ) 0.0041 sgn ( s ( t ) ) ,
u 2 ( t ) = 4.1543 Δ f ^ 2 6.8601 Δ P ^ m 2 2.0328 Δ P ^ v 2 + 1.8618 Δ P ^ tie 2 3.2704 δ ^ 2 8.0172 Δ P ^ WTG 2 7.8554 Δ P ^ e 2 1.3495 ( Δ f 2 Δ f ^ 2 ) + 0.0110 ( δ 2 δ ^ 2 ) 0.3704 s ( t ) 0.0041 sgn ( s ( t ) ) ,
where Δ f ^ i is the estimation of Δ f i , Δ P ^ mi is the estimation of Δ P mi , Δ P ^ vi is the estimation of Δ P vi , Δ P ^ tie i is the estimation of Δ P tie i , δ ^ i is the estimation of δ i , Δ P ^ e i is the estimation of Δ P e i , Δ P ^ WTG i is the estimation of Δ P WTG i .
The simulation results under the SMC are shown in Figure 5, Figure 6 and Figure 7.
As can be seen from Figure 5, the integrated SMC enables the stabilization of the load frequency deviation in both regions. Additionally, from Figure 5, it is shown that the observed value Δ f ^ i tracks the parameter Δ f i effectively. Further combined with Figure 6, the observer proposed can achieve the estimation of the unknown state system.
In SMC, we consider the LFC problem of system including EV and without EV. To compare the system performance between PID and SMC, we summarized the data from the simulation results, which is shown in Table 3, compared with the results of PID, SMC designed can reach stability in shorter time. Also, the fluctuation phenomenon is smoother. In addition, through Figure 7, it is seen that SMC reduce the stable time for system containing EVs. According to the above analysis, better control performance of the SMC and the effectiveness of EVs participation in frequency control is verified.

4.3. Optimize the SMC

In this part, the simulation of the optimized SMC will be carried out. The matrix in the controllers is selected as follows:
S 1 = 0 0 1 0 0 0 17 ,   S 2 = 0 0 1 0 0 0 18 , μ 1 = 5 , μ 2 = 9 , ν 1 = 0 . 1 , ν 2 = 0 . 1
By solving (29), the optimized sliding surfaces of region 1 and region 2 are obtained as follows:
s 1 ( t ) = 2.4957 Δ f ^ 1 0.0002 Δ P ^ m 1 12.5012 Δ P ^ v 1 + Δ P ^ ˙ v 1 0.0512 Δ P ^ tie 1 0.9015 δ ^ 1 + 0.0005 Δ P ^ WTG 1 + 17.0251 Δ P ^ e 1 + 17 Δ P ^ ˙ e 1 ,
s 2 ( t ) = 2.5102 Δ f ^ 1 0.0004 Δ P ^ m 1 + 12.4598 Δ P ^ v 1 Δ P ^ ˙ v 1 0.0450 Δ P ^ tie 1 0.8845 δ ^ 1 0.0004 Δ P ^ WTG 1 + 17.9895 Δ P ^ e 1 + 18 Δ P ^ ˙ e 1 .
The optimized sliding controllers of region 1 and region 2 are obtained as follows:
u 1 ( t ) = 0.2445 Δ f ^ 1 0.0152 Δ P ^ m 1 0.0814 Δ P ^ v 1 + 0.0054 Δ P ^ tie 1 0.0907 δ ^ 1 1.3501 Δ P ^ WTG 1 0.1491 Δ P ^ e 1 4.5103 ( Δ f 1 Δ f ^ 1 ) 0.0021 ( δ 1 δ ^ 1 ) 0.2066 s ( t ) 0.0041 sgn ( s ( t ) ) ,
u 2 ( t ) = 1.2105 Δ f ^ 2 0.0135 Δ P ^ m 2 0.1437 Δ P ^ v 2 + 0.0098 Δ P ^ tie 2 0.1592 δ ^ 2 0.7510 Δ P ^ WTG 2 0.2619 Δ P ^ e 2 2.5012 ( Δ f 2 Δ f ^ 2 ) 0.1511 ( δ 2 δ ^ 2 ) 0.3704 s ( t ) 0.0041 sgn ( s ( t ) ) ,
The simulation results under the optimized SMC are shown in Figure 8, Figure 9 and Figure 10.
In this section, we first verify that the optimized observer can observe the state of the system, which is shown in Figure 8. According to Figure 8, it is shown that the optimized SMC can enable the load frequency deviation to reach stable values. Similarly, after optimization, the control performance was verified in system with and without EVs. By comparing the simulation results in Figure 9, for system containing EVs, it can achieve stability with smaller fluctuations and shorter time.
To further verify the effectiveness of the optimized SMC, we compare the control performance of the controller before and after the optimization. Simulations are performed in system with and without EVs, the results are shown in Figure 10. To compare the results, we have summarized the data from the simulation results, as shown in Table 4.
Through comparing the simulation results of in Figure 10a,b, in the case of the system without EVs, the parameter optimized SMC can make Δ f i stable in a shorter time. Similarly, in the case of the system containing EVs, the same conclusion can be obtained from the simulation results of in Figure 10c,d. From the above analysis, the effectiveness and the superiority of the parameter optimized controller are also further verified.

5. Conclusions

This paper aims to reduce the frequency fluctuation of systems with wind power and EVs. Through the above simulation analysis, we have solved the impact of monitoring equipment limitations, disturbances, and wind power uncertainty on the system. Further, the controller parameters have been optimized by combining the Lyapunov stability principle. The subsequent conclusions are drawn from the presented work:
  • The observer designed in this paper can realize the observation of the system state. Moreover, the observer-based SMC is designed in this paper. It not only realizes the stability in a certain time, but also has better control performance.
  • Through the simulation results, the EVs has better suppression effect on the load frequency fluctuation of the power system. Thus, the effectiveness of the view that EVs as storage participate in frequency control is also verified.
  • By comparing the figures before and after optimization, it is shown that the optimized SMC can enable the system load frequency reach stable value in shorter time, which improves the robustness of the interconnected power system.
With the above analysis, this paper realizes the LFC for the system with wind power and EVs. Also, we found that the energy storage of EVs can effectively participate in the frequency control. It is noted that through the simulation results of the system with two regions interconnected, we have verified the research objectives of this paper. In the future research, simulations of multiple region interconnections may be further considered.

Author Contributions

All authors contributed to the research in this paper. Y.L. conceived the methodology and model of the research system of this article. Z.W. performed simulation verification and analyzed the data. Z.Y. finished the original draft of the article. W.Y. edited and reviewed the article. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (grant no. 61803253, 61973209).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Nomenclature of system model.
Table A1. Nomenclature of system model.
Δ f i frequency deviation D i load damping coefficient
Δ P m i generator power deviation M i inertia constant
Δ P v i control valve position deviation T g i governor time constant
Δ P tie i tie-line power exchange T ch i turbine time constant
δ i area control deviation R i governor droop characteristic
Δ P WTG i wind turbine power deviation β i frequency deviation factor
Δ P w - i power of wind Δ P d i load demand disturbance
Δ P e i EVs output power deviation α g i turbine proportionality factor
k e i EVs gain factor α e i EVs proportionality factor
T e i EVs time instants

References

  1. Ourahou, M.; Ayrir, W.; Hassouni, B.E. Review on smart grid control and reliability in presence of renewable energies: Challenges and prospects. Math. Comput. Simul. 2020, 167, 19–31. (In Chinese) [Google Scholar] [CrossRef]
  2. Zhou, X.X.; Lu, Z.X.; Liu, Y.M. Development models and key technologies of future grid in China. Proc. CSEE 2014, 23, 4999–5008. (In Chinese) [Google Scholar]
  3. Li, Y.P.; Yan, Z.; Jia, Y.B. Review of multi-area interconnected power system reliability evaluation. East China Electr. Power 2009, 37, 1848–1854. (In Chinese) [Google Scholar]
  4. Yao, J.G.; Yang, S.C.; Shan, M.H. Reflections on operation supporting system architecture for future interconnected power grid. Autom. Electr. Power Syst. 2012, 37, 52–59. (In Chinese) [Google Scholar]
  5. Bai, H.B.; Wang, R.H. Influence of the grid-connected wind farm on power quality. Proc. CSU-EPSA 2012, 24, 120–124. (In Chinese) [Google Scholar]
  6. Li, J.H.; Zuo, J.J.; Wang, S. Analysis and assessment of operation risk for power system with large-scale wind power integration. Power Syst. Technol. 2016, 42, 3503–3513. (In Chinese) [Google Scholar]
  7. Tang, X.J.; Cai, J.C.; Ma, S.Y. Effects of DFIG wind power grid-connection on frequency response of power system. Proc. CSU-EPSA 2020, 32, 37–43, 62. (In Chinese) [Google Scholar]
  8. Yang, D.Y.; Cai, G.W. Decentralized model predictive control-based load frequency control for high wind power penetrated power system. Proc. CSEE 2015, 35, 583–591. (In Chinese) [Google Scholar]
  9. Han, Y.H.; Ma, C.; Zhu, Y.Z. Load frequency control for power systems with wind energy. Control Eng. China 2018, 25, 2046–2051. (In Chinese) [Google Scholar]
  10. Yang, C.; Yao, W.; Wen, J.Y. Event triggering based robust control of load frequency for interconnected power grid with wind power. Autom. Electr. Power Syst. 2018, 42, 57–64, 207–211. (In Chinese) [Google Scholar]
  11. Yang, L.; Sun, Y.Z.; Xu, J. Adaptive load frequency control of wind power system based on online reinforcement learning. Autom. Electr. Power Syst. 2020, 44, 74–89. (In Chinese) [Google Scholar]
  12. Scarabaggio, P.; Grammatico, S.; Carli, R. Distributed demand side management with stochastic wind power forecasting. IEEE Trans. Control Syst. Technol. 2021, 1–16. [Google Scholar] [CrossRef]
  13. Huang, Y.; Liu, J.Y.; Chen, J.R. Load frequency control considering vehicle to grid. Autom. Electr. Power Syst. 2012, 36, 24–28. (In Chinese) [Google Scholar]
  14. Zhang, Q.; Zhou, L.; Zhou, L.W. Load frequency control considering charging and discharging static frequency characteristics of electric vehicles. Autom. Electr. Power Syst. 2014, 38, 74–80. (In Chinese) [Google Scholar]
  15. Pan, Y.T.; Tan, W. Load frequency control considering vehicle to grid based on linear active disturbance rejection control. Control Eng. China 2018, 25, 1226–1231. (In Chinese) [Google Scholar]
  16. Yang, J.W.; Dong, H.Z.; Liao, K. Joint optimization of load frequency control considering auxiliary frequency regulation of electric vehicles. Power Syst. Technol. 2019, 39, 200–206. (In Chinese) [Google Scholar]
  17. Shang-Guan, X.C.; He, Y.; Zhang, C. Sampled-data based discrete and fast load frequency control for power systems with wind power. Appl. Energy 2020, 259, 114202. [Google Scholar] [CrossRef]
  18. Sperstad, I.B.; Korpås, M. Energy storage scheduling in distribution systems considering wind and photovoltaic generation uncertainties. Energies 2019, 12, 1231. [Google Scholar] [CrossRef] [Green Version]
  19. Aravindh, D.; Sakthivel, R.; Kaviarasan, B. Design of observer-based non-fragile load frequency control for power systems with electric vehicles. ISA Trans. 2019, 91, 21–31. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Topology of two-region power system including wind power and electric vehicles (EVs).
Figure 1. Topology of two-region power system including wind power and electric vehicles (EVs).
Energies 14 02288 g001
Figure 2. System load frequency control (LFC) model including wind power and EVs.
Figure 2. System load frequency control (LFC) model including wind power and EVs.
Energies 14 02288 g002
Figure 3. Flow chart of Section 3.
Figure 3. Flow chart of Section 3.
Energies 14 02288 g003
Figure 4. (a) Frequency deviation Δ f i with and without EVs in region 1; (b) frequency deviation Δ f i with and without EVs in region 2.
Figure 4. (a) Frequency deviation Δ f i with and without EVs in region 1; (b) frequency deviation Δ f i with and without EVs in region 2.
Energies 14 02288 g004
Figure 5. (a) Frequency deviation Δ f i and observed frequency deviation Δ f ^ i in region 1; (b) frequency deviation Δ f i and observed frequency deviation Δ f ^ i in region 2.
Figure 5. (a) Frequency deviation Δ f i and observed frequency deviation Δ f ^ i in region 1; (b) frequency deviation Δ f i and observed frequency deviation Δ f ^ i in region 2.
Energies 14 02288 g005
Figure 6. (a) Estimated error of frequency deviation e 1 in region 1; (b) estimated error of frequency deviation e 1 in region 2.
Figure 6. (a) Estimated error of frequency deviation e 1 in region 1; (b) estimated error of frequency deviation e 1 in region 2.
Energies 14 02288 g006
Figure 7. (a) Frequency deviation Δ f i with and without EVs in region 1; (b) frequency deviation Δ f i with and without EVs in region 2.
Figure 7. (a) Frequency deviation Δ f i with and without EVs in region 1; (b) frequency deviation Δ f i with and without EVs in region 2.
Energies 14 02288 g007
Figure 8. (a) Frequency deviation Δ f i and observed frequency deviation Δ f ^ i in region 1; (b) frequency deviation Δ f i and observed frequency deviation Δ f ^ i in region 2.
Figure 8. (a) Frequency deviation Δ f i and observed frequency deviation Δ f ^ i in region 1; (b) frequency deviation Δ f i and observed frequency deviation Δ f ^ i in region 2.
Energies 14 02288 g008
Figure 9. (a) Frequency deviation Δ f i with and without EVs in region 1; (b) frequency deviation Δ f i with and without EVs in region 2.
Figure 9. (a) Frequency deviation Δ f i with and without EVs in region 1; (b) frequency deviation Δ f i with and without EVs in region 2.
Energies 14 02288 g009
Figure 10. (a) Comparison of frequency deviation Δ f i without EVs in region 1. (b) Comparison of frequency deviation Δ f i without EVs in region 2. (c) Comparison of frequency deviation Δ f i with EVs in region 1. (d) Comparison of frequency deviation Δ f i with EVs in region 2.
Figure 10. (a) Comparison of frequency deviation Δ f i without EVs in region 1. (b) Comparison of frequency deviation Δ f i without EVs in region 2. (c) Comparison of frequency deviation Δ f i with EVs in region 1. (d) Comparison of frequency deviation Δ f i with EVs in region 2.
Energies 14 02288 g010
Table 1. Regional system parameters.
Table 1. Regional system parameters.
ParametersRegion 1Region 2
Tchi0.35 s0.352 s
Tgi0.08 s0.08 s
Ri5 Hz/up5 Hz/put
βi0.4 Hz/pu0.4 Hz/pu
Di0.02 Hz/pu0.01Hz/pu
Mi0.3 pu·s0.3 pu·s
Tei1 s1 s
TWTGi1.5 s1.52 s
Tij0.2 pu/rad0.22 pu/rad
Table 2. PID controller parameters.
Table 2. PID controller parameters.
SymbolRegion 1Region 2
Kp−0.6−0.6
Ki−1−1.1
Kd−1.2−1.3
Table 3. Summary of simulation results data of sliding mode controller (SMC) and PID.
Table 3. Summary of simulation results data of sliding mode controller (SMC) and PID.
PIDSMC
Overshoot (Hz)Stable Time(s)Overshoot (Hz)Stable Time (s)
without EVsregion 10.021549.700.015245.20
region 20.021948.600.020243.50
with EVsregion 10.019047.600.013040.30
region 20.020545.300.021239.80
Table 4. Summary of simulation results data of SMC.
Table 4. Summary of simulation results data of SMC.
Before OptimizationAfter Optimization
Overshoot (Hz)Stable Time(s)Overshoot (Hz)Stable Time (s)
without EVsregion 10.015245.200.014035.40
region 20.020243.500.024237.30
with EVsregion 10.013040.300.008031.45
region 20.021239.800.023530.70
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.; Liu, Y.; Yang, Z.; Yang, W. Load Frequency Control of Multi-Region Interconnected Power Systems with Wind Power and Electric Vehicles Based on Sliding Mode Control. Energies 2021, 14, 2288. https://doi.org/10.3390/en14082288

AMA Style

Wang Z, Liu Y, Yang Z, Yang W. Load Frequency Control of Multi-Region Interconnected Power Systems with Wind Power and Electric Vehicles Based on Sliding Mode Control. Energies. 2021; 14(8):2288. https://doi.org/10.3390/en14082288

Chicago/Turabian Style

Wang, Zhenghao, Yonghui Liu, Zihao Yang, and Wanhao Yang. 2021. "Load Frequency Control of Multi-Region Interconnected Power Systems with Wind Power and Electric Vehicles Based on Sliding Mode Control" Energies 14, no. 8: 2288. https://doi.org/10.3390/en14082288

APA Style

Wang, Z., Liu, Y., Yang, Z., & Yang, W. (2021). Load Frequency Control of Multi-Region Interconnected Power Systems with Wind Power and Electric Vehicles Based on Sliding Mode Control. Energies, 14(8), 2288. https://doi.org/10.3390/en14082288

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