南水北调工程水指标的分配问题

南水北调工程水指标的分配问题摘要南水北调中线工程建成后，有限的水资源需要我们合理分配利用，本文根据2000年的统计数据，用Excel软件计算处理数据，计算出每个城市具体分配的水量。针对问题一，对生活用水、工业用水和服务综合业的用水的分配，我们首先按2000年的用水标准对各城市进行初始分配。对于剩余的水，采用简化后的公平席位分配Q值法进行再分配。针对问题二，我们将所调用的工业用水和综合服务业用水产值总和作为目标函数，建立以总产值最大为目标的线性规划模型，限制各城市用于工业增加值和综合服务业的调水量在2000年平均值的50％~150%的围.利用matlab程序求得城市获得最大经济效益的调水分配，并求得最大收益为7080亿元。关键词：南水北调ExcelQ值法线性规划经济效益1、问题重述南水北调中线工程建成后，预计2010年年调水量为110亿立米，主要用来解决京、津、冀、豫四省（市）的沿线20个大中城市的生活用水、工业用水和综合服务业的用水，分配比例分别为40%、38%、22%。这样可以改善我国中部地区的生态环境和投资环境，推动经济发展。用水指标的分配总原则是：改善区域的缺水状况、提高城市的生活水平、促进经济发展、提高用水效益、改善城市环境。根据2000年的统计数据，各城市的人口数量差异大，基本状况和经济情况也不相同，缺水程度也不同（如下）。要研究的问题是：（1）、请你综合考虑各种情况，给出2010年每个城市的调水分配指标，使得各城市的总用水情况尽量均衡；（2）、由于各城市的基本状况和自然条件不同，对相同的供水量所产生的经济效益不同，请从经济效益的角度，给出调水指标的分配案。但是，要注意到每个城市的工业和综合服务业的发展受产业规模的限制，不可能在短时间无限制的增长。2、问题分析本题属于分配问题，分配的关键在于要适当照顾各城市经济发展的均衡，在此基础上再对各城市进行生活用水、工业和综合服务业用水的分配。对于生活用水的分配，我们首先要考虑保障人民的生活用水，先按2000年人均用水量标准进行初始分配。对初始分配有剩余的情况，为了体现公平原则，采取公平席位分配中的Q值法进行再分配。但这里的Q值法要进行简化处理。对工业用水和综合服务业用水的分配，。我们可以采用类似的法进行分配。对于问题二，既要考虑经济利益的最大化，又要适当照顾各城市经济发展的均衡，因此我们采用单目标线性规划的法，将所调的工业用水和综合服务业用水产值总和作为目标函数，使其达到最大即可。约束条件中，工业用水和综合服务业用水总量均有限制，更重要的是要考虑各城市的均衡发展，在这里我们可以从2000年各城市用于工业增加值和综合服务业的调水量的平均值出发，限定该值得围在50％~150%，进而求得最大经济效益下各城市调水分配。3、模型假设1）原有各城市的供水量不变，运水过程中水量未曾损耗2）各城市的经济相互独立，互不影响3）不考虑水的循环利用4）城市人口，工业产值及综合服务业年均产值增产值不变4、符号说明2000年第i个城市的人口数2000年第i个城市的人均生活用水量2010年第i个城市生活用水初始分配量5、模型的建立与求解5.1问题一模型准备5.1.1Q值法介绍设A,B两人数分别为；分别占有和个席位,则两每个席位所代表的人数分别为和.则相对公平的席位分配案为：如果两分别占有和席,当总席位增加1席时,应该分配给A还是B.不妨设,即对A不公平,当再分配一个席位时,有以下三种情况：I.当时,这说明即使给A增加1席,仍然对A不公平,所以这一席显然应给A.II.当时,这说明给A增加1席,变为对B不公平,此时对B的相对不公平值为:（2）III.当时,这说明给B增加1席,将对A不公平,此时对A的相对不公平值为:（3）因为公平分配席位的原则是使相对不公平值尽可能小,所以如果（4）则这1席给A,反之这1席给B.由(2)(3)可知,(4)等价于（5）不难上述的第I种情况也与(5)式等价,于是我们的结论是当(5)式成立时,增加的1席应给A,反之给B.若记:则增加的1席给Q值大的一.5.1.2Q值法简化应用说明表示2010年第i个城市的总人口数，表示2000第i个城市的人均用水量，所以第i个城市的Q值为：剩余水对第i个城市的分配量=剩余水量5.2问题一模型的建立5.2.1生活用水分配模型针对问题一，在分配供水时应充分满足用户的需求，故先从2000年新增人口生活用水新增总量为初始分配量，初步分配后发现有剩余，我们将剩余的水进行合理公平的分配。因此可采用席位分配模型，我们选取简化的Q值分配模型，根据生活用水初始分配量=城市新增人口数*人均用水量，则有：先以2000年各城市人均生活用水为标准进行初始分配，计算结果如下表 城市序号 1 2 3 4 5 6 7 8 9 10 初始分配（亿立米） 0.3418 0.1598 0.0478 0.0625 0.0187 0.0319 0.1178 0.0144 0.0255 0. 城市序号 11 12 13 14 15 16 17 18 19 20 初始分配（亿立米） 0. 0.0163 0.026 0.0834 0.069 0.0346 0.0347 0.0128 0.0147 0.0541由上表可知初始分配的水的总量为1.2597亿立米。由于110亿立米的水中有40%用于生活用水，即44亿立米，所以经初始分配后剩余的水有42.7403亿立米。根据简化Q值模型公式计算各城市的Q值如下表： 城市序号 1 2 3 4 5 6 7 8 9 10 Q值 103.0185 84.912 0.4705 0.605 0.5217 1.62845 5.754 0.8156 1.0021 0.5761 城市序号 11 12 13 14 15 16 17 18 19 20 Q值 0.3208 0.4707 1.7135 2.2587 14.96 1.6064 2.8871 0.321 1.2636 3.03再根据模型计算剩余水的分配如下表所示： 城市序号 1 2 3 4 5 6 7 8 9 10 剩余水分配1 19.3154 15.9205 0.0882 0.1134 0.0978 0.3053 1.0451 0.1529 0.1879 0.108 城市序号 11 12 13 14 15 16 17 18 19 20 剩余用水 0.0602 0.0882 0.3213 0.4235 2.8049 0.3012 0.5413 0.0602 0.2369 0.568所以各城市生活用水分配如下所示： 城市序号 1 2 3 4 5 6 7 8 9 10 生活用水 19.6572 16.0803 0.136 0.176 0.1165 0.3372 1.1629 0.1673 0.2134 0.174 城市序号 11 12 13 14 15 16 17 18 19 20 生活用水 0.0881 0.1046 0.3472 0.5069 2.8739 0.3359 0.576 0. 0.2516 0.6221同理，对于工业用水及综合服务业用水的分配我们可以采用类似的法进行分配，最终得到各城市三种用水的分配案如下表所示 生活用水量 工业用水量 综合服务用水量（亿立米） 序号 城市名称 （亿立米） （亿立米） 1 北京 19.6572 4.4855 13.631 2 天津 16.0803 19.8058 7.6323 3 0.136 0.4956 0.4132 4 0.176 1.6905 0.1044 5 0.1165 2.5408 0.0639 6 0.3372 0.4053 0.0439 7 1.1629 4.1766 0.7276 8 0.1673 0.3394 0.0573 9 0.2134 4.5311 0.1023 10 0.174 0.0275 0.0628 11 0.0881 0.0039 0.0823 12 0.1046 0.0357 0.0557 13 0.3472 0.0136 0.0818 14 0.5069 0.0134 0.0521 15 2.8739 2.8894 0.8957 16 0.3359 0.0201 0.0432 17 0.576 0.0407 0.0436 18 0. 0.0356 0.0307 19 0.2516 0.0201 0.0511 20 0.6221 0.205 0.02475.3问题二的模型建立5.3.1工业用水及综合服务业用水的分配模型对于问题二，我们要从经济效益的角度考虑工业用水和综合服务业用水的分配，同时注意到每个城市工业和综合服务业发展受到产业规模限制，不可能在短时间无限制地增长，故对每个城市的调水指标都应有上下限的约束，我们指定实际用水可调度量可在该城市用水需求量的50％~150%的围波动，使得各大城市在工业和综合服务业所产生的经济效益最大,于是以工业和服务业总经济效益最大为目标得到目标函数为从而得到总经济效益最大的线性规划数学模型为：S,t根据已有的条件限制，我们利用matlab软件计算，得到如下的分配结果（表）（程序见附录）： 　 　 生活用水量 工业用水量 综合服务用水量（亿立米） 序号 城市名称 （亿立米） （亿立米） 1 北京 19.6572 2.6247 5.1769 2 天津 16.0803 28.54 2.7126 3 0.136 1.0546 1.2138 4 0.176 1.6517 0.2432 5 0.1165 2.37 1.0203 6 0.3372 1.0217 1.1054 7 1.1629 2.7976 0.7708 8 0.1673 0.9402 0.431 9 0.2134 2.4927 1.0326 10 0.174 0.801 0.6893 11 0.0881 0.7739 0.57 12 0.1046 0.8254 0.6719 13 0.3472 0.7761 1.0616 14 0.5069 0.8153 0.2927 15 2.8739 3.0299 0.9847 16 0.3359 0.8089 0.7051 17 0.576 0.806 0.6398 18 0. 0.8058 0.3647 19 0.2516 0.7933 0.5693 20 0.6221 0.9308 1.0524求得2010年各城市获得工业经济效益为：995亿元，服务经济效益为：6085亿元，，并求得最大收益为7080亿元。6、模型的推广与评价6.1模型的推广本题属于规划问题，建立的是单目标线性规划模型，这个模型不仅仅适用于水的资源配置问题，它对规划类问题的求解都可以起到指导作用。通过资源配置最优化为杠杆平衡它们之间的分配关系。决策者要通过概念抽象、关系分析可将各类影响因子放入规划模型中，可以通过相关的计算机软件得到兼顾全局的最优解。本模型涉及预测，资源调度分配及合理优化使效益最大等问题，可在于发电厂的电力分配，公司生产销售及资源调度优化配置等领域有着广泛的应用。6.2模型的评价1)在生活用水公平分配模型中我们采用了简化的Q值法，对初始分配后剩余的水按Q值比例进行分配,达到了公平的效果。但模型中缺乏对简化Q值法合理性的讨论，这使得简化法应用的根据不足。2)在工业用水和综合服务业用水的分配模型中，我们采用了单目标线性规划的法。并且充分利用了有限的水资源创造最大的经济效益，不过在分配水资源时没有考虑到每个城市之间的公平性问题。7、参考文献[1]启源，金星，叶俊.数学建模.北京：高等教育出版社，2010.[2]刁在筠，桂真，宿洁，马建华.运筹学.高等教育出版社.2010.8.附录8.1分配给工业用水的最优解的程序c=[-1/143;-1/72;-1/102;-1/96;-1/110;-1/120;-1/86;-1/131;-1/126;-1/;-1/210;-1/;-1/205;-1/;-1/88;-1/210;-1/;-1/210;-1/200;-1/];r2=0.236088;v1=4.4855;v2=19.8057;v3=0.4956;v4=1.6904;v5=2.5407;v6=0.4053;v7=4.1766;v8=0.3394;v9=4.5311;v10=0.0275;v11=0.0039;v12=0.0357;v13=0.0136;v14=0.0134;v15=2.8894;v16=0.0201;v17=0.0407;v18=0.0356;v19=0.0201;v20=0.2050;a=[10000000000000000000;01000000000000000000;00100000000000000000;00010000000000000000;00001000000000000000;00000100000000000000;00000010000000000000;00000001000000000000;00000000100000000000;00000000010000000000;00000000001000000000;00000000000100000000;00000000000010000000;00000000000001000000;00000000000000100000;00000000000000010000;00000000000000001000;00000000000000000100;00000000000000000010;00000000000000000001;-10000000000000000000;0-1000000000000000000;00-100000000000000000;000-10000000000000000;0000-1000000000000000;00000-100000000000000;000000-10000000000000;0000000-100000000000000000000-100000000000;000000000-10000000000;0000000000-1000000000;00000000000-100000000;000000000000-10000000;0000000000000-1000000;00000000000000-100000;000000000000000-10000;0000000000000000-1000;00000000000000000-100;000000000000000000-10;0000000000000000000-1;];b=[1.5*r2*v1;1.5*r2*v2;1.5*r2*v3;1.5*r2*v4;1.5*r2*v5;1.5*r2*v6;1.5*r2*v7;1.5*r2*v8;1.5*r2*v9;1.5*r2*v10;1.5*r2*v11;1.5*r2*v12;1.5*r2*v13;1.5*r2*v14;1.5*r2*v15;1.5*r2*v16;1.5*r2*v17;1.5*r2*v18;1.5*r2*v19;1.5*r2*v20;-0.5*r2*v1;-0.5*r2*v2;-0.5*r2*v3;-0.5*r2*v4;-0.5*r2*v5;-0.5*r2*v6;-0.5*r2*v7;-0.5*r2*v8;-0.5*r2*v9;-0.5*r2*v10;-0.5*r2*v11;-0.5*r2*v12;-0.5*r2*v13;-0.5*r2*v14;-0.5*r2*v15;-0.5*r2*v16;-0.5*r2*v17;-0.5*r2*v18;-0.5*r2*v19;-0.5*r2*v20];aeq=[11111111111111111111];beq=[41.8];[x,fval,exitflag,output]=linprog(c,a,b,aeq,beq,[],[])8.2分配给综合服务业用水的最优解的程序c=[-1/160;-1/140;-1/;-1/360;-1/315;-1/318;-1/235;-1/315;-1/320;-1/310;-1/320;-1/352;-1/280;-1/310;-1/220;-1/320;-1/310;-1/340;-1/280;-1/320];v1=13.6310;v2=7.6323;v3=0.4132;v4=0.1044;v5=0.0637;v6=0.0439;v7=0.7276;v8=0.0573;v9=0.1023;v10=0.0628;v11=0.0823;v12=0.0557;v13=0.0818;v14=0.0521;v15=0.0897;v16=0.0432;v17=0.0436;v18=0.0308;v19=0.0511;v20=0.0248;a=[10000000000000000000;01000000000000000000;00100000000000000000;00010000000000000000;00001000000000000000;00000100000000000000;00000010000000000000;00000001000000000000;00000000100000000000;00000000010000000000;00000000001000000000;00000000000100000000;00000000000010000000;00000000000001000000;00000000000000100000;00000000000000010000;00000000000000001000;00000000000000000100;00000000000000000010;00000000000000000001;-10000000000000000000;0-1000000000000000000;00-100000000000000000;000-10000000000000000;0000-1000000000000000;00000-100000000000000;000000-10000000000000;0000000-100000000000000000000-100000000000;000000000-10000000000;0000000000-1000000000;00000000000-100000000;000000000000-10000000;0000000000000-1000000;00000000000000-100000;000000000000000-10000;0000000000000000-1000;00000000000000000-100;000000000000000000-10;0000000000000000000-1;];b=[1.5*r1*v1;1.5*r1*v2;1.5*r1*v3;1.5*r1*v4;1.5*r1*v5;1.5*r1*v6;1.5*r1*v7;1.5*r1*v8;1.5*r1*v9;1.5*r1*v10;1.5*r1*v11;1.5*r1*v12;1.5*r1*v13;1.5*r1*v14;1.5*r1*v15;1.5*r1*v16;1.5*r1*v17;1.5*r1*v18;1.5*r1*v19;1.5*r1*v20;-0.5*r1*v1;-0.5*r1*v2;-0.5*r1*v3;-0.5*r1*v4;-0.5*r1*v5;-0.5*r1*v6;-0.5*r1*v7;-0.5*r1*v8;-0.5*r1*v9;-0.5*r1*v10;-0.5*r1*v11;-0.5*r1*v12;-0.5*r1*v13;-0.5*r1*v14;-0.5*r1*v15;-0.5*r1*v16;-0.5*r1*v17;-0.5*r1*v18;-0.5*r1*v19;-0.5*r1*v20];aeq=[11111111111111111111];beq=[24.2];[x,fval,exitflag,output]=linprog(c,a,b,aeq,beq,[],[])8.32000年的基本情况 序号 城市名称 城市人口 工业产值 综合服务业总产值 人均生活用水量（升/日） 万元综合服务业用水量（m3） 万元工业增加值用水量（m3） 总数（万人） 年自然增长率（%） 增加值（亿元） 年增长率（%） 人均产值（万元） 年增长率（%） 1 北京 1285 2.04 737 11.1 1.16 13.2 354 160 143 2 天津 682 3.03 739 11.7 0.83 12.2 209 140 72 3 56 9.15 193 10.0 0.30 10.0 245 102 4 87 5.90 268 12.5 0.23 12.1 325 360 96 5 46 5.87 480 9.8 0.22 8.6 315 110 6 78 6.12 256 7.6 0.20 7.6 318 120 7 218 5.41 464 10.2 0.44 11.8 267 235 86 8 52 4.50 10.9 0.15 12.0 165 315 131 9 81 3.69 721 10.0 0.22 11.4 230 320 126 10 83 6.61 110 6.9 0.16 11.7 320 310 11 42 8.00 36 9.6 0.18 12.3 220 320 210 12 41 6.10 97 8.7 0.14 13.5 174 352 13 72 6.01 104 10.3 0.22 8.9 160 280 205 14 128 6.92 67 8.0 0.18 9.4 250 310 15 220 5.12 310 12.9 0.53 10.2 164 220 88 16 78 6.56 72 11.1 0.17 9.2 320 210 17 90 6.61 114 8.8 0.18 8.4 155 310 18 32 6.44 106 10.0 0.12 11.3 165 340 210 19 58 4.60 83 9.0 0.15 10.3 148 280 200 20 121 5.90 211 10.4 0.13 8.8 202 320 全国平均值 - 10.7 - 9.9 0.23 7.8 219 610 288专业资料_1641650083.unknown_1641650084.unknown_1641650085.unknown_1641650082.unknown