1.一种基于炼油厂全厂调度离散时间建模方法，其特征在于，所述建模方法是一种考虑生产模式过渡过程的炼油生产过程和存储控制的离散时间调度建模方法。把整个炼油厂系统划分为原油供应、炼油生产、成品油调和交付三个部分，基于离散时间，从生产装置的运行模式、生产装置运行模式的过渡过程的角度进行建模，按以下步骤建模：
步骤1：问题描述
首先对建模对象进行分析，分析出模型中的决策变量；将一个典型的炼油厂系统划分为三部分：第一部分是原油供应，假定来自原油储罐的原油供应是充足的；第二部分是生产装置，包括常压蒸馏装置(ATM)、减压蒸馏装置(VDU)、流化催化裂化装置(FCCU)、加氢精制装置(HTU)、加氢脱硫装置(HDS)、催化重整装置(RF)、醚化装置(ETH)和甲基叔丁基醚装置(MTBE)；第三部分是成品油调和及交付，该建模对象中有一种供应原油和八种成品油，其中包括五种汽油(JIV93、JIV97、GIII90、GIII93、GIII97)和三种柴油(GIII0、GIIIM10、GIV0)，假定成品油均存放于储油罐中，以最大限度地满足订单要求同时最小化总的生产成本代价为调度优化目标。
步骤2：操作模式定义
所述步骤1中不同炼油生产装置的操作模式不同，ATM和VDU装置具有汽油模式(G)和柴油模式(D)两种操作模式；FCCU、HDS和ETH装置具有汽油—汽油模式(GG)、汽油—柴油模式(GD)、柴 油—汽油模式(DG)和柴油—柴油模式(DD)四种操作模式；HTU1和HTU2装置具有苛刻操作模式(H)和温和操作模式(M)两种操作模式；RF和MTBE都仅有一种操作模式。
步骤3：调度模型采用离散时间表述。
先不考虑模式切换过渡过程，则在每个时刻点上，生产装置的操作模式和输入物料量、用于调和的组分油油量和成品油的交付量都是确定的，yu,m,t表征生产装置u在时间间隔t内的操作模式是否为m；
在调度模型中引入模式切换过渡过程后，需额外增加的决策变量：xu,m,m',t表征生产装置u在不同的时间间隔t内是否处于从操作模式m到m'的模式切换过渡过程；Cu,m,m',t表征生产装置u在时间间隔t-1与时间间隔t之间是否有从操作模式m到m'的切换。
步骤4：问题公式化‐构建数学模型
基于离散时间表示的炼油厂全厂调度模型可构建为混合整数非线性规划(MINLP)数学模型，其具有各类必要的约束条件，包括：
A运行模式切换约束：运行模式变量约束、模式切换变量约束、过渡过程变量约束和过渡过程保持时间约束；B物料平衡及容量、组分油调和及成品油交付的约束；
步骤5：目标函数‐构建的调度模型
炼油厂调度问题的目标函数是最小化炼油厂的生产成本、物料储存成本以及订单缺货惩罚费用，目标函数的数学表达式如下：
minf ＝ minΣT(QIATM,tOPC+ΣuΣmΣm'xu,m,m',tQIu,ttOpCostu,m,m'+ΣuΣm'yu,m',t(1-Σmxu,m,m',t)QIu,tOpCostu,m')+Σtα(ΣoINVo,t+ΣocINVoc,t)+ΣlΣoβl(Rl,o-ΣtDl,o,t) (24)QIATM,t为生产装置ATM处在时间间隔t内的输入流量；
OPC为原油c的价格；
tOpCostu,m,m'为生产装置u在操作模式从m到m'的过渡过程中的操作成本；
OpCostu,m为生产装置u处于操作模式m的操作成本；
α为每个时间间隔组分油和成品油的罐存成本；
βl为每个时间间隔对订单l交付延迟的惩罚因子。
目标函数式中第一项是购买原油的成本费用，第二项和第三项是生产装置在过渡过程和稳态运行过程中操作成本，第四项是物料储存费用，第五项是订单缺货惩罚。
混合整数非线性规划模型如下：
(P0)：
min f＝ΣT(QIATM,tOPC+ΣuΣmΣm'xu,m,m',tQIu,ttOpCostu,m,m'+ΣuΣm'yu,m',t(1-Σmxu,m,m',t)QIu,tOpCostu,m')+Σtα(ΣoINVo,t+ΣocINVoc,t)+ΣlΣoβl(Rl,o-ΣtDl,o,t)
步骤6：模型线性化
以上构建的调度模型(P0)中涉及有双线性项和三线性项，双线性项是一个二进制变量与一个连续变量的乘积，三线性项是两个二进制变量与一个连续变量的乘积，可以通过引入额外的辅助变量将这些项线性化；
具体来说，步骤4中约束条件(7)和步骤5目标函数中都涉及 到相同的双线性项和
三线性项；双线性项是xu,m,m',tQIu,t，其中，xu,m,m',t是二进制变量，QIu,t是连续变量；三线性项是yu,m',t(1-Σmxu,m,m',t)QIu,t，由xu,m,m',t的定义知，xu,m,m',t＝1表示生产装置u在时间间隔t内处于从m至m'的切换过渡过程。因为生产装置在时间间隔t-1上的操作模式唯一，故Σmxu,m,m',t的值不超过1，所以(1-Σmxu,m,m',t)可以视为二进制变量，yu,m',t是二进制变量，QIu,t是连续变量。
步骤7：线性化后的约束和目标函数
如步骤5、6所述，生产装置流出口物料平衡约束和目标函数可重新书写如下：
min f'＝minΣT(QIATM,tOPC+ΣuΣmΣm'xQIu,m,m',ttOpCostu,m,m'+ΣuΣm'xyQIu,m',tOpCostu,m')+Σtα(ΣoINVo,t+ΣocINVoc,t)+ΣlΣoβl(Rl,o-ΣtDl,o,t) (24’)则最终重构的离散时间混合整数非线性规划调度模型如下：
(P1)：
min f'＝ ΣT(QIATM,tOPC+ΣuΣmΣm'xQIu,m,m',ttOpCostu,m,m'+ΣuΣm'xyQIu,m',tOpCostu,m')+Σtα(ΣoINVo,t+ΣocINVoc,t)+ΣlΣoβl(Rl,o-ΣtDl,o,t)
2.根据权利要求1所述一种基于炼油厂全厂调度离散时间建模方法，其特征在于，所述步骤1模型中的决策变量有：
a)每个生产装置上运行的操作模式的次序及其起止时间，包括yu，m，t、Cu，m，m′，t、xu，m，m′，t；
b)每个生产装置在每个时间点的产率(负荷)，包括QOu，s，t、QOu，oi，t；
c)每个时间点上用于调和的组分油的品种和数量，包括Qoc，o，t、QOoc，t；
d)每个时间点上存储或交付的成品油的品种和数量，包括Dl，o，t、QIo，t。
模型中根据外部信息可确定的参数有：
a)每个生产装置的运行操作模式Mu和相应的过渡过程；
b)每个生产装置在稳态运行时的收率Yieldu，s，m及过渡过程中的收率tYieldu，s，m，m′；
c)每个生产装置在稳态运行时的运行成本OpCostu，m及过渡过程中的运行成本
tOpCostu，m，m′；
d)每个过渡过程的时长(过渡过程持续时间)TTu，m，m′；
e)成品油的关键特性值范围，包括
f)每个订单要求的交付时间及所需成品油油量，包括Tl1、Tl2、Rl,o；
g)生产装置的最小输入流量值 和最大输入流量值
min max
h)所有储油罐的容量限度，包括 INVo 、INVo ；
i)组分油的最小调和比例值 及最大调和比例值
j)原油价格OPC；
k)物料储存成本α和订单缺货惩罚值βl；
l)调度时间跨度范围T。
3.根据权利要求1所述一种基于炼油厂全厂调度离散时间建模方法，其特征在于，所述步骤2中不同生产装置的操作模式如下：所述生产装置具有多种操作模式，如表1所示，表1 生产装置的操作模式
其中，
a)ATM和VDU
对于ATM和VDU，有两种分馏操作运行模式：汽油模式(G)和柴油模式(D)；汽油模式下装置会尽可能多的产出汽油馏分，柴油模式下装置会尽可能多的产出柴油馏分；
b)FCCU
FCCU有两个主要部分：反应部分和分馏部分，与ATM、VDU装置相似，这两个部分的操作模式也分为汽油模式和柴油模式，同样，汽油模式下装置会尽可能多的产出汽油馏分，柴油模式下装置会尽可能多的产出柴油馏分；因此将两部分相结合，FCCU共有四个操作模式，分别命名为：汽油-汽油模式(GG)，汽油-柴油模式(GD)，柴油-汽油模式(DG)，柴油-柴油模式(DD)；FCCU的模式切换过渡过程如表2所示，
表2 FCCU的模式切换过渡过程
c)HDS和ETH
对HDS来说，产出物的收率和关键性能指标都与来自FCCU的待处理物料的种类相关。
如果FCCU的操作模式发生变化，则FCCU的产出物种类会发生变化，相应地HDS的生产处理过程也要发生变化，即进行操作模式切换。将这些不同的处理过程定义成不同的操作模式，模式名称与FCCU装置的模式名称完全相同；
ETH装置的生产过程与上述相似，产出物的收率和关键性能指标都与来自HDS的待处理物料的种类相关，采用与分析HDS相同的方法定义ETH的操作模式；
d)HTU1和HTU2
对HTU1和HTU2来说，共有两种操作模式：苛刻操作模式(H)和温和操作模式(M)；与温和操作模式相比，苛刻操作模式产出的组分油有更低的含硫量和更高的十六烷值；相应的，苛刻操作模式的运行成本也更高；
e)RF和MTBE
RF和MTBE都仅有一种操作模式，假定在过渡过程中，运行成本的变化与收率的变化保持一致，采用积分后求平均的方法得到过渡过程的固定运行成本与收率；与稳态运行过程相比，过渡过程的操作成本更高而收率更低。
4.根据权利要求1所述一种基于炼油厂全厂调度离散时间建模方法，其特征在于，所述步骤4构建数学模型，是将基于离散时间表示的炼油厂全厂调度模型构建为混合整数非线性规划(MINLP)数学模 型，其中包括：
A运行模式切换约束
A.1运行模式变量约束
任何一种生产装置在任何时间只能有一种运行模式，
其中，yu，m，t＝1表示生产装置u在时间间隔t的运行模式是否为m；U为生产装置的集合；T为时间间隔的集合；
A.2模式切换变量约束
其中，Cu，m，m′，t＝1表示生产装置u在时间间隔t-1与时间间隔t之间运行模式从m切换至m′；则当Cu，m，m′，t＝1时，yu，m，(t-1)和yu，m′，t均为1；Mu为生产装置u的生产模式的集合；
A.3过渡过程变量约束；
其中，xu，m，m′，t表示生产装置u在时间间隔t内是否处于操作模式从m至m′的切换过渡过程；如果xu，m，m′，t＝1，模式切换一定发生在从时间间隔t-TTu，m，m′+1到时间间隔t的范围内；如果xu，m，m′，t＝0，在从时间间隔t-TTu，m，m′+1到时间间隔t的范围内就一定没有模式切换；如果式(4)中的m＝m′，则TTu，m，m′＝0，因此xu，m，m′，t＝0，表明，若未发生模式切换，就无过渡过程；表示为
过渡过程必须在整个调度时间范围内完成，这意味着在最后一个时间间隔内所有生产装置都不处于过渡过程；
A.4过渡过程保持时间约束
在一个过渡过程结束前，时间间隔都处于过渡过程中，不应有新的模式切换发生；
在某时间间隔内，生产装置u发生了从运行模式m到运行模式m′的模式切换，在过渡过程期间的时间间隔内不可以发生新的模式切换，直至生产装置u稳定于运行模式m′的第二个时间间隔才能发生下一次模式切换或者继续保持运行模式m′，则约束如下：
TTu，m，m′表示的是生产装置u从运行模式m到运行模式m′的过渡过程时长；如果Cu，m，m′，t＝1，表示生产装置u在时间间隔t-1与时间间隔t之间运行模式从m切换至m′，为使约束(6)成立，在时间间隔t+1到时间间隔t+TTm，m′的范围内应有Cu，m′，m′，t′＝1，如果Cu，m，m′，t＝0，约束(6)恒成立；
B物料平衡及容量、组分油调和、成品油交付约束
B.1质量平衡约束
B.1.1生产装置流出口质量平衡约束
如果一个生产装置具有一个以上运行模式，它的约束为：
式中，Yieldu,s,m'为生产装置u在操作模式为m'时端口s产出物料的收率；tYieldu,s,m,m'为生产装置u在操作模式从m到m'的过渡过程中端口s产出物料的收率；QIu,t为生产装置u在时间间隔t内的输入流量；QOu,s,t为生产装置u的端口s在时间间隔t内的输出流量；
S为生产装置输出端口的集合；
如果生产装置处于稳态运行过程，xu,m,m',t＝0，则∑m∑m'xu,m,m',tQIu,ttYieldu,s,m,m'＝0，因此，
QOu,s,t＝∑m'yu,m',t(1-∑mxu,m,m',t)QIu,tYieldu,s,m'
如果生产装置处于过渡过程，1-∑mxu,m,m',t＝0，则∑m'yu,m',t(1-∑mxu,m,m',t)QIu,tYieldu,s,m'＝0，因此
QOu,s,t＝∑m∑m'xu,m,m',tQIu,ttYieldu,s,m,m'
如果生产装置仅有一种生产运行模式，则约束(7)转变为：
B.1.2中间油品的质量平衡约束
中间油品来自每个生产装置的流出口；在时间间隔t内，对于中间油品oi，来自上游装置的流出量总和等于进入下游装置的输入量总和，约束表示如下：
式中，QOu,oi,t为时间间隔t内生产装置u的中间油品oi输出流量；QIu,oi,t为时间间隔t内生产装置u的中间油品oi输入流量；OI为中间 油品的集合；
B.1.3储罐质量平衡约束
每个储罐在时间间隔t结束时的储量等于在在时间间隔t-1结束时的储量加上时间间隔t内储罐的输入量并减去时间间隔t内储罐的输出量：
QOoc,t和QIo,t的关系是
其中，INVoc,t为时间间隔t结束时组分油oc的罐存量；INVoc,ini为组分油oc的初始罐存量；QIu,oc,t为时间间隔t内来自生产装置u的组分油oc输入流量；QOoc,t为时间间隔t内组分油oc输出流量；INVo,t为时间间隔t结束时成品油o的罐存量；INVo,ini为成品油o的初始罐存量；QIo,t为时间间隔t内成品油o的输入流量；Dl,0,t为时间间隔t内订单l的成品油o交付量；Qoc,o,t为时间间隔t内成品油o中的组分油oc 调和流量；OC为用于调和的组分油的集合；O为成品油的集合；L为订单的集合；
B.2容量约束
B.2.1生产装置的容量约束
该约束明确要求在时间间隔t内生产装置u的装载量必须满足容量的最小值和最大值要求；
其中， 为生产装置u的输入流量最小值； 为生产装置u的输入流量最大值；
B.2.2储罐的容量约束，储罐的库存量，包括组分油和成品油，都必须处于最小限值和最大限值之间，
式中， 为组分油oc的罐存容量最小值； 为组分油oc的罐存容量最大
min max
值；INVo 为成品油o的罐存容量最小值；INVo 为成品油o的罐存容量最大值；
B.3调和约束
B.3.1组分油调和比例约束
组分油有调和最大比例值以及调和最小比例值。相应的约束关系为：
式中， 为用于调和成品油o的组分油oc最小比例成分； 为用于调和成品油o
的组分油oc最大比例成分；
B.3.2成品油特性值约束
石油产品的主要特性值，包括汽油的研究法辛烷值(RON)和硫磺浓度值，柴油的十六烷值、硫磺浓度值和冷凝点因子值等必须处于要求的最大限值和最小限值范围内；其约束关系为：
其中， 通过将各项乘以ΣocQoc,o,t，该约束条件可等价转
化成线性表示：
为简化起见，该模型采用线性调和准则，即调和过程中的成品油主要特性值呈线性；上述式中， 为成品油o的特性p最小值；PROo,p,t为时间间隔t内成品油o的特性p的
值； 为成品油o的特性p最大值；PROoc,p为组分油oc的特性p的值；P为油品特性
的集合
B.4成品油交付约束
每个订单都有交付的起始时间和结束时间要求，成品油的交付既不能早于起始时间，也不能迟于截止时间；订单缺货会有惩罚值，调度时间结束时可计算总的缺货惩罚大小。因此成品油的供需约束要求为：
其中，Rl,o是订单l内成品油o的需求量；Tl1是订单l要求的开始交付时间间隔；Tl2是订单l要求的结束交付时间间隔。
5.根据权利要求1所述一种基于炼油厂全厂调度离散时间建模方法，其特征在于，所述步骤6模型线性化，是依据所述步骤4中的约束条件(7)和目标函数都涉及到相同的双线性项和三线性项；以上构建的调度模型(P0)中涉及有双线性项和三线性项，双线性项是一个二进制变量与一个连续变量的乘积，三线性项是两个二进制变量与一个连续变量的乘积，双线性项是xu,m,m',tQIu,t，其中，xu,m,m',t是二进制变量，QIu,t是连续变量；三线性项是yu,m',t(1-Σmxu,m,m',t)QIu,t，由xu,m,m',t的定义知，xu,m,m',t＝1表示生产装置u在时间间隔t内处于从m至m'的切换过渡过程；因为生产装置在时间间隔t-1上的操作模式唯一，故Σmxu,m,m',t的值不超过1，所以(1-Σmxu,m,m',t)可以视为二进制变量，yu,m',t是二进制变量，QIu,t是连续变量；通过引入辅助变量将这些项线性化，具体包括：
A对模型中的双线性项进行线性化
为实现线性化，引入两个辅助性连续变量xQIu,m,m',t和xQI1u,m,m',t以及如下的辅助约束条件：
上述约束条件(26)、(27)中的参数 是QIu,t最大值；约束条件(26)、(27)、(28)、(29)可确保，如果xu,m,m',t＝0，则xQIu,m,m',t＝0；如果xu,m,m',t＝1，则xQI1u,m,m',t＝0；因此，由上述约束条件可得，xQIu,m,m',t等价于xu,m,m',t和QIu,t的乘积；
B对模型中的三线性项进行线性化
B.1先引入辅助性二进制变量xyu,m',t表达yu,m',t(1-Σmxu,m,m',t)；相应的辅助约束条件如下：
上述约束条件(30)、(31)、(32)确保，如果yu,m',t＝0或1-Σmxu,m,m',t＝0，则xyu,m',t＝
0；约束条件(32)确保，如果yu,m',t＝1且1-Σmxu,m,m',t＝1，则xyu,m',t＝1；
B.2再引入两个辅助性连续变量xyQIu,m',t和xyQI1u,m',t，实现双线性项xyu,m',tQIu,t的线性化；相应的辅助约束条件如下：
(36)
约束条件(35)和(36)中的 与约束条件(26)和(27)中的 相同；
约束条件(35)、(36)、(37)、(38)确保，如果xyu,m',t＝0，则xyQIu,m',t＝0；如果xyu,m',t＝1，则xyQI1u,m',t＝0；因此，由上述约束条件可得，xyQIu,m',t等价于xyu,m',t和QIu,t的乘积。