有压管道充水过程水力特性分析.docVIP

有压管道充水过程水力特性分析 摘 要：本文在圣维南方程基础上，采用假想窄缝法和虚拟流量法，并运用稳定性高的Priessmann隐格式建立有压的充水数学模型。该充水数学模型能模拟有压管道从无水变为有压流的水流运动过程，其中包括漫流和明满过渡流。文中详细阐述了实现虚拟流量的方法和过程，并结合具体工程分析充水过程中水位和流量的变化规律，为有压管道的充水运行提供科学依据。 关键词：有压管道，充水过程，虚拟流量，数值模拟，水力特性 有压管道输水方式具有占地面积少、输水损失小、维护费用低、使用寿命长、过流能力大等优点而往往为许多大型的调水工程所采用，如吉林引松入长、山西引黄入晋、南水北调等工程。但输水过程中水力控制引起的非恒定流态复杂，易形成明满过渡流[4]，给管道带来危害，在有压管道充水过程[1]中尤为突出。有压输水管道往往长达数公里，由于管道初始处于放空状态，流量持续充入后会在同一管道先后经历漫流[2]、无压流、明满交替流和有压流四种复杂的过渡流态。模拟有压管道的充水过程，不仅可以预测波前[3]的流达位置及时间，还可以掌握管道在明漫交替流水位及流量变化情况。因此，分析其水力特性的变化在工程应用上具有重要的意义。 充水过程可以采用一维明渠非恒定流控制方程来描述： 运动方程 （1） 连续方程 （2） 式中：为流量； 为水深； 为管道底坡；为水面宽；为过水断面面积；为湿周； 为曼宁粗糙系数，为有压流动的波速。 充水过程中管道经无压流过渡到有压流，而明渠非恒定流方程仅能计算无压流，因此可以采用Preissmann提出的假想窄缝法[]。即当水流接触管顶后，水面宽变为常数，这样，无压流、明满流和有压流就可以用统一的微分方程求解。 漫流在充水初期水流还未到达管道末端，不能将无水部分的水深和流量都设置为零，而认为自波前开始至末端有一股虚拟流动水流，使水流在迭代计算时平稳过渡到无压流阶段if(kk/=ns+1)then; if(kk.eq.ns)then; if(h(ns)hv)then; kk=ns+1; else; kk=ns; endif; else; do i=1,ns; if(h(i)hv.or.q(i)qv)then; kk=I; exit; endif; enddo; endif; endif kk为波前位置（如果波前到达最后一个断面ns，则令kk=ns+1，表明管道全部为无压流动，无需用到此段代码）；h、q分别为水深和流量；hv、qv分别为虚拟水深和虚拟流量。 3 工程计算实例 3.1 基本资料 南水北调天津干线将于建成通水，工程首次大规模使用4.4m管径，提高了输水流量，但增大单孔箱涵容水体积的同时，在复杂的水力过渡过程中，也增加了输水运行风险。首端堰后边界满足波动方程： （） 末端保水堰边界分以下两种情况考虑当，末端堰前的波动方程： （）当，末端堰前的波动方程： （） 式中：、、首端时段初与时段末的水深流量；、、末端时段初与时段末的水深流量；为首/末端底板高程；为恒定入流流量；虚拟流量；面积；末端堰前面积；过堰流量，为堰宽，为流量系数为堰顶高程。 4 计算结果及分析 数据采用模型实测的流量系数，虚拟流量取m3/s，虚拟水深取，首端以恒定的5m3/s流量流入，模拟10+679至14+000管道的充水过程.1 充水局部过程分析 充水局部过程分析可以水流在管道中的细微变化过程，能更全面的认识和了解充水过程中流态变化趋势。 图为漫流阶段Ⅰ、Ⅱ、Ⅲ三组沿管长方向分布的水面线族，每一组水面线之间间隔100s，分别从100s至500s、1000s至1400s、2400s至3000s。由图可明显看出，三组水面线间隔由稀至密分布，因为随着水流的波前越往前流动，水体在管长方向拉长，首端给水一定的条件下，涨高相同的距离需要更长的时间。第Ⅲ组的6条水面线显示了由漫流过渡到明渠流动的变化过程，在2400s时波前到达末端，当水流遇到堰的阻挡后，水面开始上升，此时进入到明渠流动阶段，并逐渐往前影响水流运动趋势，使之升高，在末端形成的“V”字形状，而高差二十公分的“V”字形状在实际长3.3多公里的管道中很难观察得到，体现出了数值模拟的优势。图清晰显示了明渠流动阶段水流的运动趋势。图a中各水面线间隔100s绘制，管内前后水面落差较大，进入的水流基本上用于提升管道下半段水面高程，缩小上下游水位差；图b中各水面线时间间隔200s，4800s时水面大致相平，充入的水使整个管内水面上升，开始进入明渠流的