黨 冰 劉勇洪 杜吳鵬 陳可仁

1 北京市氣候中心,北京 100089 2 中國氣象局地球系統數值預報中心,北京 100081 3 災害天氣國家重點實驗室,北京 100081 4 龍源(北京)風電工程設計咨詢有限公司,北京 100034

提 要: 風電場存在著明顯的局地氣候效應,但陸/海兩種不同類型風電場是否存在局地風環(huán)境效應的明顯差異?以河北尚義陸地風電場和江蘇如東海上風電場為例,利用風電場周邊氣象站及測風塔典型年份觀測資料,開展了陸/海風電場對湍流強度(TI)和風切變指數(WSE)的影響差異初步分析研究。研究結果表明:陸/海風電場對TI 和WSE會產生顯著影響;其中,陸/海風電場對TI均為增強效應,建設后年均TI分別增加31%和37%,最大增幅分別發(fā)生在春季(47%)和冬季(49%);影響差異主要為陸地風電場TI增幅在高層明顯大于低層,夜間大于白天,而海上風電場不同高度TI增幅及日變化則較為平穩(wěn);陸/海風電場對WSE影響差異顯著,陸地風電場建設后WSE白天增加、夜晚降低,日變化明顯減小,年均降低8%,最大降幅發(fā)生在秋季(12%);海上風電場建設后WSE白天、夜晚均明顯增加,年均增加24%,最大增幅發(fā)生在春季(37%)。

引 言

由于化石能源的不可再生性,以及其在燃燒過程中排放的溫室氣體和污染氣體帶來的氣候環(huán)境問題日益嚴重,開發(fā)、利用可再生能源已成為世界各國急需解決的問題(Dincer and Acar,2015;李國慶和李曉兵,2016;陳正洪等,2018)。風能作為一種潔凈的、儲量極為豐富的可再生能源,越來越受到廣泛重視(Garrigle and Leahy,2015;Phillips,2015;葉小嶺等,2019)。我國風能資源儲量豐富,總儲量為322 600 GW,實際可開發(fā)量為25 300 GW(薛桁等,2001;梁紅等,2021)。近年來我國風電發(fā)展迅速,據全球風能協會公布的數據顯示,2021年中國風力發(fā)電裝機總量達338.31 GW,占世界總裝機容量的40.4%;其中海上風電發(fā)展迅猛,2021年新增裝機量16.9 GW,約占全球新增裝機量80%。由于風能的能量密度小,發(fā)展風能需要較大的集能面積,因此大規(guī)模建設風力發(fā)電場會改變地表粗糙度,從而改變陸表和大氣的熱交換過程;同時由于風機渦輪振動、風能轉化為電能,改變了自然界原有的能源循環(huán)模式,導致風電場對局地或區(qū)域氣候可能產生影響(Buckley et al,2005;張雪婷等,2019;蔣俊霞等,2019)。由此,大型風電場建設和運行對氣候和環(huán)境的影響受到越來越多的關注(Fitch,2015;Li et al,2018;梁紅等,2021)。

關于大型風電場建設運行對局地風環(huán)境的影響已有大量研究。Christiansen and Hasager(2006)利用SAR雷達數據發(fā)現丹麥地區(qū)海上風電場的平均風速減小8%～9%,對下風向的影響范圍達20 km;Fitch et al(2012)通過數值模擬發(fā)現從海上風電場延伸至其60 km下游區(qū)域均有風速衰減現象,風電場區(qū)域內風速衰減最多可達16%;Smith et al(2013)對美國中西部大型陸上風電場觀測研究顯示發(fā)現在葉片掃過區(qū)域及風機輪軸高度處風速明顯降低、湍流強度顯著增強,且夜間更為明顯。胡菊(2012)模擬發(fā)現我國河西走廊地區(qū)的大型風電場建成后平均風速減小了0.3 m·s-1;徐榮會(2014)根據觀測發(fā)現內蒙古地區(qū)風電場局地風速在風機的擾動下有減小趨勢;張雪婷等(2019)發(fā)現我國大別山區(qū)風電場群建設和運行使大悟縣局地風速略有減小。但國內研究主要集中在陸上風電場,對海上風電場研究較少,且關于風電場建設對一些重要指標如湍流強度和風切變指數的影響研究也很少。

湍流強度(TI)和風切變指數(WSE)是評價氣流穩(wěn)定程度的指標,對于風電場的風能評估及風電場穩(wěn)定運行具有重要意義(李鴻秀等,2014;黃林宏等,2016;龔璽等,2018;陳燕等,2019;靳雙龍等,2019)。TI不但與地表摩擦和風速切變引起的動力因子以及溫度層結引起的熱力因子密切相關,還與地理位置、地形和地表粗糙度等因素有關(賀德馨,2006;李鴻秀等,2014;劉敬樂等,2020);WSE表征風速隨高度的變化程度,主要受地表粗糙度、大氣穩(wěn)定度等因素影響(徐寶清等,2014;龔強等,2015;馬曉梅等,2020;常蕊等,2022;程雪玲等,2023)。目前利用測風塔資料開展了風電場建設前風環(huán)境特征如風速、TI或WSE等分析(李鵬和田景奎,2011;班欣等,2012),或者是利用數值模擬技術開展風電場對大氣邊界層特征影響研究(劉維維,2015;路婷婷,2017;張穎超等,2016;常蕊等,2013),但對風電場群建設后局地風速以及TI和WSE變化的觀測研究較少。由于風電場建設后,局地風速和下墊面粗糙度發(fā)生變化,且風機的擾動改變了局地大氣的穩(wěn)定性,相應改變TI和WSE(吳正人等,2014;靳雙龍等,2019;夏馨等,2022),從而改變了風電場建設前的風能評估結果以及風電場安全運行所需的風速條件(劉磊等,2012);同時由于陸/海風電場存在著下墊面及周圍環(huán)境的巨大差異,因此有必要開展二者在風電場建設后對局地風環(huán)境關鍵參數影響及對比研究,為進一步闡明陸/海風電場對局地風環(huán)境的不同影響機理提供參考。

由此,本文分別選擇河北尚義陸地風電場和江蘇如東海上風電場,利用風電場群周邊國家氣象站風速觀測資料以及質量控制后的測風塔測風資料,選定風電場建設前后典型時段,通過分析局地風速、TI和WSE的變化,開展陸/海風電場群建設對局地風場影響差異初步研究。一方面闡明典型陸/海下墊面風電場群建設后局地風環(huán)境關鍵參數的變化規(guī)律,另一方面揭示這兩種不同下墊面風電場群對風環(huán)境關鍵參數的影響差異機制,從而為今后陸/海風能資源評估、有效規(guī)避相關風險等提供科學依據。

1 研究區(qū)與數據

1.1 研究區(qū)概況

本文研究選取的陸地風電場位于中國河北張家口尚義縣(如圖1a所示),尚義風電場群從2008—2015年底陸續(xù)建成224臺風機(表1),風機輪轂高度為65～90 m,總裝機容量為450 MW;風電場所在區(qū)域屬于半干旱大陸性季風氣候區(qū),海拔高度為1375～1573 m,自然植被以草本植物為主,山谷及緩坡地帶有零星小面積成塊林地、耕地及灌木;氣候特點表現為氣溫低溫差大、雨量少而集中、光照充足,年均降水量為414 mm,年均風速為3.4 m·s-1,全年主導風為西北向(NW),次主導風為西西南(WWS)。

表1 陸/海風電場群建設概況一覽表Table 1 The general construction situation of inland/offshore wind farms

圖1 (a)尚義風電場群和(b)如東風電場群的位置、地形及周邊氣象站點分布Fig.1 The geographical location of inland/offshore wind farms and meteorological observation stations (a) Shangyi Wind Farms, (b) Rudong Wind Farms

海上風電場研究區(qū)位于中國江蘇南通如東縣(如圖1b所示),如東風電場群從2009—2015年陸續(xù)建成155臺風機,風機輪轂高度為65～90 m,總裝機容量為482 MW;風電場所在區(qū)域屬于亞熱帶海洋性季風氣候區(qū),境內地勢平坦,海拔高度為3～5 m,陸地地貌是典型的濱海平原,沿海屬黃海南部海域,沿海潮汐24小時漲落兩次;受海洋調節(jié)和季風環(huán)流影響,四季分明、氣候溫和、光照充足、雨水充沛,年降水量為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)以及氣象行業(yè)標準《風電場氣象觀測及資料審核、訂正技術規(guī)范(QX/T 74-2007)》(中國氣象局,2007)進行質量控制后的氣象觀測資料,選擇合適的陸/海風電場建設前后分析時段(表2)。尚義風電場選用建設前后各3年(2005—2007年和2016—2018年)的國家氣象站資料進行局地風速變化分析,而0704#測風塔資料由于存在不連續(xù)和缺測情況,選擇數據較為完整的2006年和2017年進行TI和WSE變化分析。如東風電場選用0865#測風塔2009年和2018年數據進行TI和WSE變化分析。

表2 陸/海風電場群建設前后選擇的分析時段表Table 2 The analysis period of inland/offshore wind farms before and after construction

2 研究方法

2.1 評估站選擇

2.1.1 陸上風電場影響站與背景站選擇

河北尚義陸地風電場群周圍60 km內包含4個 國家氣象站:尚義、商都、懷安和興和站(圖1a),其中尚義站緊鄰風電場群,可作為風電場影響站進行評估;此外,尚義0704#測風塔位于風電場群內,也可作為風電場影響站進行評估。

利用尚義站評估尚義陸地風電場局地風速變化時,除了風電場群的可能影響,氣象站氣象要素變化還受年際氣象條件波動和氣候變化影響,僅靠單一氣象站是不能識別出風電場對其單獨影響,需要借助研究站與背景站(風電場非影響區(qū))之間的對比分析來消除年際波動與氣候變化影響。由此,選擇合適的背景站是尚義陸地風電場局地風速變化評估的一個關鍵步驟。

背景站一般要求與研究站的氣候類型和地理條件相似,但又不能受風電場影響;數值模擬研究表明風電場影響風速衰減的距離在風電場下游方向可以達到30～60 km(Frandsen et al, 2009),要求背景站距風電場至少30 km。這里分別對風電場建設前(1981—2007年)的尚義站與懷安、興和、商都站的主要氣候要素和地理條件進行相似性分析,如表3所示。

表3 國家氣象站(尚義、懷安、商都、興和)氣候和地理條件相似性分析Table 3 The similarity analysis of climate and geographical conditions among four national meteorological stations

由表3分析可知,尚義站年均氣溫、風速和相對濕度均與商都站最高,且兩者海拔差最小(僅44 m),空間距離在30 km以上,且商都站在風電場主導風向(NW)上游區(qū)域,因此可選擇商都站作為背景站進行尚義陸地風電場的局地風速變化分析。

為進一步分析尚義陸上風電場的影響范圍,還需對可能受到風電場影響的其他兩個站(懷安與興和站)進行分析。由于懷安站和背景站海拔相差很大(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年年均風速資料分析發(fā)現其在2014年風速明顯下降且低于如皋站,而在1981—2013年期間如東站風速均明顯高于如皋站;根據文中表1可知,2014年前如東海上風電場已投產94個風機,2014年以后投產61個,因此如東站的風速突然降低很可能是其他外部原因造成,經調查得知2014—2015年如東氣象站周圍建立了大量高樓,致使2014年后風速下降明顯,并非風電場影響造成;而如皋站距離風電場較遠(50 km),且在海上風電場群的西南方向,與風電場群主導風向(東南—偏北風)相差甚大,幾乎不受風電場影響,因此如東站和如皋站均不適宜作為風電場影響站進行評估,即如東海上風電場沒有合適的國家氣象站進行局地風速變化評估。

如東0865#測風塔位于風電場群偏東方向約7 km 處,處于風電場海上影響范圍之內,可作為影響站進行評估。

2.2 局地風速影響評估

定義尚義陸上風電場群對局地風速的影響值為:

dV=VI-VB

(1)

ΔV=dVAF-dVBF

(2)

式中:dV為國家氣象站影響站(或可能影響站)10 m 高度風速(VI)與背景站10 m高度風速(VB)之差(單位:m·s-1);ΔV為風電場建設后(dVAF)相對于風電場建設前(dVBF)的風速變化(單位:m·s-1),即風電場對風速的局地影響。

為評估風電場群建設前后局地風速是否發(fā)生顯著變化,這里以風電場建設前影響站(或可能影響站)與背景站風速差dV的年際標準差(StddV)作為年際波動,如果風速差變化值(ΔV)的絕對值大于StddV,則表明該地風速受到風電場群影響,如果ΔV絕對值小于或等于StddV,則表明該地風速沒有受到風電場群影響。

由于如東海上風電場沒有合適的國家氣象站進行局地風速研究,雖然0865#測風塔可作為影響站進行評估,但由于該塔建設前僅有1年資料,且沒有10 m高度的測風數據,因此本文不開展如東海上風電場的局地風速影響研究。

2.3 局地湍流強度(TI)影響評估

根據國家標準《風電場風能資源評估方法(GBT 18710—2002)》(中華人民共和國國家質量監(jiān)督檢驗檢疫總局,2004),湍流強度的計算公式為:

TI=σ/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)表示風速垂直切變的強度。在這里,采用風速隨高度變化呈指數規(guī)律分布(王承煦和張源,2003;Farrugia,2003;Rehman and Al-Abbadi,2007)進行估算,即:

(4)

式中:Vn、V1分別為Zn、Z1高度上的風速,α為風切變指數WSE。

對式(4)作變換,并兩邊取對數,即得:

(5)

目前計算WSE時,不同風速資料的選擇會造成估算結果偏差明顯。彭懷午等(2010)、吳婷婷(2015)研究表明采用全部數據進行風切變計算產生的平均值代表性差、偏離系數較大,而采用去除3 m·s-1以下風速數據計算產生的平均值代表性要優(yōu)于前者、偏離系數較小。因此,本文利用測風塔資料,選取不同高度去除3 m·s-1以下且高層風速大于低層風速的數據,根據式(5)計算每兩個高度層間的WSE,再利用最小二乘法(徐寶清等,2014;陳燕等,2019;孫淑芳和張廣興,2019)將各層間的WSE進行擬合,得到該區(qū)域平均WSE,即根據式(5)可得到:

ln(Vn/V1)=αln(Zn/Z1)

(6)

本文定義風電場對WSE的局地影響為風電場建設后WSE與建設前WSE的距平百分率,該值為正值時,表明風電場建設運行使得風切變增加;該值為負值時,表明風電場建設運行使得風切變下降。

3 結果分析

3.1 陸地風電場群局地風速變化

3.1.1 年變化趨勢

就影響站尚義站而言,2005—2011年影響站和背景站(商都站)風速變化趨勢具有良好的一致性(圖2a),2011年后,背景站風速(VSD)總體呈上升趨勢,而影響站風速(VSY)則呈下降趨勢,意味著2011年后影響站受到非氣候變化因素導致其風速變化趨勢與背景站差異顯著。從兩站風速差(dV)變化來看,2012年后兩站年風速差呈增大趨勢,其中2016—2018年風速差明顯大于風電場竣工前任何一年(圖2a),該結果顯示尚義站很可能受到風電場群影響,致使風速在風電場建設前后發(fā)生明顯變化。

圖2 2005—2018年(a)尚義站和(b)興和站年平均風速及與商都站風速差變化趨勢Fig.2 The variation trend of average wind speed at (a) Shangyi Station and (b) Xinghe Station and the differences between two stations and Shangdu Station respectively during 2005-2018

就可能影響站(興和站)而言,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波動范圍內,表明風電場建設后風速降低主要發(fā)生在春季和冬季。

圖3 尚義風電場群建設前后(a)尚義站和(b)興和站風速差變化與氣候年際標準差對比Fig.3 The comparison of wind speed difference changes (ΔV) and climate interannual standard deviations (StddV) of (a) Shangyi Station and (b) Xinghe Station before and after the construction of Shangyi wind farms

就可能影響站(XH)而言,年和春、夏、秋季ΔV均在StddV波動范圍內,僅冬季ΔV略高于StddV(偏高0.04 m·s-1),表明風電場建設后興和站風速變化基本處于正常年際波動內。該結果顯示興和站受尚義風電場影響很小,即使受到影響也僅僅發(fā)生在冬季,而且這種影響非常輕微。由于興和站位于風電場偏南方向約29 km處,因此可以推測尚義陸上風電場群影響范圍很可能不超過30 km。

3.2 陸/海風電場群湍流強度變化

3.2.1 日變化

分別選擇尚義陸地測風塔30 m、70 m和如東海上測風塔15 m、80 m高度層,分析建設前后TI日變化(圖4):無論是尚義陸上風電場群,還是如東海上風電場群,建設后各高度層不同時次TI均明顯增加;這主要是由于風電場群內大范圍風機的攪動,使得空氣上下交換加強,致使陸/海風電場區(qū)TI增加明顯(Christiansen and Hasager,2006;Smith et al,2013)。就不同高度來看,低層TI明顯高于高層,這是由于在低層大氣更易受到地表粗糙度與地形影響,大氣上下交換強,因此TI更強。

圖4 (a)尚義和(b)如東風電場群建設前后不同高度湍流強度日變化對比Fig.4 The comparison of diurnal variation of TI at different heights at (a) Shangyi Station and (b) Rudong Station before and after construction of inland/offshore wind farms

尚義陸上風電場TI日變化呈單峰曲線,白天TI明顯大于夜間,這是由于日出后地表受熱,對流逐漸加強,導致大氣層結不穩(wěn)定,TI逐漸增大,至中午前后達到最大;日落后地表受熱停止,對流作用減小,大氣層結逐漸轉為中性,TI也隨之減小。與陸上風電場并不相同,如東海上風電場TI日變化不明顯,這主要是由于海面水體熱容量大,增溫和降溫過程較陸地緩慢,導致海上區(qū)域熱力對流變化小于陸地,因此海上風電場TI日變化波動較小。

對比尚義陸地和如東海上風電場群同一高度層50 m建設前后TI日變化,如圖5a所示:無論建設前后,陸地TI日均值均明顯高于海上,尤其在白天,這是由于海洋熱容量明顯高于陸地,造成海上白天近地層溫度梯度小于陸地,加之下墊面平滑,從而使得同一高度上海洋TI明顯小于陸地;到了夜晚,陸地和海水上大氣層結都比較穩(wěn)定,溫度梯度接近,因此TI比較接近。

圖5 陸/海風電場群(a)建設前后50 m高度湍流強度及(b)建設后不同高度湍流強度增幅的日變化Fig.5 The diurnal variation of (a) the TI at 50 m height before and after construction and (b) TI increasing range at different heights after construction of inland/offshore wind farms

對比陸/海風電場群建設后各典型高度層(尚義:30 m、70 m;如東:15 m、80 m)TI增幅日變化(圖5b),可以看出:就不同時次來看,陸地風電場群各高度層TI較建設前增加11%～66%,其中夜晚增幅(平均為45%)明顯大于白天(平均為20%)。這是因為陸地區(qū)域夜間大氣穩(wěn)定度高,對流本身較弱,但由于風機的攪動,使得上下空氣交換加強,致使湍流強度大幅增加;而白天地面受大氣加熱影響,大氣層結不穩(wěn)定,對流本身較強,雖然風機的攪動能更進一步加強空氣擾動,但相比于夜間增幅小得多。

就不同高度來看,陸地風電場建設后70 m高度TI平均增幅為40%,明顯高于30 m TI增幅(28%),表明TI在高層增幅大于低層,這種增幅在夜間更為明顯。這很可能一方面與70 m高度TI在夜間處于低值區(qū)且明顯小于30 m高度(圖4a)有關;另一方面由于70 m高度層更接近風機輪轂高度(平均75 m),受風機攪動影響更大,而尚義風輪平均直徑約為90 m,表明葉片最底端與地面距離約為30 m,因此30 m高度受風機影響較小。

如東海上風電場群建設后各高度TI增加29%～ 48%,但增幅沒有明顯的日變化特征。這是因為雖然海上區(qū)域風機攪動加強了空氣擾動,使TI增加,但由于TI日變化較小(圖4b),TI增幅的日變化也較小。此外,80 m高度TI平均增幅為39%,稍高于15 m高度平均增幅(36%),表明海上風機對低層和高層TI影響造成的差異遠不如陸地風機的影響明顯。

3.2.2 年和季節(jié)變化

尚義陸地和如東海上風電場群建設前后不同高度年均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來看,無論建設前后,陸地區(qū)域TI均明顯高于海上,偏高約21%,這是由于海面相對于陸地下墊面平滑、地表粗糙度小,加之海水熱容量高于陸地、近海面層溫度梯度小于陸地近地層,導致海上風電場區(qū)域TI相對于陸地風電場整體較小。

圖6 (a)尚義和(b)如東風電場群建設前后不同高度年湍流強度對比Fig.6 Comparison of annual TI at different heights at (a) Shangyi Station and (b) Rudong Station before and after construction of inland/offshore wind farms

以尚義陸地和如東海上風電場群測風塔50 m高度TI四季變化進行對比,如圖7a所示:從季節(jié)變化上看,無論建設前后,陸地風電場TI四季變化均較為明顯(建設前0.091～0.126,建設后0.121～0.155),其中夏季最大、冬季最小,春秋季介于之間;海上風電場TI季節(jié)變化建設前較為平穩(wěn)(0.080～0.096),季節(jié)變化幅度明顯小于陸地,這與已有的風電場建設前不同地形TI變化相關研究結果較為一致(李鴻秀等,2014;龍強等,2020);而風電場群建設后如東海上風電場TI季節(jié)變化加大(0.098～0.135),其中冬秋TI強于春夏,這與陸地風電場季節(jié)變化規(guī)律并不一致。

圖7 陸/海風電場群50 m高度湍流強度(a)建設前后四季變化及(b)建設后湍流強度增幅對比Fig.7 (a) The seasonal variation of TI at the 50 m height before and after construction and (b) comparison of TI increasing range after construction in different seasons

從季節(jié)TI增幅看(圖7b),尚義陸地風電場春季TI增幅(47%)最為明顯,然后依次為冬季(33%)、秋季(33%)和夏季(24%),這與前面分析的春季和冬季風速受風電場影響顯著具有一致性。而如東海上風電場冬季TI增幅最大(49%),然后依次為秋季(42%)、春季(34%)和夏季(21%)。

3.3 陸/海風電場群風切變指數變化

3.3.1 日變化

尚義陸地和如東海上風電場群建設前后,年均WSE日變化如圖8a所示:兩個風電場群WSE均表現為夜間高、白天低等日變化特征,但陸地風電場WSE日變化明顯大于海上風電場。這是由于白天對流作用強、大氣層結不穩(wěn)定,湍流活動旺盛,使不同大氣垂直層間的動量交換頻繁,風速垂直梯度較小,導致WSE較低;夜間地表溫度低,大氣層趨于穩(wěn)定,層間很少混合或不混合,風速隨高度的上升急劇增大,導致WSE較高,這與雷楊娜等(2015)相關研究結果較為一致。

圖8 陸/海風電場群(a)建設前后風切變指數及(b)建設后風切變指數變幅日變化Fig.8 (a) Diurnal variation of WSE before and after construction and (b) variable amplitude of WSE after the construction of inland/offshore wind farms

對比陸/海風電場群建設后WSE變幅的日變化分布(圖8b),可以看出:風電場建設后,陸地風電場群WSE在白天明顯增加(平均增幅為33%),夜間則降低(平均降幅為9%),WSE變幅趨勢與日WSE大小分布正相反,WSE小的時段增幅明顯,WSE大的時段降幅較大,從而使得WSE的日變化值變小(從建設前的0.199降至建設后的0.133)。這是由于陸地區(qū)域WSE日變化較為明顯,白天大氣不穩(wěn)定、風速垂直梯度較小,但由于高層風速更接近風機輪轂高度受風機擾動影響更大,風速增大幅度大于低層,使得上下層風速梯度增大,致使白天WSE大幅增加;夜間大氣穩(wěn)定、風速垂直梯度較大,風機擾動使得上下空氣交換加強、風速梯度變小,致使夜晚WSE降低。海上風電場群建設后,WSE在白天夜晚均有所增加(平均增幅為24%),但由于WSE日變化較小,WSE增幅無明顯日變化趨勢,且建設前后WSE日變化值也沒有明顯變化(分別為0.040和0.039)。

3.3.2 年和季變化

尚義陸地和如東海上風電場群建設前后,年和四季平均WSE及建設后變幅如圖9所示。

圖9 陸/海風電場群(a)建設前后風切變指數及(b)建設后風切變指數變幅對比Fag.9 (a) Annual and seasonal variation of WSE before and after the construction and (b) variable amplitude of WSE after the construction of inland/offshore wind farms

風電場建設前,陸地區(qū)域年均WSE為0.196(圖9a),明顯高于海上區(qū)域(年均WSE為0.133)。這是由于尚義處于內陸山區(qū),地形起伏和地表粗糙度都較大,使得風速隨高度變化大,因此WSE也相對較大;如東沿海地形平坦,氣流受地形和粗糙度影響較少,風速在不同高度的分布較穩(wěn)定,故WSE相對較小。陸地風電場WSE表現為冬季和秋季較大,春季和夏季較小,這是由于冬季氣溫低,太陽輻射強度較低,大氣層結穩(wěn)定度高,動量、熱量交換較少,高低層之間風速差異大,從而WSE大;夏季氣溫高,高強度的太陽輻射引起大氣層的不穩(wěn)定性增強,大氣層間的混合加強,高低層之間風速差異變小,導致WSE較小(杜燕軍和馮長青,2010);海上風電場WSE表現為冬季和春季較大,夏季和秋季較小,這很可能與如東海上風電場夏秋溫度高于冬春、大氣混合較冬春更加明顯有關。

尚義陸地風電場建設后,無論年和季節(jié)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年逐日觀測資料、區(qū)域內測風塔2006—2017年不同高度逐時測風資料以及江蘇如東海上風電場群區(qū)域測風塔2009—2018年不同高度逐時測風資料,根據風電場群建設和投產時間以及質控后可用的測風塔數據,選定風電場建設前后分析時段,通過計算建設前后局地風速差值、湍流強度TI和風切變指數WSE,對比研究陸/海風電場群建設對局地風環(huán)境關鍵參數的影響,初步得到以下主要結論。

(1)尚義陸地風電場具有明顯風速減弱效應,年均風速較建設前降低0.43 m·s-1,且風速降低主要發(fā)生在春季和冬季,風速分別降低了0.89 m·s-1和0.71 m·s-1,風電場影響范圍很可能不超過30 km。

(2)陸地風電場區(qū)TI值高于海上風電場區(qū),TI日變化和季節(jié)變化也比海上區(qū)域明顯。陸/海風電場群建設對TI均具有增強效應,建設后年均TI分別增加了31%和37%。陸地風電場區(qū)TI夜間增幅(45%)明顯大于白天增幅(20%),且高層增幅(40%)大于低層(28%),而海上風電場區(qū)TI不同高度增幅及其日變化較為平穩(wěn)。陸地風電場區(qū)TI四季增幅為24%～47%;海上風電場區(qū)TI四季增幅為21%～49%。

(3)陸/海風電場區(qū)域WSE均呈現冬季高、夏季低、夜間大于白天等特征,但陸地WSE日變化明顯大于海上。陸/海風電場建設后對WSE影響存在顯著差異,陸地風電場建設后WSE白天增加了33%、夜晚則降低了9%,從而導致WSE日變化變小,且年和四季WSE均有所降低,其中年均降低了8%,四季降低了3%～12%;海上風電場建設后WSE白天、夜晚均明顯增加,從而導致年和四季WSE均明顯增加,其中年均增加了24%,四季增加了4%～37%。

4.2 討 論

(1)風電場群對湍流強度TI影響的不確定性。本文所得結果顯示陸/海風電場群在建設前湍流強度TI年、季、日變化均存在明顯差異,這也與已有風電場建設前相關研究結果(史軍等,2016;鄭新倩等,2019)較為一致,表明本文的研究資料具有一定代表性,方法較為合理。而風電場群建設后,陸上和海上TI均發(fā)生了較為明顯的變化,凸顯了風電場群對海/陸局地風場的影響,由于海陸兩種下墊面性質的不同,這種影響存在著不同高度、白天/夜晚和季節(jié)差異,但這種差異是否受到其他條件例如年際變化或者背景氣候變化的影響還并不清楚,因此未來還需更長時間的氣象觀測資料來開展風電場群對TI長期的影響研究。

(2)風電場群對風切變指數WSE影響的不確定性。本文結果顯示風電場建設前陸/海WSE呈現冬季較高、夏季較低、白天較小、夜間較高等特點,與已有風電場建設前的相關研究結果較為一致(彭懷午等,2010;孫淑芳和張廣興,2019;陳燕等,2019);風電場建設后由于風機一定程度上改變了不同高度的風速使得WSE發(fā)生變化,但這種變化實際上是很復雜的,因為WSE不僅受地表粗糙度影響,還受地形、大氣穩(wěn)定度、海拔高度以及熱力層結的影響(李鵬和田景奎,2011),而風電場不僅影響地表粗糙度,還影響風速和大氣溫度如局地增溫降風效應(張雪婷等,2019;Liu et al,2021)。因此風電場群對WSE的影響機制未來還需要更多觀測資料進行佐證。

(3)資料分析的不確定性。由于測風塔資料限制,本文僅選取了風電場群建設前后各1年資料,而且由于陸/海測風塔資料時間差異,選取的年份也不同,因此存在資料代表性不夠充分造成的不確定性問題,從而造成陸/海風電場群對局地風環(huán)境影響差異分析存在一定不確定性。此外,不同測風塔資料實際上是受氣候背景影響并存在年際變化的,雖然TI和WSE是關于風速的相對變化指標,很大程度上已經消除了氣候背景場(含年際變化)的影響,但由于沒有周圍參考塔同高度風速資料,陸/海風電場對TI和WSE的分析與對比仍存在一定程度受背景場影響的不確定性。另外,如東測風塔位于海上風電場東南方向,考慮盛行風向的季節(jié)變換(冬季盛行風為偏北風,夏季盛行風為偏南風),測風塔受風電場影響在冬季較夏季更為明顯,這可能也是如東海上風電場對TI和WSE的影響在冬季最大、夏季最小的另外一個原因。因此,未來還需選取更多的站點以及更有代表性的資料進行詳細分析,才能消除這種由于資料不足帶來的不確定性。