CN105591395A - 一种双馈风机虚拟惯量控制后的转速恢复方法 - Google Patents
一种双馈风机虚拟惯量控制后的转速恢复方法 Download PDFInfo
- Publication number
- CN105591395A CN105591395A CN201610006983.9A CN201610006983A CN105591395A CN 105591395 A CN105591395 A CN 105591395A CN 201610006983 A CN201610006983 A CN 201610006983A CN 105591395 A CN105591395 A CN 105591395A
- Authority
- CN
- China
- Prior art keywords
- double
- blower fan
- omega
- fed blower
- fed
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000011084 recovery Methods 0.000 title claims abstract description 20
- 230000006698 induction Effects 0.000 title abstract 3
- 230000000694 effects Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 5
- 238000013016 damping Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 abstract description 7
- 230000001133 acceleration Effects 0.000 abstract description 4
- 230000005611 electricity Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- HOWHQWFXSLOJEF-MGZLOUMQSA-N systemin Chemical compound NCCCC[C@H](N)C(=O)N[C@@H](CCSC)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CC(O)=O)C(=O)OC(=O)[C@@H]1CCCN1C(=O)[C@H]1N(C(=O)[C@H](CC(O)=O)NC(=O)[C@H](CCCN=C(N)N)NC(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H]2N(CCC2)C(=O)[C@H]2N(CCC2)C(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H](CCC(N)=O)NC(=O)[C@@H](NC(=O)[C@H](C)N)C(C)C)CCC1 HOWHQWFXSLOJEF-MGZLOUMQSA-N 0.000 description 1
- 108010050014 systemin Proteins 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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/24—Arrangements for preventing or reducing oscillations of power in networks
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P9/00—Arrangements for controlling electric generators for the purpose of obtaining a desired output
- H02P9/007—Control circuits for doubly fed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P9/00—Arrangements for controlling electric generators for the purpose of obtaining a desired output
- H02P9/04—Control effected upon non-electric prime mover and dependent upon electric output value of the generator
-
- 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]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P2101/00—Special adaptation of control arrangements for generators
- H02P2101/15—Special adaptation of control arrangements for generators for wind-driven turbines
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/30—Wind power
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Control Of Positive-Displacement Air Blowers (AREA)
Abstract
本发明涉及一种双馈风机虚拟惯量控制后的转速恢复方法,属于新能源发电中的控制技术领域。首先建立风力机传动链运动模型,并将运动模型转化为功率形式;将功率形式的运动模型转化为可设计扩张状态观测器的标准形式;设计扩张状态观测器,估计风力机输出的机械功率;整定扩张状态观测器参数;在估计得到的机械功率基础上设计风机加速功率并计算风机电磁功率参考值。相比于传统双馈风机虚拟惯量的比例/比例积分转速恢复策略,所提出的控制方法避免了比例/比例积分参数难以调节的问题,可以在快速恢复转速的同时避免***频率的二次跌落。
Description
技术领域
本发明涉及一种双馈风机虚拟惯量控制后的转速恢复方法,属于新能源发电中的控制技术领域。
背景技术
风力发电近年来在全球范围内得到了迅猛发展。随着风电渗透率的不断提高,风电的波动性、不确定性和低可控性对电力***的安全稳定运行有显著的影响。现在常用的双馈风机(Double-fedInductionGenerator,DFIG)通过电力电子变换器与电网相连,具有较高的风能捕获效率以及有功无功解耦控制能力,同时风机可以变转速运行,进行最大功率跟踪控制(Maximumpowerpointtracking,MPPT),提高了风能利用率。但在MPPT控制下,风机机械***与电磁***解耦,失去了传统同步机的惯性响应能力,这将导致***惯量大幅减小,发生功率扰动后***频率的变化率和偏移量增大,甚至超出限值,不利于***的稳定运行。
虽然双馈风机转子转速与***频率相互解耦,但是运行时其旋转器件(叶片,传动轴,发电机转子等)储存了可观的动能,这点与传统同步机组并无二致。若在风机控制***中加入附加功率控制环节,则可在***频率变化时释放一定的旋转动能给电网或者从电网吸收一定的电能转化为旋转动能,从而在短时间内改变风机有功出力,模拟传统同步发电机组的惯性响应以及一次调频响应,参与电力***频率控制。这种控制称作双馈风机虚拟惯量控制。
虚拟惯量控制下,若***发生频率跌落,风机为了增加出力会释放动能。虚拟惯量控制结束后,风机转速将低于MPPT控制所对应的最优转速,这将使得风机输出的机械功率减小。为了恢复到最优转速,风机还需进行转速恢复控制。通过比例(P)或者比例积分(PI)控制器进行转速恢复的方法简单易行,但该方法参数难以调节,P参数和I参数太小转速恢复太慢,P参数太大容易在开始转速恢复的时刻引起风机功率突然下降,I参数选择过大,虽然不会引起风机功率突然下降,但仍然会因为积分参数太大导致风机功率下降过快,这都会导致***频率的二次跌落,不利于电力***的频率稳定。由此可见,基于P或者PI控制器的转速恢复策略仍然存在缺陷。
发明内容
本发明的目的是提出一种双馈风机虚拟惯量控制后的转速恢复方法,采用扩张状态观测器估计双馈风机的机械功率,并在此基础上设计加速功率曲线进行转速恢复,可有效改善传统基于比例或者比例积分控制器转速恢复策略参数难以调节、容易引起***频率二次跌落的问题。
本发明提出的双馈风机虚拟惯量控制后的转速恢复方法,包括如下步骤:
(1)建立双馈风机的传动链运动模型如下:
其中,Hm为与双馈风机的叶片、轮毂、低速传动轴、齿轮箱、高速传动轴和发电机转子等效的质块的惯性时间常数,D为质块的自阻尼系数,Hm和D均为由双馈风机厂商提供的双馈风机固有参数,ωw为双馈风机的风轮转速(弧度/s),Te为双馈风机的异步发电机的电磁转矩,Tm为双馈风机的捕获的机械转矩;
在上式的两边同时乘以ωw,得到双馈风机的传动链运动模型的功率形式:
其中,Pm是双馈风机捕获的机械功率,Pe是双馈风机输出的电磁功率,Dωw 2是双馈风机传动过程中机械功率的损耗,双馈风机捕获的机械功率Pm为:
Pm=0.5ρCp(λ,β)πR2vw 3
上式中,ρ为空气密度(kg/m3);vw为风速(m/s);R为双馈风机的风轮半径(m),πR2为双馈风机的叶片扫过的面积(m2);Cp(λ,β)是双馈风机的风能利用系数,可近似表达为
其中λ为叶尖速比,是风轮叶尖的线速度与风速之比,ωw为双馈风机的风轮转速(弧度/s),a1~8是双馈风机风能利用系数表达式中的参数,根据实际风机风能利用特性拟合得到,λi是计算过程中的中间变量,β为双馈风机叶片的桨距角;
(2)将上述功率形式的双馈风机的传动链运动模型改写成一阶标准形式如下:
(3)上述步骤(2)的一阶标准形式的双馈风机的传动链运动模型 在扩张状态后改写为如下形式:
上式中,表示双馈风机的输出机械功率PM的微分,令状态扩张状态z2=PM/Hm,输入u=Pe,输入系数b=-1/Hm,构造如下的扩张状态观测器如下:
利用上述扩张状态观测器,计算得到双馈风机的输出机械功率PM的估计值和双馈风机转速平方的估计值
其中,是机械功率估计值的微分值,是转速平方估计值的微分值,α,δ,β01,β02分别为扩张状态观测器参数,其中,参数α取值为0.5,滤波因子取值范围为控制器采样时间的500-1500倍,数β01,β02的取值范围分别为50-500和500-2000;D和Hm为双馈风机自身参数,由双馈风机生产厂家测定并提供,转速ωw可由双馈风机角度编码器测得,电磁功率Pe由功率仪测量得到;
(4)根据上述双馈风机的输出机械功率PM,利用下式计算得到双馈风机的加速功率:
其中ωw为双馈风机转速,ωw0为双馈风机需要恢复到的最优转速,ωwrec为双馈风机虚拟惯量控制后开始恢复时的风机转速,ΔPC为设定的常数,取值范围为(0.05-0.5)×双馈风机额定功率,f(T)为一个随恢复时间T变化的常数;
(5)根据上述双馈风机的加速功率和估计得到的机械功率,计算得到转速恢复过程中双馈风机的电磁功率为实现虚拟惯量控制后转速恢复。
本发明提出的双馈风机虚拟惯量控制后的转速恢复方法,其优点是:
1、本发明方法能实时估计风机难以测量的捕获和输出的机械功率,并在此基础上合理安排风机加速功率,可以在双馈风机快速恢复转速的同时避免电网中二次频率跌落以及风机恢复过慢等问题。
2、本发明方法的控制过程简单明确,仅需修改风机参考功率指令即可实施。
3、本发明方法中的扩张状态观测器,采用非线性结构,抗干扰能力强,对状态的估计较为精确,对不同工况的适应性强,无需调整参数,即可在大多数工况下取得较好的控制效果。
附图说明
图1为本发明方法中的加速功率梯形递增函数f(T)。
图2为本发明方法的流程框图。
图3为本发明方法的实施例中,简单两区域测试***示意图。
图4为实施例中***频率变化曲线。
图5为实施例中风机转速变化曲线。
具体实施方式
本发明提出的双馈风机虚拟惯量控制后的转速恢复方法,其流程框图如图2所示,包括如下步骤:
(1)建立双馈风机的传动链运动模型如下:
其中,Hm为与双馈风机的叶片、轮毂、低速传动轴、齿轮箱、高速传动轴和发电机转子等效的质块的惯性时间常数,D为质块的自阻尼系数,Hm和D均为由双馈风机厂商提供的双馈风机固有参数,ωw为双馈风机的风轮转速(弧度/s),Te为双馈风机的异步发电机的电磁转矩,Tm为双馈风机的捕获的机械转矩;
在上式的两边同时乘以ωw,得到双馈风机的传动链运动模型的功率形式:
其中,Pm是双馈风机捕获的机械功率,Pe是双馈风机输出的电磁功率,是双馈风机传动过程中机械功率的损耗,双馈风机捕获的机械功率Pm为:
Pm=0.5ρCp(λ,β)πR2vw 3
上式中,ρ为空气密度(kg/m3);vw为风速(m/s);R为双馈风机的风轮半径(m),πR2为双馈风机的叶片扫过的面积(m2);Cp(λ,β)是双馈风机的风能利用系数,可近似表达为
其中λ为叶尖速比,是风轮叶尖的线速度与风速之比,ωw为双馈风机的风轮转速(弧度/s),a1~8是双馈风机风能利用系数表达式中的参数,具体数值可以根据实际风机风能利用特性拟合得到,一组典型取值为0.5176,116.0,0.4,5.0,-21.0,0.0068,0.08,0.035,λi是计算过程中的中间变量,β为双馈风机叶片的桨距角;
(2)将上述功率形式的双馈风机的传动链运动模型改写成一阶标准形式如下:
(3)扩张状态观测器是在一般观测器的基础上,将影响***被控输出的总扰动扩张成新的状态变量,然后对***状态变量和总扰动进行估计的一种非线性观测器,其具体原理如下。
对于如下非线性不确定一阶***:
其中x是状态变量,f(x,w(t))为未知函数,w(t)为未知外扰,u为控制输入,b为控制输入系数。
如果对b可以准确估计,则可令a(t)=f(x,w(t)),表示***的总扰动。
对于上述一阶非线性不确定***,取新的状态变量为:
则新的状态方程为:
对于上式所示的***,可以构造出如下的扩张状态观测器:
其中为状态观测量,β01,β02为扩张状态观测器的参数。函数fal(ε,α,δ)的表达式为:
由于双馈风机的输出机械功率PM等于双馈风机捕获的机械功率减去传动损耗的功率,即因此上述一阶标准形式的双馈风机的传动链运动模型 在扩张状态后改写为如下形式:
上式中,表示双馈风机的输出机械功率PM的微分,令状态扩张状态(干扰)z2=PM/Hm,输入u=Pe,输入系数b=-1/Hm,构造如下的扩张状态观测器如下:
利用上述扩张状态观测器,计算得到双馈风机的输出机械功率PM的估计值和双馈风机转速平方的估计值
其中,是机械功率估计值的微分值,是转速平方估计值的微分值,α,δ,β01,β02分别为扩张状态观测器参数,其中,参数α取值为0.5,滤波因子δ为比***采样时间稍大的数以加强滤波效果,具体数值可以根据实际实验效果适当调整,取值范围为控制器采样时间的500-1500倍,参数β01主要与状态z1的估计效果有关,参数β02与扩张状态z2的估计效果有关。参数β01,β02的取值越大,估计收敛的效果越好,但太大的值在实际***中容易引起扩张状态观测器输出的振荡,参数β01,β02的取值范围分别为50-500和500-2000,一种工程中实用的整定方法如下:先选择参数使得***稳定(实际上在很大范围内参数都可以使***稳定),然后再调整参数使得估计效果既快速又精确。
在双馈风机虚拟惯量转速恢复方法中,需要获得参数D,Hm,状态以及输入Pe。其中D和Hm是都是双馈风机自身参数,同一型号的双馈风机应当这些参数基本一致,可由生产厂家测定并提供。转速ωw可由角度编码器测得,电磁功率Pe可由功率仪测量得到。
(4)根据上述双馈风机的输出机械功率PM,利用下式计算得到双馈风机的加速功率:
其中ωw为双馈风机转速,ωw0为双馈风机需要恢复到的最优转速,ωwrec为双馈风机虚拟惯量控制后开始恢复时的风机转速,ΔPC为设定的常数,取值范围为(0.05-0.5)×双馈风机额定功率,f(T)为一个随恢复时间T变化的常数,(如图1所示,横坐标为时间,纵坐标为f(T)的取值,设定m和T2的取值,使得m的取值范围为-0.2到-0.05,T1为转速恢复开始时刻,T2-T1的值取值范围为15-45秒。此策略一方面通过梯形函数f(T),使得在转速恢复的最初阶段,风机有功出力不至于降得太低;另一方面,当风机转速ωw从ωwrec恢复到ωw0的过程中,系数(ωw0-ωw)/(ωw0-ωwrec)也逐渐减小到0。因此,相对采用比例或者比例积分(控制器进行转速恢复的方法,这种转速恢复策略可以更好地避免转速恢复过程对***频率造成的不利影响。)。从图1中可以看出,f(T)从小于0的数开始,这是为了进一步减小开始转速恢复时双馈风机有功出力的跌落。
(5)根据上述双馈风机的加速功率和估计得到的机械功率,计算得到转速恢复过程中双馈风机的电磁功率为实现虚拟惯量控制后转速恢复。
以下介绍本发明方法的一个实施例:
实验***中电力***模型采用简化后的电力***有功-频率响应模型,该模型中发电机方程采用二阶摇摆方程,忽略发电机的电压励磁调节,保留完整的调速器模型;***采用直流潮流模型。所有实验中的扰动均为在母线1增加0.1p.u.的负荷。电力***结构如图3所示、参数以及扰动参数如表1所示。具体实施本发明方法的步骤如下:
(1)建立双馈风机的传动链运动模型如下:
其中,Hm为与双馈风机的叶片、轮毂、低速传动轴、齿轮箱、高速传动轴和发电机转子等效的质块的惯性时间常数,本例中为5.4s,D为质块的自阻尼系数,由于这个参数通常较小,本例中忽略此参数,设其为0。ωw为双馈风机的风轮转速(弧度/s),Te为双馈风机的异步发电机的电磁转矩,Tm为双馈风机的捕获的机械转矩;
在上式的两边同时乘以ωw,得到双馈风机的传动链运动模型的功率形式:
其中,Pm是双馈风机捕获的机械功率,Pe是双馈风机输出的电磁功率,是双馈风机传动过程中机械功率的损耗,双馈风机捕获的机械功率Pm为:
Pm=0.5ρCp(λ,β)πR2vw 3
上式中,ρ为空气密度(kg/m3);vw为风速(m/s);R为双馈风机的风轮半径(m),πR2为双馈风机的叶片扫过的面积(m2);Cp(λ,β)是双馈风机的风能利用系数,可近似表达为
其中λ为叶尖速比,是风轮叶尖的线速度与风速之比,ωw为双馈风机的风轮转速(弧度/s),a1~8是双馈风机风能利用系数表达式中的参数,具体数值可以根据实际风机风能利用特性拟合得到,本例取值为0.5176,116.0,0.4,5.0,-21.0,0.0068,0.08,0.035,λi是计算过程中的中间变量,没有实际的物理意义,β为双馈风机叶片的桨距角;
本例试验***计算全部采用标幺值,对ρ,vw,R均不取具体值,只假设这些参数刚好满足在实验风速下MPPT控制时出力为0.2p.u.,风机最优转速为1p.u.。若风速发生改变,则按照 修改对应的最大功率,按照修改对应的最优转速。
(2)将上述功率形式的双馈风机的传动链运动模型改写成一阶标准形式如下:
(3)由于双馈风机的输出机械功率PM等于双馈风机捕获的机械功率减去传动损耗的功率,即因此上述一阶标准形式的双馈风机的传动链运动模型 在扩张状态后改写为如下形式:
上式中,表示双馈风机的输出机械功率PM的微分,令状态扩张状态(干扰)z2=PM/Hm,输入u=Pe,输入系数b=-1/Hm,构造如下的扩张状态观测器如下:
函数fal(ε,α,δ)的表达式为:
利用上述扩张状态观测器,计算得到双馈风机的输出机械功率PM的估计值和双馈风机转速平方的估计值
其中,是机械功率估计值的微分值,是转速平方估计值的微分值,α,δ,β01,β02分别为扩张状态观测器参数,经过调试,参数α取值为0.5,滤波因子δ取为采样时间的1000倍,参数β01,β02的取值分别取为280和1000。参数D和Hm都是已知参数,转速ωw可由角度编码器测得,电磁功率Pe可由功率仪测量到。
4)根据上述双馈风机的输出机械功率PM,利用下式计算得到双馈风机的加速功率:
其中ωw为双馈风机转速,ωw0为双馈风机需要恢复到的最优转速,ωwrec为双馈风机虚拟惯量控制后开始恢复时的风机转速,ΔPC为设定的常数,取值为0.2倍双馈风机额定功率,f(T)为一个随恢复时间T变化的常数,如附图1所示,横坐标为时间,纵坐标为f(T)的取值,设定m和T2的取值,本例中设m的取值-0.15,T1为转速恢复开始时刻,T2-T1的值取值为34.5秒。
根据上述双馈风机的加速功率和估计得到的机械功率,计算得到转速恢复过程中双馈风机的电磁功率为实现虚拟惯量控制后转速恢复。
此外,在虚拟惯量控制Kpf=0.5,Kdf=30(以***容量基值作为基值),持续时间为20s。作为对比目标采用比例控制器,参数设为:KP=0.8。
表1简单两母线电力***参数
由附图4和附图5的结果可知,与传统的基于比例控制器的转速恢复策略相比,该方法可以在保证转速恢复速度的同时避免频率的二次跌落。
Claims (1)
1.一种双馈风机虚拟惯量控制后的转速恢复方法,其特征在于该方法包括如下步骤:
(1)建立双馈风机的传动链运动模型如下:
其中,Hm为与双馈风机的叶片、轮毂、低速传动轴、齿轮箱、高速传动轴和发电机转子等效的质块的惯性时间常数,D为质块的自阻尼系数,Hm和D均为由双馈风机厂商提供的双馈风机固有参数,ωw为双馈风机的风轮转速(弧度/s),Te为双馈风机的异步发电机的电磁转矩,Tm为双馈风机的捕获的机械转矩;
在上式的两边同时乘以ωw,得到双馈风机的传动链运动模型的功率形式:
其中,Pm是双馈风机捕获的机械功率,Pe是双馈风机输出的电磁功率,是双馈风机传动过程中机械功率的损耗,双馈风机捕获的机械功率Pm为:
Pm=0.5ρCp(λ,β)πR2vw 3
上式中,ρ为空气密度(kg/m3);vw为风速(m/s);R为双馈风机的风轮半径(m),πR2为双馈风机的叶片扫过的面积(m2);Cp(λ,β)是双馈风机的风能利用系数,可近似表达为
其中λ为叶尖速比,是风轮叶尖的线速度与风速之比,ωw为双馈风机的风轮转速(弧度/s),a1~8是双馈风机风能利用系数表达式中的参数,根据实际风机风能利用特性拟合得到,λi是计算过程中的中间变量,β为双馈风机叶片的桨距角;
(2)将上述功率形式的双馈风机的传动链运动模型改写成一阶标准形式如下:
(3)上述步骤(2)的一阶标准形式的双馈风机的传动链运动模型 在扩张状态后改写为如下形式:
上式中,表示双馈风机的输出机械功率PM的微分,令状态扩张状态z2=PM/Hm,输入u=Pe,输入系数b=-1/Hm,构造如下的扩张状态观测器如下:
利用上述扩张状态观测器,计算得到双馈风机的输出机械功率PM的估计值和双馈风机转速平方的估计值
其中,是机械功率估计值的微分值,是转速平方估计值的微分值,α,δ,β01,β02分别为扩张状态观测器参数,其中,参数α取值为0.5,滤波因子取值范围为控制器采样时间的500-1500倍,数β01,β02的取值范围分别为50-500和500-2000;D和Hm为双馈风机自身参数,由双馈风机生产厂家测定并提供,转速ωw可由双馈风机角度编码器测得,电磁功率Pe由功率仪测量得到;
(4)根据上述双馈风机的输出机械功率PM,利用下式计算得到双馈风机的加速功率:
其中ωw为双馈风机转速,ωw0为双馈风机需要恢复到的最优转速,ωwrec为双馈风机虚拟惯量控制后开始恢复时的风机转速,ΔPC为设定的常数,取值范围为(0.05-0.5)×双馈风机额定功率,f(T)为一个随恢复时间T变化的常数;
(5)根据上述双馈风机的加速功率和估计得到的机械功率,计算得到转速恢复过程中双馈风机的电磁功率为实现虚拟惯量控制后转速恢复。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610006983.9A CN105591395B (zh) | 2016-01-06 | 2016-01-06 | 一种双馈风机虚拟惯量控制后的转速恢复方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610006983.9A CN105591395B (zh) | 2016-01-06 | 2016-01-06 | 一种双馈风机虚拟惯量控制后的转速恢复方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105591395A true CN105591395A (zh) | 2016-05-18 |
CN105591395B CN105591395B (zh) | 2017-10-31 |
Family
ID=55930734
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610006983.9A Active CN105591395B (zh) | 2016-01-06 | 2016-01-06 | 一种双馈风机虚拟惯量控制后的转速恢复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105591395B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107359832A (zh) * | 2017-06-28 | 2017-11-17 | 国家电网公司 | 一种非线性变结构励磁控制方法 |
CN107453410A (zh) * | 2017-09-04 | 2017-12-08 | 上海电力学院 | 负荷扰动的双馈风机参与风柴微网调频控制方法 |
CN108011381A (zh) * | 2017-11-30 | 2018-05-08 | 沈阳工业大学 | 一种风储一体化***调频控制方法 |
CN108131247A (zh) * | 2017-12-20 | 2018-06-08 | 北京金风科创风电设备有限公司 | 用于风力发电机组的数据处理方法和装置 |
CN108512251A (zh) * | 2018-04-25 | 2018-09-07 | 华北电力科学研究院有限责任公司 | 一种风机转速恢复方法及*** |
CN109193778A (zh) * | 2018-10-12 | 2019-01-11 | 浙江运达风电股份有限公司 | 一种基于风功率估计的惯量响应转速恢复控制方法 |
CN111859650A (zh) * | 2020-07-10 | 2020-10-30 | 中国科学院电工研究所 | 基于在线式联合仿真的风电机组传动链虚拟地面试验方法 |
CN112186814A (zh) * | 2020-09-10 | 2021-01-05 | 天津大学 | 一种双馈风机有功功率输出速降控制***及方法 |
CN112523945A (zh) * | 2020-12-29 | 2021-03-19 | 重庆邮电大学 | 一种双馈风力机最大风能捕获自抗扰非线性控制方法 |
CN113574270A (zh) * | 2019-03-22 | 2021-10-29 | 西门子歌美飒可再生能源公司 | 对风力涡轮机的惯性响应功率的估计 |
CN117108444A (zh) * | 2023-06-21 | 2023-11-24 | 国网甘肃省电力公司经济技术研究院 | 一种基于惯量需求的双馈风机附加惯量的控制方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103939290A (zh) * | 2014-04-17 | 2014-07-23 | 武汉新能源接入装备与技术研究院有限公司 | 一种适用于虚拟惯量优化控制研究的模拟风电机组 |
CN104836253A (zh) * | 2015-05-19 | 2015-08-12 | 清华大学 | 双馈风机虚拟惯量的控制方法及控制装置 |
-
2016
- 2016-01-06 CN CN201610006983.9A patent/CN105591395B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103939290A (zh) * | 2014-04-17 | 2014-07-23 | 武汉新能源接入装备与技术研究院有限公司 | 一种适用于虚拟惯量优化控制研究的模拟风电机组 |
CN104836253A (zh) * | 2015-05-19 | 2015-08-12 | 清华大学 | 双馈风机虚拟惯量的控制方法及控制装置 |
Non-Patent Citations (3)
Title |
---|
ARANTXA TAPIA, GERARDO TAPIA, J. XABIER OSTOLAZA, AND JOSÉ RAMÓN: "Modeling and Control of a Wind Turbine Driven Doubly Fed Induction Generator", 《IEEE TRANSACTIONS ON ENERGY CONVERSION 》 * |
J. G. SLOOTWEG,S. W. H. DE HAAN,H. POLINDER,W. L. KLING: "General Model for Representing Variable Speed Wind Turbines in Power System Dynamics Simulations", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 * |
李辉等: "风电机组等效模型对机组暂态稳定分析结果的影响", 《中国电机工程学报》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107359832A (zh) * | 2017-06-28 | 2017-11-17 | 国家电网公司 | 一种非线性变结构励磁控制方法 |
CN107453410A (zh) * | 2017-09-04 | 2017-12-08 | 上海电力学院 | 负荷扰动的双馈风机参与风柴微网调频控制方法 |
CN107453410B (zh) * | 2017-09-04 | 2020-03-20 | 上海电力学院 | 负荷扰动的双馈风机参与风柴微网调频控制方法 |
CN108011381A (zh) * | 2017-11-30 | 2018-05-08 | 沈阳工业大学 | 一种风储一体化***调频控制方法 |
CN108131247A (zh) * | 2017-12-20 | 2018-06-08 | 北京金风科创风电设备有限公司 | 用于风力发电机组的数据处理方法和装置 |
CN108131247B (zh) * | 2017-12-20 | 2020-09-29 | 北京金风科创风电设备有限公司 | 用于风力发电机组的数据处理方法和装置 |
CN108512251B (zh) * | 2018-04-25 | 2021-01-01 | 华北电力科学研究院有限责任公司 | 一种风机转速恢复方法及*** |
CN108512251A (zh) * | 2018-04-25 | 2018-09-07 | 华北电力科学研究院有限责任公司 | 一种风机转速恢复方法及*** |
CN109193778A (zh) * | 2018-10-12 | 2019-01-11 | 浙江运达风电股份有限公司 | 一种基于风功率估计的惯量响应转速恢复控制方法 |
CN113574270A (zh) * | 2019-03-22 | 2021-10-29 | 西门子歌美飒可再生能源公司 | 对风力涡轮机的惯性响应功率的估计 |
CN111859650A (zh) * | 2020-07-10 | 2020-10-30 | 中国科学院电工研究所 | 基于在线式联合仿真的风电机组传动链虚拟地面试验方法 |
CN111859650B (zh) * | 2020-07-10 | 2023-10-20 | 中国科学院电工研究所 | 基于在线式联合仿真的风电机组传动链虚拟地面试验方法 |
CN112186814A (zh) * | 2020-09-10 | 2021-01-05 | 天津大学 | 一种双馈风机有功功率输出速降控制***及方法 |
CN112186814B (zh) * | 2020-09-10 | 2022-04-19 | 天津大学 | 一种双馈风机有功功率输出速降控制***及方法 |
CN112523945A (zh) * | 2020-12-29 | 2021-03-19 | 重庆邮电大学 | 一种双馈风力机最大风能捕获自抗扰非线性控制方法 |
CN117108444A (zh) * | 2023-06-21 | 2023-11-24 | 国网甘肃省电力公司经济技术研究院 | 一种基于惯量需求的双馈风机附加惯量的控制方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN105591395B (zh) | 2017-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105591395A (zh) | 一种双馈风机虚拟惯量控制后的转速恢复方法 | |
Muljadi et al. | Understanding inertial and frequency response of wind power plants | |
Fingersh et al. | Baseline results and future plans for the NREL controls advanced research turbine | |
CN110454328A (zh) | 一种无风速传感器下的风力发电***功率追踪方法 | |
Muljadi et al. | Fixed-speed and variable-slip wind turbines providing spinning reserves to the grid | |
CN104779642A (zh) | 一种双馈风力发电机组的频率与阻尼综合控制方法 | |
CN104675629A (zh) | 一种变速风力发电机组的最大风能捕获方法 | |
Burnham et al. | Variable rotor-resistance control of wind turbine generators | |
Nam et al. | Feedforward pitch control using wind speed estimation | |
Cho et al. | Development and experimental verification of counter-rotating dual rotor/dual generator wind turbine: Generating, yawing and furling | |
CN106762405A (zh) | 一种能够抑制风力发电机组超速的控制方法及装置 | |
CN103746628A (zh) | 一种双馈感应风力发电机转子侧换流器的控制方法 | |
Krpan et al. | The mathematical model of a wind power plant and a gas power plant | |
Beltran et al. | A combined high gain observer and high-order sliding mode controller for a DFIG-based wind turbine | |
Sahin et al. | Performance comparison of two turbine blade pitch controller design methods based on equilibrium and frozen wake assumptions | |
CN104135207B (zh) | 一种大型双馈风力发电机并网联合优化控制方法 | |
CN114876732A (zh) | 一种风电机组变桨的控制方法及装置 | |
Le-peng et al. | Simulation for strategy of maximal wind energy capture of doubly fed induction generators | |
CN113852095A (zh) | 风电机组惯量响应控制方法及*** | |
Lak et al. | Speed control for direct drive permanent magnet wind turbine | |
Licari et al. | Coordinated inertia response from permanent magnet synchronous generator (PMSG) based wind farms | |
Guo et al. | Nonlinear control of wind power generation systems | |
CN105317632A (zh) | 一种风电机组转动惯量的测量方法 | |
Han et al. | Supplementary power control of PMSG-based wind farms for system dynamic stability | |
Van et al. | Improved pitch angle control for variable-speed wind turbine system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |