南水北调工程是解决我国北方地区水资源短缺,实现社会经济可持续发展的特大型跨流域调水工程,分别从长江上、中、下游调水,经东、中、西三线输送至水资源短缺的黄淮海流域。南水北调东线工程利用江苏省已建的江水北调工程,逐步扩大调水规模并延长输水线路,工程从长江下游扬州江都水利枢纽抽引长江水,利用京杭大运河及与其平行的河道逐级提水北送至黄淮海平原东部和胶东地区,其中江苏境内受水区涉及扬州、淮安、盐城、宿迁、连云港和徐州六市,总面积3.81万km2。

随着全球区域气候的变化,将导致区域降雨量发生变化,从而对水资源供给和自然生态系统产生重大影响。因此,分析研究区域降雨变化的趋势和周期变化规律[1-4],对认识区域降雨特征、应对强降雨等极端气候以及农业生产、防汛抗旱、水资源规划等具有重要意义。南水北调工程实施以来,水文工作者进行了很多有价值的研究,但是有关降雨特征分析的研究仅南水北调中线京石段开展过[5]。目前研究降雨变化趋势常用的方法有滑动平均法[6-7]、线性估计、Mann-Kendall秩次相关检验法[8-9]等,常用的周期变化研究方法有回归周期分析法、傅立叶分析法[10]、小波分析法等。本文以南水北调东线江苏境内受水区的沂南水利分区作为研究区域,利用滑动平均法、Mann-Kendall秩次相关检验法、游程检验法[11-12]、Spearman秩次相关检验法[13-15]四种方法相结合对该地区降雨变化趋势进行分析,并应用Morlet小波变换[11-12]分析降雨的时间及周期变化特征,此研究为南水北调东线江苏境内水量调度运行管理、区域抗旱排涝调度工作的开展以及水资源开发与利用等提供参考。

1 区域概况

江苏省沂南水利分区(简称:沂南地区)是南水北调东线江苏境内水量配置模型的计算单元之一,位于新沂河以南、废黄河以北,东临黄海、西至中运河,涉及四市七县区,包括盐城市的响水县,连云港市的灌南县,淮安市的涟水县、淮阴区,宿迁市的宿豫区、泗阳县、沭阳县,区域总面积7 477 km2。区域地形总趋势是西高东低,自西南向东北倾斜,东西向平均地面比降约万分之一点三。该区域地处亚热带向暖温带过渡性气候中,又是南北气流过渡地带,有明显的海洋季风特征,全年四季分明,夏季暖热多雨,受台风侵袭频繁,冬季寒冷干燥,春秋两季昼暖夜寒,冷空气活动频繁,温差变化较大。

2 暴雨类型及特征

区域降雨相对比较集中,多突发性暴雨,最大日降雨量610.3 mm(陈港站,2000年8月30日)。降雨年际变化幅度大,多年平均面雨量为926.12 mm(1951-2012年),最大年降雨量为1 489.69 mm(2003年),最小年降雨量仅为 572.24 mm(2004年),多年平均年降雨天数99.2 d,最多年降雨天数为131 d(1954年),最少年降雨天数63 d(1957年);降雨年内分配极不均匀,汛期6月-9月降雨量比较集中,约占全年降雨量的70%,秋末的10月至次年的2月,降雨量一般只占全年降雨量的15%左右。

台风则多发生在7月到9月,尤其以8月最频繁,常因8～9级大风大潮伴随着100 mm以上降雨同时发生,极易产生突发性的洪、涝、潮等自然灾害。另龙卷风、六级以上的大风、冰雹和5月中下旬的干热风,在局部地区也会造成危害[16-17]。

3 资料选取及分析方法

3.1 资料选取

根据沂南地区雨量站网分布,并综合考虑降雨资料的代表性、可靠性以及完整性,选用区域内的宿迁闸、泗阳闸、淮阴闸等37个雨量站点1951-2012年逐日降雨数据,利用泰森多边形法[18]进行面平均雨量计算,以面平均雨量作为区域降雨趋势分析的资料,能更好更全面地分析该地区的降雨变化特征。沂南地区雨量站点分布见图1。

  

图1 沂南地区雨量站点分布Fig.1 Location of precipitation stations in Yinan Area

3.2 分析方法

式中:w0为常数,i为虚数。通过小波变换后,得到的小波变换图能够反映降雨在不同时间尺度上的周期变化[19]。小波系数实部为正的用实线表示,代表降雨偏多;小波系数实部为负的用虚线表示,代表降雨偏少。小波方差图反映了能量随时间尺度的分布,可以用来分析一个时间序列中各种尺度扰动的相对强度。对应峰值处的尺度称为该序列的主要时间尺度,即主要周期[20-22]。为了进一步确定降雨的周期变化,需通过小波方差检验,确定显著周期。小波方差的计算公式为:

(1)滑动平均法。

在一个系列{xn}中,将连续的2n或2n+1个数值取平均,然后向后逐年滑动求均值直到最后一个数据,得到的新序列{yn}比原序列{xn}更光滑,即滑动平均法。

(2)Mann-Kendall检验法。

本文对沂南地区1951-2012年面平均雨量采用5年滑动平均法,直观反映该地区降雨的变化趋势,同时采用Mann-Kendall秩次相关检验法、游程检验法、Spearman秩次相关检验法,选用置信水平α=0.05,通过比较检验统计量与临界值大小来判断降雨变化趋势的显著性。沂南地区年降雨量5年滑动平均过程线见图2,显著性检验结果见表1。

应用 Mann-Kendall方法检验降雨时间序列{xn}的变化趋势时,构造统计量S如下:

(4)Spearman秩次相关检验法。

 

式中:

 

统计量 S服从正态分布,均值为 0,方差Var(S)=n(n-1)(2n+5)/18。当n＞10时,将统计量S标准化:

 

T服从自由度为n-2的t分布,原假设为序列无趋势,则根据降雨时序的秩次相关系数计算统计量T,然后选择显著水平α,在t分布表中查出临界值tα/2,当|T|＞ tα/2时,则拒绝原假设,说明序列随时间有相依关系,从而推断序列趋势明显,否则,接受原假设,趋势不显著[13-15]。

从包含有n个数据点的数据集中随机选取K( )K=3个数据点作为初始聚类中心点，然后采用欧几里得公式计算出每个数据点到K个聚类中心点的距离，并将数据点归类到离其最近的聚类中心点所在的类当中，当所有的数据点归类完成以后，分别计算出每个类中数据点的均值，将类均值作为新的类中心，并重新进行聚类，如此不断循环，直到前后两次聚类中心点相同，则聚类结束。

(3)游程检验法。

根据已知的降雨时间序列{xn},求出相应的距平系列若将连续的0或0看作一个游程,统计序列中出现的游程总数S,其中n1为序列中x'i＞0出现的总数,n2为序列中x'i＜0出现的总数。实际分析中,如n1或n2大于15时,构造统计量U:

武陵山特色旅游小镇应在武陵山场镇的基础上优化建设，要走小而精的建设道路，避免大规模大拆大建，遵循紧凑布局和节约建设用地的原则，完善现有建设区，尽可能减少对自然生态环境的破坏。同时以此为契机，推进新型城镇化，促进城乡一体化发展。

 

构造统计量U近似服从N(0,1)标准正态分布,式中:

5）我们在现场施工时，可以运用信息化模型和现场情况进行对比，可以更加直接准确的检查施工情况和修缮保护情况。

 

根据游程检验的S分布表,选取某一显著α水平作双边检验。当|U|＞Uα/2时,拒绝原假设,序列趋势变化明显,反之不明显[11-12]。

急性阑尾炎在临床上是一种比较常见的急腹症，该疾病在各个年龄段都会出现。目前，在临床上主要的治疗方式就是手术治疗，患者在经手术治疗之后，虽然可以起到一定的治疗效果，但是术后会严重影响患者的机体，给患者带来一定的疼痛感。因此，经过不断的研究与发现，在急性阑尾炎术后对患者实施有效的护理服务，可以起到良好的效果。因此，在本次研究中，主要是对60例急性阑尾炎患者的护理方式进行深入性的研究，并进行了全面性的总结。

小波分析是一种信号的时间/尺度(频率)分析方法,具有多分辨率分析的特点,在时频两域都具有表征信号局部特征的能力,被广泛应用于水文时间序列的尺度分析、水文预测预报及水文随机模拟。Morlet小波是一维连续的,其波形与降雨的波形较为接近。函数表达式为:

 

式中:n为序列长度;di=Ri-i。如果秩次与时序相近,则较小,秩次相关系数较大,趋势性显著。通常采用t检验法检验降雨时序的趋势性是否显著,统计量T的计算公式为:

2012年8月2—4日，福建省柘荣县受第9号台风“苏拉”影响，遭遇持续强降雨，雨量达328 mm，最大1小时雨量达88 mm、最大3小时雨量达141 mm。特大暴雨造成柘荣县境内河水急剧上涨，城区局部内涝，各乡镇不同程度受灾，多处护岸、堤防和道路被冲毁，山体滑坡频发。面对特大强降雨，柘荣县及时启动县级山洪灾害防御预案。各地各部门按照预案加强工程巡查、及时转移群众，23支应急分队果断处置险情，各地防御预警铜锣长紧急鸣锣告警。8月3日23时，乍洋乡石山站监测到最大1小时降雨量达88 mm，由于山洪灾害监测预警系统及时预警，县乡两级立即启动应急预案，转移下游村民300多人，避免了人员伤亡。

 

若当U＜0,则认为有下降趋势,反之呈上升趋势。在显著水平α下,若|U|≥Uα/2,认为趋势性显著,反之趋势性不显著[8-9]。

煤矿井下断层常以大中断层和小断层为主，采掘生产中常见到的大中断层的落差一般在5 m以上，这种断层在地质勘探阶段生产准备阶段就能查明，常将其设置为采掘区或工作面的自然边界。但是小于5 m的小断层则不宜查明，若前期工作不十分仔细，会给后期工作带来不便。须在进行采区和工作面设计之初，调查清楚各断层的展布及延伸状况，避免“二次工作”支出的额外费用。

(5)Morlet小波变换。

Spearman秩次相关检验主要通过分析降雨序列{xi}与其时序i的相关性而检验降雨时序是否有趋势性。在运算时,将降雨序列按从大到小排列,用其新秩序Ri代表,计算秩次相关系数:

 

本文采用滑动平均法、Mann-Kendall秩次相关检验法、游程检验法、Spearman秩次相关检验法等4种方法来研究分析区域降雨趋势。同时应用Morlet小波变换分析该地区降雨的时间及周期变化特征。

2.2.1 开敞式景观配置模式。垂直面上层是高大乔木，下层由矮生灌木、大面积草坪或地被组成。槭树作为前景树或中层树种，小面积植物群落种植或孤植于整个空间的边缘或中心，形成外向、开敞式空间。该空间视野开阔，视线通透，多用于道路景观节点、林缘、广场等景观配置中。典型配置模式：①乔灌草型，白蜡(女贞等)+三角枫(五角枫等)—月季(牡丹等)—草坪； ②乔灌型,建始槭(三角枫等)—麻叶绣线菊；③乔草型，茶条槭(飞蛾槭、红枫、复叶槭等)—草坪(或麦冬、玉簪)。

 

4 降雨变化特征分析

4.1 降雨变化趋势分析

这则有些像笑话的故事，其实不是笑话，而是明人沈德符《万历野获》第二十二卷“海忠介抚江南”一章中，讲述的海瑞巡抚江南之时的一件“轶事”（海瑞谥号“忠介”，故史上又称他为“海忠介”）。这故事的发生有这样一个背景，那就是海瑞对于下级官员趋步至官厅，按礼谒见长官习惯性地弓腰低眉很有些看不惯，因此他下令郡县的长官“庭参”之时不得俯首；但是，他命令虽然下达了，可鉴于作为属下的那些官吏一个个十分畏惧上司，所以，依然没有人敢于仰视上级。

图2可以直观反映沂南地区年降雨量的变化情况,降雨呈现增多和减少交替变化,但总体趋势是减少了,同时通过Mann-Kendall检验,U=-0.996＜0,也定量反映出降雨呈下降趋势。并且从表1可以看出,三种检验一致表明这种下降趋势是显著的。

  

图2 沂南地区年降雨量5年滑动平均过程线Fig.2 T he 5-year moving average curve of rainfall in Yinan area

 

表1 沂南地区年降雨量变化趋势检验T ab.1 T he change trend test of annual rainfall in Yinan area

  

方法 检验统计量 显著水平 临界值 结果 趋势性Mann-Kendall秩次相关检验法 U=-0.996 α=0.05 1.96 |U|＞Uα/2 显著游程检验法 U=0.39 α=0.05 1.96 |U|＞Uα/2 显著spearman秩次相关检验法 T=0.942 α=0.05 2.00 |T|＜Tα/2 显著

4.2 降雨周期变化分析

应用Matlab的Morlet小波分析工具对沂南地区1951-2012年降雨资料进行周期特性分析,得到小波系数实部等值线图和小波方差图,具体见图3、图4。

  

图3 小波系数实部等值线图Fig.3 T he isograms of the real part of wavelet coefficient

  

图4 小波方差图Fig.4 T he variance of Morlet wavelet coefficient

小波系数实部等值线图能反映降雨系列不同时间尺度的周期变化及其在时间域中的分布,从而判断在不同时间尺度降雨的未来变化趋势。图3中,实线表示小波系数正值部分,代表降雨偏多,H表示正值中心;虚线表示小波系数负值部分,代表降雨偏少,L表示负值中心。从图中可以清楚的看到降雨演变过程中存在多时间尺度特征,大致为22～32 a,12～21 a,3～11 a的3类尺度的周期变化规律。其中,22～32 a、12～21 a两个时间尺度上,周期变化在整个分析时间中表现的比较稳定,具有全域性。在22～32 a尺度上先增后减,出现多-少准3次震荡;在12～21 a时间尺度上,出现多-少6次震荡。在3～11 a时间尺度上1954-1981年较为明显,多—少循环交替。

从图4小波方差图可以看出,不同时间尺度扰动的相对强度不同,图中存在4个较为明显的峰值,依次 27 a、15 a、10 a、4 a。其中 27 a对应着最大峰值,说明27 a左右的周期震荡最强,为沂南地区降雨变化的第一主要周期;15 a对应着第二峰值,为沂南地区降雨变化的第二主要周期;10 a对应着第三峰值,为沂南地区降雨变化的第三主要周期;4 a对应着第四峰值,为沂南地区降雨变化的第四主要周期。以上4个周期的波动控制着降雨在整个时间域内的变化特征。

根据小波方差检验结果,绘制控制区域降雨演变的第一、二周期即27 a、15 a特征时间尺度的小波系数图,见图5、图6。从主周期趋势图中可以分析出在不同的时间尺度下,区域降雨存在的平均周期及多-少变化特征不同。图5表明,在27 a时间尺度上,平均周期为18 a左右,大约经历了3个多-少转换期。降雨峰值出现在1955年、1973年、1990年和2009年。图6表明,在15 a时间尺度上,平均周期为10 a左右,大约经历了6个多-少转换期。降雨峰值出现在 1953年、1963年、1973年、1982年、1991年、2000年、2011年。

  

图5 沂南地区降雨变化的27 a特征时间尺度小波实部过程线Fig.5 The curve of the real part of wavelet coefficient for rainfall change on 27 a characteristic time scale in Yinan area

  

图6 沂南地区降雨变化的15 a特征时间尺度小波实部过程线Fig.6 The curve of the real part of wavelet coefficient for rainfall change on 15 a characteristic time scale in Yinan area

5 结论

(1)在降雨序列具有显著变化趋势时,滑动平均法所得到的曲线能够直观的反映出变化趋势,但是不能表明趋势是否显著。文中采用3种检验方法与滑动平均法相结合,既可以直观反映变化趋势,又可以定量分析检验趋势的显著性,同时三种检验方法可以对比参考。最终分析结果一致表明沂南地区1951-2012年面平均降雨量呈现增多和减少交替变化,但总体呈递减趋势,并且趋势显著。

砖子不吭声，他不知如何吭声，他从反馈的消息得知，这次赵仙童确实没有李金枝演得好，那是一锤定音的事，谁也改变不了的，除非下次再去摘。他也明白，摘梅花是多少戏人的终极追求，赵仙童例外不了。可命运女神是不会光顾每一个求梦者的，这就是现实，残酷的现实啊。

Fig.2 shows a schematic of the fiow duct test rig in the Fluid and Acoustic Engineering Laboratory(FAEL)at BUAA.Theflow is provided by a continuous wind tunnel consisting of a centrifugal fan,plenum with metal screens,and silencing ducts.

(2)通过Matlab的Morlet小波分析,沂南地区降雨存在22～32 a,12～21 a,3～11 a的3类尺度的周期变化规律。不同的时间尺度表现出不同的循环交替。小尺度周期震荡剧烈,大尺度周期震荡平稳,规律更明显,全局性更强。

(3)通过小波方差判断降雨时间序列各种时间尺度扰动的相对强度,沂南地区降雨受27 a、15 a、10 a、4 a四个时间尺度的干扰较大,其中27 a干扰强度最大,为该序列的主要时间尺度,即主要周期。在27 a时间尺度上,该地区降雨的平均周期为18 a左右,大约经历了3个多-少转换期。

参考文献(References):

[1] 潘杰.常州地区1978年-2012年降雨量时空分布特征[J].南水北调与水利科技,2015,13(2):245-248.(PAN J.Spatial and temporal characteristics of precipitation in Changzhou Area from 1978 to 2012[J].South-to-North Water Transfers and Water Science&T echnology,2015,13(2):245-248.(in Chinese))DOI:10.13476/j.cnki.nsbdqk.2015.02.012.

[2] 徐宗学,张楠.黄河流域近50年降水变化趋势分析[J].地理研究,2006,25(1):27-34.(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):27-34.(in Chinese))DOI:10.3321/j.issn:1000-0585.2006.01.004.

[3] EVANS J P,JAKEM AN A J.Development of a simple,Catchmentscale,rainfall-evaportranspiration-runoff model[J].Environmental Modelling&Software.1998,13:385-393.

[4] WAGENER T,LEES M J.Rainfall-Runoff Modeling T oolbox user mannual[J].Department of Civil and Environmental Eegineering,Imperial College London,U K.2001.

[5] 尹航,王文君,吴英杰,等.南水北调中线京石段降雨序列特性分析[J].水资源与水工程学报,2016,27(4):113-118.(YIN H,WANG W J,WU Y J,et al.Analysis of rainfall characteristics in Beijing-Shijiazhuang section of South-to-North Water Diversion[J].Journal of Water Resources&Water Engineering,2016,27(4):113-118.(in Chinese))DOI:10.11705/j.issn.1672-643X.2016.04.21.

[6] 裴益轩,郭民.滑动平均法的基本原理及应用[J].火炮发射与控制学报,2001(1):21-23.(PEI Y X,GUO M.The fundamental principle and application of sliding average method[J].GUN Launch& Control Journal,2001(1):21-23.(in Chinese))DOI:10.3969/j.issn.1673-6524.2001.01.007.

[7] 李娟,张维江,马轶.滑动平均-马尔可夫模型在降雨预测中的应用[J].水土保持研究,2005,12(6):196-198.(LI J,ZHANG W J,M A Y.Application of moving Average-Markov Chain in the prediction of precipitation[J].Research of Soil and Water Conservation,2005,12(6):196-198.(in Chinese))

[8] 康淑媛,张勃,柳景峰,等.基于Mann-Kendall法的张掖市降雨量时空分布规律分析[J].资源科学,2009,31(3):501-508.(KANG S Y,ZHANG B,LIU J F,et al.Analysis of the spatiotemporal distribution of precipitation in Zhangye City Using Mann-Kendall method[J].Resources Science,2009,31(3):501-508.(in Chinese))DOI:10.3321/j.issn:1007-7588.2009.03.022.

[9] 郭小芹,刘明春,钱莉,等.从Mann-Kendall特征看石羊河流域降雨量的演变规律[J].干旱区地理,2010,33(4):593-599.(GU O X Q,LIU M C,QIAN L,et al.Evolution regularity of precipitation in the Shiyang River Basin on Mann-Kendall features[J].Arid Land Geography,2010,33(4):593-599.(in Chinese))DOI:10.13826/j.cnki.cn65-1103/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].北京:化学工业出版社,2005.(WANG W S,DING J,LI Y Q.Hydrology wavelet analysis[M].Beijing:Chemical Industry Press,2005.(in Chinese))

[12] 丁晶,刘权授.随机水文学[M].北京:中国水利水电出版社,2008.(DING J,LIU Q S.Stochastic hydrology[M].Beijing:China Water Power Press,2008.(in Chinese))

[13] 王国庆,贾西安,陈江南,等.人类活动对水文序列的显著影响干扰点分析——以黄河中游无定河流域为例[J].西北水资源与水工程,2001,12(3):13-15.(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):13-15.(in Chinese))DOI:10.3969/j.isnn.1672-643X.2001.03.004.

[14] 曾小凡,翟建青,姜彤,等.长江流域年降雨量的空间特征和演变规律分析[J].河海大学学报(自然科学版),2008,36(6):727-732.(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):727-732.(in Chinese))DOI:10.3876/j.issn.1000-1980.2008.06.001.

[15] 杜富慧,郝振纯,陈新美,等.邯郸地区年降雨量统计特性[J].水资源保护,2009,25(6):16-20.(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):16-20.(in Chinese))DOI:10.3969/j.issn.1004-6933.2009.06.005.

[16] 徐晨光,冯新峰,郑志伟,等.江苏省降雨时空分布特征研究[J].现代农业科技,2014,20:236-239.(XU C G,FENG X F,ZHENG Z W,et al.Study on characteristics for spatial and temporal distribution in Jiangsu province[J].M odern Agricultural Sciences and T echnology,,2014,20:236-239.(in Chinese))DOI:10.3969/j.issn.1007-5739.2014.20.139.

[17] 张宇,钟平安,万新宇,等.近57年江苏沿海降水量演变特征分析[J].南水北调与水利科技,2015,13(2):198-201.(ZHA NG Y,ZH ONG P A,WAN X Y,et al.Evolution characteristics of precipitation along the coastal area of Jiangsu Province in recent 57 years[J].South-to-N orth Water T ransfers and Water Science& T echnology,2015,13(2):198-201.(in Chinese))DOI:10.13476/j.cnki.nsbdqk.2015.02.002.

[18] 芮孝芳.水文学原理[M].北京:中国水利水电出版社,2004.(RU I X F.Hydrological principles[M].Beijing:China Water Conservancy and Hydropower Press,2004.(in Chinese))

[19] 山城菊,董增川,付晓花,等.近50年来滦河流域径流演变规律研究[J].南水北调与水利科技,2013,11(4):5-8.(SHAN C J,DONG Z C,FU X H,et al.Rules of runoff variation in luanhe river basin in recent 50 years[J].South-to-North Water Transfers and Water Science& Technology,2013,11(4):5-8.(in Chinese))DOI:10.3724/SP.J.1201.2013.04005.

[20] 邱海军,曹明明,曾彬.基于小波分析的西安降水时间序列的变化特征[J].中国农业气象,2011,32(1):23-27.(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):23-27.(in Chinese))DOI:10.3963/j.isnn.1000-6362.2011.01.005.

[21] 王哲,只德国,李涛涛,等.海河流域降雨时间序列小波分析[J].海河水利,2012(3):30-32.(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):30-32.(in Chinese))DOI:10.3969/j.issn.1004-7328.2012.03.010.

[22] 许月卿,李双成,蔡运龙.基于小波分析的河北平原降水变化规律研究[J].中国科学D辑(地球科学),2004,34(12):1176-1183.(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):1176-1183.(in Chinese))DOI:10.3321/j.issn:1006-9267.2004.12.011.