包道日娜, 劉嘉文, 尚 偉, 王帥龍, 劉旭江
(內(nèi)蒙古工業(yè)大學(xué) 能源與動(dòng)力工程學(xué)院, 呼和浩特 010051)
目前,我國(guó)風(fēng)力發(fā)電技術(shù)已經(jīng)較為發(fā)達(dá),裝機(jī)總?cè)萘空紦?jù)世界市場(chǎng)份額的33.6%,早在2010年就位列世界第一,超越了整個(gè)歐盟,并且我國(guó)風(fēng)力發(fā)電裝機(jī)量仍在持續(xù)上漲[1-2]。風(fēng)力發(fā)電技術(shù)的長(zhǎng)足進(jìn)步離不開(kāi)風(fēng)電測(cè)試認(rèn)證技術(shù)的支持,但隨著風(fēng)力發(fā)電技術(shù)的日益提高,我國(guó)風(fēng)電測(cè)試認(rèn)證技術(shù)的發(fā)展顯得尤為不足[3-4]。在風(fēng)電測(cè)試認(rèn)證技術(shù)領(lǐng)域,除了風(fēng)電測(cè)試場(chǎng)的微觀選址,最重要的就是風(fēng)力機(jī)的機(jī)位布置。良好的布置方案不但可以滿(mǎn)足測(cè)試標(biāo)準(zhǔn)、縮小測(cè)試周期,還可以提高測(cè)試場(chǎng)的年發(fā)電量和風(fēng)能利用率。
在風(fēng)電測(cè)試場(chǎng)實(shí)際運(yùn)行過(guò)程中,因不同機(jī)位布置方案引發(fā)的尾流損失是造成測(cè)試場(chǎng)測(cè)試周期較長(zhǎng)和發(fā)電能力不足的一大因素,特別是在風(fēng)資源相對(duì)匱乏的區(qū)域,尾流損失對(duì)風(fēng)電測(cè)試場(chǎng)效益的影響更為嚴(yán)重。因此,通過(guò)優(yōu)化風(fēng)電測(cè)試場(chǎng)機(jī)位布置方案以減少尾流損失從而提高測(cè)試場(chǎng)效益引發(fā)了越來(lái)越廣泛的關(guān)注。國(guó)內(nèi)外學(xué)者在風(fēng)電場(chǎng)風(fēng)力機(jī)機(jī)位布置方面做了大量研究,取得了卓有成效的研究成果。文獻(xiàn)[5]提出了一種通過(guò)優(yōu)化風(fēng)電場(chǎng)主風(fēng)向上風(fēng)機(jī)塔筒高度從而降低風(fēng)電場(chǎng)尾流損失的方法,并用Wasp10.0軟件驗(yàn)證了方法的有效性。文獻(xiàn)[6]基于尾流損失影響,利用Wind PRO軟件進(jìn)行了風(fēng)電機(jī)組排布優(yōu)化,最終得到了一種尾流影響較小的風(fēng)電場(chǎng)風(fēng)力機(jī)排布方案。文獻(xiàn)[7]以高斯尾流模型為基礎(chǔ),通過(guò)粒子群優(yōu)化算法對(duì)風(fēng)電場(chǎng)風(fēng)機(jī)布置進(jìn)行了排布優(yōu)化。相對(duì)于風(fēng)電場(chǎng)機(jī)位優(yōu)化布置研究,國(guó)內(nèi)外在風(fēng)電測(cè)試場(chǎng)風(fēng)力機(jī)優(yōu)化排布方面研究較少。為了彌補(bǔ)風(fēng)電測(cè)試認(rèn)證技術(shù)方面的缺陷,我國(guó)有關(guān)單位分別在內(nèi)蒙古和山東策劃并建立了中小型野外風(fēng)電測(cè)試場(chǎng)。這一舉措為風(fēng)力機(jī)市場(chǎng)提供了強(qiáng)有力的監(jiān)督手段,填補(bǔ)了我國(guó)中小型風(fēng)電測(cè)試認(rèn)證技術(shù)的空白。同時(shí)也推動(dòng)了風(fēng)力發(fā)電產(chǎn)業(yè)的發(fā)展[8]。為了縮短風(fēng)電測(cè)試場(chǎng)測(cè)試周期并提高其經(jīng)濟(jì)效益,本文主要以東營(yíng)沿海野外中小型風(fēng)電測(cè)試場(chǎng)內(nèi)風(fēng)力機(jī)優(yōu)化布置為研究目的,探究測(cè)試場(chǎng)風(fēng)力機(jī)機(jī)位最優(yōu)布置方法。
東營(yíng)市中小型風(fēng)力機(jī)測(cè)試認(rèn)證中心位于黃河入??谌侵薜貐^(qū),地勢(shì)平坦,具備優(yōu)良的風(fēng)能資源[9]。依據(jù)IEC標(biāo)準(zhǔn)對(duì)六臺(tái)被選風(fēng)力發(fā)電機(jī)進(jìn)行合理布置,該六臺(tái)風(fēng)力機(jī)機(jī)型分別為5 kW傘形風(fēng)力機(jī)、10 kW三葉片風(fēng)力機(jī)、30 kW三葉片風(fēng)力機(jī)、100 kW 三葉片風(fēng)力機(jī)、3 kW五葉片風(fēng)力機(jī)和1.5 kW五葉片風(fēng)力機(jī)。風(fēng)力機(jī)具體參數(shù)如表1所示。
對(duì)中小型風(fēng)力機(jī)認(rèn)證測(cè)試場(chǎng)進(jìn)行全年風(fēng)資源數(shù)據(jù)觀測(cè),并經(jīng)過(guò)數(shù)據(jù)篩選得出當(dāng)?shù)孛吭嘛L(fēng)向玫瑰圖。按照其分布規(guī)律可分為春夏秋冬四季的風(fēng)向玫瑰圖。圖1為中小型風(fēng)電認(rèn)證測(cè)試場(chǎng)一年之中四個(gè)季節(jié)的風(fēng)向玫瑰圖。由圖1(a)和(b)可知,SSE(南南東)方向?yàn)榇合膬杉镜氖⑿酗L(fēng)向,NE(北東)方向?yàn)榇合膬杉镜诙黠L(fēng)向。由圖1(c)和(d)可知,秋冬季節(jié)盛行風(fēng)向在157.5°~180°之間,第二主風(fēng)向在315°~337.5°之間。由圖2全年風(fēng)向玫瑰圖可知, SSE(南南東)方向?yàn)槿晔⑿酗L(fēng)向, NE(北東)方向?yàn)槿甑诙⑿酗L(fēng)向。
表1 被測(cè)風(fēng)力機(jī)參數(shù)表
圖1 中小型風(fēng)電認(rèn)證測(cè)試場(chǎng)一年四個(gè)季節(jié)的風(fēng)向玫瑰圖
圖2 中小型風(fēng)電認(rèn)證測(cè)試場(chǎng)全年的風(fēng)向玫瑰圖
考慮到第二主風(fēng)向和盛行風(fēng)向的影響,結(jié)合IEC 61400-12-1:2005《風(fēng)力發(fā)電機(jī)組 功率特性測(cè)試》標(biāo)準(zhǔn)對(duì)給定輪轂高度、風(fēng)輪直徑的被選機(jī)型作出具體布置要求,提出“梯形布置”和“分散布置”兩種測(cè)試場(chǎng)布置方案,如圖3所示。
圖3 中小型風(fēng)電測(cè)試場(chǎng)布置方案
根據(jù)春夏兩季和秋冬兩季風(fēng)向分布規(guī)律,在WAsP軟件中分別對(duì)兩種布置方案進(jìn)行年發(fā)電量和尾流損失的計(jì)算并匯入表2。對(duì)比“分散布置”和“梯形布置”兩種方案,無(wú)論春夏兩季還是秋冬兩季,“分散布置”方案下的總尾流損失因子較小,年發(fā)電量較大。繼而結(jié)合全年測(cè)風(fēng)數(shù)據(jù)對(duì) “分散布置”方案進(jìn)行尾流損失和年發(fā)電量計(jì)算,得出風(fēng)電測(cè)試場(chǎng)在使用“分散布置”方案時(shí),尾流損失最低為0.27%,年發(fā)電量最高為142 MWh,因此將“分散布置”作為風(fēng)電測(cè)試場(chǎng)的初步布置方案,并利用數(shù)值模擬的方法深入探究尾流損失的具體原因。
依據(jù)實(shí)際情況,可利用非穩(wěn)態(tài)的連續(xù)性方程來(lái)描述來(lái)流風(fēng)[10],可由式(1)所表述:
(1)
表2 不同方案下測(cè)試場(chǎng)尾流損失及年發(fā)電量計(jì)算
式中:ρ—?dú)饬髅芏?,kg·m-3;t—時(shí)間,s;u—速度矢量在x方向上的分量;v—速度矢量在y方向上的分量;w—速度矢量在z方向上的分量;μ—?jiǎng)恿φ承韵禂?shù),Pa·s;p—壓力,Pa。
在實(shí)際情況中,流體密度往往是發(fā)生變化的,此時(shí)可用如式(2)關(guān)系進(jìn)行描述:
(2)
式中:ρ—?dú)饬髅芏龋琸g·m-3;t—時(shí)間,s;u—速度矢量在x方向上的分量;v—速度矢量在y方向上的分量;w—速度矢量在z方向上的分量;μ—?jiǎng)恿φ承韵禂?shù),Pa·s;p—壓力,Pa。
利用CFD數(shù)值模擬的方法對(duì)尾流損失的具體原因進(jìn)行驗(yàn)證分析。首先對(duì)測(cè)試場(chǎng)進(jìn)行建模并導(dǎo)入六臺(tái)風(fēng)力機(jī)模型,然后進(jìn)行網(wǎng)格劃分和設(shè)置邊界條件,經(jīng)相關(guān)網(wǎng)格無(wú)關(guān)性驗(yàn)證得出,當(dāng)網(wǎng)格數(shù)從700萬(wàn)變至900萬(wàn)時(shí),網(wǎng)格數(shù)量的變化對(duì)計(jì)算結(jié)果的影響極小,因此該網(wǎng)格數(shù)量可作為計(jì)算網(wǎng)格。最后利用CFX求解器進(jìn)行求解。其中,東南南(ESS)方向作為測(cè)試場(chǎng)盛行風(fēng)向,東北(NE)方向作為測(cè)試場(chǎng)第二主風(fēng)向,根據(jù)盛行風(fēng)向和第二主風(fēng)向的特點(diǎn),進(jìn)行了兩種不同方案的數(shù)值模擬來(lái)探究風(fēng)機(jī)的尾流速度分布規(guī)律。具體布置方案為方案一的入口為盛行風(fēng)向,出口為反方向。方案二的入口為第二主風(fēng)向,出口為反方向。如圖4所示,(a)和(b)分別為方案一和方案二模型立體圖,(c)和(d)為模型的俯視圖,圖中多邊形角點(diǎn)處為風(fēng)力機(jī)安裝位置。
圖4 測(cè)試場(chǎng)模型圖
不同方案下入口風(fēng)向都與測(cè)試場(chǎng)模型的入口邊界相垂直,入口距離前排風(fēng)力機(jī)長(zhǎng)度均為2 Df(Df為前排風(fēng)力機(jī)的風(fēng)輪直徑)。出口距離后排風(fēng)力機(jī)長(zhǎng)度為14 Db和9 Db(Db為出口處風(fēng)力機(jī)的風(fēng)輪直徑,由于方案二中后排風(fēng)力機(jī)風(fēng)輪直徑大,因此其尾流距離較長(zhǎng),為了保證整體計(jì)算的準(zhǔn)確性,因此將該方案出口距離后排風(fēng)力機(jī)長(zhǎng)度設(shè)置為9 Db)。兩側(cè)距離鄰近兩側(cè)風(fēng)力機(jī)的長(zhǎng)度分別為2 Ds1和2 Ds2(Ds1和Ds2分別為左右兩側(cè)風(fēng)力機(jī)的風(fēng)輪直徑)。測(cè)試場(chǎng)模型高為H0+Hh(H0和Hh分別為容量最大風(fēng)力機(jī)的輪轂高度和風(fēng)輪直徑)。
經(jīng)計(jì)算可得不同方案下風(fēng)力機(jī)尾流速度分布曲線(xiàn)如圖5所示。由圖5(a)可看出六臺(tái)風(fēng)機(jī)尾流速度曲線(xiàn)分布趨勢(shì)較為相似,尾流速度均是從5D處開(kāi)始恢復(fù),到18D處基本恢復(fù)至來(lái)流風(fēng)速。由圖5(b)方案二可得A機(jī)位和E機(jī)位風(fēng)力機(jī)的尾流正常發(fā)揮,速度曲線(xiàn)較為平滑,但其余B機(jī)位、C機(jī)位、D機(jī)位和F機(jī)位風(fēng)力機(jī)的尾流速度分布變化較為明顯。
圖5 不同方案下風(fēng)機(jī)尾流速度分布曲線(xiàn)
為了分析得出上述方案中影響機(jī)位之間尾流速度的原因,作出相互影響的BF機(jī)位和CD機(jī)位的布置方法簡(jiǎn)化圖,如圖6所示。
圖6 BF機(jī)位位置簡(jiǎn)圖
圖7為BF機(jī)位縱向的尾流速度分布圖,對(duì)比分析兩種方案下的速度分布曲線(xiàn)。由圖7(a)可知,B機(jī)位風(fēng)力機(jī)尾流速度曲線(xiàn)在下游7D處突然開(kāi)始變得平緩,直到10D處速度才開(kāi)始恢復(fù)上升趨勢(shì),到18D處尾流速度基本上升到來(lái)流風(fēng)速。結(jié)合圖6 BF機(jī)位位置簡(jiǎn)圖分析可知,這是因?yàn)橄嘛L(fēng)向7D到10D處F機(jī)位風(fēng)輪的旋轉(zhuǎn)對(duì)此區(qū)域內(nèi)尾流速度的恢復(fù)造成了一定影響,致使 B機(jī)位風(fēng)機(jī)尾流速度恢復(fù)在此區(qū)域內(nèi)突然變得平緩,雖然在10D 處B機(jī)位尾流速度開(kāi)始恢復(fù),但是F機(jī)位在此處的尾流發(fā)揮使得B機(jī)位尾流恢復(fù)減慢,到18D 處基本恢復(fù)至來(lái)流風(fēng)速。由圖7(b)可以看出,相較于方案二,方案一尾流速度曲線(xiàn)更為平滑,結(jié)合圖6可知,這是因?yàn)镕機(jī)位始終受到B機(jī)位的尾流干擾,不但造成了F機(jī)位對(duì)風(fēng)能的吸收比較消極,而且造成了F機(jī)位的尾流速度出現(xiàn)波動(dòng),影響了其尾流速度的恢復(fù)。
圖7 BF機(jī)位縱向尾流速度分布曲線(xiàn)圖
圖8為BF機(jī)位橫向的尾流速度分布圖。由圖8(a)可知,B機(jī)位下風(fēng)向5D處速度曲線(xiàn)較為對(duì)稱(chēng),且速度變化區(qū)域在-0.75D到0.75D之間。下風(fēng)向8D處的速度曲線(xiàn)從-1.5D處開(kāi)始出現(xiàn)波動(dòng),到輪轂中心處呈現(xiàn)上升趨勢(shì),到0.5D處速度曲線(xiàn)恢復(fù)平穩(wěn)。結(jié)合圖6分析可知,因?yàn)榇藚^(qū)域恰好處于F機(jī)位風(fēng)輪旋轉(zhuǎn)處,由于風(fēng)輪的旋轉(zhuǎn)對(duì)此區(qū)域造成擾動(dòng)使得B機(jī)位風(fēng)機(jī)尾流速度出現(xiàn)波動(dòng)。B機(jī)位下風(fēng)向12D處受F機(jī)位尾流影響速度曲線(xiàn)從-1.5D開(kāi)始呈現(xiàn)下降趨勢(shì),到-0.75D時(shí),最低速度為10 m/s。但隨后從-0.75D開(kāi)始,B機(jī)位進(jìn)入到單獨(dú)尾流區(qū)域,尾流速度開(kāi)始恢復(fù)平穩(wěn)。由圖8(b)可知,F(xiàn)機(jī)位下風(fēng)向5D、8D、12D處橫向尾流速度變化較為相似,均是從-1D處開(kāi)始下降,到輪轂中心處尾流速度降到最低,隨后呈現(xiàn)上升趨勢(shì)。但輪轂中心右側(cè)速度始終低于左側(cè),由圖6可知,這是因?yàn)镕機(jī)位下風(fēng)向0.75D到1.5D區(qū)域處于B機(jī)位尾流區(qū)域內(nèi),B機(jī)位尾流發(fā)揮對(duì)F機(jī)位尾流速度恢復(fù)造成了干擾,使其尾流速度始終低于-1.5D到-0.75D的速度。
圖8 BF機(jī)位下風(fēng)向尾流速度橫向分布曲線(xiàn)圖
沿第二主風(fēng)向CD機(jī)位位置簡(jiǎn)化圖如圖9所示。
圖9 CD機(jī)位風(fēng)力機(jī)位置簡(jiǎn)圖
圖10為CD機(jī)位尾流速度縱向分布曲線(xiàn)圖。由圖10(a)可知,相較于方案一,方案二下的尾流速度曲線(xiàn)更為平滑。結(jié)合圖9可知,這是因?yàn)镈機(jī)位處于C機(jī)位的尾流區(qū)域內(nèi),D機(jī)位風(fēng)輪的旋轉(zhuǎn)對(duì)周?chē)諝猱a(chǎn)生了擾動(dòng),因此減慢了C機(jī)位尾流的恢復(fù)。由圖10(b)可知,在5D到7D處,方案二下的尾流速度略低于方案一,這是因?yàn)镈機(jī)位位于C機(jī)位的尾流區(qū)域內(nèi),C機(jī)位的尾流發(fā)揮不但使D機(jī)位對(duì)于風(fēng)能的吸收極為被動(dòng),還造成了其尾流速度恢復(fù)減慢。但由于兩機(jī)位距離較遠(yuǎn), D機(jī)位在下風(fēng)向7D 后尾流速度開(kāi)始恢復(fù)正常。
圖10 CD機(jī)位尾流速度縱向分布曲線(xiàn)圖
CD機(jī)位尾流速度橫向曲線(xiàn)圖,如圖11所示。由圖11(a)可得,C機(jī)位下風(fēng)向5D、8D和12D處尾流速度從-1.5D開(kāi)始均產(chǎn)生了波動(dòng),且整個(gè)橫向尾流速度曲線(xiàn)為非對(duì)稱(chēng)分布,其中8D位置處的波動(dòng)最為明顯。這是因?yàn)镈機(jī)位恰好位于C機(jī)位下風(fēng)向5D和8D之間偏左側(cè)的位置,由于D機(jī)位風(fēng)輪吸收部分能量使得空氣流速降低,再加上風(fēng)輪的旋轉(zhuǎn)對(duì)周?chē)諝猱a(chǎn)生擾動(dòng),造成了此區(qū)域內(nèi)湍流強(qiáng)度較大,因此C機(jī)位在此區(qū)域內(nèi)的尾流速度波動(dòng)較大。但在距離較遠(yuǎn)的12D處,C機(jī)位受其尾流影響逐漸降低,尾流速度開(kāi)始恢復(fù)。且輪轂中心右側(cè)下風(fēng)向不同位置處的尾流速度均恢復(fù)至來(lái)流風(fēng)速。
圖11 CD機(jī)位下風(fēng)向尾流速度橫向分布曲線(xiàn)圖
由圖11(b)可知,D機(jī)位下風(fēng)向不同位置處的橫向尾流速度分布趨勢(shì)都大體相似,且為對(duì)稱(chēng)分布。尾流速度均是從-1.5D處開(kāi)始下降,到輪轂中心處降為最低,隨后開(kāi)始恢復(fù)至來(lái)流風(fēng)速。但其中5D處輪轂中心右側(cè)速度略低于左側(cè),這是因?yàn)镈機(jī)位處于C機(jī)位下風(fēng)向偏左側(cè)的位置,D機(jī)位右側(cè)尾流速度恢復(fù)受C機(jī)位尾流發(fā)揮的影響,但由于D機(jī)位風(fēng)輪直徑較大,其下風(fēng)向5D處距離已經(jīng)足夠遠(yuǎn),因此C機(jī)位的尾流速度在此位置已經(jīng)基本恢復(fù)至來(lái)流風(fēng)速,對(duì)D機(jī)位尾流速度影響極小。
在保證盛行風(fēng)向下各機(jī)型尾流正常發(fā)揮,避免沿第二主風(fēng)向下BF機(jī)位和CD機(jī)位尾流相互影響的情況下,再加入WN(西北)風(fēng)向來(lái)重新校準(zhǔn)機(jī)位布置。得到如圖12所示的測(cè)試場(chǎng)機(jī)位布置調(diào)整圖,其中,虛線(xiàn)為初步方案機(jī)位連線(xiàn),實(shí)線(xiàn)為調(diào)整后的機(jī)位連線(xiàn)。
圖12 測(cè)試場(chǎng)機(jī)位布置調(diào)整圖
為驗(yàn)證重新調(diào)整后的機(jī)位布置的可行性,利用WAsP軟件計(jì)算了調(diào)整后方案下的年發(fā)電量和尾流損失,結(jié)果如表3所示??梢钥吹剑臧l(fā)電量提高到了142.323 MWh,尾流損失降低到0.219%。
表3 測(cè)試場(chǎng)最優(yōu)布置方案下風(fēng)力機(jī)尾流損失及發(fā)電量計(jì)算統(tǒng)計(jì)
由此可以確定,調(diào)整后的方案為最優(yōu)布置方案,其中的尾流損失可能是其他風(fēng)向的影響造成的。最終得出的最優(yōu)布置方案,不僅增加了測(cè)試場(chǎng)的年發(fā)電量,而且還減少了風(fēng)機(jī)的尾流損失。且該布置方案綜合研究了盛行風(fēng)向、第二主風(fēng)向和第三風(fēng)向下機(jī)位之間的尾流影響,提高了測(cè)試場(chǎng)的經(jīng)濟(jì)效益和風(fēng)能的利用率,有效地縮短了風(fēng)力機(jī)的測(cè)試周期。
本文對(duì)不同類(lèi)型的風(fēng)力機(jī)在沿海野外中小型風(fēng)電認(rèn)證測(cè)試場(chǎng)中的機(jī)位布置進(jìn)行了研究,依據(jù)該地區(qū)四季的風(fēng)資源分布情況,對(duì)測(cè)試場(chǎng)的風(fēng)力機(jī)布置提出“梯形”和“分散”兩種布置方案,利用WAsP軟件進(jìn)行年發(fā)電量和尾流損失的計(jì)算,得到其年發(fā)電量和尾流損失,將“分散布置”方案作為初步布置方案,為探究尾流損失的具體原因。利用CFD對(duì)其進(jìn)行數(shù)值模擬計(jì)算,得到改方案下機(jī)位之間尾流相互影響的具體原因,并對(duì)其作出調(diào)整,并最終確定了最優(yōu)布置方案,最優(yōu)布置方案下年發(fā)電量從142.246 MWh增加到了142.323 MWh,尾流損失從0.273%減少到了0.219%。