陳世游,陸 海,蘇 適
(云南電網(wǎng)有限責(zé)任公司電力科學(xué)研究院工作站,云南 昆明 650217)
風(fēng)能作為清潔便利的能源,是風(fēng)電開(kāi)發(fā)利用的重要載體。過(guò)去10 年,我國(guó)風(fēng)電裝機(jī)增長(zhǎng)了47倍,2010年風(fēng)力發(fā)電累計(jì)裝機(jī)容量達(dá)到4182.7 萬(wàn)千瓦,躍居世界第一。預(yù)計(jì)到2020年,我國(guó)風(fēng)電裝機(jī)容量將至少達(dá)到1.5億千瓦,約占10%。由此可見(jiàn),風(fēng)能最具規(guī)?;_(kāi)發(fā)前景。目前國(guó)內(nèi)風(fēng)電機(jī)組能量可利用率水平較低一般為60%~70%[1],而國(guó)外已經(jīng)達(dá)到了90%以上,即使在不考慮限電的影響情況下,機(jī)組的能量可利用率仍然有較大的提升空間。隨著風(fēng)電大規(guī)模、集中式接入主網(wǎng),如何合理地評(píng)估風(fēng)電場(chǎng)運(yùn)行情況及為調(diào)度提供可靠的風(fēng)電場(chǎng)運(yùn)行信息,如理論出力,對(duì)指導(dǎo)風(fēng)電場(chǎng)優(yōu)化控制運(yùn)行具有重要意義。
隨著風(fēng)力發(fā)電的大規(guī)模應(yīng)用,對(duì)風(fēng)電場(chǎng)全場(chǎng)風(fēng)機(jī)出力研究也逐漸展開(kāi)。由于風(fēng)力發(fā)電機(jī)組的理論發(fā)電能力只與機(jī)組參數(shù)和入流的風(fēng)速有關(guān),風(fēng)機(jī)艙前風(fēng)速精準(zhǔn)測(cè)量是評(píng)估單臺(tái)及全場(chǎng)機(jī)組理論出力重要前提。然而,在針對(duì)高原地形風(fēng)電場(chǎng)多分布在山脊上和立體氣候風(fēng)資源隨海拔變化顯著的情況下,紊流、湍流和風(fēng)向迂回的相互交錯(cuò)作用,這種傳統(tǒng)的測(cè)量方法導(dǎo)致測(cè)風(fēng)塔測(cè)得的風(fēng)數(shù)據(jù)與風(fēng)機(jī)艙前風(fēng)速不一致?,F(xiàn)有文獻(xiàn)主要從風(fēng)電場(chǎng)出力建模和以精確測(cè)風(fēng)為基礎(chǔ)進(jìn)行功率曲線(xiàn)繪制來(lái)評(píng)估單機(jī)出力進(jìn)行研究,文獻(xiàn)[2]利用雷達(dá)測(cè)風(fēng)儀激光反射原理進(jìn)行風(fēng)速風(fēng)向測(cè)量精準(zhǔn)測(cè)量,利用區(qū)間分析的方法繪制出功率曲線(xiàn)風(fēng)力發(fā)電機(jī)組的功率特性曲線(xiàn),但是文獻(xiàn)給出方法只適合去評(píng)估單臺(tái)風(fēng)機(jī)出力,無(wú)法對(duì)全場(chǎng)風(fēng)機(jī)出力進(jìn)行評(píng)估。文獻(xiàn)[3]以風(fēng)電場(chǎng)的實(shí)際風(fēng)速及其發(fā)電出力的歷史實(shí)測(cè)數(shù)據(jù)為基礎(chǔ),通過(guò)對(duì)風(fēng)電場(chǎng)內(nèi)部多個(gè)機(jī)群的風(fēng)速序列進(jìn)行相關(guān)性處理,建立由多種Copula函數(shù)組成的混合Copula模型,然后將多組機(jī)群出力線(xiàn)性疊加,得到風(fēng)電場(chǎng)聚合出力。文獻(xiàn)在進(jìn)行全場(chǎng)出力評(píng)估時(shí),需要大量機(jī)群的精確風(fēng)速歷史數(shù)據(jù)。而在高原立體氣候條件下,獲得各個(gè)機(jī)群精確風(fēng)速困難性大。文獻(xiàn)[4]采用 CFD數(shù)值模擬評(píng)估復(fù)雜山地地形中的風(fēng)能分布,采用了六面體結(jié)構(gòu)化網(wǎng)格系統(tǒng),計(jì)算收斂迅速, 結(jié)果精確。文獻(xiàn)為評(píng)估全場(chǎng)風(fēng)資源提供了可靠方案,為全場(chǎng)風(fēng)機(jī)出力評(píng)估提供科學(xué)方法。文獻(xiàn)[5]提出了一種真實(shí)功率曲線(xiàn)的模型,考慮了各范圍風(fēng)速數(shù)據(jù)的正態(tài)分布。此外,還提出了一種基于蒙特卡羅的仿真技術(shù),以基于正規(guī)的模型來(lái)產(chǎn)生風(fēng)速數(shù)據(jù),為風(fēng)電場(chǎng)出力建模提供可靠數(shù)據(jù)。
現(xiàn)有研究針對(duì)高原立體氣候條件下,沒(méi)有對(duì)風(fēng)電場(chǎng)出力研究提供可行的實(shí)驗(yàn)方案,因此有待深入研究;在此,本文提出采用一臺(tái)激光雷達(dá)裝置對(duì)機(jī)組進(jìn)行風(fēng)速精準(zhǔn)測(cè)量;然后,以激光雷達(dá)的風(fēng)速實(shí)測(cè)數(shù)據(jù)為參考點(diǎn),對(duì)風(fēng)電場(chǎng)風(fēng)資源分布進(jìn)行CFD模擬推算風(fēng)電場(chǎng)內(nèi)機(jī)組的艙前風(fēng)速,最后,采用推算的風(fēng)速和風(fēng)機(jī)實(shí)際輸出功率進(jìn)行擬合,將擬合功率曲線(xiàn)與廠家提供的靜態(tài)功率曲線(xiàn)做對(duì)比,驗(yàn)證風(fēng)場(chǎng)理論出力評(píng)估方法的可行性。本文研究流程圖如圖 1所示。
激光雷達(dá)是采用激光多普勒原理,對(duì)空氣中的微粒的雷達(dá)回波進(jìn)行追蹤,從而對(duì)風(fēng)速進(jìn)行測(cè)量的裝置。與傳統(tǒng)機(jī)械式測(cè)風(fēng)相比,激光雷達(dá)可以采集全面風(fēng)數(shù)據(jù):風(fēng)速,風(fēng)向,風(fēng)切變,向變,湍流強(qiáng)度,氣流傾角等數(shù)據(jù);另外其測(cè)量原理決定了可以對(duì)機(jī)組前方自由流的風(fēng)速進(jìn)行測(cè)量,保證了測(cè)量數(shù)據(jù)的可靠性。同時(shí),在風(fēng)機(jī)尾流、大氣現(xiàn)象、復(fù)雜氣流的可視化方面也具有較強(qiáng)的優(yōu)勢(shì)[6]。
激光雷達(dá)風(fēng)速測(cè)量的示意如圖2所示。光學(xué)頭發(fā)出兩束夾角為2θ的激光束,并采用多普勒原理分別測(cè)量沿激光束的徑向風(fēng)速,一旦 R WS1和 R WS0被測(cè)出,那水平風(fēng)速可以采用如下的方法進(jìn)行重構(gòu)[7]。
徑向風(fēng)速的計(jì)算公式為:
RWS=徑向風(fēng)速=沿激光束的風(fēng)速=
其中λ為激光波長(zhǎng),f為多普勒頻移量
水平風(fēng)速的兩個(gè)分量可以采用下式得到
則風(fēng)速矢量和風(fēng)向可以通過(guò)下式得到
圖2 激光雷達(dá)測(cè)風(fēng)示意圖Fig.2 Laser radar wind measurement
CFD(Computational Fluid Dynamics)計(jì)算流體力學(xué)技術(shù)是通過(guò)計(jì)算機(jī)的幫助,根據(jù)流體力學(xué)的規(guī)律進(jìn)行模擬求解,將幾何區(qū)域分解為一個(gè)個(gè)小的網(wǎng)格,并在其中求解復(fù)雜的偏微分方程組的過(guò)程[8]。綜合計(jì)算精度,運(yùn)算時(shí)間,計(jì)算機(jī)性能,采用雷諾時(shí)均法對(duì)湍流進(jìn)行模擬,考慮地表粗糙度,本次計(jì)算采用WindSIM軟件帶有壁面函數(shù)的k-s瑞流模型,三維標(biāo)準(zhǔn)k-s模型方程為:
其中k和ε分別為湍動(dòng)能及其耗散率,kP為湍動(dòng)能生成項(xiàng),可以用如下公式表示:
在式6中,ju和ju′分別為 j方向的平均風(fēng)速和脈動(dòng)風(fēng)速。其中湍流粘度為:
式4、5和7中:Cμ,1Cε,2Cε,kσ,εσ為模型常數(shù),需要根據(jù)測(cè)風(fēng)數(shù)據(jù)對(duì)模型常數(shù)進(jìn)行修訂。
利用地球電子地形數(shù)據(jù)ASTER GDEM(先進(jìn)星載熱發(fā)射和反射輻射儀全球數(shù)字高程模型)獲得地形數(shù)據(jù);根據(jù)提供的1:2000的測(cè)繪地形圖(DWG格式),采用AutoCAD和Global mapper、Wasp Map Editor等軟件對(duì)地形數(shù)據(jù)進(jìn)行處理,地形圖處理步驟如下[9]。
(1)地形圖的外擴(kuò)
本次計(jì)算除了風(fēng)電場(chǎng)區(qū)域外,還需要包含周邊5000m的區(qū)域,外擴(kuò)區(qū)域地形圖采用 ASTER電子地圖數(shù)據(jù),該數(shù)據(jù)較能真實(shí)的反映真實(shí)的地形特征。
通過(guò)Autocad、Global mapper等軟件對(duì)原始地形數(shù)據(jù)和下載地形數(shù)據(jù)進(jìn)行合并處理,優(yōu)化邊界等高線(xiàn),保證計(jì)算收斂性。
(2)地形圖處理及生成
通過(guò) cad對(duì)地形數(shù)據(jù)中的無(wú)效信息進(jìn)行剔除,僅保留等高線(xiàn)信息;
通過(guò) Global mapper軟件對(duì)等高線(xiàn)信息進(jìn)行處理,輸出map、gws等需要的格式
依靠國(guó)家基礎(chǔ)地理信息中心研制的全球高分辨率地表覆蓋遙感制圖獲得 30米地表覆蓋數(shù)據(jù)(GlobeLand30),通過(guò)獲取GLC30粗糙度信息,截取風(fēng)電場(chǎng)區(qū)域范圍內(nèi)的粗糙度信息,導(dǎo)入 global mapper中;或者根據(jù)提供的測(cè)繪地形圖中的關(guān)于對(duì)地貌信息的描述,手動(dòng)生成地貌信息文件,導(dǎo)入global mapper文件;對(duì)導(dǎo)入的文件進(jìn)行處理,輸出需要的格式。
采用Global mapper或Wasp Map Editor對(duì)地形數(shù)據(jù)進(jìn)行處理及合并,根據(jù)所采用的軟件所需格式輸出。Windsim需要gws文件,WT需要map文件。
針對(duì)風(fēng)電場(chǎng)范圍及周?chē)?5000米的范圍進(jìn)行CFD模擬仿真。各項(xiàng)參數(shù)設(shè)置如下。
(1)水平分辨率
水平分辨率為25 m~267 m。對(duì)風(fēng)電場(chǎng)區(qū)域及測(cè)風(fēng)塔區(qū)域進(jìn)行網(wǎng)格加密,網(wǎng)格精度為25 m*25 m,風(fēng)電場(chǎng)以外的區(qū)域采用 0.3的系數(shù)逐漸擴(kuò)大,在不影響計(jì)算結(jié)果的情況下,節(jié)省計(jì)算機(jī)資源,提高計(jì)算速度。
(2)垂直網(wǎng)格數(shù)量
垂直網(wǎng)格數(shù)量采用31個(gè),并通過(guò)參數(shù)設(shè)置保證150 m以下存在10個(gè)以上網(wǎng)格,以確保模擬結(jié)果的準(zhǔn)確性。
(3)平滑設(shè)置
實(shí)際地形計(jì)算結(jié)果收斂性較好,未對(duì)地形圖進(jìn)行平滑處理
(4)森林模型
因無(wú)區(qū)域森林?jǐn)?shù)據(jù),未采用森林模型,采用粗糙度設(shè)置替代。通過(guò)以上設(shè)置,本次計(jì)算x方向網(wǎng)格數(shù)量410個(gè),y方向網(wǎng)格數(shù)量418個(gè),z方向網(wǎng)格數(shù)量31個(gè),總網(wǎng)格數(shù)量500多萬(wàn)個(gè)。
表1 風(fēng)電場(chǎng)邊界條件設(shè)置參數(shù)Tab.1 Wind farm boundary conditions setting parameters
計(jì)算風(fēng)電場(chǎng)位于中國(guó)云貴高原地區(qū),所在區(qū)域海拔高度1800~2400 m,風(fēng)電場(chǎng)內(nèi)高度落差約500 m。風(fēng)電場(chǎng)占地面積約為30 km2,境內(nèi)溝壑縱橫,群山高聳、山谷相間、 立體氣候明顯,區(qū)域內(nèi)存在大量藥用值物。風(fēng)電場(chǎng)內(nèi)共有66臺(tái)東方電氣FD82機(jī)型風(fēng)電機(jī)組,風(fēng)機(jī)輪轂高度為70 m,葉輪直徑為82 m,額定容量 1500 KW。場(chǎng)內(nèi)風(fēng)電機(jī)組的排布如圖 3所示。
本次試驗(yàn)所采用的激光雷達(dá)為 AVENT公司生產(chǎn)的 WindIris機(jī)艙式激光雷達(dá)。在實(shí)際測(cè)試時(shí),處理單元放置在風(fēng)力機(jī)機(jī)艙內(nèi)部,光學(xué)頭和通過(guò)三腳架打孔固定安裝于機(jī)艙頂部如圖4所示。
圖3 機(jī)組排布圖Fig.3 Unit layout diagram
圖4 機(jī)艙頂安裝方式Fig.4 Nacelle top installation
本次測(cè)試所采用的激光雷達(dá)主要技術(shù)參數(shù)如表2所示。
表2 W indIris技術(shù)參數(shù)表Tab.2 W indIris technical parameters
3.2.1 雷達(dá)數(shù)據(jù)處理
通過(guò)機(jī)載激光裝置,測(cè)得單臺(tái)風(fēng)機(jī)艙前不同距離處的風(fēng)速如表3所示。從表3所測(cè)的數(shù)據(jù)來(lái)看風(fēng)速在距離風(fēng)機(jī)110 m~200 m之間基本穩(wěn)定。
為進(jìn)一步清晰描述不同測(cè)風(fēng)距離測(cè)風(fēng)數(shù)據(jù)相關(guān)性,引入相關(guān)系數(shù)來(lái)刻畫(huà)變量之間關(guān)系密切程度,相關(guān)系數(shù)的取值區(qū)間在1到–1之間。1表示兩個(gè)變量完全線(xiàn)性相關(guān),–1表示兩個(gè)變量完全負(fù)相關(guān),0表示兩個(gè)變量不相關(guān)。數(shù)據(jù)越趨近于0表示相關(guān)關(guān)系越弱。測(cè)風(fēng)數(shù)據(jù)之間相關(guān)性越高,數(shù)據(jù)質(zhì)量越可靠。下式是相關(guān)系數(shù)[10]計(jì)算公式。
其中xyr表示樣本相關(guān)系數(shù),xys 表示樣本協(xié)方差,xs表示X的樣本標(biāo)準(zhǔn)差,ys表示Y的樣本標(biāo)準(zhǔn)差。下面分別是xs和ys標(biāo)準(zhǔn)差的計(jì)算公式。
按照如上公式9、10、11和12我們可以求得數(shù)據(jù)相關(guān)關(guān)系,具體數(shù)值如表4所示,從表4可以看出,測(cè)風(fēng)數(shù)據(jù)之間相關(guān)性非常好,數(shù)據(jù)質(zhì)量可靠。距離風(fēng)輪越遠(yuǎn),風(fēng)速越高。這是由于距離風(fēng)機(jī)越遠(yuǎn),海拔越低,測(cè)量點(diǎn)距離地表越遠(yuǎn)。白天熱交換劇烈,大氣極不穩(wěn)定,對(duì)流強(qiáng)烈,貼近地表位置的風(fēng)速將高于遠(yuǎn)離地表位置。夜間則相反。
3.2.2 CFD模型驗(yàn)證及通過(guò)雷達(dá)數(shù)據(jù)推算全場(chǎng)
將激光雷達(dá)測(cè)量的110 m、140 m、170 m距離處風(fēng)速均導(dǎo)入WindSim軟件中計(jì)算,模擬風(fēng)機(jī)所處上述位置處風(fēng)速,將雷達(dá)測(cè)速與推算作對(duì)比來(lái)驗(yàn)證模型準(zhǔn)確性,結(jié)果如表5所示。
表4 不同距離風(fēng)速相關(guān)性分析統(tǒng)計(jì)表Tab.4 Statistical table of correlation analysis of wind speed at different distances
表5 各個(gè)點(diǎn)處風(fēng)速推算驗(yàn)證Tab.5 Verification of wind speed at various points
110 m測(cè)量結(jié)果和推算結(jié)果偏低,230 m測(cè)量結(jié)果和推算結(jié)果偏高,其余結(jié)果均在8.34 m/s附近,偏差均小于0.1,選取200 m處(2.5D)進(jìn)行推算是可行的。因風(fēng)輪阻擋,造成實(shí)測(cè)風(fēng)速小于實(shí)際風(fēng)速,采用距離近處風(fēng)速推算結(jié)果小于距離遠(yuǎn)處推算結(jié)果,模型并未產(chǎn)生較大誤差,結(jié)果準(zhǔn)確可靠。
將云南某地區(qū)風(fēng)電場(chǎng)各個(gè)風(fēng)機(jī)進(jìn)行 CFD模擬推算,并用對(duì)機(jī)組的 SCADA/風(fēng)速儀數(shù)據(jù)進(jìn)行采集,獲得機(jī)組各個(gè)機(jī)位模擬數(shù)據(jù)。將推算的數(shù)據(jù)與SCADA/風(fēng)速儀采集的數(shù)據(jù)進(jìn)行相關(guān)性分析,結(jié)果如表 6所示。對(duì)比 J02~J05及附近機(jī)位的實(shí)測(cè)SCADA數(shù)據(jù),發(fā)現(xiàn) J02、J03、J04、J05風(fēng)速與周?chē)鷻C(jī)位風(fēng)速變化趨勢(shì)差異較大??赡苁菣C(jī)艙風(fēng)速儀損壞造成或系統(tǒng)錯(cuò)誤。從表中可以看出各個(gè)風(fēng)機(jī)推算得到的艙前風(fēng)速與 SCADA/風(fēng)速儀采集的相關(guān)性較好,模擬效果較好[11]。
3.2.3 功率曲線(xiàn)對(duì)比
采用推算的風(fēng)速和風(fēng)機(jī)實(shí)際輸出功率擬合功率曲線(xiàn),擬合結(jié)果與廠家提供的靜態(tài)功率曲線(xiàn)對(duì)比,如圖5、6所示。
從圖5、6可以看出,功率曲線(xiàn)和Cp曲線(xiàn)變化存在四個(gè)階段:
(1)在低風(fēng)速段(3 m/s~7 m/s),擬合的實(shí)際功率曲線(xiàn)略高于廠家的靜態(tài)功率曲線(xiàn),這主要是因?yàn)殪o態(tài)功率曲線(xiàn)模擬未考慮湍流強(qiáng)度的影響,而實(shí)際功率曲線(xiàn)受到湍流強(qiáng)度影響。在低風(fēng)速段,Cp為上升趨勢(shì),當(dāng)風(fēng)速?gòu)? m/s增加到4 m/s時(shí),風(fēng)能利用系數(shù)Cp從0.2左右增加到0.4左右,風(fēng)速的提升對(duì)Cp的影響非常大。因此在低風(fēng)速段,湍流強(qiáng)度越高,功率曲線(xiàn)表現(xiàn)越好。
表6 模擬結(jié)果與運(yùn)行結(jié)果相關(guān)性統(tǒng)計(jì)Tab.6 Correlation between simulation results and operating results
(2)在8 m/s附近,靜態(tài)功率曲線(xiàn)與動(dòng)態(tài)功率曲線(xiàn)基本一致。因?yàn)榇穗A段,Cp隨風(fēng)速的變化非常小,湍流強(qiáng)度對(duì)風(fēng)速的影響無(wú)法在功率曲線(xiàn)中體現(xiàn)。
(3)在額定功率附近(9 m/s~12 m/s),擬合的實(shí)際功率曲線(xiàn)低于廠家的靜態(tài)功率曲線(xiàn),這是因?yàn)镃p曲線(xiàn)開(kāi)始下滑,這與低風(fēng)速段情況正好相反,因此有上述表現(xiàn)。
圖5 功率曲線(xiàn)對(duì)比Fig.5 Pow er curve comparison
圖6 Cp 曲線(xiàn)對(duì)比Fig.6 Cp curve comparison
(4)在額定風(fēng)速以后(13 m/s~25 m/s),風(fēng)機(jī)實(shí)際輸出功率高于額定功率,這是因?yàn)轱L(fēng)機(jī)在設(shè)計(jì)時(shí)可以超發(fā)10%左右,因此高出靜態(tài)功率曲線(xiàn)。
綜上所述,采用激光雷達(dá)測(cè)風(fēng)與CFD外推全場(chǎng)風(fēng)機(jī)出力方法擬合的功率曲線(xiàn)與實(shí)際較為相符,此方法具有可行性。
本文首先從復(fù)雜地形氣候條件下,傳統(tǒng)測(cè)風(fēng)手段不能精確測(cè)得機(jī)艙前風(fēng)速入手,采用激光雷達(dá)精確測(cè)得云南曲靖某地區(qū)單個(gè)機(jī)艙前不同距離處風(fēng)速,并對(duì)測(cè)風(fēng)數(shù)據(jù)進(jìn)行相關(guān)性分析,論證了激光雷達(dá)測(cè)風(fēng)數(shù)據(jù)可靠性;然后,將激光雷達(dá)測(cè)量的110 m、140 m、170 m距離處風(fēng)速均導(dǎo)入WindSim軟件中計(jì)算,對(duì)全場(chǎng)風(fēng)速進(jìn)行模擬,并將模擬結(jié)果與測(cè)風(fēng)塔處測(cè)得風(fēng)速做對(duì)比,得出 CFD模擬方案可行性;最后,采用CFD模擬推算的風(fēng)速和風(fēng)機(jī)實(shí)際輸出功率進(jìn)行擬合得到功率曲線(xiàn),將擬合結(jié)果與廠家提供的靜態(tài)功率曲線(xiàn)對(duì)比,結(jié)果可知擬合得到功率曲線(xiàn)與廠家提供的靜態(tài)功率曲線(xiàn)較吻合,從而驗(yàn)證了以激光雷達(dá)測(cè)風(fēng)為基礎(chǔ)進(jìn)行 CFD模擬外推全場(chǎng)理論功率方案的可行性。
過(guò)去在評(píng)價(jià)機(jī)組或風(fēng)電場(chǎng)時(shí),基本都是用年利用小時(shí)數(shù),發(fā)電量等指標(biāo)來(lái)進(jìn)行。這樣做的缺點(diǎn)是沒(méi)有考慮到機(jī)組微觀選址時(shí)帶來(lái)的風(fēng)資源天然差異。我們只知道某臺(tái)機(jī)組發(fā)電量大,某臺(tái)小。但是并不知道每臺(tái)機(jī)組實(shí)際上的理論最大出力是多少,也無(wú)法有效地判斷造成這種差異的深層次原因。這時(shí),通過(guò)本文所提方案評(píng)估單機(jī)及全場(chǎng)理論出力,進(jìn)而引入“能量可利用率”這一新的評(píng)價(jià)標(biāo)準(zhǔn),解決上述評(píng)價(jià)方法缺陷,合理地評(píng)估風(fēng)電場(chǎng)運(yùn)行情況及為調(diào)度提供可靠的風(fēng)電場(chǎng)信息,指導(dǎo)風(fēng)電場(chǎng)優(yōu)化控制運(yùn)行,且對(duì)提高新能源區(qū)域消納能力具有重要意義。
[1] 馬平, 劉昌華. 風(fēng)力發(fā)電機(jī)組功率曲線(xiàn)的驗(yàn)證[J]. 可再生能源, 2008, 26(12): 82-84.
[2] 潘寧. 基于激光雷達(dá)測(cè)風(fēng)儀的風(fēng)電機(jī)組功率曲線(xiàn)測(cè)試方法研究[J]. 節(jié)能技術(shù), 2013, 31(178): 112-115.
[3] Pan Xiong, Wang Lili, et al. Modeling method of wind farm output based on mixed Copula function[J]. Automation of Electric Power Systems, 2014, 38(14): 17-22.
[4] S. Gill, B. Stephen, S. Galloway. Wind turbine condition assessment through power curve copula modeling, IEEE Transactions on Sustainable Energy[J]. 2012, 3(1): 94-101.
[5] Kim Jaejin,Baik Jongjin.A numerical study of the effects of ambient wind direction on flow and dispersion in urban street canyons using the RNG k-s turbulence model[J]. Atmospheric Environment, 2004, 38(19): 3039-3048.
[6] 肖創(chuàng)英, 汪寧渤, 陟晶, 等. 甘肅酒泉風(fēng)電出力特性分析[J]. 電力系統(tǒng)自動(dòng)化, 2010, 34(17): 64-67.
[7] 劉昊, 柳亦兵, 辛衛(wèi)東, 李若昭. 基于運(yùn)行數(shù)據(jù)的風(fēng)力發(fā)電機(jī)組功率特性分析[J]. 電網(wǎng)與清潔能源. 2009(7): 53-56.
[8] P. Mohan Kumar, M. Mohan Ram Surya, et al. Comparitive CFD analysis of darrieus wind turbine with NTU-20-V and NACA0018 airfoils, 2017 IEEE International Conference on Smart Grid and Smart Cities [C]. Singapore, 2017.
[9] 肖儀清, 李朝, 等.復(fù)雜地形風(fēng)能評(píng)估的CFD方法[J]. 華南理工大學(xué)學(xué)報(bào), 2009. 37(9): 30-35.
[10] 卞海紅, 鄭維高, 林章歲, 等. 集群風(fēng)電場(chǎng)出力統(tǒng)計(jì)指標(biāo)建模與應(yīng)用[J]. 電力自動(dòng)化設(shè)備, 2015, 35(12): 21-27.