[发明专利]一种准确插补延长设计站水文特征值的处理方法

权利要求书

1.一种准确插补延长设计站水文特征的处理方法，其特征在于包括以下步骤：

(1)简单线性相关分析

设计站自然特征实测物理数值为m个的短系列Y₁，在设计站临近地区寻找相同或相近自然特征、具有物理成因关系的参证站，参证站的自然特征实测物理数值为n个的长系列X，其中包括与Y₁时间对应的实测物理数值系列X₁和超过设计站的实测物理数值系列X₂；

建立Y₁与X₁简单线性相关关系，插补延长设计站中与参证站时间对应但缺测的物理数值Y₂；由于受简单线性相关的技术限制，无法得到完整的Y₂，得到的是被均化的y₂；

y₂＝aX+b…………………………………………………①

(2)完整线性相关分析

在Y₁与X₁线性相关关系方程中考虑偏值Δy₂：

Y₂＝y₂+Δy₂＝aX+b+Δy₂………………………………………②

式中：Y₂——设计站完整插补延长的物理数值；

Δy₂——完整插补延长物理数值的偏值；

(3)偏值计算

根据设计站、参证站对应年的实测物理数值相关数据点，以及分析确定的简单相关方程，计算出每一个相关点的偏值Δy_1j，j＝1、2、3、…、m，通过m个Δy₁值可统计出样本的算术平均值标准差S，近似代替总体的数学期望α，S近似代替总体的均方差σ；根据正态分布函数，通过随机取得Δy₂；

Δy₁＝Y₁-(aX+b)………………………………………………③

j＝1、2、3、…、m；m﹤n

(4)完整线性相关的物理数值计算

根据偏值计算设计站完整插补延长的物理数值：

Y₂＝y₂+Δy₂＝aX+b+Δy₂………………………………………⑦

Y₂∪Y₁＝Y…………………………………………………⑧

式中，Cv₂为完整长系列Y的离散系数、Cv为X的离散系数、k_i为模数、c为绝对相对误差，若c＜3％，则Cv₂与Cv符合地区水文特征，计算结束；若c≥3％，则返回重新计算偏值。