Voltage control in a distribution system using active power loads

Transcription

1 Scholars' Mine Masters Theses Student Research & Creative Works Fall 2016 Voltage control in a distribution system using active power loads Mounika Chava Follow this and additional works at: Part of the Electrical and Computer Engineering Commons Department: Recommended Citation Chava, Mounika, "Voltage control in a distribution system using active power loads" (2016). Masters Theses This Thesis - Open Access is brought to you for free and open access by Scholars' Mine. It has been accepted for inclusion in Masters Theses by an authorized administrator of Scholars' Mine. This work is protected by U. S. Copyright Law. Unauthorized use including reproduction for redistribution requires the permission of the copyright holder. For more information, please contact scholarsmine@mst.edu.

2 VOLTAGE CONTROL IN A DISTRIBUTION SYSTEM USING ACTIVE POWER LOADS by MOUNIKA CHAVA A THESIS Presented to the Faculty of the Graduate School of the MISSOURI UNIVERSITY OF SCIENCE AND TECHNOLOGY In Partial Fulfillment of the Requirements for the Degree MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 2016 Approved by Jhi-Young Joo, Advisor Mariesa L. Crow Pourya Shamsi

3 2016 Mounika Chava All rights reserved

4 iii ABSTRACT High penetration of renewable energy resources such as rooftop solar photovoltaic (PV) systems is exacerbating violations of the voltage limits in power distribution networks. Existing solutions to these over-/under-voltage issues include tap-changing transformers and shunt capacitors. However, tap changers only allow voltage changes in discrete steps and shunt capacitors cannot handle over-voltage situations. To overcome these problems, this thesis proposes voltage control in a distribution system by adjusting active power loads. The proposed solution can be used for both under- and over-voltage cases at any levels of voltage adjustment. First, the active power that needs to be adjusted at selected nodes to bring the voltages within acceptable limits is calculated for under and over-voltage cases using alternating current optimal power flow (ACOPF). ACOPF is solved for three different cases by varying the cost of active power adjustment and the marginal cost from the feeder supply, which is assumed to be the market price. Secondly, demand response on electric water heaters is implemented to achieve active power adjustments at selected nodes obtained from the ACOPF results. The problem is formulated over multiple time steps to minimize the energy costs of water heaters at a specific node subject to the dynamics of the water temperature, the energy consumption, the temperature constraints, and the voltage limits at the node. The IEEE 34-bus radial distribution system is used as a case study. By solving the same problem for different price settings, it is observed that the price structure affects the energy demand only in the under-voltage system, but not in the over-voltage system. The main contribution of this thesis is to show that by adjusting active power loads, voltages in a distribution system can be maintained within the permissible limits.

5 iv ACKNOWLEDGMENTS First, I would like to thank my advisor, Dr. Jhi Young Joo, for giving me the opportunity to work on this thesis. She has been a great support throughout my Master s program. I am grateful for the guidance she provided and for her valuable time. Second, I would also like to thank my thesis committee members, Dr. Mariesa L. Crow and Dr. Pourya Shamsi, for their guidance and taking time to review my thesis work. I would like to thank my lab mate, Maigha, for her timely help. Finally, I would like to thank my parents for their love and support, who stood by me and supported all my decisions.

6 v TABLE OF CONTENTS Page ABSTRACT... iii ACKNOWLEDGMENTS... iv LIST OF FIGURES... vi LIST OF TABLES... viii SECTION 1. INTRODUCTION BACKGROUND PROBLEM STATEMENT FORMULATION AND METHODOLOGY ACOPF FORMULATION DEMAND RESPONSE PROBLEM FORMULATION Electric Water Heater Mode DR Problem Formulation SIMULATION RESULTS ACTIVE POWER ADJUSTMENT Under-voltage System Over-voltage System DEMAND RESPONSE MODEL Under-voltage System Over-voltage System CONCLUSION BIBLIOGRAPHY VITA.40

7 vi LIST OF FIGURES Page Figure 2.1. Block diagram of water heater model... 6 Figure 3.1. IEEE 34 bus system Figure 3.2. Voltage profile of initial, under-voltage and over-voltage test system Figure 3.3. Under-voltage test system Figure 3.4. Voltage Profile of under-voltage test system Figure 3.5. Over-voltage test system Figure 3.6. Voltage Profile of over-voltage test system Figure 3.7. Electricity price settings 1, 2 and Figure 3.8. Hot water usage profiles of houses 1, 2 and 3 at node Figure 3.9. Energy consumption at node 824 a) price setting 1 b) price settings 2 and Figure Temperature profile at node 824 a) house 1 b) house Figure Temperature profile of house 3 at node Figure Hot water usage profiles of house 1 and 2 at node Figure Energy consumption at node 834 a) price setting 1 b) price settings 2 and Figure Temperature profile at node 834 a) house 1 b) house Figure Hot water usage profiles of houses 1 and 2 at node Figure Energy consumption at node 836 a) price setting 1 b) price settings 2 and Figure Temperature profile at node 836 a) house 1 b) house Figure Hot water usage profiles 1, 2 and Figure Energy consumption at node 824 a) price setting 1 b) price setting 2 and Figure Temperature profiles at node 824 a) house 1 b) house Figure Temperature profile of house 3 at node Figure Energy consumption at node 834 a) price setting 1 b) price settings 2 and Figure Temperature profiles at node 834 a) house 1 b) house

8 vii Figure Energy consumption at node 836 a) price setting 1 b) price setting 2 and Figure Temperature profiles at node 836 a) house 1 b) house

9 viii LIST OF TABLES Page Table 3.1. Total load of the system Table 3.2. Active and Reactive power limits for under and over-voltage system Table 3.3. ACOPF with adjusting active power only at the slack bus Table 3.4. ACOPF with slack bus and active power adjustment for Cases 1, 2 and Table 3.5. Cost analysis of correcting voltages in the under-voltage system Table 3.6. Active power adjustment at feeder and selected nodes Table 3.7. Comparison between cost for feeder adjustment and cost for adjustment at nodes 824, 834 and 836 for 3 cases Table 3.8. ACOPF with only slack bus generation Table 3.9. ACOPF with slack bus and active power adjustment for cases 1, 2 and Table Cost analysis of correcting voltages in the over-voltage system Table Active power adjustments in under & over-voltage system Table Simulation parameter setup Table Maximum energy consumption limit for under-voltage system Table Minimum energy consumption limits for over-voltage system... 30

10 1. INTRODUCTION 1.1 BACKGROUND In power systems, the voltage at each node is maintained between above and below 5% of the rated value. An electric power distribution system is generally designed as a radial system, which is arranged in a tree like structure with usually only one power supply source at the beginning of the feeder. Due to this structure, distribution systems are more prone to voltage issues than the transmission system. Under-voltages may be the result of 1) faults on the power system 2) capacitor bank switching off 3) high load on the system 4) loss of renewable distributed generation. Over voltages can occur because of 1) switching off high load 2) energizing a capacitor bank [2] 3) high penetration of distributed generation renewable sources particularly solar [3]. The voltage issues in a distribution system have become worse with increase in uncontrollable renewable distributed generation especially solar photovoltaic (PV). In the last few years, there has been a dramatic increase in renewable distributed generation. With the increasing PV, high PV penetration and low demand together can lead to overvoltage problems [18]. On the other hand, when there is peak load and low PV generation, voltage can drop below lower limit causing under-voltage issues [12]. There have been various studies on mitigating voltage violation problems in distribution system considering DG renewable energy resources. Reference [24] presents impact of distributed generation on voltages especially voltage sag problems. Reference [26] addressed effect of PV power variability on voltage regulation in distribution system considering 20% of PV penetration. The typical approach is to provide reactive power support from inverter units to tackle the voltage rise caused by high penetration of PV [8], [9], [10], [11] and [21]. References [19] and [21] proposed PV generation curtailment to prevent over voltages. Reference [20] suggested use of shunt reactors, shunt capacitors and transformer tap changer to prevent voltage instability. However, the tap changers only allow voltage change at discrete steps. Reactive power is adjusted at the point of PV installation to mitigate voltage rise problems [23]. The authors of [25] suggested use of onsite battery energy storage that is integrated with PV inverter to reduce the effects of

11 2 PV output variability. Reference [27] proposed use of PV over production in the LV feeder by shifting the peak loads i.e. demand side management and compared it against the PV energy curtailment. In this thesis, maintaining the bus voltages (under-voltages and over-voltages) within the acceptable limits is achieved using demand response (DR) with active power loads. The proposed solution can be used for both under and over-voltage issues and any level of voltage adjustment can be achieved. Demand response is defined as changes in electric usage by end use customers from their normal consumption patterns in response to changes in the price of electricity over time or incentive payments [5]. In this thesis, electric water heaters are used as controllable loads for demand response (DR). Water heaters contribute a large portion of household loads in the U.S. and have following advantages to implement demand response [6] Water heaters have relatively high consumption, which is up to 30% of household load compared to other home appliances [6] The heating element in an electric water heater is a resistor, which does not require reactive power support. Main aim is to adjust only active power but not reactive power 1.2 PROBLEM STATEMENT This section describes the objective of the thesis problem and steps involved in it. The main goal is to control the voltage in a distribution system by adjusting active power loads. The first step is to determine the active power adjustments needed to bring the voltages within limits. Over-voltage and under-voltage cases were setup by increasing and decreasing the system load respectively. These under-voltage and over-voltage systems are used for the rest of study and analysis. Few nodes has been selected to adjust the active power and assumed that the active power at those nodes can be adjusted. ACOPF problem (as shown in Section 2) is then formulated by modeling the load buses where the active power can be adjusted as generation buses. The ACOPF problem is modeled and solved over one time step in MATPOWER. The ACOPF results provide

12 3 information on the optimal feeder source generation and active power adjustment needed at the selected nodes to maintain voltages within acceptable limits. The second step is to achieve the active power adjustments at selected nodes obtained in step one using demand response on electric water heaters. A demand response model for electric water heaters is designed (as given in Section 2) over multiple time steps minimizing the energy costs subject to the dynamics of the water temperature, energy and the temperature limits. The DR model is solved over multiple time steps since the electric water heater cannot be on/off throughout the time and the water temperature has to be maintained within the temperature limits. The limits on energy consumption are derived from the active power adjustments obtained in the first step. DR model is formulated as a mixed integer linear programming and solved in MATLAB.

13 4 2. FORMULATION AND METHODOLOGY In this section, first ACOPF problem is formulated that determines the active power adjustments needed to bring the voltages within permissible limits. Next, a demand response model on electric water heaters is formulated to achieve the active power adjustments obtained from the ACOPF solution. 2.1 ACOPF FORMULATION The objective of ACOPF problem is to optimize active and reactive power dispatch subject to demand, transmission network, voltage constraints and active and reactive power generation output limits [7]. Active power adjustment by the loads in this work is equated as active power adjustment as generator. If the load decreases its consumption, it is equivalent as a generator producing active power at the node, and vice versa. Therefore, the load buses in the distribution system where the active power can be adjusted are treated as generation buses in the ACOPF problem. The formulation is given by minimize P gi, Q gi, V i, θ i Subject to N g C gi P gi i=1 n P gi P di = k=1 V i V k (G ik cos θ ik + B ik sin θ ik ) (1) n Q gi Q di = k=1 V i V k (G ik sin θ ik B ik cos θ ik ) (2) S ik,min S ik S ik,max (3) P gi,min P gi P gi,max for i = 1,, N g (4) Q gi,min Q gi Q gi,max for i = 1,, N g (5) V i,min V i V i,max (6) Where P gi Active power generation at bus i P di Active power demand at bus i Q gi Reactive power injection of generator at bus i

14 5 Q di Reactive power demand at bus i S ik - MVA flow on line ik V i Voltage magnitude at bus i C gi Cost function of generator at bus i G ik Conductance of line ik B ik Susceptance of line ik S ik,min, S ik,max Lower and upper MVA flow limits on line ik P gi,min, P gi,max - Lower and upper active power generation limits on line ik Q gi,min, Q gi,max - Lower and upper reactive power generation limits on line ik V i,min, V i,max Lower and upper limit of voltage magnitude at bus i Equations (1) and (2) are the power flow equations. (3) is related with the lower and upper flow on the lines. Equation 4 and 5 define active and reactive power limits of each generating unit. Equation 6 limits voltage at each bus. P gi limits for the load buses are obtained from the original P di at those nodes. ACOPF problem is modeled and solved in MATPOWER. 2.2 DEMAND RESPONSE PROBLEM FORMULATION To develop the demand response model, the electric water heater load model is first derived. Then an optimization problem is formulated for controlling water heaters with the water temperature dynamics as a constraint Electric Water Heater Model. This section presents electric water heater load model. Figure 2.1 shows the block diagram of the water heater model [14]. The water heater model parameters are classified into three parts: 1. The temperature profile that includes the ambient temperature, the inlet water temperature and the hot water temperature set points 2. The water heater characteristics including the heat resistance of the tank (R), the surface area of the tank (SA tank ) and the rated power 3. Hot water usage profile by the user

15 6 Figure 2.1. Block diagram of water heater model The water temperature in the tank is calculated as [14] T t+1 = T t(v tank V wd. Δt) V tank + T in. V wd. Δt V tank + 1gal 8.34lb. [P 3412 Btu t. η. kwh SA tank. (T t T amb ) R ]. Δt 60 min h. 1 V tank (7)

16 7 where T t V tank V wd Δt Water temperature in the tank ( ) in time slot t Volume of tank (gallons) Volume of hot water withdrawn (gallons per minute) Duration of each time slot (minutes) T in Inlet water temperature ( ) P t Electricity demand of water heating unit in time slot t (KW) SA tank Surface area of the tank (ft 2 ) T amb Ambient temperature ( ) R Heat resistance of tank (. ft 2. h/btu) η Efficiency factor DR Problem Formulation. The objective of demand response model is to minimize the total energy consumption subject to temperature dynamics of electric water heater, set point temperature and energy limits. As water temperature has to be maintained within bounds, the water heater cannot be on or off for the whole hour, hence the demand response problem has been solved for each 5-minute interval over an hour. It can be formulated as n 12 minimize C t E t Subject to i=1 t=1 T t+1 = T t(v tank V wd. Δt) V tank + T in. V wd. Δt V tank + 1gal 8.34lb. [P 3412 Btu t. η. kwh SA tank. (T t T amb ) R ]. Δt 60 min h. 1 V tank T t,min T t T t,max t (9) (8)

17 8 n i=1 E t E t,max t (10) n i=1 E t E t,min t (11) Where n is number of electric water heaters and t is the time interval. C t is electricity market price at each time step t. E t (P t /12) is energy consumption of electric water heater at each time step. At each time step t, demand for electricity of the water heater unit is calculated as P t = p r. W s (12) where p r W s Rated power of water heater (EWH) Water heater status, 0=OFF, 1=ON Equation 9 represents lower (T t,min ) and upper (T t,max ) bounds of water temperature. Equation 10 and 11 are Energy consumption limits (E t,min and E t,max ) which are determined based on the active power adjustments obtained in Section 3. Equation 12 defines the water heater status during the given interval, p r is the rated input power of electric water heater, which is constant and Ws is the solving variable that tells the status of water heater. High load and low generation together causes under-voltages in the system. Hence, the energy consumption in under-voltage system has to be decreased and in overvoltage system, energy consumption has to be increased. Maximum limit on energy consumption for under-voltage system and minimum energy consumption for overvoltage system has been imposed in problem (8)-(12). In under voltage case, maximum limit on energy consumption at each time step t is given as difference between energy consumption before DR and amount of energy consumption reduction that needs to be achieved at each node. In over voltage case, energy consumption lower bound is obtained by adding the energy consumption before DR and the amount of energy consumption that has to be increased.

18 9 For the DR optimization problem, expected hot water usage for each house and the electricity market price are given in 5-minute time intervals, with the initial water temperature known. The optimization problem is solved using a mixed integer linear programming with an equality constraint and minimum and maximum bounds. Energy consumption and water temperature at each time step t are the results of the optimization problem. To observe the price sensitivity of the loads, the same optimization problem is solved for different price settings.

19 10 3. SIMULATION RESULTS 3.1 ACTIVE POWER ADJUSTMENT The IEEE 34-node radial distribution network with generation at node 800 (slack bus) is used for the study as shown in Figure 3.1[13]. The test system is modeled after an actual feeder located in Arizona. The feeder s nominal voltage is 24.9 kv [13]. The single-phase balanced test system was modeled in MATPOWER ignoring voltage regulators between 814 and 850 and between 852 and 832 nodes. Figure 3.1. IEEE 34 bus system To study the test system further, the total load was adjusted (as shown in Table 3.1) randomly to setup a system with over voltages and under voltages. Over-voltage and under-voltage cases were derived by decreasing 13% of the load at each node and by increasing 16% of the load at each bus respectively.

20 11 Table 3.1. Total load of the system Initial system load Under -voltage system load Over-voltage system load P (MW) Q (Mvar) P (MW) Q (Mvar) P (MW) Q (Mvar) Figure 3.2 shows the voltage profile of the initial, under-voltage and over-voltage test system. From the graph, the voltage profile of the initial test system is within limits (0.95 p.u p.u.). In the under-voltage system, voltages at a few (16) nodes especially the nodes that are far from the source are below the lower limit (0.95 p.u.). In the overvoltage system, some of the node voltages are above (1.05 p.u.). After setting up the under-voltage and over voltage systems, next how much active power needs to be adjusted is determined to bring the voltages within acceptable limits of 0.95 p.u. to 1.05 p.u. The nodes that have lateral branches (nodes 824, 834, and 836) has been selected and assumed the active power at those nodes can be adjusted. V o l t a g e ( P U ) Node Initial System Under- voltage system Over-voltage system Figure 3.2. Voltage profile of initial, under-voltage and over-voltage test system

21 12 The ACOPF problem is solved for both under-voltage and over-voltage test systems, to calculate the active power adjustment to bring the voltages within limits. The problem is formulated in Section 2. Active power adjustment at nodes 824, 834 and 836 nodes is modeled as generation for solving ACOPF. The solution of the ACOPF problem gives the optimal slack bus generation and the active power that needs to be adjusted at nodes 824, 834 and 836. Maximum active power output limits for the active power adjustment at nodes 824, 834 and 836 is given as 10% of the load at those buses respectively (as shown in Table 3.2). Reactive power limits are given as zero since the objective is to adjust active power only. Bus Table 3.2. Active and Reactive power limits for under and over-voltage system Under-voltage system Over-voltage system Active Power (MW) Reactive Power (Mvar) Active Power (MW) Reactive Power (Mvar) Under-voltage System. Figure 3.3 represents the under-voltage test system and highlighted nodes represent the nodes whose voltages were below 0.95 p.u.

22 13 Figure 3.3. Under-voltage test system ACOPF problem is solved for the under-voltage system having only feeder supply, to obtain the optimal generation, or net active power production satisfying the voltage limits at each node. ACOPF results with adjusting active power only at the slack bus is given in the Table 3.3. Table 3.3. ACOPF with adjusting active power only at the slack bus Bus Slack bus Generation Active power (MW) Reactive power (Mvar) ACOPF problem is solved for the under-voltage system for below cases, to evaluate the economy of adjusting active power at the nodes (i.e., demand response), compared with the cost of adjusting active power from the feeder supply. Case 1: Cost of feeder supply ($10/MWh) < cost of active power adjustment ($50/MWh) Case 2: Cost of feeder supply ($50/MWh) > cost of active power adjustment ($10/MWh) Case 3: Cost of feeder supply ($10/MWh) = cost of active power adjustment ($10/MWh)

23 14 Table 3.4 shows the ACOPF results for all three cases. From the results, it is clear that in Case 1, slack bus is generating the total active and reactive power required by the system while the active power adjustment at nodes 824, 834 and 836 nodes is zero. This is because the cost for adjusting active power is higher than the cost of the slack bus generation. Bus Table 3.4. ACOPF with slack bus and active power adjustment for Cases 1, 2 and 3 Under voltage test system Case 1 cost of (feeder supply < Active power (MW) adjusting P) Reactive Power (Mvar) Case 2 cost of (feeder supply > Active power (MW) adjusting P) Reactive Power (Mvar) Case 3 cost of (feeder supply = adjusting P) Active power (MW) Reactive Power (Mvar) In Case 2, active power adjustment at nodes 824, 834, and 836 is using almost its full capacity (2,500 W at 824, 2,000 W at 834, and 1,800 W at 836). In total, 6,628 W needs to be adjusted to bring the voltage profile within the limits. Case 3 is similar to case 2 where active power adjustment needed at 824, 834 and 836 is its maximum active power limit (2500 W at 824, 2000 W at 834 and 1800 W at 836). To maintain voltage conditions in the under-voltage system, it requires 6,578 W active power adjustments at nodes 824, 834, and 836.

24 15 In all three cases, reactive power supplied by the slack (2.801 Mvar) is same. In Cases 2 and 3 active power adjustment needed is its maximum active power limit where as in Case 1 it is zero. This is because the cost for active power adjustment is less than or equal to slack bus generation in Cases 2 and 3. Figure 3.4. Voltage Profile of under-voltage test system shows voltage profile of initial system and improved voltage profile before adjusting and after adjusting active power at selected nodes. Voltage profile before and after adjusting active power is same. V o l t a g e ( P U ) Node Initial voltage Profile Voltage profile (only slack bus) Voltage profile (Slack bus & real power adjustment) Figure 3.4. Voltage Profile of under-voltage test system From Table 3.5, the total cost of generation remains the same before and after adjusting active power when the cost of slack bus generation is less than the cost of adjusting active power. However, in Cases 2 and 3 total cost of generation is low after adjusting active power compared to the total cost of generation before adjusting active power. The total optimal cost of generation in Case 2 is much higher than Cases 1 and 3. It is because most of the total generation is coming from slack bus and it is costly (five times the cost of active power adjustment) in Case 2. It is economical to adjust active

25 16 power to maintain the voltage condition of the system within permissible limits only if the cost of adjusting active power is less than or equal to the cost of slack bus generation. Case 1 Table 3.5. Cost analysis of correcting voltages in the under-voltage system The total cost of generation ($/hr) (cost for slack < P adjustment) Case 2 (cost for slack > P adjustment) Case 3 (cost for slack = P adjustment) Before adjusting active power After adjusting active power Over-voltage System. Figure 3.5 represents the over-voltage test system. The voltages at the highlighted nodes are above 1.05 p.u. The optimal generation of slack bus without any active power adjustment is given in Table 3.8. In an over-voltage system, voltage profile of the system can be improved by increasing the load. Producing negative active power is equivalent to increasing active power. Hence the selected load nodes are modeled as negative generation for solving ACOPF for the over voltage system in the study. Negative per-mwh cost with negative power production is a positive cost, which is adding to generation cost. Hence, it involves more cost to adjust active power using demand response compared to cost of active power adjustment at the slack bus. Table 3.6 show the ACOPF results with only feeder supply and feeder supply along with adjusting active power.

26 17 Figure 3.5. Over-voltage test system Table 3.6. Active power adjustment at feeder and selected nodes Node Only feeder adjustment (MW) Adjustment at slack and nodes (MW) ACOPF is solved for three different cases Case 1- Cost of feeder supply (10 $/MWh) < cost for active power adjustment (-50 $/MWh) Case 2- Cost of feeder supply (50 $/MWh) > cost for active power adjustment (-10 $/MWh) Case 3- Cost of feeder supply (10 $/MWh) = cost for active power adjustment (-10 $/MWh) From the Table 3.7, for all the three cases, cost of active power adjustment at the slack bus is cheaper than the active power adjustment at the nodes 824, 834 and 836. From simulations, assuming a positive per-wh cost for active power consumption (or negative active power production) did not turn out to be an economical option. Therefore,

27 18 in the over-voltage case, demand response can be economically used only when a negative cost is assigned to active power consumption (or negative active power production). In other words, a per-wh net benefit for active power consumption, which is equivalent to a negative cost of producing active power is considered. Net benefit of adjusting active power is defined as revenue minus costs. The costs include maintenance, losses and costs for implementing demand response, such as meter/controller installation. Hence, cases 1, 2 and 3 are redefined as Case 1 Cost of feeder supply (10 $/MWh) < net benefit of adjusting active power (50 $/MWh) Case 2 Cost of feeder supply (50 $/MWh) > net benefit of adjusting active power (10 $/MWh) Case 3 Cost of feeder supply (10 $/MWh) = net benefit of adjusting active power (10 $/MWh) Table 3.7. Comparison between cost for feeder adjustment and cost for adjustment at nodes 824, 834 and 836 for 3 cases Cost for Cost for only feeder adjustment adjustment at nodes 824,834 and 836 Case 1: (cost of feeder supply (10 $/MWh) < cost for active power adjustment (-50 $/MWh) Case 2: (cost of feeder supply (50 $/MWh) > cost for active power adjustment (-10 $/MWh) Case 3: (cost of feeder supply (10 $/MWh) = cost for active power adjustment (-10 $/MWh)

28 19 The ACOPF results for the over-voltage case with only slack bus generation are shown in Table 3.9. From Table 3.9, Case 2 is same as system without any active power adjustment. In Cases 1 and 3, active power adjustments needed at nodes 824, 834 and 836 are equal to its minimum active power limits.in total, 5,046 W needs to be adjusted in Cases 1 and 3 respectively, to bring the voltage profile within the limits. Table 3.8. ACOPF with only slack bus generation Bus Optimal slack bus Generation Active power (MW) Reactive power (Mvar) Bus Table 3.9. ACOPF with slack bus and active power adjustment for cases 1, 2 and 3 Over -voltage test system Case 1 (cost of feeder supply < net benefit of adjusting Active power (MW) P) Reactive Power (Mvar) Case 2 (cost of feeder supply > net benefit of adjusting Active power (MW) P) Reactive Power (Mvar) Case 3 (cost of feeder supply = net benefit of adjusting P) Active power (MW) Reactive Power (Mvar)

29 20 Figure 3.6 shows the voltage profile of initial system, which has voltage violations, and improved voltage profile after adjusting active power. The total cost of generation in cases 2 and 3 before and after adjusting active power remains same. In case 1, the total cost of generation after adjusting active power is less than the cost before adjusting active power (as shown in Table 3.10). V o l t a g e ( P U ) Node Initial volatge profile Voltage profile(slack bus only) Voltage profile ( Slack bus & real power adjustment) Figure 3.6. Voltage Profile of over-voltage test system Table 3.11 shows the active power adjustments needed at the nodes 824, 834 and 836 to bring the voltages within limits (0.95 p.u p.u.). In the under-voltage system, the load by the amount of active power is decreased as shown in Table 3.11 at each node whereas the active power at each node is increased in the over-voltage system.

30 21 Table Cost analysis of correcting voltages in the over-voltage system The total cost of generation ($/hr) Before adjusting P After adjusting P Case 1 (cost for slack < net benefit of adjusting P) Case 2 (cost for slack > net benefit of adjusting P) Case 3 (cost for slack = net benefit of adjusting P) Table Active power adjustments in under & over-voltage system Node Active power adjustment needed(w) Under-voltage system Over-voltage system DEMAND RESPONSE MODEL In this section, demand response (DR) on electric water heater is implemented to achieve the active power adjustments at nodes 824, 834 and 836 as calculated in Section 3.1. First, the number of water heaters that needs to be controlled is determined. Typical individual household peak load is 6 kw [15]. Considering the household peak load and total load at nodes 824, 834 and 836, 3 houses are assigned at node 824, 2 houses at node 834 and 836. It is assumed that each household has an electric water heater. Demand response problem is simulated in MATLAB over an hour in 5-minute intervals.

31 22 The lower (T t,min ) and upper (T t,max ) bounds of the hot water temperature are considered as 110 and 130 [17]. Energy consumption limits (E t,min and E t,max ) is determined based on the active power adjustments needed at nodes 824, 834 and 836 (as shown in Table Active power adjustments in under & over-voltage system). Table 3.12 shows the parameters setup [14] needed to run the DR optimization problem. Table Simulation parameter setup Attribute Value Inlet water temperature (T in ) 60 ( ) Ambient temperature (T amb ) 60 ( ) Tank volume ( V tank ) 40 (gallons) Heat resistance of tank (R) 20 (. ft 2. h/btu) Surface area of tank (SA tank ) 30.8 (ft 2 ) Rated power (P) 4.5 KW In order to obtain the price sensitivity of the electric water heater load, the same optimization problem is solved for three different electricity price profiles (as shown in Figure 3.7 [16] Under-voltage System. It is assumed that all the houses at nodes at 824, 834 and 836 are at their peak during the hour. Each house at a node has different water usage profiles and the initial water temperature is assumed 118. Three different water usage profiles are used for the analysis and are given as Figure 3.8. First the default case is simulated i.e., before DR when there are no constraints on the energy consumption of water heaters. The water heater status is determined according to the following rules: the water heater status is OFF when the tank water temperature is within the limits. If the tank temperature drops below the lower temperature bound, the heating coils start working at its rated power until it reaches the upper limit.

32 Electricity Price ( /kwh) Price setting 1 Price setting 2 Price setting Figure 3.7. Electricity price settings 1, 2 and 3 To solve the DR model (formulated in Section 2), E t,max has to be determined at each node as they vary with the active power adjustment needed at each node. E t,max = consumption before DR active adjustment needed at the node 12 (13) Hence, E t,max at each node is given as in Table The inputs to DR model and simulated results before and after DR at each node are arranged in following order 1. Hot water usage [17] 2. Total energy consumption 3. Water temperature profile of each electric water heater Table Maximum energy consumption limit for under-voltage system Node E t,max (kwh)

33 Hot water consumption (gal/min) 24 Hot water consumption (gal/min), house House 1 House 2 House 3 Figure 3.8. Hot water usage profiles of houses 1, 2 and 3 at node 824 Figures 3.9, 3.10 and 3.11 are the DR results at node 824. From Figure 3.9, total energy consumption over an hour at node 824 before and after DR is kWh and kwh respectively. At node 824, as the upper limit on energy consumption is kWh, during any interval, maximum number of water heaters that are allowed to operate is two. From Figure 3.9, before DR three water heaters are ON during 6 th,7 th,8 th,9 th,10 th,11 th and 12 th intervals but after DR, maximum number of water heaters ON in any of the intervals is two. Before DR when water temperature drops below the lower bound (110 ), the heating element starts working and tries to reach the upper bound (130 ). The demand response model does not require the water temperature to reach its maximum but it maintains the temperature between 110 and 130 minimizing the energy consumption. From Figure 3.10 and 3.11 before DR, water temperature dropped below 110 before

34 Water Temperature ( ) Water Temperature ( ) Energy Consumption (kwh) Energy Consumption (kwh) 25 and during the fifth interval. From then, the water heater is ON until the last interval to reach the maximum temperature limit. Energy consumption (kwh) Energy consumption (kwh) Before DR After DR Price setting 2 Price setting 3 Figure 3.9. Energy consumption at node 824 a) price setting 1 b) price settings 2 and 3 Temperature Profile, House 1 Temperature Profile, House Before DR After DR Before DR After DR Figure Temperature profile at node 824 a) house 1 b) house 2

35 Hot water consumption (gal/min) Water Temperature ( ) 26 Temperature Profile, House Before DR After DR Figure Temperature profile of house 3 at node 824 By solving the same optimization problem for different price settings, obtained different energy usage at each time step are obtained. Energy consumption for price settings 2 and 3 is obtained and shown in Figure 3.9. From Figure 3.9, total energy consumption is same for all price settings but the energy usage over a particular time interval is different. Figure 3.12 shows the hot water usage profile of house 1 and 2 at node 834. Hot water consumption (gal/min) House 1 House Figure Hot water usage profiles of house 1 and 2 at node 834

36 Energy consumption (kwh) Energy consumption (kwh) 27 By solving the DR model at node 834, obtained results as shown in Figure 3.13 and From Figure 3.13, total energy consumption at node 834 for price setting 1 before and after DR is 6.75 kwh and kwh respectively. Energy consumption (kwh) Energy consumption (kwh) Before DR After DR Price setting 2 Price setting 3 Figure Energy consumption at node 834 a) price setting 1 b) price settings 2 and Temperature ( ) Temperature ( ) Before DR After DR Before DR After DR Figure Temperature profile at node 834 a) house 1 b) house 2

37 Hot water usage (gal/min) 28 Before DR, both the water heaters are ON from 5 th interval to 12 th interval but after DR due to limit on energy consumption, any one of the water heaters is allowed to operate at time step t. From Figure 3.13, it is clear that only one water heater is on while the other water heaters are off at a given time interval. Figure 3.14 represents water temperature profile for each house at node 834 before and after DR. From Figure 3.14, before DR water temperature dropped below 110 at second and fifth intervals for electric water heater one and two respectively. From then until the last interval, water heater is ON to bring the water temperature to its maximum. After DR, the water temperatures are maintained between 110 and 130. By solving the same DR model at node 834 for different price settings 2 and 3, energy consumption of water heater is obtained as shown in Figure Total energy consumption over an hour is same for price setting one, two and three but energy usage in a specific time interval varies with the price structure. Hot water usage at node 834 is given as in Figure Hot water usage(gal/min) House 1 House 2 Figure Hot water usage profiles of houses 1 and 2 at node 836 Similar to the results at nodes 824 and 834, energy consumption and temperature profile of water heaters 1 anad 2 at node 836 is given in Figures 3.16 and From

38 Temperature ( ) Temperature ( ) Energy consumption (kwh) Energy consumption (kwh) 29 Figure 3.16, total energy consumption at node 836 before and after DR for price setting 1 is 6.375kWH and 3.37 kwh respectively. Energy consumption at node 836 for price settings 2 and 3 is as in Figure Before DR After DR Price setting 2 Price setting 3 Figure Energy consumption at node 836 a) price setting 1 b) price settings 2 and Before DR After DR Before DR After DR Figure Temperature profile at node 836 a) house 1 b) house 2

39 30 Total energy consumption remains same for all price settings but the energy usage over a time step t varies with price structure. Figure 3.17 represents water temperature profiles water heater 1 and 2. Before DR, water temperature dropped below 110 during 5 th and 4 th time intervals. From then the water heaters are on until the the last time interval trying to reach the maximum limit. After DR, temperatures are maintained between limits ( ) Over-voltage System. Over-voltages are due to high generation or less load. The load is increased by the amount of active power as calculated in Section 3.1. In over-voltage case, limits on the minimum energy demand is imposed. The initial temperature is assumed 115. Hence, E t,min at each node is calculated using below equation (given as in Table 3.14). E t,min = consumption before DR + active adjustment needed at the node 12 (14) Table Minimum energy consumption limits for over-voltage system Node E t,min (kwh) The hot water usage demand is very low in the over-voltage system compared to the demand in the under-voltage system. Three different hot water usage profiles are defined as shown in Figure As the hot water usage is very low, the energy consumption of water heaters at node 824 before DR is 0 and the water temperature is between 110 and 130 during all the intervals. In the DR model as limits on minimum energy consumption ( kwh) has been imposed throughout the whole period, this means that any one of the three water heaters has to be on at each interval. Figure 3.19 shows the energy consumption of water heaters at node 824 before and after DR. From the Figure 3.19, energy consumption before DR is zero because of low water usage. After

40 Hot water consumption (gal/min) 31 DR, energy consumption is set to its minimum limit i.e. any one of the three water heaters are on during any given interval. Hot water consumption (gal/min) Profile 1 Profile 2 Profile 3 Figure Hot water usage profiles 1, 2 and 3 The same optimization problem is solved for different price settings to observe the price sensitivity. Unlike the under-voltage system case, price does not have any effect on the energy consumption in the over-voltage system. This is because of the constraint on the minimum energy consumption of water heaters. The result of DR model is setting energy usage output at the minimum limit (shown in Figure 3.19). From Figure 3.20 and 3.21 before DR, the water temperature profile during all the intervals is between 110 and 130, it did not drop below 110 hence water heater status is off and energy consumption is zero. After DR, since the water heaters, one and two are on for 15 minutes and water heater 3 is on for 30-minutes.Hnece, water temperature of each water heater raised gradually compared to the default case.

41 Temperature ( F) Temperature ( F) Energy consumption (kwh) Energy consumption (kwh) 32 Energy consumption (kwh) Energy consumption (kwh) Before DR After DR Price setting 2 Price setting 3 Figure Energy consumption at node 824 a) price setting 1 b) price setting 2 and 3 Temperature profile, house 1 Temperature profile, house Before DR After DR Before DR After DR Figure Temperature profiles at node 824 a) house 1 b) house 2

42 Energy consumption (kwh) Energy consumption (kwh) Temperature ( F) 33 Temperature profile, house Before DR After DR Figure Temperature profile of house 3 at node 824 Figure 3.22 and 3.24 are the results of DR at node 834. Hot water usage profile 1 and 2, price settings 1 and 2 are given as inputs to DR model. From Figure 3.22, before DR the energy usage is zero. Since the water temperature did not drop below 110, water heater is off throughout the period. After DR, any one of the waters is on during all the intervals because of the limit on minimum energy consumption. Hence the water temperature of water heater 1 and 2 increased gradually (shown in Figure 3.23). Energy consumption (kwh) Energy consumption (kwh) Before DR After DR Price setting 3 Price setting 2 Figure Energy consumption at node 834 a) price setting 1 b) price settings 2 and 3

43 Energy consumption (kwh) Energy consumption (kwh) Temperature ( F) Temperature ( F) 34 Temperature profile, house 1 Temperature profile, house Before DR After DR Before DR After DR Figure Temperature profiles at node 834 a) house 1 b) house 2 Water heater is not required to be on because for the given low water usage, water temperature is within limits. Since the objective of DR model is to minimize energy costs satisfying energy and the temperature constraints, it is setting the energy usage at its minimum limit. The DR results at node 836 are shown in Figure 3.24 and Similar to the energy consumption at node 824 and 834, from Figure 3.24 energy usage at node 836 before DR is 0 and it is at its minimum limit after DR. Energy consumption (kwh) Energy consumption (kwh) Before DR After DR Price setting 3 Price setting 2 Figure Energy consumption at node 836 a) price setting 1 b) price setting 2 and 3

44 Temperature ( F) Temperature ( F) 35 Before DR, the water temperature is between the limits and it did not drop below 110. The water temperature rose gradually after DR. Energy usage for different price settings is same (shown in Figure 3.24) Temperature profile, house Temperature profile, house Before DR After DR Before DR After DR Figure Temperature profiles at node 836 a) house 1 b) house 2

45 36 4. CONCLUSION This thesis provides an approach to maintain voltages within the limits in a distribution system with active power adjustment. The ACOPF problem is solved to calculate the active power adjustments at each node to bring the violated voltages within the limits. The problem is solved for three different cases where the marginal cost of feeder supply and the cost of active power adjustment were varied. The solution to this problem i.e., the active power adjustments at eligible nodes was realized by adjusting energy consumption of electric water heaters. It can be concluded that voltages in a distribution system can be maintained within limits by adjusting active power using demand response. This is especially economical when the cost of adjusting active power with loads is lower than or equal to the cost form the feeder supply in under-voltage system. In over-voltage case, demand response is economical when cost of feeder supply is less than net benefit of adjusting P. There is much future work ahead where the assumption on nodes that has the active power adjustment capability that were made to show the proof of concept can be relaxed and systematic approach can be defined. Assumptions that were made in solving demand response on electric water heater temperature such as the same initial temperature, inlet and ambient temperature for all the houses can be relaxed and more variability can be included..

46 37 BIBLIOGRAPHY [1] Akash T. Davda, Bhupendra. R. Parekh System Impact Analysis of Renewable Distributed Generation on an Existing Radial Distribution Network 2012 IEEE Electrical Power and Energy Conference. [2] Surya Santoso Fundamentals of Electric Power Quality [3] E. Liu and J. Bebic Distribution System Voltage Performance Analysis for High- Penetration Photovoltaics [online]. [4] W. El-Khattam and M. M. A. Salama, Distributed Generation technologies, definitions and benefits in Electric Power Systems Research., [5] W. El-Khattam and M. M. A. Salama, Benefits of demand response in electricity markets and recommendations for achieving them [online]. Available: Demand_Response_in_Electricity_Markets_and_Recommendations_for_Achievi ng_them_report_to_congress.pdf. [6] Ruisheng Diao, Shuai Lu, Marcelo Elizondo, Ebony Mayhorn, Yu Zhang, Nader Samaan Electric Water Heater Modeling and Control Strategies for Demand Response 2012 IEEE Power and Energy Society General Meeting. [7] Hui Zhang, StuMIEEE, Vijay Vittal, Gerald T.Heydt, and Jaime Quintero A relaxed AC Optimal Power Flow Model Based on a Taylor Series. IEEE Innovative Smart Gris Technologies-Asia (ISGT Asia) November [8] Roozbeh Kabiri, Donald Grahame Holmes, Brendan P. McGrath, Lasantha Gunaruwan Meegahapola LV Grid Voltage Regulation Using Transformer Electronic Tap Changing, With PV Inverter Reactive Power Injection IEEE journal of emerging and selected topics in power electronics Vol.3. No.4,December [9] Mohamed M. Aly, Mamdouh Abdel-Akher, Zakaria Ziadi and Tomonobo Senjyu Voltage Stability Assessment of Photovoltaic Energy Systems with Voltage Control Capabilities Renewable Energy Research and Applications (ICRERA), 2012 International Conference Nov [10] Tzung-Lin Lee, Shih-Sian Yang Shang-Hung Hu Design of Decentralized Voltage Control for PV Inverters to Mitigate Voltage Rise in Distribution Power System without Communication The 2014 International Power Electronics Conference. [11] Juan C. Vasquez, R. A. Mastromauro, Josep M. Guerrero, Marco Liserre Voltage Support Provided by a Droop-Controlled Multifunctional Inverter IEEE Transactions on industrial electronics Vol 56, No 11, November 2009.