DG Installation in Distribution System for Minimum Loss

Similar documents
Academic Editor: João P. S. Catalão Received: 10 February 2016; Accepted: 6 April 2016; Published: 14 April 2016

Computational Analysis of IEEE 57 Bus System Using N-R Method

ECONOMIC OPERATION OF POWER SYSTEMS

Optimization Methods MIT 2.098/6.255/ Final exam

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

Information-based Feature Selection

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled

subject to A 1 x + A 2 y b x j 0, j = 1,,n 1 y j = 0 or 1, j = 1,,n 2

Unbalanced Radial Distribution System Load Flow and Voltage Profile Enhancement in the presence of Distributed Generators

IP Reference guide for integer programming formulations.

A NEW APPROACH TO SOLVE AN UNBALANCED ASSIGNMENT PROBLEM

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

SNAP Centre Workshop. Basic Algebraic Manipulation

THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS

Analytic Theory of Probabilities

Introduction to Optimization Techniques. How to Solve Equations

TR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT

An Introduction to Randomized Algorithms

CS322: Network Analysis. Problem Set 2 - Fall 2009

A numerical Technique Finite Volume Method for Solving Diffusion 2D Problem

Deterministic Model of Multipath Fading for Circular and Parabolic Reflector Patterns

Study on Coal Consumption Curve Fitting of the Thermal Power Based on Genetic Algorithm

Polynomials with Rational Roots that Differ by a Non-zero Constant. Generalities

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Open Circuit Fault Analysis of Electrical Power System Using MATLAB

NUMERICAL METHODS FOR SOLVING EQUATIONS

The improvement of the volume ratio measurement method in static expansion vacuum system

Number of fatalities X Sunday 4 Monday 6 Tuesday 2 Wednesday 0 Thursday 3 Friday 5 Saturday 8 Total 28. Day

c. Explain the basic Newsvendor model. Why is it useful for SC models? e. What additional research do you believe will be helpful in this area?

Kinetics of Complex Reactions

Reliability and Queueing

ELECTRICAL energy is continuously lost due to resistance

SEQUENCES AND SERIES

mx bx kx F t. dt IR I LI V t, Q LQ RQ V t,

Optimally Sparse SVMs

Random Walks on Discrete and Continuous Circles. by Jeffrey S. Rosenthal School of Mathematics, University of Minnesota, Minneapolis, MN, U.S.A.

MULTI-OBJECTIVE GENERATION DISPATCH USING PARTICLE SWARM OPTIMISATION WITH MULTIPLE FUEL OPTION

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Linear Regression Demystified

Session 5. (1) Principal component analysis and Karhunen-Loève transformation

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

THE EFFECT OF DISPERSED GENERATION UNITS WITH AND WITHOUT CAPACITOR BANKS ON THE SYSTEM LOSSES AND THE VOLTAGE REGULATION.

OPTIMAL PIECEWISE UNIFORM VECTOR QUANTIZATION OF THE MEMORYLESS LAPLACIAN SOURCE

Exact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method

Optimum allocation of Renewable DG sources and Synchronous capacitor simultaneously using PSO

Balancing. Rotating Components Examples of rotating components in a mechanism or a machine. (a)

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan

The Simplex algorithm: Introductory example. The Simplex algorithm: Introductory example (2)

l -State Solutions of a New Four-Parameter 1/r^2 Singular Radial Non-Conventional Potential via Asymptotic Iteration Method

Analytical solutions for multi-wave transfer matrices in layered structures

The target reliability and design working life

THE APPEARANCE OF FIBONACCI AND LUCAS NUMBERS IN THE SIMULATION OF ELECTRICAL POWER LINES SUPPLIED BY TWO SIDES

Scientific Review ISSN(e): , ISSN(p): Vol. 4, Issue. 4, pp: 26-33, 2018 URL:

Lecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables

Scheduling under Uncertainty using MILP Sensitivity Analysis

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

Example 3.3: Rainfall reported at a group of five stations (see Fig. 3.7) is as follows. Kundla. Sabli

The Local Harmonious Chromatic Problem

Appendix: The Laplace Transform

Numerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials

ROLL CUTTING PROBLEMS UNDER STOCHASTIC DEMAND

A comparative study of a system of Lotka-Voltera type of PDEs through perturbation methods

11 Correlation and Regression

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Fastest mixing Markov chain on a path

Linear Support Vector Machines

Uniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations

Polynomial Functions and Their Graphs

CHEE 221: Chemical Processes and Systems

1. pn junction under bias 2. I-Vcharacteristics

Some Properties of the Exact and Score Methods for Binomial Proportion and Sample Size Calculation

Time-Domain Representations of LTI Systems

MEI Casio Tasks for Further Pure

Upper bound for ropelength of pretzel knots

Effect of Temperature Variations on Strain Response of Polymer Bragg Grating Optical Fibers

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014

THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.

Ma 530 Introduction to Power Series

Application of Homotopy Perturbation Method for the Large Angle period of Nonlinear Oscillator

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

OPTIMIZED SOLUTION OF PRESSURE VESSEL DESIGN USING GEOMETRIC PROGRAMMING

CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS

Senvion SE Franz-Lenz-Str. 1, Osnabrück, Germany

G r a d e 1 1 P r e - C a l c u l u s M a t h e m a t i c s ( 3 0 S )

MAT 271 Project: Partial Fractions for certain rational functions

Stochastic Simulation

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

U8L1: Sec Equations of Lines in R 2

MAXIMALLY FLAT FIR FILTERS

SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

Streamfunction-Vorticity Formulation

Chapter 7: The z-transform. Chih-Wei Liu

ADVANCED SOFTWARE ENGINEERING

Linear Programming and the Simplex Method

Elementary Statistics

Assignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.

Integer Programming (IP)

Transcription:

DG Istallatio i Distributio System for Miimum Loss Aad K Padey Om Mishra Alat Saurabh Kumar EE, JSSATE EE, JSSATE EE, JSSATE EE, JSSATE oida,up oida,up oida,up oida,up Abstract: This paper proposes optimal placemet ad optimal size of DG system i the distributed system. Earlier proposal gives aalytical expressio for optimal placemet of DG that is capable of delivery real power oly. This research paper has proposed two other type of DG system also (type 3 ad type 4). 14 test bus systems with stadard IEEE data is used. Result shows optimal locatio for placemet of DG ad its optimal size with some costraits. The result required less computatio ad are verified with the help of illustrated graphs. Keywords: Optimal Locatio, Optimal Size, Power Losses, Distributed Geeratio. I. ITRODUCTIO The objective of power system is to supply electrical eergy to each cosumer with rated coditio with miimum losses ad maximum power trasfer alog with maitaiig stability of the system. There are two types of power techologies: Traditioal techologies ad o-traditioal techologies. Traditioal techologies are o-reewable techologies such as thermal power plat. o-traditioal techologies iclude fuel cells ad reewable eergy based techologies such as wid eergy, solar eergy, geothermal eergy, ocea eergy. For icreasig demad of eergy with miimum cost has led to developmet of may ew techologies. These techologies are also very importat because of growig evirometal cocer aroud the whole world. The fear of depletio of covetioal source of eergy has led to the developmet of o covetioal techologies. Amog them owadays DG system has attracted several researchers to fid ew ideas about how it ca be utilised effectively. Amog these optimal size ad optimal placemet of DG system is very importat. DG (distributed geeratio ) is small scale power statio whose rage may vary from few KW to 50 MW.They are geerally placed at the cosumer side or at the ed of distributio etwork.there are four types of DG system : 1. Type 1 DG, capable of ijectig active power. 2. Type 2 DG, capable of ijectig reactive power 3. Type 3 DG capable of ijectig active ad reactive power. 4. Type 4 DG capable of ijectig active power ad cosumig reactive power. DG system improves reliability of the system. But has also led to may problems. It has foud that if DG is placed at wrog positio i the distributio system the effectiveess of DG decreases. Therefore it is very importat for oe to decide where the DG should be used ad it should be what size. This research paper focuses o this problem oly takig ito cosideratio of all the four types of II. LITERATURE REVIEW I [1] for reducig loss approximate power flow techique was used. As this method is approximate method therefore it ca ot be used i the practical large distributio system. I [2] sigle loop method was used for solvig etwork cofiguratio. I [3] 2/3 are golde rule is used. This method is very effective for uiform load. But this is applicable to ot applicable to o-uiformly distributed load.this rule is simple ad easy to apply. I[4] approach has bee made to fid appropriate locatio to place that DG system i radial as well i loop system this method helps i miimizig losses but here optimal sizig is ot cosidered. Because i practical case differet sizes of DG are preset ad due to ecoomic reaso sometime it becomes difficult to use DG of same size so this method does ot goes well i practical case. I [5] GA method was used to fid optimal placemet ad size of DG it ca solve may problem. But it requires lots of computatioal time. I [6] it used exact loss formula for fidig optimal locatio of DG system i distributed etwork. Here the two load flow calculatio is performed. First of all losses are calculated without the placemet of DG.The the secod load flow calculatio are performed to fid the losses.but this techique gives solutio for type 1 DG oly. Other techique were also used such as dyamic programmig [7], fuzzy expert system which requires lots of aalysis ad have more computatio time [8] ad [9] TS techique is used for fidig the best locatio of placemet of capacitor oly. It does ot give solutio for iductive loads which is the most commo load i distributio system. III. METHODOLOGY Before the placemet of DG system losses of the distributed system are calculated.ow DG system are placed at differet buses ad losses are calculated with the help of exact loss formula. Amog the various cases the locatio which gives miimum loss gives the optimal positio i the distributed system of particular size. This process is repeated for may sizes withi the give costrait The aalytical expressio is give below. P L = [α ij (P i P j + Q i Q j ) + β ij (Q i P j P i Q j )] i=1 1008

Where α ij = r ij V i V j cos(δ i δ j ); β ij = r ij V i V j si(δ i δ j ); V i i r ij + x ij = Z ij P i Q i Bus voltage load agle Z bus Active Power ijectio at i bus Reactive Power ijectio at i bus o of buses 3. Type 3 DG: Here Q is positive ad the value of P DG ad Q DG will give the value of total apparet power that determies size of 4. Type 4 DG: The value of Q is egative ad here also we ca calculate the total apparet power to determie the size of DG system. With the help of above expressio we ca ot oly fid the optimal size of DG but also its type which will give miimum loss i the distributed etwork P i = P DGi P Di Q i = Q DGi Q Di = ap DGi - Q Di a=ta -1 (Q/P) Puttig the above two equatio i exact loss formula the ew equatio obtaied is: P L = [α ij ((P DGi P Di )P j + (ap DGi Q Di )Q j ) i=1 + β ij ((ap DGi Q Di )P j (P DGi P Di )Q j )] We will differetiate the followig equatio with respect to P DGi P L = 2 [α P ij (P j + aq j ) + β ij (ap j Q j )] = 0 DGi Fially o substitutio ad simplificatio we get P DGi = α ij(p Di + aq Di ) + β ij (ap Di Q Di ) X i ay i a 2 α ii + α ii Where X i = (α ij P j β ij Q j ) Y i = (α ij Q j + β ij P j ) With the help of above equatio we ca estimate the size of DG as follows- 1. Type 1 DG: For type 1 DG value of a=0 because power factor is uity ad size is: P DGi = P Di 1 [β α ii Q Di + (α ij P j β ij Q j )] i 2. Type 2 DG: I this value of a= because power factor=0: Q DGi = Q Di + 1 α i [β ii P Di (α ij Q j + β ij P j ) ] IV. ALGORITHM This research paper uses the stadard IEEE 14 test bus system ad voltage ad load agle of each buses are also preset accordig to IEEE stadards. The followig steps are followed- 1. Z bus formulatio codig has bee doe to fid the Z bus of IEEE stadards. 2. 100 MVA base is take. 3. First of all the total losses of the system is calculated. 4. There is o placemet of 5. May costraits has bee take before the aalysis of system such as size of 6. The loop is performed to fid positio of bus system where miimum losses are preset. 7. We first take miimum size of DG ad the we place at all the 13 bus except the slack bus. 8. Out of all the 13 positio, oly that positio is located which give miimum loss after placemet of 9. This positio is saved i the form of variable. 10. The process is agai repeated util the maximum limit of the size. 11. Several loopig is doe to calculate the optimum positio of differet size of 12. Out of the several size oly that size is chose which gives miimum loss i the distributio system. 13. So we get optimal size ad locatio for type 1 14. Above process is repeated for two other type of DG also (Type 3 ad Type 4 DG). 15. We get miimum losses for all three types of 16. Out of the three cases that DG type is selected which give the miimum loss. V. THE TEST SYSTEM IEEE 14 bus system is used to check the validity of proposed work ad to fid the optimal size ad locatio of The detail about the test bus system is give i the appedix sectio from where aalysis of DG system is doe. 1009

VI. RESULT The below tables shows the values of losses i per uit ad takig 100MVA base ad its variatio o varyig locatio, size ad type of DG ad their correspodig graph are give i appedix sectio. 1. TYPE 1 DG Bus o. Size (p.u.).01.898.896.894.896.895.896.8954.896.896.898.895.896.897.02.891.892.886.892.875.894.877.888.881.884.887.888.887.03.887.884.878.882.879.883.880.881.881.881.881.880.879.04.883.878.870.875.872.875.872.874.873.874.872.871.873.05.878.872.862.868.865.869.865.868.866.865.866.866.863.06.873.866.854.862.857.863.858.861.86.857.859.856.858.07.868.860.846.855.850.856.851.855.852.853.850.852.849.08.863.854.838.848.842.850.844.848.845.847.844.845.841.09.858.849.831.842.836.844.838.842.839.840.837.838.834.1.854.843.823.836.829.838.832.836.832.834.829.831.827 0.855 0.85 5 0.835 0.83 5 0.815 0.81 0.805 10, 0.854 (size MW,losses p.u.) 10, 3 10, 0.836 10, 0.838 10, 0.836 10, 0.832 10, 0.832 10, 0.834 10, 9 10, 9 10, 0.831 10, 7 10, 3 Series1 2. TYPE 3 DG Size Bus o. (p.u.).01.898.897.8958.896.896.897.898.896.896.8958.896.896.895.02.893.891.888.890.891.891.892.890.889.890.888.889.887.03.88.885.881.883.885.884.887.885.882.883.881.882.880.04.884.880.874.877.880.879.883.879.875.877.874.875.873.05.879.874.867.871.874.873.878.873.869.871.867.869.866.06.874.869.861.865.869.867.874.868.874.868.863.864.859.07.870.864.854.860.865.862.870.863.857.859.855.856.853.08.866.859.848.854.86.857.866.858.851.854.849.85.847.09.862.854.842.849.856.852.863.853.846.848.846.844.841.1.858.849.836.843.851.848.860.849.841.843.837.838.835 1010

0.92 0.9 LOSSES 10, 0.894 (size MW,losses p.u 0.88 10, 0.886 10, 0.886 10, 0.878 10, 0.87 10, 0.854 10, 2 10, 1 10, 6 10, 0.839 10, 0.834 10, 0.831 10, 3 Series1 0.8 0.78 BUS Size(p.u.) 3. TYPE 4 DG Bus o..01.898.896.894.896.894.896.893.895.895.895.895.895.894.02.893.890.885.888.884.888.883.887.887.887.887.887.886.03.887.883.881.887.881.873.874.879.879.879.879.879.878.04.882.877.867.874.865.874.864.872.872.872.872.872.871.05.878.871.858.867.856.867.855.865.865.865.864.865.863.06.873.865.850.860.848.860.847.857.858.858.857.858.855.07.868.859.842.854.839.854.838.851.851.851.85.851.848.08.864.853.834.847.831.847.830.844.844.845.843.844.841.09.860.848.826.840.823.841.822.838.838.838.837.838.834.1.855.842.818.834.815.835.814.832.832.832.831.831.828 LOSSES 0.85 0.83 10, 0.855 10, 2 10, 0.834 10, 0.835 10, 0.818 10, 0.815 (size MW,loss p.u.) 10, 0.832 10, 0.831 10, 0.831 10, 0.832 LOSSES VS POSITIO 10, 0.831 10,8 0.81 10,0.814 0.8 0.79 BUS 1011

VII.COCLUSIO The loss i the 14 bus system without the use of DG is calculated to be 90.03 MW. Usig type 1 DG the loss is calculated to be 82.36 MW ad the optimal locatio is bus o.4 ad the optimal size is of 10 MW Usig type 3 DG the loss is calculated to be 83.53 MW ad the optimal locatio is bus o.14 ad the optimal size is of 10 MW Usig type 4 DG the loss is calculated to be 81.46 MW ad the optimal locatio is bus o.8 ad the optimal size is of 10 MW From above cases type 4 DG is suggested at locatio bus o.8 of size 10 MW which gives miimum loss of 81.46 MW. VIII. REFERECES [1] M. E.Bara ad F. F.Wu, etwork recofiguratio i distributio systems for loss reductio ad load balacig, IEEE Tras. Power Del., vol. 4, o. 2, pp. 1401 1407, Apr. 1989. IX. [2] J. Y. Fa, L. Zhag, ad J. D. McDoald, Distributio etwork recofiguratio: [3] J. V. Schmill, Optimum size ad locatio of shut capacitors o distributiofeeders, IEEE Tras. Power App. Syst., vol. PAS-84, o. 9,pp. 825 832, Sep. 1965. [4] C.Wag adm. H. ehrir, Aalytical approaches for optimal placemetof distributed geeratio sources i power systems, IEEE Tras PowerSyst., vol. 19, o. 4, pp. 2068 2076, ov. 2004. [5] A. Silvestri, A. Berizzi, ad S. Buoao, Distributed geeratio plaigusig geetic algorithms, i Proc. It. Cof. Electr. Power Eg.(PowerTech Budapest), Sep., p. 257. [6]. Acharya, P. Mahat, ad. Mithulaatha, A aalytical approach for DG allocatio i primary distributio etwork, It. J. Electr. Power Eergy Syst., vol. 28, o. 10, pp. 669 678, 2006. [7] T. H. Fawzi, S. M. El-Sobki, ad M. A. Abdel-Halim, A ew approach for the applicatio of shut capacitors to the primary distributio feeders, IEEE Tras. Power App. Syst., vol. PAS-102, o. 1, pp. 10 13, Ja. 1983. [8] R. A. Gallego, A. J. Moticelli, ad R. Romero, Optimal capacitor placemet i radial distributio etworks, IEEE Tras. Power Del., vol. 16, o. 4, pp. 630 637, ov. 2001 APPEDIX Bus o. Bus code Voltage magitude Agle degrees 1. STADARD IEEE DATA Load Geerator Ijected MVAR MW MVAR Qmi Qmax 1 1 1.06 0 30.38 17.78 40-40 0 0 0 2 2 1.045 0 0 0 232 0-40 50 0 3 2 1.01 0 131.88 26.6 0 0 0 40 0 4 0 1 0 66.92 10 0 0 0 0 0 5 0 1 0 10.64 2.24 0 0 0 0 0 6 2 1.07 0 15.68 10.5 0 0-6 24 0 7 0 1 0 0 0 0 0 0 0 0 8 2 1.09 0 0 0 0 0-6 24 0 9 0 1 0 41.8 23.24 0 0 0 0 0 10 0 1 0 12.6 8.12 0 0 0 0 0 11 0 1 0 4.9 2.52 0 0 0 0 0 12 0 1 0 8.54 2.24 0 0 0 0 0 13 0 1 0 18.9 8.12 0 0 0 0 0 14 0 1 0 2 7 0 0 0 0 0 2. IEEE-14 BUS SYSTEM 1012