吳天太

(中國三峽新能源(集團(tuán))股份有限公司甘肅分公司，蘭州 730000)

0 引言

目前市場主流的風(fēng)電場發(fā)電量計算軟件WT、WAsP 等均采用年平均空氣密度作為發(fā)電量計算時的空氣密度，然而這種用年平均值替代每時每刻瞬時值的方法，未考慮空氣密度的年變化情況對風(fēng)電場發(fā)電量計算的影響，因此一些相關(guān)的風(fēng)電場發(fā)電量計算軟件開發(fā)商在考慮是否有必要使用空氣密度時間序列來計算風(fēng)電場發(fā)電量[1]。

基于此，本文主要分析空氣密度年變化情況對風(fēng)電場發(fā)電量計算的影響，旨在評估空氣密度年變化情況引起的風(fēng)電場發(fā)電量計算誤差，考慮使用空氣密度時間序列計算風(fēng)電場發(fā)電量的必要性。首先對中國空氣密度的年變化幅度的大致情況進(jìn)行分析；然后選取10 座有代表性的測風(fēng)塔，對其所在區(qū)域的空氣密度年變化情況進(jìn)行統(tǒng)計分析，探究空氣密度的區(qū)域年變化情況，并分析空氣密度的主要影響因素；最后選取其中空氣密度年變化幅度較大的1 座測風(fēng)塔數(shù)據(jù)來分析空氣密度年變化情況對風(fēng)電場發(fā)電量計算的影響。

由于風(fēng)速的年變化情況和空氣密度的年變化情況沒有固定的函數(shù)，很難從理論公式上推導(dǎo)空氣密度年變化情況對風(fēng)電場發(fā)電量的影響，因此，本文以某型號風(fēng)電機(jī)組的單機(jī)理論發(fā)電量計算為例來分析空氣密度年變化情況對風(fēng)電場發(fā)電量計算的影響。

1 中國空氣密度年變化幅度的大致情況

為探討空氣密度年變化情況對風(fēng)電場發(fā)電量計算的影響，必須分析中國空氣密度年變化的幅度，以此來分析各區(qū)域空氣密度年變化情況對風(fēng)電場發(fā)電量影響的大小，同時探究空氣密度的主要影響因素。

1.1 空氣密度的計算公式[2]

空氣密度的計算公式可表示為：

式中：ρ為空氣密度，kg/m3；P為年平均大氣壓力，Pa；T為年平均開氏溫標(biāo)空氣溫度，K，T=273.15+t(其中t為攝氏溫度)；R為氣象常數(shù)，J/(kg·K)，取287。

從式(1)可以看出，空氣密度和年平均大氣壓力成正比，和年平均開氏溫標(biāo)空氣溫度成反比。

1.2 中國代表性區(qū)域的空氣密度年變化情況分析

在中國范圍內(nèi)各種地形和氣候條件下，選取最具代表性的10 座測風(fēng)塔，統(tǒng)計分析這些測風(fēng)塔所在10 個區(qū)域的空氣密度年變化情況[3]。選取的10 座測風(fēng)塔分別位于云南省楚雄州(下文簡稱為“云南楚雄”)、福建省莆田市(下文簡稱為“福建莆田”)、陜西省渭南市(下文簡稱為“陜西渭南”)、新疆維吾爾自治區(qū)烏魯木齊市(下文簡稱為“新疆烏魯木齊”)、黑龍江省同江市(下文簡稱為“黑龍江同江”)、內(nèi)蒙古自治區(qū)烏蘭察布市(下文簡稱為“內(nèi)蒙古烏蘭察布”)、青海省海西州(下文簡稱為“青海海西”)、江蘇省揚(yáng)州市(下文簡稱為“江蘇揚(yáng)州”)、廣東省清遠(yuǎn)市(下文簡稱為“廣東清遠(yuǎn)”)、湖北省宜城市(下文簡稱為“湖北宜城”)，分析每座測風(fēng)塔某1 個完整年的數(shù)據(jù)。

各代表區(qū)域測風(fēng)塔某1 個完整年的空氣密度、空氣溫度、大氣壓力的月平均值及年變化情況分別如表1～表3 所示。表中，年平均值采用10 min 時間序列計算得到。各代表區(qū)域的空氣密度的年變化幅度如圖1 所示。

表1 各代表區(qū)域測風(fēng)塔某1 個完整年的空氣密度的月平均值及年變化情況Table 1 Monthly average value and annual variation of air density in a complete year of anemometer tower in each representative region(單位：kg/m3)

表2 各代表區(qū)域測風(fēng)塔某1 個完整年的空氣溫度的月平均值及年變化情況Table 2 Monthly average value and annual variation of air temperature in a complete year of anemometer tower in each representative region(單位：℃)

表3 各代表區(qū)域測風(fēng)塔某1 個完整年的大氣壓力的月平均值及年變化情況Table 3 Monthly average value and annual variation of atmospheric pressure in a complete year of anemometer tower in each representative region(單位：kPa)

(續(xù)表)

圖1 各代表區(qū)域的空氣密度的年變化幅度Fig. 1 Annual variation range of air density in each representative region

從表1～表3 和圖1 可以看出：在中國，空氣溫度年變化幅度較大的地區(qū)(例如內(nèi)蒙古烏蘭察布和黑龍江同江)，其空氣密度的年變化幅度也較大，空氣密度的月平均值標(biāo)準(zhǔn)差也較大，在0.050 kg/m3以上；空氣溫度的年變化幅度較小的地區(qū)，空氣密度的年變化幅度也很小，空氣密度的月平均值標(biāo)準(zhǔn)差也很小。

1.3 空氣密度變化幅度的主要影響因素

為分析空氣密度的主要影響因素，對空氣密度標(biāo)準(zhǔn)差、空氣溫度標(biāo)準(zhǔn)差、大氣壓力標(biāo)準(zhǔn)差的相關(guān)關(guān)系進(jìn)行分析。

對空氣密度標(biāo)準(zhǔn)差和空氣溫度標(biāo)準(zhǔn)差做線性相關(guān)，相關(guān)關(guān)系圖如圖2 所示。相關(guān)方程式為y=0.0046x-0.0001，決定系數(shù)r2=0.8858，相關(guān)系數(shù)r=0.941，相關(guān)性非常好。

圖2 空氣密度標(biāo)準(zhǔn)差與空氣溫度標(biāo)準(zhǔn)差的相關(guān)關(guān)系Fig. 2 Related relationship between air density standard deviation and air temperature standard deviation

對空氣密度標(biāo)準(zhǔn)差和大氣壓力標(biāo)準(zhǔn)差做線性相關(guān)，相關(guān)關(guān)系圖如圖3 所示。相關(guān)方程式為y=0.0128x+0.0364，r2=0.0144，r=0.120，基本不相關(guān)。

圖3 空氣密度標(biāo)準(zhǔn)差與大氣壓力標(biāo)準(zhǔn)差的相關(guān)關(guān)系Fig. 3 Related relationship between air density standard deviation and atmospheric pressure standard deviation

從表1～表3 可以進(jìn)一步統(tǒng)計得到10 個代表區(qū)域測風(fēng)塔測得的空氣密度年平均值為1.093 kg/m3，年平均值標(biāo)準(zhǔn)差為0.04 kg/m3，年平均值標(biāo)準(zhǔn)差占年平均值的3.66%；10 個代表區(qū)域測風(fēng)塔的空氣溫度年平均值為11.2 ℃，年平均值標(biāo)準(zhǔn)差為8.92 ℃，年平均值標(biāo)準(zhǔn)差占年平均值的79.64%；10 個代表區(qū)域測風(fēng)塔的大氣壓力年平均值為89.5 kPa，年平均值標(biāo)準(zhǔn)差為0.32 kPa，年平均值標(biāo)準(zhǔn)差占年平均值的0.36%。通過對年平均值標(biāo)準(zhǔn)差占年平均值的比率進(jìn)行比較分析，發(fā)現(xiàn)大氣壓力的年平均值標(biāo)準(zhǔn)差占年平均值的比率很小，空氣溫度的年平均值標(biāo)準(zhǔn)差占年平均值的比率較大，即大氣壓力的年變化幅度很小，空氣溫度的年變化幅度較大。

從圖2 和圖3 可以分析得出：空氣密度年變化幅度受大氣壓力年變化幅度的影響很小，受空氣溫度年變化幅度的影響較大。因此，四季空氣溫度差異大的區(qū)域，空氣密度的年變化幅度也大；四季空氣溫度差異小的區(qū)域，空氣密度的年變化幅度也小。

2 按某測風(fēng)塔得到的不同值計算發(fā)電量的結(jié)果比較

選取上文10 個區(qū)域中陜西渭南的測風(fēng)塔(記為1#測風(fēng)塔)，分別按空氣密度年平均值(見表1)、空氣密度月平均值(見表1)及空氣密度時間序列值計算發(fā)電量，然后對發(fā)電量計算結(jié)果進(jìn)行比較。由于該測風(fēng)塔所在區(qū)域的空氣密度年變化幅度較大，分析空氣密度年變化情況對風(fēng)電場發(fā)電量計算的影響有一定的典型性。年等效利用小時數(shù)為年發(fā)電量與裝機(jī)容量的比值，因此為方便比較，本文以等效利用小時數(shù)代替發(fā)電量進(jìn)行分析。

1#測風(fēng)塔的塔高為100 m，不同高度處安裝有風(fēng)速傳感器、風(fēng)向傳感器、空氣溫度傳感器、大氣壓力傳感器，傳感器和記錄儀均為美國NRG 品牌的測風(fēng)設(shè)備。

1#測風(fēng)塔的傳感器和記錄儀的配置情況如表4 所示。

表4 1#測風(fēng)塔的傳感器和記錄儀的配置情況Table 4 Configuration of sensors and recorders for 1 # wind measuring tower

2.1 按1#測風(fēng)塔年平均空氣密度計算的等效利用小時數(shù)

根據(jù)1#測風(fēng)塔數(shù)據(jù)得到的年平均空氣溫度為15.6 ℃，年平均大氣壓力為97.5 kPa，使用式(1)，可計算得到該區(qū)域的年平均空氣密度為1.174 kg/m3。

選擇浙江運(yùn)達(dá)風(fēng)電股份有限公司生產(chǎn)的WD156-3000 型號風(fēng)電機(jī)組(下文簡稱為“運(yùn)達(dá)WD156-3000 風(fēng)電機(jī)組”)，輪轂高度為100 m。在年平均空氣密度為1.174 kg/m3、中等湍流強(qiáng)度的功率曲線下，測算測風(fēng)塔單點(diǎn)的年等效利用小時數(shù)，空氣密度折減系數(shù)取78%，則按年平均空氣密度計算得到的年等效利用小時數(shù)為2301.22 h。

2.2 按1#測風(fēng)塔月平均空氣密度計算的等效利用小時數(shù)

統(tǒng)計1#測風(fēng)塔每月的空氣密度，依舊使用運(yùn)達(dá)WD156-3000 風(fēng)電機(jī)組(輪轂高度取100 m)，空氣密度折減系數(shù)取78%，按月平均空氣密度計算1#測風(fēng)塔單點(diǎn)的月等效利用小時數(shù)。

按1#測風(fēng)塔月平均空氣密度和年平均空氣密度分別計算得到的等效利用小時數(shù)結(jié)果如表5所示，月等效利用小時數(shù)對比如圖4 所示。

表5 按1#測風(fēng)塔月平均空氣密度和年平均空氣密度計算的結(jié)果Table 5 Calculated results based on monthly average air density and annual average air density of 1 # anemometer tower

圖4 按年平均空氣密度和月平均空氣密度分別計算得到的月等效利用小時數(shù)對比Fig. 4 Comparison of monthly equivalent utilization hours calculated based on annual average air density and monthly average air density

從表5 和圖4 可以看出：使用月平均空氣密度計算得到的年等效利用小時數(shù)為2293.92 h，使用年平均空氣密度計算得到的年等效利用小時數(shù)為2301.22 h，差7.3 h，使用月平均空氣密度得到的計算結(jié)果比使用年平均空氣密度得到的計算結(jié)果低0.32%，總體差異不大。使用年平均空氣密度和使用月平均空氣密度計算得到的月等效利用小時數(shù)之間的差異最大為3.83%。使用年平均空氣密度來計算月等效利用小時數(shù)時，由于夏季時的空氣密度低于年平均空氣密度，得到的結(jié)果屬于高估；冬季時的空氣密度高于年平均空氣密度，得到的結(jié)果屬于低估。

2.3 按1#測風(fēng)塔空氣密度時間序列值計算的等效利用小時數(shù)

1#測風(fēng)塔記錄儀每10 min 記錄1 組數(shù)據(jù)，包含風(fēng)速、風(fēng)向、空氣溫度、大氣壓力等數(shù)據(jù)。使用1#測風(fēng)塔每10 min 的空氣溫度和大氣壓力，計算每10 min 的空氣密度，精確到0.01 kg/m3，可得到1#測風(fēng)塔全年的空氣密度時間序列值，其中，全年空氣密度時間序列值的最小值為1.06 kg/m3，最大值為1.31 kg/m3。

依然使用運(yùn)達(dá)WD156-3000 風(fēng)電機(jī)組(輪轂高度取100 m)，按照0.01 kg/m3空氣密度間隔的功率曲線、空氣密度折減系數(shù)取78%，計算1#測風(fēng)塔單點(diǎn)的等效利用小時數(shù)序列，并統(tǒng)計年等效利用小時數(shù)。計算得到的年等效利用小時數(shù)為2297.41 h，與使用月平均空氣密度計算得到的年等效利用小時數(shù)(2293.92 h)相比高3.49 h，即高了0.152%；與使用年平均空氣密度計算得到的年等效利用小時數(shù)(2301.22 h)相比低3.81 h，即低了0.166%，差異均較小。

2.4 小結(jié)

在實際應(yīng)用中，如果要考慮風(fēng)電場發(fā)電的月度出力情況，為了精確起見，使用對應(yīng)的月平均空氣密度計算風(fēng)電場發(fā)電量更精確；若只是為了計算風(fēng)電場年發(fā)電量，使用月平均空氣密度和使用年平均空氣密度計算發(fā)電量得到的結(jié)果差異不大。

3 結(jié)論

本文在中國范圍內(nèi)各種地形和氣候條件下，選擇了10 個具有代表性的區(qū)域，對代表性區(qū)域測風(fēng)塔數(shù)據(jù)中的空氣溫度、大氣壓力、空氣密度進(jìn)行分析，并選取其中空氣密度年變化幅度較大的陜西省渭南市的1 座測風(fēng)塔分析空氣密度年變化情況對風(fēng)電場發(fā)電量計算結(jié)果的影響。分析結(jié)果如下：

1)空氣溫度年變化幅度較大的地區(qū)，空氣密度的年變化幅度也較大，空氣密度的月平均值標(biāo)準(zhǔn)差也較大，在0.050 kg/m3以上；而空氣溫度年變化幅度較小的地區(qū)，空氣密度的年變化幅度也很小，空氣密度的月平均值標(biāo)準(zhǔn)差也很小。

2)比較分析年平均值標(biāo)準(zhǔn)差占年平均值的比率，大氣壓力的年平均值標(biāo)準(zhǔn)差占年平均值的比率很小，空氣溫度的年平均值標(biāo)準(zhǔn)差占年平均值的比率較大，即大氣壓力的年變化幅度很小，空氣溫度的年變化幅度較大。

3)空氣密度年變化幅度受大氣壓力年變化幅度的影響很小，受空氣溫度年變化幅度的影響較大。四季空氣溫度差異大的區(qū)域，空氣密度的年變化幅度也大；四季空氣溫度變化小的區(qū)域，空氣密度的年變化幅度也小。

4)計算風(fēng)電場發(fā)電量時，使用月平均空氣密度計算得到的年等效利用小時數(shù)比使用年平均空氣密度計算得到的低0.32%；使用空氣密度時間序列計算得到的年等效利用小時數(shù)比使用月平均空氣密度計算得到的高0.152%，比使用年平均空氣密度計算得到的低0.166%，差異均較小。

使用年平均空氣密度和使用月平均空氣密度計算得到的月等效利用小時數(shù)之間的差異最大為3.83%。使用年平均空氣密度來計算等效利用小時數(shù)時，由于夏季時的空氣密度低于年平均空氣密度，得到的結(jié)果屬于高估；冬季時的空氣密度高于年平均空氣密度，得到的結(jié)果屬于低估。在實際應(yīng)用中，如果要考慮風(fēng)電場發(fā)電的月度出力情況，為了精確起見，使用對應(yīng)的月平均空氣密度計算風(fēng)電場發(fā)電量更精確；若只是為了計算風(fēng)電場年發(fā)電量，使用月平均空氣密度和使用年平均空氣密度計算發(fā)電量得到的結(jié)果差異不大。