2. 中国气象局公共气象服务中心, 北京 100081
2. Public Meteorological Service Center of CMA, Beijing 100081
近地层风速的垂直分布即风廓线是大气边界层研究和工程设计所需的基本气象参数[1-2],也是风电场风资源评估的关键要素。目前,近地层风廓线的获取主要依靠测风塔的梯度观测[3]。测风塔稳定性和可靠性高,但由于测风塔单点的观测代表性有限,在复杂地形上可能需要布设多个测风塔,导致成本增加[4]。因此,已有一些代替测风塔的仪器研发应用在边界层风观测中[5-8],脉冲激光风廓线仪就是其中一种。
脉冲激光风廓线仪是向大气发射激光脉冲,接收大气 (气溶胶粒子和大气分子) 的后向散射信号,通过分析发射激光的径向多普勒频移反演风速,是一种新型的测风手段。夏俊荣等[9]对北京325 m气象铁塔和激光风廓线仪WINDCUBE在63~200 m之间的6个层高度观测的10 min平均水平风速、风向进行对比,结果显示水平风速、风向相关系数分别达到0.98和0.99以上,表明该风廓线仪是一款观测结果可靠、机动性良好的低层大气风廓线测量仪器。Caadillas等[10]对海上测风平台FINO1上的脉冲激光风廓线仪WINDCUBE与该平台100 m塔上的超声风速仪、杯式风速仪的测风结果进行对比,结果显示10 min平均水平风速、风向线性拟合度均在0.99以上,但在进行数据筛选时将降雨期间的数据剔除,未分析降雨对脉冲激光风廓线仪观测的影响。Gottschall等[11]在丹麦国家实验室进行了脉冲激光风廓线仪WINDCUBE V2与杯式测风仪的对比试验,结果显示60~116 m高度层的回归斜率均小于1.5%,水平风速误差不受风向切变影响。郑启明等[12]比较激光风廓线仪与超音波风速计的风速风向发现,10 min平均风向有良好的一致性,当风速接近或大于10 m·s-1时,两种仪器的风速测量误差在10%以内,低风速的相对误差较大。
以上研究多基于对平均风速、风向的比较,对垂直气流、风的脉动特征量以及对降雨条件下激光风廓线仪的测风数据对比鲜见报道。而工程上也需要准确测量垂直速度、湍流强度、均方差等脉动特征量,同时还要考虑仪器在不同条件下的测风性能。为确定脉冲激光风廓线仪WINDCUBE V2的测风性能,针对风速、风向、湍流强度等风工程和风能资源研究关注的要素,与测风塔上杯式风速仪和超声风速仪观测结果进行了一系列的对比分析,也为确认这类遥感式测风仪的观测数据可靠性的检验方法提供参考。
1 对比观测试验介绍 1.1 试验场地描述对比试验的场地设在广东省湛江市东海岛气象观测站,满足《地面气象观测规范》[13]要求。站内设有高度为100 m的气象观测塔,脉冲激光风廓线仪WINDCUBE V2放置于观测塔东南面5 m处。
1.2 观测仪器设置及激光风廓线仪测风原理测风塔的50 m,70 m和100 m高度层安装了NRG-Symphonie型测风仪,风速量程范围为1~96 m·s-1,记录逐10 min平均风速、风向、3 s阵风风速和1 s相对10 min平均风速的标准差。测风塔65 m高度层安装了WindMaster Pro型超声波测风仪,采样频率是10 Hz,最大风速量程65 m·s-1,记录每0.1 s的u,v,w三维风速值、虚温和数据判别码。采用文献[14-15]的方法对超声测风数据进行质量控制处理。
脉冲激光风廓线仪WINDCUBE V2的设备性能参数见表 1。
|
|
表 1 脉冲激光风廓线仪参数表 Table 1 Specifications of WINDCUBE V2 Lidar |
WINDCUBE V2脉冲激光风廓线仪向空中发射激光脉冲,并接收返回信号,比较发出和返回信号频率的变化,利用多普勒频移原理算出指定高度气流径向速度。激光风廓线仪可向4个方向进行测量,扫描锥角度为28°,通过三维坐标系投影算出风的水平和垂直分量。激光风廓线仪可同时测量输出12个任意高度层的平均风、阵风和垂直速度等要素。
2 观测数据处理 2.1 对比观测数据的同步性调整为确保各种仪器的观测数据代表同一空气团的运动状况,需要对数据的同步性进行精准调整。由于NRG测风仪和脉冲激光风廓线仪的10 min数据记录所代表的观测时段有10 min的时间差①,即NRG风速仪在11:10(北京时,下同) 和脉冲激光风廓线仪在11:20记录的才是同一时间段的观测数据,因此,对数据进行10 min的时差调整。超声风速仪与脉冲激光风廓线仪也进行了相应的同步性调整。然后再利用转折天气过程中3种仪器测得的风速、风向的显著突变点进行一致性比对分析,确保3种仪器测量的观测数据达到完全同步。
① Leosphere公司.WINDCUBE V2 Lidar Remote sensor用户手册.62.
2.2 对比观测数据的有效性检验和样本筛选为保证对比观测数据的可靠性和一致性,对数据进行检验和筛选:① 依据数据有效率筛选样本。选取脉冲激光风廓线仪50 m,65 m,70 m和100 m高度层10 min数据有效率不低于80%的样本,及相应时刻的杯式测风仪和超声风速仪的数据。② 剔除环境干扰因子。由于测风塔分别位于杯式测风仪和脉冲激光风廓线仪的西面和西北面,对这两个方位的来风有阻挡作用,因此,将激光雷达受测风塔阻挡的方位附近±10°,即240°~274°和304°~332°方位角的数据剔除。最终筛选出9204个10 min样本,进行激光风廓线仪与杯式测风仪比对,2885个10 min样本,进行激光风廓线仪与超声风速仪对比。
3 特征参数计算方法以风能和风工程参数分析需求为主要目标,对风的主要统计特征参数进行对比。
① 拟合度。拟合度反映了拟合参数与实测值的吻合程度。
② 偏差。偏差是随机变量与标准值之差。文中用E表示激光风廓线仪与杯式测风仪测量值的偏差。由于风向在360°左右时会出现虚假突变,依据式 (1) 对脉冲激光风廓线仪的风向数据进行转换:
|
(1) |
式 (1) 中,θ1为杯式测风仪的风向值,θ2和θ′2分别为换算前后的脉冲激光风廓线仪的风向值。
③ 相对偏差。相对偏差在此表征两类仪器测风数据之偏差的离散度,表达式如下:
|
(2) |
式 (2) 中,x和y分别为杯式测风仪与脉冲激光风廓线仪的测量值。
④ 湍流强度。湍流强度反映了风的脉动强度,是确定结构脉动风荷载的关键参数,定义湍流强度为10 min时距的脉动风速标准差与水平平均风速的比值[16]。
⑤ 风廓线幂指数。风廓线用于描述风随高度的变化情况。工程应用中常以幂指数或对数率来量化拟合描述风廓线,目前,国内外大部分相关规范推荐幂指数形式[17-20]。
4 测风结果对比 4.1 对比观测数据说明杯式测风仪与脉冲激光风廓线仪的同步观测日期为2012年5月21日-8月16日。超声风速仪与脉冲激光风廓线仪的同步观测日期为2012年7月6-29日。对比观测期间,最大风速出现在5月21日,100 m高度层测得10 min平均最大风速为15.8 m·s-1,风向为83°;最大降水强度出现在8月5日,10 min降水量达15.5 mm,100 m高度层的风速为11.9 m·s-1,风向为225°。
本次试验中,激光风廓线仪先后对40~300 m范围内的16个高度进行观测,图 1给出激光风廓线仪测得的总体样本和降雨样本的数据有效率随高度的变化,总体样本有效率在40~155 m高度层大于97%,155 m以上随高度增加,样本有效率减小,275 m高度处为56.2%。降雨情况下数据有效率下降趋势略大一些。
|
|
| 图1 激光风廓线仪数据有效率随测量高度的变化 Fig.1 Variation of data availability of WINDCUBE V2 Lidar with measuring height | |
4.2 总体样本对比
基于式 (3),对激光风廓线仪和杯式测风仪测量的风速、风向、湍流强度、3 s阵风风速和风速标准差进行线性拟合, 并在表 2中给出拟合参数及偏差。图 2给出100 m高度层的风速、风向散点图。
|
|
| 图2 100 m高度层总体样本的风参数散点图及拟合线 (a) 风速, (b) 风向 Fig.2 Scatter plot of wind parameter and the fitting curve for all samples at 100 m height (a) wind speed, (b) wind direction | |
|
(3) |
|
|
表 2 脉冲激光风廓线仪与杯式测风仪总体样本的风参数拟合参数及偏差 Table 2 Results of regression analysis for wind parameters measured by lider and cup anemometer |
式 (3) 中,x和y分别为杯式测风仪与脉冲激光风廓线仪的测量值,a和b分别为拟合直线的斜率和截距。下文拟合均基于式 (3)。
图 2显示,两种仪器测量的风参数具有很好的一致性。由表 2可见,风参数的拟合关系式的斜率均接近1,风速、风向的拟合度均在0.99以上。脉冲激光风廓线仪的水平风速较杯式测风仪的测量值略小0.041~0.088 m·s-1。两种仪器测得的风向具有一定的固定差,且50 m和70 m,100 m显示的风向偏差值差异较大,原因可能是由于测风塔上的测风仪和激光风廓线仪在人工安装的过程中,使用指南针进行方向调整时存在一定的人为偏差。湍流强度的拟合度大于0.67,脉冲激光风廓线仪的湍流强度偏大约0.01。3 s阵风风速的拟合度大于0.96,脉冲激光风廓线仪的阵风风速偏大约0.3 m·s-1。风速标准差的拟合度大于0.79,随着高度的增加,拟合度有降低的趋势,同时偏差增大。
4.3 不同水平风速下的对比脉冲激光风廓线仪的风速测量范围为0~80 m·s-1,但否在各个风速量程均具有较好的探测性能需要进一步讨论,以确定该仪器在强风天气下的适用性。依据不同风速条件下边界层风特性参数的不同特征表现,将测风数据分为大风速和小风速两类。以杯式测风仪10 min平均风速为基础,将1~6 m·s-1的数据定为小风样本,将不低于6 m·s-1的数据定为大风样本数据进行比对分析。
图 3给出激光风廓线仪和杯式测风仪测量的大小风样本风参数的散点图及拟合线,拟合参数及偏差在表 3中给出。结合图 3和表 3可见,大风样本的风速相关较高,水平风速拟合度均在0.98以上,比小风样本的风速拟合度略大0.004~0.015。风向的拟合度均在0.99以上,大风样本的风向相关较高,且风向偏差稍小。小风和大风样本的湍流强度拟合度分别为0.60~0.75和0.76~0.84。大风样本的阵风风速拟合度大于0.93,较小风样本的拟合度偏大0.033~0.102。大风样本的风速标准差的相关较高,拟合度均在0.84以上,较小风样本的拟合度偏大0.1左右。总体而言,两种仪器测得的大风样本的风参数相关更好。
|
|
| 图3 100 m高度层大小风样本的风参数散点图及拟合线 (红点和实线:小风样本及拟合线;蓝点和虚线:大风样本及拟合线) (a) 风速, (b) 风向, (c) 湍流强度, (d) 阵风风速, (e) 风速标准差 Fig.3 Scatter plot of wind parameters and the fitting curve of samples of weak wind and strong wind observed at 100 m height (red dot and solid line denote weak wind and the linear fitting line; blue dot and dashed line denote strong wind and linear fitting line) (a) wind speed, (b) wind direction, (c) turbulence intensity, (d) gust wind speed, (e) wind speed standard deviation | |
|
|
表 3 脉冲激光风廓线仪与杯式测风仪测得的大小风样本的风参数拟合参数及偏差 Table 3 Results of regression analysis for wind parameters measured by lider and cup anemometer for samples of weak wind and strong wind |
基于50 m,70 m和100 m的测风数据,计算两种仪器测得的大风 (50 m水平风速不低于6 m·s-1) 和小风 (50 m水平风速在1~6 m·s-1) 样本的平均风廓线,并根据幂指数律进行拟合,结果显示大风和小风情况下两种仪器测得的平均风廓线幂指数值均一致。
图 4给出两种仪器观测的10 min 100 m高度平均风速相对偏差的概率分布。其中,概率分布的统计区间间隔均为1,[-0.5, 0.5]为1个区间,图 4中标示的[-6.5%, 3.5%], p=90.5%表示有90.5%的相对偏差集中在-6.5%到3.5%的范围内 (下同)。由图 4可见,大风样本水平风速相对偏差的波动范围显著较小。综合50 m,70 m和100 m 3个高度层的样本来看,90%以上的风速相对偏差值在-4.5%到2.5%之间,而90%以上小风样本的水平风速相对偏差在-11.5%到5.5%之间,大风情况下脉冲激光风廓线仪对平均风速的测量效果更好。
|
|
| 图4 100 m高度层小风 (a) 和大风 (b) 样本风速相对偏差的概率分布 Fig.4 Probability distribution of the relative deviation of wind speed at 100 m height for samples of weak wind (a) and strong wind (b) | |
另外,两种仪器观测的大风样本的风向偏差更集中在较小范围内。90%以上的大风样本的风向偏差为8.5°~14.5°,而90%以上小风样本的风向偏差集中在5.5°~18.5°。
4.4 降雨天气下的对比有降雨发生时,由于脉冲激光风廓线仪会同时接收到大气湍流回波和降水粒子的散射回波,降水粒子的垂直下落速度往往会使信号谱发生明显变化,从而造成测量误差[21-23],因此,降水往往是影响风廓线仪测风精度的重要因素之一。为了了解WINDCUBE V2在降雨情况下的观测效果,本文根据试验基地内国家气象站 (区站号为59751) 的降雨资料,选出观测期间有降雨的10 min测风样本350个、无降雨的10 min测风样本8545个进行对比分析。
图 5给出脉冲激光风廓线仪和杯式测风仪测量的有无降雨样本的风参数散点图及其拟合线,拟合参数及偏差在表 4中给出。由图 5和表 4可见,在10 min降水量小于15 mm的降雨情况下,降雨对脉冲激光风廓线仪的平均风速、风向的测量没有显著影响,线性拟合的斜率均接近于1,降雨样本的拟合度要略低于无降雨样本,但均大于0.99。降雨情况下,脉冲激光风廓线仪的偏差波动较大。降雨对脉冲激光风廓线仪湍流强度的观测没有影响,降雨情况下50 m高度层湍流强度的拟合度相对更高。降雨样本阵风风速的拟合度要较无降雨样本的拟合度略减小约3%~6%。降雨样本的风速标准差拟合度要较无降雨样本的拟合度显著减少10%以上。降雨情况下风速标准差的偏差及其波动也较大。可见,降雨对脉冲激光风廓线仪风速标准差的测量有一定影响。
|
|
| 图5 100 m高度层有无降雨样本风参数散点图及拟合线 (红点和实线:降雨样本及拟合线;蓝点和虚线:无降雨样本及拟合线) (a) 风速, (b) 风向, (c) 湍流强度, (d) 阵风风速, (e) 风速标准差 Fig.5 Scatter plot of wind parameters and the fitting curve of samples with and without rain at 100 m (red dot and solid line:samples with rain and the linear fitting line; blue dot and dashed line: samples without rain and linear fitting line) (a) wind speed, (b) wind direction, (c) turbulence intensity, (d) gust wind speed, (e) wind speed standard deviation | |
|
|
表 4 脉冲激光风廓线仪与杯式测风仪测得的有无降雨样本的风参数拟合参数及偏差 Table 4 Results of regression analysis for wind parameters measured by lider and cup anemometer for samples with and without rain |
图 6给出两种仪器测得的风速相对偏差的概率分布及其随雨强的变化。由图 6可见,无降雨样本的风速相对偏差更为集中在较小的范围内。另外,随观测高度的增加,风速相对偏差也越小。当降雨增强时风速的相对偏差也有减小趋势,对降雨时段的测风样本分析发现,降雨增强往往伴随的风速也增大。
|
|
| 图6 100 m高度层无降雨 (a) 及有降雨样本 (b) 风速相对偏差的概率分布及其随降水量变化 (c) Fig.6 Probability distribution of the relative deviation of wind speed for samples without rain (a), withrain (b) and variations of the relative deviations with rainfall intensity (c) at 100 m height | |
两种仪器测得的90%以上的风向偏差集中在6.5°~15.5°范围内,与无降雨样本的分布差异很小。可见10 min降水量小于15 mm的降雨对脉冲激光风廓线仪的风向测量几乎没有影响,这一现象仍可能与强降雨往往伴随较大的风速有关。另外,当降水增强时,风向偏差有减小的趋势。
4.5 水平风速和降雨对垂直速度测量的影响气流的垂直速度是大气风场结构研究的重要参数,潘乃先等[24]研究发现,横向风对多普勒天气雷达的垂直速度测量会造成误差。另外,有降雨时,雨滴的散射也是主要散射机制,可能对激光有较强的衰减作用,因此,本节基于脉冲激光风廓线仪与超声风速仪的对比观测资料,分析水平风速大小和降雨对脉冲激光风廓线仪WINDCUBE V2的垂直速度测量的影响。
本文的垂直速度,正值表示上升气流,负值表示下沉气流。图 7给出无降雨时,两种仪器测得的垂直速度根据式 (3) 进行线性拟合得到的拟合度随水平风速范围的变化。由图 7可见,当水平风速范围增加时,拟合度有显著降低的趋势,当水平风速为0~1.5 m·s-1时,拟合度为0.872,而当水平风速范围增大为0~6 m·s-1时,拟合度下降为0.25。可见,水平风速对脉冲激光风廓线仪的垂直速度测量有显著影响。
|
|
| 图7 脉冲激光风廓线仪与超声风速仪测得的无降雨样本的垂直速度的拟合度随水平风速的变化 Fig.7 Variation of vertical airflow velocity fitting degree of lidar against ultrasonic anemometer with the horizontal wind velocity for samples without rain | |
降雨情况下,当雨强增加时,激光风廓线仪测得的垂直速度显著小于超声风速仪的测量值,显示较大的下沉气流,当10 min降水量在2 mm以上时,垂直速度的偏差基本集中在-3.2~-1.4 m·s-1。
为了进一步细化分析降雨的影响,计算两种仪器的1 min平均水平风速和垂直速度,并结合气象站的5 min降水量计算的1 min降水量 (假设5 min内降水量平均分布),给出观测期间测得的2个较强降雨过程的1 min水平风速、垂直速度和1 min降水量的时程变化图 (图 8)。由图 8可见,在降雨前后,两种仪器的垂直速度均具有较好的一致性,但当降雨出现时,激光风廓线仪测得的垂直速度显著减小,降水越强,垂直速度下降越明显,当1 min降水量达到0.5 mm时,激光风廓线仪的垂直速度较超声风速仪的测量值最大下降4 m·s-1左右,而此时对应的水平风速则受降雨影响不大。降雨时,风廓线仪测得的下沉垂直速度是大气垂直运动和降水粒子下降速度的合成,激光风廓线仪对垂直速度的测量效果欠佳。
|
|
| 图8 1 min垂直速度、水平风速和降水量时程 Fig.8 Time series of 1 min vertical airflow and horizontal velocity and the rainfall | |
5 结论
研究表明:
1) 脉冲激光风廓线仪与杯式测风仪测量的10 min平均风速、风向的拟合度均大于0.99,3 s阵风风速的拟合度大于0.96,相关性高。100 m高度的平均风速、风向的偏差平均值分别为-0.088 m·s-1和11.7°,激光风廓线仪测得的3 s阵风风速要较杯式测风仪的测量值略大0.347~0.378 m·s-1。
2) 脉冲激光风廓线仪与杯式测风仪测得的湍流强度的拟合度为0.67~0.78,风速标准差的拟合度为0.79~0.87。
3) 大风情况下脉冲激光风廓线仪的测量效果更佳。在10 min降水量小于15 mm的情况下,降雨对激光风廓线仪的风速、风向、湍流强度和3 s阵风风速的测量没有显著影响,但降雨对风速标准差有一定影响。
4) 当水平风速增大和有降雨时,脉冲激光风廓线仪对垂直速度的测量效果欠佳。降雨时,受降雨粒子散射回波的影响,激光风廓线仪测得显著向下的垂直速度。
| [1] | Song Lili, Chen Wenchao, Wang Binglan, et al. Characteristics of wind profiles in the landing typhoon boundary layer. J Wind Eng Ind Aerodyn, 2016, 149: 77–88. DOI:10.1016/j.jweia.2015.11.008 |
| [2] | 中华人民共和国住房和城乡建设部. 建设结构荷载规范 (GB 50009-2012). 北京: 中国建筑工业出版社, 2012. |
| [3] | 张容焱, 张秀芝, 杨校生, 等. 台风莫拉克 (0908) 影响期间近地层风特性. 应用气象学报, 2012, 23, (2): 184–194. DOI:10.11898/1001-7313.20120207 |
| [4] | 王青梅, 张以谟. 气象激光雷达的发展现状. 气象科技, 2006, 34, (3): 246–249. |
| [5] | 陈元昭, 俞小鼎, 陈训来, 等. 2015年5月华南一次龙卷过程观测分析. 应用气象学报, 2016, 27, (3): 334–341. DOI:10.11898/1001-7313.20160308 |
| [6] | 王志春, 植石群, 丁凌云. 强台风纳沙 (1117) 近地层风特性观测分析. 应用气象学报, 2013, 24, (5): 595–605. DOI:10.11898/1001-7313.20130509 |
| [7] | 魏应植, 汤达章, 许健民, 等. 多普勒雷达探测"艾利"台风风场不对称结构. 应用气象学报, 2007, 18, (3): 285–294. |
| [8] | 李明华, 范绍佳, 王宝民, 等. 珠江三角洲秋季大气边界层温度和风廓线观测研究. 应用气象学报, 2008, 19, (1): 53–60. DOI:10.11898/1001-7313.20080110 |
| [9] | 夏俊荣, 王普才, 闵敏. 新型多普勒测风激光雷达Windcube的风参数观测与验证. 气候与环境研究, 2011, 6, (6): 733–741. |
| [10] | Canadillas B, Westerhellweg A, Neumann T, et al.Testing the Performance of a Ground-based Wind LiDAR System.DEWI MAGAZIN, 2011. |
| [11] | Gottschall J L, Michael S C.WINDCUBE V2 WLS7-0091 Lidar Verification Test at Høvsøre Test Site.Risø-I-3080(EN), 2010. |
| [12] | 郑启明, 罗元隆, 王军翰. 自然风场特性之实场量测与风洞模拟//第十六届全国结构风工程会议论文集. 成都: 西南交通大学出版社, 2013: 461-468. |
| [13] | 中国气象局. 地面气象观测规范. 北京: 气象出版社, 2003. |
| [14] | 宋丽莉, 陈雯超, 黄浩辉. 工程抗台风研究中风观测数据的可靠性和代表性判别. 气象科技进展, 2011, 1, (1): 33–37. |
| [15] | Song L L, Li Q S, Chen W C, et al. Wind characteristics of a strong typhoon in marine surface boundary layer. Wind and Structures, 2012, 15, (1): 1–15. DOI:10.12989/was.2012.15.1.001 |
| [16] | 盛裴轩, 毛节泰, 李建国, 等. 大气物理学. 北京: 北京大学出版社, 2003. |
| [17] | 项海帆, 林志兴, 鲍卫刚, 等. 公路桥梁抗风设计指南. 北京: 人民交通出版社, 1996. |
| [18] | Architectural Institute of Japan.Recommendations for Loads on Buildings (AIJ-RLB-2004).Tokyo, 2004. |
| [19] | ASCE.Minimum Design Loads for Buildings and Other Structures, Reston, VA, ASCE/SEI 7-05, New York, 2006. |
| [20] | Buildings Department.Code of Practice for Wind Effects in Hong Kong, Buildings Department of Hong Kong, 2004. |
| [21] | 邓闯, 阮征, 魏鸣, 等. 风廓线雷达测风精度评估. 应用气象学报, 2012, 23, (5): 523–533. DOI:10.11898/1001-7313.20120502 |
| [22] | 林晓萌, 何平, 黄兴友. 一种抑制降水对风廓线雷达水平风干扰的方法. 应用气象学报, 2015, 26, (1): 66–75. |
| [23] | 何平, 朱小燕, 阮征. 风廓线雷达探测降水过程的初步研究. 应用气象学报, 2009, 20, (4): 465–470. DOI:10.11898/1001-7313.200904011 |
| [24] | 潘乃先, 郑毅. 风对多普勒声雷达测风的某些影响. 北京大学学报 (自然科学版), 1986, 22, (1): 98–105. |
2017, 28 (3): 327-339



