Next Article in Journal
Effect of Wind Turbine Blade Rotation on Triggering Lightning: An Experimental Study
Next Article in Special Issue
Multi-Agent System Fault Protection with Topology Identification in Microgrids
Previous Article in Journal
The Desalination Process Driven by Wave Energy: A Challenge for the Future
Previous Article in Special Issue
Mitigation of the Impact of High Plug-in Electric Vehicle Penetration on Residential Distribution Grid Using Smart Charging Strategies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Privacy-Preserving Distributed Optimal Scheduling for Interconnected Microgrids

State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources, North China Electric Power University, Beijing 102206, China
*
Author to whom correspondence should be addressed.
Energies 2016, 9(12), 1031; https://doi.org/10.3390/en9121031
Submission received: 20 October 2016 / Revised: 21 November 2016 / Accepted: 29 November 2016 / Published: 7 December 2016
(This article belongs to the Special Issue Smart Microgrids: Developing the Intelligent Power Grid of Tomorrow)

Abstract

:
With the development of microgrids (MGs), interconnected operation of multiple MGs is becoming a promising strategy for the smart grid. In this paper, a privacy-preserving distributed optimal scheduling method is proposed for the interconnected microgrids (IMG) with a battery energy storage system (BESS) and renewable energy resources (RESs). The optimal scheduling problem is modeled to minimize the coalitional operation cost of the IMG, including the fuel cost of conventional distributed generators and the life loss cost of BESSs. By using the framework of the alternating direction method of multipliers (ADMM), a distributed optimal scheduling model and an iteration solution algorithm for the IMG is introduced; only the expected exchanging power (EEP) of each MG is required during the iterations. Furthermore, a privacy-preserving strategy for the sharing of the EEP among MGs is designed to work with the mechanism of the distributed algorithm. According to the security analysis, the EEP can be delivered in a cooperative and privacy-preserving way. A case study and numerical results are given in terms of the convergence of the algorithm, the comparison of the costs and the implementation efficiency.

1. Introduction

Microgrids (MGs) are self-controlled entities, which facilitate the penetration of renewable energy and distributed energy resources (DERs) for economic and reliability purposes. Generally, the MG can be operated in either grid-connected [1] or island mode [2]. With the development of MGs, a new typical approach is to consider where several MGs exchange energy with one another even when the MGs are isolated from the utility grid. In other words, there exist energy flows within a group of physically-connected MGs, but not between the MGs and the main grid. One promising application of interconnected MGs (IMG) is for the service restoration in a self-healing distribution system to improve the dynamic performance when islanding occurs and to extend the power supply during the system outages [3,4]. Another important application is to ensure the full utilization of renewable energy resources (RESs) and improve the economic benefit and reliability in isolated areas [5]. Consequently, the interconnected operation of multi-MGs is becoming a promising strategy for the smart grid.
Considering the type of MGs, the interconnection among autonomous MGs can be realized via an AC [6,7], DC [8] or AC-DC hybrid architecture [9]. As a newly-emerging topic of the smart grid, there are still limited research works focused on the related problems of IMGs. The studies can be roughly divided into two categories: (1) coordinated control and (2) energy management. For the coordinated control, the feasibility (benefits) and basic controls of IMGs via AC connections were studied in [6,7]. Moreover, in [5,9,10], considering that the AC MGs may be operated at different voltages and frequencies, hybrid AC/DC links and coordinated control schemes were studied for maintaining the autonomy of each MG during interconnections. For the interconnection of DC MGs, several control methods were proposed, such as the control strategy for two DC MGs [11], distributed hierarchical control [12] and distributed two-level tertiary control [13] for a cluster of DC MGs.
Energy management is another category of study for the IMG. To ensure the full utilization of RES and to minimize the operation cost, effective scheduling of the DERs and battery energy storage systems (BESS) [14] is highly necessary. Generally, if all of the MGs could share the information on their respective load data, generation data and parameters and grid parameters, the optimal scheduling could be easily realized based on the traditional centralized optimization [15], such as optimal power flow (OPF). However, for security considerations, it is not desirable for each MG to do so because the shared information could compromise the privacy of each MG. Thus, this is the basic driving force of the distributed optimization of IMG scheduling. For instance, a distributed convex optimization framework was developed for energy trading between islanded MGs [16]. The problem consists of several islanded MGs that exchange energy flows by means of an arbitrary topology. In [17], the authors studied the self-organization and decentralized energy management of an MG cluster islanded from the main grid after a disruptive event.
However, even with applying distributed optimization, there are still weak points in these works. First, for privacy preserving, the algorithms require the MGs to share enough data with others during the optimization. For instance, the trading prices and corresponding trading energy should be shared between the neighboring MGs in [16], and the expected power exchange of each MG should be shared with all other MGs in [17]. After long-term operation, the generation details and energy consumption patterns of one MG can be gradually inferred by other MGs using data mining [18], which could gradually affect the confidentiality and privacy [19]. In the smart grid, data perturbation and aggregation are the two main privacy-preserving approaches that have been applied to avoid the inference of sensitive information from the metering data [20]. For privacy preserving, the assumed attacks are mainly based on non-intrusive monitoring of energy usage traces. Typical solutions rely on battery-based load hiding [21,22], on noise injection [23] or on a multi-party computation cryptosystem [24,25]. However, in the distributed optimization framework, the shared data are tightly coupled with the local optimization and circulated during the execution of the scheduling horizon. In this regard, the solutions of battery-based hiding and noise injection are difficult to apply. Second, the modeling of BESS is not considered in the optimization of IMG. For the MGs consisting of wind, solar and diesel generations, the BESS is an important component to increase the utilization of RES [26,27]. With suitable sizing [26,28] and coordinated control [29], BESS could have great influence on the optimal scheduling of MGs.
To this end, this paper proposes a privacy-preserving distributed optimal scheduling method for energy management of IMG, with the consideration of BESS and RES. The main contributions of the paper are as follows: (1) based on the relationship between charging-discharging features and the battery life cost, the cost model of BESS is formulated as a quadratic polynomial; (2) by using the framework of the alternating direction method of multipliers (ADMM), a distributed optimal scheduling model and iterative algorithm is proposed; (3) based on the multi-party computation cryptosystem, a privacy-preserving strategy to protect the sharing of expected exchanging power (EEP) between MGs is designed to work with the iterative algorithm.

2. System Model

In this section, we describe the system model for the distributed optimal scheduling of IMG. In the MG, we take RESs, conventional distributed generation, and BESS into consideration. RESs, such as wind turbines (WT) and photovoltaics (PV), are non-dispatchable, and conventional distributed generation such as diesel, is dispatchable. BESS can be regarded as dispatchable resources for discharging and controllable load for charging.

2.1. Renewable Energy Resources

A RES unit, such as WT and PV, is not dispatchable, and its output power depends on the availability of the primary sources (i.e., solar irradiance or wind speed). Therefore, forecasting is required in order to consider them in the optimal operation of IMG [30]. In the optimization, we can assume that the complex power of an RES over the scheduling horizon is given, and there is no generation cost for RES [31].

2.2. Diesel Generation

Diesel generation (DG) is a dispatchable source, which plays the role of a supplement source when there is a shortage of RES power. Diesel fuel is necessary for generation, and the relationship between cost and output power can be expressed as:
C g e n ( P g e n ) = c f u e l · a g P g e n + b g P g e n 2 + c g
where c f u e l is the unit price of diesel fuel (CYN/L); P g e n is the output power of the diesel generation; a g , b g , c g are the coefficients of the fuel demand of the generation to its output power [32].

2.3. BESS

The lifetime of the battery will be depleted after a certain number of charging and discharging cycles. The relationship of the BESS’ lifetime loss can be described as:
L l o s s = A c A t o t a l
where L l o s s is the life loss percentage of BESS; A c is the effective cumulative capacity in a certain period, kWh. A t o t a l is the total effective cumulative capacity of the full life cycles, kWh. According to the manufacturing manual of the lead acid battery, it will deliver 390 Q effective kWh over its lifetime with the capacity of Q [33].
Furthermore, A c is decided by the battery’s state of charge (SOC) and the actual capacity A c , and the relationship is described as:
A c = λ S O C · A c
where λ S O C is an effective weight. When the SOC of BESS takes the value from 0 . 5 to 1 . 0 , we can approximately assume that λ S O C has a linear relation with SOC [33].
λ S O C = a · S O C + b
where a and b are two coefficients.
When BESS’ SOC changes from S O C i n i t to S O C e n d during the operation of discharging, the effective capacity is described as:
A c = S O C e n d S O C i n i t λ S O C Q d S O C
With the consideration of the relationship between discharging power P and its duration T , S O C e n d can be further described as:
S O C e n d = S O C i n i t P · T / Q
and the equivalent accumulative capacity is
A c = a 2 Q P 2 T 2 + P T ( a · S O C i n i t + b )
Integrating (2) and (7), the lifetime cost can be described as a function of P:
C B S ( P ) = I a P 2 T 2 + 2 P T Q ( a S O C i n i t + b ) 2 Q A t o t a l
where I is the investment of BESS.

2.4. Optimal Scheduling Model of MGs

Given the forecasted power of WT, PV and loads, the MG operator can schedule the dispatchable resources’ output power to minimize the total operation cost. Generally, the MG-n’s optimal scheduling problem can be formulated as:
m i n h n ( x n ) = f n ( x n ) + γ ( x n i n ) = C D G n ( x D G n ) + C B S n ( x B S n ) + γ ( x n i n )
s . t . d n P u c n x D G n x B S n = x n i n
0 x D G n P D G n m a x
P c h a r g e m a x x B S n P d i s c h a r g e m a x
S O C m i n S O C n S O C m a x
where x n = [ x D G , x B S ] is composed by the dispatchable sources’ output power, in which x D G n and x B S n are DG and BESS’s output power, respectively. x n i n is the expected exchange power of MG-n. If MG-n imports power from IMG, x n i n is positive; otherwise, if MG-n exports power to IMG, x n i n is negative. γ ( · ) is the cost of trading x n in the IMG, which can be charged by the network investor considering the investment of power lines and network losses. d n is the MG-n’s total power demand, and P u c n is the total non-dispatchable output power. Considering that the energy efficiency of BESS in charging/discharging is less than 100%, when the BESS’ scheduling power is x B S n , the actual charging power should be x B S n · η , whereas the actual discharging power should be x B S n / η . As shown in (1) and (8), the cost function of DG and BS is convex, which indicates that the cost function h n ( x ) is convex, as well.
As for the constraints, (10) is set to keep the power balance of MG; (11) ensures that the output power of DG is in the range of rated capacity; (12) is used to limit the output power of the BESS, according to the fact that the charging/discharging power should not exceed the maximum power of the bi-directional inverter; (13) sets the maximum and minimum SOC of BESS in operation.

2.5. Optimal Scheduling Model of the IMG

In the environment of interconnected operation, all participating MGs are trying to lower their total operation cost by considering the cost efficiency of DERs and scheduling in a cooperative way. For an IMG containing N MGs, its optimization problem can be formulated as [16]:
m i n n = 1 N f n ( x n ) + γ n ( x n ) s . t . n = 1 N 1 T x n = n = 1 N ( d n P u c n )
where 1 T is the unit row vector.
It is noted that in (14), all of the participating MGs are assumed to work as a coalition [34], in which they will aim to minimize the global operation cost of the IMG.
In general, there are two methods to solve this optimal scheduling problem: centralized optimization and distributed optimization. In the centralized method, all private information of each MG needs to be collected by the control center (CC), and the optimal scheduling result will be distributed to each MG after optimization, as shown in Figure 1.

3. Distributed Optimization

3.1. ADMM

ADMM was proposed by Gabay in the 1970s and has been widely used in large-scale distributed computing and optimization [35]. With the excellent performance in convergence and robustness, the method has drawn great attention in recent years. In the field of power system optimization, ADMM has also been used to solve the distributed optimization problem [36,37]. The basic form can be described as:
m i n s ( x ) + g ( z ) s . t . A x + B z = c
where x R n ; z R m ; A R p × n ; B R p × m ; c R p . Functions s and g are all convex. When function s is R n R { + } and g is R m R { + } , the algorithm can converge to the optimized solution [35].

3.2. Decentralizing the Problem

To decentralize objective Function (14), we define the following function according to the standard ADMM form [35]:
g ( z ) = 0 , z R N + z R N C = { z R N | z 1 + z 2 + + z N = 0 }
To make the equation simple, we can express x i n as:
x n i n = 1 T x n d n P u c n , n = [ 1 , 2 , , N ]
Combined with (9), Function (14) can be further described as:
m i n n = 1 N h n ( x n ) + g ( z ) s . t . x i n z = 0
where x i n = [ x 1 i n , x 1 i n , , x N i n ] indicates the EEP vector of MGs in IMG. Through variable z , objective function (18) can be solved in a distributed manner.
According to the standard iteration form of ADMM, the problem’s iteration process can be described as:
x n k + 1 = arg min x h n ( x ) + ρ 2 x n i n z n k + u n k 2 2 z k + 1 = arg min z g ( z ) + ρ 2 n N z n x n i n k + 1 u n k 2 2 u n k + 1 = u n k + x n i n k + 1 z n k + 1
where ρ > 0 is a given penalty coefficient; k is the number of iterations; u n is the Lagrange multiplier for the equation constraint of (18). In the process of updating z , all x n i n should be collected. To enhance the privacy-preserving ability and simplify the iteration, we can have:
z ¯ = 1 N n = 1 N z n
and then, (16) can be transformed as:
g ˜ ( z ¯ ) = 0 , z ¯ Q + , z ¯ Q Q = { z ¯ | z ¯ = 0 }
For any z , its update in (19) is equivalent to:
min g ˜ ( z ¯ ) + ρ 2 n = 1 N z n x n i n k + 1 u n k 2 2 s . t . z ¯ = 1 N n = 1 N z n
For any constant z ¯ , the Lagrange equation of (22) can be described as:
L ( z , μ ) = g ˜ ( z ¯ ) + ρ 2 n = 1 N z n x n i n k + 1 u n k 2 2 + μ · ( z ¯ 1 N n = 1 N z n )
where μ is the Lagrange multiplier for the constraint in (22). As (22) is a convex function, the necessary and sufficient condition is:
z L ( z , μ ) = 0 a n d z ¯ = 1 N n = 1 N z n
For any constant z ¯ , g ˜ ( z ¯ ) is a scalar, and the condition z L ( z , μ ) = 0 translates to:
ρ ( z n x n i n k + 1 u n k ) 1 N μ = 0 , n = [ 1 , 2 , , N ]
After some manipulation, (24) and (25) collectively imply:
z n x n i n k + 1 u n k = z ¯ x ¯ i n k + 1 u ¯ k
where we have used the same overline-notation for x i n k + 1 and u k :
x ¯ i n k + 1 = 1 N n = 1 N x n i n k + 1 , u ¯ k = 1 N n = 1 N u n k + 1
By employing (26) in (22), the z -update introduced in (19) can be equivalently formulated as:
z ¯ k + 1 = arg min z ¯ g ˜ z ¯ + N · ρ 2 z ¯ x ¯ i n k + 1 u ¯ k 2 2
After being simplified, it is only needed to update z ¯ , and the Lagrange multiplier u n can be expressed as:
u k + 1 = u k + x ¯ i n k + 1 z ¯ k + 1
Combining (28) and (29), replace z n k in the updating process of x n , and then, the simplified distributed iteration form is:
x n k + 1 = arg min x h n ( x ) + ρ 2 x n i n x n i n k + x ¯ i n k z ¯ k + u k 2 2 z ¯ k + 1 = arg min z ¯ g ˜ ( z ¯ ) + N · ρ 2 z ¯ x ¯ i n k + 1 u k 2 2 u k + 1 = u k + x ¯ i n k + 1 z ¯ k + 1
According to (21), the result of update z ¯ would be:
z ¯ k + 1 = 0
The final form of the distributed iteration can be described as:
x n k + 1 = arg min x h n ( x ) + ρ 2 x n i n x n i n k + x ¯ i n k + u k 2 2 u k + 1 = u k + x ¯ i n k + 1
From (32), we can know that only the average EEP in IMG is required in the optimization, and the Lagrange multiplier’s updating can be calculated within each MG.
As shown in Figure 2, by using the distributed optimization, each MG solves its own optimal scheduling and can finally reach the global optimal solution by informing its EEP to other MGs during the iterations.
Furthermore, the dual residual can be used as the stopping criteria for the iteration convergence.
r k = 1 N n = 1 N x n i n k
When the residual satisfies the following conditions, it is considered to achieve convergence, and the iteration is stopped.
r k 2 2 ε
in which ε is the convergence error. The calculation process of the distribution optimal scheduling is described in Algorithm 1.
Algorithm 1 Distributed scheduling algorithm
 1:  procedure DISPATCHING( x n , n = [ 1 , 2 , , N ] )    
 2:   Initiate u,ρ,k
 3:   while r 2 2 > ε do
 4:      x n k + 1 = arg min x h n ( x ) + ρ 2 x n i n x n i n k + x ¯ i n k + u k 2 2
 5:      x ¯ i n k + 1 = 1 N n = 1 N x n i n k + 1
 6:      u k + 1 = u k + x ¯ i n k + 1
 7:      r k + 1 = 1 N n = 1 N x n i n k + 1
 8:      k = k + 1
 9:   end while
10:  end procedure

4. Privacy-Preserving Strategy

As shown in (27) and Algorithm 1, each MG should share its own EEP. After that, the average EEP will be calculated and sent back to all of the participating MGs. Considering the privacy, each MG should share its own EEP, but does not want to show it to the other MGs. By coincidence, as the Paillier cryptosystem has the features of homomorphic addition, indistinguishability and self-blinding, it is feasible to transfer the addition of plaintext EEPs to the multiplication of ciphertext EEPs, which makes it quite suitable for privacy preserving in the sharing of EEPs [38].

4.1. Basic Theory of the Paillier Cryptosystem

The Paillier cryptosystem is widely used in a distributed environment and can protect users’ privacy in the collaboration process [39]. The basic schema of Paillier includes three parts, which are described as follows [38].

4.1.1. Key Generation

Select two prime numbers p and q, and calculate N k = p q and λ k :
λ k = l c m ( p 1 , q 1 )
where l c m means the operation of calculating the lowest common multiple. We denote Euler’s totient function ϕ ( n ) taken on N k , i.e., ϕ ( N k ) = ( p 1 ) ( q 1 ) in the present case. Recall that | Z N k 2 * | = ϕ ( N k 2 ) , then a random number g k Z N k 2 * is needed and should satisfy:
g c d ( L ( g k λ k m o d N k 2 , N k ) ) = 1
where L ( x ) = ( x 1 ) / N k and g c d is the operation of calculating the greatest common divisor. The Paillier public key and private key are < N k , g k > and λ k .

4.1.2. Encryption

Assume that m c Z N is plaintext, w Z N is a random number and E ( · ) is the encryption operation.
E ( m c m o d N k , w m o d N k ) = g k m c w N k m o d N k 2

4.1.3. Decryption

Assume that c Z N k 2 * is ciphertext and D ( · ) is the decryption operation with the private key λ k .
D ( c ) = L ( c λ k m o d N k 2 ) L ( g k λ k m o d N k 2 ) m o d N k
To secure the Paillier cryptosystem’s privacy, parameter N k and g k should be big enough. m c 1 , m c 2 , w 1 , w 2 Z N , the formula holds:
E ( m c 1 , w 1 ) E ( m c 2 , w 2 ) = E ( m c 1 + m c 2 , w 1 w 2 ) m o d N k 2
which indicates that the result of the multiple operation of two ciphertexts is equal to the encryption result of the sum of two plaintexts, and this is Paillier’s homomorphic addition characterization.
As can be seen, the plaintext should be an integer in the process of encryption, but the information exchanged might be decimals in actual applications. To make this cryptosystem useful, a constant multiple is used to amplify the plaintext.

4.2. Protocol of EEP Sharing

Firstly, an authority entity for the IMG is required, which is denoted as A U . The public key and private key of the A U are denoted as < N A U , g A U > and λ A U . The public key < N A U , g A U > is known to every MG participating in the IMG, and the private key λ A U is only known by the A U itself.
Secondly, the sharing for the EEP of every MG is implemented in a series manner. The route of information delivery must be fixed and known to every MG for each iteration.
As we know, x n i n is the private data of MG n, and its ciphertext encrypted by the public key of A U can be denoted as E ( x n i n ) , which is transmitted from MG n to n + 1 . After each step, A U is the last receiver of the ciphertext and uses its private key to decrypt and broadcast the data to all of the MGs. The detailed information exchange procedure is listed as follows (see Figure 3).
Step 1: Every MG optimizes its own EEP x n i n , encrypts the plaintext by using A U ’s public key and gets the ciphertext E ( x n i n ) .
Step 2: MG 1 transmits E ( x 1 i n ) to MG 2, and then MG 2 should calculate as:
E ( x 1 i n + x 2 i n ) = E ( x 1 i n ) E ( x 2 i n ) .
After calculation, E ( x 1 i n + x 2 i n ) is transferred to the next MG and does the same operation until the last MG.
E ( n = 1 N x n ) = n = 1 N E ( x n i n )
Step 3: The last MG N sends the cipher data to A U , and A U has the private key to perform the decryption process. Finally, the average EEP for distributed optimization is broadcast to all MGs.
x ¯ i n = 1 N n = 1 N x n i n

4.3. Security Analysis of the Protocol

We use the semi-honest security model to analyze the protocol [40]. This model assumes that the involved parties, in our case MG 1 to MG N and the A U , perform the protocol steps properly, but are also allowed to store the messages from previous steps in an attempt to deduce more information than they are entitled to. Based on the semi-honest model, we can enumerate three different attach scenarios: (1) all of the MGs are trustful, but the AU is untrustful; (2) some of the MGs are untrustful and collusive, but the AU is trustful; (3) some of the MGs are untrustful and collusive, and the AU is also untrustful.
  • For the first scenario, the sensitive data of each MG x n i n are encrypted and delivered on a fixed path. The AU can receive the ciphertext of summing EEP E ( n = 1 N x n ) . Although the plaintext can be decrypted as n = 1 N x n , the actual data of arbitrary MGs cannot be deduced.
  • For the second scenario, as the EEP of each MG is encrypted by the public key of AU E ( x n i n ) , even if there are a number of collusive MGs, the plaintext of any MGs cannot be obtained without the private key of AU.
  • For the third scenario, we assume that there are N c collusive MGs. (a) When N c = N 1 , that is only one non-collusive MG exists (denoted as MG-a), then, the MG-a is the target of the collusive MGs. In this case, if all of the collusive MGs can share their own plaintexts with the AU or set their EEPs to zero, the AU can get the privacy data of MG-a. Otherwise, the privacy of MG-a still can be guaranteed. However, this condition is difficult to achieve, since it would also affect the privacy of the collusive MGs or make it easy to be detected by MG-a. (b) When N c < N 1 , that is the number of non-collusive MGs is not less than two. Similar to the condition of N c = N 1 , the untrustful AU can only get the sum of EEPs from the non-collusive MGs (e.g. n = 1 N N c x n ) even with the worst conditions.
Therefore, the privacy-preserving protocol can guarantee the privacy of every MG under the check of the semi-honest model.

5. Case Study

5.1. Basic Data

As shown in Figure 4, an IMG composed of three MGs is used to test the proposed distributed scheduling method. The parameters of DERs and BESS are listed in Table 1 and Table 2, respectively.
In the case study, 15 min is set as the length of a time slot. That is, there are 96 time slots in a day. The optimization is solved in a distributed manner by using the short-term forecasted power of PVs, WTs and loads, and only the encrypted information of the EEP should be delivered during the process of optimization. The net load curves of the MGs and the total net load curve of the IMG are shown in Figure 5.
Apparently, these MGs have different net load profiles. The IMG has negative net load in Time Slots 1–17, 33–61 and 91–96. In these time slots, BESSs need to store the excess renewable energy. While in other time slots, BESSs or DGs should discharge or generate power to satisfy the load demand, considering their operation costs.

5.2. Result and Analysis of the Distributed Optimal Scheduling

The convergence of the distributed algorithm in each time slot is shown in Figure 6. The algorithm can achieve convergence in 40–50 iterations for almost every time slot, which indicates that the ADMM has good performance on solving this problem. Moreover, to analyze the evolution of the internal variables of three MGs, the 26th time slot in the scheduling horizon is chosen as an example, as shown in Figure 7. During the iteration process, the results show that each MG has decided its own scheduling scheme with the interaction of the shared EEP values, and all of the MGs finally approach the global optima without any changes for their solutions.
The output power curves of all of the controllable DERs are shown in Figure 8, and the curves of power exchanged via the connecting tie lines are shown in Figure 9. In Time Slots 1–17, 33–61 and 91–96, the total power of PVs and WTs is much higher than the load demand, so the BESSs are working in charging mode to store the surplus RES energy (the values are negative in Figure 8).
When the PVs and WTs cannot satisfy the load demand, the BESSs and DGs are scheduled to discharge or generate power and, thus, secure the IMG’s operation. The output power of each controllable energy resource is decided by the cost characteristics. In Figure 8, we can find out that the energy of BESSs is utilized as first priority when there is a shortage of RES (i.e., Time Slots 15–25 and 61–71). After that, the DGs are scheduled to generate the remaining shortages (i.e., Time Slots 25–35 and 71–91). Furthermore, Figure 9 shows the curves of the decided exchanging power of different tie lines in the IMG. Combined with Figure 5, Figure 8 and Figure 9, the results show that the distributed optimal scheduling is reasonable for the operation of the IMG.

5.3. Comparison with the Centralized Optimization

The operation costs of the IMG are respectively obtained by using the centralized optimization and the distributed optimization, shown in Figure 10 and Table 3. The results show that the distributed optimal scheduling method is approximately the same as the centralized optimization; there are only minor errors in several time slots, which are mainly because of the rounding errors in the calculation.

5.4. Comparison of Operation Cost between the Isolated Mode and Interconnected Mode

Assuming that all of the MGs can be operated in the isolated or interconnected mode, the operation cost of both modes are given in Table 4.
In the isolated operation mode, the results show that there is discarded RES energy in MG1 and MG3. Due to the low load requirements and limited BESS in these two MGs, the RES energy cannot be fully utilized. Finally, the discarded energy has a great influence on the operation cost of MGs. As for the comparison, in the interconnected mode, the results show that the MGs are highly complementary in RES energy consumption, which has significantly decreased the overall cost of IMG.

5.5. Operation Cost Analysis for the Different Configured Capacities of BESS

Assuming that the BESS capacity has been changed from 0.7–1.2-times the original capacity in Table 1, the sensitivity analysis of the overall operation costs with different capacities of the BESS in MGs is shown in Figure 11. The results show that the cost change tendency can be categorized as two stages. First, from 0.7–0.9, the operation cost of IMG is reduced in this stage. This is mainly due to the discarded power of RESs having been decreased with the increasing of the BESS capacity, which eventually saves the cost of DGs. Second, from 0.9–1.2, the operation cost of IMG is increased with the increasing of the BESS capacity. Considering that no more discarded RES power can be saved, the overall cost of IMG has been increased with the investment of BESS.

5.6. Analysis for the Impact of Forecasting Errors

In the optimal scheduling, the solutions are obtained based on the short-term forecasting of RES power and load. In fact, the forecasting errors may cause unbalance between the generations and loads. Thus, in the interconnected operation, we assume that all of the MGs are working in the fixed power control of interconnection tie-lines. In this regard, all of the forecasting errors of load and RES could result in the autonomous power regulation of BESS or diesel generators in each MG, to keep the power balance and possible cost deviations by the MG itself. The detail of this control model can be found in [15]. Moreover, to analyze the influence of operation cost caused by the uncertain forecasting errors, we assume that the forecasting errors are randomly located in [−10%, 10%] [30]. By using the control model [15], the final operation cost with balance regulations is shown in Table 5.
The results show that the forecasting errors have very limited influence on the total cost of IMG, in which the costs of MG1and MG3 are reduced and the cost of MG3 is increased. It is apparent that the forecasting errors have the largest impact on the cost of MG1. Considering that MG1 has the largest capacity of RES, this is reasonable in the operation.

5.7. Efficiency Analysis of the Privacy-Preserving Protocol

As the privacy-preserving strategy is implemented as a series procedure in a wide area distributed environment, it is important to test the time requirement to finish the whole process of the distributed optimization. Therefore, we construct an experimental platform consisting of several computers and 3G communication networks. After complete implementation of one day’s 96 time slot optimization, the execution time for each procedure is listed in Table 6, including the objection function solution, encryption, communication and decryption.
As shown in Figure 6, almost all of the scheduling of every time slot can be converged in around 50 iterations. Combined with the privacy-preserving protocol and the results of Table 6, the time requirement for one scheduling time slot can be finished in around 10 min. The result means that the minimum length of one time slot for the scheduling should be more than 10 min. Therefore, the method is suitable for the IMG’s optimal scheduling problem of day ahead, hour ahead and even 15 min ahead.

6. Conclusions

A privacy-preserving distributed scheduling method for IMG with battery storage and renewable energy resources is proposed. According to the case study and experiment results, the method has almost the same optimization result with the centralized optimization, and the privacy of each MG can also be guaranteed. The processes of the iteration algorithm and privacy-preserving have delayed the implementation time of optimal scheduling, but still can be accepted for the condition with the minimum time slot longer than 10 min, which is suitable for the optimal scheduling problem of IMG.
The proposed method can be extended in several ways. Firstly, one interesting problem is the involvement of flexible loads (i.e., shiftable appliances, controllable heating, ventilation and air conditioning loads, etc.). In this scenario, the MG operator can choose to adjust the load, the generation of DGs or trade with other MGs, with the comprehensive consideration of the generation cost, trading price, load characteristic and demand response cost. Secondly, another interesting problem is the time horizon of optimization. In energy management that includes multiple time periods, the intertemporal constraint for remaining energy of BESS should be further considered. This constraint may have an influence on the distributed model, which needs further study. Finally, as the participating MGs are assumed to work as a coalition, the basic principle for forming a coalition and the related re-distribution strategy of cooperative benefits by using the IMG still need further research.

Acknowledgments

This work was supported in part by the National High-tech R&D Program of China (863 Program) (No. 2014AA052001) and in part by the Fundamental Research Funds for the Central Universities (No. 2015ZD02).

Author Contributions

The paper was a collaborative effort between the authors. Nian Liu, Cheng Wang, Minyang Cheng and Jie Wang contributed collectively to the theoretical analysis, modeling, simulation and manuscript preparation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, N.; Chen, Q.; Liu, J.; Lu, X. A Heuristic Operation Strategy for Commercial Building Microgrids Containing EVs and PV System. IEEE Trans. Ind. Electron. 2015, 62, 2560–2570. [Google Scholar] [CrossRef]
  2. Oureilidis, K.O.; Bakirtzis, E.A.; Demoulias, C.S. Frequency-based control of islanded microgrid with renewable energy sources and energy storage. J. Mod. Power Syst. Clean Energy 2016, 4, 1–9. [Google Scholar] [CrossRef]
  3. Resende, F.; Gil, N.; Lopes, J. Service restoration on distribution systems using Multi-MicroGrids. Eur. Trans. Electr. Power 2011, 21, 1327–1342. [Google Scholar] [CrossRef]
  4. Shahnia, F.; Chandrasena, R.P.; Rajakaruna, S.; Ghosh, A. Primary control level of parallel distributed energy resources converters in system of multiple interconnected autonomous microgrids within self-healing networks. IET Gener. Transm. Distrib. 2014, 8, 203–222. [Google Scholar] [CrossRef]
  5. Nutkani, I.U.; Loh, P.C.; Blaabjerg, F. Distributed operation of interlinked AC microgrids with dynamic active and reactive power tuning. IEEE Trans. Ind. Appl. 2013, 49, 2188–2196. [Google Scholar] [CrossRef]
  6. Kamel, R.M.; Chaouachi, A.; Nagasaka, K. Analysis of transient dynamic response of two nearby micro-grids under three different control strategies. Low Carbon Econ. 2010, 1, 39. [Google Scholar] [CrossRef]
  7. Pashajavid, E.; Shahnia, F.; Ghosh, A. Interconnection of two neighboring autonomous microgrids based on small signal analysis. In Proceedings of the 2015 9th International Conference on Power Electronics and ECCE Asia (ICPE-ECCE Asia), Seoul, Korea, 1–5 June 2015; pp. 213–220.
  8. Werth, A.; Kitamura, N.; Tanaka, K. Conceptual Study for Open Energy Systems: Distributed Energy Network Using Interconnected DC Nanogrids. IEEE Trans. Smart Grid 2015, 6, 1621–1630. [Google Scholar] [CrossRef]
  9. Koyanagi, K.; Hida, Y.; Yokoyama, R.; Nagata, S.; Nakao, K.; Hirai, T. Electricity cluster-oriented network: A grid-independent and autonomous aggregation of micro-grids. In Proceedings of the 2010 International Symposium Modern Electric Power Systems (MEPS), Wroclaw, Poland, 20–22 September 2010; pp. 1–6.
  10. Che, L.; Shahidehpour, M.; Alabdulwahab, A.; Al-Turki, Y. Hierarchical Coordination of a Community Microgrid with AC and DC Microgrids. IEEE Trans. Smart Grid 2015, 6, 3042–3051. [Google Scholar] [CrossRef]
  11. Kumar, M.; Srivastava, S.C.; Singh, S.N.; Ramamoorty, M. Development of a control strategy for interconnection of islanded direct current microgrids. IET Renew. Power Gener. 2014, 9, 284–296. [Google Scholar] [CrossRef]
  12. Shafiee, Q.; Dragicevic, T.; Vasquez, J.C.; Guerrero, J.M. Hierarchical control for multiple dc-microgrids clusters. IEEE Trans. Energy Convers. 2014, 29, 922–933. [Google Scholar] [CrossRef]
  13. Moayedi, S.; Davoudi, A. Distributed Tertiary Control of DC Microgrid Clusters. IEEE Trans. Power Electron. 2016, 31, 1717–1733. [Google Scholar] [CrossRef]
  14. Zhao, H.; Qiuwei, W.U.; Wang, C.; Cheng, L.; Rasmussen, C.N. Fuzzy logic based coordinated control of battery energy storage system and dispatchable distributed generation for microgrid. J. Mod. Power Syst. Clean Energy 2015, 3, 422–428. [Google Scholar] [CrossRef]
  15. Li, Y.; Liu, N.; Zhang, J. Jointly optimization and distributed control for interconnected operation of autonomous microgrids. In Proceedings of the 2015 IEEE Innovative Smart Grid Technologies-Asia (ISGT ASIA), Pathumwan Bangkok, Thailand, 3–6 November 2015; pp. 1–5.
  16. Gregoratti, D.; Matamoros, J. Distributed Energy Trading: The Multiple-Microgrid Case. IEEE Trans. Ind. Electron. 2015, 62, 2551–2559. [Google Scholar] [CrossRef]
  17. He, M.; Giesselmann, M. Reliability-constrained self-organization and energy management towards a resilient microgrid cluster. In Proceedings of the 2015 IEEE Power & Energy Society Innovative Smart Grid Technologies Conference (ISGT), Washington, DC, USA, 17–20 February 2015; pp. 1–5.
  18. Das, K.; Kargupta, H.; Bhaduri, K. A local distributed peer-to-peer algorithm using multi-party optimization based privacy preservation for data mining primitive computation. In Proceedings of the 2009 IEEE Ninth International Conference on Peer-to-Peer Computing, P2P’09, Seattle, WA, USA, 9–11 September 2009; pp. 212–221.
  19. Liu, N.; Chen, J.; Zhu, L.; Zhang, J.; He, Y. A Key Management Scheme for Secure Communications of Advanced Metering Infrastructure in Smart Grid. IEEE Trans. Ind. Electron. 2013, 60, 4746–4756. [Google Scholar] [CrossRef]
  20. Rottondi, C.; Barbato, A.; Chen, L.; Verticale, G. Enabling Privacy in a Distributed Game-Theoretical Scheduling System for Domestic Appliances. IEEE Trans. Smart Grid 2016. [Google Scholar] [CrossRef] [Green Version]
  21. Chen, Z.; Wu, L. Residential Appliance DR Energy Management with Electric Privacy Protection by Online Stochastic Optimization. IEEE Trans. Smart Grid 2013, 4, 1861–1869. [Google Scholar] [CrossRef]
  22. Yang, L.; Chen, X.; Zhang, J.; Poor, H.V. Cost-Effective and Privacy-Preserving Energy Management for Smart Meters. IEEE Trans. Smart Grid 2015, 6, 486–495. [Google Scholar] [CrossRef]
  23. Acs, G.; Castelluccia, C. I Have a DREAM! (Differentially Private Smart Metering); Springer: Berlin/Heidelberg, Germany, 2011; pp. 118–132. [Google Scholar]
  24. Marmol, F.G.; Sorge, C.; Ugus, O.; Perez, G.M. Do not snoop my habits: Preserving privacy in the smart grid. IEEE Commun. Mag. 2012, 50, 166–172. [Google Scholar] [CrossRef]
  25. Badra, M.; Zeadally, S. Design and Performance Analysis of a Virtual Ring Architecture for Smart Grid Privacy. IEEE Trans. Inf. Forensics Secur. 2014, 9, 321–329. [Google Scholar] [CrossRef]
  26. Gan, L.K.; Shek, J.; Mueller, M. Hybrid wind-photovoltaic-diesel-battery system sizing tool development using empirical approach, life-cycle cost and performance analysis: A case study in Scotland. Energy Convers. Manag. 2015, 106, 479–494. [Google Scholar] [CrossRef]
  27. Zhou, N.; Liu, N.; Zhang, J.; Lei, J. Multi-Objective Optimal Sizing for Battery Storage of PV-Based Microgrid with Demand Response. Energies 2016, 9, 591. [Google Scholar] [CrossRef]
  28. Ogunjuyigbe, A.S.O.; Ayodele, T.R.; Akinola, O.A. Optimal allocation and sizing of PV/Wind/Split-diesel/ Battery hybrid energy system for minimizing life cycle cost, carbon emission and dump energy of remote residential building. Appl. Energy 2016, 171, 153–171. [Google Scholar] [CrossRef]
  29. Salazar, J.; Tadeo, F.; Prada, C.D. A microgrid library in a general simulation language. IFAC Proc. Vol. 2014, 47, 3599–3604. [Google Scholar] [CrossRef]
  30. Liu, N.; Tang, Q.; Zhang, J.; Fan, W.; Liu, J. A hybrid forecasting model with parameter optimization for short-term load forecasting of micro-grids. Appl. Energy 2014, 129, 336–345. [Google Scholar] [CrossRef]
  31. Shi, W.; Xie, X.; Chu, C.C.; Gadh, R. Distributed Optimal Energy Management in Microgrids. IEEE Trans. Smart Grid 2015, 6, 1137–1146. [Google Scholar] [CrossRef]
  32. Nguyen, T.A.; Crow, M. Stochastic Optimization of Renewable-Based Microgrid Operation Incorporating Battery Operating Cost. IEEE Trans. Power Syst. 2015, 31, 2289–2296. [Google Scholar] [CrossRef]
  33. Zhao, B.; Zhang, X.; Chen, J.; Wang, C.; Guo, L. Operation optimization of standalone microgrids considering lifetime characteristics of battery energy storage system. IEEE Trans. Sustain. Energy 2013, 4, 934–943. [Google Scholar] [CrossRef]
  34. Saad, W.; Han, Z.; Poor, H.V. Coalitional Game Theory for Cooperative Micro-Grid Distribution Networks. In Proceedings of the 2011 IEEE International Conference on Communications Workshops (ICC), Kyoto, Japan, 5–9 June 2011; pp. 1–5.
  35. Boyd, S.; Parikh, N.; Chu, E.; Peleato, B.; Eckstein, J. Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers. Found. Trends Mach. Learn. 2011, 3, 1–122. [Google Scholar] [CrossRef]
  36. Sulc, P.; Backhaus, S.; Chertkov, M. Optimal distributed control of reactive power via the alternating direction method of multipliers. IEEE Trans. Energy Convers. 2014, 29, 968–977. [Google Scholar] [CrossRef]
  37. Erseghe, T. Distributed optimal power flow using ADMM. IEEE Trans. Power Syst. 2014, 29, 2370–2380. [Google Scholar] [CrossRef]
  38. Paillier, P. Public-key cryptosystems based on composite degree residuosity classes. In Proceedings of the International Conference on the Theory and Applications of Cryptographic Techniques, Prague, Czech Republic, 2–6 May 1999; Springer: Heidelberg, Germany, 1999; pp. 223–238. [Google Scholar]
  39. Lu, R.; Liang, X.; Li, X.; Lin, X.; Shen, X.S. Eppa: An efficient and privacy-preserving aggregation scheme for secure smart grid communications. IEEE Trans. Parallel Distrib. Syst. 2012, 23, 1621–1631. [Google Scholar]
  40. Jung, T.; Li, X.Y. Collusion-Tolerable Privacy-Preserving Sum and Product Calculation without Secure Channel. IEEE Trans. Dependable Secur. Comput. 2015, 12, 45–57. [Google Scholar] [CrossRef]
Figure 1. Centralized optimal scheduling of the interconnected microgrids (IMG).
Figure 1. Centralized optimal scheduling of the interconnected microgrids (IMG).
Energies 09 01031 g001
Figure 2. Distributed optimal scheduling of the IMG.
Figure 2. Distributed optimal scheduling of the IMG.
Energies 09 01031 g002
Figure 3. Information transmission process of the privacy-preserving strategy.
Figure 3. Information transmission process of the privacy-preserving strategy.
Energies 09 01031 g003
Figure 4. Structure of the IMG.
Figure 4. Structure of the IMG.
Energies 09 01031 g004
Figure 5. Net power curves of the MGs.
Figure 5. Net power curves of the MGs.
Energies 09 01031 g005
Figure 6. Convergence of the distributed optimization algorithm.
Figure 6. Convergence of the distributed optimization algorithm.
Energies 09 01031 g006
Figure 7. Evolution of the internal variables for the 26th time slot.
Figure 7. Evolution of the internal variables for the 26th time slot.
Energies 09 01031 g007
Figure 8. Scheduling results of the controllable DERs using the distributed optimization.
Figure 8. Scheduling results of the controllable DERs using the distributed optimization.
Energies 09 01031 g008
Figure 9. Decided exchanging power among MGs.
Figure 9. Decided exchanging power among MGs.
Energies 09 01031 g009
Figure 10. Errors between the operation cost of centralized optimization and distributed optimization.
Figure 10. Errors between the operation cost of centralized optimization and distributed optimization.
Energies 09 01031 g010
Figure 11. Sensitivity analysis for the different capacities of BESS.
Figure 11. Sensitivity analysis for the different capacities of BESS.
Energies 09 01031 g011
Table 1. Parameters of the distributed energy resources (DERs). WT, wind turbine.
Table 1. Parameters of the distributed energy resources (DERs). WT, wind turbine.
DERRated Capacity (kW)DERRated Capacity (kW)
MG1.WT1500MG3.PV3800
MG1.PV1500MG1.DG1500
MG2.PV2800MG2.DG2800
MG3.WT3500MG3.DG31000
Table 2. Parameters of the battery energy storage systems (BESSs).
Table 2. Parameters of the battery energy storage systems (BESSs).
BESSRated Power (kW)Rated Energy (kWh)I (CYN)
MG1.BESS1250800800,000
MG2.BESS235010001,000,000
MG3.BESS340012001,200,000
Table 3. Comparison of the operation cost by using the centralized and distributed optimization.
Table 3. Comparison of the operation cost by using the centralized and distributed optimization.
MGsCentralized Optimization (CNY)Distributed Optimization (CNY)
MG13211.633213.14
MG26754.126747.35
MG35515.765520.88
Total15,481.5215,481.37
Table 4. Cost and discarded renewable energy resources (RES) energy in isolated and interconnected operation.
Table 4. Cost and discarded renewable energy resources (RES) energy in isolated and interconnected operation.
MGsCost of Isolated Operation (CNY)Discarded RES Energy of Isolated Operation (kWh)Cost of Interconnected Operation (CNY)Discarded RES Energy of Interconnected Operation (kWh)
MG1190.211128.132815.830
MG218,616.3005856.920
MG34622.87846.465307.430
Total23,429.381974.5915,481.370
Table 5. Operation cost of the MGs considering the forecasting errors.
Table 5. Operation cost of the MGs considering the forecasting errors.
MGsDistributed Scheduling (CNY)Considering Forecasting Error (CNY)Cost Errors (%)
MG13213.143040.24−5.38
MG26747.356758.25+0.16
MG35520.885396.10−2.2
Total15,481.3715,194.60−1.85
Table 6. Execution time of the experiment.
Table 6. Execution time of the experiment.
Scheduling (s)Encryption (s)Transmission (s)Decryption (s)
Average0.02690.00421.8870.0039
Minimize0.01460.00321.330.0036
Maximize0.06170.00983.150.0052

Share and Cite

MDPI and ACS Style

Liu, N.; Wang, C.; Cheng, M.; Wang, J. A Privacy-Preserving Distributed Optimal Scheduling for Interconnected Microgrids. Energies 2016, 9, 1031. https://doi.org/10.3390/en9121031

AMA Style

Liu N, Wang C, Cheng M, Wang J. A Privacy-Preserving Distributed Optimal Scheduling for Interconnected Microgrids. Energies. 2016; 9(12):1031. https://doi.org/10.3390/en9121031

Chicago/Turabian Style

Liu, Nian, Cheng Wang, Minyang Cheng, and Jie Wang. 2016. "A Privacy-Preserving Distributed Optimal Scheduling for Interconnected Microgrids" Energies 9, no. 12: 1031. https://doi.org/10.3390/en9121031

APA Style

Liu, N., Wang, C., Cheng, M., & Wang, J. (2016). A Privacy-Preserving Distributed Optimal Scheduling for Interconnected Microgrids. Energies, 9(12), 1031. https://doi.org/10.3390/en9121031

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