随着中国经济的快速发展,管道运输在国民经济、国防工业生产和人们的日常生活中,正发挥着越来越重要的作用,如石油,天然气,自来水管已成为国家经济发展的命脉。这些管道在建设过程中管道弯曲或因机械压力导致管道变形致使管道堵塞的情况发生。在寒冷的环境里油气输运管道会形成天然气水合物,造成管道堵塞进而带来经济损失。由于输运管道的介质多种多样,加之管道所处环境十分复杂,以及管道堵塞形式的不确定性,因此管道堵塞的检测非常困难。所以为了解决管道堵塞的检测问题寻找一种简单可行同时高效可靠的方法很重要。本文通过对堵塞管道情况的分析,建立水击方程的数学模型分析管道堵塞后水力瞬态流动的变化与特性,为以后的管道探堵工作提供理论基础。本研究利用水击方程建立管道堵塞的数学模型,引入了基于有限体积法的Godunov格式,再用Fortran对瞬态流动模型方程进行编程得到计算求解程序,得到有压管道在瞬变流条件下指定位置的压力实时值。并与ADINA数值模拟二维模型的结果进行了对比,验证了结果的正确性。通过改变计算网格的大小说明了本文计算结果与网格精度无关。考察了水锤波的传播过程几个重要时间点,通过数值比较说明本文的数值模拟方法结果与Wylie(1993)的一维水锤计算理论的结果十分接近,验证了本研究方法可以计算管道中波传播与反射的复杂特性。在计算过程中考虑了非定常摩擦力对管道流动的影响,在一维水锤控制方程中引入扩散项替代非定常摩擦力所引起的粘性作用,通过比较不同堵塞长度时有无扩散项的压强曲线,说明对堵塞长度越短的管道这种粘性作用越重要,进一步证明了在非定常管道流动中引入扩散项近似非定常摩擦力的影响是有必要的。通过改变变截面长度观察堵塞长度对水力瞬变的影响,可以看出在变截面为10m和5m时结果与Wylie(1993)的一维特征线法的结果很接近。而当变截面长度越小时Wylie的一维特征线法却无法得到实际的结果。研究发现在管道流动中瞬态压力波峰面有厚度,而且这个厚度对波反射的峰值影响很大。通过比较堵塞长度为10m和1m的传播过程得出堵塞长度越短反射量越小的结论。