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

    • 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)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
  • Optical Radar Systems And Details Thereof (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

Two-dimensional Jensen model and double-beam laser radar-based wake region fan power prediction method
Technical Field
The invention belongs to the technical field of fan control, and particularly relates to a wake flow fan power prediction method based on a two-dimensional Jensen model and a double-beam laser radar.
Background
Wind energy is one of the most promising renewable energy sources in the world as a rich and clean new energy source. With the increase of the transmission and accommodation capacity of a power grid, wind power plant owners pay more attention to the improvement of the power generation efficiency of the fans. The wake flow of the upstream fan can obviously reduce the power output of the downstream fan, so that the power output value of the downstream fan in the wake area of the upstream fan is accurately calculated, and the method has important significance for controlling and optimizing the fans in the wind power plant and improving the overall power generation efficiency of the wind power plant. At present, a method for calculating the power of a fan considering wake flow exists. However, under the influence of the yaw error angle of the upstream and downstream fans, the method is inaccurate in judging the position of the wake zone where the downstream fan is located, and is inaccurate in calculating the power output value when the downstream fan is located in the wake zone.
Disclosure of Invention
The invention aims to provide a wake region fan power prediction method based on a two-dimensional Jensen model and a double-beam laser radar.
The technical solution for realizing the purpose of the invention is as follows: a method for predicting the power of a fan in a wake flow area based on a two-dimensional Jensen model and a double-beam laser radar comprises the following specific steps:
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 whether 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 flow area;
and 8, calculating the power output value of the downstream fan.
Further, in the step 1, an SCADA system of the fan is used for collecting inflow wind speed and inflow wind direction of an upstream fan, a double-beam laser radar arranged above an upstream fan cabin is used for collecting yaw error angles, a double-beam laser radar arranged above a downstream fan cabin 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 building site selection of the wind field.
Further, in step 2, according to the two-dimensional Jensen model and the distribution distance and the distribution angle of the wind turbine, respectively calculating a wake flow radius of a cross section position where a radar wind measurement point of the downstream wind turbine is located and a wake flow radius of a cross section position where a wind wheel plane of the downstream wind turbine is located, 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·β1+θFWT-θx
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)-z0cos(α)
where Ld is the distance between the cabin links of the upstream and downstream fans, and z0For downstream fansThe distance between the wind measuring point and the radar is reached, and alpha is the included angle between the radar laser beam of the downstream fan and the central axis of the downstream fan;
the calculation formula of the wake radius of the cross section 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)
then, the calculation formula of the wake radius of the cross section position where the wind wheel plane of the downstream fan is located is as follows:
Rw＝kLw+r0
further, in step 3, it is determined 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 whether two wind measuring points are both in the wake area, the specific method is as follows:
if Ldsin (theta)L)-z0sin(α)＞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)-z0sin(α)≤Rl&Ldsin(θL)+z0sin(α)＞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)-z0sin(α)≤Rl&Ldsin(θL)+z0sin(α)≤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, z0The distance between the radar wind point and the radar is measured for the radar of the downstream fan, and alpha is the clamp between the radar laser beam of the downstream fan and the central axis of the downstream fanCorner, RlThe radius of the wake flow of the cross section where the radar wind measurement point of the downstream fan is located.
Further, in step 4, according to an included angle between a connecting line of the upstream fan and the downstream fan and a wake central axis of the upstream fan, whether the downstream fan is positioned on the left side or the right side of the wake central axis is judged, 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.
Further, in step 5, according to the position of the radar wind measurement point in the wind speed area, calculating the wind speed values measured by the radar left and right wind measurement points of the compensated downstream fan and the yaw error angle of the compensated downstream fan, 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 BDA0002665701110000031
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 BDA0002665701110000041
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)-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
(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 BDA0002665701110000042
Figure BDA0002665701110000043
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 BDA0002665701110000051
Figure BDA0002665701110000052
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 BDA0002665701110000053
further, in step 6, it is determined 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, and the specific method is as follows:
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＝r0cos(0.3CT·β1+β2)
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 wake radius of the cross section position of the downstream fan wind wheel plane, 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.
Further, according to the position of the downstream fan wind wheel plane in the wake area, calculating the equivalent inflow wind speed of the downstream fan, 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 BDA0002665701110000061
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 BDA0002665701110000062
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 BDA0002665701110000063
Figure BDA0002665701110000064
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 wind wheel radius.
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 BDA0002665701110000071
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 fan wind wheel plane and the center line of the tail flow plane, u0The inflow wind speed of an upstream fan;
the equivalent inflow wind speed calculation formula of the downstream fan is as follows:
Figure BDA0002665701110000072
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 BDA0002665701110000073
wherein,
a＝LR-r1
b＝LR+r1
further, in step 8, calculating a power output value of the downstream fan, wherein a specific formula is as follows:
Figure BDA0002665701110000074
wherein rho is air density, S is downstream fan wind wheel area, CpFor the downstream fan power utilization factor,
Figure BDA0002665701110000075
is the equivalent inflow wind speed of a downstream fan, beta2The yaw error angle of the downstream fan after compensation.
A wake flow area fan power prediction system based on a two-dimensional Jensen model and a double-beam laser radar comprises:
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.
Compared with the prior art, the invention has the following remarkable advantages: based on a two-dimensional Jensen model, the yaw error angle of the downstream fan after compensation in the wake area is accurately calculated by combining the change of the yaw error angle of the upstream fan to the position of the wake area, and the equivalent inflow wind speed of the downstream fan in the wake area of the upstream fan is calculated on the basis of the yaw error angle, so that the accurate calculation of the power output value of the wake area of the downstream fan is finally realized, and a foundation is laid for the control optimization of the wake suppression in the wind power plant.
Drawings
FIG. 1 is a flow chart of a method for predicting the fan power in the wake zone based on a two-dimensional Jensen model and a double-beam laser radar.
Fig. 2 is a schematic structural diagram of a dual-beam lidar.
FIG. 3 is a schematic diagram of the present invention for determining the wind speed zone of the wind wheel plane of the downstream wind turbine.
FIG. 4 is a schematic diagram of calculating the area of a downstream fan rotor plane in a wake area according to the present invention.
FIG. 5 is a plot of the compensated downstream fan yaw error angle of the present invention.
FIG. 6 is a graph comparing the predicted power output value of the wake region of the downstream fan with the actual power output value.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
As shown in fig. 1, the method for predicting the power of the fan in the wake area based on the two-dimensional Jensen model and the dual-beam laser radar specifically comprises the following steps:
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;
an SCADA system of the wind turbine is used for collecting inflow wind speed and inflow wind direction of an upstream wind turbine, a double-beam laser radar (shown in figure 2) arranged above an upstream wind turbine cabin is used for collecting yaw error angles, a double-beam laser radar (shown in figure 2) arranged above a downstream wind turbine cabin 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 wind turbine in the wind field are determined according to the building site selection of the wind field.
Step 2, calculating the wake flow radius of the cross section position of the radar wind measuring point of the downstream fan according to the two-dimensional Jensen model and the distribution distance and the distribution angle of the fan, and calculating the wake flow radius of the cross section position of the wind wheel plane of the downstream fan;
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·β1+θFWT-θx
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)-z0cos(α)
where Ld is the distance between the cabin links of the upstream and downstream fans, and 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;
then, the calculation formula of the wake radius of the cross section where the radar wind measuring point 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)
then, the calculation formula of the wake radius of the cross section position where the wind wheel plane of the downstream fan is located is as follows:
Rw＝kLw+r0
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 whether the two wind measuring points are both in the wake area;
if Ldsin (theta)L)-z0sin(α)＞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)-z0sin(α)≤Rl&Ldsin(θL)+z0sin(α)＞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)-z0sin(α)≤Rl&Ldsin(θL)+z0sin(α)≤RlIf so, two wind measuring points of the downstream fan radar are both in the wake zone;
step 4, 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;
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.
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;
according to the position of the radar wind measuring point in the wind speed area, calculating the wind speed values measured by the compensated radar left and right wind measuring points of the downstream fan, and dividing the values into five specific situations:
(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 BDA0002665701110000101
in the formula, CTIs the lift coefficient of the fan;
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(α)
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 BDA0002665701110000111
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)-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
(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 BDA0002665701110000112
Figure BDA0002665701110000113
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 BDA0002665701110000114
Figure BDA0002665701110000115
in the formula, rlkAnd rrkRespectively the 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 planeThe calculation formula is as follows:
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 BDA0002665701110000121
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, wherein the schematic diagram is shown in fig. 3;
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＝r0cos(0.3CT·β1+β2)
step 7, calculating the equivalent inflow wind speed of the downstream fan according to the difference that the wind wheel plane of the downstream fan is positioned in the wake flow area;
according to the position of the downstream fan wind wheel plane in the wake area, the method is divided into three specific situations:
(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 BDA0002665701110000123
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 BDA0002665701110000122
as shown in FIG. 4, let O1，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 BDA0002665701110000131
Figure BDA0002665701110000132
wherein L isRIs the distance from the downstream fan center to the upstream fan wake centerlineVertical distance of (d):
LR＝Ldsin(θL)
the wind wheel area of the downstream fan is as follows:
S＝πr0 2
the projected area of the wind wheel on the central line of the vertical wake flow is as follows:
S1＝πr1 2
wind speed u in wake zonerThe calculation formula is as follows:
Figure BDA0002665701110000133
the equivalent inflow wind speed calculation formula of the downstream fan is as follows:
Figure BDA0002665701110000134
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 BDA0002665701110000135
wherein,
a＝LR-r1
b＝LR+r1
step 8, calculating the power output value of the downstream fan;
the calculation formula is as follows:
Figure BDA0002665701110000141
the invention also provides a system for predicting the power of the fan in the wake zone based on the two-dimensional Jensen model and the double-beam laser radar, which comprises the following steps:
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.
A computer apparatus comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing a wake sector fan power prediction based on a two-dimensional Jensen model and a dual-beam lidar when executing the computer program.
A computer-readable storage medium, having stored thereon a computer program which, when executed by a processor, enables wake sector fan power prediction based on a two-dimensional Jensen model and a dual-beam lidar.
Examples
In order to verify the effectiveness of the scheme, the following simulation experiment is carried out on a wind turbine with the number of A04 in a certain wind farm. Two fans A03 and A05 are arranged near the fan A04, the azimuth angle of the fan A03 relative to the fan A04 is 35 degrees, the cabin connecting line distance is 316m, the azimuth angle of the fan A05 relative to the fan A04 is 218 degrees, and the cabin connecting line distance is 314 m; the power utilization coefficients of the three fans are all 0.55, the lift coefficients are all 0.5, the radiuses of the wind wheels are all 56.5m, the distance between a radar wind measuring point and a radar is 80m, the included angle between laser emitted by the radar and the central axis of the fan is 30 degrees, and the wake flow attenuation coefficient is 0.075.
With the fan a04 as the downstream fan, the yaw error angle of the compensated downstream fan a04 is shown in fig. 5, and it can be seen that the corrected yaw error angle of the fan a04 successfully eliminates the interference of the wake effect on the yaw error calculation, so that the distortion degree of the yaw error value in the wake area is obviously reduced. And calculating a corresponding wind speed inflow angle when the wind wheel plane of the fan A04 is in the wake zone based on the corrected yaw error angle of the downstream fan A04. In the space 360 degrees, when the fan A03 is an upstream fan of the fan A04, in a wake flow area generated by the fan A03, the wake flow radius at the wind wheel plane of the fan A04 is 68.8m, the wind direction angle of the fan A04 completely positioned in the wake flow area is 29-44 degrees, the wind direction angle of the fan A04 partially positioned in the wake flow area is 12-28 degrees, 45-62 degrees, and the wind wheel plane of the fan A04 under the rest wind direction angles is completely positioned in a natural wind speed area. In the 360 degrees of space, when fan A05 is the upstream fan of fan A04, in the wake area generated by fan A05, the wake radius at the rotor plane of fan A04 is 68.8m, the wind direction angle of the wind wheel plane of fan A04 completely positioned in the wake area is 212-225 degrees, the wind direction angle of the wind wheel plane of fan A04 partially positioned in the wake area is 193-211 degrees, 226-242 degrees, and the wind wheel plane of fan A04 under the rest wind direction angles is completely positioned in the natural wind speed area.
And calculating the equivalent inflow wind speed of the fan A04 according to the calculated wind speed inflow angle corresponding to the wake area of the wind wheel plane of the fan A04, and finally calculating to obtain the power output value of the wake area of the fan A04, as shown in FIG. 6. And comparing the calculated power output value of the wake area with the actual power to obtain the comprehensive error of the predicted power and the actual power of 4.26 percent.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory, among others. Non-volatile memory can include read-only memory (ROM), Programmable ROM (PROM), Electrically Programmable ROM (EPROM), Electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Link DRAM (SLDRAM), Rambus Direct RAM (RDRAM), direct bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent 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·β1+θFWT-θx
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·β1+β2)
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
CN112096576B (en) Method for improving annual generated energy of multiple fan arrays based on wake field optimization control
AU2017352549B2 (en) Method and device for calculating power generation of wind farm
CN111794909B (en) Sector regulation-oriented wind farm level yaw dynamic optimization method and system
CN105160060B (en) A kind of wind power plant theoretical power (horse-power) based on actual power curve matching determines method
CN108953060B (en) Wind power plant field level yaw control method based on laser radar anemoscope
AU2018327223B2 (en) Method and apparatus for self-adaption of a cut-out strategy
CN114169614B (en) Wind power plant optimal scheduling method and system based on wind turbine wake model optimization
CN111894802B (en) Fan wake zone yaw error correction method based on Jensen model and double-beam radar
KR101383792B1 (en) Method for calibration of nacelle anemometer using lidar measurements
CN112434252B (en) Wake flow calculation method considering local environment factors of wind power plant
CN113219478B (en) Method, device and storage medium for measuring and predicting wind speed and wind direction in fan impeller range
CN113536487A (en) Bivariate Gaussian function-based wake flow calculation method and device and storage medium
CN114707437A (en) Method for establishing three-dimensional full wake flow model of horizontal axis wind turbine
CN111120202B (en) Yaw angle adjusting method, device, medium and equipment of wind generating set
CN113236487B (en) Wind power plant noise control method, system, device and readable storage medium
CN110765640B (en) Method, system and equipment for calculating effective wind speed of rotor
CN109766611B (en) Wind farm power simplified prediction method considering terrain gradient
CN114294157A (en) Wind power fan yaw deviation identification method, deviation correction method, electronic device and medium
Dai et al. Study on Obtaining Real Power Curve of Wind Turbines Using SCADA Data
Hwang et al. Wind farm layout optimization using multidisciplinary model
CN118428280B (en) Wind turbine wake flow speed and turbulence intensity prediction method considering atmospheric stability
Ahmad et al. Analysis of two onshore wind farms with a dynamic farm controller
CN115443608A (en) Method for optimizing power generation capacity of photovoltaic module
CN117709942B (en) Fan field clustering distributed operation and maintenance method based on wake flow analysis

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

CI02 Correction of invention patent application
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

GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract
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