Table of Contents Author Guidelines Submit a Manuscript
International Journal of Rotating Machinery
Volume 2017, Article ID 5142825, 9 pages
Research Article

Optimal Sizing for Wind/PV/Battery System Using Fuzzy -Means Clustering with Self-Adapted Cluster Number

1School of Electrical Engineering, Wuhan University, Wuhan 430072, China
2Shandong Electric Power Research Institute, Jinan 250001, China
3State Grid Jiangsu Electric Power Maintenance Branch Company, Nanjing 210013, China

Correspondence should be addressed to Hong-Kun Chen; moc.361@zsnikhc

Received 24 April 2017; Revised 14 July 2017; Accepted 16 August 2017; Published 26 September 2017

Academic Editor: Yuanzheng Li

Copyright © 2017 Xin Liu 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.


Integrating wind generation, photovoltaic power, and battery storage to form hybrid power systems has been recognized to be promising in renewable energy development. However, considering the system complexity and uncertainty of renewable energies, such as wind and solar types, it is difficult to obtain practical solutions for these systems. In this paper, optimal sizing for a wind/PV/battery system is realized by trade-offs between technical and economic factors. Firstly, the fuzzy -means clustering algorithm was modified with self-adapted parameters to extract useful information from historical data. Furthermore, the Markov model is combined to determine the chronological system states of natural resources and load. Finally, a power balance strategy is introduced to guide the optimization process with the genetic algorithm to establish the optimal configuration with minimized cost while guaranteeing reliability and environmental factors. A case of island hybrid power system is analyzed, and the simulation results are compared with the general FCM method and chronological method to validate the effectiveness of the mentioned method.

1. Introduction

Hybrid power systems (HPS) [1, 2], especially those dependent on renewable energy generations (REGs), such as solar photovoltaic (PV), together with wind turbine generations (WTGs), have been regarded as the most promising configurations for remote areas power supply, since it is neither economical nor practical for delivering power over long distances. Although these clean energies provide significant contributions and opportunities, the unpredictable nature [3] of these resources has posed serious challenges to power systems [4, 5]. In the context of remote HPS, the greatest obstacle is to maintain power balance, because the adjustable capacity depends merely on REGs and batteries [6]. Hence, the dynamic characteristic of the wind speed and solar irradiance, together with the power management of batteries, should be investigated to obtain practical configurations for HPS.

In previous literature, various methods have been introduced for sizing optimization of HPS. The stochastic nature of the REGs has been investigated with several probabilistic and chronological methods. The autoregressive moving average (ARMA) is utilized to model the uncertainties of wind generation, photovoltaic (PV) power, and load in [3, 7, 8]. However, methods for parameters estimation of ARMA are always somewhat cumbersome. Reference [9] put forward an efficient approach for sizing optimization in a stand-alone HPS with Hybrid Big Bang-Big Crunch algorithm. Reference [10] analyzed the different results obtained by four heuristic algorithms; nevertheless the uncertainties of renewable energies have not been considered detailedly. Reference [11] suggested a method for technical and economic optimization in an isolated PV system; the solar radiation classification and power supply reliability calculation are performed hourly. However, only the cluster corresponding to the minimum solar radiation is selected, which may not be suitable in the context of HPS. The other investigations are based on gross chronological data [12, 13] where the computation time is always too unbearable. In [14], the traditional fuzzy -means (FCM) is adopted, which divides the data of wind speed, solar radiation, and load evenly. Thus the inherent characteristics of the data are handled in a somewhat arbitrary way.

This proposed methodology will be complementary to the previous studies and take a step further. First of all, time series analysis [15] was used to describe the characteristic of hourly wind and solar and load data with FCM, the function of which is to group the elements of data sets that have analogous characteristics. Considering that FCM is sensitive to the initialization number of clusters [16], a parameters self-adaptive method is introduced to optimize the initial state. Furthermore, the Markov model [17] is combined to obtain the system scenarios of HPS. Then the correlation and time dependency of data sets are maintained with the time-dependent clusters of the renewable generations and load power consumption. The optimal sizes for WTGs, PV, DG, and batteries are determined with the genetic algorithm (GA), in which a power balance strategy is designed to ensure that the capitalized and operational costs are minimized and the reliability requirements, CO2 emission, and batteries constraints are preserved at the same time.

The remaining parts of the article can be demonstrated in the following manner. The models of the components in HPS and the technique of FCM with self-adapted clustering number are introduced in Section 2. Section 3 presents the objective function and constraints for optimal sizing method in HPS. Section 4 utilizes a case of stand-alone hybrid system located in Hainan, China, to verify the advantage of the proposed methodology, where the comparison between the self-adapted FCM model and the traditional model with chronological data is analyzed. In Section 5, conclusions are summarized and the relationship between reliability and cost is discussed.

2. Models of the Components in HPS

2.1. The Components in HPS
2.1.1. WTGs Generation System

The output power of each WTG [10, 18] is obtained by (1), and the power curve of IEC 61400-12 standard is displayed in Figure 1.

Figure 1: Power curve of IEC 61400-12 standard.

The parameters and are calculated bywhere denotes the rated wind speed, and are, respectively, the cut-in wind speed and cut-off wind speed. means the rated power of WTGs.

The overall wind power can be derived with (3), and the total wind energy can be obtained with (4):where is the number of wind turbines and is the time step.

2.1.2. PV System

For each PV panel, the output power can be obtained [10] with where denotes the PV module power at time . means the rated power, means the solar irradiance, and means the referenced solar irradiance, in 1000 . means the referenced temperature on the surface of panels, which can be set to be equal to 25°C. means the temperature-coefficients of PV panels, and it can be set to be equal to . The temperature of each cell is deduced with where is the atmospheric temperature. means the solar irradiance and means the standard operation cell temperature and set to be 25°C.

The overall wind power can be obtained with (7), and the total wind energy can be calculated with (8):where the number of PV panels is given with .

2.1.3. Battery Bank and the Power Balance Strategy

To accommodate the stochastic behavior of PV and wind resources, battery banks are widely utilized for hybrid power systems. The power balance strategy is mainly based on the flexibility of batteries and diesel generations. The diesel generator in [14] is the only adjustable power generation without consideration of storage devices. Thus the power balance strategy is limited to only one pattern, namely, the diesel generators run to make up for the power shortage of renewable energies. To maximize the utilization of REG and minimization of diesel generation, a power balance strategy is illustrated in Figure 2.

Figure 2: Power balance strategy of HPS.
2.2. FCM with Self-Adaptive Cluster Number

In this section, FCM clustering is modified to identify the operation state of HPS. The calculation complexity can be significantly optimized considering the number of states will be much less than the 8760 h in the chronological methods. Traditional FCM clustering algorithm can only deal with a prescribed data set with clustering number given in advance, which is not flexible in the context of large data sets. A new validity function [19] is introduced to construct the proportion of compactness and divergence; thus the cluster number can be obtained according to the given data set.

2.2.1. FCM Clustering for HPS

The FCM clustering algorithm established by Dunn and then further improved by Bezdek has been used extensively.

The given data set is divided into clusters relating to some given criterions to optimize an objective function. The problem can be formulated as

An effective partition of a given data set should be divergent and compacted. The degree of compactness and divergence are evaluated with a clustering validity function. Hence, a new validity function is adopted inwhere is the fuzzy weighting index and greater than 1 and means the clustering number. is the subscript index of the data set , in which is each pattern. is the clustering prototype. means the statistical distance of the very data set. The central vector of the overall data is given with .

2.2.2. Clustering Procedure

Then, the FCM algorithm with self-adapted clustering number is outlined with the proposed validity function . The clustering procedure is illustrated in Figure 3.

Figure 3: Self-adapted clustering procedure.

The partition matrix is set as an initial condition; merely two local values of are required to be compared since the solution is locally minimized of the object function, which validates the effectiveness of Step in Figure 3.

In [14], the cluster number for WTGs is selected by evenly dividing the range between the cut-in and cut-off wind speed. The selection of cluster number of PV and load is the same. However, the inner uncertain nature of the wind resource may be disregarded by this means, and the accuracy of this method may be reduced significantly.

In this paper, the cluster numbers for WTGs, PV, and load are obtained via the method proposed in Section 3.1. More specifically, the wind speed , solar irradiance, and load power can be divided into , , and clusters coherently. The cluster centers , , and are the representative in this cluster, namely, the representative state of wind speed, solar radiation, and load.

2.2.3. The Markov Stochastic Process

In the analysis of a stochastic process, the Markov chain is an effective method to relate the probability of a state with the frequency of the corresponding event. The operation states are , , and , where , , and are the state indices obtained by the proposed method in Section 3.1. Take a wind farm with four Markov states as an example, shown in Figure 4. The state transfer probability and failure rate among different states are given with and , respectively.

Figure 4: Four-state Markov model.

If , , are probabilistic for states ( is the cluster number), then they should satisfy

The diagonal elements are equal to the negative-sum of each off-diagonal element at column of transition probability matrix , and the elements in other positions correspond, respectively, to .where denotes the transition number from state to , and means the number of states . Rate of departure (RD) is the modulus of the diagonal elements; the frequency () and duration for state () can be formulated as

The system states of HPS are mainly determined by the combination of WT and PV states, which also determine the state of DG and batteries.

3. Sizing Optimization

3.1. Objective Functions

In the proposed optimal sizing methodology of HPS, the main goal is to determine the amounts of each kind of DG. For a given system load, the objective can be set for the overall cost optimization.where subscript denotes the type of generations, namely, WT, PV, DG, and batteries, and contains the cost of the unit, installation, and fuel consumption of the HPS.

The fuel consumption of the DGs using fossil fuel is given withThen the cost model of a diesel generator is

For the DGs using renewable energy resources like WG and PV, the operation cost can be ignored. For the DGs using fossil fuel, the operation cost should be accumulated in the studied period. The combustion of fossil fuel will contribute to the emission of CO2 and gaseous pollutants. The ramping characteristics of diesel generation are neglected in the article since the time resolution is set to be one hour.

In (17), the fuel consumption cost can be obtained.

By means of FCM clustering, it can be reduced to where can be obtained with the proposed power balance strategy.

3.2. Constraints

On the basis of normal operation in the stand-alone HPS, in order to associate the reliability factor and environmental factors, the main constraints of the proposed methodology are as follows.

3.2.1. Power Balance in the Given Time Resolution

The foundation of the sizing optimization problem is the power balance. The power balance strategy in this paper is illustrated in Figure 2 (power balance strategy of HPS). The power balance equation iswhere means output power of each kind of DG in the given time interval , the length of which is associated with the length of the planning period and is set to be 1 h in this paper.

3.2.2. The Minimum and Maximum Scale of the DGs

where the subscripts max and min mean the maximum and minimum restraints for DG scales, respectively.

3.2.3. Battery Constraints

At time , the value of charge level for each battery bank should satisfy

The battery banks capacity () sets the maximum value of charge level () and the depth of discharge (DOD) determines the minimum value of charge level for the batteries ().

3.2.4. Reliability Index

LPSP (loss of power supply probability) is selected for its simplicity, which can be derived from the division of total time of power unbalance and the studied period with The subscript “max” signifies the maximum limits.

3.2.5. Environmental Factors

The emission of the DGs using fossil fuel can be obtained with (24), when considering it in the total period, and should be constrained by a maximum .

In this article, , , and of a 30-kW diesel generator are set to be 0.028144, 0.001728, and 0.0000017.

3.3. Using Genetic Algorithm to Get Optimal Solution

The genetic algorithm (GA) is chosen to solve the sizing problem considering its ability to obtain a globally optimal solution for optimization problems. It is inspired by the process of biological evolution, namely, crossover, mutation, and selection. The population of individual solutions is repeatedly modified with a “fitness” function, typically related to the objective function.

In this article, the optimization problem (14)~(25) is handled with GA, where the variables , , , and are linked to form the gene strings in the state variable (chromosome), and (14) is set to be the fitness function.

4. Results

4.1. Case Introduction

The data from an island in Hainan province of China are used to analyze the proposed problem. There are abundant wind and solar resources in this island.

4.2. Simulation Results and Analysis

The monthly average wind speed, solar irradiation, temperature, and load power profile consumption are shown in Figures 58.

Figure 5: Wind speed.
Figure 6: Illumination rate.
Figure 7: Temperature.
Figure 8: Load profile.

Firstly, the wind speed , illumination , and temperature are imported to the HPS model to obtain the power output of WTGs () and PV (). Then the cluster numbers for WT, PV, and load are obtained via the method proposed in Section 2.1. The cluster centers Vc, Gc, and Lc are, respectively, the representative states for wind speed, solar irradiance, and load, illustrated in Tables 13.

Table 1: Markov model of WT/35 kW.
Table 2: Markov model of PV/200 W.
Table 3: Markov model of load.

According to the simulation results from self-adaptive FCM, the optimal cluster numbers for wind WT, PV, and LD are 10, 5, and 8. Namely, each WT has 10 states, each PV has 5 states, and the overall load has 8 states. There are 80 possible states for the renewable energy generations as a whole. It can be noted here that the scenario of HPS has been significantly simplified.

The operation scenario considered here is greatly simplified to be the aggregation of clustering states. The outputs of DGs and batteries are also determined by these states, according to the power balance strategy demonstrated before.

For a new state, the probability is the multiplication of state probabilities for every individual WTG, PV, and LD. The frequency and duration time for a new scenario can be obtained likewise.

Then the results from chronology-based, traditional FCM based, and the self-adapted clustering number based GAs are compared; results are illustrated in Table 4. GA is capable of realizing global optimization but cannot guarantee it. The chronological-based method requires 8760 iterative loops, and the traditional FCM based method needs 920 (=23 4 10) iterative loops, and the proposed method needs 400 (=10 5 8) iterative loops.

Table 4: Performance comparison.

It should be noted that the wind energy is superfluous in the winter nights. The output power of PV panels is redundant in the summer daytimes.

It can be found that chronology-based method is the most time-consuming due to complicate loops. With regard to the overall cost, the proposed method still has advantages. Figure 9 shows the iteration performance for the mentioned algorithms.

Figure 9: Iteration performance comparison.

The proposed self-adapted FCM clustering model is superior to traditional method in two aspects:(1)Compared to chronological-based methods, with the reduction of investigated data set, the number of system scenarios of HPS can be significantly reduced, and the computation burden and CPU time can be significantly reduced, shown in Table 4.(2)Compared to the traditional FCM based method, the clustering numbers of data sets are inherently obtained and optimized, which increases the probability to obtain a global optimum solution. While the traditional method simply makes the partition by uniform division, the scenarios selection is somewhat arbitrary and disregards the inner stochastic characteristic of the data sets. Thus the basic FCM based method finally obtains a locally optimized solution.

In the proposed method, the benefits of the renewable energies are thought to be the reduction of CO2 and the improvement of LOSP, which has been set as constraint to the very problem. Considering the impacts of different reliability index on the investment cost, let be 30000 kg/y. The overall investment cost grows higher as the reliability request () increases, as shown in Figure 10. The impact of CO2 is similar to the LPSP index.

Figure 10: Sensitivity analysis of and cost.

Actually, the benefits of the reduction of CO2 and the improvement of LOSP are negatively related to the cost optimization procedure. We have modified the discussion on this issue in our revised manuscript.

5. Conclusion

In this paper, a novel method utilizing the self-adapted FCM clustering combined with the Markov model and GA is proposed to determine the best mix of HPS. A power balance strategy is also designed to guide the optimization process. The self-adapted FCM clustering can handle the stochastic characteristics of REGs, and the Markov model can significantly reduce the operational scenarios of REGs. The proposed method has comparable competitive overall cost, and it can be concluded that the benefits of the reduction of CO2 and the improvement of LOSP are negatively related to the cost optimization procedure.

The future work will include the following:(1)Improving the clustering model to further study the correlation among renewable resources(2)Adding the local and global control strategy to the power balance analysis process(3)Extending the proposed method to the operation plan stage of the HPS.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


  1. R.-J. Wai, S. Cheng, Y.-F. Lin, and Y.-C. Chen, “Installed capacity selection of hybrid energy generation system via improved particle-swarm-optimisation,” IET Generation, Transmission and Distribution, vol. 8, no. 4, pp. 742–752, 2014. View at Publisher · View at Google Scholar · View at Scopus
  2. J. G. McGowan and J. F. Manwell, “Hybrid wind/pv/diesel system experiences,” Renewable Energy, vol. 16, no. 1-4, pp. 928–933, 1999. View at Publisher · View at Google Scholar · View at Scopus
  3. A. Arabali, M. Ghofrani, M. Etezadi-Amoli, and M. S. Fadali, “Stochastic performance assessment and sizing for a hybrid power system of Solar/Wind/Energy Storage,” IEEE Transactions on Sustainable Energy, vol. 5, no. 2, pp. 363–371, 2014. View at Publisher · View at Google Scholar · View at Scopus
  4. R. R. Barton, “Presenting a more complete characterization of uncertainty: Can it be done,” in Proceedings of the in Proceedings of the Simulation Society research workshop (INFORMS '07 ), pp. 26–60, 2007.
  5. E. Ela and M. O'Malley, “Studying the variability and uncertainty impacts of variable generation at multiple timescales,” IEEE Transactions on Power Systems, vol. 27, no. 3, pp. 1324–1333, 2012. View at Publisher · View at Google Scholar · View at Scopus
  6. A. Keane, L. F. Ochoa, C. L. T. Borges et al., “State-of-the-art techniques and challenges ahead for distributed generation planning and optimization,” IEEE Transactions on Power Systems, vol. 28, no. 2, pp. 1493–1502, 2013. View at Publisher · View at Google Scholar · View at Scopus
  7. A. Mohamed Abd El Motaleb, S. Kazim Bekdache, and L. A. Barrios, “Optimal sizing for a hybrid power system with wind/energy storage based in stochastic environment,” Renewable and Sustainable Energy Reviews, vol. 59, pp. 1149–1158, 2016. View at Publisher · View at Google Scholar · View at Scopus
  8. E. Koutroulis, D. Kolokotsa, A. Potirakis, and K. Kalaitzakis, “Methodology for optimal sizing of stand-alone photovoltaic/wind-generator systems using genetic algorithms,” Solar Energy, vol. 80, no. 9, pp. 1072–1088, 2006. View at Publisher · View at Google Scholar · View at Scopus
  9. S. Ahmadi and S. Abdi, “Application of the Hybrid Big Bang-Big Crunch algorithm for optimal sizing of a stand-alone hybrid PV/wind/battery system,” Solar Energy, vol. 134, pp. 366–374, 2016. View at Publisher · View at Google Scholar · View at Scopus
  10. A. Maleki and F. Pourfayaz, “Optimal sizing of autonomous hybrid photovoltaic/wind/battery power system with LPSP technology by using evolutionary algorithms,” Solar Energy, vol. 115, pp. 471–483, 2015. View at Publisher · View at Google Scholar · View at Scopus
  11. K. Benmouiza, M. Tadj, and A. Cheknane, “Classification of hourly solar radiation using fuzzy c-means algorithm for optimal stand-alone PV system sizing,” International Journal of Electrical Power and Energy Systems, vol. 82, pp. 233–241, 2016. View at Publisher · View at Google Scholar · View at Scopus
  12. M. Mosadeghy, R. Yan, and T. K. Saha, “A time-dependent approach to evaluate capacity value of wind and solar PV generation,” IEEE Transactions on Sustainable Energy, vol. 7, no. 1, pp. 129–138, 2016. View at Publisher · View at Google Scholar · View at Scopus
  13. L. Xu, X. Ruan, C. Mao, B. Zhang, and Y. Luo, “An improved optimal sizing method for wind-solar-battery hybrid power system,” IEEE Transactions on Sustainable Energy, vol. 4, no. 3, pp. 774–785, 2013. View at Publisher · View at Google Scholar · View at Scopus
  14. Y.-Y. Hong and R.-C. Lian, “Optimal sizing of hybrid wind/PV/diesel generation in a stand-alone power system using markov-based genetic algorithm,” IEEE Transactions on Power Delivery, vol. 27, no. 2, pp. 640–647, 2012. View at Publisher · View at Google Scholar · View at Scopus
  15. S. N. Gopiya Naik, D. K. Khatod, and M. P. Sharma, “Analytical approach for optimal siting and sizing of distributed generation in radial distribution networks,” IET Generation, Transmission and Distribution, vol. 9, no. 3, pp. 209–220, 2015. View at Publisher · View at Google Scholar · View at Scopus
  16. X. Yejun, “Fitzwilliam College,” in Proceedings of the 10th International Conference on Computer Science & Education (ICCSE' 15), pp. 22–24, Fitzwilliam College, Cambridge University, Cambridge, UK, 2015.
  17. A. S. Dobakhshari and M. Fotuhi-Firuzabad, “A reliability model of large wind farms for power system adequacy studies,” IEEE Transactions on Energy Conversion, vol. 24, no. 3, pp. 792–801, 2009. View at Publisher · View at Google Scholar · View at Scopus
  18. L. Chen, H. Chen, Z. Shu, G. Zhang, T. Xia, and L. Ren, “Comparison of inductive and resistive SFCL to robustness improvement of a VSC-HVDC system with wind plants against DC fault,” IEEE Transactions on Applied Superconductivity, vol. 26, no. 7, pp. 1–8, 2016. View at Publisher · View at Google Scholar · View at Scopus
  19. Y. Li and F. Yu, “A new validity function for fuzzy clustering,” in Proceedings of the International Conference on Computational Intelligence and Natural Computing (CINC '09), vol. 1, pp. 462–465, IEEE, June 2009. View at Publisher · View at Google Scholar · View at Scopus