燃气管网水力计算方法

现代燃气工程结课论文题目：燃气管网水力计算姓名：王朋飞学号：S20090146教师：范慧方引言随着能源结构的不断改变，燃气开发规模和应用规模的不断扩大。城市燃气管网是现代化城市人民生活和工业生产的一种主要能源配送方式，燃气输配管网的设计和运行要求对系统进行水力计算，获取必要的参数。燃气输配管网系统由高度整体化的管网所组成，在系统内燃气压力和流量变化很大，需要通过水力计算来确定管网中每一管段的尺寸（如管径、管径）、材质等参数以及压缩机的台数功率以保证既向用户合理地供应天然气，又能降低操作管理费用。1同时，考虑在满足用户用气量的前提下，当某一条或几条管道的使用有一定的压力要求时，水力计算数据可确定在这种最大承受压力下管道各个节点的压力，从而保证管网的正常运行。另外，水力计算也用于调整各个调压阀的出口压力来适应事故工况下输送压力的要求。随着燃气事业的发展，燃气输配管网系统也日趋庞大和复杂，为了掌握燃气在管道内的运行规律，合理地确定管道系统的设计和改造方案，保证管道系统的优化运行，提高管道系统的调度管理水平，解决管网流动的动态特性，在一些比较大型的城市燃气管网的水力计算分析中，必须要依靠相关的计算分析软件进行，以减少手工量和人工误差。1燃气管网水力计算燃气是可压缩流体，一般情况下管道内燃气的流动是不稳定流，由压送机站开动压缩机不同台数的工况以及用户用气量变化的工况，这些因素都导致了燃气管道内燃气压力和流量的变化。管内燃气沿程压力下降会引起燃气密度的减小。但是在低压管道中燃气密度变化可以忽略不计。所以，除了单位时间内输气量波动大的超高压天然气长输管线要用不稳定流进行计算外，在大多数情况下，设计燃气管道时都将燃气流动按稳定流计算。此外，很多情况下，燃气管道内的流动可认为是等温的，其温度等于埋管周围土壤的温度。燃气管网按照敷设形式可分为两大类：枝状管网和环状管网。下面就分别介绍两种形式的管网的水力计算特点和方法。1.1枝状管网水力计算枝状管网水力计算特点枝状管网是由输气管段和节点组成。任何形状的枝状管网，其管段数P和节点数m的关系均符合：燃气在枝状管网中从气源至各节点只有一个固定流向，输送至某管段的燃气图1枝状管网流量分配只能由一条管道供气，流量分配方案也是唯一的，枝状管道的转输流量只有一个数值，任意管段的流量等于该管段以后（顺气流方向）所有节点流量之和，因此每一管段只有唯一的流量值，如图1所示。管段3-4的流量为：Q34二q4757879710管段4-8的流量为：此外，枝状管网中变更某一管段的直径时，不影响管段的流量分配，只导致管道终点压力的改变。因此，枝状管网水力计算中各管段只有直径d与压力降厶P两个未知数枝状管网水力计算步骤对管网的节点和管段编号。确定气流方向，从主干线末梢的节点开始，利用Qi=0的关系，求得管网各管段的计算流量。根据确定的允许压力降，计算管线单位长度的允许压力降。根据管段的计算流量及单位长度允许压力降预选管径。根据计算选定的标准管径，求摩擦阻力损失和局部阻力损失，计算总的压力降。检查计算结果。若总的压力降超出允许的精度范围，则适当变动管径，直至总压力降小于并趋近于允许值为止。摩擦阻力损失，局部阻力损失和附加压头的计算方法一、摩擦阻力损失的计算方法根据城镇燃气设计规范(GB50028-93,2003版)附录A燃气管道摩擦阻力计算：1、低压燃气管道根据燃气在管道中不同的运动状态，其单位长度的摩擦阻力损失采用下列各式计算：(1)层流状态：Re_210064扎=ReP=1.131010Q4?Tld4T0Re2100（3）湍流状态:Re35001）钢管:K68.0.25花）=6.9106（K192.2）025Q5dQd5T。2）铸铁管:：.-0.102236（15158）0284dQ-6.4106（丄5158）0284QdQdT02、中压燃气管道根据燃气在管道中不同的材质，其单位长度的摩擦阻力损失采用下列各式计算:1）钢管:一0.1時RV25dRe2）铸铁管:222子九4109吟192宠）。.25存TT0临界状态：Re=21003500=0.03565Re-10P611.8Q7104dQ2T23Q-105d、1.9106（155l23Q-105dd5T0二0.102236（丄5158）0284dQT。根据聚乙烯燃气管道工程技术规程（CJJ63-95）知，聚乙烯燃气管道单位长度的摩擦阻力计算和钢管公式一样,只是K=0.01.而钢管K=0.15.二、局部阻力损失的计算方法当燃气流经三通、弯头、变径管、阀门等管道附件时，由于几何边界的急剧改变，燃气流线的变化，必然产生额外的压力损失，称之为局部阻力损失。在进行城市燃气管网的水力计算时，管网的局部阻力损失一般不逐项计算，可按然气管道摩擦阻力损失的5%-10%进行估算。对于庭院管和室内管道及厂、站区域的燃气管道，由于管路附件较多，局部阻力损失所占比例较大，常需逐一计算。局部阻力损失，可用下式求得：厂u2.P=Z2式中P局部阻力的压力损失（Pa）；x计算管段中局部阻力系数的总和；u管段中燃气流速（m/s）；燃气的密度（kg/m3）。局部阻力损失也可用当量长度来计算，各种管件折成相同管径管段的当量长度L2，各种管件当量长度L2查燃气热力工程常用数据手册，实际工程中通常按当量长度计算局部阻力。三、附加压头的计算方法由于燃气与空气的密度不同，当管段始末段存在标高差值时，在燃气管道中将产生附加压头，其值由下式确定：卩二g（a-g）H式中P附加压头（Pa）；g重力加速度；打空气的密度（kg/Nm3）；订燃气的密度（kg/Nm3）；H管段终端和始端的标高差值（m）。计算室内燃气管道及地面标高变化相当大的室外或厂区的低压燃气管道，应考虑附加压头1.2环状管网的水力计算环状管网的水力计算特点环状管网是由一些封闭成环的输气管段与节点组成。任何形式的环状管网，其管段数P、节点数m和环数n的关系均符合下式：P=mn-1环状管网任何一个节点均可由两向或多向供气，输送至某管段的燃气同时可由一条或多条管道供气，可以有不同的流量分配方案。此外，环状管网中变更某一管段的直径时，就会引起管段流量的重新分配，并改变管段各节点的压力值。因此，环状管网水力计算中各管段有三个未知量：直径di与压力降=R以及流量Qi，这与之前枝状管网的流动特点不一样，因而水力计算也不一样。水力计算基本方程1、连续方程气体被视为连续介质，在气体流动时总是连续地充满它所占据的空间，不出现空隙。这样根据质量守恒定律，流入节点的流体质量必然等于流出节点的流体质量，又由于假设气体稳定流动，所以体积流量也相等。也可以说连接于任何节点的所有管段流量，其代数和为零。用方程表示为：迟Q=0(i=1,2,.m-1)2、压降方程(能量方程)气体在管线中流动将产生摩阻损失，长输管线的摩阻损失一般包括两部分，一是气体通过直管段所产生的沿程摩阻；二是气体通过各种阀件、管件所产生的局部摩阻。压降方程反应的就是管段流量与压力损失之间的关系，可用下面的表达式来描述：记二匚宁儿(j=1,2,.P)dj式中和1值与燃气流动状况及管道粗糙度有关，而常数Kj则与燃气性质有关。3、回路方程回路方程表示的是每个基环中各管段的压力损失代数和为零。即?.R=0（n=1,2,.n）上述的三个方程构成的联立方程组亦称为管网的基本方程。现在已有的天然气水力计算方法和计算机软件都是对这三个方程的联立求解。123环状管网的水力计算方法一般地，天然气管网稳态分析是指：己知管径D压力损失厶P、摩阻系数f以及配气站或控制点的压力和节点流量，求出各管段的流量和各节点的压力。当管道中的气体为等温流动时，可以用流量Q压力损失AP、管径D管长L和摩阻系数f来描述其流动状态。当DL、f为已知时，只有Q-P为未知数时，而Q和厶P的关系可以用压降方程来表示。根据求解的未知量不同，可以将管网分析方法分为：（1）通过天然气管网基本方程，将AP消去，以Q为未知量的计算方法，称为流量法，根据方程组的构造方法又可将其分为解环方程法和解管段方程法。3：（2）通过天然气管网基本方程，将Q消去，以P为未知量的计算方法，称为节点压力法。下面将分别介绍这三种方法。1、环方程法环方程法的基本思想是：在满足节点连续性条件下,由能量方程组成非线性方程组，引进环修正流量概念，以此修正流量为未知数，将非线性方程组线性化来求解，方程组的个数为管网基环数。解环方程的主导思想是在满足连续性方程的前提下,逐步修正管段流量减小环闭合差，从而最后满足压降方程。环方程法的典型算法有牛顿一拉夫逊（Newton-RaphSon法和哈迪一克劳斯法等。国内学者李德波（2005年）基于哈迪一克劳斯法原理，通过严密的数学推导,给出了节点流量修正量的一般表达式，从而解决了复杂管网的流量分配问题，是改进目前大城市复杂环状管网水力计算的一种有效方法。2、管段方程法管段方程法的基本思想是：将能量方程线性化处理，与节点连续方程联立形成以管段流量为未知数的线性方程组来求解，方程个数为管网管段数。由于其原理主要是将非线性的能量方程线性化，故也可称之为线性化方法。其具体过程是，先将回路方程线性化，然后结合连续性方程组成线性方程组进行迭代求解。田贯三等学者（2002年）总结了城镇燃气管网稳态分析管段方程法的基本原理及其数学推导过程，并对管网稳态分析几种常用方法进行了系统的研究与比较，在计算工作量、计算精度、收敛速度等方面做出了综合评价，为提高编程质量提供了理论依据5。3、节点方程法节点方程法的基本思想是：将管段流量通过管段压降方程转化为用管段两端的节点压力表示之后，代入节点连续方程，形成以节点压力为未知量的非线性方程组来求解，方程个数为管网的节点数（一般情况下不包括参考节点）。解节点方程法的具体过程是，先将压降方程改变形式，用各管段两端节点的压力差表示管段的压力损失，使其成为关于节点压力的表达式，然后代入连续性方程构造关于节点压力的非线性方程组，通过各种算法来求解这个非线性方程组。应用于节点方程法的典型算法有牛顿拉夫逊法、有限单元法和线性逼近法段常贵等（2000年）在对燃气管网进行水力分析的基础上，建立了通用的管网稳态分析数学模型，用简单迭代法（原理同线性逼近法）和牛顿拉夫逊法分别求解数学模型，并根据两种算法的不同特点将其结合起来。通过两个算例来比较各种算法的计算效果，提出了稳定、高效的求解模型算法。6除此之外，还有很多管网计算方法，如拟牛顿法7、蒙卡特罗法、延拓法等等，然而无论使用那种方法，都是在满足天然气管网连续性方程、压降方程和回路方程基础上求解未知管径、流量和压降。2燃气管网水力计算模拟软件随着计算机技术的发展，现有很多关于燃气管网的软件。燃气管网仿真经历了由数字仿真到可视化仿真阶段，而可视化仿真是数字模拟与科学计算可视化技术相结合的产物。这些软件不仅仅能进行水力计算，还能进行管道模拟仿真，越来越集成化。2.1单一水力计算软件编制西南石油大学袁崇明、汪玉春、李长俊等人在建立数学模型和求解模型所采用的数学方法上作了许多探讨，并且取得了一定的成果，其开发的GASFLOV和EGPN动态仿真软件已得到了初步应用，但这两套软件的建模思想均忽略了温度变化。可见，其计算所得结果对于低压燃气管网是可以接受的，相对误差较小，但对于长输管线、复杂的燃气管网以及在燃气和环境温差较大时计算所得结果将会与实际情况有较大偏差。石油大学吴长春等人与石油规划总院进行了干线输气管道（管网）稳态优化运行方面的研究，并研发了稳态优化运行软件GPNOPT22002年，石油大学又开发了西气东输稳态优化运行软件WEGPOP另外，哈尔滨工业大学、天津大学、同济大学、重庆大学的许多科研人员在燃管网方面也取得了很大的研究成果。由于受到动态仿真研究水平的限制，我国目前在管道在线控制系统仿真方面还没有进行过任何研究工作，这方面将是我国学者今后研究的一个重要方向。等（1977年）介绍了一种天然气管网稳态分析计算程序CASPR（CapabilityAndSimulationProgram），此程序基于“牛顿-拉夫逊节点方程法”的同时增加了优化节点编号、按相同比例分配进气量等功能，成功地应用于工程实际，证明是一种合理有效的天然气管网稳态分析工具。王荣和等（2000年）根据优化原理将管网稳态分析问题转化为单目标优化问题，其目标函数由节点计算气压与SCADA实测气压的平方差、管段流量与SCADA实测流量的平方差和各气源计算供气量与实际供气量的平方差这三部分的和组成，并充分考虑了各种约束条件，用遗传算法对燃气管网进行了分析，从而扩充了燃气管网系统优化设计软件WDOC软件系统的功能，在多个大中型城市的实际燃气系统中加以应用取得了良好的效应。2.2相关仿真软件介绍到20世纪末，输气管道或管网稳态模拟和优化运行技术已基本成熟。近几年来，国外的一些软件公司，如AdvanticaStoner公司和Simone公司、美国TETCO输气管道公司等都进行了输气管网非稳态下优化运行的研究，并且已经在模型和算法上取得了一定的进展。现在己有的多种天然气管网稳态分析的软件、仿真程序，虽各有特色，但究其根本，都是对连续性方程、压降方程和回路方程的联立求解。一个好的燃气管网仿真软件应至少能够实现以下功能：图形的绘制及修改；图形的识别及数据的输入；管网水力计算；管网稳动态仿真；管网负荷动态平衡；管网可靠性分析；管网遥测监控；仿真图表的输出等。(1) TGNE软件美国的SSI公司开发的TGNET(Transientgasnetwork)软件是一套功能较强的燃气长输管道的仿真软件，也是我国最早引进的燃气管网仿真软件，已在世界上45条油气管道上应用。该软件的自动选择时步功能确保了计算结果精确，且其计算速度高于实时运行速度，另外，该软件含有精确的普通管道设备模型，如球阀、泵、止回阀及调节阀。图2TGNET界面该软件可以对管道运行的瞬态水力状况进行模拟，其在线模拟系统由实时模型、预测模型和自动先行模拟等几个模拟软件组成。其中，实时模型是根据管道液（气）体流动的连续性方程、动量方程、能量方程及状态方程建立的数学模型，编制实时模拟应用程序，再根据管内流体的温度、压力、流量、密度及组分等参数完成管道的泄漏检测与定位、超低压保护、批量及组分跟踪、热值跟踪等功能。而预测模型可计算某一时刻管道的存储量、管内批量油品的位置，预测喘振条件,检测输差等。自动先行模拟软件是根据目前的实时模拟条件和设备的设定来确定将来的运行条件有何变化，预测水击的产生，做到超前保护。管道在线仿真系统的应用可有效地提高管道运行的安全性和经济性。（2）SynerGEEGas软件SynerGEEGas是天然气运输和分输网络模拟软件，用于天然气运输和分输管道系统的一个离线的液压模拟工具（网络分析和计划），可为天然气运输分输公司、设施和工程公司提供各种模拟分析，诸如满足政府在管道运行方面的立法要求，更有效、安全的管道运行、提供更多管道目前信息和高级应用。可将安全、高效的网络运行最大化，增加产出、降低运行成本、运行状况分析、潜在问题的警告等，满足法律法规的要求，将管道数据转换成运行信息等。图3SynerGEEGas软件示意图此外，还有美国Stoner公司开发了用于气体稳态管网设计的软件SWS,以及模拟长输管道动态工况的软件SPS。还有Williams管道公司开发的CBT(Computer-BasedTranining)软件，以及法国燃气公司GazdeFrance开发的ACROLOLE的软件。意大利Snamprogetti公司开发了专门用于气体管网的仿真软件DTS(DisPatchinTutoralSystem)，该软件可以通过SCADA系统接口软件，以在线仿真软件方式运行，也可以离线形式运行。1 参考文献刘伟,唐文芳等.城市燃气管网水力计算软件的开发.油气储运.2002,21(1):2 2527段长贵.燃气输配.中国建筑工业出版.2007.(3) 白建辉.天然气管网稳态分析综合方法及比较研究.西南石油大学.2007李德波.复杂环状管网水力计算方法探究.节能技术,2005;23(132):327330田贯三.城镇天然气管网水力分析数学模型与计算方法.天然气工业,2002;22:9698.3 段常贵,徐彦峰等.燃气管网的稳态分析.煤气与热力,2000;20(2):9597.4 谢伟光.牛顿法在燃气管网节点压力计算中的应用.煤气与热力,1998.5 宋汉成,蒋永明等.城市燃气管网仿真研究.