基于麻雀搜索算法考虑病患流的医院诊疗设备维护优化研究.pdf

文章编号：1007 6735(2023)05 0513 10DOI:10.13255/ki.jusst.20220317004基于麻雀搜索算法考虑病患流的医院诊疗设备维护优化研究李岱高，刘勤明，李佳翔（上海理工大学管理学院，上海200093）摘要：针对医院诊疗设备出现故障后对病人的影响不可控、产生的维修费用高等问题，提出了一种基于病患流分布不均的预防性维护决策模型。首先，针对医院诊疗设备退化速度受病患流分布影响的特点，采用病患流因子反映其人流量变化规律，并基于设备衰退演化规则构建了设备退化过程模型。其次，通过将量化的诊疗设备故障风险成本以及延迟或提前维护带来的变动成本纳入成本计算体系，提出了病患流分布不均状况下以故障风险值为约束的设备维护策略。最后，采用麻雀搜索算法进行寻优获取关键初始变量，并对模型进行数值仿真分析，结果表明，新维护策略在保障诊疗设备可靠度的同时，也能保持较低的维护成本。关键词：诊疗设备；病患流分布；故障风险成本；变动成本；麻雀搜索算法中图分类号：TP207文献标志码：AMaintenance optimization of medical equipment considering patientflow distribution based on the sparrow search algorithmLI Daigao，LIU Qinming，LI Jiaxiang(Business School,University of Shanghai for Science and Technology,Shanghai 200093,China)Abstract:Tosolvetheproblemsofuncontrollableimpactonpatientsandhighmaintenancecostsafterthefailureofhospitaldiagnosisandtreatmentequipment,apreventivemaintenancedecision-makingmodelbasedonunevendistributionofpatientflowwasproposed.Firstly,inviewofthecharacteristicsofthehospitaltreatmentequipmentdegradationrateaffectedbythedistributionofpatientflow,thepatient flow factor was used to reflect the changing law of its human flow,and the equipmentdegradation process model was constructed based on the equipment degradation evolution rule.Secondly,byincorporatingthequantifiedcostofdiagnosisandtreatmentequipmentfailureriskandthevariablecostcausedbydelayedorearlymaintenanceintothecostcalculationsystem,anequipmentmaintenancestrategythatconstrainsthefailureriskvalueundertheunevendistributionofpatientflow上 海 理 工 大 学 学 报第45卷第5期J.UniversityofShanghaiforScienceandTechnologyVol.45No.52023收稿日期：20220317基金项目：国家自然科学基金资助项目（71632008，71840003）；上海市自然科学基金资助项目（19ZR1435600）；教育部人文社会科学研究规划基金资助项目（20YJAZH068）；上海理工大学科技发展项目（2020KJFZ038）；2021 年上海市大学生创新创业训练计划项目（SH2021078）.第一作者：李岱高（1996），男，硕士研究生研究方向：设备维护与决策E-mail：通信作者：刘勤明（1984），男，特聘教授研究方向：维护调度、人工智能等E-mail：wasproposed.Finally,throughthesparrowsearchalgorithm(SSA)optimizationtoobtainthekeyinitialvariables,andnumericalsimulationanalysisofthemodel,theresultsshowthatthemaintenancestrategyproposed in this paper can ensure the reliability of the diagnosis and treatment equipment whilemaintainingalowmaintenancecost.Keywords:medical equipment;distribution of patient flow;failure risk costs;variable costs;sparrowsearch algorithm在诊疗活动中，就诊患者人流量受节假日、工作时段等因素影响呈现较大波动，如周末、春节假期尾声等就诊高峰期，经常出现病患集中就诊，导致诊疗设备连续不停长时间运转，而非就诊高峰期内设备运行时间间隔明显增加。病患流时间分布的差异导致诊疗设备运行工况不同；同时，病患流分布不均情况下设备运行时间段的收益不同，进而导致不同时刻停机维护的成本产生较大差异。因此，在设备不同病患流运行工况下，保证诊疗设备可靠度的同时降低维护成本成为亟待解决的问题。目前，医院各种类型诊疗设备的维护和保养常根据历史经验和维护手册进行1，易造成维护过剩或维护不足，增加医疗风险。中国国家医疗器械不良事件监测年度报告（2020 年）显示，2020 年共收到医疗器械不良事件报告 536055 份，相较于 2019 年增长 35.25%。对于医院诊疗设备的预防性维护，不少国内外学者都作了相应的研究。Carmignani2提出了基于综合成本的故障影响分析 FEMA 方法；AbbasgholizadehRahimi 等3用FEMA 方法对放射性诊疗设备的故障进行优先排序；Hernndez-lpez等4采用层次分析法确定医院诊疗设备维护优先级，并给出相应设备权重便于后续评估。然而，他们的研究局限于患者本身的安全性及设备对医院的重要性，忽视了故障维修成本和停机损失成本等客观成本对维护优先级的影响。张宁等5建立了基于 NSGA-的医疗设备预防性维护与病人调度联合优化模型。但文献 2-5均未考虑设备运行工况对于设备衰退速率的影响，不能较好地考虑病患流分布不均对诊疗设备维护费用的影响。研究不同工况下的设备维护，首先需要建立时变工况与设备衰退规律的集成模型。沈南燕等6建立了考虑总成本最小化的设备故障率演化模型，并结合支持向量机对模型有效性进行检验；Xia 等7提出了基于环境变化因子和维修效果的设备预防性维护模型，将外部环境变化对故障率递增的影响纳入模型，但二者并未考虑设备历史工况对设备故障率的影响；Li 等8提出了考虑故障率刷新因子的设备维护策略，用以求解定期维护和预防性维护相结合的设备故障率递增模型；Vilarinho 等9将更换部件的成本纳入预防性维护成本中，从而选择更经济的方式对设备进行维护；甘婕等10以长期平均维修费用率最小为目标建立了随机规划模型；盖京波等11考虑了设备使用时间有限条件下设备预防性维修优化模型；汤乐成等12提出了考虑交付时间、产品质量和企业信誉的联合设备维修策略。以上文献均采用了加速失效模型描述设备退化规律，但在医疗设备领域研究较少。基于上述研究中存在的问题，本文考虑诊疗设备服役过程中病患流分布不均对其退化速度的影响，构建了病患流分布影响下故障率模型；有别于其他研究将故障风险成本与时间变动成本分开研究，本文同时将维护风险、设备重要性等故障风险以及维护时间变动产生的变动成本纳入成本体系，以维护总成本最小为优化目标构建了不同病患流分布下的设备维护调整策略，以满足诊疗设备运行高可靠要求和降低维修费用的需求。最后以某型号的 CT 机服役数据为例进行仿真分析，并采用麻雀搜索算法对维护策略进行优化，验证了其合理性和有效性。1 问题描述与假设本文的研究对象主要为大型复杂且不适合长时间连续使用的诊断治疗类医疗设备，例如 X 射线机、CT 机、核磁共振仪等，简称诊疗设备。该类设备使用寿命的衰退趋势会随其连续运行时长的增加而发生变动。例如 CT 机需要连续曝光产生514上海理工大学学报2023年第45卷射线，连续使用多次长时间曝光会导致设备过热，从而影响设备寿命。而就诊病患人流越密集越容易导致诊疗设备的连续运行，就诊排队病患流时域分布不均不仅会造成诊疗设备运行工况的差异，导致设备可靠度的衰退；同时，诊疗设备在病患流高峰期进行维护还将给医院造成较大的经济损失，需要对设备维护计划进行局部修正。如果提前对诊疗设备进行维护，则会造成一定的维护成本浪费，而延迟维护又会造成较大的失效风险。因此，本文基于病患流工况影响建立诊疗设备故障率演化模型，在病患就诊高峰期考虑可靠性及经济性，对诊疗设备制定合理的维护计划。模型基本假设：a.诊疗设备以全新状态投入使用，初始役龄为 0；b.在不同病患流量影响下设备失效机理不发生改变，服从加速失效模型；c.不考虑由于备件不足引发的维修等待，设备发生的各种故障能被及时修复；md.初始设备故障率服从两参数威布尔分布，形状参数为，尺度参数为；e.每位患者在诊断设备的就诊时间为相同的固定值；f.病患流区间内病患流因子 可以表征其真实病患流水平。2 模型建立 2.1 病患流分布影响下故障率模型构建ai病患流分布密度不同导致各区间设备连续运行时长不同，从而影响医疗设备衰退速率。为更好地表征病患流分布对设备的影响，借鉴文献 13中对环境应力的等效处理方法，将病患人流量数据进行分段等效处理，由此获得不同病患流区间内等效平均人流量。定义不同病患流区间的人流量因子 反映其变化规律，即ai=i/0（1）ii0式中：为第 个病患流区间等效平均人流量；为该诊疗设备基准人流量，即在设备设计最佳运行状态下病患人流量。tn诊疗设备在第 n 个病患流区间的等效役龄 为tn=aiti（2）ti式中，为第 i 个病患流区间时长。t0n诊疗设备在不同病患人流量影响下的等效役龄如图 1 所示，设备在经历 n 个病患流区间后的等效役龄为t0n=tn+n1i=1t0i=antn1i=1ti+n1i=1t0i（3）tt0nt1t2t3a1a2a3tntiann 个区间tn1i=1.图 1 病患流区间等效役龄折算Fig.1 Equivalent service age conversion在可靠性理论中，可靠度与故障率的关系为14R(t)=exp(wt0l(t)dt)（4）t式中：R(t)表示设备可靠度；l(t)表示设备故障率；为设备在当前预防性维护周期内的服役时长。结合式(2)和式(3)，可知在经过 n 个病患流区间后，诊疗设备的可靠度为R(t,n)=expwan(tn1i=1ti)+n1i=1t0(n1)0l0(an(tn1i=1ti)+n1i=1t0i)dt（5）l0(t)式中，为设备在基准病患流工况下的故障率函数。执行预防性维护时第 n 个病患流区间内设备的故障率函数为l(t,n)=anl0antn1i=1ti+n1i=1t0i（6）虽然诊疗设备运行状态可通过维护来改善，但修复后设备衰退速率较维护前加快，文献 15使用役龄回退因子和故障率递增因子的概念用以描述这个过程，并验证了其合理性。为刻画维护活动对设备衰退过程的影响，引入役龄回退因子和故障率递增因子 描述维护过程对设备故障率的影响，第 k 次维护后设备的故障率函数为lk(t)=k1lk1(t+k1Tk1)（7）kk式中：为第 k 个预防性维护周期内回退因子，用来描述在役龄层面维护活动对设备故障率的影响；为第 k 个预防性维护周期内故障率递增因子，用来描述不完美维护状况下设备故障率随维第5期李岱高，等：基于麻雀搜索算法考虑病患流的医院诊疗设备维护优化研究515kkTk1护次数增加而变化；和取值可根据历史维护数据进行估计；表示第 k1 和第 k 个预防性维护周期时间间隔。维护策略调整总体流程图如图 2 所示。2.2 故障风险模型故障风险的定义为故障发生的概率及故障后果的严重程度，其值由两者量化的结果表示。医院诊疗设备相较于其他普通的生产设备，发生故障造成的后果不仅包含维修成本和停机损失，还包含造成病人安全事故的舆情成本和诉讼成本等。在病患流高峰时，需要对诊疗设备进行提前维护或延迟维护，因此，为辅助决策是否提前或延迟进行维护，需要事先对诊疗设备进行故障风险评估，确定其故障风险系数，使其在保证诊疗活动安全性的情况下尽量兼顾高水平的可靠度。Crisk采用故障风险成本衡量风险后果。Crisk=PdId（8）PdId式中：为设备故障发生概率；为设备故障发生的后果。诊疗设备在第 n 个病患流区间内发生故障的概率16为Pd=wt0lk(t,n)dt（9）故障后果由故障维修成本和故障风险因子构成，即Id=cdn（10）cdn式中：为单位故障维修成本；为故障风险因子，用于衡量设备发生故障后产生的风险。n提前评定医疗设备的故障风险需考虑设备发生故障带来的负面舆情、医院声誉损失、补救成本等多方面因素，参与评估的要素根据设备实际服役过程中关键因素确定4。由于涉及模糊环境下考虑多因素的影响问题，因此采用模糊综合评价法提前确定好的值。U=u1,u2,ui,unV=v1,v2,vi,vnW根据大型诊疗设备的特点分别确定因素集和评语集，据此构成评价矩阵为W=w11w12.w1mw21w22.w2m.wn1wn2.wnm（11）开始专家评估获取模糊综合评价权重集 A,得到维护点故障风险值 i借助算法工具求解优化模型,获得最优故障风险阈值 0 和设备可靠度阈值 R0获得下一个试行维护时间 tk人流量因子 ai1?Ttmax?0?RR0?CtCy?否否否否是是是是是否维护周期提前至 tkt维护周期推迟至 tky输出维护最优时间 tk结束累计运行时长 T=T+tk不进行维护调整图 2 考虑病患流分布的医院诊疗设备维护策略流程图Fig.2 Flow chart of medical equipment maintenance strategy516上海理工大学学报2023年第45卷wik模糊关系矩阵中各个元素为wik=wiknk=1wik（12）A=a1,a2,ai,anB通过层次分析法确定因素集中各影响因素的权重17，获得权重集，合成模糊综合评价结果矢量 为B=AR=(b1,b2,bi,bn)（13）i=max bi（14）Rbii式中：为模糊关系矩阵；为第 个指标的模糊综合评价结果。i由于医院诊疗设备性质特殊，发生故障的后果相较于其他行业往往更为严重，评判故障风险时需为后果最为严重的情况作准备，故直接取其中最大权重值作为该点故障风险值。2.3 变动成本模型PM(n1)PMnPM(n+1)n1，n，n+1PMnynPMntn0根据假设 c，维护活动及时开展，因此每次维护的停机时间相同。图 3 为维护时刻调整示意图，其中，分别为不考虑维护调整时的第次预防性维护时刻，为延迟维护时第 次预防性维护时刻，为维护提前时第 次预防性维护时刻。为避免停机损失过大，在就诊高峰期时需对诊疗设备进行提前或延迟维护。提前对诊断设备进行维护会造成其使用价值的浪费，而维护活动延迟进行则会产生一定的设备失效风险，严重情况下可能导致误诊或无效诊疗甚至更严重的医疗事故。为保证诊疗活动的安全性，延迟维护必须在故障风险限制范围内进行，因此，规定诊疗设备故障风险不能高于阈值。停机成本服役时长PMn1PMn+1PMntPMnPMny图 3 维护周期调整示意图Fig.3 Schematic diagram of maintenance cycle adjustmentCCtCy进行设备动态维护时的变动成本包括设备提前维护变动成本与延迟维护变动成本。提前进行诊疗设备维护可以有效保障诊疗活CL动的安全性，但会造成其使用价值产生一定的浪费，从而产生一定的经济损失。由于在诊疗活动中对诊疗设备可靠度要求较高，因此将可靠度阈值考虑到维护策略的成本建模中十分必要18。设备提前维护价值浪费成本的计算公式为CL=cr(RtR0)/(1R0)（15）crRtR0式中：为每次维护的可靠度利用价值；为提前维护时刻设备可靠度；为设备维护可靠度阈值。CY诊疗设备维护延迟将会导致设备故障风险高于预期，其惩罚成本为CY=cy(R0Ry)/R0（16）cyRy式中：为延迟维护的基本惩罚成本；为设备延迟维护时刻可靠度。Cd维护调整会使诊疗设备在不同周期内的故障率发生变动，其故障率变化引起的故障维修变动成本为Cd=Cdwt0lk(t,n)dt（17）Cdt式中：为单次故障维修成本；为设备进行维护调整时其维护时间调整量。CtCy结合式(15)、式(16)和式(17)，部件提前维护变动成本与延迟维护变动成本可由以下两式联立计算得出。Ct=CLCd（18）Cy=CY+Cd（19）J假设诊断设备在本次维护周期内进行了 次维护调整，则其由于动态维护调整引起的变动成本为C=Ji=1(Ct+(1)Cy)（20）=1=0式中，为两种调整成本的决策变量，当维护提前时取，延迟维护时取。2.4 停机成本模型iCsi对于医院诊疗设备来说，计算停机成本时仅考虑停机时间是不够的，其停机成本会随病患流动态变化。诊疗设备第 次预防性维护的停机损失成本为Csi=Tpiir+Cp（21）Tpiiir式中：为维护停机时间，包括故障维修用时和预防性维护用时；为第 个病患流区间内等效平均病人流量；为单位病人诊疗平均收益。若诊疗设备在服役周期内进行了 I 次预防性维护停机、K 次故障停机，则设备维护停机总成本第5期李岱高，等：基于麻雀搜索算法考虑病患流的医院诊疗设备维护优化研究517Cs为Cs=Ii=1Csi+Ki=1Cdi（22）2.5 诊疗设备总维护成本模型IJKCCriskCCs假设诊疗设备在服役周期内进行了 次维护停机，次维护调整，次故障停机，则其运行过程中的总成本 包括：故障风险成本、变动成本和停机损失成本，故考虑病患流分布的诊疗设备维护模型优化目标可表示为minC=Crisk+C+Css.t.0R R0（23）0R0tk通过求解优化模型获得设备确定合适的最高故障风险阈值和诊疗设备最低可靠度阈值，从而得出诊疗设备下一个最优维护时间。3 麻雀搜索算法 3.1 算法过程0R0从上文内容可看出，该优化模型涉及多个非线性函数，理论上难以直接求解，而麻雀搜索算法(sparrowsearchalgorithm,SSA)是受麻雀觅食行为和反捕食行为启发而提出的新型算法，与目标优化领域常用的遗传算法(geneticalgorithm,GA)、粒子群算法(particleswarmoptimization，PSO)相比，具有参数设置简单、寻优能力强和易于实现的优点19。因此，本文采用麻雀搜索算法求解该模型以获取故障风险阈值和诊疗设备可靠度阈值，算法具体步骤如下：P步骤 1确定种群数量，初始化麻雀种群X=x11x21.xd1x12x22.xd2.x1nx2n.xdnLubLlbmaxR2ST确定种群上界、种群下界和最大迭代次数，设定系统预警值和安全值；C(,R)f=C(,R)fiXbestXworstfgfw步骤 2将总维护成本作为每只麻雀个体的适应度函数，即，计算各麻雀个体的适应度值并进行排序，确定位置最优个体和位置最差个体，更新全局最佳适应度值和全局最差适应度值；步骤 3通过式(24)更新麻雀发现者，为整个种群寻找食物，并为加入者提供觅食的方向；Xt+1i,j=Xi,jexp(imax),R2n2Xt+1P+?Xi,jXt+1P?A+L,其他（25）SD步骤 5意识到危险的麻雀在种群中的比例为，位置由式(26)随机产生；Xt+1i,j=Xtbest+?Xti,jXtbest?,fi fgXi,j+K?Xti,jXtworst?(fi fw)+,fi=fg（26）maxfg步骤 6若当前迭代次数，转到步骤2，否则输出最优麻雀个体位置信息。3.2 基于麻雀搜索算法的维护策略模型求解nkk考虑病患流分布的诊疗设备维护是一种复杂的多变量优化问题，为将麻雀搜索算法应用于模型求解，需先在麻雀种群矩阵中合理布局变量的位置以便于模型求解。麻雀种群矩阵参数含义如图 4 所示，矩阵中每一行代表一个单独的麻雀个体。代表麻雀种群个体数，矩阵 列代表共有 维决策变量。将麻雀搜索算法应用于优化诊疗设备维护模型的具体步骤如下：步骤 1读取设置区间内模型所需的诊疗设备Xx11x21x31x22x32x2nxk1xk2xknx3nx12x1n0 取值R0 取值维护点 1第 k 维变量麻雀个体 1麻雀个体 2麻雀个体 n图 4 麻雀种群含义示意图Fig.4 Schematic diagram meaning of sparrow population518上海理工大学学报2023年第45卷所有历史数据与可读信息，并对所得数据进行预处理；NPDSDSTmaxLub，Llbnk步骤 2将参数初始化，输入麻雀种群规模、发现者比例、警戒者比例、安全值、最大迭代次数、种群上下限。初始化麻雀种群，生成 只 维麻雀种群矩阵 X。C(,R)fifgfwXbestXworst步骤 3以诊疗设备总维护成本最低作为目标函数，计算每只麻雀适应度值，排序获得最优适应度值、最差适应度值，及其对应位置和。步骤 4根据式(24)更新发现者位置，根据式(25)更新加入者位置，根据式(26)更新警戒者位置。fgfwXbestXworst步骤 5计算麻雀个体适应度值，并更新麻雀位置，重新确定，以及和。maxfgXbest步骤 6判断迭代次数是否达到最大次数，如果是则输出最优适应度值作为目标值，输出最优麻雀位置作为最优目标解；否则转到步骤4继续进行运算。算法求解流程如图 5 所示。开始输入诊疗设备服役数据,初始化麻雀种群及参数计算麻雀个体适应度值 fi,并通过排序得到 fg,在搜索区间内分别更新发现者位置、加入者位置、警戒者位置,计算更新后的个体适应度值并重新排序是否达到最大迭代次数?输出最优适应度值 fg 作为诊疗设备最低维护总成本,从相应麻雀个体位置 Xbest 中提取最优 0 取值和最优 R0 取值是否结束fw,及其位置 Xbest,Xworst图 5 SSA 算法流程图Fig.5 Flow chart of SSA algorithm 4 算例分析为验证上述模型的有效性，本文以某医院联影 UCT710 型号 CT 机的维护数据为例进行分析。根据假设 d，CT 机在基准病患流工况下的故障率函数为l0(t)=m(t)m1（27）m=2.5=500k=k3k+7k=12k+110k+1Cd=15000Cp=2000Cr=6000r=400Cy=800Tp=1 h根据联影 UCT710 型号 CT 历史失效数据分析拟合，得到威布尔分布形状参数，尺度参数，役龄回退因子、故障率递增因子。假 设 单 次 故 障 维 修 成 本，设备单次预防性维护成本，可靠度利用价值，单位病人诊疗收益，延迟维护基本惩罚成本，故障停机用时与预防性维护停机时间均取，在设备生命周期内随机生成初始维护点。U=事故风险u1,设备复杂性u2,设备可替代性u3,设备利用率u4V=高v1,较高v2,中v3,较低v4,低v5R=0.20.50.3000.20.40.20.10.10.10.60.20.100.10.50.20.20A=0.3,0.2,0.25,0.25B=0.15,0.505,0.23,0.095,0.02因 素 集，评 价 集，某 时 间 点 风 险 评 价 矩 阵，各 因 素 权 重，则由式(13)可得模糊综合评判结果，取其中隶属值最高的 0.505 作为故障风险值。4.1 算法求解比较首先，为了验证麻雀搜索算法在考虑病患流分布的医院诊疗设备预防性维护模型中的有效性和优越性，将其与普遍应用在寻优问题上的粒子群算法(PSO)和鲸鱼算法(WOA)进行对比，所有算法的共有参数保持一致，种群规模设置为 50，最大迭代次数设置为 1000。各算法的参数参照文献 20 设置，如表 1 所示。表中：W 为惯性因子；C1，C2为加速度常数；a，b 为系数。表 1 算法参数设置Tab.1 Algorithm parameter setting算法参数SSASTPDSD=0.2,=0.2,=0.8PSOWC1C2=0.729,=1.49445,=1.49445WOAab=2,=1比较结果如表 2 所示，可知 SSA 算法与 WOA算法求解结果相近，且明显优于 PSO 算法所得结果。各算法迭代如图 6 所示，可看出 PSO 算法收敛速度虽比 SSA 算法和 WOA 算法快，在迭代到50 次左右就达到了最优值，但所得最优值相较于其他两个算法偏差较大。由此可知，PSO 算法在第5期李岱高，等：基于麻雀搜索算法考虑病患流的医院诊疗设备维护优化研究519此处陷入了局部最优，SSA 算法和 WOA 算法由于较强的全局寻优能力，避免了陷入局部最优。在迭代次数相同且收敛精度相差无几的情况下，SSA 算法的收敛速度明显优于 WOA 算法。表 2 算法结果对比Tab.2 Comparison of algorithm results算法0R0CSSA0.85450.790634573.6WOA0.86380.802434591.9PSO0.87550.753335893.135 00036 00037 00038 00039 0002000600400800迭代次数1 000适应度值SSAWOAPSO图 6 各算法迭代图Fig.6 Algorithm iteration graph 4.2 结果分析0R0tk通过 SSA 算法寻优得到最佳故障风险阈值=0.8545，最佳设备可靠度阈值=0.7906，通过数值仿真求解得到诊疗设备最优维护时间 等参数如表 3 所示。表 3 SSA 算法求解结果Tab.3 Solution results of SSA algorithmtk/hCriskCSCC16474633.145297.578761.4918692.2考虑病患流分布情况下诊疗设备故障率演化状况与未考虑病患流分布（即基准病患流）情况下有一定差异。如图 7 所示，不考虑病患流分布情况下诊疗设备维护时间点延后且维护间隔增加，停机维护成本固然较低，但未考虑病患人流量因素会带来较高的故障风险成本，其区别由是否考虑病患流影响下初始维护方案的不同来反映。将考虑病患流分布记为方案 1，基准病患流分布情况记为方案 2，对比结果如表 4 所示。从表 4 中可看出，由于不考虑病患流的影响，方案 2 不产生变动成本，总停机维护成本较低。但忽略病患人流量的影响会埋下更多安全隐患，从而产生更高的故障风险成本，导致总维护成本高，达到 40649.8 元，高于考虑病患流影响下的总成本。将方案 2 中维护时间点 6391198 代入考虑病患流分布的诊疗设备维护模型中计算故障风险成本和停机维护成本，发现故障风险成本上升，产生惩罚成本 15286.5，总成本增加 13498.4。因此，方案 1 相较于方案 2 能更全面地考虑病患流的影响及其带来的故障风险。考虑病患流的影响在本文框架下有其合理性，建立诊疗设备可靠度和故障率模型时充分考虑病患流造成的影响，有利于更合理地制定设备维护策略。4.3 敏感度分析0R0rCr为探究考虑病患流分布的诊疗设备维护模型各参数和模型总维护成本之间的关系，需要分析模型对各参数的敏感度。本文采用只改变一个参数，其余参数固定不变的方法，来分析不同参数取值对总维护成本的影响。需分析参数有故障风险阈值、设备可靠度阈值、单位病人平均诊疗收益、可靠度利用价值。参数变化幅度为20%20%，具体的变化情况如表 5 所示。0.850.800.900.951.002000600400800时间/d1 0001 200可靠度基准病患流量下病患流影响下图 7 是否考虑病患流分布下设备可靠度变化Fig.7 Change of equipment reliability considering patientflow distribution表 4 是否考虑病患流影响下初始维护方案比较Tab.4 Comparison of initial maintenance schemes under the influence of patient flow方案维护时间点序列/h平均维护间隔/hai维护时刻病患流因子CriskCSCC155210844251.191.129751.34836.519985.834573.62639119844719627.223628.240649.8520上海理工大学学报2023年第45卷由表 5 结果分析可知：Crisk0R0rCr0CR0CriskR0Ca.故障风险成本对故障风险阈值和设备可靠度阈值较为敏感，对单位病人平均诊疗收益 和可靠度利用价值不敏感。当增大时，对故障风险值 的要求降低，制定诊疗设备维护策略时可接受更高的故障风险，从而带来更高的故障风险成本；同时，增大的故障风险阈值使维护过程中的维护变动更少，从而带来变动成本下降。当增大时，对设备可靠度要求更高，高可靠度使设备发生故障的风险降低，进而使故障风险成本减小；同时，高可靠度需要更多次的维护来保障，因此增大导致变动成本上升。CSr0R0Crr0R0aiCS0R0b.设备维护停机成本对单位病人平均诊疗收益 最敏感，对故障风险阈值和设备可靠度阈值较为敏感，对可靠度利用价值不敏感。增大直接影响设备停机维护时产生的病人诊疗费用损失，从而增加了停机维护成本。和的变动影响设备停机维护时间点和停机次数，但停机时病人诊疗费用损失受该时段病患流密度，即参数影响，因此，的变动趋势与和无正负相关。CCrrCrc.设备维护变动成本对可靠度利用价值较为敏感，对单位病人平均诊疗收益 不敏感。增大导致设备维护时间提前时产生更多的可靠度浪费，从而造成维护变动成本增加。C0R0CCrCrCd.设备维护总成本 对故障风险阈值、设备可靠度阈值最敏感，这两项参数直接影响诊疗设备预防性维护周期内停机维护与预防性维护次数；设备维护总成本 对可靠度利用价值较为敏感，因为在模型中参数仅影响维护时间提前时的维护变动成本，而维护变动成本仅为维护总成本的构成部分，因此对维护总成本的影响较小。5 结论针对诊疗设备这类故障危害大、事故影响坏、故障维修费用和停机损失高昂的特殊设备，将分别考虑安全性、经济性的故障风险方法和考虑设备连续运行状态的病患人流量因子纳入设备预防性维护模型，克服已有模型对于设备故障和设备连续运行对于寿命的影响认识不足的问题。研究结果表明，考量故障风险因素和病患人流量表 5 参数灵敏度分析Tab.5 Parameter sensitivity analysis参数变化幅度CriskCSCC总成本变动率/%4633.145297.578761.4918692.20020%2983.755836.4211280.7120100.887.5410%3168.444738.1410837.6418744.220.28+10%5966.835167.988461.3319596.144.84+20%6807.555311.698267.8520387.098.79R020%5883.424772.388340.7918996.591.6310%5559.775613.498557.1619730.425.55+10%3567.275536.8910030.7719134.932.37+20%3239.635447.2611470.0120156.907.84r20%4633.144237.068761.4917633.705.6610%4633.144768.728761.4918163.352.83+10%4633.145827.338761.4919221.962.83+20%4633.146259.098761.4919653.725.14Cr20%4633.145297.577584.4217515.136.2910%4633.145297.578411.0318341.741.87+10%4633.145297.579199.5619130.272.34+20%4633.145297.579800.4819731.195.56第5期李岱高，等：基于麻雀搜索算法考虑病患流的医院诊疗设备维护优化研究521因子使模型更贴合实际，根据模型制定的维护策略更加合理。本文所构建的基于病患流分布不均的预防性维护决策模型较为复杂、求解难度较高，故引入 SSA 算法工具辅助求解模型，算例分析结果表明，SSA 算法能有效提高模型求解速度与精确度。为诊疗设备制定合理有效的维护策略，对保证诊疗活动安全性和提高其使用效率都具有重大意义。本文提出的基于病患流分布不均的预防性维护决策模型能够帮助医院诊疗设备在运行期间更好地兼顾高可靠度和经济性问题，有效降低设备维护成本，减少设备停机对正常诊疗活动的影响，对诊疗设备维修具有借鉴作用。未来的研究可考虑多机系统下通过调度进一步降低维护停机成本，降低病人等待时间，使模型实用性更强。另外，鉴于医院诊疗设备是一个复杂的多部件系统，可进一步考虑系统内部件维修之间的关系，使模型更贴合实际。参考文献：韩浩.医疗设备预防性维护在医疗设备管理中的应用价值分析 J.中国设备工程,2021(7):4344.1CARMIGNANIG.Anintegratedstructuralframeworktocost-basedFMECA:thepriority-costFMECAJ.Reliability Engineering&System Safety,2009,94(4):861871.2ABBASGHOLIZADEHRAHIMIS,JAMSHIDIA,AIT-KADID,etal.PrioritizationoffailuresinradiationtherapydeliveryJ.IFAC-PapersOnLine,2016,49(12):1898 1903.3HERNNDEZ-LPEZLA,PIMENTEL-AGUILARAB,ORTIZ-POSADAS M R.An index to prioritize thepreventive maintenance of medical equipmentJ.HealthandTechnology,2020,10(2):399403.4张宁,刘勤明,叶春明,等.基于 NSGA-的医疗设备预防性维护与病人调度联合优化 J.重庆师范大学学报（自然科学版）,2021,38(1):130140.5沈南燕,武星,李静,等.自动化生产线中关键设备的预维护策略研究 J.机械工程学报,2020,56(21):231240.6XIA T B,XI L F,LEE J,et al.Optimal CBPM policy7consideringmaintenanceeffectsandenvironmentalconditionJ.The International Journal of AdvancedManufacturingTechnology,2011,56(9):11811193.LIDW,ZHANGZH,ZHONGQH,etal.Performancedeteriorationmodelingandoptimalpreventivemaintenancestrategy under schedu