CN109217374B - 一种风电电力***无功电压事前多时间尺度优化控制方法 - Google Patents

一种风电电力***无功电压事前多时间尺度优化控制方法 Download PDF

Info

Publication number
CN109217374B
CN109217374B CN201811346566.4A CN201811346566A CN109217374B CN 109217374 B CN109217374 B CN 109217374B CN 201811346566 A CN201811346566 A CN 201811346566A CN 109217374 B CN109217374 B CN 109217374B
Authority
CN
China
Prior art keywords
power
wind
reactive
output
voltage
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.)
Active
Application number
CN201811346566.4A
Other languages
English (en)
Other versions
CN109217374A (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN201811346566.4A priority Critical patent/CN109217374B/zh
Publication of CN109217374A publication Critical patent/CN109217374A/zh
Application granted granted Critical
Publication of CN109217374B publication Critical patent/CN109217374B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H02J3/386
    • 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
    • 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
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • 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
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/50Controlling the sharing of the out-of-phase component
    • 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
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003Load forecast, e.g. methods or systems for forecasting future load demand
    • 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
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A30/00Adapting or protecting infrastructure or their operation
    • 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

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明提供一种风电电力***无功电压事前多时间尺度优化控制方法,用于火电机组与风电场并网的电力***,解决现有技术中由于风速剧烈波动导致风电场无功调节能力减小而无法满足***电压稳定需求的技术问题,计及风电场有功出力对其无功可控容量的影响,兼顾风电场有功无功出力和***无功补偿装置对并网电电压的影响,通过建立事前多时间尺度优化模型协调风电场有功和无功出力、协调风电场和同步机的出力,在风速变化前主动调整风电场出力以避免风速波动对电网电压造成较大影响,在保证电网安全性的前提下,实现经济效益最大化。

Description

一种风电电力***无功电压事前多时间尺度优化控制方法
技术领域
本发明涉及风电控制领域,具体地说,涉及一种风电电力***无功电压事前多时间尺度优化控制方法。
背景技术
风力发电技术近年来在世界范围内得到了迅猛发展。风速波动具有随机性,风电机组的输出功率随风速波动变化,使电力***的有功和无功平衡面临挑战。随着可再生能源发电需求的增加,我国风电装机容量快速增长,风电正在向“大规模集中开发、远距离高压输送”方向快速发展。由于大型风电场大多位于电网末端,电网对风电并网点的电压支撑不足,大规模风电场的功率波动加剧使并网点电压问题更为突出,甚至造成***电压问题。
在风电场出口母线接入无功补偿装置是目前解决风电入网引起的无功电压问题的主要方法。但是,由于风速波动具有随机性和快速性的特点,通过投切电容器或电抗器平抑无功缺额难以满足补偿容量和响应速度的要求。静止无功补偿器可以对无功出力进行实时的平滑调节,是风电场无功补偿的常用手段。但是,风电场具有运行的多态性和输出的多样性特征,风电场中静止无功补偿器之间及其与风电机组之间的配合较为困难。特别是静止无功补偿器具有显著的滞后响应,无法快速响应风速波动,可能导致电网无功出力缺额,甚至出现电压过冲现象,造成风电机组连锁性脱网。
目前,以双馈风电机组和直驱风电机组为代表的变速风电机组已成为风力发电的主流机型。变速风电机组能够通过变流器控制***实现有功出力和无功出力的解耦控制,具有快速灵活的无功调控能力,适用于电压的动态调控。但是,目前风电场内的变速风电机组普遍采用单位功率因数运行,无功出力的控制优化值始终为零,远未发挥其可用于电压调控的潜力与优势。近些年,风电***的无功电压控制问题得到了大量的关注。但是,现有方法均是在风速波动后对***的无功出力进行调整,由于风速增大时风电场的无功极限随着有功出力的增大而减小,因此在风速剧烈波动下,风电机组的可控无功容量可能无法满足电网电压稳定的要求。
已有方法提出在风速波动前通过降低风电机组的有功出力来增大其可控无功容量,从而利用风电机组的无功出力控制平抑电压的变化。但是,风速波动下电网电压的变化源于风电机组有功出力的波动,即风电机组有功出力的调整不仅影响其可控无功容量,也会直接对电网电压产生直接的影响。忽略风电机组有功出力变化的影响,不仅难以充分挖掘风电机组的无功电压控制能力,无法最大限度地平抑电网电压的变化,还可能对电网安全性和经济性造成负面的效果。因此,亟待通过协调变速风电机组有功和无功出力实现***安全性和经济性的最大化。
发明内容
针对上述现有技术的不足,本发明提供一种风电电力***无功电压事前多时间尺度优化控制方法,用于火电机组与风电场并网的电力***,解决现有技术中由于风速剧烈波动导致风电场无功调节能力减小而无法满足***电压稳定需求的技术问题,计及风电场有功出力对其无功可控容量的影响,兼顾风电场有功无功出力和***无功补偿装置对并网电电压的影响,通过建立事前多时间尺度优化模型协调风电场有功和无功出力、协调风电场和同步机的出力,在风速变化前主动调整风电场出力以避免风速波动对电网电压造成较大影响,在保证电网安全性的前提下,实现经济效益最大化。
为了解决现有技术中的问题,本发明采用了如下的技术方案:
一种风电电力***无功电压事前多时间尺度优化控制方法,用于火电机组与风电场并网的电力***,其中,***中火电机组的数量为NG,风电场中双馈风电机组的数量为ND,直驱风电机组的数量为NP,包括如下步骤:
S101、采集风电并网点电压Uacc、风电场无功出力Qw、风电场有功出力Pw、电力***无功补偿量Qc,计算风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11,计算风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度a12
S102、获取下一时刻的风速预测信息,基于下一时刻的风速预测信息计算下一时刻风电场有功出力预测值
Figure GDA0002573357090000021
结合风电场有功出力Pw、风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11以及风电并网点电压Uacc,计算风电并网点电压预测值
Figure GDA0002573357090000022
S103、基于风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11及风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度a12,计算风速波动下平抑风电并网点电压变化的***无功需求变化量ΔQw,de;当ΔQw,de-ΔQc,max>NDQD.max+NPQP.max时,执行步骤S104,否则,认为电网无功容量足够,不采取控制措施,QD.max为双馈风电机组的无功可控范围,QP.max为直驱风电机组的无功可控范围,ΔQc,max为无功补偿装置的最大调节范围;
S104、利用包含电力***运行约束和弃风惩罚的事前长时间尺度出力优化模型,确定火电机组的开停机状态和下一时刻的电容器的无功补偿量;
S105、基于已确定的火电机组的开停机状态,利用包含电力***运行约束的事前短时间尺度出力优化模型,确定下一时刻的火电机组和风电场的有功出力、无功出力以及静止无功补偿装置的无功出力;
S106、基于下一时刻的火电机组和风电场的有功出力、无功出力,下一时刻的静止无功补偿装置的无功出力以及下一时刻的电容器的无功补偿量实现对风速波动下的电力***无功电压的事前电压控制。
优选地,步骤S101中,基于公式
Figure GDA0002573357090000031
计算风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11以及风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度a12,其中,ΔUacc、ΔUB1分别为风电并网点电压变化量和母线B1的电压变化量,ΔPw为流经升压站的有功功率变化量;ΔQw为流经升压站的无功功率变化量;a11为风电并网点电压Uacc关于风电场有功出力Pw的灵敏度,a21为母线B1的电压关于风电场有功出力Pw的灵敏度;a12为风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度,a22为母线B1的电压关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度。
优选地,步骤S102中,风电场的有功出力预测值
Figure GDA0002573357090000032
与下一时刻的风速预测信息之间的关系为:
Figure GDA0002573357090000033
其中,ρ为空气密度;R为风力机叶片半径;
Figure GDA0002573357090000034
为下一时刻的预测风速;Cpmax为最大风能利用系数;s为转差率;ND为风电场中双馈风电机组的数量;
风电并网点电压预测值
Figure GDA0002573357090000035
按照如下公式计算:
Figure GDA0002573357090000036
优选地,步骤S103中:
风速波动下平抑风电并网点电压变化的***无功需求变化量ΔQw,de按如下公式计算:
Figure GDA0002573357090000041
其中,ΔUpcc,al为风电并网点电压变化的最大允许范围,ΔPD为双馈风电机组输出有功功率的变化量,ΔPP为直驱风电机组输出有功功率的变化量;
双馈风电机组的无功可控范围QD.max按如下公式计算:
Figure GDA0002573357090000042
其中,PDs为双馈风电机组定子输出的有功功率;SDc为双馈风电机组网侧变流器容量;RSC、PSC及QSC均为中间变量,
Figure GDA0002573357090000043
Figure GDA0002573357090000044
Usm为双馈风电机组的机端电压幅值;Irmax为双馈风电机组转子允许流过的电流最大值;Rs、Xs分别为双馈风电机组的定子电阻、电抗;Xm为双馈风电机组的激磁电抗;s为转差率;
直驱风电机组的无功可控范围为QP.max按如下公式计算:
Figure GDA0002573357090000045
其中,SPc为直驱风电机组网侧变流器的容量;PP为直驱风电机组的有功出力。
优选地,设t-1时刻为当前时刻,t时刻为下一时刻,步骤S104中:
所述事前长时间尺度出力优化模型的目标函数为F1,其中,F1=min(f1+f2),f1表示事前长时间尺度出力优化模型的第一分函数,f2表示事前长时间尺度出力优化模型的第二分函数,
Figure GDA0002573357090000046
其中,t=1,2,...,Nt,Nt为长时间尺度的时段数;i=1,2,...,NG,NG为电力***中火电机组的总数目;k=1,2,...,Nw,Nw为电力***中风电机组的总数目;Bi.t表示第i个火电机组在t时刻开停机状态的0-1变量,Bi.t=max(Ui.t-Ui.t-1,0),Ui.t表示第i个火电机组在t时刻开机状态的0-1变量,Ui.t-1表示第i个火电机组在t-1时刻开机状态的0-1变量;Bi.t表示第i个火电机组在t时刻开停机状态的0-1变量;PG.i.t为长时间尺度下第i个火电机组在t时刻的有功出力;Si为第i个火电机组的启动成本;ai、bi、ci均为第i个火电机组的耗量特性参数;Cw为弃风惩罚的成本系数;
Figure GDA0002573357090000051
为长时间尺度下t时刻第k个风电机组的有功出力预测值;Pw.k.t为长时间尺度下t时刻第k个风电机组的有功出力;
事前长时间尺度出力优化模型的电力***运行约束包括:
风电并网点电压约束:
Uacc.min≤Uacc.t≤Uacc.max
其中,Uacc.t为长时间尺度下t时刻的风电并网点电压,Uacc,min为风电并网点允许的电压最小值,Uacc,max为风电并网点允许的电压最大值;
长时间尺度下t时刻的风电并网点电压Uacc.t与风电场有功出力、无功出力及无功补偿装置的无功出力的关系为:
Figure GDA0002573357090000052
其中,
Figure GDA0002573357090000053
为长时间尺度下t时刻风电并网点电压预测值,Qw.k.t为长时间尺度下t时刻第k个风电机组的无功出力,Pw.k.t-1、Qw.k.t-1分别为t-1时刻第k个风电机组的有功出力和无功出力,Qc.t、Qsvc.t分别为长时间尺度下t时刻电容器和静止无功补偿装置的无功出力;Qc.t-1、Qsvc.t-1分别为长时间尺度下t-1时刻电容器和静止无功补偿装置的无功出力;
静止无功补偿装置的无功约束:
Figure GDA0002573357090000054
其中,Qc.max是电容器的最大补偿容量;Qsvc.max是静止无功补偿装置的最大出力值;
电力***功率平衡约束:
Figure GDA0002573357090000055
其中,PL.t、QL.t分别为长时间尺度下t时刻电力***的总有功负荷和无功负荷;QG.i.t为长时间尺度下第i个火电机组在t时刻输出的无功出力;
风电场有功出力及无功出力约束:
Figure GDA0002573357090000061
其中,
Figure GDA0002573357090000062
为第k个风电机组在t时刻的有功出力预测值对应的无功输出极限值;
火电机组有功出力及无功出力约束:
Figure GDA0002573357090000063
其中,
Figure GDA0002573357090000064
分别为第i个火电机组的有功出力最大、最小值;
Figure GDA0002573357090000065
为第i个火电机组的额定视在功率;
火电机组爬坡能力约束:
Figure GDA0002573357090000066
其中,Ru.i、Rd.i分别为第i个火电机组的爬坡速率和滑坡速率;
火电机组最小开停机时间约束:
Figure GDA0002573357090000067
其中,
Figure GDA0002573357090000068
分别为第i个火电机组在t-1时刻已连续开机时间和已连续停机时间;
Figure GDA0002573357090000069
分别为第i个火电机组的最小连续开机时间和最小连续停机时间。
优选地,设t-1时刻为当前时刻,t时刻为下一时刻,步骤S105中:
事前短时间尺度出力优化模型的目标函数为F2,其中,
Figure GDA00025733570900000610
其中,i=1,2,...,NG,NG为电力***中火电机组的总数目;P′G.i.t为短时间尺度下第i个火电机组在t时刻的有功出力;ai、bi、ci均为第i个火电机组的耗量特性参数;Ui.t表示第i个火电机组在t时刻开机状态的0-1变量;
事前短时间尺度出力优化模型的电力***运行约束包括:
风电并网点电压约束:
Uacc.min≤U′acc.t≤Uacc.max
其中,U′acc.t为短时间尺度下t时刻的并网点电压,Uacc,min为风电并网点允许的电压最小值,Uacc,max为风电并网点允许的电压最大值;
静止无功补偿装置的无功约束:
0≤Q′svc.t≤Qsvc.max
其中,Q′svc.t为短时间尺度下静止无功补偿装置在t时刻的无功出力,Qsvc.max是静止无功补偿装置的最大出力值;
***功率平衡约束:
Figure GDA0002573357090000071
其中,P′L.t、Q′L.t分别为短时间尺度下t时刻电力***的总有功负荷和无功负荷;Q′G.i.t为短时间尺度下第i个火电机组在t时刻输出的无功出力;P′w.k.t、Q′w.k.t分别为短时间尺度下第k个风电机组在t时刻输出的有功出力和无功出力;
风电场有功出力和无功出力约束:
Figure GDA0002573357090000072
其中,
Figure GDA0002573357090000073
为短时间尺度下第k个风电机组在t时刻的有功出力预测值。
与现有技术相比,本发明具有以下有益效果:
1、与现有技术中风速波动后再进行电压控制不同,本发明是在风速波动前对即将发生变化的电网电压进行预测,在***的可控无功容量不能满足维持电网电压稳定所需的无功功率时,通过事前主动降低风电场的有功功率以提高风电场的无功调节能力,计及风电场有功出力对并网点电压的影响,在多时间尺度下协调***中各类装备的有功无功出力达到维持电网电压稳定的目的。
2、长时间尺度下,以火电机组发电成本最小为目标,引入风电弃风惩罚项,结合***运行约束优化计算出***中风电机组的有功、无功出力,火电机组的有功、无功出力以及电容器和静止无功补偿装置的无功功率。长时间尺度下确定的火电机组开停机状态和电容器的无功补偿量在短时间尺度下不再调整。
3、短时间尺度下,为了使电网经济效益更好,以火电机组运行成本最小为控制目标,在满足约束条件的情况下,对***中风电机组的有功、无功出力,火电机组的有功、无功出力以及静止无功补偿装置的无功功率进行最优分配,最终实现在保证电网安全性的前提下,实现经济效益最大化。
附图说明
为了使发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步的详细描述,其中:
图1为本发明的控制原理图;
图2为风电电力***无功电压事前多时间尺度优化控制方法的流程图;
图3为风电电力***的结构图实例。
具体实施方式
下面结合附图对本发明作进一步的详细说明。
如图3所示,为风电电力***的结构示意图,从左往右依次为变速风电场、高压交流输电线路及大电网,风电场经各自出口母线后,汇集到场群公共连接点B1,再经过升压站T1由高压交流输电线接入大电网,其中B2为风电场并网点。图中,B3表示连接大电网的母线。
如图1及图2所示,本发明公开了一种风电电力***无功电压事前多时间尺度优化控制方法,用于火电机组与风电场并网的电力***,其中,***中火电机组的数量为NG,风电场中双馈风电机组的数量为ND,直驱风电机组的数量为NP,包括如下步骤:
S101、采集风电并网点电压Uacc、风电场无功出力Qw、风电场有功出力Pw、电力***无功补偿量Qc,计算风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11,计算风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度a12
S102、获取下一时刻的风速预测信息,基于下一时刻的风速预测信息计算下一时刻风电场有功出力预测值
Figure GDA0002573357090000081
结合风电场有功出力Pw、风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11以及风电并网点电压Uacc,计算风电并网点电压预测值
Figure GDA0002573357090000082
S103、基于风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11及风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度a12,计算风速波动下平抑风电并网点电压变化的***无功需求变化量ΔQw,de;当ΔQw,de-ΔQc,max>NDQD.max+NPQP.max时,执行步骤S104,否则,认为电网无功容量足够,不采取控制措施,QD.max为双馈风电机组的无功可控范围,QP.max为直驱风电机组的无功可控范围,ΔQc,max为无功补偿装置的最大调节范围;
S104、利用包含电力***运行约束和弃风惩罚的事前长时间尺度出力优化模型,确定火电机组的开停机状态和下一时刻的电容器的无功补偿量;
S105、基于已确定的火电机组的开停机状态,利用包含电力***运行约束的事前短时间尺度出力优化模型,确定下一时刻的火电机组和风电场的有功出力、无功出力以及静止无功补偿装置的无功出力;
S106、基于下一时刻的火电机组和风电场的有功出力、无功出力,下一时刻的静止无功补偿装置的无功出力以及下一时刻的电容器的无功补偿量实现对风速波动下的电力***无功电压的事前电压控制。
具体实施时,步骤S101中,基于公式
Figure GDA0002573357090000091
计算风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11以及风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度a12,其中,ΔUacc、ΔUB1分别为风电并网点电压变化量和母线B1的电压变化量,ΔPw为流经升压站的有功功率变化量;ΔQw为流经升压站的无功功率变化量;a11为风电并网点电压Uacc关于风电场有功出力Pw的灵敏度,a21为母线B1的电压关于风电场有功出力Pw的灵敏度;a12为风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度,a22为母线B1的电压关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度。
具体实施时,步骤S102中,风电场的有功出力预测值
Figure GDA0002573357090000092
与下一时刻的风速预测信息之间的关系为:
Figure GDA0002573357090000093
其中,ρ为空气密度;R为风力机叶片半径;
Figure GDA0002573357090000094
为下一时刻的预测风速;Cpmax为最大风能利用系数;s为转差率;ND为风电场中双馈风电机组的数量;
风电并网点电压预测值
Figure GDA0002573357090000095
按照如下公式计算:
Figure GDA0002573357090000096
具体实施时,步骤S103中:
风速波动下平抑风电并网点电压变化的***无功需求变化量ΔQw,de按如下公式计算:
Figure GDA0002573357090000101
其中,ΔUpcc,al为风电并网点电压变化的最大允许范围,ΔPD为双馈风电机组输出有功功率的变化量,ΔPP为直驱风电机组输出有功功率的变化量;
双馈风电机组的无功可控范围QD.max按如下公式计算:
Figure GDA0002573357090000102
其中,PDs为双馈风电机组定子输出的有功功率;SDc为双馈风电机组网侧变流器容量;RSC、PSC及QSC均为中间变量,
Figure GDA0002573357090000103
Figure GDA0002573357090000104
Usm为双馈风电机组的机端电压幅值;Irmax为双馈风电机组转子允许流过的电流最大值;Rs、Xs分别为双馈风电机组的定子电阻、电抗;Xm为双馈风电机组的激磁电抗;s为转差率;
直驱风电机组的无功可控范围为QP.max按如下公式计算:
Figure GDA0002573357090000105
其中,SPc为直驱风电机组网侧变流器的容量;PP为直驱风电机组的有功出力。
具体实施时,设t-1时刻为当前时刻,t时刻为下一时刻,步骤S104中:
所述事前长时间尺度出力优化模型的目标函数为F1,其中,F1=min(f1+f2),f1表示事前长时间尺度出力优化模型的第一分函数,f2表示事前长时间尺度出力优化模型的第二分函数,
Figure GDA0002573357090000106
其中,t=1,2,...,Nt,Nt为长时间尺度的时段数;i=1,2,...,NG,NG为电力***中火电机组的总数目;k=1,2,...,Nw,Nw为电力***中风电机组的总数目;Bi.t表示第i个火电机组在t时刻开停机状态的0-1变量,Bi.t=max(Ui.t-Ui.t-1,0),Ui.t表示第i个火电机组在t时刻开机状态的0-1变量,Ui.t-1表示第i个火电机组在t-1时刻开机状态的0-1变量;Bi.t表示第i个火电机组在t时刻开停机状态的0-1变量;PG.i.t为长时间尺度下第i个火电机组在t时刻的有功出力;Si为第i个火电机组的启动成本;ai、bi、ci均为第i个火电机组的耗量特性参数;Cw为弃风惩罚的成本系数;
Figure GDA0002573357090000111
为长时间尺度下t时刻第k个风电机组的有功出力预测值;Pw.k.t为长时间尺度下t时刻第k个风电机组的有功出力;
事前长时间尺度出力优化模型的电力***运行约束包括:
风电并网点电压约束:
Uacc.min≤Uacc.t≤Uacc.max
其中,Uacc.t为长时间尺度下t时刻的风电并网点电压,Uacc,min为风电并网点允许的电压最小值,Uacc,max为风电并网点允许的电压最大值;
长时间尺度下t时刻的风电并网点电压Uacc.t与风电场有功出力、无功出力及无功补偿装置的无功出力的关系为:
Figure GDA0002573357090000112
其中,
Figure GDA0002573357090000113
为长时间尺度下t时刻风电并网点电压预测值,Qw.k.t为长时间尺度下t时刻第k个风电机组的无功出力,Pw.k.t-1、Qw.k.t-1分别为t-1时刻第k个风电机组的有功出力和无功出力,Qc.t、Qsvc.t分别为长时间尺度下t时刻电容器和静止无功补偿装置的无功出力;Qc.t-1、Qsvc.t-1分别为长时间尺度下t-1时刻电容器和静止无功补偿装置的无功出力;
静止无功补偿装置的无功约束:
Figure GDA0002573357090000114
其中,Qc.max是电容器的最大补偿容量;Qsvc.max是静止无功补偿装置的最大出力值;
电力***功率平衡约束:
Figure GDA0002573357090000115
其中,PL.t、QL.t分别为长时间尺度下t时刻电力***的总有功负荷和无功负荷;QG.i.t为长时间尺度下第i个火电机组在t时刻输出的无功出力;
风电场有功出力及无功出力约束:
Figure GDA0002573357090000121
其中,
Figure GDA0002573357090000122
为第k个风电机组在t时刻的有功出力预测值对应的无功输出极限值;
火电机组有功出力及无功出力约束:
Figure GDA0002573357090000123
其中,
Figure GDA0002573357090000124
分别为第i个火电机组的有功出力最大、最小值;
Figure GDA0002573357090000125
为第i个火电机组的额定视在功率;
火电机组爬坡能力约束:
Figure GDA0002573357090000126
其中,Ru.i、Rd.i分别为第i个火电机组的爬坡速率和滑坡速率;
火电机组最小开停机时间约束:
Figure GDA0002573357090000127
其中,
Figure GDA0002573357090000128
分别为第i个火电机组在t-1时刻已连续开机时间和已连续停机时间;
Figure GDA0002573357090000129
分别为第i个火电机组的最小连续开机时间和最小连续停机时间。
具体实施时,设t-1时刻为当前时刻,t时刻为下一时刻,步骤S105中:
事前短时间尺度出力优化模型的目标函数为F2,其中,
Figure GDA00025733570900001210
其中,i=1,2,...,NG,NG为电力***中火电机组的总数目;P′G.i.t为短时间尺度下第i个火电机组在t时刻的有功出力;ai、bi、ci均为第i个火电机组的耗量特性参数;Ui.t表示第i个火电机组在t时刻开机状态的0-1变量;
事前短时间尺度出力优化模型的电力***运行约束包括:
风电并网点电压约束:
Uacc.min≤U′acc.t≤Uacc.max
其中,U′acc.t为短时间尺度下t时刻的并网点电压,Uacc,min为风电并网点允许的电压最小值,Uacc,max为风电并网点允许的电压最大值;
静止无功补偿装置的无功约束:
0≤Q′svc.t≤Qsvc.max
其中,Q′svc.t为短时间尺度下静止无功补偿装置在t时刻的无功出力,Qsvc.max是静止无功补偿装置的最大出力值;
***功率平衡约束:
Figure GDA0002573357090000131
其中,P′L.t、Q′L.t分别为短时间尺度下t时刻电力***的总有功负荷和无功负荷;Q′G.i.t为短时间尺度下第i个火电机组在t时刻输出的无功出力;P′w.k.t、Q′w.k.t分别为短时间尺度下第k个风电机组在t时刻输出的有功出力和无功出力;
风电场有功出力和无功出力约束:
Figure GDA0002573357090000132
其中,
Figure GDA0002573357090000133
为短时间尺度下第k个风电机组在t时刻的有功出力预测值。
风电场输出的有功功率随着风速的波动而波动,可能对电压稳定带来不利影响,通过主动降低风电场的有功功率提高其无功可控范围,通过协调风电场有功无功出力和***无功补偿量能够使并网点电压维持在指定值附近,根据功率平衡,***中的火电机组出力能够平抑***有功缺额,根据协调***中各个装置的有功无功出力使***经济效益最大化。本发明中风电场以较小的功率损失,将并网点电压维持在了指定值附近,极大地提高了***的安全性。与只考虑风电场无功出力和***无功补量对并网点电压的影响进行对比,本发明所提方法更有利于***的经济性。
本发明考虑了风电场有功出力对无功可控容量的影响,兼顾了风电场有功无功出力和***无功补偿量对并网点电压的影响,以风电及负荷的滚动预测值为输入变量,通过建立事前多时间尺度优化模型协调风电场有功和无功出力、协调风电场出力和***无功补偿量,协调风电场和同步机的出力,在风速变化前通过主动调整各个装置的出力保证电压满足要求的情况下,最大限度地提高***经济性。
在发明中,事前短时间尺度出力优化模型以ΔT′(ΔT′<ΔT)为采样步长,对当前时刻的风电和负荷信息进行超短期预测,ΔT′即为短时间尺度下当前时刻和下一时刻的时间间隔;
事前长时间尺度出力优化模型每隔ΔT间隔启动,对未来MΔT时段内的风电和负荷数据进行预测”,ΔT即为长时间尺度下当前时刻和下一时刻的时间间隔,且MΔT=Nt
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管通过参照本发明的优选实施例已经对本发明进行了描述,但本领域的普通技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离所附权利要求书所限定的本发明的精神和范围。

Claims (6)

1.一种风电电力***无功电压事前多时间尺度优化控制方法,用于火电机组与风电场并网的电力***,其中,***中火电机组的数量为NG,风电场中双馈风电机组的数量为ND,直驱风电机组的数量为NP,其特征在于,包括如下步骤:
S101、采集风电并网点电压Uacc、风电场无功出力Qw、风电场有功出力Pw、电力***无功补偿量Qc,计算风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11,计算风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度a12
S102、获取下一时刻的风速预测信息,基于下一时刻的风速预测信息计算下一时刻风电场有功出力预测值
Figure FDA0002573357080000011
结合风电场有功出力Pw、风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11以及风电并网点电压Uacc,计算风电并网点电压预测值
Figure FDA0002573357080000012
S103、基于风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11及风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度a12,计算风速波动下平抑风电并网点电压变化的***无功需求变化量ΔQw,de;当ΔQw,de-ΔQc,max>NDQD.max+NPQP.max时,执行步骤S104,否则,认为电网无功容量足够,不采取控制措施,QD.max为双馈风电机组的无功可控范围,QP.max为直驱风电机组的无功可控范围,ΔQc,max为无功补偿装置的最大调节范围;
S104、利用包含电力***运行约束和弃风惩罚的事前长时间尺度出力优化模型,确定火电机组的开停机状态和下一时刻的电容器的无功补偿量;
S105、基于已确定的火电机组的开停机状态,利用包含电力***运行约束的事前短时间尺度出力优化模型,确定下一时刻的火电机组和风电场的有功出力、无功出力以及静止无功补偿装置的无功出力;
S106、基于下一时刻的火电机组和风电场的有功出力、无功出力,下一时刻的静止无功补偿装置的无功出力以及下一时刻的电容器的无功补偿量实现对风速波动下的电力***无功电压的事前电压控制。
2.如权利要求1所述的风电电力***无功电压事前多时间尺度优化控制方法,其特征在于,步骤S101中,基于公式
Figure FDA0002573357080000013
计算风电并网点电压Uacc关于风电场有功出力Pw的灵敏度a11以及风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度a12,其中,ΔUacc、ΔUB1分别为风电并网点电压变化量和母线B1的电压变化量,ΔPw为流经升压站的有功功率变化量;ΔQw为流经升压站的无功功率变化量;a11为风电并网点电压Uacc关于风电场有功出力Pw的灵敏度,a21为母线B1的电压关于风电场有功出力Pw的灵敏度;a12为风电并网点电压Uacc关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度,a22为母线B1的电压关于风电场无功出力Qw和电力***无功补偿量Qc的灵敏度。
3.如权利要求1所述的风电电力***无功电压事前多时间尺度优化控制方法,其特征在于,步骤S102中,风电场的有功出力预测值
Figure FDA0002573357080000021
与下一时刻的风速预测信息之间的关系为:
Figure FDA0002573357080000022
其中,ρ为空气密度;R为风力机叶片半径;
Figure FDA0002573357080000023
为下一时刻的预测风速;Cpmax为最大风能利用系数;s为转差率;
风电并网点电压预测值
Figure FDA0002573357080000024
按照如下公式计算:
Figure FDA0002573357080000025
4.如权利要求1所述的风电电力***无功电压事前多时间尺度优化控制方法,其特征在于,步骤S103中:
风速波动下平抑风电并网点电压变化的***无功需求变化量ΔQw,de按如下公式计算:
Figure FDA0002573357080000026
其中,ΔUpcc,al为风电并网点电压变化的最大允许范围,ΔPD为双馈风电机组输出有功功率的变化量,ΔPP为直驱风电机组输出有功功率的变化量;
双馈风电机组的无功可控范围QD.max按如下公式计算:
Figure FDA0002573357080000027
其中,PDs为双馈风电机组定子输出的有功功率;SDc为双馈风电机组网侧变流器容量;RSC、PSC及QSC均为中间变量,
Figure FDA0002573357080000028
Figure FDA0002573357080000031
Usm为双馈风电机组的机端电压幅值;Irmax为双馈风电机组转子允许流过的电流最大值;Rs、Xs分别为双馈风电机组的定子电阻、电抗;Xm为双馈风电机组的激磁电抗;s为转差率;
直驱风电机组的无功可控范围为QP.max按如下公式计算:
Figure FDA0002573357080000032
其中,SPc为直驱风电机组网侧变流器的容量;PP为直驱风电机组的有功出力。
5.如权利要求1所述的风电电力***无功电压事前多时间尺度优化控制方法,其特征在于,设t-1时刻为当前时刻,t时刻为下一时刻,步骤S104中:
所述事前长时间尺度出力优化模型的目标函数为F1,其中,F1=min(f1+f2),f1表示事前长时间尺度出力优化模型的第一分函数,f2表示事前长时间尺度出力优化模型的第二分函数,
Figure FDA0002573357080000033
其中,t=1,2,...,Nt,Nt为长时间尺度的时段数;i=1,2,...,NG,NG为电力***中火电机组的总数目;k=1,2,...,Nw,Nw为电力***中风电机组的总数目;Bi.t表示第i个火电机组在t时刻开停机状态的0-1变量,Bi.t=max(Ui.t-Ui.t-1,0),Ui.t表示第i个火电机组在t时刻开机状态的0-1变量,Ui.t-1表示第i个火电机组在t-1时刻开机状态的0-1变量;Bi.t表示第i个火电机组在t时刻开停机状态的0-1变量;PG.i.t为长时间尺度下第i个火电机组在t时刻的有功出力;Si为第i个火电机组的启动成本;ai、bi、ci均为第i个火电机组的耗量特性参数;Cw为弃风惩罚的成本系数;
Figure FDA0002573357080000034
为长时间尺度下t时刻第k个风电机组的有功出力预测值;Pw.k.t为长时间尺度下t时刻第k个风电机组的有功出力;
事前长时间尺度出力优化模型的电力***运行约束包括:
风电并网点电压约束:
Uacc.min≤Uacc.t≤Uacc.max
其中,Uacc.t为长时间尺度下t时刻的风电并网点电压,Uacc,min为风电并网点允许的电压最小值,Uacc,max为风电并网点允许的电压最大值;
长时间尺度下t时刻的风电并网点电压Uacc.t与风电场有功出力、无功出力及无功补偿装置的无功出力的关系为:
Figure FDA0002573357080000041
其中,
Figure FDA0002573357080000042
为长时间尺度下t时刻风电并网点电压预测值,Qw.k.t为长时间尺度下t时刻第k个风电机组的无功出力,Pw.k.t-1、Qw.k.t-1分别为t-1时刻第k个风电机组的有功出力和无功出力,Qc.t、Qsvc.t分别为长时间尺度下t时刻电容器和静止无功补偿装置的无功出力;Qc.t-1、Qsvc.t-1分别为长时间尺度下t-1时刻电容器和静止无功补偿装置的无功出力;
静止无功补偿装置的无功约束:
Figure FDA0002573357080000043
其中,Qc.max是电容器的最大补偿容量;Qsvc.max是静止无功补偿装置的最大出力值;
电力***功率平衡约束:
Figure FDA0002573357080000044
其中,PL.t、QL.t分别为长时间尺度下t时刻电力***的总有功负荷和无功负荷;QG.i.t为长时间尺度下第i个火电机组在t时刻输出的无功出力;
风电场有功出力及无功出力约束:
Figure FDA0002573357080000045
其中,
Figure FDA0002573357080000046
为第k个风电机组在t时刻的有功出力预测值对应的无功输出极限值;
火电机组有功出力及无功出力约束:
Figure FDA0002573357080000047
其中,
Figure FDA0002573357080000048
分别为第i个火电机组的有功出力最大、最小值;
Figure FDA0002573357080000049
为第i个火电机组的额定视在功率;
火电机组爬坡能力约束:
Figure FDA0002573357080000051
其中,Ru.i、Rd.i分别为第i个火电机组的爬坡速率和滑坡速率;
火电机组最小开停机时间约束:
Figure FDA0002573357080000052
其中,
Figure FDA0002573357080000053
分别为第i个火电机组在t-1时刻已连续开机时间和已连续停机时间;
Figure FDA0002573357080000054
分别为第i个火电机组的最小连续开机时间和最小连续停机时间。
6.如权利要求1所述的风电电力***无功电压事前多时间尺度优化控制方法,其特征在于,设t-1时刻为当前时刻,t时刻为下一时刻,步骤S105中:
事前短时间尺度出力优化模型的目标函数为F2,其中,
Figure FDA0002573357080000055
其中,i=1,2,...,NG,NG为电力***中火电机组的总数目;P′G.i.t为短时间尺度下第i个火电机组在t时刻的有功出力;ai、bi、ci均为第i个火电机组的耗量特性参数;Ui.t表示第i个火电机组在t时刻开机状态的0-1变量;
事前短时间尺度出力优化模型的电力***运行约束包括:
风电并网点电压约束:
Uacc.min≤U′acc.t≤Uacc.max
其中,U′acc.t为短时间尺度下t时刻的并网点电压,Uacc,min为风电并网点允许的电压最小值,Uacc,max为风电并网点允许的电压最大值;
静止无功补偿装置的无功约束:
0≤Q′svc.t≤Qsvc.max
其中,Q′svc.t为短时间尺度下静止无功补偿装置在t时刻的无功出力,Qsvc.max是静止无功补偿装置的最大出力值;
***功率平衡约束:
Figure FDA0002573357080000056
其中,P′L.t、Q′L.t分别为短时间尺度下t时刻电力***的总有功负荷和无功负荷;Q′G.i.t为短时间尺度下第i个火电机组在t时刻输出的无功出力;P′w.k.t、Q′w.k.t分别为短时间尺度下第k个风电机组在t时刻输出的有功出力和无功出力;
风电场有功出力和无功出力约束:
Figure FDA0002573357080000061
其中,
Figure FDA0002573357080000062
为短时间尺度下第k个风电机组在t时刻的有功出力预测值。
CN201811346566.4A 2018-11-13 2018-11-13 一种风电电力***无功电压事前多时间尺度优化控制方法 Active CN109217374B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811346566.4A CN109217374B (zh) 2018-11-13 2018-11-13 一种风电电力***无功电压事前多时间尺度优化控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811346566.4A CN109217374B (zh) 2018-11-13 2018-11-13 一种风电电力***无功电压事前多时间尺度优化控制方法

Publications (2)

Publication Number Publication Date
CN109217374A CN109217374A (zh) 2019-01-15
CN109217374B true CN109217374B (zh) 2020-08-25

Family

ID=64995510

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811346566.4A Active CN109217374B (zh) 2018-11-13 2018-11-13 一种风电电力***无功电压事前多时间尺度优化控制方法

Country Status (1)

Country Link
CN (1) CN109217374B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110635523B (zh) * 2019-08-26 2021-02-05 国电南瑞科技股份有限公司 考虑新能源有功计划影响的无功电压协调预控方法和装置
CN110673569B (zh) * 2019-10-25 2022-09-23 国网安徽省电力有限公司电力科学研究院 一种绿色能源与火力发电机组协同控制方法及装置
CN111030142B (zh) * 2020-01-02 2021-03-23 重庆大学 一种高比例风电电力***多时间尺度频率优化控制方法
CN111342446B (zh) * 2020-03-02 2023-10-27 国网天津市电力公司电力科学研究院 三维家庭智慧用能泛在调度控制方法及***
CN112039077B (zh) * 2020-08-26 2022-07-05 武汉大学 一种考虑电压优化控制约束的有功频率调整方法
CN112260327B (zh) * 2020-10-22 2022-09-30 国网河南省电力公司电力科学研究院 基于特高压交直流电网的新能源无功支撑分析方法
CN117833378A (zh) * 2024-01-09 2024-04-05 国网江苏省电力有限公司淮安供电分公司 一种无avc主站的风电场电压控制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102611118A (zh) * 2012-03-14 2012-07-25 清华大学 一种引入预测信息的风电场综合无功电压控制方法
CN108695907A (zh) * 2018-06-26 2018-10-23 四川大学 一种微电网多时间尺度优化调度方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102738833B (zh) * 2012-06-20 2014-07-09 湖北省电力公司 一种含风电的电力***多时间尺度滚动协调调度方法
JP2014176257A (ja) * 2013-03-12 2014-09-22 Tokyo Denki Univ 風力発電システム
CN103441506B (zh) * 2013-06-18 2017-05-10 国家电网公司 不同时间尺度下多目标协调分散风电场无功优化控制方法
CN108090632B (zh) * 2018-01-23 2020-06-02 南方电网科学研究院有限责任公司 基于鲁棒优化的新能源并网电力***多时间尺度调度方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102611118A (zh) * 2012-03-14 2012-07-25 清华大学 一种引入预测信息的风电场综合无功电压控制方法
CN108695907A (zh) * 2018-06-26 2018-10-23 四川大学 一种微电网多时间尺度优化调度方法

Also Published As

Publication number Publication date
CN109217374A (zh) 2019-01-15

Similar Documents

Publication Publication Date Title
CN109217374B (zh) 一种风电电力***无功电压事前多时间尺度优化控制方法
CN106786806B (zh) 一种基于模型预测控制的配电网有功无功协调调控方法
Lamsal et al. Output power smoothing control approaches for wind and photovoltaic generation systems: A review
CN102299527B (zh) 一种风电场无功功率控制方法和***
EP2232667B1 (en) Method and system for controlling wind power plants
CN109066770B (zh) 一种风电接入柔性直流输电***的控制方法和装置
CN104467030B (zh) 一种基于风电与火电联合发电的功率分配方法
CN103219732A (zh) 一种含变速恒频风电场的配电网无功电压控制方法
CN102354992A (zh) 风电场无功功率控制方法
CN103346577A (zh) 降低风电场功率损耗的风电场avc无功控制***及方法
CN103490426A (zh) 大规模风电接入的弱送端电力***电压稳定综合控制方法
CN108494020B (zh) 风速波动下风电送出***无功电压主动控制方法
CN104578086A (zh) 一种风电和光伏发电接入电网的无功电压控制方法
CN105262098A (zh) 基于风电场发电功率波动评估的敏捷自动电压控制方法
CN108711868A (zh) 一种计及孤岛运行电压安全的配电网无功优化规划方法
CN105470978B (zh) 一种静止无功补偿装置成组协调控制方法
CN109936146A (zh) 一种基于改进灵敏度算法的风电场协调优化控制方法
CN114759620A (zh) 一种风光储场站群无功协同优化调控方法、装置及***
CN108183510A (zh) 面向风速波动的双馈风电***无功功率主动控制方法
CN112332462A (zh) 考虑源-荷功率随机波动特性的双馈风力发电机组一次频率平滑调节方法
CN111027179B (zh) 一种计及辅助调频服务的双馈风电场的等值建模方法
Fang et al. An optimal control strategy for reactive power in wind farms consisting of VSCF DFIG wind turbine generator systems
CN110611332A (zh) 一种海上风电***储能装置及其控制方法
CN202178583U (zh) 一种风电场无功功率控制***
CN106972553A (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