Chapter 25
25.1 Load Classification 25-1
25.2 Modeling Applications 25-2
25.3 Load Modeling Concepts and Approaches 25-3
25.4 Load Characteristics and Models 25-3
25.5 Static Load Characteristics 25-6
Exponential Models • Polynomial Models • Combined Exponential and Polynomial Models • Comparison of Exponential and Polynomial Models • Devices Contributing to Modeling Difficulties
25.6 Load Window Modeling 25-10
References 25-11
Raymond R. Shoults
University of Texas at Arlington
Larry D. Swift
University of Texas at Arlington
The physical structure of most power systems consists of generation facilities feeding bulk power into a high-voltage bulk transmission network that in turn serves any number of distribution substations. A typical distribution substation will serve from 1 to as many as 10 feeder circuits. A typical feeder circuit may serve numerous loads of all types. A light to medium industrial customer may take service from the distribution feeder circuit primary, while a large industrial load complex may take service directly from the bulk transmission system. All other customers, including residential and commercial, are typically served from the secondary of distribution transformers that are in turn connected to a distribution feeder circuit. Figure 25.1 illustrates a representative portion of a typical configuration.
The most common classification of electrical loads follows the billing categories used by the utility companies. This classification includes residential, commercial, industrial, and other. Residential customers are domestic users, whereas commercial and industrial customers are obviously business and industrial users. Other customer classifications include municipalities, state and federal government agencies, electric cooperatives, educational institutions, etc.
Although these load classes are commonly used, they are often inadequately defined for certain types of power system studies. For example, some utilities meter apartments as individual residential customers, while others meter the entire apartment complex as a commercial customer. Thus, the common classifications overlap in the sense that characteristics of customers in one class are not unique to that class. For this reason some utilities define further subdivisions of the common classes.
A useful approach to classification of loads is by breaking down the broader classes into individual load components. This process may altogether eliminate the distinction of certain of the broader classes, but it is a tried and proven technique for many applications. The components of a particular load, be it residential, commercial, or industrial, are individually defined and modeled. These load components as a whole constitute the composite load and can be defined as a “load window.”
It is helpful to understand the applications of load modeling before discussing particular load characteristics. The applications are divided into two broad categories: static (“snap-shot” with respect to time) and dynamic (time varying). Static models are based on the steady-state method of representation in power flow networks. Thus, static load models represent load as a function of voltage magnitude. Dynamic models, on the other hand, involve an alternating solution sequence between a time-domain solution of the differential equations describing electromechanical behavior and a steady-state power flow solution based on the method of phasors. One of the important outcomes from the solution of dynamic models is the time variation of frequency. Therefore, it is altogether appropriate to include a component in the static load model that represents variation of load with frequency. The lists below include applications outside of Distribution Systems but are included because load modeling at the distribution level is the fundamental starting point.
Static applications: Models that incorporate only the voltage-dependent characteristic include the following.
Dynamic applications: Models that incorporate both the voltage- and frequency-dependent characteristics include the following.
Strictly power-flow based solutions utilize load models that include only voltage dependency characteristics. Both voltage and frequency dependency characteristics can be incorporated in load modeling for those hybrid methods that alternate between a time-domain solution and a power flow solution, such as found in Transient Stability and Dynamic Stability Analysis Programs, and Operator Training Simulators (EPRI User’s Manual, 1992; EPRI Final Report EL-5003, 1987; Kundur, 1994).
Load modeling in this section is confined to static representation of voltage and frequency dependencies. The effects of rotational inertia (electromechanical dynamics) for large rotating machines are discussed in Chapters 13 and 17 of Power System Stability and Control. Static models are justified on the basis that the transient time response of most composite loads to voltage and frequency changes is fast enough so that a steady-state response is reached very quickly.
There are essentially two approaches to load modeling: component based and measurement based. Load modeling research over the years has included both approaches (EPRI, 1981, 1984, 1985). Of the two, the component-based approach lends itself more readily to model generalization. It is generally easier to control test procedures and apply wide variations in test voltage and frequency on individual components.
The component-based approach is a “bottom-up” approach in that the different load component types comprising load are identified. Each load component type is tested to determine the relationship between real and reactive power requirements versus applied voltage and frequency. A load model, typically in polynomial or exponential form, is then developed from the respective test data. The range of validity of each model is directly related to the range over which the component was tested. For convenience, the load model is expressed on a per-unit basis (i.e., normalized with respect to rated power, rated voltage, rated frequency, rated torque if applicable, and base temperature if applicable). A composite load is approximated by combining appropriate load model types in certain proportions based on load survey information. The resulting composition is referred to as a “load window.”
The measurement approach is a “top-down” approach in that measurements are taken at either a substation level, feeder level, some load aggregation point along a feeder, or at some individual load point. Variation of frequency for this type of measurement is not usually performed unless special test arrangements can be made. Voltage is varied using a suitable means and the measured real and reactive power consumption recorded. Statistical methods are then used to determine load models. A load survey may be necessary to classify the models derived in this manner. The range of validity for this approach is directly related to the realistic range over which the tests can be conducted without damage to customers’ equipment. Both the component and measurement methods were used in the EPRI research projects EL-2036 (1981) and EL-3591 (1984, 1985). The component test method was used to characterize a number of individual load components that were in turn used in simulation studies. The measurement method was applied to an aggregate of actual loads along a portion of a feeder to verify and validate the component method.
Static load models for a number of typical load components appear in Table 25.1 and Table 25.2 (EPRI, 1984, 1985). The models for each component category were derived by computing a weighted composite from test results of two or more units per category. These component models express per-unit real power and reactive power as a function of per-unit incremental voltage and/or incremental temperature and/or per-unit incremental torque. The incremental form used and the corresponding definition of variables are outlined below:
Static Models of Typical Load Components—AC, Heat Pump, and Appliances
Load Component |
Static Component Model |
1-ϕ central air conditioner |
P = 1.0 + 0.4311 * Δ V + 0.9507 * Δ T + 2.070 * Δ V2 + 2.388 * Δ T2 − 0.900 * Δ V * Δ T |
Q = 0.3152 + 0.6636 * Δ V + 0.543 * Δ V2 + 5.422 * Δ V3 + 0.839 * Δ T2 − 1.455 * Δ V * Δ T |
|
3-ϕ central air conditioner |
P = 1.0 + 0.2693 * Δ V + 0.4879 * Δ T + 1.005 * Δ V2 − 0.188 * Δ T2 − 0.154 * Δ V * Δ T |
Q = 0.6957 + 2.3717 * Δ V + 0.0585 * Δ T + 5.81 * Δ V2 + 0.199 * Δ T2 − 0.597 * Δ V * Δ T |
|
Room air conditioner (115 V rating) |
P = 1.0 + 0.2876 * Δ V + 0.6876 * Δ T + 1.241 * Δ V2 + 0.089 * Δ T2 − 0.558 * Δ V * Δ T |
Q = 0.1485 + 0.3709 * Δ V + 1.5773 * Δ T + 1.286 * Δ V2 + 0.266 * Δ T2 − 0.438 * Δ V * Δ T |
|
Room air conditioner (208/230 V rating) |
P = 1.0 + 0.5953 * Δ V + 0.5601 * Δ T + 2.021 * Δ V2 + 0.145 * Δ T2 − 0.491 * Δ V * Δ T |
Q = 0.4968 + 2.4456 * Δ V + 0.0737 * Δ T + 8.604 * Δ V2 − 0.125 * Δ T2 − 1.293 * Δ V * Δ T |
|
3-ϕ heat pump (heating mode) |
P = 1.0 + 0.4539 * Δ V + 0.2860 * Δ T + 1.314 * Δ V2 − 0.024 * Δ V * Δ T |
Q = 0.9399 + 3.013 * Δ V − 0.1501 * Δ T + 7.460 * Δ V2 − 0.312 * Δ T2 − 0.216 * Δ V * Δ T |
|
3-ϕ heat pump (cooling mode) |
P = 1.0 + 0.2333 * Δ V + 0.5915 * Δ T + 1.362 * Δ V2 + 0.075 * Δ T2 − 0.093 * Δ V * Δ T |
Q = 0.8456 + 2.3404 * Δ V − 0.1806 * Δ T + 6.896 * Δ V2 + 0.029 * Δ T2 − 0.836 * Δ V * Δ T |
|
1-ϕ heat pump (heating mode) |
P = 1.0 + 0.3953 * Δ V + 0.3563 * Δ T + 1.679 * Δ V2 + 0.083 * Δ V * Δ T |
Q = 0.3427 + 1.9522 * Δ V − 0.0958 * Δ T + 6.458 * Δ V2 − 0.225 * Δ T2 − 0.246 * Δ V * Δ T |
|
1-ϕ heat pump (cooling mode) |
P = 1.0 + 0.3630 * Δ V + 0.7673 * Δ T + 2.101 * Δ V2 + 0.122 * Δ T2 − 0.759 * Δ V * Δ T |
Q = 0.3605 + 1.6873 * Δ V + 0.2175 * Δ T + 10.055 * Δ V2 − 0.170 * Δ T2 − 1.642 * Δ V * Δ T |
|
Refrigerator |
P = 1.0 + 1.3958 * Δ V + 9.881 * Δ V2 + 84.72 * Δ V3 + 293 * Δ V4 |
Q = 1.2507 + 4.387 * Δ V + 23.801 * Δ V2 + 1540 * Δ V3 + 555 * Δ V4 |
|
Freezer |
P = 1.0 + 1.3286 * Δ V + 12.616 * Δ V2 + 133.6 * Δ V3 + 380 * Δ V4 |
Q = 1.3810 + 4.6702 * Δ V + 27.276 * Δ V2 + 293.0 * Δ V3 + 995 * Δ V4 |
|
Washing machine |
P = 1.0 + 1.2786 * Δ V + 3.099 * Δ V2 + 5.939 * Δ V3 |
Q = 1.6388 + 4.5733 * Δ V + 12.948 * Δ V2 + 55.677 * Δ V3 |
|
Clothes dryer |
P = 1.0 − 0.1968 * Δ V − 3.6372 * Δ V2 − 28.32 * Δ V3 |
Q = 0.209 + 0.5180 * Δ V + 0.363 * Δ V2 − 4.7574 * Δ V3 |
|
Television |
P = 1.0 + 1.2471 * Δ V + 0.562 * Δ V2 |
Q = 0.2431 + 0.9830 * Δ V + 1.647 * Δ V2 |
|
Fluorescent lamp |
P = 1.0 + 0.6534 * Δ V − 1.65 * Δ V2 |
Q = −0.1535 − 0.0403 * Δ V + 2.734 * Δ V2 |
|
Mercury vapor lamp |
P = 1.0 + 0.1309 * Δ V + 0.504 * Δ V2 |
Q = −0.2524 + 2.3329 * Δ V + 7.811 * Δ V2 |
|
Sodium vapor lamp |
P = 1.0 + 0.3409 * Δ V −2.389 * Δ V2 |
Q = 0.060 + 2.2173 * Δ V + 7.620 * Δ V2 |
|
Incandescent |
P = 1.0 + 1.5209 * Δ V + 0.223 * Δ V2 |
Q = 0.0 |
|
Range with oven |
P = 1.0 + 2.1018 * Δ V + 5.876 * Δ V2 + 1.236 * Δ V3 |
Q = 0.0 |
|
Microwave oven |
P = 1.0 + 0.0974 * Δ V + 2.071 * Δ V2 |
Q = 0.2039 + 1.3130 * Δ V + 8.738 * Δ V2 |
|
Water heater |
P = 1.0 + 0.3769 * Δ V + 2.003 * Δ V2 |
Q = 0.0 |
|
Resistance heating |
P = 1.0 + 2 * Δ V + Δ V2 |
Q = 0.0 |
Static Models of Typical Load Components—Transformers and Induction Motors
Load Component |
Static Component Model |
Transformer |
|
Core loss model |
|
|
|
Where V is voltage magnitude in per unit |
|
1-ϕ motor |
P = 1.0 + 0.5179 * Δ V + 0.9122 * Δτ + 3.721 * Δ V2 + 0.350 * Δ τ 2 − 1.326 * Δ V * Δτ |
Constant torque |
Q = 0.9853 + 2.7796 * Δ V + 0.0859 * Δ τ +7.368 * Δ V2 + 0.218 * Δτ 2 − 1.799 * Δ V * Δτ |
3-ϕ motor (1–10 HP) |
P = 1.0 + 0.2250 * Δ V + 0.9281 * Δτ + 0.970 * Δ V2 + 0. 086 * Δτ 2 − 0.329 * Δ V * Δτ |
Const. torque |
Q = 0.7810 + 2.3532 * Δ V + 0.1023 * Δτ − 5.951 * Δ V2 + 0.446 * Δτ 2 − 1.48 * Δ V * Δτ |
3-ϕ motor (10 HP/above) |
P = 1.0 + 0.0199 * Δ V + 1.0463 * Δτ + 0.341 * Δ V2 + 0.116 * Δτ 2 − 0.457 * Δ V * Δτ |
Const. torque |
Q = 0.6577 + 1.2078 * Δ V + 0.3391 * Δτ + 4.097 * Δ V2 + 0.289Δτ 2 − 1.477 * Δ V * Δτ |
1-ϕ motor |
P = 1.0 + 0.7101 * Δ V + 0.9073 * Δτ + 2.13 * Δ V2 + 0.245 * Δτ 2 − 0.310 * Δ V * Δ τ |
Variable torque |
Q = 0.9727 + 2.7621 * Δ V + 0.077 * Δτ + 6.432 * Δ V2 + 0.174 * Δτ 2 − 1.412 * Δ V * Δτ |
3-ϕ motor (1–10 HP) |
P = 1.0 + 0.3122 * Δ V + 0.9286 * Δτ + 0.489 * Δ V2 + 0.081 * Δτ 2 − 0.079 * Δ V * Δτ |
Variable torque |
Q = 0.7785 + 2.3648 * Δ V + 0.1025 * Δτ + 5.706 * Δ V2 + 0.13 * Δτ 2 − 1.00 * Δ V * Δτ |
3-ϕ motor (10 HP and above) |
P = 1.0 + 0.1628 * Δ V + 1.0514 * Δτ ∠0.099 * Δ V2 + 0.107 * Δτ 2 + 0.061 * Δ V * Δτ |
Variable torque |
Q = 0.6569 + 1.2467 * Δ V + 0.3354 * Δτ + 3.685 * Δ V2 + 0.258 * Δτ 2 − 1.235 * Δ V * Δτ |
ΔV = Vact − 1.0 (incremental voltage in per unit)
ΔT = Tact − 95°F (incremental temperature for Air Conditioner model), or
= Tact − 47°F (incremental temperature for Heat Pump model)
Δτ = τact − τrated (incremental motor torque, per unit)
If ambient temperature is known, it can be used in the applicable models. If it is not known, the temperature difference, ΔT, can be set to zero. Likewise, if motor load torque is known, it can be used in the applicable models. If it is not known, the torque difference, Δτ, can be set to zero.
Based on the test results of load components and the developed real and reactive power models as presented in these tables, the following comments on the reactive power models are important.
The component models appearing in Table 25.1 and Table 25.2 can be combined and synthesized to create other more convenient models. These convenient models fall into two basic forms: exponential and polynomial.
The exponential form for both real and reactive power is expressed in Equations 25.1 and 25.2 below as a function of voltage and frequency, relative to initial conditions or base values. Note that neither temperature nor torque appear in these forms. Assumptions must be made about temperature and/or torque values when synthesizing from component models to these exponential model forms.
(25.1)
(25.2)
The per-unit models of Equations 25.1 and 25.2 are as follows.
(25.3)
(25.4)
The ratio Qo/Po can be expressed as a function of power factor (pf) where ± indicates a lagging/leading power factor, respectively.
After substituting R for Qo/Po, Equation 25.4 becomes the following.
(25.5)
Equations 25.1 and 25.2 (or 25.3 and 25.5) are valid over the voltage and frequency ranges associated with tests conducted on the individual components from which these exponential models are derived. These ranges are typically ±10% for voltage and ±2.5% for frequency. The accuracy of these models outside the test range is uncertain. However, one important factor to note is that in the extreme case of voltage approaching zero, both P and Q approach zero.
EPRI-sponsored research resulted in model parameters such as found in Table 25.3 (EPRI, 1987; Price et al., 1988). Eleven model parameters appear in this table, of which the exponents α and β and the power factor (pf) relate directly to Equations 25.3 and 25.5. The first six parameters relate to general load models, some of which include motors, and the remaining five parameters relate to nonmotor loads—typically resistive type loads. The first is load power factor (pf). Next in order (from left to right) are the exponents for the voltage (αv, αf) and frequency (βv, βf) dependencies associated with real and reactive power, respectively. Nm is the motor-load portion of the load. For example, both a refrigerator and a freezer are 80% motor load. Next in order are the power factor (pfnm) and voltage (αvnm, αfnm) and frequency (βvnm, βfnm) parameters for the nonmotor portion of the load. Since the refrigerator and freezer are 80% motor loads (i.e., Nm = 0.8), the nonmotor portion of the load must be 20%.
Parameters for Voltage and Frequency Dependencies of Static Loads
Component/Parameters |
pf |
α v |
α f |
β v |
β f |
Nm |
pfnm |
α vnm |
α fnm |
β vnm |
β fnm |
Resistance space heater |
1.0 |
2.0 |
0.0 |
0.0 |
0.0 |
0.0 |
— |
— |
— |
— |
— |
Heat pump space heater |
0.84 |
0.2 |
0.9 |
2.5 |
−1.3 |
0.9 |
1.0 |
2.0 |
0.0 |
0.0 |
0.0 |
Heat pump/central AC |
0.81 |
0.2 |
0.9 |
2.5 |
−2.7 |
1.0 |
— |
— |
— |
— |
— |
Room air conditioner |
0.75 |
0.5 |
0.6 |
2.5 |
−2.8 |
1.0 |
— |
— |
— |
— |
— |
Water heater and range |
1.0 |
2.0 |
0.0 |
0.0 |
0.0 |
0.0 |
— |
— |
— |
— |
— |
Refrigerator and freezer |
0.84 |
0.8 |
0.5 |
2.5 |
−1.4 |
0.8 |
1.0 |
2.0 |
0.0 |
0.0 |
0.0 |
Dish washer |
0.99 |
1.8 |
0.0 |
3.5 |
−1.4 |
0.8 |
1.0 |
2.0 |
0.0 |
0.0 |
0.0 |
Clothes washer |
0.65 |
0.08 |
2.9 |
1.6 |
1.8 |
1.0 |
— |
— |
— |
— |
— |
Incandescent lighting |
1.0 |
1.54 |
0.0 |
0.0 |
0.0 |
0.0 |
— |
— |
— |
— |
— |
Clothes dryer |
0.99 |
2.0 |
0.0 |
3.3 |
−2.6 |
0.2 |
1.0 |
2.0 |
0.0 |
0.0 |
0.0 |
Colored television |
0.77 |
2.0 |
0.0 |
5.2 |
−4.6 |
0.0 |
— |
— |
— |
— |
— |
Furnace fan |
0.73 |
0.08 |
2.9 |
1.6 |
1.8 |
1.0 |
— |
— |
— |
— |
— |
Commercial heat pump |
0.84 |
0.1 |
1.0 |
2.5 |
−1.3 |
0.9 |
1.0 |
2.0 |
0.0 |
0.0 |
0.0 |
Heat pump comm. AC |
0.81 |
0.1 |
1.0 |
2.5 |
−1.3 |
1.0 |
— |
— |
— |
— |
— |
Commercial central AC |
0.75 |
0.1 |
1.0 |
2.5 |
−1.3 |
1.0 |
— |
— |
— |
— |
— |
Commercial room AC |
0.75 |
0.5 |
0.6 |
2.5 |
−2.8 |
1.0 |
— |
— |
— |
— |
— |
Fluorescent lighting |
0.90 |
0.08 |
1.0 |
3.0 |
−2.8 |
0.0 |
— |
— |
— |
— |
— |
Pump, fan, (motors) |
0.87 |
0.08 |
2.9 |
1.6 |
1.8 |
1.0 |
— |
— |
— |
— |
— |
Electrolysis |
0.90 |
1.8 |
−0.3 |
2.2 |
0.6 |
0.0 |
— |
— |
— |
— |
— |
Arc furnace |
0.72 |
2.3 |
−1.0 |
1.61 |
−1.0 |
0.0 |
— |
— |
— |
— |
— |
Small industrial motors |
0.83 |
0.1 |
2.9 |
0.6 |
−1.8 |
1.0 |
— |
— |
— |
— |
— |
Industrial motors large |
0.89 |
0.05 |
1.9 |
0.5 |
1.2 |
1.0 |
— |
— |
— |
— |
— |
Agricultural H2 O pumps |
0.85 |
1.4 |
5.6 |
1.4 |
4.2 |
1.0 |
— |
— |
— |
— |
— |
Power plant auxiliaries |
0.80 |
0.08 |
2.9 |
1.6 |
1.8 |
1.0 |
— |
— |
— |
— |
— |
A polynomial form is often used in a Transient Stability program. The voltage dependency portion of the model is typically second order. If the nonlinear nature with respect to voltage is significant, the order can be increased. The frequency portion is assumed to be first order. This model is expressed as follows.
(25.6)
(25.7)
where
ao + a1 + a2 = 1
bo + b1 + b2 = 1
Dp ≡ real power frequency damping coefficient, per unit
Dq ≡ reactive power frequency damping coefficient, per unit
Δf ≡ frequency deviation from scheduled value, per unit
The per-unit form of Equations 25.6 and 25.7 is the following.
(25.8)
(25.9)
The two previous kinds of models may be combined to form a synthesized static model that offers greater flexibility in representing various load characteristics (EPRI, 1987; Price et al., 1988). The mathematical expressions for these per-unit models are the following.
(25.10)
(25.11)
where
(25.12)
(25.13)
(25.14)
The expressions for the reactive components have similar structures. Devices used for reactive power compensation are modeled separately.
The flexibility of the component models given here is sufficient to cover most modeling needs. Whenever possible, it is prudent to compare the computer model to measured data for the load.
Table 25.4 provides typical values for the frequency damping characteristic, D, that appears in Equations 25.6 through 25.9, 25.13, and 25.14 (EPRI, 1979; Warnock and Kilpatrick, 1986). Note that nearly all of the damping coefficients for reactive power are negative. This means that as frequency declines, more reactive power is required which can cause an exacerbating effect for low-voltage conditions.
Static Load Frequency Damping Characteristics
Frequency Parameters |
||
Component |
Dp |
Dq |
Three-phase central AC |
1.09818 |
−0.663828 |
Single-phase central AC |
0.994208 |
−0.307989 |
Window AC |
0.702912 |
−1.89188 |
Duct heater w/blowers |
0.528878 |
−0.140006 |
Water heater, electric cooking |
0.0 |
0.0 |
Clothes dryer |
0.0 |
−0.311885 |
Refrigerator, ice machine |
0.664158 |
−1.10252 |
Incandescent lights |
0.0 |
0.0 |
Florescent lights |
0.887964 |
−1.16844 |
Induction motor loads |
1.6 |
−0.6 |
Both models provide good representation around rated or nominal voltage. The accuracy of the exponential form deteriorates when voltage significantly exceeds its nominal value, particularly with exponents (α) greater than 1.0. The accuracy of the polynomial form deteriorates when the voltage falls significantly below its nominal value when the coefficient ao is non zero. A nonzero ao coefficient represents some portion of the load as constant power. A scheme often used in practice is to use the polynomial form, but switch to the exponential form when the voltage falls below a predetermined value.
Some load components have time-dependent characteristics that must be considered if a sequence of studies using static models is performed that represents load changing over time. Examples of such a study include Voltage Stability and Transient Stability. The devices that affect load modeling by contributing abrupt changes in load over periods of time are listed below.
Protective relays: Protective relays are notoriously difficult to model. The entire load of a substation can be tripped off line or the load on one of its distribution feeders can be tripped off line as a result of protective relay operations. At the utilization level, motors on air conditioner units and motors in many other residential, commercial, and industrial applications contain thermal and/or over-current relays whose operational behavior is difficult to predict.
Thermostatically controlled loads: Air conditioning units, space heaters, water heaters, refrigerators, and freezers are all controlled by thermostatic devices. The effects of such devices are especially troublesome to model when a distribution load is reenergized after an extended outage (cold-load pickup). The effect of such devices to cold-load pickup characteristics can be significant.
Voltage regulation devices: Voltage regulators, voltage controlled capacitor banks, and automatic LTCs on transformers exhibit time-dependent effects. These devices are present at both the bulk power and distribution system levels.
Discharge lamps (mercury vapor, sodium vapor, and fluorescent lamps): These devices exhibit time-dependent characteristics upon restart, after being extinguished by a low-voltage condition—usually about 70%–80% of rated voltage.
The static load models found in Table 25.1 and Table 25.2 can be used to define a composite load referred to as the “load window” mentioned earlier. In this scheme, a distribution substation load or one of its feeder loads is defined in as much detail as desired for the model. Using the load window scheme, any number of load windows can be defined representing various composite loads, each having as many load components as deemed necessary for accurate representation of the load. Figure 25.2 illustrates the load window concept. The width of each subwindow denotes the percentage of each load component to the total composite load.
Construction of a load window requires certain load data be available. For example, load saturation and load diversity data are needed for various classes of customers. These data allow one to (1) identify the appropriate load components to be included in a particular load window, (2) assign their relative percentage of the total load, and (3) specify the diversified total amount of load for that window. If load modeling is being used for Transient Stability or Operator Training Simulator programs, frequency dependency can be added. Let P(V) and Q(V) represent the composite load models for P and Q, respectively, with only voltage dependency (as developed using components taken from Table 25.1 and Table 25.2). Frequency dependency is easily included as illustrated below.
Table 25.5 shows six different composite loads for a summer season in the southwestern portion of the United States. This “window” serves as an example to illustrate the modeling process. Note that each column must add to 100%. The entries across from each component load for a given window type represent the percentage of that load making up the composite load.
Composition of Six Different Load Window Types
LW 1 |
LW 2 |
LW 3 |
LW 4 |
LW 5 |
LW 6 |
|
Load Window Type |
Res. 1 |
Res. 2 |
Res. 3 |
Com 1 |
Com 2 |
Indust |
Load Component |
(%) |
(%) |
(%) |
(%) |
(%) |
(%) |
3-phase central AC |
25 |
30 |
10 |
35 |
40 |
20 |
Window type AC |
5 |
0 |
20 |
0 |
0 |
0 |
Duct heater with blower |
5 |
0 |
0 |
0 |
0 |
0 |
Water heater, range top |
10 |
10 |
10 |
0 |
0 |
0 |
Clothes dryer |
10 |
10 |
10 |
0 |
0 |
0 |
Refrigerator, ice machine |
15 |
15 |
10 |
30 |
0 |
0 |
Incandescent lights |
10 |
5 |
10 |
0 |
0 |
0 |
Fluorescent lights |
20 |
30 |
30 |
25 |
30 |
10 |
Industrial (induct. motor) |
0 |
0 |
0 |
10 |
30 |
70 |
EPRI User’s Manual—Extended Transient/Midterm Stability Program Package, version 3.0, June 1992.
General Electric Company, Load modeling for power flow and transient stability computer studies, EPRI Final Report EL-5003, January 1987 (four volumes describing LOADSYN computer program).
Kundur, P., Power System Stability and Control, EPRI Power System Engineering Series, McGraw-Hill, Inc., New York, pp. 271–314, 1994.
Price, W.W., Wirgau, K.A., Murdoch, A., Mitsche, J.V., Vaahedi, E., and El-Kady, M.A., Load modeling for power flow and transient stability computer studies, IEEE Trans. Power Syst., 3(1): 180–187, February 1988.
Taylor, C.W., Power System Voltage Stability, EPRI Power System Engineering Series, McGraw-Hill, Inc., New York, pp. 67–107, 1994.
University of Texas at Arlington, Determining load characteristics for transient performances, EPRI Final Report EL-848, May 1979 (three volumes).
University of Texas at Arlington, Effect of reduced voltage on the operation and efficiency of electrical loads, EPRI Final Report EL-2036, September 1981 (two volumes).
University of Texas at Arlington, Effect of reduced voltage on the operation and efficiency of electrical loads, EPRI Final Report EL-3591, June 1984 and July 1985 (three volumes).
Warnock, V.J. and Kirkpatrick, T.L., Impact of voltage reduction on energy and demand: Phase II, IEEE Trans. Power Syst., 3(2): 92–97, May 1986.