基于PSCAD的同步发电机组自定义建模及初始化方法

基于pscad的同步发电机组自定义建模及初始化方法
技术领域
1.本发明涉及电磁暂态仿真技术领域，具体为基于pscad的同步发电机组自定义建模及初始化方法。


背景技术：

2.随着新能源产业的发展以及电网规模的不断扩大，大量电力电子设备、高压直流输电系统接入电网，这种以新能源为主体的新型电力系统对仿真提出了个更高的要求。我国电力系统传统以机电暂态仿真为主的现状，难以满足对电网精细化建模的需求。由此便催生了机电暂态仿真模型向电磁暂态仿真模型的转换，电力系统将迈入大规模电磁暂态仿真的新时代。但是相比于机电暂态仿真，电磁暂态仿真的仿真规模更小。
3.国内通用机电暂态仿真软件bpa于psasp中，同步发电机数学模型均为6阶模型，而电磁暂态仿真软件pscad（power systems computer aided design电力系统计算机辅助设计）中同步发电机数学模型则为8阶数学模型。目前的研究几乎一致认为二者的区别对转换结果的影响几乎是可以忽略的，但是电磁暂态仿真中，多机电力系统若考虑励磁系统、原动机及其调速器，系统将面临“维数灾”，给分析计算带来巨大的困难。因此，根据实际需求选择不同阶数的同步发电机模型，从而简化仿真模型，对大规模电力系统电磁暂态仿真的重要性不言而喻。
4.由于同步发电机组的惯性时间常数较大，其机电过程较长，若多台发电机在系统仿真之初直接按其实际特性投入运行，系统将难以稳定，甚至无法达到稳态。对此，针对自定义建模得到的同步发电机组提出一套初始化方法，在仿真开始时，便令该系统到达稳态也是及其必要的。


技术实现要素：

5.针对上述问题，本发明的目的在于提供一种基于pscad的同步发电机组自定义建模及初始化方法，可根据实际仿真需求，选择不同阶数的同步发电机模型，且其初始化功能保证了同步发电机组在仿真开始时便直接到达稳态，跳过了长达数十秒的启动时间，且避免了系统发散的情况。技术方案如下：一种基于pscad的同步发电机组自定义建模及初始化方法，包括以下步骤：步骤1：将同步发电机实用模型转化为状态空间表达式的形式；步骤2：根据所述状态空间表达式画出同步发电机的系统结构图；步骤3：根据所述同步发电机的系统结构图在pscad中进行自定义建模；且根据所述状态空间表达式的输入输出变量不同，分别等效为受控电压源串联定子电阻、受控电流源并联定子电阻两种形式；步骤4：对同步发电机的控制系统，即励磁系统、原动机及其调速器分别进行自定义建模；步骤5：对步骤3与步骤4所搭建的同步发电机组模型进行初始化。
6.进一步的，所述同步发电机组包括五阶实用模型同步发电机、ieee dc1a型励磁系统、以及理想水轮机与pid型调速器；五阶实用模型方程如下：式中：ra为定子电阻；ud、uq分别为发电机定子d、q轴端电压；id、iq分别为发电机定子d、q轴端电流；ef为定子励磁电动势；分别为发电机q轴暂态电动势、次暂态电动势；为d轴次暂态电动势；分别表示d轴同步电抗、暂态电抗、次暂态电抗；分别为q轴同步电抗、次暂态电抗；分别为d轴开路暂态、次暂态时间常数；为q轴开路次暂态时间常数；tj为电机惯性时间常数；tm为机械转矩；ω为转子角速度；δ为q轴领先同步旋转坐标系实轴的角度；p为微分算子；根据输入、输出变量可将上式中前五个方程变换为状态空间表达式，以为状态变量，为输入，为输出，具体状态空间表达式如下：式中，
；；；其中，。
7.更进一步的，所述步骤5具体包括：步骤51：将同步发电机组分为机电暂态部分与电磁暂态部分；所述机电暂态部分包括发电机转子运动方程、原动机及其调速器；以保证转子运动的角速度稳定在参考值附近；所述电磁暂态部分包括发电机电磁方程与励磁系统，以保证机端电压稳定在参考值附近；步骤52：对机电暂态部分和电磁暂态部分进行初始化：步骤521：确定励磁系统结构图中所需初始化的传递函数框图；步骤522：令积分环节输入为0，微分环节输出为0，与稳定状态保持相同；步骤523：令步骤521中所述各传递函数框图采用其常数增益代替，即令复频率s=0；步骤524：确定励磁系统中需要初始化的各环节输出变量，并用下标“0”表示各变量稳态值；步骤525：根据步骤522与步骤523，进一步确定步骤524中各稳态值之间的关系；
步骤526：确定输入变量参考值以及输出变量的稳态值；在五阶实用模型中，需确定ud、uq、id、iq四个变量中三个变量的数值；在已知机端电压的条件下，根据dq变换确定ud、uq；采用定子电阻上的电流，近似等于受控电流源的电流，再进行dq变换以确定id、iq；至此确定各个模块的稳态输出值；步骤527：以步骤526中所确定各数值为输入，步骤524中各稳态值为输出，将步骤525中所确定的各稳态值之间的关系通过系统结构图表示出来，并在pscad中搭建相应的框图；步骤528：通过pscad的复位功能强制各环节的输出保持在步骤527中由框图所确定的稳定值。
8.更进一步的，在仿真开始时，令各传递函数框图输出等于其稳态值，从而使同步发电机组跳过数十秒的启动时间直接到达稳态；计算得到的各环节稳态值通过采样保持器，再赋值给各传递函数框图，以便于通过对采样保持器进行时间设定，保证完成初始化后，同步发电机组真实动态特性的投入运行。
9.与现有技术相比，本发明的有益效果是：1、本发明提供了一套基于同步发电机数学模型及其控制系统结构图，进行同步发电机组自定义建模的方法，可通过选择低阶同步发电机数学模型，以简化电力系统，进而解决多机电力系统电磁暂态仿真时所面临的“维数灾”问题。
10.常用仿真软件均自带同步发电机模型供用户仿真建模使用，若进行大规模电磁暂态仿真，将面临“维数灾”，因此，大多直接将同步发电机简化为理想电压源（同样为仿真软件自带），但是理想电压源并不能反映同步发电机的动态过程。因此，一定条件下，低阶同步发电机模型便是更好的选择，但是部分仿真软件自带同步发电机并不能降低阶数，因此本发明给出了新的建模方法。
11.2、本发明提出了一套适用于自定义同步发电机组的初始化方法，避免了多台发电机直接按其实际特性投入运行系统难以稳定甚至无法达到稳定的情况。
附图说明
12.图1为五阶实用模型输出电流形式系统结构图。
13.图2为同步发电机等效结构图。
14.图3为ieee dc1a型励磁系统的系统结构图。
15.图4为理想水轮机与pid型调速器系统结构图。
16.图5为电力系统组成结构图。
17.图6为励磁系统参考值-5%阶跃时机端电压仿真波形。
18.图7为调速器参考值-2%阶跃时角速度仿真波形。
19.图8为突加20%负载状态下机端电压仿真波形。
20.图9为突加20%负载状态下角速度仿真波形。
21.图10为初始化前后机端电压仿真波形对比。
22.图11为初始化前后角速度仿真波形对比。
23.图12本发明基于pscad的同步发电机组自定义建模及初始化方法的流程图。
具体实施方式
24.下面结合附图和具体实施例对本发明做进一步详细说明。
25.实施例：由五阶实用模型同步发电机、ieee dc1a型励磁系统、以及理想水轮机与pid型调速器组成的同步发电机组。
26.步骤1、将同步发电机实用模型转化为状态空间表达式的形式。
27.常用的五阶实用模型方程如下：式中：ra为定子电阻；ud、uq分别为发电机定子d、q轴端电压；id、iq分别为发电机定子d、q轴端电流；ef为定子励磁电动势；分别为发电机q轴暂态电动势、次暂态电动势；为d轴次暂态电动势；分别表示d轴同步电抗、暂态电抗、次暂态电抗；分别为q轴同步电抗、次暂态电抗；分别为d轴开路暂态、次暂态时间常数；为q轴开路次暂态时间常数；tj为电机惯性时间常数；tm为机械转矩；ω为转子角速度；δ为q轴领先同步旋转坐标系实轴的角度。
28.根据输入、输出变量的不同，可将上式中前五个方程变换为两种状态空间表达式。第一种以为状态变量，为输入，为输出；第二种以为状态变量，为输入，为输出。
29.本实施方案仅对第一种形式进行说明，具体状态空间表达式如下：式中，
；；；；其中，。
30.步骤2、根据步骤1得到的状态空间表达式画出同步发电机的系统结构图。
31.由于两种状态空间表达式均含多个输出变量与输入变量，难以对各输出变量与输入变量之间的传递函数进行解耦，因此，这里采用电气模型的系统结构图进行自定义建模。步骤1所确定的状态空间表达式所对应的系统结构图（如图1所示）。
32.输入变量ef由励磁系统提供，ud、uq可由机端电压测量值进行dq变换得到，输出变量id、iq进行dq反变换得到受控电流源电流ia、ib、ic。故可用受控电流源并联定子电阻进行等效，具体等效发电机模型（如图2所示），图中ra、rb、rc分别为a相、b相、c相定子电阻。
33.步骤3、根据步骤2所确定的同步发电机系统结构图在pscad中进行自定义建模。
34.步骤4、对同步发电机组的控制系统，即励磁系统、原动机及其调速器进行自定义建模。
35.ieee dc1a型励磁系统的系统结构图（如图3所示），图中，vc为发电机机端电压，v
ref
为机端电压参考值；tc、tb分别为超前滞后环节的超前时间常数与滞后时间常数；ka、ta分别为放大器增益与时间常数；te、ke、se分别为励磁机的时间常数、自励系数与饱和系数；
kf、tf、vf分别为励磁系统稳定器的增益、时间常数与输出电压；e
fd
为励磁电压。
36.理想水轮机与pid型调速器组成的系统结构图（如图4所示），图中，ω为转子角速度，ω
ref
为角速度参考值；k
p
、ki、kd分别为调速系统比例环节系数、积分环节系数与微分环节系数；ty为随动系统时间常数；tw为水启动时间；tm、te分别为机械转矩与电磁转矩；b
p
为调差系数。
37.步骤5、对步骤3与步骤4所搭建的同步发电机组模型进行初始化。
38.首先，应将同步发电机组分为机电暂态部分与电磁暂态部分（如图5所示），图中，、分别表示电力系统联络线上的电流、电压；分别表示负荷电流、负荷电压。机电暂态部分由发电机转子运动方程、原动机及其调速器构成，保证转子运动的角速度稳定在参考值附近；电磁暂态部分由发电机电磁方程与励磁系统组成，保证机端电压稳定在参考值附近。之后，分别对机电暂态部分、电磁暂态部分进行初始化。
39.现以电磁暂态部分为例以对初始化方法进行说明，具体过程如下：（1）、确定系统结构图中所需初始化的传递函数框图。根据图3所示励磁系统，需要初始化的模块包括超前滞后环节，输出量为y3；放大器环节，输出量为y2；积分环节，输出量为e
fd
；励磁系统稳定器，输出量为vf。
40.（2）、令积分环节输入为0，微分环节输出为0，与稳定状态保持相同。
41.（3）、令（1）中所述各传递函数框图采用其常数增益代替，即令s=0。
42.（4）、确定（1）中所述各环节输出变量，并用下标“0”表示各变量稳态值。
43.（5）、根据（2）与（3），可进一步确定（4）中各稳态值之间的关系，如下所示：。
44.（6）、确定输入变量参考值以及输出变量的稳态值。输入变量参考值参考实际仿真需求，现分析输出变量的稳态值。由于定子励磁电动势ef为同步发电机状态空间表达式的输入变量，且与励磁电压e
fd
存在如下关系：；式中，x
ad
为d轴电枢反应电抗；rf为励磁电阻。
45.因此，若要确定e
fd0
需从同步发电机数学模型进行分析。分析图1所示运算结图，仅有3个积分环节需进行初始化。而五阶实用模型中，含有8个变量，除去3个状态变量，以及待确定的定子励磁电动势，应至少确定ud、uq、id、iq四个变量中三个变量的数值。在已知机端电压的条件下，可根据dq变换确定ud、uq。又由于定子阻值一般很小，因此可采用定子电阻上的电流，近似等于受控电流源的电流，再进行dq变换以确定id、iq。至此，便可确定各个模块的稳态输出值。
46.（7）、以（6）中所确定各数值为输入，（4）中各稳态值为输出，将（5）中所确定的各稳态值之间的关系通过系统结构图表示出来，并在pscad中搭建相应的框图。
47.（8）、通过pscad的复位功能强制各环节的输出保持在（7）中由框图所确定的稳定值。需要注意的是，计算得到的各环节稳态值应通过采样保持器，再赋值给各传递函数框图，以便于通过对采样保持器进行时间设定，保证完成初始化后，同步发电机组真实动态特性的投入运行。
48.至此，便完成了同步发电机电磁暂态部分的初始化，机电暂态部分初始化方法同电磁暂态部分。
49.现针对本实施例对本发明所提技术路线进行验证：首先，采用pscad库文件中的同步发电机，控制系统则采用自定义建模的dc1a型调速系统、理想水轮机与pid型调速系统，得到模型1。在该模型额定负载稳定运行的状态下，分别对机端电压参考值作-5%阶跃，对角速度参考值作-2%阶跃。由图6和图7可知，机端电压超调量约为3.16%，角速度超调量约为0.45%，其中，pu表示标幺值。
50.其次，将仿真模型1中的同步发电机，更改为自定义建模元件，得到全自定义发电机组仿真模型2。在稳定运行状态下，突加20%负载。由图8和图9可知，机端电压几乎无影响，角速度最大超调量约为2.5%，约6秒后，振荡范围小于
±
2%。
51.最后，对全自定义发电机组仿真模型2进行初始化。由图10和图11可知，通过初始化后，机端电压可提前约9秒进入稳态，角速度则可提前约65秒进入稳态，省去大量启动时间。
52.图12为本发明基于pscad的同步发电机组自定义建模及初始化方法的流程图。
53.综上，本发明解决了多机电力系统电磁暂态仿真时，pscad软件自带同步发电机数学模型阶数过高所面临的“维数灾”问题。由于同步发电机数学模型根据精度不同，具有不同的阶数（也就是维数），常用的同步发电机数学模型为2~6阶，本发明提供了一种实现阶数降低的方法，数学模型阶数的降低，可减少所需的存储空间，加快仿真软件的计算速度。
54.本发明还提出了一套适用于同步发电机组模型的初始化方法，避免了多台发电机直接按其实际特性投入运行系统难以稳定甚至无法达到稳定的情况。
55.最后说明的是，以上实施例仅说明发明的技术方案而非限制。若本领域技术人员在不脱离本发明的宗旨与范围所作的非实质性改变或改进，其均应涵盖在本发明的权利要求范围中。