CN114172202A - 一种基于需求侧资源主动响应的含风电互联电力系统负荷频率控制方法 - Google Patents

一种基于需求侧资源主动响应的含风电互联电力系统负荷频率控制方法 Download PDF

Info

Publication number
CN114172202A
CN114172202A CN202111519102.0A CN202111519102A CN114172202A CN 114172202 A CN114172202 A CN 114172202A CN 202111519102 A CN202111519102 A CN 202111519102A CN 114172202 A CN114172202 A CN 114172202A
Authority
CN
China
Prior art keywords
frequency
power
formula
air conditioner
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.)
Granted
Application number
CN202111519102.0A
Other languages
English (en)
Other versions
CN114172202B (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.)
Huaneng Clean Energy Research Institute
North China Electric Power University
Clean Energy Branch of Huaneng Zhejiang Energy Development Co Ltd
Original Assignee
Huaneng Clean Energy Research Institute
North China Electric Power University
Clean Energy Branch of Huaneng Zhejiang Energy Development Co Ltd
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 Huaneng Clean Energy Research Institute, North China Electric Power University, Clean Energy Branch of Huaneng Zhejiang Energy Development Co Ltd filed Critical Huaneng Clean Energy Research Institute
Priority to CN202111519102.0A priority Critical patent/CN114172202B/zh
Publication of CN114172202A publication Critical patent/CN114172202A/zh
Application granted granted Critical
Publication of CN114172202B publication Critical patent/CN114172202B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241The oscillation concerning frequency
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/26Arrangements for eliminating or reducing asymmetry in polyphase networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/28The renewable source being wind energy
    • 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
    • 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
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/50Arrangements for eliminating or reducing asymmetry in polyphase networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Air Conditioning Control Device (AREA)

Abstract

本发明公开了一种需求侧资源主动响应的含风电互联电力系统负荷频率控制方法,通过设置响应死区和温度恢复功率,建立变频空调负荷主动响应模式,并设计相适应的线性自抗扰控制策略,给出多区域互联电力系统的线性自抗扰控制器参数整定。本发明用于平抑风电随机波动性,提高电网接纳能力,利用需求侧资源的灵活可调性提升含风电互联系统的负荷频率控制品质,在高风电波动下能使频差稳定在0.2Hz的允许偏差内,同时使系统联络线功率波动影响最小化,为新能源电力系统提供了可靠的频率稳定性。

Description

一种基于需求侧资源主动响应的含风电互联电力系统负荷频 率控制方法
技术领域
本发明属于新能源发电负荷频率控制领域,涉及一种基于需求侧资源主动响应的自抗扰负荷频率控制设计方法,特别涉及一种基于需求侧资源主动响应的含风电互联电力系统负荷频率控制方法。
背景技术
近年来,环境问题日益严峻,为实现巴黎气候协议,将全球温升限制在2摄氏度以内,全球能源系统去碳化速度加快,以风电为代表的新能源发展迅速,然而,风电天然具有波动性和不确定性,以风电为代表的新能源电力规模化并网会使电力系统的能量平衡、频率稳定等面临更大挑战,挖掘并提升系统中灵活资源的主动响应能力是保障新能源电力系统安全稳定运行的重要手段。
电力系统容许频率波动的范围有限,根据GB/T15945-2008《电能质量电力系统频率允许偏差》,我国电力系统频率偏差允许值为0.2Hz,当系统容量较大时,偏差值可放宽到0.5Hz。通过设计需求侧资源的主动响应使风电场输出的有功功率波动满足并网要求,常用的需求侧资源有空调、冰箱、热水器和电动汽车等,对其经过集中或分散的控制手段,可以对电压、频率等系统信号做出响应,参与系统调节。
空调作为一种温控设备具备一定的储能能力。人体能够接受的温度舒适范围较大,对于短暂微小的温度变化敏感性较弱,因此空调负荷不仅可以在不影响用户舒适度的前提下快速响应系统调频,而且由于其体量大的优势(夏季高峰期时占尖峰负荷的30%~50%),具有巨大的调频潜力,进而据此设计空调类需求侧资源参与电力系统调频方法。
发明目的
本发明的目的是提出一种基于需求侧资源主动响应的负荷频率优化控制方法,以解决风电规模化并网造成的能量失衡,频率振荡的问题。
发明内容
本发明提供了一种基于需求侧资源主动响应的含风电互联电力系统负荷频率控制方法,包括以下步骤:
步骤1:确定含风电互联电力系统负荷频率控制模型和变频空调负荷模型;
步骤2:在频率偏差很小,无需进行频率响应的调频段,设置空调资源响应调频死区Δfd,并利用调频死区Δfd及时修正室内温度,使之恢复到初始设定,建立温度自恢复功率PT0
步骤3:对步骤2中超出死区Δfd的频段,结合空调资源自身功率限制以及用户舒适度,对参与调频响应的变频空调设置最大输出功率PAC,max和最小输出功率PAC,min约束;
步骤4:根据步骤1的空调负荷模型和步骤2-3的约束,设计空调入网响应功率PAC和对应于频率偏差的动态补偿功率ΔPC
步骤5:针对多区域互联的新能源互联电力系统负荷频率控制对象设计线性自抗扰控制模型,并整定参数;
优选地,步骤1中所述含风电互联电力系统的负荷频率控制模型包括调速器模型、非再热汽轮机模型、发电机与负载模型、风电机组的数学模型:
其中,所述调速器模型表示为表示为如式(1)-(2)所示:
Figure BDA0003407188090000021
Figure BDA0003407188090000022
式中,ΔPV为调速器的输入,Δu为控制器输入指令,Δf为频率偏差,ΔPG为调速器输出指令;
所述非再热汽轮机模型表示为如式(3)所示:
Figure BDA0003407188090000031
式中ΔPT为汽轮机输出功率变化,ΔPG为调速器输出指令;
所述发电机与负载模型表示为如式(4)所示:
Figure BDA0003407188090000032
式中,Kps=1/D,Tps=2H/fD,D是负载频率相关参数,表示为D=Pl/f;Pl是额定负载,H是惯性常数,f是额定频率;
所述风电机组的数学模型表示为如式(5)所示:
Figure BDA0003407188090000033
式中,ΔPw为风电机组的输出功率,A为风机有效扫风面积,Vw为风速,ρ为空气的密度,Cp为风功率转换系数,λ为叶尖速比,β为风机桨距角;
所述含风电的互联区域的互联电力系统的负荷频率响应表示为如式(6)所示:
Δf=GP(s)(ΔPT(s)+ΔPw(s)-ΔPD(s)) (6),
式中,ΔPD为负载变化量;
所述含风电的互联区域的互联电力系统的被控量被表示为如式(7)所示:
ACEi=ΔPtie,i+BiΔfi (7),
式中,ACEi是区域i的区域控制偏差,ΔPtie,i表示区域
Figure BDA0003407188090000034
的联络线交换功率,Bi表示区域i的频率偏差常数,Δfi表示区域i的频率偏差。
进一步优选地,步骤1中所述的变频空调数学模型表示为如式(8)-(9)所示:
Tin(t)=Tin(t-1)+ω[Tout(t)-Tin(t-1)]+μPAC(t) (8),
Figure BDA0003407188090000041
式中Tin(t)、Tout(t)分别表示t时刻空调的室内温度和室外温度;ω、μ分别为设备和周围环境特性系数,PAC(t)为t时刻空调的入网运行功率;f(t)表示空调t时刻的运行频率;fmax和fmin分别表示空调运行的最大、最小频率,空调在不停机运行时应满足fmin≤f(t)≤fmax;θmaxmin分别表示温差上下限;K为常系数。
优选地,步骤2中变频空调资源在t时刻的温度自恢复功率PT0被为如式(10)所示:
Figure BDA0003407188090000042
式中T0为初始设定温度,Tin(t)、Tout(t)分别表示t时刻空调的室内温度和室外温度,t-1时刻类推;t时刻在死区的空调输出功率为PT0,满足设定温度的同时给频率大偏差期间的空调出力留出了响应裕度,即在需求侧资源需要出力时能够主动响应功率补偿。
优选地,步骤3中对参与调频响应的变频空调设置最大输出功率和最小输出功率约束,表示为如式(11)所示:
PAC,min≤PAC≤PAC,max (11),式中PAC为空调入网功率,PAC,max,PAC,min分别为变频空调入网的最大,最小输出功率。
优选地,步骤4中根据温度自恢复模式和变频空调模型确定的空调资源入网响应与补偿功率,表示为如式(12)-(13)所示:
Figure BDA0003407188090000051
ΔPC=PAC(t)-PT0(t) (13),
式中PAC为空调入网功率,KAC为频率调节特性系数,Δf为频率偏差,ΔPC为频率偏差的动态补偿功率,PAC,max、PAC,min分别为变频空调入网的最大、最小输出功率,PT0为温度自恢复功率。
优选地,步骤5中,所述设计线性自抗扰控制模型指设计适应需求侧主动响应的线性自抗扰控制器,包括以下子步骤:
步骤S51、将一切不确定因素和误差归为总干扰进行估计,表示为如式(14)所示:
y(t)(p)=f(y(t),…,y(t)p-1,u(t),d(t))+bu(t) (14),
式中,y(t)(p)是p阶系统对象,f是系统的总扰动,LADRC通过设计扩张状态观测器LESO来估计系统的总扰动f,并通过反馈回路消除扰动对系统的影响;
步骤S52、令
Figure BDA0003407188090000053
zp+1=f,假设f可微,令
Figure BDA0003407188090000054
则所设计的扩张状态观测器LESO被表示为如式(15)所示:
Figure BDA0003407188090000052
式中各参数分别表示为如式(16)-(18)所示:
z=[z1 z2 z3…zp zp+1] (16),
Figure BDA0003407188090000061
Figure BDA0003407188090000062
进一步优选地,设计一个全阶的Luenberger状态观测器,表示为如式(19)所示:
Figure BDA0003407188090000063
式中,
Figure BDA0003407188090000064
表示为z的估计值;L0是观测器的增益向量,结构表示为如式(20)所示:
L0=[β1 β2 … βp βp+1]T (20)。
进一步优选地,将扩张状态观测器LESO的带宽与特征值之间联系起来,将扩张状态观测器LESO极点配置在-ωo处,即表示为如式(21)所示:
|sI-(Ac-L0Cc)|=sp+11sp+…+βp+1=(s+ωo)p+1 (21),
式中,通过整定ω0的数值得到LESO的增益L0;当ω0的值增大时,扩张状态观测器LESO的速度会加快,同时对噪声的敏感度提高;
当Ac-L0Cc渐进稳定时,
Figure BDA0003407188090000065
将近似等于y(t)和它的各阶微分,
Figure BDA0003407188090000066
为系统总扰动
Figure BDA0003407188090000067
的估计值;
所述自抗扰的控制率表示为如式(22)-(23所示:
Figure BDA0003407188090000068
Figure BDA0003407188090000071
式中,
Figure BDA0003407188090000072
是估计的总扰动,
Figure BDA0003407188090000073
被控对象就被视为一个p阶的积分对象,表示为如式(24)-(25所示:
y(t)(p)≈u0(t) (24),
Figure BDA0003407188090000074
式中
Figure BDA0003407188090000075
近似等于
Figure BDA0003407188090000076
所述自抗扰的控制率最进一步被表示为如式(26)所示:
Figure BDA0003407188090000077
式中,K0=[k1 k2 k3 … 1]/b,将反馈控制器参数通过带宽法全部配置在-ωc处,示为如式(27)所示:
sp+kpsp-1+…+k1=(s+ωc)p (27),
式中,只需调节参数ωc,即可得到控制器的增益参数K0
附图说明
图1为本发明的整体思想的方法流程示意图。
图2为两区域互联电力系统中风电场原始输出功率。
图3为某时刻变频空调频率响应特性。
图4为本发明所述需求侧响应参与的两区域电力系统线性自抗扰负荷频率控制模型。
图5为风电波动下本发明中所述方法与传统方法的频率对比仿真图。
图6为负载扰动下本发明中所述方法与传统方法的联络线功率对比仿真图。
具体实施方式
本发明提出一种基于需求侧资源主动响应的负荷频率优化控制方法,下面结合附图和具体实施例对本发明做进一步说明。
图1为需求侧响应参与的负荷频率控制方法流程图。此负荷频率控制方法包括变频空调负荷的主动响应模式和相适应的线性自抗扰控制策略两部分,其具体实施步骤为:
1)确定以含风电两区域互联的电力系统为实施对象,具体模型见图4。
两区域的调速器模型为:
Figure BDA0003407188090000081
Figure BDA0003407188090000082
式中ΔPV为调速器的输入,Δu为控制器输入指令,Δf为频率偏差,ΔPG为调速器输出指令。
两区域均采用的非再热汽轮机模型为:
Figure BDA0003407188090000083
式中ΔPT为汽轮机输出功率变化,ΔPG为调速器输出指令。
发电机与负载模型为:
Figure BDA0003407188090000084
式中其中Kps=1/D,Tps=2H/fD。D是负载频率相关参数,D=Pl/f。Pl是额定负载,H是惯性常数,f是额定频率。
风电机组的数学模型为:
Figure BDA0003407188090000085
式中:ΔPw为风电输出功率;A为风机有效扫风面积;Vw为风速;ρ为空气的密度;Cp为风功率转换系数;λ为叶尖速比;β为风机桨距角。
2)确定变频空调负荷数学模型为:
Tin(t)=Tin(t-1)+ω[Tout(t)-Tin(t-1)]+μPAC(t)
Figure BDA0003407188090000091
式中Tin(t),Tout(t)分别表示为t时刻空调的室内温度和室外温度;ω,μ为设备和周围环境特性系数,PAC(t)为t时刻空调的入网运行功率。f(t)表示空调t时刻的运行频率;fmax和fmin分别表示空调运行的最大、最小频率,空调在不停机运行时应满足fmin≤f(t)≤fmax;θmax,θmin分别表示温差上下限;K为常系数。
3)频率偏差很小,无需进行频率响应的调频段,设置空调资源响应调频死区Δfd=0.05Hz;并利用死区及时修正室内温度,使之恢复到初始设定,建立温度自恢复功率PT0
变频空调资源在t时刻的温度自恢复功率PT0为:
Figure BDA0003407188090000092
式中T0为初始设定温度。此时在死区的空调输出功率为PT0,满足设定温度的同时给频率大偏差期间的空调出力留出了响应裕度,即在需求侧资源需要出力时可以主动响应功率补偿。
4)对参与调频响应的变频空调设置最大输出功率和最小输出功率约束:
PAC,min≤PAC≤PAC,max
式中PAC为空调入网功率,PAC,max,PAC,min分别为变频空调入网的最大,最小输出功率。
5)根据温度自恢复模式和变频空调模型确定的空调资源入网响应与补偿功率,并得出t时刻需求侧频率响应特性,见图3所示。
Figure BDA0003407188090000101
ΔPC=PAC(t)-PT0(t)
式中PAC为空调入网功率,KAC为频率调节特性系数,Δf为频率偏差,ΔPC为频率偏差的动态补偿功率。
6)进一步地,基于需求侧空调主动响应的负荷频率响应为:
Δf=GP(s)(ΔPT(s)+ΔPc(s)-ΔPw(s)-ΔPD(s))
式中ΔPD为负载变化;ΔPC为需求侧动态补偿功率;ΔPT为汽轮机输出功率变化;ΔPw为风电振荡功率,见图2所示。
进一步地,互联电力系统的被控量为:
ACEi=ΔPtie,i+BiΔfi
式中,ACEi是区域i的区域控制偏差,ΔPtie,i表示区域
Figure BDA0003407188090000102
的联络线交换功率,Bi表示区域i的频率偏差常数,Δfi表示区域i的频率偏差。
7)将一切不确定因素和误差归为总干扰进行估计进行适应需求侧主动响应的线性自抗扰控制器设计。
y(t)(p)=f(y(t),…,y(t)p-1,u(t),d(t))+bu(t)
式中,f是系统的总扰动,LADRC通过设计扩张状态观测器(LESO)来估计系统的总扰动f,并通过反馈回路消除扰动对系统的影响。
进一步地,令
Figure BDA0003407188090000103
zp+1=f,假设f可微,令
Figure BDA0003407188090000104
那么(13)就化为如下扩张状态观测器形式:
Figure BDA0003407188090000111
式中,z=[z1 z2 z3…zp zp+1],
Figure BDA0003407188090000112
Figure BDA0003407188090000113
Cc=[1 0 … 0 0]1×(p+1)
进一步地,一个全阶的Luenberger状态观测器可以设计为下面形式:
Figure BDA0003407188090000114
式中,
Figure BDA0003407188090000115
表示为z的估计值,L0是观测器的增益向量,结构如下:
L0=[β1β2…βpβp+1]T
8)对两区域电力系统的负荷频率控制采用一种分散整定方法。具体包括以下子步骤:
步骤S1:先忽略联络线交换功率,分别对单个区域设计LADRC控制器,见图4所示中,每个区域机组模型的开环传递函数为:
Figure BDA0003407188090000116
步骤S2:由传递函数得出其相对阶次p=3,对于两区域电力系统来说,设计三阶LADRC对系统进行控制。用带宽法将观测器LESO的带宽与特征值之间联系起来,将观测器极点配置在-ωo处:
|sI-(Ac-L0Cc)|=sp+11sp+…+βp+1=(s+ωo)p+1
通过整定ω0的数值就可得到LESO的增益L0。当ω0的值增大时,观测器的速度会加快,同时也会造成LESO对噪声的敏感度提高。对于需求侧资源的主动响应,选取高ω0=100,提高对空调负荷出力的适应度。
步骤S3:对于3阶区域系统(p=3)的实施对象,且Ac-L0Cc渐进稳定时,
Figure BDA0003407188090000121
将近似等于y(t)和它的各阶微分,
Figure BDA0003407188090000122
为系统总扰动f的估计值,因此扰动的估计值在控制系统中能够快速得到并消除它对系统的影响。自抗扰的控制率如下:
Figure BDA0003407188090000123
Figure BDA0003407188090000124
式中,当LESO参数选取合适,
Figure BDA0003407188090000125
为估计的总扰动;则
Figure BDA0003407188090000126
那被控对象就可以看成一个3阶的积分对象:
y(t)(3)≈u0(t),
步骤S4:对上述对象采用如下控制规律:,
Figure BDA0003407188090000127
因为
Figure BDA0003407188090000128
近似等于y(t),
Figure BDA0003407188090000129
那么控制率最终可以表示为:
Figure BDA00034071880900001210
式中,K0=[k1 k2 k31]/b。
步骤S5:多区域电力系统设计局部负荷频率控制器时需要将每个区域的对象都乘以该区域的频率偏差系数,本发明相应的控制量增益为:
Figure BDA00034071880900001211
步骤S6:将反馈控制器参数通过带宽法全部配置在-ωc处,对于图4中三阶负荷频率控制系统:
s3+k3s2+k2s1+k1=(s+ωc)3
只需调节参数ωc,就可以得到控制器的增益参数K0。改变ωc值,对需求侧资源的响应速度和系统频率的稳定性都有影响,为了充分激发需求侧资源的主动响应特性,选取最佳ωc=3。
9)以优化的工程整定PID作为策略1,加入本发明的需求侧主动响应的PID作为策略2,本发明的全套基于需求侧主动响应的自抗扰控制作为策略3,在图4所示的两区域互联的电力系统验证优化控制效果如图5和图6所示。
本发明所设计的LADRC是一种不依赖系统模型的控制结构,仅需知道系统的相对阶P和相应的增益b,就可对系统实现控制。另外,LADRC参数的整定只有ωc和ωo两个,整定方便,设计简单,便于工程实现,并且针对相对独立的空调负荷,LADRC可以充分激发需求侧资源的主动响应特性。
本发明的有益效果在于:
挖掘了电力系统中的灵活可调度的需求侧资源,并匹配先进的调控方法实现多能源互补,设计的基于需求侧资源主动响应的LADRC策略可用于补偿新能源电力系统中经典控制器的不足,从而平抑新能源电力随机波动性,提高电网接纳能力。该方法在降低频率跌落幅度上效果更好,且大大缩短了频率恢复时间,有效地控制了联络线功率的波动,提升了风电互联系统的负荷频率控制品质,为新能源电力系统提供了可靠的频率稳定性。
综上所述,用基于需求侧资源的灵活可调性,建立变频空调负荷主动响应模式,并用适应需求侧资源主动响应的LADRC策略,进行多区域互联电力系统的LADRC参数整定,达到了平抑风功率波动的要求,提升含风电互联系统的负荷频率控制品质。
此实施案例仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (9)

1.一种基于需求侧资源主动响应的含风电互联电力系统负荷频率控制方法,其特征在于,包括以下步骤:
步骤1:确定含风电互联电力系统负荷频率控制模型和变频空调负荷模型;
步骤2:在频率偏差很小,无需进行频率响应的调频段,设置空调资源响应调频死区Δfd,并利用调频死区Δfd及时修正室内温度,使之恢复到初始设定,建立温度自恢复功率PT0
步骤3:对步骤2中超出死区Δfd的频段,结合空调资源自身功率限制以及用户舒适度,对参与调频响应的变频空调设置最大输出功率PAC,max和最小输出功率PAC,min约束;
步骤4:根据步骤1的空调负荷模型和步骤2-3的约束,设计空调入网响应功率PAC和对应于频率偏差的动态补偿功率ΔPC
步骤5:针对多区域互联的新能源互联电力系统负荷频率控制对象设计线性自抗扰控制模型,并整定参数。
2.根据权利要求1所述的含风电互联电力系统负荷频率控制方法,其特征在于,步骤1中所述确定含风电互联电力系统负荷频率控制模型包括调速器模型、非再热汽轮机模型、发电机与负载模型、风电机组的数学模型:
其中,所述调速器模型表示为表示为如式(1)-(2)所示:
Figure FDA0003407188080000011
Figure FDA0003407188080000021
式中,ΔPV为调速器的输入,Δu为控制器输入指令,Δf为频率偏差,ΔPG为调速器输出指令;
所述非再热汽轮机模型表示为如式(3)所示:
Figure FDA0003407188080000022
式中ΔPT为汽轮机输出功率变化,ΔPG为调速器输出指令;
所述发电机与负载模型表示为如式(4)所示:
Figure FDA0003407188080000023
式中,Kps=1/D,Tps=2H/fD,D是负载频率相关参数,表示为D=Pl/f;Pl是额定负载,H是惯性常数,f是额定频率;
所述风电机组的数学模型表示为如式(5)所示:
Figure FDA0003407188080000024
式中,ΔPw为风电机组的输出功率,A为风机有效扫风面积,Vw为风速,ρ为空气的密度,Cp为风功率转换系数,λ为叶尖速比,β为风机桨距角;
所述含风电互联电力系统的负荷频率响应表示为如式(6)所示:
Δf=GP(s)(ΔPT(s)+ΔPw(s)-ΔPD(s)) (6),
式中,ΔPD为负载变化量;
所述含风电互联电力系统的被控量被表示为如式(7)所示:
ACEi=ΔPtie,i+BiΔfi (7),
式中,ACEi是区域i的区域控制偏差,ΔPtie,i表示区域
Figure FDA0003407188080000025
的联络线交换功率,Bi表示区域i的频率偏差常数,Δfi表示区域i的频率偏差。
3.根据权利要求2所述的含风电互联电力系统负荷频率控制方法,其特征在于,步骤1中所述的变频空调数学模型表示为如式(8)-(9)所示:
Tin(t)=Tin(t-1)+ω[Tout(t)-Tin(t-1)]+μPAC(t) (8),
Figure FDA0003407188080000031
式中Tin(t)、Tout(t)分别表示t时刻空调的室内温度和室外温度;ω、μ分别为设备和周围环境特性系数,PAC(t)为t时刻空调的入网运行功率;f(t)表示空调t时刻的运行频率;fmax和fmin分别表示空调运行的最大、最小频率,空调在不停机运行时应满足fmin≤f(t)≤fmax
θmaxmin分别表示温差上下限;K为常系数。
4.根据权利要求1所述的含风电互联电力系统负荷频率控制方法,其特征在于,步骤2中变频空调资源在t时刻的温度自恢复功率PT0被为如式(10)所示:
Figure FDA0003407188080000032
式中T0为初始设定温度,Tin(t)、Tout(t)分别表示t时刻空调的室内温度和室外温度,t-1时刻类推;t时刻在死区的空调输出功率为PT0,满足设定温度的同时给频率大偏差期间的空调出力留出了响应裕度,即在需求侧资源需要出力时能够主动响应功率补偿。
5.根据权利要求1所述的含风电互联电力系统负荷频率控制方法,其特征在于,步骤3中对参与调频响应的变频空调设置最大输出功率和最小输出功率约束,表示为如式(11)所示:
PAC,min≤PAC≤PAC,max (11),
式中PAC为空调入网功率,PAC,max,PAC,min分别为变频空调入网的最大,最小输出功率。
6.根据权利要求1所述的含风电互联电力系统负荷频率控制方法,其特征在于,步骤4中根据温度自恢复模式和变频空调模型确定的空调资源入网响应与补偿功率,表示为如式(12)-(13)所示:
Figure FDA0003407188080000041
ΔPC=PAC(t)-PT0(t) (13),
式中PAC为空调入网功率,KAC为频率调节特性系数,Δf为频率偏差,ΔPC为频率偏差的动态补偿功率,PAC,max、PAC,min分别为变频空调入网的最大、最小输出功率,PT0为温度自恢复功率。
7.根据权利要求1所述的含风电互联电力系统负荷频率控制方法,其特征在于,步骤5中,所述设计线性自抗扰控制模型指设计适应需求侧主动响应的线性自抗扰控制器,包括以下子步骤:
步骤S51、将一切不确定因素和误差归为总干扰进行估计,表示为如式(14)所示:
y(t)(p)=f(y(t),…,y(t)p-1,u(t),d(t))+bu(t) (14),
式中,y(t)(p)是p阶系统对象,f是系统的总扰动,LADRC通过设计扩张状态观测器LESO来估计系统的总扰动f,并通过反馈回路消除扰动对系统的影响;
步骤S52、令
Figure FDA0003407188080000057
zp+1=f,假设f可微,令
Figure FDA0003407188080000056
则所设计的扩张状态观测器LESO被表示为如式(15)所示:
Figure FDA0003407188080000051
式中各参数分别表示为如式(16)-(18)所示:
z=[z1 z2 z3…zp zp+1] (16),
Figure FDA0003407188080000052
Figure FDA0003407188080000053
8.根据权利要求7所述的含风电互联电力系统负荷频率控制方法,其特征在于,进一步设计一个全阶的Luenberger扩张状态观测器,表示为如式(19)所示:
Figure FDA0003407188080000054
式中,
Figure FDA0003407188080000055
表示为z的估计值;L0是观测器的增益向量,结构表示为如式(20)所示:
L0=[β1 β2…βp βp+1]T (20)。
9.根据权利要求8所述的含风电互联电力系统负荷频率控制方法,其特征在于,将扩张状态观测器LESO的带宽与特征值之间联系起来,将扩张状态观测器LESO极点配置在-ωo处,即表示为如式(21) 所示:
|sI-(Ac-L0Cc)|=sp+11sp+…+βp+1=(s+ωo)p+1 (21),
式中,通过整定ω0的数值得到LESO的增益L0;当ω0的值增大时,扩张状态观测器LESO的速度会加快,同时对噪声的敏感度提高;
当Ac-L0Cc渐进稳定时,
Figure FDA0003407188080000061
将近似等于y(t)和它的各阶微分,
Figure FDA0003407188080000062
为系统总扰动
Figure FDA0003407188080000063
的估计值;
所述自抗扰的控制率表示为如式(22)-(23所示:
Figure FDA0003407188080000064
Figure FDA0003407188080000065
式中,
Figure FDA0003407188080000066
是估计的总扰动,
Figure FDA0003407188080000067
被控对象就被视为一个p阶的积分对象,表示为如式(24)-(25所示:
y(t)(p)≈u0(t) (24),
Figure FDA0003407188080000068
式中
Figure FDA0003407188080000069
近似等于
Figure FDA00034071880800000610
所述自抗扰的控制率最进一步被表示为如式(26)所示:
Figure FDA00034071880800000611
式中,K0=[k1 k2 k3…1]/b,将反馈控制器参数通过带宽法全部配置在-ωc处,示为如式(27)所示:
sp+kpsp-1+…+k1=(s+ωc)p (27),
式中,只需调节参数ωc,即可得到控制器的增益参数K0
CN202111519102.0A 2021-12-13 2021-12-13 一种基于需求侧资源主动响应的含风电互联电力系统负荷频率控制方法 Active CN114172202B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111519102.0A CN114172202B (zh) 2021-12-13 2021-12-13 一种基于需求侧资源主动响应的含风电互联电力系统负荷频率控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111519102.0A CN114172202B (zh) 2021-12-13 2021-12-13 一种基于需求侧资源主动响应的含风电互联电力系统负荷频率控制方法

Publications (2)

Publication Number Publication Date
CN114172202A true CN114172202A (zh) 2022-03-11
CN114172202B CN114172202B (zh) 2023-11-28

Family

ID=80485999

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111519102.0A Active CN114172202B (zh) 2021-12-13 2021-12-13 一种基于需求侧资源主动响应的含风电互联电力系统负荷频率控制方法

Country Status (1)

Country Link
CN (1) CN114172202B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114935205A (zh) * 2022-07-27 2022-08-23 国网湖北省电力有限公司营销服务中心(计量中心) 多变频空调的双模型优化的协同控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110048437A (zh) * 2019-05-06 2019-07-23 昆明理工大学 一种基于线性自抗扰的孤岛交流微电网二次频率控制方法
CN111654047A (zh) * 2020-06-12 2020-09-11 昆明理工大学 一种基于自抗扰的抽水蓄能与电化学储能联合参与电网负荷频率控制方法
CN112531792A (zh) * 2020-12-03 2021-03-19 江苏方天电力技术有限公司 一种包含储能资源互联电力系统频率控制方法和系统
CN112994024A (zh) * 2021-02-02 2021-06-18 中南大学 一种改进型温控负荷参与的负荷频率控制方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110048437A (zh) * 2019-05-06 2019-07-23 昆明理工大学 一种基于线性自抗扰的孤岛交流微电网二次频率控制方法
CN111654047A (zh) * 2020-06-12 2020-09-11 昆明理工大学 一种基于自抗扰的抽水蓄能与电化学储能联合参与电网负荷频率控制方法
CN112531792A (zh) * 2020-12-03 2021-03-19 江苏方天电力技术有限公司 一种包含储能资源互联电力系统频率控制方法和系统
CN112994024A (zh) * 2021-02-02 2021-06-18 中南大学 一种改进型温控负荷参与的负荷频率控制方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李峰;余为杰;张周胜;: "含风电的多域互联电力系统负荷频率控制方法", 水电能源科学, no. 01 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114935205A (zh) * 2022-07-27 2022-08-23 国网湖北省电力有限公司营销服务中心(计量中心) 多变频空调的双模型优化的协同控制方法

Also Published As

Publication number Publication date
CN114172202B (zh) 2023-11-28

Similar Documents

Publication Publication Date Title
CN111900742B (zh) 一种风储系统基于双层协同控制的调频方法
CN107732977B (zh) 一种基于需求响应的agc实时调度方法
Zhang et al. Deloading power coordinated distribution method for frequency regulation by wind farms considering wind speed differences
CN107482649B (zh) 一种基于分频控制的两域互联系统负荷频率控制方法
CN110492498B (zh) 一种基于双层控制的温控负荷参与电力系统一次调频方法
CN111030141A (zh) 一种基于一致性算法的源-荷协同分布式优化调控方法
CN109103927B (zh) 调速系统pid控制器参数整定方法
CN109659961B (zh) 一种基于分频控制的动态电力系统负荷频率协调方法
CN108462212B (zh) 一种新能源电力系统在多源多调控域运行方式下控制方法
CN113757931B (zh) 一种空调控制方法及系统
CN112491070B (zh) 一种储能自适应阻尼vsg控制方法
CN104578181B (zh) 一种在限电弃风工况下双馈风电机组有功出力的控制方法
WO2024022543A1 (zh) 基于一致性算法的光伏电站分区分布式频率控制方法
CN107134795B (zh) 一种电热锅炉参与电网频率调节的控制系统及控制方法
CN114172202B (zh) 一种基于需求侧资源主动响应的含风电互联电力系统负荷频率控制方法
CN115733137A (zh) 通过压缩空气储能抑制电力系统低频振荡的系统及方法
CN113809760B (zh) 一种风电场参与电网二次调频的控制方法及装置
CN114123249B (zh) 一种基于电池储能主动响应的含风电互联电力系统负荷频率控制方法
CN109961185B (zh) 考虑电力系统灵活性综合协调配置的机组组合优化方法
CN108110803B (zh) 双馈风机辅助同步发电机参与电网二次调频协调控制方法
CN116470531B (zh) 一种应用于光伏场站的频率主动支撑控制方法及装置
CN115377998A (zh) 基于柔性负荷有功-电压耦合特性的电网频率控制方法
CN108539762A (zh) 一种基于自抗扰的风力发电机系统的频率控制系统及方法
CN110994669B (zh) 一种针对光伏电站集中式逆变器的控制方法及系统
CN112260306A (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