Comparison of Wheeling Cost using Power Flow Tracing Methods in Deregulated Electric Power Industry. K.Hema Lalitha Student M.Tech MVGR College of Engineering Vizianagaram I.Kranthi Kiran Associate Professor MVGR College of Engineering Vizianagaram ABSTRACT: The power industry across the World has experienced a lot of restructuring from regulated to deregulated environment to improve the technical and commercial efficiency. The purpose of deregulation is to restructure the power industry so that power production and retail sales will be competitive, leading to Independent power producers (IPPs).In deregulated power sector, the overall power system is restructured in three main parts - generation, transmission and distribution of electricity. Electricity transmission and wheeling service pricing has become a more complex and more important task. Pricing of transmission services plays a vital role in determining whether the providing transmission services are economically beneficial to both wheelers and customers. The wheeling cost depends on how much power is being moved through a particular line. To know the power flow in this case, Bialek s method and Kirschen s method has been considered. There are several methods to determine the wheeling cost. In this paper, the power flow based MW-KM method is used to calculate the wheeling cost of IEEE 14 bus system and the results are compared in both tracing methods. Keywords:Deregulation, Power tracing, Transmission pricing, Wheeling cost. 1. INTRODUCTION: The electric power industries over the years are vertically integrated utilities which are the only electricity providers in the region. It was difficult to segregate the costs incurred in generation, transmission or distribution. So, the utilities often charged their customers on an average tariff rate depending on their aggregated cost during a period. Therefore, deregulation occurred. Electric deregulation is the process of changing rules and regulations that control the electric industry for providing customers the choice of electricity suppliers by allowing competition leading to Independent Power Producers (IPPs). In this process, some new entities are expected to appear and there are: Generation Companies (GenCos.), Transmission Companies (TransCos.), Distribution Companies (DisCos.), Independent Power producer (IPP), Independent System Operator (ISO), Power Exchange (PX), Retail Energy Service Companies (RESCos.). Pricing of transmission services plays a crucial role in determining whether the providing transmission services are economically beneficial to both the wheelers and the customers. Wheeling can be defined as the transfer of electrical power from one network to another to deliver power of and for another entity or entities. To evaluate the costs of wheeling transaction, various salient methods are there like Postage Stamp method, Contract Path method, Boundary flow method and Flow based methods. In those methods, Power Flow based MW-km method is used here which is accurate and reliable.to know the power transfer between individual generators to lines and generators to loads, there are different power flow tracing methods, but in this paper we use Bialek s Power Flow Tracing method (Node method) and Kirchen s power flow tracing method(common method) to determine the power flows in an IEEE 14 bus system. The wheeling cost is then calculated by using Power flow based MW-km method for this case. 861
2. BIALEK S POWER FLOW TRACING METHOD: Bialek s power tracing technique describes that which generators are supplying a particular load, how much use each generator is making use of a transmission line and what is each generator s contribution to the system losses. It is possible to calculate how much of a particular generators output supplies a particular load and how much of a particular load is supplied by particular generator. The main principle used in this power flow tracing method is Proportionate Sharing principle. In Bialek s method, there are two algorithms for power tracing: Upstream tracing algorithm and Downstream tracing algorithm. In this paper, only Upstream tracing is used so downstream tracing is not discussed. 2.1 Proportionate Sharing Principle: The proportional sharing principle is based on Kirchhoff s current law and is topological in nature. It assumes that the network node is a perfect mixer of incoming flows. Practically the only requirement for the input data is that Kirchhoff s current law must be satisfied for all the nodes. Figure 1: Example of Proportionate Sharing Principle Here, four lines are connected to node i, two inflows and two outflows. The total power flow through the node is P i=40+60= 100 MW of which 40% is supplied by line j-i and 60% by line k-i. Hence, the 70 MW out flowing in line i-m consist of 70 x 40/100=28 MW supplied by line j-i and 70 x 60/100=42 MW supplied by line k-i. In the same way, 30MW out flowing in line i-l 30 x 40/100=12 MW consists of supplied by line j-i and 30 x 60/100 = 18 MW supplied by line k-i. 2.2 Upstream tracing algorithm: The total flow and the inflow to the i th bus, is the sum of all the power inflows through the lines connected to that bus and the bus power injection at that bus by generators. It is given by equation 1. P = P + P for i=1,2,.n; j=1,2, n (1) where i= Sending end bus number, j= Receiving end bus number, n=number of nodes/buses, α =Set of nodes supplying directly node i (i.e. power must flow towards node i in the relevant lines), P =Line flow into node i in line j-i, P = Generation at node i If the line losses are neglected, P = P Equation 1 can further be expanded to the following given by equation 2 P = P + P (2) 862
Let C = to express the relationship between line flow and the nodal flow at the j th node, using proportional sharing principle. P = C P. Substituting this in 2 which gives equation 3 P = C P + P (3) Which on arrangement becomes P C P = P or A u P=P G P is the vector of gross nodal flows; P G is the vector of nodal generations, A u is the Upstream matrix of order (n n). The (i, j) element of A u is given by the following equation 4 1 for i = j [A ] = C = for j α 0 otherwise -1 If A exists, then P=A u P G and its i th element is given by equation 5 P = [A ] P for i=1,2,.. n (5) for i=1,2,.n; j=1,2, n (4) This equation shows the contribution of k th system generator to i th nodal power flow. A line out flow in line j-i from node i can be therefore calculated using proportional sharing principle, as given by equation 6 P = [A ] P for i=1,2,.n (6) This equation gives the contribution of k th generator to power flow in each line j-i. Finally, load demand at the i th bus, applying the proportional method, is given by equation 7 as P = [A ] P (7) P = P P P This equation shows that the contribution of the k th generator to the i th load demand and can be used to trace where the power of a particular load comes from. 3. KIRSCHEN POWER TRACING METHOD: This method is proposed by Daniel Kirschen, Ron Allan and Goran Strbac. Firstly, the active power flows are calculated using a solved power flow computation. This method arranges the buses and branches of the network into homogeneous groups. Once this organization is over, it is then possible to represent the state of the system by a directed, acyclic graph. Based on the graph and the assumption, the contribution of each generator to the branches and the loads can then be computed. From the starting of a power flow solution, first identify the busses and then find the set of busses supplied to the same generators. The method involves some new concepts as: Domains, Commons, Links. Domain of a generator is defined as set of buses that are reached by the power produced by the generator. A common is defined as a group of neighboring buses supplied by the same generators. One bus belongs to one common only. Unconnected sets of buses supplied by same generators are treated as separate commons. The number of generators that supply a common is defined as the rank of the common, and its 863
range will be between one and the total number of generators in the system. The lines that connect two different commons are defined as links. The power flows in the lines of a particular link are in the same direction, always from a common of lesser rank to common of higher rank. In this method, the assumption is that the proportion of inflow from generators to common is equal to the proportion of outflows to that common. The contribution of each generator to each load in the common can be calculated recursively by using the following relations below: F = C. F (8) I = F (9) C = (10) Where, C ij: contribution by generator i to the load and external flow of common j C ik: contribution by generator i to the load and external flow of common k F jk: flow from common j to common k through the link F ijk: flow from common j to common k through the link, coming from common i I k : internal flow of common k. 4. WHEELING COST CALCULATION USING POWER FLOW BASED MW-KM METHOD: Wheeling can be defined as the transfer of electrical power through transmission and distribution lines from one utility service area to other. Wheeling charge is the price per unit (MW-hour) that a transmission owner receives for using its transmission lines to transfer power from seller to buyer i.e., from generating stations or IPPs to consumers. Wheeling is more important for Independent Power Producers as they do not own transmission lines. They need some path to transfer the generated power. So, wheeling charge computation become more prominent in deregulated power sector. Among all wheeling cost computation methods, Power flow based MW-KM method is selected, because the actual power flows in the lines are taken into consideration which gives the actual transmission charges for every network. The step by step procedure to calculate the transmission charges is described below. The transmission charges for power flow in each line is calculated by multiplying the magnitude of power flow in each line due to each transaction (MW t,l in MW) with the length of each line (L l in km) and a predetermined weighing factor which reflects the unit rate (W l in Rs./MW/km) and summed over for all lines in the network. It is given as: MWkm = W MW, L (11) Then determine the total cost for each generator in the network. The cost allocated for each generator for power transaction i.e., the share of the total transmission network capacity cost, TC, allocated to transaction t is calculated using the formula below: Total transmission of each generator Transaction cost = total line cost Total transmission cost of all generators TC = TC,, (12) Where, TC t =cost allocated to transaction t TC = total cost of all lines in Rs. L k = length of line k in km W k =cost per MW per unit length of line k 864
MW t,k = flow in line k, due to transaction t T - set of transactions K - set of lines Total transmission network capacity cost is calculated as follows: TC = W P, L (13) Finally, calculate the per unit cost for each generation using equation below: per unit cost Rs. MW = (14) 4. CASE STUDY: IEEE 14 bus system consists of 14 buses,20 lines,2 generator buses and 11 load buses. Bus 1 is considered as reference bus. Figure-2 shows the line diagram of IEEE 14 bus system. Firstly, the load flow is done for the base case IEEE 14 bus system to know the power flows in each line and load flows through Matpower 4.1 software in MATLAB. Then, by using the Bialek s Upstream looking algorithm, the contributions of each generator in the system to power flow through each line in the network and loss in each line, contribution of each generator to each load in the system are determined using equations (1) to (7). By using equations (8) to (10 ), the contributions of each generator to each line and each load in the system are calculated in Kirschen s method. The wheeling cost is computed using equations (11) to (14). Figure 2: Line Diagram of IEEE 14 bus system 865
5. RESULTS: Table 1: Generator contribution for power flow in each line S.no. Line (MW) (MW) 1 1-2 156.8829 0 2 1-5 75.51038 0 3 2-3 58.35816 14.87942 4 2-4 44.72746 11.40404 5 2-5 33.08151 8.434702 6 3-4 20.47688 3.182258 7 4-5 57.16409 4.508562 8 4-7 24.29808 3.776101 9 4-9 13.91696 2.162798 10 5-6 40.86433 3.222992 11 6-11 6.815718 0.537559 12 6-12 7.216869 0.569198 13 6-13 16.45052 1.297461 14 7-8 0 0 15 7-9 24.29808 3.776101 16 9-10 4.524424 0.703129 17 9-14 8.158491 1.26789 18 10-11 3.520259 0.277645 19 12-13 1.496248 0.11801 20 13-14 5.231259 0.412592 Table 2: Total transmission network capacity cost Maximum power flow in the line (MW) TC = W P, L ( Rupees ) S.no Line P max W l =0.5 Rs./MW/km W l =1 Rs./MW/km 1 1-2 156.8829 1764.9325 3529.865 2 1-5 75.51038 3205.4157 6410.8314 3 2-3 73.23758 2761.0567 5522.1135 4 2-4 56.1315 1880.4051 3760.8102 5 2-5 41.51622 1376.2625 2752.5251 6 3-4 23.65914 770.10485 1540.2097 7 4-5 61.67265 493.3812 986.7624 8 4-7 28.07418 557.6935 1115.387 9 4-9 16.07976 849.57399 1699.148 10 5-6 44.08732 1055.4505 2110.9009 866
11 6-11 7.353277 278.6892 557.3784 12 6-12 7.786067 379.57077 759.14153 13 6-13 17.74798 439.26243 878.52485 14 7-8 5.62E-14 0 0 15 7-9 28.07418 588.15399 1176.308 16 9-10 5.227552 84.163595 168.32719 17 9-14 9.426381 485.45862 970.91725 18 10-11 3.797904 138.81338 277.62676 19 12-13 1.614258 61.503221 123.00644 20 13-14 5.643851 375.31609 750.63218 Total 17545 35090 Table 3: Transmission cost of each generator for power flow in each linefor W k=1rs./mw/km Kirschen method Bialek method S.no Line 1 1-2 3529.865 0 3529.865 0 2 1-5 6410.831 0 6410.831 0 3 2-3 4400.205 1121.908 4400.205 1121.908 4 2-4 2996.74 764.0705 2996.74 764.0705 5 2-5 2193.304 559.2208 2193.304 559.2208 6 3-4 1332.993 207.2163 1333.045 207.165 7 4-5 915.6414 71.12099 914.6254 72.137 8 4-7 965.3254 150.0616 965.3625 150.0245 9 4-9 1470.549 228.5995 1470.605 228.5429 10 5-6 1958.758 152.1434 1956.584 154.3168 11 6-11 517.2053 40.1731 516.6314 40.747 12 6-12 704.4263 54.71519 703.6447 55.49684 13 6-13 815.2051 63.31976 814.3005 64.22432 14 7-8 0 0 0 0 15 7-9 1018.05 158.2578 1018.089 158.2186 16 9-10 145.6808 22.64635 145.6864 22.64075 17 9-14 840.2923 130.625 840.3246 130.5927 18 10-11 257.6168 20.00997 257.3309 20.29583 19 12-13 114.1407 8.8657 114.0141 8.992353 20 13-14 696.5303 54.10188 695.7574 54.87476 Total 31283.36 3807.056 31276.95 3813.469 867
Table 4: Transmission cost of each generator for power flow in each linefor W k=0.5 Rs./MW/km Kirschen method Bialek method S.no Line 1 1-2 1764.933 0 1764.933 0 2 1-5 3205.416 0 3205.416 0 3 2-3 2200.103 560.9541 2200.103 560.9541 4 2-4 1498.37 382.0353 1498.37 382.0353 5 2-5 1096.652 279.6104 1096.652 279.6104 6 3-4 666.4967 103.6081 666.5223 103.5825 7 4-5 457.8207 35.56049 457.3127 36.0685 8 4-7 482.6627 75.0308 482.6813 75.01224 9 4-9 735.2743 114.2997 735.3025 114.2715 10 5-6 979.3788 76.07169 978.292 77.15842 11 6-11 258.6026 20.08655 258.3157 20.3735 12 6-12 352.2132 27.3576 351.8223 27.74842 13 6-13 407.6025 31.65988 407.1503 32.11216 14 7-8 0 0 0 0 15 7-9 509.0251 79.12889 509.0447 79.10931 16 9-10 72.84042 11.32318 72.84322 11.32038 17 9-14 420.1461 65.31249 420.1623 65.29633 18 10-11 128.8084 10.00499 128.6655 10.14791 19 12-13 57.07037 4.43285 57.00704 4.496176 20 13-14 348.2651 27.05094 347.8787 27.43738 Total 15641.68 1903.528 15638.47 1906.734 Table 5: Transaction cost of each generator for different values of W k Kirschen method Bialek method W k=0.5 Rs./MW/km 15641.68 1903.528 15638.47 1906.734 W k=1 Rs./MW/km 31283.36 3807.056 31276.95 3813.469 868
(Rupees/MW) Table 6: Per unit cost of each generator Kirschen method (Rupees/MW) (Rupees/MW) Bialek method (Rupees/MW) W k=0.5 Rs./MW/km 67.30694 47.5882 67.29314 47.66836 W k=1 Rs./MW/km 134.6139 95.1764 134.5863 95.33672 6. DISCUSSION ON RESULTS: Table-1 shows the contribution of each generator for power flow through each line in IEEE 14 bus system. Table-2 shows the total transmission network capacity cost for two different values of W l. The value of weighing factor depends on the type of the line (short, medium or long), rating of the line. So, it is assumed as 0.5 and 1 and the results in both cases are observed. Table-3 and table-4 gives the line-wise transmission cost of each generator in both Kirschen and Bialek method. The values clearly says that the contribution of particular generator for power flow in a particular line may vary in both methods but the total power flow in the line is equal to the sum of the contributions from and. Table-5 and Table-6 shows the Transaction cost and Per unit cost of each generator in both power flow tracing methods for different values of W l..the transaction cost of each generator in both the methods is compared. 7. CONCLUSIONS: As the electric power industry is deregulated, generation, transmission and distribution are separated. The main objective is to bring fair and open access in the system. Power tracing became important in the process of unbundling. For power tracing, two efficient methods Bialek s method and Kirschen s method are considered in this paper. These two methods are based on proportionate sharing principle. The contributions of each generator to each line and each load are calculated by Bialek s Upstream looking algorithm and Kirschen s method and the results are compared. For a system with less number of buses, both methods give similar results but as the number of buses increases the difference will be more. Power flow based MW-km method is used in this paper which gives accurate pricing. 8. REFERENCES: [1] POWER SECTOR REFORMS AND RESTRUCTURING IN INDIA -S.A.Khaparde Department of Electrical Engineering Indian Institute of Technology - Bombay, Mumbai, India. [2] M. A. Pai. Computer Techniques in Power System Analysis, Tata McGraw-Hill, New Delhi, 1979. [3] Electric Power Industry Restructuring in India: Present Scenario and Future Prospect,S.N. Singh, Senior Member, IEEE and S.C. Srivastava, Senior Member, IEEE,2004 IEEE International Conference on Electric Utility Deregulation, Restructuring and Power Technologies (DRPT2004) April 2004 Hong Kong. [4] H.H. Happ, Cost of Wheeling Methodologies, IEEE Trans. On Power System, Vol.9, No.1, Feb. 1994, pp147-156 [5] Lee, W.J., Lin, C.H., and Swift, L.D. (2001) Wheeling charge under a deregulated environment, IEEE Trans. Ind. Appl., vol. 37, 1, pp.178 183. [6] J. Park, J. Lim and J. Won, An analytical approach for transmission costs allocation in transmission system, IEEE Trans. on Power Systems, vol. 13, no. 4, pp. 1407-1412, November, 1998. 869
[7] Janusz Bialek, Topological generation and load distribution factors for supplement charge allocation in transmission open access, IEEE Trans. On power systems.vol.12, no.3,august 1997. [8] J. Bialek, D. B. Tam, Tracing The Generators Output, Opportunities and Advances in International Power Generation, 18-20th March 1996. [9] Daniel Kirschen Ron Allan Goran Strbac, Contributions of Individual Generators to Loads and Flows, IEEE Transactions on Power Systems, Vol. 12, No. 1, February 1997. [10] Comparitive Analysis of Different Wheeling Charge Methodology, T.Palanisamy, R.Loganathan, K.Anand, International Journal of Engineering Research & Technology(IJERT), Vol.3, Issue 12, December-2014. [11] M. W. Mustafa and H. Shareef, A Comparison of Electric Power Tracing Methods Used in Deregulated Power Systems FirstInternational Power and Energy Conference PECon 2006, Putrajaya,Malaysia, November 28-29, 2006, pp. 156160 [12] Satyavir Singh, Power tracing in a deregulated power system: IEEE14-Bus case, International of computertechnology&applications,vol3(3),888894,june2012. 870