CN106194582B - 基于风速测量与估计的风电系统mppt控制装置及方法 - Google Patents
基于风速测量与估计的风电系统mppt控制装置及方法 Download PDFInfo
- Publication number
- CN106194582B CN106194582B CN201610828284.2A CN201610828284A CN106194582B CN 106194582 B CN106194582 B CN 106194582B CN 201610828284 A CN201610828284 A CN 201610828284A CN 106194582 B CN106194582 B CN 106194582B
- Authority
- CN
- China
- Prior art keywords
- wind
- speed
- output
- capacitance
- connect
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 88
- 239000000523 sample Substances 0.000 claims abstract description 43
- 230000008569 process Effects 0.000 claims abstract description 20
- 230000001537 neural effect Effects 0.000 claims description 46
- 210000002569 neuron Anatomy 0.000 claims description 45
- 210000000349 chromosome Anatomy 0.000 claims description 37
- 230000002068 genetic effect Effects 0.000 claims description 24
- 230000009194 climbing Effects 0.000 claims description 22
- 238000012549 training Methods 0.000 claims description 22
- 230000004927 fusion Effects 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 7
- 238000013528 artificial neural network Methods 0.000 claims description 5
- 230000035772 mutation Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 230000005611 electricity Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 210000005036 nerve Anatomy 0.000 description 2
- 241000288673 Chiroptera Species 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000007664 blowing Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- Wind Motors (AREA)
Abstract
本发明提供的基于风速测量与估计的风电系统MPPT控制装置及方法,风电系统包括永磁直驱风力发电机、整流器、第一电容、第二电容和负载,MPPT控制装置包括MPPT控制器、风速传感器、温度传感器、气压传感器、转速传感器、电压传感器、电流传感器、DC‑DC变换器和驱动模块。本发明对风速估计信号与地面风速测量信号进行融合,可有效提高寻找最大功率点的速度和准确度,降低扰动过程中的功率损失。
Description
技术领域
本发明涉及风力发电技术领域,更为具体地,涉及一种基于风速测量与估计的风电系统MPPT控制装置及方法。
背景技术
风力发电因其蕴藏量大、清洁可再生等特性,已成为当前新能源发电的主要方向之一。与双馈型风机相比,永磁直驱型风机具有无需变速箱、能耗少、维护成本低、更能适应低风速等优点。风力发电系统的输出功率与外界环境以及负载有关,在一定的外界环境下,存在唯一的最大功率点。为了提高风电系统的发电效率,需要为风力发电系统配备最大功率跟踪(Maximum Power Point Tracking,简称MPPT)控制装置,并实施有效的MPPT控制方法。风力发电系统的最大功率跟踪方法有很多,比如叶尖速比法,爬山搜索法,功率反馈法等。
叶尖速比法:在已知有效风速的条件下,通过调节风机转速,使叶尖速比保持最佳值,即可以实现最大功率跟踪。此方法的优点为思路简单,难点在于有效风速的获取。有效风速为整个风轮面积上的风速有效值,而风速传感器只能对某一点的风速进行测量,且装设较为困难,因此,有效风速的准确测量十分不易,从而导致该方法的实用性受限。
针对有效风速难以准确测量的问题,业界提出了基于风速估计的研究思路:首先测量风机转速与功率信息,再依据风力机的功率特性建立风速方程,进而求解出风速。该类方法的缺陷在于:一是风速方程为非线性方程,直接对其求解较为困难,不易于硬件实现;二是忽略了空气密度对风力机的功率特性影响,由于在不同的气温、气压条件下,空气密度不同,造成求解出的风速与真实有效风速的误差较大,从而影响了最大功率跟踪的准确性。
爬山搜索法:对风机转速或者占空比的控制指令值以一定的步长进行扰动,然后观察风机的功率变化。如果功率增加,那么风机转速或者占空比控制指令的扰动方向不变;如果功率减小,则将风机转速或者占空比控制指令的扰动反向。爬山搜索法的优点在于不需要进行风速测量,难点在于扰动周期以及扰动步长的选择。另外,当风机容量较大时,由于系统存在惯性,扰动并不能得到及时的功率响应,致使爬山搜索法的实际控制效果并不理想。
功率反馈法:此方法是通过转速传感器测量出风机的转速,利用发电机电磁转矩直接控制发电机的有功功率来跟踪已知的最佳转速-功率特性曲线,从而实现系统的最大功率跟踪控制。该方法的难度在于获取风机的最佳转速-功率曲线,需要大量的实验和分析,过程颇为繁杂。
综上所述,现有的风电系统的最大功率跟踪控制方法均存在一定的缺陷,限制了风电系统的最大功率捕获,因此有必要加以改进。
发明内容
鉴于上述背景技术所存在的问题,本发明提供一种基于风速测量与估计的风电系统MPPT控制装置及方法。
一种基于风速测量与估计的风电系统MPPT控制装置,风电系统包括永磁直驱风力发电机、整流器、第一电容、第二电容和负载,MPPT控制装置包括MPPT控制器、风速传感器、温度传感器、气压传感器、转速传感器、电压传感器、电流传感器、DC-DC变换器和驱动模块,其中,负载的一端与第二电容的正极连接,另一端接地,而第二电容的负极接地;整流器的三相输入端与永磁直驱风力发电机的三相输出端连接,整流器的单相输出正端与第一电容的正极连接,整流器的单相输出负端与第一电容的负极连接;电压传感器的待测电压输入正端与第一电容的正极连接,电压传感器的待测电压输入负端与第一电容的负极连接;电压传感器的测量信号输出端与MPPT控制器连接;电流传感器的待测电流输入端与第一电容的正极连接,电流传感器的待测电流输出端与DC-DC变换器的输入正端连接;电流传感器的测量信号输出端与MPPT控制器连接;DC-DC变换器的输入负端与第一电容的负极连接,DC-DC变换器的输出正端与第二电容的正极连接,DC-DC变换器的输出负端接地,DC-DC变换器的脉宽调制信号输入端与驱动模块的一端连接;驱动模块的另一端与MPPT控制器连接;转速传感器的两个输入端与永磁直驱风力发电机的三相输出端中的其中两端连接,转速传感器的测量信号输出端与MPPT控制器连接;风速传感器的测量信号输出端与MPPT控制器连接;温度传感器的测量信号输出端与MPPT控制器连接;气压传感器测量信号输出端与MPPT控制器连接。
本发明还提供上述的基于风速测量与估计的风电系统MPPT控制装置的控制方法,其工作步骤为:
步骤1:通过风速传感器采集地面某一高度Z0处的风速信号V0,并按照以下公式计算风机轮毂高度Zm处的风速信号Vm:
其中,α为风切变指数;
步骤2:分别利用电压传感器、电流传感器采集风电系统中整流器输出的直流电压Vz与直流电流Iz,利用转速传感器采集风机转速ωr,计算出此刻的永磁直驱风力发电机的输出机械功率Pm为:
其中,J为永磁直驱风力发电机(1)的转动惯量;
步骤3:通过温度传感器采集空气温度T,气压传感器采集大气压力p,并将风机转速ωr、机械功率Pm同空气温度T、大气压力p输入遗传-径向基神经网络风速估计模型,遗传-径向基神经网络风速估计模型输出的是风速的估计值Vm1;
步骤4:将风速的估计值Vm1与风速信号Vm按照以下方式进行融合,得出融合后的风速估计信号为:
Vmr=w*Vm1+(1-w)*Vm
其中,w为融合系数;
步骤5:利用最佳叶尖速比法,计算出有效风速的估计值Vmr所对应的最佳风机转速预测值ωref:
其中,λopt为最佳叶尖速比,R为叶轮半径;
步骤6:通过比例积分控制方法调节风机转速ωr,使之达到最佳风机转速预测值ωref;然后以ωref为初始值,采用占空比爬山搜索法以设定的扰动步长Δd跟踪风机的最大功率;
步骤7:若占空比爬山搜索法求得的扰动前后的功率差值大于等于设定阈值Tr,说明风速发生了突变,重复步骤1至步骤7;否则,继续采用占空比爬山搜索法以设定的扰动步长Δd跟踪风机的最大功率。
所述的基于风速测量与估计的风电系统MPPT控制方法,其中的遗传-径向基神经网络风速估计模型的建立过程如下:
步骤S1:采集风速-转速-气压-温度-功率训练样本,针对不同空气温度T与大气压力p,查阅空气在不同压力和温度下的密度表获取空气密度ρ,令风速信号Vw与转速信号ωr分别在一定范围内变化,并按照如下公式获取风机的机械功率输出信号Pm:
其中:R为风轮半径,Cp为风能利用系数,ρ为空气密度,λ为叶尖速比,通过构造大量风速-转速-气压-温度-功率样本,形成训练样本集{(Vw,ωr,T,p,Pm)};
步骤S2:建立遗传-径向基神经网络风速估计模型,输入层为4个神经元,隐含层为10个神经元,输出层为1个神经元,其中,输入层的第一个神经元对应输入为转速ωr,输入层的第二个神经元对应输入为温度T,输入层的第三个神经元对应输入为气压p,输入层的第四个神经元对应输入为风机机械功率Pm,输出层神经元对应输出为风速Vw,隐含层的基函数选择为高斯函数;
步骤S3:基于训练样本集{(Vw,ωr,T,p,Pm)},利用遗传算法对径向基神经网络风速估计模型的参数{(ci,δi,wi)}进行优化,其中,ci为径向基神经网络隐含层第i个神经元的基函数的数据中心,δi为径向基神经网络隐含层第i个神经元的基函数的扩展常数,wi为径向基神经网络隐含层第i个神经元到输出层的权值,具体步骤为:
步骤a):设置染色体的种群规模M,遗传终止进化代数T,染色体适应度函数,采用实数编码方式对染色体进行编码,并进行初始化;其中染色体的适应度函数为该染色体对应的径向基神经网络在训练样本集上的输出误差绝对值之和的倒数,即:
其中,表示该染色体对应的径向基神经网络在第k个训练样本上的输出,Yk表示目标输出值,m表示训练样本的总个数;
步骤b):计算种群中每个染色体的适应度函数值,并确定适应度函数值最大的为最优秀个体;
步骤c):判断本次遗传代数是否达到遗传终止进化代数T;如果达到,优化过程结束,最优秀个体对应的径向基神经网络模型即为最终模型;如果未达到,则进行选择、自适应交叉和变异运算,产生下代种群,将遗传代数增加1次,并转至步骤b)。
本发明能够取得以下技术效果:
1、对风速估计信号与风速测量信号进行融合作为最终的风速信号,提高其可靠性;
2、当外界环境发生变化时,借助于最佳叶尖速比法可以直接将工作转速调节至最佳风机转速的附近,避免了传统爬山搜索法逐步试探的过程,具有跟踪速度快的优点,从而提高了风机的发电效率。
附图说明
通过参考以下结合附图的说明及权利要求书的内容,并且随着对本发明的更全面理解,本发明的其它目的及结果将更加明白及易于理解。在附图中:
图1为根据本发明实施例的所采用的硬件结构图;
图2为根据本发明实施例的径向基神经网络风速估计模型的结构示意图;
图3为根据本发明实施例的遗传算法优化径向基神经网络的流程图;
图4为风速为7m/s时传统占空比爬山搜索法的跟踪过程示意图;
图5为风速为7m/s时本发明所提控制方法的跟踪过程示意图;
图6为变风速工况下本发明所提控制方法的跟踪过程示意图。
在所有附图中相同的标号指示相似或相应的特征或功能。
其中的附图标记包括:永磁直驱风力发电机1、整流器2、第一电容3、第二电容4、负载5、MPPT控制器6、风速传感器7、温度传感器8、气压传感器9、转速传感器10、电压传感器11、电流传感器12、DC-DC变换器13、驱动模块14。
附图中和文本中各符号的意义为:Vm为风速信号,ωr为永磁直驱风力发电机的转速(以下简称风机转速),Pm为永磁直驱风力发电机的机械功率(简称为风机机械功率),R为风轮半径,ρ为空气密度,Cp为风能利用系数,λ为叶尖速比,Vz为直流电压,Iz为直流电流,Vm1为径向基神经网络风速估计值,ωref为风机转速预测值;ci为径向基神经网络隐含层第i个神经元基函数的数据中心,δi为径向基神经网络隐含层第i个神经元基函数的扩展常数,wi为径向基神经网络隐含层第i个神经元到输出层的权值,φi(x)为隐含层基函数。
具体实施方式
在下面的描述中,出于说明的目的,为了提供对一个或多个实施例的全面理解,阐述了许多具体细节。然而,很明显,也可以在没有这些具体细节的情况下实现这些实施例。在其它例子中,为了便于描述一个或多个实施例,公知的结构和设备以方框图的形式示出。
以下将结合附图对本发明的具体实施例进行详细描述。
如图1所示,本发明实施例提供一种基于风速测量与估计的风电系统MPPT控制装置,其中,风电系统包括:永磁直驱风力发电机1、整流器2、第一电容3、第二电容4和负载5;MPPT控制装置包括:MPPT控制器6、风速传感器7、温度传感器8、气压传感器9、转速传感器10、电压传感器11、电流传感器12、DC-DC变换器13和驱动模块14。
负载5的一端与第二电容4的正极连接,负载5的另一端接地;整流器2的三相输入端与永磁直驱风力发电机1的三相输出端连接,整流器2的单相输出正端与第一电容3的正极连接,整流器2的单相输出负端与第一电容3的负极连接;电压传感器11的待测电压输入正端与第一电容3的正极连接,电压传感器11的待测电压输入负端与第一电容3的负极连接;电压传感器11的测量信号输出端与MPPT控制器6连接;电流传感器12的待测电流输入端与第一电容3的正极连接,电流传感器12的待测电流输出端与DC-DC变换器13的输入正端连接;电流传感器12的测量信号输出端与MPPT控制器6连接;DC-DC变换器13的输入负端与第一电容3的负极连接,DC-DC变换器13的输出正端与第二电容4的正极连接,第二电容4的负极接地,DC-DC变换器13的输出负端接地,DC-DC变换器13的脉宽调制信号输入端与驱动模块14的一端连接;驱动模块14的另一端与MPPT控制器6连接;转速传感器10的两个输入端与永磁直驱风力发电机1的三相输出端中的其中两端连接,转速传感器10的测量信号输出端与MPPT控制器6连接;风速传感器7的测量信号输出端与MPPT控制器6连接;温度传感器8的测量信号输出端与MPPT控制器6连接;气压传感器9的测量信号输出端与MPPT控制器6连接。
在本实施例中,永磁直驱风力发电机1的主要参数为:风轮直径为1.3m,额定电压为12V,额定转速为1400r/min,额定功率为200W,启动风速为1m/s,额定风速为10m/s,安全风速为25m/s;最大功率跟踪装置的主要参数为:MPPT控制器6采用dsPIC33FJ06GS101单片机,DC-DC变换器13采用Boost电路,驱动模块14为MCP14E3型号,电压传感器11为LV28-P型号,电流传感器12为LA25-NP型号,气压传感器9为MS5611-01BA03型号,温度传感器8为DS18b20型号,转速传感器10选用电压过零检测式频率计,第一电容3为C1=10uF,第二电容4为C2=100uF,负载5选择16欧电阻。
上述说明了MPPT控制装置的结构,以及与风电系统的连接关系。下面对MPPT控制装置的工作原理进行说明:
步骤1:通过风速传感器7采集地面某一高度Z0处的风速信号V0,并按照以下公式计算风机轮毂高度Zm处的风速信号Vm:
其中,α为风切变指数。
步骤2:分别利用电压传感器11、电流传感器12采集风电系统中整流器2输出的直流电压Vz与直流电流Iz,利用转速传感器10采集风机转速ωr,计算出此刻的永磁直驱风力发电机1的输出机械功率Pm为:
步骤3:通过温度传感器8采集空气温度T,气压传感器9采集大气压力p,并将风机转速ωr、机械功率Pm同空气温度T、大气压力p输入遗传-径向基神经网络风速估计模型,遗传-径向基神经网络风速估计模型输出的是风速的估计值Vm1。
步骤4:将风速的估计值Vm1与风速信号Vm按照以下方式进行融合,得出融合后的风速估计信号为:
Vmr=w*Vm1+(1-w)*Vm
其中,w为融合系数。
步骤5:利用最佳叶尖速比法,计算出有效风速的估计值Vmr所对应的最佳风机转速预测值ωref:
其中,λopt为最佳叶尖速比,R为叶轮半径。
步骤6:通过比例积分控制方法调节风机转速ωr,使之达到最佳风机转速预测值ωref;然后以ωref为初始值,采用占空比爬山搜索法以设定的扰动步长Δd跟踪风机的最大功率。
步骤7:若占空比爬山搜索法求得的扰动前后的功率差值大于等于设定阈值Tr,说明风速发生了突变,重复步骤1至步骤7;否则,继续采用占空比爬山搜索法以设定的扰动步长Δd跟踪风机的最大功率。
其中,步骤3中的遗传-径向基神经网络风速估计模型的建立过程如下:
步骤S1:采集风速-转速-气压-温度-功率训练样本,针对不同空气温度T与大气压力p,查阅空气在不同压力和温度下的密度表获取空气密度ρ,令风速信号Vw与转速信号ωr分别在一定范围内变化,并按照如下公式获取风机的机械功率输出信号Pm:
其中:R为风轮半径,Cp为风能利用系数,ρ为空气密度,λ为叶尖速比,通过构造大量风速-转速-气压-温度-功率样本,形成训练样本集{(Vw,ωr,T,p,Pm)}。
训练样本的采集过程为:针对不同空气温度T与大气压力p,令风速信号Vw与转速信号ωr分别在一定范围内变化,并按照公式(1)至公式(3)获取风机的机械功率输出信号Pm。在本实施例中,针对一台额定功率200W、风轮半径R=0.65m的永磁直驱风力发电机,令空气温度T以10℃为间隔在-30℃-+40℃内变化,令大气压力p以1kPa为间隔在97kPa-102kPa内变化,共计48种工况条件,通过查阅压力-温度-空气密度表计算相应的空气密度ρ,并令风速Vw信号在4m/s-13m/s内变化,风机转速ωr在10rad/s-140rad/s内变化,计算相应的风机机械功率Pm。考虑到实际应用时能量一般单向流动,故将输出功率为负值的样本进行剔除,共采集到132000个(T,p,Vw,ωr,Pm)样本。
步骤S2:建立遗传-径向基神经网络风速估计模型,输入层为4个神经元,隐含层为10个神经元,输出层为1个神经元。其中,输入层的第一个神经元对应输入为转速ωr,输入层的第二个神经元对应输入为温度T,输入层的第三个神经元对应输入为气压p,输入层的第四个神经元对应输入为风机机械功率Pm,输出层神经元对应输出为风速Vw。隐含层的基函数选择为高斯函数。
步骤S3:基于训练样本集{(Vw,ωr,T,p,Pm)},利用遗传算法对径向基神经网络风速估计模型的参数{(ci,δi,wi)}进行优化,其中,ci为径向基神经网络隐含层第i个神经元的基函数的数据中心,δi为径向基神经网络隐含层第i个神经元的基函数的扩展常数,wi为径向基神经网络隐含层第i个神经元到输出层的权值。具体步骤为:
步骤a):设置染色体的种群规模M,遗传终止进化代数T,染色体适应度函数,采用实数编码方式对染色体进行编码,并进行初始化;其中染色体的适应度函数为该染色体对应的径向基神经网络在训练样本集上的输出误差绝对值之和的倒数,即:
其中,表示该染色体对应的径向基神经网络在第k个训练样本上的输出,Yk表示目标输出值,m表示训练样本的总个数;
步骤b):计算种群中每个染色体的适应度函数值,并确定适应度函数值最大的个体为最优秀个体;
步骤c):判断本次遗传代数是否达到遗传终止进化代数T;如果达到,优化过程结束,最优秀个体对应的径向基神经网络模型即为最终模型;如果未达到,则进行选择、自适应交叉和变异运算,产生下代种群,将遗传代数增加1次,并转至步骤b);
图2示出了根据本发明实施例的径向基神经网络风速估计模型结构。径向基神经网络由输入层,隐含层,输出层组成,本发明中隐含层神经元的基函数选择为高斯函数:
其中,X表示输入向量,n为隐含层神经元的个数,ci表示隐含层第i个神经元的基函数的数据中心,δi为隐含层第i个神经元的基函数的扩展常数,范数||X-ci||表示输入向量X与数据中心δi的距离。对于输出层具有单个神经元的径向基神经网络,其输出为:
上述公式中,m代表输出隐含层神经元的个数,wi为径向基神经网络隐含层神经元到输出层的权值。
径向基神经网络的精度主要取决于三个主要参数,数据中心ci,扩展常数δi和权值wi,为了取得高精度的径向基神经网络预测模型,本发明采用遗传算法对其三个参数进行优化。
遗传算法由4个部分组成,分别为:
第一部分
染色体编码:编码方式一般有两类,一类为二进制法,一类为实属编码法,即将数据用一个常数表示,并赋予常数一定的取值范围进行编码;
第二部分
个体适应度评价:是通过构造适应度函数来实现,适应度函数为实际中的目标函数转换而来,表示解的优劣程度,也可以表示个体从父代中遗传到子代的遗传概率。适应度越高的个体遗传概率越高。
第三部分
遗传算子:遗传算法的遗传算子共有3个,分别为选择运算,交叉运算和变异运算。选择运算即计算适应度,然后进行选择,选择形式有局部选择和比例选择等;交叉运算即在两父代染色体中选择同一交叉点,并将该点前后两段互换并形成子代染色体;变异运算即对群体中的某些基因作更改,使算法保持多样性和局部搜索能力。
第四部分
运行参数:遗传算法的运行参数有群体大小M,遗传终止进化代数T,其中种群大小是指初始化种群时的群体规模,终止进化代数是指遗传的最末代数,这些参数对遗传算法的运行结果有着极大的影响。
如图3所示,本发明利用遗传算法对径向基神经网络风速预测模型进行优化的具体步骤如下:
步骤a):设置染色体的种群规模M,遗传终止进化代数T,染色体适应度函数,采用实数编码方式对染色体进行编码,并进行初始化;种群规模大可以增加遗传算法的选择空间,但会使计算速度下降;种群规模越小可以提高计算速度,但却容易使搜寻结果陷入局部极小值中。实施例中,设置种群规模M=30,遗传终止进化代数T=1000次,染色体的适应度函数为该染色体对应的径向基神经网络在训练样本集上的输出误差绝对值之和的倒数,即:
其中,表示该染色体对应的径向基神经网络在第k个训练样本上的输出,Yk表示目标输出值,m表示训练样本的总个数;
实施例中,对径向基神经网络隐含层神经元的数据中心Ci,扩展常数δi以及权值wi三个参数一起进行编码,组成染色体串。为了提高数据中心Ci和扩展常数δi在遗传过程中同时发生变化的概率,对Ci,δi两个参数进行交叉排列,即一个数据中心与一个扩展常数组成为一组,再将权重wi放在最后排列。
若输入层神经元个数为M,隐含层神经元个数为P,输出层神经元个数为N,则权值个数为MN,数据中心和宽度为PN+P个,染色体总长度为MN+PN+P。实施例中,M=4,P=10,N=1,则染色体总长度为24,编码方式采用实数编码法,采用随机数生成的方式对各染色体进行初始化;
步骤b):计算种群中每个染色体的适应度函数值,并确定适应度函数值最大的个体为最优秀个体;
步骤c):判断本次遗传代数是否达到遗传终止进化代数T;若达到,优化过程结束,最优秀个体对应的径向基神经网络模型即为最终模型;若未达到,则进行选择、交叉和变异运算,产生下代种群,将遗传代数增加1次,并转至步骤b);
在本实施方式中,选择运算采用比例选择法,交叉运算使用3点算数交叉法,变异运算采用均匀变异法。优化过程结束后,对种群中染色体的适应度函数值大小进行比较,最大值为0.0085,其所对应的染色体所代表的参数即为优化后的径向基神经网络模型参数。
结合具体实施方式对某一时刻的MPPT控制装置的控制方法进行说明:
1、通过风速传感器采集地面某一高度Z0处的风速信号V0,并按照以下公式计算风机轮毂高度Zm处的风速信号Vm:
其中,α为风切变指数;
在本实施方式中,取Z0=1m,Zm=8m,α=0.12;通过风速传感器采集地面风速信号V0=5.19m/s,计算得风机轮毂高度处的风速信号Vm=6.66m/s;
2、分别利用电压传感器、电流传感器采集风电系统中整流器输出的直流电压Vz与直流电流Iz,利用转速传感器采集风机转速ωr,计算出此刻的风机输出机械功率Pm为:
在本实施方式中,通过电压传感器采集的整流器输出的直流电压Vz=12.6V,直流电流Iz=8.2A,通过转速传感器采集到的风机转速ωr=93.8rad/s,得出此时的风机输出机械功率Pm=106.4W;
3、通过温度传感器采集空气温度T,气压传感器采集大气压力p,并将风机转速ωr、机械功率Pm同空气温度T、大气压力p输入遗传-径向基神经网络风速估计模型,风速估计模型输出为风速的估计值Vm1;
在本实施方式中,通过温度传感器采集空气温度T=20℃,通过气压传感器采集大气压力p=101kPa,将风机转速ωr=93.8rad/s、机械功率Pm=106.4W同空气温度T=20℃、大气压力p=101kPa输入遗传-径向基神经网络风速估计模型,风速估计模型输出为风速的估计值Vm1=6.84m/s;
4、将风速的估计值Vm1与风速信号Vm按照以下方式进行融合,得出融合后的风速估计信号为:
Vmr=w*Vm1+(1-w)*Vm;
在本实施方式中,取w=0.8,得出融合后的风速估计信号Vmr=6.81m/s;
5、利用最佳叶尖速比法,计算出有效风速的估计值Vmr所对应的最佳风机转速预测值ωref:
其中,λopt为最佳叶尖速比,R为叶轮半径;
在本实施方式中,取λopt=6.4,R=0.65m,计算出有效风速的估计值Vmr所对应的最佳风机转速预测值ωref=67.1rad/s;
6、通过比例积分控制方法调节风机的转速,使之达到最佳风机转速预测值ωref;然后以ωref为初始值,采用占空比爬山搜索法以设定的扰动步长Δd跟踪风机的最大功率;在实施例中,取扰动步长Δd=0.01;
7、若占空比爬山搜索法求得的扰动前后的功率差值大于等于设定阈值Tr,说明风速发生了突变,重复步骤1)至步骤7);否则,继续采用占空比爬山搜索法以设定的扰动步长Δd跟踪风机的最大功率。实例中,取阈值Tr=15W。
将上述的MPPT控制方法写入控制芯片中,输出PWM信号驱动DC-DC变换器,就可以实现最大功率跟踪功能。
为了说明本发明的功效,在不同风速工况下对本发明所提出的最大功率控制算法进行仿真验证。其中,占空比爬山搜索法的扰动周期为1s,比例积分控制器的控制周期为0.2s。
工况1:考虑风速恒定,设定风速Vw=7m/s。仿真时间t=2s时,系统进入稳态,此时风机转速ωr=93.8rad/s,风轮机械功率Pm=96.3W。仿真时间t=4s时,采用传统占空比爬山搜索法进行最大功率跟踪控制,扰动步长Δd=0.02,如图4所示。可知,需经过连续11拍的正向扰动后,风机转速ωr=67.9rad/s,风轮机械功率Pm=122.2W,开始进入稳态过程,即系统工作点开始在最大功率点附近波动。
在同样的工况下,仿真时间t=4s时,采用发明所提控制方法进行最大功率跟踪控制,其中占空比爬山搜索过程的扰动步长Δd=0.01,如图5所示。可知,首先本发明计算出融合后的风速估计信号Vmr=6.81m/s,由最佳叶尖速比法(λopt=6.4)计算出对应的最佳转速为67.1rad/s。利用PI控制器将风机调到此最佳转速,调整时间仅为1.4s,此时风机转速ωr=67.6rad/s,风轮机械功率Pm=121.5W,然后开始进行占空比爬山搜索过程,由于此时的工作点已经很接近最大功率点,仅仅经过2拍的正向扰动便进入了稳态过程,相对于占空比爬山搜索法而言大大减少了跟踪时间。
工况2:考虑变风速工况。设定风速按照正弦规律变化,即Vw=[7+sin(0.314t+1.57)]m/s。采用本发明所提控制方法的跟踪曲线如图6所示。可知,在变风速工况下,所提控制方法能够不断地对转速进行调节,从而使得系统工作点不断地向最大功率点靠近。在风速Vw=8m/s时,风轮机械功率Pm=182.4W;在风速为Vw=6m/s时,风轮机械功率Pm=73.1W,对比工况1,在风速Vw=7m/s条件下达到的最大功率点风轮机械功率Pm=122.4W,均符合功率与风速的三次方成正比这一关系,说明了本发明所提控制方法进行最大功率跟踪的有效性。
如上参照附图以示例的方式描述了根据本发明提出的基于风速测量与估计的风电系统MPPT控制装置及方法。但是,本领域技术人员应当理解,对于上述本发明所提出的基于风速测量与估计的风电系统MPPT控制装置,还可以在不脱离本发明内容的基础上对其中的实现细节做出各种改进。因此,本发明的保护范围应当由所附的权利要求书的内容确定。
Claims (2)
1.一种基于风速测量与估计的风电系统MPPT控制方法,风电系统包括永磁直驱风力发电机(1)、整流器(2)、第一电容(3)、第二电容(4)和负载(5),其MPPT控制装置包括MPPT控制器(6)、风速传感器(7)、温度传感器(8)、气压传感器(9)、转速传感器(10)、电压传感器(11)、电流传感器(12)、DC-DC变换器(13)和驱动模块(14);其中,所述负载(5)的一端与第二电容(4)的正极连接,另一端接地;所述整流器(2)的三相输入端与所述永磁直驱风力发电机(1)的三相输出端连接,所述整流器(2)的单相输出正端与所述第一电容(3)的正极连接,所述整流器(2)的单相输出负端与所述第一电容(3)的负极连接;所述电压传感器(11)的待测电压输入正端与所述第一电容(3)的正极连接,所述电压传感器(11)的待测电压输入负端与所述第一电容(3)的负极连接;所述电压传感器(11)的测量信号输出端与所述MPPT控制器(6)连接;所述电流传感器(12)的待测电流输入端与所述第一电容(3)的正极连接,所述电流传感器(12)的待测电流输出端与所述DC-DC变换器(13)的输入正端连接;所述电流传感器(12)的测量信号输出端与所述MPPT控制器(6)连接;所述DC-DC变换器(13)的输入负端与所述第一电容(3)的负极连接,所述DC-DC变换器(13)的输出正端与所述第二电容(4)的正极连接,所述第二电容(4)的负极接地,所述DC-DC变换器(13)的输出负端接地,所述DC-DC变换器(13)的脉宽调制信号输入端与所述驱动模块(14)的一端连接;所述驱动模块(14)的另一端与所述MPPT控制器(6)连接;所述转速传感器(10)的两个输入端与所述永磁直驱风力发电机(1)的三相输出端中的其中两端连接,所述转速传感器(10)的测量信号输出端与所述MPPT控制器(6)连接;所述风速传感器(7)的测量信号输出端与所述MPPT控制器(6)连接;所述温度传感器(8)的测量信号输出端与所述MPPT控制器(6)连接;所述气压传感器(9)的测量信号输出端与所述MPPT控制器(6)连接;其特征在于,所述MPPT控制方法的工作步骤为:
步骤1:通过风速传感器(7)采集地面某一高度Z0处的风速信号V0,并按照以下公式计算风机轮毂高度Zm处的风速信号Vm:
其中,α为风切变指数;
步骤2:分别利用电压传感器(11)、电流传感器(12)采集风电系统中整流器(2)输出的直流电压Vz与直流电流Iz,利用转速传感器(10)采集风机转速ωr,计算出此刻的永磁直驱风力发电机(1)的输出机械功率Pm为:
其中,J为永磁直驱风力发电机(1)的转动惯量;
步骤3:通过温度传感器(8)采集空气温度T,气压传感器(9)采集大气压力p,并将风机转速ωr、机械功率Pm同空气温度T、大气压力p输入遗传-径向基神经网络风速估计模型,遗传-径向基神经网络风速估计模型输出的是风速的估计值Vm1;
步骤4:将风速的估计值Vm1与风速信号Vm按照以下方式进行融合,得出融合后的风速估计信号为:
Vmr=w*Vm1+(1-w)*Vm
其中,w为融合系数;
步骤5:利用最佳叶尖速比法,计算出有效风速的估计值Vmr所对应的最佳风机转速预测值ωref:
其中,λopt为最佳叶尖速比,R为叶轮半径;
步骤6:通过比例积分控制方法调节风机转速ωr,使之达到最佳风机转速预测值ωref;然后以ωref为初始值,采用占空比爬山搜索法以设定的扰动步长Δd跟踪风机的最大功率;
步骤7:若占空比爬山搜索法求得的扰动前后的功率差值大于等于设定阈值Tr,说明风速发生了突变,重复步骤1至步骤7;否则,继续采用占空比爬山搜索法以设定的扰动步长Δd跟踪风机的最大功率。
2.根据权利要求1所述的一种基于风速测量与估计的风电系统MPPT控制方法,其特征在于所述遗传-径向基神经网络风速估计模型的建立过程如下:
步骤S1:采集风速-转速-气压-温度-功率训练样本,针对不同空气温度T与大气压力p,查阅空气在不同压力和温度下的密度表获取空气密度ρ,令风速信号Vw与转速信号ωr分别在一定范围内变化,并按照如下公式获取风机的机械功率输出信号Pm:
其中:R为风轮半径,Cp为风能利用系数,ρ为空气密度,λ为叶尖速比,通过构造大量风速-转速-气压-温度-功率样本,形成训练样本集{(Vw,ωr,T,p,Pm)};
步骤S2:建立遗传-径向基神经网络风速估计模型,输入层为4个神经元,隐含层为10个神经元,输出层为1个神经元,其中,输入层的第一个神经元对应输入为转速ωr,输入层的第二个神经元对应输入为温度T,输入层的第三个神经元对应输入为气压p,输入层的第四个神经元对应输入为风机机械功率Pm,输出层神经元对应输出为风速Vw,隐含层的基函数选择为高斯函数;
步骤S3:基于训练样本集{(Vw,ωr,T,p,Pm)},利用遗传算法对径向基神经网络风速估计模型的参数{(ci,δi,wi)}进行优化,其中,ci为径向基神经网络隐含层第i个神经元的基函数的数据中心,δi为径向基神经网络隐含层第i个神经元的基函数的扩展常数,wi为径向基神经网络隐含层第i个神经元到输出层的权值,具体步骤为:
步骤a):设置染色体的种群规模M,遗传终止进化代数T,染色体适应度函数,采用实数编码方式对染色体进行编码,并进行初始化;其中染色体的适应度函数为该染色体对应的径向基神经网络在训练样本集上的输出误差绝对值之和的倒数,即:
其中,表示该染色体对应的径向基神经网络在第k个训练样本上的输出,Yk表示目标输出值,m表示训练样本的总个数;
步骤b):计算种群中每个染色体的适应度函数值,并确定适应度函数值最大的为最优秀个体;
步骤c):判断本次遗传代数是否达到遗传终止进化代数T;如果达到,优化过程结束,最优秀个体对应的径向基神经网络模型即为最终模型;如果未达到,则进行选择、自适应交叉和变异运算,产生下代种群,将遗传代数增加1次,并转至步骤b)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610828284.2A CN106194582B (zh) | 2016-09-19 | 2016-09-19 | 基于风速测量与估计的风电系统mppt控制装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610828284.2A CN106194582B (zh) | 2016-09-19 | 2016-09-19 | 基于风速测量与估计的风电系统mppt控制装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106194582A CN106194582A (zh) | 2016-12-07 |
CN106194582B true CN106194582B (zh) | 2018-09-04 |
Family
ID=58067792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610828284.2A Active CN106194582B (zh) | 2016-09-19 | 2016-09-19 | 基于风速测量与估计的风电系统mppt控制装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106194582B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107122845A (zh) * | 2017-03-20 | 2017-09-01 | 华南理工大学 | 一种基于划分微风区的改进型神经网络架空线路风速预测方法 |
CN107229736B (zh) * | 2017-06-14 | 2021-10-22 | 北京唐浩电力工程技术研究有限公司 | 一种风电场风信息估测方法 |
CN107100794B (zh) * | 2017-06-14 | 2021-10-22 | 北京唐浩电力工程技术研究有限公司 | 一种风轮机能量感知及控制方法 |
CN109931230B (zh) * | 2017-12-19 | 2020-02-28 | 北京金风科创风电设备有限公司 | 检测风力发电机组的有功功率的方法和设备 |
CN108614612B (zh) * | 2018-04-24 | 2021-02-19 | 青岛高校信息产业股份有限公司 | 太阳能光伏电池最大功率追踪方法和系统 |
CN109185053B (zh) * | 2018-08-10 | 2020-01-17 | 无锡曼克斯电子科技有限公司 | 一种风力发电机低风速发电算法和实现方法、及其控制器 |
CN109779836B (zh) * | 2018-12-20 | 2020-02-21 | 明阳智慧能源集团股份公司 | 一种基于遗传算法寻优的风力发电机组发电量提升方法 |
CN112182974B (zh) * | 2020-10-09 | 2022-04-26 | 中国人民解放军陆军工程大学 | 一种神经元电路演化设计方法 |
CN112731805B (zh) * | 2020-12-01 | 2022-04-08 | 南京航空航天大学 | 一种基于风速估计的风力发电机最大功率跟踪无传感器鲁棒控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102678451A (zh) * | 2011-06-15 | 2012-09-19 | 华锐风电科技(集团)股份有限公司 | 一种风力发电机组变桨系统控制方法 |
CN103437955A (zh) * | 2013-08-13 | 2013-12-11 | 华北电力大学(保定) | 小型永磁直驱风力发电系统最大功率跟踪装置及控制方法 |
CN104454346A (zh) * | 2014-11-09 | 2015-03-25 | 华北电力大学(保定) | 一种小型永磁直驱风力发电系统最大功率跟踪控制方法 |
CN104806450A (zh) * | 2015-03-25 | 2015-07-29 | 华北电力大学(保定) | 一种基于万有引力神经网络的风电系统mppt控制方法 |
CN206386223U (zh) * | 2016-09-19 | 2017-08-08 | 华能新能源股份有限公司辽宁分公司 | 基于风速测量与估计的风电系统mppt控制装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008138636A (ja) * | 2006-12-05 | 2008-06-19 | Tokyo Metropolitan Univ | 風力発電システム |
-
2016
- 2016-09-19 CN CN201610828284.2A patent/CN106194582B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102678451A (zh) * | 2011-06-15 | 2012-09-19 | 华锐风电科技(集团)股份有限公司 | 一种风力发电机组变桨系统控制方法 |
CN103437955A (zh) * | 2013-08-13 | 2013-12-11 | 华北电力大学(保定) | 小型永磁直驱风力发电系统最大功率跟踪装置及控制方法 |
CN104454346A (zh) * | 2014-11-09 | 2015-03-25 | 华北电力大学(保定) | 一种小型永磁直驱风力发电系统最大功率跟踪控制方法 |
CN104806450A (zh) * | 2015-03-25 | 2015-07-29 | 华北电力大学(保定) | 一种基于万有引力神经网络的风电系统mppt控制方法 |
CN206386223U (zh) * | 2016-09-19 | 2017-08-08 | 华能新能源股份有限公司辽宁分公司 | 基于风速测量与估计的风电系统mppt控制装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106194582A (zh) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106194582B (zh) | 基于风速测量与估计的风电系统mppt控制装置及方法 | |
CN104806450B (zh) | 一种基于万有引力神经网络的风电系统mppt控制方法 | |
Mozafarpoor-Khoshrodi et al. | Improvement of perturb and observe method for maximum power point tracking in wind energy conversion system using fuzzy controller | |
CN107061158B (zh) | 一种低风速下风力发电机的预测与跟踪控制方法 | |
CN107045574B (zh) | 基于svr的风力发电机组低风速段有效风速估计方法 | |
CN104675629B (zh) | 一种变速风力发电机组的最大风能捕获方法 | |
CN111476293B (zh) | 一种计及风速风向的风电场分群方法 | |
CN106058937A (zh) | 一种含直驱永磁风电场的配电网广义负荷建模方法 | |
CN206386223U (zh) | 基于风速测量与估计的风电系统mppt控制装置 | |
Sarkar et al. | A study of MPPT schemes in PMSG based wind turbine system | |
Peng et al. | A hybrid DMST model for pitch optimization and performance assessment of high-solidity straight-bladed vertical axis wind turbines | |
Liew et al. | Wind characterization by three blade savonius wind turbine using iot | |
Hossain et al. | Energy efficient wind turbine system based on fuzzy control approach | |
Fang et al. | Design of Savonius model wind turbine for power catchment | |
Bertašienė et al. | Synergies of Wind Turbine control techniques | |
MURATOĞLU et al. | Numerical analyses of a straight bladed vertical axis Darrieus wind turbine: Verification of DMS algorithm and Qblade code | |
Haydlaar et al. | Optimization of power coefficient (Cp) in variable low rated speed wind turbine using increamental Particle Swarm Optimization (IPSO) | |
Ye et al. | Dynamic equivalent modeling approach of wind power plant with PMSG-WTGs | |
CN109611268A (zh) | 一种双叶轮水平轴风力机设计优化方法 | |
Wang et al. | ANN based pitch angle controller for variable speed variable pitch wind turbine generation system | |
Yusong et al. | The control strategy and simulation of the yaw system for MW rated wind turbine | |
Benaaouinate et al. | Emulation of Wind Turbine for Standalone Wind Energy Conversion Systems | |
Pehlivan | A novel fuzzy logic pitch angle controller with genetic algorithm optimization for wind turbines | |
Porselvan et al. | Design and Simulation of VAWT based on Power Supply for EV | |
Brahimi et al. | Aerodynamic analysis and performance prediction of vawt and hawt using cardaav and qblade computer codes |
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 |