[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN110750882A - 一种考虑频率约束的风电占比极限值解析计算方法 - Google Patents

一种考虑频率约束的风电占比极限值解析计算方法 Download PDF

Info

Publication number
CN110750882A
CN110750882A CN201910933156.8A CN201910933156A CN110750882A CN 110750882 A CN110750882 A CN 110750882A CN 201910933156 A CN201910933156 A CN 201910933156A CN 110750882 A CN110750882 A CN 110750882A
Authority
CN
China
Prior art keywords
frequency
power
wind power
limit value
wind
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
Application number
CN201910933156.8A
Other languages
English (en)
Other versions
CN110750882B (zh
Inventor
唐飞
刘涤尘
卿梦琪
刘福锁
都念纯
胡本溪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201910933156.8A priority Critical patent/CN110750882B/zh
Publication of CN110750882A publication Critical patent/CN110750882A/zh
Application granted granted Critical
Publication of CN110750882B publication Critical patent/CN110750882B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明涉及风电装机容量技术,尤其涉及一种考虑频率约束的风电占比极限值解析计算方法,对系统常规机组建模,得到反映各物理过程的详细频率响应模型;随后对系统进行等值处理,得到电力系统频率响应SFR单机等值降阶模型;根据风机接入的容量及系统原有容量计算风电占比值,进一步计算风机并网后的系统参数;代入原单机等值降阶模型中,得到风机并网后的电力系统频率响应模型,得出系统频率偏差标幺后的频域表达式,计算出系统频率稳态偏差的表达式。分别考虑系统稳态频率偏差约束边界及角频率变化率约束边界求出风电占比极限值,取两者的较小值,得到实际的风电功率占比极限值。该方法求解步骤简单,降低了计算复杂度,且具有一定的计算普适性。

Description

一种考虑频率约束的风电占比极限值解析计算方法
技术领域
本发明属于风电装机容量技术领域,尤其涉及一种考虑频率约束的风电占比极限值解析计算方法。
背景技术
由于风能资源和电力需求分布存在明显的不平衡性和区域差异,我国风能资源采取的远距离、大容量的大规模开发模式,导致系统频率控制能力下降,电力系统频率失稳风险增大。风电占比极限值是表征一定技术指标前提下接入系统的最大风机装机容量与系统总发电容量的百分比。当风电容量占比较小时,风电场对整个系统带来的影响并不十分明显;但风电容量占比较大时,风能的波动给电网带来的频率稳定、节点电压波动等问题,将严重危害电力系统的安全稳定运行。因此,有必要研究频率约束下系统接纳新能源的能力,确定系统在正常运行前提下可以接受的最大风电装机容量,即研究考虑频率约束下新能源在系统中的极限占比,进而较为准确地评估大规模风电并网对系统供电质量和可靠性的影响,确定合理的风能装机容量,有利于最大程度地防止过高比例的风电接入对系统带来的波动,确保电力系统的安全稳定运行。目前,风电占比极限值的计算方法主要有两种。一是考虑风速相关性和暂态稳定约束下,将风电接入容量最大作为优化目标进行求解,但该方法求解步骤较为复杂,且忽略了频率稳定约束;二是考虑频率稳定约束条件下,采用动态仿真法,通过预设风电占比值,再根据每次仿真的频率响应偏差不断修正,直到满足频率偏差约束为止,进而得到最终风电占比极限值。该方法虽适配性高,但仍是一种试凑法,结果不具备普遍性,且计算需反复操作。随着风电渗透率持续增加,现有的风电占比极限值研究方法已经难以适应,亟需提出考虑电力系统频率约束的风电占比极限值解析计算方法。
发明内容
本发明的目的是提供一种考虑了暂态频率变化率约束及稳态频率偏差约束的边界条件,对风电占比极限值进行具体的解析表达的方法。
为实现上述目的,本发明采用的技术方案是,一种考虑频率约束的风电占比极限值解析计算方法,包括以下步骤:
步骤1、系统常规机组建模:不考虑电力系统频率的时空分布特性时,将系统中电源和负荷集中等值处理,构建火电机组模型、水电机组模型、负荷阻尼调节模型、同步发电机转子动能模型,得到反映各物理过程的详细频率响应模型;
步骤2、简化等值:根据步骤1中构建的详细频率响应模型,将系统中所有机组等值为一台再热式火电机组,并忽略调速器动态响应过程及主进汽容积和汽室时间常数,得到电力系统频率响应SFR单机等值降阶模型,并表示出模型中等值机的惯性时间常数、调差系数、阻尼系数;
步骤3、计算风电接入后系统频率响应特性:根据风机接入的容量及系统原有容量计算风电占比值,基于风电占比值进一步计算风机并网后系统的等值机惯性时间常数、阻尼系数和调速器的调差系数;将各参数代入单机等值降阶模型中,得到风机并网后的电力系统频率响应模型,进而表示出系统频率偏差标幺后的频域表达式;
步骤4、根据系统频率稳态偏差约束计算风电占比极限值:通过终值定理,计算系统频率稳态偏差的表达式;系统稳态频率偏差的边界应为ξ,根据此表达式求解出频率稳态偏差约束下风电占比极限值α1
步骤5、根据系统频率变化率约束计算风电占比极限值:设在功率缺额产生的频率扰动初期,频率变化率达到最大值,计算出此时的频率变化率;考虑约束边界的角频率变化率η,反解出受频率变化率约束的风电占比极限值α2
步骤6、计算实际的风电功率占比极限值:同时考虑系统频率稳态偏差约束及系统频率变化率约束,取步骤4及步骤5中的较小值,得到实际的风电功率占比极限值。
在上述的考虑频率约束的风电占比极限值解析计算方法中,步骤1的实现包括以下步骤:
步骤1.1、对火电机组进行建模,汽轮机模型表示为:
Figure BDA0002220798600000031
式中,GTT(s)为火电机组汽轮机传递函数,ΔPT(s)为机械功率增量,Δδ为汽门开度增量,FHP、TRH、TCH分别为高压涡轮级HP级产生功率对汽轮机总功率占比、再热器时间常数、主进汽容积和汽室时间常数、汽门开度改变量;
调速器模型可以表示为:
Figure BDA0002220798600000032
式中,GnT(s)为火电机组调速器传递函数,Δδ为汽门开度增量,Δωs为系统频率增量,TGT、RT分别为汽轮机调速器的静态调差系数和调速器时间常数;
步骤1.2、对水电机组进行建模,水轮机模型可表示为:
Figure BDA0002220798600000033
式中,GTH(s)为水电机组水轮机传递函数,ΔPTH(s)为机械功率增量,ΔY为导水叶开度增量,Tw为额定负荷时水的起动时间;
调速器模型可表示为:
式中,GnH(s)为水电机组调速器传递函数,ΔY为导水叶开度增量,Δωs为系统频率增量,RH,TGH,Thy,R,RHT分别为水轮机调速器静态调差系数,调速器时间常数,复位时间,永久下降率,暂时下降率;
步骤1.3、负荷阻尼调节模型可表示为:
ΔPL=a1Δf+2a2f0Δf...nanf0 n-1Δf
=(a1+2a2f0...nanf0 n-1)Δf
=DΔf
式中,PL为负荷功率标幺值,f为系统频率标幺值,a0,a1,a2,an为与频率变化无关负荷,与频率变化线性相关负荷,与频率变化2次方相关负荷,与频率变化n次方相关负荷占总负荷的比例系数;D为负荷阻尼系数;
步骤1.4、同步发电机组转子动能模型可表示为:
Figure BDA0002220798600000042
其中,Wki为第i台同步发电机的旋转动能,ωi为同步发电机电角速度,ωn为系统额定角频率,Hi为第i台同步发电机的惯性时间常数,SNi为第i台同步发电机的总容量;
步骤1.5、将电力系统动态功率平衡方程变换到频域,并结合步骤1.1-1.4所建立的模型,得到系统详细频率响应模型,及频率偏差的表达式:
Figure BDA0002220798600000051
式中,GTTi(s)为第i台火电机组汽轮机传递函数,ΔPTNi(s)为第i台火电机组的机械功率增量,SNi为第i台同步发电机的总容量,PLN为负荷功率标幺值,Hi为第i台同步发电机的惯性时间常数,D为负荷阻尼系数,PTNi为第i台机组的机械功率,GTHi(s)为第i台水电机组水轮机传递函数,GnT(s)为第i台火电机组调速器传递函数,GnH(s)为水电机组调速器传递函数。
在上述的考虑频率约束的风电占比极限值解析计算方法中,步骤2所述电力系统频率响应SFR单机等值降阶模型为:
Figure BDA0002220798600000052
Figure BDA0002220798600000053
式中,Δω(s)为系统频率增量的频域表达,Km为机械功率增益因子,表示机械功率转换效率;Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数;Sj为第j台风机的额定容量,n和m分别为同步机台数和风机的台数;H为等值后发电机转子惯性时间常数,D为阻尼系数,R为调速器的调差系数,TR为汽轮机再热器时间常数,FH为高压锅炉输出功率所占比例,转子所承担加速功率,ΔPd(s)为扰动功率,ΔPL为扰动功率实际大小,Si、Hi和Ri分别为第i台同步机的额定容量、惯性时间常数和调速器调差系数,系统基准容量SB为系统各台发电机额定容量之和,即
Figure BDA0002220798600000061
在上述的考虑频率约束的风电占比极限值解析计算方法中,步骤3的实现包括以下步骤:
步骤3.1、根据风机接入的容量及系统原有容量计算风电占比值;
式中,α为电力系统中风电占比,Sj为第j台风机的额定容量,n和m分别同步机和风机的台数;
步骤3.2,基于此进一步计算风机并网后系统的等值机惯性时间常数、阻尼系数和调速器的调差系数:
式中,α为电力系统中风电容量占比,Sj为第j台风机的额定容量,n和m分别同步机和风机的台数,H’为风机并网后等值机惯性时间常数、D’为风机并网后阻尼系数、R’为风机并网后调速器的调差系数,Hi为各台同步发电机转子惯性时间常数,D为阻尼系数,Ri为各台同步发电机调速器的调差系数;
步骤3.3、将各参数代入单机等值降阶模型中,得到风机并网后的电力系统频率响应模型,进而表示出系统频率偏差标幺后的频域表达式:
Figure BDA0002220798600000071
式中,Δω(s)为系统频率增量的频域表达,α为电力系统中风电容量占比,Km为机械功率增益因子,Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数;Sj为第j台风机的额定容量,H为未加入风电前系统等值后发电机转子惯性时间常数,D为负荷阻尼系数,R为调速器的调差系数,TR为汽轮机再热器时间常数,FH为高压锅炉输出功率所占比例,ΔPd(s)为扰动功率,ΔPL(s)为扰动功率实际大小。
在上述的考虑频率约束的风电占比极限值解析计算方法中,步骤4的实现包括以下步骤:
步骤4.1、通过终值定理,计算系统频率稳态偏差的表达式:
Figure BDA0002220798600000072
其中,ΔPL为扰动功率实际大小,fN为系统频率额定值,α为电力系统中风电容量占比,D为负荷阻尼系数,Km为机械功率增益因子,Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数,R为系统中全部同步发电机调速器的等效调差系数,SB为系统基准容量,等于系统各台发电机额定容量之和,即
Figure BDA0002220798600000073
步骤4.2、根据国家标准相关规定,系统稳态频率偏差的边界应为ξ,即SFD≤ξ;根据此不等式求解出频率稳态偏差约束下风电占比极限值α1
Figure BDA0002220798600000074
其中,ΔPL为扰动功率实际大小,fN为系统频率额定值,α为电力系统中风电容量占比,D为负荷阻尼系数,Km为机械功率增益因子,Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数,R为系统中全部同步发电机调速器的等效调差系数,SB为系统基准容量,等于系统各台发电机额定容量之和,即
Figure BDA0002220798600000081
ξ为系统稳态频率偏差的边界。
在上述的考虑频率约束的风电占比极限值解析计算方法中,步骤5的实现包括以下步骤:
步骤5.1、设在功率缺额产生的频率扰动初期,频率变化率达到最大值;计算出此时的频率变化率:
Figure BDA0002220798600000082
步骤5.2、考虑约束边界的角频率变化率,边界应为η,反解出受频率变化率约束的风电占比极限值α2
Figure BDA0002220798600000083
其中,η为系统角频率变化率的边界值,H为未加入风电前系统等值后发电机转子惯性时间常数,ΔPL(s)为扰动功率实际大小,SB为系统基准容量,等于系统各台发电机额定容量之和,即
Figure BDA0002220798600000084
fN为系统频率额定值。
在上述的考虑频率约束的风电占比极限值解析计算方法中,步骤6所述实际的风电功率占比极限值:
α=min(α12)。
本发明的有益效果:1.同时考虑了暂态频率变化率约束及稳态频率偏差约束的边界条件,对风电占比极限值进行具体的解析表达,确定电力系统在正常运行前提下可以接受的最大风电装机容量,最大程度地防止过高比例的风电接入对电力系统带来的波动,确保电力系统的安全稳定运行。2.提出的方法对风电占比极限值进行了具体的解析表达,具有一定的计算普适性;3.求解步骤简单,降低了计算复杂度,能较为准确地评估大规模风电并网对系统供电质量和可靠性的影响,确定合理的风能装机容量,确保电力系统的安全稳定运行。
附图说明
图1是本发明一个实施例步骤1中反映系统各物理过程的详细频率响应模型;
图2是本发明一个实施例步骤2中忽略调速器动态响应过程及主进汽容积和汽室时间常数后得到的电力系统频率响应SFR单机等值降阶模型;
图3是本发明一个实施例步骤3中考虑风电并网后的电力系统频率响应模型。
具体实施方式
下面结合附图对本发明的实施方式进行详细描述。
本实施例通过以下技术方案来实现,一种考虑频率约束的风电占比极限值解析计算方法,首先对系统常规机组建模,得到反映各物理过程的详细频率响应模型;随后对系统进行等值处理,得到电力系统频率响应SFR单机等值降阶模型;根据风机接入的容量及系统原有容量计算风电占比值,进一步计算风机并网后的系统参数;代入原单机等值降阶模型中,得到风机并网后的电力系统频率响应模型,进而表示出系统频率偏差标幺后的频域表达式,计算出系统频率稳态偏差的表达式。随后分别考虑系统稳态频率偏差约束边界及角频率变化率约束边界,分别反解出受频率变化率约束的风电占比极限值,取两者的较小值,得到实际的风电功率占比极限值。
具体实施时,一种考虑频率约束的风电占比极限值解析计算方法,利用稳态频率偏差约束及暂态频率变化率约束对极限值占比进行解析计算,包括以下步骤:
S1.系统常规机组建模:不考虑电力系统频率的时空分布特性时,将系统中电源和负荷集中等值处理,构建火电机组模型、水电机组模型、负荷阻尼调节模型、同步发电机转子动能模型,得到反映各物理过程的详细频率响应模型;如图1所示。
S2.简化等值:根据S1中构建的详细频率响应模型,将系统中所有机组等值为一台再热式火电机组,并忽略调速器动态响应过程及主进汽容积和汽室时间常数,得到电力系统频率响应SFR单机等值降阶模型,并表示出模型中等值机的惯性时间常数、调差系数、阻尼系数;如图2所示。
S3.计算风电接入后系统频率响应特性:根据风机接入的容量及系统原有容量计算风电占比值,并基于此进一步计算风机并网后系统的等值机惯性时间常数、阻尼系数和调速器的调差系数。将各参数代入原单机等值降阶模型中,得到风机并网后的电力系统频率响应模型,进而表示出系统频率偏差标幺后的频域表达式;如图3所示。
S4.根据系统频率稳态偏差约束计算风电占比极限值:通过终值定理,计算系统频率稳态偏差的表达式。根据国家标准相关规定,系统稳态频率偏差的边界应为ξ,根据此不等式求解出频率稳态偏差约束下风电占比极限值α1
S5.根据系统频率变化率约束计算风电占比极限值:认为在功率缺额产生的频率扰动初期,频率变化率达到最大值。计算出此时的频率变化率,考虑约束边界的角频率变化率,反解出受频率变化率约束的风电占比极限值α2
S6.计算实际的风电功率占比极限值:同时考虑系统频率稳态偏差约束及系统频率变化率约束,取S4及S5中的较小值,得到实际的风电功率占比极限值。
实施例:
一种考虑频率稳定约束的风电占比极限值计算方法,包括以下步骤:
1,将系统中电源和负荷集中等值处理,构建火电机组模型、水电机组模型、负荷阻尼调节模型、同步发电机转子动能模型,得到反映各物理过程的详细频率响应模型。其中电源和负荷集中等值处理的结果如下。
1.1,对火电机组进行建模,汽轮机模型可表示为:
Figure BDA0002220798600000111
式中,GTT(s)为火电机组汽轮机传递函数,ΔPT(s)为机械功率增量,Δδ为汽门开度增量,FHP、TRH、TCH分别为高压涡轮级(HP级)产生功率对汽轮机总功率占比、再热器时间常数、主进汽容积和汽室时间常数、汽门开度改变量。
同时。调速器模型可以表示为:
Figure BDA0002220798600000112
其中,GnT(s)为火电机组调速器传递函数,Δδ为汽门开度增量,Δωs为系统频率增量,TGT、RT分别为汽轮机调速器的静态调差系数和调速器时间常数。
1.2,对水电机组进行建模,水轮机模型可表示为:
其中,GTH(s)为水电机组水轮机传递函数,ΔPTH(s)为机械功率增量,ΔY为导水叶开度增量,Tw为额定负荷时水的起动时间。
调速器模型可表示为:
Figure BDA0002220798600000121
其中,GnH(s)为水电机组调速器传递函数,ΔY为导水叶开度增量,Δωs为系统频率增量,RH,TGH,Thy,R,RHT分别为水轮机调速器静态调差系数,调速器时间常数,复位时间,永久下降率,暂时下降率。
1.3,负荷阻尼调节模型可表示为:
ΔPL=a1Δf+2a2f0Δf...nanf0 n-1Δf
=(a1+2a2f0...nanf0 n-1)Δf
=DΔf
其中,PL为负荷功率标幺值,f为系统频率标幺值,a0,a1,a2,an为与频率变化无关负荷,与频率变化线性相关负荷,与频率变化2次方相关负荷,与频率变化n次方相关负荷占总负荷的比例系数。D为负荷阻尼系数。
1.4,同步发电机组转子动能模型可表示为:
Figure BDA0002220798600000122
其中,Wki为第i台同步发电机的旋转动能,ωi为同步发电机电角速度,ωn为系统额定角频率,Hi为第i台同步发电机的惯性时间常数,SNi为第i台同步发电机的总容量。
1.5,将电力系统动态功率平衡方程变换到频域,并结合步骤1.1-1.4中建立的各项模型,得到系统详细频率响应模型,及频率偏差的表达式:
式中,GTTi(s)为第i台火电机组汽轮机传递函数,ΔPTNi(s)为第i台火电机组的机械功率增量,SNi为第i台同步发电机的总容量,PLN为负荷功率标幺值,Hi为第i台同步发电机的惯性时间常数,D为负荷阻尼系数,PTNi为第i台机组的机械功率,GTHi(s)为第i台水电机组水轮机传递函数,GnT(s)为第i台火电机组调速器传递函数,GnH(s)为水电机组调速器传递函数。
2,根据步骤1中构建的详细频率响应模型,将系统中所有机组等值为一台再热式火电机组,并忽略调速器动态响应过程及主进汽容积和汽室时间常数,得到电力系统频率响应SFR单机等值降阶模型。计算出系统频率偏差表达式,并表示出模型中等值机的惯性时间常数、调差系数、阻尼系数如下:
Figure BDA0002220798600000132
其中,Δω(s)为系统频率增量的频域表达,Km为机械功率增益因子,反映机械功率转换效率。Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数。Sj为第j台风机的额定容量,n和m分别为同步机台数和风机的台数。H为等值后发电机转子惯性时间常数,D为阻尼系数,R为调速器的调差系数,TR为汽轮机再热器时间常数,FH为高压锅炉输出功率所占比例,转子所承担加速功率,ΔPd(s)为扰动功率,ΔPL为扰动功率实际大小,Si、Hi和Ri分别为第i台同步机的额定容量、惯性时间常数和调速器调差系数,系统基准容量SB为系统各台发电机额定容量之和,即
Figure BDA0002220798600000141
3,构建风机并网后的电力系统频率响应模型,计算风电接入后系统频率响应特性,进而表示出系统频率偏差标幺后的频域表达式。具体步骤为:
3.1,根据风机接入的容量及系统原有容量计算风电占比值
Figure BDA0002220798600000142
式中,α为电力系统中风电占比,Sj为第j台风机的额定容量,n和m分别同步机和风机的台数。
3.2,基于此进一步计算风机并网后系统的等值机惯性时间常数、阻尼系数和调速器的调差系数:
Figure BDA0002220798600000143
式中,α为电力系统中风电容量占比,Sj为第j台风机的额定容量,n和m分别同步机和风机的台数,H’为风机并网后等值机惯性时间常数、D’为风机并网后阻尼系数、R’为风机并网后调速器的调差系数,Hi为各台同步发电机转子惯性时间常数,D为阻尼系数,Ri为各台同步发电机调速器的调差系数。
3.3,将各参数代入原单机等值降阶模型中,得到风机并网后的电力系统频率响应模型,进而表示出系统频率偏差标幺后的频域表达式:
式中,Δω(s)为系统频率增量的频域表达,α为电力系统中风电容量占比,Km为机械功率增益因子,Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数。Sj为第j台风机的额定容量,H为未加入风电前系统等值后发电机转子惯性时间常数,D为负荷阻尼系数,R为调速器的调差系数,TR为汽轮机再热器时间常数,FH为高压锅炉输出功率所占比例,ΔPd(s)为扰动功率,ΔPL(s)为扰动功率实际大小。
4,计算系统频率稳态偏差的表达式,并根据系统频率稳态偏差约束计算风电占比极限值α1。具体步骤为:
4.1,通过终值定理,计算系统频率稳态偏差的表达式:
Figure BDA0002220798600000152
其中,ΔPL为扰动功率实际大小,fN为系统频率额定值,α为电力系统中风电容量占比,D为负荷阻尼系数,Km为机械功率增益因子,Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数,R为系统中全部同步发电机调速器的等效调差系数,SB为系统基准容量,等于系统各台发电机额定容量之和,即
Figure BDA0002220798600000153
4.2,根据国家标准相关规定,系统稳态频率偏差的边界应为ξ,即SFD≤ξ(通常,中国国家标准《电能质量:电力系统频率允许偏差》中规定稳态频率偏差需在±0.2Hz以内),根据此不等式求解出频率稳态偏差约束下风电占比极限值α1
Figure BDA0002220798600000161
其中,ΔPL为扰动功率实际大小,fN为系统频率额定值,α为电力系统中风电容量占比,D为负荷阻尼系数,Km为机械功率增益因子,Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数,R为系统中全部同步发电机调速器的等效调差系数,SB为系统基准容量,等于系统各台发电机额定容量之和,即ξ为系统稳态频率偏差的边界。
5,计算系统频率变化率的最大值,根据系统频率变化率约束计算风电占比极限值α2。具体步骤为:
5.1,认为在功率缺额产生的频率扰动初期,频率变化率达到最大值。计算出此时的频率变化率:
Figure BDA0002220798600000163
5.2,考虑约束边界的角频率变化率,边界应为η(如欧洲Entso-e技术规定,频率变化率边界值为±0.005HZ/s),反解出受频率变化率约束的风电占比极限值α2
Figure BDA0002220798600000164
其中,η为系统角频率变化率的边界值,H为未加入风电前系统等值后发电机转子惯性时间常数,ΔPL(s)为扰动功率实际大小,SB为系统基准容量,等于系统各台发电机额定容量之和,即
Figure BDA0002220798600000171
fN为系统频率额定值。
6,同时考虑系统频率稳态偏差约束及系统频率变化率约束,取步骤4中计算得到的α1及步骤5中的计算求得的α2中的较小值,得到实际的风电功率占比极限值。
α=min(α12)
本实施例考虑了暂态频率变化率约束及稳态频率偏差约束的边界条件,对风电占比极限值进行具体的解析表达,确定系统在正常运行前提下可以接受的最大风电装机容量,最大程度地防止过高比例的风电接入对系统带来的波动,确保电力系统的安全稳定运行。求解步骤简单,降低了计算复杂度,且具有一定的计算普适性。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
虽然以上结合附图描述了本发明的具体实施方式,但是本领域普通技术人员应当理解,这些仅是举例说明,可以对这些实施方式做出多种变形或修改,而不背离本发明的原理和实质。本发明的范围仅由所附权利要求书限定。

Claims (7)

1.一种考虑频率约束的风电占比极限值解析计算方法,其特征是,包括以下步骤:
步骤1、系统常规机组建模:不考虑电力系统频率的时空分布特性时,将系统中电源和负荷集中等值处理,构建火电机组模型、水电机组模型、负荷阻尼调节模型、同步发电机转子动能模型,得到反映各物理过程的详细频率响应模型;
步骤2、简化等值:根据步骤1中构建的详细频率响应模型,将系统中所有机组等值为一台再热式火电机组,并忽略调速器动态响应过程及主进汽容积和汽室时间常数,得到电力系统频率响应SFR单机等值降阶模型,并表示出模型中等值机的惯性时间常数、调差系数、阻尼系数;
步骤3、计算风电接入后系统频率响应特性:根据风机接入的容量及系统原有容量计算风电占比值,基于风电占比值进一步计算风机并网后系统的等值机惯性时间常数、阻尼系数和调速器的调差系数;将各参数代入单机等值降阶模型中,得到风机并网后的电力系统频率响应模型,进而表示出系统频率偏差标幺后的频域表达式;
步骤4、根据系统频率稳态偏差约束计算风电占比极限值:通过终值定理,计算系统频率稳态偏差的表达式;系统稳态频率偏差的边界应为ξ,根据此表达式求解出频率稳态偏差约束下风电占比极限值α1
步骤5、根据系统频率变化率约束计算风电占比极限值:设在功率缺额产生的频率扰动初期,频率变化率达到最大值,计算出此时的频率变化率;考虑约束边界的角频率变化率η,反解出受频率变化率约束的风电占比极限值α2
步骤6、计算实际的风电功率占比极限值:同时考虑系统频率稳态偏差约束及系统频率变化率约束,取步骤4及步骤5中的较小值,得到实际的风电功率占比极限值。
2.如权利要求1所述的考虑频率约束的风电占比极限值解析计算方法,其特征是,步骤1的实现包括以下步骤:
步骤1.1、对火电机组进行建模,汽轮机模型表示为:
Figure FDA0002220798590000021
式中,GTT(s)为火电机组汽轮机传递函数,ΔPT(s)为机械功率增量,Δδ为汽门开度增量,FHP、TRH、TCH分别为高压涡轮级HP级产生功率对汽轮机总功率占比、再热器时间常数、主进汽容积和汽室时间常数、汽门开度改变量;
调速器模型可以表示为:
式中,GnT(s)为火电机组调速器传递函数,Δδ为汽门开度增量,Δωs为系统频率增量,TGT、RT分别为汽轮机调速器的静态调差系数和调速器时间常数;
步骤1.2、对水电机组进行建模,水轮机模型可表示为:
Figure FDA0002220798590000023
式中,GTH(s)为水电机组水轮机传递函数,ΔPTH(s)为机械功率增量,ΔY为导水叶开度增量,Tw为额定负荷时水的起动时间;
调速器模型可表示为:
Figure FDA0002220798590000024
式中,GnH(s)为水电机组调速器传递函数,ΔY为导水叶开度增量,Δωs为系统频率增量,RH,TGH,Thy,R,RHT分别为水轮机调速器静态调差系数,调速器时间常数,复位时间,永久下降率,暂时下降率;
步骤1.3、负荷阻尼调节模型可表示为:
ΔPL=a1Δf+2a2f0Δf...nanf0 n-1Δf
=(a1+2a2f0...nanf0 n-1)Δf
=DΔf
式中,PL为负荷功率标幺值,f为系统频率标幺值,a0,a1,a2,an为与频率变化无关负荷,与频率变化线性相关负荷,与频率变化2次方相关负荷,与频率变化n次方相关负荷占总负荷的比例系数;D为负荷阻尼系数;
步骤1.4、同步发电机组转子动能模型可表示为:
Figure FDA0002220798590000031
其中,Wki为第i台同步发电机的旋转动能,ωi为同步发电机电角速度,ωn为系统额定角频率,Hi为第i台同步发电机的惯性时间常数,SNi为第i台同步发电机的总容量;
步骤1.5、将电力系统动态功率平衡方程变换到频域,并结合步骤1.1-1.4所建立的模型,得到系统详细频率响应模型,及频率偏差的表达式:
Figure FDA0002220798590000032
式中,GTTi(s)为第i台火电机组汽轮机传递函数,ΔPTNi(s)为第i台火电机组的机械功率增量,SNi为第i台同步发电机的总容量,PLN为负荷功率标幺值,Hi为第i台同步发电机的惯性时间常数,D为负荷阻尼系数,PTNi为第i台机组的机械功率,GTHi(s)为第i台水电机组水轮机传递函数,GnT(s)为第i台火电机组调速器传递函数,GnH(s)为水电机组调速器传递函数。
3.如权利要求2所述的考虑频率约束的风电占比极限值解析计算方法,其特征是,步骤2所述电力系统频率响应SFR单机等值降阶模型为:
Figure FDA0002220798590000041
Figure FDA0002220798590000042
式中,Δω(s)为系统频率增量的频域表达,Km为机械功率增益因子,表示机械功率转换效率;Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数;Sj为第j台风机的额定容量,n和m分别为同步机台数和风机的台数;H为等值后发电机转子惯性时间常数,D为阻尼系数,R为调速器的调差系数,TR为汽轮机再热器时间常数,FH为高压锅炉输出功率所占比例,转子所承担加速功率,ΔPd(s)为扰动功率,ΔPL为扰动功率实际大小,Si、Hi和Ri分别为第i台同步机的额定容量、惯性时间常数和调速器调差系数,系统基准容量SB为系统各台发电机额定容量之和,即
Figure FDA0002220798590000043
4.如权利要求3所述的考虑频率约束的风电占比极限值解析计算方法,其特征是,步骤3的实现包括以下步骤:
步骤3.1、根据风机接入的容量及系统原有容量计算风电占比值;
Figure FDA0002220798590000044
式中,α为电力系统中风电占比,Sj为第j台风机的额定容量,n和m分别同步机和风机的台数;
步骤3.2,基于此进一步计算风机并网后系统的等值机惯性时间常数、阻尼系数和调速器的调差系数:
Figure FDA0002220798590000051
式中,α为电力系统中风电容量占比,Sj为第j台风机的额定容量,n和m分别同步机和风机的台数,H’为风机并网后等值机惯性时间常数、D’为风机并网后阻尼系数、R’为风机并网后调速器的调差系数,Hi为各台同步发电机转子惯性时间常数,D为阻尼系数,Ri为各台同步发电机调速器的调差系数;
步骤3.3、将各参数代入单机等值降阶模型中,得到风机并网后的电力系统频率响应模型,进而表示出系统频率偏差标幺后的频域表达式:
Figure FDA0002220798590000052
式中,Δω(s)为系统频率增量的频域表达,α为电力系统中风电容量占比,Km为机械功率增益因子,Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数;Sj为第j台风机的额定容量,H为未加入风电前系统等值后发电机转子惯性时间常数,D为负荷阻尼系数,R为调速器的调差系数,TR为汽轮机再热器时间常数,FH为高压锅炉输出功率所占比例,ΔPd(s)为扰动功率,ΔPL(s)为扰动功率实际大小。
5.如权利要求4所述的考虑频率约束的风电占比极限值解析计算方法,其特征是,步骤4的实现包括以下步骤:
步骤4.1、通过终值定理,计算系统频率稳态偏差的表达式:
Figure FDA0002220798590000061
其中,ΔPL为扰动功率实际大小,fN为系统频率额定值,α为电力系统中风电容量占比,D为负荷阻尼系数,Km为机械功率增益因子,Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数,R为系统中全部同步发电机调速器的等效调差系数,SB为系统基准容量,等于系统各台发电机额定容量之和,即
Figure FDA0002220798590000062
步骤4.2、根据国家标准相关规定,系统稳态频率偏差的边界应为ξ,即SFD≤|ξ|;根据此不等式求解出频率稳态偏差约束下风电占比极限值α1
Figure FDA0002220798590000063
其中,ΔPL为扰动功率实际大小,fN为系统频率额定值,α为电力系统中风电容量占比,D为负荷阻尼系数,Km为机械功率增益因子,Km=PF(1-fsr),其中PF为功率因数,fsr为备用系数,R为系统中全部同步发电机调速器的等效调差系数,SB为系统基准容量,等于系统各台发电机额定容量之和,即
Figure FDA0002220798590000064
ξ为系统稳态频率偏差的边界。
6.如权利要求5所述的考虑频率约束的风电占比极限值解析计算方法,其特征是,步骤5的实现包括以下步骤:
步骤5.1、设在功率缺额产生的频率扰动初期,频率变化率达到最大值;计算出此时的频率变化率:
Figure FDA0002220798590000071
步骤5.2、考虑约束边界的角频率变化率,边界应为η,反解出受频率变化率约束的风电占比极限值α2
Figure FDA0002220798590000072
其中,η为系统角频率变化率的边界值,H为未加入风电前系统等值后发电机转子惯性时间常数,ΔPL(s)为扰动功率实际大小,SB为系统基准容量,等于系统各台发电机额定容量之和,即
Figure FDA0002220798590000073
fN为系统频率额定值。
7.如权利要求6所述的考虑频率约束的风电占比极限值解析计算方法,其特征是,步骤6所述实际的风电功率占比极限值:
α=min(α12)。
CN201910933156.8A 2019-09-29 2019-09-29 一种考虑频率约束的风电占比极限值解析计算方法 Active CN110750882B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910933156.8A CN110750882B (zh) 2019-09-29 2019-09-29 一种考虑频率约束的风电占比极限值解析计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910933156.8A CN110750882B (zh) 2019-09-29 2019-09-29 一种考虑频率约束的风电占比极限值解析计算方法

Publications (2)

Publication Number Publication Date
CN110750882A true CN110750882A (zh) 2020-02-04
CN110750882B CN110750882B (zh) 2022-09-02

Family

ID=69277364

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910933156.8A Active CN110750882B (zh) 2019-09-29 2019-09-29 一种考虑频率约束的风电占比极限值解析计算方法

Country Status (1)

Country Link
CN (1) CN110750882B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113255485A (zh) * 2021-05-13 2021-08-13 云南电网有限责任公司 一种水力发电机组并网模式的识别方法及装置
CN113688511A (zh) * 2021-08-10 2021-11-23 国网江苏省电力有限公司经济技术研究院 一种电力系统频率稳定性的评估方法
CN113839416A (zh) * 2021-09-13 2021-12-24 国网湖南省电力有限公司 基于频率响应模型的风电最大接入比例确定方法及装置
CN113922389A (zh) * 2021-10-09 2022-01-11 湖南铁路科技职业技术学院 一种考虑频率约束的电力系统改进故障分析方法
CN113991647A (zh) * 2021-10-25 2022-01-28 国网辽宁省电力有限公司经济技术研究院 一种面向频率响应容量规划的电力系统随机生产模拟方法
CN114204611A (zh) * 2021-12-08 2022-03-18 华北电力大学 一种适用所有阻尼状态的频率响应解析计算方法
CN114243805A (zh) * 2021-12-20 2022-03-25 华北电力大学 一种考虑调速器限幅的同步机系统频率响应解析计算方法
CN115085290A (zh) * 2022-08-22 2022-09-20 中国电力科学研究院有限公司 一种新能源最大接入容量的确定方法、装置及介质
CN115189397A (zh) * 2022-07-07 2022-10-14 国网湖南省电力有限公司 电力系统的新能源接入容量计算方法
CN115986766A (zh) * 2022-10-28 2023-04-18 浙江大学 一种考虑调差系数灵敏度与频率约束的机组组合方法
CN117293860A (zh) * 2023-09-28 2023-12-26 暨南大学 电力系统的频率安全指标评估方法、系统及存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001178194A (ja) * 1999-12-16 2001-06-29 Okinawa Electric Power Co Ltd 風力発電装置の出力制御方法
CN103178554A (zh) * 2013-04-09 2013-06-26 国家电网公司 基于多重约束的大区电网间交换容量极限评估方法
CN106952180A (zh) * 2017-03-02 2017-07-14 河海大学 一种建立双馈分布式风电系统低阶频率响应模型的方法
CN107240918A (zh) * 2017-07-27 2017-10-10 江苏大学 一种风电接入的电力系统等值简化方法
CN107679769A (zh) * 2017-10-25 2018-02-09 东南大学 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法
CN108832658A (zh) * 2018-06-22 2018-11-16 三峡大学 一种考虑频率约束及风电调频的风电穿透功率极限计算方法
CN108964032A (zh) * 2018-07-13 2018-12-07 华北电力大学(保定) 一种建立含风电的系统频率响应模型的方法
CN109980660A (zh) * 2019-04-17 2019-07-05 陈建华 基于发电机改进频率响应特性的鲁棒区间风电调度方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001178194A (ja) * 1999-12-16 2001-06-29 Okinawa Electric Power Co Ltd 風力発電装置の出力制御方法
CN103178554A (zh) * 2013-04-09 2013-06-26 国家电网公司 基于多重约束的大区电网间交换容量极限评估方法
CN106952180A (zh) * 2017-03-02 2017-07-14 河海大学 一种建立双馈分布式风电系统低阶频率响应模型的方法
CN107240918A (zh) * 2017-07-27 2017-10-10 江苏大学 一种风电接入的电力系统等值简化方法
CN107679769A (zh) * 2017-10-25 2018-02-09 东南大学 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法
CN108832658A (zh) * 2018-06-22 2018-11-16 三峡大学 一种考虑频率约束及风电调频的风电穿透功率极限计算方法
CN108964032A (zh) * 2018-07-13 2018-12-07 华北电力大学(保定) 一种建立含风电的系统频率响应模型的方法
CN109980660A (zh) * 2019-04-17 2019-07-05 陈建华 基于发电机改进频率响应特性的鲁棒区间风电调度方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MENGJIAO QING,ET AL.: "A two-step splitting stragedy based on laplacian eigenmap algorothm for large power grid with power wind connected", 《2019 IEEE PES INNOVATIVE SMART GRID TECHNOLOGIES ASIA》 *
肖固城等: "一种可改善失步振荡模式的风电虚拟惯性控制策略", 《电力系统保护与控制》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113255485B (zh) * 2021-05-13 2022-05-17 云南电网有限责任公司 一种水力发电机组并网模式的识别方法及装置
CN113255485A (zh) * 2021-05-13 2021-08-13 云南电网有限责任公司 一种水力发电机组并网模式的识别方法及装置
CN113688511A (zh) * 2021-08-10 2021-11-23 国网江苏省电力有限公司经济技术研究院 一种电力系统频率稳定性的评估方法
CN113688511B (zh) * 2021-08-10 2023-10-13 国网江苏省电力有限公司经济技术研究院 一种电力系统频率稳定性的评估方法
CN113839416B (zh) * 2021-09-13 2023-10-27 国网湖南省电力有限公司 基于频率响应模型的风电最大接入比例确定方法及装置
CN113839416A (zh) * 2021-09-13 2021-12-24 国网湖南省电力有限公司 基于频率响应模型的风电最大接入比例确定方法及装置
CN113922389A (zh) * 2021-10-09 2022-01-11 湖南铁路科技职业技术学院 一种考虑频率约束的电力系统改进故障分析方法
CN113922389B (zh) * 2021-10-09 2023-06-20 湖南铁路科技职业技术学院 一种考虑频率约束的电力系统改进故障分析方法
CN113991647A (zh) * 2021-10-25 2022-01-28 国网辽宁省电力有限公司经济技术研究院 一种面向频率响应容量规划的电力系统随机生产模拟方法
CN114204611A (zh) * 2021-12-08 2022-03-18 华北电力大学 一种适用所有阻尼状态的频率响应解析计算方法
CN114243805A (zh) * 2021-12-20 2022-03-25 华北电力大学 一种考虑调速器限幅的同步机系统频率响应解析计算方法
CN114243805B (zh) * 2021-12-20 2023-01-10 华北电力大学 一种考虑调速器限幅的同步机系统频率响应解析计算方法
CN115189397A (zh) * 2022-07-07 2022-10-14 国网湖南省电力有限公司 电力系统的新能源接入容量计算方法
CN115085290A (zh) * 2022-08-22 2022-09-20 中国电力科学研究院有限公司 一种新能源最大接入容量的确定方法、装置及介质
CN115085290B (zh) * 2022-08-22 2022-10-28 中国电力科学研究院有限公司 一种新能源最大接入容量的确定方法、装置及介质
CN115986766B (zh) * 2022-10-28 2023-09-29 浙江大学 一种考虑调差系数灵敏度与频率约束的机组组合方法
CN115986766A (zh) * 2022-10-28 2023-04-18 浙江大学 一种考虑调差系数灵敏度与频率约束的机组组合方法
CN117293860A (zh) * 2023-09-28 2023-12-26 暨南大学 电力系统的频率安全指标评估方法、系统及存储介质
CN117293860B (zh) * 2023-09-28 2024-04-23 暨南大学 电力系统的频率安全指标评估方法、系统及存储介质

Also Published As

Publication number Publication date
CN110750882B (zh) 2022-09-02

Similar Documents

Publication Publication Date Title
CN110750882B (zh) 一种考虑频率约束的风电占比极限值解析计算方法
CN111245014B (zh) 一种基于虚拟惯量的电力系统控制方法
CN106773666A (zh) 一种针对压水堆一回路系统的模型参数获取方法
CN110829487A (zh) 一种电力系统的频率动态预测方法
CN109449937B (zh) 一种通过系统频率稳定约束确定电网新能源承载力的方法
CN109713661B (zh) 风电场接入对多机系统故障极限切除时间影响的分析方法
CN109599865B (zh) 一种确定电力系统中新能源占比的方法
CN105186502B (zh) 基于安全域的含双馈风机电力系统暂态稳定性分析方法
CN110649596B (zh) 一种考虑系统初始状态的频率全响应解析模型
Gao et al. A fast high-precision model of the doubly-fed pumped storage unit
Tan et al. Multi-time scale model reduction strategy of variable-speed pumped storage unit grid-connected system for small-signal oscillation stability analysis
CN110445194A (zh) 一种基于水、火电不同比例一次调频能力优化方法
CN109980660B (zh) 基于发电机改进频率响应特性的鲁棒区间风电调度方法
CN109787274A (zh) 一种虚拟同步控制方法及转子侧变频器控制器
Wang et al. Development of wind-energy modeling technology and standards
CN114665506A (zh) 一种大规模风电并网情况下电力系统惯性常数的测量方法
Liu et al. Coordinated Frequency Regulation Strategy of Pumped Storage Units and Battery Energy Storage System
CN103956767A (zh) 一种考虑尾流效应的风电场并网稳定性分析方法
CN113991647B (zh) 一种面向频率响应容量规划的电力系统随机生产模拟方法
CN109657380A (zh) 一种基于扩展卡尔曼滤波的双馈风电场动态等值方法
CN116093970A (zh) 计及转速保护的双馈风机一次调频模型预测控制方法
CN104037806A (zh) 一种风力发电机组基本模型和潮流问题的研究方法
Omran et al. Grid integration of a renewable energy system: modeling and analysis
Gu et al. Research on Equivalent Modeling of PMSG-based Wind Farms using Parameter Identification method
Hu et al. Frequency characteristic analysis of power system considering deep peak load regulation and renewable energy injection

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