Load Static Models for Conservation Voltage Reduction in the Presence of Harmonics

The Conservation Voltage Reduction (CVR) is a technique that aims to achieve the decrease of power consumption as a result of voltage reduction. The customer is supplied with the lowest possible voltage level compatible with the stipulated level by the regulatory agency. International Standards ANSI C84.1-2006 and IEEE std 1250-1995 specify the range of supply voltage to electronics equipment from 0.9 to 1.05 pu of nominal voltage. To analyse the CVR effect in distribution systems with different load characteristics (residential, commercial, industrial or a combination of these), mathematical load models are used. Typically, these equipment/load models are used to analyse load aggregation without any consideration of its nonlinearity characteristics. Aiming to analyse the nonlinear characteristics and its consequences, this paper presents a discussion of the neglected variables as well as the results of a set of measurements of nonlinear loads. Different mathematical models are applied to obtain them for each load. Using these models the load aggregation is evaluated. It is presented that although the models show adequate results for individual loads, the same does not occur for aggregated models if the harmonic contribution is not considered. Consequently, to apply the load model in CVR it is necessary to consider the harmonics presence and the model has to be done using only the fundamental frequency data. The discussion about the causes is done and the models are compared with the measurements.


Introduction
Conservation Voltage Reduction (CVR) is a method of energy reduction consumption resulting from a reduction of feeder voltage [1].According to Schneider [1] while in North America there are numerous CVR systems deployed in distribution feeders, the majority of the results published are based on empirical field measurements.Consequently, it is not possible to extrapolate the results of this technology on other distribution feeders.To obtain a better evaluation, it is necessary to use load models and analyse the combination of these models in simulation processes.
The load model for CVR analysis is a mathematical representation of the relationship between the active and reactive power consumed and load feed voltage; in other words, it describes the load behavior when it is connected on power grid considering the voltage variation.Such models can be done without thermal cycles or static [1] [2] where the power at each instant, t, is expressed as function of voltage and frequency at the same instant ( ) ; or dynamic, where the power is expressed as function of voltage and frequency both at past and present time [2].
Loads without thermal cycles consume energy in a time-invariant manner, with the exception of voltage variations.Specifically, there is no control feedback loop.As an example, a light bulb will consume energy when turned on, as a function of voltage in a fixed manner.In contrast, a load with a thermal cycle, such as a water heater, will have a varying duty cycle dependent on the supply voltage, Schneider [1].
The load model, according to Price et al. [2], may bring benefits in planning studies, in electrical power systems operation, resulting in, for example, reduction of system modifications and equipment investments, and system preventing emergencies due to better operating limits.
Some researchers found in literature demonstrate the importance of load model application on Conservation Voltage Reduction (CVR) [1]- [11].In 1992, a research developed by Electric Power Research Institute (EPRI) [3] presented tests in order to verify the reduction voltage effects in the efficiency of loads existing in residential, commercial and industrial sectors.Such tests were divided in two parts.In the first part, data were collected to model energy consumption as a function of supply voltage.In the second part, the collected data were used to model the loads.In 1989, a field test was performed in Snohomish [4] in four circuits connected to three substations whose loads were predominantly residential and commercial, where there were alternating cycles with nominal voltage and low voltage every 24 hours.In each substation, power meters were installed to collect active power, reactive power and voltage in each phase of the test circuit, where the data were collected every 15 minutes.The test resulted in a reduction of energy consumption.However the result was not considered optimal because each substation at each instant of time had different loads connected.It is evident that the load model is needed for predicting energy saved as a result of applying the CVR method.In 2010, the study described in Schneider [1] presented the CVR application in some feeders, where the static loads were modeled by a polynomial model (ZIP), and the dynamic loads were modeled by an equivalent thermal parameter (ETP) model.The dynamic loads used a more complex model due to the control system in each load that determines the moment when it is energized and for how long.
None of those publications consider the nonlinear characteristics of the load.The contribution of the harmonics in the load behavior has to be considered in the analysis.This paper evaluates three model static load techniques as used in literature, without any consideration of nonlinearities and including their analyses.This paper is divided as following: Firstly three CVR model techniques were presented.These techniques can be applied considering or not the harmonics components presence.Then, in the second part the approach differences were shown.To obtain the model using the measured data, it was presented the parameter estimation technique.It is desired to maintain the physics aspects of the model; then the ellipsoid algorithm (using restrictions) and its application were presented.Finally the used methodology was shown as well as the results.The comparison has been done and the conclusions were shown.

End-User Load Modeling-Static Loads
End-user load models are mathematical functions used to describe the behavior of commercial or residential loads considering the active and reactive power as function of both variables: voltage and sometimes frequency.The static models used to represent the load are performed in two main types: by an exponential model and by polynomial models, also known as ZIP models [1] [2] [5] [7].

Exponential Model-EXP
In exponential models the relationship between the consumed power and the feed voltage is given by an exponential function, as in Equation (1) for the active power and in Equation (2) for the reactive power [2] [5] [10]- [12]: where α and β parameters describe the load behavior.
When both coefficients α and β are equal to zero, it means that the load behaves as Constant Power.In case of both α and β are equal to one, the load behaves as Constant Current; and finally if both coefficients are equal to two, the load behaves as Constant Impedance.

ZIP Model (ZIP)
In polynomial models, ZIP model, the load is represented by following some physical characteristics.Some equipment, such as electronic ballasts, have constant power consumption under voltage variation.On the other hand, others behave as a combination of two or three of the following characteristics: Constant Impedance (Z), Constant Current (I) and Constant Power (P).
The expressions for active and reactive power (P i , Q i ) of the ZIP model are presented by Equations ( 3) and ( 4) respectively.For both Equations, the coefficients have to obey the constraints presented in Equation (5).The relationship between power consumption and voltage magnitude is given by a polynomial Equation [2] p p p q q q Z I P Z I P where: i P : Active power consumption of the i th load., p q Z : Fraction of load that is Constant Impedance., p q I : Fraction of load that is Constant Current., p q P : Fraction of load that is Constant Power.

Modified ZIP Model-(ZIP_mod)
Schneider et al. [1] proposed a new ZIP model where, besides the information about the Constant Impedance, Constant Current and Constant Power, there is also information about the phase angle of each one of these components.In this case, both active and reactive Equations are coupled.The active and reactive power calculus is given in Equations ( 6) and ( 7) respectively, obeying the constraint in Equation ( 8): where: i P : Real power consumption of the i th load.
i Q : Reactive power consumption of the i th load.% Z : Fraction of load that is Constant Impedance.
% I : Fraction of load that is Constant Current.
% P : Fraction of load that is Constant Power.Z θ : Phase angle of the Constant Impedance component.I θ : Phase angle of the Constant Current component.P θ : Phase angle of the Constant Power component.
In any ZIP model, when the load is represented as Constant Impedance, it means that the power varies directly with the square of the voltage; in the case of the Constant Current, the power varies directly with the voltage; and, lastly, if the load is modeled as Constant Power, the power does not vary when the voltage varies [2].

Analysing the Harmonics Presence and Its Effect in the Model
None in all presented models is mentioned about harmonic components in any reference although its presence affects the data.Once measurements are used to obtain the parameters for each model, and as the components suffer influences from harmonics, it is necessary to discuss how to consider its influence before start using the data.

For Sinusoidal Voltage and Current
The apparent power measured by the equipment gives us two parcels called P and Q.For sinusoidal source supplying a linear load, these parcels are [13]:

For Nonsinusoidal Voltage and Current
Considering nonsinusoidal source supplying a nonlinear load, however, the apparent power assumes a different value including the harmonics effect.The measured active power P M now has also the contribution of the harmonic active power P H . Similarly the measured reactive power Q M is not only composed by Q 1 but also by a sum of other parcels, namely D V , D I and D H .This measured power should no longer be called reactive power, but nonactive power, [13].The demonstration is presented below:

Parameter Estimation
The parameter estimation is a procedure which uses some samples from measurements to calculate one or more unknown parameters.The used measurements (samples) in the estimation process are subject to errors, so that the estimated parameters also have associated errors [14].Equation (13) shows the value of the samples obtained by a measurement device.The measured value is close to the true value, but with the difference of an error.
where meas Z is the value obtained by measurement equipment; true Z is the true value of the measurement; and η is the random measurement error that serves to model the uncertainty in the measurements.
The statistical criteria used in the parameter estimation may be quoted the maximum likelihood, where the goal is to maximize the probability that the estimated state variable is the actual value of the state variable vector [14].Such estimator assumes that the probability density function (PDF) is known.Other criteria may also be used, for example, the weighted least-squares criterion that aims to minimize the sum of the weighted squares of the estimated measure deviations.In this case, it is not required that the PDF of samples or measurement errors be known, but it is assumed to have Gaussian normal distribution [14].
In the maximum likelihood criterion, it is necessary to estimate the variable x that maximizes the likelihood of measuring meas Z [14].According to the author, the maximum likelihood estimate of the unknown parameter is thus expressed as the parameter value that minimizes the sum of squares of the difference between each measured value and the actual measured value, with each squared difference weighted by the variance of the meter error, as in ( 14): where: ( ) i f x : function used to calculate the value being measured by the i th measurement.

( )
J x : measurement residual.N : number of independent measurements.
2 i σ : variance of the i th measurement.meas i Z : i th measurement.In one matrix approach, this Equation can be written as in (15): where: R: covariance matrix of the measurements, as can be seen in ( 16).

Ellipsoid Algorithm
This algorithm consists in finding a solution to the optimization problem through a succession of increasingly smaller ellipsoids, starting from the initial ellipsoid containing the optimal point [15].Given the constrained optimization problem described in (17) [12], the convergence is achieved if the objective function and the constraints are convex and if n x R ∈ [16].The algorithm starts with an ellipsoid o E centered at the o x point as in (18): , where: ( ) q ⋅ : depicts the constraints of inequality,

( )
h ⋅ : represents equality constraints described as two inequality constraints, and ( ) The ellipsoid algorithm can be described by recursive Equations ( 19), ( 20) and ( 21), which generate a sequence of k x points [12].
( ) where m is the gradient or sub-gradient of the function or of the most violated constraint, according to the rule described in ( 22) and (23) [12].
( ) Figure 1 illustrates the flowchart of the ellipsoid algorithm.Firstly, the variables that define i. the region which is contained in the solution (Q) (initial ellipsoid), ii.initial optimal solution (x opt ) and its objective function value (f opt ), and iii. the algorithm parameters (for example the tolerance for calculate the gradient numerically (tol)) are initialized.After that, the variable receives the value of the largest constraint function at x. Then the objective function value at the point x is calculated.If r is negative, the algorithm calculates the objective function gradient at point x (f x ) and checks whether the value of f k is smaller than the optimal value achieved so far.If f k is smaller than f opt (optimal function value), x opt and f opt receive the value of x k and f k , respectively, otherwise the algorithm goes to the next step.If r is bigger than or equal to zero, the algorithm calculates the gradient of the most violated constraint.In the next step, the algorithm computes the variables to generate the next ellipsoid based on the computed gradient vector.This process is repeated until the convergence is achieved.

Method to Obtain the Model
The data measurements were performed three times for each load and at different times.During the experiments, the voltage has been changed in steps.An energy analyzer was used to collect data.Every 5 minutes the input voltage was increased by five volts starting the measurement at 110 V and ending at 130 V.The nominal voltage is 127 V. Samples were collected every second.For each load, it was used at least 4500 measurements.These data were used in the following models.
For the exponential model, the optimization problem is given as in Equation ( 24) [6]: where: ( ) f x is the function described in (1).The same formulation is presented for the exponential reactive power model.
For ZIP model [2] the parameters are estimated according to the optimization problem presented in Equation (25) [6]: where: 0 T p p p x P I Z P   =   ; 0 P is the active power at nominal voltage and ( ) f x is the function described in (3).The same formulation is presented for the ZIP reactive power model [2].For the ZIP_mod model [1], the estimated parameters must satisfy simultaneously the Equations of active and reactive power, so that the optimization problem is defined as in Equation ( 26): where: 0,5 α = in order to minimize both active and reactive power functions at the same intensity; ( ) x is the Equation concerning the active power with ( )  7); and n S is the apparent power consumption at no- minal voltage.
The model proposes to minimize the active and reactive power functions simultaneously through a weighted sum, that is, both power functions are minimized with the same weight.Often the objective functions (active and reactive power) have different magnitude order.Therefore, one function can be minimized faster than the other.To solve this problem, the optimal objective function was calculated for the active power ( 1 J * ) and reactive ( 2 J * ) separately using the constraints of Equation ( 26), and then, the optimal values obtained were used to calculate the objective function of the problem in Equation ( 27): ( ) This paper uses the Ellipsoid optimization technique for parameter estimation obeying the model constraints described on Equations ( 24), ( 25) and (26).

Methodology
Knowing load mathematical models, it is possible to obtain a model of a group of loads by simply adding the single models, that is, the procedure used for linear loads (superposition theorem).In order to analyse the effect of the harmonics in the load models and in their aggregation, an experiment is carried out with six different loads, and the analysis is divided in two phases.
1) During the first phase, each equipment has been tested for different voltage values in the range of 0.9 and 1.05 pu.The measurements have been made repeatedly and the variables S, S 1 , P, P 1 , N, Q 1 have been registered.Using the two groups of variables (S 1 , P 1 , Q 1 and S, P, N) the mathematical models have been obtained by following the procedure previously mentioned.Consequently, two representations (i.only the fundamental fre-quency components (S 1 , P 1 , Q 1 ) and ii.harmonics contributions (S, P, N)) have been calculated for each of the three models (Exponential, ZIP and ZIP_mod) and for each equipment.
2) In the second phase, six pieces of equipment were connected simultaneously and the measurements for the same variables shown in phase 1 have been collected.Six models have been adjusted, two representations for the two groups of variables (i.only the fundamental frequency components (S 1 , P 1 , Q 1 ) and ii.harmonics contributions (S, P, N)) for each method (Exponential, ZIP and ZIP_mod).

Results
Following the methodology, the experiments have been done and the results are presented in this section.Table 1 contains the modeled loads and their rated power.Among the modeled loads, there are resistive loads, electronics and motors.

Phase 1-Models Obtained for Each Load
Figure 2(a) is presented to illustrate how the active powers P and P 1 varies when the voltage changes in steps.They are almost the same, as the active harmonic power P H is very low.Otherwise Figure 2(b) shows a big difference between Q M (that in fact is N) and Q 1 .Figure 3 presents similar results for Liquid Crystal Display measurements.In this case the difference between active powers P and P 1 is more evident.Using these data, the models, as a function of voltage, Exp, ZIP and ZIP_mod, are obtained.
Table 2 presents the parameters α and β for the exponential model expressed by Equation ( 7) and ( 8) for all modeled loads, considering the measured data (P, N) and the filtered data (P 1 , Q 1 ).As the active power suffers little influence of harmonics, the calculated parameters ( α , Po) are almost equal, even for models considering or not harmonics contribution.Otherwise, the nonactive power is strongly influenced by the presence of harmonics.Table 3 presents the parameters for active and nonactive power Equation considering the ZIP model.It can be seen that there are some subtle differences in active Equation parameters, mainly for the nonlinear loads like CFL, LCD and TV.On the other hand, significant differences are observed for nonactive power parameters.In fact the model parameters for nonactive power are completely different.
For Modified ZIP model, the results are presented in a different way.As the Equations of active and reactive power are coupled, the results are presented in Table 4, each one for different reference data S 1 , P 1 , Q 1 and S, P, N respectively.It can be seen that the results are strongly affected by the presence of harmonics.

Phase 2-Load Aggregation
In the second phase the modeled loads, used in phase I, are turned on in the same time, and measurements of total active and nonactive power are done.Using these data, the models for the group of equipments data are obtained (Collective Model).Table 5, Table 6 and Table 7 present the obtained models.
Finally, the model obtained by summing up the individual models is calculated (Sum of Models) by applying superposition theorem used in Equation (28).Total 1  a These values were obtained using the model built using data for all equipment turned on simultaneously; b These values were obtained using the model built using the sum of individual models.
values presents the reference data.It can be observed that the difference between obtained models (Exp, ZIP and ZIP_mod) of active power for both voltage levels is almost irrelevant, considering or not harmonics (P m or P 1 ).
Otherwise for nonactive power the model presents very significant differences.For fundamental frequency component (Q 1 ) only, the results present similar values.Nonetheless, when the harmonic components are considered (Q m , N), the results are completely different, since the superposition theorem cannot be applied.

Conclusion
This paper demonstrated the application of optimization technique for the treatment of ellipsoid three static load models widely used in the literature.This technique was used in order to obtain real models for the loads, so that the result achieved can produce physical meaning.All the mathematical formulations, as well as the methodology used to achieve the results are available in the initial sections.
For most of the loads, the three models showed similar results, except for electronic loads such as compact fluorescent lamp, monitor, computer and television.Such loads showed a reduction in power consumption with the increase in voltage level.Thus, it was evident that among the three physical models presented, only the ZIP_mod can faithfully represent the load, once the other models do not show curves with negative slope.
The literature presents these models and suggests using the superposition theorem to obtain the model which represents a group of loads by simple summing up the individual models.This paper shows that the results of the obtained models, which do not filter the data (use only the fundamental frequency), cannot be used for load aggregation.In other words, to apply this kind of load model it is necessary to work with only the filtered data of the load, otherwise the final model will present very significant errors.

iQ:
Reactive power consumption of the i th load.

oP:
Active power consumption at nominal voltage.

oQ:
Reactive power consumption at nominal voltage.

iV:
Actual terminal voltage.

iV:
Actual terminal voltage.

nS:
Apparent power consumption at nominal voltage.

Table 2 .
Exponential model parameters for active and nonactive power.

Table 8
presents a comparison among the obtained models outputs (Collective Model and Sum of Models), considering or not harmonics components (i.only the fundamental frequency components (S1, P1, Q1) and ii.harmonics contributions (S, P, N)), for two different voltage levels (110 and 130V).The line of the measured

Table 3 .
ZIP model parameters for active and nonactive power.

Table 4 .
Modified ZIP model parameters for active and nonactive power.

Table 5 .
Exponential model parameters for active and nonactive power.

Table 6 .
ZIP model parameters for active and nonactive power.

Table 7 .
Modified ZIP model parameters for active and nonactive power.

Table 8 .
Comparison among models for two different voltages and for different references.