許超宗，陳蜀江*，黃鐵成，朱 選，姚艷霞，李紅軍，郭 超

(1.新疆師范大學地理科學與旅游學院，新疆 烏魯木齊 830054；2.北京林業(yè)大學，北京 100083；3.烏魯木齊空間遙感應用研究所，新疆 烏魯木齊 830054；4.莫納什大學，澳大利亞 墨爾本 3000；5.中國林業(yè)科學研究院森林生態(tài)環(huán)境與保護研究所，北京 100083；6.西天山國家級自然保護區(qū)，新疆 伊犁 835000)

【研究意義】隨著全球氣候變化、能源危機的出現(xiàn)，可再生能源漸漸成為世界各國關注的重點[1]。風能作為可再生清潔能源，具有安全、清潔、蘊藏量大、開發(fā)技術成熟等優(yōu)點，對社會能源結(jié)構的調(diào)整和環(huán)境保護方面具有重要的意義。風速是重要氣候要素之一，它的變化對氣候變化及氣候形成具有一定的反應[2-3]，被廣泛應用于氣候變化研究、風蝕強度評估、沙塵暴預測、大氣污染評價等領域[4-5]。【前人研究進展】近些年，有關風速變化的研究逐漸增加。研究表明：從全球范圍來看，過去30～50年間中緯度大部分地區(qū)近地面風速變化呈現(xiàn)減小趨勢，風速減小的地區(qū)包括:澳大利亞、中國、歐洲、北美洲[6]。而高緯度地區(qū)(緯度>65°)南極洲和阿拉斯加風速增加大約0.005 m/s[7]。與此同時，中國大部分地區(qū)風速有減小的趨勢:我國平均風速以-0.012 m/(s·a)[8]的速率在遞減，王遵婭等[9]對中國近50年氣候變化研究表明，中國大部分地區(qū)風速均顯著減小，其中冬、春季較為明顯。WANG et al[10]研究中國1951-2000年風速變化特征顯示，我國普遍存在風速降低的現(xiàn)象，以西北地區(qū)下降的幅度最大，季節(jié)最大降幅發(fā)生在冬季。何毅等對南北疆風速變化研究表明，北疆地區(qū)平均風速呈現(xiàn)明顯下降趨勢，遞減率為-0.014 m/(s·a)[11]，20世紀90年代和50年代相比，全國平均風速減小了約16 %，減小最顯著的西北西部，減小了近29 %，幅度非常大[9]。因而對于西北地區(qū)尤其是西北西部風速變化的研究就尤極為重要，風是影響沙塵天氣最為直接、最優(yōu)相關的因子，風速的大小直接關系到沙塵天氣頻數(shù)及強度的變化[12]，風速的變化也將直接影響水循環(huán)，風能發(fā)展?jié)摿εc風速變化也休戚相關?！颈狙芯壳腥朦c】伊犁地區(qū)位于新疆維吾爾自治區(qū)西部，42°14′～44°53′N，80°09′～84°56′E，地處南北天山(高度為3000 m以上)之間，位于中緯度內(nèi)陸天山山區(qū)，氣候為大陸性溫帶半干旱氣候。地形東、南、北三面高山環(huán)繞，地勢東高西低，東窄西寬，呈喇叭型向西敞開。伊犁地區(qū)地處西風環(huán)流區(qū)，在我國西風帶上游，同時受到副熱帶高壓的影響，春、夏之交盛行西風和西北風，是我國著名的多風地區(qū)[13]。伊犁地區(qū)作為整個新疆重要的糧食基地，素有“新疆糧倉”之稱。而大風會產(chǎn)生沙塵暴災害，造成土壤風蝕、沙化，生態(tài)環(huán)境遭受破壞，危害農(nóng)業(yè)生產(chǎn)[14-15]?！緮M解決的關鍵問題】本文主要對伊犁地區(qū)年平均風速變化特征進行研究，既可揭示中高緯度西風帶控制下干旱區(qū)平均風速的變化規(guī)律，也可為當?shù)厣鷳B(tài)環(huán)境的改善、能源合理開發(fā)利用、工農(nóng)業(yè)生產(chǎn)提供科學依據(jù)。

1 材料與方法

1.1 資料來源

依據(jù)氣象數(shù)據(jù)資料的連續(xù)性及最長時段性等標準，本研究對缺測和誤測的數(shù)據(jù)進行剔除及預測(線性回歸方法預測)，經(jīng)過嚴格的質(zhì)量控制，選取質(zhì)量較好的氣象站8個(表1)。地面氣象站的平均風速實測日值資料來自中國氣象科學數(shù)據(jù)共享服務網(wǎng)站國家氣象地面基準站。為確保所選氣象站數(shù)據(jù)資料長度的均一性與穩(wěn)定性，氣象資料時間跨度均為1960年1月1日至2016年12月31日。按照氣候方法劃分四季，即3-5月為春季，6-8月為夏季，9-11月為秋季，12月至次年2月為冬季。

表1 伊犁地區(qū)氣象站點位置

1.2 研究方法

利用Excel對實測的風速數(shù)據(jù)進行處理，研究方法主要選取氣候傾向率法、6階多項式擬合曲線等數(shù)理統(tǒng)計方法；對研究區(qū)近56年風速變化進行分析，用M-K突變檢驗法、小波分析等方法對風速變化的突變性和周期性進行分析；對空間上的變化主要采用反距離空間插值進行分析。

(1)氣候傾向率法。氣候傾向率反映氣候要素的變化趨勢，一般可采用一元線性回歸方程的方法計算：

y=b+at，t=1961，1969，…，2016(t—年份序號)

(1)

式中，b為截距；a為回歸系數(shù)也稱為傾向率。當回歸系數(shù)為正值時，即a>0,表明要素y在計算時間內(nèi)是線性增加，反之則是減小的。

(2)M-K突變檢驗。Mann-Kendall突變檢驗簡稱M-K檢驗，是非參數(shù)的統(tǒng)計檢驗法，其不受異常值干擾，也無需遵循樣本分布。在統(tǒng)計中，UF大于0說明序列呈上升趨勢，小于0則相反，超出臨界線時，表明變化趨勢顯著。若UF和UB在臨界線內(nèi)交于一點，則其對應的年份即為突變的時間點。

(3)小波分析。小波分析(Wavelet Analysis)是一種具有時-頻多分辨率功能的分析方法，能清晰揭示出隱藏在時間序列中的多種變化周期，充分反映系統(tǒng)在不同時間尺度的變化趨勢，是研究氣象要素長期變化的重要工具[16-18]。

圖1 1961-2016年伊犁地區(qū)年均風速變化趨勢Fig.1 Trend of annual average wind speed in Yili area from 1961 to 2016

2 結(jié)果與分析

2.1 時間變化

2.1.1 年際變化特征及趨勢分析 近56年間，伊犁地區(qū)風速總體處于波動下降趨勢，傾向率為-0.171 m/(s·a)，與全國風速變化吻合均處于遞減狀態(tài)，遞減速率高于全國速率[0.012 m/(s·a)][19]。其中年平均風速17.76 m/s，最高風速1971年23.59 m/s，最低風速12.95 m/s(1998年)，二者相差10.64 m/s。6階多項式擬合曲線顯示，伊犁地區(qū)平均風速變化具有明顯的階段性特征，大致可以分為2個時期和5個階段：1961-1985年的相對大風時期和1986-2016年相對小風時期；1961-1964年、1974-1995年、2010-2016年的下降階段和1965-1973年、1996-2010年的上升階段，這與劉櫟杉[20]等研究的全疆風速變化具有一致性。年均風速的變化趨勢為波動下降—緩慢增長—趨于穩(wěn)定。

圖2 1961-2016年伊犁地區(qū)平均風速M-K突變檢驗Fig.2 M-K mutation test of average wind speed in Yili area from 1961 to 2016

2.1.2 突變檢驗 研究區(qū)平均風速M-K突變檢驗(圖2)表明，年均風速UF曲線多位于0值以下，風速呈波動下降趨勢。1972年UF﹥0，風速呈增加趨勢，其他年份UF均位于小于0，風速減小，并在1980年以后超過了0.01顯著性水平臨界值置信區(qū)間(U0.01=2.58)，風速達到極顯著減小。UF與UB曲線在置信區(qū)間內(nèi)有一交點1979年，經(jīng)過滑動t檢驗，此點通過了檢驗，為伊犁地區(qū)年均風速的突變時間。

2.1.3 季節(jié)變化特征 四季風速變化與全年風速變化基本一致，均表現(xiàn)為波動下降趨勢(圖3)。四季平均風速依次為22.07、18.47、16.40和14.08 m/s，風速最大為春季，冬季風速最小。傾向率變化顯示，下降最顯著的為春季-0.196 m/(s·a)，其次為秋季-0.185 m/(s·a)、夏季-0.166 m/(s·a)，冬季下降最小-0.139 m/(s·a)。6階多項式擬合曲線表明，春季、秋季和冬季平均風速變化階段具有一致性，大致可分為4個階段：2個風速增加時期和2個風速減小時期，各季節(jié)風速增加、減小的階段性又存在差異；夏季平均風速的變化較為特殊，可分為2個增加期和3個減小期。

圖3 1961-2015年伊犁四季平均風速變化趨勢Fig.3 Trends in average wind speed over the four seasons of Yili from 1961 to 2015

平均風速1960s1970s1980s1990s2001-2010春距平4.314.00-0.45-3.59-2.33c.v.0.0870.0810.1260.0590.026夏距平4.113.3-0.62-3.51-2.21c.v.0.0590.0480.1260.0470.037秋距平4.353.070.58-3.96-2.67c.v.0.0350.0720.0940.0730.041冬距平3.582.42-0.24-2.92-2.04c.v.0.1060.1210.0930.0960.069年距平4.263.24-0.12-3.52-2.38c.v.0.0420.0640.1020.0530.028

距平值和變異系數(shù)(c.v.)分別反映氣候要素相對于某個長期平均值的高低程度和離散情況[21]。從風速年代際變化情況(表2)分析，全年風速距平在1960s和1970s為正值(偏大)，1980s后為負值風速轉(zhuǎn)小，變異系數(shù)最大值出現(xiàn)在80s，最小值出現(xiàn)在2001-2010年，表明風速在1980s時期波動的幅度最為顯著。四季風速年代際變化中，春季、秋季和冬季的變化基本一致，表現(xiàn)為1960s和1970s期間風速偏大，20世紀80年代至今，風速偏小，與研究區(qū)全年風速的變化保持一致；夏季風速變化較為復雜，3個風速偏小時期和2個風速偏大期。變異系數(shù)則表現(xiàn)為：秋季變化幅度最小，其次為夏季、春季，冬季變化幅度最大，表明秋季風速較為穩(wěn)定波動不大，冬季風速離散程度最大，變化較為劇烈。春季、夏季和秋季在1980s變化較劇烈，冬季則在1970s變化幅度較劇烈。

2.1.4 周期性特征 小波系數(shù)實部等值線圖表明，研究區(qū)年平均風速的周期性變化存在2個周期震蕩(圖4)，即18和28年。18年震蕩周期貫穿于整個研究時間段，在全年平均風速的震蕩周期最為強烈，為時序變化的主要周期。28年周期在1975-2005年期間較為明顯，其余時間內(nèi)較為弱化。對四季風速周期性變化而言，春季和夏季的變化周期與全年風速變化周期吻合，均存在18、28年2種年尺度周期變化。較為特殊的是夏季在1970-2000年間存在著一個較為完整的6年次周期，且夏季的18年周期震蕩較春季強。秋季主要存在18、22、8～10年2種年尺度周期變化，其中22年年尺度周期震蕩在2010年后震蕩逐漸減弱。冬季主要存在10、18、22年3種年尺度周期變化，10年震蕩周期貫穿于整個時間段，且由強逐漸減弱；1973年以后22年震蕩周期開始出現(xiàn)。四季風速的周期性變化與全年風速的變化具有統(tǒng)一性，春、夏、秋、冬四季與全年在18年尺度的震蕩周期上存在一致性，但四季具有各自的特殊性。

2.2 空間變化

利用Arcgia反距離空間插值法結(jié)合研究區(qū)8個氣象站點對其風速在空間上的變化進行分析(圖5)表明，年平均風速在空間上變化主要由西北向東、向西遞減，具有兩個明顯的高值中心和2個低值中心，高值中心分別為尼勒克(5.78 m/s)、察布查爾(5.08 m/s)，低值中心霍城(1.70 m/s)、新源(1.83 m/s)。整個伊犁地區(qū)風速傾向率的變化均表現(xiàn)為減小的趨勢，風速傾向率空間變化與平均風速變化大體呈反方向變化，即尼勒克、察布查爾等年平均風速較高的地區(qū)，傾向率減小的速率較慢，反之平均風速較低的地區(qū)，其減小的速率較快。風速的這種空間變化與其所處的地理位置(中緯西風帶)、地形等條件具有密切的聯(lián)系。

3 討 論

伊犁地區(qū)風速的變化趨勢與全國基本吻合，但其減小的速率要高于全國的速率，這一變化與前人[11,13]對新疆風速變化的研究一致。空氣的水平運動成為風，它是氣團之間存在溫差、出現(xiàn)氣壓梯度而產(chǎn)生的氣體流動，因此風速的減小與氣候變暖也有一定的關系[22]，XUE et al[23]的研究認為中國北方地區(qū)比南海地區(qū)變暖的幅度大，主要與我國冬季風的減弱有關。首先，隨著全球氣候變暖的發(fā)生，我國新疆地區(qū)氣候的變化朝著暖濕化發(fā)展[24]，伊犁地區(qū)也不例外，因而，其風速減小與當?shù)貧夂蜃兣嬖谝欢ǖ年P聯(lián)性。特別是我國和亞歐大陸近年來的氣溫上升，使得冬季較高的地面氣溫可能降低地表氣壓，進而減弱陸地和鄰近海洋之間的氣溫和氣壓梯度，最終降低氣壓梯度力并導致風速下降[25]。其次，伊犁地區(qū)屬于中緯西風帶，且地形呈喇叭型向西敞開，易受西風環(huán)流的影響，其風速的大小會受到中緯西風強弱的影響。最后，隨著我國城市化速度及水平的加快，迅速崛起的建筑物對風速有一定的阻擋作用，會使平均風速減小。江瀅等[26]研究表明，全國城市氣象站風速減小趨勢要比鄉(xiāng)鎮(zhèn)氣象站更為明顯，說明城市化擴大影響了風速的變化。國家對新疆發(fā)展的支持與援助，使得新疆地區(qū)的經(jīng)濟發(fā)展、社會建設等各方面取得了進步，伊犁地區(qū)城市化的發(fā)展對伊犁地區(qū)風速的變化也具有重要影響。

圖4 1961-2015年伊犁風速小波分析Fig.4 Yili wind speed wavelet analysis from 1961 to 2015

圖5 伊犁地區(qū)年平均風速空間分布Fig.5 Spatial distribution of annual average wind speed in Yili area

4 結(jié) 論

(1)1961-2016年伊犁地區(qū)平均風速呈明顯的下降趨勢，其遞減速率為-0.171 m/(s·a)，高于全國遞減速率。四季變化趨勢與全年基本吻合，春季遞減速率最高，冬季最小。

(2)研究時間段內(nèi)，該區(qū)年平均風速在1979年前后出現(xiàn)了突變。此研究與劉櫟彤[20]對全疆風速變化研究的突變年份1981年基本相似，但與全國的突變點1969年不一致，說明風速的突變存在地區(qū)的空間差異性。

(3)周期性研究表明，伊犁地區(qū)平均風速存在明顯的18年震蕩周期，且貫穿于整個研究時間段，全年和四季都符合此震蕩周期。但四季風速的變化周期存在差異性，夏、秋季節(jié)基本吻合全年風速震蕩周期變化，冬、春季節(jié)較為特殊，存在早于顯著震蕩周期(18年)的周期。

(4)空間上，伊犁地區(qū)平均風速的變化趨勢主要是由西北向東、向西遞減，存在2個高值中心和2個低值中心；傾向率的空間變化與年均風速空間分布呈反方向。