基于灰色ARIMA组合模型的砂石骨料物流需求预测——以浙江省为例

卢尔赛; 张改平; 赵良; 张朝晖

doi:10.16503/j.cnki.2095-9931.2022.05.010

交通运输研究 >

2022 , Vol. 8 >Issue 5: 84 - 90

DOI: https://doi.org/10.16503/j.cnki.2095-9931.2022.05.010

基于灰色ARIMA组合模型的砂石骨料物流需求预测——以浙江省为例

卢尔赛 ^,¹ ,
张改平 ¹ ,
赵良 ^,¹ ,
张朝晖 ²

展开

1.交通运输部科学研究院，北京 100029
2.日昌升集团有限公司，浙江杭州 310002

赵良（1992—），男，内蒙古兴安盟人，博士，助理研究员，研究方向为交通运输规划与管理。E-mail: 382100283@qq.com

卢尔赛（1988—），男，北京人，工程师，研究方向为物流管理与工程、交通运输规划与管理等。E-mail: 821105817@qq.com

收稿日期: 2022-04-19

网络出版日期: 2022-11-18

基金资助

工业和信息化部政策研究项目(YCLS-JCC-2002)

收起

Logistics Demand Prediction of Sand and Gravel Aggregate Based on Grey ARIMA Combination Model: A Case Study of Zhejiang Province

LU Er-sai ^,¹ ,
ZHANG Gai-ping ¹ ,
ZHAO Liang ^,¹ ,
ZHANG Zhao-hui ²

Expand

1. China Academy of Transportation Science, Beijing 100029, China
2. Richangsheng Group Co., Ltd., Hangzhou 310002, China

Received date: 2022-04-19

Online published: 2022-11-18

Fold

摘要

为提升砂石骨料物流需求量预测精度，建立灰色ARIMA组合预测模型，并用于预测浙江省内砂石骨料物流需求。首先，分别运用灰色ARIMA组合模型和灰色GM(1,1)模型两种预测方法对浙江省水泥产量进行预测并将二者的结果进行对比，验证灰色ARIMA组合模型预测精度。结果显示，灰色ARIMA组合模型的预测精度相较于仅应用灰色GM(1,1)模型有所提升。然后，运用灰色ARIMA组合模型对浙江省2021—2025年的砂石骨料物流需求量进行预测，结果显示2021—2025年浙江省砂石骨料物流需求量呈逐年增加的趋势。

关键词： 砂石骨料; 物流需求; 需求预测; 灰色GM(1,1)模型; ARIMA模型

本文引用格式

卢尔赛 , 张改平 , 赵良 , 张朝晖 . 基于灰色ARIMA组合模型的砂石骨料物流需求预测——以浙江省为例[J]. 交通运输研究, 2022 , 8(5) : 84 -90 . DOI: 10.16503/j.cnki.2095-9931.2022.05.010

Abstract

In order to improve the prediction accuracy of the logistics demand of sand and gravel aggregates, this paper established the grey ARIMA combination prediction model, and forecasted the logistics demand of sand and gravel aggregates in Zhejiang Province with this model. Firstly, the grey ARIMA combined model and grey GM(1,1) model were used to forecast the cement output in Zhejiang Province respectively, and the results were compared to verify the prediction accuracy of the grey ARIMA combined model. The results showed that the prediction accuracy of the grey ARIMA combined model was higher than that of the grey GM(1,1) model only. Secondly, the grey ARIMA combination model was used to forecast the logistics demand of sand and gravel aggregates in Zhejiang Province from 2021 to 2025. The results showed that the logistics demand of sand and gravel aggregates in Zhejiang Province would increase year by year from 2021 to 2025.

Key words： sand and gravel aggregate; logistics demand; demand forecast; grey GM(1,1) model; ARIMA model

0 引言

砂石骨料是一种在混凝土中起骨架、填充和稳定体积作用的岩石颗粒等粒状松散材料，广泛应用于房屋建设、交通基础设施建设、市政工程建设等领域，是我国基础设施建设用量较大且不可或缺的原材料。2021年，我国砂石骨料产量居大宗物资首位，为200亿t^[1]。为了更好地保障砂石骨料供应、节约砂石骨料资源，应对砂石骨料物流需求进行精准预测。目前国内砂石骨料市场需求预测主要依靠经验，误差较大，已不适应目前经济迅速发展的形势。因此，准确地预测砂石骨料物流需求，对于节约砂石骨料资源、保障砂石骨料供应具有一定的现实意义。

目前，物流需求量预测方法有很多。国外学者对物流需求预测研究起步较早，且集中于改进单一算法的缺陷、提高单一算法的精度方面。Bruzda^[2]基于大数据分析研究出一种新型农产品物流需求量预测系统，并设计了系统的软硬件部分；Baisariyev等^[3]运用Bootstrap方法对航空备件物流需求进行了预测；Erwin等^[4]运用改进灰色模型对人道主义物流需求进行预测，以进一步优化阿姆斯特丹行动中心的物流调配过程；Ryuichi等^[5]通过建立贸易物流预测模型，基于亚太经合组织的贸易额对该组织成员国的国际物流进行了预测。国内学者对物流需求预测的研究思路与国外学者的相似，也着重于研究单一算法精度的提高。陈长英^[6]以2008—2017年的广西壮族自治区货运量为基础，运用改进的灰色-马尔可夫链模型对广西壮族自治区货运周转量进行了预测，得出该自治区物流需求量将逐年递增的结论；王子健^[7]运用BP神经网络对关中平原城市群中不同城市物流需求进行了预测，并根据未来物流需求走向对城市群内部各个城市之间的物流竞争力进行了评价；李明书等^[8]基于时间序列法选择ARIMA模型对长春市邮政物流总量进行了研究与预测，结果表明ARIMA模型对呈周期性变化的时间序列预测较为准确，可以用来进行邮政物流总量预测；谭伟华^[9]融合多元回归分析法与神经网络分析法，对江西省物流需求进行了预测；王晓平等^[10]提出了基于支持向量机模型的北京城镇农产品冷链物流需求预测方法，并验证了新构建模型的预测结果具有较高的精度，可在一定程度上为相关决策提供依据；徐慧^[11]构建了基于q阶模糊的层次分析法，对医药应急物流供应能力水平进行评价及预测，并有效解决了已有的模糊层次分析法不足以匹配模糊分析水平的问题；刘庆庆等^[12]基于2010—2019年连云港物流量数据，运用灰色GM(1, 1)模型对连云港2020—2024年物流需求量进行了预测研究，并确定了GM(1, 1)模型的精度。

综合来看，以往的研究多基于原始数据直接建模，并往往运用单一的预测方式。如果单纯运用一种模型进行预测，对数据的规律性要求较高，但现实数据往往没有很强的规律性，导致预测结果误差较大。因此，为减小预测结果的误差，有必要运用多模型组合进行预测^{[5⇓⇓⇓-9]}。组合预测是指将广泛认可的预测模型，比如灰色GM(1,1)模型、多元回归模型、马尔科夫链模型、ARIMA模型、BP神经网络模型等组合运用，一般组合方法为先运用GM(1,1)模型或多元回归模型对原始数据进行处理，再用马尔科夫链、ARIMA模型等进行预测。组合预测在使用无序的时间序列数据前,先对数据进行回归处理，再进行预测，可以在一定程度上提高预测的准确性^[6]。

因此，本文运用灰色GM(1,1)模型与ARIMA模型对砂石骨料物流需求进行组合预测。首先运用灰色ARIMA模型和GM(1,1)模型对浙江省水泥产量进行预测，然后以真实产量为基准进行对比研究，验证灰色ARIMA模型的预测精度。接下来运用灰色ARIMA模型对浙江省2021—2025年的砂石骨料需求量进行预测。

1 模型构建

1.1 问题描述

砂石骨料物流需求的增长水平受经济环境、人口规模、建设政策等可知因素以及突发重大事件、突发政策、自然灾害等无法预测的未知因素的影响，且相互之间关系复杂、变化无序、相互影响，难以做到客观量化，是一个典型的灰色系统，故可用灰色GM(1, 1)模型进行预测。但灰色预测依然存在一定的问题，即对时间序列的平滑性要求较高，如果平滑性不够高，则会产生较大误差。为此，本文在对原始时间序列进行灰色处理后，运用ARIMA模型对时间序列数据进行预测，以降低数据平滑性不足带来的误差。

1.2 灰色ARIMA组合模型

灰色ARIMA 组合模型主要由灰色预测模型和 ARIMA 模型组合而成。首先建立灰色预测模型对产量数据序列进行拟合，然后建立ARIMA 模型对灰色模型的拟合误差序列进行预测^[13-14]，最后将两种模型的预测值求和构成产量预测值。

组合模型的预测步骤如下：

（1）设训练集X⁽⁰⁾={X⁽⁰⁾(1), X⁽⁰⁾(2),…, X⁽⁰⁾(n)}，并设预测数据长度为r的观测数据集为

X (0)' = {X (0) (n + 1), X (0) (n + 2), …, X (0) (n + r)}

。对训练集序列建立灰色GM(1, 1)预测模型，得拟合序列

X^(0) (n)

，并用残差检验法对灰色预测结果进行检验^[15⇓-17]。

求

X^(0) (n)

拟合序列的相对误差序列

E 1 (0) (n) = X (0) (n) - X^(0) (n) = {E 1 (0) (1), E 1 (0) (2), …, E 1 (0) (n)}

其中，

X (0) (n)

为初始序列；

X^(0) (n)

为预测序列。

（2）对相对误差序列进行修正。找出

E 1 (0) (n)

中最小负数

E [1, m i n] (0)

，然后令

E 2 (0) (n) = E 1 (0) (n) + E [1, m i n] (0)

，得到数列

E 2 (0) (n)

^[18-19]。

（3）判断误差序列

E 2 (0) (n)

的平稳性，并通过对序列

E 2 (0) (n)

做ADF检验，确定差分阶数

d

。

（4）对差分处理后的平稳序列建立ARIMA模型，运用自相关分析法对序列的自相关和偏相关函数图进行模型阶数初步识别，参照表1，经过AIC准则判定，针对ARIMA（p, d, q）（p为自回归项系数，q为滑动平均项数，d为使数列成为平稳序列所进行的差分次数）模型确定结果最优时对应的p, q，并用该模型得到误差预测序列

E^2 (1) (n)

。

表1 ARIMA模型的相关性特征

模型	自相关函数	偏相关函数
AR(p)	拖尾	p阶截尾
MA(q)	q阶截尾	拖尾
ARMA(p, q)	拖尾	拖尾
ARIMA(p, d, q)	拖尾	拖尾

（5）对

E^2 (1) (n)

进行差分还原，得序列

E^2 (0) (n)

，再把

E^2 (0) (n)

序列还原成原始序列的预测数列，还原公式为

E^1 (0) (n) = E^2 (0) (n) + E [1, m i n] (0)

。

（6）对ARIMA模型拟合结果进行白噪声检验，若未通过检验，说明还有一些重要的信息没有提取，则回到步骤（5）对拟合模型进行重新设定，直到通过白噪声检验。

（7）将灰色模型的预测结果

X^(0) (n)

与修正后的ARIMA模型预测结果

E^2 (0) (n)

相加，得出最后组合模型的数据

Y^(0) (k)

，即

Y^(0) (k) = X^(0) (k) + E^1 (0) (k)

^{[20⇓⇓-23]}。

灰色ARIMA模型的预测流程如图1所示。

显示原图|下载原图ZIP|生成PPT

图1 灰色ARIMA模型的预测流程

2 实证分析

2.1 浙江省砂石骨料物流需求现状

由于砂石骨料需求市场没有官方统计，市场不规范，且国家规定的开采量与实际开采量之间差异较大，故本文首先梳理砂石骨料供需情况，以得到砂石骨料需求总量的实际构成。2021年，中国200亿t的砂石骨料需求量总体自给自足，进出口总量在3 000万t以下^[24]，故在计算砂石骨料总需求时，不需要考虑进出口量。根据《2021年中国砂石行业运行报告》，砂石需求、运输和供给如图2所示^[24]。

显示原图|下载原图ZIP|生成PPT

图2 砂石骨料供需关系图

从图2中可以看出，砂石骨料的总需求规模D为预拌混凝土用砂石需求量、沥青混凝土用砂石需求量和水稳层用砂石需求量的总和，即：

（1）

D = D 1 + D 2 + D 3

式（1）中：D为砂石骨料总需求量；

D 1

为预拌混凝土用砂石需求量；

D 2

为沥青混凝土用砂石需求量；

D 3

为水稳层用砂石需求量。

预拌混凝土用砂石量

D 1

可通过各地水泥产量推测。水泥砂石配比往往按照水泥∶砂石骨料 = 1∶6来进行混合，即水泥用量是砂石骨料用量的1/6。同时，根据《2021年中国砂石行业运行报告》得出，

D 1

约占市场总需求量的70%^[24]。基于上述条件，可根据水泥产量近似推算砂石骨料需求量

D

：

（2）

D = D 水 泥 × 670

式（2）中：

D 水 泥

为水泥产量。

2.2 模型精度验证及求解

根据上述分析，本文选择1990—2020年浙江省水泥产量（资料来源于1990—2020年的《浙江省统计年鉴》^[25]）来预估浙江省2021—2025年砂石骨料需求，建立灰色ARIMA组合模型进行拟合和预测。

首先，基于1990—2015年的水泥产量数据（见图3），分别运用GM(1, 1)模型与灰色ARIMA模型预测2016—2020年水泥产量数据，并将预测结果与真实产量进行对比，以验证灰色ARIMA模型的预测精度。

显示原图|下载原图ZIP|生成PPT

图3 浙江省1991—2015年水泥产量

然后对浙江省2021—2025年的水泥产量进行预测，最后根据水泥产量预测值对该省相应年份的砂石骨料需求量进行预测，以得出浙江省2021—2025年的砂石骨料物流需求量。

由图3可看出，水泥生产量出现了几个特殊的节点：在1990—1998年，水泥生产量保持平缓上升趋势；在1999—2008年间水泥生产量快速上升；自2009年之后，产量在波动中呈缓慢上升趋势，且伴随下降趋势。

运用GM(1, 1)模型对1990—2015年的水泥产量数据进行拟合，分析并预测2016—2020年的水泥生产量。GM(1, 1)模型预测值如表2所示。

表2 GM(1, 1)模型预测值

预测阶数	预测值/万t
1	14 287.7
2	14 921.7
3	15 565.9
4	16 220.6
5	16 885.8

GM(1, 1)模型预测拟合结果如图4所示。

显示原图|下载原图ZIP|生成PPT

图4 浙江省水泥产量GM（1, 1）模型结果拟合图

从拟合曲线以及预测值上看，单一的GM(1,1)模型预测误差较大，需要应用更高精度的预测模型才能更好地对浙江省水泥产量进行预测。

对水泥生产量实际值和GM(1,1)模型拟合值求残差序列，再将非负处理后的灰色残差序列代入ARIMA 模型进行ADF平稳性检验，检验结果如表3所示。

表3 ADF平稳性检验结果

期数	Q统计量	Prob值
1	21.5	0.000
2	37.1	0.000
3	47.3	0.000
4	51.5	0.000
5	52.6	0.000
6	52.7	0.000
7	52.8	0.000
8	54.7	0.000
9	58.9	0.000
10	64.9	0.000
11	72.5	0.000
12	81.2	0.000

要判断序列是否平稳，需要看ADF检验结果中的Prob值，当大多数的Prob值都大于0.05（置信水平）时序列为白噪声序列，序列平稳；反之则为非平稳序列。由表3中Prob值可知序列是非平稳的，需要对序列进行差分处理。

对序列进行一阶差分处理，即d=1，结果如表4所示。由于此时Prob值全部大于0.05，所以序列是平稳的，停止差分处理，对其进行相关性检验，结果如表5所示。

表4 ADF单位根检验

单位根检验方法	置信水平	Prob值
ADF-Fisher检验	0.6	0.7
ADF-Choi检验	0.6	0.7

表5 残差序列的自相关和偏相关函数

期数	AIC(p)	AIC(q)	Q统计量	Prob值
1	142 837.3	143 179.3	1.3	0.3
2	133 743.3	134 433.4	1.9	0.4
3	127 447.6	130 945.9	2.9	0.4
4	129 907.7	128 834.1	2.9	0.6
5	128 597.4	129 338.4	4.5	0.5
6	129 332.1	131 345.9	4.7	0.6
7	130 857.2	132 945.6	5.1	0.7
8	131 581.3	133 432.1	5.3	0.7
9	133 512.1	133 956.2	7.1	0.6
10	131 651.6	134 512.8	7.8	0.6
11	132 621.8	133 612.7	8.9	0.6
12	133 638.9	135 671.1	13.9	0.3

由于此时Prob值全部大于0.05，所以序列是平稳的，此时停止差分处理，对其相关性进行检验。通过表5可初步判断，自相关函数拖尾，偏相关函数拖尾，并且其中p与q的AIC值在p=3, q=4时达到最小，即选定最优模型为ARIMA(3, 1, 4)。

模型阶数确定后，运用ARIMA(3, 1, 4)模型对GM(1,1)模型的残差修正值进行拟合，然后将 ARIMA 预测模型差分还原后的数据与GM(1, 1)模型预测数据相加，即为灰色 ARIMA 组合模型对于2021—2025年浙江省水泥产量的最终预测结果，如图5所示。

显示原图|下载原图ZIP|生成PPT

图5 浙江省水泥产量组合模型结果拟合图

从图5可以看出，灰色ARIMA组合模型的拟合精度明显优于GM(1, 1)模型。两种模型的预测值如表6所示。

表6 GM(1, 1)模型和灰色ARIMA组合模型预测值

年份	实际值	GM(1, 1)模型		灰色ARIMA组合模型
年份	实际值	预测值 /万t	误差	预测值 /万t	误差
2016	10 796.5	12 630.798	16.99	10 713.823	0.77
2017	11 231.2	13 127.898	16.888	10 922.706	2.75
2018	12 248.3	13 631.252	11.291	11 909.98	0.02
2019	13 441.04	14 140.938	5.207	13 218.228	0.02
2020	13 272.91	14 657.037	10.428	13 639.809	0.03

由表6可知，灰色ARIMA组合模型的预测精度较高，可用于预测浙江省2021—2025年的水泥产量，预测结果如表7所示。

表7 2021—2025年浙江省水泥产量预测值

年份	灰色ARIMA模型预测值/万t
2021	13 965.22
2022	14 284.1
2023	14 653.1
2024	15 017.5
2025	15 582.3

将上述数据代入式（2）计算可得，浙江省2021—2025的砂石骨料物流需求量，如表8所示。

表8 浙江省2021—2025的砂石骨料物流需求预测值

年份	水泥需求/万t	砂石骨料物流需求/万t
2021	13 965.2	119 701.9
2022	14 284.1	122 434.9
2023	14 653.1	125 598.4
2024	15 017.5	128 721.3
2025	15 582.3	133 562.4

3 结语

本文采用灰色ARIMA组合模型对浙江省砂石骨料物流需求量进行了预测研究。通过对比GM(1, 1)模型与灰色ARIMA组合模型预测结果，验证了灰色ARIMA组合模型相较于GM(1, 1)有更高的精度，确认了灰色ARIMA模型的实用性。运用灰色ARIMA组合模型对浙江省砂石骨料需求进行预测，结果显示2021—2025年浙江省的砂石骨料需求量会继续呈现不断上升的趋势。但本文在研究中未考虑到砂石骨料物流需求量的精准求法，下一步的研究方向是运用大数据算法对砂石骨料的用量进行精确求解，以更好地对砂石骨料需求量进行预测。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

[1]	钟明然, 左清富. 我国建筑砂石骨料供求现状分析及管理对策研究[J]. 桂林航天工业学院学报, 2021, 26(4):504-509.

[2]	BRUZDA J. Quantile smoothing in supply chain and logistics forecasting[J]. International Journal of Production Economics, 2019, 208: 122-139.

[3]	BAISARIYEV M, BAKYTZHANULY A, MUKHANOVA B, et al. Demand forecasting methods for spare parts logistics for aviation: a real-world implementation of the Bootstrap method[J]. Procedia Manufacturing, 2021, 55: 500-506.

[4]	ERWIN V D L, JAN V D, MICHAEL R, et al. Demand forecasting and order planning for humanitarian logistics: an empirical assessment[J]. Journal of Operations Management, 2016, 45(1): 114-122.

[5]	RYUICHI S, TOMIHIRO W. Future forecast of trade amount and international cargo flow in the apec region: an application of trade-logistics forecasting model[J]. Asian Transport Studies, 2012, 2(2): 194-208.

[6]	陈长英. 基于改进的灰色-马尔可夫链模型的广西物流需求预测研究[J]. 西部交通科技, 2019(10):152-155.

[7]	王子健. 基于BP神经网络的关中平原城市群物流竞争力评价[J]. 经济研究导刊, 2022(7):57-59.

[8]	李明书, 刘仁云, 胡姝祺, 等. 基于时间序列分析的长春市物流研究及预测[J]. 科技创新与应用, 2020(32):34-36.

[9]	谭伟华. 基于多元回归和神经网络的江西省物流需求预测研究[D]. 南昌: 江西财经大学, 2020.

[10]	王晓平, 彭文凯, 卢怀宇, 等. 基于支持向量机模型的北京城镇农产品冷链物流需求预测[J]. 湖北农业科学, 2018, 57(15):88-94.

[11]	徐慧. 基于q阶模糊层次分析法的医药应急物流供应能力评价研究[D]. 北京: 北京交通大学, 2021.

[12]	刘庆庆, 郑琰, 孙萌, 等. 基于灰色系统理论的物流需求预测——以连云港港为例[J]. 物流工程与管理, 2022, 44(3):8-11.

[13]	马巧云, 邹相林, 郭柯楠. 基于灰色ARIMA模型的河南省物流需求预测研究[J]. 河南科学, 2019, 37(7):1171-1180.

[14]	曹培格. 基于灰色组合模型的粮食产量预测研究[D]. 郑州: 河南工业大学, 2018.

[15]	王建设, 张皓. 关于北京市砂石骨料铁路绿色供应链发展分析[J]. 物流工程与管理, 2019, 41(5):102-103.

[16]	路正南, 张志娟. 我国风电产业装机容量预测分析[J]. 科技与管理, 2011, 13(1):77-80.

[17]	刘思峰, 郭天榜, 党耀国. 灰色系统理论及其应用[M]. 北京: 科学出版社, 1999.

[18]	杨克磊, 张振宇, 和美. 应用灰色 GM(1, 1)模型的粮食产量预测研究[J]. 重庆理工大学学报, 2015(4):124-127.

[19]	周慧秋. 灰色模型GM(1, N)在东北地区粮食综合生产能力预测中的应用研究[J]. 农业技术经济, 2006(3):58-62.

[20]	李葆春, 马琦. 灰色GM(1, 1)模型在定西县粮食产量预测中的应用[J]. 甘肃农业大学学报, 2005, 40(5):660-663.

[21]	崔毅, 徐伟, 张松涛. 基于灰色马尔科夫模型的青岛冷链物流需求预测[J]. 物流科技, 2019, 42(10):78-82,87.

[22]	SURESH K K, PRIYA S R K. Forecasting sugarcane yield of tamilnadu using ARIMA Models[J]. Sugar Tech, 2011, 13(1): 23-26.

[23]	杨铁军, 杨娜, 朱春华, 等. 一种新的基于 ARIMA 模型的粮食产量预测[J]. 河南工业大学学报(自然科学版), 2015, 36(5):19-22.

[24]	中国砂石协会. 2021年中国砂石行业运行报告[EB/OL]. (2021-02-28)[2022-04-01]. https://mnr.gov.cn/dt/mtsy/202203/t20220301_2729521.html.

[25]	浙江省统计局. 1990—2021年浙江省统计年鉴[EB/OL].(2021-09-01)[2022-04-01]. http://tjj.zj.gov.cn/col/col1525563/index.html.

Options

文章导航

模态框（Modal）标题

摘要

本文引用格式

Abstract

0 引言

1 模型构建

1.1 问题描述

1.2 灰色ARIMA组合模型

表1 ARIMA模型的相关性特征

图1 灰色ARIMA模型的预测流程

2 实证分析

2.1 浙江省砂石骨料物流需求现状

图2 砂石骨料供需关系图

2.2 模型精度验证及求解

图3 浙江省1991—2015年水泥产量

表2 GM(1, 1)模型预测值

图4 浙江省水泥产量GM（1, 1）模型结果拟合图

表3 ADF平稳性检验结果

表4 ADF单位根检验

表5 残差序列的自相关和偏相关函数

图5 浙江省水泥产量组合模型结果拟合图

表6 GM(1, 1)模型和灰色ARIMA组合模型预测值

表7 2021—2025年浙江省水泥产量预测值

表8 浙江省2021—2025的砂石骨料物流需求预测值

3 结语

参考文献