1. 引言

在当前保护环境、应对气候变暖的迫切需求下，清洁、可再生的风能对环境友好，亟待进一步开发利用 [1] 。我国拥有超过1.8万公里的海岸线，海上风电潜力巨大，可开发容量达到30亿千瓦。与陆地风电不同，海上风电更少受土地占用、噪声污染等因素制约 [1] [2] [3] [4] 。我国传统能源主要分布在西北内陆，而经济重心在东南沿海，海上风电不仅可以弥补能源供应的不足，而且能够有效降低火力发电导致的环境污染，“向海争风”正成为东部沿海地区绿色低碳发展的“蓝色动力” [5] 。根据《2021年中国海洋经济统计公报》 [4] ，我国海上风电新增容量连续多年领跑全球。据国际能源署预测，2040年我国海上风电装机容量将与整个欧盟相当，减排能力将进一步提升。

近年来，随着国家和沿海省份高度重视开发潜在的海上风能，风电设备大型化，风电生产基地正在向沿海转移。然而，随着风机高度不断突破，传统的气象站点测风手段已无法满足监测需求。为了在风电场设计和运行过程中避免资金、人力、物力投入的浪费，亟待深入研究和揭示沿海地区风的时空变化规律 [6] 。国内外学者在这方面已开展了一些研究。如钱莉等 [7] 利用河西走廊东部设置的气象自动观测站收集的1971~2004年期间逐小时观测数据，发现相较于年平均风速仅2.7 m/s的北部地区，南部地区大于等于3.5 m/s，年平均有效风能时数大于5000 h，风能资源丰富。基于测风塔观测，Serhat Kucukali等 [8] 分析了2008年6月至2009年6月伊兹米特地区50米高度处收集的风数据，发现当地盛行东北偏东风向，湍流强度与高度呈现负相关。王丽娜等 [9] 通过气象站常规观测数据对大连地区风能资源进行评估时发现，大连多年平均风速为3.8 m/s，沿海的风速风能最为丰富，尤其是旅顺地区，30年平均风功率密度高达191.7 W/m²，是大连内陆地区的3~4倍。龙强等 [10] 针对放置在渤海湾北岸地区的两座测风塔的观测数据进行比对，发现海上的测风塔相对于沿岸的测风塔探测得到的平均风速更大，在110 m高度能达到7 m/s的风速；而湍流强度刚好相反，海上的湍流强度普遍小于沿岸的湍流强度，且随着高度的上升湍流强度的变化愈发平缓，70 m至110 m之间湍流强度基本稳定在0.08左右。

综上可见，以往研究人员在进行风能研究时多基于测风塔或气象站等传统观测手段，垂直高度信息有限、时间和空间分辨率低、设备维护成本昂贵、后期人员维修困难等不足，使得这方面的研究很缺乏，特别是针对电机垂直高度上的风能变化规律的研究尚罕见。而多普勒测风激光雷达可以提供秒级时间分辨率、米级空间分辨率的监测数据，能够弥补传统手段的不足。因此，本研究利用基于福建三沙和沿海两部多普勒测风激光雷达探测到的数据，对在风机典型高度100~200米的湍流强度、风功率密度等风能特性日变化和月际变化规律进行研究，以此来更具体地了解电机垂直高度上的风能变化规律，以期为沿海风能的监测与开发利用提供新的参考依据。

2. 研究方法

本研究在福建省宁德市三沙镇(经度119.9205，纬度26.8636)展开。福建省宁德市，地处我国东南沿海地区，海洋资源得天独厚，海域面积4.45万平方公里，大陆海岸线1046公里，占全省三分之一左右。海洋渔业资源丰富，是福建省乃至中国重要的海水养殖区域，约60万人从事海上养殖、捕捞和加工销售，渔业已成为当地沿海乡村振兴的重要支柱产业。由于常年受海洋季风气候影响，该地区拥有储能丰富的风资源，年平均风速可达到6 m/s [11] 。

整个观测过程共配备三台仪器，分别是两台不同型号的多普勒测风激光雷达(以下简称激光雷达)和一部安装在测风塔上的超声风温仪。Wind PrintV2000型号激光雷达放置在三沙气象站，代表内陆区域。Wind PrintS4000激光雷达安装在台风观测基地，代表沿海区域。测风塔相距三沙气象站300米左右的距离，距离观测基地几百米。

研究首先验证了激光雷达在边界层风场的探测精度，将安装在测风塔最高70米处的超声风温仪的观测数据作为真值，与激光雷达77.9米的风场数据进行比对，因高度非常接近故视为在同一高度观测。在此基础上，本研究针对湍流强度、风功率密度等风能特性日变化和月际变化规律进行分析，重点探究100米至200米风机典型高度上的风能变化特征。

从以往对福建沿海地区的研究来看，秋季风能资源丰富，其有效风能占比高，且秋季风能变异系数低，稳定性好，同时夏季是风能最丰富的时段 [12] 。据文明章等人研究 [13] ，福建沿海地区8~11月份的平均风速以及平均有效风功率密度在年变化中处于中高值阶段，其中10月份风能资源最丰富，平均风速和平均有效风功率密度分别为9.1 m/s和742.0 W/m²；朱琳等人的研究表明 [14] ，福建沿海地区8~11月份的风速和风能密度变化趋势相同，在年变化中处于中高值阶段，8月最低，11月最高。因此，本研究针对2020年8月至11月进行研究。

2.1. 多普勒测风激光雷达

多普勒测风激光雷达的原理是向大气中发射激光信号，与大气中的组分发生作用，利用多普勒频移与风速之间的关系，测量多普勒频移，进而反演计算出风场 [15] ，而反演技术是通过多普勒光束定向摆动(Doppler-Beam-Display, DBS)和速度方位显示(Velocity-Azimuth-Display, VAD)这两种扫描方式进行。多普勒测风激光雷达跟传统观测风装置相比，克服了传统装置使用寿命短、使用环境要求高、高度受限、分辨率低等缺点，因此大大提高了监测数据的测量精度和测风准确性 [16] 。

本研究采用的多普勒测风激光雷达是青岛华航环境科技有限责任公司生产的WindPrint S4000型号和WindPrint V2000型号的激光雷达(图1)。这两种型号的激光雷达都采用DBS扫描方法，即在固定角度的情况下，通过不同角度方位组合的探测来实现测量的一种技术方法。WindPrint S4000激光雷达能探测4000 m以下的径向风速、风廓线和垂直风速。时间分辨率5 s，空间分辨率15m，测风精度在0.100 m/s以内，工作温度范围为−30℃~+50℃，有着体积小、功耗低，安装方便、启动快捷的优点；WindPrint V2000激光雷达同样能探测垂直高度3000 m以下的风速与风向数据。时间分辨率3~4 s，风速测量范围0~75 m/s，工作温度范围为−40℃~60℃，具有重量轻、方便携带测量和利用率高等特点。

2.2. 超声风温仪

本研究中配备的超声风温仪是由Gill公司提供的WindMaster Pro型三维超声波风速风向计(图2)，分别在测风塔10 m、30 m、50 m和70 m高度处各分布一台，以此探测各高度处的风场情况。超声风温仪能提供声温、声速和UVW风矢量输出，可在0~65 m/s风速下正常工作，风速分辨率达到0.010 m/s，风向的探测范围为0˚~359˚。此外，从数据采集频率来看，该仪器可达到20 Hz的观测频率。

2.3. 监测过程与数据处理

为了了解观测期间福建三沙地区的气候状况，采用当地气象自动站逐小时的探测数据进行分析(图3)，发现观测期间自动站收集到的最高气温是在9月14日的11:00，气温为32.4℃，最低气温则是在10月22日凌晨1:00，为13.1℃。8、9月份平均气温为27℃，气温总体上保持在23℃~30℃范围，较为稳定；10月份平均气温为22℃，气温上下波动幅度较大。8、9月份平均相对湿度82.4%，相对湿度总体维持在70%~90%之间，10月份平均相对湿度76.1%，整体波动较大，相对湿度总体处于50%~80%区间。

为确保数据处理过程可靠，需要针对个别变量进行质量控制，例如水平风速大于50 m/s；垂直风速大于20 m/s或小于−20 m/s；风向小于0˚或大于360˚；信噪比小于7的数据都要筛选出来剔除 [17] 。为评估三沙地区的风能资源情况，本研究决定从湍流强度、风功率密度等风能特性展开分析，具体的计算公式如下：

$TI=\frac{{\sigma }_{\mu }}{\bar{\mu }}$ (1)

$k={\left(\frac{{\sigma }_{\mu }}{\bar{\mu }}\right)}^{-1.086}$ (2)

$c=\frac{\bar{\mu }}{\Gamma \left(1+\frac{1}{k}\right)}$ (3)

$p=\frac{1}{2}\rho c^3\Gamma \left(1+\frac{3}{k}\right)$ (4)

注：TI为湍流强度， ${\sigma }_{\mu }$ 为10分钟风速标准偏差， 为10分钟平均风速，k为威布尔分布的形状参数，c为威布尔分布的尺度参数，p为风功率密度， $\rho$ 为空气密度， $\Gamma$ 为gamma函数。

3. 结果与讨论

3.1. 数据情况概述

福建省宁德市三沙镇地处我国东南沿海地区，常年受海洋季风气候影响，拥有储能丰富的风资源，年平均风速可达到6 m/s，因此监测该地区的风能具有代表性 [18] 。如图4，实测数据结果显示S4000型和V2000型两部激光雷达在77.9 m高度探测到的三沙地区的风速与放置在测风塔70 m上的超声风温仪观测的风速相关性均达到0.9以上，整体的变化趋势较为一致。其中，8、9月份的风速总体处于0~8 m/s范围，偏小，且上下波动剧烈；10、11月份的平均风速集中在6~12m/s区间且无明显波动，说明10、11月份的风速相对稳定的。S4000激光雷达于9月17日22:30探测到最小风速0.149 m/s，在8月26日18:30探测到最大风速18.048 m/s。同一高度，V2000激光雷达在9月17日的23:10观测到最小风速0.247 m/s，在11月6日观测到风速峰值17.835 m/s。

进一步对计算得到的湍流强度、风功率密度作日变化平均处理时发现(图5)，三部仪器在77.9 m高度探测的湍流强度从总体上而言均是经历了0时至9时增大、9时至24时减小的日变化过程，S4000激光雷达、V2000激光雷达和超声风温仪探测到的湍流强度各自在9:10、8:10、9:00时刻是最大的，分别为0.190、0.202、0.258，在19:30、20:40、21:40三个时刻湍流强度都是最小的，仅为0.120、0.135、0.135。风功率密度的日变化也有着相似的变化趋势，7时至9时是风功率密度最小的时段，下午15:30~19:00是风功率密度最大的时段。湍流强度在9时达到极值的可能原因为，随着时间的推移地面温度不断升高大气温度也不断升高，地面与海平面上空的大气温度差异逐渐缩小，并且空气密度差异也不断缩小。不论是湍流强度，还是风功率密度，两部激光雷达和超声风温仪的计算结果都具有很强的相关性，湍流强度方面超声风温仪与两部激光雷达之间的相关性均达到了0.65以上，风功率密度方面相关性也是都超过了0.8，由此也验证了两部激光雷达对湍流强度、风功率密度的测算是真实可靠的。

3.2. 日变化

为继续探究三沙地区风能的日变化垂直分布特征，本研究选取了103.9 m、155.9 m、207.8 m三个典型高度分析整个试验期间8月至11月的湍流强度、风功率密度的日变化平均情况。

湍流强度通常随地形、天气、季节变化等因素发生改变，能有效地反映风的脉动变化特征 [19] 。湍流强度作为风资源前期评估工作中的重点评估对象，其结果对风电场选址、风机设置起着重要参考价值 [20] 。如图6，两部激光雷达在三个高度估算的湍流强度其日变化平均情况的走势较为一致。S4000激光雷达计算得到的湍流强度有明显的日变化，0时至9时湍流强度上升，9时至24时湍流强度下降，V2000激光雷达的湍流强度的日变化时间相比S4000激光雷达出现滞后，可能与两部雷达放置位置不同有关，0时至11时湍流强度呈上升趋势，11时至24时湍流强度呈下降趋势。

由表1的统计结果发现，整体的湍流强度是随着探测高度的增大而减小的，且随着高度的升高湍流强度减小的趋势也在趋于平缓，S4000激光雷达从103.9 m到155.9 m湍流强度降低13.534%，从155.9 m上升到207.8 m湍流强度降幅12.174%，同样的，V2000激光雷达计算的湍流强度从103.9 m到155.9 m减小22.819%，从155.9 m到207.8 m降幅的趋势进一步减缓，只有15.652%，这与龙强的研究结论是一致的。从探测位置来看，沿海放置的S4000激光雷达103.9 m探测高度处的湍流强度相比内陆监测点的V2000激光雷达同一高度处的湍流强度要小，可能是因为陆地地形更为复杂，下垫面更为粗糙，使得风速减小所导致的。155.9 m和207.8 m两个高度沿海和内陆监测点的湍流强度大小基本相同。从日变化平均情况来看，上午时段的湍流强度总体大于下午时段的湍流强度，S4000激光雷达在三个高度上探测到的湍流强度日变化的最小值的时刻均是19:50，V2000激光雷达探测到最小湍流强度也是稳定在19:30或19:40时刻。由此看出在100~200 m高度区间内的风资源较传统测风设备所测风资源更为丰富。

图7表现了S4000型、V2000型激光雷达湍流强度在三个垂直高度上每个月的日变化，结果表明两部激光雷达的测算结果具有较高的相关性。无论是S4000激光雷达，还是V2000激光雷达，在三个不同高度上四个月的湍流强度日变化都存在类似的走势分布，日出时段的湍流强度达到最大，日落时分的湍流强度则最小。

风场的可靠预测是风电场优化设计和运行的重要基础，而风功率密度是预测风场功率输出的关键因素 [21] 。相比湍流强度，两部激光雷达在三个高度上估算到的风功率密度的日变化平均情况整体吻合程度更高(图8)。S4000激光雷达计算到在0时至9时、19时至24时两个时间段内风功率密度随着时间的变化而减小，9时至19时随时间的增大而增大，V2000激光雷达的探测结果显示在0时到10时和19时到24小时两时段内风功率密度呈减小趋势，10时至19时内呈上升走势。

结合表2可以发现，从探测高度来看，高度越高风功率密度也更大，但随着高度的提高风功率密度增大的趋势也在减缓，S4000激光雷达在103.9 m到155.9 m探测高度内风功率密度上涨幅度达20.46%，而从155.9 m到207.8 m之间风功率密度涨幅仅为10.814%。V2000激光雷达测算的风功率密度从103.9 m高度到155.9 m高度，中间涨幅有19.207%，而从155.9 m到207.8 m高度只上涨了6.407%。可以看出风功率密度在100~200 m高度区间内是逐渐增大的，涨幅速度却趋于平缓。从探测位置来看，沿海的风功率密度整体上均低于内陆监测点的风功率密度，与王丽娜等人关于沿海与内陆监测点之间风功率密度的结论相矛盾，这可能是由于地形因素造成激光雷达受遮挡所造成的。从日变化平均情况来看，上午时段的风功率密度明显小于下午的风功率密度。两部激光雷达计算得到的最大风功率密度的时间稳定在18:50左右，风功率密度最小的时段集中在9:00~11:00之间，都与湍流强度的分布特征呈现相反的趋势。

图9展现了S4000型、V2000型激光雷达风功率密度每个月的日变化在垂直高度上的分布，计算结果显示两部激光雷达对风功率密度的测算一致性是比较好的。8、9月份的风功率密度日变化总体上没有过大的波动，基本是稳定在100 W/m²左右，10、11月份风功率密度日变化上下产生较大波动，下午时段的风功率密度相较于上午时段的风功率密度整体有50%的上涨。

3.3. 月际变化

结合月平均湍流强度图10发现，155.9米和207.8米高度上湍流强度最大值出现在9月，湍流强度

最小值出现在11月，推断可能与9月份的风速都小于11月份的水平风速有关。从探测高度看，随着高度的上升该地区的湍流强度也在不断的减小，且减小的趋势也是趋于缓和，103.9 m处每个月的平均湍流强度范围在0.12~0.15之间，155.9 m和207.8 m的平均湍流强度范围则稳定在0.09~0.12。从探测位置来看，8、9月份内陆监测点的湍流强度整体要高于沿海边的湍流强度。

结合月平均风功率密度图11看，随着高度的上升该地区的风功率密度也在不断的增大，而增大的趋势也是在减缓，在103.9 m到155.9 m之间风功率密度上涨了超过20%，而从155.9 m到207.8 m上涨幅度仅有12%左右。从探测位置来看，内陆监测点每个月的平均风功率密度都是大于沿海的这个可能是因为受地形干扰或者是受风向影响。相比8、9月份，10、11月份的风能资源更加丰富，结合季节来看，8月份属于夏季，9、10、11月份属于秋季，整体风能资源情况夏季少，秋季多，尤其是秋末时段风功率密度最盛。

4. 结论

本研究利用位于福建三沙的S4000型和V2000型多普勒测风激光雷达，分析了我国典型的东南沿海地区——福建三沙内陆和沿海两个监测点在风机典型高度100-200米的湍流强度、风功率密度的分布特征，为后续高效开发利用沿海风能提供科学参考依据。主要研究结论如下。

1) 观测期S4000多普勒测风激光雷达探测到的77.9米高度的最大风速是18.048 m/s，V2000多普勒测风激光雷达探测到的77.9米高度的最大风速是17.835 m/s。

2) 随着探测高度的上升，两部多普勒测风激光雷达测算的湍流强度逐渐减小、风功率密度逐步增大，且减小或增大的趋势趋于平缓。一般来说，为了产生更高的电能，风力涡轮机都会建在风功率密度更高的地方，考虑到成本问题，建议将风机安装高度设置在150米左右。

3) 风功率密度和湍流强度的日变化基本呈现相反的分布，日出前湍流强度最大，风功率密度最小；日落后湍流强度最小，风功率密度达到最大。8、9月份三个高度上的风功率密度基本稳定在100 W/m²左右，风能资源相对较少；10、11月份风功率密度随着高度上升而增大，具有丰富的风能储能资源。考虑到8、9月份属于夏末秋初时节，太阳能资源较多，因此建议将风能和太阳能进行结合来提供发电，而10、11月份则充分利用风能资源进行发电。

基金项目

科技部重点研发计划(2018YFC1506303，2018YFC1506305)；科技部重点研发计划国际合作专项(2017YFE0107700)；国家自然科学基金(41775065，41475060，41805088)；本研究同时受上海气象局台风科技创新团队支持。

NOTES

^*第一作者。

^#通讯作者。

参考文献