CN113569202A - 统一基准订正的河流流量全量程测量不确定度计算方法 - Google Patents

统一基准订正的河流流量全量程测量不确定度计算方法 Download PDF

Info

Publication number
CN113569202A
CN113569202A CN202111132322.8A CN202111132322A CN113569202A CN 113569202 A CN113569202 A CN 113569202A CN 202111132322 A CN202111132322 A CN 202111132322A CN 113569202 A CN113569202 A CN 113569202A
Authority
CN
China
Prior art keywords
flow
water level
full
curve
measured
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111132322.8A
Other languages
English (en)
Other versions
CN113569202B (zh
Inventor
熊明
吴琼
朱子园
许全喜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Bureau of Hydrology Changjiang Water Resources Commission
Original Assignee
Bureau of Hydrology Changjiang Water Resources Commission
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Bureau of Hydrology Changjiang Water Resources Commission filed Critical Bureau of Hydrology Changjiang Water Resources Commission
Priority to CN202111132322.8A priority Critical patent/CN113569202B/zh
Publication of CN113569202A publication Critical patent/CN113569202A/zh
Application granted granted Critical
Publication of CN113569202B publication Critical patent/CN113569202B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/002—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow wherein the flow is in an open channel

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Remote Sensing (AREA)
  • Mathematical Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Hydrology & Water Resources (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Fluid Mechanics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明提供了本发明提供了统一基准订正的河流流量全量程测量不确定度计算方法，包括全量程水位流量关系曲线及全量程水位流量关系曲线上包络线和全量程水位流量关系曲线下包络线的拟定、统一基准确定、实测流量同基准相对误差订正计算、流量标准差及总不确定度计算；应用本发明可避免因定位、水深、流速等测量仪器或方法的初始常误差对高、中、低水流量相对误差计算差异的影响；可为统一高、中、低水流量测验允许误差标准创造条件；可解决生态流量、水资源监测等中小流量测验误差评价的难题。

Description

统一基准订正的河流流量全量程测量不确定度计算方法
技术领域
本发明涉及水文测量应用技术领域，尤其涉及统一基准订正的河流流量全量程测量不确定度计算方法。
背景技术
我国河流流量测验精度采用总不确定度表达，各类水文站每年按高、中、低水计算一次总不确定度用于精度评定，并与不同水位级的允许误差比较判断是否满足要求。现有的单次流量的相对误差按
Figure 918514DEST_PATH_IMAGE001
计算，其假定各级实测流量误差为相同的正态分布，并据此计算测点标准差和总随机不确定度。由于未考虑起点距、水深、流速等测量仪器初始常误差影响，导致该类初始误差对低水相对误差计算影响较大、对高水影响较小的问题长期存在，不能制定高、中、低水统一的流量测验允许误差，即使《河流流量测验规范》等国家标准，也不得不按水位级以及资料的用途或服务对象分别确定不同水位级的流量允许误差。
如何在计算实测流量相对误差时，消除或弱化各类初始常误差的影响，使所有高、中、低水的流量误差计算保持统一的基准，最终形成统一的流量允许误差标准，是当前水文测验领域急需解决的问题。
发明内容
本发明的目的在于针对上述现有技术的不足，提供了统一基准订正的河流流量全量程测量不确定度计算方法，可解决生态流量、水资源监测的中小流量测验误差评价的难题。
为实现上述目的，本发明采用了如下技术方案：
本发明提供了统一基准订正的河流流量全量程测量不确定度计算方法，包括以下步骤：
S1、全量程水位流量关系曲线及全量程水位流量关系曲线上包络线和全量程水位流量关系曲线下包络线的拟定：根据需要整编时间段内的所有实测水位、流量点据，拟定该整编时间段的水位流量关系及拟定该整编时间段的水位流量关系的上包络线和拟定该整编时间段的水位流量关系的下包络线；
S2、统一基准确定：根据所述实测水位、所述流量点据分布，选择包括最大、最小流量中值或平均值
Figure 510033DEST_PATH_IMAGE002
及所对应的水位流量关系曲线水位
Figure 876423DEST_PATH_IMAGE003
，最高、最低水位中值或平均值
Figure 462738DEST_PATH_IMAGE003
及所对应的水位流量关系曲线流量
Figure 242475DEST_PATH_IMAGE002
作为基准点，并查读
Figure 309788DEST_PATH_IMAGE003
对应的上包络线值
Figure 593002DEST_PATH_IMAGE004
、
Figure 353147DEST_PATH_IMAGE005
对应的下包络线值
Figure 557864DEST_PATH_IMAGE006
；
S3、实测流量同基准相对误差订正计算：第
Figure 225605DEST_PATH_IMAGE007
个实测流量点
Figure 832167DEST_PATH_IMAGE008
的相对误差计算按推求的统一基准订正公式：
Figure 763214DEST_PATH_IMAGE009
其中，
Figure 475735DEST_PATH_IMAGE010
为第
Figure 150429DEST_PATH_IMAGE011
个实测流量点；
Figure 611498DEST_PATH_IMAGE012
为第
Figure 510184DEST_PATH_IMAGE013
个实测流量点相应水位
Figure 689492DEST_PATH_IMAGE014
在水位流量关系上的插值；
Figure 902299DEST_PATH_IMAGE015
第
Figure 217874DEST_PATH_IMAGE013
个实测流量点相应水位
Figure 553040DEST_PATH_IMAGE016
在水位流量下包络线的插值；
Figure 688486DEST_PATH_IMAGE017
第
Figure 767301DEST_PATH_IMAGE013
个实测流量点相应水位
Figure 934452DEST_PATH_IMAGE016
在水位流量上包络线的插值；
S4、流量标准差及总不确定度计算：根据所有所述实测流量点计算的统一基准订正的相对误差，计算全量程的流量测点标准差
Figure 174941DEST_PATH_IMAGE018
，并由
Figure 63263DEST_PATH_IMAGE019
计算各种置信水平的总随机不确定度，用于对流量测验精度的总体评价；
Figure 883451DEST_PATH_IMAGE020
为采用的实测流量点的个数。
进一步，所述整编时间段包括年、月或给定的固定时间段。
进一步，所述水位流量关系包括稳定的水位流量关系、临时曲线或单值化处理的单一线。
进一步，所述水位流量关系拟定采用包括目估适线法、基于最小二乘或最小绝对值原理的最优化方法。
进一步，所述拟定该整编时间段的水位流量关系的上包路线和所述拟定该整编时间段的水位流量关系的下包路线指位于所有水位、流量实测点据之上和所述位、流量实测点据之下的直线或曲线，采用包括目估适线法、流量独立分量置信水平90%随机不确定度确定。
进一步，所述统一基准订正系指将所有用于计算相对误差的实测流量点据，统一换算成某一标准流量位置的相对误差值。
进一步，所述全量程系指包括所有流量测点且涵盖最高最低水位。
本发明的有益效果为：可避免因定位、水深、流速等测量仪器或方法的初始常误差对高、中、低水流量相对误差计算差异的影响；可为统一高、中、低水流量测验允许误差标准创造条件；可解决生态流量、水资源监测等中小流量测验误差评价的难题。
附图说明
图1 为本发明统一基准订正的河流流量全量程测量不确定度计算方法的原理图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，下面结合附图，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
请参阅1，统一基准订正的河流流量全量程测量不确定度计算方法，其特征在于，包括以下步骤：
S1、全量程水位流量关系曲线及全量程水位流量关系曲线上包络线和全量程水位流量关系曲线下包络线的拟定：根据需要整编时间段内的所有实测水位、流量点据，拟定该整编时间段的水位流量关系及拟定该整编时间段的水位流量关系的上包络线和拟定该整编时间段的水位流量关系的下包络线；
S2、统一基准确定：根据所述实测水位、所述流量点据分布，选择包括最大、最小流量中值或平均值
Figure 970356DEST_PATH_IMAGE002
及所对应的水位流量关系曲线水位
Figure 319429DEST_PATH_IMAGE003
，最高、最低水位中值或平均值
Figure 163888DEST_PATH_IMAGE003
及所对应的水位流量关系曲线流量
Figure 584505DEST_PATH_IMAGE002
作为基准点，并查读
Figure 260337DEST_PATH_IMAGE003
对应的上包络线值
Figure 42961DEST_PATH_IMAGE021
、
Figure 905874DEST_PATH_IMAGE022
对应的下包络线值
Figure 864603DEST_PATH_IMAGE006
；
S3、实测流量同基准相对误差订正计算：第
Figure 863783DEST_PATH_IMAGE007
个实测流量点
Figure 882555DEST_PATH_IMAGE023
的相对误差计算按推求的统一基准订正公式：
Figure 701606DEST_PATH_IMAGE024
其中，
Figure 932867DEST_PATH_IMAGE025
为第
Figure 786554DEST_PATH_IMAGE011
个实测流量点；
Figure 976227DEST_PATH_IMAGE012
为第
Figure 548153DEST_PATH_IMAGE013
个实测流量点相应水位
Figure 845755DEST_PATH_IMAGE014
在水位流量关系上的插值；
Figure 553948DEST_PATH_IMAGE026
第
Figure 914522DEST_PATH_IMAGE013
个实测流量点相应水位
Figure 973745DEST_PATH_IMAGE016
在水位流量下包络线的插值；
Figure 15651DEST_PATH_IMAGE027
第
Figure 906246DEST_PATH_IMAGE013
个实测流量点相应水位
Figure 375405DEST_PATH_IMAGE016
在水位流量上包络线的插值；
S4、流量标准差及总不确定度计算：根据所有所述实测流量点计算的统一基准订正的相对误差，计算全量程的流量测点标准差
Figure 718661DEST_PATH_IMAGE028
，并由
Figure 564258DEST_PATH_IMAGE019
计算各种置信水平的总随机不确定度，用于对流量测验精度的总体评价；
Figure 247043DEST_PATH_IMAGE029
为采用的实测流量点的个数。
所述整编时间段包括年、月或给定的固定时间段。
所述水位流量关系包括稳定的水位流量关系、临时曲线或单值化处理的单一线。
所述水位流量关系拟定采用包括目估适线法、基于最小二乘或最小绝对值原理的最优化方法。
所述拟定该整编时间段的水位流量关系的上包路线和所述拟定该整编时间段的水位流量关系的下包路线指位于所有水位或流量实测点据之上和所述水位或流量实测点据之下的直线或曲线，采用包括目估适线法或流量独立分量置信水平90%随机不确定度确定。
所述统一基准订正系指将所有用于计算相对误差的实测流量点据，统一换算成某一标准流量位置的相对误差值。
所述全量程系指包括所有流量测点且涵盖最高最低水位。
以上所述实施例仅表达了本发明的实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (7)

1.统一基准订正的河流流量全量程测量不确定度计算方法，其特征在于，包括以下步骤：
S1、全量程水位流量关系曲线及全量程水位流量关系曲线上包络线和全量程水位流量关系曲线下包络线的拟定：根据需要整编时间段内的所有实测水位、流量点据，拟定该整编时间段的水位流量关系及拟定该整编时间段的水位流量关系的上包络线和拟定该整编时间段的水位流量关系的下包络线；
S2、统一基准确定：根据所述实测水位、所述流量点据分布，选择包括最大、最小流量中值或平均值
Figure 832562DEST_PATH_IMAGE001
及所对应的水位流量关系曲线水位
Figure 952964DEST_PATH_IMAGE002
，最高、最低水位中值或平均值
Figure 404805DEST_PATH_IMAGE002
及所对应的水位流量关系曲线流量
Figure 596752DEST_PATH_IMAGE001
作为基准点，并查读
Figure 348808DEST_PATH_IMAGE002
对应的上包络线值
Figure 577795DEST_PATH_IMAGE003
、
Figure 782511DEST_PATH_IMAGE004
对应的下包络线值
Figure 916165DEST_PATH_IMAGE005
；
S3、实测流量同基准相对误差订正计算：第
Figure 522727DEST_PATH_IMAGE006
个实测流量点
Figure 391457DEST_PATH_IMAGE007
的相对误差计算按推求的统一基准订正公式：
Figure 817890DEST_PATH_IMAGE008
其中，
Figure 492585DEST_PATH_IMAGE009
为第
Figure 953653DEST_PATH_IMAGE010
个实测流量点；
Figure 180235DEST_PATH_IMAGE011
为第
Figure 359544DEST_PATH_IMAGE012
个实测流量点相应水位
Figure 327279DEST_PATH_IMAGE013
在水位流量关系上的插值；
Figure 908433DEST_PATH_IMAGE014
第
Figure 181283DEST_PATH_IMAGE012
个实测流量点相应水位
Figure 582308DEST_PATH_IMAGE015
在水位流量下包络线的插值；
Figure 333226DEST_PATH_IMAGE016
第
Figure 627942DEST_PATH_IMAGE012
个实测流量点相应水位
Figure 71692DEST_PATH_IMAGE015
在水位流量上包络线的插值；
S4、流量标准差及总不确定度计算：根据所有所述实测流量点计算的统一基准订正的相对误差，计算全量程的流量测点标准差
Figure 960014DEST_PATH_IMAGE017
，并由
Figure 514623DEST_PATH_IMAGE018
计算各种置信水平的总随机不确定度，用于对流量测验精度的总体评价；其中，
Figure 801861DEST_PATH_IMAGE020
为采用的实测流量点的个数。
2.根据权利要求1所述的统一基准订正的河流流量全量程测量不确定度计算方法，其特征在于：所述整编时间段包括年、月或给定的固定时间段。
3.根据权利要求1所述的统一基准订正的河流流量全量程测量不确定度计算方法，其特征在于：所述水位流量关系包括稳定的水位流量关系、临时曲线或单值化处理的单一线。
4.根据权利要求1所述的统一基准订正的河流流量全量程测量不确定度计算方法，其特征在于：所述水位流量关系拟定采用包括目估适线法、基于最小二乘或最小绝对值原理的最优化方法。
5.根据权利要求1所述的统一基准订正的河流流量全量程测量不确定度计算方法，其特征在于：所述拟定该整编时间段的水位流量关系的上包络线和所述拟定该整编时间段的水位流量关系的下包络线指位于所有水位或流量实测点据之上和所述水位或流量实测点据之下的直线或曲线，采用包括目估适线法或流量独立分量置信水平90%随机不确定度确定。
6.根据权利要求1所述的统一基准订正的河流流量全量程测量不确定度计算方法，其特征在于：所述统一基准订正系指将所有用于计算相对误差的实测流量点据，统一换算成某一标准流量位置的相对误差值。
7.根据权利要求1所述的统一基准订正的河流流量全量程测量不确定度计算方法，其特征在于：所述全量程系指包括所有流量测点且涵盖最高最低水位。
CN202111132322.8A 2021-09-27 2021-09-27 统一基准订正的河流流量全量程测量不确定度计算方法 Active CN113569202B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111132322.8A CN113569202B (zh) 2021-09-27 2021-09-27 统一基准订正的河流流量全量程测量不确定度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111132322.8A CN113569202B (zh) 2021-09-27 2021-09-27 统一基准订正的河流流量全量程测量不确定度计算方法

Publications (2)

Publication Number Publication Date
CN113569202A true CN113569202A (zh) 2021-10-29
CN113569202B CN113569202B (zh) 2022-01-07

Family

ID=78174658

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111132322.8A Active CN113569202B (zh) 2021-09-27 2021-09-27 统一基准订正的河流流量全量程测量不确定度计算方法

Country Status (1)

Country Link
CN (1) CN113569202B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080133194A1 (en) * 2006-10-30 2008-06-05 Schlumberger Technology Corporation System and method for performing oilfield simulation operations
CN102967334A (zh) * 2012-09-26 2013-03-13 朱作行 利用对信号包络线处理测量流体流量的系统及方法
CN109635435A (zh) * 2018-12-12 2019-04-16 中山大学 一种基于贝叶斯理论的天然河道水位流量关系确定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080133194A1 (en) * 2006-10-30 2008-06-05 Schlumberger Technology Corporation System and method for performing oilfield simulation operations
CN102967334A (zh) * 2012-09-26 2013-03-13 朱作行 利用对信号包络线处理测量流体流量的系统及方法
CN109635435A (zh) * 2018-12-12 2019-04-16 中山大学 一种基于贝叶斯理论的天然河道水位流量关系确定方法

Also Published As

Publication number Publication date
CN113569202B (zh) 2022-01-07

Similar Documents

Publication Publication Date Title
CN112197749B (zh) 一种波浪浮标有效波高的交叉校准方法及装置
CN114047472B (zh) 一种智能电表的计量误差监测系统及其监测方法、装置
CN110895354A (zh) 一种基于泰森多边形动态调整的面雨量计算方法
CN115453662B (zh) 一种时间维度和空间维度结合的异常站点筛查方法
CN108779995B (zh) 用于检测环境参数的传感器和用于校准这种传感器的方法
CN111998918A (zh) 一种误差校正方法、误差校正装置及流量传感系统
CN111429365A (zh) 一种海洋内孤立波的振幅反演装置及方法
CN105572191A (zh) 一种电化学气体传感器的压力补偿方法
CN110792928A (zh) 一种基于大数据的管道泄漏诊断组合算法
CN115854999A (zh) 基于场景自适应的h-adcp断面平均流速自校正方法
CN101319940B (zh) 一种长杆铂电阻温度计的检测装置及其使用方法
CN115876288A (zh) 一种基于大数据的电子仪表故障分析方法及系统
CN113531403B (zh) 水管泄漏检测方法及装置
CN105204047B (zh) 一种卫星导航系统中观测量单个粗差的探测与修复方法
CN113569202B (zh) 统一基准订正的河流流量全量程测量不确定度计算方法
CN114777887A (zh) 一种超声波水表的零点漂移消除方法
CN114279415A (zh) 一种基于卫星测高数据计算和校准河湖水位的方法
CN110553631A (zh) 一种关于水位流量关系的水位测量系列误差分析方法
CN104061985A (zh) 一种河道压力式水位计的检定与校准方法
CN108562616A (zh) 一种天然气管道外检测方法
CN115545112B (zh) 一种大量地下水实时自动监测数据自动识别和处理的方法
CN107655539B (zh) 水位监测方法与系统
CA2661081C (en) Method for the calculation of fluid interface level
CN113157684B (zh) 一种水利海量数据的查错方法
CN102645186B (zh) 一种基于双线性插值算法的水膜厚度传感器实现方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant