Methodologies of Pallet Allocation Over Pallet Pool in Intermodal Transportation Network

行为决策是智能网联汽车实现自动控制的重要组成^[1-6]，也是该领域的研究热点. 车辆通过信号交叉口时，其行为会受到信号灯和周围车辆的影响，如果不能准确预判通行条件，得到的决策结果会导致通行效率低、闯红灯甚至引发交通事故.

目前，智能网联汽车领域有关车辆通过信号交叉口的研究，主要关注简单交通环境下的辅助驾驶速度规划，研究目的大多为降低排放和燃油消耗^[7-12]，或者提高通行效率^[13-16]，涉及到的车辆行为是一些基本规则. 在行为决策层面，现有的智能网联汽车将一般道路场景的决策方法应用到信号交叉口，也有少部分学者专门研究在信号交叉口这种特殊复杂场景下的行为决策，只考虑信号灯影响时，Park等^[17]根据信号、车速、交叉口大小和车辆与停车线的距离构建了简单的停止/通过决策算法；Zhou等^[18]提出了基于强化学习的跟驰模型来获取合适的驾驶行为；Liu等^[19]基于行驶规则，在各种转向条件下构建了不同信号灯色时的决策流程图；苏冲等^[20]将动作序列划分为不同的行为状态，建立了各种工况下的状态迁移图. 考虑周围车辆影响时，学者们主要解决各方向车辆之间的冲突，在机器学习、运筹优化等方面进行探索，Gadepally等^[21]利用隐马尔可夫模型从连续的滤波中观测中估计驾驶员行为，为行为决策提供有效依据；Beaucorps等^[22]根据人类驾驶员的录音数据建立有其他车辆存在场景下的无信号交叉口和环形路决策算法；杜明博^[23]提出基于有限状态机的驾驶场景转换模型，构建基于ID3决策树的驾驶行为决策模型；de Campos等^[24]量化每辆车为避免潜在碰撞而配置的自由度，用成本函数和局部约束条件解决车辆之间的协调问题；Alonso等^[25]根据靠近交叉口车辆的状态信息构建车辆优先通行图结构，以此来解决不同方向车辆之间的冲突；宋威龙等^[26]设计了基于模糊逻辑的驾驶员类型判别模型，根据其他车辆的行驶意图和驾驶员激进程度建立了决策规则；Makarem等^[27]提出了一种分散式模型预测控制方法，对具有预定路径的车辆引入线性二次最优控制器，加入线性约束来避免碰撞.

从以上回顾可以发现，无论是否在车联网环境下，现有智能车辆在信号交叉口的控制研究主要集中在速度控制层面，在行为决策层面只有很少的研究，鲜有对通行决策策略的研究，考虑的行驶环境也非常单一，没有把周围车辆的干扰纳入考虑范围，特别是前车对智能网联汽车的行驶空间有显著影响^[28-29]. 如果不考虑前车状态，只考虑自身状态和信号灯，会使车辆在判断通过/不通过问题上有局限性，容易产生两难区. 鉴于此，本文关注的是车联网环境下智能网联汽车在信号交叉口的行为决策问题，为了让车辆准确控制运动，从微观交通流的本质出发，基于信号时间信息，利用Gipps跟驰理论预测前方车辆的状态和判断自身通过条件，同时加入换道行为来寻求更优的决策结果，从而构建一种智能网联汽车直行通过信号交叉口的决策模型.

1.   信号交叉口行为决策模型

1.1   信号交叉口行为决策流程

当前车路协同技术还未成熟，本文研究的交通环境是非协同控制的车联网环境，通过交叉口的智能网联主车辆SV （subject vehicle）只能通过车联网获取信号信息和其他车辆的位置、速度等，并根据信号灯色和时间、前车行驶意图做出相应的决策. 考虑一般情况，将信号灯情况分为绿灯（倒计时未出现）、黄灯、红灯、绿灯倒计时，只要车辆在可识别到信号灯范围内，就启动决策；当识别到的是红灯或黄灯，则做出停车的决策；当识别到的是绿灯（倒计时未出现），则继续按照跟驰状态行驶；当识别到的是绿灯倒计时，则识别自身在当前车道是领头车还是跟随车，并结合前车行驶意图判断能否在绿灯结束前通过停止线，若不能通过，则根据换道时间模型判断换道后能否通过，再做出相应的决策. 车辆的行为决策流程如图1所示.

图  1  信号交叉口行为决策流程

Figure  1.  Behavior decision process at signalized intersection

下载: 全尺寸图片 幻灯片

1.2   前方车辆状态预测模型

预测前方车辆未来的状态，本质就是判断车辆是否能够在绿灯倒计时间内通过停止线，这取决于车辆在倒计时间内走过的距离和倒计时间结束时的位置. 在交叉口，车辆为了追求高效行驶，都保持在跟驰状态，本文采用Gipps跟驰模型^[30]来表征前后车行驶状态的相互影响.

1.2.1   Gipps跟驰模型

Gipps跟驰模型属于安全距离模型，其基本思想是寻找一个特定的跟驰距离，当前车紧急刹车时，如果后车与前车的跟驰距离小于该特定距离，就有可能发生碰撞. 在交叉口，根据为避免与其前车（车辆n-1）发生追尾所需保证的安全距离，可以得到跟驰车（车辆n）行驶的安全速度，如式（1）.

式中：T为车辆n的反应时间；b_n为车辆n的最大制动减速度；s_n车辆n的效用尺寸，其值等于车身长度；v_n(t)为车辆n的速度；x_n(t)为车辆n沿行驶方向的位置；t为时间.

Gipps跟驰模型根据运动学推导，具有明确的物理意义，能较好地体现出行驶速度因素在跟驰过程中的影响，可以真实地模拟车辆之间的跟随行为，在交通仿真中有着广泛应用，如英国交通部的SISTM模型、美国的CARSIM模型等，是交通流理论的核心内容之一.

1.2.2   车辆状态预测

以反应时间作为车辆决策间隔，领头车是位置最靠近交叉口停止线的车辆，将领头车编号为P_V1，其行为不受其他车辆的影响. 通过Gipps跟驰模型，可以计算出跟驰车辆在反应时间后的速度，再结合初始的速度和位置可以得到车辆在反应时间后的位置. 求得领头车经过每个反应时间后的速度和位置后，则可以通过递推方法依次求得其后的所有跟随车辆在每个反应时间后的速度和位置.

按照与交叉口距离从小到大的顺序对所有车辆进行编号，依次为P_V1, P_V2, ···, P_Vn, 用v_P,1, v_P,2, ···, v_P,n分别表示其速度，x_P,1, x_P,2, ···, x_P,n分别表示其沿行驶方向的位置. 设所有车辆在初始时刻t₀的位置为x_P,1(t₀), x_P,2(t₀), ···, x_{P, n}(t₀)，速度为v_P,1(t₀), v_P,2(t₀), ···, v_P,n(t₀)，从初始时刻开始的第k （k = 1, 2, ···）个反应时间后的车辆状态如下：

1） 当n = 1时，

2） 当n > 1时，

1.3   不换道通过交叉口条件

车辆SV通过车联网获取到信号灯信息后，首先需要判断能否在当前车道通过交叉口，根据其前方是否有车辆，可将SV分为领头车和跟随车.

1.3.1   SV为领头车

当SV为领头车，其行驶行为不受其他车辆的影响，设SV在当前位置为x_P,S(t₀)，速度为v_P,S(t₀)，可以计算出SV先匀加速行驶，后匀速行驶到绿灯结束时刻的位置：

式中：等号右边第2项为车辆SV匀加速行驶的距离；第3项为按最大限速值匀速行驶的距离；t_g(t₀)为当前绿灯剩余时长；v_lim为道路最大限速值；a_max为车辆的最大舒适加速度.

设x_stop为停止线位置，可得到车辆SV在绿灯结束时刻与停止线的位置差为

1.3.2   SV为跟随车

当SV为跟随车，其行为决策结果受到前方车辆行为的影响，假设SV在未来按照Gipps跟驰模型行驶，那么可以用该模型描述其行驶状态并进行预测. 车辆SV每隔一个反应时间进行一次决策，可以得到车辆走完剩余绿灯时间需要进行h次决策，如式（6）.

根据式(2)和式(3)，可以得到当前车道中车辆SV之前的每辆车在绿灯信号结束时刻的位置为x_P,n(t₀ + hT)，车辆SV在绿灯信号结束时刻的位置为x_P,S(t₀ + hT)，如图2 （图中x、y分别为车辆的纵向、横向位置），根据车辆位置和交叉口停止线位置，可以判断每辆车是否能够在绿灯结束时通过停止线.

式中：Δx_P,n为第n辆车在绿灯结束时刻与停止线的位置差.

图  2  车辆位置

Figure  2.  Position of vehicles

下载: 全尺寸图片 幻灯片

SV在绿灯结束时刻与停止线的位置差$\Delta {x_{{\rm{P,S}}}}$也可由式（7）计算得出. 根据SV为领头车和跟随车的两种情况分别判断，若$\Delta {x_{{\rm{P,S}}}} > 0$，则SV在绿灯结束时刻不能通过停止线；如果$\Delta {x_{{\rm{P,S}}}} \leqslant 0$，则SV在绿灯结束时刻能够通过停止线.

1.4   换道通过交叉口条件

1.4.1   SV换道时间模型

当SV不能在当前车道通过时，可以寻求相邻车道的换道通行方案，在对换道策略进行评估时，最主要的考虑因素就是换道带来的时间损失，故以下分析换道过程的时间计算方法，SV换道的运动过程如图3所示.

图  3  换道过程示意

Figure  3.  Schematic of lane changing process

下载: 全尺寸图片 幻灯片

参考Yang等^[31]的研究，将车辆的换道轨迹描述为多项式曲线，假设车辆换道开始时的横纵向位置均为0、导航角处于水平状态、车辆在换道结束时的纵横向位置分别为x_f和y_f，以此来计算SV换道执行时间. SV换道执行时间为换道行驶的路程与平均速度的比值，如式（8）.

式中：

式中：${\bar v}$为SV在换道过程中的平均速度，取值为初始速度v₀与目标速度v_f的平均值；v_f为SV在换道结束时的目标速度，取值为目标车道前面车辆的初始速度，若目标车道无前车，则取值为初始速度提升15%后的值；l为换道行驶的路程；a_max,S为车辆在所有不同换道轨迹下产生的最大法向加速度；x_max,f为车辆在所有换道轨迹下完成换道需要的最大纵向位移；v_L,i为目标车道前面的第i辆车在当前的速度；m为目标车道的前车数量；ω为协调经济性和舒适性的权重.

1.4.2   目标车道前车通过条件判断

若目标车道的前车数量太多，会导致：1） SV没有换道空间，即车辆不能从当前位置开始换道；2） SV换道后通过停止线的概率减小. 所以，本文认为如果目标车道上远离交叉口的最后一辆车位置在与SV平行的道路横向截面之后，则目标车道前车不能通过停止线；如果该车在与SV平行的道路横向截面之前，则目标车道有前车；如果目标车道上与SV平行的的道路横向截面之前没有车辆，则目标车道无前车.

在换道过程中，目标车道前车保持为跟驰行驶状态，对其进行编号，依次为L_V1, L_V2, ···, L_Vm. 与前述计算当前车道的车辆状态相同的思路，可以得到目标车道前面的车辆在绿灯结束时刻与停止线的位置差为

式中：x_L,m(t₀ + hT)为目标车道前面的车辆L_Vm在绿灯结束时刻的位置.

1.4.3   SV换道后通过条件判断

有了换道时间后，可以计算出换道过程经历的反应时间个数为h₁=t_S/T，从而得到SV在换道结束时刻的速度和位置如式（14）.

SV换道结束后，后续阶段的通过条件判断和前述在当前车道的方法类似，从SV换道结束后的状态开始计算，如果目标车道没有前车，SV可以加速行驶，可以根据式（4）计算其在绿灯结束时刻的位置；如果目标车道有前车，SV只能跟随前车行驶，可以根据式（3）计算其在绿灯结束时刻的位置，最后得到SV在绿灯结束时刻与停止线的位置差Δx_L,S. 如果$\Delta {x_{{\rm{L,S}}}} > 0$，则SV换道后在绿灯结束时刻不能通过停止线；如果$\Delta {x_{{\rm{L,S}}}} \leqslant 0$，则SV换道后在绿灯结束时刻可以通过停止线.

2.   仿真分析

2.1   仿真方案

给定所有车辆当前的状态数据（车辆位置、速度）、车辆固定属性（最大舒适加速度、最大制动减速度）、道路属性（最大限速值、车道宽度、交叉口停止线位置、绿灯倒计时间），可以模拟出SV在剩余绿灯时间内的行驶过程. 本文考虑一般情况，场景设置为SV前方有车辆，相邻车道前方有车辆，通过配置车辆位置和当前绿灯剩余时间参数来构建不同的交通环境，以此来评估模型的决策能力，分析车辆的行为决策过程和速度变化.

2.2   模型评价

搭建不同的场景来评价模型的决策能力，设置两个不同绿灯倒计时间10.0 s和14.0 s的场景，分别为场景1和场景2，在两个场景中道路限速值为60 km/h，停车线位置为300 m，最大舒适加速度为2 m/s²，最大制动减速度为3 m/s². 在其他设置都相同的条件下，使用基于规则的现有决策模型和本文决策模型分别仿真车辆的运行，得到车辆行驶过程中的决策结果，两个场景如图4、图5所示.

在场景1下，由对比可知：在车辆靠近信号交叉口的前阶段，现有模型和本文模型都做出减速停车的决策，说明本文模型具有和现有模型一致的基本决策能力，在后阶段越来越靠近停止线时，在现有模型做出减速停车决策的情况下，本文模型做出直行跟驰的决策，即前者决策出不能通过交叉口，而后者预测出可以通过，说明本文模型能捕获到微小的交通状态变化，有使得决策结果更加精确的优势，从而提高车辆的通过效率. 在场景2中，增加了绿灯倒计时间长度，车辆使用本文模型的决策结果就一直为直行跟驰，这也符合绿灯时间越长，车辆通过交叉口概率越大的交通规律，说明本文模型能应用到实际的车辆运行控制中.

图  4  场景1决策过程对比

Figure  4.  Comparison of decision-making process for scenario 1

下载: 全尺寸图片 幻灯片

图  5  场景2决策过程对比

Figure  5.  Comparison of decision-making process for scenario 2

下载: 全尺寸图片 幻灯片

2.3   行为决策影响因素分析

2.3.1   当前绿灯剩余时间

将当前绿灯剩余时间作为变化量，分别以20.0，19.0，···，1.0 s的剩余时间来构建场景，场景的主要参数设置见表1.

表  1  不同绿灯剩余时间的场景输入

Table  1.  Scene input information for different green light countdown time

输入信息	数值		输入信息	数值
停止线位置/m	300		SV当前位置/m	171
道路最大限速值/（km•h−1）	60		SV当前速度/（km•h−1）	29
最大舒适加速度/（m•s−2）	2		当前车道头车位置/m	200
最大制动减速度/（m•s−2）	3		相邻车道头车位置/m	290

下载: 导出CSV 

| 显示表格

选取3个典型决策过程和速度变化结果如图6所示，当前绿灯剩余时间较短时，SV既不能在当前车道通过也不能换道，每次决策只能是停车；当前绿灯剩余时间在10.0 s和13.0 s之间时，SV在前阶段计算出不能在当前车道通过停止线，故做出了换道决策，换道之后的阶段做出了直行加速决策；当前绿灯剩余时间较长时，SV就一直做出直行跟驰决策. 综合以上分析，绿灯剩余时间每增加一个单位，SV在当前绿灯周期通过交叉口的概率可提升20%，车辆面对不同的当前绿灯剩余时间时做出的决策有较大差异.

图  6  不同绿灯剩余时间的决策过程和速度变化

Figure  6.  Decision-making process and speed change of different remaining time of green light

下载: 全尺寸图片 幻灯片

2.3.2   当前车道车辆位置

将当前车道车辆位置作为变化量，设置当前车道前车和车辆SV的不同位置来构建场景，场景的主要参数设置见表2.

选取3个典型决策过程和速度变化结果如图7所示，当前车道车辆与停止线的距离大于220 m时，SV不能在绿灯倒计时间内通过，故一直做出减速停车的决策；当头车位置在160～220 m时，SV在前阶段有换道的决策，换道之后的阶段做出了直行加速的决策；当前车道车辆与停止线的距离小于160 m时，SV可以在当前车道通过停止线，故一直做出直行跟驰的决策. 由此可知，不同的前车和SV自身位置得到的决策结果有较大差异.

表  2  不同当前车道车辆位置的场景输入

Table  2.  Scene input information for different vehicle positions in current lane

输入信息	数值		输入信息	数值
当前绿灯剩余时间/s	11		最大制动减速度/ （m•s−2）	3
停止线位置/m	300		SV当前速度/（km•h−1）	29
道路最大限速值/（km•h−1）	60		当前车道头车位置
最大舒适加速度/（m•s−2）	2		相邻车道头车位置/m	290

下载: 导出CSV 

| 显示表格

图  7  不同当前车道车辆位置的决策过程和速度变化

Figure  7.  Decision process and speed change of different vehicle positions in current lane

下载: 全尺寸图片 幻灯片

3.   结　论

1） 本文以周围车辆对智能网联汽车的影响为研究重点，以行驶时间和行驶路程作为通行条件计算依据，建立智能网联汽车在信号交叉口的行为决策方法，让其利用前车状态对执行各种行为的条件进行了超前判断，并得到准确的决策结果.

2） 仿真结果表明，决策模型可以使车辆根据车联网获取的信号倒计时间、车辆位置和速度状态，对交叉口停止线附近车辆的通过、停车和换道行为之间做出准确的选择，解决了前车不确定行为带来的决策困难问题，对智能网联汽车的实际控制具有指导意义.

3） 本研究考虑了车辆直行时的情况，当车辆有转向需求时，不同车道上车辆的状态和其他进口道车辆的冲突对决策也具有不可忽略的影响，未来研究关注的是车辆在交叉口有不同转向需求时的行为决策问题.