论文部分内容阅读
近几十年来,随着黄土高原大规模水土保持措施与生态恢复工程的实施,黄河水量沙量锐减,水沙情势巨变。厘清不同时空尺度水沙对变化环境的响应不仅有利于深入理解高度治理下的河流径流输沙动态特征,而且可为水土保持规划及区域生态治理提供科学依据。本研究以黄土丘陵沟壑区西川河流域为例,基于流域逐日降雨及水沙资料,采用多种时间序列统计方法,分析流域不同时间尺度降雨、径流、输沙的演变特征,定量评估降雨及人为因素对流域水沙变化的贡献率;筛选1974-1989年与2007-2013年间112场次洪水事件,基于洪水径流深、洪水历时、洪峰流量三个表征洪水产流输沙特征的指标,综合采用聚类分析和判别分析方法划分不同洪水类型,探究不同洪水类型输沙规律及水沙关系特征。同时,提出了“相似洪水事件法”,量化人类活动对流域径流及输沙减少的贡献率。最后,采用野外人工降雨模拟试验,探究坡面尺度产流产沙对不同植被盖度的响应,量化坡面不同植被组分的减水减沙效应,揭示坡面植被恢复的水沙调控效应。主要研究结论如下:(1)西川河年际尺度径流及输沙量减少的主要驱动因素为人类活动,而降雨变化所造成的影响相对较低。西川河流域多年年降雨量并无显著的增加或减少趋势,但年径流量和输沙量均呈显著的下降趋势(P<0.05);在流域降雨量并未出现显著突变年份情况下,西川河径流量、输沙量时间序列均于1996年发生突变。人类活动对西川河径流量和输沙量减少的贡献率分别为82.0%和87.9%,而降雨变化的贡献率分别为18.0%和12.1%。(2)揭示了西川河径流及输沙量的年内变化特征及其对水土保持措施的响应。流域基准期(1974-1995)与变化期(1996-2012)月均径流量和输沙量的年内分配特征均呈现“单峰”分布,最高值出现于七月。汛期不同月份径流及输沙量在水土保持措施实施后均有明显降低,而汛期径流输沙量在年内的占比并未发生显著变化;随着流域水土保持措施的实施,流域日尺度水沙关系发生了改变;1974-2012年间,年内日最大流量、输沙率以及含沙量均呈显著下降趋势。相比于基准期,变化期历年日最大流量、输沙率以及含沙量的平均值分别下降了37.3%、42.6%、18.1%;1974-1995年间,西川河极端日流量最优分布为Ln 3分布,其5年、10年、25年、50年、100年一遇的年极端日流量分别为63.5 m3·s-1、90.7 m3·s-1、132.1 m3·s-1、168.0 m3·s-1、208.4 m3·s-1。1996-2012年间,极端日流量的最优分布函数为Wakeby分布,其5年、10年、25年、50年、100年一遇的年极端日流量分别为25.1 m3·s-1、28.6 m3·s-1、34.5 m3·s-1、41.5 m3·s-1、52.1 m3·s-1,同一级别重现期日流量水平较上一时段明显降低。(3)基于WASASED模型分析了不同降雨情景下流域产流量的变化特征。采用1988-1991年西川河流域逐日流量资料对WASASED模型进行率定与验证,日流量模拟结果的纳什系数ENS和决定系数R2均超过了0.6,表明模拟结果满足精度要求,可用于流域日尺度径流过程的模拟。保持其余输入因子不变,设置四种降雨情景(降雨量增加5%、降雨增加10%、降雨量减少5%、降雨量减少10%),定量评估不同降雨情景下流域产流量的变化。模型模拟结果表明,相比于实测降雨量条件下流域年均径流深的模拟值,当降雨增量分别为5%和10%时,流域产流量分别增加22.8%和35.2%。当降雨量分别减少5%和10%时,流域产流量分别降低13.8%和31.0%。(4)揭示不同类型洪水的水沙动态特征及洪沙输移效应。综合采用K-medoids聚类分析和判别分析方法将西川河流域112场次洪水划分为4种类型。不同洪水类型产沙量差异较大,A型、B型、C型及D型洪水对西川河输沙总量的贡献率分别为24.5%(1.48×107 t)、11.2%(6.77×106 t)、21.4%(1.29×107 t)和42.9%(2.59×107 t)。其中,A型与C型洪水次洪平均输沙量显著高于其余两类洪水,是流域治理需重点关注的洪水类型。流域主要滞回曲线类型为逆时针滞回曲线与“8”字形滞回曲线,其中,A类型洪水均呈复合型滞回曲线,C类型洪水亦多呈现复合型滞回曲线,B类型洪水以逆时针滞回曲线为主,D类型洪水多为“8”字滞回曲线。此外,相比于基准期,变化期复合型滞回曲线的占比明显降低,而逆时针滞回曲线的占比有所提高。(5)提出“相似洪水事件法”用于量化人类活动对流域径流输沙减少的贡献率。以西川河流域112场次洪水事件为样本,通过阈值设定,共筛选出八组具有相似降雨条件的洪水进行对比分析,结果表明人类活动对西川河径流量减少的贡献率介于20%与79%之间,平均贡献为47.0%,人类活动对输沙量减少的贡献为55%-95%,平均值为81.0%。(6)随坡面植被覆盖度的提升,产流率与产沙率均呈现下降趋势,草被盖度与产流率间呈线性函数关系,而草被盖度与产沙率间呈指数函数关系。不同雨强条件下,随坡面植被盖度提升,初始产流时间均逐渐延长。1.5 mm·min-1雨强下,20%、40%、60%、90%盖度坡面相比于裸土坡面的减水效应分别为24%、36%、53%和79%,减沙效应分别为37%、73%、78%和94%;2.0 mm·min-1雨强下,20%、40%、60%、90%盖度坡面相比于裸土坡面的减水效应分别为17%、28%、37%和71%,减沙效应分别为27%、67%、78%和89%,相同植被盖度的减沙效应高于其减水效应。(7)草被枯落物层对坡面径流的拦蓄作用最强,而根系对坡面减沙的贡献最大。不同雨强条件下,原状草坡的初始产流时间均最长,去除枯落物坡面的初始产流时间次之,裸土坡面初始产流时间最短。1.5 mm·min-1雨强下,枯落物层、植被茎叶以及植被根系的减流效益分别为26%、7%和12%,减沙效益分别为36%、3%和51%;2.0mm·min-1雨强下,枯落物层、植被茎叶以及植被根系的减水效益分别为17%、10%和11%,减沙效益分别为30%、7%和51%。