Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Optimal Placement of D-STATCOMs into the Radial Distribution Networks in the Presence of Distributed Generations

Babak Safari Chabok, Ahmad Ashouri*

Department of Electrical Engineering, Khodabandeh Branch, Islamic Azad University, Khodabandeh, Iran

Abstract: This paper proposes a combination of Discrete Imperialistic Competition and Nelder-Mead algorithms to solve D-STATCOMs placement optimization problem. It is assumed that optimal number, locations and sizes of D-STATCOMs are determined in radial distribution network while Distributed Generations are previously installed in it. Indeed, we focus on voltage control and concept of reactive power management in annual scheduling time interval. Objective function is defined in terms of network active power losses and voltage stability characteristic at different load levels by appropriate weighting factors. Voltage Stability Index (VSI) is applied to identify the weak buses in radial distribution network. To validate performance and effectiveness of proposed DICA-NM hybrid algorithm, simulation studies are applied on the 30-bus IEEE radial distribution test feeder. Finally, numerical results for some important network variables have been compared in three different case studies.

Keywords: Discrete Imperialistic Competition Algorithm and Nelder-Mead (DICA-NM), Distributed Generations (DGs), location and Size of D-STATCOM, Forward and Backward Power Flow, Voltage Stability Index (VSI)

Citation: Babak Safari Chabok, Ahmad Ashouri, Optimal Placement of D-STATCOMs into the Radial Distribution Networks in the Presence of Distributed Generations. American Journal of Electrical and Electronic Engineering, 4 (2). 40-48. January 2016. doi:10.12691/ajeee-4-2-1

1. Introduction

In recent years, restructuring process in power system, demand side competition promotion, in addition to generation side and deregulation into the power system management causes distribution networks to play an important role in power system planning. For above challenging situation, the important subjects in optimal planning and operation planning of electrical distribution system are loss reduction, voltage and current profile improvement which are related to voltage control and reactive power management concept. This concept includes all planning and actions to improve voltage profile and voltage stability. It consists of two steps: first step is reactive power resources planning and second step is reactive power dispatch on reactive power resources that have up to now been installed on networks. Reactors, capacitor banks and transformer equipped to Under Load Tap Changer (ULTC) are available reactive power resources in distribution network. Main goal of voltage control and reactive power management are indeed correct coordination and setting between available controlling devices in sub-transmission substation and distribution feeders. Nowadays, DG resources, because of high importance in energy generation process, are vastly applied in the electrical distribution systems. Synchronous generators are one of the most applicable types of DGs that are installed in the medium voltage systems. DGs can affect the voltage control devices and reactive power especially on Under Load Tap Changer because of its operation ability in two different modes including power factor control (PQ) and voltage control (PV) and also installation ability on different locations in distribution network. Therefore, for reliable voltage regulation in system, DG unit connections should be coordinated with other system devices. During reactive power management studies with high penetration of DGs in networks, high flexibility of DGs reactive power related to network should be considered. In spite of these, power system operators cannot use all available reactive abilities of DGs up to now due to lack of optimal management, thus distribution systems with large capacity of DGs are faced with major issues such as voltage stability and voltage collapse. To cover this gap, recently, D-FACTS devices due to control ability of voltage profile and power flow are proposed as an applicable option for voltage stability improvement. By development in modern control techniques and power electronic devices technology, use of these devices would be inevitable to improve the quality of distribution network operation. Desirable controlling methods for improvement in distribution network operation have been achieved by Custom Power Devices (CPD) application. These elements are used in distribution networks for voltage profile and power quality improvement, and loss reduction, which consist of devices such as Solid State Transfer Switch (SSTS), Distribution Static Compensator (D-STATCOM) and Unified Power Flow Controller (UPFC).D-STATCOM as a Voltage Source Converter (VSC) with shunt connection is more applied for power quality concerns. Under voltage sag conditions or lines overload, bus voltage to which D-STATCOM is connected can be regulated with proper compensation current injection. Injected current to voltage balance/regulation in desired level is estimated by bus voltage and requested reactive power. Optimal use of these devices is achieved via finding out its optimal numbers, locations, and sizing. Ref. [1] pays attention to the optimal placement of D-STATCOM in electrical distribution systems by use of particle swarm optimization (PSO) method. Ref. [2] proposes Immune Algorithm (IA) for finding the optimal location of unified power flow controller (UPFC) for the optimal power flow determination and congestion management. Ref. [3] presents a method for optimal placement of D-STATCOM using the firefly algorithm (FA) by goal of power quality enhancement in distribution network and obtained results are also compared with the particle swarm optimization and genetic algorithm. In ref. [4], optimal location and sizing of DSTATCOM for the sake of power loss reduction, and improvement of current and voltage profile in distribution networks are investigated using Immune Algorithm while considering voltage and current profile improvement, minimum cost of installation of DSTATCOM and maximum power loss reduction are integrated into an objective function through appropriate weighting factors. This paper proposes an application of DICA-NMhybrid algorithm to determine the optimal numbers, locations and sizing of D-STATCOMs in radial distribution networks in presence of Distributed Generation (DGs). Objective function is defined in two terms by appropriate weighting factors: summation of active power losses, D-STATCOMs devices' installation costs after normalization and reverse value of normalized voltage stability index (VSI) for system critical bus at three different load levels. Problem constraints consist of under/over current and voltage limits, minimum and maximum active and reactive power generation of DGs, number and size of D-STATCOMs, Voltage Stability Index (VSI).Rest of the paper is organized as follows: In section.2, Forward-Backward power flow is reviewed as an efficient study method in radial distribution network and then entering procedure of DGs and D-STATCOM into the Forward-Backward distribution load flow is presented. In section.3, Voltage Stability Index (VSI) formulation is performed. In section.4, Optimal Placement of D-STATCOMs as an optimization problem is modeled by an objective function and system constraints. In section.5, proposed DICA-NM hybrid algorithm is illustrated. Simulation studies and numerical analysis have been performed in section.6 and concluding remarks appear in the last part of paper.

2. Forward-Backward Power Flow

Because of high ratio and vast changing range of in distribution systems, especially, in radial type, common power flow methods hardly converge. Therefore in this paper forward-backward load flow is used for power flow study of radial distribution systems [5]. It is assumed that, substation bus of radial feeder is slack bus with known voltage magnitude and angle, under initial conditions that other buses' voltages are equal to slack bus, and power losses of all branches are equal to zero. Forward-backward power flow repetitive algorithm is composed of four stages: these four stages in repetition are as below:

A) Nodes Power Calculation

Considering this point that the absorbed power with constant current and impedance loads and also injected reactive power by capacitors due to change in nodes voltages will be changed during algorithm repetition, all aforementioned values should be calculated during performing repetition of algorithm, again. Only, absorbed power with constant power loads will be unchanged during sequential repetitions.

B) Backward Sweep for Summing Over All Branches

With starting from end-bus and move toward the slack bus, branches' powers are calculated according to relation (1):

[picture] | (1)

Indexes and are sending and receiving nodes of branch , respectively; is a set of branches which from one side are connected to the receiving node of branch and from the other side are connected to sending node of branch . Index is a set of sending end nodes of branches in set . Index is absorbed power by the load which is connected to receiving node of branch n. Index is power flow from node toward node through branch and is power losses in branch n which is connected to nodes and and alsoit is assumed in first iteration equal to zero.

C) Forward Sweep to Update Nodes Voltages

With starting from the branch which is connected to slack bus and move toward the end branches, currents in sending bus of branch , voltages in receiving bus of branch and active and reactive losses are calculated by relations (2), (3) and (4):

[picture] | (2)

[picture] | (3)

[picture] | (4)

Index is the current flow from node toward node through branch , and is series impedance of branch between nodes and . Indexes and are nodes voltages and of branch and is power losses of branch between nodes and .

D) Voltage Miss match

After performing previous stages in each repetition, voltage miss match is calculated for all buses by relation (5):

[picture] | (5)

Even if one of miss match values is greater than convergence criteria, all stages of algorithm should be repeated until convergence condition is achieved.

2.1. DGs Model in Forward-Backward Power Flow

One of the most common DG types is synchronous generators which are installed in the medium voltage of distribution system and it has the capability of direct connection to network. DGs, depending on operation of excitation system, can be modeled as voltage control (PV) and power factor control (PQ) [6,7]. Aforementioned types of DGs used in this paper have been modeled as power factor control (PQ). Active and reactive power generation limits in synchronous generators are determined with regard to operational curve, which considers mechanical limit of turbine, thermal limit of armature winding, etc. for generator at working set point. DG units that are controlled in form of PQ buses, as a negative load, enter the forward-backward load flow model, but DG units that are controlled in form of PV buses need additional process in power flow solution [8].

2.2. D-STATCOM Model in Forward-Backward Power Flow

A D-STATCOM can work as synchronize voltage source with variable amplitude and phase angle. Therefore, it has a capability to control its bus voltage and correct power factor. In steady state operation with heavy load or some types of short circuit, usually, D-STATCOM injects appropriate current to Point of Common Coupling (PCC) and regulates bus voltage to nominal rating or predetermined value. Generally, D-STATCOM can transfer active and reactive powers to network, simultaneously [9]. As can be seen, Figure 1 shows D-STATCOM installation in bus of radial distribution feeder. In this paper, we consider D-STATCOM application only for transferring reactive power. Many of distribution systems have a radial structure that is fed through one side. It is assumed that three phase radial distribution feeder has a balance condition.

Figure 1. D-STATCOM installation in bus of radial distribution feeder

In this investigation, D-STATCOM device is used for compensating voltage of bus of radial distribution feeder to desirable value as can be seen in Figure 1; Phasordiagram is available in Figure 2. Considering the system voltage, consequently, should be kept in quarter.

Figure 2. Phasor diagram of voltage and current of system in Figure 1 after D-STATCOM installation

According to the Figure 1 with D-STATCOM installation in bus , currents and flow in branch, simultaneously and based on Phasor diagram in Figure 2, we have:

[picture] | (6)

[picture] | (7)

By separating real and imaginary parts of equation 7 and some calculations, we have:

[picture] | (8)

[picture] | (9)

As can be found from equation block9, there are two roots for variable x and therefore two values for and are calculated. To determine the correct solution, boundary conditions in these roots have been tested as (10).

[picture] | (10)

The results show that the correct answer to equation 7 for and are calculated with appropriate separated real and imaginary equations as below (11) and (12):

[picture] | (11)

[picture] | (12)

Therefore, current and voltage Phasors of D-STATCOM in the point of common coupling with network can be calculated by relations (10), (11) and (12). Finally, injected reactive power to network by D-STATCOM to recover bus voltage to desire value can be express as below relation (13) which “*” is conjugation mark.

[picture] | (13)

To evaluate network power flow considering D-STATCOM in distribution system for any iteration of forward sweep, the magnitude of compensated bus voltage can be assumed equal to desirable value , phase angle of compensated bus, reactive power and size of D-STATCOM can currently be calculated by above equations. Voltages are up to dated and reactive power injection by D-STATCOM is applied to continue the forward sweep and with goal of determining load current in the next backward sweep. This process is repeated until power flow can converge [10].

3. Voltage Stability Index (VSI)

Based on voltage level, electrical power systems can be divided in generation, transmission and distribution systems. Modern distribution networks are faced with load demand increment, and also, experience significant changes in load demand from low to high levels during a day. In different industrial areas of distribution network, voltage collapse has been observed under critical load demand. In ref. [11], real occurrence of these phenomena have been reported which system voltage periodically collapsed and to avoid the repetition of voltage collapse, fast reactive compensation has been required. Up to now, a few investigations have been carried out about voltage stability in radial distribution networks, unlike transmission systems. In this paper, voltage stability index (VSI) proposed in ref. [12] is applied to identify weak buses of radial distribution networks due to its simplicity and enough accuracy. The node in which VSI value is minimum, related to other nodes, is more susceptible to voltage collapse. Consider Figure 3 and Figure 4 for mathematical formulation of voltage stability index. After forward-backward power flow study, voltage and current of all branches and nodes are known, therefore the values of and for node can be easily calculated and then will be determined.

Figure 3. single line diagram of 15-bus radial distribution feeder

Figure 4. equivalent circuit of radial distribution feeder for Figure 3

From ref. [6] voltage stability index can be formulated as below relation (14):

[picture] | (14)

Indexes and are total active and reactive powers fed to load and branches which are connected to receiving node of branch . For steady state operation of radial distribution networks, should be greater than zero and lower than one, for all nodes of the network. By use of voltage stability index calculation, stability level of radial distribution network can be measured and consequently if VSI is weak, then, proper operation should be performed.

4. Problem Formulation

By development in management system and distribution automation, the load demands among feeders and secondary bus of sub-transmission substation can be estimated for a day to one year via short-term and long-term forecasting methods. Therefore, by determining annual load profile which consists of load levels and its durations, D-STATCOM placement optimization problem in radial distribution network in presence of DGs can be formulated and solved by the goal of loss reduction and voltage stability improvement.

4.1. Objective Function of D-STATCOMs Placement in Radial Distribution Feeders

In this paper, optimal number, location and size of DSTATCOMs are determined, considering DGs presence in radial distribution networks. Objective Function (OF) is minimizing the sum of two important terms with appropriate weighting factor. First one is total normalization value of active power losses cost and D-STATCOMs installation and operation costs. Second term is reverse value of normalized critical bus voltage stability index in three low, medium and peak load-levels. Indeed, objective function is implicitly a cost-benefit analysis between total saving in active losses cost against D-STATCOMs installation and operation costs and simultaneously improve voltage stability characteristic of network buses whose mathematical formulation has been expressed as relation (15).

[picture] | (15)

Weighting factor W1 and W2 can be selected based on priorities of decrease in network costs and voltage stability characteristic of radial distribution system. Index is defined as relation (16).

[picture] | (16)

In above relation can be calculated by summation of distribution network losses cost without D-STATCOM presence and annual investment cost of D-STATCOM including installation and operation with regard to three low, medium and peak load levels. So, index is defined as follow relation (17). Index shows the cost of energy losses and is time duration of load level and is the ratio of time duration of load level to total time period [13]. Index is active power losses of distribution network in load level after D-STATCOM installation and operation.

[picture] | (17)

Investment cost of D-STATCOM can be extracted by mathematical formulation as relation (18).

[picture] | (18)

Index shows the annual cost of D-STATCOM in load level, is investment cost assigned to load level, is the asset rate of return and is longevity of D-STATCOM [14]. Total cost saving (TCS) is the difference between total energy losses cost before installation and total energy losses cost and annual cost of D-STATCOM after installation in different load levels which is calculated by relation (19).

[picture] | (19)

Index in OF is defined as relation (20) by dividing .

[picture] | (20)

Indexes and show respectively the weakest and the strongest network buses as voltage stability characteristic in three load levels by performing forward-backward power flow study during solution of D-STATCOM placement optimization problem.

4.2. Technical Constraints in Radial Distribution Feeders

Usually, voltage control and reactive power management in long-term scheduling have been performed via optimal placement of D-STATCOM and other reactive resources in radial distribution network to maintain network in static safe operation limits and improve voltage stability profile. The limits that should be satisfied for the optimal placement of D-STATCOM in presence of DGs in radial distribution system are as follows:

4.2.1. Permissible Range of Bus Voltage and line Current Magnitude

In steady state conditions, buses voltages and lines currents magnitude of radial distribution feeder should be in permissible range.

[picture] | (21)

[picture] | (22)

4.2.2. Reactive Power Injection Limit of D-STATCOMs

In steady state conditions, D-STATCOM reactive power injection should be in allowable range to the bus to which it is connected.

[picture] | (23)

4.2.3. Active and Reactive Power Generation Limit of DGs

In steady state conditions, active and reactive power generation constraints of DGs are modeled as relations (24), (25) and (26).

[picture] | (24)

[picture] | (25)

[picture] | (26)

Indexes and are active and reactive power generation of DG units, respectively. Indexes are minimum and maximum active power generation of DGs units, respectively. Index and are minimum and maximum reactive power generation of DGs respectively. Index is maximum apparent power of DGs.

5. Proposed DICA-NM Hybrid Algorithm

Imperialistic Competition Algorithm (ICA) is a new method in evolutionary computations that helps to find optimal solution to different optimization problems. This algorithm presents a solution to different problems by mathematical model of political-social human evolution process. From application viewpoint, this algorithm is categorized in optimization evolutionary algorithms like Genetic Algorithm (GA), Particle Swarm Optimization (PSO), Ant Colony Optimization (ACO), Simulated Annealing (SA), which lead to global optimal solution. Similar to other algorithms in this category, ICA provides a primary set of probabilistic solutions. The assimilation policy, imperialistic competition and revolution construct the basis of this algorithm. This algorithm presents the operators in a regular form as an algorithm which can help to solve complex optimization problems by imitation of social evolution, economic and political processes in countries and mathematical modeling of parts of these processes. Indeed, Algorithm sees the optimization problem' solutions in form of a country and, tries to gradually improve these solutions during repetitive process and finally get to optimal solution. Similar to other evolutionary algorithms, this algorithm also begins with a number of primary random populations which each of them has been called a country. Some of the best population elements (equal to genes in genetic algorithm) are selected as imperialist. The rest of the populations are called colonies. Imperialists, depending on their power draw these colonies to their side with a special procedure that is explained later. Total power of each empire depends on both imperialist country (as central core) and its colonies. In mathematical form, this dependency is modeled by imperial power definition as summation of total power of imperial country and a percentage of its colonies' average power. After forming primary empires, their imperialistic competitions begin between them. Each empire which cannot be success in imperial competition and increase itself power (or at least maintains its power), will be removed from imperialistic competition scene. Therefore, an imperial survival would be in absorbing and seizing the territory of adversary empires. Thus, during imperialistic competitions, gradually bigger empires gain more power and the weaker ones will be eliminated. Empires to increase their power will have to expand their territories [15]. Imperialistic competition repetitively continues and in different populations until just one empire remains that is the global optimal solution of problem [16,17]. NelderMeadalgorithm is one of the best known algorithms and also a very fast method for finding local optima in multi-dimensional unconstrained optimization problems. Unlike gradient method, this method does not use function derivatives for finding optimal solution. Nelder mead algorithm converges to local optima by constructing a unique structure. By use of this unique structure, search in different directions with high potential has been performed to optimize the objective function. Unique structure in NelderMead method is defined as a geometrical figure which is composed of pointsthat are considered as the vertices of a working simplex which is the number of variables of objective function in optimization problem. In each iteration NM algorithm starts calculating the reflection of worst point (vertex with maximum value in objective function) along the central point (vertex with mid value in objective function).Based on the calculated value in first step, NM algorithm performs expansion and contraction actions in new structure. In other words, values of objective function in each iteration has been evaluated for all vertices of created structure and highest value of objective function for each vertex of structure in each iteration replaces the worst value calculated from previous step and otherwise structure is contracted in along the best point (vertex with lowest value of objective function).This process will be repeated until desired error is achieved. Convergence speed is affected by tree parameters: α, β and γ which are coefficients of reflection, contraction and expansion, respectively [18]. It should be noted that ICA-NM hybrid algorithm operation is based on real numbers. With regard to structure of problem variables, solution is categorized in two parts: in first part with regard to maximum permissible number of D-STATCOMs, the buses numbers which D-STATCOMs should be installed on determined locations and in second part the size of D-STATCOM and DGs generation in three load levels are determined. Figure 5 shows flowchart of DICA-NM hybrid algorithm. Because buses’ numbers as problem variables of D-STATCOM placement in first part are discrete types, here ICA algorithm cannot be directly used in problem solution and therefore DICA-NMHybrid Algorithm is applied for aforementioned problem solution. Under these conditions, DICA-NM search space for problem variables in first part will be according to relation (27) for and ICA-NM search space in second part will be as relation (28) for .

[picture] | (27)

[picture] | (28)

That is used to get to the integer partand also should be limited as follows (29):

[picture] | (29)

Figure 5. flowchart of DICA-NM hybrid algorithm for D-STATCOM placement

Indexes and are the number of network buses and maximum permissible number of D-STATCOMs for installation, respectively.

6. Simulation Studies

In this section to validate the effectiveness and performance of proposed DICA-NM hybrid algorithm, optimization problem of D-STATCOM placement in presence of DGs has been solved for 30-IEEE radial test feeder [19]. As can be seen in Figure 6, there is a sub-transmission substation 63kv/20kv, therefore voltage level of radial feeder is 20kv, and apparent power of substation transformer is 500KVA. It is assumed that optimal location of Distributed Generation (DG) was determined on bus8 based on previous network studies and its nominal rating is 4.5MVA with the power factor of 0.9.

Figure 6. Single line diagram of 30-bus IEEE radial test feeder

Branches' maximum current limit is 520A. Allowable range for bus voltage magnitude is considered from 0.9 to 1.10p.u. It is assumed maximum number of D-STATCOMs which can be installed on 30-bus radial feeder is equal to 3. In addition, maximum reactive capacity of D-STATCOM is assumed 5000kvar. Load demands along the feeder are modeled as constant power type which can cause the worst case of voltage stability for radial distribution network. Table 1 shows investment parameters of D-STATCOMs and values of three load levels including light, medium and peak and their durations in a year. In addition, network load demand on each bus is available in ref. [20]. Proposed approach has been implemented on computer ASUS (Intel) with CPU processor 2.4GHz.Three different case studies have been considered on IEEE 30-bus radial test feeder as follow:

1) Radial distribution test feeder without DGs and D-STATCOM presence.

2) Radial distribution test feeder without DGs and with D-STATCOM presence.

3) Radial distribution test feeder with DGs and D-STATCOM presence.

Table 1. D-STATCOMs investment parameters in IEEE 30-bus radial test feeder

Achievement results and conclusion remarks related to bellow subjects for three different case studies are validated and compared: Voltage profile in each load-level, network losses costs of active power, D-STATCOM installation and operation costs in network, reactive power injection by D-STATCOM, active and reactive power injections by DGs, economic benefit of D-STATCOM installation and coordination between DG and D-STATCOMs belong to the radial feeder due to network losses reduction, and improvement in voltage stability characteristic of radial distribution feeder.

6.1. Case 1: Base case (Radial Distribution Feeder without D-STATCOM and DGs)

In base case study, there are no DGs and D-STATCOMs in radial test feeder. Branches power flows, active and reactive power losses and also buses voltages are determined by use of forward-backward power flow method. Total active power loss in each load-level, and critical bus which is susceptible to voltage instability have been shown in Table 2. Total active power losses cost of distribution network is equal to in a year. In addition, simulation results show that Bus14 has the lowest VSI values among all network buses in three load levels and therefore, it is more susceptible to voltage collapse than other buses along the radial distribution feeder.

Table 2. Simulation results for base case study

6.2. Case2: D-STATCOM Placement in Radial Distribution Feeder without DGs

In second case study, it is assumed that all information as compared to base case study remains unchanged. In new case study optimization problem of D-STATCOM placement without DGs in radial distribution test feeder is solved again by DICA-NM hybridalgorithm. From simulation results in case2, Bus6 is selected as optimal location of a D-STATCOM installation in the 30-bus radial distribution feeder. It is expected that optimal amount of reactive power is injected to network by D-STATCOM and thus, as compared to case1 not only network voltage and current profiles are improved but also total active power losses cost in radial distribution network is decreased. According to our expectation, in case2 the total active power losses cost is equal to 103040.89$ which shows 6879$ decrement in annual cost as compared to base case study. From simulation result in case2, Total Cost Saving (TCS) in a year is equal to4757.1$ due to only a D-STATCOM presence in radial distribution feeder. Network active power losses, D-STATCOM reactive power injection, D-STATCOM current magnitude and its phase angle and also annual cost of D-STATCOM in load level have been shown in Table 3.

Table 3. Simulation results for Case2

From simulation result in case2, the lowest voltage stability index is calculated in three load levels for Bus14 similar to base case, but, VSI values show little increment as compared to base case study.

6.3. Case 3: D-STATCOM Placement in Radial Distribution Feeder with DGs

In third case study, D-STATCOM placement optimization problem is solved again in presence of DGs. In case3 similar to case2 all network information remains unchanged as compared to case1.It is assumed that one DG is installed on Bus8 in radial distribution test feeder. Now, optimization problem of D-STATCOM placement in presence of DGs can be solved again with proposed DICA-NM hybrid algorithm. It can be seen from simulation results that Bus6 is selected again as optimal location of D-STATCOM installation and its optimal number is one similar to previous case study (case2).Under new condition total active power losses cost of distribution network is equal to , and we can see and decrements as compared to case 1 and case 2, respectively, due to reactive power management coordination between D-STATCOM and DGs. From simulation result in case3, Total Cost Saving (TCS) in a year is equal to , due to simultaneous presence of DG and D-STATCOM in radial distribution feeder. In case 3, we can calculate the values of network active power losses, D-STATCOM reactive power injection, D-STATCOM current magnitude and its phase angle, annual cost of D-STATCOM and also DGs active and reactive power injections in three load-level as the Table 4.

Table 4. Simulation results for Case3

Table 5. The comparison of all buses voltage magnitude in test feeder

Results show that total network active power losses in three load-level is decreased as compared to previous case studies using coordinated reactive power management of DGs and D-STATCOM. In three case studies values of all bus voltage magnitude are in desirable ranges between and while their comparison comes into the Table 5. In case3 again Bus14 is identified with the lowest VSI values in three load-level as critical network bus which shows that this bus is more susceptible to voltage instability relative to other buses of radial distribution test feeder. Change of voltage stability index (VSI) for critical network bus (Bus14) in each network load-level conditions including light, medium and peak have been shown in Figure 7 and for three case studies are compared. This comparison shows that VSI index for case 2 and case 3 are significantly affected in each three load-level, and these effects on VSI values are ascending especially in peak loading condition of distribution network.

Figure 7. voltage stability index comparison between three case studies

7. Conclusion

In this paper a new DICA-NM hybrid algorithm is proposed for solving optimization problem of the annual D-STATCOM placement when distributed generations have been previously installed on distribution network. Objective function is minimizing the sum of normalized active power losses, D-STATCOM installation cost and the reverse value of normalized voltage stability index. Forward-Backward method is applied for power flow study in radial distribution network. Simulation analyses are applied for three case studies. Numerical results show that utilizing coordinated D-STATCOMs and DGs as voltage control and reactive power management resources cause more reduction in active power losses and more improvement in voltage profile of network buses than conditions that D-STATCOMs have been applied without DGs presence in distribution feeder. Therefore, it can be concluded that new DICA-NM hybrid algorithm that converges to optimal solution of D-STATCOM placement problem with appropriate accuracy would be comparable to other evolutionary methods such as GA, PSO, ACO and etc.

References

[1]   | Y. Del Valle, J.C. Hernandez, G.K. Venayagamoorthy, and R.G. Harley,“ Multiple STATCOM allocation and sizing using particle swarm optimization,” In: IEEE conference power systems and exposition, pp. 1884-1891, 2006.
[2]   | S.A. Taher, M.K. Amooshahi ,“Optimal placement of UPFC in power systems using immune algorithm,” Simulation Modeling and Practice Theory, Vol. 19, No. 5, pp. 1399-412, 2011.
[3]   | M. Farhoodnea, A. Mohamed, H. Shareef, H. Zayandehroodi, “Optimum D-STATCOM placement using firefly algorithm for power quality enhancement,” IEEE 7th International Power Engineering and Optimization Conference (PEOCO), pp. 98-102, June 2013.
[4]   | S.A. Taher, and S.A. Afsari,“Optimal location and sizing of D-STATCOM in distribution systems by immune algorithm,” International Journal of Electrical Power&Energy Systems, Vol. 60, pp. 34-44, Sep. 2014.
[5]   | H. Hamouda, K. Zehar, “ Efficient Load Flow Method for Radial Distribution Feeders”, Journal of Applied Sciences, Vol. 6, No. 13, pp. 2741-2748, 2006.
[6]   | A. Losi and M. Russo, “Dispersed Generation Modeling for Object-Oriented Distribution Load Flow,” IEEE Transactions on Power Delivery, vol. 20, no. 2, pp. 1532-1540, 2005.
[7]   | H. Chen, J. Chen, D. Shi, and X. Duan, “Power Flow Study and Voltage Stability Analysis for Distribution Systems with Distributed Generation,” in Power Engineering Society General Meeting, june2006, IEEE (PP. 8-PP).
[8]   | S.M. Moghaddas-Tafreshi and E. Mashhour, “Distributed Generation Modeling for Power Flow Studies and a Tree-Phase Unbalanced Power Flow Solution for Distribution Systems Considering Distributed Generation,” Electric Power System Research, vol. 79, pp. 680-686, 2009.
[9]   | K. Somsai, T. Kulworawanichpong, “Modeling, simulation and control of DSTATCOM using ATP/EMTP,” In: 13th IEEE conference harmonics and quality of power, pp. 1-4, 2008.
[10]   | M. Hosseini, H.A. Shayanfar, M. Fotuhi, “Modeling of series and shunt distribution facts devices in distribution load flow,” Journal Electrical System, pp. 1-12, 2008.
[11]   | G. Brownell, H. Clarke, “Analysis and solutions for bulk system voltage instability,” IEEE Computer Applications in Power Systems, Vol. 2, No. 3, pp. 31-35, 1989.
[12]   | M. Chakravorty, D. Das,“Voltage stability analysis of radial distribution networks,” Electrical Power and Energy Systems, Vol. 23, pp.129-135, Elsevier Journal (2001).
[13]   | L. Dasilva, S. Carneiro, E. DeOliveira, E. DeSouzaCosta, J. RezendePereira, P. Garcia,“ A heuristic constructive algorithm for capacitor placement on distribution systems,” IEEE Transaction on Power System, vol. 23, no. 4, pp.1619-1626, Nov. 2008.
[14]   | E. Lakervi, E.J. Holmes, “Electricity distribution network design,” IET, No. 21, 1996.
[15]   | E. Atashpaz-Gargari, C. Lucas, “Imperialist Competitive Algorithm: An Algorithm for Optimization inspired by Imperialistic Competition,” IEEE Congress on Evolutionary Computation, vol. 7, pp. 4661-4667, Sept. 2007.
[16]   | L.D.S. Coelho, L. D. Afonso and P. Alotto, “A modified imperialist competitive algorithm for optimization in electromagnetic,” vol. 48, no. 2, pp.579-582, Feb. 2012.
[17]   | G. Zhang, J. Peng, “Improved ICA for Constrained optimization,” International forum on computer science technology and application, pp.204-207, 2009.
[18]   | J.A. Nelder, R. Mead, “A simplex for function minimization,” The computer journal, vol. 7, pp. 308-313, 1965
[19]   | Z.H. Wang, H. Chen, and G.A. Taylor, “Volt/Var Control in Distribution Systems Using a Time-Interval Based Approach,” IEE proceedings Generation Transmission & Distribution, vol. 150, no. 5, pp. 548-554, September2003.
[20]   | F. Viawan and D. Karlsson, “Voltage and Reactive Power Control in Systems with Synchronous Machine-Based Distributed Generation,” IEEE Transactions on Power Delivery, vol. 23, no. 2, pp. 1079-1087, 2008.