位山灌区输沙渠泥沙运移特性研究.docxVIP

? ? 位山灌区输沙渠泥沙运移特性研究 ? ? 马 强，徐立荣，董晓知，徐 晶 (济南大学 水利与环境学院，山东 济南 250022) 黄河是世界上输沙量最大和含沙量最大的河流，因此引黄灌区具有“引黄必引沙”的重要特点[1]。引黄灌区自开始引水以来，泥沙问题一直伴随着灌区发展[2]。虽然自2002年小浪底水库进行调水调沙试验，并转入正式生产运行后，黄河下游特别是山东境内的引黄条件发生了很大变化，如引黄能力下降，输沙量变小等，但是引黄灌渠排沙问题仍未减轻，如何排沙依旧是黄河下游治理的重中之重，需要重新提出科学的水沙调度和合理的工程措施，从而延长沉沙池使用寿命，减少骨干渠道淤积，更好地发挥灌区的灌溉效益。 渠道是灌区灌排系统组成的重要部分，而渠道输水输沙能力的有效实现是整个灌区正常运行的关键[3-4]。为了解决灌区的泥沙问题，提高渠道的输水输沙能力，减少泥沙淤积，保证灌区的正常运行，人们通常以控制渠道泥沙淤积为目标，建立相应水沙数学模型。一般对河道泥沙淤积的数学模型研究多采用一维水沙模型，模拟断面各水沙要素的总平均值[5]，相比之下，平面二维水沙模型能更清楚地表达河床的平面变形，因此被广泛应用于河道水沙运动研究[6-7]。 目前，对于引黄灌渠输沙渠道二维水沙的泥沙淤积研究成果较少。本文中以山东省聊城市位山灌区一支改造后的输沙干渠——东输沙渠作为研究对象，以控制输沙渠渠道泥沙淤积为目标，利用水力学模型MIKE 21具有的简单、快速、精确、可视性好等优点[8]，构建输沙渠二维水沙模型，在经过模型率定与验证后，设计典型工况对水流特性及泥沙淤积分布规律进行模拟预测，筛选保水输沙的最佳方案，为黄河下游引黄灌区干渠泥沙治理和利用提供参考。 1 研究区概况 山东省聊城市位山灌区地处黄河下游，是全国第五大灌区、山东省最大灌区[9-10]。灌区共包括四区一市五县，土地总面积为5 380 km2，总耕地面积为3 793 km2，设计灌溉面积为3 387 km2[11]。灌区南靠黄河，北临卫运河，东与山东省德州市相接，马颊河和徒骇河从灌区自西南向东北贯穿而过。灌区属黄河冲积形成的平原，地势平坦开阔，走向略微向东北方向倾斜[9，11]。 东输沙渠于2019年改造完成，全长14.5 km，担负着一干渠供水供沙的任务，具体位置如图1所示。东输沙渠下接进入东沉沙区，有沉沙池3个，已还耕1个，主要使用1#、3#池轮用，2个沉沙池长度分别为3.0、3.2 km。渠道断面被改造成梯形复式断面，为了使水流顺畅、流态稳定，渠道中心线保持不变。渠道比降仍采用原比降1/6 000，边坡坡比为1∶2.0，复式断面底宽为16.0 m，渠道开挖深度为1.5 m，浅滩宽度为2 m。 位山灌区地理位置地图从国家测绘地理信息局标准地图服务网站下载，地图审批编号为GS(2019)3333(/browse.html?picId=%224o28b0625501ad2bfc0210%22)，经过ArcGIS10.3软件数字化处理得到。图1 位山灌区输沙干渠图 2 二维水沙模型建立及验证 2.1 水动力模型理论 水力学模型MIKE 21中的水动力(HD)模块是必需选择的核心模块，是所有模拟的基础。该模块建立在基于二维数值求解方法的浅水方程基础上，在深度上集成了三向不可压缩雷诺(Reynolds)平均Navier-Stokes方程，并服从Boussinesq假设和静水压假设[12-15]。 二维非恒定浅水方程组包括连续性方程，X、Y方向动量方程，即 连续性方程 (1) X方向动量方程 (2) Y方向动量方程 (3) 该模型是由动量守恒方程和水流连续方程组成，在水平面上可以使用笛卡尔坐标。 2.2 泥沙模型理论 水力学模型MIKE 21中的输泥模块(MT)结合了多粒径级和底床分层，描述了黏聚性泥沙(淤泥或黏土)在波浪和水流作用下的冲刷、传输和沉积[16]。 1)基本方程。二维输泥模块的对流扩散控制方程为 (4) 2)泥沙沉积过程。泥沙在水体中通常处于沉积和悬浮状态，泥沙淤积是泥沙颗粒从水体沉降至床面的输移过程。当床面切应力小于淤积临界切应力时，产生泥沙淤积；反之，床底发生冲刷。细颗粒泥沙的沉降速度取决于颗粒大小、温度、泥沙浓度等。黏性泥沙的淤积的控制方程为 D=ωsρbPD， (5) 式中：D为沉积率；ωs为泥沙沉降速度，m/s；ρb为泥沙的近底质量浓度，kg/m3；PD为泥沙沉降在河床的概率函数， (6) 其中τb、τcd分别为河床剪切应力、河床淤积临界剪切应力，N/m2。 3)河床的冲刷过程。河床剪切应力τb大于临界冲刷剪切应力τce时，就会发生河床层的冲刷。床底冲刷程度表达式为 (7) 式中：SE为侵蚀率；E为床底冲刷程度，kg/(s·m2)；n为侵蚀程度。 2.3