炼油厂的选址问题.doc

井冈山大学第五届“井冈杯”数学建模竞赛B 题参赛论文（油田选址问题）彭安平（数理学院 10 计算本一 ） （会员） 联系方式：郭本阳（数理学院 10 计算本一 ） 曾燕 （电信学院 09 信息本一） 1摘要 本文是对九个油井来选择最优的炼油厂地址使其总运输费用最少问题通过给出的所对应的油井口的产量分别分析 3 个问题，列出所求的目标函数建立数学模型，借助 matlab 来实现复杂的计算，从而得到较优的方案此题的根公式： 2sdkz问题一：选点的最优问题，通过综合考虑分别在九个井口处建炼油厂建立模型，由根公式并通过 matlab 实现计算比较得出最优方案（一号油井处建炼油厂，总费用为 853520*k） ；问题二：通过随机模拟在其范围内的各点建炼油厂，分别计算出所需要的总费用，由根公式通过 matlab 编程计算找出其最佳地址求出的炼油厂建在（34.95，41.51）km 处，总费用为 952103459250267/2147483648*k；问题三：将 9 个油井随机分为 2 组为两炼油厂分别负责的油井号。

由随机模拟的方法和根公式并通过 matlab 分别计算出 4 个方案中较优的方案，然后对比得到最终的较优方案，两厂坐标为(4.0460,80.7837)和(43.8511,29.8459)总费用为 66741342533787/268435456*k最后对本文所建立的模型在实际中进行了分析，并提出了改进方向问题一中假设建立货物虚设站，利用树算法，可以进一步优化总费用（见附录五简述） 关键词: 根公式 动态随机点模拟 matlab 编程计算 虚设站2一、 问题重述在商品高度市场化的当今社会，物品时时刻刻处于流动之中资料显示，商品的平均物流成本占总成本的 36%，而其制造成本仅占总成本的 13%据2011 年 5 月 9 日 CCTV《经济半小时》报道，去年我国物流总费用占 GDP 的18%，比发达国家高出近一倍，放在市场竞争日益激烈的今天，物流管理显得日益重要某一油田在一平坦地区拥有九口油井，其年产量及位置如下表所示所有的原油都需要运输到炼油厂进行提炼现在不考虑炼油厂的建设费用，因此总的费用仅与炼油厂的位置有关在假定的单位运费与运输距离成正比的条件下，需对以下问题做出决策：一、 如果两点间的距离以折线计算，且九个井口均可作为炼油厂的候选位置，问炼油厂建在哪个井口附近（该井口到炼油厂距离以 0计）最佳，总运输费用是多少？二、 若两点间距离以直线距离计算，且该区域的任一点均可作为炼油厂的候选厂址，炼油厂应建在何处，总费用是多少？三、 若油田高层已决定在该地区见两个炼油厂，不考虑炼油厂的建造费用，仅考虑运费，两个炼油厂分别建在什么位置，各应服务于哪几个油井（假定一个油井的原油只能运往一个炼油厂） ，才能使总运费最低，总费用是多少？井号 位置（X，Y）(km) 产量（万吨）1 （22，38 ） 172 （8，13 ） 403 （4，81 ） 604 （51，32 ） 205 （38，11 ） 256 （17，12 ） 157 （81，63 ） 508 （19，45 ） 89 （62，12 ） 30请分别建立以上三个问题的数学模型，并予以求解，对你所建模型的优劣性进行评估。

3二、问题分析问题—：要求炼油厂建设在所给的 9 个油井口附近（不计厂与此油井的距离） ，以其中一个油井为中心分别向其它井作折线分别计算出以 1-9 号油井为炼油厂所需的总运输费用，从而得出较好的方案问题二：因为炼油厂可以建在所给的区域的任一点，此时利用 MATLAB 的随机函数给出随机点作为炼油厂厂址再分别计算此随机厂与 9 口油井的总运输费用经过大量随机点，就能得出较优的方案问题三：现要求建 2 个炼油厂，所以可将那 9 个油井口分成 2 组：(1,8)、(2,7) 、(3,6)、(4,5)再由大量随机点分别计算 4 种情况所需要的总运输费用，最后比较得出最优解以上问题均用 matlab 实现三、模型的假设1、当炼油厂建在某个井口附近，则不考虑该井口到炼油厂的距离；2、题中我所给的平坦地区假设为 ;]81,0[,][kmykx3、假设运输车严格按照路线行驶且炼油厂占地为零4、假设单位运费与运输距离的比例系数为 k5、假设随机点能取遍整个可行区域6、不考虑炼油厂和虚设站的建设费用四、 符号说明X 各油井横坐标的数组Y 各油井纵坐标的数组Z 各油井产量的数组m 单位运费，d 两点之间的距离，k 单位费用与运输距离的比例系数，z 产油量，s 单条路线费用S 总费用M 炼油厂 MN 炼油厂 N号油井的横坐标iX4号油井的纵坐标iY号油井的产量Zi第 次随机时油厂的横坐标ix第 次随机时油厂的纵坐标yi第 次随机时炼油厂 M 的横坐标ix1x2(i) 第 次随机时炼油厂 N 的横坐标iy1(i) 第 次随机时油厂 M 的纵坐标y2(i) 第 次随机时油厂 N 的纵坐标五、模型的建立与求解问题一： 模型建立：1） 、由单位运费与运输距离成正比可得：mkd.(1)又因为 sdz .(2) 由（1）和（2）可得 2k……（3） 。

所以以各口油井且以折线建炼油厂的各总运费公式为： 9 21()(|)(|)(|()jSikXjiYjiZj（i=1,2…..9）……(4)比较 S(i)，得出最终运费 S=min(S(i))模型求解：整理得下表：（用 matlab 编程求解，源程序见附录二）井号 位置（X，Y ）Km产量（万吨） 在 i 号井旁建炼油厂的总运费S(i)1 （22,38） 17 853520*k2 （8,13） 40 1302869*k3 （4,81） 60 1782125*k4 （51,32） 20 1031265*k55 （38,11） 25 1244205*k6 （17,12） 15 1225997*k7 （81,63） 50 1963232*k8 （19,45） 8 854100*k9 （62,12） 30 1519220*k由上表可知炼油厂应建在 1 号油井处，此时的总运输费用为 853520*k坐标图示和源程序见附录一（1） ）：问题二模型建立：1） 、在坐标一定范围内随机产生一点：（x(i),y(i)）作为炼油厂的坐标点，同理由根公式 2(3)sdkzL可得在这点建厂且以直线运输所需的运输费用公式：9221()()()jSiXjxiYjiZj2） 、一定条件下利用 MATLAB 软件进行随机模拟尽可能曾大随机数 ，使运输i费用 S 越趋向于所有可能。

再找出其中最小的 S(i)，和相应的炼油厂坐标iyx,模型求解：（源程序见附录三）通过 MTLAB 程序得ans =34.94967107903618ans =41.51413800062496ans =952103459250267/2147483648*k所以炼油厂建在（34.95，41.51）km 处，如下图所示且此时的总费用为952103459250267/2147483648*k其坐标图和坐标图源程序见附录一（2）问题三：模型建立：将 9 个油井随机分为 2 组为两炼油厂分别负责的油井号，由随机模拟的方法得6出模拟炼油厂并计算出各组中两模拟炼油厂所负责的油井的总费用，取各次模拟总费用的最小值，得出各组模拟的最小值，最后得到一个最小费用的分组并用 matlab 输出此费用的两厂坐标值即为最终所求模型求解：（源程序见附录四）整理得下表：M 负责1口，N 负责8口(h=1)M 负责2口，N 负责7口(h=2)M 负责3口，N 负责6口(h=3)M 负责4口 N 负责5口(h=4)M 负责8口，N 负责1口（h=8） （验证数据）M 负责的油井号 3 3,7 3,7,8 3,4,7,9 1,2,4,5,6,7,8,9N 负责的油井号 1,2,4,5,6,7,8,9 1,2,4,5,6,8,9 1,2,4,5,6,9 1,2,5,6,8 3M(4.0460,80.7837) M(38.9410,72.7700) M(37.5563,71.0205) M(44.7150,56.0694) M(43.9905,29.8865)两炼油厂的坐标 N(43.8511,29.8459) N(31.6494,19.2309) N(32.4835,17.8656) N(19.3524,18.8838) N(4.0187,80.9354)总运输费用 2.4863e+005*k 3.4003e+005*k 4.1257e+005*k 4.0207e+005*k 2.4864e+005*k由上表可知：建两个炼油厂时的较优方案为一个厂负责 3 号油井的加工，另一个则负责其余8 口油井的加工。

其坐标图和坐标图源程序见附录一（3）六、模型评价（1）优点：1.计算过程采用 MATLAB 数学软件，使计算比较精确简便，减少误差 ;2.合理的假设，使复杂问题简单化，抽象问题具体化;3.在模型分析中，通过多种方法对比计算，更容易找到我们的优化解决方案;4.运用了一些图形，用数形结合法来进行分析，使模型思路更清晰，更有说服力;5.本文建立模型符合实际情况，真实的考虑了各种车型、运费单价等等情况对总费用最省的影响，并且具有一般性，便于模型的实施、应用和推广;6.本文用的数学方法都比较简单易懂，方便方案的利用7.我们采用了随机性，更加具有普遍性，这样更能接近真实值2）缺点：1.虽然用了数学软件 MATLAB 编程计算，但由于自身电脑性能的限制，不能进行更加精确的计算，其中不免有些误差2.模型假设的时候可能考虑不周，难免存在一些细小问题被忽略；3.语言方面表达可能有所欠缺，影响论文7七、参考文献井冈山大学数学建模协会，MTLAB 与数学建模，第一版重庆大学数学系 数学实验 科学出版社 第一版原思聪，MATLAB 语言机械工程应用，机械工程出版社，200 年 8 月第一版优酷视频 matlab 基础教程， 8 4 51 38 17 81 19 62];B=[38 13 81 32 11 12 63 45 12];plot(A,B,'pentagram')title('油井位置 ')gtext('1 号(22,38)')gtext('2 号(8,13)')gtext('3 号(4,81)')gtext('4 号(51,32)')gtext('5 号(38,11)')gtext('6 号(17,12)')gtext('7 号(81,63)')gtext('8 号(19,45)')gtext('9 号(62,12)')gtext('炼油厂')8（2）问题二：A=[22 8 4 51 38 17 81 19 62 34.95];B=[38 13 81 32 11 12 63 45 12 41.51];plot(A,B,'pentagram')title('油井位置 ')gtext('1 号(22,38)')gtext('2 号(8,13)')gtext('3 号(4,81)')gtext('4 号(51,32)')gtext('5 号(38,11)')gtext('6 号(17,12)')gtext('7 号(81,63)')gtext('8 号(19,45)')gtext('9 号(62,12)')gtext('炼油厂（ 34.95，41.51）')9（。