劉麗影,李楊,侯淑蓮,錢序
(中國(guó)電力工程顧問(wèn)集團(tuán)東北電力設(shè)計(jì)院,長(zhǎng)春市130021)
國(guó)外核電站設(shè)計(jì)中一般都要求提供廠用電系統(tǒng)暫態(tài)穩(wěn)定分析報(bào)告,在一些國(guó)外設(shè)計(jì)項(xiàng)目的合同中也明確要求提供廠用電源切換過(guò)程的暫態(tài)分析報(bào)告。我國(guó)對(duì)于廠用電源切換的暫態(tài)穩(wěn)定判據(jù)是校驗(yàn)廠用電源切換中的最低母線電壓,在設(shè)計(jì)中不需要提供暫態(tài)穩(wěn)定分析報(bào)告,顯然跟不上國(guó)外發(fā)展步伐,因此對(duì)發(fā)電廠(尤其是核電站)進(jìn)行暫態(tài)分析研究勢(shì)在必行。
近年來(lái)國(guó)內(nèi)外學(xué)者相繼對(duì)暫態(tài)分析進(jìn)行了研究,文獻(xiàn)[1]指出機(jī)電暫態(tài)仿真是分析電力系統(tǒng)的有效工具,在系統(tǒng)設(shè)計(jì)、規(guī)劃、運(yùn)行和調(diào)度等方面有廣泛的應(yīng)用;文獻(xiàn)[2]介紹了采用天津大學(xué)研發(fā)的動(dòng)態(tài)軟件包(BUSTP)對(duì)內(nèi)蒙古達(dá)拉特發(fā)電廠1號(hào)機(jī)組的廠用電切換過(guò)程進(jìn)行仿真研究;文獻(xiàn)[3]采用PSCAD/EMTDC軟件,以某發(fā)電廠為例,分析了感應(yīng)電動(dòng)機(jī)參數(shù)對(duì)廠用電切換方式的影響;文獻(xiàn)[4]應(yīng)用DL/T 5153《火力發(fā)電廠廠用電設(shè)計(jì)技術(shù)規(guī)定》的計(jì)算公式,研究了二灘水力發(fā)電廠大功率異步電動(dòng)機(jī)起動(dòng)對(duì)廠用電系統(tǒng)電壓的影響;文獻(xiàn)[5]依托某國(guó)外發(fā)電廠(2×660 MW 空冷機(jī)組),利用 ETAP(electrical transient analysis program)軟件對(duì)發(fā)電機(jī)系統(tǒng)及單元、公用高壓廠用電系統(tǒng)、低壓廠用電系統(tǒng)以及升壓站400 kV系統(tǒng)進(jìn)行了仿真計(jì)算。國(guó)外學(xué)者除了對(duì)廠用電切換進(jìn)行研究外,還將研究成果運(yùn)用到設(shè)計(jì)工作中,暫態(tài)穩(wěn)定分析技術(shù)已經(jīng)成熟。
隨著計(jì)算機(jī)仿真技術(shù)的發(fā)展,越來(lái)越多的軟件可以實(shí)現(xiàn)系統(tǒng)暫態(tài)計(jì)算及仿真。在電氣設(shè)計(jì)中,由美國(guó)OTI公司發(fā)行的ETAP軟件具有潮流分析、短路分析、諧波分析、暫態(tài)穩(wěn)定分析、可靠性分析等功能,其界面友好,在國(guó)際上被通用。
田灣核電站1、2號(hào)機(jī)組由俄羅斯設(shè)計(jì)院(簡(jiǎn)稱俄方)設(shè)計(jì),俄方提供了1套完整的廠用電系統(tǒng)的暫態(tài)穩(wěn)定分析報(bào)告,文獻(xiàn)[6]中給出了廠用電源切換過(guò)程中的暫態(tài)仿真曲線和仿真數(shù)據(jù)。本文基于ETAP軟件,以田灣核電站1、2號(hào)機(jī)組的暫態(tài)穩(wěn)定分析報(bào)告為參考,對(duì)田灣核電站3、4號(hào)機(jī)組(2×1 000 MW)進(jìn)行暫態(tài)穩(wěn)定分析研究。
本研究采用ETAP軟件進(jìn)行計(jì)算,ETAP軟件中有潮流分析、短路分析、諧波分析、暫態(tài)穩(wěn)定分析、可靠性分析等功能模塊。ETAP暫態(tài)穩(wěn)定分析模塊旨在分析系統(tǒng)變化(或干擾)時(shí)和變化前、后電力系統(tǒng)的動(dòng)態(tài)響應(yīng)和穩(wěn)定極限。程序模擬電力系統(tǒng)的動(dòng)態(tài)特性,執(zhí)行設(shè)定的事件和動(dòng)作,求解系統(tǒng)網(wǎng)絡(luò)方程和機(jī)械微分方程,仿真系統(tǒng)和電機(jī)在時(shí)域內(nèi)的響應(yīng)。通過(guò)這些響應(yīng),可以確定系統(tǒng)的暫態(tài)特性,做出穩(wěn)態(tài)評(píng)估,及時(shí)調(diào)整繼電保護(hù)整定值,并運(yùn)用必要的補(bǔ)償或加強(qiáng)措施來(lái)提高系統(tǒng)的穩(wěn)定性。ETAP暫態(tài)穩(wěn)定性分析模塊遵從IEEE 1110、IEEE 421.5標(biāo)準(zhǔn)。
田灣核電站規(guī)劃建設(shè)8臺(tái)核電機(jī)組,一期工程2臺(tái)俄羅斯產(chǎn)的1 000 MW WWER1000型壓水堆核電機(jī)組(1、2號(hào)機(jī)組),已于2007年建成投產(chǎn);二期工程計(jì)劃再建2臺(tái)俄羅斯產(chǎn)同型(改進(jìn)型)壓水堆機(jī)組(3、4號(hào)機(jī)組)。二期工程總設(shè)計(jì)方是俄羅斯設(shè)計(jì)院,核島部分由俄方設(shè)計(jì),常規(guī)島部分由東北電力設(shè)計(jì)院設(shè)計(jì)。目前,該工程正進(jìn)行施工圖設(shè)計(jì),在初步設(shè)計(jì)中,俄方要求中方提供廠用電系統(tǒng)的暫態(tài)穩(wěn)定分析報(bào)告,以便開(kāi)展核島部分設(shè)計(jì)。
二期工程2臺(tái)機(jī)組均采用發(fā)變組單元接線,發(fā)電機(jī)出口設(shè)置斷路器,發(fā)電機(jī)與主變壓器通過(guò)離相封閉母線連接。根據(jù)接入系統(tǒng)設(shè)計(jì),二期工程2臺(tái)機(jī)組接入一期工程已建成的500 kV高壓配電裝置,二期工程新增1回500 kV出線。
依托俄羅斯產(chǎn)WWER1000(改進(jìn)型)1 000 MW級(jí)壓水堆核電機(jī)組工藝系統(tǒng)特點(diǎn)及俄方相關(guān)設(shè)計(jì)規(guī)范,根據(jù)現(xiàn)階段廠用負(fù)荷計(jì)算結(jié)果,廠用電源配置的主要原則如下:
(1)每臺(tái)機(jī)組設(shè)2臺(tái)容量為63/31.5-31.5 MVA的分裂繞組高壓廠用工作變壓器30(40)BBT01、30(40)BBT02,其電源由主變壓器低壓側(cè)和發(fā)電機(jī)出口斷路器之間引接。2臺(tái)高壓廠用變壓器均裝設(shè)有載調(diào)壓開(kāi)關(guān),以保證機(jī)組正常起停時(shí)高壓母線電壓水平的要求。
(2)每臺(tái)機(jī)組設(shè)1臺(tái)容量為70/35-35 MVA分裂繞組、帶有載調(diào)壓開(kāi)關(guān)的備用變壓器,其電源引自220 kV輔助開(kāi)關(guān)站二期工程擴(kuò)建的220 kV母線。
(3)廠用電系統(tǒng)采用6.3 kV和0.4 kV兩級(jí)電壓供電,6.3 kV系統(tǒng)采用單母線接線。根據(jù)核電站廠用負(fù)荷對(duì)供電可靠性的要求,以及俄方相關(guān)設(shè)計(jì)標(biāo)準(zhǔn),廠用電設(shè)計(jì)分為3個(gè)系統(tǒng):正常運(yùn)行供電系統(tǒng)、正常運(yùn)行可靠供電系統(tǒng)、應(yīng)急供電系統(tǒng)。圖1為3號(hào)機(jī)組廠用電原則性接線圖(4號(hào)機(jī)組與此圖相同)。
(1)發(fā)電機(jī):額定容量為1 318 MVA,額定功率為1 186.2 MW,功率因數(shù)為0.9,額定電壓為24 kV;直軸超瞬變電抗Xd″(飽和)容差為28.7%(±10%),發(fā)電機(jī)轉(zhuǎn)動(dòng)慣量為326 000 kg·m2。
(2)主變壓器(下稱主變):額定容量為3×430 MVA,電壓比為,阻抗電壓Ud=18%,接線組別為YN、d11。
(3)高壓廠用變壓器(下稱高廠變):額定容量為63/31.5-31.5 MVA,電壓比為24±8×1.5%/6.3-6.3 kV,阻抗電壓Ud=26.5%,容差為+7.5%,接線組別為 D、d0-d0。
(4)高壓廠用備用變壓器(下稱高備變):額定容量為70/35-35 MVA,電壓比為230±8×1.5%/6.3-6.3 kV,阻抗電壓Ud=22.5%,容差為+7.5%,接線組別為 YN、d11-d11。
(5)廠用電系統(tǒng)電動(dòng)機(jī):參照高廠變運(yùn)行負(fù)荷清單和高備變運(yùn)行負(fù)荷清單以及電動(dòng)機(jī)參數(shù)[7]。
圖1 3號(hào)機(jī)組廠用電系統(tǒng)原則性接線圖Fig.1 Connection scheme for auxiliary power system in No.3 unit
本文將校驗(yàn)各種廠用電源切換工況,包括機(jī)組正常運(yùn)行工況,即廠用電系統(tǒng)由發(fā)電機(jī)組通過(guò)高廠變提供電源工況和事故運(yùn)行工況,即廠用電源切換至高備變供電(由廠外220 kV備用電源系統(tǒng)供電)工況。
在廠用電源切換校驗(yàn)中,都是以高廠變(高備變)帶最大負(fù)荷為條件進(jìn)行電源切換計(jì)算的。根據(jù)電源的故障類型及廠用負(fù)荷用電的要求,將廠用電源切換分為以下4種工況。
工況1:單元機(jī)組主變或高廠變故障,主保護(hù)動(dòng)作,同時(shí)跳開(kāi)500 kV斷路器、發(fā)電機(jī)出口斷路器和正常運(yùn)行供電系統(tǒng)6.3 kV母線段工作電源進(jìn)線斷路器。供電電源自動(dòng)切換至高備變,正常運(yùn)行供電系統(tǒng)6.3 kV母線段斷電時(shí)間約0.2 s(保護(hù)動(dòng)作時(shí)間加上斷路器動(dòng)作時(shí)間)。
工況2:單元機(jī)組主變或高廠變故障,主保護(hù)拒動(dòng),后備保護(hù)動(dòng)作,同時(shí)跳開(kāi)500 kV斷路器、發(fā)電機(jī)出口斷路器和正常運(yùn)行供電系統(tǒng)6.3 kV母線段工作電源進(jìn)線斷路器。供電電源自動(dòng)切換至高備變,正常運(yùn)行供電系統(tǒng)6.3 kV母線段斷電時(shí)間約0.7 s(保護(hù)動(dòng)作時(shí)間0.6 s加上斷路器動(dòng)作時(shí)間0.1 s)。
工況3:單元機(jī)組高廠變低壓繞組故障,主保護(hù)拒動(dòng),后備保護(hù)(如高廠變后備距離保護(hù))動(dòng)作,同時(shí)跳開(kāi)500 kV斷路器、發(fā)電機(jī)出口斷路器和正常運(yùn)行供電系統(tǒng)6.3 kV母線段工作電源進(jìn)線斷路器。供電電源自動(dòng)切換至高備變。此時(shí)故障變壓器供電的2段正常運(yùn)行供電系統(tǒng)6.3 kV母線段斷電時(shí)間約0.7 s(保護(hù)動(dòng)作時(shí)間0.6 s加上斷路器動(dòng)作時(shí)間0.1 s);非故障變壓器供電的2段正常運(yùn)行供電系統(tǒng)6.3 kV母線段斷電時(shí)間約0.2 s。
工況4:500 kV系統(tǒng)短路故障,主保護(hù)或斷路器拒動(dòng)后,斷路器失靈保護(hù)又拒動(dòng)(或操縱員誤操作,500 kV系統(tǒng)失電)。正常運(yùn)行供電系統(tǒng)6.3 kV母線段電壓降低,該段工作電源進(jìn)線斷路器要等1.5 s后依靠低電壓保護(hù)動(dòng)作跳開(kāi),供電電源自動(dòng)切換至高備變。此時(shí)4段正常運(yùn)行供電系統(tǒng)6.3 kV母線段斷電時(shí)間約1.5 s。
俄方對(duì)廠用電系統(tǒng)暫態(tài)穩(wěn)定的要求為:在電源切換過(guò)程中,各段高、低壓母線電壓恢復(fù)至65%額定電壓的時(shí)間均小于5 s,各段高、低壓母線電壓恢復(fù)至80%額定電壓的時(shí)間均小于15 s。
應(yīng)用ETAP軟件建立模型,輸入數(shù)據(jù)包括同步發(fā)電機(jī)、主變、等效電網(wǎng)系統(tǒng)、高廠變、高備變、感應(yīng)電動(dòng)機(jī)等設(shè)備參數(shù)以及等效負(fù)荷、靜態(tài)負(fù)荷等。建立模型后應(yīng)用ETAP暫態(tài)穩(wěn)定分析模塊,對(duì)俄方提出的各種工況進(jìn)行仿真,結(jié)果如圖2所示。
仿真工況1:單元機(jī)組主變或高廠變故障,主保護(hù)動(dòng)作,機(jī)組停機(jī);正常運(yùn)行供電系統(tǒng)6.3 kV母線段斷電時(shí)間約為0.2 s;廠用電源切換至高備變,維持反應(yīng)堆40%功率所需的廠用負(fù)荷參加自起動(dòng)。圖2(a)為工況1時(shí),6.3 kV廠用母線電壓曲線。
圖2 不同電源切換工況時(shí)高、低壓母線電壓Fig.2 Voltage of high-voltage and low-voltage bus under power supply switch
仿真工況2:單元機(jī)組主變或高廠變故障,主保護(hù)拒動(dòng),后備保護(hù)動(dòng)作,機(jī)組停機(jī);正常運(yùn)行供電系統(tǒng)6.3 kV母線段斷電時(shí)間約為0.7 s;廠用電源切換至高備變,維持反應(yīng)堆40%功率所需的廠用負(fù)荷參加自起動(dòng)。圖2(b)為工況2時(shí),6.3 kV廠用母線電壓曲線。
仿真工況3A:單元機(jī)組高廠變低壓繞組故障,主保護(hù)拒動(dòng),后備保護(hù)動(dòng)作,機(jī)組停機(jī);故障變壓器供電的2段正常運(yùn)行供電系統(tǒng)6.3 kV母線段BBA和BBB斷電時(shí)間為0.7 s,非故障變壓器供電的2段正常運(yùn)行供電系統(tǒng)6.3 kV母線段BBC和BBD斷電時(shí)間為0.2 s;廠用電源切換至高備變,維持反應(yīng)堆40%功率所需的廠用負(fù)荷參加自起動(dòng)。圖2(c)為工況3A時(shí),6.3 kV廠用母線電壓曲線。
仿真工況3B:單元機(jī)組高廠變低壓繞組故障,主保護(hù)拒動(dòng),后備保護(hù)動(dòng)作,機(jī)組停機(jī);故障變壓器供電的2段正常運(yùn)行供電系統(tǒng)6.3 kV母線段BBC和BBD斷電時(shí)間為0.7 s,非故障變壓器供電的2段正常運(yùn)行供電系統(tǒng)6.3 kV母線段BBA和BBB斷電時(shí)間為0.2 s;廠用電源切換至高備變,維持反應(yīng)堆40%功率所需的廠用負(fù)荷參加自起動(dòng)。圖2(d)為工況3B時(shí),6.3 kV廠用母線電壓曲線。
仿真工況4:500 kV系統(tǒng)短路故障,主保護(hù)或斷路器拒動(dòng)后,斷路器失靈保護(hù)又拒動(dòng);或操縱員誤操作時(shí),500 kV系統(tǒng)失電;廠用系統(tǒng)電壓降低,低電壓保護(hù)動(dòng)作甩負(fù)荷,跳開(kāi)正常運(yùn)行供電系統(tǒng)6.3 kV母線段進(jìn)線斷路器,正常運(yùn)行供電系統(tǒng)6.3 kV母線段斷電時(shí)間約為1.5 s;廠用電源切換至高備變,維持反應(yīng)堆40%功率所需的廠用負(fù)荷參加自起動(dòng)。圖2(e)為工況4時(shí),6.3 kV廠用母線電壓曲線。
圖2(f)為工況1時(shí),0.4 kV應(yīng)急廠用母線電壓曲線。
圖2(g)為工況2時(shí),0.4 kV應(yīng)急廠用母線電壓曲線。
圖2(h)為工況3A時(shí),0.4 kV應(yīng)急廠用母線電壓曲線。
圖2(i)為工況3B時(shí),0.4 kV應(yīng)急廠用母線電壓曲線。
圖2(j)為工況4時(shí),0.4 kV應(yīng)急廠用母線電壓曲線。
國(guó)內(nèi)工程是依據(jù)DL/T 5153—2002《火力發(fā)電廠廠用電設(shè)計(jì)技術(shù)規(guī)定》[8](下稱“規(guī)定”)來(lái)校驗(yàn)成組電動(dòng)機(jī)自起動(dòng)時(shí)母線電壓是否滿足要求的。在“規(guī)定”的第5.4、5.5、5.6條中給出了計(jì)算方法,并提供了保證自起動(dòng)成功的最低母線電壓要求。
本文將ETAP軟件計(jì)算結(jié)果與應(yīng)用規(guī)定計(jì)算方法的計(jì)算結(jié)果進(jìn)行對(duì)比,對(duì)母線電壓進(jìn)行校驗(yàn)。成組電動(dòng)機(jī)自起動(dòng)時(shí)廠用母線電壓校驗(yàn)結(jié)果如表1所示。
表1 成組電動(dòng)機(jī)自起動(dòng)時(shí)廠用母線電壓校驗(yàn)Tab.1 Voltage checking of auxiliary buses during grouped motors self-starting
從表1可得,俄方提出了4種切換工況,工況1、2、3切換時(shí)間均小于0.8 s,屬于快速切換,工況4切換時(shí)間超過(guò)0.8 s,屬于慢速切換。本工程采用“規(guī)定”計(jì)算方法進(jìn)行了母線電壓校驗(yàn)計(jì)算,慢速切換中高壓母線電壓不滿足“規(guī)定”要求,成組電動(dòng)機(jī)自起動(dòng)不成功;快速切換中高、低壓母線電壓都滿足“規(guī)定”要求,成組電動(dòng)機(jī)自起動(dòng)成功[9-10]。
采用ETAP軟件對(duì)4種暫態(tài)工況進(jìn)行仿真,結(jié)果表明可滿足俄方對(duì)廠用電源的要求,即高、低壓各段母線電壓恢復(fù)至65%額定電壓的時(shí)間均小于5 s,高、低壓各段母線電壓恢復(fù)至80%額定電壓的時(shí)間均小于15 s。
采用ETAP軟件仿真,工況1中斷電0.2 s后投入備用電源,BBA段母線初始電壓較 BBB、BBC、BBD段母線電壓高,母線電壓為69.35%額定電壓;而采用“規(guī)定”方法計(jì)算得,高壓母線初始電壓為83.56%基準(zhǔn)電壓。采用ETAP軟件仿真,工況4中斷電1.5 s后投入備用電源,BBB段母線初始電壓較BBA、BBC、BBD段母線電壓高,母線電壓為59.37%額定電壓;而采用“規(guī)定”方法計(jì)算后,高壓母線初始電壓為67.50%基準(zhǔn)電壓。2種方法的計(jì)算結(jié)果相差較大。
采用ETAP軟件仿真,工況1、2、3滿足俄方對(duì)廠用電設(shè)計(jì)的要求,采用“規(guī)定”方法計(jì)算結(jié)果,也滿足成功自起動(dòng)要求,2種計(jì)算結(jié)果一致。采用ETAP軟件仿真,工況4起動(dòng)時(shí)間為6 s左右,在5 s內(nèi)達(dá)到了65%額定電壓,計(jì)算結(jié)果滿足俄方對(duì)廠用電設(shè)計(jì)的要求;而采用“規(guī)定”方法計(jì)算結(jié)果則不滿足成功自起動(dòng)要求,2種計(jì)算結(jié)果產(chǎn)生了偏差。這是由于“規(guī)定”單純地將暫態(tài)過(guò)程簡(jiǎn)化為快速切換和慢速切換2種工況,勢(shì)必造成計(jì)算結(jié)果不合理。
(1)母線電壓校驗(yàn)公式不能準(zhǔn)確反映廠用母線電壓狀況,計(jì)算結(jié)果與實(shí)際暫態(tài)穩(wěn)定過(guò)程有一定的出入,具有一定的局限性,對(duì)于核電站和承建的國(guó)外電廠不完全適用。
(2)采用ETAP軟件對(duì)俄方提出的4種暫態(tài)工況進(jìn)行仿真,仿真結(jié)果滿足俄方對(duì)廠用電源的要求。提供的ETAP暫態(tài)穩(wěn)定分析報(bào)告得到了俄方和業(yè)主的認(rèn)可,并根據(jù)計(jì)算結(jié)果開(kāi)展下一步設(shè)計(jì)工作,此仿真結(jié)果可為繼電保護(hù)調(diào)整提供理論依據(jù)。
(3)田灣核電站3、4號(hào)機(jī)組的暫態(tài)穩(wěn)定分析表明,ETAP軟件可以很好地仿真電源切換過(guò)程中的母線電壓,并提供電壓曲線上各點(diǎn)的數(shù)據(jù),據(jù)此可判斷廠用電系統(tǒng)受到的擾動(dòng)是否能達(dá)到穩(wěn)定。仿真結(jié)果是廠用電系統(tǒng)繼電保護(hù)調(diào)整的依據(jù),打破了傳統(tǒng)的實(shí)驗(yàn)及經(jīng)驗(yàn)整定模式。
(4)應(yīng)用本研究成果,在核電工程廠用電系統(tǒng)設(shè)計(jì)時(shí),可以很好地實(shí)現(xiàn)暫態(tài)穩(wěn)定分析,同時(shí)還可供火電工程參照使用,尤其對(duì)大容量空冷機(jī)組的設(shè)計(jì)更具有實(shí)用意義。
(5)本研究成果提高了核電機(jī)組廠用電系統(tǒng)供電的設(shè)計(jì)可靠性水平,實(shí)現(xiàn)了從靜態(tài)估算向暫態(tài)穩(wěn)定仿真的跨進(jìn)。
[1]李亞樓,周孝信,吳中習(xí).基于PC機(jī)群的電力系統(tǒng)機(jī)電暫態(tài)仿真并行算法[J].電網(wǎng)技術(shù),2003,27(11):6-12.
[2]段剛,余貽鑫,殷啟志.廠用電切換機(jī)電動(dòng)態(tài)過(guò)程研究[J].電網(wǎng)技術(shù),1998,22(1):61-67.
[3]任果.基于PSCAD/EMTDC的廠用電源切換研究[D].天津:天津大學(xué),2007.
[4]張曉輝,胡曉芬.淺析大功率異步電動(dòng)機(jī)起動(dòng)對(duì)二灘電廠廠用電系統(tǒng)電壓的影響[J].電力系統(tǒng)保護(hù)與控制,2009,37(16):112-114.
[5]錢序,萬(wàn)里寧,楊華.基于ETAP的發(fā)電廠電氣系統(tǒng)暫態(tài)分析[J].電力建設(shè),2012,33(3):74-78.
[6]俄羅斯設(shè)計(jì)院.技術(shù)方案的電氣技術(shù)部分:廠用電系統(tǒng)[R].俄羅斯:俄羅斯設(shè)計(jì)院,1999.
[7]東北電力設(shè)計(jì)院.基于ETAP仿真的核電機(jī)組電源切換中的暫態(tài)分析研究[R].長(zhǎng)春:東北電力設(shè)計(jì)院,2012.
[8]DL/T 5153—2002火力發(fā)電廠廠用電設(shè)計(jì)技術(shù)規(guī)定[S].北京:中國(guó)電力出版社,2002.
[9]李光琦.電力系統(tǒng)暫態(tài)分析[M].北京:中國(guó)電力出版社,2000.
[10]李發(fā)海,朱東起.電機(jī)學(xué)[M].北京:科學(xué)出版社,2007.