IEEE Personl use of this mteril is permitted Permission from IEEE must e otined for ll other uses, in ny urrent or future medi, inluding reprinting/repulishing this mteril for dvertising or promotionl purposes, reting new olletive works, for resle or redistriution to servers or lists, or reuse of ny opyrighted omponent of this work in other works rxiv:36798v [ssy] 6 Nov
Effiient Computtion of Sensitivity Coeffiients of Node Voltges nd Line Currents in Unlned Rdil Eletril Distriution Networks Konstntin Christkou, Memer, IEEE, Jen-Yves Le Boude, Fellow, IEEE, Mrio Polone, Senior Memer, IEEE, Dn-Cristin Tomozei, Memer, IEEE Astrt The prolem of optiml ontrol of power distriution systems is eoming inresingly ompelling due to the progressive penetrtion of distriuted energy resoures in this speifi lyer of the eletril infrstruture Distriution systems re, indeed, experiening signifint hnges in terms of opertion philosophies tht re often sed on optiml ontrol strtegies relying on the omputtion of linerized dependenies etween ontrolled (eg voltges, frequeny in se of islnding opertion) nd ontrol vriles (eg power injetions, trnsformers tp positions) As the implementtion of these strtegies in rel-time ontrollers imposes stringent time onstrints, the derivtion of nlytil dependeny etween ontrolled nd ontrol vriles eomes non-trivil tsk to e solved With referene to optiml voltge nd power flow ontrols, this pper ims t providing n nlytil derivtion of node voltge nd line urrent flows s funtion of the nodl power injetions nd trnsformers tp-hngers positions Compred to other pprohes presented in the literture, the one proposed here is sed on the use of the [Y] ompound mtrix of generi multi-phse rdil unlned network In order to estimte the omputtionl enefits of the proposed pproh, the relevnt improvements re lso quntified versus trditionl methods The vlidtion of the proposed method is rried out y using oth IEEE 3 nd 34 node test feeders The pper finlly shows the use of the proposed method for the prolem of optiml voltge ontrol pplied to the IEEE 34 node test feeder Index Terms Voltge/urrent sensitivity oeffiients, unlned eletril distriution networks, power systems optiml opertion, smrt grids I INTRODUCTION O ptiml ontrols of power systems re often sed on the solution of liner prolems tht link ontrol vriles to ontrolled quntities y mens of sensitivity oeffiients Typil optimiztion prolems refer to sheduling of genertors, voltge ontrol, losses redution, et So fr, these tegories of prolems hve een ommonly investigted in the domin of high voltge trnsmission networks However, during the pst yers, the inresed penetrtion of distriuted energy resoures (DERs) in power distriution systems hs rised the importne of developing optiml ontrol strtegies speifilly pplied to the opertion of these networks (eg [], [], [3], [4], [5], [6]) Within this ontext, it is worth noting tht the solution of optiml prolems eomes of interest only Konstntin Christkou, Mrio Polone, Jen-Yves Le Boude nd Dn- Cristin Tomozei (emil: konstntinhristkou@epflh, jen-yvesleoude @epflh, mriopolone@epflh, dn-ristintomozei @epflh) re with the Éole Polytehnique Fédérle de Lusnne, CH-5 Lusnne, Switzerlnd if it meets the stringent time onstrints required y reltime ontrols nd imposed y the higher dynmis of these networks ompred to the trnsmission ones Typil exmples of optiml ontrols tht re not yet deployed in tive distriution networks (ADNs) re voltge nd power flow ontrols Usully this tegory of prolems hs een ddressed in the literture y mens of liner-pprohes pplied to the dependeny etween voltges nd power flows s funtion of the power injetions (eg [4], [5], [7], [8]) The typil pproh for the solution of this lss of ontrol prolems is the use of the sensitivity oeffiients through n updted Join mtrix derived from the lod flow prolem [9], [], [], [], [3] However, from the omputtionl point of view, the min disdvntge of suh tegory of methods is tht, for every hnge in the opertion onditions of the network, n updted Join mtrix needs to e uilt on the sis of the network stte nd needs, then, to e inverted This proedure involves non-trivil omputtion onstrints for the implementtion in rel-time entrlized or deentrlized ontrollers For this reson, the uthors of [4] hve proposed the diret omputtion of voltges nd network losses sensitivity oeffiients, sed on the Guss-Seidel formultion of the lod flow prolem, y mking use of the [Z] mtrix of lned network Also, in [7] it hs een proposed the use of the [Z] mtrix long with the onstnt-urrent model for lods nd genertors In [8] the sensitivity oeffiients re proposed to e lulted strting from the network rnh urrents The pproh presented in [5], [6], [7], [8], [9] elongs to lss of methods typilly derived from iruit theory nd is sed on the use of the so-lled djoint network In order to inrese the omputtionl effiieny of this tegory of pprohes, nd to extend it to the inherent multiphse unlned onfigurtion of distriution networks, the min ontriution of this pper is to provide strightforwrd nlytil derivtion of node voltges nd line urrents sensitivities s funtion of the power injetions nd trnsformers tp-hngers positions To this end, we propose to use the so-lled [Y] ompound mtrix, whih hs the dvntge of eing sprse Compred to [7] the pproh here proposed tkes into ount the whole dmittne mtrix of the network On the other hnd the nlytil derivtion of sensitivities in [7] ws sed on the pproximted representtion of the network
3 lines where lines shunt prmeters re negleted The method presented in [8] lwys requires se-se lod flow solution nd it relies on the ssumption tht ll genertors re PV nodes (ie with fixed voltge mgnitude) Also, it does not ount for the mutul oupling etween different phse ondutors The pproh tht ppers the more generl mong the ove listed is the one proposed in [4] However, this method depends on pseudo-lod flow pproh (ie it mkes use of Guss-Seidel itertive proess with fixed numer of itertions) whih influenes the ury of the omputed oeffiients Furthermore, ompred to [4] we hve een le to: generlize the prolem formultion for generi numer of slk usses; extend the omputtion of sensitivities to tp-hngers positions (ie hnges of slk usses referene voltges); provide the proof tht the nlytil omputtion of sensitivities dmits unique solution for the se of rdil networks nd tke into ount the inherent multiphse nd unlned nture of distriution networks The struture of the pper is the following: Setion II fouses on the prolem formultion y desriing, in detil, the nlytil proedure t the se of the proposed method It lso inludes proof of uniqueness of the solution of the liner system used to lulte the sensitivity oeffiients for the se of rdil networks The sme setion lso provides omputtionl ost nlysis of the proposed method versus trditionl pprohes Setion III vlidtes the proposed method using the IEEE 3 nd 34 node test feeders Setion IV shows n pplition exmple of sensitivity oeffiients relted to the optiml voltge ontrol in unlned distriution networks tking dvntge of the possiility of omputing them for ll the phses Setion V provides the finl remrks out possile pplitions of the proposed method II PROBLEM FORMULATION A Clssil Computtion of Sensitivity Coeffiients in Power Networks In this prgrph we mke referene to lned network omposed y K usses Trditionlly, there re three proposed wys to lulte the sensitivity oeffiients of our interest The first method onsists of estimting them y series of lod flow lultions eh performed for smll vrition of single ontrol vrile (ie nodl power injetions, P l, Q l ) [4]: Ēi Ēi = Ēi P l = Ēi Q l Pi,i l = Q i,i l = Pi,i l = Q i,i l = Īij Īij = Īij P l = Īij Q l Pi,i l = Q i,i l = Pi,i l = Q i,i l = It is importnt to oserve tht line shunt prmeters re non-negligile in se of networks hrterized y the presene of oxil les These types of omponents re typil in the ontext of urn distriution networks In the rest of the pper omplex numers re denoted with r ove (eg Ē) nd omplex onjugtes with r elow (eg Ē) () where Ēi is the diret sequene phse-to-ground voltge of node i nd Īij is the diret sequene urrent flow etween nodes i nd j (i, j { K}) The seond method uses the Newton Rphson formultion of the lod flow lultion to diretly infer the voltge sensitivity oeffiients s sumtries of the inverted Join mtrix (eg [9], [], [], [], [3]): J = P Ē Q Ē P θ Q θ () It is worth oserving tht suh method does not llow to ompute the sensitivities ginst the trnsformers tp-hngers positions Additionlly, s known, the sumtrix Q is Ē usully dopted to express voltge vritions s funtion of retive power injetions when the rtio of longitudinl line resistne versus retne is negligile It is worth noting tht suh n ssumption is no longer pplile to distriution systems tht require in ddition to tke into ount tive power injetions A third method is derived from iruit theory In this method Tellegen s theorem is pplied in power networks nd the omputtion of sensitivities relies on the onept of the so-lled djoint networks (eg [5], [6], [7], [8], [9]) This pproh requires se-se lod flow solution in order to uild speifi djoint network tht needs to e solved in order to infer the desired sensitivities B Anlytil Derivtion of Voltge nd Current Sensitivity Coeffiients This susetion ontins the min nlytil development of this pper relted to the derivtion of the voltge sensitivity oeffiients 3 ) Voltge Sensitivity Coeffiients: the nlysis strts with the voltge sensitivity oeffiients To this end, we derive mthemtil expressions tht link us voltges to us tive nd retive power injetions For this purpose, K-us 3- phse generi eletril network is onsidered The following nlysis trets eh phse of the network seprtely nd, thus, it n e pplied to unlned networks As known, the equtions tht link the voltge of eh phse of the usses to the orresponding injeted urrent re in totl M = 3K nd they re given y: [Ī ] = [Ȳ ] [Ē ] (3) where [Ī ] = [Ī, Ī, Ī, ĪK, ĪK, ĪK ] T, [Ē ] = [Ē, Ē, Ē, ĒK, ĒK, ĒK ] T We denoted y,, the three network phses The [Ȳ ] mtrix is formed y using the so-lled ompound dmittne mtrix (eg []) s 3 As shown in susetion II-B the urrent sensitivities n e strightforwrdly derived from the voltge ones
4 follows: ] [Ȳ = Ȳ K Ȳ K Ȳ K Ȳ K Ȳ K Ȳ K Ȳ K Ȳ K Ȳ K In order to simplify the nottion, in wht follows we will ssume the following orrespondenes: [Ī ] = [Ī,, ĪM ] T, [Ē ] = [Ē,, ĒM ] T nd Ȳ M ] [Ȳ = Ȳ M Ȳ MM For the rest of the nlysis we will onsider the network s omposed y S slk usses nd N usses with P Q injetions, (ie {,, M} = S N, with S N = ) The P Q injetions re onsidered onstnt nd independent of the voltge In this respet, we re ssuming tht for eh seprte perturtion of nodl power injetions, the other lods/genertors do not hnge their power set points Therefore, the omputtion of the sensitivities inherently ounts for the whole response of the network in terms of vrition of oth tive nd retive power flows Suh onsequene llows to ompute the sensitivities in the lose viinity of the network stte The link etween power injetions nd us voltges reds: i S = Ē i Ȳ ij Ē j, i N (4) The derived system of equtions (4) holds for ll the phses of eh us of the network Sine the ojetive is to lulte the prtil derivtives of the voltge mgnitude over the tive nd retive power injeted in the other usses, we hve to onsider seprtely the slk us of the system As known, the ssumptions for the slk us equtions re to keep its voltge onstnt nd equl to the network rted vlue, y lso fixing its phse equl to zero Hene, for the three phses of the slk us, it holds tht: Ēi =, i S (5) At this point, y using eqution (4) s strting point one n derive losed-form mthemtil expressions to define nd quntify voltge sensitivity oeffiients with respet to tive nd retive power vritions in orrespondene of the N usses of the network To derive voltge sensitivity oeffiients, the prtil derivtives of the voltges with respet to the tive nd retive power P l nd Q l of us l N hve to e omputed The prtil derivtives with respet to tive power stisfy the following system of equtions: {i=l} = Ē i Ȳ ij Ē j Ē i Ȳ ij Ēj (6) where it hs een tken into ount tht: i S = {P i jq i } = {i=l} (7) The system of equtions (6) is not liner over omplex numers, ut it is liner with respet to Ēi, Ē i, therefore it is liner over rel numers with respet to retngulr oordintes As we show next, it hs unique solution nd n therefore e used to ompute the prtil derivtives in retngulr oordintes to redue the omputtionl effort A similr system of equtions holds for the sensitivity oeffiients with respet to the injeted retive power Q l With the sme resoning, y tking into ount tht: we otin tht: j {i=l} = Ē i i S = {P i jq i } = j {i=l} (8) Ȳ ij Ē j Ē i Ȳ ij Ēj (9) By oserving the ove liner systems of equtions (6) nd (9), we n see tht the mtrix tht needs to e inverted in order to solve the system is fixed independently of the power of the l-th us with respet to whih we wnt to ompute the prtil derivtives The only element tht hnges is the left hnd side of the equtions One Ēi, Ē i re otined, the prtil derivtives of the voltge mgnitude n e expressed s: Ēi = Ēi Re(Ē Ēi i ) () nd similr equtions hold for derivtives with respet to retive power injetions Theorem : The system of equtions (6), where l is fixed nd the unknowns re Ēi, i N, hs unique solution for every rdil eletril network The sme holds for the system of equtions (9), where the unknowns re Ēi, i N Proof: Sine the system is liner with respet to retngulr oordintes nd there re s mny unknowns s equtions, the theorem is equivlent to showing tht the orresponding homogeneous system of equtions hs only the trivil solution The homogeneous system n e written s: = i Ȳ ij Ē j Ē i Ȳ ij j, i N () where i re the unknown omplex numers, defined for i N We wnt to show tht i = for ll i N Let us onsider two eletril networks with the sme topology, ie sme [Ȳ ] mtrix, where the voltges re given In the first network, the voltges re Ē i = Ēi, i S Ē i = Ēi i, i N nd in the seond network they re Ē Ē i = Ēi, i S i = Ēi i, i N () (3)
5 Let i e the onjugte of the sored/injeted power t the ith us S in the first network, nd i in the seond Apply eqution (4) to us i N in the first S network: i = Ē S i Ȳ ij Ē j ( = (Ē i i ) Ȳ ij Ē j ) Ȳ ij (Ēj j ) j S = Ē i Ȳ ij Ē j i Ȳ ij j i Ȳ ij Ē j Ē i Ȳ ij j Similrly, for the seond network nd for ll usses i N : i = Ē S i Ȳ ij Ē j i Ȳ ij j i Ȳ ij Ē j Ē i Ȳ ij j Sutrt the lst two equtions nd otin ( ) i S S i = i Ȳ ij Ē j Ē i Ȳ ij j By eqution (), it follows tht S i for ll i N Thus i = S the two networks hve the sme tive nd retive powers t ll non slk usses nd the sme voltges t ll slk usses As disussed in [] for rdil distriution networks suh n ssumption mens tht the lod flow prolem lwys hs unique solution Therefore, it follows tht the voltge profile of these networks must e extly the sme, ie Ē i i = Ē i i for ll i N nd thus i = for ll i N ) Current Sensitivity Coeffiients: From the previous nlysis, the sensitivity oeffiients linking the power injetions to the voltge vritions re known Thus, it is strightforwrd to express the rnh urrent sensitivities with respet to the sme power injetions Assuming to represent the lines tht ompose the network y mens of π models, the urrent flow Ī ij etween nodes i nd j n e expressed s funtion of the phse-to-ground voltges of the relevnt i, j nodes s follows: Ī ij = Ȳij(Ēi Ēj) (4) where Ȳij is the generi element of [Ȳ ] mtrix etween node i nd node j Sine the voltges n e expressed s funtion of the power injetions into the network usses, the prtil derivtives of the urrent with respet to the tive nd retive power injetions in the network n e expressed s: Īij = Ēi Ȳij( Ēj ) Īij Ēi = Ȳij( Ēj ) (5) Applying the sme resoning s erlier, the rnh urrent sensitivity oeffiients with respet to n tive power P l n e omputed using the following expressions: Īij = Īij Re(Ī Īij ij ) (6) Similr expressions n e derived for the urrent oeffiients with respet to the retive power in the usses s: Īij = Īij Re(Ī Īij ij ) (7) C Sensitivity Coeffiients with respet to tp positions of trnsformers This susetion is devoted to the derivtion of nlytil expressions for the voltge sensitivity oeffiients 4 with respet to tp positions of trnsformer We ssume tht trnsformers tp-hngers re loted in orrespondene of the slk usses of the network s for distriution networks these represent the onnetions to externl trnsmission or su-trnsmission networks As onsequene, the voltge sensitivities s funtion of the tp positions re equivlent to the voltge sensitivities s funtion of the slk referene voltge We ssume tht the trnsformers voltge vritions due to tp position hnges re smll enough so tht the prtil derivtives onsidered in the following nlysis re meningful Furthermore, we ssume tht the power injetions t the network usses re onstnt nd independent of the voltge With the sme resoning s in Se II-B, the nlysis strts in eqution (4) We write Ēl = Ēl e jθ l for ll usses l For us i N the prtil derivtives with respet to the voltge mgnitude Ēk of slk us k S re onsidered: Ē i Ȳ ik e jθ k = ik Ȳ ij Ē j Ē W i Ȳ ij Wjk, (8) where W ik := ( Ēi Ēk = Ēi Ēi Ēk j θ i Ēk ) Ē i, i N We hve tken into ount tht: Ȳ ij Ē j = Ēk Ȳike jθ k (9) nd j S i S = () Ēk The derived system of equtions (8) is liner with respet to ik W nd W ik, nd hs the sme ssoited mtrix s the system in (6) Sine the resulting homogeneous system of equtions is identil to the one in (), y Theorem it hs unique solution After resolution of (8), we find tht the sensitivity oeffiients with respet to the tp position of the trnsformer t us k re given y ( ) Ēi Wik Ēk = Ēi Re () 4 Note s shown in Se II-B one the voltge sensitivities re otined the ones of urrents n e omputed diretly Ē i
6 D Computtionl Cost Anlysis for Voltge Sensitivities with respet to PQ injetions The im of this susetion is to show the omputtionl dvntge of the proposed method ompred to the lssil pproh with respet to the omputtion of voltge sensitivities s funtion of power injetions only 5 Furthermore, the two methods re pplied to the IEEE 3 nd 34 node test feeders nd ompred in terms of CPU time neessry to lulte the voltge sensitivity oeffiients We re ssuming tht: ) there re lods/injetions in ll three phses of the system nd ) the phsors of phse-to-ground voltges in ll the network re known (eg oming from stte estimtion proess []) In the following tle, Algorithm shows the steps required to lulte the voltge sensitivity oeffiients using the trditionl method nd Algorithm shows the orresponding steps using the nlytil method proposed here For the trditionl method n updted Join needs to e uilt, nd its inverse will provide the desired voltge sensitivities For the nlytil method the orresponding steps refer to invert squre mtrix of size N (s reported in Setion II-B N refers to the numer of network usses with P Q injetions) nd multiply the inverse mtrix with one olumn vetor for eh P Q us in the network Algorithm Computtion of voltge sensitivity oeffiients using the Join method : uild Join mtrix ssoited to the Newton Rphson method : invert mtrix J of size N N 3: extrt the su-mtries orresponding to the desired sensitivity oeffiients Algorithm Computtion of voltge sensitivity oeffiients using the nlytil method : uild the mtrix of the liner system of equtions : invert mtrix of size N N 3: do N multiplitions of the inverse mtrix with vetors of size N In Tle I the men CPU time neessry to lulte the voltge sensitivity oeffiients is presented for the IEEE 3 nd 34 node test feeders respetively, when itertions of the method re exeuted It n e oserved tht the nlytil pproh exhiits n improvement of performne whih is of 34 for the IEEE 3 node test feeder nd 5 for the IEEE 34 node test feeder In the sme tle the relevnt 95% onfidene intervls re lso reported for the omputtion of the oeffiients for the two enhmrk feeders One n oserve the dvntge of the proposed nlytil method s the numer of usses in the network inreses It is worth oserving tht suh n improvement depends not only on 5 As lredy pointed out in SeII-A trditionl Join sed sensitivity omputtions do not ount indeed the vritions of tp-hngers the numer of usses ut lso on the network topology (ie sprsity of the [Y] dmittne mtrix) Tle I CPU TIME NECESSARY FOR CALCULATING VOLTAGE SENSITIVITY COEFFICIENTS IN THE IEEE 3 AND THE 34 NODE TEST FEEDERS WHEN ALL PHASES OF ALL BUSSES HAVE LOADS Join Anlytil rtio 3 us feeder 88 ± 8 mse 5 ± 43 mse 34 34 us feeder 98 ± 3 mse 834 ± 59 mse 5 III NUMERICAL VALIDATION The numeril vlidtion of the proposed method for the omputtion of voltge/urrent sensitivities is performed with two different pprohes In prtiulr, s the inverse of the lod flow Join mtrix provides the voltge sensitivities, the omprison reported elow mkes referene to suh method for the voltge sensitivities only On the ontrry, s the inverse of the lod flow Join mtrix does not provide urrent sensitivity oeffiients, their ury is evluted y using numeril pproh where the lod flow prolem is solved y pplying smll injetion perturtions into given network (see Setion II-A) A similr pproh is deployed to vlidte the sensitivities with respet to tp positions of the trnsformers, ie smll perturtions of the voltge mgnitude of one phse of the slk us nd solution of the lod flow prolem Fig shows the IEEE 3 nodes test feeder implemented in the EMTP-RV simultion environment ([3], [4], [5]) dopted to perform the multiphse lod flow BUS646_6 BUS645_5 5kW 5kVAR Lod4 PI3_3 9kW 45kVAR Lod5 BUS6_ BUS684_ PI_3 BUS65_3 48kW 6kVAR Lod3 BUS65_ PI_5 PI7_ 37kW 3kVAR Lod Phse: Slk: 46kVRMSLL/_ BUS63_ PI_3 BUS67_7 PI_ PI5_ PI6_ PI_5 Lod8 86kW 58kVAR 37kW 88kVAR 45kW 4kVAR 45kW 4kVAR Lod9 PI_5 BUS68_ Lod6 BUS633_3 9kW 7kVAR Lod PI_5 Lod5 45kW 4kVAR Lod7BUS69_8 BUS634_4 PI_5 8kW 3kVAR 4kW kvar 4kW kvar BUS675_9 Lod Lod Lod3 6kW kvar 88kW 8kVAR 3kW 3kVAR Lod6 Lod4 Lod7 Figure IEEE 3 node test feeder represented in the EMTP-RV simultion environment
7 Voltge sensitivity of phse of us 8 (pu) Voltge sensitivity of phse α of us 8 (pu) 3 9 8 7 Sensitivity oeffiient x 6 3 () 6 5 6 4 4 6 Ative power genertion/sorption t phse of us 9 (kw) 34 3 3 8 6 4 3 Sensitivity oeffiient 35 4 6 4 4 6 Ative power genertion/sorption t phse of us 9 (kw) 5 5 5 5 5 5 etween Join nd Anlytil method x 7 () etween Join nd Anlytil method For the ske of revity we limit the vlidtion of the proposed method to redued numer of usses exhiiting the lrgest voltge sensitivity ginst P Q lod/injetions In prtiulr, we refer to the vrition of voltges t us 8 with respet to lod/injetion in us 9, ie Ē 8 P 9, Ē 8 P 9, Ē 8 Q, Ē 8 9 Q 9 In Fig() the voltge sensitivity of phse us 8 is shown with respet to tive power sorption nd genertion t phse of us 9 We ssume the onvention tht positive vlues of P nd Q denote power sorption, wheres negtive vlues orrespond to power genertion Fig() shows for the sme usses s Fig(), the sme sensitivity ut referring to voltge nd power elonging to different phses Additionlly, Fig () nd (d) show the voltge sensitivity of us 8 with respet to retive power sorption nd genertion t us 9 In ll these four figures the dshed line represents the reltive error etween the trditionl pproh (ie sed on the inverse of the Join mtrix) nd the nlytil method proposed here As it n e oserved, the overll errors re in the order of mgnitude of 6 In Fig3() nd Fig3() the urrent sensitivity oeffiient of phse of rnh 3 is presented with respet to tive nd retive power sorption/genertion t phse of us 3 In the sme figures, the dshed lines represent the reltive error etween the nlytil vlues nd Voltge sensitivity of phse of us 8 (pu) 8 6 4 Sensitivity oeffiient x 6 () 8 6 4 4 6 Retive power genertion/sorption t phse of us 9 (kvr) 8 6 4 etween Join nd Anlytil method Current sensitivity oeffiient of phse α of rnh 3 (pu) Sensitivity oeffiient x 4 () 6 3 4 4 6 Ative power genertion/sorption t phse α of us 3 (kw) 4 6 8 etween Anlytil nd numeril method Voltge sensitivity of phse α of us 8 (pu) 3 3 34 36 38 4 Sensitivity oeffiient x 7 (d) 4 4 6 4 4 6 6 Retive power genertion/sorption t phse of us 9 (kvr) 4 6 8 etween Join nd Anlytil method Current sensitivity oeffiient of phse α of rnh 3 (pu) 3 Sensitivity oeffiient x 4 () 6 4 4 4 6 Retive power genertion/sorption t phse α of us 3 (kvr) 4 6 8 etween Anlytil nd numeril method Figure Voltge sensitivity oeffiient of phse nd of us 8 with respet to tive nd retive power genertion/sorption t phse of us 9 Figure 3 Current sensitivity oeffiients of phse of rnh -3 with respet to power genertion/sorption t phse of node 3
Voltge sensitivity of phse α of us 7 Voltge sensitivity oeffiients (pu) Voltge sensitivity of phse α of us 7 Voltge sensitivity oeffiients (pu) Voltge sensitivity of of phse α of of us us 7 7 Voltge sensitivity oeffiients (pu) Voltge sensitivity oeffiients (pu) 8 the numeril ones Even for these oeffiients extremely low errors re otined Conerning the vlidtion of voltge sensitivities ginst tp-hnger positions, we hve mde referene to the IEEE 3 node test feeder where the slk us nd therefore the primry susttion trnsformer is pled in orrespondene of node We ssume to vry the slk us voltge of ±6% over 7 tp positions (where position refers to the network rted voltge) In Fig 4 the sensitivity of voltge in phse of us 7 is shown wrt the tp positions in phse, nd of the slk Also, in this se the differene etween the nlytilly inferred sensitivities nd the numeril omputed ones is negligile (ie in the order of mgnitude of 4 ) 5 5 5 95 85 75 65 55 45 5 3 6 4 7 8 9 5 5 35 45 55 Distne from slk (ft) () Voltge sensitivity oeffiients Ē i P3 with respet to tive power sorption t phse of node 3 s funtion of the distne from the slk us 3 5 4 3 9 8E-4 x -4 7E-4 6E-4 5E-4 4E-4 3E-4 E-4 8 Sensitivity oeffiient E-4 7 E -38-34-3-6--8-4- -6-6 4 8 6 3 34 38 Trnsformer's tp tp positions phse phse α () Voltge sensitivity oeffiient of phse of us 7 with respet to trnsformer s tp position t phse of the slk us E-5 x -5 8E-5 6E-5 8 4E-5 E-5 6 E-5 8E-6 8 4 Sensitivity oeffiient 6E-6 6 4 4E-6 E-6 E -38-34-3-6--8-4- -6-6 4 8 6 3 34 38 Trnsformer's tp positions phse () Voltge sensitivity oeffiient of phse of us 7 with respet to trnsformer s tp position t phse of the slk us - 5 5 35 45 55 3 5 6 4-3 -4-5 -6-7 7 Distne from slk (ft) 8 9 () Voltge sensitivity oeffiients Ē i P3 with respet to tive power sorption t phse of node 3 s funtion of the distne from the slk us 9 7 5 3 9 3 4 7 5 6 5 5 5 35 45 55 Distne from slk (ft) () Voltge sensitivity oeffiients Ē i Q with respet 3 to retive power sorption t phse of node 3 s funtion of the distne from the slk us 7 8 9 3 3-5 -3-35 -4 3E-4 x -4 5E-4 E-4-5 5 35 45 55-3 5 6 3 4-4 -45 5E-4-5 E-4-55 Sensitivity oeffiient -6 5 5E-5-38-34-3-6--8-4- -6-6 4 8 6 3 34 38 Trnsformer's tp positions phse -5-6 -7-8 8 7 Distne from slk (ft) 9 3 () Voltge sensitivity oeffiient of phse of us 7 with respet to trnsformer s tp position t phse of the slk us Figure 4 Voltge sensitivity oeffiient of phse of us 7 with respet to trnsformer s tp positions (d) Voltge sensitivity oeffiients Ē i Q with respet 3 to retive power sorption t phse of node 3 s funtion of the distne from the slk us Figure 5 Voltge sensitivity oeffiients with respet to power sorption t phse of us 3 s funtion of the distne from the slk us
9 It is worth oserving tht for the se of the voltge sensitivities, oeffiients tht refer to the voltge vrition s funtion of perturtion (power injetion or tphnger position) of the sme phse, show the lrgest oupling lthough non-negligile ross dependeny n e oserved etween different phses Finlly, Fig5 depits the vrition of voltge sensitivity oeffiients in ll the network with respet to tive nd retive power sorption t phse of us 3 s funtion of the distne from the slk us in feet This type of representtion llows to oserve the overll network ehvior ginst speifi PQ usses sorptions/injetions In prtiulr, we n see tht lrger sensitivities re oserved when the distne etween the onsidered voltge nd the slk us inreses Furthermore, lower, ut quntified dependeny etween oeffiients relted to different phses, n e oserved Also, s expeted, retive power hs lrger influene on voltge vritions lthough the tive power exhiits non negligile influene From the opertionl point of view it is worth oserving tht, figures s Fig5, provide to network opertors n immedite view of the response of the eletril network ginst speifi lods/injetions tht ould lso e used for losed loop ontrol or ontingeny nlysis Lod7 Lod8 Lod6 Lod7 Lod6 PI5 PI6 PI7 Lod3 PI3 Lod4 Lod5 PI4 Lod4 Slk: 49kVRMSLL/_ PI Lod34 Lod33 Lod PI PI PI3 PI5 PI6 PI7 PI8 PI9 PI BUS88_5 BUS84_3 PI4 Lod5 BUS83_6 BUS86_4 BUS85_9 BUS83_ Lod9 BUS854_7 PI8 PI9 Lod3 Lod3 PI Phse: P=3kW Q= DER_8 BUS856_8 BUS888_ BUS89_ BUS8_5 PI3 BUS8_ BUS8_ BUS86_3 BUS88_4 BUS86_9 BUS86_3 Lod35 PI3 BUS834_5 Lod3 BUS836_3 Lod5 Lod9 Lod BUS8_6 PI4 Lod6 BUS84_7 BUS85_8 BUS88_ BUS8_ PI5 PI34 PI3 Lod4 PI6 Lod BUS84_34 PI PI7 Lod DER_3 P=3kW Q= Phse: BUS8_ Lod7 Lod3 BUS858_3 BUS864_4 DER_4 P=6kW PI Q= Phse: PI8 PI3 PI33 BUS86_3 BUS84_6 BUS844_7 Lod8 BUS848_9 BUS846_8 Lod9 PI9 BUS838_33 DER_33 P=3kW Q= Figure 6 IEEE 34 node test feeder represented in the EMTP-RV simultion environment Lod Lod Lod IV APPLICATION OF THE PROPOSED PROCEDURE TO THE PROBLEM OF OPTIMAL VOLTAGE CONTROL For the pplition prt, the IEEE 34 test node feeder is onsidered s depited in Fig6 In usses 8, 3, 4 nd 33 we ssume to hve distriuted energy resoures tht the Distriution Network Opertor (DNO) n ontrol in terms of tive nd retive power Their initil operting vlues, s well s their rted power outputs, re shown in Tle II Furthermore, the DNO hs ontrol on the trnsformer s tp positions Tle II INITIAL AND MAXIMUM OPERATIONAL SET POINTS OF THE DERS AND THE TAP-CHANGERS IN THE 34 TEST NODE FEEDER P init (kw) P mx(kw) n init n min n mx DER 8 3 DER 3 6 36 36 DER 4 5 6 DER 33 5 3 The optiml ontrol prolem is formulted s liner one tking dvntge of the voltge sensitivity oeffiients The ontrolled vriles re the us node voltges nd the ontrol vriles re the tive nd retive power injetions of the DER nd the trnsformer s tp positions under the ontrol of the DNO, x = [ P DER, Q DER, n] It is importnt to stte tht, formlly, this prolem is mixed integer optimiztion prolem due to the tp positions of the trnsformers However, for resons of simpliity, the tp positions re onsidered pseudo-ontinuous vriles whih re rounded to the nerest integer one the optiml solution is rehed The ojetive of the liner optimiztion prolem relevnt to the prolem is: min x Ēi Ē () The linerized reltionship tht links us voltges with ontrol vriles is expressed in the following wy (eg [4]): Ēi = K Pi P i K Qi Q i K ni n i (3) where K Pi is the vetor of sensitivity oeffiients with respet to the tive powers of the DERs, K Qi is the vetor of sensitivity oeffiients with respet to the retive powers of the DERs nd K ni is the vetor of sensitivity oeffiients with respet to the trnsformer s tp positions The imposed onstrints on the opertionl points of the DERs nd the tp positions re the following: P DERi P DERimx (4) Q DERimin Q DERi Q DERimx n min n n mx In order to simplify the nlysis, we hve ssumed tht the DER pility urves re retngulr ones in the PQ plne The formulted linerized prolem is solved y using the liner lest squres method The method used to lulte nlytilly the sensitivity oeffiients llows us to onsider two different optimiztion senrios In the first (opt ), the opertor of the system is ssumed to ontrol the set points of the DERs onsidering tht they re injeting equl powers
Voltge (pu) Voltge (pu) Voltge (pu) into the three phses, wheres in the seond se (opt ) it is ssumed to hve more sophistited ontrol on eh of the phses independently exept for the tp-hngers positions It is worth noting tht this seond option, lthough fr from relisti implementtion, llows us to show the pility of the proposed method to del with the inherent unlned nture of distriution networks Tle III nd Tle IV show the optiml opertionl set points orresponding to these ses Tle III OPTIMAL OPERATIONAL SET POINTS OF THE DERS AND THE TAP-CHANGERS IN THE 34 TEST NODE FEEDER WHEN THE SYSTEM OPERATOR HAS CONTROL ON THEIR 3-PHASE OUTPUT P opt (kw) Q opt (kvr) n opt DER 8 3 3 DER 3 6 6 - DER 4 6 646 DER 33 3-446 Additionlly, in Fig7 the voltge profile of the usses of the system is presented in the initil nd the optiml ses The solid line in the figures shows the initil voltge profile, the solid line with the mrkers shows the first se optiml senrio (opt ) nd the dshed line represents the seond se where the DNO hs full ontrol in eh of the phses of the DERS (opt ) The offset in the grphs, oserved in the slk us, depits the optiml tp position in eh se Wht n e oserved is tht, when there is possiility to ontrol eh of the three phses of the DERs output, the optiml voltge profile is etter thn the one orresponding to ontrol of the 3-phse output of the set points of the DERs Tle IV OPTIMAL OPERATIONAL SET POINTS OF THE DERS AND THE TAP-CHANGERS IN THE 34 TEST NODE FEEDER WHEN THE SYSTEM OPERATOR HAS CONTROL ON EACH OF THE THREE PHASES INDEPENDENTLY P opt (kw) Q opt (kvr) n opt DER 8 DER 8-8856 DER 8 83 DER 3 DER 3 DER 3 DER 4 8 DER 4 965 DER 4 4 DER 33-788 DER 33 DER 33 984 V CONCLUSION In this pper we hve proposed new method for the nlytil omputtion of voltges nd urrents sensitivity oeffiients s funtion of the nodl power injetions The ontriutions of the proposed method re the following: (i) it is generlized to ount for generi numer of slk usses; (ii) it llows the omputtion of sensitivities wrt tp-hnger positions (iii) it is proved to dmit unique solution for the se of rdil networks nd (iv) it supports the omputtion of the sensitivities for generi unlned eletril network y using the [Y] ompound mtrix eing, thus, suitle for distriution systems 5 995 99 985 98 975 97 965 96 955 95 945 94 Phse initil voltge profile three phse lned injetions three phse unlned injetions 5 5 Distne from slk us (ft) () IEEE 34 node test feeder - Voltge profile of phse of the usses 5 5 995 99 985 98 975 97 Phse initil voltge profile three phse lned injetions three phse unlned injetions 5 5 Distne from slk us (ft) () IEEE 34 node test feeder - Voltge profile of phse of the usses 5 995 99 985 98 975 97 965 96 955 95 945 Phse initil voltge profile three phse lned injetions three phse unlned injetions 5 5 Distne from slk us (ft) () IEEE 34 node test feeder - Voltge profile of phse of the usses Figure 7 feeder Initil nd optimized voltge profile of the IEEE 34 node test Compred to the trditionl use of the Join lod-flow mtrix, it llows us to redue the omputtion time y lmost ftor of three, thus enling, in priniple, its implementtion in rel-time optiml ontrollers The pper hs lso vlidted the proposed method y mking referene to typil IEEE 3 nd 34 nodes distriution test feeders The former hs een used to numerilly vlidte the omputtion of the oeffiients whilst the ltter hs een used to show n pplition exmple relted to possile integrtion of the proposed method for the prolem of optiml voltge ontrol in unlned distriution systems It is worth oserving tht the proposed nlytil omputtion of voltges nd urrents sensitivities enles the redution of the omputtionl time of severl trditionl power systems
prolems involving non-negligile omputtionl efforts, suh s rel-time entrlized ontrols, ontingeny nlysis or optiml plnning REFERENCES [4] J Mhseredjin, Simultion des trnsitoires életromgnétiques dns les réseux életriques, Édition Les Tehniques de l Ingénieur, 8 [5] J Mhseredjin, S Dennetière, L Dué, B Khodkhhin, nd L Gérin-Ljoie, On new pproh for the simultion of trnsients in power systems, Eletri power systems reserh, vol 77, no, pp 54 5, 7 [] N Singh, E Kliokys, H Feldmnn, R Kussel, R Chrustowski, nd C Joorowiz, Power system modelling nd nlysis in mixed energy mngement nd distriution mngement system, Power Systems, IEEE Trnstions on, vol 3, no 3, pp 43 49, 998 [] N Jenkins, R Alln, P Crossley, D Kirshen, nd G Str, Emedded genertion, A T Johns nd D F Wrne, Eds IET,, vol 9, no 3 [3] N-G Jmes, Control nd Automtion of Eletril Power Systems Hooken, NJ: CRC Press, 6 [4] A Borghetti, M Bosetti, S Grillo, S Mssuo, C Nui, M Polone, nd F Silvestro, Short-term sheduling nd ontrol of tive distriution systems with high penetrtion of renewle resoures, Systems Journl, IEEE, vol 4, no 3, pp 33 3, [5] Q Zhou nd J Bilek, Genertion urtilment to mnge voltge onstrints in distriution networks, Genertion, Trnsmission & Distriution, IET, vol, no 3, pp 49 498, 7 [6] T Senjyu, Y Miyzto, A Yon, N Urski, nd T Funshi, Optiml distriution voltge ontrol nd oordintion with distriuted genertion, Power Delivery, IEEE Trnstions on, vol 3, no, pp 36 4, 8 [7] S Conti, S Riti, nd G Vglisindi, Voltge sensitivity nlysis in rdil mv distriution networks using onstnt urrent models, in Industril Eletronis (ISIE), IEEE Interntionl Symposium on IEEE,, pp 548 554 [8] D Khtod, V Pnt, nd J Shrm, A novel pproh for sensitivity lultions in the rdil distriution system, Power Delivery, IEEE Trnstions on, vol, no 4, pp 48 57, 6 [9] J Peshon, D Piery, W Tinney, nd O Tveit, Sensitivity in power systems, Power Apprtus nd Systems, IEEE Trnstions on, no 8, pp 687 696, 968 [] D Shirmohmmdi, H Hong, A Semlyen, nd G Luo, A ompenstion-sed power flow method for wekly meshed distriution nd trnsmission networks, Power Systems, IEEE Trnstions on, vol 3, no, pp 753 76, 988 [] A Wood nd B Wollenerg, Power genertion, opertion, nd ontrol Wiley New York, 996, vol [] R Mronto, Eletri power systems, vol, CEI, Itlin Eletrotehnil Committee, Milno, Itly, [3] M Begovi nd A Phdke, Control of voltge stility using sensitivity nlysis, Power Systems, IEEE Trnstions on, vol 7, no, pp 4 3, 99 [4] Q Zhou nd J Bilek, Simplified lultion of voltge nd loss sensitivity ftors in distriution networks, in Pro of the 6th Power Systems Computtion Conferene (PSCC8), Glsgow, Sotlnd, 8 [5] J Bndler nd M El-Kdy, A unified pproh to power system sensitivity nlysis nd plnning, prt i: Fmily of djoint systems, in Pro IEEE Int Symp Ciruits Syst, 98, pp 68 687 [6], A new method for omputerized solution of power flow equtions, Power Apprtus nd Systems, IEEE Trnstions on, no, pp, 98 [7], A unified pproh to power system sensitivity nlysis nd plnning, prt ii: Speil lss of djoint systems, in Pro IEEE Int Symp Ciruits Syst, 98, p 688 [8] L Ferreir, Tellegen s theorem nd power systems-new lod flow equtions, new solution methods, Ciruits nd Systems, IEEE Trnstions on, vol 37, no 4, pp 59 56, 99 [9] R Gurrm nd B Surmnym, Sensitivity nlysis of rdil distriution network-djoint network method, Interntionl Journl of Eletril Power & Energy Systems, vol, no 5, pp 33 36, 999 [] J Arrillg, D Brdley, nd P Bodger, Power system hrmonis, 985 [] H Ching nd M Brn, On the existene nd uniqueness of lod flow solution for rdil distriution power networks, Ciruits nd Systems, IEEE Trnstions on, vol 37, no 3, pp 4 46, 99 [] A Aur nd A Exposito, Power system stte estimtion: theory nd implementtion CRC, 4, vol 4 [3] J Mhseredjin, S Lefevre, nd X Do, A new method for time domin modelling of nonliner iruits in lrge liner networks, in Pro of th Power Systems Computtion Conferene PSCC, 993