基于MATLAB的電壓閃爍下電弧爐特性仿真研究

2012-01-15 06:03:08姚俊濤史俊霞

電子設計工程 2012年16期

關鍵詞：模型

姚俊濤，史俊霞，朱婷

（武漢大學電氣工程學院，湖北武漢 430072）

電弧爐是一種固有的非線性時變負載，它能引起很多的電能質量問題，比如：諧波和電壓閃爍等。電弧爐的操作會產生奇偶諧波電流，這些諧波電流在電網中流通時會產生諧波電壓，諧波電壓反過來又會影響其他的用戶。

閃爍是人眼在物體光照強度改變時的一種經驗感覺。光照中的最大敏感性改變的頻率范圍是5～15 Hz。由于電弧爐是一種大閃爍源，導致在相關電網中的電壓波動。因此，電弧爐的建模吸引了很多的電力系統工程師，以此來解決電弧爐帶來的相關問題。

弧的建模中最重要的問題就是弧的模擬，現在有很多模型來描述電弧[1]，幾乎都使用了平衡穩態方程，基于微分方程的時域模型也被許多人應用。其他的模型比如頻域響應、V－I特性被引進用來分析EAF特性[2]。通過比較時域和頻域模型可以發現，時域模型對EAF的研究更有幫助。以上的研究中，有許多的局限存在，比如微分方程的初始條件，三相電流的平衡位置以及用復雜的數學方程來模擬EAF等。

文中對在時域用MATLAB對EAF模型進行模擬。本文所用模型能夠體現出EAF的主要特性，并且在模擬過程中采用一種更好的近似方式而不需要去設置EAF的初始條件。同時本文所用模型能夠很好的描述EAF的不同操作條件以及它對相關電網的影響。

1 電網模型

圖1描述了一個簡單的單相電網，用它來為EAF供電。

圖1 為EAF供電的電網模型Fig.1 Electric network supplying an EAF

在圖1中系統阻抗用Zs表示，PCC總線代表通用的耦合點，AF總線表示變壓器的低電壓邊，它的阻抗是Zt，系統的整體參數如表1所示。

2 作為負載時的電弧爐建模

在文章的這部分，將電弧的電壓和電流估算來進行電弧爐建模，主要是基于參考文獻[1－3]所說的真實VIC劃分和線性分段。下面主要討論本文所用的基于VIC的EAF方程。電弧VIC模型，該模型綜合雙曲線和指數模型，表述如下：■

表1 電力系統參數及電弧爐參數Tab.1 Electric system and EAF model parameters

其中：V和i分別表示每一步的弧電壓和弧電流，Vat是電壓門限值，C與D結合是通過功率和電流來定義弧長，從而決定門限電壓，I0是電流常數，用來抵消異常電流值，是符號函數。

文中所用的VIC模型能夠準確描述時域電弧爐特性，并且可以解釋電弧爐的各種工作條件[4－5]。

3 仿真電路

如圖1所示的電力系統中的電弧爐特性通過SIMULINK/MATLAB建立仿真電路，如圖2所示。

由于電弧爐是一種非線性時變的電壓控制源，文中用function/MATLAB程序來模擬它。動力學電弧爐模型需要對弧的影響進行實時分析，通過不同的弧電導率來模擬動力學弧特性。

圖2 MATLAB仿真圖Fig.2 Complete Simulink/MATLAB file of simulation

一般而言，所有的變化都是隨機的自然的，但是為了比較，本文主要研究了兩種變化：正弦變化和隨機變化。為了研究電壓閃爍對電弧爐的影響，電壓門限值Vat也是正弦和隨機變化的，Vat被進行如下調制[7]。

1）正弦變化

其中：m是調制指數，ωf是閃爍頻率。

2）隨機變化

Vat被加載一個均值為0的隨機信號，頻率范圍為：4～14 Hz。因此有：

其中：N（t）是一個限帶的白噪聲，并且均值為0，方差為1。在Simulink/MATLAB中這兩種調制的仿真圖如圖3所示。

圖3 隨機閃爍和正弦閃爍仿真Fig.3 Simulink/MATLAB simulation of random flicker and sinusoidal flicker

4 仿真結果及分析

文中為了更好地對比穩態時和電壓閃爍發生時電弧爐的主要特性的變換，設置了以下幾種參數的比較：1）弧電流和弧電壓（Arc Current and Arc Voltage）；2）伏安特性（VIC）；3）弧電導率變換（Arc Conductance）；4）能量功耗（P 為有效功率；Q為無功功率）。

穩態弧電流和弧電壓及伏安特性仿真效果如圖4所示；穩態弧電導率及能量功耗變化仿真效果如圖5所示；正弦閃爍下的弧電流和弧電壓及能量功耗仿真效果如圖6所示；隨機閃爍下的弧電流和弧電壓及能量功耗仿真效果如圖7所示；隨機閃爍下的伏安特性及弧電導率仿真效果如圖8所示。

1）穩態特性仿真

正如文中所描述的一樣，電壓閃爍對電弧爐系統的影響可以用電壓隨著時間的變換來仿真研究。通過仿真分析可以看見，加入電壓閃爍后，弧電壓和弧電流的變換周期幾乎和穩態時是一致的，但是對于幅度則是隨著閃爍的強弱而不斷變化。伏安特性也比穩態時更加復雜凌亂多變。對于能量功耗而言，穩態時的有效功率明顯由于無功功率，加入電壓閃爍后，無功功率會出現瞬時的增強，有效功率則會出現瞬時的降低，致使系統幾乎無功率輸出（死機），同時有效功率和無功功率的變換非常劇烈，導致系統的輸出非常不穩定，但是隨著時間的推移，系統將慢慢趨于穩定，在PCC處出現如穩態時的可靠性。加入電壓閃爍后，弧電導率變得波動異常的大。

圖4 穩態弧電流和弧電壓及伏安特性仿真Fig.4 Simulation of Arc Voltage and Arc Current and VIC under steady state

圖5 穩態弧電導率及能量功耗變化仿真Fig.5 Simulation of Arc Conductance and power consumption variation under steady state

2）正弦閃爍仿真

圖6 正弦閃爍下的弧電流和弧電壓及能量功耗仿真Fig.6 simulation of Arc Voltage and Current and power consumption variation under sinusoidal flicker

3）隨機閃爍仿真

圖7 隨機閃爍下的弧電流和弧電壓及能量功耗仿真Fig.7 simulation of Arc Voltage and Current and power consumption variation under random flicker

圖8 隨機閃爍下的伏安特性及弧電導率仿真Fig.8 Simulation of VIC and Arc Conductance variation under random flicker

5 結論

通過在電壓閃爍下對電弧爐特性進行仿真分析，可以了解其特性隨電壓閃爍的變換規律，這對于電弧爐的設計以及對電網中的電壓閃爍的防御具有很好的指導作用。對于電壓閃爍引起的電網系統的各種特性參數突變，可以通過增強靜態參數補償，添加自適應的有源濾波器等來改善在PCC處的電能質量。

[1]Zheng T，Makram E B，Girgis A A.Effect of different arc furnace models on voltage distortion[C]//IEEE Transactions，InternationalConferenceon Harmonicsand Quality of Power，1998（2）:1079-1085.

[2]Golkar M A，Meschi S.MATLAB modeling of arc furnace for flicker study[C]//IEEE Conference on Industrial Technology（ICIT），2008（5）:1-6.

[3]Hooshmand R，Banejad M，Esfahani M T.A new time domain model for electric arc furnace[J].Journal of Electrical Engineering，2008，59（4）:195-202.

[4]祁碧茹，肖湘寧.用于電壓波動研究的電弧爐的模型和仿真[J].電工技術學報，2000，15（3）：31-35.QI Bi-ru，XIAO Xiang-ning.Modeling and simulation of an arc furnace for voltage fluctuation investigation[J].Transactions of China Electrotechnical Society，2000，15（3）:31-35.

[5]王育飛，潘艷霞，姜建國.基于MATLAB的交流電弧爐隨機模型與仿真[J].高電壓技術，2008，34（5）： 973-977.WANG Yu-fei，PAN Yan-xia，JIANG Jian-guo.Stochastic model of AC electric arc furnace based on MATLAB[J].High Voltage Engineering，2008，34（5）:973-977.

[6]劉小河，楊秀媛.電弧爐電氣系統諧波分析的頻域方法研究[J].中國電機工程學報，2006，26（2）：30-35.LIU Xiao-he，YANG Xiu-yuan.Research on harmonic analysis of arc furnace electric system based on frequency domain method[J].Proceedings of the Chinese Society for Electrical Engineering，2006，26（2）:30-35.