地震灾后物资分配优化模型

1、2011年南京理工大学数学建模竞赛承诺书我们仔细阅读了全国大学生数学建模的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网 上咨询等）与本队以外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的 ，如果引用别人的成果或其他公开的 资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参 考文献中明确歹0出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规 则的行为，我们愿意承担由此引起的一切后果。我们的参赛（报名）队号为：17参赛组别（研究生或本科）：本科参赛队员（先打印，后

2、签名,并留联系电话）：打印签名联系电话队员1（队长）：张杰队员2:李晶队员3:胡雪瑶3编号专用页参赛队伍的参赛队号：（请各个参赛队提前填写好）：竞赛统一编号（由竞赛组委会送至评委团前编号）：竞赛评阅编号（由竞赛评委团评阅前进行编号）：题目地震灾后物资分配优化模型摘要本文研究了在地震发生后对各灾区进行物资合理分配的问题。针对问题我们首先利用层次分析法确定了各个灾区对各种物资的相对急需程度，记为an。其次，我们根据救灾物资分配使所有灾区整体效用最高的原则， 引入符号：为，一第i个灾区对所需的第j 种物资的实际分配量，Cj 第j种物资的可供分配总量，构建了以整体效用为目标函数 的整数规划模型及约束条

3、件：MAX W W aijXj S.T. Z Xj <Cj，然后利用LINDO软件进行数据处理与求解， ijJ在求解过程中我们考虑到每个灾区对每种物资的最低保障量Mij和满意度，丁是引入满意度系数，记为XijTij以及最低保障量的约束条件得到以整体满意度为目标函数的整数规划模型及约束条件:MAXIZi =1Jzj =1XijTijS.T. M ij < Xij < Tj其中T 第i个灾区对'1 二- aij Xij + ' 2 ''3i =1j =1i =1j =1T ij第j种物资的需求量，最后利用多目标规划法对两个目标函数进行优化：MAXS

4、.T.E XijM Cj , Mi产 X产 Tij。其中人 1，人2 i习为各单目标函数修正系数,旨在修正两目标函数因目标数值不在同一数量级上产生的差异。随后我们用汶川地震的相关数据对模型进行实际算例，证明了模型是可操作的。同时，我们还讨论了当应急物资重要性不同时的改进模型并对灾情发生变化或在物资分 发过程中存在的接受意愿、歧视等主观因素进行了定性的分析，使模型更加全面、科学、 合理。最后，我们对模型进行了优缺点总结及合理推广。关键词：物资分配、层次分析法、满意度系数、整数规划、多目标规划一、问题重述近年来地震频发，其预测目前比较困难但在灾后如能及时救援，可以很大程度减少 伤亡，其中救援物资分

5、配非常关键。但在物资分配中实际存在着诸多问题，例如：智利 在震后救援物资的分配上严重不均，日本大地震中已经躲过大地震和海啸的灾民因得不 到及时的物资分配而在避难所死亡等，因此研究地震灾害后的物资分配问题是具有现实 意义的，请从以下方面来考虑物资分配的问题：1、考虑灾区、受灾者和物资的不同，建立数学模型制定分配原则并给出合理的分 配方法。2、收集各类实际数据，给出一个符合题意的数值算例。3、通过以上分析，给出你的量化优化方案及建议。在考虑本本题时，还需注意一些I可题：1、受灾者的灾情不同，对生活物资的急需程度和需求量不同。而且各地的灾情在 不断发生变化，如何优化方案应对这种变化。2、生活物资应当

6、根据受灾区域大小、受灾程度、人口密度、灾区群众需求进行分 配，保证重点，确保及时、快捷、公开、高效、公平、公正发放。严禁物资发放中的优 亲厚友、性别歧视、年龄歧视和孤残歧视行为，在保障需求的同时，避免浪费。3、数值算例最好采用实际数据，且尽量提高数据容量。4、可以考虑线性规划、整数规划、多目标规划等各类优化模型。5、分配中有很多主观因素，例如：接受意愿、歧视等，这些因素能否考虑到模型 中。二、问题分析在救灾物资的分配过程中，考虑到物资的流动性与多样性等特点，影响救灾物资的 分配效果的因素是众多的，总的来说，灾区对物资的相对急需程度是非常重要的一个因 素，如果能够按照各灾区对各种物资的相对急需程

7、度来合理分配物资，那么救灾的效果 也就越好，因此如何确定各灾区对各物资的相对急需程度及如何构造效用函数是核心问18题。根据以上问题的特点我们用层次分析法确定物资相对急需程度，用整数规划法构造 目标函数建立了基础模型。但是我们通过数学分析和进一步的考虑发现了基本模型存在 着为达到目标函数最大化而牺牲了对物资相对急需程度较小的灾区的利益的问题，于是我们引进了各灾区对各种所需物资的最低保障量作为约束条件的下线。同时我们也意识到我们忽略了各灾区对各种物资分发的满意度的研究，因为满意度是物资的实际分配量与物资的需求量的比值，满意度系数越大，所得分配的满意程度也就越大，则救灾的效 果也就越好，即满意度系数

8、与整体效用成正相关，而这也是影响救灾效果的重要因素。基于上述分析，我们定义了整体满意度函数，显然，由于定义了两个目标函数使得要求 解的问题变成了一个多目标规划的问题，故需要将目标函数转化为新的目标向量。通过 以上分析我们就找到了基础模型建立的路径并在此基础上对模型进行了合理的优化，以下是我们研究思路图：整体效用目标函 数和整体满意度目 标函数各灾区对各类物资的相对急需程度的分析判定 n多目标函数规划及约束条件的确U定最低保障量 和物资总量 以及各灾区需求量LINDO进行优化求解AHP软件进行层次分析各灾区对各物资实际分 配量具体算例比较模型改进与推广优缺点评价三、模型假设1、所有参与分配的物资

9、都是灾区急需的重要物资，不同物资之间不可替代。2、在地震发生以后，信息是及时且完全对称的，即当地民政部门已经对灾情进行 了调查统计并分析评估出了基本数据o3、分配物资的供给是小丁灾区的需求的。4、为了方便模型建立，将非整数通过数据整数化处理转换为整数来考虑。5、救灾物资在运输的过程中损坏率为零。6、不考虑实际分发中主观因素的影响，并且忽略救灾物资在运往救灾点以及救灾 物资分配时所占用的时间。7、分配的物资只能由灾民个人享用，不存在物资的转移，且每个灾区之间相互独 立。四、符号说明团：第i个灾区对第j种物资的相对急需程度Aj : I个灾区对J种物资的相对急需程度矩阵为：第i个灾区的第j种物资的实

10、际分配量Tj ：第i个灾区的第j种物资的需求量Mj ：第i个灾区的第j种物资的最低保障量Cj :第j种物资可供分配量I:灾区总个数J:物资总类数Ei:物资分配的整体效用E2：物资分配的整体满意度乳i :效用目标函数的修正权重系数%2 :整体满意度目标函数的修正权重系数j ：各种物资相互问的重要性比较后第j种物资的重要指标五、模型的建立与求解准则层定义为影响该种物资一、用层次分析法确定I个灾区对J种物资相对急需程度 将目标层设为各个灾区对某一种物资的相对急需程度，急需程度的因素，方案层为各个受灾区。通过相互比较确定各灾区对某一种物资的相对 急需程度，然后依次利用层次分析法确定各灾区对其他种类的物

11、资的相对急需程度，最 终得到各灾区对各种物资的相对急需程度，记为 aij。1、构建灾区对某种物资急需程度评价体系。灾害发生后，需要对各个灾区物资急需程度进行评价，但是对每一种物资急需程度 进行评价的因素是不一样的，这就需要建立多个层次分析体系分别确立不同的物资急需 程度。在此我们将应急物资分为五类：饮用水、帐篷、食品、医疗用品、生活用品，将 分别影响这五种物资的因素制表如下所示：“是”代表影响因素，空格代表不是影响因素。影响因素饮用水帐篷医疗用品食品生活用品人口密度是是是是是水电等基础设施破坏程度是是是房屋倒塌数量是备灾能力是是是是是田地破坏亩数是2、比较准则层中的因素对目标层的影响。在评价各

12、灾区对各种物资的急需程度的问题 中分别两两比较各因素对某种物资急需程度的影响，得到判断矩阵：ai i& i1日3 1如分析各灾区对丁帐篷的相对急需程度时，根据影响它的因素，我们可得到准则层 的三个因素分别是人口密度，房屋倒塌数量，备灾能力。通过它们的两两比较得到这三 个因素对各灾区对丁帐篷的相对急需程度的判断矩阵。同理，我们也可得到其他因素对 每个灾区对丁其他物资的相对急需程度的判断矩阵。3、比较方案层对准则层的影响。根据各灾区反馈的受灾信息以及专家根据经验和数 据的评估歹0出各灾区对丁每个准则的成对比较矩阵：准则1准则2准则3比较具有不同性质的两个因素对丁上一层因素的影响时所采用的相

13、对尺度是由在与灾 害相关的领域的资深专家根据丰富的灾害评估经验及救员经验和在灾害发生第一时间 内对灾害进行灾情评估的基础上而提供的， 这些数据是具有科学性、客观性与合理性的3、对上述判断矩阵进行一致性比较。如果各个判断矩阵的随机一致性比率 CR有一个不满足小丁 0.1的条件，即不具有 可接受的一致性，则向各个专家反馈，进行判断矩阵的调整，直到各判断矩阵都具有可 接受的一致性。4、计算组合权向量。将各矩阵得到的归一化后的权向量相乘，可借助层次分析法软件yaahp辅助计算，得出第i个灾区对每种物资的相对急需程度 a” ai2, ai3, ai4 , a5 ,最后建立I个灾区对J种物资的急需程度的矩

14、阵 Aj。二、建立有约束条件的物资分配的总效用函数通过对问题的分析与理解并结合实际情况，我们确定了物资分配的基本原则：在救 灾物资有限的情况下，对救灾物资进行最有效的分配就是使灾区的救灾物资分配达到的 整体效用最高，这就需要我们根据各灾区对各种物资的急需程度来合理分配物资，构建 效用函数。模型I ：地震发生后灾区的数量记为I,救灾物资种类记为J,根据以上分析可得 总效用函数：i jEi=二 二 a j x j( 1)i =1 j =1同时由丁救灾物资是有限的，所以我们得到了一个约束条件即为各灾区对同一种物资的 实际分配量应小丁这种物资的总量：I云 xi j < C j(2)i =1再由上

15、面得到的物资分配的总效用程度为目标函数和约束条件的整数规划模型，即可得：MAX E1S.T.' xi j - C ji =1x/ 0(3)模型:在实际的分配过程中，会存在人为因素导致部分灾民分不到物资，例如： 抢劫、歧视等。而计算机在计算过程中也会因为要得到目标函数的最大值而牺牲部分灾 民的利益，同时，各种物资的供应相对丁需求总是不足的，为了达到对物资合理科学的 分配，以保障在物资分配上每个灾区的每种物资都有最低的满意程度，我们引入了每个 灾区对所需每种物资的最低保障量和各灾区对各种物资分配的满意度系数，其中最低保障量是由相关专家根据以往丰富的数据与实际救援经验且结合地震发生后的及时反

16、馈 而得出的数据，是具有科学性、合理性与可信性。满意度系数是由每个灾区所需的每种 物资分配到的实际量与需求量的比值来表示的，用它来衡量每个灾区对每种物资的满意 程度，从而定义整体满意度函数：I Jx iiE2=W W 广(4)i =1 j =1 1 ij乂由丁每个灾区每种物资的分配量不应小丁最底保障量，所以有约束方程：X/ M ij(5)根据模型假设3可得：Xi j ° T i j(6)模型m :效用函数El从物资供给者的角度定义了救灾效果的好坏，既仅仅根据各 灾区对各种物资的相对急需程度分配给各灾区相应的物资，并不考虑物资对各灾区的实际影响；而满意度函数 E2则从物资接受者地角度定

17、义了救灾效果的好坏，能够反映出 灾民对实际救灾效果的评价。我们认为仅仅考虑其中一方对丁评价整体救灾效果都是不 全面的，应将两函数结合起来引入多目标规划，即可得：MAXij1 '、'、i -1j -1IJ3ij xij + ' 2 -i =1 j -1x ijTij(7)S.T.I' X i ji =1EC j(3)x - ijM ij(5)x - xijTij(6)其中:i=（1,2,3.I）;j=（1,2,3 ;J）Xi , A值由实际灾情发生后aj和Tj的值确定。确定原则为消除因效用函数和满意度函数目标数值不在同一数量级上产生的差异对整个目标函数的影响，即修

18、正了权重系数使二个分目标函数对丁整个目标函数的权重相等。乂考虑各修正系数分别乘以的是各自的整个分目标函数，所以可以断定 4 ,九2并不会影响最终的优化分配方案。最后根据以上建立好的优化模型m,应用LINDO软件对数据进行处理，给出此模型的一个合理的最优化的分配方案。六、数值算例1、算例描述。2008年5月12日四川发生特大地震，由丁波及范围极为广泛，受灾地区众多，所 以导致统计数据巨大，因此我们只选取了 5个具有代表性的灾区进行实例分析。这 5个 地区每个地区有人员伤亡，以及房屋和各种基础设施的不同程度的破坏，对灾区人民的 生命和生活造成了极大的威胁。地震发生后救灾总指挥部对各灾区的受灾状况和

19、主要的 应急物资的需求量进行了统计，并且救灾物资是有限的。其中主要救灾物资有五种：饮 用水、帐篷、食品、药品、生活用品，并且对这五种物资进行编号依次为1、2、3、4、5。表1是各灾区对各种物资的需求量，表 2是各灾区对所需的各种物资的最低保障量。表1各灾区对各种物资的需求量物资种类汶川北川广元彭州雅安饮用水（万箱）1.020.890.540.390.27帐篷（万顶）1.230.720.440.200.11食品（万箱）1.321.120.650.130.08药品（万盒）0.970.870.540.230.10生活用品（万件）1.100.740.580.320.09表2各灾区对所需物资的最低保障量

20、物资种类汶川北川广元彭州雅安饮用水（万箱）0.810.630.340.150.10帐篷（万顶）0.930.540.300.100.05食品（万箱）1.010.900.520.070.05药品（万盒）0.750.580.280.110.04生活用品（万件）0.820.470.250.140.032、用层次分析法确定各灾区对各种物资的相对急需程度。首先确定各灾区对饮用水的相对急需程度。将目标层设为各灾区对饮用水的相对急 需程度，准则层为人口密度、水电设施破坏程度、备灾能力，方案层为5个受灾区域。通过两两比较人口密度、水电设施破坏程度、备灾能力相对丁各个灾区对饮用水的急需 程度的影响建立判断矩阵，通

21、过yaahp软件先得出判断矩阵的随机一致性指标CR经检验CR满足小丁 0.1的要求，因此认为它有可接受的一致性。进而得到各灾点对饮用水 的相对急需程度，得到aii =0.3814 0.2964 0.1527 0.0981 0.0714。然后运用相同的方法依次可得到各灾区对帐篷的相对急需程度ai2 =0.3011 0.2456 0.20230.1587 0.0923,对食品的相对急需程度ai3 =0.4207 0.2333 0.1503 0.12310.0726,对药品的相对急需程度 ai4 =0.4538 0.2483 0.1421 0.0944 0.0615, 对 生活用品的相对急需程度 血

22、=0.3926 0.3179 0.1711 0.0756 0.0428,从而得到5个灾区对5种物资的相对急需程度，记为矩阵:'0.38140.30110.42070.45380.3926、0.29640.24560.23330.24830.3179Aij =0.15270.20230.15030.14210.17110.09810.15870.12310.09440.075610.07140.09230.07260.06150.0428 ；根据Aij和上面确定的Tj ,为了消除数量级间的差异，经过分析与计算得到知=0.7,2 =0.3.3、 已知5种救灾物资的可供分配量。饮用水为 2.

23、5万箱，帐篷为2.35万顶，食品 为2.95万箱，药品为2.00万箱，生活用品为2.00万件。4、用构建模型中的函数求解该问题。将相关数据带入到所建模型的函数中应用LINDO软件求解最优分配方案。用软件实现的 5个灾区对5种物资的实际分配量如下表 所示：表3基于效用与满意度最大化的5种物资分配明细表物资种类汶川北川广元彭州雅安饮用水（万箱）1.150.650.370.120.06帐篷（万顶）0.960.830.430.170.11食品（万箱）1.230.990.600.080.05药品（万盒）0.880.650.310.120.04生活用品（万箱）0.920.600.290.160.035、结

24、果分析。从结果可以看出，各种物资的总分配量与可供分配量是相等的，并 且每个灾区所分配到的每种物资的实际量是不超过需求量的并保证了每个灾区对每种 物资的最低保障量，以上都说明了我们计算出来的分配量是满足模型的约束条件的。各 灾区的实际分配量是可以体现出各灾区对各种物资的相对急需程度与满意度的，例如： 雅安灾区对5种救灾物资的相对急需程度较低，在分配过程中就将相对较少的物资分配 给雅安灾区，而将大量物资分配给相对需求程度较大的灾区汶川、北川，避免了需求少 者多拿需求大者拿不够的局面的出现，能够保证效用与满意度的最优。同时，我们根据 得到的数据计算出满意度系数：表4各灾区对各种物资分配的满意度系数物

25、资种类汶川北川广元彭州雅安饮用水0.930.900.840.590.55帐篷0.920.930.800.420.43食品0.900.880.920.620.63药品0.910.750.570.520.40生活用品0.820.810.520.530.33由上表可看出大部分满意度系数是大丁 0.5的，只有极个别满意度系数并没达到 0.5, 这一方面说明就整体而言各灾区对丁所分配物资的满意程度是较高的，另一方面也说明了存在着牺牲需求量小的区域的满意度而满足整体满意度最大化的问题，因此模型在有些方面还是不太准确完备的，仍需要进一步改进。总体而言运用本模型是可以达到整体 救灾效果最优的结果，可以在一定程

26、度上解决灾后物资合理分配的问题，因此模型是合 理、实用且较优的。七、模型的改进1. 讨论当各种救灾物资相互重要性不一样时对模型的改进。引入n , R可根据灾情发生后，由丁众多原因造成的第j种物资的相对重要性指标。丁是对原目标函数进行改动得:MAXi j，寸 a*i =1 j =1rj1 J x.,、广 J j=lTij(8)约束方程(3) (5) (6)不变;模型的求解建议：由丁改进后的模型加入了一个衡量某种物资相对重要性指标的系数，且是以指数形式出现，所以使目标函数无法再用线性规划的思维考虑，因此我们建 议可以利用遗传算法对此非线性规划模型进行数据处理和求解。2、讨论对灾区的灾民分类保证各个

27、灾民的效用与满意度。我们所研究的模型主要针 对的是各个灾区对丁各种救灾物资的相对急需程度，然而当物资被分配到各灾区后并不 代表着灾区中的每位灾民都能够实现他们对各种物资的需要。因此，我们认为有必要引 入不同灾民对各种物资的相对急需程度的分析， 但是具体去分析每个灾民对各种物资的 相对急需程度的工作量是巨大的，且可操作性不强，故我们提出对灾民分类的建议，即 将具有同样届性的灾民归为一类，针对每一类灾民类似地建立相应的效用函数与满意度 函数，从个人角度衡量了救灾物资的分配是否合理。3、基于灾情进一步扩大的问题，我们认为是对所建模型的再次运用。因为灾情变化的实质我们认为就是各灾区对各种物资的相对急需

28、程度的变化，核心即为确定新的各灾区对各种物资的相对急需程度，因此我们可以通过层次分析法确定新的各灾区对各种 物资的相对急需程度，进而利用所建函数求出各灾区对各种物资的最优分配量，达到解 决问题的目的。八、模型评价与推广优点1、模型的建立基于实际情况，具有一定的实用价值。2、模型对问题研究合理、科学，理论性强。3、模型适用范围广，易于推广，例如：在经济生活中对于“如何分配有限的资源 使人们获得的最大的收益”此类问题同样适用。4、模型具有简洁性，层次分析法的基本原理和基本步骤易于理解，计算也相对简 便，容易为决策者了解和掌握。缺点1、模型中定义的对各种物资的急需程度的因素不够全面。2、模型对于救灾物资的分类过于笼统宽泛，应进一步细化。3、模型对于一些人为因素只是进行了简单的定性的分析或是一些泛泛的改进想法 而并未深入定量分析，显得不够具体准确。4、层次分析法的比较与判断是粗糙的且人为主观因素的影响很大，使得决策可能 难以被众人接受。九、建议基于对以上模型的讨论，我们想给物资分配相关部门和灾民提出一些意见与建议。当一些不可抗的自然灾害发生后灾民对各种物资的需求量总是远远大于相对有限的可分配量的，因此如何将有限的资源最合理的分配是非常重要的，要做到合理分配就 必须对各灾区的实际情况进行详细、准确的了解，此时根据从灾区反馈回来的第一手信 息对物资急需程度的评定