CN106997509A - 一种不确定信息融合的应急物资分布需求预测方法 - Google Patents
一种不确定信息融合的应急物资分布需求预测方法 Download PDFInfo
- Publication number
- CN106997509A CN106997509A CN201710193250.5A CN201710193250A CN106997509A CN 106997509 A CN106997509 A CN 106997509A CN 201710193250 A CN201710193250 A CN 201710193250A CN 106997509 A CN106997509 A CN 106997509A
- Authority
- CN
- China
- Prior art keywords
- condition
- disaster
- attribute
- value
- rule
- 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.)
- Granted
Links
- 239000000463 material Substances 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 53
- 230000004927 fusion Effects 0.000 title claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims abstract description 4
- 239000000284 extract Substances 0.000 claims abstract description 4
- 238000000205 computational method Methods 0.000 claims abstract description 3
- 230000008569 process Effects 0.000 claims description 7
- 230000004069 differentiation Effects 0.000 claims description 5
- 238000011524 similarity measure Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 3
- 238000005303 weighing Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 14
- 238000007726 management method Methods 0.000 description 8
- 239000011449 brick Substances 0.000 description 6
- 238000013277 forecasting method Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 238000007418 data mining Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000003064 k means clustering Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2413—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on distances to training or reference patterns
- G06F18/24133—Distances to prototypes
- G06F18/24137—Distances to cluster centroïds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Data Mining & Analysis (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Game Theory and Decision Science (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Operations Research (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种不确定信息融合的应急物资分布需求预测方法,首先分析历史灾情数据库中各范例的基本属性,提取范例的共同特征并运用粗糙集方法进行简约,根据各条件属性对物资预测结果的影响程度计算出属性权重值;然后将简约后的灾情数据运用反馈补偿的K‑means方法将案例分类,并根据当前案例与各类中心点之间的距离判断当前案例属于哪一类,根据范例相似度的计算方法,检索出相似度最高的历史目标范例;最后根据目标范例的物资消耗以及范例属性权重值来线性推测当前范例的物资需求量。此外,根据简约后的范例还可得出当前案例的救援优先级。通过该方法,解决了救援物资需求分布预测问题,为航空应急救援指挥系统的调度决策提供依据。
Description
技术领域
本发明涉及一种应急物质需求预测技术,尤其涉及救援信息不确定性下应急事件的多源信息融合方式和信息规则挖掘方法,为准确的实施应急物质调度和救援提供依据。
背景技术
在抗灾救援及处置突发事件的各项措施中,获取准确的救援信息特别是救援物资需求信息,是应急协同调度决策的先决条件。目前,应急救援物资需求预测主要从以下几个方面开展研究:⑴应急物资需求的专家经验预测法,主要利用历史震灾分析的数据和相关理论,根据统计的震级、死亡人数、受伤人数、倒塌房屋等多种参数,提出了需求预测模型。⑵应急救援物资需求的事例推理法,首先描述和提取发生或即将发生的灾害特征,根据这些特征从历史灾害库中搜索相似案例,对比分析新旧灾害案例,对历史灾害需求进行调整,从而获取本次灾害的需求。⑶基于空间遥感和GIS技术的应急物资需求预测方法。一些学者对利用高清卫星遥感、航拍遥感影像结合GIS来快速评估地震等巨灾后损失情况,评估结果可用于应急物资需求预测的重要依据。⑷基于时间序列的应急物资需求预测方法。该方法适用于动态需求预测,目前基于该类理论的自回归移动平滑法、指数平滑法和独立同分步法已经广泛用于应急需求预测。⑸基于神经网络的应急物资需求预测方法。由于应急物资需求影响的多元随机性,传统的数学建模难以精确预测,而神经网络因为其广泛的适应能力和学习能力,在非线性系统的预测方面得到广泛的应用;⑹不确定信息条件下实时应急需求预测方法,其目标是预测给定时期内每个受灾地区的动态应急需求。
但应急救援资源的需求预测,受到诸多社会、环境等因素影响,具有很强的时效性和阶段性,同时灾情和物资需求信息具有模糊性和不确定性,灾情发生短时间内获取的信息极为有限,需要结合对灾区历史统计数据挖掘基础上进行的资源需求分布预测,以保证预测的准确性。以往的模型缺乏将历史物资需求数据挖掘规则结合到不确定信息融合应急物资需求预测中去,因而难以获得准确的应急救援物资需求。
发明内容
本发明的目的是提供一种不确定信息融合的应急物资分布需求预测方法,有效解决不确定性信息融合下的救援物资需求预测问题,可有效提高救援物资需求准确性,为航空救援调度和实施提供决策依据。
本发明公开的不确定信息融合的应急物资分布需求预测方法,包括如下步骤,
步骤1,分析历史灾情数据库中各范例的基本属性,提取范例的共同特征,将这些特征运用粗糙集方法进行简约;
步骤2,根据简约后的各条件属性对物资预测结果的影响程度计算出属性权重值;步骤3,对条件属性值和决策属性值相同的范例进行进一步简约;
步骤4,将进一步简约后的灾情数据运用K-means方法将案例分为K类,找出这些类的中心点,并根据当前案例与各类中心点之间的距离判断当前案例属于哪一类;
步骤5,根据范例相似度的计算方法,检索出相似度最高的历史目标范例;
步骤6,根据目标范例的物资消耗以及范例属性权重值来线性推测当前范例的物资需求量。
进一步的,步骤1包括:
步骤11,设S=(C,B)为历史灾情范例数据库,Cn为第n个范例,B为范例属性所组成的集合,即B=F∪D,
其中,F={f1,f2,...,fm}为灾情范例的条件属性集,即和地震有关的情景特征因素信息集,fm为第m个灾情属性的信息;D={D1,D2,...,Di}为决策属性集,即主要应急物资需求集,Di为灾区第i类物资的需求量,D0表示耐用品需求量,D1表示消耗品需求量;
步骤12,给定各条件属性阈值,范例条件属性值满足阈值要求为1,否则为0,按此规则生成0-1信息表,如果C/ind(F)=C/ind(F-{fm}),则属性fm是可以约简的,否则是不可约简的。
进一步的,步骤2包括:
步骤21,需计算出简约后范例的各条件属性对物资预测结果的影响程度,用属性权重值ωj表示。灾情范例条件属性集F={f1,f2,...,fm}的影响权重集为{ω1,ω2,...,ωm},且满足:
步骤22,n(f)表示灾情范例在条件属性为f时的取值,当n(f)在范例库C={C1,C2,...,Cn}中的取值差异比较大的时候,表明该条件属性对分类的判别影响也比较大,应赋予较高权重值;反之,当n(f)在分类中的取值差异比较小的时候,表明该条件属性对分类的判别影响也不大,应取较小的权重值;
步骤23,将灾情范例数据库中的每一个范例当做一类。范例Ci在条件属性fj下的取值n(fj)为该案例在特征因素fj下的隶属度函数并有:
均方差为:
则可求得灾情范例各条件属性的权重ωj为:
进一步的,步骤3包括:F表示灾情范例条件属性的集合,D表示灾情范例决策属性的集合,则有:rij:des(fi)→des(dj),Fi∩Dj≠φ表示由条件属性到决策属性的决策规则集;cover(rij)=|Fi∩Dj|/|B|,Fi∩Dj≠φ为灾情救援优先级决策规则rij的覆盖度。
进一步的,步骤3还可以根据进一步简约后的范例构造救援优先级决策规则,统计各规则发生的频率并计算其覆盖度,当条件覆盖度大于一定阈值时,说明该规则成立,则可得出当前案例的救援优先级。
进一步的,救援优先级决策规则覆盖度等于各规则出现的频率之和,当覆盖度大于90%时,该规则成立。
进一步的,步骤4中是所采用的反馈补偿的K-means方法,即基于反馈补偿的K值确定算法。
进一步的,基于反馈补偿K-means的历史灾情范例聚类算法,其K值确定算法步骤为:
步骤41,设K=2,计算K个聚类中心之间的欧式距离,再设K=K+1,用同样的方法计算各聚类中心之间的平均间距;
步骤42,对比分类数目为K+1时各聚类中心的平均间距相较于分类数目为K时的平均间隔是否增大,增大则再令K=K+1,重复以上步骤,否则确定K值。
进一步的,步骤5中是所采用相关系数法来进行相似度计算,相关系数法来进行相似度计算的步骤包括:
对灾情范例间的相似度采用如下相关系数进行计算,
相关系数的取值范围是[-1,1],计算出的相关系数绝对值越大,表明灾情范例Cx与灾情范例Cy之间的相关度越高。
进一步的,步骤6中,假设相似度最大的历史目标灾情范例的属性分别为P1、P2、P3和P4,耐用品和消耗品的供应量分别为N1、N2,条件属性f1,f2,f3,f4的权重值分别为ω(f1),ω(f2),ω(f3),ω(f4)当前灾情范例的对应属性分别为P1′、P2′、P3′和P4′,则当前灾情范例耐用品和消耗品的需求量N1′和N2′分别为:
N1′=(P1/P1′×ω(f1)+P2/P2′×ω(f2)+P3/P3′×ω(f3)+P4/P4′×ω(f4))×N1
N2′=(P1/P1′×ω(f1)+P2/P2′×ω(f2)+P3/P3′×ω(f3)+P4/P4′×ω(f4))×N2。
本发明具有如下有益效果如下:
1)利用粗糙集属性简约方法,属性的维数大大降低,有效减少灾害范例的冗余属性,提高运行效率,同时也在不失真的情况下节省了大量灾情数据的存储空间。
2)根据简约后的各条件属性对物资预测结果的影响程度计算出属性权重值的方法不仅克服了仅由专家主观性判断带来的误差,也提高了应急救援物资需求预测的精度。
3)本发明基于反馈补偿的K-means聚类算法和范例推理理论引入灾后应急救援物资需求预测这一领域,解决了不确定性救援信息融合条件下的预测准确性,提高了预测的准确性,为救援决策实施提供依据。
4)采用相关系数法来进行相似度计算较夹角余弦、杰卡德相似系数度量法等,相对误差更低,判断更加精准。
5)该发明可直接应用于航空应急救援指挥系统中,有效提高航空救援的效率。
6)根据本发明得到受灾点的救援优先级以及受灾点的物资需求量的预测结果,在救援调度过程中,根据救援优先级是确定救援的先后顺序,然后根据该地物资物资需求,并结合物资分配标准,进行调度。
附图说明
图1为本发明的应急物资需求预测流程图。
图2为基于反馈补偿的K值确定流程图。
图3为K=2时的范例聚类结果。
图4为K=3时的范例聚类结果。
图5为K=4时的范例聚类结果。
图6为K=5时的范例聚类结果。
具体实施方式
本发明运用基于反馈补偿K-means范例推理法预测灾区应急物资需求量的基本流程如图1所示,具体内容包括以下六个部分:
一、基于粗糙集的历史灾情范例属性简约
设S=(C,B)为历史灾情范例数据库,Cn为第n个范例,B为范例属性所组成的集合,即B=F∪D。其中,F={f1,f2,...,fm}为灾情范例的条件属性集,即和地震有关的情景特征因素(如总人口、总面积、震级、震源深度、最高烈度、受灾人数、伤亡人数、砖混比例等)信息集,fm为第m个灾情属性的信息;D={D1,D2,...,Di}为决策属性集,即主要应急物资需求集,Di为灾区第i类物资的需求量。ωi为第i个范例的救援优先级。
给定各条件属性的阈值,若范例条件属性值满足阈值要求为1,否则为0,按此规则生成0-1信息表。
如果C/ind(F)=C/ind(F-{fm}),即条件属性集F和条件属性集F-{fm}对任意范例C的决策属性无任何影响,二者造成的决策结果是一致的,则属性fm是可以约简的,否则是不可约简的。灾情范例简约是粗糙集理论的核心内容之一,它能在持续原有分类能力的情况下,去除冗余的数据和属性。
选取2008~2012年间的22个历史灾情数据作为训练样本(见表1),样本的条件属性为:总人口(人)f1,总面积(平方公里)f2,震级(级)f3,震源深度(级)f4,最高烈度(级)f5,受灾人数(人)f6,伤亡人数(人)f7,砖混比例f8;决策属性为:耐用品需求量D0,消耗品需求量D1。将样本的各条件属性按表2的规则进行离散化处理,超过阈值的属性值设为1,否则设为0,离散化处理后的结果如表3所示。
表1 2008~2012年地震灾情历史范例数据
表2 离散化处理规则
例如,获得的一条范例信息为:{总人口30000人,总面积332.9平方公里,震级6.2,震源深度12,最高烈度8,受灾人数2857,伤亡人数1489,砖混比例0.968},参照表2将此条范例进行离散化处理,得到的结果为{0,0,1,1,1,0,0,1},如表3中第一行范例数据所示。
表3 0-1离散化后的范例信息
从表3中可以看出,属性f1、f2对各范例的属性值相同,故只保留其中一个即可,在此保留属性f1;相同的,属性f3、f4、f5对各范例影响的属性值相同,保留属性f3;属性f6、f7对各范例的属性值相同,保留属性f7。属性简约后的范例信息如表4所示:
表4 简约后的范例信息
属性简约后得到新的条件属性集为F={f1,f3,f7,f8}。利用粗糙集属性简约方法,属性的维数从8维降到了4维,能够有效减少灾害范例的冗余属性,提高运行效率,同时也在不失真的情况下节省了大量灾情数据的存储空间。
二、灾情范例属性权重值的计算
在不同的决策环境下,灾情范例的各条件属性对决策结果会有不同程度的影响,在此,需计算出简约后范例的各条件属性对物资预测结果的影响程度,用属性权重值ωj表示。
令灾情范例条件属性集F={f1,f2,...,fm}的影响权重集为{ω1,ω2,...,ωm},且满足:
令n(f)表示灾情范例在条件属性为f时的取值,即,隶属度函数值。当n(f)在范例库C={C1,C2,...,Cn}中的取值差异比较大的时候,表明该条件属性对分类的判别影响也比较大,应赋予较高权重值;反之,当n(f)在分类中的取值差异比较小的时候,表明该条件属性对分类的判别影响也不大,应取较小的权重值。因此,可将灾情范例数据库中的每一个范例当做一类。范例Ci在条件属性fj下的取值n(fj)为该案例在特征因素fj下的隶属度函数nCi(fj),并有:
均方差为:
则可求得灾情范例各条件属性的权重ωj为:
将第一部分算例简约后的属性按照上述给出的权重值计算步骤求出条件属性f1、f3、f7、f8的权值分别为ω1=0.13、ω3=0.36、ω7=0.31、ω8=0.20,权重值越大,说明该条件属性对分类判别的作用越大,反之,权重值越小,说明该条件属性对分类判别的作用越小。在本发明的算例中,对分类判别影响最大的因素是震级,总人口数对分类判别的影响最小。这种权值计算方法不仅克服了仅由专家主观性判断带来的误差,也提高了应急救援物资需求预测的精度。
三、救援优先级决策规则覆盖度的计算
按照第一部分中灾情范例数据库的定义,则有:rij:des(fi)→des(dj),Fi∩Dj≠φ表示由条件属性到决策属性的决策规则集。定义:cover(rij)=|Fi∩Dj|/|B|,Fi∩Dj≠φ为灾情救援优先级决策规则rij的覆盖度,即在灾情范例数据库中,满足决策规则集rij的对象数占整个数据库对象数的比例。
对表4训练样本中条件属性值和决策属性值相同的范例进行简约所得决策规则如表5所示,运用本项目组开发的低空应急救援系统中的500条历史灾情数据构建范例数据库,统计决策规则发生的频率如下:
表5 救援优先级决策规则
救援优先级决策规则覆盖度等于各规则出现的频率之和,当覆盖度大于90%时,则该规则成立。通过表5中各规则出现的频率,计算出满足该决策的样本数占整体样本数的92.3%,即规则覆盖率为92.3%,说明该规则成立,可用于实际救援。
由表5可知,救援优先级的判断规则可由总人口、震级、伤亡人数、砖混比例这四个属性决定,例如,当范例总人口数小于10万、震级是大于等于6、伤亡人数小于3000、房屋砖混比例大于等于0.6时,根据表5类比范例C1可推断出该范例的救援优先级是2;当总人口数大于等于10万、震级是大于等于6、伤亡人数大于等于3000、房屋砖混比例大于等于0.6时,类比范例C14可推断出当前范例的救援优先级是5,通过此决策规则可对各灾区的救援的优先级进行判定。四、基于反馈补偿K-means的历史灾情范例聚类算法
聚类分析方法作为一种被广泛使用的数据挖掘算法,能以较高的效率处理大数据,获得这些数据的分布特征及分布规则,并逐步应用于各领域。基于不同思想的聚类算法被提出,主要有基于划分的算法、基于层次的算法、基于密度的算法、基于网格的算法和基于模型的算法等,例如,K-means聚类算法、层次聚类算法、SOM聚类算法以及FCM聚类算法等。这些算法都能取得不错的聚类效果,其中应用最多且算法思想较为简单的是基于划分的K-means算法。
针对历史灾情数据多而杂的特点,本发明运用计算时间短、速度快、聚类效果较好的反馈补偿K-means算法对范例进行聚类分析,获得数据的分类,为第四部分找出与当前范灾情例相似度最高的历史灾情范例提供数据支持,并能提高对历史范例数据库检索的效率。
K-means算法采用欧式距离作为两个灾情范例相似度的评价指标,其基本思想是:随机选取历史灾情范例数据库中的K个范例作为初始聚类的中心,根据范例数据库中各灾情范例到K个中心的距离将其归类到距离最小的类中,然后计算各个类中灾情范例距离的平均值,更新每个类的中心,直到聚类中心不再发生变化。它的目标是将历史灾情范例数据库中的范例分成若干个类,使得同一类内的灾情范例之间的相似度尽可能大,不同类的灾情范例之间的相似度尽可能小。但,随着K-means算法研究的深入,一些不足也逐渐暴露出来,主要包括:K值难以预先确定以及初始的聚类中心选择困难等。K值的确定以及初始聚类中心的选择始终是K-means算法的重点及难点。
本发明提出的基于反馈补偿K-means的历史灾情范例聚类算法,其K值确定算法步骤为:首先,设K=2,计算K个聚类中心之间的欧式距离,再设K=K+1,用同样的方法计算各聚类中心之间的平均间距,对比分类数目为K+1时各聚类中心的平均间距相较于分类数目为K时的平均间隔是否增大,增大则再令K=K+1,重复以上步骤,否则确定K值。此方法将灾情范例类间平均距离、类间重叠情况及类内数据的分散程度作为类内数目选择的重要参考,即K值的确定。假设为K个灾情范例聚类中心之间的平均距离,反馈补偿K值确定过程如图2所示。
在反馈补偿K-means灾情范例聚类算法中初始聚类中心点的选择规则为:第一个初始聚类中心点从历史灾情范例数据库中随机生成,第二个初始聚类中心点离第一个初始点最远,第三个初始聚类中心点离第一、第二个初始点的平均距离最远,反复如此,选出K个初始聚类中心。
对范例数据库中简约后的灾情数据运用K-means算法进行分类时,应先确定分类数,即K值。首先,设K=2,计算出2个聚类中心之间的间距(即欧式距离)为217138.9999;再设K=3,用同样的方法计算出3个聚类中心之间的平均间距为308223.8637,对比可知,K=3时聚类中心的平均间距相较于K=2的平均间隔明显增大,说明K=3时的聚类效果要比K=2时的聚类效果好;再设K=4,计算出4个聚类中心之间的平均间距为320208.71559442,对比可知,K=4时聚类中心的平均间距相较于K=3的平均间隔明显增大,说明K=4时的聚类效果要比K=3时的聚类效果更好;再设K=5,计算出5个聚类中心之间的平均间距为274362.7997,对比可知,K=5时聚类中心的平均间距相较于K=4的平均间隔明显减小,说明K=5时的聚类效果要比K=4时的聚类效果更差,分类不明显,故确定K值为4。如图3、4、5、6分别是K=2、K=3、K=4、K=5时的范例聚类结果。将类间平均距离、类间重叠情况及类内数据的分散程度作为类内数目选择的重要参考,通过对比以上实验聚类结果,可明显看出K=2时范例聚类效果太分散,K=3及K=5时的范例聚类结果中各分类重叠区域较大,当K=4时范例聚类效果较好,分类清晰可见。故选取K=4作为分类数。
五、灾情范例相似度的计算
范例相似度计算的常用方法一般包括夹角余弦计算法、杰卡德相似系数法以及相关系数法等,其中夹角余弦更多的是从方向上区分差异,而对绝对的数值不敏感;杰卡德相似系数处理的都是非对称二元变量,并没有考虑向量中潜在数值的大小,而是简单的处理为0和1;相关系数法是度量随机向量之间的相关程度的一种方法,本发明采用相关系数对灾情范例间的相似度进行计算,定义如下:
相关系数的取值范围是[-1,1],计算出的相关系数绝对值越大,表明灾情范例Cx与灾情范例Cy之间的相关度越高。
运用反馈补偿K-means方法对灾情数据分类后,计算当前灾情范例与各分类中心点的距离,判断当前案例属于哪一类,然后再将当前灾情范例与该类中每个算例进行相似度计算,找出相似度最大历史灾情范例。
从历史灾情范例中选出X=(129320,6.3,4967,0.833)作为验算该算法的当前范例,当K=4时,各分类中心点数据如表6所示:
表1 当前范例与各中心点距离表
从表6可知,当前灾情范例X与上述四个中心点的欧式距离分别为:155114.910787806,104617.455838721,31896.3656220504,376667.762803906,与第二类中心点O2的距离最近,故将当前范例视为第二类。
运用MATLAB将属于第二类的数据筛选出来,并将当前范例X=(129320,6.3,4967,0.833)与第二类中的每个范例进行相似度计算。计算结果中最大的相关系数为0.999983050183659,其对应的历史案例是C141=(136455,8.4,4418,0.783),故可参照该范例的历史需求数据对当前范例的物资需求量进行预测。
六、当前灾情范例物资需求预测
假设相似度最大的历史目标灾情范例的属性分别为P1、P2、P3和P4,耐用品和消耗品的供应量分别为N1、N2,条件属性f1,f2,f3,f4的权重值分别为ω(f1),ω(f2),ω(f3),ω(f4)当前灾情范例的对应属性分别为P1′、P2′、P3′和P4′,则当前灾情范例耐用品和消耗品的需求量N1′和N2′分别为:
N1′=(P1/P1′×ω(f1)+P2/P2′×ω(f2)+P3/P3′×ω(f3)+P4/P4′×ω(f4))×N1
N2′=(P1/P1′×ω(f1)+P2/P2′×ω(f2)+P3/P3′×ω(f3)+P4/P4′×ω(f4))×N2
通过查找范例数据库可知,历史范例C141=(136455,8.4,4418,0.783)对耐用品和消耗品的物资需求量分别为1100kg和700kg,利用线性相关算法对当前灾情范例的耐用品和消耗品的物资需求量进行预测,结果分别为:1049kg和668kg。即当前案例X=(129320,6.3,4967,0.833)对耐用品的需求量为1049kg,对消耗品的需求量为668kg。
综合实施例的步骤,可得到两个结果数据,一个是受灾点的救援优先级,另一个是受灾点的物资需求量,在救援调度过程中,救援优先级是物资分配的参考标准,通常,对救援优先级高的受灾点予以优先救援,并且会先满足该点的物资需求,然后再逐个考虑下一救援优先级受灾点的物资需求的分配。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下还可以作出若干改进,这些改进也应视为本发明的保护范围。
Claims (10)
1.一种不确定信息融合的应急物资分布需求预测方法,其特征在于:包括如下步骤,
步骤1,分析历史灾情数据库中各范例的基本属性,提取范例的共同特征,将这些特征运用粗糙集方法进行简约;
步骤2,根据简约后的各条件属性对物资预测结果的影响程度计算出属性权重值;
步骤3,对条件属性值和决策属性值相同的范例进行进一步简约;
步骤4,将进一步简约后的灾情数据运用K-means方法将案例分为K类,找出这些类的中心点,并根据当前案例与各类中心点之间的距离判断当前案例属于哪一类;
步骤5,根据范例相似度的计算方法,检索出相似度最高的历史目标范例;
步骤6,根据目标范例的物资消耗以及范例属性权重值来线性推测当前范例的物资需求量。
2.如权利要求2所述的应急物资分布需求预测方法,其特征在于:所述步骤1包括:
步骤11,设S=(C,B)为历史灾情范例数据库,Cn为第n个范例,B为范例属性所组成的集合,即B=F∪D,
其中,F={f1,f2,...,fm}为灾情范例的条件属性集,即和地震有关的情景特征因素信息集,fm为第m个灾情属性的信息;D={D1,D2,...,Di}为决策属性集,即主要应急物资需求集,Di为灾区第i类物资的需求量,D0表示耐用品需求量,D1表示消耗品需求量;
步骤12,给定各条件属性阈值,范例条件属性值满足阈值要求为1,否则为0,按此规则生成0-1信息表,如果C/ind(F)=C/ind(F-{fm}),则属性fm是可以约简的,否则是不可约简的。
3.如权利要求1所述的应急物资分布需求预测方法,其特征在于:所述步骤2包括:
步骤21,需计算出简约后范例的各条件属性对物资预测结果的影响程度,用属性权重值ωj表示。灾情范例条件属性集F={f1,f2,...,fm}的影响权重集为{ω1,ω2,...,ωm},且满足:
步骤22,n(f)表示灾情范例在条件属性为f时的取值,当n(f)在范例库C={C1,C2,...,Cn}中的取值差异比较大的时候,表明该条件属性对分类的判别影响也比较大,应赋予较高权重值;反之,当n(f)在分类中的取值差异比较小的时候,表明该条件属性对分类的判别影响也不大,应取较小的权重值;
步骤23,将灾情范例数据库中的每一个范例当做一类。范例Ci在条件属性fj下的取值n(fj)为该案例在特征因素fj下的隶属度函数并有:
均方差为:
则可求得灾情范例各条件属性的权重ωj为:
4.如权利要求1所述的应急物资分布需求预测方法,其特征在于:所述步骤3包括:
F表示灾情范例条件属性的集合,D表示灾情范例决策属性的集合,则有:
rij:des(fi)→des(dj),Fi∩Dj≠φ表示由条件属性到决策属性的决策规则集;
cover(rij)=|Fi∩Dj|/|B|,Fi∩Dj≠φ为灾情救援优先级决策规则rij的覆盖度。
5.如权利要求1所述的应急物资分布需求预测方法,其特征在于:所述步骤3包括:根据进一步简约后的范例构造救援优先级决策规则,统计各规则发生的频率并计算其覆盖度,当条件覆盖度大于一定阈值时,说明该规则成立,则可得出当前案例的救援优先级。
6.如权利要求5所述的应急物资分布需求预测方法,其特征在于:救援优先级决策规则覆盖度等于各规则出现的频率之和,当覆盖度大于90%时,该规则成立。
7.如权利要求1所述的应急物资分布需求预测方法,其特征在于:所述步骤4中是所采用的反馈补偿的K-means方法,即基于反馈补偿的K值确定算法。
8.如权利要求7所述的应急物资分布需求预测方法,其特征在于:基于反馈补偿K-means的历史灾情范例聚类算法,其K值确定算法步骤为:
步骤41,设K=2,计算K个聚类中心之间的欧式距离,再设K=K+1,用同样的方法计算各聚类中心之间的平均间距;
步骤42,对比分类数目为K+1时各聚类中心的平均间距相较于分类数目为K时的平均间隔是否增大,增大则再令K=K+1,重复以上步骤,否则确定K值。
9.如权利要求1所述的应急物资分布需求预测方法,其特征在于:所述步骤5中是所采用相关系数法来进行相似度计算,相关系数法来进行相似度计算的步骤包括:
对灾情范例间的相似度采用如下相关系数进行计算,
相关系数的取值范围是[-1,1],计算出的相关系数绝对值越大,表明灾情范例Cx与灾情范例Cy之间的相关度越高。
10.如权利要求1所述的应急物资分布需求预测方法,其特征在于:所述步骤6包括:
假设相似度最大的历史目标灾情范例的属性分别为P1、P2、P3和P4,耐用品和消耗品的供应量分别为N1、N2,条件属性f1,f2,f3,f4的权重值分别为ω(f1),ω(f2),ω(f3),ω(f4)当前灾情范例的对应属性分别为P1′、P2′、P3′和P4′,则当前灾情范例耐用品和消耗品的需求量N1′和N2′分别为:
N1′=(P1/P1′×ω(f1)+P2/P2′×ω(f2)+P3/P3′×ω(f3)+P4/P4′×ω(f4))×N1,
N2′=(P1/P1′×ω(f1)+P2/P2′×ω(f2)+P3/P3′×ω(f3)+P4/P4′×ω(f4))×N2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710193250.5A CN106997509B (zh) | 2017-03-28 | 2017-03-28 | 一种不确定信息融合的应急物资分布需求预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710193250.5A CN106997509B (zh) | 2017-03-28 | 2017-03-28 | 一种不确定信息融合的应急物资分布需求预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106997509A true CN106997509A (zh) | 2017-08-01 |
CN106997509B CN106997509B (zh) | 2021-02-05 |
Family
ID=59431724
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710193250.5A Active CN106997509B (zh) | 2017-03-28 | 2017-03-28 | 一种不确定信息融合的应急物资分布需求预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106997509B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108022051A (zh) * | 2017-12-13 | 2018-05-11 | 吉林大学 | 一种基于信息融合的城市轨道交通调度决策方法 |
CN109443419A (zh) * | 2018-08-31 | 2019-03-08 | 广州市世科高新技术有限公司 | 一种基于机器学习的整流器在线监测方法 |
JP2019067185A (ja) * | 2017-10-02 | 2019-04-25 | ヤフー株式会社 | 情報処理装置、販売制御方法および販売制御プログラム |
CN109741715A (zh) * | 2019-02-25 | 2019-05-10 | 深圳市华星光电技术有限公司 | 显示面板的补偿方法、补偿装置及存储介质 |
CN110288118A (zh) * | 2019-05-14 | 2019-09-27 | 湖北工业大学 | 一种基于无人机的空中救援响应方法 |
CN110991763A (zh) * | 2019-12-13 | 2020-04-10 | 中国民用航空飞行学院 | 一种基于指标模糊分割和topsis的通航应急救援资源需求预测方法 |
CN111080080A (zh) * | 2019-11-25 | 2020-04-28 | 桂林理工大学南宁分校 | 一种村镇地质灾害风险预估方法及系统 |
CN113112074A (zh) * | 2021-04-14 | 2021-07-13 | 许人杰 | 一种基于物联网的智慧应急资源调动系统及其方法 |
CN113344356A (zh) * | 2021-05-31 | 2021-09-03 | 烽火通信科技股份有限公司 | 一种多目标资源分配决策方法和装置 |
CN114186711A (zh) * | 2021-10-27 | 2022-03-15 | 中山大学 | 一种基于多任务时序学习的工业原料消耗量预测方法 |
CN114579712A (zh) * | 2022-05-05 | 2022-06-03 | 中科雨辰科技有限公司 | 基于动态模型的文本属性提取匹配方法 |
CN116050676A (zh) * | 2023-03-31 | 2023-05-02 | 四川互慧软件有限公司 | 一种基于突发群体事件的应急预案推荐系统 |
CN117035374A (zh) * | 2023-10-09 | 2023-11-10 | 北京数易科技有限公司 | 一种应对突发事件的兵力协同调度方法、系统及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102819772A (zh) * | 2012-08-29 | 2012-12-12 | 广东电网公司 | 电力配网建设物资需求预测方法及装置 |
CN102831489A (zh) * | 2012-08-29 | 2012-12-19 | 广东电网公司 | 电力配网建设物资需求预测方法及装置 |
CN105139021A (zh) * | 2015-07-08 | 2015-12-09 | Tcl集团股份有限公司 | 一种基于粗糙集理论实现电视用户快速分类的方法及系统 |
CN105760888A (zh) * | 2016-02-23 | 2016-07-13 | 重庆邮电大学 | 一种基于属性聚类的邻域粗糙集集成学习方法 |
-
2017
- 2017-03-28 CN CN201710193250.5A patent/CN106997509B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102819772A (zh) * | 2012-08-29 | 2012-12-12 | 广东电网公司 | 电力配网建设物资需求预测方法及装置 |
CN102831489A (zh) * | 2012-08-29 | 2012-12-19 | 广东电网公司 | 电力配网建设物资需求预测方法及装置 |
CN105139021A (zh) * | 2015-07-08 | 2015-12-09 | Tcl集团股份有限公司 | 一种基于粗糙集理论实现电视用户快速分类的方法及系统 |
CN105760888A (zh) * | 2016-02-23 | 2016-07-13 | 重庆邮电大学 | 一种基于属性聚类的邻域粗糙集集成学习方法 |
Non-Patent Citations (4)
Title |
---|
ZHI-HUA HU: "Relief demand forecasting in emergency logistics based on tolerance model", 《2010 2ND IEEE INTERNATIONAL CONFERENCE ON INFORMATION MANAGEMENT AND ENGINEERING》 * |
喻慧: "基于反馈补偿K-means的救援物资需求预测", 《航空计算技术》 * |
段在鹏: "基于模糊粗糙集的化工园区多米诺事故救援时机分析", 《INDUSTRIAL SAFETY AND ENVIRONMENTAL PROTECTION》 * |
赵小柠等: "基于范例推理的灾害性地震应急物资需求预测研究", 《中国安全科学学报》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7023078B2 (ja) | 2017-10-02 | 2022-02-21 | ヤフー株式会社 | 情報処理装置、販売制御方法および販売制御プログラム |
JP2019067185A (ja) * | 2017-10-02 | 2019-04-25 | ヤフー株式会社 | 情報処理装置、販売制御方法および販売制御プログラム |
CN108022051A (zh) * | 2017-12-13 | 2018-05-11 | 吉林大学 | 一种基于信息融合的城市轨道交通调度决策方法 |
CN109443419A (zh) * | 2018-08-31 | 2019-03-08 | 广州市世科高新技术有限公司 | 一种基于机器学习的整流器在线监测方法 |
CN109443419B (zh) * | 2018-08-31 | 2021-08-31 | 广州致新电力科技有限公司 | 一种基于机器学习的整流器在线监测方法 |
CN109741715A (zh) * | 2019-02-25 | 2019-05-10 | 深圳市华星光电技术有限公司 | 显示面板的补偿方法、补偿装置及存储介质 |
CN110288118A (zh) * | 2019-05-14 | 2019-09-27 | 湖北工业大学 | 一种基于无人机的空中救援响应方法 |
CN110288118B (zh) * | 2019-05-14 | 2022-04-01 | 湖北工业大学 | 一种基于无人机的空中救援响应方法 |
CN111080080A (zh) * | 2019-11-25 | 2020-04-28 | 桂林理工大学南宁分校 | 一种村镇地质灾害风险预估方法及系统 |
CN110991763A (zh) * | 2019-12-13 | 2020-04-10 | 中国民用航空飞行学院 | 一种基于指标模糊分割和topsis的通航应急救援资源需求预测方法 |
CN110991763B (zh) * | 2019-12-13 | 2022-12-02 | 中国民用航空飞行学院 | 一种基于指标模糊分割和topsis的通航应急救援资源需求预测方法 |
CN113112074A (zh) * | 2021-04-14 | 2021-07-13 | 许人杰 | 一种基于物联网的智慧应急资源调动系统及其方法 |
CN113344356A (zh) * | 2021-05-31 | 2021-09-03 | 烽火通信科技股份有限公司 | 一种多目标资源分配决策方法和装置 |
CN114186711A (zh) * | 2021-10-27 | 2022-03-15 | 中山大学 | 一种基于多任务时序学习的工业原料消耗量预测方法 |
CN114579712A (zh) * | 2022-05-05 | 2022-06-03 | 中科雨辰科技有限公司 | 基于动态模型的文本属性提取匹配方法 |
CN114579712B (zh) * | 2022-05-05 | 2022-07-15 | 中科雨辰科技有限公司 | 基于动态模型的文本属性提取匹配方法 |
CN116050676A (zh) * | 2023-03-31 | 2023-05-02 | 四川互慧软件有限公司 | 一种基于突发群体事件的应急预案推荐系统 |
CN116050676B (zh) * | 2023-03-31 | 2023-06-16 | 四川互慧软件有限公司 | 一种基于突发群体事件的应急预案推荐系统 |
CN117035374A (zh) * | 2023-10-09 | 2023-11-10 | 北京数易科技有限公司 | 一种应对突发事件的兵力协同调度方法、系统及介质 |
CN117035374B (zh) * | 2023-10-09 | 2024-01-02 | 北京数易科技有限公司 | 一种应对突发事件的兵力协同调度方法、系统及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106997509B (zh) | 2021-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106997509A (zh) | 一种不确定信息融合的应急物资分布需求预测方法 | |
Sheu | Dynamic relief-demand management for emergency logistics operations under large-scale disasters | |
Xiong et al. | A decision support method for design and operationalization of search and rescue in maritime emergency | |
Sarma et al. | Uncertain demand estimation with optimization of time and cost using Facebook disaster map in emergency relief operation | |
CN116108758B (zh) | 滑坡易发性评价方法 | |
CN104715318A (zh) | 面向通信网络的多维度运行风险评价方法 | |
Csikós et al. | Traffic speed prediction method for urban networks—An ANN approach | |
Hu | Multi-objective genetic algorithm for berth allocation problem considering daytime preference | |
Shao et al. | Relief demand forecasting based on intuitionistic fuzzy case-based reasoning | |
Lahdelma et al. | Stochastic multicriteria acceptability analysis (SMAA) | |
CN109934398A (zh) | 一种钻爆法施工隧道瓦斯危险等级预测方法和装置 | |
Liu et al. | Humanitarian logistics planning for natural disaster response with Bayesian information updates. | |
Khaldi et al. | Performance prediction of pharmaceutical suppliers: comparative study between DEA-ANFIS-PSO and DEA-ANFIS-GA | |
Cai et al. | Bi-level programming enabled design of an intelligent maritime search and rescue system | |
Massel et al. | The use of machine learning in situational management in relation to the tasks of the power industry | |
Sassirekha et al. | Predicting the academic progression in student’s standpoint using machine learning | |
Zhou et al. | The landing safety prediction model by integrating pattern recognition and Markov chain with flight data | |
CN105243474A (zh) | 基于时空信息的尾矿库安全风险评估方法 | |
Kai et al. | Notice of Retraction: A Novel Forecasting Model of Fuzzy Time Series Based on K-means Clustering | |
Setianto et al. | Employee Turnover Analysis Using Comparison of Decision Tree and Naive Bayes Prediction Algorithms on K-Means Clustering Algorithms at PT. AT | |
Rokhsari et al. | Urban network risk assessment using Fuzzy-AHP and TOPSIS in GIS environment | |
CN106339357A (zh) | 寻优参数化的网络信息安全检测系统和方法 | |
Hasperué et al. | Rule extraction on numeric datasets using hyper-rectangles | |
Pang et al. | Wt model & applications in loan platform customer default prediction based on decision tree algorithms | |
Ouallane et al. | Towards Intelligent Road Traffic Management Based on Neutrosophic Logic: A Brief Review |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |