CN103390091B - 一种传染病疫情优化控制方法 - Google Patents

一种传染病疫情优化控制方法 Download PDF

Info

Publication number
CN103390091B
CN103390091B CN201210139062.1A CN201210139062A CN103390091B CN 103390091 B CN103390091 B CN 103390091B CN 201210139062 A CN201210139062 A CN 201210139062A CN 103390091 B CN103390091 B CN 103390091B
Authority
CN
China
Prior art keywords
epidemic
cost
region
person
population
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201210139062.1A
Other languages
English (en)
Other versions
CN103390091A (zh
Inventor
黄顺祥
刘峰
李静
关彩虹
王新明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CHEMICAL DEFENSE COLLEGE OF PLA
Original Assignee
CHEMICAL DEFENSE COLLEGE OF PLA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CHEMICAL DEFENSE COLLEGE OF PLA filed Critical CHEMICAL DEFENSE COLLEGE OF PLA
Priority to CN201210139062.1A priority Critical patent/CN103390091B/zh
Publication of CN103390091A publication Critical patent/CN103390091A/zh
Application granted granted Critical
Publication of CN103390091B publication Critical patent/CN103390091B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Medical Treatment And Welfare Office Work (AREA)

Abstract

本发明属于流行病传播与控制领域,涉及一种传染病疫情优化控制方法。在现有技术中,决定疫情控制方案的控制参数没有进行优化控制,而在传染病疫情控制中,控制措施越严,相应的控制代价就越大,感染者人数也就会越小,如何确定控制参数的值,使得控制代价和感染者医治成本最小,这个优化控制问题急待解决。为此,本发明提出了一种传染病疫情优化控制方法,在该方法中,首先确定传染病疫情控制代价的目标函数,然后利用遗传算法优化求解了使最小传染病疫情控制代价目标函数值最小的控制参数值。该方法解决了使控制代价和感染者医治成本最小的优化控制问题,使得优化后控制代价远远小于优化前控制代价。

Description

一种传染病疫情优化控制方法
技术领域
本发明属于流行病传播与控制领域,涉及一种传染病疫情优化控制方法。
背景技术
2002年1月,严重急性呼吸综合征(SARS)在我国广东省首先爆发,并于2003年3月初输入北京,随之在北京迅速蔓延,同期也引起在世界范围内的爆发流行,其中,以我国广州、香港、北京以及台湾地区疫情尤为严重,据世界卫生组织估计,SARS对全球经济造成了约300亿美元的损失。2009年甲型H1N1流感疫情在全球的爆发,6月11日世界卫生组织宣布其预警级别为流感警戒最高级6级。甲型H1N1流感除了造成直接人员伤亡和巨额医疗费用外,对经济的间接影响、对民众心理和社会安定的危害都是非常严重的。即使普通的流感,每年就夺走了数万人的生命。传染病直接威胁人类生命和健康,在世界迅速全球化的今天,传染病正以惊人的速度蔓延,全球任何一个角落爆发的传染病都可能会引发全球的警觉。
当前的传染病模型,可分为动力学模型、统计模型和时空交互模型三类,这些模型的共同特点是参数通常是常数,为了表征传染病疫情发展跟人的控制措施有关,黄顺祥等人建立了非线性、变系数传染病预测模型,该模型引入了控制参数,使控制措施与传染病疫情发展联系起来。
在传染病疫情控制中,控制措施越严,相应的控制代价就越大,感染者人数也就会越小,怎样才能使控制代价和感染者医治成本等最小,是一个优化控制问题,目前还没有方法很好的解决这个问题。为此,本发明提出了一种传染病疫情优化控制方法,实现了传染病控制代价与控制效果的优化,为决策部门提供关键技术支持。
发明内容
为了解决传染病疫情优化控制,实现传染病疫情控制代价和感染者医治成本等最小的问题,本发明提出一种传染病疫情优化控制方法,该方法实现了传染病控制代价与控制效果的优化。
步骤一、确定传染病疫情控制代价的目标函数J的表达式:
J=AC(I,E,NI,NE,D,P,Ok)+BC(Dc,MI)+DCS,MI)(1)
其中,Ac(I,E,NI,NE,D,P,Ok)为疫情预防代价函数:
A c ( I , E , N I , N E , D , P , O k ) = Σ i = 1 N 1 λ I N I , i I i C I , i t I , i + Σ j = 1 N 2 λ E N E , j E j C E , j t E , j + Σ k = 1 N 3 r d , k D k S k t D , k + Σ l = 1 N 4 r p , l P l N S , l + Σ m = 1 N 5 O m - - - ( 2 )
I为一固定区域内发病者总人数,将该固定区域划分为N1个区域,Ii为第i个区域内发病者人数;NI为每个发病者在整个固定区域内平均亲密接触的人数,NI,i为第i个区域内每个发病者平均亲密接触的人数,1≤i≤N1,N1为大于等于1的正整数;CI,i为第i个区域内发病者亲密接触的隔离代价系数;λI为发病者的亲密接触者隔离率;tI,i为第i个区域的发病者亲密接触者的隔离时间;E为该固定区域的潜伏期者总人数,将该固定区域划分为N2个区域,Ej为第j个区域内潜伏期者人数,NE为每个潜伏期者在整个固定区域内平均亲密接触的人数,NE,j为第j个区域内每个潜伏期者平均亲密接触的人数,1≤j≤N2,N2为大于等于1的正整数;CE,j为潜伏期者亲密接触的隔离代价系数;λE为潜伏期者的亲密接触者隔离率;tE,j为第j个区域的潜伏期者亲密接触者的隔离时间;D表示整个固定区域单位面积的洗消成本,将该固定区域划分为N3个区域,Dk表示第k个区域单位面积的洗消成本,1≤k≤N3,N3为大于等于1的正整数;Sk为第k个区域被沾染面积;tD,k为第k个区域需要洗消的时间天数,设每天洗消一次;rd,k为第k个区域被沾染面积洗消的比例;P表示整个固定区域个人防护成本,Pl表示第l个区域个人防护的成本,1≤l≤N4,N4为大于等于1的正整数;NS,l为第l个区域的人数;rp,l为第l个区域人员防护的比例;Om指除前面4种行动代价之外的应急费用,比如控制流动人口付出的代价、疫苗研制、生产和使用费用,将该固定区域划分为N5个区域,1≤m≤N5,N5为大于等于1的正整数;
BC(Hc,MI)为发病者医治成本函数:
B C ( H c , M I ) = Σ r = 1 M M I , r H c , r - - - ( 3 )
其中,将整个固定的区域划分为M个区域,MI为到t时刻整个固定区域发病者总人数,MI,r为到t时刻第r个区域发病者总人数,1≤r≤M,M为大于等于1的正整数;Hc为整个区域每个感染者的平均治疗费用,Hc,r为第r个区域每个感染者的平均治疗费用;
DCs,MI)为死亡者社会代价函数:
D C ( η S , M I ) = Σ m = 1 N η S , m { int ( γ m M I , m ) + μ S [ int ( γ m M I , m ) + 1 ] α m } - - - ( 4 )
其中,ηS为整个固定区域的社会代价系统,ηS,m为第m个区域人身伤害标准,按保险公司计算方法来确定,μS为第m个区域社会代价系数,是对旅游业、交通运输业、商贸和餐饮业、会展业、文化娱乐业的综合影响指标,根据实际统计或评估数据反演,1≤m≤N,N为大于等于1的正整数;γm为第m个区域感染者平均死亡率;int(γmMI,m)为第m个区域死亡人数;αm为第m个区域死亡代价指数,其中0<αm<1;
目标函数J的表达式(1)中涉及的参数,除I、E外,其余都通过调查或反演的方法获得;根据考虑了防护和洗消措施的具有控制变量的传染病疫情控制模型,I、E是关于控制参数λEa,τEe,λIa,τIe,rda,τd,rpa,τp的函数,控制参数决定了控制方案;其中,λEa、λIa分别为对潜伏期者和发病期者的目标隔离率,τEe、τIe分别为对潜伏期者和发病期者的隔离措施生效时间,rda为目标洗消的污染面积占总污染面积的比例;rpa为该地区目标防护的人员占该地区需要防护人数,τd为洗消措施生效时间,τp为防护措施生效时间;
定义Ψ为控制变量,Ψ=[λEa,τEe,λIa,τIe,rda,τd,rpa,τp],所以变量E为控制变量Ψ的函数,变量I为控制变量Ψ的函数,即目标函数J=J(Ψ),是Ψ的函数。
步骤二、采用遗传算法对Ψ中的各控制参数进行优化求解;
在优化求解过程中,遗传算法采用的适应度函数与目标函数J成负相关的关系,且适应度函数值恒不为负;优化求解得到的Ψ中的各控制参数组成了优化的疫情控制方案。
有益效果:
本发明提出的一种传染病疫情优化控制方法,解决了使控制代价和疫情损失的优化问题,使得优化后控制代价远远小于优化前控制代价且感染者总人数减小,因此既可以大量节省疫情控制成本而取得良好的经济效益,又可通过减少感染人数而取得良好的社会效益。
附图说明
图1为本发明的优化控制前后传染病累积发病人数统计对比图。
具体实施方式
下面结合具体实施例对本发明做进一步说明。
一种传染病疫情优化控制方法,其实现的具体实施步骤如下:
步骤一、确定传染病疫情控制代价的目标函数J的表达式:
J=AC(I,E,NI,NE,D,P,Ok)+BC(Hc,MI)+DCS,MI)(1)
其中,Ac(I,E,NI,NE,D,P,Ok)为疫情预防代价函数:
A c ( I , E , N I , N E , D , P , O k ) = Σ i = 1 N 1 λ I N I , i I i C I , i t I , i + Σ j = 1 N 2 λ E N E , j E j C E , j t E , j + Σ k = 1 N 3 r d , k D k S k t D , k + Σ l = 1 N 4 r p , l P l N S , l + Σ m = 1 N 5 O m - - - ( 2 )
I为一固定区域内发病者总人数,将该固定区域划分为N1个区域,Ii为第i个区域内发病者人数;NI为每个发病者在整个固定区域内平均亲密接触的人数,NI,i为第i个区域内每个发病者平均亲密接触的人数,1≤i≤N1,N1为大于等于1的正整数;CI,i为第i个区域内发病者亲密接触的隔离代价系数;λI为发病者的亲密接触者隔离率;tI,i为第i个区域的发病者亲密接触者的隔离时间;E为该固定区域的潜伏期者总人数,将该固定区域划分为N2个区域,Ej为第j个区域内潜伏期者人数,NE为每个潜伏期者在整个固定区域内平均亲密接触的人数,NE,j为第j个区域内每个潜伏期者平均亲密接触的人数,1≤j≤N2,N2为大于等于1的正整数;CE,j为潜伏期者亲密接触的隔离代价系数;λE为潜伏期者的亲密接触者隔离率;tE,j为第j个区域的潜伏期者亲密接触者的隔离时间;D表示整个固定区域单位面积的洗消成本,将该固定区域划分为N3个区域,Dk表示第k个区域单位面积的洗消成本,1≤k≤N3,N3为大于等于1的正整数;Sk为第k个区域被沾染面积;tD,k为第k个区域需要洗消的时间天数,设每天洗消一次;rd,k为第k个区域被沾染面积洗消的比例;P表示整个固定区域个人防护成本,Pl表示第l个区域个人防护的成本,1≤l≤N4,N4为大于等于1的正整数;NS,l为第l个区域的人数;rp,l为第l个区域人员防护的比例;Om指除前面4种行动代价之外的应急费用,比如控制流动人口付出的代价、疫苗研制、生产和使用费用,将该固定区域划分为N5个区域,1≤m≤N5,N5为大于等于1的正整数。
BC(Hc,MI)为发病者医治成本函数:
B C ( H c , M I ) = Σ r = 1 M M I , r H c , r - - - ( 3 )
其中,将整个固定的区域划分为M个区域,MI为到t时刻整个固定区域发病者总人数,MI,r为到t时刻第r个区域发病者总人数,1≤r≤M,M为大于等于1的正整数;Hc为整个区域每个感染者的平均治疗费用,Hc,r为第r个区域每个感染者的平均治疗费用;
DCs,MI)为死亡者社会代价函数:
D C ( η S , M I ) = Σ m = 1 N η S , m { int ( γ m M I , m ) + μ S [ int ( γ m M I , m ) + 1 ] α m } - - - ( 4 )
其中,ηS为整个固定区域的社会代价系统,ηS,m为第m个区域人身伤害标准,按保险公司计算方法来确定,μS为第m个区域社会代价系数,是对旅游业、交通运输业、商贸和餐饮业、会展业、文化娱乐业的综合影响指标,根据实际统计或评估数据反演,1≤m≤N,N为大于等于1的正整数;γm为第m个区域感染者平均死亡率;int(γmMI,m)为第m个区域死亡人数;αm为第m个区域死亡代价指数,其中0<αm<1。
目标函数J的表达式(1)中涉及的参数,除I、E外,其余通过调查或反演的方法获得,对于I、E由传染病模型的状态方程确定,该状态方程在2012年5月7日申请的名为“一种具有控制变量的传染病疫情预测方法”(申请号:201210137662.4)的专利中公开。
下面根据传染病模型的状态方程求解I、E的值:
设进入该地区的人都为易感者,流出的人中易感者、潜伏期者、感染期者和治愈者在流出的人中所占的比例与该地区易感者、潜伏期者、感染期者和治愈者在该地区总人口中所占的比例成正比关系。令S表示在该地区易感者的总人数,E表示在该地区潜伏期者的总人数,I表示在该地区发病者的总人数;R表示在该地区移出者的总人数,其中移出者包括治愈者和死亡者,由此根据动力学方程得到传染病模型的状态方程如公式(5)、(6)、(7)、(8)所示:
dS dt = - k ( 1 - γ ) ( 1 - λ I ) SI - k ( 1 - γ ) ( 1 - λ E ) SE - αS + α - - - ( 5 )
dE dt = k ( 1 - γ ) ( 1 - λ I ) SI + k ( 1 - γ ) ( 1 - λ E ) SE - E D 1 - αE - - - ( 6 )
dI dt = E D 1 - I D 2 - - - ( 7 )
dR dt = I D 2 - αR - - - ( 8 )
其中,κ为传染系数,κ0为基本传染率,是人口密度的函数,κ1为季节影响系数,为周期,天;λE、λI分别为针对潜伏期者和发病期者的隔离率,
即将潜伏期者和发病期者在该地区内被有效隔离的比例;γ表示感染下降率;D1表示潜伏期;D2表示传染期;α为流动人口人数,具体指疫区内外的相互流动比例,其中假设发病者不流出;
传染病模型的状态方程中,潜伏期者和发病期者的隔离率λE、λI分别通过以下两个函数确定:
λ E ( t ) = λ Ea - λ E 0 π { π 2 + arctan [ v E ( t - τ Ee ) ] } + λ E 0 λ I ( t ) = λ Ia - λ I 0 π { π 2 + arctan [ v I ( t - τ Ie ) ] } + λ I 0 - - - ( 9 )
其中λEa、λIa分别为对潜伏期者和发病期者的目标隔离率,λE0、λI0分别为对潜伏期者和发病期者的初始隔离率,νE、νI分别为对潜伏期者和发病期者的隔离效率因子,λEaE0、λIaI0分别为控制方案对应的对潜伏期者和发病期者的隔离措施增加强度,τEe、τIe分别为对潜伏期者和发病期者的隔离措施生效时间,即λE(t)和λI(t)函数曲线拐点处对应的时间;
当采取防护和洗消措施,感染下降率将改变,感染下降率γ的取值由公式(10)、(11)、(12)确定:
γ(t)=μdrd(t)+μprp(t)-μdμprd(t)rp(t)(10)
r d ( t ) = r da - r d 0 π { π 2 + arctan [ v d ( t - τ d ) ] } + r d 0 - - - ( 11 )
r p ( t ) = r pa - r p 0 π { π 2 + arctan [ v p ( t - τ p ) ] } + r p 0 - - - ( 12 )
其中,μd为洗消后被消灭病毒占洗消前病毒总数的比例,μp为防护后被隔离的病毒占防护前病毒总数的比例,rd(t)为对被洗消的病毒区占总病毒区的比例,rp(t)为该地区被防护的人员占需要防护人员总人数的比率。在公式(7)中μdrd(t)为不采取防护措施时,病毒区中被洗消病毒占洗消前病毒的比例;μprp(t)为不采取洗消措施时,该地区被防护人员被隔离病毒占防护前病毒总数的比例;μdμprd(t)rp(t)指同时采取洗消和防护措施时,重复计算的同时被洗消和隔离的病毒占病毒总数的比例。
rda为目标洗消的污染面积占总污染面积的比例;rd0为初始洗消的污染面积占总污染面积的比例;νd为洗消效率因子;rda-rd0为洗消措施增加强度。τd为洗消措施生效时间,即rd(t)函数曲线拐点处对应的时间;rpa为该地区目标防护的人员占该地区需要防护人数,比例;rp0为初始防护比例;νp为防护效率因子;rpa-rp0为防护措施增加强度;τp为防护措施生效时间,即rp(t)函数曲线拐点处对应的时间。
定义Ψ为控制变量,Ψ=[λEa,τEe,λIa,τIe,rda,τd,rpa,τp],所以变量E为控制变量Ψ的函数,变量I为控制变量Ψ的函数,因此目标函数J=J(Ψ)。
下面步骤二~步骤九采用遗传算法对Ψ中的各控制参数进行优化求解;在优化求解过程中,遗传算法采用的适应度函数与目标函数J成负相关的关系,且适应度函数值恒不为负;优化求解得到的Ψ中的各控制参数组成了优化的疫情控制方案。具体来说:
步骤二、采取多参数级联编码方法对8个控制参数进行二进制编码,编码过程中将潜伏 期者隔离措施生效时间τ Ee 和发病者隔离措施生效时间τ Ie 、洗消措施生效时间τ d 、防护措施 生效时间τ p 精度有效数字位数设置为小数点后1位;潜伏期者发病者目标洗消比例τ da 、目标 防护比例r pa 精度有效数字位数设置为小数点后2位;潜伏期者目标隔离率λ Ea 和发病者目标 隔离率λ Ia 精度有效数字位数设置为小数点后4位。
步骤三、随机生成数目适当的初始种群X(0)=(x1,x2,…xk),x1、x2,…xk均为初始种群的个体,每个个体代表Ψ中的8个控制参数,k为初始种群个体的数目。
步骤四、针对步骤三生成的初始种群X(0)的个体xi通过公式(1)、(2)、(3)、(4)得到初始种群各个体的传染病疫情控制代价目标函数J(X(0,xi));其中,X(0,xi)表示第0代种群中的其中一个个体;J(X(0,xi))表示第0代种群中的其中一个个体对应的疫情控制代价。
判断J(X(0,xi))是否有个体满足使发病者模拟累计值中的最大值小于统计累计值,如果满足有,输出满足条件的种群,并从中选取出最小传染病疫情控制代价目标函数值J的个体,该个体对应的控制参数及J值即为所求,退出本算法;如果不满足终止条件,继续下面的步骤。
其中,发病者模拟累计值指的是当前第0代种群中发病者总的人数,根据Ψ的取值带入背景技术中提及的具有控制变量的传染病预测控制模型计算得到,统计累计值为实际情况中的发病人数。
定义第T代种群的个体xi得到相应措施的传染病疫情控制代价目标函数J(X(T,xi)),X(T,xi)表示第T代种群的个体xi,i的取值范围为从1到第T代种群的个体数目总数。
步骤五、用Fit(T,xi)表示适应度函数,代表传染病疫情控制代价的大小程度,将是对控制变量进行优化取舍的基础,即:
Fit(T,xi)=(1/J(X(T,xi)))m(13)
其中m≥1,X(T,xi)表示第t代种群的第T代种群的个体xi
步骤六、初始时,令T=0。应用选择算子产生中间代Xr(T)(T为迭代次数):采用轮盘赌方法作为选择算子,即根据适应度函数值大小对X(T)中的个体进行随机选择,X(T)中的个体被选中的概率正比于其适应度函数值Fit(T,xi),经过此轮盘赌方法选中的个体组成中中间代Xr(T);同时也将精英保留策略作为选择算子,保留X(T)中使Fit(T,xi)取最大值的个体。
步骤七、对Xr(T)应用两点交叉算子:按概率pc对Xr(T)中的个体进行选取,将被选中的每两个个体的基因链进行交叉,生成两个新的个体,这两个新个体与其余未被选择的个体一起组成Xe(T),其中交叉位置是随机的,pc为人为设定。
步骤八、对Xe(T)应用变异算子,得到X(T+1),具体是:Xe(T)随机选取新个体,将新个体的基因链的所有基因按概率pm进行变异,得到的所有新个体除去使Fit(T,xi)取最小值的个体,并且与步骤五中采用精英保留策略保留的个体一起组成种群X(T+1),其中pm为人为设定。
步骤九、设定X(T+1)的种群中所有个体生成的发病者模拟累计值中的最大值小于统计累计值为终止条件。判断X(T+1)的种群是否满足终止条件,如果满足终止条件,输出X(T+1)种群的个体,并从中选取出最小传染病疫情控制代价目标函数值J的个体,该个体对应的控制参数及J值即为所求;如果不满足终止条件,T递增1,返回步骤七,开始下一次迭代,直到有满足终止条件的种群为止。
其中,发病者模拟累计值指的是当前第T+1代种群中发病者总的人数,根据控制参数对疫情发展情况的影响计算得到,统计累计值为实际情况中的发病人数。用于以后的传染病疫情优化控制时,该统计累计值为国家控制传染病疫情要达到的目标发病人数。
基于上述具体的传染病疫情优化控制方法,在本实施例中,以北京2003年SARS疫情控制为例,将北京地区看作整体疫情控制区域,则建立目标函数如下:
J = λ I N I IC I t I + λ E N E EC E t E + ∫ 0 t r d DS d dt + r p PN S + M I H c + η S { int ( γM I ) + μ S [ int ( γM I ) + 1 ] α } - - - ( 14 )
其中,I为发病者;NI为每个发病者平均亲密接触的人数,设NI=6;CI为发病者亲密接触的隔离代价系数,设CI=90元/天;λI为发病者的亲密接触者隔离率;tI为与发病者亲密接触者的隔离时间;tI=15天;E为潜伏期者人数,NE为每个潜伏期者平均亲密接触的人数,设NE=6;CE为潜伏期者亲密接触的隔离代价系数,设CE=90元/天;λE为潜伏期者的亲密接触者隔离率;tE为潜伏期者亲密接触者的隔离时间;tE=15天;D表示单位面积的洗消成本,取D=0.2元/平方米;SD为被沾染面积,取SD=497平方公里;rd为被沾染面积洗消的比例,取P表示个人防护的成本,2元/人.天;Ns为总人数;rp为人员防护的比例;Hc为每个感染者的平均治疗费用,取Hc=20000元;ηS为人身伤害标准,μS为社会代价系数,根据《中国经济快讯》文献(文钊,2003)和《21世纪经济》(邹光祥,2003)报道的数据通过反演可估算μS=451.4,北京市2003年取ηS=480000元/人;χ为感染者平均死亡率,γ=0.06;α为死亡代价指数,取α=0.9。
优化前后对比:优化前控制代价为2.60×1010元,优化后控制代价为8.49×109元,
表1优化控制前后控制参数值变化对比
参数 τIe λIa τEe λea rda τd rpa τp
优化前 37.0 0.9966 43.5 0.9024 0.80 27.1 0.48 32.4
优化后 53.7 0.5939 10.4 0.9973 0.97 14.5 0.10 26.6
表2优化控制前后各代价对比

Claims (5)

1.一种传染病疫情优化控制方法,其特征在于,该方法实现的步骤如下:
步骤一、确定传染病疫情控制代价的目标函数J的表达式:
J=AC(I,E,NI,NE,D,P,Ok)+BC(Hc,MI)+DCS,MI)(1)
其中,Ac(I,E,NI,NE,D,P,Ok)为疫情预防代价函数:
A c ( I , E , N I , N E , D , P , O k ) = Σ i = 1 N 1 λ I N I , i I i C I , i t I , i + Σ j = 1 N 2 λ E N E , j E j C E , j t E , j + Σ k = 1 N 3 r d , k D k S k t D , k + Σ l = 1 N 4 r p , l P l N S , l + Σ m = 1 N 5 O m
I为一固定区域内发病者总人数,将该固定区域划分为N1个区域,Ii为第i个区域内发病者人数;NI为每个发病者在整个固定区域内平均亲密接触的人数,NI,i为第i个区域内每个发病者平均亲密接触的人数,1≤i≤N1,N1为大于等于1的正整数;CI,i为第i个区域内发病者亲密接触的隔离代价系数;λI为发病者的亲密接触者隔离率;tI,i为第i个区域的发病者亲密接触者的隔离时间;E为该固定区域的潜伏期者总人数,将该固定区域划分为N2个区域,Ej为第j个区域内潜伏期者人数,NE为每个潜伏期者在整个固定区域内平均亲密接触的人数,NE,j为第j个区域内每个潜伏期者平均亲密接触的人数,1≤j≤N2,N2为大于等于1的正整数;CE,j为潜伏期者亲密接触的隔离代价系数;λE为潜伏期者的亲密接触者隔离率;tE,j为第j个区域的潜伏期者亲密接触者的隔离时间;D表示整个固定区域单位面积的洗消成本,将该固定区域划分为N3个区域,Dk表示第k个区域单位面积的洗消成本,1≤k≤N3,N3为大于等于1的正整数;Sk为第k个区域被沾染面积;tD,k为第k个区域需要洗消的时间天数,设每天洗消一次;rd,k为第k个区域被沾染面积洗消的比例;P表示整个固定区域个人防护成本,Pl表示第l个区域个人防护的成本,1≤l≤N4,N4为大于等于1的正整数;NS,l为第l个区域的人数;rp,l为第l个区域人员防护的比例;Om指除前面4种行动代价之外的应急费用,比如控制流动人口付出的代价、疫苗研制、生产和使用费用,将该固定区域划分为N5个区域,1≤m≤N5,N5为大于等于1的正整数;
BC(Hc,MI)为发病者医治成本函数:
B C ( H c , M I ) = Σ r = 1 M M I , r H c , r - - - ( 3 )
其中,将整个固定的区域划分为M个区域,MI为到t时刻整个固定区域发病者总人数,MI,r为到t时刻第r个区域发病者总人数,1≤r≤M,M为大于等于1的正整数;Hc为整个区域每个感染者的平均治疗费用,Hc,r为第r个区域每个感染者的平均治疗费用;
DCs,MI)为死亡者社会代价函数:
其中,ηS为整个固定区域的社会代价系统,ηS,m为第m个区域人身伤害标准,按保险公司计算方法来确定,μS为第m个区域社会代价系数,是对旅游业、交通运输业、商贸和餐饮业、会展业、文化娱乐业的综合影响指标,根据实际统计或评估数据反演,1≤m≤N,N为大于等于1的正整数;γm为第m个区域感染者平均死亡率;int(γmMI,m)为第m个区域死亡人数;αm为第m个区域死亡代价指数,其中0<αm<1;
目标函数J的表达式(1)中涉及的参数,除I、E外,其余都通过调查或反演的方法获得;根据考虑了防护和洗消措施的具有控制变量的传染病疫情控制模型,I、E是关于控制参数λEaEeIaIe,rdad,rpap的函数,控制参数决定了控制方案;其中,λEa、λIa分别为对潜伏期者和发病期者的目标隔离率,τEe、τIe分别为对潜伏期者和发病期者的隔离措施生效时间,rda为目标洗消的污染面积占总污染面积的比例;rpa为该地区目标防护的人员占该地区需要防护人数,τd为洗消措施生效时间,τp为防护措施生效时间;
定义Ψ为控制变量,Ψ=[λEaEeIaIe,rdad,rpap],所以变量E为控制变量Ψ的函数,变量I为控制变量Ψ的函数,即目标函数J=J(Ψ),是Ψ的函数;
下面步骤二~步骤九采用遗传算法对Ψ中的各控制参数进行优化求解:
步骤二、采取多参数级联编码方法对8个控制参数进行二进制编码,编码过程中将潜伏期者隔离措施生效时间τEe和发病者隔离措施生效时间τIe、洗消措施生效时间τd、防护措施生效时间τp精度有效数字位数设置为小数点后1位;潜伏期者发病者目标洗消比例τda、目标防护比例rpa精度有效数字位数设置为小数点后2位;潜伏期者目标隔离率λEa和发病者目标隔离率λIa精度有效数字位数设置为小数点后4位;
步骤三、随机生成数目适当的初始种群X(0)=(x1,x2,…xk),x1、x2,…xk均为初始种群的个体,每个个体代表Ψ中的8个控制参数,k为初始种群个体的数目;
步骤四、针对步骤三生成的初始种群X(0)的个体xi通过公式(1)、(2)、(3)、(4)得到初始种群各个体的传染病疫情控制代价目标函数J(X(0,xi));其中,X(0,xi)表示第0代种群中的其中一个个体;J(X(0,xi))表示第0代种群中的其中一个个体对应的疫情控制代价;
判断J(X(0,xi))是否有个体满足使发病者模拟累计值中的最大值小于统计累计值,如果满足有,输出满足条件的种群,并从中选取出最小传染病疫情控制代价目标函数值J的个体,该个体对应的控制参数及J值即为所求,退出本算法;如果不满足终止条件,继续下面的步骤;
其中,发病者模拟累计值指的是当前第0代种群中发病者总的人数,根据Ψ的取值带入背景技术中提及的具有控制变量的传染病预测控制模型计算得到,统计累计值为实际情况中的发病人数;
定义第T代种群的个体xi得到相应措施的传染病疫情控制代价目标函数J(X(T,xi)),X(T,xi)表示第T代种群的个体xi,i的取值范围为从1到第T代种群的个体数目总数;
步骤五、用Fit(T,xi)表示适应度函数,代表传染病疫情控制代价的大小程度,将是对控制变量进行优化取舍的基础,即:
Fit(T,xi)=(1/J(X(T,xi)))m
其中m≥1,X(T,xi)表示第t代种群的第T代种群的个体xi
步骤六、初始时,令T=0;应用选择算子产生中间代Xr(T),其中T为迭代次数:采用轮盘赌方法作为选择算子,即根据适应度函数值大小对X(T)中的个体进行随机选择,X(T)中的个体被选中的概率正比于其适应度函数值Fit(T,xi),经过此轮盘赌方法选中的个体组成中中间代Xr(T);同时也将精英保留策略作为选择算子,保留X(T)中使Fit(T,xi)取最大值的个体;
步骤七、对Xr(T)应用两点交叉算子:按概率pc对Xr(T)中的个体进行选取,将被选中的每两个个体的基因链进行交叉,生成两个新的个体,这两个新个体与其余未被选择的个体一起组成Xe(T),其中交叉位置是随机的,pc为人为设定;
步骤八、对Xe(T)应用变异算子,得到X(T+1),具体是:Xe(T)随机选取新个体,将新个体的基因链的所有基因按概率pm进行变异,得到的所有新个体除去使Fit(T,xi)取最小值的个体,并且与步骤五中采用精英保留策略保留的个体一起组成种群X(T+1),其中pm为人为设定;
步骤九、设定X(T+1)的种群中所有个体生成的发病者模拟累计值中的最大值小于统计累计值为终止条件;判断X(T+1)的种群是否满足终止条件,如果满足终止条件,输出X(T+1)种群的个体,并从中选取出最小传染病疫情控制代价目标函数值J的个体,该个体对应的控制参数及J值即为所求;如果不满足终止条件,T递增1,返回步骤七,开始下一次迭代,直到有满足终止条件的种群为止;
其中,发病者模拟累计值指的是当前第T+1代种群中发病者总的人数,根据控制参数对疫情发展情况的影响计算得到,统计累计值为实际情况中的发病人数,用于以后的传染病疫情优化控制时,该统计累计值为国家控制传染病疫情要达到的目标发病人数。
2.如权利要求1所述的一种基于遗传算法的传染病疫情优化控制方法,其特征在于,将优化求解得到的控制变量Ψ代入具有控制变量的传染病疫情预测模型中,得到I、E,进而求得目标函数J的值,即优化的疫情控制方案对应的疫情控制代价。
3.如权利要求1所述的一种传染病疫情优化控制方法,其特征在于,遗传算法的适应度函数一种优选形式为取目标函数倒数的Z次方,其中Z≥1。
4.如权利要求1所述的一种传染病疫情优化控制方法,其特征在于,采用遗传算法优化参数时,选择算子采用轮盘赌方法和精英保留策略。
5.如权利要求1所述的一种传染病疫情优化控制方法,其特征在于,在遗传算法中,采取多参数级联编码方法对8个控制参数进行二进制编码,编码过程中将潜伏期者隔离措施生效时间τEe和发病者隔离措施生效时间τIe、洗消措施生效时间τd、防护措施生效时间τp精度有效数字位数设置为小数点后1位;潜伏期者发病者目标洗消比例τda、目标防护比例rpa精度有效数字位数设置为小数点后2位;潜伏期者目标隔离率λEa和发病者目标隔离率λIa精度有效数字位数设置为小数点后4位。
CN201210139062.1A 2012-05-08 2012-05-08 一种传染病疫情优化控制方法 Expired - Fee Related CN103390091B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210139062.1A CN103390091B (zh) 2012-05-08 2012-05-08 一种传染病疫情优化控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210139062.1A CN103390091B (zh) 2012-05-08 2012-05-08 一种传染病疫情优化控制方法

Publications (2)

Publication Number Publication Date
CN103390091A CN103390091A (zh) 2013-11-13
CN103390091B true CN103390091B (zh) 2016-04-06

Family

ID=49534360

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210139062.1A Expired - Fee Related CN103390091B (zh) 2012-05-08 2012-05-08 一种传染病疫情优化控制方法

Country Status (1)

Country Link
CN (1) CN103390091B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108597617B (zh) * 2018-04-11 2022-05-20 平安科技(深圳)有限公司 流行病分级预测方法及装置、计算机装置和可读存储介质
CN109192318A (zh) * 2018-07-11 2019-01-11 辽宁石油化工大学 描述传染病传播过程的简化SIS模型的建立与Laplace分析
CN112037465B (zh) * 2020-03-02 2022-02-22 电子科技大学 疫情防控方法、装置、服务器及存储介质
CN110993119B (zh) * 2020-03-04 2020-07-07 同盾控股有限公司 基于人口迁移的疫情预测方法、装置、电子设备及介质
CN111341463B (zh) * 2020-03-17 2023-08-15 腾讯科技(深圳)有限公司 疫情排查方法、装置、计算机设备和存储介质
CN114864104A (zh) * 2020-03-26 2022-08-05 北京懿医云科技有限公司 传染病疫情预测方法及装置、存储介质、电子设备
CN111370138B (zh) * 2020-04-01 2023-04-11 大连理工大学 一种面向宫颈癌防治的hpv长期最优治疗策略制定方法
CN111540478B (zh) * 2020-04-22 2023-10-10 第四范式(北京)技术有限公司 疫情推演仿真系统和仿真方法
CN111627563B (zh) * 2020-05-22 2023-04-11 种鹏云 一种交通运输对covid-19传播影响的评估方法
CN111667920B (zh) * 2020-06-04 2023-05-02 中网联金乐盟科技(北京)有限公司 一种基于定位数据的人群密切接触排查方法
CN111768875B (zh) * 2020-06-23 2023-08-01 深圳市城市公共安全技术研究院有限公司 传染病疫情预测方法、系统、装置及存储介质
CN111739658B (zh) * 2020-07-20 2021-02-02 医渡云(北京)技术有限公司 基于输入病例预测传染病趋势的方法及装置
CN112652403B (zh) * 2020-12-25 2023-07-14 中国科学技术大学 疫情预测方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101777092A (zh) * 2009-12-18 2010-07-14 中国人民解放军防化指挥工程学院 基于考虑空间非均匀性的传染病传播模型疫情预测方法
CN101794342A (zh) * 2009-09-30 2010-08-04 中国人民解放军防化指挥工程学院 基于非线性、变系数预测模型的传染病疫情预测分析方法
CN102222163A (zh) * 2011-05-25 2011-10-19 中国人民解放军防化指挥工程学院 一种具有免疫期且受季节影响的传染病疫情预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101794342A (zh) * 2009-09-30 2010-08-04 中国人民解放军防化指挥工程学院 基于非线性、变系数预测模型的传染病疫情预测分析方法
CN101777092A (zh) * 2009-12-18 2010-07-14 中国人民解放军防化指挥工程学院 基于考虑空间非均匀性的传染病传播模型疫情预测方法
CN102222163A (zh) * 2011-05-25 2011-10-19 中国人民解放军防化指挥工程学院 一种具有免疫期且受季节影响的传染病疫情预测方法

Also Published As

Publication number Publication date
CN103390091A (zh) 2013-11-13

Similar Documents

Publication Publication Date Title
CN103390091B (zh) 一种传染病疫情优化控制方法
CN103390089B (zh) 一种具有控制变量的传染病疫情预测方法
CN101794342B (zh) 一种基于考虑隔离措施的传染病疫情预测方法
Shen Logistic growth modelling of COVID-19 proliferation in China and its international implications
CN102222163B (zh) 一种具有免疫期且受季节影响的传染病疫情预测方法
González-Parra et al. Modeling the epidemic waves of AH1N1/09 influenza around the world
CN106874238A (zh) 一种威布尔型单元备件需求量的计算方法
Farthing et al. Assessing the efficacy of interventions to control indoor SARS-Cov-2 transmission: An agent-based modeling approach
Rustan et al. The Outbreak’S modeling of coronavirus (Covid-19) using the modified seir model in Indonesia
Hu et al. Optimization of COVID-19 prevention and control with low building energy consumption
Chen et al. Cost-effectiveness of influenza control measures: a dynamic transmission model-based analysis
Lin et al. Transmission dynamics of cholera with hyperinfectious and hypoinfectious vibrios: mathematical modelling and control strategies
CN103793619A (zh) 一种新型传染病空间扩散模拟方法
Yang et al. Modeling and analysis of COVID-19 based on a time delay dynamic model
Luo et al. Analysis of potential risk of COVID-19 infections in China based on a pairwise epidemic model
Sheng et al. The modeling and analysis of the COVID-19 pandemic with vaccination and isolation: a case study of Italy
Özaltın et al. Optimal distribution of the influenza vaccine
CURIONI UK domestic heat demand modelling
Ruiz Estrada et al. The Application of the 2019-nCoV Global Economic Impact Simulator (the 2019-nCoV-GEI-Simulator) in China
Snacken Control of influenza. Public health policies
Goswami et al. A mathematical model of influenza: stability and treatment
Liu et al. Epidemic trends analysis and risk estimation of 2019-nCoV outbreak
CN112420209B (zh) 一种评估活禽市场关闭措施效果数学模型的建立方法
Liu et al. A K-Means based Model towards Ebola Virus Prorogation Prediction
CN103093116B (zh) 一种电网供需风险评估方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160406