冲积河流悬移质泥沙与床沙交换机理及计算方法研究 .docx

工程建筑学相关资料ENGINEERING ARCHITECTURE RELATED INFORMATION 工程建筑学相关资料 ENGINEERING ARCHITECTURE RELATED INFORMATION PAGE PAGE 1 冲积河流悬移质泥沙与床沙交换机理及计算方法研究  摘要：本文在分析前人研究成果的基础上，从非恒定挟沙水流沙量平衡入手，建立了冲积河流悬移质泥沙与床沙交换关系表达式，并引入平衡冲淤物粒径的概念，对冲淤物粒径的计算方法进行了系统分析，提出了一整套计算非恒定挟沙水流悬移质泥沙与床沙交换的模式与方法，同时将该方法引用到数学模型中，采用黄河下游典型洪水进行了验证。 关键词：悬移质泥沙 床沙 交换机理 平衡冲淤物粒径 1 引言　　悬移质泥沙与床沙交换机理及计算方法，作为河流泥沙动力学的重要问题之一， 长期以来，被许多学者所关注，并取得了卓有成效的进展。该问题的研究，对河道演变分析、河工模型设计、河流泥沙数学模型计算以及河流泥沙动力学学科的发展都颇有意义。　　现有方法一般为将悬移质泥沙和床沙分成若干组，通过计算分组挟沙力实现悬沙与床沙的交换，从而求出悬沙与床沙级配的变化规律。其中韩其为是由悬移质级配求分组挟沙力，该方法假定含沙量级配与挟沙力级配相同，且各粒径组泥沙独自保持平衡而不相互影响，在此基础上建立了淤积和冲刷两种情况的级配计算公式[1]，但在计算过程中出现与实际河流冲淤不相符合的情况，如河床发生冲刷时，各粒径组的泥沙均发生冲刷，而淤积时各粒径组泥沙均发生淤积。HEC-6模型[2]是以床沙级配求分组挟沙力，其基本思路为：先求每一粒径组泥沙的可能挟沙力，即这一粒径组含沙量占100%的水流挟沙力，而实际挟沙力则是此粒径组泥沙在床沙中所占百分比与上述可能水流挟沙力的乘积，这种方法的前提假定与实际挟沙图形有较大出入。　　李义天[3]通过理论分析建立了输沙平衡状态下床沙质级配与床沙级配之间的关系为(1)(2)　　式中Pi， Pbi分别为床沙质与床沙级配；ωi为第i粒径组泥沙沉速；u*为摩阻流速；κ为卡门常数；m为粒径分组数；σv为垂向紊动强度，原作者取σv=u*。在已知床沙级配的条件下，根据式(1)计算出悬沙中床沙质级配后，再由全沙挟沙力乘以挟沙力级配Pi即可得出分组挟沙力，进一步采用非饱和输沙模型计算河段内的冲淤，尔后根据冲淤结果调整床沙级配作为下一时段的计算依据。　　钱意颖、曲少军、张启卫、韦直林[4]提出的分组挟沙力计算方法在数学模型计算中也颇有应用价值。但应该承认，水沙运行过程中，各粒径组泥沙是相互影响的，人为地分组后分别计算显然与实际出入较大。特别象黄河这样的多沙细沙河流，细颗粒泥沙对水流流变特性及挟沙特性影响尤为显著。此外，现有通过计算分组挟沙力来确定泥沙级配的方法多应用于恒定流模型，而天然河流一般为非恒定流，泥沙级配不仅沿程变化，而且还随时间变化，故此进一步研究悬沙与床沙交换机理有重要意义。2 冲积河流床沙与悬沙交换模式及计算方法初步研究2.1悬沙与床沙交换过程中悬沙粒径变化规律研究　　如图1所示，取1-2之间的δx微段，Q为流量；S为含沙量；A为过水面积；Ψ(d)为悬移质泥沙粒径分布函数；δt为时间微段；ρ′为冲淤物干容重；A0为横断面淤积面积；Ψc(d)为冲淤物粒径分布函数；qL 、SL、ΨL(d)分别为侧向入流单宽流量、含沙量、入流悬沙粒径分布函数。则δt时段内，δx河段上di的净流入量P1为(3)　　微段δx的水体沙量增加中粒径为di的泥沙所占的量P2为(4)　　侧向入流与冲淤使di的增加量P3为(5)　　由质量守恒原理，则P1+P3=P2，将式(3)～(5)代入，两边除以δxδt，并忽略微小量，整理得(6)?式(6)两边乘以di，并对d积分后化简，同时将悬沙平均粒径、冲淤物平均粒径、侧向入流泥沙平均粒径代入可得?(7)　　根据多元函数求导法则将上式展开，并代入泥沙连续方程，进一步化简得(8)　　采用式(7)或式(8)即可计算任一时刻任一断面的悬沙平均粒径大小。对式(6)两边乘以\$d2i\$，积分后并取悬沙粒径二阶圆心距；冲淤物粒径二阶圆心距；侧向入流泥沙粒径二阶圆心距，则化简可得?(9)或(10)　　求得进口ξd过程后，即可依据式(9)或(10)计算出任一时刻任一段面的悬沙粒径二阶圆心距。2.2 床沙与悬沙交换过程中床沙粒径变化规律研究　　取单位宽度的断面，设δt时段床沙交换层厚度为Hc；河床高程为Z；床沙粒径分布函数为Ψ(Di)；床沙平均粒径为Dcp。则对于Di的泥沙含量，瞬初x1=HcΨ(Di);瞬末；冲淤部分。由质量守恒知x2-x1=x3，同时等号两边同除以δt，可化简为?(11) 式两边乘以Di并对D积分得上 (12)根据多元函数求导