一种泵站引渠及前池水体泥沙参考浓度的确定方法技术

本发明专利技术涉及一种泵站引渠及前池水体泥沙参考浓度的确定方法，包括以下步骤：1)计算得出泵站引渠及前池水体的流速随时间变化的关系式；2)计算泵站引渠及前池水体的最大谢尔兹数Θ

Method for determining reference sediment concentration of water diversion channel and front pool of pumping station

The invention relates to a method for determining the diversion canal and water pumping station forebay sediment reference concentration, which comprises the following steps: 1) relationship between the calculated and the diversion canal pumping station forebay water flow velocity changes with time; 2) calculation and water diversion canal pumping station forebay of the largest number of shields.

【技术实现步骤摘要】
一种泵站引渠及前池水体泥沙参考浓度的确定方法
本专利技术涉及一种泵站引渠及前池水体泥沙参考浓度的确定方法，属于水利工程泵站领域。
技术介绍
在泥沙设计的数值模拟中，很关键的一步是给出泥沙浓度方程的底部边界条件，一般都在距离壁面2倍泥沙粒径高度的位置提出参考浓度C。边界条件的精度决定了数值模拟的结果精度。精准知晓包含相位差作用和考虑泥沙体积守恒的泥沙参考浓度的大小，是准确模拟非恒定流动条件下水体挟沙率和输沙率、判断水泵机组磨蚀的重要基础。现有工程中广泛使用的参考浓度都是基于恒定流动的理论推导出来的，不包含相位差作用即泥沙运动相对于水动力条件所需要的时间响应，也未能考虑泥沙体积的守恒。这导致出现的问题是：无法体现泥沙运动相对流速的相位漂移和运动的泥沙总量在流动速度减少后的相位残留，给计算带来极大的误差；泥沙浓度随着流动条件的增大也难以达到饱和状态，参考浓度随流动强度太大从而失效。
技术实现思路
针对上述问题，本专利技术的目的是提供一种在非恒定流动条件下确定相关相位差参数和泥沙体积守恒定律来获得泵站引渠及前池水体泥沙参考浓度的方法。为实现上述目的，本专利技术采用以下技术方案：一种泵站引渠及前池水体泥沙参考浓度的确定方法，包括以下步骤：1)计算得出泵站引渠及前池水体的流速随时间变化的关系式U(t)；2)计算泵站引渠及前池水体的最大谢尔兹数Θm、摩擦因子f和粗糙高度kN；3)计算泵站引渠及前池水体的泥沙沉降速度和泥沙在输沙层内的沉降时间与流动周期的比值Ψ；4)计算泵站引渠及前池水体的泥沙对流速的相位漂移ψ和泥沙相位的残留α；5)计算泵站引渠及前池水体的泥沙体积分数参考浓度。所述步骤1)中，这里采用2阶Stokes波动理论得出泵站引渠及前池水体的流速随时间变化的关系式。所述步骤2)中，采用如下的式(1)～(3)计算泵站引渠及前池泥沙的最大谢尔兹数Θm、摩擦因子f和粗糙高度kN，具体过程如下：其中U是泵站引渠及前池水体的流速，下标m代表最大值；D是泥沙粒径；T是流动周期；Θ是谢尔兹数，下标m代表最大值；f是摩擦因子；s是泥沙与水密度的比值；所在地的重力加速度g；kN是动床面粗糙高度；先假定最大谢尔兹数Θm<1，由式(3)得到粗糙高度的值kN，代入式(2)，得到摩擦因子的值f；摩擦因子的值代入式(1)得到的最大谢尔兹数的值Θm若大于1则假定不成立；根据最大谢尔兹数Θm>1，由式(3)得到粗糙高度kN＝5ΘmD代入式(2)，并联合式(1)得到从式(4)解得摩擦因子的值f，代入式(1)得最大谢尔兹数的值Θm，代入式(3)得粗糙高度的值kN。所述步骤3)中，泵站引渠及前池水体的泥沙沉降速度和泥沙在输沙层内的沉降时间与流动周期的比值Ψ，具体计算过程如下：从如下的式(5)得到泥沙沉降速度的值，代入如下的式(6)得到泥沙在输沙层内的沉降时间与流动周期的比值Ψ；其中ν是水的运动粘性系数。所述步骤4)中，泵站引渠及前池水体的泥沙对流速的相位漂移ψ和泥沙相位的残留α，具体过程如下：采用如下提供的式(7)和式(8)，计算泥沙对流速的相位漂移ψ和泥沙相位的残留αα＝exp(-0.2/Ψ)。(8)所述步骤5)中，泵站引渠及前池泥沙的泥沙体积分数参考浓度的获取方法如下：将所述步骤3)和所述步骤4)获得的泥沙对流速的相位漂移ψ、泥沙相位的残留α和泥沙沉降速度w代入下述式(9)得到泥沙体积分数参考浓度随时间变化的关系式；将式(1)得到的谢尔兹数代入式(9)，即可得到最终的泥沙参考浓度随时间变化和相应的函数图像，式(9)如下所示：其中F代表水体的挟沙能力。本专利技术由于采取以上技术方案，其具有以下优点：1、通过泥沙体积守恒的考虑，避免参考浓度随流动强度太大而失效。2、通过所确定的相关相位差参数进行设计计算，可以反应非恒定流动过程中泥沙相对水动力条件所需要的时间响应。3、本专利技术通过所确定的相关相位差参数进行设计计算，体现了泥沙运动相对流速的相位漂移。4、本专利技术通过所确定的相关相位差参数进行设计计算，体现了具有运动的泥沙总量在流动速度减少后的相位残留。附图说明图1本专利技术工作过程采用的非恒定流动速度过程，U/Um是无量纲的流速，t/T是无量纲的时间；图2是本专利技术与恒定流理论确定的参考浓度比较，流动强度Um＝0.6m/s，C是泥沙参考浓度，实线代表本专利技术，点划线代表传统方法；图3是本专利技术与恒定流理论确定的参考浓度比较，流动强度Um＝1.0m/s，实线代表本专利技术，点划线代表传统方法；图4是本专利技术最大泥沙体积参考浓度与流动强度的比较，Θm是最大谢尔兹数，Cm是最大泥沙体积参考浓度。具体实施方式下面结合附图和实施例对本专利技术进行详细的描述。本专利技术提供的一种泵站引渠及前池水体泥沙参考浓度的确定方法，包括以下步骤：1)计算得出泵站引渠及前池水体的流速随时间变化的关系式U(t)，具体可以利用2阶Stokes波动理论得出泵站引渠及前池水体流速随时间变化的关系式，因该部分属于本领域的公知方法，因此不进行具体说明。2)计算泵站引渠及前池水体的最大谢尔兹数Θm、摩擦因子f和粗糙高度kN，具体可通过下述的式(1)～(3)进行计算：式中，U是泵站引渠及前池水体的流速，下标m代表最大值；D是泥沙粒径；T是流动周期；Θ是谢尔兹数，下标m代表最大值；f是摩擦因子；s是泥沙与水密度的比值；g是重力加速度；kN是动床面粗糙高度。在具体求解时，先假定最大谢尔兹数Θm<1，由式(3)得到粗糙高度的值kN，代入式(2)，得到摩擦因子的值f；摩擦因子的值代入式(1)得到的最大谢尔兹数的值Θm若大于1则假定不成立。根据最大谢尔兹数Θm>1，由式(3)得到粗糙高度kN＝5ΘmD代入式(2)，并联合式(1)得到从式(4)解得摩擦因子的值f，代入式(1)得最大谢尔兹数的值Θm，代入式(3)得粗糙高度的值kN。3)计算泵站引渠及前池水体的泥沙沉降速度和泥沙在输沙层内的沉降时间与流动周期的比值Ψ，具体计算过程如下：由下述式(5)计算得到泵站前池水体的泥沙沉降速度w，将w的值代入式下述式(6)得泥沙在输沙层内的沉降时间与流动周期的比值Ψ，运用的式(5)和式(6)如下：式中，ν是水的运动粘性系数。4)计算泵站引渠及前池水体的泥沙对流速的相位漂移ψ和泥沙相位的残留α，具体过程如下：利用下述式(7)得到泵站前池的泥沙对流速的相位漂移ψ，由下述式(8)得泥沙相位的残留α的值：α＝exp(-0.2/Ψ)。(8)5)计算泵站引渠及前池水体的泥沙体积分数参考浓度C(t)。将谢尔兹数Θ、泥沙对流速的相位漂移ψ、泥沙相位的残留α和泥沙沉降速度w代入下述式(9)得到泥沙体积分数参考浓度随时间变化的关系式：其中F代表水体的挟沙能力。下面结合附图和实例对本专利技术进行详细的描述。以下实例用于说明本专利技术，但不用来限制本专利技术的范围。1)计算得出泵站引渠及前池水体的流速随时间变化的关系式。已知泥沙所处气压为一个标准大气压、水温为20℃、最大流动速度Um＝0.6m/s，流动周期T＝5s，泥沙粒径D＝1.0×10-4m，泥沙与水密度的比值s＝2.65，所在地的重力加速度g＝9.8m/s2，水的运动粘性系数ν＝1.0×10-6m2/s，临界谢尔兹数Θcr取0.05。采用2阶Stoke波动理论得出流速与时间的关系式为：U(t)＝0.48本文档来自技高网...

【技术保护点】
一种泵站引渠及前池水体泥沙参考浓度的确定方法，包括以下步骤：1)计算得出泵站引渠及前池水体的流速随时间变化的关系式U(t)；2)计算泵站引渠及前池水体的最大谢尔兹数Θ

【技术特征摘要】
1.一种泵站引渠及前池水体泥沙参考浓度的确定方法，包括以下步骤：1)计算得出泵站引渠及前池水体的流速随时间变化的关系式U(t)；2)计算泵站引渠及前池水体的最大谢尔兹数Θm、摩擦因子f和粗糙高度kN；3)计算泵站引渠及前池水体的泥沙沉降速度和泥沙在输沙层内的沉降时间与流动周期的比值Ψ；4)计算泵站引渠及前池水体的泥沙对流速的相位漂移ψ和泥沙相位的残留α；5)计算泵站引渠及前池水体的泥沙体积分数参考浓度C(t)。2.如权利要求1所述的一种泵站引渠及前池水体泥沙参考浓度的确定方法，其特征在于：所述步骤1)中，这里采用2阶Stokes波动理论得出泵站引渠及前池水体的流速随时间变化的关系式。3.如权利要求2所述的一种泵站引渠及前池水体泥沙参考浓度的确定方法，其特征在于：所述步骤2)中，采用如下的式(1)～(3)计算泵站引渠及前池水体的最大谢尔兹数Θm、摩擦因子f和粗糙高度kN，具体过程如下：其中U是泵站引渠及前池水体的流速，下标m代表最大值；D是泥沙粒径；T是流动周期；Θ是谢尔兹数，下标m代表最大值；f是摩擦因子；s是泥沙与水密度的比值；所在地的重力加速度g；kN是动床面粗糙高度；先假定最大谢尔兹数Θm<1，由式(3)得到粗糙高度的值kN，代入式(2)，得到摩擦因子的值f；摩擦因子的值代入式(1)得到的最大谢尔兹数的值Θm若大于1则假定不成立；根据最大谢尔兹数Θm>1，由式(3)得到粗糙高度kN＝5ΘmD代入式(2)，并联合式(1)得到从式(4)解得摩擦因子的值f，代入式(1)得最大谢尔兹数的值Θm，代入式(...

【专利技术属性】
技术研发人员：陈鑫，王福军，
申请(专利权)人：中国农业大学，
类型：发明
国别省市：北京,11