Mathematical Problems in Engineering

Volume 2017 (2017), Article ID 5812394, 9 pages

https://doi.org/10.1155/2017/5812394

## Short-Term Photovoltaic Power Generation Forecasting Based on Multivariable Grey Theory Model with Parameter Optimization

School of Computer and Information Engineering, Hubei University, Wuhan, Hubei 430062, China

Correspondence should be addressed to Wenyang Cao; moc.qq@567618875

Received 29 July 2016; Revised 6 December 2016; Accepted 26 December 2016; Published 19 January 2017

Academic Editor: Wanan Sheng

Copyright © 2017 Zhifeng Zhong 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

Owing to the environment, temperature, and so forth, photovoltaic power generation volume is always fluctuating and subsequently impacts power grid planning and operation seriously. Therefore, it is of great importance to make accurate prediction of the power generation of photovoltaic (PV) system in advance. In order to improve the prediction accuracy, in this paper, a novel particle swarm optimization algorithm based multivariable grey theory model is proposed for short-term photovoltaic power generation volume forecasting. It is highlighted that, by integrating particle swarm optimization algorithm, the prediction accuracy of grey theory model is expected to be highly improved. In addition, large amounts of real data from two separate power stations in China are being employed for model verification. The experimental results indicate that, compared with the conventional grey model, the mean relative error in the proposed model has been reduced from 7.14% to 3.53%. The real practice demonstrates that the proposed optimization model outperforms the conventional grey model from both theoretical and practical perspectives.

#### 1. Introduction

With population growth, economic development, and nuclear confidence crisis, many countries are changing the energy structure and promoting the rapid development of renewable energy. Among them, the solar energy is being largely involved due to its highest sustainable development capability. However, photovoltaic power generation suffers from apparent intermittence and volatility resulting from illumination intensity, temperature, and so forth, which would cause alteration of both steady and transient characteristics of the power system when merged with current power grid. In this sense, the grid system planning, operation, and economic analysis will be largely impacted. As such, it is of great help to make accurate power output prediction of photovoltaic power station with the aim of coordination of conventional power and photovoltaic power, timely scheduling adjustment and proper power grid operation mode arrangement in advance. With the aid of prediction, on the one hand, the adverse effects of merging with photovoltaic power will be reduced, and the operational security and reliability of power system will be increased. On the other hand, by involving solar energy resource, the spinning reserve capacity and running cost of power system will be reduced as well as greater economic and social benefits being achieved.

Currently, a number of models are being applied for photovoltaic power generation prediction. In terms of prediction theory and methodology, they can be classified into three categories: neural network based model (NN) [1–4], time series model [5–7], and time trend extrapolation model [8]. Among these models, NN benefits from high prediction accuracy; however, it suffers from complex modeling together with high requirements of data samples, complicated training of models, and high cost. Time series model has less computational load; however, its prediction accuracy is not acceptable [5]. Markov model poses high requirement for classification scope, which is largely experience dependent. Generally speaking, the wider the scope, the simpler the model and, hence, the less accurate the prediction result, and vice versa [8].

Grey model (GM) is being widely used in data prediction due to lots of advantages. The main ones are that only few samples are needed and consideration of their distribution and variation trend is not necessary. In addition, the model benefits from low computational complexity, high accuracy of short-term prediction, easy checking, and so forth [9–15]. He and Li [10] proposed an enhanced residual error modifying model for power generation prediction for 5.6 kW photovoltaic system. However, the factor that the variation of daily power generation greatly depends on the system itself, external environment, and so forth was not considered in this model. Towards this issue, Zhong et al. [15] derived a model and obtained a good prediction result. It is reasonable to apply grey theory into photovoltaic prediction, in terms of the feature of grey theory and photovoltaic system. However, the existing model is not adapted to the photovoltaic system in this paper because of the difference of limited condition and photovoltaic data we have obtained. Therefore, how to improve the grey theory and make it applicable to the actual situation in this paper and perform better in prediction is the focus of this paper. In the further study, it is found that generation of background value in grey theory is of great importance in data prediction. Following up, Zhuang [16] verified that the prediction failed in the case of using model with in background value formula if power generation fluctuated dramatically. Lin et al. [14] regenerated a novel background value formula and proposed an optimized multivariable grey model based on the formula. It is demonstrated that the proposed model performed well in road displacement prediction.

In this paper, an integrated particle swarm optimization and multivariable grey theory model is applied for ground value formula. It is expected that, by using this method, the prediction accuracy will be largely improved. Further, to verify the feasibility of the proposed model, large amounts of real data from two separate power stations are employed for verification. The experimental results demonstrate the full functionality of the proposed mathematical tool.

This paper is structured as follows. Section 2 discusses the fundamental principle of multivariable grey theory model as well as key issues of current model for forecasting. Section 3 describes the general forecasting procedure of using the proposed optimization model. The real data from power station and their prediction results by using both the proposed model and the old model are discussed in Section 4. In Section 5, conclusions are drawn.

#### 2. The Multivariable Grey Theory Model with Parameter Optimization

##### 2.1. Overview of Multivariable Grey Theory Model

The fundamental principle of multivariable grey theory is described as follows.

###### 2.1.1. Accumulative Sequence Generation

Suppose , for , is the sample sequence, where is predicting sequence and are correlative sequences. After the Accumulated Generation Operation (AGO), the sequence is expressed aswhere

###### 2.1.2. Equation Establishment

The Albino equation of model is shown as follows:where is a parameter.

The differential equation of is given bywhere is the background value, which is generated from :

###### 2.1.3. Parameter Derivation

By using the Least Square (LS) algorithm, the parameters in (4) can be obtained as where

###### 2.1.4. Prediction Formula Generation

The approximate time response of model is given by

The original data sequence can be retrieved by Inverse Accumulated Generation Operation (IAGO) when :

##### 2.2. Problems of Multivariable Grey Theory Model

Assuming the existing data sequence is denoted as , , then, according to multivariable grey theory mentioned above, the background value of can be obtained by average operation between neighboring data:where is the last information and is the latest information.

The coefficients of are the weights on old and new information. Note that the sum of the two coefficients will always be zero; the larger the value of one side is, the more important it is and the smaller the value of the other side is, and vice versa. It can be seen from (10) that is generated under the condition of equal weight between old and new information. Generally speaking, with the lack of the old and new information’s reliability, it is more likely to choose equal weights. However, accurate prediction can hardly be expected in this case [8]. From , the parameter in background value can be derived:

It can be seen that the limit of is 0.5 when is approaching zero, while it deviates by 0.5 when the absolute value of is large. Bringing (8) to (9), we can obtainwhere , ; when bringing (12) and (5) into (4), we can get

According to (8),Then, bringing (14) into (13), the following equation is made:It is well understood that both and are constants when the sample sequence is determined. In addition, we can tell that parameter in (4) is parameter dependent in (3).

By integrating (1) in the interval of , we can obtainSince the term at the right side of (16) can be served as the grey constants, (16) can be rewritten asCombining (18) and (2), we can obtain

So, the real background value equals the integration of in the interval of, derived from (19). And the background value in the simple model was generated from the neighboring average. We should not just letbe equal to 0.5, which is the huge limitation in the simple multivariable grey theory model. At the same time, since is an ascending sequence, the value of is always between 0 and 1.

It is observed that when the time gap is small and data sequence keeps flat, conventional multivariable grey theory model is feasible to some extent. However, when the data changes fast and dramatically, this model may cause a large error [12]. In this paper, a kind of particle swarm algorithm was used for background value optimization. It is expected that, via this way, the better recovered results can be obtained and less error can be made.

##### 2.3. The Particle Swarm Optimization

Particle swarm optimization (PSO) was first proposed by Kennedy and Eberhart in 1995 [17]. It is a kind of swarm intelligence algorithm and is being widely used in various disciplines as well as engineering area due to its simple structure, fast convergence, and robustness [18, 19]. The algorithm is described as follows.

Suppose there are particles in -dimensional space; first, randomly select the initial velocity and position of each particle. Then, update these two parameters by iteration that involves local extreme value and global extreme value . Finally, the velocity and position of the th time can be given bywhere is the inertia weight factor, and are the training coefficients, and and are random numbers between 0 and 1.

#### 3. The Optimization Model Based Algorithm Design

There are a number of factors that impact daily photovoltaic power generation volume. They are usually classified into two categories: systematic factors and external factors. The former include efficiency of transformation between solar energy and battery and inclination of battery panel and power capacity, while the latter consist of air temperature, solar radiation intensity, weather, evaporation, and so forth. The impact of systematic factors has been considered in historic power generation already and thus can be omitted in the following prediction in this paper, while the external factors such as solar radiation intensity and air temperature are the main concerns for daily power generation volume [20, 21], which is of high priority in our modeling. Owing to the existence of historical data in our system database, in this paper, two parameters, that is, solar radiation intensity and air temperature, are being used as correlative inputs into portfolio model for prediction of daily photovoltaic power generation.

According to correlation analysis as well as fundamentals of multivariable grey model and PSO algorithm, the procedure of the proposed algorithm is described as follows:(a)Extract all the historic power generation volume, air temperature (average for day time), and illumination intensity from archive (the data of power generation volume comes from the inverter; solar radiation intensity and air temperature come from environmental monitor).(b)Generate sample matrix and accumulative sequence subsequently.(c)Set parameters of PSO, including training factors, weight, lower and upper bounds of position and velocity, number of initial particles, and maximum number of iterations.(d)Obtain the Fitness Function of PSO, which is given by deviation between fitted values and real values of the sample sequence:(e)Initialize the position and velocity of each particle as well as its local extreme and global extreme and then calculate the level of fitness of each particle.(f)Update local extreme and global extreme of each particle according to its level of fitness.(g)Iterate and update particle’s position and velocity based on (12) and (13) every loop.(h)Iteration continues until the number of iterations exceeds the max number. Then, the particle position can be obtained in terms of global extreme.(i)The prediction is achieved by inputting to multivariable grey model.

The detailed flowchart of the proposed model is shown in Figure 1.