李 娜，吳 娟，李欣昱

(1.綏化學(xué)院 農(nóng)業(yè)與水利工程學(xué)院，黑龍江 綏化 152000；2.黑龍江工商學(xué)院，黑龍江 哈爾濱 150016)

1 研究區(qū)概況

克山農(nóng)場隸屬于黑龍江北大荒農(nóng)墾集團(tuán)，位于黑龍江西北部，訥河市與克山縣境內(nèi)。地理位置為125°07′40″～125°37′30″E，48°11′15″～48°24′07″N，海拔240～340 m。土地面積3.51萬hm2，耕地面積2.72萬hm2，人口2.7萬人(2019年)，下設(shè)5個(gè)管理區(qū)，27個(gè)居民點(diǎn)，場直科研管理區(qū)和種子加工廠。

克山農(nóng)場地處小興安嶺西麓，松嫩平原東北部，地勢(shì)丘陵漫崗，土質(zhì)肥沃，黑鈣土占土地面積的98%，適宜小麥、大豆、玉米等農(nóng)作物生長。根據(jù)農(nóng)作物品種積溫區(qū)劃，克山農(nóng)場屬于溫帶氣候第三積溫帶，春季降雨少、多風(fēng)，夏季溫度高、降雨多，秋季溫度迅速降低，初霜期在9月18日左右，冬季時(shí)間持續(xù)較長、溫度低、降雪量較大。該農(nóng)場近十年年平均氣溫達(dá)2.1 ℃，明顯高于多年年平均氣溫；極端高溫達(dá)到36.5 ℃，極端低溫在-37.6 ℃左右；全年無霜期大約為120 d。多年平均年降水量為502.5 mm，6—8月年降水量占比達(dá)到68.3%，9—11月降水量占比16.4%，3—5月降水量占比11.0%左右，12月—次年1月降水量占比約3.0%。年平均風(fēng)速3.9 m/s，常年以西北風(fēng)為主。

2 研究方法

2.1 滑動(dòng)平均法

滑動(dòng)平均法是時(shí)間序列分析中常用的一種動(dòng)態(tài)數(shù)據(jù)處理方法。通過選取合適的平均值來減小數(shù)據(jù)波動(dòng)對(duì)研究結(jié)果的影響，進(jìn)而減少實(shí)驗(yàn)誤差。更加直觀的呈現(xiàn)時(shí)間序列變化的趨勢(shì)性和階段性。

對(duì)時(shí)間序列x1，x2，…，xn的若干個(gè)前期值和后期值取平均值，得出新的序列，這種處理方法消除了異常值對(duì)序列的影響，使原序列更加光滑，如式(1)：

yt=(xt-1+xt-2+xt-3+…+xt-n)/n

(1)

式中：yt為新的時(shí)間序列；xt-1為前期實(shí)際值；xt-2，xt-3和xt-n分別為前兩期、前三期直至前n期的實(shí)際值；n為移動(dòng)平均的時(shí)期個(gè)數(shù)。

n的取值具有關(guān)鍵性作用，它決定了滑動(dòng)平均所取前期值和后期值的個(gè)數(shù)。本文取n等于5，即滑動(dòng)平均為5年滑動(dòng)平均[1]。

2.2 線性趨勢(shì)法

將樣本量為n的某氣候變量用xi表示，其所對(duì)應(yīng)時(shí)間用ti表示，建立某氣候變量xi隨時(shí)間ti變化的線性回歸方程來描述xi與ti關(guān)系，如式(2)所示。

xi=a+bti(i=1，2，…，n)

(2)

式中：xi為某氣候變量某時(shí)刻狀態(tài)；a為回歸常數(shù)；b為回歸系數(shù)；ti為某一時(shí)刻；n為樣本量。對(duì)不同時(shí)間t所觀測(cè)的氣候變量數(shù)據(jù)的回歸常數(shù)a和回歸系數(shù)b，利用最小二乘法進(jìn)行計(jì)算，如式(3)～式(6)：

(3)

(4)

(5)

(6)

2.3 Mann-Kendall非參數(shù)檢驗(yàn)

假設(shè)有一時(shí)間序列x1，x2，…，xn構(gòu)造秩序列ri，ri表示當(dāng)1≤j≤i時(shí)的樣本累積數(shù)。如式(7)、式(8)：

(7)

(8)

式中：sk為統(tǒng)計(jì)量，當(dāng)時(shí)間系列大于某一長度時(shí)，統(tǒng)計(jì)量sk服從正態(tài)分布；n為樣本量。sk均值E(st)與方差Var(st)，如式(9)、式(10)：

(9)

(10)

假定時(shí)間的序列隨機(jī)獨(dú)立，對(duì)統(tǒng)計(jì)量UFk做如下定義，如式(11)：

(11)

式中：UFk由sk、均值E(sk)和方差Var(sk)組成。當(dāng)UFk=0時(shí)，UFk呈現(xiàn)標(biāo)準(zhǔn)正態(tài)分布，根據(jù)研究需要，可以給定k不同的顯著水平，假定給定的某一顯著水平α，通過查詢正態(tài)分布表得到統(tǒng)計(jì)量臨界值，用Uα表示，當(dāng)|UFk|>Uα表明序列存在一個(gè)明顯趨勢(shì)變化，UFk在圖中表示為UF。將上式利用到反序列中，重復(fù)上述計(jì)算過程，UBk的值可以用計(jì)算值乘以-1得到，圖中UB表示UBk，令UBi=0。

進(jìn)一步利用序列中UFk和UBk的值分析時(shí)間序列x變化趨勢(shì)，其中UFk曲線與UBk曲線交點(diǎn)為時(shí)間序列發(fā)生突變的時(shí)間，以此指出突變區(qū)域。當(dāng)UFk>0時(shí)，說明序列呈上升趨勢(shì)；當(dāng)序列呈下降趨勢(shì)時(shí)，UFk< 0。

水文研究的某個(gè)氣候變量時(shí)間序列的突變點(diǎn)往往是研究的重點(diǎn)，因?yàn)橥蛔凕c(diǎn)可以說明很多問題，比如：導(dǎo)致其突變的氣候原因、地理因素以及人類活動(dòng)等，以往對(duì)水文時(shí)間序列研究僅僅是進(jìn)行趨勢(shì)性分析，但是趨勢(shì)性分析只能說明時(shí)間序列的整個(gè)變化趨勢(shì)，時(shí)間序列內(nèi)部的變化特征無法具體反映出來，因此本文采用M-K方法更加深入分析了研究區(qū)域降水量時(shí)間序列內(nèi)部的突變特征[3]。

3 降水量時(shí)間特征分析

3.1 年際變化特征

以克山農(nóng)場1991—2020年逐年、月降水量數(shù)據(jù)為依據(jù)，采用趨勢(shì)性分析和M-K突變檢驗(yàn)特征分析方法對(duì)降水量進(jìn)行時(shí)間特征分析。

(1)趨勢(shì)性分析。應(yīng)用克山農(nóng)場逐年降水量數(shù)據(jù)，繪出年降水量變化特征圖，通過分析可知：克山農(nóng)場1991—2020年年降水量整體呈下降的趨勢(shì)。應(yīng)用滑動(dòng)平均法消除異常值影響，保證結(jié)果的準(zhǔn)確性，對(duì)年降水量數(shù)據(jù)現(xiàn)進(jìn)行5 a滑動(dòng)平均處理，結(jié)果表明：近30 a來克山農(nóng)場降水量呈現(xiàn)逐年下降趨勢(shì)，大體上呈波動(dòng)狀態(tài)。如圖1所示。

圖1 年降水量總體變化特征

(2)突變性分析。給定UFk為0.05的顯著水平，采用M-K方法分析克山農(nóng)場年降水量時(shí)間序列的突變特征，M-K突變檢驗(yàn)分析結(jié)果見圖2。

圖2 年降水量M-K突變分析

克山農(nóng)場UF統(tǒng)計(jì)量、UB統(tǒng)計(jì)量曲線(UF-UB曲線)有1992年、1994年、1995年、1998年、2000年、2005年和2008年等多個(gè)交點(diǎn)，即在1992年、1994年、1995年、1998年、2000年、2005年和2008年降水量發(fā)生突變，表現(xiàn)為克山農(nóng)場年降水量1992年、1994年和1995年后呈下降的趨勢(shì)，但變化不顯著；在1998年、2000年后降水量逐年增加，在2004年達(dá)到最大值，后開始下降，在2008年之后呈下降趨勢(shì)。

3.2 年內(nèi)變化特征

選取克山農(nóng)場1991—2020年30 a克山農(nóng)場逐月降水?dāng)?shù)據(jù)為分析樣本，采用趨勢(shì)分析和M-K突變檢驗(yàn)特征分析方法對(duì)降水量進(jìn)行時(shí)間特征分析。

(1)趨勢(shì)性分析。通過降水量變化特征圖可以直觀看出：在年際變化上，近30a來克山農(nóng)場春季、冬季年際傾向率為正值，降水量均呈現(xiàn)上升的趨勢(shì)；夏季、秋季年際傾向率為負(fù)值，降水量均呈下降的趨勢(shì)。為了降低周期性變化與異常值對(duì)研究結(jié)果對(duì)影響，將數(shù)據(jù)運(yùn)用5 a滑動(dòng)平均方法進(jìn)行處理，結(jié)果表明1991—2020年春季、冬季降水量呈現(xiàn)逐年增加趨勢(shì)，夏季、秋季降水量具有相反的特征，表現(xiàn)為逐年減少的趨勢(shì)。如圖3～圖6 所示。

圖5 秋季降水量變化特征

圖6 冬季降水量變化特征

(2)突變性分析。給定UFk為0.05的顯著水平，采用M-K方法深入分析克山農(nóng)場降水量時(shí)間序列的突變特征，M-K突變檢驗(yàn)分析結(jié)果，如圖7～圖10所示。

圖7 春季降水量M-K突變分析

如圖7所示，春季UF-UB曲線有2015年和2017年2個(gè)交點(diǎn)，即在2015年和2017年降雨量發(fā)生突變，表現(xiàn)為春季的降水量在2015年后呈下降趨勢(shì)，在2017年后降水量呈上升趨勢(shì)，但是均為不顯著變化；如圖8所示，夏季UF-UB曲線有1993年、2005年、2008年和2018年4個(gè)交點(diǎn)，即在1993年、2005年、2008年和2018年降水量發(fā)生突變，具體表現(xiàn)為1993年后降水量呈上升趨勢(shì)，2005年和2008年后年降雨量均呈現(xiàn)下降趨勢(shì)，在2018年后呈現(xiàn)上升趨勢(shì)，但各年降水量均為不顯著變化；如圖9所示，秋季UF-UB曲線有2005年、2008年、2010年、2013年和2014年5個(gè)交點(diǎn)，說明秋季降水量在2005年、2008年、2010年、2013年和2014年發(fā)生突變，表現(xiàn)為降雨量在2005年、2008年和2010年后降水量均呈下降趨勢(shì)，2013年降水量呈現(xiàn)上升趨勢(shì)，2014年后呈下降趨勢(shì)，各年均為不顯著變化；如圖10所示，冬季UF-UB曲線有2018年1個(gè)交點(diǎn)，即在2008年降雨量發(fā)生突變，具體表現(xiàn)為2018年以后降雨量呈現(xiàn)上升趨勢(shì)，且變化不顯著。

圖8 夏季降水量M-K突變分析

圖9 秋季降水量M-K突變分析

圖10 冬季降水量M-K突變分析

4 結(jié) 論

(1)降水量年際變化趨勢(shì)?？松睫r(nóng)場UF-UB曲線雖有多個(gè)交點(diǎn)，但UF曲線都在非顯著變化范圍內(nèi)，表現(xiàn)為雖在各個(gè)交點(diǎn)年份降水量發(fā)生突變，但均為不顯著變化，1991—2020年間表現(xiàn)為不平穩(wěn)的波動(dòng)狀態(tài)。

(2)降水量季節(jié)變化趨勢(shì)。春季的降水量在2015年后呈下降趨勢(shì)，在2017年后降水量呈上升趨勢(shì)，但是均為不顯著變化；夏季降水量1993年后呈上升趨勢(shì)，2005年和2008年后年均呈現(xiàn)下降趨勢(shì)，在2018年后呈現(xiàn)上升趨勢(shì)，但各年降水量均為不顯著變化；秋季降水量在2005年、2008年和2010年后均呈下降趨勢(shì)，2013年降水量呈現(xiàn)上升趨勢(shì)，2014年后呈下降趨勢(shì)，各年均為不顯著變化；冬季降水量2018年以后呈現(xiàn)上升趨勢(shì)，且變化不顯著。