CN107607482A - 河流动态泥沙含量测量方法 - Google Patents
河流动态泥沙含量测量方法 Download PDFInfo
- Publication number
- CN107607482A CN107607482A CN201710772246.4A CN201710772246A CN107607482A CN 107607482 A CN107607482 A CN 107607482A CN 201710772246 A CN201710772246 A CN 201710772246A CN 107607482 A CN107607482 A CN 107607482A
- Authority
- CN
- China
- Prior art keywords
- sediment
- river
- sample
- silt
- mixed solution
- 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
Links
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了河流动态泥沙含量测量方法,通过在河流不同高度取样烘干,加水配置混合溶液,通过单色光测量标准混合溶液的吸光度,求出最优单色光波长范围,构建稳态模型,通过模型计算出河流的泥沙含量,方便快捷。
Description
技术领域
本发明属于环境监测技术领域,尤其涉及一种河流动态泥沙含量测量方法。
背景技术
随着经济社会的发展以及光谱技术的进步,光谱技术在水环境监测与水 质参数定量反演的应用日趋广泛,方法也不断发展与成熟。根据实测透射光 谱数据,可以分析研究区域的悬浮泥沙浓度与粒度分布规律,初步探究两者 之间的关系,以及光谱与悬浮泥沙浓度、粒度变化的响应关系,但是,光谱 测量泥沙仍然存在偶然误差大,测量精度低等问题。实际河流是几乎不存在 静止的泥沙沉降,河流中时刻存在向下的水流,因此,水流流速也将对泥沙 产生影响,另一方面,泥沙将可能产生起动,进一步影响泥沙沉降。泥沙的 起动指的是在一定的床面上,逐渐增大水流强度,直到使床沙从静止转入运 动,即使对单颗泥沙颗粒来说,由于所受的瞬时作用底速或拖曳力有脉动特 性,由于沙粒大小、形状、位置具有变异性,所以,额泥沙的起动是具有一 定随机性的,但在一定的特定条件下,从宏观上进行观察的话,可以从实验 学的现象观察中得出泥沙具有弱动、中动、普动等三种不同判别标准,存在 一定的宏观必然性。当水流中含有泥沙时,水流的紊动特性会发生明显的改变,从而引起泥沙的一些特性如沉速发生改变,给准确测量泥沙含量带来困 难。
发明内容
针对现有技术存在的问题,本发明提供河流动态泥沙含量测量方法,
本发明的上述技术问题主要是通过下述技术方案得以解决的:
本发明的河流动态泥沙含量测量方法,包括如下步骤:
1)将6个取样桶放入河流中取样得到泥沙标本,具体为:
将6个同规格圆筒型取样桶(圆筒直径DD、高HH)等间距叠放在一个位 于其轴心的固定杆上,每个取样桶可以沿固定杆上下方向调节其与相邻取 样桶间的间距HHH,将固定杆竖直插入河流待测点,调节取样桶间距HHH, 使得最上方第1个取样桶上边沿与第6个取样桶上边沿间距等于待测点河 流深度H的5/6,每个取样桶间距相等,最上方第1个取样桶上边沿与水面 间距为河流深度H的1/12,经过时间T小时后取出取样桶,沥干其内部的 水分,然后再用烘干机去除水分,得到烘干的泥沙标本,称重,得到质量 m1、m2、m3、m4、m5、m6(单位:g)的泥沙标本,分别贴上标签1、2、3、 4、5、6备用;
2)将经过烘干的泥沙标本,混合均匀并加水配成标准混合溶液,具体为: 将经过烘干的泥沙标本,按比例(m1/m6):(m2/m6):(m2/m6):(m4/m6): (m5/m6):(m6/m6)(单位:g)均匀混合并加水配成标准混合溶液, 以备实验测量。
3)通过单色光测量标准混合溶液的吸光度,求出最优单色光波长范围, 具体为:先摇晃均匀标准混合溶液,混合溶液质量为1000g,通过单 色光测上述混合溶液的吸光度A,单色光的波长范围在400nm到800nm 之间取值,由朗伯-比尔定律A=εCL可以得到泥沙浓度其中, L为光程长度,ε为吸光度系数,针对不同河流还可以通过实验室光 谱测量来优选其最优单色光波长,例如采用紫外可见分光光度计测量 上述标准混合溶液在400nm—800nm之间的吸光度,得到泥沙溶液的 光谱图,得到吸光度没有突变扰动的数个优化测量区间,分别计算每 个优化测量区间与标准混合溶液浓度的线性相关系数R,得到一个相 关系数最大的最优测量区间,即得到最优单色光的波长范围。
4)继续测量标准混合溶液吸光度与时间的关系A=αe-βt,t是溶液水平静 置时间,求出系数α,β。
5)由A=εCL与A=αe-βt,即可以得到静水条件下泥沙浓度C与所处深度h 关系εCL=αe-β(h/ω),其中,h为测量点到水面的垂直距离(即水的深 度),ω为泥沙沉降速度,其大小取决于泥沙的粒径,可以通过实验求 出。
6)根根据泥沙与水流的情况求出起动系数、水流对泥沙沉降系数以及泥 沙沉降系数,再通过计算得到泥沙含量,具体为:考虑河流中时刻存 在向下的水流与泥沙的起动,设泥沙起动数为D,沉降数为J,水流 流速为V,水流流速与水平面的夹角为θ,则D=αdV sinθ, 其中,αd为起动系数,βj为水流对泥沙沉降系数, γ为泥沙沉降系数,设编号6泥沙含量为C6,可以得到编号1、2、3、 4、5的泥沙含量分别为m1*C6/m6、m2*C6/m6、m3*C6/m6、m4*C6/m6、 m5*C6/m6,分别将编号1、2、3、4、5、6的泥沙含量m1*C6/m6、m2*C6/m6、 m3*C6/m6、m4*C6/m6、m5*C6/m6、C6代入上面两个式子,可以求出αd、 βj、γ;
当泥沙起动数等于沉降数,即D=J时,河流泥沙分布达到稳态分布, 由公式即可求出C(h);
令h0点的泥沙浓度等于整个河流整体的泥沙浓度 的平均值,因此我们要准确测量某点河流动态泥沙含量,只要将测量点或 取水点设置在某点深度h0处即可。
本发明的有益效果是:为了准确测量某点河流动态泥沙含量,基于泥沙含 量吸光度测量实验和泥沙沉降规律,以及实际河流中存在的水流及泥沙启动现 象,建立了一套河流动态泥沙含量准确测定方法。该方法不仅可以准确测量河 流的泥沙含量,还具有测量方法简单、实用的特点,对特定的监测点而言,通 过区间优化以及监测点位置选取,可以实现通过单波长吸光度的测量来准确测 量该点以及整个河流动态泥沙含量,只要经过一个取样点的单波长测量即可得 到河流泥沙准确含量,简化了繁杂的测量操作,降低了数据采集数量以及计算 量,尤其适合水质实时在线的动态监控与长周期、大数据的实验测量研究
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施 例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述 中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付 出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为河流动态泥沙含量的测量方法流程图。
具体实施方式
下面结合附图对本发明的优选实施例进行详细阐述,以使本发明的优点和 特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚 明确的界定。
如图1所示,河流动态泥沙含量测量方法,包括如下步骤:
(1)将6个同规格圆筒型取样桶(圆筒直径DD 5cm、高HH 20cm)等间距 叠放在一个位于其轴心的固定杆上,每个取样桶可以沿固定杆上下方向调节其 与相邻取样桶间的间距HHH,将固定杆竖直插入南京某河流待测点,待测点河流 深度H=3m,调节取样桶间距HHH为50cm,使得最上方第1个取样桶上边沿与第 6个取样桶上边沿间距等于250cm,每个取样桶间距50cm,最上方第1个取样桶 上边沿与水面间距为25cm,经过时间T=2小时后取出取样桶,沥干其内部的水 分,然后再用烘干机去除水分,得到烘干的泥沙标本,称重,得到质量m1=0.15g、 m2=0.22g、m3=0.40g、m4=0.86g、m5=1.93g、m6=4.35g的泥沙标本6组,分别 贴上标签1、2、3、4、5、6备用;
(2)分别取上述6组标本0.03g、0.05g、0.09g、0.20g、0.44g、1g均匀 混合并加水配成1000g的混合溶液,得到泥沙浓度C=0.001* (m1+m2+m3+m4+m5+m6)/m6=1.82g/L的标准混合溶液;
(3)先摇晃均匀标准混合溶液,通过单色光测上述混合溶液的吸光度 A=3.83656,单色光的波长范围在400nm到800nm之间取值,由朗伯-比尔定律 A=εCL可以得到泥沙浓度其中,L=5cm为光程长度,求出吸光度系数 ε=0.4216L/(g*cm),采用紫外可见分光光度计测量上述6种不同标签配置的泥 沙溶液在400nm—800nm之间的吸光度,得到泥沙溶液的光谱图,分析上述6组 光谱图,得到吸光度没有突变扰动的数个优化测量区间,分别计算每个优化测 量区间与6种不同标签配置的泥沙溶液浓度的线性相关系数R,得到一个相关系 数最大的最优测量区间,最终得到最优单色光的波长范围为480--520nm;
(4)继续测量标准混合溶液吸光度与时间的关系,由A=αe-βt求出系数 α=1.2621,β=0.025157。
(5)由A=εCL与A=αe-βt,即可以得到静水条件下泥沙浓度C与所处深度h 关系εCL=αe-β(h/ω),其中,h为测量点到水面的垂直距离(即水的深度),本待 测点泥沙沉降速度ω为0.262cm/s。
(6)本待测点水流流速为V=0.96m/s,水流流速与水平面的夹角θ为9.8度。 由αd=0.4635、βj=0.4632、γ=0.0005,可得到D=αdVsinθ0.07574,
(7)当泥沙起动数等于沉降数,即D=J时,河流泥沙分布达到稳态分布, 由公式即可求出C(h)=-0.00688*e^(h/0.262)。
(8)令h0点的泥沙浓度等于整个河流整体的泥沙浓度 的平均值,因此我们要准确测量本待测点河流动态泥沙含量,只要将测量取水 点设置在本待测点深度h0=2.36m处即可。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于 此,任何不经过创造性劳动想到的变化或替换,都应涵盖在本发明的保护范围 之内。因此,本发明的保护范围应该以权利要求书所限定的保护范围为准。
Claims (6)
1.一种河流动态泥沙含量测量方法,其特征在于,包括如下步骤:
1)将6个取样桶放入河流中取样得到泥沙标本;
2)将经过烘干的泥沙标本,混合均匀并加水配成标准混合溶液;
3)通过单色光测量标准混合溶液的吸光度,求出最优单色光波长范围;
4)继续测量标准混合溶液吸光度与时间的关系A=αe-βt,t是溶液水平静置时间,求出系数α,β;
5)由A=εCL与A=αe-βt,即可以得到静水条件下泥沙浓度C与所处深度h关系εCL=αe-β(h/ω),其中,h为测量点到水面的垂直距离,即水的深度,ω为泥沙沉降速度;
6)根根据泥沙与水流的情况求出起动系数、水流对泥沙沉降系数以及泥沙沉降系数,再通过计算得到泥沙含量。
2.根据权利要求1所述的一种河流动态泥沙含量测量方法,其特征在于:所述步骤1)将6个取样桶放入河流中取样得到泥沙标本具体为:将6个同规格圆筒型取样桶等间距叠放在一个位于其轴心的固定杆上,每个取样桶可以沿固定杆上下方向调节其与相邻取样桶间的间距HHH,将固定杆竖直插入河流待测点,调节取样桶间距HHH,使得最上方第1个取样桶上边沿与第6个取样桶上边沿间距等于待测点河流深度H的5/6,每个取样桶间距相等,最上方第1个取样桶上边沿与水面间距为河流深度H的1/12,经过时间T小时后取出取样桶,沥干其内部的水分,然后再用烘干机去除水分,得到烘干的泥沙标本,称重,得到质量m1、m2、m3、m4、m5、m6的泥沙标本,分别贴上标签1、2、3、4、5、6备用。
3.根据权利要求1所述的一种河流动态泥沙含量测量方法,其特征在于:所述步骤2)将经过烘干的泥沙标本,混合均匀并加水配成标准混合溶液具体为:将经过烘干的泥沙标本,按比例(m1/m6):(m2/m6):(m2/m6):(m4/m6):(m5/m6):(m6/m6)均匀混合并加水配成标准混合溶液,以备实验测量。
4.根据权利要求1所述的一种河流动态泥沙含量测量方法,其特征在于:所述步骤3)通过单色光测量标准混合溶液的吸光度,求出最优单色光波长范围具体为:先摇晃均匀标准混合溶液,通过单色光测上述混合溶液的吸光度A,单色光的波长范围在400nm到800nm之间取值,由朗伯-比尔定律A=εCL可以得到泥沙浓度其中,L为光程长度,ε为吸光度系数,针对不同河流还可以通过实验室光谱测量来优选其最优单色光波长,例如采用紫外可见分光光度计测量上述标准混合溶液在400nm—800nm之间的吸光度,得到泥沙溶液的光谱图,得到吸光度没有突变扰动的数个优化测量区间,分别计算每个优化测量区间与标准混合溶液浓度的线性相关系数R,得到一个相关系数最大的最优测量区间,即得到最优单色光的波长范围。
5.根据权利要求1所述的一种河流动态泥沙含量测量方法,其特征在于:所述步骤6)根根据泥沙与水流的情况求出起动系数、水流对泥沙沉降系数以及泥沙沉降系数,再通过计算得到泥沙含量具体为:考虑河流中时刻存在向下的水流与泥沙的起动,设泥沙起动数为D,沉降数为J,水流流速为V,水流流速与水平面的夹角为θ,则D=αdVsinθ,其中,αd为起动系数,βj为水流对泥沙沉降系数,γ为泥沙沉降系数,设编号6泥沙含量为C6,可以得到编号1、2、3、4、5的泥沙含量分别为m1*C6/m6、m2*C6/m6、m3*C6/m6、m4*C6/m6、m5*C6/m6,分别将编号1、2、3、4、5、6的泥沙含量m1*C6/m6、m2*C6/m6、m3*C6/m6、m4*C6/m6、m5*C6/m6、C6代入上面两个式子,可以求出αd、βj、γ;
当泥沙起动数等于沉降数,即D=J时,河流泥沙分布达到稳态分布,由公式即可求出C(h);
令h0点的泥沙浓度等于整个河流整体的泥沙浓度的平均值,因此我们要准确测量某点河流动态泥沙含量,只要将测量点或取水点设置在某点深度h0处即可。
6.根据权利要求1所述的一种河流动态泥沙含量测量方法,其特征在于:所述混合溶液的质量为1000g。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710772246.4A CN107607482B (zh) | 2017-08-31 | 2017-08-31 | 河流动态泥沙含量测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710772246.4A CN107607482B (zh) | 2017-08-31 | 2017-08-31 | 河流动态泥沙含量测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107607482A true CN107607482A (zh) | 2018-01-19 |
CN107607482B CN107607482B (zh) | 2020-06-05 |
Family
ID=61056755
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710772246.4A Active CN107607482B (zh) | 2017-08-31 | 2017-08-31 | 河流动态泥沙含量测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107607482B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112748079A (zh) * | 2020-12-23 | 2021-05-04 | 中山艾尚智同信息科技有限公司 | 一种测定建筑用天然砂含泥量的方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1526024A (zh) * | 2001-02-07 | 2004-09-01 | ��ʡ��ѧԺ | 光电子探测系统 |
CN101246112A (zh) * | 2008-03-20 | 2008-08-20 | 福州大学 | 水体中悬浮泥沙含量的测量方法及测量装置 |
CN102288646A (zh) * | 2011-05-16 | 2011-12-21 | 中国海洋大学 | 一种海水中泥沙浓度垂向分布的原位自动监测方法及系统 |
CN102359862A (zh) * | 2011-08-12 | 2012-02-22 | 河海大学 | 粉沙质和淤泥质海岸泥沙运动数值模拟方法 |
CN103512829A (zh) * | 2013-10-24 | 2014-01-15 | 中国水利水电科学研究院 | 一种通过浑水密度测量泥沙浓度的方法 |
CN103940769A (zh) * | 2014-04-16 | 2014-07-23 | 河海大学 | 胶体在潜流带中输运的水槽实验方法 |
WO2015124718A1 (en) * | 2014-02-20 | 2015-08-27 | Dsm Ip Assets B.V. | Phage insensitive streptococcus thermophilus |
CN105913123A (zh) * | 2016-04-12 | 2016-08-31 | 西北工业大学 | 基于自动编码器和多层感知器网络的黄河主溜光谱建模方法 |
CN106370563A (zh) * | 2016-08-26 | 2017-02-01 | 南京昊控软件技术有限公司 | 利用光学法测量模型垂向断面泥沙浓度的方法和装置 |
-
2017
- 2017-08-31 CN CN201710772246.4A patent/CN107607482B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1526024A (zh) * | 2001-02-07 | 2004-09-01 | ��ʡ��ѧԺ | 光电子探测系统 |
CN101246112A (zh) * | 2008-03-20 | 2008-08-20 | 福州大学 | 水体中悬浮泥沙含量的测量方法及测量装置 |
CN102288646A (zh) * | 2011-05-16 | 2011-12-21 | 中国海洋大学 | 一种海水中泥沙浓度垂向分布的原位自动监测方法及系统 |
CN102359862A (zh) * | 2011-08-12 | 2012-02-22 | 河海大学 | 粉沙质和淤泥质海岸泥沙运动数值模拟方法 |
CN103512829A (zh) * | 2013-10-24 | 2014-01-15 | 中国水利水电科学研究院 | 一种通过浑水密度测量泥沙浓度的方法 |
WO2015124718A1 (en) * | 2014-02-20 | 2015-08-27 | Dsm Ip Assets B.V. | Phage insensitive streptococcus thermophilus |
CN103940769A (zh) * | 2014-04-16 | 2014-07-23 | 河海大学 | 胶体在潜流带中输运的水槽实验方法 |
CN105913123A (zh) * | 2016-04-12 | 2016-08-31 | 西北工业大学 | 基于自动编码器和多层感知器网络的黄河主溜光谱建模方法 |
CN106370563A (zh) * | 2016-08-26 | 2017-02-01 | 南京昊控软件技术有限公司 | 利用光学法测量模型垂向断面泥沙浓度的方法和装置 |
Non-Patent Citations (3)
Title |
---|
方彦军 等: "含沙量测量的新进展", 《武汉水利电力大学学报》 * |
王丽萍 等: "河水泥沙含量测量方法研究", 《辽宁师范大学学报(自然科学版)》 * |
路炳军 等: "径流小区集流桶( 池) 泥沙含量快速测定方法探讨", 《水土保持通报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112748079A (zh) * | 2020-12-23 | 2021-05-04 | 中山艾尚智同信息科技有限公司 | 一种测定建筑用天然砂含泥量的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107607482B (zh) | 2020-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dobbins | Effect of turbulence on sedimentation | |
Caldwell et al. | A laboratory study of the turbulent Ekman layer | |
Li et al. | Boundary layer dynamics and drag reduction in flows of high cohesive sediment suspensions | |
Rahma et al. | Measuring flow velocity under straw mulch using the improved electrolyte tracer method | |
CN103512829A (zh) | 一种通过浑水密度测量泥沙浓度的方法 | |
Goto et al. | Turbulence estimation using fast-response thermistors attached to a free-fall vertical microstructure profiler | |
CN110659783B (zh) | 有植被群落河道横断面平均流速纵向分布预测方法 | |
Whalley et al. | An experimental investigation into spatiotemporal intermittencies in turbulent channel flow close to transition | |
Klingbeil et al. | Thickness-weighted averaging in tidal estuaries and the vertical distribution of the Eulerian residual transport | |
Termini | Turbulent mixing and dispersion mechanisms over flexible and dense vegetation | |
Salazar et al. | Time-average local thickness measurement in falling liquid film flow | |
CN107607482A (zh) | 河流动态泥沙含量测量方法 | |
Nicosia et al. | Roughness effect on the correction factor of surface velocity for rill flows | |
Coundoul et al. | Role of local flow conditions in river biofilm colonization and early growth | |
Boyer | Rotating flow over long shallow ridges | |
CN203519460U (zh) | 一种涂料粘度测定仪的智能清洗检控系统 | |
Boyer et al. | Rotating open channel flow past right circular cylinders | |
Zagumennyi et al. | Diffusion induced flows on a strip: Theoretical, numerical and laboratory modeling | |
CN103487352B (zh) | 一种涂料粘度测定仪的智能清洗检控系统及其检控方法 | |
Xue et al. | Measurement of seawater surface tension coefficient based on bubble rising behavior | |
Moussouni et al. | The effect of raindrop Kinetic energy on soil erodibility | |
Li | Sheet flow under simulated rainfall | |
Obukhov et al. | On the microstructure of atmospheric turbulence–A review of recent work in the USSR | |
Islam et al. | Establishing morpho-dynamic baseline for flow-sediment management of a tidal river in the Ganges–Brahmaputra–Meghna Delta system through field measurement | |
da Costa Araújo | Erosion mechanisms in marine sediments |
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 |