一種新型光學電流傳感器結構設計與仿真研究

萬代，鐘力生，齊飛，周恒逸，趙邈，段緒金

(1.國網湖南省電力公司電力科學研究院，湖南 長沙 410007；2.西安交通大學，陜西 西安 710049)

一種新型光學電流傳感器結構設計與仿真研究

萬代1，鐘力生2，齊飛1，周恒逸1，趙邈1，段緒金1

(1.國網湖南省電力公司電力科學研究院，湖南 長沙 410007；2.西安交通大學，陜西 西安 710049)

本文設計對稱光程反射式光學電流傳感器實驗系統，使用COMSOL多模場耦合仿真軟件對所設計的光學傳感系統的傳輸場進行模擬仿真，得到對稱光程反射式光學電流傳感器對法拉第磁光效應的響應度，其理論靈敏度高于傳統結構的靈敏度580倍。未來將基于所建立的實驗系統對傳感器的潛在應用做進一步的深入研究。

光學電流傳感器；Faraday效應；結構設計；仿真，靈敏度

電流傳感器作為電網中至關重要的環節，關系著電力系統的安全運行。近百年來，傳統電磁式電流互感器由于其成熟的測量技術、簡單的結構和不易損壞等優點，已獲得廣泛應用；隨著電網規模的日益增大和運行機制的不斷發展，由于其結構和使用條件的特殊性，它在電力系統電壓等級不斷提高的情況下普遍暴露出安全性低〔1－2〕、電磁干擾嚴重〔3〕、環境友好性差、存在磁飽和現象〔3〕影響測量精度〔4－5〕、成本高以及裝配難度大等突出問題。

傳統電流傳感器的缺陷促使人們開始探索更先進的電流傳感器。20世紀60年代出現的半導體集成電路技術、激光技術以及70年代初出現的光纖通信技術，使光學電流傳感技術迅速興起與發展。光學電流傳感器理論上具有傳統電磁式電流傳感器功能并能克服其上述缺點〔6－7〕。但是光學材料中存在的線性雙折射現象及其對電流磁場較低的敏感度一直成為制約光學電流傳感器的主要原因。

目前，對于全光纖電流傳感技術的研究主要集中于大電流測量領域，量程一般在102～106A，且普遍存在欠缺穩定性的問題〔8－10〕。為促進智能電網發展，符合智能電網提出的先進傳感技術要求，以及拓寬光學電流傳感器的應用領域，文中設計了對稱光程反射式光學電流傳感器測試系統，理論上可解決光學電流傳感器中的線性雙折射現象，對提高全光纖小電流傳感器測量靈敏度，拓展其應用領域具有指導意義。

1 理論

1.1 Faraday效應

光學電流傳感器的工作原理是基于法拉第效應，法拉第效應主要是由于外磁場對光學材料磁化后，使材料本身具有磁矩，此時在光學材料內部傳輸的偏振光的電磁場將與材料的磁矩相互作用，從而影響到光波的電場分布，體現在宏觀上即為光波的偏振態發生旋轉〔11〕。此偏振態的旋轉角度稱為法拉第偏轉角。通過準確測量法拉第偏轉角反映了光波所處的外磁場大小，于是可進一步計算出產生此磁場的電流大小。

光學材料的介電常數ε會在材料具有磁矩后發生變化，而材料所具有的磁矩與材料磁化強度有關，于是介電常數的變化可以用材料的磁化強度M來描述。將介電常數張量ε的變化用M的冪級數展開，根據張量的性質并應用昂薩格關系〔11〕εij(M)＝εij(－M)， 介電常數張量的各個分量可以表示為:

對于光學電流傳感器所選用的光學晶體，其對稱性高于四方晶系，且設定材料分布均勻，于是有ε11＝ε22， ε12＝－ε21， ε13＝ε31＝ε23＝ε32＝0。 式 (1)可變為:

式中 設定z軸為光軸方向，xy平面為光學晶體橫截面，分別沿x，y，z坐標軸方向的介電常數ε11＝ε22＝ε33＝ε； δ為磁極化系數，與外磁場相關。

偏振光在光纖中傳輸時滿足麥克斯韋方程組。

式中 B為磁感應強度； D為電位移矢量。入射線偏振光的表達式為:

式中 E0為光波電場矢量的振幅；s為光波矢方向的單位矢量；H0為光波磁場矢量的振幅。

真空中光速c的表達式可以寫為c2＝ (μ0ε0)－1。于是聯立式 (2)和式 (3)可得:

將電場強度和介電常數分別用向量和矩陣表示，代入式 (5)；且設定光波沿光軸方向傳輸，有s(0，0，1):

對式 (6)進行矩陣運算可得:

將式 (7)帶入式 (6)第1，2行系列方程可得:

式 (8)代表兩旋向相反且以不同的速度c/n＋和c/n－沿光軸方向傳播的圓偏振光，經過一段距離后，兩束光之間存在相位差，此時兩束圓偏振光合并后依舊是線偏振光，并且相對于初始時的線偏振光的偏振態有一個偏轉角。

設線偏振光沿z方向傳播l距離，此時偏振態偏轉角度為θ。若線偏振光的初始偏振態是沿x軸方向，所以tanθ＝Ey/Ex。

聯合式 (7)可得:

式中 V為光纖材料的費爾德常數/rad·A－1；H為測試電流產生的磁場。

1.2 線性雙折射效應

若光學材料的介電常數非等值分布，那么線偏振光在其中傳輸時會產生奇異現象，稱之為線性雙折射效應。設光學晶體的橫截面上 εx≠εy。解式(6)，將得到:

相應的，方程組式 (6)的解將變為:

式 (12)代表一束橢圓偏振光。由此可知，由于線性雙折射的存在，將導致入射線偏振光退偏為橢圓偏振光。

2 結構設計

設計了一種新型結構的高靈敏度光學電流傳感器實驗系統，如圖1所示，可命名為對稱光程反射式光學電流傳感器。此系統由光源、格蘭棱鏡、光纖耦合器、傳感光纖、自聚焦透鏡、對稱光程反射式傳感頭、沃拉斯頓棱鏡、光探測器組成。

在本系統中激光器發出的光由格蘭棱鏡變為線偏振光，然后在傳感頭部位發生法拉第磁光效應使光波偏振態產生變化，然后由檢偏器分解Ex和Ey分別被光探測器接收后轉換為電流信號，最后通過計算軟件對兩路信號進行運算，最終得到光波偏振態夾角準確值。

圖1 對稱螺旋嵌套型全光纖電流傳感器

以往傳統的光學傳感頭中線性雙折射非常大以及對磁場靈敏度偏低，是導致光學電流傳感器穩定性和可靠性欠缺主要原因。所以設計了一種新型的對稱光程反射式傳感頭，這種結構與傳統直線型結構相比，可提高靈敏度達800倍，從而可提高全光纖電流傳感器測量精度與穩定性。

光程反射式光學電流傳感頭由通電螺線管和管內的圓柱形石英晶體組成，螺線管內徑為30 mm、外徑為40 mm、長度為230 mm，環形圓柱中心為空洞，每個螺線管上環形導線匝數為 n1＝1 500。在螺線管內部產生集中的磁場 (與單位長度的導線匝數有關)，然后將兩面鍍反射膜的石英玻璃置于螺線管軸線上。

圖2(a)， (b)分別為光程反射式光學電流傳感頭中柱狀石英晶體的徑向截面和軸向截面的剖面圖。柱狀石英晶體的橫截面是直徑為30 mm的圓，石英晶體圓柱的軸向長度為 500 mm。如圖2(a)所示，選取100根直徑為0.1 mm的高折射率純石英晶體圓柱體 (光纖預制棒)，安置在低折射率石英玻璃圓柱中作為光通道，每個相鄰光通道之間的最短距離為2 mm，最長距離為2.82 mm。其制備過程類似于預制棒法拉制實芯型光子晶體光纖的制備過程，首先選取一根直徑為500 mm的低折射率玻璃圓柱體，使用精密鉆孔機床在玻璃棒中鉆100個直徑為1 mm的圓柱形空氣孔。從玻璃棒的徑向截面觀測，這100個空氣孔組成10×10的正方形陣列，每相鄰最近的2個空氣孔之間的距離為10 mm。從玻璃棒的軸向截面觀測，這100個圓柱形空氣孔的對稱軸與玻璃棒的對稱軸平行。然后選取100根直徑為1 mm的高折射率純石英光纖預制棒，將這些高折射率石英棒分別安置于玻璃圓柱體內的100個空氣孔中，保持固定。將玻璃棒固定在拉絲機的端部，逐漸加熱至玻璃軟化，保持加熱區均勻加熱并對玻璃棒進行勻速拉制，合理選取拉制速度將玻璃棒拉制為直徑為30 mm的圓柱體，此時拉制出的細玻璃棒結構與原結構一致，只是徑向尺寸縮小10倍。選取拉制好的玻璃棒中間最均勻的部分，截取長度為500 mm、直徑為30 mm的圓柱體，此即為多光路反射式光纖小電流傳感頭中的柱狀石英晶體。

圖2 對稱光程反射式光學電流傳感頭

由于光程反射式光學電流傳感頭的柱狀石英晶體中光通道為高折射率純石英棒，包層為低折射率玻璃，類似于光線在光纖中傳輸過程，線偏振光將在高折射率純石英光通道中傳輸。類似于折射率引導型自聚焦透鏡原理，將玻璃棒的兩端加工為折射率引導型自聚焦透鏡陣列并在兩端面處鍍上法拉第鏡式保偏反射膜，以保證線偏振光由光通道傳輸至端面處反射入相鄰的光通道中，并消除了反射相移對實驗結果的影響。根據法拉第旋光效應的非互易性，通過兩面法拉第鏡式保偏反射膜多次反射后法拉第轉角將成倍增大。

利用安培環路定理對式 (10)中外磁場與測試電流的關系進行解析，式 (10)可簡化為:

式中 n1為通電螺線管匝數；n2為傳感頭中光通道的數目。

3 仿真研究

本文使用COMSOL多模場耦合仿真軟件對系統的傳輸場進行仿真模擬。COMSOL多模場耦合仿真軟件是以有限元法為基礎，通過求解偏微分方程(單場)或偏微分方程組 (多場)來實現真實物理現象的仿真。本實驗中應用COMSOL把含時間變量的Maxwell方程在空間中轉化為差分方程。在這種差分格式中每個網格點上的電場 (或磁場)分量僅與它相鄰的磁場或 (電場)分量及上一時間步該點的場值有關。在每一時間步計算網格空間各點的電場和磁場分量，隨著時間步的推進，即能直接模擬電磁波的傳播及其與物體的相互作用過程。由于在差分格式中被模擬空間電磁性質的參量是按空間網格給出的，因此只需對相應空間點設定適當的參數，對介質的非均勻性、各向異性、色散特性和非線性等結構均能進行精確模擬〔12〕。

3.1 傳統直線型結構仿真研究

根據傳統直線型結構的光學電流傳感器應用實例、設計仿真模型，對其進行模擬。設置電纜為銅圓柱體，半徑1 cm，高度1 m，通電流1 A，在xyz坐標系中電纜軸線沿z軸，徑向截面處于xy平面，于是電流密度設為jx＝jy＝0，jz＝3 183 A/m2。按照這些參數繪制結合結構并設定模型的材料性質與電學性質；周圍半徑0.2 m、高1.4 m的圓柱體范圍內設為空氣，最外層設為PML無限元，作為計算區域。圖3所示為xy平面內y方向磁場Hy分布的仿真計算結果，圖3中右方強度軸顏色的深淺表示磁場模的大小，正負號表示磁場的方向，從圖3中可以看出電纜徑向截面的兩側磁場方向相反，且越遠離電纜磁場強度越小。由于軟件中選擇顯示場解時一次只能選擇一個參數，所以圖中未顯示x方向的磁場強度，計算結果表明Hx的分布與Hy一致。

圖3 直線型傳感頭磁場模量仿真結果

光學電流傳感器在實際測量中，光纖環以螺旋線形式纏繞在電纜上，不可能完全閉合，所以不能直接應用安培環路定理，需要對實際情形進行仿真分析。設光纖路徑為一螺旋線，其螺旋半徑為2 cm，螺距為0.125 mm，共250匝，螺旋路徑起點在電纜中點處，縱向長度為31.25 mm。輸入參數化曲線數據，將磁場仿真結果導出，并計算磁場沿光路的的矢量積分，結果為249.875 A(這里A為Hdl的單位，并非描述電流)。仿真結果與理論計算值250 A非常接近，中間有0.125 A的區別主要是由于光纖環分布在z軸原點兩側并非集中于z軸原點截面上，而螺旋路徑無法完全閉合造成。

3.2 對稱光程反射式結構仿真研究

按照對稱光程反射式光學電流傳感器模型的設計方案，多光路反射式結構的仿真模型如圖4所示。建立模型時在螺線管中嵌入均勻分布的10×10線陣列，每個光通道之間最短間距為2 mm，傳感頭長度為 3 cm，電流密度設為 jx＝－652 174 y/(x2＋y2)0.5A/m2，jy＝652 174x/(x2＋y2)0.5A/m2，jz＝0。通過計算其仿真模型，磁場沿光路的矢量積分結果為1.449×105A/m，在波長為650 nm的紅光激光器激勵下，法拉第磁光偏轉角為θ＝0.857 8 rad。根據計算結果可知，對稱光程反射式結構與直線型結構相比，靈敏度提高了580倍。

圖4 對稱光程反射式結構磁場分布仿真

在多光路反射式傳感頭的設計中，設定了傳感頭長度為3 cm，在實際制作中可能會制作出不同尺寸的傳感頭，由于傳感頭長度的不同，在一定程度上也會影響法拉第磁光偏轉角的測量。設定線路中通過1 A仿真電流，如圖5(a)為不同長度的傳感頭對測量結果的影響，由計算結果可以發現，當傳感頭軸向長度小于螺線管長度時，系統靈敏度會大幅度降低；當傳感頭長度達到50 cm以上時，曲線漸漸趨于平緩，此時系統的靈敏度趨于穩定，測量結果誤差保持在0.1%以內，由此可知傳感頭的長度至少應該設置在50 cm以上。

若傳感頭尺寸達到要求，在實驗中由于操作失誤等原因有可能使將傳感頭位置安裝錯誤，導致傳感頭中點與螺線管軸線中點產生偏差，設定線路中通過1 A仿真電流，圖5(b)為傳感頭在軸線位置上偏差對測量結果的影響，在圖中可以發現當傳感頭的位置安放錯誤時會對測量結果帶來非常大的影響，隨著位置偏差的增大，測量結果誤差將越來越大，當位置偏差保持在1 cm以內，測量結果誤差低于0.2%。

圖5 傳感區域長度與位置偏差對測量結果的影響

4 結論

該對稱光程反射式光學電流傳感器測試系統采用COMSOL Multiphysic多模場耦合模擬仿真軟件對所設計的實驗系統進行仿真計算。通過與傳統直線型光學電流傳感器對比發現，此傳感器可解決小電流準確測量的難題，具體結論:

1)對稱光程反射式光學電流傳感器靈敏度優于靈敏度傳統直線型光學電流傳感器580倍，擴展了全光纖電流傳感器的測量范圍，解決了將全光纖電流傳感器應用于小電流測量領域的技術難題。

2)對稱光程反射式光學電流傳感器由對偶式結構減小了線性雙折射對測量結果的影響。

3)對稱光程反射式光學電流傳感器結構加工精度要求較高。傳感頭的長度至少應該設置在50 cm以上；軸向位置偏差保持在1 cm以內，測量結果誤差低于0.2%。

〔1〕劉曄，王采堂，蘇彥民.電力系統適用光學電流互感器研究新進展 〔J〕.電力系統自動化，2000，17:60-64.

〔2〕 Emerging Technologies Working Group， Fiber Optic Sensors Working Group.Optical Current Transducers for Power Systems: AReview〔J〕.IEEE Transaction on Power Delivery，1994，9 (4):1 778-1 787.

〔3〕Wan Dai，Zhong Lisheng，George Chen.A New Structure of Full OpticalFiber Current Sensor〔C/CD〕. IEEE The 10th International Conference on Properties and Applications of Dielectric Materials，Bangalore，India，2012.

〔4〕Fabio Vieira Batista de Nazare，Mazare Martins Werneck.Hybrid Optoelectronic Sensor for Current and Temperature Monitoring in Overhead Transmission Lines〔J〕.IEEE Sensors Journal，2012，12(5):1 193-1 194.

〔5〕Lee Byoungho.Review of the present status of Optical Fiber Sensors〔J〕.Optical Fiber Technology，2003，9:57-79.

〔6〕Apurba Ghosh，Punya Brata Dutta Gupta，Ajit Kumar Mandal. Development of a Fiber-Optic Current Sensor with Range-Changing Facility Using Shunt Configuration〔J〕.IEEE Sensors Journal，2013，13(4):1 347-1 354.

〔7〕Li Bingjun，Li Lijing.An Overview of the Optical Current Sensor〔C〕.2012 International Conference on Computer Science and Electronics Engineering.2012，3:202-206.

〔8〕秦冉，王倩倩，楊世海.全光纖電流互感器現場運行誤差特性研究 〔J〕.江蘇電機工程，2013，32(02):7-9，14.

〔9〕陳銘明，盧樹峰，包玉樹.光學電流互感器實時誤差分析系統的設計 〔J〕.Jiangsu electrical engineering，2013，32(02):47-49，53.

〔10〕陳文中，林一，周健.數字化變電站全光纖電流互感器準確度校驗 〔J〕.華東電力，2009，37(12):2 022-2 024.

〔11〕劉公強，樂志強，沈德芳.磁光學 〔M〕.上海:上海科技出版社，2001.

〔12〕 COMSOL Corporation，幾 何建 模 用 戶 指 南， COMSOL Multiphysics操作手冊叢書，2009.

Structure design and simulation research on a new type of optical current sensor

WAN Dai1，ZHONG Lisheng2，QI Fei1，ZHOU Hengyi1，ZHAO Miao1，DUAN Xujin1

(1.State Grid Hunan Electric Power Corporation Research Institute，Changsha 410007，China；2.Xi’an Jiaotong University，Xi’an 710049，China)

A new structure of optical current sensor is designed in this paper.The COMSOL software is used to simulate the transmission field of optical sensing system.The responsiveness of symmetry optical path reflection type optical current sensor to Faraday effect can be obtained by the simulation study.The results show that the sensitivity of current sensor of symmetry optical path reflection type optical is 580 times higher than the sensitivity of traditional linear type optical current sensor.Finally the influence rule of measuring result error caused by sensing area positional deviation is analyzed.Further study about potential application of optical current sensor will proceed based on the optical current sensing system of symmetry optical path reflection type.

optical current sensor；Faraday effect；structure design；simulation；sensitivity

10.3969/j.issn.1008-0198.2015.04.003

TM452

B

1008-0198(2015)04-0009-05

萬代(1985)，男，博士，工程師，研究方向為光纖電流傳感技術、電介質與電氣絕緣技術、配電設備技術研究。

2015-06-16

鐘力生 (1962)，男，博士，教授 (博導)，研究方向為電介質與電氣絕緣技術、聚合物光電材料及器件、生物電介質及應用。

齊飛(1985)，男，工程師，研究方向為配電設備技術研究。