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 PDF

Info

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
Application number
CN202010917961.4A
Other languages
Chinese (zh)
Other versions
CN112149363B (en
Inventor
邱颖宁
李绪阳
冯延晖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Huaxin Science And Technology Co ltd
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN202010917961.4A priority Critical patent/CN112149363B/en
Publication of CN112149363A publication Critical patent/CN112149363A/en
Application granted granted Critical
Publication of CN112149363B publication Critical patent/CN112149363B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids

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

The invention discloses a wake flow area fan power prediction method based on a two-dimensional Jensen model and a double-beam laser radar, which comprises the steps of calculating the wake flow radius of the cross section position of a radar wind measuring point of a downstream fan and the wake flow radius of the cross section position of a wind wheel plane of the downstream fan; judging whether the downstream fan is positioned on the left side or the right side of the wake central axis according to the included angle between the connecting line of the upstream fan and the downstream fan and the wake central axis of the upstream fan; calculating wind speed values measured by the compensated radar left and right wind measuring points of the downstream fan and a yaw error angle of the compensated downstream fan according to the position of the radar wind measuring point in the wind speed area; calculating 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 flow area; and calculating the power output value of the downstream fan. The invention also discloses a corresponding prediction system. The method can accurately calculate the power output value of the downstream fan with the yaw error in the wake flow area.

Description

基于二维Jensen模型和双光束激光雷达的尾流区风机功率预 测方法Power prediction of wind turbines in wake region based on two-dimensional Jensen model and dual-beam lidar test method

技术领域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,根据雷达测风点处于风速区的位置,计算补偿后的下游风机雷达左、右测风点测量的风速值与补偿后的下游风机偏航误差角度;Step 5, according to the position of the radar wind measurement point in the wind speed area, calculate the wind speed value measured by the left and right wind measurement points of the downstream wind turbine after compensation and the yaw error angle of the downstream wind turbine after compensation;

步骤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 step 1, 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, and the yaw error angle is collected by using the double beam installed above the nacelle of the downstream fan. The lidar collects the wind speed values measured by the left and right wind measuring points, and determines the spacing and azimuth angle between the fans in the wind farm according to the construction site selection of the wind farm.

进一步的,步骤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·β1FWTx θ L = 0.3C T ·β 1FWTx

式中,θ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 step 5, according to the position of the radar wind measuring point in the wind speed area, 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, and the specific method is: :

(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:

Figure BDA0002665701110000031
Figure BDA0002665701110000031

式中,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:

Figure BDA0002665701110000041
Figure BDA0002665701110000041

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:

Figure BDA0002665701110000042
Figure BDA0002665701110000042

Figure BDA0002665701110000043
Figure BDA0002665701110000043

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:

Figure BDA0002665701110000051
Figure BDA0002665701110000051

Figure BDA0002665701110000052
Figure BDA0002665701110000052

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:

Figure BDA0002665701110000053
Figure BDA0002665701110000053

进一步的,步骤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·β12)r 1 =r 0 cos( 0.3CT ·β 12 )

其中,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:

Figure BDA0002665701110000061
Figure BDA0002665701110000061

其中,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:

Figure BDA0002665701110000062
Figure BDA0002665701110000062

其中,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:

Figure BDA0002665701110000063
Figure BDA0002665701110000063

Figure BDA0002665701110000064
Figure BDA0002665701110000064

其中,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:

Figure BDA0002665701110000071
Figure BDA0002665701110000071

式中,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:

Figure BDA0002665701110000072
Figure BDA0002665701110000072

其中,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:

Figure BDA0002665701110000073
Figure BDA0002665701110000073

其中,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:

Figure BDA0002665701110000074
Figure BDA0002665701110000074

其中,ρ为空气密度,S为下游风机风轮面积,Cp为下游风机功率利用系数,

Figure BDA0002665701110000075
为下游风机等效入流风速,β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,
Figure BDA0002665701110000075
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·β1FWTx θ L = 0.3C T ·β 1FWTx

式中,θ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,根据雷达测风点处于风速区的位置,计算补偿后的下游风机雷达左、右测风点测量的风速值,以及补偿后的下游风机偏航误差角度;Step 5, according to the position of the radar wind measurement point in the wind speed area, calculate the wind speed value measured by the left and right wind measurement points of the downstream wind turbine radar after compensation, and the yaw error angle of the downstream wind turbine after compensation;

根据雷达测风点处于风速区的位置,计算补偿后的下游风机雷达左、右测风点测量的风速值分为五种具体情形: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:

Figure BDA0002665701110000101
Figure BDA0002665701110000101

式中,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:

Figure BDA0002665701110000111
Figure BDA0002665701110000111

式中,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:

Figure BDA0002665701110000112
Figure BDA0002665701110000112

Figure BDA0002665701110000113
Figure BDA0002665701110000113

式中,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:

Figure BDA0002665701110000114
Figure BDA0002665701110000114

Figure BDA0002665701110000115
Figure BDA0002665701110000115

式中,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:

Figure BDA0002665701110000121
Figure BDA0002665701110000121

步骤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·β12)r 1 =r 0 cos( 0.3CT ·β 12 )

步骤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:

Figure BDA0002665701110000123
Figure BDA0002665701110000123

其中,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:

Figure BDA0002665701110000122
Figure BDA0002665701110000122

如图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:

Figure BDA0002665701110000131
Figure BDA0002665701110000131

Figure BDA0002665701110000132
Figure BDA0002665701110000132

其中,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:

Figure BDA0002665701110000133
Figure BDA0002665701110000133

下游风机等效入流风速计算公式为:The formula for calculating the equivalent inflow wind speed of the downstream fan is:

Figure BDA0002665701110000134
Figure BDA0002665701110000134

其中,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:

Figure BDA0002665701110000135
Figure BDA0002665701110000135

其中,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:

Figure BDA0002665701110000141
Figure BDA0002665701110000141

本发明还提出一种基于二维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)

1. The method for predicting the power of the fan in the wake zone based on the two-dimensional Jensen model and the double-beam laser radar is characterized by comprising the following specific steps of:
step 1, collecting inflow wind speed and inflow wind direction of an upstream fan, collecting yaw error angle of the upstream fan, collecting wind speed values measured by a left wind measurement point and a right wind measurement point of a double-beam laser radar of a downstream fan, and determining the distance and azimuth angle between the fans in a wind field;
step 2, calculating the wake flow radius of the cross section position of the radar wind measuring point of the downstream fan and the wake flow radius of the cross section position of the wind wheel plane of the downstream fan according to the two-dimensional Jensen model and the space and the azimuth angle between the fans;
step 3, judging whether two wind measuring points of the downstream fan radar are both in the natural wind speed area, one wind measuring point is in the wake area and the other wind measuring point is in the natural wind speed area, or the two wind measuring points are both in the wake area;
step 4, judging whether the downstream fan is positioned on the left side or the right side of the wake central axis according to the included angle between the connecting line of the upstream fan and the downstream fan and the wake central axis of the upstream fan;
step 5, calculating wind speed values measured by the compensated radar left and right wind measuring points of the downstream fan and a compensated yaw error angle of the downstream fan according to the position of the radar wind measuring point in the wind speed area;
step 6, judging whether the wind wheel plane of the downstream fan is completely in the natural wind speed area, partially in the wake area or completely in the wake area;
step 7, calculating 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 area;
and 8, calculating the power output value of the downstream fan.
2. The method for predicting the fan power in the wake zone based on the two-dimensional Jensen model and the double-beam laser radar as claimed in claim 1, wherein in the step 1, an SCADA system of the fan is used for collecting the inflow wind speed and the inflow wind direction of the upstream fan, a double-beam laser radar arranged above a cabin of the upstream fan is used for collecting a yaw error angle, a double-beam laser radar arranged above a cabin of the downstream fan is used for collecting wind speed values measured by a left wind measuring point and a right wind measuring point, and the distance and the azimuth angle between the fans in the wind field are determined according to the built site of the wind field.
3. The wake sector fan power prediction method based on the two-dimensional Jensen model and the dual-beam laser radar as claimed in claim 1, wherein in the step 2, according to the distance and the azimuth angle between the two-dimensional Jensen model and the fan, the wake radius of the cross section position where the radar wind measurement point of the downstream fan is located and the wake radius of the cross section position where the wind wheel plane of the downstream fan is located are respectively calculated, and the specific method is as follows:
an included angle theta between the connecting line of the upstream fan and the downstream fan and the central axis of the wake of the upstream fan is definedLComprises the following steps:
θL=0.3CT·β1FWTx
in the formula, thetaxFor incident wind angle, θFWTAzimuth angle, β, of the upstream fan relative to the downstream fan1Is the yaw error angle of the upstream fan, CTIs the lift coefficient of the fan;
defining the vertical distance L between the cross section of a radar wind measuring point of the downstream fan and the upstream fanlThe calculation formula is as follows:
Ll=Ldcos(θL)-z0 cos(α)
where Ld is the distance between the cabin links of the upstream and downstream fans, and z0The distance between the wind point and the radar is measured for the radar of the downstream fan, and alpha is the downstreamThe included angle between the radar laser beam of the fan and the central axis of the downstream fan;
the wake radius calculation formula of the cross section position where the radar wind measurement point of the downstream fan is located is as follows:
Rl=kLl+r0
where k is the wake attenuation coefficient, r0Is the wind wheel radius of the fan.
The vertical distance between the section of the downstream fan wind wheel plane and the upstream fan is defined as LwThe calculation formula is as follows:
Lw=Ldcos(θL)
the calculation formula of the wake radius of the cross section position of the downstream fan wind wheel plane is as follows:
Rw=kLw+r0
4. the method for predicting the fan power in the wake zone based on the two-dimensional Jensen model and the dual-beam laser radar as claimed in claim 1, wherein in the step 3, it is determined whether two wind measuring points of the downstream fan radar are both in the natural wind speed zone, one wind measuring point is in the wake zone and the other wind measuring point is in the natural wind speed zone, or whether two wind measuring points are both in the wake zone, and the method specifically comprises the following steps:
if Ldsin (theta)L)-z0 sin(α)>RlIf the wind speed is higher than the set wind speed, the two wind measuring points of the downstream fan radar are both in a natural wind speed area;
if Ldsin (theta)L)-z0 sin(α)≤Rl&Ldsin(θL)+z0 sin(α)>RlIf the wind speed is higher than the set wind speed, the radar of the downstream fan is positioned in the wake area at one wind measuring point, and the other wind measuring point is positioned in the natural wind speed area at the other wind measuring point;
if Ldsin (theta)L)-z0 sin(α)≤Rl&Ldsin(θL)+z0 sin(α)≤RlIf so, two wind measuring points of the downstream fan radar are both in the wake zone;
wherein Ld is the distance between the connecting lines of the engine rooms of the upstream fan and the downstream fan, and thetaLIs an included angle between the connecting line of the upstream fan and the downstream fan and the central axis of the wake flow of the upstream fan, z0Is downstreamThe distance between the wind point and the radar is measured by the radar of the fan, alpha is the included angle between the radar laser beam of the downstream fan and the central axis of the downstream fan, RlThe radius of the wake flow of the cross section where the radar wind measurement point of the downstream fan is located.
5. The method for predicting the power of the fan in the wake zone based on the two-dimensional Jensen model and the dual-beam laser radar as claimed in claim 1, wherein in the step 4, whether the downstream fan is positioned on the left side or the right side of the central wake axis of the upstream fan is judged according to an included angle between a connecting line of the upstream fan and the downstream fan and the central wake axis of the upstream fan, and the specific method comprises the following steps:
if the included angle between the connecting line of the upstream fan and the downstream fan and the wake central axis of the upstream fan is smaller than 0, the downstream fan is positioned on the left side of the wake central axis, otherwise, the downstream fan is positioned on the right side of the wake central axis.
6. The method for predicting the power of the fan in the wake zone based on the two-dimensional Jensen model and the dual-beam laser radar as claimed in claim 1, wherein in the step 5, the wind speed values measured by the compensated left and right wind measuring points of the downstream fan radar and the yaw error angle of the compensated downstream fan are calculated according to the position of the wind measuring point of the downstream fan radar in the wind speed zone, and the specific method is as follows:
(1) if the downstream fan is positioned on the left side of the central axis of the wake flow, the radar left wind measuring point is positioned in the natural wind speed area, and the right wind measuring point is positioned in the wake flow area, the wind speed compensation coefficient of the right wind measuring point is as follows:
Figure FDA0002665701100000031
in the formula, CTIs the lift coefficient of the fan, k is the wake flow attenuation coefficient, r0Is the radius of the wind wheel of the fan, LlThe vertical distance between the cross section of the radar wind measurement point of the downstream fan and the upstream fan is obtained;
rrkthe vertical distance between the position of the right wind measuring point and the central line of the wake flow plane is calculated by the following formula:
rrk=Ldsin(θL)-z0sin(α)
where Ld is the distance between the engine room connecting lines of the upstream and downstream fans, and θLIs an included angle between the connecting line of the upstream fan and the downstream fan and the central axis of the wake flow of the upstream fan, z0The distance between a radar wind point and a radar of a downstream fan is measured, and alpha is an included angle between a radar laser beam of the downstream fan and a central axis of the downstream fan;
the wind speed values measured by the compensated left and right wind measuring points of the downstream fan radar are as follows:
Va=Vlos1
Vb=Vlos2/Crw
in the formula, Vlos1,Vlos2Actual measured wind speed values of left and right wind measuring points of a radar of a downstream fan are respectively obtained;
(2) if the downstream fan is positioned on the right side of the central axis of the wake flow, the radar left wind measuring point is positioned in the wake flow area, and the radar right wind measuring point is positioned in the natural wind speed area, the wind speed compensation coefficient of the left wind measuring point is as follows:
Figure FDA0002665701100000032
in the formula, rlkThe vertical distance between the position of the left wind measuring point and the central line of the wake flow plane is calculated by the following formula:
rlk=Ldsin(θL)-z0 sin(α)
the wind speed values measured by the compensated left and right wind measuring points of the downstream fan radar are as follows:
Va=Vlos1/Clw
Vb=Vlos2
(3) if the downstream fan is positioned on the left side of the central axis of the wake flow and the radar left and right wind measuring points are positioned in the wake flow area, the wind speed compensation coefficients of the left and right wind measuring points are respectively as follows:
Figure FDA0002665701100000041
Figure FDA0002665701100000042
in the formula, rlkAnd rrkThe vertical distances between the positions of the left wind measuring point and the right wind measuring point and the central line of the wake flow plane are respectively calculated by the following formula:
rlk=Ldsin(θL)+z0sin(α)
rrk=|Ldsin(θL)-z0sin(α)|
the wind speed values measured by the compensated left and right wind measuring points of the downstream fan radar are as follows:
Va=Vlos1/Clw
Vb=Vlos2/Crw
(4) if the downstream fan is positioned on the right side of the central axis of the wake flow and the radar left and right wind measuring points are positioned in the wake flow area, the wind speed compensation coefficients of the left and right wind measuring points are respectively as follows:
Figure FDA0002665701100000043
Figure FDA0002665701100000044
in the formula, rlkAnd rrkThe vertical distances between the positions of the left wind measuring point and the right wind measuring point and the central line of the wake flow plane are respectively calculated by the following formula:
rlk=|Ldsin(θL)-z0sin(α)|
rrk=Ldsin(θL)+z0sin(α)
the wind speed values measured by the compensated left and right wind measuring points of the downstream fan radar are as follows:
Va=Vlos1/Clw
Vb=Vlos2/Crw
(5) if the left wind measuring point and the right wind measuring point of the downstream fan radar are both in the natural wind speed area, the wind speed values measured by the left wind measuring point and the right wind measuring point of the downstream fan radar after compensation are as follows:
Va=Vlos1
Vb=Vlos2
according to the difference of the wind speed areas where the downstream fans are located, calculating the yaw error angle of the compensated downstream fan, wherein the specific formula is as follows:
Figure FDA0002665701100000051
7. the method for predicting the power of the fan in the wake zone based on the two-dimensional Jensen model and the double-beam laser radar according to claim 1, wherein in the step 6, whether the wind wheel plane of the downstream fan is completely in the natural wind speed zone, partially in the wake zone or completely in the wake zone is judged, and the specific method comprises the following steps:
if Ldsin (theta)L)-r1>RwThe plane of the wind wheel of the downstream fan is completely in the natural wind speed area;
if Ldsin (theta)L)-r1≤Rw&Ldsin(θL)+r1>RwThe plane part of the wind wheel of the downstream fan is positioned in the wake area;
if Ldsin (theta)L)-r1≤Rw&Ldsin(θL)+r1≤RwThe plane of the wind wheel of the downstream fan is completely positioned in the wake area;
wherein r is1The calculation formula is that the projection length of the radius of the wind wheel on the central line of the vertical wake flow is as follows:
r1=r0 cos(0.3CT·β12)
wherein Ld is the distance between the connecting lines of the engine rooms of the upstream fan and the downstream fan, and thetaLIs an included angle between the connecting line of the upstream fan and the downstream fan and the central axis of the wake flow of the upstream fan, RwIs the tail of the cross section position of the downstream fan wind wheel planeRadius of flow, r0Is the radius of wind wheel of the fan, beta1Is an upstream fan yaw error angle, beta2For compensated yaw error angle of downstream fan, CTIs the lift coefficient of the fan.
8. The wake sector fan power prediction method based on the two-dimensional Jensen model and the dual-beam laser radar as claimed in claim 1, wherein in the step 7, the equivalent inflow wind speed of the downstream fan is calculated according to the difference that the wind wheel plane of the downstream fan is positioned in the wake sector, and the specific method is as follows:
(1) if the plane of the wind wheel of the downstream fan is completely in the natural wind speed area, the equivalent inflow wind speed of the downstream fan is as follows:
Figure FDA0002665701100000052
wherein u is0The inflow wind speed of the upstream fan.
(2) If the plane part of the downstream fan wind wheel is positioned in the wake area, the area of the downstream fan wind wheel plane in the wake area is as follows:
Figure FDA0002665701100000061
wherein R iswIs the wake radius of the cross section position of the downstream fan wind wheel plane, r1Is the projection length of the radius of the wind wheel on the central line of the vertical wake;
let O be1,O2Is the wake flow center and the downstream fan wheel center of the downstream fan wheel plane, B is any intersection point of the wake flow area edge and the downstream fan wheel edge of the downstream fan wheel plane, and theta1Is a connecting line O1O2And a connection line O1Angle of B, θ2Is a connecting line O1O2And a connection line O2The calculation formula of the included angle B is as follows:
Figure FDA0002665701100000062
Figure FDA0002665701100000063
wherein L isRIs the vertical distance of the downstream fan center from the upstream fan wake centerline:
LR=Ldsin(θL)
wherein Ld is the distance between the connecting lines of the engine rooms of the upstream fan and the downstream fan, and thetaLThe included angle between the connecting line of the upstream fan and the downstream fan and the central axis of the wake flow of the upstream fan is formed;
the wind wheel area of the downstream fan is as follows:
S=πr0 2
wherein r is0Is the radius of the wind wheel of the fan.
The projected area of the wind wheel on the central line of the vertical wake flow is as follows:
S1=πr1 2
wherein r is1The length of the projection of the radius of the wind wheel on the central line of the vertical wake flow;
wind speed u in wake zonerThe calculation formula is as follows:
Figure FDA0002665701100000064
in the formula, CTIs the lift coefficient of the fan, k is the wake flow attenuation coefficient, r0Is the radius of the wind wheel of the fan, LwThe vertical distance between the cross section of the downstream fan wind wheel plane and the upstream fan, r is the vertical distance between any point of the downstream fan wind wheel plane and the center line of the tail flow plane, u is the vertical distance between any point of the downstream fan wind wheel plane and the center line of the tail flow plane0The inflow wind speed of an upstream fan;
the equivalent inflow wind speed calculation formula of the downstream fan is as follows:
Figure FDA0002665701100000065
wherein, a and b are integral upper and lower limits, and the values are as follows:
a=LR-r1
b=Rw
(3) if the plane part of the wind wheel of the downstream fan is positioned in the wake area, the equivalent inflow wind speed calculation formula of the downstream fan is as follows:
Figure FDA0002665701100000071
wherein,
a=LR-r1
b=LR+r1
9. the method for predicting the power of the fan in the wake zone based on the two-dimensional Jensen model and the dual-beam laser radar as claimed in claim 1, wherein in the step 8, the power output value of the downstream fan is calculated by a specific formula:
Figure FDA0002665701100000072
wherein rho is air density, S is downstream fan wind wheel area, CpFor the downstream fan power utilization factor,
Figure FDA0002665701100000073
is the equivalent inflow wind speed of a downstream fan, beta2The yaw error angle of the downstream fan after compensation.
10. Two-dimensional Jensen model and double-beam laser radar-based wake flow area fan power prediction system is characterized by comprising:
the data acquisition module is used for acquiring the inflow wind speed and the inflow wind direction of the upstream fan, acquiring the yaw error angle of the upstream fan, acquiring the wind speed values measured by the left and right wind measuring points of the dual-beam laser radar of the downstream fan, and determining the distance and the azimuth angle between the fans in a wind field;
the wake flow radius calculation module is used for calculating the wake flow radius of the cross section position of the radar wind measuring point of the downstream fan and the wake flow radius of the cross section position of the wind wheel plane of the downstream fan according to the two-dimensional Jensen model and the space and the azimuth angle between the fans;
the wind measuring point position judging module is used for judging whether two wind measuring points of the downstream fan radar are both in a natural wind speed area, one wind measuring point is in a wake area, the other wind measuring point is in the natural wind speed area, or the two wind measuring points are both in the wake area;
the fan position judging module is used for judging whether the downstream fan is positioned on the left side or the right side of the wake central axis according to an included angle between the connecting line of the upstream fan and the downstream fan and the wake central axis of the upstream fan;
the compensation module is used for calculating wind speed values measured by the compensated radar left and right wind measuring points of the downstream fan and a yaw error angle of the compensated downstream fan according to the position of the radar wind measuring point in the wind speed area;
the wind wheel plane position judging module is used for judging whether the wind wheel plane of the downstream fan is completely positioned in a natural wind speed area, a partial wake area or a wake area;
the equivalent inflow wind speed calculation module is used for calculating 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 flow area;
and the power output calculation module is used for calculating the power output value of the downstream fan.
CN202010917961.4A 2020-09-03 2020-09-03 Two-dimensional Jensen model and double-beam laser radar-based wake region fan power prediction method Active CN112149363B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
袁明友等: "基于二维尾流模型的风电场产能分析", 《可再生能源》 *

Cited By (4)

* Cited by examiner, † Cited by third party
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