1. 水文短系列插补延长与问题提出
我院在编制《勐波罗河流域水资源综合规划》中,云南境内的怒江最大支流勐波罗河主流上设有上(上游站)、中(下游站)、下(某站)3个水文站,上游站仅实测有连续15年(1980~1994年)短系列径流量(Y1),下游站实测有连续53年(1961~2013年)长系列径流量。为了插补延长上游站年径流量短系列为长系列,首先根据水量平衡法对受人类活动影响较大的下游站实测流量还原为天然径流量(X)系列,然后建立上游站与下游站的同期1980~1994年径流量相关分析 [1] [2] (见图1),相关数据点群(○)呈狭窄带状分布,二者相关密切(相关系数0.88)。通过狭窄带状数据点群(○)的不同量级重心或最小二乘法可确定相关线性关系,其相关方程(1)式为简单相关:
(1)
式中:
为上游站插补延长的年径流量,X代表下游站年径流量,a、b表示系数。
按简单相关方程插补延长上游站38年(1961~1979、1995~2013年)径流量(
),
与下游站对应年径流量的关系数据点(▲)完全落在简单线性相关关系线上,这是不符合两站年径流量实际情况的(提出问题),这是目前国内工程水文分析通常采用的方法,也是技术方面一直存在的不足。
2. 插补延长对系列统计参数 [3] 的影响
简单插补延长的38年上游站年径流量属于同量级中值,38个上游站与下游站年径流量关系数据点(▲)完全落在简单线性相关关系线上,简单插补延长年径流量(
)和实测年径流量(
)组成53年径流量(y)简单长系列,其均值
(4025)基本代表总体特征,离散系数
值为0.28,小于下游站年径流量(X)长系列离散系数
值0.29 (公式略)。上游站、下游站为上、下游关系,后者流域面积是前者的近20倍,根据径流上、下游特性,前者年径流量
值应大于后者
才合理。
综合分析,上游站简单相关方程插补延长的38年径流量(
),与实测年径流量(
)组成的年径流量(y)简单长系列被均化,其离散系数
值变小。
![](//html.hanspub.org/file/6-2410637x24_hanspub.png)
Figure 1. Simple correlation analysis of annual runoff for up and down stream stations
图1. 上游站与下游站年径流量简单相关分析图
(2)
(3)
(4)
(5)
式中:n为年径流量系列长度(n = 53),
代表上游站多年平均径流量,
表示上游站对应第i年径流量,
为模数,
表示上游站年径流量系列离散系数,
。
3. 均化系列的离散化处理
上游站插补延长简单相关方程(1)式反映了上游站与下游站对应年的年径流量线性相关关系,从相关数据点群(○)的分布形状看,完整表达上游站与下游站对应年的年径流量相关关系应该是下(6)式 [4] :
(6)
式中:
为上游站完整插补延长的年径流量,
表示完整插补延长的年径流量偏值。
上游站完整插补延长年径流量(
)系列与实测年径流量(
)系列组合成完整年径流量(Y)系列,即
。
(6)式中
项属于两站年径流量相关点的偏值,偏值是由于两站流域内降水场次大小不一、时空分布不均以及下垫面条件差异等因素造成,水文分析中一般忽略偏值
项,这对插补延长与实测组成的长系列均值是没有影响的,但对离散系数Cv值有影响,且不宜忽略。
分析偏值
项的影响因素,其具有随机性,由于Y系列的随机性为皮尔逊III型分布曲线,偏值
项的随机性基本服从正态分布,反映了关系数据点群偏离线性关系线近则密、远则疏的特点。
根据上游站、下游站对应年的实测年径流量相关数据点,以及分析确定的简单相关方程(1)式,可计算出每一个相关点的偏值
,见图2,通过15个
值可统计出样本的算术平均值
(−29)、标准差S (484),
近似代替总体的数学期望
,S近似代替总体的均方差
。根据正态分布函数(11)式 [5] ,通过随机取得
[即函数
值] [6] ,见表1。
(7)
(8)
(9)
(10)
(11)
式中:符号意义同前。在上述简单相关方程的基础上,考虑其造成长系列离散系数变小,增加偏值
相后,完整相关方程为(11)式,此式反映了上游站与下游站年径流量相关的非线性关系,完整插补延长的38年上游站年径流量与下游站同期径流量数据点(◆)分布见图3,从图可看出完整插补延长数据点(◆)分布在实测相关数据点群(○)的分布范围内,分布特性与实测点数据群(○)基本一致,完整插补延长数据点群(◆)偏离简单线性相关关系线呈近则密、远则疏的水文特点。
上游站完整插补延长年径流量(
)与实测年径流量(
)组成年径流量完整长系列(Y)的离散系数
值为0.32 (公式略),比简单长系列(y)的离散系数
值(0.28)大,比下游站年径流量系列离散系数
值0.29大0.03,解决了简单相关插补延长带来的均化不利影响,符合上游站、下游站的上下游水文特性 [7] ,具有区域上的合理
![](//html.hanspub.org/file/6-2410637x63_hanspub.png)
Figure 2. Diagram of deviation of simple correlation date of annual runoff for up and down stream stations
图2. 上游站与下游站年径流量简单相关点偏值示意
![](Images/Table_Tmp.jpg)
Table 1. Interpolation and prolongation results of annual runoff for up and down stream stations
表1. 上游站与下游站年径流量简单和完整插补延长成果表
![](//html.hanspub.org/file/6-2410637x64_hanspub.png)
Figure 3. Diagram of entire correlation of annual runoff for up and down stream stations
图3. 上游站与下游站年径流量完整相关分析图
![](//html.hanspub.org/file/6-2410637x65_hanspub.png)
Figure 4. Diagram of experiential data of annual runoff for upstream station
图4. 上游站年径流量经验频率数据点频率分布图
![](//html.hanspub.org/file/6-2410637x66_hanspub.png)
Figure 5. Simple correlation analysis of annual runoff for a station with downstream station
图5. 某站与下游站年径流量简单相关分析图
性,属目前工程水文分析的创新点。
上游站年径流量的完整长系列(Y)、简单长系列(y)的均值仅相差0.04%,完整相关分析计入偏值项不影响均值稳定。上游站年径流量完整系列验频率数据点(▲)分布明显比简单系列验频率数据点(○)离散,见图4。
4. 均化系列离散处理方法的验证 [8]
下游站的下游有与之同步的某站实测流量,且经过还原的天然实测径流量(Z)长系列。假定某站类似上游站
![](//html.hanspub.org/file/6-2410637x67_hanspub.png)
Figure 6. Diagram of entire correlation of annual runoff for a station with downstream station
图6. 某站与下游站年径流量完整相关分析图
缺相应年份径流量,按上述方法进行简单和完整的相关分析插补延长某站径流系列,图5为偏值分析图,验证分析结论是:某站完整插补延长长系列均值与实际长系列均值仅相差0.2%,二者长系列的离散系数Cv值(0.23,公式略)相同,插补延长数据点(Δ)与相关数据点(□)分布均偏离简单线性相关关系线呈近则密、远则疏的水文特点,见图6。
通过对某站验证分析,说明完整相关分析插补延长方法对用简单插补延长被均化物理特征值系列具有可推广应用价值。
5. 均化系列离散处理的应用意义
1) 本文对均化系列的离散化处理,解决了线性相关插补延长多年来存在的缺陷,使插补延长系列较好的反映实际水文特征;
2) 偏值属于正态分布的分布函数数值,每选取一组,分析得的插补延长系列与上一组是有差异的,只要考虑完整长系列与简单长系列的均值差异在允许范围内(一般小于3%),以及Cv值符合地区水文特性即可;
3) 只要两系列具有物理成因的相关关系,简单线性相关关系插补延长系列被均化影响,均可考虑本法消除其影响,目前未见国内有相似的处理方法。
基金项目
国家水体污染控制与治理科技重大专项“滇池流域水资源联合调度改善湖体水质关键技术与工程示范”(2013ZX07102-006-01);云南省技术创新人才计划(2011C1092)。