CN107086606A - 一种配电网综合等效异步电动机负荷模型参数辨识方法 - Google Patents
一种配电网综合等效异步电动机负荷模型参数辨识方法 Download PDFInfo
- Publication number
- CN107086606A CN107086606A CN201710443376.3A CN201710443376A CN107086606A CN 107086606 A CN107086606 A CN 107086606A CN 201710443376 A CN201710443376 A CN 201710443376A CN 107086606 A CN107086606 A CN 107086606A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- msup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Abstract
本发明公开了一种配电网综合等效异步电动机负荷模型参数辨识方法,首先考虑配电网参数及无功补偿,建立综合异步电动机负荷模型,然后给定实测数据,计算状态变量初始值,接着建立负荷模型的参数辨识准则函数,并利用CPSO算法对目标函数进行寻优,并输出负荷模型,最后对负荷模型进行校验,并输出结果。本发明考虑了计入配网无功补偿和变压器分接头的影响,使得负荷模型更加准确,CPSO算法中引入了混沌优化搜索的概率,以对最优粒子进行混沌优化搜索,提高了算法收敛精度,有效地提高了模型参数辨识的准确性。
Description
技术领域
本发明涉及电力系统仿真技术领域,尤其涉及一种配电网综合等效异步电动机负荷模型参数辨识方法。
背景技术
随着分布式电源广泛接入中低压配电网,使得原有配电网的组成、结构和潮流流向有了较大的改变。在电力系统综合负荷建模研究中,模型参数的获取是最关键的问题之一。美国 WSCC及IEEE等机构根据负荷类型的不同推荐了相应的典型异步电动机模型参数,在我国的工程仿真计算中,曾长期取“典型参数”,但是随着电网规模的快速发展,综合负荷构成的复杂程度日益增加,“典型参数”对实际电网的适应性问题近年来得到广泛关注,且开展了大量的理论与应用研究。总体测辨法将负荷群体看成一个整体,先进行现场采集数据,确定模型结构,再根据现场采集数据利用系统辨识理论辨识出模型参数,该方法具有简单实用,建模数据直接来源于实际系统等优点。当分布式电源接入容量不大时,原有的经典负荷模型 (CLM)和考虑配电网支路的综合负荷模型(SLM)仍具有较好的描述能力,但随着分布式电源接入容量增大,原有模型已不具备准确的描述能力。
发明内容
本发明的目的在于提供一种配电网综合等效异步电动机负荷模型参数辨识方法,能够提高模型参数辨识的准确性。
为实现上述目的,本发明采用如下技术方案:
一种配电网综合等效异步电动机负荷模型参数辨识方法,依次包括以下步骤:
(1)考虑配电网参数及无功补偿,建立综合异步电动机负荷模型;
(2)给定实测数据,计算状态变量初始值;
(3)建立负荷模型的参数辨识准则函数,利用CPSO算法对目标函数进行寻优,并输出负荷模型;
(4)对负荷模型进行校验,并输出结果。
优选地,所述步骤(1)具体包括以下步骤:
①在110kV线路上增加一个变压器,在10kV母线上并联异步电动机,同时考虑变压器分接头引入基准变换以进行配电网络与异步电动机对接实现,设异步电动机在自身基准下的端电压、暂态内电势、从电网吸收的电流、功率及等值阻抗的物理量分别为: Uim、E′im、Iim、Pim、Qim、Zim,下标im为区分作用,并非变量,各个物理量相对应系统基准下的标么值分别为US、ES'、IS、PD、QD、ZS,下标S和D为区分作用,并非变量,设变压器高压侧实际运行分接头电压和低压侧额定电压分别为和与UBS对应的变压器低压侧系统基准电压为UBL,其中,UBS表示系统电压基准,下标T、TN、BS及BL均为区分作用,并非变量,则有:
K=SBS/SBM (2)
UBL=UBS/kT (3)
其中,kT*表示变压器高压侧与低压侧的非标准电压比的标幺值,kT表示变压器高压侧与低压侧的非标准电压比,K表示系统与异步电动机之基准功率比,kB表示系统电压基准与感应电动机电压基准比值,SBS表示系统功率,SBM表示感应电动机功率,下标B、BS、BM均为区分作用,并非变量;
由式(1)可知,当取时,若变压器运行于主抽头,即或取各基准电压为网络平均额定电压且变压器变比为平均额定电压比,则kT*=1,设异步电动机电压基准为则由电压等级归算及标幺值换算原理,即可推导出“基准变换”的变换矩阵如式(4)所示:
其中,各个物理量的定义如步骤(1)所述;
配电网综合等效后的线路-变压器组首、末端功率及电压平衡关系如式(5)~式(7),变压器低压侧母线L的功率平衡关系如公式(8),
PS+jQS=(PZIP+jQZIP)+(PD+jQD)-jQC (8)
其中,PS表示线路末端有功功率,QS表示线路末端无功功率,P表示线路首端有功功率, Q表示线路首端无功功率,Re表示线路的等值电阻,Xe表示线路的等值电抗,US表示异步电动机的端电压的标幺值,USx表示感应电动机在系统基准下的端电压在x轴的坐标分量,USy表示感应电动机在系统基准下的端电压在y轴的坐标分量,U表示线路首端电压,下标e、Sx 及Sy均为区分作用,并非变量;
②在10kV母线并联静态负荷,将静态负荷用多项式(ZIP)模型描述如式(9):
其中,模型参数满足ZP+IP+PP=1、ZQ+IQ+PQ=1及 (0)为 US(0)的有名值,US(0)表示异步电动机的初始端电压,PZIP(t)表示静态负荷在t时刻的有功功率,QZIP(t)表示静态负荷在t时刻的无功功率,PZIP(0)表示静态负荷的初始有功功率,QZIP(0)表示静态负荷的初始无功功率,ZP表示静态负荷的有功阻抗,ZQ表示静态负荷的无功阻抗,US(t)表示异步电动机在t时刻的端电压,kU0为已知的电压变换系数;
③在静态负荷上并联一附加动态无功补偿元件,方程式如下:
其中,QC(t)表示无功补偿元件在t时刻的无功功率,下标c为区分作用,并非变量,f 表示频率,在标幺值下XC表示电容器容抗,XC0表示角频率近似取1时的值,即正常运行时的电容器容抗,表示补偿电容器的电容,ω表示角频率;
④将异步电动机模型用三阶机电暂态微分方程描述,如式(11):
式中,e'im=e'im.x+je'im.y为异步电动机暂态电势,e′im·x及e′im·y分别表示暂态电势在x轴和 y轴的坐标分量,t表示时间,s表示异步电动机的转差率,w0表示系统的同步角频率,频率为50Hz时,w0=314.16rad/s,Td0'表示异步电动机暂态电势的衰减时间常数,下标d仅为区分作用,并非变量,Xs表示定子绕组漏抗,Rs表示定子绕组电阻,X'表示定子和转子暂态电抗,Uim=Uim.x+jUim.y表示负荷端电压,Uimx及Uimy分别表示感应电动机电压在x轴和y轴的坐标分量,Tj表示转子惯性时间常数,Xr表示转子绕组漏抗,Xm表示激磁电抗,Rr表示转子电阻,Tm表示异步电动机的机械负载功率,具体计算公式如式(12):
Tm=T0[Aωr 2+Bωr+C] (12)
其中,T0为异步电动机负载率;ωr=1-s为转子侧角速度;A、B、C为机械转矩系数,且满足A+B+C=1;
异步电动机输出电流方程为:
其中,Iim.x及Iim.y分别表示异步电动机输出电流在x轴和y轴的坐标分量;
异步电动机的功率为:
优选地,所述步骤(2)具体包括以下步骤:
①给定实测电压激励U(k)、实测有功功率P(k)及无功功率Q(k);
②计算异步电动机及静态负荷的初始有功功率,设异步电动机的动态负荷占空比为
其中,PD(0)表示感应电动机初始有功功率,Ps(0)表示变压器低压侧母线初始有功功率,下标m、D及S均为区分作用,并非变量;
由式(15)及综合异步电动机等值电路可求得异步电动机及静态负荷初始稳态有功功率为
其中,Pim(0)表示异步电动机的初始稳态有功功率,PZIP(0)表示静态负荷的初始稳态有功功率;
③计算异步电动机的初始滑差,由综合异步电动机等值电路可知异步电动机稳态等值阻抗和导纳如式(17)所示:
其中,Zim(0)表示异步电动机的稳态等值阻抗,Yim(0)表示异步电动机的稳压等值导纳,Gim(0)表示异步电动机的稳态等值电导,Bim(0)表示异步电动机的稳态等值电纳, R=Rr/S(0),Rr表示转子电阻,S(0)表示初始转差率;
结合(16),式(17)可写为:
其中,Xsm、Xrs、Xrm、Rm、XP均无实际物理意义,为替代参数,其取值如下:
其中,
从而可得关于R的一元二次方程AR2+BR+C=0 (22)
A1、B1、C1、A2、B2、C2均无实际物理意义,为替代参数,其取值如下:
其中,A、B及C的值由式(23)确定:
解式(22)可得异步电动机初始转差率如下:
④计算异步电动机及静态负荷初始无功功率;
异步电动机的初始转差率确定后,由异步电动机应满足的初始稳态条件并由式(4)可得系统基准下的异步电动机初始稳态无功功率:
再由节点功率平衡原理唯一确定静态负荷的初始稳态无功功率如式(26):
QZIP(0)=QS(0)-QD(0)=QS(0)-Qim(0)/K (26)
⑤计算异步电动机暂态电动势初始值;
由异步电动机暂态等效电路可求得暂态电动势初始值如式(27)所示,
其中初始电流依初始稳态条件按式(28)所示确定:
⑥计算表异步电动机的机械负载率T0:
其中,ωr0=1-s(0)。
优选地,为避免仿真中开方出现复数问题,需要对式(24)进行以下约束:
①A>0,B<0且B2-4AC≥0;
②额定初始滑差率取值范围为:0.015≤s≤0.07。
优选地,所述步骤(3)中,负荷模型的参数辨识准则函数如式(30)所示:
其中,x(t)表示系统状态向量,u(t)为系统输入向量,ym(k)=[Pm(k),Qm(k)]T为输入u(t)时负荷模型得到的输出响应,y(k)表示符合模型实际输出相应,k为采样开始时刻,L为采样数,α=[Re,Xe,ZP,IP,ZQ,IQ,XC,K,Km,KT*,Rs,Xs,Rr,Xr,Rm,Xm,Tj,A,B] 表示独立待辨识参数,β=[YZIP0,PP,PQ,T0,C]表示非独立待辨识参数,其中,PP=1-ZP-IP,PQ=1-PQ-IQ,C=1-A-B,YZIP0表示静态负荷初始导纳值,QZIP0表示静态负荷初始无功功率,US0表示异步电动机的初始端电压,GZIP0表示静态负荷初始电导,BZIP0表示静态负荷初始电纳。
优选地,所述负荷模型的参数辨识准则函数的计算过程如下:
①在实际变电站系统中输入电压激励U(k)、实测有功功率P(k)及无功功率Q(k),并给定独立待辨识参数初值α,由实测数据通过式(5)~(7)得到系统基准下的末端电压激励,通过式(4)变换为自身基准下的末端电压激励;
②由步骤(2)计算出异步电动机及静态负荷初始功率响应,求出状态变量初值和非独立参数,利用改进欧拉法解微分方程组(11),求得系统转差率S(k)及暂态电动势eimx(k),并带入输出方程(13)~(14),求得感应电动机的功率响应,由式(9)求得静态负荷的功率响应,再根据式(4)~(8)求得异步电动机在自身基准下的功率响应及线路末端功率响应,最终得到线路首端功率响应,即式(6);
③采用改进CPSO算法对以下目标函数寻优:
其中,各个参数的含义与权利要求5中相同,不再赘述。
优选地,利用改进CPSO算法进行目标函数寻优的过程如下:
①根据经验值对各个参数初始化,包括学习因子c1和c2,约束因子a,最大进化代数itermax, wmin和wmax,粒子数N,混沌搜索步长调节参数β和混沌搜索步数ckmax,其中,下标max及 min均为区分作用,并非变量;
②随机生成N个粒子的xi和vi,令k=0;
③按式计算wk;
④根据步骤(2)中求得的模型的有功功率计无功功率,计算粒子适应度,即式(30),并设定粒子的Pbest及种群的gbest,Pbest表示种群搜索到的最优解,gbest表示种群中最好的粒子;
⑤按式更新粒子的速度和位置,对种群中的粒子进行操作,并更新每个粒子的pi和种群的g,式中:k表示第k次迭代,c1、c2为学习因子,r1、r2是[0,1]间的随机数,a为控制速度权重的约束因子,wk为惯性权重;如果则取如果则取如果或者则重新初始化 和为vd的取值范围,下标d及id均仅为区分作用,并非变量;
⑥按下式计算混沌优化搜索的概率Pk
如果rand(0,1)≤Pk,则令d=1,对粒子g中变量进行混沌优化搜索,其余的D-1个变量保持不变,之后进入步骤⑦,否则直接进入步骤⑦,其中,rand(0,1)为[0,1]间的随机数;
⑦令k=k+1,此时若k≥L,则按照式(30)将此函数模型响应与实测功率进行计算,然后进入步骤⑧,若k<L,则返回步骤⑤;
⑧收敛判断,如果粒子g的适应值小于给定的阀值或者k>itermax,则进化过程成功结束,然后进入步骤⑨,否则求解第i此最优化问题minJ→αi,并令i=i+1,然后返回步骤①;
⑨令α=α(i-1),β=β(i-1),输出参数α、β及模型响应序列ym(k),其中k=1,2,……,L,然后对模型进行校验,并输出结果。
优选地,所述步骤⑥中,混沌优化搜索的过程如下:
①令l=0,l表示步数,随机生成D个不同轨迹的混沌变量不包括混沌迭代方程的4个不动点(0,0.25,0.5,0.75,1),其中:d表示变量的序号,l表示第l次混沌搜索;
②将按照(32)式线性映射到优化变量取值区间[ad,bd]得到ad,bd为优化变量的取值范围:
③对进行混沌搜索:
若则其中,f*为当前最优解,为当前得到的最优变量,β为常数;
④
⑤重复步骤②-④,直到一定步数内f*保持不变或者达到给定的最大搜索步数ckmax,结束寻优计算,此时的即为算法得到的最优变量,f*为得到的最优解;
⑥d←d+1,如果d=D则结束混沌搜索,否则转到步骤①对下一变量进行混沌搜索。
本发明把配电网络(含110k V及以下配电变压器)作为综合负荷模型的有机组成部分,同时计入配网无功补偿和变压器分接头的影响,在配电网综合等效的基础上,建立了“考虑配电网参数及无功补偿的综合异步电动机负荷模型”,然后初始化模型,建立辨识准则函数,并对负荷模型中待辨识的负荷参数通过输入电压和实际系统的输出响应采用CPSO算法进行寻优,大量基于现场实测负荷特性数据的建模实践表明,该模型能有效地描述综合负荷特性,在参数稳定性、泛化能力等性能上,较传统的异步电动机模型有了一定的改进,更加符合电网仿真所要求的应用环境,具有收敛速度快及辨识结果适应性能好的优点;在建立负荷模型时,考虑了计入配网无功补偿和变压器分接头的影响,更符合实际情况,使得负荷模型更加准确;利用进化类算法进行参数辨识优化可提高模型参数的准确性,把混沌优化搜索技术引入到粒子群(PSO)算法中,提出了基于混沌搜索的粒子群优化算法(CPSO),该算法具有收敛速度快的优点,有效避免了基本粒子群算法易陷入局部极小值的缺陷,同时CPSO算法中引入了混沌优化搜索的概率Pk,以对最优粒子Pbest进行混沌优化搜索,提高了算法收敛精度,有效地提高了模型参数辨识的准确性。
附图说明
图1为本发明所述综合负荷系统的等效结构示意图;
图2为本发明所述配电网综合感应电动机模型等值电路图;
图3为本发明所述基于CPSO方法综合感应电动机模型辨识流程图;
图4为所述实施例一中有功功率拟合图;
图5为所述实施例一中无功功率拟合图;
图6为所述实施例一中参数辨识适应度曲线图。
具体实施方式
如图1至图3所示,本发明所述的一种配电网综合等效异步电动机负荷模型参数辨识方法,依次包括以下步骤:
步骤一、考虑配电网参数及无功补偿,建立综合异步电动机负荷模型;
(1)在110kV线路上增加一个变压器,在10kV母线并联异步电动机,同时考虑变压器分接头引入基准变换以进行配电网络与异步电动机对接实现,设异步电动机在自身基准下的端电压、暂态内电势、从电网吸收的电流、功率及等值阻抗的物理量分别为: Uim、Ei′m、Iim、Pim、Qim、Zim,下标im为区分作用,并非变量,各个物理量相对应系统基准下的标么值分别为US、ES'、IS、PD、QD、ZS,下标S和D为区分作用,并非变量,设变压器高压侧实际运行分接头电压和低压侧额定电压分别为和与UBS对应的变压器低压侧系统基准电压为UBL,其中,UBS表示系统电压基准,下标T、TN、BS及BL均为区分作用,并非变量,则有:
K=SBS/SBM (2)
UBL=UBS/kT (3)
其中,kT*表示变压器高压侧与低压侧的非标准电压比的标幺值,kT表示变压器高压侧与低压侧的非标准电压比,K表示系统与异步电动机之基准功率比,kB表示系统电压基准与感应电动机电压基准比值,SBS表示系统功率,SBM表示感应电动机功率,下标B、BS、BM均为区分作用,并非变量;
由式(1)可知,当取时,若变压器运行于主抽头,即或取各基准电压为网络平均额定电压且变压器变比为平均额定电压比,则kT*=1,设异步电动机电压基准为则由电压等级归算及标幺值换算原理,即可推导出“基准变换”的变换矩阵如式(4)所示:
其中,各个物理量的定义如步骤(1)所述;
配电网综合等效后的线路-变压器组首、末端功率及电压平衡关系如式(5)~式(7),变压器低压侧母线L的功率平衡关系如公式(8),
PS+jQS=(PZIP+jQZIP)+(PD+jQD)-jQC (8)
其中,PS表示线路末端有功功率,QS表示线路末端无功功率,P表示线路首端有功功率, Q表示线路首端无功功率,Re表示线路的等值电阻,Xe表示线路的等值电抗,US表示异步电动机的端电压的标幺值,USx表示感应电动机在系统基准下的端电压在x轴的坐标分量,USy表示感应电动机在系统基准下的端电压在y轴的坐标分量,U表示线路首端电压,下标e、Sx 及Sy均为区分作用,并非变量;
(2)在10kV母线并联静态负荷,将静态负荷用多项式(ZIP)模型描述如式(9):
其中,模型参数满足ZP+IP+PP=1、ZQ+IQ+PQ=1及 (0)为 US(0)的有名值,US(0)表示异步电动机的初始端电压,PZIP(t)表示静态负荷在t时刻的有功功率,QZIP(t)表示静态负荷在t时刻的无功功率,PZIP(0)表示静态负荷的初始有功功率,QZIP(0)表示静态负荷的初始无功功率,ZP表示静态负荷的有功阻抗,ZQ表示静态负荷的无功阻抗,US(t)表示异步电动机在t时刻的端电压,kU0为已知的电压变换系数;
(3)在静态负荷上并联一附加动态无功补偿元件,方程式如下:
其中,QC(t)表示无功补偿元件在t时刻的无功功率,下标c为区分作用,并非变量,f 表示频率,XC表示电容器容抗,XC0表示角频率近似取1时的值,即正常运行时的电容器容抗,表示补偿电容器的电容,ω表示角频率;
(4)将异步电动机模型用三阶机电暂态微分方程描述,如式(11):
式中,e'im=e'im.x+je'im.y为异步电动机暂态电势,e′im·x及e′im·y分别表示暂态电势在x轴和 y轴的坐标分量,t表示时间,s表示异步电动机的转差率,w0表示系统的同步角频率,频率为50Hz时,w0=314.16rad/s,Td0'表示异步电动机暂态电势的衰减时间常数,下标d仅为区分作用,并非变量,Xs表示定子绕组漏抗,Rs表示定子绕组电阻,X'表示定子和转子暂态电抗,Uim=Uim.x+jUim.y表示负荷端电压,Uimx及Uimy分别表示感应电动机电压在x轴和y轴的坐标分量,Tj表示转子惯性时间常数,Xr表示转子绕组漏抗,Xm表示激磁电抗,Rr表示转子电阻,Tm表示异步电动机的机械负载功率,具体计算公式如式(12):
Tm=T0[Aωr 2+Bωr+C] (12)
其中,T0为异步电动机负载率;ωr=1-s为转子侧角速度;A、B、C为机械转矩系数,且满足A+B+C=1;
异步电动机输出电流方程为:
其中,Iim.x及Iim.y分别表示异步电动机输出电流在x轴和y轴的坐标分量;
异步电动机的功率为:
步骤二、给定实测数据,计算状态变量初始值;
(1)给定实测电压激励U(k)、实测有功功率P(k)及无功功率Q(k);
(2)计算异步电动机及静态负荷的初始有功功率,设异步电动机的动态负荷占空比为
其中,PD(0)表示感应电动机初始有功功率,PS(0)表示变压器低压侧母线初始有功功率,下标m、D及S均为区分作用,并非变量;
由式(15)及综合异步电动机等值电路可求得异步电动机及静态负荷初始稳态有功功率为
其中,Pim(0)表示异步电动机的初始稳态有功功率,PZIP(0)表示静态负荷的初始稳态有功功率;
(3)计算异步电动机的初始滑差,由综合异步电动机等值电路可知异步电动机稳态等值阻抗和导纳如式(17)所示:
其中,Zim(0)表示异步电动机的稳态等值阻抗,Yim(0)表示异步电动机的稳压等值导纳,Gim(0)表示异步电动机的稳态等值电导,Bim(0)表示异步电动机的稳态等值电纳, R=Rr/S(0),Rr表示转子电阻,S(0)表示初始转差率;
结合(16),式(17)可写为:
其中,Xsm、Xrs、Xrm、Rm、XP均无实际物理意义,为替代参数,其取值如下:
其中,
从而可得关于R的一元二次方程AR2+BR+C=0 (22)
A1、B1、C1、A2、B2、C2均无实际物理意义,为替代参数,其取值如下:
其中,A、B及C的值由式(23)确定:
解式(22)可得异步电动机初始转差率如下:
为避免仿真中开方出现复数问题,需要对式(24)进行以下约束:
(4)计算异步电动机及静态负荷初始无功功率;
异步电动机的初始转差率确定后,由异步电动机应满足的初始稳态条件并由式(4)可得系统基准下的异步电动机初始稳态无功功率:
再由节点功率平衡原理唯一确定静态负荷的初始稳态无功功率如式(26):
QZIP(0)=QS(0)-QD(0)=QS(0)-Qim(0)/K (26)
(5)计算异步电动机暂态电动势初始值;
由异步电动机暂态等效电路可求得暂态电动势初始值如式(27)所示,
其中初始电流依初始稳态条件按式(28)所示确定:
(6)计算表异步电动机的机械负载率T0:
其中,ωr0=1-s(0)。
步骤三、建立负荷模型的参数辨识准则函数,并利用CPSO算法对目标函数进行寻优,并输出负荷模型;
负荷模型的参数辨识准则函数如式(30)所示:
其中,x(t)表示系统状态向量,u(t)为系统输入向量,ym(k)=[Pm(k),Qm(k)]T为输入u(t)时负荷模型得到的输出响应,y(k)表示符合模型实际输出相应,k为采样开始时刻,L为采样数,α=[Re,Xe,ZP,IP,ZQ,IQ,XC,K,Km,KT*,Rs,XS,Rr,Xr,Rm,Xm,Tj,A,B] 表示独立待辨识参数,β=[YZIP0,PP,PQ,T0,C]表示非独立待辨识参数,其中,PP=1-ZP-IP,PQ=1-PQ-IQ,C=1-A-B,YZIP0表示静态负荷初始导纳值,QZIP0表示静态负荷初始无功功率,Us0表示异步电动机的初始端电压,GZIP0表示静态负荷初始电导,BZIP0表示静态负荷初始电纳。
负荷模型的参数辨识准则函数的计算过程如下:
(1)在模型中输入实测电压激励U(k)、实测有功功率P(k)及无功功率Q(k),并给定独立待辨识参数初值α,由实测数据通过式(5)~(7)得到系统基准下的末端电压激励,通过式(4)变换为自身基准下的末端电压激励;
(2)由步骤二计算出异步电动机及静态负荷初始功率响应,求出状态变量初值和非独立参数,利用改进欧拉法解微分方程组(11),求得系统转差率S(k)及暂态电动势eimx(k),并带入输出方程(13)~(14),求得感应电动机的功率响应,由式(9)求得静态负荷的功率响应,再根据式(4)~(8)求得异步电动机在自身基准下的功率响应及线路末端功率响应,最终得到线路首端功率响应,即式(6);
(3)采用改进CPSO算法对以下目标函数寻优:
其中,各个参数的含义与式(30)相同,不再赘述。
改进CPSO算法以PSO算法的计算流程为主体流程,对种群中最好的粒子gbest进行给定步数的混沌优化搜索,指导粒子群往最优解方向搜索,以改进SPSO算法进化后期收敛速度慢,易陷入局部极小的缺点,为了提高计算效率,在进化初期,由于PSO算法的收敛速度较快,以小概率混沌搜索最优粒子的优化变量,在进化后期,以接近1的概率调用混沌优化搜索。
利用改进CPSO算法进行目标函数寻优的过程如下:
Step1、根据经验值对各个参数初始化,包括学习因子c1和c2,约束因子a,最大进化代数itermax,wmin和wmax,粒子数N,混沌搜索步长调节参数β和混沌搜索步数ckmax,其中,下标max及min均为区分作用,并非变量;
Step2、随机生成N个粒子的xi和vi,令k=0;
Step3、按式计算wk;
Step4、根据步骤二中求得的模型的有功功率计无功功率,计算粒子适应度,即式(30),并设定粒子的Pbest及种群的gbest,Pbest表示种群搜索到的最优解,gbest表示种群中最好的粒子;
Step5、按式更新粒子的速度和位置,对种群中的粒子进行操作,并更新每个粒子的pi和种群的g,式中:k表示第k次迭代,c1、c2为学习因子,r1、r2是[0,1]间的随机数,a为控制速度权重的约束因子,wk为惯性权重;如果则取如果则取如果或者则重新初始化 和为vd的取值范围,下标d及id均仅为区分作用,并非变量;
Step6、按下式计算混沌优化搜索的概率Pk
如果rand(0,1)≤Pk,则令d=1,对粒子g中变量进行混沌优化搜索,其余的D-1个变量保持不变,之后进入Step7,否则直接进入Step7,其中,rand(0,1)为[0,1]间的随机数;
其中,混沌优化搜索的过程如下:
①令l=0,l表示步数,随机生成D个不同轨迹的混沌变量不包括混沌迭代方程的4个不动点(0,0.25,0.5,0.75,1),其中:d表示变量的序号,l表示第l次混沌搜索;
②将按照(32)式线性映射到优化变量取值区间[ad,bd]得到ad,bd为优化变量的取值范围:
③对进行混沌搜索:
若则其中:f*为当前最优解,为当前得到的最优变量,β为常数;
④
⑤重复步骤②-④,直到一定步数内f*保持不变或者达到给定的最大搜索步数ckmax,结束寻优计算,此时的即为算法得到的最优变量,f*为得到的最优解;
⑥d←d+1,如果d=D则结束混沌搜索,否则转到步骤①对下一变量进行混沌搜索。
Step7、令k=k+1,此时若k<L,则按照式(30)将此函数模型响应与实测功率进行计算,然后进入Step8,若k<L,则返回Step5;
Step8、收敛判断,如果粒子g的适应值小于给定的阀值或者k>itermax,则进化过程成功结束,然后进入Step9,否则求解第i此最优化问题minJ→αi,并令i=i+1,然后返回Step1;
Step9、令α=α(i-1),β=β(i-1),输出参数α、β及模型响应序列ym(k),其中 k=1,2,……,L。
(4)对模型进行校验,并输出结果。
对模型的参数进行校验的过程为现有技术,不再赘述。
以下结合实施例对本发明的方案进行具体描述。
实施例一
对某变电站现场实测数据进行辨识,数据特征见表1:
表1实测动态负荷记录数据特征
设静态负荷为恒阻抗特性,即ZP=1;IP=QP=0;ZQ=1;IQ=PQ=0,取则独立待辨识参数为α=[Re,Xe,XC,K,Km,RS,XS,Rr,Xr,Rm,Xm,Tj,A,B]14个参数,其中非重点参数取典型值,见表2:
表2非重点参数典型值
设定参数的搜索范围:Km为0.1~0.8;Tj为1.2~3.2s;Xs为0.1~0.4pu;Xe为0.01~0.2pu; Re为0.01~0.1pu;Tj为1.2~3.2s,辨识结果见表3:
表3参数辨识结果
编号2的数据拟合曲线见图4、图5,适应度曲线见图6,通过图4至图6可看出,六月CPSO算法辨识该模型,其拟合数据与实测数据基本一致,该算法收敛速度快,辨识结果适应性能好,验证了算法在动态负荷模型参数辨识中的有效性。
Claims (8)
1.一种配电网综合等效异步电动机负荷模型参数辨识方法,其特征在于,依次包括以下步骤:
(1)考虑配电网参数及无功补偿,建立综合异步电动机负荷模型;
(2)给定实测数据,计算状态变量初始值;
(3)建立负荷模型的参数辨识准则函数,利用CPSO算法对目标函数进行寻优,并输出负荷模型;
(4)对负荷模型进行校验,并输出结果。
2.如权利要求1所述一种配电网综合等效异步电动机负荷模型参数辨识方法,其特征在于,所述步骤(1)具体包括以下步骤:
①在110kV线路上增加一个变压器,在10kV母线上并联异步电动机,同时考虑变压器分接头引入基准变换以进行配电网络与异步电动机对接实现,设异步电动机在自身基准下的端电压、暂态内电势、从电网吸收的电流、功率及等值阻抗的物理量分别为:Uim、E′im、Iim、Pim、Qim、Zim,下标im为区分作用,并非变量,各个物理量相对应系统基准下的标么值分别为US、ES'、IS、PD、QD、ZS,下标S和D为区分作用,并非变量,设变压器高压侧实际运行分接头电压和低压侧额定电压分别为和与UBS对应的变压器低压侧系统基准电压为UBL,其中,UBS表示系统电压基准,下标T、TN、BS及BL均为区分作用,并非变量,则有:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>T</mi>
<mo>*</mo>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<mo>/</mo>
<msub>
<mi>k</mi>
<mi>B</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>U</mi>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>U</mi>
<mrow>
<mi>T</mi>
<mi>N</mi>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>B</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>B</mi>
<mi>M</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
K=SBS/SBM (2)
UBL=UBS/kT (3)
其中,kT*表示变压器高压侧与低压侧的非标准电压比的标幺值,kT表示变压器高压侧与低压侧的非标准电压比,K表示系统与异步电动机之基准功率比,kB表示系统电压基准与感应电动机电压基准比值,SBS表示系统功率,SBM表示感应电动机功率,下标B、BS、BM均为区分作用,并非变量;
由式(1)可知,当取时,若变压器运行于主抽头,即或取各基准电压为网络平均额定电压且变压器变比为平均额定电压比,则kT*=1,设异步电动机电压基准为则由电压等级归算及标幺值换算原理,即可推导出“基准变换”的变换矩阵如式(4)所示:
<mrow>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>U</mi>
<mi>S</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>E</mi>
<mi>S</mi>
<mo>&prime;</mo>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>I</mi>
<mi>S</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mi>D</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mi>S</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>k</mi>
<mrow>
<mi>T</mi>
<mo>*</mo>
</mrow>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>k</mi>
<mrow>
<mi>T</mi>
<mo>*</mo>
</mrow>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>Kk</mi>
<mrow>
<mi>T</mi>
<mo>*</mo>
</mrow>
</msub>
</mrow>
</mfrac>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mi>K</mi>
</mfrac>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mi>K</mi>
</mfrac>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>Kk</mi>
<mrow>
<mi>T</mi>
<mo>*</mo>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,各个物理量的定义如步骤(1)所述;
配电网综合等效后的线路-变压器组首、末端功率及电压平衡关系如式(5)~式(7),变压器低压侧母线L的功率平衡关系如公式(8),
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>S</mi>
</msub>
<mo>=</mo>
<mi>P</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>P</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>Q</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>/</mo>
<msup>
<mi>U</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>S</mi>
</msub>
<mo>=</mo>
<mi>Q</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>P</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>Q</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>X</mi>
<mi>e</mi>
</msub>
<mo>/</mo>
<msup>
<mi>U</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>P</mi>
<mo>=</mo>
<msub>
<mi>P</mi>
<mi>S</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>/</mo>
<msubsup>
<mi>U</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Q</mi>
<mo>=</mo>
<msub>
<mi>Q</mi>
<mi>S</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>X</mi>
<mi>e</mi>
</msub>
<mo>/</mo>
<msubsup>
<mi>U</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>S</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>U</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>PR</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>QX</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>U</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>S</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>QR</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>PX</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>U</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
PS+jQS=(PZIP+jQZIP)+(PD+jQD)-jQC (8)
其中,PS表示线路末端有功功率,Qs表示线路末端无功功率,P表示线路首端有功功率,Q表示线路首端无功功率,Re表示线路的等值电阻,Xe表示线路的等值电抗,Us表示异步电动机的端电压的标幺值,Usx表示感应电动机在系统基准下的端电压在x轴的坐标分量,Usy表示感应电动机在系统基准下的端电压在y轴的坐标分量,U表示线路首端电压,下标e、Sx及Sy均为区分作用,并非变量;
②在10kV母线并联静态负荷,将静态负荷用多项式(ZIP)模型描述如式(9):
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>Z</mi>
<mi>I</mi>
<mi>P</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>Z</mi>
<mi>I</mi>
<mi>P</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<msub>
<mi>Z</mi>
<mi>P</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>U</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>U</mi>
<mi>S</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>P</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>U</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>U</mi>
<mi>S</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>P</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>Z</mi>
<mi>I</mi>
<mi>P</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>Z</mi>
<mi>I</mi>
<mi>P</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<msub>
<mi>Z</mi>
<mi>Q</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>U</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>U</mi>
<mi>S</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>Q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>U</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>U</mi>
<mi>S</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>Q</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,模型参数满足ZP+IP+PP=1、ZQ+IQ+PQ=1及 为US(0)的有名值,US(0)表示异步电动机的初始端电压,PZIP(t)表示静态负荷在t时刻的有功功率,QZIP(t)表示静态负荷在t时刻的无功功率,PZIP(0)表示静态负荷的初始有功功率,QZIP(0)表示静态负荷的初始无功功率,ZP表示静态负荷的有功阻抗,ZQ表示静态负荷的无功阻抗,US(t)表示异步电动机在t时刻的端电压,kU0为已知的电压变换系数;
③在静态负荷上并联一附加动态无功补偿元件,方程式如下:
<mrow>
<msub>
<mi>Q</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>U</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>X</mi>
<mi>C</mi>
</msub>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>U</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>X</mi>
<mrow>
<mi>C</mi>
<mn>0</mn>
</mrow>
</msub>
</mfrac>
<mi>f</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,QC(t)表示无功补偿元件在t时刻的无功功率,下标c为区分作用,并非变量,f表示频率,在标幺值下XC表示电容器容抗,XC0表示角频率近似取1时的值,即正常运行时的电容器容抗,表示补偿电容器的电容,ω表示角频率;
④将异步电动机模型用三阶机电暂态微分方程描述,如式(11):
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msub>
<msup>
<mi>de</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<msup>
<mi>e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>{</mo>
<msub>
<msup>
<mi>e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
</mrow>
<mrow>
<msubsup>
<mi>R</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>&omega;X</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>&omega;e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>&omega;X</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>&omega;e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msub>
<msup>
<mi>de</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<msup>
<mi>e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>{</mo>
<msub>
<msup>
<mi>e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
</mrow>
<mrow>
<msubsup>
<mi>R</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>&omega;e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>&omega;X</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>&omega;e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>m</mi>
</msub>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<msub>
<msup>
<mi>e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>+</mo>
<msup>
<mi>&omega;X</mi>
<mo>&prime;</mo>
</msup>
<msub>
<msup>
<mi>e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>&omega;e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<msub>
<msup>
<mi>e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msup>
<mi>&omega;X</mi>
<mo>&prime;</mo>
</msup>
<msub>
<msup>
<mi>e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>&omega;e</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>T</mi>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>X</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>R</mi>
<mi>r</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<msub>
<mi>X</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>m</mi>
</msub>
<msub>
<mi>X</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<msub>
<mi>X</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>r</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,e'im=e'im.x+je'im.y为异步电动机暂态电势,e′im·x及e′im·y分别表示暂态电势在x轴和y轴的坐标分量,t表示时间,s表示异步电动机的转差率,w0表示系统的同步角频率,频率为50Hz时,w0=314.16rad/s,Td0'表示异步电动机暂态电势的衰减时间常数,下标d仅为区分作用,并非变量,Xs表示定子绕组漏抗,Rs表示定子绕组电阻,X'表示定子和转子暂态电抗,Uim=Uim.x+jUim.y表示负荷端电压,Uimx及Uimy分别表示感应电动机电压在x轴和y轴的坐标分量,Tj表示转子惯性时间常数,Xr表示转子绕组漏抗,Xm表示激磁电抗,Rr表示转子电阻,Tm表示异步电动机的机械负载功率,具体计算公式如式(12):
Tm=T0[Aωr 2+Bωr+C] (12)
其中,T0为异步电动机负载率;ωr=1-s为转子侧角速度;A、B、C为机械转矩系数,且满足A+B+C=1;
异步电动机输出电流方程为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>&omega;e</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>&omega;X</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>&omega;e</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msubsup>
<mi>R</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>&omega;X</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>&omega;e</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>&omega;X</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>&omega;e</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msubsup>
<mi>R</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>&omega;X</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Iim.x及Iim.y分别表示异步电动机输出电流在x轴和y轴的坐标分量;
异步电动机的功率为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
3.如权利要求2所述一种配电网综合等效异步电动机负荷模型参数辨识方法,其特征在于,所述步骤(2)具体包括以下步骤:
①给定实测电压激励U(k)、实测有功功率P(k)及无功功率Q(k);
②计算异步电动机及静态负荷的初始有功功率,设异步电动机的动态负荷占空比为
<mrow>
<msub>
<mi>K</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mi>S</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,PD(0)表示感应电动机初始有功功率,Ps(0)表示变压器低压侧母线初始有功功率,下标m、D及S均为区分作用,并非变量;
由式(15)及综合异步电动机等值电路可求得异步电动机及静态负荷初始稳态有功功率为
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>K</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>K</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>K</mi>
<mi>m</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mi>S</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>Z</mi>
<mi>I</mi>
<mi>P</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<mi>S</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>D</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>K</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mi>S</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Pim(0)表示异步电动机的初始稳态有功功率,PZIP(0)表示静态负荷的初始稳态有功功率;
③计算异步电动机的初始滑差,由综合异步电动机等值电路可知异步电动机稳态等值阻抗和导纳如式(17)所示:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mi>S</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>R</mi>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
<mo>+</mo>
<mo>(</mo>
<mi>R</mi>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>jB</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Zim(0)表示异步电动机的稳态等值阻抗,Yim(0)表示异步电动机的稳压等值导纳,Gim(0)表示异步电动机的稳态等值电导,Bim(0)表示异步电动机的稳态等值电纳,R=Rr/S(0),Rr表示转子电阻,S(0)表示初始转差率;
结合(16),式(17)可写为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mi>R</mi>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mi>R</mi>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>k</mi>
<mrow>
<mi>T</mi>
<mo>*</mo>
</mrow>
<mn>2</mn>
</msubsup>
<msub>
<mi>P</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<msub>
<mi>U</mi>
<mi>S</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<msubsup>
<mi>k</mi>
<mrow>
<mi>T</mi>
<mo>*</mo>
</mrow>
<mn>2</mn>
</msubsup>
<mfrac>
<mrow>
<mi>K</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>K</mi>
<mi>m</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mi>S</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<msub>
<mi>U</mi>
<mi>S</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>s</mi>
<mi>m</mi>
</mrow>
</msub>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<mi>R</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>X</mi>
<mrow>
<mi>r</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msub>
<msub>
<mi>X</mi>
<mi>P</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mi>R</mi>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Xsm、Xrs、Xrm、Rm、XP均无实际物理意义,为替代参数,其取值如下:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>2</mn>
<mo>&lsqb;</mo>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>X</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>X</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<msub>
<mi>X</mi>
<mi>m</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>R</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<msubsup>
<mi>X</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>X</mi>
<mrow>
<mi>r</mi>
<mi>S</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>X</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>X</mi>
<mi>P</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>X</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>X</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<msubsup>
<mi>X</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>r</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>X</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>X</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>S</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mi>P</mi>
</msub>
<mo>=</mo>
<msub>
<mi>X</mi>
<mi>r</mi>
</msub>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>r</mi>
</msub>
<msub>
<mi>X</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<msub>
<mi>X</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
4
从而可得关于R的一元二次方程AR2+BR+C=0 (22)
A1、B1、C1、A2、B2、C2均无实际物理意义,为替代参数,其取值如下:
其中,A、B及C的值由式(23)确定:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>A</mi>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>B</mi>
<mo>=</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>C</mi>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
解式(22)可得异步电动机初始转差率如下:
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>AR</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mo>-</mo>
<mi>B</mi>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<mi>B</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>4</mn>
<mi>A</mi>
<mi>C</mi>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>24</mn>
<mo>)</mo>
</mrow>
</mrow>
④计算异步电动机及静态负荷初始无功功率;
异步电动机的初始转差率确定后,由异步电动机应满足的初始稳态条件并由式(4)可得系统基准下的异步电动机初始稳态无功功率:
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mi>L</mi>
</msub>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
<mo>/</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>T</mi>
<mo>*</mo>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>25</mn>
<mo>)</mo>
</mrow>
</mrow>
再由节点功率平衡原理唯一确定静态负荷的初始稳态无功功率如式(26):
QZIP(0)=QS(0)-QD(0)=QS(0)-Qim(0)/K (26)
⑤计算异步电动机暂态电动势初始值;
由异步电动机暂态等效电路可求得暂态电动势初始值如式(27)所示,
<mrow>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>e</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>e</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>{</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
<mtd>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>27</mn>
<mo>)</mo>
</mrow>
</mrow>
其中初始电流依初始稳态条件按式(28)所示确定:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mo>.</mo>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>28</mn>
<mo>)</mo>
</mrow>
</mrow>
⑥计算表异步电动机的机械负载率T0:
<mrow>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>r</mi>
<mn>0</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msup>
<msub>
<mi>A&omega;</mi>
<mrow>
<mi>r</mi>
<mn>0</mn>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>B&omega;</mi>
<mrow>
<mi>r</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<mi>C</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>29</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,ωr0=1-s(0)。
4.如权利要求3所述一种配电网综合等效异步电动机负荷模型参数辨识方法,其特征在于:为避免仿真中开方出现复数问题,需要对式(24)进行以下约束:
①A>0,B<0且B2-4AC≥0;
②额定初始滑差率取值范围为:0.015≤s≤0.07。
5.如权利要求4所述一种配电网综合等效异步电动机负荷模型参数辨识方法,其特征在于:所述步骤(3)中,负荷模型的参数辨识准则函数如式(30)所示:
<mrow>
<mi>min</mi>
<mi> </mi>
<mi>J</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<mo>{</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>30</mn>
<mo>)</mo>
</mrow>
</mrow>
5
其中,x(t)表示系统状态向量,u(t)为系统输入向量,ym(k)=[Pm(k),Qm(k)]T为输入u(t)时负荷模型得到的输出响应,y(k)表示符合模型实际输出相应,k为采样开始时刻,L为采样数,α=[Re,Xe,ZP,IP,ZQ,IQ,XC,K,Km,KT*,RS,XS,Rr,Xr,Rm,Xm,Tj,A,B]表示独立待辨识参数,β=[YZIP0,PP,PQ,T0,C]表示非独立待辨识参数,其中,PP=1-ZP-IP,PQ=1-PQ-IQ,C=1-A-B,YZIP0表示静态负荷初始导纳值,QZIP0表示静态负荷初始无功功率,Us0表示异步电动机的初始端电压,GZIP0表示静态负荷初始电导,BZIP0表示静态负荷初始电纳。
6.如权利要求5所述一种配电网综合等效异步电动机负荷模型参数辨识方法,其特征在于,所述负荷模型的参数辨识准则函数的计算过程如下:
①在实际变电站系统中输入电压激励U(k)、实测有功功率P(k)及无功功率Q(k),并给定独立待辨识参数初值α,由实测数据通过式(5)~(7)得到系统基准下的末端电压激励,通过式(4)变换为自身基准下的末端电压激励;
②由步骤(2)计算出异步电动机及静态负荷初始功率响应,求出状态变量初值和非独立参数,利用改进欧拉法解微分方程组(11),求得系统转差率S(k)及暂态电动势eimx(k),并带入输出方程(13)~(14),求得感应电动机的功率响应,由式(9)求得静态负荷的功率响应,再根据式(4)~(8)求得异步电动机在自身基准下的功率响应及线路末端功率响应,最终得到线路首端功率响应,即式(6);
③采用改进CPSO算法对以下目标函数寻优:
<mrow>
<mi>min</mi>
<mi> </mi>
<mi>J</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<mo>{</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>30</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,各个参数的含义与权利要求5中相同,不再赘述。
7.如权利要求6所述一种配电网综合等效异步电动机负荷模型参数辨识方法,其特征在于,利用改进CPSO算法进行目标函数寻优的过程如下:
①根据经验值对各个参数初始化,包括学习因子c1和c2,约束因子a,最大进化代数itermax,wmin和wmax,粒子数N,混沌搜索步长调节参数β和混沌搜索步数ckmax,其中,下标max及min均为区分作用,并非变量;
②随机生成N个粒子的xi和vi,令k=0;
③按式计算wk;
④根据步骤(2)中求得的模型的有功功率计无功功率,计算粒子适应度,即式(30),并设定粒子的Pbest及种群的gbest,Pbest表示种群搜索到的最优解,gbest表示种群中最好的粒子;
⑤按式更新粒子的速度和位置,对种群中的粒子进行操作,并更新每个粒子的pi和种群的g,式中:k表示第k次迭代,c1、c2为学习因子,r1、r2是[0,1]间的随机数,a为控制速度权重的约束因子,wk为惯性权重;如果则取如果则取如果或者则重新初始化 和为vd的取值范围,下标d及id均仅为区分作用,并非变量;
⑥按下式计算混沌优化搜索的概率Pk
<mrow>
<msub>
<mi>P</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>31</mn>
<mo>)</mo>
</mrow>
</mrow>
如果rand(0,1)≤Pk,则令d=1,对粒子g中变量进行混沌优化搜索,其余的D-1个变量保持不变,之后进入步骤⑦,否则直接进入步骤⑦,其中,rand(0,1)为[0,1]间的随机数;
⑦令k=k+1,此时若k≥L,则按照式(30)将此函数模型响应与实测功率进行计算,然后进入步骤⑧,若k<L,则返回步骤⑤;
⑧收敛判断,如果粒子g的适应值小于给定的阀值或者k>itermax,则进化过程成功结束,然后进入步骤⑨,否则求解第i此最优化问题minJ→αi,并令i=i+1,然后返回步骤①;
⑨令α=α(i-1),β=β(i-1),输出参数α、β及模型响应序列ym(k),其中k=1,2,……,L,然后对模型进行校验,并输出结果。
8.如权利要求7所述一种配电网综合等效异步电动机负荷模型参数辨识方法,其特征在于:所述步骤⑥中,混沌优化搜索的过程如下:
①令l=0,l表示步数,随机生成D个不同轨迹的混沌变量不包括混沌迭代方程的4个不动点(0,0.25,0.5,0.75,1),其中:d表示变量的序号,l表示第l次混沌搜索;
②将按照(32)式线性映射到优化变量取值区间[ad,bd]得到ad,bd为优化变量的取值范围:
<mrow>
<msubsup>
<mi>rx</mi>
<mi>d</mi>
<mi>l</mi>
</msubsup>
<mo>&LeftArrow;</mo>
<msubsup>
<mi>a</mi>
<mi>d</mi>
<mi>l</mi>
</msubsup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mi>d</mi>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>cx</mi>
<mi>d</mi>
<mi>l</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>32</mn>
<mo>)</mo>
</mrow>
</mrow>
③对进行混沌搜索:
若则其中,f*为当前最优解,为当前得到的最优变量,β为常数;
④l←l+1,
⑤重复步骤②-④,直到一定步数内f*保持不变或者达到给定的最大搜索步数ckmax,结束寻优计算,此时的即为算法得到的最优变量,f*为得到的最优解;
⑥d←d+1,如果d=D则结束混沌搜索,否则转到步骤①对下一变量进行混沌搜索。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710443376.3A CN107086606B (zh) | 2017-06-13 | 2017-06-13 | 一种配电网综合等效异步电动机负荷模型参数辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710443376.3A CN107086606B (zh) | 2017-06-13 | 2017-06-13 | 一种配电网综合等效异步电动机负荷模型参数辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107086606A true CN107086606A (zh) | 2017-08-22 |
CN107086606B CN107086606B (zh) | 2021-07-02 |
Family
ID=59606540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710443376.3A Active CN107086606B (zh) | 2017-06-13 | 2017-06-13 | 一种配电网综合等效异步电动机负荷模型参数辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107086606B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107677960A (zh) * | 2017-09-21 | 2018-02-09 | 国网福建省电力有限公司 | 基于pmu量测的同步发电机参数辨识方法 |
CN108388120A (zh) * | 2018-02-12 | 2018-08-10 | 沈阳工业大学 | 一种用于混沌系统参数辨识的方法 |
CN108564205A (zh) * | 2018-03-27 | 2018-09-21 | 昆明理工大学 | 一种基于实测数据的负荷模型及参数辨识优化方法 |
CN108649562A (zh) * | 2018-05-04 | 2018-10-12 | 华北水利水电大学 | 一种基于苍狼算法的电力系统负荷建模参数辨识方法 |
CN109713662A (zh) * | 2018-12-20 | 2019-05-03 | 清华大学 | 一种电力系统负荷模型辨识参数向低压节点等效的方法 |
CN109830965A (zh) * | 2019-03-18 | 2019-05-31 | 重庆大学 | 基于节点电动机最大自起动容量的低压减载量的计算方法 |
CN111162719A (zh) * | 2020-02-06 | 2020-05-15 | 深圳市测力佳控制技术有限公司 | 异步电机运行参数测量方法、装置和计算机设备 |
CN111241749A (zh) * | 2020-01-13 | 2020-06-05 | 广西师范大学 | 一种基于储备池计算的永磁同步电动机混沌预测方法 |
CN111725804A (zh) * | 2020-06-10 | 2020-09-29 | 国网新疆电力有限公司 | 一种综合负荷模型参数辨识方法和计算机可读存储介质 |
CN112039051A (zh) * | 2020-06-30 | 2020-12-04 | 长沙理工大学 | 一种含双馈风力发电机接入变电站母线负荷的实时建模方法 |
CN113111298A (zh) * | 2021-04-01 | 2021-07-13 | 北京兴达智联科技有限公司 | 一种升降压变换器电路参数在线辨识的方法 |
CN113452299A (zh) * | 2021-06-23 | 2021-09-28 | 云南电网有限责任公司电力科学研究院 | 一种基于定子电流的动态负荷建模参数辨识方法 |
CN113705145A (zh) * | 2021-08-23 | 2021-11-26 | 贵州电网有限责任公司 | 一种π型传输线模型参数的在线自校正方法 |
CN115102235A (zh) * | 2022-04-29 | 2022-09-23 | 华南理工大学 | 一种基于交流调压器的户用光伏优化管理方法及系统 |
CN116504069A (zh) * | 2023-06-26 | 2023-07-28 | 中国市政工程西南设计研究总院有限公司 | 城市道路网络容量优化方法、装置、设备及可读存储介质 |
CN117911195A (zh) * | 2024-03-18 | 2024-04-19 | 国网山西省电力公司电力科学研究院 | 一种异步电动机机群等效模型参数辨识方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103793556A (zh) * | 2014-01-07 | 2014-05-14 | 广西大学 | 综合出线负荷的变电站负荷建模方法 |
-
2017
- 2017-06-13 CN CN201710443376.3A patent/CN107086606B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103793556A (zh) * | 2014-01-07 | 2014-05-14 | 广西大学 | 综合出线负荷的变电站负荷建模方法 |
Non-Patent Citations (3)
Title |
---|
李欣然等: "配电网集结等效的异步电动机综合负荷模型及其总体测辨建模", 《电工技术学报》 * |
杨俊杰等: "基于混沌搜索的粒子群优化算法", 《计算机工程与应用》 * |
王振树等: "基于混沌与量子粒子群算法相结合的负荷模型参数辨识研究", 《电工技术学报》 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107677960B (zh) * | 2017-09-21 | 2020-01-10 | 国网福建省电力有限公司 | 基于pmu量测的同步发电机参数辨识方法 |
CN107677960A (zh) * | 2017-09-21 | 2018-02-09 | 国网福建省电力有限公司 | 基于pmu量测的同步发电机参数辨识方法 |
CN108388120A (zh) * | 2018-02-12 | 2018-08-10 | 沈阳工业大学 | 一种用于混沌系统参数辨识的方法 |
CN108564205A (zh) * | 2018-03-27 | 2018-09-21 | 昆明理工大学 | 一种基于实测数据的负荷模型及参数辨识优化方法 |
CN108649562A (zh) * | 2018-05-04 | 2018-10-12 | 华北水利水电大学 | 一种基于苍狼算法的电力系统负荷建模参数辨识方法 |
CN109713662A (zh) * | 2018-12-20 | 2019-05-03 | 清华大学 | 一种电力系统负荷模型辨识参数向低压节点等效的方法 |
CN109830965A (zh) * | 2019-03-18 | 2019-05-31 | 重庆大学 | 基于节点电动机最大自起动容量的低压减载量的计算方法 |
CN109830965B (zh) * | 2019-03-18 | 2021-02-02 | 重庆大学 | 基于节点电动机最大自起动容量的低压减载量的计算方法 |
CN111241749A (zh) * | 2020-01-13 | 2020-06-05 | 广西师范大学 | 一种基于储备池计算的永磁同步电动机混沌预测方法 |
CN111162719A (zh) * | 2020-02-06 | 2020-05-15 | 深圳市测力佳控制技术有限公司 | 异步电机运行参数测量方法、装置和计算机设备 |
CN111162719B (zh) * | 2020-02-06 | 2021-12-17 | 深圳市测力佳控制技术有限公司 | 异步电机运行参数测量方法、装置和计算机设备 |
CN111725804A (zh) * | 2020-06-10 | 2020-09-29 | 国网新疆电力有限公司 | 一种综合负荷模型参数辨识方法和计算机可读存储介质 |
CN112039051A (zh) * | 2020-06-30 | 2020-12-04 | 长沙理工大学 | 一种含双馈风力发电机接入变电站母线负荷的实时建模方法 |
CN113111298A (zh) * | 2021-04-01 | 2021-07-13 | 北京兴达智联科技有限公司 | 一种升降压变换器电路参数在线辨识的方法 |
CN113111298B (zh) * | 2021-04-01 | 2023-12-15 | 北京兴达智联科技有限公司 | 一种升降压变换器电路参数在线辨识的方法 |
CN113452299A (zh) * | 2021-06-23 | 2021-09-28 | 云南电网有限责任公司电力科学研究院 | 一种基于定子电流的动态负荷建模参数辨识方法 |
CN113452299B (zh) * | 2021-06-23 | 2024-02-06 | 云南电网有限责任公司电力科学研究院 | 一种基于定子电流的动态负荷建模参数辨识方法 |
CN113705145A (zh) * | 2021-08-23 | 2021-11-26 | 贵州电网有限责任公司 | 一种π型传输线模型参数的在线自校正方法 |
CN113705145B (zh) * | 2021-08-23 | 2024-01-26 | 贵州电网有限责任公司 | 一种π型传输线模型参数的在线自校正方法 |
CN115102235A (zh) * | 2022-04-29 | 2022-09-23 | 华南理工大学 | 一种基于交流调压器的户用光伏优化管理方法及系统 |
CN115102235B (zh) * | 2022-04-29 | 2024-04-16 | 华南理工大学 | 一种基于交流调压器的户用光伏优化管理方法及系统 |
CN116504069A (zh) * | 2023-06-26 | 2023-07-28 | 中国市政工程西南设计研究总院有限公司 | 城市道路网络容量优化方法、装置、设备及可读存储介质 |
CN116504069B (zh) * | 2023-06-26 | 2023-09-05 | 中国市政工程西南设计研究总院有限公司 | 城市道路网络容量优化方法、装置、设备及可读存储介质 |
CN117911195A (zh) * | 2024-03-18 | 2024-04-19 | 国网山西省电力公司电力科学研究院 | 一种异步电动机机群等效模型参数辨识方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107086606B (zh) | 2021-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107086606A (zh) | 一种配电网综合等效异步电动机负荷模型参数辨识方法 | |
WO2022127172A1 (zh) | 一种评估海上风电场柔性直流输电系统谐振稳定性的方法 | |
CN104917174B (zh) | 一种针对主动配电网的静态电压稳定性判定方法 | |
CN103197559A (zh) | 一种改善双馈机组接入后系统小干扰稳定性的方法 | |
CN106294993B (zh) | 一种考虑逆变器电流饱和的暂态能量函数分析方法 | |
CN104079228B (zh) | 双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法 | |
CN108336751B (zh) | 一种lcc-mmc混合直流电网的机电暂态建模方法 | |
Wang et al. | Adaptive extended Kalman filter based dynamic equivalent method of PMSG wind farm cluster | |
CN111541246B (zh) | 一种电力系统交直流潮流的全纯嵌入计算方法 | |
CN103984822A (zh) | 一种基于全数字实时仿真装置的三相潮流实现方法 | |
CN107666155A (zh) | 基于Markov模型的多能互补系统随机稳定性分析方法 | |
CN110206686A (zh) | 一种用于风力发电机组的自适应最大功率跟踪控制方法 | |
CN108649562A (zh) | 一种基于苍狼算法的电力系统负荷建模参数辨识方法 | |
CN105429132A (zh) | 一种电动机负荷模型的构建方法 | |
CN104617578B (zh) | 一种含风电场电力系统的可用输电能力的获取方法 | |
CN109617488A (zh) | 一种考虑励磁电路的虚拟同步机的建模方法 | |
CN108418242A (zh) | 一种基于相似度同调的双馈风电机动态等值方法 | |
CN106294959B (zh) | 模型参考自适应控制与双馈风机降阶模型相结合的建模仿真方法 | |
CN105576651B (zh) | 一种中小水电机群混合并行动态等值法 | |
CN106406272A (zh) | 一种风电场中静止无功发生器的控制器性能测试方法 | |
CN109149583A (zh) | 有源配电网简洁前推回代潮流计算方法 | |
CN105224728A (zh) | 一种含详细发电机模型的电网暂态稳定能量函数分析方法及系统 | |
CN108258725A (zh) | 一种基于等效功角同调的双馈风电机动态等值方法 | |
CN105140957B (zh) | 基于风电场和光伏电站聚合模型的机电振荡模式估算方法 | |
CN112736927A (zh) | 一种基于概率盒理论的风电场不确定性潮流建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |