HUTECH
B GIÁO DC VÀ ÀO TO
TRNG I HC K THUT CÔNG NGH TP. HCM
PHÙNG S TIN NG DNG CÔNG C WAVELET D
BÁO PH TI IN TI HUYN ÔNG
ANH-HÀ NI
LUN VN THC S
Chuyên ngành : Thit b, mng & nhà máy đin
Mã s ngành: 60 52 50
HUTECH
TP. H CHÍ MINH, tháng 4 nm 2012
HUTECH
B GIÁO DC VÀ ÀO TO
TRNG I HC K THUT CÔNG NGH TP. HCM
HUTECH
HUTECH
CÔNG TRÌNH C HOÀN THÀNH TI
TRNG I HC K THUT CÔNG NGH TP. HCM
Cán b hng dn khoa hc :
PGS.TS. PHAN TH THANH BÌNH
(Ghi rõ h, tên, hc hàm, hc v và ch ký) Lun vn Thc s đc bo v ti Trng i hc K thut Công ngh
TP. HCM ngày 22 tháng 04 nm 2012
Thành phn Hi đng đánh giá Lun vn Thc s gm:
(Ghi rõ h, tên, hc hàm, hc v ca Hi đng chm bo v Lun vn Thc s)
1. ……………………………………………………………
2. ……………………………………………………………
3. ……………………………………………………………
4. ……………………………………………………………
5. ……………………………………………………………
NHIM V LUN VN THC S
H tên hc viên: Phùng S Tin Gii tính: Nam
Ngày, tháng, nm sinh: 05/06/1976 Ni sinh: Thái
Nguyên
Chuyên ngành: Thit b, mng & nhà máy đin MSHV:
1081031052
I- TÊN TÀI:
“NG DNG CÔNG C
WAVELET D BÁO PH TI IN TI
HUYN ÔNG ANH-HÀ NI”
II- NHIM V VÀ NI DUNG:
- C s khoa hc lý thuyt và thc nghim ca công tác d báo ph ti
đin.
- Nghiên cu và áp dng phng pháp d báo Box-Jenkins.
- Kt hp phân tích Wavelet vi phng pháp d báo Box-jenkins nhm đt đc mô h
ình d
- ng dng d báo ph ti đin ti huyn ông anh – Hà ni.
III- NGÀY GIAO NHIM V: Ngày 15 tháng 09 nm 2011
IV- NGÀY HOÀN THÀNH NHIM V: Ngày 15 tháng 03 nm 2012
V- CÁN B HNG DN:
PGS.TS.PHAN TH THANH BÌNH
CÁN B HNG DN KHOA QUN LÝ CHUYÊN
NGÀNH
(H tên và ch ký) (H tên và ch ký)
L
I CM N
Li đu tiên xin cm n c quan ni tôi công tác (Trng Cao ng Ngh
K Thut Công Ngh-B Lao đng Thng binh và Xã hi) đã to mi điu
kin v thi gian cng nh b trí sp xp công vic cho tôi phù hp đ tôi tham
gia đc khóa hc này. Xin chân thành cm n các thy cô giáo đã ging dy
cho lp Cao hc Thit b, mng & nhà máy đin khóa I vi tt c lòng bit n
và kính trng.
Xin chân thành cm n Trng i hc K thut Công ngh Thành ph
H Chí Minh, Ban giám hiu và các thy cô trong phòng sau đi hc và hp tác
quc t, khoa C-in-in t, đc bit là cô PGS.TS Phan Th Thanh Bình
Trng khoa C-in-in t, ging viên trng i hc K thut Công ngh
Thành ph H Chí Minh đã tn tình giúp đ, hng dn, cung cp tài liu cho
tôi trong quá trình thc hin lun vn này.
Cm n tt c nhng ngi bn trong lp đã k vai sát cánh cùng tôi chia
s và hoàn thin trong khóa hc này.
Cui cùng vì trình đ kin thc có hn, thi gian đu t cho vic ly s
liu lun vn và hoàn thin lun vn còn hn ch. Do đó lun vn s có sai sót,
kính mong các thy cô, các bn trong lp đóng góp đ lun vn hoàn thin.
Xin trân thành cm n Hc viên thc hin Lun vn
HUTECH
thc tin cao, đáp ng đc nhu cu phát trin tng lai ca ngành đin.
2. GII THUT CHNG TRÌNH VÀ KT QU NG DNG:
2.1 THUT TOÁN C BN CA CHNG TRÌNH
+ Phng pháp d báo Box-Jenkins da trên c s hành vi quá kh
ca chính
bin đang cn d báo.
+ Phng pháp d báo Box-Jenkins thông qua mô hình d báo
ARMA (quá
trình trung bình trt AR và t quy hi MA).
+ 2 mô hình
ARMA m rng:
- Mô hình ARIMA: Áp dng cho chui có
xu th không n đnh tng
hay gim theo thi gian, ta đa chui v dng dng loi b yu t xu th.
- Mô hình SARIMA: Áp dng cho chui có yu t
mùa vi chu k bng
4 (quý), 12 (tháng), ta đa chui v dng dng loi b yu t mùa.
2.2 CÁC BC CA PHNG PHÁP D BÁO BOX-JENKINS.
Bc 1:
nh dng mô hình (Tìm các giá tr p, d, q): ây là bc quan trng
nht đ nhn bit tt c các h mô hình ARMA và đa ra mô hình nào thích
hp và ti u nht.
-a chui v dng dng ( Xác đnh bc d ca xu th, xác đnh chu k, ly sai
phân đ loi b yu t xu th và chu k)
-Xác đnh p, q ca mô hình nh vào biu đ tng quan (Phân tích đ th trong
các trng hp đn gin)
Bc 2:
c lng các h s ca mô hình:
-Mô hình AR(p) áp dng phng pháp bình phng cc tiu hoc dung quan
-Tng hp chui d báo xp x và chi tit ra kt qu cn d báo.
2.4 KT QU D BÁO BNG PHNG PHÁP BOX-JENKINS:
Nhp chui d liu 10 ngày đ d báo t gi 217 đn gi 240 (24 gi). th
chui ph ti cn d báo (chui thc t) :
HUTECH
Ph ti có tính mùa, ta loi b tính cht mùa vì nm trong d báo có tính
trung hn và dài hn. Ph ti có tính lp li sau 24h, ta đa chui v chui
dng. Phân tích hàm tng quan đn ACF và hàm t tng quan riêng PACF,
ta thy hàm tng quan đn ACF gim đu còn hàm t tng quan riêng PACF
gim đt ngt sau đ tr p=1. Ta chn mô hình AR(1,0).
Ví d ng vi h s AR = 0.8473; C=0.0021 ta dn ti mô hình t hi
quy là
111
8473.00021.0
+=+=
ttt
yyy
jm
Kt qu sai s trung bình Mape :
tb = 6.3468%
th chui d báo và chui thc t
0 50 100 150 200 250
100
150
200
250
300
là chui chi tit, chui A
3
là chui xp x, các chui
đc xây dng nh sau:
A
3
=
å
= Mmtp
mnmn
, ,1);(
j
f
; D
n
=
å
=
n
nmnm
Mmtp , ,1);(
y
y
Ta tin hành phân tích tng chui nh đã xây dng mô hình ARIMA
cho chui d liu đu vào. Xây dng mô hình xp x và d báo cho riêng tng
chui.
Sau khi đã phân tích tng chui, ta tng hp li thành chui d báo p
n
’:
HUTECH
ABSTRACT
1. INTRODUCE ABOUT SUBJECT:
The electricity load forecasting plays an important role in ensuringsafe working
regime to be effective and economical power system. On the other hand it is
the basis to conduct strategic planning to develop future power system.
Thus the accuracy of the forecasts need special attention
.Dong Anh District of
Hanoi city where there is rapid growth of economic and social,
Parallel with that development, the demand for power users more
pressing particularly with the emergence of many industrial parks, factories and
enterprises. During the recent shortages, power fluctuations become
increasingly popular
. The variation depends on the time of the dayled
to the use of power customers were affected
.Thereforeelectricity load forecasts
in this increasingly more necessary than ever to plan
, operate an efficient grid.
How to solve manipulation by classical experience nature on computation
power load forecast is inaccurate, how it just purely based on the experiences
of the past that the experience was notalways true when applied
to continuously changing circumstancesin the future. To overcome the above
disadvantages, the writer offers analysis method of wavelet transform.
Forecasting model application of wavelet transform analysis for short-
with the available data, if not choose theappropriate model to go back to step
one.
Coefficient of the model must be different from 0, if not satisfied,
then removed from the ARMA model.
-Inspection of the rest of the remainder, if the remainder is a white noise is
satisfactory if not satisfied then repeat the process until you
find the appropriate model.
Step 4: Forecast: Forecast limited to a few cycles (short term)because of
increasing amplitude error of prediction time.
-When we calculate the predicted value predicted for anupcoming time period.
-Using the predicted values to calculate the value of the nextforecast.
-Continue until such time period specified.
2.3 STEPS OF COMBINATION WAVELET ANALYSIS
PREDICTION METHOD BOX-JENKINS.
Step 1: Analyze and reconstruct the signal by Wavelet transform:
-Split signal into approximation and detail sequences (Waveletthey choose to
comply, the new forecast accuracy).
Step 2: Perform a prediction for each sequence: Perform the stepsby the
method of Box-Jenkins forecasting
Step 3: Re-structure prediction results using wavelet transformback to
the forecast results.
Chain-sum approximation and forecast information to forecastresults.
2.4 RESULTS FORECAST BY BOX-JENKINS:
Enter the data for 10days forecast from 217 hours to 240 hours (24
hours). Graph to chain load forecasts (actual sequence):
HUTECH
The load is the season, we remove the seasonal nature in theforecast for
the medium and long term. Repeatable load after 24 hours, we put chains on
the chain stops. Analysis of single-
2.5 RESULTS FORECAST COMBINATION WITH WAVELET
TRANSFORM BOX-JENKINS METHODOLOGY:
Using wavelet analysis of type 3 with a db5-wave pattern. Pnsequence
HUTECH
data will then be broken down into four series: p
n
= D
1
+D
2
+D
3
+A
3
Sequence D1, D2, D3 is more string, the string is a sequence of
approximately A3, the chain is built as follows:
A
3
=
å
= Mmtp
mnmn
, ,1);(
j
f
; D
n
=
å
350
400
450
500
550
chuoi du bao v-duong dut net; chuoi thuc te u-duong lien netHUTECH
NG DNG CÔNG C WAVELET D BÁO PH
TI IN TI HUYN ÔNG ANH-HÀ NI
APPLICATION WAVELET TOOL FOR FORECASTING LOAD IN
THE DISTRICT OF DONG ANH-HA NOI
H tên hc viên: PHÙNG S TIN
n v công tác: Trng khoa CNTT, Trng C Ngh K Thut Công Ngh-B LTB&XH
T 59 th trn ông Anh Hà Ni, T: 04.39686791, Fax: 04.38832436
Di đng: 098.3393834, E-mail:
Cán b hng dn: PGS.TS PHAN TH THANH BÌNH
n v công tác: Trng khoa C-in-in t, Trng i Hc K Thut Công Ngh TP HCM
TÓM TT
Vic d báo ph ti đin đóng vai trò rt quan trng trong vic đm bo ch đ làm vic an toàn đt
hiu qu và tit kim ca h thng đin. Mt khác nó là c s đ tin hành hoch đnh chin lc phát
trin h thng đin trong tng lai. Do đó tính chính xác ca d báo cn phi đc bit chú trng.
Huyn ông anh thành ph Hà ni là ni có tc đ phát trin nhanh v kinh t xã hi, song song vi
s phát trin đó thì nhu cu s dng đin nng ngày càng cp thit đc bit vi s xut hin ca nhiu
khu công nghip, nhà máy, xí nghip. Trong thi gian va qua thì s thiu ht, bin đng đin nng
ngày càng ph bin. S bin đng ph thuc vào tng thi đim trong ngày dn ti quá trình s dng
đin nng ca khách hàng b nh hng nng n. Do đó d báo ph ti đin ti đây càng tr lên cn
thit hn bao gi ht đ hoch đnh, vn hành li đin mt cách hiu qu.
hình d báo áp dng phân tích bin đi
Wavelet dùng cho d báo ngn hn cho ta kt
qu d báo khá chính xác, có tính ng dng
thc tin cao, đáp ng đc nhu cu phát
trin tng lai ca ngành đin.
2.
NI DUNG
2.1 THUT TOÁN C BN CA
CHNG TRÌNH
+ Phng pháp d báo Box -Jenkins da trên
c s
hành vi quá kh
+ Phng pháp d báo Box -Jenkins thông qua
mô hình d báo
ARMA (quá trình trung bình
trt AR và t quy hi MA).
ca chính bin đang
cn d báo.
+ 2 mô hình
ARMA m rng:
- Mô hình ARIMA: Áp dng cho chui có
xu
th
không n đnh tng hay gim theo thi
gian, ta đa chui v dng dng loi b yu t
xu th.
- Mô hình SARIMA: Áp dng cho chui có
yu t
-Kim đnh tính dng ca các phn d, nu
phn d là mt nhiu trng thì đt yêu cu nu
không tha thì lp li quá trình cho đn khi tìm
đc mô hình thích hp.
Bc 4:
-Khi d báo ta tính giá tr d báo cho mt thi
đon sp ti.
D báo: D báo gii hn trong mt
vài chu k (ngn hn) vì biên đ sai s tng
nhanh theo thi gian d báo.
-S dng giá tr d báo đ tính cho giá tr d
báo thi đon tip theo.
HUTECH
-C tip tc nh vy cho đn khi ht thi đon
xác đnh.
2.3 CÁC BC CA PHÂN TÍCH
WAVELET KT HP PHNG PHÁP
D BÁO BOX-JENKINS.
Bc 1:
-Tách tín hiu thành các chui xp x và chi
tit (chn h Wavelet phù hp thì d báo mi
chính xác).
Phân tích và tái cu trúc tín hiu
qua Wavelet transform:
Bc 2:
Tin hành d báo cho tng chui:
Thc hin các bc theo phng pháp d báo
Box-Jenkins
Bc 3:
-Tng hp chui d báo xp x và chi tit ra
−−
+=+=
ttt
yyy
ϕµ
Kt qu sai s trung bình Mape :
tb = 6.3468%
th chui d báo và chui thc t
0
50 100 150 200
250
100
150
200
250
300
350
400
450
500
550
chuoi du bao v-duong dut net; chuoi thuc te u-duong lien net
3.2 KT QU D BÁO BNG BIN I
WAVELET KT HP PHNG PHÁP
BOX-JENKINS :
S dng phân tích Wavelet bc 3 vi sóng mu
là db5. Chui d liu p
ϕ
φ
= ; D
n
∑
=
n
nmnm
Mmtp , ,1);(
ψ
ψ
=
HUTECH
Ta tin hành phân tích tng chui nh đã xây
dng mô hình ARIMA cho chui d liu đu
vào. Xây d ng mô hình xp x và d báo cho
riêng tng chui.
Sau khi đã phân tích tng chui, ta tng hp
li thành chui d báo p
n
’:
p
n
’=D
1
’ +D
2
’ +D
3
(24h), vi yêu cu phi thng xuyên cp
nht giá tr thc t nu không sai s s ln.
có đc mt kt qu d báo phù hp và
kh thi, ngoài dùng yu t kinh nghim ra còn
phi chn h Wavelet phân tích phù hp vi
chui d liu vì đây là yu t quyt đnh đ
chính xác ca d báo.
Kiu d liu d báo trên đc biu din cho
đin nng, tuy nhiên nu s liu khác thì cng
có th biu din tng t nh trên nh d báo
giá vàng, d báo thi tit…
y) ta s dng
bin đi Wavelet kt hp phng pháp d báo
Box-Jenkins, còn d báo dài hn nên s dng
phng pháp d báo tng quan s cho kt
qu chính xác hn.
4.2 HNG PHÁT TRIN
đin không sát vi
thc t, mt s ngày t
ln. n
ta cn khc phc nhc đim này.
ng dng mô hình d báo ph ti ngày vào d
báo giá đin trong tng lai khi đin nng tr
thành hàng hoá, vic giao dch mua bán theo
qui lut th tr ng, d báo chng khoán, d
báo giá vàng, d báo thi tit ….
Thit k mô hình d báo ph ti đin cho h
thng vùng min.
Vic d báo ph ti đin đóng vai trò rt quan trng trong vic đm bo ch đ làm vic an
toàn đt hiu qu và tit kim ca h thng đin. Mt khác nó là c s đ tin hành hoch
đnh chin lc phát trin h thng đin trong tng lai. Do đó tính chính xác ca d báo cn
phi đc bit chú trng. Huyn ông anh thành ph Hà ni là ni có tc đ phát trin nhanh
v kinh t xã hi, song song vi s phát trin đó thì nhu cu s dng đin nng ngày càng cp
thit đc bit vi s xut hin ca nhiu khu công nghip, nhà máy, xí nghip. Trong thi gian
va qua thì s thiu ht, bin đng đin nng ngày càng ph bin. S bin đng ph thuc
vào tng thi đim trong ngày dn ti quá trình s dng đin nng ca khách hàng b nh
hng nng n. Do đó d báo ph ti đin ti đây càng tr lên cn thit hn bao gi ht đ
hoch đnh, vn hành li đin mt cách hiu qu.
ABSTRACT
The electricity load forecasting plays an important role in ensuringsafe working regime to be
effective and economical power system. On the other hand it is the basis to conduct strategic
planning to develop future power system. Thus the accuracy of the forecasts need special
attention. Dong Anh District of Hanoi city where there is rapid growth of economic and
social, Parallel with that development, the demand for power users more
pressing particularly with the emergence of many industrial parks, factories and enterprises.
During the recent shortages
, power fluctuations become increasingly popular. The
variation depends on the time of the dayled
to the use of power customers were affected
.Thereforeelectricity load forecasts in this
increasingly more necessary than ever to plan
, operate an efficient grid. HUTECH
1. GII THIU
trung bình trt AR và t quy hi MA).
+ 2 mô hình
ARMA m rng:
- Mô hình ARIMA: Áp dng cho chui có
xu th không n đnh tng hay gim theo
thi gian, ta đa chui v dng dng loi
b yu t xu th.
- Mô hình SARIMA: Áp dng cho chui
có yu t
mùa vi chu k bng 4 (quý), 12
(tháng), ta đa chui v dng dng loi b
yu t mùa.
2.2 CÁC BC CA PHNG PHÁP
D BÁO BOX-JENKINS.
Bc 1:
nh dng mô hình (Tìm các giá
tr p, d, q): ây là bc quan trng nht đ
nhn bit tt c các h mô hình ARMA và
đa ra mô hình nào thích hp và ti u
nht.
-a chui v dng dng ( Xác đnh bc d
ca xu th, xác đnh chu k, ly sai phân
đ loi b yu t xu th và chu k)
-Xác đnh p, q ca mô hình nh vào biu
đ tng quan (Phân tích đ th trong các
trng hp đn gin)
Bc 2:
c lng các h s ca mô
hình:
-Mô hình AR(p) áp dng phng pháp