作业车间调度风险的量化评估与控制方法
【专利摘要】本发明公开了一种作业车间调度风险的量化评估与控制方法，用于解决现有调度风险控制方法控制效果差的技术问题。技术方案是通过选择抗风险能力强的方案作为初始调度方案，并配合风险事件过滤应对控制策略。通过对生产环境中的不确定因素进行识别、分析并通过聚合运算构建描述生产环境不确定程度的综合量化模型。依据该量化模型，对调度方案进行量化评估，获得在没有风险应对活动的情况下调度的风险值，采用进化算法优化调度风险值，将获得的风险值最小的调度方案作为初始调度方案。当风险事件发生时，通过对风险事件的影响程度进行预评估，与设置的阈值进行比较来决定是否对风险事件进行应对控制，达到对风险事件的过滤，减少了对调度方案的变动。
【专利说明】作业车间调度风险的量化评估与控制方法
【技术领域】
[0001] 本发明涉及一种调度风险控制方法，特别是涉及一种作业车间调度风险的量化评 估与控制方法。
【背景技术】
[0002] 文献"Hybrid Evolutionary Algorithm for job scheduling under machine maintenance, Applied Soft Computing, 2013, Voll3 (3)，pl440_1447" 公开了一种基于混 合进化算法的调度风险控制方法。该方法在传统的进化算法的基础上，加入了一个基于左 移思想（left-shifting)的局部搜索启发，增强了算法的搜索能力。利用该算法，首先在理 想的生产环境下，以最大完工时间（makespan)为优化目标，获到最大完工时间最小的调度 方案作为初始调度方案。当机器故障发生时，再次利用该混合进化算法对受影响的工序进 行再调度，以使得再调度后的方案的最大完工时间与初始的最大完工时间偏差最小。该方 法采用反应式调度的方法较好的解决了风险应对问题，可以降低调度的风险水平。然而，由 于文献所述方法以具有最小最大完工时间的调度方案为初始方案,使得机器上的空闲时间 被高度压缩，调度方案抵御风险事件的能力降低。当机器故障发生时，大量的工序将会受到 机器故障的影响，加上频繁的再调度活动，初始调度方案将发生巨大的变化，从而使调度稳 定性降低，甚至使生产过程陷入混乱。再者，生产过程中有很多种风险事件，而文献只考虑 了机器故障这一种风险事件，因此在多种风险事件作用下，该方法的控制效果不能被保证。
【发明内容】
[0003] 为了克服现有调度风险控制方法控制效果差的不足，本发明提供一种作业车间调 度风险的量化评估与控制方法。该方法通过选择抗风险能力强的方案作为初始调度方案， 并配合风险事件过滤应对控制策略，达到在控制调度风险的同时，保证生产过程的稳定性。 通过对生产环境中的不确定因素进行识别、分析并通过聚合运算构建描述生产环境不确定 程度的综合量化模型。依据该量化模型，对调度方案进行量化评估，获得在没有风险应对活 动的情况下调度的风险值，进而采用进化算法优化调度风险值，将获得的风险值最小的调 度方案作为初始调度方案。这样的初始调度方案具有很强的适应能力，可以减小风险事件 的影响。当风险事件发生时，通过对风险事件的影响程度进行预评估，进而与设置的阈值进 行比较来决定是否对风险事件进行应对控制，达到对风险事件的过滤，减少了对调度方案 的变动。此外，多种风险事件同时被考虑，提高了该调度风险控制方法的适应能力。
[0004] 本发明解决其技术问题所采用的技术方案是：一种作业车间调度风险的量化评估 与控制方法，其特点是采用以下步骤：
[0005] 步骤一、（a)搜集车间生产故障报告表和车间生产管理信息系统的生产事件记录， 提取表中和记录中影响调度性能的风险事件，并记录风险事件的类型、发生时刻、影响时 长；(b)以事件ID、事件名、事件类型、发生时刻、影响时长为字段建立数据表，将记录的风 险事件信息数据逐个录入数据表中，构成风险信息数据库。
[0006] 步骤二、（a)按风险事件的类型对风险事件数据库中的风险事件进行分类，每一种 类型为一个风险因素；
[0007] (b)针对各个风险因素，依据数据库中的风险事件基本数据，以月为统计周期对各 个风险因素进行统计分析，确定各个风险因素的统计信息完备程度；
[0008] (1)风险因素的发生概率模型。
[0009] 利用风险因素在各统计周期内的统计数据，采用统计分析方法确定风险因素发生 的概率模型。通过公式（1)求得风险因素的发生概率模型。
【权利要求】
1. 一种作业车间调度风险的量化评估与控制方法，其特征在于包括以下步骤： 步骤一、（a)搜集车间生产故障报告表和车间生产管理信息系统的生产事件记录，提取 表中和记录中影响调度性能的风险事件，并记录风险事件的类型、发生时刻、影响时长；(b) 以事件ID、事件名、事件类型、发生时刻、影响时长为字段建立数据表，将记录的风险事件信 息数据逐个录入数据表中，构成风险信息数据库； 步骤二、（a)按风险事件的类型对风险事件数据库中的风险事件进行分类，每一种类型 为一个风险因素； (b)针对各个风险因素，依据数据库中的风险事件基本数据，以月为统计周期对各个风 险因素进行统计分析，确定各个风险因素的统计信息完备程度； (1) 风险因素的发生概率模型； 利用风险因素在各统计周期内的统计数据，采用统计分析方法确定风险因素发生的概 率模型；通过公式（1)求得风险因素的发生概率模型；
式中，P代表风险因素的发生概率；n代表统计周期的个数山代表风险事件在第i个 统计周期内的总持续时间；代表在第i个统计周期内正常工作的时间； 采用3-〇表征法，按照公式（2)和公式（3)分别确定m、〇的值，进而构建该风险因素 的模糊概率模型：m-3 〇，m，m+3 〇 ;
式中，n为样本总数；Xi为样本概率值；
式中，E(Xi)为样本概率均值； (2) 风险因素的程度分布模型； 通过参数估计的方法获得风险因素程度概率分布模型；机器故障的持续时间服从指数 分布，如式（4)所示，参数A未知；
式中，d为故障持续时间八为指数分布的参数； 采用极大似然估计法对式（4)中的参数A进行估计获得风险因素的程度概率分布模 型；具体步骤如下： 写出似然函数，如式（5):
式中，屯为第i次机器故障持续时间； 求出InLU)及似然方程，式（6):
解该似然方程即得到参数A的极大似然估计，式（7):
将历史统计数据值带入式（7)求得分布函数（4)的参数估计值，从而确定机器故障的 程度概率分布模型； 每个统计区间就代表一个程度等级，区间越小，则描述越精细；用式（8)计算各个程度 等级的出现概率；
式中，Pi代表风险因素程度为i时的概率是风险因素程度为i时的出现频数；m为 总的程度等级数； 采用信息扩散的方法对统计数据进行处理，具体步骤如下： 确定论域，式（9): U = {屮，u2, ? ? ?屮，? ? ?，uj (9) 式中，W代表第i个论域，即故障程度的第i个等级；n为论域的个数； 通过式（10)将单一观测样本携带的信息按正态分布分配给U中的所有论域；
式中，h为扩散系数，由公式（11)确定：
b=max(yi}j= 1, 2,. . . ,m),a=min(y^j= 1, 2,. . . ,m) 式中，b为样本最大值；a为样本最小值，m为样本总数； 由式（12)进行信息规范化处理：
式中，&由公式（13)确定：
各论域的概率由公式（14)算得：
式中，q(Ui)，Q分别由公式（15)和公式（16)计算；
(c) 对于统计数据完备度高的风险因素采用统计推理的方法获得风险因素的发生概率 模型和程度概率分布模型，对于统计数据不足但分布形式已知的风险因素采用信息扩散方 法获得风险因素的发生概率模型和程度概率分布模型，对于统计数据不足且分布形式未知 的风险因素，则采用模糊处理操作来构建风险因素的发生概率模型和程度概率分布模型； 将明确的风险因素发生概率模型Pi模糊化，使得所有风险因素的发生概率模型都为模 糊概率，如式（17):
式中，^为风险因素发生概率Pi的模糊概率； 将所有风险因素的发生概率模型做模糊求和，获得调度风险的综合发生模糊概率模 型，如式（18):
由公式（19)对该综合模糊概率发生模型去模糊化处理，即获得调度风险综合发生概 率模型P;
通过仿真方法综合区间离散型、连续型和模糊型概率分布模型，构建一个描述调度风 险综合程度分布概率模型；依据各个风险因素的程度概率分布模型进行仿真取样，进而对 样本数据进行统计分析并通过曲线拟合技术，获得调度风险综合程度概率分布模型f(d); (d) 采用风险聚合推理技术，建立综合发生概率模型和综合程度概率分布模型；结合 调度方案与调度风险综合描述模型获得调度方案风险评估模型，如式（20)所示；通过模拟 仿真的方法近似求解该调度风险评价模型； R=EP/g(d)f(d)dd,d>0 (20) 式中，g(d)为风险因素程度为d时对调度最大完工时间的影响值； (e) 基于此风险因素综合模型，利用模拟仿真技术构建调度风险量化评估模型来对调 度风险大小进行量化评估； 步骤三、（a)调度方案的风险值评估：按照调度风险的综合发生概率P及其综合程度概 率分布f(d)模拟产生风险事件，以调度最大完工时间的平均影响值作为调度风险值r，如 式（21);
式中，m代表模拟仿真的次数；n代表一次仿真中，发生的风险事件的个数；gWi)代表 事件程度为屯时对调度最大完工时间的影响值； (b)调度风险值优化：设置优化目标为最小化f(t) =MK+r，其中MK为调度方案的最大 完工时间；以该优化目标引导进化算法的搜索过程，优化并输出具有良好抗风险能力的调 度方案； 步骤四、（a)风险事件的应对决策：风险事件应对决策，负责对调度过程中的遇到的风 险事件进行过滤控制；根据风险事件对调度性能影响大小的预测评估，通过与决策机制的 阈值进行比较，来确定是否对风险事件进行应对处理； 风险事件对调度性能的影响程度评估：当风险事件发生后，根据该风险事件的程度概 率分布模型，通过模拟仿真和右移修复策略可以快速获得该风险事件对调度性能的影响程 度Et，通过式（22)计算； Et =g(dt)/MKold = (MKnew-MKold)/MKold (22) 式中，MKnOT为经受风险事件影响后的最大完工时间；MK^是在风险事件发生之前的调 度最大完工时间； 风险应对决策：设置风险应对决策不等式为式（23)，取风险应对决策的阈值为0. 001 ; 如式（23)所示，是否对风险事件做出应对，由风险事件发生的时刻和风险事件的影响程度 共同决定；
式中，△为风险应对决策的阈值；t为风险事件发生的时刻； (b)风险应对措施来源：1)生产调度研究领域中采用的应对不确定事件影响的控制措 施，如局部重调度方法、完全重调度方法；2)生产实际中采用的应对不确定事件的控制方 法，如采用并行机、提前加工； 风险应对措施的选择：通过模拟仿真，选择在该风险事件作用下，Et值最小的控制措 施，确定风险应对措施与该风险因素的对应关系；当风险事件发生时，由确定的对应关系选 择出相应的风险应对控制措施； 风险应对控制：采用选择的风险应对控制措施，对调度方案、环境做出调整，减少风险 事件对调度最大完工时间的影响； 步骤五、风险信息采集更新； (a) 信息采集更新：随着生产的进行，风险事件会不断发生，关于风险事件的信息会以 生产事故报告表的形式或被生产管理信息系统记录下来；将这些新的风险事件信息搜集起 来、并按照风险信息数据表所需格式进行数据处理；将处理后的数据信息，按照风险信息数 据库数据表的格式，逐个录入到风险信息数据库中； (b) 风险分析模型的更新：随着风险信息数据库的不断更新、完善，就需要重新进行风 险分析过程更新已有的风险分析模型来提高对风险因素的描述精度，保证风险综合评估模 型的可靠性；根据企业实际可采用定期或自适应触发的方式来选择触发风险分析过程的时 间。
【文档编号】G06Q10/04GK104408525SQ201410631041
【公开日】2015年3月11日 申请日期:2014年11月11日 优先权日:2014年11月11日
【发明者】孙树栋, 吴自高, 肖世昌, 国欢, 徐元超, 秦建伟, 杨宏安, 王军强 申请人:西北工业大学