带有最优流量信息中断期望影响的一类新格子模型.pdf

带有最优流量信息中断期望影响的一类新格子模型刘丽梅1，化存才2*(1.云南民族大学预科教育学院，云南昆明650504；2.云南师范大学数学学院，云南昆明650500)摘要：交通中断是导致交通扰动的一个重要因素.为研究最优流量信息中断的期望行为对交通流稳定性的影响，基于 Nagatanis 模型考虑最优流量中断概率的影响，提出一类新的格子流体力学模型.首先，根据线性稳定分析方法，得到新模型的线性稳定条件；其次，通过非线性稳定理论，获得该模型的 modifiedKorteweg-deVries(mKdV)方程，并根据求解方程得到的扭结反扭结孤立波解描述交通堵塞的演变过程；最后，通过仿真算例验证了理论分析结果.结果表明，考虑最优流量中断期望影响的新格子流体力学模型比 Nagatanis 模型更稳定，当添加最优流量信息中断的期望影响时，交通流的稳定性得到显著提高.关键词：交通流；格子模型；交通中断概率；mKdV 方程；数值模拟中图分类号：O175文献标志码：A文章编号：02587971(2024)01000807j+2随着经济的迅速增长以及社会现代化的不断发展，道路上的车辆数量逐年攀升，这也导致交通拥堵和交通事故日益频发.为缓解日趋复杂的交通问题，众多学者对交通流模型进行研究，交通流模型可分为宏观的交通流模型和微观的交通流模型1-7.1998 年，Nagatani 提出涉及交通流的宏观格子流体动力学模型8，这对研究交通流具有里程碑的意义.此后，学者们在此基础上对模型进行了改进，提出了许多考虑不同因素的格子流体力学模型，如李志鹏等9在 Nagatanis 模型的基础上考虑流量差的影响，并验证考虑流量差的影响能提高交通流的稳定性.田钧方等考虑密度差影响的格子流体力学模型，这对抑制交通拥堵起到积极的作用10.在现实生活中，驾驶员的个体特征也是影响交通流稳定性的重要因素之一，彭光含等在Nagatanis 模型的基础上考虑驾驶员预估效应的影响11.Kang 等12提出考虑驾驶员反应延迟影响的格子流体力学模型.之后，考虑驾驶员急性、胆量等特征的交通流模型也被学者研究.随着道路交通量的逐渐增多，道路交通环境也时常影响交通流的稳定性.特别的，由行人、事故等因素造成的交通中断是导致交通拥堵的一个重要原因.任何行驶的车辆均有可能出现中断，而交通中断往往是以不确定的概率发生.唐铁桥等13，提出考虑前车速度中断影响的车辆跟驰模型.之后，孙棣华等14在变道的 2 车道格子流体力学模型中加入中断概率，发现变道可以平均各车道的交通量，但当前格子流量中断会使得交通状况恶化.Jiang等15考虑后视以及中断共同影响的格子流体力学模型，结果表明 2 个因素都能提高交通流的稳定性.彭光含等16-18考虑格子最优流量中断的格子流体力学模型.j+1j+1j+1这些考虑中断影响的模型均能改善交通流的稳定性.在实际交通流中，交通中断概率是驾驶员对可能发生交通中断的预期行为.格子最优流量可看成是驾驶员在行驶过程中根据运行状况预估得到的优化流量，预估中断发生时最优流量变为 0，而以往的模型并未考虑格子最优流量中断的期望行为.为此，本文在 Nagatanis 格子流体力学模型的基础上，考虑格子最优流量信息中断的期望行为，提出一类新的格子流体力学模型；利用线性稳定和非线性稳定分析方法，得到模型的稳定性条件和对应的 mKdV 方程；收稿日期：2022-12-11；接受日期：2023-06-11；网络出版日期：2023-09-11基金项目：国家自然科学基金（12161095）.作者简介：刘丽梅（1994），女，云南人，硕士，研究实习员，主要研究非线性微分方程应用和计算.E-mail：.*通信作者：化存才（1964），男，云南人，博士，教授，硕士生导师，主要研究非线性微分方程应用和计算.E-mail：.云南大学学报（自然科学版），2024,46（1）:814JournalofYunnanUniversity:NaturalSciencesEditionDOI:10.7540/j.ynu.20220651通过求解 mKdV 方程，分析交通密度波传播演变的过程；最后，通过数值仿真验证理论分析的正确性.1模型的建立在 1998 年，Nagatani8提出最早的格子流体力学模型，其模型方程如下：tj+0(jvjj1vj1)=0,(1)t(jvj)=a0V(j+1)ajvj.(2)0jvj0V(j+1)V(j+1)式中：为平均密度，分别表示在第 j 格点 t 时刻的瞬时密度和瞬时速度，a 为驾驶员的敏感系数.表示最优流量，为驾驶员根据交通状况预估得到的最优速度，其形式为：V()=vmax2tanh20201c+tanh1c,(3)c0vmax其中：表示安全密度，表示平均密度，表示最大行驶速度.j+1在实际的交通中，交通中断概率是驾驶员对可能发生交通中断的预估行为，为揭示格子最优流量信息中断的期望行为对交通流稳定性的影响，本文在(2)式中添加最优流量中断的概率，提出新的格子流体力学演化方程，模型方程如下：t(jvj)=(1 pj+1)a0V(j+1)ajvj+pj+1(ajvj),(4)pj+1j+1pj+1式中：表示优化流量在第格子中断的概率，为简化计算将中断概率看成常数 p.因此，方程(4)可以写成如下形式：j(t+2)j(t+)+20(1 p)V(j+1)V(j)=0.(5)2线性稳定分析0V(0)对提出的格子流体力学模型(5)进行线性稳定性分析，研究优化流量中断概率的期望行为对交通拥堵演变的影响.在稳定状态下，假设交通流的稳定密度为，最优速度为，得到稳定状态的密度和速度分别为：j(t)=0,vj(t)=V(0).(6)yj(t)假定表示格点 j 上稳定状态密度的小扰动量，此时密度变为：j(t)=0+yj(t).(7)将(7)式代入(5)式得yj(t+2)yj(t+)+20(1 p)V(0)(yj+1yj)=0.(8)yj(t)=Aeikj+zt通过展开，方程(8)变为：e2zez+20(1 p)V(0)(eik1)=0.(9)z=z1ik+z2(ik)2+ik将代入方程(9)，可得到关于 的 1 阶、2 阶项的系数z1=20(1 p)V,z2=32z211220(1 p)V(0).(10)z2 0z2 0当,交通流是稳定的，但当时交通流将演化为非稳态流.由此，可得到交通流中性稳定曲线=1320(1 p)V(0),(11)均匀流的稳定条件可表示为：1320(1 p)V(0).(12)第46卷刘丽梅等：带有最优流量信息中断期望影响的一类新格子模型9p=0当时，上述稳定条件变成 Nagatanis 模型的稳定条件.0 0.1 0.2 0.3(c,ac)p=0根据(12)式可知，最优流量中断的概率 p 对交通流的稳定性具有重要的影响.由(11)式可画出扰动中断概率 p 在取不同值（，）时的中性稳定曲线，如图 1 中的实线.图 1 中虚线为共存曲线，由后续的非线性分析导出.每条中性稳定曲线都存在一个最高点，被称为临界点.在中性稳定线上方的区域交通流是稳定的，而在中性稳定线下方的区域交通流是非稳定的.当，即不存在扰动中断情况时，交通流的中性稳定曲线同 Nagatanis 模型，稳定区域的面积最小.之后，随着 p 的逐渐变大，交通流的临界稳定点和中性稳定线逐渐降低，交通流稳定区域的面积扩大，交通流越稳定.由此说明添加最优流量中断期望影响的格子模型比 Nagatanis 模型更稳定.同时增大最优流量扰动中断的概率 p 能提高交通流的稳定性.3非线性稳定分析(c,ac)为研究交通流密度波的传播演变行为，本文将采用还原扰动法在临界点附近导出 mKdV 方程分析交通演变的非线性现象.现对空间变量 x 和时间变量 t 引入缓慢变量 X 和 T，即X=(j+bt),T=3t,(13)0 1式中：，b 为待定参数.并且密度为：j(t)=c+R(X,T).(14)5将方程(13)(14)代入到方程(5)中，并对方程(5)进行 Taylor 展开至 项，得到如下方程2h1XR+3h22XR+4(TR+h33XR+h4XR3)+5(h5TXR+h64XR+h72XR3)=0,(15)h1=b+2c(1 p)Vh2=32b2+122c(1 p)Vh3=76b32+162c(1 p)Vh4=16V2c(1 p)h5=3bh6=58b43+1242c(1 p)Vh7=112V2C(1 p)V=dV(j)dj?j=cV=d3V(j)dj3?j=c式 中：，.此时，.=(1+2)cb=2c(1 p)V23在临界点附近令，消去和 项，方程(15)变为：4(TRk13XR+k2XR3)+5(k32XR+k44XR+k52XR3)=0.(16)k1=1272c(1 p)Vk2=16V2c(1 p)k3=122c(1 p)Vk4=1542c(1 p)Vk5=112V2c(1 p)式中：，.现进行以下变换T=k1T,R=k1k2R.(17)将方程(16)变成标准的 mKdV 方程TR3XR+XR3+MR=0,(18)M(R)=k3k12XR+k4k14XR+k5k22XR3O()式中：，忽视校验项，方程(18)变成标准的 mKdV 方程，方程的扭结解为：图1不同扰动中断概率p下密度与灵敏度的相图Fig.1The phase diagram of the density and the sensitivityunderdifferentinterruptionprobabilitiesp10云南大学学报（自然科学版）http:/第46卷R0(X,T)=CtanhC2(XCk1T).(19)式中：C 表示扭结反扭结解的传播速度，满足以下可解性条件(R0,MR0)=w+dXR0(X,T)MR0(X,T)=0.(20)MR0=MR，扭结反扭结解的传播速度 C 为：C=5k2k32k2k43k1k5.(21)由此，mKdV 方程关于密度的扭结反扭结孤子解为：j=c+k1Ck2c1tanhC2c1j+1Ck1c1t.(22)j=cAA=k1Ck2(c1)扭结反扭结孤立波解代表低密度时的自由相和高密度时的拥堵相，共存曲线是由推导出来的，其中.从图 1 中可以发现，中性稳定曲线(实线)和共存曲线(虚线)将交通流分成稳定、亚稳定、不稳定 3 个区域.随着参数 p 的逐渐变大，交通流的临界点逐渐变小，共存曲线逐渐降低，交通流不稳定区域的面积逐渐变小，交通流越稳定.即考虑最优流量中断的期望行为对改善交通流的稳定性起着积极的作用.4仿真模拟数值模拟用来验证理论分析的正确性，本节采用周期边界条件，初始条件选取如下18.j(1)=j(0)=00.1,j=N2，j(1)=j(0)=0+0.1,j=N2+1，j(1)=j(0)=0,j,N2,N2+1；(23)N=100 c=0=0.25 a=2.3 vmax=2式中：，.104p=0、0.1、0.2、0.3p=0p=0.1、0.2、0.3p=0.3图 2 为时间步长之后，扰动中断概率 p 分别取不同值下的密度波随时间的演变图，图 2(a)(d)分别对应不同的中断概率.从图 2 中可发现当时模型变成 Nagatanis 模型，此时，密度波的振幅波动较大交通流极不稳定.但从图 2(b)(d)，随着 p 逐渐增大，交通流的密度波的振幅逐渐减小，交通流越来越稳定.并且时的密度波的振幅均低于 Nagatanis 模型的振幅，即考虑最优流量信息中断的期望影响的格子模型比 Nagatanis 模型更稳定.当，满足稳定条件(12)，即使在初始条件下加入小扰动，但随着时间的推移，交通密度波的振幅最终消失，交通流趋于稳定.因此，考虑最优流量中断的概率能提高交通流的稳定性.0、0.1、0.2、0.3t=10 300p=0.1,0.2,0.3p=0p=0.3p=0p=0.3为更直观地研究最优流量中断的概率 p 对交通流密度波的影响，图 3 画出了 p 分别取不同参数下：(a)在时间步长的瞬时密度分布情况和(b)迟滞回环曲线.图 3(a)更能直观地显示，随着最优流量中断的概率 p 变大，交通流密度波的振幅逐渐变小.同时曲线的振幅均低于即 Nagatanis 模型的振幅.特别，当时，满足稳定条件(12)时，此时交通密度波的振幅保持稳定状态.图 3(b)绘制了交通密度与流量关系的迟滞回环图，当迟滞回环区域变小，流量波动的幅度小，交通流稳定.反之，迟滞回环区域大，则流量波动的幅度大，交通流不稳定.当时模型变为 Nagatanis 模型，迟滞回环的区域最大，此时交通流极其不稳定.但随着参数 p 的逐渐增大，迟滞回环的区域变小，交通流量波动的幅度越来越小，交通流越稳定.当时，迟滞回环曲线变成一个点，流量保持不变，交通流最稳定.因此，考虑最优流量中断的概率能提高交通流的稳定性.第46卷刘丽梅等：带有最优流量信息中断期望影响的一类新格子模型115结论交通中断可看成驾驶员对未发生的交通状况的预估行为，为分析最优流量中断的期望行为对交通流稳定性的影响，本文提出了一类扩展的 Nagatanis 模型.首先，通过线性稳定分析方法，得到关于最优流量中断概率 p 影响下的稳定条件，由相图得到，随着 p 的不断增大，交通稳定区域逐渐增大，并发现含有参数p 的格子模型稳定区域比 Nagatanis 模型大.其次，采用非线性分析方法，在临界稳定点附近导出 mKdV 方程，通过求解方程得到的扭结反扭结解描述交通密度波的演变过程，并画出共存曲线，进一步验证引入最优流量中断概率 p 对改善交通拥堵是有效的.最后，利用仿真算例验证了理论分析结果，进一步证明了含图2不同扰动中断概率p下密度波随时间的演化情况Fig.2Thespatiotemporalevolutionofdensitywavesfordifferentinterruptionprobabilityp图3瞬时密度分布和迟滞回环曲线图Fig.3Thedensityprofileofthedensitywaveandhysteresisloopcurvegraph12云南大学学报（自然科学版）http:/第46卷有参数 p 的格子流体力学模型比 Nagatanis 模型更能有效地抑制交通拥堵，即考虑最优流量中断的期望行为对提高交通流的稳定性具有积极的影响.本文探讨添加最优流量中断期望影响的格子模型对交通流稳定性的影响，后续可结合能耗和心理特征等方面，收集与实际相关的多交通信息研究多车道交通流，为仿真和控制提供更多的理论基础.参考文献：姜锐,吴清松,朱祚金.一种新的交通流动力学模型J.科学通报,2000,45(17):1895-1899.DOI:10.3321/j.issn:0023-074X.2000.17.019.JiangR,WuQS,ZhuZJ.AnewtrafficflowdynamicsmodelJ.ChineseScienceBulletin,2000,45(17):1895-1899.1SunYQ,GeHX,ChengRJ.Anextendedcar-followingmodelconsideringdriversmemoryandaveragespeedofpreced-ingvehicleswithcontrolstrategyJ.PhysicaA:StatisticalMechanicsandItsApplications,2019,521(C):752-761.DOI:10.1016/j.physa.2019.01.092.2翟聪,刘伟铭,黄玲.带有最优流量记忆时间差的一类新格子模型J.交通运输系统工程信息,2018,18(4):96-102.DOI:10.16097/ki.1009-6744.2018.04.015.ZhaiC,LiuWM,HuangL.AclassofnewlatticemodelwithconsiderationoftheoptimalcurrentchangeswithmemoryJ.JournalofTransportationSystemsEngineeringandInformationTechnology,2018,18(4):96-102.3李霞,汪一戈,崔洪军,等.智能网联环境下复杂异质交通流稳定性解析J.交通运输系统工程与信息,2020,20(6):114-120.DOI:10.16097/ki.1009-6744.2020.06.015.LiX,WangYG,CuiHJ,etal.StabilityanalysisofcomplexheterogeneoustrafficflowunderconnectedandautonomousenvironmentJ.JournalofTransportationSystemsEngineeringandInformationTechnology,2020,20(6):114-120.4TangTQ,HuangHJ,ZhaoSG,etal.AnextendedOVmodelwithconsiderationofdriversmemoryJ.InternationalJournalofModernPhysicsB,2009,23(5):743-752.DOI:10.1142/S0217979209051966.5袁娜,化存才.多前车速度差的车辆跟弛模型的稳定性与孤波J.物理学报,2012,61(16):16509.DOI:10.7498/aps.61.160509.YuanN,HuaCC.Analysisofthestabilityandsolitarywavesformultivelocitydifferencecar-followingmodeloftrafficflowJ.ActaPhysicaSinica,2012,61(16):16509.6沈奎,化存才.分数阶多速度差模型的孤立波J.云南大学学报(自然科学版),2019,41(S1):7-15.DOI:10.7540/j.ynu.2017.ShenK,HuaCC.SolitarywavesofafractionalmultispeeddifferentialmodelJ.JournalofYunnanUniversity(NaturalSciencesEdition),2019,41(S1):7-15.7NagataniT.ModifiedKdVequationforjammingtransitioninthecontinuummodelsoftrafficJ.PhysicaA:StatisticalMechanicsandItsApplications,1998,261(3/4):599-607.DOI:10.1016/S0378-4371(98)00347-1.8LiZP,LiXL,LiuFQ.StabilizationanalysisandmodifiedKdVequationoflatticemodelswithconsiderationofrelativecurrentJ.InternationalJournalofModernPhysicsC,2008,19(8):1163-1173.DOI:10.1142/S0129183108012868.9TianJF,YuanZZ,JiaB,etal.PhasetransitionsandtheKorteweg-deVriesequationinthedensitydifferencelatticehydro-dynamicmodeloftrafficflowJ.InternationalJournalofModernPhysicsC,2013,24(3):1350016-1350025.DOI:10.1142/S0129183113500162.10PengGH.Anewlatticemodeloftrafficflowwiththeconsiderationofindividualdifferenceofanticipationdrivingbehavi-orJ.PhysicaA:StatisticalMechanicsandItsApplications,2013,18(10):2801-2806.DOI:10.1016/sns.2013.03.007.11KangYR,SunDH.LatticehydrodynamictrafficflowmodelwithexplicitdriversphysicaldelayJ.NonlinearDynamics,2013,71(3):531-537.DOI:10.1007/s11071-012-0679-5.12Tang T Q,Huang H J,Huang S J,et al.A new car-following model with consideration of the traffic interruptionprobabilityJ.ChinesePhysicsB,2009,18(3):975-983.DOI:10.1088/1674-1056/18/3/022.13Sun D H,Zhang G,Liu W N,et al.Effect of explicit lane changing in traffic lattice hydrodynamic model withinterruptionJ.NonlinearDynamics,2016,86(1):269-282.DOI:10.1007/s11071-016-2888-9.14JiangCT,ChengRJ,GeHX.Animprovedlatticehydrodynamicmodelconsideringthe“backwardlooking”effectandthetrafficinterruptionprobabilityJ.NonlinearDynamics,2018,91(2):777-784.DOI:10.1007/s11071-017-3908-0.15PengGH,KuangH,QingL.Anewlatticemodeloftrafficflowconsideringdriversanticipationeffectofthetrafficinter-16第46卷刘丽梅等：带有最优流量信息中断期望影响的一类新格子模型13ruptionprobabilityJ.PhysicaA:StatisticalMechanicsandItsApplications,2018,507(1):374-380.DOI:10.1016/j.physa.2018.05.095.PengGH,LuWZ,HeHD.ImpactofthetrafficinterruptionprobabilityofoptimalcurrentontrafficcongestioninlatticemodelJ.PhysicaA:StatisticalMechanicsandItsApplications,2015,425(4):27-33.DOI:10.1016/j.physa.2015.01.045.17PengGH,JiaTT,KuangH,etal.EnergyconsumptioninanewlatticehydrodynamicmodelbasedonthedelayedeffectofcollaborativeinformationtransmissionunderV2XenvironmentJ.PhysicaA:StatisticalMechanicsandItsApplications,2022,585(C):126443(1-14).DOI:10.1016/j.physa.2021.126443.18AnewlatticehydrodynamicsmodelconsideringtheanticipationimpactofoptimalcurrentinterruptionLIULi-mei1，HUACun-cai2*(1.SchoolofPreparatoryEducation,YunnanMinzuUniversity,Kunming650504,Yunnan,China；2.SchoolofMathematics,YunnanNormalUniversity,Kunming650500,Yunnan,China)Abstract:Trafficinterruptionisanimportantfactorresultingintrafficcongestion.Tostudytheinfluenceofthe anticipation behavior of the optimal flow information interruption on traffic flow stability,a new latticehydrodynamics model is proposed considering the effect of optimal traffic interruption probability based onNagatanismodel.Firstly,thelinearstabilityconditionofthenewmodelisobtainedaccordingtothemethodoflinearstabilityanalysis.Secondly,themodifiedKorteweg-deVries(mKdV)equationofthemodelisderivedfromthenonlinearstabilitytheory,andtheevolutionoftrafficjamsisdescribedaccordingtothekinkantikinksolitonsolutionobtainedbysolvingtheequation.Finally,theresultsofthetheoreticalanalysisareverifiedbynumericalsimulations.TheresultsshowthatthenewlatticehydrodynamicsmodelconsideringtheanticipationeffectofoptimalcurrentinterruptionismorestablethanNagatanismodel.Thestabilityoftrafficflowcanbeeffectivelyimprovedwhentheanticipationeffectoftheoptimalcurrentinformationinterruptionisadded.Keywords:trafficflow；latticemodel；trafficinterruptionprobability；mKdVequation；numericalsimulation14云南大学学报（自然科学版）http:/第46卷