CN112149363A - Two-dimensional Jensen model and double-beam laser radar-based wake region fan power prediction method - Google Patents
Two-dimensional Jensen model and double-beam laser radar-based wake region fan power prediction method Download PDFInfo
- Publication number
- CN112149363A CN112149363A CN202010917961.4A CN202010917961A CN112149363A CN 112149363 A CN112149363 A CN 112149363A CN 202010917961 A CN202010917961 A CN 202010917961A CN 112149363 A CN112149363 A CN 112149363A
- Authority
- CN
- China
- Prior art keywords
- fan
- wind
- wake
- downstream
- downstream 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000011144 upstream manufacturing Methods 0.000 claims abstract description 109
- 238000005259 measurement Methods 0.000 claims description 77
- 238000004364 calculation method Methods 0.000 claims description 42
- 101150067085 los1 gene Proteins 0.000 description 12
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000004308 accommodation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
Description
技术领域technical field
本发明属于风机控制技术领域,具体涉及一种基于二维Jensen模型和双光束激光雷达的尾流区风机功率预测方法。The invention belongs to the technical field of fan control, and in particular relates to a method for predicting the power of a fan in a wake region based on a two-dimensional Jensen model and a double-beam laser radar.
背景技术Background technique
风能作为一种丰富、清洁的新能源,是世界上最有前途的可再生能源之一。随着电网输送和容纳能力的增加,风电场业主更加关注风机发电效率的提升。上游风机的尾流会使下游风机的功率输出出现显著降低,因此准确计算在上游风机尾流区中下游风机的功率输出值,对于风电场中风机的控制优化,提高风电场整体发电效率具有重要意义。目前,已有考虑尾流的风机功率计算方法。但是,受上下游风机偏航误差角度的影响,该方法对于下游风机所处尾流区的位置判断不准确,对于下游风机处于尾流区时计算的功率输出值并不准确。As a kind of abundant and clean new energy, wind energy is one of the most promising renewable energy sources in the world. With the increase of grid transmission and accommodation capacity, wind farm owners pay more attention to the improvement of wind turbine power generation efficiency. The wake of the upstream fan will significantly reduce the power output of the downstream fan. Therefore, the accurate calculation of the power output value of the downstream fan in the upstream fan wake area is important for the control optimization of the fans in the wind farm and improving the overall power generation efficiency of the wind farm. significance. At present, there are existing fan power calculation methods considering wake flow. However, due to the influence of the yaw error angle of the upstream and downstream fans, this method is inaccurate in judging the position of the wake region where the downstream fans are located, and the calculated power output value when the downstream fans are in the wake region is not accurate.
发明内容SUMMARY OF THE INVENTION
本发明目的在于提出一种基于二维Jensen模型和双光束激光雷达的尾流区风机功率预测方法。The purpose of the present invention is to propose a method for predicting the power of a fan in the wake region based on a two-dimensional Jensen model and a double-beam laser radar.
实现本发明目的的技术解决方案为:一种基于二维Jensen模型和双光束激光雷达的尾流区风机功率预测方法,具体步骤如下:The technical solution for realizing the purpose of the present invention is: a method for predicting the power of a blower fan in the wake region based on a two-dimensional Jensen model and a double-beam laser radar. The specific steps are as follows:
步骤1,采集上游风机的入流风速与入流风向,采集上游风机的偏航误差角度,采集下游风机双光束激光雷达左、右测风点测量的风速值,确定风场中风机之间的间距与方位角度;Step 1: Collect the inflow wind speed and inflow wind direction of the upstream fan, collect the yaw error angle of the upstream fan, collect the wind speed value measured by the left and right wind measuring points of the double-beam lidar of the downstream fan, and determine the distance between the fans in the wind farm. azimuth angle;
步骤2,根据二维Jensen模型及风机之间的间距与方位角度,计算下游风机雷达测风点所处截面位置的尾流半径,以及下游风机风轮平面所处截面位置的尾流半径;Step 2, according to the two-dimensional Jensen model and the spacing and azimuth angle between the fans, calculate the wake radius of the section where the radar wind measuring point of the downstream fan is located, and the wake radius of the section where the downstream fan rotor plane is located;
步骤3,判断下游风机雷达的两处测风点是均处于自然风速区,是一处测风点处于尾流区且另外一处测风点处于自然风速区,还是两处测风点均处于尾流区;Step 3: Determine whether the two wind measurement points of the downstream wind turbine radar are both in the natural wind speed area, whether one wind measurement point is in the wake area and the other wind measurement point is in the natural wind speed area, or both wind measurement points are in the natural wind speed area. wake area;
步骤4,根据上下游风机的连线与上游风机尾流中心轴线的夹角,判断下游风机处于尾流中心轴线的左侧还是右侧;Step 4: According to the angle between the connection line of the upstream and downstream fans and the central axis of the wake of the upstream fan, determine whether the downstream fan is on the left or the right of the central axis of the wake;
步骤5,根据雷达测风点处于风速区的位置,计算补偿后的下游风机雷达左、右测风点测量的风速值与补偿后的下游风机偏航误差角度;
步骤6,判断下游风机风轮平面是完全处于自然风速区,是部分处于尾流区,还是完全处于尾流区;Step 6, determine whether the plane of the downstream fan rotor is completely in the natural wind speed region, partially in the wake region, or completely in the wake region;
步骤7,根据下游风机风轮平面处于尾流区位置,计算下游风机的等效入流风速;Step 7: Calculate the equivalent inflow wind speed of the downstream fan according to the position of the wind wheel plane of the downstream fan in the wake region;
步骤8,计算下游风机的功率输出值。Step 8: Calculate the power output value of the downstream fan.
进一步的,步骤1中,利用风机的SCADA系统采集上游风机的入流风速与入流风向,利用安装在上游风机机舱上方的双光束激光雷达采集偏航误差角度,利用安装在下游风机机舱上方的双光束激光雷达采集左、右测风点测量的风速值,根据风场的建造选址确定风场中风机之间的间距与方位角度。Further, in
进一步的,步骤2中,根据二维Jensen模型及风机的分布距离与分布角度,分别计算下游风机雷达测风点所处截面位置的尾流半径、下游风机风轮平面所处截面位置的尾流半径,具体方法为:Further, in step 2, according to the two-dimensional Jensen model and the distribution distance and distribution angle of the fan, the wake radius at the cross-sectional position where the radar wind measurement point of the downstream fan is located, and the wake at the cross-sectional position where the downstream fan rotor plane is located are calculated respectively. Radius, the specific method is:
定义上下游风机连线与上游风机尾流中心轴线的夹角θL为:Define the angle θ L between the connection line of the upstream and downstream fans and the center axis of the wake of the upstream fan as:
θL=0.3CT·β1+θFWT-θx θ L = 0.3C T ·β 1 +θ FWT -θ x
式中,θx为入流风向角,θFWT为上游风机相对于下游风机所在的方位角,β1为上游风机的偏航误差角度,CT为风机升力系数;In the formula, θ x is the inflow wind direction angle, θ FWT is the azimuth angle of the upstream fan relative to the downstream fan, β 1 is the yaw error angle of the upstream fan, and C T is the fan lift coefficient;
定义下游风机雷达测风点所在截面与上游风机的垂直距离为Ll,计算公式为:Define the vertical distance between the section where the radar wind measurement point of the downstream fan is located and the upstream fan as L l , and the calculation formula is:
Ll=Ldcos(θL)-z0cos(α)L l =Ldcos(θ L )-z 0 cos(α)
式中,Ld为上下游风机的机舱连线距离,z0为下游风机雷达测风点与雷达的距离,α为下游风机雷达激光束与下游风机中轴线的夹角;In the formula, Ld is the distance between the nacelles of the upstream and downstream fans, z 0 is the distance between the radar wind measurement point of the downstream fan and the radar, and α is the angle between the radar laser beam of the downstream fan and the central axis of the downstream fan;
则下游风机雷达测风点所处截面位置的尾流半径计算公式为:Then the formula for calculating the wake radius at the cross-section position where the radar wind measuring point of the downstream fan is located is:
Rl=kLl+r0 R l =kL l +r 0
式中,k为尾流衰减系数,r0为风机的风轮半径。In the formula, k is the wake attenuation coefficient, and r 0 is the rotor radius of the fan.
定义下游风机风轮平面所在截面与上游风机的垂直距离为Lw,计算公式为:Define the vertical distance between the section where the rotor plane of the downstream fan is located and the upstream fan as L w , the calculation formula is:
Lw=Ldcos(θL)L w =Ldcos(θ L )
则下游风机风轮平面所处截面位置的尾流半径计算公式为:Then the formula for calculating the wake radius at the cross-sectional position where the wind wheel plane of the downstream fan is located is:
Rw=kLw+r0 R w =kL w +r 0
进一步的,步骤3中,判断下游风机雷达的两处测风点是均处于自然风速区,是一处测风点处于尾流区且另外一处测风点处于自然风速区,还是两处测风点均处于尾流区,具体方法为:Further, in step 3, it is determined whether the two wind measurement points of the downstream fan radar are both in the natural wind speed area, whether one wind measurement point is in the wake area and the other wind measurement point is in the natural wind speed area, or whether the two wind measurement points are in the natural wind speed area. The wind points are all in the wake region, and the specific methods are as follows:
若Ldsin(θL)-z0sin(α)>Rl,则下游风机雷达的两处测风点均处于自然风速区;If Ldsin(θ L )-z 0 sin(α)>R l , the two wind measurement points of the downstream wind turbine radar are in the natural wind speed area;
若Ldsin(θL)-z0sin(α)≤Rl&Ldsin(θL)+z0sin(α)>Rl,则下游风机雷达一处测风点处于尾流区,且另外一处测风点处于自然风速区;If Ldsin(θ L )-z 0 sin(α)≤R l &Ldsin(θ L )+z 0 sin(α)>R l , then one wind measurement point of the downstream wind turbine radar is in the wake region, and the other The wind measuring point is in the natural wind speed area;
若Ldsin(θL)-z0sin(α)≤Rl&Ldsin(θL)+z0sin(α)≤Rl,则下游风机雷达的两处测风点均处于尾流区;If Ldsin(θ L )-z 0 sin(α)≤R l &Ldsin(θ L )+z 0 sin(α)≤R l , then the two wind measurement points of the downstream wind turbine radar are both in the wake region;
其中,Ld为上下游风机的机舱连线距离,θL为上下游风机连线与上游风机尾流中心轴线的夹角,z0为下游风机雷达测风点与雷达的距离,α为下游风机雷达激光束与下游风机中轴线的夹角,Rl为下游风机雷达测风点所处截面位置的尾流半径。Among them, Ld is the distance between the nacelles of the upstream and downstream fans, θ L is the angle between the upstream and downstream fans and the center axis of the upstream fan wake, z 0 is the distance between the radar wind measurement point of the downstream fan and the radar, α is the downstream fan The angle between the radar laser beam and the central axis of the downstream fan, R l is the wake radius of the cross-section where the radar wind measurement point of the downstream fan is located.
进一步的,步骤4中,根据上下游风机的连线与上游风机尾流中心轴线的夹角,判断下游风机处于尾流中心轴线的左侧还是右侧,具体方法为:Further, in step 4, according to the angle between the connection line of the upstream and downstream fans and the central axis of the wake of the upstream fan, it is determined whether the downstream fan is on the left or the right of the central axis of the wake, and the specific method is as follows:
若上下游风机的连线与上游风机尾流中心轴线的夹角小于0,则下游风机处于尾流中心轴线左侧,否则处于尾流中心轴线右侧。If the angle between the connecting line of the upstream and downstream fans and the central axis of the wake of the upstream fan is less than 0, the downstream fan is on the left side of the central axis of the wake, otherwise it is on the right side of the central axis of the wake.
进一步的,步骤5中,根据雷达测风点处于风速区的位置,计算补偿后的下游风机雷达左、右测风点测量的风速值与补偿后的下游风机的偏航误差角度,具体方法为:Further, in
(1)若下游风机处于尾流中心轴线左侧,且雷达左测风点处于自然风速区、右测风点处于尾流区,则右测风点风速补偿系数为:(1) If the downstream fan is on the left side of the center axis of the wake, and the left wind measurement point of the radar is in the natural wind speed area and the right wind measurement point is in the wake area, the wind speed compensation coefficient of the right wind measurement point is:
式中,CT为风机升力系数,k为尾流衰减系数,r0为风机的风轮半径,Ll为下游风机雷达测风点所在截面与上游风机的垂直距离;where C T is the lift coefficient of the fan, k is the wake attenuation coefficient, r 0 is the radius of the fan rotor, and L l is the vertical distance between the section where the radar wind measurement point of the downstream fan is located and the upstream fan;
rrk为右测风点所处位置与尾流平面中心线的垂直距离,计算公式为:r rk is the vertical distance between the position of the right wind measurement point and the center line of the wake plane, and the calculation formula is:
rrk=Ldsin(θL)-z0sin(α)r rk =Ldsin(θ L )-z 0 sin(α)
式中,Ld为上下游风机的机舱连线距离,θL为上下游风机连线与上游风机尾流中心轴线的夹角,z0为下游风机雷达测风点与雷达的距离,α为下游风机雷达激光束与下游风机中轴线的夹角;In the formula, Ld is the distance between the nacelles of the upstream and downstream fans, θ L is the angle between the upstream and downstream fan lines and the center axis of the upstream fan wake, z 0 is the distance between the radar wind measurement point of the downstream fan and the radar, and α is the downstream The angle between the fan radar laser beam and the central axis of the downstream fan;
补偿后的下游风机雷达左、右测风点测量的风速值为:The wind speed values measured by the left and right wind measuring points of the downstream wind turbine radar after compensation are:
Va=Vlos1 V a =V los1
Vb=Vlos2/Crw V b =V los2 /C rw
式中,Vlos1,Vlos2分别为下游风机雷达左右测风点的实际测量风速值;In the formula, V los1 and V los2 are the actual measured wind speed values of the left and right wind measuring points of the downstream wind turbine radar respectively;
(2)若下游风机处于尾流中心轴线右侧,且雷达左测风点处于尾流区、右测风点处于自然风速区,则左测风点风速补偿系数为:(2) If the downstream fan is on the right side of the wake center axis, and the left wind measurement point of the radar is in the wake area and the right wind measurement point is in the natural wind speed area, the wind speed compensation coefficient of the left wind measurement point is:
rlk为左测风点所处位置与尾流平面中心线的垂直距离,计算公式为:r lk is the vertical distance between the position of the left wind measurement point and the center line of the wake plane, and the calculation formula is:
rlk=Ldsin(θL)-z0sin(α)r lk =Ldsin(θ L )-z 0 sin(α)
补偿后的下游风机雷达左、右测风点测量的风速值为:The wind speed values measured by the left and right wind measuring points of the downstream wind turbine radar after compensation are:
Va=Vlos1/Clw V a =V los1 /C lw
Vb=Vlos2 V b =V los2
(3)若下游风机处于尾流中心轴线左侧且雷达左右测风点均处于尾流区,则左右测风点风速补偿系数分别为:(3) If the downstream fan is on the left side of the center axis of the wake and the left and right wind measuring points of the radar are in the wake region, the wind speed compensation coefficients of the left and right wind measuring points are:
rlk与rrk分别为左、右测风点所处位置与尾流平面中心线的垂直距离,计算公式为:r lk and r rk are the vertical distances between the positions of the left and right wind measurement points and the center line of the wake plane, respectively. The calculation formula is:
rlk=Ldsin(θL)+z0sin(α)r lk =Ldsin(θ L )+z 0 sin(α)
rrk=|Ldsin(θL)-z0sin(α)|r rk =|Ldsin(θ L )-z 0 sin(α)|
补偿后的下游风机雷达左、右测风点测量的风速值为:The wind speed values measured by the left and right wind measuring points of the downstream wind turbine radar after compensation are:
Va=Vlos1/Clw V a =V los1 /C lw
Vb=Vlos2/Crw V b =V los2 /C rw
(4)若下游风机处于尾流中心轴线右侧且雷达左右测风点均处于尾流区,则左右测风点风速补偿系数分别为:(4) If the downstream fan is on the right side of the center axis of the wake and the left and right wind measuring points of the radar are in the wake region, the wind speed compensation coefficients of the left and right wind measuring points are:
rlk与rrk分别为左、右测风点所处位置与尾流平面中心线的垂直距离,计算公式为:r lk and r rk are the vertical distances between the positions of the left and right wind measurement points and the center line of the wake plane, respectively. The calculation formula is:
rlk=|Ldsin(θL)-z0sin(α)|r lk =|Ldsin(θ L )-z 0 sin(α)|
rrk=Ldsin(θL)+z0sin(α)r rk =Ldsin(θ L )+z 0 sin(α)
补偿后的下游风机雷达左、右测风点测量的风速值为:The wind speed values measured by the left and right wind measuring points of the downstream wind turbine radar after compensation are:
Va=Vlos1/Clw V a =V los1 /C lw
Vb=Vlos2/Crw V b =V los2 /C rw
(5)若下游风机雷达左右测风点均处于自然风速区,补偿后的下游风机雷达左、右测风点测量的风速值为:(5) If the left and right wind measurement points of the downstream fan radar are in the natural wind speed area, the wind speed measured by the left and right wind measurement points of the downstream fan radar after compensation is:
Va=Vlos1 V a =V los1
Vb=Vlos2 V b =V los2
根据下游风机所在风速区域的不同,计算补偿后的下游风机的偏航误差角度,具体公式为:According to the different wind speed areas of the downstream fans, the yaw error angle of the downstream fans after compensation is calculated. The specific formula is:
进一步的,步骤6中,判断下游风机风轮平面是完全处于自然风速区,是部分处于尾流区,还是完全处于尾流区,具体方法为:Further, in step 6, it is judged whether the plane of the wind wheel of the downstream fan is completely in the natural wind speed area, partly in the wake area, or completely in the wake area, and the specific method is as follows:
若Ldsin(θL)-r1>Rw,则下游风机风轮平面完全处于自然风速区;If Ldsin(θ L )-r 1 >R w , the plane of the downstream fan rotor is completely in the natural wind speed region;
若Ldsin(θL)-r1≤Rw&Ldsin(θL)+r1>Rw,则下游风机风轮平面部分处于尾流区;If Ldsin(θ L )-r 1 ≤R w &Ldsin(θ L )+r 1 >R w , the plane part of the downstream fan rotor is in the wake region;
若Ldsin(θL)-r1≤Rw&Ldsin(θL)+r1≤Rw,则下游风机风轮平面完全处于尾流区;If Ldsin(θ L )-r 1 ≤R w &Ldsin(θ L )+r 1 ≤R w , the plane of the downstream fan rotor is completely in the wake region;
其中,r1为风轮半径在垂直尾流中心线上的投影长度,计算公式为:Among them, r 1 is the projected length of the rotor radius on the vertical wake centerline, and the calculation formula is:
r1=r0cos(0.3CT·β1+β2)r 1 =r 0 cos( 0.3CT ·β 1 +β 2 )
其中,Ld为上下游风机的机舱连线距离,θL为上下游风机连线与上游风机尾流中心轴线的夹角,Rw为下游风机风轮平面所处截面位置的尾流半径,r0为风机风轮半径,β1为上游风机偏航误差角度,β2为补偿后的下游风机偏航误差角度,CT为风机升力系数。Among them, Ld is the distance between the nacelles of the upstream and downstream fans, θ L is the angle between the upstream and downstream fan lines and the central axis of the wake of the upstream fan, R w is the wake radius at the cross-sectional position of the plane of the downstream fan's rotor, r 0 is the radius of the fan rotor, β 1 is the yaw error angle of the upstream fan, β 2 is the yaw error angle of the downstream fan after compensation, and C T is the lift coefficient of the fan.
进一步的,根据下游风机风轮平面处于尾流区的位置,计算下游风机的等效入流风速,具体方法为:Further, according to the position of the downstream fan rotor plane in the wake region, calculate the equivalent inflow wind speed of the downstream fan, and the specific method is as follows:
(1)若下游风机风轮平面完全处于自然风速区,则下游风机等效入流风速为:(1) If the plane of the downstream fan rotor is completely in the natural wind speed zone, the equivalent inflow wind speed of the downstream fan is:
其中,u0为上游风机的入流风速。Among them, u 0 is the inflow wind speed of the upstream fan.
(2)若下游风机风轮平面部分处于尾流区,则下游风机风轮平面处于尾流区的面积为:(2) If the plane part of the downstream fan rotor is in the wake area, the area of the downstream fan rotor plane in the wake area is:
其中,Rw为下游风机风轮平面所处截面位置的尾流半径,r1是风轮半径在垂直尾流中心线上的投影长度;Among them, R w is the wake radius of the cross-sectional position of the downstream fan rotor plane, and r 1 is the projected length of the rotor radius on the vertical wake centerline;
令O1,O2为下游风机风轮平面处的尾流中心和下游风机风轮中心,B为下游风机风轮平面处的尾流区域边缘与下游风机风轮边缘的任意一个交点,θ1为连线O1O2与连线O1B的夹角,θ2为连线O1O2与连线O2B的夹角,计算公式为:Let O 1 , O 2 be the wake center at the plane of the downstream fan rotor and the center of the downstream fan rotor, B be any intersection point between the edge of the wake region at the plane of the downstream fan rotor and the edge of the downstream fan rotor, θ 1 is the included angle between the connecting line O 1 O 2 and the connecting line O 1 B, and θ 2 is the included angle between the connecting line O 1 O 2 and the connecting line O 2 B, and the calculation formula is:
其中,LR为下游风机中心距离上游风机尾流中心线的垂直距离:Among them, LR is the vertical distance from the center of the downstream fan to the centerline of the wake of the upstream fan:
LR=Ldsin(θL)L R =Ldsin(θ L )
其中,Ld为上下游风机的机舱连线距离,θL为上下游风机连线与上游风机尾流中心轴线的夹角;Among them, Ld is the distance between the nacelles of the upstream and downstream fans, and θ L is the angle between the upstream and downstream fan lines and the central axis of the upstream fan wake;
下游风机风轮面积为:The area of the downstream fan rotor is:
S=πr0 2 S=πr 0 2
其中,r0为风轮半径。Among them, r 0 is the radius of the rotor.
风轮在垂直尾流中心线上的投影面积为:The projected area of the wind rotor on the vertical wake centerline is:
S1=πr1 2 S 1 =πr 1 2
其中,r1为风轮半径在垂直尾流中心线上的投影长度;Among them, r 1 is the projected length of the rotor radius on the vertical wake centerline;
尾流区的风速ur计算公式为:The formula for calculating the wind speed ur in the wake region is:
式中,CT为风机升力系数,k为尾流衰减系数,r0为风机的风轮半径,Lw下游风机风轮平面所在截面与上游风机的垂直距离,r为风机风轮平面任意一点与尾流平面中心线的垂直距离,u0为上游风机的入流风速;In the formula, C T is the lift coefficient of the fan, k is the wake attenuation coefficient, r 0 is the rotor radius of the fan, L w is the vertical distance between the section where the plane of the fan rotor of the downstream fan is located and the upstream fan, and r is any point on the plane of the fan rotor The vertical distance from the centerline of the wake plane, u 0 is the inflow wind speed of the upstream fan;
下游风机等效入流风速计算公式为:The formula for calculating the equivalent inflow wind speed of the downstream fan is:
其中,a,b为积分上下限,取值为:Among them, a and b are the upper and lower limits of the integral, and the values are:
a=LR-r1 a=L R -r 1
b=Rw b= Rw
(3)若下游风机风轮平面部分处于尾流区,则下游风机等效入流风速计算公式为:(3) If the plane part of the rotor of the downstream fan is in the wake region, the calculation formula of the equivalent inflow wind speed of the downstream fan is:
其中,in,
a=LR-r1 a=L R -r 1
b=LR+r1 b=L R +r 1
进一步的,步骤8中,计算下游风机的功率输出值,具体公式为:Further, in step 8, the power output value of the downstream fan is calculated, and the specific formula is:
其中,ρ为空气密度,S为下游风机风轮面积,Cp为下游风机功率利用系数,为下游风机等效入流风速,β2为补偿后的下游风机偏航误差角度。Among them, ρ is the air density, S is the rotor area of the downstream fan, C p is the power utilization coefficient of the downstream fan, is the equivalent inflow wind speed of the downstream fan, and β 2 is the yaw error angle of the downstream fan after compensation.
一种基于二维Jensen模型和双光束激光雷达的尾流区风机功率预测系统,包括:A wake region fan power prediction system based on two-dimensional Jensen model and double-beam lidar, including:
数据采集模块,用于采集上游风机的入流风速与入流风向,采集上游风机的偏航误差角度,采集下游风机双光束激光雷达左、右测风点测量的风速值,确定风场中风机之间的间距与方位角度;The data acquisition module is used to collect the inflow wind speed and inflow wind direction of the upstream fan, collect the yaw error angle of the upstream fan, and collect the wind speed value measured by the left and right wind measuring points of the double-beam lidar of the downstream fan, and determine the distance between the fans in the wind farm. spacing and azimuth angle;
尾流半径计算模块,用于根据二维Jensen模型及风机之间的间距与方位角度,计算下游风机雷达测风点所处截面位置的尾流半径,以及下游风机风轮平面所处截面位置的尾流半径;The wake radius calculation module is used to calculate the wake radius at the cross-sectional position where the radar wind measurement point of the downstream fan is located, and the cross-sectional position of the downstream fan rotor plane according to the two-dimensional Jensen model and the distance and azimuth angle between the fans. wake radius;
测风点位置判断模块,用于判断下游风机雷达的两处测风点是均处于自然风速区,是一处测风点处于尾流区且另一处测风点处于自然风速区,还是两处测风点均处于尾流区;The wind measurement point position judgment module is used to judge whether the two wind measurement points of the downstream wind turbine radar are both in the natural wind speed area, whether one wind measurement point is in the wake area and the other wind measurement point is in the natural wind speed area, or whether the two wind measurement points are in the natural wind speed area. The wind measurement points are all in the wake region;
风机位置判断模块,用于根据上下游风机的连线与上游风机尾流中心轴线的夹角,判断下游风机处于尾流中心轴线的左侧还是右侧;The fan position judgment module is used to judge whether the downstream fan is on the left or right side of the wake center axis according to the angle between the connection line of the upstream and downstream fans and the center axis of the wake flow of the upstream fan;
补偿模块,用于根据雷达测风点处于风速区的位置,计算补偿后的下游风机雷达左、右测风点测量的风速值与补偿后的下游风机偏航误差角度;The compensation module is used to calculate the wind speed value measured by the left and right wind measuring points of the downstream fan after compensation and the yaw error angle of the downstream fan after compensation according to the position of the radar wind measuring point in the wind speed area;
风轮平面位置判断模块,用于判断下游风机风轮平面是完全处于自然风速区,是部分处于尾流区,还是完全处于尾流区;The wind rotor plane position judgment module is used to judge whether the wind rotor plane of the downstream fan is completely in the natural wind speed area, partly in the wake area, or completely in the wake area;
等效入流风速计算模块,用于根据下游风机风轮平面处于尾流区的位置,计算下游风机的等效入流风速;The equivalent inflow wind speed calculation module is used to calculate the equivalent inflow wind speed of the downstream fan according to the position of the downstream fan rotor plane in the wake region;
功率输出计算模块,用于计算下游风机的功率输出值。The power output calculation module is used to calculate the power output value of the downstream fan.
本发明与现有技术相比,其显著优点为:基于二维Jensen模型,结合上游风机偏航误差角度对于其尾流区位置的改变,准确计算了尾流区中补偿后的下游风机偏航误差角度,并以此为基础计算出在上游风机的尾流区中下游风机的等效入流风速,最终实现了下游风机尾流区功率输出值的准确计算,为风电场中尾流抑制的控制优化奠定基础。Compared with the prior art, the present invention has the significant advantages that: based on the two-dimensional Jensen model, combined with the change of the yaw error angle of the upstream fan to the position of the wake region, the yaw of the downstream fan after compensation in the wake region is accurately calculated Error angle, and based on this, the equivalent inflow wind speed of the downstream fans in the wake region of the upstream fan is calculated, and finally the accurate calculation of the power output value of the downstream fan wake region is realized, which is the control optimization of wake suppression in the wind farm. Lay the foundation.
附图说明Description of drawings
图1为本发明基于二维Jensen模型和双光束激光雷达的尾流区风机功率预测方法流程图。FIG. 1 is a flow chart of a method for predicting the power of a fan in the wake region based on a two-dimensional Jensen model and a dual-beam lidar according to the present invention.
图2为双光束激光雷达的结构示意图。Figure 2 is a schematic diagram of the structure of a dual-beam lidar.
图3为本发明判断下游风机风轮平面所在风速区的示意图。FIG. 3 is a schematic diagram of the present invention for judging the wind speed zone where the plane of the downstream fan rotor is located.
图4为本发明计算下游风机风轮平面处于尾流区面积的示意图。4 is a schematic diagram of calculating the area of the downstream fan rotor plane in the wake region according to the present invention.
图5为本发明补偿后的下游风机偏航误差角度图。FIG. 5 is a yaw error angle diagram of the downstream fan after compensation according to the present invention.
图6为本发明预测的下游风机尾流区功率输出值与实际功率输出值的比较图。FIG. 6 is a comparison diagram of the power output value in the wake region of the downstream fan predicted by the present invention and the actual power output value.
具体实施方式Detailed ways
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。In order to make the purpose, technical solutions and advantages of the present application more clearly understood, the present application will be described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application, but not to limit the present application.
如图1所示,基于二维Jensen模型和双光束激光雷达的尾流区风机功率预测方法,具体步骤如下:As shown in Figure 1, based on the two-dimensional Jensen model and dual-beam lidar fan power prediction method in the wake region, the specific steps are as follows:
步骤1,采集上游风机的入流风速与入流风向,采集上游风机的偏航误差角度,采集下游风机双光束激光雷达左、右测风点测量的风速值,确定风场中风机之间的间距与方位角度;Step 1: Collect the inflow wind speed and inflow wind direction of the upstream fan, collect the yaw error angle of the upstream fan, collect the wind speed value measured by the left and right wind measuring points of the double-beam lidar of the downstream fan, and determine the distance between the fans in the wind farm. azimuth angle;
利用风机的SCADA系统采集上游风机的入流风速与入流风向,利用安装在上游风机机舱上方的双光束激光雷达(如图2所示)采集偏航误差角度,利用安装在下游风机机舱上方的双光束激光雷达(如图2所示)采集左、右测风点测量的风速值,根据风场的建造选址确定风场中风机之间的间距与方位角度。The inflow wind speed and inflow wind direction of the upstream fan are collected by the SCADA system of the fan, the yaw error angle is collected by the double-beam lidar installed above the nacelle of the upstream fan (as shown in Figure 2), and the yaw error angle is collected by using the double beam installed above the nacelle of the downstream fan. The lidar (as shown in Figure 2) collects the wind speed values measured by the left and right wind measuring points, and determines the distance and azimuth angle between the fans in the wind farm according to the construction site selection of the wind farm.
步骤2,根据二维Jensen模型及风机的分布距离与分布角度,计算下游风机雷达测风点所处截面位置的尾流半径,计算下游风机风轮平面所处截面位置的尾流半径;Step 2, according to the two-dimensional Jensen model and the distribution distance and distribution angle of the fan, calculate the wake radius at the cross-sectional position where the radar wind measuring point of the downstream fan is located, and calculate the wake radius at the cross-sectional position where the downstream fan rotor plane is located;
定义上下游风机连线与上游风机尾流中心轴线的夹角θL为:Define the angle θ L between the connection line of the upstream and downstream fans and the center axis of the wake of the upstream fan as:
θL=0.3CT·β1+θFWT-θx θ L = 0.3C T ·β 1 +θ FWT -θ x
式中,θx为入流风向角,θFWT为上游风机相对于下游风机所在的方位角,β1为上游风机的偏航误差角度,CT为风机升力系数;In the formula, θ x is the inflow wind direction angle, θ FWT is the azimuth angle of the upstream fan relative to the downstream fan, β 1 is the yaw error angle of the upstream fan, and C T is the fan lift coefficient;
定义下游风机雷达测风点所在截面与上游风机的垂直距离为Ll,计算公式为:Define the vertical distance between the section where the radar wind measurement point of the downstream fan is located and the upstream fan as L l , and the calculation formula is:
Ll=Ldcos(θL)-z0cos(α)L l =Ldcos(θ L )-z 0 cos(α)
式中,Ld为上下游风机的机舱连线距离,z0为下游风机雷达测风点与雷达的距离,α为下游风机雷达激光束与下游风机中轴线的夹角;In the formula, Ld is the distance between the nacelles of the upstream and downstream fans, z 0 is the distance between the radar wind measurement point of the downstream fan and the radar, and α is the angle between the radar laser beam of the downstream fan and the central axis of the downstream fan;
则雷达测风点所处截面位置的尾流半径计算公式为:Then the formula for calculating the wake radius at the cross-section position of the radar wind measurement point is:
Rl=kLl+r0 R l =kL l +r 0
式中,k为尾流衰减系数,r0为风机的风轮半径。In the formula, k is the wake attenuation coefficient, and r 0 is the rotor radius of the fan.
定义下游风机风轮平面所在截面与上游风机的垂直距离为Lw,计算公式为:Define the vertical distance between the section where the rotor plane of the downstream fan is located and the upstream fan as L w , the calculation formula is:
Lw=Ldcos(θL)L w =Ldcos(θ L )
则下游风机风轮平面所处截面位置的尾流半径计算公式为:Then the formula for calculating the wake radius at the cross-sectional position where the wind wheel plane of the downstream fan is located is:
Rw=kLw+r0 R w =kL w +r 0
步骤3,判断下游风机雷达的两处测风点是均处于自然风速区,是一处测风点处于尾流区且另外一处测风点处于自然风速区,还是两处测风点均处于尾流区;Step 3: Determine whether the two wind measurement points of the downstream wind turbine radar are both in the natural wind speed area, whether one wind measurement point is in the wake area and the other wind measurement point is in the natural wind speed area, or both wind measurement points are in the natural wind speed area. wake area;
若Ldsin(θL)-z0sin(α)>Rl,则下游风机雷达的两处测风点均处于自然风速区;If Ldsin(θ L )-z 0 sin(α)>R l , the two wind measurement points of the downstream wind turbine radar are in the natural wind speed area;
若Ldsin(θL)-z0sin(α)≤Rl&Ldsin(θL)+z0sin(α)>Rl,则下游风机雷达一处测风点处于尾流区,且另外一处测风点处于自然风速区;If Ldsin(θ L )-z 0 sin(α)≤R l &Ldsin(θ L )+z 0 sin(α)>R l , then one wind measurement point of the downstream wind turbine radar is in the wake region, and the other The wind measuring point is in the natural wind speed area;
若Ldsin(θL)-z0sin(α)≤Rl&Ldsin(θL)+z0sin(α)≤Rl,则下游风机雷达的两处测风点均处于尾流区;If Ldsin(θ L )-z 0 sin(α)≤R l &Ldsin(θ L )+z 0 sin(α)≤R l , then the two wind measurement points of the downstream wind turbine radar are both in the wake region;
步骤4,根据上下游风机的连线与上游风机尾流中心轴线的夹角,判断下游风机处于尾流中心轴线左侧还是右侧;Step 4: According to the angle between the connection line of the upstream and downstream fans and the central axis of the wake of the upstream fan, determine whether the downstream fan is on the left or the right of the central axis of the wake;
若上下游风机的连线与上游风机尾流中心轴线的夹角小于0,则下游风机处于尾流中心轴线左侧,否则处于尾流中心轴线右侧。If the angle between the connecting line of the upstream and downstream fans and the central axis of the wake of the upstream fan is less than 0, the downstream fan is on the left side of the central axis of the wake, otherwise it is on the right side of the central axis of the wake.
步骤5,根据雷达测风点处于风速区的位置,计算补偿后的下游风机雷达左、右测风点测量的风速值,以及补偿后的下游风机偏航误差角度;
根据雷达测风点处于风速区的位置,计算补偿后的下游风机雷达左、右测风点测量的风速值分为五种具体情形:According to the position of the radar wind measuring point in the wind speed area, the wind speed values measured by the left and right wind measuring points of the downstream wind turbine after calculation and compensation are divided into five specific situations:
(1)若下游风机处于尾流中心轴线左侧,且雷达左测风点处于自然风速区、右测风点处于尾流区,则右测风点风速补偿系数为:(1) If the downstream fan is on the left side of the center axis of the wake, and the left wind measurement point of the radar is in the natural wind speed area and the right wind measurement point is in the wake area, the wind speed compensation coefficient of the right wind measurement point is:
式中,CT为风机升力系数;In the formula, C T is the lift coefficient of the fan;
rrk为右测风点所处位置与尾流平面中心线的垂直距离,计算公式为:r rk is the vertical distance between the position of the right wind measurement point and the center line of the wake plane, and the calculation formula is:
rrk=Ldsin(θL)-z0sin(α)r rk =Ldsin(θ L )-z 0 sin(α)
补偿后的下游风机雷达左、右测风点测量的风速值为:The wind speed values measured by the left and right wind measuring points of the downstream wind turbine radar after compensation are:
Va=Vlos1 V a =V los1
Vb=Vlos2/Crw V b =V los2 /C rw
式中,Vlos1,Vlos2分别为下游风机雷达左右测风点的实际测量风速值;In the formula, V los1 and V los2 are the actual measured wind speed values of the left and right wind measuring points of the downstream wind turbine radar respectively;
(2)若下游风机处于尾流中心轴线右侧,且雷达左测风点处于尾流区、右测风点处于自然风速区,则左测风点风速补偿系数为:(2) If the downstream fan is on the right side of the wake center axis, and the left wind measurement point of the radar is in the wake area and the right wind measurement point is in the natural wind speed area, the wind speed compensation coefficient of the left wind measurement point is:
式中,rlk为左测风点所处位置与尾流平面中心线的垂直距离,计算公式为:In the formula, r lk is the vertical distance between the position of the left wind measurement point and the center line of the wake plane, and the calculation formula is:
rlk=Ldsin(θL)-z0sin(α)r lk =Ldsin(θ L )-z 0 sin(α)
补偿后的下游风机雷达左、右测风点测量的风速值为:The wind speed values measured by the left and right wind measuring points of the downstream wind turbine radar after compensation are:
Va=Vlos1/Clw V a =V los1 /C lw
Vb=Vlos2 V b =V los2
(3)若下游风机处于尾流中心轴线左侧且雷达左右测风点均处于尾流区,则左右测风点风速补偿系数分别为:(3) If the downstream fan is on the left side of the center axis of the wake and the left and right wind measuring points of the radar are in the wake region, the wind speed compensation coefficients of the left and right wind measuring points are:
式中,rlk与rrk分别为左、右测风点所处位置与尾流平面中心线的垂直距离,计算公式为:where r lk and r rk are the vertical distances between the positions of the left and right wind measurement points and the centerline of the wake plane, respectively. The calculation formula is:
rlk=Ldsin(θL)+z0sin(α)r lk =Ldsin(θ L )+z 0 sin(α)
rrk=|Ldsin(θL)-z0sin(α)|r rk =|Ldsin(θ L )-z 0 sin(α)|
补偿后的下游风机雷达左、右测风点测量的风速值为:The wind speed values measured by the left and right wind measuring points of the downstream wind turbine radar after compensation are:
Va=Vlos1/Clw V a =V los1 /C lw
Vb=Vlos2/Crw V b =V los2 /C rw
(4)若下游风机处于尾流中心轴线右侧且雷达左右测风点均处于尾流区,则左右测风点风速补偿系数分别为:(4) If the downstream fan is on the right side of the center axis of the wake and the left and right wind measuring points of the radar are in the wake region, the wind speed compensation coefficients of the left and right wind measuring points are:
式中,rlk与rrk分别为左、右测风点所处位置与尾流平面中心线的垂直距离,计算公式为:where r lk and r rk are the vertical distances between the positions of the left and right wind measurement points and the centerline of the wake plane, respectively. The calculation formula is:
rlk=|Ldsin(θL)-z0sin(α)|r lk =|Ldsin(θ L )-z 0 sin(α)|
rrk=Ldsin(θL)+z0sin(α)r rk =Ldsin(θ L )+z 0 sin(α)
补偿后的下游风机雷达左、右测风点测量的风速值为:The wind speed values measured by the left and right wind measuring points of the downstream wind turbine radar after compensation are:
Va=Vlos1/Clw V a =V los1 /C lw
Vb=Vlos2/Crw V b =V los2 /C rw
(5)若下游风机雷达左右测风点均处于自然风速区,补偿后的下游风机雷达左、右测风点测量的风速值为:(5) If the left and right wind measurement points of the downstream fan radar are in the natural wind speed area, the wind speed measured by the left and right wind measurement points of the downstream fan radar after compensation is:
Va=Vlos1 V a =V los1
Vb=Vlos2 V b =V los2
根据下游风机所在风速区域的不同,计算补偿后的下游风机的偏航误差角度,具体公式为:According to the different wind speed areas of the downstream fans, the yaw error angle of the downstream fans after compensation is calculated. The specific formula is:
步骤6,判断下游风机风轮平面是完全处于自然风速区,是部分处于尾流区,还是完全处于尾流区,其示意图如图3所示;Step 6, determine whether the plane of the downstream fan rotor is completely in the natural wind speed region, partially in the wake region, or completely in the wake region, the schematic diagram is shown in Figure 3;
若Ldsin(θL)-r1>Rw,则下游风机风轮平面完全处于自然风速区;If Ldsin(θ L )-r 1 >R w , the plane of the downstream fan rotor is completely in the natural wind speed region;
若Ldsin(θL)-r1≤Rw&Ldsin(θL)+r1>Rw,则下游风机风轮平面部分处于尾流区;If Ldsin(θ L )-r 1 ≤R w &Ldsin(θ L )+r 1 >R w , the plane part of the downstream fan rotor is in the wake region;
若Ldsin(θL)-r1≤Rw&Ldsin(θL)+r1≤Rw,则下游风机风轮平面完全处于尾流区;If Ldsin(θ L )-r 1 ≤R w &Ldsin(θ L )+r 1 ≤R w , the plane of the downstream fan rotor is completely in the wake region;
其中,r1为风轮半径在垂直尾流中心线上的投影长度,计算公式为:Among them, r 1 is the projected length of the rotor radius on the vertical wake centerline, and the calculation formula is:
r1=r0cos(0.3CT·β1+β2)r 1 =r 0 cos( 0.3CT ·β 1 +β 2 )
步骤7,根据下游风机风轮平面处于尾流区的不同,计算下游风机的等效入流风速;Step 7: Calculate the equivalent inflow wind speed of the downstream fan according to the difference between the plane of the downstream fan's rotor in the wake region;
根据下游风机风轮平面处于尾流区的位置,分为三种具体情形:According to the position of the downstream fan rotor plane in the wake region, there are three specific situations:
(1)若下游风机风轮平面完全处于自然风速区,则下游风机等效入流风速为:(1) If the plane of the downstream fan rotor is completely in the natural wind speed zone, the equivalent inflow wind speed of the downstream fan is:
其中,u0为上游风机的入流风速。Among them, u 0 is the inflow wind speed of the upstream fan.
(2)若下游风机风轮平面部分处于尾流区,则下游风机风轮平面处于尾流区的面积为:(2) If the plane part of the downstream fan rotor is in the wake area, the area of the downstream fan rotor plane in the wake area is:
如图4所示,令O1,O2为下游风机风轮平面处的尾流中心和下游风机风轮中心,B为下游风机风轮平面处的尾流区域边缘与下游风机风轮边缘的任意一个交点,θ1为连线O1O2与连线O1B的夹角,θ2为连线O1O2与连线O2B的夹角,计算公式为:As shown in Figure 4, let O 1 and O 2 be the wake center at the plane of the downstream fan rotor and the center of the downstream fan rotor, and B be the distance between the edge of the wake region at the plane of the downstream fan rotor and the edge of the downstream fan rotor For any intersection, θ 1 is the angle between the line O 1 O 2 and the line O 1 B, and θ 2 is the angle between the line O 1 O 2 and the line O 2 B. The calculation formula is:
其中,LR为下游风机中心距离上游风机尾流中心线的垂直距离:Among them, LR is the vertical distance from the center of the downstream fan to the centerline of the wake of the upstream fan:
LR=Ldsin(θL)L R =Ldsin(θ L )
下游风机风轮面积为:The area of the downstream fan rotor is:
S=πr0 2 S=πr 0 2
风轮在垂直尾流中心线上的投影面积为:The projected area of the wind rotor on the vertical wake centerline is:
S1=πr1 2 S 1 =πr 1 2
尾流区的风速ur计算公式为:The formula for calculating the wind speed ur in the wake region is:
下游风机等效入流风速计算公式为:The formula for calculating the equivalent inflow wind speed of the downstream fan is:
其中,a,b为积分上下限,取值为:Among them, a and b are the upper and lower limits of the integral, and the values are:
a=LR-r1 a=L R -r 1
b=Rw b= Rw
(3)若下游风机风轮平面部分处于尾流区,则下游风机等效入流风速计算公式为:(3) If the plane part of the rotor of the downstream fan is in the wake region, the calculation formula of the equivalent inflow wind speed of the downstream fan is:
其中,in,
a=LR-r1 a=L R -r 1
b=LR+r1 b=L R +r 1
步骤8,计算下游风机的功率输出值;Step 8, calculate the power output value of the downstream fan;
计算公式为:The calculation formula is:
本发明还提出一种基于二维Jensen模型和双光束激光雷达的尾流区风机功率预测系统,包括:The present invention also proposes a wake region fan power prediction system based on the two-dimensional Jensen model and the double-beam laser radar, including:
数据采集模块,用于采集上游风机的入流风速与入流风向,采集上游风机的偏航误差角度,采集下游风机双光束激光雷达左、右测风点测量的风速值,确定风场中风机之间的间距与方位角度;The data acquisition module is used to collect the inflow wind speed and inflow wind direction of the upstream fan, collect the yaw error angle of the upstream fan, and collect the wind speed value measured by the left and right wind measuring points of the double-beam lidar of the downstream fan, and determine the distance between the fans in the wind farm. spacing and azimuth angle;
尾流半径计算模块,用于根据二维Jensen模型及风机之间的间距与方位角度,计算下游风机雷达测风点所处截面位置的尾流半径,以及下游风机风轮平面所处截面位置的尾流半径;The wake radius calculation module is used to calculate the wake radius at the cross-sectional position where the radar wind measurement point of the downstream fan is located, and the cross-sectional position of the downstream fan rotor plane according to the two-dimensional Jensen model and the distance and azimuth angle between the fans. wake radius;
测风点位置判断模块,用于判断下游风机雷达的两处测风点是均处于自然风速区,是一处测风点处于尾流区且另一处测风点处于自然风速区,还是两处测风点均处于尾流区;The wind measurement point position judgment module is used to judge whether the two wind measurement points of the downstream wind turbine radar are both in the natural wind speed area, whether one wind measurement point is in the wake area and the other wind measurement point is in the natural wind speed area, or whether the two wind measurement points are in the natural wind speed area. The wind measurement points are all in the wake region;
风机位置判断模块,用于根据上下游风机的连线与上游风机尾流中心轴线的夹角,判断下游风机处于尾流中心轴线的左侧还是右侧;The fan position judgment module is used to judge whether the downstream fan is on the left or right side of the wake center axis according to the angle between the connection line of the upstream and downstream fans and the center axis of the wake flow of the upstream fan;
补偿模块,用于根据雷达测风点处于风速区的位置,计算补偿后的下游风机雷达左、右测风点测量的风速值与补偿后的下游风机偏航误差角度;The compensation module is used to calculate the wind speed value measured by the left and right wind measuring points of the downstream fan after compensation and the yaw error angle of the downstream fan after compensation according to the position of the radar wind measuring point in the wind speed area;
风轮平面位置判断模块,用于判断下游风机风轮平面是完全处于自然风速区,是部分处于尾流区,还是完全处于尾流区;The wind rotor plane position judgment module is used to judge whether the wind rotor plane of the downstream fan is completely in the natural wind speed area, partly in the wake area, or completely in the wake area;
等效入流风速计算模块,用于根据下游风机风轮平面处于尾流区的位置,计算下游风机的等效入流风速;The equivalent inflow wind speed calculation module is used to calculate the equivalent inflow wind speed of the downstream fan according to the position of the downstream fan rotor plane in the wake region;
功率输出计算模块,用于计算下游风机的功率输出值。The power output calculation module is used to calculate the power output value of the downstream fan.
一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现基于二维Jensen模型和双光束激光雷达的尾流区风机功率预测。A computer device, comprising a memory, a processor, and a computer program stored on the memory and running on the processor, the processor implements a wake based on a two-dimensional Jensen model and a dual-beam lidar when the processor executes the computer program District fan power prediction.
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现基于二维Jensen模型和双光束激光雷达的尾流区风机功率预测。A computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, realizes power prediction of a wake region fan based on a two-dimensional Jensen model and a double-beam lidar.
实施例Example
为了验证本发明方案的有效性,在某风电场编号为A04的风机上进行如下仿真实验。风机A04附近存在两台风机A03与A05,风机A03相对于风机A04的方位角为35°,机舱连线距离为316m,风机A05相对于风机A04的方位角为218°,机舱连线距离为314m;三台风机功率利用系数均为0.55,升力系数均为0.5,风轮半径均为56.5m,雷达测风点与雷达的距离为80m,雷达发出的激光与风机中轴线的夹角为30°,尾流衰减系数为0.075。In order to verify the effectiveness of the scheme of the present invention, the following simulation experiments are carried out on a wind turbine numbered A04 in a wind farm. There are two fans A03 and A05 near the fan A04. The azimuth angle of the fan A03 relative to the fan A04 is 35°, the connecting distance of the engine room is 316m, the azimuth angle of the fan A05 relative to the fan A04 is 218°, and the connecting distance of the engine room is 314m. The power utilization coefficients of the three fans are all 0.55, the lift coefficients are all 0.5, the radius of the wind rotor is 56.5m, the distance between the radar wind measurement point and the radar is 80m, and the angle between the laser emitted by the radar and the central axis of the fan is 30° , the wake attenuation coefficient is 0.075.
以风机A04作为下游风机,补偿后的下游风机A04的偏航误差角度如图5所示,可以看到修正后的风机A04的偏航误差角成功消除了尾流效应对于偏航误差计算的干扰,使偏航误差值在尾流区的畸变程度明显降低。以修正后的下游风机A04偏航误差角度为基础,计算出在风机A04风轮平面处于尾流区时对应的风速入流角。在空间360°中,当风机A03为风机A04的上游风机时,在风机A03产生的尾流区中,风机A04风轮平面处的尾流半径为68.8m,风机A04风轮平面完全处于尾流区的入流风向角度为29°-44°,风机A04风轮平面部分处于尾流区的入流风向角度为12°-28°,45°-62°,其余入流风向角度下的风机A04风轮平面完全处于自然风速区。在空间360°中,当风机A05为风机A04的上游风机时,在风机A05产生的尾流区中,风机A04风轮平面处的尾流半径为68.8m,风机A04风轮平面完全处于尾流区的入流风向角度为212°-225°,风机A04风轮平面部分处于尾流区的入流风向角度为193°-211°,226°-242°,其余入流风向角度下的风机A04风轮平面完全处于自然风速区。Taking the fan A04 as the downstream fan, the yaw error angle of the downstream fan A04 after compensation is shown in Figure 5. It can be seen that the corrected yaw error angle of the fan A04 has successfully eliminated the interference of the wake effect on the yaw error calculation. , so that the degree of distortion of the yaw error value in the wake region is significantly reduced. Based on the corrected yaw error angle of the downstream fan A04, the corresponding wind speed inflow angle is calculated when the plane of the rotor A04 of the fan is in the wake region. In the space of 360°, when the fan A03 is the upstream fan of the fan A04, in the wake area generated by the fan A03, the radius of the wake at the plane of the fan A04 is 68.8m, and the plane of the fan A04 is completely in the wake. The inflow wind direction angle of the area is 29°-44°, the inflow wind direction angle of the fan A04 in the wake area is 12°-28°, 45°-62°, and the rest of the inflow wind direction angle is the fan A04 wind wheel plane Completely in the natural wind speed zone. In 360° of space, when fan A05 is the upstream fan of fan A04, in the wake area generated by fan A05, the radius of the wake at the plane of fan A04 is 68.8m, and the plane of fan A04 is completely in the wake. The inflow wind direction angle of the area is 212°-225°, the inflow wind direction angle of the fan A04 in the wake area is 193°-211°, 226°-242°, and the rest of the inflow wind direction angle of the fan A04 wind wheel plane Completely in the natural wind speed zone.
根据计算所得的风机A04风轮平面处于尾流区时对应的风速入流角,计算风机A04的等效入流风速,并最终计算得到风机A04尾流区功率输出值,如图6所示。将计算所得的尾流区功率输出值与实际功率进行比较,得到预测功率与实际功率的综合误差为4.26%,可以看出本发明修正方法简单有效,能够准确计算风机尾流区功率输出值。According to the calculated inflow angle of the wind speed when the plane of the wind wheel of the fan A04 is in the wake region, the equivalent inflow wind speed of the fan A04 is calculated, and the power output value of the wake region of the fan A04 is finally calculated, as shown in Figure 6. Comparing the calculated power output value in the wake region with the actual power, the comprehensive error between the predicted power and the actual power is 4.26%. It can be seen that the correction method of the present invention is simple and effective, and can accurately calculate the power output value in the wake region of the fan.
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be implemented by instructing relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage In the medium, when the computer program is executed, it may include the processes of the above-mentioned method embodiments. Wherein, any reference to memory, storage, database or other medium used in the various embodiments provided in this application may include non-volatile and/or volatile memory. Nonvolatile memory may include read only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous chain Road (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。The technical features of the above embodiments can be combined arbitrarily. In order to make the description simple, all possible combinations of the technical features in the above embodiments are not described. However, as long as there is no contradiction in the combination of these technical features It is considered to be the range described in this specification.
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。The above-mentioned embodiments only represent several embodiments of the present application, and the descriptions thereof are specific and detailed, but should not be construed as a limitation on the scope of the invention patent. It should be pointed out that for those skilled in the art, without departing from the concept of the present application, several modifications and improvements can be made, which all belong to the protection scope of the present application. Therefore, the scope of protection of the patent of the present application shall be subject to the appended claims.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010917961.4A CN112149363B (en) | 2020-09-03 | 2020-09-03 | Two-dimensional Jensen model and double-beam laser radar-based wake region fan power prediction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010917961.4A CN112149363B (en) | 2020-09-03 | 2020-09-03 | Two-dimensional Jensen model and double-beam laser radar-based wake region fan power prediction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112149363A true CN112149363A (en) | 2020-12-29 |
CN112149363B CN112149363B (en) | 2021-09-03 |
Family
ID=73889171
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010917961.4A Active CN112149363B (en) | 2020-09-03 | 2020-09-03 | Two-dimensional Jensen model and double-beam laser radar-based wake region fan power prediction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112149363B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112879220A (en) * | 2021-03-16 | 2021-06-01 | 上海电气风电集团股份有限公司 | Fan control method and system and readable storage medium |
CN113033009A (en) * | 2021-03-31 | 2021-06-25 | 西安热工研究院有限公司 | Real-time calculation method for wake flow loss of offshore wind farm in service |
CN114511158A (en) * | 2022-04-19 | 2022-05-17 | 南京理工大学 | Wind turbine power prediction method based on wake deflection effect and 2DJensen model |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020462A (en) * | 2012-12-21 | 2013-04-03 | 华北电力大学 | Wind power plant probability output power calculation method considering complex wake effect model |
CN106919731A (en) * | 2015-12-25 | 2017-07-04 | 中国电力科学研究院 | A kind of Wind turbines wake flow for different wind angles determines method |
CN108953060A (en) * | 2018-03-30 | 2018-12-07 | 浙江大学 | Wind power plant field grade Yaw control method based on laser radar anemometer |
CN109255184A (en) * | 2018-09-11 | 2019-01-22 | 河海大学 | A kind of determination method and system of the full tail flow field wind speed profile of wind energy conversion system |
US20190370418A1 (en) * | 2018-05-29 | 2019-12-05 | Beijing Goldwind Science & Creation Windpower Equipment Co., Ltd. | Method and apparatus for arranging wind turbines based on rapid accessment fluid model and wake model |
-
2020
- 2020-09-03 CN CN202010917961.4A patent/CN112149363B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020462A (en) * | 2012-12-21 | 2013-04-03 | 华北电力大学 | Wind power plant probability output power calculation method considering complex wake effect model |
CN106919731A (en) * | 2015-12-25 | 2017-07-04 | 中国电力科学研究院 | A kind of Wind turbines wake flow for different wind angles determines method |
CN108953060A (en) * | 2018-03-30 | 2018-12-07 | 浙江大学 | Wind power plant field grade Yaw control method based on laser radar anemometer |
US20190370418A1 (en) * | 2018-05-29 | 2019-12-05 | Beijing Goldwind Science & Creation Windpower Equipment Co., Ltd. | Method and apparatus for arranging wind turbines based on rapid accessment fluid model and wake model |
CN109255184A (en) * | 2018-09-11 | 2019-01-22 | 河海大学 | A kind of determination method and system of the full tail flow field wind speed profile of wind energy conversion system |
Non-Patent Citations (1)
Title |
---|
袁明友等: "基于二维尾流模型的风电场产能分析", 《可再生能源》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112879220A (en) * | 2021-03-16 | 2021-06-01 | 上海电气风电集团股份有限公司 | Fan control method and system and readable storage medium |
CN112879220B (en) * | 2021-03-16 | 2022-11-01 | 上海电气风电集团股份有限公司 | Fan control method and system and readable storage medium |
CN113033009A (en) * | 2021-03-31 | 2021-06-25 | 西安热工研究院有限公司 | Real-time calculation method for wake flow loss of offshore wind farm in service |
CN114511158A (en) * | 2022-04-19 | 2022-05-17 | 南京理工大学 | Wind turbine power prediction method based on wake deflection effect and 2DJensen model |
Also Published As
Publication number | Publication date |
---|---|
CN112149363B (en) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112149363B (en) | Two-dimensional Jensen model and double-beam laser radar-based wake region fan power prediction method | |
CN108953060B (en) | Field-level yaw control method for wind farm based on lidar anemometer | |
Hansen et al. | Wind farm noise: measurement, assessment, and control | |
CN103823979B (en) | A kind of wind power plant noise prediction method | |
CN105160060B (en) | A method for determining the theoretical power of wind farms based on actual power curve fitting | |
CN106815456A (en) | A kind of Wind turbines power characteristic evaluation method | |
CN112096576B (en) | Method for improving annual generated energy of multiple fan arrays based on wake field optimization control | |
US20240309843A1 (en) | Cooperative operation optimization control method for wind turbine groups | |
CN114511158B (en) | Wind turbine power prediction method based on wake deflection effect and 2DJensen model | |
CN113219478B (en) | Method, device and storage medium for measuring and predicting wind speed and wind direction in fan impeller range | |
CN111894802B (en) | Yaw error correction method in wind turbine wake region based on Jensen model and dual beam radar | |
CN110307121B (en) | A method for optimizing the blade angle of wind turbines | |
CN104794259A (en) | Wind power plant on-grid energy deviation computing method based on anemometer tower mutual verification | |
CN114294156A (en) | Yaw feedback control method and system for full-wind-condition gain scheduling of wind turbine generator | |
CN116050287A (en) | Modeling method and device for wake analysis of offshore floating wind turbines | |
CN108269197B (en) | Wind turbine power characteristic evaluation method and device | |
CN106014878B (en) | The test method and system of wind driven generator unit yaw system action error | |
CN115203622A (en) | Laser radar scanning path planning method for measuring wind power plant wake flow | |
CN114294157A (en) | Wind power fan yaw deviation identification method, deviation correction method, electronic device and medium | |
CN111120202A (en) | Yaw angle adjustment method, device, medium and equipment for wind turbine | |
CN119005420A (en) | Automatic arrangement system of fan photovoltaic equipment based on satellite map vision area avoidance | |
CN109766611B (en) | A Simplified Prediction Method of Wind Farm Power Considering Terrain Slope | |
CN108734601A (en) | A kind of wind power plant generated energy evaluation method and system | |
CN116720437A (en) | Modeling method, system, terminal and medium for clearance distance of wind turbine generator blade | |
Dai et al. | Study on obtaining real power curve of wind turbines using SCADA data |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20210521 Address after: 102600 block D, AVIC Plaza, No.15, Ronghua Road, economic and Technological Development Zone, Daxing District, Beijing Applicant after: Beijing Huaxin science and Technology Co.,Ltd. Address before: 210094 Xuanwu District, Jiangsu, Xiaolingwei 200, Nanjing Applicant before: NANJING University OF SCIENCE AND TECHNOLOGY |
|
TA01 | Transfer of patent application right | ||
CI02 | Correction of invention patent application |
Correction item: Applicant|Address Correct: NANJING University OF SCIENCE AND TECHNOLOGY|210094 No. 200, Xiaolingwei, Xuanwu District, Nanjing City, Jiangsu Province False: Beijing Huaxin science and Technology Co.,Ltd.|102600 Building D, AVIC Plaza, No. 15 Ronghua Road, Economic and Technological Development Zone, Daxing District, Beijing Number: 23-02 Volume: 37 |
|
CI02 | Correction of invention patent application | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20201229 Assignee: Beijing Huaxin science and Technology Co.,Ltd. Assignor: NANJING University OF SCIENCE AND TECHNOLOGY Contract record no.: X2024980001507 Denomination of invention: A power prediction method for wind turbines in wake region based on two-dimensional Jensen model and dual beam LiDAR Granted publication date: 20210903 License type: Exclusive License Record date: 20240129 |
|
EE01 | Entry into force of recordation of patent licensing contract |