第十章天然气管网稳定流水力计算分解.docVIP

10 稳定流燃气管网计算方法与模型 10.1燃气管网水力计算数学模型与方法 10.1.1 燃气管网水力计算的数学模型 用计算机进行燃气管网水力计算，首先需要把管网的信息输入到计算机中去，这就必须用数学的语言描述管网的结构，这一任务可借助图论来完成，图10-1为一简单的管网示意图。 图10-1 管网示意图 图中 1，2，…10——接点编号； （1），（2），…（10）——管段编号； Ⅰ，Ⅱ，Ⅲ——环编号； Q1，Q2，…Q10——节点流量； q1，q2，…q10——管段流量。 由图论可知，任何环状管网在管段为p，节点数数为m，环数为n的情况下，其管段数、节点数和环数存在下列关系： p=m+n-1 燃气管网供气时，在任何情况下均需满足管道压降计算公式，节点流量方程和环能量方程，其中后两个方程称为基本方程。 10.1.1.1管段压力降计算公式 j=1,2,…p (10－1) 式中 ——管段的阻力系数； ——管段的压力降； ——管段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环能量方程 对于燃气管网中任一环路均应满足压降之和为零，可用下式表示： i=1,2,…n (10－3) 式中 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 解环方程法 在满足连续方程组（10-2）的条件下，用求解各环校正流量的方法，来间接解出各管段流量的方法叫解环方程法，也就是Hardy Cross法。 对第i环列出能量方程，最初确定的管段设计流量一般不能满足能量方程，其能量可用下式表示： i=1，2，…n （10－5） 式中： Δpi——第i环的压降不闭合差。 为了使各环的压降闭合差达到允许的计算精度，保证节点流量平衡，引入环校正流量来消除各环的闭合差，对每环的闭合差引入校正流量Δqi（i=1，2，…n），则第i环的能量方程可以改写为： i=1，2，…n （10-6） 式中 ——第i环校正流量，其正负号与bij一致； ——第i环第j管段邻环校正流量，其正负号与bkj相反。 将式（10-6）括号内多项式展开为麦可劳林级数，因为Δqi，Δqk与qj相比甚小，故只取展开式的前两项。 （10-7） 将式（10-7）代入式（10-6）得： i=1，2，…n （10-8） 将所有的各环能量方程联立，则形成求解（i=1，2，…n）的线性方程组，其方程组的矩阵表示形式为： （10-9） 式中 B——由元素bij组成的环路关联矩阵； BT——矩阵的转置矩阵； Δq——由（i=1，2，…n）组成的向量； q——由管段流量qj（j=1，2，…p）组成的向量； R——由组成的对角矩阵。 计算步骤：首先确定出各管段的初始计算流量q（0），形成线性方程组（10-9），求解得出各环校正流量Δq（1），对q（0）进行校正得q（1），判断q（1）是否满足计算精度要求，未满足要求再重新形成线性方程组（10-9），求解校正流量Δq（2），对q（1）进行校正，其管网流量校正通式为：