高雄飛，劉元會(huì)，郭建青，王媛英

(1.長(zhǎng)安大學(xué)理學(xué)院，陜西西安 710064;2.長(zhǎng)安大學(xué)環(huán)境科學(xué)與工程學(xué)院，陜西西安 710051)

土壤含水量是土壤的基本物理性質(zhì)指標(biāo)之一。準(zhǔn)確掌握土壤水分的動(dòng)態(tài)變化是確定灌溉制度并適時(shí)進(jìn)行灌溉的重要依據(jù)［1］。對(duì)土壤含水量時(shí)間序列的時(shí)空演變過(guò)程進(jìn)行研究在土壤含水量預(yù)測(cè)預(yù)報(bào)中具有重要的作用。目前，處理這類長(zhǎng)期觀測(cè)資料的方法主要有回歸分析法［2］、頻譜分析法［3］、隨機(jī)微分方程法［4］等。這些方法的結(jié)果能夠在一定程度上反映土壤動(dòng)態(tài)的內(nèi)在物理機(jī)制，但是，它們均未能很好地揭示土壤含水量動(dòng)態(tài)過(guò)程的分形特征［5］。土壤含水量由于受到降雨、蒸發(fā)、日照時(shí)間等因素的影響，具有高度的非線性，對(duì)于這類非線性的具有統(tǒng)計(jì)特性的數(shù)據(jù)系列，采用重標(biāo)極差分析法［6］(rescaled range analysis，簡(jiǎn)稱R/S分析法)，可以較好地揭示土壤含水量動(dòng)態(tài)變化過(guò)程的內(nèi)在規(guī)律及其分形特征，從而有效預(yù)報(bào)或控制其變化規(guī)律。

Hurst［7］在總結(jié)尼羅河的多年水文觀測(cè)資料時(shí)提出了R/S分析法，其后Mandelbrot［8］對(duì)該方法在理論上進(jìn)行了補(bǔ)充和完善。Wallis等［9］曾利用此法對(duì)河流流量、泥漿沉積量、樹(shù)木年輪、降雨量等自然現(xiàn)象進(jìn)行研究。該方法在其他領(lǐng)域(如物理、生物學(xué)等)也得到了廣泛的應(yīng)用［10］。本文以欒城水文站2002—2008年深度為10 cm、80 cm和150 cm土壤含水量時(shí)間序列為例，采用R/S分析法對(duì)土壤含水量的動(dòng)態(tài)過(guò)程進(jìn)行分析，進(jìn)而借助分形理論對(duì)土壤含水量的分形特征進(jìn)行定量描述，期望揭示當(dāng)?shù)剞r(nóng)田土壤水分的動(dòng)態(tài)變化特征。

1 R/S分析法

1.1 R/S分析的基本原理

將時(shí)間序列{Xi=1，2，3，…，N}劃分為m個(gè)時(shí)間長(zhǎng)度為n的獨(dú)立時(shí)間序列，其中m=N/n(m為整數(shù))，子序列記為Da，其中下標(biāo)a表示第a個(gè)子序列［11］。

a.每個(gè)時(shí)間長(zhǎng)度為n的子序列Da的均值計(jì)算公式為

式中:Xˉ——子序列平均值;j——子序列下標(biāo)。

b.每個(gè)子序列Da的累積偏差為

式中:Xj，a——子序列的第j個(gè)觀測(cè)值;Xˉa——第a個(gè)子序列的平均值。

c.每個(gè)子序列Da的極差為

d.每個(gè)子序列Da的標(biāo)準(zhǔn)差為

e.每個(gè)子序列的重標(biāo)極差定義為

f.對(duì)每個(gè)子序列重復(fù)第b步到第e步的計(jì)算，得到一個(gè)重標(biāo)極差序列，計(jì)算該序列的均值:

g.將n增加到下一個(gè)更大的值，重復(fù)第a步至第f步的計(jì)算，直到，得到一個(gè)序列。假設(shè)該序列累積偏差的極差與觀測(cè)數(shù)的一個(gè)待定冪次方成正比［12］，即

其中

式中:τ——時(shí)間步長(zhǎng);C——待定常數(shù)，可以理解為反映自然現(xiàn)象固有特性的常數(shù);H——時(shí)間序列的Hurst指數(shù)。

對(duì)式(7)兩邊取對(duì)數(shù)，有

式(8)利用最小二乘法進(jìn)行求解，其斜率H即為Hurst指數(shù)。

H表示了時(shí)間變換的標(biāo)度因子，因此時(shí)間序列具有自仿射性，經(jīng)過(guò)嚴(yán)格的數(shù)學(xué)推導(dǎo)［13］，時(shí)間序列的分形維數(shù)與H的關(guān)系［14-15］如下:

式中:D——分形維數(shù)，D定量地描述了土壤含水量動(dòng)態(tài)所具有的分形特征，可以表征不同深度土壤含水量變化過(guò)程的復(fù)雜程度，D越大說(shuō)明變化過(guò)程越復(fù)雜［16］。

1.2 Hurst指數(shù)的意義

自相似性和長(zhǎng)程依賴性已成為時(shí)間序列的重要特征［17-18］，定量描述長(zhǎng)程依賴性的主要方法就是估計(jì)Hurst指數(shù)。Hurst指數(shù)是刻畫時(shí)間序列的重要參數(shù)，依據(jù)分形理論，時(shí)間序列的Hurst指數(shù)介于0～1之間。以0.5為間隔，時(shí)間序列在不同的區(qū)間表現(xiàn)出不同的特性［19］:(a)如果H=0.5，則該時(shí)間序列為獨(dú)立同分布的隨機(jī)序列，時(shí)間序列在各個(gè)時(shí)間尺度上都是相互獨(dú)立的;如果H≠0.5，則可判斷該時(shí)間序列為非隨機(jī)序列。(b)如果0.5＜H＜1，表明該時(shí)間序列是持續(xù)性的，或稱其具有長(zhǎng)期記憶性(也就是說(shuō)，若時(shí)間序列在過(guò)去是上升(下降)的，現(xiàn)在或?qū)?lái)時(shí)期將繼續(xù)上升(下降))。H越接近1，時(shí)間序列的持續(xù)性程度就越強(qiáng)，反之就越弱。如果0＜H＜0.5，表示時(shí)間序列存在反持續(xù)性，也就是說(shuō)，該時(shí)間序列在過(guò)去是下降趨勢(shì)就意味著將來(lái)是上升趨勢(shì)，反持續(xù)行為的強(qiáng)度依賴于H接近0的程度。

2 實(shí)例應(yīng)用

2.1 基本數(shù)據(jù)

本例數(shù)據(jù)引自文獻(xiàn)［20］。該文獻(xiàn)給出了2002—2008年中國(guó)科學(xué)院欒城農(nóng)業(yè)生態(tài)系統(tǒng)試驗(yàn)站內(nèi)試驗(yàn)農(nóng)田15個(gè)不同深度層次的土壤含水量實(shí)測(cè)數(shù)據(jù)、觀測(cè)點(diǎn)的地理位置及觀測(cè)方法。為了探討不同深度土壤含水量的動(dòng)態(tài)變化特征，筆者選用了深度為10 cm、80 cm和150cm共3層土壤含水量的實(shí)測(cè)數(shù)據(jù)進(jìn)行計(jì)算與分析。深度10 cm代表土壤表層，深度80 cm代表作物根系的主要發(fā)育區(qū)，深度150 cm代表作物根系發(fā)育區(qū)底部。

2.2 土壤含水量的動(dòng)態(tài)分析

將土壤含水量看成是隨機(jī)變量，即t～X(t)，利用Hurst經(jīng)驗(yàn)公式(即式(8))對(duì)其進(jìn)行計(jì)算，并將計(jì)算結(jié)果表示在雙對(duì)數(shù)坐標(biāo)中，再利用最小二乘法對(duì)點(diǎn)進(jìn)行擬合，得到的直線斜率即是H的估計(jì)值。序列的D可以通過(guò)式(9)間接求得。該時(shí)間序列長(zhǎng)度N=379，為了研究不同層次土壤含水量時(shí)間序列的動(dòng)態(tài)變化特征，分別取不同的時(shí)間尺度N=100、200、379，通過(guò)Matlab軟件將所得的一系列l(wèi)gτ、lgF(τ)在雙對(duì)數(shù)坐標(biāo)上進(jìn)行擬合，從而確定直線的斜率，結(jié)果見(jiàn)圖1、圖2和圖3。

圖1 深度為10 cm不同時(shí)間尺度土壤含水量R/S分析結(jié)果Fig.1 Results of R/S analysis of soil moisture on different time scales at depth of 10 cm

表1給出了深度為10 cm、80 cm和150 cm土壤含水量時(shí)間序列在不同時(shí)間尺度上的H和D。從表1可以看出，在不同的時(shí)間尺度下，H均大于0.5，由此可以判斷當(dāng)?shù)夭煌疃韧寥篮康臅r(shí)間序列是非隨機(jī)時(shí)間序列，時(shí)間序列存在顯著的長(zhǎng)程依賴性，而且該時(shí)間序列是持續(xù)性的，即同一深度含水量時(shí)間序列總體呈現(xiàn)下降或上升的趨勢(shì)。從圖4對(duì)3個(gè)土壤深度層次變化曲線擬合的直線可以看出，深度為10 cm處含水量呈現(xiàn)上升趨勢(shì)，深度為80 cm和150 cm處含水量呈現(xiàn)下降趨勢(shì)。產(chǎn)生這一現(xiàn)象的原因可能是隨著當(dāng)?shù)亟?jīng)濟(jì)發(fā)展速度加快，過(guò)量開(kāi)采地下水，導(dǎo)致地下水持續(xù)下降，引起深層土壤含水量呈現(xiàn)出減小的趨勢(shì)。然而，由于受到降雨補(bǔ)充或灌溉蒸發(fā)的消耗，表層土壤含水量的變化幅度較大，呈現(xiàn)上升的趨勢(shì)，這也與圖4所示的結(jié)果基本相符。

深度為10cm、80 cm和150 cm土壤含水量時(shí)間序列的平均分形維數(shù)D分別為1.262、1.123、1.047，D隨著深度的增加越來(lái)越小，由此可知，隨著土壤深度的增加，土壤含水量波動(dòng)幅度越來(lái)越小，噪聲更少。產(chǎn)生這一現(xiàn)象的原因可能是表層土壤含水量的大小受降水和蒸發(fā)的交互影響較強(qiáng)，土壤含水量變化幅度較大，波動(dòng)程度較高，故分形維數(shù)也較大;對(duì)深層土壤含水量而言，由于受降雨與蒸發(fā)的影響較弱，土壤含水量變化幅度較小，波動(dòng)程度較低，故分形維數(shù)減小。另一方面，土壤含水量時(shí)間序列的長(zhǎng)程依賴性隨著深度的增加也越來(lái)越弱，并得出當(dāng)?shù)赝寥篮康臅r(shí)間序列為分形時(shí)間序列的結(jié)論。

圖2 深度為80 cm不同時(shí)間尺度土壤含水量R/S分析結(jié)果Fig.2 Results of R/S analysis of soil moisture on different time scales at depth of 80 cm

圖3 深度為150 cm不同時(shí)間尺度土壤含水量R/S分析結(jié)果Fig.3 Results of R/S analysis of soil moisture on different time scales at depth of 150 cm

表1 不同土壤深度在不同時(shí)間尺度的相關(guān)參數(shù)Table 1 Parameters at different depths on different time scales

3 結(jié) 語(yǔ)

借助分形理論中的R/S分析法，對(duì)欒城縣2002—2008年深度為10 cm、80 cm和150 cm土壤含水量時(shí)間序列進(jìn)行動(dòng)態(tài)分析，得到如下結(jié)論:

a.在深度為10cm、80cm和150cm處土壤含水量時(shí)間序列的Hurst指數(shù)H分別為0.733、0.877、0.953，均大于0.5。不同層次的土壤含水量表現(xiàn)出長(zhǎng)程的依賴性，也就是說(shuō)，同一土壤層含水量在整個(gè)序列上總體呈現(xiàn)下降或上升的趨勢(shì)，而且該序列是非隨機(jī)序列。

圖4 不同深度土壤含水量的變化趨勢(shì)Fig.4 Trend of soil moisture dynamics at different depths

b.在10 cm、80 cm和150 cm深度，土壤含水量序列的平均分形維數(shù)D分別為1.262、1.123、1.047，其分形維數(shù)逐漸變小，說(shuō)明土壤的含水量波動(dòng)幅度也越來(lái)越小，噪聲更少，其長(zhǎng)程依賴性也越來(lái)越弱，另一方面也可以判斷出欒城區(qū)土壤含水量的時(shí)間序列為分形時(shí)間序列。

［1］李松敏，王仰仁.考慮下界面水分通量的棉田土壤水分模擬研究［J］.人民黃河，2010，32(11):89-91.(LI Songmin，WANG Yangren.Consider the interface moisture flux simulation study of soil moisture in cotton fields［J］.Yellow River，2010，32(11):89-91.(in Chinese))

［2］劉政，孫艷，李保成，等.SPSS逐步回歸分析法構(gòu)建北疆棉蚜年發(fā)生程度的預(yù)報(bào)模型［J］.安徽農(nóng)業(yè)科學(xué)，2009，37(6):2815-2817.(LIU Zheng，SUN Yan，LI Baocheng，et al.Stablishment of prediction model on year occurrence degree ofaphidsgossipier in North Xinjiang with SPSSstepwise regression analysis［J］.Journal of Anhui Agricultural Sciences，2009，37(6):2815-2817.(in Chinese))

［3］于福榮，王友賀.頻譜分析法在鄭州市降水量預(yù)報(bào)中的應(yīng)用［J］.人民長(zhǎng)江，2012，43(10):50-53.(YU Furong，WANG Youhe.Application of frequency spectrum analysis in annual precipitation forecast for Zhengzhou City ［J］.Yangtze River，2012，43(10):50-53.(in Chinese))

［4］劉宇峰，孫虎，原志華.基于小波分析的汾河河津站徑流與輸沙的多時(shí)間尺度特征［J］.地理科學(xué)，2012，32(6):764-770.(LIU Yufeng，SUN Hu，YUAN Zhihua.Multiple time scale features of run off and sediment discharge on wavelet analysis at Hejin Station of Fenhe River Basin［J］.Scientia Geographical Sinica，2012，32(6):764-770.(in Chinese))

［5］桑燕芳，王中根，劉昌明.水文時(shí)間序列分析方法研究進(jìn)展［J］.地理科學(xué)進(jìn)展，2013，32(1):20-30.(SANG Yanfang，WANG Zhonggeng，LIU Changming.Hydrological time series analysis method is reviewed［J］.Progress in Geography，2013，32(1):20-30.(in Chinese))

［6］蔡江東.重標(biāo)極差分析法預(yù)測(cè)徑流量問(wèn)題探究［J］.人民黃河，2012，34(7):25-27.(CAI Jiangdong.Approach to the problems of rescaled range analysis in predicting the runoff［J］.Yellow River.2012，34(7):25-27.(in Chinese))

［7］HURST H E.Long-term storage capacity of reservoirs［J］.Transactions of the American Society of Civil Engineers，1951，116:770-808.

［8］MANDELBROT B B.Fractional brownian motion，fractional noise and Application［J］.SIAM Rev，1968(10):425-437.

［9］WALLISJ R，MATALAS N C.Small sample properties ofHandK-estimations of the Hurst coefficient［J］.Water Resource Research，1970，6(6):1583-1594.

［10］劉志，姚運(yùn)生.鄂西及鄰區(qū)地震時(shí)空叢集的物理分形及R/S分析［J］.大地測(cè)量與地球動(dòng)力學(xué)，2012，32(1):13-16.(LIU Zhi，YAO Yunsheng.Western Hubei and its adjacent area earthquake cluster of time and space physics and fractal R/Sanalysis［J］.Journal of Geodesy and Geodynamics，2012，32(1):13-16.(in Chinese))

［11］吳克金，汪波.基于多重分形與混沌理論的金融市場(chǎng)研究［D］.天津:天津大學(xué)，2005.

［12］王斌.基于分形理論的我國(guó)股票市場(chǎng)價(jià)格波動(dòng)研究［D］.西安:西安電子科技大學(xué)，2010.

［13］屈世顯，張建華.復(fù)雜系統(tǒng)的分形理論與應(yīng)用［M］.西安:陜西人民出版社，1995:344-358.

［14］黃勇，周志芳，王錦國(guó)，等.R/S分析在地下水動(dòng)態(tài)分析中的應(yīng)用［J］.河海大學(xué)學(xué)報(bào):自然科學(xué)版，2002，30(1):83-87.(HUANGYong，ZHOU Zhifang，WANGJinguo，et al.R/Sanalysis in the application of the groundwater dynamic analysis［J］.Journal of Hohai University:Natural Sciences，2002，30(1):83-87.(in Chinese))

［15］趙晶，王乃昂.近50年來(lái)蘭州城市氣候變化的 R/S分析［J］.干旱區(qū)地理，2002，25(1):90-95.(ZHAO Jing，WANG Naiang.Nearly 50 years of Lanzhou City of R/Sanalysis of climate change［J］.Arid Land Geography，2002，25(1):90-95.(in Chinese))

［16］王小軍，張強(qiáng)，古璇清.基于分形理論的灌溉水有效利用系數(shù)空間尺度變異［J］.地理學(xué)報(bào)，2012，67(9):1201-1212.(WANG Xiaojun，ZHANG Qiang，GU Xuanqing.Fractal-based effective utilization coefficient of irrigation water space scale variability［J］.Acta Geographica Sinica，2012，67(9):1201-1212.(in Chinese))

［17］MESA O J，POVEDA G.The Hurst effect:the scale of fluctuation approach［J］.Water Resour Res，1993，29(12):3995-4002.

［18］SHAO Quanxi，LI Ming.A new trend analysis for seasonal time series with consideration of data dependence［J］.Journal of Hydrology，2011，17:104-112.

［19］MANDELBROT B B，TOQUE M S.Robust R/S analysis of long-rain serial correlation［R］.Manila:Bulletin of the International Statistical Institute，1979:69-99.

［20］孫鴻烈，胡春勝.農(nóng)田生態(tài)系統(tǒng)卷河北欒城站(1988—2008)［M］.北京:中國(guó)農(nóng)業(yè)出版社，2011:120-128.