基于FPGA 的電力系統暫態仿真器研究

顧勇，張永明，強嘉懿

（上海電機學院電氣學院，上海201306）

0 引言

新世紀以來全球經濟飛速發展，對能源的需求量日益擴大，煤、石油以及天然氣等傳統的能源的不斷減少，難以維持各國經濟的長期發展。風能、太陽能由于清潔、成本低、可再生等優點成為當前應用廣泛的可再生新能源，尤其是在發電方面得到了廣泛的大力應用，由此分布式發電技術開始得到廣泛的應用。同時由于當前電力系統的規模越來越龐大，再加上分布式發電的并網運行給電網帶來了越來越多的問題，對于電力系統的仿真的實時性和對電力系統的算法和模型的要求也在不斷地提高，針對第三方硬件而言，一些商業的電力系統仿真軟件如Dspace、RT-LAB 等，并不適用于所有的硬件平臺[1]，并且其軟件的基礎架構以及其所采用的各種器件的模型與仿真算法并不會對用戶開放，利用這些軟件完成相關的仿真并不是困難的，但是如果需要進行基于模型算法層面的深入研究以及基于第三方硬件平臺實現更大規模的、更經濟的仿真研究則是一個不可實現的事情。所以，基于自主開發的電力系統的仿真軟件成為一項基礎且重要的工作。

本文利用Python 語言編寫了大眾化適用性廣的暫態仿真軟件，可以適用于第三方FPGA 硬件平臺，并在建立電力系統的基本元器件模型，選用IEEE 4 節點模型作為實驗模型，完成系統模型的單相接地短路情況下的暫態仿真，對基于FPGA 的仿真平臺進行了驗證。

1 暫態仿真方法

在電力系統仿真領域常用的數值積分算法有：前向歐拉法、后向歐拉法、梯形積分法以及龍格庫塔法。如果對于一個穩定的微分方程一個算法會得到一個穩定的差分方程近似，則稱之為A 類穩定[2]。梯形積分法是一種具有A 類穩定性且具有二階精度的數值算法，所以可以選取梯形積分法作為電力系統電磁暫態仿真的主要算法。

如圖1 所示為一電感及其差分離散形式，其離散形式為電導與歷史項電流的并聯。以此為例，應用梯形積分法后可以得到式（1）（2）（3）的差分方程。

上世紀60 年代，Dommel 教授提出了基于節點電壓法的EMTP 方法作為電磁暫態仿真的基本框架[3]，EM?TP 方法在基爾霍夫定律和電壓電流關系的基礎上，采用數值積分方法將連續的元件離散等效為等值的諾頓電路，采用節點電壓法對線性電路進行求解。EMTP 方法計算效率高穩定性好且計算的精度高，已成為當前主流的仿真方法。節點電壓法的基本形式如式（4）所示。

由于電力系統的模型眾多且復雜，單純的利用基本的節點電壓方法不能夠對電力系統的所有模型都利用EMTP 方法分析，在前人的努力中發展出了改進節點電壓法和增廣的改進節點電壓法[4]。增廣的改進節點電壓法的表示形式如式（5）所示。

增廣的改進節點電壓法相對于改進節點電壓法而言只保留了Yn 的部分，這個方程式可以看作是一個AX=b 的系統，利用這種形式的節點方程方法可以良好地處理含受控源支路（如理想變壓器模型）和不接地的理想電壓源模型以及理想開關等器件的分析。

2 元器件建模

本文中針對電力電子開關器件的建模方法是拓撲建模法。拓撲建模法中的經典成熟的方法有二值電阻法和二值L/C 建模法。二值電阻法的優點是仿真時的精度比較高，但是針對小步長仿真的效果不理想。二值L/C 建模法可以適用于小步長仿真，即當開關閉合時建模為一個值很小的電感，當開關斷開時建模為一個值很小的電容。當選擇合適的電感、電容和仿真步長的數值時，無論開關打開還是閉合都將獲得一個不變的電導的數值，可以保證電路對應的導納矩陣是不變的。二值L/C 建模法開關等效示意圖如圖2 所示。

圖2 二值L/C開關建模示意圖

對于要求不是很高的電力系統的仿真而言，電源類器件可以直接用單相交流源進行模擬，三相電源可分解為3 個單相交流源進行模擬，然后利用相位差區分不同的相別[6]。單相交流源的模擬可以用式（6）表示。

其中，A 是幅值，ω是角頻率，θ是初始相位。

負荷及支路模型的建立要簡單得多，對于多個元件的組合電路可以通過降階處理來降低支路矩陣的階數，從而可以使得在仿真時的仿真步長內的計算量減少，利于快速仿真的實現。圖3 是RLC 支路的降階示意圖。

圖3 RLC支路降階

RLC 電路的統一表達式如式（7）所示：

表1 RLC 電路的統一參數表達

對于輸電線路較近的電力系統模型來說，輸電線路模型采用以集中參數表示的PI 等效模型時足夠精確的[5]，如圖4 所示。在某些情況下，可以忽略線路中的對地電容的影響和RL 串聯阻抗之間的耦合作用，利用多條RL 串聯阻抗來表示表示輸電線路模型。

在電力系統中，變壓器的接線方式會對暫態仿真時的結果有影響，在電力系統中常用的三相變壓器的接線方式有：Yy0 以及Yd11 兩種，如圖5 所示。

圖5 變壓器連接組別

在不考慮變壓器在公共磁路上的耦合與不對稱現象時，可以使用三個單相線性變壓器按照不同的接線方式來組成三相變壓器。圖6 所示為簡化的單相線性變壓器模型。

圖6 單相變壓器模型

3 仿真平臺設計

本仿真平臺的程序軟件基于Python 編程語言開發，模型平臺是MATLAB 2015b，配合以Vivado 2015.4版本，程序框架采用基于節點電壓法的EMTP 框架。Python 語言具有面向對象編程思維方法的特點，具有這種思維方法的編程語言具有易維護且易擴展的特點，在需求變化的時候可以快速簡單有彈性的擴展程序功能用以滿足仿真的需求，同時Python 語言中存在大量的擴展庫，簡單易學，可以降低開發的難度。針對不斷擴大的電力系統模型而言，仿真程序軟件具有良好而擴展性是一個重要的指標。

整個仿真平臺由建模模塊、模型的識別與分析仿真模塊、FPGA 硬件部分以及外部測量儀器如示波器組成。模型仿真模塊基于MATLAB/Simulink 平臺，通過程序對MATLAB/Simulink 平臺的模型進行二次開發，主程序部分識別并獲取模型的電氣信息以及拓撲構成，由主程序生成計算矩陣，并得到Verilog HDL 硬件描述語言的代碼，借助Vivado 軟件將程序燒寫進入FPGA 平臺，獲得模型的仿真波形，完成針對模型的仿真。整個仿真平臺的結構圖如圖7 所示。

圖7 仿真平臺結構圖

程序可以分為五個主要的模塊：數據輸入模塊、模型識別模塊、矩陣生成模塊、開關檢查模塊以及仿真計算模塊。程序各模塊的結構圖如圖8 所示。其中，數據輸入模塊可以基于Simulink 建立的系統模型傳遞給模型識別模塊，以便后續對模型的識別處理；模型識別模塊主要用于識別用戶基于MATLAB/Simulink 建立的系統模型，獲得系統模型的電氣信息及拓撲信息；矩陣形成模塊用于形成系統模型的節點計算矩陣，獲得節點阻抗矩陣、電流源矩陣和歷史項矩陣等；開關檢查模塊主要用于檢查系統模型中的頂管期間的通斷狀態是否發生變化，然后根據結果判斷是否需要修改節點計算矩陣；仿真計算模塊可以根據節點矩陣自動生成Verilog HDL 代碼，調用Vivado 軟件下載到FPGA 中完成系統的仿真計算。

圖8 模塊組成

程序實現的整體思想框架如圖9 所示。程序的核心模塊是接收用戶設定的對仿真的步長等系統仿真的要求信息、自動讀取識別系統模型的各種電氣信息以及拓撲信息并將這些信息保留等待后續處理運行，根據不同元器件的不同的描述方式及關鍵詞對不同的模塊進行適宜的建模方式，并對元器件進行參數初始化、暫態矩陣的形成和電力電子開關器件的狀態的判斷，將獲得的系統模型的信息以及用戶提供信息等仿真數據送入程序的暫態仿真部分，然后通過軟件完成電力系統的仿真，將結果輸出到示波器等外部設備，完成系統的仿真過程。

圖9 程序的架構

程序的計算流程如圖10 所示。

程序的實現首先要讀取模型的數據文件，獲得電力系統模型的元件參數以及拓撲信息。元件信息包括元器件的類型、名稱和基本的參數，拓撲信息包括元器件的連接關系以及各支路和各節點的編號信息。系統的初始化包括對系統中的所包含電感電容元件的諾頓等效電路的生成計算，由初始信息確定模型中的歷史量的賦值信息，完成系統參數的初始化。然后程序可以根據系統模型的初始化信息以及電氣與拓撲信息形成仿真計算所使用的增廣的改進節點法的計算矩陣。矩陣的計算過程放在FPGA 中實現，若要實現這個目的必須把計算過程用硬件描述語言來實現，本文使用的語言為Verilog HDL。程序具備自動將矩陣方程的計算流程用HDL 描述實現的功能，然后調用Vivado 軟件，將HDL 代碼寫入FPGA 中。對于模型中沒有開關類器件時，利用FPGA 對系統模型的計算就可以按照設定步長開始計算。若系統中存在開關器件，那就需要對開關的狀態進行判斷。在開始對系統進行仿真計算以后，需要在每一個計算步長的計算完成以后更新由于離散差分所產生的歷史項電流源的數值，然后更新到計算矩陣中，開始下一個時步的計算。若對系統的仿真設定了計算時間，還需要考慮是否達到總的仿真時間。

圖10 程序計算流程

4 仿真與實驗驗證

圖11 為搭建的實驗平臺，包括上位機、FPGA，示波器以及波形發生器等。IEEE 4 節點模型[7]如圖12 所示。

圖11 實驗仿真平臺

圖12 IEEE 4節點模型

暫態過程的產生可以有多種情況產生，本文選取最常見的一種暫態情況，即發生單相接地短路造成暫態過程。為了不失一般性，本文選擇A 相作為發生單相接地短路的故障相。

本例中，A 相的單相接地故障發生在母線4 上，測量的是母線3 上故障相即A 相的電壓和電流，如圖13所示，藍色的線為電壓波形，綠色的線為電流波形。同樣的選擇程序仿真與RT-LAB 的仿真作為對比。由圖13 可以看出，發生暫態時，程序的仿真波形的毛刺會顯著的增多，由于示波器的濾波的效果的影響也有一定的影響，接地短路發生時，程序的仿真波形會出現一些粗糙，FPGA 的仿真波形的大小比RTLAB 的仿真波形略大，在算法的精確度和在開關動作發生時的對數值的處理會造成波形的差別。仿真程序的電壓與電流的波形趨勢與RT-LAB 仿真情況下略有差異，但都在誤差的允許范圍以內。

5 結語

通過實驗結果可以看出，本文所基于FPGA 和自己開發的程序軟件的電網仿真實驗平臺在對于小型電力系統時的仿真能夠基本準確地得到暫態情況下的仿真表現，對于系統發生暫態情況時的電壓和電流能夠與RTLAB 的結果做到在誤差范圍基本相同。因此，所研究的電網仿真平臺可以實現基本的電網暫態的仿真任務。

圖13 IEEE 4節點模型的暫態仿真