南水北調(diào)東線典型受水區(qū)降雨特征分析

張曉松 曹命凱 丁艷霞 王少波 張艷 孫敏

摘要：根據(jù)南水北調(diào)東線典型受水區(qū)1951-2012年降雨資料，采用滑動平均法、MannKendall檢驗法、游程檢驗法、Spearman秩次相關檢驗法及Morlet小波變換分析方法，分析了江蘇省沂南水利分區(qū)降雨演變趨勢和周期變化特征。研究結(jié)果表明，近62年來，該區(qū)域降雨量呈不顯著減少趨勢，并且呈現(xiàn)27 a顯著周期變化。研究結(jié)果對該地區(qū)應對防洪抗旱、農(nóng)業(yè)生產(chǎn)以及南水北調(diào)東線可持續(xù)調(diào)水具有重要的指導意義。

關鍵詞：南水北調(diào)；沂南地區(qū)；降雨量；Morlet小波變換；降雨特征

中圖分類號：TV213文獻標志碼：A文章編號：16721683（2018）03005906

Rainfall characteristics in typical waterreceiving area of East Route of

SouthtoNorth Water Diversion Project

ZHANG Xiaosong1，CAO Mingkai1，DING Yanxia1，WANG Shaobo1，ZHANG yan1，SUN Ming2

（1.Jiangsu Surveying and Design Institute of Water Resources，Co.，Ltd.，Yangzhou 225127，China；

2.Gaoyou Water Conservancy Bureau，Yangzhou 225600，China）

Abstract：On the basis of the rainfall data from 1951 to 2012，in this paper we analyzed the rainfall evolution trend in a typical waterreceiving area of East Route of SouthtoNorth Water Diversion Project.We adopted four methods including the moving average method，MannKendall test method，Run test method，and Spearman rank test.Meanwhile we used Morlet wavelet transform to analyze the periodic variation characteristics.The results showed that the rainfall in the region showed an insignificant decreasing trend from 1951 to 2012 and had a 27year variation cycle.The results have important guiding significance for flood control，drought resisting，agricultural production，and sustainable water diversion for the East Route of SouthtoNorth Water Diversion Project.

Key words：SouthtoNorth Water Diversion；Yinan area；rainfall；Morlet wavelet transform；rainfall characteristics

南水北調(diào)工程是解決我國北方地區(qū)水資源短缺，實現(xiàn)社會經(jīng)濟可持續(xù)發(fā)展的特大型跨流域調(diào)水工程，分別從長江上、中、下游調(diào)水，經(jīng)東、中、西三線輸送至水資源短缺的黃淮海流域。南水北調(diào)東線工程利用江蘇省已建的江水北調(diào)工程，逐步擴大調(diào)水規(guī)模并延長輸水線路，工程從長江下游揚州江都水利樞紐抽引長江水，利用京杭大運河及與其平行的河道逐級提水北送至黃淮海平原東部和膠東地區(qū)，其中江蘇境內(nèi)受水區(qū)涉及揚州、淮安、鹽城、宿遷、連云港和徐州六市，總面積381萬km2。

隨著全球區(qū)域氣候的變化，將導致區(qū)域降雨量發(fā)生變化，從而對水資源供給和自然生態(tài)系統(tǒng)產(chǎn)生重大影響。因此，分析研究區(qū)域降雨變化的趨勢和周期變化規(guī)律[14]，對認識區(qū)域降雨特征、應對強降雨等極端氣候以及農(nóng)業(yè)生產(chǎn)、防汛抗旱、水資源規(guī)劃等具有重要意義。南水北調(diào)工程實施以來，水文工作者進行了很多有價值的研究，但是有關降雨特征分析的研究僅南水北調(diào)中線京石段開展過[5]。目前研究降雨變化趨勢常用的方法有滑動平均法[67]、線性估計、MannKendall秩次相關檢驗法[89]等，常用的周期變化研究方法有回歸周期分析法、傅立葉分析法[10]、小波分析法等。本文以南水北調(diào)東線江蘇境內(nèi)受水區(qū)的沂南水利分區(qū)作為研究區(qū)域，利用滑動平均法、MannKendall秩次相關檢驗法、游程檢驗法[1112]、Spearman秩次相關檢驗法[1315]四種方法相結(jié)合對該地區(qū)降雨變化趨勢進行分析，并應用Morlet小波變換[1112]分析降雨的時間及周期變化特征，此研究為南水北調(diào)東線江蘇境內(nèi)水量調(diào)度運行管理、區(qū)域抗旱排澇調(diào)度工作的開展以及水資源開發(fā)與利用等提供參考。

第16卷 總第96期·南水北調(diào)與水利科技·2018年6月張曉松等·南水北調(diào)東線典型受水區(qū)降雨特征分析1區(qū)域概況

江蘇省沂南水利分區(qū)（簡稱：沂南地區(qū)）是南水北調(diào)東線江蘇境內(nèi)水量配置模型的計算單元之一，位于新沂河以南、廢黃河以北，東臨黃海、西至中運河，涉及四市七縣區(qū)，包括鹽城市的響水縣，連云港市的灌南縣，淮安市的漣水縣、淮陰區(qū)，宿遷市的宿豫區(qū)、泗陽縣、沭陽縣，區(qū)域總面積7 477 km2。區(qū)域地形總趨勢是西高東低，自西南向東北傾斜，東西向平均地面比降約萬分之一點三。該區(qū)域地處亞熱帶向暖溫帶過渡性氣候中，又是南北氣流過渡地帶，有明顯的海洋季風特征，全年四季分明，夏季暖熱多雨，受臺風侵襲頻繁，冬季寒冷干燥，春秋兩季晝暖夜寒，冷空氣活動頻繁，溫差變化較大。

2暴雨類型及特征

區(qū)域降雨相對比較集中，多突發(fā)性暴雨，最大日降雨量6103 mm（陳港站，2000年8月30日）。降雨年際變化幅度大，多年平均面雨量為92612 mm（1951-2012年），最大年降雨量為1 48969 mm（2003年），最小年降雨量僅為57224 mm（2004年），多年平均年降雨天數(shù)992 d，最多年降雨天數(shù)為131 d（1954年），最少年降雨天數(shù)63 d（1957年）；降雨年內(nèi)分配極不均勻，汛期6月-9月降雨量比較集中，約占全年降雨量的70%，秋末的10月至次年的2月，降雨量一般只占全年降雨量的15%左右。

臺風則多發(fā)生在7月到9月，尤其以8月最頻繁，常因8～9級大風大潮伴隨著100 mm以上降雨同時發(fā)生，極易產(chǎn)生突發(fā)性的洪、澇、潮等自然災害。另龍卷風、六級以上的大風、冰雹和5月中下旬的干熱風，在局部地區(qū)也會造成危害[1617]。

3資料選取及分析方法

3.1資料選取

根據(jù)沂南地區(qū)雨量站網(wǎng)分布，并綜合考慮降雨資料的代表性、可靠性以及完整性，選用區(qū)域內(nèi)的宿遷閘、泗陽閘、淮陰閘等37個雨量站點1951-2012年逐日降雨數(shù)據(jù)，利用泰森多邊形法[18]進行面平均雨量計算，以面平均雨量作為區(qū)域降雨趨勢分析的資料，能更好更全面地分析該地區(qū)的降雨變化特征。沂南地區(qū)雨量站點分布見圖1。圖1沂南地區(qū)雨量站點分布

Fig.1Location of precipitation stations in Yinan Area

3.2分析方法

本文采用滑動平均法、MannKendall秩次相關檢驗法、游程檢驗法、Spearman秩次相關檢驗法等4種方法來研究分析區(qū)域降雨趨勢。同時應用Morlet小波變換分析該地區(qū)降雨的時間及周期變化特征。

（1）滑動平均法。

在一個系列{xn}中，將連續(xù)的2n或2n+1個數(shù)值取平均，然后向后逐年滑動求均值直到最后一個數(shù)據(jù)，得到的新序列{yn}比原序列{xn}更光滑，即滑動平均法。

（2）MannKendall檢驗法。

應用MannKendall方法檢驗降雨時間序列{xn }的變化趨勢時，構(gòu)造統(tǒng)計量S如下：

S=∑n-1k=1∑nj=k+1Sgn（xj-xk）（1）

式中：

Sgn（xj-xk）=+1（xj-xk）>0

0（xj-xk）=0

-1（xj-xk）<0）（2）

統(tǒng)計量S服從正態(tài)分布，均值為0，方差Var（S）=n（n-1）（2n+5）/18。當n>10時，將統(tǒng)計量S標準化：

U=S-1Var（S）S>0

0S=0

S+1Var（S）S<0（3）

若當U<0，則認為有下降趨勢，反之呈上升趨勢。在顯著水平α下，若|U|≥Uα/2，認為趨勢性顯著，反之趨勢性不顯著[89]。

（3）游程檢驗法。

根據(jù)已知的降雨時間序列{xn }，求出相應的距平系列{x′n}，其中x′i=xi-，=1/n ∑ni=1xi 。若將連續(xù)的x′i>0或x′i<0看作一個游程，統(tǒng)計序列中出現(xiàn)的游程總數(shù)S，其中n1為序列中x′i>0出現(xiàn)的總數(shù)，n2為序列中x′i<0出現(xiàn)的總數(shù)。實際分析中，如n1或n2大于15時，構(gòu)造統(tǒng)計量U：

U=S-E（S）D（S）（4）

構(gòu)造統(tǒng)計量U近似服從N（0，1）標準正態(tài)分布，式中：

E（S）=2n1 n2n1+n2+1（5）

D（S）=2n1 n2 [2n1 n2-（n1+n2）]（n1+n2）2 （n1+n2-1）（6）

根據(jù)游程檢驗的S分布表，選取某一顯著α水平作雙邊檢驗。當|U|>Uα/2時，拒絕原假設，序列趨勢變化明顯，反之不明顯[1112]。

（4）Spearman秩次相關檢驗法。

Spearman秩次相關檢驗主要通過分析降雨序列{xi }與其時序i的相關性而檢驗降雨時序是否有趨勢性。在運算時，將降雨序列按從大到小排列，用其新秩序Ri代表，計算秩次相關系數(shù)：

r=1-6·∑ni=1d2i/（n3-n）（7）

式中：n為序列長度；di=Ri-i。如果秩次與時序相近，則較小，秩次相關系數(shù)較大，趨勢性顯著。通常采用t檢驗法檢驗降雨時序的趨勢性是否顯著，統(tǒng)計量T的計算公式為：

T=r（n-4）/（1-r2）（8）

T服從自由度為n-2的t分布，原假設為序列無趨勢，則根據(jù)降雨時序的秩次相關系數(shù)計算統(tǒng)計量T，然后選擇顯著水平α，在t分布表中查出臨界值tα/2，當|T|>tα/2時，則拒絕原假設，說明序列隨時間有相依關系，從而推斷序列趨勢明顯，否則，接受原假設，趨勢不顯著[1315]。

（5）Morlet小波變換。

小波分析是一種信號的時間/尺度（頻率）分析方法，具有多分辨率分析的特點，在時頻兩域都具有表征信號局部特征的能力，被廣泛應用于水文時間序列的尺度分析、水文預測預報及水文隨機模擬。Morlet小波是一維連續(xù)的，其波形與降雨的波形較為接近。函數(shù)表達式為：

Φ（t）=exp（iw0 t）exp（-t22）（9）

式中：w0為常數(shù)，i為虛數(shù)。通過小波變換后，得到的小波變換圖能夠反映降雨在不同時間尺度上的周期變化[19]。小波系數(shù)實部為正的用實線表示，代表降雨偏多；小波系數(shù)實部為負的用虛線表示，代表降雨偏少。小波方差圖反映了能量隨時間尺度的分布，可以用來分析一個時間序列中各種尺度擾動的相對強度。對應峰值處的尺度稱為該序列的主要時間尺度，即主要周期[2022]。為了進一步確定降雨的周期變化，需通過小波方差檢驗，確定顯著周期。小波方差的計算公式為：

Var（α）=∫∞-∞|Wf（a，b）|2db（10）

4降雨變化特征分析

4.1降雨變化趨勢分析

本文對沂南地區(qū)1951-2012年面平均雨量采用5年滑動平均法，直觀反映該地區(qū)降雨的變化趨勢，同時采用MannKendall秩次相關檢驗法、游程檢驗法、Spearman秩次相關檢驗法，選用置信水平α=005，通過比較檢驗統(tǒng)計量與臨界值大小來判斷降雨變化趨勢的顯著性。沂南地區(qū)年降雨量5年滑動平均過程線見圖2，顯著性檢驗結(jié)果見表1。

圖2可以直觀反映沂南地區(qū)年降雨量的變化情況，降雨呈現(xiàn)增多和減少交替變化，但總體趨勢是減少了，同時通過MannKendall檢驗，U=-0996<0，也定量反映出降雨呈下降趨勢。并且從表1可以看出，三種檢驗一致表明這種下降趨勢是顯著的。

方法檢驗統(tǒng)計量顯著水平臨界值結(jié)果趨勢性MannKendall秩次相關檢驗法U=-0.996α=0.051.96|U|>Uα/2顯著游程檢驗法U=0.39α=0.051.96|U|>Uα/2顯著spearman秩次相關檢驗法T=0.942α=0.052.00|T|

應用Matlab的Morlet小波分析工具對沂南地區(qū)1951-2012年降雨資料進行周期特性分析，得到小波系數(shù)實部等值線圖和小波方差圖，具體見圖3、圖4。

小波系數(shù)實部等值線圖能反映降雨系列不同時間尺度的周期變化及其在時間域中的分布，從而判斷在不同時間尺度降雨的未來變化趨勢。圖3中，實線表示小波系數(shù)正值部分，代表降雨偏多，H表示正值中心；虛線表示小波系數(shù)負值部分，代表降雨偏少，L表示負值中心。從圖中可以清楚的看到降雨演變過程中存在多時間尺度特征，大致為22～32 a，12～21 a，3～11 a的3類尺度的周期變化規(guī)律。其中，22～32 a、12～21 a兩個時間尺度上，周期變化在整個分析時間中表現(xiàn)的比較穩(wěn)定，具有全域性。在22～32 a尺度上先增后減，出現(xiàn)多-少準3次震蕩；在12～21 a時間尺度上，出現(xiàn)多-少6次震蕩。在3～11 a時間尺度上1954-1981年較為明顯，多—少循環(huán)交替。

從圖4小波方差圖可以看出，不同時間尺度擾動的相對強度不同，圖中存在4個較為明顯的峰值，依次27 a、15 a、10 a、4 a。其中27 a對應著最大峰值，說明27 a左右的周期震蕩最強，為沂南地區(qū)降雨變化的第一主要周期；15 a對應著第二峰值，為沂南地區(qū)降雨變化的第二主要周期；10 a對應著第三峰值，為沂南地區(qū)降雨變化的第三主要周期；4 a對應著第四峰值，為沂南地區(qū)降雨變化的第四主要周期。以上4個周期的波動控制著降雨在整個時間域內(nèi)的變化特征。

根據(jù)小波方差檢驗結(jié)果，繪制控制區(qū)域降雨演變的第一、二周期即27 a、15 a特征時間尺度的小波系數(shù)圖，見圖5、圖6。從主周期趨勢圖中可以分析出在不同的時間尺度下，區(qū)域降雨存在的平均周期及多-少變化特征不同。圖5表明，在27 a時間尺度上，平均周期為18 a左右，大約經(jīng)歷了3個多-少轉(zhuǎn)換期。降雨峰值出現(xiàn)在1955年、1973年、1990年和2009年。圖6表明，在15 a時間尺度上，平均周期為10 a左右，大約經(jīng)歷了6個多-少轉(zhuǎn)換期。降雨峰值出現(xiàn)在1953年、1963年、1973年、1982年、1991年、2000年、2011年。

5結(jié)論

（1）在降雨序列具有顯著變化趨勢時，滑動平均法所得到的曲線能夠直觀的反映出變化趨勢，但是不能表明趨勢是否顯著。文中采用3種檢驗方法與滑動平均法相結(jié)合，既可以直觀反映變化趨勢，又可以定量分析檢驗趨勢的顯著性，同時三種檢驗方法可以對比參考。最終分析結(jié)果一致表明沂南地區(qū)1951-2012年面平均降雨量呈現(xiàn)增多和減少交替變化，但總體呈遞減趨勢，并且趨勢顯著。

（2）通過Matlab的Morlet小波分析，沂南地區(qū)降雨存在22～32 a，12～21 a，3～11 a的3類尺度的周期變化規(guī)律。不同的時間尺度表現(xiàn)出不同的循環(huán)交替。小尺度周期震蕩劇烈，大尺度周期震蕩平穩(wěn)，規(guī)律更明顯，全局性更強。

（3）通過小波方差判斷降雨時間序列各種時間尺度擾動的相對強度，沂南地區(qū)降雨受27 a、15 a、10 a、4 a四個時間尺度的干擾較大，其中27 a干擾強度最大，為該序列的主要時間尺度，即主要周期。在27 a時間尺度上，該地區(qū)降雨的平均周期為18 a左右，大約經(jīng)歷了3個多-少轉(zhuǎn)換期。

參考文獻（References）：

[1]潘杰.常州地區(qū)1978年-2012年降雨量時空分布特征[J].南水北調(diào)與水利科技，2015，13（2）：245248.（PAN J.Spatial and temporal characteristics of precipitation in Changzhou Area from 1978 to 2012[J].SouthtoNorth Water Transfers and Water Science & Technology，2015，13（2）：245248.（in Chinese）） DOI：10.13476/j.cnki.nsbdqk.2015.02.012.

[2]徐宗學，張楠.黃河流域近50年降水變化趨勢分析[J].地理研究，2006，25（1）：2734.（XU Z X，ZHANG N.Long term trend of precipitation in the Yellow River basin during the past 50 years[J].Geographical Research，2006，25（1）：2734.（in Chinese）） DOI：10.3321/j.issn：10000585.2006.01.004.

[3]EVANS J P，JAKEMAN A J.Development of a simple ，Catchmentscale，rainfallevaportranspirationrunoff model[J].Environmental Modelling& Software.1998，13：385393.

[4]WAGENER T ，LEES M J.RainfallRunoff Modeling Toolbox user mannual[J].Department of Civil and Environmental Eegineering，Imperial College London，UK.2001.

[5]尹航，王文君，吳英杰，等.南水北調(diào)中線京石段降雨序列特性分析[J].水資源與水工程學報，2016，27（4）：113118.（YIN H，WANG W J，WU Y J，et al.Analysis of rainfall characteristics in BeijingShijiazhuang section of SouthtoNorth Water Diversion[J].Journal of Water Resources & Water Engineering，2016，27（4）：113118 .（in Chinese）） DOI：10.11705/j.issn.1672643X.2016.04.21.

[6]裴益軒，郭民.滑動平均法的基本原理及應用[J].火炮發(fā)射與控制學報，2001（1）：2123.（PEI Y X，GUO M.The fundamental principle and application of sliding average method[J].GUN Launch & Control Journal，2001（1）：2123.（in Chinese）） DOI：10.3969/j.issn.16736524.2001.01.007.

[7]李娟，張維江，馬軼.滑動平均馬爾可夫模型在降雨預測中的應用[J].水土保持研究，2005，12（6）：196198.（LI J，ZHANG W J，MA Y.Application of moving AverageMarkov Chain in the prediction of precipitation[J].Research of Soil and Water Conservation，2005，12（6）：196198.（in Chinese））

[8]康淑媛，張勃，柳景峰，等.基于MannKendall法的張掖市降雨量時空分布規(guī)律分析[J].資源科學，2009，31（3）：501508.（KANG S Y，ZHANG B，LIU J F，et al.Analysis of the spatiotemporal distribution of precipitation in Zhangye City Using MannKendall method[J].Resources Science，2009，31（3）：501508.（in Chinese）） DOI：10.3321/j.issn：10077588.2009.03.022.

[9]郭小芹，劉明春，錢莉，等.從MannKendall特征看石羊河流域降雨量的演變規(guī)律[J].干旱區(qū)地理，2010，33（4）：593599.（GUO X Q，LIU M C，QIAN L，et al.Evolution regularity of precipitation in the Shiyang River Basin on MannKendall features[J].Arid Land Geography，2010，33（4）：593599.（in Chinese）） DOI：10.13826/j.cnki.cn651103/x.2010.04.015.

[10]肖志國.幾種水文序列周期分析方法的比較研究[D].南京：河海大學，2006.（XIAO Z G.Compasison of different analysis methods for the periodicity of hydrological time series[D].Nanjing：Hohai University，2006.（in Chinese））

[11]王文圣，丁晶，李耀清.水文小波分析[M].北京：化學工業(yè)出版社，2005.（WANG W S，DING J，LI Y Q.Hydrology wavelet analysis[M].Beijing：Chemical Industry Press，2005.（in Chinese））

[12]丁晶，劉權(quán)授.隨機水文學[M].北京：中國水利水電出版社，2008.（DING J，LIU Q S.Stochastic hydrology[M].Beijing：China Water Power Press，2008.（in Chinese））

[13]王國慶，賈西安，陳江南，等.人類活動對水文序列的顯著影響干擾點分析——以黃河中游無定河流域為例[J].西北水資源與水工程，2001，12（3）：1315.（WANG G Q，JIA X A ，CHEN J N，et al.Analysis on the transition point of hydrological series impacted by human activities—a case study of Wudinghe basin in the middle reache of the Yellow River[J].Northwest Water Resources & Water Engineering，2001，12（3）：1315.（in Chinese）） DOI：10.3969/j.isnn.1672643X.2001.03.004.

[14]曾小凡，翟建青，姜彤，等.長江流域年降雨量的空間特征和演變規(guī)律分析[J].河海大學學報（自然科學版），2008，36（6）：727732.（ZENG X F，ZHAI J Q，JIANG T，et al.Spatial characteristics and evolutional trends of annual precipitation in the Yangtze River basin[J].Journal of Hohai University（Natural Sciences），2008，36（6）：727732.（in Chinese）） DOI：10.3876/j.issn.10001980.2008.06.001.

[15]杜富慧，郝振純，陳新美，等.邯鄲地區(qū)年降雨量統(tǒng)計特性[J].水資源保護，2009，25（6）：1620.（DU F H，HAO Z C，CHEN X M，et al.Statistical characteristics of annual precipitation in Handan region[J].Water Resources Protection，2009，25（6）：1620.（in Chinese）） DOI：10.3969/j.issn.10046933.2009.06.005.

[16]徐晨光，馮新峰，鄭志偉，等.江蘇省降雨時空分布特征研究[J].現(xiàn)代農(nóng)業(yè)科技，2014，20：236239.（XU C G，F(xiàn)ENG X F，ZHENG Z W，et al.Study on characteristics for spatial and temporal distribution in Jiangsu province[J].Modern Agricultural Sciences and Technology，，2014，20：236239.（in Chinese）） DOI：10.3969/j.issn.10075739.2014.20.139.

[17]張宇，鐘平安，萬新宇，等.近57年江蘇沿海降水量演變特征分析[J].南水北調(diào)與水利科技，2015，13（2）：198201.（ZHANG Y，ZHONG P A，WAN X Y，et al.Evolution characteristics of precipitation along the coastal area of Jiangsu Province in recent 57 years[J].SouthtoNorth Water Transfers and Water Science & Technology，2015，13（2）：198201.（in Chinese）） DOI：10.13476/j.cnki.nsbdqk.2015.02.002.

[18]芮孝芳.水文學原理[M].北京：中國水利水電出版社，2004.（RUI X F.Hydrological principles[M].Beijing：China Water Conservancy and Hydropower Press，2004.（in Chinese））

[19]山城菊，董增川，付曉花，等.近50年來灤河流域徑流演變規(guī)律研究[J].南水北調(diào)與水利科技，2013，11（4）：58.（SHAN C J，DONG Z C，F(xiàn)U X H，et al.Rules of runoff variation in luanhe river basin in recent 50 years[J].SouthtoNorth Water Transfers and Water Science & Technology，2013，11（4）：58.（in Chinese）） DOI：10.3724/SP.J.1201.2013.04005.

[20]邱海軍，曹明明，曾彬.基于小波分析的西安降水時間序列的變化特征[J].中國農(nóng)業(yè)氣象，2011，32（1）：2327.（QIU H J，CAO M M，ZENG B.Wavelet analysis on the temporal series of precipitation in Xi′an[J].Chinese Journal of Agrometeorology，2011，32（1）：2327.（in Chinese）） DOI：10.3963/j.isnn.10006362.2011.01.005.

[21]王哲，只德國，李濤濤，等.海河流域降雨時間序列小波分析[J].海河水利，2012（3）：3032.（WANG Z，ZHI D G，LI T T，et al.Wavelet analysis of Haihe river basin rainfall time series[J].Haihe Water resources，2012（3）：3032.（in Chinese）） DOI：10.3969/j.issn.10047328.2012.03.010.

[22]許月卿，李雙成，蔡運龍.基于小波分析的河北平原降水變化規(guī)律研究[J].中國科學D輯（地球科學），2004，34（12）：11761183.（XU Y Q，LI S C，CAI Y L.The research of rainfall variation based on wavelet in Hebei Plain[J].Science China Ser.D Earth Sciences，2004，34（12）：11761183.（in Chinese）） DOI：10.3321/j.issn：10069267.2004.12.011.第16卷第3期