周艷華,聶冰青

(1.東北電力大學(xué) 電氣工程學(xué)院,吉林 吉林 132012； 2.國電南京自動化股份有限公司,南京 211100)

?

功率連接型數(shù)?；旌戏抡嫦到y(tǒng)接口特性分析

周艷華1,聶冰青2

(1.東北電力大學(xué) 電氣工程學(xué)院,吉林 吉林 132012； 2.國電南京自動化股份有限公司,南京 211100)

為解決功率連接型數(shù)?；旌戏抡嬷虚]環(huán)穩(wěn)定性與仿真精確性兩個關(guān)鍵問題,提出基于系統(tǒng)傳遞函數(shù)的穩(wěn)定性分析方法和精確性評估方法。介紹并比較了5種接口算法的特性,研究論證了其穩(wěn)定性與精確性的優(yōu)劣。在Matlab/Simulink中搭建了簡單電路數(shù)字仿真子系統(tǒng)與物理被試系統(tǒng)的模型,對5種接口算法的穩(wěn)定性與精確性進行了離線仿真,結(jié)果表明特定應(yīng)用中阻尼阻抗法相比于其他算法可行性更高。

功率連接型數(shù)?；旌戏抡?接口算法;穩(wěn)定性;精確性

隨著計算機技術(shù)的快速發(fā)展,數(shù)?；旌戏抡?hardware-in-the-loop, HIL)技術(shù)在電力電子及電力系統(tǒng)領(lǐng)域得到越來越廣泛的應(yīng)用。由于數(shù)?；旌戏抡婕夹g(shù)是一種先進的實驗/測試方法,因此可在虛擬環(huán)境中反復(fù)、安全、經(jīng)濟地模擬新型的物理裝置原型[1]。

目前大部分數(shù)?；旌戏抡鎽?yīng)用均以控制裝置作為物理被試系統(tǒng)。數(shù)字仿真系統(tǒng)與物理被試系統(tǒng)之間傳輸?shù)男盘枮榈凸β实目刂菩盘?變化范圍為+10～-10 V,由A/D、D/A裝置準(zhǔn)確轉(zhuǎn)換,稱為信號型數(shù)模混合仿真(control hardware—in-the-loop, CHIL)。而實際應(yīng)用中一些功率設(shè)備,如電動機也需要被測試。這種情況下的物理被試系統(tǒng)會發(fā)出或吸收有功功率,此時接口設(shè)備需包含功率放大及功率轉(zhuǎn)換裝置,稱之為功率連接型數(shù)?；旌戏抡?power hardware-in-the-loop, PHIL)[2]。

在功率連接型數(shù)?；旌戏抡嫦到y(tǒng)中,尤其是對于高功率設(shè)備,功率接口帶來的誤差(如延時和失真等)會引起嚴重不穩(wěn)定問題和不準(zhǔn)確的仿真結(jié)果。文獻[3-6]介紹了將功率連接型數(shù)?；旌戏抡孀鳛橐环N便捷工具來解決電力系統(tǒng)各種問題的初步試驗,而未研究接口算法對系統(tǒng)穩(wěn)定性的影響。文獻[7]提出了一種基于一階近似的新型接口算法,但未對該算法的穩(wěn)定性進行詳盡的理論分析。文獻[8-9]提出一種一般性的實驗步驟,然而只停留在理論研究層面,很難應(yīng)用于實踐中。以上研究均針對獨立系統(tǒng)和特定問題,對如何實現(xiàn)功率連接型數(shù)模混合仿真,尤其是接口算法的選擇缺乏一般性的分析。因此,本文論述了引起功率連接型數(shù)?；旌戏抡娌环€(wěn)定問題的根本原因,并比較了常用的5種接口算法,提出只有選擇合適的接口算法才能得到最優(yōu)仿真結(jié)果。

1 功率連接型數(shù)?；旌戏抡娌环€(wěn)定問題

閉環(huán)穩(wěn)定性是功率連接型數(shù)模混合仿真中的關(guān)鍵問題。仿真的不穩(wěn)定性不僅導(dǎo)致仿真結(jié)果錯誤,甚至?xí)斐蓽y試設(shè)備和物理被試系統(tǒng)的損壞。引起不穩(wěn)定的主要原因包括接口延時、有限帶寬、接口功率放大器產(chǎn)生的諧波等。以一個簡單的功率連接型數(shù)模混合仿真系統(tǒng)為例(如圖1所示),對其穩(wěn)定性進行分析。

圖1 數(shù)?；旌戏抡嫦到y(tǒng)不穩(wěn)定問題分析電路

原始電路為分壓式電路,PHIL電路中,負載阻抗ZL為實際的硬件電阻,而電壓源VS和阻抗ZS都是由數(shù)字仿真實現(xiàn)。為了實現(xiàn)這個接口,采用一個電壓放大器將“數(shù)字電壓V1”變?yōu)椤皩嶋H的電壓V2”,并加在負載電阻兩端。流過負載電阻的實際電流I2通過傳感器測量并反饋至數(shù)字側(cè),即數(shù)字側(cè)電流源I1。盡管原始分壓電路穩(wěn)定,但解耦后的PHIL電路并不穩(wěn)定。

假設(shè)時間tk時,V2的電壓放大器產(chǎn)生一個誤差ε,即Δv2=ε。I2=V2/ZL,則電流I2對應(yīng)的誤差為Δi2(tk)=ε/ZL。當(dāng)電流Δi2反饋回數(shù)字側(cè)時,進一步導(dǎo)致v1的更大誤差,即

Δi2(tk)=ε/ZL

V1=Vs-Zs×I1

Δv1(tk+1)=-(Zs/ZL)ε

下一仿真步長更新后的V1值包含上一步長的誤差Δv1,Δv1被進一步放大,系數(shù)為-(ZS/ZL)。如果ZS/ZL>1,誤差的幅值將不斷增大,最終導(dǎo)致系統(tǒng)失穩(wěn)。

1.1 功率連接型數(shù)?；旌戏抡娣€(wěn)定性分析

功率連接型數(shù)?；旌戏抡嫦到y(tǒng)可表示成傳遞函數(shù)的形式,如圖2所示。

圖2 傳遞函數(shù)形式的功率連接型數(shù)?；旌戏抡嫦到y(tǒng)

為保證功率連接型數(shù)?；旌戏抡嫦到y(tǒng)穩(wěn)定,數(shù)字仿真子系統(tǒng)與物理被試子系統(tǒng)都必須是穩(wěn)定的系統(tǒng),且反饋回路必須只含有負實部的極點,即開環(huán)傳遞函數(shù)GOL=T12TBT22TF需滿足奈奎斯特穩(wěn)定判據(jù)。

將上述分析應(yīng)用于圖1所示系統(tǒng),假設(shè)功率接口為理想接口,電壓放大器以一個微小延時Δt來表示,可得到以下關(guān)系:

u1=Vs,T11=1,T12=-Zs,z1=V1

r2=V2,TF=e-sΔt,u2=0,T21=0

T22=1/ZL,z2=I2,r1=I1,TB=1

代入開環(huán)傳遞函數(shù)表達式可得

ZS/ZL的值為2,e-sΔt代表當(dāng)系統(tǒng)頻率增加時一個固定的相位差。開環(huán)傳遞函數(shù)的奈奎斯特曲線為一組半徑不同的圓,無限次包含(-1,0)點,表明系統(tǒng)不穩(wěn)定。

1.2 功率連接型數(shù)?；旌戏抡嫦到y(tǒng)精確性分析

穩(wěn)定的功率連接型數(shù)?；旌舷到y(tǒng)并不能確保該系統(tǒng)在實際應(yīng)用中可行。除穩(wěn)定之外,必須保證該系統(tǒng)精確。由于缺乏判斷精確性的參考標(biāo)準(zhǔn),因此判斷一個系統(tǒng)是否精確較為困難。而評估功率連接型數(shù)?；旌戏抡嫦到y(tǒng)精確性的一個可行方法是評估考慮接口擾動時系統(tǒng)傳遞函數(shù)的響應(yīng)。

假設(shè)圖1所示系統(tǒng)中v1代表系統(tǒng)的狀態(tài)。由于絕大多數(shù)系統(tǒng)的誤差來源于接口電壓放大器,將電壓誤差視為擾動并研究擾動與v1之間的傳遞函數(shù),即

分析可知,正常帶寬范圍內(nèi)TF接近于單位增益,系統(tǒng)的開環(huán)傳遞函數(shù)GOL幅值越小,系統(tǒng)越精確。

2 接口算法

將圖1所示原始分壓電路按圖3所示系統(tǒng)進行數(shù)?；旌戏抡?。放大的信號不再是v1,而是電流i1。測量加在負載上的電壓v2并反饋至數(shù)字仿真系統(tǒng)。

此時系統(tǒng)的開環(huán)傳遞函數(shù)為

GOL的幅值為0.5,系統(tǒng)穩(wěn)定。表明采用不同結(jié)構(gòu)的接口拓撲可提高系統(tǒng)的精確性與穩(wěn)定性。為區(qū)分不同拓撲結(jié)構(gòu)的接口,引入接口算法的概念。接口算法定義了兩個重要的特征:一是傳輸信號的類型,如電流、電壓、功率等;二是傳輸信號的方式,如增益、低通濾波等。

圖3 放大接口電流的數(shù)模混合仿真系統(tǒng)

2.1 理想變壓器法(ITM)

理想變壓器法是進行功率連接型數(shù)?；旌戏抡孀顐鹘y(tǒng)、最直接的方法之一。根據(jù)放大信號的不同分為電壓型理想變壓器法(見圖1)和電流型理想變壓器法(見圖3)。假設(shè)PHIL系統(tǒng)誤差只來源于接口放大設(shè)備的延時,則基于ITM法的系統(tǒng)開環(huán)傳遞函數(shù)可表示為

其中,ZL和ZS分別為物理被試系統(tǒng)和數(shù)字仿真系統(tǒng)的等效阻抗。該等效在實際中是可行的,因為任何時刻電力系統(tǒng)都可以表示成戴維南或諾頓形式。該接口的穩(wěn)定性取決于ZS/ZL的值。

2.2 時變一階近似法(TFA)

TFA是基于假設(shè)PHIL系統(tǒng)物理被試系統(tǒng),可以視為一個一階RL電路或一階RC電路。結(jié)合以往歷史實驗數(shù)據(jù),該模型的回歸系數(shù)可以在線升級解決。在數(shù)字側(cè)進行補償以糾正接口帶來的誤差。

以RL電路為例,假設(shè)物理被試系統(tǒng)存在下述一階等式:

根據(jù)梯形近似法,電流i2的值可由上一仿真步長中的電壓v1和電流i2計算得到,即

i2(k)≈αv1(k-1)+βi2(k-1)

i2(k)=Geqv1(k-1)+Ieq

α、β可由下式計算:

RL型時變一階近似法結(jié)構(gòu)如圖4所示。

圖4 時變一階近似法接口結(jié)構(gòu)

RC型時變一階近似法與RL型特性相似。TFA算法本質(zhì)上是一種預(yù)測方法,通過建立物理被試系統(tǒng)的模型來預(yù)測系統(tǒng)狀態(tài)的變化。因此,和其他任何一種預(yù)測型方法一樣,時變一階近似法在處理非線性系統(tǒng)及高頻信號時存在缺陷。

2.3 傳輸線模型法(TLM)

輸電線模型是利用連接電感或電容對輸電線路進行等效的一種方法。由于傳輸線可進一步用Bergeron模型表示(從電路任何一端看,線路為獨立的戴維南或諾頓電路,只能看到另一端當(dāng)前的參數(shù)),因此該方法適用于解耦大規(guī)模電力系統(tǒng),且方便進行計算。

應(yīng)用傳輸線模型的功率連接型數(shù)?；旌戏抡嫦到y(tǒng)結(jié)構(gòu)如圖5所示。

圖5 傳輸線法接口結(jié)構(gòu)

將連接電容或電感以傳輸線路的Bergeron模型代替,等效電阻R1k可由下式計算得到:

式中Δt為線路行波傳輸時間,在功率連接型數(shù)?；旌戏抡嬷凶兂闪私涌谘訒r。

傳輸線模型算法的等效傳遞函數(shù)為

由于TLM是嚴格基于梯形近似的,因此它非常穩(wěn)定,然而在功率連接型數(shù)?；旌戏抡嬷泻茈y實現(xiàn)。首先,該算法以一個電阻R1k代替了連接電容或電感,在高功率的數(shù)?；旌戏抡嬷?電阻將消耗大量功率。其次,R1k的參數(shù)與連接電容或電感的參數(shù)有關(guān),當(dāng)仿真系統(tǒng)發(fā)生變化時,R1k值隨其變化,以致實現(xiàn)不夠靈活,耗費巨大成本。最后,接口延時Δt與負載狀態(tài)和信號的頻率有關(guān),假設(shè)Δt為固定值,則會導(dǎo)致仿真精度降低。

2.4 部分電路復(fù)制法(PCD)

PCD法來源于稀疏技術(shù)。文獻[10]將該方法應(yīng)用于大型電路仿真軟件SPICE。PCD法將原始電路分割成多個子電路,使用迭代法求解。部分電路復(fù)制法的結(jié)構(gòu)如圖6所示。原始電路中的連接阻抗zab在數(shù)字側(cè)和物理側(cè)被重復(fù)。

圖6 部分電路復(fù)制法接口結(jié)構(gòu)

PCD法的開環(huán)傳遞函數(shù)為

(1)

由式(1)可知,對于電阻性網(wǎng)絡(luò),部分電路復(fù)制法比理想變壓器法的穩(wěn)定性更好,原因在于PCD法更容易使開環(huán)傳遞函數(shù)的幅值小于1。假設(shè)該方法在特定的應(yīng)用中收斂,仿真精度可由足夠的迭代次數(shù)保證。但在實際仿真中,每個仿真步長只能迭代一次,因此每次迭代的誤差必須盡可能小。為保證仿真精度,連接阻抗zab的數(shù)值應(yīng)比za和zb的值大。實際應(yīng)用很難滿足上述條件,因此部分電路復(fù)制法的精度一般較低。

2.5 阻尼阻抗法(DIM)

阻尼阻抗法是介于理想變壓器法和部分電路復(fù)制法之間的一種方法,通過插入一個阻尼阻抗z*來實現(xiàn)。阻尼阻抗法的結(jié)構(gòu)如圖7所示。

圖7 阻尼阻抗法結(jié)構(gòu)

由圖7可知,當(dāng)阻尼阻抗z*等于0時,電壓v*等于v1,DIM法與PCD法相同;當(dāng)z*無窮大時,DIM法變成了ITM法。根據(jù)阻尼阻抗z*的取值不同,DIM法的穩(wěn)定性介于部分電路復(fù)制法和理想變壓器法之間。z*的最優(yōu)取值可由DIM法的開環(huán)傳遞函數(shù)求得:

(2)

由式(2)可知,當(dāng)z*=Zb時,開環(huán)傳遞函數(shù)的幅值為零,則系統(tǒng)高度穩(wěn)定。然而,很難獲得Zb的準(zhǔn)確值,因為這要求物理被試系統(tǒng)為理想模型,所以實際應(yīng)用中并不存在。但可由上一步仿真結(jié)果計算出Zb的近似值,可由仿真得到的v2′和i2的有效值計算出Zb的平均值。盡管阻抗元件的增加會帶來損耗,但阻尼阻抗法仍很大程度地改善了系統(tǒng)穩(wěn)定性和精確性。

3 仿真驗證

3.1 仿真算例

以圖8所示的簡單電路為例,在Matlab/Simulink中分別搭建了數(shù)字仿真子系統(tǒng)與物理被試子系統(tǒng)模型,進行了數(shù)模混合仿真的離線驗證,并與理論分析進行了對比。仿真步長設(shè)為60 μs,接口延時取500 μs。

圖8 用于數(shù)?；旌戏抡娴暮唵坞娐?/p>

系統(tǒng)穩(wěn)態(tài)運行時的波形如圖9所示。

圖9 RS=2 Ω時,不同接口算法的數(shù)?；旌戏抡鎸Ρ葓D及局部放大圖

由圖9可知,傳輸線法和部分電路復(fù)制法波形誤差較大,而其他算法的結(jié)果相對準(zhǔn)確。當(dāng)RS的值增加到5 Ω時,傳輸線法和理想變壓器法使系統(tǒng)仿真不穩(wěn)定。同時傳輸線法和部分電路復(fù)制法的仿真結(jié)果仍然誤差較大,而阻尼阻抗法依然保持較高的穩(wěn)定性與精度,如圖10所示。

3.2 分析對比

通過理論分析與仿真驗證,將5種接口算法的特性總結(jié)如下:

圖10RS=5 Ω時,不同接口算法下的數(shù)模混合仿真對比圖及局部放大圖

Fig.10RS=5 Ω, comparison chart and partial discharge of different interface algorithm for digital analog mixed simulation

1) 理想變壓器法作為最傳統(tǒng)、最直接的方法,精度高但穩(wěn)定性相對較低。

2) 傳輸線法和部分電路復(fù)制法具有較高的穩(wěn)定性,但其精確性過度依賴于連接元件的參數(shù)及電源電阻和負載電阻的比例。此外,傳輸線模型要求原始電路中存在連接電容或電感進行等效解耦,并在數(shù)?；旌戏抡娴臄?shù)字仿真系統(tǒng)及物理被試系統(tǒng)中加入電阻元件,在許多應(yīng)用中,電阻元件的增加限制了該方法的可行性。

3) 時變一階近似法使用預(yù)測方法以補償接口延時,但該方法穩(wěn)定性較低,且接口中存在大量噪聲,實現(xiàn)起來不夠靈活。

4) 阻尼阻抗法結(jié)合了理想變壓器法和部分電路復(fù)制法的優(yōu)點。當(dāng)z*接近物理被試系統(tǒng)的阻抗值時,該方法具有較高的穩(wěn)定性與精確性,適用于絕大多數(shù)功率連接型數(shù)?；旌戏抡嫦到y(tǒng)。

4 結(jié) 語

為保證功率連接型數(shù)?；旌戏抡娴目煽啃耘c安全性,本文提出通過奈奎斯特穩(wěn)定判據(jù)判斷系統(tǒng)穩(wěn)定性方法。詳細分析并比較了5種接口算法特性。在Matlab/Simulink中搭建了簡單電路數(shù)字仿真子系統(tǒng)與物理被試系統(tǒng)的模型,對5種接口算法的穩(wěn)定性與精確性進行了離線仿真。結(jié)果表明:盡管不存在適用于所有功率連接型數(shù)?；旌戏抡娴淖顑?yōu)接口算法,但與其他算法相比,阻尼阻抗法的穩(wěn)定性與精確性更高。

[1] 周俊,郭劍波,朱藝穎,等.特高壓交直流電網(wǎng)數(shù)?；旌蠈崟r仿真系統(tǒng)[J].電力自動化設(shè)備,2011,31(9):18-21.

ZHOU Jun,GUO Jianbo,ZHU Yiying, et al.Digital-analog hybrid real time simulation system of UHVDC and UHVAC power grid[J]. Electric Power Automation Equipment 2011,31(9):18-21.

[2] 胡昱宙,張沛超,包海龍,等.功率連接型數(shù)字物理混合仿真系統(tǒng)(一)接口算法特性[J].電力系統(tǒng)自動化,2013,7(37):36-42.

HU Yuzhou,ZHANG Peichao,BAO Hailong, et al.Power hardware-in-the-loop simulation system. part one characteristics of interface algorithms[J]. Automation of Electric Power Systems,2013,7(37:):36-42.

[3] ACHA E, ANAYA-LARA O, PARLE J, et al. Real-time simulator for power quality disturbance applications[J].Harmonics Quality Power,2000,4(3): 763-768.

[4] ZHU W, PEKAREK S, JATSKEVICH J, et al. A model-in-the-loop interface to emulate source dynamics in a zonal DC distribution system[J].IEEE Trans.Power Eletron,2005,2(20):438-445.

[5] VIEHWEIDER A, LAUSS G, FELIX L. Stabilization of power hardware-in-the-Loop simulations of electric energy syste-ms[J]. Simulation Modelling Practice and Theory, 2011,19(7):1699-1708.

[6] REN W, STEURER M, BALDWIN T L. An effective method for evaluating the accuracy of power hardware-in-the-loop simulations[J]. IEEE Transactions on Industry Applications,2009,45(4):1484-1490.

[7] 劉云,蔣衛(wèi)平,印永華,等.高壓交直流大電網(wǎng)的數(shù)?；旌蠈崟r仿真系統(tǒng)建模[J].力系統(tǒng)自動化,2008,32(12):52-56.

LIU Yun, JIANG Weiping, YIN Yonghua, et al.Modeling of analogue-digital hybrid real-time simulation system applied in the UHV AC/DC great power grid[J].Automation of Electric Power Systems,2008,32(12):52-56.

[8] 陳磊,閔勇,葉駿,等.數(shù)字物理混合仿真系統(tǒng)的建模及理論分析[J].電力系統(tǒng)自動化,2009,33(23):9-13.

CHEN Lei,MIN Yong,YE Jun, et al. Modeling and theoretical analysis of Hardware-in-the-loop Simulation[J]. Automation of Electric Power Systems,2009,33(23):9-13.

[9] DMITRIEV-ZDOROV V B. Generalized coupling as a way to improve the convergence in relaxation-based sovers[J].Proc.Euro Daceuro VHDL,1996,16(20):15-20.

(責(zé)任編輯 侯世春)

Analysis of power hardware-in-the-loop simulation system interface characteristics

ZHOU Yanhua1, NIE Bingqing2

(1.College of Electrical Engineering, Northeast Dianli University,Jilin 132012,China; 2.Guodian Nanjing Automation CO,Ltd., Nanjing 211100, China)

The closed-loop stability and the simulation accuracy are two paramount issues in power hardware-in-the-loop (HIL) simulation. In order to solve the issues, the stability analysis and accuracy evaluation based on system transfer function were proposed. Five different interface algorithms were described, and their respective characteristics with respect to the system stability were compared. Through Matlab/Simulink simulations, the models of simple circuit digital simulation subsystem and physical tested system were established to simulate the stability and accuracy of the five interface algorithms offline. The results show the feasibility of damping and impedance method is higher than that of others in specific application.

power hardware-in-the-loop(HIL) simulation; interface algorithm; stability; accuracy

2015-06-23。

周艷華(1989—)，女，碩士研究生，研究方向為電力系統(tǒng)仿真。

TM743

A

2095-6843(2016)01-0028-06