吳 文 軍

(上海電力學(xué)院， 上海 200090)

0 引 言

同塔多回高壓輸電系統(tǒng)的諧波污染，除了線路交流電的高次諧波外，電磁干擾也是一大成因。因此，任意分布導(dǎo)體散射場的分析與預(yù)測(cè)，對(duì)于同塔多回高壓輸電系統(tǒng)的設(shè)計(jì)與布局、電磁環(huán)境的研究和諧波干擾的抑制具有重要意義。其中，在任意干擾源散射作用下，復(fù)雜導(dǎo)電體近場分布的數(shù)值計(jì)算方法的研究具有十分關(guān)鍵的作用和意義，因而得到了理論研究和工程技術(shù)人員的廣泛關(guān)注和深入探討[1-3]。

任意散射源照射下，復(fù)雜導(dǎo)電體近場分布的數(shù)模計(jì)算在系統(tǒng)電磁兼容性的設(shè)計(jì)和預(yù)測(cè)方法中常見的方法是將同塔多回高壓輸電系統(tǒng)中導(dǎo)電體模擬為面元模型，構(gòu)造模擬面近似導(dǎo)電體表面，在模擬面上建立子域基電流模型，并運(yùn)用積分方程結(jié)合矩量法求解模擬面上的電流和電荷分布[4-6]。面元模型不僅可描述任意形狀的導(dǎo)電體，且適用于計(jì)算近場和表面量，但計(jì)算更為復(fù)雜，所以尋求快速、合理的求解方法成為成功實(shí)施面元模型解決電磁兼容性預(yù)測(cè)的關(guān)鍵。

目前，解決同塔多回高壓輸電系統(tǒng)非線性電磁干擾求解問題的主要途徑有：① 頻域分析法[7]。求解頻域積分-微分方程，利用高斯-牛頓法反演變換得到系統(tǒng)電磁特性的時(shí)域解，系統(tǒng)能夠快速收斂，但穩(wěn)定性不高，易發(fā)散。② 直接時(shí)域法[8-9]。直接在時(shí)域求解時(shí)域積分、差分方程，利用梯度法可確保系統(tǒng)穩(wěn)定反演，獲得精確的散射模式的數(shù)值解，但收斂速度較慢，特別是逼近極小值時(shí)更慢。③ 最小二乘估計(jì)方法[10]。結(jié)合矩量法求解積分方程，求得矩陣形式的解，然后在復(fù)頻域利用圍線積分和柯西定理，把復(fù)頻域解表示成留數(shù)項(xiàng)之和，可直接把時(shí)域解表示成復(fù)指數(shù)項(xiàng)之和。該方法計(jì)算精確度高、穩(wěn)定性好，兼有前兩種方法的優(yōu)點(diǎn)，但計(jì)算量較大，影響了實(shí)際效果與應(yīng)用范圍。

頻域分析法是分析瞬變電磁場問題的經(jīng)典方法，但在揭示瞬態(tài)響應(yīng)機(jī)理方面不如直接時(shí)域法[11-12]。直接時(shí)域法能夠提供瞬態(tài)響應(yīng)隨時(shí)間推進(jìn)的演變過程，但不能像最小二乘估算方法那樣反映出瞬態(tài)響應(yīng)域目標(biāo)特征參量間的直接關(guān)系[13-14]。在多回高壓輸電系統(tǒng)復(fù)雜電磁環(huán)境中，大多數(shù)散射的積分方程分析是在頻域，微分方程法是在時(shí)域，而時(shí)域積分方程位于兩者之間，極值分解方法通過結(jié)合矩量法求解時(shí)域積分方程，使其在確保穩(wěn)定性的同時(shí)，通過預(yù)測(cè)模型結(jié)構(gòu)獲得精確解，因而論文基于極值分解法，結(jié)合最小二乘估算策略，提出模式綜合電磁預(yù)測(cè)方法，為電磁兼容性預(yù)測(cè)面元模型的建立提供了有效的求解途徑，實(shí)現(xiàn)了同塔多回高壓輸電線路簡單高效、切實(shí)可行的電磁兼容預(yù)測(cè)與分析[15]。

1 耦合電磁干擾問題描述

針對(duì)同塔多回高壓輸電系統(tǒng)中普遍存在的非線性耦合電磁干擾，考察N(N>1)維散射向量任意分布情形下的耦合電磁干擾問題，滿足：

(1) 對(duì)同塔多回高壓輸電系統(tǒng)中散射源的幾何結(jié)構(gòu)不作限制；

(2) 對(duì)向量元素的類型不作限制，無需歸一化處理；

(3) 對(duì)同塔多回高壓輸電系統(tǒng)中散射源的放置位置不作限制。

根據(jù)電氣系統(tǒng)的電磁特性，考察電磁散射模型

F(θ,φ)=θ0Fθ(θ,φ)+φ0Fφ(θ,φ)

(1)

式中：θ0、φ0為θ、φ空間球面坐標(biāo)系中沿角度方向的單位矢量。據(jù)式(1)，總散射能量為

方便起見，定義角球面坐標(biāo)對(duì)(θ,φ)為符號(hào)Λ，對(duì)應(yīng)的角度微分為dΛ，則總散射能量可表示為

(2)

(3)

式(3)為最小二乘最優(yōu)化問題的標(biāo)準(zhǔn)方程形式。定義離散空間中的N維列向量組成矢量p，激勵(lì)值am組成列向量a，則式(3)可寫為

p=Ga

(4)

2 耦合電磁干擾的互耦與極值分解

若式(4)的Gram矩陣中存在互耦，則會(huì)令G成為病態(tài)矩陣，表明擴(kuò)展矩陣中的某些向量存在相關(guān)性。Gram矩陣中互耦程度的測(cè)量依據(jù)它的行列式秩數(shù)，當(dāng)

det[(fm,fn)]?0

由Gram矩陣的直接轉(zhuǎn)置求解式(4)時(shí)，可能會(huì)因舍入誤差導(dǎo)致無解。極值分解技術(shù)為此類情況提供了有效的解決途徑。

若A為復(fù)雜M×N矩陣，則可構(gòu)造(M+N)×(M+N)矩陣:

式中：W為Hermit矩陣。由構(gòu)成向量qj做歸一化變換，使得

Wqj=σjqj,j=1,2,…,M+N

(5)

式中：σj為特征值。W為對(duì)角陣,將qj分解為

式中：uj和vj的分別為M維和N維。式(5)可表示為

Avj=σjuj

(6)

AHuj=σjvj

(7)

式(6)兩邊左乘AH，結(jié)合式(7)，可得

(8)

同樣，式(7)兩邊左乘A，結(jié)合式(6)，可得

(9)

(10)

利用式(10)可生成N×N維矩陣INN，使得:

(11)

(12)

式中：σj為A的極值。所有極值選為正，經(jīng)歸一化變換，可得式(12)的矩陣形式,

A=UHVH

式中：U、V分別為M×M維、N×N維矩陣；H為M×N維矩陣，滿足:

式(12)通常稱為降階極值分解形式，可將U、V降為R列，定義UR、VR使得:

(13)

利用式(13)，求解最小二乘最優(yōu)化問題的標(biāo)準(zhǔn)方程式(3)，可得:

(14)

由此可向量得解

求取極值實(shí)現(xiàn)了散射向量元素與激勵(lì)參數(shù)間的相位轉(zhuǎn)換，從而降低電磁散射效能。

利用模式綜合技術(shù)，將非線性電磁耦合等效為約束條件，采用最小二乘方法，求解符合應(yīng)用需求的最優(yōu)觀測(cè)值，能夠有效抑制能量耦合的影響，確保導(dǎo)電體優(yōu)良的電磁兼容性能，推導(dǎo)可得求解結(jié)果：

(15)

3 實(shí)驗(yàn)結(jié)果分析

對(duì)110 kV同塔四回輸電線路采用垂直排列方式，線路模型如圖1所示，線路主要參數(shù)見表1所示。

參數(shù)數(shù)量桿塔高度/m58導(dǎo)線型號(hào)ACSR?720/50導(dǎo)線直徑/mm36．2導(dǎo)線直流電阻/(Ω·km-1)0．03984避雷線型號(hào)LBGJ?150?40AC避雷線直徑/mm15．75避雷線直流電阻/(Ω·km-1)0．2935避雷線接地方式中性點(diǎn)直接接地

為了驗(yàn)證所設(shè)計(jì)的電磁抗擾系統(tǒng)的穩(wěn)定性，基于dSPACE系統(tǒng)， 實(shí)時(shí)仿真同塔多回高壓輸電系統(tǒng)電磁場環(huán)境，對(duì)同塔多回高壓輸電系統(tǒng)12條線路中近域耦合電磁干擾最為強(qiáng)烈的一條線路進(jìn)行電磁兼容性能分析，辨識(shí)獲得的散射向量幅值與相角分別如圖2、3中實(shí)線所示的幅值軌跡與相角軌跡。利用模式綜合預(yù)測(cè)算法對(duì)耦合干擾進(jìn)行計(jì)算與預(yù)測(cè)，得到了如圖2、3中虛線所示的補(bǔ)償軌線。由波形分析可知：隨著干擾源向中心領(lǐng)域的移近，傳統(tǒng)最小二乘方法作用下的系統(tǒng)受耦合干擾的影響愈加嚴(yán)重，幅值與相角的最大旁瓣出現(xiàn)了明顯的干擾形變，甚至有影響主頻域波形的趨勢(shì)；應(yīng)用基于極值分解的改進(jìn)算法，不僅可使干擾作用方向上波形不發(fā)生形變，且在主頻域內(nèi)期望軌跡獲得快速響應(yīng)與跟蹤，有效補(bǔ)償了由耦合干擾引起的檢測(cè)誤差，減小了電磁干擾的耦合強(qiáng)度，提高了系統(tǒng)的運(yùn)行效能。

圖2 近域干擾幅值波形對(duì)比圖

圖3 近域干擾相角波形對(duì)比圖

4 結(jié) 語

在分析同塔多回高壓輸電線路電磁干擾耦合特性的基礎(chǔ)上，提出了一種新型的模式綜合電磁兼容預(yù)測(cè)數(shù)值分析方法。將極值分解與矩量法相結(jié)合，通過適當(dāng)配置系統(tǒng)的電磁散射方程，精確預(yù)測(cè)散射干擾的影響，補(bǔ)償耦合干擾引起的跟蹤誤差，減少耦合強(qiáng)度，增強(qiáng)系統(tǒng)抗擾能力，能夠確保系統(tǒng)穩(wěn)定運(yùn)行的同時(shí)，實(shí)現(xiàn)期望幅值、相角的快速跟隨。

由同塔多回高壓輸電線路中電磁兼容預(yù)測(cè)方法的設(shè)計(jì)過程、推導(dǎo)證明和實(shí)驗(yàn)測(cè)試可得到下列結(jié)論：

(1) 存在耦合干擾作用的實(shí)際操作環(huán)境中，采用模式綜合電磁兼容預(yù)測(cè)算法的同塔多回高壓輸電線路運(yùn)行平穩(wěn)，辨識(shí)向量可快速跟蹤期望軌跡，克服了傳統(tǒng)優(yōu)化算法中補(bǔ)償力度不足的缺點(diǎn)，系統(tǒng)穩(wěn)定性高，魯棒性強(qiáng)。

(2) 模式綜合電磁兼容算法結(jié)構(gòu)簡單，運(yùn)算量小，計(jì)算精度高，控制效果好，為分析和設(shè)計(jì)電磁干擾實(shí)時(shí)預(yù)測(cè)方案提供了有效的手段和工具。

(3) 采用模式綜合電磁兼容算法的同塔多回高壓輸電線路具有較好的控制性能和魯棒特性，算法結(jié)構(gòu)不復(fù)雜，實(shí)時(shí)控制中的計(jì)算量也不大，為實(shí)現(xiàn)在線實(shí)時(shí)控制提供了可能。

[1] 孫秋芹,李慶民,呂鑫昌. 同塔多回輸電線路并聯(lián)電抗器的參數(shù)優(yōu)化分析[J]. 中國電機(jī)工程學(xué)報(bào), 2011,31(25):50-57.

SUN Qiu-qin, LI Qing-min, LV Xin-chang. Parameter optimization analysis of shunt reactor of multi-circuit transmission on the same tower[J]. Proceedings of the CSEE, 2011,31(25):50-57.

[2] 亞哈基,阿比利,沙哈巴蒂,等. 基于大重疊域矩量法的頻率選擇表面分析[J]. IEEE 微波變電理論與技術(shù), 2010, 58(8):2175-2187.

Yahaghi A, Abiri H, Shahabadi M,etal. Large Overlapping Subdomain Method of Moments for the Analysis of Frequency Selective Surfaces[J]. IEEE Transactions on Microwave Theory and Techniques, 2010, 58(8):2175-2187.

[3] 韓軍科,楊靖波,李清華. 超/特高壓交流同塔多回輸電線路覆冰不平衡張力分析[J]. 電網(wǎng)技術(shù), 2011,35(12):34-37.

HAN Jun-ke, YANG Jing-bo, LI Qing-hua. Analysis on unbalanced tension caused by ice-coating on conductors of UHV/EHV AC multi-circuit transmission lines on the same tower[J]. Power System Technology, 2011,35(12):34-37.

[4] 顏 偉,趙 陽,王恩榮. 復(fù)雜電子系統(tǒng)輻射EMI 噪聲機(jī)理分析與抑制[J]. 中國電機(jī)工程學(xué)報(bào), 2013,33(15):173-179.

YAN Wei, ZHAO Yang, WANG En-rong. Analysis and suppression of EMI noise mechanism complex electronic systems radiation[J]. Proceedings of the CSEE, 2013,33(15):173-179.

[5] 王玉峰,鄒積巖,廖敏夫. 瞬態(tài)電磁干擾信號(hào)發(fā)生器的研制[J]. 高電壓技術(shù),2009,35(4):866-871.

WANG Yu-Feng, ZOU Ji-yan, LIAO Min-fu. Development of transient electromagnetic interference signal generator[J]. High Voltage Engineering, 2009,35(4):866-871.

[6] 凱伯勒. 未來智能電網(wǎng)的網(wǎng)格化電能質(zhì)量和電磁兼容性[J]. IEEE電磁兼容雜志, 2012(1):100-103.

Keebler P F. Meshing power quality and electromagnetic compatibility for tomorrow’s smart grid[J]. Electromagnetic Compatibility Magazine, IEEE, 2012(1):100-103.

[7] Fujii K, Harada S, Sugiura A. An estimation method for the free-space antenna factor of VHF EMI antennas[J]. IEEE Transactions on Electromagnetic Compatibility, 2005, 47(3): 627-634.

[8] 張向明,趙治華,郭 飛. 電磁兼容測(cè)試系統(tǒng)電磁干擾問題分析與解決[J]. 電工技術(shù)學(xué)報(bào),2010,25(10):14-17.

ZHANG Xiang-ming, ZHAO Zhi-hua, GUO Fei. Analysis and resolution of electromagnetic interference in emc test system[J]. Transactions of China Electrotechnical Society, 2010,25(10):14-17.

[9] 易馬澤,哈姆卡,劉洪江. 基于直接數(shù)據(jù)域最小二乘法的中心厄米特形式的兩步優(yōu)化[C]∥IEEE 系統(tǒng)與控制國際會(huì)議. 土耳其伊斯坦布爾,2010.

Yilmazer N, Bhumkar S, Hongjiang Liu. A two-step optimization of the centro-hermitian form in direct data domain least squares approach[C]. IEEE International Conference on Systems Man and Cybernetics. Istanbul, Turkey.2010.

[10] 馬 超,賈云峰. 基于正交試驗(yàn)設(shè)計(jì)的系統(tǒng)級(jí)電磁自我兼容性的檢測(cè)方法[C]. IEEE第4次國際研討會(huì)議. 中國北京, 2011.

Chao Ma, Yunfeng Jia. The method of system-level electromagnetic self-compatibility inspection based on orthogonal experimental design[C]. 2011 IEEE 4th International Symposium. Beijing, China.

[11] 薛 花,王育飛. 復(fù)雜環(huán)境電磁干擾的解耦控制方法研究[J]. 電機(jī)與控制學(xué)報(bào), 2013,17(7):76-60.

XUE Hua, WANG Yu-fei. Decoupling control method of electromagnetic interference in complex environment[J]. Electric Machines and Control, 2013,17(7):76-60.

[12] 于振偉,杰森米克斯. 基于近場耦合掃描特性化和IC遠(yuǎn)場輻射的一種改進(jìn)偶極矩模型[J]. IEEE 電磁兼容變電, 2013,55(1): 97 - 108.

Zhenwei Yu, Jason A Mix. An Improved Dipole-Moment Model Based on Near-Field Scanning for Characterizing Near-Field Coupling and Far-Field Radiation From an IC[J]. IEEE Transactions on Electromagnetic Compatibility,2013,55(1): 97 - 108.

[13] 里佐利,科斯坦佐,蒙塔納里. 遠(yuǎn)場性自激振蕩微帶天線的耦合電磁/非線性優(yōu)化[C]∥微波研討會(huì)摘要, 美國西雅圖,華盛頓州,2002: 123-126.

Rizzoli V,Costanzo A,Montanari E. Coupled electromagnetic/nonlinear optimization of self-oscillating microstrip antennas with far-field performance specifications[C]∥Microwave Symposium Digest, Seattle, WA, USA,2002: 123-126.

[14] 覃宇建,吳 超,何建國. 航天器電磁兼容預(yù)測(cè)系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)[J]. 系統(tǒng)仿真學(xué)報(bào), 2008,20(18):4859-4861.

TAN Yu-jian, WU Chao, HE Jian-guo. Design and realization of spacecraft electromagnetic compatibility prediction system[J]. Journal of System Simulation, 2011,35(12):34-37.

[15] 李 勃,黃大慶. 一種基于預(yù)測(cè)法的電磁兼容設(shè)計(jì)方法研究[J]. 儀器儀表學(xué)報(bào),2006,27(6):2499-2500.

LI Bo, HUANG Da-qing. Research of an EMC design methodology based on prediction method[J]. Chinese Journal of Scientific Instrument, 2006,27(6):2499-2500.