2974.抗旱优化方案

1、抗旱优化方案摘要：抗旱方案的优化，旨在使我们在花费最小的情况下获得最大的回报。因此，我们进行了大量的数据实验，运用excel表格，matlb，lingo等软件对我们的方案进行优化，并且运用0-1整数规划对我们的选井进行优化选择方案。手心我们考虑了三种方案的可行性，通过分析我们发现一次只打井和只铺设管道军是不可行的；于是只有采用第三种方案同时打井与铺设管道施工。在第三种方案的基础上，我们有进一步将其优化为一次性打井和逐年分次打井两套方案。在建立模型上，首先我们计算出了老四井在2010至2014年的每年产水量总和；之后进行优选打井方案，当然，管道的铺设使我们必不可少的考虑因素。最后，我们对两种方案

2、经性对比，取的最优的方案，即为逐年打井更能为这一旱区提供更多的适量保障。关键词：0-1规划 线性方程 优化设计 拟合一、问题的重述我国地域辽阔，水资源分布不均。使得我国部分地区处于半干旱或干旱状态，为了更好的解决干旱地区群众的生活用水问题。我们以一个处于我国西南地区的某个偏远贫困村为例，进行我们的分析解剖。我们所的此地年平均降水量不足20mm，是典型的缺水地区。过去村民的日常生活和农业生产用水来自于，一是每家每户自行建造的小蓄水池，用来屯积雨水，二是利用村里现有的四口水井。由于近年来环境破坏，经常是一连数月滴雨不下，这些小蓄水池的功能完全丧失。而现有的四口水井经过多年使用后，年产水量也在逐渐减

3、少，在附表1中我们知道了其9年来的产水量粗略统计数据。2009年以来，由于水井的水远远不能满足需要，不仅各种农业生产全部停止，而且大量的村民每天要被迫翻山越岭到相隔十几里外去背水来维持日常生活。 为此，政府打算着手帮助该村解决用水难的问题。从两方面考虑，一是地质专家经过勘察，在该村附近又找到了8个可供打井的位置，它们的地质构造不同，因而每个位置打井的费用和预计的年产水量也不同，详见表2，而且预计每口水井的年产水量还会以平均每年10%左右的速率减少。二是从长远考虑，可以通过铺设管道的办法从相隔20公里外的地方把河水引入该村。铺设管道的费用为（万元），其中表示每年的可供水量（万吨/年），表示管道长

4、度（公里）。铺设管道从开工到完成需要三年时间，且每年投资铺设管道的费用为万元的整数倍。要求完成之后，每年能够通过管道至少提供100万吨水。政府从2010年开始，连续三年，每年最多可提供60万元用于该村打井和铺设管道，为了保证该村从2010至2014年这五年间每年分别能至少获得150、160、170、180、190万吨水，请作出一个从2010年起三年的打井和铺设管道计划，以使整个计划的总开支尽量节省。二、问题的分析针对这一现状，我们发现在这三年内只铺设管道或只打井均是不符合现状要求的，因为其不仅达不到村民在2010年之后的用水需求而且不符合三年内铺设完管道的长远利益要求；故而，我们采取了铺设管道

5、与打井的同时施工的措施方案。我们根据具体情况建立了二种模型；分别是一次性打井和逐年分次打井两套模型方案。对此，我们首先运用matlab软件拟合出了老四井随时间的产水量变化趋势，并获得了其在2010至2014年的产水量。然后，依据题设中村民各时间段用水总量的限制要求；分别在两种模型中进行优化求解获得最小花费。而且，我们发现铺设管道的费用在要求的输水量情况下，其花费是一个固定的值；因此，为节约花费，我们需要在选打井口数以及各井的产水量和打井费用上行优化选择，是的总打井费用最低。最终我们获得：结余=总投资钱数-铺设管道用费-打尽付费用，获得最大结余。三、模型的假设与符号说明模型假设1、假设不考虑小蓄

6、水池的作用和利息的因素。2、假设该村水的供水仅有水井和管道提供。3、假设在表2的每口水井的年产水量会以固定的每年10%速率减少。4、假设不考虑其他因素的影响，铺设管道三年内完成且第4年开始供水。5、假设在所打井中实际供水量不低于预测的供水量。符号说明xi- 表示选取i号井，i=（1，2，3，4，5，6，7，8）qij- 表示第i口井第j年的产水量，i=（1，2，3，4）.j=（2010，2011，2012，2013，2014）t- -表示时间，（单位：年）q -表示管道输水量qi-表示第i年打井的产水量，i=（1，2，3）j -表示结余费用m -表示总开支m1- 表示铺设管道总开支m2- 表示

7、打井所需总开支q旧n-表示旧井在n年的总产量，n=（2010,2011,2012,2013,2014）q新n- 表示新井在n年的总产量，n=（2010,2011,2012）z -表示固有资金量四、模型的建立与求解我们首先考虑到这个村庄的老四井问题；即这个村庄以前的四口井预期的产水量的变化情况。我们依据matlb拟合出其每口井产水量的变化趋势，并预测其2010至2014年每年每口井的总产水量。4.1 1的拟合图从上拟合图，我们可万一看出1的产水量q1出于逐年递减的状态。并从matlab求解中我们得到了其产水量随时间变化的线性方程，即为： qij =-1.2017*t+33.4639；并且进一步求

8、解出了该井在2010至2014年的产水量q11、q12、q13、q14、q15，分别为：21.4469、20.2452 、19.0435 、17.8418 、16.6401。4.2 2的拟合图同样的，我们从2的2001至2009年的产水量变化的三点图中，我们发现他合乎一个一元三次方程的坐标分布规律，因此我们得出了如上的拟合图。并且获得了该井产水量随时间的变化方程即为：qij =-0.0311*t.3+0.7966*t.2-7.5771*t+28.2357；并且进一步求出了其2010至2014年的预期产水量q21、q22、q23、q24、q25分别为1.0247、0、0、0、0. 4.3 3拟合

9、图 从3的三点图，我们发现其2001至2009年产水量随时间的变化基本符合一次线性函数，故而，我们拟合并求解除了其产水量随时间的变化方程，即为：qij =-2.0850t+29.9806.并进一步求出了该井2010至2014年的预期产水量q31、q32、q33、q34、q35，分别为：9.1306、7.0456 、4.9606、2.8756、0.7906。4.4 4的拟合图由4的散点图分布，我们可以看出他基本符合产水量随时间的一元三次方程变化规律曲线。我们拟合出了上图，并且求出了他的产水量随时间变化的一元三次方程的系数，即该方程为：qij =-0.1773*t.3+3.7333*t.2-27.

10、0442*t+86.2095，并求出了该井在2010至2014年的预期产水量q41、q42、q43、q44、q45，分别为：11.7975、4.4663、0、0、0。4.5 从上面的老四井中的预测中，我们可以得到这四口老井在2010至2014年，每年总共能给该村提供的水量q旧n。如下表04-1：表04-4 老井产水量表年份1234q旧n201021.44691.02479.130611.797543.3997201120.245207.04564.466331.7571201219.043504.9606024.0041201317.841802.8756020.7174201416.6401

11、00.7906017.4307模型一 一次性打井 为了满足，2010至2014年该村的需水要求，即该村在2010至2014年的需水分别为150、160、170、180、190万吨。我们建立了这样一个模型，即一次性打井能够为其后的2010、2011、2012年的需水提供足够的保障。根据如下的8个打井位置和所需费用：表2 8个位置打井费用（万元）和当年产水量（万吨）编号12345678打井费用57546553当年产水2536321531282212我们发现要是得一次打井满足以后三年的需求量，就必须优先考虑其对第三年，也就是对2012年的供水量是否满足，我们知道2012年最低需求供水总量为170万吨

12、，那么所打新井的供水量q新3=170- q旧3=145.9959（万吨），由此我们可以得出所打新井在2010年的供水量即为q新1= q新3/（1-10%）2,得出：q新1=180.2419，也就是说，我们所打新井在2010你的供水量不得小于180.2419.那么我们该在这8口井中如何选择才能使得花费最后，得到最优解呢？我们通过运用0-1整数规划，以及lingo软件；获得了我们的最优打井数和选井方案。即为：model:min=5*x1+7*x2+5*x3+4*x4+6*x5+5*x6+5*x7+3*x8;x1+x2+x3+x4+x5+x6+x7+x8=7;25*x1+36*x2+32*x3+15

13、*x4+31*x5+28*x6+22*x7+12*x8=180.2419从运算的结果中我们发现必须打7口井才能满足要求，而且打井所需花费为36万元。这几口井分别是第1、2、3、5、6、7、8号井；求解过程见附录五。对于每年的铺设管道费用和结余费用见我们整理后的表04-2和04-3 表04-2 花费细目表年份要求（万吨）老四（万吨）新六井（万吨）供应水（万吨）打井费用（万元）管道费用（万元）铺设管道（公里）结余花费（万元）201015043.3997186229.39936243.4735201116031.7571167.4199.1570608.682201217024.0041150.66

14、174.6640557.958201318020.7174135.594156.311000201419017.4307122.0346139.465000提供100顿水管道的费用139即m1=139，m2=36，m=m1+m2=175，结余j=z-m=180-175=5.表04-3 花费细目表要求（万吨）年份老四（万吨）新六井（万吨）供应水（万吨）打井费用（万元）管道费用（万元）铺设管道（公里）结余花费（万元）150201043.3997186229.399736243.354915440160201131.7571167.4199.15710608.38728862170201224.00

15、41150.66174.66410608.38728862180201320.7174135.594156.3114000190201417.4307122.0346139.4653000提供水量（q）107铺设管道花费144铺设总公里数20.12949我们从04-2和04-3的对比中我们发现，随着我们的管道提供水量的增加，我们的结余相应减少。且当结余0时，管道最大提供水量q=107万吨。同时，由于2013年开始管道为该村提供最少不小于100万吨的水量，再加上经所提供的2013和2014年分别为156.3114和139.4653万吨，所以说他满足，该村2013和2014 年分别180和190万

16、吨水的需求。模型二 逐年分次打井一次性打井，一次性投入过大；我们采取分次打井，分流资金；希望可以能减少总投资量。我们设其第一年所打的井的产水量为q1，第二年以及第三年所打的井的产水量分别为q2和q3。那么只需求出满足2010至2012年的每年所需最低水量150、160.170的q1、q2、q3的最优打井方案，即为q1+ q2+ q3取之最小。那么这三年新井总的最小产水量为w1，他应满足如下的方程目标函数：w1= q1-+ q2+ q3;106.6003= q1;st 128.2429= q2+ q1*0.9;145.9959= q3+ q2*0.9+ q1*0.81;381=107;最后我们所

17、得结果为选择第1、2、6、7号井作为2010年的打井方案，且其满足了新井的最低要求107万吨水，花费最低为22万。总体上该方案的花费总额为打井33万元，而铺设管道随着需要供应的水量有所变化。如下表04-4，当余额j=8万元时，其供水量q=100万吨，当余额j=0时，每年可以提供q=112万吨水。这可以有我们的附excel表中所得，随着输入提供水量的变化，可以得到不同的结余额。表04-4 花费细目表年份要求（万吨）老四（万吨）新六井（万吨）供应水（万吨）打井费用（万元）管道费用（万元）铺设管道（公里）结余花费（万元）201015043.3997107150.399722385.311949482

18、01116031.7571128.3160.05715557.688347201217024.0041146.47170.47416466.430254201318020.7174131.823152.5404000201419017.4307118.6407136.0714000提供水量100铺设管道花费139铺设总公里数20.43055由图标04-4中，2013和2014年的几口井产水量总和分别为131.823万吨和118.6407万吨。我们再加上管道提供的水量，已完全满足该村2013和2014年的180万吨和190万吨水量的需求。故而，有m1=139，m2=33，结余j=z-m1-m2=

19、8该模型可以说是成立的。 对方案一和方案二的对比，我们从长远的考虑来说方案二更优，相比较来说其比第一方案更结余花费；在以后的供水中，方案二的管道能够提供的水量也最大。方案二的管道能提供最大112万吨水而方案一的管道运输最大只能提供107万吨水；而且从资金的周转上方案二也更优，因为他太能提供分流资金的作用。故而，我们的从经济和长远发展来说，我们应该优选方案二。五、模型的评价与推广模型评价：优点：1该模型简单，思路清晰，使人一目了然。2利用matlab和lingo软件使数据更加准确。3该模型铺设管道和打井同时应用，在保证水量供给正常的情况下，使得总费用最少，做到了节约。缺点：该模型中没有雨水的积蓄

20、，雨水资源得不到充分的利用。该模型中没有考虑到经济过快发展和人口剧增所带来的用水量过大。模型推广：因为模型没雨水的积蓄量，雨水资源得不到充分的利用。在生活和生产上，应该合理利用积蓄的雨水进行生产，合理的调配水源，减少生活用水用于生产上，而又能使雨水积蓄的水量得到利用，保证生产的正常。该模型还可以石油管道的铺设、道路的修建等。六、软件应用附录1t=1,2,3,4,5,6,7,8,9;q=32.2 ,31.3 ,29.7 ,28.6 ,27.5 ,26.1 ,25.3 ,23.7, 22.7;aa=polyfit(t,q,1);a=aa(1)b=aa(2)，r=polyval(aa,t);a =

21、-1.2017，b = 33.4639 t=10;11;12;13;14q=-1.2017*t+33.4639t =10，11，12，13， 14q = 21.4469， 20.2452， 19.0435，17.8418，16.6401附录2 t=1,2,3,4,5,6,7,8,9;q=21.5 15.9 11.8 8.7 6.5 4.8 3.5 2.6 2.0;aa=polyfit(t,q,3);a=aa(1)，b=aa(2)，c=aa(3)，d=aa(4)，r=polyval(aa,t);plot(t,q,k+,t,r,r);a = -0.0311，b = 0.7966，c = -7.57

22、71，d =28.2357 t=10;11;12;13;14;q=-0.0311*t.3+0.7966*t.2-7.5771*t+28.2357q = 1.0247，-0.1179，-1.7199，-3.9679，-7.0485附录3 t=1,2,3,4,5,6,7,8,9;q=27.9 25.8 23.8 21.6 19.5 17.4 15.5 13.3 11.2;aa=polyfit(t,q,1);a=aa(1)，b=aa(2)，r=polyval(aa,t);plot(t,q,k+,t,r,r);a =-2.0850，b =29.9806 t=10;11;12;13;14;q=-2.08

23、50*t+29.9806q = 9.1306，7.0456，4.9606，2.8756，0.7906附录4t=2,3,4,5,6,7,8,9;q=46.2 32.6 26.7 23.0 20.0 18.9 17.5 16.3;aa=polyfit(t,q,3);a=aa(1)，b=aa(2)，c=aa(3)，d=aa(4)r=polyval(aa,t);，plot(t,q,k+,t,r,r);a = -0.1773，b = 3.7333，c = -27.0442，d = 86.2095 t=10;11;12;13;14;q=-0.1773*t.3+3.7333*t.2-27.0442*t+86

24、.2095q = 11.7975，4.4663，-7.1001，-23.9655，-47.1937附录5model:min= q1+q2+q3;106.6003=q1;128.2429=q2+q1*0.9;145.9959=381;q10;q20;q30;gin(q1);gin(q2);gin(q3);endglobal optimal solution found. objective value: 170.0000 extended solver steps: 0 total solver iterations: 2 variable value reduced cost q1 107.0

25、000 1.000000 q2 32.00000 1.000000 q3 31.00000 1.000000 row slack or surplus dual price 1 170.0000 -1.000000 2 0.3997000 0.000000 3 0.5710000e-01 0.000000 4 0.4741000 0.000000 5 0.7700000 0.000000 6 107.0000 0.000000 7 32.00000 0.000000 8 31.00000 0.000000附录6model:min=5*x1+7*x2+4*x4+5*x6+5*x7+3*x8;x1

26、+x2+x4+x6+x7+x8=4;25*x1+36*x2+15*x4+28*x6+22*x7+12*x8107;bin(x1);bin(x2);bin(x4);bin(x6);bin(x7);bin(x8);endglobal optimal solution found. objective value: 22.00000 extended solver steps: 0 total solver iterations: 0 variable value reduced cost x1 1.000000 5.000000 x2 1.000000 7.000000 x4 0.000000 4

27、.000000 x6 1.000000 5.000000 x7 1.000000 5.000000 x8 0.000000 3.000000 row slack or surplus dual price 1 22.00000 -1.000000 2 0.000000 0.000000 3 4.000000 0.000000六、参考文献1数学建模方法与范例，寿纪麟等编，西安交通大学出版社，1993。2 数学模型基础，王树禾编著，中国科学技术大学出版社，1996。3 数学模型选谈(走向数学从书)，华罗庚，王元著，王克译，湖南教育出版社，1991。4 数学模型与数学建模，刘来福、曾文艺编，北京师范大学出版杜，1997。5 姜启源.谢金星.数学建模m.北京：高等教育出版社，2003年.6 应用数学建模 applied mathematical modelling (elsevier)