一种压缩机组智能增减控制方法与流程

1.本发明属于天然气管道压缩机组自动控制技术领域，应用在具有多台压缩机组联合运行的压气站，具体为一种压缩机组智能增减控制方法。


背景技术：

2.长距离输送天然气管道压气站场一般设置有多台压缩机组，用于对管道天然气进行增压，而压缩机组的运行数量是根据目标控制参数(例如：该压气站进、出口压力需求)，并结合当前管网的运行工况进行确定的。
3.目前，压缩机组的控制方式分为手动模式和自动模式两种情况：手动模式为：由操作人员根据当前实时工况，手动下发转速控制命令，人为调节转速；自动模式为：通过设定压缩机组进、出口压力目标控制参数，压缩机组自动调节转速。
4.无论手动模式还是自动模式，随着工况的动态变化，当机组转速接近其上限或下限，但仍未达到目标控制需求时，此时由操作人员判断是否需要增减机组以及增减哪一台机组，并远程下发启、停机命令，并增启或减停一台压缩机，以此确保压缩机组的运行数量能够满足工况需求。
5.然而，目前的这种控制方式，仍存在如下方面问题：
6.(1)压缩机组运行台数的控制需要专人值守，实时监视，依靠人为经验判断，并人工手动操作，人力成本高，且工作效率低下。
7.(2)压缩机组的数量控制不够精确，造成能耗增加、能源浪费，进而效益相应减少；同时，压缩机组工作时间的增加，不仅寿命缩短，而且维护成本增加。


技术实现要素：

8.为旨在至少在一定程度上解决精确计算压缩机组运行数量、并智能选择压缩机运行组合方式、提升机组综合运行效率以及降低运行能耗的技术问题，本发明提出了一种压缩机组智能增减控制方法。
9.为达上述目的，本发明提供了一种压缩机组智能增减控制方法，应用在具有多台压缩机组联合运行的天然气管道压气站，包括如下步骤：
10.s1：通过压缩机控制系统判断压缩机启动或停止时间，若启停机时间大于30分钟，系统自动调节平稳，并进入步骤s2，若启停机时间小于等于30分钟，系统跳出机组运行数量计算模型；
11.s2：获取天然气管道压气站与压缩机组运行相关的实时工况参数，依次建立压缩机组的运行数量计算模型、运行负荷增减机组计算模型和喘振裕度计算模型；
12.s3：根据运行数量计算模型、运行负荷增减机组计算模型和喘振裕度计算模型的计算结果，延时180秒，触发增启或减停机组命令；
13.s4：当需要增启、减停压缩机组时，根据优先级比选结果，自动触发对应机组的启停机命令。
14.优选的，在步骤s2中，所述建立压缩机组的运行数量计算模型具体为：根据压缩机对管道内天然气做功的计算模型，建立最少可启动机组的数量计算模型n，所述压缩机运行数量计算模型n的表达式如下：
[0015][0016]
其中，n代表理论计算得出的压缩机运行数量，n代表实际压缩机运行数量，p
i
代表各压缩机对管道内天然气所做的功，p
e
代表各台压缩机额定功率，η代表综合效率修正系数，一般取常数0.7～0.75。
[0017]
优选的，所述根据压缩机对管道内天然气做功的计算模型的具体表达式为：
[0018]
p＝qh
a
[0019]
其中，p代表压缩机组对管道内天然气所做的功，h
a
为绝热压头，q为天然气质量流量；
[0020]
所述天然气质量流量q的计算表达式如下：
[0021]
q＝d
r
dq
b
[0022]
其中，d
r
代表天然气相对密度，d代表空气密度(单位：kg/m3)，近似取常数1.205kg/m3，q
b
代表压缩机进口的标况体积流量，可由压缩机进口流量计测量给出；
[0023]
所述绝热压头h
a
的计算表达式如下：
[0024][0025]
其中，p
s
、p
d
、t
s
分别代表压缩机天然气进口压力、出口压力、进口温度，m
w
代表天然气摩尔质量，r代表普适气体常数，r＝8.3145j/(mol
·
k)，z
s
代表压缩机进口条件天然气压缩系数，γ代表天然气绝热指数；
[0026]
所述天然气绝热指数γ的计算表达式如下：
[0027][0028][0029]
其中，m代表天然气多变指数，p
s
、p
d
、t
s
、t
d
分别代表压缩机天然气进口压力、出口压力、进口温度、出口温度，z
s
、z
d
分别代表压缩机进口、出口天然气压缩系数；
[0030]
所述压缩机进口、出口天然气压缩系数z
s
、z
d
的计算表达式分别如下：
[0031][0032]
[0033]
其中，p
pc
、t
pc
分别代表临界压力和临界温度，p
pc
、t
pc
计算公式分别通过如下公式求出：
[0034][0035]
优选的，所述触发增启或减停机组命令具体包括：
[0036]
若n＜n
‑
1，延时180秒，触发减停机组命令；
[0037]
若n＞n，延时180秒，触发增启机组命令；
[0038]
若n
‑
1＜n＜n，说明压缩机组运行数量满足工况要求，无需增启和减停机组。
[0039]
优选的，所述运行负荷增减机组计算模型通过负荷占用比load表达，所述负荷占用比load通过如下表达式求得：
[0040][0041]
其中，load代表负荷占用比，r代表平均负荷，r
max
代表最大负荷，r
min
代表最小负荷；
[0042]
所述平均负荷r的计算通过如下表达式求得：
[0043][0044]
其中，n代表实际压缩机运行数量。
[0045]
优选的，所述触发增启或减停机组命令具体还包括：
[0046]
若load＞95％，延时180秒，触发增启机组命令；
[0047]
若load＜5％，延时180秒，触发减停机组命令。
[0048]
优选的，所述喘振裕度计算模型具体通过机组平均裕度m表达，所述机组平均裕度m通过如下表达式求得：
[0049][0050]
其中，m代表机组平均裕度，n代表实际压缩机运行数量，m
i
代表各压缩机喘振裕度。
[0051]
优选的，所述触发增启或减停机组命令具体还包括：
[0052]
若m＜5％，延时180秒，触发减停机组命令。
[0053]
优选的，所述当需要增启、减停压缩机组时，根据优先级比选结果，自动触发对应机组的启停机命令，具体包括：
[0054]
(1)采用数字表示优先级，“1”代表优先级最高，数字越大，优先级越低，依次类推，0表示不可用状态；
[0055]
(2)优先级判定采用手动、自动两种方式，其中：
[0056]
手动方式下，可人为写入优先级数字；
[0057]
自动方式下，比较周期内压缩机组累计运行时间，运行时间越长，优先级越低，反之优先级越高；
[0058]
(3)当触发增启压缩机组命令后，在停机状态下的机组中进行优先级比选，优先级越高，即数字越大，优先启动；
[0059]
(4)当触发减停压缩机组命令后，在运行状态下的机组中进行优先级比选，优先级越低，即数字越小，优先停止。
[0060]
优选的，所述压缩机控制系统包括压缩机单机控制系统、压缩机组负荷分配控制系统和站场scada控制系统，所述压缩机单机控制系统、所述压缩机组负荷分配控制系统和所述站场scada控制系统之间能够实现数据互传，获取压气站与压缩机运行相关的实时工况参数。
[0061]
相对于现有技术，上述压缩机组智能增减控制方法的有益效果为：
[0062]
1、根据压缩机组的运行数量计算模型，给出实时工况下的最优压缩机应启动数量，有效提升机组综合运行效率，降低生产能耗，节约运行成本。
[0063]
2、根据已计算出最优压缩机运行数量，智能控制机组增启减停过程，无需人员参与操作，提高生产效率，释放劳动力。
附图说明
[0064]
图1为本发明实施例中压缩机组智能增减控制方法的整体控制逻辑图；
[0065]
图2为本发明实施例中压缩机组智能增减控制方法的基本流程图。
具体实施方式
[0066]
下面将参照附图详细描述本发明的实施例，描述涉及附图时，除非另有表示，不同附图中的相同附图标记表示相同或相似的要素。要说明的是，以下示例性实施例中所描述的实施方式并不代表本发明的所有实施方式。它们仅是与如权利要求书中所详述的、本发明公开的一些方面相一致的装置和方法的例子，本发明的范围并不局限于此。在不矛盾的前提下，本发明各个实施例中的特征可以相互组合。
[0067]
参见图1
‑
2，本发明实施例提供一种压缩机组智能增减控制方法，包括如下步骤：
[0068]
s1：通过压缩机控制系统判断压缩机启动或停止时间，若启停机时间大于30分钟，系统自动调节平稳，并进入步骤s2，若启停机时间小于等于30分钟，系统跳出机组运行数量计算模型；
[0069]
s2：获取天然气管道压气站与压缩机组运行相关的实时工况参数，依次建立压缩机组的运行数量计算模型、运行负荷增减机组计算模型和喘振裕度计算模型；
[0070]
s3：根据运行数量计算模型、运行负荷增减机组计算模型和喘振裕度计算模型的计算结果，延时180秒，触发增启或减停机组命令；
[0071]
s4：当需要增启、减停压缩机组时，根据优先级比选结果，自动触发对应机组的启停机命令。
[0072]
具体的，在本发明实施例当中，某压气站设有4台压缩机，在运行压缩机2台，额定功率均为20mw，相关参数如下表1：
[0073]
表1实时天然气组分
[0074]
序号组分名称分子式实时值1甲烷ch494.0710％2乙烷c2h62.9710％3丙烷c3h80.5110％
4异丁烷c4h10_10.0780％5正丁烷c4h10_20.0710％6异戊烷c5h12_10.0180％7正戊烷c5h12_20.0210％9己烷c6h140.0670％10氮气n21.4970％11二氧化碳co20.6940％12相对密度\0.584
[0075]
其中，压缩机控制系统给出如下各台运行机组实时工况参数如下表2：
[0076]
表2机组实时工况参数
[0077][0078]
进一步的，在步骤s2中，所述建立压缩机组的运行数量计算模型具体为：根据压缩机对管道内天然气做功的计算模型，建立最少可启动机组的数量计算模型n，所述压缩机运行数量计算模型n的表达式如下：
[0079][0080]
其中，n代表理论计算得出的压缩机运行数量，n代表实际压缩机运行数量，p
i
代表各压缩机对管道内天然气所做的功，p
e
代表各台压缩机额定功率，η代表综合效率修正系数，一般取常数0.7～0.75。
[0081]
由此，在本发明实施例中，由可启动机组数量计算模型，取η＝0.70，得出n＝0.91。
[0082]
进一步的，所述根据压缩机对管道内天然气做功的计算模型的具体表达式为：
[0083]
p＝qh
a
[0084]
其中，p代表压缩机组对管道内天然气所做的功，h
a
为绝热压头，q为天然气质量流量；
[0085]
由此，在本发明实施例中，由压缩机对管道内天然气做功的计算模型：得到1#压缩机做功p1＝6.27mw，2#压缩机做功p1＝6.53mw。
[0086]
所述天然气质量流量q的计算表达式如下：
[0087]
q＝d
r
dq
b
[0088]
其中，d
r
代表天然气相对密度，d代表空气密度(单位：kg/m3)，近似取常数1.205kg/m3，q
b
代表压缩机进口的标况体积流量，可由压缩机进口流量计测量给出；
[0089]
由此，由天然气质量流量计算公式：得到1#压缩机瞬时质量流量q1＝655983.1538kg/h，2#压缩机瞬时质量流量q2＝571820.353kg/h。
[0090]
所述绝热压头h
a
的计算表达式如下：
[0091][0092]
其中，p
s
、p
d
、t
s
分别代表压缩机天然气进口压力、出口压力、进口温度，m
w
代表天然气摩尔质量，r代表普适气体常数，r＝8.3145j/(mol
·
k)，z
s
代表压缩机进口条件天然气压缩系数，γ代表天然气绝热指数；
[0093]
具体的，在本实施例当中，通过站场scada系统给出天然气组分参数，可计算得出当前天然气摩尔质量m
w
为17.0628；
[0094]
由绝热压头计算模型：得到1#压缩机绝热压头h
a1
＝32.4259kj/kg，2#压缩机绝热压头h
a2
＝41.1162kj/kg。
[0095]
天然气绝热指数γ的计算表达式如下：
[0096][0097][0098]
其中，m代表天然气多变指数，p
s
、p
d
、t
s
、t
d
分别代表压缩机天然气进口压力、出口压力、进口温度、出口温度，z
s
、z
d
分别代表压缩机进口、出口天然气压缩系数；
[0099]
在本发明实施例当中，由天然气绝热指数计算公式得到1#压缩机绝热指数γ1＝0.2796，2#压缩机绝热指数γ2＝0.2584。
[0100]
压缩机进口、出口天然气压缩系数z
s
、z
d
的计算表达式分别如下：
[0101][0102][0103]
其中，p
pc
、t
pc
分别代表临界压力和临界温度，p
pc
、t
pc
计算公式分别通过如下公式求出：
[0104][0105]
在本发明实施例当中，通过负荷分配控制系统，根据已建立模型，计算压缩机组运行数量过程为：
[0106]
(1)根据压缩机组实际做功情况，建立最少可启动机组数量计算模型，步骤如下：
[0107]
步骤1：由天然气压缩系数计算公式
[0108][0109][0110][0111]
得到1#压缩机进口压缩系数z
s1
＝0.9005，1#压缩机出口压缩系数z
d1
＝0.9052，2#压缩机进口压缩系数z
s2
＝0.8853，2#压缩机出口压缩系数z
d2
＝0.8886。
[0112]
进一步的，所述触发增启或减停机组命令具体包括：
[0113]
若n＜n
‑
1，延时180秒，触发减停机组命令；
[0114]
若n＞n，延时180秒，触发增启机组命令；
[0115]
若n
‑
1＜n＜n，说明压缩机组运行数量满足工况要求，无需增启和减停机组。
[0116]
由此，在本发明实施例当中，与实际运行数量比较，得到n＜n
‑
1，则延时180秒，触发减停机组命令。
[0117]
进一步的，所述运行负荷增减机组计算模型通过负荷占用比load表达，所述负荷占用比load通过如下表达式求得：
[0118][0119]
其中，load代表负荷占用比，r代表平均负荷，r
max
代表最大负荷，r
min
代表最小负荷；
[0120]
所述平均负荷r的计算通过如下表达式求得：
[0121][0122]
其中，n代表实际压缩机运行数量。
[0123]
进一步的，所述触发增启或减停机组命令具体还包括：
[0124]
若load＞95％，延时180秒，触发增启机组命令；
[0125]
若load＜5％，延时180秒，触发减停机组命令。
[0126]
进一步的，所述喘振裕度计算模型具体通过机组平均裕度m表达，所述机组平均裕度m通过如下表达式求得：
[0127][0128]
其中，m代表机组平均裕度，n代表实际压缩机运行数量，m
i
代表各压缩机喘振裕度。
[0129]
进一步的，所述触发增启或减停机组命令具体还包括：
[0130]
若m＜5％，延时180秒，触发减停机组命令。
[0131]
进一步的，所述当需要增启、减停压缩机组时，根据优先级比选结果，自动触发对
应机组的启停机命令，具体包括：
[0132]
(1)采用数字表示优先级，“1”代表优先级最高，数字越大，优先级越低，依次类推，0表示不可用状态；
[0133]
(2)优先级判定采用手动、自动两种方式，其中：
[0134]
手动方式下，可人为写入优先级数字；
[0135]
自动方式下，比较周期内压缩机组累计运行时间，运行时间越长，优先级越低，反之优先级越高；
[0136]
(3)当触发增启压缩机组命令后，在停机状态下的机组中进行优先级比选，优先级越高，即数字越大，优先启动；
[0137]
(4)当触发减停压缩机组命令后，在运行状态下的机组中进行优先级比选，优先级越低，即数字越小，优先停止。
[0138]
由此，在本发明实施例当中，通过优先级比选方法，选择优先停止机组，具体为：
[0139]
当前处于优先级自动模式，压缩机组运行状态正常，
[0140]
进一步的，比较两年内累计运行时间，2#机组大于1#机组，则2#机组优先停止。
[0141]
停机结束30分钟内，跳出压缩机组运行数量计算模型。
[0142]
另一方面，根据运行负荷增减机组计算模型，比较运行机组的平均负荷与最大负荷：
[0143][0144][0145]
得出load＝39.84％，5％＜load＜95％，运行负荷正常，未触发运行负荷增减机组控制逻辑。
[0146]
根据运行机组喘振裕度减停压缩机计算模型
[0147][0148]
得出m＝45.5％，m＞3％，喘振裕度正常，未触发喘振裕度减停压缩机控制逻辑。
[0149]
进一步的，所述压缩机控制系统包括压缩机单机控制系统、压缩机组负荷分配控制系统和站场scada控制系统，所述压缩机单机控制系统、所述压缩机组负荷分配控制系统和所述站场scada控制系统之间能够实现数据互传，获取压气站与压缩机运行相关的实时工况参数。
[0150]
尽管上面已经示出和描述了本发明的实施例，应当理解的是，上述实施例是示例性的，不能解释为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。