一种制药化工生产中原料结晶过程的优化调度方法
【专利摘要】本发明涉及一种制药化工生产中原料结晶过程的优化调度方法，属于化工生产过程智能优化调度【技术领域】。本发明通过确定原料结晶过程调度模型和优化目标，并使用基于贝叶斯网络的分布估计算法的优化调度方法对优化目标进行优化；其中调度模型依据每台结晶机器上所加工原料的最大完工时间来建立，优化目标为最小化最大完工时间。本发明使得化工生产中原料结晶过程的表达清晰准确，调度方法合理有效；使得算法的寻优能力得到增强。
【专利说明】-种制药化工生产中原料结晶过程的优化调度方法
【技术领域】
[0001] 本发明设及一种制药化工生产中原料结晶过程的优化调度方法，属于化工生产过 程智能优化调度【技术领域】。
【背景技术】
[0002] 随着全球经济的持续快速增长，流程工业的发展程度已成为衡量国家工业水平的 重要指标。化工生产是流程工业最重要的组成部分之一，化工生产通过对原材料进行相应 的物理和化学加工，进而实现原材料价值的升值。在化工生产中，主要包括过滤、萃取、结晶 和蒸发等基本操作，其中结晶操作应用的场合相当广。结晶过程常用于原料的初步处理和 最终产品的输出，原料的处理的速度和质量将直接对后续的深加工产生重要影响，对整个 系统结构的优化升级具有重要意义。因此，化工生产中原料结晶过程的优化调度具有重要 研究价值。
[0003] 对原材料进行结晶是常见的化工操作，由于原材料的理化性质存在差异，有些原 材料经过一次结晶即可完成加工，而有些原材料需要进行二次或多次结晶操作才能完成加 工；同时，为了防止产品间的交叉污染，每台机器在加工完一种原材料后，需要一定的设置 时间进行清洗和调整，才能继续加工另外一种原材料，而设置时间取决于原材料之间的加 工顺序；此外，实际结晶操作的加工（机器）单元往往由多台异构的并行机器构成，该些机 器的加工能力不尽相同，原材料需根据自身的理化性质、体积、质量等因素选择合适的机器 进行加工。该过程就是典型的异构并行机调度问题。与其他生产调度问题一样，异构并行 机调度问题的优化指标主要包括最大完工时间、拖期产品数、平均流经时间等，其中W最大 完工时间（makespan或CmJ使用最为广泛。异构并行机调度问题属于NP-Complete问题， 其解空间随问题规模的增大呈指数增长。因此，对异构并行机调度问题求解算法的研究具 有较高的实际和理论价值，可为相关化工生产优化系统的设计提供切实指导。
[0004] 由于原料结晶过程调度问题属NP完全范畴，传统的数学规划方法和启发式构造 性方法无法保证解的优化质量。因此，本发明设计了一种基于贝叶斯网络的分布估计算 法化ayesian statistical inference-based estimation of distribution algorithm ，BEDA)的优化调度方法对优化目标进行优化，可在较短时间内获得制药化工生产中原料结 晶过程的优化调度方法的近似最优解。
【发明内容】
[0005] 本发明所要解决的技术问题是在较短时间内获得制药化工生产中原料结晶过程 调度问题的近似最优解的问题，提供了一种制药化工生产中原料结晶过程的优化调度方 法。
[0006] 本发明的技术方案是：一种制药化工生产中原料结晶过程的优化调度方法，通过 确定原料结晶过程调度模型和优化目标，并使用基于贝叶斯网络的分布估计算法的优化调 度方法对优化目标进行优化；其中调度模型依据每台结晶机器上所加工原料的最大完工时 间来建立，优化目标为最小化最大完工时间：
[0007] Cmax= max(C 1= F i+ti) i G D
[000引 ZieDXuk=l j G D，（j，k|k = Mi=Mj)
[0009] ZjED Xyk= 1 i G D, (i, k|k = M i= M j)
[0010] Fi+ti+Su《 Fj'+B(l-Xuk) (i，j) G D，（i，j，k|k = Mi= Mj)
[0011] Fi+ti《Fj i, j G Sij
[0012] /&,'=〇 k G M
[0013] S馬,' =O'S,'馬=0 i e D, Q, k I k = Mi)
[0014] 其中，B表示一个正数且大于模型中的任意一个变量值；D表示所有有效的操作序 列；Mj.表示操作j对应的加工机器；X uk表示是在机器k上的操作j正好在操作i后边时 xyk= 1，否则等于0 ;t J表示操作j对应的加工时间；F J表示操作j的开始时间；C J表示操 作j的完成时间化表示在机器k上的虚拟操作；S U表示操作i在操作j后加工时机器上 的设置时间；n为待结晶的原料排序；操作表示的是第1种原料的第m个结晶阶段在机器 k上的加工；调度的目标为在所有待结晶的原料排序的集合n中找到一个JT%使得最大完 工时间最小。
[0015] 所述基于贝叶斯网络的分布估计算法的优化调度方法的具体步骤如下：
[0016] St巧1、编码方式；W待加工原料排序进行编码n = [JT。JI2，...，JTts];其中TS 表示所有产品的总工序数；
[0017] Step2、种群的初始化：利用步骤Stepl中的编码方式，采用随机方法产生初始化 种群，直至初始解的数量达到种群规模的要求，同时选择当代中种群中的"最优个体"作为 第一代的"历史最优个体"；其中，gen = 1 ;
[0018] St巧3、邻域生成；利用Insed的方法生成当前种群中"最优个体"的邻域；
[0019] Step4、条件概率矩阵的更新机制；将精英个体描绘在贝叶斯网络上，每出现一次 节点 Nu，p (曰 G 1，...，TS，0 G l，...，n)到节点 （曰 G 1，...，TS，0 ' G l，...，n) 的情况，便在该有向弧的权重上加一；接着根据之前选择的精英个体构造的网络，生成条件 概率矩阵；其中用于构建条件概率矩阵的精英个体是由从该代种群及"最优个体"的邻域中 适配值最好的前e%个优质个体构成，n表示原料数；
[0020] steps、产生新种群；根据轮盘赌的方法对条件概率矩阵采样形成新种群 POPeaWi(gen)中的个体；
[0021] steps、更新"历史最优个体";对比"历史最优个体"与新种群中的"最优个体"的 适配值，选择适配值较小的个体更新"历史最优个体"；
[002引 St巧7、终止条件；设定终止条件的最大迭代次数为200,如果满足，则输出"历史 最优个体"；否则转至步骤Step3,反复迭代，直到满足终止条件。
[0023] 所述种群规模设置为50, e%= 0. 6,交叉概率为0. 7。
[0024] 本发明的工作原理是；
[0025] 步骤1 ;建立结晶过程调度模型和优化目标。
[0026] 调度模型依据每台结晶机器上所加工原料的最大完工时间来建立，优化目标为最 小化最大完工时间Cmax;
[0027] Cmax= max (C i= F i+ti) i G D
[002引 ZieDXuk=l j G D, (j, k|k = M i= M J)
[0029] Zj'eDXuk=l i G D，（i，k|k = Mi= Mj)
[0030] Fi+ti+Su《F j'+B(l-Xuk) (i，_j) G D，（i，j，k|k = Mi= Mj)
[0031] Fj+tj^ F j i, j G S。'
[0032] hk=Q k G M
[003引 = 0,％/、. = 0 i e D, (i, k I k = Mi)
[0034] 式中；B表示一个很大的正数；D表示所有有效的操作序列；Mj.表示操作j对应的 加工机器；xuk表示是在机器k上的操作j正好在操作i后边时X uk= 1，否则等于0 ;t J表 示操作j对应的加工时间；。表示操作j的开始时间；Cj.表示操作j的完成时间；Rk表示在 机器k上的虚拟操作；Sy表示操作i在操作j后加工时机器上的设置时间；操作表示的是 第1种原料的第m个结晶阶段在机器k上的加工。设JT为待结晶的原料排序；调度的目标 为在所有待结晶的原料排序的集合n中找到一个JT %使得最大完工时间最小。
[0035] 步骤2;解的表达。
[0036] 编码是指建立反映结晶操作解空间结构和特征的编码排列，智能算法通过操作编 码排列实现迭代寻优。解码是指通过利用特定的单元选择规则，结合编码排列合理安排每 台机器上的原料排序，进而形成可行的调度方案。对于调度问题，经典的编码方式主要包括 基于机器的编码，基于工件的编码和基于随机键的编码等。针对化工生产中原料结晶机器 的特征，本发明提出基于原料总工序数的排序编码。
[0037] 譬如，对于n = 4的一个问题（对于n个待加工的原料数），则对于一个解31 = [2, 4, 3, 1，4, 3, 4, 4, 1，3]就是问题的一个编码排列，同时也是问题的一个可行解，该种编码 方式建立了编码排列到问题可行解的一一映射关系，简单直观，且容易实现。
[003引步骤3;种群的初始化。
[0039] 令gen = 1，利用步骤2中的编码方式，采用随机方法产生初始化种群，直至初始解 的数量达到种群规模的要求。选择当代中种群中的"最优个体"作为第一代的"历史最优个 体"。
[0040] 步骤4 ;邻域生成：
[0041] 利用Insed的方法生成当前种群中"最优个体"的邻域。
[0042] 步骤5 ;条件概率矩阵的更新机制。
[0043] 将精英个体描绘在贝叶斯网络上，每出现一次节点 Na，p (a G 1，...，TS，0 G l，...，n)到节点 （a G 1，...，TS，0 ' G l，...，n)的 情况，便在该有向弧的权重上加一；接着根据之前选择的精英个体构造的网络，生成条件概 率矩阵；其中用于构建条件概率矩阵的精英个体是由从该代种群及"最优个体"的邻域中适 配值最好的前e%个优质个体构成，n表示原料数；
[0044] 步骤6;产生新种群。
[0045] 新种群pop。。。,; (gen)中的个体根据轮盘赌的方法对条件概率矩阵采样形成。
[0046] 步骤7 ;更新历史最优个体。
[0047] 对比"历史最优个体"与新种群中的"最优个体"的适配值，更新"历史最优个体"。
[0048] 步骤8 ;设定终止条件的最大迭代次数为200。如达到设定的最大迭代次数200， 则输出"历史最优个体"；否则，令gen = gen+1，返回步骤4。
[0049] 本发明的有益效果是：
[0050] 1、提出了最小化最大完工指标下的制药化工生产中原料结晶过程的优化调度方 法，使得化工生产中原料结晶过程的表达清晰准确，调度方法合理有效；
[0化^ 2、通过贝叶斯与邸A的结合，使得邸DA具有全局与局部捜索能力；使得邸DA可W 有效克服在局部捜索上的不足；有利于充分利用优秀个体的信息来指导捜索方向，进而使 得算法的捜索宽度和深度得到合理平衡，使得算法的寻优能力得到增强。
[0化引 3、提出的基于邸DA的优化调度方法可W有效解决制药化工生产中原料结晶过程 调度问题。
【专利附图】
【附图说明】
[0053] 图1为本发明的整体流程图；
[0054] 图2为本发明中"Insert"操作示意图；
[0化5] 图3为本发明中贝叶斯网络示意图。
【具体实施方式】
[0化6] 实施例1 ;如图1-3所示，一种制药化工生产中原料结晶过程的优化调度方法，通 过确定原料结晶过程调度模型和优化目标，并使用基于贝叶斯网络的分布估计算法的优化 调度方法对优化目标进行优化；其中调度模型依据每台结晶机器上所加工原料的最大完工 时间来建立，优化目标为最小化最大完工时间：
[0057] Cmax= max(C 1= F i+ti) i G D
[005引 ZieDXuk=l j G D, (j, k|k = M i= M j)
[0059] ZjED Xyk= 1 i G D, (i, k|k = M i= M j)
[0060] Fi+ti+Su《F J+B(l-Xuk) (i，j) G D，（i，j，k|k = Mi= Mj)
[0061] Fi+ti《Fj i, j G Sij
[0062] k G M
[0063] =化'、?心二 0 i G D, (i, k I k = Mi)
[0064] 其中，B表示一个正数且大于模型中的任意一个变量值；D表示所有有效的操作序 列；Mj.表示操作j对应的加工机器；X uk表示是在机器k上的操作j正好在操作i后边时 Xyk= 1，否则等于0 j表示操作j对应的加工时间；F j表示操作j的开始时间；C j表示操 作j的完成时间化表示在机器k上的虚拟操作；S U表示操作i在操作j后加工时机器上 的设置时间；n为待结晶的原料排序；操作表示的是第1种原料的第m个结晶阶段在机器 k上的加工；调度的目标为在所有待结晶的原料排序的集合n中找到一个31%使得最大完 工时间最小。
[00化]所述基于贝叶斯网络的分布估计算法的优化调度方法的具体步骤如下：
[0066] St巧1、编码方式；W待加工原料排序进行编码n = [JTi，312，...，JTts];其中TS 表不所有广品的总工序数；
[0067] Step2、种群的初始化：利用步骤Stepl中的编码方式，采用随机方法产生初始化 种群，直至初始解的数量达到种群规模的要求，同时选择当代中种群中的"最优个体"作为 第一代的"历史最优个体"；其中，gen = 1 ;
[0068] St巧3、邻域生成；利用Insed的方法生成当前种群中"最优个体"的邻域；
[0069] Step4、条件概率矩阵的更新机制；将精英个体描绘在贝叶斯网络上，每出现一次 节点 Nu，p (曰 G 1，...，TS，0 G l，...，n)到节点 （曰 G 1，...，TS，0 ' G l，...，n) 的情况，便在该有向弧的权重上加一；接着根据之前选择的精英个体构造的网络，生成条件 概率矩阵；其中用于构建条件概率矩阵的精英个体是由从该代种群及"最优个体"的邻域中 适配值最好的前e%个优质个体构成，n表示原料数；
[0070] steps、产生新种群；根据轮盘赌的方法对条件概率矩阵采样形成新种群 POPeaWi(gen)中的个体；
[0071] steps、更新"历史最优个体";对比"历史最优个体"与新种群中的"最优个体"的 适配值，选择适配值较小的个体更新"历史最优个体"；
[007引 St巧7、终止条件；设定终止条件的最大迭代次数为200,如果满足，则输出"历史 最优个体"；否则转至步骤Step3,反复迭代，直到满足终止条件。
[0073] 实施例2 ;如图1-3所示，一种制药化工生产中原料结晶过程的优化调度方法，通 过确定原料结晶过程调度模型和优化目标，并使用基于贝叶斯网络的分布估计算法的优化 调度方法对优化目标进行优化；其中调度模型依据每台结晶机器上所加工原料的最大完工 时间来建立，优化目标为最小化最大完工时间：
[0074] Cmax= max(C 1= F i+ti) i G D
[007引 ZieDXuk=l j G D, (j, k|k = M i= M j)
[0076] ZjED Xyk= 1 i G D, (i, k|k = M i= M j)
[0077] Fi+ti+Su《F j'+B(l-Xuk) (i，_]?) G D，（i，j，k | k = M;= M j)
[007引 Fi+ti《Fj i, j G Sij
[0079] 二、、 k G M
[0080] -S'/如' =化 &鮮=0 i e D, Q, k I k = Mi)
[0081] 其中，B表示一个正数且大于模型中的任意一个变量值；D表示所有有效的操作序 列；Mj.表示操作j对应的加工机器；X uk表示是在机器k上的操作j正好在操作i后边时 xyk= 1，否则等于0 J表示操作j对应的加工时间；F J表示操作j的开始时间；C J表示操 作j的完成时间化表示在机器k上的虚拟操作；S U表示操作i在操作j后加工时机器上 的设置时间；n为待结晶的原料排序；操作表示的是第1种原料的第m个结晶阶段在机器 k上的加工；调度的目标为在所有待结晶的原料排序的集合n中找到一个JT%使得最大完 工时间最小。
[0082] 所述基于贝叶斯网络的分布估计算法的优化调度方法的具体步骤如下：
[0083] St巧1、编码方式；W待加工原料排序进行编码n = [JT。JI2，...，jTts];其中TS 表示所有产品的总工序数；
[0084] Step2、种群的初始化：利用步骤Stepl中的编码方式，采用随机方法产生初始化 种群，直至初始解的数量达到种群规模的要求，同时选择当代中种群中的"最优个体"作为 第一代的"历史最优个体"；其中，gen = 1 ;
[0085] St巧3、邻域生成；利用Insed的方法生成当前种群中"最优个体"的邻域；
[0086] Step4、条件概率矩阵的更新机制；将精英个体描绘在贝叶斯网络上，每出现一次 节点 Nu，p (曰 G 1，...，TS，0 G l，...，n)到节点 （曰 G 1，...，TS，0 ' G l，...，n) 的情况，便在该有向弧的权重上加一；接着根据之前选择的精英个体构造的网络，生成条件 概率矩阵；其中用于构建条件概率矩阵的精英个体是由从该代种群及"最优个体"的邻域中 适配值最好的前e%个优质个体构成，n表示原料数；
[0087] steps、产生新种群；根据轮盘赌的方法对条件概率矩阵采样形成新种群 POPeaWi(gen)中的个体；
[008引 steps、更新"历史最优个体";对比"历史最优个体"与新种群中的"最优个体"的 适配值，选择适配值较小的个体更新"历史最优个体"；
[0089] St巧7、终止条件；设定终止条件的最大迭代次数为200,如果满足，则输出"历史 最优个体"；否则转至步骤Step3,反复迭代，直到满足终止条件。
[0090] 所述种群规模设置为50, e% = 0. 6,交叉概率为0. 7, B取值为10000，
[0091] 具体的对比实验如下：
[0092] 采用国内某制药化工企业生产中结晶过程的实际生产数据，进一步对BEDA(如图 1)的性能进行了验证。该企业在生产中需要通过结晶的方式对制药原料进行提纯处理，现 有5种组份不同的原料（即纯度、杂质类型等不同），4台不同类型的结晶机器（即容量、结 晶方式等不同），用于生产5个批次的药品。在对5种原料进行结晶操作时，需考虑如下因 素；（1)、5种原料的纯度和杂质类型不同，使得每种原料所需要的结晶提纯次数不同；(2)、 不同原料的体积不同，故需要根据机器容量选取机器；(3)、在结晶过程中，原料产品的纯度 不断提高，使得同种原料在不同结晶阶段的溶剂选择、结晶温度、揽拌方式都有所不同，所 W对于同种产品的不同结晶阶段所需的机器不同；(4)、为了防止原料批次间污染，每台机 器在加工完一种原材料后，需要一定的设置时间进行清洗和调整，才能继续加工另外一种 原材料；(5)、原料的加工时间tim化）（表示第1种原料（表中"产品"）的第m个结晶阶段 (表中"阶段"）在机器k(表中"M1/M2/M3/M4")上的加工时间）取决于原料所处的结晶阶 段和所选择的机器，因此每种原料在各个结晶阶段所需的加工时间不同。5种原料的工序数 集合31 = [1，2, 5, 3, 4, 5, 2, 3, 5, 4, 4, 3, 5]，原料的加工约束如表1所示；
[0093] 表1实例的加工约束
[0094]
【权利要求】
1. 一种制药化工生产中原料结晶过程的优化调度方法，其特征在于：通过确定原料结 晶过程调度模型和优化目标，并使用基于贝叶斯网络的分布估计算法的优化调度方法对优 化目标进行优化；其中调度模型依据每台结晶机器上所加工原料的最大完工时间来建立， 优化目标为最小化最大完工时间：
其中，B表示一个正数且大于模型中的任意一个变量值；D表示所有有效的操作序列；Mj表示操作j对应的加工机器；Xijk表示是在机器k上的操作j正好在操作i后边时Xijk = 1，否则等于〇 4表示操作j对应的加工时间；F^表示操作j的开始时间；q表示操作j的 完成时间；Rk表示在机器k上的虚拟操作；Su表示操作i在操作j后加工时机器上的设置 时间；为待结晶的原料排序；操作表示的是第1种原料的第m个结晶阶段在机器k上的 加工；调度的目标为在所有待结晶的原料排序的集合n中找到一个31%使得最大完工时 间最小。
2. 根据权利要求1所述的制药化工生产中原料结晶过程的优化调度方法，其特征在 于：所述基于贝叶斯网络的分布估计算法的优化调度方法的具体步骤如下： St印1、编码方式：以待加工原料排序进行编码JT= [JT1,JT2,...，JTts];其中TS表示 所有广品的总工序数； Step2、种群的初始化：利用步骤Stepl中的编码方式，采用随机方法产生初始化种群， 直至初始解的数量达到种群规模的要求，同时选择当代中种群中的"最优个体"作为第一代 的"历史最优个体"；其中，gen= 1 ; Step3、邻域生成：利用Insert的方法生成当前种群中"最优个体"的邻域； Step4、条件概率矩阵的更新机制：将精英个体描绘在贝叶斯网络上，每出现一次节点Na，p(aG1，...，TS，0Gl，...，n)到节点Na+1，p,(aG1，...，TS，0，Gl，...，n)的 情况，便在该有向弧的权重上加一；接着根据之前选择的精英个体构造的网络，生成条件概 率矩阵；其中用于构建条件概率矩阵的精英个体是由从该代种群及"最优个体"的邻域中适 配值最好的前e%个优质个体构成，n表示原料数； Step5、产生新种群：根据轮盘赌的方法对条件概率矩阵采样形成新种群pop_di (gen) 中的个体； Step6、更新"历史最优个体":对比"历史最优个体"与新种群中的"最优个体"的适配 值，选择适配值较小的个体更新"历史最优个体"； Step7、终止条件：设定终止条件的最大迭代次数为200,如果满足，则输出"历史最优 个体"；否则转至步骤St印3,反复迭代，直到满足终止条件。
3. 根据专利要求1或2所述的制药化工生产中原料结晶过程的优化调度方法，其特征 在于：所述种群规模设置为50,e%= 0. 6,交叉概率为0. 7。
【文档编号】G05B19/418GK104503382SQ201410668205
【公开日】2015年4月8日 申请日期:2014年11月20日 优先权日:2014年11月20日
【发明者】钱斌, 陈少峰, 胡蓉, 张桂莲 申请人:昆明理工大学