地震灾后的物资分配

1、精选优质文档-倾情为你奉上 地震灾后的物资分配摘要通过层次分析法和模糊综合判别的模型，经过分析灾区受灾情况和受灾人员，求得不同灾区所受灾害严重权重（i表示不同的受灾地区），从而确定不同灾区的物资需求量，然后由物资供应点与受灾地区构造距离矩阵（i表示不同物资供应点，j表示不同的受灾地区），建立规划模型，由约束条件得到供需向量，从而求解目标值规划问题。以汶川地震为例，搜集相关数据（见附录6），求解得到灾后物资优化分配，根据模型得到的结论，提出部分建议。一、问题重述近年来，我们生活的地球发生了多次大地震，虽然地震的预测目前比较困难，但如果在灾后能及时援救，可以很大程度减少伤亡，其中救援物资的分配非常

2、关键。在我国汶川大地震中，由于物资调配及时，在很大程度上降低了灾害的影响。为研究地震灾害后的物资分配，考虑以下问题：1. 考虑灾区、受灾者和物资等的不同，建立数学模型制定分配原则并给出合理的分配方法。2. 收集各类实际数据，给出一个符合题意的数值算例。3. 通过以上分析，给出你的量化优化方案及建议。二、模型假设1. 物资供应量能够满足灾区物资需求量。2. 物资供应点与灾区需求点距离按照两地之间的直线距离。3. 物资供应点分配方式道路运输不行时可以采取空中运输等方式。三、符号说明1.：不同灾区的受灾程度权重。2.：综合评价向量。3.：综合评价值权重向量。4.:个物资供应点。5.: 个受灾区。6.

3、：代表每个供应点有救灾物资。7.：代表每个受灾点需要救援物资。8.：物资供应点与受灾区之间的距离矩阵。9.：评价指标向量。10.：不同的准则；四、问题分析问题1：第一问要求考虑灾区、受灾者和物资等的不同，建立数学模型制定分配原则并给出合理的分配方法。不同灾区受灾严重情况不同，所以要建立模型：按照受灾严重情况求出不同地区的受灾程度，根据受灾程度来决定受灾物资分配多少。还要考虑到救灾物资如何分配，比如，共有m个物资供应点，n个受灾区，如何将m个物资供应点的物资高效的分配到灾区。剩下的问题就是解决将有限的救灾物资投入到应急救援最急需的地方，实现最优化配置问题。可以通过建立物资集散地与受灾区的距离矩阵

4、，转化成求解最短路径问题，从而达到救灾物资的最优分配。问题2：第二问要求通过收集各类实际数据，给出一个符合题意的数值算例。以汶川地震为例，收集相关数据，受灾区距震中的距离、人口密度、受伤人数、死亡人数及经济发展程度数据，根据第一问建立的分配模型来求解灾区物资的优化分配。问题3：通过第一问所建立的模型，以及第二问的数值解给出量化分配方案和建议。五、模型建立与求解问题1：通过评估不同灾区受灾的严重程度，从而受灾程度得到不同灾区物资分配多少，可以用两种方法求解受灾严重程度：模糊评判和层次分析。下面分别介绍两种方法求解过程。求解不同灾区救灾物资需求量：模型（1）模糊综合评判1、 评价对象受灾地区：下标

5、n代表有n个不同评价对象2、 评价指标受灾因素： 下标m代表m个不同的评价指标3、 权重系数权重向量： 其中4、 价值权重综合评价向量：模型求解:即：5、 综合评价值权重归一化模型（2）层次分析法目标层：O方案层： （下标n代表有n个不同评价对象。）准则层：设要比较各准则对目标O的重要性 要由A确定对O的权向量。考察完全一致的情况： 阶一致阵的唯一非零特征根为。求解公式： 对应的权向量(特征向量)：记第2层（准则层）对第1层（目标层）的权向量为：因考虑完全一致的情况，故不需进行一致性检验。同样求第3层(方案)对第2层每一元素(准则)的权向量。构造矩阵：则第3层对第1层的组合权向量建立物资供应点

6、如何向灾区分配物资：救灾物资分配就是将有限的救灾物资投入到应急救援最急需的地方，实现最优化配置问题。建立物资集散地与受灾区的距离矩阵，求解最短路径。模型建立：符号说明：m:m个物资供应点。n：n个受灾区。：代表每个供应点有救灾物资。：代表每个受灾点需要救援物资。：物资供应点与受灾区之间的距离矩阵。约束条件： 目标函数：目标函数S表示：合考虑物资供应点与受灾地区距离和灾区物资需求，在满足物资需求的情况下，按照最近邻分配，从而得到最优值S,得到的就是供应点到灾区的物资分配情况。问题2:以汶川地震为例，选取八个灾区，它们分别是：成都市、德阳市、绵阳市、广元市、阿坝州、雅安市、眉山市、巴中市。评价指标

7、选取受灾区距震中的距离、人口密度、受伤人数、死亡人数及经济发展程度，通过以上因素结合模糊综合评判和层次分析求解受灾程度。模糊综合评判根据汶川地震的受灾情况，我们选取受灾区成都市、德阳市、绵阳市、广元市、阿坝州、雅安市、眉山市、巴中市作为方案层，分别用表示；选取受灾区距震中的距离、人口密度、受伤人数、死亡人数及经济发展程度作为评价指标，分别用表示；权重向量为：。评价因素对应矩阵为：由公式得归一化后以上程序见附录一。分配优化：选取四川经济基础比较好的成都市、绵阳市、德阳市作为救援物资集结地点，将这三个地点的物资总量当做“1”，测出三个集结点距8个灾区的距离，根据归一化后的综合评价指标权重，运用目标

8、规划模型约束条件： 目标函数：求解供需分配矩阵。下面使用lingo进行线性规划求解供需矩阵。用lingo求解，程序见附录二，结果为：Global optimal solution found. Objective value: 90.21531 Infeasibilities: 0.E-16 Total solver iterations: 11 VOLUME( Y1, Z1) 0. 0. VOLUME( Y1, Z2) 0. 0. VOLUME( Y1, Z3) 0. 1. VOLUME( Y1, Z4) 0.E-01 0. VOLUME( Y1, Z5) 0. 0. VOLUME( Y1,

9、 Z6) 0.E-02 0. VOLUME( Y1, Z7) 0.E-02 0. VOLUME( Y1, Z8) 0.E-02 0. VOLUME( Y2, Z1) 0. 123.2000 VOLUME( Y2, Z2) 0. 0. VOLUME( Y2, Z3) 0. 2. VOLUME( Y2, Z4) 0. 0. VOLUME( Y2, Z5) 0. 62.20000 VOLUME( Y2, Z6) 0. 166.1000 VOLUME( Y2, Z7) 0. 123.2000 VOLUME( Y2, Z8) 0. 9. VOLUME( Y3, Z1) 0. 213.7000 VOLUM

10、E( Y3, Z2) 0. 90.90000 VOLUME( Y3, Z3) 0. 0. VOLUME( Y3, Z4) 0.E-01 0. VOLUME( Y3, Z5) 0. 120.3000 VOLUME( Y3, Z6) 0. 210.4000 VOLUME( Y3, Z7) 0. 213.9000 VOLUME( Y3, Z8) 0. 17.00000供应点-受灾点-供应量柱状图（matlab柱状图作图程序见附录3。）模糊综合判别距离矩阵物资需求供应量关系成都德阳绵阳广元阿坝雅安眉山巴中供应量成都061.6107.8260.3222.5125.266.2288.90.55德阳61.6

11、046.6199.2223.1229.7127.8236.30.24绵阳107.846.60154.4236.9229.7174.2200.00.21需求量0.1350.35640.12050.1870.18880.00680.00270.0028层次分析法根据汶川地震的受灾情况，我们选取受灾区成都市、德阳市、绵阳市、广元市、阿坝州、雅安市、眉山市、巴中市作为方案层，分别用表示；选取受灾区距震中的距离、人口密度、受伤人数、死亡人数及经济发展程度作为 评价指标，分别用表示。说明：求解方案层对准则层、准则层对目标层程序见附录4。 求得的准则层对目标层、方案层对准则层的一致性矩阵见附录5。准则层对目

12、标层权重方案层对目标层权重-0.2044 -0.7002 -0.1242 -0.3188 -0.0217 -0.1707 -0.5186 -0.4959 -0.7725 -0.1109 -0.2043 -0.2070 -0.6361 -0.0941 -0.1067 -0.4741 -0.1476 -0.1389 -0.4478 -0.3593 -0.2795 -0.0080 -0.5610 -0.3033 -0.6873 -0.3534 -0.0772 -0.0008 -0.0163 -0.4019 -0.3136 -0.3333 -0.0003 -0.0038 -0.2130 -0.6011

13、 -0.2430 -0.0003 -0.0031 -0.4092方案层对目标层权重归一化得分配优化：同样选取四川经济基础比较好的成都市、绵阳市、德阳市作为救援物资集结地点，将这三个地点的物资总量当做“1”，测出三个集结点距8个灾区的距离，根据归一化后的综合评价指标权重，运用目标规划模型约束条件： 目标函数：求解供需分配矩阵。下面使用lingo进行线性规划求解供需矩阵。用lingo求解，结果为：Global optimal solution found. Objective value: 103.2536 Infeasibilities: 0.E-16 Total solver iteratio

14、ns: 12VOLUME( Y1, Z1) 0. 0. VOLUME( Y1, Z2) 0. 0. VOLUME( Y1, Z3) 0. 1. VOLUME( Y1, Z4) 0.E-01 0. VOLUME( Y1, Z5) 0. 0. VOLUME( Y1, Z6) 0.E-01 0. VOLUME( Y1, Z7) 0.E-01 0. VOLUME( Y1, Z8) 0.E-01 0. VOLUME( Y2, Z1) 0. 122.7000 VOLUME( Y2, Z2) 0. 0. VOLUME( Y2, Z3) 0. 1. VOLUME( Y2, Z4) 0.E-01 0. VOLU

15、ME( Y2, Z5) 0. 61.70000 VOLUME( Y2, Z6) 0. 165.6000 VOLUME( Y2, Z7) 0. 122.7000 VOLUME( Y2, Z8) 0. 8. VOLUME( Y3, Z1) 0. 213.7000 VOLUME( Y3, Z2) 0. 91.40000 VOLUME( Y3, Z3) 0. 0. VOLUME( Y3, Z4) 0.E-01 0. VOLUME( Y3, Z5) 0. 120.3000 VOLUME( Y3, Z6) 0. 210.4000 VOLUME( Y3, Z7) 0. 213.9000 VOLUME( Y3

16、, Z8) 0. 17.00000供应点-受灾点-供应量柱状图层次分析法距离矩阵 物资需求 供应量关系成都德阳绵阳广元阿坝雅安眉山巴中供应量成都061.6107.8260.3222.5125.266.2288.90.55德阳61.6046.6199.2223.1229.7127.8236.30.24绵阳107.846.60154.4236.9229.7174.22000.21需求量0.12430.22990.12520.15180.17090.0560.0570.0849问题3模型求解的准确程度很大一部分取决于受灾地区受灾情况数据的多少，及时，准确，要迅速搜集各种种类的灾害信息，及时反馈到决策

17、中心，通过模型得到各个灾区的物资需求，安排不同的物资供应点调配物资供应受灾地区。六、模型评价与改进模型评价：1、 解决问题一时运用了两种模型，通过对比易知运用层次分析法的主观性相对较小，在求解物资需求时层次分析法相对较好。2、 建立的目标函数S，是按照最近邻物资分配条件得到的，通过求解目标函数S,达到求解物资供需的合理分配 不足与改进：1、 评估各个灾区物资需求量时，层次分析法较好，因为模糊综合评价在利用评价指标的权重时主观性相对于层次分析法有更大的主观性，比较两种方法求得的最终结果可以看出，层次分析给出了更好的物资分配方式。2、 该模型考虑的灾害因素太少，增加考虑的因素可以更好的确定各个灾区

18、的物资需求，从而很好的把握好物资调配量，同时也能更好地安排供应点与灾区物资的调配，最大程度地减少受灾程度。3、 模型中简化了距离因素，考虑的是直线距离，实际情况并不是如此，即使是空运也很难距离最短飞行，根据地形的不同，肯定会更复杂，道路曲折坎坷，道路受阻情况也十分普遍，这一点在实际求解模型时要给于考虑，使模型求解的结果更符合真实情况，也更趋合理。4、 现实过程当中，灾区信息根据时间的变化而变化，要建立一个动态物资分配模型，时时调整决策。七、参考文献1、地震应急救援物资最优分配问题初探 李晋, 袁志祥(陕西省地震局, 陕西西安 )2、1.C+面向对象程序设计教程 游洪跃 清华大学出版社3、数学建

19、模与数学实验（第三版）赵静 但琦 高等教育出版社4、MATLAB 7.0 实用教程 张圣勤 机械工业出版社八、附录附录1：模糊综合评判测评判向量 cd G:load 0.txt; w=0.25 0.20 0.30 0.35 -0.10'data=X0;sum=0;for n=1:8 y(n)=data(n,:)*w; sum=sum+y(n);end a=y'sum m=y'/sum附录2：用lingo求解：（模糊综合评判求解物资分配）model :sets: gongying/y1.y3/:capacity; xuqiu/z1.z8/:demand; links(go

20、ngying,xuqiu):volume,distance;endsetsmin=sum(links:distance*volume);for(xuqiu(J): sum(gongying(I):volume(I,J)=demand(J);for(gongying(I): sum(xuqiu(J):volume(I,J)=capacity(I);data : capacity=0.55 0.24 0.21; demand=0.135 0.3564 0.1205 0.187 0.1888 0.0068 0.0027 0.0028; distance=0 61.6 107.8 260.3 222.

21、5 125.2 66.2 288.9 61.6 0 46.6 199.2 223.1 229.7 127.8 236.3 107.8 46.6 0 154.4 236.9 229.7 174.2 200;enddataend附录3：柱状图作图程序：（模糊综合评判）function zhuzhuangtub = 0.1350.11640 0.09750.18880.00680.00270.00280 0.24 0 0 0 0 0 00 0 0.12050.08950 0 0 0;bar3(b);ylabel('供应点');xlabel('受灾点');zlabel(

22、'供应量');附录4：求一致矩阵程序：（方案层对准则层）#include<iostream>#include<fstream>#include<iomanip>using namespace std;int main()float b85; float a88;int i,j;int n;ifstream infile("1.txt");if(infile.fail()cout<<"da kai shi bai"<<endl;exit(1);while(!infile.eof()

23、for(i=0;i<8;i+)for(j=0;j<5;j+)infile>>bij;infile.close();ofstream outfile("2.txt"); if(outfile.fail()cout<<"da kai shi bai"<<endl;exit(2); for(n=0;n<5;n+)for(i=0;i<8;i+)for(j=0;j<8;j+)aij=bin/bjn;for(i=0;i<8;i+)for(j=0;j<8;j+)outfile<<

24、aij<<setiosflags(ios:scientific)<<setprecision(2)<<" "cout<<aij<<endl;outfile<<endl;outfile.close();system("pause");return 0;求一致矩阵（准则层对目标层）#include<iostream>#include<fstream>#include<iomanip>using namespace std;int main()float

25、 b51; float a55;int i,j;ifstream infile("1.txt");if(infile.fail()cout<<"da kai shi bai"<<endl;exit(1);while(!infile.eof()for(i=0;i<5;i+)for(j=0;j<1;j+)infile>>bij;infile.close();ofstream outfile("2.txt"); if(outfile.fail()cout<<"da kai

26、 shi bai"<<endl;exit(2); for(i=0;i<5;i+)for(j=0;j<5;j+)aij=bi0/bj0;for(i=0;i<5;i+)for(j=0;j<5;j+)outfile<<aij<<setiosflags(ios:scientific)<<setprecision(2)<<" "cout<<aij<<endl;outfile<<endl;outfile.close();system("pause&

27、quot;);return 0;附录5：准则层对目标层矩阵：1 1.40e+000 7.00e-001 5.38e-001 1.75e+000 7.14e-001 1.00e+000 5.00e-001 3.85e-001 1.25e+000 1.43e+000 2.00e+000 1.00e+000 7.69e-001 2.50e+000 1.86e+000 2.60e+000 1.30e+000 1.00e+000 3.25e+000 5.71e-001 8.00e-001 4.00e-001 3.08e-001 1.00e+000 方案层对准则层矩阵：方案层对震中距离：1 1.20e+0

28、00 1.00e+000 4.31e-001 7.32e-001 5.78e-001 6.52e-001 3.40e-001 8.35e-001 1.00e+000 8.36e-001 3.60e-001 6.11e-001 4.83e-001 5.44e-001 2.84e-001 9.99e-001 1.20e+000 1.00e+000 4.31e-001 7.31e-001 5.78e-001 6.51e-001 3.40e-001 2.32e+000 2.78e+000 2.32e+000 1.00e+000 1.70e+000 1.34e+000 1.51e+000 7.89e-0

29、01 1.37e+000 1.64e+000 1.37e+000 5.89e-001 1.00e+000 7.91e-001 8.90e-001 4.65e-001 1.73e+000 2.07e+000 1.73e+000 7.46e-001 1.26e+000 1.00e+000 1.13e+000 5.88e-001 1.53e+000 1.84e+000 1.54e+000 6.62e-001 1.12e+000 8.88e-001 1.00e+000 5.22e-001 2.94e+000 3.52e+000 2.94e+000 1.27e+000 2.15e+000 1.70e+0

30、00 1.92e+000 1.00e+000 方案层对人口密度：1.00e+000 1.35e+000 3.38e+000 4.74e+000 8.74e+001 9.07e+000 2.10e+000 2.88e+000 7.41e-001 1.00e+000 2.50e+000 3.51e+000 6.47e+001 6.72e+000 1.56e+000 2.13e+000 2.96e-001 3.99e-001 1.00e+000 1.40e+000 2.58e+001 2.68e+000 6.22e-001 8.52e-001 2.11e-001 2.85e-001 7.13e-00

31、1 1.00e+000 1.84e+001 1.91e+000 4.43e-001 6.07e-001 1.14e-002 1.54e-002 3.87e-002 5.43e-002 1.00e+000 1.04e-001 2.41e-002 3.30e-002 1.10e-001 1.49e-001 3.73e-001 5.23e-001 9.63e+000 1.00e+000 2.32e-001 3.18e-001 4.76e-001 6.42e-001 1.61e+000 2.26e+000 4.16e+001 4.31e+000 1.00e+000 1.37e+000 3.47e-00

32、1 4.68e-001 1.17e+000 1.65e+000 3.03e+001 3.15e+000 7.30e-001 1.00e+000 方案层对死亡人数：1.00e+000 2.50e-001 1.95e-001 8.95e-001 2.21e-001 1.53e+002 4.28e+002 4.28e+002 3.99e+000 1.00e+000 7.78e-001 3.57e+000 8.83e-001 6.10e+002 1.71e+003 1.71e+003 5.13e+000 1.28e+000 1.00e+000 4.59e+000 1.13e+000 7.83e+002

33、 2.19e+003 2.19e+003 1.12e+000 2.80e-001 2.18e-001 1.00e+000 2.47e-001 1.71e+002 4.78e+002 4.78e+002 4.52e+000 1.13e+000 8.81e-001 4.04e+000 1.00e+000 6.91e+002 1.93e+003 1.93e+003 6.55e-003 1.64e-003 1.28e-003 5.86e-003 1.45e-003 1.00e+000 2.80e+000 2.80e+000 2.34e-003 5.86e-004 4.56e-004 2.09e-003

34、 5.17e-004 3.57e-001 1.00e+000 1.00e+000 2.34e-003 5.86e-004 4.56e-004 2.09e-003 5.17e-004 3.57e-001 1.00e+000 1.00e+000 方案层对受伤人数：1.00e+000 4.13e-001 3.39e+000 7.12e-001 1.05e+000 1.96e+001 8.39e+001 1.02e+002 2.42e+000 1.00e+000 8.21e+000 1.73e+000 2.55e+000 4.73e+001 2.03e+002 2.48e+002 2.95e-001 1.22e-001 1.00e+000 2.10e-001 3.10e-001 5.77e+000 2.47e+001 3.02e+001 1.40e+000 5.80e-001 4.76e+000 1.00e+000 1.48e+000 2.74e+001 1.18e+002 1.44e+002 9.51e-001 3.93e-001 3.22e+000 6.78e-001 1.00e+000 1.86e+001 7.98e+001 9.74e+001 5