CN111366747B - 一种提高电解质示踪测量侵蚀坡面水流流速的方法 - Google Patents

一种提高电解质示踪测量侵蚀坡面水流流速的方法 Download PDF

Info

Publication number
CN111366747B
CN111366747B CN202010236663.9A CN202010236663A CN111366747B CN 111366747 B CN111366747 B CN 111366747B CN 202010236663 A CN202010236663 A CN 202010236663A CN 111366747 B CN111366747 B CN 111366747B
Authority
CN
China
Prior art keywords
data
electrolyte
probe
water flow
electrolyte concentration
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
Application number
CN202010236663.9A
Other languages
English (en)
Other versions
CN111366747A (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.)
China Agricultural University
Original Assignee
China Agricultural University
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 China Agricultural University filed Critical China Agricultural University
Priority to CN202010236663.9A priority Critical patent/CN111366747B/zh
Publication of CN111366747A publication Critical patent/CN111366747A/zh
Application granted granted Critical
Publication of CN111366747B publication Critical patent/CN111366747B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/18Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance
    • G01P5/20Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance using particles entrained by a fluid stream

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

本发明属于地表水文过程技术领域。本发明涉及一种提高电解质示踪测量侵蚀坡面水流流速的方法,所述方法包括如下步骤:1)电信号数据采集;2)拖尾数据去除;3)去除拖尾数据后归一化和流速计算;采用拖尾数据去除方法,在易侵蚀坡面可以准确地测量电解质吸附条件下的薄层水流流速。该方法方便实用,工作可靠。

Description

一种提高电解质示踪测量侵蚀坡面水流流速的方法
技术领域
本发明属于地表水文过程技术领域。本发明涉及一种电解质示踪测量薄层水流流速的改进方法,尤其是涉及床面吸附和阻滞引起的电解质浓度分布拖尾现象,造成溶质运移模型无法求解的问题。
背景技术
薄层水流是指毫米级深度的坡面水流,是造成坡面水土流失的主要营力,也是输移泥沙和面源污染物的媒介。薄层水流的水动力特征及流速是认识坡面侵蚀动力过程、分析水流输沙能力和剥蚀能力的重要参数。确定坡面水流流速是建立水土流失量预报模型的关键。因此,在含沙水流条件下,改进以溶质运移理论为基础的电解质示踪技术,快速准确地测量薄层水流流速,具有重要的理论和实践意义。
静止的水中,电解质的运动仅由分子扩散作用控制,某一位置的电导率随时间的分布可用正态分布描述。同理,由于对流作用,水流在通过某一固定位置的电解质浓度(电导率)在水流方向上也符合正太分布,这个过程仅受水动力弥散系数的影响。
人工草皮、砾石、冻结土壤等侵蚀较弱的坡面,电解质示踪剂受水流中泥沙颗粒的吸附和扰动影响很小,同时坡面多孔介质对电解质吸附滞留影响也较小。在上述条件下,电解质运移满足溶质运移控制方程的理想边界条件,因此,测量过程中探针测定的电导率数据随时间的分布仅受水动力弥散作用的影响,左偏(拖尾)很小,采用已有的溶质运移模型可准确计算水流流速。然而,在易侵蚀坡面,如解冻土坡面、新耕农地坡面、陡坡坡面等,薄层水流侵蚀搬运土壤颗粒,形成高含沙水流,侵蚀造成坡面床底糙度和孔隙增加,同时由于电解质溶液密度大于清水密度,床面会吸附和滞留部分电解质。
在电导率探针布设点,含电解质示踪剂的水流经过探针,电导率从零值逐渐增加到峰值,此时间段为上升区;此后电导率从峰值逐渐降低到零值,此时间段为下降区。电解质被沉积泥沙和坡面孔隙吸附阻滞,且被吸附的电解质的量较电解质总量小。
在上升区,水流中总电解质浓度处上升阶段,水流电导率总大于被吸附电解质对应的电导率;故在上升区,电解质的吸附和阻滞对水流电导率随时间的分布影响较小,电导率可仍可准确反映溶质运移过程。坡面水流流速比床底孔隙中的流速大2-3个数量级。因此,在下降区,水流电导率从峰值快速降低至被阻滞和吸附在泥沙中的电解质对应的床面水流电导率值。此后,由于床面吸附电解质的释放,水流电导率缓慢降低并形成拖尾,此时,探针测定值无法真实地反映电解质在水流中的传导过程。
上述由于床面吸附和阻滞引起的电解质浓度分布拖尾现象,无法真实地反映示踪剂溶质在水流中由于对流和水动力弥散作用下的电解质浓度分布;采用含吸附和阻滞造成的拖尾误差的电导率分布数据计算水流流速,会造成溶质运移方程无法求解,或求解结果错误。
测定溶质在薄层水流中电导率变化,用溶质运移模型拟合电解质输移过程,从而计算出水流流速等模型参数。溶质运移控制方程如下:
Figure BDA0002431228910000021
边界条件如下:
C(x,t)=C0δ(t) x=0 方程2a
C(x,t)=0 x=∞ 方程2b
C(x,t)=0 t=0 方程2c
方程1存在如下解析解:
Figure BDA0002431228910000031
其中,C是正则化后的电解质浓度,C0是拟合得到的电解质初始浓度;x是探针布设点到马氏瓶的电解质溶液的入水口的坡长,单位为m;t是时间,单位为s;DH是水动力弥散系数,单位为m2 s-1;u是流速,单位是m s-1
将探针测量电解质浓度分布数据带入方程3,用最小二乘法,拟合参数C0,DH和u。u即为测量的坡面水流流速。
目前电解质示踪法测定薄层水流流速被国内外广泛应用,但都集中在非侵蚀坡面,较少关注侵蚀坡面,也缺乏对电解质被侵蚀泥沙吸附影响下的溶质运移过程模拟的认识,对减少电解质吸附影响的方法未能展开研究。
发明内容
本发明提供一种提高电解质示踪测量侵蚀坡面水流流速的方法,其降低电解质被吸附和阻滞影响,改进采用溶质运移示踪技术测量薄层水流流速。
本发明的目的是通过以下技术方案实现的:
一种提高电解质示踪测量侵蚀坡面水流流速的方法,所述方法包括如下步骤:
1)电信号数据采集
在侵蚀坡面上游设置一个马氏瓶6,马氏瓶6内装有高溶解度电解质的饱和溶液;
所述马氏瓶6出口设有电磁阀1,电磁阀1的出口距离坡面有一定铅垂距离;电磁阀1的开启时长为固定值,确保马氏瓶6释放确定体积电解质溶液;
一个第一探针2布设在马氏瓶6的下游的侵蚀坡面,第一探针2测定水流电导率,确定溶质运移模型脉冲边界条件;
在第一探针2的下游的侵蚀坡面布设一个第二探针3,第二探针3测定集中水流经过该点的电导率;
第一探针2、第二探针3和电磁阀1均连接至数据采集器4,数据采集器4和操作控制部件5连接;
集中水流稳定后,电磁阀1开启,释放电解质溶液,同时第一探针2和第二探针3开始测定水流电导率实时变化,测定溶质运移过程;电导率实时数据存储在数据采集器4;
2)拖尾数据去除
首先根据电解质浓度和电导率的函数关系,将第一探针2和第二探针3测定的电导率数据转换为电解质浓度数据,再对电解质浓度数据归一化;
对电解质浓度数据进行归一化处理,即各数据点的电解质浓度ci均除以本次测量过程中所有数据点的电解质浓度最大值cmax
Ci=ci/cmax
其中,i对应测量时间,ci为各数据点的电解质浓度,cmax为本次测量过程中所有数据点的电解质浓度的最大值;
去除拖尾数据的方法如下:
首先,基于对称性原则,根据归一化后的电解质浓度上升段数据预测归一化后的电解质浓度下降段测量数据,并将相同时刻预测数据与测量数据比较;如果在某时刻开始,测量值均比预测值大一个预定的比例以上,则认为从该时刻以后的测量数据为拖尾数据;该预定的比例值可根据测量工况不同调整;
3)去除拖尾数据后归一化和流速计算
去除拖尾数据点后,对剩下的数据进行归一化处理,即各数据点的电解质浓度ci均除以本次测量过程中所有数据点的电解质浓度最大值cmax
Ci=ci/cmax
其中,i对应测量时间,ci为各数据点的电解质浓度,cmax为本次测量过程中所有数据点的电解质浓度的最大值;
采用归一化后的电解质浓度Ci及方程3,采用参数拟合,即得流速、初始浓度、弥散系数:
Figure BDA0002431228910000051
其中,C是归一化后后的电解质浓度,C0是拟合得到的电解质初始浓度;x是探针布设点到马氏瓶6的电解质溶液的入水口的坡长,单位为m;t是时间,单位为s;DH是水动力弥散系数,单位为m2 s-1;u是流速,单位是m s-1
将探针测量电解质浓度分布数据带入方程3,用最小二乘法,拟合参数C0、DH和u;
u即为测量的坡面水流流速。
步骤1中,所述高溶解度电解质的饱和溶液为氯化钾饱和溶液。
步骤1中,电磁阀1的出口距离坡面的铅垂距离为10cm。
步骤1中,第一探针2与马氏瓶6的电解质溶液的入水口之间的距离为10cm。
步骤1中,根据测量条件和精度要求,第二探针3能够布置多个,能够得到多个探针位置的流速数据。
步骤2中,所述预定的比例值为30%。
本发明的有益效果在于:
采用拖尾数据去除方法,在易侵蚀坡面可以准确地测量电解质吸附条件下的薄层水流流速。该方法方便实用,工作可靠。
附图说明
图1本发明的设备布置和电信号数据采集步骤示意图;
图2为去拖尾数据电解质溶质运移过程模拟示意图;
图3为含拖尾数据的电解质浓度随时间变化过程;
图4为去除拖尾数据的电解质浓度随时间变化过程;
图5有拖尾的电解质溶质运移过程模拟;
图6去除拖尾数据的电解质溶质运移过程模拟。
附图标记:
1、电磁阀
2、第一探针
3、第二探针
4、数据采集器
5、操作控制部件
6、马氏瓶
具体实施方式
下面结合附图和实施例对本发明进行进一步说明。
本发明的一种提高电解质示踪测量侵蚀坡面水流流速的方法,包括如下步骤:
1、电信号数据采集
如图1所示,在侵蚀坡面上游设置一个马氏瓶6,马氏瓶6内装有高溶解度电解质的饱和溶液。优选地,所述高溶解度电解质的饱和溶液为氯化钾饱和溶液。
所述马氏瓶6出口设有电磁阀1。电磁阀1的出口距离坡面有一定铅垂距离。优选地,电磁阀1的出口距离坡面的铅垂距离为10cm。电磁阀1的开启时长为固定值,确保马氏瓶6释放确定体积电解质溶液。
一个第一探针2布设在马氏瓶6的下游的侵蚀坡面。优选地,第一探针2与马氏瓶6的电解质溶液的入水口之间的距离为10cm。第一探针2测定水流电导率,确定溶质运移模型脉冲边界条件。
在第一探针2的下游的侵蚀坡面布设一个第二探针3。第二探针3测定集中水流经过该点的电导率。根据测量条件和精度要求,第二探针3能够布置多个。
第一探针2和第二探针3均测定水流电导率,其中第一探针2的数据用于确定溶质运移模型脉冲边界条件,第二探针3数据用于分析溶质运移过程。如果布置多个第二探针3,即可得到多个探针位置的流速数据。
第一探针2、第二探针3和电磁阀1均连接至数据采集器4。数据采集器4和操作控制部件5连接。
集中水流稳定后,电磁阀1开启,释放电解质溶液,同时第一探针2和第二探针3开始测定水流电导率实时变化,测定溶质运移过程。电导率实时数据存储在数据采集器4。
2、拖尾数据去除
首先根据电解质浓度和电导率的函数关系,将第一探针2和第二探针3测定的电导率数据转换为电解质浓度数据,再对电解质浓度数据归一化。
对电解质浓度数据进行归一化处理,即各数据点的电解质浓度ci均除以本次测量过程中的所有数据点的电解质浓度最大值cmax
Ci=ci/cmax
其中,i对应测量时间,ci为各数据点的电解质浓度,cmax为本次测量过程中所有数据点的电解质浓度的最大值。
拖尾是电解质浓度在达到峰值后并不以和已经测得的上升段(图2的左侧)对称的趋势降低到0值,而是长时间保持在高于背景浓度。
去除拖尾数据的方法如下:
首先,基于对称性原则,根据图2左侧归一化后的电解质浓度上升段数据预测右侧归一化后的电解质浓度下降段测量数据,并将相同时刻预测数据与测量数据比较。如果在某时刻开始,测量值均比预测值大一个预定的比例以上,则认为从该时刻以后的测量数据为拖尾数据。该预定的比例值可根据测量工况不同调整。优选地,所述预定的比例值为30%。
电解质浓度为因变量,时间为自变量作图,如图2所示。下降区,三角符号对应的电解质浓度即为由于床面吸附和阻滞作用造成的电解质浓度的拖尾数据,该部分水流电解质浓度随时间缓慢降低;圆圈符号对应的点即为反映水流中由于对流和水动力弥散作用的电解质浓度变化,为真实的溶质运移过程。
归一化的电解质浓度随时间变化如图3所示,在电解质峰值过后,有一系列数据保持较高值,形成一个比较明显的拖尾。
去除拖尾后的电解质浓度如图4所示。
3、去除拖尾数据后归一化和流速计算
去除拖尾数据点后,对剩下的数据进行归一化处理,即各数据点的电解质浓度ci均除以本次测量过程中所有数据点的电解质浓度最大值cmax
Ci=ci/cmax
其中,i对应测量时间,ci为各数据点的电解质浓度,cmax为本次测量过程中所有数据点的电解质浓度的最大值。
采用归一化后的电解质浓度Ci及方程3,采用参数拟合,即得流速、初始浓度、弥散系数:
Figure BDA0002431228910000091
其中,C是归一化后后的电解质浓度,C0是拟合得到的电解质初始浓度;x是探针布设点到马氏瓶6的电解质溶液的入水口的坡长,单位为m;t是时间,单位为s;DH是水动力弥散系数,单位为m2 s-1;u是流速,单位是m s-1
将探针测量电解质浓度分布数据带入方程3,用最小二乘法,拟合参数C0、DH和u。
u即为测量的坡面水流流速。
采用方程3及包含有拖尾数据的电解质运移过程进行参数拟合,结果如图5中虚线所示,流速值为0.12m/s,而流量法计算得流速为0.31m/s,显然拖尾数据导致模拟结果出现较大偏差。
采用方程3及去除拖尾数据的电解质运移过程进行参数拟合,结果如图6中虚线所示,流速值为0.34m/s,计算结果与流量法所得流速近似。故去除拖尾后模型所得结果为正确的水流流速。

Claims (6)

1.一种提高电解质示踪测量侵蚀坡面水流流速的方法,其特征在于:所述方法包括如下步骤:
1)电信号数据采集
在侵蚀坡面上游设置一个马氏瓶(6),马氏瓶(6)内装有高溶解度电解质的饱和溶液;
所述马氏瓶(6)出口设有电磁阀(1),电磁阀(1)的出口距离坡面有一定铅垂距离;电磁阀(1)的开启时长为固定值,确保马氏瓶(6)释放确定体积电解质溶液;
一个第一探针(2)布设在马氏瓶(6)的下游的侵蚀坡面,第一探针(2)测定水流电导率,确定溶质运移模型脉冲边界条件;
在第一探针(2)的下游的侵蚀坡面布设一个第二探针(3),第二探针(3)测定集中水流经过该点的电导率;
第一探针(2)、第二探针(3)和电磁阀(1)均连接至数据采集器(4),数据采集器(4)和操作控制部件(5)连接;
集中水流稳定后,电磁阀(1)开启,释放电解质溶液,同时第一探针(2)和第二探针(3)开始测定水流电导率实时变化,测定溶质运移过程;电导率实时数据存储在数据采集器(4);
2)拖尾数据去除
首先根据电解质浓度和电导率的函数关系,将第一探针(2)和第二探针(3)测定的电导率数据转换为电解质浓度数据,再对电解质浓度数据归一化;
对电解质浓度数据进行归一化处理,即各数据点的电解质浓度ci均除以本次测量过程中所有数据点的电解质浓度最大值cmax
Ci=ci/cmax
其中,i对应测量时间,ci为各数据点的电解质浓度,cmax为本次测量过程中所有数据点的电解质浓度的最大值;
去除拖尾数据的方法如下:
首先,基于对称性原则,根据归一化后的电解质浓度上升段数据预测归一化后的电解质浓度下降段测量数据,并将相同时刻预测的测量数据与测量数据比较;如果在某时刻开始,测量值均比预测值大一个预定的比例以上,则认为从该时刻以后的测量数据为拖尾数据;该预定的比例值可根据测量工况不同调整;
3)去除拖尾数据后归一化和流速计算
去除拖尾数据点后,对剩下的数据进行归一化处理,即各数据点的电解质浓度ci均除以本次测量过程中所有数据点的电解质浓度最大值cmax
Ci=ci/cmax
其中,i对应测量时间,ci为各数据点的电解质浓度,cmax为本次测量过程中所有数据点的电解质浓度的最大值;
采用归一化后的电解质浓度Ci及方程3,采用参数拟合,即得流速、初始浓度、弥散系数:
Figure FDA0002757721870000021
其中,C是归一化后的电解质浓度,C0是拟合得到的电解质初始浓度;x是探针布设点到马氏瓶(6)的电解质溶液的入水口的坡长,单位为m;t是时间,单位为s;DH是水动力弥散系数,单位为m2 s-1;u是流速,单位是m s-1
将探针测量电解质浓度分布数据代入方程3,用最小二乘法,拟合参数C0、DH和u;
u即为测量的坡面水流流速。
2.如权利要求1所述的提高电解质示踪测量侵蚀坡面水流流速的方法,其特征在于:步骤1中,所述高溶解度电解质的饱和溶液为氯化钾饱和溶液。
3.如权利要求1所述的提高电解质示踪测量侵蚀坡面水流流速的方法,其特征在于:步骤1中,电磁阀(1)的出口距离坡面的铅垂距离为10cm。
4.如权利要求1所述的提高电解质示踪测量侵蚀坡面水流流速的方法,其特征在于:步骤1中,第一探针(2)与马氏瓶(6)的电解质溶液的入水口之间的距离为10cm。
5.如权利要求1所述的提高电解质示踪测量侵蚀坡面水流流速的方法,其特征在于:步骤1中,根据测量条件和精度要求,第二探针(3)能够布置多个,能够得到多个探针位置的流速数据。
6.如权利要求1所述的提高电解质示踪测量侵蚀坡面水流流速的方法,其特征在于:步骤2中,所述预定的比例值为30%。
CN202010236663.9A 2020-03-30 2020-03-30 一种提高电解质示踪测量侵蚀坡面水流流速的方法 Expired - Fee Related CN111366747B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010236663.9A CN111366747B (zh) 2020-03-30 2020-03-30 一种提高电解质示踪测量侵蚀坡面水流流速的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010236663.9A CN111366747B (zh) 2020-03-30 2020-03-30 一种提高电解质示踪测量侵蚀坡面水流流速的方法

Publications (2)

Publication Number Publication Date
CN111366747A CN111366747A (zh) 2020-07-03
CN111366747B true CN111366747B (zh) 2021-02-05

Family

ID=71204926

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010236663.9A Expired - Fee Related CN111366747B (zh) 2020-03-30 2020-03-30 一种提高电解质示踪测量侵蚀坡面水流流速的方法

Country Status (1)

Country Link
CN (1) CN111366747B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112462088B (zh) * 2020-11-18 2022-02-25 中国农业大学 一种盐热耦合测量坡面水流流速的方法
CN113009177B (zh) * 2021-02-09 2022-07-05 中国农业大学 一种电解质质心法测量坡面水流流速的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1356534A (zh) * 2000-12-07 2002-07-03 中国科学院南京土壤研究所 测量地表薄层径流流速的方法
CN101825646A (zh) * 2010-05-07 2010-09-08 中国农业大学 电解质示踪测量薄层水流流速的装置及方法
CN102680732A (zh) * 2012-05-24 2012-09-19 中国农业大学 坡面薄层水流流速测量方法及系统
JP5047886B2 (ja) * 2008-06-17 2012-10-10 愛知時計電機株式会社 水流計測システム及び水流計測方法
CN108593958A (zh) * 2017-12-15 2018-09-28 北京航空航天大学 同步获取气固两相流速度场的方法和装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1219207C (zh) * 2001-12-18 2005-09-14 中国科学院水利部水土保持研究所 电解质脉冲法测量水流流速和泥沙含量的方法
CN108802417B (zh) * 2018-06-28 2024-03-29 安徽理工大学 多钻孔含水层水流流速流向测定方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1356534A (zh) * 2000-12-07 2002-07-03 中国科学院南京土壤研究所 测量地表薄层径流流速的方法
JP5047886B2 (ja) * 2008-06-17 2012-10-10 愛知時計電機株式会社 水流計測システム及び水流計測方法
CN101825646A (zh) * 2010-05-07 2010-09-08 中国农业大学 电解质示踪测量薄层水流流速的装置及方法
CN102680732A (zh) * 2012-05-24 2012-09-19 中国农业大学 坡面薄层水流流速测量方法及系统
CN108593958A (zh) * 2017-12-15 2018-09-28 北京航空航天大学 同步获取气固两相流速度场的方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Velocity of water flow along saturated loess slopes under erosion effects;Yuhan Huang;《Journal of Hydrology》;20180403;第304-311页 *
降雨条件下电解质示踪正态模型的应用;史晓楠 等;《中国农业大学学报》;20101231;第15卷(第4期);第124-129页 *

Also Published As

Publication number Publication date
CN111366747A (zh) 2020-07-03

Similar Documents

Publication Publication Date Title
CN111366747B (zh) 一种提高电解质示踪测量侵蚀坡面水流流速的方法
Liu et al. Structure of air-water bubbly flow in a vertical pipe—II. Void fraction, bubble velocity and bubble size distribution
Thulasidas et al. Dispersion during bubble-train flow in capillaries
CN108254032A (zh) 河流超声波时差法流量计算方法
Zhang et al. Measurements of dust deposition velocity in a wind-tunnel experiment
Li et al. Effect of saltating sediment load on the determination of the mean velocity of overland flow
Cheng et al. Turbulent open-channel flow with upward seepage
CN109060638B (zh) 岩心渗透率确定方法及装置
CN1793881A (zh) 土壤水分、溶质运移参数自动测定装置
CN109583035B (zh) 基于元胞自动机的城市面源污染物累积冲刷过程计算方法
CN106971085A (zh) 一种泥石流沟道侵蚀量的计算方法
CN109736788A (zh) 一种判别化学驱前缘波及状态的实验方法
Shao et al. A method for directly measuring the hydraulic conductivity of unsaturated soil
Baas et al. Duration of deposition from decelerating high-density turbidity currents
CN113737710B (zh) 天然河道推移质输沙率的估算方法
CN109063370B (zh) 漫滩弯曲河道水深平均水流流向预测方法
CN107090798B (zh) 一种水库超饱和输沙状态下恢复饱和系数计算方法
Guymer et al. Longitudinal dispersion due to surcharged manhole
Piedra-Cueva et al. Erosion of a deposited layer of cohesive sediment
CN108827871B (zh) 一种管式泥沙侵蚀试验装置中泥沙表面切应力确定方法
CN116593359A (zh) 一种确定输水渠道藻沙混合物动水沉速和有效沉降距离的方法
CN113009177B (zh) 一种电解质质心法测量坡面水流流速的方法
CN105823846B (zh) 一种基于分配示踪法的过硫酸盐氧化修复分析方法
Wei et al. Experimental study of air-water interface properties in self-aerated flows
CN112462088B (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210205

CF01 Termination of patent right due to non-payment of annual fee