黄河上游典型河流水热耦合平衡方程研究.docxVIP

黄河上游典型河流水热耦合平衡方程研究 越来越激烈的人类活动严重改变了流域的底部，并影响了流域的产配流机制，导致河流等水文因素的变化，从而影响水循环过程和水量平衡。 针对河川径流变化影响因素的评估方法已有较多研究，如广泛应用的水文模型法、水文统计法等 洮河、大夏河作为黄河上游区重要的一级支流，受气候变化和人类活动双重作用，河川径流发生显著的变化 研究区域的概况和数据来源 大夏河流域概况 洮河流域位于甘肃省西部，介于101°36′— 104°20′E、34°06′—36°01′N(图1),发源于西倾山北麓，由西向东依次经过碌曲、卓尼、岷县，岷县受阻后北拐穿过九甸峡和海甸峡，经东乡、临洮等地，在永靖注入黄河。流域总面积24973 km 大夏河发源于甘南高原，河源由夏河和咯河(又名：德乌鲁河)两大支流组成(图1)。大夏河流经夏河县、合作市、临夏、东乡县注入刘家峡水库，流域面积6808 km 月径流、月径流数据和气象数据分析 本文选择洮河流域的红旗水文站和大夏河流域的折桥水文站(图1),通过黄河水文年鉴、《黄河泥沙公报》及黄土高原数据中心(http: ///hyd/inedx.jsp)获取两个水文站的月径流数据。通过中国气象数据共享网(http: //)下载获得流域范围内12个气象站点1961—2017年的基础气象数据，包括最高温度、最低温度、平均温度、日照时数、平均相对湿度、日降雨量、平均风速。采用世界粮农组织(FAO)修正的Penman-Monteith公式获取研究区的潜在蒸散发，该方法充分考虑了作物地表特征和空气动力学的参数的变化 (1) 式中，ET 学习方法 实行endell突变检验分析 Mann-Kendall(MK)趋势检验和Mann-Kendall突变检验分析是水文时间序列常用的方法，能够很好地反映河川径流变化发展规律，具有精确度高、不受异常值干扰、适用于非正态分布数据的优点。计算方法详见参考文献 径流变化归因识别 基于Budyko水热耦合平衡理论量化气候变异和人类干扰对径流变化的敏感性和贡献率，能够准确、快速识别引导径流变化的主导因素，在研究气候变化和下垫面变化对径流影响程度分析中广泛应用 径流弹性系数法 基于Budyko水热耦合平衡方程计算长时间尺度下的实际蒸散发 (2) 式中：ET表示流域年实际蒸散发( 在多年尺度上，闭合流域满足水量平衡关系，即： R=P-ET (3) 式中，R为多年平均径流深( (4) 假设P、ET (5) 定义 (6) 进一步推导，各弹性系数表达式为(?=ET (7) (8) (9) 降雨、潜在蒸散发、下垫面的径流弹性系数由(7)—(9)式求得，弹性系数的正负表示径流与该因子的正负相关关系，弹性系数的绝对值大小表示该因子对径流变化的影响程度。 径流时空变化 如果径流发生显著性变化，基于突变点，研究时期分为基准期和若干人类活动期，基准期和人类活动期的径流深变化为 (10) 式中i代表降雨P、潜在蒸散发ET 通过计算可以分别获得人类活动期降雨、潜在蒸散发、下垫面系数相对基准期引起的径流深变化为ΔR η 其中η 结果与分析 径流时空分布 洮河和大夏河流域的径流特征如表1所示，结果表明：洮河流域的多年平均径流深为177.71 mm, 是大夏河121.49 mm的1.41倍。大夏河年际径流深极值比为22.66,是洮河的5—6倍，表明大夏河年际径流深差异明显，这主要与大夏河南北地域差异性和季风性气候引起的降雨时空分布有关。采用非参数MK趋势检验流域径流深的变化趋势可知：洮河、大夏河流域的径流深均呈现减少趋势，显著性超过99%。图2为流域年径流深的线性拟合，由此可知，洮河、大夏河1961—2017年径流深减少速率分别为-1.85 mm/a、-1.36 mm/a。Mann-Kendall突变点检验发现，洮河、大夏河径流突变分别发生在1987、1985年，突变之后人类活动期径流相对变化率超过30%。 流域降雨、蒸发等气象因素直接影响着河川径流变化。采用MK趋势检验分析研究区降雨、潜在蒸散发变化特征发现(表1),洮河、大夏河年降雨量无显著的变化趋势，年降雨量变化速率分别为-1.08、0.005 mm/a(图3);而洮河和大夏河流域年潜在蒸散发显著增长( 流量变化的归因分析 蒸散发、下垫面参数的弹性系数分析 基于Budyko水热耦合平衡方程理论计算洮河流域和大夏河流域各时段水文特征以及径流对降雨量、潜在蒸散发、下垫面参数的弹性系数(表2)。分析发现：洮河、大夏河流域在人类活动期降雨量、径流深相对基准期均减少，年潜在蒸散发增加，径流系数减少，降雨转换为蒸散发的比例增加。 总体而言，洮河、大夏河流域径流