总第315期交 通 科 技
SerialNo.315
2
022第6期TransportationScience&Technology
No.6Dec.2022
DOI10.3963/j
.issn.1671 7570.2022.06.004 陕西省自然科学基金(
2022JM 307);陕西交通科技项目(22
67B);陕西高速机械化工程有限公司科研项目(KY01 22
)资助收稿日期:2022 09 14
基于时间序列分析的高速公路小修工程量预测
李许峰1 成高立1 罗雅琼2 史小丽2
(1.陕西高速机械化工程有限公司 西安 710038; 2.长安大学公路学院 西安 710064)摘 要 日常养护预算计划编制是高速公路日常养护管理的重要内容之一,小修工程量的大小决定了日常养护预算的高低,科学合理地确定小修工程量至关重要。文中通过定量化分析不同小修细目工程量与潜在影响因素之间的关系,
将高速公路小修细目分为有具体影响因素的细目和无具体影响因素的细目两类;
针对后者,采用时间序列分析法进行小修工程量预测,并依托某高速公路的3个小修细目进行验证。结果表明,时间序列分析法对于表现为随机性特征的小修细目工程量预测具有较好的测算效果。
关键词 高速公路 小修 工程量预测 时间序列分析 ARIMA模型中图分类号 U418
在公路养护市场化改革的过程中,
养护预算编制的准确性和及时性是科学制定养护投资计划的基础。对公路养护资金测算进行有针对性地定量化研究,对实现养护费用的科学化、规范化管理具有较大的理论和现实意义。目前高速公路日常养护费用估算主要是根据养护预算定额计算各单项资产养护活动的费用总和确定的,这种估价方法需要预先确定各养护活动的工程量,再分别乘以其单价,
汇总后得出总价。目前国内对高速公路养护工程量的预测主要是通过分析影响养护费用的因素,建立养护费用或工程量与影响因素的数学模型,影响因素的分
析过程常采用定性分析法[1 5]或定量分析法[6 7],
常用模型包括多元回归分析模型、BP神经网络、灰系统理论模型。部分学者着眼于探索路面使
用性能与公路养护造价之间的关联性[8
],文献[9]
用层次分析法对历年路面使用性能指数PQI重新标定,采用灰模型GM(1,1)预测短期PQI值,得出路面使用性能指数下降程度(100 PQI)与单位里程小修工程量呈线性关系。文献[10]建立路面状况指数与破损率之间的关系,并拟合工程量,对坑槽类破损效果较好。文献[11 12]采用灰理论GM(1,1
)模型预测高速公路养护成本;文献[13]使用时间序列分析的ARMA(1,1)模型对养管资金做出预测。
由于高速公路小修范围宽广、细目众多,现阶段研究更多聚焦于具有显著影响因素的小修细目(如路面小修细目),或者直接以日常养护费用这一总体性指标作为分析对象;
前者使小修预算在编制过程中出现部分细目工程量无据可依的状况,后者太过笼统,不够精细化。事实上,有些细目其小修工程量并不受具体因素的影响,而是在运营过程中自然衰减呈现出随机性的变化特征;本文以此类小修细目为研究对象,对这些随时间变化的小修工程量引用py
thon中的statsmodels库采用ARIMA模型进行时间序列分析,试图探寻此类小修细目工程量的变化规律。1 小修工程量变化特征分析
某高速公路双向路基长度256.618km,
路面面积3740827.23m2,桥梁双幅长度25355m,
桥涵长31457m,波形梁钢护栏总长536m。现有2004-2021年期间每年平均日交通量
(AADT)、标准轴载当量轴次(ESAL)、重车交通量(HT)、重车比例(HR)、特大桥比例、交通事故次数和各养护细目小修工程量数据。
对于高速公路小修工程量的逐年变化,首先会考虑是否与交通条件有关,在此借助SPSS软件运用双变量相关性分析方法研究小修工程量的影响因素。相关性分析要求数据服从正态分布,
Copyright ©博看网. All Rights Reserved.
对于不满足正态分布的数据取对数,取对数之后再进行相关性分析。以波形梁钢护拦板更换(二波)小修细目为例,对工程量与通车年限(OA)、AADT、ESAL、HT、HR进行皮尔逊相关性检验(皮尔逊系数为2个变量之间的协方差和标准差的商)。表1为相关性检验结果。
表1 相关性检验
ln工程量OAlnAADTlnESALlnHTlnHR
ln工程量1-0.219-0.235-0.26-0.276-0.107
OA10.646 0.615 0.605 -0.104
lnAADT10.971 0.939 -0.105
lnESAL10.992 0.131
lnHT10.241
lnHR1
注:“ ”表示2个变量相关系数检验的狋统计量显著性概率狆<0.01,相关性非常显著;“ ”表示狆
<0.05,相关性显著;“1”表示2个相同变量的相关系数为1,即完全相关;其余无标记数字表示狆>0.05,没有显著的相关关系。
由相关性检验表第一行可知,工程量与通车年限、年平均日交通量、标准轴载当量轴次、重车交通量、重车比例均无显著相关关系,即波形梁钢护拦板更换(二波)工程量的变化与项目记录交通条件无关。
归纳发现,标志标线的维修养护、伸缩缝维修、路肩与边坡的防护与清理、排水设施新建与养护、护栏零部件更换等小修细目并无具体的影响因素。对于这些细目养护工程量的预测,不能依赖因变量与自变量的关系建立回归模型,因此将目光放在时间序列分析上,研究工程量随时间的变化规律。
2 时间序列模型
时间序列分析通常用于对某一指标按一定顺序排列而成的数列进行研究分析,揭示数据本身蕴含的内在规律,进行短期预测。常用的时间序列模型有自回归模型(AR模型)、移动平均模型(MA模型)、自回归移动平均模型(ARMA模型)和自回归综合移动平均模型(ARIMA模型)。ARMA模型是前2种模型的组合,其三者原始时间序列均为平稳序列,ARIMA模型将非平稳序列转化为平稳序列后再进行剩余建模步骤。平稳序列即时间序列的均值和方差不发生明显变化。本文实例中无具体影响因素的小修细目工程量时间序列经平稳性检验后为非平稳序列,因此选用ARIMA(狆,犱,狇)模
型建立工程量与时间的关系。
ARIMA(狆,犱,狇)模型是一种时间序列建模方法,其建模的基本思想是对非平稳的时间序列用犱次差分(时间序列在狋时刻与(狋-1)时刻的差值)使其成为平稳序列,再用ARMA(狆,狇)模型对该平稳序列建模,之后经反变换得到原序列。式(1)为ARIME(狆,犱,狇)的一般表达式。
狔狋= 1狔狋-1+ 2狔狋-2+…+ 狆狔狋-狆+
狌狋-θ1狌狋-1-θ2狌狋-2-…-θ狇狌狋-狇(1)式中:
1
,
2
,…,
狆
为自回归系数;θ
1
,θ
2
,…,θ
狇
为移动平均系数;{狌}为误差白噪声,服从均值为0、方差为σ2的正态分布;{狔}为原始时间序列差分后的平稳序列。
3 分析步骤
3.1 数据处理
将历年小修工程量按等间隔的时间先后顺序组成时间序列,对于缺少数据的情况,采用插值法或借助移动趋势线补齐。当预测高速公路某年的小修工程量时,截取该年度以前的小修工程量作为时间序列。
3.2 平稳性检验
平稳性即要求经由样本时间序列得到的拟合
曲线在未来的一段时间内仍能顺着现有的形态“惯性”地延续下去,且序列的均值和方差不发生明显变化。对时间序列进行平稳性检验的目的是根据其平稳性建立与其相适应的时间序列模型。ADF检验即判断含有高阶序列自相关的序列是否存在单位根:如果序列平稳,就不存在单位根;
否则,就会存在单位根。所以,ADF检验的H
0假设是存在单位根,如果得到的显著性检验统计量TestStatistic小于置信度(10%,5%,1%)的临界值CriticalValues,则对应有(90%,95%,99%)的把握来拒绝原假设,即序列平稳。
3.3 数据的纯随机性检验(白噪声检验)
当经过检验确定是平稳序列之后的数据并非都能值得建立模型,只有符合动态规律性,即数据之间存在一定的相依性,历史数据对未来的发展有一定的影响才值得去建立模型并做出预测。纯随机性检验就是检验数据是否值得建立模型的一项重要步骤。在进行纯随机检验之前,先引入纯随机序列,纯随机序列是指序列值彼此之间没有任何相依性,过去的数据行为对未来的发展完全
7
1
2022年第6期李许峰等:基于时间序列分析的高速公路小修工程量预测
Copyright©博看网. All Rights Reserved.
没有影响。
从统计学角度来看,当时间序列为纯随机序列时,则表示该序列没有任何分析的价值。Py thon白噪声检验返回值lb_pvalue为基于卡方分布的狆统计量,当狆值小于0.05时表明序列为非白噪声序列。
3.4 模型定阶
狆、犱、狇是ARIME(狆,犱,狇)模型的3个关键参数;犱为差分次数,可由时间序列经平稳性检验和纯随机性检验后确定;狆、狇为阶数,即当前时间节点的值与前狆、狇个历史值相关,阶数的确定需根据自相关图(图1)和偏自相关图(图2)判断。
图1 自相关图
图2 偏自相关图
自相关函数(autocorrelationfunction,ACF)反映了狋时刻和(狋 犽)时刻小修工程量狔狋、狔狋-犽的
相关性;犽为滞后阶数,即当前时刻的值与前犽个时间点的值有关。自相关图显示沿狓轴的滞后阶数及狔轴上的自相关系数(-1~1之间)。自相关系数实际上反映的不仅是狔狋、狔狋-犽两者的相关性,还有狔狋-犽+1、狔狋-犽+2、…、狔狋-1对狔狋-犽的间接影响。
偏自相关函数(partialautocorrelationfunc tion,PACF)是剔除中间(犽-1)个变量的间接影响后,狔狋-犽与狔狋的直接相关性。长条状渲染为95%的置信区间。
通过函数的自相关图和偏自相关图,可以确定出该模型的阶数,然后进行参数估计及显著性检验。表2为ARIMA(狆,犱,狇)
模型阶数确定表。表2 犃犚犐犕犃(狆,犱,狇)
阶数确定表模型
自相关图
偏自相关图
AR(狆)
拖尾(衰减趋于0)狆阶后截尾
MA(狇)狇阶后截尾拖尾(衰减趋于0)ARMA(狆,狇)狇阶后衰减趋于0
狆阶后衰减趋于0
图1表示自相关图2阶截尾,
图2表示偏自相关图2阶拖尾。若原始序列是平稳序列,模型则为MA(2)
;若序列是经犱次差分后的平稳序列,模型则为ARIMA(犱,2)。4 工程实例
为了检验ARIMA模型预测方法的有效性,对依托工程高速公路2021年的小修工程量进行预测,并与实际值对比检验。以更换护栏螺丝、维修收水口(C25混凝土)、波形梁钢护拦板更换(二波)3个养护细目为例,表3为各年度小修细目工程量汇总表。
表3 小修工程量汇总表
年份更换
护栏
螺丝(
套)维修收水口
(C25混凝土)
/m
3波形梁钢
护拦板更换(二波)/m
200488.40200551.62
200621681.636452.95200764796.462328.882008
664818.341128.882009118164.151694.37201055422
1.272302.342011696118.881986.422012330356.962212.352013160887.25253.9620146448.753216.00201524114.346304.27201647712.154323.74201766446.3410880.47201880428.235360.002019100212.00242.792020112520.15730.002021
1330
16.80
3117.00
当预测2
021年的小修工程量时,取2020年及以前的数据作为时间序列。表4、表5为小修细目工程量时间序列平稳性检验和白噪声检验结果。
表4 时间序列平稳性检验结果表
细目名称
更换护栏螺丝(2阶)
维修收水口(C25混凝土)(2阶)
波形梁钢护拦板更换(二波)(3阶)
TestStatistic
-7.157
-3.482
-5.465
CriticalValues(1%)-4.01-4.22-4.33CriticalValues(5%)-3.10-3.19-3.23CriticalValues(10%)
-2.69
-2.73
-2.75
注:
TestStatistic为时间序列平稳性检验统计量,当小于CriticalValues(5%)的值代表置信水平为95%时的时间序列平稳。表5 时间序列白噪声检验结果表
细目名称更换护栏螺丝(2阶)
维修收水口(C25混凝土)(2阶)
波形梁钢护栏板更换(二波)(3阶)
lb_p
value0.038978
0.045194
0.043934
注:lb_p
value为基于卡方分布的狆统计量,当狆值小于0.05时表明序列为非白噪声序列。8
1李许峰等:基于时间序列分析的高速公路小修工程量预测
2022年第6期
Copyright ©博看网. All Rights Reserved.2022年高速免费日期
确定差分阶数犱后,根据犱阶差分序列的自相关图和偏自相关图,确定狆、狇值。表6为更换护栏螺丝、维修收水口(C25混凝土)、波形梁钢护栏板更换(二波)小修工程量模型系数表,表7为小修工程量预测值与实际值对比结果。
由于误差项是不确定的,所以没有具体的函数表达式,statsmodels库可直接得到预测值。
表6 模型系数表
细目模型变量coef犘>|狕|
更换护栏螺丝ARIMA(1,2,0)ar.L1
sigma2-0.5751
2.8466
0.014
0.037
维修收水口(C25混凝土)ARIMA(2,2,0)ar.L1
ar.L2
sigma2
-0.7329
-0.5373
1.0895
0.001
0.050
0.049
波形梁钢护栏板更换(二波)ARIMA(1,3,0)ar.L1
sigma2-0.6101
8.8466
0.033
0.031
注:ar.L1、ar.L2为回归模型的系数项(值为coef),sigma2为误差项犝t的方差。
表7 预测值与实际值对比表
养护细目
实际
工程量
预测
工程量
绝对值
误差/%
相对值
误差/%
更换护栏螺丝1330套1342.96套-12.96-0.97 维修收水口
(C25混凝土)16.8m315.00m31.8010.72
波形梁钢护栏
板更换(二波)3117m2627.45m489.5515.71 由表7可知,预测值与实际值的绝对差异较小,表明用时间序列模型进行预测的有效性;相对误差受到基数大小的影响,可作为辅助指标对预测结果进行判定。
5 结语
小修工程量的合理确定对日常养护预算编制具有十分重要的作用。通过对某高速公路的小修工程量数据分析,得出以下结论。
1)需采用定量化分析方法研究确定各小修细目是否具有影响因素。
2)当小修细目不具有显著的影响因素时,采用时间序列分析法预测小修工程量是一种可靠的方法。
3)不同养护细目需分别建立时间序列模型,且其模型不是一成不变的;因此需要及时更新积累的年度数据,用新数据加入后反映的趋势来预测未来年度的变化情况。
4)养护过程中数据的积累非常关键,从数据中挖掘出信息并用于决策,才能体现用数据进行养护管理的理念。
参考文献
[1] 付浩然,姚玉玲.高速公路养护工程量预测研究[J].北方交通,2016(6):79 81,85.[2] 高洪如.高速公路养护造价影响因素分析与费用控制研究[J].路基工程,2018,(增刊1):29 32.[3] 牛永亮,邱兆文.高速公路路面小修工程量预测模型研究[J].筑路机械与施工机械化,2008(1):54 56.
[4] 于英,王麒麟,田晋跃.高速公路养护成本预测模型[J].交通运输工程学报,2007(3):76 79,85.
[5] 滕伟玲,姚玉玲.高速公路小修保养工程量的预测模型[J].长安大学学报(自然科学版),2012,32(6):23 27.
[6] 张平,吴雪莲,高楠,等.高速公路交通安全设施小修保养费用测算[J].武汉理工大学学报(交通科学
与工程版),2021,45(1):117 121.
[7] 史小丽,李玉环,张明,等.高速公路路面小修保养费用影响因素量化分析[J].同济大学学报(自然科
学版),2019,47(9):1302 1309.
[8] 张玉峰,顾剑.公路养护工程造价影响因素评价研究[J].公路工程,2014,39(5):247 251.
[9] 吕朝阳,李梦晨,杨卿.基于路面使用性能指数的道路养护工程量预测[J].建设机械技术与管理,2015,28(4):61 65.
[10] 吴正莺.深圳市沥青路面日常养护工程量预测及养护决策分析[D].西安:长安大学,2012.
[11] 向红艳,徐莲怡.基于分数阶拓展算子GM(1,1)模型的高速公路养护成本预测[J].中外公路,2020,
40(1):278 282.
[12] 董香军,姚峰.灰理论GM(1.1)模型预测高速公路养护成本研究[J].山东交通科技,2014(5):18
20.
[13] 郜业强.基于ARMA模型的高速公路养护与管理资金需求预测[J].科学技术与工程,2011,11(7):
1636 1639,1648.
(下转第35页)
9
1
2022年第6期李许峰等:基于时间序列分析的高速公路小修工程量预测
Copyright©博看网. All Rights Reserved.
犃狆狆犾犻犮犪狋犻狅狀犚犲狊犲犪狉犮犺狅狀犉狌犾犾犉犾犲狓犻犫犾犲犎犻犵犺犕狅犱狌犾狌狊犃狊狆犺犪犾狋犘犪狏犲犿犲狀狋犛狋狉狌犮狋狌狉犲狌狀犱犲狉犎犻犵犺犜犲犿狆犲狉犪狋狌狉犲犪狀犱犎犲犪狏狔犔
狅犪犱犆狅狀犱犻狋犻狅狀狊犔犐犖犣犺犻狑犲犻1,犔犐犃犖犌犢狅狀犵犮狌狀1,犔犗犖犌犠犲犻狔狌1,犎犃犖犚狅狀犵犼
犻犲2
(1.No.1EngineeringCo.,Ltd.,FHECofCCCC,Beijing1
02205,China;2.SchoolofTransportationandLogisticsEngineering,WuhanUniversityofTechnology
,Wuhan430063,China)犃犫狊狋狉犪犮狋:BasedontheYaoundétoDoualaExpressway
inCameroon,consideredthecharacteristicsofhightemperatureandheavyloadtrafficunderthetropicalrainforestclimateinCameroon,themechan icalsensitivityofthefullflexibleasphaltpavementstructureisstudiedbyusingn
onlinearfiniteele mentmethod,soastooptimizethepavementstructuredesignscheme.Theresultsshowthatavoidingheavyloadandoverloadingisaneffectivemeasuretopreventdiseasesoffullflexiblehighmodulusas phaltpavement.Increasingthebasemoduluscaneffectivelyr
educethetensilestressunderthebaseandimprovetheanti fatigueperformanceoffullflexiblehighmodu
lusasphaltpavement,whilein creasingthesurfacemodulushasageneraleffectonimprovingtheperformanceoffullflexiblehighmodulusasphaltpavement.Byincreasingthethicknessofupperandlowerbasecourse,thehorizontalshearstress,verticalcompressivestressandbottomtensilestressofeachstructurelay
eroffullflexi blehighmodulusasphaltpavementcanbeeffectivelyreduced,andtheanti ruttingandanti fatigueperformancecanalsobeimproved.Accordingtotheresultsofmechanicalresp
onse,thefinalcon structionplanisthecombinationthicknessof6cmBBMupperlayer,9cmBBMlowerlayer,and10cmEMEbaseasphaltlay
er.犓犲狔狑狅狉犱狊:fullflexible;highmodulusasphaltpavement;mechanicalsensitivity;optimizationofp
avementstructure(上接第19页)
犘狉犲犱犻犮狋犻狅狀狅犳犕犻狀狅狉犚犲狆犪犻狉犙狌犪狀狋犻狋狔犅
犪狊犲犱狅狀犜犻犿犲犛犲狉犻犲狊犃狀犪犾狔狊犻狊犕犲狋犺狅犱犳狅狉犈狓狆狉犲狊狊狑犪狔
犔犐犡狌犳犲狀犵1,犆犎犈犖犌犌犪狅犾犻1,犔犝犗犢犪狇犻狅狀犵2,
犛犎犐犡犻犪狅犾犻2
(1.ShaanxiExpresswayMechanizationEngineeringC
o.,Ltd.,Xi′an710038,China;2.SchoolofHighway,Chang′anUniversity
,Xi′an710064,China)犃犫狊狋狉犪犮狋:Thedailymaintenancebudgetplanningisoneoftheimportantcontentsofthedailymainte nancemanagementofexpressways.Thesizeoftheminorrep
airworksdeterminesthelevelofthedai lymaintenancebudget,anditisveryimportanttodeterminetheamountoftheminorrepairworkssci entificallyandreasonably.Byquantitativelyanalyzingtherelationshipb
etweenthequantitiesofdiffer entminorrepairsandpotentialinfluencingfactors,theminorrepairsofexpresswayaredividedintotwocategories:thosewithspecificinfluencingfactorsandthosewithoutspecificinfluencingf
actors.Forthelatter,timeseriesanalysisisusedtopredictthequantityofminorrepairs,anditisverifiedbyrelyingonthreeminorrepairsofacertainexpressway.Theresultsshowthatthetimeseriesanalysismethodhasagoodcalculationeffectforthepredictionofthequantityofminorrepairswithrandomchar
acteristics.犓犲狔狑狅狉犱狊:expressway;minorrepair;predictionofrepairquantity;timeseriesanalysis;ARIMAmodel
5
32022年第6期
林志伟等:高温重载条件下全柔性高模量沥青路面结构的应用研究
Copyright ©博看网. All Rights Reserved.
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论