Energies 2015, 8, 2268-2289; doi:10.3390/en8032268

OPEN ACCESS

energies

ISSN 1996-1073

www.mdpi.com/journal/energies

Article

Optimal Dispatch Strategy of a Virtual Power Plant Containing Battery Switch Stations in a Unified Electricity Market

Hao Bai, Shihong Miao *, Xiaohong Ran and Chang Ye

State Key Laboratory of Advanced Electromagnetic Engineering and Technology, Huazhong University of Science and Technology, Wuhan 430074, China;

E-Mails: baihao713@gmail.com (H.B.); ranxhcsust@126.com (X.R.); yech0131@163.com (Y.C.)

* Author to whom correspondence should be addressed; E-Mail: shmiao@hust.edu.cn; Tel.: +86-27-8755-6034.

Academic Editor: Enrico Sciubba

Received: 12 October 2014 /Accepted: 16March 2015 /Published: 23 March 2015

Abstract: A virtual power plant takes advantage of interactive communication and energy management systems to optimize and coordinate the dispatch of distributed generation, interruptible loads, energy storage systems and battery switch stations, so as to integrate them as an entity to exchange energy with the power market. This paper studies the optimal dispatch strategy of a virtual power plant, based on a unified electricity market combining day-ahead trading with real-time trading. The operation models of interruptible loads, energy storage systems and battery switch stations are specifically described in the paper. The virtual power plant applies an optimal dispatch strategy to earn the maximal expected profit under some fluctuating parameters, including market price, retail price and load demand. The presented model is a nonlinear mixed-integer programming with inter-temporal constraints and is solved by the fruit fly algorithm.

Keywords: virtual power plant; distributed generation; unified electricity market; energy storage system; battery switch station; interruptible load

1. Introduction

With the graveness of environmental pollution and energy shortage, distributed generation (DG) has been a widely researched topic in energy system studies. Although DG has many merits, there are still

many problems in its application [1]. Microgrids can coordinate well technical contradictions between DG and the power grid, and also have some energy management functions. Power from DGs is used locally most of the time, so the application of DG is limited by its geographic area. It not possible to effectively use DG on a large-scale and in multiple areas and achieve the benefits of scale in the electricity market. The virtual power plant (VPP) provides new thinking and methods to solve these problems.

Figure 1 shows the basic elements in VPP. With remote control and communication technology, VPP aggregates distributed energy resources (DERs) including DGs, interruptible loads (ILs), energy storage systems (ESSs) and electric vehicles (EVs), as an integrated entity [2]. It takes part in the electricity market in the form of a "single plant" and provides convenient service to the system operator, which improves the controllability and visibility of DERs. VPP forms an appropriate and flexible portfolio for DERs to maximize expected benefits [3]. VPP realizes optimal energy allocation by exchanging energy with the grid and dispatching energy among the elements inside the VPP. As a core energy management system (EMS) in VPPs, the optimal dispatch strategy (ODS) regulates and controls the power flow between DERs inside the VPP [4,5]. ODS is aimed at profit maximization and cost minimization of VPPs [6]. The operating characteristic of ODS make EMS have two-way information flow, as EMSs input real-time information from the DERs, and simultaneously output control signals to DERs. The internal operation state is that ESSs charge or discharge, ILs join or quit operation, DGs start up or shut down and EVs charge or discharge [7-9]. The external operation state of a VPP is that the ODS makes a bidding plan for the electricity market according to the electricity market price, electricity retail price [10], forecasted load and forecasted error [11]. These papers only consider the day-ahead market, and ignore the real-time market or separate the real-time market from the day-ahead market. Based on the separate electricity market, VPP will face the dilemma of less profit.

Figure 1. Basic elements in a Virtual Power Plant.

The ODS model usually covers objective function [12] and constraint conditions, including energy balance constraints [13], battery storage constraints [14], DG operation constraints [15], IL operation constraints [16], power grid constraints [17] and the cost curve function of DERs [11,12]. Cost curves of ILs and DGs can be formulated as a quadratic polynomial of the output power. Conversely, the cost curve of an ESS is defined as a linear polynomial of the output power [13]. The energy flow loss between the VPP and the electricity market is a linear polynomial of the active power [14]. The objective function means maximum benefit and minimum cost, VPP sells electricity to earn the maximum benefit [15]. The cost contains start-up and shut-down costs of DGs, fuel costs of DGs, operation losses of ESSs, expenses of electricity purchase. The start-up cost of DGs is usually a constant [16], which does not conform to the actual operation situation, as it has a mathematical function relationship with the off time.

Reference [18] considers state of charge (SOC) to the ensure the operational safety of ESSs, but the service life of ESSs affect ageing cost, so some constraints must be considered to extend the service life of ESSs, but these constraints do not appear in the previous research [15,17]. In addition, batteries should be charged in a timely manner when they are approximately fully discharged, which is ignored in the references [11,18]. According to the dispatch strategy, ILs such as washing machines, air conditioners and water heaters, can be regulated to maintain energy balance [19]. The simplified model only limits the amount of interruption load, which cannot reflect the characteristics of ILs, so the interruption interval constraints and interruption time constraints in the IL model still need further research.

Electric vehicles are usually regarded as a battery to optimize the charge and discharge strategy [20], but the model cannot provide a visual interface for system operators. A VPP needs to obtain operation information of all EVs and make charge-discharge schedules for each EV, which increases the computational complexity and reduces ODS which is a dynamic programming project, so its solving speed and precision affect the dynamic response and economic benefit of the VPP. ODS has more constraints and numerical variables, so the solving process needs a great amount of calculations. The genetic algorithm [11], accelerated particle swarm optimization [21] and Bender's decomposition algorithm [22] have unsatisfactory solution speed, so it is necessary to choose a reasonable algorithm to improve timeliness and accuracy.

To overcome the drawbacks of existing power dispatch models, the paper proposes an optimal dispatch strategy of a VPP containing battery switch stations (BSSs) in a unified electricity market. Based on the unified electricity market, the model optimizes electricity trading in the day-ahead market to maximize the profits of selling electricity and minimize costs of purchasing electricity. DER constraints, power balance constraints, security constraints and network constraints are taken into consideration in the VPP operation. Compared to previous VPP models, the presented model expands and enriches the description of the elements in the VPP in which the state-switch of DGs, cycle life of ESSs, charging constraints of ESSs, and the interruption constraints of ILs and the described BSSs are introduced to recharge EVs and its operation details are illustrated in the paper. The model is applied to a test system to verify its validity, and the simulation results solved by the fruit fly algorithm show that VPPs provide the optimal portfolio of DERs and play a reasonable dual role in trading with the unified market, and the dispatch strategy reduces the difficulty of the power grid dispatching and improves the profitability of the VPP.

This paper is organized as follows: in the next section, we present the trading model in a unified electricity market. The third section proposes the objective function and operational constraints for the optimal dispatch strategy. The next section illustrates the fruit fly algorithm. In the fifth section, the proposed model is evaluated with a VPP case, then simulation results are presented and analyzed. Finally the conclusions are drawn in the fifth section.

2. Trading Model in a Unified Electricity Market

At present, spot transactions in the electricity market mainly use two kinds of pricing method: the market clearing price (MCP) and pay-as-bid (PAB) [23]. The MCP, also known as system marginal pricing, is also adopted in this paper. The proposed power purchase trading model and algorithm have one thing in common: day-ahead trading purchases all electricity to meet the forecasted load demand of the next day, then the Independent System Operator (ISO) buys some electricity to make up the load balance between the forecasted load and the real time load in the real-time market, where the trading model is a separate independent sequential trading model. The total electricity purchase cost is the sum of a two-part purchase cost. The model separates day-ahead trading and real-time trading, and ignores the influence caused by the day-ahead market on the real-time market, so it misses the optimal power purchase deal. To get the optimal trading scheme, the day-ahead market and real-time market should be unified to make an electricity purchase scheme with minimum cost and an electricity selling scheme with maximum benefits. Therefore the new model should combine the day-ahead market with the real-time market, breaking the constraint that the VPP should purchase or sell electricity power according to the forecasted load in the day-ahead market. The trading power in the day-ahead market is optimized to achieve the greatest gains.

A large number of single loads make up the total load in a power grid, so a single load generally represents a small proportion in total load. On the basis of the central limit theorem the total load approximately follows a normal distribution. If some single loads have a larger proportion, one can use data analysis methods, such as histograms and empirical distributions, etc., to approximately describe the total load distribution, based on historical load data.

2.1. Power Purchases in a Unified Market

Some research defines the forecasted load and forecast error as a random variable and use a normal distribution to describe their uncertainty [10,24-28]. The references [29,30] use a Gauss model to describe the load, in which it is assumed that the load at each hour follows a normal distribution with known mean and variance. In addition, some research shows how to transfer the load to a normal distribution [31-34], and the daily load can be divided into several segments to use normal distributions with different mean and variance. Douglas [34] thinks that it is a reasonable assumption to consider the load demand follows normal distributions because the system load is actually an aggregation of all end users connected to the grid. If each customer is considered a random variable, then their summed loads will approximate a normal distribution according to the central limit theorem. The paper supposes the actual load demand DLt follows a normal distribution and it can be expressed as:

DL ~ N(SL,at) (1)

where StL is the forecasted load; ot is the variance of forecasted load, that reflects the fluctuation of the actual load demand.

In a unified electricity market, the electricity purchased in the day-ahead market does not need to be equal to the forecasted load, moreover the DG and IL in the VPP are taken into consideration, then the electricity purchased in the real-time market follows a normal distribution and can be formulated as:

[bRT = SL - SL - BDA - PDa ~ N , ot)

t _ St j-tsj-t Bt G '.tp'.t

Mb _ SL ~ ¿u °IL SIL ~ BDA ~ ^ UG PDG

j_1 i_1

where i is the index of DG units, running from 1 to G ; j is the index of IL blocks, running from 1 to IL; t is the index of time periods, running from 1 to T ; B'DA is the electricity purchased in the day-ahead market; B'RT is the electricity purchased in the real-time market; S'l is the total load shedding of ILs; S^ is the load shedding of IL j; PDG is the total output of DGs; PDG is the power output of DG i; vb is the forecasted load in the real-time market. is a binary variable denoting the curtailment decision for a IL; u'g' is a binary variable denoting the generation decision for a DG. According to the MCP, the total power purchase cost Cb can be expressed as:

c - c + C

^b ~ ^DA T ^RT

CDA P'dBDA

CRT - X [BRT P'rfb (BRT )] dBR

where CDA is the cost of power purchasee in the day-ahead market; CRT is the cost of power purchases in the real-time market; p^ is the unified marginal price in the day-ahead market; p' is the adjustment marginal price in the real-time market; f (BRT) is the density function of BRT .

Differences between the forecasted load and actual load are inevitable. When B'DA cannot meet the

load demand, VPP needs to purchase extra electricity in the real-time market at a price which is higher than p'd, meanwhile when B'DA is larger than the load demand, VPP will sell redundant electricity at a price which is lower than p'd. ptr can be expressed as:

t k, brt > 0

Pr =k,BRT < 0 (5)

where p^ is the purchasing price; p^ is the selling price.

According to Equations (4) and (5), CRT can be expressed as:

CRT = X j-w [BRTptrfb (BRT dBRT

= p'bSio"[BRt fb (BrtdBRT + p'sZflB'xrft (BrtdBRT

fb Bt )--t=— e 2c' (7)

yJZKö'

(Brt "p-b )2

Further CRT can be derived from Equations (6) and (7) and expressed as:

k ( "Prs)exp(-MbV2'

t=1 v2rc

+Mb ( -P*)(-mS/at) + pr^M-b}

where o(-mS/ at) is the distribution function of a normal distribution.

2.2. Power Selling in a Unified Market

The electricity sold in the real-time market follows the normal distribution and can be formulated as:

'N , O)

s' = P' + S' - S' - S'

RT ~ DG IL DA

mS =Z uGtPG, S^visL - sl - SD

where S'DA is the electricity sold in the day-ahead market; SRT is the electricity sold in the real-time market; mS is the forecasted electricity sold in the real-time market. The total electricity sold benefits Es can be expressed as:

EDA =X ptdSD

ERT =Si_OT [SRTprrfs (SRT ^ dSR

Es = EDA + ERT

(P'rb, sRT < 0 Ip'rs, sRT ^ 0

f (srt )=-72-e

(srt-m'S )2 2a?

where EDA is the benefits of electricity sold in the day-ahead market; ERT is the benefits of electricity sold in the real-time market; f (SRT) is the density function of SRT ;When the load demand exceeds its

expectation, electricity is already inadequate after selling electricity in the day-ahead market, so the VPP must purchase electricity at the price level of p^, to meet the load demands. If electricity is also adequate after selling electricity in the day-ahead market, the VPP will continue to sell electricity at the price level p^ . According to Equations (10)-(12), ERT can be expressed as:

T 1 f 12

ERT h ( - P-)eXP(-^ +mS ( - p'rb)) -'

t t + PrbMs

3. Optimal Dispatch Strategy

3.1. Objective Function

max profit = Es + El + EBSS - Cb - SUCDG

-SDC - OC - OC - C (14)

^ ^ DG DG ESS IL

where El is the benefits of selling electricity to consumers in the VPP; EBSS is the benefits of BSS; SUCDG is the start-up cost of DGs; SDCdg is the shut-down cost of DGs; OCDG is the operating cost of DGs; OCESS is the operating cost of ESSs; CIL is the outage compensation costs for ILs.

3.2. Grid Security and Power Balance Constraints

P -pj = VWIKIcos(e„ + s, -5j) (15)

Q - Qj = MvIM sin(9 ¡y + 5 - 5j ) (16)

V ■ < v. < V (17)

min i max ^ 1 / J

Sj < Sim (18)

-»max "'] ij

where P and Pj, Qt and Qj, V and Vj are the active power, reactive power, voltage and phase angle of two terminals in a bus, respectively; 9ij is the angle of complex Y-bus matrix elements; Yij is magnitude of admittance matrix elements; Vmin and Vmax are the minimum and maximum limits on the bus voltage respectively; Sij is the line capacity; Simj ax is the line maximum capacity.

BDA + BRT + SBdisch + SEdisch + PG (19)

= SDA + SRT +SBchar +SEchar + SL - SIL

where SB^ is the discharge capacity of BSS; SE'disch is the discharge capacity of ESS; SB'char is the charge capacity of BSS; SEtch(^r is the charge capacity of ESS.

3.3. DG Model

3.3.1. Operation Constraints

pmin ^ ~Dl < pmax PDG — PDG — PDG

-pl.t+i pl.t UTTpl 1-nn i nr. — JVC/in,

lDG 1 DG — DG (21)

PG - PDG Z RDN'dg (22)

[T f-1 - MDTt ]x[ I.t-1 - I.t ]> 0 (23)

[To- -MUT]x[i ut-]>o (24)

where PDGT and PDGX are the minimum and maximum DG capacity limit for active power; RDN'DG and RUPDG are the ramp-down and ramp-up limits for a DG unit; MDT and MUT are the minimum down time and minimum up time limit for a DG unit; Toff and Ton are the number of hours for which a DG unit has been on/off; 1),, is a binary variable denoting commitment status of DG t at time t.

3.3.2. Operation Costs

The fuel cost of a DG can be formulated as a quadratic polynomial of the output power, and besides the start-up cost and shut-down cost are also covered in operation costs:

OCdg (PD'G)2 + b^DG (25)

where a'DG and b1DG are positive coefficients of the quadratic cost function.

T G ( _T°ff \

.............^ _ 1i,t_\

--DG ............" '

t=l 1=1

SUCdg =YZsUCIg ydg exp

v Tdg y

where SUCDG is the start-up cost of a DG; yDg is a binary variable denoting the start-up decision for DG i at time t. tDg is the time constant of the start-up cost for DG i.

SDCdg =X ^SDC^G (27)

t=i i=i

where SDC'dG is the shut-down cost of a DG; is a binary variable denoting the shut-down decision for DG i at time t.

3.4. IL Model

3.4.1. Operation Benefits and Costs

The benefits of selling electricity to consumers in the VPP can be expressed as:

E, =Z^t(SL -SL) (28)

where is the electricity retail price; S'IL is the total amount of load to be shed for all ILs. VPP should pay outage compensation for an IL that gets cut off.

Cd )2 + blLSiL (29)

t=1 1=1

where a'IL and b\L are positive coefficients of quadratic outage compensation function.

3.4.2. Operation Constraints

c<1 ^ rimax

SIL ^ SIL (30)

where Smx is upper limit for curtailing on interruptible load. The duration constraint for an interruption:

k+T1 'd

^Mli At < TL

._. IL

k=T (31)

t = 1,2...T- TlILd

where M'L is a binary variable denoting the interrupted status for IL i at time t; T^ is the longest duration for IL i.

The time interval constraint between two interruptions:

X(1 - ML ) At > TLn (ML1 - Mil )

t = 1,2...T - T^

where TL" is the shortest interval requirement. The interruption times constraint:

X Mil-* < Ni

ZMil+* < N

i ,max iL

¡M'L - M^ < M£-* Imi^^-1 - Ml < m'lI+*

where is the time tag when IL suffers an outage; M'l_* is the time tag when the power restoration.

The total interruption durations constraint is given by:

(34) IL enjoys

Mil At < TL

where TLmx is the longest duration for each interruption.

3.5. ESSModel

3.5.1. Operation Costs

Lead-acid and lithium-ion batteries will lose part of their available capacity after charging-discharging several times. When the available capacity drops below the capacity limit, the battery must be replaced. The battery aging/wear cost represents the operation cost of the ESS. The depth of discharge (DoD) contributes to the cycle life of an EV battery [35]. In this study, due to the lack of direct operational data on DoD, the state of charge (SOC) subtracted from 100% charge (1-SOC) is employed to represent the DoD of the BEV battery. The relationship between the cycle life and the DoD is an empirical formula:

cl = p0 (Dr)Plexp(02 (l - f)) (36)

where cl is the cycle life; D is the DoD; DR is the rated DoD; Pq=320, ^=1.703, p2= -3.59 [35].

T ESS rc

OCESS =XX (37)

t=1 i=1 Cli ,t

where rc is replacement cost for one ESS; rcjcli,t is the aging/wear cost for each discharge.

3.5.2. Operation Constraints

The SOC must be limited in a reasonable range to avoid overcharge or deep charge.

SOCmm < SOCi < SOCmax (38)

where SOC,min and SOC,max are the maximum and minimum SOC for ESS i.

The ESS only has one operating state between charge state and discharge state at the same time, the operating state constraint can be expressed as:

0 < dut + ci,t < 1 (39)

The battery should be charged in a timely manner when it's approximately fully discharged SOC, < 0.35, and the next discharge is not allowed before charging is completed.

TESS ( +1 - dht )< TESmn exp

J,SQC _ 1 - SOC

(( jSQC ^ ^

VV tsqc j )

1 - SQC™n (41)

SOCmin < SOCt < 0.35 (42)

where TE'ss"1 is the minimum time interval between the next charge and this full discharge; TSQC is DoD after the normalization; tsqc is a time constant for SOC.

3.6. BSSModel

EV is charged by charging posts in which slow charging only applies to cars parked for a longer periods of time, while quick charging has a negative effect on the lifespan of the battery. The two charging modes have a high requirement for charging time and charging spots. Battery swap mode removes battery packs that have been used in an EV [36], and installs full battery packs that have been charged in a battery switch station (BSS), so the process can be finished in a few minutes. The replacement battery packs are then charged slowly at low power grid load to ensure the battery lifespan. The battery packs in BSSs can restrain the power fluctuations caused by intermittent DGs and provide reserve capacity for the power grid [37].

As shown in Figure 2, a BSS consists of a charger controller and battery packs. The charger controller charges and discharges the battery packs under the unified control and management of the VPP. The battery packs, including empty packs and full packs, are controlled by a dispatching center. The dispatching center collects empty packs from service points, and transports them to a warehouse for recharging. Meanwhile the center dispatches full packs to service points, then rents them to EV consumers. When the power price is high, BSS may sell power to the grid by discharging, and it can charge the batteries at low price.

Figure 2. Schematic diagram of a battery switch station. 3.6.1. Operation Benefits

BSS makes a profit mainly in two ways: battery rental benefits f and selling electricity benefits f2:

EBSS = Z(f + f2 )

(1) Service points charge a battery rental fee from EV consumers, after deducting the battery depreciation cost caused by charge, so the BSS's net benefits can be expressed as:

/1 = (a - ß )SNBSs - Z

ß= rclcli,, SNi't = P IE0

v BSS 1 EV I B.

sn'bSS * N'BSS

snB'Ss f y E0 ^

— A j^BSS

where a is the battery rental fee for one battery pack; p is the battery depreciation cost for one discharge; rc is the replacement cost for one battery pack; cli,t is the cycle life of battery pack i at time t; SN'bss is the demand of battery packs at the service points in unit time; Xj is the mean power purchase price during the charging of battery pack i; nc is the charging efficiency of the battery pack; E°BSS is the rated capacity of the battery pack; PEV is the battery demand of EVs; is the number of full charged battery packs in charger i at time t.

(2) BSS earn profits by selling power to the grid at a high price and purchasing power from the grid at a lower price:

DNBss f y \

/ = £ I W*EBSs.off -XDN%S (46)

j=i V nc )

\DNgs + SN^ <

I 0 < Ei,l < Ei,l (47)

I" — BSS.off BSS.on

where DN'b'ss is the number of charge-discharge cycles for a battery pack selling power to the grid in unit time; X is the aging cost coefficient; XDN'/SS is the total depreciation cost of a battery selling power to the grid in unit time; Xc.j and Xd.j are the mean power purchase price and mean power selling price for a battery pack i; nd is the discharging efficiency of the battery pack; E'BSS off is the power sold to the grid by discharging the battery pack i; E'BtSS on is the SOC of battery pack i at time t.

3.6.2. Operation Constraints (1) Charger controller constraint:

0 < MBSS + + MBSS -< 1

where M'b'ss + is a binary variable denoting the charge decision for charger i at time t; M'b'ss- is a binary variable denoting the discharge decision for charger i at time t:

f pi,min ^ pi,t ^ pi,max I PBSS + < PBSS + < PBSS+

r}i,min ^ T)i,t ^ r)i,max (49)

IPBSS - < PBSS - < PBSS -

where PBbSs + and P'bSs- are the charge power and discharge power for charger i at time t; PBSS+ and PBsm+n are the maximum and minimum charge power for charger i; PBSS- and PBsSS- are the maximum and minimum discharge power for charger i.

f Api,min < pi,t - pi,t-1 < a pi,max I BSS + — rBSS + rBSS + — BSS+

Api,min ^ pi,t pi,t-1 ^ a pi,max

BSS - — rBSS - rBSS - — BSS -

where AP^ and APS+ are the ramp-down and ramp-up limits of the charge power for charger i; APS- and APJBSm-1 are the ramp-down and ramp-up limits of the discharge power for charger i; At is unit time.

(2) State flag constraint:

I Mi ,t _ Mi ,t-1 < ML

I IV1 BSS + 1V1 BSS + —1V1 BSS+*

IM1 ,t - M1 ,t-1 < M1 ,t (51)

llvl BSS - 1V1BSS - — 1V1 BSS -

where M1Btss+* is a flag denoting the state of charger i at time t changing to charge from discharge; M'Sis a flag denoting the state of charger i at time t changing to discharge from charge:

Mi,t < BSS.on

DOT 44 — "

1 BSS+# — „CI

Ti/fi't < 1 _ BSS .on 1V1BSS-#- 1 E 0

E't (52)

J— R

where MRSS+# is a flag denoting the charger i complete charging at time t; MRRS_# is a flag denoting the charger i completing discharge at time t.

To expand the battery service life, the online battery must complete charging before its state changes to discharge from charge, and besides the online battery must complete discharge before its state changes to charge from discharge:

\mRSS+# > mRSS+* I-MRSS_# > mRSS_*

(3) SOC constraint:

E',t _ F

RSS RSS

RSS .on

1 + (pi,t_1 _ pi,t_1 = o

' ^ y1 RSS _ 1 RSS + J1-*1 u

= E',t _( N',t_1 _ M',t_M E

~ RSS \ RSS 1V1 RSS _)

0 < F',t < E0

u — ^RSS .on ^RSS

Ni ,t =Y(M

RSS \1V1.

_ M1 ,t

RSS+# lyl RSS_#

where ERRs is the power charged by charger i at time t.

4. Fruit Fly Optimization Algorithm

The fruit fly optimization algorithm (FOA) [38] is a new global optimization method based on the food finding behavior of the fruit fly. The fruit fly has a very keen sense of sight and smell [39], its food finding process is as follows: firstly, the fruit fly smells the food source by is osphresis organ, and flies toward that location [40]; when it gets close to the food location, its sensitive vision is also used for finding food and the company's flocking, and it flies toward that direction [41]. According to the food finding characteristics of fruit fly swarm, the FOA can be divided into six steps [38]:

(1) Initialize the fruit fly swarm location (x, y), and set the maximum number of generations and

population size;

(2) Give the random direction and distance for foraging using osphresis by an individual fruit fly:

xi = x + Lr < yt = y + Lr Lr = rand (_L, L)

where L is the fixed step size of a fruit fly using osphresis for foraging.

(3) Calculate the smell concentration judgment value si for each fruit:

U =V xi + y2

U = 1/ d

where di is the distance of the food location to the origin;

(4) Calculate the smell concentration smelli for each fruit by substituting si into the smell

concentration judgment function, then find the fruit fly with maximal smell concentration:

smelli = function (si) [bestsmell, bestindex] = max (smell)

(5) Record the maximal smell concentration value and x, y coordinates. Then, the fruit flies swarm towards that location with the maximal smell concentration value by using vision:

smellbest=bestsmell l xbest = x (bestindex) .ybest = y (bestindex)

(6) End the algorithm if the maximum number of generations is reached; otherwise, repeat the implementation of steps (2)-(3).

The FOA needs few adaptable parameters and is easy to implement, it has higher optimizing precision and less calculation. In addition, the algorithm has stronger global optimization ability, and can effectively prevent premature convergence.

5. Study Case

5.1. Test System

In this section, the test system used for evaluation of the proposed model is introduced. A VPP exchanges electricity with the grid in a unified trading market. The network configuration of the VPP is shown in Figure 3. The VPP contains four DGs at bus 2, 6, 7 and 11. The characteristics and constraints of the DGs such as generation limit, cost function coefficients, minimum up/down time limit and ramping capability for reserve are shown in Table 1. The four DGs do not contain intermittent DGs, such as wind generation and photovoltaic. If DGs contain intermittent DG, the VPP also dispatches other controllable DGs, such as micro turbines, diesel generators and fuel cells, based on the forecasted output of intermittent DGs, and the dispatch strategy focuses on the power generation scheme of DGs. Therefore the absence of intermittent DGs does not affects the research idea about the dispatch strategy of the VPP.

At bus 5 and 10, the interruptible loads can be curtailed up to 20 and 30 kW. The VPP should pay an outage compensation to consumers according to CL = 0.01SfL + 1.5SfL . The rated capacity of the ESS located at bus 9 is 80 KW and of the BSS located at bus 10 is 100 KW. The forecasted load SR of VPP is shown in Figure 2a and the variance of forecasted load ot is shown in Table 2. In this case, the time cycle is 24 h.

Figure 4b-d show the retail price for the end consumers of the VPP, the marginal price in the day-ahead market as well as the selling price and purchase price in the real-time market.

Figure 3. The single line diagram of VPP test system. Table 1. The characteristic parameters of DGs.

DG P . -1 min P -1 max a b MUT MDT RUP RDN SUC SDC

DG1 0 100 0.01 10.5 2 2 10 10 8.5 8.5

DG2 0 120 0.01 8.5 1.5 1.5 12.5 12.5 10.3 10.3

DG3 0 80 0.01 9.2 1 1 1 7.5 1 7.5 7.6 7.6

DG4 0 80 0.01 12.6 4 4 14 14 15 15

Table 2. The variance of forecastes DoD load.

t 1 2 3 4 5 6 7 8 9 10 11 12

S(t) 13 12 11.5 12.5 14 15.5 17 19 18.5 20 19.5 18.5

t 13 14 15 16 17 18 19 20 21 22 23 24

S(t) 18 18 19 20 21 22 22.5 22 19.5 17 16 14

Figure 4. Cont.

Figure 4. (a) Forecasted load; (b) The EV battery demand; (c) The retail price; (d) The marginal price in the day-ahead market as well as the selling price and purchase price in the real-time market.

5.2. Simulation Result and Discussion

The paper introduces FOA to solve the optimal dispatch problem in VPP. Figure 5 shows the results of the VPP optimal dispatch problem using FOA. As shown in Figure 5a, during hours 1-5, the prices in both the day-ahead and real-time market are lower than the operation costs of DG units. Therefore the VPP tends to buy electricity from the market rather than schedule the generation of DGs, then VPP sells electricity to the end consumers and charges the ESS and battery packs in the BSS. During hours 10-15 and 18-22, the prices in the day-ahead market are greater than the production costs of the DG units. According to load demand and storage needs, the VPP dispatches DGs to operate at full capacity to make a large gain. During hours 15-18, the market price is at a low point, so DGs tend to slash or suspend production. Due to the higher cost of star-up and shut-down and the lower limit of ramp-up and ramp-down, the DG1, DG2 and DG3 select to slash production and DG4 chooses shutdown to reduce spending.

Figure 5b explains the curtailed load of the IL. During hours 15-17 and 21-24, the retail price is larger than the market price, so the VPP earns more benefits by selling electricity to all consumers. During that period, IL does need to be curtailed. Although the retail price is also at a higher level in the time period between hours 1 to 4, the VPP cuts off all the IL to charge the ESS and the empty battery packs in the BSS. The main reasons for giving a priority to BSS and ESS is the potential profit of the next period. In other times, the retail price is lower than the market price, so the VPP makes more profit by flexibly curtailing the loads.

Figure 5c illustrates the SOC change of ESSs and BSSs, during hours 1-6 and 20-24, based on the less load demand lower market price, when the VPP purchases power from the market then stores electricity in ESSs, so the ESSs are charged to full capacity. During hours 14-17, the VPP does the same thing after selling electricity to its end consumers, ESS waits a reasonable time to discharge electricity. In other time periods, the VPP sells electricity to the market for more benefits by

discharging ESS. During hours 8-11 and 16-22, EVs have more demands for battery replacement. According to the dispatch strategy, the BSS purchases low-price electricity to produce full battery packs before the load peaks. Because the battery rental price is larger than the retail price, the BSS does not sell any surplus electricity to the end consumers in the VPP rather than trade electricity with the market. In other times, the BSS develops a reasonable charge-discharge plan after meeting the battery demands of EVs, so the SOC curve has more fluctuation than the ESS.

Figure 5d shows the electricity transactions of the VPP in a unified market. The positive values denote that the VPP sells power to the energy market as a producer, while the negative values denote that the VPP purchases power from the energy market as a consumer. The dispatch strategy optimizes the electricity trading in the day-ahead market, and purchases or sells electricity to meet the difference between the forecasted load and actual load. The VPP can make more benefits in a unified market than a separate market. Table 3 illustrates the profits of the VPP in a unified market and separate market. The total net benefits are 11,483.8 monetary units in the unified market and 9091.8 monetary units in a separate market, so the model proposed in this paper can help the VPP make more benefits.

Figure 5. (a) Generation scheduling of DGs; (b) Interruptible loads; (c) Charging and discharging behaviors of BSS and ESS; (d) Electricity transaction in unification market.

Table 3. The net benefits in a unified market and separate market.

Hour Revenue in Unified Market Cost in Unified Market Net Benefits in Unified Market Net Benefits in Separate Market

1 941.7 623.1 318.6 181.6

2 729.9 507.7 222.2 165.2

3 801.9 566.0 235.9 165.9

4 863.3 649.9 213.4 163.4

5 1,008.3 707.4 300.8 180.8

6 1,390.2 1,005.1 385.1 320.1

7 1,476.4 1,005.1 471.3 341.3

8 1,751.3 1,230.3 521.0 386.0

9 1,931.9 1,313.3 618.6 492.6

10 1,954.4 1,313.3 641.2 518.2

11 1,517.4 1,042.9 474.5 328.5

12 2,055.9 1,313.3 742.6 626.6

13 2,063.2 1,218.5 844.7 753.7

14 1,848.8 1,218.5 630.3 530.3

15 1,481.4 999.7 481.8 418.8

16 1,085.0 707.4 377.5 234.5

17 1,376.1 1,064.6 311.5 238.5

18 1,356.2 1,064.6 291.6 213.6

19 1,915.7 1,360.7 555.0 472.0

20 2,088.9 1,360.1 728.8 597.8

21 2,359.0 1,467.3 891.7 755.7

22 1,993.9 1,467.3 526.6 447.6

23 1,046.6 707.4 339.2 249.2

24 1,473.2 1,113.3 359.9 309.9

Total 36,510.6 25,026.8 11,483.8 9,091.8

The presented dispatch strategy is solved by using Matlab under FOA on a laptop computer with a 2.5 GHz core processor. The average time consumption of five cycles is 356 s. The obtained average execution times are acceptable for solving the dispatch strategy of the VPP. If we change the time step from 1 h to 15 min, VPP can improve the accuracy of the load prediction, and the variance of the forecasted load will be reduced at least 5.3%. The VPP can cope better with short-term price spikes and load fluctuations to enhance the economic benefits. The new net benefits are 11,647.6 monetary unit, which is less than 11,483.8, but the execution times is 1572 s. The result shows that a little time step can increase less the benefits.

6. Conclusions

The paper presents a trading model of a VPP in a unified market combining a day-ahead market with a real-time market. Based on the forecasted load, the volume of electricity trading in the day-ahead market is optimized to minimize the costs of purchasing power and maximize the benefits of power selling. Trading in a unified market can produce more benefits than a separate market. The load expectation in load forecasting is facile, so the paper adopts a statistical analysis method of historical

data to achieve a normal distribution of the forecasted load. The trading model is based on MCP and it also can be applied to PAB.

The model proposed in the paper expands and enriches the description of elements in a VPP. Considering the start-up costs of a DG, the model is redefined to consider the influence of shut-down time. The duration constraint for an interruption, total interruption times, total interruption durations and the time intervals between two interruptions are covered in the IL model. The paper takes cycle life into consideration to emulate the aging cost of ESSs, and in addition, the minimum time interval between the next charge and full discharge is limited to protect the battery. The BSS is introduced to recharge EVs and its operation benefits and constraints are illustrated in the paper.

The proposed dispatch strategy is also tested with a VPP system including 11 buses, four DGs, two ILs as well as one ESS and one BSS. The simulation result achieved by the fruit fly algorithm shows that VPP can reasonably merge the various elements and flexibly dispatch power among elements based on a global benefits goal. VPP is a consumer or supplier of the unified market, and it switches roles according to the retail price and market price for increasing the potential profit. The proposed model can determine the optimal power dispatch strategy among elements in the VPP and a trading plan with the electricity market for the whole VPP.

Acknowledgments

The project was supported by the National Natural Science Foundation of China (Grant No. 51377068) and Science and Technology Foundation of the State Grid Corporation of China for the project in 2015 "Research and Application on Storage Technology with High Penetration of Distributed Photovoltaics".

Author Contributions

The author Hao Bai designed the research, conducted the programing, developed and validated the models, wrote and edited the paper. Shihong Miao provided some ideas on the paper, Xiaohong Ran checked the results, Chang Ye checked the entire manuscript.

Conflicts of Interest

The authors declare no conflict of interest. References

1. El-Khattam, W.; Salama, M.M.A. Distributed generation technologies, definitions and benefits. Electr. Power Syst. Res. 2004, 71, 119-128.

2. Pudjianto, D.; Ramsay, C.; Strbac, G. Virtual power plant and system integration of distributed energy resources. IETRenew. Power Gener. 2007, 1, 10-16.

3. Mashhour, E.; Moghaddas-Tafreshi, S. A review on operation of micro grids and virtual power plants in the power markets. In Proceedings of the 2nd International Conference on Adaptive Science & Technology (ICAST), Accra, Ghana, 14-16 January 2009; pp. 273-277.

4. Giuntoli, M.; Poli, D. Optimized thermal and electrical scheduling of a large scale virtual power plant in the presence of energy storages. IEEE Trans. Smart Grid2013, 4, 942-955.

5. Bakari, K.; Kling, W.L. Virtual power plants: An answer to increasing distributed generation. In Proceedings of the Innovative Smart Grid Technologies Conference Europe (ISGT Europe), Gothenburg, Sweden, 1-6 October 2010.

6. Trichakis, P.; Taylor, P.; Lyons, P.; Hair, R. Transforming Low Voltage Distribution Networks into Small Scale Energy Zones. Proc. ICE Energy 2009, 162, 37-46.

7. Petersen, M.; Bendtsen, J.; Stoustrup, J. Optimal dispatch strategy for the agile virtual power plant. In Proceedings of the American Control Conference (ACC), Montreal, QC, Canada, 27-29 June 2012; pp. 288-294.

8. Vahidinasab, V. Optimal distributed energy resources planning in a competitive electricity market: Multiobjective optimization and probabilistic design. Renew. Energy 2014, 66, 354-363.

9. Kok, K. Short-term economics of virtual power plants. In Proceedings of the 20th International Conference and Exhibition on Electricity Distribution—Part 1 (CIRED), Prague, Czech Republic, 1-4 June 2009.

10. Mashhour, E.; Moghaddas-Tafreshi, S.M. Bidding strategy of virtual power plant for participating in energy and spinning reserve markets—Part II: Numerical analysis. IEEE Trans. Power Syst. 2011, 26, 957-964.

11. Robu, V.; Kota, R.; Chalkiadakis, G.; Rogers, A.; Jennings, N.R. Cooperative virtual power plant formation using scoring rules. In Proceedings of the Twenty-Sixth Conference on Artificial Intelligence, Toronto, ON, Canada, 22-26 July 2012; pp. 370-376.

12. Pandzic, H.; Kuzle, I.; Capuder, T. Virtual power plant mid-term dispatch optimization. Appl. Energy 2013, 101, 134-141.

13. Mashhour, E.; Moghaddas-Tafreshi, S.M. Bidding strategy of virtual power plant for participating in energy and spinning reserve markets—Part I: Problem formulation. IEEE Trans. Power Syst. 2011, 26, 949-956.

14. Caldon, R.; Patria, A.R.; Turri, R. Optimal control of a distribution system with a virtual power plant. Bulk Power Syst. Dyn. Control 2004, 2004, 278-284.

15. Caldon, R.; Patria, A.R.; Turri, R. Optimisation algorithm for a virtual power plant operation. In Proceedings of the 39th International Universities Power Engineering Conference (UPEC), Bristol, UK, 8 September 2004; pp. 1058-1062.

16. Zdrilic, M.; Pandzic, H.; Kuzle, I. The mixed-integer linear optimization model of virtual power plant operatio. In Proceedings of 8th International Conference on the European Energy Market (EEM), Zagreb, Croatia, 25-27 May 2011; pp. 467-471.

17. Peik-herfeh, M.; Seifi, H.; Sheikh-El-Eslami, M.K. Two-stage approach for optimal dispatch of distributed energy resources in distribution networks considering virtual power plant concept. Int. Trans. Electr. Energy Syst. 2014, 24, 43-63.

18. Pandzic, H.; Morales, J.M.; Conejo, A.J.; Kuzle, I. Offering model for a virtual power plant based on stochastic programming. Appl. Energy 2013, 105, 282-292.

19. Mashhour, E.; Moghaddas-Tafreshi, S.M. Integration of distributed energy resources into low voltage grid: A market-based multiperiod optimization model. Electr. Power Syst. Res. 2010, 80, 473-480.

20. Arslan, O.; Karasan, O.E. Cost and emission impacts of virtual power plant formation in plug-in hybrid electric vehicle penetrated networks. Energy 2013, 60, 116-124.

21. Ivanecky, J.; Hropko, D.; Kovác, M. Accelerated particle swarm optimization for controlling virtual power plant consisting of renewable energy sources. J. Energy Power Eng. 2013, 7, 1408-1414.

22. Seyyed Mahdavi, S.; Javidi, M.H. VPP decision making in power markets using benders decomposition. Int. Trans. Electr. Energy Syst. 2014, 24, 960-975.

23. Huisman, R.; Huurman, C.; Mahieu, R. Hourly electricity prices in day-ahead markets. Energy Econ. 2007, 29, 240-248.

24. Liu, J.; Wu, X.; Yu, J. Distribution network reconfiguration considering load uncertainty and dependence. Trans. China Electrotech. Soc. 2006, 21, 54-59.

25. Bo, R.; Li, F. Probabilistic LMP forecasting considering load uncertainty. IEEE Trans. Power Syst. 2009, 24, 1279-1289.

26. Da Silva, A.L.; Arienti, V. Probabilistic load flow by a multilinear simulation algorithm. IETDigit. Lab. 1990, 137, 276-282.

27. Leite da Silva, A.M.; Ribeiro, S.M.P.; Arienti, V.; Allan, R.; do Coutto Filho, M. Probabilistic load flow techniques applied to power system expansion planning. IEEE Trans. Power Syst. 1990, 5,1047-1053.

28. Pedersen, L.; Stang, J.; Ulseth, R. Load prediction method for heat and electricity demand in buildings for the purpose of planning for mixed energy distribution systems. Energy Build. 2008, 40, 1124-1134.

29. Valenzuela, J.; Mazumdar, M.; Kapoor, A. Influence of temperature and load forecast uncertainty on estimates of power generation production costs. IEEE Trans. Power Syst. 2000, 15, 668-674.

30. Breipohl, A.M.; Lee, F.N. A stochastic load model for use in operating reserve evaluation. In Proceedings of the Probabilistic Methods Applied to Electric Power Systems, London, UK, 3-5 July 1991; pp. 123-126.

31. Charytoniuk, W.; Chen, M.-S.; Kotas, P.; van Olinda, P. Demand forecasting in power distribution systems using nonparametric probability density estimation. IEEE Trans. Power Syst. 1999, 14, 1200-1206.

32. Li, W.; Billinton, R. Effect of bus load uncertainty and correlation in composite system adequacy evaluation. IEEE Trans. Power Syst. 1991, 6, 1522-1529.

33. Seppälä, A. Statistical distribution of customer load profiles. In Proceedings of International Conference on Energy Management and Power Delivery, Singapore, 21-23 November 1995; pp. 696-701.

34. Douglas, A.P.; Breipohl, A.M.; Lee, F.N.; Adapa, R. Risk due to load forecast uncertainty in short term power system planning. IEEE Trans. Power Syst. 1998, 13, 1493-1499.

35. Zhou, C.; Qian, K.; Allan, M.; Zhou, W. Modeling of the cost of EV battery wear due to V2G application in power systems. IEEE Trans. Energy Convers. 2011, 26, 1041-1050.

36. Lombardi, P.; Heuer, M.; Styczynski, Z. Battery switch station as storage system in an autonomous power system: Optimization issue. In Proceedings of the Power and Energy Society General Meeting, Minneapolis, MN, USA, 1-6 July 2010.

37. Miao, Y.; Jiang, Q.; Cao, Y. Operation strategy for battery swap station of electric vehicles based on microgrid. Autom. Electr. Power Syst. 2012, 15, 33-38. (In Chinese)

38. Pan, W.-T. A new fruit fly optimization algorithm: Taking the financial distress model as an example. Knowl. Based Syst. 2012, 26, 69-74.

39. Li, H.; Guo, S.; Zhao, H.; Su, C.; Wang, B. Annual electric load forecasting by a least squares support vector machine with a fruit fly optimization algorithm. Energies 2012, 5, 4430-4445.

40. Li, H.-Z.; Guo, S.; Li, C.-J.; Sun, J.-Q. A hybrid annual power load forecasting model based on generalized regression neural network with fruit fly optimization algorithm. Knowl. Based Syst. 2013, 37, 378-387.

41. Wang, L.; Zheng, X.-L.; Wang, S.-Y. A novel binary fruit fly optimization algorithm for solving the multidimensional knapsack problem. Knowl. Based Syst. 2013, 48, 17-23.

© 2015 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/4.0/).

Copyright of Energies (19961073) is the property of MDPI Publishing and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.