CN102830246B - 一种测量泥沙启动流速的方法 - Google Patents

一种测量泥沙启动流速的方法 Download PDF

Info

Publication number
CN102830246B
CN102830246B CN 201210314718 CN201210314718A CN102830246B CN 102830246 B CN102830246 B CN 102830246B CN 201210314718 CN201210314718 CN 201210314718 CN 201210314718 A CN201210314718 A CN 201210314718A CN 102830246 B CN102830246 B CN 102830246B
Authority
CN
China
Prior art keywords
concentration
velocity
silt
suspended particles
water
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
CN 201210314718
Other languages
English (en)
Other versions
CN102830246A (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN 201210314718 priority Critical patent/CN102830246B/zh
Publication of CN102830246A publication Critical patent/CN102830246A/zh
Application granted granted Critical
Publication of CN102830246B publication Critical patent/CN102830246B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Abstract

本发明涉及一种测量泥沙启动流速的方法,包括如下步骤:、从流速为零开始,利用B超仪测量并记录河床边界上方悬移粒子浓度,所述粒子浓度为每帧图像所含悬移粒子的数目,用N表示;、逐步增加水流速度,测量并记录悬浮粒子浓度及浓度变化率,所述水流速度用V表示,单位cm/s,所述浓度变化率为ΔN/ΔV;、当悬浮粒子浓度超过静水中悬浮粒子浓度的3~6倍且N>30,同时悬浮粒子的浓度变化率达到100~200颗/(cm/s)时,判定此时的水流速度为泥沙启动流速,测量此时的水流速度为泥沙的启动流速。本发明方法能够对水流中悬移沙粒进行定量分析,从而建立定量、客观、一致的判定标准。

Description

一种测量泥沙启动流速的方法
技术领域
本发明属于河工测量技术领域,涉及一种泥沙启动流速的测量方法。
背景技术
在河工模型试验中,泥沙启动流速的测量对于泥沙运移规律的研究具有十分重要的意义。泥沙启动是指在一定的水流条件下,泥沙由静止转为运动的现象。然而,现有泥沙启动流速测量的方法还比较落后,主要依靠人工观察的方法,即由人工调整试验水槽中的流速,并通过人工目视观察试验水槽床面泥沙的启动情况,进而判定是否达到泥沙的启动流速。这种方法存在以下问题:1、受人为因素影响较大,不同人员、不同的习惯、不同光照条件可能得出不同临界启动流速;2、缺乏合理的、定量的、客观一致的判别依据,人为随意性很大。3、受水流透明性影响较大。在含沙水流中往往会有很多微小的尘粒,长时间悬浮在水中,尤其是流动的水中,使得水体混浊,实验人员根本无法看清床面泥沙的运动情况,因而无法判定泥沙的启动状态。
发明内容
本发明要解决的技术问题是提供一种基于仪器的测量泥沙启动流速的方法。本发明采用B超成像,自动判别泥沙启动状态和测量启动流速。
当流速低于泥沙的启动流速时,除少数微小粒子流动外,床面的沙粒是不动的,在床底附近悬浮运动的沙粒浓度很低,只有少数杂散悬浮粒子。显然,在这种情况下,河床床面也保持不变;当水流速度达到并超过泥沙临界启动流速时,床面沙粒的力学平衡条件被打破,就会发生沙粒沿床面的滚动和部分沙粒上扬运动,伴随而至的是床面附近悬浮运动粒子的浓度的显著增加。经历一定的时间后,床面的地形会发生一定的变化;进一步增大水流速度,将伴随大量沙粒的扬起和流动,运动沙粒的浓度将急剧增加,同时,河床地貌也将发生快速的变化。
显然,测量出河床边界上方悬移沙粒的数量(或浓度)及其变化情况,就可以判定泥沙的启动与否,对应流速由流速仪测得,从而测量出泥沙的临界启动流速。
由于超声波在水下有很好的传导性,而且高频超声波对水中的粒子较为敏感,借助于B型超声扫描方式,只要选择合适的超声频率,进行水下二维成像,可以获得河床界面及其上方的悬浮沙粒的的实时图像信号。依据这些粒子的光斑图像信号可以分析判定泥沙的启动状况。虽然图像中的粒子的大小可能会有一定失真,但足以反映悬移粒子的存在和粒子的浓度变化情况,为分析判定泥沙的启动状态提供足够的信息。
实验证明,在泥沙临界启动流速附近,河床界面悬移沙粒的浓度有一个显著的增加的过程,这为通过B超成像准确判定泥沙的启动状况提供了有效特征。
定义术语:
粒子浓度:每帧图像所含悬移粒子的数目,用N表示。
基础浓度:静水中每帧图像所含悬移粒子浓度,用N0表示。
水流速度:用V表示,单位cm/s。
浓度变化率:ΔN/ΔV,  用U表示。
相对浓度: N/(N0), 用S表示。
根据不同沙质多次实验结果和数据的统计分析,本发明建立了以下基于B超成像的泥沙启动流速的判定法则:
1、粒子浓度N显著增加到基础浓度N0的3~6倍以上且N>30;
2、粒子浓度的相对变化率达到100~200颗/cm/s。
   以上两个条件同时满足对应的最小流速即为泥沙的启动流速。
基于上述判定准则,本发明提供的测量泥沙启动流速的方法为:
                                                 
Figure 2012103147189100002DEST_PATH_IMAGE001
、从流速为零开始,利用B超仪测量并记录河床边界上方悬移粒子浓度,所述粒子浓度为每帧图像所含悬移粒子的数目,用N表示;
Figure 261892DEST_PATH_IMAGE002
、逐步增加水流速度,测量并记录悬浮粒子浓度及浓度变化率,所述水流速度用V表示,单位cm/s,所述浓度变化率为ΔN/ΔV;
Figure 2012103147189100002DEST_PATH_IMAGE003
、当悬浮粒子浓度超过静水中悬浮粒子浓度的3~6倍且N>30,悬浮粒子的浓度变化率达到100~200时判定此时的水流速度为泥沙启动流速,测量此时的水流速度为泥沙的启动流速。
    测量水流速度可用常用的流速仪如旋桨流速仪、ADV流速仪或电磁流速仪等。
本发明方法采用B超成像仪,对启动流速试验设施的水下河床边界和悬移沙粒直接进行成像,获取泥沙粒子运动的直观视频图像信息,通过分析床面附近悬移粒子浓度的变化特征,结合流速仪测得的流速,最终给出泥沙的启动流速。这种方法有以下优点:1、适应性较强。由于超声波在水下有良好的穿透性,混水流和清水流都能获得比较清晰的图像。2、能够对水流中悬移沙粒进行定量分析,从而建立定量、客观、一致的判定标准。3、可利用计算机自动分析处理,实现启动流速的自动测量。
附图说明
图1平均粒径为0.1mm沙粒图像。
图2平均粒径为0.25mm沙粒图像。
图3平均粒径为0.5mm沙粒成像。
图4泥沙启动前的水下图像 。
图5泥沙启动后的水下图像 。
图6粒子浓度与流速的关系曲线。
具体实施方式
实施例1,通过B超对静止水中塑料沙进行成像
 试验证实,根据所测粒子的粒度大小,采用医用B超仪,选择合适的超声波频率(2.5~7.5MHZ)和适当的信号增益,完全可以获得水中沙粒的反射图像。
图1,2,3是实际采集到的不同粒径时悬移粒子图像。可以看出,B超仪可以获得粒度小至0.1mm的沙粒的清晰图像。
实施例2通过B超对水中泥沙启动前后进行成像
根据所测粒子的粒度大小,采用医用B超仪,选择合适的超声波频率(2.5~7.5MHZ)和适当的信号增益,获得水中沙粒启动流速前后河床界面附近水下悬移粒子图像信号。如图4,5所示。图中亮线为河床边界,河床上方的亮斑为悬移沙粒的图像。由图可以看出,启动前后河床边界上方悬移粒子数量发生了明显变化。
实施例3以模型试验塑料沙为例,图6是该方法实测得到的悬移粒子浓度与流速的关系曲线。
从中可以看出:当流速小于3.11cm/s时,床面附近悬移粒子浓度很低,曲线平坦,颗粒浓度变化不大,约为每帧8~11颗左右,平均浓度为10颗/每帧(以下称为基础浓度,这是由其它一些轻质的悬浮物造成的,并非真正的悬浮沙粒)。此时,人工目视观察河床界面处没有明显沙粒的启动运动现象;当流速增加到3.21cm/s时,悬移粒子浓度由每帧10颗增加到每帧27颗,粒子浓度相对于流速的变化率(即曲线的斜率)为:170颗/cm/s,相对变化倍率(与基础浓度的比值)为2.7,目视观察河床界面处有沙粒的轻微移动;当流速增加到3.33 cm/s时,悬移粒子浓度猛增到N=47颗/每帧,粒子浓度的变化率为:168颗/cm/s,相对变化倍率为4.7,人工目视观察河床界面处有沙粒的明显移动;当流速增加到3.45 cm/s时,悬移粒子浓度猛增到76颗/每帧,粒子浓度变化率为:194颗/cm/s,相对变化倍率为7.6。此时,目视观察河床界面处有沙粒的显著移动;很显然,从上面数据和目视观察的结果来看,在流速为3.21cm/s时,发生了泥沙的轻微启动,在流速为3.45cm/s时泥沙发生了显著运动,所以泥沙临界启动流速应在3.21~3.45cm/s之间,可取临界启动流速为3.33cm/s。不确定度为±0.1cm/s,完全满足启动流速测量精度的要求。

Claims (2)

1.一种测量泥沙启动流速的方法,其特征在于包括如下步骤:
                                                
Figure 2012103147189100001DEST_PATH_IMAGE002
、从流速为零开始,利用B超仪测量并记录河床边界上方悬浮粒子浓度,所述粒子浓度为每帧图像所含悬浮粒子的数目,用N表示;
、逐步增加水流速度,测量并记录悬浮粒子浓度及浓度变化率,所述水流速度用V表示,单位cm/s,所述浓度变化率为ΔN/ΔV;
Figure 2012103147189100001DEST_PATH_IMAGE006
、当悬浮粒子浓度超过静水中悬浮粒子浓度的3~6倍且N>30,同时悬浮粒子的浓度变化率达到100~200颗/(cm/s)时,判定此时的水流速度为泥沙启动流速,测量此时的水流速度为泥沙的启动流速。
2.根据权利要求1所述的测量泥沙启动流速的方法,其特征在于,测量水流速度用旋桨流速仪、ADV流速仪或电磁流速仪。
CN 201210314718 2012-08-30 2012-08-30 一种测量泥沙启动流速的方法 Expired - Fee Related CN102830246B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210314718 CN102830246B (zh) 2012-08-30 2012-08-30 一种测量泥沙启动流速的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210314718 CN102830246B (zh) 2012-08-30 2012-08-30 一种测量泥沙启动流速的方法

Publications (2)

Publication Number Publication Date
CN102830246A CN102830246A (zh) 2012-12-19
CN102830246B true CN102830246B (zh) 2013-12-18

Family

ID=47333451

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210314718 Expired - Fee Related CN102830246B (zh) 2012-08-30 2012-08-30 一种测量泥沙启动流速的方法

Country Status (1)

Country Link
CN (1) CN102830246B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114739629B (zh) * 2022-04-13 2023-01-17 武汉大学 一种泥沙起动流速测量装置
CN115655994B (zh) * 2022-09-15 2023-09-19 浙江天禹信息科技有限公司 一种水域泥沙超声波测量方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4842774A (en) * 1987-08-07 1989-06-27 The United States Of America As Represented By The United States Department Of Energy Pyramiding tumuli waste disposal site and method of construction thereof
WO2009018597A1 (de) * 2007-08-09 2009-02-12 Technische Universität Graz Vorrichtung zum bestimmen von strömungsparametern einer partikel - fluidum - strömung
CN101533035A (zh) * 2009-04-02 2009-09-16 华东师范大学 河口海岸高浊度环境近底水沙观测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4807626B2 (ja) * 2006-12-01 2011-11-02 東京電力株式会社 画像解析システム及びコンピュータプログラム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4842774A (en) * 1987-08-07 1989-06-27 The United States Of America As Represented By The United States Department Of Energy Pyramiding tumuli waste disposal site and method of construction thereof
WO2009018597A1 (de) * 2007-08-09 2009-02-12 Technische Universität Graz Vorrichtung zum bestimmen von strömungsparametern einer partikel - fluidum - strömung
CN101533035A (zh) * 2009-04-02 2009-09-16 华东师范大学 河口海岸高浊度环境近底水沙观测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JP特开2008-140103A 2008.06.19
基于粒子运动分析的慢流速超声测量技术研究;赵小红等;《仪器仪表学报》;20080430;第29卷(第4期);845-849 *
赵小红等.基于粒子运动分析的慢流速超声测量技术研究.《仪器仪表学报》.2008,第29卷(第4期),845-849.

Also Published As

Publication number Publication date
CN102830246A (zh) 2012-12-19

Similar Documents

Publication Publication Date Title
CN104535475B (zh) 碳酸盐岩微观结构的确定方法及装置
Honkanen et al. Digital imaging measurement of dense multiphase flows in industrial processes
Florens et al. Defining the roughness sublayer and its turbulence statistics
CN106228532B (zh) 一种lcd屏幕亚像素级缺陷检测方法
CN110473260A (zh) 一种波浪视频测量装置及方法
CN103033560B (zh) 一种基于b超成像技术的低含沙量测量方法
CN102830246B (zh) 一种测量泥沙启动流速的方法
Larcher et al. Set of measurement data from flume experiments on steady uniform debris flows
CN103839281A (zh) 一种利用造影微泡测量血流速度的方法
Ziegenhein et al. Particle tracking using micro bubbles in bubbly flows
CN104458542A (zh) 非接触式砂石骨料粒度粒形检测仪及检测方法
CN102889907B (zh) 一种基于视觉检测的流量监控装置及流量监控方法
CN109272548A (zh) 一种浮选过程气泡直径的测量方法
CN107631960A (zh) 海洋物质输移扩散系数ccd自动监测实验装置及其监测方法
CN110501268A (zh) 一种基于显微图像处理的微量粉尘检测方法
CN202770854U (zh) 一种测量泥沙启动流速的装置
CN116026855A (zh) 基于视觉分析的固体颗粒检验用样品分析系统
Zou et al. B-mode ultrasound imaging measurement and 3D reconstruction of submerged topography in sediment-laden flow
CN204594856U (zh) 非接触式砂石骨料粒度粒形检测仪
Osorio-Cano et al. A method for extracting surface flow velocities and discharge volumes from video images in laboratory
Terwisscha van Scheltinga et al. Sand particle velocities over a subaqueous dune slope using high‐frequency image capturing
CN204065126U (zh) 一种岸滩表面流速粒子追踪测量系统
Abid et al. Volume and rate measurement of slowly generated gas bubbles
Gaeta et al. Non-intrusive measurements of wave-induced flow over dikes by means of a combined ultrasound doppler velocimetry and videography
CN115908271A (zh) 一种泵送混凝土坍落度在线监测方法、装置及设备

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

Granted publication date: 20131218

Termination date: 20150830

EXPY Termination of patent right or utility model