CN102680732B - 坡面薄层水流流速测量方法及系统 - Google Patents
坡面薄层水流流速测量方法及系统 Download PDFInfo
- Publication number
- CN102680732B CN102680732B CN201210165464.9A CN201210165464A CN102680732B CN 102680732 B CN102680732 B CN 102680732B CN 201210165464 A CN201210165464 A CN 201210165464A CN 102680732 B CN102680732 B CN 102680732B
- Authority
- CN
- China
- Prior art keywords
- tau
- flow
- convection current
- dispersion equation
- boundary condition
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种坡面薄层水流流速测量方法及系统,涉及边界获取技术领域。本发明以正弦函数和正态分布函数表达的边界条件实现坡面薄层水流流速测量,规避了因边界条件测量不准确而导致的流速测量精度低的问题,从而提高了坡面薄层水流流速测量的精度。另外,不需要先计算边界条件,再计算坡面薄层水流流速,而是同时计算出水流流速和边界条件,进而提高了流速测量的效率。
Description
技术领域
本发明涉及边界获取技术领域,特别涉及一种坡面薄层水流流速测量方法及系统。
背景技术
电解质脉冲模型测量坡面薄层水流流速的方法,从理论上能够在一定条件下准确预测水流流速。然而,由于设备的限制,输入边界为脉冲函数的假定在实际中很难实现,尤其对于短距离流速测量误差较大。夏卫生等人通过实验验证结果表明,脉冲模型方法在距离电解质注入点0.5m处测量的流速相对误差达到38%。电解质示踪实测边界条件法,是在距离电解质注入点约0.05m处新增一组探针,测量得到电解质浓度信号作为模型的边界条件,改进了原来模型中脉冲边界的假设。该方法分别利用正态分布函数和正弦函数形式模拟溶质输入边界,对理论的发展起到了推动性的作用,较脉冲模型而言,该方法显著提高了流速预测的精度,然而,模型采用0.05m处测量得到的函数近似地作为边界条件,仍不是发生在0m处的真实边界条件。由于该方法对边界条件在时间和距离上反应是非常敏感的,边界条件的不准确仍将导致流速测量的误差。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:如何提高坡面薄层水流流速测量的精度。
(二)技术方案
为解决上述技术问题,本发明提供了一种坡面薄层水流流速测量方法,所述方法包括:
在当前坡面薄层水流的待测截面设置感应探针,在所述感应探针上游的预定位置注入电解质,设所述预定位置和待测截面之间的距离为x;
所述感应探针采集所述待测截面的电导率,并将所述电导率换算为归一化电解质浓度;
利用所述归一化电解质浓度、距离x、及采用正弦函数和正态分布函数来表示的边界条件计算所述当前坡面薄层水流的流速。
优选地,计算所述当前坡面薄层水流的流速具体包括以下步骤:
将所述归一化电解质浓度和距离x分别代入采用正弦函数来表示边界条件的对流弥散方程w和采用正态分布函数来表示边界条件的对流弥散方程m中,对所述对流弥散方程w和对流弥散方程m进行拟合,以获得所述当前坡面薄层水流的流速。
优选地,所述对流弥散方程为:
其中, DH为水动力弥散系数,t为时间,u为所述当前坡面薄层水流的流速,f(τ)为所述边界条件,C1(x,t)为所述归一化电解质浓度。
优选地,所述对流弥散方程w为:
所述对流弥散方程m为:
其中,A、B和D为所述边界条件的待确定系数。
优选地,对所述对流弥散方程w和对流弥散方程m采用最小二乘法进行拟合。
本发明还公开了一种坡面薄层水流流速测量系统,所述系统包括:
探针设置模块,用于在当前坡面薄层水流的待测截面设置感应探针,在所述感应探针上游的预定位置注入电解质,设所述预定位置和待测截面之间的距离为x;
换算模块,用于所述感应探针采集所述待测截面的电导率,并将所述电导率换算为归一化电解质浓度;
流速计算模块,用于利用所述归一化电解质浓度、距离x、及采用正弦函数和正态分布函数来表示的边界条件计算所述当前坡面薄层水流的流速。
优选地,计算所述当前坡面薄层水流的流速具体包括以下步骤:
将所述归一化电解质浓度和距离x分别代入采用正弦函数来表示边界条件的对流弥散方程w和采用正态分布函数来表示边界条件的对流弥散方程m中,对所述对流弥散方程w和对流弥散方程m进行拟合,以获得所述当前坡面薄层水流的流速。
优选地,所述对流弥散方程为:
其中, DH为水动力弥散系数,t为时间,u为所述当前坡面薄层水流的流速,f(τ)为所述边界条件,C1(x,t)为所述归一化电解质浓度。
优选地,所述对流弥散方程w为:
所述对流弥散方程m为:
其中,A、B和D为所述边界条件的待确定系数。
优选地,对所述对流弥散方程w和对流弥散方程m采用最小二乘法进行拟合。
(三)有益效果
本发明以正弦函数和正态分布函数表达的边界条件实现坡面薄层水流流速测量,规避了因边界条件测量不准确而导致的流速测量精度低的问题,从而提高了坡面薄层水流流速测量的精度。另外,不需要先计算边界条件,再计算坡面薄层水流流速,而是同时计算出水流流速和边界条件,进而提高了流速测量的效率。
附图说明
图1是按照本发明一种实施方式的坡面薄层水流流速测量方法的流程图;
图2是按照本发明一种实施方式的坡面薄层水流流速测量系统的结构框图;
图3是按照本发明一种实施例的实验硬件结构图;
图4是在流量为12L min-1且坡度为4°的条件下,实测值和模拟值的归一化浓度示意图;
图5是在流量为12L min-1且坡度为8°的条件下,实测值和模拟值的归一化浓度示意图;
图6是在流量为12L min-1且坡度为12°的条件下,实测值和模拟值的归一化浓度示意图;
图7是在流量为24L min-1且坡度为4°的条件下,实测值和模拟值的归一化浓度示意图;
图8是在流量为24L min-1且坡度为8°的条件下,实测值和模拟值的归一化浓度示意图;
图9是在流量为24L min-1且坡度为12°的条件下,实测值和模拟值的归一化浓度示意图;
图10是在流量为48L min-1且坡度为4°的条件下,实测值和模拟值的归一化浓度示意图;
图11是在流量为48L min-1且坡度为8°的条件下,实测值和模拟值的归一化浓度示意图;
图12是在流量为48L min-1且坡度为12°的条件下,实测值和模拟值的归一化浓度示意图;
图13是边界条件通过正态分布函数和正弦函数表示时对应的流速结果比较图;
图14是边界条件通过正态分布函数和正弦函数表示时对应的确定性系数(即R2)结果比较图;
图15是本发明的方法和实测边界条件法分别计算的流速值的对比图;
图16是本发明的方法和实测边界条件法分别计算的确定性系数的对比图;
图17是本发明的方法和实测边界条件法分别计算的确定性系数的统计均值对比图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
图1是按照本发明一种实施方式的坡面薄层水流流速测量方法的流程图;参照图1,所述方法包括:
在当前坡面薄层水流的待测截面设置感应探针,在所述感应探针上游的预定位置注入电解质,设所述预定位置和待测截面之间的距离为x;
所述感应探针采集所述待测截面的电导率,并将所述电导率换算为归一化电解质浓度;
利用所述归一化电解质浓度、距离x、及采用正弦函数和正态分布函数来表示的边界条件计算所述当前坡面薄层水流的流速。
优选地,计算所述当前坡面薄层水流的流速具体包括以下步骤:
将所述归一化电解质浓度和距离x分别代入边界条件采用正弦函数来表示的对流弥散方程w和边界条件采用正态分布函数来表示的对流弥散方程m中,对所述对流弥散方程w和对流弥散方程m进行拟合(本实施方式中,对所述对流弥散方程w和对流弥散方程m采用最小二乘法进行拟合),以获得所述当前坡面薄层水流的流速。
描述溶质在一维稳态水流中的对流弥散方程表达式为:
式中,h为水流深度;w为水流宽度;C为归一化的溶液浓度;x为所述预定位置和待测截面之间的距离,其单位为m;u为水流流速,其单位为m/s;t为时间,其单位为s;DH为水动力弥散系数,其单位为m2/s。
不考虑降雨及入渗影响条件下,流量为常数,水流流速的变化可以忽略,因此:
整理得:
hw=constant (3)
其中,Q0是流量。
结合公式(1)和公式(3),得到:
将边界条件假设为脉冲函数,初始条件和边界条件由下式给出:
C(x,t)/C0=δ(t=t0) x=0(4a)
C(x,t)=0 x=∞ (4b)
C(x,t)=0 t=0 (4c)
因此,可以计算出上面初边值问题的解析解为:
在实际试验操作中电解质一般采用盐溶液,由于盐溶液注入设备的限制,注入溶液的时间是达不到理论意义上的脉冲函数效果的。如果真实的边界条件函数形式以f(t)表示,相应的输出信号可以用真实输入信号与脉冲响应的卷积来计算,优选地,所述对流弥散方程为:
其中, DH为水动力弥散系数,t为时间,u为所述当前坡面薄层水流的流速,f(τ)为所述边界条件,C1(x,t)为所述归一化电解质浓度。
边界条件仍采用正弦函数和正态分布函数两种形式来表达,分别为:
式中,A、B、D为边界条件的确定系数。
优选地,所述对流弥散方程w为:
所述对流弥散方程m为:
其中,A、B和D为所述边界条件的待确定系数。
图2是按照本发明一种实施方式的坡面薄层水流流速测量系统的结构框图;参照图2,所述系统包括:
探针设置模块,用于在当前坡面薄层水流的待测截面设置感应探针,在所述感应探针上游的预定位置注入电解质,设所述预定位置和待测截面之间的距离为x;
换算模块,用于所述感应探针采集所述待测截面的电导率,并将所述电导率换算为归一化电解质浓度;
流速计算模块,用于利用所述归一化电解质浓度、距离x、及采用正弦函数和正态分布函数来表示的边界条件计算所述当前坡面薄层水流的流速。
优选地,计算所述当前坡面薄层水流的流速具体包括以下步骤:
将所述归一化电解质浓度和距离x分别代入采用正弦函数来表示边界条件的对流弥散方程w和采用正态分布函数来表示边界条件的对流弥散方程m中,对所述对流弥散方程w和对流弥散方程m进行拟合,以获得所述当前坡面薄层水流的流速。
优选地,所述对流弥散方程为:
其中, DH为水动力弥散系数,t为时间,u为所述当前坡面薄层水流的流速,f(τ)为所述边界条件,C1(x,t)为所述归一化电解质浓度。
优选地,所述对流弥散方程w为:
所述对流弥散方程m为:
其中,A、B和D为所述边界条件的待确定系数。
优选地,对所述对流弥散方程w和对流弥散方程m采用最小二乘法进行拟合。
实施例1
参照图3,本实施例中的实验硬件结构由5部分组成,分别为水槽C、电解质脉冲发生器A、水量注入单元(未示出,在图中的箭头方向按一定流量注入水)、电导率探针B(即感应探针)、以及数据采集与存储设备D。试验水槽长4m,宽15cm,高50cm。当水流达到稳定后,从水槽的上部由电解质脉冲发生器(位于所述预定位置)注入高浓度的KCl溶液(盐溶液),感应探针(位于所述待测截面)监测的电导率通过数据采集系统由计算机输出,数据采集频率为每秒100个点。
实验设计了3组流量(Q=12,24,48L min-1),3组坡度(S=4,8,12°)和5组预定位置和待测截面之间的距离(x=0.3,0.6,0.9,1.2,1.5m)的组合。
图4~12是在不同流量和不同坡度,实测值和模拟值的归一化溶度示意图,其中,图中最左侧的纵坐标为实测值和模拟值对应的归一化浓度,图中最右侧的纵坐标为边界条件对应的归一化浓度,计算的水流流速和确定性系数列于表1中。
由图4~12可见,不同实验条件下模型较好地拟合了不同位置测得的电解质浓度变化趋势,随着测量距离的增加,曲线峰值逐渐降低。随着流量和坡度的增加,穿透时间逐渐减少。
表1
在不同距离x测得的流速值变化不显著,而随着流量和坡度的增加,流速呈增加趋势,这一结论是合理的,亦说明了本发明方法的有效性。所有测量位置处模型拟合的确定性系数均大于0.9,表明模型能够较好地拟合不同实验工况条件下的实测结果。
为了便于参数对比,图13和图14分别点绘了分别利用正态函数和正弦函数作为边界条件的模型参数的对比结果。
图13和图14中,水流流速和确定性系数的数据点均紧密围绕1:1的直线,说明两种边界条件形式计算的流速结果是一致的。对两组数据分别利用T检验进行分析,表明利用正态函数和正弦函数两种函数形式作为虚拟边界计算的流速和确定性系数均没有显著差异。这一结果表明正弦函数和正态函数均能够准确地描述真实的边界条件,并进一步证明了本发明方法的合理性。
分别利用本发明的方法和实测边界条件法对实测数据进行拟合,并计算水流流速和确定性系数。图15和图16对比分析了两种方法在不同实验条件和测量位置情况下计算的流速值和计算精度。
由图15和图16可见,不同实验条件下两种方法测量的流速值呈一致性变化规律。大部分空心点(4°)位于1:1线以下,实心点(8°)分布在1:1直线上,而十字标号(12°)位于1:1线上部。表明在水槽坡度较小时,本发明的方法计算的流速值低于实测边界条件法计算的流速值,当坡度到达8°时,两种方法测量的流速值相当,当坡度继续增大时,本发明的方法计算的流速值要高于实测边界条件法计算的流速值。拟合的确定性系数分布在1:1直线附近,统计均值分析结果(图17)表明本发明的流速测量精度要高于实测边界条件。
以上实施方式仅用于说明本发明,而并非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型,因此所有等同的技术方案也属于本发明的范畴,本发明的专利保护范围应由权利要求限定。
Claims (4)
1.一种坡面薄层水流流速测量方法,其特征在于,所述方法包括:
在当前坡面薄层水流的待测截面设置感应探针,在所述感应探针上游的预定位置注入电解质,设所述预定位置和待测截面之间的距离为x;
所述感应探针采集所述待测截面的电导率,并将所述电导率换算为归一化电解质浓度;
利用所述归一化电解质浓度、距离x、及采用正弦函数和正态分布函数来表示的边界条件计算所述当前坡面薄层水流的流速;
其中,计算所述当前坡面薄层水流的流速具体包括以下步骤:
将所述归一化电解质浓度和距离x分别代入采用正弦函数来表示边界条件的对流弥散方程w和采用正态分布函数来表示边界条件的对流弥散方程m中,对所述对流弥散方程w和对流弥散方程m进行拟合,以获得所述当前坡面薄层水流的流速;
其中,所述对流弥散方程为:
其中, DH为水动力弥散系数,t为时间,u为所述当前坡面薄层水流的流速,f(τ)为所述边界条件,C1(x,t)为所述归一化电解质浓度;
其中,所述对流弥散方程w为:
所述对流弥散方程m为:
其中,A、B和D为所述边界条件的待确定系数。
2.如权利要求1所述的方法,其特征在于,对所述对流弥散方程w和对流弥散方程m采用最小二乘法进行拟合。
3.一种坡面薄层水流流速测量系统,其特征在于,所述系统包括:
探针设置模块,用于在当前坡面薄层水流的待测截面设置感应探针,在所述感应探针上游的预定位置注入电解质,设所述预定位置和待测截面之间的距离为x;
换算模块,用于所述感应探针采集所述待测截面的电导率,并将所述电导率换算为归一化电解质浓度;
流速计算模块,用于利用所述归一化电解质浓度、距离x、及采用正弦函数和正态分布函数来表示的边界条件计算所述当前坡面薄层水流的流速;
其中,计算所述当前坡面薄层水流的流速具体包括以下步骤:
将所述归一化电解质浓度和距离x分别代入采用正弦函数来表示边界条件的对流弥散方程w和采用正态分布函数来表示边界条件的对流弥散方程m中,对所述对流弥散方程w和对流弥散方程m进行拟合,以获得所述当前坡面薄层水流的流速;
其中,所述对流弥散方程为:
其中, DH为水动力弥散系数,t为时间,u为所述当前坡面薄层水流的流速,f(τ)为所述边界条件,C1(x,t)为所述归一化电解质浓度;
其中,所述对流弥散方程w为:
所述对流弥散方程m为:
其中,A、B和D为所述边界条件的待确定系数。
4.如权利要求3所述的系统,其特征在于,对所述对流弥散方程w和对流弥散方程m采用最小二乘法进行拟合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210165464.9A CN102680732B (zh) | 2012-05-24 | 2012-05-24 | 坡面薄层水流流速测量方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210165464.9A CN102680732B (zh) | 2012-05-24 | 2012-05-24 | 坡面薄层水流流速测量方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102680732A CN102680732A (zh) | 2012-09-19 |
CN102680732B true CN102680732B (zh) | 2014-04-16 |
Family
ID=46812973
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210165464.9A Expired - Fee Related CN102680732B (zh) | 2012-05-24 | 2012-05-24 | 坡面薄层水流流速测量方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102680732B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104535794B (zh) * | 2015-01-16 | 2018-03-20 | 中国水利水电科学研究院 | 一种薄层水流流速测量系统和方法 |
CN105911306B (zh) * | 2016-05-10 | 2019-01-11 | 河海大学 | 流速自动化测量装置及测量方法 |
CN107290129B (zh) * | 2017-05-04 | 2019-09-27 | 中国水利水电科学研究院 | 一种坡面水文试验流量场观测系统及方法 |
CN108759949B (zh) * | 2018-05-30 | 2020-08-11 | 河海大学 | 测量浅水区地下水出流量的装置及其操作方法 |
CN110441546B (zh) * | 2019-09-04 | 2021-03-26 | 中国电建集团中南勘测设计研究院有限公司 | 一种泄水建筑物水流流速测量装置及方法 |
CN110823194A (zh) * | 2019-11-18 | 2020-02-21 | 中国农业大学 | 一种坡面流水深测量方法及装置 |
CN111366747B (zh) * | 2020-03-30 | 2021-02-05 | 中国农业大学 | 一种提高电解质示踪测量侵蚀坡面水流流速的方法 |
CN112462088B (zh) * | 2020-11-18 | 2022-02-25 | 中国农业大学 | 一种盐热耦合测量坡面水流流速的方法 |
CN112633554B (zh) * | 2020-11-27 | 2023-07-25 | 北京林业大学 | 坡面薄层水流流速修正系数预测方法及装置 |
CN113009177B (zh) * | 2021-02-09 | 2022-07-05 | 中国农业大学 | 一种电解质质心法测量坡面水流流速的方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3906998A1 (de) * | 1989-03-04 | 1990-09-06 | Horst Prof Dr Frankenberger | Verfahren zur stroemungsmessung von waessrigen, elektrolytischen fluessigkeiten und vorrichtung hierzu |
CN1219207C (zh) * | 2001-12-18 | 2005-09-14 | 中国科学院水利部水土保持研究所 | 电解质脉冲法测量水流流速和泥沙含量的方法 |
CN101825646B (zh) * | 2010-05-07 | 2012-11-07 | 中国农业大学 | 电解质示踪测量薄层水流流速的方法 |
-
2012
- 2012-05-24 CN CN201210165464.9A patent/CN102680732B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102680732A (zh) | 2012-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102680732B (zh) | 坡面薄层水流流速测量方法及系统 | |
Zhang et al. | A new vertical coordinate system for a 3D unstructured-grid model | |
CN101825646B (zh) | 电解质示踪测量薄层水流流速的方法 | |
CN103487234A (zh) | 基于多尺度排列熵的气液两相流型动力学表征及识别方法 | |
CN103185808A (zh) | 一种基于光电技术的气泡测速系统 | |
CN109185211A (zh) | 一种基于计算流体力学瞬时计算的泵站机组压力脉动预测方法 | |
CN105547957A (zh) | 土壤下渗湿润峰吸力和降雨初损量计算方法 | |
CN106443614B (zh) | 高超声速目标测加速度方法 | |
CN103245765B (zh) | 一种海洋参数测定系统 | |
CN111089627A (zh) | 基于面积流速法的多分辨率明渠自动测流装置 | |
CN103575979B (zh) | 一种数字化测量交流电频率的方法 | |
CN104036120A (zh) | 一种测定材料与构件疲劳s-n曲线性能的单点法 | |
Wang et al. | Three-dimensional modelling of sediment transport under tidal bores in the Qiantang estuary, China | |
CN107270983A (zh) | 基于瞬时流量的数据搭建超声波流量模型的方法 | |
CN219368806U (zh) | 一种基于超声波时差法的多通道明渠流速测量装置 | |
Xu et al. | Effects of large-scale embankments on the hydrodynamics and salinity in the Oujiang River Estuary, China | |
Sakphrom et al. | Simplified stream discharge estimation for hydrological application based on NB-IoT deployment | |
CN104570104A (zh) | 一种基于两步法avf的纵横波地震品质因子提取方法 | |
CN105512492B (zh) | 潮汐流能发电机输出功率的概率建模方法 | |
CN104268334A (zh) | 一种污染物空间展布的验证方法 | |
CN210464571U (zh) | 渠道流量测量设备 | |
CN211373724U (zh) | 基于面积流速法的多分辨率明渠自动测流装置 | |
CN204666303U (zh) | 一种新型高温高压升泡仪 | |
CN105095669A (zh) | 水流冲刷条件下坡面地表糙率及入渗参数的估算方法 | |
CN103076080B (zh) | 基于局部压力测量二维流路内t-s波的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140416 Termination date: 20210524 |