混合交通流参数之间关系模型标定
- 格式:pdf
- 大小:256.95 KB
- 文档页数:4
第32卷 第3期2008年6月武汉理工大学学报(交通科学与工程版)Jou rnal of W uhan U n iversity of T echno logy(T ran spo rtati on Science &Engineering )V o l .32 N o.3June 2008混合交通流参数之间关系模型标定3 收稿日期:2008201215 熊烈强:男,45岁,教授,主要研究领域为流体力学、智能运输系统 3中国博士后基金项目资助(批准号:20040350533);湖北省教育厅项目资助(批准号:2004D 008)熊烈强 刘海岷(武汉工业学院机械工程系 武汉 430023)摘要:提出了混合交通流参数之间关系模型的标定和超车换道率的计算方法,采用广深高速公路、广佛高速公路和沪宁高速公路的实测数据对模型进行了验证,并计算了各数据点的超车换道率.结果表明,流量、密度数据的相关系数都接近1,超车换道率介于-1和1之间.关键词:混合交通流;动力学模型;超车换道率;参数关系中图法分类号:U 491.20 引 言文献[1]提出的混合交通流流量、密度和超车换道率之间的关系模型为q =u f [m k +(1-m )k f e n]k f e n ≤k <k ″0<n ≤1u f [m k +(1-m )k f e n ]k f e n≤k <k ′e n-1≤n ≤014m u f k j e n k ′e n≤k ≤0.5k j e n-1≤n ≤0m u f k -k2k je -n0.5k j e n <k ≤k j e n-1≤n ≤0(1)式中:q 为干道流量;k 为密度;u f 为自由流速度,取道路的设计时速;k f 为以自由流速度行驶时的最大密度,由实测数据回归得到;m 为波速系数,0<m <1;k j 为速度为0时的阻塞密度,取111.1pcu(km ・ln )[2];n 为超车换道率,等于超车换道流量r 与干道流量q 的比值;k ′=k j4-(1-m )mk f .由式(1)得到的关系曲线如图1所示. 本文根据实测数据研究模型的标定及超车换道率的识别.图1 流量、密度和超车换道率之间关系u f 2自由流行驶;n =12超车换道率为1;n =02超车换道率为0;n =-12超车换道率为-11 数据来源数据来源于1998~1999年间在珠江三角洲地区进行的交通调查资料中广深高速公路、广佛高速公路的实测数据和国家“九五”攻关项目“公路通行能力研究”中沪宁高速公路1999年的实测数据.4条高速公路的设计时速均为120km h .表1为广深高速公路2车道的实测数据,表2为广佛高速公路2车道的实测数据;表3为沪宁高速公路2车道的实测数据.所有数据都是以实测5m in 流量(以5m in 作为一个时段)为基数乘以12,化为小时流量,并考虑了车辆折算系数的影响作出的.速度可由流量q 除以密度k 得到.表1 广深高速公路2车道数据序号q [pcu・(h・ln)-1]k [pcu・(km・ln)-1]n序号q[pcu・(h・ln)-1]k [pcu・(km・ln)-1]n1900.7999-0.78138019 4604.3889 0.458925 21101.0300-0.84960020 4804.6747 0.344875 31300.8735 0.31603321 5005.0015 0.129212 41501.3441-0.30734022 5305.3944-0.028080 51701.5946-0.42495023 5505.3858 0.434169 61901.7802-0.30685024 5705.7413 0.180211 72101.8996-0.01557025 5906.0149 0.056056 82302.1397-0.07347026 6106.0624 0.397194 92502.3665-0.09729027 6306.3425 0.286269 102702.6219-0.20660028 6506.7554-0.173660 112902.8003-0.09053029 7007.0958 0.305253 123102.9119 0.16711930 7407.6304 0.103546 133303.2220-0.05119031 8008.3065 0.052897 143503.4000 0.04994232 8308.9199-1.121250 153803.6658 0.18722633 8709.5562-1.074180 164003.9254 0.09655534 90010.2951-1.040280 174204.1416 0.10094235 95010.8559-0.986210 184404.3116 0.20496736104013.0293-0.895700 注:带下划线的数据是绝对值大于1,都是测量中变异的数据点.表2 广佛高速公路2车道数据序号q [pcu・(h・ln)-1]k [pcu・(km・ln)-1]n序号q[pcu・(h・ln)-1]k [pcu・(km・ln)-1]n1900.9168-0.00637022110013.8104-0.67047 2 1001.1148-0.33517023115016.1426-0.62602 3 1501.8061-0.63257024120016.0772-0.58346 4 2002.3384-0.00969025125016.4214-0.54263 5 2503.0255-0.19437026130016.9447-0.50341 6 3003.4932 0.43423827135018.3499-0.46567 7 3504.1716 0.34498728140021.3447-0.42931 8 4004.6431 0.75001529145021.7914-0.39421 9 4505.7070-0.85346030150021.9651-0.36031 10 5006.5020 #NUM!①31155021.8926-0.32752 11 5506.5947-0.72257832160023.4742-0.29577 12 6007.6258-0.69053033165023.8992-0.26500 13 6508.1057 0.15904934170024.4921-0.23515 14 7008.9755-1.322750②35175027.6899-0.20616 15 7509.2994 0.45010536180029.0980-0.17799 16 80010.0718 0.05491437185030.1794-0.15059 17 85010.9536-0.92830038190026.4919-0.12392 18 90012.2816-0.87114039195034.6236-0.09795 19 95012.2328-0.81707040200037.7786-0.07263 20100013.4300-0.76578041205034.4075-0.04794 21105014.2412-0.71699042212036.8696-0.01436 注:①#NUM!是计算结果为无穷大;②带下划线的数据是绝对值大于1,都是测量中变异的数据点.2 数据拟合由图1可知,在低密度k f e-1≤k≤14k je-1-1-mmk f e-1(2)范围内,流量2密度曲线是一组直线.因此,可用该密度范围内的数据,标定m和k f,步骤如下.1)初选k≤14k je-1范围内的数据,采用M athcad[3]中的线性拟合函数进行拟合,得到m和k f e n.・265・武汉理工大学学报(交通科学与工程版)2008年 第32卷表3 沪宁高速公路2车道数据序号q [pcu・(h・ln)-1]k [pcu・(km・ln)-1]n序号q[pcu・(h・ln)-1]k [pcu・(km・ln)-1]n11702.0736-0.769990163604.5837-1.495100 21902.0588 0.427753173704.5427-0.114740 32102.7190#NUM!①183804.6047 0.143584 42202.7747-1.423310②193904.6295 0.446815 52302.5413 0.509215204004.9219-0.082100 62402.7812 0.205424214105.0837-0.233200 72703.4148-1.353210②224205.1065 0.189608 82803.3593-0.017420234305.2327 0.198378 92903.3986 0.288706244405.6314-1.834110②103003.6554-0.183460254505.7166-1.104300②113103.6020 0.437257264605.8015-0.684570 123203.8548 0.060010274705.7488 0.196905 133303.8146 0.545849284806.0553-0.653730 143403.9645 0.496909295005.7754 0.967797 153504.2101 0.170001 注:①#NUM!是计算结果为无穷大;②带下划线的数据是绝对值大于1,都是测量中变异的数据点. 2)认为拟合曲线是n=0的曲线,则由k f e n可得到k f.3)将m和k f代入式(2),得到新的密度范围. 4)若所选用的数据均在新的密度范围内,则所得到的m和k f满足要求;否则,去掉超出新的密度范围的数据点,重新拟合,得到新的m和k f,返回步骤3. 按上述步骤标定的m和k f结果如表4所列.表4 4组数据标定m和k f回归参数值广深2车道广佛2车道沪宁2车道m0.76410.64520.6464 k f [pcu・(km・ln)-1]1.28530.44960.4655 采用M athcad中的corr(x,y)函数对数据分别进行速度与密度、流量与密度相关性分析,结果如表5所列.表5 数据相关系数广深2车道广佛2车道沪宁2车道corr(u,k)-0.963-0.949-0.512corr(q,k) 0.998 0.999 0.992 从表5中可知,广深高速公路和广佛高速公路的速度、密度数据的相关系数较高,速度、密度基本符合直线关系;沪宁高速公路的速度、密度数据的相关系数较低.因此,将Green sh ields的速度2密度线性模型[425]应用于沪宁高速公路时,误差较大.而4组流量、密度数据的相关系数都接近1,符合直线关系,符合式(1)的流量2密度关系.3 超车换道率识别按式(1)计算各数据点的n.当密度k<k f e-1时,交通流处于稀薄流,本文的理论模型没有涉及,不计算其n值.当密度处于k f e-1≤k≤14k je-1-1-mmk f e-1内时,由式(1)知,按下式计算n值n=lnq-m u f k(1-m)u f k f(3) 当密度处于14k je-1-1-mmk f e-1<k< 12k je-1范围内时,若q≥m2u f k jm k j-4(1-m)k fk,仍按式(3)计算n值;否则,按下式计算n值n=ln4qm u f k j(4) 当密度处于12k je-1≤k≤14k j-1-mmk f范围内时,若q≥m2uf k jm k j-4(1-m)k fk,则按式(3)计算n值;若12m u f k≤q<m2u f k jm k j-4(1-m)k fk,则按式(4)计算n值;否则按下式计算n值.n=-lnm u f k-qm u f k2k j(5) 当密度处于k>14k j-1-mmk f范围内时,若q≥12m u f k,则按式(4)计算n值;否则,按式(5)计算n值.n值的计算结果见表1~表3.4 超车换道率分析 以广佛高速公路2车道(超车道)为例,分析n 值的意义.・365・ 第3期熊烈强,等:混合交通流参数之间关系模型标定1)在低密度范围k f e -1≤k ≤14k j e -1-1-mmk f e-1(对应0.1654≤k ≤10.1279p cu(km ・ln ))内,n 值有“+”有“-”,表明既有车辆进入车道又有车辆驶离车道,并且车辆的进出是随机的.2)在密度范围14k j e -1-1-m mk f e -1<k <14k j -1-m mk f (对应10.1279<k <27.5305p cu(km ・ln ))内,理论上n 值应有“+”有“-”,但实测数据显示均为“-”,因为此时车间距不足100m ,车辆很难插入进来.密度较低时,n 值接近“-1”,表明大部分车辆回到了行车道;随着密度增加,n 值的绝对值越来越小,因为此时行车道密度也较高,回到行车道的车辆也越来越少.3)在密度范围k ≥14k j -1-mmk f (对应k ≥27.5305pcu (km ・ln ))内,n 值逐渐接近“0”,表明密度增加,超车换道越来越困难.4)结合图1,有2条直线值得注意.第1条是连接点(k =14k j e-1-1-mmk f e-1,q =14m u f k j e -1)和点(k =14k j -1-mmk f ,q =14m u f k j )的直线:q =m 2u f k j m k j -4(1-m )k fk ,该直线通过q 2k 坐标平面的原点,即在直线上的速度相等,;第2条是连接点(k =0.5k j e -1,q =0.25m u f k j e -1)和点(k =0.5k j ,q =0.25m u f k j )的直线:q =0.5m u f k ,该直线也通过q 2k 坐标平面的原点,在直线上的速度相等,u =0.5m u f .5)表1至3中绝大部分数据满足本文的理论模型.但标注#NUM !和带下划线的数据不满足.表明模型对变异的数据非常敏感,如广佛高速公路2车道的10号数据(q =500pcu(h ・ln ),k =6.502p cu(km ・ln ))是变异点,若将密度改为k =6.362p cu(km ・ln ),则n =-0.9466,为非变异点.从另一个角度看,按本文的模型处理数据可以判断数据的变异性.参考文献[1]熊烈强,王 富,李 杰,等.混合交通流参数之间的关系[J ].华中科技大学学报:自然科学版,2005,33(7):97299.[2]熊烈强,殷燕芳,王 富,等.匝道连接处交通流动力学的理论、模型及其应用[J ].武汉理工大学学报:交通科学与工程版,2006,30(2):2612264.[3]郑桂水.M athcad 2000实用教程[M ].北京:国防工业出版社,2000.[4]熊烈强,李 杰,严新平,等.适应IT S 的交通流参数之间的关系[J ].武汉理工大学学报:交通科学与工程版,2004,28(2):1712173.[5]A po sto lo sKo tsialo s ,M arko s Papageo rgi ou ,M o rgan M angeas ,et al.Coo rdinatedandintegrated contro l of mo to r w ay netw o rk s via non 2linearop ti m alcontro l [J ].T ranspo rtati onR esearch Part C ,2002,10:65284.Calib rati on of the R elati on sh i p Am ong Param eters of M ixed T raffic F lowX iong L ieq i ang L iu Ha i m i n(D ep a rt m en t of M echan ica l E ng ineering ,W uhan P oly techn ic U n iversity ,W uhan 430023)AbstractGiven are the m ethod calib rating the relati on sh i p am ong param eters of m ixed traffic flow and calcu lating the overtak ing rate .T he su rveyed data of Guangzhou 2Shenzhen ,Guangzhou Fo shan and Shanghai N an jing freew ay in Ch ina are adop ted to test the relati on sh i p s .T he resu lts indicate that the co rrelati on coefficien t betw een flow and den sity app roaches 1,and the overtak ing rate is betw een -1and 1.Key words :m ixed traffic flow ;dynam ical m odel ;overtak ing rate ;relati on sh i p of param eters・465・武汉理工大学学报(交通科学与工程版)2008年 第32卷。