CN110513281A - 基于负荷预测的压缩空气系统优化控制方法 - Google Patents
基于负荷预测的压缩空气系统优化控制方法 Download PDFInfo
- Publication number
- CN110513281A CN110513281A CN201910618029.9A CN201910618029A CN110513281A CN 110513281 A CN110513281 A CN 110513281A CN 201910618029 A CN201910618029 A CN 201910618029A CN 110513281 A CN110513281 A CN 110513281A
- Authority
- CN
- China
- Prior art keywords
- compressor
- stage
- layer
- hidden layer
- follows
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04B—POSITIVE-DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS
- F04B49/00—Control, e.g. of pump delivery, or pump pressure of, or safety measures for, machines, pumps, or pumping installations, not otherwise provided for, or of interest apart from, groups F04B1/00 - F04B47/00
- F04B49/06—Control using electricity
- F04B49/065—Control using electricity and making use of computers
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及压缩空气技术领域,尤其涉及一种基于负荷预测的压缩空气系统优化控制方法。包括以下步骤:一,压缩空气的建模及参数辨识:二,利用BP神经网络建立短期供气负荷预测模型,以一周的负荷数据对模型进行训练,预测接下来24小时内每30min的供气负荷;三、以压缩机组全天运行费用最小为目标,综合考虑压缩机组的运行能耗和启动能耗,以基准运行的限制条件为约束,建立优化控制模型并采用动态规划法进行求解。本发明综合考虑压缩机组的运行能耗和启动能耗,结合末端用气需求、排气压力、压缩机最小运行时间和最小停机时间约束条件,建立压缩机组全天运行费用的数学模型,并利用动态规划法进行求解,得到使机组全天费用最小。
Description
技术领域
本发明涉及压缩空气技术领域,尤其涉及一种基于负荷预测的压缩空气系 统优化控制方法。
背景技术
压缩空气系统占据了全国工业用电总量6.5%左右,已是我国工业领域的重 点节能对象。其在工业领域几乎是随处可见,从小型的食品、纺织企业,到大 型的钢铁石化工厂都可以看到其身影,可以说压缩空气对于工业生产已是不可 或缺。然而与其被广泛使用相对应的是系统的高耗能。
从压缩空气系统来看,工业压缩机消耗了整个系统耗电的96%,根据全生命 周期成本理论分析,压缩空气系统全生命周期成本由设备初始投资成本、维护 成本、运行成本三部分组成,以10年作为压缩空气系统的生命周期,设备初始 投资成本和维护成本分别占压缩空气系统全生命周期成本的16%、6%,而运行成 本则高达78%。
随着节能减排的深入推动以及企业为提高经济效益对压缩空气系统节能的 现实需求,需要采取合适的技术对系统进行节能优化,以降低能源消耗。
发明内容
本发明的目的在于克服上述技术的不足,而提供基于负荷预测的压缩空气 系统优化控制方法。
本发明为实现上述目的,采用以下技术方案:一种基于负荷预测的压缩空 气系统优化控制方法,其特征在于:包括以下步骤:
步骤一,压缩空气的建模及参数辨识:
步骤二,利用BP神经网络建立短期供气负荷预测模型,以一周的负荷数据 对模型进行训练,预测接下来24小时内每30min的供气负荷;
步骤三、以压缩机组全天运行费用最小为目标,综合考虑压缩机组的运行 能耗和启动能耗,以基准运行的限制条件为约束,建立优化控制模型并采用动 态规划法进行求解。
优选地,步骤一中,具体步骤如下:
第一步,建立压缩比模型:
气体压缩以等温过程耗功最少,绝热过程耗功最多,多变压缩过程则是介 于二者之间,因此应尽量使压缩过程中气体的温度变化最小,同时因为压缩机 的容积效率是随着单级压缩比的增大而减小的,所以通常采用分级压缩、中间 冷却的方法;
选择不同的压缩比时,整个压缩过程消耗的功不一样,对于z级压缩,各 级压缩比相等时耗功最小,即:
式中,ε为单级压缩的压缩比;p1为进气压力,MPa;pz+1为压缩终了的压 力,MPa;z为压缩机级数;
第二步,建立实际排气量模型:
空压机的理论容积流量与很多因素有关,多级压缩时活塞式压缩机理论排 气量为:
式中,Ql为理论排气量,m3/min;D为一级气缸的直径,m;S为活塞的行 程,m;t为同级气缸的数量;N为曲轴转速,r/min;
从理论上讲,压缩机的排气量就是进入一级气缸的气量,进入的多,排出 的也就多;但在实际情况下,压缩机的排气量会受到进气压力、气缸的余隙容 积、级间泄漏等因素的影响,气缸有效容积减小,实际排气量并不能达到理论 排气量。这些实际因素对实际排气量的影响可用容积效率λ来表征;
压缩机的容积效率为容积系数、压力系数、温度系数以及气密系数的乘积, 即:
λ=λvλpλtλl
式中,λ为容积效率,λv为容积系数,λp为压力系数,λt为温度系数,λl为 气密系数;
活塞式空压机实际排气量Q等于理论排气量与容积效率的乘积,即:
Q=λvλpλtλlQl=λQl
式中,Q为实际排气量,m3/min;Ql为理论排气量,m3/min;
第三步,建立压缩机功率模型:
气体被压缩的过程是非常迅速的,与外界的热传递很少,因此实际的压缩 过程均趋于绝热。多级压缩时,压缩机的绝热功率可按下式计算:
式中,Pad——绝热功率,kW;
pi、pd——进气、排气压力,bar;
κ——气体的绝热指数,对于氮气,κ=1.4;
z——压缩级数。
由于机械损失(如摩擦损失、泄漏损失)的存在,压缩机实际工作时的功 率要大于绝热功率,二者之间的关系可通过绝热效率来衡量,绝热效率为压缩 机的绝热功率与轴功率的比值。因此实际功率为:
式中,P为实际功率,kW;η为压缩机的绝热效率。
优选地,步骤二中,具体步骤如下:
BP神经网络的正向传播和学习过程包括以下步骤:
第一步,输入层有n个神经单元,隐含层有p个神经单元,输出层有q个 神经单元,共有m对学习模式。首先定义以下变量:
输入模式向量希望输出向量隐含层各神经单元的输入激活向量隐含层各神经单元输出向量输出层各神经单元的输入激活向量输出实际 值向量输入层至隐含层的权值Wij,i为目的神经单元,j为源 神经单元;隐含层至输出层的权值Vjt;隐含层各神经单元阈值θj;输出层各神 经单元阈值γt;激活函数f(·);其中i=1,2,···,n;j=1,2,···,p; t=1,2,···,q;k=1,2,···,m;
第二步,输入信号正向传播
隐含层各神经单元的激活值为:
带入激活函数f(x)中,得到隐含层j单元的输出值为:
同理可得输出层第t个神经单元的激活值为lt为:
输出层第t个神经单元的实际输出值为ct,则:
ct=f(lt) (t=1,2,L,q)
第三步,误差逆向传播
第二步的输入信号正向传播过程可以得到神经网络的实际输出值,当实际 输出值的误差较大时,就需要对传递过程进行校正;所谓校正,就是误差的反 向传播,计算过程反过来;
输出层的校正误差为:
式中:t=1,2,L,q;k=1,2,L,m;为希望输出;为实际输出;为输 出层函数的导数;
隐含层各神经单元的校正误差为:
式中:j=1,2,L,p;k=1,2,L,m;
当输出层和隐含层的误差校正量求得后,则可逐层逆向调整输出层 到隐含层、隐含层到输入层的权值;
输出层至隐含层的权值和输出层的阈值校正量为:
式中:ΔVjt为隐含层至输出层的权值校正量;Δγt为输出层各单元阈值校正 量;为隐含层第j神经单元的输出;为输出层第t单元的校正误差;α为隐 含层至输出层的学习系数,0<α<1;
同理,隐含层至输入层的权值和隐含层各单元阈值的校正量为:
式中:ΔWij为输入层至隐含层的权值校正量;Δθj为隐含层各单元阈值校正 量;为输入层第i单元的输入;为隐含层第j神经单元的校正误差;β为输 入层至隐含层的学习系数,0<β<1;
第四步,循环训练学习
实际上就是对上述正向和逆向过程的重复。为使实际输出值的误差最小, BP神经网络的每一组模式一般都需要训练多次,才能使神经网络记住这一训练 模式;
第五步,学习结果判别
每次学习训练结束后都需要对学习结果判别,以检验实际输出值的误差是 否已经小于最小误差ε或者学习次数达到设定次数N,如果是的话则可以结束 学习,否则继续学习过程。在第(2)步中误差修正是针对某一组的学习模式, 对于全部学习模式,全局误差E为:
优选地,步骤三中,具体步骤如下:设压缩机组共有N台压缩机,将一天 的调度周期分为T个阶段,每个阶段的对应序号为t(t=1,2,L,T),t阶段的末 端用气量需求为QD,t,第i台压缩机在阶段t的产气量为Qi,t;
实际运行时,压缩机组总的产气量必须要满足末端的用气需求,即:
其中Qi,t为第i台压缩机的产气量;
由于储气罐的允许工作压力要求在Pmin-PmaxMPa,所以压缩机的排气压力也 存在限制条件:
1.6≤pd≤2.4
根据工艺要求,压缩机的最小连续运行时间和最小连续停机时间均为 30min;
以Ui,t表示压缩机i在时间段t的工作状态,Ui,t=0表示没有运行,Ui,t=1表 示运行,即:
则第i台压缩机的启动能耗表示为:
Ui,t(1-Ui,t-1)Si i=1,2,L,N
压缩机的启停约束可表示为:
式中:T1为压缩机的最小连续停机时间,T2为最小连续运行时间;
综合考虑压缩机的运行能耗和启动能耗,则压缩机组的优化控制就是找到 最优的机组运行方案,使得总费用最小,压缩机组的总费用可表示为:
式中Pi,t为第i台压缩机在时段t运行的用电量,f为阶梯电价;
综上所述,压缩机组的优化控制模型为:
目标函数:
约束条件:
(1)末端用气需求约束
(2)排气压力约束
1.6≤pd≤2.4
(3)最小连续运行时间和最小连续停机时间约束
(4)压缩机状态约束
压缩机组共有N台相同型号的压缩机,以N=4为例,设在t阶段用气需求 为QD,t,它可能由4台压缩机排列组合形成的24=16种运行方案的某一种来满足, 即从第t个阶段到第t+1个阶段允许的决策有16种,为便于求解,将压缩机组 运行方式按顺序编号;
根据最优性原理,建立递推方程如下:
Ft+1,j=min(Ft,x+Sx+Cj(QD,t))
式中Ft+1,j表示从第1个阶段到第t+1个阶段采用第j种机组运行方式的能 耗,Ft,x为t阶段压缩机组采用的最优允许方案——第x种的总能耗,Sx为t阶 段到t+1阶段机组的启动能耗,Cj(QD,t)为t+1阶段采用第j种机组运行方式来 满足用气需求QD,t的运行能耗。0<j,x<2N,0<t<T,N为压缩机台数,N=4;t为 一个调度周期内的阶段数,t=48;
在进行迭代优化过程中,阶段t可能会出现多种满足约束条件的压缩机组运 行方式,这时就需要逐个比较,然后选取总能耗最小的运行方式,进行寻优过 程;
模型迭代寻优过程如下:
(1)选择状态变量:t+1阶段的状态变量就是从种压缩机组运行方式中满足 约束条件的运行方式;
(2)第t阶段的第x状态转移到t阶段的j状态的决策变量需要满足约束 条件;
(3)按式(1)递推;
(4)重复计算过程,在完成T-1阶段的计算后,求解minFT,x对应的运行方 式x后,根据状态转移记录的各阶段的最优路径,从而得到整个周期上压缩机 组的最优运行方式。
本发明的有益效果是:本发明综合考虑压缩机组的运行能耗和启动能耗, 以全天费用最小为目标,结合末端用气需求、排气压力、压缩机最小运行时间 和最小停机时间约束条件,建立压缩机组全天运行费用的数学模型,并利用动 态规划法进行求解,得到使机组全天费用最小。
具体实施方式
下面以较佳实施例详细说明本发明的具体实施方式。
步骤一,压缩空气的建模及参数辨识
压缩空气系统主要包括生产、调压和输送三个环节,其中生产环节主要由 压缩机(组)构成,输送环节主要包括储气罐、调压、管网和阀门等,调压环 节主要通过调压站实现。系统基本工作流程为:气体在压缩机中进行压缩,然 后以较高的压力输出,经过调压站调压后以需求压力通过管网输送至用气末端。 当产气大于用气,气体将储存至储气罐,相反,当用气大于产气,储存在储气 罐中的气体将进行补充。
本发明中压缩空气系统均为活塞式压缩机,以活塞式压缩机为对象,建立 压缩机压缩比、实际排气量、实际功率以及储气罐充放气特性模型。具体步骤 如下:
第一步,建立压缩比模型。气体压缩以等温过程耗功最少,绝热过程耗功 最多,多变压缩过程则是介于二者之间,因此应尽量使压缩过程中气体的温度 变化最小,同时因为压缩机的容积效率是随着单级压缩比的增大而减小的,所 以通常采用分级压缩、中间冷却的方法。
选择不同的压缩比时,整个压缩过程消耗的功不一样。对于z级压缩,各 级压缩比相等时耗功最小,即:
式中,ε为单级压缩的压缩比;p1为进气压力,MPa;pz+1为压缩终了的压 力,MPa;z为压缩机级数。
第二步,建立实际排气量模型。空压机的理论容积流量与很多因素有关, 例如活塞直径,活塞的行程,气缸的数量等,多级压缩时活塞式压缩机理论排 气量为:
式中,Ql为理论排气量,m3/min;D为一级气缸的直径,m;S为活塞的行 程,m;t为同级气缸的数量;N为曲轴转速,r/min。
从理论上讲,压缩机的排气量就是进入一级气缸的气量,进入的多,排出 的也就多。但在实际情况下,压缩机的排气量会受到进气压力、气缸的余隙容 积、级间泄漏等因素的影响,气缸有效容积减小,实际排气量并不能达到理论 排气量。这些实际因素对实际排气量的影响可用容积效率λ来表征。
压缩机的容积效率为容积系数、压力系数、温度系数以及气密系数的乘积, 即:
λ=λvλpλtλl
式中,λ为容积效率,λv为容积系数,λp为压力系数,λt为温度系数,λl为 气密系数。
活塞式空压机实际排气量Q等于理论排气量与容积效率的乘积,即:
Q=λvλpλtλlQl=λQl
式中,Q为实际排气量,m3/min;Ql为理论排气量,m3/min。
第三步,建立压缩机功率模型。
气体被压缩的过程是非常迅速的,与外界的热传递很少,因此实际的压缩 过程均趋于绝热。多级压缩时,压缩机的绝热功率可按下式计算:
式中,Pad——绝热功率,kW;
pi、pd——进气、排气压力,bar;
κ——气体的绝热指数,对于氮气,κ=1.4;
z——压缩级数。
由于机械损失(如摩擦损失、泄漏损失)的存在,压缩机实际工作时的功 率要大于绝热功率,二者之间的关系可通过绝热效率来衡量,绝热效率为压缩 机的绝热功率与轴功率的比值。因此实际功率为:
式中,P为实际功率,kW;η为压缩机的绝热效率。
步骤二,利用BP神经网络建立短期供气负荷预测模型,以一周的负荷数据 对模型进行训练,预测接下来24小时内每30min的供气负荷。
BP神经网络训练学习时,首先是正向预测,输入层接受外部输入的信号(历 史负荷数据),由各节点传递给直接相连的、对应的隐含层各节点,隐含层将信 号处理后传递给输出层,当预测结果与实际值的误差大于规定时,则将误差反 向传递:预测误差经隐含层反传回输入层并且分摊给各节点,将获得的各层单 元的误差信号作为修正各单元连接权的依据。重复正向、反向过程直至预测结 果符合要求为止,也是BP神经网络不断学习训练的过程。BP神经网络的正向传 播和学习过程包括以下步骤:
第一步,输入层有n个神经单元,隐含层有p个神经单元,输出层有q个 神经单元,共有m对学习模式。首先定义以下变量:
输入模式向量希望输出向量隐含层各神经单元的输入激活向量隐含层各神经单元输出向量 输出层各神经单元的输入激活向量输出实际 值向量输入层至隐含层的权值Wij,i为目的神经单元,j为源 神经单元;隐含层至输出层的权值Vjt;隐含层各神经单元阈值θj;输出层各神 经单元阈值γt;激活函数f(·)。其中i=1,2,···,n;j=1,2,···,p; t=1,2,···,q;k=1,2,···,m。
第二步,输入信号正向传播
隐含层各神经单元的激活值为:
带入激活函数f(x)中,得到隐含层j单元的输出值为:
同理可得输出层第t个神经单元的激活值为lt为:
输出层第t个神经单元的实际输出值为ct,则:
ct=f(lt) (t=1,2,L,q)
第三步,误差逆向传播
第二步的输入信号正向传播过程可以得到神经网络的实际输出值,当实际 输出值的误差较大时,就需要对传递过程进行校正。所谓校正,就是误差的反 向传播,计算过程反过来。
输出层的校正误差为:
式中:t=1,2,L,q;k=1,2,L,m;为希望输出;为实际输出;为输 出层函数的导数。
隐含层各神经单元的校正误差为:
式中:j=1,2,L,p;k=1,2,L,m。
当输出层和隐含层的误差校正量求得后,则可逐层逆向调整输出层 到隐含层、隐含层到输入层的权值。
输出层至隐含层的权值和输出层的阈值校正量为:
式中:ΔVjt为隐含层至输出层的权值校正量;Δγt为输出层各单元阈值校正 量;为隐含层第j神经单元的输出;为输出层第t单元的校正误差;α为隐 含层至输出层的学习系数,0<α<1。
同理,隐含层至输入层的权值和隐含层各单元阈值的校正量为:
式中:ΔWij为输入层至隐含层的权值校正量;Δθj为隐含层各单元阈值校正 量;为输入层第i单元的输入;为隐含层第j神经单元的校正误差;β为输 入层至隐含层的学习系数,0<β<1。
第四步,循环训练学习
实际上就是对上述正向和逆向过程的重复。为使实际输出值的误差最小, BP神经网络的每一组模式一般都需要训练多次,才能使神经网络记住这一训练 模式。
第五步,学习结果判别
每次学习训练结束后都需要对学习结果判别,以检验实际输出值的误差是 否已经小于最小误差ε或者学习次数达到设定次数N,如果是的话则可以结束 学习,否则继续学习过程。在第(2)步中误差修正是针对某一组的学习模式, 对于全部学习模式,全局误差E为:
步骤三,以压缩机组全天运行费用最小为目标,综合考虑压缩机组的运行 能耗和启动能耗,以基准运行的限制条件为约束,建立优化控制模型并采用动 态规划法进行求解。
设压缩机组共有N台压缩机,将一天的调度周期分为T个阶段,每个阶段 的对应序号为t(t=1,2,L,T),t阶段的末端用气量需求为QD,t,第i台压缩机在 阶段t的产气量为Qi,t。
实际运行时,压缩机组总的产气量必须要满足末端的用气需求,即:
其中Qi,t为第i台压缩机的产气量。
由于储气罐的允许工作压力要求在Pmin-PmaxMPa,所以压缩机的排气压力也 存在限制条件:
1.6≤pd≤2.4
如果压缩机启停过于频繁,不仅会对压缩机的运行调度带来困难,也会造 成压缩机的疲劳损坏、对电路造成冲击,因此压缩机频繁启停往往是不允许的。 为避免这个问题,通常规定压缩机的最小连续停机时间和最小连续运行时间。 最小连续停机时间即压缩机停机后,至少要连续停机的时间;最小连续运行时 间即压缩机开始运行后,最小的连续运行时间。本发明研究的压缩机组,根据 工艺要求,压缩机的最小连续运行时间和最小连续停机时间均为30min。
以Ui,t表示压缩机i在时间段t的工作状态,Ui,t=0表示没有运行,Ui,t=1表 示运行,即:
则第i台压缩机的启动能耗表示为:
Ui,t(1-Ui,t-1)Si i=1,2,L,N
压缩机的启停约束可表示为:
式中:T1为压缩机的最小连续停机时间,T2为最小连续运行时间。
综合考虑压缩机的运行能耗和启动能耗,则压缩机组的优化控制就是找到 最优的机组运行方案,使得总费用最小,压缩机组的总费用可表示为:
式中Pi,t为第i台压缩机在时段t运行的用电量,f为阶梯电价。
综上所述,压缩机组的优化控制模型为:
目标函数:
约束条件:
(1)末端用气需求约束
(2)排气压力约束
1.6≤pd≤2.4
(3)最小连续运行时间和最小连续停机时间约束
(4)压缩机状态约束
压缩机组共有N台相同型号的压缩机,以N=4为例,设在t阶段用气需求 为QD,t,它可能由4台压缩机排列组合形成的24=16种运行方案的某一种来满足, 即从第t个阶段到第t+1个阶段允许的决策有16种,为便于求解,将压缩机组 运行方式按顺序编号。
根据最优性原理,建立递推方程如下:
Ft+1,j=min(Ft,x+Sx+Cj(QD,t))
式中Ft+1,j表示从第1个阶段到第t+1个阶段采用第j种机组运行方式的能 耗,Ft,x为t阶段压缩机组采用的最优允许方案——第x种的总能耗,Sx为t阶 段到t+1阶段机组的启动能耗,Cj(QD,t)为t+1阶段采用第j种机组运行方式来 满足用气需求QD,t的运行能耗。0<j,x<2N,0<t<T,N为压缩机台数,N=4;t为 一个调度周期内的阶段数,t=48。
在进行迭代优化过程中,阶段t可能会出现多种满足约束条件的压缩机组运 行方式,这时就需要逐个比较,然后选取总能耗最小的运行方式,进行寻优过 程。
模型迭代寻优过程如下:
(1)选择状态变量:t+1阶段的状态变量就是从种压缩机组运行方式中满 足约束条件的运行方式;
(2)第t阶段的第x状态转移到t阶段的j状态的决策变量需要满足约束 条件;
(3)按式(1)递推;
(4)重复计算过程,在完成T-1阶段的计算后,求解minFT,x对应的运行方 式x后,根据状态转移记录的各阶段的最优路径,从而得到整个周期上压缩机 组的最优运行方式。
本发明综合考虑压缩机组的运行能耗和启动能耗,以全天费用最小为目标, 结合末端用气需求、排气压力、压缩机最小运行时间和最小停机时间约束条件, 建立压缩机组全天运行费用的数学模型,并利用动态规划法进行求解,得到使 机组全天费用最小的优化控制方法。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通 技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰, 这些改进和润饰也应视为本发明的保护范围。
Claims (4)
1.一种基于负荷预测的压缩空气系统优化控制方法,其特征在于:包括以下步骤:
步骤一,压缩空气的建模及参数辨识:
步骤二,利用BP神经网络建立短期供气负荷预测模型,以一周的负荷数据对模型进行训练,预测接下来24小时内每30min的供气负荷;
步骤三、以压缩机组全天运行费用最小为目标,综合考虑压缩机组的运行能耗和启动能耗,以基准运行的限制条件为约束,建立优化控制模型并采用动态规划法进行求解。
2.根据权利要求1所述的基于负荷预测的压缩空气系统优化控制方法,其特征在于,步骤一中,具体步骤如下:
第一步,建立压缩比模型:
气体压缩以等温过程耗功最少,绝热过程耗功最多,多变压缩过程则是介于二者之间,因此应尽量使压缩过程中气体的温度变化最小,同时因为压缩机的容积效率是随着单级压缩比的增大而减小的,所以通常采用分级压缩、中间冷却的方法;
选择不同的压缩比时,整个压缩过程消耗的功不一样,对于z级压缩,各级压缩比相等时耗功最小,即:
式中,ε为单级压缩的压缩比;p1为进气压力,MPa;pz+1为压缩终了的压力,MPa;z为压缩机级数;
第二步,建立实际排气量模型:
空压机的理论容积流量与很多因素有关,多级压缩时活塞式压缩机理论排气量为:
式中,Ql为理论排气量,m3/min;D为一级气缸的直径,m;S为活塞的行程,m;t为同级气缸的数量;N为曲轴转速,r/min;
从理论上讲,压缩机的排气量就是进入一级气缸的气量,进入的多,排出的也就多;但在实际情况下,压缩机的排气量会受到进气压力、气缸的余隙容积、级间泄漏等因素的影响,气缸有效容积减小,实际排气量并不能达到理论排气量。这些实际因素对实际排气量的影响可用容积效率λ来表征;
压缩机的容积效率为容积系数、压力系数、温度系数以及气密系数的乘积,即:
λ=λvλpλtλl
式中,λ为容积效率,λv为容积系数,λp为压力系数,λt为温度系数,λl为气密系数;
活塞式空压机实际排气量Q等于理论排气量与容积效率的乘积,即:
Q=λvλpλtλlQl=λQl
式中,Q为实际排气量,m3/min;Ql为理论排气量,m3/min;
第三步,建立压缩机功率模型:
气体被压缩的过程是非常迅速的,与外界的热传递很少,因此实际的压缩过程均趋于绝热。多级压缩时,压缩机的绝热功率可按下式计算:
式中,Pad——绝热功率,kW;
pi、pd——进气、排气压力,bar;
κ——气体的绝热指数,对于氮气,κ=1.4;
z——压缩级数。
由于机械损失(如摩擦损失、泄漏损失)的存在,压缩机实际工作时的功率要大于绝热功率,二者之间的关系可通过绝热效率来衡量,绝热效率为压缩机的绝热功率与轴功率的比值。因此实际功率为:
式中,P为实际功率,kW;η为压缩机的绝热效率。
3.根据权利要求1所述的基于负荷预测的压缩空气系统优化控制方法,其特征在于,步骤二中,具体步骤如下:
BP神经网络的正向传播和学习过程包括以下步骤:
第一步,输入层有n个神经单元,隐含层有p个神经单元,输出层有q个神经单元,共有m对学习模式。首先定义以下变量:
输入模式向量希望输出向量隐含层各神经单元的输入激活向量隐含层各神经单元输出向量输出层各神经单元的输入激活向量输出实际值向量输入层至隐含层的权值Wij,i为目的神经单元,j为源神经单元;隐含层至输出层的权值Vjt;隐含层各神经单元阈值θj;输出层各神经单元阈值γt;激活函数f(·);其中i=1,2,···,n;j=1,2,···,p;t=1,2,···,q;k=1,2,···,m;
第二步,输入信号正向传播
隐含层各神经单元的激活值为:
带入激活函数f(x)中,得到隐含层j单元的输出值为:
同理可得输出层第t个神经单元的激活值为lt为:
输出层第t个神经单元的实际输出值为ct,则:
ct=f(lt) (t=1,2,L,q)
第三步,误差逆向传播
第二步的输入信号正向传播过程可以得到神经网络的实际输出值,当实际输出值的误差较大时,就需要对传递过程进行校正;所谓校正,就是误差的反向传播,计算过程反过来;
输出层的校正误差为:
式中:t=1,2,L,q;k=1,2,L,m;为希望输出;为实际输出;为输出层函数的导数;
隐含层各神经单元的校正误差为:
式中:j=1,2,L,p;k=1,2,L,m;
当输出层和隐含层的误差校正量求得后,则可逐层逆向调整输出层到隐含层、隐含层到输入层的权值;
输出层至隐含层的权值和输出层的阈值校正量为:
式中:ΔVjt为隐含层至输出层的权值校正量;Δγt为输出层各单元阈值校正量;为隐含层第j神经单元的输出;为输出层第t单元的校正误差;α为隐含层至输出层的学习系数,0<α<1;
同理,隐含层至输入层的权值和隐含层各单元阈值的校正量为:
式中:ΔWij为输入层至隐含层的权值校正量;Δθj为隐含层各单元阈值校正量;为输入层第i单元的输入;为隐含层第j神经单元的校正误差;β为输入层至隐含层的学习系数,0<β<1;
第四步,循环训练学习
实际上就是对上述正向和逆向过程的重复。为使实际输出值的误差最小,BP神经网络的每一组模式一般都需要训练多次,才能使神经网络记住这一训练模式;
第五步,学习结果判别
每次学习训练结束后都需要对学习结果判别,以检验实际输出值的误差是否已经小于最小误差ε或者学习次数达到设定次数N,如果是的话则可以结束学习,否则继续学习过程。在第(2)步中误差修正是针对某一组的学习模式,对于全部学习模式,全局误差E为:
。
4.根据权利要求1所述的基于负荷预测的压缩空气系统优化控制方法,其特征在于,步骤三中,具体步骤如下:设压缩机组共有N台压缩机,将一天的调度周期分为T个阶段,每个阶段的对应序号为t(t=1,2,L,T),t阶段的末端用气量需求为QD,t,第i台压缩机在阶段t的产气量为Qi,t;
实际运行时,压缩机组总的产气量必须要满足末端的用气需求,即:
其中Qi,t为第i台压缩机的产气量;
由于储气罐的允许工作压力要求在Pmin-PmaxMPa,所以压缩机的排气压力也存在限制条件:
1.6≤pd≤2.4
根据工艺要求,压缩机的最小连续运行时间和最小连续停机时间均为30min;
以Ui,t表示压缩机i在时间段t的工作状态,Ui,t=0表示没有运行,Ui,t=1表示运行,即:
则第i台压缩机的启动能耗表示为:
Ui,t(1-Ui,t-1)Si i=1,2,L,N
压缩机的启停约束可表示为:
式中:T1为压缩机的最小连续停机时间,T2为最小连续运行时间;
综合考虑压缩机的运行能耗和启动能耗,则压缩机组的优化控制就是找到最优的机组运行方案,使得总费用最小,压缩机组的总费用可表示为:
式中Pi,t为第i台压缩机在时段t运行的用电量,f为阶梯电价;
综上所述,压缩机组的优化控制模型为:
目标函数:
约束条件:
(1)末端用气需求约束
(2)排气压力约束
1.6≤pd≤2.4
(3)最小连续运行时间和最小连续停机时间约束
(4)压缩机状态约束
压缩机组共有N台相同型号的压缩机,以N=4为例,设在t阶段用气需求为QD,t,它可能由4台压缩机排列组合形成的24=16种运行方案的某一种来满足,即从第t个阶段到第t+1个阶段允许的决策有16种,为便于求解,将压缩机组运行方式按顺序编号;
根据最优性原理,建立递推方程如下:
Ft+1,j=min(Ft,x+Sx+Cj(QD,t))
式中Ft+1,j表示从第1个阶段到第t+1个阶段采用第j种机组运行方式的能耗,Ft,x为t阶段压缩机组采用的最优允许方案——第x种的总能耗,Sx为t阶段到t+1阶段机组的启动能耗,Cj(QD,t)为t+1阶段采用第j种机组运行方式来满足用气需求QD,t的运行能耗。0<j,x<2N,0<t<T,N为压缩机台数,N=4;t为一个调度周期内的阶段数,t=48;
在进行迭代优化过程中,阶段t可能会出现多种满足约束条件的压缩机组运行方式,这时就需要逐个比较,然后选取总能耗最小的运行方式,进行寻优过程;
模型迭代寻优过程如下:
(1)选择状态变量:t+1阶段的状态变量就是从种压缩机组运行方式中满足约束条件的运行方式;
(2)第t阶段的第x状态转移到t阶段的j状态的决策变量需要满足约束条件;
(3)按式(1)递推;
(4)重复计算过程,在完成T-1阶段的计算后,求解minFT,x对应的运行方式x后,根据状态转移记录的各阶段的最优路径,从而得到整个周期上压缩机组的最优运行方式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910618029.9A CN110513281A (zh) | 2019-07-10 | 2019-07-10 | 基于负荷预测的压缩空气系统优化控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910618029.9A CN110513281A (zh) | 2019-07-10 | 2019-07-10 | 基于负荷预测的压缩空气系统优化控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110513281A true CN110513281A (zh) | 2019-11-29 |
Family
ID=68622565
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910618029.9A Pending CN110513281A (zh) | 2019-07-10 | 2019-07-10 | 基于负荷预测的压缩空气系统优化控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110513281A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112539161A (zh) * | 2020-12-09 | 2021-03-23 | 上海安悦节能技术有限公司 | 一种基于数据驱动的空压机系统控制方法 |
CN113391656A (zh) * | 2021-06-07 | 2021-09-14 | 国家石油天然气管网集团有限公司西气东输分公司 | 一种压缩机组智能增减控制方法 |
CN113464415A (zh) * | 2021-09-03 | 2021-10-01 | 蘑菇物联技术(深圳)有限公司 | 控制空压站的空压机的方法、设备和计算机存储介质 |
CN114576554A (zh) * | 2022-05-06 | 2022-06-03 | 浙江浙能航天氢能技术有限公司 | 一种基于神经网络的氢气压缩模块控制方法及装置 |
CN114841075A (zh) * | 2022-05-19 | 2022-08-02 | 华北电力大学 | 一种电-气综合能源系统的最优能流计算建模方法 |
WO2023082427A1 (zh) * | 2021-11-11 | 2023-05-19 | 烟台杰瑞石油服务集团股份有限公司 | 故障诊断方法、诊断模型构建方法、装置、设备和介质 |
CN116578074A (zh) * | 2023-07-14 | 2023-08-11 | 德耐尔能源装备有限公司 | 一种集装箱氮气压缩机组的集中监测控制方法及系统 |
CN117311430A (zh) * | 2023-11-29 | 2023-12-29 | 广东鑫钻节能科技股份有限公司 | 一种节能的数字能源空压站控制系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102032443A (zh) * | 2009-09-27 | 2011-04-27 | 中国石油化工股份有限公司 | 气柜负荷预测设备和方法以及瓦斯调度系统和方法 |
CN108734649A (zh) * | 2017-04-24 | 2018-11-02 | 英特尔公司 | 神经网络训练机构 |
-
2019
- 2019-07-10 CN CN201910618029.9A patent/CN110513281A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102032443A (zh) * | 2009-09-27 | 2011-04-27 | 中国石油化工股份有限公司 | 气柜负荷预测设备和方法以及瓦斯调度系统和方法 |
CN108734649A (zh) * | 2017-04-24 | 2018-11-02 | 英特尔公司 | 神经网络训练机构 |
Non-Patent Citations (6)
Title |
---|
徐少明等: "《空气压缩机实用技术》", 30 April 1994, 机械工业出版社 * |
李丽侦: "基于神经网络的压缩空气用能预测模型研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
李艳: "基于BP神经网络的本体映射方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
李长春: "空压机群节能监控系统在电厂的应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
白文翔: "活塞式空压机在山地物探爆破孔施工中的应用研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
陆良福: "《炼油过程及设备》", 31 May 1993, 中国石化出版社 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112539161A (zh) * | 2020-12-09 | 2021-03-23 | 上海安悦节能技术有限公司 | 一种基于数据驱动的空压机系统控制方法 |
CN113391656A (zh) * | 2021-06-07 | 2021-09-14 | 国家石油天然气管网集团有限公司西气东输分公司 | 一种压缩机组智能增减控制方法 |
CN113391656B (zh) * | 2021-06-07 | 2022-06-28 | 国家石油天然气管网集团有限公司西气东输分公司 | 一种压缩机组智能增减控制方法 |
CN113464415A (zh) * | 2021-09-03 | 2021-10-01 | 蘑菇物联技术(深圳)有限公司 | 控制空压站的空压机的方法、设备和计算机存储介质 |
CN113464415B (zh) * | 2021-09-03 | 2021-12-03 | 蘑菇物联技术(深圳)有限公司 | 控制空压站的空压机的方法、设备和计算机存储介质 |
WO2023082427A1 (zh) * | 2021-11-11 | 2023-05-19 | 烟台杰瑞石油服务集团股份有限公司 | 故障诊断方法、诊断模型构建方法、装置、设备和介质 |
US11886178B2 (en) | 2021-11-11 | 2024-01-30 | Yantai Jereh Oilfield Services Group Co., Ltd. | Fault diagnosis method, method for building fault diagnosis model, equipment, device and medium |
CN114576554B (zh) * | 2022-05-06 | 2022-08-02 | 浙江浙能航天氢能技术有限公司 | 一种基于神经网络的氢气压缩模块控制方法及装置 |
CN114576554A (zh) * | 2022-05-06 | 2022-06-03 | 浙江浙能航天氢能技术有限公司 | 一种基于神经网络的氢气压缩模块控制方法及装置 |
CN114841075A (zh) * | 2022-05-19 | 2022-08-02 | 华北电力大学 | 一种电-气综合能源系统的最优能流计算建模方法 |
CN116578074A (zh) * | 2023-07-14 | 2023-08-11 | 德耐尔能源装备有限公司 | 一种集装箱氮气压缩机组的集中监测控制方法及系统 |
CN116578074B (zh) * | 2023-07-14 | 2023-09-26 | 德耐尔能源装备有限公司 | 一种集装箱氮气压缩机组的集中监测控制方法及系统 |
CN117311430A (zh) * | 2023-11-29 | 2023-12-29 | 广东鑫钻节能科技股份有限公司 | 一种节能的数字能源空压站控制系统 |
CN117311430B (zh) * | 2023-11-29 | 2024-02-20 | 广东鑫钻节能科技股份有限公司 | 一种节能的数字能源空压站控制系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110513281A (zh) | 基于负荷预测的压缩空气系统优化控制方法 | |
Li et al. | A hybrid optimization-based scheduling strategy for combined cooling, heating, and power system with thermal energy storage | |
WO2021189739A1 (zh) | 一种基于生产计划的钢铁企业氧气负荷预测方法 | |
CN111256294A (zh) | 一种基于模型预测的冷水机组组合运行优化控制方法 | |
CN101155995A (zh) | 使多个压缩机机组的运行最优化的方法和用于此方法的装置 | |
CN102734890B (zh) | 中央空调冷冻水的模糊控制方法、装置及中央空调系统 | |
CN102242868B (zh) | 一种工业装置蒸汽管网优化运行方法 | |
CN113032935B (zh) | 一种大型并联天然气管网优化运行模型及求解方法 | |
CN113129164B (zh) | 一种天然气管网天然气流量压力调度决策指标的计算方法 | |
CN113537644B (zh) | 一种多空压站动态协同优化调控系统及方法 | |
CN111832138B (zh) | 一种区域管网拓扑优化方法 | |
CN109917656A (zh) | 基于工艺介质多温度目标的循环冷却水最小压差节能控制系统及方法 | |
CN116542370A (zh) | 一种考虑碳捕集与碳交易的园区低碳经济运行方法 | |
CN111340305A (zh) | 一种建筑运行能耗预测方法 | |
CN112132481A (zh) | 一种区域能效综合评价方法 | |
CN108734419A (zh) | 一种基于知识迁移的高炉煤气调度系统建模方法 | |
CN112365129A (zh) | 基于交叉超效率ccr模型的综合能源系统综合效率评价方法 | |
CN114925928B (zh) | 一种智能动态调控的地下综合体浅层地热能利用、存储方法及其系统 | |
CN116777154A (zh) | 基于分布式近端策略梯度算法的综合能源系统调度方法 | |
Mei et al. | Energy efficiency prediction of screw chillers on BP neural network optimized by improved genetic algorithm | |
CN113420397B (zh) | 储能电站动态运行调控方法 | |
CN107045583B (zh) | 一种海洋平台发电优化和二氧化碳减排方法 | |
CN114781750A (zh) | 一种低输量原油管道间歇输送优化运行模型及求解方法 | |
CN113555907A (zh) | 一种考虑非理想通信情况的分布式多能系统滚动优化控制方法 | |
Fadnes et al. | Artificial neural network model for predicting CO2 heat pump behaviour in domestic hot water and space heating systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191129 |