用于异步电动机负荷的通用电磁暂态仿真解耦方法和系统技术方案

本发明专利技术公开了一种用于异步电动机负荷的通用电磁暂态仿真解耦方法和系统，包括：确定异步电动机解耦模型的单相简化等值电路；通过所述单相简化等值电路，获得异步电动机负荷

【技术实现步骤摘要】
用于异步电动机负荷的通用电磁暂态仿真解耦方法和系统


[0001]本专利技术涉及电磁暂态仿真领域，具体是涉及一种用于异步电动机负荷的通用电磁暂态仿真解耦方法和系统。

技术介绍

[0002]负荷模型对研究大规模电力系统的暂态稳定性和电压稳定性起着相当重要的作用。从系统角度看，电动机负荷占电力系统总负荷比例超过50％，因此在电磁暂态仿真中通常用电动机负荷模型来表示动态负荷模型。当系统中包含大量电动机负荷时，其动态特性会显著影响仿真结果。为提高仿真效率，在电力系统仿真中，异步电动机通常使用集总参数模型(忽略了电机的几何复杂性、涡流、饱和以及磁滞效应)，包括相域模型(phase domain，PD)，VBR(voltage
‑
behind
‑
reactance)模型和d
‑
q模型。
[0003](1)在PD模型中，电机变量均在abc坐标系下表示，定子和转子电路与外部网络直接相连，电机变量会和外部网络同时求解，具有较好的数值稳定性，但由于PD模型的电感矩阵参数会随着转子位置的变化而改变，因此PD模型的计算复杂度较高。
[0004](2)在VBR模型中，定子变量用abc坐标表示，而转子变量用d
‑
q坐标表示，电机等效为一个与RL支路串联的电压源。VBR模型的定子可以直接与外部网络相连，电机的状态变量和网络方程可以同时求解，改善了数值稳定性，然而这是以时变的RL参数为代价的。
[0005](3)将Park变换应用于PD模型，得到d
>‑
q模型，其优点是在d
‑
q坐标系下可以得到恒定的电感矩阵，仿真效率较高。异步电动机模型中应用最广泛的d
‑
q模型，d
‑
q模型在PSCAD/EMTDC，EMTP
‑
RV，Matlab/Simulink等商业电磁暂态仿真软件中广泛采用。
[0006]与本专利技术相关的现有技术一
[0007]诺顿电流源法，该方法是商业软件PSCAD/EMTDC的采用电机
‑
电网接口方法。
[0008]诺顿电流源法将异步电动机的模型等效为诺顿电流源与电网相连。这种类型的电动机
‑
电网接口如图2所示。注入电流I
m
(t)是利用上一个时间步长的异步电动机定子端电压来计算的，这意味着在计算异步电动机变量时有一个时间步长的延迟。因此可能导致数值不稳定的问题，特别是当异步电动机处在接近开路的状态。为了提高数值稳定性，防止异步电动机完全开路，在电动机
‑
网络接口中引入了一个端口特征阻抗r
″
。然后，这个增加的电阻的影响通过注入到异步电动机定子端口的额外电流I
c
(t)源来补偿。
[0009]技术一缺点：由于注入电流I
m
(t)是根据上一个时步的定子端电压计算得到，因此注入电流I
m
(t)的计算存在一个步长的延时。数学上，这种延时的引入相当于数值积分中的显式方法，这一因素可能会影响仿真的数值稳定性。为了保证接口的数值稳定性和精度，该方法需要对注入电流I
m
(t)的求解进行多次迭代计算，即需要多次求解电动机
‑
电网整体方程。当系统规模很大时，电动机
‑
电网整体方程阶数很高，多次求解电动机
‑
电网整体方程的效率将十分低下。
[0010]与本专利技术相关的现有技术二
[0011]网络迭代法，该方法是商业软件EMTP
‑
RV的采用电机
‑
电网接口方法。
[0012]在EMTP
‑
RV中，异步电动机模型与网络方程采用牛顿迭代法联立求解。此外，在EMTP
‑
RV异步电动机模型中，利用速度和电压循环迭代获得每个时间步长的电动机方程的收敛解。这两个循环的收敛判据是由用户指定的误差偏差来控制。
[0013]技术二的缺点：异步电动机和网络方程在每个时间步长的联立迭代求解通常需要更多的计算资源，当网络方程规模增加时，整体的仿真效率将显著降低。

技术实现思路

[0014]为解决含大量异步电动机负荷的电力系统的电磁暂态仿真的全网方程求解迭代次数过多，导致仿真效率低，接口数值稳定性差的问题，本专利技术提供一种用于异步电动机负荷的通用电磁暂态仿真解耦方法和系统，所述方法，包括：
[0015]确定异步电动机解耦模型的单相等值电路；
[0016]通过合并所述单相等值电路中的异步电动机阻抗，以及使用集中参数表示分布式参数解耦元件，简化所述单相等值电路，获得异步电动机解耦模型的单相简化等值电路；
[0017]通过所述单相简化等值电路，获得异步电动机负荷
‑
电网电磁暂态仿真接口；
[0018]基于所述异步电动机负荷
‑
电网电磁暂态仿真接口和d
‑
q模型，对异步电动动机负荷电磁暂态仿真，在当前仿真时刻，根据系统节点电压更新解耦元件历史电流源；将定子三相电压和转子三相电压经Park变换，获得d
‑
q参考系下的定子电压和转子电压；根据所述定子电压和转子电压，求解异步电动机状态方程，获得d
‑
q参考系下的定子电流和转子电流；
[0019]通过将所述定子电流通过Park反变换获得定子三相电流，将所述定子三相电流传递至电流源，求解异步电动机定子侧接口电路，获得异步电动机定子侧电压；
[0020]在所述异步电动机状态方程满足误差条件时，基于解耦元件特征电导，更新系统的节点电导矩阵，基于解耦元件定子侧历史电流源和补偿电容历史电流源，更新系统的节点电压方程右端项电流；
[0021]根据所述节点电导矩阵和节点电压方程右端项电流，获得系统的节点电压；通过所述节点电压，更新系统补偿电容历史电流源；
[0022]根据预设步长，更新当前仿真时刻，再次对异步电动动机负荷电磁暂态仿真，在当前仿真时刻大于仿真时间时，结束仿真。
[0023]进一步的，确定异步电动机解耦模型的单相等值电路，包括：
[0024]确定用于构建异步电动机解耦模型的输入数据和相关基准值；
[0025]通过所述输入数据和相关基准值，以及异步电动机原始模型的单相电路拓扑，确定异步电动机解耦模型的单相等值电路。
[0026]进一步的，通过合并异步电动机阻抗，以及使用集中参数表示分布式参数解耦元件，简化所述单相等值电路，获得异步电动机解耦模型的单相简化等值电路，包括：
[0027]通过以下公式合并异步电动机阻抗：
[0028][0029]其中，Z
im
为异步电动机等值阻抗，z
s
为定子等值阻抗，z
r
为转子等值阻抗，z
m
为励磁阻抗；
[0030]将分布式参数解耦元件用π型本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种用于异步电动机负荷的通用电磁暂态仿真解耦方法，其特征在于，包括：确定异步电动机解耦模型的单相等值电路；通过合并所述单相等值电路中的异步电动机阻抗，以及使用集中参数表示分布式参数解耦元件，简化所述单相等值电路，获得异步电动机解耦模型的单相简化等值电路；通过所述单相简化等值电路，获得异步电动机负荷
‑
电网电磁暂态仿真接口；基于所述异步电动机负荷
‑
电网电磁暂态仿真接口和d
‑
q模型，对异步电动动机负荷电磁暂态仿真，在当前仿真时刻，根据系统节点电压更新解耦元件历史电流源；将定子三相电压和转子三相电压经Park变换，获得d
‑
q参考系下的定子电压和转子电压；根据所述定子电压和转子电压，求解异步电动机状态方程，获得d
‑
q参考系下的定子电流和转子电流；通过将所述定子电流通过Park反变换获得定子三相电流，将所述定子三相电流传递至电流源，求解异步电动机定子侧接口电路，获得异步电动机定子侧电压；在所述异步电动机状态方程满足误差条件时，基于解耦元件特征电导，更新系统的节点电导矩阵，基于解耦元件定子侧历史电流源和补偿电容历史电流源，更新系统的节点电压方程右端项电流；根据所述节点电导矩阵和节点电压方程右端项电流，获得系统的节点电压；通过所述节点电压，更新系统补偿电容历史电流源；根据预设步长，更新当前仿真时刻，再次对异步电动动机负荷电磁暂态仿真，在当前仿真时刻大于仿真时间时，结束仿真。2.根据权利要求1所述的方法，其特征在于，确定异步电动机解耦模型的单相等值电路，包括：确定用于构建异步电动机解耦模型的输入数据和相关基准值；通过所述输入数据和相关基准值，以及异步电动机原始模型的单相电路拓扑，确定异步电动机解耦模型的单相等值电路。3.根据权利要求1所述的方法，其特征在于，通过合并异步电动机阻抗，以及使用集中参数表示分布式参数解耦元件，简化所述单相等值电路，获得异步电动机解耦模型的单相简化等值电路，包括：通过以下公式合并异步电动机阻抗：其中，Z
im
为异步电动机等值阻抗，z
s
为定子等值阻抗，z
r
为转子等值阻抗，z
m
为励磁阻抗；将分布式参数解耦元件用π型集中参数表示，其中解耦元件等值阻抗X
d
和解耦元件等值电纳B
d
的计算公式为：
其中，z
c
为解耦元件特征阻抗，γ为传输线的传播常数，j为虚数单位，ω为同步角速度，L
d
为解耦元件电感，C
d
为解耦元件电容；通过上述方式获得异步电动机解耦模型的单相简化等值电路。4.根据权利要求1所述的方法，其特征在于，所述异步电动机负荷
‑
电网电磁暂态仿真接口，具体包括：1个异步电动机方程求解模块和3个电路模块；所述3个电路模块，分别为异步电动机定子侧接口电路，异步电动机电网侧接口电路和外部电网；异步电动机电网侧接口电路与外部电网直接相连，相连节点为k；异步电动机定子侧接口电路包含1个节点，该节点为m；将异步电动机电网侧接口电路的解耦元件特征电导G
c
和解耦元件电网侧历史电流源I
k
，与异步电动机定子侧接口电路的G
c
和解耦元件定子侧历史电流源I
m
合称为解耦元件。5.根据权利要求1所述的方法，其特征在于，基于所述异步电动机负荷
‑
电网电磁暂态仿真接口和d
‑
q模型，对异步电动动机负荷电磁暂态仿真，在当前仿真时刻，根据系统节点电压更新解耦元件历史电流源，包括：开始对异步电动动机负荷电磁暂态仿真，令仿真时刻t＝0；将仿真时刻向前推进一个步长，更新仿真时刻：t＝t+Δt，将迭代次数置零，令j＝0；根据t
‑
Δt时刻求解得到的系统节点电压U，更新解耦元件历史电流源，具体为：I
m
(t
‑
Δt)＝I
m
(t
‑
2Δt)+α(
‑
2G
c
u
k
(t
‑
Δt)+I
k
(t
‑
2Δt)
‑
I
m
(t
‑
2Δt))I
k
(t
‑
Δt)＝I
k
(t
‑
2Δt)+α(
‑
2G
c
u
m
(t
‑
Δt)+I
m
(t
‑
2Δt)
‑
I
k
(t
‑
2Δt))其中，u
k
是节点k电压，u
m
是节点m电压，α∈[0,1]是插值系数。6.根据权利要求1所述的方法，其特征在于，将定子三相电压和转子三相电压经Park变换，获得d
‑
q参考系下的定子电压和转子电压，包括：令迭代次数j＝j+1，将上一次迭代获得的定子三相电压v
abcs
(t)和转子三相电压v
abcr
(t)，经Park变...

【专利技术属性】
技术研发人员：刘文焯，庞博涵，连攀杰，郁舒雁，许克，汤涌，
申请(专利权)人：中国电力科学研究院有限公司，
类型：发明
国别省市：