CN114154296A - Fully developed wind power plant equivalent roughness calculation method considering atmospheric stability - Google Patents

Fully developed wind power plant equivalent roughness calculation method considering atmospheric stability Download PDF

Info

Publication number
CN114154296A
CN114154296A CN202111242620.2A CN202111242620A CN114154296A CN 114154296 A CN114154296 A CN 114154296A CN 202111242620 A CN202111242620 A CN 202111242620A CN 114154296 A CN114154296 A CN 114154296A
Authority
CN
China
Prior art keywords
wind
power plant
wind power
equivalent roughness
atmospheric stability
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.)
Pending
Application number
CN202111242620.2A
Other languages
Chinese (zh)
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.)
Three Gorges New Energy Dachaidan Wind Power Co ltd
North China Electric Power University
Original Assignee
Three Gorges New Energy Dachaidan Wind Power Co ltd
North China Electric Power University
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 Three Gorges New Energy Dachaidan Wind Power Co ltd, North China Electric Power University filed Critical Three Gorges New Energy Dachaidan Wind Power Co ltd
Priority to CN202111242620.2A priority Critical patent/CN114154296A/en
Publication of CN114154296A publication Critical patent/CN114154296A/en
Pending legal-status Critical Current

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/06—Wind turbines or wind farms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Wind Motors (AREA)

Abstract

The invention relates to the technical field of wind power plant equivalent modeling and wind power plant micro site selection, and particularly provides a wind power plant equivalent roughness calculation method considering atmospheric stability and aiming at solving the problem that the existing method is not perfect in calculation of wind power plant equivalent roughness. To this end, the method of the invention comprises: acquiring parameters of a wind turbine generator; obtaining an influence factor of interaction of atmospheric stability and a wind power plant on equivalent roughness; acquiring a dimensionless wake flow additional vortex-viscosity coefficient of the wind power plant; and calculating the equivalent roughness of the wind power plant according to the parameters of the wind turbine generator, the influence factors and the dimensionless wake flow additional vortex-viscosity coefficient of the wind power plant. According to the method, the influence of the atmospheric stability is considered, and the reaction of the wind power plant to the atmospheric stability is considered, so that the equivalent roughness of the boundary layer of the wind power plant under different conditions can be quickly, simply and accurately calculated, and support is provided for the optimization design of the large wind power plant.

Description

Fully developed wind power plant equivalent roughness calculation method considering atmospheric stability
Technical Field
The invention relates to the technical field of fully developing wind power plant equivalent modeling and wind power plant micro site selection, and particularly provides a wind power plant equivalent roughness calculation method considering atmospheric stability.
Background
Roughness is one measure of the amount of friction of the ground to the wind, and in general, the flatter and smoother the ground, the less roughness. In a logarithmic wind profile, roughness represents the height at which the wind speed is zero. If the wind turbine generator is regarded as a rough element, the roughness is increased due to the existence of the wind power plant, and the roughness at the moment is called wind power plant equivalent roughness. The wind power plant equivalent roughness is actually equivalent to a special 'terrain' in the wind power plant, and the wind power plant equivalent roughness is influenced by the wind turbine generator arrangement distance, the wind turbine diameter, the hub height, the running state of the wind turbine generator, the original ground roughness and the original atmospheric stability. For a large wind power plant, the wind power plant equivalent roughness is accurately predicted, and the method has very important significance for micro site selection and power prediction of the wind power plant and improvement of the economic benefit of the wind power plant.
The calculation of the equivalent roughness of the wind power plant is very complex, but for a large wind power plant, most of the downstream flow state of the wind power plant tends to be stable, and the wind power plant can be approximately regarded as a fully developed wind power plant. This provides convenience for establishing the wind power plant equivalent roughness model. The strong interaction between a large wind power plant and an atmospheric boundary layer is an important factor influencing the accurate evaluation of roughness, and the interaction can be described as follows: on one hand, the complex wake effect of the wind power plant can destroy the original momentum balance of the atmospheric boundary layer, obviously influence the vertical momentum transport of the atmospheric boundary layer and change the stability of the atmospheric boundary layer. On the other hand, the atmospheric stability can change the original turbulence characteristic of the atmospheric boundary layer, and the wake effect of the wind turbine generator is adversely affected. The two are mutually influenced and continuously iterated, and finally a balanced state is achieved.
Atmospheric stability is not always maintained in neutral conditions during actual wind farm operation. From the time scale of the hour scale, the atmospheric stability is in the unstable to stable to unstable state respectively from the noon of the first day to midnight to the noon of the second day. Atmospheric stability at a location typically deviates from neutrality, even on a long-term statistical average. Therefore, the equivalent roughness of the wind power plant under different atmospheric stability degrees needs to be evaluated.
In commercial wind resource analysis software WAsP, a simple Lettau model calculated based on wind farm visualization parameters (wind wheel diameter, hub height, average wind turbine generator footprint) is adopted.
The E & F model divides the wind power plant atmospheric boundary layer under the neutral condition into two stress layers along the vertical direction, and accordingly two logarithmic law average speed expressions are obtained.
The Calaf (top down) model considers three stress layer structures on the basis of the E & F model, and the speed distribution of the Calaf (top down) model is more consistent with the large vortex simulation result of the fully developed wind power plant.
The P & R model considers the atmospheric stability correction on the basis of the E & F model, but does not consider the reaction of the wind power plant to the atmospheric stability and three stress layer structures of a wind power plant boundary layer.
The calculation of the equivalent roughness of the wind power plant by the model is not perfect.
Disclosure of Invention
The invention aims to solve the technical problems, namely, the problem that the existing model is not perfect in calculating the equivalent roughness of the wind power plant is solved.
In a first aspect, the invention provides a wind farm equivalent roughness calculation method taking atmospheric stability into account. The method comprises the following steps:
acquiring parameters of a wind turbine generator;
obtaining an influence factor of interaction of atmospheric stability and a wind power plant on equivalent roughness;
acquiring a dimensionless wake flow additional vortex-viscosity coefficient of the wind power plant;
and calculating the equivalent roughness of the wind power plant according to the parameters of the wind turbine generator, the influence factors and the dimensionless wake flow additional vortex-viscosity coefficient of the wind power plant.
In a preferred embodiment of the wind farm equivalent roughness calculating method, the step of calculating the equivalent roughness of the wind farm according to the parameters of the wind turbine, the influence factors and the dimensionless wake flow additional vortex-viscosity coefficient of the wind farm specifically includes calculating the equivalent roughness of the wind farm according to the following equation:
Figure BDA0003320024250000031
wherein D is the diameter of the wind wheel of the wind turbine generator, zhThe hub height of the wind turbine generator;
w(Lhi) Is the influence factor of the interaction of the atmospheric stability above the height of the hub and the wind power plant on the equivalent roughness, w (L)lo) The influence factor of the interaction of the atmospheric stability below the height of the hub and the wind power plant on the equivalent roughness is shown;
Figure BDA0003320024250000032
attaching vortex viscosity coefficients to dimensionless wake flows of the wind farm,
Figure BDA0003320024250000033
CTis the thrust coefficient of the wind wheel, SxAnd SyThe parameters are dimensionless parameters of the flow direction and the span-wise distance of the wind turbine generator relative to the diameter of the wind wheel.
In a preferred embodiment of the wind farm equivalent roughness calculating method, the wind farm equivalent roughness is calculated by using the wind farm equivalent roughness and the atmospheric stability above the hub height as an influence factor w (L) of the interaction between the wind farm equivalent roughness and the atmospheric stabilityhi) Determined by the following equation:
Figure BDA0003320024250000034
wherein the content of the first and second substances,
Figure BDA0003320024250000035
and
Figure BDA0003320024250000036
respectively, are atmospheric stability correction functions at different heights above the hub height.
In a preferred embodiment of the wind farm equivalent roughness calculating method, the wind farm equivalent roughness is calculated by using the wind farm equivalent roughness and the atmospheric stability below the hub height as an influence factor w (L)lo) Determined by the following equation:
Figure BDA0003320024250000037
wherein the content of the first and second substances,
Figure BDA0003320024250000041
and
Figure BDA0003320024250000042
respectively, the atmospheric stability correction functions at different heights below the hub height.
In a preferred embodiment of the wind farm equivalent roughness calculating method, a dimensionless wake flow of the wind farm is added with a vortex viscosity coefficient
Figure BDA0003320024250000043
Determined by the following equation:
Figure BDA0003320024250000044
where κ is the Karman constant, u*Is the friction speed;
Figure BDA0003320024250000045
(zh) is the wind speed averaged over time and space for the hub height plane.
In a preferred embodiment of the wind farm equivalent roughness calculation method described above, the wind farm is a fully developed area of a large wind farm.
In a second aspect, the present invention provides a method of fully developing a wind farm generated power estimate, the method comprising the steps of:
calculating the equivalent roughness of the wind power plant according to the wind power plant equivalent roughness calculation method considering the atmospheric stability;
determining the wind speed for fully developing the height of the hub of the wind power plant according to the calculated equivalent roughness;
and according to the wind speed of the hub height of the wind power plant, performing linear interpolation on the wind speed-power curve of the wind turbine generator to estimate the generating power of the wind power plant.
In a preferred embodiment of the above wind farm generated power estimation method, the wind speed at which the hub height of the wind farm is fully developed is determined by the following equation:
Figure BDA0003320024250000046
wherein z is0，hiTo an equivalent roughness, u*hiD is the diameter of the wind wheel of the wind turbine generator, z for fully developing the friction speed of the wind power plant above the hub heighthHub height of wind turbine, w (L)hi) Is an influence factor of the interaction of the atmospheric stability above the height of the hub and the wind power plant on the equivalent roughness,
Figure BDA0003320024250000047
vortex viscosity coefficients are added to the dimensionless wake of the wind farm, and κ is the karman constant.
In a preferred embodiment of the wind farm generated power estimation method, the friction speed u above the height of the hub inside the wind farm is obtained*hiDetermined by the following equation:
Figure BDA0003320024250000051
wherein u is*inflowIs the inflow friction velocity, delta is the boundary layer height, Z0，loRepresenting the roughness of the ground, Z0，hiRepresenting wind farm equivalent roughness, LinflowLength of Monin-obuff, L, of inflow profilehiIs the length psi of the Morin-obufhoff above the hub heightmIs a dimensionless stability parameter.
In a third aspect, the present invention provides a wind farm micro-siting method, comprising the steps of:
carrying out primary micro site selection on the wind power plant;
estimating the generation power of the wind turbines of a fully developed wind park according to the method of claim 7 or 8, based on the parameters of the preliminary micro-siting;
calculating the sum of the generated power of all the wind turbine generators, judging whether the sum of the generated power meets the design requirement, if so, selecting the site to pass, and if not, readjusting the site selection parameter until the sum of the generated power meets the design requirement.
In an embodiment of the method for estimating the generated power of the wind farm, the step of performing preliminary micro-site selection on the wind farm includes: and (4) performing micro site selection on the wind power plant by adopting a Jensen wake flow model and a wake flow square sum superposition model.
The invention achieves the following beneficial effects:
according to the wind power plant equivalent roughness model, the atmospheric stability, the interaction with the wind power plant and the three-stress-layer structure are considered, a more complete wind power plant equivalent roughness model is established, the equivalent roughness of a wind power plant boundary layer under different conditions can be rapidly, simply and accurately calculated through the model, and support is provided for the optimization design of a large wind power plant.
Drawings
Preferred embodiments of the present invention are described below with reference to the accompanying drawings, in which:
FIG. 1 is a flow chart of main steps of a wind power plant equivalent roughness calculation method considering atmospheric stability of the invention;
FIG. 2 is a schematic distribution of an inflow wind profile taking into account atmospheric stability;
FIG. 3 is a comparison graph of wind farm equivalent roughness model calculation results and LES results for different atmospheric stability;
FIG. 4 is a relative error diagram of wind farm equivalent roughness model calculation results and LES results with different atmospheric stability.
Detailed Description
For convenience of describing the technical scheme of the present invention, technical terms involved in the present invention are first explained as follows:
atmospheric stability: that is, the degree of atmospheric thermal stability refers to the ability of an air micelle to maintain its original state after the air micelle is disturbed. General atmospheres can be divided into three categories, stable, neutral and unstable, among which: unstable atmospheres may also be called convective atmospheres. The atmospheric stability varies continuously from unstable to neutral to unstable.
Unstable atmosphere: after the disturbance, if the air micelles accelerate away from their original position, they are referred to as unstable atmosphere, or convective atmosphere.
And (3) stabilizing the atmosphere: if the disturbance tends to return to the original position, the atmosphere is called a stable atmosphere.
Neutral atmosphere: after the disturbance, the air is called neutral atmosphere if it neither accelerates away nor returns to its original position.
Large-scale wind power plants: the wind power station with large flow direction distance is not specially specified, and the general large wind power base can be regarded as a large wind power station.
Fully developing a wind power plant: also called an infinite wind farm, refers to the latter half of a large wind farm. The main characteristics are two: 1) with the increase of the depth of the wind power plant, the power of the wind turbine tends to be stable. 2) The energy source is primarily the downward momentum transport at the top of the boundary layer.
Wind profile: the variation curve of the wind speed along with the height from the ground under a certain roughness of the ground.
Roughness of the ground: refers to the "average height" of the obstacle above the ground.
The method of the present invention will be described in detail below with reference to the accompanying drawings.
As shown in FIG. 1, the wind farm equivalent roughness calculation method considering the atmospheric stability comprises the following steps:
s100, acquiring parameters of the wind turbine generator;
s200, obtaining an influence factor of interaction of atmospheric stability and a wind power plant on equivalent roughness;
s300, acquiring a dimensionless wake flow additional vortex-viscosity coefficient of the wind power plant;
s400, calculating the equivalent roughness of the wind power plant according to the parameters of the wind turbine generator, the influence factors and the dimensionless wake flow additional vortex viscosity coefficient of the wind power plant.
The method not only considers the influence of the atmospheric stability, but also considers the adverse effect of the wind power plant on the atmospheric stability, so that the equivalent roughness of the boundary layer of the wind power plant under different conditions can be quickly, simply and accurately calculated, and support is provided for the optimal design of the large wind power plant.
Specifically, the influence of the atmospheric stability on the boundary layer wind profile is considered at first, a top down boundary layer model is supplemented, the coupling effect of the wind power plant on the atmospheric boundary layer stability is considered, and an equivalent roughness model suitable for the wind power plant with different atmospheric stability is deduced. The method can quickly, simply and accurately calculate the equivalent roughness of the boundary layer of the wind power plant under different conditions, and provides reference for the optimization design of the large wind power plant.
The concrete description is as follows:
the intervals of the wind generating sets in the wind power plant along the flow direction and the spreading direction are respectively SxD and SyD, D is the diameter of the wind wheel, SxAnd SyRespectively, the dimensionless parameters of the flow direction and the span-wise distance of the wind turbine generator set relative to the diameter of the wind wheel, and the height of the hub is zh。
Assuming that the wind power plant has three stress layers, the average flow velocity under different vertical heights can be obtained, and the specific calculation is as follows:
step 1: the internal stress balance (including ground friction stress and retardation stress generated by a wind turbine generator) of a boundary layer of the wind power plant is fully developed, and the equivalent friction stress of the wind power plant is equal to the sum of the ground viscous bottom layer friction stress and the wind field resistance in unit area:
Figure BDA0003320024250000081
wherein ρ air density; u. of*loIs the ground friction speed; u. of*hiIs the wind farm friction speed;
Figure BDA0003320024250000082
Figure BDA0003320024250000083
CTis the thrust coefficient of the wind wheel;
Figure BDA0003320024250000084
(zh) represents the wind speed over the time and space average of the hub height plane.
Step 2: passing through two different heights z according to non-interfering meteorological data1、z2(z1＜z2) Calculating the gradient checking number according to the temperature difference and the wind speed difference, and solving the formula as follows:
Figure BDA0003320024250000085
wherein g is the acceleration of gravity (m/s)2)，
Figure BDA0003320024250000086
Is the mean absolute temperature of the gas layer, Δ T and
Figure BDA0003320024250000087
respectively two heights z of boundary layer1、z2Temperature difference and flow direction wind speed difference betweendIs the dry adiabatic desuperheating rate.
And step 3: calculating the length L of the inflow wind profile according to the calculation result of the step 2inflowThe calculation formula is as follows:
Figure BDA0003320024250000088
wherein the content of the first and second substances,
Figure BDA0003320024250000089
representing the average geometric height.
And 4, step 4: introducing dimensionless atmospheric stability correction of wind speed gradient according to the similar theory of the Morin-obufh:
Figure BDA0003320024250000091
wherein, k is a karman constant, k is 0.4; z is the height from the ground;
Figure BDA0003320024250000092
(z) represents the wind speed averaged over time and space in a plane of height z;
Figure BDA0003320024250000093
and
Figure BDA0003320024250000094
is a dimensionless stability parameter, and the two have the following relationship:
Figure BDA0003320024250000095
u*by friction speed, for an inflow profile u*＝u*inflowFor fully developed wind farm interior, u is above the hub height*＝u*hiU below the hub height*＝u*lo；
L is the length of Monin-obufh, and L is L for inflow wind profileinflowFor fully developed wind power plant interior, L is equal to L above the hub heighthiL is less than the height of the hublo。
Z0Roughness is substituted, and Z is taken below the height of the hub in the fully developed wind power plant0＝Z0，loRepresenting the roughness of the ground; above the hub height Z0，hiAnd representing the equivalent roughness of the wind power plant.
Unless otherwise specified, u hereinafter*L and Z0Should be so treated.
And 5: to formula (4) from z0To the z integral and neglecting the small term available inflow profile:
Figure BDA0003320024250000096
step 6: according to the meteorological data in the step 2, calculating the inflow friction speed u through the logarithmic wind profile corrected by the atmospheric stability*inflowAnd ground roughness Z0，loThe solving formula is:
Figure BDA0003320024250000097
Figure BDA0003320024250000098
note that the inflow friction velocity u may be obtained by fitting the actually measured wind profile by the least square method using the formula (5)*inflowSurface roughness Z0，loAnd the atmospheric stability L which is not interfered by the wind power fieldinflow。
And 7: stability correction function psi of wind speedmThe following steps can be taken:
Figure BDA0003320024250000101
wherein the content of the first and second substances,
Figure BDA0003320024250000102
ψ0＝-ln an+31/2bnan 1/3π/6 (8.2)
an＝0.33，bn＝0.41，am＝6.1，bm＝2.5 (8.3)
and 8: the wind power plant is assumed to have small influence on ground heat flux, and the ground temperature changes little. The surface heat flux can be calculated according to the similar theory of moxin-obufh:
Figure BDA0003320024250000103
wherein, thetasIs the ground temperature.
And step 9: for the part except the tail current layer of the fully developed wind power plant, according to the result of the large vortex simulation, (1)
Figure BDA0003320024250000104
Or zh+ D/4 < z ≦ δ) integrated (from z) according to equation (4)0Integration to z) gives:
Figure BDA0003320024250000105
Figure BDA0003320024250000111
wherein the content of the first and second substances,
Figure BDA0003320024250000112
respectively representing the atmospheric stability below the hub height and above the hub height inside the wind farm. δ represents the inner boundary layer height of the wind farm.
Step 10: equation (4) can be expressed as:
Figure BDA0003320024250000113
wherein the equivalent vortex viscosity coefficient vT＝κzu*。
Step 11: for fully developed wind farm wake horizon (z)h-3D/4≤z≤zh+ D/4), the phenomenon of speed reduction and turbulence kinetic energy enhancement occurs due to the action of the wind wheel. Equivalent vortex-viscous system vTNeed to be from vT＝κzu*Increase to vT＝(κzu*+vw). Equation (4) can be varied as:
Figure BDA0003320024250000114
wherein v iswIs the wake additional vortex viscosity coefficient.
Step 12: defining dimensionless wake additional vortex viscosity coefficient
Figure BDA0003320024250000115
The formula (13) can be modified:
Figure BDA0003320024250000116
step 13: the turbulence level increase of the wake layer is caused by the momentum loss of the wind wheel
Figure BDA0003320024250000117
Is proportional, thereby estimating a turbulent velocity scale of
Figure BDA0003320024250000118
The wake length scale is the wind wheel diameter D, so the wake additional vortex-viscosity coefficient can be estimated as:
Figure BDA0003320024250000119
step 14: the dimensionless wake additional vortex viscosity coefficient is obtained by equation (15):
Figure BDA0003320024250000121
step 15: to obtain a reaction with CftDirectly related values, let D be zhObtained from (4)
Figure BDA0003320024250000122
Figure BDA0003320024250000123
Wherein z ish100m, and z 0m 1 m. Equation (16) can be simplified as:
Figure BDA0003320024250000124
step 16: equation (14) can be obtained by performing an indefinite integration:
Figure BDA0003320024250000125
and step 17: will be expressed as the formula (18)
Figure BDA0003320024250000126
The integral constant C is determined in conjunction with equations (10), (11), respectively, and can be obtained:
Figure BDA0003320024250000127
Figure BDA0003320024250000128
step 18: according to the continuity, the formulae (19) and (20) are represented by the formula (z) ═ zhThe time-average wind speeds are equal, and the following can be obtained:
Figure BDA0003320024250000129
step 19: to simplify the formula, define w:
Figure BDA00033200242500001210
Figure BDA0003320024250000131
step 20: equation (21) can be simplified as:
Figure BDA0003320024250000132
step 21: bringing the formula (24) into the formula (1) to obtain wind field equivalent roughness analytical models under different atmospheric stability:
Figure BDA0003320024250000133
the independent unknown parameters contained in the calculation formula (25) of the equivalent roughness of the wind power plant comprise u*lo，u*hiAnd Z0，hi. The solving process is as follows:
obtaining an equation according to the equation (5) and the equation (11) that the boundary layer height δ wind speed is equal in the wind farm:
Figure BDA0003320024250000134
the value of the boundary layer height delta in the wind power plant is as follows: the stable atmospheric boundary layer (SBL) was taken to be 600 meters, the neutral atmospheric boundary layer (NBL) was taken to be 850 meters, and the unstable atmospheric boundary layer (CBL) was taken to be 1100 meters. And the height delta of the boundary layer in the wind power plant can also be measured by measurement means such as direct observation, ground-based remote sensing, space-based remote sensing and the like.
Three independent equations (1), (24) and (26) are combined, and u is calculated by an iterative method*lo，u*hiAnd Z0，hi。
Step 22: changing z to zhThe wind speed of the hub in the wind power plant can be fully developed by the drive-in type (20):
Figure BDA0003320024250000141
the effectiveness of the process of the invention is illustrated by the following specific examples.
The embodiment adopts the large eddy simulation data to calculate and verify the accuracy of the model. In all the examples, the hub height is equal to the rotor diameter, i.e. zhD100 m. The different atmospheric stability large vortex simulation example settings are shown in table 1. A periodic boundary condition is adopted to simulate an infinite wind power field, and different atmospheric stabilities are simulated by changing ground heat flux. And acquiring the equivalent roughness of the wind power plant according to the wind speed of 2 Zh.
Model calculation:
the method comprises the following steps: the inflow atmospheric stability parameter L is obtained according to the methodinflow；
Step two: according to the equations (6) and (7), the inflow friction speed u is calculated*inflowAnd ground roughness Z0，lo；
Step three: calculating the ground heat flux q according to equation (9)s；
For the convenience of verification, the parameters are adjusted to be consistent with the large vortex simulation, but the actual measurement data can still be adopted in practical application.
Step four: joint type (1), (24) and (26), and iterative computation u*lo、u*hiAnd Z0，hi。
The comparison of the wind power plant equivalent roughness model calculation result and the LES result and the relative error are shown in the figures 2 and 3.
Table 1 different atmospheric stability large eddy simulation example settings
Figure BDA0003320024250000142
As can be seen from the comparison results, the model calculation substantially coincided with the large vortex simulation results.
The top down model is developed by considering the atmospheric stability, and the analytic model for calculating the equivalent roughness of the fully developed wind power plant is deduced by considering the adverse effect of the wind power plant on the atmospheric stability.
The method has the innovation points that: according to the model, the correction of the atmospheric stability on the wind speed profile is considered on the basis of a top down boundary layer model, so that the influence of the atmospheric stability is introduced, and the equivalent roughness and the equivalent friction speed of the wind power plant are predicted. The method is closer to the running state of the actual wind power plant under the non-neutral atmosphere, and the expression is simple and easy to calculate.
So far, the technical solutions of the present invention have been described in connection with the preferred embodiments shown in the drawings, but it is easily understood by those skilled in the art that the scope of the present invention is obviously not limited to these specific embodiments. Equivalent changes or substitutions of related technical features can be made by those skilled in the art without departing from the principle of the invention, and the technical scheme after the changes or substitutions can fall into the protection scope of the invention.

Claims (10)

1. A fully-developed wind power plant equivalent roughness calculation method considering atmospheric stability is characterized by comprising the following steps:
acquiring parameters of a wind turbine generator;
obtaining an influence factor of interaction of atmospheric stability and a wind power plant on equivalent roughness;
acquiring a dimensionless wake flow additional vortex-viscosity coefficient of the wind power plant;
and calculating the equivalent roughness of the wind power plant according to the parameters of the wind turbine generator, the influence factors and the dimensionless wake flow additional vortex-viscosity coefficient of the wind power plant.
2. The fully developed wind farm equivalent roughness calculation method according to claim 1, wherein the step of calculating the equivalent roughness of the wind farm according to the parameters of the wind turbines, the impact factors and the dimensionless wake additional vortex viscosity coefficients specifically comprises calculating the equivalent roughness of the wind farm according to the following equation:
Figure FDA0003320024240000011
wherein D is the diameter of the wind wheel of the wind turbine generator, zhThe hub height of the wind turbine generator;
w(Lhi) Is the influence factor of the interaction of the atmospheric stability above the height of the hub and the wind power plant on the equivalent roughness, w (L)lo) The influence factor of the interaction of the atmospheric stability below the height of the hub and the wind power plant on the equivalent roughness is shown;
Figure FDA0003320024240000012
attaching vortex viscosity coefficients to dimensionless wake flows of the wind farm,
Figure FDA0003320024240000013
CTis the thrust coefficient of the wind wheel, SxAnd SyThe parameters are dimensionless parameters of the flow direction and the span-wise distance of the wind turbine generator relative to the diameter of the wind wheel.
3. The fully developed wind farm equivalent roughness calculation method according to claim 2, characterized in that the factor of influence w (L) of the interaction of the atmospheric stability above the hub height and the wind farm on the equivalent roughnesshi) Determined by the following equation:
Figure FDA0003320024240000014
wherein the content of the first and second substances,
Figure FDA0003320024240000015
and
Figure FDA0003320024240000016
respectively, are atmospheric stability correction functions at different heights above the hub height.
4. The fully developed wind farm equivalent roughness calculation method according to claim 2, characterized in that the influence factor w (L) of the interaction of the atmospheric stability below the hub height and the wind farm on the equivalent roughnesslo) By the followingThe equation determines:
Figure FDA0003320024240000021
wherein the content of the first and second substances,
Figure FDA0003320024240000022
and
Figure FDA0003320024240000023
respectively, the atmospheric stability correction functions at different heights below the hub height.
5. Method for fully developed wind farm equivalent roughness calculation according to any of the claims 2 to 4, characterized in that dimensionless wake additional vortex viscosity coefficients of the wind farm
Figure FDA0003320024240000024
Determined by the following equation:
Figure FDA0003320024240000025
wherein κ is the karman constant; u. of*Is the friction speed;
Figure FDA0003320024240000026
the wind speed is averaged over time and space for the hub altitude plane.
6. The fully developed wind farm equivalent roughness calculation method according to claim 1, wherein the wind farm is a fully developed area of a large wind farm.
7. A fully developed wind power plant generated power estimation method is characterized by comprising the following steps:
calculating an equivalent roughness for a fully developed wind farm according to the method of any one of claims 1 to 6;
determining the wind speed for fully developing the height of the hub of the wind power plant according to the calculated equivalent roughness;
and according to the wind speed of the hub height of the wind power plant, performing linear interpolation on the wind speed-power curve of the wind turbine generator to estimate the generating power of the wind power plant.
8. The fully developed wind farm generated power estimation method according to claim 7, wherein the wind speed for the fully developed wind farm hub height is determined by the following equation:
Figure FDA0003320024240000031
wherein z is0，hiTo an equivalent roughness, u*hiD is the diameter of the wind wheel of the wind turbine generator, z for fully developing the friction speed of the wind power plant above the hub heighthHub height of wind turbine, w (L)hi) Is an influence factor of the interaction of the atmospheric stability above the height of the hub and the wind power plant on the equivalent roughness,
Figure FDA0003320024240000032
vortex viscosity coefficients are added to the dimensionless wake of the wind farm, and κ is the karman constant.
9. A micro-site selection method for fully developing a wind power plant is characterized by comprising the following steps:
carrying out primary micro site selection on the wind power plant;
estimating the generation power of the wind turbines of a fully developed wind park according to the method of claim 7 or 8, based on the parameters of the preliminary micro-siting;
calculating the sum of the generated power of all the wind turbine generators, judging whether the sum of the generated power meets the design requirement, if so, selecting the site to pass, and if not, readjusting the site selection parameter until the sum of the generated power meets the design requirement.
10. The method for micro-siting a fully developed wind farm according to claim 9, characterized in that the step of "preliminary micro-siting a wind farm" comprises in particular: and (4) performing micro site selection on the wind power plant by adopting a Jensen wake flow model and wake flow square sum superposition.
CN202111242620.2A 2021-10-25 2021-10-25 Fully developed wind power plant equivalent roughness calculation method considering atmospheric stability Pending CN114154296A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111242620.2A CN114154296A (en) 2021-10-25 2021-10-25 Fully developed wind power plant equivalent roughness calculation method considering atmospheric stability

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111242620.2A CN114154296A (en) 2021-10-25 2021-10-25 Fully developed wind power plant equivalent roughness calculation method considering atmospheric stability

Publications (1)

Publication Number Publication Date
CN114154296A true CN114154296A (en) 2022-03-08

Family

ID=80458752

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111242620.2A Pending CN114154296A (en) 2021-10-25 2021-10-25 Fully developed wind power plant equivalent roughness calculation method considering atmospheric stability

Country Status (1)

Country Link
CN (1) CN114154296A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117291015A (en) * 2023-09-07 2023-12-26 宁夏回族自治区电力设计院有限公司 Equivalent roughness-based stock wind power plant modeling method and system

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117291015A (en) * 2023-09-07 2023-12-26 宁夏回族自治区电力设计院有限公司 Equivalent roughness-based stock wind power plant modeling method and system

Similar Documents

Publication Publication Date Title
Bartl et al. Blind test comparison of the performance and wake flow between two in-line wind turbines exposed to different turbulent inflow conditions
Wu et al. Large-eddy simulation of wind-turbine wakes: evaluation of turbine parametrisations
Allaerts et al. Sensitivity and feedback of wind-farm-induced gravity waves
Ren et al. Numerical simulation and prediction of spatial wind field under complex terrain
Edwards et al. PIV measurements and CFD simulation of the performance and flow physics and of a small‐scale vertical axis wind turbine
CN106548414B (en) Method for calculating power generation capacity of offshore wind farm
CN108717593A (en) A kind of microcosmic structure generated energy appraisal procedure based on wind wheel face equivalent wind speed
Apsley et al. Unsteady thrust on an oscillating wind turbine: Comparison of blade-element momentum theory with actuator-line CFD
CN112784509A (en) Wind power generator power calculation method based on improved Gaussian distribution wake flow model
CN110321632B (en) Method for calculating equivalent roughness of fully developed wind power plant
CN115943255A (en) System and method for measuring turbulence of wind flow in complex terrain with LiDAR
CN116341410B (en) Typical terrain wind field simulation method and system considering different stabilities
CN115062563B (en) Wind turbine three-dimensional wake flow wind speed calculation method based on space change
Cathelain et al. Calibration of a super-Gaussian wake model with a focus on near-wake characteristics
Smaïli et al. On the rotor effects upon nacelle anemometry for wind turbines
CN114154296A (en) Fully developed wind power plant equivalent roughness calculation method considering atmospheric stability
Odemark Wakes behind wind turbines-Studies on tip vortex evolution and stability
Wei et al. Power-generation enhancements and upstream flow properties of turbines in unsteady inflow conditions
CN112580891A (en) Method for predicting turbulent vertical momentum flux of boundary layer of wind power plant and evaluating generated energy
Polster Comprehensive comparison of analytical wind turbine wake models with wind tunnel measurements and wake model application on performance modelling of a downstream turbine
CN114139346A (en) Fully developed wind power plant equivalent roughness calculation method
Jackson Application of Reynolds stress model using direct modeling and actuator disk approaches for a small-scale wind turbine
Burmester et al. Determination of scaled wind turbine rotor characteristics from three dimensional RANS calculations
Polster et al. Experimental validation of analytical wake and downstream turbine performance modelling
Chivaee Large eddy simulation of turbulent flows in wind energy

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