Abstract

This paper intends to enumerate the impact of distributed generation (DG) on distribution system in terms of active as well as reactive power loss reduction and improved voltage stability. The novelty of the method proposed in this paper is the simple and effective way of sizing and siting of DG in a distribution system by using two-port Z-bus parameters. The validity of the method is verified by comparing the results with already published methods. Comparative study presented has shown that the proposed method leads existing methods in terms of its simplicity, undemanding calculation procedures, and less computational efforts and so does the time. The method is implemented on IEEE 69-bus test radial distribution system and results show significant reduction in distribution power losses with improved voltage profile of the system. Simulation is carried out in MATLAB environment for execution of the proposed algorithm.

1. Introduction

With the development of economy, load demands in distribution networks are sharply increasing. Hence, the distribution networks are operating more close to the voltage instability boundaries. The decline of voltage stability margin is one of the important factors which restrict the increase in load served by distribution companies [1]. Therefore, it is necessary to consider voltage stability with the integration of DG units in distribution systems. The insertion of DG presents opportunities as it leads many technical as well as economical benefits along with voltage stability attained due to reduction of line currents. Reduction of line currents also results in the line loss reduction. Planning issues, regulatory framework, and the availability of resources limit distribution network operators (DNOs) and developers in their ability to accommodate distributed generation but governments are incentivizing low carbon technologies, as a means of meeting environmental targets and increasing energy security. This momentum can be harnessed by DNOs to bring network operational benefits through improved voltage profile and lower line losses delivered by investment in DGs [2]. The main hurdles for the DNOs are the implementation and reliability of the loss minimization technique; however researchers are trying hard to locate better techniques to exploit all possible benefits of DG [2].

This paper mainly focuses on the investigation of simple and efficient analytical approach for siting and sizing of DG for insertion in distribution system and evaluation of its performance in terms of system loss reduction and voltage profile correction. Two-port Z-bus parameters are utilized to find the expressions for designing DGs to be installed in distribution network. Optimal location and capacity of DG are calculated by the proposed method and results are compared with well-known exhaustive load flow (ELF) and results generated by software voltage stability and optimization (VS&OP) package [3]. Results show that proposed method follows similar prototype of those methods but with much less computational efforts and time.

Many methods of DG allocation, available in the literature, require complicated equations to solve which may further require calculations of many subcoefficients and rigorous iterative steps. Thus all those methods [412] gratuitously make the process time consuming and tedious especially for large systems. The proposed method directly utilizes the two-port Z-bus equations and converts them into power form to get both types of losses (active and reactive power losses) in the distribution system.

2. Proposed Methodology

This paper proposes an analytical method based on two-port Z-bus equations. It is observed that Z-bus equations in the modified form can give the optimal size of DG to be installed at each location and the proportional losses.

2.1. Two-Port Z-Bus Equations

Consider a distribution system as shown in Figure 1.

Two-port Z-bus equations for above system are The negative sign indicates the opposite direction of as shown in Figure 1. Hence for a radial distribution system the two-port feeder equations can be written as where , , , and are column matrices of load and generator voltages and currents, respectively. and are generator and load buses, respectively.

The ratio of source to load voltage at light or no load condition can be written as [] & [] are corresponding partitioned portions of network matrix. The above relation gives a transmission voltage drop (TVD) which is similar to the ILF matrix investigated in [13]. The TVD matrix will be of the dimension .

And

The TVD matrix thus formed gives the factor by which generator voltage is reduced to load voltage due to the drop across the impedance present in the power flow path from generator to respective load. Thus it can give the power contribution of each generator to each load as per the impedance encountered in the path so as to attain the total demand with maximum efficiency. The implication of TVD matrix is explained in the next section.

2.2. Significance of TVD Matrix

It is assumed that the 7-bus sample radial distribution system as shown in Figure 2 has load A to E = 10 MW each which is fed by substation (S.S.) and DG. Each section of the line is of the same length, that is, 100 Km with line parameters in per unit per 100 km being = 10 and = 25. To demonstrate the utility of TVD matrix three cases with different line conditions are considered.

Case 1. Each section of the line has the same line parameters ( and per unit).

Case 2. Values of line parameters between sections 2 and 3 are doubled.

Case 3. Values of line parameters between sections 3 and 6 are doubled.

For the cases mentioned above three TVD matrices are evaluated and are given as

Case 1. Since load A is at 100 Km from S.S. and 300 Km from DG, that is, three times of 100 Km, the corresponding factors of (TVD) matrix are 0.75 and 0.25, respectively. Similarly load E is at 300 Km from DG and 100 Km from S.S.; the corresponding TVD factors are 0.25 and 0.75, respectively. Whereas the loads B, C, and D are at the same distance from both S.S. and DG, the corresponding TVD factors are 0.50 for each.

Case 2. Load A is at the same distance from S.S. but the DG is now 400 Km away from it; the corresponding TVD factors for this condition are 0.80 and 0.20, respectively. Similarly load E is at the same distance from DG whereas S.S. is now 400 Km away from it, and the corresponding TVD factors are 0.20 and 0.80, respectively. The loads B, C, and D are at 300, 400, and 500 Km, respectively, from S.S. and 200, 300, and 400 Km, respectively, from DG. In other terms it can be said that these loads are at 60% of distance from S.S. and at 40% of distance from DG and the corresponding TVD matrices for each of them are 0.40 and 0.60, respectively.

Case 3. Load E is at the same distance from DG but the S.S. is now 400 Km away from it; the corresponding TVD factors for this condition are 0.80 and 0.20, respectively. Similarly load A is at the same distance from S.S whereas DG is now 400 Km away from it, and the corresponding TVD factors are 0.20 and 0.80, respectively. The loads B, C, and D are at 200, 300, and 400 Km, respectively, from S.S. but at 300, 400, and 500 Km, respectively, from DG. In other terms it can be said that these loads are at 40% of distance from S.S. and at 60% of distance from DG and the corresponding TVD matrices for each of them are 0.60 and 0.40, respectively.

From the analysis of three cases and the corresponding TVD factors it can be observed that the TVD matrix gives the accurate idea of the proportion of impedance encountered in the path between respective load and the generator feeding that load. Thus TVD factor can play a decisive role in recommending the proportion of power contribution from each generator present in the system to the loads which will give the minimum power losses. The detailed analysis is furnished in Tables 1 and 2.

Thus from Table 2 it is clear that if DG and S.S. are scheduled according to TVD factor, the voltage profile improvement as well as power loss reduction is significant whereas the load sharing with different values will result in reduced power loss minimization percentage as well as declined voltage profile. Hence TVD factor can be considered as the decisive factor for DG capacity planning.

2.3. Optimal Loss

Optimal loss calculation assortment of equations can be developed as given below.

Multiplying both sides by , (2) will be rewritten as Thus, Pre-multiplying equation (8) by , and rewritten as; Thus, (9) can be written in power form as where where Corresponding partitioned portion of matrix. Corresponding partitioned portion of matrix. Column matrix of generator voltages. Column matrix of load currents. unity matrix of dimension , with all diagonal elements replaced by respective generator voltages.

Thus, where is total number of generators and is the total number of load buses.

Equation (13) will give the power contribution of each generator to the system load. So it may also be used for calculation of DG capacity. Equation (14) will give total active power losses supplied by available sources in the system. This equation is used for loss calculation of the system with DG installed at any particular load bus and, thus, it can be used for finding the optimal location for DG insertion.

2.4. Optimal DG Capacity

As explained in Section 2.2. transmission voltage drop obtained by (4) plays an important role in calculating the approximate value of DG capacity. Equation (13) can be used for developing the equation for DG capacity evaluation which is given as

From the matrix given by (15) the power contribution of the generators to each load present in the system can be obtained which is actually similar to the T-index value explained in [14]. Here each row of the matrix gives equation for evaluating the power supplied by respective generator to loads present in the system and each individual term of the equation indicates the fraction of load contributed by it. Thus the equation explored at each row will give the maximum capacity of the respective generator.

Now if DG is installed at th load bus then that load bus is considered as the generator bus so the total number of generator buses is increased by one and the number of load buses is decreased by one. So the modified matrix with new dimension as will be Thus the desired capacity of the DG is the summation of terms of power contribution of DG to each load where is the total number of generator buses and is the total number of load buses.

2.5. Optimal Location

Steps to find the optimum location are given below.

Step 1. Calculate the system losses without DG by using (14).

Step 2. Put DG at every load bus and calculate the losses by updating (14) and corresponding DG capacity by using (17).

Step 3. Compare the losses with the base case (without DG).

Step 4. Prepare the priority list in ascending order of losses calculated for each location of DG.

Step 5. Top ranking bus is considered as the optimum location for DG installation.

Step 6. The capacity of DG calculated for that location is considered as the optimum capacity of DG at that location. From (17) corresponding optimal DG capacity at optimal location can be evaluated.

3. Flowchart for Proposed Method

See Figure 3.

4. Results and Analysis

4.1. Test System

The proposed method is tested on 69-bus radial distribution system with a total real and reactive load of 3.8 MW and 2.69 MVAr, respectively [15]. An analytical software tool has been developed in MATLAB environment to execute the proposed algorithm.

4.2. Assumptions

(1) Designing of DG is performed at peak load only. (2) Maximum active and reactive power limit of DG for different test systems is assumed to be equal to the total active load of the system. (3) The lower and upper voltage thresholds for DG are set at 0.95–1.05 pu.

4.3. Simulation Results

Table 3 and Figure 4 show the comparative analysis of three methods of loss minimization for 69-bus system. Results show that the proposed method follows the similar trend line as that of the results generated by VS&OP package as well as ELF. The priority list shown in Table 3 indicates that bus 61 is the best location for installing the DG. The next suitable locations are buses number 60, 62, 63, and so forth. Other two methods can calculate only active power losses whereas proposed method gives both active and reactive loss reductions. Figure 5 indicates the optimal capacity of DG at corresponding location which will give the minimum possible power losses. Hence it is observed that if DG is installed at bus 61 with capacity of 2.03 MWAtts and 1.22 MVAR recommended by proposed method then maximum loss reduction is possible. Since DG installation at suitable location reduces the power flows of many lines in the system which in turn reduces the overall loss and improves the corresponding node voltages, it is seen that DG insertion at the appropriate location can improve the voltage profile of the system which can be observed from Figure 6. Table 4 gives the summarized results of all the three methods which indicate that the proposed method shows enhanced optimality of solution. Table 5 indicates the voltage profile of the system for case (i) without DG and case (ii) with DG installed at optimum location, that is, bus 61. It indicates that the proposed method suggests the location as well as capacity of DG which will give improved voltage profile as compared to other methods.

5. Conclusion

This paper presents a novel method which uses two-port Z-bus equations for finding the optimal location and optimal size at that location. It avoids the time consuming and cumbersome iterative approach for handling the undemanding problem of designing the new DG to be installed in radial distribution system. The proposed method for designing DG to install in distribution system is tested on 69-bus test distribution system and results are verified by comparing with well-known exhaustive load flow solutions and results generated by already published VS&PO package. The proposed method provides an efficient tool for distribution system loss calculation, for loss minimization, and for designing the DG to install in system to attain enhanced voltage stability. The method is applied for single DG allocation but it can also be implemented for multiple DG allocation.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.