救护中心建立问题的研究
摘要
本文对某小镇建立两个救护中心,使应对突发事件总的响应时间最少的问题进行了
分析,并建立了数学模型进行了求解。
在假设(I)的前提下,即需要救护的事件集中在每个街区的中心。考虑到街区数目
不是很多,本文采用穷举法进行了最优解的搜索。即先任意选取两点作为救护中心的位
置,然后计算其他街区到这两个救护中心的总响应时间,总响应时间最少的旧最优的方
案。同时为了考虑障碍区域和水塘,本文首先对那些设置救护中心需要穿越障碍区域和
水塘的点进行了剔除,然后在利用计算机一一穷举。
在假设(Ⅱ)的前提下,需要救护的事件均匀分布在街道上,在计算总响应时间时,
本文把整个街道的事件发生频率集中在街道的中心位置处进行计算。同时本文证明了当
救护中心仍设立在街角处时所需的总响应时间是最少的,这样仍可以按照假设(I)中的
穷举方法求出救护中心设立的最优位置。
关键词:穷举法;剔除;街道中心;街角
一. 问题的重述
某小镇开始计划建立两个救护中心,把救护站、消防队和派出所结合在一起。
图 1 指出每个长方形街区所发生的需要救护事件的次数,北边的 L 形区域是障碍,
而南边的长方形区域是浅水池,救护车辆驶过一条南北向的街道平均花 15 秒,而
救护车辆驶过一条东西向的街道平均花 20 秒,请确定这两个救护中心的位置,使
得总响应时间最少。
(1)假定需要救护的事件集中在每个街区的中心,救护中心位于街角处。
(2)假定需要救护的事件沿包围每个街区的街道上均匀分布,救护中心可位于街道
的任何地方。
图 1 小镇的街区分布图
二.问题分析
对于假设(I)的情况,要建立救助站的位置,使总的响应时间最短。在考虑障碍区
域的情况下,可以首先把那些建立救护站需要穿过障碍区域的点剔除掉 ,然后可以考
虑穷举法利用计算机求出最佳的建立救护中心的位置。
对于假设(Ⅱ)的情况,由于突发事件是均匀分布在每条街道上的,可以利用每条街
道的中心点位置来作为这整条街道突发事件的频率集中点。同时可以证明:在街角处设
置救护中心是所需总响应时间最短的。这样建立救护中心的位置使有限的,仍可以考虑
穷举法利用计算机进行求解。
1
T(X ,Y ,X ,Y )
2
2
1
t(X ,Y , , )
n i j
m
P
( , )
P i
j
i
j
三.符号说明
建立两个救护中心所需最少时间的函数
第(i,j)个小区到( X ,Ym
街区的事故发生频率矩阵
n )个救护站的时间函数
坐标中第(i,j)个区域的应急事件次数
每个街区街角的横坐标
每个街区街角的纵坐标
四.模型假设
1.在同一时刻各街区及街区内不会同时发生两件应急事件;
2.应急事件集中在街区中心,而应急设施在街角处;
3.障碍区及浅水塘区应急事件次数为零;
4.应急车辆的响应时间只考虑在街道上行驶的时间,不考虑转弯的时间;
5.假设应急需求沿包围每个街区的街道是均与分布的,而应急措施可以位于街道的
任何地方。
5.1 问题一的模型建立与求解
五.模型的建立与求解
在假定(I)的前提下,应急需求集中在每个街区的中心,本文进一步假定应急车辆
只要到达该街区四个街角中最近的一个,就认为到达率该街区。则有假定(I)每个应急
设施选在街角处,可能的位置只有 6×11=66 个,两个应急设施的位置的可能组合至多
只有 66×65/2=2145 个。考虑到这个数目不是很大,可以利用计算机进行穷举,对每种
组合一一算出所对应的总响应时间,依次比较得出最小的响应时间及其对应的选址方
案,具体的算法是:
建立直角坐标系,一该镇的东北角为原点,以从东到西为 X 轴的正方向,以从北向
南为 Y 轴的正方向,在南北和东西方向分别以一个街区的长作为单位长,则街角的坐标
(X,Y) 是 满 足 条 件 0
4,0
(
i
塘的影响则应急车辆从设在(X,Y)点的应急到一 (
i
的 整 数 。 而 每 个 街 区 的 中 心 的 坐 标 具 有 形 式
的整数。如果不考虑障碍和水
,其中 i,j 满足的条件:0
为中心的街区的行驶时
9
j
0.5,
j
0.5)
0.5)
0.5,
5,0
10
X
Y
j
i
间等于:
(
t X Y i
, , )
j
,
j 为以 (
i
记 ( , )
P i
0.5
0.5) 15
20(
X i
(
20
X
i
0.5,
0.5)
j
为中心的街区的事故发生频率(即在图上该街区所标的数
0.5) 15(
Y
(
j
0.5) 17.5
0.5
Y
0.5)
j
字)。用矩阵 P 表示,则 P 为
P
5 2 2 1 5 0 3 2 4 2
2 3 3 3 3 4 1 3 0 4
4 3 3 0 3 4 0 0 0 0
1 2 0 0 4 3 2 2 0 1
3 3 2 5 3 2 1 0 3 3
则如果应急设施设在 1
1
(
X Y ,
)
,
(
X Y 这两点,则该设置方案的总响应时间为
)
,
2
2
,
T X Y X Y
2
(
,
,
1
1
2
)
4
9
i
0
j
1
( , ) min{ (
P i
j
,
t X Y t X Y
2
), (
,
1
1
2
)}
2,Y Y 分别独立地取遍 0—10,依次对 1
)
,
让 1X 取遍 0—4, 2X 取遍 1X —4, 1
X Y X Y
2
(
T X Y X Y 。通过比较可得出所有这些总响应
(
)
,
,
,
,
,
1
2
1
2
2
的每一个值计算出对应的总响应时间 1
时间的最小值。
以上计算的过程没考虑障碍区域和水塘的情况,对于考虑障碍区域和水塘的情况,
见图 2.
图 2.障碍区域
只要不把救护站设在图 2 所标记的 A、B、C、D、E 等点,则救护的路线就不会穿越
障碍区域。所以在计算中可以事先把这些点剔除。同理对于水塘障碍区域见图 3.只要剔
除 F、G 点即可保证救护车的路线不会穿越水塘。
图 3. 水塘区域
按照以上计算过程利用计算机求解可得在 M 和 N 处建立救护中心,所需的最少时间
为 5065 秒。见图 4.
图 4.救护中心的位置
按照以上计算过程利用计算机求解可得在 M 和 N 处建立救护中心,所需的最少时间
为 5065 秒。即两个救护中心分别建立在(2,2)和(2,6)两点。见图 4.(求解代码见附件)
5.2 模型二的建立与求解
在假定(Ⅱ)的情况下,由于允许应急设施在街道上任何位置,这就有无穷多种可能,
不能直接用计算机穷举,不过,本文可以证明:应急设施仍应设在街角处,才能使总的
响应时间最短。
对已选定的两个应急设施的位置 A 和 B。首先,本文将街道上所有的点的集合划分
成两个责任区 VA 和 VB ,分别由 A,B 进行救助:街道上的点 P 如果由 A 点去救助比由 B
点去救助的时间更短,就将 P 划进 A 的责任区 VA ,反之就划进 VB 。为了叙述方便,本
文将每个长方形街区的死条边中的每一条称为一条“街道”,街道的一段称为“街段”。
没挑街道中属于 VA 的点与属于 VB 的点各组成一个街段。分别成为 A 的或 B 的责任段。
一条街道最多被分成两个责任段,责任段只有有限多条。对于每个应急设施,分别计算
出它的每个责任段的总响应时间,将总响应时间求和就得到这个设施的责任区的总响应
时间,将两个责任区各自的总响应使劲相加就得到这一选址方案的总响应时间。
下面需要知道:任一设施 A 到它的一个责任区段 EF 的总响应时间怎样计算。按假
定(Ⅱ),街区出现事故的频率平均分布在它周围的四条街道上,每条街段的事故发生频
率与它的长度成正比。将应急车辆每秒钟行驶的路程作为长度单位,则当街区事故频率
为 p,街段的长度为 t 是,这一街段的事故频率为 p×t/70,70 是街区的周长,即车辆
绕街区行驶一周需 70 秒,在大多数情况下,一条街段同时与两个街区相邻,两个街区
的事故它都有份,它的事故频率应为(p+q)×t/70,p,q 分别是这两个街区的事故的总频
率。事实上,由于 EF 的每一小段的事故发生频率只与这一小段的长度有关,换句话说:
频率密度是常数,只要求出 EF 到 A 的平均行驶时间 T,再乘以 EF 的总的事故频率就行
了 。 当 A 设 在 街 角 处 时 , 平 均 行 驶 时 间 也 就 是 A 到 EF 的 中 点 M 的 行 驶 时 间
秒,这里(X,Y),( 1m , 2m )分别是 A,M 的坐标,将 MAT 乘以
MAT
EF 的事故频率,就得到 EF 的总响应时间。若应急设施 A 不是设在街角处,而是设在某
X m
1
Y m
15
2
20
条街道 CD 的两个端点 C、D 之间,则可能出现这样的情况:从 A 出发到 EF 中的某些点
的最短救助路线应该向 C 方向行驶,而到另一些点去则应向 D 方向行驶。这时,平均时
间就不等于 A 到 EF 中点 M 的时间 MAT ,而是比 MAT 小。在这样的情况下 EF 可以分成两段
EG、GF,从 A 到其中一段上的所有的点的最短救助路线应向 C 方向行驶,而到另一段上
的所有的点则应向 D 方向行驶。分别计算 EG、GF 的事故发生频率 EGP , GFP ,将这两个
频率分别集中在 EG、GF 各自的中点 1M , 2M ,就可分别算出 EG、GF 的总响应时间,再
将他们相加就得到 EF 的总响应时间。
下面证明:最短的总响应时间必可又这在街角处的应急设施 A、B 来实现,假定已
选择两个应急设施 A、B 的位置使总响应时间最短,且至少有一个设施(比如 A)不是设在
街角处,而是设在某一条街道 CD 的两个端点 C、D 之间。本文证明:可以吧这个设施从
A 移到 C 或 D,使总响应时间不增加。证明的主要想法是:将设施迁移到街角后,它到
某些街段缩短了一段路程,同时到另外某些街段增加同样长的一段路程。如果路程缩短
的那些街段的事故总频率大于路程增加的那些街段的事故总频率,则总的响应时间缩短
了,设施位置得到优化,说明原来的位置不是最优。设 P 是 A 的责任区 AV 内需要救助
的一点。从 A 出发到 P,有两种可能的最短救助路线 AP:一种是沿 AC、经由 C 点到 P;
另一种是沿 AD、经由 D 点到 P。凡是 AP 属于前一种情况的,这样的点 P 组成的集合记
作 CU ;凡是 AP 属于后一种情况的,这样的点 P 组成的集合记作 DU 。这样就将 A 的责
任区按最短救助路线出发时的两个不同方向分成了两个区域。比较 CU , DU 这两个区域
各自的事故总频率 CP , DP 的大小。如果 CP 比 DP 大,就将设施从 A 移到 C,向 CU 靠拢;
反之,当 DP 比 CP 大时,将设施由 A 迁到 D 去靠近 DU ;当 CP = DP 时将设施任意迁到 C 或
D 都可以。本文将证明:将设施经过这样的迁移后,总响应时间只可能减少,不可能增
P
加。因此假如迁移前的方案最优,迁移后一定还是最优(事实上,当 C
P
的方案一定比原来最优,说明原来的不可能最优)。不妨先假定 C
C 点为了便于比较迁移前后的总响应时间的变化情况,本文先作下面两个假设:
P 时,迁移后
P ,设施从 A 迁到
D
D
(1) 应急设施从 A 搬迁到 C 后,两个旧的责任区 AV , BV 先仍分别由 C 和 B 负责救
助,暂时不变。如果在这样不改变责任区的情形下都能证明总响应时间不增加,则
再进一步合理调整 C、B 的责任区还可能进一步缩短(至少不会增加)总响应时间。
(2) 搬迁后从新设施 C 到旧区域 CU 中的任何一点 P 的救助路线为:从 C 出发离开
CD,沿原先 A 的旧的救助路线到 P。从 C 到旧区域 DU 的任何一点 P 的救助路线为:
从 C 出发沿 CD(经过 A)到 D,再沿原先 A 的旧的救助路线到 P。
设应急车辆从 A 到 C 的行驶时间为 T。则按(2)的行驶路线, CU 的点到设施的路
程减少了 AC,行驶时间减少 T,总响应时间减少 CP T ; DU 的点则相反,路程都增
加了 AC,行驶时间都增加 T,总响应时间增加 DP T 。由于 C
P
,总
响应时间减少量超过(或等于)增加量,总的响应时间是减少了(或不改变),设施搬
P T P T
D
P , C
D
迁后的位置比原来更优,至少同样优。
假设(2)的路线不一定是最短路线,如果再进一步选择最短路线,则还有可能
进一步缩短新设施方案的总响应时间,更加说明其优越性。假设(I)的责任区的划
分不一定是合理的,可以再进行调整,将街道上的每一点划给离他最近的设施的责
任区,这样又可能再减少新设施方案的总响应时间,再一次增加它的优越程度。这
样就证明了新设施比就设施更优,或者同样优。
因此在题目假定(Ⅱ)之下,仍可设应急设施设在街角处,于是与假定(I)的情
况 类 似 地 可 用 计 算 机 穷 举 算 出 答 案 来 。 对 任 一 对 候 选 的 应 急 设 施 位 置
2X ,Y ),求出每一条街道 CD 的总响应时间,将所有街道的总响应时间
A( 1
相加就得到这一选址方案的总响应时间。进行比较就可得出最短的总响应时间及相
1X ,Y ),B( 2
对的选址方案。CD 的总响应时间的计算方法已在前面讲过。并且由于设施都设在街
角处,只要将 CD 分成两个责任段 CE、ED,将这两个责任段的事故频率分别集中在