大凌河中下游水文站月径流分析.docxVIP

大凌河中下游水文站月径流分析 1 辽宁主要河流水系 大凌河位于辽宁省西部。南支起源于建昌县黑山以北的水泉沟，西支起源于河北省平泉县的泉子沟。它们流经沙左县大城市东南部的大阳iqa。流经朝阳、北票、义县、凌海等县市，向渤海注入凌海市南部的南魁河和南井子。流域面积23 263 km2, 辽宁省境内流域面积20 285 km2, 河长435 km, 主要支流有老虎山河、牤牛河及西河等。流域内多属荒山秃岭, 水土流失严重, 河流含沙量大, 气候干旱, 暴雨集中, 洪旱灾害较为频繁。朝阳水文站位于大凌河的中游, 控制流域面积10 236 km2。 2 流量变化的特征分析 2.1 分析年度变化的特征分析 2.1.1 模型:单峰形,20世纪2020年代 计算朝阳站各年代及多年平均月径流量结果为:径流量年内分配不均匀, 各年代月径流量年内变化情况均呈现为单峰形, 1～5月径流量小且变化平缓, 6月份开始迅速增加, 8月份达到最大, 10月份径流量开始减少;20世纪60～70年代较多年平均峰值偏前, 80～90年代峰值明显减少;60～70年代线形偏尖瘦, 80～90年代峰形较宽平。各年代6～9月的汛期径流量占年径流量的比值为60%～66%。 2.1.2 合成向量和合成总效应的公式 按径流量年内分配集中的程度和最大径流量出现的时间计算, 取1 a内各月的径流量作为向量, 月径流量的大小作为向量的长度, 所处的月份作为向量的方向。1～12月, 每月的方位角分别为15°, 45°, 75°, …, 345°。1～12月逐个累加各个向量, 合成的分向量的大小表示分向量之和的总效应, 合成向量的方向表示总效应的方向, 其计算公式为: Rx=12∑i=1risinθi(1)Rx=∑i=112risinθi(1) Ry=12∑i=1ricosθi(2)Ry=∑i=112ricosθi(2) R=√R2x+R2y(3)R=R2x+R2y???????√(3) Cn=R/12∑i=1ri(4)Cn=R/∑i=112ri(4) D=arctan (Ry/Rx) (5) 式中:ri、θi为月径流量的大小和方向;Rx、Ry、分别为X、Y方向合成向量和合成总向量;Cn和D分别为集中度和集中期。 经计算, 朝阳站径流在20世纪50年代、60～70年代、80～90年代、21世纪初的集中期分别为227°、231°、225°、229°, 集中度分别为0.49、0.48、0.44、0.38, 根据集中期的变化范围可推算出其最大径流出现在8月中旬至8月下旬, 这与实际情况相符合。 2.2 年龄变化的特征分析 2.2.1 年际变化程度 径流年际变化的总体特征常用变差系数Cv或年际极值比来表示。Cv和年际极值比反映了一个地区径流的相对变化程度, 它们的值越大, 表示径流的年际变化越剧烈;反之表示径流年际变化平缓。对朝阳站1954～2008年的径流量进行统计分析发现, 年最大径流量为23.44亿m3, 最小为0.7 952亿m3, 最大值与最小值的比值为29.5, 年径流变差系数Cv为0.725, 说明该地区径流年际变化较大。 2.2.2 丰水段时段 按照径流的距平百分率P来划分径流的丰枯等级, 计算表达式为: p=(wi-ˉw)/ˉw×100%(6) 式中:wi为某年的径流量;ˉw为多年平均径流量。P≥20%为丰水;10%≤P20%为偏丰;-10%≤P10%为平水;-20%≤P-10%为偏枯;P -20%为枯水。 经分析计算, 朝阳站55 a径流资料, 丰水年16 a, 偏丰年2 a, 平水年8 a, 偏枯水年为3 a, 枯水年26 a。变化过程中, 丰水段历时较短, 通常只有1～5 a;而枯水时段历时很长, 最长为10 a。在20世纪60 年代中期以前, 流域处于丰水期时段; 60年代中期以后, 流域基本处于枯水或偏枯期。 2.3 水文时间序列kjt秩次相关检验 首先采用累积滤波器对朝阳站径流进行趋势分析: S=n∑i=1Ρi/nˉΡ(n=1?2?3?n)(7) 式中:S为累计平均值;Pi为时间 (年) 序列;ˉΡ为时间序列的平均值;n为序列的长度。 计算得出朝阳站年径流累积平均曲线见图1。 朝阳站年径流虽年份不同呈现波动性变化, 但总体上呈明显下降趋势。 然后定量分析径流年际变化, 采用水文时间序列Kendall秩次相关检验方法。方法为:对序列x1,x2, …,xn(n为样本数) , 所有对偶观测值 (xi,xj,ji)xixj中出现的个数Pi, 顺序 (i,j) 的子集是 (i=1,j=2, 3, …,n) , (i=2,j=3, 4, …,n) , (i=n-1,j=n) , 则有: {U=τ/[var(τ)]1/2τ=4∑Ρi/n(n-1)-1var(τ)=2(2