多风电场并网时安全约束机组组合的混合整数规划解法

,,,(,510640):,,,,2IEEE118,:;;;;:2009209215;:2009212203(NCET20820207);(109128)0,,[122],,,(SCUC)(SCED)[325]:,,;,[6210],,[7]SCUC,,,[8]SCUC[9],[10],2:,,,;(MIP),()MIP[11212],:;,,,;,,MIP,,MIP[13],[14]MIP[13]MIP[5]SCUCMIP;,,2,MIP1MIPSCUC1.1(GGDF),:P=Bj-jxjLjPjJ(1)923452010310Vol.34No.5Mar.10,2010:P=B,j;xjj;Ljj;J[14]GGDF,,(1),3MIPGGDFSCUC,,21.2,:mingGTtNTF(g,t)+Su(g,t)(2):GT;NT;FSu1:s.t.F(g,t)=(agp2(g,t)+bgp(g,t)+cg)I(g,t)gGT;tNT(3)Su(g,t)K(g,mg)I(g,t)-mgh=1I(g,t-h)Su(g,t)0;gGT;tNT;Pmg(4)pgI(g,t)p(g,t)+r(g,t)pgI(g,t)gGT;tNT(5)gGTr(g,t)R(t)tNT(6)p(g,t)-Rd(g)p(g,t+1)p(g,t)+Ru(g)gGT;t=1,2,,NT-1(7)Tu(g)-1h=0I(g,t+h)(I(g,t)-I(g,t-1))min(Tu(g),NT-t+1)gGT;tNT(8)Td(g)-1h=0(1-I(g,t+h))(I(g,t-1)-I(g,t))min(Td(g),NT-t+1)gGT;tNT(9)gGTp(g,t)+gGWp(g,t)=D(t)tNT(10)gGTgjp(g,t)+gGWgjp(g,t)LjjJ;tNT(11):I(g,t)tg;p(g,t)tg;K(g,mg)gmg;(5);pgpg;(6);r(g,t)gt,R(t)t;(7);Rd(g)Ru(g)g;(8)(9);Tu(g)Td(g)g;(10);GW;D(t)t;(11);gjgj2:pgI(g,t)p(s,g,t)pgI(g,t)PtNT;PsSt;gGT(12)gGTp(s,g,t)+gGWp(s,g,t)=D(t)(13)gGTgjp(s,g,t)+gGWgjp(s,g,t)LjjJ(14)-16Rd(g)p(s,g,t)-p(g,t)16Ru(g)gGT(15):Stt;p(s,g,t)tgs;(12)s,I(g,t);(13)s;(14)s;(15),p(s,g,t)-p(g,t),g,1/6,10min:3,[16]I(g,t),MIP,1.31,X0t,si,,032010,34(5)siXi,si,tSts1Fig.1Feasiblesolutionsundertheconstraintsofadjustmenttime,,,2St,,,,,2Fig.2Curveoftheforecastedwindpoweroutputanditsprobabilityintervalcomparedwithrealvalues2St2,t110MW120MW,(40MW,180MW)(60MW,180MW)3,tSt,,,,[10],St,,n(n3),2nn,n,St,32Fig.3Spaceofwindgenerationconsideringtwowindfarms,,St,,,(12)(15)p(s,g,t)(gGTGW),I(g,t)p(g,t),(12)(15):Ax=bx0(16)x=[,xGW,]T,xGWp(s,g,t)(gGW),xxxGWU={x1,x2,,x2n},1,2,,2n,i[0,1]ki=1i=1,2ni=1ixiU,UU2nx,U,{xGW1,xGW2,,xGW2n},U,(16)[17],UA,UAn,,,,,13,3IEEE118,5418691,3733MW12h14h15h,motor.ece.iit.edu/data/SCUC_118test.xls[10],124h,CPLEX10.11Table1Hourlyforecastedgenerationofwindpower/MW26597377921207.270.062.1100100.02204.672.161.2103102.53253.377.766.0111104.04257.160.966.0113107.55212.881.262.411695.06164.188.255.512690.07215.297.362.113992.58140.5109.251.016097.59126.6120.448.018085.010115.5128.145.018379.011173.4134.457.019280.012147.1138.651.019890.013170.3138.659.719897.514269.2130.266.0170105.015185.2105.060.015095.016181.5116.251.016687.517131.3138.645.019880.018204.3138.643.519875.019122.4138.642.019570.020133.5138.645.019275.021105.6138.640.519677.522100.8138.642.019880.023126.7130.245.018685.024146.1113.448.016280.01,26,59,9230%,73,7750%,,,31:(1)St,4,10,102:GGDF,10000,10St3:GGDFSt,4,10,1059,26,92,77,731,,443Fig.4Performancecurvesofthreemethods4(a)1,3,4,1310,,210,,,4,34(b)34(a),23GGDF10000,24h240000,,(15)4(c),,13100%,53210,2,,4232010,34(5)4(d),3,,4,2,4(c)5,[18]5,4h28.9%3253.8s,730747,96.32%,AA155Fig.5Penetrationlevelsofwindpowerwith5windfarmsineachhour,:MIP;;31185,20%4GGDF,MIP,2MIP,,,()[1]..,2003,27(8):84289.LEIYazhou.Studiesonwindfarmintegrationintopowersystem.AutomationofElectricPowerSystems,2003,27(8):84289.[2],LIGHTBODYG..,2005,29(10):125.LEIYazhou,LIGHTBODYG.Windenergyandelectricitymarket.AutomationofElectricPowerSystems,2005,29(10):125.[3]KUSIAKAndrew,ZHENGHaiyang,SONGZhe.Short2termpredictionofwindfarmpower:adataminingapproach.IEEETransonPowerSystems,2009,24(1):1252136.[4],,..,2008,32(10):84288.QIAOJiageng,XUFei,LUZongxiang.Optimizationanalysismodelofgrid2connectedwindcapacitybasedondependentchanceprogramming.AutomationofElectricPowerSystems,2008,32(10):84288.[5],,..,2006,30(2):22226.CHENHaiyan,CHENJinfu,DUANXianzhong.Fuzzymodelingandoptimizationalgorithmondynamiceconomicdispatchinwindpowerintegratedsystem.AutomationofElectricPowerSystems,2006,30(2):22226.[6]TUOHYA,MEIBOMP,DENNYE,etal.Unitcommitmentforsystemswithsignificantwindpenetration.IEEETransonPowerSystems,2009,24(2):5922601.[7]WULei,SHAHIDEHPOURM,LITao.Stochasticsecurity2constrainedunitcommitment.IEEETransonPowerSystems,2007,22(2):8002811.[8]RUIZPA,PHILBRICKCR,ZAKE,etal.Uncertaintymanagementintheunitcommitmentproblem.IEEETransonPowerSystems,2009,24(2):6422651.[9]BOUFFARDF,GALIANAFD.Stochasticsecurityforoperationsplanningwithsignificantwindpowergeneration.IEEETransonPowerSystems,2008,23(2):3062316.[10]WANGJH,SHAHIDEHPOURM,LIZY.Security2constrainedunitcommitmentwithvolatilewindpowergeneration.IEEETransonPowerSystems,2008,23(3):131921327.[11]LITao,SHAHIDEHPOURM.Price2basedunitcommitment:acaseofLagrangianrelaxationversusmixedintegerprogramming.IEEETransonPowerSystems,2005,20(4):201522025.[12]STREIFFERTD,PHILBRIEKR,OTTA.Amixedintegerprogrammingsolutionformarketclearingandreliabilityanalysis//ProceedingsofIEEEPowerEngineeringSocietyGeneralMeeting,June12216,2005,SanFrancisco,CA,USA:272422731.[13]BIXBYRE,FENOLONM,GUZonghao,etal.MIP:theoryandpracticelosingthegap//POWELLMJD,SCHOLTESS.Systemmodellingandoptimization:methods,theoryandapplications.Dordrecht,TheNetherlands:Kluwer,2000:19249.[14]DUPACOV#J,GR;WE2KUSKAN,R;MISCHW.Scenarioreductioninstochasticprogramming:anapproachusingprobabilitymetrics.MathematicalProgramming,2003,95(3):4932511.[15]NGWY.Generalizedgenerationdistributionfactorsforpowersystemsecurityevaluations.IEEETr