Optimal Distributed Generator Allocation Method Considering Voltage Control Cost

Up till now, the high penetration of intermittent distributed generation (DG) has posed great challenges to the planning and operation of the grid. To achieve the best balance between economic cost and acceptable capacity of DG, this paper proposes a new integrated planning method of the active distribution network while considering voltage control cost. Firstly, characteristics of decentralized and centralized voltage control methods were analyzed. The technical frameworks, voltage control strategies and economical models of different voltage control systems were put forward. Then, an integrated planning model with objectives to minimize the comprehensive cost and maximize clean energy utilization under the constraint of maintaining acceptable voltage was implemented. Simulations were conducted using the Multi-objective Differential Evolution Algorithm (MODE). IEEE 33-bus test systems were employed to verify the effectiveness of the proposed method. The results demonstrate that the proposed approach is able to connect larger distributed generators and decrease the economic cost of Distribution Network Operators while maintaining voltage within the statutory limits.


Introduction
Fossil fuels exhaustion and potentially environment problems have seen a growing increase in the usage of clean energy [1].The connection of distributed wind generation (DWG), photovoltaic energy (PV) and run-of-river small hydropower (SHP) to distribution networks can improve the reliability and flexibility of distribution network.However, the output of distributed generation (DG) is usually influenced by climate environment, which is intermittent, uncertain and fluctuant.Therefore, power system operation will be affected [2,3], and the impact is closely related to the location and capacity of DGs [4][5][6].With uncertain variation of DG output and load, it may lead to some serious problems, such as bidirectional horizontal load flows, greater voltage fluctuation and serious voltage quality problems [7,8].Voltage rise is a significant constraint for increasing the share of sustainable sources.Therefore, power engineers nowadays are facing new challenges in both the planning and operation of power distribution systems.
There have been many studies on the negative impacts of DG integration.Microgrid can provide some solutions to DG integration by the means of self-control, protection and management.However, the complex equipment and high cost restrict its large-scale popularization [9].In the 2008 CIGRE, the theme of operation and development of active distribution system (ADS) was proposed by C6.11 project team of distribution and distributed generation special committee (C6), aiming at achieving effective management of power flow by controlling DG.In addition, many valuable studies on optimal allocation of DGs in ADS have been done.Harrison and Wallace [10] proposed a planning method based on optimal power flow algorithm, in which increasing the capacity of DG was taken as the objective function.Borges and Martins [11] adopted the bi-level scenario programming given that ADS has the ability to control voltage.Zhang et al. [12] took the uncertainties of DG and load into account, and chance-constrained programming method was adopted to establish bi-level optimal allocation model.
The studies above lay the foundation for further work, but some important issues to study still remain.The planning of ADS, which has flexible network structure and integrated control system, is more complex than that of traditional distribution network.When research is conducted on the economic analysis of ADS planning, the costs of its information, control and communication systems are ignored.Recently, the concept of cyber physical system for power grid (power CPS) has become a research subject of great concern worldwide.The cyber physical system (CPS) is designed to realize the interoperability and deep integration of physical and cyber systems, so that it can obtain better operating effect and performance level beyond the traditional application system [13].Besides the conventional physical power system, power CPS emphasizes that the impacts of related information, control and communication systems should be analyzed for it has strong relationships with the operation situation of the grid.With the rapid development of power CPS, the control system and its supporting systems (information system, communication system and automation system) will become more and more complex, and the total cost is closely related to the complexity of them.In addition, the control cost and ability of different control strategies are different.Therefore, to give a reasonable planning scheme of DG allocation, the cost of relevant control system and the control effect of control strategies should be included [14].
The purpose of this study is to provide a novel DG allocation methodology that considers both the cost and effects of the voltage control systems.The basic idea is to optimize the installed capacity of DG units considering both the benefits of Distribution Network Operators (DNOs) and the penetration level of DG units.Meanwhile, two categories of voltage control approaches aimed at minimizing the voltage deviation from a reference value are proposed and summarized.Their technical frameworks, specific voltage control strategies and economical models are analyzed.
The main contributions and differences can be summarized as follows: (1) Compared to other related works, an integration planning methodology combines DG units' optimization and voltage control is proposed; (2) two main categories of voltage control strategies with DG units involved actively are proposed and the cost and control effect of control systems are analyzed in detail; (3) in the model of DG capacity optimization, the cost of relevant control and communication systems are included, which is a non-ignorable part with the development of ADS; (4) the volatilities, uncertainties and errors of consumption and generation are considered and stochastic models of DG generation and load demand are put forward in the case study; (5) the control effect of voltage regulation is included in the simulation, which gives the operators a visualized display of the control strategy ; and (6) the influence of certain voltage control strategy on the optimal acceptable capacity of DG is shown, which can help make a decision on which control strategy to take in a specific situation.
The structure of this paper will be broken down as follows.The voltage control system in distribution network will first be shown, followed by the method of DG capacity optimization.Then, the case study carried out on IEEE 33-bus test systems is performed.Finally, the conclusions of this paper will be drawn.

The Operation Mode of DG
The interfaces of DG are generally classified into three forms [15]: synchronous generators, asynchronous generators, and DC/AC or AC/DC converter.There are usually two types of DG operation modes including the constant power factor and constant voltage modes, which can be taken as PQ bus and PV bus in power flow calculation, respectively.
(1) Power Factor Control Mode (PFC) In power factor control mode, the P/Q ratio of a generator is kept constant, with the reactive power following the variation of real power.Traditionally, in order to ensure the availability of DG unit's full real power output, the power factor of DG remains (near) uniform.The generator bus will serve as a PQ bus.
(2) Voltage Control Mode (VC) In voltage control mode, automatic voltage regulator (AVR) of generator ensures voltage constant by changing excitation [15].When the output of DGs increases, AVR needs to regulate the field current of generator to keep voltage constant.Under this situation, the reactive power output will be decreased or be absorbed.This control strategy is relatively complex to deal with.The generator bus will serve as a PV bus.However, the reactive power after modified by equations may exceed limits.In this case, the PV bus should be converted to a PQ bus.

Category of Voltage Control in ADS
In conventional distribution networks, some data are insufficient and difficult to measure.A nine-zone diagram control strategy is used by VQC devices to control voltage [16].With the development of ADS, DGs provide many controllable elements in voltage control and new control strategies.Now, some researches propose that centralized and decentralized voltage control can be used in ADS to control voltage in a more active manner [17].In [18], the centralized voltage control strategy is defined which uses a wide range of communication system to coordinate different devices to control voltage.Although the control effect is good, the control system is relatively complex and high-cost.In [19], decentralized control strategy is realized by controlling DGs, capacitors, on-load tap changer (OLTC) and other devices locally to ensure the voltage of monitoring points within limits.The coordination of this strategy is relevantly poor, and some devices may be activated frequently.However, it is cost-effective because a communication system is not required.Moreover, it can not only minimize the power losses but also increase the generation capacity thanks to its flexibility [20].

Controllable Elements in Voltage Control System
Controllable elements participating in ADS voltage control: OLTC, secondary capacitor in substation C, line capacitor C 1 and DG units.The operation of the tap changer is limited to its tapping limits and capacity.The action times of OLTC and C are limited.There are several modes of voltage control with DG involved: Power Factor Control (PFC), Voltage Control (VC), Power Factor-Voltage Control (PFC-VC) and Generation Curtailment.PFC depends on a certain limit of generation connected to the system.VC is disruptive to the network devices such as OLTCs.PFC-VC method combines the behavior of the generator's operation in two modes namely, PFC and VC.Generation Curtailment is the last resort if other methods are not successful.

Decentralized Voltage Control System
In decentralized voltage control method, local information is used to control voltage at a particular bus independently.This method is widely adopted in China where measurement, optimization and communication methods are limited [21].
(1) Characteristics of the Method This is a voltage control strategy based on self-information to improve overall network performance.Thus, it would not require extensive deployment of sensors and communications equipment.Moreover, it can have positive effects on both power losses decrease and generation capacity increase thanks to its flexibility [20].It is applicable in remote mountain areas, but it cannot achieve global optimization.
(2) Economical Model Only local information is needed, so decentralized voltage control system is composed of locally automatic controller and corresponding communication equipment.After analyzing its economic characteristics, the construction cost of the voltage control system composes the cost of its controller and the construction cost of its communication system.It can be described as follows: C tx " where C vo.de is the construction cost of decentralized voltage control system; C local is the cost of local automatic controller; C tx is the construction cost of communication system; λ 1 , . . ., λ 5 are, respectivel,y the length of fiber, the number of EPON-OLT, EPON-ONU, GPRS terminal and integrated network management equipment; β 1 , . . ., β 5 are, respectively, comprehensive unit price of fiber, EPON-OLT, EPON-ONU, GPRS terminal and integrated network management equipment; and C 0 is the cost of construction and management. (

3) Control Strategy
To make the system more complete, this paper involves day-ahead optimization in decentralized voltage control system.The strategy is divided into day-ahead and real-time scale.
Day-ahead control adopts static optimization method.Based on load forecast and DG output prediction, the upper and lower voltage limits of OLTC and operation voltage of secondary capacitor in substation are computed [22].After considering operation times, the results of optimization are sent down to dispatchers.
Real-time control adopts decentralized voltage control mode.The dynamic adjustment features of DG can be used to adjust the power factor of generator.Prediction deviations of load and coming water volume can be balanced.Both the day-ahead and real-time control can be achieved without communications.Thus, this decentralized voltage control system needs no communication systems.
To make it clear, the simulation process of this strategy is shown in Figure 1: Sustainability 2016, 8,193 Only local information is needed, so decentralized voltage control system is composed of locally automatic controller and corresponding communication equipment.After analyzing its economic characteristics, the construction cost of the voltage control system composes the cost of its controller and the construction cost of its communication system.It can be described as follows: where Cvo.de is the construction cost of decentralized voltage control system; Clocal is the cost of local automatic controller; Ctx is the construction cost of communication system; 1, …, 5 are, respectivel,y the length of fiber, the number of EPON-OLT, EPON-ONU, GPRS terminal and integrated network management equipment; 1, …, 5 are, respectively, comprehensive unit price of fiber, EPON-OLT, EPON-ONU, GPRS terminal and integrated network management equipment; and C0 is the cost of construction and management. (

3) Control Strategy
To make the system more complete, this paper involves day-ahead optimization in decentralized voltage control system.The strategy is divided into day-ahead and real-time scale.
Day-ahead control adopts static optimization method.Based on load forecast and DG output prediction, the upper and lower voltage limits of OLTC and operation voltage of secondary capacitor in substation are computed [22].After considering operation times, the results of optimization are sent down to dispatchers.
Real-time control adopts decentralized voltage control mode.The dynamic adjustment features of DG can be used to adjust the power factor of generator.Prediction deviations of load and coming water volume can be balanced.Both the day-ahead and real-time control can be achieved without communications.Thus, this decentralized voltage control system needs no communication systems.
To make it clear, the simulation process of this strategy is shown in Figure 1:  ‚ Day-ahead optimization: Make a plan for the distribution network according to the data of typical day; and obtain the action sequence of OLTC, capacitor banks and other devices based on the results of optimization.

‚
Real-time decentralized voltage control: According to the stochastic models of DG generation and load, start up the voltage control system when voltage of measured node exceeds limits and adjust power factor of DGs successively for under-excited operation.

‚
Stop voltage control operation when the power factor of the last DG unit reaches cosφ min (capacitive) but voltage remains unqualified, which means this voltage control strategy unable to adjust the voltage to normal level.

Centralized Voltage Control System
Centralized voltage control is the most effective way to manage and control the operation of the whole ADS.With the development of communication technology, SCADA (Supervisory Control And Data Acquisition)-centered centralized voltage control can optimize operation based on global information.
(1) Characteristics of the Method This method is based on optimization theory.Thus, it can make full use of various voltage control devices to regulate globally, which is able to achieve overall optimization.However, the application of centralized control strategies to the existing networks faces several drawbacks: in addition to the heavy investments necessary for devices and control systems, all centralized approaches require a highly reliable communication channel through the overall distribution network [23].
(2) Economical Model Centralized voltage control system is composed of Active Distribution Network Management System (ADMS), area coordination controller and locally automatic controller [24].After analyzing its economic characteristics, the construction cost of the voltage control system composes ADMS cost, area coordination controller cost, locally automatic controller cost and its communication system construction cost.It can be described as follows: where C vo.de is the construction cost of centralized voltage control system, C ADMS is the ADMS cost, C area is the area coordination controller cost, C local is the cost of locally automatic controller, and C tx is the construction cost of communication system, which is illustrated in Equation ( 2). (

3) Control Strategy
From the perspective of mathematics, centralized voltage control by optimization is a multi-objective non-linear programming problem [25], which can be described as follows: where u d is discrete control variable, x is dependent variable that can only be controlled indirectly, and u c is continuous control variable.Usually, dependent variable x is node voltage or angle.Discrete control variable u d is switching variable such as position of taps, parallel capacitor and reactor.Continuous control variable u c is the output of DGs.Objective function is network losses minimization in consideration of voltage quality and other economic indexes.It can be illustrated with Interior Point Method: Equality constraints are power flow equations of every node.Inequality constraints include physical constraints of network components and capacity constraints of controllable resources: voltage within the statutory limits; active power and reactive power constraints of controllable resources; and transformer ratio constraint.

Capacity Optimization of DG Considering Voltage Control
The connection of DGs is considered as a solution for environmental pollution, global warming and the rapid depletion of fossil fuels.However, the high penetration of DGs could cause unexpected voltage variation, which is a significant constraint for increasing the share of sustainable sources.By voltage regulation ancillary service with DGs involved, it could avoid the DG units' disconnections due to the infringement of voltage regulatory limits as much as possible [26].Therefore, ADS planning should not only keep cost efficient but also maximize DG generation capacity with voltage acceptable.The control ability, demand on communication system and cost of different voltage control system are different.Therefore, besides capacity and location of DGs, decision variables should also include the category of voltage control system.

Objective Functions
On the one hand, considering the profits of DNOs, economy should be considered in ADS planning.On the other hand, generation capacity of DGs should be as large as possible with voltage acceptable.Therefore, this paper takes minimal comprehensive cost and maximal clean energy generation ratio as objectives.DGs capacity optimization model with active voltage control is solved.

Objective Function 1: Minimizing Comprehensive Cost
Comprehensive cost brought to DNOs should be fully considered, which should be as low as possible.Annual comprehensive cost of DG is comprised of construction cost and operation cost.Construction cost includes connection cost of DG units and construction cost of voltage control systems.The main characteristics of DGs project are its high construction cost in prophase and low maintenance cost in operation.Therefore, power losses cost is the only factor considered in the operation cost [27].Then the following equations illustrate the objective function [28].
where C is annual comprehensive cost; C co is discounted construction cost, C op is operation cost; C i.tr is the connection cost of the ith DG unit or means DGs' newly-built lines cost; N DG is the number of DGs; Tl, T de and T ce are the life of newly built lines, decentralized and centralized voltage control systems, respectively; α stands for bank interest rates; C vo.de is the construction cost of decentralized voltage control system, which is shown in Equation ( 1); C vo.ce is the construction cost of centralized voltage control system, which is shown in Equation (3); x s is a binary variable determining the category of control system; ∆P t loss is the total network losses at time t; T is the total time of dispatch; and λ is residential electricity prices.

Objective Function 2: Maximizing Clean Energy Generation Ratio
If the capacity of DGs is overlarge, it may cause voltage unacceptable and out of control.Therefore, when optimizing the capacity of DGs, the actual amount of DG output should be increased as much as possible on the basis of the existing voltage control system, namely to maximize the clean energy generation ratio.According to the definition of clean energy generation ratio in [29], this objective function can be written as: where ξ is the ratio of clean energy generation, which is the ratio of DGs output to the sum of active load and network losses, and represents the utilization ratio of DGs.T is the total time of dispatch.The dispatching interval is 15 minutes.N G is the number of DGs, P DGn (t) is the output of the nth DG at the time of t, N is the total number of nodes, P Lm (t) is the active load of node m at the time of t, N B is the total number of distribution network branches, and P loss j ptq is the total real power losses on branch j at the time of t.

Constraints Constraint 1: Constraints of Voltage Qualified Rate
Since the output of most DGs is influenced by environment, being intermittent may cause voltage fluctuations.In order to ensure power quality after DG connected to the grid, it is required that voltage qualified rate meet certain requirements after on-line voltage control.
where η t is the voltage qualified rate at time t, δ is the lower limit of voltage qualified rate, N is measured node at time t, N i is the node in feasible region, U i is voltage of node i, and U min and U max are the lower and upper limits of voltage, respectively.

Constraint 2: Constraints of DGs' Annual Comprehensive Cost
Considering the profits of DNOs, DGs' annual comprehensive cost is asked to be acceptable.
where C min and C max are the lower and upper limit of comprehensive cost, respectively.
where n is the number of nodes, P t i , P t iDG , P t Li , Q t i , Q t iDG , Q t Li , and U t i are active power output of the main source, active power output of DG, active power demand, reactive power output of infinite system, reactive power output of DG, reactive power demand, and voltage of node i at time t, respectively.G ij , B ij and θ ij are the conductance, susceptance and voltage phase angle difference between node i and node j.
where ω DG min and ω DG max are the upper and lower limits of DG installed capacity, ω DG i is the installed capacity of the ith DG unit, P t iDG and Q t iDG are active power and reactive power of the ith DG in t period, respectively, δ t iDG is the maximal available DG active power, µ t is the efficiency of DG unit in t period, and Q DGmin and Q DGmax represent the maximum value of reactive power that the converter is able to absorb and/or inject into distribution network, respectively.

Multi-Objective Differential Evolution Algorithm
With two objective functions, DGs capacity optimization is a multi-objective problem.However, to solve the multi-objective optimization problems, it needs to meet two or more objectives.Sometimes the multiple goals might contradict each other.Therefore, in solving multi-objective optimization problems, Pareto solution set is usually used.Here, a new intelligent optimization algorithm called differential evolution algorithm (DE), based on population optimization is used.This algorithm has faster convergence rate, fewer adjustable parameters, simple operation and strong robustness [30].Similar to the standard DE, MODE also includes population initialization, crossover, mutation, selection and other operations.However, being different from DE just based on objective function value, population evolution of MODE is based on the fast non-dominated sorting and the calculation of crowding.
(1) Population Initialization Using the classic uniform random initialization method: where x 0 ij is dimension j of initial individual i. rand is uniformly distributed random number between [0,1].x U j and x L j are the upper and lower limits of variable j, respectively.j = 1, 2, ..., D. D is the dimension of the optimization issues.
(2) Mutation Operation When the difference vector is added to another individual vector selected randomly, the mutated vector is generated.For each target vector x t i , the mutation operation is shown in Equation (13).
where r 1 , r 2 , r 3 P{1, 2, . . ., NP} are different integers, and r 1 , r 2 , r 3 are different from current target vector index i.Thus, the size of population NP ě 4. F is a scaling factor.
(3) Crossover Operation For target vector individual x t i , it will go on crossover operation with the mutated vector u t`1 i , generating trial individuals u t`1 i .To ensure the evolution of individual x t i by random selection, there is at least one bit of u t`1 i that is contributed by v t`1 i .As for other bits, crossover probability factor CR can determine which bit of u t`1 i is contributed by x t i .The crossover equation is shown in Equation (14).
where rand(j) P [0,1] is uniformly distributed random number, randn(i) P [1, 2, . . ., D] is dimension variable index selected randomly, which ensures that at least one dimension variable of trial vector is contributed by the mutated vector.
(4) Selection Operation In accordance with a certain strategy, select individuals from parent generation to second generation.Using a search strategy called "greedy", the test subject u t`1 i is competed with x t i .When the adaptation degree of u t`1 i is better than x t i , u t`1 i can be selected as second generation.Otherwise, x t i is selected directly as second generation.Taking the minimal optimization as an example, the equation for selection operation is as follow: (5) Non-Dominated Ranking The individuals of population are ranked based on non-dominated relationship.A fast non-dominated ranking strategy of typical NSGA-II is used in the case.

(6) Calculation of Congestion Degree
When new population is generated, usually individuals of high level and small aggregation density are reserved to participate in evolution.

(7) Shear Operation
After adding every individual of fronts F 1 -F j to the new offspring in turn, if the number exceeds N, they are ranked according to congestion degree.According to the distance, add individuals selected from F j to the new offspring until the number of individual reaches N.

Optimization Based on Multi-Objective Differential Evolution Algorithm
The flow chart of sizing optimization of DGs considering active voltage control strategy is depicted in Figure 2: 8,193 After adding every individual of fronts F1-Fj to the new offspring in turn, if the number exceeds N, they are ranked according to congestion degree.According to the distance, add individuals selected from Fj to the new offspring until the number of individual reaches N.

Optimization Based on Multi-Objective Differential Evolution Algorithm
The flow chart of sizing optimization of DGs considering active voltage control strategy is depicted in Figure 2:  Specific steps are as follows: ‚ Input the network parameters, initialize the parameters and population, and then conduct the mutation and crossover operations.

‚
According to the data of typical day, determine the day-ahead optimal dispatch schedule.
According to the actual load and output of DG, simulate the real-time voltage regulation with DG participation.

‚
Taking minimal DGs' comprehensive cost and maximal clean energy generation ratio as multi-objective functions and voltage qualified rate as constraint, the optimal compromise solution is solved by intelligent algorithm, and then the most appropriate DG capacity is selected.

Case Study
The simulation is carried on IEEE-33 system, which is a 12.66 kV radial distribution system.The system has one supply point and 32 nodes.The total substation loads of the initial configuration are 3.715 MW and 2.3 Mvar.More parameters of the system can be found in [31].In renewable SHP is the most widely used and flexible for power generation in China, especially in remote areas.According to the situation of actual installation described in [32], four small hydropower plants are connected at different positions along the feeder.The gird topology and installed SHP plants are shown in Figure 3.The parameter values used in MODE are as follows [33]: Population size is 100, the maximum number of iterations is 2000, the maximum and minimum factor ratios are 0.9 and 0.1, respectively, and the maximum and minimum cross ratios are 0.8 and 0.3, respectively.
Sustainability 2016, 8,193 connected at different positions along the feeder.The gird topology and installed SHP plants are shown in Figure 3.The parameter values used in MODE are as follows [33]: Population size is 100, the maximum number of iterations is 2000, the maximum and minimum factor ratios are 0.9 and 0.1, respectively, and the maximum and minimum cross ratios are 0.8 and 0.3, respectively.

Stochastic Modeling of SHP Generation and Load
The generated power of SHP plants depends on the specific weight of the water, the river inflow, the height of the waterfall and the efficiencies of the electric generator [34].However, the river inflow is modeled as a stochastic process, and rainfall has hysteresis and cumulative effects [35].It is difficult to describe the stochastic distribution of SHP generation by specific probability model.Water resource is usually abundant in wet season, which makes the power generation of SHP plants very large.In this case, the system voltage tends to be higher than the statutory limits, which is the focus of our study.Therefore, the output characteristics of SHP in wet season are analyzed here.In wet season, the SHP units usually operate at closely to the rated condition with the maximal efficiency approaching 85% [36].The efficiency of SHP plant is defined as follows: where  t is efficiency of SHP plant during t period,  SHP is the installed capacity of SHP plant, t SHP P is active power of SHP plant during t period.The daily output of SHP plant follows Gamma distribution [35].Its probability distribution is illustrated as follows: where PSHP is active power of SHP plant, (•) is Gamma equation, a = 0.1055, b = −0.0102.

Stochastic Modeling of SHP Generation and Load
The generated power of SHP plants depends on the specific weight of the water, the river inflow, the height of the waterfall and the efficiencies of the electric generator [34].However, the river inflow is modeled as a stochastic process, and rainfall has hysteresis and cumulative effects [35].It is difficult to describe the stochastic distribution of SHP generation by specific probability model.Water resource is usually abundant in wet season, which makes the power generation of SHP plants very large.In this case, the system voltage tends to be higher than the statutory limits, which is the focus of our study.Therefore, the output characteristics of SHP in wet season are analyzed here.In wet season, the SHP units usually operate at closely to the rated condition with the maximal efficiency approaching 85% [36].The efficiency of SHP plant is defined as follows: where µ t is efficiency of SHP plant during t period, ω SHP is the installed capacity of SHP plant, P t SHP is active power of SHP plant during t period.The daily output of SHP plant follows Gamma distribution [35].Its probability distribution is illustrated as follows: where P SHP is active power of SHP plant, Γ(¨) is Gamma equation, a = 0.1055, b = ´0.0102.
As for the stochastic modeling of load, it can be illustrated as follows: where P t load and Q t load are the load at time t, P t 0 and Q t 0 are the basic load at time t which can be obtained by the data of typical day considering the daily variation characteristics of load, and α t and β t are the possible fluctuations of load, which follow Normal distribution [35].Load data on the typical day are shown in Figure 4.In the simulation, parameters of the day-ahead optimal dispatch schedule adopt the parameters of the typical day.To give consideration to both accuracy and the amount of calculation, the typical day is used to represent the operation of the grid throughout the year [37].The initial power factor of every SHP is 1.0.

Simulation Analysis of Multi-Objective Functions Capacity Optimization
In this scenario, SHP plants are mostly located in remote mountainous area.On the one hand, it is very difficult to lay wire communication lines.On the other hand, wireless communication is lack of stability and less effective.The centralized control method needs great additional investment for the construction of a control central and some communication infrastructures.The cost of communication system here is very high.Centralized voltage control scheme is unsuitable to adopt, which needs global information.Thus, centralized voltage control is abandoned and decentralized voltage control is adopted, which needs no communication system.
To make the problem simple, in simulation, minimal network losses and maximal clean energy generation ratio are taken as the multi-objective functions.The Pareto front solution set solved by MODE is shown in Figure 5.In the simulation, parameters of the day-ahead optimal dispatch schedule adopt the parameters of the typical day.To give consideration to both accuracy and the amount of calculation, the typical day is used to represent the operation of the grid throughout the year [37].The initial power factor of every SHP is 1.0.

Simulation Analysis of Multi-Objective Functions Capacity Optimization
In this scenario, SHP plants are mostly located in remote mountainous area.On the one hand, it is very difficult to lay wire communication lines.On the other hand, wireless communication is lack of stability and less effective.The centralized control method needs great additional investment for the construction of a control central and some communication infrastructures.The cost of communication system here is very high.Centralized voltage control scheme is unsuitable to adopt, which needs global information.Thus, centralized voltage control is abandoned and decentralized voltage control is adopted, which needs no communication system.
To make the problem simple, in simulation, minimal network losses and maximal clean energy generation ratio are taken as the multi-objective functions.The Pareto front solution set solved by MODE is shown in Figure 5.
for the construction of a control central and some communication infrastructures.The cost of communication system here is very high.Centralized voltage control scheme is unsuitable to adopt, which needs global information.Thus, centralized voltage control is abandoned and decentralized voltage control is adopted, which needs no communication system.
To make the problem simple, in simulation, minimal network losses and maximal clean energy generation ratio are taken as the multi-objective functions.The Pareto front solution set solved by MODE is shown in Figure 5.All of the solutions presented by the Pareto fronts are valid thresholds for the optimization problem [38].The choice of one of them represents the optimal tradeoff between active losses minimization and renewable energy usage maximization.Considering that the network losses should be within acceptable range, the optimal compromise solution is selected.Results show that the optimal total capacity of SHPs is 7.8261 MW.The capacities of SHPs connected to node 18, node 22, node 25 and node 33 are 0.6682 MW, 2.7836 MW, 2.8936 MW and 1.4807 MW, respectively.Voltage qualified rate at any monitoring time during the day is 1.Network losses of the whole day are 4.526 MWh, and network losses rate is 3.29%.Clean energy generation ratio is 4.451.Power flow runs reversely back in this scenario.The voltage regulation system is started, but not frequently.Power losses during 24 h on typical day are shown in Table 1.To show the reduction of total cost, a comparison of different study cases is conducted when the centralized voltage control approach is used instead.The capacities of SHPs connected to node 18, node 22, node 25 and node 33 are 0.6682 MW, 2.7836 MW, 2.8936 MW and 1.4807 MW, respectively.Results show that, with centralized voltage control, voltage qualified rate at any monitoring time during the day is 1.The network losses of the whole day are 4.0274 MWh and the network losses rate is 2.364%.To analyze the costs clearly, Table 3 shows the result of the optimal scheme.Compared with the centralized method, decentralized voltage control has a great reduction in the total cost, especially in the part of the construction cost of communication system.Although the centralized approach reduces the operation cost, it is too small to ignore compared to the cost of the control system.

Influence of Voltage Control System on DG Capacity
To analyze the influence of voltage control system on DG capacity, the optimal capacity of DGs without voltage control can be calculated in the same way.In the process of simulation, active voltage control is abandoned.Results show that the optimal total capacity of SHPs is 6.397 MW.The capacities of SHPs connected to node 18, node 22, node 25 and node 33 are 0.5397 MW, 2.2384 MW, 2.7493 MW and 0.8696 MW, respectively.Network losses of the whole day are 3.1243 MWh, and network losses rate is 1.46%.Clean energy generation ratio is 3.6424.
From Table 4, with the help of improving the voltage profile, voltage control system increases the capacity of DGs.Although network losses increase either, network losses rate is still acceptable.Meanwhile, the construction cost of centralized voltage control is regarded as infinite compared with decentralized.

The Control Effect of Proposed Voltage Control Method
To prove the ability of the proposed decentralized voltage control method, a further simulation was carried out with daily variations of loads and generations.
From Figures 6 and 7 the proposed voltage control strategy reduces the voltage rise effectively and the voltage profile rapidly drops as a consequence.The effect of proposed voltage control method is illustrated.Once the voltage of the monitoring node exceeds specified limits, the voltage control system begins.It is confirmed that the proposed technique contributes to voltage adjustment.With voltage profile adjusted, larger DGs are permitted to be connected to distribution networks.However, the ability to control voltage is limited with the voltage profile worsen as shown in Figure 8.In this case, the voltage at, or around, 4:00 p.m. is still over the limit with control action, as the algorithm cannot find a valid solution.

Relevance between the Capacity of DG and Load, Circuit Structure
According to the calculation results, distribution of SHP capacity is shown in Figure 9, where SHP 1, SHP 2, SHP 3 and SHP 4 stand for SHP plants connected to node 18, node 22, node 25 and node 33, respectively.

Relevance between the Capacity of DG and Load, Circuit Structure
According to the calculation results, distribution of SHP capacity is shown in Figure 9, where SHP 1, SHP 2, SHP 3 and SHP 4 stand for SHP plants connected to node 18, node 22, node 25 and node 33, respectively.From Figure 8, capacities of SHP 2 and SHP 3 are larger.The reason is that they are at the load center and closer to the system bus with less influence on voltage.As SHP 1 and SHP 4 are far away from the system bus, the load is lighter and the load relevance is weaker, easily making network voltage out of range.Compared to SHP 3, the range of capacity distribution and the maximal achievable capacity of SHP 2 are larger.

Conclusions
This study has presented a new ADS planning approach for promoting renewable energy usage and the benefits of DNOs while keeping voltage profile acceptable.Especially, the costs of specific control system and its supporting systems are included in the optimization model, and the effect of voltage control is included in the simulation.Based on the theoretical analysis and simulation results above, this paper draws the following conclusions: (1) The cost of control system and its supporting systems are included in the objectives of planning.Therefore, it takes the impacts of control system on the costs of operation and construction into account.The control effect of voltage regulation is included in the simulation.Therefore, the control ability can be verified instead of being estimated roughly.Having made improvements in the above two aspects, the precision of planning can be improved.(2) Different control systems have different influences on the planning and operation of the grid.Both the decentralized and centralized approaches can reduce voltage rises and increase the acceptable capacities of DG units to a certain degree, and the effect of the latter is better.However, the centralized approach means a great investment in related costs.With the power grid becoming smarter, more automatic and complicated, the cost of these systems will account for a large share of the total cost.If the cost saving is the priority, adopting the decentralized approach is suggested.If the control effect or DG penetration is the priority, adopting the centralized approach is suggested.From Figure 8, capacities of SHP 2 and SHP 3 are larger.The reason is that they are at the load center and closer to the system bus with less influence on voltage.As SHP 1 and SHP 4 are far away from the system bus, the load is lighter and the load relevance is weaker, easily making network voltage out of range.Compared to SHP 3, the range of capacity distribution and the maximal achievable capacity of SHP 2 are larger.

Conclusions
This study has presented a new ADS planning approach for promoting renewable energy usage and the benefits of DNOs while keeping voltage profile acceptable.Especially, the costs of specific control system and its supporting systems are included in the optimization model, and the effect of voltage control is included in the simulation.Based on the theoretical analysis and simulation results above, this paper draws the following conclusions: (1) The cost of control system and its supporting systems are included in the objectives of planning.
Therefore, it takes the impacts of control system on the costs of operation and construction into account.The control effect of voltage regulation is included in the simulation.Therefore, the control ability can be verified instead of being estimated roughly.Having made improvements in the above two aspects, the precision of planning can be improved.(2) Different control systems have different influences on the planning and operation of the grid.Both the decentralized and centralized approaches can reduce voltage rises and increase the acceptable capacities of DG units to a certain degree, and the effect of the latter is better.However, the centralized approach means a great investment in related costs.With the power grid becoming smarter, more automatic and complicated, the cost of these systems will account for a large share of the total cost.If the cost saving is the priority, adopting the decentralized approach is suggested.If the control effect or DG penetration is the priority, adopting the centralized approach is suggested.(3) The proposed approach allows DNOs to obtain benefits by inducing the comprehensive cost and maximizes the usage of renewable energy.The algorithm of MODE can compute the optimal capacity of DG units.
(4) In the absence of a widespread communication channel, decentralized voltage control method provides an effective solution to mitigate voltage problem.The simulation results show that the proposed voltage control method helps improve voltage to some extent, and DG capacity can be increased by 12.88%.(5) Compared with the traditional voltage control methods such as the installation of additional reactive power supply, the proposed voltage control from DGs strategy has more potential.
Traditionally, it is difficult to determine the optimal location of reactive power controllers because the configuration of the distribution system may be changed in the future.Furthermore, the setting costs for the installation of additional reactive power compensator is not beneficial for power utilities.The case study proves the effectiveness and advantages of the proposed method.(6) The optimal capacity of DG near the system bus is relatively larger.The optimal capacity of DG near heavy loads and with better load relevance is also relatively larger.The basic load demand of each node in IEEE-33 system used in the case study is shown in Table A3.Pd means the active load of each node.Qd means the reactive load of each node.In Section 4.4, the specific voltage of node 18 during the whole day is shown as follows.The scenario that the proposed voltage control strategy reduces the voltage rise effectively is shown in Table A4.The scenario that the proposed voltage control strategy cannot reduce the voltage rise effectively is shown in Table A5.V1 means the voltage before control.V2 means the voltage after control.

Figure 1 .
Figure 1.Simulation flow chart of decentralized voltage control system.

Figure 1 .
Figure 1.Simulation flow chart of decentralized voltage control system.

Figure 2 .
Figure 2. Flow chart of sizing optimization of distributed generation (DG) considering voltage control system.Specific steps are as follows:  Input the network parameters, initialize the parameters and population, and then conduct the mutation and crossover operations. According to the data of typical day, determine the day-ahead optimal dispatch schedule.According to the actual load and output of DG, simulate the real-time voltage regulation with DG participation.

Figure 2 .
Figure 2. Flow chart of sizing optimization of distributed generation (DG) considering voltage control system.

Figure 7 .
Figure 7. Daily voltage qualified rate with (solid line) and without (dotted line) the control action.

Figure 8 .
Figure 8. Daily voltage profile# with (solid line) and without (dotted line) the control action.

Figure 6 .Figure 6 .
Figure 6.Daily voltage profile with (solid line) and without (dotted line) the control action.

Figure 7 .
Figure 7. Daily voltage qualified rate with (solid line) and without (dotted line) the control action.

Figure 8 .
Figure 8. Daily voltage profile# with (solid line) and without (dotted line) the control action.

Figure 6 .
Figure 6.Daily voltage profile with (solid line) and without (dotted line) the control action.

Figure 7 .
Figure 7. Daily voltage qualified rate with (solid line) and without (dotted line) the control action.

Figure 8 .
Figure 8. Daily voltage profile# with (solid line) and without (dotted line) the control action.

Table 1 .
[39]r losses profile based on optimal capacity considering voltage control.Average electricity sales price in a certain region is 0.48 RMB/kWh.To accept the newly-built four SHP plants, electricity lines use LGJ-185.The total length of lines is 30 km.The unit price is RMB 130 million per km.The Comprehensive unit prices of communication equipment in distribution network are shown in Table2[39].

Table 2 .
Comprehensive unit prices (k RMB) of communication equipment in distribution network Note: The unit k RMB means thousand RMB.

Table 3 .
Cost analysis of the optimal scheme with different control systems.

Table 4 .
Influence of voltage control system on DG capacity.
Daily voltage qualified rate with (solid line) and without (dotted line) the control action.

Table A3 .
Load demand of IEEE 33 bus system.

Table A4 .
Voltage profile with effective regulation.

Table A5 .
Voltage profile without effective regulation.