孫淑芳，張廣興

（1.烏魯木齊市氣象局，新疆 烏魯木齊830002；2.中國(guó)氣象局烏魯木齊沙漠氣象研究所，新疆 烏魯木齊830002）

20 世紀(jì)50 年代以來(lái)，由于航空的發(fā)展、高層建筑的修建和城市大氣污染的發(fā)生，風(fēng)隨高度的廓線直接影響近地面層生活的人們，地學(xué)界開(kāi)始重視大氣邊界層的研究。最具劃時(shí)代意義的工作是Monin和Obukhov[1]1954 年提出了Monin-Obukhov 相似性理論，建立了近地層湍流統(tǒng)計(jì)量和平均量之間的聯(lián)系。Wyngaard[2]1971 年提出了局地自由對(duì)流近似，補(bǔ)充了近地面層相似理論。Dyer[3]等1982 年利用1976年在澳大利亞國(guó)際湍流實(shí)驗(yàn)（ITCE）的數(shù)據(jù)對(duì)近地面層相似理論進(jìn)行了完善，使得該理論具有了極好的應(yīng)用價(jià)值。Neuwstadt[4]1984 年研究了穩(wěn)定的夜間邊界層的湍流結(jié)構(gòu)。Shao[5]1990 年研究了水平均勻邊界層的局部相似關(guān)系，進(jìn)一步探究了全邊界層的相似性理論。我國(guó)的胡隱樵等[6]開(kāi)展了野外實(shí)驗(yàn)，驗(yàn)證了局地相似性理論，并建立了各種局地相似性理論之間的關(guān)系。張強(qiáng)[7]等對(duì)局地相似性理論在非均勻下墊面近地面層的適應(yīng)性做了一些研究。隨著研究的深入相似性理論也成了大氣邊界層氣象學(xué)中最主要的分析和研究手段之一。而相似性理論最重要的應(yīng)用就是近地面層平均風(fēng)速廓線分布，由此導(dǎo)出的風(fēng)切變指數(shù)在高層建筑物的荷載，風(fēng)電渦輪機(jī)的結(jié)構(gòu)和輸電線路的選址、架設(shè)運(yùn)維方面得到了廣泛應(yīng)用，也是本文撰寫的理論基礎(chǔ)。有了風(fēng)切變指數(shù)，就可把一般氣象站標(biāo)準(zhǔn)10 m 高度測(cè)風(fēng)推算出高架設(shè)施高度的風(fēng)速；另一方面目前的WRF 模式采用σ坐標(biāo)，模式底部一般為30～60 m，就需要利用邊界層Monin-Obukhov 相似性理論的方法插值到10 m 風(fēng)向桿標(biāo)準(zhǔn)高度，風(fēng)切變指數(shù)就是風(fēng)廓線的風(fēng)速高度關(guān)系式，目前，區(qū)域數(shù)值模式對(duì)大風(fēng)區(qū)的10 m 風(fēng)速預(yù)報(bào)偏小[9]，部分原因是模式中適合大部分區(qū)域的風(fēng)廓線插值方案不能代表大風(fēng)區(qū)的風(fēng)廓線。

由此可見(jiàn)，α 是一個(gè)具有實(shí)際預(yù)報(bào)和工程建設(shè)應(yīng)用價(jià)值的重要參數(shù)，因此，往往在重要的區(qū)域建設(shè)鐵塔來(lái)觀測(cè)近地層風(fēng)廓線并計(jì)算風(fēng)切變指數(shù)，這方面各地做了不少工作。王志春等[8]計(jì)算了廣東省湛江市位于沿海徐聞?wù)镜娘L(fēng)切變指數(shù)。杜燕軍和馮長(zhǎng)青[10]研究了內(nèi)蒙古地區(qū)不同地形地貌以及地面粗糙度條件下的風(fēng)切變指數(shù)；彭懷午等[11]分析了內(nèi)蒙古地區(qū)風(fēng)切變指數(shù)的日和月時(shí)間變化；賀志明等[12]指出鄱陽(yáng)湖區(qū)風(fēng)切變指數(shù)較平原地區(qū)明顯偏小。楊興華等開(kāi)展的近地層、邊界層研究為本文提供了很好的參考[13-16]。金莉莉等[17]分析了烏魯木齊市近地層每座塔各層間風(fēng)切變指數(shù)以及主要特征，為本文提供了地理位置最近的參考范例。

中亞是寒潮冷空氣的關(guān)鍵區(qū)，冷空氣東移南下受天山阻攔易形成天山準(zhǔn)靜止鋒[18]，但冷空氣在天山山脈的博格達(dá)山與巴里坤山之間的七角井埡口可以翻山，在埡口南面的蘭新鐵路紅旗坎站至了墩站間123 km 寬度的背風(fēng)坡坡腳形成大風(fēng)區(qū)域，即通常所稱的“百里風(fēng)區(qū)”，地表為風(fēng)蝕戈壁地貌，最大風(fēng)速可達(dá)64 m·s-1，是12 級(jí)風(fēng)底限值的近兩倍[19]。該風(fēng)區(qū)具有風(fēng)速高、風(fēng)期長(zhǎng)、季節(jié)性強(qiáng)、風(fēng)向穩(wěn)定、起風(fēng)速度快等規(guī)律[20]，也有大風(fēng)爆發(fā)迅速，低空風(fēng)速?gòu)?qiáng)勁的特點(diǎn)[21]。在百里風(fēng)區(qū)這樣的大風(fēng)區(qū)研究風(fēng)切變指數(shù)無(wú)論是預(yù)報(bào)理論還是工程實(shí)用都有意義，但目前研究尚少。

1 資料與方法

1.1 資料

百里風(fēng)區(qū)十三間房的測(cè)風(fēng)塔上測(cè)風(fēng)儀器分別安裝在10、30、50、70、100 m 共5 層。2009 年5 月（建塔并開(kāi)始正常觀測(cè)）至2011 年1 月（停止觀測(cè)）20 個(gè)月逐時(shí)風(fēng)速資料，每層共14 000 個(gè)樣本，在統(tǒng)計(jì)上屬于大樣本，統(tǒng)計(jì)量具有代表性。取逐時(shí)10 min 平均風(fēng)速資料，對(duì)測(cè)風(fēng)塔原始數(shù)據(jù)，首先進(jìn)行各層時(shí)間序列上一致性判斷，再進(jìn)行各層間的一致性判斷。然后再用十三間房氣象站（43°13′N，91°44′E；海拔高度721.4 m）經(jīng)過(guò)中國(guó)氣象局審核的資料進(jìn)行極值比對(duì)和合理性判斷，以此對(duì)所用資料進(jìn)行了質(zhì)量控制。10 min 平均風(fēng)速以時(shí)為統(tǒng)計(jì)單位，給出各層平均風(fēng)速，避免反復(fù)平均導(dǎo)致精度降低。

1.2 計(jì)算方法

在近地層中，風(fēng)速隨高度變化顯著。造成風(fēng)在近地層中垂直變化的原因有動(dòng)力因素和熱力因素，前者主要來(lái)源于地面的摩擦效應(yīng)，即地面的粗糙度，后者主要是近地層大氣垂直穩(wěn)定度所致。當(dāng)大氣層結(jié)為中性時(shí)，湍流將完全依靠動(dòng)力原因來(lái)發(fā)展，這時(shí)風(fēng)速隨高度變化服從普朗特[22]經(jīng)驗(yàn)公式：

式中，k 為馮卡門（von Karman）常數(shù)，k=0.4；u*為摩擦速度，τ 為表面剪切應(yīng)力，ρ 為空氣密度。

圖1 測(cè)風(fēng)塔位置和測(cè)場(chǎng)地形地貌

近地層風(fēng)速的垂直分布主要取決于地表粗糙度和低層大氣的層結(jié)狀態(tài)。在中性大氣層結(jié)下，對(duì)數(shù)和冪指數(shù)方程都可以較好地描述風(fēng)速的垂直廓線，實(shí)測(cè)數(shù)據(jù)檢驗(yàn)結(jié)果表明，冪指數(shù)公式比對(duì)數(shù)公式能更精確地?cái)M合風(fēng)速的垂直廓線，我國(guó)新修訂的《建筑結(jié)構(gòu)設(shè)計(jì)規(guī)范》[23]和《風(fēng)電場(chǎng)風(fēng)能資源評(píng)估方法》[24]（GB／T18710—2002）均推薦使用冪指數(shù)公式，其表達(dá)式為：

由公式（4）變形得到：

公式（6）用來(lái)計(jì)算各層間的風(fēng)切變指數(shù)α。

2 結(jié)果與分析

2.1 整層風(fēng)切變指數(shù)

圖2 十三間房（#31005）測(cè)風(fēng)塔風(fēng)廓線圖

風(fēng)切變指數(shù)不僅受周圍地形粗糙度的影響，更多的是受大氣層穩(wěn)定度影響。大氣層穩(wěn)定度隨著季節(jié)和天氣過(guò)程相關(guān)的氣象條件而變化。相較于其他地區(qū)，本文研究區(qū)域，晴天日數(shù)多，太陽(yáng)輻射是最主要的因子。一日之內(nèi)，隨著日出，太陽(yáng)高度角的增大，大氣層穩(wěn)定度由夜間的穩(wěn)定變?yōu)橹行院筒环€(wěn)定，湍流增強(qiáng)，上下層大氣混合，上下層風(fēng)速大小差別減少，根據(jù)公式（3）風(fēng)切變指數(shù)α 則變小。圖3 顯示一天內(nèi)風(fēng)速≥3 m·的風(fēng)切變指數(shù)和太陽(yáng)輻射強(qiáng)度在時(shí)間上同步，這與沿海地區(qū)的結(jié)論較為一致[27，28]，也與金莉莉在烏魯木齊的研究結(jié)論一致[17]。8：00 左右日出后，隨著太陽(yáng)高度角增大α 緩慢變小，14：00 α 最小，一直維持最小值到20：00 太陽(yáng)落山，可以看到α 迅速增大，這和戈壁的地理環(huán)境相關(guān)。而大風(fēng)期間的α 值沒(méi)有明顯的日變化，因?yàn)?，該地是大風(fēng)區(qū)，大風(fēng)的發(fā)生頻率較高，大風(fēng)及以上風(fēng)速湍流很強(qiáng)，其強(qiáng)度遠(yuǎn)超太陽(yáng)輻射引起的湍流日變化，進(jìn)而導(dǎo)致α 值日變化不顯著。

圖3 風(fēng)切變指數(shù)日變化

圖4 給出了兩個(gè)風(fēng)速段的風(fēng)切變指數(shù)α 各月的值，代表了季節(jié)變化。可以看出，對(duì)于啟動(dòng)風(fēng)速段，3—10 月α 維持較小值，究其原因，夏季氣溫高、日照時(shí)間長(zhǎng)、湍流強(qiáng)、α 值小，物理意義與上一節(jié)分析的白天α 值小的原因一致。春秋，該地地處風(fēng)區(qū)風(fēng)口，寒潮強(qiáng)冷空氣活動(dòng)頻繁，風(fēng)力強(qiáng)勁，使得這兩季風(fēng)切變指數(shù)α 值較小。1—2 月和11—12 月是該地的冬半年，天氣寒冷，太陽(yáng)高度角變小，日照時(shí)間變短，湍流變?nèi)?，同時(shí)該地為吐鄯托盆地，冬季盆地特有的冷湖效應(yīng)使得邊界層內(nèi)100 m 塔高這一層大氣較為穩(wěn)定，因此α 值較大。在一年內(nèi)，該地風(fēng)切變指數(shù)α 最低的月份不僅和太陽(yáng)輻射最高的月份吻合，也與當(dāng)?shù)氐睦淇諝饣顒?dòng)相關(guān)，可見(jiàn)大風(fēng)區(qū)的風(fēng)切變指數(shù)α 要比一般沿海地區(qū)復(fù)雜。本文計(jì)算的風(fēng)切變指數(shù)α 的月和日變化規(guī)律與郭曉寧和保廣裕[29]在青海省格爾木站分析的太陽(yáng)輻射特征有較好的反向相關(guān)性，說(shuō)明本文的推論是正確的。

大風(fēng)風(fēng)速段風(fēng)切變指數(shù)α 值各月差異較小，但也可以看出冬季較大，其他三季較小，最小在8 月，其原因與上一節(jié)分析的一致，因?yàn)榘倮镲L(fēng)區(qū)的大風(fēng)天氣過(guò)程是高頻率事件且風(fēng)速很大，湍流很強(qiáng)，打破了太陽(yáng)輻射季節(jié)的影響，使得風(fēng)切變指數(shù)α 季節(jié)特征不顯著。

2.2 層間風(fēng)切變指數(shù)分布特征

圖4 風(fēng)切變指數(shù)季節(jié)變化

實(shí)際的風(fēng)電場(chǎng)與輸電線路以及其他如高鐵經(jīng)過(guò)風(fēng)區(qū)的設(shè)施建設(shè)與運(yùn)維，不僅要獲得總的風(fēng)切變指數(shù)α 值，以方便用常規(guī)的10 m 觀測(cè)的風(fēng)速計(jì)算出其他高度的風(fēng)速。更需要知道各層間的α 值，以便了解風(fēng)切變強(qiáng)度，避開(kāi)風(fēng)速?gòu)?qiáng)切變層。由表1 可見(jiàn)，啟動(dòng)風(fēng)速段的各層間的α 值，最小出現(xiàn)在100～70 m間，最大出現(xiàn)在50～30 m；大風(fēng)風(fēng)速段，最小出現(xiàn)在70～50 m，無(wú)獨(dú)有偶，最大也出現(xiàn)在50～30 m。這不是巧合，說(shuō)明十三間房站，50～30 m 有一個(gè)風(fēng)速切變，也就是這兩層風(fēng)速差異較其他各層要大。

表1 十三間房測(cè)風(fēng)塔各層間風(fēng)切變指數(shù)

夜間和白天的各層間風(fēng)切變指數(shù)α 值平均后繪制圖5。由圖5a 可見(jiàn)白天各層間的α 值沿著高度逐步減小，說(shuō)明白天由于太陽(yáng)輻射使得大氣層趨于中性和不穩(wěn)定，湍流加強(qiáng)，上下層混合使得層間風(fēng)速差異減少，相對(duì)應(yīng)α 值各層較為均一，風(fēng)切變不顯著。而夜間，在50～30 m 間α 值出現(xiàn)最大值，表明在這一層有一個(gè)風(fēng)速切變。

由圖5b 可見(jiàn)，冬季的風(fēng)切變指數(shù)α 極值位于70～50 m，主要是因?yàn)榕璧氐亩灸鏈貙虞^厚，使得70～50 m 和100～70 m 兩個(gè)層間α 值均較大。而夏季的極值在50～30 m。與圖5a 白天的極值出現(xiàn)原因相近，不再贅述。需要注意的是50～30 m 這一層，無(wú)論是日變化還是季節(jié)變化均存在一個(gè)拐點(diǎn)。

3 結(jié)論

利用吐鄯托盆地百里風(fēng)區(qū)的十三間房測(cè)風(fēng)塔僅有的一年半觀測(cè)資料，分析了風(fēng)區(qū)山盆地形和戈壁地貌條件下，啟動(dòng)風(fēng)速段和大風(fēng)風(fēng)速段，整層和層間風(fēng)切變指數(shù)α 的分布規(guī)律以及日變化和季節(jié)變化特征，得到以下結(jié)論：

（1）啟動(dòng)風(fēng)速段整層風(fēng)切變指數(shù)α 為0.084 5，大風(fēng)風(fēng)速段整層風(fēng)切變指數(shù)α 為0.036 4。2 個(gè)風(fēng)速段的各個(gè)層間風(fēng)切變指數(shù)α 最大值均出現(xiàn)在50～30 m，分別是0.109 6 和0.056 9。風(fēng)切變指數(shù)α 與風(fēng)速成反比關(guān)系，即風(fēng)速較小時(shí)風(fēng)切變指數(shù)α 較大，風(fēng)速較大時(shí)風(fēng)切變指數(shù)α 較小，較好體現(xiàn)了風(fēng)速越大垂直湍流交換越強(qiáng)，上下風(fēng)速更趨一致，進(jìn)而風(fēng)切變指數(shù)α 越小。

圖5 層間風(fēng)切變指數(shù)變化

（2）啟動(dòng)風(fēng)速段的風(fēng)切變指數(shù)α 白天較小，夜間較大。隨著日出風(fēng)切變指數(shù)α 緩慢變小，直至當(dāng)?shù)貢r(shí)間正午，達(dá)到最小值。日落后α 迅速變大，在日落2 h 達(dá)到夜間的較大值。一日內(nèi)對(duì)應(yīng)于日出日落，α 日變化曲線接近一個(gè)正弦波形。各層間風(fēng)切變指數(shù)α 白天較為均一，夜晚在50～30 m 有一個(gè)極大值拐點(diǎn)，對(duì)應(yīng)夜晚在該層有一個(gè)風(fēng)速切變。大風(fēng)風(fēng)速段幾乎沒(méi)有日變化。

（3）風(fēng)切變指數(shù)α 無(wú)論啟動(dòng)風(fēng)速段還是大風(fēng)風(fēng)速段均呈現(xiàn)出冬季大，夏季小的季節(jié)特征，但啟動(dòng)風(fēng)速段季節(jié)特征更顯著。各層間的風(fēng)切變指數(shù)α 值，冬季的最大值在70～50 m，而夏季的極值在50～30 m，且存在顯著的拐點(diǎn)。