《第十章天然气管网稳定流水力计算》由会员分享，可在线阅读，更多相关《第十章天然气管网稳定流水力计算（16页珍藏版）》请在金锄头文库上搜索。
1、-10 稳定流燃气管网计算方法与模型10.1燃气管网水力计算数学模型与方法10.1.1 燃气管网水力计算的数学模型用计算机进展燃气管网水力计算，首先需要把管网的信息输入到计算机中去，这就必须用数学的语言描述管网的构造，这一任务可借助图论来完成，图10-1为一简单的管网示意图。1234567812610934587Q6Q8Q7Q4Q1Q2Q3Q5q1q2q6q10q9q3q4q5q7q8IIIIII图10-1 管网示意图图1管网示意图图中 1，2，10接点编号； 1，2，10管段编号；，环编号； Q1，Q2，Q10节点流量； q1，q2，q10管段流量。由图论可知，任何环状管网在管段为p，节点数
2、数为m，环数为n的情况下，其管段数、节点数和环数存在以下关系： p=m+n-1燃气管网供气时，在任何情况下均需满足管道压降计算公式，节点流量方程和环能量方程，其中后两个方程称为根本方程。10.1.1.1管段压力降计算公式 j=1,2,p (101)式中 管段的阻力系数；管段的压力降；管段j的流量；常数。 可列出p个管段压降计算公式。10.1.1.2节点流量连续方程对燃气管网任一节点i均满足流量平衡，可用下式表示： i=1,2,m 10-2式中 管段j与节点i的关联元素，管段j与节点i关联，且是管段的起点，管段j与节点i关联，且是管段的终点，管段j与节点i不关联。 可建立个独立的方程10.1.1
3、.3环能量方程对于燃气管网中任一环路均应满足压降之和为零，可用下式表示： i=1,2,n (103)式中 bij管网环路与管段的关联元素，bij=1管段j在第i个环中，且管段j的方向与环的方向一致，bij= -1，管段j在第i个环中，且管段的方向与环的方向相反，bij=0，管段j不在第i个环中。 可建立n个独立的环能量方程。10.1.2三种计算方法总之，对于一个管网，当管径时，每条管段有压降和流量两个未知数，共有2p个未知数，可列出的方程数为： 10-4这样未知数与方程的个数相等，可以进展求解，方程组为非线性的，直接求解困难，一般可通过以下三种方法求解。10.1.2.1 解环方程法 在满足连续
4、方程组10-2的条件下，用求解各环校正流量的方法，来间接解出各管段流量的方法叫解环方程法，也就是Hardy Cross法。对第i环列出能量方程，最初确定的管段设计流量一般不能满足能量方程，其能量可用下式表示： i=1，2，n 105式中： pi第i环的压降不闭合差。为了使各环的压降闭合差到达允许的计算精度，保证节点流量平衡，引入环校正流量来消除各环的闭合差，对每环的闭合差引入校正流量qii=1，2，n，则第i环的能量方程可以改写为： i=1，2，n 10-6式中 第i环校正流量，其正负号与bij一致；第i环第j管段邻环校正流量，其正负号与bkj相反。将式10-6括号内多项式展开为麦可劳林级数，
5、因为qi，qk与qj相比甚小，故只取展开式的前两项。 10-7将式10-7代入式10-6得： i=1，2，n10-8将所有的各环能量方程联立，则形成求解i=1，2，n的线性方程组，其方程组的矩阵表示形式为：10-9式中 B由元素bij组成的环路关联矩阵；BT矩阵的转置矩阵；q由i=1，2，n组成的向量； q由管段流量qjj=1，2，p组成的向量； R由组成的对角矩阵。计算步骤：首先确定出各管段的初始计算流量q0，形成线性方程组10-9，求解得出各环校正流量q1，对q0进展校正得q1，判断q1是否满足计算精度要求，未满足要求再重新形成线性方程组10-9，求解校正流量q2，对q1进展校正，其管网流
6、量校正通式为： 10-10进展循环校正，直到第l次校正后的j=1，2，p，满足精度要求为止，在计算过程中，把控制每环的压降残差来到达计算精度要求，转化为控制各管段前后两次修正后的管段流量差满足一定的计算要求为止，迭代完毕后，根据管段压降计算式算出各管段压降，从给定压力的基准点推算出各节点压力等参数。不考虑邻环的影响时，称为单一回路法；考虑邻环影响时，称为联立回路法。10.1.2.2解节点方程法 以节点连续方程为根底，把方程中的管段流量通过管段压降计算公式，转化为用管段两端的节点压力表示，这样连续方程转化为满足能量方程，以节点压力为变量的方程组，通过求解方程组便可得各节点压力，此法称为节点法。
7、节点发按其解法分为有限元节点法和联立节点法。 A有限元节点法 对燃气管网进展水力计算，要求满足以下三个方程组 a节点流量连续方程组Aq+Q=0 b管段压力降方程组ATP=p c管段流量方程组q=Cp由上述三式可得，求解节点压力的方程组。 ACATP+Q=0 10-11式中 A由元素aij组成的节点关联矩阵； C由元素组成的节点对角矩阵； P节点压力向量； Q节点流量向量； q管段流量向量；p管段压降向量； AT矩阵的转置矩阵。 计算步骤：首先初设管段流量q0，形成方程组10-11，求解节点压力p1，计算出q1；q1不满足要求进展修正，再形成方程组1011进展逐次逼近，直到第K+1次的qK+1与
8、qK差的绝对值满足计算精度要求为止。 B联立节点法联立节点法也称为牛顿拉普森法，求解节点方程的数学模型为： i=1，2，m 10-12将上式按台劳级数展开，为了简化计算，取一次项来逼近。 i=1，2，m 10-13这就是联立节点法所求解的方程组，实际计算中应写成下述迭代形式。 i=1，2，m 10-14其方程组的矩阵表示形式为： 10-15式中 pi，p1第i管段的起点压力和终点压力；pi节点i的压力修正值； Rk由形成的对角矩阵；第k+1次求解方程组的解向量。求解出节点压力修正值后，按下式进展修正。 10-16 按给出的迭代形式进展迭代求解，最后解得满足计算精度要求即可。10.1.2.3解管
9、段方程法将节点连续方程和环能量方程联立形成有p个独立方程的方程组，其个数为管网管段数，将其转化为以管段流量为变量的方程组。由于能量方程为非线性方程，难以直接求解，因此可以通过线性化进展迭代逼近，其数学过程表示如下： (10-17)设则3-17式可表示为: (10-18) 式中 C由元素cij组成的矩阵;当im-1时,im+1时,l=i-m+1;j=1,2,p。计算步骤：首先设各管段初始流量，形成线性化方程组10-18，求解q1，当q1不满足计算精度要求时，对管段流量按下式修正后进展迭代求解。 10-19式中 形成线性化能量方程系数的管段计算流量；流量修正系数，取00.5。 当迭代到满足精度要求
10、，计算出管网其它参数，输出计算结果即可。10.2三种算法的比较与评价10.2.1 方程组矩阵的性质三种算法方程组的系数矩阵均为稀疏矩阵，其中解环方程法和解节点方程法的系数矩阵为正定对称矩阵，解管段方程法的系数矩阵为一般大型稀疏矩阵，三种算法的系数矩阵均可压缩一维贮存，节省计算机内存。解环方程法和解节点方程法的方程组易形成，编程简单，解管段方程法的方程组形成较复杂，编程难度大，且占内存多。10.2.2计算工作量三种算法的方程个数分别为管网的环数、节点数减1和管段数，所以三种算法的工作量依次为：解环方程法最小，解节点方程法居中，解管段方程法最大。10.2.3 对计算初值的要求解环方程法需设管段流量
11、的初值，要求管段流量初值必须满足节点连续方程。有限元节点法是求解各节点压力，需初设各管段流量，对管段流量初值要求不高；联立节点法是求解各节点压力修正值，需初设各节点压力，对各节点压力初值要求较高。解管段方程法需初设各管段流量初值，对管段流量初值要求不高。10.2.4收敛速度与计算精度解环方程法中的联立回路法收敛性好，计算精度较高，单一回路法收敛速度慢，计算精度低。有限元节点法和联立节点法的收敛性和计算精度都比较好，单在使用联立节点法时如节点压力初值选取不当，则有限元节点法在计算精度和收敛速度上优于联立节点法。节点法在平差时，如遇到大管径低摩阻的管段，收敛速度和计算精度都降低。这是因为各管段的流
12、量是由两端的压差求得，当大管径管段压差很小,对节点压力影响的灵敏度与相邻管段不同时，所求得管段流量难以满足计算精度要求。解管段方程法收敛速度快,计算精度高，可以用它的计算结果作为标准，来衡量其它两种算法的精度。10.2.5原始数据准备工作量三种方法均需输入节点数、管段数、每条管段的起点号、终点号、管长、管径。解环方程法和解管段方程法需要环的信息，输入环与管段的关联矩阵。因此节点法的原始数据准备工作两最少。综上所述，在一般情况下进展燃气管网平差应优先选用节点法或解环方程法，在计算精度要求高时，应选用解管段方程法。10.3大管径低摩阻不收敛问题的研究在燃气管网水力计算过程中，当燃气管网中的*管段的
13、管径过大或流量过小时，存在不收敛或收敛精度低，影响水力计算的准确性。国内外有关文献均报导过，但未找出原因。本工程通过深入多年的研究，找出了原因，并提出了解决的方法。 不收敛的原因在燃气管网水力计算过程中，国标GB50028-93规定根据燃气在管道中不同的运动状态，其摩擦阻力系数一般按以下格式计算：10.3.1.1 低压燃气管道1层流状态 2临界状态 错误！无效。3湍流状态 钢管 铸铁管 式中 雷诺数；管道内外表的当量绝对粗糙度，对于钢管取0.2mm；燃气管道的摩擦阻力系数；运动粘度m2/s；管道内径mm；燃气管道的计算流量m3/h。10.3.1.2高中压燃气管道国标GB50028-93规定高中压管道的流态按湍流考虑。钢管 铸铁管 而根据燃气设
