1. 引言
近年来,天然气作为一种清洁能源走进了千家万户,其在我国能源结构中所占比例不断攀升,但仍远落后于发达国家。未来几年,我国的天然气工业将会有飞跃式的发展。目前,从霍尔果斯入境的中亚天然气经过西气东输二线送达华东和华南;从云南入境的缅甸天然气经中缅天然气管道送达云贵;而华北地区的用气紧张,也将会随着陕京四线的建成而解决。
早在上世纪60年代,美国就开始对输气管道优化运行进行研究,当时仅限于处理稳态问题。随着仿真技术的成熟,国外学者在输气管道(网)稳态运行优化的研究上已经有较为成熟的模型,并将很多模型应用到实际管道运行中。纵观近几年的国外文献,可以发现从事天然气管道运行优化的学者们大多将研究重点放在天然气长输管道非稳态运行优化上。截至到2011年为止,以全线能耗最小为目标函数的天然气长输管道非稳态运行优化模型仍然没有成熟的方案,而且对非稳态优化的求解方法也有待进一步探索[1] [2] [3] [4] [5]。
上世纪90年代我国才开始对输气管道运行优化进行研究,虽然在理论上与国外有一定差距,但是许多国内学者仍然针对我国实际管线建立了数学模型,在模型求解的方法上做了很多的研究和探索,并开发了不少实用性软件。2004年,北京华油天然气公司和中国石油大学(北京)合作开发了陕京线稳态优化运行软件SJGPOPT,计算结果表明陕京线有较大的节能降耗潜力。2008年杨义在其硕士毕业论文中研究了中国主干输气管网的稳态优化运行问题,其将管网拆成单根管道,应用枚举法和顺序动态规划法依次求解每条管道的优化运行方案,最后得到管网的优化运行方案。基于上述理论,其开发了输气管网稳态优化运行软件SimuOpt2007,并用陕京一线和西气东输一线作为算例,对算法和软件进行了验证[6] [7] [8] [9]。
本文针对单条天然气管线开发能耗计算软件,用户可以随意增加站场,并在站场之间加入阀室,输入压气站压缩机基本参数,就可以计算出能耗。同时,该软件还可以优化能耗,继而计算出该管线的能耗最优值。
2. 软件介绍
2.1. 数学模型
2.1.1. 问题分析
在起点进气量给定后,输气成本主要取决于该系统压缩机的总能耗费用,因此选择总能耗费用最小作为该系统稳态运行优化的最优准则。稳态运行优化问题可以分为两个层次,即管线级优化和站内优化[10],管线级优化可确定各个压气站最优进/出站压力,站级优化是在已知压气站进/出站压力和进站流量的前提下确定站内最优开机方案。在优化问题求解过程中,两个优化层次相互影响,管线级优化需要站内优化得到在确定压气站进/出站压力和进站流量下的最优压气站能耗费用;而站级优化只有在管线级优化给定进/出站压力和进站流量的情况下方可进行。
2.1.2. 管线级优化[10]
在稳态运行工况下,当气源供气流量、供气压力、供气温度;各分输点分输量;沿线地温和总传热系数等参数确定后,对任意压气站的一个出站状态(包括压力和温度),都能唯一确定下一站的进站状态。该出站状态还可结合该压气站的进站状态和进站流量确定该站的多种运行方案。下面建立管线级优化顺序形式的动态规划模型[11]。
如图1所示,建立动态规划模型时设压气站数目为n。
Figure 1. Schematic diagram of dynamic programming model
图1. 动态规划模型示意图
将第k − 1个压气站出口到第k个压气站出口的输气过程视为对应问题的第k阶段,第k阶段的状态变量xk (即该阶段的起点状态)为第k − 1站的出站压力(k ≥ 2),决策变量dk为第k站的压比,阶段效应为第k站的能耗费用,则可建立针对该管线各压气站压比的的动态规划模型:
状态变量:
k = 2~n + 1
决策变量:
k = 1~n
状态演变方程:
k = 1~n
阶段效应:
k = 1~n
指标函数:
k = 1~n
最优指标函数:
k = 1~n
函数递推方程:
k = 1~n
初始条件:
x1 = Pd0 (气源的供气压力)
f0(x1) = 0
其中,Pdk为第k站的出站压力;
Dkr(xk+1)为第k阶段相应于其终点状态xk+1的允许决策集合;
P1,k(xk+1):从多阶段过程起点x1到达第k阶段终点状态xk+1的前部子过程的策略,它是由前k个阶段的决策值构成的序列;
Ck(dk, xk+1)为第k个压气站的能耗费用函数。
2.1.3. 站级优化[10]
该系统所有压气站的压缩机均是并联方式,对于并联压缩机的流量分配问题,有流量均分和流量高低分配组合两种方式,据文献[10]介绍,对于同型号机组并联的情况,如果不考虑每套压缩机流量上/下限约束,理论上可以证明流量均分是最优的。
2.2. 压缩机特性工艺计算[11]
陕京二线进行优化计算时需要确定以下数据,单套压缩机的最高与最低压比;最高与最低转速;单套压缩机在入口状态下的可行流量范围;压气站功率;压气站能耗费用。由此可见,该系统优化的核心就是压气站优化,而压气站优化的关键是计算压气站在指定工况下的能耗费用。
(1) 压缩机组能耗费用计算公式
压缩机组的能耗费用计算分两步,先求出压缩机组的输入功率,然后结合驱动机的效率及能耗单价得到能耗费用值。该系统包括燃驱离心压缩机组和电驱离心压缩机组,故只考虑这两种压缩机组能耗费用的计算。
如果驱动机为电动机,则压缩机组的能耗费用为:
(1)
其中,Cost为单套压缩机组的能耗费用,元/天;
Ns为压缩机组输入功率,kW;
ηd为电机效率;
Pd为电价,元/(kW·h)。
(2) 压缩机特性曲线
在回归和应用压缩机特性方程时,需要注意两个问题:
① 在压缩机的各种特性方程中,压缩机的流量通常是指在其入口状态(包括温度和压力)下的体积流量。如果特性曲线图上的流量是标准状态(101,325 Pa, 20℃)下的体积流量,则在回归特性方程时应首先应将标准体积流量换算为实际入口状态下的体积流量。
② 特性方程是在压缩机特性曲线的测试条件下回归得到的,当利用这些特性方程计算压缩机实际工艺参数时,需要根据相似原理将测定条件下的特性方程换算成实际条件下的特性方程。
研究压缩机特性,先根据需要将标准状况下的流量换算成压缩机进口状态(进站压力与进站温度)下的体积流量,换算公式为:
(2)
其中,Q1、T1、P1和Z1分别为压缩机进口状态下的流量、温度、压力及压缩因子,单位分别为m3/d、K、MPa;
Q0为标况下的流量,单位为Nm3/d。
然后将Q1化为在转速n0下的对应流量:
(3)
对单套压缩机组,需要先求得其输入功率,如公式(4)所示。
(4)
其中,H为多变能头,kJ/kg;
G为质量流量,kg/s;
Η为压缩机的效率。
从上式可知,需要回归压缩机的压头–流量曲线和效率–流量曲线。
① 压头–流量特性方程
如果所给特性曲线图上的流量是标准状况下的值,需要先利用公式(2)将标准流量换算成压缩机进口状态下的流量。然后选择对应转速n0的压头–流量曲线进行取点回归,得到压头–流量特性方程为:
(5)
其中,b0、b1和b2为系数;
H0为压缩机进口状态下的体积流量Q1对应的压头,m。
设压缩机的实际转速为n,标准状况下的体积流量为Q0。为了利用公式(5)计算压缩机的压头,首先需要用公式(2)将Q0换算成压缩机进口状态下的体积流量Q1,然后利用(3)将Q1换算成转速n0下的相应流量Q,并将得到的Q代入公式(5)得到H0。最后将H0代入式(6)得到压缩机的实际压头:
(6)
② 效率–流量特性方程
如果效率–流量特性曲线图上的流量是标准状况下的值,同样需要利用公式(2)将标准流量换算成压缩机进口状态下的流量。然后选择某一转速下的效率–流量曲线进行取点回归,得到对应转速n0的效率–流量特性方程为:
(7)
其中,a0、a1和a2为系数;
η0为压缩机进口状态下体积流量Q1对应的效率。
设压缩机的实际转速为n,标准状况下的体积流量为Q0。为了利用公式(7)计算压缩机的效率,首先需要用公式(2)将Q0换算成压缩机进口条件下的体积流量,然后利用(3)将Q1换算成在转速n0下的相应流量Q,将得到的Q代入效率–流量特性方程(7),得到效率η。
3. 举例
本文以陕京二线为例,如表1所示,各压气站内压缩机在其额定转速下回归得到的曲线方程。
Table 1. Characteristic equation of Shaanjing second line compressor
表1. 陕京二线压缩机特性方程
|
榆林 |
临县 |
压头–流量 |
H = 64.9 + 9.1613Q − 0.2132Q2 |
H = 60.058 + 1.8436Q − 0.0467Q2 |
效率–流量 |
Η = −0.0836 + 0.0645Q − 0.0011Q2 |
η = −0.2539 + 0.0329Q − 0.0045Q2 |
滞止–转速 |
Qmax = 1.053 + 0.029n + 10 − 7n2 |
Qmax = 1.492 + 0.0044n + 10−7n2 |
喘振–转速 |
Qmin = 9.4383 − 0.008n + 2×10 − 7n2 |
Qmin = 8.717 − 0.0006n + 2 × 10−7n2 |
转速–滞止 |
n = 901.58 + 263.72Q − 1.0262Q2 |
n = 60.888 + 204.6Q − 0.4812Q2 |
转速–喘振 |
n = −2376.6 + 887.48Q − 14.062Q2 |
n = −1151.5 + 630.73Q − 8.158Q2 |
|
阳曲 |
石家庄 |
压头–流量 |
H = 57.274 + 2Q − 0.0488Q2 |
H = 59.072 + 1.9086Q − 0.0476Q2 |
效率–流量 |
η = −0.2037 + 0.03923Q − 0.0005Q2 |
η = − 0.2593 + 0.0328Q − 0.00045Q2 |
滞止–转速 |
Qmax = 3.653 + 0.004n + 2 × 10−7n2 |
Qmax = 0.4066 + 0.0051n + 10−7n2 |
喘振–转速 |
Qmin = 9.3548 − 0.0009n + 3 × 10−7n2 |
Qmin = 8.6737 − 0.0007n + 3 × 10−7n2 |
转速–滞止 |
n = −259.55 + 204.57Q − 0.6457Q2 |
n = 86.806 + 183.95Q − 0.3721Q2 |
转速–喘振 |
n = −1040.1 + 582.09Q − 7.5623Q2 |
n = − 1095.5 + 590.37Q − 7.6374Q2 |
3.1. 相关工艺参数
(1) 标准状态为:压力1.01325 × 105 Pa (绝对压力),温度20℃。
(2) 年工作天数为350天。
(3) 管道的管内壁绝对当量粗糙度取10 μm。
(4) 各压气站内压降取0.2 MPa,其中压气站入口压降取0.15 MPa (主要考虑分离设备、计量设备及管路的摩阻损失),出口压降取0.05 MPa (主要考虑管路空冷器的摩阻损失)。各分输站站内压降取0.05 MPa。
(5) 在本案例分析中,沿线温度取各地的年平均地温;沿线总传热系数均取1.75 W/(m2∙K)。
(6) 在本论文案例分析中,各站最高出站压力为10 MPa (绝对压力);榆林压气站出站温度上限取60℃,其余各压气站出站温度上限取55℃。
(7) 长庆来气压力为5 MPa,进站温度10℃;通州最低允许进站压力4.2 MPa。
(8) 各管段的输气效率统一取0.94。
(9) 电价每度1.00元,电机效率95%。
(10) 气量在安平站分输一半给联络线。
3.2. 优化案例分析
3.2.1. 固定进气量(3700万方/天)
Table 2. Operation optimization plan (37 million cubic meters/day)
表2. 运行优化方案(3700万方/天)
|
进站流量 |
进站压力 |
出站压力 |
压比 |
开机 |
转速 |
进站温度 |
出站温度 |
能耗 |
单位 |
万方/天 |
MPa |
MPa |
|
套 |
rmp |
℃ |
℃ |
万元/天 |
榆林 |
3700 |
5 |
9.9 |
1.96 |
3 |
8110 |
10 |
60 |
85.16 |
兴县 |
3700 |
8.28 |
9.8 |
1.187 |
2 |
5544 |
27.4 |
41.8 |
20.46 |
阳曲 |
3700 |
7.67 |
8.5 |
1.109 |
2 |
4669 |
19.8 |
29.3 |
13.02 |
石家庄 |
3700 |
6.26 |
6.26 |
1.000 |
0 |
0 |
17.4 |
17.4 |
0.00 |
通州站 |
1850 |
4.32 |
|
|
|
|
|
|
|
总计 |
|
|
|
|
7 |
|
|
|
118.63 |
由上表2可以得出:该软件在设定进气量的情况下,可以较快计算出全线的优化运行工况。
3.2.2. 改变进气量
设进气量由3100万方/天增加到4300万方/天,步长定为200万方/天。
Table 3. Changes in intake volume (increased from 31 million cubic meters per day to 43 million cubic meters per day)
表3. 改变进气量(由3100万方/天增加到4300万方/天)
进气量 |
3100 |
3300 |
3500 |
3700 |
3900 |
4100 |
4300 |
能耗 |
77.36 |
86.43 |
98.98 |
118.63 |
136.04 |
156.04 |
180.54 |
将表3中的进气量为横坐标,以能耗为纵坐标,得到下图2
Figure 2. Energy consumption change chart (increased from 31 million cubic meters per day to 43 million cubic meters per day)
图2. 能耗变化图(由3100万方/天增加到4300万方/天)
将由上图可以得出:1) 随着输量的增加,全线运行能耗费用增加;2) 与实际运行费用相比,随着数量增加,节能效果越来越不明显;3) 增输相同的气量,所需运行费用随起点输量的增加而增加。
3.2.3. 管线效益分析
假设该管线的管输费用单价为0.1元/方、0.15元/方和0.2元/方。将其结合3.2.2中的表可得如下三表(表4~6)数据。
Table 4. Energy consumption and profit comparison table (0.1 unit price)
表4. 能耗和利润对比表(0.1单价)
进气量 |
3100 |
3300 |
3500 |
3700 |
3900 |
4100 |
4300 |
能耗 |
77.36 |
86.43 |
98.98 |
118.63 |
136.04 |
156.04 |
180.54 |
利润 |
232.64 |
243.57 |
251.02 |
251.37 |
253.96 |
253.96 |
249.46 |
Table 5. Energy consumption and profit comparison table (0.15 unit price)
表5. 能耗和利润对比表(0.15单价)
进气量 |
3100 |
3300 |
3500 |
3700 |
3900 |
4100 |
4300 |
能耗 |
77.36 |
86.43 |
98.98 |
118.63 |
136.04 |
156.04 |
180.54 |
利润 |
387.64 |
408.57 |
426.02 |
436.37 |
448.96 |
458.96 |
464.46 |
Table 6. Energy consumption and profit comparison table (0.2 unit price)
表6. 能耗和利润对比表(0.2单价)
进气量 |
3100 |
3300 |
3500 |
3700 |
3900 |
4100 |
4300 |
能耗 |
77.36 |
86.43 |
98.98 |
118.63 |
136.04 |
156.04 |
180.54 |
利润 |
542.64 |
573.57 |
601.02 |
621.37 |
643.96 |
663.96 |
679.46 |
由上述三表可以得出:只要该管线的管输费用单价不小于0.1元/方,则输量大利润大;且随着管输费用单价的增加,输量越大利润越大。
4. 结论
(1) 该软件在设定进气量的情况下,可以较快计算出全线的优化运行工况;
(2) 1) 随着输量的增加,全线运行能耗费用增加;2) 与实际运行费用相比,随着数量增加,节能效果越来越不明显;3) 增输相同的气量,所需运行费用随起点输量的增加而增加;
(3) 只要该管线的管输费用单价不小于0.1元/方,则输量大利润大;且随着管输费用单价的增加,输量越大利润越大。