Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 135673, 8 pages

http://dx.doi.org/10.1155/2015/135673

## Stochastic Optimal Dispatch of Virtual Power Plant considering Correlation of Distributed Generations

^{1}School of Electrical Engineering, Southeast University, Nanjing 210096, China^{2}Department of Mathematics and Research Center for Complex Systems and Network Sciences, Southeast University, Nanjing 210096, China^{3}Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia

Received 22 October 2014; Revised 3 March 2015; Accepted 5 March 2015

Academic Editor: Hector Puebla

Copyright © 2015 Jie Yu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Virtual power plant (VPP) is an aggregation of multiple distributed generations, energy storage, and controllable loads. Affected by natural conditions, the uncontrollable distributed generations within VPP, such as wind and photovoltaic generations, are extremely random and relative. Considering the randomness and its correlation of uncontrollable distributed generations, this paper constructs the chance constraints stochastic optimal dispatch of VPP including stochastic variables and its random correlation. The probability distributions of independent wind and photovoltaic generations are described by empirical distribution functions, and their joint probability density model is established by Frank-copula function. And then, sample average approximation (SAA) is applied to convert the chance constrained stochastic optimization model into a deterministic optimization model. Simulation cases are calculated based on the AIMMS. Simulation results of this paper mathematic model are compared with the results of deterministic optimization model without stochastic variables and stochastic optimization considering stochastic variables but not random correlation. Furthermore, this paper analyzes how SAA sampling frequency and the confidence level influence the results of stochastic optimization. The numerical example results show the effectiveness of the stochastic optimal dispatch of VPP considering the randomness and its correlations of distributed generations.

#### 1. Introduction

The continuing growth of energy consumption and environmental pollution spurs the application of distributed generations. But with more and more distributed generations connecting into the power grid, it became a serious problem to effectively manage a number of diffuse and different distributed generations, which maybe impede the green energy application. Virtual power plant (VPP) provides a practical and feasible solution to operational management of massive distributed generations, so as to make full use of distributed generations and increase the reliability of power supply. Conceptually, VPP is the collection of multiple decentralized distributed generations, energy storage, and some loads, with the characteristic of controllability like conventional power plants. Currently, the VPP is still in theoretical research and early stage of the pilot phase. Early VPP project [1] is mainly concentrated in Europe and the United States, whose research content focuses on interactive information architecture [2], operation control architecture [3], auxiliary services [4], and VPP energy management optimization scheduling.

VPP, whose essence is energy management schedule, should realize the optimal generation coordination among distributed generations, controllable loads, and storage devices within itself [5]. VPP contains a number of distributed generations with different types and different features, so the VPP energy management becomes more complicated than traditional power plants. From the view of optimization, the VPP optimal dispatch mainly focuses on the VPP costs minimum or profits maximum. The uncontrolled distributed generations within the VPP, such as wind and photovoltaic, show strong random features affected by natural conditions. The uncertainties need to be given full consideration in the process of VPP optimization scheduling. Vasirani et al. [6] consider electric vehicles to balance the uncertainty of wind generation in the VPP dispatch, which participates in short-term power market. The VPP studied by Shropshire et al. [7] contains a coastal wind farm and a small medium scaled nuclear power plant, so the nuclear power would effectively stabilize the volatility of wind generation. Houwing et al. [8] use deterministic power generation, such as microgas turbine, to burden the volatility of wind generation. From the perspective of demand side, Mohammadi et al. [9] apply flexible loads to counteract the effects of wind generation uncertainties in a VPP. Peikherfeh et al. [10], thinking of the uncertainty caused by distributed generation prediction error, deal with it by VPP spare capacity schedule. Other researchers construct VPP stochastic optimal dispatch, which previously focuses on VPP expectation [11], chance constrained [12], and two-stage stochastic optimization model [13–15]. Pandžić et al. [13], considering the uncertainty of distributed generations and electricity prices in the power market, build VPP two-stage stochastic optimization model with the goal of VPP profits maximum. A little similar VPP stochastic optimization model is built by Handschin et al. [14], aiming at VPP, costs minimum with the constraint conditions that not only contains electrical balance but also considers thermal equilibrium. Moghaddam et al. [15] provide a day-ahead optimal schedule for VPP systems containing wind generations and hydropower units.

Within the same region, distributed generations, consisting of VPP, theoretically have output power correlations due to the geographical proximity. But seldom research results considering correlations of distributed generations are published till now. This paper analyzes the problem of VPP optimal dispatch based on the correlations between wind and photovoltaic generations. The probability distributions of independent wind and photovoltaic generations are described by empirical distribution functions, and their joint probability density models are established with Frank-copula function according to their joint histogram and Squared Euclidean Distance. Containing random variables and their correlations, VPP stochastic optimization model is built. Based on sample average approximation (SAA), VPP chance constrained stochastic optimization model is converted into a deterministic optimization model. Numerical cases show effective results.

#### 2. The Randomness and Correlations in VPP

The wind and photovoltaic generations in VPP have strong randomness and correlations affected by natural conditions. The probability distributions of independent wind and photovoltaic generations could be described by empirical distribution functions, and copula function could be applied to quantitatively define the corrections of them. The empirical distribution functions are shown as follows:

Copula theory was proposed by Sklar in 1959, who pointed that any -dimensional joint cumulative distribution function can be divided into marginal cumulative distribution and a copula function. The marginal cumulative distribution describes the distribution of the variables, and the copula function describes the correlations of the variables, listed as follows:where is the joint probability distribution function of . are the probability distribution function of single variable . denotes the copula function.

Normality-copula, -copula, and Archimedes-copula are the most frequently used type of copula function.(1)Normality-copula function is as follows: where is the symmetric positive definite matrix whose diagonal elements are 1. is the standard normal distribution function whose correlation coefficient matrix is , whose marginal distribution is standard normal distributions. is the inverse function of standard normal distribution function.(2)-copula function is as follows: where is the symmetric positive definite matrix whose diagonal elements are all ones. is the distribution function, whose correlation coefficient matrix is and degrees of freedom are . is the inverse function of distribution function whose degrees of freedom are .(3)Archimedes-copula function is as follows: where , as a convex decreasing function, is the generator of Archimedes function. is the inverse function of . Among the different types of Archimedes-copula function, Gumbel-copula, Clayton-copula, and Frank-copula are the most frequently applied.(a)Gumbel-copula is as follows:(b)Clayton-copula is as follows:(c)Frank-copula is as follows: where , are the empirical distribution functions of random variables and indicates the connection parameters of the joint distribution.

Different copula function can describe different correlation of the random variables, so we need to select the optimal copula function to describe the correlation of specific variables. Squared Euclidean Distance is applied to select the optimal copula function. Firstly, the Squared Euclidean Distance of all kinds of theoretical copula and empirical copula function are calculated. The smallest distance theoretical copula is the optimal copula to describe the correlations of variables. The expression of Squared Euclidean Distance is listed as follows:where expresses the empirical copula of variables. denotes the theoretical copula of them.

#### 3. VPP Stochastic Optimal Dispatch Mathematic Model

The VPP contains a number of uncontrollable generations, such as wind and photovoltaic and controllable generations, such as fuel cell and diesel generators, as well as batteries as energy storage devices. The aim of VPP dispatch model is to maximize the profits of VPP, considering the power grid constraints, the operating constraints of generations, and the energy storage constraints. In view of stochastic generation feature of the uncontrollable generations (wind and photovoltaic), VPP dispatch model would be constructed as stochastic optimal mathematic model as follows.

##### 3.1. Objective Function

Consider where is the power energy price from VPP selling to the power grid. is the load energy price bought by the customers from VPP. means the plan generation that VPP submits to the power grid control center. represents interior loads in VPP. means the output generation of distributed generator in the VPP. is the maintenance costs of distributed generator in the VPP. indicates the fuel costs function of distributed generator . The expressions of fuel cells and diesels are as follows.

The costs function of fuel cell is as follows:where denotes the fuel costs of fuel cell. expresses the price of natural gas per unit volume. is the abbreviation of the calorific value of natural gas. means the output generation of fuel cell.

The costs function of the diesel is as follows:where denotes the fuel costs of diesel. , , and are the coefficients of cost function. means the output generation of the diesel.

##### 3.2. Constraints

The generation plan constraint is formulated in (13). The VPP, undertaking power generation like traditional power plants, outputs power generation according to the generation plan from the power grid control center. VPP outputs would be determined by the total generation of all distributed generator (in the VPP) subtracting the load (in the VPP):

The upper and lower limits of each controllable distributed generation are as follows:, represent the upper and lower limits, respectively.

The operating constraints of uncontrollable distributed generation are as follows:where , represent the upper limits of wind and photovoltaic generations respectively, whose output power generation meets the copula joint probability distributed denotes the probability values of inequality. is the confidence level, which is given in advance.

Power energy storage constraints are as follows:where , represent the charge and discharge power of the battery, respectively. , denote the upper and lower limits of charging power. , denote the upper and lower limits of discharging power.

#### 4. Solution Methods and Process

Sample average approximation (SAA) is an effective method based on Monte Carlo simulation to solve stochastic optimization. This method used sample average function to approximately substitute the expected value of the function and specifically use Monte Carlo to generate sample values. Firstly, the average values of the corresponding function are calculated based on the sample values. And then it was determined whether the mean function satisfies the probability, so as to convert the uncertain chance constraints of probability factors into certain constraints. The main solving steps of this method are listed as follows.

Sample groups of wind and photovoltaic outputs power from copula joint probability distributed and shown as .

Calculate the number of samples which satisfies the constraints used indicator function . When the value of the samples satisfies the constraints, is 1. Otherwise, is 0:

Calculate the probability that the samples satisfied the constraints. The probability is written as :

Convert the original uncertainty constraint into certainty constraints:

Judge whether satisfies . If it is satisfied, SAA method is applied to obtain the average outputs power of wind and photovoltaic generations. describes the expected values of wind and photovoltaic generations. Otherwise, return to to recalculate.

#### 5. Numerical Cases

##### 5.1. Initial Data

In this paper’s numerical cases, the VPP contains a wind farm, a photovoltaic station, a fuel cell, a diesel generator, and the storage device—battery. The related data information is listed in Tables 1 and 2.