鄱阳湖区风电场风电功率预报研究

风力发电作为一种无污染可再生的能源,已逐渐成为许多国家能源战略可持续发展的重要组成部分。风电场风能预报是风力发电开发中的关键技术问题。为研究鄱阳湖区风力发电预报技术,采用中尺度模式WRF和微尺度模块CALMET对鄱阳湖区长岭风电场进行了200 m水平分辨率风能预报,并根据长岭机组理论功率曲线表和实测数据拟合出理论和实际发电机组功率曲线模型及平均有功功率与发电量模型。根据WRF+CALMET模式预报风速及建立的发电机组功率曲线模型和平均有功功率与发电量模型,预报了长岭风电场发电量。结果表明:长岭风电场23座风机逐小时风速预报值与观测值相关系数为0.42~0.61,均方根误差为2.59~3.68,相对误差为-13.7%~17.4%;对整个风场,预报风速与观测风速的相关系数为0.55,均方根误差为2.8,相对误差为-4.79%。实测发电量值高于预报值,平均偏大39.7 kW,相对误差为-12.6%,预报值与实测值相关性较好,相关系数达到0.52。总体来说,根据中尺度数值模式预报的风速结合风功率、发电量模型预测出的发电量与实测值较为接近,但各月差异性较大。     

Simulation of Turbulent Flow Inside and Above Wind Farms: Model Validation and Layout Effects

A recently-developed large-eddy simulation framework is validated and used to investigate turbulent flow within and above wind farms under neutral conditions. Two different layouts are considered, consisting of thirty wind turbines occupying the same total area and arranged in aligned and staggered configurations, respectively. The subgrid-scale (SGS) turbulent stress is parametrized using a tuning-free Lagrangian scale-dependent dynamic SGS model. The turbine-induced forces are modelled using two types of actuator-disk models: (a) the ‘standard’ actuator-disk model (ADM-NR), which calculates only the thrust force based on one-dimensional momentum theory and distributes it uniformly over the rotor area; and (b) the actuator-disk model with rotation (ADM-R), which uses blade-element momentum theory to calculate the lift and drag forces (that produce both thrust and rotation), and distributes them over the rotor disk based on the local blade and flow characteristics. Validation is performed by comparing simulation results with turbulence measurements collected with hot-wire anemometry inside and above an aligned model wind farm placed in a boundary-layer wind tunnel. In general, the ADM-R model yields improved predictions compared with the ADM-NR in the wakes of all the wind turbines, where including turbine-induced flow rotation and accounting for the non-uniformity of the turbine-induced forces in the ADM-R appear to be important. Another advantage of the ADM-R model is that, unlike the ADM-NR, it does not require a priori specification of the thrust coefficient (which varies within a wind farm). Finally, comparison of simulations of flow through both aligned and staggered wind farms shows important effects of farm layout on the flow structure and wind-turbine performance. For the limited-size wind farms considered in this study, the lateral interaction between cumulated wakes is stronger in the staggered case, which results in a farm wake that is more homogeneous in the spanwise direction, thus resembling more an internal boundary layer. Inside the staggered farm, the relatively longer separation between consecutive downwind turbines allows the wakes to recover more, exposing the turbines to higher local wind speeds (leading to higher turbine efficiency) and lower turbulence intensity levels (leading to lower fatigue loads), compared with the aligned farm. Above the wind farms, the area-averaged velocity profile is found to be logarithmic, with an effective wind-farm aerodynamic roughness that is larger for the staggered case.     

大气边界层剪切湍流统计特性的风洞实验及其层次相似律

对风洞模拟的大气边界层的剪切湍流场的速度脉动时间序列信号以及某一高度位置上的小波变换分解重构的信号进行了系统的统计分析.结果表明, 在非均匀下垫面大气边界层的不同垂直高度上的信号, She-Leveque湍流层次相似律都成立, 层次相似指数β和最高激发态的奇异标度指数γ随垂直高度的变化而变化, 并给出不同尺度的分解重构信号的一些层次结构特点.     

用于风电场功率预测的逐时风速预报

为满足风电场精细化风速预报的需要,利用MM5模式的格点输出数值预报产品以及福建沿海地区两座海上测风塔的逐时气象资料,采用逐步回归的统计方法,将MM5数值预报产品与测风塔实测气象资料要素指标值共同引入回归方程,拟合两座测风塔特定高度的逐小时风速预报方程,详细介绍了该方法的实现过程.通过一年的预报效果指标以及15日的独立样本效果检验,结果表明该方法对于逐小时风速的预报有一定的应用价值.     

The Transport Time pdf for a Simple Turbulent Flow

In a turbulent fluid, the time a particle needs to travel from a point source to the observation point, can be considered as a random variable. It is shown that the probability density function (pdf) for this random variable is determined by the Lagrangian particle position pdf. The characteristics of the transport time pdf are discussed for the simple case of a turbulent fluid moving with a constant mean velocity.     

Large-Eddy Simulation of Turbulent Flow Over a Rough Surface

A family of wall models is proposed that exhibits moresatisfactory performance than previousmodels for the large-eddy simulation (LES) of the turbulentboundary layer over a rough surface.The time and horizontally averaged statistics such asmean vertical profiles of windvelocity, Reynolds stress, turbulent intensities, turbulentkinetic energy and alsospectra are compared with wind-tunnel experimental data.The purpose of the present study is to obtain simulatedturbulent flows that are comparable with wind-tunnelmeasurements for use as the wind environment for thenumerical prediction by LES of source dispersion in theneutral atmospheric boundary layer.     

垂直网格计算频散性的研究

在线性斜压原始方程组的基础上, 从频率和群速方面讨论了几种垂直跳点网格和时间—垂直跳点网格的计算频散性, 并指出了出现错误群速的垂直尺度范围。以便为原始方程大气模式选取垂直网格提供指导。     

Wind Tunnel And Field Measurements Of Turbulent Flow In Forests. Part I: Uniformly Thinned Stands

Many forest management methods alterstand density uniformly. The effectsof such a change on the wind andturbulence regimes in the forest arecritical to a number of processes governingthe stability of the stand and itsmicroclimate. We measured wind speed andturbulence statistics with a Dantec tri-axialhot-film probe in model forests of variousdensities (31–333 trees m^-2), created byremoving whole trees in a regular pattern in awind tunnel, and compared them with similarmeasurements made with propeller anemometers insimilarly thinned plots (156–625 trees ha^-1)within a Sitka spruce stand in Scotland. The results agree well, in general, with measurements made inother such studies with diverse canopy types.The systematic variations with density and verticalleaf-area distribution (which differed betweenwind-tunnel and field trees) in our work can explainmuch of the variability shown in scaled profiles ofbasic turbulence statistics reported in theliterature. The wind tunnel and field results are shown to be in good agreement overalldespite the difference in vertical leaf-areadistribution. Within-canopy and isolated-treedrag coefficients in the wind tunnel showthat tree-scale shelter effects increase astree density increases. The measurements indicatethat turbulence in the canopy is dominated bylarge-scale structures with dimensions of the sameorder as the height of the canopy as found inother studies but suggest that inter-tree spacing also modulates the size of these structures. These structures are associated with the sweeps that dominatemomentum exchange in the canopy and it is thisfact that allows the tri-axial probe to operate sowell despite the relatively narrow range of anglesin which the wind vector is correctly measured. Theratio of streamwise periodicity of these structuresto vorticity thickness varies systematically withtree density in the range 2.7–5.1, which spans theexpected range of 3.5–5 found in a laboratorymixing-layer, suggesting that tree spacing imposes another relevant length scale. This test andothers show that the results are in agreement withthe idea that canopy turbulence resembles that of a mixing layer even though they disagree with, and challenge the linear relationship between, streamwise periodicity andshear length scale presented recently in theliterature. The measurements are also in goodoverall agreement with simple drag models presented recently by other researchers.     

Experimental Study on the Wake Meandering Within a Scale Model Wind Farm Subject to a Wind-Tunnel Flow Simulating an Atmospheric Boundary Layer

The phenomenon of meandering of the wind-turbine wake comprises the motion of the wake as a whole in both horizontal and vertical directions as it is advected downstream. The oscillatory motion of the wake is a crucial factor in wind farms, because it increases the fatigue loads, and, in particular, the yaw loads on downstream turbines. To address this phenomenon, experimental investigations are carried out in a wind-tunnel flow simulating an atmospheric boundary layer with the Coriolis effect neglected. A \(3 \times 3\) scaled wind farm composed of three-bladed rotating wind-turbine models is subject to a neutral boundary layer over a slightly-rough surface, i.e. corresponding to offshore conditions. Particle-image-velocimetry measurements are performed in a horizontal plane at hub height in the wakes of the three wind turbines occupying the wind-farm centreline. These measurements allow determination of the wake centrelines, with spectral analysis indicating the characteristic wavelength of the wake-meandering phenomenon. In addition, measurements with hot-wire anemometry are performed along a vertical line in the wakes of the same wind turbines, with both techniques revealing the presence of wake meandering behind all three turbines. The spectral analysis performed with the spatial and temporal signals obtained from these two measurement techniques indicates a Strouhal number of \(\approx 0.20 - 0.22\) based on the characteristic wake-meandering frequency, the rotor diameter and the flow speed at hub height.     

Mean and Turbulent Flow Statistics in a Trellised Agricultural Canopy

Flow physics is investigated in a two-dimensional trellised agricultural canopy to examine that architecture’s unique signature on turbulent transport. Analysis of meteorological data from an Oregon vineyard demonstrates that the canopy strongly influences the flow by channelling the mean flow into the vine-row direction regardless of the above-canopy wind direction. Additionally, other flow statistics in the canopy sub-layer show a dependance on the difference between the above-canopy wind direction and the vine-row direction. This includes an increase in the canopy displacement height and a decrease in the canopy-top shear length scale as the above-canopy flow rotates from row-parallel towards row-orthogonal. Distinct wind-direction-based variations are also observed in the components of the stress tensor, turbulent kinetic energy budget, and the energy spectra. Although spectral results suggest that sonic anemometry is insufficient for resolving all of the important scales of motion within the canopy, the energy spectra peaks still exhibit dependencies on the canopy and the wind direction. These variations demonstrate that the trellised-canopy’s effect on the flow during periods when the flow is row-aligned is similar to that seen by sparse canopies, and during periods when the flow is row-orthogonal, the effect is similar to that seen by dense canopies.     

不同强风样本湍流特性参数的计算分析

在结构风工程中, 风湍流统计参数计算的正确与否直接影响到风荷载的计算精度。在实际风参数计算与分析中, 多选用风速较大的样本资料, 但过分强调大风可能产生不合理的计算结果。利用超声风速仪瞬时风速观测资料, 分别划分成相对强风和持续强风样本, 计算并比较其湍流统计特性参数, 发现湍流统计参数特性值 (湍流度、阵风因子、摩擦速度等) 与风速大小并不能很好匹配, 有时风速不大但其湍流特性值却很大, 反之也然。研究表明:选取的湍流风资料样本或统计方法不同, 都会影响风特性参数的计算结果, 进而影响到风荷载计算的精度。这一结果对于提高结构风工程中风参数计算与设计的科学性和合理性具有现实意义。     

几种垂直跳点网格计算频散性的比较研究

为证明两种新的垂直跳点网格--LZ网格和LY网格的有效性,从频率、群速和出现错误群速的垂直尺度范围等方面与现有的几种垂直跳点网格和时间--垂直跳点网格进行了比较.结果表明:LZ网格计算频散性与目前使用的最广泛的并且计算频散性也是最佳的CP网格相当,而LY网格次之,从而为模式设计者选取垂直网格提供了又一种选择.     

山东即墨风电场风力资源分析

以即墨气象站风速风向为参考,利用其30年风速资料确定风电场代表年,根据铁塔气象观测资料分析即墨风电场的风力资源状况;依据“风能资源评价技术规定”对即墨风电场的风速风向进行了检验修正。分析了测风塔各高度风速日、年变化,风向频率年月变化。计算分析了测风塔不同高度平均风速、风功率密度、有效小时数等风能参数。计算了不同风电机型60~70 m轮毂高度年发电量及其随高度变化规律;估算了风电场年发电量。得出即墨风电场达到3级风场标准,具有开发价值。     

大规模风电场建成后对风能资源影响的研究

考虑了千万千瓦级风电基地建成后风电机群对近地面层风速的影响,采用Frandsen研究了大规模风电场内部风速损失时所使用的方法,在内边界层已经充分发展成新边界层的区域内,对轮毂高度65m处风速Uh进行了计算。结果表明,风电场建成后研究区内,风速Uh与未建场时的65m风速U0相比变小,存在风速损失,该风速损失随着U0的增大而减小,与风电机的推力系数CT性质有关;大规模风电场建成后,Uh在3～20m.s-1范围内的平均风功率密度与未建场时U0在此范围内的平均风功率密度相比损失约为58.45%,这与建场地区建场前65m处风速值大小以及各风速值出现的概率有关。     

基于风速高相关分区的风电场风速预报订正

为了提高风电场风速预报和功率预测的精度和准确率,并考虑风机测风数据的不稳定因素,以多年服务的内蒙古中部某风力发电场A为研究区,在勘察风电场地形及风机布局后,按照季节、风向进行风机间风速时空相关性分析,划分出风机轮毂高度风速高相关为典型特征的风机网格分类片区,采用卡尔曼滤波方法,通过直接和间接两种订正方案,分别进行风机片区风速订正。结果表明:风速高相关风机片区的划分,对于提高风电场风速预报及功率预测精度和准确率具有一定作用,利用风电场区测风塔梯度观测风速,对风机片区进行间接订正,可有效改善数值模式预报风速,15个片区类型下相关系数由0.18~0.72提高至0.67~0.91,误差绝对值由1.6~2.9 m·s-1降低至1.0~1.5 m·s-1。     

Large-eddy Simulation of Turbulent Flow Across a Forest Edge. Part I: Flow Statistics

The statistics of turbulent flow across a forest edge have been examined using large-eddy simulation, and results compared with field and wind-tunnel observations. The moorland-to-forest transition is characterized by flow deceleration in the streamwise direction, upward distortion of the mean flow, formation of a high pressure zone immediately in front of the edge, suppression of the standard deviations and covariance of velocity components, and enhancement of velocity skewnesses. For the selected forest density, it is observed that the maximum distortion angle is about 8 degrees from the horizontal. Instead of approaching a downwind equilibrium state in a monotonic manner, turbulence (standard deviations and covariances of velocity components) and mean streamwise velocity undershoot in the transition zone behind the edge. Evolution of flow statistics clearly reveals the growth of an internal boundary layer, and the establishment of an equilibrium layer downwind of the edge. It is evident that lower-order moments generally adjust more quickly over the new rough surface than do higher-order moments. We also show that the streamwise velocity standard deviation at canopy height starts its recovery over the rough surface sooner than does the vertical velocity standard deviation, but completes full adjustment later than the latter. Despite the limited domain size upstream of the edge, large-eddy simulation has successfully reproduced turbulent statistics in good agreement with field and wind-tunnel measurements.     

Direct Numerical Simulation of a Turbulent Channel Flow with Forchheimer Drag

Boundary-Layer Meteorology - We characterize the turbulent flow, using direct numerical simulations (DNS), within a closed channel between two parallel walls with a canopy of constant areal density...     

A Boundary-Layer Scaling for Turbulent Katabatic Flow

Scaling relationships are proposed for the turbulent katabatic flow of a stably stratified fluid down a homogeneously cooled planar slope—the turbulent analogue of a Prandtl-type slope flow. The \(\Pi \) Theorem predicts that such flows are controlled by three non-dimensional parameters: the slope angle, the Prandtl number, and a Reynolds number defined in terms of the surface thermal forcing (surface buoyancy or surface buoyancy flux), Brunt-Väisälä frequency, slope angle, and molecular viscosity and diffusivity coefficients. However, by exploiting the structure of the governing differential equations in a boundary-layer form, scaled equations are deduced that involve only two non-dimensional parameters: the Prandtl number and a modified Reynolds number. In the proposed scaling framework, the slope angle does not appear as an independent governing parameter, but merely acts as a stretching factor in the scales for the dependent and independent variables, and appears in the Reynolds number. Based on the boundary-layer analysis, we hypothesize that the full katabatic-flow problem is largely controlled by two rather than three parameters. Preliminary tests of the scaling hypothesis using data from direct numerical simulations provide encouraging results.