馬宜軍-VFTO下頻變變壓器繞組過電壓時(shí)域算法研究

分類號： 學(xué)校代碼：10079 密級：華北電力大學(xué)碩 士 學(xué) 位 論 文題 目：陡波前過電壓下頻變變壓器繞組過電壓時(shí)域算法研究英文題目：ime-domain Calculation of ransformer indingsOvervoltage under VFTO Considering Frequency-dependentParameters研究生姓名：馬宜軍 專業(yè)：電工理論與新技術(shù)研究方向：電網(wǎng)絡(luò)理論及其在電力系統(tǒng)中的應(yīng)用導(dǎo)師姓名：梁貴書、董華英 職稱：教授2006年12月22日聲 明本人鄭重聲明：此處所提交的碩士學(xué)位論文《陡波前過電壓下頻變變壓器繞組過電壓時(shí)域算法研究，是本人在華北電力大學(xué)攻讀碩士學(xué)位期間，在導(dǎo)師指導(dǎo)下進(jìn)行的研究工作和取得的研究成果據(jù)本人所知除了文中特別加以標(biāo)注和致謝之處外論文中不包含其他人已經(jīng)發(fā)表或撰寫過的研究成果也不包含為獲得華北電力大學(xué)或其他教育機(jī)構(gòu)的學(xué)位或證書而使用過的材料與我一同工作的同志對本研究所做的任何貢獻(xiàn)均已在論文中作了明確的說明并表示了謝意。學(xué)位論文作者簽名： 日 期： 本人完全了解華北電力大學(xué)有關(guān)保留使用學(xué)位論文的規(guī)定即①學(xué)校有權(quán)保管、并向有關(guān)部門送交學(xué)位論文的原件與復(fù)印件②學(xué)?？梢圆捎糜坝】s印或其它復(fù)制手段復(fù)制并保存學(xué)位論文③學(xué)?？稍试S學(xué)位論文被查閱或借閱④學(xué)?？梢詫W(xué)術(shù)交流為目的,復(fù)制贈(zèng)送和交換學(xué)位論文⑤同意學(xué)?？梢杂貌煌绞皆诓煌襟w上發(fā)表傳播學(xué)位論文的全部或部分內(nèi)容。(涉密的學(xué)位論文在解密后遵守此規(guī)定)作者簽名： 導(dǎo)師簽名 日 期： 日 期： 華北電力大學(xué)碩士學(xué)位論文目錄華北電力大學(xué)碩士學(xué)位論文目錄目 錄中文摘要英文摘要第一章緒論11.1陡波前過電壓（VFTO）11.1.1VFTO的產(chǎn)生11.1.2VFTO的波形特點(diǎn)與產(chǎn)生機(jī)理11.1.3VFTO對變壓器絕緣的影響21.2本課題國內(nèi)外研究概況31.2.1變壓器模型的建立31.2.2多導(dǎo)體傳輸線方程的數(shù)值解法31.3本文完成的主要工作5第二章預(yù)備知識62.1多導(dǎo)體傳輸線的基本知識62.2變壓器分布參數(shù)模型的建立62.3變壓器線圈MTL模型基本參數(shù)的計(jì)算8第三章VFTO下變壓器線圈的空間離散模型103.1緊湊有限差分法103.2變壓器多導(dǎo)體傳輸線模型的空間離散123.3矢量匹配和遞歸卷積133.4稀疏矩陣技術(shù)16第四章頻變變壓器的線圈時(shí)域仿真184.1直接求解微分方程仿真電位184.1.1精細(xì)積分法184.1.2仿真與實(shí)驗(yàn)驗(yàn)證194.2微分方程離散為代數(shù)方程仿真電位254.2.1反向差分法254.2.2BDF應(yīng)用于變壓器線圈284.2.3仿真驗(yàn)證與實(shí)例計(jì)算294.3變壓器線圈餅間并聯(lián)MOV限制VFTO幅值36第五章基于krylov子空間的時(shí)域降階算法40iiii5.1基于krylov子空間的代數(shù)方程的降階405.1.1子空間簡述405.1.2Krylov子空間415.1.3Arnoldi算法425.1.4Lanczos算法445.2仿真驗(yàn)證46第六章結(jié)論與展望47參考文獻(xiàn)48致謝52攻讀學(xué)位期間發(fā)表的學(xué)術(shù)論文和參加科研情況53華北電力大學(xué)碩士學(xué)位論文華北電力大學(xué)碩士學(xué)位論文PAGEPAGE1第一章緒論1.1陡前過電壓（VFO）1.1.1VFO氣體絕緣變電站(GasInsulatedSubstation簡稱GIS),又稱為全封閉組合電器[1]，GIS是將高壓電器放于接地的金屬殼內(nèi)，以高壓氣體為主要絕緣的電站。所用的絕緣氣體主要是SF6。GIS具有結(jié)構(gòu)緊湊、體積小、占地少、運(yùn)行可靠、維護(hù)工作量少、對環(huán)境污染小等優(yōu)點(diǎn)，為城市和人口稠密地區(qū)以及大型水電工程建設(shè)節(jié)省用地創(chuàng)造了條件，所以在電力系統(tǒng)中得到越來越廣泛的應(yīng)用。在GIS變電站中由于斷路器隔離開關(guān)以及接地開關(guān)操作或帶電線路對地閃絡(luò)，甚至雷電波入侵，都可能會在GIS內(nèi)部產(chǎn)生一個(gè)上升速度極快(幾到幾十ns)的電壓陡波這個(gè)電壓陡波沿著GIS管道傳播遇到波阻抗發(fā)生改變就會發(fā)生反射和折射所有在GIS中產(chǎn)生的多次反射和折射的各行波分量疊加在一起就會形成波頭很陡、頻率高達(dá)幾到幾十MHz甚至上百M(fèi)Hz數(shù)量級的陡波前過電壓，又稱特快速暫態(tài)過電壓(eryFastransientOvervoltage，簡稱VFTO)[2]。陡波前過電壓主要包括內(nèi)部陡波前過電壓、瞬態(tài)外殼電壓和外部陡波前過電壓三種在GIS內(nèi)高壓導(dǎo)體或管道和外殼之間產(chǎn)生的VFFO稱內(nèi)部陡波過電壓出現(xiàn)在GIS外部及GIS以外設(shè)備上或設(shè)備內(nèi)的陡波前過電壓稱之為外部陡波前過電壓；在GIS外殼和地之間產(chǎn)生的陡波前過電壓稱為瞬態(tài)外殼電壓[3]。1.1.2VFO的波形特點(diǎn)與產(chǎn)生機(jī)理VFTO波形包括多個(gè)頻率段，波形比較復(fù)雜，通常由四個(gè)分量組成[4]：(1)階躍電壓；(2)由于GIS內(nèi)母線管道(即電暈屏蔽彎管等)波阻抗的多次微弱變化形成的極高頻范圍(最高達(dá)100MHz)；(3)由于GIS母線管道和電纜末端或架空線終端處波阻抗的顯著變化而引起的反射形成的高頻范圍(最高達(dá)30MHz)；(4)由于外部設(shè)備的大電容(如電容式互感器或輸電線載波系統(tǒng)的耦合電容)引起的諧振產(chǎn)生的低頻范圍(0.1～5MHz)。因此，內(nèi)部陡波前過電壓的波形取決于GIS的內(nèi)部結(jié)構(gòu)和外部配置。此外，由于陡波前過電壓的行波特性，其波形隨位置不同可能有很大的變化(在某些情況下，1米的距離就會造成顯著的變化[4])通常情況下并不是因?yàn)閂FTO幅值特別高而是其高頻振蕩的電壓波對變壓器的絕緣和系統(tǒng)的正常運(yùn)行造成危害[5]。1.1.3VFO對變壓器絕緣的影響當(dāng)GIS進(jìn)行斷路器隔離開關(guān)等開關(guān)操作時(shí)產(chǎn)生的陡波前過電壓極有可能引起GIS內(nèi)部設(shè)備及與GIS相連的電力系統(tǒng)高壓設(shè)備的絕緣事故[6]。決定該類絕緣事故的主要因素有[7,8]：(1)開關(guān)在操作期間的動(dòng)態(tài)絕緣特性；(2)高壓設(shè)備的絕緣介質(zhì)在VFTO作用下的絕緣耐受能力。通過對VFTO的產(chǎn)生機(jī)理、波形及傳播途徑的分析可知，VFTO的特性與高壓設(shè)備耐壓試驗(yàn)所采用的雷電沖擊波和操作沖擊波的特性有很大不同。盡管制造商已掌握了對高壓設(shè)備(如變壓器)在標(biāo)準(zhǔn)雷電沖擊波和操作波作用下的各種絕緣數(shù)據(jù)，但是，由于VFTO的高陡度及其高頻振蕩特性，它會怎樣影響絕緣強(qiáng)度，絕緣的擊穿強(qiáng)度是多少，絕緣裕度應(yīng)當(dāng)怎樣選取等等問題，制造者和設(shè)計(jì)者目前還不十分清楚[9]。在各種電力設(shè)備中，變壓器受VFTO影響非常大，因?yàn)閂FTO電壓波頻率高，當(dāng)幅值不是特別高時(shí)，甚至不足以使變壓器入口避雷器動(dòng)作的情況下，就可能侵入變壓器造成匝間絕緣的破壞因此不論是與GIS直接連接的還是非直接連接的變壓器，其內(nèi)部的絕緣都會受到很大的威脅。但是當(dāng)變壓器不是直接和開關(guān)相連接的時(shí)候，VFTO會有一定削弱。由陡波前過電壓引起的絕緣擊穿事故已在許多國家電網(wǎng)中出現(xiàn)，如美國電力系統(tǒng)(ACP)、加拿大及中國廣東核電站均發(fā)生過由于此類因素引起的超高壓變壓器絕緣擊穿事故[3]對于VFTO造成的變壓器內(nèi)部的絕緣擊穿，其原因就目前的研究而言，可以歸納為以下兩個(gè)方面[10,1]：(1)在開關(guān)觸頭擊穿瞬間產(chǎn)生的VFTO到達(dá)變壓器時(shí)，在變壓器端部加上了一個(gè)陡波頭波，對直接連接的變壓器，其上升時(shí)間可能只有數(shù)十ns，遠(yuǎn)遠(yuǎn)低于雷電沖擊波試驗(yàn)時(shí)的波頭上升時(shí)間(約為1.2μs)，該陡波的波頭在變壓器繞組上造成極不均勻的匝間電壓分布，大部分電壓降落在靠近入波端的一小部分線圈或?qū)w上，電位梯度極大，危害較大。對非直接連接的變壓器，因?yàn)榻?jīng)過了其它的設(shè)備，陡波頭趨于平緩，其作用與雷電沖擊波相近。(2)VFTO中含有的振蕩諧波的頻率與變壓器中的若干固有振蕩頻率匹配從而引起諧振產(chǎn)生幅值很高的高頻諧振過電壓導(dǎo)致繞組與鐵芯以及匝間的絕緣破壞。由于變壓器的絕緣設(shè)計(jì)只考慮了雷電沖擊波以及操作波作用下引起的沖擊情況，故匝間絕緣在高頻振蕩情況下顯得很脆弱。1.2本題國內(nèi)外研究概況1.2.1變壓器模型的建立變壓器繞組建模的方法有場的方法[12]、場-路結(jié)合的方法[13]和路的方法[14-17]三大類。但由于變壓器線圈的纏繞結(jié)構(gòu)復(fù)雜，采用前兩種方法計(jì)算量龐大。將“場”問題簡化為“路”來求解，可有效地減少計(jì)算量。實(shí)驗(yàn)和仿真結(jié)果的吻合表明，從工程應(yīng)用的角度來說這種簡化是合理、可以接受的。因此，實(shí)際中主要采用路的方法。采用路的方法建立的模型又分為以下幾類：(1)低頻和中頻下繞組模型低頻和中頻下一般采用集中參數(shù)模型。文獻(xiàn)[3]基于變壓器單餅或雙餅為單元建立了等值集中電路模型，利用該模型仿真雷電波侵入變壓器的暫態(tài)過程的時(shí)候可以滿足需要。如果用此模型仿真變壓器VFTO作用下暫態(tài)電壓分布存在以下問題，首先仿真精度無法滿足要求；其次只能計(jì)算線餅首末端的過電壓值，對于線匝之間的電壓分布，不能通過該電路模型計(jì)算獲得[8]。此外，還有很多集中參數(shù)的變壓器模型，例如：EMTP下的BCTRAN和TRELEG這兩個(gè)標(biāo)準(zhǔn)的變壓器模型[14]。(2)高頻下繞組模型頻率高的時(shí)候，線匝長度相對電磁波的波長已經(jīng)不能忽略，因此需要采用分布參數(shù)理論建立變壓器繞組的高頻模型[14-21]建立變壓器繞組的電路高頻模型主要有兩種方法：基于多端口網(wǎng)絡(luò)理論建立變壓器線圈的集中參數(shù)電路模型[14-16]；通過分割變壓器線圈為若干單元來建立詳細(xì)的內(nèi)部模型[17-21]。文獻(xiàn)[17]采用集中參數(shù)元件構(gòu)造模型在最容易損壞的前幾匝線圈,每匝線圈為一單元構(gòu)造模型,其余線圈以兩餅為一單元但是文中實(shí)驗(yàn)驗(yàn)證的輸入波采用波前時(shí)間為μs級的雷電波,而VFTO下波前時(shí)間達(dá)到ns級；文獻(xiàn)[18]采用分布參數(shù)與集中參數(shù)相結(jié)合的混合模型，模型中集中電路部分的元件參數(shù)不易確定；文獻(xiàn)[19]對殼式變壓器線圈采用單導(dǎo)體傳輸線和多導(dǎo)體傳輸線相結(jié)合的模型，該模型僅考慮了餅內(nèi)互感，未考慮餅間互感；實(shí)踐表明，該模型應(yīng)用于其它類型變壓器線圈時(shí)，仿真計(jì)算結(jié)果與實(shí)際測量存在較大誤差[17]文獻(xiàn)[20-21]將變壓器的每一匝線圈看成一條傳輸線建立了繞組的全多導(dǎo)體傳輸線模型可以詳細(xì)反映變壓器線圈的電磁過程。1.2.2多導(dǎo)體傳輸線方程的數(shù)值解法由于VFTO頻率很高變壓器線圈參數(shù)呈現(xiàn)明顯的頻變效應(yīng)[19]因此一般采用頻域分析方法對多導(dǎo)體傳輸線模型進(jìn)行分析[19-21]即借助快速傅立葉變換首先求其頻域響應(yīng)，然后轉(zhuǎn)換到時(shí)域。這種頻域分析方法對于大型變壓器計(jì)算量龐大，且難以處理非線性問題。為了解決這一問題，本文采用時(shí)域方法分析頻變線圈。其基本思想是直接在時(shí)域用數(shù)值方法將傳輸線的偏微分方程通過空間變量的離散，轉(zhuǎn)化為常微分方程求解。求解多導(dǎo)體傳輸線方程的時(shí)域數(shù)值分析方法主要有以下幾種：(1)時(shí)域有限差分法(FDTD)[22]，采用具有二階精度的中心差分近似替代偏微分方程中的空間和時(shí)間導(dǎo)數(shù)，即dpl)=dl

pl+Δ/2)?pl?Δ/2)+ο(Δ2)Δ其中，Δ是步長間隔并且要構(gòu)造差分網(wǎng)格用以實(shí)現(xiàn)中心差分近似這是FDTD的關(guān)鍵步驟。該方法簡單直接，但是為了得到理想的精度，差分點(diǎn)數(shù)非常密集，降低了計(jì)算效率。(2)特征法[24]，將具有時(shí)間t和空間坐標(biāo)x的偏微分方程轉(zhuǎn)換到兩組特征線族上以常微分方程的形式求解，具有較高的計(jì)算效率。該方法最基本特點(diǎn)是將電壓電流分解為入射波和反射波，且各具有波速Φ,這對于單根傳輸線的方程容易做到。對于單根線，傳輸線越簡單，該方法優(yōu)勢越明顯，計(jì)算均勻無耗線時(shí)，其計(jì)算效率遠(yuǎn)遠(yuǎn)高于其他方法。當(dāng)傳輸線較為復(fù)雜時(shí)，例如對于單根有耗線或非均勻線，相對優(yōu)勢就減少了。當(dāng)傳輸線參數(shù)高度頻變時(shí)，該方法就難以處理了。在實(shí)際的系統(tǒng)中，大多數(shù)都是多導(dǎo)體耦合線，因?yàn)樘卣鞣ū举|(zhì)是將電壓和電流分解為入射波和反射波，因此需要對多導(dǎo)體耦合線解耦，對無耗耦合線容易做到。對于有耗多導(dǎo)體耦合線就比較復(fù)雜，電報(bào)方程右側(cè)有兩項(xiàng)，兩個(gè)分布參數(shù)矩陣不相同，方程組去耦時(shí)要求對兩個(gè)分布參數(shù)矩陣同時(shí)對角化，在一般情況下比較困難，只有某些特定的多導(dǎo)體耦合系統(tǒng)能做到，這就在很大程度上限制了特征法的應(yīng)用范圍。(3)微分求積法(DQ)[23]，主要思想是將某點(diǎn)對座標(biāo)的微分算子以該座標(biāo)全部定義域中一系列離散點(diǎn)的函數(shù)值加權(quán)逼近，將微分方程化為常微分方程或代數(shù)方程求解。其主要環(huán)節(jié)在于利用一組試驗(yàn)函數(shù)確定微分求積逼近式中的線性加權(quán)全系數(shù)。該方法得到的方程組所包含的微分方程數(shù)可明顯少于前述由差分而轉(zhuǎn)化成的微分方程組的方程數(shù)，但是DQ方法中試驗(yàn)函數(shù)的選擇至關(guān)重要，試驗(yàn)函數(shù)有很多種類，沒有任何試驗(yàn)函數(shù)適用所有仿真對象，當(dāng)所選試驗(yàn)函數(shù)不適合仿真對象時(shí)，就會嚴(yán)重影響精度。(4)緊湊有限差分法[2],主要思想是采用四階精度的插值公式αf(x)

+αf(x)

+αf(x)

=fi12?fi121 2 1xi1

xi

xi1 x對時(shí)域電報(bào)方程空間離散化為微分方程式中f代表電壓或者電流變量α1,α2為根據(jù)采用的插值公式?jīng)Q定的常數(shù)。該方法離散點(diǎn)僅是傳統(tǒng)離散方法的1/3，極大地提高了計(jì)算效率。1.3本完成的主要工作本文主要研究VFTO下變壓器繞組的傳輸線高頻電路模型的時(shí)域快速求解方法，確定變壓器在VFTO作用下繞組內(nèi)部的電位分布。這對于改進(jìn)變壓器的絕緣設(shè)計(jì)從而提高電力系統(tǒng)的穩(wěn)定運(yùn)行將起到重要的作用。本文的具體工作如下：(1)在時(shí)域計(jì)算頻變參數(shù)是難點(diǎn)，本文采用矢量匹配結(jié)合遞歸卷積在處理時(shí)域頻變問題，在保證高的計(jì)算精度情況下，避免了常規(guī)卷積積分的巨大計(jì)算量。使用緊湊有限差分法在空間離散電報(bào)方程使用精細(xì)積分法求解微分方程得到了VFTO下頻變變壓器繞組的時(shí)域響應(yīng)。(2)基于(1)中的理論，采用了反向差分法求解微分－代數(shù)方程組，仿真了內(nèi)屏蔽式變壓器線圈VFTO下的響應(yīng)，極大地減少了計(jì)算量。利用時(shí)域算法可以處理非線性的優(yōu)點(diǎn)，仿真了MOV做變壓器的內(nèi)保護(hù)時(shí)線圈的電位分布,仿真結(jié)果表明，線圈諧振電壓幅值得到了有效的抑制。(3)研究了Krylov子空間降階方法，以降低求解大規(guī)模代數(shù)方程的計(jì)算量，并且減少內(nèi)存的使用仿真結(jié)果表明在保證精度的基礎(chǔ)上方程的階數(shù)大幅降低。第二章預(yù)備知識2.1多體傳輸線的基本知識當(dāng)信號頻率提高時(shí)，導(dǎo)體的長度相對電磁波的波長不能忽略，就需要應(yīng)用傳輸線理論[25]研究問題傳輸線是約束電磁波沿著規(guī)定方向傳輸能量和信息的系統(tǒng)傳輸線的幾何長度l與電磁波的工作波長λ之比值lλ稱為傳輸線的電長度，通常把lλ>.05的傳輸線稱為長線。這時(shí)傳輸線導(dǎo)體上存在的損耗電阻，兩導(dǎo)體間介質(zhì)損耗產(chǎn)生的電導(dǎo)、傳輸線的自感以及兩導(dǎo)線間的互電容，這些量沿線分布，因而這些量稱為分布參數(shù)。在傳輸線方程的建立以及分析中，最基本的假設(shè)是所分析的傳輸線為TEM結(jié)構(gòu),即電場和磁場的方向?yàn)闄M向。傳輸線的數(shù)學(xué)表達(dá)形式為電報(bào)方程，其時(shí)域形式為??u=?i?L?i??x ?t?? ??i=?u?Cu??x ?t其中，R,G,L,C分別為傳輸線單位長度的電阻電導(dǎo)電感電容如果初始電壓、電流值為零，頻域電報(bào)方程可寫成下列形式?d?dx?

=?(R+jωL)=?Z?d=?(G+jωC)=?dx其中Z=R+ωL為單位長度的串聯(lián)阻抗，Y=G+ωC為單位長度的并聯(lián)導(dǎo)納。2.2變器分布參數(shù)模型的建立為了建立變壓器線圈的多導(dǎo)體傳輸線(MTL,MulticondutorransmissionLines)模型，對變壓器線圈進(jìn)行如下理想化假設(shè)[26]：(1)認(rèn)為線圈的平均直徑遠(yuǎn)大于其幅向尺寸(繞組的徑向?qū)挾?因而可以忽略線匝彎曲的影響，且所有線匝的長度近似相等(取平均值，即平均匝長)。(2)認(rèn)為線餅間連線及電壓源引線充分短，因而可以忽略它們對電磁場分布的影響。(3)認(rèn)為線圈的平均匝長大于所分析線圈的截面尺寸，因而認(rèn)為電磁波沿線匝傳播過程在同一子午面上是瞬時(shí)建立起的電磁場分布，即忽略電磁波沿軸向、徑向的延時(shí)效應(yīng)。根據(jù)上面的假設(shè)，可將線圈在線端處沿子午面剖開，將線匝展成直線，每一線匝成一“傳輸線對于連續(xù)式線圈模型這些傳輸線按線圈繞制關(guān)系首尾相連(如圖2-1所示)。為了方便，以下仍將“傳輸線”稱為線匝。規(guī)定線匝按電氣聯(lián)接順序編號。當(dāng)線匝連接順序改變時(shí)（線圈類型改變，線匝編號隨之改變。這樣，線匝的邊界條件為第i根線末端的電壓電流分別等于第i＋1根線首端的電壓電（i＝1，2，…，N－1；第一根線首端接電壓源，第N根線末端或接地、或接負(fù)載阻抗、或懸空，邊界條件總數(shù)為2N個(gè)。US

IS

IR)

URUS()

IS(2)

IR(2)

UR(2)US(N?)

IS(N?)

IR(N?)

UR(N?)US(N)

IS(N)

IR(N)

UR(N)圖2-1線圈首位相連示意圖沿用傳統(tǒng)的建立MTL方程的方法，可得?U(x,t)=?I(x,t)?LI(x,t)? x t?

(2-1)?I(x,t)=GU(x,t)?CU(x,t)x t其中U(x,t)和I(x,t)是沿線分布的N×1的電壓、電流列向量．R，L，C，G分別是單位長度上的電阻、電感、電容和電導(dǎo)矩陣。將上式變換到頻域可得頻域MTL方程為?dU(x,s)=?(R+sL)I(x,s)=?Z(s)I(x,s)? dx?

(2-2)?dI(x,s)=?G+s)U(x,s)=?Y(s)U(x,s)? dxZ和Y分別是N×N串聯(lián)阻抗矩陣和并聯(lián)導(dǎo)納矩陣。對于圖2.1所示的連續(xù)式線圈的MTL模型，其邊界條件可以描述為iri)=isi+)uri)=usi+)

i=,,",N?1us)=ust)（變壓器入口VFTO波形）上式中下標(biāo)r代表線匝末端，s代表線匝首端。最后加上線圈最后一匝末端的邊界條件，根據(jù)線圈末端邊界條件的不同(可能是直接接地，懸空，經(jīng)過阻抗接地等)，添加不同的邊界條件方程(可能是電壓為零電流為零或電壓和電流滿足一定的約束條件)。求解MTL上的電壓電流分布情況，通?？煞譃橄铝腥齻€(gè)步驟[27]：(1)從靜電場的角度出發(fā)，計(jì)算傳輸線分布參數(shù)模型的單位長度參數(shù)。(2)從電路角度建立MTL方程并添加邊界條件；(3)求解第二步中在邊界條件約束下的MTL方程。2.3變器線圈MTL模型基本參數(shù)的計(jì)算(1)電容參數(shù)的確定文獻(xiàn)[27]中使用有限元軟件FEMM計(jì)算靜電場中變壓器線匝間儲存的能量計(jì)算時(shí)候激勵(lì)導(dǎo)體加1V電壓，其余導(dǎo)體電勢為零，可以由如下公式計(jì)算電容1N12∑ijΔij2

=Wji12其中Wj為第j個(gè)激勵(lì)導(dǎo)體產(chǎn)生的能量，i為與激勵(lì)導(dǎo)體之間有電容的導(dǎo)體編號。因?yàn)棣j值為1所以當(dāng)有足夠數(shù)目的類似上述方程聯(lián)立就可以計(jì)算出匝間餅間電容，從而形成線圈的電容矩陣C。(2)電阻參數(shù)的確定計(jì)算電阻參數(shù)，必須考慮高頻下的集膚效應(yīng)。高頻時(shí)單位長度電阻為[25]R=0+S其中，0表示靜態(tài)情況下的電阻，S是由于集膚效應(yīng)導(dǎo)致的頻變電阻。fs=式中

δ(1+d2)

(2-3)δ=1

1d22(1+d2)1，d2分別為矩形導(dǎo)體橫截面的長和寬，σ為導(dǎo)體的電導(dǎo)率，f為對應(yīng)的頻率。為了把頻域的線圈模型轉(zhuǎn)化到時(shí)域，需要對頻變電阻做如下處理，把f代入s表達(dá)式(2-3)中

=ωπs sR= f = 1s s

ωR′ωσδ(1+d2) σδ(1+d2) π其中把s=jω代入上式

′=

1 1δ(1+d2) πR=R′

ω=R′ jω=

2R′

jω=

2R′

s(1?j)s s s j

1+js 2 s(3)電感參數(shù)的確定同電阻參數(shù)類似，高頻時(shí)電感參數(shù)也是由兩部分構(gòu)成的[25]。L=e+i其中，e是低頻情況下的電感，i是在高頻時(shí)由于集膚效應(yīng)造成的電感。e的計(jì)算公式為L=εrC1e c2i其中，C為(1)中所計(jì)算的電容矩陣，εr為絕緣材料的相對介電常數(shù)，c為光在真空中的速度，取c=3×8/s。L其表達(dá)式如下[28]：iL=si ω其中，s為電阻的頻變部分，ω為相應(yīng)的頻率。線圈的串聯(lián)阻抗為Z=R+jωL=0+S+jω(e+i)=0+jωe+S+jωi=0+jωe+S+jS把S的表達(dá)式(2-3)代入上式，可以得到Z=0+se+

2′s上式右側(cè)前兩項(xiàng)可以使用反拉氏變換直接化到時(shí)域，第三項(xiàng)通過矢量匹配和遞歸卷積處理，具體方法將在以后章節(jié)詳細(xì)敘述。(4)電導(dǎo)參數(shù)的確定電導(dǎo)G可以由下式?jīng)Q定[28]G=σCε其中，σ是導(dǎo)體電導(dǎo)率；ε是周圍介質(zhì)介電常數(shù)，C是電容。根據(jù)文獻(xiàn)[25]，當(dāng)頻率低于吉赫茲的時(shí)候，電導(dǎo)受頻變因素的影響可以忽略，VFTO信號最高在百z，所以本文不考慮電導(dǎo)的頻變。第三章VFTO下變壓器線圈的空間離散模型本章以連續(xù)式線圈為例，詳細(xì)敘述了VFTO下頻變線圈的空間離散。本章首先使用緊湊有限差分法(CompactFiniteDifference,CFD)[29-30]對變壓器繞組的多導(dǎo)體傳輸線模型進(jìn)行空間離散運(yùn)用矢量匹配法[31]處理頻變參數(shù)并用遞歸卷積[32-34]處理由此引起的卷積項(xiàng)，從而把帶有卷積積分的偏微分方程化為狀態(tài)方程；由于矢量匹配可以較為準(zhǔn)確地逼近頻變參數(shù)緊湊有限差分法只需取傳統(tǒng)分段方法1/3的段數(shù)，空間離散點(diǎn)少；遞歸卷積只需考慮前一步的值，避免了普通卷積積分大量的計(jì)算，因此，基于本章方法的時(shí)域仿真算法不僅可以方便地處理非線性，而且可有效地減少計(jì)算量。3.1緊有限差分法變壓器多導(dǎo)體模型在時(shí)域的數(shù)學(xué)形式是偏微分方程組CFD就是用來空間離散這個(gè)偏微分方程組，把模型數(shù)學(xué)形式化作微分方程組。為了敘述的簡便，以單導(dǎo)體情況為例進(jìn)行方法的說明。設(shè)單根傳輸線的長度為l，如圖3-1所示。沿空間將其等分為M段，每段長度為Δx=l

M。u(x)表示在x=(n?12(Δx)(n=,,??,M)處的電壓i(x)表示在x=n(Δx)(n=,,??,M)處的電流。i0u0

xu12

i1u212 i2

iM?2uM?32iM1uM12iMuM圖3-1傳輸線的分段傳輸線的時(shí)域電報(bào)方程為?u(x,t)=?Li(x,t)?i(x,t)? x t?

(3-1)?i(x,t)=Cu(x,t)?u(x,t)x t式中，R,G,L和C分別為傳輸線單位長度的電阻、電導(dǎo)、電感和電容,u和i分別為線上的電壓和電流,x為長度上空間坐標(biāo)，t為時(shí)間變量。使用下式對電報(bào)方程進(jìn)行空間離散[29]αf(x)

+αf(x)

+αf(x)

=fn12?

fn121 2 1

(3-2)x n1

x n

x n1 Δx[30]式中，f代表u(x)或者i(x)；α1,α2為根據(jù)采用的插值公式?jīng)Q定的常數(shù)

。當(dāng)采用四階精度的插值公式時(shí)，α1=1/24，α2=1/12。使用式(3-2)離散式(3-1)中第一個(gè)方程，令f(x)=u(x,t)，式(3-2)可化為αu(x,t)

+αu(x,t)

+αu(x,t)

u ?u=n12 n121 2 1x n1

x n

x n1 x把式(3-1)中第一個(gè)方程等號右側(cè)代入上式等號左側(cè)，經(jīng)過化簡α

+Ldn1＋α

i+Ldn＋α

+Ldn1

u ?u=－n12 n12

(3-3a)1( n1

dt) 2( n

dt) 1(

n1 dt) x同理，令f(x)=i(x,t)，可以利用式(3-2)離散方程式(3-1)中第二個(gè)方程α

du+C n+32＋α

u

du+C n12＋α

du+C n12)1 n+32 dt

2 n12 dt

1 n12 dt

(3-3b)=－n1?n(n=,??,M?2)x由于式(3-3)中，每點(diǎn)的值和其前后點(diǎn)的值相關(guān)，但是首端沒有前一點(diǎn)，末端沒有后一點(diǎn)因此空間離散方程(3-3)對傳輸線兩端不適用本文對兩端采用二階精度的離散公式[30]。對于首端有?α(Ri

+Ld0)?α(RI

+Ld1)=12?u0? 4 0 dt?

1 1 dt Δx

(3-4)?α

Gu

du+C 12)?αGu

du i?i+C 32)=1 0對于末端有

3 12 dt

1 32 dt Δx?α

(Ri

+LdM)?α(Ri

+LdM1)=uM?uM1/2? 4 M dt

1 M1 dt Δx? (3-5)??α

Gu

+C M12)?αu

u i+C M?32)=M

M13 M12 dt

1 M?32 dt x其中，α3=1?α1,α4=1/2?α1。經(jīng)過上述空間離散得到2M+1個(gè)方程，其矩陣形式為Pdx+Qx+f=0dt

(3-6)其中 x=?12 u32

"uM1/2 0 1

TT"M ;Tf=[0 " 0

u0

" uM];?3 1 ?? ??1

a2 % ?? %%% 0 ?? ?? % a2 1 ?? 1 3 ?P=? ?? 4? 0 1

1 ?2 % ?? ?? % % % ?? ?? % 2? 1

1?4?3 1

?1 1 ?? ??1

2 % %% 0 ?? %%% %% ?? ?? %2 1

%% ?? 1 3

?1 1?? ?Q=?1??1% 0

4 1 0 ?1 2 % ?? ?? %% %%% ?? %% %%% ?? ?? 0 % 1

% 2

1?? ?? ?1

1 4?式中，i=C?Δx?αi，i=L?Δx?αi,i=G?Δx?αi,i=R?x?αi由于矩陣P逆矩陣存在[35]，則由方程(3-6)得

(i=1,2,3,4)。dx=x+f (3-7)P其中A=?P?Q，fP

dt P=?Pf。CFD和傳統(tǒng)差分方法相比，采用的點(diǎn)數(shù)僅為傳統(tǒng)方法的1/3；并且可以保證系統(tǒng)降階時(shí)的無源性。3.2變器多導(dǎo)體傳輸線模型的空間離散連續(xù)式變壓器線圈模型是由首尾相連的N根多導(dǎo)體傳輸線構(gòu)成[27]如圖2-1所示。多導(dǎo)體傳輸線模型的方程為?u(x,t)=?Li(x,t)?Ri(x,t)? x t?

(3-8)?i(x,t)=?Cu(x,t)?u(x,t)x t上式只是把式(3-1)中的

R,L,G,C變?yōu)镹×N矩陣，i(x)，u(x)變?yōu)镹維列向量。T1 2 N1 N其中 u(x)=[1(x,u2(x)???uN1(x,uNT1 2 N1 N

，i(x)=[i(x,i(x)???i (x,i(x)T。首尾相連的邊界條件為ijM

=(j)0,

ujM

=(j)0

(j=,2???N?);運(yùn)用緊湊有限差分法,由式(3-8)可得與式(3-6)相似的下列方程?d?+??+?=0 (3-9)dt上式中各量的意義與式(3-6)中相同只是，?和?均為(2M+)N維方陣，?和?均為(2M+)N維列向量。將電流邊界條件代入式(3-9)，合并同類項(xiàng)，消去未知數(shù)1M,2M,??,(N)M；將電壓邊界條件代入式(3-9)，把第2MN+1～(2M+)N?1個(gè)方程對應(yīng)加到第MN+2~(M+)N個(gè)方程，然后劃去第2MN+1～(2M+)N?1個(gè)方程。這樣，?向量中未知數(shù)1M,u2M???，(N)M被消去。至此，式(3-9)化為2MN+1維的方程組。?中10,uNM為已知的邊界條件。上述對方程的變換，相當(dāng)于對矩陣?和?把第2MN+1～(2M+)N?1列（行）對應(yīng)加到第MN+2~(M+)N（行然后劃去第2MN+1～(2M+)N?1（行。列向量?和x的處理方式與?和?中對行的處理類似。由于?可逆[35]，故方程(3-9)可變?yōu)橄铝袠?biāo)準(zhǔn)狀態(tài)方程由于?可逆，故方程(3-9)可變?yōu)橄铝袠?biāo)準(zhǔn)狀態(tài)方程d?=??+?

(3-10)dt P式中各項(xiàng)的意義與方程(3-7)相同。上述處理方法同樣適用于其它類型的變壓器線圈。3.3矢匹配和遞歸卷積VFTO作用下變壓器繞組傳輸線模型的分布參數(shù)呈現(xiàn)明顯的頻變特性。引入矢量匹配[1](ectorFitting)和遞歸卷積[2,33](ecursiveConvolution)相結(jié)合的方法處理頻變參數(shù)下面以集膚效應(yīng)導(dǎo)致的頻變效應(yīng)[25]為例進(jìn)行說明該方法對其他因素引起的參數(shù)隨頻率變化同樣適用。對頻變參數(shù)傳輸線，頻域中的電報(bào)方程為?dU(x,s)=?Z(s)I(x,s)? dx?

(3-1)?dI(x,s)=Y(s)U(x,s)? dx其中Z(s)=R(s)+sL(s)，Y(s)=G(s)+sC(s)。考慮集膚效應(yīng)時(shí)，單位長度電阻可表示為式(3-1)中第一式可寫為

R(s)=0+1sdU(x,s)=?R

+R s+sL

+sL(s))I(x,s)

(3-12)(dx 0 1 0 1(使用矢量匹配，上式中 s，1(s)化為如下形式N kf(s)≈∑ i i1s?i式中ki,i分別為第i個(gè)極點(diǎn)和相應(yīng)的留數(shù)。將上式代入式(3-12)得dU(x,s) ?R R

kR L s

kj

?I(x,s)n m∑ ∑=?? +

i + + ?

(3-13)dx ?

0 1 s?p 0 s?p?? i1 Ri

j1

j?式中n,m分別代表

s,L(s)匹配的階數(shù)，k,p分別表示相應(yīng)的極點(diǎn)和留數(shù)。對式(3-13)取拉氏反變換u(x,t)=?Ri(x,t)?L

i(x,t)?λ?λ

(3-14)x其中1,2均為卷積項(xiàng)

0 0 t 1 2?λ=λ(x,t)=

n tp(tτ)k e ix, d?1 1??

∑i1

Ri0

(3-15)?λ=λ

(x,t)=

m t p(tτ)Ri ( τ)τj ( Ri ( τ)τj ( )2 2 ∑j1

j0j對式(3-1)中第二個(gè)方程可類似處理。將具有N匝的變壓器線圈使用N根耦合頻變的多導(dǎo)體傳輸線系統(tǒng)建模把每根導(dǎo)體分為M段，得到形式如式(3-9)，傳輸線上分布電壓和電流空間離散方程?d?+??+?+?+?=0

(3-16)dt 1 212其中，?,?是由形如(3-14)的附加項(xiàng)（卷積項(xiàng)）形成的2MN+1維列向量12意義和(3-9)中相同。對任意空間點(diǎn)x，式(3-15)可以做如下變化：n m1(x,t)=∑i(x,t),2=∑2j(x,t)i1若時(shí)間步長為τ，在時(shí)刻n有如下遞推關(guān)系

j1?i(x,n+)=e?

pRτ

i(x,n)+di(x,n)+di2i(x,n+)

(3-17)2j(x,n+)=e

pLτ

2j(x,n)+dji(x,n)+dj2i(x,n+)k ? epRτ?1?

k ?epRτ?1 ?其中 di1=

Ri

?eRτ?

?，di2=

Ri?

?1?pi?

pτ ?

pi?

pτ ?d =k

?pτ?eLj?

epLτ?1??，d

?epτ?1=k?

??1?,j1 Lj?

pLτ ?

j2 Lj? pτ ?? ? ? Lj ?遞歸卷積關(guān)鍵是式(3-17)，對(3-17)中一式進(jìn)行證明，二式類似可證明。由式(3-15)n1 Ri( )n1 Ri( )λ(x,t

)=k

e i(x,t)dti n1 Ri0(=kRi(

t∫np(tτ)( )∫np(tτ)( )0

t∫n1p(tτ∫n1p(tτ)( ) )tn

(3-18)n1 Rin1 Ri( )

ptτ=eiλ

(x,t

)+k

e i(x,t)dtni nn令

Rit?(x,t

)=tn1eRi(tτ)i(x,t)dt

(3-19)ni n1 tn利用線性插值逼近傳輸線上分布電流關(guān)于時(shí)間t的函數(shù)，即假定電流在各時(shí)間點(diǎn)之間是線性的，則有如下結(jié)論i(x,t)=i(x,n)+i(x,n+)?i(x,n)(t?t)τ n

(3-20)把式(3-20)代入式(3-19)t+

ptτ?(x,t )=

e i(x,τ)dti n1

t?t+ pt?tn1n1

tn+

?tdt??ixt

+?epτ

n1t+ ptn1tet

t?tndt??ixtn1 Ri( )nn1 Ri( )n

Ri(n1) 1

? (,n) ?

Ri(n

) ? (

,n1)∫Ri?n τ ? ? n∫Ri

τ ?(3-21)使用分步積分法，式(3-21)可化為? 1?

epRτ?1?

1?epRτ?1 ?i(x,tn1)=

epRτ

? ?i(x,tn)+ ?1

?i(x,tn1)

(3-22)+ ? ? ? ?pi?

pτ ?

pi?

pτ ?將式(3-22)代入式(3-18)，即可得到式(3-17)中第一式。式(3-17)等號右邊的第三項(xiàng)可與方程中的同類未知量項(xiàng)合并，其余兩項(xiàng)可看作傳輸線上隨時(shí)間、空間變化的已知附加分布電源項(xiàng)，可與方程中的已知項(xiàng)合并。這樣，方程(3-16)可轉(zhuǎn)化為標(biāo)準(zhǔn)的狀態(tài)方程。3.4稀矩陣技術(shù)稀疏矩陣是指含有大量零元素，非零元素在全部矩陣元素中占比例很小的矩陣?？梢钥紤]利用稀疏矩陣內(nèi)含有大量零元素的特點(diǎn)，采取措施只存儲非零元素和只對非零元素進(jìn)行運(yùn)算減少存儲量和運(yùn)算量這一類方法稱之為稀疏矩陣技術(shù)[36]。電路方程中，除了個(gè)別的，大多數(shù)是稀疏矩陣。而且隨著電路規(guī)模的增加，矩陣的階數(shù)增高稀疏程度也越大例如10個(gè)節(jié)點(diǎn)的電路系數(shù)矩陣的非零元素可能占50%；而100個(gè)節(jié)點(diǎn)的典型電路中，非零元素僅占5%左右。據(jù)統(tǒng)計(jì)，對多數(shù)實(shí)際電路，其系數(shù)矩陣中的非零元素?cái)?shù)目在4n~6n之間(n為矩陣的階數(shù))。稀疏矩陣技術(shù)應(yīng)用主要優(yōu)點(diǎn)：(1)減少存儲量只存儲非零元素節(jié)省存儲空間如果矩陣是對稱的還可以進(jìn)一步節(jié)省存儲空間。(2)提高計(jì)算速度只進(jìn)行非零元素的運(yùn)算避免無效的零運(yùn)算采用合理的數(shù)據(jù)存儲技術(shù)，使運(yùn)算過程中對數(shù)據(jù)信息的分類、檢索、插入和刪除等處理盡可能方便，以提高數(shù)據(jù)檢索效率，減少運(yùn)算時(shí)間和存儲量。以上兩個(gè)問題的解決是相互依存的又可能是相互矛盾的往往需要綜合權(quán)衡，在存儲量、計(jì)算時(shí)間之間采取適當(dāng)?shù)恼壑苑桨?。為了感性認(rèn)識稀疏矩陣技術(shù)，在這里簡單介紹一種比較簡單的稀疏矩陣存儲方法：線性表存儲方法，對于稀疏矩陣T，可以使用3個(gè)向量α,β,γ分別存儲T中非零元素的橫坐標(biāo)，縱坐標(biāo)，數(shù)值。這樣就節(jié)省了大量的存儲空間。選擇存儲方案不僅考慮節(jié)省存儲空間，還要考慮存取時(shí)間，檢索、查找方便等因素，如果運(yùn)算過程中，稀疏矩陣是變化的，就要考慮刪除、添加、改變非零元素的效率，即存取效率問題。本文算法中，式(3-6),(3-16)中矩陣P、Q，基本的電阻矩陣R,電容矩陣C，電導(dǎo)矩陣G，均為稀疏矩陣。使用稀疏矩陣存儲這些矩陣，并且涉及這些空間離散以后稀疏矩陣非零元素比例為15%~20%采用稀疏矩陣以后稀疏矩陣存儲空間節(jié)省了70%。稀疏矩陣是一個(gè)計(jì)算方法，在本文以后章節(jié)，涉及大型稀疏矩陣，均要采用稀疏矩陣技術(shù)，用來節(jié)省存儲空間和計(jì)算量。應(yīng)用稀疏矩陣應(yīng)注意以下問題：(1)對一個(gè)稀疏矩陣進(jìn)行某些運(yùn)算，例如求逆、平方等。結(jié)果很可能不再是稀疏矩陣，應(yīng)立即轉(zhuǎn)為普通矩陣存儲格式。否則占用存儲空間不但不會減少，還會增加很多，并且由于存取數(shù)據(jù)的原因，計(jì)算速度也會大幅下降。(2)稀疏矩陣只對非零元素進(jìn)行計(jì)算，但是某些計(jì)算，零元素也會產(chǎn)生數(shù)值，例如：對稀疏矩陣T,求eT,其中零元素會產(chǎn)生數(shù)值，應(yīng)避免錯(cuò)誤。第四章頻變變壓器的線圈時(shí)域仿真基于第三章的理論，本章分別使用了精細(xì)積分法和反向差分法，在時(shí)域仿真了VFTO作用下變壓器線圈的電位分布仿真與實(shí)驗(yàn)結(jié)果接近證明了理論的正確性。然后對拉西瓦750kV晉東南1000kV變壓器進(jìn)行了VFTO作用下的電位分布計(jì)算，對其絕緣設(shè)計(jì)進(jìn)行指導(dǎo)。文中還計(jì)算了不同上升沿時(shí)間最大匝間電壓的變化規(guī)律，以找出VFTO對絕緣影響的最主要因素本章最后仿真了用氧化鋅非線性電阻做變壓器內(nèi)保護(hù)時(shí)變壓器的電壓分布，發(fā)揮了時(shí)域方法處理非線性問題的優(yōu)點(diǎn)，取得了良好的效果。4.1直求解微分方程仿真電位4.1.1精細(xì)積分法CFD把高頻變壓器模型數(shù)學(xué)形式化作微分方程組，而精細(xì)積分法[37,39](PreciseIntegrationMethod,PIM)則是用來求解這個(gè)微分方程組，由此得出變壓器線圈在VFTO下過電壓分布。式(3-7)的解析解為[39]xt)=et

x(0)+

ftf0e

A(t?ξ)P

ξ)dξ以τ為時(shí)間步長，則t時(shí)刻的解為τ

τ ξx(n+)=e

x(n)+0

e fPt?ξ)ξ

(4-1)在時(shí)間步長內(nèi)，設(shè)非齊次項(xiàng)fp是線性的或可以線性化的，即可以表示為fPt)=0+1t?tn)

(4-2)其中0=

fPtn,1=

fP;tn將式(4-2)代入式(4-1)可得迭代公式x(n+)=eτ?x(n)+τ?f

(n)?+τ?f

(n+)

(4-3)2 p 2 p該算法的關(guān)鍵在于計(jì)算eτ。令T=eτ，m=2N，N=2,Δt=τ/m，則eA?Δt≈1+A?Δt+1[A?Δt2+1[A?Δt3+1[A?Δt42 3! 4!=1+a因此2T=e2

Aτ

=?e

A?t

N2N=1+a2N

2N1[ ]=1+[ ]

2N1[ ]×1+[ ]這種分解共進(jìn)行N次。而

1+a]×1+a]=1+2a+a×a精細(xì)計(jì)算可用下列語句實(shí)現(xiàn)

for（j＝0；j<N，j++）a=2a+a×aendT=1+a這里最后一步再加單位矩陣l，能有效防止嚴(yán)重喪失有效數(shù)字。PIM的優(yōu)點(diǎn)為(1)Δt的值極小，eA?Δt近似到了10-30絕大多數(shù)計(jì)算機(jī)已經(jīng)作為舍入誤差，保證了計(jì)算的高精度；(2)eτ只需要計(jì)算一次，可以一直使用；時(shí)間步長τ可以取較大的數(shù)值，提高了計(jì)算效率。4.1.2仿真與實(shí)驗(yàn)驗(yàn)證為了驗(yàn)證所提方法的可行性和準(zhǔn)確性，首先對連續(xù)式變壓器線圈模型進(jìn)行了計(jì)算和測量。線圈模型的基本參數(shù)見下表[27]表4-1變壓器線圈基本參數(shù)參數(shù)量參數(shù)值線餅數(shù)18每餅匝數(shù)10導(dǎo)體寬度(mm)6.95導(dǎo)體高度(mm)1.2平均匝長(m)1.4828絕緣紙厚度(mm)3.00絕緣相對介電常數(shù)2.2使用本文敘述方法，VFTO下線圈數(shù)學(xué)模型在頻域?yàn)槿缦码妶?bào)方程[27]?dU(x,s) ?R=? +RR

s+sL

+sL

1?I(x,s)? dx

?0 1 0 1 s?? ? ???dI(x,s)=?(sC? dx

0+0

)U(x,s)其中，0,1,0,1,C0,0為由線圈模型單位長度的電氣參數(shù)常數(shù)矩陣，β為常量。合并同類項(xiàng)，上式可化為?dU(x,s)=?R

+sL

+?

s)I(x,s)? dx?

0 0 1(1 1 1?dI(x,s)=?(s?+G(1 1 1

)U(x,s)dx 0式中 ?＝R＋L，本文在H～0MHz頻率范圍內(nèi)采用對數(shù)平均分布的104個(gè)頻點(diǎn)使用10階實(shí)數(shù)極點(diǎn)和留數(shù)進(jìn)行矢量匹配，匹配前后均方根誤差為5×3。 s匹配后的幅頻特性曲線如圖4-1所示。圖4-1 s匹配前、后的幅頻特性曲線在仿真計(jì)算及試驗(yàn)中，線圈首端直接輸入電壓，末端直接接地。輸入電壓波形的上升沿時(shí)間為20ns，如圖4-2。計(jì)算中每根傳輸線分為5段，時(shí)間步長為5ns。圖4-3分別給出了第2、4、6、8、10、1214、16餅?zāi)┒穗妷旱姆抡嬗?jì)算波形和測量波形；表4-2為測量與仿真輸出電壓峰值的對比；表4-3為測量與仿真輸出電壓的諧振頻率及對應(yīng)的幅值由圖4-3表4-2和表4-3可見本文仿真的結(jié)果與實(shí)際測量結(jié)果基本吻合。誤差主要是由于變壓器線圈模型電氣參數(shù)的誤差所致。在采用完全傳輸線模型和相同計(jì)算（CPU為PⅢ533的條件下本文方法的仿真時(shí)間為1807s，比文獻(xiàn)[20]中頻域方法的2592s減少了30%的計(jì)算量。圖4-2仿真及實(shí)驗(yàn)中輸入電壓波形(a)第2餅輸出波形 (b)第4餅輸出波形(c)第6餅輸出波形 (d)第8餅輸出波形(e)第10餅輸出波形 (f)第12餅輸出波形(g)第14餅輸出波形 (h)第16餅輸出波形圖4-3計(jì)算和測量時(shí)域結(jié)果的比較表4-2計(jì)算和測量結(jié)果的峰值比較餅號測量(電壓/V)仿真(電壓/V)2100.0199.40468.3165.84664.6260.34856.8155.491048.4451.171246.3150.951442.5641.251640.6341.05表4-3測量與仿真結(jié)果頻譜比較餅號測量仿真諧振頻率(MHz)頻譜×10-5V/Hz諧振頻率(MHz)頻譜×10-5V/Hz20.503.280.493.640.903.550.894.411.303.521.313.241.702.981.752.152.101.812.211.3810.00.8110.040.6610.81.051.10.8512.070.5812.050.5740.605.10.585.280.974.560.915.251.312.881.313.471.701.351.750.9810.000.5510.040.491.900.4912.050.3813.200.3413.120.3360.586.280.566.170.953.140.923.291.600.711.590.782.040.732.010.769.400.639.310.7110.810.5210.900.581.200.331.230.2980.608.310.587.430.971.320.911.701.302.151.312.761.701.131.771.1410.100.5410.080.4710.800.5310.950.561.900.4712.050.47100.607.150.606.940.971.810.931.981.291.301.271.201．711.261.761.1710.000.5410.130.481.890.451.940.46120.606.150.606.430.972.880.933.101.301.431.311.411.711.291.761.2410.000.5510.090.481.900.461.940.48140.605.10.585.290.974.560.915.301.312.871.323.451.701.351.770.9810.00.5510.040.491.900.4912.070.38160.605.510.585.590.974.620.914.881.312.771.322.921.701.191.741.0510.00.5310.090.491.900.4812.030.41本文還仿真了不考慮頻變時(shí)的結(jié)果，在圖4-4中，虛線是不考慮頻變時(shí)的各餅輸出電壓，實(shí)線是考慮頻變時(shí)的各餅輸出電壓。從圖中可以看到，未考慮頻變時(shí)，電壓峰值明顯比考慮頻變時(shí)高。主要原因是，考慮頻變時(shí)，由于集膚效應(yīng)，電流集中在導(dǎo)線表層，阻抗要大的多，所以電壓波衰減比較快。(a)第2餅輸出波形 (b)第4餅輸出波形(c)第8餅輸出波形 (d)第14餅輸出波形圖4-4考慮頻變前、后計(jì)算結(jié)果的比較4.2微方程離散為代數(shù)方程仿真電位由于內(nèi)屏蔽式線圈和連續(xù)式線圈邊界條件的不同，導(dǎo)致PIM的數(shù)值精度變差，因此引入了基于反向差分代替微分代數(shù)方程仿真內(nèi)屏蔽式變壓器線圈的方法。由于內(nèi)屏蔽式變壓器線圈中屏蔽線末端都是懸空的，所以屏蔽線末端邊界條件是Ib=0(pb代表屏蔽線末端的編號)。從第三章方程(3-6)可以看到，對于懸空屏蔽線附加向量f中的電壓邊界條件是未知的；如果在方程(3-6)中添加懸空線電流邊界條件，P矩陣不可逆就無法使用PIM求解狀態(tài)方程為了解決上述問題本節(jié)把方程(3-6)首先改寫為另一種形式然后使用反向差分法(ackwardDierentiationForulas,BDF)進(jìn)行仿真。由于本方法不再對邊界條件限制，故可以應(yīng)用到任何形式的線圈。4.2.1反向差分法(1)CFD方程的改寫方程組(3-3)、(3-4)、(3-5)可以改寫做如下形式Pdx+Qx+f=0dt方程(4-1)中各矩陣和向量的意義和(3-6)中已經(jīng)不同。其中

(4-4)?0 3 1 0 ?? ?? 1

a2 % ?? %%% ?? ?? % a2 1 0 ?? 1 3 ?P=? ??# # 4 1 ?? 0 1

2 % ?? ?? % % % ?? ???0 0

% 21

1?4? 3 1

1 1 ?? ?? 1

2 % %% 0 ?? %%% 0

%% ?? ?? %2 1

%% ?? 1 3

1 1?? ?Q=??1 1 0

4 1 ?? %%

1 2 % ?? ?? %% %%% ?? %% %%% ?? ?? 0 %% %2

1?? ?? ?11

1 4?x=?u0f=[0

12"

320],

" uM1/2 uM

I0 1

T" IM ,矩陣PQ的維數(shù)是(2M+)N×(2M+3)N,也就是方程數(shù)目為(2M+)N未知數(shù)數(shù)目為(2M+3)N。注意方程(4-4)中并未添加邊界條件，因此附加向量f為零向量。對于首尾相連的線匝，有如下邊界條件ijM

=(j+)0，ujM

=(j+)0(j=,2???N?)對于末端懸空的線匝，邊界條件為添加邊界條件有兩種方法：

ipb=0。(a)使用類似第三章中處理連續(xù)線圈邊界條件的方法進(jìn)行矩陣變換添加首尾相連的線匝的邊界條件；對末端懸空的線匝，寫做新的方程，列入方程組。(b)把首尾相連的線和末端懸空線均寫做新的方程，列入方程組。無論采用那種方法，都要把變壓器線圈首、末端邊界條件列方程寫入方程組。第一種方法可以節(jié)省存儲空間和計(jì)算量。本文采用第一種方法進(jìn)行仿真。添加邊界條件以后方程和未知數(shù)個(gè)數(shù)均為(2M+)N+2(Npb+)其中Npb為屏蔽線的懸空端數(shù)目。(2)化微分方程為代數(shù)方程[39]對于方程(4-4)，添加邊界條件以后，微分項(xiàng)的系數(shù)矩陣是不可逆的，因此，需要化作代數(shù)方程求解，以避免求逆問題。對于形式如下的方程稱為微分代數(shù)方程f(x,,t)=0； (4-5)其中未知數(shù)個(gè)數(shù)和方程個(gè)數(shù)相等。相對傳統(tǒng)所述顯式狀態(tài)方程=h(x,t)

(4-6)避免了把方程首先轉(zhuǎn)化為(4-6)式的形式，這樣就可以避免直接求逆。假設(shè)方程(4-5)在t=tn，t=tn1，…，t=tn?k時(shí)的解x(t)均已經(jīng)求得，其中步長i=ti1?ti不一定一致。使用i表示x(i)，則方程(4-5)在t=tn1時(shí)的解n1必須滿足f(n1,(n+),tn1)=0

(4-7)使用反向差分公式,在任何給定精度范圍在t=tn1時(shí)，用n1和k個(gè)過去值n,n1,",n?k1表示時(shí)間導(dǎo)數(shù)的現(xiàn)在值(tn1)如果設(shè)n1表示(tn1)的近似值則k階BDF表示為1k1n1=? ∑αixn?ig(n1)hi=0

(4-8)其中α0,α1,",αk為常數(shù)h為時(shí)間步長式(4-8)中只有n1是未知的因此表示為其函數(shù)g(n1)。將(4-8)代入式(4-7)f(n1,g(n1),tn1)=0

(4-9)上式成為以未知變量n1表示的一個(gè)代數(shù)方程組，可以是線性或者非線性的，本文中高頻變壓器線圈模型的方程是線性的，但是本文方法也可以計(jì)算非線性模型。BDF的階數(shù)k取值范圍為3-6k取6的時(shí)候精度已經(jīng)非常高了足以滿足絕大多數(shù)問題的需要。本文仿真過程階數(shù)可以設(shè)定，仿真結(jié)果證明，采用4階精度時(shí)已經(jīng)可以滿足精度要求當(dāng)變步長時(shí)BDF系數(shù)是與步長有關(guān)的一個(gè)方程組本文中步長是固定的，式(4-8)中系數(shù)α可由如下公式求解?1 1 1 1 "

1??α0?

?0?? ?? ? ???01 2 3 "

k??α1?

?1?2?01 222

2 "

k2??α?

?0?? ??

?=??

(4-10)?# # # # " #??#??# # # # " #??#?

?#??#?? ?? ? ??01 2k

k " kk α 0? ??k? ??上式左邊系數(shù)矩陣是范德蒙德(andermonde)矩陣定步長情況下α0,α1,",αk的值與離散時(shí)間和步長無關(guān)。(3)反向差分表示反向差分公式使用以上算法，執(zhí)行(4-8)式時(shí)，必須存儲前k步的值n,n1,",n?k1，文獻(xiàn)[40]證明存儲反向差分有利于減少舍入誤差，令則式(4-8)化為

n?in?i?n?i,i=,,",k1k1

(4-1)n1=?

∑?iΔxn?i

(4-12)hi=0式中式(4-8)化為

i?i∑αij=0kpn1=n+∑iΔxn?ip

(4-13)i1式中ii1+∑γjj1文獻(xiàn)[39]中證明了反向差分表示反向差分公式是等價(jià)的。PAGEPAGE284.2.2BDF應(yīng)用于變壓器線圈本章4.1中得到的方程組，為了敘述簡便沒有添加頻變因素，現(xiàn)在使用矢量匹和配遞歸卷積，可以得到考慮頻變時(shí)變壓器線圈的數(shù)學(xué)模型。頻變的變壓器線圈可表示為??+??+?

+?=0

(4-14)t 00其中?，?為相關(guān)系數(shù)矩陣，?為未知電壓電流列向量，?為邊界條件，?為頻變0產(chǎn)生的附加項(xiàng)，式中各項(xiàng)和第三章中相比，意義相同，但是內(nèi)容已經(jīng)不同。應(yīng)用改進(jìn)后的BDF使用式(4-12)逼近式(4-14)中的導(dǎo)數(shù)項(xiàng)化簡后在時(shí)刻n+1,可以表示為?(

??a???

n1+1)

+??

n1+f0

+?=0

(4-15)0 1上式中，a=α/τ，?為前k部值的加權(quán)項(xiàng)，其余各相意義和式(4-14)中相同。0 1對式(4-15)進(jìn)一步化簡? ?(a?+?)?n1=?(??1+f0+?)? ?

(4-16)顯然，方程(4-16)的解為

? 1

?? ? ?n1=?(aP+?)

?(P?1+f0+λ)

(4-17)對于方程(4-16)，如果階數(shù)比較低，可以直接求逆。如果階數(shù)比較高，本文使用了LU分解法[36]來求解。方程(4-16)可以簡寫為代數(shù)方程形式Ax=bL、U分解后有

A=LU其中，L為下三角矩陣，U為上三角矩陣。則原方程可以化為令則原方程化為

LUx=bx=y

(4-18)Ly=b (4-19)使用追趕三角法求解方程(4-19)得到y(tǒng)值，再使用一次追趕三角法求解方程(4-18)得到x值，即解得了原方程的解。本文使用LU分解法主要原因在于，分解后可以保證矩陣L、U的稀疏性，使用稀疏矩陣技術(shù)，節(jié)省計(jì)算量和存儲空間。4.2.3仿真驗(yàn)證與實(shí)例計(jì)算(1)仿真驗(yàn)證由于本節(jié)的算法同時(shí)可以求解本章第一節(jié)試驗(yàn)中的連續(xù)式線圈，所以首先用上一節(jié)的實(shí)驗(yàn)檢驗(yàn)本節(jié)算法的正確性。仿真與測量結(jié)果比較見圖4-5。由圖可以看到仿真與測量結(jié)果在幅值與波形都很接近，造成誤差的主要原因是變壓器線圈的復(fù)雜性和建立模型做了一些簡化。相對本章第一節(jié)的方法，本節(jié)方法使矩陣維數(shù)增大，計(jì)算機(jī)資源占有增多。圖4-5仿真與測量結(jié)果比較(2)實(shí)例驗(yàn)證天威集團(tuán)使用變壓器線圈沖擊電壓計(jì)算軟件仿真雷電波下線圈過電壓分布經(jīng)過與實(shí)測比較，該軟件結(jié)果與實(shí)際接近，精度足夠指導(dǎo)實(shí)際生產(chǎn)。該軟件對每餅線圈使用集中參數(shù)模型進(jìn)行建模無法仿真VFTO高頻信號作用下的線圈中過電壓分布應(yīng)用本文算法計(jì)算拉西瓦750kV主變雷電波作用下線圈過電壓分布和已有軟件結(jié)果比較，以證明本文算法的正確性。輸入波為雙指數(shù)函數(shù)疊加成的標(biāo)準(zhǔn)實(shí)驗(yàn)雷電波，波前時(shí)間為1.2μs，波尾時(shí)間為50μs，峰值為100kV。即??t

?t?()ut=A?eτ1?eτ2()? ?? ?

()其中，τ1為波尾時(shí)間，τ2為波前時(shí)間，A=104kV.本文算法和軟件得到的一些餅的末端對地電壓如表4-4所示表4-4軟件和本文方法各餅峰值比較餅號軟件電壓峰值(電壓/V)本文方法電壓峰值(電壓/V)兩者誤差百分比293.5292.101.5184％487.6086.121.6895％1071.0782.2115.6747％1871.6073.192.2207％2469.8361.921.3275％3458.3655.435.0206％4433.6635.164.4563％5414.7114.014.7587％表4-4給出了前幾餅和中間兩者誤差比較大的一些餅的末端電位從表中數(shù)據(jù)可以看出在VFO最容易造成損害的前幾餅兩者非常接近中間的一些餅誤差比較大，但是O一般不會對這些餅造成損害。誤差主要原因是本文算法未考慮副邊的影響。在O作用時(shí)間比較短副邊影響小雷電波作用時(shí)間比較長副邊影響就會造成比較大的誤差。(3)上升沿時(shí)間對匝間電壓的影響由文獻(xiàn)[5]VFTO頻率對變壓器的匝間絕緣影響相當(dāng)大本文在不改變輸入波幅值的情況下，僅改變電壓波的頻率。當(dāng)頻率提高時(shí)候，上升沿時(shí)間相應(yīng)變短。以拉西瓦750kV主變?yōu)槔幼饔迷谧儔浩魅肟谔幍腣FTO電壓保持幅值不變波前時(shí)間分別為20ns,60ns,100ns.圖4-6中給出了不同波前時(shí)間情況下最大匝間電壓的波形最大匝間電壓峰值的絕對值分別為19.2kV18.7kV16.1kV由圖可以看到，保持幅值不變上升沿時(shí)間相應(yīng)變短，匝間電壓絕對值越大，對匝間絕緣影響越大。因此頻率越高，對匝間絕緣的影響也越大。圖4-6波前時(shí)間為20ns、60ns、100ns時(shí)最大匝間電壓的波形(4)實(shí)例計(jì)算本文對設(shè)計(jì)中的晉東南1000kV主變高壓線圈VFTO下電壓分布進(jìn)行了計(jì)算。該變壓器高壓線圈屏蔽段示意圖如圖4-7所示?？傇褦?shù)為842。按照設(shè)計(jì)要求，GIS開關(guān)投切時(shí)可能產(chǎn)生波頭很陡的暫態(tài)過電壓(VFTO)(暫按2.5倍設(shè)備最高運(yùn)行電壓考慮)，因此施加在變壓器高壓繞組首端的暫態(tài)過電壓最大幅值按2.5*100/3＝1588kV考慮本文輸入波的波前時(shí)間為20ns是情況比較嚴(yán)重的時(shí)候，如圖4-8(a)所示。圖4-8(b)-4-8(m)給出了變壓器線圈匝間電壓，可以看出最高匝間電壓出現(xiàn)在第1餅1號線與屏蔽線之間峰值為384.8kV從圖中還可以看出電壓差比較大的匝主要出現(xiàn)在線圈前幾匝,也就是電壓降主要在前幾匝上是絕緣保護(hù)的重點(diǎn)而以后的線匝壓降就比較低了。圖4-8(n)-4-8(r)給出了變壓器前幾匝線圈對地電壓從圖中可以看出電壓變化比較大的匝主要出現(xiàn)在線圈前幾匝，這也是前幾匝線圈匝間電壓比較大的原因。圖4-8(s)-4-8(v)給出了變壓器線圈前幾餅的餅間電壓從圖中可以看出電壓差比較大的餅主要出現(xiàn)在線圈前幾餅，峰值為445.3kV。圖4-7 高壓線圈屏蔽段結(jié)構(gòu)示意圖(a)線圈入口處VFO電壓波形波前時(shí)間為2ns,峰值為5kV

(b)最大匝間電壓峰值為348k,出現(xiàn)在第1線圈與屏蔽線之間(c)第1餅2號線與其前一匝屏蔽線的電壓差，峰值絕對值為224kV

(d)第1餅2號線與其后面屏蔽線的電壓差，峰值絕對值為2.5kV(e)第1餅3號線與其前面屏蔽線的電壓差，峰值絕對值為2.3kV

(f)第1餅3號線與其后面屏蔽線的電壓差，峰值絕對值為1.3kV(g)第1餅4號線與其前面屏蔽線的電壓差，峰值絕對值為1.1kV

(h)第1餅6號工作線與其前面屏蔽線的電壓差，峰值絕對值為103kV(i)第1餅6號工作線與7號工作線的電壓差，峰值絕對值為83.1kV

(j)第2餅8號線與其前面屏蔽線的電壓差，峰值絕對值為78kV(k)第3餅6號工作線與7號工作線的電壓差，峰值絕對值為88.4kV

(l)第5餅7號工作線與8號工作線的電壓差，峰值絕對值為26.1kV()第1餅3號工作與4號工線的電壓差，峰值絕對值為1.kV

(n)第1第1匝屏線的電壓電壓峰值為186kV(o)第1餅2號工作線的電壓電壓峰值為18.V

(p)第1第2匝屏線的電壓電壓峰值為148kV(q)第1餅6號工作線的電壓電壓峰值為16.V

(r)第2餅2號工作線電壓電壓峰值為10.V(s)第1餅第2餅最電壓值，位置在第1匝作線，峰值為453kV

(t)第2餅第3餅最電壓值，位置在3餅8工作線，峰值為350kV(u)第34餅間最大電壓值位置在4餅134工作線，峰值為319kV

(v)第45餅間最大電壓值位置在5餅5號工線，峰值為787kV圖4-8變壓器入口電壓、匝間電壓、對地電壓、餅間電壓4.3變器線圈餅間并聯(lián)MOV限制VFO幅值文獻(xiàn)[41][42]分別通過實(shí)驗(yàn)和計(jì)算研究了在變壓器每餅上跨接金屬氧化物可變電阻(MentalOxidearisto,MOV)來限制變壓器內(nèi)部的電壓諧振，取得了令人滿意的效果。但是文中把每餅線圈用集中參數(shù)元件建模，這種模型只適用于雷電波侵入的情況，無法滿足研究VFTO作用下的電壓分布的需要。本文使用多導(dǎo)體傳輸線對變壓器線圈建模，在拉西瓦水電站750kV主變高壓線圈餅間跨接MOV。仿真VFTO作用下的線圈電位分布因?yàn)镸OV為非線性元件頻域方法無法處理而本文使用的時(shí)域方法可以很容易處理含非線性元件的問題。本文中的MOV是氧化鋅非線性電阻，MOV的U-I特性[43]可以用下式表示U=Iα式中α為非線性系數(shù)，取值為0.01-0.05。其中k=CH/Sα

(4-21)H為厚度，S為截面積，C為材料常數(shù)。也就是k與幾何尺寸，材料常數(shù)，非線性系數(shù)有關(guān)。MOV做變壓器線圈內(nèi)保護(hù)原理可以用如圖4-9的U-I特性圖說明當(dāng)作用在MOV兩端的電壓差小于閥值時(shí)，MOV工作在小電流區(qū)，如圖4-9所示的1區(qū)，流過MOV的電流小于1mA,這時(shí)MOV呈現(xiàn)很大的電阻，幾百兆歐級，連接MOV的支路近似于開路。在過電壓作用下，當(dāng)作用在MOV兩端的電壓差超過閥值時(shí)，MOV的工作點(diǎn)要發(fā)生躍變，進(jìn)入大電流區(qū)．如圖4-9中2區(qū)所示．在大電流區(qū)，MOV呈現(xiàn)很小的電阻，流過電流可達(dá)百安或千安級，并伴隨吸收過電壓能量，實(shí)際上其電壓變化不大這樣對VFTO過電壓來說。MOV可以將過電壓限制在一定的范圍內(nèi)，對共振過電壓．MOV能起阻尼作用，抑制局部共振的幅值。當(dāng)電流持續(xù)增大，MOV過流能力達(dá)到飽和時(shí)，電流值隨電壓變化比2區(qū)顯著變小，如圖4-9所示3區(qū)，進(jìn)入3區(qū)臨界點(diǎn)電流大約是10kA左右。U1 2 3I圖4-9 MOV的U-I特性曲線圖4-10給出了拉西瓦水電站主變高壓線圈前幾餅線匝分布簡圖以及MOV的連接方式.該變壓器線圈詳細(xì)結(jié)構(gòu)類似圖4-7中所示，也是內(nèi)屏蔽式，為了清楚明了，圖4-10中略去了屏蔽線文獻(xiàn)[44]使用EMTP仿真了操作過電壓,在變壓器入口處的波形如圖4-1所示圖4-10變壓器線圈中MOV連接簡圖圖4-1變壓器入口處操作過電壓本文使用公式(4-21)構(gòu)造了一個(gè)壓控非線性電阻(MOV)，MOV閥值電壓均為1.5kV,非線性系數(shù)取0.03。在狀態(tài)方程(4-14)中添加MOV的VCR方程，代入接MOV的線匝端部的KCL方程這樣就得到了接MOV時(shí)的狀態(tài)方程.新狀態(tài)方程未知變量中增加了三個(gè)MOV的電流為未知量。計(jì)算過程中，在未知量里把MOV電流寫在一起，以便使用矩陣分塊計(jì)算，加快計(jì)算速度。圖4-12給出了變壓器第10匝入口處對地電壓,該處是由諧振引起的對地電壓最高的線匝,由圖可見,MOV抑制諧振電壓效果明顯.由圖4-13可以看到餅間電壓差比較大主要在前四餅之間,接MOV前后相比餅間電壓有了明顯的降低電壓峰值限制在了MOV閥值電壓內(nèi).圖4-14所示為各MOV中的電流曲線,電流峰值為2.1kA，出現(xiàn)在編號為1的非線性電阻處.圖4-12加MOV前后第10匝入口端電壓圖4-13加避雷器前后1-2、2-3、3-4餅間電壓圖4-14避雷器上的電流第五章基于krylov子空間的時(shí)域降階算法求解大規(guī)模矩陣問題包括的線性方程組求解是科學(xué)工程計(jì)算中的重大課題。本文涉及了高頻下的變壓器線圈模型，按照第四章的方法，多導(dǎo)體耦合線模型最后形成一個(gè)高階的代數(shù)方程組最高達(dá)到七千多階計(jì)算量非常龐大幾乎已經(jīng)到了PC機(jī)的處理極限。因此引入了階數(shù)縮減技術(shù)(OrderReduction)[45-48]，來節(jié)省存儲空間和降低計(jì)算量。最近幾年,數(shù)學(xué)領(lǐng)域關(guān)于階數(shù)縮減方面的研究工作取得了許多重大進(jìn)展，許多成果已經(jīng)應(yīng)用到了工程方面，例如微電子學(xué)科。在過去一段時(shí)間里，電子和電力學(xué)科研究中也使用了很多階數(shù)縮減的方法。第一類是漸進(jìn)波形估計(jì)法(WE)又稱Pade逼近[45]來實(shí)現(xiàn)對大型線性集中網(wǎng)絡(luò)的階數(shù)縮減原理是尋求一種形式簡單而特性又和原網(wǎng)絡(luò)函數(shù)逼近的函數(shù)。這種方法適用范圍廣，但是，它不能保證階數(shù)縮減模型的無源性第二類是基于Krylov子空間的階數(shù)縮減方法[46-50]這類方法是使用低階矩陣的特征值來等效高階矩陣的特征值可以保證系統(tǒng)的無源性例如Arnoldi,Lanczos算法。本文使用Arnoldi,Lanczos算法來處理本文的高階模型。求解線性方程組問題的Krylov子空間方法可以追溯到上世紀(jì)50年代初，此后Lanczos,Arnoldi,Paige,Saad[51-52]等人做了許多的工作。人們做了大量的理論分析和數(shù)值試驗(yàn),充分認(rèn)識到Krylov子空間方法是求解大型線性方程組和大型矩陣特征值問題的一類最有效的方法。5.1基于krylov子空間的代數(shù)方程的降階5.1.1空間簡述為便于理解，首先簡述子空間映射在解大型線性方程中的應(yīng)用。設(shè)一個(gè)線性方程組以下面的矩陣方程表示式中A∈Rn×n

為系數(shù)矩陣、x∈Rn

Ax=b (5-1)為待求變量列向量、b∈Rn為常系數(shù)列向量。令m是n維空間中m個(gè)線性無關(guān)的n維向量m=[1,v2,"vm]×mm∈RKm=span(1,v2,"vm)則稱Km是以上述m個(gè)n維向量作為基底在n維空間中展成的m維子空間

[24]利用向子空間投影的概念可求上述方程(5-1)的近似解。設(shè)x(m)表示利用子空間投影的近似解，此時(shí)可得到如下的等式m其中X(m)∈K，因此可以得到m

x(m)=b (5-2)令 (m) (m)

(Ax

(m)

?b)

⊥vj

j=,,",m

(5-3)x =mym其中V為上述子空間的一個(gè)基，式中y(m)為一個(gè)m維向量。m則可將方程(5-3)改寫為mVT(m

my

(m)

?b)