龔璽 徐金閣

摘要利用宜昌市1961—2020年逐日2 min平均風(fēng)速資料，采用氣候傾向法、MannKendall檢驗法和小波分析法，從不同時間尺度分析平均風(fēng)速的變化趨勢、突變和周期特征，揭示三峽局部地區(qū)長時間序列風(fēng)速變化規(guī)律。結(jié)果表明：近60年來宜昌市年平均風(fēng)速有增大的趨勢，遞增率為0.084 m/（s·10 a）；四季平均風(fēng)速均呈略增大趨勢，四季氣候傾向率相近。8月平均風(fēng)速增大最快，6月平均風(fēng)速增大最慢。MannKendall突變檢驗分析得出宜昌市年平均風(fēng)速在1971、1996和2013年發(fā)生突變。小波分析結(jié)果表明，年平均風(fēng)速變化有周期性規(guī)律，其第一主周期為32年。

關(guān)鍵詞風(fēng)速；變化特征；氣候傾向；MannKendall檢驗；小波分析；宜昌市

中圖分類號S162? 文獻(xiàn)標(biāo)識碼A? 文章編號05176611（2024）07020504

doi：10.3969/j.issn.05176611.2024.07.048

Study on the Change Characteristics of Wind Speed in Yichang City from 1961 to 2020

GONG Xi1，XU Jin-ge2

（1.National Meteorological Information Center， Beijing 100081;2.China Energy Construction Group Hunan Electric Power Design Institute Co.， Ltd.， Changsha， Hunan 410007）

AbstractUsing the daily 2-minute average wind speed data from 1961 to 2020 in Yichang City， the climate tendency method， Mann-Kendall test method and wavelet analysis method were used to analyze the change trend， mutation， and periodic characteristics of average wind speed at different time scales， revealing the long-term wind speed variation patterns in the local area of the Three Gorges area.The results showed that the annual mean wind speed in Yichang City had increased in recent 60 years with a decreasing rate of 0.084 m/（s·10 a）， and the four seasons mean wind speed had a slightly increasing trend with a similar climate tendency rate. The increase of average wind speed was the fastest in August and the slowest in June. The Mann-Kendall mutation test showed that the annual mean wind speed in Yichang City had a sudden change in 1971，1996 and 2013. Wavelet analysis showed that the annual mean wind speed changes periodically， and the first main period of the annual mean wind speed change was 32 years.

Key wordsWind speed;Change characteristic;Climate tendency;Mann-Kendall test;Wavelet analysis;Yichang City

長江三峽工程是世界上最大的水利樞紐之一，它建成后對國民經(jīng)濟(jì)產(chǎn)生巨大效益，同時也對局地的天氣氣候產(chǎn)生了影響。武慧鈴等［13］通過分析三峽水庫及其周邊氣象站的降水?dāng)?shù)據(jù)，發(fā)現(xiàn)三峽水庫蓄水改變了下墊面，水庫年均氣溫呈弱的下降趨勢，年降水量與降水日數(shù)增加，年均相對濕度減小，年平均風(fēng)速呈增加趨勢，蒸發(fā)量有所下降。封瑞雪等［45］利用重慶多年氣象站資料，分析了各氣象要素的年際變化，尤其是蓄水前后的差異情況。陳祥義等［6］利用三峽庫區(qū)及周邊氣象站的逐月降水?dāng)?shù)據(jù)，對三峽庫區(qū)1951─2012年年均、雨季、旱季降水量與降水日數(shù)的變化趨勢進(jìn)行了研究，并對三峽庫區(qū)內(nèi)年降水量和降水日數(shù)的空間分布情況進(jìn)行了分析。黃治勇等［7］研究了長江三峽庫區(qū)持續(xù)12 h以上和連續(xù)3 d以上極端大霧天氣氣候變化特征，探討了庫區(qū)蓄水后大霧天氣氣候變化的原因。張?zhí)煊畹龋?］利用1998—2016年TRMM和CMORPH衛(wèi)星降水資料和三峽庫區(qū)氣象觀測站數(shù)據(jù)，分析了三峽庫區(qū)局地降水變化。董新寧等［9］利用區(qū)域氣候模式RegCM3和NCEP（NNRP1）再分析資料模擬了2006年重慶三峽區(qū)域的氣候情況。

三峽庫區(qū)建成后除了提供豐富的水能資源以外，也存在風(fēng)能資源的應(yīng)用。前人對三峽庫區(qū)的降水、氣溫、相對濕度等要素進(jìn)行了比較詳盡的分析［513］，但是對風(fēng)的研究還較少。筆者利用宜昌市國家基本氣象站1961—2020年逐日2 min平均風(fēng)速資料，采用氣候傾向法、MannKendall檢驗法和小波分析法，分析平均風(fēng)速在不同月份的變化趨勢、突變和周期特征，揭示三峽地區(qū)長時間序列風(fēng)速變化規(guī)律，為三峽地區(qū)氣候變化提供參考。

1資料與方法

1.1數(shù)據(jù)來源

宜昌地處湖北省西南部，長江三峽起始地，三峽大壩、葛洲壩所在地。該研究使用的數(shù)據(jù)為宜昌市國家基本氣象站（30°42′N，111°18′E）1961—2020年逐日2 min平均風(fēng)速的觀測資料，對缺測數(shù)據(jù)進(jìn)行了補齊［14］。

1.2研究方法

1.2.1氣候傾向法。

假設(shè)y表示氣候要素，用x表示時間序列：

y=β0+β1x+ε（1）

式中：β0為回歸常數(shù)；β1為回歸系數(shù)。

1.2.2MannKendall檢驗。

對于n個樣本量的時間序列x，構(gòu)造一秩序列：

Sk=ki=1ri（k=2，3，…，n）（2）

其中，

ri=±1xi>xj0xi≤xj（j=1，2，…，i）（3）

定義統(tǒng)計量：

UFk=SkE（Sk）Var（Sk）（k=1，2，…，n）（4）

UF1=0時符合標(biāo)準(zhǔn)正態(tài)分布，E（Sk）和Var（Sk）是累計數(shù)Sk的均值和方差，在x1，x2，…，xn相互獨立且連續(xù)分布時，均值和方差可由下式計算得出：

E（Sk）=n（n+1）4（5）

Var（Sk）=n（n1）（2n+5）72（6）

給定顯著性水平α，查正態(tài)分布表可得Uα/2，若|UFi|>|Uα/2|，則在顯著水平α下序列具有顯著的趨勢變化，如α=0.05，置信區(qū)間臨界值uα/2=±1.96。將時間序列x逆序，重復(fù)上述計算過程，同時使

UBk=UFkk=n+1k（k=1，2，…，n）（7）

繪制UFk和UBk曲線圖，若UFk>0，則序列在第k時刻呈上升趨勢，反之，則呈下降趨勢；當(dāng)UFk超過臨界值，表明上升或下降趨勢顯著。若UFk和UBk 2條曲線出現(xiàn)交點，且交點位于置信區(qū)間內(nèi)，那么交點對應(yīng)時刻即為突變開始的時刻。其中變化趨勢的大小可用傾斜度β（也稱為趨勢系數(shù)）表示，計算式如下：

β=medianxixjijj

式中，n為樣本數(shù)量，若β>0，則說明樣本呈上升或增大趨勢，反之，樣本呈下降或減小趨勢。β的絕對值即為增減率。

1.2.3小波分析。

將任意L2（R）空間中的函數(shù)f（t）在小波基函數(shù)下展開：

Wf（a，b）≤f（t）（12）

Ψa，b（t）≥1|a|∫∞t=∞f（t）Ψtbadt（13）

Wf（a，b）為小波變換系數(shù)，與頻率參數(shù)a和時間參數(shù)b有關(guān)。將（13）式離散化為：

Wf（a，b）=1|a|Δtf（kΔt）ΨkΔtba（14）

將復(fù)值Morlet小波作為母函數(shù)，如下式：

Ψ（t）=eiωtet2/22eω2t2/4（15）

Morlet小波變換包含模、實部、虛部和幅值等變量，模的大小表示變量信息在不同時間尺度上的強弱，實部表示變量在不同時間尺度上的分布和位相兩方面信息。小波系數(shù)反映了變量在不同時間尺度下的變化幅度，其絕對值越大，說明在對應(yīng)的時間尺度下的變化越顯著，小波系數(shù)峰值對應(yīng)的時間尺度值即為變量變化的主要周期［15］。

2結(jié)果與分析

2.1趨勢分析

2.1.1平均風(fēng)速的年際變化趨勢。

從1961—2020年宜昌市年平均風(fēng)速變化趨勢（圖1）可以看出，近60年來宜昌市年平均風(fēng)速為1.3 m/s，與張靜等［3］對三峽庫區(qū)年平均風(fēng)速的研究成果一致。

1990年以前年平均風(fēng)速變化起伏明顯，20世紀(jì)60年代呈不明顯減小趨勢，70年代基本呈顯著增大趨勢，自1977年開始減小，持續(xù)至1989年。最大值出現(xiàn)在1976年，為1.8 m/s；1971年年平均風(fēng)速顯著增大。1990年以后年平均風(fēng)速基本波動中呈上升趨勢，20世紀(jì)90年代年平均風(fēng)速顯著增大，2000—2005年顯著減小，2006年略有增大后，呈緩慢減小至2013年，2014年顯著增大，后呈緩慢減小趨勢。

采用氣候傾向法對1961—2020年的年平均風(fēng)速進(jìn)行趨勢分析，最小二乘法得到回歸方程，由趨勢圖和回歸方程（圖1）可知，近60年來宜昌市年平均風(fēng)速總體呈略增大趨勢，氣候傾向率為0.084 m/（s·10 a）（相關(guān)系數(shù)為0.446，通過α=0.01的顯著性檢驗），尤其是自2006年開始年平均風(fēng)速的變化總體轉(zhuǎn)為上升趨勢，除大氣環(huán)流變化影響［16］，宜昌市年平均風(fēng)速變化與三峽水庫蓄水后水面粗糙度減小有關(guān)［3］。

2.1.2平均風(fēng)速的季節(jié)變化趨勢。

采用氣候傾向法對1961—2020年宜昌市各季節(jié)平均風(fēng)速的變化趨勢進(jìn)行回歸分析，計算得各季節(jié)氣候傾向率，結(jié)果見表1。由表1可知，四季平均風(fēng)速均呈略增大趨勢，四季氣候傾向率相近，春季平均風(fēng)速以0.089 m/（s·10 a）遞增最快，秋季以0.081 m/（s·10 a）遞增最慢。

2.1.3平均風(fēng)速的月變化趨勢。

分析1961—2020年宜昌市各月平均風(fēng)速的變化趨勢，結(jié)果發(fā)現(xiàn)（圖2），4月平均風(fēng)速最大，為1.5 m/s，10月平均風(fēng)速最小，為1.1 m/s。對各月平均風(fēng)速采用氣候傾向法進(jìn)行回歸分析，結(jié)果見表2。由表2可知，近60年宜昌市各月平均風(fēng)速均呈增大趨勢；其中，8月平均風(fēng)速增大最快，氣候傾向率為0.099 m/（s·10 a），6月平均風(fēng)速增大最慢，氣候傾向率為0.066 m/（s·10 a）。

2.1.4平均風(fēng)速的日變化趨勢。

從1961—2020年宜昌市近地面風(fēng)速日均值變化特征（圖3）可以看出，總體而言，宜昌市日均風(fēng)速先隨日序增加而增加，大約在第90天達(dá)到最大值，而后逐漸平穩(wěn)，之后又隨著日序增加而減小，在約第300天降到最小值，最后日均風(fēng)速基本穩(wěn)定在1.2 m/s左右。與1991—2020年氣候標(biāo)準(zhǔn)值相比，60年的日均值略小于近30年的氣候標(biāo)準(zhǔn)值，說明近30年宜昌市近地面風(fēng)速略有增大。

2.2突變檢驗分析

2.2.1年平均風(fēng)速。

采用MannKendall檢驗法對宜昌市1961—2020年年平均風(fēng)速序列進(jìn)行突變檢驗，結(jié)果見圖4。MannKendall檢驗法的UF和UB兩變量曲線交點位于1971、1996和2013年，且以上交點位于uα=±1.96的置信區(qū)間。

2.2.2季節(jié)平均風(fēng)速。

采用MannKendall檢驗法對1961—2020年宜昌市各季節(jié)平均風(fēng)速序列進(jìn)行突變檢驗，結(jié)果發(fā)現(xiàn)（圖5），MannKendall檢驗法秋季UF和UB兩變量曲線位于uα=±1.96的置信區(qū)間的交點有5個，其他3個季節(jié)均為2個；春季、秋季平均風(fēng)速的突變點均出現(xiàn)于1973和2013年，夏季出現(xiàn)于1971和2013年，冬季出現(xiàn)于1973和2011年，說明20世紀(jì)70年代夏季平均風(fēng)速的突變早于其他3個季節(jié)，21世紀(jì)冬季平均風(fēng)速突變早于其他3個季節(jié)。

2.3周期分析

通過小波分析法對1961—2020年宜昌市年平均風(fēng)速序列進(jìn)行周期性分析，由時頻分布（圖6a）可知，宜昌市年平均風(fēng)速在10～15、30～35年的時間尺度上存在周期性。

由小波方差（圖6b）可知，32年的時間尺度為年平均風(fēng)速變化的第一主周期。32年的時間尺度上，年平均風(fēng)速表現(xiàn)為顯著的“增—減—增”的變化規(guī)律，小波系數(shù)在1961—1982年為正值，1983—2004年為負(fù)值，2005—2020年為正值，且后一正值中心數(shù)值大于前一正值中心數(shù)值，表明年平均風(fēng)速在進(jìn)入21世紀(jì)20年間增大趨勢較20世紀(jì)60—70年代更顯著。

3結(jié)論

該研究利用宜昌市1961—2020年的風(fēng)速觀測資料，采用氣候傾向法、MannKendall檢驗法和小波分析法，分析了平均風(fēng)速的變化趨勢、突變和周期性特征，揭示了三峽地區(qū)風(fēng)速變化規(guī)律，得到的結(jié)論如下：

（1）近60年來宜昌市年平均風(fēng)速總體為增大趨勢，氣候傾向率為0.084 m/（s·10 a），相關(guān)系數(shù)為0.446，通過α=0.01的顯著性檢驗。

（2）四季平均風(fēng)速均呈略增大趨勢，四季氣候傾向率相近。8月平均風(fēng)速增大最快，6月平均風(fēng)速增大最慢。

（3）MannKendall突變檢驗分析得出宜昌市年平均風(fēng)速在1971、1996和2013年發(fā)生突變，春季、秋季平均風(fēng)速的突變點均出現(xiàn)于1973和2013年，夏季出現(xiàn)于1971和2013年，冬季出現(xiàn)于1973和2011年。

（4）小波分析結(jié)果表明，宜昌市年平均風(fēng)速在10～15、30～35年的時間尺度上存在周期性變化；32年的時間尺度為年平均風(fēng)速變化的第一主周期。

參考文獻(xiàn)

［1］? 武慧鈴，周建中，田夢琦，等.三峽水庫蓄水前后氣候變化分析［J］.水力發(fā)電，2021，47（5）：30-35.

［2］ 符坤，張六一，任強.蓄水前后三峽庫區(qū)氣候時空變化特征［J］.環(huán)境影響評價，2018，40（3）：82-86，96.

［3］ 張靜，劉增進(jìn)，肖偉華，等.三峽水庫蓄水后庫區(qū)氣候要素變化趨勢分析［J］.人民長江，2019，50（3）：113-116，165.

［4］ 封瑞雪，劉軍旗，姚夢輝，等.三峽水庫蓄水前后重慶氣候變化分析［J］.長江流域資源與環(huán)境，2019，28（4）：994-1002.

［5］ 周小英，謝世友，任偉.1955—2014年三峽庫區(qū)降水特征分析：以重慶市萬州區(qū)為例［J］.西南大學(xué)學(xué)報（自然科學(xué)版），2017，39（10）：102-108.

［6］ 陳祥義，肖文發(fā)，黃志霖，等.1951—2012 年三峽庫區(qū)降水時空變化研究［J］.生態(tài)環(huán)境學(xué)報，2015，24（8）：1310-1315.

［7］ 黃治勇，牛奔，葉麗梅，等.長江三峽庫區(qū)極端大霧天氣的氣候變化特征［J］.長江流域資源與環(huán)境，2012，21（5）：646-652.

［8］ 張?zhí)煊?，桂術(shù)，楊若文，等.TRMM和CMORPH衛(wèi)星資料對三峽庫區(qū)降水的評估分析［J］.氣象，2020，46（8）：1098-1112.

［9］ 董新寧，程炳巖，王勇，等.三峽區(qū)域氣候變化模擬研究［J］.西南大學(xué)學(xué)報（自然科學(xué)版），2012，34（1）：111-116.

［10］? 陳鮮艷，張強，葉殿秀，等.三峽庫區(qū)局地氣候變化［J］.長江流域資源與環(huán)境，2009，18（1）：47-51.

［11］ 陳鮮艷，張強．長江三峽局地氣候監(jiān)測（1961—2007年）［Ｍ］．北京：氣象出版社，2009．

［12］ 陳鮮艷，廖要明，張強．長江三峽工程生態(tài)與環(huán)境監(jiān)測系統(tǒng)：三峽氣候及影響因子研究［Ｍ］．北京：氣象出版社，2011．

［13］ 陳鮮艷，宋連春，郭占峰，等.長江三峽庫區(qū)和上游氣候變化特點及其影響［J］.長江流域資源與環(huán)境，2013，22（11）：1466-1471.

［14］ 徐金閣，龔璽，戴昌明.近60年洞庭湖南部地區(qū)風(fēng)速變化特征研究［J］.氣象與環(huán)境科學(xué)，2022，45（3）：83-88.

［15］ 魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù)［M］.2版.北京：氣象出版社，2007.

［16］ 劉燕.氣象站遷站前后主要氣象要素的差異性分析及對預(yù)報服務(wù)工作的影響［J］.氣象與環(huán)境科學(xué)，2017，40（2）：138-143.

作者簡介龔璽（1990—），女，江蘇泰州人，工程師，博士，從事風(fēng)能資源數(shù)值模擬研究。