隨著科技的發展、醫學的進步，人口有逐漸成長的現象，但需要在有限的土地上求生活空間，所以建築物有向上發展的趨勢，加上政府近年來推動都市更新的計畫，都市中存在著越來越多的高層建築。高層建築除了受到地震力的影響，對風力的敏感程度也會增加，因此耐風設計對結構物的反應是高層建築設計裡重要的一環。在風工程研究中，無論是風洞實驗、數值模擬，最主要的目的都是在實場上的耐風設計及風環境評估，例如相鄰建築物間風力的相互作用，與建築物周遭環境的行人風場，皆與生活品質與居住安全有著密不可分的關係，而台灣地區颱風所造成之財物損失，位居台灣所有天然災害的第一位，解決的方法唯有對風場行為有完整的了解，才能以工程的手段來降低風災對工程設施的破壞，進而提高人類居住環境的品質。因此，風工程儼然成為一門重要的課題。但實場量測儀器昂貴，所量測數據非常有限，因此可以利用風洞實驗所做出來的結果去補足在實場中無法量測到的資料數據，做為設計及評估上的參考。本研究以風洞模擬實驗為主分作兩個部分，一為符合標準地況自然風場之風洞模擬，在實場量測中我們不易得知大氣邊界層真實的紊流特性，透過風洞實驗的結果，我們能更了解該地況的邊界層特性。但如果要進行一個真實地況的大氣邊界層模擬，在製作地況縮尺模型上需要耗費大量的人力、物力。本研究希望能將真實地物，以一簡單的方法、公式模擬成均勻方塊的粗
糙元素，代替原實場地形來進行風洞實驗，以節省許多的經費與時間，此為
符合標準地況自然風場之風洞模擬。二為複雜地形自然風場之風洞模擬，由於台灣地窄人稠，各地形皆相鄰，許多建物皆建築在山坡、山腳上，大氣邊界層特性實為複雜，亦與各國不同，本研究最終以台灣北部淡水北投一帶大屯山系山腳一處俗稱小坪頂的山區平地，為綜合型的複雜地形場址，加上前述的符合標準地況自然風場之風洞模擬，以進行了解台灣本地的大氣邊界層特性，進而建立出適合台灣本地的耐風設計規範及風洞實驗技術。

The simulation of turbulence boundary layers in this study can be divided into two parts. The first part is simulation of flow field developed over terrains that compliance with the definitions of the terrain categories in the building wind code. Since current field measurements can only provide reliable mean wind speed profile. Properly simulated turbulence boundary layer is the best and only viable means to acquire turbulence features of the atmospheric boundary layers. Besides simulating flow fields via scaled terrain models, the thesis further studied the procedure to develop equivalent floor roughness elements to generate turbulent boundary layers over urban, suburban and open terrains. The second part of this thesis is the wind tunnel simulation of wind field over complex terrain. A 1/5000 mountainous terrain model centered at Xiaopingding, tamsui in the northeastern corner of Taiwan was built. Wind tunnel measurements were then compared with field measurement of Lidar and numerical simulation by LES. The initial result indicates that wind tunnel simulation of rugged terrain has its limitations. In real life, hills and mountains are not homogenous terrains. The shape and geographic lay out of high mountains would significantly distort the approaching flow. This geographic effect can be extremely difficult properly simulated in a fixed sized wind tunnel test section. More works need to be done in this subject.

目錄
第1章、 緒論	1
1.1. 前言	1
1.2. 研究動機與目的	1
1.3. 研究方法	2
1.4. 研究流程	3
1.5. 論文架構	4
第2章、 文獻回顧	5
2.1. 大氣邊界層流場之風洞模擬	5
2.1.1. 雷諾數效應	6
2.1.2. 風洞實驗之阻塞效應	7
2.2. 粗糙元素	8
2.3. 粗糙元素的間隔對流場的影響	8
2.4. 流場發展長度的限制	10
2.5. 粗糙長度	10
2.5.1. 實場對於粗糙長度的估算公式	11
2.5.2. 粗糙元素對於粗糙長度的估算公式	11
2.6. 零平面高度	13
2.7. 粗糙長度與α值的關係	14
2.8. 本國規範	15
第3章、 理論背景	16
3.1. 大氣邊界層紊流特性	16
3.2. 風速剖面	17
3.3. 紊流強度	19
3.4. 紊流長度尺度	20
3.5. 紊流頻譜	21
3.6. 空間相關函數	23
3.7. 陣風因子	24
第4章、 實驗儀器與量測方法	25
4.1. 風洞實驗室	25
4.2. 大氣邊界層流場之模擬	26
4.3. 風洞實驗測量儀器	27
4.3.1. 測量儀器-熱膜探針	27
4.3.2. 測量儀器- 4孔式壓力感測器	28
4.4. 實驗數據及採樣分析	29
4.5. 可攜式風速遙測系統-光達	31
4.5.1. 光達系統	31
4.5.2. 數據採樣	32
第5章、 前期實驗研究	34
5.1. 自然風場之實場量測	34
5.1.1. 符合標準地況自然風場之實場量測	34
5.1.2. 複雜地形自然風場之實場量測	40
5.2. 大氣邊界層之風洞模擬	43
5.2.1. 都市地況測點之風洞縮尺模擬	43
5.2.2. 設計粗糙元素	43
5.2.3. 鄉鎮地況與開闊地況測點	46
5.3. 複雜地形之數值模擬	49
5.3.1. 地形數值資料	49
5.3.2. 複雜地形繪製	50
5.3.3. 數值模擬	51
5.3.4. 風洞試驗與數值模擬比較	53
5.3.5. 後續研究	54
第6章、 符合標準地況自然風場之風洞模擬	55
6.1. 實驗設計	55
6.2. 都市地況	58
6.2.1. 量測內容	58
6.2.2. 模擬修正及結果數據比較	61
6.3. 鄉鎮地況	64
6.3.1. 量測內容	64
6.3.2. 模擬修正及結果數據比較	68
6.4. 開闊地況	71
6.4.1. 量測內容	71
6.4.2. 模擬修正及結果數據比較	75
6.5. 小結	78
第7章、 複雜地形自然風場之風洞模擬	79
7.1. 第一階段	81
7.1.1. 實驗設置	81
7.1.2. 實驗結果及數據比對	82
7.2. 第二階段	84
7.2.1. 實驗設置	84
7.2.2. 實驗結果及數據比對	87
7.3. 第三階段	95
7.3.1. 實驗設置	95
7.3.2. 實驗結果及數據比對	96
7.4. 討論	98
7.5. 小結	98
第8章、 結論與建議	99
8.1. 結論	99
8.2. 建議	101
參考文獻	102
附錄A	106
附錄B	111
附錄C	114
附錄D	115

 
圖目錄
圖 2-1 Isolated roughness flow	8
圖 2-2 Wake interference flow	8
圖 2-3 Skimming flow	8
圖 2-4 順風向紊流強度與表面積密度關係圖	9
圖 2-5 順風向紊流強度剖面	9
圖 2-6 發展長度示意圖	10
圖 2-7 粗糙長度（z0）與α值的關係圖	14
圖 3-1 邊界成形成示意圖（1）	16
圖 3-2 邊界成形成示意圖（2）	16
圖 3-3 m與z0關係圖	21
圖 4-1 淡江大學一、二號風洞實驗室平、立面圖	25
圖 4-2 擾流板示意圖	26
圖 4-3 決定擾流板之高度與寬度之經驗曲線	26
圖 4-4 Cobra probe全體圖	28
圖 4-5 Cobra probe探頭風向分量說明	29
圖 4-6 Cobra probe探頭風向角度說明	29
圖 4-7 不同類型的光達掃描策略	31
圖 4-8 VAD法之示意圖	32
圖 4-9 宜蘭光復國小十分鐘平均風速剖面圖	33
圖 4-10 宜蘭新南國小十分鐘平均風速剖面圖	33
圖 5-1 A地況量測位置	35
圖 5-2 B地況量測位置	36
圖 5-3 C地況量測位置	36
圖 5-4 A地況每小時平均所得之風速剖面（a）季風（b）颱風	38
圖 5-5 B地況每小時平均所得之風速剖面（a）季風（b）颱風	39
圖 5-6 C地況每小時平均所得之風速剖面（a）季風（b）颱風	39
圖 5-7 複雜地形量測位置	40
圖 5-8 複雜地形東北季風SNR訊噪比	41
圖 5-9 複雜地形東北季風 風速剖面	41
圖 5-10 複雜地形東北季風 風速歷時	42
圖 5-11 複雜地形東北季風 風向歷時	42
圖 5-12 都市地況測點之平均風速、紊流強度與紊流尺度剖面	44
圖 5-13 都市地況測點之平均風速、紊流強度與紊流尺度剖面	45
圖 5-14 均勻粗糙元素模擬之風速頻譜 (a)100m (b)240m (c)450m)	46
圖 5-15 鄉鎮地況測點之平均風速、紊流強度與紊流尺度剖面	47
圖 5-16 鄉鎮地況測點之平均風速、紊流強度與紊流尺度剖面	47
圖 5-17 開闊地況測點之平均風速、紊流強度與紊流尺度剖面	48
圖 5-18 淡水小坪頂位置	49
圖 5-19 網格繪製示意圖	51
圖 5-20 風洞試驗與數值模擬 (a)西北風 (b)東北風	53
圖 5-21 逼近流之風速與紊流強度剖面	54
圖 6-1 粗糙元素轉換、符號示意圖	57
圖 6-2 都市地況千分之一數值地形圖 模擬範圍	58
圖 6-3 都市地況縮尺模型	58
圖 6-4 擾流板尺寸配置圖（一）	59
圖 6-5 擾流板尺寸配置圖（二）	59
圖 6-6 擾流板尺寸配置圖（三）	59
圖 6-7 無因次化風速剖面比對	59
圖 6-8 都市地況均勻粗糙元素	61
圖 6-9 都市地況之平均風速、紊流強度與紊流尺度剖面	62
圖 6-10 都市地況之平均風速、紊流強度與紊流尺度剖面	62
圖 6-11 都市地況縮尺模擬之風速頻譜 (a)z=90m (b)285m (c)480m	63
圖 6-12 均勻粗糙元素模擬之風速頻譜 (a)z=90m (b)285m (c)480m	63
圖 6-13 都市地況測點之空間相關性	63
圖 6-14 鄉鎮地況千分之一數值地形圖 模擬範圍	64
圖 6-15 鄉鎮地況縮尺模型	64
圖 6-16 擾流板尺寸配置（一）	65
圖 6-17 擾流板尺寸配置（二）	65
圖 6-18 擾流板尺寸配置（三）	66
圖 6-19 無因次化風速剖面比對	66
圖 6-20 紊流尺度剖面比對	66
圖 6-21 鄉鎮地況均勻粗糙元素	68
圖 6-22 鄉鎮地況之平均風速、紊流強度與紊流尺度剖面	69
圖 6-23 鄉鎮地況之平均風速、紊流強度與紊流尺度剖面	69
圖 6-24 鄉鎮地況縮尺模擬之風速頻譜 (a)z=60m (b)210m (c)360m	70
圖 6-25 均勻粗糙元素模擬之風速頻譜 (a)z=60m (b)210m (c)360m	70
圖 6-26 鄉鎮地況測點之空間相關性	70
圖 6-27 開闊地況千分之一數值地形圖 模擬範圍	71
圖 6-28 開闊地況縮尺模型	71
圖 6-29 擾流板尺寸配置（一）	72
圖 6-30 擾流板尺寸配置（二）	72
圖 6-31 擾流板尺寸配置（三）	73
圖 6-32 無因次化風速剖面比對	73
圖 6-33 紊流強度剖面比對	73
圖 6-34 開闊地況均勻粗糙元素	74
圖 6-35 開闊地況之平均風速、紊流強度與紊流尺度剖面	76
圖 6-36 開闊地況之平均風速、紊流強度與紊流尺度剖面	76
圖 6-37 縮尺模擬(有河堤)之風速頻譜 (a)z=60m (b)165m (c)270m	76
圖 6-38 縮尺模擬(無河堤)之風速頻譜 (a)z=60m (b)165m (c)270m	77
圖 6-39 均勻粗糙元素模擬之風速頻譜 (a)z=60m (b)165m (c)270m	77
圖 6-40 開闊地況之空間相關性	77
圖 7-1 複雜地形量測位置地形圖	79
圖 7-2 複雜地形模擬範圍、測點及大屯山系相關位置示意圖	80
圖 7-3 複雜地形1/5000縮尺模型	81
圖 7-4 複雜地形剖面圖	81
圖 7-5  2003年~2013年 淡水氣象站風向風花圖	82
圖 7-6 東北來風向 風速剖面、紊流強度、紊流尺度	82
圖 7-7 東北來風向 總體紊流強度、風向角、風攻角剖面	82
圖 7-8 西北來風向 風速剖面、紊流強度、紊流尺度	83
圖 7-9 西北來風向 總體紊流強度、風向角、風攻角剖面	83
圖 7-10 西北風、東北風 風速剖面比對	83
圖 7-11 延伸地形、過渡段配置示意圖	84
圖 7-12 北北東（22.5°）過渡段實際安裝示意圖	85
圖 7-13 鼓浪嶼氣象站近十年各月風花圖	85
圖 7-14 東北（45°）地表粗糙度改善示意圖	86
圖 7-15 配置一、二 風速、紊流強度、紊流尺度、總體紊流強度、風向角、風攻角剖面	87
圖 7-16 配置三、四 風速、紊流強度、紊流尺度、總體紊流強度、風向角、風攻角剖面	88
圖 7-17 配置一、三 風速、紊流強度、紊流尺度、總體紊流強度、風向角、風攻角剖面	89
圖 7-18 配置二、四 風速、紊流強度、紊流尺度、總體紊流強度、風向角、風攻角剖面	90
圖 7-19 配置二、五 風速、紊流強度、紊流尺度、總體紊流強度、風向角、風攻角剖面	91
圖 7-20 東北風 風洞實驗與數值模擬風速剖面比對	92
圖 7-21 空風洞風速發展剖面	93
圖 7-22 空風洞A點無側板 風速、紊流強度、紊流尺度 風向角、風攻角剖面	93
圖 7-23 空風洞B點無側板 風速、紊流強度、紊流尺度 風向角、風攻角剖面	94
圖 7-24 空風洞B點模型設置示意圖	95
圖 7-25 地形模型加6cm端版設置示意圖	95
圖 7-26 端版前緣無因次化風速剖面圖 （●：距前緣40 cm  ○：60 cm  ▲：80 cm  △：100 cm）	96
圖 7-27 端版前緣總體紊流強度、紊流尺度、風向角、風攻角剖面 （●：距前緣40 cm  ○：60 cm  ▲：80 cm  △：100 cm）	96
圖 7-28 第三階段 風速、紊流強度、紊流尺度 風向角、風攻角剖面	97
 
表目錄
表 3-1 Simiu 粗糙長度與地形對照表	17
表 3-2 John D.Holmes 粗糙長度與地形對照表	17
表 3-3 α、δ與地形分類對照表	18
表 3-4 地表粗糙長度與β關係表	19
表 3-5 同時間長度之c（t）值	24
表 5-1 實場監測所得之平均風速剖面	39
表 5-2 淡江大學台北校園測點風場模擬之地表粗糙元素	44
表 6-1 都市地況之各項參數資料修正	61
表 6-2 鄉鎮地況之各項參數資料修正	68
表 6-3 開闊地況之各項參數資料修正	75

[1]Armitt, J. & Counihan, J., “The Simulation of the Atmospheric Environment”, 2 (1968) 49-71.
[2]Counihan, J., “An Improved Method of Simulation Atmospheric Boundary Layer”, Atmospheric Environment, 4 (1970) 159-275.
[3]Counihan, J., “Further Measurements in a Simulated Atmospheric Boundary Layer”﹐Atmospheric Environment, 4 (1970) 159-275. 
[4]Counihan, J., “Simulation of an Adiabatic Urban Boundary Layer in a Wind Tunnel”﹐Atmospheric Environment, 7 (1973) 673-689. 
[5]Standen, N. M., “A Spire Array for Generating Thick Turbulent Shear Layers for Natural Wind Simulation in Wind Tunnels”, Rep. LTR-LA-94, National Aeronautical Establishment, Ottawa, Canada, 1972.
[6]Barret, R. V., “A Versatile Compact Wind Tunnel for Industrial Aerodynamics”, Technical note, Atmospheric Environment, 6 (1972) 491-495.
[7]Cook, N. J., “On Simulating the lower Third of the Urban Adiabatic Boundary Layer in a Wind Tunnel”, Atmospheric Environment, 7 (1973) 691-705.
[8]Cermak, J. E. & Peterka, J. A.﹐1974﹐ “Simulation of Atmospheric Flows in Short Wind Tunnel Test Sections”, Center for Building Technology, IAT, National Bureau of Standards Washington, D.C., June. 1974.
[9]Simiu, E., Scanlan, R. H., “Wind Effects on Structures”, 2nd edit., John Wiley & Sons, 1986.
[10]Nakamura, Y. & Ohya, Y., “The effects of turbulence on the mean flow past two dimensional rectangular cylinders”, J. of Fluid. Mech., 149 (1984) 255-273.
[11]Hoerner, S. F., “Fluid Dynamic Drag”, published by the author, 148 Busteed Drive. Midland Park. N. J., 1965.
[12]Townsend, A. A., “The Structure of Turbulent Shear Flow”, Cambridge Univ. Press, 1956.
[13]Jensen, M., 1958, “The Model Law for Phenomena in Natural Wind”, Ingeioen International Edition, 4 (1958) 121-123.
[14]Whitbread, R. E., “Model Simulation of Wind Effects on Structures”, Proceeding of the Conference on Wind Effects on Buildings and Structures, (1963) 284-306.
[15]Biggs, J. M., “Wind Load on Truss Bridges”, ASCE, (1954) 879.
[16]Hunt, A., “Wind Tunnel Measurement of Surface Pressure on Cubic Building Models at Several Scales”, Journal of Wind Engineering and Industrial Aerodynamics, 10 (1982) 137-163.
[17]Wieringa, J., “Representative roughness parameters for homogeneous terrain”, Boundary Layer Meteorology, 63 (1993) 323-363.
[18]Yunqiu Jia & Sill, B. L. & Reinhold, T. A., “Effects of surface roughness element spacing on boundary-layer velocity profile parameters”, Journal of Wind Engineering and Industrial Aerodynamics., 73 (1998) 215-230.
[19]Reiji Kimura & Kyoichi Otsuki & Makio Kamichika, “Relationships between the zero-plance displacement and the roughness length over sorghum and alfalfa canopies”, J. Agric. Meteorol., 55 (1999) 15-24.
[20]Macdonald, R. W. & Griffiths, R. F. & Hall, D. J., “An improved method for the estimation of surface roughness of obstacle arrays”, Atmospheric Environment, 32 (1998) 1857-1864.
[21]Bottema, M., “Roughness parameters over regular rough surfaces: Experimental requirements and model validation”, Journal of Wind Engineering and Industrial Aerodynamics, 64 (1996) 249-265.
[22]Molion, L. C. B. & Moore, C. J., “Estimating the zero-plane displacement for tall vegetation using a mass conservation method”, Boundary-Layer Meteorology, 26 (1983) 115-125.
[23]Tieleman, H. W., “Roughness estimation for wind-load simulation experiments”, Journal of Wind Engineering and Industrial Aerodynamics, 91 (2003) 1163-1173.
[24]Yan MENG & Masahiro MATSUI, “Characteristics of the vertical wind profile in neutrally atmospheric boundary layers”, Wind Engineers, JAWE, 65 (1995) 1-15.
[25]內政部建研所，建築物耐風設計規範及解說，建築出版社，2007。
[26]Holmes, John D., “Wind Loading of Structures”, 1st edit., Spon Press, 2001.
[27]Davenport, A. G., “The Spectrum of Horizontal Gustiness Near the Ground in High Winds”, Journal of Royal Meteorological. Society., 87 (1961) 194-211.
[28]Kaimal, J. C. et al., “Spectral Characteristic of Surface layer Turbulence ”, Journal of Royal Meteorological. Society., 98 (1972) 563-589.
[29]von Karman, T. “Progress in the Statistical Theory of Turbulence”,.Proceedings of the National Academy of Science of the United States of America., 34 (1948) 530-539.
[30]Shiau, B. S., “Velocity Spectra and Turbulence Statistics at the Northeastern Coast of Taiwan under High-wind condition”, Journal of Wind Engineering and Industrial Aerodynamics, 88 (2000) 139-151.
[31]Irwin, H.P.A.H., “The design of spires for wind simulation”, Journal of Wind Engineering and Industrial Aerodynamics., 7 (1981) 361-366.
[32]Browning, K. A. & Wexler, R. “The determination of Kinematic Properties of a wind field using Doppler radar”, Journal of applied meteorology, 7 (1968) 105-113.
[33]曾鈺婷，大氣邊界層特性之實場監測及風洞實驗模擬，淡江大學土木工程學系碩士學位論文，2011。
[34]郭庭瑋，不同地形之邊界層特性之風洞模擬，淡江大學土木工程學系碩士班畢業論文，2011。
[35]胡琮普，地表粗糙元素對於模擬大氣邊界層之影響，淡江大學土木工程學系碩士班畢業論文，2012。
[36]內政部建研所委託研究報告，標準地況自然風場實場監測與特性分析(三)，2011。