管道线路布置的优化设计(DOC)

管道线路布置的优化设计

摘要

管道运输是输送石油的一个重要途径，设计合理的管线铺设方案，不仅可以节省铺设的费用，还可以减少后期运输的成本，提高经济效益。本文针对题目中给出的不同情况，设计了不同情况输油管线的详细方案。

针对问题一：根据两个炼油厂到铁路线距离和两个炼油厂间的不同距离以及共用管线与非共用管线的两种不同情况，对不共用管线时进行B A B A w w w w ≠=,的分析，对共用管线时进行S B A S B A w w w w w w ≠===,，S B A w w w ≠≠的分析。最终可以将模型归纳为：运用轴对称定理建立的非线性优化模型。在模型检验中运用费马点对模型进行检验，可以证明该模型的正确性。

针对问题二：在已经确定了两个炼油厂的地点一个在郊区一个在城区的情况下，由于在城区的管道铺设还需增加拆迁和工程补偿等附加费，首先按照各级公司的各项数据运用matlab 进行Topsis 综合评价法分析，得到甲，乙等级公司的评估可信度之比为1:0.426，从而得到拆迁和工程补偿等附加费用的期望值为54.23=P W 万元/千米，。然后根据问题一中的模型三，运用Lingo 编程的方法，得出将火车站建在点（4.427985,0) 时费用最少为502.6264万元，此时的城郊结合处坐标为(15，6.547257)，无共用管线，两厂管道交汇处坐标为(4.427985，0.4435016)。在用Lingo 求解得到费用最小的线路后，控制变量x ，保持y 和1y 的条件不变，对x 进行灵敏度分析，可以总结出如下结论：当x 的值大于4.427985时，随着x 值得增大，y 和1y 的值都在小幅度的减小，以此来保证费用较小。

针对问题三：根据题中给出的数据，可以将火车站分为建立在城区和郊区两种情况，根据通用模型三，运用Lingo 编程的方法，将已知数据代入，得到将火车站建在郊区坐标为（5.323864,0）时费用最少为458.6181万元。为了检验计算的准确性，利matlab 编程进行模拟，得到最小总费用为523.6968万元，火车站的坐标点位（14.9820,0），共用管道的坐标为（14.9820,0.0772）.由此可得，我们建立的模型是可行的。

针对论文的实际情况，对论文的优缺点做了评价，文章最后还给出了其他的方法，以用于参考。

关键词：轴对称定理 非线性优化模型 费马点 T o p s i s 综合评价法

1、问题重述

某油田计划在铁路线一侧建造两家炼油厂，同时在铁路线上增建一个车站，用来运送成品油，成品油由输油管线运往火车站。

问题一，要针对两炼油厂到铁路线距离和两炼油厂间距离的各种不同情形提出设计方案，并要考虑是否使用共用管线以及共用管线费用与非共用管线费用相同或不同的情形。

问题二，对于一个具体实例，在所有管线的铺设费用均为每千米9.5万元时，但铺设在城区的那部分管线，还需增加拆迁和工程补偿等附加费用。为对此项附加费用进行估计，聘请了三家工程咨询公司进行了估算。三家工程咨询公司的资质分别为甲级、乙级、乙级，附加费用的估价分别为每千米24万元、21万元、25万元。要求给出管线布置方案及相应的费用。

问题三，在与问题二相同的实例中，为进一步节省费用，可以根据炼油厂的生产能力，选用相适应的油管，这时的管线铺设费用将分别降为输送A 厂成品油的每千米7.6万元，输送B 厂成品油的每千米8.0万元，共用管线费用为每千米11.2万元，拆迁等附加费用同上。要求给出管线最佳布置方案及相应的费用。

2、问题分析

本题中，某油田计划在铁路一侧建造两家炼油厂，同时在铁路上增建一个车站，用于运送成品油。

针对问题一，根据两个炼油厂到铁路线距离和两个炼油厂间的不同距离以及共用管线与非共用管线的两种不同情况，在不共用管线时设一个炼油厂为A,一个炼油厂为B ，从炼油厂A 到火车站的管线费用为A w 万元/千米，从炼油厂B 到火车站的管线费用为B w 万元/千米，对不共用管线时进行B A B A w w w w ≠=,的分析；在共用管线时，设共用管线的费用为S w 万元/千米，对共用管线时进行,,,S B A S B A S B A w w w w w w w w w ≠≠===的分析。最终可以将模型归纳为：运用轴对称定理建立的非线性优化模型。在模型检验中运用费马点对模型进行检验，可以证明该模型的正确性。

针对问题二，在给定数据的条件下铺设管线时，因为所有管线均为9.5万元/千米，因此不考虑共用和非共用管线的价格不同的情况，但因为在城区的管线需增加拆迁和工程补偿等附加费用，设计学院聘请了三家不同资质的咨询公司，得到3个估价，因此需要对着3个估价进行数据处理，得到一个加权后的附加费用估计值。确定在郊区和在城区的管道长度，就能得到费用的最低值。

针对问题三，对于问题二这个实例，为了进一步节省费用，选用相适应的的油管，这就存在有共用管道和没有共用管道的两种情况，对模型进行比较，就可以得到最优解，最后利用Matlab 编写模拟程序进行模拟，得出模拟值与其比较，检验计算的准确性。

3、模型假设与符号说明

3.1模型假设

假设一：城区和郊区地形良好，管道在城区与郊区都能直线铺设；

假设二：在铺设管道过程中，不考虑由于河流、山坡等障碍而增加费用；

假设三：共用管道与非共用管道接口处的长度忽略不计；

假设四：管道铺设在边界线上时不算入拆迁和工程补偿等附加费用；

假设五：不考虑由于在铺设管道时造成的意外事故所赔偿的费用；

假设六：管道铺设后不会对周围的环境造成污染；

3.2符号说明

符号解释

A炼油厂A的地点

B炼油厂B的地点A

C铁路上的火车站点

C 共用管线的起点

b炼油厂A点的纵坐标

1

a炼油厂B的横坐标

2

b炼油厂B的纵坐标

2

x铁路上的火车站点横坐标

y铁路上的火车站点纵坐标

Z管线建设的总费用

w通往炼油厂A的管线每千米的费用

A

w通往炼油厂B的管线每千米的费用

B

w共用管线每千米的费用

S

w最终估计的附加费用

p

p公司可信度估计权重

4、模型的准备

在已经确定了两个炼油厂的地点一个在郊区一个在城区的情况下，由于在城区的管道铺设还需增加拆迁和工程补偿等附加费用，对此项附加费用进行估计。

在得到三家工程咨询公司的估算价格之后，由于三家公司的资质情况和估算费用结果不经相同，如下表所示：

因公司一具有甲级资质，公司二和公司三具有乙级资质，我们在查阅中国工程咨询公司资格认定[6]方法后，按照各级公司的各项数据运用Matlab 进行Topsis 法分析，得到二者的与最优方案接近程度比为1:0.426为我们将附加费用按此照权重进行再次估计：

332211E P E P E P W P ++= 根据公司不同赋予不同权重：

%;23

%;23%,54321===P P P 计算可得拆迁和工程补偿等附加费用的期望值为54.23=P W （万元/千米）。

5、模型的建立于求解

５.1问题一：

5.1.1.1模型一的建立： 不共用管线的情况

首先考虑没有共用管线的方案，管线铺设总费用主要包括两部分： 首先设点A 的坐标为),0(1b A ，点B 的坐标为),(22b a B ，A 厂到车站的管线铺设费用A w 万元每千米、B 厂到车站管线铺设费用B w 。在没有共用管线情况下，我们应该将车站建铁路线上，不妨令该点的坐标为),(y x C 。因为不共线时C 点在轴上面，即C 点的坐标为)0,(x C ，所以我们有2222212)(_,b x a BC b x AC +-=+=，从而总建设费用为：

图（5.1.1）

工程咨询公司 公司一 公司二 公司三 附加费用（万元/千米） 24 21 25

O

)

,(y x C

)

,(22b a B

)

,0(1b A

2

222212)(b x a w b x w Z B A +-++=

那我们的问题转化为求点C 的位置（求x 的值），使得管线建设费用最少。为此我们可以得到模型一：

2

2222

12)(b x a w b x w MinZ B A +-++=

(5.1.1)

??????

?>≤≤≥>0,0,0,0,..2221B A w w a x a b a t s 对于模型一：

2

2

222

12)(b x a w b x w MinZ B A +-++= (5.1.2)

???????>≤≤≥>0,0,0,0,..22

21B A w w a x a b a t s

下面确定车站位置使得总建设费用最小。

利用费尔马极值原理，考虑其导数

2

2

2

222

1

2

)()(b x a x a w b x x w dx

dz

B

A +---+= (5.1.3)

讨论：

Step1、当B A w w =，即两炼油厂的每千米的管线铺设费用相等，则(5.2)式可转化为：

2

2

2

222

1

2

)()(b x a x a w b x x w dx

dz

A

A +---+=

令0=dx dz

，即：

0)()(22

222212=+---+b x a x a b x x

(5.1.4)

当21b b =时，如图5.1.2所示。

由于每千米的管线建设费用相同，且不考虑拆迁和工程补偿等附加费用，只需建设车站使得A 厂和B 厂到车站的管线总长度最小即可。由光的反射定理，可知将火车站建在直线'AA 与x 轴的交点处（其中'A 点是A 点关于x 轴的对称点）。

解得：

???

?

??

?

-=

+=

)()(21222121b a b a x b a b a x 讨论这两个解：

在没有共用管线的情况下，且炼油厂A 与炼油厂B 的每千米管线建设费相等时，即B A w w =。

图(5.1.3)

如图5.1.3所示：由于炼油厂A 与B 的管道建设费用相等，即A w =B w ，则连接车站到炼油厂A 的距离CA 与车站到炼油厂B 的距离CB 之和最短，为最优。在平面直角坐标系中，以x 轴为对称轴，找一点),0('a A -使得点'A 与点A 关于x 轴对称。连接点B A '，

图(5.1.2)

O

)0,(x C

),(22b a B

),0(1b A ),0(1'b A -

O

)0,(x C

),(22b a B

),0(1b A ),0(1'b A -

与x 轴相交于点C 。故有：+=CA BA '

CB ,根据两点间直线距离最短可得出点C 到点

A 的距离加上点C 到点

B 的距离为最短。故点

C 必须在X 轴上。

根据两炼油厂都在铁路一侧，有0>a 且0>b ，首先考虑第一个解：

22121)

(0a b a b

a x ≤+=

≤，故1x 有意义。

将1x 代入(5.1.1)式，则最省的费用为：

)(212b b a w Z A ++=

其中车站应建在点)0,(1x C 处。 其次考虑第二个解2x ，

特别的，当21b b =时，点C 的坐标为)0,2

(

2

a ； 当2

1b b ≠时，

(a)、如果21b b >，则22122)(a b a b

a x >-=

，

(b)、如果21b b <，则0)

(2122<-=

b a b

a x 。

故解2x 无意义。 Step2、当B A w w ≠，则

2222()A B Z w x a w l x b =++-+

其导数为：

2

2

2

222

1

2

)()(b x a x a w b x x w dx

dz

A

A +---+=

(5.1.5)

令0,dz

dx

=可得x 的符号解。

5.1.2.1模型二的建立： 共用管线的情况

a.当两炼油厂的连线所在直线垂直于铁路线时，如图5.1.4所示（A ，B 上下位置可互换）：

A

(5．1.4)

B

O 铁路线

明显的，从远油厂铺设非共用管道到近油厂，再从近油厂铺设共用管道到增建的车站，费用最少（注：通过相关资料查得安全距离大致应在1000米以上,所以 AB 的距离

m AB 1000≥ ）。

此时，总费用为： OB w AB w Z S A ?+?= 当S A w w =时：

)(OB AB w Z A +?=

当S A w w ≠时： OB w AB w Z S A ?+?=

b.当两炼油厂的连线所在直线不垂直于铁路线时，建立模型二，如图5.1.6所示（其中21,b b 大小不确定）： 由于两种管道的单位造价相同，要使总铺设费用最小，则应使铺设距离最短，为求得最小距离，则必须确定两种管道的结合点O 。 于是问题转化为在四边形ABDO 内找一点'C ，使'''CC BC AC ++最小。据实际情况判断，'C 应界定在

矩形AEDO 内。

图（5.1.5）

)

0,(x C

)

,('y x C

)

0,(2a D

d

)

,(22b a B O

)

,0(1b A

E

铁路线

如图5.1.5所示，设车站建在点)0,(x C 处，由于两种管道的单位造价相同，要使总铺设费用最小，则应使铺设距离最短，为求得最小距离，则必须确定两种管道的结合点

'C ,连接尺寸CC ’垂直于铁路线，炼油厂A 与炼油厂B 的管线相交于点),(y x C 处，使

'''CC BC AC ++最小。

炼油厂A 的管线建设费用为A w ，炼油厂B 的管线建设费用为B w ，共用管线建设费用为S w 。则总建设费用为：

y w b y x a w b y x w Z s B A +-+-+-+=2222212)()()(

在此情况下，我们的问题实际上可以转化为模型二：

??

?≤≤≤≤+-+-+-+=)

,min(0,0.)()()(2122222212b b y a x t s y w b y x a w b y x w MinZ s B A

对于模型二：

??

?≤≤≤≤-+-+-+=)

,min(0,0.)()()(2122222212b b y a x t s b y x a w b y x w MinZ B A

讨论:

Step1、当s B A w w w ==时，即所有管道价格一样

则y w b y x a w b y x w Z A A A +-+-+-+=2222212)()()(

其对x 的偏导数：

???

?

?

?-+--+

-+=??22222222)()()()(b y a x a x b y x x w x z A 对y 的偏导数

???

? ??+-+--+-+-=??1)()()()()(222222121b y a x b y b y x b y w y z A

令

x z

??=0，y

z ??=0 有方程组：

????

?????=???

? ??+-+--+-+-=????

?

?-+--+

-+01)()()()()(0)()()()(2

2222212

122222222b y a x b y b y x b y w b y a x a x b y x x w A A 解方程组得：

???

???

?-+-=-+=)(2163)(2

3221212121b b a y b b a x ???

????++=-+=)(2163)(232212

21222b b a y b b a x

把其结果代入目标函数有：

()

(

)

2

2122139156

39156a

b b

c Z a b b c

Z +-=--=

Step2、当S B A w w w ≠=时.即共管线，且共用管线费用与非共用管线费用不同， 其目标函数为：

y w b y a x w b y x w Z S A A +-+-+-+=2222212)()()(

其对x 的偏导数：

???

?

??-+--+

-+=??22222212)()()()(b y a x a x b y x x w x z A 对y 的偏导数

s

A w b y a x b y b y x b y w y z

+???

?

??-+--+

-+-=??222222121)()()()()(

令

x z

??=0，y

z ??=0 有方程组：

?????????=+???

?

??-+--+-+-=????

?

?-+--+

-+0

)()()()()(0)()()()(2

2222212

122222212s

A A w b y a x b y b y x b y w b y a x a x b y x x w 利用matlab 解方程组得：

????

????

?

?

?? ??-+----

=??? ??-+---+=2222221221

2

2222212212142)82)(()4(4142)82)((2)(21S A S A S A S S A S A S w w w a w w b b w w a y w w w a w w b b a b b a x

????????

???? ??------

=??

? ??-----+=2222221222

22222212212242)82)(()4(4142)82)((2)(21S A S A S A S S A S A S w w w a w w b b w w a y w w w a w w b b a b b a x

Step3、当B A S w w w ≠≠，即所有费用都不相同的情况下，管道费用为：

y w b y a x w b y x w Z S B A +-+-+-+=2222212)()()(

对它求x 的偏导数：

2

22

222

12

)

()()()

(b y a x a x w b y x x w x

z B A -+--+

-+=??

求对y 的偏导数：

2

2

2

2

()()()

()()

A B w y a w y b z w

y

x y a x l y b --?=+

+?+--+-

s B A w b y a x a x w b y x x w x

z +-+--+

-+=??2

22

222

12

)

()()()

(

令0z x ?=?，0z y ?=? 有方程组：00z

x

z y

??=??????=???，对于给定的数值，可以求出(),x y ,

然后可以求出目标的最优值Z 。

5.1.3模型三的建立

综合以上的模型一与模型二可以建立一个综合模型三：

??

?≤≤≤≤+-+-+-+=)

,min(0,0.)()()(2122222212b b y a x t s y w b y x a w b y x w MinZ s B A

当0=y 时：2

2222

12)(b x a w b x w MinZ B A +-++=

当0≠y 时：y w b y x a w b y x w MinZ s B A +-+-+-+=2

222212)()()(

根据上述思路编写matlab 程序，流程图如下所示。

Start

初始化，给坐标及费用赋值

是否溢出？利用轴对称求最短距离原理，计算当前费用cost_cur

冒号法生成序列

在虚线（即城郊分界线）上，以0.01为步长，按y 坐标从0到8生成一个序列。每个y 值对应虚线上的一个坐标点。

是否比cost_min 小记录当前E 点的y 坐标，更

新cost_min 的值

Yes

Yes

NO

NO

END

注：后面的各程序不同之处在于此，即通过两点求会合点的算法不同

5.２问题二：

5.2.1模型四的建立

在此种情况下两炼油厂位置分别处于郊区和城区，所以我们分别讨论车站建在郊区和车站建在城区的情况，分别计算出两种方案管线所需要的费用。最终确定方案的费用Z 为车站位于郊区的费用1Z 和车站位于城区的费用2Z 二者的最小值。按照问题1所简历的模型，我们以铁路为x 轴，A 点与铁路的垂线为y 轴做出如下图：

图（5.21）

当火车站位于郊区时（如上图a ），位于城区的管道线路就只有BD 段，：

()P B B S A W W BD W DC W CC W AC Z +?+?+?+?='''1

当火车站位于城区时（如上图b ），位于城区的管道有'

BC 、'

CC 和'

DC 三段，

()()()P B P S P A A W W BC W W CC W W DC W AD Z +?++?++?+?='''2 所以对于模型四为

()()()()???+?++?++?+?=+?+?+?+?==P B P S P A A P B C B A W W BC W W CC W W DC W AD Z W W BD W DC W CC W AC Z MinZ ''''''2

1

5.2.2模型四求解：

由题目可求得4.23=p (万元/千米)。根据题意，)(3千米=a ，8()b =千米，

15()c =千米，20()l =千米，)/(5.9千米万元===S B A w w w 。带入题目可得 当火车站建在郊区时：

()()()()())

54.235.9(815255.9153min 2

122212

21+?-+-+???

? ?

?-+-+

+-+=y x y y y y x Z ??

???≤≤≥≤≤8

0015

0..1y y x t s 当火车站建在城区时：

()()()()()()

54.235.9825155.9315min 222

122

122+???

?

?

?-+-+

+-+-+?-+=y x y y y x y Z 郊区

)

,(b l B

)

,0(1a A 城区

郊区

城区

)

,0(1a A )

0,(x C )

,('y x C

)

,(1y c D

)

,(b l B

)

,(1y c D )

0,(x C )

,('y x C

O

a

O

b

L

L

??

???≤≤≥≤≤8

002515..1y y x t s 运用Lingo 软件编程求解（程序见附录）

可以得到以下两种情况下的结果如下表：

车站建在郊区 车站建在城市 总费用（万元） 502.6264 568.4405 城郊结合处坐标（x,y ） (15，6.547257) (15,0) 铁路坐标(x,y) (4.427985,0) (15,0) 两厂管道交汇处坐标（x,y ） (4.427985，0.4435016) 无共用管线

则有 502.6264)568.4405, 502.6264

min()min(2,1===Z Z Z 所以最优方案总花费为502.6264万元，车站建在距离坐标远点4.427985千米的火车

道上。，两油厂的管线在坐标(4.427985，0.4435016)处交接，然后使用公共管线垂直到达铁路线上的车站。

用matlab 做出其准确图形如下：

5

10

1520

25

30

-1

01234

56789

10炼油厂A 炼油厂B

管线的城郊交界处

两炼油厂的管线交界处车站位置

O 千米

千米

5.3问题三

5.3.1 模型五的建立

在问题二中，我们已经建立了近似B A s

w w w ≠≠情况下的一般模型。本问题和问

题二相比，正好是管线铺设费用的数值不一样了的情况，根据问题2模型具有的推广性质，我们完全可以用问题二的模型来解决问题三。

()()()()???

+?++?++?+?=+?+?+?+?==P B P S P A A P B C B A W W BC W W CC W W DC W AD Z W W BD W DC W CC W AC Z MinZ ''''''2

1

5.3.2模型五的求解

根据题意，我们将数据)(3千米=a ，8()b =千米，15()c =千米，20()l =千米，

)/(6.7千米万=A W ，()千米万/8=B W ，()千米万/2.11=S W ，千米）（万/54.23=P W 代入模型得

将火车站建立在郊区时：

()()

)

54.238(815252.118

)15()(6.7)3(2

12

21221+?-+-+

?+?-+-+?-+=y y x y y y x MinZ

??

?

??≤≤≥≤≤8

00150..1y y x t s 将火车站建立在城区时：

()())

54.232.11()54.238()8()25(54.236.7)()15(6.73152

22122

122+?++?-+-+

+?-+-+?-+=y y x y y x y MinZ

??

?

??≤≤≥≤≤8

002515..1y y x t s

用Lingo 编程求解（程序见附件）可解出结果如下：

车站建在郊区 车站建在城市 总费用（万元）

458.6181

520.1667

管线的城郊处坐标（x,y ） （15，6.561749） （15，0.2040007E-07） 铁路坐标(x,y)

（5.323864,0） （15,0） 两厂管道交汇处坐标（x,y ）

无共用管线

无共用管线

所以最优方案总花费为458.6181万元，车站建在距离坐标远点5.323864千米的火车道上。，两厂无共用管道。

用matlab 做出其准确图形如下：

5

10

1520

25

30

-1

01

234

56

789

10炼油厂A 炼油厂B

管线的城郊交界处

车站位置

O 千米

千米

6、模型结果分析与检验

对于问题一，当管线费用相同我们出了用以上方法证明之外，我们也可以根据费马点来来证明此点为管线交接的最优点有且仅有此点三角形的内角都小于120°的情况：首先证明AA',BB',CC'三条线交于一点。 设P 为线段CC'和BB'的交点。注意到三角形AC C'和三角形BAB'是全等的，三角形AC C'可以看做是

三角形AB B'以A 点为轴心顺时针旋转60度得到的，所以角∠PB C'等于60度，和∠AB C '相等。因此，

P C B A '...四点共圆。同样地，可以证明P C B A .'..四点共圆。于是：

∠=APB ∠0120=APC 从而∠0

120=CPB 。于是可以得出：P C B A .'.'.四点共圆，即

∠

=PB A '∠060'=CB A ,∠=PA A '∠+APB ∠0060120+=APB

P A A '..三点共线。也就是说''.'.AA BB CC 三条线交于一点。接下来证明交点P 就是到三个顶点距离之和最小的点。

在线段'AA 上选择一点Q ，使得PC QP =。由于∠060=QPC ，所以等腰三角形PQC 是正三角形。于是∠=PCB ∠'QCA 。同时PC QC =、C A' BC =，于是可以得出三角形

BPC 和三角形QC A'是全等三角形。所以PB QA'=。综上可得出： AA' PC PB PA =++对

于平面上另外一个点P'，以C P'为底边，向下作正三角形C Q'P'。运用类似以上的推理可以证明三角形C BP'和三角形C Q'A'是全等三角形。因此也有：

A Q' + Q'P' + AP' = C P' +

B P' +A P'平面上两点之间以直线长度最短。因此

PC. PB PA AA' Y A'Q' Q'P' AP' C P' B P' A P'++=++=++ 也就是说，点P 是平面上到点

C B A ..距离的和最短的一点。最后证明唯一性。

如果有另外一点P'使得PC PB PA C P' B P' A P'++=++，那么A Q' Q'P' AP' AA'++=因此点P'和Q''也在线段AA'之上。依照P'和Q'的定义，可以推出

∠

=B AP '∠0120'=C AP

因此P'也是AA'BB'CC'三条线的交

点。因此P'点也就是P 点。因此点P 是唯一的。一内角大于120°的情况。 如右图，∠BAC 大于120°，P 为三角形内一点。以BA 为底边，向上作正三角形BAF ；以PA 为底边，向上作正三角形PAQ 。于是三角形AQF 和三角形APB 是全等三角形。

PB FQ =。所以PC QP FQ PC PB PA ++=++. 延长FA 交QC 于D 点，则

QC FQ PC QP FQ +>++DC FD DC QD FQ +>++=AC FA DC AD FA +>++=AC AB +=. 即AC AB PC PB PA +>++.

所以A 点到三顶点的距离比三角形内任意一点到三顶点的距离都小，即A 点为费马点。

对于问题2的模型，我们使用的是赋权值进行机选取得加权平均数的方法来确定估计值，但是在实际的工程中，人们常使用其它经济学方法来进行分析讨论，例如乐观决策法（好中求好，附近费用为21万元/千米），悲观决策法（小中取大,附加费用为25万元/千米）等方法来确定附近费用的值，但由于存在较大的误差。所以我们选用了加权

平均数来确定附加费用。有较强的说服力。在用Lingo 求解得到费用最小的线路后，控制变量x ，保持Y 和1Y 的条件不变，对x 进行灵敏度分析，可以总结出如下结论：当x 的值超过我们求出来的费用最小的x 值，即4.427985时，随着x 值得增大，Y 和1Y 的值都在小幅度的减小，以此来保证费用较小。 对于问题三的模型，虽然对于假设车站在城区的情况下计算出之后两厂有2CM 的共用管道距离，但是在现实生活中，考虑到链接管道时需要的额外费用，此距离是可以忽略不计的，所以我们可以认为在问题三的第二种情况下是没有共用管道的，所以模型仍有其的合理性。

7、模型推广与改进

在此模型中，我们的假设条件在现实使用中有一定困难，在实际运用上的提高还有待提高。同时模型可以应用到在铁路线一侧建立多个炼油厂的设计中，若只在铁路线上建设一个车站，则依然设共用管线与非共用管线的交点，若建立多个车站就不适用了，

应建立更好的模型。另外，考虑炼油厂建设在铁路线的两侧，则不需要考虑共用管线的问题，即只需假设出车站的位置即可。

8、模型的优缺点

模型的优点：模型的可操作性强，运用我们提供的模型，在直接修改数据的情况下就可以计算出全局最优解，具有一定的普遍性；运用的模型和理论都很浅显易懂，且具有较强的说服力。模型的缺点：为了建模的方便，对一些非主要的因素做了合理的假设，这就造成结果上的一定误差；运用得方法较单一，可比性不强。

参考文献

[1]谢金星、薛毅编著；优化建模与LINDO/LINGO软件,北京：清华大学出版社,2005年7月。

[2]姜启源、谢金星、叶俊；数学建模(第三版),北京：高等教育出版社,2004年。

[3]赵静、但琦；数学建模与数学实验,高等教育出版社施普林格出版社,2000年。

[4]姜启源等编著；大学数学实验,北京：清华大学出版社,2005年2月。

[5] 皮耶·德·费马；费马点，

https://www.360docs.net/doc/052813531.html,/link?url=hHHIAEWyqrCe0oLEGVQyi-qNek7y3DcT_BMKfgC CGNiILdB_gl1mZAVJIVvGGJJj

[6]《中华人民共和国国家发展和改革委员会令第29号—工程咨询单位资格认定办法》，https://www.360docs.net/doc/052813531.html,/zcfb/zcfbl/zcfbl2005/t20050602_6090.htm，2005年

附录

问题二模型及结果：

TOPISIS

a=[500 60 18 9 6 5 2 10

200 30 9 5 3 5 2 8

50 15 3 2 2 5 1 5];

[m,n]=size(a);

b=[];

for i=1:m

for j=1:n

b(i,j)=a(i,j)/ sqrt(sum(a(:,j).^2));

end

end

c=[];d=[];

for j=1:n

for i=1:1

c(i,j)=max(b(:,j)) ;

d(i,j)=min(b(:,j));

end

end

e=[];f=[];g=[];

for i=1:m

e(i,1)=sqrt(sum(b(i,:)-c)^2);

f(i,1)=sqrt(sum(b(i,:)-d)^2);

g(i,1)=f(i,1)/(e(i,1)+f(i,1));

end

g

Z:

1

min=(@sqrt(x^2+(3-y)^2)+y+@sqrt((y1-y)^2+(15-x)^2))*9.5+(@sqrt((25-15)^2+(8 -y1)^2)*(9.5+23.54));

x>=0;

x<=15;

y>=0;

y1>=0;

y1<=8;

Global optimal solution found.

Objective value: 502.6264

Objective bound: 502.6264

Infeasibilities: 0.000000

Extended solver steps: 767

Total solver iterations: 56971

Variable Value Reduced Cost

X 4.427985 0.000000

Y 0.4435016 -0.3828085E-08

Y1 6.547257 0.1127801E-07

Row Slack or Surplus Dual Price

1 502.6264 -1.000000

2 4.427985 0.000000

3 10.57201 0.000000

4 0.4435016 0.000000

5 6.547257 0.000000

6 1.452743 0.000000

Z:

2

min=@sqrt(15^2+(3-y1)^2)*9.5+(@sqrt((x-15)^2+(y1-y)^2)+y+@sqrt((25-x)^2+(8-y)^2))*(9.5+23.54);

x>=15;

x<=25;

y>=0;

y1>=0;

y1<=8;

Global optimal solution found.

Objective value: 568.4405

Objective bound: 568.4405

Infeasibilities: 0.000000

Extended solver steps: 410

Total solver iterations: 69054

Variable Value Reduced Cost

Y1 0.000000 -1.863103

X 15.00000 0.000000

Y 0.000000 12.40008

Row Slack or Surplus Dual Price

1 568.4405 -1.000000

2 0.000000 0.000000

3 10.00000 0.000000

4 0.000000 0.000000