Next Article in Journal
Effect of Surface Thermal Resistance on the Simulation Accuracy of the Condensation Risk Assessment for a High-Performance Window
Next Article in Special Issue
A Network Flow Model for Price-Responsive Control of Deferrable Load Profiles
Previous Article in Journal
Bayesian Energy Measurement and Verification Analysis
Previous Article in Special Issue
On Conceptual Structuration and Coupling Methods of Co-Simulation Frameworks in Cyber-Physical Energy System Validation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Allocation of Reserve Policies for a Multiple-Cell Based Power System

1
State Key Laboratory of Alternate Electrical Power Systems with Renewable Energy Sources, North China Electric Power University, Beijing 102206, China
2
Global Energy Interconnection Research Institute Europe GmbH, 10117 Berlin, Germany
3
Center for Electrical Power and Energy, DK2800 Lyngby, Denmark
*
Author to whom correspondence should be addressed.
Energies 2018, 11(2), 381; https://doi.org/10.3390/en11020381
Submission received: 6 December 2017 / Revised: 25 January 2018 / Accepted: 30 January 2018 / Published: 7 February 2018
(This article belongs to the Special Issue Methods and Concepts for Designing and Validating Smart Grid Systems)

Abstract

:
This paper applies a robust optimization technique for coordinating reserve allocations in multiple-cell based power systems. The linear decision rules (LDR)-based policies were implemented to achieve the reserve robustness, and consist of a nominal power schedule with a series of linear modifications. The LDR method can effectively adapt the participation factors of reserve providers to respond to system imbalance signals. The policies considered the covariance of historic system imbalance signals to reduce the overall reserve cost. When applying this method to the cell-based power system for a certain horizon, the influence of different time resolutions on policy-making is also investigated, which presents guidance for its practical application. The main results illustrate that: (a) the LDR-based method shows better performance, by producing smaller reserve costs compared to the costs given by a reference method; and (b) the cost index decreases with increased time intervals, however, longer intervals might result in insufficient reserves, due to low time resolution. On the other hand, shorter time intervals require heavy computational time. Thus, it is important to choose a proper time interval in real time operation to make a trade off.

1. Introduction

With increasing concerns regarding global warming and environment pollution, there has been a worldwide movement in the promotion of renewable technologies for electricity generation and for reducing the greenhouse-gas emissions. Many distributed generation units, including wind turbines [1], photovoltaic generators [2,3], fuel cells and fuel cell/gas/steam powered combined heat and power systems [4], are being used and connected to the power systems. However, a large penetration of intermittent and variable generation introduces operational challenges to the power system. To accommodate the variability of the renewables, it is important to harness the flexibility of the newly introduced units, such as batteries and electric vehicles, especially in the distribution network level [5,6,7]. Enabling this flexibility comes at the cost of an increasingly complex control system, characterized by many state and decision variables [8].
The problems encountered in the power systems have received much attention, and various efforts have been made to address the problems. These range from developing proper control schemes for individual component operation such as hierarchical aggregation method [9,10] to radical rethinking of system operations [11,12,13,14]. Traditionally, the system is kept secure by distinguishing the role between transmission system operators (TSO) and distribution system operators (DSO). For example, the TSO centrally controls a few big power plants through a supervisory control and data acquisition system. The DSO centrally manages the status of key devices, such as breakers, reference setting points of on/off load tap changers, capacity banks, etc. However, it is impossible for TSO to control the large number of distributed energy resources found today, as the grid control systems are centralized by design, and do not yet actively integrate distributed energy resources into the operation on a meaningful scale.
To address the aforementioned problem, the European FP7 project ELECTRA IRP proposes and develops a Web-of-Cells (WoC) architecture for operating the future power system [11,12,13,14]. In this approach, the power systems operation is divided into connected cells, each responsible for their own balancing and voltage control, thus establishing a robust, decentralized horizontal decomposition as opposed to the conventional centralized and vertical system operation. The WoC concept reformulates the control architecture of electric power systems to accommodate the challenges of fully distributed generation, reduced inertia, storage integration and flexible demand. Cells, which are defined as non-overlapping topological subsets of a power system, are associated with a scale-independent operational responsibility to contribute to system operation and stability. The operating state, including power exchanges and reserve parameters, can then be continuously optimized by coordination across cells.
In this paper, we study how joint energy and reserve scheduling, which are provided by flexible load units such as storage units and electric vehicles (EVs) that have limited power, energy and specific use patterns, could be operated more efficiently in the WoC architecture based power system. To do so, a robust power system reserve allocation approach [15,16] combining a predictive dispatch with optimal control policies in the form of linear decision rules (LDR) [17] is proposed. LDR concept is used in operations research field, where current states, past data or future predictions are combined linearly to make an operational decision. In [15], LDR-based reserve policies consisting of a nominal power schedule with a series of planned linear modification is proposed to accommodate fluctuating renewable energy resources. These policies are time-coupled, which exploits the temporal correlation of these prediction errors. The study showed that LDR-based reserve policies can reduce reserve operation cost compared to existing standard reserve operation method. In [16], the authors proposed an adjustable robust optimization approach to account for the uncertainty of renewable energy sources in optimal power flow. The optimized solution has two part: (1) the base-point generation is calculated to serve the forecast load which is not balanced by RES; and (2) the generation control using participation factors ensures a feasible solution for all realizations of RES output within a prescribed uncertainty set. However, both papers only applied the LDR method for one power system. Compared to stochastic programming, robust optimization method only requires knowledge of the range of variation of the uncertain parameters as opposed to an accurate specification of the uncertain parameter in stochastic programming. Therefore, robust optimization has been gaining popularity for decision making under uncertainty. In [18], the authors applied LDR-based reserve policies to the Web-of-Cells [19]. Firstly, the study shows that the method works fine for a single cell operation, i.e., the power and energy curve of batteries are within the capacity for any realized RES output. Then, three ad hoc cooperation strategies of web-of-cells are studied and compared using the LDR method. The three cooperation strategies include: (a) no cooperation between cells; (b) full cooperation between all cells; and (c) in between these two extreme cases. The study shows that Strategy (a) has a clear disadvantage over other two cooperation strategies.
Building forth on the previously developed work [18], this paper has two advancements: (1) It develops a model that adapts the application of the LDR-based reserve policies to multiple-cell based power systems rather than single-cell based power systems. Based on the proposed model, cross-cell reserve allocation, indicating the cooperation scheme among cells, can be determined. The results show that the involvement of the cross-cell reserve depends on the availability of reserve resources in the local cell. (2) To facilitate the real time operation in the real system, the effects of different time intervals on the LDR control policies are investigated in this paper. The investigation of the effects is made considering energy curves, power curves, and cost index of each discussed case. The remainder of the paper is organized as follows. In Section 2, the optimization formulation for one and multiple cell based power system is proposed, given the basic power system model. Comprehensive case studies are performed in Section 3. Conclusions are drawn in Section 4.

2. Methodology

The methodology used for solving the robust optimal reserve allocation problem is introduced in this section. The methodology is based on LDR, and can be applied to the operation of power systems, which could be single-cell based or multiple-cell based. A single-cell based power system can be considered as an isolated system, which could be an autonomous microgrid or a regional network. A multiple-cell based power system consists of two or more cells linked via interconnectors between each other. In the following subsections, the LDR based methodology is firstly proposed. Then, the optimization problem regarding single-cell and multiple-cell based power systems are, respectively, formulated using proposed methodology.

2.1. Basic Power System Model

A power system with various participants connected to a power grid is considered. The participants of a power system could be production units, loads, or storage units, which either inject power into or extract power from a node in the network. They are categorized into two types in terms of power injection: participants with inelastic power injection and those with elastic power injection. The inelastic power flows indicate the power flows of the participants cannot be regulated by control signals. These participants could be certain loads and renewable generators, such as wind turbines and PV panels. Regarding a participant i of this kind, the power injection into or extraction from a network y i can be modeled as:
y i = r i + G i δ
where r i R T indicates a nominal prediction of the power injection or extraction; T is the divided discrete time steps of a planning time horizon, over which electricity can be traded on intra-day markets [20]; δ R N δ T is the random forecast error vector; G i is a linear function used for mapping the uncertainty δ to power flows; and N δ is the number of elements in the uncertainty vector at a given time. If the power flows of the participants can be perfectly predicted, the prediction error δ will be zero.
The elastic power injection indicates that the power flows of the participants can be influenced by control signals. In other words, the flexibility of these participants can be exploited and used to mitigate the disturbance in the network. This can be achieved using the control signals determined by the results of proper optimization. According to [15], the elastic power injections can be modeled as C j x j , where x j R n j T is a vector of future states of participant j, n j is the state dimension, and C j R T × n j T is the stacked output matrix used for selecting the needed element of the state vector x j . The future state vector x j is given as:
x j = A j x 0 j + B j u j
where x 0 j is the current state of the participant j, u j R T is the control input to the elastic power participant j for balancing the system, and A j and B j are the corresponding stacked state transition matrices.

2.1.1. Constraints for Production Units

Different participants need to be limited by corresponding constraints. Regarding a production unit at period t, which could be an inelastic or elastic power unit, the upper and lower bounds of the power are imposed by the following constraints:
P i min y i P i max
P j min C j x j P j max

2.1.2. Constraints for Storage Units

Regarding a storage unit j, which could be a battery, an electric vehicle, etc., whose power injection is elastic, the storage unit’s power and energy constraints need to be imposed:
P j min C j x j P j max
E j min E j , t 1 + C j x j Δ t E j max
As all participants are connected to a network, the sum of all the inelastic power injection y i and elastic power injection C j x j has to be zero at all times t = 1 , , T . This can be achieved using an equality constraint:
i = 1 N inelas y i + j = 1 N elas C j x j = 0
where N inelas is the number of the inelastic power participants in a network, and N elas is the number of the elastic power participants.

2.2. Linear Decision Rule Based Robust Optimization of Reserve Allocation

Inserting Equations (1) and (2) into Equation (7), the following equation can be obtained:
i = 1 N inelas ( r i + G i δ ) + j = 1 N elas C j ( A j x 0 j + B j u j ) = 0
As mentioned, to keep power balance in a network, the power injection or extraction of the inelastic power units, i = 1 N inelas ( r i + G i δ ) , must be balanced by the power contribution of the elastic power units, j = 1 N elas C j ( A j x 0 j + B j u j ) , at any point in time. Regarding elastic power units, u j is the control input that can regulate the power flow of the elastic power participant j. According to the LDR method, control input signal u j can be expressed to policies of the affine form:
u j = D j δ + e j
where u j is described by a nominal schedule e j R T plus a linear variation D j R T × T , the nominal schedule e is mainly used for balancing the nominal prediction of the inelastic power flows r, and D is the dynamic response to the prediction errors δ . The matrix D defines a map from the uncertainty into the realization of the power contribution of the elastic power units. In order for the use of future disturbances to be causal, D j takes the lower-triangular form.

2.3. One Cell-Based Reserve Allocation Model

In a one-cell based power system, the reserve allocation optimization problem is formulated as a cost function as follows:
min E { j ϕ ( α j P j ( δ ) T P j ( δ ) + β j P j ( δ ) + γ j I ¯ ) }
where j ϕ indicates the elastic power participant; P i is the power contribution of the participant; α j , β j , and γ j are stacked vectors of quadratic, linear, and constant coefficients of the cost function of the power contribution, respectively; and I ¯ is a vector of all of them. This optimization problem is subject to system constraints, power constraints, and energy constraints, as mentioned in the previous section. The power contribution of the elastic power participant can be given as:
P j ( δ ) = C j x j ( δ )
with the help of Equation (2), P j ( δ ) can be further written as:
P j ( δ ) = C j ( A j x 0 j + B j u j )
According to policies of the affine form, as given in Equation (9), the entire optimization problem in Equation (10) becomes tractable due to the restricted variety of candidate u j .
Substituting Equation (12) into the objective function, i.e. Equation (10), the following form can be obtained:
min E { j ϕ ( α j ( C j ( A j x 0 j + B j u j ) ) T ( C j ( A j x 0 j + B j u j ) ) + β j C j ( A j x 0 j + B j ( D j δ + e j ) ) + γ j I ¯ ) }
where u j can be further expressed as policies of the affine form, as shown in Equation (9); together with an assumption that E [ δ ] = 0 , the optimization problem can be written as:
min j ϕ α j ( a j T a j + D j T b j T b j D j , E [ δ δ T ] ) + β j a j + γ j I ¯
with
a j = C j A j x 0 j + C j B j e j
b j = C j B j
where X , Y is the trace of product X Y ; regarding the reformulation of the equality and inequality constraints, the approach is similar to the one presented in [15,18].

2.4. Multiple Cells-Based Reserve Allocation Model

Compared to the single cell-based power system, the robust optimization of reserve allocation using LDR in multiple cells-based is more complicated. An example of a three-cell power system is shown in Figure 1. Each cell has its respective generations, loads, and storage units. Three cell are connected with tie lines, as depicted in the Figure 1. Regarding the resources with flexibilities, such as electric vehicles and batteries, it is assumed that a portfolio of all flexible energy-constrained resources in one cell can be represented by an aggregator as a single unit.
In the WoC concept, local imbalances within a cell are supposed to be solved locally. Cells are not necessarily supposed to be self-sufficient, but are required to have the balancing capability provided by elastic power units for mitigating deviations from a given schedule. The basic concept of reserve allocation in multi-cell-based power systems is that the local reserve provided by local elastic power units is prioritized for handling of local imbalances, which is reflected in the cost function of each available power unit. As long as local reserves can handle the local imbalances, no reserve is needed from other cells. Cross-cell reserve allocation happens when the local resource cannot handle the local imbalance.
Cells are managed by so-called Cell System Operators (CSOs), whose roles incorporate the tasks of traditional Distribution and Transmission System Operators (DSOs and TSOs, respectively) [21]. A dedicated Cell Controller (CC) performs monitoring and automated balancing tasks. Each cell is assigned one CC, which is managed by one CSO using bi-directional communication, as shown in Figure 2. A CSO, on the other hand, may be responsible for more than one CC to allow for flexible adaption to present-day grid partitioning and management schemes.
Among other tasks, CSOs procure reserves and send the schedules to the CCs, which automatically carry out balancing tasks around the given setpoints and trajectories. As indicated in Figure 2, information and measurements from each cell controller, including the prediction of the generation of the renewables, availability of the elastic power units as well as their power and energy constraints, information of loads, etc., are transferred to the cell operator, based on which the cell operator can utilize the proposed optimization to distribute the control signals, i.e., the D and e presented in Equation (9), for each cell controller to allocate the reserve in each cell.
The reserve optimization problem for multiple-cell based power system can be generally formulated as:
min E { l ϕ cell j ϕ l m ϕ cell , l m ( α l , l j P l , l j ( δ ) T P l , l j ( δ ) + β l , l j P l , l j ( δ ) + γ l , l j I ¯ + v l , m j ( α l , m j P l , m j ( δ ) T P l , m j ( δ ) + β l , m j P l , m j ( δ ) + γ l , m j I ¯ ) ) }
where l and m indicate the cell indices; j is the participant involved in the reserve allocation; ϕ cell and ϕ l are the sets of cell and participants in cell l, respectively; P l , l is the allocated resource in cell l that is reserved for the imbalance in cell l; and P l , m is resource that is reserved from cell l to cell m, also called the cross-cell reserve service. Finally, v l , m j is a binary decision variable used for determining the involvement of the cross-cell reserve allocation:
v l , m j = 0 P l , l j ( δ max ) > = P max j 1 P l , l j ( δ max ) < P max j
Similar to Equation (12), P l , l and P l , m can be expressed as:
P l , l j ( δ ) = C l , l j ( A l , l j x 0 , l , l j + B l , l j u l , l j )
P l , m j ( δ ) = C l , m j ( A l , m j x 0 , l , m j + B l , m j u l , m j )
According to the LDR method, u l , l j and u l , m j can be expressed as:
u l , l j = D l , l j δ + e l , l j
u l , m j = D l , m j δ + e l , m j
where participant j’s power schedule in local cell u l , l j is determined by a nominal schedule e l , l j , a linear variation D l , l j and predication errors δ . Similarly, expression of u l , m j can be obtained. Together with the assumption that E [ δ ] = 0 , the optimization problem can be written as:
min l ϕ cell j ϕ l m ϕ cell , l m ( ( α l , l j ( ( a l , l j ) T a l , l j + ( D l , l j ) T ( b l , l j ) T b l , l j D l , l j , E [ δ δ T ] ) + β l , l j a l , l j + γ l , l j I ¯ ) + v l , m j ( α l , m j ( ( a l , m j ) T a l , m j + ( D l , m j ) T ( b l , m j ) T b l , m j D l , m j , E [ δ δ T ] ) + β l , m j a l , m j + γ l , m j I ¯ ) )
with
a l , l j = C l , l j A l , l j x 0 , l , l j + C l , l j B l , l j e l , l j
b l , l j = C l , l j B l , l j
a l , m j = C l , m j A l , m j x 0 , l , m j + C l , m j B l , m j e l , m j
b l , m j = C l , m j B l , m j
Equation (23) establishes an overall optimization for a multiple-cell based power system. To solve the one-cell and three-cell optimization problems, YALMIP, a toolbox in MATLAB, is used [22,23].

3. Case Studies

In this section, several case studies are carried out. The setup is based on a configuration of SYSLAB, which is a laboratory testing facility for Smart Grid concepts located at the Risø campus of the Technical University of Denmark.

3.1. SYSLAB System

Figure 3 shows the SYSLAB facility [24,25], which is a 400 V three-phase grid designed for studying advanced grid control and communication concepts. The facility has 16 busbars and 116 automated coupling points. A wide range of distributed energy resources, such as wind turbines, solar panels, electric vehicles, etc., can be remotely operated via a distributed monitoring and control platform. Because of the very flexible connection interface of the busbars, various topologies of the system can be configured and operated. Division of the cells in the system can also be flexible and different from the one shown in Figure 3.

3.2. One Cell-Based Simulation

Figure 4 shows one-cell based isolated system, which comprises a battery, a solar PV, an EV, and a mobile load. The mobile load is used for imitating a typical residential electric load profile. The imbalance in the system is mainly caused by the difference between the fluctuated PV generation and the electric load demand. In this system, the inelastic power unit is the PV panel, whose power generation cannot be regulated by control signals. The EV and the battery act as elastic power units, which can provide reserve service due to their controllability of the power flow. The simulation is carried out for 30 min from 11:30 to 12:00. The time interval is chosen to be 3 min, and the horizon length therefore is 10.
From the simulation, the policies for the reserve allocation, i.e., the matrix D introduced in Section 2, of the battery and the EV are presented in Figure 5. As described in Section 2.2, D is a lower-triangular matrix that defines a map from the uncertainty to the realization of power contribution of a reserve participant. The matrix is visualized in Figure 5 for both the battery and the EV. According to Equation (9), the dynamic response to the prediction errors at time instant t is determined by [ D j ] t , 0 , [ D j ] t , 1 , , [ D j ] t , t , which are the elements in row t in the figure. The battery is allocated with more reserve than the reserve allocation contributed by the EV. This is because the battery has a wider power range and larger capacity. Furthermore, the reserve cost of the policy based scheme over the simulation period is compared to that of a flexible-rate scheme. Flexible rate scheme is commonly used for reserve optimization [15,26]. In this paper, the results given by flexible rate scheme act as a reference case. The differences between the two schemes are summarized as follows:
  • Flexible-rate reserves [15]: D j is a diagonal matrix. This indicates the best possible response to uncertainty without time coupling. The previous uncertainty therefore has no impact on the present operation because the causality of the uncertainty is omitted. The optimization is over the elements of e j and the diagonal parts of D j .
  • Policy-based reserves: Compared to the above scheme, this scheme considers the time coupling by taking D j as the lower-triangular form. It allows full exploitation of the information that will be available at each time step when the reserve is deployed.
The comparison results regarding the two schemes are given in Table 1. It is shown that the cost index given by the policy-based reserve is smaller than that given by the flexible-rate reserve. This is due to the full exploitation of the covariance matrix using the lower-triangular form of D. From the viewpoint of optimization, the feasibility region of D is larger for policy-based reserves, which gives the optimization more opportunity to find better results.

3.3. Three Cells-Based Simulation

A three-cell based isolated system is presented in Figure 6. Cell 1 has a vanadium battery described in Table 2 and an Aircon wind turbine that has a maximum power generation of 9.8 kW. Cell 2 comprises an EV and a solar panel. Cell 3 only has a mobile load in it. The basic idea of the co-operation is that the local imbalance has priority to be handled by using local reserve, while the cross-cell reserve is allocated only when there is a need. To demonstrate the co-operation of multiple cells, two scenarios are developed by setting the maximum power availability of EV. A three-hour simulation is carried out. The horizon length is set to 12, and the time interval is 15 min.

3.3.1. Scenario 1

In this scenario, P min and P max of the EV are set to −15 kW and 15 kW, respectively. With this power availability, the imbalance in Cell 2 caused by the fluctuated power generation of the PV panel can be handled locally. Similarly, Cell 1 can also handle its local imbalance due to the large power capacity of the vanadium battery. Therefore, it is expected that there will be no cross-cell reserve allocation between Cell 1 and Cell 2. Cell 3, however, has only an inelastic power unit, which is the mobile load that represents the residential electric load demand. As it is assumed that the load profile can be precisely predicted, there is no reserve allocation between Cell 3 and other cells, but only nominal schedule e 13 and e 23 exist for supporting the residential load. Figure 7 shows the simulation results of Scenario 1. Because of the sufficient flexibility in each cell, the cross-cell reserve allocation does not exist, which leads to D 12 = 0 and D 21 = 0 , as depicted in the figure. The imbalances inside Cell 1 and Cell 2 are all handled in their own cell; it can be seen that D 11 = I and D 22 = I , indicating that the predication errors in the two cells are fully compensated using the vanadium battery and the EV in their respective cells.

3.3.2. Scenario 2

To create the need for cross-cell reserve allocation, in Scenario 2, the P min and P max of the EV in Cell 2 are set to −5 kW and 5 kW, respectively. In this case, the local reserve in Cell 2 cannot handle the fluctuated power generation of the PV. Therefore, reserve allocation from Cell 1 to Cell 2 is necessary. In other words, the elastic resource in Cell 1 needs to be utilized to support the compensation of the predication errors in Cell 2. For that purpose, the vanadium battery in Cell 1 is divided virtually into two parts. One part, corresponding to D 11 , is used for allocating the reserve for the imbalance in Cell 1. The other part, corresponding to D 12 , is to reserve the imbalance in Cell 2. This reserve might be the power flow across Cells 1 and 2 via the tie line 1–2. The simulation results are presented in Figure 8. In Cell 1, the local reserve is sufficient for handling the local imbalance, the D matrix for Cell 1 is D 11 = I . The imbalance in Cell 2, as expected, is covered by local reserve in Cell 2 and some reserve from Cell 1, as can be seen in Figure 8 that D 12 + D 22 = I .

3.4. Impact of Time Interval

According to Equations (14) and (23), the optimization formulations depend on the covariances of the historic imbalance signals E [ δ δ T ] . The covariances of the historic imbalance further depend on the time interval that is used in the simulation and even real time operation. Three cases are carried out based on the one-cell system shown in Figure 4. The configurations of the three cases are given in Table 1. The simulation duration is set to 30 min for all three cases. A simulation with time interval of 3 min is introduced in Section 3.2, while the policy results of a simulation with time intervals of 2 min and 5 min are given in Table 1. Then, a comparison is made based on the time intervals that are less and more than 3 min. Combining the results presented in Section 3.2, the cost indices of the three cases are presented in Table 1. It is noted that the cost index decreases with the increased time interval. This is because the calculation of the covariance of historic imbalance is determined by the time interval. Furthermore, the effects of the time interval on power and energy curves are depicted in Figure 9 based on the renewable inputs and load profile in a specific day. Only the battery curves are presented as an example. It is shown that the power and energy curves contain more information when the time interval is smaller. A longer time interval will smooth the prediction errors of the historic data. This will finally reduce the needed reserve and cost. Longer time intervals can reduce the computation time. However, due to the average operation of the prediction errors, time intervals that are too long might result in inaccurate optimization results, i.e., insufficient reserves. On the other hand, shorter time intervals can make the results more reliable, but the computation time will increase accordingly. In real time operation, it is important to choose a proper time interval to make a trade off.
Table 1 presents the comparison of the cost index given by two schemes. It is shown that the cost index of flexible-rate reserve is higher than that of policy-based reserve for all three cases. The cost saving is important because the presented results are obtained from the 30 min simulation. The cost saving of the policy-based reserve will further increase in accordance with the increasing scale and the operation duration of the system.

4. Discussion and Conclusions

Linear decision rule-based control policy for coordinating reserve allocation to the ELECTRA Web-of-Cells system architecture has been developed, implemented, demonstrated, and applied to a three-cell based power system. It has been demonstrated that control policies based on linear decision rules are well suited for integration of distributed energy resources which have flexibility in balancing control, and the robust allocation method is applicable to realistic imbalance signals. Furthermore, it is concluded that a proper time interval selection is important in the linear decision rules-based application. We note that this method relies on several inputs such as: (a) the covariance information of historical data; and (b) the predicted base power production and the prediction error bound of the inelastic power injection, which needs to be improved in the near future. In addition, when adapting the method into real cases, the cost function of elastic power units that represents the balancing service provision should be carefully designed.

Acknowledgments

The research leading to these results has received funding from the European Union Seventh Framework Programme [FP7/2007–2013] under grant agreement No.PADSGCI609687, and includes support from the ELECTRA REX researcher exchange programme. Besides, to accomplish the research, Tian Lan has also received funding from Global Energy Interconnection Research Institute Europe GmbH under the project of “Development of key equipments and engineering technology for large-scale offshore windfarms VSC-HVDC grid connection”.

Author Contributions

Junjie Hu, Tian Lan and Kai Heussen conceived the research idea and designed the simulation setup; Tian Lan and Junjie Hu mainly wrote the paper; amd Kai Heussen, Mattia Marinelli, Alexander Prostejovsky proofread the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Indices

iIndex of inelastic participant.
jIndex of elastic participant.
l,mIndex of cell.

Variable and Parameters

δ Random forecast error vector.
e j Participant j’s nominal elastic power.
r i Nominal prediction of power injection or extraction of participant i.
u j Stacked vector of participant j’s future control inputs.
x j Stacked vector of participant j’s future states.
x 0 j Vector of participant j’s current states.
y i Power injection or extraction of inelastic participant i.
A j Stacked state transition matrix for participant j.
B j Stacked state transition matrix for participant j.
C j Stacked output matrix for participant j.
D j Matrix adjusting power in response to δ .
G i Map from uncertainty to inelastic power injection.
TLength of time horizon in steps.

References

  1. Pineda, I.; Wilkes, J. Wind in Power: 2014, European Statistics; Technical Report; WindEurope: Brussels, Belgium, 2015; pp. 1–12. [Google Scholar]
  2. Agency, I.E. Renewable Energy Medium-Term Market Report 2015; Market Analysis and Forecasts to 2020; Technical Report; International Energy Agency: Paris, France, 2015. [Google Scholar]
  3. Dabra, V.; Paliwal, K.K.; Sharma, P.; Kumar, N. Optimization of photovoltaic power system: A comparative study. Prot. Control Mod. Power Syst. 2017, 2, 3. [Google Scholar] [CrossRef]
  4. Lan, T.; Strunz, K. Multiphysics Transients Modeling of Solid Oxide Fuel Cells: Methodology of Circuit Equivalents and Use in EMTP-Type Power System Simulation. IEEE Trans. Energy Convers. 2017, 32, 1309–1321. [Google Scholar] [CrossRef]
  5. Hu, J.; Yang, G.; Bindner, H.W.; Xue, Y. Application of Network-Constrained Transactive Control to Electric Vehicle Charging for Secure Grid Operation. IEEE Trans. Sustain. Energy 2017, 8, 505–515. [Google Scholar] [CrossRef]
  6. Hiskens, I.; Callaway, D. Achieving controllability of plug-in electric vehicles. In Proceedings of the 2009 IEEE Vehicle Power and Propulsion Conference (VPPC), Dearborn, MI, USA, 7–10 September 2009; IEEE: Piscataway, NJ, USA, 2009; pp. 1215–1220. [Google Scholar]
  7. Knezović, K.; Martinenas, S.; Andersen, P.B.; Zecchino, A.; Marinelli, M. Enhancing the Role of Electric Vehicles in the Power Grid: Field Validation of Multiple Ancillary Services. IEEE Trans. Transp. Electr. 2017, 3, 201–209. [Google Scholar] [CrossRef]
  8. De Martini, P.; Mani Chandy, K.; Fromer, N. Grid 2020 Towards a Policy of Renewable and Distributed Energy Resources; Technical Report; The Resnick Sustainability Institute at Caltech: Pasadena, CA, USA, 2012. [Google Scholar]
  9. Palizban, O.; Kauhaniemi, K.; Guerrero, J.M. Microgrids in active network management—Part I: Hierarchical control, energy storage, virtual power plants, and market participation. Renew. Sustain. Energy Rev. 2014, 36, 428–439. [Google Scholar] [CrossRef]
  10. Bessa, R.J.; Matos, M.A. Economic and technical management of an aggregation agent for electric vehicles: A literature survey. Eur. Trans. Electr. Power 2012, 22, 334–350. [Google Scholar] [CrossRef]
  11. ELECTRA IRP. European Liaison on Electricity Committed towards Research Activity Integrated Research Programm. Available online: http://www.electrairp.eu/ (accessed on 6 February 2018).
  12. Marinelli, M.; Pertl, M.; Rezkalla, M.M.; Kosmecki, M.; Canevese, S.; Obushevs, A.; Morch, A.Z. The Pan-European Reference Grid Developed in the ELECTRA Project for Deriving Innovative Observability Concepts in the Web-of-Cells Framework. In Proceedings of the 51st International Universities Power Engineering Conference (UPEC), Coimbra, Portugal, 6–9 September 2016; IEEE: Piscataway, NJ, USA, 2016. [Google Scholar]
  13. D’hulst, R.; Fernández, J.M.; Rikos, E.; Kolodziej, D.; Heussen, K.; Geibelk, D.; Temiz, A.; Caerts, C. Voltage and frequency control for future power systems: The ELECTRA IRP proposal. In Proceedings of the 2015 International Symposium on Smart Electric Distribution Systems and Technologies (EDST), Vienna, Austria, 8–11 September 2015; pp. 245–250. [Google Scholar]
  14. Morch, A.Z.; Jakobsen, S.H.; Visscher, K.; Marinelli, M. Future control architecture and emerging observability needs. In Proceedings of the 2015 IEEE 5th International Conference on Power Engineering, Energy and Electrical Drives (POWERENG), Riga, Latvia, 11–13 May 2015; pp. 234–238. [Google Scholar]
  15. Warrington, J.; Goulart, P.; Mariéthoz, S.; Morari, M. Policy-based reserves for power systems. IEEE Trans. Power Syst. 2013, 28, 4427–4437. [Google Scholar] [CrossRef]
  16. Jabr, R.A. Adjustable Robust OPF With Renewable Energy Sources. IEEE Trans. Power Syst. 2013, 28, 4742–4751. [Google Scholar] [CrossRef]
  17. Bertsimas, D.; Brown, D.B.; Caramanis, C. Theory and applications of robust optimization. SIAM Rev. 2011, 53, 464–501. [Google Scholar] [CrossRef]
  18. Hu, J.; Heussen, K.; Claessens, B.; Sun, L.; D’Hulst, R. Toward coordinated robust allocation of reserve policies for a cell-based power system. In Proceedings of the 2016 IEEE PES Innovative Smart Grid Technologies Conference Europe (ISGT-Europe), Ljubljana, Slovenia, 9–12 October 2016; pp. 1–6. [Google Scholar]
  19. Caerts, C.; Rikos, E.; Syed, M.; Guillo Sansano, E.; Merino-Fernández, J.; Rodriguez Seco, E.; Evenblij, B.; Rezkalla, M.M.N.; Kosmecki, M.; Temiz, A.; et al. Description of the Detailed Functional Architecture of the Frequency and Voltage Control Solution (Functional and Information Layer); Technical Report; Technical University of Denmark: Lyngby, Denmark, 2017. [Google Scholar]
  20. Stoft, S. Power System Economics: Designing Markets for Electricity, 1st ed.; Wiley-IEEE Press: New York, NY, USA, 2002. [Google Scholar]
  21. Brunner, H.; Tornelli, C.; Cabiati, M. The Web of Cells Control Architecture for Operating Future Power Systems; Technical Report; ELECTRA IRP—Public Deliverables: WP5 Increased Observability; European Union: Brussels, Belgium, 2018; under review. [Google Scholar]
  22. Lofberg, J.Y. YALMIP: A toolbox for modeling and optimization in MATLAB. In Proceedings of the International Symposium on Computer Aided Control Systems Design (CACSD), New Orleans, LA, USA, 2–4 September 2004; IEEE: Piscataway, NJ, USA, 2004; pp. 284–289. [Google Scholar]
  23. Optimization, Gurobi. Inc. Gurobi Optimizer Reference Manual; Technical Report; Gurobi Optimization, Inc.: Houston, TX, USA, 2017. [Google Scholar]
  24. Oliver, G.; Bindner, H. Building a test platform for agents in power system control: Experience from SYSLAB. In Proceedings of the International Conference on Intelligent Systems Applications to Power Systems (ISAP), Toki Messe, Niigata, Japan, 5–8 November 2007; IEEE: Piscataway, NJ, USA, 2007; pp. 1–5. [Google Scholar]
  25. Prostejovsky, A.M.; Marinelli, M.; Rezkalla, M.; Syed, M.H.; Guillo-Sansano, E. Tuningless Load Frequency Control Through Active Engagement of Distributed Resources. IEEE Trans. Power Syst. 2017. [Google Scholar] [CrossRef]
  26. Energinet. Ancillary Services to Be Delivered in Denmark Tender Conditions; Technical Report; Energinet: Erritsø, Denmark, 2017. [Google Scholar]
Figure 1. An example of a three-cell small scale power system.
Figure 1. An example of a three-cell small scale power system.
Energies 11 00381 g001
Figure 2. Control system architecture.
Figure 2. Control system architecture.
Energies 11 00381 g002
Figure 3. SYSLAB layouts for the two test cases. Devices and lines used in the one-cell case are marked in red. Additional components for the three-cell case are indicated in orange.
Figure 3. SYSLAB layouts for the two test cases. Devices and lines used in the one-cell case are marked in red. Additional components for the three-cell case are indicated in orange.
Energies 11 00381 g003
Figure 4. One-cell based system.
Figure 4. One-cell based system.
Energies 11 00381 g004
Figure 5. Policy of reserve allocation for one-cell system with time interval of 3 min.
Figure 5. Policy of reserve allocation for one-cell system with time interval of 3 min.
Energies 11 00381 g005
Figure 6. Three-cell based system.
Figure 6. Three-cell based system.
Energies 11 00381 g006
Figure 7. Policy of reserve allocation for three-cell power system in Scenario 1.
Figure 7. Policy of reserve allocation for three-cell power system in Scenario 1.
Energies 11 00381 g007
Figure 8. Policy of reserve allocation for three-cell power system in Scenario 2.
Figure 8. Policy of reserve allocation for three-cell power system in Scenario 2.
Energies 11 00381 g008
Figure 9. Energy and power curves of battery.
Figure 9. Energy and power curves of battery.
Energies 11 00381 g009
Table 1. Comparison of cost index regarding different time intervals.
Table 1. Comparison of cost index regarding different time intervals.
CaseSimulation DurationTime IntervalHorizon LengthCost Index of Policy-Based ReserveCost Index of Flexible-Rate Reserve
130 min2 min1530.5130.62
230 min3 min1030.2030.29
330 min5 min629.9730.04
Table 2. Properties of devices used in the one-cell system.
Table 2. Properties of devices used in the one-cell system.
DeviceTest Case P nom (kW) P min (kW) P max (kW)Description
Solar1, 3 Cell10.10.010.1Orientation
az. 180 , el. 40
Battery1, 3 Cell0.0−15.015Vanadium redox flow type
190 kWh, initial state of charge is 50%
EV1, 3 Cell0.0−2.02.0Bidirectional charger
20 kWh, initial state of charge is 50%
Mob. Load1, 3 Cell−33.0−33.00.0Thyristor-contr.
Aircon3-Cell9.80.09.8Wind turbine

Share and Cite

MDPI and ACS Style

Hu, J.; Lan, T.; Heussen, K.; Marinelli, M.; Prostejovsky, A.; Lei, X. Robust Allocation of Reserve Policies for a Multiple-Cell Based Power System. Energies 2018, 11, 381. https://doi.org/10.3390/en11020381

AMA Style

Hu J, Lan T, Heussen K, Marinelli M, Prostejovsky A, Lei X. Robust Allocation of Reserve Policies for a Multiple-Cell Based Power System. Energies. 2018; 11(2):381. https://doi.org/10.3390/en11020381

Chicago/Turabian Style

Hu, Junjie, Tian Lan, Kai Heussen, Mattia Marinelli, Alexander Prostejovsky, and Xianzhang Lei. 2018. "Robust Allocation of Reserve Policies for a Multiple-Cell Based Power System" Energies 11, no. 2: 381. https://doi.org/10.3390/en11020381

APA Style

Hu, J., Lan, T., Heussen, K., Marinelli, M., Prostejovsky, A., & Lei, X. (2018). Robust Allocation of Reserve Policies for a Multiple-Cell Based Power System. Energies, 11(2), 381. https://doi.org/10.3390/en11020381

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