CN116378898B - 一种风电系统桨距角控制方法、装置、系统及存储介质 - Google Patents
一种风电系统桨距角控制方法、装置、系统及存储介质 Download PDFInfo
- Publication number
- CN116378898B CN116378898B CN202310654046.4A CN202310654046A CN116378898B CN 116378898 B CN116378898 B CN 116378898B CN 202310654046 A CN202310654046 A CN 202310654046A CN 116378898 B CN116378898 B CN 116378898B
- Authority
- CN
- China
- Prior art keywords
- wind
- speed
- pitch angle
- wind power
- fan
- 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 48
- 238000010248 power generation Methods 0.000 claims abstract description 70
- 238000013178 mathematical model Methods 0.000 claims abstract description 13
- 238000004590 computer program Methods 0.000 claims description 10
- 238000010276 construction Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 20
- 238000010586 diagram Methods 0.000 description 14
- 230000008859 change Effects 0.000 description 13
- 230000003044 adaptive effect Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 230000001276 controlling effect Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000009466 transformation Effects 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
-
- 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
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/022—Adjusting aerodynamic properties of the blades
- F03D7/0224—Adjusting blade pitch
-
- 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
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/04—Automatic control; Regulation
- F03D7/042—Automatic control; Regulation by means of an electrical or electronic controller
- F03D7/043—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic
- F03D7/045—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic with model-based controls
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/327—Rotor or generator speeds
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/328—Blade pitch angle
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/40—Type of control system
- F05B2270/402—Type of control system passive or reactive, e.g. using large wind vanes
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/70—Type of control algorithm
- F05B2270/702—Type of control algorithm differential
-
- 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)
- 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)
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种风电系统桨距角控制方法、装置、系统及存储介质,属于风力发电技术领域,方法包括:根据预获取的风力发电系统的数学模型,构建获取风机转速非线性微分方程;根据所述风机转速非线性微分方程,对风力发电系统中未知量风能利用系数进行估计,获取估计参数;根据预构建的风能利用系数的非线性模型,利用所述估计参数,求解获取桨距角,并将所述桨距角反馈给风力发电系统,实现风力发电系统桨距角控制。该方法能够通过对风能利用系数进行估计,实现风力发电系统桨距角控制。
Description
技术领域
本发明涉及一种风电系统桨距角控制方法、装置、系统及存储介质,属于风力发电技术领域。
背景技术
为推动社会生产向资源节约型转型,发展清洁能源,提高资源利用率成为了当务之急。风能作为绿色清洁能源,在我国分布广泛,因此提高风能利用率,加强风能技术开发成为了当前新能源开发的重中之重。
变桨距技术是相较于定桨距技术的改进。定桨距即叶片的桨距角固定,不会随风速的变化而变化,风轮的风能利用率完全由气动性来调节。当风轮转速在额定转速以下时,定桨距技术拥有较好的控制效果,当风速过大,风轮转速超过额定值时,由于没有功率反馈系统,定桨距控制无法维持功率的恒定,系统的功率会持续下降,风能利用系数降低,系统的载荷也会大大增加,降低了使用寿命。变桨距控制通过增加变桨距角执行机构,当叶片转速大于额定值时,桨距角执行机构启动,控制桨距角缓慢变化,使得系统的功率能维持在额定功率附近,并且拥有较高的风能利用系数。独立变桨距就是在统一变桨控制的基础上,给每个桨距角额外附加一个不同的桨距角信号,以优化风轮的空气动力学特性,达到降低不平衡载荷的目的。
风电机是高度非线性化的模型,传统的PID控制对其控制效果有限,无法实现桨距角的连续控制。自适应控制是针对非线性模型的一种有效控制方法,通过设计合适的控制律,在保证系统稳定的同时,提高系统的精度和动态特性。文献(韩兵,马杰.基于模型参考自适应的大型风电机组独立变桨控制方法.分布式能源,2021,6(05):26-32.)提出了基于模型参考的自适应方法,仿真结果表明风机的输出功率能保持稳定,现有的大多文献中的控制方案大多是对桨距角进行参数估计,通过调节控制器参数进而对桨距角直接进行控制,这类方案在一定程度上增加了观测器数目,提高了系统的复杂度和设计成本。
发明内容
本发明的目的在于提供一种风电系统桨距角控制方法、装置、系统及存储介质,通过对风能利用系数进行估计,实现风力发电系统桨距角控制。
为达到上述目的,本发明提供如下技术方案:
第一方面,本发明提供一种风电系统桨距角控制方法,包括:
根据预获取的风力发电系统的数学模型,构建获取风机转速非线性微分方程;
根据所述风机转速非线性微分方程,对风力发电系统中未知量风能利用系数进行估计,获取估计参数;
根据预构建的风能利用系数的非线性模型,利用所述估计参数,求解获取桨距角,并将所述桨距角反馈给风力发电系统,实现风力发电系统桨距角控制。
结合第一方面,进一步的,所述风力发电系统的数学模型的表达式如公式(1)所示:
(1)
公式(1)中,为风轮转动惯量,/>为风机转速,/>为风机转速的导数,/>为风机机械转矩,/>为齿轮比,/>为电机转矩,/>为空气密度,/>为叶轮长度,/>为风速,/>为叶尖速比,为风能利用系数。
结合第一方面,进一步的,所述风机转速非线性微分方程的表达式如公式(2)所示:
(2)
公式(2)中,为风机转速的导数,/>为常数项,/>,/>为空气密度,/>为叶轮长度,/>为风轮转动惯量,/>为风速,/>为叶尖速比,/>为风能利用系数,/>为齿轮比,为电机转矩。
结合第一方面,进一步的,根据所述风机转速非线性微分方程,对风力发电系统中未知量风能利用系数进行估计,获取估计参数包括:
引入估计参数和风力发电系统的控制律,对所述风机转速非线性微分方程进行改写,获取风机转速非线性微分方程改写模型;
定义转速误差和估计误差,并构建获取关于转速误差和估计误差的李雅普诺夫函数;
根据所述风机转速非线性微分方程改写模型和李雅普诺夫函数,基于李雅普诺夫第二定理,结合等价确定性原则,求解获取风力发电系统的控制律;
根据所述风力发电系统的控制律,求解获取估计参数。
结合第一方面,进一步的,所述风机转速非线性微分方程改写模型的表达式如公式(3)所示:
(3)
公式(3)中,为风机转速的导数,/>为常数项,/>,/>为空气密度,/>为叶轮长度,/>为风轮转动惯量,/>为风速,/>为叶尖速比,/>为估计参数,/>为齿轮比,/>为电机转矩,/>为风力发电系统的控制律;
所述关于转速误差和估计误差的李雅普诺夫函数的表达式如公式(4)所示:
(4)
公式(4)中,为关于转速误差和估计误差的李雅普诺夫函数,/>为转速误差,/>,/>为额定转速,/>为风机转速,/>为估计误差,/>,/>为风能利用系数,/>为估计参数;
所述风力发电系统的控制律的表达式如公式(5)所示:
(5)
公式(5)中,为额定转速的导数,/>为任意大于零的常数;
所述估计参数的表达式如公式(6)所示:
(6)
公式(6)中,为风机工作时间。
结合第一方面,进一步的,所述风能利用系数的非线性模型的表达式如公式(7)所示:
(7)
公式(7)中,为风能利用系数,/>为叶尖速比,/>为桨距角,/>为关于桨距角/>和叶尖速比/>的函数,/>为转速误差,/>为风机转速,/>为叶轮长度,/>为风速。
结合第一方面,进一步的,根据预构建的风能利用系数的非线性模型,利用所述估计参数,求解获取桨距角包括:
将所述风能利用系数的非线性模型的表达式中的风能利用系数用求解获取的估计参数/>代替,获取仅包含桨距角/>一个未知数的非线性方程组;
对所述仅包含桨距角一个未知数的非线性方程组进行求解,获取桨距角/>。
第二方面,本发明提供一种风电系统桨距角控制装置,包括:
风机转速非线性微分方程构建模块:用于根据预获取的风力发电系统的数学模型,构建获取风机转速非线性微分方程;
风能利用系数估计模块:用于根据所述风机转速非线性微分方程,对风力发电系统中未知量风能利用系数进行估计,获取估计参数;
桨距角控制模块:用于根据预构建的风能利用系数的非线性模型,利用所述估计参数,求解获取桨距角,并将所述桨距角反馈给风力发电系统,实现风力发电系统桨距角控制。
第三方面,本发明提供一种系统,包括处理器及存储介质;
所述存储介质用于存储指令;
所述处理器用于根据所述指令进行操作以执行根据第一方面任一项所述方法的步骤。
第四方面,本发明提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现第一方面任一项所述方法的步骤。
与现有技术相比,本发明的有益效果是:
本发明提供的风电系统桨距角控制方法,基于风力发电系统的数学模型,构建获取风机转速非线性微分方程,通过直接对风能利用系数进行参数估计,求解获取桨距角,从而达到对桨距角进行间接估计的效果,使估计参数踪上真实值,在控制效果上更为精确,相对于传统的PID控制,不仅能够维持功率稳定,而且能够提高风力发电系统的稳定性。传统的变桨距控制大多对桨距角进行参数估计,但是桨距角作为风机的物理参数,其变化曲线并不能直接反应风机性能的优良,相反,风能利用系数是评估风机性能的一个重要指标,本发明通过对风能利用系数进行估计,其变化曲线能够直接反应当前风力发电系统的性能。另外,本发明还能够通过调节控制器参数,优化风能利用系数,使其达到最优值。
附图说明
图1是本发明实施例提供的风电系统桨距角控制方法流程图;
图2是本发明实施例提供的风轮转速的变化示意图;
图3是本发明实施例提供的额定风速9m/s时PID控制下风力发电系统功率和风能利用系数的变化示意图;
图4是本发明实施例提供的额定风速9m/s时自适应控制下风力发电系统功率和风能利用系数的变化示意图;
图5是本发明实施例提供的额定风速11m/s时PID控制下风力发电系统功率和风能利用系数的变化示意图;
图6是本发明实施例提供的额定风速11m/s时自适应控制下风力发电系统功率和风能利用系数的变化示意图。
具体实施方式
下面结合具体实施方式对本专利的技术方案作进一步详细地说明。
下面详细描述本专利的实施例,实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本专利,而不能理解为对本专利的限制。在不冲突的情况下,本申请实施例以及实施例中的技术特征可以相互组合。
实施例1
图1是本发明实施例一提供的一种风电系统桨距角控制方法流程图,本流程图仅仅示出了本实施例方法的逻辑顺序,在互不冲突的前提下,在本发明其它可能的实施例中,可以以不同于图1所示的顺序完成所示出或描述的步骤。
本实施例提供的风电系统桨距角控制方法可应用于终端,可以由风电系统桨距角控制装置来执行,该装置可以由软件和/或硬件的方式实现,该装置可以集成在终端中,例如:任一具备通信功能的平板电脑或计算机设备。参见图1,本实施例的方法具体包括如下步骤:
步骤一:根据预获取的风力发电系统的数学模型,构建获取风机转速非线性微分方程;
风力发电系统的数学模型的表达式如公式(1)所示:
(1)
公式(1)中,为风轮转动惯量,/>为风机转速,/>为风机转速的导数,/>为风机机械转矩,/>为齿轮比,/>为电机转矩,/>为空气密度,/>为叶轮长度,/>为风速,/>为叶尖速比,为风能利用系数。
风机转速非线性微分方程的表达式如公式(2)所示:
(2)
公式(2)中,为风机转速的导数,/>为常数项,/>,/>为空气密度,/>为叶轮长度,/>为风轮转动惯量,/>为风速,/>为叶尖速比,/>为风能利用系数,/>为齿轮比,为电机转矩。
步骤二:根据风机转速非线性微分方程,对风力发电系统中未知量风能利用系数进行估计,获取估计参数;
风能利用系数为高度非线性方程,无法将桨距角线性化,若要通过估计桨距角的方法对风力发电系统进行控制,必须对风能利用系数进行线性化处理。为避免对桨距角进行估计时造成估计参数的数目增加,本实施例中,直接对风能利用系数进行估计。
根据风机转速非线性微分方程,对风力发电系统中未知量风能利用系数进行估计,获取估计参数括如下步骤:
步骤A:引入估计参数和风力发电系统的控制律,对风机转速非线性微分方程进行改写,获取风机转速非线性微分方程改写模型;
风机转速非线性微分方程改写模型的表达式如公式(3)所示:
(3)
公式(3)中,为风机转速的导数,/>为常数项,/>,/>为空气密度,/>为叶轮长度,/>为风轮转动惯量,/>为风速,/>为叶尖速比,/>为估计参数,/>为齿轮比,/>为电机转矩,/>为风力发电系统的控制律。
步骤B:定义转速误差和估计误差,并构建获取关于转速误差和估计误差的李雅普诺夫函数;
关于转速误差和估计误差的李雅普诺夫函数的表达式如公式(4)所示:
(4)
公式(4)中,为关于转速误差和估计误差的李雅普诺夫函数,/>为转速误差,/>,/>为额定转速,/>为风机转速,/>为估计误差,/>,/>为风能利用系数,/>为估计参数。
步骤C:根据风机转速非线性微分方程改写模型和李雅普诺夫函数,基于李雅普诺夫第二定理,结合等价确定性原则,求解获取风力发电系统的控制律;
①对李雅普诺夫函数进行求导,获取李雅普诺夫导函数:
其中,为李雅普诺夫导函数,/>为转速误差,/>为转速误差的导数,/>为估计误差,为估计误差的导数,/>为额定转速的导数。
②基于李雅普诺夫第二定理,令,结合等价确定性原则,得到:
进而得到风力发电系统的控制律,风力发电系统的控制律的表达式如公式(5)所示:
(5)
公式(5)中,为任意大于零的常数。
步骤D:根据风力发电系统的控制律,求解获取估计参数;
(ⅰ)将风力发电系统的控制律代入李雅普诺夫导函数,得到:
其中,为估计参数的导数。
(ⅱ)为使成立,令:
进而得到:
(ⅲ)对估计参数的导数求积分,得到估计参数,估计参数的表达式如公式(6)所示:
(6)
公式(6)中,为风机工作时间。
步骤三:根据预构建的风能利用系数的非线性模型,利用估计参数,求解获取桨距角,并将桨距角反馈给风力发电系统,实现风力发电系统桨距角控制;
风能利用系数的非线性模型的表达式如公式(7)所示:
(7)
公式(7)中,为风能利用系数,/>为叶尖速比,/>为桨距角,/>为关于桨距角/>和叶尖速比/>的函数,/>为转速误差,/>为风机转速,/>为叶轮长度,/>为风速。
根据预构建的风能利用系数的非线性模型,利用估计参数,求解获取桨距角包括如下步骤:
步骤a:将风能利用系数的非线性模型的表达式中的风能利用系数用求解获取的估计参数/>代替,获取仅包含桨距角/>一个未知数的非线性方程组;
步骤b:对仅包含桨距角一个未知数的非线性方程组进行求解,获取桨距角/>。
为验证本实施例提供的风电系统桨距角控制方法的有效性,将本实施例提供的风电系统桨距角控制方法在Matlab/simulink平台进行建模仿真,相关参数如表1所示。
表1 风力发电机组参数
本实施例提供的风电系统桨距角控制方法,基于转速的功率跟踪控制,因此,当转速超过额定转速时,若系统没有功率反馈结构,功率会持续下降。如图2所示,为风轮转速的变化示意图,如图3所示,为PID控制下风力发电系统功率和风能利用系数的变化示意图,如图4所示,为自适应控制下风力发电系统功率和风能利用系数的变化示意图,由图2-4可知,当转速超过额定转速时,自适应控制器仍然能够稳定功率,且在临界时刻拥有较好的动态特性。
为验证不同风速情况下,自适应控制器对功率的控制效果,将表1中额定风速设置为11m/s。如图5所示,为额定风速11m/s时PID控制下风力发电系统功率和风能利用系数的变化示意图,如图6所示,为额定风速11m/s时自适应控制下风力发电系统功率和风能利用系数的变化示意图,由图5-6可知,自适应控制器仍然能对功率进行有效控制。
本实施例提供的风电系统桨距角控制方法,基于风力发电系统的数学模型,构建获取风机转速非线性微分方程,通过直接对风能利用系数进行参数估计,求解获取桨距角,从而达到对桨距角进行间接估计的效果,使估计参数踪上真实值,在控制效果上更为精确,相对于传统的PID控制,不仅能够维持功率稳定,而且能够提高风力发电系统的稳定性。传统的变桨距控制大多对桨距角进行参数估计,但是桨距角作为风机的物理参数,其变化曲线并不能直接反应风机性能的优良,相反,风能利用系数是评估风机性能的一个重要指标,本发明通过对风能利用系数进行估计,其变化曲线能够直接反应当前风力发电系统的性能。另外,本发明还能够通过调节控制器参数,优化风能利用系数,使其达到最优值。
实施例2
本实施例提供一种风电系统桨距角控制装置,包括:
风机转速非线性微分方程构建模块:用于根据预获取的风力发电系统的数学模型,构建获取风机转速非线性微分方程;
风能利用系数估计模块:用于根据风机转速非线性微分方程,对风力发电系统中未知量风能利用系数进行估计,获取估计参数;
桨距角控制模块:用于根据预构建的风能利用系数的非线性模型,利用估计参数,求解获取桨距角,并将桨距角反馈给风力发电系统,实现风力发电系统桨距角控制。
本发明实施例所提供的风电系统桨距角控制装置可执行本发明任意实施例所提供的风电系统桨距角控制方法,具备执行方法相应的功能模块和有益效果。
实施例3
本实施例提供一种系统,包括处理器及存储介质;
存储介质用于存储指令;
处理器用于根据指令进行操作以执行实施例一中方法的步骤。
实施例4
本实施例提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现实施例一中方法的步骤。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (7)
1.一种风电系统桨距角控制方法,其特征在于,包括:
根据预获取的风力发电系统的数学模型,构建获取风机转速非线性微分方程;
根据所述风机转速非线性微分方程,对风力发电系统中未知量风能利用系数进行估计,获取估计参数;
根据预构建的风能利用系数的非线性模型,利用所述估计参数,求解获取桨距角,并将所述桨距角反馈给风力发电系统,实现风力发电系统桨距角控制;
其中,所述风力发电系统的数学模型的表达式如公式(1)所示:
(1)
公式(1)中,为风轮转动惯量,/>为风机转速,/>为风机转速的导数,/>为风机机械转矩,/>为齿轮比,/>为电机转矩,/>为空气密度,/>为叶轮长度,/>为风速,/>为叶尖速比,/>为风能利用系数;
所述风机转速非线性微分方程的表达式如公式(2)所示:
(2)
公式(2)中,为风机转速的导数,/>为常数项,/>,/>为空气密度,/>为叶轮长度,/>为风轮转动惯量,/>为风速,/>为叶尖速比,/>为风能利用系数,/>为齿轮比,/>为电机转矩;
所述风能利用系数的非线性模型的表达式如公式(7)所示:
(7)
公式(7)中,为风能利用系数,/>为叶尖速比,/>为桨距角,/>为关于桨距角/>和叶尖速比/>的函数,/>为转速误差,/>为风机转速,/>为叶轮长度,/>为风速。
2.根据权利要求1所述的风电系统桨距角控制方法,其特征在于,根据所述风机转速非线性微分方程,对风力发电系统中未知量风能利用系数进行估计,获取估计参数包括:
引入估计参数和风力发电系统的控制律,对所述风机转速非线性微分方程进行改写,获取风机转速非线性微分方程改写模型;
定义转速误差和估计误差,并构建获取关于转速误差和估计误差的李雅普诺夫函数;
根据所述风机转速非线性微分方程改写模型和李雅普诺夫函数,基于李雅普诺夫第二定理,结合等价确定性原则,求解获取风力发电系统的控制律;
根据所述风力发电系统的控制律,求解获取估计参数。
3.根据权利要求2所述的风电系统桨距角控制方法,其特征在于,所述风机转速非线性微分方程改写模型的表达式如公式(3)所示:
(3)
公式(3)中,为风机转速的导数,/>为常数项,/>,/>为空气密度,/>为叶轮长度,/>为风轮转动惯量,/>为风速,/>为叶尖速比,/>为估计参数,/>为齿轮比,/>为电机转矩,/>为风力发电系统的控制律;
所述关于转速误差和估计误差的李雅普诺夫函数的表达式如公式(4)所示:
(4)
公式(4)中,为关于转速误差和估计误差的李雅普诺夫函数,/>为转速误差,,/>为额定转速,/>为风机转速,/>为估计误差,/>,/>为风能利用系数,/>为估计参数;
所述风力发电系统的控制律的表达式如公式(5)所示:
(5)
公式(5)中,为额定转速的导数,/>为任意大于零的常数;
所述估计参数的表达式如公式(6)所示:
(6)
公式(6)中,为风机工作时间。
4.根据权利要求1所述的风电系统桨距角控制方法,其特征在于,根据预构建的风能利用系数的非线性模型,利用所述估计参数,求解获取桨距角包括:
将所述风能利用系数的非线性模型的表达式中的风能利用系数用求解获取的估计参数/>代替,获取仅包含桨距角/>一个未知数的非线性方程组;
对所述仅包含桨距角一个未知数的非线性方程组进行求解,获取桨距角/>。
5.一种风电系统桨距角控制装置,其特征在于,包括:
风机转速非线性微分方程构建模块:用于根据预获取的风力发电系统的数学模型,构建获取风机转速非线性微分方程;
风能利用系数估计模块:用于根据所述风机转速非线性微分方程,对风力发电系统中未知量风能利用系数进行估计,获取估计参数;
桨距角控制模块:用于根据预构建的风能利用系数的非线性模型,利用所述估计参数,求解获取桨距角,并将所述桨距角反馈给风力发电系统,实现风力发电系统桨距角控制;
其中,所述风力发电系统的数学模型的表达式如公式(1)所示:
(1)
公式(1)中,为风轮转动惯量,/>为风机转速,/>为风机转速的导数,/>为风机机械转矩,/>为齿轮比,/>为电机转矩,/>为空气密度,/>为叶轮长度,/>为风速,/>为叶尖速比,/>为风能利用系数;
所述风机转速非线性微分方程的表达式如公式(2)所示:
(2)
公式(2)中,为风机转速的导数,/>为常数项,/>,/>为空气密度,/>为叶轮长度,/>为风轮转动惯量,/>为风速,/>为叶尖速比,/>为风能利用系数,/>为齿轮比,/>为电机转矩;
所述风能利用系数的非线性模型的表达式如公式(7)所示:
(7)
公式(7)中,为风能利用系数,/>为叶尖速比,/>为桨距角,/>为关于桨距角/>和叶尖速比/>的函数,/>为转速误差,/>为风机转速,/>为叶轮长度,/>为风速。
6.一种系统,其特征在于,包括处理器及存储介质;
所述存储介质用于存储指令;
所述处理器用于根据所述指令进行操作以执行根据权利要求1~4任一项所述方法的步骤。
7.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现权利要求1~4任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310654046.4A CN116378898B (zh) | 2023-06-05 | 2023-06-05 | 一种风电系统桨距角控制方法、装置、系统及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310654046.4A CN116378898B (zh) | 2023-06-05 | 2023-06-05 | 一种风电系统桨距角控制方法、装置、系统及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116378898A CN116378898A (zh) | 2023-07-04 |
CN116378898B true CN116378898B (zh) | 2023-09-12 |
Family
ID=86977284
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310654046.4A Active CN116378898B (zh) | 2023-06-05 | 2023-06-05 | 一种风电系统桨距角控制方法、装置、系统及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116378898B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103291543A (zh) * | 2013-06-20 | 2013-09-11 | 上海电力学院 | 基于滑模控制理论的设计风机变桨距控制器方法 |
CN110454329A (zh) * | 2019-09-04 | 2019-11-15 | 风脉能源(武汉)股份有限公司 | 一种风电机组桨距角控制方法 |
CN115573857A (zh) * | 2022-05-27 | 2023-01-06 | 北京华能新锐控制技术有限公司 | 基于时滞adrc的风电机组变桨距控制方法及控制装置 |
-
2023
- 2023-06-05 CN CN202310654046.4A patent/CN116378898B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103291543A (zh) * | 2013-06-20 | 2013-09-11 | 上海电力学院 | 基于滑模控制理论的设计风机变桨距控制器方法 |
CN110454329A (zh) * | 2019-09-04 | 2019-11-15 | 风脉能源(武汉)股份有限公司 | 一种风电机组桨距角控制方法 |
CN115573857A (zh) * | 2022-05-27 | 2023-01-06 | 北京华能新锐控制技术有限公司 | 基于时滞adrc的风电机组变桨距控制方法及控制装置 |
Non-Patent Citations (1)
Title |
---|
大型风机的独立变桨控制方法;韩兵 等;《电力系统保护与控制》;第第44卷卷(第第2期期);1-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116378898A (zh) | 2023-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102828910B (zh) | 风力发电机组的变桨控制方法、装置和风力发电机组 | |
CN105649875B (zh) | 风力发电机组的变桨控制方法及装置 | |
CN107947228B (zh) | 基于Markov理论的含风电电力系统随机稳定性分析方法 | |
CN113931793B (zh) | 一种串列式双风轮风电机组仿真的变桨环节等值方法及系统 | |
CN113883008B (zh) | 一种可抑制多扰动因素的风机模糊自适应变桨距控制方法 | |
CN113048017B (zh) | 基于内模控制的风电机组有功功率控制优化方法及系统 | |
CN104314755A (zh) | 基于ipso的dfig变桨距ladrc方法及系统 | |
CN106294959B (zh) | 模型参考自适应控制与双馈风机降阶模型相结合的建模仿真方法 | |
CN108757306B (zh) | 一种风力发电机控制方法、装置以及风力发电机 | |
CN109488526B (zh) | 基于比例-极限学习机稳态估计的变桨距控制方法 | |
CN108614416B (zh) | 光伏虚拟同步发电机控制参数的整定方法及系统 | |
CN116378898B (zh) | 一种风电系统桨距角控制方法、装置、系统及存储介质 | |
CN116131288A (zh) | 一种考虑风光波动的综合能源系统频率控制方法及系统 | |
CN113987802B (zh) | 一种水泵水轮机特性系数的计算方法及系统 | |
CN103199524A (zh) | 一种适应多种分布式电源接入的潮流计算方法 | |
Serhoud et al. | Sensorless sliding power control of doubly fed induction wind generator based on MRAS observer | |
CN105986962B (zh) | 一种风力发电机组的最大风能捕获方法 | |
CN112486019A (zh) | 不确定风力发电机系统的最大功率跟踪模糊控制方法 | |
CN111786395A (zh) | 一种风电机组出力调节方法及装置 | |
CN105320793A (zh) | 一种从动力学角度评估风电机组变浆控制模型的方法 | |
CN115528954B (zh) | 双馈风力发电系统控制参数设置方法和系统 | |
CN109185055B (zh) | 风力发电机组快速起动方法、装置及计算机可读存储介质 | |
CN118100326B (zh) | 基于数据驱动抑制风电场功率波动的有功功率控制方法 | |
CN111082454B (zh) | 一种用于获取电力系统的振荡风险量化指标的方法及系统 | |
CN117145697A (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 |