小浪底水库坝区扰动泥沙运行规律研究与扰沙方案计算.pdfVIP

第十二届中国海岸工程学术讨论会论文集 173 小浪底水库坝区扰动泥沙运行规律 研究及扰沙方案计算 茹玉英1，张俊华1，陈大宏2，屈孟浩1 (1．黄河水利科学研究院，河南郑州450003：2．武汉大学，湖北武汉430072) 摘要：采用实体模型试验资料及三维数学模型，对小浪底水库坝区流场进行了分析计算，探讨了坝区扰动泥沙的运动规律， 给出了运动方程，对各种扰动方案及排沙效果进行了计算分析，为2004年黄河上进行的大空间尺度的“基于人工扰动的调 水调沙试验”的实施方案提供参考。 关键词：小浪底；坝区；扰动泥沙；排沙比：流场；实体模型 1小浪底坝区流场分布分析计算 1．1实体模型的坝区流场分布分析 20世纪80年代，为了研究小浪底枢纽建成投入运用后，坝区的水流流态和淤积形态，黄河水利科学 研究院、南京水利科学研究院和中国水利科学研究院，根据各自的条件，分别进行了小浪底枢纽工程整体 悬沙模型试验，黄科院还进行了概化模型试验。 其中黄科院和中国水科院进行的整体模型试验的几何比尺是相同的，彳，=彳H=100，南京水利科学 研究院进行整体模型试验几何比尺为五，=名Ⅳ=80。 黄科院的概化模型在宽lm，高1m，长20m的水槽内进行，概化模型的几何比尺与黄科院的整体模 型相同彳，=五Ⅳ=100，模型沙亦为迥料沙，泄水排沙洞的尺寸、高程亦与整体模型相同。 通过对黄科院模型试验资料的分析，得到以下认识： 000 1)在水位较高的条件下，若上游来水流量很小，仅靠排沙洞和泄洪洞泄流，在总泄量接近3 m3／s 时，坝前流速可达1．0m／s以上，但沿程衰减很快。在距坝100m处，180m高程的流速仅0．4m／s，在 000 距坝l m处，180111高程的流速接近于零，说明在小浪底枢纽具体运用条件下，底孔吸流长度很短， 不可能超过1000m。 2)小浪底枢纽在坝前水位240 m条件下，打开排沙洞泄水(Q600m3／s)，有可能在坝前形成一个 长约200～500m的冲刷漏斗，即通过排沙底孔的吸流作用，可以把淤在坝前至200,500m区域内的泥沙 排出库外。 1．2三维数模计算坝区流场分布 1)数学模型及模型参数 流体的紊流运动可以用雷诺平均方程组(RANS)来描述： V．U=0 (1) (2) Ot j 上述方程中U为流体速度、P为压力、p为流体密度、“为流体动力黏性系数、m为紊流黏性系数，6t。是克 罗内克尔符号，k为紊动能。计算分析采用了标准k-e模型： k+V．(pUk)-V．[c肛+等，V七]=P—p￡ (3) (4) 去p￡+V·cpu￡，一V·[cp+等，V￡]=c。·昙P-C2P譬 其中： (5) 174 第卜二届中国海岸工程学术讨论会论文集 fir=C∥p等 (6) 上述方程中￡为紊动能耗散率，6t为紊动能普朗特数，o￡为紊动能耗散率普朗特数。标准k-E模型参数 在计算中采用： C。=1．44，C2=1．92，C。，=0．09，仃t=1．0，o。=1．3 上述偏微分方程组采用有限体积法进行离散。 2)流场分布