适用于电磁暂态仿真的电动机负荷暂态过程计算方法技术

本发明专利技术涉及电网暂态仿真计算的技术领域，更具体地，涉及适用于电磁暂态仿真的电动机负荷暂态过程计算方法，快速计算电力系统对称和不对称故障电压跌落下感应电动机负荷暂态特性：(1)先求解转子的一阶机械暂态方程，得到扰动后转子的转速和滑差变化；(2)求解转子相量形式的暂态电势微分方程，得到暂态电势的变化；(3)获得扰动后电动机的电流和功率响应。本发明专利技术与PSCAD/EMTDC电磁暂态仿真结果的比较，验证了算法的有效性。与PSD/BPA机电暂态仿真结果的比较，体现了算法在不对称电压跌落下的计算优势。解析法显式给出了受扰后感应电动机各机械、电气参量的计算表达式，除可用于快速评估感应电动机与电压跌落的相互影响、还可应用于电力系统稳定性分析中。

Calculation method of motor load transient process suitable for electromagnetic transient simulation

The invention relates to the technical field of power grid transient simulation calculation, more specifically, to the calculation method of motor load transient process suitable for electromagnetic transient simulation, which can quickly calculate the load transient characteristics of induction motor under symmetrical and asymmetric fault voltage drop of power system: (1) first solve the first-order mechanical transient equation of rotor to obtain the change of rotor speed and slip after disturbance; (2) The transient electromotive force differential equation in the form of rotor phasor is solved to obtain the change of transient electromotive force; (3) the current and power response of the motor after disturbance are obtained. The comparison between the invention and PSCAD / EMTDC electromagnetic transient simulation results verifies the effectiveness of the algorithm. Compared with the results of PSD / BPA electromechanical transient simulation, the algorithm has the advantage of calculation under asymmetric voltage drop. The analytical method can be used not only to evaluate the interaction between induction motor and voltage drop, but also to analyze the stability of power system.

【技术实现步骤摘要】
适用于电磁暂态仿真的电动机负荷暂态过程计算方法
本专利技术涉及电网暂态仿真计算的
，更具体地，涉及适用于电磁暂态仿真的电动机负荷暂态过程计算方法。
技术介绍
随着“西电东送”战略的实施和深化，我国在“长三角”和“珠三角”等负荷中心均已形成典型的受端系统结构。受端系统负荷重、感应电动机负荷比例高，电网发生故障及恢复过程中，感应电动机由于电磁转矩下降，滑差增大，常会吸收大量的暂态功率，导致系统暂态电压难以恢复。而各项研究均表明电动机暂态特性和电压跌落相互影响的重要性，并先后提出有试验法、时域仿真法与解析法等三类研究方法。试验法借助电压跌落发生器产生各种类别的电压跌落波形，记录电动机的输出响应，进而分析电压跌落与电动机的相互影响。时域仿真法借助电力系统电磁暂态或机电暂态仿真程序，考虑较为精确的电动机暂态模型，通过数值计算对问题进行探讨。二者的优点在于结果真实可信，但为揭示某一因素的影响，需要多次试验或仿真，除繁琐、费时外，在分析和解释试验和仿真现象时，都显得力所不及。而解析法需要借助电路和电机分析理论，显式求解电动机暂态响应和系统电压间的交互作用，是最根本、最透彻的研究方法，但难度较大，且不能准确求得受扰后电动机的全部运行和状态变量；仅能考虑恒转矩机械负载，也会造成电压稳定性分析结果偏悲观。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供一种适用于电磁暂态仿真的电动机负荷暂态过程计算方法，在将外部网络等效为简单两节点系统的基础上，提出不对称故障电压跌落下精确计算感应电动机暂态响应的解析计算方法，具有较好的有效性和准确性。为解决上述技术问题，本专利技术采用的技术方案是：提供适用于电磁暂态仿真的电动机负荷暂态过程计算方法，电动机为感应电动机并设于供电系统中；其特征在于，所述计算方法包括以下步骤：S10.建立感应电动机转子的一阶机械暂态模型并求解，得到电源发生对称或不对称故障电压跌落时感应电动机的转速ωm和滑差s变化；感应电动机转子的一阶机械暂态模型表示为下式：式(1)中，ωm为感应电动机的转速，Te和Tm分别是感应电动机的电磁转矩和机械负载转矩，Tj为电动机惯性时间常数；式(2)中，ωs为感应电动机的同步转速，s为滑差；S20.根据步骤S10求得的感应电动机的转速ωm和滑差s变化，求解感应电动机转子向量形式的暂态电势微分方程得到暂态电势的变化；S30.根据步骤S20中求得的暂态电势的变化获得扰动后感应电动机的电流相应和功率响应。本专利技术的适用于电磁暂态仿真的电动机负荷暂态过程计算方法，优选地，步骤S10中，感应电动机转子的一阶机械暂态模型为忽略定子绕组和转子绕组的暂态等值电路模型，所述暂态等值电路包括串联连接的定子电阻Rs、定子电抗Xs以及激磁电抗Xm，所述激磁电抗Xm的两端并联连接有单笼模型转子电抗Xr以及单笼模型转子电阻Rr。优选地，步骤S10中：不对称故障电压跌落期间，感应电动机转子的运动方程表示为：式(3)中，K1、K2为常数；Vs1、Vs2分别表示外施定子电压的正序分量、负序分量；电压跌落清除后，转子的运动方程不包含负序电磁转矩Te2，感应电动机转子运动方程表示为：优选地，感应电动机的初始转速为ωn，从电压跌落发生时刻t0开始不断下降，到故障清除时刻t1，感应电动机转速降至ω'm；在电压跌落清除后，感应电动机转速不断上升恢复至感应电动机初始转速ωn。优选地，根据式(3)求解故障期间感应电动机转速ωm-dur：式中，t为感应电动机工作时刻；电压跌落清除后，根据式(4)求解感应电动机转速度ωm-after：优选地，步骤S20按以下步骤进行：S21.当电动机定子端施加正序电压时，以向量形式表示外施定子电压如下式：式(7)中，和分别为转子正序暂态电势和定子正序电流相量，X为转子开路电抗，X'表示转子短路电抗，表示定子开路、转子回路暂态时间常数；以向量形式表示外施定子电流如下式：式(9)中，Rs为定子电阻；将式(8)代入式(7)中，得到下式：S22.当电动机定子端施加负序电压时，流过负序电流建立负向旋转磁场，所述负向旋转磁场与转子旋转方向相反，负序转差率为(2-s)，转子负序暂态电势微分方程的相量形式可表示为：式(10)中，为转子负序暂态电势；S23.整理式(9)和(10)后，得：令：S24.将步骤S10计算得到的滑差s代入，则KE1、KE2和BE均为已知量，可将式(12)简化为：S25.求解式(13)，不对称故障期间电动机转子暂态电势的解析表达式为：式(14)中，和为感应电动机正、负序暂态电势初值，计算表达式为：式(15)中Vs10表示外施定子电压正序分量的初始值，Is10表示外施定子电流正序分量的初始值；S26.电压跌落清除后恢复过程中，定子电压中不含负序电压分量，电动机的暂态电势仅含正序分量，解析表达式为：式(16)中，为电压跌落清除时刻电动机的正序暂态电势。优选地，步骤S30按以下步骤进行：S31.将步骤S25中计算得到的故障期间电动机转子暂态电势和的d轴分量E'd1-dur、E'd2-dur和q轴分量E'q1-dur、E'q2-dur分别代入式(17)，得电压跌落期间电动机d、q轴的正序定子电流Ids1、Iqs1和负序定子电流Ids2、Iqs2：S32.电压跌落期间，感应电动机消耗的有功Pd和无功Qd计算如式(18)：式中，Pd1、Pd2为正、负序有功功率，Qd1、Qd2为正、负序无功功率，负序电流产生的功率与正序电流产生的功率方向相反；S33.电压跌落清除后，暂态电势的d轴分量E'd1-after和q轴分量E'q1-after分别代入式(17)，可得电压跌落期间电动机d、q轴的正序定子电流Ids1、Iqs1，电动机消耗的有功Pd和无功Qd计算如式(19)：求得电动机的功率消耗(Pd+jQd)后，定子电流正序分量负序分量计算如式(20)：与现有技术相比，本专利技术的有益效果是：本专利技术的计算方法，计算模型均为代数表达式，与电力系统机电和电磁暂态仿真联立求解微分-代数方程组相比，在计算速度上具有明显优势；本专利技术的计算方法，计入了不对称故障时负序分量对电动机暂态响应的影响，计算精度高于只考虑正序基波分量的机电暂态仿真；本专利技术的计算方法，显式计算电动机机械和电气参量，易于分析和解释结果与现象。附图说明图1为适用于电磁暂态仿真的电动机负荷暂态过程计算方法的流程图；图2为实施例一的暂态等值电路的示意图；图3为实施例一的单机带动感应电动机的供电系统示意图；图4为电压跌落时感应电动机转速本文档来自技高网...

【技术保护点】
1.适用于电磁暂态仿真的电动机负荷暂态过程计算方法，所述电动机为感应电动机并设于供电系统中；其特征在于，所述计算方法包括以下步骤：/nS10.建立感应电动机转子的一阶机械暂态模型并求解，得到电源发生对称或不对称故障电压跌落时感应电动机的转速ω

【技术特征摘要】
1.适用于电磁暂态仿真的电动机负荷暂态过程计算方法，所述电动机为感应电动机并设于供电系统中；其特征在于，所述计算方法包括以下步骤：
S10.建立感应电动机转子的一阶机械暂态模型并求解，得到电源发生对称或不对称故障电压跌落时感应电动机的转速ωm和滑差s变化；感应电动机转子的一阶机械暂态模型表示为下式：



式(1)中，ωm为感应电动机的转速，Te和Tm分别是感应电动机的电磁转矩和机械负载转矩，Tj为电动机惯性时间常数；



式(2)中，ωs为感应电动机的同步转速，s为滑差；
S20.根据步骤S10求得的感应电动机的转速ωm和滑差s变化，求解感应电动机转子向量形式的暂态电势微分方程得到暂态电势的变化；
S30.根据步骤S20中求得的暂态电势的变化获得扰动后感应电动机的电流相应和功率响应。


2.根据权利要求1所述的适用于电磁暂态仿真的电动机负荷暂态过程计算方法，其特征在于，步骤S10中，感应电动机转子的一阶机械暂态模型为忽略定子绕组和转子绕组的暂态等值电路模型，所述暂态等值电路包括串联连接的定子电阻Rs、定子电抗Xs以及激磁电抗Xm，所述激磁电抗Xm的两端并联连接有单笼模型转子电抗Xr以及单笼模型转子电阻Rr。


3.根据权利要求2所述的适用于电磁暂态仿真的电动机负荷暂态过程计算方法，其特征在于，步骤S10中：
不对称故障电压跌落期间，感应电动机转子的运动方程表示为：



式(3)中，K1、K2为常数；Vs1、Vs2分别表示外施定子电压的正序分量、负序分量；
电压跌落清除后，转子的运动方程不包含负序电磁转矩Te2，感应电动机转子运动方程表示为：





4.根据权利要求3所述的适用于电磁暂态仿真的电动机负荷暂态过程计算方法，其特征在于，感应电动机的初始转速为ωn，从电压跌落发生时刻t0开始不断下降，到故障清除时刻t1，感应电动机转速降至ω'm；在电压跌落清除后，感应电动机转速不断上升恢复至感应电动机初始转速ωn。


5.根据权利要求4所述的适用于电磁暂态仿真的电动机负荷暂态过程计算方法，其特征在于，根据式(3)求解故障期间感应电动机转速ωm-dur：



式中，t为感应电动机工作时刻；
电压跌落清除后，根据式(4)求解感应电动机转速度ωm-after：





6.根据权利要求1至5任一项所述的适用于电磁暂态仿真的电动机负荷暂态过程计算方法，其特征在于，步骤S20按以下步骤进行：
S21.当电动机定子端施加正序电压时，以向量形式表示外...

【专利技术属性】
技术研发人员：廖卫平，伍建炜，张艳，温健锋，黄练栋，麦炳灿，颜金佑，刘海光，周锐，谢锟，
申请(专利权)人：广东电网有限责任公司，广东电网有限责任公司江门供电局，
类型：发明
国别省市：广东;44