文峪河径流演变规律的分析.docVIP

文峪河径流演变规律分析 摘要：在分析文峪河流域特征和流域生态环境变化的基础上，应用HHT变换方法和非趋势波动分析法等方法，分析文峪河径流演变规律。结果表明：文峪河水库坝址处的天然径流演化与文峪河水库坝址以上流域的降水、气温等气象因素的周期变化规律是相吻合的，受气候因子的影响较为强烈，且表现出与气象要素相同或相近的周期变化特征。结合EMD分解得到的趋势走向，可推知未来文峪河年径流将存在一个持续的缩减趋势，径流缩减主要是受到区域降雨减少和气温上升的影响，其中以降雨偏低为主。在气候等因素变化的前提下，必须加强地表水资源的保护和调控。 关键词：文峪河 径流演变 分析 1文峪河概况及流域特征 1.1 概况 文峪河山西省中部西侧属黄河流域汾河水系，是汾河的最大支流。主河道全长158.6公里，上游河段93.9公里，中下游河段64.7公里。流域总面积4112.4平方公里。文峪河是山西省吕梁地区的重要水源地之一，承担了交城、文水、汾阳和孝义平川四县市、近50万人口的供水任务，在整个吕梁区域水安全保障体系中占有极其重要的战略地位。 1.2流域特征 1.2.1 气候特征 文峪河流域属暖温带大陆性半干旱气候，四季分明，春季多风、少雨、干燥；夏季多暴雨、炎热；秋季旱涝多变、早凉；冬季少雪、寒冷。极端最高气温39.9℃，最低极端气温-27.4℃。无霜期160-183天，最大冻土深0.8米。平均降水量523.9毫米。降水量呈现年际变化大、年内分配不均、空间垂直分布变化大等特点。整体呈现减少的趋势，常出现旱涝年相间。年内分配差异大，最大日降水和最大月降水，大多出现在6-9月，个别年份出现在5月或10月。汛期4个月降水量集中，占全年降水量的50%-88%。 1.2.2 水文特征 流域多年平均水资源总量为3.83亿m3，其中地表水为2.34亿m3，地下水为2.67亿m3。其中，多年平均径流量1952年～1996年间，实测年平均径流量1.741亿m3，最大为4.78亿m3（1964年），最小为0.596亿m3（1987年）。20世纪70年代以后，径流量呈下降趋势。结冰期为11月～次年3月。企业节省费用、降低成本、躲避环境监管，次将含水污泥用高压水泵绕过在线监控点偷排至文峪河，造成文峪河水质恶化、水体污染，影响了民众正常的生产生活 3.1数据来源 文中采用文峪河水库坝址处1952年-2007年共56年的天然年、月径流为基础数据，数据来源于山西省水文局 3.2 研究方法Hilbert-Huang变换方法HHT（Hilbert-Huang Transform）是一种适用于非线性和非平稳性数据的处理方法，主要分为经验模态分解法（EMD）及Hilbert变换两部分。 （1）EMD方法 EMD是美国N.E. Huang提出的信号分解算法，此算法的目的在于把复杂信号分解为简单的单分量信号的组合，即将性能不好的信号分解为一组性能较好的具有瞬时频率的固有模态函数（IMF）和具有单调性或一个极值点的残余项Rn。计算步骤： 原时间序列x(t)的各局部极值然后用三次样条函数进行插值，得到原时间序列x(t)的上包络线序列值max(t)和下包络线序列值min(t)。每个时刻max(t)和min(t)的瞬时平均值m(t)： (1) 用原序列x(t)减去瞬时平均值m(t)，得到新数据系列： (2) 得到第一个固有模态函数I1(t)后，用原序列减去I1(t)，得到剩余值序列r1(t)： () 将r1(t)作为一个新的原序列，以上步骤，依次提取I(t)、I3(t) …In(t)，直至没有固有模态函数被提取出来。 分解各分量合并，得原序列x(t)： () 把一个数据分解成为固有模态函数组和残余量之和。 (5) 式中，P表示取积分的主值。 解析信号z(t)可定义如下： (6) 其中， (7) (8) 式（7），（8）分别为解析信号z(t)的瞬时振幅和瞬时相位。在此基础上，可以把瞬时频率ω(t)定义为： (9) 通过上式定义的瞬时频率是时间的单值函数，在任何时间里有唯一确定的瞬时频率。为了使瞬时频率具有实际意义，必须满足在任何信号经变换分解后，得到的各个