Optimized Planning of Power Source Capacity in Microgrid, Considering Combinations of Energy Storage Devices

Since renewable energy resource is universally accepted as a promising method to solve the global energy problem, optimal planning and utilization of various distributed generators (DG) and energy storage (ES) devices deserve special concern. ES devices possess various characteristics in power density, energy density, response speed (switching speed) and lifetime. Besides, as different load types have various requirements on power supply reliability according to their importance, coordinated planning with consideration of reasonable matching between power source and load can efficiently improve power supply reliability and economic efficiency via a customized power supply and compensation strategy. This paper focuses on optimization of power source capacity in microgrid and a coordinated planning strategy is proposed with integrated consideration of characteristics of DG, ES and load. An index named additional compensation ratio (ACR) for balancing economic efficiency and reliability is proposed and considered in the strategy. The objective function which aims to minimize life cycle cost (LCC) is established considering economic efficiency, reliability and environmental conservation. The proposed planning strategy and optimizing model is calculated and verified through case study of an autonomy microgrid.


Introduction
With the increasing demand on energy resources, human society is confronted with significant issues including fossil resource exhaustion, climate change and pollutant emission.In recent decades, renewable energy resources (RER) such as wind and solar energy attracted world-wide attention for their merits in energy conservation as well as in reducing pollutant emission.It is highly recognized that utilization of RER will be one of the most promising methods to solve the global energy problem and achieve sustainable development [1,2].
The distributed integration of wind and solar power helps to reduce loss from long-distance transmission and is a significant method of RER utilization, but the intermittency and fluctuation of wind and solar energy will inevitably influence power supply reliability.Besides, the imbalance between the amount of generation and the load demand may result in power system instability [3].Hence energy storage (ES) devices, including lead-acid batteries, Li-ion batteries, super capacitor (SC), flywheel storage and superconducting magnetic energy storage (SMES) can be utilized for output power smoothing.A small power system consisting of power converter-based generation, ES devices, small classic synchronous generation, and various types of loads forms a microgrid (MG) [4].While a number of research works on MG have been done on various aspects, including control strategy [5][6][7][8][9][10], energy management systems [11][12][13][14][15][16][17], marketing policies [18][19][20][21] and system design [22][23][24][25][26][27], this paper mainly focuses on the optimal sizing of Distributed Energy Resources (DERs), especially the capacity of ES devices in MG.
Power source optimization of MG is basically done by minimizing economic cost while the reliability is guaranteed [28,29].A number of existing works are done with consideration of complementarity among various power sources as well as source-load coordination.In [25], both reserve sizing problem and loss of load probability (LOLP) index are integrated into the existing storage sizing problem and a probabilistic model is established to strive for a trade-off between reserve cost, storage cost and selling profit.The relationship between storage capacity, reserve capacity and LOLP index is obtained for decision-making on balancing financial cost and reliability.The planning framework proposed in [26] takes into consideration the various characteristics of ES as well as the availability of different renewable energies.The case illustrated in that work adopts hour-level data and emphasized the influence on optimization result caused by variation of diesel generation capacities.The authors of [27] propose a method for coordinated sizing of energy storage (ES) and diesel generators in an isolated microgrid based on discrete Fourier transform.The power fluctuation is decomposed into high-frequency component and low-frequency component, then ES and diesel generators are assigned to these two parts of power correspondingly for power smoothing according to their respond characteristics.The authors of [30] proposed an optimal sizing method of a Vanadium Redox battery in MG and the cost-effective solution is obtained based on the optimal ratings for the battery in both short-term (one-day) and long-term (entire year) scenarios proposed in the paper.
ES devices in the above papers basically involve battery energy storage, while other types of ESs, as well as hybrid energy storage systems (HESS), are discussed in some other works.In [31], the authors proposed a super capacitor energy controller (SCEC) and presented a method for selection of the SCEC and filter parameters as well as precise sizing of the super capacitor for a given application.Battery lifetime is extended by diverging the high-frequency power variations to the super capacitor, and the size of the super capacitor is significantly lower than the one estimated for the traditional super capacitor voltage controller.The authors of [32] investigated the mathematical model and the topology of a HESS based on superconducting magnetic energy storage (SMES) and battery, and a novel system-level control strategy for reasonable and effective power allocation between SMES and battery is proposed.In [33], the authors presented a methodology to size flywheel energy storage for grid-level applications using an optimal control law, which performs trade-offs between minimizing grid power consumption and flywheel energy storage size.The characteristics of flywheel-such as high efficiency, high power density and high self-discharge rate-are considered in the methodology, and case study in this work indicates that the flywheel reduced peak grid power by 8%.
In these works, load is normally considered as a whole and power fluctuation is decomposed into high-frequency component and low-frequency component for power dispatch among ES devices in a HESS.However, different consumers are not equally important and have various requirements on power supply reliability.The authors of [34,35] proposed the classification of load based on importance and the corresponding requirements on power recovery time.For instance, power outage longer than 200 ms is not acceptable for extremely important load (level 1 load); hence, fast-respond energy compensation devices shall be utilized to ensure uninterrupted power supply.In this aspect, although characteristics including power/energy density, efficiency and self-discharge rate are already discussed in existing works on power source sizing of microgrid, the respond speed or switching time of various compensation devices, including ES and diesel generator, are seldom considered.
In order to obtain a reasonable planning scheme with integrated consideration of both power source and load characteristics, analysis and discussion on ES respond speed and load classification is done in this paper.The paper then focuses on the source-load matching mechanism with analysis on economic loss caused by power outage, and a capacity planning strategy considering combinations of ES devices is proposed and an ACR index is used for balancing economic efficiency and reliability.With differentiated allocation of capacity of ES devices according to load importance, economic efficiency can be improved by slightly sacrificing the supply reliability of level 2 and level 3 load, which are less important.
The rest of this paper is organized as follows.In Section 2, modeling of distributed energy resources as well as major characteristics of ES devices are introduced.In Section 3, the source-load matching mechanism and capacity planning strategy are proposed.The objective function based on life cycle cost (LCC) theory is established in Section 4 and case study on an autonomy microgrid is introduced in Section 5.

Modeling of Distributed Generators
The model of a wind turbine (WT) is derived from its power output curve, as illustrated in Figure 1.Though the curve of each type of wind turbine is slightly different from others', they can be generally represented as a piecewise function [36].
where P WT is the output power of wind turbine (WT), P r is the rated power of WT, v r is the rated wind speed of WT, v i is the cut-in wind speed of WT, v o is the cut-out wind speed of WT, a 1 ~a4 are coefficients.
The output power of a photovoltaic (PV) panel array is determined by variables including P r , R c , R r , T c , T r and k [37].
where P PV is the output power of PV panel array, P r is the rated power of a single PV unit under standard conditions, R c is the current sunlight intensity, R r is the rated sunlight intensity under standard conditions, T c is the current temperature, T r is the rated temperature under standard conditions, k is the power temperature coefficient, n PV is the number of units in the PV panel array.The rest of this paper is organized as follows.In Section 2, modeling of distributed energy resources as well as major characteristics of ES devices are introduced.In Section 3, the source-load matching mechanism and capacity planning strategy are proposed.The objective function based on life cycle cost (LCC) theory is established in Section 4 and case study on an autonomy microgrid is introduced in Section 5.

Modeling of Distributed Generators
The model of a wind turbine (WT) is derived from its power output curve, as illustrated in Figure 1.Though the curve of each type of wind turbine is slightly different from others', they can be generally represented as a piecewise function [36].
where PWT is the output power of wind turbine (WT), Pr is the rated power of WT, vr is the rated wind speed of WT, vi is the cut-in wind speed of WT, vo is the cut-out wind speed of WT, a1~a4 are coefficients.
The output power of a photovoltaic (PV) panel array is determined by variables including Pr, Rc, Rr, Tc, Tr and k [37].
where PPV is the output power of PV panel array, Pr is the rated power of a single PV unit under standard conditions, Rc is the current sunlight intensity, Rr is the rated sunlight intensity under

Modeling and Characteristic Analysis of Energy Storage Devices and Back-Up Source
ES devices and back-up source are used to compensate energy shortage caused by insufficient distributed generator (DG) generation, hence relative analysis is significant for further study on source-load coordination mechanism.Table 1 shows the characteristics of both energy-type and power-type ES devices, including energy density, power density, investment cost, approximate cycle times and response speed.Through observation, it is obvious that two types of ES devices are complementary to each other, and a hybrid storage system can provide enhanced power supply capability.Besides, although flywheel storage and SMES have outstanding performance, they are relatively expensive and are still in the experimental stages, which made them unsuitable for practical application.In this paper, the HESS consists of a lead-acid battery, Li-ion battery and super capacitor.

Energy state model of ES device
The operation status of all ES devices can be generally described via the following formula [38,39]: where E ES (t) is the remaining energy of the ES device at moment t, σ is the self-discharge rate, P ES is the charge/discharge power during time interval t, η is the charge/discharge efficiency, ∆t is the duration of time interval.The energy state of the ES device is determined by E ES at the start of a time interval and charged/discharged power during that period.For different ES devices, the self-discharge rate can vary from 0.1%-0.3%/d(batteries) to 20%-40%/d (super capacitor).

Lead-acid battery lifetime model
The lifetime of a lead-acid battery is relatively short and is significantly influenced by operation conditions compared to super capacitors and a Li-ion battery.While the lifetime of super capacitor (SC) can extend as long as 20 years and a Li-ion battery can be charged/discharged for over 5000 cycles, substitution of a lead-acid battery may be required, which makes the lifetime calculation model of batteries essential for optimization of MG source planning.
An existing lifetime calculation model [37] is applied to predict the lifetime of lead-acid batteries.The cycle time of lead-acid batteries is related to depth of discharge (DOD), as illustrated in Figure 2. Through observation, it is obvious that two types of ES devices are complementary to each other, and a hybrid storage system can provide enhanced power supply capability.Besides, although flywheel storage and SMES have outstanding performance, they are relatively expensive and are still in the experimental stages, which made them unsuitable for practical application.In this paper, the HESS consists of a lead-acid battery, Li-ion battery and super capacitor.

Energy state model of ES device
The operation status of all ES devices can be generally described via the following formula [38,39]: where EES (t) is the remaining energy of the ES device at moment t, σ is the self-discharge rate, PES (t) is the charge/discharge power during time interval t, η is the charge/discharge efficiency, Δt is the duration of time interval.The energy state of the ES device is determined by EES at the start of a time interval and charged/discharged power during that period.For different ES devices, the self-discharge rate can vary from 0.1%-0.3%/d(batteries) to 20%-40%/d (super capacitor).

Lead-acid battery lifetime model
The lifetime of a lead-acid battery is relatively short and is significantly influenced by operation conditions compared to super capacitors and a Li-ion battery.While the lifetime of super capacitor (SC) can extend as long as 20 years and a Li-ion battery can be charged/discharged for over 5000 cycles, substitution of a lead-acid battery may be required, which makes the lifetime calculation model of batteries essential for optimization of MG source planning.
An existing lifetime calculation model [37] is applied to predict the lifetime of lead-acid batteries.The cycle time of lead-acid batteries is related to depth of discharge (DOD), as illustrated in Figure 2. The calculation formula can be expressed as: where DN is the depth of discharge, N is the equivalent cycle times under DN, a1~a5 are coefficients.The calculation formula can be expressed as: where D N is the depth of discharge, N is the equivalent cycle times under D N , a 1 ~a5 are coefficients.
The cost of a single charge/discharge cycle is represented as: where C cap-B is the Capital investment of batteries, C cycle is the cost of a single charge/discharge cycle.
Related studies [40,41] indicate that the cumulative lifetime of a battery is influenced by the state of charge (SOC) under which the battery is operating.Besides, a battery may not experience a fixed and whole charge-discharge cycle.Hence, the operation cost model for both the charging and discharging processes, considering the influence caused by SOC, is given in [37].

• Charging process
An influence factor λ c is introduced to represent the influence on operation cost caused by SOC: where S c-start is the SOC at the start of the charging process, S c-end is the SOC at the end of the charging process, k c is the adjustment coefficient to fix the final value of λ c .The cost of a single charging process is then shown by the following equation: where C cl is the cost of a single charging process, λ c is the influence factor representing SOC's influence on C cl .

• Discharging process
An influence factor λ d is introduced to represent the influence on operation cost caused by SOC: where S d-start is the SOC at the start of the discharging process, S d-end is the SOC at the end of the discharging process, k d is the adjustment coefficient to fix the final value of λ d .
The cost of a single discharging process is then represented as: where C dl is the cost of a single charging process, λ d is the influence factor representing SOC's influence on C dl .

Fuel consumption model of a diesel generator
The fuel consumption model of a diesel generator is obtained through fitting the fuel consumption curve [42]: where C Diesel is the amount of consumed fuel, P G is the generated power of diesel generator, P r is the rated power of diesel generator, A and B are coefficients.

Load Character Analysis
Load in a microgrid consists of different types including: industrial load; commercial load; hospitals; schools; residential load, etc. Different types of load possess unique characteristics and various requirements on power supply reliability.Basically, load can be divided into three categories according to the importance grade.The classification rule of each load grade is regulated in [34], as illustrated in Table 2.The analysis on the power supply reliability requirement of each load grade is essential for the study on the power-source coordination mechanism.Herein, the time limit of power interruption (t interrupt ) is used to measure the reliability requirement of load.This index describes the maximum acceptable power interruption time during a power outage failure.The t interrupt varies among specific types of each load grade.For instance, t interrupt of electrical equipment utilized in the electronics manufacturing industry is illustrated in Table 3.For simplification, t interrupt of a certain type of load is selected as the minimum value among the t interrupt of each type of electrical equipment.For instance, t interrupt of the electronics manufacturing industry is selected as 200 ms.t interrupt of major types of level 1 load is illustrated in Table 4.
As illustrated above, t interrupt of major types of level 1 load ranges from several hundred milliseconds up to 1 min.According to the characteristics of ES devices analyzed in Section 2.2, fast responding ES devices, i.e., super capacitor, shall be considered for improving the power supply reliability.
As regards level 2 and level 3 load, the time limit of power interruption is not specified in existing manuals, regulations and national standards.Hence, the power supply reliability of these loads is restricted by the annual power outage limit, as well as by considering energy shortage compensation during life cycle cost (LCC) calculation.A sector customer damage function (SCDF) to calculate power outage loss is proposed in [43], and the economic loss of a typical load is illustrated in Figure 3.For level 1 load, i.e. industrial load, the economic loss will rise significantly within several minutes after the outage.Hence, the interrupt time shall be limited within the minute level to avoid massive economic loss.For level 2 and level 3 load, the curve increases approximately linearly; this part of the loss is calculated and added to the life cycle cost of the project.

Power Source Capacity Planning Strategy Considering Source-Load Coordination
According to the discussion in Section 3.1, the energy shortage of level 1 load shall be fully compensated at all costs.In order to achieve uninterrupted power supply even if power generation is not sufficient, ES devices and backup sources shall be coordinated according to their response speed, discharge durations and other characteristics.Hence, super capacitors with fast response speed are utilized for smoothing short-time fluctuation or power supply during the initial stage of a long-term energy shortage, while a battery energy storage system (BESS) will successively be utilized for subsequent supply.In extreme cases with intensively long shortage duration, backup sources such as a diesel generator will operate as a stable power source until renewable energy generation is restored.
As for level 2 load, the requirement on power supply reliability is not as strict compared to level 1 load; thus, a short-time outage at the initial stage of a shortage period is acceptable in most cases and investment on SC can then be reduced.When it comes to level 3 load, BESS capacity can also be reduced since level 3 load basically has no requirement on power supply reliability.The principle of source-load coordination is illustrated in Figure 4.  BESS in this research consists of a Li-ion battery and a lead-acid battery.Compared to a lead-acid battery, a Li-ion battery has advantages including: (a) relatively high energy density; (b) longer lifetime at high depth of discharge (DOD); (c) high power receiving capability to realize high current charging in a short time, while the charging current of a lead-acid battery is usually limited within 0.25 C to avoid capacity degradation.Nevertheless, the unit cost of a Li-ion battery can be 5~7 times higher than a lead-acid battery, which is the major factor restricting its application.A reasonable portion of a Li-ion battery in BESS, which is charged and discharged with priority, will significantly reduce the utilization frequency of a lead-acid battery and extend the overall lifetime of BESS.
Through surveying typical daily load and climate parameters including temperature, sunlight intensity and wind speed, the curves of load and generated power of wind turbine (WT) and photovoltaic (PV) can be obtained.The difference of generated power and load is then defined as ∆P, and the area enclosed by the ∆P curve and the axis is ∆E i .Whether ∆E i is a positive or negative value determines the operation status of ES devices (charging or discharging) during this time interval.t B-switch and t D-switch represents the response speed, or switching time, for BESS and diesel generator respectively.
The overall compensation strategy is arranged as follows: 1. Charging process For a certain time interval with a duration of t i , if ∆E i > 0, ES devices will be charged.The redundant energy generated by WT and PV shall firstly be used to charge the SC in order to guarantee energy reserve of fast response ES for uninterrupted supply of level 1 load.BESS will then be charged if energy still remains, i.e., ∆E i -(E SC-max -E SC-available-i ) > 0 (Figure 5a).

Discharging process
For a certain time interval t i , if ∆E i < 0, ES devices will discharge to supplement the energy shortage of load.The time intervals are classified into three categories: For t i < t B-switch , SC is utilized for supplementary and the available capacity of SC is used for supplying all level 1 load and a part of level 2 and level 3 load.
For t B-switch <t i < t D-switch , SC and BESS will both be utilized.SC will be switched in immediately and discharge until t B-switch , when BESS is ready to switch in, then BESS is used for supplying all level 1 and level 2 load as well as a part of level 3 load.
For t i > t D-switch , SC will be switched in immediately and discharge until t B-switch , when BESS is ready to switch in.Then the batteries will discharge until t D-switch , when diesel generators are ready to be switched in, then a diesel generator is used for supplying all level 1, level 2 and level 3 load (Figure 5b).A flowchart of the proposed strategy is illustrated in Figure 6.A flowchart of the proposed strategy is illustrated in Figure 6.Further explanation on discharge strategy is illustrated in Figure 7.For each time interval, the energy shortage is divided into three stages distributed along the time axis in Figure 7a.Each stage of energy shortage consists of forced-compensated part, optional compensated part and abandoned part, as illustrated in Figure 7b.Further explanation on discharge strategy is illustrated in Figure 7.For each time interval, the energy shortage is divided into three stages distributed along the time axis in Figure 7a.Each stage of energy shortage consists of forced-compensated part, optional compensated part and abandoned part, as illustrated in Figure 7b.Forced-compensated part is the part of energy shortage which must be supplied by energy storage (ES) considering the supply requirement in Table 2. On this basis, additionally compensated energy shortage, which is an optional compensated part for ES, helps to increase the supply reliability.The uncompensated part of energy shortage is considered the abandoned part.Herein, an index named "additional compensation ratio" (ACR) is proposed for balancing economic efficiency and reliability.It represents the ratio of optional compensated part for each stage of energy shortage.Optimal ACR is obtained through calculation of energy shortage compensation cost Cl.
l-1 l 2 (1 ) where En represents each stage of energy shortage in a time interval; Eb is the forced-compensated part of En and shall be compensated at all costs.Uncompensated energy El-b in this part will be punished with c2, which has an extremely high value.The abandoned part of energy shortage is considered and calculated with c1, which is relatively small.The components of each part of energy shortage stages are illustrated in Table 5.

Objective Function
In order to obtain a scientific and economic efficient scheme of a microgrid planning, an integrated consideration of both long-term costs-including operation as well as maintenance costand short-term construction investment is included in this research.The objective function of MG source capacity optimization is set to minimize the economic cost considering life cycle cost (LCC) of  Forced-compensated part is the part of energy shortage which must be supplied by energy storage (ES) considering the supply requirement in Table 2. On this basis, additionally compensated energy shortage, which is an optional compensated part for ES, helps to increase the supply reliability.The uncompensated part of energy shortage is considered the abandoned part.Herein, an index named "additional compensation ratio" (ACR) is proposed for balancing economic efficiency and reliability.It represents the ratio of optional compensated part for each stage of energy shortage.Optimal ACR is obtained through calculation of energy shortage compensation cost C l .
where E n represents each stage of energy shortage in a time interval; E b is the forced-compensated part of E n and shall be compensated at all costs.Uncompensated energy E l-b in this part will be punished with c 2 , which has an extremely high value.The abandoned part of energy shortage is considered and calculated with c 1 , which is relatively small.The components of each part of energy shortage stages are illustrated in Table 5.

Objective Function
In order to obtain a scientific and economic efficient scheme of a microgrid planning, an integrated consideration of both long-term costs-including operation as well as maintenance cost-and short-term construction investment is included in this research.The objective function of MG source capacity optimization is set to minimize the economic cost considering life cycle cost (LCC) of the project.In this model, the LCC optimization theory is applied to take into consideration each part of the cost incurred during the entire duration of the MG project, including the capital investment C cap , the operation and maintenance cost C om , the recycling profit of scrapped equipment C r , the pollutant emission compensation cost C p and the energy shortage compensation cost C l , i.e., where N i is the amount of i-th DG/ES unit, c i is the unit price of i-th equipment, k i represents the replacement times of i-th equipment, L i is the service time of i-th equipment, r is the discount rate, c om-i is the annual operation and maintenance cost of i-th equipment, j max is the service lifetime of the project, c r-i is the recycling profit of equipment N p-i is the annual emission amount of i-th pollutant, c p-i is the unit cost of i-th pollutant compensation, c l-i is the annual energy shortage compensation cost of level i-th load.

1.
Constrains of output power of PV, WT and diesel generator: For PV and WT, the limits (especially for minimum values) are associated with the primary source.

2.
Constrains of charge and discharge power of ES devices: where P c-max , P c-min are the maximum and minimum charging power of ES devices.

3.
Constrains of power supply reliability, where the unserved time of level i-th load (T l-i ) shall not exceed its corresponding limit: Constrains of ES discharging time span, for each discharging process, the discharging time span of SC and battery shall be limited according to the proposed strategy: Different ES devices and backup sources possess various switching time characters.As for super capacitor, the switching time can be as short as microseconds (set as 20 ms in this paper).Meanwhile, BESS has a longer switching time compared to super capacitor, set as 1 min in this paper, via integrated consideration of component character and power quality requirement.For diesel generators, a warm-up process is required before supplying regulated power to the load, normally 40% of the rated power.This time varies according to the ambient temperature, ranging from 4 to 8 min, and can be up to 20~30 min under extreme conditions (i.e., −40 • C temperature).The switching time of the diesel generator is set as 5 min in this paper.It is worth noting that the switching time of all power source devices relies on parameters of the control system and power-electronic converters, hence t switch can be adjusted and modified accordingly in specific applications and this research aims to propose a generic approach.

Load Character and DER Data
To verify the proposed optimization strategy and method, a microgrid project with a 20-year service lifetime is selected as an example.Load in this microgrid project consists of level 1 load (industrial load, government department, hospital, telecommunication facility), level 2 load (schools, commercial load), level 3 load (residential load).With a general consideration of power balance in the system, renewable energy resource applied in the grid consists of 14 × 30 kW wind turbines and 48 kWh photovoltaic panels.The typical daily load curves are illustrated in Figure 8.The renewable energy generation curve is illustrated in Figure 9 according to the environment data illustrated in Table 7. commercial load), level 3 load (residential load).With a general consideration of power balance in the system, renewable energy resource applied in the grid consists of 14 × 30 kW wind turbines and 48 kWh photovoltaic panels.The typical daily load curves are illustrated in Figure 8.The renewable energy generation curve is illustrated in Figure 9 according to the environment data illustrated in Table 7.Based on the illustrated data, C om of a typical day can be obtained, and annual costs for operation and maintenance are calculated as 365C om .Annual pollutant compensation and energy shortage compensation cost is obtained similarly.The recycling profit C r is calculated as 5% of the capital investment cost for each type of device, and the scrap value at the ending year of lifetime is discounted into the value at the initial year with r = 3%.

ES Device Capacity Optimization
Particle Swarm Optimization (PSO) algorithm is applied to solve this problem.PSO algorithm is an evolutionary algorithm proposed by Kennedy J. and Eberhart R. in 1995 [45], and attracted extensive attention for its advantages in easy application, high precision and convergence speed.In this case, algorithm parameters are set as: iterations I max = 1000; maximum speed V max = 50; population size m = 20.I max is set to avoid the iterative process being trapped into a loop, and the iterations will stop when no improvement (larger than 10 −5 ) has occurred in the last 100 iterations [46].
The optimized LCC result of this planning scheme is $1,594,127 with ACR = 0.4818; specific results of optimized capacity are shown in Table 8, LCC costs of each device are shown in Table 9 and calculation results of each LCC component are shown in Table 10.The optimizing curve of LCC is illustrated in Figure 10.ACR is the ratio of optional compensated part for each stage of energy shortage.It reflects the improvement on supply reliability while basic requirements on power supply for level 1, level 2 and level 3 load are satisfied.ACR is the ratio of optional compensated part for each stage of energy shortage.It reflects the improvement on supply reliability while basic requirements on power supply for level 1, level 2 and level 3 load are satisfied.
When ACR = 0, I, IV, VII in Figure 7 will be compensated while energy shortage of II, III, VI will be abandoned and calculated into compensation cost in Cl, as illustrated in Figure 11.The calculation result is shown in Figure 12.When ACR = 0, I, IV, VII in Figure 7 will be compensated while energy shortage of II, III, VI will be abandoned and calculated into compensation cost in C l , as illustrated in Figure 11.ACR is the ratio of optional compensated part for each stage of energy shortage.It reflects the improvement on supply reliability while basic requirements on power supply for level 1, level 2 and level 3 load are satisfied.
When ACR = 0, I, IV, VII in Figure 7 will be compensated while energy shortage of II, III, VI will be abandoned and calculated into compensation cost in Cl, as illustrated in Figure 11.The calculation result is shown in Figure 12.ACR is the ratio of optional compensated part for each stage of energy shortage.It reflects the improvement on supply reliability while basic requirements on power supply for level 1, level 2 and level 3 load are satisfied.
When ACR = 0, I, IV, VII in Figure 7 will be compensated while energy shortage of II, III, VI will be abandoned and calculated into compensation cost in Cl, as illustrated in Figure 11.The calculation result is shown in Figure 12.As illustrated, although the investment cost of ES devices slightly decreased (C capSC decreased by 4% and C capBESS decreased by 22%) compared with the optimal planning scheme at ACR = 0.4818, the compensation cost C l significantly increased from $15,746 to $30,526 by 94%.The LCC result increased from $1,594,127 to $1,608,205.
When ACR = 1, all parts of energy shortage in Figure 7 (I~VII) will be totally compensated (III, VI are reduced to zero), as illustrated in Figure 13.The compensation cost C l is theoretically zero.As illustrated, although the investment cost of ES devices slightly decreased (CcapSC decreased by 4% and CcapBESS decreased by 22%) compared with the optimal planning scheme at ACR = 0.4818, the compensation cost Cl significantly increased from $15,746 to $30,526 94%.The LCC result increased from $1,594,127 to $1,608,205.
When ACR = 1, all parts of energy shortage in Figure 7 (I~VII) will be totally compensated (III, VI are reduced to zero), as illustrated in Figure 13.The compensation cost Cl is theoretically zero.The calculation result is shown in Figure 14.As illustrated, the compensation cost Cl is reduced to zero as expected.However, the investment costs of SC and BESS increase approximately 5.8 and 16 times than those of the optimal planning scheme at ACR = 0.4818 respectively.The LCC result increases to $1,774,287.
Based on the above discussion, the proposed ACR is the key to balancing economic efficiency and supply reliability.The relationship between ACR and optimal cost is illustrated in Figure 15 through calculation and fitting.The calculation result is shown in Figure 14.
As illustrated, although the investment cost of ES devices slightly decreased (CcapSC decreased by 4% and CcapBESS decreased by 22%) compared with the optimal planning scheme at ACR = 0.4818, the compensation cost Cl significantly increased from $15,746 to $30,526 by 94%.The LCC result increased from $1,594,127 to $1,608,205.
When ACR = 1, all parts of energy shortage in Figure 7 (I~VII) will be totally compensated (III, VI are reduced to zero), as illustrated in Figure 13.The compensation cost Cl is theoretically zero.The calculation result is shown in Figure 14.As illustrated, the compensation cost Cl is reduced to zero as expected.However, the investment costs of SC and BESS increase approximately 5.8 and 16 times than those of the optimal planning scheme at ACR = 0.4818 respectively.The LCC result increases to $1,774,287.
Based on the above discussion, the proposed ACR is the key to balancing economic efficiency and supply reliability.The relationship between ACR and optimal cost is illustrated in Figure 15 through calculation and fitting.As illustrated, the compensation cost C l is reduced to zero as expected.However, the investment costs of SC and BESS increase approximately 5.8 and 16 times than those of the optimal planning scheme at ACR = 0.4818 respectively.The LCC result increases to $1,774,287.
Based on the above discussion, the proposed ACR is the key to balancing economic efficiency and supply reliability.The relationship between ACR and optimal cost is illustrated in Figure 15 through calculation and fitting.
Through observation, the optimal planning scheme is obtained at ACR = 0.4818, and the cost will increase whether the ACR value is too small or too high.The decreasing tendency on the left side occurs because the exceeded compensation cost is reduced by reasonable investment on ES devices, while excessive investment would lead to over-sizing of ESS and the cost rises on the right side of the curve.Moreover, when ACR exceeds 0.5, the cost will increase sharply.This phenomenon indicates that improving reliability gradually becomes harder.For instance, improving ACR from 0.8 to 0.9 requires a $42,083 increase in the cost, while only $23,636 is needed for improving ACR from 0.6 to 0.7.
It is easy to achieve good performance, but hard to be perfect.This characteristic can provide guidance for decision-makers to obtain optimal planning schemes under different conditions with particular requirements on reliability and economic efficiency.Through observation, the optimal planning scheme is obtained at ACR = 0.4818, and the cost will increase whether the ACR value is too small or too high.The decreasing tendency on the left side occurs because the exceeded compensation cost is reduced by reasonable investment on ES devices, while excessive investment would lead to over-sizing of ESS and the cost rises on the right side of the curve.Moreover, when ACR exceeds 0.5, the cost will increase sharply.This phenomenon indicates that improving reliability gradually becomes harder.For instance, improving ACR from 0.8 to 0.9 requires a $42,083 increase in the cost, while only $23,636 is needed for improving ACR from 0.6 to 0.7.It is easy to achieve good performance, but hard to be perfect.This characteristic can provide guidance for decision-makers to obtain optimal planning schemes under different conditions with particular requirements on reliability and economic efficiency.

Conclusions
In this paper, the character analysis as well as modeling of typical DG and ES devices are completed with special emphasis on ES characters.Load character is considered to obtain a reasonable coordination strategy considering combinations of ES devices and diesel generator, followed by establishing the objective function based on life cycle cost (LCC) calculation.The proposed planning strategy and optimizing model are verified through optimized planning of the energy source capacity of an autonomy microgrid, and is able to ensure power supply for an important load, while maintaining a reasonable economic cost.The relationship between economic cost and the proposed ACR index, which reflects supply reliability condition of the system, will help decision-makers to obtain optimal planning schemes under different conditions with particular requirement on supply reliability.Hopefully, this study will provide suggestions for decision making in future power source capacity planning schemes of autonomy microgrid on isolated islands or in remote areas, and the power source coordination strategy can also be used as a reference in DG integration in power systems as well as smart grids.Further work on microgrid structure planning and evaluation methods is already on the agenda; all of these research works will eventually form an integrated, flexible and coordinated system of a microgrid planning method.

Figure 2 .
Figure 2. Cycle time of lead-acid batteries.DOD, depth of discharge.

Figure 3 .
Figure 3. Economic loss curve of typical load after power outage.

Figure 5 .
Figure 5. Operation status and coordination strategy of power sources in different conditions.(a) Illustration of charging process; (b) Illustration of discharging process.SC, super capacitor.

Figure 6 .
Figure 6.Flowchart of the planning strategy.

Figure 5 .
Figure 5. Operation status and coordination strategy of power sources in different conditions.(a) Illustration of charging process; (b) Illustration of discharging process.SC, super capacitor.

Figure 6 .
Figure 6.Flowchart of the planning strategy.

Figure 7 .
Figure 7. Illustration of discharge strategy with additional compensation ratio (ACR) considered.(a) Stages of energy shortage; (b) Three parts of an energy shortage stage.ACR, additional compensation ratio.

Figure 7 .
Figure 7. Illustration of discharge strategy with additional compensation ratio (ACR) considered.(a) Stages of energy shortage; (b) Three parts of an energy shortage stage.ACR, additional compensation ratio.

Figure 15 .
Figure 15.Relationship between ACR and optimal cost.

Table 1 .
Parameters of various types of energy storage (ES) devices.

Table 1 .
Parameters of various types of energy storage (ES) devices.

Table 2 .
Classification of load based on their importance.

Table 3 .
t interrupt of electrical equipment in the electronics manufacturing industry.

Table 4 .
t interrupt of major types of level 1 load.

Table 5 .
Component of each stage of energy shortage.

Table 5 .
Component of each stage of energy shortage.

Table 9 .
Life cycle costs (LCC) of each device.

Table 10 .
Calculation results of each LCC component.