DL∕T 1961-2019 火电厂流量测量不确定度计算方法(电力).pdf

《DL∕T 1961-2019 火电厂流量测量不确定度计算方法(电力).pdf》由会员分享，可在线阅读，更多相关《DL∕T 1961-2019 火电厂流量测量不确定度计算方法(电力).pdf（33页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、DL/T 1961 2019 I ICS 27.060.30 J 75 备案号：-20 中华人民共和国电力行业标准 DL/T 1961 2019 火力发电厂流量测量不确定度计算方法 Calculation method of uncertainty in flow measurement of power plant 火力发电厂锅炉耐火材料 火力发电厂锅炉耐火材料 2019-06-04发布 2019-10-01实施 国家能源局 发 布 DL/T 1961 2019 I 目次 前言.II 1 范围.1 2 规范性引用文件.1 3 术语和定义.1 4 参数符号.2 5 不确定度评定的一般性原则.3

2、 6 独立测量中的不确定度的计算方法.4 7 线性拟合和不确定度计算方法.5 8 非线性拟合和不确定度计算方法.10 附录 A(规范性附录)一般函数标准差的计算.13 附录 B(资料性附录)明渠校准的算例.14 附录 C(资料性附录)封闭管道流量校准的不确定度算例.19 附录 D(规范性附录)回归方法.26 附录 E(资料性附录)正交多项式曲线拟合.29 DL/T 1961 2019 II 前言 本标准按照GB/T 1.1-2009标准化工作导则第1部分：标准的结构和编写给出的规则起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别这些专利的责任。本标准由中国电力企业联合会提出

3、。本标准由电力行业电站汽轮机标准化技术委员会归口。本标准起草单位：西安热工研究院有限公司、西安西热节能技术有限公司。本标准主要起草人：高登攀、曾立飞、杨荣组、祁文玉、张永海、余小兵、谷伟伟、朱蓬勃、石 慧、王 汀、高 庆，薛朝囡、穆祺伟。本标准为首次制订。本标准在执行过程中的意见或建议反馈至中国电力企业联合会标准化管理中心（北京市白广路二条1号，100761）DL/T 1961 2019 1 火力发电厂流量测量不确定度计算方法 1 范围 本标准规定了火电厂流量测量在校准和使用过程中的线性和非线性两类关系的不确定度计算方法。本标准适用于各类封闭管路或明渠的流量测量。2 规范性引用文件 下列文件对

4、于本标准的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用于本标准。凡是不注日期的引用文件，其最新版本（包括所有的修改单）适用于本标准。GB/T27759-2011 流体流量测量不确定度评定程序 ISO772 水力学测量词汇与符号 ISO-TR-7066-1 流量测量装置校准和使用方法不确定度的估计第一部分：线性校准关系 ISO-TR-7066-2 流量测量装置校准和使用方法不确定度的估计第二部分：非线性校准关系 3 术语和定义 3.1 校准图校准图calibration graph 流量计的某些响应参数和流量函数为坐标所得的点绘曲线。3.2 置信限置信限 confidence li

5、mits 观测值或计算值的（置信）上限和下限，由样本统计量所构造的总体参数的估计区间。3.3 相关系数相关系数 correlation coefficient 表示两个变量线性关系程度的指标。3.4 方差方差 variance 表征数据分布，衡量变量值与期望值之间的偏离程度，定义为()()2()1iVar xxxn=3.5 协方差协方差 covariance 两个变量之间的总体误差的期望，定义为()cov(,)()(1)iix yxxyyn=3.6 测量误差测量误差 error of measurement DL/T 1961 2019 2 测量值与真实值之间的差值，包含系统误差和随机误差。3

6、.7 随机误差随机误差 random error 测量结果与同一待测量的大量重复测量的平均结果之差。3.8 系统误差系统误差 systematic error 在重复性条件下，对同一被测量进行无限多次测量所得结果的平均值与被测量的真值之差。3.9 粗大误差粗大误差 spurious error 超出规定条件下的误差。3.10 残差残差 residual 指实际观察值与估计值之间的差。3.11 样本样本(试验试验)标准差标准差 sample（experimental）standard deviation 同一被测量 n 次测量结果平均值的分散性度量，定义为()0.52()()1is xxxn=3

7、.12 不确定度不确定度 uncertainty 与测量结果关联的一个参数，用于表征合理赋予被测量值的分散性。3.13 独立测量独立测量 individual measurement 仅需测量单一变量的测量过程。3.14 无偏拟合无偏拟合 unbiased fitting 拟合曲线能真实的反应数据的分布趋势，用于评估拟合程度的优良特性。4 参数符号 4.1 参数符号定义 本标准中的符号及测量值单位应符合表 1 的规定。表 1 符号及说明 符号 名称 a 校准曲线在坐标轴上的截距 b 拟合直线的斜率 c 加权最小二乘法中的系数 DL/T 1961 2019 3 Cov 变量的协方差 eR 变量的

8、随机不确定度 eS 变量的系统不确定度 gj 第 j 个正交多项式的系数 Q 流量 r 相关系数 s 标准偏差 sR 最佳拟合直线上各点的标准差 t 学生分布 wi 加权最小二乘法第 i 个加权因子 x 自变量 y 因变量 U 总不确定度 UADD 根据加法模型获得的不确定度，置信概率 95%99%URSS 根据平方和的根的模型获得的不确定度，置信概率 95%自变量 x 的标准差与因变量 y 的标准差的比值 测量值与计算值的偏差 样本平均值 样本标准差 灵敏度系数 4.2 参数符号上、下标定义 本标准参数符号的上标、下标定义应符合表 2 的要求 表 2 参数符号上标、下标及定义()变量的平均值

9、()由拟合曲线所确定的预测值 i 变量的第 i 个值 ij 矩阵第 j 列的第 i 个值 5 不确定度评定的一般性原则 a)不确定度评估：评估基本数学关系式的形式和评估拟合曲线的不确定度。b)独立测量中的随机不确定度与系统不确定度按第 6 节给出的方法计算。c)两个变量之一的误差忽略不计时，以及自变量的误差相对因变量的误差可以忽略不计时，应使用 7.3 至 7.8 节所描述的方法进行线性拟合和不确定度计算，所用的基本方程形式如式（1）：y=a+bx（1）DL/T 1961 2019 4 式中：x自变量；a，b拟合直线的系数。自变量 x 的误差较小时，在评估过程中可将 x 设定为一个预定值，应使

10、用 7.1 节的计算方法。e)变量 y 独立于变量 x 时，即拟合直线与 x 轴平行，这种情况应采用第 7.5 节描述的方法评估其不确定度。f)当变量关系为曲线时，则宜考虑将其转化为线性关系。若不能转化则使用 8 节所描述的非线性拟合方法。6 独立测量中的不确定度的计算方法 6.1 随机不确定度的计算 6.1.1 对每个变量列出不确定度的来源并制作表格，表格应包括所有测量中的误差，并分别列出随机误差项和系统误差项。6.1.2 对于直接测量所确定的变量值可以采用式（2）或式（3）计算标准差。s=(xi x)2/(n 1)0.5（2）或者 s(x)=nxi2(xi)2/n(n 1)0.5（3）6.

11、1.3 将标准差计算结果代入式（4），可得到被测量 x 的某个特定 t 值的随机不确定度。eR(x)=ts(x)（4）6.1.4 为了提高评估的准确度，必须采取足够多的数据点。6.1.5 变量为多个独立测量结果的和或差时，根据下式计算总体的标准差。s(x)=s(xi)20.5（5）之后代入式（4）计算不确定度。6.1.6 对更复杂的函数，如因变量和自变量之间是积或商等较复杂的关系时，标准差应按照附录 A 计算，之后再代入式（4）计算相应的不确定度。6.2 系统不确定度的计算 6.2.1 系统误差中的无法确定的分量应根据以往有效校准数据和历史记录等进行评估。6.2.2 当变量为多个独立测量结果的

12、和，且各分量的正负不确定时，系统误差用下式计算：()0.52,eSS iie=（6）DL/T 1961 2019 5 6.2.3 对于更复杂的函数关系，系统误差限应用附录 A 中给出的方法，用2,eS i替换公式中的方差项。6.2.4 当确定了所有不确定度的来源，并且每个变量都合成了随机不确定度和系统不确定度，则完成了评估过程。6.2.5 明渠类和封闭管道中流量测量的不确定度评估算例参见附录 B 和附录 C。7 线性拟合和不确定度计算方法 7.1 校准图表的线性判断 为了确定拟合值与测量值能否无偏拟合，应通过观察测量结果和拟合直线的残差，采用下述方法获得近似直线：a)数据按照 x 或 y 的升

13、序排列，两个变量的平均值由（7）式确定。;iixx nyy n=（7）将数据等分为独立的三组，算出两端的两组数据平均值，分别记为x 1，y 1和x 3，y 3，则近似直线的斜率为 3131()()byyxx=（8）由于直线必须通过总体平均值x 和y，因此直线方程可以由下式获得：()()iiyyb xx=（9）用y bx =a进行替换，得到 y 的最佳估计值 iiyabx=+（10）最后，可用下式确定残差 ()()()iiiiiyyyyabx=（11）另一种线性判断方法，可以采用 7.3.2 描述的最小二乘法中的相关系数进行判断，其误差限同样按照式(11)计算。b)试验数据分布的初步判断。将计算

14、的残差结果按升序排列，并在标准概率纸上描绘成累积频率曲线。如果数据点基本处于一条直线附近，且没有明显弯曲，则可认为数据点近似为正态分布。c)检查异常大或异常小的残差，因为它们的存在会严重影响最终拟合直线的定位并不可避免地增大不确定度。为了识别这些异常值，可采用 GB/T27759-2011 的附录 D 中描述的方法。在对所有相关因素进行仔细分析的基础上，可适当剔除残差异常的观测数据，然后重新拟合直线和计算残差。d)描绘残差y相对自变量 x 的曲线，此曲线和估计值y 进行比较，若曲线满足下列任一情况：a)数学关系恰当合理；b)拟合过程正确实施；c)方差随 x 没有明显变化。则上述描绘的点应处于均

15、一宽度的水平区间内(见图 1 a)。残差分布也可能出现以下三种偏离理想模式的情形：1)分布带出现明显的向上或向下的弯曲(见图 1 b)，表明变量间的关系为曲线而非直线。2)分布带依然保持水平，但逐渐变宽或者收窄(见图 1 c)，这说明整个测量范围内的方差不是常数，DL/T 1961 2019 6 应在拟合过程中进行加权计算。3)分布带表现为向上或向下的直线带(见图 1 d)，表明在拟合过程或 y 的相关计算存在误差。图 1 残差、自变量、预测值关系曲线 7.2 校准曲线的线性化方法 7.2.1 如果根据第 7.1 节线性判断表明校准图为曲线形式，则考虑按照以下两种方法进行线性化处理：a)第一种

16、方法是变量置换，仅适用于变量之间存在数学函数关系性质，转换的形式取决于数学表达式本身。例如，对于明渠校准，水位与流量的关系可以表达为：0()bQc hh=+（12）式中：h测量的水位，单位为 m；h0零流量水位的基准面修正值，单位为 m；c系数 b指数 采用对数变换将上式改写为：0lnln*ln()Qcbhh=+（13）上式实现了数据线性化处理。b)第二种方法，将校准曲线分为几部分，若每个部分可线性化处理，则校准曲线依旧可以线性化。此方法需满足两个条件：第一，每部分曲线应基于相类似数量的观测点；第二，曲线的每一部分必须与DL/T 1961 2019 7 相邻部分共用 2 到 3 个点。7.2.

17、2 完成线性化后，应重复 7.1 所述的校准图标的线性判断。7.3 最佳线性拟合 7.3.1 根据 7.1 两个标准差的计算值，按照式（14）计算比值，()/()s ys x=（14）当比值20 时，应采用 7.3.2 给出的最小二乘法拟合直线。当比值小于 20 时，若变量能按 7.3 至 7.4 节的要求设定至预定值，仍可采用最小二乘拟合直线的方法。7.3.2 最小二乘法 a)若自变量的误差与因变量的误差相比可忽略不计，采用最小二乘法进行校准直线的拟合，可用下式计算：iiyabx=+（15）b)斜率 b 按照式（16）计算：()()()2iiibxxyyxx=（16）截距按照式（17）计算：

18、aybx=（17）表征 x 与 y 关系程度的相关系数 r 可按式（18）计算：()()()()0.522iiiirxxyyxxyy=（18）4)为完成拟合过程，拟合直线的标准差可按式（19）和式（20）计算：()()0.522iiryyn=（19）()()0.52a2iiybxn=（20）sR也可按式（21）计算：()0.52s1Rs yr=（21）上式中 ()()0.52s1iyyyn=（22）当采用式（21）时，需要足够多的数据点以减小舍入误差。5)b，r 和 s(y)也可通过以下方程计算：()22iiiiiibnx yxynxx=（23）DL/T 1961 2019 8 ()()()0

19、.522222iiiiiiiiirnx yxynxxnyyy=（24）()()0.522()(1)iinyys yn n=（25）对于 a 通过式(17)计算，这里仍需要由足够多的数据点以减小舍入误差的影响。6)当校准曲线由两段或多段组成，应确定这些区间的交界点。相邻两个区间的方程可按式（26）计算：y1=a1+b1x 和y2=a2+b2x（26）则在交界点有 y1=y2，x 的公值可根据式（27）计算：x=(a1 a2)/(b2 b1)（27）y 值可通过将 x 值代入式(26)来计算。7.4 最佳加权曲线的拟合 7.4.1 当 y 的方差随着 x 值变化，则需使用加权回归分析方法。此时 a

20、 和 b 通过式（28）和式（29）进行计算：iiiic ybc xan=（28）和 ()()()22iiiiiiiiiiic xc yc xbc x yc xnn=（29）通过式(22)计算 yi的标准差，通过式(19)或式(21)计算拟合直线的标准差。7.4.2 当已知各个变量的方差时，权重系数 ci可按式（30）计算：wi=1/var yi;ci=wi/w （30）其他情况下，权重系数的计算方法如下：a)按照第 7.1 所述，根据拟合直线用式（11）计算(yi)；b)以(yi)及对应的 xi值为坐标作图；c)采用本标准第 8 节给出的方法对数据进行非线性拟合；d)通过该拟合曲线得到 2(

21、yi)；e)用 2(yi)替代式(30)中的 yi计算得到 ci 的值。7.5 y 独立于 x 时的拟合直线确定方法 7.5.1 校准曲线斜率为零时，即在 x 的变化范围内 y 为常数，校准曲线变成一条水平直线，校准曲线简化为 yi的平均值。即 y=yi/n（31）7.5.2 通过试验进一步测试拟合曲线的斜率是否为零，按照式（32）计算，ts(b)（32）式中，s(b)按下式计算：DL/T 1961 2019 9 s(b)=sR/n 10.5s(x)（33）当 0 包含在式(32)给出的限制范围内时，说明拟合直线是水平的。7.6 随机不度定度的计算 7.6.1 拟合直线在 x=xk处的随机不确

22、定度按下式计算：eR(y)=tsR(1/n)+(xk x)2/(xi x)20.5（34）对于单次测量，x=xk时的不确定度为 eR(yk)=tsR1+(1/n)+(xk x)2/(xi x)20.5（35）7.7 系统不确定度计算 7.7.1 根据第 6 节给出的原则，计算校准过程中的系统不确定度，用式（36）计算：eS=es,i20.5（36）当变量的值基于诸如积或商等更复杂的函数关系时，则可采用附录 A 中方法。7.7.2 随机不确定度和系统不确定度可用以下两种公式合成：UADD=es+ts（37）或者 URSS=es2+(ts)20.5（38）7.8 测量结果的不确定度评估 7.8.1

23、 由拟合直线上的定位不准确形成的附加不确定度，以及由数据的增加或减少引入的其它附加不确定度，应使用 GB/T27759-2011 给出的方法进行评估。7.8.2 用URSS模型合成随机不确定度和系统不确定度，则测量过程中总的不确定度表示为URSS(y)，可按下式计算：URSS（y）=URSS2(y 0)+URSS2(y c)0.5（39）式中：URSS(y 0)总的附加不确定度 URSS(y c)由校准图表引起的不确定度 附加不确定度对流量总不确定度的影响由校准图表的特性决定。7.8.3 当拟合直线的斜率为 0 时，流量由流量计的输出函数与独立于流量的系数乘积得到，此时不存在附加不确定度，式(

24、39)简化为：URSS(y)=URSS(y c)（40）7.8.4 当拟合直线的斜率不等于 0 时，需用迭代法计算流量。为进行迭代，需用拟合系数初始估计值来获得流量初始计算值，然后使用计算的流量得到更加准确的校准系数，重复这一过程，直到流量的预估值不再发生较大变化。在这种情况下，流量计的任何测量误差都会引入所使用系数的误差，因此，总的不确定度应根据式(39)计算。DL/T 1961 2019 10 7.8.5 当使用条件与校准条件不同时，如系统测量条件、流体、装置等不同时，还会进一步增加测量的不确定度。在这种情况下，需要评估每种情形的置信限。8 非线性拟合和不确定度计算方法 8.1 非线性拟合

25、方法基本原则 8.1.1 如不能使用线性拟合方法，应通过多项式法建立变量之间的非线性校准曲线。例如二次多项式，形式为：y =b0+b1x+b2x2（41）多项式的一般表达形式为：y =b0+b1x+bjxj+bmxm 或 y =bjxjmj=0（42）8.1.2 利用最小二乘法计算系数bj，利用下式使校准曲线与数据点偏差的平方和最小：2i 1()niiyy=式中：y i为x=xi时由式(42)计算的预测值。8.1.3 多项式中最高次幂可根据以往经验获得，否则通过 8.3 中方法获得。8.1.4 拟合多项式项数过多可能导致曲线震荡，此种情况可将 x 的区间分为若干可采用线性或者低阶多项式拟合的区

26、间。8.1.5 对一个或两个变量进行适当变换也可实现线性或低阶多项式拟合。例如，将自变量转换为它的倒数 1/x 可以使原始数据线性化。8.1.6 当给定数据xi的随机不确定度er(x)相对于yi的随机不确定度er(y)不能忽略时，将不再适用最小二乘法。当校准曲线的斜率小于er(y)/er(x)的 1/5 时，可以认为该方法有效。当拟合曲线斜率超出这一范围时，本部分所描述的数学处理方法将不再适用。因此，在实际校准过程中，如果待拟合变量不满足上述条件，则本部分所描述的方法均不再适用。8.1.7 如果在拟合前对某个变量进行了变换，则不确定度与新变量相关。由于变量变换导致随机不确定度er(y)在 x

27、范围内不能作为常数看待，应需要采用加权最小二乘法进行拟合。8.2 非线性拟合计算方法 8.2.1 本标准所描述的拟合直线的方法也被称为线性或简单线性回归。与之类似，拟合多项式可以称为多项式或曲线回归，它是多重线性回归的一种特殊形式。数据回归处理的计算方法见附录 D。8.2.2 作为回归方法的替代，可采用附录 E 描述的正交多项式法，这一方法尤其适用于事先不知道拟合次数的情形。8.2.3 当 x 不是平均分布时，可采用有限差分方法快速预测适当的数据拟合多项式的次数，并计算多项式的系数，有限差分法可参考 ISO-TR-7066-2 中附录 E，此方法的不确定的计算已超出了本标准的范围。8.3 最优

28、拟合次数的选择 DL/T 1961 2019 11 8.3.1 最优拟合次数的确定原则：进一步增加拟合次数，当拟合结果没有显著改善时的多项式最高次数为最优拟合次数。对于每个拟合次数，应采用下式计算拟合偏差的标准差sr：sr2=(yi y i)2/(n m 1)ni=1（43）式中：y ix=xi时由拟合多项式(式 42)计算的预测值。sr2等价于本标准的sr2(y,x)项。8.3.2 拟合多项式的次数 m 应远小于数据点的个数 n。8.3.3 如果数据可以用次数为 m 的多项式较好的拟合，则当次数达到 m 时，sr会显著减小，此后sr几乎保持不变。若sr的变化不显著，应使用其它显著性检查方法确

29、定最优拟合次数或寻求较为显著的目标来确定最优的多项式次数。8.3.4 多项式次数从m-1增加到m的过程中，如果新系数bm明显不等于零，比如bm+t95s(bm)和bmt95s(bm)(bm的置信概率为 95%)不包括 0，则认为次数 m 的增加显著改善了曲线拟合效果。这一条件可以表述为：bms(bm)t95 式中：t95自由度=n m 1，置信概率为 95%时学生分布的 t 值。t95作为自由度的指数函数，可以根据下面的经验公式计算：t95=1.95+2.36 +3.2 2+5.2 3.84（44）8.3.5 对于正交多项式的系数gm的计算参见附录 E：gms(gm)t95 s2(bm)和s2

30、(gm)的系数方差的计算式参见附录 D 和附录 E。8.3.6 在第一次增加多项式次数对拟合结果没有改善时，应再增加一次多项式的次数以检验多项式次数对拟合结果是否产生明显变化。8.3.7 多项式最高次数对拟合结果的提高达到置信水平 95%时，可认为是最优次数。在选择这一次数作为待拟合数据的最优表达式之前，应考虑曲线的预期性形状，需要拟合的区间，拟合的精度等因素的影响。尽量避免拟合多项式的形式过于复杂。画出数据点和可能的拟合曲线。可以更直观的展示数据点的真实关系，给出合理的多项式形式。8.4 非线性拟合过程中的不确定度计算 8.4.1 在 95%的置信水平下，拟合预测值y 的随机不确定度由下式给

31、定：(y)=t95s(y)（45）上式中，s(y)是y 的方差s2(y)的平方根。附录 D 和附录 E 给出了s2(y)的定义，通常s2(y)可以用最高拟合次数为2m的多项式函数计算。计算s2(y)所选择的数据应足够多，以避免因减法运算引起的误差。8.4.2 对于 y 值，其 95%随机置信限为：y er(y)DL/T 1961 2019 12 校准系数的不确定度为：e(y c)=er2(y)+es2(y)0.5 （46）式中，es2(y)是y 的系统不确定度。8.4.3 如果因变量是已经过转换，上述所有的不确定度计算应根据转换后的参数计算。DL/T 1961 2019 13 附录附录 A(规

32、范性附录规范性附录)一般函数标准差的计算一般函数标准差的计算 A.1 如果总的方差是两个或两个以上构成变量的积或商，则不能采用公式（5），而应采用与通用函数相对应的更为复杂的标准差的表达式。式中的微分项(x xi)与 GB/T27759-2011 中用于合成灵敏度的系数(R Yi)的意义相同。A.2 如果 X=f(x1,x2xn)，其中 f 表示任意函数，则：()()()()()()()()()()()22112212121313112 .+2,.,nnnVarXXxVar xXxVar xXxXxCov x xXxXxCov x xXxXxCov xx=+（A.1）如果所有涉及高阶导数的项均

33、可以忽略且协方差为 0，例如，变量是独立的，则式(A.1)简化为只有第一行。A.3 作为例子，考虑明渠上流量测量表的流量方程：,ic iiiQb d v=（A.2）式中，bc，d 和 v 为自变量，应用式(A.1)的第一行，则有：()()()()()()222,222,=iiiic ic iiiiiiiiic ic iic iiVarQQbVar bQdVar dQvVar vd vVar bb vVar db dVar v=+（A.3）DL/T 1961 2019 14 附录附录 B(资料性附录资料性附录)明渠校准的算例明渠校准的算例 B.1 使用符号 h0 零流量时，级的基准修正，单位为米

34、；h 测量级，单位为米 c 系数 b 指数 Q 流量，单位为立方米每秒 B.2 表 B.1 中数据给定了级-流量的关系，计算校准方程和测量点相对于最佳拟合直线的标准差(sR)和随机不确定度 eR(Q)的关系。B.3 明渠的流量测量通常采用速度-面积法进行校准，水位与流量的关系表达式为：0()bQc hh=+（B.1）将上式改写为对数形式：0lnlnln()Qcbhh=+（B.2）进一步用下式进行替代：0ln()lny lnhhxQca+=；（B.3）则原式简化为如式(1)所示的线性方程：yabx=+（B.4）B.4 对于此类校准，水位的测量误差远小于流量测量误差，式(14)中 大于 20，则可

35、以采用本部分第 7.3.2的最小二乘法进行拟合。将表 B.1 中的数据依次代入式(23)和式(17)，计算得到校准曲线的斜率为：232(2.9337)93.7855(15.5798)32(35.5093)(15.5798)=1.5301b=（B.5）截距为：ln2.9308 1.5301(0.4869)3.6757c=（B.6）因此 ln3.6757 1.5301ln(0.115)Qh=（B.7）或 1.530139.479(0.115)Qh=（B.8）DL/T 1961 2019 15 表 B.1 采用最小二乘法计算水位-流量曲线的典型数据 观测值 编号 Q m3/s 水位（h）m(h+h0

36、)m lnQi(yi)Ln(h+h0)(xi)xiyi xi2 1 2.463 0.272 0.157 0.9014-1.8515-1.6689 3.4280 2 2.352 0.273 0.158 0.8437-1.8452-1.5568 3.4048 3 2.923 0.303 0.188 1,0726-1.6713-1.7926 2.7932 4 3,242 0.307 0.192 1.1762-1.6502-1.9410 2.7232 5 3,841 0.334 0.219 1.3457-1.5187-2.0437 2.3064 6 4.995 0.374 0.259 1.6084-1

37、.3509-2.1728 1.8248 7 5.410 0.393 0.278 1.6882-1.2801-2.1611 1.6386 8 5.422 0.394 0.279 1.6905-1.2765-2.1579 1.6294 9 5.883 0.402 0.287 1.7721-1.2483-2.2121 1.5582 10 6.154 0.410 0.295 1.8171-1.2208-2.2183 1.4904 11 7.376 0.463 0.348 1.9982-1.0556-2.1093 1.1143 12 9.832 0.520 0.405 2.285-0.9039-2.06

38、60 0.8170 13 11.321 0.548 0.433 2.4266-0.8370-2.0311 0.7006 14 12.372 0.576 0.461 2.5154-0.7744-1.9479 0.5997 15 11.825 0.580 0.465 2.4702-0.7657-1.8914 0.5863 16 13.826 0.616 0.501 2.6266-0.6911-1.8152 0.4776 17 14.102 0.626 0.511 2.6463-0.6714-1.7767 0.4508 18 19.020 0.721 0.606 2.9455-0.5009-1.47

39、54 0.2609 19 19.790 0.739 0.624 2.9852-0.4716-1.4078 0.2224 20 20.280 0.747 0.632 3.0096-0.4589-1.3811 0.2106 21 21.204 0.796 0.681 3.0542-0.3842-1.1734 0.1476 22 23.996 0.846 0.731 3.1779-0.3133-0.9956 0.0982 23 36.242 1.041 0.926 3.5902-0.0769-0.2761 0.0059 24 54.591 1.340 1.225 3.9999 0.2029 0.81

40、16 0.0412 25 67.327 1.526 1.411 4.4096 0.3443 1.4494 0.1185 26 79.050 1.761 1.646 4.3701 0.4983 2.1776 0.2483 27 110.783 2.010 1.895 4.7076 0.6392 3.0091 0.4086 28 162.814 2.632 2.517 5.0926 0.9231 4.7010 0.8521 29 227.600 3.265 3.150 5.4276 1.1474 6.2276 1.3165 30 228.800 3.280 3.165 5.4328 1.1524

41、6.2597 1.3276 31 228.500 3.306 3.191 5.4315 1.1603 6.3022 1346.3 32 236.600 3.340 3.225 5.4664 1.1709 6.4006 1.3710 总计 93.7855-15.5798-2.9337 35.5093 注:数据修正 h0=-0.115m B.5 测量点相对于最佳拟合直线的标准方差由式(19)定义。即：DL/T 1961 2019 16 ()0.5lnln(2)RisQQn=（B.9）代入表 B.2 中的数据，得：()0.50.02918 300.031Rs=（B.10）B.6 lnQ的在点(h+h

42、0)k计算得到的随机不确定度可按式(34)计算，即：0.52002000.52ln()ln()1(ln)100ln()-ln()=6.3 0.03125ln(0.115)0.486927.9238kRRkhhhheQtsnhhhhh+=+（B.11）同样地，lnQ的随机不确定度可按式(35)计算，即 0.52002000.52ln()ln()1(ln)1+100ln()-ln()=6.3 1.03125ln(0.115)0.486927.9238kRRkhhhheQtsnhhhhh+=+(B.12)B.7 在观测点(h+h0)k计算流量的eR(lnQ)按式(B.11)计算，将计算结果绘制于水位

43、-流量曲线的两侧，以给出流量计算的系统置信限，最小的宽度应该在ln(h+h0)位置。代入表 B.2 中的中第一组数据得：()0.52(ln)=6.3 0.03125-1.8515+0.486927.9238=1.97%ReQ+相似的，对于第 18 组数据：()0.52(ln)=6.3 0.03125-0.5009+0.486927.9238=1.11%ReQ+对于第 32 组的数据：()0.52(ln)=6.3 0.03125-1.1709+0.486927.9238=2.27%ReQ+这三组数据和其它汇总数据在表 B.2 最后一列给出。DL/T 1961 2019 17 图 B.1 基于表

44、B.1 水位-流量曲线 B.8 未转换流量的非对称误差限可按下式计算 置信限高于 95%时：()100 e1z 置信限低于 95%时：()100 1-ez(B.13)式中，z 为式(B.8)右侧不包括系数 100 的数值。以 NO.1 行的观测数据为例，(概率)高于 95%时的置信限为：()0.0197100 e1100(1.0199 1)1.99%=(概率)低于 95%时的置信限为：()()()0.01970.0197100 1-e100 11100 1-0.9804=1.96%e=表 B.2 计算sR和eRln(Q)所需的数据 观测值 编号(h+h0)ln(h+h0)(xi)(xi-x)Q

45、i lnQi(yi)Q ln(Q)（y）(yi-y)2 eR(lnQ)1 0.157-1.8515 1.8621 2.463 0.9014 2.323 0.8428 0.00342 1.97 2 0.158-1.8452 1.8450 2.325 0.8437 2.345 0.8523 0.00007 1.96 3 0.188 1.6713 1.4028 2.923 1.0726 3.060 1.1184 0.00209 1.80 4 0.192-1.6502 1.3533 3.242 1.1762 3.160 1.1506 0.00065 1.78 5 0.219-1.5187 1.0646

46、 3.841 1.3457 3.865 1.3520 0.00003 1.66 6 0.259-1.3509 0.7465 4.995 1.6084 4.996 1.6086 1.00000 1.52 7 0.278-1.2801 0.6292 5.410 1.6882 5.568 1.7170 0.00083 1.46 8 0.279-1.2765 0.6235 5.422 1.6905 5.598 1.7224 0.00101 1.46 9 0.287-1.2483 0.5797 5.883 1.7721 5.846 1.7658 0.00004 1.44 DL/T 1961 2019 1

47、8 10 0.295-1.2208 0.5386 6.154 1.8171 6.097 1.8078 1.00008 1.42 11 0.348-1.0556 0.3234 7.376 1.9982 7.851 2.0606 0.00389 1.30 12 0.405-0.9039 0.1739 9.832 2.2856 9.902 2.2927 0.00005 1.22 13 0.433-0.8370 0.1226 11.321 2.4266 10.968 2.3950 0.00099 1.19 14 0.461-0.7744 0.0826 12.372 2.5154 12.072 2.49

48、09 0.00060 1.16 15 0.465-0.7657 0.0777 11.825 2.4702 12.233 2.5041 0.00115 1.16 16 0.501-0.6911 0.0417 13.826 2.6266 13.711 2.6182 0.00007 1.14 17 0.511-0.6714 0.0340 14.102 2.6463 14.132 2.6484 0.00000 1.14 18 0.606-0.5009 0.0002 19.020 2.9455 18.345 2.9094 0.00130 1.11 19 0.624-0.4716 0.0002 19.97

49、0 2.9852 19.185 2.9541 0.00096 1.11 20 0.632-0.4589 0.0008 20.280 3.096 19.563 2.9736 0.00129 1.11 21 0.681-0.3842 0.0105 21.204 3.0542 21.931 3.0879 0.00113 1.12 22 0.731-0.3133 0.0301 23.996 3.1779 24.442 3.1963 0.00033 1.13 23 0.926-0.076 0.1681 36.242 3.5902 35.098 3.5581 0.00102 1.22 24 1.225 0

50、.2029 0.4758 54.591 3.9999 53.855 3,9863 0.00018 1.38 25 1,411 0.3443 0.6909 67.327 4.2096 66.859 4.2026 0.00004 1.49 26 1.646 0.4983 0.9706 79.050 4.3701 84.631 4.4383 0.00465 1.62 27 1.895 0.6392 1.2681 110.783 4.7076 104.989 4.6538 0.00288 1.74 28 2.517 0.9231 1.9881 182.814 5.0926 162.095 5.0882