复杂电力系统电气故障电磁暂态数字计算方法研究

2712Vol.27No.12200312PowerSystemTechnologyDec.20031000-3673200312-0031-06TM731A1,2,216500512150001ANEWELECTROMAGNETICTRANSIENTSSIMULATIONMETHODFORELECTRICFAULTSINCOMPLEXPOWERSYSTEMSHUHong-chun1,SIDa-jun2,CHENXue-yun21.KunmingUniversityofScienceandTechnology,Kunming650051,YunnanProvince,China;2.HarbinInstituteofTechnology,Harbin150001,HeilongjiangProvince,ChinaABSTRACT:Thestudyondigitalsimulationofelectro-magnetictransientsisaneverlastingtopicintheresearchofpowersystem.Especially,afterimprovedbyDommel,thetime-domainBergeronmodelissuccessfullyappliedinEMTP,butthepre-processingbeforesimulatingtransientscausedbyvariousfaultsandoperationsisquitetroublesome.Anewelectromagnetictransientssimulationmethodwhichcanbeversatilelyusedtocalculatemostoftheelectricfaultsinpowersystemispresented.UnlikeEMTP,itisnotnecessaryforthepresentedmethodtocalculatetheinitialvaluesofthesystemagainwhenthestructureorparameterofthesystemischanged.Meanwhile,thefaultsoccurredatanylocationofsingle-anddouble-circuitlinesincludingseriescompensatedlines,busfaults,operationofbreakersandopen-circuitfaultscanbesimulatedbythismethod.Inaddition,thefault-start-anglecanbesetinthemethod.Becausethenewmethodcanconvenientlysimulatenotonlyinrushcurrents,overvoltagesandharmonicscomponentsatdifferentfaultlocations,butalsodevelopingfaults,itissignificantinthesimulationoflineprotectionverificationandovervoltagecomputation.Theresultsfromlotsoftestsshowthatthemethodisaccurateandfast.KEYWORDS:Electromagnetictransients;Digitalsimulation;Complexpowersystem;ElectricfaultsBergeronDommelEMTPEMTPEMTPEMTP──1EMTP[1]EMTP[2]2000B20298E0409M99E006G2002E0025M32PowerSystemTechnologyVol.27No.12[3,4][5,6]22.11(a)lRxLxCx1(b)(Bergeron)[7]KNiKN(t)iNK(t)uK(t)uN(t)(a)ZKNiKN(t)iNK(t)IK(t-t)IN(t-t)uK(t)uN(t)(b)1Fig.1TheequivalentcircuitandsinglephasetransmissionlinetZc(1)⎪⎩⎪⎨⎧+=⋅==4//cxxxcRZZRlRCLZ(1)IK(t-t)IN(t-t)(2)[]}[]{[]}⎪⎪⎩⎪⎪⎨⎧-+-⋅-+-+-⋅+-=--+-⋅-+-+-⋅+-=-)()()1()()()1(5.0)()()()1()]()()[1{(5.0)(ttttttttttthituYhthituYhtIthituYhthituYhtINKNKNKNKNKNKNK(2)(2)⎩⎨⎧+-=+=)25.0/()25.0()25.0/(1cccRZRZhRZY(3)(2)iKN(t)iNK(t)1⎩⎨⎧-+⋅=-+⋅=)()()()()()(tttItuYtitItuYtiNNNKKKKN(4)(4)(2)⎪⎪⎪⎩⎪⎪⎪⎨⎧-++----++--⋅+-=--++----++--⋅+-=-(5))]2()1()2()1[(0.5)]()1()()1[()1(5.0)()]2()1()2()1[(0.5)]()1()()1[()1(5.0)(tttttttttttIhtIhhtuhtuhYhtItIhtIhhtuhtuhYhtIKNKNNNKNKK(5)TKarenbauerTT-1⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=211121111T⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=-101011111311T2.22(a)2(b)2RLLIL(t-Δt)RL=tLΔ2(6)271233)]()([1)()(ttuttuRttittImkLkmLΔ--Δ-+Δ-=Δ-(7)ikm(t)kmLuk(t)um(t)uL(t)(a)uk(t)um(t)kmikm(t)IL(t-Δt)RL(b)2Fig.2Theequivalentcircuitofinductanceforcalculation(7)IL(t-tΔ)ikm(t-tΔ)uk(t-tΔ)um(t-tΔ)2(b)(8))()]()([1)(ttItutuRtiLmkLkmΔ-+-=(8)(6)(8))]()([2)2()(ttuttuRttIttImkLLLΔ--Δ-+Δ-=Δ-(9)2.33kmCuk(t)um(t)uC(t)ikm(t)(a)uk(t)um(t)kmikm(t)IC(t-Δt)RC(b)3Fig.3Theequivalentcircuitofcapacitorforcalculation3RCIC(t-tΔ)CRC=Ct2Δ(10))]()([1)()(ttuttuRttittImkCkmCΔ--Δ--Δ--=Δ-(11))]()([2)2()(ttuttuRttIttImkCCCΔ--Δ--Δ--=Δ-(12)2.44[8,9]MOVMOVMOVMOVMOVImax(Imax)XcIMOVIMOV4Fig.4Theseriescompensatedcapacitorprotectionscheme4XcIIMOVMOVInMOVIpl=k×In(k=2~3)2~3IplMOVVplVpl=2×Xc×IplMOV()a1PmaxMOV/VVII=(13)Va30~502.534PowerSystemTechnologyVol.27No.123fU*∠qfU*∠-=ΔqqqΔEMTP4[3]C++C++1234567855(a)(b)(c)5Fig.5Partialinterfacesofthesimulationprogram27123556C=100μFVpl=340kVa=30°S1S2S3S41l1l2l3l4l5l62S1S2S4l1l2l6CMS3l5l4l3NQR6Fig.6Thesystemusedforsimulation1Tab.1Parametersofsystems/kV/()/H/H/Ω/ΩS133500.198640.168649.1396.139S2265-300.20730.1467712.208.56S3305-200.14330.123648.427.32S4286-250.33740.2537623.6420.482Tab.2ParametersoflinesΩ/kmΩ/kmmH/kmmH/kmpF/kmpF/kmkml10.2156580.0296062.154970.85892157610.514967.6150l22156580296062.1549785892157610.514967.6150l30.2233640.0243351.963550.9890518131.513122100l40.2233640.0243351.963550.9890518131.513122120l50.2233640.0243351.963550.9890518131.513122160l60.2233640.0243351.963550.9890518131.513122135Cl20.02s1C33Nl2789Al2NAN60kml40.02sA090Nl4101112l3MABC0.01s0.015s0.02sl3M7l2NFig.7Zero-sequencevoltageandwaveformonlinel2atN-terminal8AFig.8A-Phasevoltagewaveformonseriescompensationcapacitor9l2NAFig.9A-Phasecurrentwaveformonlinel4atN-terminal100ANl4Fig.10Voltageandcurrentwaveformsonlinel4atN-terminalafterA-Gfaultwhenfaultstartangleis036PowerSystemTechnologyVol.27No.121190ANl4Fig.11Voltageandcurrentwaveformsonlinel4atN-terminalafterA-Gfaultwhenfaultstartangleis9012l3MFig.12Voltageandcurrentwaveformsonlinel3atM-terminalcausedbyclosingbreakers612EMTPEMTP3——[1]MichelisInitializationandSteady-StateSolution[Z]EMTPSummerCourseLEUVEN1987[2]JohnsATAggarwalRKetalDigitalsimulationoffaultedEHVtransmissionlineswithparticularreferencetovery-high-speedprotection[J]IEEProc1976123(4)353-359[3]ShuHongchunSiDajunGaoFengetalAnewnumericalmethodforcomputingelectromagnetictransientsofHVtransmissionlinefaultsandcircuitbreakeroperations[J]AutomationofElectricPowerSystems200024(20)23-2648[4]ShuHongchunSiDajunGeYaozhongAstudyofelectromagnetictransientdigitalsimulationondevelopingfaultsonHVtransmissionlines[J]TransactionsofChinaElectro-technicalSociety200015(3)70-74[5]Wan