I nte rna t io na l J o urna l o f   E lect rica l a nd   Co m p ute E ng in ee ring   ( I J E CE )   Vo l.   8 ,   No .   2 A p r il   201 8 ,   p p .   741 ~ 7 4 8   I SS N:  2088 - 8708 DOI : 1 0 . 1 1 5 9 1 / i j ec e . v 8 i 2 . p p 7 4 1 - 7 4 8           741       J o ur na l ho m ep a g e h ttp : //ia e s co r e . co m/ jo u r n a ls /in d ex . p h p / I JE C E   Rev iew   o u nder  F requen cy  Lo a d Shedding  P ro g ra m   o f   K o so v o   Po w er Sys te m   b a sed  o n ENTS O - Require m e nts       G a z m en d K a ba s hi 1 S k e nd e K a ba s hi 2   1 KO S TT - Ko so v o   T ra n s m issio n   S y ste m   a n d   M a rk e Op e ra to r,   P risti n a ,   Ko so v o   2 F a c u lt y   o f   M a th e m a ti c a a n d   Na tu ra S c ien c e s,  M o th e T h e re sa   S tr .   Nr.  5 ,   Un iv e rsity   o f   P rish ti n a ,   P r ish ti n a ,   K o so v o       Art icle  I nfo     AB ST RAC T   A r ticle  his to r y:   R ec eiv ed   No v   1 0 ,   2 0 1 7   R ev i s ed   Ma r   1 2 ,   2 0 1 8   A cc ep ted   Ma r   2 0 ,   2 0 1 8     Un d e r - f re q u e n c y   lo a d   sh e d d in g   (UFL S is  d e sig n e d   to   p r o tec t h e   p o w e r   s y ste m   w h e n   th e   f re q u e n c y   d ro p b e lo w   g iv e n   th re sh o l d b y   s w i tch in g   o f c e rtain   a m o u n ts  o f   th e   lo a d   a im in g   th u to   b a lan c e   g e n e ra ti o n   a n d   lo a d .   T h is   p a p e p re se n ts  a   re v ie w   o f   th e   e x isti n g   UF L S   (Un d e F re q u e n c y   L o a d   S h e d d in g p r o g ra m   in   c o m p li a n c e   w it h   re c e n tl y   re v ise d   P o li c e - 5   o Op e ra ti o n a Ha n d b o o k   o f   ENT S O - e .   T h e   p ro p o se d   re v iew   o f   th e   c u rre n t   UFL S   p ro g ra m   f o Ko s o v o   P o w e S y ste m   h a c o n sid e re d   th e   m a in   sta n d a rd re q u irem e n ts  a n d   g u id e li n e f o UFL S   se b y   EN T S O - E.   T h is  w o rk   e x a m in e   s y ste m   p e r f o r m a n c e   b y   c o n d u c ti n g   d y n a m ic  si m u latio n o f   UF LS   sc h e m e s   su b jec t o   d if fe re n t   im b a lan c e b e tw e e n   lo a d   a n d   g e n e ra ti o n ,   a n d   i n c lu d e s   th re e   p o w e s y st e m   islan d   m o d e     sc e n a rio w it h   d if f e r e n e q u iv a len i n e rti a   o f   th e   sy ste m ,   re sp e c ti v e l y   d iff e r e n si z e   o f   th e   sy ste m s.  W it h   a im   to   d e f in e   th e   b e st p r o g ra m   o f   UF L S ,   w h ich   f it s to   th e   Ko s o v o   P o w e S y ste m   f re q u e n c y   b e h a v io r,   tw o   d if f e r e n UFL S   p ro g ra m a re   a n a l y z e d   a n d   re su lt a re   c o m p a re d .   T h e   p ro p o se d   p r o g ra m   is  tes ted   u sin g   a   larg e   sc a le  P S S /E   m o d e l   w h ich   re p re se n ts  in terc o n n e c ted   p o w e s y ste m   a re a   o f   S o u th e a st E u ro p e .   K ey w o r d :   E NT SO -   Fre q u en c y   P SS /E  d y n a m ic  s i m u latio n   S y s te m   i n er tia    UFLS    Co p y rig h ©   2 0 1 8   In stit u te o A d v a n c e d   E n g i n e e rin g   a n d   S c ien c e   Al rig h ts  re se rv e d .   C o r r e s p o nd ing   A uth o r :   Sk e n d er   Kab ash i,    Facu lt y   o f   Ma t h e m atica l a n d   Natu r al  Scie n ce s ,     Mo th er   T h er esa  Str.   Nr .   5 ,     Un i v er s it y   o f   P r is h t in a,     P r is h tin a,   Ko s o v o .   E m ail: s k e n d er . k ab ash i@ u n i - p r . ed u       1.   I NT RO D UCT I O N     Un d er - f r eq u e n c y   lo ad   s h ed d i n g   ( U F L S)  i s   d ef i n ed   as  c o o r d in ated   s et  o f   co n tr o ls   u s in g   u n d er   f r eq u en c y   r ela y s   w h ic h   r esu l ts   in   th d ec r ea s o f   elec tr ical  lo ad s   in   th p o w er   s y s te m ,   w it h   ai m   to   r ec o v er   th s y s te m   f r eq u en c y .   L o ad   s h ed d in g   as   th e   las r eso r to   av o id   m aj o r   p o w er   s y s te m   b r e ak d o w n   h as   b ee n   u tili ze d   f o r   lo n g   t i m e.   I t   is   m ai n l y   tr i g g er ed   b y   u n d er - f r eq u en c y   o r   u n d er - v o ltag e   p r o tectio n   r ela y s   a n d   ac tu ated   b y   d i s tr ib u tio n   s y s t e m   cir c u it  b r ea k er s .   P r o p er   d esig n   o f   lo ad   s h ed d i n g   s c h e m es  w h ich   in cl u d p r o p er   s ettin g s   o f   u n d er - f r eq u en c y   p r o tectio n   is   m o s r ele v a n is s u to   en s u r s m o o th   lo ad   r elief ,   in   s i tu at io n s   w h er th p o w er   s y s te m   o t h e r w i s w o u ld   g o   u n s tab le.   T h cu r r en r ev is ed   p r o g r a m   r eq u ir e m en ts   f o r   UF L S   ar p r esen ted   in   Op er atio n al  H an d b o o k   o f   E NT SO - E /P o lice - 5 .     E ac h   T SO sh all  i m p le m en t t h e   E NT SO - E   R C E   g e n er al  U FLS sc h e m a s   f o llo w ed :   Fre q u en c y   in   t h r an g 4 9 . 0   to   4 8 . 0   Hz:   a.   A t   leas a n   a m o u n t   o f   d e m a n d   co r r esp o n d in g   to   5 o f   t h e   to tal  lo ad   s h all   b d is co n n ec ted   at    4 9 . 0   Hz.   b.   I n   to tal,   an   a m o u n o f   d e m an d   co r r esp o n d in g   to   4 5 +/ -   7 o f   th to tal  lo ad   s h all  b d is co n n ec ted   b et w ee n   4 9 . 0   an d   4 8 . 0   Hz.     Evaluation Warning : The document was created with Spire.PDF for Python.
                        I SS N :   2 0 8 8 - 8708   I n t J   E lec  &   C o m p   E n g ,   Vo l.  8 ,   No .   2 A p r il   201 8   :   7 4 1     7 4 8   742   T h UFL S sc h e m s h all  b i m p le m en ted   s tep w is ta k in g   i n t o   ac co u n t f o llo w i n g   co n s id er a tio n s :   a.   T h n u m b er   o f   d is co n n ec t io n   s tep s   s h all  b m i n i m u m   6   ( i n c lu d in g   t h s tep   tr ig g er ed   at  4 9 . 0   Hz) ,   b.   Fo r   ea ch   s tep ,   an   am o u n o f   d e m an d   co r r esp o n d in g   to   1 0 o f   to tal  lo ad   s h all  b d is co n n e cted   at   m ax i m u m .   c.   A d d itio n al  d f /d t f u n ctio n   i n   U FLS r ela y s   i s   allo w ed   in   t h r an g 4 9 . 8     4 9 . 0   Hz.   d.   No   in ten t io n al  ti m d ela y   s h al l b s et  in   UF L S r ela y s .   e.   Ma x i m u m   d is co n n ec t io n   d ela y   s h all  b 1 5 0   m s   i n cl u d in g   b r ea k er s   o p er atio n   ti m [ 1 ] .   T h ex is ti n g   U F L p r o g r a m   w h ic h   i s   o p er atio n al  i n   Ko s o v o   P o w er   S y s te m ,   w as  e s tab li s h ed   b ased   o n   p r ev io u s   r ec o m m e n d atio n   o f   E NT SO - E ,   i n cl u d in g   o n l y   f o u r   s tep s   w it h   d i f f er e n a m o u n o f   lo ad   d is co n n ec t io n s   a n d   w ith   ad d iti o n al  ti m d ela y   in c lu d ed   in   u n d er - f r eq u e n c y   r ela y s   [ 2 ] .   T h m ai n   f ac to r s   th at  E NT SO - E   h a s   in i tiated   th U FLS r ev i e w   ar e:   a.   S y s te m   i m p le m e n tatio n   e n s u r es  th e   af f ec ti v it y   o f   UF L S:  it   m ea n s   m i n i m al  n ec es s ar y   s h ed d in g   o f   lo ad ,   b.   C o m p en s ate  d is co n n ec tio n   o f   d is p er s ed   g en er atio n   at  u n f a v o r ab le  f r eq u en cies,   c.   Av o id   o v er   f r eq u en c y   ( o v er co m p e n s atio n ) ,   o v er v o lta g an d   p o w er   tr an s ie n ts   t h at  ca n   lead   to   an   ad d itio n al  lo s s   o f   g en er atio n .     1 . 1 .   T heo re t ica B a ck g ro un d o f   un der  F re qu e ncy   L o a S hed din g   P ro t ec t io n ( UF L S)   E ac h   p ar o f   th p o w er   s y s te m   ca n   b u n b ala n ce d   if   t h lo ad   ex ce ed   th g e n er atio n .   S u c h   ch a n g e   b et w ee n   t h g e n er ated   p o w er   an d   th lo ad   p o w er   i s   k n o w n   as  g e n er ato r   i m b ala n ce .   I n   s u ch   c ase s   g e n er ato r s   au to m at icall y   w ill  in cr ea s e   th m ec h an ica p o w er   o n   th t u r b in t h r o u g h   th p r i m ar y   r eg u la tio n   ( s p i n n in g   r eser v e,   if   an y ) ,   s o   th at  th f r eq u en c y   p o in w o u ld   b k ep n ea r   to   th n o m i n al  v a lu o f   5 0   Hz.   Ho w ev er ,   i f   th g en er ato r s   r ea c h   t h r e g u l atio n   li m it s   w h er ea s   t h s y s te m   r e m ain s   i n   t h i s lan d   m o d o f   w o r k   t h e n   t h e   r ap id   f r eq u en c y   d ec li n w i ll  o cc u r .   I n   t h ese   ca s es   i f   t h u n d er - f r eq u en c y   p r o tectio n s   ar n o ac ti v ated   f o r   au to m at ic  s w itc h   o f f   o f   t h lo ad ,   th en   u n d er - f r eq u e n c y   p r o tectio n   ( 4 7 . 5   Hz)   o f   g e n er ato r s   w ill   ca r r y   o u t h eir   tr ip p in g   o f f ,   av o id i n g   p o s s ib le   d am a g o n   t h g e n er ato r   s y s t e m .     T h w o r k   o f   g e n er ato r s   w it h   t h f r eq u e n c y   b elo w   t h cr itica v alu e s   af f ec t s   n e g ati v el y   th g en er ato r .   Su c h   ef f ec in   t h ca s es  o f   h y d r o - g en er ato r s ,   e v e n   w ith   d r o p   u p   t o   1 0 %,  is   n o as  g r ea as  in   th ca s es  o f   th er m a g e n er ato r s   as  t h latt er   ar s en s it iv e v e n   f o r   t h v alu e s   ab o v 5 %,  as  is   il lu s tr ated   in   F i g u r 1 .   I n   o r d er   to   av o id   to tal  d ec lin o f   th s y s te m s   is o lated   p ar t,  t h u n d er - f r eq u e n c y   p r o tectio n s   a r ap p lied   in   ev er y   p o w er   s y s te m   i n   o r d er   to   ac tiv ate  au to m atic  lo ad - s h ed d in g   i. e.   b y   s w itc h i n g   o f f   ce r tai n   a m o u n t s   o f   th p o w er   s y s te m   lo ad     ai m in g   t h u s   to   b alan ce   g e n er atio n   a n d   lo ad   [ 3 ] .           Fig u r 1 . Stea m   t u r b in o f f - f r eq u en c y   li m its       T h d y n a m ic  b e h av io u r   o f   a n   in d iv id u al  s y n ch r o n o u s   t u r b in e - g e n er ato r     ca n   b d escr ib ed   u s i n g   t h s w i n g   eq u atio n   o f   r o tatin g   m as s :     ) ( 2 2 ei mi ni n i i P P S f dt df H                 ( 1 )     w h er   is   t h f r eq u e n c y   o f   g en er ato r   n   is   th n o m i n al  f r eq u en c y ,   m   is   t h m ec h a n ica p o w er   o f   tu r b in e - g en er ato r   ,   an d   e   is   t h elec t r ical  p o w er   o f   g e n er ato r     [ 4 ] .     5 0 4 9 4 8 4 7 5 1 5 2 F r e q u e n c y   ( H z ) 1 s 6 s 6 0 s 1 0 m i n 1 0 0 m i n t i m e 4 6 5 3 C o n t i n o u s   o p e r a t i o n Evaluation Warning : The document was created with Spire.PDF for Python.
I n t J   E lec  &   C o m p   E n g     I SS N:  2088 - 8708       R ev iew   o f u n d er F r eq u en cy   L o a d   S h ed d i n g   P r o g r a o f Ko s o vo   P o w er … .   ( Ga z men d   K a b a s h i)   743   T h is   n o n - li n ea r   d i f f er e n tial  eq u atio n   d e f i n es  t h f r eq u e n c y   ch an g d ep en d i n g   o n   t h d i f f er en ce   b et w ee n   th e   elec tr ical  an d   m ec h a n ical  p o wer s ,   r esp ec tiv el y ,   f r o m   t h u n b alan ce   b et w ee n   th g e n er atio n   an d   th lo ad .   T h e   r ate  o f   ch an g o f   f r eq u e n c y   i s   d ef in ed   b y   t h i m b ala n ce   an d   in er tia  o f   t h tu r b i n e - g e n er ato r .     T h in er tia  co n s ta n i s   g i v en   in   s ec o n d s   a n d   it   ca n   b i n ter p r eted   as  th e   ti m t h at  e n er g y   s to r ed   in   r o tatin g   p ar ts   o f   tu r b in e - g e n er ato r   is   ab le  to   s u p p l y   lo ad   eq u al  to   th r ated   ap p a r en p o w er   o f   t h tu r b in e - g en er ato r .   S y s te m   in er t ia  co n s tan p la y s   s ig n i f ica n r o le  o n   f r eq u e n c y   b eh a v io u r   af ter   d is tu r b an ce   o cc u r s   i n   p o w er   s y s te m .   T h g r ea ter   th e   in e r tia,   th le s s   ac ce ler atio n   w il l b o b s er v ed   an d   th les s   w il l b th f r eq u en c y   d ev iatio n .   T h i n er tia  co n s tan t     d escr ib es t h i n er tia  o f   a n   i n d iv id u al  t u r b in e - g e n er ato r   an d   d if f er s   f r o m   t y p e   o f   g e n er ato r s .   I n er tia  co n s ta n f o r   in d iv id u al  g e n er ato r s   is   r ep r esen ted   b y   eq u a tio n :     n n S J H 2 2 1                     ( 2 )     W h er e:     J   is   th m o m en t o f   in er tia  o f   g en er ato r   an d   tu r b i n [ k g · m 2 ] ,     n   is   t h r ated   m ec h a n ical  a n g u lar   v elo cit y   o f   t h r o to r   [ r ad / s ] ,     n   is   th r ated   ap p ar en t p o w er   o f   th g en er ato r   [ V A ].     R ea i n ter co n n ec ted   p o w er   s y s te m s   u s u all y   o p er ates  i n   l ar g s ca le   s y n ch r o n o u s   zo n e,   s u c h   a   E NT SO - E   an d   co n s is t   o f   l ar g n u m b er   o f   g en er ato r s   a n d   th n e t w o r k   w h ic h   co n n ec t s   th e   g e n er ato r s   to   lo ad s .   Hen ce ,   w h e n   an   i m b ala n ce   in   t h s y s te m   ar is e s ,   f r eq u e n c y   i s   n o u n i f o r m   t h r o u g h o u t th s y s te m .     T h to tal  eq u iv alen s y s te m   in er tia  co n s tan t is  g iv e n   b y :     n ni n ni i nn n n nn n n n s y s S S H S S S S H S H S H H 1 1 1 1 2 1 2 2 1 1 ....... .......           ( 3 )       W h er e:   S ni   is   t h r ated   ap p ar en t p o w er   o f   th g e n er ato r   i   an d   H i   is   t h in er tia  co n s ta n t o f   tu r b in e - g e n er ato r   i   Fo r   ex a m p le  f o r   Ko s o v o   P o w er   S y s te m   th to tal  eq u i v al en s y s te m   in er tia  co n s tan is H sys =2 . 6   s ec o n d s   T h eq u iv ale n s y s te m   i n er tia  co n s ta n v ar iate  b ased   o n   th e   n u m b er   o f   g e n er atio n s   i n   o p er at io n   a n d   th eir   c h ar ac ter is tic s .   T h lo wer   th i n er tia,   t h h i g h   ac ce le r atio n   w ill   b o b s er v ed   an d   t h h ig h   w ill  b t h e   f r eq u en c y   d ev ia tio n .     Fo r   w i n d   an d   s o lar   s o u r ce s   n o   ef f ec ts   o f   t h eir   i n er tia  ar co n s id er ed   s in ce   n o n - s y n c h r o n o u s l y   co n n ec ted   g e n er atio n   u n its ,   li k m o d er n   w i n d   tu r b in g en er ato r s ,   P s o lar   u n it s   ar c o n n ec ted   v ia  p o w er   co n v er ter s .   Fo r   w i n d   t u r b in es  th eir   r o tatio n al  s p ee d   is   i s o lated   f r o m   t h s y s te m   f r eq u e n c y .   T h ey   t h er e f o r d o   n o d eliv er   n at u r al  i n er tial   r esp o n s a n d   d o   n o co n tr i b u te  to   th e   s y s te m   i n er tia.   Mo r R E S,  an d   le s s   co n v e n tio n al  g en er atio n   th s y s te m   w ill  lo s g r ad u a ll y   t h i n er tia  w h ic h   w ill  h a v n e g ati v i m p ac o n   s y s te m   s tab ilit y .   T h s y s te m   g e n er atio n   i m b ala n ce   is   ca lc u lated   u s in g   eq u at io n     100 ) %( l o s s e s l o a d l o s s e s l o a d P P G P P i m b a l a c e             ( 4 )     I is   v er y   d if f ic u lt  to   d e f i n th m o s p r o b ab le  lev el  o f   s y s te m   g e n er atio n   i m b ala n ce   to   b e   en co u n ter ed   p ar ticu lar l y   w it h   h i g h   m e s h ed   s y s te m s   h a v i n g   s tr o n g   in ter co n n ec t io n s   as  ar th i n d iv id u al   p o w er   s y s te m s   o f   SEE .   T h is   ap p lies   p ar ticu lar l y   i n   cr itical   ca s es  s u ch   a s   s u d d en   d is co n n ec tio n s   f r o m   th e   in ter co n n ec tio n .   I n   t h ese  e v e n ts   ce r tai n   ele m e n ts   o f   t h s y s te m   ar o v er lo ad ed   to   th cr itical  v al u e s ,   w h ich   d ep en d   m a in l y   o n   t h le v el  o f   s y s te m   o v er lo ad   at  th ti m e   o f   o cc u r r en ce .   I f   s u ch   o v er l o ad   elem e n v alu e s   ex ce ed   th r a n g e   o f   n o r m al  o v er lo ad   lev els  t h e n   t h cu r r en v alu e s   e n ter s   i n   t h r an g o f   o v er - c u r r en t   p r o tectio n   ac tu atio n   ( I >> )   w it h   a   ti m d ela y   b ei n g   s h o r ter   t h an   t h ti m e   o f   f r eq u e n c y   d ec lin in   t h le v el  o f   u n d er - f r eq u e n c y   p r o tectio n s   a ctu atio n .   T h er ef o r th a n al y s is   i s   f o c u s ed   in   ca s es  w h er e   s y s te m   is   e n tire l y   is o lated   f r o m   t h r est o f   s y n ch r o n o u s   zo n e.   Evaluation Warning : The document was created with Spire.PDF for Python.
                        I SS N :   2 0 8 8 - 8708   I n t J   E lec  &   C o m p   E n g ,   Vo l.  8 ,   No .   2 A p r il   201 8   :   7 4 1     7 4 8   744   T h n atu r o f   th s y s te m   f r eq u en c y   c h a n g e s   w as  a n al y ze d   in   ter m s   o f   t h s y s te m   g en er ato r   i m b alan ce   le v el.   Fi g u r 2   p r esen t s   f o u r   ca s es  o f   f r eq u en c y   d ec li n e s   w ith   d i f f er en g en er atio n   i m b ala n ce s   w it h o u t a c tin g   o f   U F L S.           Fig u r 2 .   S y s te m   f r eq u e n c y   r e s p o n s f r o m   d if f er en v al u o f   g en er ato r   i m b ala n ce s         2.   SE L E C T I O O F   T H E   P O WE SY ST E M   I SL AND  M O DE   O P E RA T I O R E G A RDIN G   UF L S   I n   o r d er   to   d ef in th UF L p r o g r a m   f o r   Ko s o v o   P o w er   S y s te m ,   b ased   o n   E NT SO - E   r eq u ir e m e n t s   an d   r ec o m m e n d atio n s ,   t h r ee   d if f er en s i m u latio n   ca s es  w er e   cr ea ted   s o   th at  th b eh av io r   o f   f r eq u e n c y   ca n   b s tu d ied   in   d i f f er en t i s lan d   m o d s ce n ar io s :   a.   Scen ar io - 1   ( r ef er en t) L ar g p ar o f   So u th - E ast  E u r o p ea n   P o w er   s y s te m s   i s o lated   f r o m   th r est  o f   E NT SO - E   b.   Scen ar io - 2 :   Fo u r   P o w er   S y s t e m s   ( Gr ee ce ,   A lb an ia,   Ma ce d o n ia  an d   Ko s o v o )   is o lated   f r o m   t h r est o f     E NT SO -   c.   Scen ar io - 3 : K o s o v o   P o w er   S y s te m   i s o lated     Fo r   ea ch   s ce n ar io ,   d if f er en i m b a lan ce   le v els  ar s i m u lated 5 0 %,  4 0 %,  3 0 an d   2 0 %.   I n   Fig u r 3   is   s h o w n   th f ir s t sce n ar io   w h e n   SEE   p o w er   s y s te m s   r e m ai n   is o lated   f r o m   th r e s t o f   E NT S O - E .   Sce n ar io - 2   an d   3   co n s id er s   s m a ller   p ar t o f   p o w er   s y s te m s   o p er atin g   i n   is la n d   m o d e.           Fig u r 3 .   So u th - E as t E u r o p ea n   P o w er   s y s te m s   i s o lated   f r o m   th r est o f   E NT SO - E   ( Sce n ar i o - 1)       4 0 %     3 0 % 2 0 % 5 0 %     F o r :   5 0 %   i m b a l a n c e ;     d f / d t   = 2 . 5 5   H z / s               4 0 %   i m b a l a n c e ;     d f / d t   = 2 . 2 5   H z / s             3 0 %   i m b a l a n c e ;     d f / d t   = 2 . 0   H z / s             2 0 %   i m b a l a n c e ;     d f / d t   = 1 . 5 5   H z / s   R O M A N I A B U L G A R I A S E R B I A K o s o v a M O N T E N E G R O B O S I N I A   &   H E R C E G O V I N A A L B A N I A M A C E D O N I A G R E E C E C R O A T I A S L O V E N I A X X X X I T A L Y T U R K Y X U K R A I N E U K R A I N E X X X X X X X X X X X X I T A L Y A U S T R I A H U N G A R I A X X M o l d a i a Evaluation Warning : The document was created with Spire.PDF for Python.
I n t J   E lec  &   C o m p   E n g     I SS N:  2088 - 8708       R ev iew   o f u n d er F r eq u en cy   L o a d   S h ed d i n g   P r o g r a o f Ko s o vo   P o w er … .   ( Ga z men d   K a b a s h i)   745   B y   m ea n s   o f   s y s te m   s i m u lat i o n s ,   th f o llo w i n g   r ec o m m en d atio n s   f o r   S y s te m   Op er ato r   o f   Ko s o v o   w il l b d eter m in ed   f o r :   a.   Op ti m al  to tal  s h ed d in g   lo ad   in   p er ce n tag o f   r ef er e n ce   lo ad ,     b.   Op ti m al  f r eq u en c y   s tep p in g   f o r   s y s te m   w it h   d is p er s ed   f r eq u en c y   r ela y s   i m p le m e n ted ,   c.   Op ti m al  n u m b er   o f   lo ad   s h ed d in g   s ta g es i n   p er ce n tag o f   r ef er en ce   lo ad .       3.   SIMUL A T IO N  M O D E L   T h p o w er   f lo w   an d   d y n a m ic  s i m u lat io n s   ar ca r r ied   o u o n   lar g e   s ca le   P SS /E  m o d el  h a v in g   o v er   1 4   in ter co n n ec ted   ar ea s   r ep r e s en ti n g   g i v e n   co u n t r ies  o f   So u th ea s E u r o p ( SEE )   p o w er   s y s te m s ,   in cl u d i n g   th Ko s o v o   P o w er   S y s te m   [ 5 ] .   T h lo ad   m o d els  w ith   v o lta g an d   f r eq u e n c y   d ep en d en ce   h as  b ee n   in cl u d ed   in   th d esi g n   f o r   UF L S i n   o r d er   t o   ac h iev ac cu r ate  ac ti v p o w er   i m b alan ce   est i m a tio n .   Usi n g   lar g s ca le  m o d el  t h s i m u lta n eo u s   o p en i n g   o f   i n ter c o n n ec tio n   lin e s   w as  s i m u lated   w it h   t h e   g en er ato r   u n i ts   u n d er   o p er atio n   a f f ec ted   b y   d i f f er e n lev el  o f   g e n er atio n   i m b ala n ce s ,   co n s id er in g   th r ee   s ce n ar io s   o f   i s o lated   s y s te m s   i n   w i n ter   r eg i m o f   o p er atio n .     T h f o cu s   o f   m o n ito r i n g   is   t h s y s te m   f r eq u e n c y   o f   Ko s o v o   P o w er   S y s te m   in   o r d er   to   d ef in e   t h o p tim a p r o g r a m   f o r   UFLS.   T h f r eq u en c y   r esp o n s o f   th s y s te m   i s   an a l y ze d   w it h   r e s p ec o f   E NT SO - r ec en tl y   u p d ates r eg ar d in g   UF L S,  w h ic h   is   d e s cr ib ed   in   P o lice  5 ,   o f   Op er atio n al  Han d b o o k .     3 . 1 .   L o a d sh ed din g   re la y   m o del i ng   T h m ax i m u m   n u m b er   o f   s tep s   i m p le m e n ted   i n to   t h m o d el  is   eq u al   to   8 .   Fo r   ea ch   s tep ,   th r es h o ld   an d   p er ce n tag o f   d is co n n ec ted   lo ad   is   ass o ciate d ,   j o in tl y   w it h   d ela y   t h at   r e p r esen th r ela y   co m p u tatio n al  ti m e   an d   t i m e   t o   ex ec u te   an d   o p en   t h lo ad   f e ed er   cir cu it b r ea k er .     D is co n n ec tio n   d ela y   ca u s ed   b y   cir cu it  b r ea k er   s w i tch i n g   ti m an d   r ela y   o p er atio n   ti m ( w it h o u ad d itio n al  ti m d ela y )   is   tak en   to   b 1 0 0   m s .   W ith   ai m   to   d ef in t h b est  p r o g r a m   o f   UF L S,  w h ic h   f its   t o   th Ko s o v o   P o w er   S y s te m   f r eq u en c y   ch ar ac ter is tic s ,   t w o     d i f f er e n UFLS  p r o g r a m s   ar an a l y ze d   an d   r esu lt s   ar co m p ar ed :   a.   6   s tep   UFL w i th   6 x 8 =4 8 % o f   to tal  lo ad   s h ed   s tar ti n g   f r o m   4 9 . 0 0   Hz  an d   en d in g   at  4 8 . 0 0   HZ   b.   8   s tep   UFL w i th   8 x 6 =4 8 % o f   to tal  lo ad   s h ed   s tar ti n g   f r o m   4 9 . 0 0   Hz  an d   en d in g   at  4 8 . 0 0   HZ   T h f r eq u en c y   r an g 4 9 . 0 0 - 4 8 . 0 0   Hz   w h er th U F L s h o u ld   ac t,  is   m an d ato r y   f o r   ea ch   T SO - o p er atin g   i n   s y n ch r o n o u s   zo n o f   E NT SO - E .   T h f ir s s tep   o f   lo ad   s h ed d i n g   is   f i x ed   at  4 9 . 0   Hz  T h r ea s o n   is   to   r eser v a   r an g ( 1   Hz)   b et wee n   5 0   Hz   an d   4 9   Hz   w h er p r i m ar y   r eser v e   is   tr y in g   to   r ec o v er   t h e f f ec f r o m   th p o w er   d e f icit.  T h is   r an g o f   1   Hz  s h o u ld   co n tr o th u n d er   f r eq u en c y   tr an s ien b y   lo a d s   s h ed d in g .   B elo w   th is   f r eq u en c y   th er e   is   ce r t ain   m ar g in   ( ar o u n d   0 . 5   Hz)   w h er g en er ati n g   u n its   ca n   o p er ate  an d   p r o b ab ly   r ec o v er   w it h o u tr ip .   Ne w   E NT SO - E   g r id   co d e   f o r   n e w   g en er atio n   u n its   t y p is   m o r r ig o r o u s l y   as   p r ev io u s   n atio n al  g r id   co d es  r eg ar d in g   ti m d u r atio n   o f   i s la n d   o p er atio n   m o d o f   g en er at o r s ,   f o r   ex a m p le  t h e   g en er ato r s   s h o u ld   o p er ate  at  le ast 3 0   m in   i n   u n d er   f r eq u e n c y   r eg i m 4 8 . 0   Hz - 4 7 . 5   Hz.   [ 6 ]       4.   SI M UL AT I O R E S UL T S   Fig u r 4   s h o w s   s i m u latio n   r esu lt  o f   s y s te m   f r eq u e n c y   r esp o n s ( 4 0 0   k b u s   b ar   at  SS   Ko s o v o   B )   af ter   i s o latio n   o f   SEE   P o w er   s y s te m   f r o m   t h r est   o f   E NT SO - E   n et w o r k ,   w it h   5 0   %   g e n e r atio n   i m b a lan ce .   I n   th is   ca s th s ix   s tep   UF L p r o g r a m   is   co n s id er ed .   A ll  s i x   s tep s   ar ac tiv ated   w it h   6 x 8 % =4 8 o f   to tal  lo ad   s h ed .   C o n s id er i n g   th e   ti m o f   1 0 0   m s   b et w ee n   u n d er - f r eq u e n c y   r ela y   tr ig g er in g   a n d   cir cu i t b r ea k er   o p er atin g ,   it  ca n   b s h o w n   t h at  f o r   ar o u n d   2 . 2   s ec o n d s   th e f f ec t s   o f   UF L S   o n   f r eq u e n c y   r ec o v er y   ar v i s ib le  an d   s u cc e s s f u l.  T h in ter v al s   b et w ee n   t w o   l o ad   s h ed d in g   b lo ck s   ar e   0 . 2   Hz.   T h m i n i m u m   f r eq u en c y   d ec li n is   4 7 . 9   Hz,   w it h   0 . 4   Hz  p o s itiv m ar g i n   f r o m   cr it ical  s y s te m   f r eq u en c y   o f   4 7 . 5   Hz  [ 7 ] .   Du r in g   t h f r eq u e n c y   d ec lin in   t h i s   ca s all  p r i m ar y   p o w er   r eser v e s   f r o m   s y n c h r o n o u s   g en er ato r s   lo ca ted   in   e f f ec ted   ar ea   w il b au to m at icall y   ac t iv ated   tr y i n g   to   r ec o v er   th s y s te m   f r eq u e n c y .   Us u all y   i n   s u c h   cir cu m s ta n ce s   e f f ec ts   o f   p r im ar y   p o w er   r eg u latio n   i s   n o c h ar ac ter ized   as   f ast   r esp o n s e   b ased   o n   w ell - k n o w n   s p ee d   g o v er n o r   ch ar ac ter is tic s   o f   s tea m   t u r b i n g en er ato r s .   I ca n   b co n s i d er ed   th at  v i s i b le  e f f ec ts   ca n   b o b s er v ed   af ter   1   s ec o n d   f r o m   t h i n itial  tr i g g e r   o f   s p ee d   g o v er n o r .     I n   th i s   v er y   s h o r ti m f r a m an d   a f t er   th to tal  s y s te m   in er tia  co n s tan t p la y s   s i g n i f i ca n t r o le  o n   f r eq u e n c y   b e h av i o r   af ter   d is tu r b an ce   o cc u r s .   Fig u r 5   s h o w s   s i m u la tio n   r esu lt   o f   s y s te m   f r eq u en c y   r e s p o n s w i th   s a m s y s te m   co n d itio n   a s   p r ev io u s   o n e,   b u w it h   d if f er en U FLS  p r o g r a m   w h ich   co n s id er   ei g h s tep s   w it h   8 x 6 % =4 8 o f   to tal  lo ad   s h ed .   T h in ter v als  b et w ee n   t w o   lo ad   s h ed d in g   b lo ck s   ar 0 . 1 4 2 8   Hz,   m o r clo s to   ea c h   o th er   t h a n   s i x   s tep   UFLS  p r o g r a m .   T h m i n i m u m   f r eq u en c y   d ec lin i s   4 7 . 8 5   Hz ,   w it h   0 . 0 5 Hz  h ig h er   th a n   s i x   s tep   UFL S.   Fig u r 6   s h o w s   f r eq u e n c y   r es p o n s d if f er e n ce s   b et w ee n   t wo   UFL p r o g r a m s   w it h   s a m s i m u lat io n   co n d itio n   co n s id er in g   5 0   % g e n er atio n   i m b ala n ce .     Evaluation Warning : The document was created with Spire.PDF for Python.
                        I SS N :   2 0 8 8 - 8708   I n t J   E lec  &   C o m p   E n g ,   Vo l.  8 ,   No .   2 A p r il   201 8   :   7 4 1     7 4 8   746   As  is   s h o w n   in   Fi g u r 6 ,   d i f f e r en ce s   b et w ee n   t w o   U F L p r o g r a m s   ar r elati v el y   s m all,   e v en   w ca n   b o b s er v ed   th at  8   s tep s   UFLS   s h o w s   s o m e h o w   b etter   r es u lts   b u t n o t so   s i g n if ica n tl y .   Si m u latio n   w it h   th s a m i m b alan ce   co n d itio n   is   p er f o r m ed   also   f o r   s m al ler   s ize  o f   is o lat ed   p o w er   s y s te m s ,   s ce n ar io - 2   ( f o u r   r e m ai n in g   p o w er   s y s te m s )   an d   s ce n ar io - 3   w it h   o n l y   Ko s o v o   P o w er   S y s te m   i n   o p er atio n .   Fig u r 7   s h o w s   s i m u latio n   r e s u l ts   o f   f r eq u e n c y   r esp o n s f o r   t h r ee   s ce n ar io s   o f   d if f er en s ize  o f   p o w er   s y s te m s   r e m a in i n g   in   o p er atio n   in   is la n d   m o d e.           Fig u r 4 . S y s te m   f r eq u e n c y   r es p o n s w ith   5 0 g en er atio n   i m b alan ce   an d   s i x   s tep   UF L S a ct iv atio n     ( Scen ar io - 1)           Fig u r 5 .   S y s te m   f r eq u e n c y   r e s p o n s w it h   5 0 % g e n er atio n   i m b a lan ce   a n d   eig h t step   UF L S a ctiv atio n   ( Scen a r io - 1)       L ar g er   t h p o w er   s y s te m   i s ,   th lo w er   f r eq u e n c y   e x c u r s i o n   w ill  o cc u r s .   T h f r eq u e n c y   r ate  o f   ch an g f o r   d if f er en s ize  o f   is o lated   p o w er   s y s te m s   an d   d if f er en lev el  o f   g en e r atio n   i m b a lan ce s   ar s h o w n   in   T ab le  1 .   A s   is   s h o w n   i n   th tab le  th f r eq u e n c y   r ate  o f   c h an g is   lo w er   w h en   t h s ize   o f   is o lated   p o w er     S E C I   R E G I O N A L   M O D E L   -   W I N T E R   M A X I M U M   2 0 2 0 3 R D   W E D N E S D A Y   I N   J A N U A R Y   2 0 2 0   A T   1 9 : 0 0   C E T T i m e   ( s e c o n d s ) 6 5 . 6 5 . 2 4 . 8 4 . 4 4 3 . 6 3 . 2 2 . 8 2 . 4 2 1 . 6 1 . 2 0 . 8 0 . 4 0 1 2   -   F R E Q   4 7 1 0 0 5   [ _ K O B S P 1           4 0 0 . 0 0 ]   :   S E E 4 0 % 0 - 0 . 0 0 5 - 0 . 0 1 - 0 . 0 1 5 - 0 . 0 2 - 0 . 0 2 5 - 0 . 0 3 - 0 . 0 3 5 - 0 . 0 4 - 0 . 0 4 5 0 4 9   H z 4 8   H z I f   f < 4 7 . 5   H z     g e n e r a t o r   t r i p   a f t e r   2 0   s 4 8 . 8   H z   s t e p   2 4 8 . 6   H z   s t e p   3 4 8 . 4   H z   s t e p   4 4 8 . 2   H z   s t e p   5 4 9   H z   - s t e p   1 4 8   H z   s t e p   6 5 0   H z R e l a y   a c t i v a t i o n C i r c u i t   b r e a k e r   o p e n F r e q u e n c y   r a n g e   o f   6   s t e p     l o a d   s h e d   a c t i v a t i o n 5 0 %   g e n e r a t i o n   i m b a l a n c e     S i x   s t e p   a c t i v a t i o n   o f   U F L C   f o r   5 0 %   g e n e r a t i o n   i m b a l a n c e F r e q u e n c y   d e v i a t i o n   ( p u . ) 6   s t e p s   e a c h   8 %   o f   l o a d   s h e d f m i n = 4 7 . 9   H z   2 0 % 3 0 % 4 0 % 5 0 % f m i n = 4 7 . 8 5   H z U F L S   R e l a y   a c t i v a t i o n 4 9 . 0   H z 4 8 . 8 5   H z 4 8 . 7 1   H z 4 8 . 5 7   H z 4 8 . 4 3   H z 4 8 . 2 8   H z 4 8 . 0 5   H z 4 8 . 3 5   H z 4 8 . 6 5   H z 4 8 . 1 4   H z 4 8   H z Evaluation Warning : The document was created with Spire.PDF for Python.
I n t J   E lec  &   C o m p   E n g     I SS N:  2088 - 8708       R ev iew   o f u n d er F r eq u en cy   L o a d   S h ed d i n g   P r o g r a o f Ko s o vo   P o w er … .   ( Ga z men d   K a b a s h i)   747   s y s te m   i s   lar g er .   T h er ar t w o   f ac to r s   i n f lu e n ci n g   th v a lu o f   f r eq u e n c y   r ate  o f   ch a n g e:  th a v ailab le  o f   p r im ar y   r eser v a n d   th i n er tia   co n s tan t o f   t h s y s te m .   I t c an   b ex p ec ted   th at  th lar g s y s t e m   w ill  h a v m o r e   av ailab le  p r i m ar y   p o w er   r eser v es  a n d   also   lar g er   to tal  eq u i v alen i n er tia  co n s ta n b ec au s h i g h   n u m b er   o f   an d   h i g h   d i v er s it y   o f   s y n ch r o n o u s   g e n er atio n   t y p es  w it h   d if f er en t p ar a m eter s   [ 8 ] .           Fig u r 6 .   Fre q u en c y   r e s p o n s co m p ar is o n   b et w ee n   s ix   a n d   eig h t step   UF L f o r   5 0 % g en er atio n   i m b ala n ce   ( s ce n ar io - 1)           Fig u r e   7 .   Fre q u en c y   r e s p o n s f o r   d if f er e n t size  o f   is o lated   p o w er   s y s te m s   f o r   s ix   s tep   UF L S p r o g r a m .         T ab le  1 .   T h f r eq u en c y   r ate  o f   ch an g f o r   d if f er en t size  o f   i s o lated   p o w er   s y s te m s   a n d   d if f er en t i m b ala n ce s       F r e q u e n c y   r a t e   o f   c h a n g e   I mb a l a n c e   2 0 %   3 0 %   4 0 %   5 0 %   I s l a n d   mo d e   d f / d t     [ H z / s]   d f / d t     [ H z / s]   d f / d t     [ H z / s]   d f / d t     [ H z / s]   S EE   1 . 5 5   2 . 0   2 . 2 5   2 . 5 5   G r ,   A l ,   M k ,   K s   1 . 2   1 . 6 5   2 . 1 5   2 . 4 5   K o s o v o   ( K s )   1 . 1 5   1 . 8   2 . 5 5   2 . 9           C o m p a r i s o n   b e t w e e n   6   a n d   8   s t e p   U F L S   f o r   5 0 %   i m b a l a n c e 8   s t e p   U F L S 6   s t e p   U F L S     C o m p a r i s o n   o f   f r e q u e n c y   r e s p o n s e   f o r   d i f f e r e n t   s i z e   o f   i s o l a t e d   p o w e r   s y s t e m s   K S - i s o l a t e d S E E - G R , A L , M K , K S Evaluation Warning : The document was created with Spire.PDF for Python.
                        I SS N :   2 0 8 8 - 8708   I n t J   E lec  &   C o m p   E n g ,   Vo l.  8 ,   No .   2 A p r il   201 8   :   7 4 1     7 4 8   748   5.   CO NCLU SI O N   C o m p ar is o n   o f   s i m u latio n   r es u lts   b et w ee n   6   s tep   U FLS  a n d   8   s tep   UFL S   w it h   4 8 o f   to tal  lo ad   s h ed ,   s h o w s   n o   s i g n if ican t   d i f f er en ce s   w h en   t h e y   ap p l y   i n   So u t h ea s P o w er   S y s te m   r e g io n   o f   E NT SO - E.  E v en   w it h   t h ac tiv atio n   o f   all  s tep s ,   th en d   f r eq u e n c y   i s   n ev er   f u ll y   r ec o v er ed .   T h 8   s tep   UFL s h o ws   s o m b etter   p er f o r m a n ce   r eg a r d in g   m i n i m u m   f r eq u e n c y   d ec a y   an d   e n d   f r eq u e n c y .   F u r th er m o r 8   s tep   UF L S   h as  s o m e   s i g n if ican t   d is ad v a n t ag es r elate d   to   th e     n ar r o w   ti m in ter v als     b et w ee n   t w o   lo ad   s h ed d in g   b lo ck s   i n   co m p ar is o n   w it h   6   s tep   U F L S.  W h e n   w co n s id er   t h at  lar g n u m b er   o f   cir cu i b r ea k er s   i n   d is tr ib u tio n   ne t w o r k   h a s   d if f er en b r ea k er s   o p er atio n   tim d u to   th t ec h n o lo g y   an d   th a g e,   th er e   is   p o s s ib ilit y   o f   s i m u lta n eo u s   tr ig g er i n g   o f   t w o   b lo ck s   w h e n   p o w er   s y s t e m   ex p er ien ce s   f r eq u en c y   d ec a y .   W h ile  r ela y   m an u f ac t u r er s   p r o v id s et  p o in ts   i n   in cr e m en ts   o f   0 . 0 1 Hz ,   th er is   ty p ical l y   1 0 0 - 1 5 0   m s   o f   ti m d ela y ,   in cl u d in g   r ela y   an d   cir cu it  b r ea k er   o p er atin g   ti m e s ,   f r o m   th r ela y   tr ig g er in g   w h e n   t h f r eq u e n c y   r ea ch   s etti n g   p o in to   th m o m e n o f   lo ad   s h ed .   Du r in g   th i s   d ela y ,   th f r eq u en c y   co n ti n u e s   t o   d r o p ,   th at  is ,   th f r eq u en c y   at  w h ic h   lo ad   ac t u a ll y   s h ed s   i s   b elo w   t h s et  p o i n t.  I f   t h s p ac in g   b et w ee n   t h e   s h ed d in g   s ta g es   is   to o   clo s e,   th lo ad   s h ed   in itiated   b y   s tag co u ld   ac tu all y   s h ed   th lo ad   w h ile  t h f r eq u en c y   i s   in   p h a s o f   r ec o v er y   cr ea ti n g   an   u n ac ce p t ab le   en d   f r eq u en c y   w h ic h   ca n   b to o   h ig h .   T h er ef o r th r ec o m m e n d atio n   f o r   th S y s te m   Op er ato r   o f   Ko s o v o   is   to   ad o p t 6   s tep s ,   6 x 8 =4 8 % lo ad   s h ed   U F L S p r o g r a m ,   w h ic h   f u ll y   co m p lie s   w it h   E NT SO - E   r eq u ir e m en ts .   T h 6   s tep   UFL S p r o g r a m   is   p r esen ted   in   T ab le   2   b elo w :       T ab le  2 .   Six   s tep ,   6 x 8 % UF L S   s ettin g s   U F L S   S t e p   1   S t e p   2   S t e p   3   S t e p   4   S t e p   5   S t e p   6   L o a d   sh e d     8%   8%   8%   8%   8%   8%   I n i t i a l   f r e q u e n c y   4 9 . 0   H z   4 8 . 8   H z   4 8 . 6   H z   4 8 . 4   H z   4 8 . 2   H z   4 8 . 0   H z       RE F E R E NC E S     [1 ]   h tt p s:/ /www . e n tso e . e u /f il e a d m in /u se r_ u p lo a d /_ li b ra ry /p u b li c a ti o n s /en tso e /Op e ra ti o n _ Ha n d b o o k / P o l ic y _ 3 _ f in a l. p d f   [2 ]   G rid   Co d     KO S T T ,   ww w . k o stt. c o m   [3 ]   L o a d   S h e d d in g ,   L o a d   Re st o ra ti o n   a n d   G e n e ra to P ro tec ti o n   Us in g   S o li d - sta te an d   El e c tro m e c h a n ica ww w . G Ei n d u strial. c o m /p m   [4 ]   Ku n d u r,   P .   P o w e S y ste m   S tab il it y   a n d   Co n tro l.   Ne w   Yo rk M c G r a w - Hill   ,   1 9 9 4 .   [5 ]   T ra n s m issio n s Ne t w o rk   De v e lo p m e n P lan   ( 2 0 1 7 - 2 0 2 6 w ww . k o st t. c o m   [6 ]   h tt p s:/ /www . e n tso e . e u /m a jo r - p ro j e c ts/n e tw o rk - c o d e - d e v e lo p m e n t/ re q u irem e n ts - f o r - g e n e ra to rs/P a g e s/d e f a u lt . a sp x   [7 ]   P r o tec ti v e   Re la y in g - P rin c i p le an d   A p p li c a ti o n ,   th ir d   E d it i o n ,   J.L e wis Blac k b u rn ,   T h o m a s J.Do m in   CRC P re ss ,   Ne w   Yo rk ,   IS BN 1 - 5 7 4 4 4 - 7 1 6 - 5 ,   p p . 5 4 1 - 638   [8 ]   h tt p s:/ /www . e n tso e . e u /Do c u m e n ts/P u b l ica ti o n s/ S OC/No rd ic/No rd i c _ re p o rt _ F u tu re _ S y ste m _ In e rti a . p d f # s e a rc h = F u t u re % 2 0 sy ste m % 2 0 in e rti a       B I O G RAP H I E S   O F   AUTH O RS        G a z m e n d   K a b a s h i   re c e iv e d   h is  e lec tri c a e n g in e e rin g   De g re e   f ro m   Un iv e rsit y   o P rish ti n a ,   Ko so v o   in   2 0 0 3 .   T h e   P h . D.  d e g re e   re c e iv e d   in   2 0 1 3   in   t h e   sa m e   F a c u lt y .   S in c e   2 0 0 3   h e   is  w it h   T ra n s m issio n   S y st e m   a n d   M a rk e Op e ra to o f   Ko so v o   (KO S T T w o rk in g   o n   lo n g   term   d e v e lo p m e n o tran sm issio n   n e tw o rk .   He   h a a u th o re d   o v e 2 0   p a p e rs  in   p e e r - re v ie w e d   jo u rn a ls.   His  re se a rc h   in tere sts a re   in   tran sie n sta b il it y   o f   p o w e s y ste m s,  a n d   W in d   P o we e n e rg y .             S k e n d e r   K a b a s h i   is   A ss o c iate   P ro f e ss o o f   P h y sic s   a n d   h e a d   a De p a rtm e n o P h y sic s   o f   F a c u lt y   o Na tu ra S c ien c e s,  Un iv e rsity   o P r ish ti n a .   Ho l d P h in   En v ir o n m e n tal  P h y sic a n d   Re n e wa b le  En e rg y   o n   Jo z e f   S te f a n   In stit u te  in   L ju b lj a n a ,   S lo v e n ia  in   2 0 0 8 .   P ro f .   Ka b a sh is  a c ti v e   in   b o t h   n u c lea p h y sic a n d   in   re se a rc h   o f   re n e w a b le  e n e rg y   so u rc e s.   He   h a a u th o re d   o v e 3 0   p a p e rs  i n   p e e r - re v ie w e d   jo u rn a ls.   Its  sp e c if ic   sc ien ti f ic   re se a rc h   in tere st   is  th e   D y n a m ic   M o d e li n g   th e   En v iro n m e n tal  S y ste m s p a rti c u larly   m o d e li n g   th e   En e rg y   a n d   En v ir o n m e n t,   a s we ll   a s th e   s y ste m a ti c   m e a su re m e n ts  re late d   to   t h e   e n v ir o n m e n tal  p ro tec ti o n   a n d   im p a c o f   re n e w a b le en e rg y   o n   su sta in a b le  e n e rg y   p ro d u c ti o n   a n d ,   M o d e li n g   o En e rg y   S y st e m s   re late d   to   Cli m a te  Ch a n g e   f o S u sta in a b le  De v e lo p m e n t.       Evaluation Warning : The document was created with Spire.PDF for Python.