CN104343628B - A kind of Wind turbines variable pitch control method containing dead-zone compensator - Google Patents
A kind of Wind turbines variable pitch control method containing dead-zone compensator Download PDFInfo
- Publication number
- CN104343628B CN104343628B CN201410432782.6A CN201410432782A CN104343628B CN 104343628 B CN104343628 B CN 104343628B CN 201410432782 A CN201410432782 A CN 201410432782A CN 104343628 B CN104343628 B CN 104343628B
- Authority
- CN
- China
- Prior art keywords
- pitch
- wind turbine
- dead zone
- wind
- output
- 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.)
- Expired - Fee Related
Links
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
- 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/70—Type of control algorithm
- F05B2270/709—Type of control algorithm with neural networks
-
- 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)
- Wind Motors (AREA)
Abstract
本发明公开了一种含有死区补偿器的风电机组变桨距控制方法,属于风电机组控制系统领域,其包括建立风电机组的数学模型,得到变桨距控制器输出的桨距角期望值与变桨距执行机构输出的桨距角实际值的数学表达式、设计变桨距执行机构的死区补偿器,得到桨距角补偿值的数学表达式、测量风轮机周边的风速和风力发电机输出的有功功率、计算出桨距角期望值并分别传送至变桨距执行机构和变桨距执行机构的死区补偿器、计算出桨距角实际值和桨距角补偿值并传送给风轮机的内部调相器、风轮机的内部调相器对风轮机进行调节。该方法很好的解决了风电机组运行中,变桨距执行机构死区对整个系统的干扰和影响,提升风电机组的运行状态,保证功率输出的稳定性。
The invention discloses a pitch control method of a wind turbine with a dead zone compensator, which belongs to the field of wind turbine control systems, and includes establishing a mathematical model of the wind turbine and obtaining the expected value of the pitch angle output by the pitch controller and the variable pitch angle. The mathematical expression of the actual value of the pitch angle output by the pitch actuator, the design of the dead zone compensator of the pitch actuator, the mathematical expression of the pitch angle compensation value, the measurement of the wind speed around the wind turbine and the output of the wind turbine The active power of the pitch angle is calculated and sent to the pitch actuator and the dead zone compensator of the pitch actuator respectively, the actual value of the pitch angle and the compensation value of the pitch angle are calculated and sent to the wind turbine The internal phase regulator, the internal phase regulator of the wind turbine, regulates the wind turbine. This method well solves the interference and influence of the dead zone of the variable pitch actuator on the entire system during the operation of the wind turbine, improves the operating state of the wind turbine, and ensures the stability of power output.
Description
技术领域technical field
本发明属于风电机组控制系统领域,具体涉及一种含有死区补偿器的风电机组变桨距控制方法。The invention belongs to the field of wind turbine control systems, and in particular relates to a pitch control method for wind turbines with a dead zone compensator.
背景技术Background technique
风能,作为一种清洁的可再生能源,越来越受到世界各国的重视,并得到快速的发展。对于风能的利用,主要依靠于风电机组来进行能量转换。风电机组控制系统是机组正常运行的核心,其控制技术是风电机组的关键技术之一,与风电机组的其他部分关系密切,其精确的控制、完善的功能将直接影响整个机组的安全与效率。Wind energy, as a clean and renewable energy, has been paid more and more attention by countries all over the world, and has developed rapidly. For the utilization of wind energy, it mainly relies on wind turbines for energy conversion. The wind turbine control system is the core of the normal operation of the wind turbine. Its control technology is one of the key technologies of the wind turbine and is closely related to other parts of the wind turbine. Its precise control and perfect functions will directly affect the safety and efficiency of the entire wind turbine.
变桨距控制技术作为风电机组的一种主流控制技术已经得到了广泛的应用,然而,由于风能具有很强的随机性、间歇性和不可调控等性质,再加上变桨距执行机构存在死区以及其他诸多不确定因素,使得变桨距控制系统具有参数非线性、参数时变性、滞后性等特点,特别是变桨距执行机构死区所带来的影响,使得变桨距系统的控制调节存在偏差,造成风电机组输出功率不稳定。Pitch control technology has been widely used as a mainstream control technology of wind turbines. However, due to the strong randomness, intermittent and uncontrollable properties of wind energy, and the deadlock of the pitch actuator zone and many other uncertain factors make the pitch control system have the characteristics of parameter nonlinearity, parameter time-varying, hysteresis, etc., especially the influence brought by the dead zone of the pitch actuator, making the control of the pitch system There is a deviation in the adjustment, which causes the output power of the wind turbine to be unstable.
发明内容Contents of the invention
针对现有技术存在的不足,本发明提供一种含有死区补偿器的风电机组变桨距控制方法。Aiming at the deficiencies in the prior art, the present invention provides a pitch control method for a wind turbine with a dead zone compensator.
本发明的技术方案:Technical scheme of the present invention:
步骤1:建立风电机组的数学模型,得到变桨距控制器输出的桨距角期望值βr与变桨距执行机构输出的桨距角实际值β的数学表达式:Step 1: Establish a mathematical model of the wind turbine, and obtain the mathematical expressions of the expected value of the pitch angle β r output by the pitch controller and the actual value of the pitch angle β output by the pitch actuator:
式中:τ为时间常数;s为复频率,是傅里叶变换的补充。Where: τ is the time constant; s is the complex frequency, which is the complement of Fourier transform.
同时得到风轮机捕获风功率能力的数学表达式:At the same time, the mathematical expression of the wind turbine's ability to capture wind power is obtained:
式中:Cp(λ,β)为风能利用系数。Where: C p (λ, β) is the utilization coefficient of wind energy.
步骤2:根据步骤1得到的变桨距控制器输出的桨距角期望值βr与变桨距执行机构输出的桨距角实际值β的数学表达式,设计变桨距执行机构的死区补偿器,得到桨距角补偿值Δβr的数学表达式,按如下步骤进行:Step 2: According to the mathematical expression of the pitch angle expected value β r output by the pitch controller and the actual pitch angle β output by the pitch actuator obtained in step 1, design the dead zone compensation of the pitch actuator device to obtain the mathematical expression of the pitch angle compensation value Δβ r , proceed as follows:
步骤2.1:建立变桨距执行机构死区的数学模型;Step 2.1: Establish a mathematical model of the dead zone of the variable pitch actuator;
本发明中所使用的变桨距执行机构死区的数学模型如下所示:The mathematical model of the dead zone of the variable pitch actuator used in the present invention is as follows:
式中:D′(u)表示关于变桨距执行机构死区的函数表达式;br为离开变桨距执行机构死区上限风速大小,m/s;bl为进入变桨距执行机构死区下限风速大小,m/s;u为变桨距执行机构死区的输入;gr(u)为当u≥br时变桨距执行机构死区函数表达式;gl(u)为当u≤bl时变桨距执行机构死区函数表达式;In the formula: D′(u) represents the functional expression about the dead zone of the pitch actuator; b r is the upper limit of the wind speed leaving the dead zone of the pitch actuator, m/s; b l is the wind speed entering the pitch actuator The wind speed of the lower limit of the dead zone, m/s; u is the input of the dead zone of the variable pitch actuator; g r (u) is the dead zone function expression of the variable pitch actuator when u≥b r ; g l (u) is the dead zone function expression of the variable pitch actuator when u≤b l ;
步骤2.2:对非线性动态函数f(x)的估计;Step 2.2: Estimation of the nonlinear dynamic function f(x);
f(x)=W1 *Tσ(V1 TX1)+ε1 (4)f(x)=W 1 *T σ(V 1 T X 1 )+ε 1 (4)
式中:X1=[x1 x2…xn 1]T∈Rn+1;V1为第一神经网络输入层与隐含层之间的权矩阵;W1 *为第一神经网络隐含层与输出层之间的理想权矩阵;ε1为第一神经网络重构误差;σ为激活函数,具体形式为 In the formula: X 1 =[x 1 x 2 …x n 1] T ∈ R n+1 ; V 1 is the weight matrix between the input layer and the hidden layer of the first neural network; W 1 * is the weight matrix of the first neural network The ideal weight matrix between the hidden layer and the output layer; ε 1 is the reconstruction error of the first neural network; σ is the activation function, and the specific form is
考虑到阈值作用,本发明的激活函数选为其中,本发明中有l为隐含层节点数,z代表整个系统的广义输出;而实际的第一神经网络的输出为Considering the threshold effect, the activation function of the present invention is selected as Among them, the present invention has l is the number of hidden layer nodes, z represents the generalized output of the whole system; and the actual output of the first neural network is
式中,W1为第一神经网络隐含层与输出层之间的实际权矩阵;则第一神经网络估计误差可以表示为In the formula, W 1 is the actual weight matrix between the hidden layer and the output layer of the first neural network; then the estimation error of the first neural network can be expressed as
式中:为神经网络权值估计误差;In the formula: Estimated error for neural network weights;
步骤2.3:设计变桨距执行机构死区的补偿器,得到桨距角补偿值Δβr的数学表达式;Step 2.3: Design the compensator for the dead zone of the variable pitch actuator, and obtain the mathematical expression of the pitch angle compensation value Δβr ;
在补偿器设计中利用第二神经网络的输出抵消由变桨距执行机构死区引起的系统调节误差η(u),其中第二神经网络输入为uc;Utilizing the output of the second neural network in compensator design Offset the system adjustment error η(u) caused by the dead zone of the pitch actuator, where the second neural network input is u c ;
利用第二神经网络得到桨距角补偿值Δβr的广义数学表达式为The generalized mathematical expression of the pitch angle compensation value Δβ r obtained by using the second neural network is
式中:W2为第二神经网络隐含层与输出层之间的实际权矩阵;X2=[x1 x2…xn 1]T∈Rn+1;V2为第二神经网络输入层与隐含层之间的权矩阵;ε2为第二神经网络重构误差。In the formula: W 2 is the actual weight matrix between the hidden layer and the output layer of the second neural network; X 2 =[x 1 x 2 …x n 1] T ∈ R n+1 ; V 2 is the second neural network The weight matrix between the input layer and the hidden layer; ε 2 is the reconstruction error of the second neural network.
步骤3:利用风速传感器测量风轮机周边的风速v,利用霍尔电压、电流传感器分别测得风力发电机输出的电压和电流,两者相乘得到风力发电机输出的有功功率P。Step 3: Use the wind speed sensor to measure the wind speed v around the wind turbine, use the Hall voltage and current sensors to measure the output voltage and current of the wind turbine, and multiply the two to obtain the active power P output by the wind turbine.
步骤4:变桨距控制器计算出桨距角期望值βr并分别传送至变桨距执行机构和变桨距执行机构的死区补偿器;Step 4: The pitch controller calculates the pitch angle expected value β r and sends them to the pitch actuator and the dead zone compensator of the pitch actuator respectively;
步骤5:变桨距执行机构计算出桨距角实际值β并传送给风轮机的内部调相器,同时,变桨距执行机构的死区补偿器计算出桨距角补偿值Δβr并传送给风轮机的内部调相器;Step 5: The pitch actuator calculates the actual value of the pitch angle β and transmits it to the internal phase modulator of the wind turbine. At the same time, the dead zone compensator of the pitch actuator calculates the pitch angle compensation value Δβ r and transmits it to Internal phase modulators for wind turbines;
步骤6:根据接收到的桨距角值,即桨距角实际值β和桨距角补偿值Δβr,风轮机的内部调相器对风轮机进行调节。Step 6: According to the received pitch angle value, that is, the actual pitch angle value β and the pitch angle compensation value Δβ r , the internal phase regulator of the wind turbine adjusts the wind turbine.
有益效果:Beneficial effect:
变桨距执行机构中存在死区并且风电机组受其很大影响的风险是客观存在的,因此构建一套完善的专用于风电机组的含有死区补偿器的变桨距控制方法具有迫切性和必要性。本发明的变桨距控制方法解决了以下几个存在的问题:There is an objective risk that there is a dead zone in the pitch actuator and the wind turbine is greatly affected by it. Therefore, it is urgent and necessary to construct a complete pitch control method with a dead zone compensator for wind turbines. necessity. The pitch control method of the present invention solves the following problems:
(1)考虑到风电机组中的数据量具有模糊性、随机性、不确定性和冗余性等特点,利用本发明的变桨距控制方法,成功克服了目前现有技术只能将变桨距执行机构死区看作理想状态而无法去实际解决的缺点;(1) Considering that the amount of data in wind turbines has the characteristics of ambiguity, randomness, uncertainty, and redundancy, the pitch control method of the present invention successfully overcomes the current existing technology that can only change pitch The shortcoming that the distance from the dead zone of the actuator is regarded as an ideal state and cannot be practically solved;
(2)提高了风电机组的整体性能及变桨距系统的控制精度,基本消除了变桨距执行机构死区对风电机组的影响。专门用于风电机组的含有死区补偿器的变桨距控制方法的使用,很好的解决了风电机组运行中,变桨距执行机构死区对整个系统的干扰和影响,能够很好的提升风电机组的运行状态,保证功率输出的稳定性。(2) The overall performance of the wind turbine and the control accuracy of the pitch control system are improved, and the influence of the dead zone of the pitch actuator on the wind turbine is basically eliminated. The use of the variable pitch control method including the dead zone compensator specially used for wind turbines solves the interference and influence of the dead zone of the variable pitch actuator on the entire system during the operation of the wind turbine, and can be well improved. The running status of wind turbines ensures the stability of power output.
附图说明Description of drawings
图1为本发明一种实施方式的风电机组变桨距控制结构示意图;Fig. 1 is a schematic diagram of a pitch control structure of a wind turbine in an embodiment of the present invention;
图2为本发明一种实施方式的变桨距执行机构的死区区域示意图;Fig. 2 is a schematic diagram of a dead zone area of a pitch actuator according to an embodiment of the present invention;
图3为本发明一种实施方式的含有死区补偿器的风电机组变桨距控制方法流程图;Fig. 3 is a flow chart of a pitch control method for wind turbines including a dead zone compensator according to an embodiment of the present invention;
图4为本发明一种实施方式的含有死区补偿器的风电机组变桨距控制结构示意图;Fig. 4 is a schematic diagram of a pitch control structure of a wind turbine with a dead zone compensator according to an embodiment of the present invention;
图5为本发明一种实施方式的风轮机转速仿真图;Fig. 5 is a wind turbine speed simulation diagram of an embodiment of the present invention;
图6为本发明一种实施方式的风力发电机输出功率仿真图。Fig. 6 is a simulation diagram of output power of a wind power generator according to an embodiment of the present invention.
具体实施方式detailed description
为详细说明本发明的技术内容、所达成目的及效果,下面结合附图及具体实施方式作进一步详细说明。In order to describe the technical content, purpose and effect of the present invention in detail, the following will be further described in detail in conjunction with the accompanying drawings and specific implementation methods.
如附图1所示,为一个未加死区补偿器的风电机组变桨距控制结构示意图,包括风轮机、传动系统、风力发电机和变桨距系统,变桨距系统又包括变桨距控制器和变桨距执行机构。风轮机是风电机组的原动件,用于捕获风能;传动系统将风轮机捕获的风能传递给风力发电机;风力发电机经过电磁转换,输出有功功率,变桨距系统则是通过调节风轮机的桨距角来调节风能捕获大小,从而来保证风电机组能够输出恒定的功率值,符合系统要求。其中,变桨距控制器根据接收的风力发电机输出的有功功率P和风轮机周边的风速v,利用风电机组的数学模型计算出应该调节的桨距角期望值βr,并将该值传送给变桨距执行机构来执行。然而,由于变桨距执行机构存在死区,如附图2所示,在风速v高于额定风速v1小于切出风速v2时,变桨距系统开始工作,也就是附图2中的区域2。在这段区域内,由于变桨距执行机构存在死区,导致最终变桨距执行机构执行时所传送给风轮机的桨距角值并不是βr,而是产生一定偏差的桨距角值实际值β,也就是变桨距执行机构实际输出桨距角值,所产生的偏差量大小为β′=βr-β。由于存在一定偏差,故而导致风轮机的内部调相器并未将风轮机调节到所应调节的部位,从而影响了风力发电机的有功功率P的输出。As shown in Figure 1, it is a schematic diagram of the pitch control structure of a wind turbine without a dead zone compensator, including a wind turbine, a transmission system, a wind generator, and a pitch control system, and the pitch control system includes a pitch control system. Controller and pitch actuator. The wind turbine is the original moving part of the wind turbine, which is used to capture wind energy; the transmission system transmits the wind energy captured by the wind turbine to the wind generator; the wind generator outputs active power through electromagnetic conversion, and the pitch system adjusts the wind turbine The pitch angle is used to adjust the wind energy capture, so as to ensure that the wind turbine can output a constant power value, which meets the system requirements. Among them, the pitch controller calculates the expected value of the pitch angle β r that should be adjusted according to the received active power P output by the wind turbine and the wind speed v around the wind turbine using the mathematical model of the wind turbine, and sends this value to the variable Pitch actuator to perform. However, due to the dead zone of the pitch actuator, as shown in Figure 2 , when the wind speed v is higher than the rated wind speed v1 and less than the cut-out wind speed v2, the pitch system starts to work, that is, in Figure 2 area 2. In this area, due to the dead zone of the pitch actuator, the pitch angle value transmitted to the wind turbine when the final pitch actuator is executed is not β r , but a pitch angle value with a certain deviation The actual value β, that is, the actual output pitch angle value of the variable pitch actuator, produces a deviation of β′=β r −β. Due to a certain deviation, the internal phase regulator of the wind turbine does not adjust the wind turbine to the position that should be adjusted, thereby affecting the output of the active power P of the wind turbine.
基于此,在本发明中通过设计一个变桨距执行机构的死区补偿器产生一个桨距角补偿值Δβr来抵消桨距角偏差量β′,使得最终传送到风轮机内部调相器的桨距角值尽量达到桨距角期望值βr。Based on this, in the present invention, a pitch angle compensation value Δβ r is generated by designing a dead zone compensator of a variable pitch actuator to offset the pitch angle deviation β', so that the final value transmitted to the internal phase modulator of the wind turbine The pitch angle value should reach the expected value β r of the pitch angle as much as possible.
本实施方式的含有死区补偿器的风电机组变桨距控制方法,如图3所示,包括如下具体步骤:The wind turbine pitch control method containing the dead zone compensator in this embodiment, as shown in Figure 3, includes the following specific steps:
步骤1:建立风电机组的数学模型,得到变桨距控制器输出的桨距角期望值βr的数学表达式、变桨距执行机构输出的桨距角实际值β的数学表达式;Step 1: Establish a mathematical model of the wind turbine, obtain the mathematical expression of the pitch angle expected value β r output by the pitch controller, and the mathematical expression of the actual pitch angle β output by the pitch actuator;
步骤1.1:建立风轮机的数学模型;Step 1.1: Establish a mathematical model of the wind turbine;
风轮机是将风能转化为机械能的装置,根据贝茨理论,风轮机从风能中捕获的功率和机械转矩为:A wind turbine is a device that converts wind energy into mechanical energy. According to Bates theory, the power and mechanical torque captured by a wind turbine from wind energy are:
式中:Pr为风轮机捕获的功率,w;Tr为风轮的机械转矩,N·m;λ为叶尖速比;ρ为空气的密度,kg/m3;R为风轮的半径,m;v为风速,m/s;β为桨距角实际值,rad;CT(λ,β)为机械转矩系数;Cp(λ,β)为风能利用系数;In the formula: P r is the power captured by the wind turbine, w; T r is the mechanical torque of the wind rotor, N m; λ is the tip speed ratio; ρ is the density of the air, kg/m 3 ; R is the wind rotor The radius of , m; v is the wind speed, m/s; β is the actual value of the pitch angle, rad; C T (λ, β) is the mechanical torque coefficient; C p (λ, β) is the wind energy utilization coefficient;
CT(λ,β)与Cp(λ,β)的关系可由下列表达式表示出来:The relationship between C T (λ, β) and C p (λ, β) can be expressed by the following expression:
Cp(λ,β)=λCT(λ,β) (10)C p (λ, β) = λC T (λ, β) (10)
风能利用系数Cp(λ,β)代表了风轮机捕获风功率的能力,本发明采用如下公式进行近似:The wind energy utilization coefficient Cp (λ, β) represents the ability of the wind turbine to capture wind power, and the present invention adopts the following formula to approximate:
步骤1.2:建立风力发电机的数学模型;Step 1.2: Establish a mathematical model of the wind turbine;
根据风轮机的数学模型,可以得出风轮机捕获的风能大小。风轮机作为原动件将所捕获的风能Pr传送给传动系统,传动系统再将风能转化为机械能传送给风力发电机,最终传动系统传送给风力发电机的转矩为Tg,风力发电机再根据传送过来的机械能通过电磁转换转化为所需要的电能,为了便于计算,本实施方式建立的风力发电机的数学模型如下所示:According to the mathematical model of the wind turbine, the amount of wind energy captured by the wind turbine can be obtained. The wind turbine, as the original moving part, transmits the captured wind energy P r to the transmission system, and the transmission system converts the wind energy into mechanical energy and transmits it to the wind turbine. The final torque transmitted by the transmission system to the wind turbine is T g . Then, according to the transmitted mechanical energy, it is converted into the required electric energy through electromagnetic conversion. For the convenience of calculation, the mathematical model of the wind power generator established in this embodiment is as follows:
式中:Jr为风轮转动惯量,kg·m2;Jg为风力发电机转动惯量,kg·m2;ωr为风轮机实际转速,r/min;vi为减速比;Te为风力发电机电磁转矩,N·m;TD为风力发电机阻力矩,N·m;C1,C2,C3为阻尼系数;In the formula: J r is the moment of inertia of the wind rotor, kg m 2 ; J g is the moment of inertia of the wind turbine, kg m 2 ; ω r is the actual speed of the wind turbine, r/min; v i is the reduction ratio; T e is the electromagnetic torque of the wind turbine, N m; T D is the resistance torque of the wind turbine, N m; C 1 , C 2 , C 3 are the damping coefficients;
步骤1.3:建立变桨距控制器的数学模型,得到变桨距控制器输出的桨距角期望值βr与变桨距执行机构输出的桨距角实际值β的数学表达式;Step 1.3: establish the mathematical model of the variable pitch controller, and obtain the mathematical expression of the pitch angle expected value β r output by the variable pitch controller and the actual value β of the pitch angle output by the variable pitch actuator;
通过建立的风轮机和风力发电机的数学模型,能够得出本实施方式中风轮机所捕获的风能和风力发电机所转换输出的电能。在附图2区域2中,需要通过变桨距系统来调节桨距角从而保证风力发电机输出功率恒定。对于变桨距系统,包括变桨距控制器和变桨距执行机构两部分,本发明需要通过建立变桨距控制器的数学模型,得到变桨距控制器输出的桨距角期望值βr与变桨距执行机构输出的桨距角实际值β的数学表达式。Through the established mathematical models of the wind turbine and the wind generator, the wind energy captured by the wind turbine and the electrical energy converted and output by the wind generator in this embodiment can be obtained. In area 2 of Fig. 2, it is necessary to adjust the pitch angle through the pitch control system to ensure the constant output power of the wind turbine. For the variable pitch system, including the two parts of the pitch controller and the pitch actuator, the present invention needs to obtain the pitch angle expected value β r output by the pitch controller by establishing the mathematical model of the pitch controller and The mathematical expression of the actual value β of the pitch angle output by the pitch actuator.
本实施方式建立的变桨距控制器的数学模型如下所示:The mathematical model of the pitch controller established in this embodiment is as follows:
式中:τ为时间常数;βr为桨距角期望值,rad;where: τ is the time constant; β r is the expected value of the pitch angle, rad;
得到的变桨距控制器输出的桨距角期望值βr的数学表达式、变桨距执行机构输出的桨距角实际值β的数学表达式,为共同的一个关系式,如式(1)所示:The obtained mathematical expression of the pitch angle expected value β r output by the pitch controller and the mathematical expression of the actual pitch angle β output by the pitch actuator are a common relational expression, such as formula (1) Shown:
式中:s为复频率,是傅里叶变换的补充。In the formula: s is the complex frequency, which is the complement of Fourier transform.
步骤2:根据步骤1的结果,设计变桨距执行机构的死区补偿器,得到桨距角补偿值Δβr的数学表达式;Step 2: According to the result of step 1, design the dead zone compensator of the variable pitch actuator, and obtain the mathematical expression of the pitch angle compensation value Δβ r ;
步骤2.1:建立变桨距执行机构死区的数学模型;Step 2.1: Establish a mathematical model of the dead zone of the variable pitch actuator;
本实施方式中所使用的变桨距执行机构死区的数学模型如下所示:The mathematical model of the dead zone of the variable pitch actuator used in this embodiment is as follows:
式中:D′(u)表示关于变桨距执行机构死区的函数表达式;br为离开变桨距执行机构死区上限风速大小,m/s;bl为进入变桨距执行机构死区下限风速大小,m/s;u为变桨距执行机构死区的输入;gr(u)为当u≥br时变桨距执行机构死区函数表达式;gl(u)为当u≤bl时变桨距执行机构死区函数表达式;In the formula: D′(u) represents the functional expression about the dead zone of the pitch actuator; b r is the upper limit of the wind speed leaving the dead zone of the pitch actuator, m/s; b l is the wind speed entering the pitch actuator The wind speed of the lower limit of the dead zone, m/s; u is the input of the dead zone of the variable pitch actuator; g r (u) is the dead zone function expression of the variable pitch actuator when u≥b r ; g l (u) is the dead zone function expression of the variable pitch actuator when u≤b l ;
步骤2.2:对非线性动态函数f(x)的估计;Step 2.2: Estimation of the nonlinear dynamic function f(x);
为了设计变桨距执行机构的死区补偿器,本发明必须对被控系统进行动态估计,即针对被控系统的非线性动态函数f(x)进行估计,然后将得到的估计值输入到死区补偿器作为一个补偿参考值。根据神经网络函数近似理论,本实施方式中,通过一个简单的第一神经网络来对被控系统的非线性动态函数进行估计。假设这个第一神经网络的隐含层与输入层之间的权值随机给定后不再调节,则有:In order to design the dead zone compensator of the variable pitch actuator, the present invention must dynamically estimate the controlled system, that is, estimate the nonlinear dynamic function f(x) of the controlled system, and then input the obtained estimated value into the dead zone zone compensator as a compensation reference. According to the neural network function approximation theory, in this embodiment, a simple first neural network is used to estimate the nonlinear dynamic function of the controlled system. Assuming that the weights between the hidden layer and the input layer of the first neural network are randomly given and no longer adjusted, then:
f(x)=W1 *Tσ(V1 TX1)+ε1 (4)f(x)=W 1 *T σ(V 1 T X 1 )+ε1 (4)
式中:X1=[x1 x2…xn1]T∈Rn+1;V1为第一神经网络输入层与隐含层之间的权矩阵;W1 *为第一神经网络隐含层与输出层之间的理想权矩阵;ε1为第一神经网络重构误差;σ为激活函数,具体形式为 In the formula: X 1 =[x 1 x 2 …x n 1] T ∈ R n+1 ; V 1 is the weight matrix between the input layer and the hidden layer of the first neural network; W 1 * is the weight matrix of the first neural network The ideal weight matrix between the hidden layer and the output layer; ε 1 is the reconstruction error of the first neural network; σ is the activation function, and the specific form is
考虑到阈值作用,本实施方式中的激活函数选为其中,l为隐含层节点数,z代表整个系统的广义输出;而实际的第一神经网络的输出为Considering the threshold effect, the activation function in this embodiment is selected as in, l is the number of hidden layer nodes, z represents the generalized output of the whole system; and the actual output of the first neural network is
式中,W1为第一神经网络隐含层与输出层之间的实际权矩阵;In the formula, W 1 is the actual weight matrix between the hidden layer and the output layer of the first neural network;
则第一神经网络估计误差可以表示为Then the estimation error of the first neural network can be expressed as
式中:为神经网络权值估计误差;In the formula: Estimated error for neural network weights;
步骤2.3:设计变桨距执行机构死区的补偿器,得到桨距角补偿值Δβr的数学表达式;Step 2.3: Design the compensator for the dead zone of the variable pitch actuator, and obtain the mathematical expression of the pitch angle compensation value Δβr ;
在变桨距执行机构的死区补偿器设计中利用第二神经网络的输出抵消由变桨距执行,机构死区引起的系统调节误差η(u),其中第二神经网络输入为uc;Utilizing the output of the second neural network in the dead zone compensator design of the pitch actuator Offset the system adjustment error η(u) caused by the pitch adjustment and the dead zone of the mechanism, wherein the input of the second neural network is u c ;
利用第二神经网络得到桨距角补偿值Δβr的广义数学表达式为The generalized mathematical expression of the pitch angle compensation value Δβ r obtained by using the second neural network is
式中:W2为第二神经网络隐含层与输出层之间的实际权矩阵;X2=[x1 x2…xn 1]T∈Rn+1;V2为第二神经网络输入层与隐含层之间的权矩阵;ε2为第二神经网络重构误差。In the formula: W 2 is the actual weight matrix between the hidden layer and the output layer of the second neural network; X 2 =[x 1 x 2 …x n 1] T ∈ R n+1 ; V 2 is the second neural network The weight matrix between the input layer and the hidden layer; ε 2 is the reconstruction error of the second neural network.
如图4所示,本发明在变桨距执行机构处加了一个死区补偿器来弥补变桨距执行机构死区所带来的系统偏差。如图4所示,变桨距控制器除了要将所得出的桨距角期望值βr传送给变桨距执行机构外,还要传送给变桨距执行机构死区补偿器,与此同时,变桨距执行机构也要将其最终所输出的桨距角实际值β传送给变桨距执行机构死区补偿器,死区补偿器则通过上述内部神经网络算法以及误差和偏差设定值,最终给出所应桨距角补偿值Δβr,并将该桨距角补偿值传送给风轮机内部调相器,风轮机内部调相器则同时接收到来自变桨距执行机构死区补偿器的桨距角补偿值Δβr和变桨距执行机构的桨距角实际值β。风轮机的内部调相器再根据所接收到的桨距角值对风轮机进行调节。As shown in Figure 4, the present invention adds a dead zone compensator at the pitch actuator to compensate for the system deviation caused by the dead zone of the pitch actuator. As shown in Figure 4, the pitch controller not only transmits the obtained pitch angle expected value β r to the pitch actuator, but also transmits it to the dead zone compensator of the pitch actuator. At the same time, The pitch actuator should also transmit the actual value of the pitch angle β that it finally outputs to the dead zone compensator of the pitch actuator, and the dead zone compensator uses the above-mentioned internal neural network algorithm and the error and deviation set value, Finally, the corresponding pitch angle compensation value Δβ r is given, and the pitch angle compensation value is sent to the internal phase modulator of the wind turbine. The pitch angle compensation value Δβ r and the actual value β of the pitch angle of the pitch actuator. The internal phase modulator of the wind turbine then adjusts the wind turbine according to the received pitch angle value.
步骤3:测量现场风轮机周边的风速v和风力发电机输出的有功功率P并传送至变桨距控制器;Step 3: Measure the wind speed v around the on-site wind turbine and the active power P output by the wind turbine and send it to the pitch controller;
本实施方式在风力发电机组中加入了风速传感器和功率传感器,如附图4所示。利用风速传感器测量风轮机周边风速v,这里的风速传感器型号采用上海微川精密仪器有限公司生产的EE65风速变送器,其测量范围为0-20m/s,输出信号为4-20mA/0-10V,精度为±0.2m/s,工作环境为-25-50℃,输出信号幅值与风速大小成正比。In this embodiment, a wind speed sensor and a power sensor are added to the wind power generating set, as shown in FIG. 4 . Use the wind speed sensor to measure the wind speed v around the wind turbine. The wind speed sensor here adopts the EE65 wind speed transmitter produced by Shanghai Weichuan Precision Instrument Co., Ltd. The measurement range is 0-20m/s, and the output signal is 4-20mA/0- 10V, the accuracy is ±0.2m/s, the working environment is -25-50℃, the output signal amplitude is proportional to the wind speed.
风力发电机的功率信号经功率传感器输入到变桨距控制器,最终完成对风力发电机组的额定功率的控制。这里的功率传感器分别采用霍尔电流传感器和霍尔电压传感器检测电流和电压,根据所检测的电流电压便获得了风力发电机的输出有功功率P。本发明中,电压传感器型号为瑞士莱姆公司(LEM)生产的LV 200-AW/2,原边额定有效电流为20mA,采用±15-24V电源供电,精度为±0.5%;电流传感器型号为瑞士莱姆公司(LEM)生产的LA28-NP,原边额定有效值电流为25mA,采用±15V电源供电,精度为±0.5%,原边电流测量范围为0~±36A。The power signal of the wind turbine is input to the pitch controller through the power sensor, and finally completes the control of the rated power of the wind turbine. The power sensor here uses a Hall current sensor and a Hall voltage sensor to detect current and voltage, and the output active power P of the wind turbine is obtained according to the detected current and voltage. In the present invention, the voltage sensor model is LV 200-AW/2 produced by Swiss Lem Company (LEM), the rated effective current of the primary side is 20mA, and the ±15-24V power supply is used for power supply, and the accuracy is ±0.5%; the current sensor model is The LA28-NP produced by LEM of Switzerland has a rated RMS current of 25mA on the primary side, is powered by a ±15V power supply, and has an accuracy of ±0.5%. The measuring range of the primary side current is 0 to ±36A.
通过以上传感器测量现场风轮机周边的风速v和风力发电机输出的有功功率P并传送至变桨距控制器。最终所测得的其中一组数据及实际风力发电系统参数如下表所示:The wind speed v around the wind turbine and the active power P output by the wind generator are measured by the above sensors and sent to the pitch controller. One of the final measured data and the actual wind power system parameters are shown in the following table:
表1测试风力发电系统所设定的参数Table 1 The parameters set for testing the wind power generation system
步骤4:变桨距控制器计算出桨距角期望值βr并分别传送至变桨距执行机构和变桨距执行机构的死区补偿器;Step 4: The pitch controller calculates the pitch angle expected value β r and sends them to the pitch actuator and the dead zone compensator of the pitch actuator respectively;
在本实施方式中,利用仿真软件MATLAB/Simulink对风电机组的变桨距执行机构进行仿真。In this embodiment, simulation software MATLAB/Simulink is used to simulate the pitch-variable actuator of the wind turbine.
本实施方式的仿真过程是基于一个额定功率为600KW的变桨距风电机组上实现的。将第一神经网络和第二神经网络的结构分别选作为2-21-1(即,输入层数为2,隐含层数为21和输出层数为1)和3-11-1(即,输入层数为3,隐含层数为11和输出层数为1),权值的初始值W1和设置值W2分别为0.1和0.1,其他的在这两个网络的权重的初始值设定为0.01。激活函数设定为单极性S型函数 The simulation process of this embodiment is realized on the basis of a variable-pitch wind turbine with a rated power of 600KW. The structures of the first neural network and the second neural network are respectively selected as 2-21-1 (that is, the number of input layers is 2, the number of hidden layers is 21 and the number of output layers is 1) and 3-11-1 (that is, , the number of input layers is 3, the number of hidden layers is 11 and the number of output layers is 1), the initial value W 1 and the setting value W 2 of the weight are 0.1 and 0.1 respectively, and the other initial weights of the two networks The value is set to 0.01. The activation function is set to a unipolar sigmoid function
将步骤3中表1中的数据设定为所仿真的风电机组的参数,将步骤3中传感器所测得的数据作为输入值直接传送给变桨距控制器,则变桨距控制器利用步骤1中所建立的数学模型计算出桨距角期望值βr,本实施方式的仿真过程中,并不会把桨距角期望值βr的具体数值显示出来,而是分别传送至变桨距执行机构和变桨距执行机构的死区补偿器;Set the data in Table 1 in step 3 as the parameters of the simulated wind turbine, and directly transmit the data measured by the sensor in step 3 to the pitch controller as the input value, then the pitch controller uses the step The mathematical model established in 1 calculates the expected value of the pitch angle β r . In the simulation process of this embodiment, the specific value of the expected value of the pitch angle β r is not displayed, but is transmitted to the pitch actuator respectively and the dead zone compensator of the pitch actuator;
步骤5:变桨距执行机构计算出桨距角实际值β并传送给风轮机的内部调相器,同时,变桨距执行机构的死区补偿器计算出桨距角补偿值Δβr也传送给风轮机的内部调相器;Step 5: The pitch actuator calculates the actual value of the pitch angle β and transmits it to the internal phase modulator of the wind turbine. At the same time, the pitch angle compensation value Δβ r calculated by the dead zone compensator of the pitch actuator is also transmitted Internal phase modulators for wind turbines;
在本实施方式的MATLAB/Simulink仿真中,所搭建的仿真模型将步骤4计算得出的桨距角期望值βr代入到步骤1中所得到的变桨距控制器输出的桨距角期望值βr的数学表达式与变桨距执行机构输出的桨距角实际值β的数学表达式,得出变桨距执行机构桨距角实际值β并传送给风轮机的内部调相器,同时,变桨距执行机构的死区补偿器根据步骤2中的数学表达式利用步骤3中表1的参数值以及步骤4所得到的桨距角期望值βr计算出桨距角补偿值Δβr并传送给风轮机的内部调相器;In the MATLAB/Simulink simulation of this embodiment, the simulation model built substitutes the expected value of pitch angle β r calculated in step 4 into the expected value of pitch angle β r output by the pitch controller obtained in step 1 The mathematical expression and the mathematical expression of the actual pitch angle value β output by the variable pitch actuator are obtained, and the actual value β of the pitch angle of the pitch actuator is obtained and sent to the internal phase modulator of the wind turbine. At the same time, the variable The dead zone compensator of the pitch actuator calculates the pitch angle compensation value Δβ r by using the parameter values in Table 1 in step 3 and the expected pitch angle value β r obtained in step 4 according to the mathematical expression in step 2 and sends it to Internal phase modulators of wind turbines;
步骤6:根据接收到的桨距角值,即桨距角实际值β和桨距角补偿值Δβr,风轮机的内部调相器对风轮机进行调节。Step 6: According to the received pitch angle value, that is, the actual pitch angle value β and the pitch angle compensation value Δβ r , the internal phase regulator of the wind turbine adjusts the wind turbine.
附图5和附图6分别为风轮机转速和风力发电机输出的有功功率的仿真图。如图5所示,风轮机的转速基本保持在额定转速附近,变化平稳,没有出现大的波动。从附图6可以看出,有功功率输出已经基本恒定,基本维持在600KW,上下波动也是在极小范围内。这两个仿真图很好的说明了本发明的含有死区补偿器的风电机组的变桨距控制方法具有很好的效果,基本上抵消了变桨距执行机构死区所产生的偏差对风电机组的影响。Accompanying drawing 5 and accompanying drawing 6 are the simulation diagrams of the rotational speed of the wind turbine and the active power output by the wind generator, respectively. As shown in Figure 5, the speed of the wind turbine is basically maintained near the rated speed, and the change is stable without large fluctuations. It can be seen from Figure 6 that the active power output has been basically constant, basically maintained at 600KW, and the fluctuation is also within a very small range. These two simulation diagrams have well illustrated that the pitch control method of the wind turbine containing the dead zone compensator of the present invention has a good effect, and basically offsets the deviation caused by the dead zone of the pitch variable actuator on the wind power. impact on the unit.
虽然以上描述了本发明的具体实施方式,但是本领域内的熟练的技术人员应当理解,这些仅是举例说明,可以对这些实施方式做出多种变更或修改,而不背离本发明的原理和实质。本发明的范围仅由所附权利要求书限定。Although the specific embodiments of the present invention have been described above, those skilled in the art should understand that these are only examples, and various changes or modifications can be made to these embodiments without departing from the principles and principles of the present invention. substance. The scope of the invention is limited only by the appended claims.
Claims (4)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410432782.6A CN104343628B (en) | 2014-08-27 | 2014-08-27 | A kind of Wind turbines variable pitch control method containing dead-zone compensator |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410432782.6A CN104343628B (en) | 2014-08-27 | 2014-08-27 | A kind of Wind turbines variable pitch control method containing dead-zone compensator |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104343628A CN104343628A (en) | 2015-02-11 |
| CN104343628B true CN104343628B (en) | 2017-04-05 |
Family
ID=52499935
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410432782.6A Expired - Fee Related CN104343628B (en) | 2014-08-27 | 2014-08-27 | A kind of Wind turbines variable pitch control method containing dead-zone compensator |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104343628B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106150899B (en) * | 2015-04-21 | 2019-04-16 | 兰州交通大学 | A power optimization control method for front-end speed-regulated wind turbines |
| CN114000974B (en) * | 2021-10-31 | 2023-11-07 | 重庆科凯前卫风电设备有限责任公司 | Redundancy control method based on variable pitch variable speed wind generating set |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8353667B2 (en) * | 2008-08-27 | 2013-01-15 | General Electric Company | Method and apparatus for adjusting a yaw angle of a wind turbine |
| CN102182633B (en) * | 2011-06-07 | 2013-07-24 | 浙江运达风电股份有限公司 | Independent paddle change control method and device of large-scale wind generator system |
| EP2620639B1 (en) * | 2012-01-30 | 2016-01-27 | ALSTOM Renewable Technologies | A method for dampening oscillations in a wind turbine |
| CN102828910B (en) * | 2012-09-12 | 2014-09-10 | 华锐风电科技(集团)股份有限公司 | Wind generating set as well as paddle changing control method and device thereof |
| CN103016266B (en) * | 2013-01-11 | 2014-10-29 | 华北电力大学 | Wind generating set variable-pitch control method combining fuzzy feed-forward with linear active disturbance rejection |
-
2014
- 2014-08-27 CN CN201410432782.6A patent/CN104343628B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN104343628A (en) | 2015-02-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105179164B (en) | Wind-energy changing system sliding-mode control and device based on T-S fuzzy models | |
| CN106014877B (en) | Multi-fault diagnosis and fault-tolerant control of wind power generation system | |
| CN102136036B (en) | Double-feed wind power station equivalent modeling method applied to analysis on small signal stability of power system | |
| CN110566403A (en) | Wind power generation T-S fuzzy robust scheduling fault-tolerant control method | |
| CN108334672A (en) | Variable Speed Wind Power Generator maximal wind-energy capture method based on effective wind speed estimation | |
| CN105068424A (en) | Kaplan turbine adjusting system dynamic model suitable for electric power system analysis | |
| CN105545595A (en) | Wind turbine feedback linearization power control method based on radial basis function neural network | |
| WO2012163979A2 (en) | Assessment of power systems | |
| CN101807798B (en) | Section data integration method for power system safety and stability online analysis | |
| Wu et al. | Adaptive fault diagnosis and active tolerant control for wind energy conversion system | |
| CN103259289A (en) | Method for obtaining available transfer capability of wind power system with DFIG based on optimal power flow | |
| CN102879199A (en) | Full-condition wind turbine simulator | |
| CN115249980A (en) | New energy station impedance identification method and system driven by data and knowledge in combined mode | |
| CN104343628B (en) | A kind of Wind turbines variable pitch control method containing dead-zone compensator | |
| CN104993489B (en) | Determination Method of Optimal Load Shedding Amount in Isolated Grid Emergency | |
| CN110929405B (en) | An electro-pneumatic dynamic analysis method considering wind turbines and gas turbines | |
| Rao et al. | Wideband impedance online identification of wind farms based on combined data-driven and knowledge-driven | |
| CN109450310B (en) | A H∞ Robust Control Method for Suppressing Disturbances of Wind Turbines | |
| CN103207082B (en) | The test system of paddle change system of wind turbines and method of testing | |
| Benaaouinate et al. | Development of a useful wind turbine emulator based on permanent magnet DC motor | |
| CN106940959A (en) | The Megawatt fan analogy method observed based on acceleration | |
| CN106532785A (en) | Load frequency control method of considering novel cloud model | |
| Zhao et al. | Event-Triggered H-Infinity Pitch Control for Floating Offshore Wind Turbines | |
| CN105226635B (en) | Double-feedback type wind generator analogue system | |
| CN107086603A (en) | A Stochastic Fuzzy Continuous Power Flow Method for Power Systems Containing DFIG |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170405 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |