6节点天然气潮流计算程序 （matlab程序）

0.代码链接

6节点天然气潮流计算程序（matlab程序）资源-CSDN文库

1.简述

  天然气潮流计算理论

由于天然气涉及到流体的运动方程，直接计算非常复杂，因此需要提前做出一些假设来简化计算，经过研究，适当的假设对结果影响很小，因此本文对天然气系统做出如下假设： 

1）管道内的天然气为等温气体，气体温度不会随着传输距离发生改变；

2）各个管道对天然气的摩擦系数相同； 

3）各个管道的气体压缩因子相同。

计算步骤如下： 

（1）由用户负荷求得微型燃气轮机的功率𝑃g，求出其消耗的天然气流量； 

（2）根据天然气节点压力与节点间流量的关系式求出前一个节点的节点压力； 

（3）如果节点中存在压缩机，则可以求得压缩机消耗的等效天然气流量，将其作为一个负荷加入到潮流计算中； 

（4）求出节点的节点流量方程； 

（5）重复步骤2)-4），直至求解出整个天然气网络潮流。 

2.部分代码

%% 求解天然气潮流
%主循环
for p=3:1e-5:5.5
    Hk=Bk*fp(1,6)*((fp(2,6)/p)^(Zk*0.3/1.3)-1);
    f6=1.07e-7*Hk; 
    fp(1,5)=f6+fp(1,6);
    fp(2,5)=Pi(fp(1,5),p,Pline(5,4),1);%求解出5节点压力值
    fp(2,3)=Pi(fp(1,5),fp(2,5),Pline(4,4),1);%求解出3节点压力值
    fp(2,4)=Pi(fp(1,4),fp(2,3),Pline(3,4),2);%求解出4节点压力值
    fp(1,3)=fp(1,4)+fp(1,5);%3节点总流量等于4、5节点流量之和
    fp(2,2)=Pi(fp(1,3),fp(2,3),Pline(2,4),1);%求解出2节点压力值
    fp(1,2)=fp(1,3)+15; %2节点流量为3节点流量+2节点负荷
    fp(1,1)=Pi(fp(1,2),fp(2,2),Pline(1,4),1);%求解出1节点压力值
    
    wucha=5-fp(1,1);
    if (wucha<1e-5 )
    break;
    end
end

.。。。。。

3.运行结果

4.参考文献

[1]刘健锋. 电-气区域综合能源系统能流计算及其收敛性调整[D].湘潭大学,2021.DOI:10.27426/d.cnki.gxtdu.2021.001665.

[2]张萌. 电-气综合能源系统多能流计算与最优能流分析[D].山东大学,2022.DOI:10.27272/d.cnki.gshdu.2022.002028.