長方體建筑屋頂風能利用效能研究

作者簡介：袁行飛(1972)，女，教授，博士生導師，主要從事大跨空間結構研究，（Email）yuanxf@zju.edu.cn。摘要：屋頂安裝小型風力發電機是近幾年來新發展的一種風能利用形式。以年平均風速為對象，研究建筑屋頂風能利用的效能，提出衡量建筑風能利用效能的6個性能指標。利用CFD數值模擬方法分析了長方體建筑平屋頂對風能的利用效果，考察了不同迎風角下控制點的性能指標，探討屋頂風力機安裝的最佳位置和排列形式。對長方體建筑屋頂形式進行分析，對比了4種不同屋頂形式的風能利用效能。

關鍵詞：風能；建筑環境；建筑屋頂；利用效能

中圖分類號：TU023；TK89文獻標志碼：A文章編號：16744764（2012）03011707

Wind Power Utilization Efficiency on Roof of Rectangular Buildings

YUAN Xingfei， ZHANG Yu

(Space Structures Research Center， Zhejiang University， Hangzhou 310058， P. R. China)

Abstract:Installing small wind turbines on the roof of buildings is a new way of wind power utilization in recently years. Based on the annual mean wind speed， the efficiency of wind power utilization on roofs of buildings was studied and the evaluation indexes for the utilization efficiency were proposed. With the help of CFD numerical analysis， the analysis on wind power utilization efficiency on flat roofs of rectangular buildings was carried out. The performance indexes of reference points under different windward angles were investigated. The best position and arrangement of wind turbines on the roof were also discussed. Furthermore， the utilization efficiency of four roofs with different shapes were analyzed and compared with that of flat roof.

Key words:wind power; building environment; building roofs; utilization efficiency

大力開發風能是解決中國能源短缺、實現可持續發展戰略的重要途徑之一［1］。建筑環境中的風能利用具有免于輸送的優點，所產生的電能可直接用于建筑本身，為綠色建筑的發展提供了一種新的思路，對實現建筑可持續化發展、緩解能源與環境矛盾具有重要而深遠的意義。

近年來一些學者對建筑環境中的風能利用技術進行了探索。研究主要著眼于可行性分析［24］、建筑風環境模擬［57］、建筑風力集中器設計［89］以及適宜建筑環境的風力發電機開發［1011］等方面。與沙漠、曠野、近海相比，建筑環境中的風場相對復雜、紊流加劇、風速降低，但局部區域也存在較大風速，如高層建筑頂部的小急流、建筑物開洞部位的穿堂風、相鄰建筑通道的夾道效應等。為提高風能利用效率，風力機通常安裝在上述部位。根據建筑中風力機的放置位置，風能利用建筑大致可分為：頂部安裝型、空洞安裝型和通道安裝型。其中頂部安裝型應用較為廣泛。本文采用CFD數值模擬方法分析了長方體平屋頂建筑的風能利用效能，研究了不同迎風角對建筑風能集結效果的影響，探討了風力機的最佳安裝位置，并對屋頂形式進行了優化分析，提出了有利于風能集結效果的屋頂形式，為實際工程應用提供理論基礎和技術支撐。1建筑風能利用效能評價指標

風力機的功率與風速三次方成正比，因此風力機的安裝位置應選取風速相對較大的位置，同時避免高湍流。為降低成本，還應盡量降低風力機的安裝高度。由此提出以下評價屋頂風能利用效能的指標：〖=D(〗袁行飛，等：長方體建筑屋頂風能利用效能研究〖=〗

1）風速增大系數Cv=V/V0－1，衡量建筑風能集結效果最主要的參數，其中V為高度z處的實際風速，V0為高度z處未受建筑撓動時的風速。Cv越大，表明建筑對風能的集結效果越好，N個位置的總風速增大系數為Cv總=ΣNi=1Cvi。

2）湍流強度I，衡量建筑風場的湍流程度。風場的湍流會減少風機輸出功率，引起極端荷載，最終削弱破壞風力機，因此風力機的安裝位置應避免高湍流。

3）屋頂紊流厚度δ，是屋頂風速變化梯度較大區域和高湍流區域兩者高度的最大值，用來評估風力機的安裝高度。

4）實際風速V，衡量具體位置的風能情況，評價其風能利用的可行性，同時為適宜風力機的選擇提供依據。

5）風速傾斜角γ，水平軸風力機尚不能利用垂直方向的風速進行發電；對于垂直軸風力機，風速傾斜角在一定程度上可提高風力機的功率［12］。

6）行人高度處風速V2，max，衡量建筑風環境的舒適度。為滿足建筑風環境舒適度，要求V2，max≤5 m/s。2建筑風環境數值模擬假定和方法

2.1數值模擬基本假定

1）流體不可壓縮；2)流體為Newton流體，忽略粘性耗散；3)流體在固壁上無滑移；4)流體各向同性；5)流動為穩態。

2.2計算模型和方法

1）幾何建模及網格劃分

幾何建模和網格劃分在FLUENT 6.3前處理軟件Gambit中完成。計算流域如圖1， H為建筑物的高度。采用四面體單元進行非結構化網格離散，網格尺寸由內往外逐步增大。

2）湍流模型

考慮到在鈍體繞流中的精度，采用可實現的kε湍流模型（Realizable kε）。采用非平衡的壁面函數法處理近壁面的湍流狀態。該法能考慮壓力梯度和偏離平衡假設帶來的影響，對具有像環繞、分離、再附、撞擊等復雜流動的計算具有較高的精度。

圖1計算流域示意圖

3）邊界條件

進流面采用速度入口邊界，出口采用壓力出口邊界，計算區域上壁面及兩側壁面采用對稱邊界，相當于設置光滑壁面，其剪切應力為0；建筑物表面和地面采用無滑移的壁面條件（wall），沿壁面切向流體速度為0。

入口邊界采用UDF接口輸入，風速采用滿足指數律的風剖面函數：高度Z處的風速U－(z)=U10－(Z10)a；湍動能k=32(U(z)－×I)2；耗散率ε=0.090.75×k3/2/l，湍流強度Iuz=I0ZHT－α－0.05，Z=30 m時的湍流強度 Iu(30)可由地貌系數α得出：A類地貌的I0=0.081 1，B類地貌的I0=0095 5，C類地貌的I0=0.109 3，D類地貌的I0=0.116 3；湍流積分長度l=100×(Z/30)0.5。取C類地貌，U－10=3 m/s，α=0.22。

4）計算方法

計算采用3D單精度，分離式求解器，空氣模型選用不可壓縮的常密度空氣模型，對流項的離散采用精度較高的二階迎風格式，速度壓力耦合采用SIMPLEC算法。計算中使用自適應網格技術（Adapted Grid），每步迭代均對壓力梯度大的區域進行網格加密，以便更準確地預測復雜的流動。3長方體平屋頂建筑風能利用效果及風力機最佳位置3.1迎風角φ=00時屋頂風能利用效果

長方體平屋頂建筑基本模型：長L=50 m，寬W=20 m，高H=30 m，如圖2，h為屋頂以上高度，由于模型的對稱性，取屋頂前沿中間和后沿各2個參考位置，分別為A（-W/2， 0， H+h）、D（W/2，-L/2， H+h）、B（0， 0， H+h）、E（0，-L/2， H+h）、C（W/2， 0， H+h）、F（W/2，-L/2， H+h），分析屋頂的風能集結效果。

圖2長方體平屋頂建筑模型通過CFD數值模擬，可得出各指標參數的變化情況。圖3為建筑物中間豎向剖面（y=0 m）和邊線豎向剖面（y=-25 m）的風速V等值線圖。圖4為屋頂風速增大系數Cv隨屋頂以上高度h的變化情況。

圖3建筑豎向剖面風速等值線圖

由圖3可見：風速V在屋頂附近的變化梯度較大，且隨h的增加而增大，并趨于穩定。當h>3 m時，前沿點的風速增大梯度趨于穩定；當h>8 m時，中線點的風速增大梯度趨于穩定；當h>12 m時，后沿點的風速增大梯度趨于穩定。屋頂湍流厚度δ沿來流方向呈增厚趨勢。屋頂湍流強度隨h的增大先增大后減小，在h=5 m內較大，高湍流區出現在建筑屋頂的中前方。

圖4φ=0時風速增大系數Cv

由圖4可見，前沿A點的風速增大系數在h=5 m處達到最大值0.085，對應的實際風速為4.29 m/s，湍流強度為0.22；D點在h=3.5 m處達到最大值0079，對應的實際風速為4.24 m/s，湍流強度為023。A點的風速傾斜角最大值可達70°，D點的風速傾斜角最大值可達46°；在h=5 m處，A點的風速傾斜角仍有28°，D點為17°。中線B點的風速增大系數在h=11 m處達到最大值0.093，對應的實際風速為4.47 m/s，湍流強度為0.18；E點在h=10 m處達到最大值0.075，對應的實際風速為4.37 m/s，湍流強度0.16。后沿C點的風速增大系數在h=19.5 m處達到最大值0.060，對應的實際風速為451 m/s，湍流強度為0.13；F點在h=14.5 m處達到最大值0.055，對應的實際風速為4.39 m/s，湍流強度為0.14。中線點和后沿點的風速傾斜角均在0°附近波動。

綜上所述，沿來流風方向，屋頂紊流厚度逐漸增厚；風速增大系數、湍流強度和實際風速隨h的變化均是前沿點先達到最大值，接著是中線點，之后是后沿點；各參考位置風速增大系數最值所對應的風速均大于3 m/s，處于有效風速范圍內；前沿點的風向變化比較大，中線點和后沿點的風向變化不大。

3.2不同迎風角下屋頂風能利用效果

考慮建筑對稱性，分別在建筑屋頂的前沿、中線和后沿增加一個參考位置：D0、E0、F0與位置D、E、F對稱。由分析知，為提高風力機的發電效能，風力機應在年主風向垂直的方向上盡量多的排放。假定風力機沿建筑長度方向排布，分別分析h=5 m處前沿點（A、D、D0），h=10 m處中線點（B、E、E0）和h=15 m處后沿點（C、F、F0）在不同來流風向（000、2250、4500、67.50、9000）的風能利用效果。分析結果見表1。表1不同迎風角下各個參考點的性能參數

迎風角前沿(h=5 m)中線(h=10 m)后沿(h=15 m)ADD0BEE0CFF0V2，max/

（m·s-1）0.00Cv0.0850.0720.0720.0860.0740.0740.0480.0540.054V/（m·s-1）4.2904.2404.2404.4204.3704.3704.3804.4004.400I0.2200.2000.2000.1900.1600.1600.1500.1300.130cosγ0.8770.9530.9530.9840.9940.9940.9960.9980.996ΣCv0.2290.2340.1563.622.50Cv0.0510.038-0.0130.0650.061-0.0080.0470.0530.013V/（m·s-1）4.1504.1003.9004.3404.3204.0404.3704.4104.230I0.2200.1600.1900.1800.1300.1800.1500.1200.150cosγ0.9230.9710.9210.9900.9970.9730.9980.9990.992ΣCv0.0760.1180.1133.3450Cv0.0640.018-0.0040.0460.061-0.0010.0540.0540.041V/（m·s-1）4.2004.0203.9404.2604.3104.0704.4004.4004.350I0.240.1700.1700.1900.1600.2000.1600.1500.160cosγ0.9590.9761.0000.9990.9930.9941.0000.9970.998ΣCv0.0780.1060.1493.867.50Cv0.001-0.005-0.0650.0290.012-0.0320.0390.0130.008V/（m·s-1）3.9603.9303.7004.1904.1203.9404.3404.2304.210I0.1500.1400.1400.1400.1300.1300.1300.1200.120cosγ0.9840.9780.9821.0000.9861.0001.0000.9940.999ΣCv-0.0690.0100.0603.5900Cv-0.1500.053-0.229-0.0830.028-0.175-0.0020.013-0.045V/（m·s-1）3.3604.1603.0503.7304.1803.3604.1704.2303.990I0.1800.2100.1600.1800.1800.1800.1600.1500.150cosγ1.0000.9830.9991.0000.9871.0001.0000.9961.000ΣCv-0.326-0.230-0.0343.5

通過分析可見，在（0.00，450）之間，前沿，中線和后沿的風速增大系數合值均隨著φ的增大先減小后增大，最大值出現迎風角φ=0.00時；在（45.00，90.00）之間，風速增大系數合值急劇減小；當φ=900時，達到最小值。

3.3屋頂風力機安裝位置

屋頂風力機的安裝位置應選取風速相對較大，且變化較平穩的位置，以減小風力機葉輪的脈動，同時避免高湍流。屋頂安裝風力機的高度宜滿足h≥δ+D/2，D為水平軸風力機的直徑，對于垂直軸風力機的安裝高度宜滿足h≥δ。

由前述分析知，對于平屋頂不同位置來說，前沿點的風能集結效果最好——風速比達到最值時的高度最低，且風速變化梯度趨于平穩時的高度也最低，同時避開了高湍流；由于建筑屋面對風速的阻礙作用，建筑物屋頂上的風速在同一高度下，由前沿到后沿呈下降趨勢，隨著高度的增加，同一高度處，屋面的風速趨于一致。

建筑物屋頂安裝風力機的紊流，除了來自建筑物的影響之外，還有一個重要的因素：風力機尾流的相互影響。為了盡量減少風力機尾流的相互影響，風力機的排列間距應滿足一定的要求［12］。因此風力機應在年主風向垂直的方向上盡量多的排放，對于安裝風力機的建筑來說，應盡可能使建筑的長度方向垂直于年主風向方向，即來流風向控制在（-45°，45°），以提高風力機的發電效能。

假設年主風向垂直于建筑長度方向，當風機尺寸D>W/8時，風力機沿前沿單排布置的風能利用效率最佳，考慮到風速傾斜角的問題，應優先選擇升力型垂直軸風力機［12］，當安裝水平軸風力機時，可傾斜一定角度，安裝高度大于5 m；亦可沿中線單排布置，安裝高度大于10 m，其利用效能較好，還可減少結構不均勻受力；當風機尺寸D＜W/8，風力機沿寬度方向可多排布置，其中前后排有高差排列的風能利用效能較前沿單排布置和中線單排布置方式好，既可以有效地利用風能，降低安裝高度，又可以減小風力機之間的尾流影響。4長方體建筑不同屋頂形式分析

4.1不同屋頂形式的幾何模型

目前，在將風力機引入建筑物的技術中，平屋頂上直接安裝風力機是一種最直接最易實現的方式，因此本文在分析其對風能集結效果的基礎上，將其作為分析其他屋頂形式風能集結效果的參考依據，如圖7為4種不同的屋頂形式，建筑物的長度L、寬度W和總高H均不變。其中圖7(a)為前高單坡型屋頂；圖7(b)為后高單坡型屋頂；圖7(c)為圓弧雙坡型屋頂；圖7(d)為梯形雙坡型屋頂，屋頂面的寬度l=5 m。各屋頂形式的屋頂傾斜高度分別用Ma、Mb、Mc、Md表示。取各屋頂形式的屋頂傾斜高度分別為2.5、5.0、7.5、10.0 m與平屋頂（M=0）進行比較分析。

圖7屋頂不同形式

4.2不同形式屋頂的風能利用效能

分析迎風角φ=0時，各屋頂形式的風能利用效能。圖8為各屋頂形式屋頂傾斜高度等于5 m時，建筑中間豎直剖面的風速V等值線圖和湍流強度I等值線圖。由圖8可見，前高單坡型屋頂前沿點的紊流厚度與平屋頂的基本相同，中線點和后沿點的紊流厚度有所增大，同一屋頂高度處，風速有所減小；后高單坡型屋頂前沿點的紊流厚度較平屋頂的高，中線點尤其是后沿點的紊流厚度有所減小；雙坡型屋頂中線點的紊流厚度較平屋頂有所減小，屋頂以上同一h下，圓弧雙坡型屋頂中線點的風速較平屋頂的大。

圖8各屋頂形式建筑y=0 m豎直剖面風速等值線圖(m/s)

圖9為不同屋頂傾斜高度下，各形式屋頂的風速增大系數Cv隨屋頂以上高度h的變化情況。相比平屋頂形式，以上4種屋頂形式，只有圓弧雙坡型屋頂和梯形雙坡型屋頂的風速增大系數最值有所提高，其中圓弧屋頂提高的幅度最大。隨著屋頂傾斜高度的增加，圓弧屋頂的風速增大系數最值越來越大，且達到最值時所對應的高度h也越來越小；其他形式屋頂的風速增大系數最值越來越小，但最值所對應的高度h越來越小。

綜上，當迎風角φ=0時，前單坡屋頂適宜沿前沿單排安放風力機，中線和后沿的紊流厚度較大，不宜安放風力機；雖然后單坡屋頂的風速增大系數較平屋頂有所減小，但中線和后沿的紊流厚度明顯減

圖9不同屋頂高度下各屋頂形式特征位置風速

增大系數Cv變化情況

小，風力機沿前沿、中線和后沿安放的利用效能相當，當建筑寬度允許的情況下，可以利用建筑本身前后的高差實現風力機的多排布置；圓弧雙屋頂和梯形雙屋頂沿中線安放風力機的風能利用效能較平屋頂有所提高，由于只可單排安放風力機，適宜建筑寬度較小的建筑。5結論

通過CFD數值模擬分析了長方體建筑屋頂風能利用效果，比較了不同迎風角、不同屋頂形式對風能集結效果的影響，得出以下結論：

1）長方體平屋頂建筑的屋頂對風能有一定的集結效果，其集結效果受來流風向影響較大。

2）長方體平屋頂建筑沿屋頂前沿點安裝風力機的風能利用效果最佳——相同風速下，其安裝高度最低，同時應將來流風向角控制在（-450，450）的范圍內。當建筑寬度較大時，長方體平屋頂建筑屋頂可前后有高差地多排安放風力機，既減小風力機尾流的影響，同時也可增加屋頂風力機的安裝數目，提高建筑屋頂的總發電量。

3）幾種不同屋頂形式下，后單坡屋頂宜發展前后有高差多排布置風力機；當屋頂單排布置風力機時，梯形雙坡屋頂和圓弧雙坡屋頂沿中線單排布置風力機的效能最佳。

4）通過分析得知建筑屋頂對風能具有一定的集結效果，但建筑屋頂風場比較復雜，紊流很大，對風力機性能的影響不容忽視，風力機能否在建筑中得到廣泛應用，還需進一步的研究。

參考文獻：

［1］石定環. 可再生能源與可持續發展［J］.中國科技產業， 2008(1):1518.

SHI Dinghua. Renewable energy and sustainable development［J］.Science Technology Industry of China， 2008(1): 1518.

［2］Taylar D.Using buildings to harvest wind energy［J］.Building Research Information， 1998， 26(3): 199202.

［3］Mertens S. Wind energy in urban areas: concentrator effects for wind turbines close to buildings［J］. Refocus， 2002， 3(2): 2224.

［4］李秋勝， 李永貴， 陳伏彬， 等. 風能發電在超高層建筑中的應用研究［C］//第十四屆全國風工程學術會議論文集， 2009:801805.

LI Qiusheng， LI Yonggui， CHEN Fubin， et al. Application research of energy generation by wind on highrise buildings［C］//The 14th proceeding of national wind engineering， 2009， 801805.

［5］Lin L U， Yan K I P．Investigation on the feasibility and enhancement methods of wind power utilization in highrise buildings of Hong Kong［J］．Renewable Sustainable Energy Reviews， 2009， 13 (2):450461.

［6］姜瑜君，桑建國，張伯寅.高層建筑的風環境評估［J］. 北京大學學報:自然科學版，2006，42(1)：6873.

JIANG Yujun， SANG Jianguo， ZHANG Boyin. An evaluation on wind environment around skyscrapers［J］. Acta Scientiarum Naturalium Universitatis Pekinensis， 2006， 42(1):6873.

［7］蘇國， 陳水福. 復雜體型高層建筑表面風壓周圍風環境數值模擬［J］. 工程力學， 2006， 23(8): 152157.

SU Guo， CHEN Shuifu.Numerical simulation of wind pressures and wind environment around a complex shaped highrise building［J］. Engineering Mechanics， 2006， 23(8): 144149.

［8］Mertens S. The energy yield of roof mounted wind turbines［J］.Wind Engineering， 2003， 27(6):507517.

［9］潘雷，陳寶明.建筑環境中風能利用的研究［D］. 濟南: 山東建筑大學，2006.

［10］Holdsworth B. Options for microwind generation: part1［J］. Renewable Energy Focus， 2009，10 (2): 6063.

［11］王乃粒.城市風力發電［J］．世界科學，2005(1):3940.

WANG Naili.Urban wind power generation［J］. World Science， 2005 (1):3940.

［12］陳坤，賀德馨.風力機尾流數學模型及尾流對風力機性能的影響研究［J］. 流體力學實驗與測量， 2003，17(1):8487.

CHEN Kun， HE Dexin.The study for the wake model of wind turbine and the wake effects on the performance of wind turbine［J］. Experiments and Measurements in Fluid Mechanics， 2003， 17(1):8487.

（編輯王秀玲）