1. 引言
行人过街是城市交通体系中重要的组成部分,在城市交通出行中占有很大比例。在大型信号交叉口行人二次过街路段,以步行交通量大,过街设施复杂,成为影响城市交通系统运行的一个重要因素。在城市平面交叉口的研究中,规模较大的混合交通下的交叉口内行人在一个有效绿灯时间内无法到达对面,容易产生安全隐患并且降低机动车流的运行效率[1]。通过信号控制行人二次过街,减少行人与机动车流的冲突,降低了事故隐患,但增加了行人过街时间。信号交叉口行人二次过街是城市复杂交通系统中的重要组成部分,为了探究行人二次过街的影响因素,提高行人二次过街时过街效率保障行人过街安全,在过去的研究中不同学者运用不同的模型与方法来探索研究影响行人过街的因素,分析行人的过街特征。
通过行人行为力学分析,构建行人动力模型,能够深入解析行人过街的交通特性,促进信号交叉口交通系统高效和安全的运行。许多学者在研究行人特征时提出了一系列行人流模型,包括气体动力学模型、元胞自动机模型和社会力模型等。其中,社会力模型受到了越来越多的关注,在微观行人流领域得到了广泛的应用。Helbing等人提出了社会力模型的概念[2],将人内部动机进行量化,使内心驱动转化为力,从力学方面考虑心理学问题,将力划分为三种力,自身驱动力、人与人之间的相互作用力以及障碍物对人产生的作用力。此后许多学者从压力的生理反应[3]、复杂环境中的应用[4]、心理学的原理解释[5]等多方面对社会力的理论进行了研究。
除此之外,许多学者还对社会力模型进行了改进,Sticco等运用真实视频数据对社会力模型的两个参数摩擦系数和身体力系数进行了优化[6]。Jiang等使用环境驱动的社会力模型来探索主体执行复杂的个体运动[7]。李娟等人根据实际行人过街情况,对行人的相对速度、期望速度和行人与障碍物之间的作用力进行了改进[8]。
在20多年的发展中,社会力模型不仅在疏散人群、消防救援等方面得到广泛应用,在交通领域也得到了部分应用,例如考虑行人过街时周边行人的干扰对模型进行了改进及仿真[9]。运用社会力模型对平衡车与行人混合交通的情况下搜集的数据进行校准[10]、对期望速度进行了分析和验证[11]等。但在行人过街方面,目前许多的研究主要是利用调查问卷进行对行人的调查与定性分析,例如运用调查问卷搜集在十字路口影响行人闯红灯过街的行为[12]。Avinash等人在印度混合交通的情况下通过现场调查及数据搜集分析得出年轻人的过街的步行速度最快并且表明城市特征也影响行人步行速度[13]。行人作为人行横道处的主要交通群体,且在过路过程时与其他群体有着交互行为[14]。因此,为了研究人行横道处多群体行为的演化规律,可以运用博弈论建立模型来分析行为的演化趋势以及影响因素。
综上所述,关于行人过街的行为特性分析,目前研究多集中在心理领域的定性描述,量化计算较少。为此,本研究聚焦信号交叉口行人二次过街行为,基于博弈论改进的社会力模型,解析二次过街的影响因素,为行人过街安全及提高过街效率提供更多的参考建议。
2. 改进后的社会力模型
行人二次过街主要运用于车流量较大车道过宽的十字路口,在道路中心位置设置安全岛,在运用原始社会力模拟行人二次过街后,认为在原模型基础上还要考虑行人二次过街的特性以及绿灯倒计时对行人的影响,所以本文对模型进行了改进,在利用视频数据分析研究行人二次过街的特性之后,发现绿灯倒计时对行人的步速以及运动轨迹都有较大的影响,所以在原始社会力模型的基础上增加了信号灯对行人的影响力。
改进后的社会力模型主要分为四部分:行人自驱力
;周边行人对被观察者产生的作用力
;斑马线边界线对行人产生排斥力
;绿灯时间对人产生的作用力
。如公式(1)所示。图1表明了各个分力的方向以及合力的方向。
Figure 1. Components of pedestrian social forces
图1. 行人社会力的各个分力
(1)
式中,
为行人i的社会力;
为表示行人i的自驱力;
为前方的行人j对行人i的吸引力和排斥力;
为边界线对行人i的排斥力;
为绿灯倒计时对行人i的吸引力;M为行人的质量;
为加速度矢量。
2.1. 自驱力
自身驱动力是由行人运动中受到目标吸引所产生的。在进行计算时,为了更符合行人过街特性增加了方向系数将不同方向的速度进行折算,如公式(2)以及(3)。
(2)
(3)
式中,
为t时刻时期望速度的方向向量;
为实际速度;
为期望速度;
为松弛时间;
为方向系数。
2.2. 行人的相互作用力
行人之间的相互作用力如公式(4)所示。其中sign函数用来表示行人的相互作用力是同向行人的吸引力还是对向行人的排斥力,如公式(5)所示。此外,由于行人之间的作用力受对方的速度和过街时间影响较大,为更符合行人二次过街A、B侧对时间认知不同(即A侧行人对时间的关注程度更高)的特征增加了时间系数。
(4)
式中,
为相用力强度系数;
为吸引力作用范围;
为行人视野的扇形半径;
为时间系数。
(5)
在行人行走过程中,过街时行人会形成椭圆力场与目标行人i有关,如公式(6)。
(6)
式中,
为椭圆的短半轴;
为目标行人i的实际行走方向;
为行人j的实际行走方向;
为行人j的当前实际速度;
为行人j的步幅宽度。
2.3. 边界线排斥力
在过街时当行人密度到达一定程度时,行人会超过斑马线边界向两边外扩,此时斑马线边界线对行人产生一股拉力且力场呈梯度
下降与边界线的力半径
及行人的作用力半径R有关,边界线形成的排斥力呈高斯分布,如公式(7)所示。
(7)
式中,
为边界线力的半径;R为行人受力范围半径。
2.4. 绿灯倒计时影响力
1) 绿灯倒计时影响力公式
在考虑行人二次过街的特性之后,增加了绿灯倒计时对行人过街的影响力,A侧行人受绿灯时长的影响大且倒计时越短对行人的影响越大,而绿灯倒计时
对于B侧行人的影响较,因此对于A、B侧行人进行了绿灯倒计时影响分级。并且绿灯倒计时对行人的影响主要以加速度的形式呈现。除此之外,增加了关于行人速度差的线性函数
以及关于时间的函数
,以及对于两侧过街行人的需求侧重不同,在关于时间和速度的函数运用完全信息动态博弈得出各自的权重系数a和b,且根据行人的特性及行走方向以及距离还设置了力的有效系数,进行区分不同类型的行人如公式(8)所示。
(8)
式中,
为绿灯倒计时;
为行人i的力场半径;
为绿灯作用距离;a,b为对时间和速度关注度概率;
为绿灯倒计时力场方向;
为行人i和j的相对速度。
2) 运用完全信息动态博弈设置的权重系数
在设置信号灯力的时,根据行人二次过街的特性,运用完全信息动态博弈,对时间和速度关注度概率设置权重,根据视频数据将A、B侧行人对时间和速度的关注度的不同设定速度矩阵中5个单位为零收益,时间矩阵中6个单位为零收益,根据两侧的零收益来计算出A、B侧行人各个博弈情况下的收益,如表1所示。
Table 1. Velocity matrix gains of pedestrians on both sides of A and B
表1. A、B两侧行人的速度矩阵收益
|
改变 |
不改变 |
改变 |
(−1, 0) |
(−2, 0.5) |
不改变 |
(0.5, −0.5) |
(1, 1) |
当A、B两侧行人行走方向都改变时,A侧行人追求尽快通过人行横道,设定A侧行人的距离值为6;B侧行人追求最短的通行距离,更快的通行速度,轨迹线进行小范围的改变,则B侧的距离值为5;所以该博弈的收益值为(−1, 0)。当A侧行人行走方向改变、B侧行人不变时,因B侧行人不改变行走轨迹,A侧行人要在原有轨迹上做出改变来避免与B侧行人发生冲突,因此取A侧的距离值为7,B侧不改变方向轨迹线对比都改变时的距离值要小,取B侧的距离值为4.5;所以该博弈的收益值为(−2, 0.5)。A侧行人不改变方向时,取距离值为4.5,B侧行人改变且相比于两侧都改变时B侧行人要改变的幅度要大,因此B侧的距离值为5.5;所以该博弈的收益值为(0.5, −0.5)。A、B两侧行人的行走方向都不改变时,两侧都是寻求最短轨迹因此两者的距离值都最小取为4,所以该博弈的收益值为(1, 1)。
由于A、B侧行人对时间的追求与对时间的认知都有所差异,并且在信号配时上,A侧绿灯时间短、过街距离长;B侧过街距离短、所以相比于A侧更追求过街舒适性,在对视频资料的行人过街的时间数据分析之后确定A、B侧行人时间收益矩阵,如表2所示。
Table 2. Velocity matrix gains of pedestrians on both sides of A and B
表2. A、B两侧行人的速度矩阵收益
|
改变 |
不改变 |
改变 |
(2, 1) |
(1.5, 0) |
不改变 |
(0, −1.5) |
(−1, −2.5) |
A、B两侧行人行走方向都改变时,A侧行人更追求过街速度,过街时间少,B侧行人追求舒适,轨迹改变较小且花费时间要比A侧多,所以将A侧时间值设为4,B侧时间值为5;A、B侧以支付6个单位的时间为零收益,所以该博弈的收益值为(2, 1)。A侧行人行走方向改变、B侧不变时,A侧行人为避让B侧行人在时间上会比两侧都改变要长,B侧通行时间长,但对时间的追求小于A侧,且B侧不改变方向时间花费更长,所以该博弈的收益值为(1.5, 0)。A侧行人行走方向不变、B侧改变时,A侧过街时间会更长,B侧因为时间充足但追求过街舒适度,过街时间花费更长,所以该博弈的收益值为(0, −1.5)。A、B两侧行人行走方向都不改变时,由于对时间认知的差异且双方都不让步,过街花费的时间会更久,所以该博弈的收益值为(−1, −2.5)。
计算出所有可能性的博弈收益之后,在时间和距离综合考虑下的收益找到纳什均衡,分别有8个不同的P值如表3所示,用来进行计算。
Table 3. Explanation of the variable P
表3. 变量P的解释说明
变量名称 |
变量解释 |
P1 |
当A侧行人和B侧行人都选择变方向时A侧行人的收益 |
P2 |
当A侧行人选择变方向,B侧行人不改变方向时A侧行人的收益 |
P3 |
当A侧行人选择不变方向而B选择变方向时A侧行人的收益 |
P4 |
当A、B侧行人都选择不变方向时A侧行人的收益 |
P5 |
当A和B都选择变方向时B的收益 |
P6 |
当A选择变方向而B选择不变方向时B的收益 |
P7 |
当A选择不变方向而B选择变方向时B的收益 |
P8 |
当A和B都选择不变方向时B的收益 |
同时,本文也引入权重系数b来当作B侧行人在考虑过街距离上的占比,这样就可以得到B侧行人在距离和时间同时考虑下的综合收益,用公式描述如(9)和(10)所示。
(9)
(10)
在行人二次过街时行人间的速度差和时间差都有较大关系,所以按照所测量的数据进行分析,得出了单调递增的速度差函数和单调递减的时间差函数。公式如(11)和(12)所示。
(11)
(12)
式中,
为行人i和j的速度差;
为行人i和j时间差;
为速度系数;
为时间系数。
行人所受的社会力会随着绿灯倒计时越来越短而逐渐增大,且绿灯倒计时只影响力的大小而不会改变力的方向。
3. 结果与讨论
3.1. 数据收集与处理
在模型改进前选取淄博市南京路与共青团西路交叉口的北进口东西向的行人二次过街区域为研究目标,且过街以年轻人为主,斑马线全长32 m、宽9 m,中间设有安全岛4 m,如图2所示,将二次过街分为两个阶段,这次主要研究AB段行人过街特性。
Figure 2. Schematic diagram of pedestrians crossing the street twice
图2. 行人二次过街示意图
在所要研究的交叉口行人二次过街中,利用无人机收集了10个工作日晚高峰时段的视频数据,共提取到276个信号周期的有效数据。运用Kinovea软件对目标行人进行了追踪,得到行人的速度、轨迹等。本文设置了4个横截面分别是1、2、3、4将路分成了五个路段,分别是路段S0、S1、S2、S3、S4,如图3所示。
Figure 3. The location of the cross section for data collection
图3. 数据收集横截面位置
其中在1横截面(大约4.5 m处)时,A侧后部分行人处于启动加速状态,B侧行人开始启动状态,如图4(a)所示。此时行人前方视野开阔且空间较为空旷,行人此时主要受自驱力和同向行人的吸引力。在S1至S2的路段,A侧行人慢慢加速,B侧行人开始启动并向A侧前进。
在2横截面(大约7.1 m处)在A侧与B侧“首位行人”交汇时的区段,如图4(b)所示,此时A侧行人速度开始下降,增加了对向行人对A侧行人的阻力。
在3横截面(大约10.2 m处)时,两侧行人完全进入了拥挤状态如图4(c)所示,A侧行人速度开始放缓,且增加了改变路线的次数。S3至S4的路段,拥挤的人群开始消散,行人的速度开始提升,随着绿灯倒计时的时间越来越小,行人的速度也越来越快。
在4横截面(13.5 m处)时,两侧行人拥挤状态即将消散。A侧行人马上到达安全岛,由于信号控制不均衡,导致A侧行人行走速度越来越大。
Figure 4. Pedestrian status on the sides of cross section A and B
图4. 横截面A、B侧行人状态
图5说明了信号交叉口A、B侧行人过街信号,从图中可以看出所调查的信号交叉口行人过街信号为两相位。在行人二次过街的第一阶段A侧为全红,B侧有效绿灯时间为51 s,在行人二次过街的第二阶段A、B侧行人有效绿灯时间均为33 s。
Figure 5. Effective green time for pedestrians crossing the street
图5. 行人过街的有效绿灯时间
在四个横截面我们用视频比对以及参考其他资料之后,结合路口特点,在每个横截面上分别搜集了行人的年龄、性别、步速、运动轨迹等多种数据。首先将收集到的数据按照横截面进行分类,在S1路段,主要受内心驱动的行人的实际速度大多在0.9~1.2 m/s。在S2路段由于感受到绿灯倒计时的压力及对向行人的阻力行人的实际速度大多在1.4 m/s~1.7 m/s之间,如图6所示。
(a) S1横截面 (b) S2横截面
(c) S3横截面 (d) S4横截面
Figure 6. Four cross-sectional pedestrian speed profiles
图6. 四个横截面行人速度分布图
拥挤状态下行人的横截面处的瞬时速度方向各异,处理计算时较为复杂,所以设置了6个速度方向以及8个区间
来进行数据的计算,如图7所示。由于y轴在实际场景中为南北方向所以不设置速度方向,当行人的行走方向在b区间时,都按2方向进行处理,有便于计算速度在x轴与y轴的分量。以此类推c区间速度方向都按照方向3进行计算,f区间的行人速度的方向都以速度方向5进行计算,g区间的速度方向都以速度方向6进行计算,以上每个区间的所占圆的1/8,且每个扇形的弧度均为45˚。
两区间的行人速度方向都归为速度方向4进行计算,
两区间按速度方向1计算。
Figure 7. Partition decision diagram for different speed directions
图7. 不同速度方向分区判定图
3.2. 结果讨论
3.2.1. 自驱力
根据视频中提取的各项数据以及对行人的调查与访谈,从A到B侧行人的平均期望速度
。行人的平均速度为
。计算后得到图8的曲线,行人在过街时的期望速度和自驱力与实际速度呈反比,即实际速度越大,自驱力就越小。从图中可以清楚的看到自驱力的变化,在开始过街时慢慢加速,此时内心对于过街目标明确,速度一直增长,从S1到S2时,当看到对向行人移动时,A侧行人的速度加快上升逼近于期望速度,所以自驱力开始减小。在2横截面处,由于与对向行人的距离拉近,斑马线上开始形成拥挤状态,A侧行人的实际速度有所下降,与期望速度的差值开始增大,自驱力也开始增大,并在3横截面的拥挤中心实际速度最小时,自驱力达到一个高峰值。在3横截面之后,拥挤的人群逐渐消散,实际速度增大,自驱力减小。
Figure 8. Trend graph of pedestrian self-drive force on side A
图8. A侧行人自驱力趋势图
3.2.2. 行人相互作用力
在开始过街时,行人的实际速度比期望速度要低,经过松弛时间达到期望速度。在S2段由于周围影响力增大,不同侧的行人受到的影响也不同,例如在行人二次过街时,影响行人过街的有同行的同伴,
Figure 9. Trend graph of pedestrian interaction force on side A
图9. A侧行人相互作用力趋势图
周围的同向行人以及对向行人对行人的影响力。
本文将行人过街时受到其他行人影响的作用力强度A分为三个等级,第一等级为距离特别近的同伴的此时将系数设为1、第二等级为对向行人的压力按照计算设为0.85,第三等级为同行的陌生人的压力设为0.5。以及在二次过街时AB侧行人对时间认知度不同,所以统计计算出A侧行人过街的平均过街时间为8~10 s,B侧为10~12 s。
在图9中我可以看出,在1横截面与向行人交叉之前,A侧行人属于起步状态,人与人之间的距离较小,状态较为稳定此时相互之间的作用力也较为平稳,在S1到S2时,人群逐渐散开人与人之间的距离拉大,行人与行人之间的相互作用力也开始减小,在S2至S3时,对向行人使A侧行人的行走空间缩小,并且对A侧行人有阻力所以合力进一步减小,从S3至S4对向行人减少,A侧行人聚集至安全岛,相互作用力逐渐上升。
3.2.3. 边界力
在行人过街时,斑马线附近并无较大障碍物,所以边界力对行人影响不大,根据计算行人流位置靠近斑马线边界线的行人流受边界力影响较大,而靠近中间的行人受到的边界力相对较小,所以根据计算采用定值作为行人二次过街中斑马线的边界力,如图10所示。
Figure 10. Trend map of boundary force on side A
图10. A侧边界力趋势图
Figure 11. Trend chart of the green countdown force on side A
图11. A侧绿灯倒计时作用力趋势图
3.2.4. 信号灯作用力
作为在行人二次过街中,导致行人速度波动变化最大的力是由绿灯倒计时所导致,比如在过街绿灯后期刚开始过街的行人可能会选择跑步前进。而在快要接近安全岛或已经在安全岛内的行人想要一次性通过该路口,尤其是在A侧的行人对时间的关注度较大,绿灯倒计时就成了人们后半段加速的主要原因。在A、B侧分阶段进行计算并且每个路段的绿灯倒计时力的影响系数的数值不同。A侧行人对于时间的看重且A侧行人在斑马线的前半段对于时间的关注度相比于后半段稍微低一些,在S3至S4时,绿灯倒计时越来越小且由于信号灯设置的不合理导致行人内心想要过街的心情变得更加焦虑所以力的走向呈上升趋势。如图11所示。
3.3. 模型验证
根据实际情况运用动能定理改进后的社会力模型与通过动能定理计算出的实际情况以及原始社会力模型进行量化之后得出如图12所示的趋势对比图,可以直观的看出改进后的社会力模型更符合行人二次过街的特性,在前半部分行人处于启动和加速状态,此时主要受自驱力的影响且行人的行走空间较大,所以原模型、改进后的模型与实际情况走向大致相同。随着行人密度的增大导致斑马线开始出现拥挤状况,所以后半部分开始人群到达拥挤状态,行人行走空间减小,到图中最低点时为过街人群最拥挤的状态,然后人群开始消散,并且A侧行人受绿灯倒计时的影响步速明显增大。所以图中可以看出,在后半段原模型的趋势与改进后模型、实际情况的走向相差较大,由于改进后的模型考虑了实际绿灯倒计时对行人的影响,所以与实际情况更为接近。但在图中最低点,人群到达最拥挤的状态,此时行人过街处于复杂交通状态,影响行人过街的因素增加,所以此点的原模型数据与改进后的模型数据都与实际情况相差较大。
Figure 12. Comparison of the original model, the improved model and the actual data
图12. 原模型、改进后模型与实际数据的对比图
4. 结论
本研究运用博弈论分析了实际情况下的行人二次过街两侧行人的过街需求,在A侧行人过街时对于时间的考虑比重较大且后期绿灯倒计时对其步速和情绪的影响较大,所以根据不同侧行人过街的不同需求,对原始社会力模型进行了改进,在考虑了饱和状态下两侧行人的过街需求以及信号配时不均衡的情况,得出了以下结论:
1) 对改进后的社会力的各个分力(自驱力、行人相互作用力以及绿灯倒计时力)进行了量化,将得出的改进后的模型与原模型和实际情况进行对比验证,结果表明改进后的模型能够更准确地反应行人过街特性。
2) 运用博弈论对流量不均衡的二次过街的行人对时间和速度的关注进行了量化,得出A侧行人对时间认知更强烈。
3) 改进后的社会力模型更能准确描述行人二次过街,前半段行走过程中原始社会力、改进后的社会力模型与实际情况下力的走向大致相同,在后半段行走路过程中考虑了绿灯倒计时影响的改进后的社会力模型更符合实际情况。
以上结论在未来为交叉口设置行人二次过街信号和改进行人过街设施设置提供参考。我们通过量化行人二次过街影响因素以及由于实际条件限制,只选取了一个实际路况进行了验证,模型还有很大的空间进行改进,未来我们将研究更多复杂交通状况下的行人过街,为了提高行人通行效率和通行安全,更适用于多种情况下的行人过街。