发明名称 一种泥沙群体沉速的测量方法及测量装置
摘要 本发明涉及一种泥沙群体沉速测量方法采用的测量装置,属于河工模型试验领域,该装置包括沉降筒、多个测压设备和计算机三部分;沉降筒为垂直放置的任意截面形状的柱状容器,边壁沿垂向均匀布置多个测压孔,测压孔的个数不少于3个;测压设备垂直安装在测压孔中,并与边壁密封,测压设备的感压元件的中心与测压孔中心齐平;各测压设备与计算机相连,在计算机的控制下同步测量压强,并将测量结果实时传回计算机进行分析和存储。本方法通过测量压强的变化过程计算出泥沙浓度的变化过程,进而计算出泥沙的群体沉降速度。本发明能实现自动、无干扰、准确地测量含沙水体中泥沙群体沉降速度。
申请公布号 CN103529237B 申请公布日期 2016.06.29
申请号 CN201310499361.0 申请日期 2013.10.22
申请人 清华大学 发明人 陈启刚;李文杰;李丹勋;钟强;王兴奎
分类号 G01P5/00(2006.01)I;G01N15/04(2006.01)I 主分类号 G01P5/00(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 廖元秋
主权项 一种泥沙群体沉速测量方法,其特征在于,该方法采用的测量装置包括沉降筒、多个测压设备和计算机三部分;沉降筒为垂直放置的任意截面形状的柱状容器,边壁沿垂向均匀布置多个测压孔;测压设备垂直安装在测压孔中,并与边壁密封,测压设备的感压元件的中心与测压孔中心齐平;各测压设备与计算机相连;该方法包括以下步骤:1)将待测沙浆盛入沉降筒,静置沙浆待液面恢复平静后,通过计算机开启各测压设备开始按设定的频率同步测量浑水压强,各测压设备将实时测量的压强传回计算机;2)根据各测压设备测量的压强计算沙浆在任意时刻t的平均含沙量:设测压设备的个数为N,各测压设备(2)及所在的水平面分别编号为n,n=1,2,…N,水平面n距离沙浆液面的高差为H<sub>n</sub>,n=1,2,…N,测压设备n在任意时刻t的压强读数为P<sub>n</sub>(t),n=1,2,…N,各测压设备测量压强的频率为f=1/Δt,其中Δt为测压时间间隔;根据静水压强的定义,水平面n与液面之间的沙浆在任意时刻t的体积平均容重<img file="FDA0000895822900000011.GIF" wi="110" he="79" />为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><msub><mi>&gamma;</mi><mi>n</mi></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><msub><mi>H</mi><mi>n</mi></msub></mfrac><mo>,</mo><mrow><mo>(</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000895822900000012.GIF" wi="1220" he="151" /></maths>根据质量含沙量的定义,水平面n与液面之间的沙浆的平均含沙量<img file="FDA0000895822900000013.GIF" wi="111" he="78" />与体积平均容重<img file="FDA0000895822900000014.GIF" wi="113" he="77" />之间满足如下关系:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><msub><mi>S</mi><mi>n</mi></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&gamma;</mi><mi>s</mi></msub><mfrac><mrow><mover><msub><mi>&gamma;</mi><mi>n</mi></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>&gamma;</mi></mrow><mrow><msub><mi>&gamma;</mi><mi>s</mi></msub><mo>-</mo><mi>&gamma;</mi></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000895822900000015.GIF" wi="1262" he="158" /></maths>式中γ<sub>s</sub>为待测泥沙的干容重,γ为水的容重;3)根据式(2)求得所有水平面与液面之间的沙浆的平均含沙量,得到水平面n与水平面n‑1之间的沙浆在时刻t的平均含沙量为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><msub><mi>S</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi></mrow></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>H</mi><mi>n</mi></msub><mover><msub><mi>S</mi><mi>n</mi></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mover><msub><mi>S</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>H</mi><mi>n</mi></msub><mo>-</mo><msub><mi>H</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000895822900000016.GIF" wi="1358" he="151" /></maths>将<img file="FDA0000895822900000017.GIF" wi="177" he="85" />取为与液面高差为(H<sub>n</sub>+H<sub>n‑1</sub>)/2的平面处的含沙量,插值得到水平面n处的含沙量S<sub>n</sub>(t);<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>S</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mover><msub><mi>S</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi></mrow></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mover><msub><mi>S</mi><mrow><mi>n</mi><mo>,</mo><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>,</mo><mrow><mo>(</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000895822900000021.GIF" wi="1230" he="142" /></maths>4)根据输沙率的定义,时刻t‑Δt至t+Δt之间通过水平面n的泥沙量除以沉降筒截面面积A及时间即为输沙率g<sub>n</sub>(t),计算如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>g</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>(</mo><mover><msub><mi>S</mi><mi>n</mi></msub><mo>&OverBar;</mo></mover><mo>(</mo><mi>t</mi><mo>+</mo><mi>&Delta;</mi><mi>t</mi><mo>)</mo><mo>-</mo><mover><msub><mi>S</mi><mi>n</mi></msub><mo>&OverBar;</mo></mover><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;</mi><mi>t</mi><mo>)</mo><mo>)</mo><msub><mi>H</mi><mi>n</mi></msub><mi>A</mi></mrow><mrow><mn>2</mn><mi>A</mi><mi>&Delta;</mi><mi>t</mi></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000895822900000022.GIF" wi="1342" he="159" /></maths>将g<sub>n</sub>(t)取为水平面n在时刻t的瞬时输沙率;5)根据式(4)和式(5)得到任意时刻t距液面高差为H<sub>n</sub>处的泥沙的群体沉降速度ω<sub>n</sub>(t)为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>g</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>S</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000895822900000023.GIF" wi="1157" he="175" /></maths>
地址 100084 北京市海淀区清华园1号