Optimization of subway platform's passenger flow guidance methodbased on AnyLogic
-
摘要: 为避免地铁车厢客流过度集中,引导站台乘客分散乘车,文章以地铁各车厢客流均衡为目标,以均衡度标准差为指标,综合考虑车厢现有人数、候车人数及预计下车人数,优化地铁站台客流引导方式。通过AnyLogic仿真软件,模拟广州地铁珠江新城5号线站台条件,对比了客流量范围在12万人/日~60万人/日条件下,不同引导方式对车厢乘客分布不均衡度的影响。结果表明,地铁站客流量达到24万人/日时,考虑预计下车人数的站台客流引导方式下的车厢客流均衡度标准差明显提升;与不考虑该站下车人数的引导方式相比,客流均衡度标准差提升幅度最大达到25.6%,相比人工引导方式的提升幅度达53.6%,有效提升了地铁精细化服务能力。
-
关键词:
- 城市轨道交通 /
- 地铁站台 /
- 客流引导 /
- AnyLogic仿真 /
- 均衡度标准差
Abstract: In order to avoid excessive concentration of passenger flow in subway carriages and guide passengers on the platform to disperse and board trains, this paper aimed to balance the passenger flow in each carriage of the subway, used the standard deviation of balance degree as an indicator, comprehensively considered the existing number of passengers in the carriage, the number of waiting passengers, and the expected number of alighting passengers, to optimize the subway platform passenger flow guidance method. Through AnyLogic simulation software, the paper simulated the platform conditions of the Pearl River New Town Line 5 of Guangzhou Metro, and compared the influence of different guidance modes on the uneven distribution of passengers in the carriage under the conditions of passenger flow ranging from 120 000 people/day to 600 000 people/day. The results show that when the passenger flow of the subway station reaches 240 000 people/day, the standard deviation of the passenger flow balance of the carriages under the platform passenger flow guidance method considering the expected number of passengers getting off significantly improves. Compared with the guidance method that does not consider the number of passengers getting off at the station, the standard deviation of passenger flow balance has increased by a maximum of 25.6%, and compared to the manual guidance method, the improvement has reached 53.6%. The optimized guidance method effectively improves the refined service capacity of the subway. -
随着我国铁路建设的不断推进,高速列车在运行中能耗问题的重要性愈发彰显。列车在运行过程中有4种工况,能源消耗主要集中在牵引和制动这两种工况。列车运行曲线对应不同工况组合下的列车能耗,因此研究能耗最小的列车运行曲线具有重要的现实意义。
目前,众多学者对能耗优化问题进行了研究,文献[1]将列车节能问题转换为一定约束条件下的寻优问题;文献[2]利用极大值原理求解列车转换点;文献[3]提出通过划分坡道求解列车转换原则和运行位置的关系;文献[4]引入遗传算法,优化列车运行曲线的惰行距离,通过调整惰行点实现节能效果;文献[5]基于遗传算法,讨论列车运行曲线优化及晚点恢复等问题;文献[6]利用遗传算法解决列车运行调整问题;文献[7]利用遗传算法对列车运行过程中最优的组合序列进行求解。由于遗传算法的固有缺陷,一些学者开始研究针对列车节能问题的算法改进。文献[8]利用改进后变长算子的遗传算法,对城轨列车节能运行问题进行优化;文献[9]通过在算法中加入准点调整和局部搜索的引导机制对传统遗传算法进行改进。
本文在前人研究的基础上,通过建立约束条件下的列车运行模型和对遗传算法中算子的改进,增强算法的全局搜索能力和收敛能力,同时,引入自适应机制,并通过Matlab仿真,验证了改进后算法的有效性。
1 列车运行节能模型及策略
1.1 列车运行节能模型
列车在运行过程中受到牵引力、制动力、运行阻力的共同作用[10],根据牛顿第二定律,对运行列车进行受力分析可得,列车在正常行驶过程中的动力学方程为:
$$\frac{{dv(x)}}{{dt}} = \frac{{F(v) - B(v)}}{M} - {w_0}(v) - {w_f}(x)$$ (1) 其中,
$v$ 是列车当前运行速度;$t$ 是列车行驶时间;$F(v)$ 为列车在当前速度下的牵引力;$B(v)$ 为列车在当前速度下的制动力;${w_0}(v)$ 为列车在当前速度下的单位运行基本阻力,计算方式如式(2)所示;${w_f}(x)$ 为列车当前位置运行单位附加阻力;$M$ 为列车总质量;$x$ 为列车当前运行里程数。$${w_0}(v) = a + bv + c{v^2}$$ (2) 其中,
$a、b、c$ 均为常数,一般根据运行车型来取值。列车运行单位附加阻力
${w_f}(x)$ 包含单位坡道附加阻力${w_p}(x)$ 、单位曲线附加阻力${w_q}(x)$ 、单位隧道附加阻力${w_s}({x})$ 。由于列车长度较长,计算阻力时将其视为单一质点并不可靠,因此考虑将列车的每节车厢视为单质点,来计算列车单位附加阻力,而列车所受的总附加阻力为各个车厢所受附加阻力之和,如公式(3)所示:$$\begin{array}{l} {w_f}(x) = \displaystyle\sum\limits_{n = 1}^m {[{w_p}({x_n}) + {w_q}({x_n}) + {w_s}({x_n})]} \\ \left\{ \begin{array}{l} {w_p}\left( {{x_n}} \right) = p\left( {{x_n}} \right)\\ {w_q}\left( {{x_n}} \right) = \dfrac{{600}}{{R\left( {{x_n}} \right)}}\\ {w_s}\left( {{x_n}} \right) = 0.000\;13 \cdot {L_s}\left( {{x_n}} \right) \end{array} \right. \end{array}$$ (3) 其中,
${w_p}({x_n})$ 为第$n$ 节列车在$x$ 处的单位坡道附加阻力;$p({x_n})$ 是该位置的坡道千分数;${w_q}({x_n})$ 是第$n$ 节列车在x处所受到的单位曲线附加阻力;$R({x_n})$ 是第$n$ 节列车在$x$ 位置处的曲线半径;${w_s}({x_n})$ 是第$n$ 节列车在$x$ 处的单位隧道附加阻力;${L_s}({x_n})$ 是第$n$ 节列车在$x$ 位置处的轨道长度。线路附加阻力在同一位置是固定的,为方便计算,考虑提前处理线路附加阻力,并将其存储起来,后续使用时根据车头与列车各车厢相对距离判断所处区间,从而得到附加阻力。
将列车运行能耗最小作为优化目标,同时需要满足准时、定点等安全约束条件,建立列车节能运行数学模型[11]如下:
$$\begin{array}{l} \min E = \min \displaystyle\int_0^T {F(t) \cdot v} (t)dt\\ s.t.\left\{ \begin{array}{l} v = \dfrac{{ds}}{{dt}},\\ v(0) = v(S) = 0,\\ v(x) \leqslant {v_{\lim }}(x),\\ \displaystyle\int_0^S {v{{(x)}^{ - 1}}dx} = T,\\ t(0) = 0,t(S) = T \end{array} \right. \end{array}$$ (4) 其中,
$E$ 是列车从出发到终点所消耗的能量,本文只考虑列车牵引能耗;$F(t)$ 是$t$ 时刻的牵引力;$v(t)$ 是$t$ 时刻列车的速度;$S$ 是出发站到终点站的距离;$T$ 是列车总运行时间;${v_{\lim }}(x)$ 是当前路段的最大速度限制。1.2 列车运行节能策略
列车在同一路线中采取不同的运行策略,运行能耗也不同。本文将(xi , Mqi)称为工况转换点,其中,xi为列车工况切换点的位置;Mqi为对应工况。根据列车节能操纵原则[12],列车起动时以最大牵引力起动,制动时以最大制动力停车,有利于节能。如图1所示,列车运行工况可分为牵引(
${{\rm{SQ}}}$ )、匀速(${{\rm{SY}}}$ )、惰性(${{\rm{SD}}}$ )和制动(${{\rm{SZ}}}$ ),列车以${\rm{SQ}}$ 工况从${{\rm{s_1}}}$ 处加速行驶至${{\rm{s_2}}}$ ,切换成${\rm{SY}}$ 工况匀速行驶至${{\rm{s_3}}}$ ,再切换为${\rm{SD}}$ 工况至${{\rm{s_4}}}$ 处,接着继续以${\rm{SY}}$ 工况行驶,到${{\rm{s_5}}}$ 转换为${\rm{SZ}}$ 工况直至终点。$[({{\rm{s}}_{\rm{1}}},{\rm{SQ}}),({{\rm{s}}_2},{\rm{SY}}),$ $({{\rm{s}}_{\rm{3}}},{\rm{SD}}),({{\rm{s}}_{\rm{4}}},{\rm{SY}}), $ $ ({{\rm{s}}_{\rm{5}}},{\rm{SZ}})]$ 称为列车的一个操纵序列可行解。实际运行中存在着无数种工况转换点的组合,因此需要找出所有可行解中能耗最低的运行策略。由本文设置的列车操纵序列推导出的列车运行曲线未必全符合运行条件,当出现停车位置或运行时间超过一定范围的不可行解时,则不予保留。本文停车精确度取10 m,运行时间精确度取60 s。
1.3 制动点设置
列车制动时按最大减速度
${a_{\max }}$ 制动,制动曲线从停车点开始反算。如图2所示,当列车运行曲线与停车制动曲线相交时,可切换为制动工况,因此只要制动曲线落于${S_{QMNL}}$ 内且与$x$ 轴相交,则可满足停车精确度要求。2 改进遗传算法求解
数值法求解列车节能问题模型较复杂,而遗传算法适用于解决复杂情况的寻优问题[13],但传统的遗传算法存在易早熟和收敛性差的问题,因此本文考虑改进算法中的遗传算子,并引入进化逆转算子。
2.1 算法概述
(1)编码
列车运行策略是由一系列上文讨论过的工况组成,将工况转换点
$({{{x_i}}},{{{M_{qi}}}})$ 做为基因,运行区间内所有工况转换点的集合做为染色体。本文采用实数编码。染色体结构如图3所示,其中$ l $ 为染色体长度。(2)种群初始化
根据线路情况和列车运行工况设置原则,随机生成原始工况转换点位置,初始种群规模为n。本文设置初始种群规模为100。
(3)适应值计算
列车节能优化的目标函数为
$\min E$ ,为方便计算,取目标函数的倒数为适应度值$fitness = \dfrac{1}{{\min E}}$ 。适应度越高的染色体越接近理想求解。(4)收敛条件
以是否达到最大进化次数作为收敛条件,如果没有达到条件,则继续进化,否则退出。本文设置最大进化代数为200。
2.2 自适应遗传进化算子
采用轮盘赌法选择染色体以保证在进化过程中适应度较大的个体有较大的概率,进入下一代种群。
2.2.1 改进交叉算子
通常,在算法初期普遍种群适应度值较低,染色体质量较差,后期染色体适应度值和质量都会相应提升,因此,考虑将普通遗传算法中交叉算子的概率改变为随着迭代次数变化的动态值:
$$P_c^m = \left\{ \begin{array}{l} {P_{c,\max }} \cdot {e^{ - \frac{m}{M}}},{\text{若}}{P_{c,\max }} \cdot {e^{ - \frac{m}{M}}} < {P_{c,\min }}\\ {P_{c,\min }},{\text{其他}} \end{array} \right.$$ (5) 其中,
$m$ 表示当前进化代数;$M$ 为总进化代数;${P_{c,\max }}$ ,${P_{c,\min }}$ 分别是预先设定的最大交叉概率和最小交叉概率。交叉概率随着进化代数的增加而变小,以便初期能够得到更多的新染色体,增大解空间,提高算法全局搜索能力,后期交叉概率相应减小,保证算法稳定收敛。2.2.2 改进变异算子
实际运行中,一般加速转匀速的工况转换点不会轻易变化,只有
${{({{{x}}}_{2},\;{{{x}}}_{3},\;{{{x}}}_{7},\;{{{x}}}_{8})}}$ 这类工况转换位置可以改变,而传统的变异算子随机性过强,无法保留上一代进化过程中的重要信息,因此针对这类关键基因点,考虑引入3代评价函数来逼近搜索,得到最佳值。如图4所示。具体操作为,记录下每个关键转换点的更新方向,称为标志位,以及更新后的适应度值。标志位前进为1,后退为−1,若当前标志位连续3代适应度值都增加,则继续在该方向上搜索;若3代中连续2代适应度值减小,则更新标志位方向,即从第1代开始的位置反方向变异搜索。其他变异规则如表1所示。
表 1 3代逼近搜索规则第1代 + + + + – – – 第2代 + + – – + + – 第3代 + – + – – + 下次变异方向 保持当前标志位方向 回朔 保持当前标志位方向 更新标志位方向 回朔 保持当前标志位方向 更新标志位方向 变异算子通过选择染色体中某个基因进行变异,以增强算法的局部搜索能力,提高种群的多样性,防止出现早熟现象。变异概率应随着迭代次数增多而趋小,且适应度较好的个体变异概率也应较小,以保留优良特性,适应度较低个体的变异概率应较大。
$$\qquad\qquad\qquad\qquad\qquad P_b^m = \left\{ \begin{array}{l} {P_{b,\max }} \cdot {e^{ - \left\{ \left(\left| {1 - \frac{{fitness({X_i})}}{{fitnes{s_{\max }}}}} \right|\right) \cdot \alpha \right\} }}, {\text{若}}{P_{b,\max }} \cdot {e^{ - \left\{ \left(\left| {1 - \frac{{fitness({X_i})}}{{fitnes{s_{\max }}}}} \right|\right) \cdot \alpha \right\} }} > {P_{b,\min }}\\ {P_{b,\min }},{\text{其他}} \end{array} \right.$$ (6) 其中,
$fitnes{s_{\max }}$ 为当前种群最大适应度;$fitness({X_i})$ 为当前个体的适应度值;${P_{b,\max }}$ 为预设的最大变异概率;${P_{b,\min }}$ 为预设的最小变异概率;$\alpha $ 为系数调节因子。2.2.3 改进逆转算子
传统遗传算法易陷入局部搜索,且传统交叉算子的随机性易使父代的优秀基因在进化过程中消失。本文结合上述引导机制,为增大解空间,考虑引入逆转算子[14],经逆转操作后,适应度值变大的种群参与下一代进化,适应度趋近或更低的种群则舍弃。算法初期,整体种群适应度较差,多样性较低,逆转概率可以稍大,以增大解空间;算法后期,整体适应度较高,此时逆转操作的概率应当较低以保留优良个体。逆转概率设置如下:
$$\qquad\qquad\qquad\qquad\qquad P_r^m = \left\{ \begin{array}{l} {P_{r,\max }} \cdot {e^{ - \left\{ \left(\left| {1 - \frac{{fitness({X_i})}}{{fitnes{s_{\max }}}}} \right|\right) \cdot \frac{m}{M} \cdot \beta \right\} }}, {\text{若}}{P_{r,\max }} \cdot {e^{ - \left\{ \left(\left| {1 - \frac{{fitness({X_i})}}{{fitnes{s_{\max }}}}} \right|\right) \cdot \frac{m}{M} \cdot \beta \right\} }} > {P_{r,\min }}\\ {P_{r,\min }},{\text{其他}} \end{array} \right.$$ (7) 其中,
${P_{r,\max }}$ 为预设的最大逆转概率;${P_{r,\min }}$ 为预设的最小逆转概率;$\beta $ 为系数调节因子。考虑本文所设置的染色体情况,进行逆转操作时,需要对染色体进行如下处理,工况转换点中运行工况应与逆转前和其位置最相近点的工况保持一致,以防过多不可行解的出现。
3 仿真实例
为验证本文算法的有效性,选定运行线路为合福高铁(合肥—福州)中,合肥至无为的线路,仿真车型为CRH380B。列车共8节编组,线路总长度100.59 km,限速300 km/h,在66~68 km路段限速250 km/h,图定运行时间为1424 s。列车参数如表2所示。列车牵引制动曲线如图5所示。
表 2 列车参数参数名称 参数值 列车总重 495 t 列车长度 200.3 m 最大运行速度 300 km/h 单位阻力 (0.0257+0.050 7 v+0.000 505 57 v2)kN
v为列车当前运行时速定员 510 人 将列车参数及线路信息输入列车运行仿真模型,对算法进行初始化。经过仿真分析,比较正常行驶下的列车运行曲线及根据本文提出的算法优化后的列车运行曲线,如图6所示,可以看出,优化后的列车运行曲线明显增加了惰性区段长度。计算可得优化前后的列车运行数据变化,如表3所示。
表 3 改进前后列车运行数据对比分析评价指标 优化前指标 优化后指标 性能比较 运行时间 1430 s 1459 s +2% 能耗 6.87×10 9 J 6.13×10 9 J −10.7% 由表3可知,对高速列车运行曲线进行优化后,列车运行时间变为1459 s,满足列车准点性要求,而能耗由6.87×109 J降为6.13×109 J,降低了10.7%的能源消耗,本论文提出的算法在满足运行时间指标的情况下有效降低了列车运行中的能耗。
4 结束语
(1)本文基于既有的列车节能操纵策略,建立了以列车安全性、准时性为指标的列车优化模型,讨论了列车停车制动点的设置,为后续列车能耗计算提供基础;
(2)在遗传算法的迭代过程中,提出一种3代逼近搜索的寻优引导机制,解决传统遗传算法求解速度较慢的问题。同时为了扩大解空间,引入逆转算子,进一步提高了算法搜索能力;
(3)以合福高铁区间线路数据和CRH380B型列车为基础进行仿真,结果表明,采用本文提出的方法,列车运行能耗降低10.7%,同时满足列车停车的安全性和准点性,具有较好的理论参考价值。
-
表 1 AnyLogic模型的相关参数设置
模块 参数名称 参数值 TrainSouce 新车厢 carindex==0?new RailCar():carindex%2==0?new PassengerCar2():new PassengerCar1() 车厢长度 carindex==0?21:19 巡航速度 150/35(m/s) Delay 开门等候时间 42.02 s TrainMoveTo 巡航速度 200/45.53(km/h) Dropoff 下车比例 triangular(0, 1, 0.7) Pedsource(S) 到达速率 All*0.132 Pedsource(M) 到达速率 All*0.369 Pedsource(All) 舒适速度 uniform(0.5, 1) 初始速度 uniform(0.3, 0.7) 直径 uniform(0.2, 0.3) PedEnter 舒适速度 uniform(0.7, 1.2) 初始速度 uniform(0.5, 0.9) 直径 uniform(0.1, 0.2) secSelectUnfollow 动态值1 self.TYPE_PROBABILITIES(0.1, 0.1, 0.1, 0.1, 0.3,0.3) 动态值2 self.TYPE_PROBABILITIES(0.1, 0.1, 0.3, 0.3, 0.1,0.1) 动态值3 self.TYPE_PROBABILITIES(0.3, 0.3, 0.1, 0.1, 0.1,0.1) SelectOutput 乘客选择由站台中部通道离开概率(M) 0.47 乘客选择由站台两侧通道离开概率(S) 0.265 表 2 采用引导方式1时乘客分布不均衡状况
客流量/
(万人/日)计数/次 均衡度标准差 偏差 均值置信 总客流量/
人次12 147674 24.574 8.143 0.042 3628992.526 24 199761 48.433 15.824 0.069 9674977.691 36 210348 72.225 23.866 0.100 15192431.946 48 - - - - - 60 - - - - - 表 3 采用引导方式2时乘客分布不均衡状况
客流量/
(万人/日)计数/次 均衡度标准差 偏差 均值置信 总客流量/
人次12 40491 15.442 4.813 0.047 625260.697 24 134002 30.223 9.240 0.049 4049875.470 36 73084 42.361 12.963 0.094 3095923.827 48 111060 50.647 16.528 0.097 5624844.113 60 182196 64.214 21.855 0.100 11699623.509 表 4 采用引导方式3时乘客分布不均衡状况
客流量/
(万人/日)计数/次 均衡度标准差 偏差 均值置信 总客流量/
人次12 35141 11.519 4.409 0.046 404773.036 24 113375 22.488 8.567 0.050 2549572.184 36 68578 33.950 12.802 0.096 2328243.152 48 97391 43.842 15.812 0.099 4269783.179 60 180316 61.658 20.711 0.099 11117996.994 表 5 相对已有引导方式提升幅度
客流量/
(万人/日)引导方式1
均衡度标准差引导方式2
均衡度标准差引导方式3
均衡度标准差方式3相对方式1的提升幅度 方式3相对方式2的提升幅度 12 24.574 15.442 11.519 53.1% 25.4% 24 48.433 30.223 22.488 53.6% 25.6% 36 71.023 42.361 33.950 52.2% 19.9% 48 - 50.502 43.672 - 13.5% 60 - 64.135 62.323 - 2.8% -
[1] 广州地铁. 2017,2018,2019,2020,2021,2022年年报[EB/OL] 2023-06-28]. https://www.gzmtr.com/ygwm/gsgk/qynb/.
[2] 田洪金,殷 玲. 一种地铁站台门乘客引导系统的设计与实现[J]. 物联网技术,2022,12(5):75-76,79. [3] 王世轶,祁 勇,潘永剑,等. 一种地铁车厢拥挤度智能监测系统的研究[J]. 电子设计工程,2022,30(17):55-59. [4] 庞梦圆,杜彦君,厉超吉,等. 基于图像采集技术的地铁站台客流引导系统研究[J]. 都市快轨交通,2017,30(2):37-41. [5] Çapalar J, Nemec A, Zahradnik C, et al. Optimization of passenger distribution at metro stations through a guidance system[C]//Proceedings of the 16th International Conference on Computer Aided Systems Theory, 19-24 February, 2017, Las Palmas de Gran Canaria, Spain. Cham: Springer, 2018. 397-404.
[6] Ding H, Di Y R, Zheng X Y, et al. Passenger arrival distribution model and riding guidance on an urban rail transit platform[J]. Physica A:Statistical Mechanics and its Applications, 2021, 571(2): 125847. DOI: 10.1016/j.physa.2021.125847
[7] Zheng Z X, Wang H, Liu W M, et al. Toward real-time congestion measurement of passenger flow on platform screen doors based on surveillance videos analysis[J]. Physica A: Statistical Mechanics and its Applications, 2023, 612(4): 128474. DOI: 10.1016/j.physa.2023.128474
[8] 深圳地铁试点推出车厢拥挤度智能显示系统[J]. 都市快轨交通,2019,32(3):18. [9] 沙永杰,纪 雁,陈婉婷. 新加坡公共交通规划与管理[J]. 国际城市规划,2021,36(1):142-147. [10] 方晨晨,周继彪,董 升,等. 基于BP神经网络的地铁车厢拥挤度预测方法[J]. 交通信息与安全,2018,36(6):47-53. [11] 徐玉萍,江睿智. 交通运输效率、交通结构与经济发展[J]. 企业经济,2021,40(5):142-150. [12] 刘光通,孔德文,姚丽亚,等. 基于Wi-Fi嗅探技术的站台引导员效用及布置研究[J]. 交通信息与安全,2020,38(1):84-91. [13] 广州地铁官方网站. 广州地铁-列车介绍[EB/OL]. (2022-12-31)[2023-06-28]. https://cs.gzmtr.com/ckfw/lcjs_2020/. [14] 莫俊泓,石红国. 基于仿真的珠江新城站改造前后客流组织分析[J]. 综合运输,2021,43(12):119-126. -
期刊类型引用(3)
1. 刘文斌,马子彧,张春来,徐晓磊. 基于GIS技术的铁路网能力利用分析系统设计与应用. 铁路计算机应用. 2025(03): 38-43 . 本站查看
2. 黄方圆,徐泽儒,郭禅. 基于GIS技术的多模式航运运输里程精细化分析. 中国航务周刊. 2024(48): 59-61 . 百度学术
3. 杨文澈,刘茂朕,史晓磊,董卓皇,陈亚勋. 铁路货车里程管理及辅助分析应用系统设计及应用. 铁路计算机应用. 2024(12): 59-64 . 本站查看
其他类型引用(0)