同軸電極的負(fù)電暈特里切爾脈沖特性分析

何旺齡　何俊佳　張　錦　萬(wàn)保權(quán)

(1.華中科技大學(xué)電氣與電子工程學(xué)院　武漢　430074 2.中國(guó)電力科學(xué)研究院(武漢)　武漢　430074)

同軸電極的負(fù)電暈特里切爾脈沖特性分析

何旺齡1何俊佳1張錦2萬(wàn)保權(quán)2

(1.華中科技大學(xué)電氣與電子工程學(xué)院武漢430074 2.中國(guó)電力科學(xué)研究院(武漢)武漢430074)

摘要利用流體動(dòng)力學(xué)模型，對(duì)大氣壓條件下空氣間隙的同軸電極結(jié)構(gòu)進(jìn)行了負(fù)極性電暈放電特性研究，獲得了一組規(guī)則的特里切爾脈沖。除第一個(gè)脈沖外，其他脈沖的幅值均在1.5～2.5 mA。然后，對(duì)第一個(gè)脈沖(稱為第一脈沖)不同發(fā)展階段的空間電場(chǎng)分布和粒子密度分布進(jìn)行了描述，并討論了第一脈沖與第二脈沖之間的死區(qū)時(shí)間內(nèi)不同時(shí)刻負(fù)離子密度的變化過(guò)程及其對(duì)陰極表面電場(chǎng)的影響作用，最后，對(duì)第二脈沖發(fā)展前后空間電場(chǎng)分布、電子密度分布和負(fù)離子密度分布及其變化規(guī)律進(jìn)行了進(jìn)一步的描述。仿真表明，隨著負(fù)離子在電場(chǎng)力的作用下逐漸向外遷移，導(dǎo)線表面電場(chǎng)強(qiáng)度開(kāi)始恢復(fù)，但其恢復(fù)速度在負(fù)離子遠(yuǎn)離導(dǎo)線的過(guò)程中逐漸變緩。同時(shí)，當(dāng)導(dǎo)線表面電場(chǎng)強(qiáng)度恢復(fù)至低于初始拉普拉斯場(chǎng)強(qiáng)的某一特定值時(shí)，新的脈沖就開(kāi)始發(fā)生，第二脈沖發(fā)生時(shí)空間電場(chǎng)、電子密度及負(fù)離子密度的最大值均小于第一脈沖。

關(guān)鍵詞：負(fù)電暈特里切爾脈沖空間電場(chǎng)分布粒子密度

0引言

我國(guó)正在大力發(fā)展特高壓交直流輸電線路[1，2]，隨著輸電線路電壓等級(jí)的提高，導(dǎo)線周圍的電暈放電現(xiàn)象變得愈加明顯，強(qiáng)烈的電暈放電現(xiàn)象會(huì)引起線路的絕緣劣化和電暈損失增大等問(wèn)題，影響線路的輸送安全和輸送效率[3-6]。同時(shí)，也會(huì)導(dǎo)致一些環(huán)境效應(yīng)，如電磁干擾和可聽(tīng)噪聲等[7-9]。此外，電暈放電也被人們用來(lái)進(jìn)行空氣凈化和污水處理等工業(yè)應(yīng)用[10]。而特里切爾脈沖放電是電暈放電的重要組成部分，對(duì)解釋電暈放電現(xiàn)象起著至關(guān)重要的作用。因此，有必要對(duì)負(fù)極性特里切爾脈沖進(jìn)行分析研究。

特里切爾于1938年首次發(fā)現(xiàn)，在負(fù)極性電暈放電條件下，存在一種規(guī)則的電暈電流脈沖族，通常表現(xiàn)為上升時(shí)間較快(小到ns級(jí))，脈沖幅值較小，持續(xù)時(shí)間較短[11]。此后，人們對(duì)特里切爾脈沖進(jìn)行了大量的理論研究和實(shí)驗(yàn)研究。在實(shí)驗(yàn)研究方面，文獻(xiàn)[12]對(duì)不同氣體條件下的特里切爾脈沖進(jìn)行了測(cè)量，認(rèn)為特里切爾脈沖只存在于電負(fù)性氣體之中。文獻(xiàn)[13]在空氣條件下，采用針板電極結(jié)構(gòu)對(duì)特里切爾脈沖進(jìn)行了大量的實(shí)驗(yàn)測(cè)量研究，建立了脈沖頻率、平均電暈電流與施加電壓、針尖的曲率半徑以及針板距離之間的相互關(guān)系。文獻(xiàn)[14]發(fā)現(xiàn)在給定電極條件下，特里切爾脈沖的形成過(guò)程與電壓、氣壓以及外電路并聯(lián)電容并無(wú)關(guān)系。

在研究氣體放電機(jī)理方面，由于空氣放電等離子體測(cè)量技術(shù)的限制，許多放電特征參數(shù)還不能完全通過(guò)實(shí)驗(yàn)來(lái)獲取，所以國(guó)內(nèi)外很多學(xué)者利用數(shù)學(xué)手段對(duì)特里切爾脈沖進(jìn)行建模分析和研究。近年來(lái)，P.Sattari等[15，16]利用針板結(jié)構(gòu)，提出一種新的二維特里切爾脈沖仿真模型，討論了不同施加電壓條件下特里切爾脈沖頻率的變化規(guī)律，并與W.L.Lama等的實(shí)測(cè)數(shù)據(jù)進(jìn)行了對(duì)比分析。T.N.Tran等[17]提出了基于軸對(duì)稱的有限元模型對(duì)負(fù)極性電暈放電進(jìn)行仿真分析，對(duì)不同施加電壓及不同二次發(fā)射系數(shù)對(duì)特里切爾脈沖頻率及幅值的影響進(jìn)行了分析研究，并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比分析。尹晗等[18]建立了二維時(shí)域模型對(duì)線板結(jié)構(gòu)的特里切爾脈沖進(jìn)行仿真分析，其中利用有限元方法進(jìn)行泊松方程求解，利用有限體積法進(jìn)行對(duì)流方程求解。伍飛飛等[19，20]基于流體動(dòng)力學(xué)模型，提出了一種考慮不同化學(xué)碰撞反應(yīng)和光電離的混合數(shù)值模型，并對(duì)針板結(jié)構(gòu)下的特里切爾脈沖進(jìn)行了仿真分析。

上述學(xué)者在進(jìn)行負(fù)極性電暈放電的數(shù)值計(jì)算時(shí)，對(duì)特里切爾脈沖族中第一脈沖的微觀過(guò)程進(jìn)行了建模和仿真分析，并對(duì)整個(gè)脈沖族的宏觀特性(如脈沖幅值和脈沖頻率)進(jìn)行了討論。基于前人的研究成果，本文采用流體模型，利用有限元方法，對(duì)特里切爾脈沖發(fā)展過(guò)程中的3種帶電粒子進(jìn)行了建模仿真，討論了第一脈沖不同發(fā)展階段的空間電場(chǎng)分布和粒子密度分布，并對(duì)死區(qū)時(shí)間內(nèi)空間電場(chǎng)和負(fù)離子密度的變化規(guī)律進(jìn)行了細(xì)致的描述，最后討論了第二個(gè)脈沖起始前后電場(chǎng)、電子密度和負(fù)離子密度的變化過(guò)程。

1數(shù)值模型建立

1.1控制方程

本文基于流體動(dòng)力學(xué)模型，對(duì)負(fù)極性電暈放電過(guò)程中電子、正離子和負(fù)離子的對(duì)流擴(kuò)散運(yùn)動(dòng)及產(chǎn)生和消失過(guò)程進(jìn)行了數(shù)學(xué)建模，并與考慮空間電荷影響的泊松方程進(jìn)行耦合求解，從而模擬特里切爾脈沖的發(fā)展過(guò)程。泊松方程和連續(xù)性方程見(jiàn)式(1)～式(4)。

泊松方程

(1)

連續(xù)性方程

(2)

(3)

(4)

式中，Φ為電勢(shì)，V；t為時(shí)間，s；Ne、Np、Nn分別為電子、正離子和負(fù)離子的密度，1/m3；e為單位電荷量，1.69×10-19C；ε0為真空介電常數(shù)，8.854×10-12F/m；εr為相對(duì)介電常數(shù)，取值為1；We、Wp、Wn分別為電子、正離子和負(fù)離子的運(yùn)動(dòng)速度矢量，m/s；α為空氣在標(biāo)準(zhǔn)大氣壓條件下的湯森電離系數(shù)，1/m；η為吸附系數(shù)，1/m；βep為電子-正離子復(fù)合率；βnp為正負(fù)離子復(fù)合率；De為電子擴(kuò)散系數(shù)。

相關(guān)參數(shù)的具體計(jì)算方法如下所示[21]。

(5)

(6)

(7)

Wp=μpE，Wn=-μnE

(8)

βep=βnp=2.2×10-12m3/s

(9)

De=1.28×10-3m2/s

(10)

式中，E為空間電場(chǎng)強(qiáng)度矢量，V/m；δ為相對(duì)空氣密度，取值為1；μp、μn分別為正離子和負(fù)離子遷移率，μp=1.4×10-4m2/(V·s)，μn=1.8×10-4m2/(V·s)。

1.2計(jì)算模型

圖1為負(fù)極性電暈放電模型，其中，導(dǎo)線半徑為0.05 cm，同軸電極外半徑為20 cm，導(dǎo)線施加電壓為-30 kV，外電極接地。考慮到整個(gè)區(qū)域計(jì)算量大，而導(dǎo)線電離層很薄的情況[21]，本文對(duì)計(jì)算區(qū)域進(jìn)行了簡(jiǎn)化，在導(dǎo)線周圍設(shè)置了圓形人工邊界，邊界半徑為子導(dǎo)線半徑的20倍，即1 cm。人工邊界的電位通過(guò)靜電場(chǎng)計(jì)算設(shè)置為-15 kV。

圖1　負(fù)極性電暈放電模型Fig.1　Schematic diagram of negative corona discharge

電子二次發(fā)射是維持負(fù)極性電暈放電進(jìn)一步發(fā)展的重要機(jī)理，而目前對(duì)于電子二次發(fā)射過(guò)程有不同的描述機(jī)理，包括陰極表面場(chǎng)強(qiáng)較大導(dǎo)致的場(chǎng)致發(fā)射、放電等離子體中光子輻射產(chǎn)生的電子發(fā)射以及正離子碰撞陰極表面產(chǎn)生的電子發(fā)射[17]。文獻(xiàn)[17，18]均認(rèn)為正離子碰撞陰極是負(fù)極性電暈放電持續(xù)向前發(fā)展的關(guān)鍵原因。因此，本文在仿真時(shí)僅考慮正離子碰撞陰極產(chǎn)生二次電子這一過(guò)程。

電暈電流計(jì)算采用文獻(xiàn)[17]提出的方法，即電暈電流為導(dǎo)線表面的帶電粒子通量與導(dǎo)線表面感應(yīng)電荷變化量的環(huán)路積分。

(11)

式中，Is為電暈電流，A；n為陰極表面的法向量；dS為陰極表面的面積微元；σca為陰極表面的電荷密度，當(dāng)陰極為金屬電極時(shí)，其計(jì)算公式為

σca=-n(ε0εrEcathode)

(12)

式中，Ecathode為陰極表面的電場(chǎng)強(qiáng)度矢量，V/m。

1.3仿真模型及邊界條件

本文進(jìn)行仿真計(jì)算時(shí)，具體的網(wǎng)格剖分結(jié)構(gòu)圖如圖2和圖3所示。在導(dǎo)線表面附近區(qū)域及放電通道可能發(fā)展到達(dá)區(qū)域中，電子密度和電場(chǎng)強(qiáng)度等物理量相對(duì)于其他區(qū)域具有較大的數(shù)值及梯度，因此該部分區(qū)域的網(wǎng)格劃分非常細(xì)。在計(jì)算過(guò)程中，模型總共劃分為161 080個(gè)網(wǎng)格，平均網(wǎng)格質(zhì)量為0.982 8。

圖2　計(jì)算幾何結(jié)構(gòu)網(wǎng)格剖分圖(總體)Fig.2　Meshed simulation geometry(overall)

圖3　計(jì)算幾何結(jié)構(gòu)網(wǎng)格剖分圖(細(xì)化)Fig.3　Meshed simulation geometry(enlarged)

同時(shí)，對(duì)于流體模型來(lái)說(shuō)，邊界條件的設(shè)置至關(guān)重要，本文仿真時(shí)，具體的邊界條件設(shè)置如表1所示。其中，γ為電子的二次發(fā)射系數(shù)，取值為0.01[18]；De為電子擴(kuò)散系數(shù)，取值為1.28×10-3m2/s;Dp和Dn分別為正離子和負(fù)離子擴(kuò)散系數(shù)，取值為0[18]。

表1　邊界條件

2仿真結(jié)果及討論

本文主要基于Comsol Multiphysics有限元分析軟件的流體模塊進(jìn)行建模計(jì)算，計(jì)算的環(huán)境條件為300 K和1個(gè)標(biāo)準(zhǔn)大氣壓。在下列各仿真圖中，電流以及電場(chǎng)幅值均為絕對(duì)值。

2.1特里切爾脈沖仿真波形

圖4為一組特里切爾電流脈沖，在500 μs的時(shí)間內(nèi)呈現(xiàn)出80個(gè)連續(xù)脈沖。圖5為脈沖族的局部放大圖。圖6為單個(gè)電暈電流脈沖。

由圖4～圖6可見(jiàn)，仿真的電流脈沖符合特里切爾脈沖的基本特性，除第一脈沖外，其他脈沖的幅值和頻率均比較規(guī)則，脈沖峰值在1.5～2.5 mA，與文獻(xiàn)[21]中測(cè)量的脈沖幅值基本一致。單個(gè)脈沖呈較快的上升時(shí)間和短的持續(xù)時(shí)間，在一段死區(qū)時(shí)間后，新的脈沖重新發(fā)生。

圖4　仿真的特里切爾脈沖Fig.4　The simulated Trichel pulses

圖5　特里切爾脈沖局部放大圖Fig.5　The partial enlarged Trichel pulses

圖6　第一個(gè)特里切爾脈沖Fig.6　The first Trichel pulse

通過(guò)圖4發(fā)現(xiàn)，第一脈沖的幅值較大，其后的死區(qū)時(shí)間也長(zhǎng)于后續(xù)脈沖中較規(guī)則的死區(qū)時(shí)間，這種現(xiàn)象在文獻(xiàn)[15，18]中均有描述。引起這一現(xiàn)象的主要原因是：仿真開(kāi)始時(shí)，由于沒(méi)有空間電荷的抑制作用，第一脈沖發(fā)生時(shí)的導(dǎo)線表面場(chǎng)強(qiáng)要大于后續(xù)脈沖，從而導(dǎo)致第一脈沖發(fā)生時(shí)電離區(qū)的電離程度比較強(qiáng)烈，產(chǎn)生的空間電子和正負(fù)離子密度較大，從而使得第一脈沖幅值較大。由于第一脈沖產(chǎn)生的離子較多，導(dǎo)致脈沖發(fā)生后，空間離子對(duì)電場(chǎng)的抑制作用比后續(xù)離子強(qiáng)烈，空間電場(chǎng)減小程度較大，離子的移動(dòng)速度緩慢，導(dǎo)致導(dǎo)線表面電場(chǎng)恢復(fù)至起暈場(chǎng)強(qiáng)的時(shí)間被拉長(zhǎng)，從而使得整個(gè)死區(qū)時(shí)間變長(zhǎng)。

2.2第一脈沖的空間電場(chǎng)分布特性

為了更好地理解特里切爾脈沖的發(fā)展過(guò)程，本文對(duì)第一脈沖不同發(fā)展階段的空間電場(chǎng)分布進(jìn)行了分析討論。如圖6所示，t1～t5所對(duì)應(yīng)的時(shí)刻分別為16 ns、24 ns、100 ns、300 ns和400 ns。不同時(shí)刻的空間電場(chǎng)分布如圖7所示。

圖7　第一脈沖不同時(shí)刻電場(chǎng)分布Fig.7　Electric field distribution at different times during the first Trichel pulse

可以看出，在0時(shí)刻，電場(chǎng)分布為不受空間電荷影響的拉普拉斯(Laplace)場(chǎng)。隨著負(fù)極性電暈脈沖的起始與發(fā)展，碰撞電離產(chǎn)生的大量正離子向陰極方向移動(dòng)，使得正離子與陰極之間的區(qū)域場(chǎng)強(qiáng)發(fā)生畸變。正離子在碰撞陰極之后發(fā)射出的二次電子加劇了陰極附近的電離程度，從而在該區(qū)域產(chǎn)生了大量的帶電粒子，形成了包含電子云和正負(fù)離子云的等離子體區(qū)域。更加靠近陰極的正離子云使得陰極與等離子體之間的電場(chǎng)強(qiáng)度急劇增大，這一現(xiàn)象可從t1、t2時(shí)刻的電場(chǎng)分布圖看出。由于受空間電荷場(chǎng)反方向影響，等離子體通道內(nèi)的區(qū)域場(chǎng)強(qiáng)很小。

此后，隨著等離子體頭部越來(lái)越靠近陰極，高場(chǎng)強(qiáng)的電離區(qū)域變得非常小，最后，沒(méi)有足夠的高場(chǎng)強(qiáng)區(qū)域進(jìn)行碰撞電離，脈沖停止發(fā)展，正離子被陰極逐漸吸收，陰極表面及附近的場(chǎng)強(qiáng)受到負(fù)離子和電子的抑制作用開(kāi)始下降。

2.3第一脈沖的空間粒子分布特性

電子在脈沖放電過(guò)程中起至關(guān)重要的媒介作用，由于其運(yùn)動(dòng)速度較快，在與中性分子發(fā)生碰撞時(shí)，承擔(dān)能量傳遞的功能。而空間分布的負(fù)離子對(duì)導(dǎo)線表面電場(chǎng)有抑制作用。因此，下面對(duì)兩種帶電粒子在第一脈沖發(fā)展不同階段的空間密度分布進(jìn)行了描述。

由圖8和圖9可以看出，在脈沖起始、等離子體頭部向陰極發(fā)展的過(guò)程中，形成了大量的電子和負(fù)離子，兩者的密度分別達(dá)到1019m-3和1018m-3數(shù)量級(jí)。隨后在等離子體停止向前發(fā)展，逐漸消散的過(guò)程中，電子和負(fù)離子密度開(kāi)始下降，同時(shí)由于電場(chǎng)力的作用，兩者向遠(yuǎn)離導(dǎo)線的方向遷移。

圖8　第一脈沖不同時(shí)刻空間電子密度分布Fig.8　Electron density distribution at different times during the first Trichel pulse

圖9　第一脈沖不同時(shí)刻空間負(fù)離子密度分布Fig.9　Negative ion density distribution at different times during the first Trichel pulse

2.4死區(qū)時(shí)間內(nèi)的帶電粒子特性

圖10和圖11分別為第一個(gè)死區(qū)時(shí)間內(nèi)不同時(shí)刻的負(fù)離子和空間電場(chǎng)分布。由圖10可以看出，隨著時(shí)間的增加，負(fù)離子逐漸向外擴(kuò)散，由于在擴(kuò)散過(guò)程中與空間分布正離子的復(fù)合作用，負(fù)離子密度逐漸降低。從2 μs發(fā)展到8 μs時(shí)，最大負(fù)離子密度由2.03×1017m-3下降到4.5×1016m-3。通過(guò)圖11可以看出，隨著負(fù)離子向外遷移和復(fù)合，其對(duì)陰極表面及附近區(qū)域的電場(chǎng)抑制作用開(kāi)始減弱，該區(qū)域的電場(chǎng)逐漸上升。

圖10　死區(qū)時(shí)間不同時(shí)刻空間負(fù)離子二維分布(1/m3)Fig.10　Negative ion density distribution at different times during the dead time(1/m3)

圖11　死區(qū)時(shí)間不同時(shí)刻空間電場(chǎng)分布Fig.11　 Electric field distribution at different times during the dead time

圖12為第二個(gè)脈沖發(fā)展前不同時(shí)刻的空間電場(chǎng)變化圖。可以看出，在t=0.1 μs時(shí)，第一脈沖逐漸停止發(fā)展，陰極表面的場(chǎng)強(qiáng)開(kāi)始下降。而后，隨著負(fù)離子的遷移和復(fù)合，該區(qū)域電場(chǎng)逐漸恢復(fù)。當(dāng)負(fù)離子距離導(dǎo)線表面越來(lái)越遠(yuǎn)時(shí)，場(chǎng)強(qiáng)的恢復(fù)速度也逐漸放緩。可以看出，在t=10 μs與t=20 μs時(shí)，陰極表面及附近區(qū)域的場(chǎng)強(qiáng)分布幾乎重合，表面場(chǎng)強(qiáng)僅從6.32×106V/m恢復(fù)至6.37×106V/m。這主要是由于隨著負(fù)離子向外擴(kuò)散，其所處空間電場(chǎng)強(qiáng)度越來(lái)越小，由式(8)可知，當(dāng)電場(chǎng)減小時(shí)，負(fù)離子運(yùn)動(dòng)速率降低，從而導(dǎo)致場(chǎng)強(qiáng)恢復(fù)速度變慢。這一過(guò)程持續(xù)到導(dǎo)線表面場(chǎng)強(qiáng)達(dá)到起暈場(chǎng)強(qiáng)，第二次脈沖放電開(kāi)始。

圖12　不同時(shí)刻電場(chǎng)分布Fig.12　Electric field distribution at different times

2.5新的脈沖起始及發(fā)展

圖13～圖15分別為第二脈沖起始、發(fā)展及終止的不同時(shí)刻空間電場(chǎng)分布和電子、負(fù)離子密度分布。由圖13可以看出，當(dāng)新的脈沖起始時(shí)，陰極表面場(chǎng)強(qiáng)并未恢復(fù)至初始外加電壓條件下的Laplace場(chǎng)。在36.9 ～37.1 μs時(shí)，二次脈沖發(fā)生，碰撞電離形成的等離子體向陰極發(fā)展，陰極附近區(qū)域的電子密度和負(fù)離子密度急劇增大，該區(qū)域的場(chǎng)強(qiáng)也發(fā)生較大的畸變。

圖13　第二個(gè)脈沖發(fā)展不同時(shí)刻空間電場(chǎng)分布Fig.13　Electric field distribution at different times during the second pulse

圖14　第二個(gè)脈沖發(fā)展不同時(shí)刻空間電子密度分布Fig.14　Electron density distribution at different times during the second pulse

圖15　第二個(gè)脈沖發(fā)展不同時(shí)刻空間負(fù)離子密度分布Fig.15　Negative ion density distribution at different times during the second pulse

隨著二次脈沖結(jié)束，空間電子逐漸消散，負(fù)離子在電場(chǎng)力的作用下開(kāi)始向外遷移，導(dǎo)線表面電場(chǎng)也逐漸恢復(fù)。一段死區(qū)時(shí)間后，新的脈沖開(kāi)始重復(fù)發(fā)生。

3結(jié)論

本文建立了基于有限元的流體動(dòng)力學(xué)模型，對(duì)同軸電極結(jié)構(gòu)的負(fù)極性特里切爾脈沖進(jìn)行了數(shù)值仿真。通過(guò)計(jì)算，在500 μs的時(shí)間長(zhǎng)度內(nèi)，總共仿真出80個(gè)連續(xù)規(guī)則的特里切爾電流脈沖，除第一個(gè)脈沖外，脈沖幅值均在1.5～2.5 mA之間。

通過(guò)描述第一個(gè)脈沖不同發(fā)展階段的電場(chǎng)、電子和負(fù)離子密度變化情況，可以得出，當(dāng)碰撞電離產(chǎn)生的等離子體向陰極發(fā)展時(shí)，等離子體頭部電離程度非常強(qiáng)烈，場(chǎng)強(qiáng)、電子和負(fù)離子密度較大，其中，電子密度比負(fù)離子密度高出一個(gè)數(shù)量級(jí)。第一脈沖發(fā)展結(jié)束后，負(fù)離子逐漸向外擴(kuò)散，導(dǎo)線表面的場(chǎng)強(qiáng)開(kāi)始恢復(fù)。該過(guò)程一直持續(xù)到導(dǎo)線表面場(chǎng)強(qiáng)達(dá)到起暈場(chǎng)強(qiáng)，新的電暈脈沖發(fā)生。通過(guò)觀察新脈沖發(fā)展前后導(dǎo)線表面場(chǎng)強(qiáng)的變化可知，后續(xù)脈沖在起始時(shí)，導(dǎo)線表面場(chǎng)強(qiáng)不能恢復(fù)至無(wú)空間電荷影響的初始Laplace場(chǎng)強(qiáng)。同時(shí)，后續(xù)脈沖發(fā)生時(shí)空間電子密度和負(fù)離子密度要小于第一脈沖。

參考文獻(xiàn)

[1]劉振亞.特高壓電網(wǎng)[M].北京：中國(guó)經(jīng)濟(jì)出版社，2006.

[2]Huang Daochun，Shu Yinbiao，Ruan Jiangjun，et al.Ultra-high-voltage transmission in China：developments，current status，and future prospects[J].Proceeding of the IEEE，2009，97(3)：553-554.

[3]Maruvada P S.Corona performance of high-voltage transmission lines[M].Baldock：Research Studies Press Ltd.，2000.

[4]劉陽(yáng)，崔翔，盧鐵兵，等.濕度對(duì)絞線正極性電暈電流脈沖及其無(wú)線電干擾影響的試驗(yàn)[J].電工技術(shù)學(xué)報(bào)，2015，30(3)：79-89.

Liu Yang，Cui Xiang，Lu Tiebing，et al.Experimental investigation of humidity impact on positive corona current pulse and its radio interference of stranded conductors[J].Transactions of China Electrotechnical Society，2015，30(3)：79-89.

[5]喬驥，鄒軍，袁健生，等.采用有限差分求解高壓直流輸電線路空間離子流場(chǎng)的新方法[J].電工技術(shù)學(xué)報(bào)，2015，30(6)：85-91.

Qiao Ji，Zou Jun，Yuan Jiansheng，et al.A new finite difference based approach for calculating ion flow field of HVDC transmission lnes[J].Transactions of China Electrotechnical Society，2015，30(6)：85-91.

[6]胡琴，吳執(zhí)，舒立春，等.交流電場(chǎng)下水滴對(duì)導(dǎo)線電暈特性的影響[J].電工技術(shù)學(xué)報(bào)，2015，30(18)：237-245.

Hu Qin，Wu Zhi，Shu Lichun，et al.The crona caracteristics of the cnductor atached with wter dops under AC electrical feld[J].Transactions of China Electrotechnical Society，2015，30(18)：237-245.

[7]謝莉，趙錄興，陸家榆，等.有限長(zhǎng)高壓直流輸電線路無(wú)線電干擾電磁場(chǎng)的計(jì)算方法研究[J].電工技術(shù)學(xué)報(bào)，2016，31(1)：96-102.

Xie Li，Zhao Luxing，Lu Jiayu，et al.Study on the clculation mthod for rdio interference electromagnetic field of HVDC transmission line with finite length[J].Transactions of China Electrotechnical Society，2016，31(1)：96-102.

[8]劉云鵬，朱雷，律方成，等.特高壓電暈籠直流分裂導(dǎo)線正極性電暈起始特性分析[J].電工技術(shù)學(xué)報(bào)，2013，28(1)：73-79.

Liu Yunpeng，Zhu Lei，Lü Fangcheng，et al.Analysis of the positive corona onset characteristic of the bundle conductors in the UHV corona cage[J].Transactions of China Electrotechnical Society，2013，28(1)：73-79.

[9]何旺齡，萬(wàn)保權(quán)，劉震寰，等.特高壓交流試驗(yàn)線段的無(wú)線電干擾特性研究及長(zhǎng)線路預(yù)測(cè)[J].高電壓技術(shù)，2014，40(12)：3953-3960.

He Wangling，Wan Baoquan，Liu Zhenhuan，et al.Radio interference characteristics of UHVAC test lines and prediction for long lines[J].High Voltage Engineering，2014，40(12)：3953-3960.

[10]Lin Lie，Zhang Peng，Wang Yongqing.Analysis of instabilities in non-equilibrium plasmas[J].Chinese Physics Letters，2004，21(10)：1993-1996.

[11]Trichel G W.The mechanism of the negative point to plane corona near onset[J].Physics Review，1938，54(12)：1078-1084.

[12]Loeb L B.The mechanism of the trichel pulses of short time duration in air[J].Physical Review，1952，86(2)：256-257.

[13]Lama W L，Gallo C F.Systematic study of the electrical characteristics of the Trichel current pulses from negative needle-to-plane coronas[J].Journal Applied Physics，1974，45(1)：103-13.

[14]歐陽(yáng)吉庭，張子亮，張宇，等.空氣負(fù)電暈Trichel 脈沖特性的實(shí)驗(yàn)研究[J].高電壓技術(shù)，2014，40(4)：1194-1200.

Ouyang Jiting，Zhang Ziliang，Zhang Yu，et al.Experimental study on the characteristics of negative-corona Trichel pulses in air[J].High Voltage Engineering，2014，40(4)：1194-1200.

[15]Sattari P，Castle G S P，Adamiak K，et al.Numerical simulation of Trichel pulses in a negative corona discharge in air[J].IEEE Transactions on Industry Applications，2011，47(4)：1935-1943.

[16]Sattari P，Castle G S P，Adamiak K，et al.Trichel pulse characteristics—negative corona discharge in air[J].Journal of Physics D：Applied Physics，2011，44(15)：155502-155512.

[17]Tran T N，Golosnoy I O，Lewin P L，et al.Numerical modelling of negative discharges in air with experimental validation[J].Journal of Physics D：Applied Physics，2011，44(1)：015203.

[18]Yin Han，Zhang Bo，He Jinliang，et al.Modeling of Trichel pulses in the negative corona on a line-to-plane geometry[J].IEEE Transactions on Magnetics，2014，50(2)：473-476.

[19]伍飛飛，廖瑞金，楊麗君，等.棒-板電極直流負(fù)電暈放電特里切爾脈沖的微觀過(guò)程分析[J].物理學(xué)報(bào)，2013(11)：348-357.Wu Feifei，Liao Ruijin，Yang Lijun，et al.Analysis on micro process of DC negative Trichel pulses in a bar-plane electrode[J].Acta Physica Sinica，2013(11)：348-357.

[20]廖瑞金，伍飛飛，劉康淋，等.棒-板電極直流負(fù)電暈放電脈沖過(guò)程中的電子特性研究[J].電工技術(shù)學(xué)報(bào)，2015，30(10)：319-329.

Liao Ruijin，Wu Feifei，Liu Kanglin，et al.Simulation of characteristics of electrons during a pulse cycle in Bar-Plate DC negative corona discharge[J].Transactions of China Electrotechnical Society，2015，30(10)： 319-329.

[21]尹晗.直流輸電線路高頻電暈電流與無(wú)線電干擾的轉(zhuǎn)換關(guān)系[D].北京：清華大學(xué)，2014.

Characteristics of Negative Corona Trichel Pulses in a Coaxial Electrode System

He Wangling1He Junjia1Zhang Jin2Wan Baoquan2

(1.School of Electrical and Electronic EngineeringHuazhong University of Science and Technology Wuhan430074China 2.China Electric Power Research Institute(Wuhan)Wuhan430074China)

AbstractIn this paper,the characteristics of negative corona discharge in air at atmospheric pressure are studied by using the hydrodynamic drift-diffusion model with a coaxial electrode system,and a set of regular Trichel pulses is obtained.It is found that the pulse amplitude is between 1.5 mA and 2.5 mA except the first one.Then the distributions of the electric field,the ion densities at different stages during the first pulse are plotted.The distribution of the negative ion density and its influence on the conductor surface electric field are discussed during the dead time between the first pulse and the second pulse.After that,the distributions of the electric field,the electron and negative ion densities around the second pulse are plotted.The simulation results show that the electric field at the conductor surface increases when the negative ions remove from cathode under the influence of electric field force;and the recovery velocity of electric field slows down when the negative ions move away.Besides,before the electric field at the conductor surface recovers to its original Laplacian field strength,the next new pulse begins to happen.It is also found that the maximum values of the electric field,the electron and negative ion densities are reduced during the second pulse comparing with those in the first pulse.

Keywords：Negative corona，Trichel pulse，electric field distribution，ion density

收稿日期2015-04-15改稿日期2015-07-05

作者簡(jiǎn)介E-mail：wanglinghe88@gmail.com(通信作者) E-mail：h-junjia@vip.sina.com

中圖分類號(hào)：TM723

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展(973)計(jì)劃資助項(xiàng)目(2011CB209402)。

何旺齡男，1988年生，博士研究生，研究方向?yàn)榻涣鬏旊娋€路的電暈放電與電磁環(huán)境。

何俊佳男，1968年生，教授，博士生導(dǎo)師，研究方向?yàn)殡娀‰娊佑|、高電壓與絕緣技術(shù)以及脈沖功率技術(shù)。