利用城市雨水回灌地下水防治地面沉降技术方法探讨

利用城市雨水回灌地下水防治地面沉降技术方法探讨 Vo. 26 No. 4l第 26 卷第 4 期 学灾 害 JOURNAL OF CATASTROPHOLOGY 2011 10 Oct. 2011年 月 利用冲击力信号判断泥石流颗粒垂向分选的 * 试验研究 21，2 2 2 ，，，勇杨红娟韦方强胡凯衡洪 ( 1. ，610041;中国科学院山地灾害与地表过程重点实验室四川 成都 ?，610041)中国科学院水利部成都山地灾害与环境研究所四川 成都 2. : ，摘 要分别提出根据冲击力数据的波动强度和峰值情况进行泥石流颗粒垂向分选研究的系统方法并配置了3 3 3 ( 2 095 kg / m、2 008 kg / m960 kg / m) 。1 三组不同容重和 的粘性泥石流样品开展泥石流冲击试验两种方法的分 3 ，2 095 kg / m，析结果基本一致即容重为 的泥石流分选不明显其他两组泥石流呈现出正粒序分布且容重越小分 。，，，选越显著量纲分析表明开展的试验粘滞力在运动中起主导作用颗粒之间作用力较小不容易发生反粒序 ，，分选这与通过冲击力分析颗粒垂向分选的结果一致因此提出的利用冲击力信号判断泥石流颗粒垂向分选的 。系统方法具有适用性 : ; ; ; 关键词泥石流颗粒分选冲击力试验 : P642.2 3: A: 1000 ) 811X(2011)04 ) 0029 ) 06中图分类号文献标志码文章编号 ，泥石流是发生于山坡与沟谷中的含有小至粘野外测量根据距离沟床不同位置处的冲击信号,1,，的尖峰初步分析了颗粒垂向分选状况但是该分 、。粒大至巨石的固液两相流由于重力及颗粒 。，析没有考虑流速的影响为此本文进一步开展 ，之间的相互作用泥石流运动中可能发生固体颗 ，泥石流冲击力试验结合流速资料提出利用冲击 ，粒的垂向分选因此颗粒垂向分选研究是泥石流 。力信号判断泥石流颗粒垂向分选的系统方法 ，动力学研究的重要内容对于认识泥石流的复杂 。运动机理具有重要意义 在泥石流堆积体的不同位置取样并分析样品 1 试验概况的颗粒级配是研究泥石流颗粒垂向分选的常用方 ,2,，Major法如 了美国俄勒冈州森林观测站大 1. 1 ( ) 试验装置和测量系统型泥石流流槽试验堆积物的反粒序上粗下细现 ,3, 本研究利用可变坡钢结构实验槽开展粘性泥 ，象王裕宜等分析了云南蒋家沟某支沟泥石流 。1 ，石流冲击力测试试验装置如图 所示由四部分 。堆积扇的剖面特征堆积体是泥石流停积后形成 ，的因此这种方法得到的颗粒垂向分选信息不能 : 、 、 。构成泥 浆 池实 验 槽传感器支架和回 收 池。完全反映运动时的状况对于不含粘土颗粒的泥 600 cm，30 cm，40 cm，实验槽长 宽 高 实验槽坡( ) ，石流即水石流可通过影像资料观测颗粒在运 。，度可以无级变动实验槽上端连接泥浆池末端 ,4, ，动中的分选情况如王兆印开展的泥石流龙头 ， 。连接回收 池为实验提供物料供应和回收利用。，运动试验研究对于含有粘粒的泥石流流体的 ,5, 2 m ，在距离实验槽末端 位置设立传感器支架供安。Savage 不透明性阻碍了该方法的应用等通过 。装冲击力传感器使用 在实验槽末端设置水平隔板并在回收样品的容器 ，中设置竖直隔板研究干颗粒的垂向分选一些学 ,6 ) 7,者进一步利用该方法研究颗粒与牛顿体混合后 ，的分选状况该方法存在的问题是隔板对流场扰 动较大并且由于大石块的破坏作用难以用于野外 ,8, 。测量胡凯衡等在蒋家沟开展了泥石流冲击力 图 1 试验装置示意 2011 ) 05 ) 23* 收稿日期::( 40771026) ; 基金项目 国家自然科学基金项目中国科学院成都山地灾害与环境研究所青年百人团队项目 :( 1982 )) ，，，，，, E-mai: yanghj@ imde. ac. cnl作者简介 杨红娟 女河南襄城人博士助理研究员主要从事泥石流动力学研究 。，测量系统包括冲击力测量和泥位测量两部分其中一次试验获取的泥深和冲击压强过程线其 、7. 4 cm 冲击力测量系统包括冲击力传感器组数据采集 中冲击压强过程线为 处的传感器获取的数 。8 1 cm 。，6 s，器和计算机传感器组由 个直径 的圆形应 据从泥深过程看流动约持续 期间泥深先 ，，10. 3 cm。变式传感器组成这些传感器被固定在一个中空 增大后减小最大值为 龙头刚到达冲击 ，，7. 4 cm 的支架的迎水面受力面与立柱迎水面位于同一 力测量断面时泥深较小处传感器受力为 。8 1. 4 cm、2.9 cm、，，平面个传感器分别距离槽底 零之后该传感器开始受力在泥深最大时冲击 4. 4 cm、5.9 cm、7.4 cm、8.9 cm、10.4 cm ，， 。和力达到最大随着泥深减 小冲击力逐渐归 零 11. 9 cm。24 ，数据采集器共有 个通道采样频率为这些试验结果为研究泥石流颗粒垂向分选提供了 4 000 Hz。，。采集器与计算机相连冲击力数据可在 试验数据支持 。计算机中同步记录和显示使用激光泥位计测量 ，0. 1 cm，支架处泥深变化过程测量精度为 采样 30 Hz，。频率为 安装在传感器前上方 1. 2 试验和试验结果 试验所用材料来自云南东川蒋家沟泥石流原状 ，1 cm cm ，1 。堆积物经 过筛处理去除大于 的颗粒 2，(, 0. 004 mm、) 试验材料颗粒级配曲线见图 粘粒 ( 0. 004: 0. 062 mm)、 ( 0. 062: 2 mm) 粉砂砂和砾 ( 2 :16 mm ) 8. 5% 、18. 6% 、石所占比例分别为 28. 6% 44. 3% 。和 图 3 某测次试验泥深和 7. 4 cm 处的冲击压强过程 2 颗粒垂向分选分析 2. 1 利用动压强偏差分析颗粒垂向分选 ，以两相流的观点来看泥石流可分为液相和 ，固相两部分其中液相由水和较细的固体颗粒构 ，。，成固相由较粗的固体颗粒构成与此相应泥 图 2 固体颗粒粒径级配曲线 。石流的动压强可分为液相动压强和固相动压强 ，、、试验分三组其容重实验槽坡降重复次,9, ，液相在空间上的分布较为均匀其动压强比较稳1。数和秒表所测龙头平均流速列于表 康志成等，定固相动压强因粒径分布范围较大且空间上不 、、将泥石流分为稀性泥石流亚粘性泥石流粘性 。，连续而不稳定因此相同粒径级配下动压强波 、，泥石流高粘性泥石流等类型其中粘性泥石流 3 ，动强度与固相体积浓度呈正相关关系相同固相 1 950 : 2 300 kg / m。，容重为 据此判断三组样品 体积浓度下动压强波动强度与固相粒径呈正相关 ，III 均为粘性泥石流其中第 组接近粘性泥石流的 。1 容重下限根据表 所列同组各测次龙头平均流速 ,10, 。关系张玉萍 利用三种粒径的颗粒分别配置不，。范围可以看出测量重复性较好 ，同体积浓度的水石流开展冲击试验利用动压力 1 表 泥石流冲击力试验 ，标准偏差表示其波动强度试验结果与该论断相 。，符若泥石流在运动中发生颗粒垂向分选粗颗 容重 / 实验槽龙头平均粒集中的部位与其他部位相比固相体积浓度和粒 分组重复次数3 ( kg / m) 坡降 / %流速 / ( m / s) ，径均会有所增加因此动压强波动强度的垂向分 2 09514.4 2. 33 : 2. 597?。布可以反映颗粒垂向分选状况 2 008 13.1 2. 57 : 2. 69 4 ? 利用动压强波动强度判断颗粒垂向分选状况 ? 1 960 13.1 6 2. 54 : 2. 90 ，时需要考虑信号噪声和流速的影响同样粒径级 ，配和固相体积浓度的泥石流动压强波动强度与 每次试验可以通过测量系统获取泥石流运动 、。噪声强度流速均为正相关关系削减噪声要求 ，，过程中的两组数据即冲击力和泥深数据其中 。3 冲击力除以传感器面积得到冲击压强图 显示了 314 期杨红娟，等: 利用冲击力信号判断泥石流颗粒垂向分选的试验研究 p; ，。小波降噪后的泥石流流速的变化情况记为 对原始冲击压强数据进行降噪处理冲击压强信 k1 ，号的频谱分析表明信号能量集中在低频部分基 冲击压强对应的动压强则反映了流速和颗粒分选 于快速傅里叶变换的低通滤波在使信号变光滑的 ，p。P的综合影响记为 定义相对动压强 为 k2 r ,10, ，同时滤掉了有用信号因此本文参考张玉萍的 p k2 。P= 。( 3)研究利用小波分析进行降噪以某一测次的冲击 r p k1，4 、压强为例图 显示了原始信号低通滤波后的信 : P，，式中为无量纲量消除了流速影响因此可 r 。0 号和小波降噪后的信号图中 时刻前泥石流尚未 P通过 波动强度的垂向分布判断泥石流颗粒垂向 r ，到达测量断面该时段的信号强度可代表噪声强 。分选状况 ，304 Pa，度该时段原始信号的标准偏差为 小波降 60 Pa，， 噪后为 因此小波降噪可以大幅度削弱噪声3 ，图 显示试验中泥深先增加后减小并且在最 。同时保留有用信号 ，1 s 大值处比较稳定因此取平均泥深最大的 时段 ，进行数据分析期间完全淹没在泥石流中的传感 5 。PSTDP，P器有 个记 的标准偏差为 代表 的 r r r ，波动强度分别计算同组各测量高度不同测次泥 1 s STDP，2 深最大 时段内的 的平均值结果如表 r 5 。5 STDP和图 所示第?组 个测量高度 平均值的 r 0. 031，标准偏差为 小于大多数测量高度不同测次 STDP; ，的标准偏差表明颗粒垂向分选不明显第 r II III 0. 063 0. 109，组和第 组则分别为 和 大于大多 STDP，数测量高度不同测次 的标准偏差表明颗粒 r 。，垂向分选较明显从平均值来看第? 组和第? STDP2. 9 cm 4. 4cm组 的最大值分别出现在 和 r ，7. 4 cm ，处最小值均出现在 处因此可以初步判 ( ) 。断为正粒序上细下粗分布 图 4 原始冲击压强信号与降噪处理后的信号 ,11,p通过静水压公式计算静压强 s ( 1)p= g( h )z ) cos，ρθ s c 32 : ( kg / m ) ;= 9. 81 m / s g 式中ρ为泥石流容重为 c ; h ( m) ; z 重力加速度为泥深为传感器中心距离 ( m) ; 。槽底的高度θ为实验槽倾 角冲击压强减去 p，静压强即为动压强 为 k p= p )p ，( 2) k s : p ( Pa) 。式中为冲击压强 低通滤波后的冲击压强对应的动压强反映了STDP5 图 同组各测次 平均值的垂向分布 r STDP2的统计结果 各测量高度 表 r 分组指标1. 4 cm2. 9 cm4. 4 cm5. 9 cm7. 4 cm标准偏差 0. 4830. 5090. 4710. 4390. 433?平均值0. 031 0. 060 0. 048 0. 021 0. 034 0. 040 标准偏差 0. 479 0. 507 0. 453 0. 397 0. 352 平均值 0. 063?0. 097 0. 023 0. 026 0. 031 0. 020 标准偏差 0. 502 0. 595 0. 660 0. 477 0. 379 平均值 ?0. 109标准偏差 0. 108 0. 073 0. 146 0. 046 0. 044 2. 2 、利用动压强尖峰数分析颗粒垂向分选泥石流冲击力概化为锯齿型脉冲矩形脉冲和尖,12, ，章书成等通过对蒋家沟泥石流的观测将 ， ，峰型脉冲其中尖峰型脉冲由大石块撞击引起 因此可通过尖峰型脉冲的数量来判断大颗粒的垂 。向分布情况 ,13,何思明等根据能量平衡原理考察了泥石流 。中大石块横向冲击拦挡坝的问题假设静力条件 ( ) F 下大石块与拦挡坝接触面上的法向应力冲击力 与法向变形量 δ 之间满足如下关系 k F = c， ( 4)δ : c k 。式中和 为材料的特性参数可以推导出泥石 F流大石块最大冲击力 为 max k 2 k + 1 u( k + 1)m 图 6 某测次 2. 9 cm 处泥深最大 1 s 内的 N 过程 c F= c ，( 5)max , , 2c ， 。峰太少 统计结果缺乏代表性 综合考虑两项因 : m ( kg) 。式中为石块质量对于本试验所用冲击' = 5 000。N，首先计算各测次泥深最 素本研究取 c ，k = 1，( 5 ) 力传感器由式可知大颗粒的最大冲击 '1s NN，大 内的 过程并统计大于 的尖峰个数c c p压强 为 max 然 槡 p= ( 6) 1s ( 1 m) ?u，后结合 内的平均流速计算单位距离内的 max cmS ，尖峰个数最后统计出同组各测量高度不同测次 2 : S ( m) 。式中为传感器面积由此可知大颗粒的，7 。1. 4 cm 尖峰数的平均值如图 所示各组 处的 ，最大冲击压强与流速成正比同时与颗粒质量正 ，:尖峰数均明显小于其它位置可能原因有两个 。Npu 。相关记 为 与 的比值 c max ，?该位置距离槽底很近传感器前方因立柱的阻 p max槡 N=。( 7)= ，挡作用而停积有泥石流体对大颗粒撞击具有缓 c cmu S ; ，冲作用?该处剪切速率很高促使大颗粒向上 c S ，N由于各传感器参数 和 相同因此 只与 c 。I ，运动第 组其它四个高度的平均尖峰数接近第 II III ，组和第 组则呈现随高度增加而减小的趋势 ，mN有关某传感器在单位距离内测得的 大于一 c I II ，表明第 组的泥石流体大颗粒分布均匀第 组和，定阈值的尖峰数越多该传感器处大于一定粒径 III ，第 组因重力作用大颗粒出现沉降泥石流容重 ， 的颗粒数 也 越 多据此可以判断大颗粒的空间 ，。I III 越小沉降越明显第 组到第 组泥石流容重 。分布 ，，递减因颗粒级配不变大颗粒含量必然相应减 p，取泥石流的动压强为小波降噪对应的 根 k2 ，，小但是统计结果与此相反这可能是因为容重 ,14,，2 mm 据蒋家沟泥石流研究成果水和 以下固体 ，减小时紊动增强颗粒冲击力因叠加到紊动作用 。，力之上而有所增加 p。颗粒构成液相其动压强 与一般流体相同 kf 2 ( 8)p= 0. 5Cu ，ρ kf f vf 3: ( kg / m ) ; C 式中ρ为液相容重为液相体积浓 f vf ,11, :p。u 度流速 利用低通滤波对应动压强 计算 k1 ( 9)u = 2p/。ρk1 c 槡 p:则固相对应的动压强 为 kc ( 10)p= p)p 。 kc k2 kf ( 7 ) pN，6 根据式和 可以计算出 图 是某测 kc c 1 s N。 次泥深最大 内的 过程只有大于一定阈值c ' ，N的尖峰才对应一个较大的颗粒即需要统计单 c :位长度内满足以下条件的信号个数 ii i ) 1 i i + 1 '' 图 7 同组各测次单位距离内 N, N 的尖峰平均数的 c , ; , ; ( 11)NN Nc N N ,N ， c c c ccc 垂i ) 1 i i+ 1 : 、、i ) 1、i i + 1 NNN式中分别代表 和 时 cc c向分布 N。 刻的 值c 2. 3 泥石流内部作用力的量纲分析 '、泥石流运动阻力包括粘滞力颗粒碰撞力和 ，取值过小尖峰可能是多个小颗粒同时冲 Nc ' ，。，; N( 11 ) 摩擦力等颗粒垂向粒径分布与主导作用力有关击传感器的结果取值过大满足式的尖 c 334 期杨红娟，等: 利用冲击力信号判断泥石流颗粒垂向分选的试验研究 ,15, Iverson 在前人研究基础上了三个无量纲数，值情况进行泥石流颗粒垂向分选研究的系统方法3 3，及相应阈值表征这些力的相对大小包括表征碰 ( 2 095 kg / m、2 008 kg / m并配置了三组不同容重 3Bagnold N、撞力和粘滞力相对大小的 数 表征碰 Bag 960 kg / m ) 1 和 的粘性泥石流样品开展泥石流冲击 。，Savage N、试验两种方法的分析结果基本一致即容重为 撞力和摩擦力相对大小的 数 表征摩擦 Sav 3N， 力和粘 滞 力相对大小的摩擦数 其 表 达 式 分 2 095 kg / m ，f 的泥石流分选不明显其他两组泥石 :别为 。流呈现出正粒序分布且容重越小分选越显著量 2 ，纲分析表明本文开展的试验粘滞力在运动中起 1 ) C D ργ vf s =，( 12)N? Bag ，，主导作用颗粒之间作用力较小不容易发生反 C μ vff，粒序分选这与通过冲击力分析颗粒垂向分选的 2 2 D ργ s N= ( 13)，，结果一致因此本文提出的利用冲击力信号判断 Sav ( ρ) ) ghtanρ s f。泥石流颗粒垂向分选的系统方法具有适用性该 1 )C ( ) ) ghtanρρ vf s f ( 14)，N=? f C μ γ vff ，，方法会受紊动影响因此只能在同一流态下进行 3: g / m ; D ，2 700 k式中ρ为颗粒密度本文取 为 s 。不同流态之间的计算结果不具有可比性 ，5. 7 mm;固相中值粒径根据粒径级配曲线为 γ ，u为剪切速率通过龙头平均流速 和泥深计算 surge :参考文献 ( = u/ h) ; ， 得到γ 为颗粒内摩擦角根据康志surge ,9, ,1,费祥俊，舒安平, 泥石流运动机理与灾害防治,M,, 北京: 清29. 8?;， 成等的实验结果取 μ为浆体粘滞系 数 f2004: 12,，华大学出版社 。通过液相流变试验得到 ,2, Major J J, Experimenta studies of deposition at a debris fow fumelll利用分析颗粒垂向分选时段的平均泥深计算 ,R,, U S Geological Survey Fact Sheet，19942: 8 , Bagnold 、Savage 得到三组试验各测次的 数数和摩 ,3,王裕宜，詹钱登，韩文亮，等, 泥石流堆积层理结构的分析研，3。I III ，擦数其范围列于表 从第 组到第 组随着 ,J,, ，2001，15( 3) : 68 ) 71,究水土保持学报 , J,, ,王兆印泥石流龙头运动的实验研究及能量理论水利学 ，，泥石流容重减小液相粘滞系数显著减小剪切 ,4,2001( 3) : 18 ) 26, ，报，NNIII 。速率有所增加故 和 逐渐增加第 组Bag f Savage SB ，Lun C K K, Particle size segregation in inclined chute ,5,flow of dry cohesionless granular solids ,J,, Journal of Fluid Me- 因 ,15,chanics，1988，189: 311 ) 335, NIverson ，。剪切速率较大大于前两组根据 界 Sav Vallance J W，Savage SB , Particle segregation in granular flows ，N, 200 ， 定的阈值表明碰撞作用超过粘滞作用Bag down chutes ,C,/ / RosatoA ，Blackmore D, Segregation in granu- ,6,lar flows ( International Union of Theoretical and Applied Mechan- N, 0. 1 ，N,2 表明碰撞作用超过摩擦作用Sav f ics Symposium) ,Dor d recht: K luwer，2005: 31 ) 51, 000 ，表明摩擦作用大于粘滞作用由此判断三组Zanuttigh B，Ghilardi P, Segregation proces s of water-granular mixtures released down a steephu cte ,J,, Journal of Hydrology， ，试验 中粘滞力在泥石流运动中均起主导作用 ,7,2010，391: 175) 187, ,16, 。其次为 颗粒碰撞力目前主要通过运动筛分和，，，, ,J,,胡凯衡韦方强洪勇等泥石流冲击力的野外测量 25( s1) : 2813) 2819,，岩石力学与工程学报 颗粒碰 ，，，, ,M,, :,8,康志成李焯芬马蔼乃等中国泥石流研究北京 ,17, ，撞理论解释颗粒物质运动中反粒序的形成本，2004: 21 ) 96,科学出版社 研究开展的试验中颗粒作用力在泥石流运动中不 , ,D,, : 张玉萍泥石流冲击信号识别方法研究重庆重庆交 ,9,2009, ，通大学，。， 起主导作用故没有出现反粒序分布与此同时Yang H，Wei F，Hu K，et al, Measuring the internal velocity of ，，泥石流容重越大液相粘滞系数越小颗粒沉降 ,10,debris flows using impact pressure edtecting in the flume experi- ，。速度越大因此正粒序分布越显著 ment, J,, Journal of Mountain Science，2011，8( 2) : 1 ) 9,章 3无量纲数计算结果表 ,11,, ,C,/ / 泥石流冲击力及其测试中国科学院 兰，书成袁建模, 中国科学院兰州冰川冻土研究所集刊 第 剪切速州冰川冻土研究所液相粘滞系N NN分组 BagSavf) 14 ( ) , : 1985:， 号中国泥石流研究专辑北京科学出版社 / ( Pas)?率 / s 数 269 ) 274, ,12, 24. 0 :0. 29 :2. 6 :0. 11 :22. 6 : ?, J,,，，,何思明吴永沈均泥石流大块石冲击力的简化计算 26. 50. 333. 00. 1325. 72009，18( 5) : 51 ) 56,，自然灾害学报 25. 6 : 0. 13 : 5. 4 : 0. 11 : 43. 2 : , M,,，，,?王裕宜詹钱登严璧玉泥石流体结构和流变特性 27. 30. 156. 10. 1354. 1: ，2001: 5 ) 7,长沙湖南科学技术出版社 26. 8 : 0. 07 : 10. 4 : 0. 12 : 58. 6 : ,13,Iverson R M, The physics of debris flows ,J,, Reviews of Geo- ? 33. 80. 0812. 80. 2287. 9 physics，1997，35( 3) : 245) 296, ,14,Mideleton G V, Experimental studies related to problems of flysch 3 结论 ,15, ,16,本文分别提出根据冲击力数据的波动强度和峰 sedimentation ,C,/ / Lajoie J, Flysch Sedimentology in NorthBagnold R A, Experiments on a rgavity-free dispersion of large,17,America, Toronto: B usiness and Economics Service，1970: 253 solid spheresin a Newtonian fluid under shea,r J,, Proc Roy ) 272, Soc，1954，22549: ) 63, Expeimental Study on Vetical Soting of Paticles in Debis rrrrr Flow with Impact Signals 1，2 2 2 2Yang Hongjuan，Wei Fangqiang，Hu Kaihengand HongY ong ( 1, Key Laboratory of Mountain Hazards and Surface Process，Chinese Academy of Sciences， Chengdu 610041，China; 2, Institute of Mountain Hazards and Environment，Chinese Academy of Sciences， Chengdu 610041， China)Abstract: Studying vertical sorting of p articles is important for understanding the complex mechanism of debris flow, In this paper wer opposed two methods jtuod ge vertical sorting of particles in debris flow, Onei s based on the standard deviation of impact signals，and the otheris based on the peaksin im pact signals, We carried out 3 3 3 three groups oflu mfe tests with debris densities of 2 095 kg / m，2 008 kg / mand 1 960 kg / m，respectively,Results of the two methods sareimil ar, There is no significant vertical sorting with the first group of edbris flows,There is normal sorting with the other two groups，andn ortmhea l sorting becomes moresig nificant for debris flows wth smaer densty, Dmensona anayss ndcates that vscous stress s domnant in debrs fow movement, s illiiilliiiiiiilIt idifficult for thes maller stress between aprticles to make a reverssoert ing, This is consistent with the results using mpact sgnas, Therefore，the methodsropo psed in ths paper areapp cabe to udge vertca sortng of partces in iilililjiliildebris flow, Key words: debris flow; particle sorting; impact pressure; te st檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪 ( 28 )上接第 页 Spatial and Tempoal Evolution Chaacteistics of Exteme rrrr Rainfalls in the Pearl River Basin 1，2 1，2 1，2 3PengJ unta，Zhang Qang，Chen aohongand Chen ongqn iiXiYi ( 1, Sun Yat-sen University，Department of Water Resources and Environment，Guangzhou 510275， China; 2, Sun Yat-sen University，Guangdong University Key Laboratory of Water Cycle and Security in South China，Guangzhou 510275，China; 3, The Chinese University of Hong Kong，Department of Geography and Resource Management，Shatiian : Three extremera infall indexes as the m aximum consecutive rainy days M( CP) ，the Abstract maxmum i China) number ofco nsecutive rainy days ( C) and related maximum precipitation intensity ( CI) and theZ index MDMare defined and analyzed based on thea ildy precipitation data of42 rain gauging stations in the period of 1960: 2005, Pearl River basin, The annual and summerM CP，MCD，MCI and Z values are high in the Pearl River The resuts ndcate that: ( 1) h gher / ower ，，I andZ vaues are foundin the east / west parts liiilMCPMCDMCldelta，Gun，Base andC hongzuo，mpyng hgher food rsk in theser egons, No evdent changes of ，，iliiiliiliiiMCPMCD of the MCI and Z values can be observedin winter; ( 2) Nos ignificant abruptb ehaviors can bei dentified in thesev ariables;( 3) Thes patial distribution of precipitation extremesis in close relations with that ofto pography, Those areasi twh hgher probabty of precptaton extremes arues uay surrounded twh mountans and mayea sy trgger mountanous iiliiiilliiilii flash floods, In this sense，there sults of this study should be of greatsc ientific and practical merits in scientific water resources management and enhancement mofit ighuatimonan to natural hazards, Key words: precipitation extremes; M ann-Kendall analysis; spatio-temporal changes; thPeea rl River basin