应急设施的优化选址问题(数学建模)

..第17讲应急设施的优化选址问题问题（AMCM-86B题）里奥兰翘镇迄今还没有自己的应急设施。1986年该镇得到了建立两个应急设施的拨款，每个设施都把救护站、消防队和警察所合在一起。图17-1指出了1985年每个长方形街区发生应急事件的次数。在北边的L形状的区域是一个障碍，而在南边的长方形区域是一个有浅水池塘的公园。应急车辆驶过一条南北向的街道平均要花15秒，而通过一条东西向的街道平均花20秒。你的任务是确定这两个应急设施的位置，使得总响应时间最少。图17-11985年里奥兰翘每个长方街区应急事件的数目（I）假定需求集中在每个街区的中心，而应急设施位于街角处。（II）假定需求是沿包围每个街区的街道上平均分布的，而应急设施可位于街道的任何地方。§1若干假设1、图17-1所标出的1985年每个长方形街区应急事件的次数具有典型代表性，能够反映该街区应急事件出现的概率的大小。2、应急车辆的响应时间只考虑在街道上行驶时间，其他因纱（如转弯时间等）可以忽略不计。3、两个应急设施的功能完全相同。在应急事件出现时，只要从离事件发生地点最近的应急设施派出应急车辆即可。4、执行任何一次应急任务的车辆都从某一个应急设施出发，完成任务后回到原设施。不出现从一个应急事件点直接到另一事件点的情况。（这是因为，每一个地点发生事件的概率都很小，两个地点同时发生事故的概率就更是小得可以忽略不计）。..§2假定（I）下的模在假定（I）下，应急需求集中在每个街区中心。我们可以进一步假定应急车辆只要到达该街区四个街角中最近的一个，就认为到达了该街区，可以开始工作了。按假定（I），每个应急设施选在街角处，可能的位置只有6×11=66个。两个应急设施的位置的可能的组合至多只有66×65/2=2145个。这个数目对计算机来说并不大，可用计算机进行穷举，对每种组合一一算出所对应的总响应时间，依次比较得出最小的响应时间及对应的选址方案。具体算法是：建立直角坐标系，以该镇的西北角为原点，从北到南为X-轴正方向，从西到东为Y-轴正方向，在南北、东西方向上分别以一个街区的长作为单位长，则街角的坐标),(YX是满足条件50,100YX的整数。而每个街区中心的坐标具有形式)5.0,5.0(ji，其中ji,是满足条件：40,90ji的整数。如果不考虑障碍和水塘的影响，同应急车辆从设在),(YX点的应急设施到以)5.0,5.0(ji为中心的街区的行驶时间等于)5.05.0(20)5.05.0(15),,,(jYiXjiYXt)5.17)5.0(20)5.0((15jYiX秒记),(jip为以)5.0,5.0(ji为中心的街区的事故发生频率（即在图上该街区所标的数字）。如果应急设施设在),(),,(2211YXYX这两点，总不妨设21XX，则该设置方案的总响应时间为),,,(2211YXYXT90402211)},,,(),,,,(min{),(ijjiYXtjiYXtjip让1X取遍0—10，2X取遍101X，21,YY分别独立地取遍0—4。依次对四数组),,,(2211YXYX的每一个值算出对应的总响应时间的最小值及对应的四数组。以上算法不难用计算机编程实现。由于数组的个数不算多（只有两千多个），计算机可很快得出答案。答案是：两个应急设施分别设在点（2，3），（6，3）时最优。这是在不考虑L形障碍区域和水塘的影响的假定下得出的最优解，但从这两个点到..任何街区都可避开L形障碍区域和水塘，故它们也就是原题所需的最优选址。§2假定（II）下的模型在假定（II）下，由于允许应急设施设在街道上任何位置，这就有无穷多种可能位置，不能直接用计算机穷举。不过，我们可证明：应急设施仍应设在街角处，才能使总响应时间最少。对已选定的两个应急设施的位置A和B，我们先来看总响应时间怎样计算。首先，我们将街道上所有的点的集合划分成两个责任区BAVV,，分别由BA,进行救助：街道上的点P如果由A点去救助比由B点去救助的路程更近，就将P划进A的责任区AV，反之就划进BV，为叙述方便，我们将每个长方形街区的四条边中的每一条称为一条“街道”，街道的一段称为“街段”。每条街道中属于AV的点与属于BV的点各组成一个街段，分别称为A的或B的“责任段”。一条街道最多被分成两个责任段（也有可能整条街道属于同一个责任区，因而本身就是一个责任段），责任地段只有有限多条，对每个应急设施，我们分别算出它的每个责任段的总响应时间，将这些总响应时间求和就得到这个设施的责任区的总响应时间。将两个责任区各自的总响应时间相加就得到这一选址方案的总响应时间。下面需要知道：任一设施A到它的一个责任段EF的总响应时间怎样计算。按假定（II），街区出现事故的频率平均分布在它周围的四条街道上，每条街段的事故发生频率与它的长度成正比。将应急车辆每秒钟行驶的路程作为长度单位，则当街区事故频率为p、街段的长度为t时，这一街段的事故频率为70,70/tp是街区的周长，即车辆绕街区行驶一周需70秒。在大多数情况下，一条街段同时与两个街区相邻，两个街区的事故它都有份，它的事故频率应为qptqp、,70/)(分别是两个街区的事故的总频率（即原题图上标出的数）。当然可以用积分的方法。即插入分点将责任段EF分成许多微小街段i，对每一小段i按其长度计算出它的事故发生频率iikdsp，其中ids是i的长度，k是与i无关（但与EF的选取有关）的常数。取应急车辆人A到i中任意一点的行驶时间iT作为A到i的时间，则微小街段i的响应时间近似地等于iidsT。对这些微小的响应时间求和即得到EF的总响应时间的近似值。让每个0ids，求和变成求积分即可。但在这里，问题比较简单，可以不用积分。事实上，由于EF的每一小段的事故..发生频率只与这一小段的长度有关，换句话说：频率密度是常数，只要求出EF到A的平均行驶时间T，再乘以EF的总的事故频率就行了。当A设在街角处时，平均行驶时间也就是A到EF的中点M的行驶时间212015mYmXTMA秒，这里),(),,(21mmYX分别是MA,的坐标，而且不考虑障碍和水塘的影响。将MAT乘以EF的事故频率，就得到EF的总响应时间。换句话说，就是将EF的事故频率EFP集中到M点，认为M按频率EFP发生事故，而EF的其他点都不发生事故。这样不会改变EF的总响应时间，却便于计算，如果应急设施A不是设在街角处，而是设在某条街道CD的两个端点DC、之间，则可能出现这样的情况：从A出发到EF中的某些点的最短救助路线应向C方向行驶，崦到另一些点去则应向D方向行驶。这时，平均时间就不等于A到EF中点M的时间AMT，而是比AMT小。在这样的情况下EF可以分成两段GFEG、，从A到其中一段（比如EG）上的所有的点的最短救助路线应向C方向行驶，而到另一段（比如GF）上的所有的点的点则应向D方向行驶。分别计算GFEG、的事故发生频率GFEGPP,，将这两个频率分别集中在GFEG、各自的中点21,MM，就可分别算出GFEG、的总响应时间，再将它们相加就得到EF的总响应时间。下面证明：最短的总响应时间必可由设在街角处的应急设施BA、来实现。假定已选择两个应急设施BA、的位置使总响应时间最短，且至少有一个设施（比如A）不是设在街角处，而是设在某一条街道CD的两个端点DC、之间。我们证明：可以把这个设施从A移到C或D，使总响应时间不增加，（而且很可能减少）。证明的主要想法是：将设施迁移到街角后，它到某些街段缩短了一段路程，同时到另外某些街段增加同样长的一段路程。如果路程缩短的那些街段的事故总频率大于路程增加的那些街段的事故总频率，则总响应时间缩短了，设施位置得到优化，说明原来的位置不是最优。先考虑与街道CD相邻的街区，也就是与急救站A相邻的街区。要使总响应时间最少，两个急救站BA,的位置显然不应当靠得太近。因此，可以假定与A相邻的街区周界上所有的点到A的路程都小于它们到B的路程，因而都应当由A负责救助。这个街区的事故频率p均匀分布在街区的周界上。我们指出：救助这个街区的事故频率p均匀分布在街区的周界上。我们指出：救助这个街区的事故的总响应时间与A在CD上的位置选取无关。事实上，无论A处于街道CD上哪一个位置，总存在一点A将街区周界分成路程相等的两段，第一段由A经C到A，第二段由A经D到A，每一段的总行驶时间是7/2=35秒，事故..总频率是2/p。由A出发去救助每一段上各点的平均行驶时间等于35/2秒，因而两段的总响应时间为2)2/35()2/(p秒，确实与A点位置的选取无关。因此，在讨论A在CD上的位置选取时，不需考虑到CD相邻的街区的事故的影响，不妨暂时假定这样的街区的事故频率为0，特别是街道CD上不发生事故，不需要救助。设P是A的责任区AV内需要救助的任一点，从A出发到P，有两种可能的最短救助路线AP：一种是沿AC、经由C点到P，另一种是沿AD、经由D点到P。凡是AP属于前一种情况的，这样的点P组成的集合记作CU；凡是AP属于后一种情况的，这样的点P组成的集合记作DU。这样就将A的责任区按最短救助路线出发时的两个不同方向分成了两个区域（各由一些街段组成）。比较DCUU,这两个区域各自的事故总频率DCPP,的大小。如果CP比DP大，我们就将设施从A移到C，向CU靠拢（同时远离DU）；反之，当DP比CP大时，将设施由A迁到D去靠近DU（同时远离CU）；当DCPP时将设施任意迁到C或D都可以。我们证明：将设施经过这样的迁移后，总响应时间只可能减少，不可能增加。因此，假如迁移前的方案最优，迁移后一定还是最优（事实上，当DCPP时，迁移后的方案一定比原来更优，说明原来不可能最优）。不妨先假定DCPP，设施从A迁到C点。（CDPP的情况同理）。为了便于比较迁移前后的总响应时间的变化情况，我们先作下面两个假设（其中所说的“旧”是指设施迁移前的情况，而“新”则是指迁移后的情况）：（1）应急设施从A搬透到C后，两个旧的责任区BAVV,先仍分别由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，总响应时间减少DCUTP;的点则相反，路程都增加AC，行驶时间都增加T，总响应时间增加TPD。由于TPTPPPDCDC,，总响应时间减少量超过（或等于）增加量，总的效果是减少了（或不改变）总响应时间，设施搬迁后的位置比原来更优，至少同样优。..假设（2）的路线不一定是最短路线。如果再进一步选择最短路线，则还有可能进一步缩短新设施方案的总响应时间，更加说明其优越性，假设（1）的责任区的贡分不一定是合理的，可以再进行调整，将街道上的每一点划给离它最近的设施的责任区，这样又可能再减少新设施方案的总响应时间，再一次增加它的优越程度，这样就证明了新设施比旧设施更优，或同样优。因此，在假定（II）下，仍可设应急设施设在街角处。于是与假定（I）的情况类似地可用计算机穷举算出答案来，对任一对候选的应急设施位置),(),,(2211YXBYXA，（坐标为整数），求出每一条街道CD的总响应时间，将所有街道的总响应时间相加就得到这一选址方案的总响应时间。进行比较就可得出最短的总响应时间及对应的选址方案。CD的总响应时间的计算方法已在前面讲过。并且由于设施都设在街角处，只要将CD分成两个责任段（在多数情形下实际上只有一个责任段）EDCE,，将这两个责任段的事故频率分别集中在它们各自的中点计算就可以了。计算结果：应急设施以设在点（2,3）,(7,3)时最优