CN107633110A - 一种负荷频率特性参数的确定方法和装置 - Google Patents
一种负荷频率特性参数的确定方法和装置 Download PDFInfo
- Publication number
- CN107633110A CN107633110A CN201710684435.6A CN201710684435A CN107633110A CN 107633110 A CN107633110 A CN 107633110A CN 201710684435 A CN201710684435 A CN 201710684435A CN 107633110 A CN107633110 A CN 107633110A
- Authority
- CN
- China
- Prior art keywords
- mrow
- load
- msub
- active power
- static load
- 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.)
- Pending
Links
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供了一种负荷频率特性参数的确定方法和装置,先计算单个设备类型中静态负荷的有功功率和单个设备类型中动态负荷的有功功率;然后计算静态负荷的综合有功功率和动态负荷的综合有功功率;最后确定静态负荷频率特性参数和动态负荷频率特性参数,静态负荷频率特性参数包括静态负荷的有功功率频率因子和无功功率频率因子,动态负荷频率特性参数包括动态负荷的机械转矩系数。本发明提供的负荷频率特性参数的确定方法确定的静态负荷频率特性参数和动态负荷频率特性参数能够提高电力系统的仿真精度,且符合实际物理电网,考虑负荷动态频率响应过程对电网稳定特性及控制策略的影响,为大型交直流混联电网安全稳运行与控制决策提供基础。
Description
技术领域
本发明涉及电力系统的安全稳定分析技术,具体涉及一种负荷频率特性参数的确定方法和装置。
背景技术
随着中国特高压交直流混联电网构建和形成,大规模潮流远距离穿越,多回大容量直流在送受端电网的集中接入,致使大电网电压稳定问题和频率稳定问题突出,负荷的电压无功特性和频率特性对大电网(即送端电网和受端电网)的安全稳定特性影响日益凸显,对大电网的安全稳定控制产生很大影响,负荷模型在大电网仿真中更为重要。
近年来发生的一系列停电事故也进一步突出了建立准确负荷模型的重要性。比如2003年发生在美国和加拿大的“8.14”大停电,2006年发生在欧洲的“11.4”事故等,在最后的调查总结报告中,均提到实际系统或者仿真系统所使用的电气元件模型不够准确,难以仿真分析事故特性。而且在1996年WSCC对“8.10”事故进行了数字仿真和模型校核比对,当修改了相关参数(包括负荷模型和参数)时,功率的仿真结果与故障录波数据基本一致,由此可见负荷模型对数字仿真结果有很大的影响。
在实际运行的电力系统中,通过实测可得频率动态过程曲线,但系统仿真结果与实测频率动态过程曲线有时会存在较大差异。1996年,美国西部协调委员会(WSCC)的事故分析报告中指出,采用不同的负荷模型进行仿真,将得到不同甚至截然相反的分析结果,这使人们认识到负荷模型对仿真计算的影响和重要性。
当系统发生故障造成功率不平衡时,频率会随之发生变化,尤其是在一些独立电网或者微网中,故障时频率变化往往较大,而电网的频率特性取决于负荷频率特性,因此,考虑频率特性的负荷模型结构及参数对正确认识微网或独立电网的系统频率动态特性十分重要。中国电网当前仿真采用的负荷模型及参数大多都是基于上世纪80年代左右的事故仿真确定的(局部电网陆续有所调整)。但是随着科技发展和产业结构的巨变,电网负荷构成、特性均发生巨大变化,尤其随着跨区混联电网发展,当前负荷模型参数仿真精度与现实存在较大偏差,导致当前负荷模型参数无法准确描述负荷频率特性。
发明内容
为了克服上述现有技术的不足,本发明提供一种负荷频率特性参数的确定方法和装置,先计算单个设备类型中静态负荷的有功功率和单个设备类型中动态负荷的有功功率;然后根据单个设备类型中静态负荷的有功功率计算静态负荷的综合有功功率,并根据单个设备类型中动态负荷的有功功率计算动态负荷的综合有功功率;最后根据静态负荷的综合有功功率确定静态负荷频率特性参数,并根据动态负荷的综合有功功率确定动态负荷频率特性参数,静态负荷频率特性参数和动态负荷频率特性参数能够提高电力系统的仿真精度。
为了实现上述发明目的,本发明采取如下技术方案:
本发明提供一种负荷频率特性参数的确定方法,所述负荷包括静态负荷和动态负荷,所述确定方法包括:
计算单个设备类型中静态负荷的有功功率和单个设备类型中动态负荷的有功功率;
根据单个设备类型中静态负荷的有功功率计算静态负荷的综合有功功率,并根据单个设备类型中动态负荷的有功功率计算动态负荷的综合有功功率;
根据静态负荷的综合有功功率确定静态负荷频率特性参数,并根据动态负荷的综合有功功率确定动态负荷频率特性参数;
所述静态负荷频率特性参数包括静态负荷的有功功率频率因子和无功功率频率因子;
所述动态负荷频率特性参数包括动态负荷的机械转矩系数。
所述单个设备类型中静态负荷的有功功率按下式计算:
PSi=Ni×NSi×PS
其中,Ni表示设备类型i的有功功率占负荷站点有功功率的百分比,Nsi为设备类型i中静态负荷的有功功率百分比,PS表示静态负荷的有功功率初值。
所述单个设备类型中动态负荷的有功功率按下式计算:
PRi=Ni×NMi×P
其中,NMi表示设备类型i中动态负荷的有功功率百分比,Nsi为设备类型i中静态负荷的有功功率百分比,P表示负荷站点的有功功率。
所述静态负荷的综合有功功率按下式计算:
其中,PSa表示静态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
所述动态负荷的综合有功功率按下式计算:
其中,PRa表示动态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
所述根据静态负荷的综合有功功率确定静态负荷频率特性参数包括:
根据静态负荷的综合有功功率确定如下式的静态负荷的有功功率频率因子:
其中,LDP表示静态负荷的有功功率频率因子,Pfi表示设备类型i中静态负荷的有功功率-频率特性系数;
根据静态负荷的综合有功功率确定如下式的静态负荷的无功功率频率因子:
其中,LDQ表示静态负荷的无功功率频率因子,Qfi表示设备类型i中静态负荷的无功功率-频率特性系数。
所述根据动态负荷的综合有功功率确定动态负荷频率特性参数包括:
根据动态负荷的综合有功功率确定如下式的动态负荷的机械转矩系数:
C=1-A-B
其中,A、B、C均表示动态负荷的机械转矩系数,Ai、Bi均表示设备类型i中动态负荷的机械转矩系数。
本发明提供一种负荷频率特性参数的确定装置,包括:
第一计算模块,用于计算单个设备类型中静态负荷的有功功率和单个设备类型中动态负荷的有功功率;
第二计算模块,用于根据单个设备类型中静态负荷的有功功率计算静态负荷的综合有功功率,并根据单个设备类型中动态负荷的有功功率计算动态负荷的综合有功功率;
确定模块,用于根据静态负荷的综合有功功率确定静态负荷频率特性参数,并根据动态负荷的综合有功功率确定动态负荷频率特性参数;
所述静态负荷频率特性参数包括静态负荷的有功功率频率因子和无功功率频率因子;
所述动态负荷频率特性参数包括动态负荷的机械转矩系数。
所述第一计算模块按下式计算单个设备类型中静态负荷的有功功率:
PSi=Ni×NSi×PS
其中,Ni表示设备类型i的有功功率占负荷站点有功功率的百分比,Nsi为设备类型i中静态负荷的有功功率百分比,PS表示静态负荷的有功功率初值。
所述第一计算模块按下式计算单个设备类型中动态负荷的有功功率:
PRi=Ni×NMi×P
其中,NMi表示设备类型i中动态负荷的有功功率百分比,Nsi为设备类型i中静态负荷的有功功率百分比,P表示负荷站点的有功功率。
所述第二计算模块按下式计算静态负荷的综合有功功率:
其中,PSa表示静态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
所述第二计算模块按下式计算动态负荷的综合有功功率:
其中,PRa表示动态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
所述确定模块具体用于:
根据静态负荷的综合有功功率确定如下式的静态负荷的有功功率频率因子:
其中,LDP表示静态负荷的有功功率频率因子,Pfi表示设备类型i中静态负荷的有功功率-频率特性系数;
根据静态负荷的综合有功功率确定如下式的静态负荷的无功功率频率因子:
其中,LDQ表示静态负荷的无功功率频率因子,Qfi表示设备类型i中静态负荷的无功功率-频率特性系数。
所述确定模块具体用于:
根据动态负荷的综合有功功率确定如下式的动态负荷的机械转矩系数:
C=1-A-B
其中,A、B、C均表示动态负荷的机械转矩系数,Ai、Bi均表示设备类型i中动态负荷的机械转矩系数。
与最接近的现有技术相比,本发明提供的技术方案具有以下有益效果:
本发明提供的负荷频率特性参数的确定方法中,先计算单个设备类型中静态负荷的有功功率和单个设备类型中动态负荷的有功功率;然后根据单个设备类型中静态负荷的有功功率计算静态负荷的综合有功功率,并根据单个设备类型中动态负荷的有功功率计算动态负荷的综合有功功率;最后根据静态负荷的综合有功功率确定静态负荷频率特性参数,并根据动态负荷的综合有功功率确定动态负荷频率特性参数;
本发明提供的负荷频率特性参数的确定方法确定的静态负荷频率特性参数和动态负荷频率特性参数能够提高电力系统的仿真精度,且符合实际物理电网,考虑负荷动态频率响应过程对电网稳定特性及控制策略的影响,为大型交直流混联电网安全稳运行与控制决策提供基础。
附图说明
图1是本发明实施例中负荷频率特性参数的确定方法流程图;
图2是本发明实施例中220kV变电站配电网络系统图;
图3是本发明实施例中仿真系统示意图;
图4是本发明实施例中系统的频率变化曲线图;
图5是本发明实施例中城西220kV负荷站点有功功率曲线图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
本发明实施例提供了一种负荷频率特性参数的确定方法,其中的负荷包括静态负荷和动态负荷,于是负荷频率特性参数就包括静态负荷频率特性参数和动态负荷频率特性参数。本发明实施例提供的负荷频率特性参数确定方法流程图如图1所示,具体过程如下:
S101:计算单个设备类型中静态负荷的有功功率和单个设备类型中动态负荷的有功功率;
S102:根据S101计算得到的单个设备类型中静态负荷的有功功率计算静态负荷的综合有功功率,并根据S101计算得到的单个设备类型中动态负荷的有功功率计算动态负荷的综合有功功率;
S103:根据S102计算得到的静态负荷的综合有功功率确定静态负荷频率特性参数,并根据S102计算得到的动态负荷的综合有功功率确定动态负荷频率特性参数;其中静态负荷频率特性参数包括静态负荷的有功功率频率因子和无功功率频率因子,动态负荷频率特性参数包括动态负荷的机械转矩系数。
上述S101中,单个设备类型中静态负荷的有功功率按下式计算:
PSi=Ni×NSi×PS
其中,Ni表示设备类型i的有功功率占负荷站点有功功率的百分比,Nsi为设备类型i中静态负荷的有功功率百分比,PS表示静态负荷的有功功率初值。
上述S101中,单个设备类型中动态负荷的有功功率按下式计算:
PRi=Ni×NMi×P
其中,NMi表示设备类型i中动态负荷的有功功率百分比,Nsi为设备类型i中静态负荷的有功功率百分比,P表示负荷站点的有功功率。
上述S102中,静态负荷的综合有功功率按下式计算:
其中,PSa表示静态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
上述S102中,动态负荷的综合有功功率按下式计算:
其中,PRa表示动态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
上述S103中,根据静态负荷的综合有功功率确定静态负荷频率特性参数具体过程如下:
1)根据静态负荷的综合有功功率确定如下式的静态负荷的有功功率频率因子:
其中,LDP表示静态负荷的有功功率频率因子,Pfi表示设备类型i中静态负荷的有功功率-频率特性系数;
2)根据静态负荷的综合有功功率确定如下式的静态负荷的无功功率频率因子:
其中,LDQ表示静态负荷的无功功率频率因子,Qfi表示设备类型i中静态负荷的无功功率-频率特性系数。
本发明实施例进一步根据上述S103得到的静态负荷的有功功率频率因子LDP和静态负荷的无功功率频率因子LDQ建立如下式的静态负荷模型:
其中,P1表示静态负荷的有功功率,PS表示静态负荷的有功功率初值,PZ表示恒定阻抗有功负荷比例,PI表示恒定电流有功负荷比例,PP表示恒定功率有功负荷比例;Q1表示静态负荷的无功功率,QS表示静态负荷的无功功率初值,QZ表示恒定阻抗无功负荷比例,QI表示恒定电流无功负荷比例,QP表示恒定功率无功负荷比例;V表示负荷母线实际电压幅值,V0为负荷母线电压幅值初值;Δf表示频率变化量。
上述S103中,根据动态负荷的综合有功功率确定如下式的动态负荷频率特性参数中的动态负荷的机械转矩系数:
C=1-A-B
其中,A、B、C均表示动态负荷的机械转矩系数,Ai、Bi均表示设备类型i中动态负荷的机械转矩系数。
本发明实施例进一步根据上述S103得到的动态负荷的机械转矩系数A、B、C建立如下式的动态负荷模型:
TE=Pe/ωr
上述动态负荷取电动机,其中,ωr表示电动机的实时转速,Pe表示电动机的电磁功率,TE表示电动机的电磁转矩,TM表示电动机的,T'0表示电动机的时间常数,E'd表示电动机的d轴电势,E'q表示电动机的q轴电势,Id表示电动机的d轴电流,Iq表示电动机的q轴电流,ωb表示电动机的基准转速,X表示电动机的等值电抗,X'表示电动机的次暂态电抗。
采用实际电网中某220kV变电站当前采用的负荷频率特性参数以及采用本发明实施例提供的确定方法得到的负荷频率特性参数和原系统(包括某变负荷区的110kV、35kV配电网络、无功补偿及110kV、35kV、10kV、6kV负荷节点的系统,如图2所示)进行仿真对比,验证本发明实施例提供的负荷频率特性参数确定方法的有效性。
采用如图3所示的仿真系统,Bus 1-Bus 4均为母线,一台同步发电机组通过双回线路向城西变和Bus 4供电,某变的有功负荷为167MW,Bus 4的有功负荷为40MW。
仿真条件:在仿真系统运行0.1秒时,Bus 4节点增加40MW有功负荷。
分别将图2所示的城西220kV变电站110kV及其以下的系统、等值SLM模型和华东现有负荷模型接于图3所示的负荷母线上进行仿真,得到系统的频率变化曲线图和城西220kV负荷站点有功功率曲线图分别如图4和图5所示。对比分析频率变化曲线和有功功率曲线,可以看到采用SLM模型与详细系统的仿真曲线的拟合效果明显好于采用现有负荷模型的负荷频率特性参数。因此与当前负荷模型的负荷频率特性参数相比,采用本发明实施例提供的SLM模型能够更好地描述负荷频率特性,且克服了当前负荷模型参数无法准确描述负荷频率特性的缺点,使故障后仿真计算中的系统特性更逼近真实的系统行为,提高了电力系统仿真计算的可信度,为电力系统的科学规划和安全稳定运行提供了有力保障。
基于同一发明构思,本发明实施例还提供了负荷频率特性参数的确定装置,这些设备解决问题的原理与负荷频率特性参数的确定装置方法相似,本发明实施例提供的负荷频率特性参数的确定装置主要包括第一计算模块、第二计算模块和确定模块,下面分别介绍上述三个模块的功能:
其中的第一计算模块,主要用于计算单个设备类型中静态负荷的有功功率和单个设备类型中动态负荷的有功功率;
其中的第二计算模块,主要用于根据单个设备类型中静态负荷的有功功率计算静态负荷的综合有功功率,并根据单个设备类型中动态负荷的有功功率计算动态负荷的综合有功功率;
其中的确定模块,主要用于根据静态负荷的综合有功功率确定静态负荷频率特性参数,并根据动态负荷的综合有功功率确定动态负荷频率特性参数;其中的静态负荷频率特性参数包括静态负荷的有功功率频率因子和无功功率频率因子,动态负荷频率特性参数包括动态负荷的机械转矩系数。
上述的第一计算模块按下式计算单个设备类型中静态负荷的有功功率:
PSi=Ni×NSi×PS
其中,Ni表示设备类型i的有功功率占负荷站点有功功率的百分比,Nsi为设备类型i中静态负荷的有功功率百分比,PS表示静态负荷的有功功率初值。
上述的第一计算模块进一步按下式计算单个设备类型中动态负荷的有功功率:
PRi=Ni×NMi×P
其中,NMi表示设备类型i中动态负荷的有功功率百分比,Nsi为设备类型i中静态负荷的有功功率百分比,P表示负荷站点的有功功率。
上述的第二计算模块按下式计算静态负荷的综合有功功率:
其中,PSa表示静态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
上述的第二计算模块进一步按下式计算动态负荷的综合有功功率:
其中,PRa表示动态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
上述的确定模块根据静态负荷的综合有功功率确定如下式的静态负荷的有功功率频率因子:
其中,LDP表示静态负荷的有功功率频率因子,Pfi表示设备类型i中静态负荷的有功功率-频率特性系数;
上述的确定模块进一步根据静态负荷的综合有功功率确定如下式的静态负荷的无功功率频率因子:
其中,LDQ表示静态负荷的无功功率频率因子,Qfi表示设备类型i中静态负荷的无功功率-频率特性系数。
上述的静态负荷的有功功率频率因子LDP和静态负荷的无功功率频率因子LDQ具体用于建立如下式的静态负荷模型:
其中,P1表示静态负荷的有功功率,PS表示静态负荷的有功功率初值,PZ表示恒定阻抗有功负荷比例,PI表示恒定电流有功负荷比例,PP表示恒定功率有功负荷比例;Q1表示静态负荷的无功功率,QS表示静态负荷的无功功率初值,QZ表示恒定阻抗无功负荷比例,QI表示恒定电流无功负荷比例,QP表示恒定功率无功负荷比例;V表示负荷母线实际电压幅值,V0为负荷母线电压幅值初值;Δf表示频率变化量。
上述的确定模块还需要根据动态负荷的综合有功功率确定如下式的动态负荷的机械转矩系数:
C=1-A-B
其中,A、B、C均表示动态负荷的机械转矩系数,Ai、Bi均表示设备类型i中动态负荷的机械转矩系数。
上述的动态负荷的机械转矩系数A、B、C具体应用于建立如下式的动态负荷模型:
TE=Pe/ωr
上述动态负荷取电动机,其中,ωr表示电动机的实时转速,Pe表示电动机的电磁功率,TE表示电动机的电磁转矩,TM表示电动机的,T'0表示电动机的时间常数,E'd表示电动机的d轴电势,E'q表示电动机的q轴电势,Id表示电动机的d轴电流,Iq表示电动机的q轴电流,ωb表示电动机的基准转速,X表示电动机的等值电抗,X'表示电动机的次暂态电抗。
为了描述的方便,以上所述装置的各部分以功能分为各种模块或单元分别描述。当然,在实施本申请时可以把各模块或单元的功能在同一个或多个软件或硬件中实现。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。
Claims (14)
1.一种负荷频率特性参数的确定方法,所述负荷包括静态负荷和动态负荷,其特征在于,所述确定方法包括:
计算单个设备类型中静态负荷的有功功率和单个设备类型中动态负荷的有功功率;
根据单个设备类型中静态负荷的有功功率计算静态负荷的综合有功功率,并根据单个设备类型中动态负荷的有功功率计算动态负荷的综合有功功率;
根据静态负荷的综合有功功率确定静态负荷频率特性参数,并根据动态负荷的综合有功功率确定动态负荷频率特性参数;
所述静态负荷频率特性参数包括静态负荷的有功功率频率因子和无功功率频率因子;
所述动态负荷频率特性参数包括动态负荷的机械转矩系数。
2.根据权利要求1所述的负荷频率特性参数的确定方法,其特征在于,所述单个设备类型中静态负荷的有功功率按下式计算:
PSi=Ni×NSi×PS
其中,Ni表示设备类型i的有功功率占负荷站点有功功率的百分比,Nsi为设备类型i中静态负荷的有功功率百分比,PS表示静态负荷的有功功率初值。
3.根据权利要求1所述的负荷频率特性参数的确定方法,其特征在于,所述单个设备类型中动态负荷的有功功率按下式计算:
PRi=Ni×NMi×P
其中,NMi表示设备类型i中动态负荷的有功功率百分比,Nsi为设备类型i中静态负荷的有功功率百分比,P表示负荷站点的有功功率。
4.根据权利要求2所述的负荷频率特性参数的确定方法,其特征在于,所述静态负荷的综合有功功率按下式计算:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
其中,PSa表示静态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
5.根据权利要求3所述的负荷频率特性参数的确定方法,其特征在于,所述动态负荷的综合有功功率按下式计算:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
其中,PRa表示动态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
6.根据权利要求4所述的负荷频率特性参数的确定方法,其特征在于,所述根据静态负荷的综合有功功率确定静态负荷频率特性参数包括:
根据静态负荷的综合有功功率确定如下式的静态负荷的有功功率频率因子:
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>D</mi>
<mi>P</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
</mrow>
其中,LDP表示静态负荷的有功功率频率因子,Pfi表示设备类型i中静态负荷的有功功率-频率特性系数;
根据静态负荷的综合有功功率确定如下式的静态负荷的无功功率频率因子:
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>D</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
</mrow>
其中,LDQ表示静态负荷的无功功率频率因子,Qfi表示设备类型i中静态负荷的无功功率-频率特性系数。
7.根据权利要求5所述的负荷频率特性参数的确定方法,其特征在于,所述根据动态负荷的综合有功功率确定动态负荷频率特性参数包括:
根据动态负荷的综合有功功率确定如下式的动态负荷的机械转矩系数:
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
</mrow>
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
</mrow>
C=1-A-B
其中,A、B、C均表示动态负荷的机械转矩系数,Ai、Bi均表示设备类型i中动态负荷的机械转矩系数。
8.一种负荷频率特性参数的确定装置,其特征在于,包括:
第一计算模块,用于计算单个设备类型中静态负荷的有功功率和单个设备类型中动态负荷的有功功率;
第二计算模块,用于根据单个设备类型中静态负荷的有功功率计算静态负荷的综合有功功率,并根据单个设备类型中动态负荷的有功功率计算动态负荷的综合有功功率;
确定模块,用于根据静态负荷的综合有功功率确定静态负荷频率特性参数,并根据动态负荷的综合有功功率确定动态负荷频率特性参数;
所述静态负荷频率特性参数包括静态负荷的有功功率频率因子和无功功率频率因子;
所述动态负荷频率特性参数包括动态负荷的机械转矩系数。
9.根据权利要求8所述的负荷频率特性参数的确定装置,其特征在于,所述第一计算模块按下式计算单个设备类型中静态负荷的有功功率:
PSi=Ni×NSi×PS
其中,Ni表示设备类型i的有功功率占负荷站点有功功率的百分比,Nsi为设备类型i中静态负荷的有功功率百分比,PS表示静态负荷的有功功率初值。
10.根据权利要求8所述的负荷频率特性参数的确定装置,其特征在于,所述第一计算模块按下式计算单个设备类型中动态负荷的有功功率:
PRi=Ni×NMi×P
其中,NMi表示设备类型i中动态负荷的有功功率百分比,Nsi为设备类型i中静态负荷的有功功率百分比,P表示负荷站点的有功功率。
11.根据权利要求9所述的负荷频率特性参数的确定装置,其特征在于,所述第二计算模块按下式计算静态负荷的综合有功功率:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
其中,PSa表示静态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
12.根据权利要求10所述的负荷频率特性参数的确定装置,其特征在于,所述第二计算模块按下式计算动态负荷的综合有功功率:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
其中,PRa表示动态负荷的综合有功功率,i=1,...,k,k为负荷中包含的设备类型个数。
13.根据权利要求11所述的负荷频率特性参数的确定装置,其特征在于,所述确定模块具体用于:
根据静态负荷的综合有功功率确定如下式的静态负荷的有功功率频率因子:
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>D</mi>
<mi>P</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
</mrow>
其中,LDP表示静态负荷的有功功率频率因子,Pfi表示设备类型i中静态负荷的有功功率-频率特性系数;
根据静态负荷的综合有功功率确定如下式的静态负荷的无功功率频率因子:
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>D</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>S</mi>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
</mrow>
其中,LDQ表示静态负荷的无功功率频率因子,Qfi表示设备类型i中静态负荷的无功功率-频率特性系数。
14.根据权利要求12所述的负荷频率特性参数的确定装置,其特征在于,所述确定模块具体用于:
根据动态负荷的综合有功功率确定如下式的动态负荷的机械转矩系数:
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
</mrow>
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>R</mi>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
</mrow>
C=1-A-B
其中,A、B、C均表示动态负荷的机械转矩系数,Ai、Bi均表示设备类型i中动态负荷的机械转矩系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710684435.6A CN107633110A (zh) | 2017-08-11 | 2017-08-11 | 一种负荷频率特性参数的确定方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710684435.6A CN107633110A (zh) | 2017-08-11 | 2017-08-11 | 一种负荷频率特性参数的确定方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107633110A true CN107633110A (zh) | 2018-01-26 |
Family
ID=61099471
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710684435.6A Pending CN107633110A (zh) | 2017-08-11 | 2017-08-11 | 一种负荷频率特性参数的确定方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107633110A (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105224812A (zh) * | 2015-10-21 | 2016-01-06 | 中国电力科学研究院 | 一种负荷模型中的静态负荷频率因子聚合方法 |
CN105429132A (zh) * | 2015-12-07 | 2016-03-23 | 中国电力科学研究院 | 一种电动机负荷模型的构建方法 |
CN105429131A (zh) * | 2015-12-07 | 2016-03-23 | 中国电力科学研究院 | 一种考虑负荷频率特性的负荷模型构建方法 |
CN105552917A (zh) * | 2015-12-07 | 2016-05-04 | 中国电力科学研究院 | 一种考虑负荷低电压释放特性的负荷模型构建方法 |
US20170033560A1 (en) * | 2015-07-31 | 2017-02-02 | Enphase Energy, Inc. | Method and apparatus for controlling a load shedding device |
-
2017
- 2017-08-11 CN CN201710684435.6A patent/CN107633110A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170033560A1 (en) * | 2015-07-31 | 2017-02-02 | Enphase Energy, Inc. | Method and apparatus for controlling a load shedding device |
CN105224812A (zh) * | 2015-10-21 | 2016-01-06 | 中国电力科学研究院 | 一种负荷模型中的静态负荷频率因子聚合方法 |
CN105429132A (zh) * | 2015-12-07 | 2016-03-23 | 中国电力科学研究院 | 一种电动机负荷模型的构建方法 |
CN105429131A (zh) * | 2015-12-07 | 2016-03-23 | 中国电力科学研究院 | 一种考虑负荷频率特性的负荷模型构建方法 |
CN105552917A (zh) * | 2015-12-07 | 2016-05-04 | 中国电力科学研究院 | 一种考虑负荷低电压释放特性的负荷模型构建方法 |
Non-Patent Citations (1)
Title |
---|
张孝 等: "考虑频率特性的负荷模型研究及应用", 《河海大学学报(自然科学版)》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | An oscillatory stability criterion based on the unified $ dq $-frame impedance network model for power systems with high-penetration renewables | |
CN102682358B (zh) | 一种评估风电并网规模与电网网架适应性的规划仿真方法 | |
CN104362622B (zh) | 一种带小水电机组地区电网的负荷模型和建模方法 | |
CN104201674B (zh) | 一种考虑负荷低电压释放特性的综合负荷模型的建模方法 | |
Xu et al. | Small-signal modeling and analysis of DC-link dynamics in type-IV wind turbine system | |
CN104992017B (zh) | 计及风电功率随机波动的广义负荷建模方法与装置 | |
CN107132772A (zh) | 交直流电网的实时仿真系统及方法 | |
CN107069800A (zh) | 一种双馈风电场短路电流模型的建立方法 | |
CN105429131B (zh) | 一种考虑负荷频率特性的负荷模型构建方法 | |
CN105429132A (zh) | 一种电动机负荷模型的构建方法 | |
Zhang et al. | Mechanism analysis of subsynchronous torsional interaction with PMSG-based WTs and LCC-HVDC | |
CN105224812B (zh) | 一种负荷模型中的静态负荷频率因子聚合方法 | |
Liu et al. | Impedance modeling of DFIG wind farms with various rotor speeds and frequency coupling | |
Trevisan et al. | Analysis of low frequency interactions of DFIG wind turbine systems in series compensated grids | |
CN115358079A (zh) | 风电场场站实时仿真模型的构建方法和阻抗特性评估方法 | |
CN111130135B (zh) | 一种适用于高比例新能源接入下的电力系统惯量计算方法 | |
Suvorov et al. | Comprehensive validation of transient stability calculations in electric power systems and hardware-software tool for its implementation | |
Li et al. | Assessing high-order harmonic resonance in locomotive-network based on the impedance method | |
Liu et al. | Highly-scalable transmission and distribution dynamic co-simulation with 10,000+ grid-following and grid-forming inverters | |
CN106655210A (zh) | 一种电网无功补偿方法 | |
Yin et al. | Identification of control parameters of the permanent magnetic synchronous generator using least square method | |
CN107633110A (zh) | 一种负荷频率特性参数的确定方法和装置 | |
CN101483343A (zh) | 以电网不对称运行和频率偏移信息辨识综合负荷等值马达方法 | |
CN111769571B (zh) | 基于光伏虚拟同步发电机的调频方法及系统 | |
Liu et al. | Construction of detailed energy function of double‐fed induction generator for small signal stability analysis of power systems with a high penetration of double‐fed induction generator |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180126 |
|
RJ01 | Rejection of invention patent application after publication |