两电平桥式换流器的定导纳建模与实时仿真方法与流程

本发明涉及电力系统，特别是一种两电平桥式换流器的定导纳建模与实时仿真方法。

背景技术：

电磁暂态仿真是电力系统仿真的重要组成部分。其基本理论与方法由加拿大的h.wdommel于20世纪60年代末提出。针对不同类型的应用,电磁暂态仿真可分为离线仿真和实时仿真。通常来说,离线仿真工具的计算时间要远多于所研究暂态现象的持续时间。而在面向对时间有严格要求的应用场景时，实时仿真通过软硬件平台相互配合，保证了仿真器内部时钟与现实世界时钟的精确同步，可以为各种电力系统保护与控制装置提供高度模拟现场实际的测试环境。

随着越来越多的电力电子设备引入电力系统，电力电子开关的高频离散特性给电力电子设备的建模与实时仿真的实现带来了巨大挑战。目前，在电磁暂态仿真中，电力电子开关主要采用的建模方法可以分为以下两种：

1)二值电阻建模，即开关导通时用小电阻等效，关断时用大电阻等效；

2)基于电感/电容等效的定导纳建模，即开关导通时用小电感等效，关断时用小电容等效；

二值电阻建模在开关状态变化时将引起开关支路的导纳突变，因此每次开关动作都需要重新形成导纳矩阵，效率低下，难以满足实时性要求，更多地应用于离线的电磁暂态仿真工具中，如pscad-emtdc，matlab的simpowersimulationtoolbox，emtp系列仿真软件等；

基于电感/电容等效的定导纳建模可以通过合理的参数设置，使小电感和小电容的等效导纳相等，避免了因开关动作导致的导纳矩阵改变，极大地提高了仿真效率，被应用于实时数字仿真仪rtds的小步长模型库中。然而，受电感和电容物理特性的制约，这种定导纳模型在开关动作后存在明显暂态误差，仿真得到的换流器功率损耗远大于实际情况，严重影响了仿真精度，这一现象被称为虚拟功率损耗问题；

技术实现要素：

针对已有两电平桥式换流器建模方法的不足，本发明的目的在于提供一种专门用于两电平桥式换流器的定导纳建模与仿真方法，该方法在避免因开关动作导致的导纳矩阵改变的同时，解决了传统建模仿真方法在两电平桥式换流器实时仿真时的虚拟功率损耗问题，极大地提高了仿真精度。

本发明的技术解决方案如下：

一种两电平桥式换流器的定导纳建模与实时仿真方法，其特点在于该方法，包括下列步骤：

1)将两电平桥式换流器及其所在电路中的支路和节点分别进行编号，其中接地节点的编号为0；

2)将各个电阻支路、电感支路、电容支路和开关支路，分别用一个等效导纳和一个历史电流源并联的等效模型替换，独立电压源支路用一个等效导纳和一个等效电流源并联的等效模型替换，各类支路的等效导纳计算公式如下：

电阻支路的等效导纳其中，r是电阻支路的电阻值；

电感支路的等效导纳其中，l是电感支路的电感值，δt是实时仿真的时间步长；

电容支路的等效导纳其中，c是电容支路的电容值；

开关支路的等效导纳其中，cdc是两电平桥式换流器直流侧的电容，lac是两电平桥式换流器交流侧的电感；

独立电压源支路的等效导纳其中，rs是独立电压源支路的内电阻；

3)根据各个支路的等效导纳，计算被仿真电路的节点导纳矩阵(导纳矩阵的计算方式可参考电路理论相关的书籍)；

4)如果当前仿真时刻t是仿真初始时刻，则所有支路的历史电流源大小取零，独立电压源支路的等效电流另外计算，如果不是仿真初始时刻，则根据上一仿真时刻各个支路的支路电压和支路电流，计算出当前仿真时刻各个支路的历史电流源大小，独立电压源支路的等效电流另外计算，各类支路的历史电流和等效电流的计算公式如下：

电阻支路的历史电流ih_r＝0；

电感支路的历史电流ih_l(t)＝il(t-δt)

其中，il(t-δt)是上一仿真时刻电感支路的支路电流；

电容支路的历史电流：ih_c(t)＝-yb_cuc(t-δt)

其中，uc(t-δt)是电容支路上一仿真时刻的支路电压，yb_c是电容支路的等效导纳；

导通时，开关支路的历史电流：

关断时，开关支路的历史电流：或者采用

导通时，开关支路的历史电流

关断时，开关支路的历史电流其中,usw(t-δt)是上一仿真时刻开关支路的支路电压，isw(t-δt)是上一仿真时刻开关支路的支路电流，yb_sw是开关支路的等效导纳；

独立电压源支路的等效电流源：其中,vs(t)是独立电压源支路的内电势大小，rs是独立电压源支路的内电阻大小；

5)根据当前仿真时刻的历史电流源和等效电流源的大小，计算流入各个节点的注入电流大小；

6)根据流入各个节点的注入电流，结合节点导纳矩阵，计算出各个节点的电压大小；

7)根据各个节点的电压，结合各个支路的等效导纳，计算出各个支路电压和支路电流；

8)如果没有到达最后的仿真时刻，令下一仿真时刻为t+δt，并返回步骤4)，当到达最后的仿真时刻，结束。

本发明的技术效果:

1)本发明在两电平桥式换流器中的开关支路等效模型在导通状态下和关断状态下采用相同的等效导纳，避免了仿真过程中因开关状态改变而重新形成导纳矩阵的操作，使得开关动作后的暂态误差快速衰减。

2)采用本发明方法的实时仿真波形比传统的基于电感电容等效法更贴近理想开关波形，极大地提高了两电平桥式换流器的实时仿真精度。对于理想开关组成的换流器而言，没有虚拟功率损耗，而在传统的基于电感/电容等效的仿真方法下，两电平换流器的虚拟功率损耗随着开关频率的增大而增大，如图3所示，在100khz时高达40％以上，和实际严重不符。在本发明方法下，两电平换流器的虚拟功率损耗基本不随开关频率变化，始终维持在接近零的水平，更加贴近理想开关组成的换流器。

附图说明

图1是支路的等效模型示意图；

图2是本发明方法与传统实时仿真方法结果对比；

图3是两种方法在不同换流器开关频率下的虚拟功率损耗率；

图4是含单相两电平桥式换流器的简单电路；

图5是本发明两电平桥式换流器的定导纳建模与实时仿真方法流程图。

具体实施方式

为便于理解，下面将以一个如图4所示的含单相两电平桥式换流器的简单电路为例，对本发明进行阐述，但不应以此限制本发明的保护范围。

在具体实现该两电平桥式换流器的实时仿真时，本发明采用硬件平台：pxie机箱内分别装有美国国家仪器(ni)公司的pxie-8135(pxie控制器)和pxie-7975r(fpga模块)，pxie控制器主要负责换流器控制系统的仿真，fpga模块主要负责两电平桥式换流器电路部分的仿真，它们之间通过pxie总线进行通信。此外，pxie控制器还可以通过以太网与上位机进行通信，在上位机上显示实时仿真波形，fpga模块可以通过i/o端口与外置控制器和示波器相连，进行硬件在环仿真。

上位机、pxie控制器和fpga模块中的程序都是通过美国国家仪器(ni)公司的labview开发环境进行统一编程。通过labview编程，上位机中的程序完成与pxie控制器的通讯以及仿真波形显示等功能；pxie控制器中的程序完成与上位机的通信、从fpga模块中读取和写入数据、模拟换流器的控制系统等功能。上述程序不属于本发明保护范围，且美国国家仪器(ni)公司官网上提供相关的程序范例，故不再详细描述。而fpga模块是本发明的具体实施载体，通过labview编程，参见图5，图5是本发明两电平桥式换流器的定导纳建模与实时仿真方法流程图。

本发明两电平桥式换流器的定导纳建模与实时仿真方法实施例的步骤如下：

(1)将两电平桥式换流器及其所在电路中的支路和节点分别进行编号，如图所4示；

(2)将各个电阻支路、电感支路、电容支路以及开关支路，分别用一个等效导纳和一个历史电流源并联的如图1所示的等效模型替换，独立电压源支路用一个等效导纳和一个等效电流源并联的等效模型替换，仿真步长取1μs，则各个支路的等效导纳如下：

支路1(独立电压源支路)

支路2(电容支路)

支路3(电容支路)

支路4(开关支路)

支路5(开关支路)

支路6(电感支路)

支路7(电阻支路)计算得到的各个支路的等效导纳值也标注在图4b中；

(3)根据各个支路的等效导纳，计算被仿真电路的节点导纳矩阵yn

(4.0)当前仿真时刻t＝0μs，是仿真初始时刻，则所有支路的历史电流源大小取零，独立电压源支路的等效电流另外计算，即：

支路1(独立电压源支路)

支路2(电容支路)ih_cdc1＝0a

支路3(电容支路)ih_cdc2＝0a

支路4(开关支路)ih_sw1＝0a

支路5(开关支路)ih_sw2＝0a

支路6(电感支路)ih_l＝0a

支路7(电阻支路)ih_r＝0a

(5.0)根据当前仿真时刻的历史电流源和等效电流源的大小，计算流入各个节点的注入电流大小(流入为正，流出为负)：

节点1的注入电流in1＝7500a

节点2的注入电流in2＝-7500a

节点3的注入电流in3＝0a

节点4的注入电流in4＝0a

(6.0)已知流入各个节点的注入电流，结合节点导纳矩阵，根据节点电压方程ynvn＝in，计算出各个节点的电压大小:

节点1的电压vn1＝3.7125v

节点2的电压vn2＝-3.7125v

节点3的电压vn3＝0v

节点4的电压vn4＝0v

(7.0)根据各个节点的电压，结合各个支路的等效导纳，计算出各个支路电压和支路电流：

支路1(独立电压源支路)vb_vs＝7.4250v,ib_vs＝7.4257e3a

支路2(电容支路)vb_cdc1＝3.7125v,ib_cdc1＝7.4250e3a

支路3(电容支路)vb_cdc2＝3.7125v,ib_cdc2＝7.4250e3a

支路4(开关支路)vb_sw1＝3.7125v,ib_sw1＝0.7425a

支路5(开关支路)vb_sw2＝3.7125v,ib_sw2＝0.7425a

支路6(电感支路)vb_lac＝0v,ib_lac＝0a

支路7(电阻支路)vb_r＝0v,ib_r＝0a

(8.0)当前仿真时刻t＝0μs，没有到达最后的仿真时刻，返回步骤(4)，从而进入下一仿真时刻t＝1μs；

(4.1)当前仿真时刻t＝1μs，不是仿真初始时刻，则根据各个支路上一仿真时刻的支路电压和支路电流，计算出当前仿真时刻各个支路的历史电流源大小，各类支路的历史电流计算公式如下：

支路1(独立电压源支路)

支路2(电容支路)ih_cdc1＝-7.4250e3a

支路3(电容支路)ih_cdc2＝-7.4250e3a

支路4(导通的开关支路)ih_sw1＝2.5351a

支路5(关断的开关支路)ih_sw2＝-0.4349a

支路6(电感支路)ih_l＝0a

支路7(电阻支路)ih_r＝0a；

(5.1)根据当前仿真时刻的历史电流源和等效电流源的大小，计算流入各个节点的注入电流大小；

节点1的注入电流in1＝1.4922e4a

节点2的注入电流in2＝-1.4925e4a

节点3的注入电流in3＝2.9700a

节点4的注入电流in4＝0a；

(6.1)根据流入各个节点的注入电流，结合节点导纳矩阵，计算出各个节点的电压大小：

节点1的电压vn1＝7.3874v

节点2的电压vn2＝-7.3874v

节点3的电压vn3＝7.4213v

节点4的电压vn4＝0.0148v；

(7.1)根据各个节点的电压，结合各个支路的等效导纳，计算出各个支路电压和支路电流：

支路1(独立电压源支路)vb_vs＝14.7747v,ib_vs＝7.3523e3a

支路2(电容支路)vb_cdc1＝7.3874v,ib_cdc1＝7.3497e3a

支路3(电容支路)vb_cdc2＝7.3874v,ib_cdc2＝7.3497e3a

支路4(开关支路)vb_sw1＝-0.0339v,ib_sw1＝2.5283a

支路5(开关支路)vb_sw2＝14.8087v,ib_sw2＝2.5268a

支路6(电感支路)vb_lac＝7.4065v,ib_lac＝0.0015a

支路7(电阻支路)vb_r＝0.0148v,ib_r＝0.0015a；

(8.1)当前仿真时刻t＝1μs，没有到达最后的仿真时刻，返回步骤(4)，从而进入下一仿真时刻t＝2μs；

(4.2)当前仿真时刻t＝2μs，不是仿真初始时刻，则根据各个支路上一仿真时刻的支路电压和支路电流，计算出当前仿真时刻各个支路的历史电流源大小……

如上重复步骤(4)～(8),直到达到最后的仿真时刻，结束仿真程序。

实施结果请参见图2、图3，图2是本发明方法与传统实时仿真方法结果对比，图3是两种方法在不同换流器开关频率下的虚拟功率损耗率。如图2所示，采用本发明方法的实时仿真波形比传统的基于电感电容等效法更贴近理想开关波形，极大地提高了两电平桥式换流器的实时仿真精度。对于理想开关组成的换流器而言，没有虚拟功率损耗，而在传统的基于电感/电容等效的仿真方法下，两电平换流器的虚拟功率损耗随着开关频率的增大而增大，如图3所示，在100khz时高达40％以上，和实际严重不符。在本发明方法下，两电平换流器的虚拟功率损耗基本不随开关频率变化，始终维持在接近零的水平，更加贴近理想开关组成的换流器。