ISSN: 1693-6
930
167
Sim
u
lasi Optim
a
si Aliran Daya
….(Subi
yanto).
SIMULASI OPTI
MASI ALIRAN DAYA SISTEM TENAGA
LISTRI
K SEBAGAI PENDEKATAN EFISIENSI BIAYA
OPERASI
Subiy
a
nto
Tekni
k
E
l
e
k
t
r
o,
Fakult
a
s
T
e
kni
k
,
Univ
e
r
sit
a
s
Neg
e
ri
S
e
mara
ng
Gedu
ng E6 L
t.2 Kampus S
e
ka
ran G
unu
ngpati Sema
rang
Email:
biyant@lyco
s
.com
,
biyantote_u
n
nes@yah
oo.com
Abs
t
rak
Dalam
si
ste
m
tenaga
setiap p
e
m
bangkit m
e
m
p
u
n
yai
bia
y
a
b
ahan
ba
ka
r
m
a
sing-
m
a
sing dan
berad
a pa
da jara
k be
ban yang ti
dak
sam
a
dari pu
sat. Kapa
sitas
sel
u
ru
h
pem
bang
kit
haru
s
lebih
da
ri
keb
u
tuhan
be
ban
dan
rugi-rugi. Dalam
si
stem
ten
aga
terinterko
ne
ksi perl
u
dila
kuka
n penj
ad
wala
n pen
yal
u
ran d
a
ya aktif dan rea
k
tif m
a
sing-m
a
si
ng
pem
bang
kit untuk
m
e
m
i
nim
u
m
k
an
bia
y
a ope
ra
si. Tulisa
n
ini
m
e
m
berikan
su
atu sim
u
lasi
tentang
opti
m
isasi ali
r
a
n
da
ya
si
ste
m
tenaga li
strik
yan
g
m
e
rup
a
kan
su
atu teknik u
n
tuk
m
e
m
i
nim
a
lkan biaya o
perasi si
stem
te
naga. Da
ri
si
m
u
lasi ini diharap
ka
n dap
at m
e
m
berika
n
gam
bara
n
n
y
ata ba
gaim
ana p
engam
bi
lan
keputu
s
a
n
dal
am
ope
rasi
si
stem
tenag
a. Meto
de
aliran
da
ya yan
g
dig
una
kan adal
ah m
e
tode
Ne
wton-Rap
s
on
yan
g
di
kom
b
inasi den
gan
optim
isasi p
e
rsam
aan
koordi
na
si d
an iter
a
s
i l
a
m
bda. Ne
wton-Ra
psho
n m
e
m
berikan
pen
yele
saia
n
aliran
da
ya
yan
g
jum
l
ah itera
s
i
sed
i
kit se
dan
g i
t
erasi
lam
b
d
a
m
e
m
berikan
pen
yele
saia
n
optim
isasi da
ric fung
si ob
yektif biaya pe
m
bangkitan.
1. PEN
DA
HU
LU
AN
Optimisa
si al
iran d
a
ya m
e
rup
a
kan
sal
ah satu ma
salah d
a
lam a
nalisa si
stem
tenaga
yang be
rpe
r
an pentin
g d
a
lam an
alisa
pere
n
canaa
n siste
m
ten
aga bai
k dal
am pen
gada
an
sistem
yang
baru m
aup
un p
enge
mb
anga
n
siste
m
yang tel
a
h ad
a. Opti
misa
si ali
r
an
daya
seb
agai
suat
u studi
si
ste
m
tenaga ya
ng mem
beri
k
an ba
nyak inf
o
rma
s
i yan
g
antara l
a
in b
e
rup
a
sud
u
t fasa te
gang
an tiap
b
u
s d
a
lam
sist
em, be
sar
da
ya pemba
ng
kitan dan
beb
a
n
aktif ma
upu
n
rea
k
tif pada
tiap bus d
a
n
informa
s
i lain. Aliran day
a dapat juga
dipakai untu
k
memp
erol
e
h
kon
d
isi a
w
al
pada a
nalisa
ke
stabila
n.
Dua la
ng
kah
utama perhi
tungan o
p
timisa
si alira
n
d
a
ya sistem t
enag
a listri
k
adala
h
perhitu
nga
n
aliran
daya
d
an o
p
timisa
si
biaya
o
p
e
r
a
s
i
pemb
ang
kit seb
agai pem
beri daya.
Bi
aya
baha
n ba
kar
adala
h
faktor utama dala
m
stasi
un pe
mbang
kit yan
g
mengg
una
kan b
aha
n b
a
ka
r
fosil perl
u
di
minimisasi m
e
lalui pem
be
bana
n ekono
mis.
2.
ALIR
AN
DA
YA DEN
G
A
N
NEWTO
N
-R
APSHO
N
Tiap bu
s dal
am si
stem te
naga li
strik m
e
libatk
an em
pat besaran
yaitu:
Daya nyata (P),
Daya rea
k
tif (Q), teg
ang
a
n
(V) d
an
su
dut fasa te
ga
ngan
(
). Dala
m penyele
s
ai
an aliran d
a
ya
dua da
ri emp
a
t besa
r
an di
atas ditentu
k
an besa
r
nya
dan dua si
sanya dihitung
selama p
r
o
s
es
penyele
s
aia
n
aliran daya.
Prosedu
r iterasi pe
nyele
s
a
i
an alira
n
day
a deng
an met
ode Newto
n
-Rap
sh
on ad
a
l
ah:
Untu
k bu
s-b
u
s
beb
an, dim
ana
sch
i
P
dan
sch
i
Q
ditentukan, be
sar tegan
gan
dan sudut fasa
disetin
g
sam
a
deng
an nila
i-nilai bu
s sl
a
ck yaitu
0
,
1
0
i
V
dan
0
,
0
0
i
. Untuk b
u
s te
gang
an
terko
n
trol ata
u
bus g
ene
rat
o
r, diman
a
i
V
dan
sch
i
P
ditentu
k
an
, sudut fasa
n
y
a diseting
sa
ma
deng
an sudut
fasa bu
s sla
c
k yaitu
0
,
0
0
i
.
Untu
k bu
s bu
s beb
an,
)
(
k
i
P
dan
)
(
k
i
Q
dihitung d
e
ngan p
e
rsam
aan
Evaluation Warning : The document was created with Spire.PDF for Python.
ISSN: 16
93-6
930
TELKOM
NIKA
Vol. 3, No. 3, Dese
mbe
r
2005 : 167 - 176
168
j
i
ij
n
j
ij
j
i
i
Y
V
V
P
cos
1
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...
..(1)
j
i
ij
n
j
ij
j
i
i
Y
V
V
Q
sin
1
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...
..(2)
dan
)
(
k
i
P
=
sch
i
P
-
)
(
k
i
P
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...
..(3)
)
(
k
i
Q
=
sch
i
Q
-
)
(
k
i
Q
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...
..(4)
1.
Bus terkont
ro
l teganga
n
)
(
k
i
P
dan
)
(
k
i
P
masin
g
-masin
g dihitu
ng den
gan p
e
rsama
an (1
)
dan (3
).
2.
Elemen-elem
en matri
ks ja
cobi
an (
J
1
, J
2
, J
3
, dan
J
4
) dari
V
J
J
J
J
Q
P
4
3
2
1
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...
..(5)
elemen
-ele
m
en
J
1
adalah
j
i
ij
i
j
ij
j
i
i
i
Y
V
V
P
sin
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...
..(6)
j
i
Y
V
V
P
j
i
ij
ij
j
i
j
i
sin
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...
..(7)
elemen
-ele
m
en J2 ad
alah
j
i
ij
i
j
ij
j
ii
ii
i
i
i
Y
V
Y
V
V
P
cos
cos
2
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...
..(8)
j
i
Y
V
V
P
j
i
ij
ij
i
j
i
cos
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...
..(9)
elemen
-ele
m
en
J
3
adalah
j
i
ij
i
j
ij
j
i
i
i
Y
V
V
Q
cos
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...(10)
j
i
Y
V
V
P
j
i
ij
ij
j
i
j
i
cos
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...(11)
elemen
-ele
m
en
J
4
adalah
j
i
ij
i
j
ij
j
ii
ii
i
i
i
Y
V
Y
V
V
Q
sin
sin
2
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...(12)
Evaluation Warning : The document was created with Spire.PDF for Python.
TELKOM
NIKA
ISSN:
1693-6
930
Sim
u
lasi Optim
a
si Aliran Da
ya….(S
ubi
yanto
)
.
169
j
i
Y
V
V
P
j
i
ij
ij
i
j
i
sin
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...(13)
3.
Selesai
k
a
n
p
e
rsama
an si
multan (5
)de
ngan fa
ktori
s
asi da
n elemi
nasi G
a
u
ssi
a
n
.
4.
Besa
r dan
su
dut fasa tega
ngan di
hitung
denga
n
)
(
)
(
)
1
(
k
i
k
i
k
i
V
V
V
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...(14)
)
(
)
(
)
1
(
k
i
k
i
k
i
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...(15)
5.
Proses dil
anj
utkan
samp
ai
)
(
)
(
toleransi
P
k
i
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...(16)
)
(
)
(
toleransi
Q
k
i
....
....
...
...
....
...
.
....
....
...
...
....
...
.
....
....
...
...(17)
deng
an,
sch
i
P
dan
sch
i
Q
:
daya aktif da
n rea
k
tif dala
m
peru
n
it
i
V
dan
i
:
besar da
n su
dut fasa tega
ngan b
u
s i.
ij
Y
dan
ij
:
besar da
n su
dut fasa admi
t
ansi bu
s i ke
bus j.
ii
Y
dan
ii
:
besar da
n su
dut fasa admi
t
ansi sendi
ri bus i.
3. PEMBEB
A
N
AN
EK
ON
O
M
IS
Sistem tena
g
a
listri
k
yan
g
melayani
sua
t
u beba
n list
r
i
k
de
nga
n me
ngab
aikan ru
gi pad
a
salu
ran tran
smisi dipe
rlihat
kan p
ada g
a
m
bar 1.
Gamba
r
1. Pe
mbang
kit The
r
mis Seb
anya
k
n Bero
pera
s
i Melayani B
eban mel
a
lui
jaring
an Tran
smisi.
Sistem terse
but terdiri da
ri n pembang
ki
t thermis yan
g
dihubu
ng
ka
n denga
n su
atu bus
tunggal yang mensuplai beban li
st
rik
sebesar
D mega watt. Masukan setiap unit pembangkit
ditunju
k
an
ol
eh
Ci, yaitu ti
ngkat biaya
baha
n b
a
kar
dari
unit
pem
bang
kit. Kelu
aran
setiap
u
n
it
adala
h
Pi, merup
a
kan da
ya listrik yan
g
diban
gkit
ka
n oleh
setiap
unit pemba
n
g
kit. Total tingkat
biaya dari sistem terseb
ut adalah ju
mlah
dari bi
aya bahan
bakar ma
sin
g
-ma
s
in
g un
it
pemba
ng
kit.
Fung
si biaya
masal
ah pe
m
beba
nan e
k
o
nomis di
defini
s
ikan seba
gai
berikut:
Evaluation Warning : The document was created with Spire.PDF for Python.
ISSN: 16
93-6
930
TELKOM
NIKA
Vol. 3, No. 3, Dese
mbe
r
2005 : 167 - 176
170
)
(
i
i
i
P
C
C
.......................................................
(18)
,
)
(
2
i
i
i
i
i
i
P
c
P
b
a
P
C
.......................................................
(19)
deng
an
C
i
(P
i
)
: biaya generator ke
-i
P
i
: daya output gene
rato
r ke
-i (MW).
a
i
, b
i
, c
i
: koefisie
n bia
y
a generator
ke-i.
Dalam Mi
nimi
sa
si total biaya prod
uksi terseb
ut haru
s
d
i
penu
hi 2 ke
kanga
n:
1) Keseim
ban
ga
n
Daya
i
i
L
P
P
D
0
.......................................................
(20)
deng
an D a
d
a
lah be
ban to
tal, dan P
L
adalah ru
gi tran
smisi m
a
si
ng-masin
g
dala
m
MW.
2)
Batas ma
ksi
mum dan min
i
mum Daya
(max)
(min)
i
i
i
P
P
P
.......................................................
(21)
deng
an
P
i(min)
: Daya pemba
ng
kita
n minimum p
e
mban
gkit ke
-i (M
W).
P
i(ma
x
)
: Daya pemba
ng
kita
n maksimum
pemba
ng
kit ke-i (M
W).
Operasi o
p
timal (biaya o
p
e
ra
si minimal
)
dapat di
rum
u
skan
seb
a
g
a
i beri
k
ut
1) Minimum
k
an
biaya
ope
ra
si
)
(
i
i
i
P
C
yaitu fungsi tujuan
2)
Mempe
r
timba
ngkan ke
ka
n
gan ke
setimb
anga
n
d
a
ya dan bata
s
p
e
m
bang
kitan
(Persamaa
n
20 dan 2
1
).
Pembeb
ana
n
optimal
sem
ua unit p
e
mb
ang
kit yang b
e
rop
e
rasi
aka
n
dicapai ji
ka
sem
ua
unit memiliki
pertam
bah
an
biaya baha
n bakar yang
sama, yaitu:
i
i
P
C
atau
1
1
P
C
,
2
2
P
C
,
3
3
P
C
, ….,
n
n
P
C
................................................
(22)
deng
an n ad
a
l
ah jumlah u
n
i
t pembang
kit
.
Selain de
nga
n penyam
aa
n pert
a
mba
h
an biaya
bah
an ba
ka
r (pe
r
sa
maa
n
koo
r
dina
si
)
adala
h
deng
a
n
metode iterasi lamb
da (p
ertamb
aha
n biaya bah
an
bakar).
1. Nilai
dari
λ
di
tentuka
n
terle
b
ih dah
ulu at
au nilai a
w
al (aca
k).
2.
Dihitun
g
kelu
aran
tiap
pe
mbang
kit
(P
i
) den
gan
pe
rsamaan
koo
r
d
i
nasi
da
n
rugi
tran
smi
s
i
(P
L
) den
gan f
o
rmul
a Kron’
s loss
i
i
i
c
b
P
2
...........................................................
(23)
n
i
n
j
n
i
i
i
j
ij
i
L
B
P
B
P
B
P
P
11
1
00
0
...........................................................
(24)
deng
an n = ju
mlah pem
ban
gkit
3. Jika
L
i
i
k
P
D
P
P
)
(
(nilai toleran
s
i yan
g
ditetapka
n
), maka pe
n
y
elesai
an
optimal dicap
a
i. Jika tida
k
kerj
akan lan
g
ka
h 4 beri
k
u
t
4. Update
nilai
λ
,
i
k
k
c
P
2
1
)
(
)
(
...........................................................
(25)
)
(
)
(
)
1
(
k
k
k
...........................................................
(26)
5.
Kemudia
n
ula
ngi 2.
Jika b
eba
n
unit ke-j
me
nca
pai
bata
s
P
j,maks
atau
P
i,min
maka
b
eban
bata
s
i
n
i a
k
an
meru
pa
kan b
eban tetap u
n
i
t j, beban sisanya
didi
strib
u
si
kan di
anta
r
a unit lainny
a.
Evaluation Warning : The document was created with Spire.PDF for Python.
TELKOM
NIKA
ISSN:
1693-6
930
Sim
u
lasi Optim
a
si Aliran Da
ya….(S
ubi
yanto
)
.
171
4.
SIMULASI OPTIMISASI ALIRAN
DAY
A
SISTEM TENAGA LISTRIK
Simulasi dila
ku
kan d
eng
a
n
baha
sa p
e
m
rog
r
am
an
MATLAB 5.3. Sistem tena
ga listri
k
yang
disimul
a
si
ka
n
terdi
r
i
dari
28 b
u
s
sep
r
ti pada
G
a
m
bar 2
be
riku
t. Bus1
seb
agai
bu
s
sla
c
k
tegang
an dite
ntuka
n
V
1
= 1,025
0
0
per unit.
Gamba
r
2. Di
agra
m
se
gari
s
si
stem tena
ga listri
k
26 b
u
s.
Den
gan d
a
ta beba
n dan d
a
t
a pemban
gki
t
an sep
e
rti Ta
bel 1 dan Ta
b
e
l 2beri
k
ut.
Tabel 1. Data
beban
D
a
ta
Be
ba
n
Bu
s Be
b
an
Bu
s
B
e
b
an
No. MW
MVar
No.
M
W
M
Var
1 51
41
15
70
31
2 22
15
16
55
27
3 64
50
17
78
38
4 25
10
18
153
67
5 50
30
19
75
15
6 76
29
20
48
27
7 0
0
21
46
23
8 0
0
22
45
22
9 89
50
23
25
12
10
0
0
24
54
27
11
25
15
25
28
13
12
89
48
26
40
20
13
31
15
27
28
13
14
24
12
28
40
20
Evaluation Warning : The document was created with Spire.PDF for Python.
ISSN: 16
93-6
930
TELKOM
NIKA
Vol. 3, No. 3, Dese
mbe
r
2005 : 167 - 176
172
Tabel 2. Data
pemban
gkita
n
Bus Mag.
Gen.
Batas
MVar
No.
Tegan
gan
MW
Min.
Max.
1 1,025
2 1,02
79
40
250
3 1,025
20
40
150
4 1,05
100
40
80
5 1,045
300
40
160
26
1,015
60
15
50
28
1,015
60
15
50
Biaya operasi
masing
-ma
s
i
ng gen
erato
r
dalam
$/h, Pi dalam M
W
a
dalah
seb
aga
i berikut:
,
0070
,
0
0
,
7
240
)
(
2
1
1
1
1
P
P
P
C
,
0095
,
0
0
,
10
200
)
(
2
2
2
2
2
P
P
P
C
,
0090
,
0
5
,
8
220
)
(
2
3
3
3
3
P
P
P
C
,
0090
,
0
0
,
11
200
)
(
2
4
4
4
4
P
P
P
C
,
0080
,
0
5
,
10
220
)
(
2
5
5
5
5
P
P
P
C
,
0075
,
0
0
,
12
190
)
(
2
26
26
26
26
P
P
P
C
,
0075
,
0
0
,
12
190
)
(
2
28
28
28
28
P
P
P
C
Batas Pemba
ngkita
n
masi
ng-m
a
si
ng pe
mbang
kit ada
lah se
bag
ai b
e
rikut (Tab
el 3):
Tabel 3. Bata
s pemb
ang
kit
an daya a
k
tif
Gen.
Min. MW
Max. MW
1 100
500
2 50
200
3 80
300
4 50
150
5 50
200
26
50
120
28
50
120
Kapasito
r
shu
n
t dan tap set
t
ing transfo
rm
ator adal
ah sebag
ai beri
k
u
t
:
Tabel 4. Kap
a
sitor
shu
n
t dan tap setting
transfo
rmato
r
Kapasito
r
Shunt
Tap
Tran
sformator
Bus No.
M
var
Pada
Tap Setting
1 4
2--3
0,96
4 2
2--13
0
,96
5 5
3--13
1
,017
6 2
4--8
1,05
11
1,5
4--12
1
,05
12
2
6--19
0
,95
15
0,5
7--9
0,95
19
5
Evaluation Warning : The document was created with Spire.PDF for Python.
TELKOM
NIKA
ISSN:
1693-6
930
Sim
u
lasi Optim
a
si Aliran Da
ya….(S
ubi
yanto
)
.
173
Data saluran
dan tran
sfo
r
mator teridi
ri dari rei
s
tan
s
i,
reakta
nsi d
a
n
seten
gah suse
ptan
si
kap
a
sitif dala
m
per unit pa
da ba
sis 1
00
MVA sebag
ai
berikut (Tab
el 5):
Tabel 5. Data
saluran da
n tran
sform
a
tor
Ha
sil p
e
rhitu
ngan
optimi
s
asi
alira
n
d
a
ya de
ngan
ba
ntuan
pro
g
ra
m Matlab
5.3
(poto
nga
n o
u
tput
ekse
ku
si) ad
alah seba
gi b
e
rikut:
Pen
y
elesaian Aliran Da
ya De
nga
n Metode Ne
wto
n
-Raphson
Ke
salahan Da
y
a
M
a
ximum = 8.050
1e-008
Pada Iterasi ke = 9
Bus Teganga
n Sudut -
---
--B
eban--
-
-
-
-
-
Pe
mbangkitan-- I
n
jeksi
No. Mag.
D
e
rajat MW
Mvar MW
Mvar Mvar
1 1.025 0
.
000 51.000
4
1
.000 727.7
19
224.087 4.
00
0
2 1.020
-0.
939 22.000
1
5
.000 79.00
0
124.609 0.
000
3 1.035
-4.
237 64.000
5
0
.000 20.00
0
61.420 0.0
0
0
4 1.050
-3.
614 25.000
1
0
.000 100.0
00
46.823 2.0
0
0
5 1.045 1
.
076 50.000
3
0
.000 300.0
00
124.333 5.
00
0
6 0.999
-2.
626 76.000
2
9
.000 0.00
0
0.000 2.00
0
7 0.994
-3.
238 0.000
0
.
000 0.000
0.000 0.00
0
8 0.997
-3.
329 0.000
0
.
000 0.000
0.000 0.00
0
9 1.009
-5.
429 89.000
5
0
.000 0.00
0
0.000 3.00
0
10 0.990
-5.
603 0.000
0
.
000 0.000
0.000 0.00
0
11 0.997
-3.
264 25.000
1
5
.000 0.00
0
0.000 1.50
0
12 0.993
-4.
727 89.000
4
8
.000 0.00
0
0.000 2.00
0
13 1.014
-4.
457 31.000
1
5
.000 0.00
0
0.000 0.00
0
14 1.000
-5.
071 24.000
1
2
.000 0.00
0
0.000 0.00
0
15 0.992
-5.
575 70.000
3
1
.000 0.00
0
0.000 0.50
0
16 0.983
-5.
944 55.000
2
7
.000 0.00
0
0.000 0.00
0
17 0.987
-5.
143 78.000
3
8
.000 0.00
0
0.000 0.00
0
18 1.007
-1.
899 153.00
0
67.000 0.0
00
0.000 0.00
0
19 1.004
-6.
454 75.000
1
5
.000 0.00
0
0.000 5.00
0
20 0.982
-6.
078 48.000
2
7
.000 0.00
0
0.000 0.00
0
21 0.977
-5.
977 46.000
2
3
.000 0.00
0
0.000 0.00
0
22 0.979
-6.
489 45.000
2
2
.000 0.00
0
0.000 0.00
0
23 0.977
-7.
141 25.000
1
2
.000 0.00
0
0.000 0.00
0
24 0.969
-7.
436 54.000
2
7
.000 0.00
0
0.000 0.00
0
25 0.975
-6.
825 28.000
1
3
.000 0.00
0
0.000 0.00
0
26 1.015
-1.
836 40.000
2
0
.000 60.00
0
32.640 0.0
0
0
Evaluation Warning : The document was created with Spire.PDF for Python.
ISSN: 16
93-6
930
TELKOM
NIKA
Vol. 3, No. 3, Dese
mbe
r
2005 : 167 - 176
174
27 0.975
-6.
530 28.000
1
3
.000 0.00
0
0.000 0.00
0
28 0.985
-6.
021 40.000
2
0
.000 60.00
0
39.525 0.0
0
0
Total
1331.000 6
70.
000 1346.71
9
653.437 25.
00
0
B =
0.0015 0.0
0
16 0.0010
-0.
0000 -0.0
003 -
0
.0001 -0
.0012
0.0016 0.0
0
45 0.0054
0.
0002 -0.0
007 -
0
.0001 -0
.0020
0.0010 0.0
0
54 0.0309
0.
0000 -0.0
017 -
0
.0013 -0
.0032
-0.0000
0.00
02 0.0000
0.
0028 -0.0
005 -
0
.0008 -0
.0003
-0.0003
-0.00
07 -0.001
7 -0.
0005 0.008
6 -
0
.0000 -0
.0013
-0.0001
-0.00
01 -0.001
3 -0.
0008 -0.0
000
0.0178 -0
.0020
-0.0012
-0.00
20 -0.003
2 -0.
0003 -0.0
013 -
0
.0020 0.01
66
B0 =
-0.0002
-0.00
06 0.0067
0.
0000 0.000
1 -
0
.0011 -0
.0006
B00 =
0.0055
Total rugi sistem = 15.7121 MW
Total Bia
y
a Pem
bangkitan = 178
37.949 $/h
Incremental bia
y
a pengiriman da
ya
(sistem lambda) = 13.69805
4 $/
MWh
Pen
y
alura
n
Da
ya
Optimal :
460.5427
164.2100
187.1736
146.5788
186.4737
95.2026
114.6383
Nilai absolut kesalahan da
ya n
y
at
a pada bus slack, dpslack = 2.6718 pu
Setelah itera
s
i ke-4 p
ada p
r
ose
s
optimi
s
a
s
i (pe
m
be
ban
an ekonomi
s
) diperol
eh :
B =
0.0018 0.0
0
13 0.0008
-0.
0000 -0.0
005 -
0
.0001 -0
.0012
0.0013 0.0
0
14 0.0010
0.
0001 -0.0
005 -
0
.0000 -0
.0010
0.0008 0.0
0
10 0.0032
0.
0000 -0.0
009 -
0
.0005 -0
.0011
-0.0000
0.00
01 0.0000
0.
0024 -0.0
005 -
0
.0007 -0
.0003
-0.0005
-0.00
05 -0.000
9 -0.
0005 0.013
4 -
0
.0001 -0
.0013
-0.0001
-0.00
00 -0.000
5 -0.
0007 -0.0
001
0.0154 -0
.0016
-0.0012
-0.00
10 -0.001
1 -0.
0003 -0.0
013 -
0
.0016 0.01
25
B0 =
1.0e-003 *
-0.3053
-0.08
94 0.7354
0.
0419 0.399
9 -
0
.6373 -0
.2544
B00 =
0.0055
Total rugi sistem = 12.3531 MW
Incremental bia
y
a pengiriman da
ya
(sistem lambda) = 13.44765
7 $/
MWh
Pen
y
alura
n
Da
ya
Optimal :
441.4958
168.6529
258.6314
Evaluation Warning : The document was created with Spire.PDF for Python.
TELKOM
NIKA
ISSN:
1693-6
930
Sim
u
lasi Optim
a
si Aliran Da
ya….(S
ubi
yanto
)
.
175
133.4328
159.8849
82.2755
98.9705
Nilai absolut kesalahan da
ya n
y
at
a pada bus slack, dpslack = 0.0006 pu
Pen
y
elesaian Aliran Da
y
a
De
ngan Met
ode Ne
w
t
o
n
-Ra
p
h
son
Ke
salahan Da
y
a
M
a
ximum = 2.210
04e-005
Pada Iterasi ke = 2
Bus Teganga
n Sudut -
---
--B
eban--
-
-
-
-
-
Pe
mbangkitan-- I
n
jeksi
No. Mag.
D
e
rajat MW
Mvar MW
Mvar Mvar
1 1.025 0
.
000 51.000
4
1
.000 441.4
32
251.850 4.
00
0
2 1.020
-0.
193 22.000
1
5
.000 168.5
55
57.198 0.0
0
0
3 1.045
-0.
606 64.000
5
0
.000 258.5
16
78.682 0.0
0
0
4 1.050
-2.
011 25.000
1
0
.000 133.3
46
32.499 2.0
0
0
5 1.045
-1.
498 50.000
3
0
.000 160.4
32
143.741 5.
00
0
6 1.001
-2.
824 76.000
2
9
.000 0.00
0
0.000 2.00
0
7 0.995
-2.
320 0.000
0
.
000 0.000
0.000 0.00
0
8 0.998
-2.
187 0.000
0
.
000 0.000
0.000 0.00
0
9 1.011
-4.
224 89.000
5
0
.000 0.00
0
0.000 3.00
0
10 0.991
-4.
056 0.000
0
.
000 0.000
0.000 0.00
0
11 0.998
-2.
844 25.000
1
5
.000 0.00
0
0.000 1.50
0
12 0.994
-3.
131 89.000
4
8
.000 0.00
0
0.000 2.00
0
13 1.022
-1.
211 31.000
1
5
.000 0.00
0
0.000 0.00
0
14 1.008
-2.
357 24.000
1
2
.000 0.00
0
0.000 0.00
0
15 0.999
-3.
115 70.000
3
1
.000 0.00
0
0.000 0.50
0
16 0.990
-3.
816 55.000
2
7
.000 0.00
0
0.000 0.00
0
17 0.982
-4.
293 78.000
3
8
.000 0.00
0
0.000 0.00
0
18 1.007
-1.
860 153.00
0
67.000 0.0
00
0.000 0.00
0
19 1.005
-5.
924 75.000
1
5
.000 0.00
0
0.000 5.00
0
20 0.984
-4.
358 48.000
2
7
.000 0.00
0
0.000 0.00
0
21 0.977
-5.
090 46.000
2
3
.000 0.00
0
0.000 0.00
0
22 0.980
-5.
028 45.000
2
2
.000 0.00
0
0.000 0.00
0
23 0.978
-6.
182 25.000
1
2
.000 0.00
0
0.000 0.00
0
24 0.969
-6.
415 54.000
2
7
.000 0.00
0
0.000 0.00
0
25 0.975
-6.
122 28.000
1
3
.000 0.00
0
0.000 0.00
0
26 1.015
-0.
471 40.000
2
0
.000 82.20
6
28.505 0.0
0
0
27 0.976
-4.
488 28.000
1
3
.000 0.00
0
0.000 0.00
0
28 0.985
-3.
449 40.000
2
0
.000 98.87
9
27.825 0.0
0
0
Total
1331.000 6
70.
000 1343.36
6
620.299 25.
00
0
Total Bia
y
a Pem
bangkitan = 164
82.508 $/h
Pada kondi
si
awal op
erasi total biaya
pemba
ng
kita
n adala
h
178
37.949 $/h d
an total
biaya pemb
a
ngkita
n
setel
ah dila
ku
kan
optim
isa
s
i adala
h
1648
2.508 $/h. b
e
ra
rti tiap jam
diperoleh p
e
n
gura
nga
n bia
y
a $1355.4
4
1
.
Sehingga ki
ra-ki
r
a d
a
lam
satu tahu
n
24 h x 365 x
1355.4
41 $/h
= $ 1187
366
3.160.
5. PENUT
UP
Dari Pem
bah
asa
n
diata
s
d
i
samp
aikan b
eberapa
ke
si
mpulan
seb
a
gai beri
k
ut:
1.
Dari
ha
sil
si
mulasi
diata
s
28
bu
s
dipe
roleh
keuntu
ngan
$
118
7
3663.1
60
seti
ap tah
un ji
ka
dilakukan o
p
timisa
si ali
r
an
daya d
a
rip
a
da tida
k. Se
hingg
a jika S
i
stem inte
rko
neksi
Ja
wa-
Bali-NTB ya
ng terdiri
da
ri baya
k
b
u
s
a
k
an
lebih
iperoleh
ke
untung
an ji
ka dila
ku
kan
optimisasi.
2.
Dalam
studi
anali
s
is
siste
m
tenaga tid
a
k
cu
kup h
a
n
y
a denga
n ali
r
an d
a
ya, perlu dilakukan
optimisa
s
i u
n
t
uk men
e
kan
biaya ope
ra
si si
stem ten
aga d
an dip
e
role
h keunt
unga
n yang
lebih be
sa
r. Penerapa
n JST Hopfield
untuk p
r
o
s
e
s
optimisa
s
i u
n
tuk pe
nelitia
n lebih lanj
ut
diupaya
k
an
selai
n
men
d
apatkan bi
aya ope
ra
si yang mini
mal
juga
ketep
a
tan anta
r
a
keb
u
tuha
n da
ya denga
n total pemba
ng
ki
tan.
Evaluation Warning : The document was created with Spire.PDF for Python.
ISSN: 16
93-6
930
TELKOM
NIKA
Vol. 3, No. 3, Dese
mbe
r
2005 : 167 - 176
176
Tulisa
n
ini h
anya suatu g
a
mba
r
an
ke
cil optim
isa
s
i
sistem
tenag
a. Untu
k si
st
em yang
besar dan
p
e
rma
s
al
ahan
yang l
ebih
kompl
e
ks
se
perti koo
r
di
n
a
si pemb
ang
kit
hid
r
o
-
thermal
haru
s
dila
ku
kan penye
k
le
saian
Unit Comm
itm
ent
.
DAF
TA
R PU
STAK
A
[1]
Hadi Saa
dat, 1999 , “
Pow
e
r System An
alysis”
M
c
G
r
aw-Hill, Singapore.
[2]
Muhamm
ad
Arif Adha, 1996,
“Pro
gr
am Simulasi Operasi e
k
onomis Pe
mbangkit
Deng
an M
e
tode Progra
m
Dinamis”
Tugas
Ak
hir, UNDIP, Semarang.
[3] Subiyanto,
2
003,
“Op
t
imisasi Pemba
ngkitan
Ten
a
ga Lis
t
rik
M
e
nggunaka
n Jaring
a
n
Syaraf Tiru
a
n
Hopfield
Adaptif”
Tes
i
s
UGM, Yogyak
arta.
[4] Sulasn
o,
199
3,
“
A
nalisis Sistem Tenaga Listrik”
Satya Wacana, Semarang.
[5]
Unggul Satriatama, 1996,
“Progr
am
Simulasi An
alisis Aliran
Da
ya
Den
g
an M
e
tode
New
t
on Rap
h
son
”
Tuga
s Akhi
r, UNDI
P, Semarang.
[6]
Kamal Idris,
Wlliam
D. Stevenson,
Jr, 1996,
“Analisis Si
stem
Tena
g
a
Lis
t
rik”
Erlang
ga, Ja
karta.
[7]
Woo
d
A. J. and Woll
enb
erg B. F. , 1994,
“Power
G
e
nera
tion, O
p
eratio
n and
Control”
John Willey & Sons, Singapore.
Evaluation Warning : The document was created with Spire.PDF for Python.