修复状态不确定的巡回修理任务多目标动态调度
刘 彦1， 陈春良1， 陈伟龙2， 范江波3
(1. 陆军装甲兵学院装备保障与再制造系， 北京100072； 2. 中国卫星海上测控部， 江苏 江阴 214431， 3. 66058部队， 天津 301700)
摘 要： 针对战时装备巡回修理任务重、修理时间有限，而现有战时装备维修任务调度研究存在调度目标单一、约束简单的问题，开展了装备巡回修理任务多目标动态调度研究。提出了复杂约束条件下的装备巡回修理任务调度军事问题，考虑待修装备修复状态的不确定性，在修理能力和修理时间限制的基础上引入了修理时间窗、非遍历、巡回修理组修理能力变化等约束，构建了以修竣装备总数、修竣装备重要度总和、获得的二次作战总时间最大为调度目标的多目标动态调度模型。设计了基于改进非支配排序遗传算法(Non-dominated Sorting Genetic Algorithm II，NSGA-II)进行模型求解，并通过示例仿真与分析，验证了模型和算法的合理性及有效性。
关键词： 巡回修理；多目标动态调度；修复状态不确定性；修理时间窗
现代战争作战节奏越来越快，作战进程不断缩短，作战空间逐渐变广，使得战损装备多、战时装备保障时间缩短，维修任务繁重，对战时装备维修的时效性提出了更高的要求。巡回修理作为战时装备维修的主要修理方式之一，通过派遣多个巡回修理组对分布在战场的众多待修装备实施抢救修理，实现作战部队战斗力的保持与恢复。
面对众多随机出现的不同损伤程度的待修装备，如何在有限的修理时间内，综合权衡待修装备位置、修理工作量、重要度、各巡回修理组的修理能力及其变化等复杂因素，合理安排巡回修理任务，解决“先修谁”“由谁修”“何时修”“怎么修”的战时装备维修任务调度问题，从而实现维修效益全局最优化，是具有重大意义且亟待解决的军事难题。
由于战场态势复杂、约束条件多、研究难度大，众多学者对装备维修任务调度进行了深入探索，取得了一定成果，具有一定的指导意义。如：张立民等[1]分析了战时装备维修任务调度的特点，建立了基于最大保障时间的维修任务静态调度模型；王正元等[2]以尽快恢复装备战斗力为目标，探究了动态维修任务调度的优化方法；陈春良等[3]考虑时间约束及负载能力约束的影响，将装备维修任务调度问题转化为车辆路径问题，并设计了基于改进最大最小蚂蚁算法(Max-Min Ant System,MMAS)的维修任务规划方法，实现了单修理组的维修任务调度；陈伟龙等[4]提出了进攻作战抢修任务动态调度问题，将其抽象为动态车辆路径问题，以获得的二次作战时间最大为调度目标构造了调度模型，并设计了变体遗传算法(Genetic Algorithm,GA)进行模型求解。
从模型抽象的角度来看，装备维修任务调度问题主要可抽象为车间作业调度问题[1]、资源约束项目调度问题[5]以及车辆路径问题[4，6]，其中前二者适用于基地级修理及定点修理，后者更适用于伴随修理及巡回修理。随着装备维修任务调度研究的不断深入，调度目标根据任务需求(如维修效益[7]、装备重要度[8]、二次作战时间[9]等)不断丰富，但大多数研究仅限于单目标调度；另一方面，约束条件根据战场情况不断细化，但主要聚焦在传统约束，如考虑装备维修流程[10]、维修专业[11]、维修力量休息[5]、不确定性[12-13]等。陈春良等[9]考虑战场抢修时修理能力的限制，提出了非遍历型调度，符合战时修理实际，具有一定指导意义；陈伟龙等[12]虽然考虑了待修装备修复状态，但仅仅停留在对完成任务数的影响，并未考虑修复状态对维修时间及装备重要度的影响。
将设计图纸变成现实是室内设计最直接的体现，此过程需要相应的施工单位来完成。通常在设计方案完善后，会有很多承包单位进行竞标，发标方会综合考虑竞标单位的实际情况，然后按照自身标准选出实力最强的单位，委托其完成商场建造[2]。在后期施工中，还要做好现场监管、质量评估等多项工作。
总体来看，现阶段装备维修任务调度研究尚未针对修理方式的差异进行研究，也未对巡回修理中如何合理规划诸多实际约束下维修任务以及动态处理不确定性信息问题进行深入研究。现有研究大多数是以最大保障时间、最长二次作战时间等作为单一调度目标进行任务调度，难以全面考虑维修任务需求，且约束条件相对较为理想化，忽略了修复状态的不确定性及其对维修时间和装备重要度的影响，也未考虑时间窗约束对维修任务调度的影响。
周边的各地区都想发展当地经济，但许多项目扎堆建设，缺乏管理规划。洞庭湖地区虽然各种资源丰富，文化底蕴雄厚宜打造文化之旅，旅游相关配套设施正在建设，但招牌口碑未形成集群，对游客吸引力欠缺。同时当地居民缺乏科学保护观念，开发同时造成了生态的破坏，加之受周围河湖气候变化和人类活动等的影响，低枯水位期提前，年入湖水量大幅减少，萎缩日益明显；化学肥料排放增加，“三农”问题突出，湖体自净能力降低，水质呈中营养状态，居民饮水安全受到严重威胁；同时基础设施滞后，地方财力薄弱，镇化水平不高，使湖区变成相对落后地区。
针对以上问题，笔者以战时装备巡回修理任务调度问题为研究对象，考虑待修装备修复状态的不确定性及其对维修时间和装备重要度的影响，引入待修装备修理时间窗约束及非遍历约束，考虑待修装备的不确定性，构建了多个装备巡回修理组任务动态调度模型，设计了改进的NSGA-Ⅱ(Non-dominated Sorting Genetic Algorithm Ⅱ)算法对模型进行求解，并通过示例验证了模型及算法的有效性及合理性。
由表3和表4可知，各因素对提取辣椒红色素含量影响的主次顺序为：粒度>酶解时间>SSL的添加量>料液比>酶解温度>酶添加量，而其最优水平组合为A2B3C3D1E1F2，即酶添加量0.8 %(W/W)、酶解时间2 h、酶解温度55 ℃、料液比1∶5(g/mL)、SSL的添加量0.6%(W/W)，辣椒粒度为粉碎后过100目筛，此时，辣椒红色素的含量为31.25 mg/g。
1 基本描述
1.1 问题描述
随着作战进程的推进及作战任务的遂行，受敌方火力打击的影响，我方作战装备不可避免地会在不同时间、不同地点出现不同程度的损伤。巡回修理力量编成多个巡回修理组，根据待修装备损伤及分布情况，有计划地对分布在战场各地域的待修装备进行巡回修理，既可相对有规律地遂行保障，又可对突然出现的待修装备做出快速反应。在一体化指挥信息平台的支撑下，各待修装备的位置、预计修理时间以及装备重要度等信息可通过相应手段获知。各巡回修理组从初始地域出发，前往各待修装备所在地域实施巡回修理，完成各自修理任务后不返回初始地域，而是根据保障指挥员的安排执行后续任务。在实施修理任务的过程中，保障指挥员会根据不断更新的待修装备及各巡回修理组信息，对各巡回修理组的修理计划进行动态调度。
在待修装备不断出现，而修理时间及修理能力有限的前提下，需考虑待修装备修理时间、重要度、修复状态的不确定性、巡回修理组修理能力的变化等因素的影响，从而为各巡回修理组合理分配维修任务，确定各巡回修理组的各待修装备的修理顺序和修复状态，并根据不断出现的待修装备及诸多不确定因素对维修任务进行动态调整，使整体的巡回修理效果达到最优化，最有利于战斗力的快速恢复及作战任务的完成。
1.2 假设描述
为了简化问题，进行如下假设：
1) 参与装备巡回修理任务调度的各待修装备均在巡回修理组的修理能力范围内，不考虑备品备件短缺情况；
姨太太。金铨有两个，金凤举子承父业，有一个晚香。姨太太一种是不问世事，安分度日。一种是趁年轻邀宠，却也是在夹缝中生存的可怜人。
2) 研究对象为合成旅巡回修理力量，编成多个巡回修理组对所属部队实施全域巡回修理；
3) 在修理任务实施过程中，各巡回修理组独立完成各自修理任务，相互独立，不存在相互支援；
4) 各待修装备的位置、预计修理时间、修理时间窗、装备重要度等相关信息已通过技术侦察获知；
5) 各修理组修理过程中不会由于任务调整而中止当前正在执行的任务；
近年来，随着我国对环保要求的提高，特别是2014年9月国家发改委、环保部、国家能源局联合发布了《关于印发〈煤电节能减排升级与改造行动计划（2014——2020年）〉的通知》，并正式将新建燃煤发电机组污染物推到了大众的视野之内，这也直接拉开中国燃煤电厂实施超净排放的序幕，湿法脱硫出口烟气携带的雾滴、气溶胶等污染物受到前所未有的关注。湿式电除尘器布置在湿法脱硫吸收塔出口、烟囱入口，可对吸收塔出口的污染物进行 “终端拦截”，进一步降低烟囱出口排放，达到“超净排放”的标准。
6) 各巡回修理组完成修理任务后不返回初始地域，而是等待保障指挥员的后续指挥调度；
7) 待修装备修复后直接归建作战部队并参与作战，忽略归建时间；
8) 不考虑敌方火力打击对待修装备造成的二次伤害。
2 模型构建
2.1 参数定义
为方便研究，对相关参数定义如下：
1)T sta为战斗开始时刻，T beg为巡回修理开始时刻，T end为战斗结束时刻。
2)E 为待修装备集合，且|E |=n ，n 为待修装备数量，H 为始发点集合。
3)R 为巡回修理组集合，且|R |=m ，m 为巡回修理组数量，
为第k 个巡回修理组初始位置，记坐标为
4)t i 为待修装备i 出现的时刻，(x i ,y i )为待修装备i 的位置坐标，其修理时间窗为[t i ,T i up]，其中T i up为修理窗时间上限。
考察修正后最终模型的主要拟合度指标，依据Amos Output中的Model Fit指标显示，显著性概率值P=0.132>0.05,卡方自由度比2/df=1.477<2.00，近似误差均方根RMSEA=0.049<0.08,拟合优度指数GFI=0.980>0.9，校正拟合优度指数AGFI=0.936,基准化适合度指标NFI=0.971>0.9，比较适配度指标CFI=0.990>0.9，综合考虑各项拟合指标，该结构模型整体拟合情况良好。
5)S 1为待修装备修复至“能正常作战”状态，其受损功能得到完全修复；S 2为待修装备修复至“能应急作战”状态，其受损功能并未完全修复；
分别为待修装备i 修复至S 1、S 2状态时的装备重要度。
2.4.4 非遍历性
7)T i 为待修装备i 的修竣时刻，T i c为计划修复待修装备i 所耗时间，
为实际修复待修装备i 所耗时间，T i S1 、T i S2 分别为将待修装备i 修复至S 1、S 2状态时的预计修理工时。
8)p k 为第k 个巡回修理组的修理能力，即每小时能提供的人时。
为第k 个巡回修理组执行第s 次修理任务时的修理效率,β 为修理效率衰减系数。
10)α 为惩罚系数，用以表征时间窗对装备重要度的影响。
为0-1变量，当待修装备i 由巡回修理组k 修复，则
否则
（c）急性妊娠期脂肪肝/HELLP综合征导致的肝衰竭：建议立即终止妊娠，如果终止妊娠后病情仍继续进展，需考虑人工肝和肝移植治疗（Ⅲ）。
为0-1变量，当巡回修理组k 由待修装备i 驶往待修装备j ，则
否则
13)O =[o 1,o 2,…,o m ],为维修任务调度问题的一个可行解，o k 为第k 个巡回修理组的一个规划修理路径，代表该巡回修理组开展修理任务的装备编号序列；|o k |为该规划序列中的元素个数，即该规划路径中的装备数；
为该规划路径中第l 台装备。
14)q k 为第k 个巡回修理组修理路径中的截点，表示该组最终实际修理的最后一台装备，且有
(1) 遍历时，q k =|o k |且
(2) 非遍历时，q k +1≤|o k |且
（1）“一带一路”沿线各国物流绩效水平的改善会拓宽中国机械运输设备出口的种类，增加中国机械运输设备出口的数量，提高中国机械运输设备出口的价格。毫无疑问，无论种类、数量还是价格都有利于中国机械运输设备出口总量的增长。因此，中国应该积极推进与 “一带一路”沿线各国的基础设施投资合作，改善贸易伙伴国的物流绩效水平，与 “一带一路”沿线各国实现设施联通、贸易畅通。
1）将药品审评审批、质量标准管理和药品监管（例如：药品注册管理、药品说明书的管理[9]等）有机结合起来，将药品质量标准的发布、修订、勘误等信息及时在国家标准管理部门网站公布，以便相关企业、药检机构和社会公众查询、执行和使用，切实提高质量标准的可及性。
2.2 多目标参数确定
实施战时装备巡回修理旨在最大可能地修复待修装备，使其能尽快返回战场继续参与战斗，以发挥修竣装备对装备体系的贡献，使作战部队的战斗力能得到最大程度的恢复。因此，可从以下3个方面来综合衡量装备巡回修理任务调度的优劣。
1) 修竣装备总数。修竣装备总数是整个战斗过程中各巡回修理组修复的待修装备数量总和，它反映了巡回修理的快慢程度，直接影响装备的参战率，其表达式为
(1)
2) 修竣装备总重要度。修竣装备总重要度是修竣装备重要度的总和，反映了所修竣装备对整个装备体系的贡献度，是巡回修理重要性的反映，其表达式为
(2)
结合巡回修理特点，考虑各巡回修理组修复各待修装备的修理时间关系、遍历与非遍历特点、待修装备修竣与否，确定相关约束条件，构建装备巡回修理任务调度模型如下:
(3)
以上3个参数从不同侧面反映巡回修理效果，且相互影响、相互制约，需要构建多目标优化模型来实现全局的综合权衡。
2.3 调度模型建立
3) 二次作战总时间。二次作战总时间是指各修竣装备从修竣时刻至战斗结束时刻的时间长度的总和，反映了装备修竣后返回战场参战的有效时长，是修理及时程度的体现，其表达式为
(4)
(5)
(6)
其中，式(4)-(6)表示面向装备巡回修理的任务调度目标分别为修竣装备总数最大、修竣装备总重要度最大、获得的二次作战总时间最长；式(7)表示巡回修理组从初始位置前往并修竣该组计划内第一台待修装备的时间关系；式(8)表示巡回修理组修理计划内相邻2台待修装备的修竣时刻的约束关系；式(9)表示修理时间以及修理能力的约束关系；式(10)表示巡回修理组的遍历与非遍历约束；式(11)表示每个待修装备至多仅由一个巡回修理组修理一次；式(12)、(13)表示任意一条弧线的终点(起点)待修装备有且仅有一个起点(终点)待修装备与之相连。
2.4 调度模型分析
2.4.1 多目标分析
装备巡回修理任务调度模型从修理数量、修理对象重要程度以及修理及时性3个方面权衡维修任务调度方案的优劣，能够克服传统维修任务调度单纯追求某一目标而导致调度方案在其他需求方面存在较大偏差的问题，从而能够更加全面地权衡维修任务调度方案的优劣。
Pareto解集[14]通过支配关系判断调度方案的优劣，能够实现各调度目标间的均衡，是处理多目标调度问题的较好思路。通过在多个目标中协调平衡，可获得一组可接受解(即Pareto最优解集)，从而增加了保障指挥员的决策余地。保障指挥员可根据战场实际需求及决策偏好从Pareto最优解集中选择满意的调度方案。
设定满意方案的决策策略为通过对Pareto最优解集中的解进行规范化处理，采用加权法进行排序优选，从而选择满意解，即
(14)
式中：F ′为规范化后的解，
分别为F 1、F 2、F 3的规范化解；μ 1、μ 2、μ 3分别为保障指挥员的偏好权重系数，可根据战场需求及态势进行调整。
2.4.2 修复状态的不确定性
巡回修理相对于伴随修理而言，修理开始的时间相对滞后、待修装备相对较多且损伤程度更严重。由于修理时间的限制，为了使待修装备尽快得到修复并返回战场参战，待修装备的修复状态并不全是完全修复，而有可能是仅恢复了部分功能，即待修装备的修复状态分为“能正常作战”S 1和“能应急作战”S 2。
待修装备修复至不同状态对维修任务调度均产生较大影响，主要是所需的修理工时和修复后获得的装备重要度不同。因此,在制定调度方案时，不但要确定各待修装备的修理顺序，还需确定各待修装备的修复状态，这使得维修任务调度问题更为复杂，求解难度更大，但同时也更贴近实际，具有更强的实用性。
2.4.3 修理时间窗
总之，在英美治安法庭模式下，立法设置了两套体系完全独立和明显区别的程序，即小额程序与普通程序，［5］前者简便、快捷、非正式、强调本人参与不依赖律师，成本低廉，能够快速解决小额纷争，但却蕴含着高度风险性和事后审查的困难;后者程序复杂、时间冗长、技术性强、依赖律师、成本高昂，但具有程序的安定性和结果的可预测性。两套系统以及两种程序的区分度极为显著，当事人可以根据自己的情势理性决定是否选择小额程序。
(15)
式中：
为经修理时间窗惩罚后的待修装备i 修复至状态S p 的装备重要度。
为待修装备i 与j 之间的距离，且d ij =d ji ，考虑到地形因素的影响，引入迂回系数
为巡回修理组的平均机动速度。
巡回修理组负责全域的巡回修理，待修装备多、修理难度大，而受修理时间及修理能力的限制，难以对散布在战场的各待修装备实现全部修复。在调度初期，参与调度的待修装备较少，在任务规划时一般可实现遍历修理(即对全部待修装备进行修理)，但随着作战进程的推进，待修装备不断增多，无法对参与调度的众多待修装备实现全部修复，需采用非遍历型调度。
方案2，在实验室栽植收割后的小麦茎秆，按田间种植行距进行模拟，栽培分布图如图5所示，单位cm。用该测量系统对模拟区域进行反复测量，设定一定的电机转速使滑块在导轨上匀速运行，正反向连续重复测试10次。
遍历型维修任务调度的约束条件为
且q k =|o k |，非遍历型维修任务调度将约束条件松弛为
且q +1≤|o k |。巡回修理装备维修任务调度既包含遍历型调度(前期)，又包含非遍历调度(后期)，在更加符合战时巡回修理实际的同时，也大大增加了模型求解难度，使求解算法难以收敛。
2.4.5 修理能力变化
目前肾病内科住院患者逐年增多,多数患者病情严重,且住院期间存在较大的安全隐患,这些安全隐患可能来自患者自身、护士、环境、用药等,通过针对不同因素强化管理和安全防范,可有效规避和减少相关不良事件的发生,保障患者安全,提高患者生存质量[4-5]。本研究中,对照组用常规护理管理方法,观察组在对照组基础上采用安全隐患排查和积极干预。结果显示,观察组护理满意度高于对照组,P<0.05;观察组护士安全管理意识、安全防范完善性优于对照组,P<0.05;观察组安全事件发生率低于对照组,P<0.05。
场上，艾尔一拳狠狠地打到威尔的下巴上。我往回缩了一下。艾瑞 克站在房间另一边，一脸假笑地看着艾尔，不断用手翻动他的眉环。
(16)
式中：
为执行第s (s ≥2)次修理任务时的修理效率，且ξ 1=1。
待修装备的修理时间不但与所承担修理的巡回修理组有关，还受其修理顺序的影响，导致装备维修任务调度更加复杂。
修理人员疲劳所导致的修理能力变化可通过修理效率来度量，其计算公式为
2.4.6 动态驱动策略分析
战时各巡回修理组的修理能力存在差异，且受战场环境等诸多因素影响，各巡回修理组的修理能力会随着修理任务的遂行而逐渐发生变化。主要体现在随着修理任务的实施，修理人员会产生疲劳，从而影响修理效率，使巡回修理组的修理能力发生变化。
在处理装备维修任务动态调度问题时，引入滚动时域思想[11]，根据调度需求将整个任务剖面转换为一系列离散时间点的静态调度问题，通过多次维修任务调度，实现动态问题静态处理。而动态驱动策略就是设定某一驱动条件用以判断该时刻是否需要对维修任务进行再次调度，是动态调度的基础。
设定面向装备巡回修理任务动态调度的动态驱动策略为：当新出现的待修装备数量达到一定阈值
时，则根据该时刻各巡回修理组信息、待修装备信息进行一次重调度。在重调度时刻将巡回修理组正在修理或刚修竣的待修装备记为关键点。
该动态调度策略在每次重调度时，更新巡回修理组及参与调度的待修装备信息，消除待修装备的实际修理时间与计划修理时间差所产生的调度误差，在实现维修任务动态调度的同时也提升了动态调度的可靠性。
3 模型求解
装备巡回修理任务调度模型是一个存在诸多约束的多目标调度模型，针对模型特点，笔者结合Pareto解集思想，设计了改进遗传算法来求解装备巡回修理任务多目标动态调度问题。
3.1 Pareto最优解集构建
传统多目标优化方法通过线性加权或目标规划等方法处理多目标调度问题，其本质是将多目标问题转化为一个或一系列单目标优化问题进行求解，难以解决多个优化目标之间的冲突关系，存在依赖先验知识、难以处理Pareto最优前端非凸等问题。NSGA-II算法作为最优秀的多目标进化算法之一，在保证种群多样性及保护种群优良个体的同时降低了计算复杂度[16]。因此，笔者通过NSGA-II算法的精英策略，采用非支配排序方法并结合拥挤度比较算子，获得Pareto最优解集，为保障指挥员提供决策依据。保障指挥员可根据战场态势及实际需求，依托决策偏好及决策策略，在Pareto最优解集中选择最满意解。
3.2 编码与解码设计
考虑装备巡回修理任务调度具有多个巡回修理组、非遍历性以及待修装备修复状态不确定性等特点，设计3段式编码以实现相关约束：前段采用顺序编码，用于表示巡回修理组的任务分工；中段采用“1-2”整数编码，用于表示待修装备的修复状态，1表示修复状态为S 1，2表示修复状态为S 2；后段采用整数编码，表示断点位置，用于划分各巡回修理组。该编码方式能实现多个巡回修理组的非遍历约束以及待修装备修复状态不确定性约束，且染色体与解一一对应，避免了遗传操作中不可行解的产生，大大提高了算法收敛速度。解码时，选取目标函数作为适应度函数，通过相关约束，计算求得截点的信息及相应适应值，从而实现解码。
编码及解码示例如图1、2所示，其中染色体X =(4,8,2,6,3,10,7,5,1,9,1,1,2,1,2,2,1,2,1,1,3,7)，n =10，m =3。经解码可知，共有3个巡回修理组负责10台待修装备的修理任务，任务计划安排为：巡回修理组1的任务安排及修理顺序为4-8-2，修复状态分别为S 1、S 1、S 2；巡回修理组2的任务安排及修理顺序为6-3-10-7，修复状态分别为S 1、S 2、S 2、S 1；巡回修理组3的任务安排及修理顺序为5-1-9，修复状态分别为S 2、S 1、S 1。截点信息为(1，5，8)，表明受诸多约束的影响，部分待修装备并未来得及修复，最终实际修理任务完成情况为：巡回修理组1修竣待修装备4，修复状态为S 1；巡回修理组2修竣待修装备6-3，修复状态为S 1、S 2；巡回修理组3修竣待修装备5，修复状态为S 2。适应值分别为F 1、F 2、F 3。
图1 编码示例
图2 解码示例
3.3 遗传算子设计
3.3.1 父染色体选择
根据NSGA-II算法中的非支配排序和个体间拥挤距离得到比较算子，采用BTS(Binary Tournament Selection)从上一代染色体中选取20%的个体作为父染色体进行后续交叉变异操作。
3.3.2 交叉变异
由于染色体采取3段式编码，各段编码的编码方式及其实际意义不同，其交叉变异无法采用传统方式进行。为了增大种群多样性，提高收敛速度，针对编码特点，确定“前段只进行变异操作，中段和后段可进行交叉及变异操作”思路，设计了相应遗传算子，对每条父染色体进行如下遗传操作：1)前段、后段均采用随机更新操作；2)前段采取随机更新，后段采取倒置操作；3)前段采用倒置操作，后段采用随机更新操作；4)前段采用倒置操作，后段采用滑动平移操作；5)前段不采取操作，后段采用随机更新操作。
3.4 算法流程设计
模型求解步骤如下：
1) 初始化相关参数。主要是确定种群规模pop_size、最大迭代次数num_gen等，并随机产生初始种群P 0。
2) 对种群P 0中任一染色体进行解码，计算其适应值和截点信息，得到初始化可行解种群O ，并记为pop_chrom。
3) 根据NSGA-Ⅱ算法对可行解种群O 进行快速非支配排序，计算非支配集个体间拥挤距离，并令gen=1。
由于修理能力的限制，难以实现所有待修装备在第一时间均得到修复。随着作战阶段、作战任务的变化，当超过某一时刻待修装备仍未得到修复，该待修装备的重要度将会受到影响。因此借鉴车辆路径问题中的时间窗[15]概念，引入待修装备修理时间窗，对超出修理时间窗而未得到修复的待修装备重要度进行惩罚，其惩罚公式如下：
4) 根据BTS从pop_chrom中随机选出数量规模为pop_size的父代染色体种群parent_chrom。通过遗传算子进行遗传操作，产生子代染色体种群offspring_chrom，并计算offspring_chrom中任一染色体的截点信息及其适应值。
5) 采用NSGA-Ⅱ算法对pop_chrom及offspring_chrom进行快速非支配排序，计算非支配集个体间拥挤距离，根据BTS从pop_chrom及offspring_chrom中筛选出规模为pop_size的新一代染色体pop_chrom，从而实现父代优秀个体基因的精英保留。
6) 判断是否gen<num_gen。若成立，gen=gen+1，转入步骤4)继续寻优；否则，停止迭代，获得本次维修任务调度的Pareto最优解集，并依据决策策略从中选出满意解，转入步骤7)。
7) 判断重调度驱动策略是否满足。若满足，转入步骤1)进行重调度；否则运算终止，输出调度结果。
4 示例仿真与分析
4.1 示例仿真
某合成旅受上级指示执行机动进攻任务，受敌方火力打击，经过一段时间后，陆续出现待修装备，在一体化指挥信息网络的支撑下，经技术侦察，各待修装备相关信息均已获知，保障指挥员根据已知信息，派遣3个巡回修理组对全旅实施全域保障，各巡回修理组遂行巡回修理任务，负责修理120 min内能完成的待修装备。
设作战开始时刻为
由于战时保障的延迟性，作战开始第50 min开始实施巡回修理，此时各巡回修理组位置为(-2.1 km，0.2 km) 、(0.5 km，0.6 km)、(3.7 km，0.3 km)。待修装备位置信息、修理时间窗、恢复至不同状态的预期修理时间等相关信息如表1所示。
战斗初期(t ≤300 min)，需要更多、更重要的装备及时参战，取μ 1=0.2，μ 2=0.6，μ 3=0.2，战斗后期(t >300 min)，需要提供更多的作战时间，取μ 1=0.3，μ 2=0.2，μ 3=0.5。
根据巡回修理组的实际机动时间、修理时间以及调度策略对装备维修任务进行动态调度，得到调度规划结果如表2所示。
表1 待修装备信息
表2 调度规划结果
注：
为各组实际调度目标值;
为整体实际目标值;
为关键点的计划修竣时刻;
为关键点的实际修竣时刻。
此次装备巡回修理任务调度共经历了6次调度过程，在MATLAB软件平台用时分别为11.45、10.87、9.68、9.75、8.32、9.06 s。最终装备巡回修理路径规划如图3所示。
图3 装备巡回修理路径规划
4.2 结果分析
由表2可以看出：
1) 从巡回修理组开始实施巡回修理至战斗结束，共修复装备26台，获得的装备重要度总和为12.17，二次作战时间为6 099 min，约101 h，对恢复部队战斗力提供了较有力的保障，间接证明了战时巡回修理的重要性。
2) 调度时间均在12 s以内，满足战时装备维修任务调度的实效性要求，同时也证明了所构模型和算法的可行性。
3) 在动态驱动策略的影响下，一共进行了6次调度，实现了待修装备不断出现情况下维修任务动态调整，其中前5次调度为遍历型调度，第6次调度为非遍历性调度；待修装备20、28、29、30未纳入最终的巡回修理任务规划中，这主要是由于战时修理时间的限制导致修理任务无法全部完成，反映了巡回修理的非遍历性。
4) 在第105 min进行第2次调度时，待修装备8纳入巡回修理组3的修理计划中，并安排在巡回修理组3的第3顺位，但在第206 min进行第3次调度时，巡回修理组3还在修理第2顺位的待修装备7，待修装备8参与重调度，并最终将其规划到巡回修理组3的第4顺位，并在第335 min得到修竣，类似的还有装备9、16、20、22。这是动态调度的体现，反映了每个巡回修理组的任务序列对其自身而言并不一定是最优的，但是对于整体的巡回修理效益是最优的。这也间接证明了通过动态调度能够更有利于全局的最优化，体现了动态调度的优势。
5) 在第206 min进行第3次调度时，待修装备9被纳入巡回修理组3的第5修理顺位，且计划修复状态为S 1，但是在第4、5、6次重调度时，待修装备9均参与了重调度，且最终在巡回修理组1的第7修理顺位得到修复，且修复状态为S 2。这反映了在重调度时刻，会根据该时刻的全局信息对各巡回修理组的任务序列及各待修装备的修复状态进行动态调整，以实现全局的最优化，从侧面反映了调度模型的合理性。
6) 前3次调度，待修装备的修复状态以S 1为主，因为前期待修装备较少，预计的修理时间相对较充裕，而μ 2取值较大，将装备修复至S 1状态能获得更多的重要度，且前期修竣装备参战时间相对较长，更高的重要度有利于发挥作战效能。后3次调度，待修装备的修复状态多为S 2，因为此时参与调度的待修装备数量多，且在战斗中后期，μ 3取值较大，希望能获得更多装备二次作战时间，因此将待修装备修复至S 2状态既能节约修理时间以便完成更多修理任务，又能迅速使其返回战场获得更多的二次作战时间。
由图3可以看出：巡回修理组1的规划路径中出现了折线3-6-5以及15-21-9-25-26，前者是因为装备5相较于装备6，修复到相同状态的装备重要度高且修理时间短，因此巡回修理组1在修竣装备3后，舍弃较近的待修装备6，而先修理重要度高且更易修复的待修装备5，这满足“先修重要装备”以及“先修易修装备”的战时修理要求，也从侧面反映了调度模型的合理性。而后者是因为该规划是多次重调度的综合结果，装备15在第4次调度后得到修竣，装备21、9是在第5次调度后得到修竣，而装备25、26是在第6次调度才纳入调度计划并最终得到修竣，这也从侧面反映了动态调度的核心是使整体修理效益最优化，而不是追求某一巡回修理组的修理效益最优。
5 结论
合理高效的装备巡回修理任务动态调度模型，可以为战时保障指挥员的维修决策提供数学支撑，并能根据战时不断变化的维修需求作出及时有效的调整，大大降低了决策工作量和人为决策风险。笔者在提出巡回修理装备维修任务动态调度军事问题的基础上，构建了维修任务的多目标动态调度模型，考虑待修装备修复状态的不确定性，引入待修装备修理时间窗、非遍历性以及修理能力的变化等约束条件，使模型更具现实意义并增强了模型的合理性，通过制定动态调度策略，实现了维修计划的动态调整，设计了NSGA-Ⅱ改进遗传算法求解了模型，并通过示例验证了模型及算法的合理性。下一步将对定点修理中装备维修任务动态调度问题展开研究。
参考文献：
[1] 张立民,彭乐,刘敬虎,等.求解装备保障任务调度问题的改进粒子群算法[J].计算机工程与设计,2012,33(8):3105-3110.
[2] 王正元,朱昱,宋建社,等.动态维修任务调度的优化方法[J].机械工程学报,2008,44(1):92-97.
[3] 陈春良,昝翔,张仕新,等.基于改进MMAS的装备维修任务路径规划方法[J].系统工程与电子技术,2017,39(12):2716-2720.
[4] 陈伟龙,陈春良,史宪铭,等.基于变体GA的进攻作战抢修任务动态调度[J].系统工程与电子技术,2017,39(3):577-583.
[5] 吕学志,陈乐,尹健,等.考虑休息的维修任务调度模型及其求解算法[J].兵工学报,2014,35(12):2116-2123.
[6] PILLAC V,GENDREAU M,GUERET C,et al.A review of dynamic vehicle routing problems[J].European journal of operational research,2013,225(1):1-11.
[7] 朱昱,宋建社,王正元.一种基于最大保障时间的战时装备维修任务调度[J].系统工程与电子技术,2007,29(11):1900-1903.
[8] 昝翔,陈春良,张仕新,等.多约束条件下战时装备维修任务分配方法[J].兵工学报,2017,38(8):1603-1609.
[9] 陈春良,陈伟龙,陈康柱,等.考虑非遍历的抢修任务多目标动态调度[J].兵工学报,2017,38(8):1593-1602.
[10] 朱昱,宋建社,曹继平,等.一种考虑装备维修流程的多维修任务调度[J].系统工程与电子技术,2008,30(7):1366-1369.
[11] 王正元,严小琴,朱昱,等.一种考虑专业的动态维修任务调度的优化方法[J].兵工学报,2009,30(2):252-256.
[12] 陈伟龙,陈春良,陈康柱,等.考虑不确定性的进攻作战抢修任务动态调度[J].兵工学报,2017,38(5):1011-1019.
[13] 曾斌,姚路,胡炜,等.考虑不确定因素影响的保障任务调度算法[J].系统工程与电子技术,2016,38(3):595-601.
[14] 汪开普,张则强,邹宾森,等.模糊作业时间的拆卸线平衡Pareto多目标优化[J].计算机工程与运用,2018,54(1):257-263.
[15] TANG Y L,CAI Y G,YANG Q J.Improved ant colony optimization for multi-depot heterogeneous vehicle routing problem with soft time windows[J].Journal of Southeast University(English Edition),2015,31(1):94-99.
[16] DEB K,PRATAP A,AGARWAL S,et al.A fast and elitist multi-objective genetic algorithm:NSGA-Ⅱ[J].IEEE transaction on evolutionary computation,2002,6(2):182-197.
Multi -objective Dynamic Scheduling with Circuit Repair Task under Uncertain Repair Status
LIU Yan1, CHEN Chun-liang1, CHEN Wei-long2, FAN Jiang-bo3
(1. Equipment Support and Remanufacturing Department, Army Academy of Armored Forces, Beijing 100072, China;2. China Satellite Maritime Tracking and Control Department, Jiangyin 214431, China;3. Troop No.66058 of PLA, Tianjin 301700, China)
Abstract ：To deal with the problem of heavy tasks and limited repair time of equipment circuit repair in wartime, single scheduling target and simple constraint of the existing wartime equipment maintenance task scheduling research, the multi-objective dynamic scheduling study of equipment circuit repair task is carried out. Circuit repair tasks scheduling with the complex constraints, which is a military problem, is proposed. Considering the uncertainty of the repair status of the equipment to be repaired, the constraints of repair time window, non-traversal, and variety of repair capability are introduced besides repair capability and repair time limit. A multi-objective dynamic scheduling model is established regarding the total number of repaired equipment, the sum of repaired equipment importance degree, and the maximum total time of the second operation as scheduling targets. An improved multi-objective genetic algorithm based on Non-dominated Sorting Genetic Algorithm Ⅱ (NSGA-Ⅱ) is designed to solve the model. Through the case simulation and analysis, the rationality and effectiveness of the model and algorithm are verified.
Keywords :circuit repair; multi-objective dynamic scheduling; uncertainty of the repair status; repair time window
中图分类号： E92
文献标志码： A
DOI： 10.3969/j.issn.1672-1497.2019.01.004
文章编号： 1672-1497(2019)01-0016-09
收稿日期： 2019-01-05
作者简介： 刘 彦(1994-)，男，博士研究生。
(责任编辑:王生凤)
标签：巡回修理论文; 多目标动态调度论文; 修复状态不确定性论文; 修理时间窗论文; 陆军装甲兵学院装备保障与再制造系论文; 中国卫星海上测控部论文; 66058部队论文;