三相變壓器建模及仿真及MATLAB仿真

1、三相變壓器建模及仿真及MATLAB仿真XXXXXXX 學(xué)院課程設(shè)計(jì)報(bào)告課程名稱: 系 部: 專業(yè)班級(jí): 學(xué)生姓名: 指導(dǎo)教師: 完成時(shí)間: 報(bào)告成績:評(píng)閱意見:學(xué)疣救學(xué)工作部制目錄摘要第一章變壓器介紹41.1 蚯器的磁化雌 41.2 變壓器保護(hù)41.3 勵(lì)筋甬流 7第二章變壓器基本原理92.1 變壓器工作原理 92.2 三相變壓器的等效電路及聯(lián)結(jié)組 10第三章變壓器仿真的方法 111. 1基于基本勵(lì)磁曲線的靜態(tài)模型 113. 2基于暫態(tài)磁化特性曲線的動(dòng)態(tài)模型 134. 3非線性時(shí)域等效電路模型 14第四章三相變壓器的仿真 165. 1三相變壓器仿真的數(shù)學(xué)模型 166. 2電源電壓的描述 207

2、. 3鐵心動(dòng)態(tài)磁化過程簡述 21第五章變壓器MATLAB仿真研究255.1仿真長線路末端電壓升高 255. 2仿真三相變壓器T2的勵(lì)磁涌流 288. 3三相變壓器仿真模型圖 345.4變壓器仿真波形分析 36結(jié)論 40參考文獻(xiàn)41摘要在電力變壓器差動(dòng)保護(hù)中，勵(lì)磁涌流和內(nèi)部故障電流的判別一直是一個(gè)關(guān) 鍵問題。文章闡述了勵(lì)磁涌流的產(chǎn)生及其特性，利用MATLAB對(duì)變壓器的勵(lì)磁 涌流、內(nèi)部故障和外部故障進(jìn)行仿真，對(duì)實(shí)驗(yàn)的數(shù)據(jù)波形分析，以此來區(qū)分故 障和涌流，目的是減少空載合閘產(chǎn)生的勵(lì)磁涌流對(duì)變壓器差動(dòng)保護(hù)的影響，提 高保護(hù)的靈敏性。本文在Matlab的編程環(huán)境下，分析了當(dāng)前的變壓器仿真的方法。在單相情

3、 況下，分析了在飽和和不飽和的勵(lì)磁涌流現(xiàn)象，和單相勵(lì)磁涌流的特征。在三 相情況下，在用分段擬和加曲線壓縮法的基礎(chǔ)上，分別用兩條修正的反正切函 數(shù)，和兩條修正的反正切函數(shù)加上兩段模擬飽和情況的直線兩種方法建立了 Ydll> YndlK YnyO和YyO四種最常用接線方式下三相變壓器的數(shù)學(xué)仿真模型， 并在Matlab下仿真實(shí)現(xiàn)。通過對(duì)三相勵(lì)磁涌流和磁滯回環(huán)波形分析，三相勵(lì)磁 涌流的特征分析，總結(jié)出影響三相變壓器勵(lì)磁涌流地主要因素。最后，分析了 兩種方法的優(yōu)劣，建立比較完善的變壓器仿真模型。關(guān)鍵字：變壓器；差動(dòng)保護(hù)；勵(lì)磁涌流；內(nèi)部故障；外部故障；波形分析；仿真; 數(shù)學(xué)模型-6 -第一章變壓器介

4、紹1.1 變壓器的磁化特性初始磁化曲線當(dāng)電流從0逐漸增加，線圈中的磁場強(qiáng)度H也隨之增加，這樣就可 以測(cè)出若干組B,H值。以H為橫坐標(biāo)，B為縱坐標(biāo)，畫出B隨H的變 化曲線，這條曲線稱為初始磁化曲線。當(dāng)H增大到某一值后，B幾乎不 再變化，這時(shí)鐵磁材料的磁化狀態(tài)為磁飽和狀態(tài)。此時(shí)的磁感應(yīng)強(qiáng)度Bs 叫做飽和磁感應(yīng)強(qiáng)度。這種磁化曲線一般如下圖中曲線所示：OH1.2 變壓器保護(hù)電力變壓器是電力系統(tǒng)中大量使用的重要電氣設(shè)備，他的故障給供電可靠 性和系統(tǒng)的正常運(yùn)行帶來嚴(yán)重的后果，同時(shí)大容量變壓器也是非常貴重的元 件，因此，必須根據(jù)變壓器的容量和重要程度裝設(shè)性能良好的、動(dòng)作可靠的保 護(hù)元件。電力變壓器的故障分為

5、內(nèi)部和外部兩種故障。內(nèi)部故障指變壓器油箱 里面發(fā)生的各種故障，主要靠瓦斯和差動(dòng)保護(hù)動(dòng)作切除變壓器；外部故障指油 箱外部絕緣套管及其引出線上發(fā)生的各種故障，一般情況下由差動(dòng)保護(hù)動(dòng)作切 除變壓器。速動(dòng)保護(hù)（瓦斯和差動(dòng)）無延時(shí)動(dòng)作切除故障變壓器，設(shè)備是否損 壞主要取決于變壓器的動(dòng)穩(wěn)定性。而在變壓器各側(cè) 母線及其相連間隔的引出設(shè) 備故障時(shí)，若故障設(shè)備未配保護(hù)（如低壓側(cè)母線保護(hù)）或保護(hù)拒動(dòng)時(shí)，則只能 靠變壓器后備保護(hù)動(dòng)作跳開相應(yīng)開關(guān)使變壓器脫離故障。因后備保護(hù)帶延時(shí)動(dòng) 作，所以變壓器必然要承受一定時(shí)間段內(nèi)的區(qū)外故障造成的過電流，在此時(shí)間 段內(nèi)變 壓器是否損壞主要取決于變壓器的熱穩(wěn)定性。因此，變壓器后備保

6、護(hù)的 定值整定與變壓器自身的熱穩(wěn)定要求之間存在著必然的聯(lián)系。1）瓦斯保護(hù)對(duì)變壓器油箱內(nèi)部的各種故障及油面的降低應(yīng)裝設(shè)瓦斯保護(hù)。容量為 800KVA及以上的油浸式變壓器，對(duì)于容量為400KVA及以上的車間內(nèi)油浸式變 壓器，勻應(yīng)裝設(shè)瓦斯保護(hù)。當(dāng)油箱內(nèi)部故障產(chǎn)生輕微瓦斯或油面下降時(shí)，；保 護(hù)裝置應(yīng)瞬間動(dòng)作于信號(hào)：當(dāng)產(chǎn)生大量瓦斯時(shí)，瓦斯保護(hù)宜動(dòng)作于斷開變壓器 各電源側(cè)斷路器。對(duì)于高壓側(cè)未裝設(shè)斷路器的線路變壓器組，未采取使瓦斯保 護(hù)能切除變壓器內(nèi)部故障的技術(shù)措施時(shí)瓦斯保護(hù)可僅動(dòng)作與信號(hào)。2）縱差保護(hù)或電流速斷保護(hù)容量在10000KVA及以上的變壓器應(yīng)裝設(shè)縱差保護(hù)，用以反應(yīng)變壓器 內(nèi)部繞組、絕緣套管及引出

7、線相間短路、中性點(diǎn)直接接地電網(wǎng)側(cè)繞組和 引出線的接地短路以及繞組匝間短路。3）過流保護(hù)變壓器的過流保護(hù)用作外部短路及變壓器內(nèi)部短路的后備保護(hù)。4）零序過流保護(hù)變壓器中性點(diǎn)直接接地或經(jīng)放電間隙接地時(shí)，應(yīng)補(bǔ)充裝設(shè)零序過流 保護(hù)。用以提高保護(hù)在單相接地時(shí)的靈敏度。零序過流保護(hù)主要用作 外部電網(wǎng)接地短路的后備保護(hù)。5）過負(fù)荷保護(hù)變壓器過負(fù)荷時(shí)，應(yīng)利用過負(fù)荷保護(hù)發(fā)出信號(hào)，在無人值班的變電 所內(nèi)可將其作用于跳閘或自動(dòng)切除一部分負(fù)荷。靈敏度高、結(jié)構(gòu)簡單，并能反應(yīng)變壓器油面內(nèi)部各種類型的故障。 特別是當(dāng)繞組短路匝數(shù)很少時(shí)，故障點(diǎn)的循環(huán)電流雖然很大，可能造成 嚴(yán)重的過熱，但反應(yīng)在外部電流的變化卻很小，各種反應(yīng)電

8、流量的保護(hù) 都難以動(dòng)作，因此瓦斯保護(hù)對(duì)保護(hù)這種故障有特殊的優(yōu)越性。7）縱聯(lián)差動(dòng)保護(hù)差動(dòng)保護(hù)是一種依據(jù)被保護(hù)電氣設(shè)備進(jìn)出線兩端電流差值的變化構(gòu)成的對(duì) 電氣設(shè)備的保護(hù)裝置，一般分為縱聯(lián)差動(dòng)保護(hù)和橫聯(lián)差動(dòng)保護(hù)。變壓器的差動(dòng) 保護(hù)屬縱聯(lián)差動(dòng)保護(hù)，橫聯(lián)差動(dòng)保護(hù)則常用于變電所母線等設(shè)備的保護(hù)。縱差動(dòng)保護(hù)是變壓器的電氣主保護(hù)，由于變壓器在電力系統(tǒng)中占有重要 地位，縱差動(dòng)保護(hù)必須滿足如下要求：(1)能反應(yīng)保護(hù)區(qū)內(nèi)各種相間和接地短路故障。(2) 動(dòng)作速度快，一般動(dòng)作時(shí)間不能大于30mso(3)在變壓器空載合閘或外部故障切除后電壓恢復(fù)期間產(chǎn)生勵(lì)磁涌流時(shí)不應(yīng) 誤動(dòng)作。(4)在變壓器過勵(lì)磁時(shí)，縱差動(dòng)保護(hù)不應(yīng)該動(dòng)作。(

9、5)發(fā)生外部故障時(shí)電流互感器飽和應(yīng)可靠不動(dòng)作。(6)保護(hù)區(qū)內(nèi)故障時(shí)，電流互感器飽和，縱差動(dòng)保護(hù)不應(yīng)拒動(dòng)或延時(shí)動(dòng)作。(7)保護(hù)區(qū)內(nèi)發(fā)生短路故障，在短路電流中含有諧波分量時(shí)，縱差動(dòng)保 護(hù)不應(yīng)拒動(dòng)或延時(shí)動(dòng)作。變壓器縱差保護(hù)的原理要求變壓器在正常運(yùn)行和縱差保護(hù)區(qū)(縱差保 護(hù)區(qū)為電流互感器TAK TA2之間的范圍)外故障時(shí)，流入差動(dòng)繼電器中的電流為零，保 證縱差保護(hù)不動(dòng)作。但由于變壓器高壓側(cè)和低壓側(cè)的額定電流不同，因 此，為了保證縱差保護(hù)的正確工作，就須適當(dāng)選擇兩側(cè)電流互感器的變 比，使得正常運(yùn)行和外部故障時(shí)，兩個(gè)電流相等。差動(dòng) 保護(hù)的原理接線圖：霞品嬲蹴瞬贏瓣流分布減小縱聯(lián)差動(dòng)保護(hù)的不平衡電流的措施:

10、1）保證電流互感器在外部最大短路電流流過時(shí)能滿足10%誤差曲線的要求。2）減小電流互感器二次回路負(fù)載阻抗以降低穩(wěn)態(tài)不平衡電流。3）可在差流回路中接入具有速飽和特性的中間變流器以降低暫態(tài)不平衡電流。 為保護(hù)縱聯(lián)差動(dòng)保護(hù)的選擇性，差動(dòng)保護(hù)的動(dòng)作電流必須躲開可能出現(xiàn)的最大 不平衡電流。而變壓器的勵(lì)磁電流是縱差動(dòng)保護(hù)不平衡電流產(chǎn)生原因之一, 特別是空載合閘時(shí)產(chǎn)生的很大的勵(lì)磁涌流會(huì)嚴(yán)重影響保護(hù)的靈敏性。1.3勵(lì)磁涌流勵(lì)磁涌流產(chǎn)生的機(jī)理變壓器是基于電磁感應(yīng)原理的一種靜止元件。在電能-磁能-電能能量的轉(zhuǎn) 換過程中，它必須首先建立一定的磁場，而在建立磁場的過程中，變壓器繞組 中就會(huì)產(chǎn)生一定的勵(lì)磁 電流。當(dāng)空載

11、變壓器稔態(tài)運(yùn)行時(shí)，勵(lì)磁電流很小，僅為 額定電流的0.35%10%。但當(dāng)變壓器空載合閘時(shí)，由于變壓器鐵芯剩磁的影響 以及合閘初相角的隨機(jī)性會(huì)使鐵芯磁通趨于飽和，從而產(chǎn)生幅值很大的勵(lì)磁涌 流。當(dāng)變壓器在電壓過零點(diǎn)合閘時(shí)，由于鐵芯中磁通最大，鐵芯嚴(yán)重飽和，因 此產(chǎn)生最大的勵(lì)磁電流，其峰值最大可達(dá)額定電流的68倍。如果在合閘瞬 間，電壓正好達(dá)到最大值時(shí)，則磁通的瞬間值正好為零，即在鐵芯里一開始就 建立了穩(wěn)態(tài)磁通。在這種情況下，變 壓器不會(huì)產(chǎn)生勵(lì)磁涌流。勵(lì)磁涌流的特點(diǎn)D勵(lì)磁涌流往往含有大量高次諧波分量（以二次諧波為主），使涌流波形偏于 時(shí)間軸的一側(cè)，波形含有間斷角為j-2）勵(lì)磁涌流的衰減常數(shù)與鐵芯的飽和

12、程度有關(guān)，飽和越深，電抗越小，衰減越 快。因此，在開始瞬間衰減很快，以后逐漸減慢，經(jīng)0.51s后其值不超過 （0.250.5） In,3）變壓器的容量越大，涌流的幅度越大，持續(xù)的時(shí)間越長。對(duì)于容量小的變壓 器衰減得快，約幾個(gè)周波即達(dá)到穩(wěn)定，大型變壓器衰減得慢，全部衰減持續(xù)時(shí) 間可達(dá)幾十秒。勵(lì)磁涌流的危害空載合閘產(chǎn)生的很大的勵(lì)磁涌流可能會(huì)引起繼電保護(hù)裝置的誤動(dòng)作，誘發(fā)操作 過電壓，損壞電氣設(shè)備，造成電網(wǎng)電壓和頻率的波動(dòng)；勵(lì)磁涌流包含的大量諧 波也會(huì)對(duì)電能質(zhì)量造成嚴(yán)重的污染。因此對(duì)變壓器勵(lì)磁涌流的仿真有著重要的 意義。第二章變壓器基本原理2.1 變壓器工作原理變壓器是一種靜止的電器，用于將一種形式

13、的交流電能改變成另一種形式 的交流電能，其形式的改變是多種多樣的。既可以改變電壓、電流；也可以改 變等效阻抗或電源相數(shù)、頻率等。以單相為例，研究變壓器臺(tái)變壓器的示意圖。它由鐵芯和線圈組成。接電源 的原邊線圈成為初級(jí)線圈；接負(fù)載的副邊線圈稱次級(jí)線圈。設(shè)原、副邊線圈匝 數(shù)分別為例、°?。根據(jù)電磁感應(yīng)現(xiàn)象，電能可從原邊輸送到副邊，但原、副邊 具有不同的電壓和電流。變壓器內(nèi)部的磁場分布的情況是非常復(fù)雜的，但是我們總可以把它們折算為 等效的兩部分磁通。其中一部分磁通。沿鐵芯閉合，同時(shí)與原、副繞組相鏈，是 變壓器能量變換和傳遞的主要因素，稱為主磁通或互感磁通；另一部分磁通必$ 主要是通過非磁性介

14、質(zhì)（空氣或油），它僅與原繞組全部相鏈（只與原繞組部分 匝數(shù)相鏈的露刺痛已折算為全部原繞組相鏈而數(shù)值減少的等效磁通），故稱它為 原繞組的漏磁通。根據(jù)電磁感應(yīng)定律，當(dāng)磁通族和九隨時(shí)間變化時(shí)，分別在它們所交鏈的繞 組內(nèi)感應(yīng)電動(dòng)勢(shì):dt 她 dt(2.1)“一吟式中與、立是主磁通在原、副繞組所感應(yīng)的電動(dòng)勢(shì)瞬時(shí)值；氣是原繞組漏磁通 在原邊感應(yīng)的電動(dòng)勢(shì)瞬時(shí)值。冬=佻k=曳I產(chǎn)L所以，與 % ,設(shè)變壓器的變比為q，則0產(chǎn)破2 J k。所以利用變壓器可以在傳輸電能的同時(shí)改變其電壓和電流。2.2 三相變壓器的等效電路及聯(lián)結(jié)組現(xiàn)在電力系統(tǒng)都采用三相制，所以實(shí)際上使用得最廣泛的是三相變壓。從 運(yùn)行原理來看，三相變壓

15、器在對(duì)稱負(fù)載下運(yùn)行時(shí)，各相的電壓、電流大小相等, 相位彼此互差120 ,故可任取一相分析，即三相問題可簡化為單相問題。根據(jù)變壓器原、副繞組電動(dòng)勢(shì)的相位關(guān)系，把變壓器繞組的連接分成各種 不同組號(hào)稱為繞組的連接組。在不同的連接組下，三相變壓器的等效電路略有 不同?，F(xiàn)以Ydll連接組為例，做三相等效電路等效電路圖如圖15所示。Yd11連接組三相等效電路-9 -在三相變壓器中，用大寫字母A、B、C表示高壓繞組的手段，用X、Y、Z 表示高壓繞組的末端；低壓繞組首、末端則應(yīng)用對(duì)應(yīng)的小寫字母”、b、c和X、 y、z表示。星形連接的中點(diǎn)用字母0表示。不論原繞組或副繞組，我國主要采用星形和三角形兩種連接方式。為

16、了形象地表示原、副邊電動(dòng)勢(shì)相位地關(guān)系，采用所謂的時(shí)鐘表示法：即 把高壓繞組的電動(dòng)勢(shì)向量作為時(shí)鐘的長針并指向12,低壓繞組的電動(dòng)勢(shì)相量作 為時(shí)鐘的短針，其所指數(shù)字作為單相變壓器連接組的組好。在我國生產(chǎn)的變壓器中，以Ydll、YndlK YnyO. YyO （n表示中性點(diǎn)接地） 四種連接組為主。第三章 變壓器仿真的方法從20世紀(jì)60年代開始，人們就花費(fèi)大量的精力去解決變壓器的計(jì)算機(jī)模 型問題。由于變壓器的非線性特性，這被證明是困難的課題。不像線性系統(tǒng)一 樣，沒有一般的解決方案可以解決非線性方程。即便是數(shù)字式的解決方案，也 只能很困難的解決某一類的非線性方程，在穩(wěn)定的狀態(tài)下，存在好的變壓器模 型。然

17、而，在瞬變的狀態(tài)下，還沒有完全令人滿意的變壓器模型。變壓器的性能主要取決于其鐵心的磁化特性，即鐵心的磁滯回環(huán)，因此對(duì) 鐵心磁滯回環(huán)的擬合是最基礎(chǔ)、最重要的工作。在變壓器特性的數(shù)值仿真計(jì)算 中，對(duì)磁滯回環(huán)的擬合提出了以下幾個(gè)要求：具有較高的精度；在大范圍 內(nèi)不分段，具有光滑性，否則會(huì)引起變壓器特性仿真計(jì)算過程的不穩(wěn)定；具 有稠密性，因?yàn)橥ㄟ^實(shí)驗(yàn)只能得到有限條磁化曲線，而變壓器仿真中需要知道B 一H平面中的任意一條曲線。由于鐵心的飽和特性、磁滯現(xiàn)象等非線性因素的影 響，很難用數(shù)學(xué)模型精確地描述鐵心的動(dòng)態(tài)磁化過程。變壓器通過鐵心磁場作用建立一次側(cè)和二次側(cè)的電磁聯(lián)系。因此變壓器暫 態(tài)建模的關(guān)鍵是對(duì)鐵心

18、動(dòng)態(tài)磁化過程的數(shù)學(xué)描述。根據(jù)對(duì)磁化特性曲線描述的-io -不同，現(xiàn)有研究用的變壓器模型大致有下列4種：（1）模型A基于基本勵(lì)磁曲線的靜態(tài)模型；（2）模型B基于暫態(tài)勵(lì)磁特性曲線的動(dòng)態(tài)模型；（3）模型C一一基于暫態(tài)勵(lì)磁特性曲線的非線性時(shí)域等效電路模型;（4）模型D一基于ANN的變斜率BP算法創(chuàng)建的模型。3.1基于基本勵(lì)磁曲線的靜態(tài)模型基于基本勵(lì)磁曲線的變壓器模型只考慮飽和引起的非線性，即采用如圖3 一1所示的基本磁化曲線作為變壓器暫態(tài)工作特性曲線進(jìn)行二次側(cè)電流的計(jì)算。 等效電路圖如圖3-2所示。/Wb圖3-1基本勵(lì)磁曲線3-2靜態(tài)模型由磁通守恒和KCL定律可以得到以下基本方程組:'乂a-，

19、。）=心2-力 _ _ -力(3.1)式中4為一次側(cè)電流；'。為勵(lì)磁電流；2為二次側(cè)電流；。為主磁通；M、 生為一、二次側(cè)匝數(shù)；&、為二次側(cè)負(fù)載。由方程組（31）中的第一和第三個(gè)方程得到-nJ 0,將“ 嘰 出代入方程組第二方程，整理可得:NRI + NL 幺-N1R,i0 /1 N/ 史 + N】Ldt 7 dtu力。(3.2)-43 -.=竺_ £ N、dB因° = 8$和° M,故有山。/ 代入式（3.2）可得:f J > Y(3.3)*=(N悶+N吟- N網(wǎng)°)/(N；s清：十 %也)用四階龍格一庫塔法或隱式梯形公式就可以求

20、解一階常微分方程式（3-2） 或式（3-3）,從而建立了變壓器仿真數(shù)學(xué)模型。3. 2基于暫態(tài)磁化特性曲線的動(dòng)態(tài)模型這類變壓器模型建立在對(duì)動(dòng)態(tài)磁化特性曲線的數(shù)學(xué)描述之上。暫態(tài)磁化特性曲線° =的描述，最常用的是采用極限回環(huán)壓縮法。即假定鐵心磁化曲線的主磁滯回環(huán)和次磁滯回環(huán)具有相似性，由主磁滯回環(huán)壓縮 生成次磁滯回環(huán)。例如用反正切函數(shù)擬合主磁滯回環(huán)，其表達(dá)式為：（/0） = a arctan h（i0 ± C） + /7z0八人式中a、和°為常數(shù)。在上升軌跡修可和下降軌跡倍可的轉(zhuǎn)折點(diǎn)將主磁滯回環(huán)K =1""。| 一（編"一 £）

21、/|"arctan/?（F" ±C）±g 按壓縮系數(shù)L2_2向直線g3三巧（32壓縮生成次級(jí)回環(huán)的下降支或上升支。圖33所示為動(dòng)態(tài)磁化特性曲線，其中工&）,力&）為極限磁滯回環(huán),（。"二0"）為轉(zhuǎn)折點(diǎn)（假設(shè),從上升變成下降），則力&）為經(jīng)過該點(diǎn)的次級(jí)回 環(huán)下降支。在力”。）形成的回環(huán)內(nèi)的人&）部分為經(jīng)過該轉(zhuǎn)折點(diǎn)的暫態(tài) 磁化軌跡。q/wb 0.61-0.6 - I &/A-3-2-101233-3局部磁滯回環(huán)軌跡3. 3非線性時(shí)域等效電路模型該模型用幾個(gè)電路元件分別模擬造成變壓器非線性的因素。因?yàn)?/p>

22、引起變壓器非線性的主要因素有飽和、渦流和磁滯，所以用三個(gè)電路 元件模擬這些因素，并將各元件流過的電流線性疊加，得到勵(lì)磁電流。其表達(dá) 式為：10 =，,+.（35）式中im為磁化電流；i力為磁滯電流；為渦流電流。因?yàn)樘蕹似渌绊懸蛩囟鴨为?dú)進(jìn)行考慮，故"可以用無磁滯曲線（基本磁 化曲線）來表示，這是一個(gè)僅僅與磁鏈有關(guān)的表達(dá)式。其表達(dá)式可以表示為= 呻（36）磁滯是由交變電流產(chǎn)生，其大小和電壓以及頻率有關(guān)。但實(shí)驗(yàn)表明，在 50Hz到400Hz內(nèi)，磁滯隨頻率的變化而改變得很小，故頻率的影響一般用一個(gè)常數(shù)表示。磁滯電流部分的表達(dá)式為:(3.7)j 、尸-i 絲dt其中夕為斯坦梅茨(Stein

23、metz)系數(shù)，由鐵磁材料的特性決定。除設(shè)定為 在50Hz下的一個(gè)常數(shù)。渦流電流和磁通、磁通變化率以及頻率有關(guān)。但是在電流頻率不超過400Hz 的情況下，渦流電流不會(huì)因頻率改變而顯著變化。因此可以不考慮頻率變化對(duì) 渦流的影響，表達(dá)式如下：(3.8)綜合以上各式可得到考慮了飽和、磁滯和渦流影響的變壓器勵(lì)磁電流暫 態(tài)數(shù)學(xué)模型，其表達(dá)式為：i。)+(勺嚕廠+二)務(wù)令 心二G =與(華產(chǎn)2 +廠則有iojJvdt+R;%(3.9)其等效電路如圖34所示。3-4非線性時(shí)域等效電路模型第四章三相變壓器的仿真電力系統(tǒng)中的變壓器通常是三相的，而三相變壓器的磁路結(jié)構(gòu)型式、繞組 接線方式（Y結(jié)、D結(jié)）、中性點(diǎn)接地

24、與否等多種因素對(duì)勵(lì)磁涌流的大小和波形有 著較大影響，故本文僅對(duì)電力系統(tǒng)中最常見的Ydll、Yndll> YnyO、YyO （n表 示中性點(diǎn)接地）接線的三相三柱心式變壓器進(jìn)行仿真研究。為簡化分析，在研 究變壓器空載合閘哲態(tài)過程時(shí)忽略鐵心的損耗，認(rèn)為勵(lì)磁支路為純電感支路。4. 1三相變壓器仿真的數(shù)學(xué)模型首先對(duì)各種不同連接組情況下，根據(jù)電路原理的基礎(chǔ)知識(shí)，建立三相 變壓器的數(shù)學(xué)模型。4.1.1三相變壓器Yd11連接組模式圖11為Ydll接線的變壓器的三相接線圖和單相等效電路。圖4-1 Yd11接線得變壓器空載合閘時(shí)三項(xiàng)接線圖和單相等效電路當(dāng)Y側(cè)空載合閘后其暫態(tài)方程如下：di與=a+機(jī)+w）9+

25、w n+，atdif =（G + 八吊 + W + 乙沖 + （4. Datdi/=a+）（ +（4+qat式中，5為Y側(cè)中性點(diǎn)電壓，其它符號(hào)見圖31?？紤]到一次為Y接線，二次為D接線，所以有：。+ h + 4 = 0（4. 2）(4. 3)式（4.3）,化簡、+%+金=3（d* + 缶，0）而Ua+Ub+Uc=o,將式（4.1）三式相加并計(jì)及式（4.2）、 得：liN = 一出）寸 + rDiD）（4.4）又由單相等效電路可知： 汕；（4.5）4 = + L式（4.5）三式相加得到：認(rèn)（4.6）d(b N；SdB, dinw dim, e ="=_2 = Kdt l/m dtdtK

26、imj(j = ahc) (4.7)式中in"電流i*的導(dǎo)數(shù)。將式(4.4). (4.5)、 組漏抗近似相等(n=r»(4.6)、(4.7)代入方程式(4.1),計(jì)及一、二次繞LfLd)，經(jīng)化簡得：3ua3%3u34+2,rsrs3M+2mamb34+3( + 24-A-L,3，+3K/,+24L,LL,3，+3K + 2zJ若忽略系統(tǒng)阻抗，即r.=0, L=0, U=0,則上式可化簡為:O' 0(4.8)式中：K, = (j = a,b,c)動(dòng)態(tài)感應(yīng)系數(shù)4，4。電源內(nèi)部等值正序電感與零序電感S , la > lh > lc變壓器鐵心截面積與各相磁路長度

27、L，Na, Nb，N,次繞組漏抗和各相匝數(shù)5 3;電流乙，十，一的導(dǎo)致4.1.2三相變壓器Ynd11連接組模式Y(jié)ndll接線的三相變壓器Yn側(cè)空載合閘時(shí)，其暫態(tài)方程為:%=a+也+w+4冷+(Ao - 4) atdjuh=a+r Vb+區(qū)+4)+(Ao L) at%=a+ri)< +(4+4)牛+(4o-4) atdiD dt 嘎 dt ，魚+ eh (4. 9)clt考慮到一次為Yn接線，二次為D接線，所以:(4. 10)繪= 3（L）霽 +Nd）(4. 11)又4+%+/=。，則式（4.9）三式相加得:。=Wo+4）牛+a+町+七窯+江 同樣將式（4. 5）三式相加得：1。=；（/+

28、強(qiáng)+乙）+ 7。(4. 13)(4.12)將式（4.5）、（4.7）、（4.13）代入方程式（4.9）并聯(lián)立式（4.12）,若不 計(jì)系統(tǒng)阻抗且認(rèn)為變壓器一、二次繞組漏抗相等，則可得其空載合閘狀態(tài)方程:%ubUc0100rimambme0。+勺064.1.3三相變壓器YnyO連接組模式bnc(4.14)YnyO接線的三相變壓器Yn側(cè)空載合閘時(shí)，其暫態(tài)方程與Yndll接線一樣, 如式（4.9）所示。因?yàn)橐淮螢閅n接線，二次為y接線，所以"+ib+ic=3i°；Jd= 0因而由單相等效電路可得:(4.15)me(4.16)乙+i力+乙=3，。(4.17).18)同樣，根據(jù)類似的推

29、導(dǎo)過程并計(jì)及前述各假設(shè)，可得該接線三相變壓器空 載合閘狀態(tài)方程：L + K,004.1.4三相變壓器YyO連接組模式Y(jié)yO接線的三相變壓器空載合閘時(shí)，其暫態(tài)方程與Yndll接線一樣，同樣如 式(4.1)所示?？紤]到一次為Y接線，二次為y接線，貝II：(4.19)從而可得:me又u.+Ub+%=0,將式(4.1)三式相加并計(jì)及式(4.19),化簡得:(4. 20)3匹=_(，+與+，)同樣，根據(jù)類似的推導(dǎo)過程，可得YyO接線得三相變壓器空載合閘狀態(tài)方 程為：me%+ 2(Ka-K.Kb+ 2(*nia34+2("tnc(4. 21)至此，式(4.8)、(4. 14). (4.18)、(

30、4. 21)和各相動(dòng)態(tài)磁化曲線鳥二"4,)及H, 二(Im/ N/),(7=4,c)構(gòu)成了 Ydll、Yndll> Yny0> YyO 接線三相變壓 一次側(cè)空載合閘的基本方程。4. 2電源電壓的描述根據(jù)前述假設(shè)，電源電壓u (相電壓)可用式(4. 22)描述。ua = (Um /1.73)*sin(69t+ a)< 叫=(U3L73)sin(Gt+o + 120 )(4. 22)uc = (Uni /1.73)«sin(d>t+ a -120 )式中，Um為電源線電壓峰值，取1.1倍額定電壓。a為A相空載合閘初相 角。在用Matlab仿真得過程中，a

31、的設(shè)定并非是一個(gè)可以輸入的變量，如果需 要改變初相角，可在程序內(nèi)部直接改變相電壓u。4. 3鐵心動(dòng)態(tài)磁化過程簡述根據(jù)試驗(yàn)所得到變壓相鐵心磁化曲線數(shù)據(jù)分段擬和其極限磁滯回環(huán)是我們 的基本原理。由試驗(yàn)所得到的數(shù)據(jù)可以幫助我們界定程序中一些參數(shù)，而如何 選擇界定函數(shù)將很大程度上影響試驗(yàn)仿真得結(jié)果。在這次的試驗(yàn)計(jì)劃中，我們 將選擇兩種方式(即選擇不同的函數(shù)逼近)進(jìn)行仿真，然后分別討論兩種方案 的優(yōu)劣，得出最佳的方案。第一種是比較簡單的模式，基本上不考慮曲線進(jìn)入飽和區(qū)的情況(盡管飽 和區(qū)是不可回避的問題，但這樣做亦不失其合理性，這一點(diǎn)將在后面被討論 到。)，采用兩條修正的反正切函數(shù)做為極限磁滯回環(huán)。然后

32、，對(duì)于主區(qū)間內(nèi)的 動(dòng)態(tài)磁滯回環(huán)，根據(jù)不同的轉(zhuǎn)折點(diǎn)和運(yùn)行趨勢(shì)對(duì)極限磁滯回環(huán)向極限磁滯回環(huán) 擬合。極限磁滯回環(huán)的數(shù)學(xué)描述由于和第二種情況相近，只是將第二種方式的 飽和區(qū)考慮在外，所以具體方法將不再贅述，可以參考4. 1.3.1部分。對(duì)于暫 態(tài)局部磁滯回環(huán)的描述，具體方法可以參考4. L 3. 2部分。對(duì)于剩磁的處理的 處理，具體方法可以參考4.1.3. 3。第二種是比較復(fù)雜得模式，需要在第一種的情況下考慮飽和的問題。這種 方法不但描述了鐵心的飽和特性，而且能夠反映鐵心的磁滯特性?；驹硎? 首先，格局試驗(yàn)所得的變壓器鐵心磁化曲線數(shù)據(jù)分段擬合其極限磁滯回環(huán)：(1)對(duì)于未飽和時(shí)主區(qū)間內(nèi)的兩條極限磁滯

33、回環(huán)，采用修正的反正切函數(shù)加以擬合;(2)對(duì)于飽和后主區(qū)間外的磁化曲線，認(rèn)為其已進(jìn)入線性可逆區(qū)(直線段), 采用兩條平行的直線段加以描述。然后，對(duì)于主區(qū)間內(nèi)的動(dòng)態(tài)磁滯回環(huán)，我們 根據(jù)其不同的轉(zhuǎn)折點(diǎn)和運(yùn)行趨勢(shì)對(duì)極限磁滯回環(huán)向飽和后的兩條平行直線進(jìn)行 壓縮，就可得變壓器鐵心實(shí)際運(yùn)行的動(dòng)態(tài)磁化軌跡。下面，就以第二種方法為例，詳細(xì)的解釋一下極限磁滯回環(huán)的描述、暫態(tài) 局部磁滯回環(huán)的描述及剩磁的處理等等問題。4. 3.1極限磁滯回環(huán)的數(shù)學(xué)描述A.主區(qū)間內(nèi)【一及 HJ極限磁滯回環(huán)可用下式所示的修正反正切函數(shù)表示。 a*p = tan_1 /3(H-C) + /H(4. 23)a/3 = tan-1 隊(duì) H

34、+ C) + /H(4. 24)式中，參數(shù)a、4、/、C可根據(jù)實(shí)測(cè)磁滯回環(huán)數(shù)據(jù)由非線性的曲線擬合 程序求得。因此：dB 1 p=5 + y >dH a 1 + 4("一0丁(4. 25)dBp ,(舊+夕(“ + C)于+/(4. 26)B.飽和后的磁化曲線(|H>Hz)擬合為兩條平行的直線段。 當(dāng)H>Hz時(shí)，b(h)= b：+號(hào)工)(4. 27)當(dāng)*旦時(shí)，B _ BHH工(4. 28)因而,dB B -及dH H、一H.$z(4. 29)4. 3. 2暫態(tài)局部磁滯回環(huán)的描述由于鐵心材料電磁性能的復(fù)雜性，對(duì)動(dòng)態(tài)局部磁滯回環(huán)的精確仿真是比較 困難的。但因極限磁滯回環(huán)已

35、描述了磁滯的基本輪廓，故根據(jù)不同轉(zhuǎn)折點(diǎn)對(duì)其 進(jìn)行壓縮就可近似模擬動(dòng)態(tài)磁化過程中的某一段上升軌跡和下降軌跡。它分兩 種情況模擬。dB/dH<0,運(yùn)行點(diǎn)下降軌跡由于極限磁滯回環(huán)左側(cè)描述了減磁過程，將極限磁滯回環(huán)左側(cè)回線在縱軸 方向按比例地朝直線3 =。壓縮，可得一簇下降曲線。對(duì)通過某一轉(zhuǎn)折點(diǎn)（H（0>,B<0>）的運(yùn)行點(diǎn)下降軌跡可由左極限磁滯回環(huán)按壓 縮系數(shù)KX向直線8 =?！?J）/。壓縮而得（如圖42所示）。此處KX=（B（°） -b：）/（b-b2）則通過該點(diǎn)的下降軌跡為：。那-(加。，-為4B = -.2_ x tan-,/7(H+C) + - +-I /

36、H- (4. 30)小小加(。)+ 0 +曰L2屋 2；dBTh因此：(4. 31)廣夕”a tan-1(H(0) + C)+ y 1 + /?(/+C) a圖42局部磁滯回環(huán)的模擬dB/dH>0,運(yùn)行點(diǎn)上升軌跡同樣，將極限磁滯回環(huán)右側(cè)向直線8 = (yH + J)/a方向壓縮，可得通過轉(zhuǎn) 折點(diǎn)(H(o), B)的運(yùn)行點(diǎn)上升軌跡遜-(“。)+勺"B = r=px tan_1 p(H -C) + yH + (4. 32)a tan-'/7(H<(,)-C)-y *-22)i+夕(c)f+?dB 。那-(加°)+勺 ±xdH atai,("

37、;-Cj f第五章變壓器MATLAB仿真研究5.1 仿真長線路末端電壓升高5.1.1 仿真模型如圖：5.1.2 仿真參數(shù)介紹及波形模型窗口參數(shù)如下圖:&P'AmeteryJ«uh/C00938 (Actn<) .X1Y8SZtUg Q3SolerPat a Z*port/xpoxt DiicnortxcfHacdviKt lApltMnta.Keffr«rK>n< Sunilttioti Tm«， Codt C*n«raticfi ICCenerat&oct S3CSt»rt tiM： 0.0Stop

38、tiM： 1.0Solwr opiicnsTyp*:、“4八，p*j S«lwrsIix rt«p fixeMtolUlativ* t«l«rtr»c«:IIlin rttp »z«:wtoAbsolut# tol«rarkc«:aut。Initial «t«p ,皿：MoSMp* pr”rv"g：PnibUall.$olv*r r=6 Mth»d: RobuitThree-Phase Source 參數(shù)如下圖：。由twRvfTiiP“，"&qu

39、ot;， LoU F1o5fe-tor-pk>xe r» elC“o (V)s SMeSFhue mtlt of Zimm A (tetxret) iI 0Ffi«W7 51)3 MInt «m>l ewxct ic<t: T1.Spec仃 ii9edance ufin( rhoH*circuit levelSeuee« r”ftw (%u)， 0.8929Swirce Lntctence QD:Bar，soly, (Vr*» pK-ph.)z 500e3|l- M 匚fiiwL匚 Jtto 1 3Multimeter 的參數(shù)

40、:其中，Us_phl_gnd代表Scope中的實(shí)線，Ur_phl_gnd代表虛線。Powergui 的參數(shù)：將 Simulation type 選為 Continuous,將 Lond flow frequency 改為 50Hz 即可。Scope的波形如下：（長度為300km）將 Distributed Parameters Line 參數(shù)中的 Line Length 改為 500km，則 Scope的波形為:改為1000km,波形為:可見，分布參數(shù)導(dǎo)線長度越長，其末端電壓Us升高越明顯。5. 2仿真三相變壓器T2的勵(lì)磁涌流仿真模型如圖:仿真參數(shù)介紹及波形模型窗口參數(shù)不變。由勵(lì)磁涌流的特性可

41、知：當(dāng)變壓器在電壓過零點(diǎn)合閘時(shí)，產(chǎn) 生最大的勵(lì)磁電流；當(dāng)變壓器在電壓最大值時(shí)合閘，不會(huì)產(chǎn)生勵(lì)磁電流。因此先 仿真三相變壓器的電壓。三相電源參數(shù)、Powergui參數(shù)不變。分布參數(shù)導(dǎo)線長度設(shè)為300km。QF3參數(shù)如下圖：P>xa»etersiMti&l ms of breik4ri |el”ed=)Sntchu< of phue A Snt chine ©f Phud B Svit chine ©f phise CTr ant iti<on(力Id愉Externil control af rv&tehin< UmsSreak

42、erx refinance Kenottfu) 0. 04BSnubbers rexkrt ancc Xp (Ghw*)I4Srm雙打 CQuie 8 (F"M) infMore- 般 fapodl :三相變壓器T2的參數(shù)如下圖:萬用表選擇測(cè)量的量由上到下依次為Uag_w2: T2、Ubg_w2： T2、Ucg_w2: T2。Demux參數(shù)的輸出量設(shè)為3 o則示波器的波形為：由上圖可得：Uag_w2: T2、Ubg_w2： T2和Ucg_w2: T2過零點(diǎn)時(shí)間可分別為8 02、0. 0267、0. 0234。為峰值的時(shí)間可分別為 0.025、0. 0317、0. 0384。改變QF3

43、的參數(shù)，如下圖：上圖表示在0.02秒時(shí)斷路器閉合。萬用表選擇測(cè)量的量由上到下依次為lexjA:T2、Iexc_B: T2、Iexc_C: T2o其他模塊參數(shù)不變。則示波器的波形為：將QF3的Transition times分別改為：0. 0267,則示波器波形為:改為0.0234,則波形如下圖:從上述波形可以看出，變壓器在某一相的電壓過零點(diǎn)合閘時(shí)，此相產(chǎn)生最大的 勵(lì)磁電流約為1200A,且經(jīng)過0. 2s左右衰減至穩(wěn)態(tài)運(yùn)行時(shí)的勵(lì)磁電流，峰值約 為 20Ao 將 QF3 的 Transition times 分別改為:為 025、0.0317、0.0384, 則示波器的波形分別為：過渡時(shí)間為0.

44、025s過渡時(shí)間為0.0317sia»40Te ©fHet 0過渡時(shí)間為0. 0384s從上述波形可以看出，當(dāng)變壓器在某一相電壓峰值時(shí)合閘，此相不會(huì)產(chǎn)生勵(lì)磁涌流, 而其他兩相則一定會(huì)產(chǎn)生勵(lì)磁涌流。5. 3三相變壓器仿真模型仿真三相變壓器外部故障仿真模型如圖:仿真三相變壓器T3的勵(lì)磁涌流仿真模型如圖:物4切P麗 £scope.仿真三相變壓器T3的內(nèi)部故障仿真T3相間短路（AB相）的模型如圖:3cTrams>2m m ，。 。 ，&e,一；加仿真T3匝間短路的模型如圖:Te七WE2 r2C-nnasa %？25.4變壓器仿真波形分析5.4.1 對(duì)勵(lì)磁涌流

45、進(jìn)行FFT分析選擇變壓器T2的勵(lì)磁涌流波形圖進(jìn)行FFT分析，如下圖:圖5-1分析了 input 1的lexjA: T2的波形。從0. 02s開始分析，分析 兩個(gè)周波。以柱形圖顯示，橫軸坐標(biāo)為頻率，最大值為300Hzo.<<r in圖5-2分析了 input 2： Iexc_B: T2的波形。從0.02s開始分析，分析兩個(gè)周波。Signal to arulyze. Dspitr leiected wgm £ DecMry fFT wnSowTmrw(OFFT mlyrtiFundifwrtal (50Hz) = 6109. THD= 188 13%Marmoftic o&l

46、t;(kr80W4020 (-T£PEg 一。«.6 yAvrtatte tigruls-aSccotOeiairut：O9U83$jgnal nutfvMr T.FFT vmMcmStan tm(» 0 s3，WRMti freou«y (to)FFT i*ettmg$6«pey UyW 二 e«r (relKh< a*BoM山.|tFrtoutr)Gy am KBnnxic order.Wax Fr«w«ncy (Hz) w .9? 4 L 010M以柱形圖顯示，橫軸坐標(biāo)為頻率，最大值為300Hz&

47、;e Edit %«b Insert look Qeiktop Window fcfeip'd j k » 4 0 口目口圖5 3分析了 input 3： Iexc_C: T2的波形。從0.02s開始分析，分析兩個(gè) 周波。以柱形圖顯示，橫軸坐標(biāo)為諧波次數(shù)，最大值為6。由此可得：勵(lì)磁涌流含有大量高次諧波分量，且以二次諧波為主。5.4.2對(duì)外部故進(jìn)行FFT任選一個(gè)短路電流（圖的Iaw3: T2）進(jìn)行FFT分析，如下圖:far Lda yi<wloch QrUdopBeipsignXsSWfCtMr©aFFT wiodow，«4f “eyeM 2

48、rendtrnrui WQtftcy MU Ud， R Q SMito Mwtyx. 3 Dtty UMM 41Oae*y FFT wM&wSlKty 沖M »eydM FFT wvito* (m rM) 2分3”, Srwr111M crMCAc toFFT MttngsU«x.r«MtKv3)-45 -圖54分析了 Iag_w3: T2的波形。從0.2s開始分析，分析兩個(gè)周波。以 柱形圖顯示，橫軸坐標(biāo)為頻率，最大值為150Hz。由圖可知：外部故障時(shí)的不正常運(yùn)行電流是標(biāo)準(zhǔn)的正弦波，不含有二次諧波。5.4. 3對(duì)內(nèi)部出亍FFT任一個(gè)短路電流（圖的Iw2: Saturable Transformer）進(jìn)行FFT分析，如下圖:b口目 AyatoUe W ScnvdB/eSnilto analyzeFFTvnnddw5untrai(