Indonesi an  Journa of El ect ri cal Engineer ing  an d  Comp ut er  Scie nce   Vo l.   23 ,  No.   1 ,   Ju ly   20 21 pp.   23 ~ 31   IS S N: 25 02 - 4752, DO I: 10 .11 591/ijeecs .v 23 .i 1 . pp 23 - 31          23       Journ al h om e page http: // ij eecs.i aesc or e.c om   Mont e Car l o anal ysis for sola r PV i mp ac t as sess m ent  in MV  distribu tion netw or ks       Dili ni Almei d a 1 Jagadees Pasupule ti 2 S ha n gari  K. R aveendr an 3 M. Re ya s udin B as ir  K han 4   1, 2, 3 Coll ege  of   E ngine er ing ,   Univ ersit i   Te n aga Na sional   (UN ITE N ),   Putra jay a ,   Mal a y si   4 Manipa l   Interna ti onal Unive rsi t y Nil ai,  Neg eri Sem bil an ,   Mal a y s ia       Art ic le  In f o     ABSTR A CT   Art ic le  history:   Re cei ved   Ma r   9 , 2 021   Re vised  Jun   7 ,   2021   Accepte J un  1 8 , 202 1       The   rap id  pen etrati on  of  solar  photovol taic  (PV sy stems   in  distri buti o n   net works   has  impos ed  var ious  implic at ions   on  net work  oper at ions .   The ref or e,   i is  i m per at ive   to  con sider  the   stoch asti na ture   of  PV   gene ra ti on   and  lo ad  deman to  add ress  th oper at ion al   c hal l enge in   fut ure   PV - ric h   distri buti on  n etw orks.  Thi pap er  proposes  M onte   Car lo  base d   proba bilis ti c   fra m ework  for  a ss essing  the   impact  of  PV   pen et ra ti on  on   m edi um   volt a g e   (MV distri buti o net works .   The   unce rt ai nt ie s a ss oci a te wi th  PV   insta llati o n   c apaci t y   and  it l oca t ion,   as  well   as  the   ti m e - var y i ng  nat ure   of  PV   gene ra ti on   and  loa demand  were   conside red   for  the   implementa ti on  of  the  proba bil ist i c   fra m ework.   ca se  stud y   w as   per form ed  for   t y pical   MV   distri but ion   net work  in  Malay s ia,  demons tra ti ng   th eff e ct iv ene ss   of  M onte   Car l o   ana l y sis  in  evalu at ing  the   pot enti al   PV   impact in  the   future .   total  of  1000   Monte  Carl si m ula ti ons  were   conduc t ed  to  acc ura tel y   ide nt if y   t he  infl uen ce  of  PV   pene tra t i on  on  volt ag p rofil es  and  n et w ork  losses.  Besi des,   seve r al   ke y   m et r ic were   used  to  q uant if y   th tec hnic a per form a nce   of  the  distri buti on  n etw ork.   The   r esults   rev ea l ed  th at   t he  wors rep er cu ss ion  of  high  solar  PV   penetr at ion   on  t y p ic a l   Malay si an  MV   distri but ion  n etw orks  is  the  viol ation  of  the   upper   volt a g st at utor y   l imit,   wh ic is  li k ely   to  o cc ur  be y ond   70%  penetration   le v el .   Ke yw or d s :   Mon te  Ca rlo  a naly sis   MV d ist ri buti on  netw ork    PV   pe netrati on   Syst e m  p ower  loss   Vo lt age  r ise   This   is an  open   acc ess arti cl e   un der  the  CC  B Y - SA   l ic ense .     Corres pond in Aut h or :   Dili ni A lm ei da   Coll ege  of   En gi neer in g   Un i ver sit i Te na ga Nasi onal  (UN ITE N)   Jal an  I KR AM - UNITE N,   43 000 Ka j a ng, S el ango r,  Mal ay sia   Em a il dili nialme ida2@gm ail.co m       1.   INTROD U CTION   Ov e t he  la st  de cade,   m any  countries   ha ve  i niti at ed  the   pr ocess   of  set ti ng  up  acce ss   to   transm issi on   and   distrib utio netw orks  to  li ber al iz their  el ect rical   sys t e m [1 ] This  has  bee ass oc ia te with  subs ta ntial   grow t in   the   par ti ci patio of  distri bu te gen e rato rs  of   va rio us   te c hnologies.   O wing  to   the  pro po s e dev el op m ents  in  the  reducti on  of  ca rbo f oo t pr i nts  an rising  te m per at ur es,   m ajo r   sh ar of  distr ibu te gen e rato rs  is  e xp ect e to   be   fo c us e on  re ne wab le   e ne rg y   [2 ] .   Re centl y,   the  la r ge - scal integ rati on  of  s olar   PV   syst e m ha sho wn  re m ark able  e xpa ns io due  t t he  e nhance m od ule  pe rform ance,  ec onom ie of  scal e,  an rese arch  an de velop m ent  relat iv to  the   ot her   ren e wa ble  ene rg s ources D ur i ng   t he  per i od  of   2010  to  2019,  the  com po und  annual  gro w th   rate  of   so la PV   has  c ro s se over  40%  a nd  it glo bal  in sta ll ed   capaci ty   is  pr oject ed  to  r each   8519   G W   by  2050  [ 3].  In   Ma la ysi a,  the  gove rn m ent  has  intende to  a dd   a annual  s olar   ca pacit of  200  M W   in   Pe nin s ula  Ma la ysi a nd  50 M i S abah  th rou gh  t he  im p lem enta ti on   of   the lar ge - scal so la r (LS S) p r ogram  [ 4].    Evaluation Warning : The document was created with Spire.PDF for Python.
                          IS S N :   2502 - 4752   Ind on esi a J  E le c Eng &  Co m Sci,   Vo l.   23 , N o.   1 Ju ly   20 21 23   -   3 1   24   Ty pical ly tradit ion al   powe syst em are  co nf i gure as  passive   netw orks  t hat   co ns ide a   un i directi onal  p owe flo fee d: f r om  central iz ed  gen e rati on  to  con s um ers.  Thus,  the r a pid pr olife rati on of P s yst e m po ses  new   te ch nical   chall enges  on  the  operati on  and   perform ance  of   these  netw orks,   since  the hav e   no bee desi gned   to   facil it at dis per se ge ner at io n   [ 5] I par ti cula r,   MV  distrib utio netw orks   to  wh ic la rg e - scal s olar  plants  a re  i nteg rated  dir e ct ly   or   by  a ggreg at io of  lo vo lt a ge  in sta ll at ion are   a lread y   exp e rienci ng  te chn ic al   c halle ng e in  a reas  wh e re  high  P c on ce ntrati ons  e xist  [ 6 ] A m on these  te chn ic a l   chall enges,  the   vo lt age  rise h a been   ide ntifie as  the  m os c riti cal   chall enge  as  it s   i m p li ca ti on di rectl a ff ect   the p e rfo rm ance of th e  d ist ri buti on n et wor k [ 7 ] [ 8 ].   Ow i ng   to  the  s tochasti natu r of   the  outp ut,   PV   un it are  c on si der e unc ontr ollable As  resu lt the   var yi ng  ge ner a ti on   an load  dem and   co uld   c ause  over volt age  pro blem i MV  distribu ti on   netw orks.   Durin per i od w hen  low  dem and   a nd  hi gh  ge ner at ion   c oin ci des,     excess  ge ner a ti on   fl ow s   to w ard t he  s ubsta ti on  form ing   rev e rs powe flo ws   [ 9 ] C on se que ntly the  ordi na ry  be hav i our  of   t he  distri bu t ion   netw ork  va ries,  le ading   t i ncrea sed  volt age  l evels.  Be sides ,   vo lt age  vio la ti on c ould  oc cu with  f ur t her  el evati on   in  th PV  pen et rati on  le vel  in  the   net work.  Noneth el ess,  the  se ve rity   of   the se  i m pacts  dep en ds   on   seve ral  factors,   includi ng  netw ork  t opology,  geog raphic  loc at ion   a nd   siz of  instal le P syst em s,  an weathe c onditi ons   [ 10 ] T he  i ncr e ased  act ive   po wer  loss   of  t he   netw ork acc om pan ie by  hi gh  s olar  P i nt egr at io n,   is   an oth e chall enge  that  cou l affe ct   the  perf or m ance  of   the  distr ibu ti on  net wor k   [ 11 ] .   Us ually the  netw ork  los s   represe nts  U - s ha pe  tra j ec tory  a gainst  t he  P pe netra ti on   le vel  [ 1 2 ] It  is  e xpect ed  to   dec rease   unde r   m od erate  PV   pen et rati on,  as   the  inte gr at io of   PV   syst em red uces  th net  po wer  fl ow  withi the  netw or [1 3 ] Howe ver,   at   hi gh  P pe netrati on  l ev el s,  the   inc rease rev e rse  pow er  fl ow  c ou l r esult  in  hi gh  ne twork   losses  [1 4 ]   Re centl y,  c on si der a ble  num ber   of  rese arch  ef forts  ha ve  bee unde rtake to   in ve sti gate  the   po te ntial   te chn ic al   chall eng es   associat ed  with  high  pen et r at ion of  PV   in  distrib utio netw orks.     Mo st  of  these  PV   im pact  assessm ents   hav bee f oc us e on  assess ing   the  e ff ect   of   gr i d - c onnec te so la PV   on  the   syst e m   and   qu antify ing   t he  a ntici pated  im pacts  as  functi on   of   t he  pe net rati on   le vel,  co ns ide rin the  s ta tutes   and   sta ndar ds   The  m et ho ds   util iz ed  in  prev iou st ud ie t analy ze,  an quantify   the  so l ar  P im pacts  cou l be  m ai nly  ca te gorized  as  determ inist ic   and  pro ba bili sti ty pes.   In  the  dete rm i nisti appr oac h,   t he   un ce rtai nties  li nk e t P ge ner at io a re  id entifi ed  a nd  e va lua te us in s cenari o - base analy sis.  T he  worst - case  scena rio   analy sis  is  on e   of  t he  frequ e ntly   ad op te decisi on - m aking   te c hniqu e in   the  e ven of   un ce rtai nties.  Be sides,  this  m et ho does  not  co ntem plate  the  pr ob a bili sti aspect  of  the  syst em as  al the   inpu data  a nd  the  e ve nts  are  kn own  in   adv a nce O ne   of   t he  ke lim i ta ti on of  this  m et ho is  the  ov e resti m at ed  resu lt wh ic cou l dri ve  unnece ssary  gr i rei nfor cem ents.  Howe ver ,   in  the  prob a bili sti appr oach,  the  r andom   pr oba bili ty   of   an  eve nt   is  con side red  t assist   the  syst e m   un certai nties.  Mo nte  Ca rlo  i s   on of   t he  wi de ly   e m plo ye pro bab il ist ic   te chn i qu e that  gen e rates  re pe at ed  ra ndom   s a m pling   to  a ddress  the   un ce rtai nties i n t he  syst em .   Ma ny  determ i nisti ty pe  stu dies  hav e   bee re porte in   the  li te ratu re  t assess  t he  P V   pe netrati on   i m pacts  on   M distri bu ti on  netw orks.   I [ 1 5 ] determ i nisti ty pe  i m pact  assessm ent   has  been   pe rfor m ed  to   exam ine  the  ne twork  l os ses  and   volt age  iss ues  relat e to  PV   pen et rati on  in  MV  distrib ution  net works   unde r   diff e re nt  so la r   var ia bili ty   con diti on s sim i la stud has  be en  co nducted  in  [1 6 ] - [ 1 8 ]   to  evaluate  the  im pact   of   P pe netra ti on   on  distrib ution   netw ork s'   vo lt age  profi le   and   power  losses.  In   [1 9 ] the  auth ors  have   stud ie the  e f fect  of  P se lf - co nsum ption  on  the  volt age   le vels  a nd   syst e m   losses  of  S pa nish   MV   distrib ution   network,  co ns ide r ing   six  ge ogra ph ic   distri bu ti ons;  rand om di sp ers ed low,   m od erate,  strong,  an extrem e.  In   [ 20 ] the  im pact   of   P ge nerat ion   on  the  volt age  pro file s   and   powe lo sses  in  distri buti on   n et w orks  has  been  assesse d,  co ns ide rin di ff ere nt  P in ver te r   capaci ti es,  instal le at   dif fer e nt  loca ti ons .   Re fer e nces  [ 21 ] - [ 2 5 ]   hav pro posed  M on t Ca rlo  base pro bab il ist ic   i m pact  assess m ents  to  dete rm ine  the  influ e nce  of  P V/D pe netrat ion   on  MV d ist rib ut ion  n et w orks I m os of   these  publica ti on s the  aut hor hav e   ta ken   into  acc ount  the  num ber   of   P syst em s,   the  instal led   capaci ty and  the  locat ion   w it hin   the  distribu ti on   netw ork  as  un certai nties  ass ociat ed  with  s olar  P gen e r at ion In   [ 2 6 ] stud has  b e en  pe r form ed  to  assess  the  i m pact  of   PV   syst em s   on  ru ral  MV  di stribu ti on  net w ork  by  util iz ing   Mon te   Ca r lo  based   te c hn i qu to  represe nt the i nterm ittency  o f  PV pro duct io n an loa ds.   Wh il a   fe st ud ie s   ha ve  pro po s ed   Mo nte  Ca rlo  a naly sis  f or   so la P i m pact  assess m ent  in     M V   distrib ution   ne tworks the  stu dies  that  consi der   detai le an reali sti network   m od el with  high  res ol ution,   tim e - series  representat io ns   ar lim it ed.   The refor e this  pa per   pro po se and   a dopts  Mon te   Ca rlo  base pro bab il ist ic   f ram ewo rk  to  inv est igate   t he   pote ntial   te chn ic al   im pacts  of  s olar  P V   pe netrati on  on  MV   distrib ution n et works  wh il e ca te ring to  t he unc ertai nties rela te to  s olar  P V gen e rati on.    The  rem ai nd er  of  the   pa pe is  orga nized   as  fo ll ows:  Sect ion  outl ines  the   re sea rch  m et hod  includi ng  the   pro po se M on te   Ca rlo - base pro ba bili sti fr am ewo rk t he   i m pact  m e tric s a nd  the   case   stu dy .   T he  sim ulati on   re su lt ob ta ined  from   the  ap plica ti on   of  the  pro pose m et ho dolo gy   are  prese nted  a nd   discusse d i n Se ct ion 3 Finall y, the c on cl us i on s  are  dra wn in S ect io 4 .   Evaluation Warning : The document was created with Spire.PDF for Python.
Ind on esi a J  E le c Eng &  Co m Sci     IS S N:  25 02 - 4752       Mo nte Carl o a na ly sis f or  s olar PV i mpact  asse ssme nt in  M V d ist rib utio n netw or ks   ( Dili ni Almei da )   25   2.   RESEA R CH MET HO D     This  sect io presents  t he  pro po s ed  im pact  assessm ent  fr am ewo r a nd   t he  im pact  m e t rics  us e t analy ze an d q ua ntify t he  im pact o f  so la P V pen et rati on on M distri buti on  netw orks.     2 . 1.      O ver vie w  o f the   i mp ac a ssessme nt   f ramew ork   The  pro posed   Mo nte  Ca rlo   base pro ba bili sti i m pact   assessm ent  f ram ewo rk  is  capab le   of   inv est igati ng  t he  po te ntial   te chn ic al   im pacts  of  s olar  PV  pe netrati on  on  MV   distri buti on   netw orks  wh il e   addressi ng   t he   uncertai nties  and  inter m it ten nat ur of  s olar  P gener at ion T he  flo wch a rt  of  the   i m pac t   assessm ent  fr am ewo r is  sho wn   in  Fig ur 1.  For  giv e MV  distrib utio net work,  this   fr am ewo r ev al uates  the possi ble im pacts  of   so la r  P V pen et rati on in  te rm s o f  sev e ral m et rics.       S t a r t M od e l l i n g M V  ne t w or k m ode l T opo l ogi c a l  a nd  e l e c t r i c a l  da t a  of  t he   ne t w or k D a i l y l oa d a nd P V   ge ne r a t i on p r of i l e s I m p a c t   a s s e s s m e n t R un t hr e e - pha s e  t i m e - s e r i e s   pow e r  f l ow  s i m ul a t i ons A not he r  P V   de pl oym e nt  s c e na r i o t a na l ys e ? R e c or d t he  s i m ul a t i on r e s ul t s E nd A na l ys e  t he  r e s ul t s Y e s No M on t e   C ar l I m p l e m e n t at i on R a ndo m l y gene r a t e  P V  c a pa c i t y a nd l oc a t i on  I nt e gr a t e  P V  s ys t e m s  t o t he  ne t w or k S e l e c t  t he  P V  pe ne t r a t i on l e ve l A not he r  P V   P e ne t r a t i on l e ve l   t a na l ys e ? No Y e s     Figur 1.  Flo w char of the  pro po s ed  Mo nte C arlo  s olar  P i m pact assess m ent fram ewo r k   Evaluation Warning : The document was created with Spire.PDF for Python.
                          IS S N :   2502 - 4752   Ind on esi a J  E le c Eng &  Co m Sci,   Vo l.   23 , N o.   1 Ju ly   20 21 23   -   3 1   26   Firstl y, the to polo gical  and el ect rical  d at a of  the select ed  M distri buti on   netw ork  is  us e to  de velo detai le netw ork  m od el   c onnected  with  so l ar  P syst em s.  In  or der   t c o nduct  reali sti analy sis,  different   load  pr of il es  and   ir rad ia nce  prof il es  are  e m plo ye to  incorp or at real - tim e   var ia ti on s   in  PV   ge ner at ion   an load  dem and As  the  first  ste to wards  the  Mote   Ca rlo  im plem entat ion la rg num ber   of  PV   de ploym ent  scenari o are  consi der e by   rando m ly   a llo cat ing   differe nt  PV   capaci ti es  acro ss  the  netw ork.   T herea fter,   three - phase,  ti m e - series,  m inu te - by - m inu te daily   po we f low  sim ulati on are  perform e to  inv e sti gate  the   po s sible  s olar  PV   im pacts  under   dif fer e nt  pe netrati on  le ve ls.  The  m ai ste ps   of  si ng l si m ulati on   unde a   par ti cula r pene trat ion  le vel ar e su m m arized b el ow.     Allocat ion  of   l oad  co nsum pti on  pr of il es  f or   reside ntial c om m ercial and   industrial   loa ds  co nnect ed  t the  netw ork.     Ra ndom   al locat ion   of  PV  ca pacit ie acr os s   the  netw ork  unti the  de fine le vel  of  penet rati on   is  reac hed  (r a ndom ly   gen erated  P ca pa ci ti es  are  assig ned   t ra ndom l gen e rated  nodes).  T he  P ge ner at io prof i le s   are all ocated  to  each PV  syst em  r and om ly  co nn ect e to  the   netw ork.     Power  flo si m ul a ti on are   perform ed  and  the  re quired   r esults  ( vo lt ag e s,  cu rr e nts,  a nd  power  flo ws are   extracte d an d r ecorde f or lat er  processin g u nd e t he  im pact assessm ent.   The  ab ove  pro cess  is  rep eat e unti def in ed  num ber   of   scenari os   are  s i m ulate unde each  P p enet rati on  le ve l.  Finall y,  the   resu lt of  al powe flo sim ulati on s   are   an al yz ed  an qua ntifie us i ng  s ever al   i m pact  m et rics  to  asse ss the  potenti al  so la P im pacts o n t he dist rib utio n netw ork.     2 . 2    Im pa c m etrics   The  f ollo wing   m et rics  are  adopted  t anal yz and   qu a ntify  the  key  te chn ic al   im pacts   of   diff e re nt   so la r PV  p e net rati on levels  on MV  d ist rib ution   netw orks .     2.2.1.  M ax im u m nodal  volt age   This m et ric ass esses the inf l ue nce of  so la P pe netrati on  on  t he  volt age  prof il es of the  netw ork.  F or  each  p ow e flo sim ulati on t he  m axi m u m   no dal   volt age  re corde in   the   ne twork   is  e xtra ct ed  to   exam ine  th e   var ia ti on  of  t he  m axi m u m  v oltage  unde r diffe ren PV pe netr at ion  lev el s.      2.2.2.  V olt ag e   v iola tio ns   This  m et ric  assesses  th pote ntial   vo lt age   rise  pro blem s   of  the   net w ork  due  to   hi gh  s olar  P V   pen et rati on  in  t he  f uture.   F or   each  powe fl ow  sim ulati on t he  daily   vo lt ag pro file of  al nodes   co nnect ed  t the  net wor ar analy zed  f or   adh e re nce  to   the  local   sta nda rd s Accord i ng  to  the   sta tuto r lim it sp eci fied  by   the  Ma la ysi an  el ect ric  util it y,    Tena ga    Nasio nal    Be r had    ( T NB)  for  vo lt ag es  in  MV   net w orks,  the  m axim u m   and   m ini m u m   vo lt age  le vels  are  1.0 an 0. 95pu,  res pecti ve ly Su bs e quen tl y,  the  per cent age  of  node th at   do  no c om ply  with  pe rm issi ble   volt age  li m it s   is  extracte t e xam ine  the  possible  volt age  rise   issues   unde r   diff e re nt P V p enetrati on leve ls.      2.2.3.  Sy s tem   po w er lo sses   This  m et ric  assesses  the   im pact  of  so la r   PV   pe netrati on  on  the   total   act ive  power  loss  of  th e   netw ork.   F or   each  powe fl ow   sim ulati on the  daily   aver age  act ive  powe loss  of  the  entire  netw ork  is  ob ta ine t de te rm ine  the  va riat ion   of   t he  syst e m   po we loss  at   di ff e re nt  so la P pe netrati on  le vel s.  Th e   equ at io n use t cal culat e t he sy stem  p ow e l os s is e xpre sse in  (1)     Av er ag e powe r l os s =   1  ( )  0   (1)     wh e re,      is t he  to ta l act ive pow e loss  of t he ne twork .     2.3.    Ca se  stu dy   The  propose Mon te   Ca rl im pact  assess m ent  fram ewo rk  was  a pp li ed   to  ty pical   Ma la ysi an  M V   distrib ution  ne twork   to  e xa m ine  the  f utur risk  of   so la P i nteg rati on.  Fi gure  pr ese nts  t he  one - li ne   diag ram   of   th e   te st  net work.  The  netw ork   is   th ree - phase  rad ia syst em   with  t wo  distri bu ti on   fee der s   serv i ng  range  of  resi den ti al com m ercial an industrial   load as   dep ic te i Fi gure   2.  T he  fee der a re  sup plied  by  a   30 M VA   distri bu ti on  tra ns f or m er  with  nom i nal  volt age  rat ing of   33 / 11kV.  Fee ders  a nd   of  the  ne twork  com pr ise   fiftee a nd  f ourtee node s,  res pec ti vely It  wa pr es um ed  that  the  m axi m u m   load  dem and   of  the   netw ork  t o be  6M W.  T he bas ic  n et w ork par a m et ers  of the  te st sy stem  are  su m m arized in  Table  1.     Evaluation Warning : The document was created with Spire.PDF for Python.
Ind on esi a J  E le c Eng &  Co m Sci     IS S N:  25 02 - 4752       Mo nte Carl o a na ly sis f or  s olar PV i mpact  asse ssme nt in  M V d ist rib utio n netw or ks   ( Dili ni Almei da )   27   2                             3                         4                                 5                           6                             7                             8                           9   1 0                     1 1                       1 2                         1 3 1 4                 1 5                         1 6 1 7                         1 8                         1 9                     2 0                         2 1                         2 2                       2 3                         2 4         2 5                       2 6                       2 7                         2 8 2 9                         3 0 Co mme r c i a l                                                 R e s i d e n ti a l                                                   I n d u s tr i a l 1 3 3 k V 1 1 k V     Figure  2 O ne - l ine d ia gr am  o f   the test  n et wor k       Table  1.   Ba sic   netw ork param et ers  of t he  te st  syst e m   Variable   Valu e   Cap acity   3 0 MVA    Vo ltag e   3 3 /1 1 k V   Nu m b e o f  f eeders   2   Nu m b e o f  no d es   30   Load   co n n ected   Feed er  1   3 .47 MW   Feed er  2   2 .53 MW       The  netw ork   par am et ers  rel at ed  to   the  i m ped ances  a nd  the   le ngths   of  each   li ne  s egm ent  wer acqu i red  f ro m   ty pical   Ma la ysi an  MV  dis tribu ti on  netw orks.  T he  nor m al iz ed  PV   ge ner at io a nd   load   consum ption   pro file adopte in  this  stud are  presente in  Fig ur 3.   A sh ow in  th figure,  resid entia consum ers  usual ly   ref le ct   low  co nsum ption   duri ng   t he   dayt i m e,  where  PV   ge nerat ion   is  at   it s   peak.  Ind us tria co nsum ers  pr ese nt  high  c on s umpti on   duri ng   t he   dayt i m and   ancho red   at   l ow e c on s um ption   durin the  ni gh t.  I co ntra st,  com m ercia con s um ers  de pict  sta ble  con s um ption  durin the  da yt i m e.   More ov e r,   P ge ner at io prof il co ncernin sunn cl i m atic  con diti on   w as  ad op t ed  to  in vestig at e   the   m axi m u m  p os s ible im pact du e to  high s olar PV  pe netrati on .           Figure  3 N orm al iz ed  PV   ge ner at io a nd lo ad  c onsu m ption   prof il es       Firstl y,  detai le netw ork  m od el   of   the  te st   syst e m   was  dev el oped  i th op e distri buti on   syst em   si m ulator   ( O pe nDSS)  sim u l at ion   platfo rm   interface wi th  MAT LAB  so ft war e I orde to   exam i ne  th e   influ e nce  of  s olar  PV  pe netr at ion   on  the   ne twork m ulti ple  PV  de plo y m ent  scena rios  we re  c onside r ed  by  Evaluation Warning : The document was created with Spire.PDF for Python.
                          IS S N :   2502 - 4752   Ind on esi a J  E le c Eng &  Co m Sci,   Vo l.   23 , N o.   1 Ju ly   20 21 23   -   3 1   28   integrati ng  te rand om l gen erate P c apacit ie at   ra ndom ly   gen er at ed  locat io ns  unde di ff e re nt  PV   pen et rati on   le vels.  I t he  re al iz at ion   of  P instal la ti on   capaci ty   for  e ach  pen et rati on  le vel,  t he  ra ndom ly  gen e rated  P capaci ti es  were  assigne unt il   the  cor res pond i ng   total   P V   capaci ty   in  ea ch  pe netrati on   le vel   was  reach ed The  m axi m u m   instal le s olar   PV  cap aci ty   was  li m i te to   12M at   100%   P pe netrati on   le vel. T he belo sel ect io ns   w ere c ounted  pri or to t he power  f lo sim ulati on s.     Sele ct ion   of th e PV pe netrati on level  (0 - 10 0%  in  ste ps   of 10%)     Ra ndom  selec t ion   of P i ns ta ll at ion  ca pacit y ( 0.1 - 2M W,  th r ee - phase  PV s yst e m s).     Ra ndom  selec t ion   of P l ocati on .   Ther ea fter th ree - ph a se,  ti m e - series,  m i nu te - by - m inu t e,  da il pow er  fl ow  sim ulati on wer e   perform ed  consi der in the  s unny  day  P gen e rati on  pro file   and   the  l oa co nsum ption   pro file s   co ve rin reside ntial co m m ercial and   industrial   loa ds   pr ese nted  i Fig ure  3.  T he  le vel  of   pe netrati on  wa el evated   from   0%  to  100%   in  ste ps   of   10%,  a nd  total   of   1000  Mo nte  Ca rl sim ulati on (10 sim ulati o ns   pe r   pen et rati on   le ve l)  wer c ondu ct ed  for  uni que   PV   de plo ym e nt  scena rios.   T he  sim ulatio r esults  wer a na ly zed  and  qua ntifie us i ng  the  de fine im pact  m et rics  to  ass ess  tw key  te chn ic al   c halle ng e of   so la P V   pen et rati on vo lt age r ise  a nd s yst e m  p ow e l os ses  on M V d ist ribu ti on  net works.       3.   RESU LT S  AN D  D I SCUS S ION   Fo eac pow er  flo sim ul at ion the  m axim u m   vo lt age,  the  pe rcen t age  of  no des  with  volt age   vio la ti ons,  a nd  the  daily   avera ge  po wer   l oss  values  we re  com pu te d.   T he reafter the  m e an  an th +/ -   one   sta nd a rd   de via ti on   of  eac 100  sim ulati on per   pe net ra ti on   le vel  were  cal culat ed.   Figures  4,   5,   and   gr a phic al ly   i ll us trat the  var ia ti on   of  the  m a xim u m   vo lt age the  per ce ntag of   node with  volt age  vio la ti on s,   and   the  da il aver age  pow er  loss  agains the  PV   penet rati on   le vel,   resp ect ively .   In   each  gra p h,  the  distrib utions  of   com pu te val ues  we re  re flec te by  scat te re points  unde each  pe netrati on  le vel.  Be side s,  the   m ini m u m  an t he  m axi m u m  v al ues of  t he dis tribu ti on  were  pr ese nted  usi ng a  square  m ark .           Figure  4 Va riat ion   of the m axi m u m  v oltage  aga inst t he  P V pen et rati on lev el       As  sho wn   in  Figure  4,  the  m axi m u m   no da vo lt age  rec orde f or   eac powe flo sim ula ti on   was   gr a dual ly   raised  with   the  i ncr ease pe ne trat ion   le vel.   Howe ver,  sim il ar  to  the  sc enar i with ou PV   integrati on,  t he   m axi m u m   no da vo lt age w as r ep or te near   t he  s ub sta ti on  a the  10% p ene trat ion   le vel T his  is   m ai nly  du to  the  fact  that  th low  P generati on   at   10%  pe netrati on   le vel  was  no su f fici ent  en ou gh   t com pen sat f or  the  vo lt ag dro du t the  higher  l oad   co ns um ption.  T he   u ppe volt ag lim i (i.e.,   1. 05pu)   vio la ti ons  were  sta rted  at   70 PV   pe netr at ion   le vel  an con ti nued  unti the  m axi m u m   le vel  con sid ered   in   the  stud y.   Acc ordin to  the  sim ula ti on   resu lt s,  there  we re  no  volt age  lim it  vio la ti on un ti the  PV   pen et r at ion   le vel  was  reac hed   60%.  Eve th ough  fe no des  wer ob s er ved   with  an  up per   volt age  li m it   vio la tio at   a   le vel  of  70%   P pen et r at ion,  the  num ber   of  vio l at ions  has  su r ge with  f ur t her   i ncr ease   in  pen et rati on   le vel.   Con se quently , th pe rce ntage  of  no des  with v oltage   vi olati on s   w as  inc rea sed  at   a   sig nifi cant  rate  a fter  7 0%   as   dep ic te d i n Fi gure  5.   Evaluation Warning : The document was created with Spire.PDF for Python.
Ind on esi a J  E le c Eng &  Co m Sci     IS S N:  25 02 - 4752       Mo nte Carl o a na ly sis f or  s olar PV i mpact  asse ssme nt in  M V d ist rib utio n netw or ks   ( Dili ni Almei da )   29       Figure  5 Va riat ion   of the  pe rc entage  of  node s w it h v oltage  vio la ti ons a gainst the  P V pene trat ion  le ve l           Figure  6 Va riat ion   of the  daily  av era ge p ower loss  a gainst  the P V pen et ra ti on  level       As  it   c ou l be   seen   in   Fig ure  6,  U - traj ec tory  was  de picte f or  the   vari at ion   of  the  da il aver ag e   powe loss  ag ai ns the  PV   pe netrati on  le ve l.  The  m ean  network   loss  w as  decr ease up  to  20%  pe ne trat ion  le vel  as  the  inf lux   of   P re du ced  the  net  power   flo withi the  net wor k.  However,  it   was  raise ab ove  the   base case  n et w ork  l os s at a  50 % p e netrati on l evel owin t t he  inc rease d re ver se  po wer fl ow.    By   exam ining   the  sim ulati on   res ults,  it   was   re vealed   that  vo lt age   rise   is  t he  worst  rep e rcu ssi on  of   high  so la PV   pen et rati on   on  the  netw ork.   B esi des,   it   was  r e m ark ed  t hat  the  occ urre nce  of   vo lt ag vio l at ions   do e not  s olely   rely   on  the   tot al   instal le PV  capaci ty   of  th netw ork  a t he  re su lt ob ta ined  from   the  Mon te   Ca rlo  analy s is  hav s how no ti ceable   distrib ution   at   higher  P pe netra ti on   le vels  (fo r   each  si m ulatio at   a   certai pe netra ti on   le vel,  the  instal le capac it was  retai ned   as  co ns ta nt).   T her e fore,  it   cou l be  c onc lud e that  the  siz a nd  locat io of  PV   syst e m conne ct ed   to  t he  netw ork  play   i m po rtant  r oles  in  a ssessi ng  the   i m pact  of   PV   pen et rati on   on  distrib ution   net works.  F ur t her,  the  stoc hastic   natu re  of  so la PV   ge ner at io a nd   load  c onsu m ption   c ou l hav e   substanti al   im pact  on   the  oc currence   of  th ese  volt age  vi ol a ti on s.  In   par ti cular ,   the  occurre nce   of   volt age  vi ol at ion co uld   increase  wh e m os PV   syst e m s   are  locat ed  at   residen ti al   nodes wh e re  co nsum ption   is  lo duri ng   high  P V   gen e rati on  pe rio ds Co nvers el y,  the  load  c on s um ption   of   bo t com m ercial   a nd   in dustria nodes  re flect go od  co rr e la ti on   with  P ge ner at io n,  wh ic c ou l te nd   to   m ini m iz e the occ urren ce  a nd the m agn it ude  of volt age  vio l at ion s.       4.   CONCL US I O NS   In   this  pa per Mon te   Ca rl based   pro ba bili sti fr a m ework   was  propos ed  to  asses th so la P V   i m pacts  on  M distrib utio netw orks  w hile  cat erin to   t he  un ce rtai nties  relat ed   to  s ol ar  P ge ner at i on.  T he  case  stud pre sented  has  il lustrate the  ef f ect iveness  of   Mon te   Ca rlo  a naly sis  in  evaluati ng   t he  pote ntial   te chn ic al   chall eng e on  ty pi cal   Ma la ys ia MV  di strib ution   netw ork.   T he  m axi m u m   n od al   vo lt age volt age   Evaluation Warning : The document was created with Spire.PDF for Python.
                          IS S N :   2502 - 4752   Ind on esi a J  E le c Eng &  Co m Sci,   Vo l.   23 , N o.   1 Ju ly   20 21 23   -   3 1   30   vio la ti ons,  an syst e m   po wer  losses  we re  use as  key  m e tric to  qu a ntif the  te chn ic al   per f orm ance  of   the   distrib ution   net work.  T he  sim ulati on   resu lt rev eal e that  the  w orst  rep e r cussion  of   high   so la r   P pen e trat ion  on   ty pical   Ma l ay sia MV  distribu ti on  netw orks  is  the  vio l at ion   of   the   up per   volt age  sta tutor li m it wh ic is   li kely   to  occur   beyo nd  70%  pen et rati on   le ve l.  More over,  it   was  identifi ed  that  the  siz and  locat io of   P V   syst e m con ne ct ed  to  t h network  play   c ru c ia ro le in   ass essing  the  im pact  of   PV  pe ne trat ion   on  distr ibu ti on   netw orks.  F urt her m or e,  it   was  obser ve that  the   sto chasti nat ur e   of  s olar  PV  ge ner at io a nd  lo a consum ption   c ou l ha ve  a substanti al   i m pact  on   t he  occ urr ence o vo lt a g e   vio la ti ons.  T hus,  it   is  im per at ive  to   consi der   t he  stocha sti beh a vi our  of  PV   outpu a nd   l oad   dem and to geth er  with  the  fee der   c har act e risti cs  in   the p la nn i ng st age,  t o defi ne  t he  li m it s f or P instal la ti ons  in or der  t o pr e ven f utu re  ove rvolta ge pr oble m s.       ACKN OWLE DGE MENTS   The  aut hors  w ou l li ke  to  th ank   Un i ver sit Tena ga  Nasi onal   (U N ITE N),  Ma la ysi fo pro vid i ng   th e   BOLD  Re sea r ch  G ran (RJ O10 517844/0 71)  to  do   this  r esearch  a nd  ackno wled ge  th publica ti on   su pp or t   thr ough J 5100 50002 - BO LD REFR ESH2 025 - CE NTRE  O F EX C ELLE N CE from  the iRM C of   U NI T E N .       REFERE NCE S   [1]   A.  Colm ena r - Sa ntos,  C.   Rei no - Rio,   D.   Borge - Diez ,   and   E .   Co ll ado - Fern ánde z ,   Distribut ed  g e ner ation:  r eview  of  fac tors  tha c an  cont ribu te   m ost  to  ac hie v sce nar io  of  DG   unit embedde in  the   ne distribution  net works ,   Re newab le   and   Sustainabl e   Ener gy  Revi ews ,   vol .   59 ,   pp .   1130 1 148,   2016 ,   doi 1 0. 1016/j.rser. 201 6. 01. 023 .   [2]   E.   Karuna r at hn e,   J.  Pasupulet i ,   J.  Eka na y a ke ,   and  D.  Alm ei da,   Optimal  Placem ent   and  Sizi ng  of  DG i Distribut ion  N et works   U sing  MLPS Alg orit hm , ”  En ergie s ,   vol .   13,   no.   23,   p .   6 185,   2020 ,   do i:  10. 3390/e n1323 6185 .   [3]   REN21,   Renew abl es  2020  Global   Stat us  Rep ort, ”  Par is,  Fran ce ,   2020 Acce ssed:  June  24,   2021.  [ On li ne ] .   Av ai la ble:  http s:/ /www . re n21.net/ gsr - 2020/ .   [4]   SEDA ,   Annual   Repor 2018, ”  Sus ta ina ble   Energ y   Deve lopment  Authority   (SEDA ),   Malay s ia,  2019.   Acc essed :   June  2021.   [Onli ne] .   Available:  h tt p://ww w.seda . gov. m y /pdfdow nloa d/ annua l - r ep ort - 2018/ ? wpdm dl= 2080   [5]   E.   Karun ara thn e ,   J.  Pasupule ti ,   J.  Eka n a y ak e,   an D.  Alm ei da,  Network  loss  red uct ion  and  voltage  improvem en t   b y   opti m al   pl ace m ent   and  sizi ng  of  distri bute ge ner at ors  with  act ive   and  react iv power  inj ection  using  fine - tune PS O,”   Indone sia Journal  of   E l e ct rical  Eng ine er ing  and  Comput er  Scienc e   ( IJE E CS) ,   vol.   21 ,   no.   2,   pp .   646 - 656 ,   2021 ,   doi 10 . 11 591/i jeec s.v21 . i 2. pp647 - 656 .   [6]   C.   Long ,   J.  W u,   L.   Thomas,  and   N.  Jenkins,   Optimal  ope ration   of  soft  open   po int in  m edi um   volt ag e le c trica distri buti on  n etw orks  with  distri bute g ene r at ion , ”  App lied   Ene rgy ,   vo l.  184,   pp.   427 437,   2016 ,   d oi:  10. 1016/j.a p ene r g y . 2016. 10 . 031 .   [7]   D.  W .   Alm ei da ,   A.  H.  M.   S.  M.  S.  Abe y sing he,   and  J.  B .   E kanay ak e,  Anal y sis  of   rooftop   solar  impa ct on  distri buti on   ne tworks,”   Ce yl on   Jo urnal  of  S ci en ce ,   vol. 48, no. 2, p p 103 112 ,   201 9 ,   doi 10 . 4038/cjs.v48i 2. 7614 .   [8]   D.  Alm ei da,   J.  Pasupulet i ,   J.  Ekana y ake,  and  E.   Karuna rat hn e,   Miti gation  of  over voltage  due  to   high  pene tra t io n   of  solar  photovol ta i cs  using  smart   inve r te rs  vol t/ var   con trol,”   I ndonesian  Journal  of  El e ct rica l   Engi nee ring  an Computer  Scien ce   ( IJEECS) ,   vol .   19 ,   no .   3 ,   pp .   1 259 1266,   2020 ,   doi:  10. 11591 / ijeec s.v19 . i3 . pp12 59 - 1266 .   [9]   D.  Alm ei da ,   S.  Abe y singhe ,   M.   P.  Ek anay ak e,  R.   I.  Godaliy a d da,   J.   Ek anay ak e,   and  J.  Pasupu le ti,  Gene r al i zed   appr oac to  asse ss   and  cha racte r ise  the   impa ct   o solar  PV   on  LV  net works , ”  In te rnational   Journal  of  Elec tri cal  Powe and  Ene r gy  Syst ems ,   vol .   121,   p .   106058 ,   2020 ,   doi 10 . 10 16/j . i je p es. 2020. 106058 .   [10]   A.  S.  N.  Huda  and  R.   Ži v anov ,   La rg e - sca le  int egr a ti on  of  distri bute g enerat ion  in to  distribution  net works :   Stud y   obj ec t ive s ,   rev i ew  of  m odel and  computa t iona too ls,”   Re n ewabl and  Sust ainabl En ergy   Re v ie ws ,   vol .   76 ,   pp.   974 988 ,   20 17 ,   doi 10 . 1016 /j . rser . 2017. 03 . 0 69 .   [11]   D.  Alm ei da,   J.  P asupule ti,  and  J .   Eka na y a ke ,   Perform anc ev al ua ti on  of  PV   pene t rat ion  at   diff ere n locati ons  in  a   LV  distri buti on  net work  conne c t ed  with  a off - loa ta cha ng in tra nsform er, ”  Indone sian  Journal  of  El ectric al   Engi ne ering  and   Computer  Sci en ce   ( IJE ECS) ,   vol .   21,   no.   2 ,   pp.   9 87 - 993,   2020 ,   d oi:   10. 11591/ ijee cs. v21. i2 . pp987 - 993 .   [12]   V.  H.  M .   Que zada,   J.  R .   Abbad ,   and  T .   G .   S .   R om án,   As sessm ent   of  en erg y   d i stribut ion  losses  for  inc r ea sing  pene tr at ion  of  di stribut ed  generation,”   in   IEEE  Tr ansacti ons  on  Powe Syste ms ,   vol.   21,   no.   2,   pp .   533 - 540,   Ma y   2006,   doi 10 . 11 09/T PW RS . 2006. 873115 .   [13]   M.  Thoms on  and  D.  G.  Infie ld ,   Network  Pow er - Flow  Anal y sis  for  High  Penet r at ion  of  Distribu te Gene r at ion ,   in   IE EE   Tr ansacti ons  on  P ower  Syste ms ,   vol.   22 ,   n o.   3,   pp.   1 157 - 1162,   Au g.   2007,   doi:  10. 1109/T PW RS . 2007. 901290 .   [14]   G.  v ar,   A .   m ez - Expósit o,   A.  Arcos - Varga s,   and  M.  Rodrígu e z - Montañé s,   Influe nc of  roof to PV   gene ration  on  net   demand ,   losses  and  netw ork  conge stion s:  ca se  stud y , ”  Inte rnat ional   Journal  of  E lect rical   Pow er  and   Ene rgy  S yste ms ,   vol.   106 ,   pp .   68 86,   2019 ,   doi 10 . 1016/j.ijepe s.20 18. 09. 013 .   [15]   V.  Annathur ai,  C.   K.  Gan,   M.  R.   A.  Ghani,   a nd  K.  A.  Baha r in ,   Im pac ts  of  solar  var ia b il i t y   on  distri bu ti o n   net works   per for m anc e,”  Int ernati onal  Journal  o Applied   Engi n ee ring  Re search ,   vol .   12 ,   no .   7 ,   pp.   1151 115 5 ,   2017.   Evaluation Warning : The document was created with Spire.PDF for Python.
Ind on esi a J  E le c Eng &  Co m Sci     IS S N:  25 02 - 4752       Mo nte Carl o a na ly sis f or  s olar PV i mpact  asse ssme nt in  M V d ist rib utio n netw or ks   ( Dili ni Almei da )   31   [16]   A.  Te m iz   et   al. ,   Inve stiga ti on  of   MV  Distribut ion  Networks  with   High - Penet ratio Distribut ed   PV s:  Study   for  an   Urban  Area ,   in   Ene rgy  Proce d ia ,   2017 ,   vol .   141 ,   pp.   517 524 ,   doi 10. 1016 /j.e g y p ro. 2017. 11 . 069 .   [17]   M.  Hashem ina m in,   V.  G.  Age lidis   and  A.  Heid ari ,   Im pac stu d y   of  high   PV   pene tr at ion  in  l ow  and  m edi um - volt ag net work when  considering  res ide n tial  and  industrial  l oad  profile,   20 13  Inte rnationa l   Confe renc on   Re newab le   E nergy   Re sear ch  and  A ppli cations  (ICRE R A) ,   2013,   pp.   3 47 - 352,   doi:  10. 1109/ICRE R A.2013. 6749778 .   [18]   N.  Mourad  and  B.   Moham ed,   Im pac of  inc rea sed  distri bu ted  photovol ta i gene ra ti on  on  rad ial  distri but io net works ,   2016   Inte rnational   C onfe renc on  Control,   Dec ision  and  Information  Technol ogie ( CoDIT) ,   2016,   pp.   292 - 295,   doi 10 . 1109/CoDIT. 20 16. 7593576 .   [19]   J.  Garc ia - Vil la lo bos,  P.  Egui a,   E.   Torr es  and  A.  Et xega r ai,  Im pac of   photovol ta i self - cons um pti on  on  po wer   losses  and  volta ge  le v el of  M distri buti on  n et works ,   2017  IEE PE Inno vat i ve   Smar Gr id  Technol og ie s   Confe renc E urope  ( ISGT - Euro pe) ,   2017,   p p.   1 - 6,   doi 10 . 1109/I SG TE urope . 201 7. 8260161 .   [20]   Y.  Ta ng,   R .   Burgos,  C.   Li   and  D .   Boro y evich,   Im pac of  PV   inve rte g ene ra ti on  on  volt age   prof ile  and  power  loss  in  m edi um   volt age   distri buti on  sy stems ,   2016  IEE 43rd  Phot o vol taic  Speciali s ts  Confe renc ( PV SC) ,   2016,   pp.   1869 - 1874,   doi : 10.1109/ PV SC . 2016. 7749946 .   [21]   E.   Zi o ,   M.  Del f ant i ,   L.   Giorg i ,   V.  Olivi e ri,   a nd  G.  Sansavini ,   Monte  Carlo  sim ula ti on - ba sed  proba bil ist i assess m ent   of  D pene tration  in   m edi um   volt age   distri buti on  n etw orks,”   Inte rnat ional   Journal  o El e ct rica Power   and  Ene rgy   Syst ems ,   vol .   64 ,   pp .   852 860,   2015 ,   doi 10 . 1016/j.ijepe s.2014. 0 8. 00 4 .   [22]   N.  Srisaen  an A.  Sangsw a ng,   Eff e ct o PV   Grid - Connec t ed  S y s tem   Loc ation  o Distr ibut i on  S y stem,   AP CC AS  2006  -   2006  IEE Asia  Pacifi Con fe ren ce   on  Circui ts  and  Syste ms ,   2006,   pp.   852 - 855,   do i :   10. 1109/APC CAS . 2006. 34217 5 .   [23]   R.   Vill a faf i la,  S.   Galc er an,   B .   B ak - Jensen,   P.  C hen,   Z .   Chen  an S.  Sorensen,   Probabil isti assess m ent   of  wind   power  produc ti o on  volt age   pr ofil in  distri bu t ion  net works ,   2 007  9th  Inte rnational   Confe ren c on  El ectric al  Powe r Quali ty   a nd  Util isat ion ,   2 007,   pp .   1 - 5 ,   doi 10. 1109 /E PQ U.2007. 4424123 .   [24]   M.  Begovi c ,   A.   Prege lj,  A.  Ro hat gi  and  D.  N ovosel,   Im pac t   of  ren ewa b le   distri bute d   generat ion  on  pow e r   s y stems ,   Proceedi ngs of   the   34 t Annual   Hawai Inte rnationa C onfe renc on  Sys te m Sc i enc es ,   20 01,   pp.   6 54 - 663,   doi:   10 . 1109/HI CS S.2001. 926265 .   [25]   M.  Delf anti,   D .   Fala bre tt i ,   and  M.  Merlo,   Dispersed  generatio impact   on  dist ribut ion  n et work   losses,”   Elec tri Powe r Sy st ems  Re search ,   vol .   9 7,   pp .   10 - 18 ,   20 13 ,   doi 10 . 1016 /j . epsr. 2012. 11 . 0 18 .   [26]   Q.  Tra n,   M.  Cong  Pha m ,   L.   Pare nt  and  K.  Sousa ,   Inte gra ti on  o PV   Sy stems   int Grid:  From   I m pac Anal y s is  to  Soluti ons,   2018   IEEE  In te rnati onal  Confe renc on  Env ironment  and  El ectri cal   Engi ne ering   and  2018  IEE Industrial  and  Comm erc ial   Powe Syste ms   Europe  ( EE EIC  I & CPS  Eu rope) ,   2018,   pp.   1 - 6,   doi :   10. 1109/EEEIC.2018.8494400 .     Evaluation Warning : The document was created with Spire.PDF for Python.