毁灭性灾难中救灾物资分配方案的数学模型研究
№.1・158・
陕西科技大学学报
JOURNAl。OFSHAANXIUNIVERSITYOF.SCIENCE&TECHNOI.OGY
Feb.2009
Vol-27
’文章编号:1000-5811{2009)01-0158—06
毁灭性灾难中救灾物资分配方案的数学模型研究
谭秋月1,孙平安2,韩希超3
(1.武夷学院经济与数学系,福建武夷山354300;2.武夷学院信息网络与教育技术中心,福建武夷山354300;3.辽宁科技大学计算机科学与工程学院,辽宁鞍山114051)
摘要:针对自然灾害发生后灾区救灾物资数量有限的情况,利用群组聚类的思想制定出了救灾物资详细的分配方案、相应的分配原则和物资分配策略,并给出了相应的数学模型或数学描述,通过一个简单直观的数值算例进一步验iz-'f该模型的合理性和可行性.关键词:救灾物资;群组聚类;分配模型;数值算例中图分类号:0224
文献标识码:A
0引言
当某一地区受灾后,根据灾区灾害发生程度的不同,将灾情分为毁灭性(Destroyed)和非毁灭性(None-Destroyed)两种.假设某一灾区共有N名受灾群众,现有一批救灾物资要发放给这些受灾者;物资共有M种,每种物资的数量有限;各受灾者的灾情不同,对每种物资的急需程度和需求量不同.我们要研究的问题是:(1)制定分配原则并给出合理的分配方案,(2)试给出一个符合假设情况的数值算例.1模型假设及符号说明1.1模型假设[1.2]
(1)假设灾区有N名受灾者,且每一名受灾者均需要救灾物资;(2)假设灾区的每名受灾者之间相互独立,不考虑其家庭因素等;(3)假设灾区受灾者所经历的受灾灾种相同;(4)假设所运来的救灾物资为M种,且每种救灾物资均不足;(5)假设所运来的M种救灾物资均为对受灾者可分配的救灾物资;(6)假设对救灾物资一律实行无偿发放.1.2符号说明
N:灾区受灾者总人数;M:运送到灾区的救灾物资总类数;6,:运往灾区的第歹种救灾物资(以下简称物资歹)的总数量,其中歹=1,2,…,M;B:运往灾区的M种救灾物资的总数量矩阵,即B=(岛)。×M;q:物资J的剩余可分配数量,其中歹一1,2’..・,M;qd:第i名受灾者(以下简称受灾者i)对物资j的需求量;Q:
gll921
q12
‘‘‘
qlM
q2M
受灾者对各种救灾物资的需求量矩阵,即:Q=(qd)N×一r=
qzz
…
;K:按急需程度对受灾者
qMqN2
…qMf
进行分类得到的群组数目;G:毁灭性灾害模型中对受灾者分类得到的第志个群组;皖:非毁灭性灾害模
型中对物资歹分配时将受灾者分类得到的第k个群组;蹰:对物资歹把受灾者分组后统计得到的第k组的受
・收稿日期;2008一lo~16
作者简介:谭秋月(1980--),女。陕西省杨陵区人.助教,在读硕士生,研究方向:图论及其应用
第1期
。!!!!!!!!!—!!!
f‘Oll
COlZ
谭秋月等:毁灭性灾难中救灾物资分配方案的数学模型研究
一.
====================II==
・159・
II!!!!!
灾者人数;叫。:受灾者i对物资j;的急需程度;W:各受灾者对救灾物资急需程度矩阵,即:W=(哟)NxM=
…∞1^f
。‘‘
l∞21∞22
WZM
卜
i
CONZ
。。‘
;^:毁灭性灾害模型中群组G内各受灾者对物资_f分配方案的约束函数,其中k
10/诅4
LOJNl
=1,2,…,K;厶:非毁灭性灾害模型中受灾者i对物资歹分配方案的约束函数,其中i一1,2,…,N;粕:受灾者i所分配到的第J种救灾物资的数量;X:题意所求救灾物资的分配方案,即:X=(奶)N州=
…Z1M
‘‘’X2.M
…X.MVI
2模型的综述
灾情发生后,根据灾情发生的灾害程度,将灾情分为毁灭性灾种和非毁灭性灾种.分别选取这两种灾情状况下能够描述受灾者灾情及对救灾物资的急需程度和需求量的某几个属性值作为已知条件,利用群组聚类的思想决定两种分配方案:
(1)毁灭性灾种,体现受灾者对救灾物资的紧迫性,将灾区受灾者按照其对救灾物资需求的某一种或者几种属性值进行聚类分组,然后对各群组受灾者实现对救灾物资的平均分配.
(2)非毁灭性灾种,对其现有灾情信息的某一种或者几种属性值经综合分析后分组,引人竞争机制,对每种救灾物资按照对每名受灾者边加权修正优先级边分配的原则进行分配,从而使得物资分配更合理,分配方法更科学.
3灾情状况及救灾模型的建立3.1毁灭性灾情状况概述E3,4]
当灾区所受灾种对受灾者来说为某一毁灭性灾种且运送来的该批救灾物资均为紧急救灾物资时,那么首先应该考虑救灾物资分配的紧迫性.对每一名受灾者均无法考虑其受灾程度、受灾损失以及个人经济等因素时,利用聚类的思想,对受灾者按照某一种或者几种属性值进行综合分析,分成几个不同的受灾群组.比如对受灾者按年龄阶段划分为儿童、青年、壮年、老人等群组,再根据分类的结果进行不同受灾群组的优先级(用急需程度表示)排序,上例群组划分中的优先级应为:儿童>老人>青年>壮年.对不同群组受灾者赋予不同的物资需求量,但同一群组受灾者的需求量应相同;然后按照群组优先级从高到低的顺序完成对该物资的分配,当分配到某一群组,而该群组的物资需求总量小于剩余物资量时,应考虑将部分剩余救灾物资先平均分给该群组的每一名受灾者.当剩余的救灾物资总数小于该群组受灾者人数时,应对剩余的救灾物资按照该群组中的个人灾情程度中的某一属性值排序,然后按照由高到低的排序顺序将物资分配到该群组部分受灾者手中以完成该种救灾物资的分配,接下来进行下一种救灾物资分配,直到将所有的物资分配完毕.
3.2毁灭性救灾模型的建立睁7J
在毁灭性灾难中,各受灾者对M种救灾物资均有需求,将N名受灾者按其某一种或者几种属性值进
行综合分析,分成K个群组,记为G1,G2,…,GK,对群组内的各受灾者,同样按照其某一种或者几种属性值进行综合分析,利用综合属性值进行内部排序,用受灾者的序号表示顺序,序号越低则表示在群组内部对物资的急需程度越高,即:
・160・
!!!!!!!!!!!!!!!I==
陕两科技大学学报
I
I
I=。=
I!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
第27卷
如果受灾者i∈G,行∈G,并且,l≤i<以≤N,则在群组G内部,受灾者i对物资歹的急需程度要比受灾者n对物资.『的急需程度高.
救灾人员根据灾区实际情况合理赋予不同群组受灾者以不同的物资需求量,但同一群组受灾者的需求量应相同.设受灾者i对物资J的需求量为qo,则灾区内的N名受灾者对现有救灾物资的需求量为已知.
统计这K个群组中受灾者的人数,分别为:s。,sz,…,蹰,应满足:
K
∑s产N
^=1
(1)
由上述条件,建立以群组内部最公平为约束的数学模型,则灾情状况(1)中可构造出对每个群组G分配物资J的救灾物资分配方法模型为:
min:^一max{ma。,x{qli—q),qo一劫一1)一min{碑!p({劬一zd},qo一奶一1)
~
∑x口≤岛
i=l
z#≤q#
,27d≥0且为整数
口i∈GK七一1,2,…,Ki一1,2,…,NZ一1,2,…,N
(2)
歹=1,2,…,M
将所有受灾者分组以后,下一阶段就需要进行排序操作,以确定对各群组的紧急程度及进行分配的优先次序.在受灾者的群组分类分析中,已经对原始指标数据进行了无量纲化处理,可以将各指标作为群组排序中模糊评估的评估准则,而将原始数据无量纲化所得的数据作为各评估准则的权值,并以各指标权值的高低作为评估群组优先级的依据.根据综合权值的高低对各受灾者予以排序,即可得到各群组的优先
在模型中,进行救灾物资分配时,采取对救灾物资进行轮流分配的方式对某一种救灾物资进行集中分配,即分配完第1种救灾物资之后再进行第2种救灾物资的分配,依次循环下去,直到把所有可分配的救灾物资分配完毕为止.
,
.对物资J进行分配时,按照群组急需程度从高到底的排序顺序依次分配,当对群组G进行物资分配时,比群组G的优先级高的各个群组均已经按照赋予的需求量完成分配并且该救灾物资仍有剩余,同样记救灾物资J的剩余可分配量为cf.
灾情状况1模型中对群组G发放物资j的分配方法步骤如下:首先检验群组G的物资需求总量是否小于等于剩余物资量fj,即:
如果式(3)条件成立,则按照群组G中各受灾者对物资歹的需求量进行分配,调整分配方案中的相关
、
数值铂和救灾物资歹的剩余可分配量c,的值,然后对具有下一优先级的群组G抖,进行该物资分配.
如果式(3)条件不成立,为了保证使得每个群组内部分配后满足约束条件,即:^一max{max{qlj一
嘞},qo一嘞一1)一rain{min{qo一嘞),劬~哪一1)最小,计算可得:
㈦碧。靠
㈤
第1期潭秋月等:毁灭性灾难中救灾物资分配方案的数学模型研究
・161・
然后对群组G中的各受灾者分别分配,.。个该物资,对于仍然剩余的rz个该物资,一定满足rz<&,再对群组G内部的各受灾者按照序号顺序将剩余的,.:个该物资分配给急需程度最高的前r2个受灾者,这样在分配完物资歹时也确保了所得的厂茸的值为最小值,从而完成了救灾物资歹的分配;然后对第歹+1种救灾物资进行如上分配,直到所有的M种救灾物资全部完成分配为止.
4
毁灭性灾难救灾模型数值算例的求解对于毁灭性灾难,我们建立如下算例:
算例l:某一灾区有lo名受灾群众,现有一批救灾物资要发放给这些受灾者.物资共有5种,每种物资
数量有限,详细数据见表1所示.
值进行分类,在此我们按照年龄段将他们分成3个群组,其中百蘸丽可—了——r—1——了1
1磊焉趸_百广—i广—磊——西、
1~18岁分为一组,记为G,,60岁以上分为一组,记为G2,18~60岁分为一组,记为G3.
按照对救灾物资的急需程度,可知G,组的受灾者对物资的急需程度最高,将其急需程度记为A;G2组对物资的急需程度次之,将其急需程度记为B;最后把G3组急需程度记为C.救灾人员按照已经分好的群组对每种救灾物资赋予每个群组内的受灾者以相同的需求量,不同群组间需求量不一定相同,总的趋势是,物资急需程度越高的群组所赋予的需求量越少.算例中的相关数据详见表2所示.
按照上述数据,可得到:N_--10,M一5,B一(31
19
33
13
对于该灾区的10名受灾者,按照其某一种或者几种属性
表1各种救灾物资总量
表2各受灾者对每种物资的需求量
及所属群组的急需程度表
受灾者物资1物资2物资3物资4物资5急需程度
18),Gl一{Ⅱ2,口g},Gz一{口3,口6,
alo),G3={口1,n4,口5,a7,a8),¥1j=2,s≈----3,s3j=5,
52455
4134434413
6343343334
3236636623
6443343344
Q一(qo)10X5一
45524
下面开始对物资l进行分配:由于f1=bl一31>0,而∑q^2sk一2×2+4×3+5×5=41>31,
^;1
则首先转到对急需程度最高的群组G1内的各受灾者进行分配.
检验q,㈣一2×2—4<fl,则对群组G1内的受灾者按照其需求量进行分配,即:zz・=q21=2,X9l=q。,一2,同时调整c。一31—2×2=27>0后,转到对下一个优先群组G2内的受灾者进行剩余物资分配.
检验q:。S。一4×3—12<f1,同样对群组G2内的受灾者按照其需求量进行分配,即:z。。=q6。=z1。1—4,调整cl=27—4×3—15>0,转到对下一群组G3内的受灾者进行剩余物资分配.
再次检验q。。ss一5×5—25>C・,此时剩余物资总量已经不能满足对群组G3中所有受灾者按其需求
・162・
陕两科技大学学报第27卷
量分配,由式(4)得到:.|r】2[詈]2[訾]。3
量分配,由式(4)得到:<‘
Lss-J
L5J
,于是得到对群组Ga中各受灾者每人分配3个物,于是得到对群组中各受灾者每人分配3个物
【r2=cl—rls3=15~3×5=0
资2即可,即:zll=z41一z5l=z7l=z8l=3,此时b1=0,厂ll=o,厂2l=o,^1=0.
同理,接下来进行物资2的分配,按照上述步骤,可得到物资2对群组G,和G2的受灾者分配结果为:z22一z。2—1,z。。=z。。一z。。2=3,此时f2—8>o,再转到对下一群组G3进行剩余物资分配.
检验q32s。一4×5=20>f。,此时剩余物资总量已经不能满足对群组G3中所有受灾者按其需求量分
配,同样由式(4)得到:j,.・2[詈]一[詈]一1
【r2=cz—r1s3—8—1×5=3
,于是得到对群组G3中各受灾者每人分配1个物
资2后,将会出现5名受灾者分配剩余3个物资2的情形,此时考虑群组G3中的所有5名受灾者的灾情程度,按某一种或者几种属性值综合灾情急需程度排序,将剩余的这3个物资2分配给排序最高的前3名受灾者即可.(注意:在同一群组内的受灾者,设受灾者i序号越小,其在同群组中的综合灾情急需程度越高.)即:在程序中,将剩余的3个物资2分配给受灾者1、4和5,结果为:z・z一础z—X52—2,z,2=X82=1,得到b2=0,f12—0,^2—0,^2—1.
如此下去,分别得到对物资3、物资4及物资5的分配结果及各自分配后的效用值如下:
z23
2z9323,z332X632
X103
24,X132
X43
2
23
X53。X732.,T83
得到fl。=0,厂23=0,^s=0,则
zz4
2勘4=2,z342z64。X10423,z142X442
X54
2
X74
2.T84。0
得到^。=0,厂2。=0,则
,T25
2
X95
24,z3524,-T652..T10523,.Z'152,3。452
3755
2
X75
2zss。0
得到^5—0;^5一1.
综上所述,可得到该算例最终的分配方案为:
32433
2132231113
3343343334
O23OO3O023
O440O3O043
X=(劫)10X5=
43324
从上述算例可以看出,此分配方案既保证了对救灾物资急需程度较高的受灾者优先得到物资,又可以在群组内部达到最大公平性,即分配结束后各群组内的各受灾者对物资的需求量和其所分配到的该物资数量之差的最大值与最小值的差额最小.
对于物资l,能够满足群组G,和G2的需求,但对于继续程度较低的群组G3则不能满足所有受灾者的需求,但剩余可分配救灾物资恰好能平均分配给该组的受灾者,使得^。=o;对于物资2,同样可以满足群组G。和G2的需求,但对于继续程度较低的群组G3也不能满足需求并且不能对该群组受灾者将剩余救灾物资平均分配,在约束条件的限制下,由题意得^。=1,最大程度地保证了群组内受灾者的公平;物资3的分配情况和物资1类似,最终^。=0;相比之下,物资4和物资5的可分配量较少,均达不到所需救灾物资
第1期谭秋月等:毁灭性灾难巾救灾物资分配方案的数学模型研究
!!!!!!!!!!!!!!11=:========================================m==Im!
・163・
的一半,按照群组优先的原则分配,只能分配给前两个继续程度相对较高的群组的受灾者,并且物资4恰好能将剩余可分配量满足群组G2的所有受灾者需求,得到^。一o;而物资5只能满足群组G2的部分受灾者,最终得到.厂25=1.
以上模型及数值算例的验证结果表明该模型可以在上述灾情状况下实现将运送到某一灾区的救灾物资向灾区受灾者及时合理的分配,此模型是合理、可行的.
5
模型的评价及改进
本文根据灾情状况建立了相应的数学模型,其通用性较强,所建立的数学模型能够简洁直观地刻画出
所要研究的问题,并用所能够描述的救灾者灾情、对救灾物资的急需程度和需求量的某几个属性值作为已知条件,在模型中体现救灾物资分配的详细原则.
对于毁灭性灾难救灾模型的改进,可以将该模型中的各受灾者针对不同的救灾物资所对应的灾情属性因素考虑得更全面一些,对不同的属性值在灾情中的影响给予一个合适的权值,然后将某几个属性综合加权之后将受灾者分成几个更合理的群组,而不是如算例1中所描述的那样只针对受灾者的年龄阶段将受灾者分成3个群组,这样会使救灾物资的分配方案模型更加合理化.
参考文献
[1]王杏.救灾物流中的物资调运模犁研究[D].北京:北京交通大学硕士学位论文,2007.[2]国家自然灾害救助应急预案IN].北京:新华社.2006-01一10.
[33王楠。刘勇。曾敏刚.自然灾害应急物流的物资分配策略研究[J].中国物流学术前沿报,2006,(11):518—525.[4]李阳。李聚轩。滕立新.大规模灾害救灾物流系统研究[J].科技导报.2005,25(7):64—67.[5]赵静,但琦.数学建模与数学实验(第2版)[M].北京:高等教育出版社。2003.[6]刁在筠,刘桂真。宿洁,等.运筹学(第三版)[M].北京:高等教育出版社。2007.[7]朱德通.最优化模型与实验[M].上海:同济大学出版社,2003.
MATHEMATICALMoDEL
OFDISTRIBUTIoN
METHoDoFRELIEFMATERIAL
TANQiu-yuel,SUNPing-anz,HANXi—eha03
(1.DepartmentofMathematics.WuyiUniversity,Wuyishan
cation
354300。China;2.InformationNetwork&Edu—
TechnologyCenter,WuyiUniversity,Wuyishan
354300,China;3.DepartmentofComputerScience
114051。China)
andEngineeringTechnology。LiaoningUniversityofScience{5LTechnology,Anshan
Abstract:Afterthenaturaldisasterhappened,becausetheamountofreliefmaterialthatiStransporting
to
a
certaindisaster
on
area
islimited,thepaperresearchesforthedifferentdistri-
butionschemesbased
thetwosituationsofdisasteraccordingtothecalamitydegree.Ca-
at
lamitydegree1istruculencewhilecalamitydegree2non-truculence.Aiming
gree
calamityde-
1,theauthor
sets
downthedistributionscheme,relevantdistributionmethodprinciple
anddistributionstrategybygroupclustering
method.What’Smore,a
to
correct
mathematical
modelandthemodel.
a
straightnumericalexample
are
given
validatetherationalityandfeasibilityof
Keywords:reliefmaterial;groupclusteringmethod;allocationmodel;numericalexample
毁灭性灾难中救灾物资分配方案的数学模型研究
作者:作者单位:
谭秋月, 孙平安, 韩希超, TAN Qiu-yue, SUN Ping-an, HAN Xi-chao
谭秋月,TAN Qiu-yue(武夷学院经济与数学系,福建,武夷山,354300), 孙平安,SUN Ping-an(武夷学院信息网络与教育技术中心,福建,武夷山,354300), 韩希超,HAN Xi-chao(辽宁科技大学计算机科学与工程学院,辽宁,鞍山,114051)
陕西科技大学学报(自然科学版)
JOURNAL OF SHAANXI UNIVERSITY OF SCIENCE AND TECHNOLOGY(NATURAL SCIENCE EDITION)2009,27(1)
刊名:英文刊名:年,卷(期):
参考文献(7条)
1. 朱德通 最优化模型与实验 20032. 刁在筠;刘桂真;宿洁 运筹学 20073. 赵静;但琦 数学建模与数学实验 2003
4. 李阳;李聚轩;滕立新 大规模灾害救灾物流系统研究[期刊论文]-科技导报 2005(07)5. 王楠;刘勇;曾敏刚 自然灾害应急物流的物资分配策略研究 2006(11)6. 国家自然灾害救助应急预案 2006
7. 王杏 救灾物流中的物资调运模型研究 2007
本文链接:http://d.g.wanfangdata.com.cn/Periodical_xbqgyxyxb200901036.aspx