Load Forecasting
Load Forecasting
Load Forecasting
LOAD FORECASTING
ØCLASSIFICATION OF LOAD.
ØCHARACTERISTIC OF LOAD.
ØAPPROACHES TO LOAD FORECASTING.
ØFORECASTING METHODOLOGY.
1.EXTRAPOLATION .
2.CORRELATION.
REQUIREMENTS FOR FORECAST
ØCATEGORIZATION OF CONSUMPTION SECTORS.
ØCOLLECTION & COMPILATION OF PAST DATA.
ØDECOMPOSITION OF DATA.
ØVALIDATION OF PAST DATA & ANALYSIS.
ØSELECTION OF SUITABLE FORECAST
TECHNIQUE.
ØCOLLECTION OF DATA OF ECONOMIC VARIABLES
FOR ECONOMETRIC MODELLING.
KEY ISSUES
ØENERGY SHORTAGE DUE TO RESTRICTIONS &
UNSCHEDULED CUTS.
ØENERGY EFFICIENCY & CONSERVATION
MEASURES.
ØSPECIAL IMPETUS ON RURAL ELECTRIFICATION &
IRRIGATION.
ØTHRUST TO ECONOMICALLY WEAKER STATES &
REGIONS
Ø
FACTORS
INFLUENCING ELECTRICITY DEMAND
ØDEMOGRAPHIC GROWTH.
ØGEOGRAPHIC VARIATIONS.
ØMETEOROLOGICAL VARIATIONS.
ØNATIONAL POLICIES FOR SECTORAL GROWTH
§NATIONAL ELECTRICITY POLICY
§INFRASTRUCTURE DEVELOPMENT PLANS.
ØECONOMIC GROWTH.
ØDIVERSITY FACTORS-INTER STATE, INTER
REGIONAL.
ØINTERNATIONAL ECONOMY LIKE FUEL /OIL PRICE,
ETC.
ELECTRICITY DEMAND FORECASTING:
A NECESSITY
Deterministic
Probabilistic
Stochastic
EXTRAPOLATION
Extrapolation techniques involve fitting trend
curves to basic historical data adjusted to
reflect the growth trend .
Forecast obtain by trend curve
Deterministic extrapolation
Probabilistic Extrapolation
Uncertainty of extrapolated results is quantified using
statistical mean & variance.
Uncertainty arises in historical data & analytical
model.
Stochatic Extrapolation Model
Residential Sales
Forecasts
Population Method
- Major factors
customers, Builders.
- Rapid growth is an extremely important step
Commercial Sales
Forecast
Commercial establishments are usually
service oriented
Hence growth patterns are close to
residential sales
By using method of extrapolation ratios of
commercial to residential sales in to
the future & then multiply by residential
sales forecast
Another approach to extrapolate
historical commercial sales.
Industrial sales forecast
Industrial sales are tied very closely to the overall economy
& overall economy not stable over selected periods.
Approaches in view
1. Multiply forecasted production levels by forecasted
energy consumption per unit of production
2. Multiply forecasted number of industrial workers by
forecasted energy consumption per worker.
ε(η(t) η(t+Ɣ))=0 Ɣ≠0 property 3
= ση Ɣ=0
σ 2 = Variance large σ 2indicates noise vary about its
η η
mean
σ 2 =0 means η=ň
η
Σ defines 99% confidence interval
η
Random weekly peak historical data gives values of
coefficients on the fitting functions.
That is done by minimizes squared difference
between weekly peak demand l (t) & true peak
demand ξ(t)
n
function J= ∑1/2W2(t)[l(t)-ξ(t)]2
t=1
J=1/2[L(N)-a’R(N)]W(N)W’(N)[L(N)-a’R(N)]’
Where, L(N)=1*N Vector whose elements are l(t)
a’= 1*n Fitting function coefficient vector
R(N)=n*N fitting function matrix whose columns
are f(t)
By taking partial derivative
â= [L(N)W(N)W’(N)R’(N)][R(N)W(N)W’(N)R’(N)]-1
Letting, F(N)=[L(N)W(N)W’(N)R’(N)]
G(N)=)[R(N)W(N)W’(N)R’(N)]
â= F(N) G(N)-1
Recursive formula F(N)
Fik =∑N fi(t)ᵦN-t fk(t)
t=0
= fi(N)fk(N)+ ∑N-1 fi(t)ᵦN-1-t fk(t)
But ∑N-1 fi(t)ᵦN-1-t fk(t) = Fik (N-1)
t=1
Thus recursive relationship for the ikth element of F
Fik = fi(1) fk(1)
Initial values of Fik for fitting function describing 3 bus
system
F11 = f12(1) =0.1
F12 = f1(1) f2(1)= 0.1
F21 = f2(1) f1(1)=0.1
F22 = f22(1) =0.1
Recursive formula G(N)
Gn =∑N l(t)ᵦN-t fn(t)
t=1
= l n(N)fn(N)+ ᵦ ∑N-1 l (t)ᵦN-1-t fn(t)
t=1
= l n(N)fn(N)+ ᵦ Gn(N-1)
Mean & Variance
Coefficients
The Variance of a (N) qualifies the range over
1
which the elements of a(N) varies as data are
simultaneously processed.
Determination of coefficient variance
con(â(N))= ε([(â(N)-ẫ(N)] [â(N)-ẫ(N)]’)
ẫ(N) mean value of (â(N)
â(N) = [εL(N)] [L(N)W(N)W’(N)R’(N)F(N)-1 ]
since [εŋ(N)] =0
ẫ(N) = F(N)-1 R(N)W(N)W’(N)ξ’(N)
con(â(N))= ε(â(N)â’(N)-ẫ(N)â’(N)+ẫ(N)ẫ’(N))
= ση2 F(N)-1 R(N)W(N)W’(N)W(N)W’(N)R(N)F(N)-1
FORECAST VARIANCE
variance qualifies uncertainty of forecast
ε(ξ(t))= â’(N)f(t)
var(ξ(t))= ε[ξ(t)ξ’(t)][ξ(t) ξ’(t)]
= f’(t)con[â’(N)f(t)]
= σξ2