徐云聰 張 逸 林才華 劉必杰

基于電氣特性-物理參數(shù)耦合的交流電弧爐模型

徐云聰1張 逸1林才華2劉必杰3

（1. 福州大學(xué)電氣工程與自動化學(xué)院 福州 350108 2. 國網(wǎng)福建省電力有限公司電力科學(xué)研究院 福州 350007 3. 國網(wǎng)福建省電力有限公司寧德供電公司 寧德 352100）

現(xiàn)有電弧爐模型主要關(guān)注電弧內(nèi)部與外部特性的聯(lián)系，割裂了電弧內(nèi)部微觀下電氣與物理的耦合影響，難以體現(xiàn)電弧的電熱物理本質(zhì)，導(dǎo)致仿真精度有所欠缺。針對上述問題，該文提出一種改進的新型交流電弧爐模型。首先，在電子連續(xù)性方程的基礎(chǔ)上，構(gòu)建電弧的電氣特性模型；其次，根據(jù)弧柱磁壓縮力平衡方程和瞬時能量平衡方程分別構(gòu)建電弧的壓強、溫度等物理參數(shù)模型，推導(dǎo)得到電弧內(nèi)部的電氣特性與物理參數(shù)之間的耦合機理，得到兩者的耦合模型；然后，基于電弧爐變壓器二次側(cè)額定電流和最大電弧溫度之間的關(guān)系，利用麻雀搜索算法對電弧溫度曲線進行辨識，確定物理模型的各個參數(shù)；最后，依據(jù)50 t、55 t煉鋼交流電弧爐的現(xiàn)場實測數(shù)據(jù)，在Matlab/Simulink平臺上搭建模型并驗證模型的正確性。仿真結(jié)果表明，所提模型更符合交流電弧爐實際運行過程中的物理本質(zhì)，能更精確地反映其諧波、電壓波動特性。

交流電弧爐 電弧溫度 麻雀算法 電氣特性 物理參數(shù)

0 引言

電弧爐以其較低的投資費用、較優(yōu)的冶煉質(zhì)量以及較低的碳排放等優(yōu)勢，在國內(nèi)外冶金企業(yè)中得到廣泛的應(yīng)用[1]。但電弧爐的冶煉過程中隨機性強，具有強時變非線性的特點[2]，產(chǎn)生的大量諧波會加劇電力系統(tǒng)設(shè)備的熱效應(yīng)，降低其使用壽命；電壓偏差及波動會使得變電站的保護頻繁啟動，造成安全隱患[3-4]。因此，交流電弧爐的精確建模對研究電弧爐在電網(wǎng)產(chǎn)生的電能質(zhì)量影響，制定相關(guān)治理措施具有重要意義。

電弧爐的非線性主要是由于電弧的隨機性導(dǎo)致。隨機理論是表征電弧隨機性的有效建模方法之一。文獻[5-7]將帶通高斯白噪聲或帶通白噪聲疊加在電弧半徑、電弧電壓或電弧弧長等靜態(tài)電弧參數(shù)上，以此得到電弧爐冶煉過程中的非線性時變特征，但隨機理論僅是機械地在靜態(tài)電弧參數(shù)上疊加隨機信號，難以解釋其內(nèi)在的機理過程以及外在特性。后續(xù)研究表明，電弧爐的電壓、電流波動與混沌系統(tǒng)的外在表現(xiàn)較為接近，應(yīng)用混沌理論建立的模型比隨機理論更接近電弧現(xiàn)象的本質(zhì)[8]。文獻[9]首次使用混沌動力學(xué)描述電弧爐的非線性特性，在此基礎(chǔ)上，文獻[10-13]采用蔡氏電路產(chǎn)生的低頻混沌電路信號對電弧半徑或電弧電壓進行調(diào)制，得到具有雙吸引子混沌現(xiàn)象的電弧爐模型，但該方法需要大量的電壓電流實測數(shù)據(jù)進行參數(shù)辨識，工程適用性較差。文獻[14]將實際的電弧爐特性曲線分段線性化來模擬交流電弧爐的波動特性，但該擬合方法較為粗糙，忽略了電弧的熱慣性。文獻[15]基于電弧的物理特性進行了理論上的初步探討，在考慮電弧弧柱熱慣性的基礎(chǔ)上將弧柱電阻等效為線性時變電阻，但該方法將電弧溫度視作調(diào)制的倍頻余弦信號，難以準(zhǔn)確反映溫度對電弧弧柱的影響。文獻[16]基于電子連續(xù)性方程構(gòu)建電弧的確定性模型，并考慮電弧半徑、弧長對電子密度的影響，得到交流電弧爐外部特性與內(nèi)部微觀特性有機結(jié)合的通用模型，該模型雖在電弧微觀上有一定深入，但其將溫度、壓強等參數(shù)視作常數(shù)，沒有考慮電弧參數(shù)變化對其物理過程的影響，難以準(zhǔn)確反映冶煉過程中電弧的時變特性與其電熱物理本質(zhì)。

現(xiàn)有研究表明，電弧的物理狀態(tài)不是靜止的，電弧溫度與電弧電流有密切聯(lián)系，其會隨著電流的變化呈現(xiàn)典型的倍頻交變特性，電弧的壓強大小也會隨著電流以及弧長的變化呈現(xiàn)相關(guān)性特征[17]，同時，電弧溫度與壓強的變化也會反過來影響電子密度進而反映到電弧電壓、電流上[18]，即電弧的電氣量與物理參數(shù)之間呈現(xiàn)出耦合影響[19]。上述模型主要將電弧視作非線性電氣元件，便于和電弧爐電氣系統(tǒng)聯(lián)合求解。但是難以反映電弧內(nèi)部粒子之間的微觀動態(tài)形成機理，不能完全描述電弧的時變特征，割裂了電弧弧柱中電壓、電流與電弧溫度、壓強等參數(shù)之間的耦合影響，與實際有一定偏離，導(dǎo)致構(gòu)建的仿真模型未能準(zhǔn)確反映其電能質(zhì)量特性。

本文中，將電弧的電壓、電流定義為電氣特性，將表征電弧狀態(tài)的溫度、壓強定義為物理參數(shù)，并針對以上問題，提出一種基于電氣特性-物理參數(shù)耦合的交流電弧爐模型，揭示了電弧爐冶煉過程中的電壓、電流等電氣特性與電弧的溫度、壓強等物理參數(shù)的動態(tài)耦合關(guān)系。首先，基于電弧的電子連續(xù)性方程構(gòu)建電氣特性模型，考慮電弧爐冶煉過程中電氣量對物理過程中電弧溫度、壓強等物理參數(shù)影響的機理；其次，計及電弧的磁壓縮與瞬時能量平衡過程中電弧溫度、壓強的變化對電弧電氣特性的影響，構(gòu)建電氣特性-物理參數(shù)耦合的電弧爐模型；進一步地，基于實際電弧爐變壓器參數(shù)構(gòu)建電弧溫度曲線，并使用麻雀搜索算法對電弧溫度進行辨識，以確定模型參數(shù)；最后，以福建某50 t、55 t交流電弧爐實測數(shù)據(jù)為例，在Matlab/Simulink軟件中，驗證了模型的正確性。

1 電弧爐電氣特性-物理參數(shù)耦合機理

本文所提出的交流電弧爐通用模型由電氣特性模型與物理參數(shù)模型兩部分構(gòu)成。電氣特性模型通過電弧內(nèi)部微觀特性來表征電弧爐的強時變非線性，與此同時，電弧的電氣量會對電弧的溫度、壓強等物理參數(shù)產(chǎn)生影響，構(gòu)建物理參數(shù)模型得到冶煉過程中電弧溫度、壓強的累積效應(yīng)反映到電弧電氣量的變化特性，考慮兩者的密切聯(lián)系與相互影響，從而建立電弧電氣特性與物理參數(shù)之間的耦合模型。

1.1 電氣特性模型

對于任何電弧等離子體，產(chǎn)生的電流可由來自電子的電流和來自離子的電流兩個基本項組成[18]。而在電弧放電過程中，電子的遷移速度遠大于離子的遷移速度。因此，在電中性條件下[20]，弧柱電導(dǎo)率可近似認(rèn)為主要由電子密度和電子遷移率決定。

為研究電弧放電的最基本特性，可忽略電子和離子的擴散，利用電子連續(xù)性方程求解電子密度[16]，即

根據(jù)歐姆定律的微分形式，可以得到

式中，為電弧弧柱的平均電流密度；為電弧弧長；為電弧電流；為電弧半徑，其變化規(guī)律可表示為

式中，c為陰極斑點半徑；為電弧的軸向坐標(biāo)。

對電弧弧柱中的電場強度軸向方向進行積分得到電弧電壓為

聯(lián)立式（1）～式（8）可得到電弧的電氣特性方程為

現(xiàn)有研究表明，電弧爐的弧長在實際運行過程中會呈現(xiàn)出周期性、準(zhǔn)周期性或隨機性的波動[16]，可用式（11）進行模擬。

1.2 物理參數(shù)模型

電弧現(xiàn)象在本質(zhì)上是電物理與熱物理的一個綜合過程，并且在大多情況下，熱物理過程起著決定性的作用[22]。基于瞬時能量平衡方程[20]可表征出交流電弧放電過程中的物理特性，該式適用于電弧弧柱的所有體積元。

在電弧環(huán)境中，電弧弧柱的磁場會對弧柱產(chǎn)生一個壓縮力，導(dǎo)致弧柱半徑收縮，從而影響電弧電阻。表明電弧屬于彈性導(dǎo)體，其形狀會受到環(huán)境壓強的影響[17]。對弧柱構(gòu)建力平衡方程以研究壓強對電弧的影響[23]。

依據(jù)安培環(huán)路定律，可對磁感應(yīng)強度進行求 解，有

將式（14）代入式（13），并對式（13）兩邊從弧柱邊緣向弧柱內(nèi)部積分，可得到磁壓縮對弧柱徑向距離處的壓強為

式中，()為弧柱內(nèi)徑向距離處的壓強；0為弧柱邊緣處壓強。

將式（7）的平均電流密度代入式（15），整理得到磁壓縮附加壓強與電弧電流、弧長的關(guān)系為

當(dāng)徑向距離為零時，即可得到弧軸處的附加壓強為

對于大電流電弧，對流引起的能量損失占電弧能量平衡的5%左右，其散失能量的90%是通過輻射散熱完成，這部分能量對電弧爐作業(yè)具有支配作用[24-26]。因此，本文關(guān)注最主要的輻射散熱為

式中，1、2為待確定參數(shù)，其大小與電弧爐的容量有關(guān)；為玻耳茲曼常數(shù)；m為受激原子的高能級平均值。

聯(lián)立式（7）、式（12）～式（18），得到交流電弧的物理參數(shù)方程為

1.3 耦合模型

電子連續(xù)性方程體現(xiàn)的電場特性與歐姆定律組成電弧的電氣特性模型，瞬時能量平衡方程與磁場產(chǎn)生的磁壓縮力平衡方程構(gòu)成電弧的物理參數(shù)模型，其中的能量轉(zhuǎn)移會非常強烈地影響電場和磁場，而電磁場又通過電弧電導(dǎo)進行加熱，作為能量源反過來影響溫度場[19]。

電弧的電氣特性模型與物理參數(shù)模型較為完整地反映了等離子體中溫度場、電場以及磁場之間重要耦合效應(yīng)，正是這些耦合效應(yīng)導(dǎo)致了電弧的形成。圖1給出了電氣特性模型與物理參數(shù)模型之間的耦合關(guān)系。

1.4 邊界條件

模型建立后，為了獲得數(shù)值解，必須給定邊界條件。現(xiàn)有研究表明，電弧在正常工作狀態(tài)下呈現(xiàn)為典型的“長鐘型（鐘罩型）”[27]，如圖2所示。因此本文模型的邊界條件建立在以下計算區(qū)域。圖2中，為等離子體電弧在空氣中的導(dǎo)電區(qū)域，其中為陰極斑區(qū)域；為電??；為陽極表面。為了簡化計算，認(rèn)為正常工作狀態(tài)下電弧是穩(wěn)定的，維持為長鐘型，且無偏弧情況。

由此可得到計算區(qū)域的邊界條件見表1[28]。

表1 電弧等離子體的邊界條件

Tab.1 Boundary conditions of arc plasma

2 電弧爐仿真模型及參數(shù)確定

2.1 電弧爐仿真模型

本文所提出的電氣特性-物理參數(shù)耦合模型以文獻[16]的電子連續(xù)性方程為基礎(chǔ)構(gòu)建電氣特性模型，以弧長變化模型來表征電弧的不確定性，以確定性模型表征電弧的非線性，如圖3所示；以電弧的磁壓縮力平衡方程和瞬時能量平衡方程構(gòu)建物理參數(shù)模型，以此反映溫度、壓強對電弧確定性模型的影響，如圖4所示。

圖3 電氣特性模型

圖4 物理參數(shù)模型

2.2 基于麻雀搜索算法的參數(shù)辨識

耦合模型中包含電子密度、電弧電壓、電弧電流、溫度、弧長和反映熱慣性影響等因素的多個參數(shù)，若想準(zhǔn)確確定這些參數(shù)，需要大量的測試數(shù)據(jù)，由于電弧爐工作的特殊性和已有測試手段的局限性，難以準(zhǔn)確獲取模型的所有參數(shù)。因此，只能從工程應(yīng)用的角度對模型中的各個參數(shù)進行近似估算及辨識，以確定電弧爐模型的各個參數(shù)。

本文使用的電氣特性模型參數(shù)可參照文獻[16]進行估算，而在物理參數(shù)模型中的溫度場參數(shù)、1和2還未確定。針對待確定的參數(shù)，本文使用麻雀搜索算法[29]（Sparrow Search Algorithm, SSA）對電弧溫度進行參數(shù)辨識。麻雀搜索算法（SSA）是一種新的元啟發(fā)式算法，相較于目前的灰狼優(yōu)化算法、粒子群算法、引力搜索算法等[30-34]其他算法，具有更強的逃離局部極值能力和全局搜索能力。它模擬了麻雀群體在抗捕食的過程中，為獲取食物個體位置不斷逼近食物的過程。發(fā)現(xiàn)者位置的更新方式可表示為

對于麻雀中的加入者，其位置更新方式為

處于群體邊緣處的麻雀作為警戒者，其位置更新方式為

現(xiàn)有文獻證明，在不同大小的電流下，弧柱溫度的變化曲線是類似的[26]。因此本文基于Bowman在電流為8 A、15 A、25 A時對碳電極交流電弧溫度測量得到的若干結(jié)果[32]，得到真實情況下的電弧溫度曲線特征，結(jié)合實際電弧爐變壓器二次側(cè)的額定電流與文獻[28]中不同電弧電流大小下電弧最大溫度的計算結(jié)果，構(gòu)建符合電弧爐實際工況下的電弧溫度曲線，并使用SSA算法對所構(gòu)建的電弧溫度曲線進行辨識，以使得模型的電弧溫度仿真計算結(jié)果與實際情況相近。本文使用文獻[15]所提出的余弦函數(shù)來近似設(shè)定需要逼近的電弧溫度曲線。

電弧電流與最大電弧溫度H之間的對應(yīng)關(guān)系見表2[28]，基于麻雀算法的參數(shù)辨識流程如圖5所示。

表2 最大電弧溫度與電流的對應(yīng)關(guān)系

Tab.2 The relationship between maximum arc temperature and current

3 實例分析

本文以福建的50 t和55 t煉鋼電弧爐實測數(shù)據(jù)進行算例驗證，以驗證模型的正確性與適用性。鑒于電弧爐僅在熔化期具有嚴(yán)重的電能質(zhì)量問題且對電弧爐設(shè)備的考核也多在熔化期進行，因此本文算例重點針對電弧爐熔化期進行驗證。且出于設(shè)備穩(wěn)定生產(chǎn)要求以及人員安全性考慮，一般難以對爐變低壓進行測試，因此本文使用電弧爐變壓器高壓側(cè)數(shù)據(jù)并基于變壓器諧波傳遞特性將其折算至低壓側(cè)[36-37]。

圖5 基于麻雀算法的電弧溫度辨識流程

3.1 50 t煉鋼交流電弧爐仿真分析

本文以福建某50 t煉鋼交流電弧爐負荷為例，圖6為其供電系統(tǒng)的等效電路。鋼廠的主變壓器T1將110 kV高壓電網(wǎng)電壓轉(zhuǎn)換為35 kV；T2為電弧爐變壓器，其設(shè)備型號為HSSPZ-22000/35；d為短網(wǎng)阻抗，d=0.4 mW；d為短網(wǎng)的電抗（自感和互感之和），d=3 mW。電弧爐變壓器的參數(shù)見表3。

圖6 50 t交流電弧爐供電系統(tǒng)等效電路

表3 HSSPZ-22000/35電弧爐變壓器參數(shù)

Tab.3 Arc furnace transformer parameters of HSSPZ-22000/35

圖7 基于SSA的電弧溫度辨識曲線

圖8 50 t交流電弧爐仿真電路

在電弧爐實際運行過程中，由于電弧的隨機性，電弧爐的每個交流電周期的伏安特性曲線并不一致，圖9為仿真得到的50 t煉鋼交流電弧爐的特性曲線。從圖9中可以看出，考慮了電弧電流、弧長變化以及物理參數(shù)的耦合影響后，所建立的模型能夠反映電弧電流和電弧電壓之間呈現(xiàn)的高度強時變非線性特征，與電弧爐實際冶煉過程中的特性基本一致。

圖9 50 t交流電弧爐U-I特性曲線

圖10為仿真所得的50 t交流電弧爐的電弧電壓和電弧電流。由圖10可以看出，由于電弧溫度的交變特性以及壓強的變化影響，電弧電流和電弧電壓呈現(xiàn)明顯的畸變特征，并且電弧電壓還存在明顯的隨機波動。綜上所述，本文所建立的電弧爐模型能夠更切合實際地體現(xiàn)電弧爐冶煉過程中的高度非線性和隨機特性。

圖10 50 t交流電弧爐仿真結(jié)果

1）諧波含量的對比分析

50 t交流電弧爐的實際電弧電壓、電流數(shù)據(jù)通過故障錄波器（采樣頻率12.6 kHz）采集得到。并在此基礎(chǔ)上，利用Matlab/Simulink軟件分別搭建文獻[16]模型和本文方法模型，并使用軟件的快速傅里葉變換（Fast Fourier Transform, FFT）分析功能對電弧爐電弧電壓及電弧電流的各次諧波含量進行計算，表4及表5給出了兩種方法的仿真諧波含量與實測的對比結(jié)果。

表4 50 t交流電弧爐的諧波電壓分析結(jié)果

Tab.4 Harmonic voltage of 50 t AC electric arc furnace

表5 50 t交流電弧爐的諧波電流分析結(jié)果

Tab.5 Harmonics current of 50 t AC electric arc furnace

根據(jù)兩種方法的仿真結(jié)果與實測的對比，從表4和表5可知，本文所提出的交流電弧爐模型具有更高的精度。各次的諧波電壓畸變率誤差值不超過0.55%，諧波電壓總畸變率誤差值僅為0.15%。各次諧波電流畸變率的仿真精度在整體上更具優(yōu)勢，諧波電流總畸變率誤差值不超過0.17%。

2）電壓波動對比及分析

根據(jù)GB/T 12326-2008《電能質(zhì)量 電壓波動和閃變》，電壓波動定義為電壓方均根曲線上相鄰兩個極值電壓之差。根據(jù)以上定義，電壓波動的計算公式為

根據(jù)式（24）將仿真模型中電壓波動狀況與實測測量值進行對比，見表6。

表6 50 t交流電弧爐電壓波動狀況

Tab.6 Voltage fluctuation of 50 t AC electric arc furnace

由表6可知，本文方法與文獻[16]方法對電壓波動的仿真精度基本一致，但本文所建立的模型對諧波和電壓波動的仿真效果上誤差均較小。

圖11為電弧電流及電弧溫度的仿真結(jié)果，可以看出：交流電弧的溫度隨電流的變化而變化，呈現(xiàn)交變特性，且為二倍頻的關(guān)系。由于電弧的熱慣性，當(dāng)電弧電流達到最大值max時，電弧溫度滯后一定時間才達到最大值max。圖12為電弧壓強的仿真結(jié)果，可以看出：電弧弧柱的壓強并不是一直處于一個大氣壓下，而是會在1～1.7個大氣壓范圍內(nèi)波動。

圖11 電弧電流及電弧溫度仿真結(jié)果

3.2 55 t煉鋼交流電弧爐仿真分析

為進一步驗證模型在不同型號電弧爐的適用性，本文使用福建某五金鑄造的55 t煉鋼交流電弧爐為例進行驗證。圖13為該負荷的供電方式，電弧爐變壓器型號為HSSPZ-38000/35，d=0.4 mW，d= 3 mW，電弧爐變壓器的參數(shù)見表7。鑒于該負荷加裝了靜止無功補償器（Static Var Generator, SVG）治理裝置，在此僅驗證諧波電流的仿真效果。

圖12 電弧壓強變化曲線仿真結(jié)果

圖13 55 t交流電弧爐供電示意圖

表7 HSSPZ-38000/35電弧爐變壓器參數(shù)

Tab.7 Arc furnace transformer parameters of HSSPZ-38000/35

仿真所得55 t交流電弧爐的諧波電流曲線如圖15所示。其中，實際電弧電流通過德維創(chuàng)電能質(zhì)量分析儀（采樣頻率為20 kHz）采集得到。根據(jù)文獻[16]方法和本文方法構(gòu)建55 t交流電弧爐模型，分別對實測和仿真的諧波電流各次含有率及總畸變率進行計算，對比結(jié)果見表8。

圖14 55 t煉鋼交流電弧爐仿真電路

圖15 55 t交流電弧爐電弧電流仿真結(jié)果

表8 55 t交流電弧爐的電流諧波狀況

Tab.8 Current harmonics of 55 t AC electric arc furnace

根據(jù)表8結(jié)果可知，本文構(gòu)建的模型可適應(yīng)不同型號容量的電弧爐仿真，并且使用本文方法在仿真精度上具有更小的誤差。

4 結(jié)論

本文提出一種基于電弧的電氣特性與物理參數(shù)耦合的交流電弧爐模型，并通過50 t及55 t煉鋼交流電弧爐實測數(shù)據(jù)進行驗證，得到以下結(jié)論：

1）揭示了電弧微觀特性下壓強、電弧溫度等物理參數(shù)與電壓、電流等電氣特性的耦合機理，充分考慮了電弧的熱慣性，在此基礎(chǔ)上構(gòu)建的電弧爐模型，更貼近交流電弧爐實際運行過程中的特征，對電弧爐的諧波電壓、電流以及電壓波動等電能質(zhì)量特性的仿真精度有更進一步的提高，且能應(yīng)用于不同型號電弧爐的仿真研究，具有較好的適用性。

2）本文方法在確定模型時僅需爐變低壓側(cè)額定電流參數(shù)以及爐變高壓側(cè)數(shù)周波錄波電流，即可確定模型的全部參數(shù)，便于從工程角度確定電弧溫度，具有較強的工程實用性和可移植性。

[1] 朱榮, 魏光升, 張洪金. 近零碳排電弧爐煉鋼工藝技術(shù)研究及展望[J]. 鋼鐵, 2022, 57(10): 1-9.

Zhu Rong, Wei Guangsheng, Zhang Hongjin. Research and prospect of EAF steelmaking with near-zero carbon emissions[J]. Iron and Steel, 2022, 57(10): 1-9.

[2] Collantes Bellido R, Gomez T, Identification and modelling of a three phase arc furnace for voltage disturbance simulation[J].IEEE Transactions on Power Delivery, 1997, 12(4): 1812-1817.

[3] 林海雪. 電弧爐的有功功率沖擊對發(fā)電機組的影響[J]. 中國電機工程學(xué)報, 2014, 34(增刊1): 232-238.

Lin Haixue. Influence of active power impact caused by EAF on the generating sets[J]. Proceedings of the CSEE, 2014, 34(S1): 232-238.

[4] 林才華, 張逸, 黃道姍, 等. 低頻間諧波檢測分析與治理技術(shù)研究綜述[J]. 供用電, 2019, 36(9): 59-65.

Lin Caihua, Zhang Yi, Huang Daoshang, et al. A summary of research on low frequency inter-harmonic detection analysis and governance technology[J]. Distribution and Utilization, 2019, 36(9): 59-65.

[5] 王晶, 束洪春, 林敏, 等. 用于動態(tài)電能質(zhì)量分析的交流電弧爐的建模與仿真[J]. 電工技術(shù)學(xué)報, 2003, 18(3): 53-58.

Wang Jing, Shu Hongchun, Lin Ming, et al. Modeling and simulation of AC arc furnace for dynamic power quality studies[J]. Transactions of China Electro- technical Society, 2003, 18(3): 53-58.

[6] 王育飛, 潘艷霞, 姜建國. 基于MATLAB的交流電弧爐隨機模型與仿真[J]. 高電壓技術(shù), 2008, 34(5): 973-977.

Wang Yufei, Pan Yanxia, Jiang Jianguo. Stochastic model of AC electric arc furnace based on MATLAB[J]. High Voltage Engineering, 2008, 34(5): 973-977.

[7] 王琰, 毛志忠, 李妍, 等. 用于電壓波動研究的交流電弧爐電弧模型[J]. 電網(wǎng)技術(shù), 2010, 34(1): 36-40.

Wang Yan, Mao Zhizhong, Li Yan, et al. An electric arc model of AC electric arc furnace for research on voltage fluctuation[J]. Power System Technology, 2010, 34(1): 36-40.

[8] KingP E, Ochs T L, Hartman A D. Chaotic responses in electric arc furnaces[J]. Journal of Applied Physics, 1994, 76(4): 2059-2065.

[9] O'Neill-Carrillo E, Heydt GT, Kostelich EJ, et al. Nonlinear deterministic modeling of highly varying loads[J]. IEEE Transactions on Power Delivery, 1999, 14(2): 537-542.

[10] 寧元中, 梁穎, 吳昊. 電弧爐的混合仿真模型[J]. 四川大學(xué)學(xué)報(工程科學(xué)版), 2005, 37(1): 85-89.

Ning Yuanzhong, Liang Ying, Wu Hao. Circuit model of arc furnace for mixed-discipline simulation[J]. Advanced Engineering Sciences, 2005, 37(1): 85-89.

[11] Ozgun O, Abur A. Flicker study using a novel arc furnace model[J]. IEEE Transactions on Power Delivery, 2002, 17(4): 1158-1163.

[12] 王育飛, 姜建國. 用于電能質(zhì)量研究的新型交流電弧爐混沌模型[J]. 中國電機工程學(xué)報, 2008, 28(10): 106-110.

Wang Yufei, Jiang Jianguo. A novel chaotic model of AC electric arc furnace for power quality study[J]. Proceedings of the CSEE, 2008, 28(10): 106-110.

[13] 林才華, 張逸, 邵振國, 等. 用于低頻非平穩(wěn)間諧波研究的超高功率電弧爐模型[J]. 中國電力, 2020, 53(11): 1-8.

Lin Caihua, Zhang Yi, Shao Zhenguo, et al. An ultra-high-power electric arc furnace model for low- frequency non-stationary inter-harmonics studies[J]. Electric Power, 2020, 53(11): 1-8.

[14] 祁碧茹, 肖湘寧. 用于電壓波動研究的電弧爐的模型和仿真[J]. 電工技術(shù)學(xué)報, 2000, 15(3): 31-35.

Qi Biru, Xiao Xiangning. Modeling and simulation of an arc furnace for voltage fluctuation investigation[J]. Transactions of China Electrotechnical Society, 2000, 15(3): 31-35.

[15] 劉小河, 程少庚, 蘇文成. 電弧爐電氣系統(tǒng)的諧波分析研究[J]. 電工技術(shù)學(xué)報, 1994, 9(2): 21-26, 10.

Liu Xiaohe, Cheng Shaogeng, Su Wencheng. Harmonic analysis of electric systems of arc furnace[J]. Transactions of China Electrotechnical Society, 1994, 9(2): 21-26, 10.

[16] 張逸, 林才華, 邵振國, 等. 基于電子連續(xù)性方程的新型交流電弧爐通用模型[J]. 中國電機工程學(xué)報, 2021, 41(21): 7425-7434.

Zhang Yi, Lin Caihua, Shao Zhenguo, et al. A novel general model of AC electric arc furnace based on electron continuity equation[J]. Proceedings of the CSEE, 2021, 41(21): 7425-7434.

[17] 曹振. 電弧爐電弧模型的研究[D]. 遼寧: 東北大學(xué), 2009.

[18] Underwood T, Roy S, Glaz B. Physics based lumped element circuit model for nanosecond pulsed dielectric barrier discharges[J]. Journal of Applied Physics, 2013, 113(8): 083301.

[19] 徐龍, 李秋菊, 馬文會, 等. 三維非軸對稱直流電弧射流過程數(shù)值模擬[J]. 昆明理工大學(xué)學(xué)報(自然科學(xué)版), 2021, 46(4): 18-26.

Xu Long, Li Qiuju, Ma Wenhui, et al. Numerical simulation of 3D non-axisymmetric DC arc jet process[J]. Journal of Kunming University of Science and Technology(Natural Sciences), 2021, 46(4): 18-26.

[20] 過增元, 趙文華. 電弧和熱等離子體[M]. 北京: 科學(xué)出版社, 1986.

[21] Lowke J J. Simple theory of free-burning arcs[J]. Journal of Physics D: Applied Physics, 1979, 12(11): 1873-1886.

[22] 趙永秀, 劉樹林, 馬一博. 爆炸性試驗電感電路分?jǐn)喾烹娞匦苑治雠c建模[J]. 煤炭學(xué)報, 2015, 40(7): 1698-1704.

Zhao Yongxiu, Liu Shulin, Ma Yibo. Analysis and modeling of inductor-disconnected-discharged charac- teristics based on explosive test[J]. Journal of China Coal Society, 2015, 40(7): 1698-1704.

[23] Abid F, Niayesh K, Jonsson E, et al. Arc voltage characteristics in ultrahigh-pressure nitrogen including supercritical region[J]. IEEE Transactions on Plasma Science, 2018, 46(1): 187-193.

[24] Latham D J. A channel model for long arcs in air[J]. The Physics of Fluids, 1980, 23(8): 1710.

[25] 于豐, 毛志忠. 電弧爐系統(tǒng)電弧模型的建立與參數(shù)辨識[J]. 東北大學(xué)學(xué)報(自然科學(xué)版), 2013, 34(2): 178-181.

Yu Feng, Ma Zhizhong. An arc model for arc furnace system and parameter identification[J]. Journal of Northeastern University (Natural Science), 2013, 34(2): 178-181.

[26] 王其平, 電器電弧理論[M]. 北京: 機械工業(yè)出版社, 1991.

[27] 姚聰林, 朱紅春, 姜周華, 等. 電弧爐內(nèi)長電弧等離子體的數(shù)值模擬[J]. 工程科學(xué)學(xué)報, 2020, 42(增刊1): 60-67.

Yao Conglin, Zhu Hongchun, Jiang Zhouhua, et al. Numerical simulation of a long arc plasma in an electric arc furnace[J]. Chinese Journal of Engin- eering, 2020, 42(S1): 60-67.

[28] 王豐華. 電弧爐建模研究及其應(yīng)用[D]. 上海: 上海交通大學(xué), 2006.

[29] Xue J K, Shen B. A novel swarm intelligence optimization approach: sparrow search algorithm[J]. Systems Science and Control Engineering, 2020, 8(1): 22-34.

[30] Zhang C L, Ding S F. A stochastic configuration network based on chaotic sparrow search algorithm[J]. Knowledge-Based Systems, 2020, 8(1): 22-34.

[31] 薛建凱. 一種新型的群智能優(yōu)化技術(shù)的研究與應(yīng)用[D]. 上海: 東華大學(xué), 2020.

[32] 李婕, 楊淑英, 謝震, 等. 基于有效信息迭代快速粒子群優(yōu)化算法的永磁同步電機參數(shù)在線辨識[J]. 電工技術(shù)學(xué)報, 2022, 37(18): 4604-4613.

Li Jie, Yang Shuying, Xie Zhen, et al. Online para- meter identification of permanent magnet syn- chronous motor based on fast particle swarm optimization algorithm with effective information iterated[J]. Transactions of China Electrotechnical Society, 2022, 37(18): 4604-4613.

[33] 劉廷章, 張子恒. 基于循環(huán)迭代粒子群算法對電動汽車無線充電系統(tǒng)的聯(lián)合參數(shù)在線辨識[J]. 電工技術(shù)學(xué)報, 2022, 37(18): 4548-4564.

Liu Tingzhang, Zhang Ziheng. Online identification of joint parameters of electric vehicle wireless charging system based on loop-iteration particle swarm optimization[J]. Transactions of China Elec- trotechnical Society, 2022, 37(18): 4548-4564.

[34] 葛磊蛟, 趙康, 孫永輝, 等. 基于孿生網(wǎng)絡(luò)和長短時記憶網(wǎng)絡(luò)結(jié)合的配電網(wǎng)短期負荷預(yù)測[J]. 電力系統(tǒng)自動化, 2021, 45(23): 41-50.

Ge Leijiao, Zhao Kang, Sun Yonghui, et al. Short- term load forecasting of distribution network based on combination of siamese network and long short-term memory network[J]. Automation of Electric Power System, 2021, 45(23): 41-50.

[35] Bowman B, Edels H. Radial temperature measure- ments of alternating current arcs[J]. Journal of Physics D: Applied Physics, 1969, 2(1): 53-63.

[36] 賈靜然, 段曉波, 盧錦玲, 等. 變壓器的諧波傳遞特性研究[J]. 中國電力, 2017, 50(1): 92-96, 157.

Jia Jingran, Duan Xiaobo, Lu Jinling, et al. Research on harmonic transfer characteristics of transformer[J]. Electric Power, 2017, 50(1): 92-96, 157.

[37] 胡雪凱, 張乾, 田廣, 等. 雙繞組變壓器諧波傳遞特性研究[J]. 河北電力技術(shù), 2021, 40(3): 62-66,71.

Hu Xuekai, Zhang Qian, Tian Guang, et al. Research on harmonic transfer characteristics of double winding transformer[J]. Heibei Electric Power, 2021, 40(3): 62-66, 71.

AC Arc Furnace Model Based on Coupling of Electrical and Physical Parameters

1123

（1. School of Electrical Engineering and Automation Fuzhou University Fuzhou 350108 China 2. Fujian Electrical Power Research Institute Fuzhou 350007 China 3. State Grid Ningde Electric Power Supply Co. Ltd Ningde 352100 China）

Accurate modeling of electric arc furnaces (EAF) is important for studying the impacton power quality of EAF and developing relevant management measures in the power grid. However, the existing research mainly treats the electric arc as a nonlinear electrical component, neglecting the coupling of electric and physical parameters in the electric arc. Thus, reflecting the electro thermal physical nature of the arc is challenging, resulting in a deficiency in simulation accuracy.

This paper proposes an improved AC arc furnace model. The voltage and current of the arc are defined as electrical characteristics, while the temperature and pressure that characterize the arc state are defined as physical parameters. The coupling mechanism between the electrical and physical parameters is determined based on the electronic continuity equation, Ohm's Law, the magnetic compression force balance equation, and the instantaneous energy balance equation. A general model is constructed that reflects harmonics and voltage fluctuations of AC electric arc furnaces more accurately. Its mathematical expressions are given as follows:

Measured data from two steel mills’ 50t and 55t electric arc furnaces are used to verify the model. The results indicate that the simulated AC arc temperature varies with the arc current, exhibiting an alternating characteristic and a double-frequency relationship. The thermal inertia of the arc induces a lad in the peak temperature relative to the peak current. The pressure in the arc column fluctuates around atmospheric pressure. Furthermore, considering the arc’s electrical-physical coupling characteristics, the simulation accuracy of the electrical characteristics, such as harmonic voltage, harmonic current, and voltage fluctuations, has been further improved compared to existing methods.

AC electric arc furnace, arc temperature, sparrow search algorithm, electrical characteristic, physical parameters

10.19595/j.cnki.1000-6753.tces.230034

TF748.41

國家自然科學(xué)基金資助項目（51777035）。

2023-01-09

2023-07-29

徐云聰 男，1999年生，碩士研究生，研究方向為電能質(zhì)量分析。E-mail: 1360397798@qq.com

張 逸 男，1984年生，博士，副教授，研究方向為電能質(zhì)量、主動配電網(wǎng)以及電力數(shù)據(jù)分析等。E-mail: zhangyi@fzu.edu.cn（通信作者）

（編輯 郭麗軍）