2010年全国大学生数学建模C题优秀论文

2024年8月28日发(作者：)

论文来源:无忧数模网

输油管的布置

摘要

“输油管的布置”数学建模的目的是设计最优化的路线，建立一条费用最省的输油

管线路，但是不同于普遍的最短路径问题，该题需要考虑多种情况，例如，城区和郊区

费用的不同，采用共用管线和非公用管线价格的不同等等。我们基于最短路径模型，对

于题目实际情况进行研究和分析，对三个问题都设计了合适的数学模型做出了相应的解

答和处理。

问题一：此问只需考虑两个加油站和铁路之间位置的关系，根据位置的不同设计相

应的模型，我们基于光的传播原理，设计了一种改进的最短路径模型，在不考虑共用管

线价格差异的情况下，只考虑如何设计最短的路线，因此只需一个未知变量便可以列出

最短路径函数；在考虑到共用管线价格差异的情况下，则需要建立2个未知变量，如果

带入已知常量，可以解出变量的值。

问题二：此问给出了两个加油站的具体位置，并且增加了城区和郊区的特殊情况，

我们进一步改进数学模型，将输油管路线横跨两个不同的区域考虑为光在两种不同介质

中传播的情况，输油管在城区和郊区的铺设将不会是直线方式，我们将其考虑为光在不

同介质中传播发生了折射。在郊区的路线依然可以采用问题一的改进最短路径模型，基

于该模型，我们只需设计2个变量就可以列出最低费用函数，利用Matlab和VC++ 都可

以解出最小值，并且我们经过多次验证和求解，将路径精度控制到米，费用精度控制到

元。

问题三：该问的解答方法和问题二类似，但是由于A管线、B管线、共用管线三者

的价格均不一样，我们利用问题二中设计的数学模型，以铁路为横坐标，城郊交汇为纵

坐标建立坐标轴，增加了一个变量，建立了最低费用函数，并且利用VC++解出了最低费

用和路径坐标。

关键字： 改进的最短路径 光的传播 Matlab 数学模型

1

输油管的布置

一、问题的重述

某油田计划在铁路线一侧建造两家炼油厂，同时在铁路线上增建一个车站，用来运

送成品油。由于这种模式具有一定的普遍性，油田设计院希望建立管线建设费用最省的

一般数学模型与方法。

利用模型分析管线布置和管线费用的情况，具体问题如下：

1. 针对两炼油厂到铁路线距离和两炼油厂间距离的各种不同情形，提出你的设计

方案。在方案设计时，若有共用管线，应考虑共用管线费用与非共用管线费用相同或不

同的情形。

2. 设计院目前需对一更为复杂的情形进行具体的设计。两炼油厂的具体位置由附

图所示，其中A厂位于郊区（图中的I区域），B厂位于城区（图中的II区域），两个

区域的分界线用图中的虚线表示。图中各字母表示的距离（单位：千米）分别为

a

= 5，

b

= 8，

c

= 15，

l

= 20。

若所有管线的铺设费用均为每千米7.2万元。 铺设在城区的管线还需增加拆迁和

工程补偿等附加费用，为对此项附加费用进行估计，聘请三家工程咨询公司（其中公司

一具有甲级资质，公司二和公司三具有乙级资质）进行了估算。估算结果如下表所示：

工程咨询公司

附加费用（万元/千米）

公司一

21

公司二

24

公司三

20

请为设计院给出管线布置方案及相应的费用。

3. 在该实际问题中，为进一步节省费用，可以根据炼油厂的生产能力，选用相适

应的油管。这时的管线铺设费用将分别降为输送A厂成品油的每千米5.6万元，输送B

厂成品油的每千米6.0万元，共用管线费用为每千米7.2万元，拆迁等附加费用同上。

请给出管线最佳布置方案及相应的费用。

2

二、模型假设

1、管道均以直线段铺设，不考虑地形影响。

2、不考虑管道的接头处费用。

3、不考虑施工之中的意外情况，所有工作均可顺利进行。

4、共用管线的价格如果和非公用管线不一致，则共用管线价格大于任意一条非公用管

线价格，小于两条非公用管线价格之和。

3

三、符号说明

h：共用管道的高度（问题一中b）

h1:共用管道高度

h2:管线与分界线的交点到B厂与铁路平行线的距离

w：方案的经费

a：A厂到铁路的距离

b：B厂到铁路的距离

c:A厂到城郊分界线的距离

l:A、B两厂之间的铁路长度

x：A厂离共用管道的距离（问题一中的c）

y：共用管道的高度（问题一中的c）

m：共用管道的费用（问题一）

n：非共用管道费用（问题一）

y1：为o点的纵坐标

y2：为o1点的纵坐标

x1：为o点的横坐标

x2：为o1点的横坐标

L: 为管线总长度（问题一中的b）

4

四、问题分析

问题一：要考虑有和没有共用管线，还要考虑共用管线与非共用管线费用相同和不

同两种情况。同时还要考虑两个工厂是否在铁路的同一侧，如果两个工厂在铁路的同一

侧那么一定没有共用管线。 不在铁路的同一侧那么就要考虑有和没有共用管线这个问

题。计算共用管线的长度时，用光学原理，把一个工厂当作光源发射一束光经过一个平

面的反射通过另一个工厂，这样能够保证路线最短。这个平面与铁路的距离即为共用管

线的长度。同时与这个平面的交点就是两厂的管线的交点。当共用管线与非共用管线费

用不相同时可以通过建立方程组来解答。

当共用管线与非共用管线费用不相同时要建立方程组来计算其最小费用从而来确定

方案的可行性，共用管线与非共用管线长度作为变量来控制总费用，那么我们就可以列

出一个方程组，从而在变量的约束条件下可以确定最小费用。

问题二：把这个问题分两部分来考虑，即市区和郊区分两个部分，火车站建立在郊

区费用要小得多，郊区共用管线与非共用管线的费用相同所以可以用最短路径的方法来

考虑，同时又要求费用最小，可以解出最低费用及对应的铺设线路。

问题三：通过建立坐标系设两个点的坐标，同时也是表达管线的长度，然后再与各

自的费用之积确定总的费用，从而算出两点的坐标值。即确定了管线的路线。

5

五、模型的建立与求解

5.1关于问题1的模型建立与求解

对于管线布置的分析，分为两种情况：

1. 两厂分别在铁路的两侧如下图：

那么连接两厂A、B与铁路的交点C即为火车站的位置。

2. 当两厂位于铁路的同一侧时,此时要分有公用管线与没有公用管线两种情况。

a.当没有公用管线时，此时找出两厂与铁路交点连线的最近路线即可，如图：

过铁路作A厂的对称点A’，连接A’B与铁路交于一点C，该点C即为火车站的位置。

b.当有共用管线时又要分为共线管线费用与非共线管线费用相同与不同两种情况：当共

线管线与非共线管线相同时，费用为m万元/千米如图所示：

6

假设共线管线的长度为h，A厂到铁路的距离为a，B厂到铁路的距离为b，则总的管线

长度为：

L[(ah)(bh)]

2

l

2

h

(0hb)

则总费用：

W1Lm

c.当共线管线与非共线管线不同时，共用管线费用为m万元/千米 ，非共用管线费用为

n万元/千米，如图所示：

总费用为：

W(my(ay)

2

x

2

(lx)

2

(by)

2

)n

其中

0xl

0ymax(a,b)

实际的费用可以根据已知道的常量a、b、l再结合x、y的取值范围可以得出最小费

用。

5.2关于问题2的模型建立与求解

因为在城区和郊区铁路管线的费用相同，但城区要增加拆迁和工程补偿等费用，因

此城区和郊区要分为两部分来考虑。我们考虑三家咨询公司给出的三个方案，我们考虑

到甲级资质和乙级资质的评估准确性，首先排除掉公司二的预算，对于公司一和公司三

的预算，我们将分别求出最小费用，考察两者的差别。

1．假设共用管线在郊区把该模型看作是一束光从B点发射在分界处G点发生了折射，

把左边的问题看作是最短路径问题，如图所示：

7

设共用管线的长度为h1，G点到O2B的距离为h2。在区域Ⅱ中即BG段每千米的费

用为：20+7.2=27.2万元。

由以上分析数据可得如下关系式：

总费用： W1（

最小

）=

((8h

1

h

2

5h

1

)

2

15

2

h

1

)7.227.2h

2

25

（式1）

参数

h

1

的取值范围：

0h

1

8

（式2）

参数

h

2

的取值范围：

0h

2

8

（式3）

利用Matlab将式（1）（2）（3）联立关系式绘图：

2

8

用Microsoft Visual C++ 6.0解:

W1(最小)= 275.13404万元

运行结果：

在这种情况下采用公司一的预算，只需要在上式中将27.2增加为28.2即可，计算得

到总费用：280.177831万元

运行结果：

2．假设共用管线在城区同理，如图所示：

由以上分析数据可得如下关系式：

9

总费用：

W2（最小）=

27.2(h

2

[8h

1

2h

2

]

2

25)5.615

2

(5h

1

)

2

（式1）参

数

h

1

的取值范围：

0h

1

8

（式2）

参数

h

2

的取值范围：

0h

2

8

（式3）

用Microsoft Visual C++ 6.0解得

W

2（

最小

）

= 355.25587

运行结果：

显然W1（

最小

）

最小

）方案一费用少于方案二，因此舍掉这种方案。

最终求得的结果为，如果采用一咨询公司的估算价格，则最终费用为275.134304

万元，如果采用三咨询公司的估算价格，则最终费用为280.177831万元，考虑到公司

一具有高级资质，因此我们采用公司一的价格方案，将最终预算设为280.177831万元，

但是实际铺设管道的价格有可能在两种估算价格之间。

5.3关于问题3的模型建立与求解

1、O点为B管线与分界线的交点，O1点为A管与B管的交点，如下图建立坐标轴，采

用公司三的估算费用，

总费用等于各段路线的长度与各段费用的积为：

W5.6AO16.0OO1(620.0)OB7.2O1O2

坐标法解答，A01，OO1，OB，如图：O(x1，y1),O1(x2，y2)

10

由以上分析数据可得如下关系式：

A厂到管道交接点O1的长度：

AO

1

=

(x

2

15)(y

2

5)

2

（式1）

管道交点O1到B厂与城郊分界线交点O的长度：

OO

1=

x

2

(y

1

y

2

)

2

（式2）

B厂到交点O的长度：

OB=

25(8y

1

)

2

（式3）

铁路站点O2到交叉管道O1的长度：

O

1

O

2=

y

2

参数

x

2

的取值范围：

15x

2

0

（式5）

（式

2

4）

参数

y

1

的取值范围：

0y

1

8

（式6）

参数

y

2

的取值范围：

0y

2

8

（式7）

总费用：

W5.6(x

2

15)(y

2

5)

2

6.0x

2

(y

1

y

2

)

2

26.025(8y

1

)

2

7.2y

2

11

2

由以上式子利用Microsoft Visual C++ 6.0软件求得最小经费：

W3（最小值）= 244.386494万元。

运行结果：

在同种情况下，用公司一的预算费的总费用：

运行结果：

当火车站建在市区费用太高同二题中的方案二故不选用那种方案。因此类似于问题

二，我们采取公司一的估算价格，最终预算为249.468791万元。但是公司三的价格也

具有一定参考性，实际铺设管道价格应为244.386494万元到249.468791万元之间。

12

六、模型的评价与应用

从实际的生活出发输油管道是石油生产过程中的重要环节，是石油工业的动脉。在

石油的生产过程中，至始至终都离不开输油管道。我们可以把石油的生产过程简单的表

示为：

油→ 计量站→ 井联合站→ 转油站→ 矿场油库→ 炼油厂→ 用户

从油井出来的油气通过管道输送到计量站，经过计量后又由管道输送往联合站，在

联合站生产出合格的原油，合格原油通过管道和转油站输到矿场油库或外输到管道首

站，通过长输原油管道输到炼油厂加工精练，生产出各种产品，通过成品油管道或铁路、

公路、水路将各种产品送往用户，其中成品油管道就需要用到管道的布置设计。

优点：

模型使问题由复杂变简单，方便运输，提高输油效率，规划线路。管线布置和规划

及相应的费用减到最小，在不同的环境下用这种环境中的最优模型，方便快捷，节约开

支，使实际问题更加精确。同时对于题目的三个问题都设计了合适的模型，并且当给出

具体数值的时候能够给出足够精确的解，具有一定的普遍性。

缺点：

该模型在提出的时候将部分因素没有考虑进来，例如管线接头处的费用，以及工作

工程中的一些意外情况等等，使得该模型在实际应用中会缺少精确性。

应用：

模型在实际运用中，不仅仅可以用在成品油运输管布置，还可运用到原油输送和污

水处理，电线电缆的布置还有公路铁路的修建等一些列的线路布置问题。

13

七、参考文献

【1】 赵静 但琦 《数学建模与数学实验》第三版 22-29页，178-194页 高等教

育出版社 2008年1月

【2】 曹戈 《MATLAB教程及实训》 37-60页 机械工业出版社 2008年5月

【3】 邬学军 周凯《数学建模竞赛辅导教程》 73-96页 浙江大学出版社 2009年1月

14

附录：

问题2程序

1、按照公司三的评估总费用为：

#include

#include

void main()

{

double h1,h2,w;

double a,b;

double min = 10000;

for(h1=0;h1<=8;h1+=0.001)

for(h2=0;h2<=8;h2+=0.001)

{

if(h1+h2>8)

continue;

w=27.2*sqrt(25+h2*h2)+(sqrt((5-h1+8-h1-h2)*(5-h1+8-h1-h2)+225)+h1)*7.2;

if(min>w)

{

min=w;

a=h1;

b=h2;

}

}

printf("%f n",min);

printf("%f %f n",a,b);

}

15

运行结果：

2、按照公司一评估总费用为：

问题3程序：

1、按照公司三评估总费用为：

#include"stdio.h"

#include"math.h"

main()

{ double x2,y2,y1,w,a,b,c,min=1000;

for(x2=-15;x2<=0;x2+=1)

for(y1=0;y1<=8;y1+=1)

for(y2=0;y2<=8;y2+=1)

{

w=5.6*sqrt((x2+15)*(x2+15)+(y2-5)*(y2-5))+6.0*sqrt((x2*x2)+(y1-y2)

*(y1-y2))+26.0*sqrt(25+(8-y1)*(8-y1))+7.2*y2;

if(min>w)

{

min=w;

a=x2;

b=y2;

c=y1;

}

}

printf("min=%f nx2=%f ny2=%f ny1=%fn",min,a,b,c);

}

16

运行结果：

将循环因子的步长降低之后，进一步精确求解值为：

#include"stdio.h"

#include"math.h"

main()

{ double x2,y2,y1,w,a,b,c,min=1000;

for(x2=-9;x2<=-7;x2+=0.001)

for(y1=6;y1<=8;y1+=0.001)

for(y2=0;y2<=1;y2+=0.001)

{

w=5.6*sqrt((x2+15)*(x2+15)+(y2-5)*(y2-5))+6.0*sqrt((x2*x2)+(y1-y2)

*(y1-y2))+26.0*sqrt(25+(8-y1)*(8-y1))+7.2*y2;

if(min>w)

{

min=w;

a=x2;

b=y2;

c=y1;

}

}

printf("min=%f nx2=%f ny2=%f ny1=%fn",min,a,b,c);

}

2、按照公司一的评估的总费用为：

17

18