2010年数学建模c题输油管的布置解析

1、12012 高教社杯全国大学生数学建模竞赛输油管的布置摘要能源的运输线路关系到国家的经济发展，本文根据问题的条件和要求，针对两炼油 厂到铁路线距离和两炼油厂间距离的各种不同情形建立最优化模型。通过分析，将炼油 厂、车站、铁路线之间的距离作为未知常量，列出费用优化模型，完整地解决了问题。针对第一问：首先画出两炼油厂及车站的位置关系图，通过对问题的分析，在位置 关系图的基础上采用分步设计的思路，设计出了输油管道及车站的通用方案图。利用通 用方案图，设定能够表示非共用管道交汇点位置及火车站建设点位置的变量x、y，依据几何知识建立费用最小方案模型：W=p(J(a-y)7 + J(b y)2+(c-x)

2、2) + y，利用lingo软件编写程序，从而求解出任意情况下的费用最小方案。针对问题二：首先分析三家公司对附加费用的不同预测及自身的资质，我们采用加权平均的方法计算出合理的附加费用法，再由第一问的模型建立最优化模型：W = R (J(x2+(a-y)2+ J(b -d - y)2+(c-x)2) + P?y +卩3孙+(|_c)2通过ling软件编程从而求解出设计方案，该方案计算的费用为283.20万。方案如图所示：针对问题三：首先比较第三问与第二问，得出第三问与第二问的区别在于输油管道 费用不再是固定的值。改进第二问中的模型，建立第三问的最优化模型：min W = R!+ R2L2+ By

3、 + P3L3=PHE+Q-y)2+ P2j(b_d_y)2+(c_x)2+卩2丫+皿2+ (l -c)2代入数据从而得出了最优方案。方案计算的费用为252.47万关键词：lin go最优化模型加权平均值一.问题重述1问题的重述某油田计划在铁路线一侧建造两家炼油厂，同时在铁路线上增建一个车站，用来运 送成品油。由于这种模式具有一定的普遍性，油田设计院希望建立管线建设费用最省的 一般数学模型与方法。2提出问题：2（1）针对两炼油厂到铁路线距离和两炼油厂间距离的各种不同情形，提出设计方 案。在方案设计时，若有共用管线，应考虑共用管线费用与非共用管线费用相同或不同 的情形。（2）设计院目前需对一更为

4、复杂的情形进行具体的设计。两炼油厂的具体位置由附图所示，其中A厂位于郊区（图中的I区域），B厂位于城区（图中的II区域），两 个区域的分界线用图中的虚线表示。图中各字母表示的距离（单位：千米）分别为a = 5,b = 8，c = 15，l = 20。ILI)若所有管线的铺设费用均为每千米7.2万元。铺设在城区的管线还需增加拆迁和工程补偿等附加费用，为对此项附加费用进行估计，聘请三家工程咨询公司（其中公司一具有甲级资质，公司二和公司三具有乙级资质）进行了估算。估算结果如下表所示:工程咨询公司公司一公司二公司三附加费用（万元/千米）212420请为设计院给出管线布置方案及相应的费用。（3）在该实际

5、问题中，为进一步节省费用，可以根据炼油厂的生产能力，选用相适应的油管。这时的管线铺设费用将分别降为输送A厂成品油的每千米5.6万元，输送B厂成品油的每千米6.0万元，共用管线费用为每千米7.2万元，拆迁等附加费用同上。 然后给出管线最佳布置方案及相应的费用。二.问题分析（1）针对问题一，由题意可知，此问未提供任何与解题有关的已知的数据，仅给 出有两个炼油厂需要铺设通往火车站的管道，火车站也是未知待定的。要想设计出合理 的方案，就需要画出草图，以此形象的显示设计方案所涉及的不同情形， 未知数据先用 字母表示出来。 设计方案的总体思路是采用分步设计，首先架设非公用管道，其次架设 共用管道，最后确定

6、火车站的位置，。根据设计过程中总结的规律，建立最优化模型， 解决问题一提出的问题。（2）针对问题二，通过比较问题一和问题二的题设条件可知，问题二给出了影响总费用的因素，每千米的铺设费用7.2万元和附加费用。因三家工程咨询公司评估的附 加费用不同（其3中公司一21万元/千米，公司二24万元/千米,公司三20万元/千米）， 并且资质也不同（其中公司一具有甲级资质，公司二和公司三具有乙级资质）。所以我们采用加权平均的方法对评估出来的三种不同费用进行处理，求出更为合理的附加费用值。在设计费用最优化模型时，考虑到城市的拆迁费用，我们将输油管穿过城区的部分 用L3表示。然后分共线和不共线两种情况进行比较分

7、析：当A厂,B厂输油管有共用管线时，从B厂出发架设到共用管线时的非共用管线分两部分考虑L2,L3,共用管线用丫表示。未知变量均用字母表示出来，建立最优化模型，在lingo中输入目标函数与约束 条件，导出结果目标函数值；当A厂,B厂到车站的输油管没用共用管线时，由模型可知， 此时Y=0b同理，运用lingo软件，在共线模型的基础上给约束条件中的丫赋值为0， 计算出此时的目标函数值。通过比较A厂，B厂到车站是否使用共用管线的两种情况所 需费用，得出最优方案。（3）针对问题三：问题三在问题二的基础上考虑到实际问题，依据炼油厂的生产能力不同，选用相适应的输油管。这时的管线铺设费用将分别降为：输送A厂成

8、品油的每千米5.6万元，输送B厂成品油的每千米6.0万元，共用管线费用为每千米7.2万元.根据问题二共线和不共线两种情况下建立的模型，修改在lingo程序里的每千米的管道 费用值，即可导出结果，通过比较A厂，B厂到车站是否共用管线时的费用，确定最佳二.问题假设1、铺设的管道所经区域均为平整无阻隔的平地，不存在影响管道铺设的地形，如 河流，山谷等。2、铺设管道的路线由若干条直线段构成，不存在曲线。3、铁路线不存在拐弯或曲线的情况，可看做是一条直线四、名词解释及符号说明1.名词解释：1.2.2.符号说明W方案总费用R非共用管线单位长度的费用P2共用管线单位长度的费用p1由炼油厂A导出的非共用管线的

9、费用p2由炼油厂B导出的非共用管线的费用4P3公司i评估的附加费用P三家公司评估的加权平均值国i公司i的权重a炼油厂A到铁路线的垂直距离b炼油厂B到铁路线的垂直距离c炼油厂A、B在以铁路线为水平线的竖直距离Li由炼油厂A导出的非共用管线的长L2由炼油厂B导出的非共用管线的长1L.建模前共用管线的长的表示y建模后共用管线的长的表示L3在U区域内的非共用管线的长T非共用管线的交汇点D车站的建设位置x车站距离炼油厂A在铁路线上的垂点之间的距离五、 模型的建立及求解针对问题一：由题知本文是根据两个炼油厂的之间的位置及与铁路线的距离设定方案，方案中要确定输油管道的铺设路线及火车站的位置。方案要使铺设管道

10、的费用最省。问题一中未提供任何数据，说明解决第一问要解决两炼油厂在任意位置下的费用最优化 问题。由题意知铺设管线分为有共用管线和无共用管线两种情况，下面就对这两种情况 单独讨论：（1）有共用管线情况下的最优方案：分析问题知，要想费用最省的方案，就要知晓这几个数据：两炼油厂距离铁路线各自的垂直距离；两炼油厂以铁路线为水平线的水平距离，单位长度共用管线的费用R及非共用管线的费用P2。如图一：5A、B为两个炼油厂，a、b分别为A、B两个炼油厂距离铁路线的垂直距离，c为以铁 路线为水平线的水平距离。由于第一问中未提供任何数据，所以先把这些变量看作是常 量。知道这几个数据的值我们才能进一步设计方案。对于

11、有共用管线的情况下确定的方案，一定包含有图二中的信息：图二旷11AJL3.rL/y1D车站-C-图中设非公用管线和公共管线单位长度的费用分别为R、P2，Li、L2分别表示从A、B炼油厂导出的非共用管线的长度，T点表示两非共用管线的交汇处，从T点导出公用 管线连接到车站D，L表示共用管线的长度。在设计方案的时候我们先确定了T点，如图三：图三再从T点导出公用管线连接到车站D，从而确定L,由于铁路线上任一点均可作为 火车站，所以有无数条线可作为L，如图四：图四611AF1i、bDID2LU12车站CD3D4由于点到直线的垂线段是该点到直线上任意一点的线段的最短直线，所以当D点为点T在铁路线上的垂点时

12、L最小，共用管线的铺设费用最低。所以我们在设计输油管道 此时我们重新定义y表示L。图五旷1A丿- x=simple(x)x =1/2*(p2A2*b-4*p1A2*b-p2A2*a+4*p1A2*a+p2*c*(4*p1A2-p2A2)A(1/2)/p2/(4*p1A2-p2A2F(1/2)-1/2*(p2A2*b-4*p1A2*b-p2A2*a+4*p1A2*a-p2*c*(4*p1A2-p2A2)A(1/2)/p2/(4*p1A2-p2A2)A(1/2) y=simple(y)y =1/2*(4*p1A2*a-p2A2*a-p2A2*b+4*p1A2*b-p2*c*(4*p1A2-p2A2

13、)A(1/2)/(4*p1A2-p2A2)1/2*(4*p1A2*a-p2A2*a-p2A2*b+4*p1A2*b+p2*c*(4*p1A2-p2A2)A(1/2)/(4*p1A2-p2A2)syms a b c x yW=(a-y)A2+xA2)A(1/2)+(b-y)A2+(c-x)A2)A(1/2)+y; jacobia n( W,x,y) ans =1/(a-y)A2+xA2)A(1/2)*x+1/2/(b-y)A2+(c-x)A2)A(1/2)*(-2*c+2*x),1/2/(a-y)A2+xA2)A(1/2)*(2*y-2*a)+1/2/(b-y)A2+(c-x)A2)A(1/2)

14、*(2*y-2*b)+120syms a b c x y f1=(1/(a-y)A2+xA2)A(1/2)*x+1/2/(b-y)A2+(c-x)A2)A(1/2)*(-2*c+2*x)=0); f2=(1/2/(a-y)A2+xA2)A(1/2)*(2*y-2*a)+1/2/(b-y)A2+(c-x)A2)A(1/2)*(2*y-2*b)+仁O); x,y=solve(f1,f2,x,y)x =1/2*(-3*bA2+6*a*b+cA2-3*aA2-3*(-b+a+1/3*3A(1/2)*c)*b+3*a*(-b+a+1/3*3A(1/2)*c)/c1/2*(-3*bA2+6*a*b+cA2-3*aA2-3*(-b+a-1/3*3A(1/2)*c)*b+3*a*