2. 中国气象局地球系统数值预报中心, 北京 100081;
3. 灾害天气国家重点实验室,北京 100081;
4. 龙源(北京)风电工程设计咨询有限公司,北京 100034
2. CMA Earth System Modeling and Prediction Centre, Beijing 100081;
3. State Key Laboratory of Severe Weather, Beijing 100081;
4. Longyuan (Beijing) Wind Power Engineering Design Consulting Co., Ltd., Beijing 100034
由于化石能源的不可再生性, 以及其在燃烧过程中排放的温室气体和污染气体带来的气候环境问题日益严重,开发、利用可再生能源已成为世界各国急需解决的问题(Dincer and Acar, 2015;李国庆和李晓兵,2016;陈正洪等,2018)。风能作为一种洁净的、储量极为丰富的可再生能源,越来越受到广泛重视(Garrigle and Leahy, 2015;Phillips,2015;叶小岭等,2019)。我国风能资源储量丰富,总储量为322 600 GW,实际可开发量为25 300 GW(薛桁等,2001;梁红等,2021)。近年来我国风电发展迅速,据全球风能协会公布的数据显示,2021年中国风力发电装机总量达338.31 GW,占世界总装机容量的40.4%;其中海上风电发展迅猛,2021年新增装机量16.9 GW,约占全球新增装机量80%。由于风能的能量密度小,发展风能需要较大的集能面积,因此大规模建设风力发电场会改变地表粗糙度,从而改变陆表和大气的热交换过程; 同时由于风机涡轮振动、风能转化为电能,改变了自然界原有的能源循环模式,导致风电场对局地或区域气候可能产生影响(Buckley et al,2005;张雪婷等,2019;蒋俊霞等,2019)。由此,大型风电场建设和运行对气候和环境的影响受到越来越多的关注(Fitch,2015;Li et al,2018;梁红等,2021)。
关于大型风电场建设运行对局地风环境的影响已有大量研究。Christiansen and Hasager(2006)利用SAR雷达数据发现丹麦地区海上风电场的平均风速减小8%~9%,对下风向的影响范围达20 km;Fitch et al(2012)通过数值模拟发现从海上风电场延伸至其60 km下游区域均有风速衰减现象,风电场区域内风速衰减最多可达16%;Smith et al(2013)对美国中西部大型陆上风电场观测研究显示发现在叶片扫过区域及风机轮轴高度处风速明显降低、湍流强度显著增强,且夜间更为明显。胡菊(2012)模拟发现我国河西走廊地区的大型风电场建成后平均风速减小了0.3 m·s-1;徐荣会(2014)根据观测发现内蒙古地区风电场局地风速在风机的扰动下有减小趋势;张雪婷等(2019)发现我国大别山区风电场群建设和运行使大悟县局地风速略有减小。但国内研究主要集中在陆上风电场,对海上风电场研究较少,且关于风电场建设对一些重要指标如湍流强度和风切变指数的影响研究也很少。
湍流强度(TI)和风切变指数(WSE)是评价气流稳定程度的指标,对于风电场的风能评估及风电场稳定运行具有重要意义(李鸿秀等,2014;黄林宏等,2016;龚玺等,2018;陈燕等,2019;靳双龙等,2019)。TI不但与地表摩擦和风速切变引起的动力因子以及温度层结引起的热力因子密切相关,还与地理位置、地形和地表粗糙度等因素有关(贺德馨,2006;李鸿秀等,2014;刘敬乐等,2020);WSE表征风速随高度的变化程度,主要受地表粗糙度、大气稳定度等因素影响(徐宝清等,2014;龚强等,2015;马晓梅等,2020;常蕊等,2022;程雪玲等,2023)。目前利用测风塔资料开展了风电场建设前风环境特征如风速、TI或WSE等分析(李鹏和田景奎,2011;班欣等,2012),或者是利用数值模拟技术开展风电场对大气边界层特征影响研究(刘维维,2015;路婷婷,2017;张颖超等,2016;常蕊等,2013),但对风电场群建设后局地风速以及TI和WSE变化的观测研究较少。由于风电场建设后,局地风速和下垫面粗糙度发生变化,且风机的扰动改变了局地大气的稳定性,相应改变TI和WSE(吴正人等,2014;靳双龙等,2019;夏馨等,2022),从而改变了风电场建设前的风能评估结果以及风电场安全运行所需的风速条件(刘磊等,2012);同时由于陆/海风电场存在着下垫面及周围环境的巨大差异,因此有必要开展二者在风电场建设后对局地风环境关键参数影响及对比研究,为进一步阐明陆/海风电场对局地风环境的不同影响机理提供参考。
由此,本文分别选择河北尚义陆地风电场和江苏如东海上风电场,利用风电场群周边国家气象站风速观测资料以及质量控制后的测风塔测风资料,选定风电场建设前后典型时段,通过分析局地风速、TI和WSE的变化,开展陆/海风电场群建设对局地风场影响差异初步研究。一方面阐明典型陆/海下垫面风电场群建设后局地风环境关键参数的变化规律,另一方面揭示这两种不同下垫面风电场群对风环境关键参数的影响差异机制,从而为今后陆/海风能资源评估、有效规避相关风险等提供科学依据。
1 研究区与数据 1.1 研究区概况本文研究选取的陆地风电场位于中国河北张家口尚义县(如图 1a所示),尚义风电场群从2008—2015年底陆续建成224台风机(表 1),风机轮毂高度为65~90 m,总装机容量为450 MW;风电场所在区域属于半干旱大陆性季风气候区,海拔高度为1375~1573 m,自然植被以草本植物为主,山谷及缓坡地带有零星小面积成块林地、耕地及灌木;气候特点表现为气温低温差大、雨量少而集中、光照充足,年均降水量为414 mm,年均风速为3.4 m·s-1,全年主导风为西北向(NW),次主导风为西西南(WWS)。
海上风电场研究区位于中国江苏南通如东县(如图 1b所示),如东风电场群从2009—2015年陆续建成155台风机,风机轮毂高度为65~90 m,总装机容量为482 MW;风电场所在区域属于亚热带海洋性季风气候区,境内地势平坦,海拔高度为3~5 m,陆地地貌是典型的滨海平原,沿海属黄海南部海域,沿海潮汐24小时涨落两次;受海洋调节和季风环流影响,四季分明、气候温和、光照充足、雨水充沛,年降水量为1100 mm,年均风速为3.2 m·s-1,夏季盛行东南风,冬季盛行偏北风。
1.2 数据资料及分析时段河北尚义陆地风电场群气象资料包括:(1)国家基本气象站(尚义、商都、怀安和兴和站)1981—2018年风速、气温、相对湿度逐日观测数据;(2)尚义0704#测风塔2006—2017年30、50、70 m高度逐时风速观测数据。
江苏如东海上风电场群气象资料包括:(1)国家基本气象站(如东、如皋站)1981—2018年逐日风速观测数据;(2)如东0865#测风塔2009—2018年15、50、80 m高度逐时风速观测数据。
根据风电场群建设和投产时间(表 1)以及气象行业标准《风电场气象观测及资料审核、订正技术规范(QX/T 74-2007)》(中国气象局,2007)进行质量控制后的气象观测资料,选择合适的陆/海风电场建设前后分析时段(表 2)。尚义风电场选用建设前后各3年(2005—2007年和2016—2018年)的国家气象站资料进行局地风速变化分析,而0704#测风塔资料由于存在不连续和缺测情况,选择数据较为完整的2006年和2017年进行TI和WSE变化分析。如东风电场选用0865#测风塔2009年和2018年数据进行TI和WSE变化分析。
河北尚义陆地风电场群周围60 km内包含4个国家气象站:尚义、商都、怀安和兴和站(图 1a),其中尚义站紧邻风电场群,可作为风电场影响站进行评估;此外,尚义0704#测风塔位于风电场群内,也可作为风电场影响站进行评估。
利用尚义站评估尚义陆地风电场局地风速变化时,除了风电场群的可能影响,气象站气象要素变化还受年际气象条件波动和气候变化影响,仅靠单一气象站是不能识别出风电场对其单独影响,需要借助研究站与背景站(风电场非影响区)之间的对比分析来消除年际波动与气候变化影响。由此,选择合适的背景站是尚义陆地风电场局地风速变化评估的一个关键步骤。
背景站一般要求与研究站的气候类型和地理条件相似,但又不能受风电场影响;数值模拟研究表明风电场影响风速衰减的距离在风电场下游方向可以达到30~60 km(Frandsen et al, 2009),要求背景站距风电场至少30 km。这里分别对风电场建设前(1981—2007年)的尚义站与怀安、兴和、商都站的主要气候要素和地理条件进行相似性分析,如表 3所示。
由表 3分析可知,尚义站年均气温、风速和相对湿度均与商都站最高,且两者海拔差最小(仅44 m),空间距离在30 km以上,且商都站在风电场主导风向(NW)上游区域,因此可选择商都站作为背景站进行尚义陆地风电场的局地风速变化分析。
为进一步分析尚义陆上风电场的影响范围,还需对可能受到风电场影响的其他两个站(怀安与兴和站)进行分析。由于怀安站和背景站海拔相差很大(583 m),空间距离也最远(129 km),且怀安站与背景站风速相关系数仅为0.65(表 3),以商都站作为背景站并不能很好分析怀安站风速变化;而且怀安站与尚义风电场群相距约60.7 km,已超出风电场最大可能影响范围(60 km)(Fitch et al, 2012), 因此,本文不考虑怀安站作为影响站进行分析。除了选择尚义站作为影响站外,还选择了距离风电场群约29 km、年均风速与背景站相关系数为0.86的兴和站作为可能影响站进行分析。
2.1.2 海上风电场影响站选择距离如东海上风电场群最近的国家气象站为如东和如皋站,分别距风电场约21 km和50 km(图 1b)。虽然如东站可能在海上风电场影响范围内,但是对如东站1981—2018年年均风速资料分析发现其在2014年风速明显下降且低于如皋站,而在1981—2013年期间如东站风速均明显高于如皋站;根据文中表 1可知,2014年前如东海上风电场已投产94个风机,2014年以后投产61个,因此如东站的风速突然降低很可能是其他外部原因造成,经调查得知2014—2015年如东气象站周围建立了大量高楼,致使2014年后风速下降明显,并非风电场影响造成;而如皋站距离风电场较远(50 km),且在海上风电场群的西南方向,与风电场群主导风向(东南—偏北风)相差甚大,几乎不受风电场影响,因此如东站和如皋站均不适宜作为风电场影响站进行评估,即如东海上风电场没有合适的国家气象站进行局地风速变化评估。
如东0865#测风塔位于风电场群偏东方向约7 km处,处于风电场海上影响范围之内,可作为影响站进行评估。
2.2 局地风速影响评估定义尚义陆上风电场群对局地风速的影响值为:
$ \mathrm{d} V=V_{\mathrm{I}}-V_{\mathrm{B}} $ | (1) |
$ \Delta V=\mathrm{d} V_{\mathrm{AF}}-\mathrm{d} V_{\mathrm{BF}} $ | (2) |
式中:dV为国家气象站影响站(或可能影响站)10 m高度风速(VI)与背景站10 m高度风速(VB)之差(单位:m·s-1);ΔV为风电场建设后(dVAF)相对于风电场建设前(dVBF)的风速变化(单位:m·s-1),即风电场对风速的局地影响。
为评估风电场群建设前后局地风速是否发生显著变化,这里以风电场建设前影响站(或可能影响站)与背景站风速差dV的年际标准差(StddV)作为年际波动,如果风速差变化值(ΔV)的绝对值大于StddV,则表明该地风速受到风电场群影响,如果ΔV绝对值小于或等于StddV,则表明该地风速没有受到风电场群影响。
由于如东海上风电场没有合适的国家气象站进行局地风速研究,虽然0865#测风塔可作为影响站进行评估,但由于该塔建设前仅有1年资料,且没有10 m高度的测风数据,因此本文不开展如东海上风电场的局地风速影响研究。
2.3 局地湍流强度(TI)影响评估根据国家标准《风电场风能资源评估方法(GBT 18710—2002)》(中华人民共和国国家质量监督检验检疫总局,2004),湍流强度的计算公式为:
$ \mathrm{TI}=\sigma / V $ | (3) |
式中: TI为湍流强度,σ为平均风速标准差,V为平均风速。
本文定义风电场对TI的局地影响为风电场建设后TI与建设前TI的距平百分率。该值为正值时,表明风电场建设运行使得TI增加;该值为负值时,表明风电场建设运行使得TI降低。
针对尚义陆上风电场,分别计算尚义0704#测风塔风电场建设前后30、50和70 m高度TI变化,开展陆上风电场对局地TI的影响评估;针对如东海上风电场,分别计算如东0865#测风塔风电场建设前后15、50和80 m高度TI变化,开展海上风电场对局地TI的影响评估。
2.4 局地风切变指数(WSE)影响评估风切变指数(WSE)表示风速垂直切变的强度。在这里,采用风速随高度变化呈指数规律分布(王承煦和张源,2003;Farrugia,2003;Rehman and Al-Abbadi, 2007)进行估算,即:
$ V_{\mathrm{n}} / V_1=\left(Z_{\mathrm{n}} / Z_1\right)^a $ | (4) |
式中:Vn、V1分别为Zn、Z1高度上的风速,α为风切变指数WSE。
对式(4)作变换,并两边取对数,即得:
$ \alpha=\frac{\ln \left(V_{\mathrm{n}} / V_1\right)}{\ln \left(Z_{\mathrm{n}} / Z_1\right)} $ | (5) |
目前计算WSE时,不同风速资料的选择会造成估算结果偏差明显。彭怀午等(2010)、吴婷婷(2015)研究表明采用全部数据进行风切变计算产生的平均值代表性差、偏离系数较大,而采用去除3 m·s-1以下风速数据计算产生的平均值代表性要优于前者、偏离系数较小。因此,本文利用测风塔资料,选取不同高度去除3 m·s-1以下且高层风速大于低层风速的数据,根据式(5)计算每两个高度层间的WSE,再利用最小二乘法(徐宝清等,2014;陈燕等,2019;孙淑芳和张广兴,2019)将各层间的WSE进行拟合,得到该区域平均WSE,即根据式(5)可得到:
$ \ln \left(V_{\mathrm{n}} / V_1\right)=\alpha \ln \left(Z_{\mathrm{n}} / Z_1\right) $ | (6) |
令
本文定义风电场对WSE的局地影响为风电场建设后WSE与建设前WSE的距平百分率,该值为正值时,表明风电场建设运行使得风切变增加;该值为负值时,表明风电场建设运行使得风切变下降。
3 结果分析 3.1 陆地风电场群局地风速变化 3.1.1 年变化趋势就影响站尚义站而言,2005—2011年影响站和背景站(商都站)风速变化趋势具有良好的一致性(图 2a),2011年后,背景站风速(VSD)总体呈上升趋势,而影响站风速(VSY)则呈下降趋势,意味着2011年后影响站受到非气候变化因素导致其风速变化趋势与背景站差异显著。从两站风速差(dV)变化来看,2012年后两站年风速差呈增大趋势,其中2016—2018年风速差明显大于风电场竣工前任何一年(图 2a),该结果显示尚义站很可能受到风电场群影响,致使风速在风电场建设前后发生明显变化。
就可能影响站(兴和站)而言,2005—2018年VXH与VSD变化趋势具有较好一致性(图 2b),其中2014—2016年两站风速差较大;2017—2018年,风速差回归到2012年以前水平,这与尚义风电场建设期(2008—2015年)时间匹配并不一致。该结果显示,风电场建设前后兴和站风速变化并不显著,兴和站受尚义风电场群影响的可能性很小。
3.1.2 局地风速变化图 3表示尚义站(SY)与兴和站(XH)年和四季风速差变化值(ΔV)与气候年际标准差(StddV)。就影响站(SY)而言,风电场年均ΔV为-0.43 m·s-1,绝对值高于StddV(0.26 m·s-1),表明受风电场群影响,SY年均风速明显降低;就四季平均来看,ΔV在春季变幅最大(-0.89 m·s-1),其次是冬季(-0.71 m·s-1),且绝对值均明显高于相应时段的StddV,而夏季和秋季的ΔV在StddV波动范围内,表明风电场建设后风速降低主要发生在春季和冬季。
就可能影响站(XH)而言,年和春、夏、秋季ΔV均在StddV波动范围内,仅冬季ΔV略高于StddV(偏高0.04 m·s-1),表明风电场建设后兴和站风速变化基本处于正常年际波动内。该结果显示兴和站受尚义风电场影响很小,即使受到影响也仅仅发生在冬季,而且这种影响非常轻微。由于兴和站位于风电场偏南方向约29 km处,因此可以推测尚义陆上风电场群影响范围很可能不超过30 km。
3.2 陆/海风电场群湍流强度变化 3.2.1 日变化分别选择尚义陆地测风塔30 m、70 m和如东海上测风塔15 m、80 m高度层,分析建设前后TI日变化(图 4):无论是尚义陆上风电场群,还是如东海上风电场群,建设后各高度层不同时次TI均明显增加;这主要是由于风电场群内大范围风机的搅动,使得空气上下交换加强,致使陆/海风电场区TI增加明显(Christiansen and Hasager, 2006;Smith et al,2013)。就不同高度来看,低层TI明显高于高层,这是由于在低层大气更易受到地表粗糙度与地形影响,大气上下交换强,因此TI更强。
尚义陆上风电场TI日变化呈单峰曲线,白天TI明显大于夜间,这是由于日出后地表受热,对流逐渐加强,导致大气层结不稳定,TI逐渐增大,至中午前后达到最大;日落后地表受热停止,对流作用减小,大气层结逐渐转为中性,TI也随之减小。与陆上风电场并不相同,如东海上风电场TI日变化不明显,这主要是由于海面水体热容量大,增温和降温过程较陆地缓慢,导致海上区域热力对流变化小于陆地,因此海上风电场TI日变化波动较小。
对比尚义陆地和如东海上风电场群同一高度层50 m建设前后TI日变化, 如图 5a所示:无论建设前后,陆地TI日均值均明显高于海上,尤其在白天,这是由于海洋热容量明显高于陆地,造成海上白天近地层温度梯度小于陆地,加之下垫面平滑,从而使得同一高度上海洋TI明显小于陆地;到了夜晚,陆地和海水上大气层结都比较稳定,温度梯度接近,因此TI比较接近。
对比陆/海风电场群建设后各典型高度层(尚义:30 m、70 m;如东:15 m、80 m)TI增幅日变化(图 5b),可以看出:就不同时次来看,陆地风电场群各高度层TI较建设前增加11%~66%,其中夜晚增幅(平均为45%)明显大于白天(平均为20%)。这是因为陆地区域夜间大气稳定度高,对流本身较弱,但由于风机的搅动,使得上下空气交换加强,致使湍流强度大幅增加;而白天地面受大气加热影响,大气层结不稳定,对流本身较强,虽然风机的搅动能更进一步加强空气扰动,但相比于夜间增幅小得多。
就不同高度来看,陆地风电场建设后70 m高度TI平均增幅为40%,明显高于30 m TI增幅(28%),表明TI在高层增幅大于低层,这种增幅在夜间更为明显。这很可能一方面与70 m高度TI在夜间处于低值区且明显小于30 m高度(图 4a)有关;另一方面由于70 m高度层更接近风机轮毂高度(平均75 m),受风机搅动影响更大,而尚义风轮平均直径约为90 m,表明叶片最底端与地面距离约为30 m,因此30 m高度受风机影响较小。
如东海上风电场群建设后各高度TI增加29%~ 48%,但增幅没有明显的日变化特征。这是因为虽然海上区域风机搅动加强了空气扰动,使TI增加,但由于TI日变化较小(图 4b),TI增幅的日变化也较小。此外,80 m高度TI平均增幅为39%,稍高于15 m高度平均增幅(36%),表明海上风机对低层和高层TI影响造成的差异远不如陆地风机的影响明显。
3.2.2 年和季节变化尚义陆地和如东海上风电场群建设前后不同高度年均TI对比如图 6所示,与TI日变化相似,陆/海风电场群年均TI不同高度均较建设前明显增加,其中尚义陆地各高度建设前年均TI为0.099~0.114,建设后增至0.131~0.144,平均增加了31%;如东海上各高度建设前年均TI为0.077~0.117,建设后增至0.107~0.159,平均增加了37%,TI增加原因与前面日变化相同。对比陆/海风电场群同一高度50 m来看,无论建设前后,陆地区域TI均明显高于海上,偏高约21%,这是由于海面相对于陆地下垫面平滑、地表粗糙度小,加之海水热容量高于陆地、近海面层温度梯度小于陆地近地层,导致海上风电场区域TI相对于陆地风电场整体较小。
以尚义陆地和如东海上风电场群测风塔50 m高度TI四季变化进行对比,如图 7a所示:从季节变化上看,无论建设前后,陆地风电场TI四季变化均较为明显(建设前0.091~0.126,建设后0.121~0.155),其中夏季最大、冬季最小,春秋季介于之间;海上风电场TI季节变化建设前较为平稳(0.080~0.096),季节变化幅度明显小于陆地,这与已有的风电场建设前不同地形TI变化相关研究结果较为一致(李鸿秀等,2014;龙强等,2020);而风电场群建设后如东海上风电场TI季节变化加大(0.098~0.135),其中冬秋TI强于春夏,这与陆地风电场季节变化规律并不一致。
从季节TI增幅看(图 7b),尚义陆地风电场春季TI增幅(47%)最为明显,然后依次为冬季(33%)、秋季(33%)和夏季(24%),这与前面分析的春季和冬季风速受风电场影响显著具有一致性。而如东海上风电场冬季TI增幅最大(49%),然后依次为秋季(42%)、春季(34%)和夏季(21%)。
3.3 陆/海风电场群风切变指数变化 3.3.1 日变化尚义陆地和如东海上风电场群建设前后,年均WSE日变化如图 8a所示:两个风电场群WSE均表现为夜间高、白天低等日变化特征,但陆地风电场WSE日变化明显大于海上风电场。这是由于白天对流作用强、大气层结不稳定,湍流活动旺盛,使不同大气垂直层间的动量交换频繁,风速垂直梯度较小,导致WSE较低;夜间地表温度低,大气层趋于稳定,层间很少混合或不混合,风速随高度的上升急剧增大,导致WSE较高,这与雷杨娜等(2015)相关研究结果较为一致。
对比陆/海风电场群建设后WSE变幅的日变化分布(图 8b),可以看出:风电场建设后,陆地风电场群WSE在白天明显增加(平均增幅为33%),夜间则降低(平均降幅为9%),WSE变幅趋势与日WSE大小分布正相反,WSE小的时段增幅明显,WSE大的时段降幅较大,从而使得WSE的日变化值变小(从建设前的0.199降至建设后的0.133)。这是由于陆地区域WSE日变化较为明显, 白天大气不稳定、风速垂直梯度较小,但由于高层风速更接近风机轮毂高度受风机扰动影响更大,风速增大幅度大于低层,使得上下层风速梯度增大,致使白天WSE大幅增加;夜间大气稳定、风速垂直梯度较大,风机扰动使得上下空气交换加强、风速梯度变小,致使夜晚WSE降低。海上风电场群建设后,WSE在白天夜晚均有所增加(平均增幅为24%),但由于WSE日变化较小,WSE增幅无明显日变化趋势,且建设前后WSE日变化值也没有明显变化(分别为0.040和0.039)。
3.3.2 年和季变化尚义陆地和如东海上风电场群建设前后,年和四季平均WSE及建设后变幅如图 9所示。
风电场建设前,陆地区域年均WSE为0.196(图 9a),明显高于海上区域(年均WSE为0.133)。这是由于尚义处于内陆山区,地形起伏和地表粗糙度都较大,使得风速随高度变化大,因此WSE也相对较大;如东沿海地形平坦,气流受地形和粗糙度影响较少,风速在不同高度的分布较稳定,故WSE相对较小。陆地风电场WSE表现为冬季和秋季较大,春季和夏季较小,这是由于冬季气温低,太阳辐射强度较低,大气层结稳定度高,动量、热量交换较少,高低层之间风速差异大,从而WSE大;夏季气温高,高强度的太阳辐射引起大气层的不稳定性增强,大气层间的混合加强,高低层之间风速差异变小,导致WSE较小(杜燕军和冯长青,2010);海上风电场WSE表现为冬季和春季较大,夏季和秋季较小,这很可能与如东海上风电场夏秋温度高于冬春、大气混合较冬春更加明显有关。
尚义陆地风电场建设后,无论年和季节WSE均有所下降(图 9b),其中年均降低8%,四季降低3%~12%,这主要是由于风机的影响造成夜间WSE降低所致(图 8b),虽然夜间WSE降低幅度不如白天WSE增幅,但夜间WSE明显大于白天WSE(图 8a),从而造成年、季WSE有所降低,与尚义风电场相反,如东海上风电场建设后无论年和季WSE均明显增加,其中年均增加了24%,各季增加了4%~37%。这主要是由于海表粗糙度较小且远离风机,低层风速受到影响较小,高层风速更接近风机轮毂高度,受风机扰动影响更大,风速增大幅度大于低层,风速垂直梯度增大,致使WSE增加。
此外,陆/海风电场均对夏季WSE影响较小,这可能是夏季大气层间混合本身较强,而风机扰动对这种混合影响较小,使得WSE变化并不显著。
4 结论与讨论 4.1 结论本研究利用河北尚义陆地风电场群周边4个国家级气象站1981—2018年逐日观测资料、区域内测风塔2006—2017年不同高度逐时测风资料以及江苏如东海上风电场群区域测风塔2009—2018年不同高度逐时测风资料,根据风电场群建设和投产时间以及质控后可用的测风塔数据,选定风电场建设前后分析时段,通过计算建设前后局地风速差值、湍流强度TI和风切变指数WSE,对比研究陆/海风电场群建设对局地风环境关键参数的影响,初步得到以下主要结论。
(1) 尚义陆地风电场具有明显风速减弱效应,年均风速较建设前降低0.43 m·s-1,且风速降低主要发生在春季和冬季,风速分别降低了0.89 m·s-1和0.71 m·s-1,风电场影响范围很可能不超过30 km。
(2) 陆地风电场区TI值高于海上风电场区,TI日变化和季节变化也比海上区域明显。陆/海风电场群建设对TI均具有增强效应,建设后年均TI分别增加了31%和37%。陆地风电场区TI夜间增幅(45%)明显大于白天增幅(20%),且高层增幅(40%)大于低层(28%),而海上风电场区TI不同高度增幅及其日变化较为平稳。陆地风电场区TI四季增幅为24%~47%;海上风电场区TI四季增幅为21%~49%。
(3) 陆/海风电场区域WSE均呈现冬季高、夏季低、夜间大于白天等特征,但陆地WSE日变化明显大于海上。陆/海风电场建设后对WSE影响存在显著差异,陆地风电场建设后WSE白天增加了33%、夜晚则降低了9%,从而导致WSE日变化变小,且年和四季WSE均有所降低,其中年均降低了8%,四季降低了3%~12%;海上风电场建设后WSE白天、夜晚均明显增加,从而导致年和四季WSE均明显增加,其中年均增加了24%,四季增加了4%~37%。
4.2 讨论(1) 风电场群对湍流强度TI影响的不确定性。本文所得结果显示陆/海风电场群在建设前湍流强度TI年、季、日变化均存在明显差异,这也与已有风电场建设前相关研究结果(史军等,2016;郑新倩等,2019)较为一致,表明本文的研究资料具有一定代表性,方法较为合理。而风电场群建设后,陆上和海上TI均发生了较为明显的变化,凸显了风电场群对海/陆局地风场的影响,由于海陆两种下垫面性质的不同,这种影响存在着不同高度、白天/夜晚和季节差异,但这种差异是否受到其他条件例如年际变化或者背景气候变化的影响还并不清楚,因此未来还需更长时间的气象观测资料来开展风电场群对TI长期的影响研究。
(2) 风电场群对风切变指数WSE影响的不确定性。本文结果显示风电场建设前陆/海WSE呈现冬季较高、夏季较低、白天较小、夜间较高等特点,与已有风电场建设前的相关研究结果较为一致(彭怀午等,2010;孙淑芳和张广兴,2019;陈燕等,2019);风电场建设后由于风机一定程度上改变了不同高度的风速使得WSE发生变化,但这种变化实际上是很复杂的,因为WSE不仅受地表粗糙度影响,还受地形、大气稳定度、海拔高度以及热力层结的影响(李鹏和田景奎,2011),而风电场不仅影响地表粗糙度,还影响风速和大气温度如局地增温降风效应(张雪婷等,2019;Liu et al, 2021)。因此风电场群对WSE的影响机制未来还需要更多观测资料进行佐证。
(3) 资料分析的不确定性。由于测风塔资料限制,本文仅选取了风电场群建设前后各1年资料,而且由于陆/海测风塔资料时间差异,选取的年份也不同,因此存在资料代表性不够充分造成的不确定性问题,从而造成陆/海风电场群对局地风环境影响差异分析存在一定不确定性。此外,不同测风塔资料实际上是受气候背景影响并存在年际变化的,虽然TI和WSE是关于风速的相对变化指标,很大程度上已经消除了气候背景场(含年际变化)的影响,但由于没有周围参考塔同高度风速资料,陆/海风电场对TI和WSE的分析与对比仍存在一定程度受背景场影响的不确定性。另外,如东测风塔位于海上风电场东南方向,考虑盛行风向的季节变换(冬季盛行风为偏北风,夏季盛行风为偏南风),测风塔受风电场影响在冬季较夏季更为明显,这可能也是如东海上风电场对TI和WSE的影响在冬季最大、夏季最小的另外一个原因。因此,未来还需选取更多的站点以及更有代表性的资料进行详细分析,才能消除这种由于资料不足带来的不确定性。
班欣, 冯还岭, 祁欣, 等, 2012. 连云港沿海近地层湍流强度特征[J]. 气象科技, 40(2): 285-292. Ban X, Feng H L, Qi X, et al, 2012. Analysis of surface layer turbulence intensity along Lianyun-gang coastland[J]. Meteor Sci Technol, 40(2): 285-292 (in Chinese).
|
常蕊, 朱蓉, 柳艳香, 等, 2013. 基于均生函数的风电场风速短临预报模型[J]. 气象, 39(2): 226-233. Chang R, Zhu R, Liu Y X, et al, 2013. Nowcasting model of wind speed based on mean generating function for wind farms[J]. Meteor Mon, 39(2): 226-233 (in Chinese).
|
常蕊, 朱蓉, 赵大军, 2022. 登陆台风影响下离地300 m高度内的强风特征[J]. 大气科学, 46(5): 1071-1086. Chang R, Zhu R, Zhao D J, 2022. Strong wind characteristics of the lower boundary layer (0-300 m) during the landfall of a typhoon[J]. Chin J Atmos Sci, 46(5): 1071-1086 (in Chinese).
|
陈燕, 张宁, 许遐祯, 等, 2019. 江苏沿海近地层强风风切变指数特征研究[J]. 高原气象, 38(5): 1069-1081. Chen Y, Zhang N, Xu X Z, et al, 2019. The wind shear exponent in the near-surface strong wind in the coastal areas of Jiangsu Province[J]. Plateau Meteor, 38(5): 1069-1081 (in Chinese).
|
陈正洪, 何飞, 崔杨, 等, 2018. 近20年来风电场(群)对气候的影响研究进展[J]. 气候变化研究进展, 14(4): 381-391. Chen Z H, He F, Cui Y, et al, 2018. Advances in research of influence on climate of the group of wind farms in past 20 years[J]. Climate Change Res, 14(4): 381-391 (in Chinese).
|
程雪玲, 贺园园, 朱蓉, 2023. 复杂地形对风速廓线的影响[J]. 大气科学. Cheng X L, He Y Y, Zhu R, 2023. Effect of complex terrain on the wind speed profile[J]. Chin J Atmos Sci (in Chinese). DOI:10.3878/j.jssn.1006-9895.2301.22181
|
杜燕军, 冯长青, 2010. 风切变指数在风电场风资源评估中的应用[J]. 电网与清洁能源, 26(5): 62-66. Du Y J, Feng C Q, 2010. Application of wind shear index in the assessment of wind resources of wind farm[J]. Power Syst Clean Energy, 26(5): 62-66 (in Chinese).
|
龚强, 汪宏宇, 朱玲, 等, 2015. 辽宁省近地层风切变特征研究[J]. 自然资源学报, 30(9): 1560-1569. Gong Q, Wang H Y, Zhu L, et al, 2015. Study on the near surface wind shear characteristics in Liaoning Province[J]. J Natural Resources, 30(9): 1560-1569 (in Chinese).
|
龚玺, 朱蓉, 李泽椿, 2018. 我国不同下垫面的近地层风切变指数研究[J]. 气象, 44(9): 1160-1168. Gong X, Zhu R, Li Z C, 2018. Study of near surface wind shear exponents of different regions in China[J]. Meteor Mon, 44(9): 1160-1168 (in Chinese).
|
贺德馨, 2006. 风工程与工业空气动力学[M]. 北京: 国防工业出版社: 29-48. He D X, 2006. Wind Engineering and Industrial Aerodynamics[M].
Beijing: National Defense Industry Press: 29-48 (in Chinese).
|
胡菊, 2012. 大型风电场建设对区域气候影响的数值模拟研究[D]. 兰州: 兰州大学. Hu J, 2012. Numerical simulation research on impact of large-scale wind farms on regional climate[D]. Lanzhou: Lanzhou University(in Chinese).
|
黄林宏, 宋丽莉, 李刚, 等, 2016. IEC的风机选型参数估算方法在中国的适用性研究[J]. 气象, 42(12): 1522-1530. Huang L H, Song L L, Li G, et al, 2016. Applicability of estimation methods of IEC WTGS selecting parameters in China[J]. Meteor Mon, 42(12): 1522-1530 (in Chinese).
|
蒋俊霞, 杨丽薇, 李振朝, 等, 2019. 风电场对气候环境的影响研究进展[J]. 地球科学进展, 34(10): 1038-1049. Jiang J X, Yang L W, Li Z C, et al, 2019. Progress in the research on the impact of wind farms on climate and environment[J]. Adv Earth Sci, 34(10): 1038-1049 (in Chinese).
|
靳双龙, 刘晓琳, 李鸿秀, 等, 2019. 不同下垫面垂直风切变特征对比[J]. 兰州大学学报(自然科学版), 55(3): 340-346. Jin S L, Liu X L, Li H X, et al, 2019. A contrast study on vertical wind shear characteristics of different underlying surfaces[J]. J Lanzhou Univ(Nat Sci), 55(3): 340-346 (in Chinese).
|
雷杨娜, 孙娴, 张侠, 等, 2015. 陕北黄土高原复杂地形风速廓线特性[J]. 可再生能源, 33(9): 1345-1350. Lei Y N, Sun X, Zhang X, et al, 2015. Speed profile characteristics of wind over complex terrain in loess plateau in Northern Shaanxi[J]. Renewable Energy Resour, 33(9): 1345-1350 (in Chinese).
|
李国庆, 李晓兵, 2016. 风电场对环境的影响研究进展[J]. 地理科学进展, 35(8): 1017-1026. Li G Q, Li X B, 2016. Research progress of wind farm impact on the environment[J]. Prog Geogr, 35(8): 1017-1026 (in Chinese).
|
李鸿秀, 朱瑞兆, 王蕊, 等, 2014. 不同地形风电场湍流强度日变化和年变化分析[J]. 太阳能学报, 35(11): 2327-2333. Li H X, Zhu R Z, Wang R, et al, 2014. Wind turbulence intensity of daily and yearly variation analysis in different kinds of terrain[J]. Acta Energ Sol Sin, 35(11): 2327-2333 (in Chinese).
|
李鹏, 田景奎, 2011. 不同下垫面近地层风速廓线特征[J]. 资源科学, 33(10): 2005-2010. Li P, Tian J K, 2011. Characteristics of surface layer wind speed profiles over different underlying surfaces[J]. Resour Sci, 33(10): 2005-2010 (in Chinese).
|
梁红, 魏科, 马骄, 2021. 我国西北大规模太阳能与风能发电场建设产生的可能气候效应[J]. 气候与环境研究, 26(2): 123-141. Liang H, Wei K, Ma J, 2021. Climate effect assessment of ideal large-scale solar and wind power farms in Northwest China[J]. Climatic Environ Res, 26(2): 123-141 (in Chinese).
|
刘敬乐, 姚青, 蔡子颖, 等, 2020. 基于天津255 m气象塔对近年天津近地面风和气温变化特征的研究[J]. 气象, 46(9): 1235-1244. Liu J L, Yao Q, Cai Z Y, et al, 2020. Analysis of temperature and wind variation characteristics based on the 255 m meteorological tower in Tianjin[J]. Meteor Mon, 46(9): 1235-1244 (in Chinese).
|
刘磊, 高晓清, 陈伯龙, 等, 2012. 大规模风电场建成后对风能资源影响的研究[J]. 高原气象, 31(4): 1139-1144. Liu L, Gao X Q, Chen B L, et al, 2012. Preliminary estimates of wind energy resources deficit in large wind farm[J]. Plateau Meteor, 31(4): 1139-1144 (in Chinese).
|
刘维维, 2015. 风力机运行对大气边界层影响的数值模拟研究[D]. 保定: 华北电力大学. Liu W W, 2015. Numerical simulation of the influence of the wind turbine on the atmosphere boundary layer[D]. Baoding: North China Electric Power University(in Chinese).
|
龙强, 王锋, 王畅, 等, 2020. 渤海湾北岸海陆风及湍流强度特征分析[J]. 应用海洋学学报, 39(3): 303-311. Long Q, Wang F, Wang C, et al, 2020. Characteristics of sea-land breeze and turbulence intensity on the north shore of Bohai Bay[J]. J Appl Oceanogr, 39(3): 303-311 (in Chinese).
|
路婷婷, 2017. 基于大涡模拟的风力机对区域大气边界层特性影响研究[D]. 保定: 华北电力大学. Lu T T, 2017. Study on the influence of wind turbine on the characteristics of regional atmospheric boundary layer based on large eddy simulation[D]. Bao-ding: North China Electric Power University(in Chinese).
|
马晓梅, 王博, 刘永前, 2020. 一种基于大气稳定度的风资源评估方法[J]. 可再生能源, 38(1): 47-52. Ma X M, Wang B, Liu Y Q, 2020. A wind resource assessment method based on atmospheric stability[J]. Renew Energy Resour, 38(1): 47-52 (in Chinese).
|
彭怀午, 冯长青, 包紫光, 2010. 风资源评价中风切变指数的研究[J]. 可再生能源, 28(1): 21-23, 28. Peng H W, Feng C Q, Bao Z G, 2010. Study on the wind shear exponent for wind resource assessment[J]. Renewable Energy Resour, 28(1): 21-23, 28 (in Chinese).
|
史军, 穆海振, 徐家良, 等, 2016. 上海近海海域低层风特性分析[J]. 长江流域资源与环境, 25(8): 1184-1190. Shi J, Mu H Z, Xu J L, et al, 2016. Characteristics of low-level wind in Shanghai offshore[J]. Resour Environ Yangtze Basin, 25(8): 1184-1190 (in Chinese).
|
孙淑芳, 张广兴, 2019. 新疆百里风区近地层垂直风切变指数特征[J]. 沙漠与绿洲气象, 13(3): 99-104. Sun S F, Zhang G X, 2019. Characteristics of vertical wind shear exponent in the atmospheric surface layer in the 100-kilometer wind area of Xinjiang[J]. Desert Oasis Meteor, 13(3): 99-104 (in Chinese).
|
王承煦, 张源, 2003. 风力发电[M]. 北京: 中国电力出版社. Wang C X, Zhang Y, 2003. Wind Power Generation[M].
Beijing: China Electric Power Press (in Chinese).
|
吴婷婷, 2015. 风电场风能资源评估中重要参数的计算与应用[D]. 呼和浩特: 内蒙古工业大学. Wu T T, 2015. Research on important parameters calculation and application of wind energy resource assessment for wind farm[D]. Huhhot: Inner Mongolia University of Technology(in Chinese).
|
吴正人, 刘维维, 王松岭, 2014. 风力发电对局地气候的潜在影响分析[J]. 中国电力, 47(6): 101-105. Wu Z R, Liu W W, Wang S L, 2014. Analysis of potential impacts of wind power generation on local climate[J]. Electric Power, 47(6): 101-105 (in Chinese).
|
夏馨, 余晔, 董龙翔, 等, 2022. 风电场建设前后近地面湍流强度变化特征[J]. 高原气象, 41(4): 1062-1073. Xia X, Yu Y, Dong L X, et al, 2022. Characteristics of near surface turbulence intensity before and after wind farm construction[J]. Plateau Meteor, 41(4): 1062-1073 (in Chinese).
|
徐宝清, 吴婷婷, 李文慧, 2014. 风能风切变指数计算方法的比选研究[J]. 农业工程学报, 30(16): 188-194. Xu B Q, Wu T T, Li W H, 2014. Screening of calculation methods for wind shear exponent[J]. Trans Chin Soc Agric Eng, 30(16): 188-194 (in Chinese).
|
徐荣会, 2014. 干旱区风电场对局地微气象环境的影响研究——以苏尼特右旗朱日和风电场为例[D]. 呼和浩特: 内蒙古农业大学. Xu R H, 2014. The impact of wind farm on local micro meteorological environment in arid region—taking Zhurihe wind farm as an example[D]. Huhhot: Inner Mongolia University of Technology(in Chinese).
|
薛桁, 朱瑞兆, 杨振斌, 等, 2001. 中国风能资源贮量估算[J]. 太阳能学报, 22(2): 167-170. Xue H, Zhu R Z, Yang Z B, et al, 2001. Assessment of wind energy reserves in China[J]. Acta Energiae Sol Sin, 22(2): 167-170 (in Chinese).
|
叶小岭, 支兴亮, 邓华, 2019. 基于WRF模式的风电场短期风速集成预报方法研究[J]. 气象, 45(1): 88-98. Ye X L, Zhi X L, Deng H, 2019. Integrated forecasting method research of short-term wind speed in wind power plant based on WRF model[J]. Meteor Mon, 45(1): 88-98 (in Chinese).
|
张雪婷, 李金鑫, 陈正洪, 等, 2019. 我国大别山区风电场群对局地气候影响研究初探——以湖北大悟为例[J]. 气象科技进展, 9(2): 56-61. Zhang X T, Li J X, Chen Z H, et al, 2019. A preliminary exploration on the influence of wind farms on local climate in Dabie Mountainous areas in China: a case study of Dawu in Hubei Province[J]. Adv Meteor Sci Technol, 9(2): 56-61 (in Chinese).
|
张颖超, 肖寅, 邓华, 2016. 基于ELM的风电场短期风速订正技术研究[J]. 气象, 42(4): 466-471. Zhang Y C, Xiao Y, Deng H, 2016. Modification technology research of short term wind speed in wind farm based on ELM method[J]. Meteor Mon, 42(4): 466-471 (in Chinese).
|
郑新倩, 杨帆, 李超凡, 等, 2019. 巴丹吉林沙漠北缘拐子湖流沙下垫面近地层湍流强度和陆面过程特征[J]. 中国沙漠, 39(4): 103-112. Zheng X Q, Yang F, Li C F, et al, 2019. The turbulence intensity of surface layer and land surface processes over Guaizi Lake shifting sandy land on the northern margin of Badain Jaran Desert[J]. J Desert Res, 39(4): 103-112 (in Chinese).
|
中国气象局, 2007. 风电场气象观测及资料审核、订正技术规范: QX/T 74—2007[S]. 北京: 中国标准出版社. China Meteorological Administration, 2007. Regulations for data inspection and correction of wind power plant meteorological observation: QX/T 74-2007[S]. Beijing: Standards Press of China(in Chinese).
|
中华人民共和国国家质量监督检验检疫总局, 2004. 风电场风能资源评估方法: GB/T 18710—2002[S]. 北京: 中国标准出版社. General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, 2004. Methodology of wind energy resource assessment for wind farm: GB/T 18710-2002[S]. Beijing: Standards Press of China(in Chinese).
|
Buckley C, Scott N, Snodin S, et al, 2005. Review of Impacts of High Wind Penetration in Electricity Networks[M]. Australia: Garrad Hassan Pacific Pty Ltd.
|
Christiansen M B, Hasager C B, 2006. Using airborne and satellite SAR for wake mapping offshore[J]. Wind Energy, 9(5): 437-455. DOI:10.1002/we.196
|
Dincer I, Acar C, 2015. A review on clean energy solutions for better sustainability[J]. Int J Energy Res, 39(5): 585-606. DOI:10.1002/er.3329
|
Farrugia R N, 2003. The wind shear exponent in a Mediterranean island climate[J]. Renewable Energy, 28(4): 647-653. DOI:10.1016/S0960-1481(02)00066-6
|
Fitch A C, 2015. Climate impacts of large-scale wind farms as Parameterized in a global climate model[J]. J Climate, 28(15): 6160-6180. DOI:10.1175/JCLI-D-14-00245.1
|
Fitch A C, Olson J B, Lundquist J K, et al, 2012. Local and mesoscale impacts of wind farms as parameterized in a mesoscale NWP model[J]. Mon Wea Rev, 140(9): 3017-3038. DOI:10.1175/MWR-D-11-00352.1
|
Frandsen S T, Jørgensen H E, Barthelmie R, et al, 2009. The making of a second-generation wind farm efficiency model complex[J]. Wind Energy, 12(5): 445-458. DOI:10.1002/we.351
|
Garrigle E V M, Leahy P G, 2015. Quantifying the value of improved wind energy forecasts in a pool-based electricity market[J]. Renewable Energy, 80: 517-524. DOI:10.1016/j.renene.2015.02.023
|
Li Y, Kalnay E, Motesharrei S, et al, 2018. Climate model shows large-scale wind and solar farms in the Sahara increase rain and vegetation[J]. Science, 361(6406): 1019-1022. DOI:10.1126/science.aar5629
|
Liu Y H, Dang B, Xu Y M, et al, 2021. An observational study on the local climate effect of the Shangyi wind farm in Hebei Province[J]. Adv Atmos Sci, 38(11): 1905-1919. DOI:10.1007/s00376-021-0290-0
|
Phillips J, 2015. A quantitative-based evaluation of the environmental impact and sustainability of a proposed onshore wind farm in the United Kingdom[J]. Renewable Sustainable Energy Rev, 49: 1261-1270. DOI:10.1016/j.rser.2015.04.179
|
Rehman S, Al-Abbadi N M, 2007. Wind shear coefficients and energy yield for Dhahran, Saudi Arabia[J]. Renewable Energy, 32(5): 738-749. DOI:10.1016/j.renene.2006.03.014
|
Smith C M, Barthelmie R J, Pryor S C, 2013. In situ observations of the influence of a large onshore wind farm on near-surface temperature, turbulence intensity and wind speed profiles[J]. Environ Res Lett, 8(3): 034006. DOI:10.1088/1748-9326/8/3/034006
|