# Optimal Allocation and Sizing of Capacitors for Distribution Systems Reinforcement Based on Minimum Life Cycle Cost and Considering Uncertainties

V. V. Thang*, N. D. Minh
Thainguyen University of Technology (TNUT), Thainguyen, Vietnam

#### Article Metrics

0
##### Total Statistics:

Full-Text HTML Views: 2863
Abstract HTML Views: 696
##### Unique Statistics:

Full-Text HTML Views: 1034
Abstract HTML Views: 411

open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC-BY 4.0), a copy of which is available at: https://creativecommons.org/licenses/by/4.0/legalcode. This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

* Address correspondence to this author at the Thainguyen University of Technology - Electric Power System Thainguyen, Viet Nam Electric Power System Vietnam, Tel: 0084-915176569; E-mail: thangvvhtd@tnut.edu.vn

## Abstract

### Background:

The power and voltage losses have been greatly concerned in the distribution system planning and operation because low normal voltage, high impedance and load density of the system lead to increases in losses and decreases in bus voltages.

### Method:

Therefore, an optimization model for the distribution system reinforcement improving efficiency of distribution systems by integrating capacitors is proposed. The optimal allocation and sizing of capacitors are determined concurrently with optimal upgrading sizing and timeframe of feeders and substations. The objective function of model is to minimize life cycle cost over reinforcement scheme including investment and operation costs of equipment in distribution systems (including feeders, substations and capacitors) and costs for purchasing energy from the electric market.

### Conclusion:

The typical load curves for each day and season with uncertainties modeled by a Probability Distribution Function (PDF) and Time of Use (TOU) price applied to enhance accuracy of calculation results and suitability for real operation conditions. General Algebraic Modeling System (GAMS) is applied to undertake calculations in a test system.

Keyword: Distribution System, Optimization, Capacitor, Life cycle cost, GAMS, Uncertainties.

## 1. INTRODUCTION

Distribution systems serve as connections between distribution substations and retail customers with low normal voltage, high feeder impedance and load density [1, 2]. Hence, there will be a rise in power and voltage losses that are particularly considered in planning, designing and operating problems of distribution systems.

Capacitors are able to reduce power losses and delay upgrading of feeders and source substations because of reduction in transmission power on devices. Many new technologies in capacitor designs have been introduced with less leakage powers and the price recently. Additionally, the investment in construction and operation is low-cost and installation places are also unlimited [3, 4]. Hence, capacitors can be easily integrated in distribution systems and improve the economic and technical criteria. Many researches have been conducted to choose locations and capacity of capacitors to reduce power and voltage losses as well as operation cost in the distribution systems [5, 6]. The objective function to minimize power losses and voltage limits on load buses are introduced in [7, 8]. Furthermore, costs of devices are also considered in the objective function including total cost of energy losses and capacitor investment in order to enhance accuracy of calculation results [9-11]. Constraints on capacity of equipment (such as feeders and transformers) and bus voltages are utilized to ensure the requirements of operations.

In above researches, capacitors are often considered in problems of seeking optimal setting and sizing in distribution systems. It is assumed that the structures and parameters of the systems remain unchanged. However, the distribution systems are always reformed and upgraded by increasing the section of feeders and the capacity of transformer substations to meet developments of loads. In other words, the peak load is utilized to calculate the operation of distribution systems that demand strongly changes by time in day and season per year. Therefore, the error of calculation can be significantly big, which causes the over compensation in time of low demand and the increase in power loss at this time.

Moreover, the load of the distribution systems depends on demands of customers and weather that have arbitrary characteristics. Hence, the normal probability density function (PDF) is chosen for modeling the uncertainties of the load at each bus [13, 14]. In researches [15-18], the system’s loads are assumed to follow the hourly load shape divided into different levels using a clustering technique with various probabilities. Therefore, the probability forecasting on load distribution becomes a significant issue that notably affect the operation and reinforcement of the distribution systems. However, the optimal compensation calculation of capacitor at the same time with the upgrading of feeders and transformer substations has not been mentioned in above researches. Uncertainties of loads are not validated because the demand in each hour is always set to be unchanged.

In recent years, the electricity price usually determined by Time-of-Use (TOU) is the simplest dynamic price. The objective of TOU price is to encourage the reduction of energy consumption during peak hours and the electricity price varies each day. TOU price is currently based on three periods of energy use as shown in [19-22] as follows: i) Peak - the highest demand; ii) Intermediate - the moderate demand; iii) Base - the lowest demand.

The mathematical optimization is widely applied in planning and reinforcement of the distribution systems with a variety of new models and methods [23]. Georgilall etal. and Tung etal. [24, 25] present models for the planning of the distribution systems with concurrent participation of distributed generators or capacitors. An optimal model in [26] is implemented with the objective function minimizing cost for investment feeder, transformer substation, distributed generator, capacitor and energy consumption from the electrical market in planning period. Similarly [27], introduces two stage planning models aim at decreasing the quantity and time of calculation. The objective function that seeks minimum value of investment and operation cost is used in associated with constraints in two stage in order to ensure the operation of system.

Besides, the objective function based on life cycle cost is also widely utilized in the problems of planning and operation of the distribution systems. The model with the objective function minimizing the life cycle cost in entire distribution system planning period is introduced in [28-29]. These models aim at finding the best distribution system planning scheme to maximize overall benefits and costs in the life cycle of the system. Constraints on maximum allowable capacity of equipment (feeder and transformer) and balanced power flow of system are considered to guarantee the economic and technical indicators. Similarly, the life cycle cost is persistently used in objective functions in [30-31] to determine the optimal transformer capacity. The model integrates the dynamic change of money with time, environmental cost, preventive test cost, sensitivity coefficient and the parameters of transformers.

In this research, an optimal model combining the uncertainties of demand described by a Probability Distribution Function (PDF) will be proposed. The optimal displacement, sizing and installation period of capacitors are simultaneously determined with optimal sizing and timeframe of feeders and transformer substations. The objective function is minimum life cycle cost of the planning scheme as well as the technical constraints are utilized to guarantee operation of the distribution systems and capacitor. The electricity price in the system and the electrical market is TOU price.

Next sections of the paper are organized as follows: an uncertainty modeling of demand is introduced in Section 2. A mathematical model with objective function and constraints are presented in Section 3. Section 4 shows calculated results from the IEEE 9-bus test distribution system. Finally, the conclusion is demonstrated in Section 5.

## 2. MODEL OF LOAD UNCERTAINTIES

The uncertainties of load at each bus are stochastic parameters often modeled by the normal PDF as shown in equation (1) with mean of the distribution μ, standard deviation σ and variance σ2 [13-16].

 (1)

In other words, there are a number of states of load in each time segment. The selection of the number of states plays an important role because a small number of states may affect badly accuracy of calculations while a large number increases the complexity of the problem and calculation time. Hence, in this research, the requirement of the distribution systems is hypothesized to follow the hourly load characteristics of the IEEE-RTS using a clustering technique [17-18]. In each time segment, the load is divided into ten levels (states) with different probabilities to guarantee a reasonable balance between accuracy and speed as shown in Table 1.

Table 1. Load data in states.
1 100 0.01 6 58.5 0.163
2 85.3 0.056 7 51 0.163
3 77.4 0.1057 8 45.1 0.0912
4 71.3 0.1654 9 40.6 0.0473
5 65 0.1654 10 35.1 0.033

## 3. THE MATHEMATICAL MODEL

The uncertainty of the electrical demand considered in this research causes an increase in the complication of the problem. The mathematical model must take an uncertain load including many states into account with the different values of load and probabilities in each state. The proposed model simultaneously examines cost for upgrading of feeders, transformer substations and investment capacitors. Furthermore, the active and reactive power flow constraints in each state are proposed to calculate power flow and losses of system in all operation states with high accuracy. The model applies binary variables to represent the cost characteristic as well as investment and upgrading decision of equipment. Details of the model are presented as follows:

### 3.1. Objective Function

In the study, an objective function minimizing the life cycle cost of the investment project during calculation period is proposed as written in equation (2).

 (2)

where: the life cycle cost is calculated at the base year by discount rate r; the upgrading cost of feeders CFt; the upgrading cost of transformer substations CTt; the investment cost of capacitors CCt; the loss cost in the capacitor itself CLt and the energy purchased cost CEt and residual value of equipment at the end of the reinforcement period.

The upgrading cost of feeders including fixed capital cost CF0 and the variable capital cost CF is presented in equation (3) with the binary variable αij.t to represent the cost characteristic and decision variable for feeder upgrading Fij, t.

 (3)

where N is total number of system buses; T is total number of planning year; i and j are indices for all buses in the system; t is index for examination year.

Similarly, the upgrading cost of transformer substation CTt including fixed capital cost CT0 and variable capital cost CT is presented in equation (4) with a binary variable γi.t to perform the cost characteristic, transformer upgrading decision variable ΔSTi, t-NT. NT is total number of transformer substation buses.

 (4)

The equation (5) represents a new investment cost of capacitor in year t including fixed capital cost CC0 and variable capital cost CC with the binary variable βi,t to represent the cost characteristic and decision variable for capacitor investment Qi, tC [12]. Qi, tC is affected by voltage and determined as (10). NC is total number of capacitor installation buses.

 (5)

The loss cost in the capacitor itself integrated uncertainties in year t are shown in equation (6) including decision variable for investment capacitor Qi, tC.

 (6)

where kC is the active power loss factor of the capacitor itself; ds is the total number of days in each season; NH is total number of hours in each day; h is index of each hour and Ns is total seasons in a year.

The energy purchased cost CEt considering uncertainties of loads is presented in equation (7) with active power Pi, t, s, h, g and reactive power Qi, t, s, h, g purchased through transformer substation.

 (7)

where PhP, PhQ are active power purchased price at hour h; Kp is increase coefficient of electrical price in year t; The load probability in state g; λg; NG is total number of states; g is index of each state.

The life cycle of each equipment is different from each other and the equipment is not concurrently invested so at the end of the calculation period the residual value of them is presented in equation (8). Hence, the residual value is the current value calculated at the base year in objective function.

 (8)

where tinvF, TLF are installation time and life cycle of feeder; tinvT, TLT are installation time and life cycle of transformer substation; tinvC, TLC are installation time and life cycle of capacitor.

### 3.2. Constraints

The technical constraints are utilized to guarantee the operability of the distribution systems and capacitor including alternating current (AC) power flow, feeder upgrading section, transformer substation upgrading capacity, limitation of bus voltage and capacitor size.

#### 3.2.1. Constraints for Power Flow

Considering uncertainties in this research will lead to calculating the power flow for every state of load. Therefore, an AC nonlinear power flow model supplemented a state index g is written in equation (9). The effects of reactive power on the calculation of power and voltage losses are interested in the model because the accuracy of calculation results is improved.

 (9)

where PDi,t,s,h,g, QDi,s,t,h,g are active and reactive power of load i; Ui, s, t,h,g, δi,s,t,h,g are module and angle of bus voltage; Yij, t, θij, t are magnitude and angle of admittance matrix element formulated by the feeder impedances in year t; QCi,s,t,h,g is reactive power of capacitor at each state as equation (10).

 (10)

#### 3.2.2. Limits on Dimension of Capacitors

The participation of capacitor changes the power flow and affects the technical and economic indicators of distribution systems. To ensure annual installation of capacity corresponding to the equipment parameters, the size of capacitor invested at each bus in each year needs to follow constraints in (11). Additionally, the planning period is short-term; therefore, at each load bus, only a capacitor system is chosen in overall the plan horizon to reduce the installation cost as shown in equation (12).

 (11)
 (12)

where Qmax, Qmin are maximum and minimum capacity limit of capacitor; βi,t is binary variable.

#### 3.2.3. Upgrading Section Constraints of Feeder

The power carrying capability of feeders is imposed by thermal limits taken into the consideration of new feeder investments. The constraints for feeder upgrading are represented in equation (13) with capacity supplemented at each year, ΔSij,tF. ΔSij,tF is set when the capacity value is greater than the capacity limit used at year t-1. Thus, the feeder capacity needs to be upgraded at year t to meet in order to supply power of the load as equation (14).

 (13)
 (14)

where ΔSmaxF, ΔSminF is maximum power running in feeder at year t; ΔSmaxF are maximum and minimum capacity limit added of feeders; Sij, t*F is standard size of feeder and αij,t is binary variable.

#### 3.2.4. Addition Capacity Constraints for Substation

In order to use the maximum of existing substations capacity and satisfy upgrading power corresponding to the equipment parameters. The constraint (15) is utilized to set substation sizes added in year t, ΔSi,tS, with binary variable γi,t.

 (15)

where SmaxS, ΔSminS are capacity limit added of transformer substation; Si,t*S is standard size and Si,tmax is maximum capacity running in the transformer at year t.

#### 3.2.5. Constraints on Bus Voltages

The distribution systems are connected to the transmission systems by transformer substations so the voltages at substation bus are usually stabilized and assumed to be constant. At load buses, the voltages are usually changed according to parameters of the equipment and the value of the loads. Thus, the limits on load bus voltage are shown in equation (16).

 (16)

where NL is total number of load buses; Umin, Umax are maximum and minimum limit of bus voltage.

The proposed model is a mixed integer nonlinear programming (MINLP) calculated in GAMS environment with BONMIN solver to find out an optimal solution [36].

## 4. CALCULATIONS AND DISCUSSIONS

### 4.1. Parameters of Distribution Systems and Assumptions in Analysis

The IEEE 9-bus radial structure with 22kV connected to grid by transformer substation is utilized to investigate the proposed model in Fig. (1) [32, 33]. Parameters changed to match the problem are presented in appendix. The total active power and reactive power at the base year are 15,024.2kW and 10,641.8kVAr, respectively. The hourly load demand of four seasons is computed by fact of peak demand as shown in Fig. (2). The uncertainties of load at each bus are stochastic parameters as introduced in section 2.

 Fig. (1). Diagram of the test distribution system.

 Fig. (2). The hourly load demand for four seasons.

The following assumptions are utilized in this research including:

• The calculation period is selected to be 5 years with the annual developing rate of demand 10% per year. The typical characteristics of demand at the entire load are set to be the same as in Fig. (2). The electrical energy price is TOU price as shown in Fig. (3) [20, 21, 34].
• The fixed and variable upgrading costs of transformer substation are 0.2M$and 0.05M$/MVA, respectively. Similarly, the upgrading costs of 22kV feeders consist of 0.15M$/km and 0.001M$/MVA.km. The assumption lifetime of transformer substation and feeder is 20 years. The areas for upgrading of substation transformers and feeders are assumed to be unlimited.
• The fixed and variable investment cost of capacitors assumed in this research are 1000$for each installed system and 20$/kVAr as shown in [7, 11, 12]. The active power loss factor of the capacitor itself kc is determined 0.02W/kVAr and the life of capacitor are 20 years [35]. The capacitors can be selected to install at all of load locations because of quick installation and small spaces of occupation.
• The voltage at load buses is allowed to change from 0.9pu to 1.1pu in order to guarantee the operation of equipment. At the substation bus, the voltage usually is stabilized and assumed to be 1.05pu.
 Fig. (3). Energy prices.

Table 2. Calculated Cases.
Case Certainty Uncertainty Capacitor
1 x
2 x
3 x x
4 x x

### 4.2. Results

The efficiency of proposed model is investigated by four cases as listed in Table 2. The plan of distribution systems is undertaken to determine the upgrading capacity and time of feeders, substations with fixed and stochastic load in each hour in case 1 and case 2, respectively. The influence of capacitor to distribution system planning is calculated in case 3 that system parameters are similar case 1. In case 4, the stochastic demand is simultaneously considered with participation of capacitor in planning problem.

The calculation results with the above assumption parameters determined upgrading time and capacity of feeders and substations as well as investment time and capacity of capacitor as shown in Table 3. All of four cases, the substation capacity is selected 10MVA but the upgrading time is different. The transformer substation is upgraded at 4th year in case 1 and 2 while case 3 and 4 is deferred to 5th year because of the capacitor reducing power flow in this case. Similarly, in cases 1 and 2, feeders 1-2 and 1-3 need to be upgraded 5th year with capacity 8MVA while remaining feeders are not upgraded during computation period. Because of the support of capacitor in case 3 and case 4, the feeder 1-3 is deferred and is not upgraded in both cases while the only feeder 1-2 is upgraded at case 4.

Table 3. The investment and upgrading decided of equipment.
No Equipment Bus Base parameter Upgraded or Invested Time (year)
1 2 3 4 5 1 2 3 4 5 1 2 3 4 5
Case 1, 2 Case 3 Case 4
feeder i-j (MVA)
1-2 6.67 8 8
1-3 6.67 8
substation i (MVA)
1 10 10 10 10
3 Capacity invested of
Capacitor i (MVAr)
2 0 0.5
3 0.5 0.5
6 0 0.5
7 0 0.5 0.5
Table 4. The investment and upgrading decided of equipment.
`
No Cost Case 1 Case 2 Case 3 Case 4
1 Total life cycle cost (M$) 32.46 20.61 32.22 20.03 2 Feeder and substation upgrading cost (M$) 1.97 1.97 0.53 0.98
3 O&M and electrical energy cost (M$) 32.36 20.03 32.16 19.98 4Investment capacitor cost (M$) 0 0 0.09 0.04

In this research, the effect of investment project is evaluated by life cycle cost during examination period. The calculation results in Table 4 show that the efficiency of distribution systems is always enhanced with participation of capacitor. The life cycle cost of project in case 3 in comparison with case 1, the fixed load parameters, reduces 0.24M$, the corresponding 0.74%. Similarly, the life cycle cost of project for investment the capacitor in case 4 compared with case 2, the uncertainties load parameters, decreases 0.58M$ in proportion to 1.8%. The cost for investment and upgrading equipment also reduces from 1.97M$of case 2 to 1.02M$ of case 4. Moreover, stochastic load with many states and different probabilities considered leads to reduce the electrical energy purchased from grid. The comparison between case 1 and case 2 shows that the energy cost reduces to 12.33M$, the corresponding 38.10%. Similarly, the energy cost also decreases to 12.18M$ in proportion to 37.87% when the case 4 is compared with case 3.

 Fig. (4). Minimum voltage at buses.

The maximum power loss is determined at state that the demand is maximization present as shown in Fig. (5). In cases 1, 2 and 4, this value is 4.56% at the 1st planning year and increases 0.51% per year at next year depending on the rise of load demand. However, at the 5th year, the maximum power loss is 6.43% and only increases 0.33% compared with the 4th year in case 1 and case 2. This decrease is made due to feeders 1-2 and 2-3 upgraded decrease the resistors of the distribution system. In case 4, this loss is only 6.23% and increases 0.13% compared with the 4th year because the capacitor invested reduces power flow of feeders. Similarly, the capacitor installed in case 3 supported power loss at just 1st year with decrease 0.21% in comparison with case 1. The maximum difference of maximum power loss is 0.42% at the 5th year.

 Fig. (5). Maximum power loss.

Moreover, the electrical energy loss always reduces during calculation period when considering to capacitor and uncertainties of loads as represented in Fig. (6). The electrical energy loss in case 3, the capacitor invested, is reduced 0.11% at just the 1st year in proportion to 7609.47kWh. This value increases 0.02% per year in 2nd, 3rd and 4th year. Nevertheless, it only is 0.13% in the 5th year because the feeders upgraded of case 1 decrease the resistors of the distribution system. Similarly, the capacitor installed at 4th year in case 4 decreases 0.05% of energy loss in proportion to 2150.05kWh compared with case 2. In 5th year, the capacity of capacitor is added 0.5MVAr so the electrical energy loss reduces 0.1%, the corresponding 4738.10kWh. Besides, the electrical energy loss in case 2 in comparison with case 1 always reduces from 1.02% to 1.78%. The maximum difference when the comparison between case 4 and case 3 is 1.66%. This result is made because of the important effect of the stochastic demand with less probabilities and load demand fact.

 Fig. (6). Total energy loss of distribution system.

## CONCLUSION

This paper presents a generalized formulation of distribution system reinforcement integrating the capacitor and uncertainties of loads. The problem is scheduling the upgrading equipment of distribution systems (feeders and transformer substation) together with installation capacitor to minimize life cycle cost of the investment project. The proposed model is based on mixed integer nonlinear programming (MINLP) with binary variables utilized to perform the nonlinear cost characteristics and installation size of the equipment. The proposed formulation is verified by the IEEE 9-bus radial structure with uncertainties of loads and TOU prices. The result of simulation shows that a significant saving of life cycle cost can be achieved in comparison with traditional schemes. Only upgrading of the feeders and transformer substation is demanded. The reduction of power loss and electrical energy loss as well as upgrading delay of equipment in during calculating period could be seen as an effect of capacitor and stochastic load to distribution system reinforcement. Moreover, the voltage at all buses is always ensured in technical requirements in all operation modes.

### APPENDIX

See Table (5 and 6).

No Bus PD0 (kW) QD0 (kVAr)
1 1 - -
2 2 2551.0 1857.7
3 3 2497.7 1935.1
4 4 2656.8 1616.6
5 5 2570.4 1542.5
6 6 1425.6 1110.2
7 7 1182.2 926.9
8 8 1041.8 798.2
9 9 1098.7 854.8
Total 15024.2 10641.8
Table 6. Data of Feeder Parameters.
No Bus i - Bus j Smax.ij
(MVA)
Lij
(km)
Rfij
(W)
Xfij
(W)
1 1-2 6.67 3.3 2.55 1.42
2 1-3 6.67 6.2 4.79 2.66
3 1-4 6.67 5.3 4.10 2.27
4 1-5 6.67 4.5 3.48 1.93
5 2-6 6.67 2.7 2.09 1.16
6 3-7 6.67 3.2 2.47 1.37
7 4-8 6.67 4.2 3.25 1.80
8 5-9 6.67 2.2 1.70 0.94

Not applicable.

### CONFLICT OF INTEREST

The authors declare no conflict of interest, financial or otherwise.

## ACKNOWLEDGEMENTS

The study is supported by Thai Nguyen University of Technology (TNUT), Thai Nguyen, Viet Nam.

## REFERENCES

 [1] John D. McDonald, Bartosz Wojszczyk, Byron Flynn, and Ilia Voloh, Distribution Systems, Substations, and Integration of Distributed Generation In: Encyclopedia Sustainability Sci Technol. Springer 2013. [2] Lee Willis H. Power Distribution Planning Reference Book. In: 2004. crc press 2004 [3] Jayalakshmi M, Balasubramanian K. Simple Capacitors to Supercapacitors-An Overview Inter J Electr Sci 2008; 3: 1196-217. [4] Li J, Cheng X, Shashurin A, Keidar M. Review of Electrochemical Capacitors Based on Carbon Nanotubes and Graphene. Graphene 2012; 1(01): 1. [5] Carlisle JC, El-Keib AA. A review of capacitor placement techniques on distribution feeders In: . Pro.: the 29th Southeastern Symposium on System Theory 1997; pp. 359-65. [6] Aman MM, Jasmon GB, Bakar AH, Mokhlis H. M. Karimi, “Optimum shunt capacitor placement in distribution system: A review and comparative study. Renew Sustain Energy Rev 2014; 30: 429-39. [7] Masoum MA, Ladjevardi M. Optimal Placement, Replacement and Sizing of Capacitor Banks in Distorted Distribution Networks by GA. IEEE Trans Power Deliv 2004; 19(4): 1794-801. [8] Dixit M, Kundu P, Jariwala HR. Optimal Allocation and Sizing of Shunt Capacitor in Distribution systemsfor Power Loss Minimization IEEE Students Con. on Electrical. In: . IEEE Students' Confrence on Electrical, Electronics and Computer Scienc 2016; pp. 1-6. [9] Eajal AA, El-Hawary ME. Optimal Capacitor Placement and Sizing in Unbalanced Distribution Systems With Harmonics Consideration Using Particle Swarm Optimization. IEEE Trans Power Deliv 2010; 25(3): 1734-41. [10] Abou El-Ela A A. Optimal Sitting and Sizing of Capacitors for Voltage Enhancement of Distribution Systems In: 50th International Universities Power Engineering Conference. 2015; pp. 1-6. [11] Tilakul K, Buasri P. Capacitor Location and Size Determination to Reduce Power Losses of a Distribution Feeder in Lao PDR Int J Comput Electron Eng 2012; 4(1): 32. [12] Dixit M, Kundu P, Jariwala HR. Incorporation of distributed generation and shunt capacitor in radial distribution system for techno-economic benefits. Eng Sci Technol Int J 2017; 20: 482-93. [13] Zhang T, Song X, Meng X, Yu J, Chen X. A probabilistic load model based on chi-square method for distribution network In: . 2nd IET Renewable Power Generation Conference (RPG 2013) 2013. [14] Liu Z, Wen F, Ledwich G, Member S. Optimal Siting and Sizing of Distributed Generators in Distribution Systems Considering Uncertainties. IEEE Trans Power Deliv 2011; 26(4): 2541-51. [15] Atwa YM, El-Saadany EF, Salama MM, Seethapathy R. Optimal Renewable Resources Mix for Distribution systemsEnergy Loss Minimization IEEE Transactions on Power Sytems 2010; 25(1): 360-70. [16] Neyestani N, Damavandi MY, Shafie-khah M, Chicco G, Catalão JP. Stochastic Modeling of Multienergy Carriers Dependencies in Smart Local Networks With Distributed Energy Resources. IEEE Trans Smart Grid 2015; 6(4): 1748-62. [17] Pinheiro JM, Dornellas CR, Melo AC. Probing the new IEEE reliability test system (RTS-96): HL-II assessment. IEEE Trans Power Syst 1998; 13: 171-6. [18] Y.M. Atwa, E.F. El-Saadany. Probabilistic approach for optimal allocation of wind based distributed generation in distribution systems. ET Renew Power Gener 2011; 5(1): 79-88. [19] Wong S, Bhattacharya K, Fuller JD. Electric power distribution systemsdesign and planning in a deregulated environment. IET Gener Transm Distrib 2009; 3(12): 1061-78. [20] Meng F-L, Zeng X-J. An Optimal Real-time Pricing Algorithm for the Smart Grid: A Bi-level Programming Approach", In: 2013 Imperial College Computing Student Workshop, 2013. [21] Shao S, Zhang T. "Impact of TOU Rates on Distribution Load Shapes in a Smart Grid with PHEV Penetration." In: IEEE/PES Transmission & Distribution Conference & Exposition (T&D), 2010. [22] Thanhtung HA, Zhang Y, Thang VV, Huang J. Energy hub modeling to minimize residential energy costs considering solar energy and BESS J Mod Power Syst Clean Energy 2017; 5(3): 389-99. [23] Aman MM, Jasmon GB, Bakar AHA. Optimum Capacitor Placement and Sizing for Distribution System Based on an Improved Voltage Stability Index, N.3 Int Rev Elect Eng (IREE) 2012; 7: 4622-30. [24] Georgilakis PS, Hatziargyriou ND. A review of power distribution planning in the modern power systems era: Models, methods and future research. Electr Power Syst Res 2015; 121: 89-100. [25] Tung NS, Chakravorty S. Optimized Power Distribution Planning A Review Int J Electron Elect Eng 2014; 2 [26] Algarni AA, Bhattacharya K. "A Novel Approach to Disco Planning in Electricity Markets: Mathematical Model", In: Power Systems Conference and Exposition, pp. 1-5, 2009. [27] Wong S, Bhattacharya K, Fuller JD. Electric power distribution systemsdesign and planning in a deregulated environment. IET Gen Trans Distribut 2009; 3(12): 1061-78. [28] Su H, Zhang J, Liang Z, Niu S. "Power Distribution Network Planning Optimization Based on Life Cycle Cost", In: China International Conference on Electricity Distribution, pp. 1-8, 2010. [29] Liu L, Cheng H, Ma Z, Zhu Z, Zhang J, Yao L. "Life Cycle Cost Estimate of Power System Planning", In: International Conference on Power System Technology (POWERCON 2010), pp. 1-8, 2010. [30] Yi-Song JR, Long-Li JQ. "Comparative study Of 10kV distribution network structure based on Life Cycle Cost analysis", In: China International Conference on Electricity Distribution (CICED 2012), pp. 1-5, 2012. [31] Qiao G-H, Zheng R, Lei S, Bo Z, Xing J-G, Zhan X-L. "A New Life Cycle Cost Model of Power Transformer and Its Comprehensive Sensitivity Analysis", In: International Conference on Power System Technology (POWERCON 2014), pp. 1342-1348, 2014. [32] Walid El-Khattam YG, Salama M M A. An Integrated Distributed Generation Optimization Model for Distribution systemsPlanning. IEEE Trans Power Syst 2005; 20(5): 1158-65. [33] Khalkhali K, Abapour S, Moghaddas-Tafreshi SM, Abapour M. Application of data envelopment analysis theorem in plug-in hybrid electric vehicle charging station planning. IET Gener Transm Distrib 2015; 9(7): 666-76. [34] Brahman F, Honarmand M, Jadid S. Optimal electrical and thermal energy management of a residential energy hub, integrating demand response and energy storage system. Energy Build 2015; 90: 65-75. [35] Solutions for Power factor Correction at Medium Voltage CIRCUTOR, S.A. http://www.circutor.com [36] Rosenthal RE. A user’s guide. Washington: Development Corporation 2015.