黄河上游近60年来水沙变化及其影响因子.docxVIP

黄河上游近60年来水沙变化及其影响因子 黄河上游（龙羊峡以上）面积约占黄河总面积的15.4%，多年平均流量仅占总流量的38.5%。平均流量模式约为总流模式的2.5倍。排水沟的丰和旱季变化不仅影响黄河上游生态环境的发展，而且直接影响黄河下游资源的变化。因此, 深入研究黄河上游区水资源的变化机理是源区生态保护和全流域经济可持续发展的关键。 在青藏高原大江河径流量演变研究方面, 施雅风、赖祖铭等认为, 20世纪正处于气候暖干化、水资源萎缩过程中, 江河流域的年径流量自60年代中期以来呈现下降趋势, 黄河唐乃亥以上流域如果气温升高2℃, 即便降水量增加20%, 该流域的夏季径流量也将减少3%;黄河流域天然径流量是由周期性和随机性成分组成的, 夏秋季降水量对同期及落后半年之内的径流有显著贡献;黄河兰州东乡以上4个流域径流有减少的趋势, 但程度各不相同, 90年代径流年内分配特征出现了较大的变化, 突出表现在汛期径流量的减少。 黄河上游地区20世纪90年代以来的暖干化趋势, 使得三江源区冰川退缩、冻土退化、湖泊萎缩、草场退化、土地荒漠化、水土流失、生物多样性遭受破坏[9,10,11,12,13,14,15,16,17], 而这些生态环境问题都与水循环和水资源有关。利用该地的水文、气象资料, 着重分析不同季节径流量的变化规律及其与影响因素的关系, 以期为黄河上游地区的经济建设和生态保护提供依据。 1 数据和方法 1.1 逐月流量资料 利用黄河上游地区(32°—36° N, 96°—103° E) 11个气象台站 (兴海、同德、河南、泽库、玛沁、甘德、达日、久治、玛多、玛曲、若尔盖) 1956—2005年1—12月逐月降水、气温和3个水文站 (达日、玛曲、唐乃亥) 同期的逐月流量资料 (资料由青海省气候资料中心和水文水资源勘测局提供) 。 1.2 计算方法 1 黄河上游地区时段降水量、加权系数分布 Ρ=W1×Ρ1+W2×Ρ2+W3×Ρ3(1)P=W1×P1+W2×P2+W3×P3(1) 式中:P1、P2、P3和W1、W2、W3分别为黄河上游地区1、2、3区的时段平均降水量和平均加权系数。黄河源头至达日段为1区, 达日以下至玛曲段为2区, 玛曲以下至唐乃亥段为3区。分区的原理和方法见文献。 2 计算表面蒸发e的公式如下所示 E=3100Ρm3100+1.8Ρ2mexp(1-34.4Τm235.0+Τm)(2) 式中:Pm为月降水量,Tm为月平均气温。 3 序列值 υ=σˉx,其中:{ˉx=1nn∑i=1xi;σ=√1nn∑i=1(xi-ˉx)2(3) 式中:xi为样本第i年序列值;ˉx为序列的多年平均值;σ为标准差。 4 计算q0模型系数的公式如下所示 Q0=QiQj(4) 式中:Q0为某时段流量 (或降水量、蒸发量等)Qi与另一时段流量 (或降水量、蒸发量等)Qj的比值。 5 回归模型的建立 y′(t)=b0+b1t,其中b1n∑i=1(yi-ˉy)(ti-ˉt)n∑i=1(ti-ˉt)2(5) 计算冬季、雨季和过渡季节降水和河川流量的历史序列, 利用计算的要素时间序列, 以时间为自变量, 要素为因变量, 建立一元回归 (方程略) 或N阶回归曲线模拟方程。计算模拟序列与原序列的相关系数、气候的变化倾向率以及各年代的平均值, 比较其变化的特征。用气候诊断方法分析河川流量的历史演变规律和年代际变化特征, 探讨影响河川流量的因子, 以及其因子所起的作用。 2 年平均流量的年际变化 依据黄河源区不同地段0 ℃积温的开始日期和结束日期、雨季开始日期和结束日期以及达日、玛曲、唐乃亥水文站月际流量演变特征, 将黄河源区流量分成4个季节来讨论其变化的特点, 其中, 5—10月和12月至翌年3月分别为雨季和枯季, 4月和11月为过渡季节。 计算得出, 1956—2005年枯季、雨季和过渡季节4月及11月均方差的变差系数分别为0.21、0.30、0.24、0.31, 枯季变差系数最小, 过渡季节的11月变差系数最大。这说明枯季和过渡季节4月流量年际之间的波动幅度小, 比较稳定, 而雨季和过渡季节11月流量年际之间的波动幅度大, 不稳定。0.21～0.31的变差系数也说明, 黄河上游的干流是以雨水、冰川融水、季节性融雪水和地下水为主要补给源的季节性河流。 从唐乃亥年平均流量年代际变化 (图2b) 可以看出, 20世纪50—80年代年平均流量呈波动性的上升趋势, 90年代至21世纪前5年年平均流量呈下降趋势, 年平均流量基本经历了一个“枯→丰→枯”的历史变化过程。从10 a尺度看, 90年代是年平均流量的最枯时段, 最丰的时段是80年代。以1971—2000年黄河上游地区年流量的30 a均值为基准, 地表径流20世纪90年代累计减少了320.16亿m3, 21世纪